The bacterial pangenome as a new tool for analysing pathogenic bacteria

Comparative analysis of Streptococcus suis genomes identifies novel candidate virulence-associated genes in North American isolates

Streptococcus suis is a significant economic and welfare concern in the swine industry. Pan-genome analysis provides an in-silico approach for the discovery of genes involved in pathogenesis in bacterial pathogens. In this study, we performed pan-genome analysis of 208 S. suis isolates classified into the pathogenic, possibly opportunistic, and commensal pathotypes to identify novel candidate virulence-associated genes (VAGs) of S. suis. Using chi-square tests and LASSO regression models, three accessory pan-genes corresponding to S. suis strain P1/7 markers SSU_RS09525, SSU_RS09155, and SSU_RS03100 (>95% identity) were identified as having a significant association with the pathogenic pathotype. The proposed novel SSU_RS09525 + /SSU_RS09155 + /SSU_RS03100 + genotype identified 96% of the pathogenic pathotype strains, suggesting a novel genotyping scheme for predicting the pathogenicity of S. suis isolates in North America. In addition, mobile genetic elements carrying antimicrobial resistance genes (ARGs) and VAGs were identified but did not appear to play a major role in the spread of ARGs and VAGs.

Pathogenicity of Shiga toxin-producing Escherichia coli (STEC) from wildlife: Should we care?

Shiga toxin-producing Escherichia coli (STEC) is one of the most frequent bacterial agents associated with food-borne outbreaks in Europe. In humans, the infection can lead to life-threatening diseases. Domestic and wild animals can harbor STEC, and ruminants are the main STEC reservoirs, although asymptomatic. In the present study we have characterized STEC from wildlife (wild boar (n = 56), red deer (n = 101), red fox (n = 37) and otter (n = 92)). Cultivable STEC (n = 52) were isolated from 17% (n = 49) of the faecal samples. All the isolates were non-O157 STEC encoding stx1 (n = 2; 4%) and/or stx2 genes (n = 51; 98%). Only one strain (2%) isolated from red fox had an antibiotic resistant phenotype. However, when the normalized resistance interpretation of epidemiological cutoffs (NRI ECOFFs) were used, 23% (n = 12) of the strains were non-wildtype to at least one of the antibiotics tested. After analysis by pulsed-field gel electrophoresis (PFGE), 20 strains were selected for whole genome sequencing and belonged to the following serotypes: O27:H30 (n = 15), O146:H28 (n = 2), O146:H21 (n = 1), O178:H19 (n = 1), and O103:H2 (n = 1). In addition to stx, all strains encode several virulence factors such as toxins, adhesins, fimbriae and secretion systems, among others. All sequenced genomes carried several mobile genetic elements (MGEs), such as prophages and/or plasmids. The core genome and the phylogenetic analysis showed close evolutionary relationships between some of the STEC recovered from wildlife and strains of clinical origin, highlighting their pathogenic potential. Overall, our results show the zoonotic potential of STEC strains originating from wildlife, highlighting the importance of monitoring their genomic characteristics following a One Health perspective, in which the health of humans is related to the health of animals, and the environment.

Comparative Genomics of Clostridium baratii Reveals Strain-Level Diversity in Toxin Abundance

Clostridium baratii strains are rare opportunistic pathogens associated with botulism intoxication. They have been isolated from foods, soil and be carried asymptomatically or cause botulism outbreaks. Is not taxonomically related to Clostridium botulinum, but some strains are equipped with BoNT/F7 cluster. Despite their relationship with diseases, our knowledge regarding the genomic features and phylogenetic characteristics is limited. We analyzed the pangenome of C. baratii to understand the diversity and genomic features of this species. We compared existing genomes in public databases, metagenomes, and one newly sequenced strain isolated from an asymptomatic subject. The pangenome was open, indicating it comprises genetically diverse organisms. The core genome contained 28.49% of the total genes of the pangenome. Profiling virulence factors confirmed the presence of phospholipase C in some strains, a toxin capable of disrupting eukaryotic cell membranes. Furthermore, the genomic analysis indicated significant horizontal gene transfer (HGT) events as defined by the presence of prophage genomes. Seven strains were equipped with BoNT/F7 cluster. The active site was conserved in all strains, identifying a missing 7-aa region upstream of the active site in C. baratii genomes. This analysis could be important to advance our knowledge regarding opportunistic clostridia and better understand their contribution to disease.

Single-gene long-read sequencing illuminates Escherichia coli strain dynamics in the human intestinal microbiome

Gut microbiome is of major interest due to its close relationship to health and disease. Bacteria usually vary in gene content, leading to functional variations within species, so resolution higher than species-level methods is needed for ecological and clinical relevance. We design a protocol to identify strains in selected species with high discrimination and in high numbers by amplicon sequencing of the flagellin gene. We apply the protocol to fecal samples from a human diet trial, targeting Escherichia coli. Across the 119 samples from 16 individuals, there are 1,532 amplicon sequence variants (ASVs), but only 32 ASVs are dominant in one or more fecal samples, despite frequent dominant strain turnover. Major strains in an intestine are found to be commonly accompanied by a large number of satellite cells, and many are identified as potential extraintestinal pathogens. The protocol could be used to track epidemics or investigate the intra- or inter-host diversity of pathogens.

Comparative pangenomics: analysis of 12 microbial pathogen pangenomes reveals conserved global structures of genetic and functional diversity

BACKGROUND

With the exponential growth of publicly available genome sequences, pangenome analyses have provided increasingly complete pictures of genetic diversity for many microbial species. However, relatively few studies have scaled beyond single pangenomes to compare global genetic diversity both within and across different species. We present here several methods for "comparative pangenomics" that can be used to contextualize multi-pangenome scale genetic diversity with gene function for multiple species at multiple resolutions: pangenome shape, genes, sequence variants, and positions within variants.

RESULTS

Applied to 12,676 genomes across 12 microbial pathogenic species, we observed several shared resolution-specific patterns of genetic diversity: First, pangenome openness is associated with species' phylogenetic placement. Second, relationships between gene function and frequency are conserved across species, with core genomes enriched for metabolic and ribosomal genes and accessory genomes for trafficking, secretion, and defense-associated genes. Third, genes in core genomes with the highest sequence diversity are functionally diverse. Finally, certain protein domains are consistently mutation enriched across multiple species, especially among aminoacyl-tRNA synthetases where the extent of a domain's mutation enrichment is strongly function-dependent.

CONCLUSIONS

These results illustrate the value of each resolution at uncovering distinct aspects in the relationship between genetic and functional diversity across multiple species. With the continued growth of the number of sequenced genomes, these methods will reveal additional universal patterns of genetic diversity at the pangenome scale.

Comparative Genomics Reveals Insights into Induction of Violacein Biosynthesis and Adaptive Evolution in Janthinobacterium

Violacein has different bioactive properties conferring distinct selective advantages, such as defense from predation and interspecific competition. Adaptation of Janthinobacterium to diverse habitats likely leads to variation in violacein production among phylogenetically closely related species inhabiting different environments, yet genomic mechanisms and the influence of adaptive evolution underpinning violacein biosynthesis in Janthinobacterium are not clear. In this study, we performed genome sequencing, comparative genomic analysis, and phenotypic characterization to investigate genomic factors regulating violacein production in nine Janthinobacterium strains, including a type strain from soil and eight strains we isolated from terrestrial subsurface sediment and groundwater. Results show that although all nine Janthinobacterium strains are phylogenetically closely related and contain genes essential for violacein biosynthesis, they vary in carbon usage and violacein production. Sediment and groundwater strains are weak violacein producers and possess far fewer secondary metabolite biosynthesis genes, indicating genome adaptation compared to soil strains. Further examination suggests that quorum sensing (QS) may play an important role in regulating violacein in Janthinobacterium: the strains exhibiting strong potential in violacein production possess both N-acyl-homoserine lactone (AHL) QS and Janthinobacterium QS (JQS) systems in their genomes, while weaker violacein-producing strains harbor only the JQS system. Preliminary tests of spent media of two Janthinobacterium strains possessing both AHL QS and JQS systems support the potential role of AHLs in inducing violacein production in Janthinobacterium. Overall, results from this study reveal potential genomic mechanisms involved in violacein biosynthesis in Janthinobacterium and provide insights into evolution of Janthinobacterium for adaptation to oligotrophic terrestrial subsurface environment.

IMPORTANCE Phylogenetically closely related bacteria can thrive in diverse environmental habitats due to adaptive evolution. Genomic changes resulting from adaptive evolution lead to variations in cellular function, metabolism, and secondary metabolite biosynthesis. The most well-known secondary metabolite produced by Janthinobacterium is the purple-violet pigment violacein. To date, the mechanisms of induction of violacein biosynthesis in Janthinobacterium is not clear. Comparative genome analysis of closely related Janthinobacterium strains isolated from different environmental habitats not only reveals potential mechanisms involved in induction of violacein production by Janthinobacterium but also provides insights into the survival strategy of Janthinobacterium for adaptation to oligotrophic terrestrial subsurface environment.

