Ginkgo biloba is an economically valuable tree worldwide. The species has nearly become extinct during the Quaternary, which has likely resulted in reduction of its genetic variability. The genetic variability is now conserved in few natural populations in China and a number of cultivars that are, however, derived from a few ancient trees, helping the species survive in China through medieval times. Despite the recent interest in ginkgo, however, detailed knowledge of its genetic diversity, conserved in cultivated trees and cultivars, has remained poor. This limits efficient conservation of its diversity as well as efficient use of the existing germplasm resources. Here we performed genotyping-by-sequencing (GBS) on 102 cultivated germplasms of ginkgo collected to explore their genetic structure, kinship, and inbreeding prediction. For the first time in ginkgo, a genome-wide association analysis study (GWAS) was used to attempt gene mapping of seed traits. The results showed that most of the germplasms did not show any obvious genetic relationship. The size of the ginkgo germplasm population expanded significantly around 1500 years ago during the Sui and Tang dynasties. Classification of seed cultivars based on a phylogenetic perspective does not support the current classification criteria based on phenotype. Twenty-four candidate genes were localized after performing GWAS on the seed traits. Overall, this study reveals the genetic basis of ginkgo seed traits and provides insights into its cultivation history. These findings will facilitate the conservation and utilization of the domesticated germplasms of this living fossil plant.

Co Innovation Center for Sustainable Forestry in Southern China Nanjing Forestry University Nanjing 210037 China

Coconut Research Institute of Chinese Academy of Tropical Agricultural Science Wenchang Hainan 571339 China

Department of Botany and Zoology Faculty of Science Masaryk University Koltlářská 2 Brno 61137 Czech Republic

Institute of Grassland Flowers and Ecology Beijing Academy of Agriculture and Forestry Sciences Beijing 100097 China

Zobrazit více v PubMed

Tang CQ, Yang Y, Ohsawa Met al. . Evidence for the persistence of wild Ginkgo biloba (Ginkgoaceae) populations in the Dalou Mountains, southwestern China. Am J Bot. 2012;99:1408–14 PubMed

Zhao Y, Paule J, Fu Cet al. . Out of China: distribution history of Ginkgo biloba L. Taxon. 2010;59:495–504

Zhao YP, Fan G, Yin PPet al. . Resequencing 545 ginkgo genomes across the world reveals the evolutionary history of the living fossil. Nat Commun. 2019;10:4201. PubMed PMC

Forest F, Moat J, Baloch Eet al. . Gymnosperms on the EDGE. Sci Rep. 2018;8:6053. PubMed PMC

Purohit VK, Phondani PC, Dhyani Det al. . Ginkgo biloba a living fossil: need conservation initiatives. Natl Acad Sci Lett. 2007;30:31–3

Sadaf HM, Bibi Y, Arshad Met al. . Status of maiden hair tree – Ginkgo biloba; living fossils becoming endangered. Nova J Med Biol Sci. 2014;2:1–5

Crane PR. Ginkgo: The Tree that Time Forgot. London: Yale University Press, 2013

Tredici P. Ginkgos and people – a thousand years of interaction. Arnoldia (Boston). 1991;51:2–15

Gong W, Zeng Z, Chen YYet al. . Glacial refugia of Ginkgo biloba and human impact on its genetic diversity: evidence from chloroplast DNA. J Integr Plant Biol. 2008;50:368–74 PubMed

Crane PR. An evolutionary and cultural biography of ginkgo. Plants People Planet. 2019;1:32–7

Zhou Q, Mu KM, Ni ZXet al. . Analysis of genetic diversity of ancient ginkgo populations using SSR markers. Ind Crop Prod. 2020;145:111942

Fan XX, Shen L, Zhang Xet al. . Assessing genetic diversity of Ginkgo biloba L. (Ginkgoaceae) populations from China by RAPD markers. Biochem Genet. 2004;42:269–78 PubMed

Elshire RJ, Glaubitz JC, Sun Qet al. . A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011;6:e19379. PubMed PMC

Chung YS, Choi SC, Jun T-Het al. . Genotyping-by-sequencing: a promising tool for plant genetics research and breeding. Hortic Environ Biotechnol. 2017;58:425–31

