Pan-plastome approach empowers the assessment of genetic variation in cultivated Capsicum species

Abstract

Pepper species (Capsicum spp.) are widely used as food, spice, decoration, and medicine. Despite the recent old-world culinary impact, more than 50 commercially recognized pod types have been recorded worldwide from three taxonomic complexes (A, B, and P). The current study aimed to apply a pan-plastome approach to resolve the plastomic boundaries among those complexes and identify effective loci for the taxonomical resolution and molecular identification of the studied species/varieties. High-resolution pan-plastomes of five species and two varieties were assembled and compared from 321 accessions. Phyloplastomic and network analyses clarified the taxonomic position of the studied species/varieties and revealed a pronounced number of accessions to be the rare and endemic species, C. galapagoense, that were mistakenly labeled as C. annuum var. glabriusculum among others. Similarly, some NCBI-deposited plastomes were clustered differently from their labels. The rpl23-trnI intergenic spacer contained a 44 bp tandem repeat that, in addition to other InDels, was capable of discriminating the investigated Capsicum species/varieties. The rps16-trnQ/rbcL-accD/ycf3-trnS gene set was determined to be sufficiently polymorphic to retrieve the complete phyloplastomic signal among the studied Capsicum spp. The pan-plastome approach was shown to be useful in resolving the taxonomical complexes, settling the incomplete lineage sorting conflict and developing a molecular marker set for Capsicum spp. identification.

References

[1] Wang, H., Li, X. & Song, J. Vegetable genetic resources in China. Hortic. Plant J. 4, 83-88 (2018).

[2] Qin, C. 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).

[3] Kumar, S., Kumar, R. & Singh, J. in Handbook of herbs and spices 299-312 (Elsevier, Cambridge, England, 2006).

[4] Liu, W. Y., Kang, W. H. & Kang, B. C. in Genetics, Genomics and Breeding of Peppers and Eggplants Genetics, Genomics and Breeding of Crop Plants (eds B. C. Kang & C. Kole) Ch. 1, 1-13 (Taylor & Francis, Clemson, South Carolina, USA, 2013).

[5] Eshbaugh, W. H. in Peppers: botany, production and uses (ed Vincent, M. Russo) 14-28 (CAB International Publishing, New York, USA, 2012).

[6] Pickersgill, B. The genus Capsicum: a multidisciplinary approach to the taxonomy of cultivated and wild plants. Biol. Zent. 107, 381-389 (1988).

[7] Onus, A. N. & Pickersgill, B. Unilateral incompatibility in Capsicum (Solanaceae): occurrence and taxonomic distribution. Ann. Bot. 94, 289-295 (2004).

[8] Walsh, Brian M. & Hoot, Sara B. Phylogenetic relationships of Capsicum (Solanaceae) using DNA sequences from two noncoding regions: the chloroplast atpB-rbcL spacer region and nuclear waxy introns. Int. J. Plant Sci. 162, 1409-1418 (2001).

[9] Moscone, E. A. et al. in VI International Solanaceae Conference: Genomics Meets Biodiversity. 745. 137-170 (Madison, Wisconsin, USA, 2007).

[10] Carrizo García, C. et al. Phylogenetic relationships, diversification and expansion of chili peppers (Capsicum, Solanaceae). Ann. Bot. 118, 35-51 (2016).

[11] D’Agostino, N. et al. The complete plastome sequences of eleven Capsicum genotypes: insights into DNA variation and molecular evolution. Genes 9, 503 (2018).

[12] Group, C. P. W. et al. A DNA barcode for land plants. Proc. Natl Acad. Sci. USA 106, 12794-12797 (2009).

[13] Parks, M., Cronn, R. & Liston, A. Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes. BMC Biol. 7, 84 (2009).

[14] Jo, Y. D. et al. Complete sequencing and comparative analyses of the pepper (Capsicum annuum L.) plastome revealed high frequency of tandem repeats and large insertion/deletions on pepper plastome. Plant Cell Rep. 30, 217-229 (2011).

[15] Raveendar, S. et al. The complete chloroplast genome sequence of Korean landrace “Subicho” pepper (Capsicum annuum var. annuum). Plant Breed. Biotechnol. 3, 88-94 (2015).

[16] Shin, M.-J. et al. in 7th World Summit on Plant Genomics conference. (Bangkok, Thailand, 2017).

[17] Ou, L. et al. Pan-genome of cultivated pepper (Capsicum) and its use in gene presence-absence variation analyses. New Phytol. 220, 360-363 (2018).

[18] Poczai, P. & Hyvönen, J. Discovery of novel plastid phenylalanine (trn F) pseudogenes defines a distinctive clade in Solanaceae. Springer Plus 2, 459 (2013).

[19] Raveendar, S. et al. The complete chloroplast genome of Capsicum annuum var. glabriusculum using Illumina sequencing. Molecules 20, 13080-13088 (2015).

[20] Kim, T. S. et al. Complete chloroplast genome sequence of Capsicum baccatum var. baccatum. Mol. Breed. 36, 110 (2016).

[21] Shim, D. et al. The complete chloroplast genome of Capsicum frutescens (Solanaceae). Appl. Plant Sci. 4, 1600002 (2016).

