Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multilocus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Microbiological and molecular studies on a multidrug-resistant Pseudomonas aeruginosa from a liver transplant patient with urinary tract infection in Egypt

BACKGROUND

Pseudomonas aeruginosa is an opportunistic pathogen responsible for complicated UTIs and exhibits high antibiotic resistance, leading to increased mortality rates, especially in cases of multidrug-resistant strains. This study aimed to investigate the antibiotic susceptibility patterns and genomic characterization of XDR strains identified in end-stage liver disease patients who underwent liver transplants.

METHODS

In this study, a number of 30 individuals who underwent liver transplants were registered. Ninety urine and 60 wound site swab samples were collected and processed for culturing, identification, and antimicrobial sensitivity. Extensively drug-resistant strain EMARA01 was confirmed through Sanger sequencing and was then processed for whole genome sequencing to characterize the genomic pattern. Sequencing data were processed for de novo assembly using various tools and databases, including genome annotation, serotype identification, virulence factor genes, and antimicrobial resistance gene. Pangenome analysis of randomly selected 147 reference strains and EMAR01 sequenced strain was performed using the Bacterial Pan Genome Analysis (BPGA) software.

RESULTS

Of these total examined samples, nosocomial infection due to P. aeruginosa was detected in twelve patients' samples. AST analysis showed that P. aeruginosa strains exhibit resistance to tobramycin, erythromycin, and gentamicin, followed by piperacillin and ofloxacin, and no strains exhibit resistance to meropenem and imipenem. The CARD database identified 59 AMR genes similar to the EMAR01 strain genome and mostly belong to the family involved in the resistance-nodulation-cell division (RND) antibiotic efflux pump. Five genes; nalC, nalD, MexR, MexA, and MexB, exhibit resistance to 14 classes of antibiotics, while two AMR; CpxR, and OprM, exhibit resistance to 15 classes of drugs. Pangenome analysis revealed that the pan-genome remained open, suggesting the potential for acquiring accessory and unique genes. Notably, the genes predominantly involved in amino acid transport metabolism were identified using the KEGG database.

CONCLUSIONS

This study provides valuable insights into the antimicrobial resistance profile, genetic features, and genomic evolution of P. aeruginosa strains causing UTIs in liver transplant patients. The findings emphasize the significance of comprehending AMR mechanisms and genetic diversity in P. aeruginosa for developing effective treatment strategies and infection control measures.

Biodegradation of phthalate acid esters and whole-genome analysis of a novel Streptomyces sp. FZ201 isolated from natural habitats

Di-n-butyl phthalate (DBP) is one of the most extensively used phthalic acid esters (PAEs) and is considered to be an emerging, globally concerning pollutant. The genus Streptomyces holds promise as a degrader of various organic pollutants, but PAE biodegradation mechanisms by Streptomyces species remain unsolved. In this study, a novel PAE-degrading Streptomyces sp. FZ201 isolated from natural habitats efficiently degraded various PAEs. FZ201 had strong resilience against DBP and exhibited immediate degradation, with kinetics adhering to a first-order model. The comprehensive biodegradation of DBP involves de-esterification, β-oxidation, trans-esterification, and aromatic ring cleavage. FZ201 contains numerous catabolic genes that potentially facilitate PAE biodegradation. The DBP metabolic pathway was reconstructed by genome annotation and intermediate identification. Streptomyces species have an open pangenome with substantial genome expansion events during the evolutionary process, enabling extensive genetic diversity and highly plastic genomes within the Streptomyces genus. FZ201 had a diverse array of highly expressed genes associated with the degradation of PAEs, potentially contributing significantly to its adaptive advantage and efficiency of PAE degradation. Thus, FZ201 is a promising candidate for remediating highly PAE-contaminated environments. These findings enhance our preliminary understanding of the molecular mechanisms employed by Streptomyces for the removal of PAEs.

Tracking and characterization of a novel conjugative transposon identified by shotgun transposon mutagenesis

The horizontal transfer of mobile genetic elements (MGEs) is an essential process determining the functional and genomic diversity of bacterial populations. MGEs facilitate the exchange of fitness determinant genes like antibiotic resistance and virulence factors. Various computational methods exist to identify potential MGEs, but confirming their ability to transfer requires additional experimental approaches. Here, we apply a transposon (Tn) mutagenesis technique for confirming mobilization without the need for targeted mutations. Using this method, we identified two MGEs, including a previously known conjugative transposon (CTn) called BoCTn found in Bacteroides ovatus and a novel CTn, PvCTn, identified in Phocaeicola vulgatus. In addition, Tn mutagenesis and subsequent genetic deletion enabled our characterization of a helix-turn-helix motif gene, BVU3433 which negatively regulates the conjugation efficiency of PvCTn in vitro. Furthermore, our transcriptomics data revealed that BVU3433 plays a crucial role in the repression of PvCTn genes, including genes involved in forming complete conjugation machinery [Type IV Secretion System (T4SS)]. Finally, analysis of individual strain genomes and community metagenomes identified the widespread prevalence of PvCTn-like elements with putative BVU3433 homologs among human gut-associated bacteria. In summary, this Tn mutagenesis mobilization method (TMMM) enables observation of transfer events in vitro and can ultimately be applied in vivo to identify a broader diversity of functional MGEs that may underly the transfer of important fitness determinants.

Comparative genomic analysis of Pseudomonas aeruginosa strains susceptible and resistant to carbapenems and aztreonam isolated from patients with healthcare-associated infections in a Mexican hospital

Pseudomonas aeruginosa (PA) is an important opportunistic pathogen that causes different infections on immunocompromised patients. Within PA accessory genome, differences in virulence, antibiotic resistance and biofilm formation have been described between strains, leading to the emergence of multidrug-resistant strains. The genome sequences of 17 strains isolated from patients with healthcare-associated infections in a Mexican hospital were genomically and phylogenetically analyzed and antibiotic resistance genes, virulence genes, and biofilm formation genes were detected. Fifteen of the 17 strains were resistant to at least two of the carbapenems meropenem, imipenem, and the monobactam aztreonam. The antibiotic resistance (mexA, mexB, and oprM) and the biofilm formation (pslA and pslD) genes were detected in all strains. Differences were found between strains in accessory genome size. The strains had different sequence types, and seven strains had sequence types associated with global high risk epidemic PA clones. All strains were represented in two groups among PA global strains. In the 17 strains, horizontally acquired resistance genes to aminoglycosides and beta-lactams were found, mainly, and between 230 and 240 genes that encode virulence factors. The strains under study were variable in terms of their accessory genome, antibiotic resistance, and virulence genes. With these characteristics, we provide information about the genomic diversity of clinically relevant PA strains.

Comparative genomics of Brucella abortus and Brucella melitensis unravels the gene sharing, virulence factors and SNP diversity among the standard, vaccine and field strains

Brucella abortus and Brucella melitensis are the primary etiological agents of brucellosis in large and small ruminants, respectively. There are limited comparative genomic studies involving Brucella strains that explore the relatedness among both species. In this study, we involved strains (n=44) representing standard, vaccine and Indian field origin for pangenome, single nucleotide polymorphism (SNP) and phylogenetic analysis. Both species shared a common gene pool representing 2884 genes out of a total 3244 genes. SNP-based phylogenetic analysis indicated higher SNP diversity among B. melitensis (3824) strains in comparison to B. abortus (540) strains, and a clear demarcation was identified between standard/vaccine and field strains. The analysis for virulence genes revealed that virB3, virB7, ricA, virB5, ipx5, wbkC, wbkB, and acpXL genes were highly conserved in most of the Brucella strains. Interestingly, virB10 gene was found to have high variability among the B. abortus strains. The cgMLST analysis revealed distinct sequence types for the standard/vaccine and field strains. B. abortus strains from north-eastern India fall within similar sequence type differing from other strains. In conclusion, the analysis revealed a highly shared core genome among two Brucella species. SNP analysis revealed B. melitensis strains exhibit high diversity as compared to B. abortus strains. Strains with absence or high polymorphism of virulence genes can be exploited for the development of novel vaccine candidates effective against both B. abortus and B. melitensis.

Delving into the lifestyle of Sundarban Wetland resident, biofilm producing, halotolerant Salinicoccus roseus: a comparative genomics-based intervention

BACKGROUND

Microbial community played an essential role in ecosystem processes, be it mangrove wetland or other intertidal ecologies. Several enzymatic activities like hydrolases are effective ecological indicators of soil microbial function. So far, little is known on halophilic bacterial contribution and function on a genomic viewpoint of Indian Sundarban Wetland. Considering the above mentioned issues, the aims of this study was to understand the life style, metabolic functionalities and genomic features of the isolated bacterium, Salinicoccus roseus strain RF1H. A comparative genome-based study of S. roseus has not been reported yet. Henceforth, we have considered the inclusion of the intra-species genome comparison of S. roseus to gain insight into the high degree of variation in the genome of strain RF1H among others.

RESULTS

Salinicoccus roseus strain RF1H is a pink-red pigmented, Gram-positive and non-motile cocci. The bacterium exhibited high salt tolerance (up to 15% NaCl), antibiotic resistance, biofilm formation and secretion of extracellular hydrolytic enzymes. The circular genome was approximately 2.62978 Mb in size, encoding 574 predicted genes with GC content 49.5%. Presence of genomic elements (prophages, transposable elements, CRISPR-Cas system) represented bacterial virulence and multidrug-resistance. Furthermore, genes associated with salt tolerance, temperature adaptation and DNA repair system were distributed in 17 genomic islands. Genes related to hydrocarbon degradation manifested metabolic capability of the bacterium for potential biotechnological applications. A comparative pangenome analysis revealed two-component response regulator, modified C4-dicarboxylate transport system and osmotic stress regulated ATP-binding proteins. Presence of genes encoding arginine decarboxylase (ADC) enzyme being involved in biofilm formation was reported from the genome. In silico study revealed the protein is thermostable and made up with ~ 415 amino acids, and hydrophilic in nature. Three motifs appeared to be evolutionary conserved in all Salinicoccus sequences.

