Lateral gene transfer in Mycobacterium avium subspecies paratuberculosis

Lateral gene transfer is an integral part of genome evolution in most bacteria. Bacteria can readily change the contents of their genomes to increase adaptability to ever-changing surroundings and to generate evolutionary novelty. Here, we report instances of lateral gene transfer in Mycobacterium avium subsp. paratuberculosis, a pathogenic bacteria that causes Johne's disease in cattle. A set of 275 genes are identified that are likely to have been recently acquired by lateral gene transfer. The analysis indicated that 53 of the 275 genes were acquired after the divergence of M. avium subsp. paratuberculosis from M. avium subsp. avium, whereas the remaining 222 genes were possibly acquired by a common ancestor of M. avium subsp. paratuberculosis and M. avium subsp. avium after its divergence from the ancestor of M. tuberculosis complex. Many of the acquired genes were from proteobacteria or soil dwelling actinobacteria. Prominent among the predicted laterally transferred genes is the gene rsbR, a possible regulator of sigma factor, and the genes designated MAP3614 and MAP3757, which are similar to genes in eukaryotes. The results of this study suggest that like most other bacteria, lateral gene transfers seem to be a common feature in M. avium subsp. paratuberculosis and that the proteobacteria contribute most of these genetic exchanges.     

Innovation from reduction: gene loss, domain loss and sequence divergence in genome evolution

Analyses of genome sequences have revealed a surprisingly variable distribution of genes, reflecting the generation of novel genes, lateral gene transfer and gene loss. The impact of gene loss on organisms has been difficult to examine, but the loss of protein coding genes, the loss of domains within proteins and the divergence of genes have made surprising contributions to the differences among organisms. This paper reviews surveys of gene loss and divergence in fungal and archaeal genomes that indicate suites of functionally related genes tend to undergo loss and divergence. Instances of fungal gene loss highlighted here suggest that specific cellular systems have changed, such as Ca 2+ biology in Saccharomyces cerevisiae and peroxisome function in Schizosaccharomyces pombe. Analyses of loss and divergence can provide specific predictions regarding protein-protein interactions, and the relationship between networks of protein interactions and loss may form a part of a parametric model of genome evolution.     

Simulated evolution of protein-protein interaction networks with realistic topology

We model the evolution of eukaryotic protein-protein interaction (PPI) networks. In our model, PPI networks evolve by two known biological mechanisms: (1) Gene duplication, which is followed by rapid diversification of duplicate interactions. (2) Neofunctionalization, in which a mutation leads to a new interaction with some other protein. Since many interactions are due to simple surface compatibility, we hypothesize there is an increased likelihood of interacting with other proteins in the target protein's neighborhood. We find good agreement of the model on 10 different network properties compared to high-confidence experimental PPI networks in yeast, fruit flies, and humans. Key findings are: (1) PPI networks evolve modular structures, with no need to invoke particular selection pressures. (2) Proteins in cells have on average about 6 degrees of separation, similar to some social networks, such as human-communication and actor networks. (3) Unlike social networks, which have a shrinking diameter (degree of maximum separation) over time, PPI networks are predicted to grow in diameter. (4) The model indicates that evolutionarily old proteins should have higher connectivities and be more centrally embedded in their networks. This suggests a way in which present-day proteomics data could provide insights into biological evolution.     

Transfer of rpl22 to the nucleus greatly preceded its loss from the chloroplast and involved the gain of an intron.

Most chloroplast and mitochondrial proteins are encoded by nuclear genes that once resided in the organellar genomes. Transfer of most of these genes appears to have occurred soon after the endosymbiotic origin of organelles, and so little is known about the process. Our efforts to understand how chloroplast genes are functionally transferred to the nuclear genome have led us to discover the most recent evolutionary gene transfer yet described. The gene rpl22, encoding chloroplast ribosomal protein CL22, is present in the chloroplast genome of all plants examined except legumes, while a functional copy of rpl22 is located in the nucleus of the legume pea. The nuclear rpl22 gene has acquired two additional domains relative to its chloroplast ancestor: an exon encoding a putative N-terminal transit peptide, followed by an intron which separates this first exon from the evolutionarily conserved, chloroplast-derived portion of the gene. This gene structure suggests that the transferred region may have acquired its transit peptide by a form of exon shuffling. Surprisingly, phylogenetic analysis shows that rpl22 was transferred to the nucleus in a common ancestor of all flowering plants, at least 100 million years preceding its loss from the legume chloroplast lineage.     