Pangenomic Analysis of Dickeya dianthicola Strains Related to the Outbreak of Blackleg and Soft Rot of Potato in the United States

Dickeya dianthicola has caused an outbreak of blackleg and soft rot of potato in the eastern half of the United States since 2015. To investigate genetic diversity of the pathogen, a comparative analysis was conducted on genomes of D. dianthicola strains. Whole genomes of 16 strains from the United States outbreak were assembled and compared with 16 previously sequenced genomes of D. dianthicola isolated from potato or carnation. Among the 32 strains, eight distinct clades were distinguished based on phylogenomic analysis. The outbreak strains were grouped into three clades, with the majority of the strains in clade I. Clade I strains were unique and homogeneous, suggesting a recent incursion of this strain into potato production from alternative hosts or environmental sources. The pangenome of the 32 strains contained 6,693 genes, 3,377 of which were core genes. By screening primary protein subunits associated with virulence from all U.S. strains, we found that many virulence-related gene clusters, such as plant cell wall degrading enzyme genes, flagellar and chemotaxis related genes, two-component regulatory genes, and type I/II/III secretion system genes, were highly conserved but that type IV and type VI secretion system genes varied. The clade I strains encoded two clusters of type IV secretion systems, whereas the clade II and III strains encoded only one cluster. Clade I and II strains encoded one more VgrG/PAAR spike protein than did clade III. Thus, we predicted that the presence of additional virulence-related genes may have enabled the unique clade I strain to become predominant in the U.S. outbreak.

Pan-genome Analyses of 3 Strains of Inonotus obliquus and Prediction of Polysaccharide and Terpenoid Genes

Inonotus obliquus is a rare, edible and medicinal fungus that is widely used as a remedy for various diseases. Its main bioactive substances are polysaccharides and terpenoids. In this study, we characterized and investigated the pan-genome of three strains of I. obliquus. The genome sizes of JL01, HE, and NBRC8681 were 32.04, 29.04, and 31.78 Mb, respectively. There were 6 543 core gene families and 6 197 accessory gene families among the three strains, with 14 polysaccharide-related core gene families and seven accessory gene families. For terpenoids, there were 13 core gene families and 17 accessory gene families. Pan-genome sequencing of I. obliquus has improved our understanding of biological characteristics related to the biosynthesis of polysaccharides and terpenoids at the molecular level, which in turn will enable us to increase the production of polysaccharides and terpenoids by this mushroom.

PangenomeNet: a pan-genome-based network reveals functional modules on antimicrobial resistome for Escherichia coli strains

BACKGROUND

Discerning genes crucial to antimicrobial resistance (AMR) mechanisms is becoming more and more important to accurately and swiftly identify AMR pathogenic strains. Pangenome-wide association studies (e.g. Scoary) identified numerous putative AMR genes. However, only a tiny proportion of the putative resistance genes are annotated by AMR databases or Gene Ontology. In addition, many putative resistance genes are of unknown function (termed hypothetical proteins). An annotation tool is crucially needed in order to reveal the functional organization of the resistome and expand our knowledge of the AMR gene repertoire.

RESULTS

We developed an approach (PangenomeNet) for building co-functional networks from pan-genomes to infer functions for hypothetical genes. Using Escherichia coli as an example, we demonstrated that it is possible to build co-functional network from its pan-genome using co-inheritance, domain-sharing, and protein-protein-interaction information. The investigation of the network revealed that it fits the characteristics of biological networks and can be used for functional inferences. The subgraph consisting of putative meropenem resistance genes consists of clusters of stress response genes and resistance gene acquisition pathways. Resistome subgraphs also demonstrate drug-specific AMR genes such as beta-lactamase, as well as functional roles shared among multiple classes of drugs, mostly in the stress-related pathways.

CONCLUSIONS

By demonstrating the idea of pan-genome-based co-functional network on the E. coli species, we showed that the network can infer functional roles of the genes, including those without functional annotations, and provides holistic views on the putative antimicrobial resistomes. We hope that the pan-genome network idea can help formulate hypothesis for targeted experimental works.

Genome Sequence Analysis of Clostridium chauvoei Strains of European Origin and Evaluation of Typing Options for Outbreak Investigations

Black quarter caused by Clostridium (C.) chauvoei is an important bacterial disease that affects cattle and sheep with high mortality. A comparative genomics analysis of 64 C. chauvoei strains, most of European origin and a few of non-European and unknown origin, was performed. The pangenome analysis showed limited new gene acquisition for the species. The accessory genome involved prophages and genomic islands, with variations in gene composition observed in a few strains. This limited accessory genome may indicate that the species replicates only in the host or that an active CRISPR/Cas system provides immunity to foreign genetic elements. All strains contained a CRISPR type I-B system and it was confirmed that the unique spacer sequences therein can be used to differentiate strains. Homologous recombination events, which may have contributed to the evolution of this pathogen, were less frequent compared to other related species from the genus. Pangenome single nucleotide polymorphism (SNP) based phylogeny and clustering indicate diverse clusters related to geographical origin. Interestingly the identified SNPs were mostly non-synonymous. The study demonstrates the possibility of the existence of polymorphic populations in one host, based on strain variability observed for strains from the same animal and strains from different animals of one outbreak. The study also demonstrates that new outbreak strains are mostly related to earlier outbreak strains from the same farm/region. This indicates the last common ancestor strain from one farm can be crucial to understand the genetic changes and epidemiology occurring at farm level. Known virulence factors for the species were highly conserved among the strains. Genetic elements involved in Nicotinamide adenine dinucleotide (NAD) precursor synthesis (via nadA, nadB, and nadC metabolic pathway) which are known as potential anti-virulence loci are completely absent in C. chauvoei compared to the partial inactivation in C. septicum. A novel core-genome MLST based typing method was compared to sequence typing based on CRISPR spacers to evaluate the usefulness of the methods for outbreak investigations.

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.

Comparative Genomic Analysis of Novel Bifidobacterium longum subsp. longum Strains Reveals Functional Divergence in the Human Gut Microbiota

Bifidobacterium longum subsp. longum is a prevalent group in the human gut microbiome. Its persistence in the intestinal microbial community suggests a close host-microbe relationship according to age. The subspecies adaptations are related to metabolic capabilities and genomic and functional diversity. In this study, 154 genomes from public databases and four new Chilean isolates were genomically compared through an in silico approach to identify genomic divergence in genes associated with carbohydrate consumption and their possible adaptations to different human intestinal niches. The pangenome of the subspecies was open, which correlates with its remarkable ability to colonize several niches. The new genomes homogenously clustered within subspecies longum, as observed in phylogenetic analysis. B. longum SC664 was different at the sequence level but not in its functions. COG analysis revealed that carbohydrate use is variable among longum subspecies. Glycosyl hydrolases participating in human milk oligosaccharide use were found in certain infant and adult genomes. Predictive genomic analysis revealed that B. longum M12 contained an HMO cluster associated with the use of fucosylated HMOs but only endowed with a GH95, being able to grow in 2-fucosyllactose as the sole carbon source. This study identifies novel genomes with distinct adaptations to HMOs and highlights the plasticity of B. longum subsp. longum to colonize the human gut microbiota.

Machine Learning Prediction of Resistance to Subinhibitory Antimicrobial Concentrations from Escherichia coli Genomes

