Genome-scale phylogenetic analysis finds extensive gene transfer among fungi

Although the role of lateral gene transfer is well recognized in the evolution of bacteria, it is generally assumed that it has had less influence among eukaryotes. To explore this hypothesis, we compare the dynamics of genome evolution in two groups of organisms: cyanobacteria and fungi. Ancestral genomes are inferred in both clades using two types of methods: first, Count, a gene tree unaware method that models gene duplications, gains and losses to explain the observed numbers of genes present in a genome; second, ALE, a more recent gene tree-aware method that reconciles gene trees with a species tree using a model of gene duplication, loss and transfer. We compare their merits and their ability to quantify the role of transfers, and assess the impact of taxonomic sampling on their inferences. We present what we believe is compelling evidence that gene transfer plays a significant role in the evolution of fungi.     

SHOT: a web server for the construction of genome phylogenies

With the increasing availability of genome sequences, new methods are being proposed that exploit information from complete genomes to classify species in a phylogeny. Here we present SHOT, a web server for the classification of genomes on the basis of shared gene content or the conservation of gene order that reflects the dominant, phylogenetic signal in these genomic properties. In general, the genome trees are consistent with classical gene-based phylogenies, although some interesting exceptions indicate massive horizontal gene transfer. SHOT is a useful tool for analysing the tree of life from a genomic point of view. It is available at http://www.Bork.EMBL-Heidelberg.de/SHOT.     

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.     

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.     

Run-Off Replication of Host-Adaptability Genes Is Associated with Gene Transfer Agents in the Genome of Mouse-Infecting Bartonella grahamii

The genus Bartonella comprises facultative intracellular bacteria adapted to mammals, including previously recognized and emerging human pathogens. We report the 2,341,328 bp genome sequence of Bartonella grahamii, one of the most prevalent Bartonella species in wild rodents. Comparative genomics revealed that rodent-associated Bartonella species have higher copy numbers of genes for putative host-adaptability factors than the related human-specific pathogens. Many of these gene clusters are located in a highly dynamic region of 461 kb. Using hybridization to a microarray designed for the B. grahamii genome, we observed a massive, putatively phage-derived run-off replication of this region. We also identified a novel gene transfer agent, which packages the bacterial genome, with an over-representation of the amplified DNA, in 14 kb pieces. This is the first observation associating the products of run-off replication with a gene transfer agent. Because of the high concentration of gene clusters for host-adaptation proteins in the amplified region, and since the genes encoding the gene transfer agent and the phage origin are well conserved in Bartonella, we hypothesize that these systems are driven by selection. We propose that the coupling of run-off replication with gene transfer agents promotes diversification and rapid spread of host-adaptability factors, facilitating host shifts in Bartonella.     

Exploring the genomes: From Arabidopsis to crops

Model systems have played a crucial role for understanding biological processes at genetic, molecular and systems levels. Arabidopsis thaliana is one of the best studied model species for higher plants. Large genomic resources and mutant collections made Arabidopsis an excellent source for functional and comparative genomics. Rice and Brachypodium have a great potential to become model systems for grasses. Given the agronomic importance of grass crops, it is an attractive strategy to apply knowledge from Arabidopsis to grasses. Despite many efforts successful reports are sparse. Knowledge transfer should generally work best between orthologous genes that share functionality and a common ancestor. In higher plants, however, recent genome projects revealed an active and rapid evolution of genome structure, which challenges the concept of one-to-one orthologous mates between two species. In this study, we estimated on the example of protein families that are involved in redox related processes, the impact of gene expansions on the success rate for a knowledge transfer from Arabidopsis to the grass species rice, sorghum and Brachypodium. The sparse synteny between dicot and monocot plants due to frequent rearrangements, translocations and gene losses strongly impairs and reduces the number of orthologs detectable by positional conservation. To address the limitations of sparse synteny and expanded gene families, we applied for the detection of orthologs in this study orthoMCL, a sequence-based approach that allows to group closely related paralogs into one orthologous gene cluster. For a total of 49 out of 170 Arabidopsis genes we could identify conserved copy numbers between the dicot model and the grass annotations whereas approximately one third (34.7%, 59 genes) of the selected Arabidopsis genes lack an assignment to any of the grass genome annotations. The remaining 62 Arabidopsis genes represent groups that are considerably biased in their copy numbers between Arabidopsis and all or most of the three grass genomes.     

