Monitoring the long-term molecular epidemiology of the pneumococcus and detection of potential 'vaccine escape' strains

Background

While the pneumococcal protein conjugate vaccines reduce the incidence in invasive pneumococcal disease (IPD), serotype replacement remains a major concern. Thus, serotype-independent protection with vaccines targeting virulence genes, such as PspA, have been pursued. PspA is comprised of diverse clades that arose through recombination. Therefore, multi-locus sequence typing (MLST)-defined clones could conceivably include strains from multiple PspA clades. As a result, a method is needed which can both monitor the long-term epidemiology of the pneumococcus among a large number of isolates, and analyze vaccine-candidate genes, such as pspA, for mutations and recombination events that could result in ‘vaccine escape’ strains.

Methodology

We developed a resequencing array consisting of five conserved and six variable genes to characterize 72 pneumococcal strains. The phylogenetic analysis of the 11 concatenated genes was performed with the MrBayes program, the single nucleotide polymorphism (SNP) analysis with the DNA Sequence Polymorphism program (DnaSP), and the recombination event analysis with the recombination detection package (RDP).

Results

The phylogenetic analysis correlated with MLST, and identified clonal strains with unique PspA clades. The DnaSP analysis correlated with the serotype-specific diversity detected using MLST. Serotypes associated with more than one ST complex had a larger degree of sequence polymorphism than a serotype associated with one ST complex. The RDP analysis confirmed the high frequency of recombination events in the pspA gene.

Conclusions

The phylogenetic tree correlated with MLST, and detected multiple PspA clades among clonal strains. The genetic diversity of the strains and the frequency of recombination events in the mosaic gene, pspA were accurately assessed using the DnaSP and RDP programs, respectively. These data provide proof-of-concept that resequencing arrays could play an important role within research and clinical laboratories in both monitoring the molecular epidemiology of the pneumococcus and detecting ‘vaccine escape’ strains among vaccine-candidate genes.     

Microdiversity of extracellular enzyme genes among sequenced prokaryotic genomes

Understanding the relationship between prokaryotic traits and phylogeny is important for predicting and modeling ecological processes. Microbial extracellular enzymes have a pivotal role in nutrient cycling and the decomposition of organic matter, yet little is known about the phylogenetic distribution of genes encoding these enzymes. In this study, we analyzed 3058 annotated prokaryotic genomes to determine which taxa have the genetic potential to produce alkaline phosphatase, chitinase and β-N-acetyl-glucosaminidase enzymes. We then evaluated the relationship between the genetic potential for enzyme production and 16S rRNA phylogeny using the consenTRAIT algorithm, which calculated the phylogenetic depth and corresponding 16S rRNA sequence identity of clades of potential enzyme producers. Nearly half (49.2%) of the genomes analyzed were found to be capable of extracellular enzyme production, and these were non-randomly distributed across most prokaryotic phyla. On average, clades of potential enzyme-producing organisms had a maximum phylogenetic depth of 0.008004–0.009780, though individual clades varied broadly in both size and depth. These values correspond to a minimum 16S rRNA sequence identity of 98.04–98.40%. The distribution pattern we found is an indication of microdiversity, the occurrence of ecologically or physiologically distinct populations within phylogenetically related groups. Additionally, we found positive correlations among the genes encoding different extracellular enzymes. Our results suggest that the capacity to produce extracellular enzymes varies at relatively fine-scale phylogenetic resolution. This variation is consistent with other traits that require a small number of genes and provides insight into the relationship between taxonomy and traits that may be useful for predicting ecological function.     

Overlapped Sequence Types (STs) and Serogroups of Avian Pathogenic (APEC) and Human Extra-Intestinal Pathogenic (ExPEC) Escherichia coli Isolated in Brazil

