Bacterial Genomic Islands: Organization,Function, and Evolutionary Role

Data on the structural organization and evolutionary role of specific bacterial DNA regions known as genomic islands are reviewed. Emphasis is placed on the most extensively studied genomic islands, pathogenicity islands (PAIs), which are present in the chromosome of Gram-negative and Gram-positive pathogenic bacteria and absent from related nonpathogenic strains. PAIs are long DNA regions that harbor virulence genes and often differ in GC content from the remainder of the bacterial genome. Many PAI occur in the tRNA gene loci, which provide a convenient target for foreign gene insertion. Some PAI are highly homologous to each other and contain sequences similar to ISs, phage att sites, and plasmid ori sites, along with functional or defective integrase and transposase genes, suggesting horizontal transfer of PAI among bacteria.     

An evolutionary analysis of genome expansion and pathogenicity in Escherichia coli

Background

There are several studies describing loss of genes through reductive evolution in microbes, but how selective forces are associated with genome expansion due to horizontal gene transfer (HGT) has not received similar attention. The aim of this study was therefore to examine how selective pressures influence genome expansion in 53 fully sequenced and assembled Escherichia coli strains. We also explored potential connections between genome expansion and the attainment of virulence factors. This was performed using estimations of several genomic parameters such as AT content, genomic drift (measured using relative entropy), genome size and estimated HGT size, which were subsequently compared to analogous parameters computed from the core genome consisting of 1729 genes common to the 53 E. coli strains. Moreover, we analyzed how selective pressures (quantified using relative entropy and dN/dS), acting on the E. coli core genome, influenced lineage and phylogroup formation.

Results

Hierarchical clustering of dS and dN estimations from the E. coli core genome resulted in phylogenetic trees with topologies in agreement with known E. coli taxonomy and phylogroups. High values of dS, compared to dN, indicate that the E. coli core genome has been subjected to substantial purifying selection over time; significantly more than the non-core part of the genome (p<0.001). This is further supported by a linear association between strain-wise dS and dN values (β = 26.94 ± 0.44, R²~0.98, p<0.001). The non-core part of the genome was also significantly more AT-rich (p<0.001) than the core genome and E. coli genome size correlated with estimated HGT size (p<0.001). In addition, genome size (p<0.001), AT content (p<0.001) as well as estimated HGT size (p<0.005) were all associated with the presence of virulence factors, suggesting that pathogenicity traits in E. coli are largely attained through HGT. No associations were found between selective pressures operating on the E. coli core genome, as estimated using relative entropy, and genome size (p~0.98).

Conclusions

On a larger time frame, genome expansion in E. coli, which is significantly associated with the acquisition of virulence factors, appears to be independent of selective forces operating on the core genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-882) contains supplementary material, which is available to authorized users.     

Life on Arginine for Mycoplasma hominis: Clues from Its Minimal Genome and Comparison with Other Human Urogenital Mycoplasmas

Mycoplasma hominis is an opportunistic human mycoplasma. Two other pathogenic human species, M. genitalium and Ureaplasma parvum, reside within the same natural niche as M. hominis: the urogenital tract. These three species have overlapping, but distinct, pathogenic roles. They have minimal genomes and, thus, reduced metabolic capabilities characterized by distinct energy-generating pathways. Analysis of the M. hominis PG21 genome sequence revealed that it is the second smallest genome among self-replicating free living organisms (665,445 bp, 537 coding sequences (CDSs)). Five clusters of genes were predicted to have undergone horizontal gene transfer (HGT) between M. hominis and the phylogenetically distant U. parvum species. We reconstructed M. hominis metabolic pathways from the predicted genes, with particular emphasis on energy-generating pathways. The Embden–Meyerhoff–Parnas pathway was incomplete, with a single enzyme absent. We identified the three proteins constituting the arginine dihydrolase pathway. This pathway was found essential to promote growth in vivo. The predicted presence of dimethylarginine dimethylaminohydrolase suggested that arginine catabolism is more complex than initially described. This enzyme may have been acquired by HGT from non-mollicute bacteria. Comparison of the three minimal mollicute genomes showed that 247 CDSs were common to all three genomes, whereas 220 CDSs were specific to M. hominis, 172 CDSs were specific to M. genitalium, and 280 CDSs were specific to U. parvum. Within these species-specific genes, two major sets of genes could be identified: one including genes involved in various energy-generating pathways, depending on the energy source used (glucose, urea, or arginine) and another involved in cytadherence and virulence. Therefore, a minimal mycoplasma cell, not including cytadherence and virulence-related genes, could be envisaged containing a core genome (247 genes), plus a set of genes required for providing energy. For M. hominis, this set would include 247+9 genes, resulting in a theoretical minimal genome of 256 genes.     

