Duplication is more common among laterally transferred genes than among indigenous genes

Horizontally transferred glycoside hydrolase 26 may aid hemipteran insects in plant tissue digestion

Glycoside hydrolases are enzymes that break down complex carbohydrates into simple sugars by catalyzing the hydrolysis of glycosidic bonds. There have been multiple instances of adaptive horizontal gene transfer of genes belonging to various glycoside hydrolase families from microbes to insects, as glycoside hydrolases can metabolize constituents of the carbohydrate-rich plant cell wall. In this study, we characterize the horizontal transfer of a gene from the glycoside hydrolase family 26 (GH26) from bacteria to insects of the order Hemiptera. Our phylogenies trace the horizontal gene transfer to the common ancestor of the superfamilies Pentatomoidea and Lygaeoidea, which include stink bugs and seed bugs. After horizontal transfer, the gene was assimilated into the insect genome as indicated by the gain of an intron, and a eukaryotic signal peptide. Subsequently, the gene has undergone independent losses and expansions in copy number in multiple lineages, suggesting an adaptive role of GH26s in some insects. Finally, we measured tissue-level gene expression of multiple stink bugs and the large milkweed bug using publicly available RNA-seq datasets. We found that the GH26 genes are highly expressed in tissues associated with plant digestion, especially in the principal salivary glands of the stink bugs. Our results are consistent with the hypothesis that this horizontally transferred GH26 was co-opted by the insect to aid in plant tissue digestion and that this HGT event was likely adaptive.

Duplicate Gene Expression and Possible Mechanisms of Paralog Retention During Bacterial Genome Expansion

Gene duplication contributes to the evolution of expression and the origin of new genes, but the relative importance of different patterns of duplicate gene expression and mechanisms of retention remains debated and particularly poorly understood in bacteria. Here, we investigated gene expression patterns for two lab strains of the cyanobacterium Acaryochloris marina with expanding genomes that contain about 10-fold more gene duplicates compared with most bacteria. Strikingly, we observed a generally stoichiometric pattern of greater combined duplicate transcript dosage with increased gene copy number, in contrast to the prevalence of expression reduction reported for many eukaryotes. We conclude that increased transcript dosage is likely an important mechanism of initial duplicate retention in these bacteria and may persist over long periods of evolutionary time. However, we also observed that paralog expression can diverge rapidly, including possible functional partitioning, for which different copies were respectively more highly expressed in at least one condition. Divergence may be promoted by the physical separation of most Acaryochloris duplicates on different genetic elements. In addition, expression pattern for ancestrally shared duplicates could differ between strains, emphasizing that duplicate expression fate need not be deterministic. We further observed evidence for context-dependent transcript dosage, where the aggregate expression of duplicates was either greater or lower than their single-copy homolog depending on physiological state. Finally, we illustrate how these different expression patterns of duplicated genes impact Acaryochloris biology for the innovation of a novel light-harvesting apparatus and for the regulation of recA paralogs in response to environmental change.

Differential Expression of Two Copies of the irmA Gene in the Enteroaggregative E. coli Strain 042

Gene duplications significantly impact the gene repertoires of both eukaryotic and prokaryotic microorganisms. The genomes of pathogenic Escherichia coli strains share a group of duplicated genes whose function is mostly unknown. The irmA gene is one of the duplicates encoded in several pathogenic E. coli strains. The function of its gene product was investigated in the uropathogenic E. coli strain CFT073, which contains a single functional copy. The IrmA protein structure mimics that of human interleukin receptors and likely plays a role during infection. The enteroaggregative E. coli strain 042 contains two functional copies of the irmA gene. In the present work, we investigated their biological roles. The irmA_4509 allele is expressed under several growth conditions. Its expression is modulated by the global regulators OxyR and Hha, with optimal expression at 37°C and under nutritional stress conditions. Expression of the irmA_2244 allele can only be detected when the irmA_4509 allele is knocked out. Differences in the promoter regions of both alleles account for their differential expression. Our results show that under several environmental conditions, the expression of the IrmA protein in strain 042 is dictated by the irmA_4509 allele. The irmA_2244 allele appears to play a backup role to ensure IrmA expression when the irmA_4509 allele loses its function. IMPORTANCE Gene duplications occur in prokaryotic genomes at a detectable frequency. In many instances, the biological function of the duplicates is unknown, and hence, the significance of the presence of multiple copies of these genes remains unclear. In pathogenic E. coli isolates, the irmA gene can be present either as a single copy or in two or more copies. We focused our work on studying why a different pathogenic E. coli strain encodes two functional copies of the irmA gene. We show that under several environmental conditions, one of the alleles dictates IrmA expression, and the second remains silent. The latter allele is only expressed when the former is silenced. The presence of more than one functional copy of the irmA gene in some pathogenic E. coli strains can result in sufficient expression of this virulence factor during the infection process.

Gene Duplications in the Genomes of Staphylococci and Enterococci

Gene duplications are a feature of bacterial genomes. In the present work we analyze the extent of gene duplications in the genomes of three microorganisms that belong to the Firmicutes phylum and that are etiologic agents of several nosocomial infections: Staphylococcus aureus, Enterococcus faecium, and Enterococcus faecalis. In all three groups, there is an irregular distribution of duplications in the genomes of the strains analyzed. Whereas in some of the strains duplications are scarce, hundreds of duplications are present in others. In all three species, mobile DNA accounts for a large percentage of the duplicated genes: phage DNA in S. aureus, and plasmid DNA in the enterococci. Duplicates also include core genes. In all three species, a reduced group of genes is duplicated in all strains analyzed. Duplication of the deoC and rpmG genes is a hallmark of S. aureus genomes. Duplication of the gene encoding the PTS IIB subunit is detected in all enterococci genomes. In E. faecalis it is remarkable that the genomes of some strains encode duplicates of the prgB and prgU genes. They belong to the prgABCU cluster, which responds to the presence of the peptide pheromone cCF10 by expressing the surface adhesins PrgA, PrgB, and PrgC.

How enzyme promiscuity and horizontal gene transfer contribute to metabolic innovation

Promiscuity is the coincidental ability of an enzyme to catalyze its native reaction and additional reactions that are not biological functions in the same active site. Promiscuity plays a central role in enzyme evolution and is thus a useful property for protein and metabolic engineering. This review examines enzyme evolution holistically, beginning with evaluating biochemical support for four enzyme evolution models. As expected, there is strong biochemical support for the subfunctionalization and innovation-amplification-divergence models, in which promiscuity is a central feature. In many cases, however, enzyme evolution is more complex than the models indicate, suggesting much is yet to be learned about selective pressures on enzyme function. A complete understanding of enzyme evolution must also explain the ability of metabolic networks to integrate new enzyme activities. Hidden within metabolic networks are underground metabolic pathways constructed from promiscuous activities. We discuss efforts to determine the diversity and pervasiveness of underground metabolism. Remarkably, several studies have discovered that some metabolic defects can be repaired via multiple underground routes. In prokaryotes, metabolic innovation is driven by connecting enzymes acquired by horizontal gene transfer (HGT) into the metabolic network. Thus, we end the review by discussing how the combination of promiscuity and HGT contribute to evolution of metabolism in prokaryotes. Future studies investigating the contribution of promiscuity to enzyme and metabolic evolution will need to integrate deeper probes into the influence of evolution on protein biophysics, enzymology, and metabolism with more complex and realistic evolutionary models. ENZYMES: lactate dehydrogenase (EC 1.1.1.27), malate dehydrogenase (EC 1.1.1.37), OSBS (EC 4.2.1.113), HisA (EC 5.3.1.16), TrpF, PriA (EC 5.3.1.24), R-mandelonitrile lyase (EC 4.1.2.10), Maleylacetate reductase (EC 1.3.1.32).

Gene duplications in the E. coli genome: common themes among pathotypes

Background