Avian pathogenic Escherichia coli (APEC) strains belong to a category that is associated with colibacillosis, a serious illness in the poultry industry worldwide. Additionally, some APEC groups have recently been described as potential zoonotic agents. In this work, we compared APEC strains with extraintestinal pathogenic E. coli (ExPEC) strains isolated from clinical cases of humans with extra-intestinal diseases such as urinary tract infections (UTI) and bacteremia. PCR results showed that genes usually found in the ColV plasmid (tsh, iucA, iss, and hlyF) were associated with APEC strains while fyuA, irp-2, fepC sitD_chrom, fimH, crl, csgA, afa, iha, sat, hlyA, hra, cnf1, kpsMTII, clpV _Sakai and malX were associated with human ExPEC. Both categories shared nine serogroups (O2, O6, O7, O8, O11, O19, O25, O73 and O153) and seven sequence types (ST10, ST88, ST93, ST117, ST131, ST155, ST359, ST648 and ST1011). Interestingly, ST95, which is associated with the zoonotic potential of APEC and is spread in avian E. coli of North America and Europe, was not detected among 76 APEC strains. When the strains were clustered based on the presence of virulence genes, most ExPEC strains (71.7%) were contained in one cluster while most APEC strains (63.2%) segregated to another. In general, the strains showed distinct genetic and fingerprint patterns, but avian and human strains of ST359, or ST23 clonal complex (CC), presented more than 70% of similarity by PFGE. The results demonstrate that some “zoonotic-related” STs (ST117, ST131, ST10CC, ST23CC) are present in Brazil. Also, the presence of moderate fingerprint similarities between ST359 E. coli of avian and human origin indicates that strains of this ST are candidates for having zoonotic potential.     

Intestine and Environment of the Chicken as Reservoirs for Extraintestinal Pathogenic Escherichia coli Strains with Zoonotic Potential

Although research has increasingly focused on the pathogenesis of avian pathogenic Escherichia coli (APEC) infections and the “APEC pathotype” itself, little is known about the reservoirs of these bacteria. We therefore compared outbreak strains isolated from diseased chickens (n = 121) with nonoutbreak strains, including fecal E. coli strains from clinically healthy chickens (n = 211) and strains from their environment (n = 35) by determining their virulence gene profiles, phylogenetic backgrounds, responses to chicken serum, and in vivo pathogenicities in a chicken infection model. In general, by examining 46 different virulence-associated genes we were able to distinguish the three groups of avian strains, but some specific fecal and environmental isolates had a virulence gene profile that was indistinguishable from that determined for outbreak strains. In addition, a substantial number of phylogenetic EcoR group B2 strains, which are known to include potent human and animal extraintestinal pathogenic E. coli (ExPEC) strains, were identified among the APEC strains (44.5%) as well as among the fecal E. coli strains from clinically healthy chickens (23.2%). Comparably high percentages (79.2 to 89.3%) of serum-resistant strains were identified for all three groups of strains tested, bringing into question the usefulness of this phenotype as a principal marker for extraintestinal virulence. Intratracheal infection of 5-week-old chickens corroborated the pathogenicity of a number of nonoutbreak strains. Multilocus sequence typing data revealed that most strains that were virulent in chicken infection experiments belonged to sequence types that are almost exclusively associated with extraintestinal diseases not only in birds but also in humans, like septicemia, urinary tract infection, and newborn meningitis, supporting the hypothesis that not the ecohabitat but the phylogeny of E. coli strains determines virulence. These data provide strong evidence for an avian intestinal reservoir hypothesis which could be used to develop intestinal intervention strategies. These strains pose a zoonotic risk because either they could be transferred directly from birds to humans or they could serve as a genetic pool for ExPEC strains.     

Population Structure and Antimicrobial Resistance Profiles of Streptococcus suis Serotype 2 Sequence Type 25 Strains

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

Multi-locus sequence typing of enteroaggregative Escherichia coli isolates from Nigerian children uncovers multiple lineages

Background

Enteroaggregative Escherichia coli (EAEC) are defined by their stacked-brick adherence pattern to human epithelial cells. There is no all-encompassing genetic marker for EAEC. The category is commonly implicated in diarrhea but research is hampered by perplexing heterogeneity.

Methodology/Principal Findings

