Genome comparison of bacterial pathogens.

Bacterial pathogens are being sequenced at an increasing rate. To many microbiologists, it appears that there simply is not enough time to digest all the information suddenly available. In this chapter we present several tools for comparison of sequenced pathogenic genomes, and discuss differences between pathogens and non-pathogens. The presented tools allow comparison of large numbers of genomes in a hypothesis-driven manner. Visualization of the results is very important for clear presentation of the results and various ways of graphical representation are introduced.

The genus burkholderia: analysis of 56 genomic sequences.

The genus Burkholderia consists of a number of very diverse species, both in terms of lifestyle (which varies from category B pathogens to apathogenic soil bacteria and plant colonizers) and their genetic contents. We have used 56 publicly available genomes to explore the genomic diversity within this genus, including genome sequences that are not completely finished, but are available from the NCBI database. Defining the pan- and core genomes of species results in insights in the conserved and variable fraction of genomes, and can verify (or question) historic, taxonomic groupings. We find only several hundred genes that are conserved across all Burkholderia genomes, whilst there are more than 40,000 gene families in the Burkholderia pan-genome. A BLAST matrix visualizes the fraction of conserved genes in pairwise comparisons. A BLAST atlas shows which genes are actually conserved in a number of genomes, located and visualized with reference to a chosen genome. Genomic islands are common in many Burkholderia genomes, and most of these can be readily visualized by DNA structural properties of the chromosome. Trees that are based on relatedness of gene family content yield different results depending on what genes are analyzed. Some of the differences can be explained by errors in incomplete genome sequences, but, as our data illustrate, the outcome of phylogenetic trees depends on the type of genes that are analyzed.

Distribution and characterization of staphylococcal interspersed repeat units (SIRUs) and potential use for strain differentiation.

Variable-number tandem repeats (VNTRs) have been shown to be a powerful tool in the determination of evolutionary relationships and population genetics of bacteria. The sequencing of a number of Staphylococcus aureus genomes has allowed the identification of novel VNTR sequences in S. aureus, which are similar to those used in the study of the evolution of Mycobacterium tuberculosis clades. Seven VNTRs, termed staphylococcal interspersed repeat units (SIRUs), distributed around the genome are described, occurring in both unique and multiple sites, and varying in length from 48 to 159 bp. Variations in copy numbers were observed in all loci, within both the sequenced genomes and the UK epidemic methicillin-resistant S. aureus (EMRSA) isolates. Clonally related UK EMRSA isolates were clustered using SIRUs, which provided a greater degree of discrimination than multi-locus sequence typing, indicating that VNTRs may be a more appropriate evolutionary marker for studying transmission events and the geographical spread of S. aureus clades.

Complete DNA sequence of the linear mitochondrial genome of the pathogenic yeast Candida parapsilosis.

The complete sequence of the mitochondrial DNA of the opportunistic yeast pathogen Candida parapsilosis was determined. The mitochondrial genome is represented by linear DNA molecules terminating with tandem repeats of a 738-bp unit. The number of repeats varies, thus generating a population of linear DNA molecules that are heterogeneous in size. The length of the shortest molecules is 30,922 bp, whereas the longer molecules have expanded terminal tandem arrays (nx738 bp). The mitochondrial genome is highly compact, with less than 8% of the sequence corresponding to non-coding intergenic spacers. In silico analysis predicted genes encoding fourteen protein subunits of complexes of the respiratory chain and ATP synthase, rRNAs of the large and small subunits of the mitochondrial ribosome, and twenty-four transfer RNAs. These genes are organized into two transcription units. In addition, six intronic ORFs coding for homologues of RNA maturase, reverse transcriptase and DNA endonucleases were identified. In contrast to its overall molecular architecture, the coding sequences of the linear mitochondrial DNA of C. parapsilosis are highly similar to their counterparts in the circular mitochondrial genome of its close relative C. albicans. The complete sequence has implications for both mitochondrial DNA replication and the evolution of linear DNA genomes.

Sequence analysis of three mitochondrial DNA molecules reveals interesting differences among Saccharomyces yeasts.

