Search

Abstract 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.

Abstract Dupuytren's contracture (DC) is the most common inherited connective tissue disease of humans and is hypothesized to be associated with aberrant wound healing of the palmar fascia. Fibroblasts and myofibroblasts are believed to play an important role in the genesis of DC and the fibroproliferation and contraction that are hallmarks of this disease. This study compares the gene expression profiles of fibroblasts isolated from DC patients and controls in an attempt to identify key genes whose regulation might be significantly altered in fibroblasts found within the palmar fascia of Dupuytren's patients. Total RNA isolated from diseased palmar fascia (DC) and normal palmar fascia (obtained during carpal tunnel release; 6 samples per group) was subjected to quantitative analyses using two different microarray platforms (GE Code Link™ and Illumina™) to identify and validate differentially expressed genes. The data obtained was analyzed using The Significance Analysis of Microarrays (SAM) software through which we identified 69 and 40 differentially regulated gene transcripts using the CodeLink™ and Illumina™ platforms, respectively. The CodeLink™ platform identified 18 upregulated and 51 downregulated genes. Using the Illumina™ platform, 40 genes were identified as downregulated, eleven of which were identified by both platforms. Quantitative RT-PCR confirmed the downregulation of three high-interest candidate genes which are all components of the extracellular matrix: proteoglycan 4 (PRG4), fibulin-1 (FBLN-1) transcript variant D, and type XV collagen alpha 1 chain. Overall, our study has identified a variety of candidate genes that may be involved in the pathophysiology of Dupuytren's contracture and may ultimately serve as attractive molecular targets for alternative therapies.

There is extensive genomic diversity among Streptococcus pneumoniae isolates. Approximately half of the comprehensive set of genes in the species (the supragenome or pangenome) is present in all the isolates (core set), and the remaining is unevenly distributed among strains (distributed set). The Streptococcus pneumoniae Supragenome Hybridization (SpSGH) array provides coverage for an extensive set of genes and polymorphisms encountered within this species, capturing this genomic diversity. Further, the capture is quantitative. In this manner, the SpSGH array allows for both genomic and transcriptomic analyses of diverse S. pneumoniae isolates on a single platform. In this unit, we present the SpSGH array, and describe in detail its design and implementation for both genomic and transcriptomic analyses. The methodology can be applied to construction and modification of SpSGH array platforms, as well to other bacterial species as long as multiple whole-genome sequences are available that collectively capture the vast majority of the species supragenome., (Copyright © 2015 John Wiley & Sons, Inc.)

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.

Two multidrug resistant strains of Streptococcus pneumoniae - SV35-T23 (capsular type 23F) and SV36-T3 (capsular type 3) were recovered from the nasopharynx of two adult patients during an outbreak of pneumococcal disease in a New York hospital in 1996. Both strains belonged to the pandemic lineage PMEN1 but they differed strikingly in virulence when tested in the mouse model of IP infection: as few as 1000 CFU of SV36 killed all mice within 24 hours after inoculation while SV35-T23 was avirulent.Whole genome sequencing (WGS) of the two isolates was performed (i) to test if these two isolates belonging to the same clonal type and recovered from an identical epidemiological scenario only differed in their capsular genes? and (ii) to test if the vast difference in virulence between the strains was mostly - or exclusively - due to the type III capsule. WGS demonstrated extensive differences between the two isolates including over 2500 single nucleotide polymorphisms in core genes and also differences in 36 genetic determinants: 25 of which were unique to SV35-T23 and 11 unique to strain SV36-T3. Nineteen of these differences were capsular genes and 9 bacteriocin genes.Using genetic transformation in the laboratory, the capsular region of SV35-T23 was replaced by the type 3 capsular genes from SV36-T3 to generate the recombinant SV35-T3* which was as virulent as the parental strain SV36-T3* in the murine model and the type 3 capsule was the major virulence factor in the chinchilla model as well. On the other hand, a careful comparison of strains SV36-T3 and the laboratory constructed SV35-T3* in the chinchilla model suggested that some additional determinants present in SV36 but not in the laboratory recombinant may also contribute to the progression of middle ear disease. The nature of this determinants remains to be identified.

