Chromosome-level genome assembly and population genomics unveil strigolactone-regulated growth adaptation in the mycoheterotrophic orchid Gastrodia elata

Abstract

Mycoheterotrophic plants rely entirely on fungal symbionts for nutrients, yet the role of intraspecific genomic variation in shaping symbiotic adaptation remains unclear. Gastrodia elata is a mycoheterotrophic orchid with multiple cultivated varieties. Here, we generated a chromosome-level genome of G. elata Bl. f. glauca. Comparative genomic analyses with published G. elata assemblies revealed extensive intraspecific variation, characterized by transposon-mediated inversions occurring in 26% of syntenic regions. Notably, these regions frequently harbored orphan genes. Population genomic analysis of 150 individuals identified three genetically distinct clades: two cultivated (Clades E and G) and one hybrid (Clade I). Transcriptomic profiling uncovered clade-specific expression patterns in symbiosis-related genes, particularly within strigolactone signaling pathways. Molecular dynamics simulations and protein interaction assays demonstrated that polymorphisms in the M domain of the suppressor protein DWARF53 (GeD53) modulate strigolactone signaling by altering the stability of its interaction with the receptor (GeD14). Specifically, a Clade G-specific haplotype enhanced signaling through stabilized protein interactions, thereby influencing tuber development genes, whereas GeD14 variants had minimal functional impact. Further co-expression networks identified LOL5 , RNP1 , and MTHD as downstream effectors correlating with clade-specific tuber phenotypes and carbohydrate allocation. These findings demonstrate how intraspecific variation in strigolactone signaling components drives functional divergence in G. elata , providing both mechanistic insights into mycoheterotrophic adaptation and genomic resources for future research.

References

[1] Jhu MY, Sinha NR . Parasitic plants: an overview of mechanisms by which plants perceive and respond to parasites. Annu Rev Plant Biol. 2022; 73: 433-55

[2] Miura C, Furui Y, Yamamoto T. et al. Autoactivation of mycorrhizal symbiosis signaling through gibberellin deactivation in orchid seed germination. Plant Physiol. 2024; 194: 546-63

[3] Gomes SIF, Merckx VSFT, Kehl J. et al. Mycoheterotrophic plants living on arbuscular mycorrhizal fungi are generally enriched in 13C, 15N and 2H isotopes. J Ecol. 2020; 108: 1250-61

[4] Tešitel J, Li AR, Knotková K. et al. The bright side of parasitic plants: what are they good for? Plant Physiol. 2021; 185: 1309-24

[5] Vogel A, Schwacke R, Denton AK. et al. Footprints of parasitism in the genome of the parasitic flowering plant Cuscuta campestris. Nat Commun. 2018;9:2515

[6] Sun G, Xu Y, Liu H. et al. Large-scale gene losses underlie the genome evolution of parasitic plant Cuscuta australis. Nat Commun. 2018; 9: 2683

[7] Cai L, Arnold BJ, Xi Z. et al. Deeply altered genome architecture in the endoparasitic flowering plant Sapria himalayana Griff. (Rafflesiaceae). Curr Biol. 2021; 31: 1002-1011.e9

[8] Yuan Y, Jin X, Liu J. et al. The Gastrodia elata genome provides insights into plant adaptation to heterotrophy. Nat Commun. 2018; 9: 1615

[9] Li MH, Liu KW, Li Z. et al. Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy. Nat Plants. 2022; 8: 373-88

[10] Aquino B, Bradley JM, Lumba S . On the outside looking in: roles of endogenous and exogenous strigolactones. Plant J. 2021; 105: 322-34

[11] Hua Z, Teng X, Huang J. et al. The Armillaria response to Gastrodia elata is partially mediated by strigolactone-induced changes in reactive oxygen species. Microbiol Res. 2024; 278: 127536

[12] Wang Y, Shahid MQ . Genome sequencing and resequencing identified three horizontal gene transfers and uncovered the genetic mechanism on the intraspecies adaptive evolution of Gastrodia elata Blume. Front Plant Sci. 2023; 13: 1035157