Bacterial genomic islands: organization, function, and role in evolution

Data on the structural organization and evolutionary role of specific bacterial DNA regions known as genomic islands are reviewed. Emphasis is placed on the most extensively studied genomic islands, pathogenicity islands (PAIs), which are present in the chromosome of Gram-negative and Gram-positive pathogenic bacteria and absent from related nonpathogenic strains. PAIs are extended DNA regions that harbor virulence genes and often differ in GC content from the remainder of the bacterial genome. Many PAI occur in the tRNA genes, which provide a convenient target for foreign gene insertion. Some PAI are highly homologous to each other and contain sequences similar to ISs, phage att sites, and plasmid ori sites, along with functional or defective integrase and transposase genes, suggesting horizontal transfer of PAI among bacteria.     

A novel method to generate unmarked gene deletions in the intracellular pathogen Rhodococcus equi using 5-fluorocytosine conditional lethality

A novel method to efficiently generate unmarked in-frame gene deletions in Rhodococcus equi was developed, exploiting the cytotoxic effect of 5-fluorocytosine (5-FC) by the action of cytosine deaminase (CD) and uracil phosphoribosyltransferase (UPRT) enzymes. The opportunistic, intracellular pathogen R. equi is resistant to high concentrations of 5-FC. Introduction of Escherichia coli genes encoding CD and UPRT conferred conditional lethality to R. equi cells incubated with 5-FC. To exemplify the use of the codA::upp cassette as counter-selectable marker, an unmarked in-frame gene deletion mutant of R. equi was constructed. The supA and supB genes, part of a putative cholesterol catabolic gene cluster, were efficiently deleted from the R. equi wild-type genome. Phenotypic analysis of the generated ΔsupAB mutant confirmed that supAB are essential for growth of R. equi on cholesterol. Macrophage survival assays revealed that the ΔsupAB mutant is able to survive and proliferate in macrophages comparable to wild type. Thus, cholesterol metabolism does not appear to be essential for macrophage survival of R. equi. The CD-UPRT based 5-FC counter-selection may become a useful asset in the generation of unmarked in-frame gene deletions in other actinobacteria as well, as actinobacteria generally appear to be 5-FC resistant and 5-FU sensitive.     

Conjugal Transfer of a Virulence Plasmid in the Opportunistic Intracellular Actinomycete Rhodococcus equi

Rhodococcus equi is a facultative intracellular, Gram-positive, soilborne actinomycete which can cause severe pyogranulomatous pneumonia with abscessation in young horses (foals) and in immunocompromised people, such as persons with AIDS. All strains of R. equi isolated from foals and approximately a third isolated from humans contain a large, ∼81-kb plasmid which is essential for the intramacrophage growth of the organism and for virulence in foals and murine in vivo model systems. We found that the entire virulence plasmid could be transferred from plasmid-containing strains of R. equi (donor) to plasmid-free R. equi strains (recipient) at a high frequency and that plasmid transmission reestablished the capacity for intracellular growth in macrophages. Plasmid transfer required living cells and cell-to-cell contact and was unaffected by the presence of DNase, factors pointing to conjugation as the major means of genetic transfer. Deletion of a putative relaxase-encoding gene, traA, located in the proposed conjugative region of the plasmid, abolished plasmid transfer. Reversion of the traA mutation restored plasmid transmissibility. Finally, plasmid transmission to other Rhodococcus species and some additional related organisms was demonstrated. This is the first study showing a virulence plasmid transfer in R. equi, and it establishes a mechanism by which the virulence plasmid can move among bacteria in the soil.     