The Journey to smORFland

The genome sequences completed so far contain more than 20 000 genes with unknown function and no similarity to genes in other genomes. The origin and evolution of the orphan genes is an enigma. Here, we discuss the suggestion that some orphan genes may represent pseudogenes or short fragments of genes that were functional in the genome of a common ancestor. These may be the remains of unsuccessful duplication or horizontal gene transfer events, in which the acquired sequences have entered the fragmentation process and thereby lost their similarity to genes in other species. This scenario is supported by a recent case study of orphan genes in several closely related species of Rickettsia, where full-length ancestral genes were reconstructed from sets of short, overlapping orphan genes. One of these was found to display similarity to genes encoding proteins with ankyrin-repeat domains.     

Eukaryotic origin of glyceraldehyde-3-phosphate dehydrogenase genes in Clostridium thermocellum and Clostridium cellulolyticum genomes and putative fates of the exogenous gene in the subsequent genome evolution

Although lateral gene transfer (LGT) events have been frequently documented in the evolution of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), no eukaryote-to-prokaryote transfer has been reported so far. Here we describe the first case of the GAPDH gene transfer from a eukaryote to a subset of Clostridium species (Bacteria, Firmicutes). A series of phylogenetic analyses of GAPDH homologues revealed that Clostridium thermocellum and Clostridium cellulolyticum homologues have the evolutionary affinity to the eukaryotic homologues, rather than to those of bacterial species closely related to the two Clostridium species in the organismal phylogeny. These results suggest that the GAPDH genes in the two Clostridium species are of eukaryotic origin, which is the first reported case of eukaryote-to-bacterium GAPDH gene transfer. Since a previously published 16S ribosomal DNA phylogeny and our GAPDH phylogeny commonly suggest an intimate evolutionary relationship between C. thermocellum and C. cellulolyticum, a common ancestor of the two species likely acquired the eukaryotic GAPDH gene. In the C. cellulolyticum genome, the exogenous GAPDH gene was physically separated from other glycolytic genes, suggesting that this gene organization was likely achieved by a random insertion of the laterally transferred gene. On the other hand, in the C. thermocellum genome, the laterally transferred GAPDH gene clusters with other bacterial glycolytic genes. We discuss possible scenarios for the evolutionarily chimeric glycolytic gene cluster in the C. thermocellum genome.     

Gene gain and gene loss in streptococcus: is it driven by habitat?

