Comparative genomic analyses of seventeen Streptococcus pneumoniae strains: insights into the pneumococcal supragenome

On the origin of Mycobacterium ulcerans, the causative agent of Buruli ulcer

Background

Mycobacterium ulcerans is an unusual bacterial pathogen with elusive origins. While closely related to the aquatic dwelling M. marinum, M. ulcerans has evolved the ability to produce the immunosuppressive polyketide toxin mycolactone and cause the neglected tropical disease Buruli ulcer. Other mycolactone-producing mycobacteria (MPM) have been identified in fish and frogs and given distinct species designations (M. pseudoshottsii, M. shinshuense, M. liflandii and M. marinum), however the evolution of M. ulcerans and its relationship to other MPM has not been defined. Here we report the comparative analysis of whole genome sequences from 30 MPM and five M. marinum.

Results

A high-resolution phylogeny based on genome-wide single nucleotide polymorphisms (SNPs) showed that M. ulcerans and all other MPM represent a single clonal group that evolved from a common M. marinum progenitor. The emergence of the MPM was driven by the acquisition of the pMUM plasmid encoding genes for the biosynthesis of mycolactones. This change was accompanied by the loss of at least 185 genes, with a significant overrepresentation of genes associated with cell wall functions. Cell wall associated genes also showed evidence of substantial adaptive selection, suggesting cell wall remodeling has been critical for the survival of MPM. Fine-grain analysis of the MPM complex revealed at least three distinct lineages, one of which comprised a highly clonal group, responsible for Buruli ulcer in Africa and Australia. This indicates relatively recent transfer of M. ulcerans between these continents, which represent the vast majority of the global Buruli ulcer burden. Our data provide SNPs and gene sequences that can differentiate M. ulcerans lineages, suitable for use in the diagnosis and surveillance of Buruli ulcer.

Conclusions

M. ulcerans and all mycolactone-producing mycobacteria are specialized variants of a common Mycobacterium marinum progenitor that have adapted to live in restricted environments. Examination of genes lost or retained and now under selective pressure suggests these environments might be aerobic, and extracellular, where slow growth, production of an immune suppressor, cell wall remodeling, loss or modification of cell wall antigens, and biofilm-forming ability provide a survival advantage. These insights will guide our efforts to find the elusive reservoir(s) of M. ulcerans and to understand transmission of Buruli ulcer.

A neutral theory of genome evolution and the frequency distribution of genes

Background

The gene composition of bacteria of the same species can differ significantly between isolates. Variability in gene composition can be summarized in terms of gene frequency distributions, in which individual genes are ranked according to the frequency of genomes in which they appear. Empirical gene frequency distributions possess a U-shape, such that there are many rare genes, some genes of intermediate occurrence, and many common genes. It would seem that U-shaped gene frequency distributions can be used to infer the essentiality and/or importance of a gene to a species. Here, we ask: can U-shaped gene frequency distributions, instead, arise generically via neutral processes of genome evolution?

Results

We introduce a neutral model of genome evolution which combines birth-death processes at the organismal level with gene uptake and loss at the genomic level. This model predicts that gene frequency distributions possess a characteristic U-shape even in the absence of selective forces driving genome and population structure. We compare the model predictions to empirical gene frequency distributions from 6 multiply sequenced species of bacterial pathogens. We fit the model with constant population size to data, matching U-shape distributions albeit without matching all quantitative features of the distribution. We find stronger model fits in the case where we consider exponentially growing populations. We also show that two alternative models which contain a "rigid" and "flexible" core component of genomes provide strong fits to gene frequency distributions.

Conclusions

The analysis of neutral models of genome evolution suggests that U-shaped gene frequency distributions provide less information than previously suggested regarding gene essentiality. We discuss the need for additional theory and genomic level information to disentangle the roles of evolutionary mechanisms operating within and amongst individuals in driving the dynamics of gene distributions.

Comparative genomic analyses of 17 clinical isolates of Gardnerella vaginalis provide evidence of multiple genetically isolated clades consistent with subspeciation into genovars

Gardnerella vaginalis is associated with a spectrum of clinical conditions, suggesting high degrees of genetic heterogeneity among stains. Seventeen G. vaginalis isolates were subjected to a battery of comparative genomic analyses to determine their level of relatedness. For each measure, the degree of difference among the G. vaginalis strains was the highest observed among 23 pathogenic bacterial species for which at least eight genomes are available. Genome sizes ranged from 1.491 to 1.716 Mb; GC contents ranged from 41.18% to 43.40%; and the core genome, consisting of only 746 genes, makes up only 51.6% of each strain's genome on average and accounts for only 27% of the species supragenome. Neighbor-grouping analyses, using both distributed gene possession data and core gene allelic data, each identified two major sets of strains, each of which is composed of two groups. Each of the four groups has its own characteristic genome size, GC ratio, and greatly expanded core gene content, making the genomic diversity of each group within the range for other bacterial species. To test whether these 4 groups corresponded to genetically isolated clades, we inferred the phylogeny of each distributed gene that was present in at least two strains and absent in at least two strains; this analysis identified frequent homologous recombination within groups but not between groups or sets. G. vaginalis appears to include four nonrecombining groups/clades of organisms with distinct gene pools and genomic properties, which may confer distinct ecological properties. Consequently, it may be appropriate to treat these four groups as separate species.

Unique core genomes of the bacterial family vibrionaceae: insights into niche adaptation and speciation

Background

The criteria for defining bacterial species and even the concept of bacterial species itself are under debate, and the discussion is apparently intensifying as more genome sequence data is becoming available. However, it is still unclear how the new advances in genomics should be used most efficiently to address this question. In this study we identify genes that are common to any group of genomes in our dataset, to determine whether genes specific to a particular taxon exist and to investigate their potential role in adaptation of bacteria to their specific niche. These genes were named unique core genes. Additionally, we investigate the existence and importance of unique core genes that are found in isolates of phylogenetically non-coherent groups. These groups of isolates, that share a genetic feature without sharing a closest common ancestor, are termed genophyletic groups.

Results

The bacterial family Vibrionaceae was used as the model, and we compiled and compared genome sequences of 64 different isolates. Using the software orthoMCL we determined clusters of homologous genes among the investigated genome sequences. We used multilocus sequence analysis to build a host phylogeny and mapped the numbers of unique core genes of all distinct groups of isolates onto the tree. The results show that unique core genes are more likely to be found in monophyletic groups of isolates. Genophyletic groups of isolates, in contrast, are less common especially for large groups of isolate. The subsequent annotation of unique core genes that are present in genophyletic groups indicate a high degree of horizontally transferred genes. Finally, the annotation of the unique core genes of Vibrio cholerae revealed genes involved in aerotaxis and biosynthesis of the iron-chelator vibriobactin.

Conclusion

The presented work indicates that genes specific for any taxon inside the bacterial family Vibrionaceae exist. These unique core genes encode conserved metabolic functions that can shed light on the adaptation of a species to its ecological niche. Additionally, our study suggests that unique core genes can be used to aid classification of bacteria and contribute to a bacterial species definition on a genomic level. Furthermore, these genes may be of importance in clinical diagnostics and drug development.

Pangenomic study of Corynebacterium diphtheriae that provides insights into the genomic diversity of pathogenic isolates from cases of classical diphtheria, endocarditis, and pneumonia

Corynebacterium diphtheriae is one of the most prominent human pathogens and the causative agent of the communicable disease diphtheria. The genomes of 12 strains isolated from patients with classical diphtheria, endocarditis, and pneumonia were completely sequenced and annotated. Including the genome of C. diphtheriae NCTC 13129, we herewith present a comprehensive comparative analysis of 13 strains and the first characterization of the pangenome of the species C. diphtheriae. Comparative genomics showed extensive synteny and revealed a core genome consisting of 1,632 conserved genes. The pangenome currently comprises 4,786 protein-coding regions and increases at an average of 65 unique genes per newly sequenced strain. Analysis of prophages carrying the diphtheria toxin gene tox revealed that the toxoid vaccine producer C. diphtheriae Park-Williams no. 8 has been lysogenized by two copies of the ω(tox)(+) phage, whereas C. diphtheriae 31A harbors a hitherto-unknown tox(+) corynephage. DNA binding sites of the tox-controlling regulator DtxR were detected by genome-wide motif searches. Comparative content analysis showed that the DtxR regulons exhibit marked differences due to gene gain, gene loss, partial gene deletion, and DtxR binding site depletion. Most predicted pathogenicity islands of C. diphtheriae revealed characteristics of horizontal gene transfer. The majority of these islands encode subunits of adhesive pili, which can play important roles in adhesion of C. diphtheriae to different host tissues. All sequenced isolates contain at least two pilus gene clusters. It appears that variation in the distributed genome is a common strategy of C. diphtheriae to establish differences in host-pathogen interactions.

