We recently showed that 37/600 (6.2%) invasive infections with group B Streptococcus (GBS) in Toronto, Ontario, Canada, were caused by serotype IV strains. We report a relatively high level of genetic diversity in 37 invasive strains of this emerging GBS serotype. Multilocus sequence typing identified 6 sequence types (STs) that belonged to 3 clonal complexes. Most isolates were ST-459 (19/37, 51%) and ST-452 (11/37, 30%), but we also identified ST-291, ST-3, ST-196, and a novel ST-682. We detected further diversity by performing whole-genome single-nucleotide polymorphism analysis and found evidence of recombination events contributing to variation in some serotype IV GBS strains. We also evaluated antimicrobial drug resistance and found that ST-459 strains were resistant to clindamycin and erythromycin, whereas strains of other STs were, for the most part, susceptible to these antimicrobial drugs.

Clostridium difficile infection (CDI) is intricately linked to the health of the gastrointestinal tract and its indigenous microbiota. In this study, we assessed whether fecal excretion of host DNA is associated with CDI development. Assuming that shedding of epithelial cell increases in the inflamed intestine, we used human DNA excretion as a marker of intestinal insult. Whole-genome shotgun sequencing was employed to quantify host DNA excretion and evaluate bacterial content in fecal samples collected from patients with incipient CDI, hospitalized controls, and healthy subjects. Human DNA excretion was significantly increased in patients admitted to the hospital for a gastrointestinal ailment, as well as prior to an episode of CDI. In multivariable analyses, human read abundance was independently associated with CDI development. Host DNA proportions were negatively correlated with intestinal microbiota diversity. Enterococcus and Escherichia were enriched in patients excreting high quantities of human DNA, while Ruminococcus and Odoribacter were depleted. These findings suggest that intestinal inflammation can occur prior to CDI development and may influence patient susceptibility to CDI. The quantification of human DNA in feces could serve as a simple and noninvasive approach to assess bowel inflammation and identify patients at risk of CDI.

Pseudomonas aeruginosa establishes life-long chronic airway infections in cystic fibrosis (CF) patients. As the disease progresses, P. aeruginosa pathoadaptive variants are distinguished from the initially acquired strain. However, the genetic basis and the biology of host-bacteria interactions leading to a persistent lifestyle of P. aeruginosa are not understood. As a model system to study long term and persistent CF infections, the P. aeruginosa RP73, isolated 16.9 years after the onset of airways colonization from a CF patient, was investigated. Comparisons with strains RP1, isolated at the onset of the colonization, and clonal RP45, isolated 7 years before RP73 were carried out to better characterize genomic evolution of P. aeruginosa in the context of CF pathogenicity.

The O-specific antigen (OSA) in Pseudomonas aeruginosa lipopolysaccharide is highly varied by sugar identity, side chains, and bond between O-repeats. These differences classified P. aeruginosa into 20 distinct serotypes. In the past few decades, O12 has emerged as the predominant serotype in clinical settings and outbreaks. These serotype O12 isolates exhibit high levels of resistance to various classes of antibiotics. Here, we explore how the P. aeruginosa OSA biosynthesis gene clusters evolve in the population by investigating the association between the phylogenetic relationships among 83 P. aeruginosa strains and their serotypes. While most serotypes were closely linked to the core genome phylogeny, we observed horizontal exchange of OSA biosynthesis genes among phylogenetically distinct P. aeruginosa strains. Specifically, we identified a "serotype island" ranging from 62 kb to 185 kb containing the P. aeruginosa O12 OSA gene cluster, an antibiotic resistance determinant (gyrA(C248T)), and other genes that have been transferred between P. aeruginosa strains with distinct core genome architectures. We showed that these genes were likely acquired from an O12 serotype strain that is closely related to P. aeruginosa PA7. Acquisition and recombination of the "serotype island" resulted in displacement of the native OSA gene cluster and expression of the O12 serotype in the recipients. Serotype switching by recombination has apparently occurred multiple times involving bacteria of various genomic backgrounds. In conclusion, serotype switching in combination with acquisition of an antibiotic resistance determinant most likely contributed to the dissemination of the O12 serotype in clinical settings.

