A Prolonged Outbreak of Human Adenovirus A31 (HAdV-A31) Infection on a Pediatric Hematopoietic Stem Cell Transplantation Ward with Whole Genome Sequencing Evidence of International Linkages.

A surge in hematopoietic stem cell transplantation (HSCT) human adenovirus A31 (HAdV-A31) infections was initially observed in late 2014/2015 at SickKids (SK) Hospital, Toronto, Canada. In response, enhanced laboratory monitoring for all adenovirus infections was conducted. Positive samples underwent genotyping, viral culture, and, in selected cases, whole-genome sequencing (WGS). HAdV-A31 specimens/DNA obtained from four international pediatric HSCT centers also underwent WGS. During the SK outbreak period (27 October 2014 to 31 October 2018), 17/20 HAdV-A31 isolates formed a distinct clade with 0 to 8 mutations between the closest neighbors. Surveillance before and after the outbreak detected six additional HAdV-A31 HSCT cases; three of the four sequenced cases clustered within the outbreak clade. Two SK outbreak isolates were identical to sequences from two patients in an outbreak in England. Three SK non-outbreak sequences also had high sequence similarity to strains from three international centers. Environmental PCR testing of the HSCT ward showed significant adenovirus contamination. Despite intense infection control efforts, we observed re-occurrence of infection with the outbreak strain. Severe but nonfatal infection was observed more commonly with HAdV-A31 compared to other genotypes, except HAdV-C1. Our findings strongly implicate nosocomial spread of HAdV-A31 over 10 years on a HSCT unit and demonstrate the value of WGS in defining and mapping the outbreak. Close linkages among strains in different countries suggest international dissemination, though the mechanism is undetermined. This large, extended outbreak emphasizes the pre-eminent role of HAdV-A31 in causing intractable pediatric HSCT outbreaks of severe illness worldwide.

Neptune: a bioinformatics tool for rapid discovery of genomic variation in bacterial populations.

The ready availability of vast amounts of genomic sequence data has created the need to rethink comparative genomics algorithms using 'big data' approaches. Neptune is an efficient system for rapidly locating differentially abundant genomic content in bacterial populations using an exact k-mer matching strategy, while accommodating k-mer mismatches. Neptune's loci discovery process identifies sequences that are sufficiently common to a group of target sequences and sufficiently absent from non-targets using probabilistic models. Neptune uses parallel computing to efficiently identify and extract these loci from draft genome assemblies without requiring multiple sequence alignments or other computationally expensive comparative sequence analyses. Tests on simulated and real datasets showed that Neptune rapidly identifies regions that are both sensitive and specific. We demonstrate that this system can identify trait-specific loci from different bacterial lineages. Neptune is broadly applicable for comparative bacterial analyses, yet will particularly benefit pathogenomic applications, owing to efficient and sensitive discovery of differentially abundant genomic loci. The software is available for download at: http://github.com/phac-nml/neptune.

Interactome analysis of longitudinal pharyngeal infection of cynomolgus macaques by group A Streptococcus.

Relatively little is understood about the dynamics of global host-pathogen transcriptome changes that occur during bacterial infection of mucosal surfaces. To test the hypothesis that group A Streptococcus (GAS) infection of the oropharynx provokes a distinct host transcriptome response, we performed genome-wide transcriptome analysis using a nonhuman primate model of experimental pharyngitis. We also identified host and pathogen biological processes and individual host and pathogen gene pairs with correlated patterns of expression, suggesting interaction. For this study, 509 host genes and seven biological pathways were differentially expressed throughout the entire 32-day infection cycle. GAS infection produced an initial widespread significant decrease in expression of many host genes, including those involved in cytokine production, vesicle formation, metabolism, and signal transduction. This repression lasted until day 4, at which time a large increase in expression of host genes was observed, including those involved in protein translation, antigen presentation, and GTP-mediated signaling. The interactome analysis identified 73 host and pathogen gene pairs with correlated expression levels. We discovered significant correlations between transcripts of GAS genes involved in hyaluronic capsule production and host endocytic vesicle formation, GAS GTPases and host fibrinolytic genes, and GAS response to interaction with neutrophils. We also identified a strong signal, suggesting interaction between host gammadelta T cells and genes in the GAS mevalonic acid synthesis pathway responsible for production of isopentenyl-pyrophosphate, a short-chain phospholipid that stimulates these T cells. Taken together, our results are unique in providing a comprehensive understanding of the host-pathogen interactome during mucosal infection by a bacterial pathogen.