A functional genomics approach to establish the complement of carbohydrate transporters in Streptococcus pneumoniae

The aerotolerant anaerobe Streptococcus pneumoniae is part of the normal nasopharyngeal microbiota of humans and one of the most important invasive pathogens. A genomic survey allowed establishing the occurrence of twenty-one phosphotransferase systems, seven carbohydrate uptake ABC transporters, one sodium:solute symporter and a permease, underlining an exceptionally high capacity for uptake of carbohydrate substrates. Despite high genomic variability, combined phenotypic and genomic analysis of twenty sequenced strains did assign the substrate specificity only to two uptake systems. Systematic analysis of mutants for most carbohydrate transporters enabled us to assign a phenotype and substrate specificity to twenty-three transport systems. For five putative transporters for galactose, pentoses, ribonucleosides and sulphated glycans activity was inferred, but not experimentally confirmed and only one transport system remains with an unknown substrate and lack of any functional annotation. Using a metabolic approach, 80% of the thirty-two fermentable carbon substrates were assigned to the corresponding transporter. The complexity and robustness of sugar uptake is underlined by the finding that many transporters have multiple substrates, and many sugars are transported by more than one system. The present work permits to draw a functional map of the complete arsenal of carbohydrate utilisation proteins of pneumococci, allows re-annotation of genomic data and might serve as a reference for related species. These data provide tools for specific investigation of the roles of the different carbon substrates on pneumococcal physiology in the host during carriage and invasive infection.

The infinitely many genes model for the distributed genome of bacteria

The distributed genome hypothesis states that the gene pool of a bacterial taxon is much more complex than that found in a single individual genome. However, the possible fitness advantage, why such genomic diversity is maintained, whether this variation is largely adaptive or neutral, and why these distinct individuals can coexist, remains poorly understood. Here, we present the infinitely many genes (IMG) model, which is a quantitative, evolutionary model for the distributed genome. It is based on a genealogy of individual genomes and the possibility of gene gain (from an unbounded reservoir of novel genes, e.g., by horizontal gene transfer from distant taxa) and gene loss, for example, by pseudogenization and deletion of genes, during reproduction. By implementing these mechanisms, the IMG model differs from existing concepts for the distributed genome, which cannot differentiate between neutral evolution and adaptation as drivers of the observed genomic diversity. Using the IMG model, we tested whether the distributed genome of 22 full genomes of picocyanobacteria (Prochlorococcus and Synechococcus) shows signs of adaptation or neutrality. We calculated the effective population size of Prochlorococcus at 1.01 × 10¹¹ and predicted 18 distinct clades for this population, only six of which have been isolated and cultured thus far. We predicted that the Prochlorococcus pangenome contains 57,792 genes and found that the evolution of the distributed genome of Prochlorococcus was possibly neutral, whereas that of Synechococcus and the combined sample shows a clear deviation from neutrality.

Differential occurrence of Streptococcus pneumoniae serotype 11E between asymptomatic carriage and invasive pneumococcal disease isolates reflects a unique model of pathogen microevolution

BACKGROUND

Streptococcus pneumoniae is a commensal colonizer of the human nasopharynx (NP) that causes disease after evasion of host defenses and dissemination. Pneumococcal strains expressing the newly identified serotype 11E arise from antigenically similar 11A progenitors by genetic inactivation of the O-acetyltransferase gene wcjE. Each 11E strain contains a distinct mutation to wcjE, suggesting that 11E strains are not transmitted among hosts despite their recovery from multiple patients with pneumococcal disease. We investigated whether the presumed lack of transmission of serotype 11E is consistent with its inability to survive in the NP.

METHODS

More than 400 pneumococcal carriage, middle ear, conjunctiva, and blood isolates, serotyped as 11A by Quellung reaction, were reexamined for reactivity to 11A- and 11E-specific antibodies. We confirmed serotyping of isolates with sequencing of wcjE alleles.

RESULTS

Serotype 11E strains were statistically more likely to occur among blood (4 of 15), conjunctiva (1 of 14), or middle ear (2 of 21) isolates than among carriage isolates (2 of 355). All 11E isolates contained unique mutations that putatively decrease wcjE expression.

CONCLUSIONS

The lack of a circulating 11E clone and the increased occurrence of 11E strains among disease isolates supports the idea that serotype 11E independently arises during infection after initial colonization with a serotype 11A progenitor. Factors encountered in the NP likely contribute to relative rarity of 11E among carriage isolates, whereas selective pressures in deeper tissues possibly promote 11E emergence. These findings illustrate a novel model of microevolution that transpires during the span of a single encounter with serotype 11A, highlighting the adaptability of bacterial pathogens within hosts.

Pneumococcal gene complex involved in resistance to extracellular oxidative stress

Streptococcus pneumoniae is a gram-positive bacterium which is a member of the normal human nasopharyngeal flora but can also cause serious disease such as pneumonia, bacteremia, and meningitis. Throughout its life cycle, S. pneumoniae is exposed to significant oxidative stress derived from endogenously produced hydrogen peroxide (H(2)O(2)) and from the host through the oxidative burst. How S. pneumoniae, an aerotolerant anaerobic bacterium that lacks catalase, protects itself against hydrogen peroxide stress is still unclear. Bioinformatic analysis of its genome identified a hypothetical open reading frame belonging to the thiol-specific antioxidant (TlpA/TSA) family, located in an operon consisting of three open reading frames. For all four strains tested, deletion of the gene resulted in an approximately 10-fold reduction in survival when strains were exposed to external peroxide stress. However, no role for this gene in survival of internal superoxide stress was observed. Mutagenesis and complementation analysis demonstrated that all three genes are necessary and sufficient for protection against oxidative stress. Interestingly, in a competitive index mouse pneumonia model, deletion of the operon had no impact shortly after infection but was detrimental during the later stages of disease. Thus, we have identified a gene complex involved in the protection of S. pneumoniae against external oxidative stress, which plays an important role during invasive disease.

In silico analysis of sequenced strains of Clostridium difficile reveals a related set of conjugative transposons carrying a variety of accessory genes

The human gut pathogen Clostridium difficile contains many conjugative transposons that have an array of accessory genes. In the current study, recently sequenced genomes were analyzed to identify new putative conjugative transposons. Eleven new elements in 5 C. difficile strains were identified and all had a similar structure to the previously described elements CTn1, CTn5 and CTn7 in C. difficile strain 630. Each element identified did however contain a new set of accessory genes compared with those previously reported; including those predicted to encode ABC transporters, a toxin/antitoxin system and multiple antibiotic resistance genes.

Differences in genotype and virulence among four multidrug-resistant Streptococcus pneumoniae isolates belonging to the PMEN1 clone

We report on the comparative genomics and characterization of the virulence phenotypes of four S. pneumoniae strains that belong to the multidrug resistant clone PMEN1 (Spain^23F ST81). Strains SV35-T23 and SV36-T3 were recovered in 1996 from the nasopharynx of patients at an AIDS hospice in New York. Strain SV36-T3 expressed capsule type 3 which is unusual for this clone and represents the product of an in vivo capsular switch event. A third PMEN1 isolate – PN4595-T23 – was recovered in 1996 from the nasopharynx of a child attending day care in Portugal, and a fourth strain – ATCC700669 – was originally isolated from a patient with pneumococcal disease in Spain in 1984. We compared the genomes among four PMEN1 strains and 47 previously sequenced pneumococcal isolates for gene possession differences and allelic variations within core genes. In contrast to the 47 strains – representing a variety of clonal types – the four PMEN1 strains grouped closely together, demonstrating high genomic conservation within this lineage relative to the rest of the species. In the four PMEN1 strains allelic and gene possession differences were clustered into 18 genomic regions including the capsule, the blp bacteriocins, erythromycin resistance, the MM1-2008 prophage and multiple cell wall anchored proteins. In spite of their genomic similarity, the high resolution chinchilla model was able to detect variations in virulence properties of the PMEN1 strains highlighting how small genic or allelic variation can lead to significant changes in pathogenicity and making this set of strains ideal for the identification of novel virulence determinants.

Contribution of serotype and genetic background to virulence of serotype 3 and serogroup 11 pneumococcal isolates

