Mathematical modeling of evolution of horizontally transferred genes

We describe a stochastic birth-and-death model of evolution of horizontally transferred genes in microbial populations. The model is a generalization of the stochastic model described by Berg and Kurland and includes five parameters: the rate of mutational inactivation, selection coefficient, invasion rate (i.e., rate of arrival of a novel sequence from outside of the recipient population), within-population horizontal transmission ("infection") rate, and population size. The model of Berg and Kurland included four parameters, namely, mutational inactivation, selection coefficient, population size, and "infection." However, the effect of "infection" was disregarded in the interpretation of the results, and the overall conclusion was that horizontally acquired sequences can be fixed in a population only when they confer a substantial selective advantage onto the recipient and therefore are subject to strong positive selection. Analysis of the present model in different domains of parameter values shows that, as long as the rate of within-population horizontal transmission is comparable to the mutational inactivation rate and there is even a low rate of invasion, horizontally acquired sequences can be fixed in the population or at least persist for a long time in a substantial fraction of individuals in the population even when they are neutral or slightly deleterious. The available biological data strongly suggest that intense within-population and even between-populations gene flows are realistic for at least some prokaryotic species and environments. Therefore, our modeling results are compatible with the notion of a pivotal role of horizontal gene transfer in the evolution of prokaryotes.     

Identification and categorization of horizontally transferred genes in prokaryotic genomes

Horizontal gene transfer （HGT）, a process through which genomes acquire genetic materials from distantly related organisms, is believed to be one of the major forces in prokaryotic genome evolution.However, systematic investigation is still scarce to clarify two basic issues about HGT： （1） what types of genes are transferred; and （2） what influence HGT events over the organization and evolution of biological pathways. Genome-scale investigations of these two issues will advance the systematical understanding of HGT in the context of prokaryotic genome evolution. Having investigated 82 genomes, we constructed an HGT database across broad evolutionary timescales. We identified four function categories containing a high proportion of horizontally transferred genes： cell envelope, energy metabolism, regulatory functions, and transport/binding proteins. Such biased function distribution indicates that HGT is not completely random;instead, it is under high selective pressure, required by function restraints in organisms. Furthermore, we mapped the transferred genes onto the connectivity structure map of organism-specific pathways listed in Kyoto Encyclopedia of Genes and Genomes （KEGG）. Our results suggest that recruitment of transferred genes into pathways is also selectively constrained because of the tuned interaction between original pathway members. Pathway organization structures still conserve well through evolution even with the recruitment of horizontally transferred genes. Interestingly, in pathways whose organization were significantly affected by HGT events, the operon-like arrangement of transferred genes was found to be prevalent. Such results suggest that operon plays an essential and directional role in the integration of alien genes into pathways.     

Limitations of compositional approach to identifying horizontally transferred genes

Genes with atypical G+C content and pattern of codon usage in a certain genome are possibly of exotic origin, and this idea has been applied to identify horizontal events. In this way, it was postulated that a total of 755 genes in the E. coli genome are relics of horizontal events after the divergence of E. coli from the Salmonella lineage 100 million years ago (Lawrence and Ochman, 1998). In this paper we propose a new way to study sequence composition more thoroughly. We found that although the 755 genes differ in composition from other genes in the E. coli genome, the difference is minor. If we accepted that these genes are horizontally transferred, then (1) it would be more likely that they were transferred from genomes evolutionarily closely related to E. coli; but (2) the dating method used by Lawrence and Ochman (1997, 1998) largely underestimated the average age of introduced sequences in the E. coli genome, in particular, most of the 755 genes should be introduced into E. coli before, instead of after, the divergence of E. coli from the Salmonella lineage. Our study reveals that atypical G+C content and pattern of codon usage are not reliable indicators of horizontal gene transfer events. Received: 27 September 2000 / Accepted: 9 April 2001     

A fragment of chloroplast DNA was transferred horizontally,probably from non-eudicots,to mitochondrial genome of <Emphasis Type="Italic">Phaseolus</Emphasis>