Bacillus pseudofirmus OF4 is an extreme but facultative alkaliphile that grows non-fermentatively in a pH range from 7.5 to above 11.4 and can withstand large sudden increases in external pH. It is a model organism for studies of bioenergetics at high pH, at which energy demands are higher than at neutral pH because both cytoplasmic pH homeostasis and ATP synthesis require more energy. The alkaliphile also tolerates a cytoplasmic pH > 9.0 at external pH values at which the pH homeostasis capacity is exceeded, and manages other stresses that are exacerbated at alkaline pH, e.g. sodium, oxidative and cell wall stresses. The genome of B. pseudofirmus OF4 includes two plasmids that are lost from some mutants without viability loss. The plasmids may provide a reservoir of mobile elements that promote adaptive chromosomal rearrangements under particular environmental conditions. The genome also reveals a more acidic pI profile for proteins exposed on the outer surface than found in neutralophiles. A large array of transporters and regulatory genes are predicted to protect the alkaliphile from its overlapping stresses. In addition, unanticipated metabolic versatility was observed, which could ensure requisite energy for alkaliphily under diverse conditions., (© 2011 Society for Applied Microbiology and Blackwell Publishing Ltd.)

In many macroorganisms, the ultimate source of potent biologically active natural products has remained elusive due to an inability to identify and culture the producing symbiotic microorganisms. As a model system for developing a meta-omic approach to identify and characterize natural product pathways from invertebrate-derived microbial consortia, we chose to investigate the ET-743 (Yondelis) biosynthetic pathway. This molecule is an approved anticancer agent obtained in low abundance (10(-4)-10(-5) % w/w) from the tunicate Ecteinascidia turbinata and is generated in suitable quantities for clinical use by a lengthy semisynthetic process. On the basis of structural similarities to three bacterial secondary metabolites, we hypothesized that ET-743 is the product of a marine bacterial symbiont. Using metagenomic sequencing of total DNA from the tunicate/microbial consortium, we targeted and assembled a 35 kb contig containing 25 genes that comprise the core of the NRPS biosynthetic pathway for this valuable anticancer agent. Rigorous sequence analysis based on codon usage of two large unlinked contigs suggests that Candidatus Endoecteinascidia frumentensis produces the ET-743 metabolite. Subsequent metaproteomic analysis confirmed expression of three key biosynthetic proteins. Moreover, the predicted activity of an enzyme for assembly of the tetrahydroisoquinoline core of ET-743 was verified in vitro. This work provides a foundation for direct production of the drug and new analogues through metabolic engineering. We expect that the interdisciplinary approach described is applicable to diverse host-symbiont systems that generate valuable natural products for drug discovery and development.

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.

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.

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.

Cellulosilyticum ruminicola H1 is a newly described bacterium isolated from yak (Bos grunniens) rumen and is characterized by its ability to grow on a variety of hemicelluloses and degrade cellulosic materials. In this study, we performed the whole-genome sequencing of C. ruminicola H1 and observed a comprehensive set of genes encoding the enzymes essential for hydrolyzing plant cell wall. The corresponding enzymatic activities were also determined in strain H1; these included endoglucanases, cellobiohydrolases, xylanases, mannanase, pectinases, and feruloyl esterases and acetyl esterases to break the interbridge cross-link, as well as the enzymes that degrade the glycosidic bonds. This bacterium appears to produce polymer hydrolases that act on both soluble and crystal celluloses. Approximately half of the cellulytic activities, including cellobiohydrolase (50%), feruloyl esterase (45%), and one third of xylanase (31%) and endoglucanase (36%) activities were bound to cellulosic fibers. However, only a minority of mannase (6.78%) and pectinase (1.76%) activities were fiber associated. Strain H1 seems to degrade the plant-derived polysaccharides by producing individual fibrolytic enzymes, whereas the majority of polysaccharide hydrolases contain carbohydrate-binding module. Cellulosome or cellulosomelike protein complex was never isolated from this bacterium. Thus, the fibrolytic enzyme production of strain H1 may represent a different strategy in cellulase organization used by most of other ruminal microbes, but it applies the fungal mode of cellulose production.

Objectives: We sought to construct and partially characterize complementary DNA (cDNA) libraries prepared from the middle ear mucosa (MEM) of chinchillas to better understand pathogenic aspects of infection and inflammation, particularly with respect to leukotriene biogenesis and response., Methods: Chinchilla MEM was harvested from controls and after middle ear inoculation with nontypeable Haemophilus influenzae. RNA was extracted to generate cDNA libraries. Randomly selected clones were subjected to sequence analysis to characterize the libraries and to provide DNA sequence for phylogenetic analyses. Reverse transcription-polymerase chain reaction of the RNA pools was used to generate cDNA sequences corresponding to genes associated with leukotriene biosynthesis and metabolism., Results: Sequence analysis of 921 randomly selected clones from the uninfected MEM cDNA library produced approximately 250,000 nucleotides of almost entirely novel sequence data. Searches of the GenBank database with the Basic Local Alignment Search Tool provided for identification of 515 unique genes expressed in the MEM and not previously described in chinchillas. In almost all cases, the chinchilla cDNA sequences displayed much greater homology to human or other primate genes than with rodent species. Genes associated with leukotriene metabolism were present in both normal and infected MEM., Conclusions: Based on both phylogenetic comparisons and gene expression similarities with humans, chinchilla MEM appears to be an excellent model for the study of middle ear inflammation and infection. The higher degree of sequence similarity between chinchillas and humans compared to chinchillas and rodents was unexpected. The cDNA libraries from normal and infected chinchilla MEM will serve as useful molecular tools in the study of otitis media and should yield important information with respect to middle ear pathogenesis.

