Skip to content
naulon

No retraction, correction or concern found when last checked on Oct 29, 2026.

nQuire: a statistical framework for ploidy estimation using next generation sequencing

Clemens L. Weiß, Marina Pais, Liliana M. Cano and 2 more

BMC Bioinformatics · 2018

doi:10.1101/143537

Abstract

Background. Intraspecific variation in ploidy occurs in a wide range of species including pathogenic and nonpathogenic eukaryotes such as yeasts and oomycetes. Ploidy can be inferred indirectly - without measuring DNA content - from experiments using next-generation sequencing (NGS). We present nQuire, a statistical framework that distinguishes between diploids, triploids and tetraploids using NGS. The command-line tool models the distribution of base frequencies at variable sites using a Gaussian Mixture Model, and uses maximum likelihood to select the most plausible ploidy model. nQuire handles large genomes at high coverage efficiently and uses standard input file formats.

Results. We demonstrate the utility of nQuire analyzing individual samples of the pathogenic oomycete Phytophthora infestans and the Baker’s yeast Saccharomyces cerevisiae. Using these organisms we show the dependence between reliability of the ploidy assignment and sequencing depth. Additionally, we employ normalized maximized log- likelihoods generated by nQuire to ascertain ploidy level in a population of samples with ploidy heterogeneity. Using these normalized values we cluster samples in three dimensions using multivariate Gaussian mixtures. The cluster assignments retrieved from a S. cerevisiae population recovered the true ploidy level in over 96% of samples. Finally, we show that nQuire can be used regionally to identify chromosomal aneuploidies.

Conclusions. nQuire provides a statistical framework to study organisms with intraspecific variation in ploidy. nQuire is likely to be useful in epidemiological studies of pathogens, artificial selection experiments, and for historical or ancient samples where intact nuclei are not preserved. It is implemented as a stand-alone Linux command line tool in the C programming language and is available at https://github.com/clwgg/nQuireunder the MIT license.

Electronic supplementary material. The online version of this article (10.1186/s12859-018-2128-z) contains supplementary material, which is available to authorized users.

References 17

  1. 1.

    Selmecki AM, Maruvka YE, Richmond PA, Guillet M, Shoresh N, Sorenson AL, De S, Kishony R, Michor F, Dowell R, Pellman D. Polyploidy can drive rapid adaptation in yeast. Nature. 2015; 519(7543):349–52. 10.1038/nature14187.

    Not in our library10.1038/nature14187

  2. 2.

    Venkataram S, Dunn B, Li Y, Agarwala A, Chang J, Ebel E. R, Geiler-Samerotte K, Hérissant L, Blundell JR, Levy SF, Fisher DS, Sherlock G, Petrov DA. Development of a comprehensive Genotype-to-Fitness map of Adaptation-Driving mutations in yeast. Cell. 2016; 166(6):1585–9622. 10.1016/j.cell.2016.08.002.

    Not in our library10.1016/j.cell.2016.08.002

  3. 3.

    Comai L. The advantages and disadvantages of being polyploid. Nat Rev Genet. 2005; 6(11):836–46. 10.1038/nrg1711.

    Not in our library10.1038/nrg1711

  4. 4.
  5. 5.
  6. 6.

    Li Y, Shen H, Zhou Q, Qian K, van der Lee T, Huang S. Changing ploidy as a strategy: The irish potato famine pathogen shifts ploidy in relation to its sexuality. 2017; 30(1):45–52. 10.1094/MPMI-08-16-0156-R.

    Not in our library10.1094/mpmi-08-16-0156-r

  7. 7.

    Otto SP, Whitton J. Polyploid incidence and evolution. Annu Rev Genet. 2000; 34:401–37. 10.1146/annurev.genet.34.1.401.

    Not in our library10.1146/annurev.genet.34.1.401

  8. 8.

    Dirihan S, Terho P, Helander M, Saikkonen K. Efficient analysis of ploidy levels in plant evolutionary ecology. Caryologia. 2013; 66(3):251–6. 10.1080/00087114.2013.849414.

    Not in our library

  9. 9.

    Margarido GRA, Heckerman D. ConPADE: genome assembly ploidy estimation from next-generation sequencing data. PLoS Comput Biol. 2015; 11(4):1004229. 10.1371/journal.pcbi.1004229.

    Not in our library10.1371/journal.pcbi.1004229

  10. 10.

    Augusto Corrêa Dos Santos R, Goldman GH, Riaño-Pachón DM. ploidyNGS: visually exploring ploidy with next generation sequencing data. Bioinformatics. 2017; 33(16):2575–6. 10.1093/bioinformatics/btx204.

    Not in our library10.1093/bioinformatics/btx204

  11. 11.

    Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J, Banks E, Garimella KV, Altshuler D, Gabriel S, DePristo MA. From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013; 43:11–10133. 10.1002/0471250953.bi1110s43.

    Not in our library10.1002/0471250953.bi1110s43

  12. 12.

    Gompert Z, Mock KE. Detection of individual ploidy levels with genotyping-by-sequencing (GBS) analysis. Mol Ecol Resour. 2017. 10.1111/1755-0998.12657.

    Not in our library10.1111/1755-0998.12657

  13. 13.

    Scrucca L, Fop M, Murphy TB, Raftery AE. mclust 5: Clustering, classification and density estimation using gaussian finite mixture models. R J. 2016;8(1):289-317.

    Not in our library

  14. 14.

    R Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2008. http://www.R-project.org.

    Not in our library

  15. 15.

    Meyer M, Kircher M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc. 2010; 2010(6):5448. 10.1101/pdb.prot5448.

    Not in our library10.1101/pdb.prot5448

  16. 16.
  17. 17.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup. The sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25(16):2078–9. 10.1093/bioinformatics/btp352.

    Not in our library10.1093/bioinformatics/btp352

Cited by papers we hold 2

Reads

Read by people on this page

None yet.

Checked or read by AI through naulon

None yet.

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.

Are you an author?

Is this your paper?

Claim it to add your note and be paid when agents cite it. Free.