The evolutionary relationships between the two bacteria Escherichia coli and Haemophilus influenzae and their putative last common ancestor

We have tried to approach the nature of the last common ancestor to Haemophilus influenzae and Escherichia coli and to determine how each bacterium could have diverged from this putative organism. The approach used was exhaustive analysis of the homologous proteins coded by genes present in these bacteria, using as criteria for sequence relatedness an alignment of at least 80 amino acid residues and a PAM distance (number of accepted point mutations per 100 residues separating two sequences) below 250. Evolutionarily significant similarities were found between 1,345 H. influenzae proteins (85% of the total genome) and 3,058 E. coli. proteins (75% of the total genome), many of them belonging to families of various sizes (from 666 doublets to 35 large groups of more than 10 members). Nearly all the genes found by this approach to be duplicated in both bacteria were already duplicated in their last common ancestor. This was deduced from (1) the comparison of the respective distributions of evolutionary distances between orthologs (genes separated only by speciation events) and paralogs (genes duplicated in the same genome) and (2) the analysis of the phylogenetic trees reconstructed for each family of paralogs containing at least two members belonging to each bacterium. The distributions of the different categories of homologs show a significant loss of paralogous genes in H. influenzae (reduction proportional to the genome size), of many sequences which are still present in one copy in E. coli, and of some entire gene families. Phylogenetic trees also confirmed this recent loss of paralogous genes in H. influenzae. Thus, the genome size of the last common ancestor of these two bacteria would have been close to that of present-day E. coli, and the evolution of H. influenzae toward a parasitic life led to an important decrease in its genome size by some mechanism of streamlining. During this recent evolution, the memory of the gene order present in the last common ancestor has been blurred, but a few short conserved chromosomal fragments can still be detected in present-day E. coli and H. influenzae.      

Site-specific recombination between cloned attP and attB sites from the Haemophilus influenzae bacteriophage HP1 propagated in recombination-deficient Escherichia coli.

Plasmids were constructed which contain both attP and attB DNA segments derived from the insertion sites of the lysogenic bacteriophage HP1 and its host, Haemophilus influenzae. Similar plasmids containing the two junction segments (attL and attR regions) between the phage genome and the lysogenic host chromosome were also prepared. The formation of recombinant dimer plasmids was observed when attP-attB plasmids were propagated in Escherichia coli HB101 (recA), while plasmids containing the junction segments did not form recombinant dimers. Deletion of the phage DNA segment adjacent to the attP site from the attP-attB constructions eliminated detectable recombination, suggesting that this sequence contains the gene encoding the HP1 integrase. No plasmid recombination was observed in strains of E. coli defective in integration host factor. This suggests that integration host factor is important in the expression or activity of the system which produces the site-specific recombination of sequences derived from HP1 and H. influenzae. Further, it suggests that a protein functionally analogous to E. coli integration host factor may be present in H. influenzae.     

A peroxide/ascorbate-inducible catalase from Haemophilus influenzae is homologous to the Escherichia coli katE gene product.

Bacterial catalases are induced by exposure to peroxide (e.g., Escherichia coli katG) or entry into stationary phase (e.g., E. coli katE). To study regulatory systems in Haemophilus influenzae, we complemented an E. coli rpoS mutant, which is unable to induce katE in stationary phase, with a plasmid library of H. influenzae Rd- chromosomal DNA. Nineteen complementing clones with a catalase-positive phenotype were obtained and characterized after screening about 10(5) transformants. All carried the same structural gene for an H. influenzae catalase. The DNA sequence of this gene, called hktE, encodes a 508-amino-acid polypeptide with strong homology to eukaryotic catalases and E. coli katE. However, hktE is regulated like E. coli katG, with catalase activity increasing 10-fold and hktE mRNA levels increasing 4-fold upon exposure to ascorbic acid, which serves to generate hydrogen peroxide. Mutations in the known global regulatory genes of H. influenzae--crp, cya, and sxy--do not affect the inducibility of hktE. The hktE gene maps to a 225-kb segment of the H. influenzae chromosome in a region encoding resistance to spectinomycin.     

