Einat Hazkani-Covo

Directed networks reveal genomic barriers and DNA repair bypasses to lateral gene transfer among prokaryotes

Lateral gene transfer (LGT) plays a major role in prokaryote evolution with only a few genes that are resistant to it; yet the nature and magnitude of barriers to lateral transfer are still debated. Here, we implement directed networks to investigate donor-recipient events of recent lateral gene transfer among 657 sequenced prokaryote genomes. For 2,129,548 genes investigated, we detected 446,854 recent lateral gene transfer events through nucleotide pattern analysis. Among these, donor-recipient relationships could be specified through phylogenetic reconstruction for 7% of the pairs, yielding 32,028 polarized recent gene acquisition events, which constitute the edges of our directed networks. We find that the frequency of recent LGT is linearly correlated both with genome sequence similarity and with proteome similarity of donor-recipient pairs. Genome sequence similarity accounts for 25% of the variation in gene-transfer frequency, with proteome similarity adding only 1% to the variability explained. The range of donor-recipient GC content similarity within the network is extremely narrow, with 86% of the LGTs occurring between donor-recipient pairs having ≤5% difference in GC content. Hence, genome sequence similarity and GC content similarity are strong barriers to LGT in prokaryotes. But they are not insurmountable, as we detected 1530 recent transfers between distantly related genomes. The directed network revealed that recipient genomes of distant transfers encode proteins of nonhomologous end-joining (NHEJ; a DNA repair mechanism) far more frequently than the recipient lacking that mechanism. This implicates NHEJ in genes spread across distantly related prokaryotes through bypassing the donor-recipient sequence similarity barrier.

Endosymbiotic origin and differential loss of eukaryotic genes

Chloroplasts arose from cyanobacteria, mitochondria arose from proteobacteria. Both organelles have conserved their prokaryotic biochemistry, but their genomes are reduced, and most organelle proteins are encoded in the nucleus. Endosymbiotic theory posits that bacterial genes in eukaryotic genomes entered the eukaryotic lineage via organelle ancestors. It predicts episodic influx of prokaryotic genes into the eukaryotic lineage, with acquisition corresponding to endosymbiotic events. Eukaryotic genome sequences, however, increasingly implicate lateral gene transfer, both from prokaryotes to eukaryotes and among eukaryotes, as a source of gene content variation in eukaryotic genomes, which predicts continuous, lineage-specific acquisition of prokaryotic genes in divergent eukaryotic groups. Here we discriminate between these two alternatives by clustering and phylogenetic analysis of eukaryotic gene families having prokaryotic homologues. Our results indicate (1) that gene transfer from bacteria to eukaryotes is episodic, as revealed by gene distributions, and coincides with major evolutionary transitions at the origin of chloroplasts and mitochondria; (2) that gene inheritance in eukaryotes is vertical, as revealed by extensive topological comparison, sparse gene distributions stemming from differential loss; and (3) that continuous, lineage-specific lateral gene transfer, although it sometimes occurs, does not contribute to long-term gene content evolution in eukaryotic genomes.

Evolutionary dynamics of large numts in the human genome: rarity of independent insertions and abundance of post-insertion duplications.

We determined the phylogenetic positions of 82 large nuclear pseudogenes of mitochondrial origin (numts) within the human genome. For each numt, two possibilities pertaining to its origin were considered: (1) independent insertion from the mitochondria into the nucleus, or (2) genomic duplication subsequent to the insertion. A significant increase in the rate of numt accumulation is seen after the divergence of Platyrrhini (New World monkeys) from the Catarrhini (Old World monkeys, apes and humans). By using pairwise phylogenetic analyses, we were able to demonstrate that this peak in numt accumulation is mostly the result of duplication of preexisting nuclear numts rather than the result of an increase in mitochondrial-sequence insertion. In fact, only about a third of all the numt repertoire in the human nuclear genome is due to insertions of mitochondrial sequences, the rest originated as duplications of preexisting numts. Hence, we conclude that numt insertion occurs at a much lower rate than previously reported. As expected under the assumption that genomic duplications occur at rates that are uninfluenced by content, older numts were found to be duplicated more times than recently inserted ones.

Endosymbiotic gene transfer from prokaryotic pangenomes: Inherited chimerism in eukaryotes

