Comparative population genomics dissects the genetic basis of seven domestication traits in jujube

Abstract

Jujube (Ziziphus jujuba Mill.) is an important perennial fruit tree with a range of interesting horticultural traits. It was domesticated from wild jujube (Ziziphus acidojujuba), but the genomic variation dynamics and genetic changes underlying its horticultural traits during domestication are poorly understood. Here, we report a comprehensive genome variation map based on the resequencing of 350 accessions, including wild, semi-wild and cultivated jujube plants, at a >15× depth. Using the combination of a genome-wide association study (GWAS) and selective sweep analysis, we identified several candidate genes potentially involved in regulating seven domestication traits in jujube. For fruit shape and kernel shape, we integrated the GWAS approach with transcriptome profiling data, expression analysis and the transgenic validation of a candidate gene to identify a causal gene, ZjFS3, which encodes an ethylene-responsive transcription factor. Similarly, we identified a candidate gene for bearing-shoot length and the number of leaves per bearing shoot and two candidate genes for the seed-setting rate using GWAS. In the selective sweep analysis, we also discovered several putative genes for the presence of prickles on bearing shoots and the postharvest shelf life of fleshy fruits. This study outlines the genetic basis of jujube domestication and evolution and provides a rich genomic resource for mining other horticulturally important genes in jujube.

References

[1] Qu, Z. & Wang, Y. Chinese Fruit Trees Record-Chinese Jujube (China Forestry Publishing House, 1993).

[2] Liu, M., Liu, P. & Liu, G. Advances of research on germplasm resources of Chinese jujube. Acta Hortic. 993, 15-20 (2013).

[3] Liu, M. J. Chinese jujube: botany and horticulture. Hort. Revi 32, 229-298 (2006).

[4] Huang, J. et al. The jujube genome provides insights into genome evolution and the domestication of sweetness/acidity taste in fruit trees. PLoS Genet. 12, e1006433 (2016).

[5] Wang, S. et al. The OsSPL16-GW7 regulatory module determines grain shape and simultaneously improves rice yield and grain quality. Nat. Genet. 47, 949-954 (2015).

[6] Huang, K. et al. WIDE AND THICK GRAIN 1, which encodes an otubain-like protease with deubiquitination activity, influences grain size and shape in rice. Plant J. 91, 849-860 (2017).

[7] Zhao, D. S. et al. GS9 acts as a transcriptional activator to regulate rice grain shape and appearance quality. Nat. Commun. 9, 1240 (2018).

[8] Liu, J., Van. Eck., J., Cong, B. & Tanksley, S. D. A new class of regulatory genes underlying the cause of pear-shaped tomato fruit. Proc. Natl Acad. Sci. USA 99, 13302-13306 (2002).

[9] Wu, S., Xiao, H., Cabrera, A., Meulia, T. & van der Knaap, E. SUN regulates vegetative and reproductive organ shape by changing cell division patterns. Plant Physiol. 157, 1175-1186 (2011).

[10] Zhao, J. et al. A functional allele of CsFUL1 regulates fruit length through repressing CsSUP and inhibiting auxin transport in cucumber. Plant Cell 31, 1289-1307 (2019).

[11] Liu, M. et al. The complex jujube genome provides insights into fruit tree biology. Nat. Commun. 5, 5315 (2014).

[12] Kellogg, A. A., Branaman, T. J., Jones, N. M., Little, C. Z. & Swanson, J. D. Morphological studies of developing Rubus prickles suggest that they are modifed glandular trichomes. Botany 89, 217-226 (2011).

[13] Ma, Z. Y., Wen, J., Ickert-Bond, S. M., Chen, L. Q. & Liu, X. Q. Morphology, structure, and ontogeny of trichomes of the grape genus (Vitis, Vitaceae). Front. Plant Sci. 7, 704 (2016).

[14] Miller, A. J. & Gross, B. L. From forest to field: perennial fruit crop domestication. Am. J. Bot. 98, 389-414 (2011).

[15] Brummell, D. A. et al. Modification of expansin protein abundance in tomato fruit alters softening and cell wall polymer metabolism during ripening. Plant Cell 11, 2203-2216 (1999).

[16] Brummell, D. A. & Harpster, M. H. Cell wall metabolismin fruit softening and quality and its manipulation in transgenic plants. Plant Mol. Biol. 47, 311-340 (2001).

[17] Goulao, L. F. & Oliveira, C. M. Cell wall modification during fruit ripening: when a fruit is not the fruit. Trends Food Sci. Tech. 19, 4-25 (2008).

[18] Mercado, J. A., Pliego-Alfaro, F. & Quesada, M. A. in Breeding for Fruit Quality (eds Jenks, M. A. & Bebeli, P. J.) 81-104 (Oxford: JohnWiley & Sons, 2011).

[19] Sheehy, R. E., Kramer, M. & Hiatt, W. R. Reduction of polygalacturonase activity in tomato fruit by antisense RNA. Proc. Natl Acad. Sci. USA 85, 8805-8809 (1988).

[20] Smith, C. et al. Antisense RNA inhibition of polygalacturonase gene expression in transgenic tomatoes. Nature 334, 724-726 (1988).

[21] Jime´nez-Bermu´dez, S. et al. Manipulation of strawberry fruit softening by antisense expression of a pectate lyase gene. Plant Physiol. 128, 751-759 (2002).

[22] Quesada, M. A. et al. Antisensedown-regulation of the FaPG1 gene reveals an unexpected central role for polygalacturonase in strawberry fruit softening. Plant Physiol. 150, 1022-1032 (2009).