Molecular Characterization of Rhodococcus equi from Horse-Breeding Farms by Means of Multiplex PCR for the vap Gene Family

This study evaluated the molecular characteristics of Rhodococcus equi isolates obtained from horses by a multiplex PCR assay that amplifies the vap gene family (vapA, -B, -C, -D, -E, -F, -G, and -H). A total of 180 R. equi isolates were studied from four different sources, namely healthy horse feces (112), soil (12), stalls (23), and clinical isolates (33) from horse-breeding farms. The technique was performed and confirmed by sequencing of amplified vap gene family controls. Thirty-two (17.8%) of the R. equi isolates were positive for the vapA gene and carried at least three other vap genes. All 147 isolates from equine feces, stalls, and soil failed to demonstrate any genes associated with virulence-inducing proteins. About 32 (97.0%) out of the 33 clinical equine isolates tested positive for the multiplex PCR assay for the vap gene family. They demonstrated six molecular profiles: 100% featured the vapA, vapD, and vapG genes, 86.6% vapF, 76.6% vapH, 43.3% vapC, 36.6% vapE, and none vapB. The most frequent molecular profile was vap A, -D, -F, G, and -H, where this profile was present in 37.5% of the strains. Moreover, there was no molecular epidemiological pattern for R. equi isolates that uniquely mapped to each horse-breeding farm studied. Our proposed technique allows the identification of eight members of the vap gene family (vapA, B, -C, -D, -E, -F, -G, and -H). It is a practical and efficient method of conducting clinical and epidemiological studies on R. equi isolates.     

Comparative Analysis of the Full Genome of Helicobacter pylori Isolate Sahul64 Identifies Genes of High Divergence

Isolates of Helicobacter pylori can be classified phylogeographically. High genetic diversity and rapid microevolution are a hallmark of H. pylori genomes, a phenomenon that is proposed to play a functional role in persistence and colonization of diverse human populations. To provide further genomic evidence in the lineage of H. pylori and to further characterize diverse strains of this pathogen in different human populations, we report the finished genome sequence of Sahul64, an H. pylori strain isolated from an indigenous Australian. Our analysis identified genes that were highly divergent compared to the 38 publically available genomes, which include genes involved in the biosynthesis and modification of lipopolysaccharide, putative prophage genes, restriction modification components, and hypothetical genes. Furthermore, the virulence-associated vacA locus is a pseudogene and the cag pathogenicity island (cagPAI) is not present. However, the genome does contain a gene cluster associated with pathogenicity, including dupA. Our analysis found that with the addition of Sahul64 to the 38 genomes, the core genome content of H. pylori is reduced by approximately 14% (∼170 genes) and the pan-genome has expanded from 2,070 to 2,238 genes. We have identified three putative horizontally acquired regions, including one that is likely to have been acquired from the closely related Helicobacter cetorum prior to speciation. Our results suggest that Sahul64, with the absence of cagPAI, highly divergent cell envelope proteins, and a predicted nontransportable VacA protein, could be more highly adapted to ancient indigenous Australian people but with lower virulence potential compared to other sequenced and cagPAI-positive H. pylori strains.     

Internally Controlled Real-Time PCR Method for Quantitative Species-Specific Detection and vapA Genotyping of Rhodococcus equi

