Evolutionary stability of DNA uptake signal sequences in the Pasteurellaceae

Multiple interspecies recombination events documented by whole-genome sequencing in multidrug-resistant Haemophilus influenzae clinical isolates

Introduction

Haemophilus influenzae (Hi) was long known as an easy-to-treat bacterium, but increasing resistance against beta-lactams and other critically important antibiotics is now a growing concern. We describe here the whole-genome sequencing (WGS) analysis of three non-typeable Hi isolates received in 2018-2019 by the Belgian National Reference Centre (NRC) for Haemophilus influenzae, as they presented an unusual multi-resistant profile.

Methods

All three isolates were sequenced by WGS and mapped to the reference isolate Hi Rd KW20. Shorten uptake signal sequences (USSs) known to be associated with homologous recombination were sought in ftsI, murE and murF genes, and inner partial sequences were compared against the blast nucleotide database to look for similarity with other Haemophilus species. Their antimicrobial resistance (AMR) genotype was studied. Core-genome multilocus sequence typing (MLST) was performed on the NTHi database pubMLST to place our isolates in the actual worldwide epidemiology.

Results

The isolates also harboured interspecies recombination patterns in the murF-murE-ftsI region involved in cell wall synthesis. The three isolates were multidrug resistant and two of them were also resistant to amoxicillin-clavulanic acid and showed a reduced susceptibility to meropenem. All three isolates belonged to the MLST clonal complex (CC) 422, and WGS revealed that the three were very similar. They harboured mobile genetic elements (carrying blaTEM-1B, mefA and msrD genes associated with resistance), mutations in gyrA and parC linked to fluoroquinolone resistance as well as remodelling events in ompP2 that might be related to lower carbapenem susceptibility.

Conclusion

The Hi evolution towards antimicrobial multiresistance (AMR) is a complex and poorly understood phenomenon, although probably linked to a large degree to the presence of USSs and exchange within the family Pasteurellaceae. To better understand the respective roles of clonal expansion, horizontal gene transfers, spontaneous mutations and interspecies genetic rearrangements in shaping Hi AMR, both analysis of Hi communities over time within individuals and worldwide monitoring of non-typeable Hi causing infections should be conducted.

Basic Characterization of Natural Transformation in Avibacterium paragallinarum

Avibacterium paragallinarum is the pathogen involved in infectious coryza (IC), an acute infectious upper respiratory disease in chickens. The prevalence of IC has increased in China in recent years. There is a lack of reliable and effective procedures for gene manipulation, which has limited the research on the bacterial genetics and pathogenesis of A. paragallinarum. Natural transformation has been developed as a method of gene manipulation in Pasteurellaceae by the introduction of foreign genes or DNA fragments into bacterial cells, but there has been no report on natural transformation in A. paragallinarum. In this study, we analyzed the existence of homologous genetic factors and competence proteins underlying natural transformation in A. paragallinarum and established a method for transformation in it. Through bioinformatics analysis, we identified 16 homologs of Haemophilus influenzae competence proteins in A. paragallinarum. We found that the uptake signal sequence (USS) was overrepresented in the genome of A. paragallinarum (1,537 to 1,641 copies of the core sequence ACCGCACTT). We then constructed a plasmid, pEA-KU, that carries the USS and a plasmid, pEA-K, without the USS. These plasmids can be transferred via natural transformation into naturally competent strains of A. paragallinarum. Significantly, the plasmid that carries USS showed a higher transformation efficiency. In summary, our results demonstrate that A. paragallinarum has the ability to undergo natural transformation. These findings should prove to be a valuable tool for gene manipulation in A. paragallinarum. IMPORTANCE Natural transformation is an important mechanism for bacteria to acquire exogenous DNA molecules during the process of evolution. Additionally, it can also be used as a method to introduce foreign genes into bacteria under laboratory conditions. Natural transformation does not require equipment such as an electroporation apparatus. It is easy to perform and is similar to gene transfer under natural conditions. However, there have been no reports on natural transformation in Avibacterium paragallinarum. In this study, we analyzed the presence of homologous genetic factors and competence proteins underlying natural transformation in A. paragallinarum. Our results indicate that natural competence could be induced in A. paragallinarum serovars A, B, and C. Furthermore, the method that we established to transform plasmids into naturally competent A. paragallinarum strains was stable and efficient.

A cross-species genome-wide analysis of sequences similar to those involved in DNA uptake bias in the Pastuerellaceae and Neisseriaceae families of pathogenic bacteria

Acquiring new DNA allows the emergence of drug resistance in bacteria. Some Pasteurellaceae and Neisseriaceae species preferentially take up specific sequence tags. The study of such sequences is therefore relevant. They are over-represented in the genomes of the corresponding species. I found similar sequences to be present only in, but not in all, the genomes of the Pasteurellaceae and Neisseriaceae families. The genomic densities of these sequences are different both between species and between families. Interestingly, the family whose genomes harbor more of such sequences also shows more sequence types. A phylogenetic analysis allowed inferring the possible ancestral Neisseriacean sequence and a nucleotide-by-nucleotide analysis allowed inferring the potential ancestral Pasteurellacean sequence based on its genomic footprint. The method used for this work could be applied to other sequences, including transcription factor binding and repeated DNAs.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