Endosymbiotic theory in eukaryotic-cell evolution rests upon a foundation of three cornerstone partners—the plastid (a cyanobacterium), the mitochondrion (a proteobacterium), and its host (an archaeon)—and carries a corollary that, over time, the majority of genes once present in the organelle genomes were relinquished to the chromosomes of the host (endosymbiotic gene transfer). However, notwithstanding eukaryote-specific gene inventions, single-gene phylogenies have never traced eukaryotic genes to three single prokaryotic sources, an issue that hinges crucially upon factors influencing phylogenetic inference. In the age of genomes, single-gene trees, once used to test the predictions of endosymbiotic theory, now spawn new theories that stand to eventually replace endosymbiotic theory with descriptive, gene tree-based variants featuring supernumerary symbionts: prokaryotic partners distinct from the cornerstone trio and whose existence is inferred solely from single-gene trees. We reason that the endosymbiotic ancestors of mitochondria and chloroplasts brought into the eukaryotic—and plant and algal—lineage a genome-sized sample of genes from the proteobacterial and cyanobacterial pangenomes of their respective day and that, even if molecular phylogeny were artifact-free, sampling prokaryotic pangenomes through endosymbiotic gene transfer would lead to inherited chimerism. Recombination in prokaryotes (transduction, conjugation, transformation) differs from recombination in eukaryotes (sex). Prokaryotic recombination leads to pangenomes, and eukaryotic recombination leads to vertical inheritance. Viewed from the perspective of endosymbiotic theory, the critical transition at the eukaryote origin that allowed escape from Muller’s ratchet—the origin of eukaryotic recombination, or sex—might have required surprisingly little evolutionary innovation.

Evolution of multicellularity in Metazoa: comparative analysis of the subcellular localization of proteins in Saccharomyces, Drosophila and Caenorhabditis

A comparison of the subcellular assignments of proteins between the unicellular Saccharomyces cerevisiae and the multicellular Drosophila melanogaster and Caenorhabditis elegans was performed using a computational tool for the prediction of subcellular localization. Nine subcellular compartments were studied: (1) extracellular domain, (2) cell membrane, (3) cytoplasm, (4) endoplasmic reticulum, (5) Golgi apparatus, (6) lysosome, (7) peroxisome, (8) mitochondria, and (9) nucleus. The transition to multicellularity was found to be characterized by an increase in the total number of proteins encoded by the genome. Interestingly, this increase is distributed unevenly among the subcellular compartments. That is, a disproportionate increase in the number of proteins in the extracellular domain, the cell membrane, and the cytoplasm is observed in multicellular organisms, while no such increase is seen in other subcellular compartments.

The Evolutionary History of Prosaposin: Two Successive Tandem-Duplication Events Gave Rise to the Four Saposin Domains in Vertebrates

Abstract. Prosaposin is a multifunctional protein encoded by a single-copy gene. It contains four saposin domains (A, B, C, and D) occurring as tandem repeats connected by linker sequences. Because the saposin domains are similar to one another, it is deduced that they were created by sequential duplications of an ancestral domain. There are two types of evolutionary scenarios that may explain the creation of the four-domain gene: (1) two rounds of tandem internal gene duplication and (2) three rounds of duplications. An evolutionary and phylogenetic analysis of saposin DNA and amino acid sequences from human, mouse, rat, chicken, and zebrafish indicates that the first evolutionary scenario is the most likely. Accordingly, an ancestral saposin-unit duplication produced a two-domain gene, which, subsequently, underwent a second complete tandem duplication to give rise to the present four-domain structure of the prosaposin gene.

Quantifying the Number of Independent Organelle DNA Insertions in Genome Evolution and Human Health

Fragments of organelle genomes are often found as insertions in nuclear DNA. These fragments of mitochondrial DNA (numts) and plastid DNA (nupts) are ubiquitous components of eukaryotic genomes. They are, however, often edited out during the genome assembly process, leading to systematic underestimation of their frequency. Numts and nupts, once inserted, can become further fragmented through subsequent insertion of mobile elements or other recombinational events that disrupt the continuity of the inserted sequence relative to the genuine organelle DNA copy. Because numts and nupts are typically identified through sequence comparison tools such as BLAST, disruption of insertions into smaller fragments can lead to systematic overestimation of numt and nupt frequencies. Accurate identification of numts and nupts is important, however, both for better understanding of their role during evolution, and for monitoring their increasingly evident role in human disease. Human populations are polymorphic for 141 numt loci, five numts are causal to genetic disease, and cancer genomic studies are revealing an abundance of numts associated with tumor progression. Here, we report investigation of salient parameters involved in obtaining accurate estimates of numt and nupt numbers in genome sequence data. Numts and nupts from 44 sequenced eukaryotic genomes reveal lineage-specific differences in the number, relative age and frequency of insertional events as well as lineage-specific dynamics of their postinsertional fragmentation. Our findings outline the main technical parameters influencing accurate identification and frequency estimation of numts in genomic studies pertinent to both evolution and human health.