We developed a novel quantitative real-time PCR (Q-PCR) method for the soil actinomycete Rhodococcus equi, an important horse pathogen and emerging human pathogen. Species-specific quantification was achieved by targeting the chromosomal monocopy gene choE, universally conserved in R. equi. The choE Q-PCR included an internal amplification control (IAC) for identification of false negatives. A second Q-PCR targeted the virulence plasmid gene vapA, carried by most horse isolates but infrequently found in isolates from other sources. The choE-IAC and vapA assays were 100% sensitive and specific as determined using 178 R. equi isolates, 77 nontarget bacteria, and a panel of 60 R. equi isolates with known vapA⁺ and vapA-negative (including vapB⁺) plasmid genotypes. The vapA⁺ frequency among isolate types was as follows: horse, 85%; human, 20%; bovine and pig, 0%; others, 27%. The choE-IAC Q-PCR could detect up to one genome equivalent using R. equi DNA or 100 bacteria/ml using DNA extracted from artificially contaminated horse bronchoalveolar lavage (BAL) fluid. Quantification was linear over a 6-log dynamic range down to ≈10 target molecules (or 1,000 CFU/ml BAL fluid) with PCR efficiency E of >0.94. The vapA assay had similar performance but appeared unsuitable for accurate (vapA⁺) R. equi quantification due to variability in target gene or plasmid copy number (1 to 9). The dual-reaction Q-PCR system here reported offers a useful tool to both medical and veterinary diagnostic laboratories for the quantitative detection of R. equi and (optional) vapA⁺ “horse-pathogenic” genotype determination.     

A Global Overview of the Genetic and Functional Diversity in the Helicobacter pylori cag Pathogenicity Island

The Helicobacter pylori cag pathogenicity island (cagPAI) encodes a type IV secretion system. Humans infected with cagPAI–carrying H. pylori are at increased risk for sequelae such as gastric cancer. Housekeeping genes in H. pylori show considerable genetic diversity; but the diversity of virulence factors such as the cagPAI, which transports the bacterial oncogene CagA into host cells, has not been systematically investigated. Here we compared the complete cagPAI sequences for 38 representative isolates from all known H. pylori biogeographic populations. Their gene content and gene order were highly conserved. The phylogeny of most cagPAI genes was similar to that of housekeeping genes, indicating that the cagPAI was probably acquired only once by H. pylori, and its genetic diversity reflects the isolation by distance that has shaped this bacterial species since modern humans migrated out of Africa. Most isolates induced IL-8 release in gastric epithelial cells, indicating that the function of the Cag secretion system has been conserved despite some genetic rearrangements. More than one third of cagPAI genes, in particular those encoding cell-surface exposed proteins, showed signatures of diversifying (Darwinian) selection at more than 5% of codons. Several unknown gene products predicted to be under Darwinian selection are also likely to be secreted proteins (e.g. HP0522, HP0535). One of these, HP0535, is predicted to code for either a new secreted candidate effector protein or a protein which interacts with CagA because it contains two genetic lineages, similar to cagA. Our study provides a resource that can guide future research on the biological roles and host interactions of cagPAI proteins, including several whose function is still unknown.     

Distribution of pathogenicity island markers and virulence factors in new phylogenetic groups of uropathogenic <Emphasis Type="Italic">Escherichia coli</Emphasis> isolates

The present study was aimed at investigating the relationship between the new Clermont’s phylogenetic groups, virulence factors, and pathogenicity island markers (PAIs) among uropathogenic Escherichia coli (UPEC) in Iran. This cross-sectional study was carried out on 140 UPEC isolates collected from patients with urinary tract infections in Bushehr, Iran. All isolates were subjected to phylogenetic typing using a new quadruplex-PCR method. The presence of PAI markers and virulence factors in UPEC strains was evaluated by multiplex PCR. The most predominant virulence gene was fimH (85%), followed by iucC (61.4%), papC (38.6%), hlyA (22.1%), cnf-1 (18.6%), afa (10.7%), papG and neuC (each 9.3%), ibeA (3.6%), and sfa/foc (0.7%). The most common phylogenetic group was related to B2 (39.3%), and the least common to A (0.7%). The most prevalent PAI marker was PAI IV536 (77.14%), while markers for PAI III536 (13.57%), PAI IIJ96 (12.86%), and PAI II536 (12.14%) were the least frequent among the UPEC strains. Meanwhile, the PAI IJ96 marker was not detected. There was a significant association between the phylogenetic group B2 and all the studied virulence genes and PAI markers. To our knowledge, this is the first study to compare the relationship between new phylogenetic groups, virulence genes and PAI markers in UPEC strains in Iran. The phylogenetic group B2 was predominantly represented among the studied virulence genes and PAI markers, indicating the preference of particular strains to carry virulence genes.     