Gene duplication underlies a significant proportion of gene functional diversity and genome complexity in both eukaryotes and prokaryotes. Although several reports in the literature described the duplication of specific genes in E. coli, a detailed analysis of the extent of gene duplications in this microorganism is needed.

Results

The genomes of the E. coli enteroaggregative strain 042 and other pathogenic strains contain duplications of the gene that codes for the global regulator Hha. To determine whether the presence of additional copies of the hha gene correlates with the presence of other genes, we performed a comparative genomic analysis between E. coli strains with and without hha duplications. The results showed that strains harboring additional copies of the hha gene also encode the yeeR irmA (aec69) gene cluster, which, in turn, is also duplicated in strain 042 and several other strains. The identification of these duplications prompted us to obtain a global map of gene duplications, first in strain 042 and later in other E. coli genomes.

Duplications in the genomes of the enteroaggregative strain 042, the uropathogenic strain CFT073 and the enterohemorrhagic strain O145:H28 have been identified by a BLASTp protein similarity search. This algorithm was also used to evaluate the distribution of the identified duplicates among the genomes of a set of 28 representative E. coli strains. Despite the high genomic diversity of E. coli strains, we identified several duplicates in the genomes of almost all studied pathogenic strains. Most duplicated genes have no known function. Transcriptomic analysis also showed that most of these duplications are regulated by the H-NS/Hha proteins.

Conclusions

Several duplicated genes are widely distributed among pathogenic E. coli strains. In addition, some duplicated genes are present only in specific pathotypes, and others are strain specific. This gene duplication analysis shows novel relationships between E. coli pathotypes and suggests that newly identified genes that are duplicated in a high percentage of pathogenic E. coli isolates may play a role in virulence. Our study also shows a relationship between the duplication of genes encoding regulators and genes encoding their targets.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5683-4) contains supplementary material, which is available to authorized users.

Comparative proteome analysis reveals pathogen specific outer membrane proteins of Leptospira

Proteomes of pathogenic Leptospira interrogans and L. borgpetersenii and the saprophytic L. biflexa were filtered through computational tools to identify Outer Membrane Proteins (OMPs) that satisfy the required biophysical parameters for their presence on the outer membrane. A total of 133, 130, and 144 OMPs were identified in L. interrogans, L. borgpetersenii, and L. biflexa, respectively, which forms approximately 4% of proteomes. A holistic analysis of transporting and pathogenic characteristics of OMPs together with Clusters of Orthologous Groups (COGs) among the OMPs and their distribution across 3 species was made and put forward a set of 21 candidate OMPs specific to pathogenic leptospires. It is also found that proteins homologous to the candidate OMPs were also present in other pathogenic species of leptospires. Six OMPs from L. interrogans and 2 from L. borgpetersenii observed to have similar COGs while those were not found in any intermediate or saprophytic forms. These OMPs appears to have role in infection and pathogenesis and useful for anti-leptospiral strategies.

Selection Effects on the Positioning of Genes and Gene Structures from the Interplay of Replication and Transcription in Bacterial Genomes

Bacterial chromosomes are partly shaped by the functional requirements for efficient replication, which lead to strand bias as commonly characterized by the excess of guanines over cytosines in the leading strand. Gene structures are also highly organized within bacterial genomes as a result of such functional constraints, displaying characteristic positioning and structuring along the genome. Here we analyze the gene structures in completely sequenced bacterial chromosomes to observe the positional constraints on gene orientation, length, and codon usage with regard to the positions of replication origin and terminus. Selection on these gene features is different in regions surrounding the terminus of replication from the rest of the genome, but the selection could be either positive or negative depending on the species, and these positional effects are partly attributed to the A-T enrichment near the terminus. Characteristic gene structuring relative to the position of replication origin and terminus is commonly observed among most bacterial species with circular chromosomes, and therefore we argue that the highly organized gene positioning as well as the strand bias should be considered for genomics studies of bacteria.

Phylogenomic networks reveal limited phylogenetic range of lateral gene transfer by transduction

Bacteriophages are recognized DNA vectors and transduction is considered as a common mechanism of lateral gene transfer (LGT) during microbial evolution. Anecdotal events of phage-mediated gene transfer were studied extensively, however, a coherent evolutionary viewpoint of LGT by transduction, its extent and characteristics, is still lacking. Here we report a large-scale evolutionary reconstruction of transduction events in 3982 genomes. We inferred 17 158 recent transduction events linking donors, phages and recipients into a phylogenomic transduction network view. We find that LGT by transduction is mostly restricted to closely related donors and recipients. Furthermore, a substantial number of the transduction events (9%) are best described as gene duplications that are mediated by mobile DNA vectors. We propose to distinguish this type of paralogy by the term autology. A comparison of donor and recipient genomes revealed that genome similarity is a superior predictor of species connectivity in the network in comparison to common habitat. This indicates that genetic similarity, rather than ecological opportunity, is a driver of successful transduction during microbial evolution. A striking difference in the connectivity pattern of donors and recipients shows that while lysogenic interactions are highly species-specific, the host range for lytic phage infections can be much wider, serving to connect dense clusters of closely related species. Our results thus demonstrate that DNA transfer via transduction occurs within the context of phage–host specificity, but that this tight constraint can be breached, on rare occasions, to produce long-range LGTs of profound evolutionary consequences.

Gene-Family Extension Measures and Correlations

The existence of multiple copies of genes is a well-known phenomenon. A gene family is a set of sufficiently similar genes, formed by gene duplication. In earlier works conducted on a limited number of completely sequenced and annotated genomes it was found that size of gene family and size of genome are positively correlated. Additionally, it was found that several atypical microbes deviated from the observed general trend. In this study, we reexamined these associations on a larger dataset consisting of 1484 prokaryotic genomes and using several ranking approaches. We applied ranking methods in such a way that genomes with lower numbers of gene copies would have lower rank. Until now only simple ranking methods were used; we applied the Kemeny optimal aggregation approach as well. Regression and correlation analysis were utilized in order to accurately quantify and characterize the relationships between measures of paralog indices and genome size. In addition, boxplot analysis was employed as a method for outlier detection. We found that, in general, all paralog indexes positively correlate with an increase of genome size. As expected, different groups of atypical prokaryotic genomes were found for different types of paralog quantities. Mycoplasmataceae and Halobacteria appeared to be among the most interesting candidates for further research of evolution through gene duplication.

Whole genome sequencing revealed host adaptation-focused genomic plasticity of pathogenic Leptospira

Leptospirosis, caused by pathogenic Leptospira spp., has recently been recognized as an emerging infectious disease worldwide. Despite its severity and global importance, knowledge about the molecular pathogenesis and virulence evolution of Leptospira spp. remains limited. Here we sequenced and analyzed 102 isolates representing global sources. A high genomic variability were observed among different Leptospira species, which was attributed to massive gene gain and loss events allowing for adaptation to specific niche conditions and changing host environments. Horizontal gene transfer and gene duplication allowed the stepwise acquisition of virulence factors in pathogenic Leptospira evolved from a recent common ancestor. More importantly, the abundant expansion of specific virulence-related protein families, such as metalloproteases-associated paralogs, were exclusively identified in pathogenic species, reflecting the importance of these protein families in the pathogenesis of leptospirosis. Our observations also indicated that positive selection played a crucial role on this bacteria adaptation to hosts. These novel findings may lead to greater understanding of the global diversity and virulence evolution of Leptospira spp.

The Dynamics of Lateral Gene Transfer in Genus Leishmania - A Route for Adaptation and Species Diversification

Background

The genome of Leishmania major harbours a comparably high proportion of genes of prokaryote origin, acquired by lateral gene transfer (LGT). Some of these are present in closely related trypanosomatids, while some are detected in Leishmania only. We have evaluated the impact and destiny of LGT in genus Leishmania.

Methodology/Principal Findings