The capsular serotype has long been associated with the virulence of Streptococcus pneumoniae. Here we present an in-depth study of phenotypic and genetic differences between serotype 3 and serogroup 11 S. pneumoniae clinical isolates from both the general and indigenous populations of Australia. Both serotypes/groups included clonally unrelated strains with differences in well-known polymorphic virulence genes, such as nanA and pspA, as demonstrated by multilocus sequence typing and Western blot analysis. Nonetheless, the serotype 3 strains were consistently and significantly more virulent in mice than the serogroup 11 strains. Despite extensive genomic analysis, noncapsular genes common to one serotype/group but not the other were not identified. Nevertheless, following the conversion of a serotype 11A isolate to serotype 3 and subsequent analysis in an intranasal infection model, it was evident that both capsular and noncapsular factors determine the virulence phenotype in mice. However, it appears that these noncapsular factors vary from strain to strain.

PGAP: pan-genomes analysis pipeline

Summary: With the rapid development of DNA sequencing technology, increasing bacteria genome data enable the biologists to dig the evolutionary and genetic information of prokaryotic species from pan-genome sight. Therefore, the high-efficiency pipelines for pan-genome analysis are mostly needed. We have developed a new pan-genome analysis pipeline (PGAP), which can perform five analytic functions with only one command, including cluster analysis of functional genes, pan-genome profile analysis, genetic variation analysis of functional genes, species evolution analysis and function enrichment analysis of gene clusters. PGAP's performance has been evaluated on 11 Streptococcus pyogenes strains.

Availability:PGAP is developed with Perl script on the Linux Platform and the package is freely available from http://pgap.sf.net.

Contact:junyu@big.ac.cn; xiaojingfa@big.ac.cn

Supplementary information:Supplementary data are available at Bioinformatics online.

Whole-genome comparison of two Campylobacter jejuni isolates of the same sequence type reveals multiple loci of different ancestral lineage

Campylobacter jejuni ST-474 is the most important human enteric pathogen in New Zealand, and yet this genotype is rarely found elsewhere in the world. Insight into the evolution of this organism was gained by a whole genome comparison of two ST-474, flaA SVR-14 isolates and other available C. jejuni isolates and genomes. The two isolates were collected from different sources, human (H22082) and retail poultry (P110b), at the same time and from the same geographical location. Solexa sequencing of each isolate resulted in 1.659 Mb (H22082) and 1.656 Mb (P110b) of assembled sequences within 28 (H22082) and 29 (P110b) contigs. We analysed 1502 genes for which we had sequences within both ST-474 isolates and within at least one of 11 C. jejuni reference genomes. Although 94.5% of genes were identical between the two ST-474 isolates, we identified 83 genes that differed by at least one nucleotide, including 55 genes with non-synonymous substitutions. These covered 101 kb and contained 672 point differences. We inferred that 22 (3.3%) of these differences were due to mutation and 650 (96.7%) were imported via recombination. Our analysis estimated 38 recombinant breakpoints within these 83 genes, which correspond to recombination events affecting at least 19 loci regions and gives a tract length estimate of 2 kb. This includes a 12 kb region displaying non-homologous recombination in one of the ST-474 genomes, with the insertion of two genes, including ykgC, a putative oxidoreductase, and a conserved hypothetical protein of unknown function. Furthermore, our analysis indicates that the source of this recombined DNA is more likely to have come from C. jejuni strains that are more closely related to ST-474. This suggests that the rates of recombination and mutation are similar in order of magnitude, but that recombination has been much more important for generating divergence between the two ST-474 isolates.

Effect of new alleles of the histidine kinase gene ciaH on the activity of the response regulator CiaR in Streptococcus pneumoniae R6

The two-component regulatory system CiaRH of Streptococcus pneumoniae affects β-lactam susceptibility, autolysis, bacteriocin production, competence development, host colonization and virulence. The system was discovered in a screen for S. pneumoniae R6 mutants resistant to the β-lactam antibiotic cefotaxime. A mutation in the histidine kinase gene ciaH led to this phenotype by enhancing CiaR-mediated gene expression. Additional mutations in ciaH have been described in other spontaneous β-lactam-resistant mutants of S. pneumoniae R6, but their influence on CiaR-mediated gene regulation has not been determined. Likewise, altered ciaH alleles are present in clinical S. pneumoniae isolates, none of which had been characterized. These novel ciaH variants were introduced into S. pneumoniae R6 to measure their ability to activate CiaR-dependent regulation. The ciaH alleles from spontaneous mutants obtained in the laboratory increased the activity of CiaR-dependent promoters between four- and 26-fold, while variants from clinical strains were less effective, with a threefold activation at most. Accordingly, phenotypes associated with a hyperactive CiaRH system, β-lactam resistance, and prevention of competence development, were far more pronounced in the laboratory mutants. Amino acid changes affecting CiaH function were positioned throughout the protein. Five of the most activating changes are located close to the conserved histidine and one in the extracytoplasmic sensor domain. The characterization of new alleles of ciaH expands the spectrum of CiaH variants, which may help to elucidate signal transduction of this important regulatory system. Our study also demonstrates that ciaH alleles overstimulating CiaR regulon expression are present in clinical isolates of S. pneumoniae.

Ecological and temporal constraints in the evolution of bacterial genomes

Studies on the experimental evolution of microorganisms, on their in vivo evolution (mainly in the case of bacteria producing chronic infections), as well as the availability of multiple full genomic sequences, are placing bacteria in the playground of evolutionary studies. In the present article we review the differential contribution to the evolution of bacterial genomes that processes such as gene modification, gene acquisition and gene loss may have when bacteria colonize different habitats that present characteristic ecological features. In particular, we review how the different processes contribute to evolution in microbial communities, in free-living bacteria or in bacteria living in isolation. In addition, we discuss the temporal constraints in the evolution of bacterial genomes, considering bacterial evolution from the perspective of processes of short-sighted evolution and punctual acquisition of evolutionary novelties followed by long stasis periods.

Population genetics and evolution of the pan-genome of Streptococcus pneumoniae

The genetic variability in bacterial species is much larger than in other kingdoms of life. The gene content between pairs of isolates can diverge by as much as 30% in species like Escherichia coli or Streptococcus pneumoniae. This unexpected finding led to the introduction of the concept of the pan-genome, the set of genes that can be found in a given bacterial species. The genome of any isolate is thus composed by a "core genome" shared by all strains and characteristic of the species, and a "dispensable genome" that accounts for many of the phenotypic differences between strains. The pan-genome is usually much larger than the genome of any single isolate and, given the ability of many bacteria to exchange genetic material with the environment, constitutes a reservoir that could enhance their ability to survive in a mutating environment. To understand the evolution of the pan-genome of an important pathogen and its interactions with the commensal microbial flora, we have analyzed the genomes of 44 strains of Streptococcus pneumoniae, one of the most important causes of microbial diseases in humans. Despite evidence of extensive homologous recombination, the S. pneumoniae phylogenetic tree reconstructed from polymorphisms in the core genome identified major groups of genetically related strains. With the exception of serotype 1, the tree correlated poorly with capsular serotype, geographical site of isolation and disease outcome. The distribution of dispensable genes was consistent with phylogeny, although horizontal gene transfer events attenuated this correlation in the case of ancient lineages. Homologous recombination, involving short stretches of DNA, was the dominant evolutionary process of the core genome of S. pneumoniae. Genetic exchange with related species sharing the same ecological niche was the main mechanism of evolution of S. pneumonia; and S. mitis was the main reservoir of genetic diversity of S. pneumoniae. The pan-genome of S. pneumoniae increased logarithmically with the number of strains and linearly with the variability of the sample, suggesting that acquired genes accumulate proportionately to the age of clones.

Effects of Streptococcus pneumoniae strain background on complement resistance

Background

Immunity to infections caused by Streptococcus pneumoniae is dependent on complement. There are wide variations in sensitivity to complement between S. pneumoniae strains that could affect their ability to cause invasive infections. Although capsular serotype is one important factor causing differences in complement resistance between strains, there is also considerable other genetic variation between S. pneumoniae strains that may affect complement-mediated immunity. We have therefore investigated whether genetically distinct S. pneumoniae strains with the same capsular serotype vary in their sensitivity to complement mediated immunity.

Methodology and Principal Findings

C3b/iC3b deposition and neutrophil association were measured using flow cytometry assays for S. pneumoniae strains with different genetic backgrounds for each of eight capsular serotypes. For some capsular serotypes there was marked variation in C3b/iC3b deposition between different strains that was independent of capsule thickness and correlated closely to susceptibility to neutrophil association. C3b/iC3b deposition results also correlated weakly with the degree of IgG binding to each strain. However, the binding of C1q (the first component of the classical pathway) correlated more closely with C3b/iC3b deposition, and large differences remained in complement sensitivity between strains with the same capsular serotype in sera in which IgG had been cleaved with IdeS.