Bacterial genomes can evolve either by gene gain, gene loss, mutating existing genes, and/or by duplication of existing genes. Recent studies have clearly demonstrated that the acquisition of new genes by lateral gene transfer (LGT) is a predominant force in bacterial evolution. To better understand the significance of LGT, we employed a comparative genomics approach to model species-specific and intraspecies gene insertions/deletions (ins/del among 12 sequenced streptococcal genomes using a maximum likelihood method. This study indicates that the rate of gene ins/del is higher on the external branches and varies dramatically for each species. We have analyzed here some of the experimentally characterized species-specific genes that have been acquired by LGT and conclude that at least a portion of these genes have a role in adaptation.     

Probing the extent of randomness in protein interaction networks

Protein-protein interaction (PPI) networks are commonly explored for the identification of distinctive biological traits, such as pathways, modules, and functional motifs. In this respect, understanding the underlying network structure is vital to assess the significance of any discovered features. We recently demonstrated that PPI networks show degree-weighted behavior, whereby the probability of interaction between two proteins is generally proportional to the product of their numbers of interacting partners or degrees. It was surmised that degree-weighted behavior is a characteristic of randomness. We expand upon these findings by developing a random, degree-weighted, network model and show that eight PPI networks determined from single high-throughput (HT) experiments have global and local properties that are consistent with this model. The apparent random connectivity in HT PPI networks is counter-intuitive with respect to their observed degree distributions; however, we resolve this discrepancy by introducing a non-network-based model for the evolution of protein degrees or "binding affinities." This mechanism is based on duplication and random mutation, for which the degree distribution converges to a steady state that is identical to one obtained by averaging over the eight HT PPI networks. The results imply that the degrees and connectivities incorporated in HT PPI networks are characteristic of unbiased interactions between proteins that have varying individual binding affinities. These findings corroborate the observation that curated and high-confidence PPI networks are distinct from HT PPI networks and not consistent with a random connectivity. These results provide an avenue to discern indiscriminate organizations in biological networks and suggest caution in the analysis of curated and high-confidence networks.     

Preferential Duplication of Intermodular Hub Genes: An Evolutionary Signature in Eukaryotes Genome Networks

Whole genome protein-protein association networks are not random and their topological properties stem from genome evolution mechanisms. In fact, more connected, but less clustered proteins are related to genes that, in general, present more paralogs as compared to other genes, indicating frequent previous gene duplication episodes. On the other hand, genes related to conserved biological functions present few or no paralogs and yield proteins that are highly connected and clustered. These general network characteristics must have an evolutionary explanation. Considering data from STRING database, we present here experimental evidence that, more than not being scale free, protein degree distributions of organisms present an increased probability for high degree nodes. Furthermore, based on this experimental evidence, we propose a simulation model for genome evolution, where genes in a network are either acquired de novo using a preferential attachment rule, or duplicated with a probability that linearly grows with gene degree and decreases with its clustering coefficient. For the first time a model yields results that simultaneously describe different topological distributions. Also, this model correctly predicts that, to produce protein-protein association networks with number of links and number of nodes in the observed range for Eukaryotes, it is necessary 90% of gene duplication and 10% of de novo gene acquisition. This scenario implies a universal mechanism for genome evolution.     

The role of laterally transferred genes in adaptive evolution

Background

Bacterial genomes develop new mechanisms to tide them over the imposing conditions they encounter during the course of their evolution. Acquisition of new genes by lateral gene transfer may be one of the dominant ways of adaptation in bacterial genome evolution. Lateral gene transfer provides the bacterial genome with a new set of genes that help it to explore and adapt to new ecological niches.

Methods

A maximum likelihood analysis was done on the five sequenced corynebacterial genomes to model the rates of gene insertions/deletions at various depths of the phylogeny.

Results

The study shows that most of the laterally acquired genes are transient and the inferred rates of gene movement are higher on the external branches of the phylogeny and decrease as the phylogenetic depth increases. The newly acquired genes are under relaxed selection and evolve faster than their older counterparts. Analysis of some of the functionally characterised LGTs in each species has indicated that they may have a possible adaptive role.

Conclusion

The five Corynebacterial genomes sequenced to date have evolved by acquiring between 8 – 14% of their genomes by LGT and some of these genes may have a role in adaptation.

     

Evolution of the enzymes of the citric acid cycle and the glyoxylate cycle of higher plants. A case study of endosymbiotic gene transfer.

The citric acid or tricarboxylic acid cycle is a central element of higher-plant carbon metabolism which provides, among other things, electrons for oxidative phosphorylation in the inner mitochondrial membrane, intermediates for amino-acid biosynthesis, and oxaloacetate for gluconeogenesis from succinate derived from fatty acids via the glyoxylate cycle in glyoxysomes. The tricarboxylic acid cycle is a typical mitochondrial pathway and is widespread among alpha-proteobacteria, the group of eubacteria as defined under rRNA systematics from which mitochondria arose. Most of the enzymes of the tricarboxylic acid cycle are encoded in the nucleus in higher eukaryotes, and several have been previously shown to branch with their homologues from alpha-proteobacteria, indicating that the eukaryotic nuclear genes were acquired from the mitochondrial genome during the course of evolution. Here, we investigate the individual evolutionary histories of all of the enzymes of the tricarboxylic acid cycle and the glyoxylate cycle using protein maximum likelihood phylogenies, focusing on the evolutionary origin of the nuclear-encoded proteins in higher plants. The results indicate that about half of the proteins involved in this eukaryotic pathway are most similar to their alpha-proteobacterial homologues, whereas the remainder are most similar to eubacterial, but not specifically alpha-proteobacterial, homologues. A consideration of (a) the process of lateral gene transfer among free-living prokaryotes and (b) the mechanistics of endosymbiotic (symbiont-to-host) gene transfer reveals that it is unrealistic to expect all nuclear genes that were acquired from the alpha-proteobacterial ancestor of mitochondria to branch specifically with their homologues encoded in the genomes of contemporary alpha-proteobacteria. Rather, even if molecular phylogenetics were to work perfectly (which it does not), then some nuclear-encoded proteins that were acquired from the alpha-proteobacterial ancestor of mitochondria should, in phylogenetic trees, branch with homologues that are no longer found in most alpha-proteobacterial genomes, and some should reside on long branches that reveal affinity to eubacterial rather than archaebacterial homologues, but no particular affinity for any specific eubacterial donor.     

Lateral gene transfer, lineage-specific gene expansion and the evolution of Nucleo Cytoplasmic Large DNA viruses

Nucleo Cytoplasmic Large DNA viruses (NCLDVs) are a diverse group that infects a wide range of eukaryotic hosts (for example, vertebrates, insects, protists,…) and also show a huge range in genome size (between 100 kb and 1.2 Mb). Here I review some recent results that shed light on the origin and genome evolution of these viruses. Current data suggests that NCLDVs could have originated from a simple and ancient viral ancestor with a small subset of 30-35 genes encoding replication and structural proteins. Subsequent lateral gene transfer of both cellular genes and diverse families of Mobile Genetic Elements, followed by massive lineage-specific gene duplications is probably responsible for the huge diversity of genome size and composition found in extant NCLDVs.     

Segregate or cooperate- a study of the interaction between two species of Dictyostelium

Background

Plastids have inherited their own genomes from a single cyanobacterial ancestor, but the majority of cyanobacterial genes, once retained in the ancestral plastid genome, have been lost or transferred into the eukaryotic host nuclear genome via endosymbiotic gene transfer. Although previous studies showed that cyanobacterial gnd genes, which encode 6-phosphogluconate dehydrogenase, are present in several plastid-lacking protists as well as primary and secondary plastid-containing phototrophic eukaryotes, the evolutionary paths of these genes remain elusive.

Results

Here we show an extended phylogenetic analysis including novel gnd gene sequences from Excavata and Glaucophyta. Our analysis demonstrated the patchy distribution of the excavate genes in the gnd gene phylogeny. The Diplonema gene was related to cytosol-type genes in red algae and Opisthokonta, while heterolobosean genes occupied basal phylogenetic positions with plastid-type red algal genes within the monophyletic eukaryotic group that is sister to cyanobacterial genes. Statistical tests based on exhaustive maximum likelihood analyses strongly rejected that heterolobosean gnd genes were derived from a secondary plastid of green lineage. In addition, the cyanobacterial gnd genes from phototrophic and phagotrophic species in Euglenida were robustly monophyletic with Stramenopiles, and this monophyletic clade was moderately separated from those of red algae. These data suggest that these secondary phototrophic groups might have acquired the cyanobacterial genes independently of secondary endosymbioses.

Conclusion

We propose an evolutionary scenario in which plastid-lacking Excavata acquired cyanobacterial gnd genes via eukaryote-to-eukaryote lateral gene transfer or primary endosymbiotic gene transfer early in eukaryotic evolution, and then lost either their pre-existing or cyanobacterial gene.     

Reconstructing evolution: gene transfer from plastids to the nucleus

During evolution, the genomes of eukaryotic cells have undergone major restructuring to meet the new regulatory challenges associated with compartmentalization of the genetic material in the nucleus and the organelles acquired by endosymbiosis (mitochondria and plastids). Restructuring involved the loss of dispensable or redundant genes and the massive translocation of genes from the ancestral organelles to the nucleus. Genomics and bioinformatic data suggest that the process of DNA transfer from organelles to the nucleus still continues, providing raw material for evolutionary tinkering in the nuclear genome. Recent reconstruction of these events in the laboratory has provided a unique tool to observe genome evolution in real time and to study the molecular mechanisms by which plastid genes are converted into functional nuclear genes. Here, we summarize current knowledge about plastid-to-nuclear gene transfer in the context of genome evolution and discuss new insights gained from experiments that recapitulate endosymbiotic gene transfer in the laboratory.     

Bacterial Protein Interaction Networks: Connectivity is Ruled by Gene Conservation,Essentiality and Function

BackgroundProtein-protein interaction (PPI) networks are the backbone of all processes in living cells. In this work, we relate conservation, essentiality and functional repertoire of a gene to the connectivity k (i.e. the number of interactions, links) of the corresponding protein in the PPI network.MethodsOn a set of 42 bacterial genomes of different sizes, and with reasonably separated evolutionary trajectories, we investigate three issues: i) whether the distribution of connectivities changes between PPI subnetworks of essential and nonessential genes; ii) how gene conservation, measured both by the evolutionary retention index (ERI) and by evolutionary pressures, is related to the connectivity of the corresponding protein; iii) how PPI connectivities are modulated by evolutionary and functional relationships, as represented by the Clusters of Orthologous Genes (COGs).ResultsWe show that conservation, essentiality and functional specialisation of genes constrain the connectivity of the corresponding proteins in bacterial PPI networks. In particular, we isolated a core of highly connected proteins (connectivities k≥40), which is ubiquitous among the species considered here, though mostly visible in the degree distributions of bacteria with small genomes (less than 1000 genes).ConclusionThe genes that support this highly connected core are conserved, essential and, in most cases, belong to the COG cluster J, related to ribosomal functions and the processing of genetic information.     

