Mycobacterial species as case-study of comparative genome analysis.

The genus Mycobacterium represents more than 120 species including important pathogens of human and cause major public health problems and illnesses. Further, with more than 100 genome sequences from this genus, comparative genome analysis can provide new insights for better understanding the evolutionary events of these species and improving drugs, vaccines, and diagnostics tools for controlling Mycobacterial diseases. In this present study we aim to outline a comparative genome analysis of fourteen Mycobacterial genomes: M. avium subsp. paratuberculosis K­10, M. bovis AF2122/97, M. bovis BCG str. Pasteur 1173P2, M. leprae Br4923, M. marinum M, M. sp. KMS, M. sp. MCS, M. tuberculosis CDC1551, M. tuberculosis F11, M. tuberculosis H37Ra, M. tuberculosis H37Rv, M. tuberculosis KZN 1435 , M. ulcerans Agy99,and M. vanbaalenii PYR­1, For this purpose a comparison has been done based on their length of genomes, GC content, number of genes in different data bases (Genbank, Refseq, and Prodigal). The BLAST matrix of these genomes has been figured to give a lot of information about the similarity between species in a simple scheme. As a result of multiple genome analysis, the pan and core genome have been defined for twelve Mycobacterial species. We have also introduced the genome atlas of the reference strain M. tuberculosis H37Rv which can give a good overview of this genome. And for examining the phylogenetic relationships among these bacteria, a phylogenic tree has been constructed from 16S rRNA gene for tuberculosis and non tuberculosis Mycobacteria to understand the evolutionary events of these species.

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.

On the total number of genes and their length distribution in complete microbial genomes.

In sequenced microbial genomes, some of the annotated genes are actually not protein-coding genes, but rather open reading frames that occur by chance. Therefore, the number of annotated genes is higher than the actual number of genes for most of these microbes. Comparison of the length distribution of the annotated genes with the length distribution of those matching a known protein reveals that too many short genes are annotated in many genomes. Here we estimate the true number of protein-coding genes for sequenced genomes. Although it is often claimed that Escherichia coli has about 4300 genes, we show that it probably has only approximately 3800 genes, and that a similar discrepancy exists for almost all published genomes.

Genome organisation and chromatin structure in Escherichia coli.

We have analysed the complete sequence of the Escherichia coli K12 isolate MG1655 genome for chromatin-associated protein binding sites, and compared the predicted location of predicted sites with experimental expression data from 'DNA chip' experiments. Of the dozen proteins associated with chromatin in E. coli, only three have been shown to have significant binding preferences: integration host factor (IHF) has the strongest binding site preference, and FIS sites show a weak consensus, and there is no clear consensus site for binding of the H-NS protein. Using hidden Markov models (HMMs), we predict the location of 608 IHF sites, scattered throughout the genome. A subset of the IHF sites associated with repeats tends to be clustered around the origin of replication. We estimate there could be roughly 6000 FIS sites in E. coli, and the sites tend to be localised in two regions flanking the replication termini. We also show that the regions upstream of genes regulated by H-NS are more curved and have a higher AT content than regions upstream of other genes. These regions in general would also be localised near the replication terminus.

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.

Critical role of Lyn kinase in inhibition of neutrophil apoptosis by granulocyte-macrophage colony-stimulating factor.

The signal pathway for control of apoptosis in human neutrophils is currently unknown. In this study, we provide the first evidence that a Src family tyrosine kinase, Lyn, plays a key role in inhibition polymorphonuclear (PMN) cell death. Several nuclear proteins associated with apoptosis, i.e., p53, cdc2, and Rb, were absent from PMN. Bcl-2, known to inhibit apoptosis, was also not expressed. Programmed cell death that rapidly occurred in PMN could be arrested by granulocyte-macrophage CSF (GM-CSF), but this activation did not induce p53, cdc2, retinoblastoma, or Bcl-2 expression. Instead, GM-CSF produced a rapid activation of Lyn and Hck, but not Fgr, tyrosine phosphorylation within 1 min. Co-immunoprecipitation studies indicated that only Lyn, but not Hck, was physically coupled to GM-CSF receptor. By histologic assessment and evaluation of DNA fragmentation, only antisense Lyn, but not antisense Hck or antisense Fgr, could reverse the cell survival advantage provided by GM-CSF. Therefore, the physical coupling of Lyn to GM-CSF receptor and its early activation are required for inhibition or delay of apoptosis in PMN.

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.

Fiber-coupler-based measurement of lenslet focal lengths.

In a technique reminiscent of the classic autocollimation method, we demonstrate that focal lengths of micro-optic lenses may be measured precisely by the use of a mirror and a 2 × 1 fiber coupler.

Interleukin 15 induction of lymphokine-activated killer cell function against autologous tumor cells in melanoma patient lymphocytes by a CD18-dependent, perforin-related mechanism.

Interleukin 15 (IL-15) is a novel cytokine that shares no homology with IL-2, but it requires the use of beta and gamma chains of the IL-2 receptor complex for binding and signaling. In vitro studies have shown induction of CTL and lymphokine-activated killer (LAK) cell activity in peripheral blood mononuclear cells (PBMCs) from normal donors by IL-15 against known tumor targets. The present study attempts to define the role of IL-15 in generating LAK activity from melanoma patient lymphocytes. PBMCs of patients newly diagnosed with metastatic melanoma were incubated with different doses of recombinant human IL-15 and tested against autologous tumor cells, LAK sensitive cell lines (i.e., FMEX and Daudi), as well as the natural killer-sensitive cell line K562, in a 15-h 51Cr release assay. The effect of IL-15 was found to be both time and dose dependent, with peak activity detected after 2 or 3 days of culture with 100 ng/ml of this cytokine. LAK and not CTL activity in patient PBMCs was detected by the inability of mAbs against CD4, CD8, and MHC class I to effectively block lysis of autologous tumor and FMEX melanoma cells. In addition, interaction via the CD18 adhesion molecule was shown to be critical in IL-15-induced LAK-mediated lysis of autologous tumor cells. Finally, incubation of patient PBMCs with IL-15 for 6 h resulted in the up-regulation of perforin mRNA transcription. These findings suggest that LAK activity can be generated from melanoma patient PBMCs in the presence of IL-15 to lyse autologous tumor cells in a non-MHC-restricted manner. This new cytokine may play an important role in antitumor immunity with a possible use for cancer immunotherapy.