The complete sequences of mitochondrial DNA (mtDNA) from the two budding yeasts Saccharomyces castellii and Saccharomyces servazzii, consisting of 25 753 and 30 782 bp, respectively, were analysed and compared to Saccharomyces cerevisiae mtDNA. While some of the traits are very similar among Saccharomyces yeasts, others have highly diverged. The two mtDNAs are much more compact than that of S.cerevisiae and contain fewer introns and intergenic sequences, although they have almost the same coding potential. A few genes contain group I introns, but group II introns, otherwise found in S.cerevisiae mtDNA, are not present. Surprisingly, four genes (ATP6, COX2, COX3 and COB) in the mtDNA of S.servazzii contain, in total, five +1 frameshifts. mtDNAs of S.castellii, S.servazzii and S.cerevisiae contain all genes on the same strand, except for one tRNA gene. On the other hand, the gene order is very different. Several gene rearrangements have taken place upon separation of the Saccharomyces lineages, and even a part of the transcription units have not been preserved. It seems that the mechanism(s) involved in the generation of the rearrangements has had to ensure that all genes stayed encoded by the same DNA strand.

DNA microarray analysis of fim mutations in Escherichia coli.

Bacterial adhesion is often mediated by complex polymeric surface structures referred to as fimbriae. Type 1 fimbriae of Escherichia coli represent the archetypical and best characterised fimbrial system. These adhesive organelles mediate binding to D-mannose and are directly associated with virulence in the urinary tract. A typical type 1 fimbriated bacterium has up to 500 fimbriae on its surface, with each fimbria consisting of approximately 1000 individual subunits. This equates to approximately 8% of the total cellular protein and is potentially a significant resource drain for the cell. Here we have used DNA microarray analysis to examine the molecular events involved in response to fimbrial gene expression in E. coli K-12. Observed differential expression levels of the fim genes were in good agreement with our current knowledge of the stoichiometry of type 1 fimbriae. Changes in fim expression correlated directly with alterations in colony morphology. Deletion of the entire fim gene cluster resulted in the converse expression of another surface protein Antigen 43 (Ag43). Specific deletion of the fimH gene did not affect expression of other fim genes or Ag43, but did dramatically reduce the number of fimbriae expressed on the cell surface. The use of high-resolution oligonucleotide arrays for defining points of transcription initiation and termination is also demonstrated.

A DNA structural atlas for Escherichia coli.