whole-genome sequence of livestock-associated st398 methicillin-resistant staphylococcus aureus Isolated from Humans in Canada.

Despite reports of high colonization rates of ST398 livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) among pigs and pig farmers, the incidence of LA-MRSA infection in the general population in Canada appears to be rare in comparison to that in some European countries. In this study, the complete genome sequence of a Canadian representative LA-MRSA isolate (08BA02176) from a human postoperative surgical site infection was acquired and compared to the sequenced genome of an LA-MRSA isolate (S0385) from Europe to identify genetic traits that may explain differences in the success of these particular strains in some locales.

Human infection with a triple-reassortant swine influenza A(H1N1) virus containing the hemagglutinin and neuraminidase genes of seasonal influenza virus.

A reassortant influenza A(H1N1) virus of swine origin distinct from the pandemic H1N1 2009 strain was isolated from 3 patients, all of whom worked at the same large hog operation in Saskatchewan, Canada. The genomic composition of the isolates has not been previously reported, to our knowledge, and was the product of a genetic reassortment between seasonal H1N1 and triple-reassortant influenza virus that emerged in North American swine during the late 1990s. The neuraminidase and hemagglutinin genes of A/Saskatchewan/5350/2009, A/Saskatchewan/5351/2009, and A/Saskatchewan/5131/2009 were derived from human H1N1 virus and were closely related to those of A/Brisbane/59/2007.

Comparative genomics of ST5 and ST30 methicillin-resistant Staphylococcus aureus sequential isolates recovered from paediatric patients with cystic fibrosis.

Staphylococcus aureus chronic airway infection in patients with cystic fibrosis (CF) allows this pathogen to adapt over time in response to different selection pressures. We have previously shown that the main sequence types related to community-acquired methicillin-resistant S. aureus (MRSA) infections in Argentina - ST5 and ST30 - are also frequently isolated from the sputum of patients with CF, but in these patients they usually display multi-drug antimicrobial resistance. In this study, we sequenced the genomes of MRSA from four paediatric CF patients with the goal of identifying mutations among sequential isolates, especially those possibly related to antimicrobial resistance and virulence, which might contribute to the adaptation of the pathogen in the airways of patients with CF. Our results revealed genetic differences in sequential MRSA strains isolated from patients with CF in both their core and accessory genomes. Although the genetic adaptation of S. aureus was distinct in different hosts, we detected independent mutations in thyA, htrA, rpsJ and gyrA - which are known to have crucial roles in S. aureus virulence and antimicrobial resistance - in isolates recovered from multiple patients. Moreover, we identified allelic variants that were detected in all of the isolates recovered after a certain time point; these non-synonymous mutations were in genes associated with antimicrobial resistance, virulence, iron scavenging and oxidative stress resistance. In conclusion, our results provide evidence of genetic variability among sequential MRSA isolates that could be implicated in the adaptation of these strains during chronic CF airway infection.

Livestock-associated methicillin-resistant Staphylococcus aureus sequence type 398 in humans, Canada.

Rates of colonization with livestock-associated methicillin-resistant Staphylococcus aureus (MRSA) sequence type 398 have been high for pigs and pig farmers in Canada, but prevalence rates for the general human population are unknown. In this study, 5 LA-MRSA isolates, 4 of which were obtained from skin and soft tissue infections, were identified from 3,687 tested MRSA isolates from persons in Manitoba and Saskatchewan, Canada. Further molecular characterization determined that these isolates all contained staphylococcal cassette chromosome (SCC) mecV, were negative for Panton-Valentine leukocidin, and were closely related by macrorestriction analysis with the restriction enzyme Cfr91. The complete DNA sequence of the SCCmec region from the isolate showed a novel subtype of SCCmecV harboring clustered regularly interspaced short palindromic repeats and associated genes. Although prevalence of livestock-associated MRSA seems to be low for the general population in Canada, recent emergence of infections resulting from this strain is of public health concern.