Adult invasive disease caused by Group B Streptococcus (GBS) is increasing worldwide. Whole-genome sequencing (WGS) now permits rapid identification of recombination events, a phenomenon that occurs frequently in GBS. Using WGS, we described that strain NGBS375, a capsular serotype V GBS isolate of sequence type (ST)297, has an ST1 genomic background but has acquired approximately 300 kbp of genetic material likely from an ST17 strain. Here, we examined the virulence of this strain in an in vivo model of GBS adult invasive infection. The mosaic ST297 strain showed intermediate virulence, causing significantly less systemic infection and reduced mortality than a more virulent, serotype V ST1 isolate. Bacteremia induced by the ST297 strain was similar to that induced by a serotype III ST17 strain, which was the least virulent under the conditions tested. Yet, under normalized bacteremia levels, the in vivo intrinsic capacity to induce the production of pro-inflammatory cytokines was similar between the ST297 strain and the virulent ST1 strain. Thus, the diminished virulence of the mosaic strain may be due to reduced capacity to disseminate or multiply in blood during a systemic infection which could be mediated by regulatory factors contained in the recombined region.

The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use.

Salmonella enterica is a bacterial species that is a major cause of illness in humans and food-producing animals. S. enterica exhibits considerable inter-serovar diversity, as evidenced by the large number of host adapted serovars that have been identified. The development of methods to assess genome diversity in S. enterica will help to further define the limits of diversity in this foodborne pathogen. Thus, we evaluated a PCR assay, which targets prophage integrase genes, as a rapid method to investigate S. enterica genome diversity. To evaluate the PCR prophage integrase assay, 49 isolates of S. enterica were selected, including 19 clinical isolates from clonal serovars (Enteritidis and Heidelberg) that commonly cause human illness, and 30 isolates from food-associated Salmonella serovars that rarely cause human illness. The number of integrase genes identified by the PCR assay was compared to the number of integrase genes within intact prophages identified by whole genome sequencing and phage finding program PHASTER. The PCR assay identified a total of 147 prophage integrase genes within the 49 S. enterica genomes (79 integrase genes in the food-associated Salmonella isolates, 50 integrase genes in S. Enteritidis, and 18 integrase genes in S. Heidelberg). In comparison, whole genome sequencing and PHASTER identified a total of 75 prophage integrase genes within 102 intact prophages in the 49 S. enterica genomes (44 integrase genes in the food-associated Salmonella isolates, 21 integrase genes in S. Enteritidis, and 9 integrase genes in S. Heidelberg). Collectively, both the PCR assay and PHASTER identified the presence of a large diversity of prophage integrase genes in the food-associated isolates compared to the clinical isolates, thus indicating a high degree of diversity in the food-associated isolates, and confirming the clonal nature of S. Enteritidis and S. Heidelberg. Moreover, PHASTER revealed a diversity of 29 different types of prophages and 23 different integrase genes within the food-associated isolates, but only identified four different phages and integrase genes within clonal isolates of S. Enteritidis and S. Heidelberg. These results demonstrate the potential usefulness of PCR based detection of prophage integrase genes as a rapid indicator of genome diversity in S. enterica.

Leishmania donovani is responsible for visceral leishmaniasis, a neglected and lethal parasitic disease with limited treatment options and no vaccine. The study of L. donovani has been hindered by the lack of a high-quality reference genome and this can impact experimental outcomes including the identification of virulence genes, drug targets and vaccine development. We therefore generated a complete genome assembly by deep sequencing using a combination of second generation (Illumina) and third generation (PacBio) sequencing technologies. Compared to the current L. donovani assembly, the genome assembly reported within resulted in the closure over 2,000 gaps, the extension of several chromosomes up to telomeric repeats and the re-annotation of close to 15% of protein coding genes and the annotation of hundreds of non-coding RNA genes. It was possible to correctly assemble the highly repetitive A2 and Amastin virulence gene clusters. A comparative sequence analysis using the improved reference genome confirmed 70 published and identified 15 novel genomic differences between closely related visceral and atypical cutaneous disease-causing L. donovani strains providing a more complete map of genes associated with virulence and visceral organ tropism. Bioinformatic tools including protein variation effect analyzer and basic local alignment search tool were used to prioritize a list of potential virulence genes based on mutation severity, gene conservation and function. This complete genome assembly and novel information on virulence factors will support the identification of new drug targets and the development of a vaccine for L. donovani.

