High-quality chromosome-scale genomes facilitate effective identification of large structural variations in hot and sweet peppers

Abstract

Pepper (Capsicum annuum) is an important vegetable crop that has been subjected to intensive breeding, resulting in limited genetic diversity, especially for sweet peppers. Previous studies have reported pepper draft genome assemblies using short read sequencing, but their capture of the extent of large structural variants (SVs), such as presence–absence variants (PAVs), inversions, and copy-number variants (CNVs) in the complex pepper genome falls short. In this study, we sequenced the genomes of representative sweet and hot pepper accessions by long-read and/or linked-read methods and advanced scaffolding technologies. First, we developed a high-quality reference genome for the sweet pepper cultivar ‘Dempsey’ and then used the reference genome to identify SVs in 11 other pepper accessions and constructed a graph-based pan-genome for pepper. We annotated an average of 42 972 gene families in each pepper accession, defining a set of 19 662 core and 23 115 non-core gene families. The new pepper pan-genome includes informative variants, 222 159 PAVs, 12 322 CNVs, and 16 032 inversions. Pan-genome analysis revealed PAVs associated with important agricultural traits, including potyvirus resistance, fruit color, pungency, and pepper fruit orientation. Comparatively, a large number of genes are affected by PAVs, which is positively correlated with the high frequency of transposable elements (TEs), indicating TEs play a key role in shaping the genomic landscape of peppers. The datasets presented herein provide a powerful new genomic resource for genetic analysis and genome-assisted breeding for pepper improvement.

References

[1] Hill TA, Chunthawodtiporn J, Ashrafi H et al. Regions underlying population structure and the genomics of organ size determination in Capsicum annuum . Plant Genome. 2017; 10: 1-14.

[2] Jeong HB, Kang MY, Jung A et al. Single-molecule real-time sequencing reveals diverse allelic variations in carotenoid biosynthetic genes in pepper (capsicum spp.) . Plant Biotechnol J. 2019; 17: 1081-93.

[3] Tao Y, Jordan DR, Mace ES . Crop genomics goes beyond a single reference genome. Trends Plant Sci. 2019; 24: 1072-4.

[4] Tao Y, Zhao X, Mace E et al. Exploring and exploiting pan-genomics for crop improvement. Mol Plant. 2019; 12: 156-69.

[5] Tranchant-Dubreuil C, Rouard M, Sabot F . Plant pangenome: impacts on phenotypes and evolution. Annu Plant Rev online. 2019; 2: 1-25.

[6] Jayakodi M, Padmarasu S, Haberer G et al. The barley pan-genome reveals the hidden legacy of mutation breeding. Nature. 2020; 588: 284-9.

[7] Tettelin H, Masignani V, Cieslewicz MJ et al. Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial ‘pan-genome’. Proc Natl Acad Sci U S A. 2005; 102: 13950-5.

[8] Golicz AA, Batley J, Edwards D . Towards plant pangenomics. Plant Biotechnol J. 2016; 14: 1099-105.

[9] Golicz AA, Bayer PE, Barker GC et al. The pangenome of an agronomically important crop plant Brassica oleracea . Nat Commun. 2016; 7: 1-8.

[10] Bayer PE, Golicz AA, Scheben A et al. Plant pan-genomes are the new reference. Nat Plants. 2020; 6: 914-20.

[11] Li YH, Zhou G, Ma J et al. De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits. Nat Biotechnol. 2014; 32: 1045-52.

[12] Gao L, Gonda I, Sun H et al. The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor. Nat Genet. 2019; 51: 1044- 51.

[13] Song JM, Guan Z, Hu J et al. Eight high-quality genomes reveal pan-genome architecture and ecotype differentiation of Brassica napus . Nat Plants. 2020; 6: 34-45.

[14] Hirsch CN, Foerster JM, Johnson JM et al. Insights into the maize pan-genome and pan-transcriptome. Plant Cell. 2014; 26: 121-35.

[15] Gordon SP, Contreras-Moreira B, Woods DP et al. Extensive gene content variation in the Brachypodium distachyon pan-genome correlates with population structure. Nat Commun. 2017; 8: 2184.

[16] Sun C, Hu Z, Zheng T et al. RPAN: Rice pan-genome browser for ~3000 rice genomes. Nucleic Acids Res. 2017; 45: 597-605.

[17] Zhao Q, Feng Q, Lu H et al. Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice. Nat Genet. 2018; 50: 278-84.

