Chromosome-scale assembly and population diversity analyses provide insights into the evolution of Sapindus mukorossi

Abstract

Sapindus mukorossi is an environmentally friendly plant and renewable energy source whose fruit has been widely used for biomedicine, biodiesel, and biological chemicals due to its richness in saponin and oil contents. Here, we report the first chromosome-scale genome assembly of S. mukorossi (covering ∼391 Mb with a scaffold N50 of 24.66 Mb) and characterize its genetic architecture and evolution by resequencing 104 S. mukorossi accessions. Population genetic analyses showed that genetic diversity in the southwestern distribution area was relatively higher than that in the northeastern distribution area. Gene flow events indicated that southwest species may be the donor population for the distribution areas in China. Genome-wide selective sweep analysis showed that a large number of genes are involved in defense responses, growth and development, including SmRPS2 , SmRPS4 , SmRPS7 , SmNAC2 , SmNAC23 , SmNAC102 , SmWRKY6 , SmWRKY26 , and SmWRKY33 . We also identified several candidate genes controlling six agronomic traits by genome-wide association studies, including SmPCBP2 , SmbHLH1 , SmCSLD1 , SmPP2C , SmLRR-RKs , and SmAHP . Our study not only provides a rich genomic resource for further basic research on Sapindaceae woody trees but also identifies several economically significant genes for genomics-enabled improvements in molecular breeding.

References

[1] Goyal S, Kumar D, Menaria G et al. Medicinal plants of the genus Sapindus (Sapindaceae) - a review of their botany, phytochemistry, biological activity and traditional uses . J Drug Deliv Ther. 2014; 4: 7-20.

[2] Li SZ . The Compendium of Materia Medica . People’s Medical Publishing House; 1982.

[3] Zhu S, Fang PJ . Beijing: People’s Medical Publishing House, 1959.

[4] Shao WH . Geographic variation of saponins contents in Sapindus mukorossi peels from different habitats. Bull Bot Res. 2012; 32: 627-31.

[5] Rupeshkumar G, Surekha SK, Balu CA et al. Study of functional properties of Sapindus mukorossi as a potential biosurfactant. Indian J Sci Technol. 2011; 4: 530-3.

[6] Huang SM, Wang JW, Du MH et al. Fatty acid composition analysis of Sapindus mukorossi Gaerth. seed oil. China Oils Fats. 2009; 34: 74-6.

[7] Sun S, Ke X, Cui L et al. Enzymatic epoxidation of Sapindus mukorossi seed oil by perstearic acid optimized using response surface methodology. Ind Crops Prod. 2011; 33: 676-82.

[8] Kanchanapoom T, Kasai R, Yamasaki K . Acetylated triterpenes saponins from the Thai medicinal plant. Biol Pharm Bull. 2001; 49: 1195-7.

[9] Guo Y, Jian-ping XIE, Ai-hua LIU et al. Overview of the pharmacological action of Sapindus mukorossi . J Path Biol. 2011; 6: 873-4.

[10] Shao WH, Yue HF, Jiang JM et al. Study on genetic variation of seedling growth from different provenances of Sapindus mukorossi . J ZheJiang For Sci Technol. 2012; 32: 21-5.

[11] Fan HH, Yao XM, Tang XH et al. A study on superior individual comprehensive appraisal of Sapindus mukurossi . J Fujian For Sci Technol 2014; 41: 1-6.

[12] Huang S, Pan WM . Good landscape trees: Sapindus mukorossi Gaertn. Garden. 2004; 9: 42.

[13] Zhang FL . Study on tissue culture for soapberry. J Mountain Agr Biol. 2005; 24: 119-23.

[14] Wei FY, Fang C . Enzyme-assisted aqueous extraction of Sapindus saponins . Appl Chem Ind. 2010; 39: 1149-51.

[15] Kanchanapoom T, Kasai R, Yamasaki K . Acetylated triterpene saponins from the Thai medicinal plant. Chem Pharm Bull (Tokyo). 2001; 49: 1195-7.

[16] Diao SF, Shao WH, Chen T et al. Study on variation in phenotypic traits of fruit and seed of Sapindus mukorossi in seedling plantation. J Northwest A&F Univ. 2016; 29: 176-82.

[17] Liu J, Liu S, Xu Y et al. Screening of Sapindus germplasm resources in China based on agro-morphological traits. J For Res. 2021; 21: 1350-8.

