Gaps and complex structurally variant loci in phased genome assemblies

• Published in: Genome research Vol. 33; no. 4; pp. 496 - 510
• Main Authors: Porubsky, David, Vollger, Mitchell R., Harvey, William T., Rozanski, Allison N., Ebert, Peter, Hickey, Glenn, Hasenfeld, Patrick, Sanders, Ashley D., Stober, Catherine, Korbel, Jan O., Paten, Benedict, Marschall, Tobias, Eichler, Evan E.
• Format: Journal Article
• Language: English
• Published: United States Cold Spring Harbor Laboratory Press 01.04.2023
• Subjects: Chromosomes; DNA, Satellite - genetics; Genomes; Haplotypes; Humans; Polymorphism, Genetic; Satellite DNA; Segmental Duplications, Genomic; Sequence Analysis, DNA
• ISSN: 1088-9051; 1549-5469; 1549-5469
• DOI: 10.1101/gr.277334.122

There has been tremendous progress in phased genome assembly production by combining long-read data with parental information or linked-read data. Nevertheless, a typical phased genome assembly generated by trio-hifiasm still generates more than 140 gaps. We perform a detailed analysis of gaps, assembly breaks, and misorientations from 182 haploid assemblies obtained from a diversity panel of 77 unique human samples. Although trio-based approaches using HiFi are the current gold standard, chromosome-wide phasing accuracy is comparable when using Strand-seq instead of parental data. Importantly, the majority of assembly gaps cluster near the largest and most identical repeats (including segmental duplications [35.4%], satellite DNA [22.3%], or regions enriched in GA/AT-rich DNA [27.4%]). Consequently, 1513 protein-coding genes overlap assembly gaps in at least one haplotype, and 231 are recurrently disrupted or missing from five or more haplotypes. Furthermore, we estimate that 6–7 Mbp of DNA are misorientated per haplotype irrespective of whether trio-free or trio-based approaches are used. Of these misorientations, 81% correspond to bona fide large inversion polymorphisms in the human species, most of which are flanked by large segmental duplications. We also identify large-scale alignment discontinuities consistent with 11.9 Mbp of deletions and 161.4 Mbp of insertions per haploid genome. Although 99% of this variation corresponds to satellite DNA, we identify 230 regions of euchromatic DNA with frequent expansions and contractions, nearly half of which overlap with 197 protein-coding genes. Such variable and incompletely assembled regions are important targets for future algorithmic development and pangenome representation.