The structural arrangements of bacterial chromosomes vary widely between closely related species and can result in significant phenotypic outcomes. The appearance of large-scale chromosomal inversions that are symmetric relative to markers for the origin of replication (OriC) has been previously observed; however, the overall prevalence of replication-associated structural rearrangements (RASRs) in bacteria and their causal mechanisms are currently unknown. Here, we systematically identify the locations of RASRs in species with multiple complete-sequenced genomes and investigate potential mediating biological mechanisms. We found that 247 of 313 species contained sequences with at least one large (>50 Kb) inversion in their sequence comparisons, and the aggregated inversion distances away from symmetry were normally distributed with a mean of zero. Many inversions that were offset from dnaA were found to be centered on a different marker for the OriC Instances of flanking repeats provide evidence that breaks formed during the replication process could be repaired to opposing positions. We also found a strong relationship between the later stages of replication and the range in distance variation from symmetry.

Strains of serotype 2 Streptococcus suis are responsible for swine and human infections. Different serotype 2 genetic backgrounds have been defined using multilocus sequence typing (MLST). However, little is known about the genetic diversity within each MLST sequence type (ST). Here, we used whole-genome sequencing to test the hypothesis that S. suis serotype 2 strains of the ST25 lineage are genetically heterogeneous. We evaluated 51 serotype 2 ST25 S. suis strains isolated from diseased pigs and humans in Canada, the United States of America, and Thailand. Whole-genome sequencing revealed numerous large-scale rearrangements in the ST25 genome, compared to the genomes of ST1 and ST28 S. suis strains, which result, among other changes, in disruption of a pilus island locus. We report that recombination and lateral gene transfer contribute to ST25 genetic diversity. Phylogenetic analysis identified two main and distinct Thai and North American clades grouping most strains investigated. These clades also possessed distinct patterns of antimicrobial resistance genes, which correlated with acquisition of different integrative and conjugative elements (ICEs). Some of these ICEs were found to be integrated at a recombination hot spot, previously identified as the site of integration of the 89K pathogenicity island in serotype 2 ST7 S. suis strains. Our results highlight the limitations of MLST for phylogenetic analysis of S. suis, and the importance of lateral gene transfer and recombination as drivers of diversity in this swine pathogen and zoonotic agent.

We describe a genome reference of the African green monkey or vervet (Chlorocebus aethiops). This member of the Old World monkey (OWM) superfamily is uniquely valuable for genetic investigations of simian immunodeficiency virus (SIV), for which it is the most abundant natural host species, and of a wide range of health-related phenotypes assessed in Caribbean vervets (C. a. sabaeus), whose numbers have expanded dramatically since Europeans introduced small numbers of their ancestors from West Africa during the colonial era. We use the reference to characterize the genomic relationship between vervets and other primates, the intra-generic phylogeny of vervet subspecies, and genome-wide structural variations of a pedigreed C. a. sabaeus population. Through comparative analyses with human and rhesus macaque, we characterize at high resolution the unique chromosomal fission events that differentiate the vervets and their close relatives from most other catarrhine primates, in whom karyotype is highly conserved. We also provide a summary of transposable elements and contrast these with the rhesus macaque and human. Analysis of sequenced genomes representing each of the main vervet subspecies supports previously hypothesized relationships between these populations, which range across most of sub-Saharan Africa, while uncovering high levels of genetic diversity within each. Sequence-based analyses of major histocompatibility complex (MHC) polymorphisms reveal extremely low diversity in Caribbean C. a. sabaeus vervets, compared to vervets from putatively ancestral West African regions. In the C. a. sabaeus research population, we discover the first structural variations that are, in some cases, predicted to have a deleterious effect; future studies will determine the phenotypic impact of these variations.