CONCLUSION

The first report of whole genome analysis of Salinicoccus roseus strain RF1H provided information of metabolic functionalities, biofilm formation, resistance mechanism and adaptation strategies to thrive in climate-change induced vulnerable spot like Sundarban. Comparative genome analysis highlighted the unique genome content that contributed the strain's adaptability. The biomolecules produced during metabolism are important sources of compounds with potential beneficial applications in pharmaceuticals.

Local and Environmental Reservoirs of Salmonella enterica After Hurricane Florence Flooding

In many regions of the world, including the United States, human and animal fecal genetic markers have been found in flood waters. In this study, we use high-resolution whole genomic sequencing to examine the origin and distribution of Salmonella enterica after the 2018 Hurricane Florence flooding. We specifically asked whether S. enterica isolated from water samples collected near swine farms in North Carolina shortly after Hurricane Florence had evidence of swine origin. To investigate this, we isolated and fully sequenced 18 independent S. enterica strains from 10 locations (five flooded and five unflooded). We found that all strains have extremely similar chromosomes with only five single nucleotide polymorphisms (SNPs) and possessed two plasmids assigned bioinformatically to the incompatibility groups IncFIB and IncFII. The chromosomal core genome and the IncFIB plasmid are most closely related to environmental Salmonella strains isolated previously from the southeastern US. In contrast, the IncFII plasmid was found in environmental S. enterica strains whose genomes were more divergent, suggesting the IncFII plasmid is more promiscuous than the IncFIB type. We identified 65 antibiotic resistance genes (ARGs) in each of our 18 S. enterica isolates. All ARGs were located on the Salmonella chromosome, similar to other previously characterized environmental isolates. All isolates with different SNPs were resistant to a panel of commonly used antibiotics. These results highlight the importance of environmental sources of antibiotic-resistant S. enterica after extreme flood events.

Genome-wide screens reveal shared and strain-specific genes that facilitate enteric colonization by Klebsiella pneumoniae

Gastrointestinal (GI) colonization by Klebsiella pneumoniae is a risk factor for subsequent infection as well as transmission to other patients. Additionally, colonization is achieved by many strain types that exhibit high diversity in genetic content. Thus, we aimed to study strain-specific requirements for K. pneumoniae GI colonization by applying transposon insertion sequencing to three classical clinical strains: a carbapenem-resistant strain, an extended-spectrum beta-lactamase-producing strain, and a non-epidemic antibiotic-susceptible strain. The transposon insertion libraries were screened in a murine model of GI colonization. At 3 days post-inoculation, 27 genes were required by all three strains for colonization. Isogenic deletion mutants for three genes/operons (acrA, carAB, and tatABCD) confirmed colonization defects in each of the three strains. Additionally, deletion of acrA reduced bile tolerance in vitro, while complementation restored both bile tolerance in vitro and colonization ability in vivo. Transposon insertion sequencing suggested that some genes were more important for the colonization of one strain than the others. For example, deletion of the sucrose porin-encoding gene scrY resulted in a colonization defect in the carbapenemase-producing strain but not in the extended-spectrum beta-lactamase producer or the antibiotic-susceptible strain. These findings demonstrate that classical K. pneumoniae strains use both shared and strain-specific strategies to colonize the mouse GI tract. IMPORTANCE Klebsiella pneumoniae is a common cause of difficult-to-treat infections due to its propensity to express resistance to many antibiotics. For example, carbapenem-resistant K. pneumoniae has been named an urgent threat by the United States Centers for Disease Control and Prevention. Gastrointestinal colonization in patients with K. pneumoniae has been linked to subsequent infection, making it a key process to control in the prevention of multidrug-resistant infections. However, the bacterial factors which contribute to K. pneumoniae colonization are not well understood. Additionally, individual strains exhibit large amounts of genetic diversity, begging the question of whether some colonization factors are strain dependent. This study identifies the enteric colonization factors of three classical strains using transposon mutant screens to define a core colonization program for K. pneumoniae as well as detecting strain-to-strain differences in colonization strategies.

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multi-locus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Genome-wide screens reveal shared and strain-specific genes that facilitate enteric colonization by Klebsiella pneumoniae

Gastrointestinal (GI) colonization by Klebsiella pneumoniae is a risk factor for subsequent infection as well as transmission to other patients. Additionally, colonization is achieved by many strain types that exhibit high diversity in genetic content. Thus, we aimed to study strain-specific requirements for K. pneumoniae GI colonization by applying transposon insertion sequencing to three classical clinical strains: a carbapenem-resistant strain, an extended-spectrum beta-lactamase producing strain, and a non-epidemic antibiotic-susceptible strain. The transposon insertion libraries were screened in a murine model of GI colonization. At three days post-inoculation, 27 genes were required by all three strains for colonization. Isogenic deletion mutants for three genes/operons (acrA, carAB, tatABCD) confirmed colonization defects in each of the three strains. Additionally, deletion of acrA reduced bile tolerance in vitro, while complementation restored both bile tolerance in vitro and colonization ability in vivo. Transposon insertion sequencing suggested that some genes were more important for colonization of one strain than the others. For example, deletion of the sucrose porin-encoding gene scrY resulted in a colonization defect in the carbapenemase-producing strain but not in the extended-spectrum beta-lactamase producer or the antibiotic-susceptible strain. These findings demonstrate that classical K. pneumoniae strains use both shared and strain-specific strategies to colonize the mouse GI tract.

IMPORTANCE

Klebsiella pneumoniae is a common cause of difficult-to-treat infections due to its propensity to express resistance to many antibiotics. For example, carbapenem-resistant K. pneumoniae (CR-Kp) has been named an urgent threat by the United States Centers for Disease Control and Prevention. Gastrointestinal colonization of patients with K. pneumoniae has been linked to subsequent infection, making it a key process to control in prevention of multidrug-resistant infections. However, the bacterial factors which contribute to K. pneumoniae colonization are not well understood. Additionally, individual strains exhibit large amounts of genetic diversity, begging the question of whether some colonization factors are strain-dependent. This study identifies the enteric colonization factors of 3 classical strains using transposon mutant screens to define a core colonization program for K. pneumoniae as well as detecting strain-to-strain differences in colonization strategies.

CRISPR Comparison Toolkit: Rapid Identification, Visualization, and Analysis of CRISPR Array Diversity

CRISPR-Cas systems provide immunity against mobile genetic elements (MGEs) through sequence-specific targeting by spacer sequences encoded in CRISPR arrays. Spacers are highly variable between microbial strains and can be acquired rapidly, making them well suited for use in strain typing of closely related organisms. However, no tools are currently available to automate the process of reconstructing strain histories using CRISPR spacers. We therefore developed the CRISPR Comparison Toolkit (CCTK) to enable analyses of array relationships. The CCTK includes tools to identify arrays, analyze relationships between arrays using CRISPRdiff and CRISPRtree, and predict targets of spacers. CRISPRdiff visualizes arrays and highlights the similarities between them. CRISPRtree infers a phylogenetic tree from array relationships and presents a hypothesis of the evolutionary history of the arrays. The CCTK unifies several CRISPR analysis tools into a single command line application, including the first tool to infer phylogenies from array relationships.

Compendium-Wide Analysis of Pseudomonas aeruginosa Core and Accessory Genes Reveals Transcriptional Patterns across Strains PAO1 and PA14

Pseudomonas aeruginosa is an opportunistic pathogen that causes difficult-to-treat infections. Two well-studied divergent P. aeruginosa strain types, PAO1 and PA14, have significant genomic heterogeneity, including diverse accessory genes present in only some strains. Genome content comparisons find core genes that are conserved across both PAO1 and PA14 strains and accessory genes that are present in only a subset of PAO1 and PA14 strains. Here, we use recently assembled transcriptome compendia of publicly available P. aeruginosa RNA sequencing (RNA-seq) samples to create two smaller compendia consisting of only strain PAO1 or strain PA14 samples with each aligned to their cognate reference genome. We confirmed strain annotations and identified other samples for inclusion by assessing each sample's median expression of PAO1-only or PA14-only accessory genes. We then compared the patterns of core gene expression in each strain. To do so, we developed a method by which we analyzed genes in terms of which genes showed similar expression patterns across strain types. We found that some core genes had consistent correlated expression patterns across both compendia, while others were less stable in an interstrain comparison. For each accessory gene, we also determined core genes with correlated expression patterns. We found that stable core genes had fewer coexpressed neighbors that were accessory genes. Overall, this approach for analyzing expression patterns across strain types can be extended to other groups of genes, like phage genes, or applied for analyzing patterns beyond groups of strains, such as samples with different traits, to reveal a deeper understanding of regulation. IMPORTANCE Pseudomonas aeruginosa is a ubiquitous pathogen. There is much diversity among P. aeruginosa strains, including two divergent but well-studied strains, PAO1 and PA14. Understanding how these different strain-level traits manifest is important for identifying targets that regulate different traits of interest. With the availability of thousands of PAO1 and PA14 samples, we created two strain-specific RNA-seq compendia where each one contains hundreds of samples from PAO1 or PA14 strains and used them to compare the expression patterns of core genes that are conserved in both strain types and to determine which core genes have expression patterns that are similar to those of accessory genes that are unique to one strain or the other using an approach that we developed. We found a subset of core genes with different transcriptional patterns across PAO1 and PA14 strains and identified those core genes with expression patterns similar to those of strain-specific accessory genes.