[18] Mahar KS, Palni LMS, Ranade SA et al. Molecular analyses of genetic variation and phylogenetic relationship in Indian soap nut (Sapindus L.) and closely related taxa of the family Sapindaceae . Meta Gene. 2017; 13: 50-6.

[19] Hong L, Bai ME, Zhang JZ et al. Genetic diversity and relationship analysis of Sapindus mukorossi germplasm by ISSR and SRAP molecular markers. J Zhejiang Agr Sci. 2013; 5: 71-3.

[20] Diao SF, Shao WH, Chen T et al. Genetic diversity of Sapindus mukorossi natural populations in China based on ISSR. Eur For Inst Proc. 2016; 29: 176-82.

[21] Sun CW, Wang L, Liu J et al. Genetic structure and biogeographic divergence among Sapindus species: an inter-simple sequence repeat-based study of germplasms in China. Ind Crops Prod. 2018; 118: 1-10.

[22] Jiang CC, Lu XK, Ye XF . Molecular genetic diversity analysis of energy plant Sapindus mukorossi germplasm resources. Mol Plant Breed. 2016; 14: 2888-95.

[23] Zeng LX, Chen S, Yang YY et al. The resistance evaluation and their utilization of the germplasm resources from international rice bacterial blight nursery (IRBBN) in South China. J Plant Genet Res. 2013; 14: 930-5.

[24] Guo M, Zhang Z, Li S et al. Genomic analyses of diverse wild and cultivated accessions provide insights into the evolutionary history of jujube. Plant Biotechnol J. 2021; 19: 517-31.

[25] Zhou Z, Jiang Y, Wang Z et al. Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nat Biotechnol. 2015; 33: 408-14.

[26] Varshney RK, Thudi M, Roorkiwal M et al. Resequencing of 429 chickpea accessions from 45 countries provides insights into genome diversity, domestication and agronomic traits. Nat Genet. 2019; 51: 857-64.

[27] Li Y, Cao K, Zhu G et al. Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history. Genome Biol. 2019; 20: 36-54.

[28] Chen JH, Huang Y, Brachi B et al. Genome-wide analysis of cushion willow provides insights into alpine plant divergence in a biodiversity hotspot. Nat Commun. 2019; 10: 5230-42.

[29] Wang XC, Feng H, Chang Y et al. Population sequencing enhances understanding of tea plant evolution. Nat Commun. 2020; 11: 4447-57.

[30] Lin YL, Min J, Lai R et al. Genome-wide sequencing of longan (Dimocarpus longan Lour.) provides insights into molecular basis of its polyphenol-rich characteristics . GigaScience. 2017; 6: 1-14. https://doi.org/10.1093/gigascience/gix023.

[31] Yang J, Wariss HM, Tao L et al. De novo genome assembly of the endangered Acer yangbiense, a plant species with extremely small populations endemic to Yunnan province, China. GigaScience. 2019; 8. https://doi.org/10.1093/gigascience/giz085.

[32] Huang SX, Diang J, Deng D et al. Draft genome of the kiwifruit Actinidia chinensis . Nat Commun. 2013; 4: 1055-73.

[33] Ming R, VanBuren R, Wai CM et al. The pineapple genome and the evolution of CAM photosynthesis. Nat Genet. 2015; 47: 1435-42.

[34] Jaillon O, Aury JM, Noel B et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007; 449: 463-7.

[35] Ma Q, Feng K, Yang W et al. Identification and characterization of nucleotide variations in the genome of Ziziphus jujuba (Rhamnaceae) by next generation sequencing. Mol Biol Rep. 2014; 41: 3219-23.

[36] Xu Q, Chen LL, Ruan X et al. The draft genome of sweet orange (Citrus sinensis) . Nat Genet. 2013; 45: 59-66.

[37] Wu J, Wang Z, Shi Z et al. The genome of the pear (Pyrus bretschneideri Rehd.) . Genome Res. 2013; 23: 396-408.

[38] Riccardo V, Chagne D, Merriman TR et al. The genome of the domesticated apple (Malus × domestica Borkh.) . Nat Genet. 2010; 42: 833-9.

[39] Petit JR, Jouzel J, Raynaud D et al. Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica. Nature. 1999; 399: 429-36.

[40] Mackey D, Holt BF, Wiig A et al. RIN4 interacts with Pseudomonas syringae type III effector molecules and is required for RPM1-mediated resistance in Arabidopsis . Cell. 2002; 108: 743-54.

[41] Zhang ZB, Wu Y, Gao M et al. Disruption of PAMP-induced MAP kinase cascade by a Pseudomonas syringae effector activates plant immunity mediated by the NBLRR protein SUMM2. Cell Host Microbe. 2012; 11: 253-63.

