Multiplex sequencing of 1.5 Mb of the Mycobacterium leprae genome

Inteins as Drug Targets and Therapeutic Tools

Multidrug-resistant pathogens are of significant concern in recent years. Hence new antifungal and anti-bacterial drug targets are urgently needed before the situation goes beyond control. Inteins are polypeptides that self-splice from exteins without the need for cofactors or external energy, resulting in joining of extein fragments. Inteins are present in many organisms, including human pathogens such as Mycobacterium tuberculosis, Cryptococcus neoformans, C. gattii, and Aspergillus fumigatus. Because intein elements are not present in human genes, they are attractive drug targets to develop antifungals and antibiotics. Thus far, a few inhibitors of intein splicing have been reported. Metal-ions such as Zn²⁺ and Cu²⁺, and platinum-containing compound cisplatin inhibit intein splicing in M. tuberculosis and C. neoformans by binding to the active site cysteines. A small-molecule inhibitor 6G-318S and its derivative 6G-319S are found to inhibit intein splicing in C. neoformans and C. gattii with a MIC in nanomolar concentrations. Inteins have also been used in many other applications. Intein can be used in activating a protein inside a cell using small molecules. Moreover, split intein can be used to deliver large genes in experimental gene therapy and to kill selected species in a mixed population of microbes by taking advantage of the toxin-antitoxin system. Furthermore, split inteins are used in synthesizing cyclic peptides and in developing cell culture model to study infectious viruses including SARS-CoV-2 in the biosafety level (BSL) 2 facility. This mini-review discusses the recent research developments of inteins in drug discovery and therapeutic research.

Qualitative and quantitative proteomic analysis of Vitamin C induced changes in Mycobacterium smegmatis

Vitamin C is a critical dietary nutrient in human which has a wide range of regulatory effects on gene expression and physiology of Mycobacterium tuberculosis that leads to a dormant drug-tolerant phenotype. In the presence of iron, vitamin C shows a high bactericidal activity even in the drug resistant phenotype of M. tuberculosis. The regulatory mechanisms underlying vitamin C induced adaptations are largely unknown due to lack of functional genomics data in this field. In this study, we attempt to characterize the direct effect of vitamin C treatment on the physiology of actively growing Mycobacterium smegmatis. The study chose M. smegmatis as it is a fast-growing bacterium and a non-pathogenic model system which shares many physiological features with the pathogenic M. tuberculosis including dormancy and its regulation. The proteomic adaptation of M. smegmatis on vitamin C treatment demonstrates the important changes in cellular and metabolic process such as reversal of tricarboxylic acid cycle, decrease in ATP synthesis, decrease in iron acquisition and storage, and induction of dormancy regulators WhiB3, PhoP, and Lsr2.

Identification of functional candidates amongst hypothetical proteins of Mycobacterium leprae Br4923, a causative agent of leprosy

Mycobacterium leprae is an intracellular obligate parasite that causes leprosy in humans, and it leads to the destruction of peripheral nerves and skin deformation. Here, we report an extensive analysis of the hypothetical proteins (HPs) from M. leprae strain Br4923, assigning their functions to better understand the mechanism of pathogenesis and to search for potential therapeutic interventions. The genome of M. leprae encodes 1604 proteins, of which the functions of 632 are not known (HPs). In this paper, we predicted the probable functions of 312 HPs. First, we classified all HPs into families and subfamilies on the basis of sequence similarity, followed by domain assignment, which provides many clues for their possible function. However, the functions of 320 proteins were not predicted because of low sequence similarity with proteins of known function. Annotated HPs were categorized into enzymes, binding proteins, transporters, and proteins involved in cellular processes. We found several novel proteins whose functions were unknown for M. leprae. These proteins have a requisite association with bacterial virulence and pathogenicity. Finally, our sequence-based analysis will be helpful for further validation and the search for potential drug targets while developing effective drugs to cure leprosy.

Activity, specificity and structure of I-Bth0305I: a representative of a new homing endonuclease family