A Structural Systems Biology Approach to High-Risk CG23 Klebsiella pneumoniae

Klebsiella pneumoniae is a leading cause of antibiotic-resistant-associated deaths in the world. Here, we report the deposition of 14 structures of enzymes from both the core and accessory genomes of sequence type 23 (ST23) K1 hypervirulent K. pneumoniae.

Integrated conjugative plasmid drives high frequency chromosomal gene transfer in Sulfolobus islandicus

Gene transfer in crenarchaea has been observed within natural and experimental populations of Sulfolobus. However, the molecular factors that govern how gene transfer and recombination manifest themselves in these populations is still unknown. In this study, we examine a plasmid-mediated mechanism of gene transfer in S. islandicus that results in localized high frequency recombination within the chromosome. Through chromosomal marker exchange assays with defined donors and recipients, we find that while bidirectional exchange occurs among all cells, those possessing the integrated conjugative plasmid, pM164, mobilize a nearby locus at a significantly higher frequency when compared to a more distal marker. We establish that traG is essential for this phenotype and that high frequency recombination can be replicated in transconjugants after plasmid transfer. Mapping recombinants through genomic analysis, we establish the distribution of recombinant tracts with decreasing frequency at increasing distance from pM164. We suggest the bias in transfer is a result of an Hfr (high frequency recombination)-like conjugation mechanism in this strain. In addition, we find recombinants containing distal non-selected recombination events, potentially mediated by a different host-encoded marker exchange (ME) mechanism.

Diversity and prevalence of type VI secretion system effectors in clinical Pseudomonas aeruginosa isolates

Pseudomonas aeruginosa is an opportunistic pathogen and a major driver of morbidity and mortality in people with Cystic Fibrosis (CF). The Type VI secretion system (T6SS) is a molecular nanomachine that translocates effectors across the bacterial membrane into target cells or the extracellular environment enabling intermicrobial interaction. P. aeruginosa encodes three T6SS clusters, the H1-, H2- and H3-T6SS, and numerous orphan islands. Genetic diversity of T6SS-associated effectors in P. aeruginosa has been noted in reference strains but has yet to be explored in clinical isolates. Here, we perform a comprehensive bioinformatic analysis of the pangenome and T6SS effector genes in 52 high-quality clinical P. aeruginosa genomes isolated from CF patients and housed in the Personalised Approach to P. aeruginosa strain repository. We confirm that the clinical CF isolate pangenome is open and principally made up of accessory and unique genes that may provide strain-specific advantages. We observed genetic variability in some effector/immunity encoding genes and show that several well-characterised vgrG and PAAR islands are absent from numerous isolates. Our analysis shows clear evidence of disruption to T6SS genomic loci through transposon, prophage, and mobile genetic element insertions. We identified an orphan vgrG island in P. aeruginosa strain PAK and five clinical isolates using in silico analysis which we denote vgrG7, predicting a gene within this cluster to encode a Tle2 lipase family effector. Close comparison of T6SS loci in clinical isolates compared to reference P. aeruginosa strain PAO1 revealed the presence of genes encoding eight new T6SS effectors with the following putative functions: cytidine deaminase, lipase, metallopeptidase, NADase, and pyocin. Finally, the prevalence of characterised and putative T6SS effectors were assessed in 532 publicly available P. aeruginosa genomes, which suggests the existence of accessory effectors. Our in silico study of the P. aeruginosa T6SS exposes a level of genetic diversity at T6SS genomic loci not seen to date within P. aeruginosa, particularly in CF isolates. As understanding the effector repertoire is key to identifying the targets of T6SSs and its efficacy, this comprehensive analysis provides a path for future experimental characterisation of these mediators of intermicrobial competition and host manipulation.

Overcoming Barriers to Preventing and Treating P. aeruginosa Infections Using AAV Vectored Immunoprophylaxis

Pseudomonas aeruginosa is a bacterial pathogen of global concern and is responsible for 10-15% of nosocomial infections worldwide. This opportunistic bacterial pathogen is known to cause serious complications in immunocompromised patients and is notably the leading cause of morbidity and mortality in patients suffering from cystic fibrosis. Currently, the only line of defense against P. aeruginosa infections is antibiotic treatment. Due to the acquired and adaptive resistance mechanisms of this pathogen, the prevalence of multidrug resistant P. aeruginosa strains has increased, presenting a major problem in healthcare settings. To date, there are no approved licensed vaccines to protect against P. aeruginosa infections, prompting the urgent need alternative treatment options. An alternative to traditional vaccines is vectored immunoprophylaxis (VIP), which utilizes a safe and effective adeno-associated virus (AAV) gene therapy vector to produce sustained levels of therapeutic monoclonal antibodies (mAbs) in vivo from a single intramuscular injection. In this review, we will provide an overview of P. aeruginosa biology and key mechanisms of pathogenesis, discuss current and emerging treatment strategies for P. aeruginosa infections and highlight AAV-VIP as a promising novel therapeutic platform.

Surveillance of Enterobacter cloacae complex colonization and comparative analysis of different typing methods on a neonatal intensive care unit in Germany

Background

Enterobacter cloacae complex is a group of common opportunistic pathogens on neonatal intensive care units. Active microbiological screening to guide empirical antimicrobial treatment or to detect transmission events is recommended in high-risk preterm neonates. A rise in colonization with E. cloacae complex was observed in a German perinatal centre. The aim of this study was to evaluate the performance of different typing techniques using whole genome sequencing (WGS) as a reference.

Methods

Enterobacter cloacae complex isolates from clinical and screening specimens with an epidemiological link to the neonatal intensive care units were further assessed. Identification and antibiotic susceptibility testing was performed by a combination of VITEK2 (bioMérieux) and MALDI-TOF (Bruker Daltonics), followed by RAPD/rep-PCR and PFGE (XbaI). Retrospectively, all isolates were analyzed by Fourier-transform infrared (FTIR) spectroscopy (IR Biotyper, Bruker Daltonics). Whole genome sequencing with SNP-based clustering was used as the reference method. Furthermore, resistome analysis, sequence type and species identification were derived from the WGS data. Transmission analysis was based on epidemiological and typing data.

Results

Between September 2017 and March 2018 32 mostly preterm neonates were found to be colonized with E. cloacae complex and 32 isolates from 24 patients were available for further typing. RAPD/rep-PCR and PFGE showed good concordance with WGS whereas FTIR displayed mediocre results [adjusted rand index (ARI) = 0.436]. A polyclonal increase and two dominant and overlapping clonal clusters of two different E. hormaechei subspecies were detected. Overall, four different species were identified. Genotyping confirmed third-generation cephalosporin resistance development in isolates of the same patient. During the six-month period several infection prevention interventions were performed and no E. cloacae complex isolates were observed during the following months.

Conclusions

Interpretation of the microbiological results alone to detect transmission events is often challenging and bacterial typing is of utmost importance to implement targeted infection control measures in an epidemic occurrence of E. cloacae complex. WGS is the most discriminatory method. However, traditional methods such as PFGE or RAPD/rep-PCR can provide reliable and quicker results in many settings. Furthermore, research is needed to quickly identify E. cloacae complex to the species level in the microbiological laboratory.

Pseudomonas aeruginosa reference strains PAO1 and PA14: A genomic, phenotypic, and therapeutic review

Pseudomonas aeruginosa is a ubiquitous, motile, gram-negative bacterium that has been recently identified as a multi-drug resistant pathogen in critical need of novel therapeutics. Of the approximately 5,000 strains, PAO1 and PA14 are common laboratory reference strains, modeling moderately and hyper-virulent phenotypes, respectively. PAO1 and PA14 have been instrumental in facilitating the discovery of novel drug targets, testing novel therapeutics, and supplying critical genomic information on the bacterium. While the two strains have contributed to a wide breadth of knowledge on the natural behaviors and therapeutic susceptibilities of P. aeruginosa, they have demonstrated significant deviations from observations in human infections. Many of these deviations are related to experimental inconsistencies in laboratory strain environment that complicate and, at times, terminate translation from laboratory results to clinical applications. This review aims to provide a comparative analysis of the two strains and potential methods to improve their clinical relevance.

Genome characterization of a uropathogenic Pseudomonas aeruginosa isolate PA_HN002 with cyclic di-GMP-dependent hyper-biofilm production

Pseudomonas aeruginosa can cause various types of infections and is one of the most ubiquitous antibiotic-resistant pathogens found in healthcare settings. It is capable of adapting to adverse conditions by transforming its motile lifestyle to a sessile biofilm lifestyle, which induces a steady state of chronic infection. However, mechanisms triggering the lifestyle transition of P. aeruginosa strains with clinical significance are not very clear. In this study, we reported a recently isolated uropathogenic hyper-biofilm producer PA_HN002 and characterized its genome to explore genetic factors that may promote its transition into the biofilm lifestyle. We first showed that high intracellular c-di-GMP content in PA_HN002 gave rise to its attenuated motilities and extraordinary strong biofilm. Reducing the intracellular c-di-GMP content by overexpressing phosphodiesterases (PDEs) such as BifA or W909_14950 converted the biofilm and motility phenotypes. Whole genome sequencing and comprehensive analysis of all the c-di-GMP metabolizing enzymes led to the identification of multiple mutations within PDEs. Gene expression assays further indicated that the shifted expression profile of c-di-GMP metabolizing enzymes in PA_HN002 might mainly contribute to its elevated production of intracellular c-di-GMP and enhanced biofilm formation. Moreover, mobile genetic elements which might interfere the endogenous regulatory network of c-di-GMP metabolism in PA_HN002 were analyzed. This study showed a reprogrammed expression profile of c-di-GMP metabolizing enzymes which may promote the pathoadaption of clinical P. aeruginosa into biofilm producers.