Gene transfer and genome plasticity in Thermotoga maritima, a model hyperthermophilic species

The genome sequence of the hyperthermophilic bacterium Thermotoga maritima MSB8 presents evidence for lateral gene transfer events between bacterial and archaeal species. To estimate the extent of genomic diversity across the order Thermotogales, a comparative genomic hybridization study was initiated to compare nine Thermotoga strains to the sequenced T. maritima MSB8. Many differences could be associated with substrate utilization patterns, which are most likely a reflection of the environmental niche that these individual species occupy. A detailed analysis of some of the predicted variable regions demonstrates many examples of the deletion/insertion of complete cassettes of genes and of gene rearrangements and insertions of DNA within genes, with the C or N terminus being retained. Although the mechanism for gene transfer in this lineage remains to be elucidated, this analysis suggests possible associations with repetitive elements and highlights the possible benefits of rampant genetic exchange to these species.     

Variation in macronuclear genome content of three ciliates with extensive chromosomal fragmentation: a preliminary analysis

The genome architecture of ciliates, including features such as nuclear dualism and large-scale genome rearrangements, impacts gene and genome evolution in these organisms. To better understand the structure of macronuclear chromosomes in ciliates with extensively processed chromosomes, a sample of complete macronuclear chromosomes was sequenced from three ciliate species: Metopus es (Class [Cl]: Armophorea), Nyctotherus ovalis (Cl: Armophorea), and Chilodonella uncinata (Cl: Phyllopharyngea). By cloning whole macronuclear chromosomes into a plasmid vector, we generated nine clones from each of M. es and C. uncinata, and 37 clones from N. ovalis. Analysis of these macronuclear chromosomes provides insight into the evolution of genome features such as chromosome content, gene structure, and genetic code. Phylogenetic patterns can be found in telomere structure and codon usage, which are both more similar in M. es and N. ovalis than C. uncinata. In addition, we provide evidence of lateral transfer of a bacterial endo-beta-mannanase gene onto a M. es chromosome and report the discovery of a 42-bp conserved sequence motif within N. ovalis untranslated regions.     

Cloning and characterization of a bovine adeno-associated virus

To better understand the relationship between primate adeno-associated viruses (AAVs) and those of other mammals, we have cloned and sequenced the genome of an AAV found as a contaminant in two isolates of bovine adenovirus that was reported to be serologically distinct from primate AAVs. The bovine AAV (BAAV) genome has 4,693 bp, and its organization is similar to that of other AAV isolates. The left-hand open reading frame (ORF) and both inverted terminal repeats (ITRs) have the highest homology with the rep ORF and ITRs of AAV serotype 5 (AAV-5) (89 and 96%, respectively). However, the right-hand ORF was only 55% identical to the AAV-5 capsid ORF; it had the highest homology with the capsid ORF of AAV-4 (76%). By comparing the BAAV cap sequence with a model of an AAV-4 capsid, we mapped the regions of BAAV VP1 that are divergent from AAV-4. These regions are located on the outside of the capsid and are partially located in exposed loops. BAAV was not neutralized by antisera raised against recombinant AAV-2, AAV-4, or AAV-5, and it demonstrated a unique cell tropism profile in four human cancer cell lines, suggesting that BAAV might have transduction activity distinct from that of other isolates. A murine model of salivary gland gene transfer was used to evaluate the in vivo performance of recombinant BAAV. Recombinant BAAV-mediated gene transfer was 11 times more efficient than that with AAV-2. Overall, these data suggest that vectors based on BAAV could be useful for gene transfer applications.     

Modeling genome evolution with a diffusion approximation of a birth-and-death process