Novel family of putative homing endonuclease genes was recently discovered during analyses of metagenomic and genomic sequence data. One such protein is encoded within a group I intron that resides in the recA gene of the Bacillus thuringiensis 0305ϕ8–36 bacteriophage. Named I-Bth0305I, the endonuclease cleaves a DNA target in the uninterrupted recA gene at a position immediately adjacent to the intron insertion site. The enzyme displays a multidomain, homodimeric architecture and footprints a DNA region of ∼60 bp. Its highest specificity corresponds to a 14-bp pseudopalindromic sequence that is directly centered across the DNA cleavage site. Unlike many homing endonucleases, the specificity profile of the enzyme is evenly distributed across much of its target site, such that few single base pair substitutions cause a significant decrease in cleavage activity. A crystal structure of its C-terminal domain confirms a nuclease fold that is homologous to very short patch repair (Vsr) endonucleases. The domain architecture and DNA recognition profile displayed by I-Bth0305I, which is the prototype of a homing lineage that we term the ‘EDxHD’ family, are distinct from previously characterized homing endonucleases.

Mutational analysis of active-site residues in the Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease: Asp(122) and Asp(193) are crucial to the double-stranded DNA cleavage activity whereas Asp(218) is not

Mycobacterium leprae recA harbors an in-frame insertion sequence that encodes an intein homing endonuclease (PI-MleI). Most inteins (intein endonucleases) possess two conserved LAGLIDADG (DOD) motifs at their active center. A common feature of LAGLIDADG-type homing endonucleases is that they recognize and cleave the same or very similar DNA sequences. However, PI-MleI is distinctive from other members of the family of LAGLIDADG-type HEases for its modular structure with functionally separable domains for DNA-binding and cleavage, each with distinct sequence preferences. Sequence alignment analyses of PI-MleI revealed three putative LAGLIDADG motifs; however, there is conflicting bioinformatics data in regard to their identity and specific location within the intein polypeptide. To resolve this conflict and to determine the active-site residues essential for DNA target site recognition and double-stranded DNA cleavage, we performed site-directed mutagenesis of presumptive catalytic residues in the LAGLIDADG motifs. Analysis of target DNA recognition and kinetic parameters of the wild-type PI-MleI and its variants disclosed that the two amino acid residues, Asp(122) (in Block C) and Asp(193) (in functional Block E), are crucial to the double-stranded DNA endonuclease activity, whereas Asp(218) (in pseudo-Block E) is not. However, despite the reduced catalytic activity, the PI-MleI variants, like the wild-type PI-MleI, generated a footprint of the same length around the insertion site. The D122T variant showed significantly reduced catalytic activity, and D122A and D193A mutations although failed to affect their DNA-binding affinities, but abolished the double-stranded DNA cleavage activity. On the other hand, D122C variant showed approximately twofold higher double-stranded DNA cleavage activity, compared with the wild-type PI-MleI. These results provide compelling evidence that Asp(122) and Asp(193) in DOD motif I and II, respectively, are bona fide active-site residues essential for DNA cleavage activity. The implications of these results are discussed in this report.

Important role of the nucleotide excision repair pathway in Mycobacterium smegmatis in conferring protection against commonly encountered DNA-damaging agents

Mycobacteria are an important group of human pathogens. Although the DNA repair mechanisms in mycobacteria are not well understood, these are vital for the pathogen's persistence in the host macrophages. In this study, we generated a null mutation in the uvrB gene of Mycobacterium smegmatis to allow us to compare the significance of the nucleotide excision repair (NER) pathway with two important base excision repair pathways, initiated by uracil DNA glycosylase (Ung) and formamidopyrimidine DNA glycosylase (Fpg or MutM), in an isogenic strain background. The strain deficient in NER was the most sensitive to commonly encountered DNA-damaging agents such as UV, low pH, reactive oxygen species, hypoxia, and was also sensitive to acidified nitrite. Taken together with previous observations on NER-deficient M. tuberculosis, these results suggest that NER is an important DNA repair pathway in mycobacteria.

Recent developments in mycobacterial research

Tuberculosis remains a major health problem in the world, which is compounded further by the alarmingly high rate of M. tuberculosis infections in AIDS patients. Thus, there is an urgent need to advance our understanding of the mycobacterium to develop new drugs. The extraordinary recent developments in mycobacterial genetic research, particularly in genomics will greatly facilitate this goal. The knowledge of the entire genome sequence of M. tuberculosis will help in designing new chemotherapeutic and immunotherapeutic interventions. This review highlights recent developments in genomics, mycobacterial genetics, novel vaccine strategies, and our understanding of tuberculous dormancy.

Conserved indels in protein sequences that are characteristic of the phylum Actinobacteria