To study the dynamics and fate of LGTs we have performed phylogenetic, as well as nucleotide and amino acid composition analyses within orthologous groups of LGTs detected in Leishmania. A set of universal trypanosomatid LGTs was added as a reference group. Both groups of LGTs have, to some extent, ameliorated to resemble the recipient genomes. However, while virtually all of the universal trypanosomatid LGTs are distributed and conserved in the entire genus Leishmania, the LGTs uniquely present in genus Leishmania are more prone to gene loss and display faster rates of evolution. Furthermore, a PCR based approach has been employed to ascertain the presence of a set of twenty LGTs uniquely present in genus Leishmania, and three universal trypanosomatid LGTs, in ten additional strains of Leishmania. Evolutionary rates and predicted expression levels of these LGTs have also been estimated. Ten of the twenty LGTs are distributed and conserved in all species investigated, while the remainder have been subjected to modifications, or undergone pseudogenization, degradation or loss in one or more species.

Conclusions/Significance

LGTs unique to the genus Leishmania have been acquired after the divergence of Leishmania from the other trypanosomatids, and are evolving faster than their recipient genomes. This implies that LGT in genus Leishmania is a continuous and dynamic process contributing to species differentiation and speciation. This study also highlights the importance of carefully evaluating these dynamic genes, e.g. as LGTs have been suggested as potential drug targets.

Leishmania parasites cause leishmaniasis, a neglected tropical disease, estimated to threaten 350 million people in 88 countries worldwide according to the WHO. The genome of Leishmania major harbours a number of genes, which have been proposed as acquired by lateral gene transfer (LGT) from a broad variety of prokaryote donors. Such genes may prove beneficial for the parasites, e.g. by promoting survival of the parasite in new environments. We have studied orthologs to LGTs previously detected uniquely in L. major as well as LGTs shared also by other trypanosomatids. The universal trypanosomatid LGTs are more conserved within genus Leishmania, as compared to LGTs that are exclusively present in genus Leishmania. One possible explanation to this observation is that these have resided in their host genomes for a longer time period. This explanation strengthens the hypothesis that the LGTs unique to genus Leishmaina were acquired after the divergence from the trypanosomes, rather than before the divergence, and then subsequently lost from the trypanosome lineage. An in-depth analysis of a subset of the LGTs, which are present only in genus Leishmania showed that LGT within genus Leishmania is a dynamic process. LGTs, providing beneficial capabilities to the parasite, are demonstrated to become conserved throughout generic diversification, hence contributing to species differentiation, while LGTs of limited use are degraded and lost.

The evolutionary dynamics of tRNA-gene copy number and codon-use in E. coli

Background

The introduction of foreign DNA by Lateral Gene Transfer (LGT) can quickly and drastically alter genome composition. Problems can arise if the genes introduced by LGT use codons that are not suited to the host’s translational machinery. Here we investigate compensatory adaptation of E. coli in response to the introduction of large volumes of codons that are rarely used by the host genome.

Results

We analyze genome sequences from the E. coli/Shigella complex, and find that certain tRNA genes are present in multiple copies in two pathogenic Shigella and O157:H7 subgroups of E. coli. Furthermore, we show that the codons that correspond to these multi-copy number tRNA genes are enriched in the high copy number Selfish Genetic Elements (SGE’s) in Shigella and laterally introduced genes in O157:H7. We analyze the duplicate copies and find evidence for the selective retention of tRNA genes introduced by LGT in response to the changed codon content of the genome.

Conclusion

These data support a model where the relatively rapid influx of LGT genes and SGE’s introduces a large number of genes maladapted to the host’s translational machinery. Under these conditions, it becomes advantageous for the host to retain tRNA genes that are required for the incorporation of amino acids at these codons. Subsequently, the increased number of copies of these specific tRNA genes adjusts the cellular tRNA pool to the demands set by global shifts in codon usage.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0441-y) contains supplementary material, which is available to authorized users.

Tackling Drug Resistant Infection Outbreaks of Global Pandemic Escherichia coli ST131 Using Evolutionary and Epidemiological Genomics

High-throughput molecular screening is required to investigate the origin and diffusion of antimicrobial resistance in pathogen outbreaks. The most frequent cause of human infection is Escherichia coli, which is dominated by sequence type 131 (ST131)-a set of rapidly radiating pandemic clones. The highly infectious clades of ST131 originated firstly by a mutation enhancing conjugation and adhesion. Secondly, single-nucleotide polymorphisms occurred enabling fluoroquinolone-resistance, which is near-fixed in all ST131. Thirdly, broader resistance through beta-lactamases has been gained and lost frequently, symptomatic of conflicting environmental selective effects. This flexible approach to gene exchange is worrying and supports the proposition that ST131 will develop an even wider range of plasmid and chromosomal elements promoting antimicrobial resistance. To stop ST131, deep genome sequencing is required to understand the origin, evolution and spread of antimicrobial resistance genes. Phylogenetic methods that decipher past events can predict future patterns of virulence and transmission based on genetic signatures of adaptation and gene exchange. Both the effect of partial antimicrobial exposure and cell dormancy caused by variation in gene expression may accelerate the development of resistance. High-throughput sequencing can decode measurable evolution of cell populations within patients associated with systems-wide changes in gene expression during treatments. A multi-faceted approach can enhance assessment of antimicrobial resistance in E. coli ST131 by examining transmission dynamics between hosts to achieve a goal of pre-empting resistance before it emerges by optimising antimicrobial treatment protocols.

A recently transferred cluster of bacterial genes in Trichomonas vaginalis--lateral gene transfer and the fate of acquired genes

Background

Lateral Gene Transfer (LGT) has recently gained recognition as an important contributor to some eukaryote proteomes, but the mechanisms of acquisition and fixation in eukaryotic genomes are still uncertain. A previously defined norm for LGTs in microbial eukaryotes states that the majority are genes involved in metabolism, the LGTs are typically localized one by one, surrounded by vertically inherited genes on the chromosome, and phylogenetics shows that a broad collection of bacterial lineages have contributed to the transferome.

Results

A unique 34 kbp long fragment with 27 clustered genes (TvLF) of prokaryote origin was identified in the sequenced genome of the protozoan parasite Trichomonas vaginalis. Using a PCR based approach we confirmed the presence of the orthologous fragment in four additional T. vaginalis strains. Detailed sequence analyses unambiguously suggest that TvLF is the result of one single, recent LGT event. The proposed donor is a close relative to the firmicute bacterium Peptoniphilus harei. High nucleotide sequence similarity between T. vaginalis strains, as well as to P. harei, and the absence of homologs in other Trichomonas species, suggests that the transfer event took place after the radiation of the genus Trichomonas. Some genes have undergone pseudogenization and degradation, indicating that they may not be retained in the future. Functional annotations reveal that genes involved in informational processes are particularly prone to degradation.

Conclusions

We conclude that, although the majority of eukaryote LGTs are single gene occurrences, they may be acquired in clusters of several genes that are subsequently cleansed of evolutionarily less advantageous genes.

Safety in numbers: multiple occurrences of highly similar homologs among Azotobacter vinelandii carbohydrate metabolism proteins probably confer adaptive benefits

BACKGROUND

Gene duplication and horizontal gene transfer are common processes in bacterial and archaeal genomes, and are generally assumed to result in either diversification or loss of the redundant gene copies. However, a recent analysis of the genome of the soil bacterium Azotobacter vinelandii DJ revealed an abundance of highly similar homologs among carbohydrate metabolism genes. In many cases these multiple genes did not appear to be the result of recent duplications, or to function only as a means of stimulating expression by increasing gene dosage, as the homologs were located in varying functional genetic contexts. Based on these initial findings we here report in-depth bioinformatic analyses focusing specifically on highly similar intra-genome homologs, or synologs, among carbohydrate metabolism genes, as well as an analysis of the general occurrence of very similar synologs in prokaryotes.

RESULTS