The Transition from a Phytopathogenic Smut Ancestor to an Anamorphic Biocontrol Agent Deciphered by Comparative Whole-Genome Analysis

Pseudozyma flocculosa is related to the model plant pathogen Ustilago maydis yet is not a phytopathogen but rather a biocontrol agent of powdery mildews; this relationship makes it unique for the study of the evolution of plant pathogenicity factors. The P. flocculosa genome of ∼23 Mb includes 6877 predicted protein coding genes. Genome features, including hallmarks of pathogenicity, are very similar in P. flocculosa and U. maydis, Sporisorium reilianum, and Ustilago hordei. Furthermore, P. flocculosa, a strict anamorph, revealed conserved and seemingly intact mating-type and meiosis loci typical of Ustilaginales. By contrast, we observed the loss of a specific subset of candidate secreted effector proteins reported to influence virulence in U. maydis as the singular divergence that could explain its nonpathogenic nature. These results suggest that P. flocculosa could have once been a virulent smut fungus that lost the specific effectors necessary for host compatibility. Interestingly, the biocontrol agent appears to have acquired genes encoding secreted proteins not found in the compared Ustilaginales, including necrosis-inducing-Phytophthora-protein- and Lysin-motif- containing proteins believed to have direct relevance to its lifestyle. The genome sequence should contribute to new insights into the subtle genetic differences that can lead to drastic changes in fungal pathogen lifestyles.     

A Real-Time Impedance Based Method to Assess Rhodococcus equi Virulence

Rhodococcus equi is a facultative intracellular pathogen of macrophages and the causative agent of foal pneumonia. R. equi virulence is usually assessed by analyzing intracellular growth in macrophages by enumeration of bacteria following cell lysis, which is time consuming and does not allow for a high throughput analysis. This paper describes the use of an impedance based real-time method to characterize proliferation of R. equi in macrophages, using virulent and attenuated strains lacking the vapA gene or virulence plasmid. Image analysis suggested that the time-dependent cell response profile (TCRP) is governed by cell size and roundness as well as cytoxicity of infecting R. equi strains. The amplitude and inflection point of the resulting TCRP were dependent on the multiplicity of infection as well as virulence of the infecting strain, thus distinguishing between virulent and attenuated strains.     

Plasmid Profiles of Virulent Rhodococcus equi Strains Isolated from Infected Foals in Poland

Rhodococcus equi is an important bacterial pathogen in foals up to 6 months old, widespread in horse farms all over the world. It was found that only virulent R. equi strains expressing 15–17 kDa virulence-associated protein (VapA) and having large virulence plasmid of 85–90 kb containing vapA gene are pathogenic for horses. To date, 12 plasmid types have been reported in VapA positive strains from horses. There are no data concerning plasmid types of Polish field R. equi strains isolated from horses and horse farm environment. The aim of the study is to determine plasmid profiles of virulent R. equi strains isolated in Poland from dead foals as well as from soil samples taken from horse breeding farms. Plasmid profiles of 10 clinical strains derived from 8 farms and 11 environmental strains from 3 farms, confirmed as virulent by PCR, were compared with 12 reference strains containing the known plasmid size and type. Plasmid DNAs were analysed by digestion with the restriction endonucleases BamHI, EcoRI, EcoT22I, and HindIII for detailed comparison and estimation of plasmid sizes. The results of RFLP analysis revealed that all except one isolates used in the study are classified as VapA 85 kb type I plasmid. One strain harboured VapA 87 kb type I plasmid. This is the first report of plasmid types of Polish field R. equi strains. The results of our preliminary investigations on horse farms located in central and eastern Poland indicate that the virulent R. equi strains thus far isolated from diseased foals and horse farms environment represent a highly uniform plasmid pattern.     

Protein interaction mapping on a functional shotgun sequence of Rickettsia sibirica