Detection of Antimicrobial Resistance Using Proteomics and the Comprehensive Antibiotic Resistance Database: A Case Study.

PURPOSE: Antimicrobial resistance (AMR), especially multidrug resistance, is one of the most serious global threats facing public health. The authors proof-of-concept study assessing the suitability of shotgun proteomics as an additional approach to whole-genome sequencing (WGS) for detecting AMR determinants. EXPERIMENTAL DESIGN: Previously published shotgun proteomics and WGS data on four isolates of Campylobacter jejuni are used to perform AMR detection by searching the Comprehensive Antibiotic Resistance Database, and their detection ability relative to genomics screening and traditional phenotypic testing measured by minimum inhibitory concentration is assessed. RESULTS: Both genomic and proteomic approaches identify the wild-type and variant molecular determinants responsible for resistance to tetracycline and ciprofloxacin, in agreement with phenotypic testing. In contrast, the genomic method identifies the presence of the ß-lactamase gene, blaOXA-61 , in three isolates. However, its corresponding protein product is detected in only a single isolate, consistent with results obtained from phenotypic testing.

The effect of the Salmonella genomic island 1 on in vitro global gene expression in Salmonella enterica serovar Typhimurium LT2.

A Salmonella genomic island 1 (SGI1) isogenic strain pair was constructed using Salmonella enterica serovar Typhimurium LT2 (ST LT2). Real-time quantitative reverse transcriptase PCR revealed detectable mRNA transcripts for all 44 putative ORFs encoded within the SGI1. The highest levels of transcripts observed in SGI1 encoded ORFs were found in genes conferring antibiotic resistance to ampicillin, streptomycin/spectinomycin, and sulphonamides. Abundant mRNA transcripts, relative to gapA, were also noted for one putative regulatory ORF and seven ORFs of unknown function encoded within SGI1, whose products could represent factors contributing to increases in virulence and/or fitness of the organism. DNA microarray analysis revealed the differential expression of known factors that contribute to virulence in many pathogens. Twenty-two chromosomal genes were significantly upregulated in ST LT2 harboring SGI1, which included increased expression of iron and sialic acid utilization genes. Decreased expression was noted for 15 genes in ST LT2 harboring SGI1, including genes involved in chemotaxis and motility. This is the first report examining gene expression within the SGI1, as well as its potential effect on global gene expression, and sets the foundation for future studies involving the effect of SGI1 in other Salmonella spp.

Food Safety in the Age of Next Generation Sequencing, Bioinformatics, and Open Data Access.

Public health labs and food regulatory agencies globally are embracing whole genome sequencing (WGS) as a revolutionary new method that is positioned to replace numerous existing diagnostic and microbial typing technologies with a single new target: the microbial draft genome. The ability to cheaply generate large amounts of microbial genome sequence data, combined with emerging policies of food regulatory and public health institutions making their microbial sequences increasingly available and public, has served to open up the field to the general scientific community. This open data access policy shift has resulted in a proliferation of data being deposited into sequence repositories and of novel bioinformatics software designed to analyze these vast datasets. There also has been a more recent drive for improved data sharing to achieve more effective global surveillance, public health and food safety. Such developments have heightened the need for enhanced analytical systems in order to process and interpret this new type of data in a timely fashion. In this review we outline the emergence of genomics, bioinformatics and open data in the context of food safety. We also survey major efforts to translate genomics and bioinformatics technologies out of the research lab and into routine use in modern food safety labs. We conclude by discussing the challenges and opportunities that remain, including those expected to play a major role in the future of food safety science.

Comparison of Sample Preparation Methods Used for the Next-Generation Sequencing of Mycobacterium tuberculosis.