Approximately 900 bacterial and archaeal genomes were analysed for the occurrence of synologs, both in general and among carbohydrate metabolism genes specifically. This showed that large numbers of highly similar synologs among carbohydrate metabolism genes are very rare in bacterial and archaeal genomes, and that the A. vinelandii DJ genome contains an unusually large amount of such synologs. The majority of these synologs were found to be non-tandemly organized and localized in varying but metabolically relevant genomic contexts. The same observation was made for other genomes harbouring high levels of such synologs. It was also shown that highly similar synologs generally constitute a very small fraction of the protein-coding genes in prokaryotic genomes. The overall synolog fraction of the A. vinelandii DJ genome was well above the data set average, but not nearly as remarkable as the levels observed when only carbohydrate metabolism synologs were considered.

CONCLUSIONS

Large numbers of highly similar synologs are rare in bacterial and archaeal genomes, both in general and among carbohydrate metabolism genes. However, A. vinelandii and several other soil bacteria harbour large numbers of highly similar carbohydrate metabolism synologs which seem not to result from recent duplication or transfer events. These genes may confer adaptive benefits with respect to certain lifestyles and environmental factors, most likely due to increased regulatory flexibility and/or increased gene dosage.

Copy number change: evolving views on gene amplification

The rapid pace of genomic sequence analysis is increasing the awareness of intrinsically dynamic genetic landscapes. Gene duplication and amplification (GDA) contribute to adaptation and evolution by allowing DNA regions to expand and contract in an accordion-like fashion. This process affects diverse aspects of bacterial infection, including antibiotic resistance and host-pathogen interactions. In this review, microbial GDA is discussed, primarily using recent bacterial examples that demonstrate medical and evolutionary consequences. Interplay between GDA and horizontal gene transfer further impact evolutionary trajectories. Complementing the discovery of gene duplication in clinical and environmental settings, experimental evolution provides a powerful method to document genetic change over time. New methods for GDA detection highlight both its importance and its potential application for genetic engineering, synthetic biology and biotechnology.

Real-time evolution of new genes by innovation, amplification, and divergence

Gene duplications allow evolution of genes with new functions. Here, we describe the innovation-amplification-divergence (IAD) model in which the new function appears before duplication and functionally distinct new genes evolve under continuous selection. One example fitting this model is a preexisting parental gene in Salmonella enterica that has low levels of two distinct activities. This gene is amplified to a high copy number, and the amplified gene copies accumulate mutations that provide enzymatic specialization of different copies and faster growth. Selection maintains the initial amplification and beneficial mutant alleles but is relaxed for other less improved gene copies, allowing their loss. This rapid process, completed in fewer than 3000 generations, shows the efficacy of the IAD model and allows the study of gene evolution in real time.

The ecology of bacterial genes and the survival of the new

Much of the observed variation among closely related bacterial genomes is attributable to gains and losses of genes that are acquired horizontally as well as to gene duplications and larger amplifications. The genomic flexibility that results from these mechanisms certainly contributes to the ability of bacteria to survive and adapt in varying environmental challenges. However, the duplicability and transferability of individual genes imply that natural selection should operate, not only at the organismal level, but also at the level of the gene. Genes can be considered semiautonomous entities that possess specific functional niches and evolutionary dynamics. The evolution of bacterial genes should respond both to selective pressures that favor competition, mostly among orthologs or paralogs that may occupy the same functional niches, and cooperation, with the majority of other genes coexisting in a given genome. The relative importance of either type of selection is likely to vary among different types of genes, based on the functional niches they cover and on the tightness of their association with specific organismal lineages. The frequent availability of new functional niches caused by environmental changes and biotic evolution should enable the constant diversification of gene families and the survival of new lineages of genes.

Evolution of cariogenic character in Streptococcus mutans: horizontal transmission of glycosyl hydrolase family 70 genes

Acquisition of the ability to produce polysaccharides from sucrose, i.e. the gtf gene encoding glucosyltransferase (GTF), is the key evolutionary event enabling dental biofilm formation by streptococci. To clarify the ancestry of streptococcal GTFs, time of its occurrence, and order of specific events, we investigated the distribution of GTFs among bacteria by phylogenetic analysis of the glycosyl hydrolase family 70 enzymes. We found that streptococcal GTFs were derived from other lactic acid bacteria such as Lactobacillus and Leuconostoc, and propose the following evolutionary model: horizontal gene transfer via transposons occurred when streptococci encountered lactic acid bacteria contained in fermented food. Intra-genomic gene duplication occurred by a secondary selection pressure such as consumption of refined sugar. Our findings concerning this evolution in Streptococcus mutans provide an important background for studies of the relationship between the historical spread of dental caries and anthropological factors.

Ancient gene transfer from algae to animals: mechanisms and evolutionary significance

Background

Horizontal gene transfer (HGT) is traditionally considered to be rare in multicellular eukaryotes such as animals. Recently, many genes of miscellaneous algal origins were discovered in choanoflagellates. Considering that choanoflagellates are the existing closest relatives of animals, we speculated that ancient HGT might have occurred in the unicellular ancestor of animals and affected the long-term evolution of animals.

Results

Through genome screening, phylogenetic and domain analyses, we identified 14 gene families, including 92 genes, in the tunicate Ciona intestinalis that are likely derived from miscellaneous photosynthetic eukaryotes. Almost all of these gene families are distributed in diverse animals, suggesting that they were mostly acquired by the common ancestor of animals. Their miscellaneous origins also suggest that these genes are not derived from a particular algal endosymbiont. In addition, most genes identified in our analyses are functionally related to molecule transport, cellular regulation and methylation signaling, suggesting that the acquisition of these genes might have facilitated the intercellular communication in the ancestral animal.

Conclusions

Our findings provide additional evidence that algal genes in aplastidic eukaryotes are not exclusively derived from historical plastids and thus important for interpreting the evolution of eukaryotic photosynthesis. Most importantly, our data represent the first evidence that more anciently acquired genes might exist in animals and that ancient HGT events have played an important role in animal evolution.

The ancestral SgrS RNA discriminates horizontally acquired Salmonella mRNAs through a single G-U wobble pair

SgrS RNA is a model for the large class of Hfq-associated small RNAs that act to posttranscriptionally regulate bacterial mRNAs. The function of SgrS is well-characterized in nonpathogenic Escherichia coli, where it was originally shown to counteract glucose-phosphate stress by acting as a repressor of the ptsG mRNA, which encodes the major glucose transporter. We have discovered additional SgrS targets in Salmonella Typhimurium, a pathogen related to E. coli that recently acquired one-quarter of all genes by horizontal gene transfer. We show that the conserved short seed region of SgrS that recognizes ptsG was recruited to target the Salmonella-specific sopD mRNA of a secreted virulence protein. The SgrS-sopD interaction is exceptionally selective; we find that sopD2 mRNA, whose gene arose from sopD duplication during Salmonella evolution, is deaf to SgrS because of a nonproductive G-U pair in the potential SgrS-sopD2 RNA duplex vs. G-C in SgrS-sopD. In other words, SgrS discriminates the two virulence factor mRNAs at the level of a single hydrogen bond. Our study suggests that bacterial pathogens use their large suites of conserved Hfq-associated regulators to integrate horizontally acquired genes into existing posttranscriptional networks, just as conserved transcription factors are recruited to tame foreign genes at the DNA level. The results graphically illustrate the importance of the seed regions of bacterial small RNAs to select new targets with high fidelity and suggest that target predictions must consider all or none decisions by individual seed nucleotides.

Pathogen-origin horizontally transferred genes contribute to the evolution of Lepidopteran insects

Background

Horizontal gene transfer (HGT), a source of genetic variation, is generally considered to facilitate hosts' adaptability to environments. However, convincing evidence supporting the significant contribution of the transferred genes to the evolution of metazoan recipients is rare.

Results