A genome-led study on the pathogenesis of Fusobacterium necrophorum infections

Fusobacterium necrophorum causes a range of mild to life threatening infections and there is uncertainty in terms of diagnosis and treatment due to the lack of knowledge on their pathogenic mechanisms. This study characterised genomes of F. necrophorum to compare their virulence factors and investigate potential infection markers. 27 isolates of F. necrophorum from patients with pharyngotonsillitis were subjected to whole genome sequencing and compared with 42 genomes published in the NCBI database. Phylogenomics, pangemome, pan-GWAS and virulome were analysed to study strain variations with reference to virulence factors. Core genome based phylogenomic tree exhibited three clades of which Clade A belonged to F. necrophorum subsp necrophorum, clades B and C were F. necrophorum subsp funduliforme. Pan-GWAS and Pan-Virulome suggest some marker genes associated with clinical sources of isolation that needs further validation. Our study highlights some interesting features of the pathogenesis of F. necrophorum infections. Although the animal isolate genomes had some marker genes, the genomes of human isolates did not exhibit clear correlation to their clinical sources of isolation. This prompts to think of other mechanisms such as co-infections or host factors that can be involved in the pathogenesis.

Phylogenetic analysis and accessory genome diversity reveal insight into the evolutionary history of Streptococcus dysgalactiae

Streptococcus dysgalactiae (SD) is capable of infecting both humans and animals and causing a wide range of invasive and non-invasive infections. With two subspecies, the taxonomic status of subspecies of SD remains controversial. Subspecies equisimilis (SDSE) is an important human pathogen, while subspecies dysgalactiae (SDSD) has been considered a strictly animal pathogen; however, occasional human infections by this subspecies have been reported in the last few years. Moreover, the differences between the adaptation of SDSD within humans and other animals are still unknown. In this work, we provide a phylogenomic analysis based on the single-copy core genome of 106 isolates from both the subspecies and different infected hosts (animal and human hosts). The accessory genome of this species was also analyzed for screening of genes that could be specifically involved with adaptation to different hosts. Additionally, we searched putatively adaptive traits among prophage regions to infer the importance of transduction in the adaptation of SD to different hosts. Core genome phylogenetic relationships segregate all human SDSE in a single cluster separated from animal SD isolates. The subgroup of bovine SDSD evolved from this later clade and harbors a specialized accessory genome characterized by the presence of specific virulence determinants (e.g., cspZ) and carbohydrate metabolic functions (e.g., fructose operon). Together, our results indicate a host-specific SD and the existence of an SDSD group that causes human–animal cluster infections may be due to opportunistic infections, and that the exact incidence of SDSD human infections may be underestimated due to failures in identification based on the hemolytic patterns. However, more detailed research into the isolation of human SD is needed to assess whether it is a carrier phenomenon or whether the species can be permanently integrated into the human microbiome, making it ready to cause opportunistic infections.

Harnessing diversity of Lactococcus lactis from raw goat milk: Design of an indigenous starter for the production of Rocamadour, a French PDO cheese

Twenty-four strains of Lactococcus lactis isolated from raw goat milk collected in the Rocamadour PDO area were analysed by MLST typing and phenotypic characterisation. The strains were combined to design an indigenous starter for the production of Rocamadour PDO cheese. The strains were divided into three classes based on their technological properties: acidifying and proteolytic strains in class I (12/24 strains), slightly acidifying and non-proteolytic strains in class II (2/24 strains), and non-acidifying and non-proteolytic strains in class III (10/24 strains). Interestingly, all but three strains (21/24) produced diacetyl/acetoin despite not having citrate metabolism genes, as would classically be expected for the production of these aroma compounds. Three strains (EIP07A, EIP13D, and EIP20B) were selected for the indigenous starter based on the following inclusion/exclusion criteria: (i) no negative interactions between included strains, (ii) ability to metabolize lactose and at least one strain with the prtP gene and/or capable of producing diacetyl/acetoin, and (iii) selected strains derived from different farms to maximise genetic and phenotypic diversity. Despite consisting exclusively of L. lactis strains, the designed indigenous starter allowed reproducible cheese production with performances similar to those obtained with an industrial starter and with the sensory qualities expected of Rocamadour PDO cheese.

Intra-Species Genomic Variation in the Pine Pathogen Fusarium circinatum

Fusarium circinatum is an important global pathogen of pine trees. Genome plasticity has been observed in different isolates of the fungus, but no genome comparisons are available. To address this gap, we sequenced and assembled to chromosome level five isolates of F. circinatum. These genomes were analysed together with previously published genomes of F. circinatum isolates, FSP34 and KS17. Multi-sample variant calling identified a total of 461,683 micro variants (SNPs and small indels) and a total of 1828 macro structural variants of which 1717 were copy number variants and 111 were inversions. The variant density was higher on the sub-telomeric regions of chromosomes. Variant annotation revealed that genes involved in transcription, transport, metabolism and transmembrane proteins were overrepresented in gene sets that were affected by high impact variants. A core genome representing genomic elements that were conserved in all the isolates and a non-redundant pangenome representing all genomic elements is presented. Whole genome alignments showed that an average of 93% of the genomic elements were present in all isolates. The results of this study reveal that some genomic elements are not conserved within the isolates and some variants are high impact. The described genome-scale variations will help to inform novel disease management strategies against the pathogen.

Genomic, morphological, and biochemical analyses of a multi-metal resistant but multi-drug susceptible strain of Bordetella petrii from hospital soil

Contamination of soil by antibiotics and heavy metals originating from hospital facilities has emerged as a major cause for the development of resistant microbes. We collected soil samples surrounding a hospital effluent and measured the resistance of bacterial isolates against multiple antibiotics and heavy metals. One strain BMCSI 3 was found to be sensitive to all tested antibiotics. However, it was resistant to many heavy metals and metalloids like cadmium, chromium, copper, mercury, arsenic, and others. This strain was motile and potentially spore-forming. Whole-genome shotgun assembly of BMCSI 3 produced 4.95 Mb genome with 4,638 protein-coding genes. The taxonomic and phylogenetic analysis revealed it, to be a Bordetella petrii strain. Multiple genomic islands carrying mobile genetic elements; coding for heavy metal resistant genes, response regulators or transcription factors, transporters, and multi-drug efflux pumps were identified from the genome. A comparative genomic analysis of BMCSI 3 with annotated genomes of other free-living B. petrii revealed the presence of multiple transposable elements and several genes involved in stress response and metabolism. This study provides insights into how genomic reorganization and plasticity results in evolution of heavy metals resistance by acquiring genes from its natural environment.

Pangenomics in Microbial and Crop Research: Progress, Applications, and Perspectives

Advances in sequencing technologies and bioinformatics tools have fueled a renewed interest in whole genome sequencing efforts in many organisms. The growing availability of multiple genome sequences has advanced our understanding of the within-species diversity, in the form of a pangenome. Pangenomics has opened new avenues for future research such as allowing dissection of complex molecular mechanisms and increased confidence in genome mapping. To comprehensively capture the genetic diversity for improving plant performance, the pangenome concept is further extended from species to genus level by the inclusion of wild species, constituting a super-pangenome. Characterization of pangenome has implications for both basic and applied research. The concept of pangenome has transformed the way biological questions are addressed. From understanding evolution and adaptation to elucidating host–pathogen interactions, finding novel genes or breeding targets to aid crop improvement to design effective vaccines for human prophylaxis, the increasing availability of the pangenome has revolutionized several aspects of biological research. The future availability of high-resolution pangenomes based on reference-level near-complete genome assemblies would greatly improve our ability to address complex biological problems.

Natural Saccharomyces cerevisiae Strain Reveals Peculiar Genomic Traits for Starch-to-Bioethanol Production: the Design of an Amylolytic Consolidated Bioprocessing Yeast

Natural yeast with superior fermentative traits can serve as a platform for the development of recombinant strains that can be used to improve the sustainability of bioethanol production from starch. This process will benefit from a consolidated bioprocessing (CBP) approach where an engineered strain producing amylases directly converts starch into ethanol. The yeast Saccharomyces cerevisiae L20, previously selected as outperforming the benchmark yeast Ethanol Red, was here subjected to a comparative genomic investigation using a dataset of industrial S. cerevisiae strains. Along with Ethanol Red, strain L20 was then engineered for the expression of α-amylase amyA and glucoamylase glaA genes from Aspergillus tubingensis by employing two different approaches (delta integration and CRISPR/Cas9). A correlation between the number of integrated copies and the hydrolytic abilities of the recombinants was investigated. L20 demonstrated important traits for the construction of a proficient CBP yeast. Despite showing a close relatedness to commercial wine yeast and the benchmark Ethanol Red, a unique profile of gene copy number variations (CNVs) was found in L20, mainly encoding membrane transporters and secretion pathway proteins but also the fermentative metabolism. Moreover, the genome annotation disclosed seven open reading frames (ORFs) in L20 that are absent in the reference S288C genome. Genome engineering was successfully implemented for amylase production. However, with equal amylase gene copies, L20 proved its proficiency as a good enzyme secretor by exhibiting a markedly higher amylolytic activity than Ethanol Red, in compliance to the findings of the genomic exploration. The recombinant L20 dT8 exhibited the highest amylolytic activity and produced more than 4 g/L of ethanol from 2% starch in a CBP setting without the addition of supplementary enzymes. Based on the performance of this strain, an amylase/glucoamylase ratio of 1:2.5 was suggested as baseline for further improvement of the CBP ability. Overall, L20 showed important traits for the future construction of a proficient CBP yeast. As such, this work shows that natural S. cerevisiae strains can be used for the expression of foreign secreted enzymes, paving the way to strain improvement for the starch-to-bioethanol route.

