Previous Article | Next Article 
Journal of Bacteriology, August 2001, p. 4626-4635, Vol. 183, No. 15
0021-9193/01/$04.00+0 DOI: 10.1128/JB.183.15.000-000.2001
Copyright © 2001, American Society for Microbiology. All rights reserved.
Evolution of an Autotransporter: Domain Shuffling
and Lateral Transfer from Pathogenic Haemophilus to
Neisseria
Jeamelia
Davis,1
Arnold L.
Smith,1
William R.
Hughes,2 and
Miriam
Golomb1,*
Division of Biological Sciences and
Department of Molecular Microbiology and Immunology, School of
Medicine, University of Missouri
Columbia, Columbia,
Missouri,1 and The Genome Center,
University of Washington, Seattle, Washington2
Received 15 February 2001/Accepted 20 April 2001
 |
ABSTRACT |
The genomes of pathogenic Haemophilus influenzae
strains are larger than that of Rd KW20 (Rd), the nonpathogenic
laboratory strain whose genome has been sequenced. To identify
potential virulence genes, we examined genes possessed by Int1, an
invasive nonencapsulated isolate from a meningitis patient, but absent from Rd. Int1 was found to have a novel gene termed lav,
predicted to encode a member of the AIDA-I/VirG/PerT family of
virulence-associated autotransporters (ATs). Associated with
lav are multiple repeats of the tetranucleotide GCAA,
implicated in translational phase variation of surface molecules.
Laterally acquired by H. influenzae, lav is restricted in
distribution to a few pathogenic strains, including H. influenzae biotype aegyptius and Brazilian purpuric fever
isolates. The DNA sequence of lav is surprisingly similar to that of a gene previously described for Neisseria
meningitidis. Sequence comparisons suggest that lav
was transferred relatively recently from Haemophilus to
Neisseria, shortly before the divergence of N. meningitidis and Neisseria gonorrhoeae. Segments of
lav predicted to encode passenger and
-domains differ
sharply in G+C base content, supporting the idea that AT genes have
evolved by fusing domains which originated in different genomes.
Homology and base sequence comparisons suggest that a novel biotype
aegyptius AT arose by swapping an unrelated sequence for the passenger
domain of lav. The unusually mobile lav locus
joins a growing list of genes transferred from H. influenzae to Neisseria. Frequent gene exchange
suggests a common pool of hypervariable contingency genes and may help
to explain the origin of invasiveness in certain respiratory pathogens.
 |
INTRODUCTION |
Horizontal acquisition of virulence-determining
genes has accelerated the evolution of bacterial pathogens
(35). The recent availability of complete microbial genome
sequences facilitates investigation of the history of virulence
genes. Of particular interest are contingency loci, which encode
phase-variable surface molecules involved in host tissue interactions
(32). Contingency gene expression is commonly modulated by
slipped-strand mispairing acting on simple sequence repeats within a
control region. Although distantly related, the respiratory commensals
and pathogens Haemophilus influenzae and Neisseria
meningitidis each use tetranucleotide repeat sequences in
translational phase variation of surface proteins (23, 37,
40). Such similarity prompts the question of whether repeat-regulated genes were present in a common ancestor, evolved independently, or had a recent common origin in laterally transferred DNA.
Several important human pathogens (Streptococcus pneumoniae, H. influenzae, and N. meningitidis) coexist in the
respiratory flora of healthy individuals. Kroll et al., who first
reported natural transfer of H. influenzae sequences to
N. meningitidis, have suggested that genetic exchange
between respiratory pathogens may spark the emergence of new invasive
strains (25). They have identified three sequences,
originally from Haemophilus and flagged by
Haemophilus-specific uptake sequences (hUSs), which now
reside in the meningococcal genome.
H. influenzae, a small, gram-negative bacterium, causes
otitis media and bronchitis as well as invasive disease (meningitis and
septicemia) (50). Invasive H. influenzae
disease is usually caused by encapsulated strains of H. influenzae serotype b (Hib); near-universal immunization with
conjugate Hib vaccine has largely eliminated Haemophilus
meningitis from developed countries (3). Nonencapsulated
(nontypeable) H. influenzae (NTHi), against which the Hib
vaccine provides no protection, remains an important cause of
respiratory infections (18). Although NTHi rarely causes invasive infection in immunocompetent hosts, sporadic exceptions warn
of potential virulence. During the 1980s an outbreak of highly lethal
septicemia among children in rural Brazil, termed Brazilian purpuric
fever (BPF), was traced to a single NTHi clone related to H. influenzae biotype aegyptius strains previously associated only
with conjunctivitis (7). In 1994, an NTHi strain (Int1; also called R2866) was isolated in the United States from the blood of
a previously healthy, Hib-immunized child who developed meningitis
(34).
H. influenzae is readily transformed with
Haemophilus DNA; uptake requires the hUS, consisting of a
conserved 9-bp core within an extended 29-bp consensus sequence. The
H. influenzae strain Rd KW20 (hereafter referred to as Rd)
genome, which has been sequenced (17), contains 1,465 hUSs. While most of these occur as single copies, 17% occur as pairs
of inverted repeats in a stem-loop configuration. Stem-loop hUS pairs
tend to occur at the ends of transcription units and may function as
transcription terminators (44). At least one horizontally
transferred island has been found inserted within a stem-loop hUS pair,
suggesting that paired hUSs may play an additional role as potential
targets for the insertion of virulence-associated genes, possibly
mediated by phage integrases (10, 30).
N. meningitidis is a gram-negative bacterium
phylogenetically distant from H. influenzae; N. meningitidis
is classified in the
subdivision and H. influenzae is
classified in the
subdivision of the proteobacteria. In sub-Saharan
Africa, N. meningitidis serogroup A causes epidemic
meningitis, whereas in the developed world serogroups B and C cause
endemic invasive disease (29, 47). Natural transformation
is the primary means of genetic exchange among neisseriae; it is so
frequent that lineage boundaries are blurred, requiring populations to
be depicted as networks rather than as distinct clades
(45). Transformation employs a neisseria-specific 10-bp
uptake sequence (nUS) distinct from the hUS (14).
To identify NTHi pathogenicity genes, we initiated whole genome
comparisons between Int1 and Rd, a nonpathogenic laboratory strain.
Here we report the identification of lav, a candidate pathogenicity gene whose inferred translation product is homologous to
virulence-associated autotransporters (AT) and whose DNA sequence is
similar to NMB1527 (also called orf2 and nmrep3)
and NMA1725 of N. meningitidis (24, 36, 37,
48). We present evidence that lav arose by fusion of
segments from different organisms and has recently been transferred
from Haemophilus to Neisseria.
 |
MATERIALS AND METHODS |
Bacteria.
H. influenzae was grown as previously
described (30, 31); strains are described in Table
1.
Differential cloning of Int1 sequences.
An M13mp18 library
was prepared from a partial Sau3A digest of Int1 DNA size
fractionated to yield fragments between 0.75 and 1.5 kb. Bacterial
clones containing recombinant plasmids were arrayed on nylon filters
(1,200 colonies each). Replicate filters were probed at high stringency
with digoxigenin-labeled Rd and Int1 genomic DNAs, using M13 clones
from Pseudomonas aeruginosa, a G+C-rich organism as
background controls. Clones hybridizing with Int1 but not with Rd were
picked, and inserts were initially sequenced in one direction with a
universal primer.
DNA isolation and PCR.
H. influenzae DNA
isolation and long PCR were done as previously described
(31). The holB-tmk region was
amplified with forward primer 5'-CTTTAATCAGCACAGCATGATGCC,
corresponding to H. influenzae Rd nucleotides 477248 to 477272, and reverse primer 5'-TGTTCAAGCTGATATTGAAAGTGC,
corresponding to H. influenzae nucleotides 477396 to
477373 (17). DNA from BPF isolates R2140 and R2141 was
amplified using the same reverse primer but with forward primer 5'-CACATTTTCAAACTGGCTTGAC.
DNA blotting.
The genomic blots, hybridization, procedures
and stringent wash conditions used have been previously described
(31). The lav probe was a gel-purified 1.4-kb
PCR fragment made using Int1 lav internal primers
5'-GCTTTTGGCTGTTGATTACG and 5'-CCGCCCATTAAGCCAACGG, that was then digoxigenin labeled.
DNA sequencing.
PCR fragments were gel purified and
sequenced directly or following subcloning into a phagemid vector; both
strands were sequenced as previously described (31).
Sequences were analyzed using BLASTX, and BLASTN (1, 2)
and the GCG and Omiga packages.
Nucleotide sequence accession number.
The sequences reported
in this paper have been deposited in GenBank and assigned accession no.
AF385403 (lav) and AF385404 (las).
 |