In this study, based on sequence data accumulated to date, we used a unified method consisting of similarity search and phylogenetic analysis to detect horizontally transferred genes (HTGs) between prokaryotes and five insect species including Drosophila melanogaster, Anopheles gambiae, Bombyx mori, Tribolium castaneum and Apis mellifera. Unexpectedly, the candidate HTGs were not detected in D. melanogaster, An. gambiae and T. castaneum, and 79 genes in Ap. mellifera sieved by the same method were considered as contamination based on other information. Consequently, 14 types of 22 HTGs were detected only in the silkworm. Additionally, 13 types of the detected silkworm HTGs share homologous sequences in species of other Lepidopteran superfamilies, suggesting that the majority of these HTGs were derived from ancient transfer events before the radiation of Ditrysia clade. On the basis of phylogenetic topologies and BLAST search results, donor bacteria of these genes were inferred, respectively. At least half of the predicted donor organisms may be entomopathogenic bacteria. The predicted biochemical functions of these genes include four categories: glycosyl hydrolase family, oxidoreductase family, amino acid metabolism, and others.

Conclusions

The products of HTGs detected in this study may take part in comprehensive physiological metabolism. These genes potentially contributed to functional innovation and adaptability of Lepidopteran hosts in their ancient lineages associated with the diversification of angiosperms. Importantly, our results imply that pathogens may be advantageous to the subsistence and prosperity of hosts through effective HGT events at a large evolutionary scale.

Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches

Horizontal genetic transfer (HGT) has played an important role in bacterial evolution at least since the origins of the bacterial divisions, and HGT still facilitates the origins of bacterial diversity, including diversity based on antibiotic resistance. Adaptive HGT is aided by unique features of genetic exchange in bacteria such as the promiscuity of genetic exchange and the shortness of segments transferred. Genetic exchange rates are limited by the genetic and ecological similarity of organisms. Adaptive transfer of genes is limited to those that can be transferred as a functional unit, provide a niche-transcending adaptation, and are compatible with the architecture and physiology of other organisms. Horizontally transferred adaptations may bring about fitness costs, and natural selection may ameliorate these costs. The origins of ecological diversity can be analyzed by comparing the genomes of recently divergent, ecologically distinct populations, which can be discovered as sequence clusters. Such genome comparisons demonstrate the importance of HGT in ecological diversification. Newly divergent populations cannot be discovered as sequence clusters when their ecological differences are coded by plasmids, as is often the case for antibiotic resistance; the discovery of such populations requires a screen for plasmid-coded functions. This paper reviews the features of bacterial genetics that allow HGT, the similarities between organisms that foster HGT between them, the limits to the kinds of adaptations that can be transferred, and amelioration of fitness costs associated with HGT; the paper also reviews approaches to discover the origins of new, ecologically distinct bacterial populations and the role that HGT plays in their founding.

Dynamics of gene duplication in the genomes of chlorophyll d-producing cyanobacteria: implications for the ecological niche

Gene duplication may be an important mechanism for the evolution of new functions and for the adaptive modulation of gene expression via dosage effects. Here, we analyzed the fate of gene duplicates for two strains of a novel group of cyanobacteria (genus Acaryochloris) that produces the far-red light absorbing chlorophyll d as its main photosynthetic pigment. The genomes of both strains contain an unusually high number of gene duplicates for bacteria. As has been observed for eukaryotic genomes, we find that the demography of gene duplicates can be well modeled by a birth-death process. Most duplicated Acaryochloris genes are of comparatively recent origin, are strain-specific, and tend to be located on different genetic elements. Analyses of selection on duplicates of different divergence classes suggest that a minority of paralogs exhibit near neutral evolutionary dynamics immediately following duplication but that most duplicate pairs (including those which have been retained for long periods) are under strong purifying selection against amino acid change. The likelihood of duplicate retention varied among gene functional classes, and the pronounced differences between strains in the pool of retained recent duplicates likely reflects differences in the nutrient status and other characteristics of their respective environments. We conclude that most duplicates are quickly purged from Acaryochloris genomes and that those which are retained likely make important contributions to organism ecology by conferring fitness benefits via gene dosage effects. The mechanism of enhanced duplication may involve homologous recombination between genetic elements mediated by paralogous copies of recA.

Lateral transfer of genes and gene fragments in Staphylococcus extends beyond mobile elements

The widespread presence of antibiotic resistance and virulence among Staphylococcus isolates has been attributed in part to lateral genetic transfer (LGT), but little is known about the broader extent of LGT within this genus. Here we report the first systematic study of the modularity of genetic transfer among 13 Staphylococcus genomes covering four distinct named species. Using a topology-based phylogenetic approach, we found, among 1,354 sets of homologous genes examined, strong evidence of LGT in 368 (27.1%) gene sets, and weaker evidence in another 259 (19.1%). Within-gene and whole-gene transfer contribute almost equally to the topological discordance of these gene sets against a reference phylogeny. Comparing genetic transfer in single-copy and in multicopy gene sets, we observed a higher frequency of LGT in the latter, and a substantial functional bias in cases of whole-gene transfer (little such bias was observed in cases of fragmentary genetic transfer). We found evidence that lateral transfer, particularly of entire genes, impacts not only functions related to antibiotic, drug, and heavy-metal resistance, as well as membrane transport, but also core informational and metabolic functions not associated with mobile elements. Although patterns of sequence similarity support the cohesion of recognized species, LGT within S. aureus appears frequently to disrupt clonal complexes. Our results demonstrate that LGT and gene duplication play important parts in functional innovation in staphylococcal genomes.

Horizontal transfer, not duplication, drives the expansion of protein families in prokaryotes

Gene duplication followed by neo- or sub-functionalization deeply impacts the evolution of protein families and is regarded as the main source of adaptive functional novelty in eukaryotes. While there is ample evidence of adaptive gene duplication in prokaryotes, it is not clear whether duplication outweighs the contribution of horizontal gene transfer in the expansion of protein families. We analyzed closely related prokaryote strains or species with small genomes (Helicobacter, Neisseria, Streptococcus, Sulfolobus), average-sized genomes (Bacillus, Enterobacteriaceae), and large genomes (Pseudomonas, Bradyrhizobiaceae) to untangle the effects of duplication and horizontal transfer. After removing the effects of transposable elements and phages, we show that the vast majority of expansions of protein families are due to transfer, even among large genomes. Transferred genes--xenologs--persist longer in prokaryotic lineages possibly due to a higher/longer adaptive role. On the other hand, duplicated genes--paralogs--are expressed more, and, when persistent, they evolve slower. This suggests that gene transfer and gene duplication have very different roles in shaping the evolution of biological systems: transfer allows the acquisition of new functions and duplication leads to higher gene dosage. Accordingly, we show that paralogs share most protein-protein interactions and genetic regulators, whereas xenologs share very few of them. Prokaryotes invented most of life's biochemical diversity. Therefore, the study of the evolution of biology systems should explicitly account for the predominant role of horizontal gene transfer in the diversification of protein families.

Inferring bacterial genome flux while considering truncated genes

Bacterial gene content variation during the course of evolution has been widely acknowledged and its pattern has been actively modeled in recent years. Gene truncation or gene pseudogenization also plays an important role in shaping bacterial genome content. Truncated genes could also arise from small-scale lateral gene transfer events. Unfortunately, the information of truncated genes has not been considered in any existing mathematical models on gene content variation. In this study, we developed a model to incorporate truncated genes. Maximum-likelihood estimates (MLEs) of the new model reveal fast rates of gene insertions/deletions on recent branches, suggesting a fast turnover of many recently transferred genes. The estimates also suggest that many truncated genes are in the process of being eliminated from the genome. Furthermore, we demonstrate that the ignorance of truncated genes in the estimation does not lead to a systematic bias but rather has a more complicated effect. Analysis using the new model not only provides more accurate estimates on gene gains/losses (or insertions/deletions), but also reduces any concern of a systematic bias from applying simplified models to bacterial genome evolution. Although not a primary purpose, the model incorporating truncated genes could be potentially used for phylogeny reconstruction using gene family content.

Emergence and evolution of modern molecular functions inferred from phylogenomic analysis of ontological data

