Diversity and Evolution in the Genome of Clostridium difficile.

Clostridium difficile infection (CDI) is the leading cause of antimicrobial and health care-associated diarrhea in humans, presenting a significant burden to global health care systems. In the last 2 decades, PCR- and sequence-based techniques, particularly whole-genome sequencing (WGS), have significantly furthered our knowledge of the genetic diversity, evolution, epidemiology, and pathogenicity of this once enigmatic pathogen. C. difficile is taxonomically distinct from many other well-known clostridia, with a diverse population structure comprising hundreds of strain types spread across at least 6 phylogenetic clades. The C. difficile species is defined by a large diverse pangenome with extreme levels of evolutionary plasticity that has been shaped over long time periods by gene flux and recombination, often between divergent lineages. These evolutionary events are in response to environmental and anthropogenic activities and have led to the rapid emergence and worldwide dissemination of virulent clonal lineages. Moreover, genome analysis of large clinically relevant data sets has improved our understanding of CDI outbreaks, transmission, and recurrence. The epidemiology of CDI has changed dramatically over the last 15 years, and CDI may have a foodborne or zoonotic etiology. The WGS era promises to continue to redefine our view of this significant pathogen.

Novel Moraxella catarrhalis prophages display hyperconserved non-structural genes despite their genomic diversity.

BACKGROUND: Moraxella catarrhalis is an important pathogen that often causes otitis media in children, a disease that is not currently vaccine preventable. Asymptomatic colonisation of the human upper respiratory tract is common and lack of clearance by the immune system is likely due to the emergence of seroresistant genetic lineages. No active bacteriophages or prophages have been described in this species. This study was undertaken to identify and categorise prophages in M. catarrhalis, their genetic diversity and the relationship of such diversity with the host-species phylogeny. RESULTS: This study presents a comparative analysis of 32 putative prophages identified in 95 phylogenetically variable, newly sequenced M. catarrhalis genomes. The prophages were genotypically classified into four diverse clades. The genetic synteny of each clade is similar to the group 1 phage family Siphoviridae, however, they form genotypic clusters that are distinct from other members of this family. No core genetic sequences exist across the 32 prophages despite clades 2, 3, and 4 sharing the most sequence identity. The analysis of non-structural prophage genes (coding the integrase, and terminase), and portal gene showed that the respective genes were identical for clades 2, 3, and 4, but unique for clade 1. Empirical analysis calculated that these genes are unexpectedly hyperconserved, under purifying selection, suggesting a tightly regulated functional role. As such, it is improbable that the prophages are decaying remnants but stable components of a fluctuating, flexible and unpredictable system ultimately maintained by functional constraints on non-structural and packaging genes. Additionally, the plate encoding genes were well conserved across all four prophage clades, and the tail fibre genes, commonly responsible for receptor recognition, were clustered into three major groups distributed across the prophage clades. A pan-genome of 283,622 bp was identified, and the prophages were mapped onto the diverse M. catarrhalis multi-locus sequence type (MLST) backbone. CONCLUSION: This study has provided the first evidence of putatively mobile prophages in M. catarrhalis, identifying a diverse and fluctuating system dependent on the hyperconservation of a few key, non-structural genes. Some prophages harbour virulence-related genes, and potentially influence the physiology and virulence of M. catarrhalis. Importantly our data will provide supporting information on the identification of novel prophages in other species by adding greater weight to the identification of non-structural genes.

Evaluation of the Cepheid Xpert C. difficile/Epi and meridian bioscience illumigene C. difficile assays for detecting Clostridium difficile ribotype 033 strains.

Clostridium difficile PCR ribotype 033 (RT033) is found in the gastrointestinal tracts of production animals and, occasionally, humans. The illumigene C. difficile assay (Meridian Bioscience, Inc.) failed to detect any of 52 C. difficile RT033 isolates, while all strains signaled positive for the binary toxin genes but were reported as negative for C. difficile by the Xpert C. difficile/Epi assay (Cepheid).