Diversity patterns of bacteriophages infecting Aggregatibacter and Haemophilus species across clades and niches

Aggregatibacter and Haemophilus species are relevant human commensals and opportunistic pathogens. Consequently, their bacteriophages may have significant impact on human microbial ecology and pathologies. Our aim was to reveal the prevalence and diversity of bacteriophages infecting Aggregatibacter and Haemophilus species that colonize the human body. Genome mining with comparative genomics, screening of clinical isolates, and profiling of metagenomes allowed characterization of 346 phages grouped in 52 clusters and 18 superclusters. Less than 10% of the identified phage clusters were represented by previously characterized phages. Prophage diversity patterns varied significantly for different phage types, host clades, and environmental niches. A more diverse phage community lysogenizes Haemophilus influenzae and Haemophilus parainfluenzae strains than Aggregatibacter actinomycetemcomitans and “Haemophilus ducreyi”. Co-infections occurred more often in “H. ducreyi”. Phages from Aggregatibacter actinomycetemcomitans preferably lysogenized strains of specific serotype. Prophage patterns shared by subspecies clades of different bacterial species suggest similar ecoevolutionary drivers. Changes in frequencies of DNA uptake signal sequences and guanine–cytosine content reflect phage-host long-term coevolution. Aggregatibacter and Haemophilus phages were prevalent at multiple oral sites. Together, these findings should help exploring the ecoevolutionary forces shaping virus-host interactions in the human microbiome. Putative lytic phages, especially phiKZ-like, may provide new therapeutic options.

Efficient construction of Haemophilus parasuis mutants based on natural transformation

Studies on virulence factors and pathogenecity of Haemophilus parasuis have long been hindered by a lack of a consistent system for genetic manipulation. In this study, competence was induced by transferring H. parasuis from rich medium to starvation medium media-IV (M-IV) and iscR gene deficient mutants of H. parasuis were generated efficiently. Transformation frequency varied from 4.1 × 10^-5 to 1.1 × 10^-8 when using circular plasmid, and increased to about 2- to 31-fold when transformed using linearized plasmid. Allele replacement occurred efficiently in 6 strains, which are transformable using both circular and linearized pTRU, but not in another 2 strains which could only be transformed using linearized plasmid. The iscR mutants were stable for at least 20 passages in vitro. Haemophilus parasuis strains vary extensively in natural transformation efficiency and the method established here allows for transformation of a larger spectrum of strains with an easily accessed plasmid. This provides important tools for genetic manipulation of H. parasuis.

Natural competence and the evolution of DNA uptake specificity

Many bacteria are naturally competent, able to actively transport environmental DNA fragments across their cell envelope and into their cytoplasm. Because incoming DNA fragments can recombine with and replace homologous segments of the chromosome, competence provides cells with a potent mechanism of horizontal gene transfer as well as access to the nutrients in extracellular DNA. This review starts with an introductory overview of competence and continues with a detailed consideration of the DNA uptake specificity of competent proteobacteria in the Pasteurellaceae and Neisseriaceae. Species in these distantly related families exhibit strong preferences for genomic DNA from close relatives, a self-specificity arising from the combined effects of biases in the uptake machinery and genomic overrepresentation of the sequences this machinery prefers. Other competent species tested lack obvious uptake bias or uptake sequences, suggesting that strong convergent evolutionary forces have acted on these two families. Recent results show that uptake sequences have multiple "dialects," with clades within each family preferring distinct sequence variants and having corresponding variants enriched in their genomes. Although the genomic consensus uptake sequences are 12 and 29 to 34 bp, uptake assays have found that only central cores of 3 to 4 bp, conserved across dialects, are crucial for uptake. The other bases, which differ between dialects, make weaker individual contributions but have important cooperative interactions. Together, these results make predictions about the mechanism of DNA uptake across the outer membrane, supporting a model for the evolutionary accumulation and stability of uptake sequences and suggesting that uptake biases may be more widespread than currently thought.

Could DNA uptake be a side effect of bacterial adhesion and twitching motility?

DNA acquisition promotes the spread of resistance to antibiotics and virulence among bacteria. It is also linked to several natural phenomena including recombination, genome dynamics, adaptation and speciation. Horizontal DNA transfer between bacteria occurs via conjugation, transduction or competence for natural transformation by DNA uptake. Among these, competence is the only mechanism of transformation initiated and entirely controlled by the chromosome of the recipient bacteria. While the molecular mechanisms allowing the uptake of extracellular DNA are increasingly characterized, the function of competence for natural transformation by DNA uptake, the selective advantage maintaining it and the reasons why bacteria take up DNA in the first place are still debated. In this synthesis, I review some of the literature and discuss the four hypotheses on how and why do bacteria take up DNA. I argue that DNA uptake by bacteria is an accidental by-product of bacterial adhesion and twitching motility. Adhesion and motility are generally increased in stressful conditions, which may explain why bacteria increase DNA uptake in these conditions. In addition to its fundamental scientific relevance, the new hypothesis suggested here has significant clinical implications and finds further support from the fact that antibiotics sometimes fail to eliminate the targeted bacterium while inevitably causing stress to others. The widespread misuse of antibiotics may thus not only be selecting for resistant strains, but may also be causing bacteria to take up more DNA with the consequent increase in the chances of acquiring drug resistance and virulence-a scenario in full concordance with the previously reported induction of competence genes by antibiotics in Streptococcus pneumoniae and Legionella pneumophila.