We have performed a computational analysis of DNA structural features in 18 fully sequenced prokaryotic genomes using models for DNA curvature, DNA flexibility, and DNA stability. The structural values that are computed for the Escherichia coli chromosome are significantly different from (and generally more extreme than) that expected from the nucleotide composition. To aid this analysis, we have constructed tools that plot structural measures for all positions in a long DNA sequence (e.g. an entire chromosome) in the form of color-coded wheels (http://www.cbs.dtu. dk/services/GenomeAtlas/). We find that these "structural atlases" are useful for the discovery of interesting features that may then be investigated in more depth using statistical methods. From investigation of the E. coli structural atlas, we discovered a genome-wide trend, where an extended region encompassing the terminus displays a high of level curvature, a low level of flexibility, and a low degree of helix stability. The same situation is found in the distantly related Gram-positive bacterium Bacillus subtilis, suggesting that the phenomenon is biologically relevant. Based on a search for long DNA segments where all the independent structural measures agree, we have found a set of 20 regions with identical and very extreme structural properties. Due to their strong inherent curvature, we suggest that these may function as topological domain boundaries by efficiently organizing plectonemically supercoiled DNA. Interestingly, we find that in practically all the investigated eubacterial and archaeal genomes, there is a trend for promoter DNA being more curved, less flexible, and less stable than DNA in coding regions and in intergenic DNA without promoters. This trend is present regardless of the absolute levels of the structural parameters, and we suggest that this may be related to the requirement for helix unwinding during initiation of transcription, or perhaps to the previously observed location of promoters at the apex of plectonemically supercoiled DNA. We have also analyzed the structural similarities between groups of genes by clustering all RNA and protein-encoding genes in E. coli, based on the average structural parameters. We find that most ribosomal genes (protein-encoding as well as rRNA genes) cluster together, and we suggest that DNA structure may play a role in the transcription of these highly expressed genes.

A classification of possible routes of Darwinian evolution.

A classification of four possible routes of Darwinian evolution is presented. These are serial direct evolution, parallel direct evolution, elimination of functional redundancy, and adoption from a different function. This classification provides a conceptual framework within which to investigate the accessibility by Darwinian evolution of complex biological structures.

Structural analysis of DNA sequence: evidence for lateral gene transfer in Thermotoga maritima.

The recently published complete DNA sequence of the bacterium Thermotoga maritima provides evidence, based on protein sequence conservation, for lateral gene transfer between Archaea and Bacteria. We introduce a new method of periodicity analysis of DNA sequences, based on structural parameters, which brings independent evidence for the lateral gene transfer in the genome of T.maritima. The structural analysis relates the Archaea-like DNA sequences to the genome of Pyrococcus horikoshii. Analysis of 24 complete genomic DNA sequences shows different periodicity patterns for organisms of different origin. The typical genomic periodicity for Bacteria is 11 bp whilst it is 10 bp for Archaea. Eukaryotes have more complex spectra but the dominant period in the yeast Saccharomyces cerevisiae is 10.2 bp. These periodicities are most likely reflective of differences in chromatin structure.

Visualization of pathogenicity regions in bacteria.

We show here how pathogenicity islands can be analysed using GenomeAtlases, which is a method for visualising repeats, DNA structural characteristics, and base composition of chromosomes and plasmids. We have applied this method to the E. coli plasmid pO157, and the Y. pestis plasmid pPCP1. In both cases pathogenic genes were shown to differ in A + T content and structural properties. Furthermore, examination of an antibiotic resistance gene cluster from S. typhimurium showed that the same was true for genes encoding antibiotic resistance.

Environmental influences on DNA curvature.

DNA curvature plays an important role in many biological processes. To study environmental influences on DNA curvature we compared the anomalous migration on polyacrylamide gels of ligation ladders of 11 specifically-designed oligonucleotides. At low temperatures (25 degrees C and below) most of the sequences exhibited a degree of anomalous migration. Increased temperature had a significant effect on the anomalous migration (curvature) of some sequences but limited effects on others; at 50 degrees C only 1 sequence migrated anomalously. Mg2+ had a strong influence on the migration of certain sequences, whilst spermine enhanced the anomalous migration of a different set of sequences. Sequences with a GGC motif exhibited greater curvature than predicted by the presently-used angles for the nearest-neighbour wedge model and are especially sensitive to Mg2+. The data have implications for models for DNA curvature and for environmentally-sensitive DNA conformations in the regulation of gene expression.

Three views of microbial genomes.

We describe here GenomeAtlases as a method for visualising three different aspects of complete microbial chromosomes: repeats, DNA structural characteristics, and base composition. We have applied this method to all publicly available genomes, and find a general strand preference of global repeats. The atlas for the Mycoplasma genitalium genome is presented as an example, and results from all three views are consistent with known characteristics of the genome.

DNA binding is not sufficient for H-NS-mediated repression of proU expression.

H-NS is a major component of bacterial chromatin and influences the expression of many genes. H-NS has been shown to exhibit a binding preference for certain AT-rich curved DNA elements in vitro. In this study we have addressed the factors that determine the specificity of H-NS action in vitro and in vivo. In bandshift studies, H-NS showed a slight binding preference for all curved sequences tested whether GC-based or AT-based; the specific architecture of the curve also influenced H-NS binding. In filter retention assays little difference in affinity could be detected for any sequence tested, including the downstream regulatory element (DRE) a downstream curved DNA element required for H-NS to repress transcription of the Salmonella typhimurium proU operon in vivo. A Kd of 1-2 microM was estimated for binding of H-NS to each of these sequences. In vivo, the distance between the proU promoter and the DRE, their relative orientations on the face of the DNA helix, and translation of the DRE had no major effect on proU regulation. None of the synthetic curved sequences tested could functionally replace the DRE in vivo. These data show that differential binding to curved DNA cannot account for the specificity of H-NS action in vivo. Furthermore, binding of H-NS to DNA per se is insufficient to repress the proU promoter. Thus, the DRE does not simply act as an H-NS binding site but must have a more specific role in mediating H-NS regulation of proU transcription.

Formation of a combined H-DNA/open TATA box structure in the promoter sequence of the human Na,K-ATPase alpha2 gene.

Structural variation of DNA within the promoter of the human Na, K-ATPase alpha2 gene, which contains a 35-base pair (bp) homopyrimidine.homopurine (Py.Pu) tract adjacent to a TATA box has been studied. The Py.Pu tract contains a 26-bp quasi-mirror repeat sequence with a potential for intramolecular triplex formation. As analyzed by two-dimensional agarose gel electrophoresis, a plasmid containing 151 bp of the promoter sequence including the 35-bp Py.Pu tract undergoes structural transitions under moderately acidic pH. Chemical probing with chloroacetaldehyde, dimethyl sulfate, and potassium permanganate is consistent with the formation of triplex DNA within the Py.Pu tract at native superhelical density as isolated from Escherichia coli. Chemical probing was used to determine a supercoil dependence for the formation of this combined unwound structure. At the superhelical density sufficient to locally unwind DNA, an H-y3 isomer of intermolecular triplex likely forms. However, at higher superhelical tension an H-y5 structure forms in the Py.Pu tract, and with increasing supercoiling the local DNA unwinding extends into the abutting TATA box. The H-y5/open TATA box combination structure might be favorable at higher superhelical densities since it relaxes more supercoils. The possible involvement of the H-y5/open TATA box structure in transcription is discussed.

The chromatin-associated protein H-NS.

H-NS is a major component of chromatin in enteric bacteria. H-NS plays a structural role in organising the chromosome, and influences DNA rearrangements as well as the expression of many genes. The biochemical and functional characteristics of H-NS are distinct from those of 'typical' DNA-binding proteins and much remains to be learned about the mechanism(s) by which H-NS acts. In this article we review our current understanding of the role of H-NS, and describe possible models by which H-NS might influence DNA structure and gene expression.

The chromatin-associated protein H-NS alters DNA topology in vitro.

H-NS is one of the two most abundant proteins in the bacterial nucleoid and influences the expression of a number of genes. We have studied the interaction of H-NS with DNA; purified H-NS was demonstrated to constrain negative DNA supercoils in vitro. This provides support for the hypothesis that H-NS influences transcription via changes in DNA topology, and is evidence of a structural role for H-NS in bacterial chromatin. The effects of H-NS on topology were only observed at sub-saturating concentrations of the protein. In addition, a preferred binding site on DNA was identified by DNase I footprinting at sub-saturating H-NS concentrations. This site corresponded to a curved sequence element which we previously showed, by in vivo studies, to be a site at which H-NS influences transcription of the proU operon. When present in saturating concentrations, H-NS did not constrain supercoils and bound to DNA in a sequence-independent fashion, covering all DNA molecules from end to end, suggesting that H-NS may form distinct complexes with DNA at different H-NS:DNA ratios. The data presented here provide direct support for the hypothesis that H-NS acts at specific sites to influence DNA topology and, hence, transcription.

Environmental influences on the in vivo level of intramolecular triplex DNA in Escherichia coli.

We describe an assay for detecting intramolecular triple-stranded DNA in living cells. The assay involves quantitative analysis of the differential photobinding of 4,5',8-trimethylpsoralen to 5'TA and 5'AT dinucleotides in a region of plasmid DNA containing the triplex-forming sequence (GA)7TA(GA)7. Psoralen photobinds to the central 5'TA within the (GA)7TA(GA)7 sequence in duplex DNA but not when the sequence exists as an intramolecular triplex structure. The reactivity of photobinding sites in regions flanking the intramolecular triplex-forming sequence, those that comprise the duplex-triplex junctions, either increase or decrease upon formation of the intramolecular triplex structure from duplex DNA. The pattern of trimethylpsoralen reactivity provides an indication of the conformation of the (GA)7TA(GA)7 sequence in plasmid DNA. The formation of the intramolecular triplex structure is dependent on both pH and negative superhelical density in vitro. The fraction of the (GA)7TA(GA)7 sequence that existed as an intramolecular triplex structure in vivo was dependent on the level of DNA supercoiling in vivo and changes in growth conditions which influence the intracellular pH. The Hy3 isomer was detected in living Escherichia coli cells.