Vavien: an algorithm for prioritizing candidate disease genes based on topological similarity of proteins in interaction networks

Genome-wide linkage and association studies have demonstrated promise in identifying genetic factors that influence health and disease. An important challenge is to narrow down the set of candidate genes that are implicated by these analyses. Protein-protein interaction (PPI) networks are useful in extracting the functional relationships between known disease and candidate genes, based on the principle that products of genes implicated in similar diseases are likely to exhibit significant connectivity/proximity. Information flow?based methods are shown to be very effective in prioritizing candidate disease genes. In this article, we utilize the topology of PPI networks to infer functional information in the context of disease association. Our approach is based on the assumption that PPI networks are organized into recurrent schemes that underlie the mechanisms of cooperation among different proteins. We hypothesize that proteins associated with similar diseases would exhibit similar topological characteristics in PPI networks. Utilizing the location of a protein in the network with respect to other proteins (i.e., the "topological profile" of the proteins), we develop a novel measure to assess the topological similarity of proteins in a PPI network. We then use this measure to prioritize candidate disease genes based on the topological similarity of their products and the products of known disease genes. We test the resulting algorithm, Vavien, via systematic experimental studies using an integrated human PPI network and the Online Mendelian Inheritance in Man (OMIM) database. Vavien outperforms other network-based prioritization algorithms as shown in the results and is available at www.diseasegenes.org.     