Defining the DNA uptake specificity of naturally competent Haemophilus influenzae cells

Some naturally competent bacteria exhibit both a strong preference for DNA fragments containing specific ‘uptake sequences’ and dramatic overrepresentation of these sequences in their genomes. Uptake sequences are often assumed to directly reflect the specificity of the DNA uptake machinery, but the actual specificity has not been well characterized for any bacterium. We produced a detailed analysis of Haemophilus influenzae’s uptake specificity, using Illumina sequencing of degenerate uptake sequences in fragments recovered from competent cells. This identified an uptake motif with the same consensus as the motif overrepresented in the genome, with a 9 bp core (AAGTGCGGT) and two short flanking T-rich tracts. Only four core bases (GCGG) were critical for uptake, suggesting that these make strong specific contacts with the uptake machinery. Other core bases had weaker roles when considered individually, as did the T-tracts, but interaction effects between these were also determinants of uptake. The properties of genomic uptake sequences are also constrained by mutational biases and selective forces acting on USSs with coding and termination functions. Our findings define constraints on gene transfer by natural transformation and suggest how the DNA uptake machinery overcomes the physical constraints imposed by stiff highly charged DNA molecules.

Prosthetic valve endocarditis caused by Pasteurella in a penicillin allergic patient: challenges in diagnosis and treatment

Pasteurella multocida is a rare cause of infective endocarditis with only a few cases described. This report involves a 38-year-old penicillin-allergic patient in an immunocompromised state with several co-morbidities. Two molecular microbiological techniques, 16S rRNA sequencing and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry were used to confirm the species identification as P. multocida. Previous reports in the literature are also reviewed.

Bacterial DNA uptake sequences can accumulate by molecular drive alone

Uptake signal sequences are DNA motifs that promote DNA uptake by competent bacteria in the family Pasteurellaceae and the genus Neisseria. The genomes of these bacteria contain many copies of their canonical uptake sequence (often >100-fold overrepresentation), so the bias of the uptake machinery causes cells to prefer DNA derived from close relatives over DNA from other sources. However, the molecular and evolutionary forces responsible for the abundance of uptake sequences in these genomes are not well understood, and their presence is not easily explained by any of the current models of the evolution of competence. Here we describe use of a computer simulation model to thoroughly evaluate the simplest explanation for uptake sequences, that they accumulate in genomes by a form of molecular drive generated by biased DNA uptake and evolutionarily neutral (i.e., unselected) recombination. In parallel we used an unbiased search algorithm to characterize genomic uptake sequences and DNA uptake assays to refine the Haemophilus influenzae uptake specificity. These analyses showed that biased uptake and neutral recombination are sufficient to drive uptake sequences to high densities, with the spacings, stabilities, and strong consensus typical of uptake sequences in real genomes. This result greatly simplifies testing of hypotheses about the benefits of DNA uptake, because it explains how genomes could have passively accumulated sequences matching the bias of their uptake machineries.

Coevolution of DNA uptake sequences and bacterial proteomes

Dramatic examples of repeated sequences occur in the genomes of some naturally competent bacteria, which contain hundreds or thousands of copies of short motifs called DNA uptake signal sequences. Here, we analyze the evolutionary interactions between coding-region uptake sequences and the proteomes of Haemophilus influenzae, Actinobacillus pleuropneumoniae, and Neisseria meningitidis. In all three genomes, uptake sequence accumulation in coding sequences has approximately doubled the frequencies of those tripeptides specified by each species’ uptake sequence. The presence of uptake sequences in particular reading frames correlated most strongly with the use of preferred codons at degenerately coded positions, but the density of uptake sequences correlated only poorly with protein functional category. Genes lacking homologs in related genomes also lacked uptake sequences, strengthening the evidence that uptake sequences do not drive lateral gene transfer between distant relatives but instead accumulate after genes have been transferred. Comparison of the uptake sequence-encoded peptides of H. influenzae and N. meningitidis proteins with their homologs from related bacteria without uptake sequences indicated that uptake sequences were also preferentially located in poorly conserved genes and at poorly conserved amino acids. With few exceptions, amino acids at positions encoded by uptake sequences were as well conserved as other amino acids, suggesting that extant uptake sequences impose little or no constraint on coding for protein function. However, this state is likely to be achieved at a substantial cost because of the selective deaths required to eliminate maladaptive mutations that improve uptake sequences.

