HiRise Scaffolding Report


August 03, 2022

Theobroma cacao

Lyndel Meinhardt

USDAARS Sustainable Perennial Crops Lab


Contents

Overview

Assembly Total Length (bp) N50 L50 N90 L90
Input Assembly 452,075,119 31,634,083 6 16,245,998 13
Dovetail HiRise Assembly 452,076,199 29,395,232 7 13,674,979 16

Contiguity Metrics

Input Assembly Dovetail HiRise Assembly
Largest scaffold 50,529,124 45,398,999
Number of scaffolds 473 480
Number of scaffolds > 1kbp 473 480
Number of gaps 0 11
Number of N's per 100 kbp 0.00 0.24

BUSCO

Assembly Complete BUSCOs (C) Complete and single-copy BUSCOs (S) Complete and duplicated BUSCOs (D) Fragmented BUSCOs (F) Missing BUSCOs (M) Total BUSCO groups searched
Input Assembly 251 (98.43%) 238 (93.33%) 13 3 1 255
Dovetail HiRise Assembly 251 (98.43%) 238 (93.33%) 13 3 1 255

Pair Size Distribution

HiRise Scaffolding Information

Number of joins made by HiRise 11
Number of breaks made to input assembly by HiRise 4
Read-pairs 35,450,557

Materials and Methods

Dovetail Omni-C Library Preparation and Sequencing

For each Dovetail Omni-C library, chromatin was fixed in place with formaldehyde in the nucleus and then extracted. Fixed chromatin was digested with DNAse I, chromatin ends were repaired and ligated to a biotinylated bridge adapter followed by proximity ligation of adapter containing ends. After proximity ligation, crosslinks were reversed and the DNA purified. Purified DNA was treated to remove biotin that was not internal to ligated fragments. Sequencing libraries were generated using NEBNext Ultra enzymes and Illumina-compatible adapters. Biotin-containing fragments were isolated using streptavidin beads before PCR enrichment of each library. The library was sequenced on an Illumina HiSeqX platform to produce a approximately 30x sequence coverage. Then HiRise used MQ>50 reads for scaffolding (see "read-pair" above for figures).

Scaffolding the Assembly with HiRise

The input de novo assembly and Dovetail OmniC library reads were used as input data for HiRise, a software pipeline designed specifically for using proximity ligation data to scaffold genome assemblies (Putnam et al, 2016). Dovetail OmniC library sequences were aligned to the draft input assembly using bwa (https://github.com/lh3/bwa). The separations of Dovetail OmniC read pairs mapped within draft scaffolds were analyzed by HiRise to produce a likelihood model for genomic distance between read pairs, and the model was used to identify and break putative misjoins, to score prospective joins, and make joins above a threshold.


References

1. Putnam NH, O'Connell BL, Stites JC, Rice BJ, Blanchette M, Calef R, Troll CJ, Fields A, Hartley PD, Sugnet CW, Haussler D, Rokhsar DS, Green RE. Genome Research. 2016 Mar;26(3):342-50.