The mitochondrial genomes of some Phaseolus species contain a fragment of chloroplast trnA gene intron, named pvs-trnA for its location within the Phaseolus vulgaris sterility sequence (pvs). The purpose of this study was to determine the type of transfer (intracellular or horizontal) that gave rise to pvs-trnA. Using a PCR approach we could not find the respective portion of the trnA gene as a part of pvs outside the Phaseolus genus. However, a BLAST search revealed longer fragments of trnA present in the mitochondrial genomes of some Citrus species, Helianthus annuus and Zea mays. Basing on the identity or near-identity between these mitochondrial sequences and their chloroplast counterparts we concluded that they had relocated from chloroplasts to mitochondria via recent, independent, intracellular DNA transfers. In contrast, pvs-trnA displayed a relatively higher sequence divergence when compared with its chloroplast counterpart from Phaseolus vulgaris. Alignment of pvs-trnA with corresponding trnAfragments from 35 plant species as well as phylogenetic analysis revealed that pvs-trnA grouped with non-eudicot sequences and was well separated from all Fabalessequences. In conclusion, we propose that pvs-trnA arose via horizontal transfer of a trnA intron fragment from chloroplast of a non-eudicot plant to Phaseolus mitochondria. This is the first example of horizontal transfer of a chloroplast sequence to the mitochondrial genome in higher plants.     

Very little intron gain in Entamoeba histolytica genes laterally transferred from prokaryotes

The evolution of spliceosomal introns remains intensely debated. We studied 96 Entamoeba histolytica genes previously identified as having been laterally transferred from prokaryotes, which were presumably intronless at the time of transfer. Ninety out of the 96 are also present in the reptile parasite Entamoeba invadens, indicating lateral transfer before the species' divergence approximately 50 MYA. We find only 2 introns, both shared with E. invadens. Thus, no intron gains have occurred in approximately 50 Myr, implying a very low rate of intron gain of less than one gain per gene per approximately 4.5 billion years. Nine other predicted introns are due to annotation errors reflecting apparent mistakes in the E. histolytica genome assembly. These results underscore the massive differences in intron gain rates through evolution.     

Laterally transferred genomic islands in Xanthomonadales related to pathogenicity and primary metabolism.

Lateral gene transfer (LGT) is considered as one of the drivers in bacterial genome evolution, usually associated with increased fitness and/or changes in behavior, especially if one considers pathogenic vs. non-pathogenic bacterial groups. The genomes of two phytopathogens, Xanthomonas campestris pv. campestris and Xanthomonas axonopodis pv. citri, were previously inspected for genome islands originating from LGT events, and, in this work, potentially early and late LGT events were identified according to their altered nucleotide composition. The biological role of the islands was also assessed, and pathogenicity, virulence and secondary metabolism pathways were functions highly represented, especially in islands that were found to be recently transferred. However, old islands are composed of a high proportion of genes related to cell primary metabolic functions. These old islands, normally undetected by traditional atypical composition analysis, but confirmed as product of LGT by atypical phylogenetic reconstruction, reveal the role of LGT events by replacing core metabolic genes normally inherited by vertical processes.     

家蚕chi、gluE和fruA基因与微生物相应基因的同源性及基因水平转移初探

通过对家蚕(Bombyxm mori)大规模EST的分析,发现家蚕的chi、gluE和fruA基因分别与微生物的相应基因存在高度的氨基酸序列同源性,且进化关系很近,但与线虫(Caenorhabditis elegans)、果蝇(Drosophila melanogaster)、按蚊(Anopheles gamble)以及家蚕近缘昆虫的类似基因之间的相似性却非常低。这表明它们可能分别与微生物的同源基因具有共同的祖先,即微生物的基因水平转移给了家蚕,进化途径不属于垂直遗传。     

Identification of horizontally transferred genes in the genus Colletotrichum reveals a steady tempo of bacterial to fungal gene transfer

Background

Horizontal gene transfer (HGT) is the stable transmission of genetic material between organisms by means other than vertical inheritance. HGT has an important role in the evolution of prokaryotes but is relatively rare in eukaryotes. HGT has been shown to contribute to virulence in eukaryotic pathogens. We studied the importance of HGT in plant pathogenic fungi by identifying horizontally transferred genes in the genomes of three members of the genus Colletotrichum.

Results

We identified eleven HGT events from bacteria into members of the genus Colletotrichum or their ancestors. The HGT events include genes involved in amino acid, lipid and sugar metabolism as well as lytic enzymes. Additionally, the putative minimal dates of transference were calculated using a time calibrated phylogenetic tree. This analysis reveals a constant flux of genes from bacteria to fungi throughout the evolution of subphylum Pezizomycotina.

Conclusions

Genes that are typically transferred by HGT are those that are constantly subject to gene duplication and gene loss. The functions of some of these genes suggest roles in niche adaptation and virulence. We found no evidence of a burst of HGT events coinciding with major geological events. In contrast, HGT appears to be a constant, albeit rare phenomenon in the Pezizomycotina, occurring at a steady rate during their evolution.

