Microarray-based comparative genomic profiling of reference strains and selected Canadian field isolates of Actinobacillus pleuropneumoniae.

BACKGROUND: Actinobacillus pleuropneumoniae, the causative agent of porcine pleuropneumonia, is a highly contagious respiratory pathogen that causes severe losses to the swine industry worldwide. Current commercially-available vaccines are of limited value because they do not induce cross-serovar immunity and do not prevent development of the carrier state. Microarray-based comparative genomic hybridizations (M-CGH) were used to estimate whole genomic diversity of representative Actinobacillus pleuropneumoniae strains. Our goal was to identify conserved genes, especially those predicted to encode outer membrane proteins and lipoproteins because of their potential for the development of more effective vaccines. RESULTS: Using hierarchical clustering, our M-CGH results showed that the majority of the genes in the genome of the serovar 5 A. pleuropneumoniae L20 strain were conserved in the reference strains of all 15 serovars and in representative field isolates. Fifty-eight conserved genes predicted to encode for outer membrane proteins or lipoproteins were identified. As well, there were several clusters of diverged or absent genes including those associated with capsule biosynthesis, toxin production as well as genes typically associated with mobile elements. CONCLUSION: Although A. pleuropneumoniae strains are essentially clonal, M-CGH analysis of the reference strains of the fifteen serovars and representative field isolates revealed several classes of genes that were divergent or absent. Not surprisingly, these included genes associated with capsule biosynthesis as the capsule is associated with sero-specificity. Several of the conserved genes were identified as candidates for vaccine development, and we conclude that M-CGH is a valuable tool for reverse vaccinology.

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.

Comparative genomics profiling of clinical isolates of Aeromonas salmonicida using DNA microarrays.

BACKGROUND: Aeromonas salmonicida has been isolated from numerous fish species and shows wide variation in virulence and pathogenicity. As part of a larger research program to identify virulence genes and candidates for vaccine development, a DNA microarray was constructed using a subset of 2024 genes from the draft genome sequence of A. salmonicida subsp. salmonicida strain A449. The microarray included genes encoding known virulence-associated factors in A. salmonicida and homologs of virulence genes of other pathogens. We used microarray-based comparative genomic hybridizations (M-CGH) to compare selected A. salmonicida sub-species and other Aeromonas species from different hosts and geographic locations. RESULTS: Results showed variable carriage of virulence-associated genes and generally increased variation in gene content across sub-species and species boundaries. The greatest variation was observed among genes associated with plasmids and transposons. There was little correlation between geographic region and degree of variation for all isolates tested. CONCLUSION: We have used the M-CGH technique to identify subsets of conserved genes from amongst this set of A. salmonicida virulence genes for further investigation as potential vaccine candidates. Unlike other bacterial characterization methods that use a small number of gene or DNA-based functions, M-CGH examines thousands of genes and/or whole genomes and thus is a more comprehensive analytical tool for veterinary or even human health research.

A new approach for the analysis of bacterial microarray-based Comparative Genomic Hybridization: insights from an empirical study.

BACKGROUND: Microarray-based Comparative Genomic Hybridization (M-CGH) has been used to characterize the extensive intraspecies genetic diversity found in bacteria at the whole-genome level. Although conventional microarray analytical procedures have proved adequate in handling M-CGH data, data interpretation using these methods is based on a continuous character model in which gene divergence and gene absence form a spectrum of decreasing gene conservation levels. However, whereas gene divergence may yet be accompanied by retention in gene function, gene absence invariably leads to loss of function. This distinction, if ignored, leads to a loss in the information to be gained from M-CGH data. We present here results from experiments in which two genome-sequenced strains of C. jejuni were compared against each other using M-CGH. Because the gene content of both strains was known a priori, we were able to closely examine the effects of sequence divergence and gene absence on M-CGH data in order to define analytical parameters for M-CGH data interpretation. This would facilitate the examination of the relative effects of sequence divergence or gene absence in comparative genomics analyses of multiple strains of any species for which genome sequence data and a DNA microarray are available. RESULTS: As a first step towards improving the analysis of M-CGH data, we estimated the degree of experimental error in a series of experiments in which identical samples were compared against each other by M-CGH. This variance estimate was used to validate a Log Ratio-based methodology for identification of outliers in M-CGH data. We compared two genome strains by M-CGH to examine the effect of probe/target identity on the Log Ratios of signal intensities using prior knowledge of gene divergence and gene absence to establish Log Ratio thresholds for the identification of absent and conserved genes. CONCLUSION: The results from this empirical study validate the Log Ratio thresholds that have been used in other studies to establish gene divergence/absence. Moreover, the analytical framework presented here enhances the information content derived from M-CGH data by shifting the focus from divergent/absent gene detection to accurate detection of conserved and absent genes. This approach closely aligns the technical limitations of M-CGH analysis with practical limitations on the biological interpretation of comparative genomics data.

Large-scale comparative genomics meta-analysis of Campylobacter jejuni isolates reveals low level of genome plasticity.

We have used comparative genomic hybridization (CGH) on a full-genome Campylobacter jejuni microarray to examine genome-wide gene conservation patterns among 51 strains isolated from food and clinical sources. These data have been integrated with data from three previous C. jejuni CGH studies to perform a meta-analysis that included 97 strains from the four separate data sets. Although many genes were found to be divergent across multiple strains (n = 350), many genes (n = 249) were uniquely variable in single strains. Thus, the strains in each data set comprise strains with a unique genetic diversity not found in the strains in the other data sets. Despite the large increase in the collective number of variable C. jejuni genes (n = 599) found in the meta-analysis data set, nearly half of these (n = 276) mapped to previously defined variable loci, and it therefore appears that large regions of the C. jejuni genome are genetically stable. A detailed analysis of the microarray data revealed that divergent genes could be differentiated on the basis of the amplitudes of their differential microarray signals. Of 599 variable genes, 122 could be classified as highly divergent on the basis of CGH data. Nearly all highly divergent genes (117 of 122) had divergent neighbors and showed high levels of intraspecies variability. The approach outlined here has enabled us to distinguish global trends of gene conservation in C. jejuni and has enabled us to define this group of genes as a robust set of variable markers that can become the cornerstone of a new generation of genotyping methods that use genome-wide C. jejuni gene variability data.

Genome-wide expression analyses of Campylobacter jejuni NCTC11168 reveals coordinate regulation of motility and virulence by flhA.

We examined two variants of the genome-sequenced strain, Campylobacter jejuni NCTC11168, which show marked differences in their virulence properties including colonization of poultry, invasion of Caco-2 cells, and motility. Transcript profiles obtained from whole genome DNA microarrays and proteome analyses demonstrated that these differences are reflected in late flagellar structural components and in virulence factors including those involved in flagellar glycosylation and cytolethal distending toxin production. We identified putative sigma(28) and sigma(54) promoters for many of the affected genes and found that greater differences in expression were observed for sigma(28)-controlled genes. Inactivation of the gene encoding sigma(28), fliA, resulted in an unexpected increase in transcripts with sigma(54) promoters, as well as decreased transcription of sigma(28)-regulated genes. This was unlike the transcription profile observed for the attenuated C. jejuni variant, suggesting that the reduced virulence of this organism was not entirely due to impaired function of sigma(28). However, inactivation of flhA, an important component of the flagellar export apparatus, resulted in expression patterns similar to that of the attenuated variant. These findings indicate that the flagellar regulatory system plays an important role in campylobacter pathogenesis and that flhA is a key element involved in the coordinate regulation of late flagellar genes and of virulence factors in C. jejuni.