1. Evolutionary Biology
  2. Genetics and Genomics

High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox

  1. Rémi Allio  Is a corresponding author
  2. Marie-Ka Tilak
  3. Celine Scornavacca
  4. Nico L Avenant
  5. Andrew C Kitchener
  6. Erwan Corre
  7. Benoit Nabholz
  8. Frédéric Delsuc  Is a corresponding author
  1. CNRS, IRD, EPHE, Université de Montpellier, France
  2. National Museum and Centre for Environmental Management, University of the Free State, Bloemfontein, South Africa, France
  3. National Museums Scotland, United Kingdom
  4. CNRS, Sorbonne Université, FR2424, ABiMS, Station Biologique de Roscoff, 29680 Roscoff, France, France
  5. UMR 5554 ISEM (Université de Montpellier-CNRS-IRD-EPHE), France
  6. CNRS - Université de Montpellier, France
https://doi.org/10.7554/eLife.63167
Share this article
Cite this article
  1. Rémi Allio
  2. Marie-Ka Tilak
  3. Celine Scornavacca
  4. Nico L Avenant
  5. Andrew C Kitchener
  6. Erwan Corre
  7. Benoit Nabholz
  8. Frédéric Delsuc
(2021)
High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox
eLife 10:e63167.
https://doi.org/10.7554/eLife.63167
  2. Download BibTeX
  3. Download .RIS

Abstract

In a context of ongoing biodiversity erosion, obtaining genomic resources from wildlife is essential for conservation. The thousands of yearly mammalian roadkill provide a useful source material for genomic surveys. To illustrate the potential of this underexploited resource, we used roadkill samples to study the genomic diversity of the bat-eared fox (Otocyon megalotis) and the aardwolf (Proteles cristatus), both having subspecies with similar disjunct distributions in Eastern and Southern Africa. First, we obtained reference genomes with high contiguity and gene completeness by combining Nanopore long reads and Illumina short reads. Then, we showed that the two subspecies of aardwolf might warrant species status (P. cristatus and P. septentrionalis) by comparing their genome-wide genetic differentiation to pairs of well-defined species across Carnivora with a new Genetic Differentiation index (GDi) based on only a few resequenced individuals. Finally, we obtained a genome-scale Carnivora phylogeny including the new aardwolf species.

Data availability

Genome assemblies, associated Illumina and Nanopore sequence reads, and mitogenomes have been submitted to the National Center for Biotechnology Information (NCBI) and will be available after publication under BioProject number PRJNA681015. The full analytical pipeline, phylogenetic datasets (mitogenomic and genomic), corresponding trees, and other supplementary materials are available from zenodo.org (DOI: 10.5281/zenodo.4479226).

The following data sets were generated
The following previously published data sets were used
    1. Liu S
    2. Lorenzen ED
    3. Fumagalli M
    4. Li B
    5. Harris K
    6. Xiong Z
    7. Zhou L
    8. Korneliussen TS
    9. Somel M
    10. Babbitt C
    11. et al.
    (2014) Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears.
    PB43 : SRR942203, SRR942290, SRR942298; PB28: SRR942211, SRR942287, SRR942295; Brown Bear: SRR935591, SRR935625, SRR935627.
    https://www.ncbi.nlm.nih.gov

Article and author information

Author details

  1. Rémi Allio

    Institut des Sciences de l Evolution de Montpellier (ISEM), CNRS, IRD, EPHE, Université de Montpellier, Montpellier, France
    For correspondence
    rem.allio@yahoo.fr
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3885-5410

  2. Marie-Ka Tilak

    Institut des Sciences de l Evolution de Montpellier (ISEM), CNRS, IRD, EPHE, Université de Montpellier, Montpellier, France
    Competing interests
    The authors declare that no competing interests exist.

  3. Celine Scornavacca

    Institut des Sciences de l Evolution de Montpellier (ISEM), CNRS, IRD, EPHE, Université de Montpellier, France, CNRS, IRD, EPHE, Université de Montpellier, Montpellier, France
    Competing interests
    The authors declare that no competing interests exist.

  4. Nico L Avenant

    Department of Mammalogy, National Museum and Centre for Environmental Management, University of the Free State, Bloemfontein, South Africa, Montpellier, France
    Competing interests
    The authors declare that no competing interests exist.

  5. Andrew C Kitchener

    Department of Natural Sciences, National Museums Scotland, Edinburgh, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.

  6. Erwan Corre

    Informatics and Bioinformatics, CNRS, Sorbonne Université, FR2424, ABiMS, Station Biologique de Roscoff, 29680 Roscoff, France, Montpellier, France
    Competing interests
    The authors declare that no competing interests exist.

  7. Benoit Nabholz

    UMR 5554 ISEM (Université de Montpellier-CNRS-IRD-EPHE), Montpellier, France
    Competing interests
    The authors declare that no competing interests exist.

  8. Frédéric Delsuc

    Institut des Sciences de l'Evolution, CNRS - Université de Montpellier, Montpellier, France
    For correspondence
    frederic.delsuc@umontpellier.fr
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6501-6287