Gram-positive bacteria with a high G+C content are currently recognized as a distinct phylum, Actinobacteria, on the basis of their branching in 16S rRNA trees. Except for an insert in the 23S rRNA, there are no unique biochemical or molecular characteristics known at present that can distinguish this group from all other bacteria. In this work, three conserved indels (i.e. inserts or deletions) are described in three widely distributed proteins that are distinctive characteristics of the Actinobacteria and are not found in any other groups of bacteria. The identified signatures are a 2 aa deletion in cytochrome-c oxidase subunit 1 (Cox1), a 4 aa insert in CTP synthetase and a 5 aa insert in glutamyl-tRNA synthetase (GluRS). Additionally, the actinobacterial specificity of the large insert in the 23S rRNA was also tested. Using primers designed for conserved regions flanking these signatures, fragments of most of these genes were amplified from 23 actinobacterial species, covering many different families and orders, for which no sequence information was previously available. All the 61 sequenced fragments, except two in GluRS, were found to contain the indicated signatures. The presence of these signatures in various species from 20 families within this phylum provides evidence that they are likely distinctive characteristics of the entire phylum, which were introduced in a common ancestor of this group. The absence of all four of these signatures in Symbiobacterium thermophilum suggests that this species, which is distantly related to other actinobacteria in 16S rRNA and CTP synthetase trees, may not be a part of the phylum Actinobacteria. The identified signatures provide novel molecular means for defining and circumscribing the phylum Actinobacteria. Functional studies on them should prove helpful in understanding novel biochemical and physiological characteristics of this group of bacteria.

Characterization of theStreptomyces coelicolor A3(2)wblE gene, encoding a homologue of the sporulation transcription factor

The gene corresponding to the recently identified whiB-paralogous gene wblE in S. coelicolor was found after sequencing the downstream region of the stress-response sporulation-specific sigma-factor gene, sigH, in S. coelicolor A3(2). Sequence analysis has revealed an ORF exhibiting high similarity to sporulation transcription factors WhiB and WhiD. A stable null mutant of the wblE gene was obtained by integrative transformation, via double cross-over. Disruption of the S. coelicolor wblE gene appeared to have no obvious effect on growth, morphology, differentiation, and production of the pigmented antibiotics actinorhodin and undecylprodigiosin. Expression of the wblE gene was investigated during differentiation by S1 nuclease mapping, using RNA prepared from S. coelicolor A3(2) and its isogenic sigF and sigH mutants grown to various developmental stages. A single promoter was identified upstream of the wblE coding region. The wblEp promoter was induced at the beginning of aerial mycelium formation and its activity decreased later in differentiation. No differences in expression of the wblEp promoter were detected in S. coelicolor A3(2) mutants in sigF and sigH genes for sporulation-specific sigma factors. Sequence of the wblEp promoter showed partial similarity to the consensus sequence of the extracytoplasmic sigma factors.

Generalized gap model for bacterial artificial chromosome clone fingerprint mapping and shotgun sequencing

We develop an extension to the Lander-Waterman theory for characterizing gaps in bacterial artificial chromosome fingerprint mapping and shotgun sequencing projects. It supports a larger set of descriptive statistics and is applicable to a wider range of project parameters. We show that previous assertions regarding inconsistency of the Lander-Waterman theory at higher coverages are incorrect and that another well-known but ostensibly different model is in fact the same. The apparent paradox of infinite island lengths is resolved. Several applications are shown, including evolution of the probability density function, calculation of closure probabilities, and development of a probabilistic method for computing stopping points in bacterial artificial chromosome shotgun sequencing.

Strategies for the systematic sequencing of complex genomes

Recent spectacular advances in the technologies and strategies for DNA sequencing have profoundly accelerated the detailed analysis of genomes from myriad organisms. The past few years alone have seen the publication of near-complete or draft versions of the genome sequence of several well-studied, multicellular organisms - most notably, the human. As well as providing data of fundamental biological significance, these landmark accomplishments have yielded important strategic insights that are guiding current and future genome-sequencing projects.

Introduction of a carboxyl group in the first transmembrane helix of Escherichia coli F1Fo ATPase subunit c and cytoplasmic pH regulation