Conclusions

These data demonstrate that bacterial factors independent of the capsule and recognition by IgG have strong effects on the susceptibility of S. pneumoniae to complement, and could therefore potentially account for some of the differences in virulence between strains.

A new variant of the capsule 3 cluster occurs in Streptococcus pneumoniae from deceased wild chimpanzees

The presence of new Streptococcus pneumoniae clones in dead wild chimpanzees from the Taï National Park, Côte d'Ivoire, with previous respiratory problems has been demonstrated recently by DNA sequence analysis from samples obtained from the deceased apes. In order to broadenour understanding on the relatedness of these pneumococcal clones to those from humans, the gene locus responsible for biosynthesis of the capsule polysaccharide (CPS) has now been characterized. DNA sequence analysis of PCR fragments identified a cluster named cps3(Taï) containing the four genes typical for serotype 3 CPS, but lacking a 5'-region of ≥2 kb which is degenerated in other cps3 loci and not required for type 3 biosynthesis. CPS3 is composed of a simple disaccharide repeat unit comprising glucose and glucuronic acid (GlcUA). The two genes ugd responsible for GlcUA synthesis and wchE encoding the type 3 synthase are essential for CPS3 biosynthesis, whereas both, galU and the 3'-truncated gene pgm are not required due to the presence of homologues elsewhere in the genome. The DNA sequence of cps3(Taï) diverged considerably from those of other cps3 loci. Also, the gene pgm(Taï) represents a full length version with a nonsense mutation at codon 179. The two genes ugd(Taï) and wchE(Taï) including the promoter region were transformed into a nonencapsulated laboratory strain S. pneumoniae R6. Transformants which expressed type 3 capsule polysaccharide were readily obtained, documenting that the gene products are functional. In summary, the data indicate that cps3(Taï) evolved independent from other cps3 loci, suggesting the presence of specialized serotype 3 S. pneumoniae clones endemic to the Taï National Park area.

The adaptive potential during nasopharyngeal colonisation of Streptococcus pneumoniae

Adaptation to host defences and antimicrobials is critical for Streptococcus pneumoniae (the pneumococcus) during colonisation of the nasopharynx--its only ecological habitat. The pneumococcus is highly transformable with the genome between different strains varying widely in both gene content and gene sequence. Thus, mixed strains colonising together will expand the genetic reservoir--"supragenome" for this highly transformable microorganism, increasing its adaptive potential. The extent of the phenotypic and genotypic diversity of strains co-colonising in the nasopharynx was determined. In contrast to most carriage studies, which characterise single colonies, a systematic analysis of up to 20 colonies per colonisation was undertaken in Tanzanian children for 12 months. The serotype was determined by conventional serology and confirmed by DNA-based methods. The antibiotype for penicillin and co-trimoxazole was determined from the minimum inhibitory concentration determined by E-test. As representative of the genotype of strains the sequence types (STs) was determined by multilocus sequence typing (MLST). Of 61 colonisation events studied, seven (11.5%) had strains expressing multiple serotypes, with a maximum of five serotypes detected. Four colonisation events also had co-colonisation of penicillin and/or co-trimoxazole susceptible and non-susceptible pneumococci. Sequence typing revealed that 58% were unique to our cohort. Simultaneous colonisation of up to six STs with two expressing serotype 6B was seen. Re-isolation of either the same or different strains of serotype 6B was seen. Genetically related single-locus and double-locus variants were identified in our cohort that differed by multiple nucleotides. Multiple colony characterisation revealed phenotypic and genetic evidence of microevolution and a greater diversity of pneumococcal strains colonising together than previously observed, thus increasing the potential to adapt in response to the host environment during colonisation.

Conserved mutations in the pneumococcal bacteriocin transporter gene, blpA, result in a complex population consisting of producers and cheaters

All fully sequenced strains of Streptococcus pneumoniae possess a version of the blp locus, which is responsible for bacteriocin production and immunity. Activation of the blp locus is stimulated by accumulation of the peptide pheromone, BlpC, following its secretion by the ABC transporter, BlpA. The blp locus is characterized by significant diversity in blpC type and in the region of the locus containing putative bacteriocin and immunity genes. In addition, the blpA gene can represent a single large open reading frame or be divided into several smaller fragments due to the presence of frameshift mutations. In this study, we use a collection of strains with blp-dependent inhibition and immunity to define the genetic changes that bring about phenotypic differences in bacteriocin production or immunity. We demonstrate that alterations in blpA, blpC, and bacteriocin/immunity content likely play an important role in competitive interactions between pneumococcal strains. Importantly, strains with a highly conserved frameshift mutation in blpA are unable to secrete bacteriocins or BlpC, but retain the ability to respond to exogenous peptide pheromone produced by cocolonizing strains, stimulating blp-mediated immunity. These “cheater” strains can only coexist with bacteriocin-producing strains that secrete their cognate BlpC and share the same immunity proteins. The variable outcome of these interactions helps to explain the heterogeneity of the blp pheromone, bacteriocin, and immunity protein content.

Streptococcus pneumoniae resides in a polymicrobial environment and competes for limited resources by the elaboration of small antimicrobial peptides called bacteriocins. A conserved cluster of genes in the S. pneumoniae genome is involved in the production of bacteriocins and their associated protective immunity proteins through secretion of a signaling pheromone. In this study, we show that a significant number of strains have lost the ability to secrete bacteriocins and signaling pheromones due to a specific mutation in a dedicated transporter protein. Because the regulatory and immunity portion of the locus is retained, these “cheater” strains can survive in the face of invasion from a bacteriocin-producing strain without the cost of bacteriocin secretion. The outcome of such interactions depends on each strain’s repertoire of pheromone, immunity protein, and bacteriocin genes, such that intrastrain competition drives the diversity in bacteriocin, immunity protein, and pheromone content.

Transformation of natural genetic variation into Haemophilus influenzae genomes

Many bacteria are able to efficiently bind and take up double-stranded DNA fragments, and the resulting natural transformation shapes bacterial genomes, transmits antibiotic resistance, and allows escape from immune surveillance. The genomes of many competent pathogens show evidence of extensive historical recombination between lineages, but the actual recombination events have not been well characterized. We used DNA from a clinical isolate of Haemophilus influenzae to transform competent cells of a laboratory strain. To identify which of the ∼40,000 polymorphic differences had recombined into the genomes of four transformed clones, their genomes and their donor and recipient parents were deep sequenced to high coverage. Each clone was found to contain ∼1000 donor polymorphisms in 3–6 contiguous runs (8.1±4.5 kb in length) that collectively comprised ∼1–3% of each transformed chromosome. Seven donor-specific insertions and deletions were also acquired as parts of larger donor segments, but the presence of other structural variation flanking 12 of 32 recombination breakpoints suggested that these often disrupt the progress of recombination events. This is the first genome-wide analysis of chromosomes directly transformed with DNA from a divergent genotype, connecting experimental studies of transformation with the high levels of natural genetic variation found in isolates of the same species.

The ability of bacteria to acquire genetic information from their relatives—called natural competence—poses a major health risk, since recombination between pathogenic bacterial lineages can help bacteria develop resistance to antibiotics and adapt to host defenses. In this study we transformed competent cells of the human pathogen Haemophilus influenzae with genomic DNA from a divergent clinical isolate and used deep sequencing to identify the recombination events in four transformed chromosomes. The results show that transformation of single competent cells is more extensive than expected, and suggests that transformation can be used as a tool to map traits that vary between clinical isolates.

Designing the next generation of vaccines for global public health

Vaccine research and development are experiencing a renaissance of interest from the global scientific community. There are four major reasons for this: (1) the lack of efficacious treatment for many devastating infections; (2) the emergence of multidrug resistant bacteria; (3) the need for improving the safety of the more traditional licensed vaccines; and finally, (4) the great promise for innovative vaccine design and research with convergence of omics sciences, such as genomics, proteomics, immunomics, and vaccinology. Our first project based on omics was initiated in 2000 and was termed reverse vaccinology. At that time, antigen identification was mainly based on bioinformatic analysis of a singular genome. Since then, omics-guided approaches have been applied to its full potential in several proof-of-concept studies in the industry, with the first reverse vaccinology-derived vaccine now in late stage clinical trials and several vaccines developed by omics in preclinical studies. In the meantime, vaccine discovery and development has been further improved with the support of proteomics, functional genomics, comparative genomics, structural biology, and most recently vaccinomics. We illustrate in this review how omics biotechnologies and integrative biology are expected to accelerate the identification of vaccine candidates against difficult pathogens for which traditional vaccine development has thus far been failing, and how research will provide safer vaccines and improved formulations for immunocompromised patients in the near future. Finally, we present a discussion to situate omics-guided rational vaccine design in the broader context of global public health and how it can benefit citizens in both developed and developing countries.