Lateral gene transfer of O1 serogroup encoding genes of Vibrio cholerae

In Gram-negative bacteria, the O-antigen-encoding genes may be transferred between lineages, although mechanisms are not fully understood. To assess possible lateral gene transfer (LGT), 21 Argentinean Vibrio cholerae O-group 1 (O1) isolates were examined using multilocus sequence typing (MLST) to determine the genetic relatedness of housekeeping genes and genes from the O1 gene cluster. MSLT analysis revealed that 4.4% of the nucleotides in the seven housekeeping loci were variable, with six distinct genetic lineages identified among O1 isolates. In contrast, MLST analysis of the eight loci from the O1 serogroup region revealed that 0.24% of the 4943 nucleotides were variable. A putative breakpoint was identified in the JUMPstart sequence. Nine conserved nucleotides differed by a single nucleotide from a DNA uptake signal sequence (USS) also found in Pastuerellaceae. Our data indicate that genes in the O1 biogenesis region are closely related even in distinct genetic lineages, indicative of LGT, with a putative DNA USS identified at the defined boundary for the DNA exchange.

An intragenic distribution bias of DNA uptake sequences in Pasteurellaceae and Neisseriae

Most sequenced strains from Pasteurellaceae and Neisseriae contain hundreds to thousands of uptake sequence (US) motifs in their genome, which are associated with natural competence for DNA uptake. The mechanism of their recognition is still unclear, and I searched for intragenic location patterns of these motifs for clues about their distribution. In all cases, one orientation of the US has a higher occurrence in the reading frame, and in all Pasteurellaceae, the US and the reverse complement motifs are biased towards the gene termini. These findings could help design experimental set-ups to study preferential DNA uptake, thereby further unravelling the phenomenon of natural competence.

The impact of the neisserial DNA uptake sequences on genome evolution and stability

A study of the origin and distribution of the abundant short DNA uptake sequence (DUS) in six genomes of Neisseria suggests that transformation and recombination are tightly linked in evolution and that recombination has a key role in the establishment of DUS.

Background

Efficient natural transformation in Neisseria requires the presence of short DNA uptake sequences (DUSs). Doubts remain whether DUSs propagate by pure selfish molecular drive or are selected for 'safe sex' among conspecifics.

Results

Six neisserial genomes were aligned to identify gene conversion fragments, DUS distribution, spacing, and conservation. We found a strong link between recombination and DUS: DUS spacing matches the size of conversion fragments; genomes with shorter conversion fragments have more DUSs and more conserved DUSs; and conversion fragments are enriched in DUSs. Many recent and singly occurring DUSs exhibit too high divergence with homologous sequences in other genomes to have arisen by point mutation, suggesting their appearance by recombination. DUSs are over-represented in the core genome, under-represented in regions under diversification, and absent in both recently acquired genes and recently lost core genes. This suggests that DUSs are implicated in genome stability rather than in generating adaptive variation. DUS elements are most frequent in the permissive locations of the core genome but are themselves highly conserved, undergoing mutation selection balance and/or molecular drive. Similar preliminary results were found for the functionally analogous uptake signal sequence in Pasteurellaceae.

Conclusion

As do many other pathogens, Neisseria and Pasteurellaceae have hyperdynamic genomes that generate deleterious mutations by intrachromosomal recombination and by transient hypermutation. The results presented here suggest that transformation in Neisseria and Pasteurellaceae allows them to counteract the deleterious effects of genome instability in the core genome. Thus, rather than promoting hypervariation, bacterial sex could be regenerative.

Local Renyi entropic profiles of DNA sequences

Background

In a recent report the authors presented a new measure of continuous entropy for DNA sequences, which allows the estimation of their randomness level. The definition therein explored was based on the Rényi entropy of probability density estimation (pdf) using the Parzen's window method and applied to Chaos Game Representation/Universal Sequence Maps (CGR/USM). Subsequent work proposed a fractal pdf kernel as a more exact solution for the iterated map representation. This report extends the concepts of continuous entropy by defining DNA sequence entropic profiles using the new pdf estimations to refine the density estimation of motifs.

Results

The new methodology enables two results. On the one hand it shows that the entropic profiles are directly related with the statistical significance of motifs, allowing the study of under and over-representation of segments. On the other hand, by spanning the parameters of the kernel function it is possible to extract important information about the scale of each conserved DNA region. The computational applications, developed in Matlab m-code, the corresponding binary executables and additional material and examples are made publicly available at .

Conclusion

The ability to detect local conservation from a scale-independent representation of symbolic sequences is particularly relevant for biological applications where conserved motifs occur in multiple, overlapping scales, with significant future applications in the recognition of foreign genomic material and inference of motif structures.

Genome dynamics of short oligonucleotides: the example of bacterial DNA uptake enhancing sequences