The multicopy subunit c of the H(+)-transporting F1Fo ATP synthase of Escherichia coli folds across the membrane as a hairpin of two hydrophobic alpha helices. The subunits interact in a front-to-back fashion, forming an oligomeric ring with helix 1 packing in the interior and helix 2 at the periphery. A conserved carboxyl, Asp(61) in E. coli, centered in the second transmembrane helix is essential for H+ transport. A second carboxylic acid in the first transmembrane helix is found at a position equivalent to Ile28 in several bacteria, some the cause of serious infectious disease. This side chain has been predicted to pack proximal to the essential carboxyl in helix 2. It appears that in some of these bacteria the primary function of the enzyme is H+ pumping for cytoplasmic pH regulation. In this study, Ile28 was changed to Asp and Glu. Both mutants were functional. However, unlike the wild type, the mutants showed pH-dependent ATPase-coupled H+ pumping and passive H+ transport through Fo. The results indicate that the presence of a second carboxylate enables regulation of enzyme function in response to cytoplasmic pH and that the ion binding pocket is aqueous accessible. The presence of a single carboxyl at position 28, in mutants I28D/D61G and I28E/D61G, did not support growth on a succinate carbon source. However, I28E/D61G was functional in ATPase-coupled H+ transport. This result indicates that the side chain at position 28 is part of the ion binding pocket.

Massive gene decay in the leprosy bacillus

Leprosy, a chronic human neurological disease, results from infection with the obligate intracellular pathogen Mycobacterium leprae, a close relative of the tubercle bacillus. Mycobacterium leprae has the longest doubling time of all known bacteria and has thwarted every effort at culture in the laboratory. Comparing the 3.27-megabase (Mb) genome sequence of an armadillo-derived Indian isolate of the leprosy bacillus with that of Mycobacterium tuberculosis (4.41 Mb) provides clear explanations for these properties and reveals an extreme case of reductive evolution. Less than half of the genome contains functional genes but pseudogenes, with intact counterparts in M. tuberculosis, abound. Genome downsizing and the current mosaic arrangement appear to have resulted from extensive recombination events between dispersed repetitive sequences. Gene deletion and decay have eliminated many important metabolic activities including siderophore production, part of the oxidative and most of the microaerophilic and anaerobic respiratory chains, and numerous catabolic systems and their regulatory circuits.

From genes to sequences to antibiotics: prospects for future developments from microbial genomics

The entire genome sequences for a number of microbial organisms have become available over the past few years. This knowledge is the beginning point for understanding the fundamental principles of bacterial structure and function. The prospects for gain in knowledge from genomics are discussed in this report.

A bacterial genome in flux: the twelve linear and nine circular extrachromosomal DNAs in an infectious isolate of the Lyme disease spirochete Borrelia burgdorferi

We have determined that Borrelia burgdorferi strain B31 MI carries 21 extrachromosomal DNA elements, the largest number known for any bacterium. Among these are 12 linear and nine circular plasmids, whose sequences total 610 694 bp. We report here the nucleotide sequence of three linear and seven circular plasmids (comprising 290 546 bp) in this infectious isolate. This completes the genome sequencing project for this organism; its genome size is 1 521 419 bp (plus about 2000 bp of undetermined telomeric sequences). Analysis of the sequence implies that there has been extensive and sometimes rather recent DNA rearrangement among a number of the linear plasmids. Many of these events appear to have been mediated by recombinational processes that formed duplications. These many regions of similarity are reflected in the fact that most plasmid genes are members of one of the genome's 161 paralogous gene families; 107 of these gene families, which vary in size from two to 41 members, contain at least one plasmid gene. These rearrangements appear to have contributed to a surprisingly large number of apparently non-functional pseudogenes, a very unusual feature for a prokaryotic genome. The presence of these damaged genes suggests that some of the plasmids may be in a period of rapid evolution. The sequence predicts 535 plasmid genes >/=300 bp in length that may be intact and 167 apparently mutationally damaged and/or unexpressed genes (pseudogenes). The large majority, over 90%, of genes on these plasmids have no convincing similarity to genes outside Borrelia, suggesting that they perform specialized functions.

Multiple paralogous genes related to the Streptomyces coelicolor developmental regulatory gene whiB are present in Streptomyces and other actinomycetes

The whiB sporulation gene of Streptomyces coelicolor was shown [Davis, N. K. & Chater, K. F. (1992). Mol Gen Genet 232, 351-358] to encode a small, cysteine-rich putative transcription factor unlike any that had been described previously. The large database of DNA sequences of mycobacteria (like Streptomyces, members of the Actinomycetales) has revealed a family of genes encoding proteins related to WhiB. Mycobacterium tuberculosis contains at least six such genes (whiB homologues in mycobacteria: whmA-F) and a likely seventh, whmG. Using conserved features of Whm proteins, a PCR-based approach led to the discovery that S. coelicolor A3(2) contains several similar genes. Cloning and sequencing of these whiB-like (wbI) genes revealed likely orthologues of four of the whm genes of M. tuberculosis. In all, S. coelicolor contains at least five wbI genes in addition to whiB itself. All five were shown by RT-PCR to be transcribed. A Southern blotting survey using each wbI gene as a probe showed that nearly all of a series of representatives of ten actinomycete genera (including morphologically simple organisms) contain close homologues of several wbI genes, suggesting that the ancient progenitor of all these organisms already contained a family of such genes, which have not been found in any other organisms.

