Reference quality assembly of the 3.5-Gb genome of Capsicum annuum from a single linked-read library

Abstract

Linked-Read sequencing technology has recently been employed successfully for de novo assembly of human genomes, however, the utility of this technology for complex plant genomes is unproven. We evaluated the technology for this purpose by sequencing the 3.5-gigabase (Gb) diploid pepper (Capsicum annuum) genome with a single Linked-Read library. Plant genomes, including pepper, are characterized by long, highly similar repetitive sequences. Accordingly, significant effort is used to ensure that the sequenced plant is highly homozygous and the resulting assembly is a haploid consensus. With a phased assembly approach, we targeted a heterozygous F₁ derived from a wide cross to assess the ability to derive both haplotypes and characterize a pungency gene with a large insertion/deletion. The Supernova software generated a highly ordered, more contiguous sequence assembly than all currently available C. annuum reference genomes. Over 83% of the final assembly was anchored and oriented using four publicly available de novo linkage maps. A comparison of the annotation of conserved eukaryotic genes indicated the completeness of assembly. The validity of the phased assembly is further demonstrated with the complete recovery of both 2.5-Kb insertion/deletion haplotypes of the PUN1 locus in the F₁ sample that represents pungent and nonpungent peppers, as well as nearly full recovery of the BUSCO2 gene set within each of the two haplotypes. The most contiguous pepper genome assembly to date has been generated which demonstrates that Linked-Read library technology provides a tool to de novo assemble complex highly repetitive heterozygous plant genomes. This technology can provide an opportunity to cost-effectively develop high-quality genome assemblies for other complex plants and compare structural and gene differences through accurate haplotype reconstruction.

References

[1] Genome 10K Community of Scientists . Genome 10K: a proposal to obtain whole-genome sequence for 10,000 vertebrate species. J. Hered. 100, 659-674 (2009).

[2] i5K Consortium . The i5K Initiative: advancing arthropod genomics for knowledge, human health, agriculture, and the environment. J. Hered. 104, 595-600 (2013).

[3] Zhang, G., Rahbek, C., Graves, G. R., Lei, F., Jarvis, E. D. & Gilbert, M. T. Genomics: bird sequencing project takes off. Nature 522, 34 (2015).

[4] Kyrpides, N. C., Hugenholtz, P., Eisen, J. A., Woyke, T., Goker, M. & Parker, C. T. et al. Genomic encyclopedia of bacteria and archaea: sequencing a myriad of type strains. PLoS Biol. 12, e1001920 (2014).

[5] McCouch, S., Baute, G. J., Bradeen, J., Bramel, P., Bretting, P. K. & Buckler, E. et al. Agriculture: feeding the future. Nature 499, 23-24 (2013).

[6] Marouli, E., Graff, M., Medina-Gomez, C., Lo, K. S., Wood, A. R. & Kjaer, T. R. et al. Rare and low-frequency coding variants alter human adult height. Nature 542, 186-190 (2017).

[7] Yang, J., Bakshi, A., Zhu, Z., Hemani, G., Vinkhuyzen, A. A. E. & Lee, S. H. et al. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat. Genet. 47, 1114 (2015).

[8] Sanger, F., Nicklen, S. & Coulson, A. R. DNA sequencing with chain-terminating inhibitors. Proc Natl Acad. Sci. USA 74, 5463-5467 (1977).

[9] Bevan, M. W., Uauy, C., Wulff, B. B., Zhou, J., Krasileva, K. & Clark, M. D. Genomic innovation for crop improvement. Nature 543, 346-354 (2017).

[10] Stewart, C. Jr., Kang, B. C., Liu, K., Mazourek, M., Moore, S. L. & Yoo, E. Y. et al. The Pun1 gene for pungency in pepper encodes a putative acyltransferase. Plant J. 42, 675-688 (2005).

[11] Kim, S., Park, M., Yeom, S. I., Kim, Y. M., Lee, J. M. & Lee, H. A. et al. Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat. Genet. 46, 270-278 (2014).

[12] Qin, C., Yu, C., Shen, Y., Fang, X., Chen, L. & Min, J. et al. Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization. Proc. Natl Acad. Sci. USA 111, 5135-5140 (2014).

[13] Hill, T., Ashrafi, H., Chin-Wo, S. R., Stoffel, K., Truco, M. J. & Kozik, A. et al. Ultra-high density, transcript-based genetic maps of pepper define recombination in the genome and synteny among related species. G3 5, 2341-2355 (2015).

[14] Hulse-Kemp, A. M., Ashrafi, H., Plieske, J., Lemm, J., Stoffel, K. & Hill, T. et al. A HapMap leads to a Capsicum annuum SNP infinium array: a new tool for pepper breeding. Hortic. Res. 3, 16036 (2016).