Electronic supplementary material

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

A computational tool for the genomic identification of regions of unusual compositional properties and its utilization in the detection of horizontally transferred sequences

Similarity Plot (S-plot) is a Windows-based application for large-scale comparisons and 2-dimensional visualization of compositional similarities between genomic sequences. This application combines 2 approaches widely used in genomics: window analysis of statistical characteristics along genomes and dot-plot visual representation. S-plot is effective in identifying highly similar regions between genomes as well as regions with unusual compositional properties (RUCPs) within a single genome, which may be indicative of horizontal gene transfer or of locus-specific selective forces. We use S-plot to identify regions that may have originated through horizontal gene transfer through a 2-step approach, by first comparing a genomic sequence to itself and, subsequently, comparing it to the genomic sequence of a closely related taxon. Moreover, by comparing these suspect sequences to one another, we can estimate a minimum number of sources for these putative xenologous sequences. We illustrate the uses of S-plot in a comparison involving Escherichia coli K12 and E. coli O157:H7. In O157:H7, we found 145 regions that have most probably originated through horizontal gene transfer. By using S-plot to compare each of these regions with 277 completely sequenced prokaryotic genomes, 1 sequence was found to have similar compositional properties to the Yersinia pseudotuberculosis genome, indicating a transfer from a Yersinia or Yersinia relative. Based upon our analysis of RUCPs in O157:H7, we infer that there were at least 53 sources of horizontally transferred sequences.     

Direct gene transfer to protoplasts: Fate of the transferred genes

Direct gene transfer to protoplasts is one of several methods developed for the production of transgenic plants. This method utilizes the efficient uptake of DNA from the surrounding medium by protoplasts (cell wall-less plant cells). Where a suitable protoplast system exists large numbers of transformant clones can be efficiently produced and often regenerated to normal fertile plants. This review concentrates on the fate of the DNA which is taken up into the protoplasts. Particular emphasis is given to the factors which can influence the integration and form of the transferred DNA, the expression of transferred genes, and the inheritance in further generations of those genes. The information available suggests (1) that DNA is taken up by a large proportion of the cells in a transformation mixture, (2) that this DNA forms complexes sometimes involving carrier DNA, (3) that fewer cells actually take up DNA into the nucleus, and (4) that the complex may be rearranged and/or amplified and then integrated into the genome. If the DNA is arranged in such a way that a gene can be expressed it does so in a normal manner and is stably inherited both mitotically and meiotically.     

An integrative method for identifying the over-annotated protein-coding genes in microbial genomes

The falsely annotated protein-coding genes have been deemed one of the major causes accounting for the annotating errors in public databases. Although many filtering approaches have been designed for the over-annotated protein-coding genes, some are questionable due to the resultant increase in false negative. Furthermore, there is no webserver or software specifically devised for the problem of over-annotation. In this study, we propose an integrative algorithm for detecting the over-annotated protein-coding genes in microorganisms. Overall, an average accuracy of 99.94% is achieved over 61 microbial genomes. The extremely high accuracy indicates that the presented algorithm is efficient to differentiate the protein-coding genes from the non-coding open reading frames. Abundant analyses show that the predicting results are reliable and the integrative algorithm is robust and convenient. Our analysis also indicates that the over-annotated protein-coding genes can cause the false positive of horizontal gene transfers detection. The webserver of the proposed algorithm can be freely accessible from www.cbi.seu.edu.cn/RPGM.     

A likelihood framework to measure horizontal gene transfer