Insights into the evolutionary process of genome degradation

Studies of noncoding and pseudogene sequence diversity, particularly in Rickettsia, have begun to reveal the basic principles of genome degradation in microorganisms. Increasingly, studies of genes and genomes suggest that there has been an extensive amount of horizontal gene transfer among microorganisms. As this inflow of genetic material does not seem generally to have resulted in genome size expansions, however, degenerative processes must be at the very least as widespread as horizontal gene transfer. The basic principles of gene degradation and elimination that are being explored in Rickettsia are likely to be of major importance for our understanding of how microbial genomes evolve.

Helicase-defective RuvB(D113E) promotes RuvAB-mediated branch migration in vitro

In Escherichia coli, the RuvA and RuvB proteins interact at Holliday junctions to promote branch migration leading to the formation of heteroduplex DNA. RuvA provides junction-binding specificity and RuvB drives ATP-dependent branch migration. Since RuvB contains sequence motifs characteristic of a DNA helicase and RuvAB exhibit helicase activity in vitro, we have analysed the role of DNA unwinding in relation to branch migration. A mutant RuvB protein, RuvB(D113E), mutated in helicase motif II (the DExx box), has been purified to homogeneity. The mutant protein forms hexameric rings on DNA similar to those formed by wild-type protein and promotes branch migration in the presence of RuvA. However, RuvB(D113E) exhibits reduced ATPase activity and is severely compromised in its DNA helicase activity. Models for RuvAB-mediated branch migration that invoke only limited DNA unwinding activity are proposed.

Learning from the genome sequence of Mycobacterium tuberculosis H37Rv

