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.     

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.     

致病菌基因组岛的研究进展

细菌基因组岛是细菌基因组上的特定区域,和水平基因转移相关,具有一定的结构特点,常携带致病、耐药及与适应性等功能相关的基因。通过基因组岛在细菌间的移动,可以造成相关基因在细菌间的传播,在细菌生存和致病等过程中具有重要作用。目前已经可通过生物信息和分子生物学实验等方法对基因组岛进行预测和验证。通过对致病菌基因组岛的研究,可以阐释细菌致病性和耐药等重要功能的获得,对疾病进行溯源,在传染病预防控制中具有重要意义。     

The recombinase IntA is required for excision of esp-containing ICEEfm1 in Enterococcus faecium

Comparative genome analysis of Enterococcus faecium recently revealed that a genomic island containing the esp gene, referred to as the esp-containing pathogenicity island (esp PAI), can be transferred by conjugation and contains a partial Tn916-like element and an integrase gene, intA. Here, we characterize the role of intA in the excision of the esp PAI. An intA insertion-deletion mutant in E. faecium E1162 (E1162ΔintA) was constructed and in trans complemented with wild-type intA (E1162ΔintA::pEF30). Circular intermediates (CI) of excised esp PAI were determined using inverse PCR analysis on purified chromosomal DNA from strains E1162, E1162Δesp, E1162ΔintA, and E1162ΔintA::pEF30. In E1162 and E1162Δesp, CI of the esp PAI were detected. No CI were detected in E1162ΔintA, while in the complemented strain E1162ΔintA::pEF30 CI formation was restored, indicating that intA is essential for excision and subsequent mobilization of the esp-containing genomic island in E. faecium. Based on the fact that this island can be mobilized and is self-transmissible, we propose to change the name of the esp PAI to ICEEfm1.     

Identification of compositionally distinct regions in genomes using the centroid method

MOTIVATION: It is known that most genomic regions of special interest, e.g. horizontally acquired sequences, genomic islands, etc. have distinct word (m-mer) compositions. Most of the earlier work along this direction, addressed di- and tri-nucleotide compositions. We present an approach that can be applied to analyze compositions of any given word size. The method, called the centroid approach, can reveal compositionally distinct regions in genomic sequences for any given word size. RESULTS: We applied our method to 50 bacterial genomes and demonstrated its ability to identify embedded sequences of varying lengths from distantly related organisms. We also investigated the genetic makeup of the regions identified as compositionally distinct by our method, for four organisms from our dataset. Pathogenicity island (PAI) components and genes encoding strain-specific proteins are all frequently seen to be constituents of these regions. AVAILABILITY: Program is available on request from the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Tests for gene clustering.

Comparing chromosomal gene order in two or more related species is an important approach to studying the forces that guide genome organization and evolution. Linked clusters of similar genes found in related genomes are often used to support arguments of evolutionary relatedness or functional selection. However, as the gene order and the gene complement of sister genomes diverge progressively due to large scale rearrangements, horizontal gene transfer, gene duplication and gene loss, it becomes increasingly difficult to determine whether observed similarities in local genomic structure are indeed remnants of common ancestral gene order, or are merely coincidences. A rigorous comparative genomics requires principled methods for distinguishing chance commonalities, within or between genomes, from genuine historical or functional relationships. In this paper, we construct tests for significant groupings against null hypotheses of random gene order, taking incomplete clusters, multiple genomes, and gene families into account. We consider both the significance of individual clusters of prespecified genes and the overall degree of clustering in whole genomes.     

细菌基因组进化的分子策略

论述了细菌基因组进化的 4个分子策略 :点突变 ,基因组内重排 ,基因水平转移 ,基因缺失。从经典的达尔文进化论角度探讨了细菌基因组进化与表型进化的关系。     

PIPS: pathogenicity island prediction software

The adaptability of pathogenic bacteria to hosts is influenced by the genomic plasticity of the bacteria, which can be increased by such mechanisms as horizontal gene transfer. Pathogenicity islands play a major role in this type of gene transfer because they are large, horizontally acquired regions that harbor clusters of virulence genes that mediate the adhesion, colonization, invasion, immune system evasion, and toxigenic properties of the acceptor organism. Currently, pathogenicity islands are mainly identified in silico based on various characteristic features: (1) deviations in codon usage, G+C content or dinucleotide frequency and (2) insertion sequences and/or tRNA genetic flanking regions together with transposase coding genes. Several computational techniques for identifying pathogenicity islands exist. However, most of these techniques are only directed at the detection of horizontally transferred genes and/or the absence of certain genomic regions of the pathogenic bacterium in closely related non-pathogenic species. Here, we present a novel software suite designed for the prediction of pathogenicity islands (pathogenicity island prediction software, or PIPS). In contrast to other existing tools, our approach is capable of utilizing multiple features for pathogenicity island detection in an integrative manner. We show that PIPS provides better accuracy than other available software packages. As an example, we used PIPS to study the veterinary pathogen Corynebacterium pseudotuberculosis, in which we identified seven putative pathogenicity islands.     