The biological processes that characterize the phenotypes of a living system are embodied in the function of molecules and hold the key to evolutionary history, delimiting natural selection and change. These processes and functions provide direct insight into the emergence, development, and organization of cellular life. However, detailed molecular functions make up a network-like hierarchy of relationships that tells little of evolutionary links between structure and function in biology. For example, Gene Ontology terms represent widely-used vocabularies of processes and functions with evolutionary relationships that are implicit but not defined. Here, we uncover patterns of global evolutionary history in ontological terms associated with the sequence of 38 genomes. These patterns unfold the metabolic origins of modern molecular functions and major biological transitions in evolution toward complex life. Phylogenies reveal the primordial appearance of hydrolases and transferases, with ATPase, GTPase, and helicase activities being the most ancient. This indicates that ancient catalysts were crucial for binding and transport, the emergence of nucleic acids and protein biopolymers, and the communication of primordial cells with the environment. Finally, the history of biological processes showed that cellular biopolymer metabolic processes preceded biopolymer biosynthesis and essential processes related to macromolecular formation, directly challenging the existence of an RNA world. Phylogenomic systematization of biological function takes the structure and function paradigm to a completely new level of abstraction, demonstrating a "metabolic first" origin of life. The approach uncovers patterns in the morphing of function that are unprecedented and necessary for systematic views in biology.

Compensatory gene amplification restores fitness after inter-species gene replacements

Genes introduced by gene replacements and other types of horizontal gene transfer (HGT) represent a significant presence in many archaeal and eubacterial genomes. Most alien genes are likely to be neutral or deleterious upon arrival and their long-term persistence may require a mechanism that improves their selective contribution. To examine the fate of inter-species gene replacements, we exchanged three native S. typhimurium genes encoding ribosomal proteins with orthologues from various other microbes. The results show that replacement of each of these three genes reduces fitness to such an extent that it would provide an effective barrier against inter-species gene replacements in eubacterial populations. However, these fitness defects could be partially ameliorated by gene amplification that augmented the dosage of the heterologous proteins. This suggests that suboptimal expression is a common fitness constraint for inter-species gene replacements, with fitness costs conferred by either a lower expression level of the alien protein compared with the native protein or a requirement for an increased amount of the alien protein to maintain proper function. Our findings can explain the observation that duplicated genes are over-represented among horizontally transferred genes, and suggest a potential coupling between compensatory gene amplification after HGT and the evolution of new genes.

Versatility of choline metabolism and choline-binding proteins in Streptococcus pneumoniae and commensal streptococci

The pneumococcal choline-containing teichoic acids are targeted by cholinebinding proteins (CBPs), major surface components implicated in the interaction with host cells and bacterial cell physiology. CBPs also occur in closely related commensal species, Streptococcus oralis and Streptococcus mitis, and many strains of these species contain choline in their cell wall. Physiologically relevant CBPs including cell wall lytic enzymes are highly conserved between Streptococcus pneumoniae and S. mitis. In contrast, the virulence-associated CBPs, CbpA, PspA and PcpA, are S. pneumoniae specific and are thus relevant for the characteristic properties of this species.

On the nature of fur evolution: a phylogenetic approach in Actinobacteria

Background

An understanding of the evolution of global transcription regulators is essential for comprehending the complex networks of cellular metabolism that have developed among related organisms. The fur gene encodes one of those regulators – the ferric uptake regulator Fur – widely distributed among bacteria and known to regulate different genes committed to varied metabolic pathways. On the other hand, members of the Actinobacteria comprise an ecologically diverse group of bacteria able to inhabit various natural environments, and for which relatively little is currently understood concerning transcriptional regulation.

Results

BLAST analyses revealed the presence of more than one fur homologue in most members of the Actinobacteria whose genomes have been fully sequenced. We propose a model to explain the evolutionary history of fur within this well-known bacterial phylum: the postulated scenario includes one duplication event from a primitive regulator, which probably had a broad range of co-factors and DNA-binding sites. This duplication predated the appearance of the last common ancestor of the Actinobacteria, while six other duplications occurred later within specific groups of organisms, particularly in two genera: Frankia and Streptomyces. The resulting paralogues maintained main biochemical properties, but became specialised for regulating specific functions, coordinating different metal ions and binding to unique DNA sequences. The presence of syntenic regions surrounding the different fur orthologues supports the proposed model, as do the evolutionary distances and topology of phylogenetic trees built using both Neighbor-Joining and Maximum-Likelihood methods.

Conclusion

The proposed fur evolutionary model, which includes one general duplication and two in-genus duplications followed by divergence and specialization, explains the presence and diversity of fur genes within the Actinobacteria. Although a few rare horizontal gene transfer events have been reported, the model is consistent with the view of gene duplication as a main force of microbial genomes evolution. The parallel study of Fur phylogeny across diverse organisms offers a solid base to guide functional studies and allows the comparison between response mechanisms in relation with the surrounding environment. The survey of regulators among related genomes provides a relevant tool for understanding the evolution of one of the first lines of cellular adaptability, control of DNA transcription.

Evolution of new gene functions: simulation and analysis of the amplification model

Creation of new genes and functions is a central feature of evolution. Duplication of existing genes has long been assumed to be the source of new genes, but the precise mechanism has remained unclear. One suggestion is that new genes are created via temporary amplifications, which simultaneously increase both the selective advantage of weak, pre-existing secondary functions and the target for optimizing mutations. This paper examines the amplification model by formalizing it into a mathematical framework. This framework is used to perform stochastic (Monte Carlo) simulations. In addition, experimental data from Salmonella typhimurium LT2 are used to support the modelling, by providing estimates for parameter values. The results show that amplification of tandem repeats is likely to contribute to creation of new genes in nature.

The Lactococcus lactis FabF fatty acid synthetic enzyme can functionally replace both the FabB and FabF proteins of Escherichia coli and the FabH protein of Lactococcus lactis

The genome of Lactococcus lactis encodes a single long chain 3-ketoacyl-acyl carrier protein synthase. This is in contrast to its close relative, Enterococcus faecalis, and to Escherichia coli, both of which have two such enzymes. In E. faecalis and E. coli, one of the two long chain synthases (FabO and FabB, respectively) has a role in unsaturated fatty acid synthesis that cannot be satisfied by FabF, the other long chain synthase. Since L. lactis has only a single long chain 3-ketoacyl-acyl carrier protein synthase (annotated as FabF), it seemed likely that this enzyme must function both in unsaturated fatty acid synthesis and in elongation of short chain acyl carrier protein substrates to the C18 fatty acids found in the cellular phospholipids. We report that this is the case. Expression of L. lactis FabF can functionally replace both FabB and FabF in E. coli, although it does not restore thermal regulation of phospholipid fatty acid composition to E. coli fabF mutant strains. The lack of thermal regulation was predictable because wild-type L. lactis was found not to show any significant change in fatty acid composition with growth temperature. We also report that overproduction of L. lactis FabF allows growth of an L. lactis mutant strain that lacks the FabH short chain 3-ketoacyl-acyl carrier protein synthase. The strain tested was a derivative (called the fabH bypass strain) of the original fabH deletion strain that had acquired the ability to grow when supplemented with octanoate. Upon introduction of a FabF overexpression plasmid into this strain, growth proceeded normally in the absence of fatty acid supplementation. Moreover, this strain had a normal rate of fatty acid synthesis and a normal fatty acid composition. Both the fabH bypass strain that overproduced FabF and the wild type strain incorporated much less exogenous octanoate into long chain phospholipid fatty acids than did the fabH bypass strain. Incorporation of octanoate and decanoate labeled with deuterium showed that these acids were incorporated intact as the distal methyl and methylene groups of the long chain fatty acids.

Uncovering rate variation of lateral gene transfer during bacterial genome evolution

BACKGROUND