Among the many bacteria naturally competent for transformation by DNA uptake-a phenomenon with significant clinical and financial implications- Pasteurellaceae and Neisseriaceae species preferentially take up DNA containing specific short sequences. The genomic overrepresentation of these DNA uptake enhancing sequences (DUES) causes preferential uptake of conspecific DNA, but the function(s) behind this overrepresentation and its evolution are still a matter for discovery. Here I analyze DUES genome dynamics and evolution and test the validity of the results to other selectively constrained oligonucleotides. I use statistical methods and computer simulations to examine DUESs accumulation in Haemophilus influenzae and Neisseria gonorrhoeae genomes. I analyze DUESs sequence and nucleotide frequencies, as well as those of all their mismatched forms, and prove the dependence of DUESs genomic overrepresentation on their preferential uptake by quantifying and correlating both characteristics. I then argue that mutation, uptake bias, and weak selection against DUESs in less constrained parts of the genome combined are sufficient enough to cause DUESs accumulation in susceptible parts of the genome with no need for other DUES function. The distribution of overrepresentation values across sequences with different mismatch loads compared to the DUES suggests a gradual yet not linear molecular drive of DNA sequences depending on their similarity to the DUES. Other genomically overrepresented sequences, both pro- and eukaryotic, show similar distribution of frequencies suggesting that the molecular drive reported above applies to other frequent oligonucleotides. Rare oligonucleotides, however, seem to be gradually drawn to genomic underrepresentation, thus, suggesting a molecular drag. To my knowledge this work provides the first clear evidence of the gradual evolution of selectively constrained oligonucleotides, including repeated, palindromic and protein/transcription factor-binding DNAs.

Rapid, accurate, computational discovery of Rho-independent transcription terminators illuminates their relationship to DNA uptake

Using a novel computational method, an extensive collection of predicted Rho-independent transcription terminators is derived from 343 prokaryotes, offering insight into their relationship to DNA uptake

Background

In many prokaryotes, transcription of DNA to RNA is terminated by a thymine-rich stretch of DNA following a hairpin loop. Detecting such Rho-independent transcription terminators can shed light on the organization of bacterial genomes and can improve genome annotation. Previous computational methods to predict Rho-independent terminators have been slow or limited in the organisms they consider.

Results

We describe TransTermHP, a new computational method to rapidly and accurately detect Rho-independent transcription terminators. We predict the locations of terminators in 343 prokaryotic genomes, representing the largest collection of predictions available. In Bacillus subtilis, we can detect 93% of known terminators with a false positive rate of just 6%, comparable to the best-known methods. Outside the Firmicutes division, we find that Rho-independent termination plays a large role in the Neisseria and Vibrio genera, the Pasteurellaceae (including the Haemophilus genus) and several other species. In Neisseria and Pasteurellaceae, terminator hairpins are frequently formed by closely spaced, complementary instances of exogenous DNA uptake signal sequences. We quantify the propensity for terminators to include these sequences. In the process, we provide the first discussion of potential uptake signals in Haemophilus ducreyi and Mannheimia succiniciproducens, and we discuss the preference for a particular configuration of uptake signal sequences within terminators.

Conclusion

Our new fast and accurate method for detecting transcription terminators has allowed us to identify and analyze terminators in many new genomes and to identify DNA uptake signal sequences in several species where they have not been previously reported. Our software and predictions are freely available.

tfoX (sxy)-dependent transformation of Aggregatibacter (Actinobacillus) actinomycetemcomitans

tfoX (sxy) is a regulatory gene needed to turn on competence genes. Aggregatibacter (Actinobacillus) actinomycetemcomitans has a tfoX gene that is important for transformation. We cloned this gene on an IncQ plasmid downstream of the inducible tac promoter. When this plasmid was resident in cells of A. actinomycetemcomitans and tfoX was induced, the cells became competent for transformation. Several strains of A. actinomycetemcomitans, including different serotypes, as well as rough (adherent) and isogenic smooth (nonadherent) forms were tested. Only our two serotype f strains failed to be transformed. With the other strains, we could easily get transformants with extrachromosomal plasmid DNA when closed circular, replicative plasmid carrying an uptake signal sequence (USS) was used. When a replicative plasmid carrying a USS and cloned DNA from the chromosome of A. actinomycetemcomitans was linearized by digestion with a restriction endonuclease or when genomic DNA was used directly, the outcome was allelic exchange. To facilitate allelic exchange, we constructed a suicide plasmid (pMB78) that does not replicate in A. actinomycetemcomitans and carries a region with two inverted copies of a USS. This vector gave allelic exchange in the presence of cloned and induced tfoX easily and without digestion. Using transposon insertions in cloned katA DNA, we found that as little as 78 bp of homology at one of the ends was sufficient for that end to participate in allelic exchange. The cloning and induction of tfoX makes it possible to transform nearly any strain of A. actinomycetemcomitans, and allelic exchange has proven to be important for site-directed mutagenesis.

Horizontal gene transfer of ftsI, encoding penicillin-binding protein 3, in Haemophilus influenzae