The evolutionary relationships between the two bacteria Escherichia coli and Haemophilus influenzae and their putative last common ancestor

We have tried to approach the nature of the last common ancestor to Haemophilus influenzae and Escherichia coli and to determine how each bacterium could have diverged from this putative organism. The approach used was exhaustive analysis of the homologous proteins coded by genes present in these bacteria, using as criteria for sequence relatedness an alignment of at least 80 amino acid residues and a PAM distance (number of accepted point mutations per 100 residues separating two sequences) below 250. Evolutionarily significant similarities were found between 1,345 H. influenzae proteins (85% of the total genome) and 3,058 E. coli. proteins (75% of the total genome), many of them belonging to families of various sizes (from 666 doublets to 35 large groups of more than 10 members). Nearly all the genes found by this approach to be duplicated in both bacteria were already duplicated in their last common ancestor. This was deduced from (1) the comparison of the respective distributions of evolutionary distances between orthologs (genes separated only by speciation events) and paralogs (genes duplicated in the same genome) and (2) the analysis of the phylogenetic trees reconstructed for each family of paralogs containing at least two members belonging to each bacterium. The distributions of the different categories of homologs show a significant loss of paralogous genes in H. influenzae (reduction proportional to the genome size), of many sequences which are still present in one copy in E. coli, and of some entire gene families. Phylogenetic trees also confirmed this recent loss of paralogous genes in H. influenzae. Thus, the genome size of the last common ancestor of these two bacteria would have been close to that of present-day E. coli, and the evolution of H. influenzae toward a parasitic life led to an important decrease in its genome size by some mechanism of streamlining. During this recent evolution, the memory of the gene order present in the last common ancestor has been blurred, but a few short conserved chromosomal fragments can still be detected in present-day E. coli and H. influenzae.