Gene duplication followed by divergence is an important mechanism that leads to molecular innovation. Divergence of paralogous genes can be achieved at functional and regulatory levels. Whereas regulatory divergence at the transcriptional level is well documented, little is known about divergence of posttranslational modifications (PTMs). Protein phosphorylation, one of the most important PTMs, has recently been shown to be an important determinant of the retention of paralogous genes. Here we test whether gains and losses of phosphorylated amino acids after gene duplication may specifically modify the regulation of these duplicated proteins. We show that when phosphosites are lost in one paralog, transitions from phosphorylated serines and threonines are significantly biased toward negatively charged amino acids, which can mimic their phosphorylated status in a constitutive manner. Our analyses support the hypothesis that divergence between paralogs can be generated by a loss of the posttranslational regulatory control on a function rather than by the complete loss of the function itself. Surprisingly, these favoured transitions cannot be reached by single mutational steps, which suggests that the function of a phosphosite needs to be completely abolished before it is restored through substitution by these phosphomimetic residues. We conclude by discussing how gene duplication could facilitate the transitions between phosphorylated and phosphomimetic amino acids.

The draft genome sequence of the ascidian Ciona intestinalis, along with associated gene models, has been a valuable research resource. However, recently accumulated expressed sequence tag (EST)/cDNA data have revealed numerous inconsistencies with the gene models due in part to intrinsic limitations in gene prediction programs and in part to the fragmented nature of the assembly.

Bacterial colonization of the urogenital tract is limited by innate defenses, including the production of antimicrobial peptides (AMPs). Uropathogenic Escherichia coli (UPEC) resist AMP-killing to cause a range of urinary tract infections (UTIs) including asymptomatic bacteriuria, cystitis, pyelonephritis, and sepsis. UPEC strains have high genomic diversity and encode numerous virulence factors that differentiate them from non-UTI-causing strains, including ompT. As OmpT homologs cleave and inactivate AMPs, we hypothesized that UPEC strains from patients with symptomatic UTIs have high OmpT protease activity. Therefore, we measured OmpT activity in 58 clinical E. coli isolates. While heterogeneous OmpT activities were observed, OmpT activity was significantly greater in UPEC strains isolated from patients with symptomatic infections. Unexpectedly, UPEC strains exhibiting the greatest protease activities harbored an additional ompT-like gene called arlC (ompTp). The presence of two OmpT-like proteases in some UPEC isolates led us to compare the substrate specificities of OmpT-like proteases found in E. coli. While all three cleaved AMPs, cleavage efficiency varied on the basis of AMP size and secondary structure. Our findings suggest the presence of ArlC and OmpT in the same UPEC isolate may confer a fitness advantage by expanding the range of target substrates.

Non-typhoidal Salmonella is a leading cause of foodborne illness worldwide. Prompt and accurate identification of the sources of Salmonella responsible for disease outbreaks is crucial to minimize infections and eliminate ongoing sources of contamination. Current subtyping tools including single nucleotide polymorphism (SNP) typing may be inadequate, in some instances, to provide the required discrimination among epidemiologically unrelated Salmonella strains. Prophage genes represent the majority of the accessory genes in bacteria genomes and have potential to be used as high discrimination markers in Salmonella. In this study, the prophage sequence diversity in different Salmonella serovars and genetically related strains was investigated. Using whole genome sequences of 1,760 isolates of S. enterica representing 151 Salmonella serovars and 66 closely related bacteria, prophage sequences were identified from assembled contigs using PHASTER. We detected 154 different prophages in S. enterica genomes. Prophage sequences were highly variable among S. enterica serovars with a median ± interquartile range (IQR) of 5 ± 3 prophage regions per genome. While some prophage sequences were highly conserved among the strains of specific serovars, few regions were lineage specific. Therefore, strains belonging to each serovar could be clustered separately based on their prophage content. Analysis of S. Enteritidis isolates from seven outbreaks generated distinct prophage profiles for each outbreak. Taken altogether, the diversity of the prophage sequences correlates with genome diversity. Prophage repertoires provide an additional marker for differentiating S. enterica subtypes during foodborne outbreaks.

