Horizontal gene transfer in Histophilus somni and its role in the evolution of pathogenic strain 2336, as determined by comparative genomic analyses
- Shivakumara Siddaramappa1, 2,
- Jean F Challacombe2,
- Alison J Duncan1,
- Allison F Gillaspy3,
- Matthew Carson3,
- Jenny Gipson3,
- Joshua Orvis3,
- Jeremy Zaitshik3,
- Gentry Barnes3,
- David Bruce2,
- Olga Chertkov2,
- J Chris Detter2,
- Cliff S Han2,
- Roxanne Tapia2,
- Linda S Thompson2,
- David W Dyer3 and
- Thomas J Inzana1Email author
© Siddaramappa et al; licensee BioMed Central Ltd. 2011
Received: 15 June 2011
Accepted: 23 November 2011
Published: 23 November 2011
Pneumonia and myocarditis are the most commonly reported diseases due to Histophilus somni, an opportunistic pathogen of the reproductive and respiratory tracts of cattle. Thus far only a few genes involved in metabolic and virulence functions have been identified and characterized in H. somni using traditional methods. Analyses of the genome sequences of several Pasteurellaceae species have provided insights into their biology and evolution. In view of the economic and ecological importance of H. somni, the genome sequence of pneumonia strain 2336 has been determined and compared to that of commensal strain 129Pt and other members of the Pasteurellaceae.
The chromosome of strain 2336 (2,263,857 bp) contained 1,980 protein coding genes, whereas the chromosome of strain 129Pt (2,007,700 bp) contained only 1,792 protein coding genes. Although the chromosomes of the two strains differ in size, their average GC content, gene density (total number of genes predicted on the chromosome), and percentage of sequence (number of genes) that encodes proteins were similar. The chromosomes of these strains also contained a number of discrete prophage regions and genomic islands. One of the genomic islands in strain 2336 contained genes putatively involved in copper, zinc, and tetracycline resistance. Using the genome sequence data and comparative analyses with other members of the Pasteurellaceae, several H. somni genes that may encode proteins involved in virulence (e.g., filamentous haemaggutinins, adhesins, and polysaccharide biosynthesis/modification enzymes) were identified. The two strains contained a total of 17 ORFs that encode putative glycosyltransferases and some of these ORFs had characteristic simple sequence repeats within them. Most of the genes/loci common to both the strains were located in different regions of the two chromosomes and occurred in opposite orientations, indicating genome rearrangement since their divergence from a common ancestor.
Since the genome of strain 129Pt was ~256,000 bp smaller than that of strain 2336, these genomes provide yet another paradigm for studying evolutionary gene loss and/or gain in regard to virulence repertoire and pathogenic ability. Analyses of the complete genome sequences revealed that bacteriophage- and transposon-mediated horizontal gene transfer had occurred at several loci in the chromosomes of strains 2336 and 129Pt. It appears that these mobile genetic elements have played a major role in creating genomic diversity and phenotypic variability among the two H. somni strains.
Histophilus somni is a commensal or opportunistic pathogen of the reproductive and respiratory tracts of cattle. H. somni was initially identified as the etiologic agent of bovine thrombotic meningoencephalitis (TME), but also causes bovine shipping fever pneumonia, either independently or in association with Mannheimia haemolytica and Pasteurella multocida. Pneumonia and myocarditis are currently the most commonly reported diseases due to H. somni. Infections resulting in abortion, infertility, arthritis, septicemia, and mastitis can also be caused by H. somni with varying degrees of frequency and severity in cattle . Similar disease conditions associated with strains of H. somni have been described in sheep . Relatively less pathogenic and/or avirulent variants of H. somni have also been isolated from cattle, most frequently from the mucosal surfaces of the genital tract .
Numerous in vitro and in vivo studies during the pre-genomic era have shed light on the differences in virulence properties between H. somni pathogenic isolates from sick animals and serum-sensitive commensal isolates from the genital tract . However, thus far only a few genes involved in lipooligosaccharide (LOS) biosynthesis and serum-resistance have been identified in H. somni using DNA/DNA and DNA/protein comparisons [5–7]. H. somni pneumonia strain 2336 and preputial strain 129Pt have been comprehensively characterized phenotypically and have been analyzed in several comparative studies [8–10]. However, a comprehensive understanding of the genetic basis that determines the phenotypic variability among H. somni stains is necessary to gain further insights into their pathogenicity.
Comparative (in silico) analysis of bacterial genomes is a powerful tool for the prediction and/or identification of biochemical differences, virulence attributes, pathogenic ability, and adaptive evolution among related species/strains . Among the Pasteurellaceae, the genomes of one or more species pathogenic to humans or animals from the genera Actinobacillus, Haemophilus, Mannheimia, Pasteurella, and others have been sequenced. The availability of these genome sequences has facilitated whole genome comparisons that have provided insights into the physiology and pathogenic evolution of the corresponding bacteria [12, 13].
Horizontal gene transfer (HGT: defined as the "acquisition of new genes either directly by transformation with naked DNA, transduction with phages, or the uptake of plasmids or chromosomal fragments by conjugation") plays a critical role in driving the evolution of pathogenic bacteria . Reduction in genome size (referred to as reductive evolution) can occur as a result of continuous loss of genetic material due to gene deletion and/or mutation followed by DNA erosion . Previous analyses by biochemical and pulsed field gel electrophoresis indicated that H. somni strains 2336 and 129Pt have common ancestry, but are non-clonal [16, 17]. The following mechanisms may have engendered the genetic differences between these strains: (i) only one strain acquired genes by HGT while the other one did not; (ii) only one strain lost genes by deletion/mutation and underwent 'reductive evolution'; (iii) both strains independently and continuously acquired and lost genes, and the net loss or gain of genes is a determinant of their divergent evolution; (iv) gene convergence and the accumulation of synonymous and/or nonsynonymous nucleotide substitutions occurred across the genomes of the two strains.
The rationale for the present study was to determine, using whole genome sequencing and comparative genomics, the mechanisms responsible for genetic variability between the two strains. It was also envisaged that a comparative genomics and bioinformatics approach would facilitate identification of H. somni genes putatively involved in virulence and pathogenesis.
Genomic DNA (2 mg) from H. somni strain 2336 was purified using the Puregene protocol (Gentra Systems, Minneapolis, MN). The shotgun sequencing phase for this genome required ~35,200 sequence reads to reach 8-fold coverage . Library construction, template preparation, sequencing, assembly, and data analyses were performed as described previously [19, 20]. The sequence data assembled with Phred-Phrap were viewed using Consed to assess data quality and design closure experiments. Consed was also used to identify putative repeat regions so that the problems associated with assembling these regions could be resolved by way of combinatorial PCR experiments to isolate the repeat sequences on PCR amplicons. The location and exact sequence of each repeat was confirmed by isolating PCR fragments that contained each repeat in its entirety, followed by primer walking across the PCR product.
For initial gap closure, Single Primer Amplification of Contig Ends (SPACE), which is similar to the single-primer PCR procedure for rapid identification of transposon insertion sites, was used . Additional primers were designed, as necessary, to verify the correct assembly of contigs by confirmatory PCR. Simultaneously, a fosmid library was constructed for scaffolding purposes using the vector pCC1fos (Epicentre Biotechnologies, Madison, WI) with 40 kb inserts. Sequencing of the fosmids was necessary to close gaps across sequences that occur more than once in the genome, such as those of insertion sequences and ribosomal genes. Gaps that were not closed by SPACE-walking were closed using the sequence of H. somni strain 129Pt as a scaffold and the reads were assembled with parallel phrap (High Performance Software, LLC). Gap closure at this stage was also facilitated by AUTOFINISH . Possible mis-assemblies were corrected with Dupfinisher  or transposon bombing of bridging clones  using an EZ::TN™ kit (Epicentre). The National Human Genome Research Institute standards for the Human Genome Project (1 error per 10, 000 assembled bases) were followed for H. somni to obtain sufficient quality genomic data.
Final automated annotation of the genome of strain 2336 was performed at the Oak Ridge National Laboratory using methods similar to those used to annotate the strain 129Pt genome . Briefly, protein domains were identified by comparing each predicted protein against a Hidden Markov Model protein family database . To estimate the number of proteins specific to each strain, the Smith-Waterman algorithm  was used to compare all predicted proteins from strain 129Pt against those from strain 2336 and vice-versa. Proteins deemed to be specific to each strain were compared against the NCBI non-redundant protein database to determine whether they were hypothetical or conserved hypothetical. The translated ORF was named a hypothetical protein if there was less than 25% identity or an aligned region was less than 25% of the predicted protein length. Prediction of the number of subsystems and pairwise BLAST comparisons of protein sets within strains 2336 and 129Pt were carried out with the Rapid Annotation using Subsystems Technology (RAST), which is a fully automated, prokaryotic genome annotation service . This platform identifies tRNA and rRNA genes using the tools tRNAscan-SE and "search_for_rnas", respectively .
Multiple genome comparisons were performed using the 'progressive alignment' option available in the program MAUVE version 2.3.0 [28, 29]. Default scoring and parameters were used for generating the alignment. A synteny plot was generated using the program NUCmer, which creates a dot plot based on the number of identical alignments between two genomes . Prophage regions (PRs) were identified using Prophinder http://aclame.ulb.ac.be/Tools/Prophinder/, an algorithm that combines similarity searches, statistical detection of phage-gene enriched regions, and genomic context for prophage prediction . Identification and annotation of genomic islets/islands (GIs) other than prophages were performed based on sequence composition bias and comparative genomic analysis . Briefly, differences in the GC content, the occurrence of 'cornerstone genes' (e.g., transposases), and/or a continuous stretch of genes encoding hypothetical proteins were used as reference points for detection of GIs. Insertion sequences (ISs) were identified by whole genome BLASTX analysis of strains 2336 and 129Pt using the IS finder http://www-is.biotoul.fr/. Gene acquisition and loss among the two strains was determined by comparing gene order, orientation of genes (forward/reverse), GC content of genes (percent above or below whole genome average), features of intergenic regions (e.g., remnants of IS elements, integration sites), and the similarity of proteins encoded by genes at a locus of interest (> 90% identity at the predicted protein level). Putative horizontally transferred genes (HTGs; defined as genes whose greatest homology, based on BLASTP scores, is to genes from a more distant phylogenetic group than to genes from the same or a close phylogenetic group as the query genome) were compiled using the integrated microbial genomes (IMG) system http://img.jgi.doe.gov. This system uses not only the best hit (i.e., the homolog with the highest bitscore), but also all the matches that have bitscores equal to or greater than that of the best hit to identify putative HTGs . DNA and protein sequences were aligned using the ClustalW http://www.ebi.ac.uk/Tools/clustalw2/index.html and BOXSHADE http://www.ch.embnet.org/software/BOX_form.html programs as described previously .
Properties of the chromosomes
Characteristics of the genomes of H.somni strains
Number of protein coding genes
Overall coding density
*Number of subsystems
5S ribosomal RNA genes
16S ribosomal RNA genes
23S ribosomal RNA genes
Number of tRNA genes
pHS129 (5178 bp) 
GenBank accession number
Comparison of prophage regions and genomic islands
Characteristics of the prophage regions and genomic islands of H.somni strains 2336 and 129Pt
Number of genes
243787 bp to 255490 bp
759194 bp to 783660 bp
1001039 bp to 1049028 bp
1301649 bp to 1329985 bp
1877815 bp to 1902856 bp
1963256 bp to 1992078 bp
2124957 bp to 2143800 bp
464154 bp to 472854 bp
567352 bp to 574742 bp
1524138 bp to 1529684 bp
1569138 bp to 1610766 bp
1955573 bp to 1961650 bp
1962205 bp to 1966694 bp
1969857 bp to 1974596 bp
Manual curation indicated that strains 2336 and 129Pt contained 3 and 6 GIs, comprising a total of 72,709 bp and 36,947 bp, respectively (Table 2). GI I of strain 2336 contained 17 ORFs of unknown function and a similar sequence was not found in strain 129Pt or other members of the Pasteurellaceae. GI II of strain 2336 was the longest, had a higher GC content (~43%) than the overall genome, and contained 16 ORFs of unknown function. The genome of P. multocida strain Pm70 contained several short sequences that had homology to this region (e.g., 3 segments totaling 5,433 bp in the 2,174,533 bp to 2,187,234 bp region, 82-99% nucleotide identity, E-value = 0). GI III of strain 2336 contained 4 ORFs encoding putative transposases ([GenBank:HSM_1871, [GenBank:HSM_1874, [GenBank:HSM_1883, and [GenBank:HSM_1887] and 8 ORFs of unknown function. Strain 129Pt lacked an analogous GI, but contained homologs of the genes encoding putative transposases. A 4,640 bp sequence from Escherichia coli strain 49 [GenBank:U23723] had homology to GI III of strain 2336 (72% identity across 3,270 bp, E-value = 0).
GI I of strain 129Pt contained an ORF encoding a resolvase/integrase-like protein [GenBank:HS_0445] and 6 ORFs of unknown function. A similar sequence was not found in strain 2336, but H. parasuis strain SH0165 contained several short sequences with homology to this region (e.g., 3 segments totaling 2,040 bp, 68-79% nucleotide identity, E-values = 0 to 1e-06). GI II of strain 129Pt contained an ORF encoding a putative phage DNA primase-like protein [GenBank:HS_0533] and 9 ORFs of unknown function. A similar sequence was not found in strain 2336 or other members of the Pasteurellaceae. GI III of strain 129Pt contained ORFs encoding a putative phage terminase protein [GenBank:HS_1334], a prophage regulatory element [GenBank:HS_1335], an integrase [GenBank:HS_1337], and 4 ORFs of unknown function. A similar sequence was not found in strain 2336, but Aggregatibacter aphrophilus strain NJ8700 contained several short sequences that had homology to this region (e.g., 2 segments totaling 1,334 bp in the 1,818,930 bp to 1,821,098 bp region, 74-80% nucleotide identity, E-values = 0 to 4e-35). All predicted ORFs from GIs IV and V of strain 129Pt were of unknown function and strain 2336 contained several short sequences that had homology to these regions (e.g., 7 segments totaling 2,184 bp in the 1,965,862 bp to 1,990,890 bp region of GI II, 67-92% nucleotide identity, E-values = 9e-42 to 1e-10). The chromosome of strain 2336 had no regions of homology to GI VI of strain 129Pt, but H. parasuis strain SH0165 contained several short sequences that had homology to this region (e.g., 3 segments totaling 1,573 bp, 68% nucleotide identity, E-values = 2e-51 to 2e-13).
Comparison of insertion sequences
Insertion sequence finder indicated that strains 2336 and 129Pt contained several IS elements distributed throughout the chromosomes. Not surprisingly, some of these IS elements were found within the PRs and/or GIs described above. Insertion sequence 1016 (IS1016), consisting of a transposase [217 amino acids (aa)] flanked by 18-29 bp of terminal inverted repeats, is a member of the IS1595 superfamily [35, 36]. Strain 2336 contained 4 full-length ([GenBank:HSM_0851], [GenBank:HSM_1211], [GenBank:HSM_1267], and [GenBank:HSM_1883]) and three truncated copies of IS1016. GI III of strain 2336 included [GenBank:HSM_1883] and one of the three truncated copies of IS1016. Strain 129Pt contained 5 full-length copies of IS1016 ([GenBank:HS_0324], [GenBank:HS_0642], [GenBank:HS_1116], [GenBank:HS_1679], and [GenBank:HS_1709]). A gene in strain 129Pt (licA; [GenBank:HS_1461]), encoding a choline kinase that is necessary for the synthesis of phosphorylcholine , appeared to be interrupted due to an insertion-excision event involving IS1016, since two partial homologs of the transposase gene occurred upstream of licA.
Strain 2336 contained 4 full-length members of the family IS200/IS605 ([GenBank:HSM_0223] was near PR I and [GenBank:HSM_1626] and [GenBank:HSM_1631] were near GI I). Strain 129Pt contained 3 full-length ([GenBank:HS_0486], [GenBank:HS_0583], and [GenBank:HS_1390]) and 6 truncated ([GenBank:HS_0224], [GenBank:HS_0666], [GenBank:HS_0667], [GenBank:HS_0668], [GenBank:HS_0716], and [GenBank:HS_0717]) members of the family IS200/IS605. PR I of strain 129Pt included [GenBank:HS_1390]. Strain 2336 contained a full-length member of the family IS30 ([GenBank:HSM_1680]) and a truncated member of the family IS1182/IS5 ([GenBank:HSM_1887]). Strain 2336 contained a full-length member of the family IS481 ([GenBank:HSM_0451]), whereas strain 129Pt contained a truncated copy of the same ORF ([GenBank:HS_1518]). Strain 2336 also contained truncated members of the family IS3 ([GenBank:HSM_0531] and [GenBank:HSM_0532]). The closest homologs of H. somni IS200/IS605 and IS1595 elements were found in M. haemolytica (e.g., [GenBank:MHA_0340], 92% identity, E-value = 0) and H. influenzae (e.g., [GenBank:CGSHiII_08436], 86% identity, E-value = 0). However, the closest homologs of the IS30 element from strain 2336 were found in A. pleuropneumoniae (e.g., [GenBank:APP7_0397], 99% identity, E-value = 0) and H. parasuis (e.g., [GenBank:HAPS_1772], 99% identity, E-value = 8e-180).
BLAST comparison of protein sets
Strains 2336 and 129Pt contained 1550 predicted protein coding genes in common (bidirectional best hits, at least 90% identity at the predicted protein level). In strain 2336, 440 ORFs could not be assigned a function based on BLAST analysis and were therefore annotated as encoding hypothetical or conserved hypothetical proteins. In strain 129Pt, 429 ORFs were annotated as encoding hypothetical or conserved hypothetical proteins. Among hypothetical proteins that were common to both strains, 30 did not have homologs outside the genus. Pairwise BLAST comparisons indicated that strain 2336 contained 311 putative protein coding genes with no homologs in strain 129Pt (additional file 1). Within this subset, proteins encoded by 302 genes had at least 51 aa and 9 genes ([GenBank:HSM_0528], [GenBank:HSM_0530], [GenBank:HSM_0532], [GenBank:HSM_0603], [GenBank:HSM_1185], [GenBank:HSM_1483], [GenBank:HSM_1636], [GenBank:HSM_1742], and [GenBank:HSM_1743]) had 38-50 aa. Strain 129Pt contained 165 putative protein coding genes with no homologs in strain 2336 (additional file 2). Within this subset, proteins encoded by all 165 genes had at least 51 aa.
In both strains, a vast majority of putative HTGs appeared to have had their origins among members of gammaproteobacteria (mostly within Pasteurellales and Enterobacteriales). Putative HTGs with possible origins among members of betaproteobacteria (27 in strain 2336, 11 in strain 129Pt) and alphaproteobacteria (1 in strain 2336, 6 in strain 129Pt) were also identified. Among HTGs identified were those encoding proteins putatively involved in virulence (e.g., filamentous hemagglutinins, proteases, and antibiotic resistance regulators). A complete list of these genes is available at the 'Organism Details' sections for strains 2336 and 129Pt within IMG. Other strain-specific genes identified encoded DNA methylases (7 in strain 2336, none in strain 129Pt), transposases (8 in strain 2336, 1 in strain 129Pt), ABC transporters (5 in strain 2336, none in strain 129Pt), ATPases (4 in strain 2336, none in strain 129Pt), transcriptional regulators (14 in strain 2336, 3 in strain 129Pt), kinases (2 in strain 2336, 1 in strain 129Pt), and several proteins related to bacteriophage functions (e.g., of the 10 integrase/resolvase-related genes found in strain 129Pt, six have no homologs in strain 2336 and of the 7 integrase/resolvase-related genes found in strain 2336, 2 have no homologs in strain 129Pt). Excluding intergenic regions, the total length of sequence that was associated with specific genes in strain 2336 was 254,052 bp (~11% of the genome), and was 98,016 bp (~5% of the genome) in specific genes of strain 129Pt.
Identification of genes encoding polysaccharide biosynthesis/modification enzymes
Putative glycosyltransferase genes of H.somni strains
Strain 2336 locus tag, protein,1(SSR, if present)n
Strain 129Pt locus tag, protein,1(SSR, if present)n
SWISS-PROT accession number
Closest Pasteurellaceae homolog Bacterium, locus tag, protein size,2identity
Closest GT homolog (Outside Pasteurellaceae) locus tag, protein,2identity
[GenBank:HSM_0148], 260 aa, (G)9
[GenBank:HS_0275], 259 aa, (G)6
Acinetobacter johnsonii, [GenBank:EEY94992], 281 aa, 44%
[GenBank:HSM_0164], 522 aa, (G)9
Actinobacillus succinogenes, [GenBank:Asuc_0521], 522 aa, 63
Neisseria mucosa, [GenBank:EFC88975], 517 aa, 53%
3[GenBank:HSM_0398], 346 aa
[GenBank:HS_1612], 346 aa
Mannheimia, succiniciproducens [GenBank:MS2260], 346 aa, 78%
Yersinia pseudotuberculosis, [GenBank:YPTB0054], 354 aa, 64%
3[GenBank:HSM_0399], 321 aa
[GenBank:HS_1611], 321 aa
Providencia rettgeri, [GenBank:EFE51471], [GenBank:MS2259], 321 aa, 67% 322 aa, 63%
4[GenBank:HSM_0856], 314 aa
Rhodopirellula baltica, [GenBank:RB9243], 326 aa, 35%
5[GenBank:HSM_0975a], 157 aa, (CAGT)18
5[GenBank:HS_0636], 354 aa, (CAGT)29
Haemophilus influenzae, [GenBank:NTHI0365], 312 aa, 38%
Neisseria meningitidis, [GenBank:AF355193_3], 311 aa, 35%
[GenBank:HSM_0975], 231 aa, (AGA)3
Pasteurella dagmatis, [GenBank:EEX51329], 334 aa, 40%
Neisseria lactamica, [GenBank:AAN08512], 251aa, 32%
6[GenBank:HSM_0976], 63 aa
Pasteurella multocida, [GenBank:PM1140], 337 aa, 61%
Streptococcus parasanguinis, [GenBank:ACF35266], 297 aa, 59%
7[GenBank:HSM_0977], 277 aa, (GAGA)11
Aggregatibacter actinomycetemcomitans GenBank:D11S_0429], 282 aa, 43%
Neisseria meningitidis, [GenBank:NMC1901], 275 aa, 41%
8[GenBank:HSM_0978], 271 aa, (GA)20
[GenBank:HS_0637], 158 aa
Aggregatibacter actinomycetemcomitans [GenBank:D11S_0429], 282 aa, 43%
Neisseria meningitidis, [GenBank:CAX50883], 280 aa, 39%
9[GenBank:HSM_0979], 293 aa, (CAAT)33
[GenBank:HS_0638], 297 aa, (CAAT)36
Actinobacillus succinogenes, [GenBank:Asuc_0526], 265 aa, 52%
Photobacterium profundum, [GenBank:PBPRA0217], 253 aa, 41%
[GenBank:HSM_1426], 300 aa
Haemophilus influenzae, [GenBank:HI0871], 306 aa, 42%
Fusobacterium nucleatum, [GenBank:EDK88896], 326 aa, 38%
[GenBank:HSM_1552], 196 aa, (G)9
[GenBank:HS_1067], 194 aa, (G)6
Edwardsiella ictaluri, [GenBank:NT01EI_2547], 247 aa, 27%
10[GenBank:HSM_1794], 311 aa
Mannheimia haemolytica, [GenBank:MHA_0104], 312 aa, 47%
Streptococcus agalactiae [GenBank:AAR29926], 299 aa, 39%
[GenBank:HSM_2001], 242 aa, (GGT)3
[GenBank:HS_0116], 242 aa, (GGT)3
Aggregatibacter aphrophilus, [GenBank:NT05HA_1706], 240 aa, 61%
[GenBank:EER56028], 241 aa, 55%
[GenBank:HS_0291], 256 aa
Mannheimia succiniciproducens [GenBank:MS0438], 253 aa, 76%
Vibrio shilonii [GenBank:EDL53925], 259 aa, 62%
[GenBank:HS_0292], 348 aa
Actinobacillus minor NM305, [GenBank:EER47697], 340 aa, 57%
Vibrio parahaemolyticus, [GenBank:EED27646], 341 aa, 39%
The l ipooligo saccharide b iosynthesis (lob) gene cluster consisting of lob1 and lob2ABCD ORFs encoding glycosyltransferases involved in attaching the outer core glycoses of the LOS was previously identified in strain 738, which is an LOS phase variant of strain 2336 [5, 7]. Strain 129Pt encoded full-length homologs of lob1 and lob2D, but only the 5' ends of lob2A and lob2C and lacked lob2B. Strain 2336 contained full-length homologs of lob1 and lob2ABC, but a truncated homolog of lob2D (table 3). The variations in the lob loci in the genomes of strains 2336 and 129Pt correlate with the differences in the structures of the LOS of these strains, as determined by NMR spectroscopy and mass spectrometry [16, 38]. In addition, both strains contained ORFs encoding a phosphoheptose isomerase ([GenBank:HSM_0840] and [GenBank:HS_1238], gmhA), D,D-heptose 1,7-bisphosphate phosphatase ([GenBank:HSM_0572] and [GenBank:HS_1532], gmhB), ADP-L-glycero-D-mannoheptose-6-epimerase ([GenBank:HSM_0397] and [GenBank:HS_1613], rfaD), bifunctional heptose 7-phosphate kinase/heptose 1-phosphate adenyltransferase ([GenBank:HSM_0925] and [GenBank:HS_0576, rfaE), and bifunctional N-acetylglucosamine-1-phosphate uridyltransferase/glucosamine-1-phosphate acetyltransferase ([GenBank:HSM_0204] and [GenBank:HS_0333], glmU). These five genes/enzymes were predicted to be involved in LOS core biosynthesis in strains 2336 and 129Pt.
Strain 2336 also contained a locus that encodes proteins putatively involved in exopolysaccharide (EPS) and/or LOS biosynthesis ([GenBank:HSM_1061] or csrA, carbon storage regulator; [GenBank:HSM_1062] or manB phosphomannomutase; [GenBank:HSM_1063] or galU, UTP-glucose-1-phosphate uridylyltransferase). Strain 129Pt contained a similar locus ([GenBank:HS_1117] to [GenBank:HS_1119]) that had an IS1016 element ([GenBank:HS_1116], transposase) adjacent to it. Furthermore, both strains contained ORFs encoding a putative UDP-glucose 4-epimerase ([GenBank:HSM_1256] and [GenBank:HS_0789], galE), phosphoglucomutase ([GenBank:HSM_1832] and [GenBank:HS_1670], pgmB), and dTDP-glucose 4,6-dehydratase ([GenBank:HSM_1118] and [GenBank:HS_0707], rmlB). These three genes/enzymes were predicted to be involved in polysaccharide biosynthesis/modification.
Acquisition and loss of genes encoding filamentous hemagglutinins
Locus III of strain 2336 was 14,066 bp (GC content of 37%) and contained genes encoding FhaB ([GenBank:HSM_1489]) and FhaC ([GenBank:HSM_1490]). The fhaB (12,288 bp) of this locus was the second largest gene in the genome and the largest among the 12 homologs. This gene encoded a putative protein homologous to the high molecular weight immunoglobulin-binding protein of H. somni and the large supernatant proteins (Lsp1 and Lsp2) of H. ducreyi that have been previously described . No transposon or phage regions were apparent in this locus. Strain 129Pt had a locus containing genes that flank locus III of strain 2336, but did not contain the FhaB and FhaC homologs (Figure 5, locus III). Locus IV of strain 2336 was 21,841 bp (GC content of 38%) and contained genes encoding six FhaB homologs ([GenBank:HSM_1638], [GenBank:HSM_1641], [GenBank:HSM_1643], [GenBank:HSM_1646], [GenBank:HSM_1647], [GenBank:HSM_1651]). A truncated gene encoding FhaC ([GenBank:HSM_1653], 98 aa) was also found in this locus. A partial homolog of [GenBank:HSM_1647] was found in strain 129Pt ([GenBank:HS_0540], 153 aa), which appeared to be the only region in the chromosome of strain 129Pt with a sequence related to the fhaB genes (Figure 5, locus IV).
Acquisition of a gene encoding subtilisin-like protease
H. somni strain 2336 GI III contained a gene encoding a putative subtilisin-like serine protease ([GenBank:HSM_1889], 730 aa; Figure 5, locus V, arrow (gene) 10). A BLASTP search revealed that the non-redundant GenBank database contained several proteins that were homologous to [GenBank:HSM_1889], with the closest relative being an E. coli protein ([GenBank:AAA64865], 75% identity, E-value = 0). The NCBI protein clusters database lists [GenBank:HSM_1889] as one of the 16 members of the [GenBank:CLSK923564] cluster (peptidase S8 and S53, subtilisin, kexin, sedolisin). Among the proteins of this cluster, homologs from Ralstonia eutropha JMP134 ([GenBank:Reut_C6419], 45% identity, E-value = 1e-178), Pseudomonas fluorescens Pf0-1 ([GenBank:Pfl01_5697], 35% identity, E-value = 2e-116), and Pseudomonas syringae pv. phaseolicola 1448A ([GenBank:PSPPH_0180], 35% identity, E-value = 1e-115) are listed in the prokaryotic subtilase database as subtilases of the D-H-S family.
Comparison of genes encoding transferrin-binding proteins
A comparison of the genomes of H. somni strains 2336 and 129Pt revealed that both strains contained genes encoding TbpA ([GenBank:HS_0449] and [GenBank:HSM_0750]) and TbpB ([GenBank:HS_0448] and [GenBank:HSM_0749]). In strain 129Pt, GI I was identified immediately upstream of the tbp locus. Furthermore, homologs of H. somni strain 649 TbpA2 were present in strain 129Pt ([GenBank:HS_0582], heme uptake protein) and strain 2336 ([GenBank:HSM_0931], TonB-dependent receptor; [GenBank:HSM_0932], TonB-dependent receptor plug; and [GenBank:HSM_1988], TonB-dependent lactoferrin and transferrin receptor). Pairwise BLAST analysis identified several other putative genes encoding proteins that may also be involved in iron transport. In strain 129Pt these included [GenBank:HS_0069], iron-regulated outer membrane protein; [GenBank:HS_0181], TonB-dependent outer membrane receptor; [GenBank:HS_0728], hemin receptor outer membrane protein; and [GenBank:HS_1306], iron-regulated outer membrane protein. In strain 2336 they included [GenBank:HSM_0047], TonB-dependent receptor; [GenBank:HSM_1168], TonB-dependent hemoglobin/transferrin/lactoferrin family receptor; [GenBank:HSM_1176], outer membrane hemin receptor protein, and [GenBank:HSM_1962], TonB-dependent receptor. Two of the genes, [GenBank:HS_1306] and [GenBank:HSM_1168], were strain-specific. Among others, HS_0582 was associated with a transposase ([GenBank:HS_0583]) and [GenBank:HSM_1168] was found near PR IV.
Overabundance of genes encoding adhesin-like proteins
Strain 129Pt contained sixteen genes randomly distributed throughout the chromosome ([GenBank:HS_0209], [GenBank:HS_0383], [GenBank:HS_0478], [GenBank:HS_0589], [GenBank:HS_0602], [GenBank:HS_0790], [GenBank:HS_1058], [GenBank:HS_1085], [GenBank:HS_1154], [GenBank:HS_1185], [GenBank:HS_1234], [GenBank:HS_1512], [GenBank:HS_1543], [GenBank:HS_1563], [GenBank:HS_1616], and [GenBank:HS_1632]) with homology to genes that encode proteins of the Yersinia adhesin (YadA) superfamily . Strain 2336 also contained sixteen genes randomly distributed throughout the chromosome ([GenBank:HSM_0077], [GenBank:HSM_0338], [GenBank:HSM_0346], [GenBank:HSM_0377], [GenBank:HSM_0394], [GenBank:HSM_0708], [GenBank:HSM_0844], [GenBank:HSM_0938], [GenBank:HSM_0953], [GenBank:HSM_1022], [GenBank:HSM_1212], [GenBank:HSM_1257], [GenBank:HSM_1484], [GenBank:HSM_1542], [GenBank:HSM_1571], [GenBank:HSM_1793]) with homology to genes that encode proteins of the YadA superfamily. Adhesin-encoding genes [GenBank:HS_0209] (15,431 bp) and [GenBank:HSM_1257] (13,970 bp) were the largest among all protein coding genes predicted in the chromosomes of strains 129Pt and 2336, respectively. Within this gene repertoire, [GenBank:HSM_0938] (388 aa) and [GenBank:HS_0589] (386 aa) were 85% similar to each other at the predicted protein level (75% identity, Score = 535 bits, E-value = 1e-156) and were associated with genes encoding putative ABC transporters ([GenBank:HSM_0939-0943] and [GenBank:HS_0590-0594], respectively, 99% identity, E-value = 0 to 5e-151). In addition, [GenBank:HSM_0077] (4063 aa) and [GenBank:HS_0209] (5143 aa) were 65% similar to each other at the predicted protein level (53% identity, Score = 1009 bits, E-value = 0) and were also associated with genes encoding putative ABC transporters ([GenBank:HSM_0079-0081] and [GenBank:HS_0210-0212], respectively, 99% identity, E-value = 2e-146 to 3e-132). Furthermore, although homologs of the H. influenzae fimbrial gene cluster (hifABCDE) were absent in both strains of H. somni, homologs of type IV pili genes (pilABCD) occurred in strains 2336 and 129Pt (e.g., [GenBank:HSM_0123] and [GenBank:HS_0250], [GenBank:HSM_0217] and [GenBank:HS_1430, [GenBank:HSM_0755 and [GenBank:HS_0457], and [GenBank:HSM_0756] and [GenBank:HS_0458], 40-69% identity, E-values = 1e-168 to 1e-38). In addition, both strains contained a gene that encoded a putative pseudopilin (pulG, [GenBank:HS_0264] and [GenBank:HSM_0137], 54% identity to NTHI1109, E-value = 9e-51) that may facilitate type II secretion.
Identification of genes encoding transcriptional regulators and drug/metal resistance
Pairwise BLAST comparisons indicated that strain 2336 contained 14 genes encoding transcriptional regulators with no homologs in strain 129Pt. These included two regulators each that belong to the TetR ([GenBank:HSM_1191] and [GenBank:HSM_1734]) and MerR ([GenBank:HSM_1728] and [GenBank:HSM_1741]) families, one each that belonged to the LysR ([GenBank:HSM_0806]), OmpR ([GenBank:HSM_0817]), and MarR ([GenBank:HSM_1737]) families, and a member of an unassigned family ([GenBank:HSM_0489]; [GenBank:COG2865K]). Among these, members of the TetR, MerR, and MarR families were found within GI II. Strain 2336 also contained a gene adjacent to tetR encoding a tetracycline resistance antiporter (tetH; [GenBank:HSM_1735]). The NCBI protein clusters database lists [GenBank:HSM_1734] as a member of the PRK13756 cluster (TetR; tetracycline repressor of Leuconostoc citreum[GenBank:COG1309K]). The closest homologs of [GenBank:HSM_1734] and [GenBank:HSM_1735] were found in P. multocida (e.g., [GenBank:AAC43249] and [GenBank:AAC43250] encoding a repressor protein and a tetracycline resistance protein, respectively, 99% identity, E-values = 0 to 1e-116). In addition to [GenBank:HSM_1737], strain 2336 also contained [GenBank:HSM_1736] (encoding small multidrug resistance protein). [GenBank:HSM_1191] is listed as one of the members of the [GenBank:CLSK391246] cluster (TetR family transcriptional regulator). This regulator gene was associated with two genes encoding an ABC-type transport system ([GenBank:HSM_1192] and [GenBank:HSM_1193]). The closest homologs of [GenBank:HSM_1191] were found among members of Firmicutes, Spirochaetes, and Fusobacteria (e.g., [GenBank:EFM38698] encoding a putative TetR family transcriptional regulator in Eubacterium yurii subsp. margaretiae, 99% identity, E-value = 1e-104), but not the Pasteurellaceae. The two MerR homologs of strain 2336 had only 39% identity to each other (E-value = 7e-23) and both were conserved in members of the Pasteurellaceae (e.g., [GenBank:HAPS_1832] from H. parasuis and [GenBank:PM1941] from P. multocida, 85-99%% identity, E-values = 4e-57 to 4e-70). However, the MerR homologs have been included in different protein clusters ([GenBank:HSM_1728] in [GenBank:CLSK892364] and [GenBank:HSM_1741] in [GenBank:CLSK2299246]) and appear to regulate different functions. Since [GenBank:HSM_1728] was associated with genes involved in copper metabolism ([GenBank:HSM_1729] encodes a metal binding protein, [GenBank:HSM_1730] encodes a multicopper oxidase, and [GenBank:HSM_1731] encodes a copper-translocating ATPase), it may encode a copper efflux regulator, CueR. However, [GenBank:HSM_1741] may be involved in zinc homeostasis since it was associated with zinc-responsive genes ([GenBank:HSM_1740] encodes a zinc efflux protein and [GenBank:HSM_1739] encodes a zinc-dependent hydrolase). The NCBI protein clusters database lists [GenBank:HSM_0806] as one of the members of the [GenBank:CLSK797597] cluster (LysR family transcriptional regulator, [GenBank:COG0583K]). The closest homologs of [GenBank:HSM_0806] were found in Pasteurella dagmatis ([GenBank:EEX49428], 65% identity, E-value = 5e-113), H. influenzae ([GenBank:NTHI0936], 59% identity, E-value = 1e-101), and several species of Streptococci (e.g., [GenBank:SUB1393] and [GenBank:SPy_1634], 50-52% identity, E-values = 2e-86 to 1e-83). In strain 2336, [GenBank:HSM_0817] encoded a putative OmpR family response regulator whereas [GenBank:HSM_1124] encoded a putative ArcA family response regulator. These response regulators have very low homology to each other (33% identity, E-value = 7e-34). However, their amino-termini contained two conserved aspartic acid residues (Asp-11 and Asp-56) that may serve as phosphate acceptors and their carboxy-termini contained several conserved residues that may be involved in DNA binding (data not shown). Strains 2336 and 129Pt also contained genes that encoded putative histidine kinases ([GenBank:HSM_0483] and [GenBank:HS_1510], [GenBank:HSM_0727] and [GenBank:HS_0402], and [GenBank:HSM_1378 and [GenBank:HS_0900]). Whereas [GenBank:HSM_1378] and [GenBank:HS_0900] had an ORF downstream that encodes a putative cognate response regulator ([GenBank:HSM_1379] and [GenBank:HS_0901]), [GenBank:HSM_0727] and [GenBank:HS_0402] had an ORF downstream that encodes a transcriptional regulator with an N-terminal XRE-type helix-turn-helix domain ([GenBank:HSM_0728] and [GenBank:HS_0403]). Strain 2336 contained an additional gene that encoded a putative histidine kinase ([GenBank:HSM_0824]) that had a homolog in P. multocida strain Pm70 ([GenBank:PM1380], 58% identity, E-value = 0), but not in strain 129Pt.
Identification of genes encoding restriction-modification enzymes
Type I and Type II restriction-modification (RM) systems comprise the most frequently encountered DNA modifying enzymes among eubacteria. Strain 2336 contained ORFs that putatively encode proteins of the Type I RM system ([GenBank:HSM_1615], hsdR, site-specific deoxyribonuclease; [GenBank:HSM_1617], hsdS, DNA specificity subunit; [GenBank:HSM_1619], hsdM, DNA methylation subunit). Although strain 129Pt lacked hsdM, it contained a full-length hsdR ([GenBank:HS_0559]) and two truncated hsdS ([GenBank:HS_0554] and [GenBank:HS_0556]) orthologs. Strain 2336 also contained an operon putatively encoding enzymes of the Type II RM system ([GenBank:HSM_0801], M.hsoI, DNA-cytosine methyltransferase; [GenBank:HSM_0802], R.hsoI, restriction endonuclease). At the predicted protein level, M.HsoI has 68% identity to M.HinP1I ([GenBank:AAW33810], E-value = 2e-121) and R.HsoI has 72% identity to R.HinP1I ([GenBank:AAW33811], E-value = 3e-105) of H. influenzae P1. Strain 129Pt lacked M.hsoI and R.hsoI but contained an operon putatively encoding a Bcg I-like RM system ([GenBank:HS_0430], site-specific DNA-methyltransferase; [GenBank:HS_0431], restriction enzyme; [GenBank:HS_0432], restriction enzyme; [GenBank:HS_0433], hypothetical protein) that is absent in strain 2336.
Identification of other genes relevant to metabolism and survival
Both H. somni strains contained ORFs encoding proteins putatively involved in the transport and metabolism of fucose ([GenBank:HSM_0580] and [GenBank:HS_1451] to [GenBank:HSM_0585] and [GenBank:HS_1446]), mannose ([GenBank:HSM_0956] and [GenBank:HS_0605] to [GenBank:HSM_0960] and [GenBank:HS_0609]), xylose ([GenBank:HSM_0933] and[GenBank:HS_0584] to [GenBank:HSM_0937] and [GenBank:HS_0588]), galactitol ([GenBank:HSM_1030] and [GenBank:HS_1146] to [GenBank:HSM_1036] and [GenBank:/HS_1140]), mannitol ([GenBank:HSM_0825] and [GenBank:HS_1252] to [GenBank:HSM_0827] and [GenBank:HS_1250]), and myo-inositol ([GenBank:HSM_0425] and [GenBank:HS_1586] to [GenBank:HSM_0435] and [GenBank:HS_1576]). Furthermore, strain 2336 contained a locus with three ORFs encoding proteins putatively involved in galactose utilization (galTKM; [GenBank:HSM_0108] to [GenBank:HSM_0110]), whereas strain 129Pt contained only galM ([GenBank:HS_0236]) and truncated galK ([GenBank:HS_0235]) orthologs. However, both strains have galE, and strain 129Pt has galactose in its LOS (38), indicating that it can metabolize galactose.
Strain 2336-specific genes encoding proteins putatively involved in carbohydrate metabolism included [GenBank:HSM_0818] (ribokinase-like domain-containing protein), [GenBank:HSM_0819] (ketose-bisphosphate aldolase), [GenBank:HSM_0820] (sugar-related regulatory protein), [GenBank:HSM_0821] (sugar binding protein), [GenBank:HSM_0822] (monosaccharide-transporting ATPase), and [GenBank:HSM_0823] (ABC type sugar transporter). Strain 2336 also contained ORFs encoding a cyclase family protein and a membrane-spanning protein ([GenBank:HSM_1907] and [GenBank:HSM_1908]). These ORFs appeared to form an operon with ORFs encoding proteins putatively involved in butyrate/pyruvate metabolism ([GenBank:HSM_1903] to [GenBank:HSM_1906]). Furthermore, strain 2336 contained an ORF that encoded an oligopeptide permease ABC transporter homolog ([GenBank:HSM_0695]). The putative Opp protein of strain 2336 was related to OppB proteins involved in peptide transport in E. coli and Bacillus subtilis (e.g., [GenBank:AAC74326] and [GenBank:OPPB_BACSU], respectively, 46-76%% identity, E-values = 2e-138 to 1e-72).
GI I of strain 129Pt contained a gene encoding a putative phosphotransferase system protein ([GenBank:HS_0437], cellobiose-specific IIC component) whose homologs are found in M. haemolytica ([GenBank:MHA_2539], 94% identity, E-value = 0) and H. parasuis ([GenBank:HAPS_1535], 86% identity, E-value = 0). Strain 129Pt also contained a gene encoding a virulence-associated protein E ([GenBank:HS_0427]). This protein was related to the VirE proteins encoded by genes found on Staphylococcus phage phi2958PVL ([GenBank:BAG74398], 31% identity, E-value = 4e-36) and Enterococcus phage phiFL4A ([GenBank:ACZ64175], 33% identity, E-value = 3e-31). Strain 129Pt contained genes ([GenBank:HS_0631] and [GenBank:0632]) encoding a putative serine/threonine protein kinase-phosphatase pair. Strain 129Pt also contained loci encoding proteins putatively involved in thiamine biosynthesis (thiMDE; [GenBank:HS_0090-0092]), tellurite resistance (terDEZ; [GenBank:HS_0633-0635]), and lysine degradation (cadBA; [GenBank:HS_1006 and [GenBank:HS_1007]).
Genetic events such as deletions, duplications, insertions, and inversions are relatively common in bacterial chromosomes as a result of bacteriophage infection, integration and excision of plasmids, transpositions, and/or replication-mediated translocations [40, 41]. In addition, different prophages embedded within a single chromosome can contain similar genes encoding integration and structural functions, and it is not uncommon for these genes to undergo homologous recombination. One of the consequences of such homologous recombination is the rearrangement of the host chromosome . These events are known to be the precursors of evolution and can bring about a significant change in the number, linear order, and orientation of genes on the circular chromosomes of different strains/species of closely related bacteria .
The presence of PRs in the chromosomes of strains 2336 and 129Pt was a notable feature since the number and diversity of genes associated with these PRs far exceeded those described in H. influenzae strains Rd KW20 and 86-028NP [44, 45]. The difference in the size of the chromosomes of strains 2336 and 129Pt was partly due to PRs and associated genes. Similar observations have been made in other bacteria wherein prophage-associated sequences constitute a large portion of strain-specific DNA . Although one of the functions of RM systems is to afford protection against bacteriophage attack (the "cellular defense hypothesis"), it is interesting to note that both strains contain several prophage-like sequences despite the presence of genes encoding putative RM systems in their chromosomes. The lack of ORFs encoding HsdM, M.HsoI, and R.HsoI in strain 129Pt indicates that these systems are not absolutely essential for cell survival. Their absence may also partially explain the relative ease with which this strain can be transformed in the laboratory.
Biosynthesis of polysaccharides requires a multitude of GTs, which catalyze the transfer of sugars from an activated donor to an acceptor molecule and are usually specific for the glycosidic linkages created . Intra and interspecies divergence of genes encoding GTs are not uncommon. Phase-variable LOS is an important virulence factor of pathogenic strains of H. somni. Phase-variation of H. somni LOS has been shown to be due to the presence of SSRs in genes that encode GTs and enzymes involved in assembling non-glycose LOS components such as phosphorylcholine [5, 7, 37]. The genes lob1, lob2AB, and lob2D contain SSRs either just before the start codons or within the open reading frame . In addition, [GenBank:HSM_0148], [GenBank:HSM_0164], [GenBank:HSM_0975], and [GenBank:HSM_1552] also contain SSRs that may be responsible for LOS phase variation, but require additional experimental investigation. Most H. somni strains also produce a biofilm-associated EPS consisting primarily of mannose and galactose . Although characterization of some of the genes involved in the biosynthesis and/or modification of H. somni LOS/EPS has been determined, the identification of several more has been facilitated by comparative genome analyses [12, 13, 37]. It is likely that some of the observed variations among genes encoding GTs in H. somni and other Pasteurellaceae members is due to recombination events and/or selective pressure. Furthermore, variation in the composition and structure of the LOSs of strains 2336 and 129Pt may, in part, be due to different GT genes they have acquired or lost. In view of this, the role of new GT genes putatively involved in LOS biosynthesis and phase variation that have been identified in this study needs to be investigated.
Several types of two-partner secretion (TPS) pathways have been identified and characterized in Gram-negative bacteria . Filamentous hemagglutinins (Fha), consisting of a membrane-anchored protein (FhaC), which is involved in the activation/secretion of the cognate hemagglutinin/adhesin (FhaB), are prototypes of two-partner virulence systems. Homologs of FhaB and FhaC that possibly play a role in pathogenesis have been found in several members of the genera Bordetella, Haemophilus, Proteus, and Pasteurella, [49–53]. Among the 4 loci containing fha homologs in strain 2336, locus II appeared to be an acquisition mediated by a transposon and locus I appeared to be an acquisition due to homologous recombination. It appears that strain 129Pt has lost an fhaB homolog due to bacteriophage excision (Locus IV). It is possible that fhaB homologs in locus I of strain 2336 are paralogs, as are the fhaB homologs in locus IV. Together, these genes represent a large collection of fhaB and fhaC homologs in a single genome. The presence of multiple fhaB homologs in strain 2336 may, in part, be responsible for the serum resistance of this strain, in contrast to strain 129Pt, which does not contain full length fhaB homologs and is serum-sensitive . Structural and functional studies of N-terminal fragments of FHA from B. pertussis and HpmA from P. mirabilis have indicated that the proteins form a right-handed parallel β-helix [52–54]. Several residues that mediate the interaction of B. pertusis FHA with its cognate FhaC, and facilitate secretion, have also been identified . Although H. somni FhaB homologs have some features in common with FHA of B. pertussis and HpmA of P. mirabilis, a distinct region is the direct repeat 2 Fic domain in the FhaB homolog in locus III of strain 2336 that has been shown to induce cytotoxicity in human HeLa cells, bovine turbinate cells, and bovine alveolar type 2 cells [55, 56]. However, the secretion determinants of these proteins and the role of FhaC proteins in their secretion remain unknown.
Subtilases ([GenBank:COG1404]; subtilisin-like serine proteases) are a large superfamily of functionally diverse endo- and exo-peptidases that occur in prokaryotes and eukaryotes . Bacterial subtilisins may have a role in pathogenesis besides facilitating protein degradation and nutrient acquisition . However, subtilisin-like serine proteases from members of the Pasteurellaceae have not been characterized previously. The presence of a gene encoding a putative subtilase whose homologs were not found in other members of the Pasteurellaceae was yet another example of HGT in strain 2336. Although H. ducreyi strain 35000HP contains genes ([GenBank:HD1094] and [GenBank:HD1278]) encoding serine proteases that belong to the D-H-S family, they are unrelated to each other and to [GenBank:HSM_188]. In Agrobacterium tumefaciens, genes encoding AAA-ATPase and subtilisin-like serine protease have been shown to be functionally related and this pair has been proposed to constitute a toxin-antitoxin system that contributes to stability of plasmid pTiC58 . A conjugative megaplasmid-encoded subtilase has been shown to be a virulence factor in E. coli and it has been suggested this toxin activates a V-ATPase in Vero cells [60, 61]. Whether the ATPase-subtilisin pair identified in this study is transcriptionally and functionally coupled, and whether the protease gene contributes to the pathogenicity of strain 2336, has yet to be determined.
The ability to acquire and metabolize iron is an important determinant of bacterial survival and adaptability. In some bacteria, genes that facilitate iron uptake have been shown to be acquired by horizontal transfer . Several members of the Pasteurellaceae possess special outer membrane protein (OMP) receptors consisting of two unrelated transferrin-binding proteins, TbpA and TbpB, which facilitate acquisition of transferrin-bound iron from their hosts . H. somni strain 649 has a TbpA-TbpB receptor system that acquires iron only from bovine transferrin, and a second, probably redundant, TbpA2 receptor that can acquire iron from bovine, caprine, or ovine transferrins . From genomic analyses, it is apparent that H. somni strains possess multiple genes for the acquisition of iron. Horizontal transfer and clustering of genes related to iron metabolism is indicative of enrichment and adaptation of pathogenic H. somni to different niches within its natural host. Furthermore, products of one or more of these genes may facilitate binding to transferrins/lactoferrins of different host species and such a gene repertoire could enhance the ability of this bacterium to survive in a variety of ruminants. Studies using a mouse model have suggested the role of bovine transferrin and lactoferrin in increasing the virulence of strain 2336 .
Many bacterial pathogens contain surface proteins that facilitate adhesion to and/or invasion of the host mucosal barriers [66, 67]. Some of these proteins may also be involved in bacterial aggregation to form biofilms and their evasion of the host's innate immune system [67, 68]. At least three major categories of bacterial adherence proteins have been identified, two of which are hair-like structures called pili and non-pilus-associated proteins called adhesins. The plasmid-encoded YadA is a prototype non-pilus-associated protein that has been well characterized . Mutation or deletion of the yadA homolog can reduce virulence in pathogenic bacteria [70, 71]. Large adhesins in some bacteria are associated with ABC transporters and may be involved in biofilm formation . A transposon mutagenesis approach has implicated several genes, including those encoding filamentous haemagglutinins, in H. somni biofilm formation . H. somni also contains genes putatively involved in adhesin synthesis/transport, pilus formation, and quorum sensing (e.g., luxS), but their role in facilitating biofilm formation remains to be investigated.
Veterinarians use several antibiotics to treat H. somni infections and feedlot cattle enterprises frequently rely on tetracyclines for prophylaxis as well as growth promotion [1, 74, 75]. Although not common, H. somni resistance to tetracycline has been reported [76, 77]. Copper and zinc are often included in commercial cattle diets to achieve optimal growth and reproduction. Emergence of copper/zinc resistance in bacteria of animal origin has been documented and attributed to the excessive presence of these metals in livestock feed . Furthermore, the occurrence of genes related to metal and antibiotic resistance on integrative/conjugative elements and their horizontal co-transfer has been noted previously [79, 80]. In view of these observations, it was not surprising to find a GI containing genes putatively involved in copper, zinc, and tetracycline resistance in strain 2336.
Transcriptional regulators play crucial roles in bacterial functions and they have been classified into a number of families . The homologs of [GenBank:HSM_0806] (LysR, [NCBI:CLSK797597] cluster) in H. influenzae and P. dagmatis are associated with genes encoding proteins involved in fatty acid metabolism (e.g., acetyl-CoA acetyltransferase, 3-oxoacid CoA- transferase, and fatty acid transporters). Therefore, this cluster may represent a novel class of metabolic regulators within the LysR family. Most members of the [NCBI:PRK13756] cluster are involved in regulation of antibiotic resistance genes . Homologs of [GenBank:HSM_1734] and [GenBank:HSM_1735] among members of the Pasteurellaceae encode tetracycline resistance and are associated with mobile genetic elements [82–86]. Homologs of [GenBank:HSM_1736] and [GenBank:HSM_1737] in other bacteria are known to be horizontally transferred and may mediate resistance to antibiotics . Furthermore, homologs of [GenBank:HSM_1192] and [GenBank:HSM_1193] are predicted to be involved in multidrug resistance [88, 89]. In summary, it appears that strain 2336 contained at least three different systems related to antibiotic resistance. Although the functional role of these genes remains to be established, their similarity to metal/antibiotic resistance genes associated with mobile genetic elements in other members of the Pasteurellaceae is clinically significant.
From genome comparisons, it appears that there is no correlation between chromosome size and the number of tRNA genes (the genomes of H. somni 2336, H. somni 129Pt, H. influenzae 86-028NP, H. ducreyi 35000HP, and P. multocida Pm70 contain 49, 49, 58, 46, and 57 tRNA genes, respectively). Whether the lower number of tRNA genes found in H. somni strains is due to disruptive integration of bacteriophages into tRNA genes (as in 'bacteriophage disruption of tRNA genes in Lactobacillus johnsonii' ) or is a result of compensatory gene loss in lieu of acquisition of new genes (as in 'genome reduction in pathogenic and symbiotic bacteria' ) is unknown. Nevertheless, comparison of the chromosomes of strains 129Pt and 2336 bolsters the proposition that prophages and transposons have played a major role in creating genomic diversity and phenotypic variability in the two strains. It is also apparent that strains 2336 and 129Pt have independently and intermittently acquired and lost genes since their divergence from a common ancestor, and that the net gain in strain 129Pt is less than the net gain in strain 2336.
H. somni strain 2336 contains a larger chromosome when compared to other Haemophilus and Histophilus strains whose genome sequences are available. Several regions that resemble the pathogenicity islands of other virulent bacteria are present in strain 2336. There is evidence to suggest that most of these regions were acquired by HGT mechanisms, whereas similar regions were not found in the commensal strain 129Pt. Although previous studies have discovered the genetic basis for some of the phenotypic dissimilarities between strains 2336 and 129Pt, complete genome sequence analyses have provided a comprehensive account of innate and acquired genetic traits. Furthermore, comparisons of the genomes of strains 2336 and 129Pt have contributed to our understanding of the biology and pathogenic evolution of these bacteria. The post-genomic era for H. somni poses new challenges and opportunities in terms of functional characterization of genes and deciphering their roles in colonization, survival, and pathogenesis. Continued analyses of the genomes of H. somni strains and comparison with newly sequenced genomes of other bacteria should enhance the current knowledge on virulence mechanisms. Nevertheless, the results from this study are expected to facilitate the development of improved diagnostic tests for and vaccines against H. somni.
This work was performed under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program and by the University of California, Lawrence Livermore National Laboratory under Contract No. DE-AC02-05CH11231, Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344, and Los Alamos National Laboratory under contract No. DE-AC02-06NA25396. This work was supported by U.S. Department of Agriculture (Cooperative State Research, Education, and Extension Service, Initiative for Future Agriculture and Food Systems) grants 2001-52100-11314, 2003-35204-13637, 2007-35204-18338 to TJI, the U.S. Department of Energy under contract no. W-7405-ENG-36, and funds from the Virginia Agricultural Experiment Station. We thank Gretchen Berg, Dr. Michael Howard, and Dr. Shaadi Elswaifi for technical assistance.
- Inzana TJ: The Haemophilus somnus complex. Current Veterinary Therapy 4. Edited by: Howard JL, Smith RA. 1999, Philadelphia, PA: W.B. Saunders Co., 358-361. 4
- Corbeil LB: Histophilus somni host-parasite relationships. Anim Health Res Rev. 2007, 8: 151-160. 10.1017/S1466252307001417.PubMedView Article
- Widders PR, Dorrance LA, Yarnall M, Corbeil LB: Immunoglobulin-binding activity among pathogenic and carrier isolates of Haemophilus somnus. Infect Immun. 1989, 57: 639-642.PubMed CentralPubMed
- Corbeil L, Gogolewski R, Stephens L, Inzana T: Haemophilus somnus: antigen analysis and immune responses. Haemophilus, Actinobacillus, and Pasteurella\. Edited by: Donachie W, Lainson F, Hodgson J. 1995, New York: Plenum Press, 63-73.View Article
- Wu Y, McQuiston JH, Cox A, Pack TD, Inzana TJ: Molecular cloning and mutagenesis of a DNA locus involved in lipooligosaccharide biosynthesis in Haemophilus somnus. Infect Immun. 2000, 68: 310-319. 10.1128/IAI.68.1.310-319.2000.PubMed CentralPubMedView Article
- Cole SP, Guiney DG, Corbeil LB: Molecular analysis of a gene encoding a serum-resistance-associated 76 kDa surface antigen of Haemophilus somnus. J Gen Microbiol. 1993, 139: 2135-2143.PubMedView Article
- McQuiston JH, McQuiston JR, Cox AD, Wu Y, Boyle SM, Inzana TJ: Characterization of a DNA region containing 5'-(CAAT)(n)-3' DNA sequences involved in lipooligosaccharide biosynthesis in Haemophilus somnus. Microb Pathog. 2000, 28: 301-312. 10.1006/mpat.1999.0351.PubMedView Article
- Inzana TJ, Gogolewski RP, Corbeil LB: Phenotypic phase variation in Haemophilus somnus lipooligosaccharide during bovine pneumonia and after in vitro passage. Infect Immun. 1992, 60: 2943-2951.PubMed CentralPubMed
- Inzana TJ, Glindemann G, Cox AD, Wakarchuk W, Howard MD: Incorporation of N-acetylneuraminic acid into Haemophilus somnus lipooligosaccharide (LOS): enhancement of resistance to serum and reduction of LOS antibody binding. Infect Immun. 2002, 70: 4870-4879. 10.1128/IAI.70.9.4870-4879.2002.PubMed CentralPubMedView Article
- Corbeil LB, Bastida-Corcuera FD, Beveridge TJ: Haemophilus somnus immunoglobulin binding proteins and surface fibrils. Infect Immun. 1997, 65: 4250-4257.PubMed CentralPubMed
- Fraser CM, Eisen J, Fleischmann RD, Ketchum KA, Peterson S: Comparative genomics and understanding of microbial biology. Emerg Infect Dis. 2000, 6: 505-512. 10.3201/eid0605.000510.PubMed CentralPubMedView Article
- Challacombe JF, Inzana TJ: Comparative genomics of Pasteurellaceae. Pasteurellaceae: Biology, Genomics, and Molecular Aspects\. Edited by: Kuhnert P, Christensen H. 2008, Norfolk, UK: Caister Academic Press, 53-77.
- Challacombe JF, Duncan AJ, Brettin TS, Bruce D, Chertkov O, Detter JC, Han CS, Misra M, Richardson P, Tapia R, Thayer N, Xie G, Inzana TJ: Complete genome sequence of Haemophilus somnus (Histophilus somni) strain 129Pt and comparison to Haemophilus ducreyi 35000HP and Haemophilus influenzae Rd. J Bacteriol. 2007, 189: 1890-1898. 10.1128/JB.01422-06.PubMed CentralPubMedView Article
- Hacker J, Blum-Oehler G, Muhldorfer I, Tschape H: Pathogenicity islands of virulent bacteria: structure, function and impact on microbial evolution. Mol Microbiol. 1997, 23: 1089-1097. 10.1046/j.1365-2958.1997.3101672.x.PubMedView Article
- Olson M: When less is more: gene loss as an engine of evolutionary change. Am J Hum Genet. 1999, 64: 18-23. 10.1086/302219.PubMed CentralPubMedView Article
- St Michael F, Li J, Howard MD, Duncan AJ, Inzana TJ, Cox AD: Structural analysis of the oligosaccharide of Histophilus somni (Haemophilus somnus) strain 2336 and identification of several lipooligosaccharide biosynthesis gene homologues. Carbohydr Res. 2005, 340: 665-672. 10.1016/j.carres.2004.12.029.PubMedView Article
- Corbeil LB, Blau K, Prieur DJ, Ward AC: Serum susceptibility of Haemophilus somnus from bovine clinical cases and carriers. J Clin Microbiol. 1985, 22 (2): 192-198.PubMed CentralPubMed
- Lander ES, Waterman MS: Genomic mapping by fingerprinting random clones: a mathematical analysis. Genomics. 1988, 2: 231-239. 10.1016/0888-7543(88)90007-9.PubMedView Article
- Chissoe SL, Wang YF, Clifton SW, Ma N, Sun HJ, Lobsinger JS, Kenton SM, White JD, Roe BA: Strategies for rapid and accurate DNA sequencing. Methods: Companion Methods Enzymol. 1991, 3: 55-65. 10.1016/S1046-2023(05)80164-1.View Article
- Munson RS, Harrison A, Gillaspy A, Ray WC, Carson M, Armbruster D, Gipson J, Gipson M, Johnson L, Lewis L, Dyer DW, Bakaletz LO: Partial analysis of the genomes of two nontypeable Haemophilus influenzae otitis media isolates. Infect Immun. 2004, 72: 3002-3010. 10.1128/IAI.72.5.3002-3010.2004.PubMed CentralPubMedView Article
- Karlyshev AV, Pallen MJ, Wren BW: Single-primer PCR procedure for rapid identification of transposon insertion sites. Biotechniques. 2000, 28: 1078-PubMed
- Gordon D, Desmarais C, Green P: Automated finishing with autofinish. Genome Res. 2001, 11: 614-625. 10.1101/gr.171401.PubMed CentralPubMedView Article
- Han CS, Chain P: Finishing repeat regions automatically with Dupfinisher. Proc 2006 Int Conf Bioinformatics Comput Biol\. Edited by: Arabnia HR, and H. Valafar. 2006, Las Vegas, NV.: CSREA Press, 141-146.
- Anonymous: Transposon-based strategies for efficient DNA sequencing and functional genomics (available online at http://www.epibio.com/forum_issue.asp. EPICENTRE Forum. 2004, 11: 22-23.
- Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, Eddy SR, Sonnhammer ELL, Bateman A: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, D247-D251. 34 Database
- Smith TF, Waterman MS: Identification of common molecular subsequences. J Mol Biol. 1981, 147: 195-197. 10.1016/0022-2836(81)90087-5.PubMedView Article
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008, 9: 75-10.1186/1471-2164-9-75.PubMed CentralPubMedView Article
- Darling AC, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14: 1394-1403. 10.1101/gr.2289704.PubMed CentralPubMedView Article
- Darling AE, Treangen TJ, Messeguer X, Perna NT: Analyzing patterns of microbial evolution using the mauve genome alignment system. Methods Mol Biol. 2007, 396: 135-152. 10.1007/978-1-59745-515-2_10.PubMedView Article
- Delcher AL, Phillippy A, Carlton J, Salzberg SL: Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res. 2002, 30: 2478-2483. 10.1093/nar/30.11.2478.PubMed CentralPubMedView Article
- Lima-Mendez G, Van Helden J, Toussaint A, Leplae R: Prophinder: a computational tool for prophage prediction in prokaryotic genomes. Bioinform. 2008, 24: 863-865. 10.1093/bioinformatics/btn043.View Article
- Langille MG, Hsiao WW, Brinkman FS: Detecting genomic islands using bioinformatics approaches. Nature Rev Microbiol. 2010, 8: 373-382. 10.1038/nrmicro2350.View Article
- Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Grechkin Y, Ratner A, Anderson I, Lykidis A, Mavromatis K, Ivanova NN, Kyrpides NC: The integrated microbial genomes system: an expanding comparative analysis resource. Nucleic Acids Res. 2010, D382-D390. 38 Database
- Siddaramappa S, Duncan AJ, Brettin T, Inzana TJ: Comparative analyses of two cryptic plasmids from Haemophilus somnus (Histophilus somni). Plasmid. 2006, 55: 227-234. 10.1016/j.plasmid.2005.11.004.PubMedView Article
- Feschotte C: Merlin, a new superfamily of DNA transposons identified in diverse animal genomes and related to bacterial IS1016 insertion sequences. Mol Biol Evol. 2004, 21: 1769-1780. 10.1093/molbev/msh188.PubMedView Article
- Siguier P, Gagnevin L, Chandler M: The new IS1595 family, its relation to IS1 and the frontier between insertion sequences and transposons. Res Microbiol. 2009, 160: 232-241. 10.1016/j.resmic.2009.02.003.PubMedView Article
- Elswaifi SF, St Michael F, Sreenivas A, Cox A, Carman GM, Inzana TJ: Molecular characterization of phosphorylcholine expression on the lipooligosaccharide of Histophilus somni. Microb Pathog. 2009, 47: 223-230. 10.1016/j.micpath.2009.08.001.PubMed CentralPubMedView Article
- St Michael F, Howard MD, Li J, Duncan AJ, Inzana TJ, Cox AD: Structural analysis of the lipooligosaccharide from the commensal Haemophilus somnus genome strain 129Pt. Carbohydr Res. 2004, 339: 529-535. 10.1016/j.carres.2003.10.032.PubMedView Article
- Ward CK, Lumbley SR, Latimer JL, Cope LD, Hansen EJ: Haemophilus ducreyi secretes a filamentous hemagglutinin-like protein. J Bacteriol. 1998, 180: 6013-6022.PubMed CentralPubMed
- Tillier ER, Collins RA: Genome rearrangement by replication-directed translocation. Nat Genet. 2000, 26: 195-197. 10.1038/79918.PubMedView Article
- Brussow H, Canchaya C, Hardt WD: Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol Mol Biol Rev. 2004, 68: 560-602. 10.1128/MMBR.68.3.560-602.2004.PubMed CentralPubMedView Article
- Canchaya C, Fournous G, Brussow H: The impact of prophages on bacterial chromosomes. Mol Microbiol. 2004, 53: 9-18. 10.1111/j.1365-2958.2004.04113.x.PubMedView Article
- Moret BME, Bader DA, Warnow T: High-performance algorithm engineering for computational phylogenetics. J Supercomput. 2002, 22: 99-111. 10.1023/A:1014362705613.View Article
- Harrison A, Dyer DW, Gillaspy A, Ray WC, Mungur R, Carson MB, Zhong H, Gipson J, Gipson M, Johnson LS, Lewis L, Bakaletz LO, Munson RS: Genomic sequence of an otitis media isolate of nontypeable Haemophilus influenzae: comparative study with H. influenzae serotype d, strain KW20. J Bacteriol. 2005, 187: 4627-4636. 10.1128/JB.187.13.4627-4636.2005.PubMed CentralPubMedView Article
- Fleischmann RD, Adams MD, White O, Clayton RA, Kirkness EF, Kerlavage AR, Bult CJ, Tomb JF, Dougherty BA, Merrick JM, Sutton G, FitzHugh W, Fields C, Gocayne JD, Scott J, Shirley R, Liu L-I, Glodek A, Kelly JM, Weidman JF, Phillips CA, Spriggs T, Hadblom E, Cotton MD, Utterback TR, Hanna MC, Nguyen DT, Saudek DM, Brandon RC, Fine LD, et al: Whole-genome random sequencing and assembly of Haemophilus influenzae Rd. Science. 1995, 269: 496-512. 10.1126/science.7542800.PubMedView Article
- Lundborg M, Modhukur V, Widmalm G: Glycosyltransferase functions of E. coli O-antigens. Glycobiol. 20: 366-368.
- Inzana TJ, Swords WE, Sandal I, Siddaramappa S: Lipopolysaccharides, biofilms and quorum sensing in Pasteurellaceae. Pasteurellaceae: Biology, Genomics and Molecular Aspects\. Edited by: Kuhnert P, Christensen H. 2008, Norfolk, UK: Caister Academic Press, 177-195.
- Jacob-Dubuisson F, Locht C, Antoine R: Two-partner secretion in Gram-negative bacteria: a thrifty, specific pathway for large virulence proteins. Mol Microbiol. 2001, 40: 306-313. 10.1046/j.1365-2958.2001.02278.x.PubMedView Article
- Fuller TE, Kennedy MJ, Lowery DE: Identification of Pasteurella multocida virulence genes in a septicemic mouse model using signature-tagged mutagenesis. Microb Pathog. 2000, 29: 25-38. 10.1006/mpat.2000.0365.PubMedView Article
- Ward CK, Mock JR, Hansen EJ: The LspB protein is involved in the secretion of the LspA1 and LspA2 proteins by Haemophilus ducreyi. Infect Immun. 2004, 72: 1874-1884. 10.1128/IAI.72.4.1874-1884.2004.PubMed CentralPubMedView Article
- May BJ, Zhang Q, Li LL, Paustian ML, Whittam TS, Kapur V: Complete genomic sequence of Pasteurella multocida, Pm70. Proc Natl Acad Sci USA. 2001, 98: 3460-3465. 10.1073/pnas.051634598.PubMed CentralPubMedView Article
- Clantin B, Hodak H, Willery E, Locht C, Jacob-Dubuisson F, Villeret V: The crystal structure of filamentous hemagglutinin secretion domain and its implications for the two-partner secretion pathway. Proc Natl Acad Sci USA. 2004, 101: 6194-6199. 10.1073/pnas.0400291101.PubMed CentralPubMedView Article
- Weaver TM, Hocking JM, Bailey LJ, Wawrzyn GT, Howard DR, Sikkink LA, Ramirez-Alvarado M, Thompson JR: Structural and functional studies of truncated hemolysin A from Proteus mirabilis. J Biol Chem. 2009, 284: 22297-22309. 10.1074/jbc.M109.014431.PubMed CentralPubMedView Article
- Hodak H, Clantin B, Willery E, Villeret V, Locht C, Jacob-Dubuisson F: Secretion signal of the filamentous haemagglutinin, a model two-partner secretion substrate. Mol Microbiol. 2006, 61: 368-382. 10.1111/j.1365-2958.2006.05242.x.PubMedView Article
- Worby CA, Mattoo S, Kruger RP, Corbeil LB, Koller A, Mendez JC, Zekarias B, Lazar C, Dixon JE: The fic domain: regulation of cell signaling by adenylylation. Mol Cell. 2009, 34: 93-103. 10.1016/j.molcel.2009.03.008.PubMed CentralPubMedView Article
- Zekarias B, Mattoo S, Worby C, Lehmann J, Rosenbusch RF, Corbeil LB: Histophilus somni IbpA DR2/Fic in virulence and immunoprotection at the natural host alveolar epithelial barrier. Infect Immun. 2010, 78: 1850-1858. 10.1128/IAI.01277-09.PubMed CentralPubMedView Article
- Siezen RJ, Leunissen JA: Subtilases: the superfamily of subtilisin-like serine proteases. Protein Sci. 1997, 6: 501-523.PubMed CentralPubMedView Article
- Siezen RJ, Renckens B, Boekhorst J: Evolution of prokaryotic subtilases: genome-wide analysis reveals novel subfamilies with different catalytic residues. Proteins. 2007, 67: 681-694. 10.1002/prot.21290.PubMedView Article
- Yamamoto S, Kiyokawa K, Tanaka K, Moriguchi K, Suzuki K: Novel toxin-antitoxin system composed of serine protease and AAA-ATPase homologues determines the high level of stability and incompatibility of the tumor-inducing plasmid pTiC58. J Bacteriol. 2009, 191: 4656-4666. 10.1128/JB.00124-09.PubMed CentralPubMedView Article
- Paton AW, Srimanote P, Talbot UM, Wang H, Paton JC: A new family of potent AB(5) cytotoxins produced by Shiga toxigenic Escherichia coli. J Exp Med. 2004, 200: 35-46. 10.1084/jem.20040392.PubMed CentralPubMedView Article
- Yahiro K, Morinaga N, Satoh M, Matsuura G, Tomonaga T, Nomura F, Moss J, Noda M: Identification and characterization of receptors for vacuolating activity of subtilase cytotoxin. Mol Microbiol. 2006, 62: 480-490. 10.1111/j.1365-2958.2006.05379.x.PubMedView Article
- Dobrindt U, Hochhut B, Hentschel U, Hacker J: Genomic islands in pathogenic and environmental microorganisms. Nat Rev Microbiol. 2004, 2: 414-424. 10.1038/nrmicro884.PubMedView Article
- Cornelissen CN: Transferrin-iron uptake by Gram-negative bacteria. Front Biosci. 2003, 8: d836-847. 10.2741/1076.PubMedView Article
- Ekins A, Bahrami F, Sijercic A, Maret D, Niven DF: Haemophilus somnus possesses two systems for acquisition of transferrin-bound iron. J Bacteriol. 2004, 186: 4407-4411. 10.1128/JB.186.13.4407-4411.2004.PubMed CentralPubMedView Article
- Geertsema RS, Kimball RA, Corbeil LB: Bovine plasma proteins increase virulence of Haemophilus somnus in mice. Microb Pathog. 2007, 42: 22-28. 10.1016/j.micpath.2006.10.001.PubMedView Article
- Hultgren SJ, Abraham S, Caparon M, Falk P, St Geme JW, Normark S: Pilus and nonpilus bacterial adhesins: assembly and function in cell recognition. Cell. 1993, 73: 887-901. 10.1016/0092-8674(93)90269-V.PubMedView Article
- Kline KA, Falker S, Dahlberg S, Normark S, Henriques-Normark B: Bacterial adhesins in host-microbe interactions. Cell Host Microbe. 2009, 5: 580-592. 10.1016/j.chom.2009.05.011.PubMedView Article
- Jurcisek JA, Bakaletz LO: Biofilms formed by nontypeable Haemophilus influenzae in vivo contain both double-stranded DNA and type IV pilin protein. J Bacteriol. 2007, 189: 3868-3875. 10.1128/JB.01935-06.PubMed CentralPubMedView Article
- Hoiczyk E, Roggenkamp A, Reichenbecher M, Lupas A, Heesemann J: Structure and sequence analysis of Yersinia YadA and Moraxella UspAs reveal a novel class of adhesins. EMBO J. 2000, 19: 5989-5999. 10.1093/emboj/19.22.5989.PubMed CentralPubMedView Article
- Heise T, Dersch P: Identification of a domain in Yersinia virulence factor YadA that is crucial for extracellular matrix-specific cell adhesion and uptake. Proc Natl Acad Sci USA. 2006, 103: 3375-3380. 10.1073/pnas.0507749103.PubMed CentralPubMedView Article
- Lawrenz MB, Lenz JD, Miller VL: A novel autotransporter adhesin is required for efficient colonization during bubonic plague. Infect Immun. 2009, 77: 317-326. 10.1128/IAI.01206-08.PubMed CentralPubMedView Article
- Gerlach RG, Hensel M: Protein secretion systems and adhesins: the molecular armory of Gram-negative pathogens. Int J Med Microbiol. 2007, 297: 401-415. 10.1016/j.ijmm.2007.03.017.PubMedView Article
- Sandal I, Shao JQ, Annadata S, Apicella MA, Boye M, Jensen TK, Saunders GK, Inzana TJ: Histophilus somni biofilm formation in cardiopulmonary tissue of the bovine host following respiratory challenge. Microbes Infect. 2009, 11: 254-263. 10.1016/j.micinf.2008.11.011.PubMedView Article
- Smith MS, Yang RK, Knapp CW, Niu Y, Peak N, Hanfelt MM, Galland JC, Graham DW: Quantification of tetracycline resistance genes in feedlot lagoons by real-time PCR. Appl Environ Microbiol. 2004, 70: 7372-7377. 10.1128/AEM.70.12.7372-7377.2004.PubMed CentralPubMedView Article
- Van Donkersgoed J, Janzen ED, Potter AA, Harland RJ: The occurrence of Haemophilus somnus in feedlot calves and its control by postarrival prophylactic mass medication. Can Vet J. 1994, 35: 573-580.PubMed CentralPubMed
- Watts JL, Yancey RJ, Salmon SA, Case CA: A 4-year survey of antimicrobial susceptibility trends for isolates from cattle with bovine respiratory disease in North America. J Clin Microbiol. 1994, 32: 725-731.PubMed CentralPubMed
- Mohd-Zain Z, Turner SL, Cerdeno-Tarraga AM, Lilley AK, Inzana TJ, Duncan AJ, Harding RM, Hood DW, Peto TE, Crook DW: Transferable antibiotic resistance elements in Haemophilus influenzae share a common evolutionary origin with a diverse family of syntenic genomic islands. J Bacteriol. 2004, 186: 8114-8122. 10.1128/JB.186.23.8114-8122.2004.PubMed CentralPubMedView Article
- Hasman H, Franke S, Rensing C: Resistance to metals used in agriculture production. Antimicrobial resistance in bacteria of animal origin\. Edited by: Aarestrup F. 2006, Washington, DC.: ASM Press
- Gilmour MW, Thomson NR, Sanders M, Parkhill J, Taylor DE: The complete nucleotide sequence of the resistance plasmid R478: defining the backbone components of incompatibility group H conjugative plasmids through comparative genomics. Plasmid. 2004, 52: 182-202. 10.1016/j.plasmid.2004.06.006.PubMedView Article
- Baker-Austin C, Wright MS, Stepanauskas R, McArthur JV: Co-selection of antibiotic and metal resistance. Trends Microbiol. 2006, 14: 176-182. 10.1016/j.tim.2006.02.006.PubMedView Article
- Ramos JL, Martinez-Bueno M, Molina-Henares AJ, Teran W, Watanabe K, Zhang X, Gallegos MT, Brennan R, Tobes R: The TetR family of transcriptional repressors. Microbiol Mol Biol Rev. 2005, 69: 326-356. 10.1128/MMBR.69.2.326-356.2005.PubMed CentralPubMedView Article
- Blanco M, Kadlec K, Gutierrez Martin CB, de la Fuente AJ, Schwarz S, Navas J: Nucleotide sequence and transfer properties of two novel types of Actinobacillus pleuropneumoniae plasmids carrying the tetracycline resistance gene tet(H). J Antimicrob Chemother. 2007, 60: 864-867. 10.1093/jac/dkm293.PubMedView Article
- Kehrenberg C, Werckenthin C, Schwarz S: Tn5706, a transposon-like element from Pasteurella multocida mediating tetracycline resistance. Antimicrob Agents Chemother. 1998, 42: 2116-2118.PubMed CentralPubMed
- Hansen LM, McMurry LM, Levy SB, Hirsh DC: A new tetracycline resistance determinant, Tet H, from Pasteurella multocida specifying active efflux of tetracycline. Antimicrob Agents Chemother. 1993, 37: 2699-2705.PubMed CentralPubMedView Article
- Kehrenberg C, Schwarz S: Identification of a truncated, but functionally active tet(H) tetracycline resistance gene in Pasteurella aerogenes and Pasteurella multocida. FEMS Microbiol Lett. 2000, 188: 191-195. 10.1111/j.1574-6968.2000.tb09192.x.PubMedView Article
- Kehrenberg C, Salmon SA, Watts JL, Schwarz S: Tetracycline resistance genes in isolates of Pasteurella multocida, Mannheimia haemolytica, Mannheimia glucosida and Mannheimia varigena from bovine and swine respiratory disease: intergeneric spread of the tet(H) plasmid pMHT1. J Antimicrob Chemother. 2001, 48: 631-640. 10.1093/jac/48.5.631.PubMedView Article
- Bay DC, Turner RJ: Diversity and evolution of the small multidrug resistance protein family. BMC Evol Biol. 2009, 9: 140-10.1186/1471-2148-9-140.PubMed CentralPubMedView Article
- Saier MH, Paulsen IT, Sliwinski MK, Pao SS, Skurray RA, Nikaido H: Evolutionary origins of multidrug and drug-specific efflux pumps in bacteria. FASEB J. 1998, 12: 265-274.PubMed
- Young J, Holland IB: ABC transporters: bacterial exporters-revisited five years on. Biochim Biophys Acta. 1999, 1461: 177-200. 10.1016/S0005-2736(99)00158-3.PubMedView Article
- Ventura M, Canchaya C, Pridmore D, Berger B, Brussow H: Integration and distribution of Lactobacillus johnsonii prophages. J Bacteriol. 2003, 185: 4603-4608. 10.1128/JB.185.15.4603-4608.2003.PubMed CentralPubMedView Article
- Hacker J, Hentschel U, Dobrindt U: Prokaryotic chromosomes and disease. Science. 2003, 301: 790-793. 10.1126/science.1086802.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.