Frequent germplasm exchanges drive the high genetic diversity of Chinese-cultivated common apricot germplasm

Abstract

The genetic diversity of germplasm is critical for exploring genetic and phenotypic resources and has important implications for crop-breeding sustainability and improvement. However, little is known about the factors that shape and maintain genetic diversity. Here, we assembled a high-quality chromosome-level reference of the Chinese common apricot ‘Yinxiangbai’, and we resequenced 180 apricot accessions that cover four major ecogeographical groups in China and other accessions from occidental countries. We concluded that Chinese-cultivated common apricot germplasms possessed much higher genetic diversity than those cultivated in Western countries. We also detected seven migration events among different apricot groups, where 27% of the genome was identified as being introgressed. Remarkably, we demonstrated that these introgressed regions drove the current high level of germplasm diversity in Chinese-cultivated common apricots by introducing different genes related to distinct phenotypes from different cultivated groups. Our results highlight the consideration that introgressed regions may provide an important reservoir of genetic resources that can be used to sustain modern breeding programs.

References

[1] Flint-Garcia, S. A. Genetics and consequences of crop domestication. J. Agric. Food Chem. 61, 8267-8276 (2013).

[2] Migicovsky, Z., Warschefsky, E., Klein, L. L. & Miller, A. J. Using living germplasm collections to characterize, improve, and conserve woody perennials. Crop Sci. 59, 2365-2380 (2019).

[3] Powell, W. & Barsby, T. Germplasm diversity and genetics to drive plant breeding for Africa. Successful Agricultural Innovation in Emerging Economies: New Genetic Technologies for Global Food Production (eds Bennett, D. & Jennings, R.) (Cambridge University Press, 2013).

[4] Bourguiba, H. et al. Genetic structure of a worldwide germplasm collection of Prunus armeniaca L. Reveals three major diffusion routes for varieties coming from the species’ center of origin. Front. Plant Sci. 11, 638 (2020).

[5] Faust, M., Surányi, D. & Nyujtó, F. Origin and Dissemination of Apricot. Hortic. Rev. 22, 225-260 (1998).

[6] Zhang, J. Y. & Zhang, Z. China fruit-plant monographs, apricot flora. 17-46 (China Forestry Press, 2003).

[7] Badenes, M. L. & Byrne, D. H. Fruit Breeding, Handbook of Plant Breeding 8. 415 (Springer Science+Business Media, 2012).

[8] Zhang, Q. & Liu, W. Advances of the apricot resources collection, evaluation and germplasm enhancement. Acta Hortic. Sin. 45, 1642-1660 (2018). in Chinese.

[9] Zhang, Q.-P. et al. Genetic diversity and relationships of common apricot (Prunus armeniaca L.) in China based on simple sequence repeat (SSR) markers. Genet. Resour. Crop Ev. 61, 357-368 (2014).

[10] Jiang, F. et al. The apricot (Prunus armeniaca L.) genome elucidates Rosaceae evolution and beta-carotenoid synthesis. Hortic. Res. 6, 128 (2019).

[11] 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).

[12] Verde, I. et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nat. Genet. 45, 487-494 (2013).

[13] Verde, I. et al. The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity. BMC Genomics 18, 225 (2017).

[14] Zhang, Q. et al. The genome of Prunus mume. Nat. Commun. 3, 1318 (2012).

[15] Shirasawa, K. et al. The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding. DNA Res. 24, 499-508 (2017).

[16] Liu, W. et al. Kernel-using apricot resources and its utilization. Acta Hortic. 966, 189-191 (2012).

[17] Cao, K. et al. Comparative population genomics reveals the domestication history of the peach, Prunus persica, and human influences on perennial fruit crops. Genome Biol. 15, 415 (2014).

[18] Wu, J. et al. Diversification and independent domestication of Asian and European pears. Genome Biol. 19, 77 (2018).

[19] Liu, S. et al. The complex evolutionary history of apricots: Species divergence, gene flow and multiple domestication events. Mol. Ecol. 28, 5299-5314 (2019).

[20] Zhebentyayeva, T. N., Reighard, G. L., Gorina, V. M. & Abbott, A. G. Simple sequence repeat (SSR) analysis for assessment of genetic variability in apricot germplasm. Theor. Appl. Genet. 106, 435-444 (2003).

[21] He, T.-M. et al. Using SSR markers to determine the population genetic structure of wild apricot (Prunus armeniaca L.) in the Ily valley of West China. Genet. Resour. Crop Ev. 54, 563-572 (2007).