Cell surface proteins in S. pneumoniae, S. mitis and S. oralis

BACKGROUND AND OBJECTIVES

Streptococcus pneumoniae, a major human pathogen, is closely related to the commensal species S. mitis and S. oralis. S. pneumoniae surface proteins are implicated in virulence and host interaction of this species, but many of them have recently been detected in S. mitis B6 in silico. We tested for the presence of such genes usinga set of eight S. mitis and eleven S. oralis strains from different geographic locations.

MATERIALS AND METHODS

An oligonucleotide microarray was designed based on the genomes of S. pneumoniae R6 and TIGR4 as well as S. mitis B6 to include 63 cell surface proteins. The S. pneumoniae genes encoding neuraminidases, hyaluronidase and pneumolysin were also included. In addition to comparative genomic hybridization experiments, homologues were identified in silico in the genome of S. oralis Uo5.

RESULTS AND CONCLUSIONS

The results document that many S. pneumoniae related surface proteins are ubiquitously present among the Mitis group of streptococci. All 19 samples hybridized with the pavA probe representing a gene important for adherence and invasion of S. pneumoniae. Only eight genes were not recognized in any strain, including the S. pneumoniae PcpC gene as the only virulence gene of the S. pneumoniae core genome.The fact that only 12 out of 26 genes present in the S. oralis Uo5 genome could be detected by microarray analysis confirms the sequence variation of surface components.

A variable region within the genome of Streptococcus pneumoniae contributes to strain-strain variation in virulence

The bacterial factors responsible for the variation in invasive potential between different clones and serotypes of Streptococcus pneumoniae are largely unknown. Therefore, the isolation of rare serotype 1 carriage strains in Indigenous Australian communities provided a unique opportunity to compare the genomes of non-invasive and invasive isolates of the same serotype in order to identify such factors. The human virulence status of non-invasive, intermediately virulent and highly virulent serotype 1 isolates was reflected in mice and showed that whilst both human non-invasive and highly virulent isolates were able to colonize the murine nasopharynx equally, only the human highly virulent isolates were able to invade and survive in the murine lungs and blood. Genomic sequencing comparisons between these isolates identified 8 regions >1 kb in size that were specific to only the highly virulent isolates, and included a version of the pneumococcal pathogenicity island 1 variable region (PPI-1v), phage-associated adherence factors, transporters and metabolic enzymes. In particular, a phage-associated endolysin, a putative iron/lead permease and an operon within PPI-1v exhibited niche-specific changes in expression that suggest important roles for these genes in the lungs and blood. Moreover, in vivo competition between pneumococci carrying PPI-1v derivatives representing the two identified versions of the region showed that the version of PPI-1v in the highly virulent isolates was more competitive than the version from the less virulent isolates in the nasopharyngeal tissue, blood and lungs. This study is the first to perform genomic comparisons between serotype 1 isolates with distinct virulence profiles that correlate between mice and humans, and has highlighted the important role that hypervariable genomic loci, such as PPI-1v, play in pneumococcal disease. The findings of this study have important implications for understanding the processes that drive progression from colonization to invasive disease and will help direct the development of novel therapeutic strategies.

Streptococcus pneumoniae clonal complex 199: genetic diversity and tissue-specific virulence

Streptococcus pneumoniae is an important cause of otitis media and invasive disease. Since introduction of the heptavalent pneumococcal conjugate vaccine, there has been an increase in replacement disease due to serotype 19A clonal complex (CC)199 isolates. The goals of this study were to 1) describe genetic diversity among nineteen CC199 isolates from carriage, middle ear, blood, and cerebrospinal fluid, 2) compare CC199 19A (n = 3) and 15B/C (n = 2) isolates in the chinchilla model for pneumococcal disease, and 3) identify accessory genes associated with tissue-specific disease among a larger collection of S. pneumoniae isolates. CC199 isolates were analyzed by comparative genome hybridization. One hundred and twenty-seven genes were variably present. The CC199 phylogeny split into two main clades, one comprised predominantly of carriage isolates and another of disease isolates. Ability to colonize and cause disease did not differ by serotype in the chinchilla model. However, isolates from the disease clade were associated with faster time to bacteremia compared to carriage clade isolates. One 19A isolate exhibited hypervirulence. Twelve tissue-specific genes/regions were identified by correspondence analysis. After screening a diverse collection of 326 isolates, spr0282 was associated with carriage. Four genes/regions, SP0163, SP0463, SPN05002 and RD8a were associated with middle ear isolates. SPN05002 also associated with blood and CSF, while RD8a associated with blood isolates. The hypervirulent isolate's genome was sequenced using the Solexa paired-end sequencing platform and compared to that of a reference serotype 19A isolate, revealing the presence of a novel 20 kb region with sequence similarity to bacteriophage genes. Genetic factors other than serotype may modulate virulence potential in CC199. These studies have implications for the long-term effectiveness of conjugate vaccines. Ideally, future vaccines would target common proteins to effectively reduce carriage and disease in the vaccinated population.

Comparative supragenomic analyses among the pathogens Staphylococcus aureus, Streptococcus pneumoniae, and Haemophilus influenzae using a modification of the finite supragenome model

BACKGROUND

Staphylococcus aureus is associated with a spectrum of symbiotic relationships with its human host from carriage to sepsis and is frequently associated with nosocomial and community-acquired infections, thus the differential gene content among strains is of interest.

RESULTS

We sequenced three clinical strains and combined these data with 13 publically available human isolates and one bovine strain for comparative genomic analyses. All genomes were annotated using RAST, and then their gene similarities and differences were delineated. Gene clustering yielded 3,155 orthologous gene clusters, of which 2,266 were core, 755 were distributed, and 134 were unique. Individual genomes contained between 2,524 and 2,648 genes. Gene-content comparisons among all possible S. aureus strain pairs (n = 136) revealed a mean difference of 296 genes and a maximum difference of 476 genes. We developed a revised version of our finite supragenome model to estimate the size of the S. aureus supragenome (3,221 genes, with 2,245 core genes), and compared it with those of Haemophilus influenzae and Streptococcus pneumoniae. There was excellent agreement between RAST's annotations and our CDS clustering procedure providing for high fidelity metabolomic subsystem analyses to extend our comparative genomic characterization of these strains.

CONCLUSIONS

Using a multi-species comparative supragenomic analysis enabled by an improved version of our finite supragenome model we provide data and an interpretation explaining the relatively larger core genome of S. aureus compared to other opportunistic nasopharyngeal pathogens. In addition, we provide independent validation for the efficiency and effectiveness of our orthologous gene clustering algorithm.

Proteome and membrane fatty acid analyses on Oligotropha carboxidovorans OM5 grown under chemolithoautotrophic and heterotrophic conditions

Oligotropha carboxidovorans OM5 T. (DSM 1227, ATCC 49405) is a chemolithoautotrophic bacterium able to utilize CO and H₂ to derive energy for fixation of CO₂. Thus, it is capable of growth using syngas, which is a mixture of varying amounts of CO and H₂ generated by organic waste gasification. O. carboxidovorans is capable also of heterotrophic growth in standard bacteriologic media. Here we characterize how the O. carboxidovorans proteome adapts to different lifestyles of chemolithoautotrophy and heterotrophy. Fatty acid methyl ester (FAME) analysis of O. carboxidovorans grown with acetate or with syngas showed that the bacterium changes membrane fatty acid composition. Quantitative shotgun proteomic analysis of O. carboxidovorans grown in the presence of acetate and syngas showed production of proteins encoded on the megaplasmid for assimilating CO and H₂ as well as proteins encoded on the chromosome that might have contributed to fatty acid and acetate metabolism. We found that adaptation to chemolithoautotrophic growth involved adaptations in cell envelope, oxidative homeostasis, and metabolic pathways such as glyoxylate shunt and amino acid/cofactor biosynthetic enzymes.

The role of biofilms: are we hitting the right target?

BACKGROUND

Chronic infections affect 17 million people yearly, and approximately 550,000 people die each year from, or with, their chronic infections. Acute and chornic infection differences are well known to clinicians, but the role of bacteria in producing these clinical differences remains poorly understood.

METHODS

This review relies on basic science, clinical studies, and a general review of the medical biofilm literature. The basic science studies are level A and B quality of evidence. The clinical studies are mainly retrospective cohort (level B) and case studies (level C). The biofilm literature includes reviews with varying levels of evidence. All articles have been peer reviewed and meet the standard of evidence-based medicine.