Differences in the Accessory Genomes and Methylomes of Strains of Streptococcus equi subsp. equi and of Streptococcus equi subsp. zooepidemicus Obtained from the Respiratory Tract of Horses from Texas

Streptococcus equi subsp. equi (SEE) is a host-restricted equine pathogen considered to have evolved from Streptococcus equi subsp. zooepidemicus (SEZ). SEZ is promiscuous in host range and is commonly recovered from horses as a commensal. Comparison of a single strain each of SEE and SEZ using whole-genome sequencing, supplemented by PCR of selected genes in additional SEE and SEZ strains, was used to characterize the evolution of SEE. But the known genetic variability of SEZ warrants comparison of the whole genomes of multiple SEE and SEZ strains. To fill this knowledge gap, we utilized whole-genome sequencing to characterize the accessory genome elements (AGEs; i.e., elements present in some SEE strains but absent in SEZ or vice versa) and methylomes of 50 SEE and 50 SEZ isolates from Texas. Consistent with previous findings, AGEs consistently found in all SEE isolates were primarily from mobile genetic elements that might contribute to host restriction or pathogenesis of SEE. Fewer AGEs were identified in SEZ because of the greater genomic variability among these isolates. The global methylation patterns of SEE isolates were more consistent than those of the SEZ isolates. Among homologous genes of SEE and SEZ, differential methylation was identified only in genes of SEE encoding proteins with functions of quorum sensing, exopeptidase activity, and transitional metal ion binding. Our results indicate that effects of genetic mobile elements in SEE and differential methylation of genes shared by SEE and SEZ might contribute to the host specificity of SEE.

IMPORTANCE Strangles, caused by the host-specific bacterium Streptococcus equi subsp. equi (SEE), is the most commonly diagnosed infectious disease of horses worldwide. Its ancestor, Streptococcus equi subsp. zooepidemicus (SEZ), is frequently isolated from a wide array of hosts, including horses and humans. A comparison of the genomes of a single strain of SEE and SEZ has been reported, but sequencing of further isolates has revealed variability among SEZ strains. Thus, the importance of this study is that it characterizes genomic and methylomic differences of multiple SEE and SEZ isolates from a common geographic region (viz., Texas). Our results affirm many of the previously described differences between the genomes of SEE and SEZ, including the role of mobile genetic elements in contributing to host restriction. We also provide the first characterization of the global methylome of Streptococcus equi and evidence that differential methylation might contribute to the host restriction of SEE.

Comparative genomics and proteomics of Eubacterium maltosivorans: functional identification of trimethylamine methyltransferases and bacterial microcompartments in a human intestinal bacterium with a versatile lifestyle

Eubacterium maltosivorans YI^T is a human intestinal isolate capable of acetogenic, propionogenic and butyrogenic growth. Its 4.3-Mb genome sequence contains coding sequences for 4227 proteins, including 41 different methyltransferases. Comparative proteomics of strain YI^T showed the Wood-Ljungdahl pathway proteins to be actively produced during homoacetogenic growth on H₂ and CO₂ while butyrogenic growth on a mixture of lactate and acetate significantly upregulated the production of proteins encoded by the recently identified lctABCDEF cluster and accessory proteins. Growth on H₂ and CO₂ unexpectedly induced the production of two related trimethylamine methyltransferases. Moreover, a set of 16 different trimethylamine methyltransferases together with proteins for bacterial microcompartments were produced during growth and deamination of the quaternary amines, betaine, carnitine and choline. Growth of strain YI^T on 1,2-propanediol generated propionate with propanol and induced the formation of bacterial microcompartments that were also prominently visible in betaine-grown cells. The present study demonstrates that E. maltosivorans is highly versatile in converting low-energy fermentation end-products in the human gut into butyrate and propionate whilst being capable of preventing the formation of the undesired trimethylamine by converting betaine and other quaternary amines in bacterial microcompartments into acetate and butyrate.

Bacillus velezensis: A Treasure House of Bioactive Compounds of Medicinal, Biocontrol and Environmental Importance

Bacillus velezensis gram-positive bacterium, is frequently isolated from diverse niches mainly soil, water, plant roots, and fermented foods. B. velezensis is ubiquitous, non-pathogenic and endospore forming. Being frequently isolated from diverse plant holobionts it is considered host adapted microorganism and recognized of high economic importance given its ability to promote plant growth under diverse biotic and abiotic stress conditions. Additionally, the species suppress many plant diseases, including bacterial, oomycete, and fungal diseases. It is also able after plant host root colonization to induce unique physiological situation of host plant called primed state. Primed host plants are able to respond more rapidly and/or effectively to biotic or abiotic stress. Moreover, B. velezenis have the ability to resist diverse environmental stresses and help host plants to cope with, including metal and xenobiotic stresses. Within species B. velezensis strains have unique abilities allowing them to adopt different life styles. Strain level abilities knowledge is warranted and could be inferred using the ever-expanding new genomes list available in genomes databases. Pangenome analysis and subsequent identification of core, accessory and unique genomes is actually of paramount importance to decipher species full metabolic capacities and fitness across diverse environmental conditions shaping its life style. Despite the crucial importance of the pan genome, its assessment among large number of strains remains sparse and systematic studies still needed. Extensive knowledge of the pan genome is needed to translate genome sequencing efforts into developing more efficient biocontrol agents and bio-fertilizers. In this study, a genome survey of B. velezensis allowed us to (a) highlight B. velezensis species boundaries and show that Bacillus suffers taxonomic imprecision that blurs the debate over species pangenome; (b) identify drivers of their successful acquisition of specific life styles and colonization of new niches; (c) describe strategies they use to promote plant growth and development; (d) reveal the unlocked strain specific orphan secondary metabolite gene clusters (biosynthetic clusters with corresponding metabolites unknown) that product identification is still awaiting to amend our knowledge of their putative role in suppression of pathogens and plant growth promotion, and (e) to describe a dynamic pangenome with a secondary metabolite rich accessory genome.

Genomic and Metabolic Characteristics of the Pathogenicity in Pseudomonas aeruginosa

In recent years, the effectiveness of antimicrobials in the treatment of Pseudomonas aeruginosa infections has gradually decreased. This pathogen can be observed in several clinical cases, such as pneumonia, urinary tract infections, sepsis, in immunocompromised hosts, such as neutropenic cancer, burns, and AIDS patients. Furthermore, Pseudomonas aeruginosa causes diseases in both livestock and pets. The highly flexible and versatile genome of P. aeruginosa allows it to have a high rate of pathogenicity. The numerous secreted virulence factors, resulting from its numerous secretion systems, the multi-resistance to different classes of antibiotics, and the ability to produce biofilms are pathogenicity factors that cause numerous problems in the fight against P. aeruginosa infections and that must be better understood for an effective treatment. Infections by P. aeruginosa represent, therefore, a major health problem and, as resistance genes can be disseminated between the microbiotas associated with humans, animals, and the environment, this issue needs be addressed on the basis of an One Health approach. This review intends to bring together and describe in detail the molecular and metabolic pathways in P. aeruginosa's pathogenesis, to contribute for the development of a more targeted therapy against this pathogen.

Phylogenetic Analysis Reveals Distinct Evolutionary Trajectories of the Fluoroquinolones-Resistant Escherichia coli ST1193 From Fuzhou, China

Escherichia coli (E. coli) ST1193 is an emerging fluoroquinolones-resistant and virulent lineage. Large gaps remain in our understanding of the evolutionary processes and differences of this lineage. Therefore, we used 76 E. coli ST1193 genomes to detect strain-level genetic diversity and phylogeny of this lineage globally. All E. coli ST1193 possessed fimH64, filCH5, and fumC14. There was 94.7% of isolates classified as O-type O75. There was 9.33% of E. coli ST1193 that possessed K5 capsular, while 90.67% of isolates possessed K1 capsular. The core genome analysis revealed that all isolates were divided into two phylogenetic clades (clade A and B). Clade A included 25 non-Chinese E. coli ST1193, and clade B contained all isolates collected from Fuzhou, China, respectively. The results of comparative genomics indicated Indels were identified in 150 clade-specific genes, which were enriched into the biological process and molecular function. Accessory genome phylogenetic tree showed a high degree of correlation between accessory genome clusters and core genome clades. There was significant difference in antibiotic resistance genes (ARGs) [bla _CTX-M-55 , bla _TEM-1 , sul2, tet(B), tet(R), APH(6)-Id, and AAC(3)-IId], virulence factors (cia, neuC, gad, and traT), and plasmid replicon types (IncQ1, Col156, and IncB/O/K/Z) between clade A (non-Chinese isolates) and clade B (Chinese isolates) (p < 0.05). Further analysis of the genetic environments of bla _CTX-M-55 demonstrated that the flanking contexts of bla _CTX-M-55 were diverse. In conclusion, our results reveal the distinct evolutionary trajectories of the spread of E. coli ST1193 in Fuzhou, China and non-China regions. This supports both global transmission and localized lineage expansion of this lineage following specific introductions into a geographic locality.

Antimicrobial Resistance and Genomic Characterization of Six New Sequence Types in Multidrug-Resistant Pseudomonas aeruginosa Clinical Isolates from Pakistan