Cloning and expression in Escherichia coli of Haemophilus influenzae fimbrial genes establishes adherence to oropharyngeal epithelial cells.

In this report the first example of functional expression of a fimbrial gene cluster of a non-enteric human pathogen in Escherichia coli is described. This is shown for Haemophilus influenzae fimbriae which mediate adherence to oropharyngeal epithelial cells. A genomic library of H.influenzae type b, strain 770235f+bo, was constructed using a cosmid vector and screened with a synthetic oligonucleotide probe derived from the N-terminal sequence of the fimbrial subunit of H.influenzae. Four cosmid clones were found which hybridized to this oligonucleotide probe. Escherichia coli strains harbouring these clones expressed the H.influenzae fimbriae at their cell surface, as was demonstrated in a whole-cell ELISA and by immunogold electron microscopy using a monoclonal antibody specific for the H.influenzae fimbriae. Surface expression could be maintained during subcloning until a minimal H.influenzae DNA insert of approximately 8.1 kb was obtained. Escherichia coli strains harbouring the 8.1 kb H. influenzae DNA were able to cause a mannose-resistant adherence to oropharyngeal epithelial cells and a mannose-resistant haemagglutination of human AnWj-positive erythrocytes. The nucleotide sequence of hifA, the gene encoding the major fimbrial subunit, was determined. The predicted amino acid sequence shows a significant homology with a number of E.coli fimbrial subunits.     

Haemophilus influenzae Rd lacks a stringently conserved fatty acid biosynthetic enzyme and thermal control of membrane lipid composition

The organization of the fatty acid synthetic genes of Haemophilus influenzae Rd is remarkably similar to that of the paradigm organism, Escherichia coli K-12, except that no homologue of the E. coli fabF gene is present. This finding is unexpected, since fabF is very widely distributed among bacteria and is thought to be the generic 3-ketoacyl-acyl carrier protein (ACP) synthase active on long-chain-length substrates. However, H. influenzae Rd contains a homologue of the E. coli fabB gene, which encodes a 3-ketoacyl-ACP synthase required for unsaturated fatty acid synthesis, and it seemed possible that the H. influenzae FabB homologue might have acquired the functions of FabF. E. coli mutants lacking fabF function are unable to regulate the compositions of membrane phospholipids in response to growth temperature. We report in vivo evidence that the enzyme encoded by the H. influenzae fabB gene has properties essentially identical to those of E. coli FabB and lacks FabF activity. Therefore, H. influenzae grows without FabF function. Moreover, as predicted from studies of the E. coli fabF mutants, H. influenzae is unable to change the fatty acid compositions of its membrane phospholipids with growth temperature. We also demonstrate that the fabB gene of Vibrio cholerae El Tor N16961 does not contain a frameshift mutation as was previously reported.     

Cloning of genes encoding a 15,000-dalton peptidoglycan-associated outer membrane lipoprotein and an antigenically related 15,000-dalton protein from Haemophilus influenzae.

We have cloned and expressed in Escherichia coli a gene encoding a 15,000-apparent-molecular-weight peptidoglycan-associated outer membrane lipoprotein (PAL) of Haemophilus influenzae. The nucleotide sequence of this gene encodes an open reading frame of 153 codons with a predicted mature protein of 134 amino acids. The amino acid composition and sequence of the predicted mature protein agree with the chemically determined composition and partial amino acid sequence of PAL purified from H. influenzae outer membranes. We have also identified a second gene from H. influenzae that encodes a second 15,000-apparent-molecular-weight protein which is recognized by antiserum against PAL. This protein has been shown to be a lipoprotein. The nucleotide sequence of this gene encodes an open reading frame of 154 codons with a predicted mature protein of 136 amino acids and has limited sequence homology with that of the gene encoding PAL. Southern hybridization analysis indicates that both genes exist as single copies in H. influenzae chromosomal DNA. Both genes encode polypeptides which have amino-terminal sequences similar to those of reported membrane signal peptides and are associated primarily with the outer membrane when expressed in E. coli.     