ChIP-seq and transcriptome analysis of the OmpR regulon of Salmonella enterica serovars Typhi and Typhimurium reveals accessory genes implicated in host colonization.

OmpR is a multifunctional DNA binding regulator with orthologues in many enteric bacteria that exhibits classical regulator activity as well as nucleoid-associated protein-like characteristics. In the enteric pathogen Salmonella enterica, using chromatin immunoprecipitation of OmpR:FLAG and nucleotide sequencing, 43 putative OmpR binding sites were identified in S. enterica serovar Typhi, 22 of which were associated with OmpR-regulated genes. Mutation of a sequence motif (TGTWACAW) that was associated with the putative OmpR binding sites abrogated binding of OmpR:6×His to the tviA upstream region. A core set of 31 orthologous genes were found to exhibit OmpR-dependent expression in both S. Typhi and S. Typhimurium. S. Typhimurium-encoded orthologues of two divergently transcribed OmpR-regulated operons (SL1068-71 and SL1066-67) had a putative OmpR binding site in the inter-operon region in S. Typhi, and were characterized using in vitro and in vivo assays. These operons are widely distributed within S. enterica but absent from the closely related Escherichia coli. SL1066 and SL1067 were required for growth on N-acetylmuramic acid as a sole carbon source. SL1068-71 exhibited sequence similarity to sialic acid uptake systems and contributed to colonization of the ileum and caecum in the streptomycin-pretreated mouse model of colitis.

Sequence-based analysis uncovers an abundance of non-coding RNA in the total transcriptome of Mycobacterium tuberculosis.

RNA sequencing provides a new perspective on the genome of Mycobacterium tuberculosis by revealing an extensive presence of non-coding RNA, including long 5' and 3' untranslated regions, antisense transcripts, and intergenic small RNA (sRNA) molecules. More than a quarter of all sequence reads mapping outside of ribosomal RNA genes represent non-coding RNA, and the density of reads mapping to intergenic regions was more than two-fold higher than that mapping to annotated coding sequences. Selected sRNAs were found at increased abundance in stationary phase cultures and accumulated to remarkably high levels in the lungs of chronically infected mice, indicating a potential contribution to pathogenesis. The ability of tubercle bacilli to adapt to changing environments within the host is critical to their ability to cause disease and to persist during drug treatment; it is likely that novel post-transcriptional regulatory networks will play an important role in these adaptive responses.

Simultaneous assay of every Salmonella Typhi gene using one million transposon mutants.

Very high-throughput sequencing technologies need to be matched by high-throughput functional studies if we are to make full use of the current explosion in genome sequences. We have generated a very large bacterial mutant pool, consisting of an estimated 1.1 million transposon mutants and we have used genomic DNA from this mutant pool, and Illumina nucleotide sequencing to prime from the transposon and sequence into the adjacent target DNA. With this method, which we have called TraDIS (transposon directed insertion-site sequencing), we have been able to map 370,000 unique transposon insertion sites to the Salmonella enterica serovar Typhi chromosome. The unprecedented density and resolution of mapped insertion sites, an average of one every 13 base pairs, has allowed us to assay simultaneously every gene in the genome for essentiality and generate a genome-wide list of candidate essential genes. In addition, the semiquantitative nature of the assay allowed us to identify genes that are advantageous and those that are disadvantageous for growth under standard laboratory conditions. Comparison of the mutant pool following growth in the presence or absence of ox bile enabled every gene to be assayed for its contribution toward bile tolerance, a trait required of any enteric bacterium and for carriage of S. Typhi in the gall bladder. This screen validated our hypothesis that we can simultaneously assay every gene in the genome to identify niche-specific essential genes.

A strand-specific RNA-Seq analysis of the transcriptome of the typhoid bacillus Salmonella typhi.