Pseudomonas aeruginosa (P. aeruginosa) is a major bacterial pathogen associated with a variety of infections with high mortality rates. Most of the clinical P. aeruginosa isolates belong to a limited number of genetic subgroups characterized by multiple housekeeping genes’ sequences (usually 5–7) through the Multi-Locus Sequence Typing (MLST) scheme. The emergence and dissemination of novel multidrug-resistant (MDR) sequence types (ST) in P. aeruginosa pose serious clinical concerns. We performed whole-genome sequencing on a cohort (n = 160) of MDR P. aeruginosa isolates collected from a tertiary care hospital lab in Pakistan and found six isolates belonging to six unique MLST allelic profiles. The genomes were submitted to the PubMLST database and new ST numbers (ST3493, ST3494, ST3472, ST3489, ST3491, and ST3492) were assigned to the respective allele combinations. MLST and core-genome-based phylogenetic analysis confirmed the divergence of these isolates and positioned them in separate branches. Analysis of the resistome of the new STs isolates revealed the presence of genes blaOXA-50, blaPAO, blaPDC, blaVIM-2, aph(3′)-IIb, aac(6′)-II, aac(3)-Id, fosA, catB7, dfrB2, crpP, merP and a number of missense and frame-shift mutations in chromosomal genes conferring resistance to various antipseudomonal antibiotics. The exoS, exoT, pvdE, rhlI, rhlR, lasA, lasB, lasI, and lasR genes were the most prevalent virulence-related genes among the new ST isolates. The different genotypic features revealed the adaptation of these new clones to a variety of infections by various mutations in genes affecting antimicrobial resistance, quorum sensing and biofilm formation. Close monitoring of these antibiotic-resistant pathogens and surveillance mechanisms needs to be adopted to reduce their spread to the healthcare facilities of Pakistan. We believe that these strains can be used as reference strains for future comparative analysis of isolates belonging to the same STs.

Genomic Features Associated with the Degree of Phenotypic Resistance to Carbapenems in Carbapenem-Resistant Klebsiella pneumoniae

Carbapenem-resistant Klebsiella pneumoniae strains cause severe infections that are difficult to treat. The production of carbapenemases such as the K. pneumoniae carbapenemase (KPC) is a common mechanism by which these strains resist killing by the carbapenems. However, the degree of phenotypic carbapenem resistance (MIC) may differ markedly between isolates with similar carbapenemase genes, suggesting that our understanding of the underlying mechanisms of carbapenem resistance remains incomplete. To address this problem, we determined the whole-genome sequences of 166 K. pneumoniae clinical isolates resistant to meropenem, imipenem, or ertapenem. Multiple linear regression analysis of this collection of largely bla_KPC-3-containing sequence type 258 (ST258) isolates indicated that bla_KPC copy number and some outer membrane porin gene mutations were associated with higher MICs to carbapenems. A trend toward higher MICs was also observed with those bla_KPC genes carried by the d isoform of Tn4401. In contrast, ompK37 mutations were associated with lower carbapenem MICs, and extended spectrum β-lactamase genes were not associated with higher or lower MICs in carbapenem-resistant K. pneumoniae. A machine learning approach based on the whole-genome sequences of these isolates did not result in a substantial improvement in prediction of isolates with high or low MICs. These results build upon previous findings suggesting that multiple factors influence the overall carbapenem resistance levels in carbapenem-resistant K. pneumoniae isolates.

IMPORTANCEKlebsiella pneumoniae can cause severe infections in the blood, urinary tract, and lungs. Resistance to carbapenems in K. pneumoniae is an urgent public health threat, since it can make these isolates difficult to treat. While individual contributors to carbapenem resistance in K. pneumoniae have been studied, few reports explore their combined effects in clinical isolates. We sequenced 166 clinical carbapenem-resistant K. pneumoniae isolates to evaluate the contribution of known genes to carbapenem MICs and to try to identify novel genes associated with higher carbapenem MICs. The bla_KPC copy number and some outer membrane porin gene mutations were associated with higher carbapenem MICs. In contrast, mutations in one specific porin, ompK37, were associated with lower carbapenem MICs. Machine learning did not result in a substantial improvement in the prediction of carbapenem resistance nor did it identify novel genes associated with carbapenem resistance. These findings enhance our understanding of the many contributors to carbapenem resistance in K. pneumoniae.

Phylogenetic systematics of Butyrivibrio and Pseudobutyrivibrio genomes illustrate vast taxonomic diversity, open genomes and an abundance of carbohydrate-active enzyme family isoforms

Butyrivibrio and Pseudobutyrivibrio dominate in anaerobic gastrointestinal microbiomes, particularly the rumen, where they play a key role in harvesting dietary energy. Within these genera, five rumen species have been classified (Butyrivibrio fibrisolvens, Butyrivibrio hungatei, Butyrivibrio proteoclasticus, Pseudobutyrivibrio ruminis and Pseudobutyrivibrio xylanivorans) and more recently an additional Butyrivibrio sp. group was added. Given the recent increase in available genomes, we re-investigated the phylogenetic systematics and evolution of Butyrivibrio and Pseudobutyrivibrio. Across 71 genomes, we show using 16S rDNA and 40 gene marker phylogenetic trees that the current six species designations (P. ruminis, P. xylanivorans, B. fibrisolvens, Butyrivibrio sp., B. hungatei and B. proteclasticus) are found. However, pangenome analysis showed vast genomic variation and a high abundance of accessory genes (91.50–99.34 %), compared with core genes (0.66–8.50 %), within these six taxonomic groups, suggesting incorrectly assigned taxonomy. Subsequent pangenome accessory genomes under varying core gene cut-offs (%) and average nucleotide identity (ANI) analysis suggest the existence of 42 species within 32 genera. Pangenome analysis of those that still group within B. fibrisolvens, B. hungatei and P. ruminis, based on revised ANI phylogeny, also showed possession of very open genomes, illustrating the diversity that exists even within these groups. All strains of both Butyrivibrio and Pseudobutyrivibrio also shared a broad range of clusters of orthologous genes (COGs) (870), indicating recent evolution from a common ancestor. We also demonstrate that the carbohydrate-active enzymes (CAZymes) predominantly belong to glycosyl hydrolase (GH)2, 3, 5, 13 and 43, with numerous within family isoforms apparent, likely facilitating metabolic plasticity and resilience under dietary perturbations. This study provides a major advancement in our functional and evolutionary understanding of these important anaerobic bacteria.

Rhizospheric and endophytic Pseudomonas aeruginosa in edible vegetable plants share molecular and metabolic traits with clinical isolates

AIM

Pseudomonas aeruginosa, a leading opportunistic pathogen causing hospital-acquired infections, is also commonly found in agricultural settings. However, there are minimal attempts to examine the molecular and functional attributes shared by agricultural and clinical strains of P. aeruginosa. This study investigates the presence of P. aeruginosa in edible vegetable plants (including salad vegetables) and analyses the evolutionary and metabolic relatedness of the agricultural and clinical strains.

METHODS AND RESULTS

Eighteen rhizospheric and endophytic P. aeruginosa strains were isolated from cucumber, tomato, eggplant, and chili directly from the farms. The identity of these strains was confirmed using biochemical and molecular assays. The genetic and metabolic traits of these plant-associated P. aeruginosa isolates were compared with clinical strains. DNA fingerprinting and 16S rDNA-based phylogenetic analyses revealed that the plant- and human-associated strains are evolutionarily related. Both agricultural and clinical isolates possessed plant-beneficial properties, including mineral solubilization to release essential nutrients (phosphorous, potassium, and zinc), ammonification, and the ability to release extracellular pyocyanin, siderophore, and indole-3 acetic acid.

CONCLUSION

These findings suggest that rhizospheric and endophytic P. aeruginosa strains are genetically and functionally analogous to the clinical isolates. In addition, the genotypic and phenotypic traits do not correlate with plant sources or ecosystems.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study reconfirms that edible plants are the potential source for human and animal transmission of P. aeruginosa.

Diversity in the composition of the accessory genome of Mexican Pseudomonas aeruginosa strains

BACKGROUND

Pseudomonas aeruginosa is an important opportunistic pathogen especially in nosocomial infections due to its easy adaptation to different environments; this characteristic is due to the great genetic diversity that presents its genome. In addition, it is considered a pathogen of critical priority due to the high antimicrobial resistance.

OBJECTIVES

The aim of this study was to characterize the mobile genetic elements present in the chromosome of six Mexican P. aeruginosa strains isolated from adults with pneumonia and children with bacteremia.

METHODS

The genomic DNA of six P. aeruginosa strains were isolated and sequenced using PacBio RS-II platform. They were annotated using Prokaryotic Genome Annotation Pipeline and manually curated and analyzed for the presence of mobile genetic elements, antibiotic resistances genes, efflux pumps and virulence factors using several bioinformatics programs and databases.

RESULTS

The global analysis of the strains chromosomes showed a novel chromosomal rearrangement in two strains, possibly mediated by subsequent recombination and inversion events. They have a high content of mobile genetic elements: 21 genomic islands, four new islets, four different integrative conjugative elements, 28 different prophages, one CRISPR-Cas arrangements, and one class 1 integron. The acquisition of antimicrobials resistance genes into these elements are in concordance with their phenotype of multi-drug resistance.

CONCLUSION

The accessory genome increased the ability of the strains to adapt or survive to the hospital environment, promote genomic plasticity and chromosomal rearrangements, which may affect the expression or functionality of the gene and might influence the clinical outcome, having an impact on the treatment.

Whole-genome analysis of Klebsiella pneumoniae from bovine mastitis milk in the U.S