[22] Park, H. S., Lee, J., Lee, S. C., Yang, T. J. & Yoon, J. B. The complete chloroplast genome sequence of Capsicum chinense Jacq. (Solanaceae). Mitochondrial DNA Part B 1, 164-165 (2016).

[23] Raveendar, S. et al. Complete chloroplast genome sequencing and genetic relationship analysis of Capsicum chinense jacq.. Plant Breed. Biotechnol. 5, 261-268 (2017).

[24] Pereira-Dias, L., Vilanova, S., Fita, A., Prohens, J. & Rodríguez-Burruezo, A. Genetic diversity, population structure, and relationships in a collection of pepper (Capsicum spp.) landraces from the Spanish centre of diversity revealed by genotyping-by-sequencing (GBS). Hortic. Res. 6, 54 (2019).

[25] Zeng, F. C., Gao, C. W. & Gao, L. Z. The complete chloroplast genome sequence of American bird pepper (Capsicum annuum var. glabriusculum). Mitochondrial DNA Part A 27, 724-726 (2016).

[26] Hunziker, A. T. Genera Solanacearum: the genera of Solanaceae illustrated, arranged according to a new system. (Konigstein, Germany, Koeltz Scientific Books, 2001).

[27] Vaughn, J. N., Chaluvadi, S. R., Rangan, L. & Bennetzen, J. L. Whole plastome sequences from five ginger species facilitate marker development and define limits to barcode methodology. PloS ONE 9, e108581 (2014).

[28] Sancho, R. et al. Comparative plastome genomics and phylogenomics of Brachypodium: flowering time signatures, introgression and recombination in recently diverged ecotypes. New Phytol. 218, 1631-1644 (2018).

[29] Daniell, H. et al. Complete chloroplast genome sequences of Solanum bulbocastanum, Solanum lycopersicum and comparative analyses with other Solanaceae genomes. Theor. Appl. Genet. 112, 1503 (2006).

[30] Logacheva, M. D., Schelkunov, M. I., Shtratnikova, V. Y., Matveeva, M. V. & Penin, A. A. Comparative analysis of plastid genomes of non-photosynthetic Ericaceae and their photosynthetic relatives. Sci. Rep. 6, 30042 (2016).

[31] Jarret, R. L. DNA barcoding in a crop genebank: the Capsicum annuum species complex. Open Biol. J. 1, 35-42 (2008).

[32] Osawa, S., Ohama, T., Jukes, T. H., Watanabe, K. & Yokoyama, S. Evolution of the mitochondrial genetic code II. Reassignment of codon AUA from isoleucine to methionine. J. Mol. Evol. 29, 373-380 (1989).

[33] Turmel, M., Lopes dos Santos, A., Otis, C., Sergerie, R. & Lemieux, C. Tracing the evolution of the plastome and mitogenome in the Chloropicophyceae uncovered convergent tRNA gene losses and a variant plastid genetic code. Genome Biol. Evol. 11, 1275-1292 (2019).

[34] Muir, G. & Filatov, D. A selective sweep in the chloroplast DNA of dioecious Silene (Section Elisanthe). Genetics 177, 1239-1247 (2007).

[35] Scally, A. et al. Insights into hominid evolution from the gorilla genome sequence. Nature 483, 169 (2012).

[36] Edwards, S. V. Is a new and general theory of molecular systematics emerging?. Evol. Int. J. Org. Evol. 63, 1-19 (2009).

[37] Kearse, M. et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647-1649 (2012).

[38] Tillich, M. et al. GeSeq-versatile and accurate annotation of organelle genomes. Nucleic acids Res. 45, W6-W11 (2017).

[39] Lowe, T. M. & Chan, P. P. tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 44, W54-W57 (2016).

[40] Katoh, K. & Standley, D. M. MAFFT Multiple Sequence Alignment Software Version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772-780 (2013).

[41] Huelsenbeck, J. P. & Ronquist, F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17, 754-755 (2001).

[42] Masters, B. C., Fan, V. & Ross, H. A. Species Delimitation-a Geneious plugin for the exploration of species boundaries. Mol. Ecol. Resour. 11, 154-157 (2011).

[43] Leigh, J. W. & Bryant, D. POPART: full-feature software for haplotype network construction. Methods Ecol. Evol. 6, 1110-1116 (2015).

[44] Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2-approximately maximum-likelihood trees for large alignments. PLoS ONE 5, e9490 (2010).

[45] Soria-Carrasco, V., Talavera, G., Igea, J. & Castresana, J. The K tree score: quantification of differences in the relative branch length and topology of phylogenetic trees. Bioinformatics 23, 2954-2956 (2007).

[46] Martin, D. P., Murrell, B., Golden, M., Khoosal, A. & Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. https://doi.org/10.1093/ve/vev003 (2015).

[47] Benson, G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573-580 (1999).

[48] Lohse, M., Drechsel, O., Kahlau, S. & Bock, R. OrganellarGenomeDRAW-a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 41, W575-W581 (2013).

[49] Rozas, J. et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299-3302 (2017).

[50] Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639-1645 (2009).