Evolutionary conservation of bacterial operons: does transcriptional connectivity matter?

In the literature, it has been frequently suggested that the connectivity of a protein, i.e., the number of proteins with which it interacts, is inversely correlated with the rate of evolution. We attempted to extrapolate from proteins to operons by testing the hypothesis that operons with high transcriptional connectivity, i.e., operons that are controlled through interactions with many transcription factors, are evolutionarily more conserved at the structure and sequence levels than low-connectivity operons. With Escherichia coli used as reference, two structural- and two sequence-conservation measures were determined for 82 groups of homologous operons from 30 completely-sequenced bacterial genomes. In E. coli, large operons tend to be regulated by more transcription factors than either smaller operons or single genes. Large E. coli operons that are regulated by single transcription factors were found to be regulated by activators more frequently than by repressors. Levels of sequence conservation and structural conservation of operons were found to be independent of each other, i.e., structurally conserved operons may be divergent in sequence, and vice versa. Transcriptional connectivity was found to influence neither sequence nor structural conservation of operons. Although this finding seems to contradict the situation in genes, a critical review of the literature indicates that although gene connectivity is frequently touted as a factor in determining rates of evolution, only a very small fraction of the variability in degrees of evolutionary conservation is explainable by this factor.

Failure to Recover Major Events of Gene Flux in Real Biological Data Due to Method Misapplication

Abstract In prokaryotes, known mechanisms of lateral gene transfer (transformation, transduction, conjugation, and gene transfer agents) generate new combinations of genes among chromosomes during evolution. In eukaryotes, whose host lineage is descended from archaea, lateral gene transfer from organelles to the nucleus occurs at endosymbiotic events. Recent genome analyses studying gene distributions have uncovered evidence for sporadic, discontinuous events of gene transfer from bacteria to archaea during evolution. Other studies have used traditional models designed to investigate gene family size evolution (Count) to support claims that gene transfer to archaea was continuous during evolution, rather than involving occasional periodic mass gene influx events. Here, we show that the methodology used in analyses favoring continuous gene transfers to archaea was misapplied in other studies and does not recover known events of single simultaneous origin for many genes followed by differential loss in real data: plastid genomes. Using the same software and the same settings, we reanalyzed presence/absence pattern data for proteins encoded in plastid genomes and for eukaryotic protein families acquired from plastids. Contrary to expectations under a plastid origin model, we found that the methodology employed inferred that gene acquisitions occurred uniformly across the plant tree. Sometimes as many as nine different acquisitions by plastid DNA were inferred for the same protein family. That is, the methodology that recovered gradual and continuous lateral gene transfer among lineages for archaea obtains the same result for plastids, even though it is known that massive gains followed by gradual differential loss is the true evolutionary process that generated plastid gene distribution data. Our findings caution against the use of models designed to study gene family size evolution for investigating gene transfer processes, especially when transfers involving more than one gene per event are possible.

Count does not recover major events of gene flux in real biological data

In prokaryotes, known mechanisms of lateral gene transfer (transformation, transduction, conjugation and gene transfer agents) generate new combinations of genes among chromosomes during evolution. In eukaryotes, whose host lineage is descended from archaea, lateral gene transfer from organelles to the nucleus occurs at endosymbiotic events. Recent genome analyses studying gene distributions have uncovered evidence for sporadic, discontinuous events of gene transfer from bacteria to archaea during evolution. Other studies have used traditional birth-and-death phylogenetic models to investigate prokaryote genome evolution to claim that gene transfer to archaea was continuous during evolution, rather than involving occasional periodic mass gene influx events. Here we test the ability of Count, a birth-and-death based program, to recover known events of mass acquisition and differential loss using plastid genomes and eukaryotic protein families that were acquired from plastids. Count showed a strong bias towards reconstructed histories having gene acquisitions distributed uniformly across the tree. Sometimes as many as nine different acquisitions by plastid DNA were inferred for the same protein family. That is, Count recovered gradual and continuous lateral gene transfer among lineages, even when massive gains followed by gradual differential loss is the true evolutionary process that generated the gene distribution data.