Dairy cattle mastitis has long been one of the most common and costly diseases in the dairy industry worldwide, due to its significant impact on milk production and animal welfare. Among all mastitis causing bacterial pathogens, Klebsiella pneumoniae causes the largest milk loss. To better understand the genomic features of this population, 180 K. pneumoniae strains isolated from dairy cattle mastitis milk in 11 U.S. states were sequenced. The phylogenetic analysis classified all mastitis-causing K. pneumoniae into two major phylogroups, with exclusive predominance in phylogroup KpI. Analysis of more than 61 sequence types, 51 capsular types and 12 lipopolysaccharide O-antigen types revealed great genomic diversity of this K. pneumoniae population. Approximately 100 gene units in accessory genomes were detected with significantly higher prevalence in bovine mastitis strains, compared to human-sourced or dairy environmental strains. The most notable genes were identified associated with ferric citrate uptake, lactose fermentation and resistance to heavy metals. The acquired antimicrobial resistance genes were identified in sporadic mastitis strains. This comprehensive genomic epidemiological study provides insights for a better understanding of the virulence of mastitis-causing K. pneumoniae strains and may lead to the development of novel diagnostic tools and preventive strategies.

Prevalence of IncFIB Plasmids Found among Salmonella enterica Serovar Schwarzengrund Isolates from Animal Sources in Taiwan Using Whole-Genome Sequencing

Salmonella enterica serovar Schwarzengrund is one of the most frequently isolated Salmonella serotypes responsible for human and poultry infections in Taiwan, and it has raised public health concerns. To better facilitate the understanding of transmission patterns and the dynamics of epidemics, sharing molecular data on pathogen profiles is urgently needed. The objectives of the current study were to determine and establish baseline data of S. enterica serovar Schwarzengrund isolates from 23 epidemiologically unrelated sources from year 2000 to 2018 and examine their phenotypic and genotypic characteristics. Genomic DNA of the Salmonella isolates was extracted and subjected to whole-genome sequencing using an Illumina platform. Results showed that all selected isolates exhibited multidrug resistance, and six of those were resistant to ciprofloxacin phenotypically. Genotypically, these isolates carried genes resistant to aminoglycoside (100%), phenicol (91.3%), β-lactams (69.5%), folate pathway antagonist (100%), tetracycline (82.6%), and fluoroquinolone (4.3%). Moreover, these isolates harbor integrons with five different gene cassettes identified for the first time, which are associated with resistance to trimethoprim, streptomycin, tetracycline, sulfonamide, chloramphenicol, and gentamicin. Furthermore, prevalence of IncFIB plasmid was found among studied isolates, which may increase its ability to colonize the chicken cecum and cause extra-intestinal disease. Salmonella pathogenicity islands SPI-1 to SPI-5, SPI-13, and SPI-14, as well as C63PI locus, were also detected in all isolates. This study demonstrated that a considerable high antimicrobial resistance with high virulence levels of Salmonella were found from animal sources. Sharing data on these pathogen profiles can not only help increase the reproducibility and accessibility of genomic analysis but can also support surveillance and epidemiological investigations for salmonellosis in the region.

Differences in the genome, methylome, and transcriptome do not differentiate isolates of Streptococcus equi subsp. equi from horses with acute clinical signs from isolates of inapparent carriers

Streptococcus equi subsp. equi (SEE) is a host-restricted bacterium that causes the common infectious upper respiratory disease known as strangles in horses. Perpetuation of SEE infection appears attributable to inapparent carrier horses because it neither persists long-term in the environment nor infects other host mammals or vectors, and infection results in short-lived immunity. Whether pathogen factors enable SEE to remain in horses without causing clinical signs remains poorly understood. Thus, our objective was to use next-generation sequencing technologies to characterize the genome, methylome, and transcriptome of isolates of SEE from horses with acute clinical strangles and inapparent carrier horses—including isolates recovered from individual horses sampled repeatedly—to assess pathogen-associated changes that might reflect specific adaptions of SEE to the host that contribute to inapparent carriage. The accessory genome elements and methylome of SEE isolates from Sweden and Pennsylvania revealed no significant or consistent differences between acute clinical and inapparent carrier isolates of SEE. RNA sequencing of SEE isolates from Pennsylvania demonstrated no genes that were differentially expressed between acute clinical and inapparent carrier isolates of SEE. The absence of specific, consistent changes in the accessory genomes, methylomes, and transcriptomes of acute clinical and inapparent carrier isolates of SEE indicates that adaptations of SEE to the host are unlikely to explain the carrier state of SEE. Efforts to understand the carrier state of SEE should instead focus on host factors.

Pan-Mitogenomics Approach Discovers Diversity and Dynamism in the Prominent Brown Rot Fungal Pathogens

Monilinia fructicola and Monilinia laxa species are the most destructive and economically devastating fungal plant pathogens causing brown rot disease on stone and pome fruits worldwide. Mitochondrial genomes (mitogenomes) play critical roles influencing the mechanisms and directions of the evolution of fungal pathogens. The pan-mitogenomics approach predicts core and accessory regions of the mitochondrial genomes and explains the gain or loss of variation within and between species. The present study is a fungal pan-mitogenome of M. fructicola (N = 8) and M. laxa (N = 8) species. The completely sequenced and annotated mitogenomes showed high variability in size within and between the species. The mitogenomes of M. laxa were larger, ranging from 178,351 to 179,780bp, than the mitogenomes of M. fructicola, ranging from 158,607 to 167,838bp. However, size variation within the species showed that M. fructicola isolates were more variable in the size range than M. laxa isolates. All the mitogenomes included conserved mitochondrial genes, as well as variable regions including different mobile introns encoding homing endonucleases or maturase, non-coding introns, and repetitive elements. The linear model analysis supported the hypothesis that the mitogenome size expansion is due to presence of variable (accessory) regions. Gene synteny was mostly conserved among all samples, with the exception for order of the rps3 in the mitogenome of one isolate. The mitogenomes presented AT richness; however, A/T and G/C skew varied among the mitochondrial genes. The purifying selection was detected in almost all the protein-coding genes (PCGs) between the species. However, cytochrome b was the only gene showing a positive selection signal among the total samples. Combined datasets of amino acid sequences of 14 core mitochondrial PCGs and rps3 obtained from this study together with published mitochondrial genome sequences from some other species from Heliotales were used to infer a maximum likelihood (ML) phylogenetic tree. ML tree indicated that both Monilinia species highly diverged from each other as well as some other fungal species from the same order. Mitogenomes harbor much information about the evolution of fungal plant pathogens, which could be useful to predict pathogenic life strategies.

Expanding of ST11 Carbapenemase-Producing Klebsiella pneumoniae Subclones in a Chinese Hospital, Shenzhen, China

Background

ST11 is the most prevalent sequence type of clinical Klebsiella pneumoniae in China.

Methods

We investigated the characteristics of the ST11 subclones using core genome multi-locus sequence typing (cgMLST). Ninety-three carbapenemase-producing K. pneumoniae isolates were collected at Shenzhen People's Hospital. Then, whole-genome sequencing and cgMLST were used to discriminate apparent subclones within the ST11 group.

Results

We analyzed the prevalence and genetic relationships of these subclones. ST11 and K. pneumoniae carbapenemase (KPC-2) were the predominant genotype and carbapenemase, respectively, in the clinical carbapenemase-producing K. pneumoniae strains. cgMLST scheme genotyping divided the ST11 group into two clades across seven complex types (CTs). CT1313 was the most prevalent subclone. The deletion of galF and a high frequency of SNPs in genes associated with the stress- and SOS-responses were found in CT1291 and CT2405 over time, respectively.

Conclusion

Our results indicated that the subclones of the ST11 group had different patterns of prevalence. Highly discriminatory genotyping techniques, such as cgMLST scheme, should be used in further molecular epidemiology investigations.

Pangenome analyses of LuxS-coding genes and enzymatic repertoires in cocoa-related lactic acid bacteria

Lactobacillaceae presents potential for interspecific Quorum Sensing (QS) in spontaneous cocoa fermentation, correlated with high abundance of luxS. Three Brazilian isolates from cocoa fermentation were characterized by Whole Genome Sequencing and luxS gene was surveyed in their genomes, in comparison with public databases. They were classified as Lactiplantibacillus plantarum, Limosilactobacillus fermentum and Pediococcus acidilactici. LuxS genes were conserved in core genomes of the novel isolates, but in some non-cocoa related Lactic Acid Bacteria (LAB) it was accessory and plasmid-borne. The conservation and horizontal acquisition of luxS reinforces that QS is determinant for bacterial adaptation in several environments, especially taking into account the luxS has been correlated with modulation of bacteriocin production, stress tolerance and biofilm formation. Therefore, in this paper, new clade and species-specific primers were designed for future application for screening of luxS gene in LAB to evaluate the adaptive potential to diverse food fermentations.

Large-scale sequencing and comparative analysis of oenological Saccharomyces cerevisiae strains supported by nanopore refinement of key genomes

Saccharomyces cerevisiae has long been part of human activities related to the production of food and wine. The industrial demand for fermented beverages with well-defined and stable characteristics boosted the isolation and selection of strains conferring a distinctive aroma profile to the final product. To uncover variants characterizing oenological strains, the sequencing of 65 new S. cerevisiae isolates, and the comparison with other 503 publicly available genomes were performed. A hybrid approach based on short Illumina and long Oxford Nanopore reads allowed the in-depth investigation of eleven genomes and the identification of putative laterally transferred regions and structural variants. A comparative analysis between clusters of strains belonging to different datasets allowed the identification of novel relevant genetic features including single nucleotide polymorphisms, insertions and structural variants. Detection of oenological single nucleotide variants shed light on the existence of different levels of modulation for the mevalonate pathway relevant for the biosynthesis of aromatic compounds.

