Benjamin Janto

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

Characterization and modeling of the Haemophilus influenzae core and supragenomes based on the complete genomic sequences of Rd and 12 clinical nontypeable strains

BackgroundThe 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.ResultsWe 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.ConclusionA 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.

Insights into the genome of large sulfur bacteria revealed by analysis of single filaments.

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.

Meta-omic characterization of the marine invertebrate microbial consortium that produces the chemotherapeutic natural product ET-743.

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.

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.

Identification of differentially expressed genes in fibroblasts derived from patients with Dupuytren's Contracture

Dupuytrens 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 Dupuytrens 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 Dupuytrens contracture and may ultimately serve as attractive molecular targets for alternative therapies.

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

Extensive Genomic Plasticity in Pseudomonas aeruginosa Revealed by Identification and Distribution Studies of Novel Genes among Clinical Isolates

ABSTRACT 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 hosts 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, ∼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.

Characterization, Distribution, and Expression of Novel Genes among Eight Clinical Isolates of Streptococcus pneumoniae

ABSTRACT 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 ∼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.

Cellulosilyticum ruminicola, a Newly Described Rumen Bacterium That Possesses Redundant Fibrolytic-Protein-Encoding Genes and Degrades Lignocellulose with Multiple Carbohydrate- Borne Fibrolytic Enzymes

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