To identify key EAEC lineages, we applied multilocus sequence typing to 126 E. coli isolates from a Nigerian case-control study that showed aggregative adherence in the HEp-2 adherence assay, and 24 other EAEC strains from diverse locations. EAEC largely belonged to the A, B1 and D phylogenetic groups and only 7 (4.6%) isolates were in the B2 cluster. As many as 96 sequence types (STs) were identified but 60 (40%) of the EAEC strains belong to or are double locus variants of STs 10, 31, and 394. The remainder did not belong to predominant complexes. The most common ST complex, with predicted ancestor ST10, included 32 (21.3%) of the isolates. Significant age-related distribution suggests that weaned children in Nigeria are at risk for diarrhea from of ST10-complex EAEC. Phylogenetic group D EAEC strains, predominantly from ST31- and ST394 complexes, represented 38 (25.3%) of all isolates, include genome-sequenced strain 042, and possessed conserved chromosomal loci.

Conclusions/Significance

We have developed a molecular phylogenetic framework, which demonstrates that although grouped by a shared phenotype, the category of ‘EAEC’ encompasses multiple pathogenic lineages. Principal among isolates from Nigeria were ST10-complex EAEC that were associated with diarrhea in children over one year and ECOR D strains that share horizontally acquired loci.     

The 100 kDa haem:haemopexin-binding protein of Haemophilus Influenzae: structure and localization

All Haemophilus influenzae strains have an absolute requirement for exogenously supplied haem for aerobic growth. A majority of strains of H. influenzae type b (Hib) produce a 100 kDa protein which binds haem: haemopexin complexes. This 100 kDa haem:haemopexin binding protein, designated HxuA, was originally detected on the Hib cell surface. Monoclonal antibody (mAb)-based analyses revealed that the HxuA protein was also present in soluble form in Hib culture supernatants. This soluble HxuA protein exhibited haem:haemopexin-binding activity in a direct binding assay. Nucleotide sequence analysis of the hxuA gene from Hib strain DL42, together with N-terminal amino acid analysis of HxuA protein purified from Hib culture supernatant, revealed that this protein was synthesized as a 101 kDa precursor with a leader peptide that was removed to yield a 99kDa protein. Southern blot analysis of chromosomal DNA from four Hib and four non-typeable H. influenzae (NTHI) strains detected the presence of a single band in each strain that hybridized a Hib hxuA gene probe. Subsequent analysis of these NTHI strains showed that all four strains released into culture supernatant a haem:haemopexin-binding protein that migrated in SDS-PAGE at a rate similar or identical to that of the Hib HxuA protein. A Hib hxuA mutant was used to screen an NTHI genomic DNA library and an NTHI gene was cloned that complemented the mutation in this Hib strain. Nucleotide sequence analysis of this NTHI gene revealed that it encoded a protein with 87% identity to the Hib HxuA protein. The expression of HxuA by both Hib and NTHI strains indicates that this particular haem acquisition system is conserved among H. influenzae strains.     

The Tryptophanase Gene Cluster of Haemophilus influenzae Type b: Evidence for Horizontal Gene Transfer

Among strains of Haemophilus influenzae, the ability to catabolize tryptophan (as detected by indole production) varies and is correlated with pathogenicity. Tryptophan catabolism is widespread (70 to 75%) among harmless respiratory isolates but is nearly universal (94 to 100%) among strains causing serious disease, including meningitis. As a first step in investigating the relationship between tryptophan catabolism and virulence, we have identified genes in pathogenic H. influenzae which are homologous to the tryptophanase (tna) operon of Escherichia coli. The tna genes are located on a 3.1-kb fragment between nlpD and mutS in the H. influenzae type b (Eagan) genome, are flanked by 43-bp direct repeats of an uptake signal sequence downstream from nlpD, and appear to have been inserted as a mobile unit within this sequence. The organization of this insertion is reminiscent of pathogenicity islands. The tna cluster is found at the same map location in all indole-positive strains of H. influenzae surveyed and is absent from reference type d and e genomes. In contrast to H. influenzae, most other Haemophilus species lack tna genes. Phylogenetic comparisons suggest that the tna cluster was acquired by intergeneric lateral transfer, either by H. influenzae or a recent ancestor, and that E. coli may have acquired its tnaA gene from a related source. Genomes of virulent H. influenzae resemble those of pathogenic enterics in having an island of laterally transferred DNA next to mutS.     

Identification of a second family of high-molecular-weight adhesion proteins expressed by non-typable Haemophilus influenzae