[18] Zhou Y, Chebotarov D, Kudrna D et al. A platinum standard pan-genome resource that represents the population structure of Asian rice. Sci Data. 2020; 7: 1-11.

[19] Walkowiak S, Gao L, Monat C et al. Multiple wheat genomes reveal global variation in modern breeding. Nature. 2020; 588: 277-83.

[20] Alonso-Blanco C, Andrade J, Becker C et al. 1,135 genomes reveal the global pattern of polymorphism in Arabidopsis thaliana . Cell. 2016; 166: 481-91.

[21] Pickersgill B . Cytogenetics and Evolution of Capsicum L. In: Tsuchiya T, Gupta PKBT-D in PG and B, eds. Chromosome Engineering in Plants.Amsterdam: Elsevier, 1991, 139-60.

[22] Paran I, Van Der Knaap E . Genetic and molecular regulation of fruit and plant domestication traits in tomato and pepper. J Exp Bot. 2007; 58: 3841-52.

[23] Lee HY, Ro NY, Jeong HJ et al. Genetic diversity and population structure analysis to construct a core collection from a large capsicum germplasm . BMC Genet. 2016; 17: 1-13.

[24] Simko I, Jia M, Venkatesh J et al. Genomics and marker-assisted improvement of vegetable crops. CRC Crit Rev Plant Sci. 2021; 40: 303-65.

[25] Kim S, Park M, Yeom SI et al. Genome sequence of the hot pepper provides insights into the evolution of pungency in capsicum species. Nat Genet. 2014; 46: 270-8.

[26] Qin C, Yu C, Shen Y et al. Whole-genome sequencing of cultivated and wild peppers provides insights into capsicum domestication and specialization. Proc Natl Acad Sci U S A. 2014; 111: 5135-40.

[27] Kim S, Park J, Yeom SI et al. New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication. Genome Biol. 2017; 18: 1-11.

[28] Hulse-Kemp AM, Maheshwari S, Stoffel K et al. Reference quality assembly of the 3.5-Gb genome of Capsicum annuum from a single linked-read library. Hortic Res. 2018; 5: 4.

[29] Acquadro A, Barchi L, Portis E et al. Whole genome resequencing of four Italian sweet pepper landraces provides insights on sequence variation in genes of agronomic value. Sci Rep. 2020; 10: 1-16.

[30] Han K, Jeong HJ, Yang HB et al. An ultra-high-density bin map facilitates high-throughput QTL mapping of horticultural traits in pepper (Capsicum annuum) . DNA Res. 2016; 23: 81-91.

[31] Jo YD, Ha Y, Lee JH et al. Fine mapping of restorer-of-fertility in pepper (Capsicum annuum L.) identified a candidate gene encoding a pentatricopeptide repeat (PPR)-containing protein. Theor Appl Genet. 2016; 129: 2003-17.

[32] Changkwian A, Venkatesh J, Lee JH et al. Physical localization of the root-knot nematode (Meloidogyne incognita) resistance locus Me7 in pepper (Capsicum annuum) . Front Plant Sci. 2019; 10: 886.

[33] Ou L, Li D, Lv J et al. Pan-genome of cultivated pepper (capsicum) and its use in gene presence-absence variation analyses . New Phytol. 2018; 220: 360-3.

[34] Hill T, Ashrafi H, Chin-Wo SR et al. Ultra-high density, transcript-based genetic maps of pepper define recombination in the genome and synteny among related species. G3 (Bethesda). 2015; 5: 2341-55.

[35] Hulse-Kemp AM, Ashrafi H, Plieske J et al. A HapMap leads to a Capsicum annuum SNP infinium array: a new tool for pepper breeding. Hortic Res. 2016; 3: 1-10.

[36] Solomon AM, Kim TG, Han K et al. Fine mapping and candidate gene identification for the CapUp locus controlling fruit orientation in pepper (capsicum spp.) . Front Plant Sci. 2021; 12: 1-12.

[37] Liao Y, Wang J, Zhu Z et al. The 3D architecture of the pepper genome and its relationship to function and evolution. Nat Commun. 2022; 13: 3479.

[38] Garrison E, Sirén J, Novak AM et al. Variation graph toolkit improves read mapping by representing genetic variation in the reference. Nat Biotechnol. 2018; 36: 875-9.

[39] Li H, Wang S, Chai S et al. Graph-based pan-genome reveals structural and sequence variations related to agronomic traits and domestication in cucumber. Nat Commun. 2022; 13: 1-14.