Escherichia coli is an important cause of bacterial infections worldwide, with multidrug-resistant strains incurring substantial costs on human lives. Besides therapeutic concentrations of antimicrobials in health care settings, the presence of subinhibitory antimicrobial residues in the environment and in clinics selects for antimicrobial resistance (AMR), but the underlying genetic repertoire is less well understood. Here, we used machine learning to predict the population doubling time and cell growth yield of 1,407 genetically diverse E. coli strains expanding under exposure to three subinhibitory concentrations of six classes of antimicrobials from single-nucleotide genetic variants, accessory gene variation, and the presence of known AMR genes. We predicted cell growth yields in the held-out test data with an average correlation (Spearman's ρ) of 0.63 (0.36 to 0.81 across concentrations) and cell doubling times with an average correlation of 0.59 (0.32 to 0.92 across concentrations), with moderate increases in sample size unlikely to improve predictions further. This finding points to the remaining missing heritability of growth under antimicrobial exposure being explained by effects that are too rare or weak to be captured unless sample size is dramatically increased, or by effects other than those conferred by the presence of individual single-nucleotide polymorphisms (SNPs) and genes. Predictions based on whole-genome information were generally superior to those based only on known AMR genes and were accurate for AMR resistance at therapeutic concentrations. We pinpointed genes and SNPs determining the predicted growth and thereby recapitulated many known AMR determinants. Finally, we estimated the effect sizes of resistance genes across the entire collection of strains, disclosing the growth effects for known resistance genes in each individual strain. Our results underscore the potential of predictive modeling of growth patterns from genomic data under subinhibitory concentrations of antimicrobials, although the remaining missing heritability poses a challenge for achieving the accuracy and precision required for clinical use. IMPORTANCE Predicting bacterial growth from genome sequences is important for a rapid characterization of strains in clinical diagnostics and to disclose candidate novel targets for anti-infective drugs. Previous studies have dissected the relationship between bacterial growth and genotype in mutant libraries for laboratory strains, yet no study so far has examined the predictive power of genome sequence in natural strains. In this study, we used a high-throughput phenotypic assay to measure the growth of a systematic collection of natural Escherichia coli strains and then employed machine learning models to predict bacterial growth from genomic data under nontherapeutic subinhibitory concentrations of antimicrobials that are common in nonclinical settings. We found a moderate to strong correlation between predicted and actual values for the different collected data sets. Moreover, we observed that the known resistance genes are still effective at sublethal concentrations, pointing to clinical implications of these concentrations.

Comparative genome analysis of Salmonella enterica serovar Gallinarum biovars Pullorum and Gallinarum decodes strain specific genes

Salmonella enterica serovar Gallinarum biovar Pullorum (bvP) and biovar Gallinarum (bvG) are the etiological agents of pullorum disease (PD) and fowl typhoid (FT) respectively, which cause huge economic losses to poultry industry especially in developing countries including India. Vaccination and biosecurity measures are currently being employed to control and reduce the S. Gallinarum infections. High endemicity, poor implementation of hygiene and lack of effective vaccines pose challenges in prevention and control of disease in intensively maintained poultry flocks. Comparative genome analysis unravels similarities and dissimilarities thus facilitating identification of genomic features that aids in pathogenesis, niche adaptation and in tracing of evolutionary history. The present investigation was carried out to assess the genotypic differences amongst S.enterica serovar Gallinarum strains including Indian strain S. Gallinarum Sal40 VTCCBAA614. The comparative genome analysis revealed an open pan-genome consisting of 5091 coding sequence (CDS) with 3270 CDS belonging to core-genome, 1254 CDS to dispensable genome and strain specific genes i.e. singletons ranging from 3 to 102 amongst the analyzed strains. Moreover, the investigated strains exhibited diversity in genomic features such as virulence factors, genomic islands, prophage regions, toxin-antitoxin cassettes, and acquired antimicrobial resistance genes. Core genome identified in the study can give important leads in the direction of design of rapid and reliable diagnostics, and vaccine design for effective infection control as well as eradication. Additionally, the identified genetic differences among the S. enterica serovar Gallinarum strains could be used for bacterial typing, structure based inhibitor development by future experimental investigations on the data generated.

Comparative Genomic Analyses of Flavobacterium psychrophilum Isolates Reveals New Putative Genetic Determinants of Virulence Traits

The fish pathogen Flavobacterium psychrophilum is currently one of the main pathogenic bacteria hampering the productivity of salmonid farming worldwide. Although putative virulence determinants have been identified, the genetic basis for variation in virulence of F. psychrophilum is not fully understood. In this study, we analyzed whole-genome sequences of a collection of 25 F. psychrophilum isolates from Baltic Sea countries and compared genomic information with a previous determination of their virulence in juvenile rainbow trout. The results revealed a conserved population of F. psychrophilum that were consistently present across the Baltic Sea countries, with no clear association between genomic repertoire, phylogenomic, or gene distribution and virulence traits. However, analysis of the entire genome of four F. psychrophilum isolates by hybrid assembly provided an unprecedented resolution for discriminating even highly related isolates. The results showed that isolates with different virulence phenotypes harbored genetic variances on a number of consecutive leucine-rich repeat (LRR) proteins, repetitive motifs in gliding motility-associated protein, and the insertion of transposable elements into intergenic and genic regions. Thus, these findings provide novel insights into the genetic variation of these elements and their putative role in the modulation of F. psychrophilum virulence.

Genome based phylogeny and virulence factor analysis of mastitis causing Escherichia coli isolated from Indian cattle

Mastitis is a highly infectious disease prevalent in dairy cattle and it is majorly caused by Escherichia coli (E. coli). The objective of present study is to investigate the occurrence of virulence genes, antimicrobial susceptibility and comparative analysis of E. coli (IVRI KOL CP4 and CM IVRI KOL-1) isolates from mastitis infected animal. Whole-genome sequencing (WGS) was performed using a PacBio RS II system and de novo assembled using Hierarchical Genome Assembly Process (HGAP3). Bacterial Pan Genome Analysis Pipeline (BPGA) was used for pangenome analysis. A set of 50 E. coli isolates were used for comparative analysis (48 collected from the database and 2 reference sequences). Core genes were further concatenated for phylogenetic analyses. In silico analysis was performed for antibiotic resistance and virulence gene identification. Both of the E. coli isolates carried many resistance genes including, b-lactamase, quinolones, rifampicin, macrolide, aminoglycoside and phenicols resistance. We detected 39 virulence genes in IVRI KOL CP4 and 52 in CM IVRI KOL-1 which include toxins, adhesions, invasins, secretion machineries or iron acquisition system. High prevalence of mastitis strains belongs to phylogroups A, although few isolates were also assigned to phylogenetic groups B1 and B2. In conclusion, the present study reported the presence of genes involved in Adherence, Iron acquisition, secretion system and toxins which shown to be crucial in MPEC pathogenicity. This is the first whole genome analysis of MPEC strains to be carried out in Indian isolate to highlights the spread of resistance and virulence genes in food animals.

Pan-proteome profiling of emerging and re-emerging zoonotic pathogen Orientia tsutsugamushi for getting insight into microbial pathogenesis

With the occurrence and evolution of antibiotic and multidrug resistance in bacteria most of the existing remedies are becoming ineffective. The pan-proteome exploration of the bacterial pathogens helps to identify the wide spectrum therapeutic targets which will be effective against all strains in a species. The current study is focused on the pan-proteome profiling of zoonotic pathogen Orientia tsutsugamushi (Ott) for the identification of potential therapeutic targets. The pan-proteome of Ott is estimated to be extensive in nature that has 1429 protein clusters, out of which 694 were core, 391 were accessory, and 344 were unique. It was revealed that 622 proteins were essential, 222 proteins were virulent factors, and 42 proteins were involved in antibiotic resistance. The potential therapeutic targets were further classified into eleven broad classes among which gene expression and regulation, transport, and metabolism were dominant. The biological interactome analysis of therapeutic targets revealed that an ample amount of interactions were present among the proteins involved in DNA replication, ribosome assembly, cellwall metabolism, cell division, and antimicrobial resistance. The predicted therapeutic targets from the pan-proteome of Ott are involved in various biological processes, virulence, and antibiotic resistance; hence envisioned as potential candidates for drug discovery to combat scrub typhus.

Genomic diversity of Serratia proteamaculans and Serratia liquefaciens predominant in seafood products and spoilage potential analyses

Serratia sp. cause food losses and waste due to spoilage; it is noteworthy that they represent a dominant population in seafood. The main spoilage associated species comprise S. liquefaciens, S. grimesii, S. proteamaculans and S. quinivorans, also known as S. liquefaciens-like strains. These species are difficult to discriminate since classical 16S rRNA gene-based sequences do not possess sufficient resolution. In this study, a phylogeny based on the short-length luxS gene was able to speciate 47 Serratia isolates from seafood, with S. proteamaculans being the main species from fresh salmon and tuna, cold-smoked salmon, and cooked shrimp while S. liquefaciens was only found in cold-smoked salmon. The genome of the first S. proteamaculans strain isolated from the seafood matrix (CD3406 strain) was sequenced. Pangenome analyses of S. proteamaculans and S. liquefaciens indicated high adaptation potential. Biosynthetic pathways involved in antimicrobial compounds production and in the main seafood spoilage compounds were also identified. The genetic equipment highlighted in this study contributed to gain further insights into the predominance of Serratia in seafood products and their capacity to spoil.

Zetaproteobacteria Pan-Genome Reveals Candidate Gene Cluster for Twisted Stalk Biosynthesis and Export

Twisted stalks are morphologically unique bacterial extracellular organo-metallic structures containing Fe(III) oxyhydroxides that are produced by microaerophilic Fe(II)-oxidizers belonging to the Betaproteobacteria and Zetaproteobacteria. Understanding the underlying genetic and physiological mechanisms of stalk formation is of great interest based on their potential as novel biogenic nanomaterials and their relevance as putative biomarkers for microbial Fe(II) oxidation on ancient Earth. Despite the recognition of these special biominerals for over 150 years, the genetic foundation for the stalk phenotype has remained unresolved. Here we present a candidate gene cluster for the biosynthesis and secretion of the stalk organic matrix that we identified with a trait-based analyses of a pan-genome comprising 16 Zetaproteobacteria isolate genomes. The “stalk formation in Zetaproteobacteria” (sfz) cluster comprises six genes (sfz1-sfz6), of which sfz1 and sfz2 were predicted with functions in exopolysaccharide synthesis, regulation, and export, sfz4 and sfz6 with functions in cell wall synthesis manipulation and carbohydrate hydrolysis, and sfz3 and sfz5 with unknown functions. The stalk-forming Betaproteobacteria Ferriphaselus R-1 and OYT-1, as well as dread-forming Zetaproteobacteria Mariprofundus aestuarium CP-5 and Mariprofundus ferrinatatus CP-8 contain distant sfz gene homologs, whereas stalk-less Zetaproteobacteria and Betaproteobacteria lack the entire gene cluster. Our pan-genome analysis further revealed a significant enrichment of clusters of orthologous groups (COGs) across all Zetaproteobacteria isolate genomes that are associated with the regulation of a switch between sessile and motile growth controlled by the intracellular signaling molecule c-di-GMP. Potential interactions between stalk-former unique transcription factor genes, sfz genes, and c-di-GMP point toward a c-di-GMP regulated surface attachment function of stalks during sessile growth.

Bacillus pumilus Group Comparative Genomics: Toward Pangenome Features, Diversity, and Marine Environmental Adaptation

Background

Members of the Bacillus pumilus group (abbreviated as the Bp group) are quite diverse and ubiquitous in marine environments, but little is known about correlation with their terrestrial counterparts. In this study, 16 marine strains that we had isolated before were sequenced and comparative genome analyses were performed with a total of 52 Bp group strains. The analyses included 20 marine isolates (which included the 16 new strains) and 32 terrestrial isolates, and their evolutionary relationships, differentiation, and environmental adaptation.

Results

Phylogenomic analysis revealed that the marine Bp group strains were grouped into three species: B. pumilus, B. altitudinis and B. safensis. All the three share a common ancestor. However, members of B. altitudinis were observed to cluster independently, separating from the other two, thus diverging from the others. Consistent with the universal nature of genes involved in the functioning of the translational machinery, the genes related to translation were enriched in the core genome. Functional genomic analyses revealed that the marine-derived and the terrestrial strains showed differences in certain hypothetical proteins, transcriptional regulators, K⁺ transporter (TrK) and ABC transporters. However, species differences showed the precedence of environmental adaptation discrepancies. In each species, land specific genes were found with possible functions that likely facilitate survival in diverse terrestrial niches, while marine bacteria were enriched with genes of unknown functions and those related to transcription, phage defense, DNA recombination and repair.

Conclusion

Our results indicated that the Bp isolates show distinct genomic features even as they share a common core. The marine and land isolates did not evolve independently; the transition between marine and non-marine habitats might have occurred multiple times. The lineage exhibited a priority effect over the niche in driving their dispersal. Certain intra-species niche specific genes could be related to a strains adaptation to its respective marine or terrestrial environment(s). In summary, this report describes the systematic evolution of 52 Bp group strains and will facilitate future studies toward understanding their ecological role and adaptation to marine and/or terrestrial environments.

Pan and Core Genome Analysis of 183 Mycobacterium tuberculosis Strains Revealed a High Inter-Species Diversity among the Human Adapted Strains

Tuberculosis (TB) is an airborne communicable disease with high morbidity and mortality rates, especially in developing countries. The causal agents of TB belong to the complex Mycobacterium tuberculosis (MTBc), which is composed of different human and animal TB associated species. Some animal associated species have zoonotic potential and add to the burden of TB management. The BCG ("Bacillus Calmette-Guérin") vaccine is widely used for the prevention against TB, but its use is limited in immunocompromised patients and animals due to the adverse effects and disseminated life-threatening complications. In this study, we aimed to carry out a comparative genome analysis between the human adapted species including BCG vaccine strains to identify and pinpoint the conserved genes related to the virulence across all the species, which could add a new value for vaccine development. For this purpose, the sequences of 183 Mycobacterium tuberculosis (MTB) strains were retrieved from the freely available WGS dataset at NCBI. The species included: 168 sensu stricto MTB species with other human MTB complex associated strains: M. tuberculosis var. africanum (3), M. tuberculosis var. bovis (2 draft genomes) and 10 BCG species, which enabled the analysis of core genome which contains the conserved genes and some virulence factor determinants. Further, a phylogenetic tree was constructed including the genomes of human (183); animals MTB adapted strains (6) and the environmental Mycobacterium strain "M. canettii". Our results showed that the core genome consists of 1166 conserved genes among these species, which represents a small portion of the pangenome (7036 genes). The remaining genes in the pangenome (5870) are accessory genes, adding a high inter-species diversity. Further, the core genome includes several virulence-associated genes and this could explain the rare infectiousness potential of some attenuated vaccine strains in some patients. This study reveals that low number of conserved genes in human adapted MTBc species and high inter-species diversity of the pan-genome could be considered for vaccine candidate development.

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.

Genome and Pangenome Analysis of Lactobacillus hilgardii FLUB-A New Strain Isolated from Mead

The production of mead holds great value for the Polish liquor industry, which is why the bacterium that spoils mead has become an object of concern and scientific interest. This article describes, for the first time, Lactobacillus hilgardii FLUB newly isolated from mead, as a mead spoilage bacteria. Whole genome sequencing of L. hilgardii FLUB revealed a 3 Mbp chromosome and five plasmids, which is the largest reported genome of this species. An extensive phylogenetic analysis and digital DNA-DNA hybridization confirmed the membership of the strain in the L. hilgardii species. The genome of L. hilgardii FLUB encodes 3043 genes, 2871 of which are protein coding sequences, 79 code for RNA, and 93 are pseudogenes. L. hilgardii FLUB possesses three clustered regularly interspaced short palindromic repeats (CRISPR), eight genomic islands (44,155 bp to 6345 bp), and three (two intact and one incomplete) prophage regions. For the first time, the characteristics of the genome of this species were described and a pangenomic analysis was performed. The concept of the pangenome was used not only to establish the genetic repertoire of this species, but primarily to highlight the unique characteristics of L. hilgardii FLUB. The core of the genome of L. hilgardii is centered around genes related to the storage and processing of genetic information, as well as to carbohydrate and amino acid metabolism. Strains with such a genetic constitution can effectively adapt to environmental changes. L. hilgardii FLUB is distinguished by an extensive cluster of metabolic genes, arsenic detoxification genes, and unique surface layer proteins. Variants of MRS broth with ethanol (10-20%), glucose (2-25%), and fructose (2-24%) were prepared to test the strain's growth preferences using Bioscreen C and the PYTHON script. L. hilgardii FLUB was found to be more resistant than a reference strain to high concentrations of alcohol (18%) and sugars (25%). It exhibited greater preference for fructose than glucose, which suggests it has a fructophilic nature. Comparative genomic analysis supported by experimental research imitating the conditions of alcoholic beverages confirmed the niche specialization of L. hilgardii FLUB to the mead environment.

The Origin of Niches and Species in the Bacterial World

Niches are spaces for the biological units of selection, from cells to complex communities. In a broad sense, "species" are biological units of individuation. Niches do not exist without individual organisms, and every organism has a niche. We use "niche" in the Hutchinsonian sense as an abstraction of a multidimensional environmental space characterized by a variety of conditions, both biotic and abiotic, whose quantitative ranges determine the positive or negative growth rates of the microbial individual, typically a species, but also parts of the communities of species contained in this space. Microbial organisms ("species") constantly diversify, and such diversification (radiation) depends on the possibility of opening up unexploited or insufficiently exploited niches. Niche exploitation frequently implies "niche construction," as the colonized niche evolves with time, giving rise to new potential subniches, thereby influencing the selection of a series of new variants in the progeny. The evolution of niches and organisms is the result of reciprocal interacting processes that form a single unified process. Centrifugal microbial diversification expands the limits of the species' niches while a centripetal or cohesive process occurs simultaneously, mediated by horizontal gene transfers and recombinatorial events, condensing all of the information recovered during the diversifying specialization into "novel organisms" (possible future species), thereby creating a more complex niche, where the selfishness of the new organism(s) establishes a "homeostatic power" limiting the niche's variation. Once the niche's full carrying capacity has been reached, reproductive isolation occurs, as no foreign organisms can outcompete the established population/community, thereby facilitating speciation. In the case of individualization-speciation of the microbiota, its contribution to the animal' gut structure is a type of "niche construction," the result of crosstalk between the niche (host) and microorganism(s). Lastly, there is a parallelism between the hierarchy of niches and that of microbial individuals. The increasing anthropogenic effects on the biosphere (such as globalization) might reduce the diversity of niches and bacterial individuals, with the potential emergence of highly transmissible multispecialists (which are eventually deleterious) resulting from the homogenization of the microbiosphere, a possibility that should be explored and prevented.

Genetic characterisation of a subset of Campylobacter jejuni isolates from clinical and poultry sources in Ireland

Campylobacter spp. is a significant and prevalent public health hazard globally. Campylobacter jejuni is the most frequently recovered species from human cases and poultry are considered the most important reservoir for its transmission to humans. In this study, 30 Campylobacter jejuni isolates were selected from clinical (n = 15) and broiler (n = 15) sources from a larger cohort, based on source, virulence, and antimicrobial resistance profiles. The objective of this study was to further characterise the genomes of these isolates including MLST types, population structure, pan-genome, as well as virulence and antimicrobial resistance determinants. A total of 18 sequence types and 12 clonal complexes were identified. The most common clonal complex was ST-45, which was found in both clinical and broiler samples. We characterised the biological functions that were associated with the core and accessory genomes of the isolates in this study. No significant difference in the prevalence of virulence or antimicrobial resistance determinants was observed between clinical and broiler isolates, although genes associated with severe illness such as neuABC, wlaN and cstIII were only detected in clinical isolates. The ubiquity of virulence factors associated with motility, invasion and cytolethal distending toxin (CDT) synthesis in both clinical and broiler C. jejuni genomes and genetic similarities between groups of broiler and clinical C. jejuni reaffirm that C. jejuni from poultry remains a significant threat to public health.

Pangenomics of the Symbiotic Rhizobiales. Core and Accessory Functions Across a Group Endowed with High Levels of Genomic Plasticity

Pangenome analyses reveal major clues on evolutionary instances and critical genome core conservation. The order Rhizobiales encompasses several families with rather disparate ecological attitudes. Among them, Rhizobiaceae, Bradyrhizobiaceae, Phyllobacteriacreae and Xanthobacteriaceae, include members proficient in mutualistic symbioses with plants based on the bacterial conversion of N₂ into ammonia (nitrogen-fixation). The pangenome of 12 nitrogen-fixing plant symbionts of the Rhizobiales was analyzed yielding total 37,364 loci, with a core genome constituting 700 genes. The percentage of core genes averaged 10.2% over single genomes, and between 5% to 7% were found to be plasmid-associated. The comparison between a representative reference genome and the core genome subset, showed the core genome highly enriched in genes for macromolecule metabolism, ribosomal constituents and overall translation machinery, while membrane/periplasm-associated genes, and transport domains resulted under-represented. The analysis of protein functions revealed that between 1.7% and 4.9% of core proteins could putatively have different functions.

Phenotype and multi-omics comparison of Staphylococcus and Streptococcus uncovers pathogenic traits and predicts zoonotic potential

BACKGROUND

Staphylococcus and Streptococcus species can cause many different diseases, ranging from mild skin infections to life-threatening necrotizing fasciitis. Both genera consist of commensal species that colonize the skin and nose of humans and animals, and of which some can display a pathogenic phenotype.

RESULTS

We compared 235 Staphylococcus and 315 Streptococcus genomes based on their protein domain content. We show the relationships between protein persistence and essentiality by integrating essentiality predictions from two metabolic models and essentiality measurements from six large-scale transposon mutagenesis experiments. We identified clusters of strains within species based on proteins associated to similar biological processes. We built Random Forest classifiers that predicted the zoonotic potential. Furthermore, we identified shared attributes between of Staphylococcus aureus and Streptococcus pyogenes that allow them to cause necrotizing fasciitis.

CONCLUSIONS

Differences observed in clustering of strains based on functional groups of proteins correlate with phenotypes such as host tropism, capability to infect multiple hosts and drug resistance. Our method provides a solid basis towards large-scale prediction of phenotypes based on genomic information.

Comprehensive genomic analysis reveals virulence factors and antibiotic resistance genes in Pantoea agglomerans KM1, a potential opportunistic pathogen

Pantoea agglomerans is a Gram-negative facultative anaerobic bacillus causing a wide range of opportunistic infections in humans including septicemia, pneumonia, septic arthritis, wound infections and meningitis. To date, the determinants of virulence, antibiotic resistance, metabolic features conferring survival and host-associated pathogenic potential of this bacterium remain largely underexplored. In this study, we sequenced and assembled the whole-genome of P. agglomerans KM1 isolated from kimchi in South Korea. The genome contained one circular chromosome of 4,039,945 bp, 3 mega plasmids, and 2 prophages. The phage-derived genes encoded integrase, lysozyme and terminase. Six CRISPR loci were identified within the bacterial chromosome. Further in-depth analysis showed that the genome contained 13 antibiotic resistance genes conferring resistance to clinically important antibiotics such as penicillin G, bacitracin, rifampicin, vancomycin, and fosfomycin. Genes involved in adaptations to environmental stress were also identified which included factors providing resistance to osmotic lysis, oxidative stress, as well as heat and cold shock. The genomic analysis of virulence factors led to identification of a type VI secretion system, hemolysin, filamentous hemagglutinin, and genes involved in iron uptake and sequestration. Finally, the data provided here show that, the KM1 isolate exerted strong immunostimulatory properties on RAW 264.7 macrophages in vitro. Stimulated cells produced Nitric Oxide (NO) and pro-inflammatory cytokines TNF-α, IL-6 and the anti-inflammatory cytokine IL-10. The upstream signaling for production of TNF-α, IL-6, IL-10, and NO depended on TLR4 and TLR1/2. While production of TNF-α, IL-6 and NO involved solely activation of the NF-κB, IL-10 secretion was largely dependent on NF-κB and to a lesser extent on MAPK Kinases. Taken together, the analysis of the whole-genome and immunostimulatory properties provided in-depth characterization of the P. agglomerans KM1 isolate shedding a new light on determinants of virulence that drive its interactions with the environment, other microorganisms and eukaryotic hosts

Comparative Genomics, from the Annotated Genome to Valuable Biological Information: A Case Study

High availability of fast, cheap, and high-throughput next generation sequencing techniques resulted in acquisition of numerous de novo sequenced and assembled bacterial genomes. It rapidly became clear that digging out useful biological information from such a huge amount of data presents a considerable challenge. In this chapter we share our experience with utilization of several handy open source comparative genomic tools. All of them were applied in the studies focused on revealing inter- and intraspecies variation in pectinolytic plant pathogenic bacteria classified to Dickeya solani and Pectobacterium parmentieri. As the described software performed well on the species within the Pectobacteriaceae family, it presumably may be readily utilized on some closely related taxa from the Enterobacteriaceae family. First of all, implementation of various annotation software is discussed and compared. Then, tools computing whole genome comparisons including generation of circular juxtapositions of multiple sequences, revealing the order of synteny blocks or calculation of ANI or Tetra values are presented. Besides, web servers intended either for functional annotation of the genes of interest or for detection of genomic islands, plasmids, prophages, CRISPR/Cas are described. Last but not least, utilization of the software designed for pangenome studies and the further downstream analyses is explained. The presented work not only summarizes broad possibilities assured by the comparative genomic approach but also provides a user-friendly guide that might be easily followed by nonbioinformaticians interested in undertaking similar studies.

Population Biology and Comparative Genomics of Campylobacter Species

The zoonotic pathogen Campylobacter is the leading cause for bacterial foodborne infections in humans. Campylobacters are most commonly transmitted via the consumption of undercooked poultry meat or raw milk products. The decreasing costs of whole genome sequencing enabled large genome-based analyses of the evolution and population structure of this pathogen, as well as the development of novel high-throughput molecular typing methods. Here, we review the evolutionary development and the population diversity of the two most clinically relevant Campylobacter species; C. jejuni and C. coli. The state-of-the-art phylogenetic studies showed clustering of C. jejuni lineages into host specialists and generalists with coexisting lifestyles in chicken and livestock-associated hosts, as well as the separation of C. coli isolates of riparian origin (waterfowl, water) from C. coli isolated from clinical and farm-related samples. We will give an overview of recombination between both species and the potential impact of horizontal gene transfer on host adaptation in Campylobacter. Additionally, this review briefly places the current knowledge of the population structure of other Campylobacter species such as C. lari, C. concisus and C. upsaliensis into perspective. We also provide an overview of how molecular typing methods such as multilocus sequence typing (MLST) and whole genome MLST have been used to detect and trace Campylobacter outbreaks along the food chain.

Building genomes to understand biology

Genetic manipulation is one of the central strategies that biologists use to investigate the molecular underpinnings of life and its diversity. Thus, advances in genetic manipulation usually lead to a deeper understanding of biological systems. During the last decade, the construction of chromosomes, known as synthetic genomics, has emerged as a novel approach to genetic manipulation. By facilitating complex modifications to chromosome content and structure, synthetic genomics opens new opportunities for studying biology through genetic manipulation. Here, we discuss different classes of genetic manipulation that are enabled by synthetic genomics, as well as biological problems they each can help solve.

Reconstructing organisms in silico: genome-scale models and their emerging applications

Escherichia coli is considered to be the best-known microorganism given the large number of published studies detailing its genes, its genome and the biochemical functions of its molecular components. This vast literature has been systematically assembled into a reconstruction of the biochemical reaction networks that underlie E. coli's functions, a process which is now being applied to an increasing number of microorganisms. Genome-scale reconstructed networks are organized and systematized knowledge bases that have multiple uses, including conversion into computational models that interpret and predict phenotypic states and the consequences of environmental and genetic perturbations. These genome-scale models (GEMs) now enable us to develop pan-genome analyses that provide mechanistic insights, detail the selection pressures on proteome allocation and address stress phenotypes. In this Review, we first discuss the overall development of GEMs and their applications. Next, we review the evolution of the most complete GEM that has been developed to date: the E. coli GEM. Finally, we explore three emerging areas in genome-scale modelling of microbial phenotypes: collections of strain-specific models, metabolic and macromolecular expression models, and simulation of stress responses.

New Insights Into DAEC and EAEC Pathogenesis and Phylogeny

Diarrheagenic E. coli can be separated into six distinct pathotypes, with enteroaggregative (EAEC) and diffusely-adherent E. coli (DAEC) among the least characterized. To gain additional insights into these two pathotypes we performed whole genome sequencing of ten DAEC, nine EAEC strains, isolated from Mexican children with diarrhea, and one EAEC plus one commensal E. coli strains isolated from an adult with diarrhea and a healthy child, respectively. These genome sequences were compared to 85 E. coli genomes available in public databases. The EAEC and DAEC strains segregated into multiple different clades; however, six clades were heavily or exclusively comprised of EAEC and DAEC strains, suggesting a phylogenetic relationship between these two pathotypes. EAEC strains harbored the typical virulence factors under control of the activator AggR, but also several toxins, bacteriocins, and other virulence factors. DAEC strains harbored several iron-scavenging systems, toxins, adhesins, and complement resistance or Immune system evasion factors that suggest a pathogenic paradigm for this poorly understood pathotype. Several virulence factors for both EAEC and DAEC were associated with clinical presentations, not only suggesting the importance of these factors, but also potentially indicating opportunities for intervention. Our studies provide new insights into two distinct but related diarrheagenic organisms.

Genomic Characteristics of Invasive Streptococcus pneumoniae Serotype 1 in New Caledonia Prior to the Introduction of PCV13

Streptococcus pneumoniae serotype 1 is a common cause of global invasive pneumococcal disease. In New Caledonia, serotype 1 is the most prevalent serotype and led to two major outbreaks reported in the 2000s. The pneumococcal conjugate vaccine 13 (PCV13) was introduced into the vaccination routine, intending to prevent the expansion of serotype 1 in New Caledonia. Aiming to provide a baseline for monitoring the post-PCV13 changes, we performed a whole-genome sequence analysis on 67 serotype 1 isolates collected prior to the PCV13 introduction. To highlight the S. pneumoniae serotype 1 population structure, we performed a multilocus sequence typing (MLST) analysis revealing that NC serotype 1 consisted of 2 sequence types: ST3717 and the highly dominant ST306. Both sequence types harbored the same resistance genes to beta-lactams, macrolide, streptogramin B, fluoroquinolone, and lincosamide antibiotics. We have also identified 36 virulence genes that were ubiquitous to all the isolates. Among these virulence genes, the pneumolysin sequence presented an allelic profile associated with disease outbreaks and reduced hemolytic activity. Moreover, recombination hotspots were identified in 4 virulence genes and more notably in the cps locus (cps2L), potentially leading to capsular switching, a major mechanism of the emergence of nonvaccine types. In summary, this study represents the first overview of the genomic characteristics of S. pneumoniae serotype 1 in New Caledonia prior to the introduction of PCV13. This preliminary description represents a baseline to assess the impact of PCV13 on serotype 1 population structure and genomic diversity.

Phylogenomic Analysis Supports Two Possible Origins for Latin American Strains of Vibrio parahaemolyticus Associated with Acute Hepatopancreatic Necrosis Disease (AHPND)

Acute hepatopancreatic necrosis disease (AHPND) is a severe disease affecting recently stocked cultured shrimps. The disease is mainly caused by V. parahaemolyticus that harbors the pVA1 plasmid; this plasmid contains the pirA and pirB genes, which encode a delta-endotoxin. AHPND originated in China in 2009 and has since spread to several other Asian countries and recently to Latin America (2013). Many Asian strains have been sequenced, and their sequences are publicly accessible in scientific databases, but only four strains from Latin America have been reported. In this study, we analyzed nine pVA1-harboring V. parahaemolyticus sequences from strains isolated in Mexico along with the 38 previously available pVA1-harboring V. parahaemolyticus sequences and the reference strain RIMD 2210633. The studied sequences were clustered into three phylogenetic clades (Latin American, Malaysian, and Cosmopolitan) through pangenomic and phylogenomic analysis. The nucleotide sequence alignment of the pVA1 plasmids harbored by the Asian and Latin American strains confirmed that the main structural difference in the plasmid between the Asian and Latin American strains is the absence of the Tn3 transposon in the Asian strains; in addition, some deletions in the pirAB region were found in two of the Latin American strains. Our study represents the most robust and inclusive phylogenomic analysis of pVA1-harboring V. parahaemolyticus conducted to date and provides insight into the epidemiology of AHPND. In addition, this study highlights that disease diagnosis through the detection of the pirA and pirB genes is an inadequate approach due to the instability of these genes.

Comparative Genomic Study of Lactobacillus jensenii and the Newly Defined Lactobacillus mulieris Species Identifies Species-Specific Functionality

Lactobacillus species play a key role in the health of the urinary tract. For instance, Lactobacillus crispatus and L. jensenii have been found to inhibit uropathogenic Escherichia coli growth. While L. crispatus is typically found only within the microbiota of women without lower urinary tract symptoms (LUTS), L. jensenii has been found in the microbiota of women both with and without LUTS. With the recent introduction of the new species Lactobacillus mulieris, several strains of L. jensenii were reclassified as L. mulieris based upon gene marker and average nucleotide identity. We took a phylogenomic and comparative genomic approach to ascertain the genetic determinants of these two species. Looking at a larger data set, we identified additional L. mulieris strains, including one distinct from other members of the species—L. mulieris UMB7784. Furthermore, we identified unique loci in each species that may have clinical implications.

Lactobacilli are dominant members of the “healthy” female urogenital microbiota. One of these species, Lactobacillus jensenii, is routinely identified in the urinary microbiota of women both with and without urinary tract symptoms. In March 2020, the new bacterial species Lactobacillus mulieris was introduced, and phylogenetic and average nucleotide identity analysis identified eight L. jensenii strains that should be classified as members of the L. mulieris species. This prompted our phylogenomic study of all publicly available L. jensenii and L. mulieris genome sequences. While there is little variation in the 16S rRNA gene sequences, the core genome shows a clear distinction between genomes of the two species. We find eight additional strains of the species L. mulieris among these genomes. Furthermore, one strain, currently classified as L. mulieris UMB7784, is distinct from both L. jensenii and L. mulieris strains. As part of our comparative genomic study, we also investigated the genetic content that distinguishes these two species. Unique to the L. jensenii genomes are several genes related to catabolism of disaccharides. In contrast, L. mulieris genomes encode several cell surface and secreted proteins that are not found within the L. jensenii genomes. These L. jensenii-specific and L. mulieris-specific loci provide insight into phenotypic differences of these two species.

IMPORTANCELactobacillus species play a key role in the health of the urinary tract. For instance, Lactobacillus crispatus and L. jensenii have been found to inhibit uropathogenic Escherichia coli growth. While L. crispatus is typically found only within the microbiota of women without lower urinary tract symptoms (LUTS), L. jensenii has been found in the microbiota of women both with and without LUTS. With the recent introduction of the new species Lactobacillus mulieris, several strains of L. jensenii were reclassified as L. mulieris based upon gene marker and average nucleotide identity. We took a phylogenomic and comparative genomic approach to ascertain the genetic determinants of these two species. Looking at a larger data set, we identified additional L. mulieris strains, including one distinct from other members of the species—L. mulieris UMB7784. Furthermore, we identified unique loci in each species that may have clinical implications.

Comparative Genomics of Microbacterium Species to Reveal Diversity, Potential for Secondary Metabolites and Heavy Metal Resistance

Microbacterium species have been isolated from a wide range of hosts and environments, including heavy metal-contaminated sites. Here, we present a comprehensive analysis on the phylogenetic distribution and the genetic potential of 70 Microbacterium belonging to 20 different species isolated from heavy metal-contaminated and non-contaminated sites with particular attention to secondary metabolites gene clusters. The analyzed Microbacterium species are divided in three main functional clades. They share a small core genome (331 gene families covering basic functions) pointing to high genetic diversity. The most common secondary metabolite gene clusters encode pathways for the production of terpenoids, type III polyketide synthases and non-ribosomal peptide synthetases, potentially responsible of the synthesis of siderophore-like compounds. In vitro tests showed that many Microbacterium strains produce siderophores, ACC deaminase, auxins (IAA) and are able to solubilize phosphate. Microbacterium isolates from heavy metal contaminated sites are on average more resistant to heavy metals and harbor more genes related to metal homeostasis (e.g., metalloregulators). On the other hand, the ability to increase the metal mobility in a contaminated soil through the secretion of specific molecules seems to be widespread among all. Despite the widespread capacity of strains to mobilize several metals, plants inoculated with selected Microbacterium isolates showed only slightly increased iron concentrations, whereas concentrations of zinc, cadmium and lead were decreased.

The rise and fall of genes: origins and functions of plant pathogen pangenomes

Plant pathogens can rapidly overcome resistance of their hosts by mutating key pathogenicity genes encoding for effectors. Pathogen adaptation is fuelled by extensive genetic variability in populations and different strains may not share the same set of genes. Recently, such an intra-specific variation in gene content became formalized as pangenomes distinguishing core genes (i.e. shared) and accessory genes (i.e. lineage or strain-specific). Across pathogens species, key effectors tend to be part of the rapidly evolving accessory genome. Here, we show how the construction and analysis of pathogen pangenomes provide deep insights into the dynamic host adaptation process. We also discuss how pangenomes should ideally be built and how geography, niche and lifestyle likely determine pangenome sizes.

Genomic Variation and Arsenic Tolerance Emerged as Niche Specific Adaptations by Different Exiguobacterium Strains Isolated From the Extreme Salar de Huasco Environment in Chilean - Altiplano

Polyextremophilic bacteria can thrive in environments with multiple stressors such as the Salar de Huasco (SH). Microbial communities in SH are exposed to low atmospheric pressure, high UV radiation, wide temperature ranges, salinity gradient and the presence of toxic compounds such as arsenic (As). In this work we focus on arsenic stress as one of the main adverse factors in SH and bacteria that belong to the Exiguobacterium genus due to their plasticity and ubiquity. Therefore, our aim was to shed light on the effect of niche conditions pressure (particularly arsenic), on the adaptation and divergence (at genotypic and phenotypic levels) of Exiguobacterium strains from five different SH sites. Also, to capture greater diversity in this genus, we use as outgroup five As(III) sensitive strains isolated from Easter Island (Chile) and The Great Salt Lake (United States). For this, samples were obtained from five different SH sites under an arsenic gradient (9 to 321 mg/kg: sediment) and isolated and sequenced the genomes of 14 Exiguobacterium strains, which had different arsenic tolerance levels. Then, we used comparative genomic analysis to assess the genomic divergence of these strains and their association with phenotypic differences such as arsenic tolerance levels and the ability to resist poly-stress. Phylogenetic analysis showed that SH strains share a common ancestor. Consequently, populations were separated and structured in different SH microenvironments, giving rise to multiple coexisting lineages. Hence, this genotypic variability is also evidenced by the COG (Clusters of Orthologous Groups) composition and the size of their accessory genomes. Interestingly, these observations correlate with physiological traits such as growth patterns, gene expression, and enzyme activity related to arsenic response and/or tolerance. Therefore, Exiguobacterium strains from SH are adapted to physiologically overcome the contrasting environmental conditions, like the arsenic present in their habitat.

Comparative genomics and pangenome-oriented studies reveal high homogeneity of the agronomically relevant enterobacterial plant pathogen Dickeya solani

Background

Dickeya solani is an important plant pathogenic bacterium causing severe losses in European potato production. This species draws a lot of attention due to its remarkable virulence, great devastating potential and easier spread in contrast to other Dickeya spp. In view of a high need for extensive studies on economically important soft rot Pectobacteriaceae, we performed a comparative genomics analysis on D. solani strains to search for genetic foundations that would explain the differences in the observed virulence levels within the D. solani population.

Results

High quality assemblies of 8 de novo sequenced D. solani genomes have been obtained. Whole-sequence comparison, ANIb, ANIm, Tetra and pangenome-oriented analyses performed on these genomes and the sequences of 14 additional strains revealed an exceptionally high level of homogeneity among the studied genetic material of D. solani strains. With the use of 22 genomes, the pangenome of D. solani, comprising 84.7% core, 7.2% accessory and 8.1% unique genes, has been almost completely determined, suggesting the presence of a nearly closed pangenome structure. Attribution of the genes included in the D. solani pangenome fractions to functional COG categories showed that higher percentages of accessory and unique pangenome parts in contrast to the core section are encountered in phage/mobile elements- and transcription- associated groups with the genome of RNS 05.1.2A strain having the most significant impact. Also, the first D. solani large-scale genome-wide phylogeny computed on concatenated core gene alignments is herein reported.

Conclusions

The almost closed status of D. solani pangenome achieved in this work points to the fact that the unique gene pool of this species should no longer expand. Such a feature is characteristic of taxa whose representatives either occupy isolated ecological niches or lack efficient mechanisms for gene exchange and recombination, which seems rational concerning a strictly pathogenic species with clonal population structure. Finally, no obvious correlations between the geographical origin of D. solani strains and their phylogeny were found, which might reflect the specificity of the international seed potato market.

Interest of bacterial pangenome analyses in clinical microbiology

Thanks to the progress and decreasing costs in genome sequencing technologies, more than 250,000 bacterial genomes are currently available in public databases, covering most, if not all, of the major human-associated phylogenetic groups of these microorganisms, pathogenic or not. In addition, for many of them, sequences from several strains of a given species are available, thus enabling to evaluate their genetic diversity and study their evolution. In addition, the significant cost reduction of bacterial whole genome sequencing as well as the rapid increase in the number of available bacterial genomes have prompted the development of pangenomic software tools. The study of bacterial pangenome has many applications in clinical microbiology. It can unveil the pathogenic potential and ability of bacteria to resist antimicrobials as well identify specific sequences and predict antigenic epitopes that allow molecular or serologic assays and vaccines to be designed. Bacterial pangenome constitutes a powerful method for understanding the history of human bacteria and relating these findings to diagnosis in clinical microbiology laboratories in order to optimize patient management.

Impacts of local population history and ecology on the evolution of a globally dispersed pathogen

BACKGROUND

Pathogens with a global distribution face diverse biotic and abiotic conditions across populations. Moreover, the ecological and evolutionary history of each population is unique. Xylella fastidiosa is a xylem-dwelling bacterium infecting multiple plant hosts, often with detrimental effects. As a group, X. fastidiosa is divided into distinct subspecies with allopatric historical distributions and patterns of multiple introductions from numerous source populations. The capacity of X. fastidiosa to successfully colonize and cause disease in naïve plant hosts varies among subspecies, and potentially, among populations. Within Central America (i.e. Costa Rica) two X. fastidiosa subspecies coexist: the native subsp. fastidiosa and the introduced subsp. pauca. Using whole genome sequences, the patterns of gene gain/loss, genomic introgression, and genetic diversity were characterized within Costa Rica and contrasted to other X. fastidiosa populations.

RESULTS

Within Costa Rica, accessory and core genome analyses showed a highly malleable genome with numerous intra- and inter-subspecific gain/loss events. Likewise, variable levels of inter-subspecific introgression were found within and between both coexisting subspecies; nonetheless, the direction of donor/recipient subspecies to the recombinant segments varied. Some strains appeared to recombine more frequently than others; however, no group of genes or gene functions were overrepresented within recombinant segments. Finally, the patterns of genetic diversity of subsp. fastidiosa in Costa Rica were consistent with those of other native populations (i.e. subsp. pauca in Brazil).

CONCLUSIONS

Overall, this study shows the importance of characterizing local evolutionary and ecological history in the context of world-wide pathogen distribution.

A Pangenome Approach for Discerning Species-Unique Gene Markers for Identifications of Streptococcus pneumoniae and Streptococcus pseudopneumoniae

Correct identifications of isolates and strains of the Mitis-Group of the genus Streptococcus are particularly difficult, due to high genetic similarity, resulting from horizontal gene transfer and homologous recombination, and unreliable phenotypic and genotypic biomarkers for differentiating the species. Streptococcus pneumoniae and Streptococcus pseudopneumoniae are the most closely related species of the clade. In this study, publicly-available genome sequences for Streptococcus pneumoniae and S. pseudopneumoniae were analyzed, using a pangenomic approach, to find candidates for species-unique gene markers; ten species-unique genes for S. pneumoniae and nine for S. pseudopneumoniae were identified. These species-unique gene marker candidates were verified by PCR assays for identifying S. pneumoniae and S. pseudopneumoniae strains isolated from clinical samples. All determined species-level unique gene markers for S. pneumoniae were detected in all S. pneumoniae clinical isolates, whereas fewer of the unique S. pseudopneumoniae gene markers were present in more than 95% of the clinical isolates. In parallel, taxonomic identifications of the clinical isolates were confirmed, using conventional optochin sensitivity testing, targeted PCR-detection for the “Xisco” gene, as well as genomic ANIb similarity analyses for the genome sequences of selected strains. Using mass spectrometry-proteomics, species-specific peptide matches were observed for four of the S. pneumoniae gene markers and for three of the S. pseudopneumoniae gene markers. Application of multiple species-level unique biomarkers of S. pneumoniae and S. pseudopneumoniae, is proposed as a protocol for the routine clinical laboratory for improved, reliable differentiation, and identification of these pathogenic and commensal species.

Disentangling the impact of environmental and phylogenetic constraints on prokaryotic within-species diversity

Microbial organisms inhabit virtually all environments and encompass a vast biological diversity. The pangenome concept aims to facilitate an understanding of diversity within defined phylogenetic groups. Hence, pangenomes are increasingly used to characterize the strain diversity of prokaryotic species. To understand the interdependence of pangenome features (such as the number of core and accessory genes) and to study the impact of environmental and phylogenetic constraints on the evolution of conspecific strains, we computed pangenomes for 155 phylogenetically diverse species (from ten phyla) using 7,000 high-quality genomes to each of which the respective habitats were assigned. Species habitat ubiquity was associated with several pangenome features. In particular, core-genome size was more important for ubiquity than accessory genome size. In general, environmental preferences had a stronger impact on pangenome evolution than phylogenetic inertia. Environmental preferences explained up to 49% of the variance for pangenome features, compared with 18% by phylogenetic inertia. This observation was robust when the dataset was extended to 10,100 species (59 phyla). The importance of environmental preferences was further accentuated by convergent evolution of pangenome features in a given habitat type across different phylogenetic clades. For example, the soil environment promotes expansion of pangenome size, while host-associated habitats lead to its reduction. Taken together, we explored the global principles of pangenome evolution, quantified the influence of habitat, and phylogenetic inertia on the evolution of pangenomes and identified criteria governing species ubiquity and habitat specificity.

Genomic analysis of Chromobacterium haemolyticum: insights into the species resistome, virulence determinants and genome plasticity

The increasing number of Chromobacterium haemolyticum human infection reports, especially in tropical regions and connected with environmental sources, resulted in an urge to better describe this species. This study aimed to characterize the C. haemolyticum resistome, virulence determinants and genetic platforms related with genome plasticity. A comparative genomic analysis was conducted between clinical C. haemolyticum genomes publicly available and the genome of an environmental isolate obtained in this study. The pangenome of C. haemolyticum was calculated and a total of 3378 core genes were predicted in its core genome, corresponding to 51.7% of the pangenome. Genetic determinants putatively encoding resistance to beta-lactams, fosfomycin, aminoglycosides and trimethoprim were predicted in all genomes, possibly constituting the intrinsic resistome of this species. In terms of resistance to beta-lactams, 4 genes were predicted encoding beta-lactamases of classes A, C and D. Moreover, the analysis of Chromobacterium genomes and C. haemolyticum environmental isolates reinforced the role of this genus as progenitor of the bla_KPC gene. Putative virulence factors (VFs) were predicted in all genomes, related to adherence, toxins production, colonization and cell invasion. Secretion systems, including type III, were detected. A significant number of transposases and genomic islands were predicted in C. haemolyticum, in some cases above the average reported for Gram-negative bacterial genomes. We conclude that C. haemolyticum strains, including those of environmental origin, present a noteworthy collection of antibiotic resistance genes and VFs. Furthermore, sequences related to gene mobility and genome plasticity suggest high adaptability potential and a possible role as disseminator of antibiotic resistance.

Pan-genomics in the human genome era

Since the early days of the genome era, the scientific community has relied on a single 'reference' genome for each species, which is used as the basis for a wide range of genetic analyses, including studies of variation within and across species. As sequencing costs have dropped, thousands of new genomes have been sequenced, and scientists have come to realize that a single reference genome is inadequate for many purposes. By sampling a diverse set of individuals, one can begin to assemble a pan-genome: a collection of all the DNA sequences that occur in a species. Here we review efforts to create pan-genomes for a range of species, from bacteria to humans, and we further consider the computational methods that have been proposed in order to capture, interpret and compare pan-genome data. As scientists continue to survey and catalogue the genomic variation across human populations and begin to assemble a human pan-genome, these efforts will increase our power to connect variation to human diversity, disease and beyond.

The Genomic Context for the Evolution and Transmission of Community-Associated Staphylococcus aureus ST59 Through the Food Chain

Sequence type 59 (ST59) is a predominant clonal lineage of community-acquired, methicillin-resistant Staphylococcus aureus (CA-MRSA) in Asia. Despite its increasing clinical relevance in China, the evolution and geographic expansion of ST59 has been relatively uncared for. Previous study has shown that ST59 was the predominant clone in food-related MRSA in China. This study compared the genomes of 87 clonal complex (CC) 59 S. aureus isolates sourced from food chain and infection cases to reconstruct the molecular evolution and geographical spread of ST59. Accordingly, three major sub-clades of ST59 were identified and these did not correlate with isolation source or location. Phylogenetic analysis estimated that ST59 in mainland China diverged from a most common recent ancestor around 1974, and most of the cases of cross-country transmission occurred between 1987 and 2000. Notably, two recent events of cross-country transmission through the food chain were observed, the isolates from these events diverged within relatively short time intervals. These isolates also showed high similarity in terms of their core genome, accessory genes, and antibiotic resistance patterns. These findings provide a valuable insight into the potential route of ST59 expansion in China and indicate a need for robust food chain surveillance to prevent the spread of this pathogen.

A comparative genomics approach identifies contact-dependent growth inhibition as a virulence determinant

Emerging evidence suggests the Pseudomonas aeruginosa accessory genome is enriched with uncharacterized virulence genes. Identification and characterization of such genes may reveal novel pathogenic mechanisms used by particularly virulent isolates. Here, we utilized a mouse bacteremia model to quantify the virulence of 100 individual P. aeruginosa bloodstream isolates and performed whole-genome sequencing to identify accessory genomic elements correlated with increased bacterial virulence. From this work, we identified a specific contact-dependent growth inhibition (CDI) system enriched among highly virulent P. aeruginosa isolates. CDI systems contain a large exoprotein (CdiA) with a C-terminal toxin (CT) domain that can vary between different isolates within a species. Prior work has revealed that delivery of a CdiA-CT domain upon direct cell-to-cell contact can inhibit replication of a susceptible target bacterium. Aside from mediating interbacterial competition, we observed our virulence-associated CdiA-CT domain to promote toxicity against mammalian cells in culture and lethality during mouse bacteremia. Structural and functional studies revealed this CdiA-CT domain to have in vitro tRNase activity, and mutations that abrogated this tRNAse activity in vitro also attenuated virulence. Furthermore, CdiA contributed to virulence in mice even in the absence of contact-dependent signaling. Overall, our findings indicate that this P. aeruginosa CDI system functions as both an interbacterial inhibition system and a bacterial virulence factor against a mammalian host. These findings provide an impetus for continued studies into the complex role of CDI systems in P. aeruginosa pathogenesis.

Genomic insight into the origins and evolution of symbiosis genes in Phaseolus vulgaris microsymbionts

Background

Phaseolus vulgaris (common bean) microsymbionts belonging to the bacterial genera Rhizobium, Bradyrhizobium, and Ensifer (Sinorhizobium) have been isolated across the globe. Individual symbiosis genes (e.g., nodC) of these rhizobia can be different within each genus and among distinct genera. Little information is available about the symbiotic structure of indigenous Rhizobium strains nodulating introduced bean plants or the emergence of a symbiotic ability to associate with bean plants in Bradyrhizobium and Ensifer strains. Here, we sequenced the genomes of 29 representative bean microsymbionts (21 Rhizobium, four Ensifer, and four Bradyrhizobium) and compared them with closely related reference strains to estimate the origins of symbiosis genes among these Chinese bean microsymbionts.

Results

Comparative genomics demonstrated horizontal gene transfer exclusively at the plasmid level, leading to expanded diversity of bean-nodulating Rhizobium strains. Analysis of vertically transferred genes uncovered 191 (out of the 2654) single-copy core genes with phylogenies strictly consistent with the taxonomic status of bacterial species, but none were found on symbiosis plasmids. A common symbiotic region was wholly conserved within the Rhizobium genus yet different from those of the other two genera. A single strain of Ensifer and two Bradyrhizobium strains shared similar gene content with soybean microsymbionts in both chromosomes and symbiotic regions.

Conclusions

The 19 native bean Rhizobium microsymbionts were assigned to four defined species and six putative novel species. The symbiosis genes of R. phaseoli, R. sophoriradicis, and R. esperanzae strains that originated from Mexican bean-nodulating strains were possibly introduced alongside bean seeds. R. anhuiense strains displayed distinct host ranges, indicating transition into bean microsymbionts. Among the six putative novel species exclusive to China, horizontal transfer of symbiosis genes suggested symbiosis with other indigenous legumes and loss of originally symbiotic regions or non-symbionts before the introduction of common bean into China. Genome data for Ensifer and Bradyrhizobium strains indicated symbiotic compatibility between microsymbionts of common bean and other hosts such as soybean.