Horizontal gene transfer has been identified in only a small number of genes in Haemophilus influenzae, an organism which is naturally competent for transformation. This report provides evidence for the genetic transfer of the ftsI gene, which encodes penicillin-binding protein 3, in H. influenzae. Mosaic structures of the ftsI gene were found in several clinical isolates of H. influenzae. To identify the origin of the mosaic sequence, complete sequences of the corresponding gene from seven type strains of Haemophilus species were determined. Comparison of these sequences with mosaic regions identified a homologous recombination of the ftsI gene between H. influenzae and Haemophilus haemolyticus. Subsequently, ampicillin-resistant H. influenzae strains harboring identical ftsI sequences were genotyped by pulsed-field gel electrophoresis (PFGE). Divergent PFGE patterns among beta-lactamase-nonproducing ampicillin-resistant (BLNAR) strains from different hospitals indicated the potential for the genetic transfer of the mutated ftsI gene between these isolates. Moreover, transfer of the ftsI gene from BLNAR strains to beta-lactamase-nonproducing ampicillin-susceptible (BLNAS) H. influenzae strains was evaluated in vitro. Coincubation of a BLNAS strain (a rifampin-resistant mutant of strain Rd) and BLNAR strains resulted in the emergence of rifampin- and cefdinir-resistant clones at frequencies of 5.1 x 10(-7) to 1.5 x 10(-6). Characterization of these doubly resistant mutants by DNA sequencing of the ftsI gene, susceptibility testing, and genotyping by PFGE revealed that the ftsI genes of BLNAR strains had transferred to BLNAS strains during coincubation. In conclusion, horizontal transfer of the ftsI gene in H. influenzae can occur in an intraspecies and an interspecies manner.

Genomic distribution and functions of uptake signal sequences in Actinobacillus actinomycetemcomitans