[22] Pedryc, A. et al. Genetic diversity of apricot revealed by a set of SSR markers from linkage group G1. Sci. Hortic.-Amst. 121, 19-26 (2009).

[23] Yuan, Z. et al. Population genetic structure in apricot (Armeniaca Mill.) revealed by fluorescent-AFLP markers. Acta Hortic. Sin. 35, 319-328 (2008). In Chinese.

[24] Kole, C. Genome Mapping and Molecular Breeding in Plants in Fruits and Nuts. (Springer-Verlag Berlin Heidelberg, 2007).

[25] Ruiz, D., Egea, J., Tomas-Barberan, F. A. & Gil, M. I. Carotenoids from new apricot (Prunus armeniaca L.) varieties and their relationship with flesh and skin color. J. Agric. Food Chem. 53, 6368-6374 (2005).

[26] Chen, Y. et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. GigaSci 7, 1-6 (2018).

[27] Mascher, M. et al. A chromosome conformation capture ordered sequence of the barley genome. Nature 544, 427-433 (2017).

[28] Servant, N. et al. HiC-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 16, 259 (2015).

[29] Gambino, G., Perrone, I. & Gribaudo, I. A rapid and effective method for RNA extraction from different tissues of grapevine and other woody plants. Phytochem. Anal. 19, 520-525 (2008).

[30] Marcais, G. & Kingsford, C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 27, 764-770 (2011).

[31] Doležel, J., Greilhuber, J. & Suda, J. Estimation of nuclear DNA content in plants using flow cytometry. Nat. Protoc. 2, 2233-2244 (2007).

[32] Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722-736 (2017).

[33] Vaser, R., Sovic, I., Nagarajan, N. & Sikic, M. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 27, 737-746 (2017).

[34] Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963 (2014).

[35] Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357-359 (2012).

[36] Durand, N. C. et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 3, 95-98 (2016).

[37] Dudchenko, O. et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science 356, 92-95 (2017).

[38] Tarailo-Graovac, M. & Chen, N. Using RepeatMasker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinforma. 4, 10 (2009).

[39] Bao, W., Kojima, K. K. & Kohany, O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob. DNA 6, 11 (2015).

[40] Sánchez-Pérez, R. et al. Mutation of a bHLH transcription factor allowed almond domestication. Science 364, 1095-1098 (2019).

[41] Baek, S. et al. Draft genome sequence of wild Prunus yedoensis reveals massive inter-specific hybridization between sympatric flowering cherries. Genome Biol. 19, 127 (2018).

[42] Velasco, R. et al. The genome of the domesticated apple (Malus x domestica Borkh.). Nat. Genet. 42, 833-839 (2010).

[43] Arabidopsis Genome, I. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 408, 796-815 (2000).

[44] Kim, D., Langmead, B. & Salzberg, S. L. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357-360 (2015).

[45] Pertea, M., Kim, D., Pertea, G. M., Leek, J. T. & Salzberg, S. L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650-1667 (2016).

[46] Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290-295 (2015).

[47] Haas, B. J. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments. Genome Biol. 9, R7 (2008).

[48] Conesa, A. et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674-3676 (2005).

[49] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

[50] McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297-1303 (2010).

[51] Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655-1664 (2009).

[52] Tajima, F. Evolutionary relationship of DNA sequences in finite populations. Genetics 105, 437-460 (1983).

[53] Nei, M. Analysis of gene diversity in subdivided populations. Proc. Natl Acad. Sci. USA 70, 3321-3323 (1973).

[54] Zeng, X. et al. Origin and evolution of qingke barley in Tibet. Nat. Commun. 9, 5433 (2018).

[55] Patterson, N. et al. Ancient admixture in human history. Genetics 192, 1065-1093 (2012).

[56] Martin, S. H., Davey, J. W. & Jiggins, C. D. Evaluating the use of ABBA-BABA statistics to locate introgressed loci. Mol. Biol. Evol. 32, 244-257 (2015).

[57] Yu, Y. et al. Genome re-sequencing reveals the evolutionary history of peach fruit edibility. Nat. Commun. 9, 5404 (2018).

[58] Pickrell, J. K. & Pritchard, J. K. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 8, e1002967 (2012).

[59] Hochberg, Y. & Benjamini, Y. More powerful procedures for multiple significance testing. Stat. Med. 9, 811-818 (1990).