Drug and surgical-based therapies in type 2 diabetes are associated with altered gut microbiota architecture. Here we investigated the role of the gut microbiome in improved glucose homeostasis following bariatric surgery.

Over the past decade, there has been a rising interest in Achromobacter sp., an emerging opportunistic pathogen responsible for nosocomial and cystic fibrosis lung infections. Species of this genus are ubiquitous in the environment, can outcompete resident microbiota, and are resistant to commonly used disinfectants as well as antibiotics. Nevertheless, the Achromobacter genus suffers from difficulties in diagnosis, unresolved taxonomy and limited understanding of how it adapts to the cystic fibrosis lung, not to mention other host environments. The goals of this first genus-wide comparative genomics study were to clarify the taxonomy of this genus and identify genomic features associated with pathogenicity and host adaptation. This was done with a widely applicable approach based on pan-genome analysis. First, using all publicly available genomes, a combination of phylogenetic analysis based on 1,780 conserved genes with average nucleotide identity and accessory genome composition allowed the identification of a largely clinical lineage composed of Achromobacter xylosoxidans, Achromobacter insuavis, Achromobacter dolens, and Achromobacter ruhlandii. Within this lineage, we identified 35 positively selected genes involved in metabolism, regulation and efflux-mediated antibiotic resistance. Second, resistome analysis showed that this clinical lineage carried additional antibiotic resistance genes compared with other isolates. Finally, we identified putative mobile elements that contribute 53% of the genus's resistome and support horizontal gene transfer between Achromobacter and other ecologically similar genera. This study provides strong phylogenetic and pan-genomic bases to motivate further research on Achromobacter, and contributes to the understanding of opportunistic pathogen evolution.

Carbapenemase-producing Enterobacterales, including KPC-2 producers, have become a major clinical problem. During an outbreak in Quebec City, Canada, KPC-2-producing Klebsiella michiganensis and Citrobacter farmeri were isolated from a patient six weeks apart. We determined their complete genome sequences. Both isolates carried nearly identical IncN2 plasmids with blaKPC-2 on a Tn4401b element. Both strains also carried IncP1 plasmids, but that of C. farmeri did not carry a Beta-lactamase gene, whereas that of K. michiganensis carried a second copy of blaKPC-2 on Tn4401b. These results suggest recent plasmid transfer between the two species and a recent transposition event.

Long diagnostic wait times hinder international efforts to address antibiotic resistance in M. tuberculosis. Pathogen whole genome sequencing, coupled with statistical and machine learning models, offers a promising solution. However, generalizability and clinical adoption have been limited by a lack of interpretability, especially in deep learning methods. Here, we present two deep convolutional neural networks that predict antibiotic resistance phenotypes of M. tuberculosis isolates: a multi-drug CNN (MD-CNN), that predicts resistance to 13 antibiotics based on 18 genomic loci, with AUCs 82.6-99.5% and higher sensitivity than state-of-the-art methods; and a set of 13 single-drug CNNs (SD-CNN) with AUCs 80.1-97.1% and higher specificity than the previous state-of-the-art. Using saliency methods to evaluate the contribution of input sequence features to the SD-CNN predictions, we identify 18 sites in the genome not previously associated with resistance. The CNN models permit functional variant discovery, biologically meaningful interpretation, and clinical applicability.

Non-typhoidal Salmonella enterica serovars, associated with different foods including poultry products, are important causes of bacterial gastroenteritis worldwide. The colonization of the chicken gut by S. enterica could result in the contamination of the environment and food chain. The aim of this study was to compare the genomes of 25 S. enterica serovars isolated from broiler chicken farms to assess their intra- and inter-genetic variability, with a focus on virulence and antibiotic resistance characteristics.