RESULTS

Acute infections are associated with planktonic bacteria and must be diagnosed rapidly and accurately to prevent tissue damage and/or death. In contrast, biofilm behavior pursues a more parasitic course by producing sustained host hyperinflammation, with the biofilm feeding on plasma exudate. Chronic infections vacillate over long periods of time, responding only partially to antibiotics and reemerging once the antibiotics are withdrawn. Chronic wounds exhibit similar clinical behavior seen in other chronic infections and are associated with biofilm phenotype bacteria on their surface. Biofilm infections, such as chronic wounds, cannot be adequately diagnosed with current clinical cultures; therefore, molecular methods are necessary.

CONCLUSIONS

Biofilm phenotype bacteria require multiple concurrent strategies, including débridement and targeted antibiofilm agents. Biofilm phenotype bacteria predominate on the surface of wounds, and biofilm-based management improves wound healing outcomes, indicating that biofilm is the right target for managing the bioburden barrier of chronic wounds.

Application of nanomaterials to arrays for infectious disease diagnosis

As our understanding of infectious agents increases it has become possible to identify the biological components of individual strains with greater accuracy. This ability is aiding clinical microbiologists to control both existing and emerging infectious disease problems. Microarray technologies are playing an increasingly important role in supplying relevant test data. Control of the nanoscale structures of the materials used in arrays, including the supports, reagents and analytes, is now essential. New techniques for the manufacture of nanostructured supports is now allowing rapid advances in the fields of multianalyte testing and high-volume DNA sequencing.

Complete genome sequence of a serotype 11A, ST62 Streptococcus pneumoniae invasive isolate

Background

Streptococcus pneumoniae is an important human pathogen representing a major cause of morbidity and mortality worldwide. We sequenced the genome of a serotype 11A, ST62 S. pneumoniae invasive isolate (AP200), that was erythromycin-resistant due to the presence of the erm(TR) determinant, and carried out analysis of the genome organization and comparison with other pneumococcal genomes.

Results

The genome sequence of S. pneumoniae AP200 is 2,130,580 base pair in length. The genome carries 2216 coding sequences (CDS), 56 tRNA, and 12 rRNA genes. Of the CDSs, 72.9% have a predicted biological known function. AP200 contains the pilus islet 2 and, although its phenotype corresponds to serotype 11A, it contains an 11D capsular locus. Chromosomal rearrangements resulting from a large inversion across the replication axis, and horizontal gene transfer events were observed. The chromosomal inversion is likely implicated in the rebalance of the chromosomal architecture affected by the insertions of two large exogenous elements, the erm(TR)-carrying Tn1806 and a functional prophage designated ϕSpn_200. Tn1806 is 52,457 bp in size and comprises 49 ORFs. Comparative analysis of Tn1806 revealed the presence of a similar genetic element or part of it in related species such as Streptococcus pyogenes and also in the anaerobic species Finegoldia magna, Anaerococcus prevotii and Clostridium difficile. The genome of ϕSpn_200 is 35,989 bp in size and is organized in 47 ORFs grouped into five functional modules. Prophages similar to ϕSpn_200 were found in pneumococci and in other streptococcal species, showing a high degree of exchange of functional modules. ϕSpn_200 viral particles have morphologic characteristics typical of the Siphoviridae family and are capable of infecting a pneumococcal recipient strain.

Conclusions

The sequence of S. pneumoniae AP200 chromosome revealed a dynamic genome, characterized by chromosomal rearrangements and horizontal gene transfers. The overall diversity of AP200 is driven mainly by the presence of the exogenous elements Tn1806 and ϕSpn_200 that show large gene exchanges with other genetic elements of different bacterial species. These genetic elements likely provide AP200 with additional genes, such as those conferring antibiotic-resistance, promoting its adaptation to the environment.

Analysis of invasiveness of pneumococcal serotypes and clones circulating in Portugal before widespread use of conjugate vaccines reveals heterogeneous behavior of clones expressing the same serotype

To estimate the invasive disease potential of serotypes and clones circulating in Portugal before extensive use of the seven-valent pneumococcal conjugate vaccine, we analyzed 475 invasive isolates recovered from children and adults and 769 carriage isolates recovered from children between 2001 and 2003. Isolates were serotyped and genotyped by pulsed-field gel electrophoresis, and a selection of isolates were also characterized by multilocus sequence typing. We found that the diversities of serotypes and genotypes of pneumococci responsible for invasive infections and carriage were identical and that most carried clones could also be detected as causes of invasive disease. Their ability to do so, however, varied substantially. Serotypes 1, 3, 4, 5, 7F, 8, 9N, 9L, 12B, 14, 18C, and 20 were found to have an enhanced propensity to cause invasive disease, while serotypes 6A, 6B, 11A, 15B/C, 16F, 19F, 23F, 34, 35F, and 37 were associated with carriage. In addition, significant differences in invasive disease potential between clones sharing the same serotype were found among several serotypes, namely, 3, 6A, 6B, 11A, 14, 19A, 19F, 22F, 23F, 34, and NT. This heterogeneous behavior of the clones was found irrespective of the serotype's overall invasive disease potential. Our results highlight the importance of the genetic background when analyzing the invasive disease potential of certain serotypes and provide an important baseline for its monitoring following conjugate vaccine use. Continuous surveillance should be maintained, and current research should focus on uncovering the genetic determinants that contribute to the heterogeneity of invasive disease potential of clones sharing the same serotype.

Comparative analysis and supragenome modeling of twelve Moraxella catarrhalis clinical isolates

Background

M. catarrhalis is a gram-negative, gamma-proteobacterium and an opportunistic human pathogen associated with otitis media (OM) and exacerbations of chronic obstructive pulmonary disease (COPD). With direct and indirect costs for treating these conditions annually exceeding $33 billion in the United States alone, and nearly ubiquitous resistance to beta-lactam antibiotics among M. catarrhalis clinical isolates, a greater understanding of this pathogen's genome and its variability among isolates is needed.

Results

The genomic sequences of ten geographically and phenotypically diverse clinical isolates of M. catarrhalis were determined and analyzed together with two publicly available genomes. These twelve genomes were subjected to detailed comparative and predictive analyses aimed at characterizing the supragenome and understanding the metabolic and pathogenic potential of this species. A total of 2383 gene clusters were identified, of which 1755 are core with the remaining 628 clusters unevenly distributed among the twelve isolates. These findings are consistent with the distributed genome hypothesis (DGH), which posits that the species genome possesses a far greater number of genes than any single isolate. Multiple and pair-wise whole genome alignments highlight limited chromosomal re-arrangement.

Conclusions

M. catarrhalis gene content and chromosomal organization data, although supportive of the DGH, show modest overall genic diversity. These findings are in stark contrast with the reported heterogeneity of the species as a whole, as wells as to other bacterial pathogens mediating OM and COPD, providing important insight into M. catarrhalis pathogenesis that will aid in the development of novel therapeutic regimens.

Genomic fluidity: an integrative view of gene diversity within microbial populations

Background

The dual concepts of pan and core genomes have been widely adopted as means to assess the distribution of gene families within microbial species and genera. The core genome is the set of genes shared by a group of organisms; the pan genome is the set of all genes seen in any of these organisms. A variety of methods have provided drastically different estimates of the sizes of pan and core genomes from sequenced representatives of the same groups of bacteria.

Results

We use a combination of mathematical, statistical and computational methods to show that current predictions of pan and core genome sizes may have no correspondence to true values. Pan and core genome size estimates are problematic because they depend on the estimation of the occurrence of rare genes and genomes, respectively, which are difficult to estimate precisely because they are rare. Instead, we introduce and evaluate a robust metric - genomic fluidity - to categorize the gene-level similarity among groups of sequenced isolates. Genomic fluidity is a measure of the dissimilarity of genomes evaluated at the gene level.

Conclusions

The genomic fluidity of a population can be estimated accurately given a small number of sequenced genomes. Further, the genomic fluidity of groups of organisms can be compared robustly despite variation in algorithms used to identify genes and their homologs. As such, we recommend that genomic fluidity be used in place of pan and core genome size estimates when assessing gene diversity within genomes of a species or a group of closely related organisms.

SNP discovery by transcriptome pyrosequencing

