{"id":608715,"date":"2023-02-16T20:49:21","date_gmt":"2023-02-17T02:49:21","guid":{"rendered":"https:\/\/news.sellorbuyhomefast.com\/index.php\/2023\/02\/16\/telomere-to-telomere-assembly-of-diploid-chromosomes-with-verkko\/"},"modified":"2023-02-16T20:49:21","modified_gmt":"2023-02-17T02:49:21","slug":"telomere-to-telomere-assembly-of-diploid-chromosomes-with-verkko","status":"publish","type":"post","link":"https:\/\/newsycanuse.com\/index.php\/2023\/02\/16\/telomere-to-telomere-assembly-of-diploid-chromosomes-with-verkko\/","title":{"rendered":"Telomere-to-telomere assembly of diploid chromosomes with Verkko"},"content":{"rendered":"

Science & Nature <\/p>\n

\n
\n

Data availability<\/h2>\n

No new data were generated for this study. All assemblies generated in this paper are archived at Zenodo78<\/a><\/sup> and we have provided convenient links to download both data and assemblies79<\/a><\/sup>. The data are also hosted in public databases: A. thaliana<\/i> PRJCA005809, H. axyridis<\/i> PRJEB45202<\/a>, CHM13 PRJNA559484<\/a>, HG002 SAMN03283347 and the HPRC AWS bucket80<\/a><\/sup>.<\/p>\n<\/div>\n

\n

Code availability<\/h2>\n

Verkko code is available from GitHub81<\/a><\/sup> and all code used for the paper is archived at Zenodo78<\/a><\/sup>.<\/p>\n<\/div>\n

\n

References<\/h2>\n
\n
    \n
  1. \n

    Logsdon, G. A., Vollger, M. R. & Eichler, E. E. Long-read human genome sequencing and its applications. Nat. Rev. Genet.<\/i> 21<\/b>, 597\u2013614 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  2. \n

    Wenger, A. M. et al. Accurate circular consensus long-read sequencing improves variant detection and assembly of a human genome. Nat. Biotechnol.<\/i> 37<\/b>, 1155\u20131162 (2019).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  3. \n

    Jain, M. et al. Nanopore sequencing and assembly of a human genome with ultra-long reads. Nat. Biotechnol.<\/i> 36<\/b>, 338\u2013345 (2018).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  4. \n

    Shafin, K. et al. Nanopore sequencing and the Shasta toolkit enable efficient de novo assembly of eleven human genomes. Nat. Biotechnol.<\/i> 38<\/b>, 1044\u20131053 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  5. \n

    Nagarajan, N. & Pop, M. Sequencing and genome assembly using next-generation technologies. Methods Mol. Biol.<\/i> 673<\/b>, 1\u201317 (2010).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  6. \n

    Koren, S. & Phillippy, A. M. One chromosome, one contig: complete microbial genomes from long-read sequencing and assembly. Curr. Opin. Microbiol.<\/i> 23C<\/b>, 110\u2013120 (2014).<\/p>\n


    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  7. \n

    Kolmogorov, M., Yuan, J., Lin, Y. & Pevzner, P. A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol.<\/i> 37<\/b>, 540\u2013546 (2019).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  8. \n

    Nurk, S. et al. HiCanu: accurate assembly of segmental duplications, satellites, and allelic variants from high-fidelity long reads. Genome Res.<\/i> https:\/\/doi.org\/10.1101\/gr.263566.120<\/a> (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  9. \n

    Cheng, H., Concepcion, G. T., Feng, X., Zhang, H. & Li, H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat. Methods<\/i> 18<\/b>, 170\u2013175 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  10. \n

    Rhie, A. et al. Towards complete and error-free genome assemblies of all vertebrate species. Nature<\/i> 592<\/b>, 737\u2013746 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  11. \n

    Jarvis, E. D. et al. Semi-automated assembly of high-quality diploid human reference genomes. Nature<\/i> 611<\/b>, 519\u2013531 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  12. \n

    Miga, K. H. et al. Telomere-to-telomere assembly of a complete human X chromosome. Nature<\/i> 585<\/b>, 79\u201384 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  13. \n

    Logsdon, G. A. et al. The structure, function and evolution of a complete human chromosome 8. Nature<\/i> 593<\/b>, 101\u2013107 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  14. \n

    Nurk, S. et al. The complete sequence of a human genome. Science<\/i> 376<\/b>, 44\u201353 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  15. \n

    Bankevich, A., Bzikadze, A. V., Kolmogorov, M., Antipov, D. & Pevzner, P. A. Multiplex de Bruijn graphs enable genome assembly from long, high-fidelity reads. Nat. Biotechnol.<\/i> 40<\/b>, 1075\u20131081 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  16. \n

    Schwartz, D. C. et al. Ordered restriction maps of Saccharomyces cerevisiae<\/i> chromosomes constructed by optical mapping. Science<\/i> 262<\/b>, 110\u2013114 (1993).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  17. \n

    Ghareghani, M. et al. Strand-seq enables reliable separation of long reads by chromosome via expectation maximization. Bioinformatics<\/i> 34<\/b>, i115\u2013i123 (2018).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  18. \n

    Porubsky, D. et al. Fully phased human genome assembly without parental data using single-cell strand sequencing and long reads. Nat. Biotechnol.<\/i> 39<\/b>, 302\u2013308 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  19. \n

    O\u2019Neill, K. et al. Assembling draft genomes using contiBAIT. Bioinformatics<\/i> 33<\/b>, 2737\u20132739 (2017).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  20. \n

    Burton, J. N. et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat. Biotechnol.<\/i> 31<\/b>, 1119\u20131125 (2013).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  21. \n

    Dudchenko, Olga et al. De novo assembly of the Aedes aegypti<\/i> genome using Hi-C yields chromosome-length scaffolds. Science<\/i> 356<\/b>, 92\u201395 (2017).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  22. \n

    Ghurye, J. et al. Integrating Hi-C links with assembly graphs for chromosome-scale assembly. PLoS Comput. Biol.<\/i> 15<\/b>, e1007273 (2019).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  23. \n

    Howe, K. et al. Significantly improving the quality of genome assemblies through curation. GigaScience<\/i> 10<\/b>, giaa153 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  24. \n

    Edge, P., Bafna, V. & Bansal, V. HapCUT2: robust and accurate haplotype assembly for diverse sequencing technologies. Genome Res.<\/i> 27<\/b>, 801\u2013812 (2017).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  25. \n

    Antipov, D., Korobeynikov, A., McLean, J. S. & Pevzner, P. A. hybridSPAdes: an algorithm for hybrid assembly of short and long reads. Bioinformatics<\/i> 32<\/b>, 1009\u20131015 (2016).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  26. \n

    Wick, R. R., Judd, L. M., Gorrie, C. L. & Holt, K. E. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol.<\/i> 13<\/b>, e1005595 (2017).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  27. \n

    Di Genova, A., Buena-Atienza, E., Ossowski, S. & Sagot, M.-F. Efficient hybrid de novo assembly of human genomes with WENGAN. Nat. Biotechnol.<\/i> 39<\/b>, 422\u2013430 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  28. \n

    Koren, S. et al. De novo assembly of haplotype-resolved genomes with trio binning. Nat. Biotechnol.<\/i> 36<\/b>, 1174\u20131182 (2018).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  29. \n

    Falconer, E. et al. DNA template strand sequencing of single-cells maps genomic rearrangements at high resolution. Nat. Methods<\/i> 9<\/b>, 1107\u20131112 (2012).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  30. \n

    Sanders, A. D., Falconer, E., Hills, M., Spierings, D. C. J. & Lansdorp, P. M. Single-cell template strand sequencing by Strand-seq enables the characterization of individual homologs. Nat. Protoc.<\/i> 12<\/b>, 1151\u20131176 (2017).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  31. \n

    Miller, J. R. et al. Aggressive assembly of pyrosequencing reads with mates. Bioinformatics<\/i> 24<\/b>, 2818\u20132824 (2008).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  32. \n

    Idury, R. M. & Waterman, M. S. A new algorithm for DNA sequence assembly. J. Comput. Biol.<\/i> 2<\/b>, 291\u2013306 (1995).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  33. \n

    Wang, B. et al. High-quality Arabidopsis thaliana genome assembly with Nanopore and HiFi long reads. Genomics Proteomics Bioinformatics<\/i> 20<\/b>, 4\u201313 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  34. \n

    Cheng, H. et al. Haplotype-resolved assembly of diploid genomes without parental data. Nat. Biotechnol.<\/i> 40<\/b>, 1332\u20131335 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  35. \n

    Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics<\/i> 29<\/b>, 1072\u20131075 (2013).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  36. \n

    Mikheenko, A., Bzikadze, A. V., Gurevich, A., Miga, K. H. & Pevzner, P. A. TandemTools: mapping long reads and assessing\/improving assembly quality in extra-long tandem repeats. Bioinformatics<\/i> 36<\/b>, i75\u2013i83 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  37. \n

    Mc Cartney, A. M. et al. Chasing perfection: validation and polishing strategies for telomere-to-telomere genome assemblies. Nat. Methods<\/i> 19<\/b>, 687\u2013695 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  38. \n

    Li, H. New strategies to improve minimap2 alignment accuracy. Bioinformatics<\/i> 37<\/b>, 4572\u20134574 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  39. \n

    Boyes, D. et al. The genome sequence of the harlequin ladybird, Harmonia axyridis<\/i> (Pallas, 1773). Wellcome Open Res.<\/i> 7<\/b>, 177 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  40. \n

    Garg, S. et al. Chromosome-scale, haplotype-resolved assembly of human genomes. Nat. Biotechnol.<\/i> 39<\/b>, 309\u2013312 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  41. \n

    Chin, C.-S. & Khalak, A. Human genome assembly in 100 minutes. Preprint at bioRxiv<\/i> https:\/\/doi.org\/10.1101\/705616<\/a> (2019).<\/p>\n<\/li>\n

  42. \n

    Zook, J. M. et al. Extensive sequencing of seven human genomes to characterize benchmark reference materials. Sci. Data<\/i> 3<\/b>, 1\u201326 (2016).<\/p>\n

    Article<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  43. \n

    Wang, T. et al. The Human Pangenome Project: a global resource to map genomic diversity. Nature<\/i> 604<\/b>, 437\u2013446 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  44. \n

    Baid, G. et al. DeepConsensus improves the accuracy of sequences with a gap-aware sequence transformer. Nat. Biotechnol.<\/i> https:\/\/doi.org\/10.1038\/s41587-022-01435-7<\/a> (2022).<\/p>\n<\/li>\n

  45. \n

    Rhie, A. et al. The complete sequence of a human Y chromosome. Preprint at bioRxiv<\/i> https:\/\/doi.org\/10.1101\/2022.12.01.518724<\/a> (2022).<\/p>\n<\/li>\n

  46. \n

    Ewing, B. & Green, P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res.<\/i> 8<\/b>, 186\u2013194 (1998).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  47. \n

    Guan, D. et al. Identifying and removing haplotypic duplication in primary genome assemblies. Bioinformatics<\/i> 36<\/b>, 2896\u20132898 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  48. \n

    Porubsky, D. et al. Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders. Cell<\/i> 185<\/b>, 1986\u20132005.e26 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  49. \n

    Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows\u2013Wheeler transform. Bioinformatics<\/i> 25<\/b>, 1754\u20131760 (2009).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  50. \n

    Li, H. et al. The Sequence Alignment\/Map format and SAMtools. Bioinformatics<\/i> 25<\/b>, 2078\u20132079 (2009).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  51. \n

    Porubsky, D. et al. breakpointR: an R\/Bioconductor package to localize strand state changes in Strand-seq data. Bioinformatics<\/i> 36<\/b>, 1260\u20131261 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  52. \n

    Mohajeri, K. et al. Interchromosomal core duplicons drive both evolutionary instability and disease susceptibility of the Chromosome 8p23.1 region. Genome Res.<\/i> 26<\/b>, 1453\u20131467 (2016).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  53. \n

    McNulty, S. M. & Sullivan, B. A. Alpha satellite DNA biology: finding function in the recesses of the genome. Chromosome Res.<\/i> 26<\/b>, 115\u2013138 (2018).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  54. \n

    Mahtani, M. M. & Willard, H. F. Pulsed-field gel analysis of \u03b1-satellite DNA at the human X chromosome centromere: high-frequency polymorphisms and array size estimate. Genomics<\/i> 7<\/b>, 607\u2013613 (1990).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  55. \n

    Wevrick, R. & Willard, H. F. Physical map of the centromeric region of human chromosome 7: relationship between two distinct alpha satellite arrays. Nucleic Acids Res.<\/i> 19<\/b>, 2295\u20132301 (1991).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  56. \n

    Waye, J. S. & Willard, H. F. Chromosome specificity of satellite DNAs: short- and long-range organization of a diverged dimeric subset of human alpha satellite from chromosome 3. Chromosoma<\/i> 97<\/b>, 475\u2013480 (1989).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  57. \n

    Waye, J. S. et al. Chromosome-specific alpha satellite DNA from human chromosome 1: hierarchical structure and genomic organization of a polymorphic domain spanning several hundred kilobase pairs of centromeric DNA. Genomics<\/i> 1<\/b>, 43\u201351 (1987).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  58. \n

    Willard, H. F. et al. Detection of restriction fragment length polymorphisms at the centromeres of human chromosomes by using chromosome-specific alpha satellite DNA probes: implications for development of centromere-based genetic linkage maps. Proc. Natl Acad. Sci. USA<\/i> 83<\/b>, 5611\u20135615 (1986).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  59. \n

    Wevrick, R. & Willard, H. F. Long-range organization of tandem arrays of alpha satellite DNA at the centromeres of human chromosomes: high-frequency array-length polymorphism and meiotic stability. Proc. Natl Acad. Sci. USA<\/i> 86<\/b>, 9394\u20139398 (1989).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  60. \n

    de Lima, L. G. et al. PCR amplicons identify widespread copy number variation in human centromeric arrays and instability in cancer. Cell Genomics<\/i> 1<\/b>, 100064 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  61. \n

    KeyGene. Maize B73 Oxford Nanopore duplex sequence data release. https:\/\/www.keygene.com\/news-events\/maize-b73-oxford-nanopore-duplex-sequence-data-release\/<\/a> (2022).<\/p>\n<\/li>\n

  62. \n

    Altemose, N. et al. Complete genomic and epigenetic maps of human centromeres. Science<\/i> https:\/\/doi.org\/10.1126\/science.abl4178<\/a> (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  63. \n

    Langley, S. A., Miga, K. H., Karpen, G. H. & Langley, C. H. Haplotypes spanning centromeric regions reveal persistence of large blocks of archaic DNA. eLife<\/i> 8<\/b>, e42989 (2019).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  64. \n

    Jain, C., Koren, S., Dilthey, A., Phillippy, A. M. & Aluru, S. A fast adaptive algorithm for computing whole-genome homology maps. Bioinformatics<\/i> 34<\/b>, i748\u2013i756 (2018).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  65. \n

    Quinlan, A. R. & Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics<\/i> 26<\/b>, 841\u2013842 (2010).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  66. \n

    Rhie, A., Walenz, B. P., Koren, S. & Phillippy, A. M. Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies. Genome Biol.<\/i> 21<\/b>, 245 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  67. \n

    Rautiainen, M. & Marschall, T. MBG: minimizer-based sparse de Bruijn graph construction. Bioinformatics<\/i> 37<\/b>, 2476\u20132478 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  68. \n

    Rautiainen, M. & Marschall, T. GraphAligner: rapid and versatile sequence-to-graph alignment. Genome Biol.<\/i> 21<\/b>, 253 (2020).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  69. \n

    Wick, R. R., Schultz, M. B., Zobel, J. & Holt, K. E. Bandage: interactive visualization of de novo genome assemblies. Bioinformatics<\/i> 31<\/b>, 3350\u20133352 (2015).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  70. \n

    Vollger, M. R., Kerpedjiev, P., Phillippy, A. M. & Eichler, E. E. StainedGlass: interactive visualization of massive tandem repeat structures with identity heatmaps. Bioinformatics<\/i> 38<\/b>, 2049\u20132051 (2022).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  71. \n

    Onodera, T., Sadakane, K. & Shibuya, T. in Algorithms in Bioinformatics<\/i> (eds Darling, A. & Stoye, J.) 338\u2013348 (Springer Berlin Heidelberg, 2013).<\/p>\n<\/li>\n

  72. \n

    Roberts, M., Hayes, W., Hunt, B. R., Mount, S. M. & Yorke, J. A. Reducing storage requirements for biological sequence comparison. Bioinformatics<\/i> 20<\/b>, 3363\u20133369 (2004).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  73. \n

    Edgar, R. Syncmers are more sensitive than minimizers for selecting conserved k<\/i>\u2011mers in biological sequences. PeerJ<\/i> 9<\/b>, e10805 (2021).<\/p>\n

    Article<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  74. \n

    Ferragina, P. & Manzini, G. Indexing compressed text. J. ACM<\/i> 52<\/b>, 552\u2013581 (2005).<\/p>\n

    Article<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  75. \n

    Myers, G. A fast bit-vector algorithm for approximate string matching based on dynamic programming. J. ACM<\/i> 46<\/b>, 395\u2013415 (1999).<\/p>\n

    Article<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  76. \n

    Li, H. Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences. Bioinformatics<\/i> 32<\/b>, 2103\u20132110 (2016).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  77. \n

    Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k<\/i>-mer weighting and repeat separation. Genome Res.<\/i> 27<\/b>, 722\u2013736 (2017).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  78. \n

    Koren, S. Verkko beta2 source and assemblies evaluated in manuscript. Zenodo https:\/\/doi.org\/10.5281\/zenodo.6618379<\/a> (2022).<\/p>\n<\/li>\n

  79. \n

    Koren, S. verkko publication readme. GitHub https:\/\/github.com\/marbl\/verkko\/blob\/master\/paper\/README.md<\/a> (2022).<\/p>\n<\/li>\n

  80. \n

    HPRC HG002 public data. Amazon S3 https:\/\/s3-us-west-2.amazonaws.com\/human-pangenomics\/index.html?prefix=<\/a> (2022).<\/p>\n<\/li>\n

  81. \n

    Koren, S. verkko repository. GitHub https:\/\/github.com\/marbl\/verkko\/<\/a> (2022).<\/p>\n<\/li>\n

  82. \n

    Robinson, J. T. et al. Integrative genomics viewer. Nat. Biotechnol.<\/i> 29<\/b>, 24\u201326 (2011).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  83. \n

    Vollger, M. R. et al. Long-read sequence and assembly of segmental duplications. Nat. Methods<\/i> 16<\/b>, 88\u201394 (2019).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  84. \n

    Smith George, P. Evolution of repeated DNA sequences by unequal crossover. Science<\/i> 191<\/b>, 528\u2013535 (1976).<\/p>\n

    Article<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  85. \n

    Alkan, C., Eichler, E. E., Bailey, J. A., Sahinalp, S. C. & T\u00fcz\u00fcn, E. The role of unequal crossover in alpha-satellite DNA evolution: a computational analysis. J. Comput. Biol.<\/i> 11<\/b>, 933\u2013944 (2004).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  86. \n

    Alkan, C., Bailey, J. A., Eichler, E. E., Sahinalp, S. C. & Tuzun, E. An algorithmic analysis of the role of unequal crossover in alpha-satellite DNA evolution. Genome Inform.<\/i> 13<\/b>, 93\u2013102 (2002).<\/p>\n

    CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n

  87. \n

    Schindelhauer, D. & Schwarz, T. Evidence for a fast, intrachromosomal conversion mechanism from mapping of nucleotide variants within a homogeneous \u03b1-satellite DNA array. Genome Res.<\/i> 12<\/b>, 1815\u20131826 (2002).<\/p>\n

    Article<\/a>\u00a0
    \n     CAS<\/a>\u00a0
    \n     PubMed<\/a>\u00a0
    \n     PubMed Central<\/a>\u00a0
    \n
    \n Google Scholar<\/a>\u00a0\n <\/p>\n<\/li>\n<\/ol>\n

    Download references<\/a><\/p>\n<\/div>\n<\/div>\n

    \n

    Acknowledgements<\/h2>\n

    This work was supported, in part, by the Intramural Research Program of the National Human Genome Research Institute, National Institutes of Health (M.R., S.N., A.R., B.P.W., A.M.P. and S.K.) as well as by grants from the US National Institutes of Health (NIH grant nos. HG010169 and HG002385 to E.E.E.) and the National Institute of General Medical Sciences (NIGMS grant no. 1F32GM134558 to G.A.L.). E.E.E. is an investigator of the Howard Hughes Medical Institute. This work utilized the computational resources of the NIH HPC Biowulf cluster (https:\/\/hpc.nih.gov<\/a>).<\/p>\n<\/div>\n

    \n

    Author information<\/h2>\n
    \n

    Author notes<\/span><\/p>\n

      \n
    1. \n

      Sergey Nurk<\/p>\n

      Present address: Oxford Nanopore Technologies, Oxford, UK<\/p>\n<\/li>\n<\/ol>\n

      Authors and Affiliations<\/h3>\n
        \n
      1. \n

        Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA<\/p>\n

        Mikko Rautiainen,\u00a0Sergey Nurk,\u00a0Brian P. Walenz,\u00a0Arang Rhie,\u00a0Adam M. Phillippy\u00a0&\u00a0Sergey Koren<\/p>\n<\/li>\n

      2. \n

        Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA<\/p>\n

        Glennis A. Logsdon,\u00a0David Porubsky\u00a0&\u00a0Evan E. Eichler<\/p>\n<\/li>\n

      3. \n

        Howard Hughes Medical Institute, University of Washington, Seattle, WA, USA<\/p>\n

        Evan E. Eichler<\/p>\n<\/li>\n<\/ol>\n

        Contributions<\/h3>\n

        M.R., S.N., B.P.W. and S.K. were responsible for the methods and software development. G.A.L., D.P., A.R. and S.K. were responsible for data analysis and validation. E.E.E. and A.M.P. provided resources. M.R., S.N., A.M.P. and S.K. wrote the first draft of the manuscript. M.R., S.N., G.A.L., D.P., A.M.P. and S.K. prepared the figures. M.R., S.N., B.P.W., A.M.P. and S.K. edited the manuscript with the assistance of all authors. E.E.E., A.M.P. and S.K. supervised the study. M.R., S.N., A.M.P. and S.K. conceptualized the study.<\/p>\n

        Corresponding authors<\/h3>\n

        Correspondence to
        \n Adam M. Phillippy<\/a> or Sergey Koren<\/a>.<\/p>\n<\/div>\n<\/div>\n

        \n

        Ethics declarations<\/h2>\n
        \n

        Competing interests<\/h3>\n

        E.E.E. is on the scientific advisory board of DNAnexus, Inc. S.K. has received travel funds to speak at events hosted by Oxford Nanopore Technologies. S.N. is an employee of Oxford Nanopore Technologies. The remaining authors declare no competing interests.<\/p>\n<\/p><\/div>\n<\/div>\n

        \n

        Peer review<\/h2>\n
        \n

        Peer review information<\/h3>\n

        Nature Biotechnology<\/i> thanks Rayan Chikhi, Anton Korobeynikov and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.<\/p>\n<\/p><\/div>\n<\/div>\n

        \n

        Additional information<\/h2>\n

        Publisher\u2019s note<\/b> Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.<\/p>\n<\/div>\n

        \n

        Extended data<\/h2>\n
        \n
        \n

        Extended Data Fig. 1 A. thaliana<\/i> chromosome unitigs in Verkko (left) vs published assembly chromosomes evaluated by VerityMap (right).<\/a><\/h3>\n

        From top to bottom, Chr1, Chr2, Chr3, Chr4, and Chr5. VerityMap compares the spacing of unique k-mers within the HiFi reads to the spacing observed in the assembly. Whenever there is a disagreement, the plot shows a spike at the discrepant location. The x<\/i>-axis indicates the coordinates along the assembly contig or scaffold while the y<\/i>-axis shows the fraction of disagreeing reads (0\u2013100%). A disagreement greater than 50% is likely not a heterozygous variant but a true error in the assembly. The BED file produced by VerityMap also indicates the size of the discrepancy, estimated from the difference in k-mer spacing between the reads and the assembly.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 2 Verkko CHM13 assembly sub-graphs.<\/a><\/h3>\n

        A<\/b>. The remaining unresolved regions in CHM13 chromosomes 5, 9 and 16, visualized using Bandage69<\/a><\/sup>, with the correct resolution marked in red paths. Left: Chr5 has a spurious edge causing a cycle, and three spurious low-coverage nodes which were not removed by bubble popping since they are a part of the cycle. Middle: Chr9 has a spurious edge. Right: Chr16 has two spurious edges, and one missing edge (dashed red curve). The spurious non-genomic edges are caused by noisy ONT alignments switching between highly similar repeats in the LA graph, while the missing edge is caused by low HiFi coverage. B<\/b>. rDNA cluster mixing in CHM13 chromosomes 13, 14, and 21, visualized using Bandage69<\/a><\/sup>. Each chromosome has a separate rDNA tangle. There are two cross-chromosomal connections by erroneous low coverage (<4x) nodes circled in red. For all three chromosomes, the remainder of the p and q arms are contained in the long unitigs shown.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 3 VerityMap discrepant reads plot for CHM13 HiFi and ONT unitigs assembled by Verkko (left) and CHM13 v1.114<\/a><\/sup> (right).<\/a><\/h3>\n

        A<\/b>. The assemblies for Chromosome 4. The Verkko assembly has no regions where a large fraction of reads are deviated even though QUAST marks an error at approximately 52\u2009Mb. This corresponds to a position in the reference with a large fraction of deviated reads and an estimated 19\u2009kb discrepancy. B<\/b>. same for Chromosome 17. There are no regions with a large fraction (>50%) of discrepant reads in the Verkko assembly despite QUAST reporting an error at approximately 25\u2009Mb on the reference. This corresponds to an approximately 3\u2009kb discrepancy identified by VerityMap in CHM13 v1.1.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 4 Merqury66<\/a><\/sup> haplotype blob plots.<\/a><\/h3>\n

        A<\/b>. HG002 downsampled Verkko B<\/b>. HG002 downsampled DeepConsensus HiFi Verkko and C<\/b>. HG002 full-coverage Verkko assemblies. The Hi-C phased assembly is on the left and the trio-phased assembly is on the right. Each contig\/scaffold is a circle on the plot, with the size scaled based on contig\/scaffold length. The x<\/i>-axis shows the number of maternal markers while the y<\/i>-axis shows the number of paternal markers. Contigs which lie along either the x-axis or y-axis show no haplotype errors and are consistently maternal or paternal. Contigs which mixed haplotypes would appear along the diagonal but are not observed in these plots, indicating an accurately phased assembly.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 5 IGV82<\/a><\/sup> views of a recently published HG002 diploid assembly of paternal Chromosome 10 11<\/a><\/sup> (top) and the Verkko full-coverage trio assembly of the same chromosome (bottom).<\/a><\/h3>\n

        The tracks show the maternal (red) and paternal (blue) markers. The centromere location is shown in gray. The published assembly has extensive switching within the centromere array, indicated by the presence of maternal markers and the absence of paternal markers. In contrast, the Verkko assembly centromere shows only paternal markers. The Verkko paternal centromere array is shorter but shows no signs of mis-assembly (Extended Data Fig. 8<\/a>) indicating the larger array in the published assembly is likely due to the incorrect insertion of maternal sequence. Overall, the Verkko assembly is more continuous, with 0 gaps vs 4, and a lower hamming error rate, 0.03%, versus 1.98% compared to the published assembly.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 6 Strand-seq validation of the full-coverage Verkko trio assembly and HPRC manually curated assembly11<\/a><\/sup>.<\/a><\/h3>\n

        The maternal haplotype is shown along the top row and the paternal along the bottom row. Leftmost:<\/b> alignment-based scaffold assignment to the maternal haplotype (top) and paternal haplotype (bottom) for the full-coverage Verkko assembly. Almost all chromosomes are a single color, indicating that Verkko scaffolds resolved most chromosomes end-to-end. The only exceptions are in the acrocentrics, where some of the scaffolds could not be assigned due to low mappability and maternal Chromosome 6 and paternal Chromosomes 5 which are each composed of two large scaffolds. Over 99.7% of the scaffold bases could be assigned to chromosomes. Middle:<\/b> the cluster assignment for the maternal haplotype (top) and paternal haplotype (bottom) based on Strand-seq data for the full-coverage Verkko assembly. Here, cluster ID is assigned to each 200\u2009kb window in a scaffold. In case of large scale chromosomal mis-joins, we expect to see multiple colors in a chromosome. The Verkko assembly is consistent with scaffolds all representing a single chromosome bin. Once again, >99.7% of the scaffold bases can be assigned using Strand-seq. Only 2 and 4\u2009Mb of sequence not scaffolded by Verkko could be assigned to the maternal and paternal haplotypes, respectively. Right:<\/b> The cluster assignment for the maternal haplotype (top) and the paternal haplotype (bottom) based on Strand-seq data for the HPRC manually curated assembly. Here, cluster ID is assigned to each 200\u2009kb window in a scaffold. In case of large scale chromosomal mis-joins, we expect to see multiple colors in a chromosome. A smaller fraction of contigs (and a slightly lower fraction of bases) was assigned than for the Verkko assembly, despite the combination of technologies and manual curation. This may be due to shorter contigs from unresolved repeats which are resolved through Verkko\u2019s ONT integration. There is also visible chromosome mixing within the acrocentric chromosomes unlike in the Verkko result.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 7 Strand-seq structural variant analysis for Verkko full-coverage assembly.<\/a><\/h3>\n

        The states assigned to each scaffold in the paternal (A<\/b>) and maternal (B<\/b>) for the full-coverage Verkko trio assembly. Strand-seq reads aligned to each assembly are genotype based on their directionality into three possible strand states. Crick-Crick (\u2018cc\u2019) state in which both homologs in Strand-seq data map in direct orientation and thus such regions are consistent with Strand-seq directional information. Watson-Watson (\u2018ww\u2019) state in which both homologs in Strand-seq data map in inverted orientation and are indicative of assembly misorientation or unresolved homozygous inversion. Lastly, there are a few (<1% of bases) Watson-Crick (\u2018wc\u2019) where there is a mixture of Watson and Crick reads and such regions are indicative of heterozygous inversions between haplotypes or low-mappability regions for short Strand-seq reads. C<\/b>. The size of the heterozygous inversion versus the count of inversions of that size in the maternal and paternal haplotypes of the full-coverage Verkko trio assembly. These regions have confident Strand-seq alignments and normal copy number so these regions indicate potential true heterozygous variation between the haplotypes. D<\/b>. Strand-seq alignments to the reference Chromosome Y before it was corrected (top) and full-coverage Verkko trio Chromosome Y assembly (bottom). Each plot shows Strand-seq directional read coverage reported as binned (bin size: 10,000, step size: 1,000) read counts represented as vertical bars above (teal; Crick read counts) and below (orange; Watson read counts) the midline. The top plot shows an inversion (dashed line) where directly oriented reads (Crick; teal) switch to inversely oriented reads (Watson, orange) and then back to directly oriented reads. The Verkko assembly in contrast is consistent with only Crick reads present in the same location (dashed line).<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 8 Full-coverage Verkko trio assemblies of chromosome 1 (a), 3 (b), 4 (c), 11 (d), 9 (e), 10 (f), 16 (g), and 18 (h) centromeric regions in the HG002 genome.<\/a><\/h3>\n

        Both maternal and paternal haplotypes are shown, with repeat element annotation generated by RepeatMasker (cite:1. Smit, A., Hubley, R. & Green, P. Smit, AFA, Hubley, R & Green, P. RepeatMasker Open-4.0. 2013\u20132015. RepeatMasker Open-4.0. http:\/\/www.repeatmasker.org<\/a> (2013)) shown on top, followed by PacBio HiFi coverage, ONT coverage, and StainedGlass70<\/a><\/sup> plots. As with the Chromosome 19 centromeres (Fig. 4<\/a>), the maternal and paternal haplotypes show large-scale structural variation, with alpha-satellite HOR arrays sizes varying by tens to hundreds of kb. Sites with discrepant HiFi mappings (low coverage or high coverage) are marked with an asterisk. There are few sites in the centromeres, and the artifacts are localized and often inconsistent between ONT and HiFi alignments, indicating the assembly is overall of high quality. To further validate assembly accuracy, we intersected centromere array locations with VerityMap errors and found that in all but four cases (two on the Chr1 paternal centromere, Chr9 paternal centromere, and Chr10 maternal centromere), the errors were short (\u22641\u2009kb) or lower frequency (\u226450% of the reads). VerityMap also identified one issue, with \u226550% of reads deviating in the Chr4 maternal centromere. However, this was not visible in the NucFreq 37<\/a>,83<\/a><\/sup> plots above, and the region only had a total of three mapped reads.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 9 Comparison of the HG002 maternal and paternal full-coverage Verkko trio assemblies for the centromeric regions of chromosomes 1 (a), 3 (b), 4 (c), 9 (d), 10 (e), 11 (f), 16 (g), 18 (h), and 19 (i) in the HG002 genome.<\/a><\/h3>\n

        The plots show the similarity between the two haplotypes, with the maternal haplotype on the y-axis and the paternal on the x-axis. The centromeric regions show varying \u0251-satellite HOR array sizes and sequence identity between the two haplotypes, consistent with earlier reports that indicate that centromeric HOR arrays often expand and contract due to their repetitive nature and their propensity for unequal crossing over84<\/a>,85<\/a>,86<\/a><\/sup> and gene conversion87<\/a><\/sup> events. For Chromosome 19, as in Fig. 4<\/a>, the tracks show the repeat annotations and read coverages. The triangles show the self-similarity within each haplotype for comparison.<\/p>\n<\/div>\n

        \n

        Extended Data Fig. 10 Examples of haplotype scaffolding by Rukki in the HG002 genome.<\/a><\/h3>\n

        The nodes are colored according to their haplotype assignments. Nodes with at least 100 total markers where 90% of the markers agree are colored: red for maternal, blue for paternal. Nodes with less than 100 markers are colored gray for unassigned. The haplotype paths are marked with solid curves with dotted curves for gaps. (A<\/b>) A well behaved genomic region consisting of phased heterozygous bubbles, homozygous nodes, and spurious nodes caused by sequencing errors. Where possible, Rukki connects the nodes attributed to the same haplotype across the homozygous regions, producing two phased unitigs without gaps. (B<\/b>) A tangle within one haplotype. Rukki scaffolds across the tangle (dotted line), reporting an estimated size of the tangled region. (C<\/b>) A gap in the paternal haplotype. Rukki uses haplotype assignments and the topology of the graph to scaffold across the gap (dotted line), and estimates the size of the gap based on the size of the paired haplotype.<\/p>\n<\/div>\n<\/div>\n<\/div>\n

        \n

        Supplementary information<\/h2>\n<\/div>\n