Revealing the genetic regulation of wood traits and secondary cell wall development in Ginkgo biloba: an integrated analysis from the perspectives of GWAS, TWAS, and WGCNA

Abstract

Understanding the genetic and regulatory mechanisms underlying wood traits and secondary cell wall (SCW) development in Ginkgo biloba is crucial for improving wood quality. We identified key genes related to wood traits and SCW development through integrated genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and weighted gene co-expression network analysis (WGCNA). Cellulose biosynthesis in the SCW is catalyzed by the CesA4–CesA7–CesA8 complex encoded by GbCesA4, GbCesA7, and GbCesA8A/ 8B. These CesA genes form a co-expression network with TUBA/ TUBB and EG, indicating coordination among cellulose synthesis, cytoskeletal guidance, and cell wall remodeling. Additionally, loss of function of GbCesA8B caused only a slight reduction in cellulose content, supporting potential functional redundancy between GbCesA8A and GbCesA8B. For hemicellulose biosynthesis, GbCSLA9A/ 9B and IRX9/ IRX14 were major contributors to mannan/glucomannan and xylan synthesis, respectively, and formed a co-expression network with UXS, UXE, IRX7, GXMT, and URGT, spanning nucleotide sugar supply, transport, and polymer elongation and modification. Moreover, MYB46 may regulate mannan/glucomannan biosynthesis in the SCW by activating CSLA9 transcription. For lignin biosynthesis, TWAS identified multiple genes involved in phenylalanine biosynthesis, phenylpropanoid metabolism, and lignin monomer polymerization, including ADT/ PDT, PAL, and PER, as well as MYB91 and several bHLH genes that may positively regulate lignin accumulation. Furthermore, several transcription factors potentially involved in SCW development were identified, including GATA9 as a putative positive regulator, WRKY12 and HB15 as potential negative regulators, and ELF6, which may facilitate tracheid expansion. Our findings provide valuable insights into the genetic regulation of wood traits and SCW development in Ginkgo.

Acknowledgements

This study was supported by the Jiangsu Provincial Key Research and Development Program (BE2022373), the Independent Research Project of State Key Laboratory of Tree Genetics and Breeding (SKLTGBNJ2024–005), and the Graduate Research and Innovation Projects of Jiangsu Province (KYCX23_1245).

Author contributions

T.G., G.W., F.F., and F.C. designed the research; T.G., J.S., X.L., and Y.C. performed the experiments; T.G. and J.S. performed the upstream and downstream analysis of GWAS, TWAS, and WGCNA; T.G. and J.S. wrote the manuscript; and J.G., F.F., F.C., and G.W. revised the manuscript.

Data availability