Isolation of Haemophilus influenzae genes that suppress Escherichia coli polA mutations

Haemophilus influenzae was found to produce a DNA polymerase that was similar to polymerase I of Escherichia coli. E. coli polA mutants were used as backgrounds for the selection of H. influenzae polA suppressor genes. Six different H. influenzae fragments were isolated that could suppress E. coli polA mutations. None of the suppressors appeared to encode the H. influenzae equivalent of the E. coli polA gene. One type of clone, represented by pGW41, caused a polymerase I activity to appear in a suppressed polA1 mutant. Plasmids from the pGW41 class contained two genes (pol-2 and pol-3) that were both required for polA suppression. Mutated nonsuppressing derivatives of the pGW41 class were used to create H. influenzae mutants that were deficient in polymerase I.     

Haemophilus influenzae immunoglobulin A1 protease genes: cloning by plasmid integration-excision, comparative analyses, and localization of secretion determinants.

Many bacteria which establish infections after invasion at human mucosal surfaces produce enzymes which cleave immunoglobulin A (IgA), the primary immunoglobulin involved with protection at these sites. Bacterial species such as Haemophilus influenzae which produce IgA1 proteases secrete this enzyme into their environment. However, when the gene encoding this protein was isolated from H. influenzae serotype d and introduced into Escherichia coli, the activity was not secreted into the medium but was localized in the periplasmic space. In this study, the IgA1 protease gene (iga) from an H. influenzae serotype c strain was isolated and the gene from the serotype d strain was reisolated. The IgA1 proteases produced in E. coli from these genes were secreted into the growth medium. A sequence linked to the carboxyl terminus of the iga gene but not present in the original clone was shown to be necessary to achieve normal secretion. Tn5 mutagenesis of the additional carboxyl-terminal region was used to define a 75- to 100-kilodalton coding region required for complete secretion of IgA1 protease but nonessential for protease activity. The iga genes were isolated by a plasmid integration-excision procedure. In this method a derivative of plasmid pBR322 containing a portion of the protease gene and the kanamycin resistance determinant of Tn5 was introduced into H. influenzae by transformation. The kanamycin resistance gene was expressed in H. influenzae, but since pBR322 derivatives are unable to replicate in this organism, kanamycin-resistant transformants arose by integration of the plasmid into the Haemophilus chromosome by homologous recombination. The plasmid, together with the adjoining DNA encoding IgA1 protease, was then excised from the chromosome with DNA restriction enzymes, religated, and reintroduced into E. coli. Comparisons between the H. influenzae protease genes were initiated which are useful in locating functional domains of these enzymes.     

Identification of the Chi site of Haemophilus influenzae as several sequences related to the Escherichia coli Chi site

The Escherichia coli Chi site 5'-GCTGGTGG-3' modulates the activity of the powerful dsDNA exonuclease and helicase RecBCD. Genome sequence analyses revealed that Chi is frequent on the chromosome and oriented with respect to replication on the E . coli genome. Chi is also present much more frequently than predicted statistically for a random 8-mer sequence. Although it is assumed that Chi is ubiquitous, there is virtually no proof that its features are conserved in other microorganisms. We therefore identified and analysed the Chi sequence of an organism for which the full genome sequence was available, Haemophilus influenzae . The biological test we used is based on our finding that rolling circle plasmids provide a specific substrate for RecBCD analogues in different microorganisms. Unexpectedly, several related sequences, corresponding to 5'-GNTGGTGG-3' and 5'-G(G/C)TGGAGG-3', showed Chi activity. As in E . coli , the H . influenzae Chi sites are frequent on the genome, which is in keeping with the need for frequent Chi sites for dsDNA break repair of chromosomal DNA. Although statistically over-represented, this feature is less marked than that of the E . coli Chi site. In contrast to E . coli , the H . influenzae Chi motifs are only slightly oriented with respect to the replication strand. Thus, although Chi appears to have a highly conserved biological role in attenuating exonuclease activity, its sequence characteristics and statistical representation on the genome may differ according to the particular features of the host.     