We suggest a likelihood-based approach to estimate an overall rate of horizontal gene transfer (HGT) in a simplified setting. To this end, we assume that the number of occurring HGT events within a given time interval follows a Poisson process. To obtain estimates for the rate of HGT, we simulate the distribution of tree topologies for different numbers of HGT events on a clocklike species tree. Using these simulated distributions, we estimate an HGT rate for a collection of gene trees representing a set of taxa. As an illustrative example, we use the "Clusters of Orthologous Groups of proteins" (COGs). We also perform a correction of the estimated rate taking into account the inaccuracies due to gene tree reconstructions. The results suggest a corrected HGT rate of about 0.36 per gene and unit time, in other words 11 HGT events have occurred on average among the 44 taxa of the COG species tree. A software package to estimate an HGT rate is available online (http://www.cibiv.at/software/hgt/).     

Genome-wide identification of novel genomic islands that contribute to Salmonella virulence in mouse systemic infection

Salmonella pathogenicity islands are inserted into the genome by horizontal gene transfer and are required for expression of full virulence. Here, we performed tRNA scanning of the genome of Salmonella enterica serovar Typhimurium and compared it with that of nonpathogenic Escherichia coli in order to identify genomic islands that contribute to Salmonella virulence. Using deletion analysis, we identified four genomic islands that are required for virulence in the mouse infection model. One of the newly identified pathogenicity islands was the pheV- tRNA-located genomic island, which is comprised of 26 126 bp, and encodes 22 putative genes, including STM3117–STM3138. We also showed that the pheV tRNA-located genomic island is widely distributed among different nontyphoid Salmonella serovars. Furthermore, genes including STM3118–STM3121 were identified as novel virulence-associated genes within the pheV- tRNA-located genomic island. These results indicate that a Salmonella -specific pheV- tRNA genomic island is involved in Salmonella pathogenesis among the nontyphoid Salmonella serovars.     

Successful lateral transfer requires codon usage compatibility between foreign genes and recipient genomes

We present evidence supporting the notion that codon usage (CU) compatibility between foreign genes and recipient genomes is an important prerequisite to assess the selective advantage of imported functions, and therefore to increase the fixation probability of horizontal gene transfer (HGT) events. This contrasts with the current tendency in research to predict recent HGTs in prokaryotes by assuming that acquired genes generally display poor CU. By looking at the CU level (poor, typical, or rich) exhibited by putative xenologs still resembling their original CU, we found that most alien genes predominantly present typical CU immediately upon introgression, thereby suggesting that the role of CU amelioration in HGT has been overemphasized. In our strategy, we first scanned a representative set of 103 complete prokaryotic genomes for all pairs of candidate xenologs (exported/imported genes) displaying similar CU. We applied additional filtering criteria, including phylogenetic validations, to enhance the reliability of our predictions. Our approach makes no assumptions about the CU of foreign genes being typical or atypical within the recipient genome, thus providing a novel unbiased framework to study the evolutionary dynamics of HGT.     

Frequent,independent transfers of a catabolic gene from bacteria to contrasted filamentous eukaryotes

Even genetically distant prokaryotes can exchange genes between them, and these horizontal gene transfer events play a central role in adaptation and evolution. While this was long thought to be restricted to prokaryotes, certain eukaryotes have acquired genes of bacterial origin. However, gene acquisitions in eukaryotes are thought to be much less important in magnitude than in prokaryotes. Here, we describe the complex evolutionary history of a bacterial catabolic gene that has been transferred repeatedly from different bacterial phyla to stramenopiles and fungi. Indeed, phylogenomic analysis pointed to multiple acquisitions of the gene in these filamentous eukaryotes—as many as 15 different events for 65 microeukaryotes. Furthermore, once transferred, this gene acquired introns and was found expressed in mRNA databases for most recipients. Our results show that effective inter-domain transfers and subsequent adaptation of a prokaryotic gene in eukaryotic cells can happen at an unprecedented magnitude.     

Evidence for transfer of antibiotic-resistance genes in soil populations of streptomycetes

Phylogenetic analysis was used to evaluate the hypothesis of gene transfer in streptomycetes, many of which are antibiotic producers. The diversity and possible origins of streptomycin-resistance genes was investigated for a population of Streptomyces strains isolated from a site in Brazil where antibiotic production had previously been implicated. The analysis provides compelling evidence for the transfer of these genes. Examination of other Streptomyces -type strains also reveals a scattered distribution of streptomycin producers with respect to the overall phylogeny. These results suggest that horizontal gene transfer may be an important factor in the evolution of antibiotic genes in streptomycetes.     

Complex histories of genes encoding 3-hydroxy-3-methylglutaryl-CoenzymeA reductase

The mevalonate pathway for the synthesis of isoprenoids can be found in organisms from all domains of life. It has been previously demonstrated that the first gene specific to that pathway, which encodes the enzyme 3-hydroxy-3-methylglutaryl-CoenzymeA reductase (HMGR), has been transferred between domains by lateral gene transfer on several occasions. Here we look within the domain Bacteria at lateral acquisition of HMGR, whether as a single gene or as part of a mevalonate pathway cluster. We observe a complex history of multiple transfer events probably reflecting the fact that HMGR could be beneficial in a variety of physiological and genetic contexts. We demonstrate that even in Vibrio species, where HMGR is not clustered with other genes to form an operon or a metabolic cluster, it is under strong purifying selection.