Ancient lateral gene transfer in the evolution of Bdellovibrio bacteriovorus

The recently sequenced genome of the predatory delta-proteobacterium Bdellovibrio bacteriovorus provides many insights into its metabolism and evolution. Because its genes are reasonably uniform in G+C content, it was suggested that B. bacteriovorus actively resists recombination with foreign DNA and horizontal transfer of DNA from other bacteria. To investigate this further, we carried out a variety of phylogenetic and comparative genomics analyses using data from >200 microbial genomes, including several published delta-proteobacteria. Although there might be little evidence for the extensive recent transfer of genes, we demonstrate that ancient lateral gene acquisition has shaped the B. bacteriovorus genome to a great extent.     

A "chimera" theory on the origin of dicyemid mesozoans: evolution driven by frequent lateral gene transfer from host to parasite

The phylogenetic status of the enigmatic dicyemid mesozoans is still uncertain. Are they primitive multicellular organisms or degenerate triploblastic animals? Presently, the latter view is accepted. A phylogenetic analysis of 18S rDNA sequences placed dicyemids within the animal clade, and this was supported by the discovery of a Hox-type gene with a lophotrochozoan signature sequence. This molecular information suggests that dicyemid mesozoans evolved from an ancestral animal degenerately. Considering their extreme simplicity, which is probably due to parasitism, they might have come from an early embryo via a radical transformation, i.e. neoteny. Irrespective of this molecular information, dicyemid mesozoans retain many protistan-like or extremely primitive features, such as tubular mitochondrial cristae, endocytic ability from the outer surface, and the absence of collagenous tissue, while they do not share noticeable synapomorphy with animals. In addition, the 5S rRNA phylogeny suggests a somewhat closer kinship with protozoan ciliates than with animals. If we accept this clear contradiction, dicyemids should be regarded as a chimera of animals and protistans. Here, we discuss the traditional theory of extreme degeneration via parasitism, and then propose a new "chimera" theory in which dicyemid mesozoans are exposed to a continual flow of genetic information via eating host tissues from the outer surface by endocytosis. Consequently, many of their intrinsic genes have been replaced by host-derived genes through lateral gene transfer (LGT), implying that LGT is a key driving force in the evolution of dicyemid mesozoans.     

Development of a gene transfer system for curing of plasmids in the marine fish pathogen Vibrio salmonicida.

All reported natural isolates of the marine fish pathogen Vibrio salmonicida contain plasmids, and in another marine fish pathogen, Vibrio anguillarum, it has been shown that a plasmid is important for expression of virulence by the organism. To study the function of the plasmids in V. salmonicida, we developed a gene transfer system based on the plasmid RSF1010 replicon. The gene transfer system was used to construct a plasmid-free strain, and this strain was found to behave similarly to the wild type in a fish pathogenicity test based on intraperitoneal injection of the bacteria. We were unable to detect any other phenotypic differences between the two strains. It could therefore be concluded that at least in the V. salmonicida strain tested, extrachromosomal DNA is not required for expression of virulence.     

Development of a gene transfer system for curing of plasmids in the marine fish pathogen Vibrio salmonicida.

All reported natural isolates of the marine fish pathogen Vibrio salmonicida contain plasmids, and in another marine fish pathogen, Vibrio anguillarum, it has been shown that a plasmid is important for expression of virulence by the organism. To study the function of the plasmids in V. salmonicida, we developed a gene transfer system based on the plasmid RSF1010 replicon. The gene transfer system was used to construct a plasmid-free strain, and this strain was found to behave similarly to the wild type in a fish pathogenicity test based on intraperitoneal injection of the bacteria. We were unable to detect any other phenotypic differences between the two strains. It could therefore be concluded that at least in the V. salmonicida strain tested, extrachromosomal DNA is not required for expression of virulence.     

The role of lateral gene transfer in the evolution of isoprenoid biosynthesis pathways