Horizontal transfer of the plant virulence gene, nec1, and flanking sequences among genetically distinct Streptomyces strains in the Diastatochromogenes cluster

Evidence for the horizontal transfer of a pathogenicity island (PAI) carrying the virulence gene nec1 and flanking sequences among Streptomyces strains in the Diastatochromogenes cluster is presented. Plant-pathogenic, thaxtomin-producing Streptomyces strains, previously classified as S. scabiei based on the conventionally used phenotypic characteristics, were found to be genetically distinct from the type strain of S. scabiei based on DNA relatedness and 16S rDNA sequence analysis. Pairwise DNA-DNA hybridizations between some of these strains and the S. scabiei type strain were as low as 36%, a value much below what is conventionally accepted for species identity (70%). The sequence of the nec1 gene, however, was identical in all the S. scabiei and S. scabiei-like strains tested, irrespective of their DNA relatedness to the type strain of S. scabiei, their geographic origin, or the isolation host. Furthermore, a 26-kb DNA fragment including and flanking nec1 was also conserved among these strains based on restriction and Southern analyses. These data indicate that the etiology of potato scab is more complex than previously recognized; this result has important implications for potato scab management strategies. Previous research has suggested that horizontal transfer of a PAI was the mechanism for evolution of pathogenicity in S. acidiscabies and S. turgidiscabies, species that lie outside of the Diastatochromogenes cluster. Data presented here support this model and indicate that PAI transfer also has occurred frequently in species closely related to S. scabiei.     

A systematic method to identify genomic islands and its applications in analyzing the genomes of Corynebacterium glutamicum and Vibrio vulnificus CMCP6 chromosome I

MOTIVATION: Some genomic islands contain horizontally transferred genes, which play critical roles in altering the genotypes and phenotypes of organisms, and horizontal gene transfer has been recognized as a universal event throughout bacterial evolution. A windowless method to display the distribution of genomic GC content, the cumulative GC profile, is proposed to identify genomic islands in genomes whose complete genome sequences are available. Two new indices are proposed to assess the codon usage bias and amino acid usage bias in genomic islands. RESULTS: A 211 kb genomic island (CGGI-1) has been identified in the genome of Corynebacterium glutamicum, and three genomic islands VVGI-1, VVGI-2 and VVGI-3, with lengths 167, 40 and 33 kb, respectively, have been identified in the genome of Vibrio vulnificus CMCP6 chromosome I. The CGGI-1 is flanked by two approximately 500 bp direct repeats, and utilizes a Val-tRNA as the integration site. For the VVGI-1 and VVGI-2, each has an integrase gene at 5' junction. All the identified genomic islands show unusual GC content, codon usage and amino acid usage, compared with the rest of the genomes. In addition, it is found that genomic islands are fairly homogenous in terms of GC content variation. An index, h, to quantify the homogeneity of GC content for genomic islands is proposed, and it is shown that h is less than 0.1 for all the genomic islands analyzed. The cumulative GC profile, as well as various indices to assess the codon usage bias, amino acid usage bias and homogeneity of the genomic islands, will be useful in the analysis of other genomes. AVAILABILITY: Programs used in this work and numerical results are available upon request.     

The novel kasugamycin 2'-N-acetyltransferase gene aac(2')-IIa, carried by the IncP island, confers kasugamycin resistance to rice-pathogenic bacteria