Mycobacterium tuberculosis, the scourge of humanity, is one of the most successful and scientifically challenging pathogens of all time. To catalyse the conception of new prophylactic and therapeutic interventions against tuberculosis, and to enhance our understanding of the biology of the tubercle bacillus, the complete genome sequence of the most widely used strain, H37Rv, has been determined. Bioinformatic analysis led to the identification of approximately 4000 genes in the 4.41 Mb genome sequence and provided fresh insight into the biochemistry, physiology. genetics and immunology of this much-feared bacterium. Genomic information is centralised in TubercuList (http://www.pasteur.fr/Bio/TubercuList/).

A gonococcal porA pseudogene: implications for understanding the evolution and pathogenicity of Neisseria gonorrhoeae

Members of the genus Neisseria, including the human pathogens Neisseria meningitidis and Neisseria gonorrhoeae, express at least one member of a family of related porins. N. meningitidis is the only species known to express a second porin, the meningococcal serosubtyping antigen PorA, the most divergent member of this family. Unexpectedly, a porA gene was identified in the gonococcal genome. Both the gonococcal and meningococcal porA loci were adjacent to a homologue of the Escherichia coli greA gene, although the IS1106 element downstream of porA in some meningococci was absent in the gonococcus. Almost identical porA loci were present in four unrelated gonococcal isolates and clinical specimens from patients with gonorrhoea. Lack of PorA expression in the gonococcus resulted from mutations in the promoter region, which prevented transcription, and frameshift mutations in the coding region of the porA gene. Hybridization and amplification experiments, showing the absence of a porA gene in seven other Neisseria species, suggested that porA was acquired by a common ancestor of the gonococcus and meningococcus but inactivated in the gonococcus on speciation. This implies that, while advantageous during colonization of the upper respiratory tract, this protein has no function in, or hinders, colonization of the urogenital tract.

Comparative mycobacterial genomics

Genomics is providing us with a mass of information about the biochemistry, physiology and pathogenesis of Mycobacterium tuberculosis and Mycobacterium leprae. Comparison of the two genome sequences is mutually enriching and indicates that the M. leprae genome appears to have undergone shrinkage and large-scale gene inactivation, which may account for the exceptionally slow growth of this organism.

Novel approaches to the discovery of antimicrobial agents

With the completion of numerous bacterial genome sequences, the discovery of antibacterial drugs has fully entered the genomic era. The strategies for effectively using genomic information for target identification, target characterization, screen development and compound evaluation are emerging, and have greatly increased the number of antibacterial targets available for screening. Fortunately, simultaneous efforts in improving miniaturization, robotics and database tools are underway so that the potential of genomics can be realized.

Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence

Countless millions of people have died from tuberculosis, a chronic infectious disease caused by the tubercle bacillus. The complete genome sequence of the best-characterized strain of Mycobacterium tuberculosis, H37Rv, has been determined and analysed in order to improve our understanding of the biology of this slow-growing pathogen and to help the conception of new prophylactic and therapeutic interventions. The genome comprises 4,411,529 base pairs, contains around 4,000 genes, and has a very high guanine + cytosine content that is reflected in the biased amino-acid content of the proteins. M. tuberculosis differs radically from other bacteria in that a very large portion of its coding capacity is devoted to the production of enzymes involved in lipogenesis and lipolysis, and to two new families of glycine-rich proteins with a repetitive structure that may represent a source of antigenic variation.

Genomics and computational molecular biology

There has been a dramatic increase in the number of completely sequenced bacterial genomes during the past two years as a result of the efforts both of public genome agencies and the pharmaceutical industry. The availability of completely sequenced genomes permits more systematic analyses of genes, evolution and genome function than was otherwise possible. Using computational methods - which are used to identify genes and their functions including statistics, sequence similarity, motifs, profiles, protein folds and probabilistic models - it is possible to develop characteristic genome signatures, assign functions to genes, identify pathogenic genes, identify metabolic pathways, develop diagnostic probes and discover potential drug-binding sites. All of these directions are critical to understanding bacterial growth, pathogenicity and host-pathogen interactions.

Complete genome sequence of Methanobacterium thermoautotrophicum deltaH: functional analysis and comparative genomics

The complete 1,751,377-bp sequence of the genome of the thermophilic archaeon Methanobacterium thermoautotrophicum deltaH has been determined by a whole-genome shotgun sequencing approach. A total of 1,855 open reading frames (ORFs) have been identified that appear to encode polypeptides, 844 (46%) of which have been assigned putative functions based on their similarities to database sequences with assigned functions. A total of 514 (28%) of the ORF-encoded polypeptides are related to sequences with unknown functions, and 496 (27%) have little or no homology to sequences in public databases. Comparisons with Eucarya-, Bacteria-, and Archaea-specific databases reveal that 1,013 of the putative gene products (54%) are most similar to polypeptide sequences described previously for other organisms in the domain Archaea. Comparisons with the Methanococcus jannaschii genome data underline the extensive divergence that has occurred between these two methanogens; only 352 (19%) of M. thermoautotrophicum ORFs encode sequences that are >50% identical to M. jannaschii polypeptides, and there is little conservation in the relative locations of orthologous genes. When the M. thermoautotrophicum ORFs are compared to sequences from only the eucaryal and bacterial domains, 786 (42%) are more similar to bacterial sequences and 241 (13%) are more similar to eucaryal sequences. The bacterial domain-like gene products include the majority of those predicted to be involved in cofactor and small molecule biosyntheses, intermediary metabolism, transport, nitrogen fixation, regulatory functions, and interactions with the environment. Most proteins predicted to be involved in DNA metabolism, transcription, and translation are more similar to eucaryal sequences. Gene structure and organization have features that are typical of the Bacteria, including genes that encode polypeptides closely related to eucaryal proteins. There are 24 polypeptides that could form two-component sensor kinase-response regulator systems and homologs of the bacterial Hsp70-response proteins DnaK and DnaJ, which are notably absent in M. jannaschii. DNA replication initiation and chromosome packaging in M. thermoautotrophicum are predicted to have eucaryal features, based on the presence of two Cdc6 homologs and three histones; however, the presence of an ftsZ gene indicates a bacterial type of cell division initiation. The DNA polymerases include an X-family repair type and an unusual archaeal B type formed by two separate polypeptides. The DNA-dependent RNA polymerase (RNAP) subunits A', A", B', B" and H are encoded in a typical archaeal RNAP operon, although a second A' subunit-encoding gene is present at a remote location. There are two rRNA operons, and 39 tRNA genes are dispersed around the genome, although most of these occur in clusters. Three of the tRNA genes have introns, including the tRNAPro (GGG) gene, which contains a second intron at an unprecedented location. There is no selenocysteinyl-tRNA gene nor evidence for classically organized IS elements, prophages, or plasmids. The genome contains one intein and two extended repeats (3.6 and 8.6 kb) that are members of a family with 18 representatives in the M. jannaschii genome.