Draft Genome Sequences of 27 Northern Maine Clinical Isolates

We report the draft genome sequences of 27 common pathogens collected from a northern Maine hospital in 2017. These were sequenced in order to determine temporal and biogeographical patterns of antibiotic gene distribution. A total of 908 antibiotic resistance genes, 848 insertion sequence elements, and 57 plasmids were identified.

We report the draft genome sequences of 27 common pathogens collected from a northern Maine hospital in 2017. These were sequenced in order to determine temporal and biogeographical patterns of antibiotic gene distribution. A total of 908 antibiotic resistance genes, 848 insertion sequence elements, and 57 plasmids were identified.

Genomic context of the two integrons of ST-111 Pseudomonas aeruginosa AG1: A VIM-2-carrying old-acquaintance and a novel IMP-18-carrying integron

Pseudomonas aeruginosa is an opportunist and versatile organism responsible for infections mainly in immunocompromised hosts. This pathogen has high intrinsic resistance to most antimicrobials. P. aeruginosa AG1 (PaeAG1) is a Costa Rican high-risk ST-111 strain with resistance to multiple antibiotics, including carbapenems, due to the activity of VIM-2 and IMP-18 metallo-β-lactamases (MBLs). These genes are harbored in two class 1 integrons located inone out of the 57 PaeAG1 genomic islands. However, the genomic context associated to these determinants in PaeAG1 and other P. aeruginosa strains is unclear. Thus, we first assessed the transcriptional activity of VIM-2 and IMP-18 genes when exposed to imipenem (a carbapenem) by RT-qPCR. To select related genomes to PaeAG1, we implemented a pan-genome analysis to define and up-date the phylogenetic relationship among complete P. aeruginosa genomes. We also studied the PaeAG1 genomic islands content in the related strains and finally we described the architecture and possible evolutionary steps of the genomic regions around the VIM-2- and IMP-18-carrying integrons. Expression of VIM-2 and IMP-18 genes was demonstrated to be induced after imipenem exposure. In a subsequent comparative genomics analysis with 211 strains, the P. aeruginosa pan-genome revealed that complete genome sequences are able to separate clones by MLST profile, including a clear ST-111 cluster with PaeAG1. The PaeAG1 genomic islands were found to define a diverse presence/absence pattern among related genomes. Finally, landscape reconstruction of genomic regions showed that VIM-2-carrying integron (In59-like) is an old-acquaintance element harbored in the same known region found in other two ST-111 strains. Also, PaeAG1 has an exclusive genomic region containing a novel IMP-18-carrying integron (registered as In1666), with an arrangement never reported before. Altogether, we provide new insights about the genomic determinants associated with the resistance to carbapenems in this high-risk P. aeruginosa using comparative genomics.

Estimation of pathogenic potential of an environmental Pseudomonas aeruginosa isolate using comparative genomics

The isolation and sequencing of new strains of Pseudomonas aeruginosa created an extensive dataset of closed genomes. Many of the publicly available genomes are only used in their original publication while additional in silico information, based on comparison to previously published genomes, is not being explored. In this study, we defined and investigated the genome of the environmental isolate P. aeruginosa KRP1 and compared it to more than 100 publicly available closed P. aeruginosa genomes. By using different genomic island prediction programs, we could identify a total of 17 genomic islands and 8 genomic islets, marking the majority of the accessory genome that covers ~ 12% of the total genome. Based on intra-strain comparisons, we are able to predict the pathogenic potential of this environmental isolate. It shares a substantial amount of genomic information with the highly virulent PSE9 and LESB58 strains. For both of these, the increased virulence has been directly linked to their accessory genome before. Hence, the integrated use of previously published data can help to minimize expensive and time consuming wetlab work to determine the pathogenetic potential.

Identification of New Helicobacter pylori Subpopulations in Native Americans and Mestizos From Peru

Region-specific Helicobacter pylori subpopulations have been identified. It is proposed that the hspAmerind subpopulation is being displaced from the Americans by an hpEurope population following the conquest. Our study aimed to describe the genomes and methylomes of H. pylori isolates from distinct Peruvian communities: 23 strains collected from three groups of Native Americans (Asháninkas [ASHA, n = 9], Shimaas [SHIM, n = 5] from Amazonas, and Punos from the Andean highlands [PUNO, n = 9]) and 9 modern mestizos from Lima (LIM). Closed genomes and DNA modification calls were obtained using SMRT/PacBio sequencing. We performed evolutionary analyses and evaluated genomic/epigenomic differences among strain groups. We also evaluated human genome-wide data from 74 individuals from the selected Native communities (including the 23 H. pylori strains donors) to compare host and bacterial backgrounds. There were varying degrees of hspAmerind ancestry in all strains, ranging from 7% in LIM to 99% in SHIM. We identified three H. pylori subpopulations corresponding to each of the Native groups and a novel hspEuropePeru which evolved in the modern mestizos. The divergence of the indigenous H. pylori strains recapitulated the genetic structure of Native Americans. Phylogenetic profiling showed that Orthogroups in the indigenous strains seem to have evolved differentially toward epigenomic regulation and chromosome maintenance, whereas OGs in the modern mestizo (LIM) seem to have evolved toward virulence and adherence. The prevalence of cagA⁺/vacA s1i1m1 genotype was similar across populations (p = 0.32): 89% in ASHA, 67% in PUNO, 56% in LIM and 40% in SHIM. Both cagA and vacA sequences showed that LIM strains were genetically differentiated (p < 0.001) as compared to indigenous strains. We identified 642 R-M systems with 39% of the associated genes located in the core genome. We found 692 methylation motifs, including 254 population-specific sequences not previously described. In Peru, hspAmerind is not extinct, with traces found even in a heavily admixed mestizo population. Notably, our study identified three new hspAmerind subpopulations, one per Native group; and a new subpopulation among mestizos that we named hspEuropePeru. This subpopulation seems to have more virulence-related elements than hspAmerind. Purifying selection driven by variable host immune response may have shaped the evolution of Peruvian subpopulations, potentially impacting disease outcomes.

Metagenomic and Metatranscriptomic Analyses Revealed Uncultured Bacteroidales Populations as the Dominant Proteolytic Amino Acid Degraders in Anaerobic Digesters

Current understanding of amino acid (AA) degraders in anaerobic digesters is mainly based on cultured species, whereas microorganisms that play important roles in a complex microbial community remain poorly characterized. This study investigated short-term enrichments degrading single AAs using metagenomics and metatranscriptomics. Metagenomic analysis revealed that populations related to cultured AA degraders had an abundance <2.5% of the sequences. In contrast, metagenomic-assembled bins related to uncultured Bacteroidales collectively accounted for >35% of the sequences. Phylogenetic analyses suggested that these Bacteroidales populations represented a yet-to-be characterized family lineage, i.e., Bacteroidetes vadinHA17. The bins possessed the genetic capacity related to protein degradation, including surface adhesion (3–7 genes), secreted peptidase (52–77 genes), and polypeptide-specific transporters (2–5 genes). Furthermore, metatranscriptomics revealed that these Bacteroidales populations expressed the complete metabolic pathways for degrading 16 to 17 types of AAs in enrichments fed with respective substrates. These characteristics were distinct from cultured AA degraders including Acidaminobacter and Peptoclostridium, suggesting the uncultured Bacteroidales were the major protein-hydrolyzing and AA-degrading populations. These uncultured Bacteroidales were further found to be dominant and active in full-scale anaerobic digesters, indicating their important ecological roles in the native habitats. “Candidatus Aminobacteroidaceae” was proposed to represent the previously uncharted family Bacteroidetes vadinHA17.

Analysis of the Molecular Diversity Among Cronobacter Species Isolated From Filth Flies Using Targeted PCR, Pan Genomic DNA Microarray, and Whole Genome Sequencing Analyses

Cronobacter species are opportunistic pathogens capable of causing life-threatening infections in humans, with serious complications arising in neonates, infants, immuno-compromised individuals, and elderly adults. The genus is comprised of seven species: Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite a multiplicity of genomic data for the genus, little is known about likely transmission vectors. Using DNA microarray analysis, in parallel with whole genome sequencing, and targeted PCR analyses, the total gene content of two C. malonaticus, three C. turicensis, and 14 C. sakazaki isolated from various filth flies was assessed. Phylogenetic relatedness among these and other strains obtained during surveillance and outbreak investigations were comparatively assessed. Specifically, microarray analysis (MA) demonstrated its utility to cluster strains according to species-specific and sequence type (ST) phylogenetic relatedness, and that the fly strains clustered among strains obtained from clinical, food and environmental sources from United States, Europe, and Southeast Asia. This combinatorial approach was useful in data mining for virulence factor genes, and phage genes and gene clusters. In addition, results of plasmidotyping were in agreement with the species identity for each strain as determined by species-specific PCR assays, MA, and whole genome sequencing. Microarray and BLAST analyses of Cronobacter fly sequence datasets were corroborative and showed that the presence and absence of virulence factors followed species and ST evolutionary lines even though such genes were orthologous. Additionally, zebrafish infectivity studies showed that these pathotypes were as virulent to zebrafish embryos as other clinical strains. In summary, these findings support a striking phylogeny amongst fly, clinical, and surveillance strains isolated during 2010–2015, suggesting that flies are capable vectors for transmission of virulent Cronobacter spp.; they continue to circulate among United States and European populations, environments, and that this “pattern of circulation” has continued over decades.