MOTIVATION: In our previous studies, we developed discrete-space birth, death and innovation models (BDIMs) of genome evolution. These models explain the origin of the characteristic Pareto distribution of paralogous gene family sizes in genomes, and model parameters that provide for the evolution of these distributions within a realistic time frame have been identified. However, extracting the temporal dynamics of genome evolution from discrete-space BDIM was not technically feasible. We were interested in obtaining dynamic portraits of the genome evolution process by developing a diffusion approximation of BDIM. RESULTS: The diffusion version of BDIM belongs to a class of continuous-state models whose dynamics is described by the Fokker-Plank equation and the stationary solution could be any specified Pareto function. The diffusion models have time-dependent solutions of a special kind, namely, generalized self-similar solutions, which describe the transition from one stationary distribution of the system to another; this provides for the possibility of examining the temporal dynamics of genome evolution. Analysis of the generalized self-similar solutions of the diffusion BDIM reveals a biphasic curve of genome growth in which the initial, relatively short, self-accelerating phase is followed by a prolonged phase of slow deceleration. This evolutionary dynamics was observed both when genome growth started from zero and proceeded via innovation (a potential model of primordial evolution), and when evolution proceeded from one stationary state to another. In biological terms, this regime of evolution can be tentatively interpreted as a punctuated-equilibrium-like phenomenon whereby evolutionary transitions are accompanied by rapid gene amplification and innovation, followed by slow relaxation to a new stationary state.     

Prokaryotic gene clusters: a rich toolbox for synthetic biology

Bacteria construct elaborate nanostructures, obtain nutrients and energy from diverse sources, synthesize complex molecules, and implement signal processing to react to their environment. These complex phenotypes require the coordinated action of multiple genes, which are often encoded in a contiguous region of the genome, referred to as a gene cluster. Gene clusters sometimes contain all of the genes necessary and sufficient for a particular function. As an evolutionary mechanism, gene clusters facilitate the horizontal transfer of the complete function between species. Here, we review recent work on a number of clusters whose functions are relevant to biotechnology. Engineering these clusters has been hindered by their regulatory complexity, the need to balance the expression of many genes, and a lack of tools to design and manipulate DNA at this scale. Advances in synthetic biology will enable the large-scale bottom-up engineering of the clusters to optimize their functions, wake up cryptic clusters, or to transfer them between organisms. Understanding and manipulating gene clusters will move towards an era of genome engineering, where multiple functions can be "mixed-and-matched" to create a designer organism.     

Genomic imprinting of a placental lactogen gene in Peromyscus

The mammalian genome contains over 30 genes whose expression is dependent upon their parent-of-origin. Of these imprinted genes the majority are involved in regulating the rate of fetal growth. In this report we show that in the deer mouse Peromyscusthe placental lactogen-1-variant ( pPl1-v) gene is paternally expressed throughout fetal development, whereas the linked and closely related pPl1gene is expressed in a biallelic manner. Neither the more distantly related pPl2Agene, nor the Mus Pl1gene displays any preferential expression of the paternal allele, suggesting that the acquisition of imprinting of pPl1-v is a relatively recent event in evolution. Although pPl1 expression is temporally mis-regulated in the dysplastic placentae of hybrids between two Peromyscus species, its over-expression cannot account for the aberrant phenotypes of these placentae. We argue that the species-specific imprinting of pPl1-v, encoding a growth factor that regulates nutrient transfer from mothers to their offspring, is consistent with the parent-offspring conflict model that has been proposed to explain the evolution of genomic imprinting.     

Genomic Signatures of Recent Adaptation in a Wild Bumblebee

Environmental changes threaten insect pollinators, creating risks for agriculture and ecosystem stability. Despite their importance, we know little about how wild insects respond to environmental pressures. To understand the genomic bases of adaptation in an ecologically important pollinator, we analyzed genomes of Bombus terrestris bumblebees collected across Great Britain. We reveal extensive genetic diversity within this population, and strong signatures of recent adaptation throughout the genome affecting key processes including neurobiology and wing development. We also discover unusual features of the genome, including a region containing 53 genes that lacks genetic diversity in many bee species, and a horizontal gene transfer from a Wolbachia bacteria. Overall, the genetic diversity we observe and how it is distributed throughout the genome and the population should support the resilience of this important pollinator species to ongoing and future selective pressures. Applying our approach to more species should help understand how they can differ in their adaptive potential, and to develop conservation strategies for those most at risk.     

Application of comparative genomics to narrow-leafed lupin (Lupinus angustifolius L.) using sequence information from soybean and Arabidopsis.