Haemophilus influenzae UvrA: overexpression, purification, and in cell complementation

UvrA protein is a major component of ABC endonuclease complex involved in nucleotide excision repair (NER) mechanism. Although NER system is best characterized in Escherichia coli, not much information is available in Haemophilus influenzae. However, based on amino acid homology, uvrA ORF has been identified on H. influenzae genome [gene identification No. HI0249, Science 269 (1995) 496]. H. influenzae Rd uvrA ORF was cloned and overexpressed in E. coli. The expressed UvrA protein was purified using a two-step column chromatography protocol to a single band of expected molecular weight (104 kDa) and characterized for its ATPase and DNA binding activity. In addition, when H. influenzae uvrA was introduced in E. coli uvrA mutant strain AB1886, its UV resistance was restored to near wild type level.     

A Haemophilus influenzae gene that encodes a membrane bound 3-deoxy-D-manno-octulosonic acid (Kdo) kinase. Possible involvement of kdo phosphorylation in bacterial virulence.

The lipopolysaccharide of Haemophilus influenzae contains a single 3-deoxy-D-manno-octulosonic acid (Kdo) residue derivatized with either a phosphate or an ethanolamine pyrophosphate moiety at the 4-OH position. In previous studies, we identified a kinase unique to H. influenzae extracts that phosphorylates Kdo-lipid IV(A), a key precursor of lipopolysaccharide in this organism. We have now identified the gene encoding the Kdo kinase by using an expression cloning approach. A cosmid library containing random DNA fragments from H. influenzae strain Rd was constructed in Escherichia coli. Extracts of 472 colonies containing individual hybrid cosmids were assayed for Kdo kinase activity. A single hybrid cosmid directing expression of the kinase was found. The kinase gene was identified by activity assays, sub-cloning, and DNA sequencing. When the putative kinase gene was expressed in E. coli behind a T7 promoter, massive overproduction of kinase activity was achieved ( approximately 8000-fold higher than in H. influenzae membranes). The catalytic properties and the product generated by the overexpressed kinase, assayed with Kdo-lipid IV(A) as the substrate, were the same as observed with H. influenzae membranes. Unexpectedly, the kinase gene was identical to a previously characterized open reading frame (orfZ), which had been shown to be important for establishing bacteremia in an infant rat model (Hood, D. W., Deadman, M. E., Allen, T., Masoud, H., Martin, A., Brisson, J. R., Fleischmann, R., Venter, J. C., Richards, J. C., and Moxon, E. R. (1996) Mol. Microbiol. 22, 951-965). However, based solely on the genome sequence of H. influenzae Rd, no biochemical function had been assigned to the product of orfZ, which we now designate kdkA ("Kdo kinase A"). Although Kdo phosphorylation may be critical for bacterial virulence of H. influenzae, it does not appear to be required for growth.     

Utilization of enterobactin and other exogenous iron sources by Haemophilus influenzae, H. parainfluenzae and H. paraphrophilus

The ability of Haemophilus influenzae, H. parainfluenzae and H. paraphrophilus to utilize iron complexes, iron-proteins and exogenous microbial siderophores was evaluated. In a plate bioassay, all three species used not only ferric nitrate but also the iron chelates ferric citrate, ferric nitrilotriacetate and ferric 2,3-dihydroxybenzoate. Each Haemophilus species examined also used haemin, haemoglobin and haem-albumin as iron sources although only H. influenzae could acquire iron from transferrin or from haemoglobin complexed with haptoglobin. None of the haemophili obtained iron from ferritin or lactoferrin or from the microbial siderophores aerobactin or desferrioxamine B. However, the phenolate siderophore enterobactin supplied iron to both H. parainfluenzae and H. paraphrophilus, and DNA isolated from both organisms hybridized with a DNA probe prepared from the Escherichia coli ferric enterobactin receptor gene fepA. In addition, a monospecific polyclonal antiserum raised against the E. coli 81 kDa ferric enterobactin receptor (FepA) recognized an iron-repressible outer membrane protein (OMP) in H. parainfluenzae of between 80 and 82 kDa (depending on the strain). This anti-FepA serum did not cross-react with any of the OMPs of H. paraphrophilus or H. influenzae. The OMPs of each Haemophilus species were also probed with antisera raised against the 74 kDa Cir or 74 kDa IutA (aerobactin receptor) proteins of E. coli. Apart from one H. parainfluenzae strain (NCTC 10665), in which an OMP of about 80 kDa cross-reacted with the anti-IutA sera, no cross-reactivity was observed between Cir, IutA and the OMPs of H. influenzae, H. parainfluenzae or H. paraphrophilus.(ABSTRACT TRUNCATED AT 250 WORDS)     