Single nucleotide polymorphisms (SNPs) are single base differences between haplotypes. SNPs are abundant in many species and valuable as markers for genetic map construction, modern molecular breeding programs, and quantitative genetic studies. SNPs are readily mined from genomic DNA or cDNA sequence obtained from individuals having two or more distinct genotypes. While automated Sanger sequencing has become less expensive over time, it is still costly to acquire deep Sanger sequence from several genotypes. "Next-generation" DNA sequencing technologies that utilize new chemistries and massively parallel approaches have enabled DNA sequences to be acquired at extremely high depths of coverage faster and for less cost than traditional sequencing. One such method is represented by the Roche/454 Life Sciences GS-FLX Titanium Series, which currently uses pyrosequencing to produce up to 400-600 million bases of DNA sequence/run (>1 million reads, ~400 bp/read). This chapter discusses the use of high-throughput pyrosequencing for SNP discovery by focusing on 454 sequencing of maize cDNA, the development of a computational pipeline for polymorphism detection, and the subsequent identification of over 7,000 putative SNPs between Mo17 and B73 maize. In addition, alternative alignment and polymorphism detection strategies that implement Illumina short reads, data processing and visualization tools, and reduced representation techniques that reduce the sequencing of repeat DNA, thus enabling efficient analysis of genome sequence, are discussed.

Omics, Bioinformatics, and Infectious Disease Research

Bioinformatics is basically the study of informatic processes in biotic systems. Actually what constitutes bioinformatics is not entirely clear and arguably varies depending on who tries to define it. This chapter discusses the considerable progress in infectious diseases research that has been made in recent years using various “omics” case studies. Bioinformatics is tasked with making sense of it, mining it, storing it, disseminating it, and ensuring valid biological conclusions can be drawn from it. This chapter discusses the current state of play of bioinformatics related to genomics and transcriptomics, briefs metagenomics that finds use in infectious disease research as well as the random sequencing of genomes from a variety of organisms. This chapter explains the various possibilities of pan-genome, transcriptional reshaping and also enormous progress of proteomics study. Bioinformatic algorithms and tools are crucial tools in analyzing the data. The chapter also attempts to provide some details on the various problems and solution in bioinformatics that current-day scientists face while concentrating on second-generation sequencing strategies.

Analysis of ultra low genome conservation in Clostridium difficile

Microarray-based comparative genome hybridisations (CGH) and genome sequencing of Clostridium difficile isolates have shown that the genomes of this species are highly variable. To further characterize their genome variation, we employed integration of data from CGH, genome sequencing and putative cellular pathways. Transcontinental strain comparison using CGH data confirmed the emergence of a human-specific hypervirulent cluster. However, there was no correlation between total toxin production and hypervirulent phenotype, indicating the possibility of involvement of additional factors towards hypervirulence. Calculation of C. difficile core and pan genome size using CGH and sequence data estimated that the core genome is composed of 947 to 1,033 genes and a pan genome comprised of 9,640 genes. The reconstruction, annotation and analysis of cellular pathways revealed highly conserved pathways despite large genome variation. However, few pathways such as tetrahydrofolate biosynthesis were found to be variable and could be contributing to adaptation towards virulence such as antibiotic resistance.

Structure and dynamics of the pan-genome of Streptococcus pneumoniae and closely related species

Background

Streptococcus pneumoniae is one of the most important causes of microbial diseases in humans. The genomes of 44 diverse strains of S. pneumoniae were analyzed and compared with strains of non-pathogenic streptococci of the Mitis group.

Results

Despite evidence of extensive recombination, the S. pneumoniae phylogenetic tree revealed six major lineages. With the exception of serotype 1, the tree correlated poorly with capsular serotype, geographical site of isolation and disease outcome. The distribution of dispensable genes - genes present in more than one strain but not in all strains - was consistent with phylogeny, although horizontal gene transfer events attenuated this correlation in the case of ancient lineages. Homologous recombination, involving short stretches of DNA, was the dominant evolutionary process of the core genome of S. pneumoniae. Genetic exchange occurred both within and across the borders of the species, and S. mitis was the main reservoir of genetic diversity of S. pneumoniae. The pan-genome size of S. pneumoniae increased logarithmically with the number of strains and linearly with the number of polymorphic sites of the sampled genomes, suggesting that acquired genes accumulate proportionately to the age of clones. Most genes associated with pathogenicity were shared by all S. pneumoniae strains, but were also present in S. mitis, S. oralis and S. infantis, indicating that these genes are not sufficient to determine virulence.

Conclusions

Genetic exchange with related species sharing the same ecological niche is the main mechanism of evolution of S. pneumoniae. The open pan-genome guarantees the species a quick and economical response to diverse environments.

Pediatric vaccines on the horizon

Vaccines have saved the lives of millions of children and continue to be essential interventions to control infectious diseases among people of all ages. The list of recommended vaccines for children has expanded in recent years; however, many viral, bacterial and parasitic infections remain a major cause of morbidity and mortality in children. Improved vaccines to prevent Streptococcus pneumoniae and Neisseria meningitidis infections in children will soon be available. Recent scientific advances are being applied to design new childhood vaccines affording enhanced efficacy, safety and tolerability. Financial barriers and other obstacles to adequate vaccine access need to be eliminated to assure coverage for all children and adolescents.

SOLiD sequencing of four Vibrio vulnificus genomes enables comparative genomic analysis and identification of candidate clade-specific virulence genes

BACKGROUND

Vibrio vulnificus is the leading cause of reported death from consumption of seafood in the United States. Despite several decades of research on molecular pathogenesis, much remains to be learned about the mechanisms of virulence of this opportunistic bacterial pathogen. The two complete and annotated genomic DNA sequences of V. vulnificus belong to strains of clade 2, which is the predominant clade among clinical strains. Clade 2 strains generally possess higher virulence potential in animal models of disease compared with clade 1, which predominates among environmental strains. SOLiD sequencing of four V. vulnificus strains representing different clades (1 and 2) and biotypes (1 and 2) was used for comparative genomic analysis.

RESULTS

Greater than 4,100,000 bases were sequenced of each strain, yielding approximately 100-fold coverage for each of the four genomes. Although the read lengths of SOLiD genomic sequencing were only 35 nt, we were able to make significant conclusions about the unique and shared sequences among the genomes, including identification of single nucleotide polymorphisms. Comparative analysis of the newly sequenced genomes to the existing reference genomes enabled the identification of 3,459 core V. vulnificus genes shared among all six strains and 80 clade 2-specific genes. We identified 523,161 SNPs among the six genomes.

CONCLUSIONS

We were able to glean much information about the genomic content of each strain using next generation sequencing. Flp pili, GGDEF proteins, and genomic island XII were identified as possible virulence factors because of their presence in virulent sequenced strains. Genomic comparisons also point toward the involvement of sialic acid catabolism in pathogenesis.

Complete genome and comparative analysis of the chemolithoautotrophic bacterium Oligotropha carboxidovorans OM5

Background

Oligotropha carboxidovorans OM5 T. (DSM 1227, ATCC 49405) is a chemolithoautotrophic bacterium capable of utilizing CO (carbon monoxide) and fixing CO₂ (carbon dioxide). We previously published the draft genome of this organism and recently submitted the complete genome sequence to GenBank.

Results

The genome sequence of the chemolithoautotrophic bacterium Oligotropha carboxidovorans OM5 consists of a 3.74-Mb chromosome and a 133-kb megaplasmid that contains the genes responsible for utilization of carbon monoxide, carbon dioxide, and hydrogen. To our knowledge, this strain is the first one to be sequenced in the genus Oligotropha, the closest fully sequenced relatives being Bradyrhizobium sp. BTAi and USDA110 and Nitrobacter hamburgiensis X14. Analysis of the O. carboxidovorans genome reveals potential links between plasmid-encoded chemolithoautotrophy and chromosomally-encoded lipid metabolism. Comparative analysis of O. carboxidovorans with closely related species revealed differences in metabolic pathways, particularly in carbohydrate and lipid metabolism, as well as transport pathways.

Conclusion

Oligotropha, Bradyrhizobium sp and Nitrobacter hamburgiensis X14 are phylogenetically proximal. Although there is significant conservation of genome organization between the species, there are major differences in many metabolic pathways that reflect the adaptive strategies unique to each species.

Probing the pan-genome of Listeria monocytogenes: new insights into intraspecific niche expansion and genomic diversification

Background

Bacterial pathogens often show significant intraspecific variations in ecological fitness, host preference and pathogenic potential to cause infectious disease. The species of Listeria monocytogenes, a facultative intracellular pathogen and the causative agent of human listeriosis, consists of at least three distinct genetic lineages. Two of these lineages predominantly cause human sporadic and epidemic infections, whereas the third lineage has never been implicated in human disease outbreaks despite its overall conservation of many known virulence factors.

Results