High-density, strand-specific cDNA sequencing (ssRNA-seq) was used to analyze the transcriptome of Salmonella enterica serovar Typhi (S. Typhi). By mapping sequence data to the entire S. Typhi genome, we analyzed the transcriptome in a strand-specific manner and further defined transcribed regions encoded within prophages, pseudogenes, previously un-annotated, and 3'- or 5'-untranslated regions (UTR). An additional 40 novel candidate non-coding RNAs were identified beyond those previously annotated. Proteomic analysis was combined with transcriptome data to confirm and refine the annotation of a number of hpothetical genes. ssRNA-seq was also combined with microarray and proteome analysis to further define the S. Typhi OmpR regulon and identify novel OmpR regulated transcripts. Thus, ssRNA-seq provides a novel and powerful approach to the characterization of the bacterial transcriptome.

Quantitative RNA-seq analysis of the Campylobacter jejuni transcriptome.

Campylobacter jejuni is the most common bacterial cause of foodborne disease in the developed world. Its general physiology and biochemistry, as well as the mechanisms enabling it to colonize and cause disease in various hosts, are not well understood, and new approaches are required to understand its basic biology. High-throughput sequencing technologies provide unprecedented opportunities for functional genomic research. Recent studies have shown that direct Illumina sequencing of cDNA (RNA-seq) is a useful technique for the quantitative and qualitative examination of transcriptomes. In this study we report RNA-seq analyses of the transcriptomes of C. jejuni (NCTC11168) and its rpoN mutant. This has allowed the identification of hitherto unknown transcriptional units, and further defines the regulon that is dependent on rpoN for expression. The analysis of the NCTC11168 transcriptome was supplemented by additional proteomic analysis using liquid chromatography-MS. The transcriptomic and proteomic datasets represent an important resource for the Campylobacter research community.

A simple method for directional transcriptome sequencing using Illumina technology.

High-throughput sequencing of cDNA has been used to study eukaryotic transcription on a genome-wide scale to single base pair resolution. In order to compensate for the high ribonuclease activity in bacterial cells, we have devised an equivalent technique optimized for studying complete prokaryotic transcriptomes that minimizes the manipulation of the RNA sample. This new approach uses Illumina technology to sequence single-stranded (ss) cDNA, generating information on both the direction and level of transcription throughout the genome. The protocol, and associated data analysis programs, are freely available from http://www.sanger.ac.uk/Projects/Pathogens/Transcriptome/. We have successfully applied this method to the bacterial pathogens Salmonella bongori and Streptococcus pneumoniae and the yeast Schizosaccharomyces pombe. This method enables experimental validation of genetic features predicted in silico and allows the easy identification of novel transcripts throughout the genome. We also show that there is a high correlation between the level of gene expression calculated from ss-cDNA and double-stranded-cDNA sequencing, indicting that ss-cDNA sequencing is both robust and appropriate for use in quantitative studies of transcription. Hence, this simple method should prove a useful tool in aiding genome annotation and gene expression studies in both prokaryotes and eukaryotes.

Differences in the population structure of Neisseria meningitidis in two Australian states: Victoria and Western Australia.

Neisseria meningitidis is the causative agent of invasive meningococcal disease (IMD). A recombinant vaccine called Bexsero® incorporates four subcapsular antigens (fHbp, NHBA, NadA and PorA) which are used to assign a Bexsero® antigen sequence type (BAST) to each meningococcal strain. The vaccine elicits an immune response against combinations of variants of these antigens which have been grouped into specific BAST profiles that have been shown to have different distributions within geographical locations thus potentially affecting the efficacy of the vaccine. In this study, invasive meningococcal disease isolates from the western seaboard of Australia (Western Australia; WA) were compared to those from the south-eastern seaboard (Victoria; VIC) from 2008 to 2012. Whole-genome sequencing (WGS) of 131 meningococci from VIC and 70 meningococci from WA were analysed for MLST, FetA and BAST profiling. Serogroup B predominated in both jurisdictions and a total of 10 MLST clonal complexes (cc) were shared by both states. Isolates belonging to cc22, cc103 and cc1157 were unique to VIC whilst isolates from cc60 and cc212 were unique to WA. Clonal complex 41/44 represented one-third of the meningococcal population in each state but the predominant ST was locally different: ST-6058 in VIC and ST-146 in WA. Of the 108 BAST profiles identified in this collection, only 9 BASTs were simultaneously observed in both states. A significantly larger proportion of isolates in VIC harboured alleles for the NHBA-2 peptide and fHbp-1, antigenic variants predicted to be covered by the Bexsero® vaccine. The estimate for vaccine coverage in WA (47.1% [95% CI: 41.1-53.1%]) was significantly lower than that in VIC (66.4% [95% CI: 62.3-70.5%]). In conclusion, the antigenic structure of meningococci causing invasive disease in two geographically distinct states of Australia differed significantly during the study period which may affect vaccine effectiveness and highlights the need for representative surveillance when predicting potential impact of meningococcal B vaccines.