Functional characterization of the DNA mismatch binding protein MutS from Haemophilus influenzae

This investigation demonstrates DNA mismatch repair activity in Haemophilus influenzae cell free extracts. The mutS gene as well as purified protein of H. influenzae restored repair activity in complementation assays performed with mutS deficient Escherichia coli strain. The difference in affinity for GT and AC mismatched bases by H. influenzae MutS was reflected in the efficiency with which these DNA heteroduplexes were repaired in vitro, with GT being repaired well and AC the least. Unlike E. coli MutS, the H. influenzae homolog failed to give protein-DNA complex with homoduplex DNA. Interestingly, MutS was found to bind single-stranded DNA but with lesser affinity as compared to heteroduplex DNA. Apart from the nucleotide- and DNA-mediated conformational transitions, as monitored by circular dichroism and limited proteolysis, our data suggest a functional role when H. influenzae MutS encounters single-stranded DNA during exonucleolytic step of DNA repair process. We propose that, conformational changes in H. influenzae MutS not only modulate mismatch recognition but also trigger some of the down stream processes involved in the DNA mismatch repair process.     

Mobilization of Haemophilus influenzae chromosomal markers by an Escherichia coli F'' factor.

Filter matings between E. coli K-12 strains carrying an F'::Tn5,Tn9 factor with H. influenzae Rd strains gave rise to kanamycin-chloramphenicol-resistant H. influenzae strains at a frequency of approximately 10(-6). Transfer of the F' factor to H. influenzae was verified by expression of unselected markers in H. influenzae (lac+ or cotransfer of the nonselected antibiotic resistance), physical presence of a high-molecular-weight plasmid in recipient H. influenzae cells, and detection by Southern hybridization analysis of DNA sequences specific for the F' factor replication and partition functions in recipient H. influenzae cells. H. influenzae (F' Tn5,Tn9) strains were capable of transferring kanamycin and chloramphenicol resistances to other H. influenzae strains and were capable of mobilizing H. influenzae chromosomal markers at a low frequency. Insertion of a Tn5 element in the H. influenzae genome near the novobiocin resistance gene increased the frequency of transfer of novobiocin resistance about 30-fold. Transfer of other chromosomal markers also increased, although to a lesser extent, and ordered transfer of chromosomal markers could be demonstrated. Gene transfer was insensitive to DNase I, and transfer of chromosomal (but not F' factor) markers was dependent on the H. influenzae rec-1 and rec-2 gene functions.     

Global versus local regulatory roles for Lrp-related proteins: Haemophilus influenzae as a case study

Lrp (leucine-responsive regulatory protein) plays a global regulatory role in Escherichia coli, affecting expression of dozens of operons. Numerous lrp-related genes have been identified in different bacteria and archaea, including asnC, an E. coli gene that was the first reported member of this family. Pairwise comparisons of amino acid sequences of the corresponding proteins shows an average sequence identity of only 29% for the vast majority of comparisons. By contrast, Lrp-related proteins from enteric bacteria show more than 97% amino acid identity. Is the global regulatory role associated with E. coli Lrp limited to enteric bacteria? To probe this question we investigated LrfB, an Lrp-related protein from Haemophilus influenzae that shares 75% sequence identity with E. coli Lrp (highest sequence identity among 42 sequences compared). A strain of H. influenzae having an lrfB null allele grew at the wild-type growth rate but with a filamentous morphology. A comparison of two-dimensional (2D) electrophoretic patterns of proteins from parent and mutant strains showed only two differences (comparable studies with lrp(+) and lrp E. coli strains by others showed 20 differences). The abundance of LrfB in H. influenzae, estimated by Western blotting experiments, was about 130 dimers per cell (compared to 3,000 dimers per E. coli cell). LrfB expressed in E. coli replaced Lrp as a repressor of the lrp gene but acted only to a limited extent as an activator of the ilvIH operon. Thus, although LrfB resembles Lrp sufficiently to perform some of its functions, its low abundance is consonant with a more local role in regulating but a few genes, a view consistent with the results of the 2D electrophoretic analysis. We speculate that an Lrp having a global regulatory role evolved to help enteric bacteria adapt to their ecological niches and that it is unlikely that Lrp-related proteins in other organisms have a broad regulatory function.     

Environmental and genetic regulation of the phosphorylcholine epitope of Haemophilus influenzae lipooligosaccharide

In response to environmental signals in the host, bacterial pathogens express factors required during infection and repress those that interfere with specific stages of this process. Signalling pathways controlling virulence factors of the human respiratory pathogen, Haemophilus influenzae, are predominantly unknown. The lipooligosaccharide (LOS) outer core represents a prototypical virulence trait of H. influenzae that enhances virulence but also provides targets for innate and adaptive immunity. We report regulation of the display of the virulence-associated phosphorylcholine (PC) epitope on the LOS in response to environmental conditions. PC display is optimal under microaerobic conditions and markedly decreased under conditions of high culture aeration. Gene expression analysis using a DNA microarray was performed to begin to define the metabolic state of the cell under these conditions and to identify genes potentially involved in PC epitope modulation. Global gene expression profiling detected changes in redox responsive genes and in genes of carbohydrate metabolism. The effects on carbohydrate metabolism led us to examine the role of the putative H. influenzae homologue of csrA, a regulator of glycolysis and gluconeogenesis in Escherichia coli. A mutant containing an in-frame deletion of the H. influenzae csrA gene showed increased PC epitope levels under aerobic conditions. Furthermore, deletion of csrA elevated mRNA expression of galU, an essential virulence gene that is critical in generating sugar precursors needed for polysaccharide formation and LOS outer core synthesis. Growth conditions predicted to alter the redox state of the culture modulated the PC epitope and galU expression as well. The results are consistent with a multifactorial mechanism of control of LOS-PC epitope display involving csrA and environmental signals that coordinately regulate biosynthetic and metabolic genes controlling the LOS structure.     

Comparison of synonymous codon distribution patterns of bacteriophage and host genomes.

Synonymous codon usage patterns of bacteriophage and host genomes were compared. Two indexes, G + C base composition of a gene (fgc) and fraction of translationally optimal codons of the gene (fop), were used in the comparison. Synonymous codon usage data of all the coding sequences on a genome are represented as a cloud of points in the plane of fop vs. fgc. The Escherichia coli coding sequences appear to exhibit two phases, "rising" and "flat" phases. Genes that are essential for survival and are thought to be native are located in the flat phase, while foreign-type genes from prophages and transposons are found in the rising phase with a slope of nearly unity in the fgc vs. fop plot. Synonymous codon distribution patterns of genes from temperate phages P4, P2, N15 and lambda are similar to the pattern of E. coli rising phase genes. In contrast, genes from the virulent phage T7 or T4, for which a phage-encoded DNA polymerase is identified, fall in a linear curve with a slope of nearly zero in the fop vs. fgc plane. These results may suggest that the G + C contents for T7, T4 and E. coli flat phase genes are subject to the directional mutation pressure and are determined by the DNA polymerase used in the replication. There is significant variation in the fop values of the phage genes, suggesting an adjustment to gene expression level. Similar analyses of codon distribution patterns were carried out for Haemophilus influenzae, Bacillus subtilis, Mycobacterium tuberculosis and their phages with complete genomic sequences available.