[23] Atkinson, R. G. et al. Down-regulation of POLYGALACTURONASE1 alters firm- ness, tensile strength and water loss in apple (Malus × domestica) fruit. BMC Plant Biol. 12, 129-142 (2012).

[24] Duan, N. et al. Genome re-sequencing reveals the history of apple and sup- ports a two-stage model for fruit enlargement. Nat. Commun. 8, 249 (2017).

[25] Liang, Z. et al. Whole-genome resequencing of 472 Vitis accessions for grapevine diversity and demographic history analyses. Nat. Commun. 10, 1190 (2019).

[26] Wang, X. et al. Genomic analyses of primitive, wild and cultivated citrus provide insights into asexual reproduction. Nat. Genet. 49, 765-772 (2017).

[27] Wang, L. et al. Genome of wild mandarin and domestication history of mandarin. Mol. Plant 11, 1024-1037 (2018).

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

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

[30] Li, Y. et al. Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history. Genome Biol. 20, 36 (2019).

[31] Cao, K. et al. Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach. Plant Biotechnol. J. 17, 1954-1970 (2019).

[32] Chen, H., Patterson, N. & Reich, D. Population differentiation as a test for selective sweeps. Genome Res. 20, 393-402 (2010).

[33] Yang, X. et al. Rice ferredoxin-dependent glutamate synthase regulates nitrogen-carbon metabolomes and is genetically gifferentiated between japonica and indica subspecies. Mol. Plant 9, 1520-1534 (2016).

[34] Bi, Z. et al. ES7, encoding a ferredoxin-dependent glutamate synthase, func- tions in nitrogen metabolism and impacts leaf senescence in rice. Plant Sci. 259, 24-34 (2017).

[35] Kasili, R. et al. BRANCHLESS TRICHOMES links cell shape and cell cycle control in Arabidopsis trichomes. Development 138, 2379-2388 (2011).

[36] Peterson, K. M. et al. Arabidopsis homeodomain-leucine zipper IV proteins promote stomatal development and ectopically induce stomata beyond the epidermis. Development 140, 1924-1935 (2013).

[37] Berg, M., Rogers, R., Muralla, R. & Meinke, D. Requirement of aminoacyl-tRNA synthetases for gametogenesis and embryo development in Arabidopsis. Plant J. 44, 866-878 (2005).

[38] Wang, T. et al. A receptor heteromer mediates the male perception of female attractants in plants. Nature 531, 241-244 (2016).

[39] Durrant, W. E., Wang, S. & Dong, X. Arabidopsis SNI1 and RAD51D regulate both gene transcription and DNA recombination during the defense response. Proc. Natl Acad. Sci. USA 104, 4223-4227 (2007).

[40] Lukowitz, W., Roeder, A., Parmenter, D. & Somerville, C. A MAPKK kinase gene regulates extra-embryonic cell fate in Arabidopsis. Cell 116, 109-119 (2004).

[41] Hao, J., Jin, Z., Wang, Y., Zhang, B. & Li., D. Study on embryo develop- ment and somatic embryogenesis of jujube. J. Mol. Cell Biol. 39, 423-430 (2006).

[42] Wang, Y. in Study on Pollination Biology of Chinese Jujube[D]. 18-30 (Agri- cultural University of Hebei, Hebei, 2008).

[43] Li, D. et al. Observation of embryo development and abortion in Chinese jujube (Ziziphus jujuba Mill.). Acta Agriculturae Boreal.-Occidentalis Sin. 25, 1379-1385 (2016).

[44] Xu, J., Hua, K. & Lang, Z. Genome editing for horticultural crop improvement. Hortic. Res. 6, 113 (2019).

[45] Chang, L., Wu, S. & Tian, L. Effective genome editing and identification of a regiospecific gallic acid 4-O-glycosyltransferase in pomegranate (Punica granatum L.). Hortic. Res. 6, 123 (2019).

[46] Liu, M. & Wang, M. Germplasm Resources of Chinese Jujube (China Forestry Publishing House, 2009).

[47] Li, D. Descriptors and Data Standard for Jujube (Ziziphus jujuba Mill.) (China Agriculture Press, 2006).

[48] Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078-2079 (2009).

[49] Vilella, A. J. et al. EnsemblCompara Gene Trees: complete, duplication-aware phylogenetic trees in vertebrates. Genome Res. 19, 327-335 (2009).

[50] Yang, J. A., Lee, S. H., Goddard, M. E. & Visscher, P. M. GCTA: a tool for genome- wide complex trait analysis. Am. J. Hum. Genet. 88, 76-82 (2011).

[51] Barrett, J. C., Fry, B., Maller, J. & Daly, M. J. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263-265 (2005).

[52] Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156-2158 (2011).

[53] Zhou, X. & Stephens, M. Genome-wide efficient mixed-model analysis for association studies. Nat. Genet. 44, 821-824 (2012).

[54] Li, M. X., Yeung, J. M., Cherny, S. S. & Sham, P. C. Evaluating the effective numbers of independent tests and significant p-value thresholds in com- mercial genotyping arrays and public imputation reference datasets. Hum. Genet. 131, 747-756 (2012).

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

[56] Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real time quantitative PCR and the 2^−ΔΔCT method . Methods 25, 402-408 (2011).

[57] Ning, J., Li, X., Hicks, L. M. & Xiong, L. A Raf-like MAPKKK gene DSM1 mediates drought resistance through reactive oxygen species scavenging in rice. Plant Physiol. 152, 876-890 (2010).