We previously reported that two surface-exposed high-molecular-weight proteins, HMW1 and HMW2, expressed by a prototypic strain of non-typable Haemophilus influenzae (NTHI), mediate attachment to human epithelial cells. These proteins are members of a family of highly immunogenic proteins common to 70–75% of NTHI strains. NTHI strains that lack HMW1/ HMW2-like proteins remain capable of efficient attachment to cultured human epithelial cells, suggesting the existence of additional adhesion molecules. We reasoned that characterization of high-molecular-weight immunogenic proteins from an HMW1/HMW2-deficient strain might identify additional adhesion proteins. A genomic library was prepared in λEMBL3 with chromosomal DNA from non-typable Haemophilus strain 11, a strain that lacks HMW1/HMW2-like proteins. The library was screened immunologically with convalescent serum from a child naturally infected with strain 11, and phage clones expressing high-molecular-weight recombinant proteins were identified by Western blot analysis. One clone was identified that expressed a protein with an apparent molecular mass greater than 200 kDa. Transformation of non-adherent Escherichia coli strain DH5α with plasmids containing the genetic locus encoding this protein gave rise to E. colitransformants that adhered avidly to Chang conjunctival cells. Subcloning and mutagenesis studies localized the DNA conferring the adherence phenotype to a 4.8 kbp fragment, and nucleotide sequence analysis further localized the gene encoding the adhesion protein to a 3.3 kbp open reading frame predicted to encode a protein of 114kDa. The gene was designated hia for Haemophilus influenzae adhesin. Southern analysis revealed an hia homologue in 13 of 15 HMW1/HMW2-deficient non-typable H. influenzae strains. In contrast, the hia gene was not present in any of 23 non-typable H. influenzae strains which expressed HMW1/HMW2-like proteins. Identification of this second family of high-molecular-weight adhesion proteins suggests the possibility of developing vaccines based upon a combination of HMW1/HMW2-like proteins and Hia-like proteins which would be protective against disease caused by most or all non-typable H. influenzae     

Molecular Evolution of Mycoplasma capricolum subsp. capripneumoniae Strains, Based on Polymorphisms in the 16S rRNA Genes

Mycoplasma capricolum subsp. capripneumoniae belongs to the so-called Mycoplasma mycoides cluster and is the causal agent of contagious caprine pleuropneumonia (CCPP). All members of the M. mycoides cluster have two rRNA operons. The sequences of the 16S rRNA genes of both rRNA operons from 20 strains of M. capricolum subsp. capripneumoniae of different geographical origins in Africa and Asia were determined. Nucleotide differences which were present in only one of the two operons (polymorphisms) were detected in 24 positions. The polymorphisms were not randomly distributed in the 16S rRNA genes, and some of them were found in regions of low evolutionary variability. Interestingly, 11 polymorphisms were found in all the M. capricolum subsp. capripneumoniae strains, thus defining a putative ancestor. A sequence length difference between the 16S rRNA genes in a poly(A) region and 12 additional polymorphisms were found in only one or some of the strains. A phylogenetic tree was constructed by comparative analysis of the polymorphisms, and this tree revealed two distinct lines of descent. The nucleotide substitution rate of strains within line II was up to 50% higher than within line I. A tree was also constructed from individual operonal 16S rRNA sequences, and the sequences of the two operons were found to form two distinct clades. The topologies of both clades were strikingly similar, which supports the use of 16S rRNA sequence data from homologous operons for phylogenetic studies. The strain-specific polymorphism patterns of the 16S rRNA genes of M. capricolum subsp. capripneumoniae may be used as epidemiological markers for CCPP.     

Molecular Characterization of the Recently Emerged Poultry Pathogen Ornithobacterium rhinotracheale by Multilocus Sequence Typing