Kasugamycin (KSM), a unique aminoglycoside antibiotic, has been used in agriculture for many years to control not only rice blast caused by the fungus Magnaporthe grisea but also rice bacterial grain and seedling rot or rice bacterial brown stripe caused by Burkholderia glumae or Acidovorax avenae subsp. avenae, respectively. Since both bacterial pathogens are seed-borne and cause serious injury to rice seedlings, the emergence of KSM-resistant B. glumae and A. avenae isolates highlights the urgent need to understand the mechanism of resistance to KSM. Here, we identified a novel gene, aac(2')-IIa, encoding a KSM 2'-N-acetyltransferase from both KSM-resistant pathogens but not from KSM-sensitive bacteria. AAC(2')-IIa inactivates KSM, although it reveals no cross-resistance to other aminoglycosides. The aac(2')-IIa gene from B. glumae strain 5091 was identified within the IncP genomic island inserted into the bacterial chromosome, indicating the acquisition of this gene by horizontal gene transfer. Although excision activity of the IncP island and conjugational gene transfer was not detected under the conditions tested, circular intermediates containing the aac(2')-IIa gene were detected. These results indicate that the aac(2')-IIa gene had been integrated into the IncP island of a donor bacterial species. Molecular detection of the aac(2')-IIa gene could distinguish whether isolates are resistant or susceptible to KSM. This may contribute to the production of uninfected rice seeds and lead to the effective control of these pathogens by KSM.     

A new experimental approach for studying bacterial genomic island evolution identifies island genes with bacterial host-specific expression patterns

Background  

Genomic islands are regions of bacterial genomes that have been acquired by horizontal transfer and often contain blocks of genes that function together for specific processes. Recently, it has become clear that the impact of genomic islands on the evolution of different bacterial species is significant and represents a major force in establishing bacterial genomic variation. However, the study of genomic island evolution has been mostly performed at the sequence level using computer software or hybridization analysis to compare different bacterial genomic sequences. We describe here a novel experimental approach to study the evolution of species-specific bacterial genomic islands that identifies island genes that have evolved in such a way that they are differentially-expressed depending on the bacterial host background into which they are transferred.     

Evidence for natural horizontal transfer of tetQ between bacteria that normally colonize humans and bacteria that normally colonize livestock.

Though numerous studies have shown that gene transfer occurs between distantly related bacterial genera under laboratory conditions, the frequency and breadth of horizontal transfer events in nature remain unknown. Previous evidence for natural intergeneric transfers came from studies of genes in human pathogens, bacteria that colonize the same host. We present evidence that natural transfer of a tetracycline resistance gene, tetQ, has occurred between bacterial genera that normally colonize different hosts. A DNA sequence comparative approach was taken to examine the extent of horizontal tetQ dissemination between species of Bacteroides, the predominant genus of the human colonic microflora, and between species of Bacteroides and of the distantly related genus Prevotella, a predominant genus of the microflora of the rumens and intestinal tracts of farm animals. Virtually identical tetQ sequences were found in a number of isolate pairs differing in taxonomy and geographic origin, indicating that extensive natural gene transmission has occurred. Among the exchange events indicated by the evidence was the very recent transfer of an allele of tetQ usually found in Prevotella spp. to a Bacteroides fragilis strain.     

Cell evolution: gene transfer agents and the origin of mitochondria

Recently, α-proteobacteria have been shown to possess virus-like gene transfer agents that facilitate high frequency gene transfer in natural environments between distantly related lineages. This system could have driven the genomic integration of the mitochondrial progenitor and its proto-eukaryote host and contributed to the evolutionary mosaic of genes seen in modern-day prokaryotic and eukaryotic genomes.     

病原菌毒力岛研究进展

毒力岛作为基因组岛的一种亚类,是细菌染色体上具有特定结构和功能特征的可移动基因大片段,经基因水平转移（转导、接合或转化）获得,可使细菌基因组进化在短期内发生“量的飞跃”,直接或间接增强细菌的生态适应性,与病原菌的致病性密切相关。毒力岛存在于多种动植物病原细菌中,对于细菌的毒力变异、遗传进化甚至新病原亚种形成有重要意义。简要综述了病原菌毒力岛的研究进展,介绍了毒力岛的结构、功能特征及其在病原菌进化中作用。     

Optimal gene partition into operons correlates with gene functional order

Gene arrangement into operons varies between bacterial species. Genes in a given system can be on one operon in some organisms and on several operons in other organisms. Existing theories explain why genes that work together should be on the same operon, since this allows for advantageous lateral gene transfer and accurate stoichiometry. But what causes the frequent separation into multiple operons of co-regulated genes that act together in a pathway? Here we suggest that separation is due to benefits made possible by differential regulation of each operon. We present a simple mathematical model for the optimal distribution of genes into operons based on a balance of the cost of operons and the benefit of regulation that provides 'just-when-needed' temporal order. The analysis predicts that genes are arranged such that genes on the same operon do not skip functional steps in the pathway. This prediction is supported by genomic data from 137 bacterial genomes. Our work suggests that gene arrangement is not only the result of random historical drift, genome re-arrangement and gene transfer, but has elements that are solutions of an evolutionary optimization problem. Thus gene functional order may be inferred by analyzing the operon structure across different genomes.     

Genomic fingerprints produced by PCR with consensus tRNA gene primers.

The polymerase chain reaction using only a single 'consensus' tRNA gene primer, or a pair of primers facing outward from tRNA genes, amplifies a set of DNA fragments in bacterial, plant and animal genomic DNAs. Presumably, these PCR fingerprints are mainly derived from the regions between closely linked tRNA genes. The pattern of the PCR products is determined by which genomes and which primer(s) are used. Genomic fingerprints are largely conserved within a species and, in bacteria, most products in the fingerprint are conserved between closely related species. Thus, PCR with tRNA gene consensus primers helps to identify species and genera.     

Role of Intraspecies Recombination in the Spread of Pathogenicity Islands within the Escherichia coli Species

Horizontal gene transfer is a key step in the evolution of bacterial pathogens. Besides phages and plasmids, pathogenicity islands (PAIs) are subjected to horizontal transfer. The transfer mechanisms of PAIs within a certain bacterial species or between different species are still not well understood. This study is focused on the High-Pathogenicity Island (HPI), which is a PAI widely spread among extraintestinal pathogenic Escherichia coli and serves as a model for horizontal transfer of PAIs in general. We applied a phylogenetic approach using multilocus sequence typing on HPI-positive and -negative natural E. coli isolates representative of the species diversity to infer the mechanism of horizontal HPI transfer within the E. coli species. In each strain, the partial nucleotide sequences of 6 HPI–encoded genes and 6 housekeeping genes of the genomic backbone, as well as DNA fragments immediately upstream and downstream of the HPI were compared. This revealed that the HPI is not solely vertically transmitted, but that recombination of large DNA fragments beyond the HPI plays a major role in the spread of the HPI within E. coli species. In support of the results of the phylogenetic analyses, we experimentally demonstrated that HPI can be transferred between different E. coli strains by F-plasmid mediated mobilization. Sequencing of the chromosomal DNA regions immediately upstream and downstream of the HPI in the recipient strain indicated that the HPI was transferred and integrated together with HPI–flanking DNA regions of the donor strain. The results of this study demonstrate for the first time that conjugative transfer and homologous DNA recombination play a major role in horizontal transfer of a pathogenicity island within the species E. coli.     

Discordant phylogenies within the rrn loci of Rhizobia

It is evident from complete genome sequencing results that lateral gene transfer and recombination are essential components in the evolutionary process of bacterial genomes. Since this has important implications for bacterial systematics, the primary objective of this study was to compare estimated evolutionary relationships among a representative set of alpha-Proteobacteria by sequencing analysis of three loci within their rrn operons. Tree topologies generated with 16S rRNA gene sequences were significantly different from corresponding trees assembled with 23S rRNA gene and internally transcribed space region sequences. Besides the incongruence in tree topologies, evidence that distinct segments along the 16S rRNA gene sequences of bacteria currently classified within the genera Bradyrhizobium, Mesorhizobium and Sinorhizobium have a reticulate evolutionary history was also obtained. Our data have important implications for bacterial taxonomy, because currently most taxonomic decisions are based on comparative 16S rRNA gene sequence analysis. Since phylogenetic placement based on 16S rRNA gene sequence divergence perhaps is questionable, we suggest that the proposals of bacterial nomenclature or changes in their taxonomy that have been made may not necessarily be warranted. Accordingly, a more conservative approach should be taken in the future, in which taxonomic decisions are based on the analysis of a wider variety of loci and comparative analytical methods are used to estimate phylogenetic relationships among the genomes under consideration.     

Analysis of genomic characters reveals that four distinct gene clusters are correlated with different functions in Burkholderia cenocepacia AU 1054

Possessing three circular chromosomes is a distinct genomic characteristic of Burkholderia cenocepacia AU 1054, a clinically important pathogen in cystic fibrosis. In this study, base composition, codon usage and functional role category were analyzed in the B. cenocepacia AU 1054 genome. Although no bias in the base and codon usage was detected between any two chromosomes, function differences did exist in the genes of each chromosome. Similar base composition and differential functional role categories indicated that genes on these three chromosomes were relatively stable and that a proper division of labor was established. Based on variations in the base or codon usage, four small gene clusters were observed in all of the genes. Multivariate analysis revealed that protein hydrophobicity played a predominant role in shaping base usage bias, while horizontal gene transfer and the gene expression level were the two most important factors that affected the codon usage bias. Interestingly, we also found that these gene clusters were correlated with different biological functions: (i) 45 pyrimidine-leading-codon preferred genes were predominantly involved in regulatory function; (ii) most drug resistance-related genes involved in 826 genes that coding for hydrophobic proteins; (iii) most of the 111 horizontal transfer genes were responsible for genomic plasticity; and (iv) 73 highly expressed genes (predicted by their codon adaptation index values) showed environmental adaptation to cystic fibrosis. Our results showed that genes with base or codon usage bias were affected by mutational pressure and natural selection, and their functions could contribute to drug assistance and transmissible activity in B. cenocepacia.