Lateral gene transfer (LGT) is a major force in microbial genome evolution. Here, we present an overview of lateral transfers affecting genes involved in isopentenyl diphosphate (IPP) synthesis. Two alternative metabolic pathways can synthesize this universal precursor of isoprenoids, the 1-deoxy-D-xylulose 5-phosphate (DOXP) pathway and the mevalonate (MVA) pathway. We have surveyed recent genomic data and the biochemical literature to determine the distribution of the genes composing these pathways within the bacterial domain. The scattered distribution observed is incompatible with a simple scheme of vertical transmission. LGT (among and between bacteria, archaea and eukaryotes) more parsimoniously explains many features of this pattern. This alternative scenario is supported by phylogenetic analyses, which unambiguously confirm several cases of lateral transfer. Available biochemical data allow the formulation of hypotheses about selective pressures favouring transfer. The phylogenetic diversity of the organisms involved and the range of possible causes and effects of these transfer events make the IPP biosynthetic pathways an ideal system for studying the evolutionary role of LGT.     

Diverse genomic species and evidences of symbiotic gene lateral transfer detected among the rhizobia associated with Astragalus species grown in the temperate regions of China

Based on the analyses of ribosomal DNA and housekeeping genes, a total of 118 bacterial isolates obtained from 13 Astragalus species grown in the temperate region of China were identified as 19 genomic species of Mesorhizobium, Rhizobium, Sinorhizobium and Bradyrhizobium, two of them being putatively new species. Phylogenetic comparison of symbiotic genes (nodC and nifH) and housekeeping genes showed that the symbiotic genes of the Astragalus rhizobia were maintained by both vertical and horizontal transfer. The results demonstrated that the Astragalus species were very promiscuous hosts for rhizobia and that their rhizobia had very diverse genomic and symbiotic gene backgrounds.     

Interactions among genomic structure, function, and evolution revealed by comprehensive analysis of the Arabidopsis thaliana genome

The genome in a higher organism consists of a number of types of nucleotide sequence-specialized components, with each having tens of thousands of members or elements. It is crucial for our understanding of how a genome as an entity is organized, functions, and evolves to determine how these components are organized in the genome and how they relate with each other; however, no such knowledge is available. Here, we report a comprehensive analysis of the organization and interaction of all 40 components constituting the genome of the plant model species, Arabidopsis thaliana, at the whole-genome and chromosome levels. The 40 components include (i) 6 genome structural components consisting of GC%, genes, retrotransposons, DNA transposons, simple repeats, and low complex repeats; (ii) 3 evolutionarily critical features consisting of recombination rate, nucleotide substitutions, and nucleotide insertions/deletions; and (iii) 31 categories of genes with different functions and numbers of functions. We show that the distributions of 39 of the 40 components of the genome (excepting GC%) deviate significantly from the random distribution model and different types of the genome components are significantly correlated. These results remained to be true even when the genomic regions, such as centromeric regions, where transposable and repeat elements are abundant were excluded from the analyses. These findings suggest that DNA molecules contained in the Arabidopsis genome are each organized and structured from their constituting components in an unambiguous manner and that different types of the components that constitute or characterize the genome interact. The analysis also showed that each chromosome consists of a similar set of the components at similar densities, suggesting that the unique organization and interaction pattern of the components in each chromosome may represent, at least in part, the identity of a chromosome or a genome at the genome level, thus partly accounting for the phenotypic variation among different species. The data also provide comprehensive and new insights into many phenomena significant in genome biology, with which we particularly discuss the variation of genetic recombination. The variation of genetic recombination rate along a chromosomal arm is shaped, not only by the distribution of simple repeats, retrotransposons, DNA transposons, and nucleotide substitutions, but also by the functions of genes contained, especially those with multiple functions, suggesting that variation of genetic recombination along a chromosomal arm is the result of interactions among the components constituting local genome structure, function, and evolution.     

The genome of Methanosarcina mazei: evidence for lateral gene transfer between bacteria and archaea

The Archaeon Methanosarcina mazei and related species are of great ecological importance as they are the only organisms fermenting acetate, methylamines and methanol to methane, carbon dioxide and ammonia (in case of methylamines). Since acetate is the precursor of 60% of the methane produced on earth these organisms contribute significantly to the production of this greenhouse gas, e.g. in rice paddies. The 4,096,345 base pairs circular chromosome of M. mazei is more than twice as large as the genomes of the methanogenic Archaea currently completely sequenced (Bult et al., 1996; Smith et al., 1997). 3,371 open reading frames (ORFs) were identified. Based on currently available sequence data 376 of these ORFs are Methanosarcina-specific and 1,043 ORFs find their closest homologue in the bacterial domain. 544 of these ORFs reach significant similarity values only in the bacterial domain. They include 56 of the 102 transposases, and proteins involved in gluconeogenesis, proline biosynthesis, transport processes, DNA-repair, environmental sensing, gene regulation, and stress response. Striking examples are the occurrence of the bacterial GroEL/GroES chaperone system and the presence of tetrahydrofolate-dependent enzymes. These findings might indicate that lateral gene transfer has played an important evolutionary role in forging the physiology of this metabolically versatile methanogen.     