Riangwong K, Wanchana S, Aesomnuk Wet al. . Mining and validation of novel genotyping-by-sequencing (GBS)-based simple sequence repeats (SSRs) and their application for the estimation of the genetic diversity and population structure of coconuts (Cocos nucifera L.) in Thailand. Hortic Res. 2020;7:156. PubMed PMC

Pavan S, Delvento C, Mazzeo Ret al. . Almond diversity and homozygosity define structure, kinship, inbreeding, and linkage disequilibrium in cultivated germplasm, and reveal genomic associations with nut and seed weight. Hortic Res. 2021;8:15. PubMed PMC

Wang X, Bao K, Reddy UKet al. . The USDA cucumber (Cucumis sativus L.) collection: genetic diversity, population structure, genome-wide association studies, and core collection development. Hortic Res. 2018;5:64. PubMed PMC

Blanca J, Canizares J, Cordero Let al. . Variation revealed by SNP genotyping and morphology provides insight into the origin of the tomato. PLoS One. 2012;7:e48198. PubMed PMC

Migicovsky Z, Gardner KM, Richards Cet al. . Genomic consequences of apple improvement. Hort Res. 2021;8:9 PubMed PMC

Fu PN, Wu W, Lai GTet al. . Identifying Plasmopara viticola resistance loci in grapevine (Vitis amurensis) via genotyping-by-sequencing-based QTL mapping. Plant Physiol Biochem. 2020;154:75–84 PubMed

Liu H, Wang X, Wang Get al. . The nearly complete genome of Ginkgo biloba illuminates gymnosperm evolution. Nat Plants. 2021;7:748–56 PubMed

Xu X, Liu X, Ge Set al. . Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes. Nat Biotechnol. 2012;30:105–11 PubMed

Lam HM, Xu X, Liu Xet al. . Resequencing of 31 wild and cultivated soybean genomes identifies patterns of genetic diversity and selection. Nat Genet. 2010;42:1053–9 PubMed

Zhang MY, Xue C, Hu Het al. . Genome-wide association studies provide insights into the genetic determination of fruit traits of pear. Nat Commun. 2021;12:1144. PubMed PMC

Duan N, Bai Y, Sun Het al. . Genome re-sequencing reveals the history of apple and supports a two-stage model for fruit enlargement. Nat Commun. 2017;8:1–11 PubMed PMC

Forutan M, Ansari Mahyari S, Baes Cet al. . Inbreeding and runs of homozygosity before and after genomic selection in North American Holstein cattle. BMC Genomics. 2018;19:98. PubMed PMC

Pemberton TJ, Absher D, Feldman MWet al. . Genomic patterns of homozygosity in worldwide human populations. Am J Hum Genet. 2012;91:275–92 PubMed PMC

Hu Y, Šmarda P, Liu Get al. . High-depth transcriptome reveals differences in natural haploid Ginkgo biloba L. due to the effect of reduced gene dosage. Int J Mol Sci. 2022;23:8958. PubMed PMC

Terhorst J, Kamm JA, Song YS. Robust and scalable inference of population history from hundreds of unphased whole genomes. Nat Genet. 2017;49:303–9 PubMed PMC

Hamrick JL, Godt MJW. Allozyme diversity in cultivated crops. Crop Sci. 2002;37:26–30

Šmarda P, Horová L, Knápek Oet al. . Multiple haploids, triploids, and tetraploids found in modern-day ‘living fossil’ Ginkgo biloba. Hortic Res. 2018;5:55. PubMed PMC

Šmarda P, Veselý P, Šmerda Jet al. . Polyploidy in a ‘living fossil’ Ginkgo biloba. New Phytol. 2016;212:11–4 PubMed

De Vladar HP, Santos M, Szathmary E. Grand views of evolution. Trends Ecol Evol. 2017;32:324–34 PubMed

Kumar S, Deng CH, Hunt Met al. . Homozygosity mapping reveals population history and trait architecture in self-incompatible pear (Pyrus spp.). Front Plant Sci. 2021;11:590846. PubMed PMC

Yang XP, Luo ZL, Todd Jet al. . Genome-wide association study of multiple yield traits in a diversity panel of polyploid sugarcane (Saccharum spp.). Plant. Genome. 2020;13:e20006. PubMed

Liu S, Zhang Y, Feng Qet al. . Tomato AUXIN RESPONSE FACTOR 5 regulates fruit set and development via the mediation of auxin and gibberellin signaling. Sci Rep. 2018;8:2971. PubMed PMC