Actinobacillus actinomycetemcomitans is naturally competent for transformation, with a transformation system similar to that of Haemophilus influenzae that preferentially takes up DNA bearing uptake signal sequences (USS) with the same 9-base USS core. This study examined the function of the extended 29-base USS, which comprises a highly conserved 1st region (containing the 9-base core) and 2nd and 3rd semi-conserved AT-rich regions, in transformation of A. actinomycetemcomitans. Transformation frequency was not affected by either location (in middle or at 5' end) or quantity (one or two) of USS in donor DNA. Relative transformation efficiencies (in comparison to the positive control) were 28-67 % for linear DNA with single-base mutations in the USS 1st region, and 47 % and 73 %, respectively, for linear DNA with USS that contained either a non-consensus 2nd or a non-consensus 3rd region. Plasmids with a stand-alone 1st or a stand-alone 2nd-3rd region exhibited 21 % and 6 % relative transformation efficiencies, respectively. It was also noted that A. actinomycetemcomitans and H. influenzae were similar in the frequencies and distribution patterns of USS in their genomes. In conclusion, all three regions of the extended 29-base USS are required for optimum transformation in A. actinomycetemcomitans.

The genome sequence of Mannheimia haemolytica A1: insights into virulence, natural competence, and Pasteurellaceae phylogeny

The draft genome sequence of Mannheimia haemolytica A1, the causative agent of bovine respiratory disease complex (BRDC), is presented. Strain ATCC BAA-410, isolated from the lung of a calf with BRDC, was the DNA source. The annotated genome includes 2,839 coding sequences, 1,966 of which were assigned a function and 436 of which are unique to M. haemolytica. Through genome annotation many features of interest were identified, including bacteriophages and genes related to virulence, natural competence, and transcriptional regulation. In addition to previously described virulence factors, M. haemolytica encodes adhesins, including the filamentous hemagglutinin FhaB and two trimeric autotransporter adhesins. Two dual-function immunoglobulin-protease/adhesins are also present, as is a third immunoglobulin protease. Genes related to iron acquisition and drug resistance were identified and are likely important for survival in the host and virulence. Analysis of the genome indicates that M. haemolytica is naturally competent, as genes for natural competence and DNA uptake signal sequences (USS) are present. Comparison of competence loci and USS in other species in the family Pasteurellaceae indicates that M. haemolytica, Actinobacillus pleuropneumoniae, and Haemophilus ducreyi form a lineage distinct from other Pasteurellaceae. This observation was supported by a phylogenetic analysis using sequences of predicted housekeeping genes.

Evolution of competence and DNA uptake specificity in the Pasteurellaceae

BACKGROUND

Many bacteria can take up DNA, but the evolutionary history and function of natural competence and transformation remain obscure. The sporadic distribution of competence suggests it is frequently lost and/or gained, but this has not been examined in an explicitly phylogenetic context. Additional insight may come from the sequence specificity of uptake by species such as Haemophilus influenzae, where a 9 bp uptake signal sequence (USS) repeat is both highly overrepresented in the genome and needed for efficient DNA uptake. We used the distribution of competence genes and DNA uptake specificity in H. influenzae's family, the Pasteurellaceae, to examine the ancestry of competence.

RESULTS

A phylogeny of the Pasteurellaceae based on 12 protein coding genes from species with sequenced genomes shows two strongly supported subclades: the Hin subclade (H. influenzae, Actinobacillus actinomycetemcomitans, Pasteurella multocida, Mannheimia succiniciproducens, and H. somnus), and the Apl subclade (A. pleuropneumoniae, M. haemolytica, and H. ducreyi). All species contained homologues of all known H. influenzae competence genes, consistent with an ancestral origin of competence. Competence gene defects were identified in three species (H. somnus, H. ducreyi and M. haemolytica); each appeared to be of recent origin. The assumption that USS arise by mutation rather than copying was first confirmed using alignments of H. influenzae proteins with distant homologues. Abundant USS-like repeats were found in all eight Pasteurellacean genomes; the repeat consensuses of species in the Hin subclade were identical to that of H. influenzae (AAGTGCGGT), whereas members of the Apl subclade shared the consensus ACAAGCGGT. All species' USSs had the strong consensus and flanking AT-rich repeats of H. influenzae USSs. DNA uptake and competition experiments demonstrated that the Apl-type repeat is a true USS distinct from the Hin-type USS: A. pleuropneumoniae preferentially takes up DNA fragments containing the Apl-type USS over both H. influenzae and unrelated DNAs, and H. influenzae prefers its own USS over the Apl type.

CONCLUSION

Competence and DNA uptake specificity are ancestral properties of the Pasteurellaceae, with divergent USSs and uptake specificity distinguishing only the two major subclades. The conservation of most competence genes over the approximately 350 million year history of the family suggests that lineages that lose competence may be evolutionary dead ends.

Genotypic diversity of Haemophilus parasuis field strains

Haemophilus parasuis is the cause of Glässer's disease and other clinical disorders in pigs. It can also be isolated from the upper respiratory tracts of healthy pigs, and isolates can have significant differences in virulence. In this work, a partial sequence from the 60-kDa heat shock protein (Hsp60) gene was assessed as an epidemiological marker. We analyzed partial sequences of hsp60 and 16S rRNA genes from 103 strains of H. parasuis and other related species to obtain a better classification of the strains and examine the correlation with virulence. The results were compared with those obtained by enterobacterial repetitive intergenic consensus PCR. Our results showed that hsp60 is a reliable marker for epidemiological studies of H. parasuis and that the analysis of its sequence is a better approach than fingerprinting methods. Furthermore, the analysis of the hsp60 and 16S rRNA gene sequences revealed the presence of a separate lineage of virulent strains and indicated the occurrence of lateral gene transfer among H. parasuis and Actinobacillus strains.

Variations in gene organization and DNA uptake signal sequence in the folP region between commensal and pathogenic Neisseria species

Background

Horizontal gene transfer is an important source of genetic variation among Neisseria species and has contributed to the spread of resistance to penicillin and sulfonamide drugs in the pathogen Neisseria meningitidis. Sulfonamide resistance in Neisseria meningitidis is mediated by altered chromosomal folP genes. At least some folP alleles conferring resistance have been horizontally acquired from other species, presumably from commensal Neisseriae. In this work, the DNA sequence surrounding folP in commensal Neisseria species was determined and compared to corresponding regions in pathogenic Neisseriae, in order to elucidate the potential for inter-species DNA transfer within this region.

Results

The upstream region of folP displayed differences in gene order between species, including an insertion of a complete Correia element in Neisseria lactamica and an inversion of a larger genomic segment in Neisseria sicca, Neisseria subflava and Neisseria mucosa. The latter species also had DNA uptake signal sequences (DUS) in this region that were one base different from the DUS in pathogenic Neisseriae. Another interesting finding was evidence of a horizontal transfer event from Neisseria lactamica or Neisseria cinerea that introduced a novel folP allele to the meningococcal population.

Conclusion

Genetic recombination events immediately upstream of folP and horizontal transfer have resulted in sequence differences in the folP region between the Neisseria species. This variability could be a consequence of the selective pressure on this region exerted by the use of sulfonamide drugs.

High level of sequence diversity in the 16S rRNA genes of Haemophilus influenzae isolates is useful for molecular subtyping

A molecular typing method based on the 16S rRNA sequence diversity was developed for Haemophilus influenzae isolates. A total of 330 H. influenzae isolates were analyzed, representing a diverse collection of U.S. isolates. We found a high level of 16S rRNA sequence heterogeneity (up to 2.73%) and observed an exclusive correlation between 16S types and serotypes (a to f); no 16S type was found in more than one serotype. Similarly, no multilocus sequence typing (MLST) sequence type (ST) was found in more than one serotype. Our 16S typing and MLST results are in agreement with those of previous studies showing that serotypable H. influenzae isolates behave as highly clonal populations and emphasize the lack of clonality of nontypable (NT) H. influenzae isolates. There was not a 1:1 correlation between 16S types and STs, but all H. influenzae serotypable isolates clustered similarly. This correlation was not observed for NT H. influenzae; the two methods clustered NT H. influenzae isolates differently. 16S rRNA gene sequencing alone provides a level of discrimination similar to that obtained with the analysis of seven genes for MLST. We demonstrated that 16S typing is an additional and complementary approach to MLST, particularly for NT H. influenzae isolates, and is potentially useful for outbreak investigation.

Evaluation of the current models for the evolution of bacterial DNA uptake signal sequences

Current opinion considers two main hypotheses for the evolutionary origin of uptake signal sequences in bacteria: one model regards the uptake signal sequence (USS) as the result of biased gene conversion, whereas the second model views the USS as a molecular tag that evolved as an adaptation. In this article, we present various computational models that implement specific versions of those hypotheses. Those models show that the two hypothesis are not necessarily as opposed to each other as may appear at first glance.

Identification, distribution, and expression of novel genes in 10 clinical isolates of nontypeable Haemophilus influenzae

We hypothesize that Haemophilus influenzae, as a species, possesses a much greater number of genes than that found in any single H. influenzae genome. This supragenome is distributed throughout naturally occurring infectious populations, and new strains arise through autocompetence and autotransformation systems. The effect is that H. influenzae populations can readily adapt to environmental stressors. The supragenome hypothesis predicts that significant differences exist between and among the genomes of individual infectious strains of nontypeable H. influenzae (NTHi). To test this prediction, we obtained 10 low-passage NTHi clinical isolates from the middle ear effusions of patients with chronic otitis media. DNA sequencing was performed with 771 clones chosen at random from a pooled genomic library. Homology searching demonstrated that approximately 10% of these clones were novel compared to the H. influenzae Rd KW20 genome, and most of them did not match any DNA sequence in GenBank. Amino acid homology searches using hypothetical translations of the open reading frames revealed homologies to a variety of proteins, including bacterial virulence factors not previously identified in the NTHi isolates. The distribution and expression of 53 of these genes among the 10 strains were determined by PCR- and reverse transcription PCR-based analyses. These unique genes were nonuniformly distributed among the 10 isolates, and transcription of these genes in planktonic cultures was detected in 50% (177 of 352) of the occurrences. All of the novel sequences were transcribed in one or more of the NTHi isolates. Seventeen percent (9 of 53) of the novel genes were identified in all 10 NTHi strains, with each of the remaining 44 being present in only a subset of the strains. These genic distribution analyses were more effective as a strain discrimination tool than either multilocus sequence typing or 23S ribosomal gene typing methods.

Genome rearrangement distances and gene order phylogeny in gamma-Proteobacteria

Genome rearrangements have been studied in 30 gamma-proteobacterial complete genomes by comparing the order of a reduced set of genes on the chromosome. This set included those genes fulfilling several characteristics, the main ones being that an ortholog was present in every genome and that none of them had been acquired by horizontal gene transfer. Genome rearrangement distances were estimated based on either the number of breakpoints or the minimal number of inversions separating two genomes. Breakpoint and inversion distances were highly correlated, indicating that inversions were the main type of rearrangement event in gamma-Proteobacteria. In general, the progressive increase in sequence-based distances between genome pairs was associated with the increase in their rearrangement-based distances but with several groups of distances not following this pattern. Compared with free-living enteric bacteria, the lineages of Pasteurellaceae were evolving, on average, to relatively higher rates of between 2.02 and 1.64, while the endosymbiotic bacterial lineages of Buchnera aphidicola and Wigglesworthia glossinidia were evolving at moderately higher rates of 1.38 and 1.35, respectively. Because we know that the rearrangement rate in the Bu. aphidicola lineage was close to zero during the last 100-150 Myr of evolution, we deduced that a much higher rate took place in the first period of lineage evolution after the divergence of the Escherichia coli lineage. On the other hand, the lineage of the endosymbiont Blochmannia floridanus did present an almost identical rate to free-living enteric bacteria, indicating that the increase in the genome rearrangement rate is not a general change associated with bacterial endosymbiosis. Phylogenetic reconstruction based on rearrangement distances showed a different topology from the one inferred by sequence information. This topology broke the proposed monophyly of the three endosymbiotic lineages and placed Bl. floridanus as a closer relative to E. coli than Yersinia pestis. These results indicate that the phylogeny of these insect endosymbionts is still an open question that will require the development of specific phylogenetic methods to confirm whether the sisterhood of the three endosymbiotic lineages is real or a consequence of a long-branch attraction phenomenon.

Evolution of DNA uptake signal sequences

The DNA of some naturally competent species of bacteria contains a large number of evenly distributed copies of a short sequence. This highly overrepresented sequence is believed to be an uptake signal sequence (USS) that helps bacteria to take up DNA selectively from (dead) members of their own species. For some time it has been assumed that the USS evolved in order to enable bacteria to distinguish between conspecific and nonconspecific DNA fragments (the preference-first hypothesis). Recently, Redfield suggested that this hypothesis is not in fact realistic, as it would require biologically implausible group selection. In this article we present a model designed to demonstrate the emergence of similar USSs in a population of simulated evolving agents. We use this model to examine the conditions under which a USS will emerge in a preference-first scenario.