[42] Bu QY, Jiang H, Li CB et al. Role of the Arabidopsis thaliana NAC transcription factors ANAC019 and ANAC055 in regulating jasmonic acid-signaled defense responses. Cell Res. 2008; 18: 756-67.

[43] Sun L, Zhang H, Li D et al. Functions of rice NAC transcriptional factors, ONAC122 and ONAC131, in defense responses against Magnaporthe grisea . Plant Mol Biol. 2013; 81: 41-56.

[44] Xue T, Zheng XH, Chen D et al. A high-quality genome provides insights into the new taxonomic status and genomic characteristics of Cladopus chinensis (Podostemaceae) . Hortic Res. 2020; 7: 46.

[45] Li JB, Luan YS, Jin H . The tomato SlWRKY gene plays an important role in the regulation of defense responses in tobacco. Biochem Biophys Res Commun. 2012; 427: 671-6.

[46] Cheng J, Frishman WH, Aronow WS . Updates on cytochrome P450-mediated cardiovascular drug interactions. Dis Mon. 2010; 56: 163-79.

[47] Wickham S, West MB, Cook PF et al. Gamma-glutamyl compounds: substrate specificity of gamma-glutamyl transpeptidase enzymes. Anal Biochem. 2011; 414: 208-14.

[48] Yue M, Szostkiewicz I, Korte A et al. Regulators of PP2C phosphatase activity function as abscisic acid sensors. Science. 2009; 324: 1064-8.

[49] Shiota H, Ko S, Wada S et al. A carrot G-box binding factor-type basic region/leucine zipper factor DcBZ1 is involved in abscisic acid signal transduction in somatic embryogenesis. Plant Physiol Biochem. 2008; 46: 550-8.

[50] Qin Y, Xue B, Liu C et al. NLRX1 mediates MAVS degradation to attenuate hepatitis C virus-induced innate immune response through PCBP2. J Virol. 2017; 91: 1264-81.

[51] Yan HF, Zhou H, Luo H et al. Characterization of full-length transcriptome in Saccharum officinarum and molecular insights into tiller development. BMC Plant Biol. 2021; 21: 228.

[52] Xin XF, Kvitko B, He SY . Pseudomonas syringae: what it takes to be a pathogen . Nat Rev Microbiol. 2018; 16: 316-28.

[53] Pan C, Tian K, Ban Q et al. Genome-wide analysis of the biosynthesis and deactivation of gibberellin-dioxygenases gene family in Camellia sinensis (L.) O. Kuntze. Genes. 2017; 8: 235.

[54] Li F, Wei CY, Qiao C et al. Three duplication events and variable molecular evolution characteristics involved in multiple GGPS genes of six Solanaceae species. J Genet. 2016; 95: 453-7.

[55] Li C, Zhu Y, Guo X et al. Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng C. A. Meyer. BMC Genomics 2013; 14: 245.

[56] Guo J, Zhou YJ, Hillwig ML et al. CYP76AH1 catalyzes turnover of miltiradiene in tanshinones biosynthesis and enables heterologous production of ferruginol in yeasts. Proc Natl Acad Sci USA. 2013; 110: 12108-13.

[57] Li Y, Gou J, Chen F et al. Comparative transcriptome analysis identifies putative genes involved in the biosynthesis of xanthanolides in Xanthium strumarium L. Front Plant Sci. 2016; 7: 1317.

[58] Torii KU . Leucine-rich repeat receptor kinases in plants: structure, function, and signal transduction pathways. Int Rev Cytol. 2004; 234: 1-46.

[59] Saxena IM, Brown MR . Biochemistry and molecular biology of cellulose biosynthesis in plants: prospects for genetic engineering. Bioeng Mol Biol Plant Pathw. 2008; 1: 135-60.

[60] David S, Dean D . Elevating the vitamin E content of plants through metabolic engineering. Science. 1998; 282: 2098-100.

[61] Mertens J, Pollier J, Bossche RV et al. The bHLH transcription factors TSAR1 and TSAR2 regulate triterpene saponin biosynthesis in Medicago truncatula . Plant Physiol. 2016; 170: 194-210.

[62] Yu F, Zhao X, Chai J et al. Chromosome-specific painting unveils chromosomal fusions and distinct allopolyploid species in the saccharum complex. New Phytol. 2022; 233: 1953-65.

[63] Huang Y, Ding W, Zhang M et al. The formation and evolution of centromeric satellite repeats in Saccharum species. Plant J. 2021; 106: 616-29.