Funding

H2020 European Research Council (ERC‐2015‐CoG‐683257)

  • Frédéric Delsuc

Agence Nationale de la Recherche (ANR‐10‐LABX‐25‐01)

  • Rémi Allio
  • Marie-Ka Tilak
  • Celine Scornavacca
  • Benoit Nabholz
  • Frédéric Delsuc

Agence Nationale de la Recherche (ANR‐10‐LABX‐0004)

  • Rémi Allio
  • Marie-Ka Tilak
  • Celine Scornavacca
  • Benoit Nabholz
  • Frédéric Delsuc

Agence Nationale de la Recherche (ANR-11-INBS-0013)

  • Erwan Corre

National Research Foundation (86321)

  • Nico L Avenant

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

Copyright

© 2021, Allio 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,133
    views
  • 239
    downloads
  • 18
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

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. Rémi Allio
  2. Marie-Ka Tilak
  3. Celine Scornavacca
  4. Nico L Avenant
  5. Andrew C Kitchener
  6. Erwan Corre
  7. Benoit Nabholz
  8. Frédéric Delsuc
(2021)
High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox
eLife 10:e63167.
https://doi.org/10.7554/eLife.63167

Share this article

https://doi.org/10.7554/eLife.63167

Categories and tags

Further reading

    1. Chromosomes and Gene Expression
    2. Evolutionary Biology

    The emergence and evolution of gene expression in genome regions replete with regulatory motifs

    Timothy Fuqua, Yiqiao Sun, Andreas Wagner
    Research Article

    Gene regulation is essential for life and controlled by regulatory DNA. Mutations can modify the activity of regulatory DNA, and also create new regulatory DNA, a process called regulatory emergence. Non-regulatory and regulatory DNA contain motifs to which transcription factors may bind. In prokaryotes, gene expression requires a stretch of DNA called a promoter, which contains two motifs called –10 and –35 boxes. However, these motifs may occur in both promoters and non-promoter DNA in multiple copies. They have been implicated in some studies to improve promoter activity, and in others to repress it. Here, we ask whether the presence of such motifs in different genetic sequences influences promoter evolution and emergence. To understand whether and how promoter motifs influence promoter emergence and evolution, we start from 50 ‘promoter islands’, DNA sequences enriched with –10 and –35 boxes. We mutagenize these starting ‘parent’ sequences, and measure gene expression driven by 240,000 of the resulting mutants. We find that the probability that mutations create an active promoter varies more than 200-fold, and is not correlated with the number of promoter motifs. For parent sequences without promoter activity, mutations created over 1500 new –10 and –35 boxes at unique positions in the library, but only ~0.3% of these resulted in de-novo promoter activity. Only ~13% of all –10 and –35 boxes contribute to de-novo promoter activity. For parent sequences with promoter activity, mutations created new –10 and –35 boxes in 11 specific positions that partially overlap with preexisting ones to modulate expression. We also find that –10 and –35 boxes do not repress promoter activity. Overall, our work demonstrates how promoter motifs influence promoter emergence and evolution. It has implications for predicting and understanding regulatory evolution, de novo genes, and phenotypic evolution.

    1. Evolutionary Biology

    Estimating dispersal rates and locating genetic ancestors with genome-wide genealogies

    Matthew Osmond, Graham Coop
    Research Article Updated

    Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation-by-distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here, we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60 km2/generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.

    1. Evolutionary Biology

    Discarded sequencing reads uncover natural variation in pest resistance in Thlaspi arvense

    Dario Galanti, Jun Hee Jung ... Oliver Bossdorf
    Research Article

    Understanding the genomic basis of natural variation in plant pest resistance is an important goal in plant science, but it usually requires large and labor-intensive phenotyping experiments. Here, we explored the possibility that non-target reads from plant DNA sequencing can serve as phenotyping proxies for addressing such questions. We used data from a whole-genome and -epigenome sequencing study of 207 natural lines of field pennycress (Thlaspi arvense) that were grown in a common environment and spontaneously colonized by aphids, mildew, and other microbes. We found that the numbers of non-target reads assigned to the pest species differed between populations, had significant SNP-based heritability, and were associated with climate of origin and baseline glucosinolate contents. Specifically, pennycress lines from cold and thermally fluctuating habitats, presumably less favorable to aphids, showed higher aphid DNA load, i.e., decreased aphid resistance. Genome-wide association analyses identified genetic variants at known defense genes but also novel genomic regions associated with variation in aphid and mildew DNA load. Moreover, we found several differentially methylated regions associated with pathogen loads, in particular differential methylation at transposons and hypomethylation in the promoter of a gene involved in stomatal closure, likely induced by pathogens. Our study provides first insights into the defense mechanisms of Thlaspi arvense, a rising crop and model species, and demonstrates that non-target whole-genome sequencing reads, usually discarded, can be leveraged to estimate intensities of plant biotic interactions. With rapidly increasing numbers of large sequencing datasets worldwide, this approach should have broad application in fundamental and applied research.