Large scale genome arrangement, such as whole gene insertion/deletion, plays an important role in bacterial genome evolution. Various methods have been employed to study the dynamic process of gene insertions and deletions, such as parsimony methods and maximum likelihood methods. Previous maximum likelihood studies have assumed that the rate of gene insertions/deletions is constant over different genes. This assumption is unrealistic. For instance, it has been shown that informational genes are less likely to be laterally transferred than non-informational genes. However, how much of the variation in gene transfer rates is due to the difference between informational genes and non-informational genes is unclear. In this study, a Gamma-distribution was incorporated in the likelihood estimation by considering rate variation for gene insertions/deletions between genes. This makes it possible to address whether a difference between informational genes and non-informational genes is the main contributor to rate variation of lateral gene transfers.

RESULTS

The results show that models incorporating rate variation fit the data better than do constant rate models in many phylogenetic groups. Even though informational genes are less likely to be laterally transferred than non-informational genes, the degree of rate variation for insertions/deletions did not change dramatically and remained high even when informational genes were excluded from the study. This suggests that the variation in rate of insertions/deletions is not due mainly to the simple difference between informational genes and non-informational genes. Among genes that are not classified as informational and among the informational genes themselves, there are still large differences in the rates that these genes are inserted and deleted.

CONCLUSION

While the difference in informational gene rates contributes to rate variation, it is only a small fraction of the variation present; instead, a substantial amount of rate variation for insertions/deletions remains among both informational genes and among non-informational genes.

Ohno's dilemma: evolution of new genes under continuous selection