[64] Guillaume C, Rugen M, Bockmayr A et al. A model-based method for investigating bioenergetic processes in autotrophically growing eukaryotic microalgae: application to the green algae Chlamydomonas reinhardtii . Biotechnol Prog. 2011; 27: 631-40.

[65] Chen DH, Ronald PC . A rapid DNA minipreparation method suitable for AFLP and other PCR applications. Plant Mol Biol Report. 1999; 17: 53-7.

[66] Korbel JO, Lee C . Genome assembly and haplotyping with Hi-C. Nat Biotechnol. 2013; 31: 1099-101.

[67] Camacho C, Coulouris G, Avagyan V et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10: 421.

[68] Langdon WB . Performance of genetic programming optimised Bowtie2 on genome comparison and analytic testing (GCAT) benchmarks. BioData Min. 2015; 8: 1-7.

[69] Flynn JM, Hubley R, Goubert C et al. RepeatModeler2 for automated genomic discovery of transposable element families. Proc Natl Acad Sci USA. 2020; 117: 9451-7.

[70] Tarailo-Graovac M, Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009; 25: 1-14.

[71] Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Current Protoc Bioinformatics. 2004; 4: 10.1-10.14.

[72] Xu Z, Wang H . LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007; 35: W265-8.

[73] Lowe TM, Chan PP . tRNAscan-SE on-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016; 44: W54-7.

[74] Nawrocki EP, Kolbe DL, Eddy SR . Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009; 25: 1335-7.

[75] Kondrashov N, Pusic A, Stumpf CR et al. Ribosome-mediated specificity in Hox mRNA translation and vertebrate tissue patterning. Cell. 2011; 145: 383-97.

[76] Hoff KJ, Lange S, Lomsadze A et al. BRAKER1: unsupervised RNA-Seq-based genome annotation with GeneMark-ET and AUGUSTUS. Bioinformatics. 2016; 32: 767-9.

[77] Stanke M, Keller O, Gunduz I et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006; 34: W435-9.

[78] Ashburner M, Ball CA, Blake JA et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000; 25: 25-9.

[79] Kanehisa M, Goto S, Sato Y et al. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014; 42: D199-205.

[80] Tatusov RL, Galperin MY, Natale DA et al. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000; 28: 33-6.

[81] Boeckmann B, Bairoch A, Apweiler R et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003; 31: 365-70.

[82] Jensen L, Julien P, Kuhn M et al. eggNOG: automated construction and annotation of orthologous groups of genes. Nucleic Acids Res. 2007; 36: D250-4.

[83] Liu K, Randal LC, Tandy W . RAxML and FastTree: comparing two methods for large-scale maximum likelihood phylogeny estimation. PLoS One. 2011; 6: e27731.

[84] Puttick MN . MCMCtreeR: functions to prepare MCMCtree analyses and visualize posterior ages on trees. Bioinformatics. 2019; 35: 5321-2.

[85] Xu Y, Hou H, Liu Q et al. Removal behavior research of orthophosphate by CaFe-layered double hydroxides. Desalin Water Treat. 2016; 57: 7918-25.

[86] 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.

[87] Zhang Z, Li J, Zhao XQ et al. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genom Proteom Bioinf. 2006; 4: 259-63.

[88] Emms DM, Kelly S . OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015; 16: 157.

[89] Edgar RC . MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32: 1792-7.

[90] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20.

[91] Li H, Durbin R . Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[92] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[93] Cingolani P, Platts A, Wang LL et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118. Fly. 2012; 6: 80-92.

[94] Purcell S, Neale B, Todd-Brown K et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81: 559-75.

[95] Robinson MD, McCarthy DJ, Smyth GK . edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010; 26: 139-40.

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

[97] Beerenwinkel N, Rahnenfuhrer J, Kaiser R . Mtreemix: a software package for learning and using mixture models of mutagenetic trees. Bioinformatics. 2005; 21: 2106-7.

[98] Zhang C, Dong SS, Xu JY et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2019; 35: 1786-8.

[99] Danecek P, Auton A, Abecassis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[100] Vilella AJ, Blanco-Garcia A, Hutter S et al. VariScan: analysis of evolutionary patterns from large-scale DNA sequence polymorphism data. Bioinformatics. 2005; 21: 2791-3.

[101] Degiorgio M, Huber CD, Hubisz MJ et al. SWEEPFINDER2: increased sensitivity, robustness, and flexibility. Bioinformatics. 2016; 32: 1895-7.