The completion of genome-sequencing initiatives for model plants and EST databases for major crop species provides a large resource for gaining fundamental knowledge of complex gene interactions and the functional significance of proteins. There are increasingly numerous opportunities to transfer this information to other plant species with uncharacterized genomes and make advances in genome analysis, gene expression, and predicted protein function. In this study, we have used DNA sequences from soybean and Arabidopsis to determine the feasibility of applying comparative genomics to narrow-leafed lupin. We have used transcribed sequences from soybean and showed that a high proportion cross hybridize to lupin DNA, identifying similar genes and providing landmarks for estimating the degree of chromosomal synteny between species. To further investigate comparative relationships in this study, a detailed analysis of three lupin genes and comparison of orthologs from soybean and Arabidopsis shows that, in some cases, gene structure and expression are highly conserved and their proteins may have similar function. In other cases, genes show variation in expression profiles indicating alternative functions across species. The advantages and limitation of using soybean and Arabidopsis sequences for comparative genomics in lupins are discussed.     

Patterns of Genomic Integration of Nuclear Chloroplast DNA Fragments in Plant Species

The transfer of organelle DNA fragments to the nuclear genome is frequently observed in eukaryotes. These transfers are thought to play an important role in gene and genome evolution of eukaryotes. In plants, such transfers occur from plastid to nuclear [nuclear plastid DNAs (NUPTs)] and mitochondrial to nuclear (nuclear mitochondrial DNAs) genomes. The amount and genomic organization of organelle DNA fragments have been studied in model plant species, such as Arabidopsis thaliana and rice. At present, publicly available genomic data can be used to conduct such studies in non-model plants. In this study, we analysed the amount and genomic organization of NUPTs in 17 plant species for which genome sequences are available. The amount and distribution of NUPTs varied among the species. We also estimated the distribution of NUPTs according to the time of integration (relative age) by conducting sequence similarity analysis between NUPTs and the plastid genome. The age distributions suggested that the present genomic constitutions of NUPTs could be explained by the combination of the rapidly eliminated deleterious parts and few but constantly existing less deleterious parts.     

Trends and barriers to lateral gene transfer in prokaryotes

Gene acquisition by lateral gene transfer (LGT) is an important mechanism for natural variation among prokaryotes. Laboratory experiments show that protein-coding genes can be laterally transferred extremely fast among microbial cells, inherited to most of their descendants, and adapt to a new regulatory regime within a short time. Recent advance in the phylogenetic analysis of microbial genomes using networks approach reveals a substantial impact of LGT during microbial genome evolution. Phylogenomic networks of LGT among prokaryotes reconstructed from completely sequenced genomes uncover barriers to LGT in multiple levels. Here we discuss the kinds of barriers to gene acquisition in nature including physical barriers for gene transfer between cells, genomic barriers for the integration of acquired DNA, and functional barriers for the acquisition of new genes.     

Identification of a Ligand on the Wip1 Bacteriophage Highly Specific for a Receptor on Bacillus anthracis

Tectiviridae is a family of tailless bacteriophages with Gram-negative and Gram-positive hosts. The family model PRD1 and its close relatives all infect a broad range of enterobacteria by recognizing a plasmid-encoded conjugal transfer complex as a receptor. In contrast, tectiviruses with Gram-positive hosts are highly specific to only a few hosts within the same bacterial species. The cellular determinants that account for the observed specificity remain unknown. Here we present the genome sequence of Wip1, a tectivirus that infects the pathogen Bacillus anthracis. The Wip1 genome is related to other tectiviruses with Gram-positive hosts, notably, AP50, but displays some interesting differences in its genome organization. We identified Wip1 candidate genes for the viral spike complex, the structure located at the capsid vertices and involved in host receptor binding. Phage adsorption and inhibition tests were combined with immunofluorescence microscopy to show that the Wip1 gene product p23 is a receptor binding protein. His-p23 also formed a stable complex with p24, a Wip1 protein of unknown function, suggesting that the latter is involved with p23 in host cell recognition. The narrow host range of phage Wip1 and the identification of p23 as a receptor binding protein offer a new range of suitable tools for the rapid identification of B. anthracis.     