New genes with novel functions arise by duplication and divergence, but the process poses a problem. After duplication, an extra gene copy must rise to sufficiently high frequency in the population and remain free of common inactivating lesions long enough to acquire the rare mutations that provide a new selectable function. Maintaining a duplicated gene by selection for the original function would restrict the freedom to diverge. (We refer to this problem as Ohno's dilemma). A model is described by which selection continuously favors both maintenance of the duplicate copy and divergence of that copy from the parent gene. Before duplication, the original gene has a trace side activity (the innovation) in addition to its original function. When an altered ecological niche makes the minor innovation valuable, selection favors increases in its level (the amplification), which is most frequently conferred by increased dosage of the parent gene. Selection for the amplified minor function maintains the extra copies and raises the frequency of the amplification in the population. The same selection favors mutational improvement of any of the extra copies, which are not constrained to maintain their original function (the divergence). The rate of mutations (per genome) that improve the new function is increased by the multiplicity of target copies within a genome. Improvement of some copies relaxes selection on others and allows their loss by mutation (becoming pseudogenes). Ultimately one of the extra copies is able to provide all of the new activity.

Selection effects on the positioning of genes and gene structures from the interplay of replication and transcription in bacterial genomes

Bacterial chromosomes are partly shaped by the functional requirements for efficient replication, which lead to strand bias as commonly characterized by the excess of guanines over cytosines in the leading strand. Gene structures are also highly organized within bacterial genomes as a result of such functional constraints, displaying characteristic positioning and structuring along the genome. Here we analyze the gene structures in completely sequenced bacterial chromosomes to observe the positional constraints on gene orientation, length, and codon usage with regard to the positions of replication origin and terminus. Selection on these gene features is different in regions surrounding the terminus of replication from the rest of the genome, but the selection could be either positive or negative depending on the species, and these positional effects are partly attributed to the A-T enrichment near the terminus. Characteristic gene structuring relative to the position of replication origin and terminus is commonly observed among most bacterial species with circular chromosomes, and therefore we argue that the highly organized gene positioning as well as the strand bias should be considered for genomics studies of bacteria.

Horizontal gene transfer in chromalveolates

BACKGROUND

Horizontal gene transfer (HGT), the non-genealogical transfer of genetic material between different organisms, is considered a potentially important mechanism of genome evolution in eukaryotes. Using phylogenomic analyses of expressed sequence tag (EST) data generated from a clonal cell line of a free living dinoflagellate alga Karenia brevis, we investigated the impact of HGT on genome evolution in unicellular chromalveolate protists.

RESULTS

We identified 16 proteins that have originated in chromalveolates through ancient HGTs before the divergence of the genera Karenia and Karlodinium and one protein that was derived through a more recent HGT. Detailed analysis of the phylogeny and distribution of identified proteins demonstrates that eight have resulted from independent HGTs in several eukaryotic lineages.

CONCLUSION

Recurring intra- and interdomain gene exchange provides an important source of genetic novelty not only in parasitic taxa as previously demonstrated but as we show here, also in free-living protists. Investigating the tempo and mode of evolution of horizontally transferred genes in protists will therefore advance our understanding of mechanisms of adaptation in eukaryotes.

Evolution of symbiotic bacteria in the distal human intestine

The adult human intestine contains trillions of bacteria, representing hundreds of species and thousands of subspecies. Little is known about the selective pressures that have shaped and are shaping this community's component species, which are dominated by members of the Bacteroidetes and Firmicutes divisions. To examine how the intestinal environment affects microbial genome evolution, we have sequenced the genomes of two members of the normal distal human gut microbiota, Bacteroides vulgatus and Bacteroides distasonis, and by comparison with the few other sequenced gut and non-gut Bacteroidetes, analyzed their niche and habitat adaptations. The results show that lateral gene transfer, mobile elements, and gene amplification have played important roles in affecting the ability of gut-dwelling Bacteroidetes to vary their cell surface, sense their environment, and harvest nutrient resources present in the distal intestine. Our findings show that these processes have been a driving force in the adaptation of Bacteroidetes to the distal gut environment, and emphasize the importance of considering the evolution of humans from an additional perspective, namely the evolution of our microbiomes.

The total number of microbes that colonize the surfaces of our adult bodies is thought to be ten times greater than the total number of our human cells. Our microbial partners provide us with certain features that we have not had to evolve on our own. In this sense, we should consider ourselves to be a supraorganism whose genetic landscape includes both our own genome as well as the genomes of our resident microbes, and whose physiologic features are a synthesis of human and microbial metabolic traits. The largest collection of microbes resides in our gut, which harbors trillions of bacteria, representing hundreds of species, most falling into two groups—the Bacteroidetes and the Firmicutes. We have sequenced the genomes of two human gut-dwelling Bacteroidetes, and compared their genomes to the genomes of other bacteria that live both inside and outside of our bodies. Our results illustrate that adaptation to the gut habitat is a dynamic process that includes acquisition of genes from other microorganisms. These findings emphasize the importance of including the evolution of “our” microbial genomes when considering the evolution of humans.

Human microbiome evolution was explored by comparing human gut Bacteroidete genomic sequences to available data; common modes of evolution were revealed that have enabled these gut-dwelling microbes to adapt to their environments.

Erosion of interaction networks in reduced and degraded genomes

Unlike eukaryotes, which often recruit duplicated genes into existing protein-protein interaction (PPI) networks, the low levels of gene duplication coupled with the high probability of lateral transfer of novel genes alters the manner in which PPI networks can evolve in bacteria. By inferring the PPIs present in the ancestor to contemporary Gammaproteobacteria, we were able to trace the changes in gene repertoires, and their consequences on PPI network evolution, in several bacterial lineages that have independently undergone reductions in genome size and genome contents. As genomes degrade, virtually all multi-partner proteins have lost interactors; however, the overall average number of connections increases due to the preferential elimination of proteins that interact with only one other protein partner. We also studied the effect of lateral gene transfer on PPI network evolution by analyzing the connectivity of genes that have been gained along the Escherichia coli lineage, as well as those acquired genes subsequently silenced in Shigella flexneri, since diverging from the gammaproteobacterial ancestor. The situation in PPI networks, in which newly acquired genes preferentially attach to the hubs of the network, contrasts that observed in metabolic networks, which evolve by the peripheral gain and loss of genes, and in regulatory networks, in which high connectivity increases the propensity of loss.

Both α-haemolysin determinants contribute to full virulence of uropathogenic Escherichia coli strain 536

Uropathogenic Escherichia coli strain 536 possesses two intact copies of the alpha-haemolysin determinant localised on distinct pathogenicity islands. The coding regions of the two hlyCABD operons are conserved; however, upstream sequences are entirely dissimilar. Consequently, expression of the encoded toxin molecules in vitro is highly different. On the other hand, the contribution of the individual determinants to the strain's virulence is the same. Isogenic mutants lacking individual hly determinants have a similar increase in LD50 value in a mouse model of urinary tract infection. Mouse lung toxicity as well as in vitro assays reveals a significant decrease in acute cytotoxicity of both mutants in comparison to the parent wild-type strain; however, the two hly mutants do not significantly differ from each other in these respects. Single channel recordings show no difference in electrophysiological characteristics of the pores formed by the individual HlyA molecules on synthetic planar lipid membranes. Nor do the paralogues have any target cell preference in an in vitro cytotoxicity assay. Our data suggest that the two hly paralogues encode identical toxin functions; however, due to different regulation of expression, they participate at distinct stages of the infectious process. Interestingly, the unrelated uropathogenic E. coli strain J96 shares the same two hly alleles, suggesting that acquisition of the two paralogues accorded a selective evolutionary advantage.

Testing of a whole genome PCR scanning approach to identify genomic variability in four different species of lactic acid bacteria

Genomes can be markedly heterogeneous in conspecific bacterial strains. Genome sequences can be used to analyze genome plasticity via a PCR(2) (plasticity of chromosome revealed by PCR) approach. Small-sized chromosomes can indeed be fully amplified by long-range PCR with a set of primers designed using a reference strain and then applied to several other strains. Analysis of the resulting patterns can reveal genome plasticity. GenoFrag, a software package for the design of primers optimized for PCR(2) [N. Ben Zakour, M. Gautier, R. Andonov, D. Lavenier, M.F. Cochet, P. Veber, A. Sorokin, Y. Le Loir, GenoFrag: Software to design primers optimized for whole genome scanning by long-range PCR amplification, Nucleic Acids Res. 32 (2004) 17-24] was developed for the analysis of bacterial genome plasticity by whole genome amplification in approximately 10-kb-long fragments. By applying GenoFrag, we provide herewith evidence that genome plasticity can be analyzed in lactic acid bacteria using a PCR(2) approach. The genome sequences of Lactococcus lactis IL1403, Lactobacillus plantarum WCFS1, Lactobacillus bulgaricus ATCC11842 and Bifidobacterium longum NCC2705 were used to design four sets of primers. Each set was evaluated in silico to check that it ensured optimum coverage of the bacterial chromosome. To validate the primers generated by GenoFrag, a subset of primers was successfully used in LR-PCR experiments on genomic DNA from four L. bulgaricus strains.

The microbial pan-genome

A decade after the beginning of the genomic era, the question of how genomics can describe a bacterial species has not been fully addressed. Experimental data have shown that in some species new genes are discovered even after sequencing the genomes of several strains. Mathematical modeling predicts that new genes will be discovered even after sequencing hundreds of genomes per species. Therefore, a bacterial species can be described by its pan-genome, which is composed of a "core genome" containing genes present in all strains, and a "dispensable genome" containing genes present in two or more strains and genes unique to single strains. Given that the number of unique genes is vast, the pan-genome of a bacterial species might be orders of magnitude larger than any single genome.

Horizontal gene transfer, genome innovation and evolution

To what extent is the tree of life the best representation of the evolutionary history of microorganisms? Recent work has shown that, among sets of prokaryotic genomes in which most homologous genes show extremely low sequence divergence, gene content can vary enormously, implying that those genes that are variably present or absent are frequently horizontally transferred. Traditionally, successful horizontal gene transfer was assumed to provide a selective advantage to either the host or the gene itself, but could horizontally transferred genes be neutral or nearly neutral? We suggest that for many prokaryotes, the boundaries between species are fuzzy, and therefore the principles of population genetics must be broadened so that they can be applied to higher taxonomic categories.

Horizontal transfer of two operons coding for hydrogenases between bacteria and archaea

Using a phylogenetic approach, we discovered three putative horizontal transfers between bacterial and archaeal species involving large clusters of genes. One transfer involves an operon of 13 genes, called mbx, which probably was transferred into the genome of Thermotoga maritima from a species belonging or close to the Pyrococcus genus. The two others implied an operon of six genes, called ech, transferred independently to the genomes of Thermoanaerobacter tengcongensis and Desulfovibrio gigas, from a species belonging or close to the Methanosarcina genus. All these transfers affected operons coding for multisubunit membrane-bound (NiFe) hydrogenases involved in the energy metabolism of the donor genomes. The functionality of the transferred operons has not been experimentally demonstrated for T. maritima, whereas in D. gigas and T. tengcongensis the encoded multisubunit hydrogenase could have a role in energy conservation. This report adds several cases of horizontal gene transfers among hydrogenases already described.

An adaptive radiation model for the origin of new gene functions

Current models for the evolution of new gene functions after gene duplication presume that the duplication events themselves have no effect on fitness. But those duplications that result in new gene functions are likely to be positively selected from their inception. The evolution of new function may start with the amplification of an existing gene with some level of preadaptation for that function, followed by a period of competitive evolution among the gene copies, resulting in the preservation of the most effective variant and the 'pseudogenization' and eventual loss of the rest.

Selection for gene clustering by tandem duplication

In prokaryotic genomes, related genes are frequently clustered in operons and higher-order arrangements that reflect functional context. Organization emerges despite rearrangements that constantly shuffle gene and operon order. Evidence is presented that the tandem duplication of related genes acts as a driving evolutionary force in the origin and maintenance of clusters. Gene amplification can be viewed as a dynamic and reversible regulatory mechanism that facilitates adaptation to variable environments. Clustered genes confer selective benefits via their ability to be coamplified. During evolution, rearrangements that bring together related genes can be selected if they increase the fitness of the organism in which they reside. Similarly, the benefits of gene amplification can prevent the dispersal of existing clusters. Examples of frequent and spontaneous amplification of large genomic fragments are provided. The possibility is raised that tandem gene duplication works in concert with horizontal gene transfer as interrelated evolutionary forces for gene clustering.

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.

Inter-genomic displacement via lateral gene transfer of bacterial trp operons in an overall context of vertical genealogy

Background

The growing conviction that lateral gene transfer plays a significant role in prokaryote genealogy opens up a need for comprehensive evaluations of gene-enzyme systems on a case-by-case basis. Genes of tryptophan biosynthesis are frequently organized as whole-pathway operons, an attribute that is expected to facilitate multi-gene transfer in a single step. We have asked whether events of lateral gene transfer are sufficient to have obscured our ability to track the vertical genealogy that underpins tryptophan biosynthesis.

Results

In 47 complete-genome Bacteria, the genes encoding the seven catalytic domains that participate in primary tryptophan biosynthesis were distinguished from any paralogs or xenologs engaged in other specialized functions. A reliable list of orthologs with carefully ascertained functional roles has thus been assembled and should be valuable as an annotation resource. The protein domains associated with primary tryptophan biosynthesis were then concatenated, yielding single amino-acid sequence strings that represent the entire tryptophan pathway. Lateral gene transfer of several whole-pathway trp operons was demonstrated by use of phylogenetic analysis. Lateral gene transfer of partial-pathway trp operons was also shown, with newly recruited genes functioning either in primary biosynthesis (rarely) or specialized metabolism (more frequently).

Conclusions

(i) Concatenated tryptophan protein trees are congruent with 16S rRNA subtrees provided that the genomes represented are of sufficiently close phylogenetic spacing. There are currently seven tryptophan congruency groups in the Bacteria. Recognition of a succession of others can be expected in the near future, but ultimately these should coalesce to a single grouping that parallels the 16S rRNA tree (except for cases of lateral gene transfer). (ii) The vertical trace of evolution for tryptophan biosynthesis can be deduced. The daunting complexities engendered by paralogy, xenology, and idiosyncrasies of nomenclature at this point in time have necessitated an expert-assisted manual effort to achieve a correct analysis. Once recognized and sorted out, paralogy and xenology can be viewed as features that enrich evolutionary histories.