Complex fitness landscape shapes variation in a hyperpolymorphic species

  1. Anastasia V Stolyarova  Is a corresponding author
  2. Tatiana V Neretina
  3. Elena A Zvyagina
  4. Anna V Fedotova
  5. Alexey Kondrashov
  6. Georgii A. Bazykin  Is a corresponding author
  1. Skolkovo Institute of Science and Technology, Russian Federation
  2. Lomonosov Moscow State University, Russian Federation
  3. University of Michigan-Ann Arbor, United States
  4. Russian Academy of Sciences, Russian Federation

Abstract

It is natural to assume that patterns of genetic variation in hyperpolymorphic species can reveal large-scale properties of the fitness landscape that are hard to detect by studying species with ordinary levels of genetic variation. Here, we study such patterns in a fungus Schizophyllum commune, the most polymorphic species known. Throughout the genome, short-range linkage disequilibrium (LD) caused by attraction of minor alleles is higher between pairs of nonsynonymous than of synonymous variants. This effect is especially pronounced for pairs of sites that are located within the same gene, especially if a large fraction of the gene is covered by haploblocks, genome segments where the gene pool consists of two highly divergent haplotypes, which is a signature of balancing selection. Haploblocks are usually shorter than 1000 nucleotides, and collectively cover about 10% of the S. commune genome. LD tends to be substantially higher for pairs of nonsynonymous variants encoding amino acids that interact within the protein. There is a substantial correlation between LDs at the same pairs of nonsynonymous mutations in the USA and the Russian populations. These patterns indicate that selection in S. commune involves positive epistasis due to compensatory interactions between nonsynonymous alleles. When less polymorphic species are studied, analogous patterns can be detected only through interspecific comparisons.

Data availability

Whole-genome alignment of 55 genomes of S. commune is available at https://makarich.fbb.msu.ru/astolyarova/schizophyllum_data/. Genome sequence data are deposited at DDBJ/ENA/GenBank under accession numbers JAGVRL000000000-JAGVSI000000000, BioProject PRJNA720428. Sequencing data are deposited at SRA with accession numbers SRR14467839-SRR14467862.

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Anastasia V Stolyarova

    Skolkovo Institute of Science and Technology, Moscow, Russian Federation
    For correspondence
    anastasia.v.stolyarova@gmail.com
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6546-3052
  2. Tatiana V Neretina

    Biological Faculty, Lomonosov Moscow State University, Moscow, Russian Federation
    Competing interests
    The authors declare that no competing interests exist.
  3. Elena A Zvyagina

    Biological Faculty, Lomonosov Moscow State University, Moscow, Russian Federation
    Competing interests
    The authors declare that no competing interests exist.
  4. Anna V Fedotova

    Skolkovo Institute of Science and Technology, Moscow, Russian Federation
    Competing interests
    The authors declare that no competing interests exist.
  5. Alexey Kondrashov

    Department of Ecology and Evolutionary Biology, University of Michigan-Ann Arbor, Ann Arbor, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Georgii A. Bazykin

    Institute for Information Transmission Problems (Kharkevich Institute), Russian Academy of Sciences, Moscow, Russian Federation
    For correspondence
    gbazykin@iitp.ru
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2334-2751

Funding

No external funding was received for this work.

Copyright

© 2022, Stolyarova 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

  • 1,521
    views
  • 208
    downloads
  • 8
    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. Anastasia V Stolyarova
  2. Tatiana V Neretina
  3. Elena A Zvyagina
  4. Anna V Fedotova
  5. Alexey Kondrashov
  6. Georgii A. Bazykin
(2022)
Complex fitness landscape shapes variation in a hyperpolymorphic species
eLife 11:e76073.
https://doi.org/10.7554/eLife.76073

Share this article

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

Further reading

    1. Evolutionary Biology
    2. Neuroscience
    Gregor Belušič
    Insight

    The first complete 3D reconstruction of the compound eye of a minute wasp species sheds light on the nuts and bolts of size reduction.

    1. Evolutionary Biology
    2. Genetics and Genomics
    Julie N Chuong, Nadav Ben Nun ... David Gresham
    Research Article

    Copy number variants (CNVs) are an important source of genetic variation underlying rapid adaptation and genome evolution. Whereas point mutation rates vary with genomic location and local DNA features, the role of genome architecture in the formation and evolutionary dynamics of CNVs is poorly understood. Previously, we found the GAP1 gene in Saccharomyces cerevisiae undergoes frequent amplification and selection in glutamine-limitation. The gene is flanked by two long terminal repeats (LTRs) and proximate to an origin of DNA replication (autonomously replicating sequence, ARS), which likely promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we evolved engineered strains lacking either the adjacent LTRs, ARS, or all elements in glutamine-limited chemostats. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. Removal of local DNA elements significantly impacts the fitness effect of GAP1 CNVs and the rate of adaptation. In 177 CNV lineages, across all four strains, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, distal ones mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination can mediate gene amplification following de novo retrotransposon events. Our study reveals that template switching during DNA replication is a prevalent source of adaptive CNVs.