Dupuytren's contracture (DC) is the most common inherited connective tissue disease of humans and is hypothesized to be associated with aberrant wound healing of the palmar fascia. Fibroblasts and myofibroblasts are believed to play an important role in the genesis of DC and the fibroproliferation and contraction that are hallmarks of this disease. This study compares the gene expression profiles of fibroblasts isolated from DC patients and controls in an attempt to identify key genes whose regulation might be significantly altered in fibroblasts found within the palmar fascia of Dupuytren's patients. Total RNA isolated from diseased palmar fascia (DC) and normal palmar fascia (obtained during carpal tunnel release; 6 samples per group) was subjected to quantitative analyses using two different microarray platforms (GE Code Linktrade mark and Illuminatrade mark) to identify and validate differentially expressed genes. The data obtained was analyzed using The Significance Analysis of Microarrays (SAM) software through which we identified 69 and 40 differentially regulated gene transcripts using the CodeLinktrade mark and Illuminatrade mark platforms, respectively. The CodeLinktrade mark platform identified 18 upregulated and 51 downregulated genes. Using the Illuminatrade mark platform, 40 genes were identified as downregulated, eleven of which were identified by both platforms. Quantitative RT-PCR confirmed the downregulation of three high-interest candidate genes which are all components of the extracellular matrix: proteoglycan 4 (PRG4), fibulin-1 (FBLN-1) transcript variant D, and type XV collagen alpha 1 chain. Overall, our study has identified a variety of candidate genes that may be involved in the pathophysiology of Dupuytren's contracture and may ultimately serve as attractive molecular targets for alternative therapies.

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 ( approximately 3,000) that are represented in the S. pneumoniae population at frequencies of >or=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.

Marine sediments are frequently covered by mats of the filamentous Beggiatoa and other large nitrate-storing bacteria that oxidize hydrogen sulfide using either oxygen or nitrate, which they store in intracellular vacuoles. Despite their conspicuous metabolic properties and their biogeochemical importance, little is known about their genetic repertoire because of the lack of pure cultures. Here, we present a unique approach to access the genome of single filaments of Beggiatoa by combining whole genome amplification, pyrosequencing, and optical genome mapping. Sequence assemblies were incomplete and yielded average contig sizes of approximately 1 kb. Pathways for sulfur oxidation, nitrate and oxygen respiration, and CO2 fixation confirm the chemolithoautotrophic physiology of Beggiatoa. In addition, Beggiatoa potentially utilize inorganic sulfur compounds and dimethyl sulfoxide as electron acceptors. We propose a mechanism of vacuolar nitrate accumulation that is linked to proton translocation by vacuolar-type ATPases. Comparative genomics indicates substantial horizontal gene transfer of storage, metabolic, and gliding capabilities between Beggiatoa and cyanobacteria. These capabilities enable Beggiatoa to overcome non-overlapping availabilities of electron donors and acceptors while gliding between oxic and sulfidic zones. The first look into the genome of these filamentous sulfur-oxidizing bacteria substantially deepens the understanding of their evolution and their contribution to sulfur and nitrogen cycling in marine sediments.

Background: The distributed genome hypothesis (DGH) posits that chronic bacterial pathogens utilize polyclonal infection and reassortment of genic characters to ensure persistence in the face of adaptive host defenses. Studies based on random sequencing of multiple strain libraries suggested that free-living bacterial species possess a supragenome that is much larger than the genome of any single bacterium., Results: We derived high depth genomic coverage of nine nontypeable Haemophilus influenzae (NTHi) clinical isolates, bringing to 13 the number of sequenced NTHi genomes. Clustering identified 2,786 genes, of which 1,461 were common to all strains, with each of the remaining 1,328 found in a subset of strains; the number of clusters ranged from 1,686 to 1,878 per strain. Genic differences of between 96 and 585 were identified per strain pair. Comparisons of each of the NTHi strains with the Rd strain revealed between 107 and 158 insertions and 100 and 213 deletions per genome. The mean insertion and deletion sizes were 1,356 and 1,020 base-pairs, respectively, with mean maximum insertions and deletions of 26,977 and 37,299 base-pairs. This relatively large number of small rearrangements among strains is in keeping with what is known about the transformation mechanisms in this naturally competent pathogen., Conclusion: A finite supragenome model was developed to explain the distribution of genes among strains. The model predicts that the NTHi supragenome contains between 4,425 and 6,052 genes with most uncertainty regarding the number of rare genes, those that have a frequency of <0.1 among strains; collectively, these results support the DGH.