[13] Bae EK, An C, Kang MJ. et al. Chromosome-level genome assembly of the fully mycoheterotrophic orchid Gastrodia elata. G3 (Bethesda). 2022; 12: jkab433

[14] Xu Y, Lei Y, Su Z. et al. A chromosome-scale Gastrodia elata genome and large-scale comparative genomic analysis indicate convergent evolution by gene loss in mycoheterotrophic and parasitic plants. Plant J. 2021; 108: 1609-23

[15] Zhou X, Chen X . The taxonomy of domestic Gastrodia species. Acta Bot Yunnanica. 1983; 5: 361-8

[16] Song B, Ning W, Wei D. et al. Plant genome resequencing and population genomics: current status and future prospects. Mol Plant. 2023; 16: 1252-68

[17] Alves MF, Pinheiro F, Nunes CEP. et al. Reproductive development and genetic structure of the mycoheterotrophic orchid Pogoniopsis schenckii Cogn. BMC Plant Biol. 2021; 21: 332

[18] Minasiewicz J, Krawczyk E, Znaniecka J. et al. Weak population spatial genetic structure and low infraspecific specificity for fungal partners in the rare mycoheterotrophic orchid Epipogium aphyllum. J Plant Res. 2022; 135: 275-93

[19] Chen S, Wang X, Wang Y. et al. Improved de novo Assembly of the Achlorophyllous Orchid Gastrodia elata. Front. Genet. 2020; 11: 580568

[20] Wang Y, Shahid MQ . Genome sequencing and resequencing identified three horizontal gene transfers and uncovered the genetic mechanism on the intraspecies adaptive evolution of Gastrodia elata Blume. Front. Plant Sci. 2023; 13: 1035157

[21] Zhang L, Wu S, Chang X. et al. The ancient wave of polyploidization events in flowering plants and their facilitated adaptation to environmental stress. Plant Cell Environ. 2020; 43: 2847-56

[22] Zhang GQ, Liu KW, Li Z. et al. The Apostasia genome and the evolution of orchids. Nature. 2017; 549: 379-83

[23] Harkess A, Zhou J, Xu C. et al. The Asparagus genome sheds light on the origin and evolution of a young Y chromosome. Nat Commun. 2017; 8: 1279

[24] Mashiguchi K, Seto Y, Yamaguchi S . Strigolactone biosynthesis, transport and perception. Plant J. 2021; 105: 335-50

[25] Zhou F, Lin Q, Zhu L. et al. D14-SCFD3-dependent degradation of D53 regulates strigolactone signalling. Nature. 2013; 504: 406-10

[26] Xu C, He C . The rice OsLOL2 gene encodes a zinc finger protein involved in rice growth and disease resistance. Mol Gen Genomics. 2007; 278: 85-94

[27] Arite T, Umehara M, Ishikawa S. et al. d14, a Strigolactone-insensitive mutant of rice, shows an accelerated outgrowth of tillers. Plant Cell Physiol. 2009; 50: 1416-24

[28] Wang Z, Zhao X, Wang B. et al. Overexpression of an Arabidopsis heterogeneous nuclear ribonucleoprotein gene, AtRNP1, affects plant growth and reduces plant tolerance to drought and salt stresses. Biochem Biophys Res Commun. 2016; 472: 353-9

[29] Guan X, Okazaki Y, Lithio A. et al. Discovery and characterization of the 3-hydroxyacyl-ACP dehydratase component of the plant mitochondrial fatty acid synthase system. Plant Physiol. 2017; 173: 2010-28

[30] Xiang D, Yang H, Venglat P. et al. POPCORN functions in the auxin pathway to regulate embryonic body plan and meristem organization in Arabidopsis. Plant Cell. 2011; 23: 4348-67

[31] Hung SH, Elliott GI, Ramkumar TR. et al. Structural basis of Qng1-mediated salvage of the micronutrient queuine from queuosine-5'-monophosphate as the biological substrate. Nucleic Acids Res. 2023; 51: 935-51