RESULTS |
Identification of lav.
We constructed a bank of
clones from the pathogenic NTH: strain Int1 and screened them by
hybridization against whole genomic DNA from strain Int1 and, as a
comparison, from nonpathogenic strain Rd. This differential screening
was intended to identify genes specifically present in Int1 which
potentially encode pathogenicity-related proteins. Of an initial 1,200 clones prepared from Int1 DNA, ~10% were novel relative to Rd.
Sequence 122 bore the 3' end of an open reading frame (ORF) homologous
to a family of AT, which includes VirG of Shigella flexneri
(27), AIDA-I of Escherichia coli
(4), VapA/nmrep2 of N. meningitidis
(37), and PerT of Bordetella pertussis
(15). These virulence factors are outer membrane proteins involved in adhesion, invasion, intercellular spread, or immune evasion
(21, 28). Each consists of a carboxyl-terminal
-barrel domain, which forms a pore in the outer membrane; a linker peptide; a
passenger effector domain, which is secreted through the pore; and an
N-terminal signal sequence. The partial Int1 ORF was ~90% identical
in DNA sequence to orf2/nmrep3/NMB1525 of N. meningitidis strain MC58 (serogroup B), previously identified by
its homology to other AT (24, 37).
Downstream from the presumed AT gene, sequence 122 contained a junction
with 93% DNA homology to a gene in strain Rd, the
3' end of
H. influenzae gene HI 0456 (
tmk), encoding thymidylate
kinase. In Rd, the gene immediately downstream from
tmk is
H. influenzae gene HI 0455 (
holB), which encodes
the

subunit of
DNA polymerase III. The initiator codon of
holB overlaps the terminator
codon of
tmk; a
paired hUS in (+/

) inverted repeat configuration
extends across the
region of overlap (Fig.
1A).