A chromosome-scale Gastrodia elata genome and large-scale comparative genomic analysis indicate convergent evolution by gene loss in mycoheterotrophic and parasitic plants

Mycoheterotrophic and parasitic plants are heterotrophic and parasitize on fungi and plants, respectively, to obtain nutrients. Large-scale comparative genomics analysis has not been conducted in mycoheterotrophic or parasitic plants or between these two groups of parasites. We assembled a chromosome-level genome of the fully mycoheterotrophic plant Gastrodia elata (Orchidaceae) and performed comparative genomic analyses on the genomes of G. elata and four orchids (initial mycoheterotrophs), three parasitic plants (Cuscuta australis, Striga asiatica, and Sapria himalayana), and 36 autotrophs from various angiosperm lineages. It was found that while in the hemiparasite S. asiatica and initial mycoheterotrophic orchids, approximately 4–5% of the conserved orthogroups were lost, the fully heterotrophic G. elata and C. australis both lost approximately 10% of the conserved orthogroups, indicating that increased heterotrophy is positively associated with gene loss. Importantly, many genes that are essential for autotrophs, including those involved in photosynthesis, the circadian clock, flowering time regulation, immunity, nutrient uptake, and root and leaf development, were convergently lost in both G. elata and C. australis. The high-quality genome of G. elata will facilitate future studies on the physiology, ecology, and evolution of mycoheterotrophic plants, and our findings highlight the critical role of gene loss in the evolution of plants with heterotrophic lifestyles.     

Gene duplication,transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.     

Archaeology and evolution of transfer RNA genes in the Escherichia coli genome

Transfer RNA genes tend to be presented in multiple copies in the genomes of most organisms, from bacteria to eukaryotes. The evolution and genomic structure of tRNA genes has been a somewhat neglected area of molecular evolution. Escherichia coli, the first phylogenetic species for which more than two different strains have been sequenced, provides an invaluable framework to study the evolution of tRNA genes. In this work, a detailed analysis of the tRNA structure of the genomes of Escherichia coli strains K12, CFT073, and O157:H7, Shigella flexneri 2a 301, and Salmonella typhimurium LT2 was carried out. A phylogenetic analysis of these organisms was completed, and an archaeological map depicting the main events in the evolution of tRNA genes was drawn. It is shown that duplications, deletions, and horizontal gene transfers are the main factors driving tRNA evolution in these genomes. On average, 0.64 tRNA insertions/duplications occur every million years (Myr) per genome per lineage, while deletions occur at the slower rate of 0.30 per million years per genome per lineage. This work provides a first genomic glance at the problem of tRNA evolution as a repetitive process, and the relationship of this mechanism to genome evolution and codon usage is discussed.     

Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis.

Myosin heavy chain genes (MYHs) are the most important functional domains of myosins, which are highly conserved throughout evolution. The human genome contains 15 MYHs, whereas the corresponding number in teleost appears to be much higher. Although teleosts comprise more than one-half of all vertebrate species, our knowledge of MYHs in teleosts is rather limited. A comprehensive analysis of the torafugu (Takifugu rubripes) genome database enabled us to detect at least 28 MYHs, almost twice as many as in humans. RT-PCR revealed that at least 16 torafugu MYH representatives (5 fast skeletal, 3 cardiac, 2 slow skeletal, 1 superfast, 2 smooth, and 3 nonmuscle types) are actually transcribed. Among these, MYH(M743-2) and MYH(M5) of fast and slow skeletal types, respectively, are expressed during development of torafugu embryos. Syntenic analysis reveals that torafugu fast skeletal MYHs are distributed across five genomic regions, three of which form clusters. Interestingly, while human fast skeletal MYHs form one cluster, its syntenic region in torafugu is duplicated, although each locus contains just a single MYH in torafugu. The results of the syntenic analysis were further confirmed by corresponding analysis of MYHs based on databases from Tetraodon, zebrafish, and medaka genomes. Phylogenetic analysis suggests that fast skeletal MYHs evolved independently in teleosts and tetrapods after fast skeletal MYHs had diverged from four ancestral MYHs.     