The raw sequencing data of whole-genome resequencing and transcriptome sequencing were uploaded to the National Genomics Data Center (NGDC, https://ngdc.cncb.ac.cn/) with the accession code PRJCA047587 and PRJCA048147.

Conflicts of interest statement

The authors declare no conflicts of interest.

Supplementary material

Supplementary material is available at Horticulture Research online.

References

[1] Zhang J, Nieminen K, Serra JAA. et al. The formation of wood and its control. Curr Opin Plant Biol. 2014; 17: 56-63

[2] Xu W, Cheng H, Zhu S. et al. Functional understanding of secondary cell wall cellulose synthases in Populus trichocarpa via the Cas9/gRNA-induced gene knockouts. New Phytol. 2021; 231: 1478-95

[3] Suzuki S, Li L, Sun Y-H. et al. The cellulose synthase gene superfamily and biochemical functions of xylem-specific cellulose synthase-like genes in Populus trichocarpa. Plant Physiol. 2006; 142: 1233-45

[4] Jiang P-F, Lin X-Y, Bian X-Y. et al. Ectopic expression of Populus MYB10 promotes secondary cell wall thickening and inhibits anthocyanin accumulation. Plant Physiol Biochem. 2022; 172: 24-32

[5] Tang F, Jiao B, Zhang M. et al. PtoMYB031, the R2R3 MYB transcription factor involved in secondary cell wall biosynthesis in poplar. Front Plant Sci. 2024; 14: 1341245

[6] Sun L, Li P, Cui X. et al. PtoNAC90 is a positive regulator for secondary xylem development in Populus. Int J Biol Macromol. 2025; 320: 145803

[7] Li C, Wang X, Ran L. et al. PtoMYB92 is a transcriptional activator of the lignin biosynthetic pathway during secondary cell wall formation in Populus tomentosa. Plant Cell Physiol. 2015; 56: 2436-46

[8] Kim M-H, Cho J-S, Bae E-K. et al. PtrMYB120 functions as a positive regulator of both anthocyanin and lignin biosynthetic pathway in a hybrid poplar. Tree Physiol. 2021; 41: 2409-23

[9] Zhang S, Li C, Cui B. et al. TCP23-WRKY15 module negatively regulates lignin deposition and xylem development of wood formation in Populus. Int J Biol Macromol. 2025; 306: 141656

[10] Quan M, Zhou J, Khan M. et al. The allelic regulatory network of PtobHLH125 regulates secondary cell wall formation by affecting lignin biosynthesis in Populus tomentosa. Ind Crop Prod. 2025; 233: 121476

[11] Sjöstrand B, Moser C, Jäppinen A. et al. Fiber length characteristics of Norway spruce (Picea abies) trees with fast height growth. Holzforschung. 2025; 79: 404-17

[12] Li M, Wu W, Mo Y. et al. Genome annotation, comparative genomics and transcriptomic analysis of Eucalyptus cloeziana reveal insights into genome evolution and wood formation in eucalyptus. Ind Crop Prod. 2025; 235: 121680

[13] Yang J, Luo Y, Bian L. et al. A Chinese fir ClWND2B-ClMYB3 transcriptional cascade promotes lignin accumulation during wood formation in poplar. Ind Crop Prod. 2025; 229: 121045

[14] Guo F, Yu W, Fu F. et al. Ginkgo biloba wood transcriptome reveals critical genes for secondary cell wall formation and transcription factors involved in lignin biosynthesis. Ind Crop Prod. 2024; 216: 118736

[15] Guo F, Yu W, Fu F. et al. Physiological, transcriptome and metabolome analyses provide molecular insights to seasonal development in Ginkgo biloba xylem. Ind Crop Prod. 2024; 208: 117930

[16] Gao T, Shang J, Guo J. et al. Whole-transcriptome sequencing reveals the regulatory roles of non-coding RNAs during xylem development in Ginkgo biloba L. Ind Crop Prod. 2025; 235: 121775

[17] Risch N, Merikangas K . The future of genetic studies of complex human diseases. Science. 1996; 273: 1516-7

[18] Klein RJ, Zeiss C, Chew EY. et al. Complement factor H polymorphism in age-related macular degeneration. Science. 2005; 308: 385-9

[19] Huang X, Wei X, Sang T. et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet. 2010; 42: 961-7

[20] Liu R, Gong J, Xiao X. et al. GWAS analysis and QTL identification of fiber quality traits and yield components in upland cotton using enriched high-density SNP markers. Front Plant Sci. 2018; 9: 1067

[21] Yu J, Hui Y, Chen J. et al. Whole-genome resequencing of 240 Gossypium barbadense accessions reveals genetic variation and genes associated with fiber strength and lint percentage. Theor Appl Genet. 2021; 134: 3249-61

[22] Zhou Y, Tang H, Cheng M-P. et al. Genome-wide association study for pre-harvest sprouting resistance in a large germplasm collection of Chinese wheat landraces. Front Plant Sci. 2017; 08: 08

[23] Zhang J, Yang Y, Zheng K. et al. Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus. New Phytol. 2018; 220: 502-16

[24] Fahrenkrog AM, Neves LG, Resende MFR. et al. Genome-wide association study reveals putative regulators of bioenergy traits in Populus deltoides. New Phytol. 2017; 213: 799-811

[25] Lamara M, Raherison E, Lenz P. et al. Genetic architecture of wood properties based on association analysis and co-expression networks in white spruce. New Phytol. 2016; 210: 240-55

[26] Cappa EP, El-Kassaby YA, Garcia MN. et al. Impacts of population structure and analytical models in genome-wide association studies of complex traits in forest trees: a case study in Eucalyptus globulus. PLoS One. 2013; 8: e81267

[27] Uchiyama K, Iwata H, Moriguchi Y. et al. Demonstration of genome-wide association studies for identifying markers for wood property and male strobili traits in Cryptomeria japonica. PLoS One. 2013; 8: e79866

[28] Gusev A, Ko A, Shi H. et al. Integrative approaches for large-scale transcriptome-wide association studies. Nat Genet. 2016; 48: 245-52

[29] Tang S, Zhao H, Lu S. et al. Genome- and transcriptome-wide association studies provide insights into the genetic basis of natural variation of seed oil content in Brassica napus. Mol Plant. 2021; 14: 470-87

[30] Li Z, Wang P, You C. et al. Combined GWAS and eQTL analysis uncovers a genetic regulatory network orchestrating the initiation of secondary cell wall development in cotton. New Phytol. 2020; 226: 1738-52

[31] Hu Y, Yu Z, Gao X. et al. Genetic diversity, population structure, and genome-wide association analysis of Ginkgo cultivars. Hortic Res. 2023; 10: uhad136

[32] Zhou Q, Shen X, Li Y . Variation pattern and genome-wide association study of leaf phenotypic traits among ancient Ginkgo biloba L. populations. Forests. 2022; 13: 1764

[33] Wang S, Li Z, Guo X. et al. Comparative analysis of microsatellite, SNP, and InDel markers in four Rhododendron species based on RNA-seq. Breed Sci. 2018; 68: 536-44

[34] Arbeithuber B, Betancourt AJ, Ebner T. et al. Crossovers are associated with mutation and biased gene conversion at recombination hotspots. Proc Natl Acad Sci USA. 2015; 112: 2109-14

[35] Shen J-C, Rideout WM, Jones PA . The rate of hydrolytic deamination of 5-methylcytosine in double-stranded DNA. Nucl Acids Res. 1994; 22: 972-6

[36] Neddermann P, Jiricny J . Efficient removal of uracil from G.U mispairs by the mismatch-specific thymine DNA glycosylase from HeLa cells. Proc Natl Acad Sci USA. 1994; 91: 1642-6

[37] Kim W-C, Kim J-Y, Ko J-H. et al. Transcription factor MYB46 is an obligate component of the transcriptional regulatory complex for functional expression of secondary wall-associated cellulose synthases in Arabidopsis thaliana. J Plant Physiol. 2013; 170: 1374-8

[38] Atanassov II, Pittman JK, Turner SR . Elucidating the mechanisms of assembly and subunit interaction of the cellulose synthase complex of Arabidopsis secondary cell walls. J Biol Chem. 2009; 284: 3833-41

[39] Nemhauser JL, Torii KU . Plant synthetic biology for molecular engineering of signalling and development. Nat Plants. 2016; 2: 16010

[40] Klose H, Günl M, Usadel B. et al. Cell wall modification in tobacco by differential targeting of recombinant endoglucanase from Trichoderma reesei. BMC Plant Biol. 2015; 15: 54

[41] Wang H, Avci U, Nakashima J. et al. Mutation of WRKY transcription factors initiates pith secondary wall formation and increases stem biomass in dicotyledonous plants. Proc Natl Acad Sci. 2010; 107: 22338-43

[42] Yu Y, Hu R, Wang H. et al. MlWRKY12, a novel Miscanthus transcription factor, participates in pith secondary cell wall formation and promotes flowering. Plant Sci. 2013; 212: 1-9

[43] Kim K, Kang A-H, Jung H. et al. The activation of the NAC transcription factor by poplar PtrGATA9 facilitates the secondary cell wall deposition of interfascicular fiber cell in Arabidopsis. Plant Biotechnol Rep. 2023; 17: 665-76

[44] Wei H, Song Z, Xie Y. et al. High temperature inhibits vascular development via the PIF4-miR166-HB15 module in Arabidopsis. Curr Biol. 2023; 33: 3203-3214.e4

[45] Abbas M, Peszlen I, Shi R. et al. Involvement of CesA4, CesA7-A/B and CesA8-A/B in secondary wall formation in Populus trichocarpa wood. Tree Physiol. 2020; 40: 73-89

[46] Davis J, Brandizzi F, Liepman AH. et al. Arabidopsis mannan synthase CSLA9 and glucan synthase CSLC4 have opposite orientations in the Golgi membrane: Hemicellulosic glycan synthase topology. Plant J. 2010; 64: 1028-37

[47] Zeng W, Lampugnani ER, Picard KL. et al. Asparagus IRX9, IRX10, and IRX14A are components of an active xylan backbone synthase complex that forms in the Golgi apparatus. Plant Physiol. 2016; 171: 93-109

[48] Du Q, Pan W, Tian J. et al. The UDP-glucuronate decarboxylase gene family in Populus: structure, expression, and association genetics. PLoS One. 2013; 8: e60880

[49] Botticella E, Savatin DV, Sestili F . The triple jags of dietary fibers in cereals: how biotechnology is longing for high Fiber-Grains. Front Plant Sci. 2021; 12: 745579

[50] Chong S-L, Virkki L, Maaheimo H. et al. O-acetylation of glucuronoxylan in Arabidopsis thaliana wild type and its change in xylan biosynthesis mutants. Glycobiology. 2014; 24: 494-506

[51] Urbanowicz BR, Peña MJ, Ratnaparkhe S. et al. 4-O-methylation of glucuronic acid in Arabidopsis glucuronoxylan is catalyzed by a domain of unknown function family 579 protein. Proc Natl Acad Sci USA. 2012; 109: 14253-8

[52] Rautengarten C, Ebert B, Moreno I. et al. The Golgi localized bifunctional UDP-rhamnose/UDP-galactose transporter family of Arabidopsis. Proc Natl Acad Sci USA. 2014; 111: 11563-8

[53] Zhang J, Xie M, Tuskan GA. et al. Recent advances in the transcriptional regulation of secondary cell wall biosynthesis in the woody plants. Front Plant Sci. 2018; 9: 1535

[54] Kim W-C, Reca I-B, Kim Y. et al. Transcription factors that directly regulate the expression of CSLA9 encoding mannan synthase in Arabidopsis thaliana. Plant Mol Biol. 2014; 84: 577-87

[55] Lynch JH, Dudareva N . Aromatic amino acids: a complex network ripe for future exploration. Trends Plant Sci. 2020; 25: 670-81

[56] Maeda H, Dudareva N . The shikimate pathway and aromatic amino acid biosynthesis in plants. Annu Rev Plant Biol. 2012; 63: 73-105

[57] Cao Y, Li X, Jiang L . Integrative analysis of the core fruit lignification toolbox in pear reveals targets for fruit quality bioengineering. Biomolecules. 2019; 9: 504

[58] Xie G, Feng Y, Chen Y. et al. Effects of 1-methylcyclopropene (1-MCP) and ethylene on postharvest lignification of common beans (Phaseolus vulgaris L). ACS Omega. 2020; 5: 8659-66

[59] Pratyusha DS, Sarada DVL . MYB transcription factors-master regulators of phenylpropanoid biosynthesis and diverse developmental and stress responses. Plant Cell Rep. 2022; 41: 2245-60

[60] Xu M, He S, Zhang H. et al. Gibberellin regulates the synthesis of stone cells in ‘Nanguo’ pear via the PUMYB91-PUERF023 module. Physiol Plant. 2025; 177: e70074

[61] Song C, Guo Y, Shen W. et al. PagUNE12 encodes a basic helix-loop-helix transcription factor that regulates the development of secondary vascular tissue in poplar. Plant Physiol. 2023; 192: 1046-62

[62] Gao Z, Sun W, Wang J. et al. GhbHLH18 negatively regulates fiber strength and length by enhancing lignin biosynthesis in cotton fibers. Plant Sci. 2019; 286: 7-16

[63] Chang J, Luo H, Guo T. et al. Transcriptome-based analysis identifies the key biosynthetic genes and regulators responsible for lignification in harvested Tsai Tai. Postharvest Biol Technol. 2024; 216: 113083

[64] Wang X, Chai X, Gao B. et al. Multi-omics analysis reveals the mechanism of bHLH130 responding to low-nitrogen stress of apple rootstock. Plant Physiol. 2023; 191: 1305-23

[65] Guo Y, Qi Y, Yang G. et al. A genome-wide identification of miPEPs in hybrid poplar reveals regulatory functions of miPEP166i in adventitious root elongation. Ind Crop Prod. 2024; 210: 118150

[66] Ding Z, Friml J . Auxin regulates distal stem cell differentiation in Arabidopsis roots. Proc Natl Acad Sci USA. 2010; 107: 12046-51

[67] Xu C, Chang P, Guo S. et al. Transcriptional activation by WRKY23 and depression by removal of bHLH041 coordinately establish callus pluripotency in Arabidopsis regeneration. Plant Cell. 2023; 36: 158-73

[68] Aida M, Beis D, Heidstra R. et al. The PLETHORA genes mediate patterning of the Arabidopsis root stem cell niche. Cell. 2004; 119: 109-20

[69] Galinha C, Hofhuis H, Luijten M. et al. PLETHORA proteins as dose-dependent master regulators of Arabidopsis root development. Nature. 2007; 449: 1053-7

[70] Minic Z, Rihouey C, Do CT. et al. Purification and characterization of enzymes exhibiting β-D-Xylosidase activities in stem tissues of Arabidopsis. Plant Physiol. 2004; 135: 867-78

[71] Foreman J, Demidchik V, Bothwell JHF. et al. Reactive oxygen species produced by NADPH oxidase regulate plant cell growth. Nature. 2003; 422: 442-6

[72] Kwak JM, Mori IC, Pei ZM. et al. NADPH oxidase AtrbohD and AtrbohF genes function in ROS-dependent ABA signaling in Arabidopsis. EMBO J. 2003; 22: 2623-33

[73] Qin Y-M, Hu C-Y, Zhu Y-X . The ascorbate peroxidase regulated by H₂O₂ and ethylene is involved in cotton fiber cell elongation by modulating ROS homeostasis . Plant Signal Behav. 2008; 3: 194-6

[74] Chen X, Guo W, Liu B. et al. Molecular mechanisms of fiber differential development between G. barbadense and G. hirsutum revealed by genetical genomics. PLoS One. 2012; 7: e30056

[75] Keyzor C, Mermaz B, Trigazis E. et al. Histone demethylases ELF6 and JMJ13 antagonistically regulate self-fertility in Arabidopsis. Front Plant Sci. 2021; 12: 640135

[76] Franklin GL . Preparation of thin sections of synthetic resins and wood-resin composites, and a new macerating method for wood. Nature. 1945; 155: 51-1

[77] Hu K, Huang Y, Fei B. et al. Investigation of the multilayered structure and microfibril angle of different types of bamboo cell walls at the micro/nano level using a LC-PolScope imaging system. Cellulose. 2017; 24: 4611-25

[78] Segal L, Creely JJ, Martin AE. et al. An empirical method for estimating the degree of crystallinity of native cellulose using the X-ray diffractometer. Text Res J. 1959; 29: 786-94

[79] Liu H, Wang X, Wang G. et al. The nearly complete genome of Ginkgo biloba illuminates gymnosperm evolution. Nat Plants. 2021; 7: 748-56

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

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

[82] Van der Auwera GA, O’Connor BD . Genomics in the Cloud: Using Docker, GATK, and WDL in Terra. 1st ed. Sebastopol, CA: O’Reilly Media; 2020:

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

[84] Kim D, Langmead B, Salzberg SL . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60

[85] Putri GH, Anders S, Pyl PT. et al. Analysing high-throughput sequencing data in Python with HTSeq 2.0. Bioinformatics. 2022; 38: 2943-5

[86] Liang M, An B, Deng T. et al. Incorporating genome-wide and transcriptome-wide association studies to identify genetic elements of longissimus dorsi muscle in Huaxi cattle. Front Genet. 2023; 13: 982433

[87] Shannon P, Markiel A, Ozier O. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13: 2498-504