Ornithobacterium rhinotracheale (ORT) is an economically important bacterial pathogen of turkeys and chickens worldwide. Since its first detection, a variety of typing methods have been used to gain basic knowledge about the bacterial population structure, an issue that still needs to be addressed. Serological characterization revealed at least 18 different serotypes (A-R) with ORT of serotype A to be predominate among poultry. This study aimed to establish a multilocus sequence typing (MLST) scheme for ORT that could easily be used by other laboratories and allows for worldwide comparison of sequence data. For this purpose, 87 ORT strains from different poultry hosts, geographical origins, years of isolation and serotypes were included in the analysis to identify correlations. Fourteen different sequence types (ST) were found. The most common ST1 was identified in 40 ORT strains from turkeys and chickens on 4 continents and in 3 different European countries. Together with ST9, both STs represented over three quarters (77%) of ORT strains used in the MLST analysis and included strains of frequently cross-reacting ORT serotypes A, E and I. Nine STs were only represented by one ORT strain and might indicate possible avian host, disease or serotype-specific relationships. In contrast, discrepancies between serotype and phylogenetic relatedness were clearly demonstrated by ORT strains that belonged to identical serotypes but differed in their ST. The overall identified low genetic diversity among strains isolated from turkeys and chickens independent of host and geographical origins suggests that ORT has only recently been introduced into domestic poultry and dispersed worldwide.     

Complex Population Structure and Virulence Differences among Serotype 2 Streptococcus suis Strains Belonging to Sequence Type 28

Streptococcus suis is a major swine pathogen and a zoonotic agent. Serotype 2 strains are the most frequently associated with disease. However, not all serotype 2 lineages are considered virulent. Indeed, sequence type (ST) 28 serotype 2 S. suis strains have been described as a homogeneous group of low virulence. However, ST28 strains are often isolated from diseased swine in some countries, and at least four human ST28 cases have been reported. Here, we used whole-genome sequencing and animal infection models to test the hypothesis that the ST28 lineage comprises strains of different genetic backgrounds and different virulence. We used 50 S. suis ST28 strains isolated in Canada, the United States and Japan from diseased pigs, and one ST28 strain from a human case isolated in Thailand. We report a complex population structure among the 51 ST28 strains. Diversity resulted from variable gene content, recombination events and numerous genome-wide polymorphisms not attributable to recombination. Phylogenetic analysis using core genome single-nucleotide polymorphisms revealed four discrete clades with strong geographic structure, and a fifth clade formed by US, Thai and Japanese strains. When tested in experimental animal models, strains from this latter clade were significantly more virulent than a Canadian ST28 reference strain, and a closely related Canadian strain. Our results highlight the limitations of MLST for both phylogenetic analysis and virulence prediction and raise concerns about the possible emergence of ST28 strains in human clinical cases.     

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

The distributed-genome hypothesis (DGH) states that pathogenic bacteria possess a supragenome that is much larger than the genome of any single bacterium and that these pathogens utilize genetic recombination and a large, noncore set of genes as a means of diversity generation. We sequenced the genomes of eight nasopharyngeal strains of Streptococcus pneumoniae isolated from pediatric patients with upper respiratory symptoms and performed quantitative genomic analyses among these and nine publicly available pneumococcal strains. Coding sequences from all strains were grouped into 3,170 orthologous gene clusters, of which 1,454 (46%) were conserved among all 17 strains. The majority of the gene clusters, 1,716 (54%), were not found in all strains. Genic differences per strain pair ranged from 35 to 629 orthologous clusters, with each strain's genome containing between 21 and 32% noncore genes. The distribution of the orthologous clusters per genome for the 17 strains was entered into the finite-supragenome model, which predicted that (i) the S. pneumoniae supragenome contains more than 5,000 orthologous clusters and (ii) 99% of the orthologous clusters (~3,000) that are represented in the S. pneumoniae population at frequencies of ≥0.1 can be identified if 33 representative genomes are sequenced. These extensive genic diversity data support the DGH and provide a basis for understanding the great differences in clinical phenotype associated with various pneumococcal strains. When these findings are taken together with previous studies that demonstrated the presence of a supragenome for Streptococcus agalactiae and Haemophilus influenzae, it appears that the possession of a distributed genome is a common host interaction strategy.     

Pseudomonas syringae naturally lacking the canonical type III secretion system are ubiquitous in nonagricultural habitats,are phylogenetically diverse and can be pathogenic