Gene and genome duplications in vertebrates: the one-to-four (-to-eight in fish) rule and the evolution of novel gene functions

One important mechanism for functional innovation during evolution is the duplication of genes and entire genomes. Evidence is accumulating that during the evolution of vertebrates from early deuterostome ancestors entire genomes were duplicated through two rounds of duplications (the 'one-to-two-to-four' rule). The first genome duplication in chordate evolution might predate the Cambrian explosion. The second genome duplication possibly dates back to the early Devonian. Recent data suggest that later in the Devonian, the fish genome was duplicated for a third time to produce up to eight copies of the original deuterostome genome. This last duplication took place after the two major radiations of jawed vertebrate life, the ray-finned fish (Actinopterygia) and the sarcopterygian lineage, diverged. Therefore the sarcopterygian fish, which includes the coelacanth, lungfish and all land vertebrates such as amphibians, reptiles, birds and mammals, tend to have only half the number of genes compared with actinopterygian fish. Although many duplicated genes turned into pseudogenes, or even 'junk' DNA, many others evolved new functions particularly during development. The increased genetic complexity of fish might reflect their evolutionary success and diversity.     

Environmental distribution of the trichloroethene reductive dehalogenase gene (tceA) suggests lateral gene transfer among Dehalococcoides

The trichloroethene reductive dehalogenase gene (tceA) of Dehalococcoides spp. was detected in 12 of 21 trichloroethene-to-ethene dechlorinating enrichment cultures established from aquifer and sediment samples collected from diverse geographic locations in the USA. Analysis of the tceA chromosomal regions indicated that the tceA genes shared greater than 95% sequence identity, and all shared identical tceAB spacer sequences and tceB genes downstream of tceA. A putative transposable element (PTE) was present 1077 bp downstream of the tceB stop codon in three of eight chromosomal regions analyzed. Sequence identity was interrupted downstream of tceB and upstream or downstream of the PTE, suggesting that intrachromosomal or interchromosomal transfer of tceAB had occurred.     

Early evolutionary history and genomic features of gene duplicates in the human genome

Background

Human gene duplicates have been the focus of intense research since the development of array-based and targeted next-generation sequencing approaches in the last decade. These studies have primarily concentrated on determining the extant copy-number variation from a population-genomic perspective but lack a robust evolutionary framework to elucidate the early structural and genomic characteristics of gene duplicates at emergence and their subsequent evolution with increasing age.

Results

We analyzed 184 gene duplicate pairs comprising small gene families in the draft human genome with 10 % or less synonymous sequence divergence. Human gene duplicates primarily originate from DNA-mediated events, taking up genomic residence as intrachromosomal copies in direct or inverse orientation. The distribution of paralogs on autosomes follows random expectations in contrast to their significant enrichment on the sex chromosomes. Furthermore, human gene duplicates exhibit a skewed gradient of distribution along the chromosomal length with significant clustering in pericentromeric regions. Surprisingly, despite the large average length of human genes, the majority of extant duplicates (83 %) are complete duplicates, wherein the entire ORF of the ancestral copy was duplicated. The preponderance of complete duplicates is in accord with an extremely large median duplication span of 36 kb, which enhances the probability of capturing ancestral ORFs in their entirety. With increasing evolutionary age, human paralogs exhibit declines in (i) the frequency of intrachromosomal paralogs, and (ii) the proportion of complete duplicates. These changes may reflect lower survival rates of certain classes of duplicates and/or the role of purifying selection. Duplications arising from RNA-mediated events comprise a small fraction (11.4 %) of all human paralogs and are more numerous in older evolutionary cohorts of duplicates.

Conclusions

The degree of structural resemblance, genomic location and duplication span appear to influence the long-term maintenance of paralogs in the human genome. The median duplication span in the human genome far exceeds that in C. elegans and yeast and likely contributes to the high prevalence of complete duplicates relative to structurally heterogeneous duplicates (partial and chimeric). The relative roles of regulatory sequence versus exon-intron structure changes in the acquisition of novel function by human paralogs remains to be determined.

Electronic supplementary material

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