Evolutionary genomics of a temperate bacteriophage in an obligate intracellular bacteria (Wolbachia)

Genome evolution of bacteria is usually influenced by ecology, such that bacteria with a free-living stage have large genomes and high rates of horizontal gene transfer, while obligate intracellular bacteria have small genomes with typically low amounts of gene exchange. However, recent studies indicate that obligate intracellular species that host-switch frequently harbor agents of horizontal transfer such as mobile elements. For example, the temperate double-stranded DNA bacteriophage WO in Wolbachia persistently transfers between bacterial coinfections in the same host. Here we show that despite the phage's rampant mobility between coinfections, the prophage's genome displays features of constraint related to its intracellular niche. First, there is always at least one intact prophage WO and usually several degenerate, independently-acquired WO prophages in each Wolbachia genome. Second, while the prophage genomes are modular in composition with genes of similar function grouping together, the modules are generally not interchangeable with other unrelated phages and thus do not evolve by the Modular Theory. Third, there is an unusual core genome that strictly consists of head and baseplate genes; other gene modules are frequently deleted. Fourth, the prophage recombinases are diverse and there is no conserved integration sequence. Finally, the molecular evolutionary forces acting on prophage WO are point mutation, intragenic recombination, deletion, and purifying selection. Taken together, these analyses indicate that while lateral transfer of phage WO is pervasive between Wolbachia with occasional new gene uptake, constraints of the intracellular niche obstruct extensive mixture between WO and the global phage population. Although the Modular Theory has long been considered the paradigm of temperate bacteriophage evolution in free-living bacteria, it appears irrelevant in phages of obligate intracellular bacteria.     

Substitution of the gene for chloroplast RPS16 was assisted by generation of a dual targeting signal

Organelle (mitochondria and chloroplasts in plants) genomes lost a large number of genes after endosymbiosis occurred. Even after this major gene loss, organelle genomes still lose their own genes, even those that are essential, via gene transfer to the nucleus and gene substitution of either different organelle origin or de novo genes. Gene transfer and substitution events are important processes in the evolution of the eukaryotic cell. Gene loss is an ongoing process in the mitochondria and chloroplasts of higher plants. The gene for ribosomal protein S16 (rps16) is encoded in the chloroplast genome of most higher plants but not in Medicago truncatula and Populus alba. Here, we show that these 2 species have compensated for loss of the rps16 from the chloroplast genome by having a mitochondrial rps16 that can target the chloroplasts as well as mitochondria. Furthermore, in Arabidopsis thaliana, Lycopersicon esculentum, and Oryza sativa, whose chloroplast genomes encode the rps16, we show that the product of the mitochondrial rps16 has dual targeting ability. These results suggest that the dual targeting of RPS16 to the mitochondria and chloroplasts emerged before the divergence of monocots and dicots (140-150 MYA). The gene substitution of the chloroplast rps16 by the nuclear-encoded rps16 in higher plants is the first report about ongoing gene substitution by dual targeting and provides evidence for an intermediate stage in the formation of this heterogeneous organelle.     

Measuring genome divergence in bacteria: a case study using chlamydian data

We have studied the relative contribution of inversions, transpositions, deletions, and nucleotide substitutions to the evolution of Chlamydia trachomatis and Chlamydia pneumoniae. The minimal number of rearrangement events required for converting the gene order structure of one genome into that of the other was estimated to 59 +/- 6 events, including 13% inversions, 38% short inversions, and 49% transpositions. In contrast to previous findings, no examples of horizontal gene transfer subsequent to species divergence were identified, nor any evidence for an excessive number of tandem gene duplications. A statistical model was used to compare nucleotide frequencies for a set of genes uniquely present in one species to a set of orthologous genes present in both species. The two data sets were not significantly different, which is indicative of a low frequency of horizontal gene transfer events. This is based on the assumption that a foreign gene of different nucleotide content will not have become completely ameliorated, as verified by simulations of the amelioration rate at twofold and fourfold degenerate codon sites. The frequencies of nucleotide substitutions at twofold and fourfold degenerate sites, deletions, inversions, and translocations were estimated to 1.42, 0.62, 0.18, 0.01, and 0.01 per site, respectively.