The type III secretion system (T3SS) is an important virulence factor of pathogenic bacteria, but the natural occurrence of variants of bacterial plant pathogens with deficiencies in their T3SS raises questions about the significance of the T3SS for fitness. Previous work on T3SS-deficient plant pathogenic bacteria has focused on strains from plants or plant debris. Here we have characterized T3SS-deficient strains of Pseudomonas syringae from plant and nonplant substrates in pristine nonagricultural contexts, many of which represent recently described clades not yet found associated with crop plants. Strains incapable of inducing a hypersensitive reaction (HR⁻) in tobacco were detected in 65% of 126 samples from headwaters of rivers (mountain creeks and lakes), snowpack, epilithic biofilms, wild plants and leaf litter and constituted 2 to 100% of the P. syringae population associated with each sample. All HR⁻ strains lacked at least one gene in the canonical hrp/hrc locus or the associated conserved effector locus, but most lacked all six of the genes tested (hrcC, hrpL, hrpK1, avrE1 and hrpW1) and represented several disparate phylogenetic clades. Although most HR⁻ strains were incapable of causing symptoms on cantaloupe seedlings as expected, strains in the recently described TA-002 clade caused severe symptoms in spite of the absence of any of the six conserved genes of the canonical T3SS according to PCR and Southern blot assays. The phylogenetic context of the T3SS variants we observed provides insight into the evolutionary history of P. syringae as a pathogen and as an environmental saprophyte.     

Is It Easy to Be Urban? Convergent Success in Urban Habitats among Lineages of a Widespread Native Ant

The most rapidly expanding habitat globally is the urban habitat, yet the origin and life histories of the populations of native species that inhabit this habitat remain poorly understood. We use DNA barcoding of the COI gene in the widespread native pest ant Tapinoma sessile to test two hypotheses regarding the origin of urban populations and traits associated with their success. First, we determine if urban samples of T. sessile have a single origin from natural populations by looking at patterns of haplotype clustering from across their range. Second, we examine whether polygynous colony structure – a trait associated with invasion success – is correlated with urban environments, by studying the lineage dependence of colony structure. Our phylogenetic analysis of 49 samples identified four well supported geographic clades. Within clades, Kimura-2 parameter pairwise genetic distances revealed <2.3% variation; however, between clade genetic distances were 7.5–10.0%, suggesting the possibility of the presence of cryptic species. Our results indicate that T. sessile has successfully colonized urban environments multiple times. Additionally, polygynous colony structure is a highly plastic trait across habitat, clade, and haplotype. In short, T. sessile has colonized urban habitats repeatedly and appears to do so using life history strategies already present in more natural populations. Whether similar results hold for other species found in urban habitats has scarcely begun to be considered.     

Validation of Mixed-Genome Microarrays as a Method for Genetic Discrimination

Comparative genomic hybridizations have been used to examine genetic relationships among bacteria. The microarrays used in these experiments may have open reading frames from one or more reference strains (whole-genome microarrays), or they may be composed of random DNA fragments from a large number of strains (mixed-genome microarrays [MGMs]). In this work both experimental and virtual arrays are analyzed to assess the validity of genetic inferences from these experiments with a focus on MGMs. Empirical data are analyzed from an Enterococcus MGM, while a virtual MGM is constructed in silico using sequenced genomes (Streptococcus). On average, a small MGM is capable of correctly deriving phylogenetic relationships between seven species of Enterococcus with accuracies of 100% (n = 100 probes) and 95% (n = 46 probes); more probes are required for intraspecific differentiation. Compared to multilocus sequence methods and whole-genome microarrays, MGMs provide additional discrimination between closely related strains and offer the possibility of identifying unique strain or lineage markers. Representational bias can have mixed effects. Microarrays composed of probes from a single genome can be used to derive phylogenetic relationships, although branch length can be exaggerated for the reference strain. We describe a case where disproportional representation of different strains used to construct an MGM can result in inaccurate phylogenetic inferences, and we illustrate an algorithm that is capable of correcting this type of bias. The bias correction algorithm automatically provides bootstrap confidence values and can provide multiple bias-corrected trees with high confidence values.     

Genotypic characteristics of Haemophilus influenzae isolates from pediatric pneumonia patients in Chengdu city, Sichuan, China