Acquisition of the capsule locus by horizontal gene transfer in Neisseria meningitidis is often accompanied by the loss of UDP-GalNAc synthesis.

Pathogenic meningococci have acquired a 24 kb capsule synthesis island (cps) by horizontal gene transfer which consists of a synthetic locus and associated capsule transport genes flanked by repetitive Regions D and D'. Regions D and D' contain an intact gene encoding a UDP-galactose epimerase (galE1) and a truncated remnant (galE2), respectively. In this study, GalE protein alleles were shown to be either mono-functional, synthesising UDP-galactose (UDP-Gal), or bi-functional, synthesising UDP-Gal and UDP-galactosamine (UDP-GalNAc). Meningococci possessing a capsule null locus (cnl) typically possessed a single bi-functional galE. Separation of functionality between galE1 and galE2 alleles in meningococcal isolates was retained for all serogroups except serogroup E which has a synthetic requirement for UDP-GalNAc. The truncated galE2 remnant in Region D' was also phylogenetically related to the bi-functional galE of the cnl locus suggesting common ancestry. A model is proposed in which the illegitimate recombination of the cps island into the galE allele of the cnl locus results in the formation of Region D' containing the truncated galE2 locus and the capture of the cps island en bloc. The retention of the duplicated Regions D and D' enables inversion of the synthetic locus within the cps island during bacterial growth.

Quantum changes in Helicobacter pylori gene expression accompany host-adaptation.

Helicobacter pylori is a highly successful gastric pathogen. High genomic plasticity allows its adaptation to changing host environments. Complete genomes of H. pylori clinical isolate UM032 and its mice-adapted serial derivatives 298 and 299, generated using both PacBio RS and Illumina MiSeq sequencing technologies, were compared to identify novel elements responsible for host-adaptation. The acquisition of a jhp0562-like allele, which encodes for a galactosyltransferase, was identified in the mice-adapted strains. Our analysis implies a new ß-1,4-galactosyltransferase role for this enzyme, essential for Ley antigen expression. Intragenomic recombination between babA and babB genes was also observed. Further, we expanded on the list of candidate genes whose expression patterns have been mediated by upstream homopolymer-length alterations to facilitate host adaption. Importantly, greater than four-fold reduction of mRNA levels was demonstrated in five genes. Among the down-regulated genes, three encode for outer membrane proteins, including BabA, BabB and HopD. As expected, a substantial reduction in BabA protein abundance was detected in mice-adapted strains 298 and 299 via Western analysis. Our results suggest that the expression of Ley antigen and reduced outer membrane protein expressions may facilitate H. pylori colonisation of mouse gastric epithelium.

Temporal Changes in BEXSERO® Antigen Sequence Type Associated with Genetic Lineages of Neisseria meningitidis over a 15-Year Period in Western Australia.

