Identical sequences found in distant genomes reveal frequent horizontal transfer across the bacterial domain

  1. Michael Sheinman  Is a corresponding author
  2. Ksenia Arkhipova
  3. Peter F Arndt
  4. Bas Dutilh
  5. Rutger Hermsen  Is a corresponding author
  6. Florian Massip  Is a corresponding author
  1. Netherlands Cancer Institute, Netherlands
  2. Utrecht University, Netherlands
  3. Max Planck Institute for Molecular Genetics, Germany
  4. Berlin Institute for Medical Systems Biology, Max Delbrueck Center for Molecular Medicine,, Germany

Abstract

Horizontal Gene Transfer (HGT) is an essential force in microbial evolution. Despite detailed studies on a variety of systems, a global picture of HGT in the microbial world is still missing. Here, we exploit that HGT creates long identical DNA sequences in the genomes of distant species, which can be found efficiently using alignment-free methods. Our pairwise analysis of 93 481 bacterial genomes identified 138 273 HGT events. We developed a model to explain their statistical properties as well as estimate the transfer rate between pairs of taxa. This reveals that long-distance HGT is frequent: our results indicate that HGT between species from different phyla has occurred in at least 8% of the species. Finally, our results confirm that the function of sequences strongly impacts their transfer rate, which varies by more than 3 orders of magnitude between different functional categories. Overall, we provide a comprehensive view of HGT, illuminating a fundamental process driving bacterial evolution.

Data availability

Results of the analysis are provided as supplementary files

The following previously published data sets were used

Article and author information

Author details

  1. Michael Sheinman

    Division of Molecular Carcinogenesis, Netherlands Cancer Institute, Amsterdam, Netherlands
    For correspondence
    mishashe@gmail.com
    Competing interests
    The authors declare that no competing interests exist.
  2. Ksenia Arkhipova

    Theoretical Biology and Bioinformatics, Utrecht University, Utrecht, Netherlands
    Competing interests
    The authors declare that no competing interests exist.
  3. Peter F Arndt

    Department of Computational Molecular Biology, Max Planck Institute for Molecular Genetics, Berlin, Germany
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1762-9836
  4. Bas Dutilh

    Theoretical Biology and Bioinformatics, Utrecht University, Utrecht, Netherlands
    Competing interests
    The authors declare that no competing interests exist.
  5. Rutger Hermsen

    Science faculty, Biology Department, Utrecht University, Utrecht, Netherlands
    For correspondence
    r.hermsen@uu.nl
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4633-4877
  6. Florian Massip

    Evolutionary and Cancer Genomics, Berlin Institute for Medical Systems Biology, Max Delbrueck Center for Molecular Medicine,, Berlin, Germany
    For correspondence
    florian.massip@mdc-berlin.de
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5855-0935

Funding

Netherlands Organisation for Scientific Research (Vidi grant 864.14.004)

  • Ksenia Arkhipova
  • Bas Dutilh

H2020 European Research Council (Consolidator Grant 865694: DiversiPHI)

  • Bas Dutilh

Fondation pour la Recherche Médicale (SPE201803005264)

  • Florian Massip

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Richard A Neher, University of Basel, Switzerland

Publication history

  1. Received: September 3, 2020
  2. Accepted: June 13, 2021
  3. Accepted Manuscript published: June 14, 2021 (version 1)
  4. Version of Record published: July 9, 2021 (version 2)

Copyright

© 2021, Sheinman et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 2,589
    Page views
  • 318
    Downloads
  • 2
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Michael Sheinman
  2. Ksenia Arkhipova
  3. Peter F Arndt
  4. Bas Dutilh
  5. Rutger Hermsen
  6. Florian Massip
(2021)
Identical sequences found in distant genomes reveal frequent horizontal transfer across the bacterial domain
eLife 10:e62719.
https://doi.org/10.7554/eLife.62719

Further reading

    1. Computational and Systems Biology
    Felix Proulx-Giraldeau, Jan M Skotheim, Paul François
    Research Article

    Cell size is controlled to be within a specific range to support physiological function. To control their size, cells use diverse mechanisms ranging from ‘sizers’, in which differences in cell size are compensated for in a single cell division cycle, to ‘adders’, in which a constant amount of cell growth occurs in each cell cycle. This diversity raises the question why a particular cell would implement one rather than another mechanism? To address this question, we performed a series of simulations evolving cell size control networks. The size control mechanism that evolved was influenced by both cell cycle structure and specific selection pressures. Moreover, evolved networks recapitulated known size control properties of naturally occurring networks. If the mechanism is based on a G1 size control and an S/G2/M timer, as found for budding yeast and some human cells, adders likely evolve. But, if the G1 phase is significantly longer than the S/G2/M phase, as is often the case in mammalian cells in vivo, sizers become more likely. Sizers also evolve when the cell cycle structure is inverted so that G1 is a timer, while S/G2/M performs size control, as is the case for the fission yeast S. pombe. For some size control networks, cell size consistently decreases in each cycle until a burst of cell cycle inhibitor drives an extended G1 phase much like the cell division cycle of the green algae Chlamydomonas. That these size control networks evolved such self-organized criticality shows how the evolution of complex systems can drive the emergence of critical processes.

    1. Computational and Systems Biology
    2. Neuroscience
    Kiri Choi, Won Kyu Kim, Changbong Hyeon
    Research Article

    The projection neurons (PNs), reconstructed from electron microscope (EM) images of the Drosophila olfactory system, offer a detailed view of neuronal anatomy, providing glimpses into information flow in the brain. About 150 uPNs constituting 58 glomeruli in the antennal lobe (AL) are bundled together in the axonal extension, routing the olfactory signal received at AL to mushroom body (MB) calyx and lateral horn (LH). Here we quantify the neuronal organization in terms of the inter-PN distances and examine its relationship with the odor types sensed by Drosophila. The homotypic uPNs that constitute glomeruli are tightly bundled and stereotyped in position throughout the neuropils, even though the glomerular PN organization in AL is no longer sustained in the higher brain center. Instead, odor-type dependent clusters consisting of multiple homotypes innervate the MB calyx and LH. Pheromone-encoding and hygro/thermo-sensing homotypes are spatially segregated in MB calyx, whereas two distinct clusters of food-related homotypes are found in LH in addition to the segregation of pheromone-encoding and hygro/thermo-sensing homotypes. We find that there are statistically significant associations between the spatial organization among a group of homotypic uPNs and certain stereotyped olfactory responses. Additionally, the signals from some of the tightly bundled homotypes converge to a specific group of lateral horn neurons (LHNs), which indicates that homotype (or odor type) specific integration of signals occurs at the synaptic interface between PNs and LHNs. Our findings suggest that before neural computation in the inner brain, some of the olfactory information are already encoded in the spatial organization of uPNs, illuminating that a certain degree of labeled-line strategy is at work in the Drosophila olfactory system.