Here we compare the genomes of 26 L. monocytogenes strains representing the three lineages based on both in silico comparative genomic analysis and high-density, pan-genomic DNA array hybridizations. We uncover 86 genes and 8 small regulatory RNAs that likely make L. monocytogenes lineages differ in carbohydrate utilization and stress resistance during their residence in natural habitats and passage through the host gastrointestinal tract. We also identify 2,330 to 2,456 core genes that define this species along with an open pan-genome pool that contains more than 4,052 genes. Phylogenomic reconstructions based on 3,560 homologous groups allowed robust estimation of phylogenetic relatedness among L. monocytogenes strains.

Conclusions

Our pan-genome approach enables accurate co-analysis of DNA sequence and hybridization array data for both core gene estimation and phylogenomics. Application of our method to the pan-genome of L. monocytogenes sheds new insights into the intraspecific niche expansion and evolution of this important foodborne pathogen.

Generation of genic diversity among Streptococcus pneumoniae strains via horizontal gene transfer during a chronic polyclonal pediatric infection

Although there is tremendous interest in understanding the evolutionary roles of horizontal gene transfer (HGT) processes that occur during chronic polyclonal infections, to date there have been few studies that directly address this topic. We have characterized multiple HGT events that most likely occurred during polyclonal infection among nasopharyngeal strains of Streptococcus pneumoniae recovered from a child suffering from chronic upper respiratory and middle-ear infections. Whole genome sequencing and comparative genomics were performed on six isolates collected during symptomatic episodes over a period of seven months. From these comparisons we determined that five of the isolates were genetically highly similar and likely represented a dominant lineage. We analyzed all genic and allelic differences among all six isolates and found that all differences tended to occur within contiguous genomic blocks, suggestive of strain evolution by homologous recombination. From these analyses we identified three strains (two of which were recovered on two different occasions) that appear to have been derived sequentially, one from the next, each by multiple recombination events. We also identified a fourth strain that contains many of the genomic segments that differentiate the three highly related strains from one another, and have hypothesized that this fourth strain may have served as a donor multiple times in the evolution of the dominant strain line. The variations among the parent, daughter, and grand-daughter recombinant strains collectively cover greater than seven percent of the genome and are grouped into 23 chromosomal clusters. While capturing in vivo HGT, these data support the distributed genome hypothesis and suggest that a single competence event in pneumococci can result in the replacement of DNA at multiple non-adjacent loci.

Bacterial infections have long been studied using Koch's postulates wherein the paradigm is that a single clone leads to a given infection. Over the past decade, it has become clear that chronic bacterial infections often do not fit this paradigm. Instead these are associated with the presence of multiple strains or species (polyclonal) of bacteria that are organized into highly structured communities, termed biofilms, which can persist in the body and are recalcitrant to antibiotic treatment. In addition, there is extensive evidence that bacteria can incorporate genes from neighboring bacteria into their own genomes. This process can produce new strains and is known as horizontal gene transfer. In this study, we investigated for the first time, the tempo and relevance of gene transfer among bacterial strains of Streptococcus pneumoniae during a naturally occurring chronic childhood infection. We identified extensive gene transfer among multiple infecting strains, by sequencing of isolates recovered sequentially over a seven-month period. This gene transfer may serve as a counterpoint to the host's adaptive immune response and help explain the phenomenon of bacterial persistence, since, as occurs with some chronic viral and parasitic infections, the immune system may become overwhelmed by a set of related strains.

Comparative genomic characterization of Actinobacillus pleuropneumoniae

The Gram-negative bacterium Actinobacillus pleuropneumoniae is the etiologic agent of porcine contagious pleuropneumoniae, a lethal respiratory infectious disease causing great economic losses in the swine industry worldwide. In order to better interpret the genetic background of serotypic diversity, nine genomes of A. pleuropneumoniae reference strains of serovars 1, 2, 4, 6, 9, 10, 11, 12, and 13 were sequenced by using rapid high-throughput approach. Based on 12 genomes of corresponding serovar reference strains including three publicly available complete genomes (serovars 3, 5b, and 7) of this bacterium, we performed a comprehensive analysis of comparative genomics and first reported a global genomic characterization for this pathogen. Clustering of 26,012 predicted protein-coding genes showed that the pan genome of A. pleuropneumoniae consists of 3,303 gene clusters, which contain 1,709 core genome genes, 822 distributed genes, and 772 strain-specific genes. The genome components involved in the biogenesis of capsular polysaccharide and lipopolysaccharide O antigen relative to serovar diversity were compared, and their genetic diversity was depicted. Our findings shed more light on genomic features associated with serovar diversity of A. pleuropneumoniae and provide broader insight into both pathogenesis research and clinical/epidemiological application against the severe disease caused by this swine pathogen.

Genomic comparisons of Brucella spp. and closely related bacteria using base compositional and proteome based methods

Background

Classification of bacteria within the genus Brucella has been difficult due in part to considerable genomic homogeneity between the different species and biovars, in spite of clear differences in phenotypes. Therefore, many different methods have been used to assess Brucella taxonomy. In the current work, we examine 32 sequenced genomes from genus Brucella representing the six classical species, as well as more recently described species, using bioinformatical methods. Comparisons were made at the level of genomic DNA using oligonucleotide based methods (Markov chain based genomic signatures, genomic codon and amino acid frequencies based comparisons) and proteomes (all-against-all BLAST protein comparisons and pan-genomic analyses).

Results

We found that the oligonucleotide based methods gave different results compared to that of the proteome based methods. Differences were also found between the oligonucleotide based methods used. Whilst the Markov chain based genomic signatures grouped the different species in genus Brucella according to host preference, the codon and amino acid frequencies based methods reflected small differences between the Brucella species. Only minor differences could be detected between all genera included in this study using the codon and amino acid frequencies based methods.

Proteome comparisons were found to be in strong accordance with current Brucella taxonomy indicating a remarkable association between gene gain or loss on one hand and mutations in marker genes on the other. The proteome based methods found greater similarity between Brucella species and Ochrobactrum species than between species within genus Agrobacterium compared to each other. In other words, proteome comparisons of species within genus Agrobacterium were found to be more diverse than proteome comparisons between species in genus Brucella and genus Ochrobactrum. Pan-genomic analyses indicated that uptake of DNA from outside genus Brucella appears to be limited.

Conclusions

While both the proteome based methods and the Markov chain based genomic signatures were able to reflect environmental diversity between the different species and strains of genus Brucella, the genomic codon and amino acid frequencies based comparisons were not found adequate for such comparisons. The proteome comparison based phylogenies of the species in genus Brucella showed a surprising consistency with current Brucella taxonomy.

The distributed genome hypothesis as a rubric for understanding evolution in situ during chronic bacterial biofilm infectious processes

Most chronic infectious disease processes associated with bacteria are characterized by the formation of a biofilm that provides for bacterial attachment to the host tissue or the implanted medical device. The biofilm protects the bacteria from the host's adaptive immune response as well as predation by phagocytic cells. However, the most insidious aspect of biofilm biology from the host's point of view is that the biofilm provides an ideal setting for bacterial horizontal gene transfer (HGT). HGT provides for large-scale genome content changes in situ during the chronic infectious process. Obviously, for HGT processes to result in the reassortment of alleles and genes among bacterial strains, the infection must be polyclonal (polymicrobial) in nature. In this review, we marshal the evidence that all of the factors are present in biofilm infections to support HGT that results in the ongoing production of novel strains with unique combinations of genic characteristics and that the continual production of large numbers of novel, but related bacterial strains leads to persistence. This concept of an infecting population of bacteria undergoing mutagenesis to produce a 'cloud' of similar strains to confuse and overwhelm the host's immune system parallels genetic diversity strategies used by viral and parasitic pathogens.

Unravelling the ecology of antibiotic resistant bacteria in the nasopharynx

To study the dynamics and diversity of pneumococcal carriage and antibiotic resistance, a more thorough and systematic approach has been employed compared with routine surveillance of serotype and anti-biotic resistance. Up to ten pneumococcal isolates from pernasal (nose) and oropharyngeal (throat) sites are isolated and characterised. Our carriage studies have revealed a diverse community of pneumococci with multiple strains colonising the nasopharynx of children. In Tanzanian children less than 6 years of age, up to six serotypes and up to six different antibiotic sensitivities (as distinguished by at least a fourfold difference in the minimum inhibitory concentration) have been found. Serotyping by the Quelling reaction is prone to inaccuracy and requires expensive serological reagents. To improve the accuracy and reduce the costs, an alternative capsular typing DNA-based method has been developed. This chapter will describe the methods we have employed with emphasis on the capsular typing method.