Neisseria meningitidis is the causative agent of invasive meningococcal disease (IMD). The BEXSERO® vaccine which is used to prevent serogroup B disease is composed of four sub-capsular protein antigens supplemented with an outer membrane vesicle. Since the sub-capsular protein antigens are variably expressed and antigenically variable amongst meningococcal isolates, vaccine coverage can be estimated by the meningococcal antigen typing system (MATS) which measures the propensity of the strain to be killed by vaccinated sera. Whole genome sequencing (WGS) which identifies the alleles of the antigens that may be recognised by the antibody response could represent, in future, an alternative estimate of coverage. In this study, WGS of 278 meningococcal isolates responsible for 62% of IMD in Western Australia from 2000-2014 were analysed for association of genetic lineage (sequence type [ST], clonal complex [cc]) with BEXSERO® antigen sequence type (BAST) and MATS to predict the annual vaccine coverage. A hyper-endemic period of IMD between 2000-05 was caused by cc41/44 with the major sequence type of ST-146 which was not predicted by MATS or BAST to be covered by the vaccine. An increase in serogroup diversity was observed between 2010-14 with the emergence of cc11 serogroup W in the adolescent population and cc23 serogroup Y in the elderly. BASTs were statistically associated with clonal complex although individual antigens underwent antigenic drift from the major type. BAST and MATS predicted an annual range of 44-91% vaccine coverage. Periods of low vaccine coverage in years post-2005 were not a result of the resurgence of cc41/44:ST-146 but were characterised by increased diversity of clonal complexes expressing BASTs which were not predicted by MATS to be covered by the vaccine. The driving force behind the diversity of the clonal complex and BAST during these periods of low vaccine coverage is unknown, but could be due to immune selection and inter-strain competition with carriage of non-disease causing meningococci.

Fucoidans Disrupt Adherence of Helicobacter pylori to AGS Cells In Vitro.

Fucoidans are complex sulphated polysaccharides derived from abundant and edible marine algae. Helicobacter pylori is a stomach pathogen that persists in the hostile milieu of the human stomach unless treated with antibiotics. This study aims to provide preliminary data to determine, in vitro, if fucoidans can inhibit the growth of H. pylori and its ability to adhere to gastric epithelial cells (AGS). We analysed the activity of three different fucoidan preparations (Fucus A, Fucus B, and Undaria extracts). Bacterial growth was not arrested or inhibited by the fucoidan preparations supplemented into culture media. All fucoidans, when supplemented into tissue culture media at 1000 µg mL(-1), were toxic to AGS cells and reduced the viable cell count significantly. Fucoidan preparations at 100 µg mL(-1) were shown to significantly reduce the number of adherent H. pylori. These in vitro findings provide the basis for further studies on the clinical use of sulphated polysaccharides as complementary therapeutic agents.

Choosing a benchtop sequencing machine to characterise Helicobacter pylori genomes.

The fully annotated genome sequence of the European strain, 26695 was first published in 1997 and, in 1999, it was directly compared to the USA isolate J99, promoting two standard laboratory isolates for Helicobacter pylori (H. pylori) research. With the genomic scaffolds available from these important genomes and the advent of benchtop high-throughput sequencing technology, a bacterial genome can now be sequenced within a few days. We sequenced and analysed strains J99 and 26695 using the benchtop-sequencing machines Ion Torrent PGM and the Illumina MiSeq Nextera and Nextera XT methodologies. Using publically available algorithms, we analysed the raw data and interrogated both genomes by mapping the data and by de novo assembly. We compared the accuracy of the coding sequence assemblies to the originally published sequences. With the Ion Torrent PGM, we found an inherently high-error rate in the raw sequence data. Using the Illumina MiSeq, we found significantly more non-covered nucleotides when using the less expensive Illumina Nextera XT compared with the Illumina Nextera library creation method. We found the most accurate de novo assemblies using the Nextera technology, however, extracting an accurate multi-locus sequence type was inconsistent compared to the Ion Torrent PGM. We found the cagPAI failed to assemble onto a single contig in all technologies but was more accurate using the Nextera. Our results indicate the Illumina MiSeq Nextera method is the most accurate for de novo whole genome sequencing of H. pylori.