[32] Guo X, Hu X, Li J. et al. The Sapria himalayana genome provides new insights into the lifestyle of endoparasitic plants. BMC Biol. 2023; 21: 134

[33] Merckx V, Freudenstein JV . Evolution of mycoheterotrophy in plants: a phylogenetic perspective. New Phytol. 2010; 185: 605-9

[34] Wong DCJ, Wang Z, Perkins J. et al. The road less taken: Dihydroflavonol 4-reductase inactivation and delphinidin anthocyanin loss underpins a natural intraspecific flower colour variation. Mol Ecol. 2024; 34: e17334

[35] Bennetzen JL, Wang H . The contributions of transposable elements to the structure, function, and evolution of plant genomes. Annu Rev Plant Biol. 2014; 65: 505-30

[36] Yue Y, Yu R, Fan Y . Characterization of two monoterpene synthases involved in floral scent formation in Hedychium coronarium. Planta. 2014; 240: 745-62

[37] Zhang Y, Liu T, Nan T. et al. Characteristics and functions of volatile organic compounds in the tripartite symbiotic system of Gastrodia elata-Armillaria gallica-Rahnella aceris HPDA25. Plant Signal Behav. 2024; 19: 2399426

[38] Gutjahr C, Gobbato E, Choi J. et al. Rice perception of symbiotic arbuscular mycorrhizal fungi requires the karrikin receptor complex. Science. 2015; 350: 1521-4

[39] Zhao LH, Zhou XE, Yi W. et al. Destabilization of strigolactone receptor DWARF14 by binding of ligand and E3-ligase signaling effector DWARF3. Cell Res. 2015; 25: 1219-36

[40] Liu S, Wang J, Song B. et al. Conformational dynamics of the D53-D3-D14 complex in strigolactone signaling. Plant Cell Physiol. 2023; 64: 1046-56

[41] Wang L, Wang B, Yu H. et al. Transcriptional regulation of strigolactone signalling in Arabidopsis. Nature. 2020; 583: 277-81

[42] Jiang L, Liu X, Xiong G. et al. DWARF 53 acts as a repressor of strigolactone signalling in rice. Nature. 2013; 504: 401-5

[43] Papaleo E, Saladino G, Lambrughi M. et al. The role of protein loops and linkers in conformational dynamics and allostery. Chem Rev. 2016; 116: 6391-423

[44] Chen A, Martin G, Viljoen A. et al. Identification of a QTL conferring resistance to the subtropical race 4 of Fusarium oxysporum f. sp. cubense in Calcutta 4 (Musa acuminata ssp. burmannica). Hortic Res. 2026; uhag001

[45] Kadium VR, Pilli R, Svyantek A. et al. Dissecting the genetic basis of relevant fruit quality traits in interspecific grapevines (Vitis spp.). Hortic Res. 2026; 13: uhaf353

[46] Zhan HD, Zhou HY, Sui YP. et al. The rhizome of Gastrodia elata Blume-an ethnopharmacological review. J Ethnopharmacol. 2016; 189: 361-85

[47] Ye H, Hou Q, Lv H. et al. D53 represses rice blast resistance by directly targeting phenylalanine ammonia lyases. J Integr Plant Biol. 2024; 66: 1827-30

[48] Rogers SO, Bendich AJ . Extraction of total cellular DNA from plants, algae and fungi. In: Gelvin SB, Schilperoort RA, eds. Plant Molecular Biology Manual.Springer Netherlands: Dordrecht, 1994, 183-90

[49] Chin CS, Peluso P, Sedlazeck FJ. et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat Methods. 2016; 13: 1050-4

[50] Walker BJ, Abeel T, Shea T. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014; 9: e112963

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

[52] Dudchenko O, Batra SS, Omer AD. et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science. 2017; 356: 92-5

[53] Durand NC, Robinson JT, Shamim MS. et al. Juicebox provides a visualization system for Hi-C contact maps with unlimited zoom. Cell Syst. 2016; 3: 99-101

[54] Ou S, Su W, Liao Y. et al. Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline. Genome Biol. 2019; 20: 275

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