View larger version (49K):
[in this window]
[in a new window]
|
FIG. 1.
(A) Arrangement of lav homologs in H. influenzae and N. meningitidis MC58 genomes. Neisserial
genes are depicted as filled arrows (indicating orientation); H. influenzae-derived genes are depicted as open arrows. Hairpin
symbols represent paired hUSs in stem-loop configuration. Relative to
that in H. influenzae-Rd, the Int1 lav island is
inserted within a stem-loop hUS pair, resulting in partial duplication
of the hUS pair flanking the site of insertion within an Int1 ancestor.
Fragments of H. influenzae holB and tmk in
Neisseria are indicated in parentheses. (B)
Repeat-containing regions of lav homologs in Int1 (hi
lav), N. meningitidis MC58 (nm) (24,
37), N. gonorrhoeae, (ng) and biotype aegyptius (hae)
reference strain R1967. (C) Inferred amino acid sequences of regions of
Lav homologs, aligned with MEGALIGN (DNAStar). aa 1 to 40 are conserved
sequences at the N terminus, including the predicted signal sequence
cleavage site (arrow). aa 121 to 156 represent a less conserved region
within the presumptive passenger domain, including a sequence bracketed
by conserved cysteine residues (asterisks).
|
|
Rd-specific primers to the 5' end of
holB and the 3' end of
tmk were used to amplify Int1 genomic DNA by PCR. The
predicted
125-bp fragment was obtained with Rd and a Hib (Eagan)
template,
but a 2.4-kb fragment was obtained with Int1, including 2.3 kb
of novel sequence (Fig.
1A). The novel DNA is located within the
(+/

) hUS pair at the
holB-tmk junction, 18 bp of which are
now
found duplicated at either end of the island. The insert contains
a
single, 2,078-bp potential ORF termed
lav (for like a VirG).
The putative initiating ATG is in a different reading frame from
the
downstream coding region; 10 bp downstream are 19 copies of
the
tetranucleotide GCAA. A gain of one GCAA repeat creates a
2,082-bp ORF,
potentially encoding a 693-amino-acid (aa) Lav polypeptide
(pI = 9.3) with a 54-aa signal sequence. Lav is homologous at
its carboxyl
terminus to VirG and other
AT.
Horizontal transfer of lav.
The complete DNA
sequence of H. influenzae lav was 89% similar to
orf2/nmrep3 of N. meningitidis MC58. The
meningococcal gene is located in a region containing rfaF
(encoding a heptosyl transferase involved in lipo-oligosaccharide
biosynthesis) between a gene (orf1) whose inferred product
is a small protein B homolog and a gene (dld) encoding a
D-lactate dehydrogenase (4, 18, 48). N. meningitidis serotype B orf2 is part of a 2.45-kb
island of similarity to H. influenzae (Fig. 1A). Four GCAA
repeats follow the initiator codon of orf2/nmrep3, placing
its translated product out of reading frame; loss of a GCAA would
permit translation of a 678-aa protein. A homologous gene (NMA1725)
95% identical to N. meningitidis serotype B. orf2, but with three GCAA repeats, is found at the same site
in the genome of N. meningitidis serogroup A strain Z2491
(36).
The unfinished
Neisseria gonorrhoeae genome
(
http://www.genome.ou.edu) was searched for homology with
N. meningitidis orf2.
A gonococcal sequence 91% identical
to
N. meningitidis orf2 and
89% identical to
H. influenzae lav was found on contig 140, at
the same chromosomal
site as in
N. meningitidis. Endpoints of
the island of
homology with
H. influenzae are identical

exactly
57 bp
downstream from
orf1 and 35 bp downstream from
dld
in the
N. meningitidis and
N. gonorrhoeae genomes.
N. gonorrhoeae orf2 is preceded by
two GCAA repeats after the initiating codon (Fig.
1B), is also out of
reading frame, and contains an in-frame stop
codon and a frameshift
mutation. Inferred proteins encoded by
lav homologs contain
highly conserved regions interspersed with
divergent regions (Fig.
1C).
The N-terminal region surrounding
the signal sequence is highly
conserved, as is the C-terminal
end of the

-domain; the passenger
domain contains two conserved
cysteine residues which can potentially
form a short loop of variable
sequence (Fig.
1C).
Three considerations indicate that transfer from
Haemophilus
to
Neisseria has occurred. (i) The neisserial island of
similarity
to
H. influenzae includes a sequence upstream of
orf2 which matches
hUS consensus at all 29 bp and a paired
(+/

) uptake sequence
downstream of
orf2 which matches
consensus at all but 1 bp. These
hUSs align with those in Int1. (ii)
The
N. gonorrhoeae island
includes 116 bp from
H. influenzae holB and 225 bp from
H. influenzae tmk in
the intergenic spaces flanking
orf2. The
N. meningitidis island also includes 116 bp from
H. influenzae
holB, but it has
only 77 bp from
H. influenzae tmk. The
77-bp
tmk sequence aligns
with 225 bp of the Rd sequence but
with a 148-bp internal deletion.
The presence in the meningococcal
genome of nonfunctional bits
of adjacent
H. influenzae genes
is highly suggestive of the direction
of transfer. (iii) The G+C
content of
N. meningitidis serotype
B
orf2 is
40.0%, and that of
N. gonorrhoeae is 39.7%. These values
are much closer to the
H. influenzae genomic average of
38.2%
than to the neisserial average of 51.8%.
Recent transfer of lav from Haemophilus to
Neisseria.
Genetic similarity suggests that the
Haemophilus lav gene and the Neisseria orf2 gene
diverged relatively recently, not long before orf2 genes
diverged in Neisseria. As a more reliable indicator of
genetic distance, we compared noncoding, intergenic DNA sequence, which
was expected to be subject to fewer selective constraints than coding
sequence. The 256-bp sequence immediately upstream of lav
can be aligned for all four bacterial genomes (Fig. 2a). It includes 19 codons of holB, but the remaining 206 bp in
H. influenzae (and all 256 bp in Neisseria) are
noncoding. The H. influenzae sequence is 94.5% identical to
that of N. gonorrhoeae and 93.8% identical to that of
N. meningitidis, whereas sequences from the two neisserial
species are 96.9% identical. The distance between the sequences for
H. influenzae and Neisseria is approximately twice that between the sequences for neisserial species (Fig. 2b).
Similar relative distances (but higher absolute divergences) were found
for the lav region including the
-barrel plus the linker
(Fig. 2b). In comparison, the DNA sequence of an evolutionarily conserved regulatory gene, rpoD (encoding the major sigma
factor), is only 47% identical between H. influenzae and
Neisseria, yet rpoD similarity between the two
Neisseria sequences is 97% (Fig. 2b).

View larger version (23K):
[in this window]
[in a new window]
|
FIG. 2.
(a) Alignment (MEGALIGN) of DNA sequence 256 bp upstream
of lav homologs. Strain designations are as described for
Fig. 1, except that NM-B is MC58 and NM-A is Z2491. The initiating
codon (complementary strand) of holB is indicated (left
arrow), and the sequence includes 19 codons at the N terminus-encoding
portion of of holB. The start of lav is also
indicated (right arrow). An 18-bp sequence from an hUS is underlined;
this sequence has been duplicated from a paired hUS at the
tmk end of the Int1 lav island. (b) Phylograms of
aligned sequences from the 256-bp upstream sequence (5' to
lav), from the -barrel region (last 385 codons, including
the predicted linker region), and from rpoD (17, 36,
48). Phylogenies were derived with the ClustalV program of
MEGALIGN, using a combination of the unweighted pair group method using
arithmetic averages algorithm and the neighbor-joining method
(39, 46). (c) H. influenzae lav is anomalously
similar to its N. meningitidis homologs relative to pairwise
comparisons between housekeeping genes shared by H. influenzae and N. meningitidis. Linear regression of
inferred amino acid sequence identity on DNA identity of 12 housekeeping genes (open squares) in N. meningitidis
serotype B and H. influenzae Rd is shown. Filled circles
(arrows) indicate Int1 lav compared to homologs in N. meningitidis serotype B and N. meningitidis serotype A,
in ascending order of DNA similarity. Genes compared, in ascending
order of DNA similarity, were tpiA, lig (lig-1 in
N. meningitidis serotype B) dnaG, mutS, gyrA, rpoD,
cysK, aroG, rpoB, recA, sucD, and ilvD
(ilvD-1 in N. meningitidis serotype B). The
sequences were compared using ClustalV (MEGALIGN).
|
|
The atypically high similarity of
lav homologs, relative to
that of other genes shared by
Haemophilus and
Neisseria, is evident
in Fig.
2c, which shows a comparison
of 12 pairs of evolutionarily
conserved housekeeping genes shared by
these species to provide
a baseline. The mean DNA similarity between
Haemophilus and
Neisseria was 47.3% (standard
deviation, 8.7%; range, 35.0 to 60.2%). The
values for the
lav homologs lie above the range of DNA similarity
represented for these highly conserved genes. With respect to
DNA
similarity, the Lav amino acid sequence is more divergent
than that of
or typical housekeeping genes, presumably because
lav is
exposed to a different type of selection (Fig.
2c).
Chimeric origin of lav as evidence for domain
shuffling.
A scan of G+C content across the lsi-dld
region revealed additional detail (Fig. 3). The island
of similarity between H. influenzae and neisserial
chromosomes coincides with a dip in G+C content to a level approaching
the H. influenzae genomic average; the profile is similar in
detail for all three species, again indicating a recent common origin.
G+C content is not uniform across lav but rather rises
sharply at 850 ± 50 bp after the GCAA repeat region. Thus, Int1
lav is composed of two distinct segments, a 5' segment of
283 codons averaging 35.1% G+C and a 3' segment of 376 codons
averaging 40.8% G+C.

View larger version (31K):
[in this window]
[in a new window]
|
FIG. 3.
G+C content of regions containing lav
homologs, scanned with a sliding window of 100 bp, in 25-bp increments.
hi, H. influenzae Int1; nm, N. meningitidis
serotype A Z2491; ng, N. gonorrhoeae. The interval on the
abscissa corresponds to bp 113 to 5147 of N. gonorrhoeae
contig 140 aligned with N. meningitidis serotype A interval
1660274 to 1654800 (which has a deletion within tmk) by
creating a 148-bp gap in the N. meningitidis serotype A
sequence starting at bp 4033 and the entire Int1 sequence between
tmk and holB.
|
|
The two segments defined by G+C content approximately coincide with
predicted boundaries of the two major structural and functional
domains
of an AT. The closest homolog to Lav with functionally
characterized
domains is the 1,286-aa
E. coli AIDA-I preprotein
(
31). Amino acid sequence homology begins near the
processing
site at AIDA-I Leu
840 and includes the linker
region and the 14-strand

-barrel, which
begins at AIDA-I
Ala
1002. On the basis of alignment with AIDA-I, the

-domain of Int1
lav is predicted to extend from codon 385 to 659 codon and the
linker is predicted to extend from codon ~290 to
codon 384. The
low-G+C segment of
lav (first 283 codons)
thus corresponds roughly
to the predicted passenger domain, and the
high-G+C segment (last
376 codons) corresponds to the

-barrel domain
plus
linker.
The two segments of
lav have evolved at markedly different
rates. The role of selection in fixing mutations can be assessed
by
comparing silent-substitution rates (
Ks = synonymous changes
per 100 synonymous sites) to nonsilent-substitution
rates (
Ka = nonsynonymous changes per 100 nonsynonymous sites) (
32); the
Ks/Ka ratio is usually >10 for
conserved housekeeping genes. Silent-substitution
rates in the C
terminus-encoding segment are consistent with recent
common origin of
these genes (Table
2). The
Ks
values derived
from a comparison of Int1 and
Neisseria are
only 1.2- to 1.5-fold
higher than those derived from a comparison of
N. meningitidis and
N. gonorrhoeae. However,
nonsynonymous-substitution rates
differ substantially between the
segments encoding the N and C
termini within each lineage. The
Ks/Ka ratio is

4 for the segment
encoding the C terminus, but it is

1 for the segment encoding
the N
terminus-segment (Table
2). These differences suggest that
the two
major segments of the gene are subject to different evolutionary
constraints.
Distribution of lav.
Various H. influenzae strains and clinical isolates, including NTHi and
representatives from each capsular serotype, were examined by PCR and
genomic Southern analyses for the presence of an insert at the
holB-tmk junction. The junction was amplified in 26 strains,
including 8 from encapsulated lineages (Rd, Hib [Eagan], and
reference type strains of serotypes a to f), 9 respiratory NTHi
isolates, and 9 NTHi isolates from patients with invasive disease
(Table 1). A 2.4-kb fragment containing a lav-related sequence was found at the holB-tmk site in six strains,
representing just two clades of NTHi. Southern analysis confirmed the
absence of a lav sequence in the other strains or at other
chromosomal sites. Strains having lav-related islands
included Int1 and two respiratory isolates previously identified as
clonally related to Int1 (T. Mutangadura and M. Golomb, unpublished
data), 1128 (isolated from a patient with otitis media)
(13) and AAr160 (isolated from a tracheal aspirate)
(12). Partial sequencing revealed lav genes at
identical sites in both isolates. The sequence of the AAr160
lav island was identical to that of Int1 for 600 bp,
starting at position 1 in Fig. 2a, with the exception of the number of
GCAA repeats and a G
T polymorphism at Int1 position 382, which
conservatively changes an L to an F within the lav-encoded signal sequence. A 2.4-kb island was found at the same chromosomal site
in R1967, the reference biotype aegyptius strain, and in R2140 and
R2141, clonally related biotype aegyptius strains isolated from
patients with BPF (Table 3).
Domain shuffling and the emergence of a novel biotype
aegyptius AT gene.
The holB-tmk island of R1967
(biotype aegyptius) had endpoints identical to those in Int1. The
island contained a single potential ORF of 2,102 bp, at a position
aligned with Int1 lav. Downstream from the predicted
initiating ATG were 25 GCAA repeats, placing the coding sequence out of
reading frame. Gain of a single repeat would allow translation of a
702-aa polypeptide with a predicted 61-aa signal polypeptide. The
C-terminal
-domain of this polypeptide is strongly homologous to
Lav. We refer to the putative biotype aegyptius gene at this site as
las. Partial sequencing of the R2140 and R2141 islands (BPF)
revealed las homologs (Table 3).
Like
lav, biotype aegyptius
las is a mosaic of
low- and high-G+C segments (Fig.
4a). Although the
region immediately upstream
from biotype aegyptius
las and
the
las 
-domain are highly homologous
(~90%) to
corresponding regions in
lav (Fig.
2), the N
terminus-encoding
portion of the gene, following the GCAA repeats, has
no significant
homology to
lav or to any sequence in the
databanks at either
the DNA or amino acid sequence level (Fig.
4b).
Dispersed motifs
within the inferred N-terminal region of Las suggest a
distant
common ancestry with
H. influenzae Lav; thus, both
proteins share
the sequence SLWEPR(W/F)NS at corresponding sites (aa
222 to 230
in Int1 and aa. 226 to 234 in biotype aegyptius). The
junction
with the high-homology segment (Fig.
4a, a) is near the G+C
transition
point, which coincides with that in
lav (~850
bp or 283 codons).
High amino acid sequence similarity starts at the
position corresponding
to
lav codon 277 (biotype aegyptius
codon 280). Evidently the
lav and
las genes
diverged from an ancestral
H. influenzae sequence
by
recombinational fusion of a novel passenger domain to the
H. influenzae 
-domain near the predicted passenger-linker
boundary.

View larger version (42K):
[in this window]
[in a new window]
|
FIG. 4.
(a) G+C bias of biotype aegyptius las codons,
computed using CodonPreference with a 25-codon window. The abscissa
shows nucleotide position in las following the GCAA repeat
region. The horizontal line indicates average H. influenzae
codon usage. (b) Homology between H. influenzae lav and
H. influenzae las, using a 10-bp window, as a function of
nucleotide position.
|
|
 |
DISCUSSION |
Lateral acquisition of lav by NTHi.
The restricted
distribution of lav could be explained either by transfer
from another species or by inheritance of the gene from an H. influenzae ancestor followed by loss within multiple lineages.
Genomic comparisons favor lateral transfer. In Rd, holB and
tmk overlap by 4 bp, with the terminating TGA codon of
tmk overlapping the initiating ATG codon of holB.
A diverse set of bacterial species have similar arrangements (Fig.
5A). Discussion In Pasteurella, sister genus
to Haemophilus (38), the genes are contiguous
but not overlapping
(http://www.cbc.umn.edu/ResearchProjects/AGAC/Pm/pmhome.html). The tmk and holB genes are also contiguous
in Buchnera, a gram-negative endocellular parasite of aphids
that has a subset of the genes found in E. coli.
Buchnera is believed to have diverged before the common
ancestor of H. influenzae and E. coli appeared
(41). The contiguity of holB and tmk
in other gram-negative bacteria and even in Mycoplasma
(phylogenetically closer to gram-positive bacteria) suggests that it is
the ancestral arrangement for eubacteria. Thus, lav was most
likely inserted into the chromosome of an Int1 ancestor.

View larger version (14K):
[in this window]
[in a new window]
|
FIG. 5.
(A) Arrangement of tmk and holB
in various genomes. Nucleotides at the junction of the two genes are
shown. In H. influenzae Rd and E. coli, the
reading frame of holB is offset 2 relative to that of
tmk; in all others shown, the relative reading frame of
holB is offset 1. In P. multocida
(http://www.cbc.umn.edu/ResearchProjects/AGAC/Pm/pmhome.html), the two
genes are separated by 4 bp. GenBank accession numbers and sources are
as follows: Rd, reference 17 E. coli, reference
5; P. multocida, contig 186, nucleotides 22421 to 23201; Yersinia pestis, accession no. AF065312
(11); Mycoplasma pneumoniae, accession no.
AE000016 (22); and Mycoplasma fermentans,
accession no. AF100324 (8). (B) Proposed scenario for
origin of lav homologs in various organisms. The original
form of the gene in H. influenzae is arbitrarily shown as
Int1 lav but could have been either a lav or a
las ancestor. nm, N. meningitidis; ng, N. gonorrhoeae; hi, H. influenzae; Hae, biotype
aegyptius.
|
|
The present distribution of
lav requires at least three
interspecies transfers (Fig.
5B). The extrageneric origin of
lav is
supported by its limited distribution within
Haemophilus, by G+C
content analysis, and by analysis of
junctions with adjoining
genes. Although the
lav segment
does not differ in G+C content
from the
H. influenzae
average, it is actually composed of two
segments, one atypically low
and one atypically high in G+C content.
This suggests that the gene
originated in another, unidentified,
organism by fusion between
segments from different species. The
chimeric gene was then laterally
acquired by
H. influenzae-and
inserted within a stem-loop
hUS pair originally present in single
copy in a chromosome resembling
that of Rd. The stem-loop pair
was partially duplicated during this
integration, resulting in
direct repeats flanking the
lav
island in an ancestor of Int1.
Insertion of virulence-associated genes
within hUS pairs has previously
been noted for the
tna
cluster (
30); the targeting of regions
of high potential
secondary structure suggests that paired hUSs,
like tRNA genes, may
serve as targets for phage integrases in
Haemophilus.
Subsequently, the
lav island with adjoining bits of
holB and
tmk was incorporated by
Neisseria, most likely by transformation.
Perhaps the
fragment recombined within a region of weak, fortuitous
homology,
though none has been identified. (The two hUS pairs,
located internally
within the transferred fragment, presumably
played no role in
neisserial uptake of
Haemophilus DNA. Their
presence in
Neisseria merely reflects the recent origin of this
fragment
in
Haemophilus.) The strongest evidence indicating polarity
of transfer is the presence of pseudogene fragments of the adjoining
Haemophilus holB and
tmk genes. (Functional,
endogenous
holB and
tmk genes are located at
different chromosomal sites in
Neisseria.)
Notably, the
usual criteria for direction of transfer (G+C content,
within-species
distribution, and presence of signature sequences
like the hUS) are
less reliable and may be deceptive at times.
For instance, the G+C
content of the Lav-encoding part of the
H. influenzae island
is species typical not because of its origin
within
H. influenzae but because it is an artifact of averaging
across two
dissimilar segments. In general, detailed sequence
analysis is needed
to infer direction of lateral
transfer.
The appearance of
lav in both
N. meningitidis and
N. gonorrhoeae could be explained either by preexistence in
a common ancestor
or by acquisition by one lineage followed by
homologous transformation
into the other (Fig.
5B). A third
possibility, independent acquisition
by
N. meningitidis and
N. gonorrhoeae from two different
H. influenzae lineages, is much less likely because the junction of neisserial
sequence with the
H. influenzae island is identical, to the
base
pair, in
N. gonorrhoeae and
N. meningitidis.
The presence of
lav (
orf2) in all
N. meningitidis strains surveyed to date (
37)
is
consistent with its preexistence in an
N. meningitidis
ancestor,
though it could also be explained by transmission across
N. meningitidis lineages and strong selective
advantage.
DNA sequence similarity attests that the meningococcus and the
gonococcus have diverged relatively recently, although no time
frame
estimates have been published. Most genes are shared by
both pathogens
and are on average

98% identical (
47,
49).
If emergence
of
N. gonorrhoeae as a sexually transmitted pathogen
required relatively high human population densities, divergence
may
have occurred as recently as the Neolithic period, 10,000
years ago.
Supposing that
lav preexisted in the common ancestor
of
N. gonorrhoeae and
N. meningitidis, it must have
emerged shortly
before species divergence, since it is only slightly
less similar
between
H. influenzae and
Neisseria
than between neisserial species.
Once in
Neisseria, lav
evolved rapidly, especially within its
passenger domain. Nonsilent- and
silent-mutation rates are approximately
equal within the passenger
domains, suggesting either selective
neutrality or (more likely)
diversifying selection acting on variable
regions interspersed within a
more conserved framework, as suggested
by the data presented in Fig.
1C.
In
Haemophilus, phase variation by repeat slippage is
associated with surface molecules exposed to immune surveillance
(
23,
51), and slippage rates increase with the number of
repeats
(
20). The 19 GCAA repeats of Int1
lav
and the 25 GCAA repeats
of biotype aegyptius
las should be
relatively unstable in vivo,
consistent with an actively expressed
gene. We have detected phase
variants ranging from 18 to 22 repeats in
laboratory populations
(data not shown). Although the Int1
lav variant described here
is out of reading frame, in-frame
variants would arise readily
in natural populations and are represented
by two clinical isolates
(Table
3). Compared to
Haemophilus
orf2 genes, the
N. meningitidis orf2 genes have fewer
GCAA repeats; among 41
N. meningitidis isolates
possessing
orf2, the number of repeats ranged from 1 to 12 (
37).
A mosaic origin for
lav was inferred from a G+C content
transition at the boundary of its presumed passenger domain with the
linker and

-barrel domains. Similarly, the junction of nonhomology
between
lav and
las coincides with the G+C
transition and inferred
domain boundaries of both genes. On the basis
of quite different
evidence (discordance between phylogenies based on
individual
domains), Loveless and Saier have proposed that AT proteins
evolve
by domain shuffling (
28). A functionally novel AT
can arise
by linking a new passenger activity to a generic

-barrel
pore.
Our analysis provides independent evidence for the combinatorial
origin and subsequent reshuffling of at least one AT
protein.
Within the small sample of
H. influenzae strains examined,
lav was restricted to a few NTHi strains with unusual
virulence
potential. Although Int1 lacks the type b polysaccharide
capsule
which protects against lysis by human serum, it is more serum
resistant than most NTHi strains (B. Williams and A. L. Smith,
unpublished data). Int1, AAr160, and 1128 constitute a clade by
various
criteria, including the results of pulsed-field gel electrophoresis
after restriction with rare-cutting enzymes (Mutangadura and Golomb,
unpublished). All three strains have the phage HP2, otherwise
present
in only a small minority of NTHi strains (Williams and
Smith,
unpublished). The biotype aegyptius is associated with
unusual
virulence because it includes BPF
isolates.
As biotype aegyptius strains and Int1 belong to different phylogenetic
subgroups, it is unlikely that they inherited
lav from
a
common ancestor. Rather, it is likely that the first
H. influenzae clade to acquire the gene passed it to one or more
other clades
by transformation and homologous recombination within
flanking
DNA. Once a laterally transferred fragment has been acquired
by
a population of naturally transformable bacteria, it can readily
be
assimilated into the species by co-opting linked homologous
sequence
and uptake signals. Interstrain and interspecies transfer
implies a
shared selective advantage in certain host
environments.
This is the second report of gene transfer from
Haemophilus
to
Neisseria (
25) and the first of a gene
belonging to a family
widely associated with virulence. Gene flow from
Neisseria to
Haemophilus has not yet been
documented and may be rarer than
gene flow in the opposite direction.
Natural transformation is
more frequent and promiscuous in
Neisseria than in
Haemophilus,
and the
requirement for the nUS can be bypassed experimentally
(
6). In contrast, the requirement of
Haemophilus transformation
for the hUS is stringent
(
42). Many genes in the
N. meningitidis genome
appear to be of external origin, as judged by anomalous
G+C content; in
contrast, the
H. influenzae Rd genome contains
fewer islands
of unusual G+C content (
16,
17,
36). One way
to discover a
sequence of external origin is to search a genome
for heterologous
uptake sequences. When Kroll et al. (
25) searched
the Rd
genome with the 10-bp nUS, its only occurrence was within
an ORF
believed to have been transferred from
Haemophilus to
Neisseria.
When the small Int1-specific library was searched
with the nUS,
a single hit was found. It consisted of two oppositely
oriented
nUS sequences in the

/+ configuration, separated by 71 bp,
within
a slightly longer region of dyad symmetry:
T
TTCAGACGGCATN
67AT
GCCGTCTGAAA
(inverted repeat nUSs are underlined and N signifies any base).
The probability that two nUSs will occur by chance within 100
bp of
each other is <10
4. Inverse PCR identified a >2.5-kb
island without significant
homology to Rd (data not shown). The two
nUSs are within a 1,431-bp
ORF that is 65% identical in deduced amino
acid sequence to the
Vibrio cholerae ORF VC1769 product, a
presumed methyltransferase
(HsdM) subunit of a type I restriction
endonuclease system (
20).
G+C content of this Int1 ORF is
46.1%, a value midway between
H. influenzae and neisserial
averages. The
Vibrio ORF is 43% G+C
and part of an island
of anomalously low G+C content. Despite
possession of a neisserial
signature, this ORF did not match sequence
entries for
N. meningitidis serotype A,
N. meningitidis serotype
B, or
N. gonorrhoeae. We conclude that the Int1 nUS-bearing
ORF,
although laterally derived, is unlikely to be from
Neisseria.
Humans are the exclusive natural host of
H. influenzae and
the meningococcus. Virulence genes are likely to be recent
specializations
for competitive adaptation to a human host, and their
most obvious
sources are other human pathogens and commensals.
Repeat-regulated
genes such as
lav are well equipped for
horizontal mobility, since
their regulatory system is self-contained.
In contrast to selfish
operons, which must be transferred as a unit
(
26), contingency
genes are controlled without the need
for accessory proteins or
for extensive compatibility between
regulatory elements from different
species. The emergence of invasive
NTHi strains, such as those
causing BPF and meningitis in
immunocompetent children, suggests
that contingency gene transfer
between dangerous respiratory pathogens
may be an ongoing
process.
 |
ACKNOWLEDGMENTS |
We are grateful to Robert Munson, Janet Gilsdorf, and
Loek van Alphen for provision of bacterial strains, to Michael Calcutt for helpful discussions, and to Stephen Lory for advice and help in
constructing the M13 library. We thank Arnie Kas, University of
Washington, for an unpublished G+C scan program. Incompletely sequenced
microbial genomes were accessed on the World Wide Web from the
Gonococcal Sequencing Project at Oklahoma University and the
Pasteurella multocida Sequencing Project at the University of Minnesota.
This work was supported by NIH grant 5RO1 HG01475 to Maynard
Olson, University of Washington, by NIH grant AI 44002 to A.L.S., and
by University of Missouri Research Board grants to M.G. and A.L.S.
 |
FOOTNOTES |
*
Corresponding author. Mailing address: 110 Tucker Hall,
Division of Biological Sciences, University of Missouri, Columbia, MO
65211. Phone: (573) 882-9628. Fax: (573) 882-0123. E-mail: GolombM{at}missouri.edu.
 |
REFERENCES |
| 1.
|
Altschul, S. F.
1993.
A protein alignment scoring system sensitive at all evolutionary distances.
J. Mol. Evol.
36:290-300[CrossRef][Medline].
|
| 2.
|
Altschul, S. F.,
W. Gish,
W. Miller,
E. W. Myers, and D. J. Lipman.
1990.
Basic local alignment search tool.
J. Mol. Biol.
215:403-410[CrossRef][Medline].
|
| 3.
|
Barbour, M. L.
1996.
Conjugate vaccines and the carriage of Haemophilus influenzae type b: Emerg.
Infect. Dis.
2:176-182.
|
| 4.
|
Benz, I., and M. A. Schmidt.
1992.
AIDA-I, the adhesin involved in diffuse adherence of the diarrhoeagenic Escherichia coli strain 2787 (O126:H27), is synthesized via a precursor molecule.
Mol. Microbiol.
6:1539-1546[CrossRef][Medline].
|
| 5.
|
Blattner, F. R.,
G. Plunkett III,
C. A. Bloch,
N. T. Perna,
V. Burland,
M. Riley,
J. Collado-Vides,
J. D. Glasner,
C. K. Rode,
G. Mayhew,
J. Gregor,
N. W. Davis,
H. Kirkpatrick,
M. Goeden,
D. Rose,
R. Mau, and Y. Shao.
1997.
The complete genome sequence of Escherichia coli K-12.
Science
277:1453-1474[Abstract/Free Full Text].
|
| 6.
|
Boyle-Vavra, S., and H. S. Seifert.
1996.
Uptake-sequence-independent DNA transformation exists in Neisseria gonorrhoeae.
Microbiology
142:2839-2845[Abstract/Free Full Text].
|
| 7.
|
Brenner, D. J.,
L. W. Mayer,
G. M. Carlone,
L. H. Harrison,
W. F. Bibb,
M. C. de Cunto Brandileone,
F. O. Sottnek,
K. Iriono,
M. W. Reeves,
S. J. M. Swenson,
K. A. Birkness,
R. S. Weyant,
S. F. Berkley,
T. C. Woods,
A. G. Steigerwalt,
P. A. D. Grimont,
R. M. McKinney,
D. W. Fleming,
L. L. Gheesling,
R. C. Cooksey,
R. J. Arko,
C. V. Broome, and The Brazilian Purpuric Fever Study Group.
1988.
Biochemical, genetic, and epidemiologic characterization of Haemophilus influenzae biogroup aegyptius (Haemophilus aegyptius) strains associated with Brazilian purpuric fever.
J. Clin. Microbiol.
26:1524-1534[Abstract/Free Full Text].
|
| 8.
|
Calcutt, M. J.,
M. F. Kim,
A. B. Karpas,
P. F. Mühlradt, and K. S. Wise.
1999.
Differential posttranslational processing confers intraspecies variation of a major surface lipoprotein and a macrophage-activating lipopeptide of Mycoplasma fermentans.
Infect. Immun.
67:760-771[Abstract/Free Full Text].
|
| 9.
|
Catlin, B. W.,
J. W. Bendler, and S. H. Goodgal.
1972.
The type b capsulation locus of Haemophilus influenzae: map location and size.
J. Gen. Microbiol.
70:411-422[Abstract/Free Full Text].
|
| 10.
|
Cheetham, B. F., and M. E. Katz.
1995.
A role for bacteriophages in the evolution and transfer of bacterial virulence determinants.
Mol. Microbiol.
18:201-208[CrossRef][Medline].
|
| 11.
|
Chenal-Francisque, V.,
L. Tourneux,
E. Carniel,
P. Christova,
I. L. de la Sierra,
O. Barzu, and A. M. Gilles.
1999.
The highly similar TMP kinases of Yersinia pestis and Escherichia coli differ markedly in their AZTMP phosphorylating activity.
Eur. J. Biochem.
265:112-119[Medline].
|
| 12.
|
Clemans, D. L.,
C. F. Marrs,
M. Patel,
M. Duncan, and J. R. Gilsdorf.
1998.
Comparative analysis of Haemophilus influenzae hifA (pilin) genes.
Infect. Immun.
66:656-663[Abstract/Free Full Text].
|
| 13.
|
Coleman, T.,
S. Grass, and R. Munson, Jr.
1991.
Molecular cloning, expression, and sequence of the pilin gene from nontypeable Haemophilus influenzae M37.
Infect. Immun.
59:1716-1722[Abstract/Free Full Text].
|
| 14.
|
Elkins, C.,
C. E. Thomas,
H. S. Seifert, and P. F. Sparling.
1991.
Species-specific uptake of DNA by gonococci is mediated by a 10-base-pair sequence.
J. Bacteriol.
173:3911-3913[Abstract/Free Full Text].
|
| 15.
|
Everest, P.,
J. Li,
G. Douce,
I. Charles,
J. De Azavedo,
S. Chatfield,
G. Dougan, and M. Roberts.
1996.
Role of the Bordetella pertussis P.69/pertactin protein and the P.69/pertactin RGD motif in the adherence to and invasion of mammalian cells.
Microbiology
142:3261-3268[Abstract/Free Full Text].
|
| 16.
|
Faguy, D. M.
2000.
The controlled chaos of shifty pathogens.
Curr. Biol.
10:R498-R501[CrossRef][Medline].
|
| 17.
|
Fleischmann, R. D.,
M. D. Adams,
O. White,
R. A. Clayton,
E. F. Kirkness,
A. R. Kerlavage,
C. J. Bult,
J.-F. Tomb,
B. A. Dougherty,
J. M. Merrick,
K. McKenney,
G. Sutton,
W. FitzHugh,
C. Fields,
J. D. Gocayne,
J. Scott,
R. Shirley,
L.-I. Liu,
A. Glodek,
J. M. Kelley,
J. F. Weidmann,
C. A. Phillips,
T. Spriggs,
E. Hedblom,
M. D. Cotton,
T. R. Utterback,
M. C. Hanna,
L. D. Nguyen,
D. M. Saudek,
R. C. Brandon,
L. D. Fine,
J. L. Fritchman,
N. S. M. Geoghagen,
C. L. Gnehm,
L. A. McDonald,
K. V. Small,
C. M. Fraser,
H. O. Smith, and J. C. Venter.
1995.
Whole-genome random sequencing and assembly of Haemophilus influenzae Rd.
Science
269:496-512[Abstract/Free Full Text].
|
| 18.
|
Foxwell, A. R.,
J. M. Kyd, and A. W. Cripps.
1998.
Nontypeable Haemophilus influenzae: pathogenesis and prevention.
Microbiol. Mol. Biol. Rev.
62:294-308[Abstract/Free Full Text].
|
| 19.
|
Geluk, F.,
P. P. Eijk,
S. M. van Ham,
H. M. Jansen, and L. van Alphen.
1998.
The fimbria gene cluster of nonencapsulated Haemophilus influenzae.
Infect. Immun.
66:406-417[Abstract/Free Full Text].
|
| 20.
|
Heidelberg, J. F.,
J. A. Eisen,
W. C. Nelson,
R. A. Clayton,
M. L. Gwinn,
R. J. Dodson,
D. H. Haft,
E. K. Hickey,
J. D. Peterson,
L. Umayam,
S. R. Gill,
K. E. Nelson,
T. D. Read,
H. Tettelin,
D. Richardson,
M. D. Ermolaeva,
J. Vamathevan,
S. Bass,
H. Qin,
I. Dragoi,
P. Sellers,
L. McDonald,
T. Utterback,
R. D. Fleishmann,
W. C. Nierman, and O. White.
2000.
DNA sequence of both chromosomes of the cholera pathogen Vibrio cholerae.
Nature
406:477-483[CrossRef][Medline].
|
| 21.
|
Henderson, I. R.,
F. Navarro-Garcia, and J. P. Nataro.
1998.
The great escape: structure and function of the autotransporter proteins.
Trends Microbiol.
6:370-378[CrossRef][Medline].
|
| 22.
|
Himmelreich, R.,
H. Hilbert,
H. Plagens,
E. Pirkl,
B. C. Li, and R. Herrmann.
1996.
Complete sequence analysis of the genome of the bacterium Mycoplasma pneumoniae.
Nucleic Acids Res.
24:4420-4449[Abstract/Free Full Text].
|
| 23.
|
Hood, D. W.,
M. E. Deadman,
M. P. Jennings,
M. Bisceric,
R. D. Fleischmann,
J. C. Venter, and E. R. Moxon.
1996.
DNA repeats identify novel virulence genes in Haemophilus influenzae.
Proc. Natl. Acad. Sci. USA
93:1121-1125.
|
| 24.
|
Jennings, M. P.,
M. Bisercic,
K. L. R. Dunn,
M. Virji,
A. Martin,
K. E. Wilks,
J. C. Richards, and E. R. Moxon.
1995.
Cloning and molecular analysis of the lsi1 (rfaF) gene of Neisseria meningitidis which encodes a heptosyl-2-transferase involved in LPS biosynthesis: evaluation of surface exposed carbohydrates in LPS mediated toxicity for human endothelial cells.
Microb. Pathog.
19:391-407[CrossRef][Medline].
|
| 25.
|
Kroll, J. S.,
K. E. Wilks,
J. L. Farrant, and P. R. Langford.
1998.
Natural genetic exchange between Haemophilus and Neisseria: intergeneric transfer of chromosomal genes between major human pathogens.
Proc. Natl. Acad. Sci. USA
95:12381-12385[Abstract/Free Full Text].
|
| 26.
|
Lawrence, J. G., and J. R. Roth.
1996.
Selfish operons: horizontal transfer may drive the evolution of gene clusters.
Genetics
143:1843-1860[Abstract].
|
| 27.
|
Lett, M.-C.,
C. Sasakawa,
N. Okada,
T. Sakai,
S. Makino,
M. Yamada,
K. Komatsu, and M. Yoshikawa.
1989.
virG, a plasmid-coded virulence gene of Shigella flexneri: identification of the virG protein and determination of the complete coding sequence.
J. Bacteriol.
171:353-359[Abstract/Free Full Text].
|
| 28.
|
Loveless, B. J., and M. H. Saier.
1997.
A novel family of channel-forming, autotransporting, bacterial virulence factors.
Mol. Membr. Biol.
14:113-123[Medline].
|
| 29.
|
Maiden, M. C. J., and I. M. Feavers.
1995.
Population genetics and global epidemiology of the human pathogen Neisseria meningitidis, p. 269-293.
In
S. Baumberg, J. P. W. Young, E. M. H. Wellington, and J. R. Saunders (ed.), Population genetics of bacteria. Cambridge University Press, Cambridge, United Kingdom.
|
| 30.
|
Martin, K.,
G. Morlin,
A. Smith,
A. Nordyke,
A. Eisenstark, and M. Golomb.
1998.
The tryptophanase gene cluster of Haemophilus influenzae type b: evidence for horizontal gene transfer.
J. Bacteriol.
180:107-118[Abstract/Free Full Text].
|
| 31.
|
Mhlanga-Mutangadura, T.,
G. Morlin,
A. L. Smith,
A. Eisenstark, and M. Golomb.
1998.
Evolution of the major pilus gene cluster of Haemophilus influenzae.
J. Bacteriol.
180:4693-4703[Abstract/Free Full Text].
|
| 32.
|
Moxon, E. R.,
P. B. Rainey,
M. A. Nowak, and R. E. Lenski.
1994.
Adaptive evolution of highly mutable loci in pathogenic bacteria.
Curr. Biol.
4:24-33[CrossRef][Medline].
|
| 33.
|
Mühlemann, K.,
M. Balz,
S. Aebi, and K. Schopfer.
1996.
Molecular characteristics of Haemophilus influenzae causing invasive disease during the period of vaccination in Switzerland: analysis of strains isolated between 1986 and 1993.
J. Clin. Microbiol.
34:560-563[Abstract].
|
| 34.
|
Nizet, V.,
K. F. Colina,
J. R. Almquist,
C. E. Rubens, and A. L. Smith.
1996.
A virulent nonencapsulated Haemophilus influenzae.
J. Infect. Dis.
173:180-186[Medline].
|
| 35.
|
Ochman, H.
1997.
How Salmonella became a pathogen.
Trends Microbiol.
9:343-349.
|
| 36.
|
Parkhill, J.,
M. Achtman,
K. D. James,
S. D. Bentley,
C. Churcher,
S. R. Klee,
G. Morelli,
D. Basham,
D. Brown,
T. Chillingworth,
R. M. Davies,
P. Davis,
K. Devlin,
T. Feltwell,
N. Hamlin,
S. Holroyd,
K. Jagels,
S. Leather,
S. Moule,
K. Mungall,
M. A. Quail,
M. A. Rajandream,
K. M. Rutherford,
M. Simmonds,
J. Skelton,
S. Whitehead,
B. G. Spratt, and B. G. Barrell.
2000.
Complete DNA sequence of a serogroup A strain of Neisseria meningitidis Z2491.
Nature
404:502-506[CrossRef][Medline].
|
| 37.
|
Peak, I. R.,
M. P. Jennings,
D. W. Hood, and E. R. Moxon.
1999.
Tetranucleotide repeats identify novel virulence determinant homologues in Neisseria meningitidis.
Microb. Pathog.
26:13-23[CrossRef][Medline].
|
| 38.
|
Pohl, S.
1981.
DNA relatedness in HPA, p. 245-253.
In
M. Kilian, W. Frederiksen, and E. L. Biberstein (ed.), Haemophilus, Actinobacillus, and Pasteurella. Academic Press, New York, N.Y.
|
| 39.
|
Saitou, N., and M. Nei.
1987.
The neighbor-joining method: a new method for reconstructing phylogenetic trees.
Mol. Biol. Evol.
4:406-425[Abstract].
|
| 40.
|
Saunders, N. J.,
A. C. Jeffries,
J. F. Peden,
D. W. Hood,
H. Tettelin,
R. Rappuoli, and E. R. Moxon.
2000.
Repeat-associated phase variable genes in the complete genome sequence of Neisseria meningitidis strain MC58.
Mol. Microbiol.
37:207-215[CrossRef][Medline].
|
| 41.
|
Shigenobu, S.,
H. Watanabe,
M. Hattori,
Y. Sakaki, and H. Ishikawa.
2000.
Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS
Nature
407:81-86[CrossRef][Medline].
|
| 42.
|
Sisco, K. L., and H. O. Smith.
1979.
Sequence-specific DNA uptake in Haemophilus transformation.
Proc. Natl. Acad. Sci. USA
76:972-976[Abstract/Free Full Text].
|
| 43.
|
Smith, A. L.,
D. H. Smith,
D. R. Averill, Jr.,
J. Marino, and E. R. Moxon.
1973.
Production of Haemophilus influenzae b meningitis in infant rats by intraperitoneal inoculation.
Infect. Immun.
8:278-290[Abstract/Free Full Text].
|
| 44.
|
Smith, H. O.,
J.-F. Tomb,
B. A. Dougherty,
R. D. Fleischmann, and J. C. Venter.
1995.
Frequency and distribution of DNA uptake signal sequences in the Haemophilus influenzae Rd genome.
Science
269:538-540[Abstract/Free Full Text].
|
| 45.
|
Smith, J. M.,
N. H. Smith,
M. O'Rourke, and B. G. Spratt.
1993.
How clonal are bacteria?
Proc. Natl. Acad. Sci. USA
90:4384-4388[Abstract/Free Full Text].
|
| 46.
|
Sneath, P. H. A., and R. R. Sokal.
1973.
Numerical taxonomy. W. H.
Freeman, San Francisco, Calif.
|
| 47.
|
Spratt, B. G.,
N. H. Smith,
J. Zhous,
M. O'Rourke, and E. Feil.
1995.
The population genetics of the pathogenic Neisseria, p. 269-160.
In
S. Baumberg, J. P. W. Young, E. M. H. Wellington, and J. R. Saunders (ed.), Population genetics of bacteria. University of Cambridge, Cambridge, United Kingdom.
|
| 48.
|
Tettelin, H.,
N. J. Saunders,
J. Heidelberg,
A. C. Jeffries,
K. E. Nelson,
J. A. Eisen,
K. A. Ketchum,
D. W. Hood,
J. F. Peden,
R. J. Dodson,
W. C. Nelson,
M. L. Gwinn,
R. DeBoy,
J. D. Peterson,
E. K. Hickey,
D. H. Haft,
S. L. Salzberg,
O. White,
R. D. Fleischmann,
B. A. Dougherty,
T. Mason,
A. Ciecko,
D. S. Parksey,
E. Blair,
H. Cittone,
E. B. Clark,
M. D. Cotton,
T. R. Utterback,
H. Khouri,
H. Qin,
J. Vamathevan,
J. Gill,
V. Scarlato,
V. Masignani,
M. Pizza,
G. Grandi,
L. Sun,
H. O. Smith,
C. M. Fraser,
E. R. Moxon,
R. Rappuoli, and R. J. C. Vente.
2000.
Complete genome sequence of Neisseria meningitidis serogroup B strain MC58.
Science
287:1809-1815[Abstract/Free Full Text].
|
| 49.
|
Tinsley, C. R., and X. Nassif.
1996.
Analysis of the genetic differences between Neisseria meningitidis and Neisseria gonorrhoeae: two closely related bacteria expressing two different pathogenicities.
Proc. Natl. Acad. Sci. USA
93:11109-11114[Abstract/Free Full Text].
|
| 50.
|
Turk, D. C.
1982.
Clinical importance of Haemophilus influenzae 1981, p. 3-9.
In
S. H. Sell, and P. F. Wright (ed.), Haemophilus influenzae epidemiology, immunology, and prevention of disease. Elsevier/North-Holland Publishing Co., New York, N.Y.
|
| 51.
|
Weiser, J. N.
2000.
The generation of diversity by Haemophilus influenzae.
Trends Microbiol.
8:433-435[CrossRef][Medline].
|
Journal of Bacteriology, August 2001, p. 4626-4635, Vol. 183, No. 15
0021-9193/01/$04.00+0 DOI: 10.1128/JB.183.15.000-000.2001
Copyright © 2001, American Society for Microbiology. All rights reserved.
This article has been cited by other articles:
-
Satola, S. W., Napier, B., Farley, M. M.
(2008). Association of IS1016 with the hia Adhesin Gene and Biotypes V and I in Invasive Nontypeable Haemophilus influenzae. Infect. Immun.
76: 5221-5227
[Abstract]
[Full Text]
-
Brockmeyer, J., Bielaszewska, M., Fruth, A., Bonn, M. L., Mellmann, A., Humpf, H.-U., Karch, H.
(2007). Subtypes of the Plasmid-Encoded Serine Protease EspP in Shiga Toxin-Producing Escherichia coli: Distribution, Secretion, and Proteolytic Activity. Appl. Environ. Microbiol.
73: 6351-6359
[Abstract]
[Full Text]
-
Fernaays, M. M., Lesse, A. J., Cai, X., Murphy, T. F.
(2006). Characterization of igaB, a Second Immunoglobulin A1 Protease Gene in Nontypeable Haemophilus influenzae.. Infect. Immun.
74: 5860-5870
[Abstract]
[Full Text]
-
Fernaays, M. M., Lesse, A. J., Sethi, S., Cai, X., Murphy, T. F.
(2006). Differential Genome Contents of Nontypeable Haemophilus influenzae Strains from Adults with Chronic Obstructive Pulmonary Disease.. Infect. Immun.
74: 3366-3374
[Abstract]
[Full Text]
-
Erwin, A. L., Nelson, K. L., Mhlanga-Mutangadura, T., Bonthuis, P. J., Geelhood, J. L., Morlin, G., Unrath, W. C. T., Campos, J., Crook, D. W., Farley, M. M., Henderson, F. W., Jacobs, R. F., Muhlemann, K., Satola, S. W., van Alphen, L., Golomb, M., Smith, A. L.
(2005). Characterization of Genetic and Phenotypic Diversity of Invasive Nontypeable Haemophilus influenzae. Infect. Immun.
73: 5853-5863
[Abstract]
[Full Text]
-
Harrison, A., Dyer, D. W., Gillaspy, A., Ray, W. C., Mungur, R., Carson, M. B., Zhong, H., Gipson, J., Gipson, M., Johnson, L. S., Lewis, L., Bakaletz, L. O., Munson, R. S. Jr.
(2005). Genomic Sequence of an Otitis Media Isolate of Nontypeable Haemophilus influenzae: Comparative Study with H. influenzae Serotype d, Strain KW20. J. Bacteriol.
187: 4627-4636
[Abstract]
[Full Text]
-
Gladitz, J., Shen, K., Antalis, P., Hu, F. Z., Post, J. C., Ehrlich, G. D.
(2005). Codon usage comparison of novel genes in clinical isolates of Haemophilus influenzae. Nucleic Acids Res
33: 3644-3658
[Abstract]
[Full Text]
-
Shen, K., Antalis, P., Gladitz, J., Sayeed, S., Ahmed, A., Yu, S., Hayes, J., Johnson, S., Dice, B., Dopico, R., Keefe, R., Janto, B., Chong, W., Goodwin, J., Wadowsky, R. M., Erdos, G., Post, J. C., Ehrlich, G. D., Hu, F. Z.
(2005). Identification, Distribution, and Expression of Novel Genes in 10 Clinical Isolates of Nontypeable Haemophilus influenzae. Infect. Immun.
73: 3479-3491
[Abstract]
[Full Text]
-
Henderson, I. R., Navarro-Garcia, F., Desvaux, M., Fernandez, R. C., Ala'Aldeen, D.
(2004). Type V Protein Secretion Pathway: the Autotransporter Story. Microbiol. Mol. Biol. Rev.
68: 692-744
[Abstract]
[Full Text]
-
Munson, R. S. Jr., Harrison, A., Gillaspy, A., Ray, W. C., Carson, M., Armbruster, D., Gipson, J., Gipson, M., Johnson, L., Lewis, L., Dyer, D. W., Bakaletz, L. O.
(2004). Partial Analysis of the Genomes of Two Nontypeable Haemophilus influenzae Otitis Media Isolates. Infect. Immun.
72: 3002-3010
[Abstract]
[Full Text]
-
Leyton, D. L., Sloan, J., Hill, R. E., Doughty, S., Hartland, E. L.
(2003). Transfer Region of pO113 from Enterohemorrhagic Escherichia coli: Similarity with R64 and Identification of a Novel Plasmid-Encoded Autotransporter, EpeA. Infect. Immun.
71: 6307-6319
[Abstract]
[Full Text]
-
Kinsella, R. J., Fitzpatrick, D. A., Creevey, C. J., McInerney, J. O.
(2003). Fatty acid biosynthesis in Mycobacterium tuberculosis: Lateral gene transfer, adaptive evolution, and gene duplication. Proc. Natl. Acad. Sci. USA
100: 10320-10325
[Abstract]
[Full Text]
-
Rose, J. E., Meyer, D. H., Fives-Taylor, P. M.
(2003). Aae, an Autotransporter Involved in Adhesion of Actinobacillus actinomycetemcomitans to Epithelial Cells. Infect. Immun.
71: 2384-2393
[Abstract]
[Full Text]
-
Williams, B. J., Golomb, M., Phillips, T., Brownlee, J., Olson, M. V., Smith, A. L.
(2002). Bacteriophage HP2 of Haemophilus influenzae. J. Bacteriol.
184: 6893-6905
[Abstract]
[Full Text]