The advent and widespread application of next-generation sequencing (NGS) technologies to the study of microbial genomes has led to a substantial increase in the number of studies in which whole genome sequencing (WGS) is applied to the analysis of microbial genomic epidemiology. However, microorganisms such as Mycobacterium tuberculosis (MTB) present unique problems for sequencing and downstream analysis based on their unique physiology and the composition of their genomes. In this study, we compare the quality of sequence data generated using the Nextera and TruSeq isolate preparation kits for library construction prior to Illumina sequencing-by-synthesis. Our results confirm that MTB NGS data quality is highly dependent on the purity of the DNA sample submitted for sequencing and its guanine-cytosine content (or GC-content). Our data additionally demonstrate that the choice of library preparation method plays an important role in mitigating downstream sequencing quality issues. Importantly for MTB, the Illumina TruSeq library preparation kit produces more uniform data quality than the Nextera XT method, regardless of the quality of the input DNA. Furthermore, specific genomic sequence motifs are commonly missed by the Nextera XT method, as are regions of especially high GC-content relative to the rest of the MTB genome. As coverage bias is highly undesirable, this study illustrates the importance of appropriate protocol selection when performing NGS studies in order to ensure that sound inferences can be made regarding mycobacterial genomes.

A Study of the Infant Nasal Microbiome Development over the First Year of Life and in Relation to Their Primary Adult Caregivers Using cpn60 Universal Target (UT) as a Phylogenetic Marker.

Whereas the infant gut microbiome is the subject of intense study, relatively little is known regarding the nares microbiome in newborns and during early life. This study aimed to survey the typical composition and diversity of human anterior nare microflora for developing infants over time, and to explore how these correlate to their primary caregivers. Single nare swabs were collected at five time points over a one-year period for each subject from infant-caregiver pairs. Our study comprised of 50 infants (recruited at 2 weeks, post delivery) and their 50 primary caregivers. Applying the chaperonin-60 (cpn60) universal target (UT) amplicon as our molecular barcoding marker to census survey the microbial communities, we longitudinally surveyed infant nares microbiota at 5 time points over the course of the first year of life. The inter- and intra-subject diversity was catalogued and compared, both longitudinally and relative to their adult primary caregivers. Although within-subject variability over time and inter-subject variability were both observed, the assessment detected only one or two predominant genera for individual infant samples, belonging mainly to phyla Actinobacteria, Firmicutes, and Proteobacteria. Consistent with previously observed microbial population dynamics in other body sites, the diversity of nares microflora increased over the first year of life and infants showed differential operational taxonomic units (OTUs) relative to their matched primary caregiver. The collected evidence also support that both temporal and seasonal changes occur with respect to carriage of potentially pathogenic bacteria (PPBs), which may influence host predisposition to infection. This pilot study surveying paired infant/caregiver nare microbiomes provides novel longitudinal diversity information that is pertinent to better understanding nare microbiome development in infants.

Chlamydia trachomatis virulence factor CT135 is stable in vivo but highly polymorphic in vitro.

Chlamydia trachomatis is an important human pathogen causing both ocular and sexually transmitted disease. Recently, we identified CT135 as an important virulence determinant in a mouse infection model. Results from CEL 1 digestion assays and sequencing analyses indicated that CT135 was much more polymorphic in high in vitro passage reference serovars than it was in clinical strains that had undergone limited passaging. Herein, we used targeted next-generation sequencing of the CT134-135 locus, from reference strains and clinical isolates, enabling accurate discovery of single nucleotide polymorphisms and other population genetic variations. Our results indicate that CT134 is stable in all C. trachomatis serovars examined. In contrast, CT135 is highly polymorphic in high-passaged reference ocular and non-LGV genital serovars, with the majority of the mutations resulting in gene disruption. In low-passaged ocular clinical isolates, CT135 was frequently disrupted, whereas in genital clinical isolates CT135 was intact in almost all instances. When a serovar K isolate, with an intact CT134 and CT135, was subjected to serial passage in vitro CT134 remained invariable, while numerous gene interrupting mutations rapidly accumulated in CT135. Collectively, our data indicate that, for genital serovars, CT135 is under strong positive selection in vivo, and negative selection in vitro.

A new pyrosequencing assay for rapid detection and genotyping of Shiga toxin, intimin and O157-specific rfbE genes of Escherichia coli.