[56] Smit AF, Hubley R, Green P . RepeatMasker Open-4.0. 2013-2015. http://www.repeatmasker.org.

[57] Kim D, Paggi JM, Park C. et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37: 907-15

[58] Pertea M, Pertea GM, Antonescu CM. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5

[59] Grabherr MG, Haas BJ, Yassour M. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29: 644-52

[60] Haas BJ, Delcher AL, Mount SM. et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003; 31: 5654-66

[61] Brůna T, Hoff KJ, Lomsadze A. et al. BRAKER2: automatic eukaryotic genome annotation with GeneMark-EP+ and AUGUSTUS supported by a protein database. NAR Genom Bioinform. 2021; 3: lqaa108

[62] Lomsadze A, Burns PD, Borodovsky M . Integration of mapped RNA-Seq reads into automatic training of eukaryotic gene finding algorithm. Nucleic Acids Res. 2014; 42: e119-9

[63] She R, Chu JSC, Wang K. et al. genBlastA: enabling BLAST to identify homologous gene sequences. Genome Res. 2009; 19: 143-9

[64] Birney E, Clamp M, Durbin R . GeneWise and Genomewise. Genome Res. 2004; 14: 988-95

[65] Keilwagen J, Wenk M, Erickson JL. et al. Using intron position conservation for homology-based gene prediction. Nucleic Acids Res. 2016; 44: e89-9

[66] Haas BJ, Salzberg SL, Zhu W. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 2008; 9: R7

[67] Li H. Protein-to-genome alignment with miniprot. Bioinformatics. 2023; 39: btad014

[68] Goel M, Sun H, Jiao WB. et al. SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies. Genome Biol. 2019; 20: 277

[69] Li L, Stoeckert CJ, Roos DS . OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003; 13: 2178-89

[70] Ranwez V, Douzery EJP, Cambon C. et al. MACSE v2: toolkit for the alignment of coding sequences accounting for frameshifts and stop codons. Mol Biol Evol. 2018; 35: 2582-4

[71] Yang Z . PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007; 24: 1586-91

[72] Sun P, Jiao B, Yang Y. et al. WGDI: a user-friendly toolkit for evolutionary analyses of whole-genome duplications and ancestral karyotypes. Mol Plant. 2022; 15: 1841-51

[73] Mendes FK, Vanderpool D, Fulton B. et al. CAFE 5 models variation in evolutionary rates among gene families. Bioinformatics. 2020; 36: 5516-8

[74] Wu T, Hu E, Xu S. et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation. 2021; 2: 100141

[75] Chen S, Zhou Y, Chen Y. et al. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018; 34: i884-90

[76] McKenna A, Hanna M, Banks E. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303

[77] Chang CC, Chow CC, Tellier LC. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015; 4: 7

[78] Manichaikul A, Mychaleckyj JC, Rich SS. et al. Robust relationship inference in genome-wide association studies. Bioinformatics. 2010; 26: 2867-73

[79] 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

[80] Xu L, He W, Tai S. et al. VCF2Dis: an ultra-fast and efficient tool to calculate pairwise genetic distance and construct population phylogeny from VCF files. GigaScience. 2025; 14: giaf032

[81] 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

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

[83] Korunes KL, Samuk K . pixy: unbiased estimation of nucleotide diversity and divergence in the presence of missing data. Mol Ecol Resour. 2021; 21: 1359-68

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

[85] Dobin A, Davis CA, Schlesinger F. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013; 29: 15-21

[86] Anders S, Pyl PT, Huber W . HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015; 31: 166-9

[87] Love MI, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550

[88] Chen H, Zou Y, Shang Y. et al. Firefly luciferase complementation imaging assay for protein-protein interactions in plants. Plant Physiol. 2008; 146: 323-4

[89] Abramson J, Adler J, Dunger J. et al. Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature. 2024; 630: 493-500

[90] Case DA, Aktulga HM, Belfon K. et al. AmberTools. J Chem Inf Model. 2023; 63: 6183-91