Protein interaction maps can reveal novel pathways and functional complexes, allowing ‘guilt by association’ annotation of uncharacterized proteins. To address the need for large-scale protein interaction analyses, a bacterial two-hybrid system was coupled with a whole genome shotgun sequencing approach for microbial genome analysis. We report the first large-scale proteomics study using this system, integrating de novo genome sequencing with functional interaction mapping and annotation in a high-throughput format. We apply the approach by shotgun sequencing and annotating the genome of Rickettsia sibirica strain 246, an obligate intracellular human pathogen among the Spotted Fever Group rickettsiae. The bacteria invade endothelial cells and cause lysis after large amounts of progeny have accumulated. Little is known about specific Rickettsial virulence factors and their mode of pathogenicity. Analysis of the combined genomic sequence and protein–protein interaction data for a set of virulence related Type IV secretion system (T4SS) proteins revealed over 250 interactions and will provide insight into the mechanism of Rickettsial pathogenicity.     

VapA of Rhodococcus equi binds phosphatidic acid

Rhodococcus equi is a multihost, facultative intracellular bacterial pathogen that primarily causes pneumonia in foals less than six months in age and immunocompromised people. Previous studies determined that the major virulence determinant of R. equi is the surface bound virulence associated protein A (VapA). The presence of VapA inhibits the maturation of R. equi‐containing phagosomes and promotes intracellular bacterial survival, as determined by the inability of vapA deletion mutants to replicate in host macrophages. While the mechanism of action of VapA remains elusive, we show that soluble recombinant VapA^32‐189 both rescues the intramacrophage replication defect of a wild type R. equi strain lacking the vapA gene and enhances the persistence of nonpathogenic Escherichia coli in macrophages. During macrophage infection, VapA was observed at both the bacterial surface and at the membrane of the host‐derived R. equi containing vacuole, thus providing an opportunity for VapA to interact with host constituents and promote alterations in phagolysosomal function. In support of the observed host membrane binding activity of VapA, we also found that rVapA^32‐189 interacted specifically with liposomes containing phosphatidic acid in vitro. Collectively, these data demonstrate a lipid binding property of VapA, which may be required for its function during intracellular infection.     

Additional Evidence That the Polymerase Subunits Contribute to the Viral Replication and the Virulence of H5N1 Avian Influenza Virus Isolates in Mice

Genetically similar H5N1 viruses circulating in the avian reservoir exhibit different levels of pathogenicity in mice. In this study, we characterized two highly pathogenic H5N1 avian isolates—A/Hunan/316/2005 (HN05), which is highly pathogenic in mice, and A/Hubei/489/2004 (HB04), which is nonpathogenic. In mammalian cells, HN05 replicates more efficiently than HB04, although both viruses have similar growth kinetics in avian cells. We used reverse genetics to generate recombinant H5N1 strains containing genes from HN05 and HB04 and examined their virulence. HN05 genes encoding the polymerase complex determine pathogenicity and viral replication ability both in vitro and in vivo. The PB2 subunit plays an important role in enhancing viral replication, and the PB1 and PA subunits contribute mainly to pathogenicity in mice. These results can be used to elucidate host-range expansion and the molecular basis of the high virulence of H5N1 viruses in mammalian species.     

Horizontal gene transfer from a flowering plant to the insular pine Pinus canariensis (Chr. Sm. Ex DC in Buch)

Horizontal gene transfer (HGT) is viewed as very common in the plant mitochondrial (mt) genome, but, to date, only one case of HGT has been found in gymnosperms. Here we report a new case of HGT, in which a mt nad5-1 fragment was transferred from an angiosperm to Pinus canariensis. Quantitative assay and sequence analyses showed that the foreign nad5-1 is located in the mt genome of P. canariensis and is nonfunctional. An extensive survey in the genus Pinus revealed that the angiosperm-derived nad5-1 is restricted to P. canariensis and present across the species'' range. Molecular dating based on chloroplast DNA suggested that the HGT event occurred in the late Miocene after P. canariensis split from its closest relatives, and that the foreign copy became fixed in P. canariensis owing to drift during its colonization of the Canary Islands. The mechanism of this HGT is unclear but it was probably achieved through either direct cell–cell contact or external vectors. Our discovery provides evidence for an important role of HGT in plant mt genome evolution.