[40] Stewart C, Kang BC, Liu K et al. The Pun1 gene for pungency in pepper encodes a putative acyltransferase. Plant J. 2005; 42: 675-88.

[41] Ruffel S, Gallois JL, Moury B et al. Simultaneous mutations in translation initiation factors elF4E and elF(iso)4E are required to prevent pepper veinal mottle virus infection of pepper . J Gen Virol. 2006; 87: 2089-98.

[42] Elitzur T, Nahum H, Borovsky Y et al. Co-ordinated regulation of flowering time, plant architecture and growth by FASCICULATE: the pepper orthologue of SELF PRUNING . J Exp Bot. 2009; 60: 869-80.

[43] Tripodi P, Rabanus-Wallace MT, Barchi L et al. Global range expansion history of pepper (capsicum spp.) revealed by over 10,000 genebank accessions . Proc Natl Acad Sci U S A. 2021; 118: 1-9.

[44] Meyers LA, Levin DA . On the abundance of polyploids in flowering plants. Evolution (N Y). 2006; 60: 1198.

[45] Wellenreuther M, Mérot C, Berdan E et al. Going beyond SNPs: the role of structural genomic variants in adaptive evolution and species diversification. Mol Ecol. 2019; 28: 1203-9.

[46] Khan AW, Garg V, Roorkiwal M et al. Super-pangenome by integrating the wild side of a species for accelerated crop improvement. Trends Plant Sci. 2020; 25: 148-58.

[47] Michael TP, Van Buren R . Building near-complete plant genomes. Curr Opin Plant Biol. 2020; 54: 26-33.

[48] Sun Y, Shang L, Zhu QH et al. Twenty years of plant genome sequencing: achievements and challenges. Trends Plant Sci. 2022; 27: 391-401.

[49] Hurgobin B, Golicz AA, Bayer PE et al. Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus . Plant Biotechnol J. 2018; 16: 1265-74.

[50] Schatz MC, Maron LG, Stein JC et al. Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica . Genome Biol. 2014; 15: 506.

[51] Imaishi H, Matsuo S, Swai E et al. CYP78A1 preferentially expressed in developing inflorescences of Zea mays encoded a cytochrome P450-dependent lauric acid 12-monooxygenase. Biosci Biotechnol Biochem. 2000; 64: 1696-701.

[52] Neupane S, Ma Q, Mathew FM et al. Evolutionary divergence of TNL disease-resistant proteins in soybean (Glycine max) and common bean (Phaseolus vulgaris). Biochem Genet. 2018; 56: 397-422.

[53] Liu L, Tong H, Xiao Y et al. Activation of big Grain1 significantly improves grain size by regulating auxin transport in rice. Proc Natl Acad Sci U S A. 2015; 112: 11102-7.

[54] Zhang H, Zhao X, Ding X et al. Preparation of megabase-size DNA from plant nuclei. Plant J. 1995; 7: 175-84.

[55] Chin CS, Peluso P, Sedlazeck FJ et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat Methods. 2016; 13: 1050-4.

[56] Tang H, Zhang X, Miao C et al. ALLMAPS: robust scaffold ordering based on multiple maps. Genome Biol. 2015; 16: 1-15.

[57] Lu F, Romay MC, Glaubitz JC et al. High-resolution genetic mapping of maize pan-genome sequence anchors. Nat Commun. 2015; 6: 6914.

[58] Alonge M, Soyk S, Ramakrishnan S et al. RaGOO: fast and accurate reference-guided scaffolding of draft genomes. Genome Biol. 2019; 20: 224.

[59] Chernomor O, Von Haeseler A, Minh BQ . Terrace aware data structure for phylogenomic inference from supermatrices. Syst Biol. 2016; 65: 997-1008.

[60] Alexander DH, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64.

[61] Emms DM, Kelly S . OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019; 20: 238.

[62] Chen C, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194-202.

[63] Tian T, Liu Y, Yan H et al. AgriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017; 45: W122-9.

[64] Wang Y, Tang H, Debarry JD et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49-14.

[65] Carver TJ, Rutherford KM, Berriman M et al. ACT: the Artemis comparison tool. Bioinformatics. 2005; 21: 3422-3.

[66] Wang J, Zhang Z . GAPIT version 3: boosting power and accuracy for genomic association and prediction. Genomics Proteomics Bioinformatics. 2021; 19: 629-40.