Luna CM, Pastori GM, Driscoll Set al. . Drought controls on H2O2 accumulation, catalase (CAT) activity and CAT gene expression in wheat. J Exp Bot. 2005;56:417–23 PubMed

Boateng ID. A critical review of current technologies used to reduce ginkgotoxin, ginkgotoxin-5′-glucoside, ginkgolic acid, allergic glycoprotein, and cyanide in Ginkgo biloba L. seed. Food Chem. 2022;382:132408 PubMed

Muller BSF, Almeida JE, Lima BMet al. . Independent and joint-GWAS for growth traits in Eucalyptus by assembling genome-wide data for 3373 individuals across four breeding populations. New Phytol. 2019;221:818–33 PubMed

Minamikawa MF, Takada N, Terakami Set al. . Genome-wide association study and genomic prediction using parental and breeding populations of Japanese pear (Pyrus pyrifolia Nakai). Sci Rep. 2018;8:11994. PubMed PMC

Lee SJ, Ban SH, Kim GHet al. . Identification of potential gene-associated major traits using GBS-GWAS for Korean apple germplasm collections. Plant Breed. 2017;136:977–86

Yu BY, Boyle K, Zhang WTet al. . Multi-trait and multi-environment QTL analysis reveals the impact of seed colour on seed composition traits in Brassica napus. Mol Breed. 2016;36:111.

Doyle JJ. Isolation of plant DNA from fresh tissue. Focus. 1990;12:13–5

Chen SF, Zhou YQ, Chen YRet al. . Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:i884–90 PubMed PMC

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

Van der Auwera GA, Carneiro MO, Hartl Cet al. . From FastQ data to high-confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;43:11.10.11-11.10.33 PubMed PMC

Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164–4 PubMed PMC

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

Chang CC, Chow CC, Tellier Let al. . Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7. PubMed PMC

Francis RM. POPHELPER: an R package and web app to analyse and visualize population structure. Mol Ecol Resour. 2017;17:27–32 PubMed

Price AL, Patterson NJ, Plenge RMet al. . Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–9 PubMed

Yang JA, Lee SH, Goddard MEet al. . GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88:76–82 PubMed PMC

Zhang C, Dong SS, Xu JYet al. . PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2019;35:1786–8 PubMed

Purcell S, Neale B, Todd-Brown Ket al. . PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75 PubMed PMC

Danecek P, Auton A, Abecasis Get al. . The variant call format and VCFtools. Bioinformatics. 2011;27:2156–8 PubMed PMC

Bates D, Mächler M, Bolker Bet al. . Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48

Massey FJ Jr. The Kolmogorov-Smirnov test for goodness of fit. J Am Stat Assoc. 1951;46:68–78

Villanueva RAM, Chen ZJ. ggplot2: Elegant Graphics for Data Analysis. Measurement: Interdisciplinary Research and Perspectives. 2019;17:160–167

Jiang L, Zheng Z, Qi Tet al. . A resource-efficient tool for mixed model association analysis of large-scale data. Nat Genet. 2019;51:1749–55 PubMed

Mbatchou J, Barnard L, Backman Jet al. . Computationally efficient whole-genome regression for quantitative and binary traits. Nat Genet. 2021;53:1097–103 PubMed

Dey R, Schmidt EM, Abecasis GRet al. . A fast and accurate algorithm to test for binary phenotypes and its application to PheWAS. Am J Hum Genet. 2017;101:37–49 PubMed PMC

Zhou W, Nielsen JB, Fritsche LGet al. . Efficiently controlling for case-control imbalance and sample relatedness in large-scale genetic association studies. Nat Genet. 2018;50:1335–41 PubMed PMC

Chen H, Huffman JE, Brody JAet al. . Efficient variant set mixed model association tests for continuous and binary traits in large-scale whole-genome sequencing studies. Am J Hum Genet. 2019;104:260–74 PubMed PMC

Wu MC, Lee S, Cai Tet al. . Rare-variant association testing for sequencing data with the sequence kernel association test. Am J Hum Genet. 2011;89:82–93 PubMed PMC

Lee S, Wu MC, Lin X. Optimal tests for rare variant effects in sequencing association studies. Biostatistics. 2012;13:762–75 PubMed PMC