The distributed genome hypothesis (DGH) states that each strain within a bacterial species receives a unique distribution of genes from a population-based supragenome that is many times larger than the genome of any given strain. The observations that natural infecting populations are often polyclonal and that most chronic bacterial pathogens have highly developed mechanisms for horizontal gene transfer suggested the DGH and provided the means and the mechanisms to explain how chronic infections persist in the face of a mammalian host's adaptive defense mechanisms. Having previously established the validity of the DGH for obligate pathogens, we wished to evaluate its applicability to an opportunistic bacterial pathogen. This was accomplished by construction and analysis of a highly redundant pooled genomic library containing approximately 216,000 functional clones that was constructed from 12 low-passage clinical isolates of Pseudomonas aeruginosa, 6 otorrheic isolates and 6 from other body sites. Sequence analysis of 3,214 randomly picked clones (mean insert size, approximately 1.4 kb) from this library demonstrated that 348 (10.8%) of the clones were unique with respect to all genomic sequences of the P. aeruginosa prototype strain, PAO1. Hypothetical translations of the open reading frames within these unique sequences demonstrated protein homologies to a number of bacterial virulence factors and other proteins not previously identified in P. aeruginosa. PCR and reverse transcription-PCR-based assays were performed to analyze the distribution and expression patterns of a 70-open reading frame subset of these sequences among 11 of the clinical strains. These sequences were unevenly distributed among the clinical isolates, with nearly half (34/70) of the novel sequences being present in only one or two of the individual strains. Expression profiling revealed that a vast majority of these sequences are expressed, strongly suggesting they encode functional proteins.

Eight low-passage-number Streptococcus pneumoniae clinical isolates, each of a different serotype and a different multilocus sequence type, were obtained from pediatric participants in a pneumococcal vaccine trial. Comparative genomic analyses were performed with these strains and two S. pneumoniae reference strains. Individual genomic libraries were constructed for each of the eight clinical isolates, with an average insert size of approximately 1 kb. A total of 73,728 clones were picked for arraying, providing more than four times genomic coverage per strain. A subset of 4,793 clones were sequenced, for which homology searches revealed that 750 (15.6%) of the sequences were unique with respect to the TIGR4 reference genome and 263 (5.5%) clones were unrelated to any available streptococcal sequence. Hypothetical translations of the open reading frames identified within these novel sequences showed homologies to a variety of proteins, including bacterial virulence factors not previously identified in S. pneumoniae. The distribution and expression patterns of 58 of these novel sequences among the eight clinical isolates were analyzed by PCR- and reverse transcriptase PCR-based analyses, respectively. These unique sequences were nonuniformly distributed among the eight isolates, and transcription of these genes in planktonic cultures was detected in 81% (172/212) of their genic occurrences. All 58 novel sequences were transcribed in one or more of the clinical strains, suggesting that they all correspond to functional genes. Sixty-five percent (38/58) of these sequences were found in 50% or less of the clinical strains, indicating a significant degree of genomic plasticity among natural isolates.

We hypothesize that Haemophilus influenzae, as a species, possesses a much greater number of genes than that found in any single H. influenzae genome. This supragenome is distributed throughout naturally occurring infectious populations, and new strains arise through autocompetence and autotransformation systems. The effect is that H. influenzae populations can readily adapt to environmental stressors. The supragenome hypothesis predicts that significant differences exist between and among the genomes of individual infectious strains of nontypeable H. influenzae (NTHi). To test this prediction, we obtained 10 low-passage NTHi clinical isolates from the middle ear effusions of patients with chronic otitis media. DNA sequencing was performed with 771 clones chosen at random from a pooled genomic library. Homology searching demonstrated that approximately 10% of these clones were novel compared to the H. influenzae Rd KW20 genome, and most of them did not match any DNA sequence in GenBank. Amino acid homology searches using hypothetical translations of the open reading frames revealed homologies to a variety of proteins, including bacterial virulence factors not previously identified in the NTHi isolates. The distribution and expression of 53 of these genes among the 10 strains were determined by PCR- and reverse transcription PCR-based analyses. These unique genes were nonuniformly distributed among the 10 isolates, and transcription of these genes in planktonic cultures was detected in 50% (177 of 352) of the occurrences. All of the novel sequences were transcribed in one or more of the NTHi isolates. Seventeen percent (9 of 53) of the novel genes were identified in all 10 NTHi strains, with each of the remaining 44 being present in only a subset of the strains. These genic distribution analyses were more effective as a strain discrimination tool than either multilocus sequence typing or 23S ribosomal gene typing methods.