Shiga toxin (stx)-producing Escherichia coli (STEC) contamination in food and water is one of the most recognized concerns and a major financial burden in human hygiene control worldwide. Rapid and highly reliable methods of detecting and identifying STEC causing gastroenteric illnesses are crucial to prevent foodborne outbreaks. A number of tests have been developed and commercialized to detect STEC using molecular microbiology techniques. Most of these are designed to identify virulence factors such as Shiga toxin and intimin as well as E. coli O and H antigen serotype specific genes. In order to screen pathogenic STEC without relying on O:H serotyping, we developed a rapid detection and genotyping assay for STEC virulence genes using a PCR-pyrosequencing application. We adapted the PyroMark Q24 Pyrosequencing platform for subtyping 4 major virulence genes, Shiga toxin 1 and 2 (stx1 and stx2), intimin (eae) and O157-antigen gene cluster target rfbE, using Single Nucleotide Polymorphism (SNP) analysis. A total of 224 E. coli strains including isolates from Canadian environment, food and clinical cases were examined. Based on the multiple alignment analysis of 30-80 base nucleotide pyrogram reads, three alleles of the Shiga toxin 1a gene (stx1a) (stx1a-I, stx1a-II, stx1a-III) were identified. Results of the stx1, stx2, eae and rfbE genotyping revealed that each group of O:H serotype shares distinctive characteristics that could be associated with the virulence of each genotype. O157:H7/NM carries stx1a-II (94%), stx2a (82%), λ/γ1-eae (100%) and rfbE type-H7/NM (100%). Whereas isolates of the "Top-6" serotypes (O26, O45, O103, O111, O121, O145) had a high incidence of stx1a-I (90%) and stx2a (100%). stx1a-III (60%) was only observed in non Top-7 (Top-6 plus O157) STEC and Shigella spp. The entire assay, from extracting DNA from colonies on a plate to the generation of sequence information, can be completed in 5h. The method of profiling these 4 STEC pathogenic genotypes as demonstrated in this paper is rapid, easily performed, informative and cost-effective, and thus has a potential to be deployed in the food industry for the routine screening of potentially pathogenic STEC isolates.

A putative iron-regulated TonB-dependent receptor of Mannheimia (Pasteurella) haemolytica A1: possible mechanism for phase variation.

A recombinant plasmid that codes for a novel iron receptor protein (Irp) of Mannheimia (Pasteurella) haemolytica A1 was isolated by the partial complementation of an Escherichia coli fur mutant. The deduced amino acid sequence of Irp exhibited characteristics typical of TonB-dependent receptors. These include: a TonB-box at the N-terminal; a 50 amino acid region homologous to the "plug" domain of the E. coli FhuA and FepA receptors; and a C-terminal TonB-dependent signature which likely functions as an outer membrane anchoring domain. Previously uncharacterized Irp homologues were detected by BLAST analysis of available databases and incomplete microbial genomes. When the irp homologues from Neisseria gonorrhoeae and N. meningitidis were cloned by PCR and expressed in E. coli, novel proteins of the predicted size (84kDa) were detected in cell lysates, demonstrating that these are functional genes. The M. haemolytica A1 irp gene undergoes phase variation at a nucleotide region which contain the sequence AAAAAAATTAAAA (7A-2T-4A) flanked by a short inverted repeat. Site-specific mutagenesis of the 7A-2T-4A sequence as well as replacement of the inverted repeats resulted in a stable construct that expressed the Irp protein without phase variation. The expression of irp in M. haemolytica A1 was regulated by iron concentrations and most likely a Fur homologue, consistent with the proposed function of Irp in iron metabolism. The irp genes may represent contingency loci that play a role in iron acquisition during infection.

Draft Genome Sequence of Streptococcus pyogenes Strain 06BA18369, a Human Pathogen Associated with Skin and Soft Tissue Infections in Northern Canada.

We report the draft sequence of Streptococcus pyogenes 06BA18369 (emm type 41.2, sequence type 579 [ST579]), isolated from a skin and soft tissue infection (SSTI) mixed with Staphylococcus aureus. This genome provides insight into the genetic composition of S. pyogenes strains associated with mixed SSTIs.

Draft Genome Sequence of Methicillin-Susceptible Staphylococcus aureus Strain 06BA18369, a Pathogen Associated with Skin and Soft Tissue Infections in Northern Saskatchewan, Canada.

Here, we announce the draft sequence of a representative methicillin-susceptible Staphylococcus aureus (MSSA) isolate (06BA18369) whose strain type (spa type t311) was commonly isolated from skin and soft tissue coinfections with Streptococcus pyogenes. This strain sequence provides insight into a highly successful community-associated MSSA strain type.