Two hundred and seventy-three Haemophilus influenzae strains isolated from pediatric pneumonia patients in China were studied. We used Multilocus Sequence Typing (MLST) to analyze genotypic characteristics. All strains were biotyped and serotyped. Relatedness and patterns of genes among isolates were determined by the analysis of MLST and eBURST. H. influenzae primarily causes acute pneumonia in children under 1 year old. Nontypeable H. influenzae was responsible for most cases of pediatric pneumonia. All 273 strains were classified into eight biotypes. They mostly belonged to the I, II, and III biotypes (17.6%, 43.6%, and 22.7%, respectively). 62 strains (22.7%) produced b-lactamase. We found 28 novel alleles. Fifty different STs were found by MLST, of which 39 were novel. These were ST477 through ST508 and ST521 through ST527. Group 17 and predicted founders 503 were new groups in this study. No STs correlated with strains from Korea, which is adjacent to China. The H. influenzae strains from China appeared to have heterogeneous ST types patterns which may be the reason no outbreaks or epidemics of H. influenzae infections have occurred in Chengdu city, Sichuan, China.     

Identification of Genes and Genomic Islands Correlated with High Pathogenicity in Streptococcus suis Using Whole Genome Tilling Microarrays

Streptococcus suis is an important zoonotic pathogen that can cause meningitis and sepsis in both pigs and humans. Infections in humans have been sporadic worldwide but two severe outbreaks occurred in China in recent years, while infections in pigs are a major problem in the swine industry. Some S. suis strains are more pathogenic than others with 2 sequence types (ST), ST1 and ST7, being well recognized as highly pathogenic. We analyzed 31 isolates from 23 serotypes and 25 STs by NimbleGen tiling microarray using the genome of a high pathogenicity (HP) ST1 strain, GZ1, as reference and a new algorithm to detect gene content difference. The number of genes absent in a strain ranged from 49 to 225 with a total of 632 genes absent in at least one strain, while 1346 genes were found to be invariably present in all strains as the core genome of S. suis, accounting for 68% of the GZ1 genome. The majority of genes are located in chromosomal blocks with two or more contiguous genes. Sixty two blocks are absent in two or more strains and defined as regions of difference (RDs), among which 26 are putative genomic islands (GIs). Clustering and statistical analyses revealed that 8 RDs including 6 putative GIs and 21 genes within these RDs are significantly associated with HP. Three RDs encode known virulence related factors including the extracellular factor, the capsular polysaccharide and a SrtF pilus. The strains were divided into 5 groups based on population genetic analysis of multilocus sequence typing data and the distribution of the RDs among the groups revealed gain and loss of RDs in different groups. Our study elucidated the gene content diversity of S. suis and identified genes that potentially promote HP.     

Characterization of Non-Starter Lactic Acid Bacteria from Italian Ewe Cheeses Based on Phenotypic, Genotypic, and Cell Wall Protein Analyses

Non-starter lactic acid bacteria (NSLAB) were isolated from 12 Italian ewe cheeses representing six different types of cheese, which in several cases were produced by different manufacturers. A total of 400 presumptive Lactobacillus isolates were obtained, and 123 isolates and 10 type strains were subjected to phenotypic, genetic, and cell wall protein characterization analyses. Phenotypically, the cheese isolates included 32% Lactobacillus plantarum isolates, 15% L. brevis isolates, 12% L. paracasei subsp. paracasei isolates, 9% L. curvatus isolates, 6% L. fermentum isolates, 6% L. casei subsp. casei isolates, 5% L. pentosus isolates, 3% L. casei subsp. pseudoplantarum isolates, and 1% L. rhamnosus isolates. Eleven percent of the isolates were not phenotypically identified. Although a randomly amplified polymorphic DNA (RAPD) analysis based on three primers and clustering by the unweighted pair group method with arithmetic average (UPGMA) was useful for partially differentiating the 10 type strains, it did not provide a species-specific DNA band or a combination of bands which permitted complete separation of all the species considered. In contrast, sodium dodecyl sulfate-polyacrylamide gel electrophoresis cell wall protein profiles clustered by UPGMA were species specific and resolved the NSLAB. The only exceptions were isolates phenotypically identified as L. plantarum and L. pentosus or as L. casei subsp. casei and L. paracasei subsp. paracasei, which were grouped together. Based on protein profiles, Italian ewe cheeses frequently contained four different species and 3 to 16 strains. In general, the cheeses produced from raw ewe milk contained a larger number of more diverse strains than the cheeses produced from pasteurized milk. The same cheese produced in different factories contained different species, as well as strains that belonged to the same species but grouped in different RAPD clusters.