[15] Han, K., Jeong, H. J., Yang, H. B., Kang, S. M., Kwon, J. K. & Kim, S. et al. An ultra-high-density bin map facilitates high-throughput QTL mapping of horticultural traits in pepper (Capsicum annuum). DNA Res. 23, 81-91 (2016).

[16] Weisenfeld, N. I., Kumar, V., Shah, P., Church, D. M. & Jaffe, D. B. Direct determination of diploid genome sequences. Genome Res. 27, 757-767 (2017).

[17] Stoffel, K., van Leeuwen, H., Kozik, A., Caldwell, D., Ashrafi, H. & Cui, X. et al. Development and application of a 6.5 million feature Affymetrix Genechip(R) for massively parallel discovery of single position polymorphisms in lettuce (Lactuca spp.). BMC Genomics 13, 185 (2012).

[18] Zheng, G. X., Lau, B. T., Schnall-Levin, M., Jarosz, M., Bell, J. M. & Hindson, C. M. et al. Haplotyping germline and cancer genomes with high-throughput linked-read sequencing. Nat. Biotechnol. 34, 303-311 (2016).

[19] Wu, T. D. & Watanabe, C. K. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 21, 1859-1875 (2005).

[20] Tang, H., Zhang, X., Miao, C., Zhang, J., Ming, R. & Schnable, J. C. et al. ALLMAPS: robust scaffold ordering based on multiple maps. Genome Biol. 16, 3 (2015).

[21] Wu, F., Eannetta, N. T., Xu, Y., Durrett, R., Mazourek, M. & Jahn, M. M. et al. A COSII genetic map of the pepper genome provides a detailed picture of synteny with tomato and new insights into recent chromosome evolution in the genus Capsicum. Theor. Appl. Genet. 118, 1279-1293 (2009).

[22] Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072-1075 (2013).

[23] Kurtz, S., Phillippy, A., Delcher, A. L., Smoot, M., Shumway, M. & Antonescu, C. et al. Versatile and open software for comparing large genomes. Genome Biol. 5, R12 (2004).

[24] Nattestad, M. & Schatz, M. C. Assemblytics: a web analytics tool for the detection of variants from an assembly. Bioinformatics 32, 3021-3023 (2016).

[25] Simao, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[26] Chin, C. S., Peluso, P., Sedlazeck, F. J., Nattestad, M., Concepcion, G. T. & Clum, A. et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat. Methods 13, 1050-1054 (2016).

[27] Sato, S., Tabata, S., Hirakawa, H., Asamizu, E., Shirasawa, K. & Isobe, S. et al. The tomato genome sequence provides insights into fleshy fruit evolution. Nature 485, 635-641 (2012).

[28] Xu, X., Pan, S. K., Cheng, S. F., Zhang, B., Mu, D. S. & Ni, P. X. et al. Genome sequence and analysis of the tuber crop potato. Nature 475, 189-U194 (2011).

[29] Sierro, N., Battey, J. N., Ouadi, S., Bakaher, N., Bovet, L. & Willig, A. et al. The tobacco genome sequence and its comparison with those of tomato and potato. Nat. Commun. 5, 3833 (2014).

[30] Hirakawa, H., Shirasawa, K., Miyatake, K., Nunome, T., Negoro, S. & Ohyama, A. et al. Draft genome sequence of eggplant (Solanum melongena L.): the representative solanum species indigenous to the old world. DNA Res. 21, 649-660 (2014).

[31] Bombarely, A., Moser, M., Amrad, A., Bao, M., Bapaume, L. & Barry, C. S. et al. Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida. Nat. Plants 2, 16074 (2016).

[32] Iorizzo, M., Ellison, S., Senalik, D., Zeng, P., Satapoomin, P. & Huang, J. et al. A high-quality carrot genome assembly provides new insights into carotenoid accumulation and asterid genome evolution. Nat. Genet. 48, 657-666 (2016).

[33] Pickersgill, B. A. Cytogenetics and evolution of Capsicum L.. Chromosome engineering in plants: genetics, breeding, evolution, part B.Elsevier, Amsterdam. 1991, pp 139-160.

[34] Rinaldi, R., Van Deynze, A., Portis, E., Rotino, G. L., Toppino, L. & Hill, T. et al. New insights on eggplant/tomato/pepper synteny and identification of eggplant and pepper orthologous QTL. Front Plant Sci. 7, 1031 (2016).

[35] Edgar, R. C. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5, 113 (2004).

[36] Moscone, E. A., Lambrou, M., Hunziker, A. T. & Ehrorfer, F. Giemsa C-banded karyotypes in Capsicum (Solanaceae). Plant Syst. Evol. 186, 213-229 (1993).