Characterization of a stable, metronidazole-resistant Clostridium difficile clinical isolate.

BACKGROUND: Clostridium difficile are gram-positive, spore forming anaerobic bacteria that are the leading cause of healthcare-associated diarrhea, usually associated with antibiotic usage. Metronidazole is currently the first-line treatment for mild to moderate C. difficile diarrhea however recurrence occurs at rates of 15-35%. There are few reports of C. difficile metronidazole resistance in the literature, and when observed, the phenotype has been transient and lost after storage or exposure of the bacteria to freeze/thaw cycles. Owing to the unstable nature of the resistance phenotype in the laboratory, clinical significance and understanding of the resistance mechanisms is lacking. METHODOLOGY/PRINCIPAL FINDINGS: Genotypic and phenotypic characterization was performed on a metronidazole resistant clinical isolate of C. difficile. Whole-genome sequencing was used to identify potential genetic contributions to the phenotypic variation observed with molecular and bacteriological techniques. Phenotypic observations of the metronidazole resistant strain revealed aberrant growth in broth and elongated cell morphology relative to a metronidazole-susceptible, wild type NAP1 strain. Comparative genomic analysis revealed single nucleotide polymorphism (SNP) level variation within genes affecting core metabolic pathways such as electron transport, iron utilization and energy production. CONCLUSIONS/SIGNIFICANCE: This is the first characterization of stable, metronidazole resistance in a C. difficile isolate. The study provides an in-depth genomic and phenotypic analysis of this strain and provides a foundation for future studies to elucidate mechanisms conferring metronidazole resistance in C. difficile that have not been previously described.

Analysis of the transcriptome of group A Streptococcus in mouse soft tissue infection.

Molecular mechanisms mediating group A Streptococcus (GAS)-host interactions remain poorly understood but are crucial for diagnostic, therapeutic, and vaccine development. An optimized high-density microarray was used to analyze the transcriptome of GAS during experimental mouse soft tissue infection. The transcriptome of a wild-type serotype M1 GAS strain and an isogenic transcriptional regulator knockout mutant (covR) also were compared. Array datasets were verified by quantitative real-time reverse transcriptase-polymerase chain reaction and in situ immunohistochemistry. The results unambiguously demonstrate that coordinated expression of proven and putative GAS virulence factors is directed toward overwhelming innate host defenses leading to severe cellular damage. We also identified adaptive metabolic responses triggered by nutrient signals and hypoxic/acidic conditions in the host, likely facilitating pathogen persistence and proliferation in soft tissues. Key discoveries included that oxidative stress genes, virulence genes, genes related to amino acid and maltodextrin utilization, and several two-component transcriptional regulators were highly expressed in vivo. This study is the first global analysis of the GAS transcriptome during invasive infection. Coupled with parallel analysis of the covR mutant strain, novel insights have been made into the regulation of GAS virulence in vivo, resulting in new avenues for targeted therapeutic and vaccine research.

Group A Streptococcus transcriptome dynamics during growth in human blood reveals bacterial adaptive and survival strategies.

The molecular basis for bacterial responses to host signals during natural infections is poorly understood. The gram-positive bacterial pathogen group A Streptococcus (GAS) causes human mucosal, skin, and life-threatening systemic infections. During the transition from a throat or skin infection to an invasive infection, GAS must adapt to changing environments and host factors. To better understand how GAS adapts, we used transcript profiling and functional analysis to investigate the transcriptome of a wild-type serotype M1 GAS strain in human blood. Global changes in GAS gene expression occur rapidly in response to human blood exposure. Increased transcription was observed for many genes that likely enhance bacterial survival, including those encoding superantigens and host-evasion proteins regulated by a multiple gene activator called Mga. GAS also coordinately expressed genes involved in proteolysis, transport, and catabolism of oligopeptides to obtain amino acids in this protein-rich host environment. Comparison of the transcriptome of the wild-type strain to that of an isogenic deletion mutant (DeltacovR) mutated in the two-component regulatory system designated CovR-CovS reinforced the hypothesis that CovR-CovS has an important role linking key biosynthetic, catabolic, and virulence functions during transcriptome restructuring. Taken together, the data provide crucial insights into strategies used by pathogenic bacteria for thwarting host defenses and surviving in human blood.