Genome-Wide Identification, Evolution, and Expression Analysis of LBD Transcription Factor Family in Bread Wheat (Triticum aestivum L.)

Genome-wide identification and functional characterization of LBD gene family in four Cymbidium species (Orchidaceae) and potential regulatory role of CsiLBD27 in floral development of Cymbidium sinense

BACKGROUND

The Lateral Organ Boundaries Domain (LBD) gene family, encompassing plant-specific LOB domain proteins, plays essential roles in various aspects of plant growth and development, and has continuously diversified its functions across numerous species. However, studies on LBD genes in the Orchidaceae family remain limited. To our knowledge, this is the first systematic investigation of the LBD gene family in Cymbidium, a genus that exhibits remarkable species and trait diversity within the Orchidaceae.

RESULTS

In this study, we identified 122 LBD genes within the genomes of four Cymbidium species, distributed on 20 chromosomes. These genes were classified into class I (109 members, including 16 in subclass Ic/d) and class II (13 members), with protein lengths ranging from 94-477 amino acids. Promoter sequences of CymLBD genes revealed various cis-elements significant for light, hormonal, biotic, and abiotic stress responses. Transcriptomic analysis revealed tissue- and stage-specific expression of CymLBD genes in Cymbidium, and RT-qPCR and yeast one-hybrid assays indicated that CsiLBD27 may regulate floral patterning by directly binding to the CsiSEP3 promoter and activating its transcription. Under ABA treatment, the genes CsiLBD13, CsiLBD19, and CsiLBD21 displayed tissue-specific expression changes, suggesting hormone-responsive regulation. In Cymbidium ensifolium, 16 CenLBD genes were differentially expressed, while in Cymbidium mannii, 9 CmaLBD gene expression exist obvious circadian rhythm. GO and KEGG enrichment of 1074 and 399 predicted target genes, respectively, indicating a diverse range of functions for LBDs.

CONCLUSIONS

Collectively, this study provides the first comprehensive insight into the evolutionary dynamics, regulatory mechanisms, and functional roles of LBD genes in Cymbidium. These findings offer a valuable genetic resource for understanding floral and vegetative development in Orchidaceae and uncover potential novel functions of LBD genes.

Comprehensive Characterization and Functional Analysis of the Lateral Organ Boundaries Domain Gene Family in Rice: Evolution, Expression, and Stress Response

In this study, the LBD (Lateral Organ Boundaries Domain) gene family, a group of plant-specific transcription factors critical for plant growth and development as well as metabolic regulation, was comprehensively characterized in rice. We identified 36 LBD genes using multi-source genomic data and systematically classified them into Class I (31 genes) and Class II (5 genes). Analysis of their physicochemical properties revealed significant variations in amino acid length, molecular weight, isoelectric points, and hydropathicity. Motif analysis identified conserved LOB domains and other motifs potentially linked to functional diversity. Cis-acting element analysis indicated the involvement of these genes in various biological processes, including light response, hormone signaling, and stress response. Expression profiling demonstrated tissue-specific expression patterns, with several genes, such as XM_015770711.2, XM_015776632.2, and XM_015792766.2, showing relatively high expression in rice roots, implying their important role in root development. Transcriptome data further supported the involvement of specific genes in responses to phytohormones such as jasmonic acid (JA) and abscisic acid (ABA), as well as environmental stresses like cold and drought. Notably, XM_015770711.2, XM_015776632.2, and XM_015772758.2 may contribute to the regulation of rice environmental adaptability by mediating ABA and JA signaling pathways, respectively. In conclusion, this study identified members of the LBD gene family through the screening of two rice gene databases, and performed a comprehensive analysis of their physicochemical properties, evolutionary relationships, and expression profiles under various conditions. These findings provided valuable insights for further functional studies of LBD genes. Moreover, this study provides a foundation for targeting LBD genes to enhance stress resilience (e.g., drought/cold tolerance) and root architecture optimization. The LBD gene family possesses dual values in both stress resistance regulation and developmental optimization. The construction of its multidimensional functional map lays the theoretical and resource foundation for the precise design of high-yield and stress-resistant varieties.

Gibberellin Regulates LBD38-1 Responses to Xanthomonas arboricola pv. juglandis Infection in Walnut Bacterial Blight Pathogenesis

BACKGROUND

Plant responses to biotic and abiotic stresses are complex processes. Previous studies have shown that the LBD gene family plays important roles in plant growth and development as well as in plant defense against biotic and abiotic stresses. The expression of LBD genes was investigated in walnuts under biotic and abiotic stresses, revealing that LBD38-1 may be a key gene in the plant stress response. This study provides new insights into the roles of LBD genes in plant responses to biotic stress.

RESULTS

Forty-nine members of the JrLBD gene family were identified in the walnut genome and classified into six subfamilies. Comparative homology analysis through phylogenetic trees revealed that the presence of Group I-a and Group VI plays an important role in resistance to stressors. The expression of walnut LBD genes under cold-temperature, high-temperature, mechanical damage, and biotic stresses was analyzed via transcriptome sequencing, and the expression of JrLBD38-1 in the Group VI subfamily was particularly prominent. According to transcriptome profile analysis, JrLBD38-1 is highly expressed in different tissues of walnuts, suggesting that it plays a regulatory role in the growth and development of different tissues. The function of the Gibberellin (GA) response element in the JrLBD38-1 promoter was further analyzed and verified. These findings confirmed that GA regulated JrLBD38-1 expression changes during Xanthomonas arboricola pv. juglandis infestation of walnut leaves.

CONCLUSION

Forty-nine walnut JrLBDs were identified and classified into six subfamilies. JrLBD38-1 has GA-inducible expression, is regulated by GA under pathogenic bacterial stress, and is involved in the response to biotic stress. This function of JrLBD38-1 provides new insights into walnut disease resistance mechanisms.

Genome-Wide Analysis Elucidates the Roles of AhLBD Genes in Different Abiotic Stresses and Growth and Development Stages in the Peanut (Arachis hypogea L.)

The lateral organ boundaries domain (LBD) genes, as the plant-specific transcription factor family, play a crucial role in controlling plant architecture and stress tolerance. However, the functions of AhLBD genes in the peanut plant (Arachis hypogea L.) remain unclear. In this study, 73 AhLBDs were identified in the peanut plant and divided into three groups by phylogenetic tree analysis. Gene structure and conserved protein motif analysis supported the evolutionary conservation of AhLBDs. Tandem and segment duplications contributed to the expansion of AhLBDs. The evolutionary relationship analysis of LBD gene family between A. hypogaea and four other species indicated that the peanut plant had a close relationship with the soybean plant. AhLBDs played a very important role in response to growth and development as well as abiotic stress. Furthermore, gene expression profiling and real-time quantitative qRT-PCR analysis showed that AhLBD16, AhLBD33, AhLBD67, and AhLBD72 were candidate genes for salt stress, while AhLBD24, AhLBD33, AhLBD35, AhLBD52, AhLBD67, and AhLBD71 were candidate genes for drought stress. Our subcellular localization experiment revealed that AhLBD24, AhLBD33, AhLBD67, and AhLBD71 were located in the nucleus. Heterologous overexpression of AhLBD33 and AhLBD67 in Arabidopsis significantly enhanced tolerance to salt stress. Our results provide a theoretical basis and candidate genes for studying the molecular mechanism for abiotic stress in the peanut plant.

Identification and Characterization of LBD Gene Family in Pseudoroegneria libanotica Reveals Functions of PseLBD1 and PseLBD12 in Response to Abiotic Stress

The lateral organ boundaries domain (LBD) plays a vital role as a transcriptional coactivator within plants, serving as an indispensable function in growth, development, and stress response. In a previous study, we found that the LBD genes of Pseudoroegneria libanotica (a maternal donor for three-quarter of perennial Triticeae species with good stress resistance, holds great significance in exploring its response mechanisms to abiotic stress for the Triticeae tribe) might be involved in responding to drought stress. Therefore, we further identified the LBD gene family in this study. A total of 29 PseLBDs were identified. Among them, 24 were categorized into subclass I, while 5 fell into subclass II. The identification of cis-acting elements reveals the extensive involvement of PseLBDs in various biological processes in P. libanotica. Collinearity analysis indicates that 86% of PseLBDs were single-copy genes and have undergone a single whole-genome duplication event. Transcriptomic differential expression analysis of PseLBDs under drought stress reveals that the most likely candidates for responding to abiotic stress were PseLBD1 and PseLBD12. They have been demonstrated to respond to drought, salt, heavy metal, and heat stress in yeast. Furthermore, it is plausible that functional divergence might have occurred among their orthologous genes in wheat. This study not only establishes a foundation for a deeper understanding of the biological roles of PseLBDs in P. libanotica but also unveils novel potential genes for enhancing the genetic background of crops within Triticeae crops, such as wheat.

Characterization and expression profiles of the ZmLBD gene family in Zea mays

BACKGROUND

The Lateral Organ Boundaries Domain (LBD) gene family is a family of plant-specific transcription factors (TFs) that are widely involved in processes such as lateral organ formation, stress response, and nutrient metabolism. However, the function of LBD genes in maize remains poorly understood.

METHODS AND RESULTS

In this study, a total of 49 ZmLBD genes were identified at the genome-wide level of maize, they were classified into nine branches based on phylogenetic relationships, and all of them were predicted to be nuclear localized. The 49 ZmLBD genes formed eight pairs of segmental duplicates, and members of the same branches' members had similar gene structure and conserved motif composition. The promoters of ZmLBD genes contain multiple types of cis-acting elements. In addition, by constructing the regulatory network of ZmLBD and other genes and miRNAs, 12 and 22 ZmLBDs were found to be involved in the gene regulatory network and miRNA regulatory network, respectively. The expression pattern analysis suggests that ZmLBD genes may be involved in different biological pathways, and drought stress induced the expressions of two inbred lines.

CONCLUSIONS

The findings enhance our comprehension of the potential roles of the ZmLBD gene family in maize growth and development, which is pivotal for genetic enhancement and breeding efforts pertaining to this significant crop.

Genome-Wide Analysis of the Lateral Organ Boundaries Domain (LBD) Gene Family in Sweet Potato (Ipomoea batatas)

The LBD family is a plant-specific transcription factor family that plays an important role in a variety of biological processes. However, the function of IbLBD genes in sweet potato remains unclear. In this study, we identified a total of 53 IbLBD genes in sweet potato. Genetic structure showed that most of the IbLBD genes contained only two exons. Following the phylogenetic investigation, the IbLBD gene family was separated into Class I (45 members) and Class II (8) members. Both classes of proteins contained relatively conservative Motif1 and Motif2 domains. The chromosomal locations, gene duplications, promoters, PPI network, and GO annotation of the sweet potato LBD genes were also investigated. Furthermore, gene expression profiling and real-time quantitative PCR analysis showed that the expression of 12 IbLBD genes altered in six separate tissues and under various abiotic stresses. The IbLBD genes belonging to Class I were mostly expressed in the primary root, the pencil root, and the leaves of sweet potatoes, while the genes belonging to Class II were primarily expressed in the various sweet potato roots. The IbLBD genes belonging to Class I were mostly expressed in the primary root, the pencil root, and the leaves of sweet potatoes, while the genes belonging to Class II were primarily expressed in the fibrous root, pencil root, and tuber root.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Genome-Wide Identification of LBD Genes in Foxtail Millet (Setaria italica) and Functional Characterization of SiLBD21

Plant-specific lateral organ boundaries domain (LBD) proteins play important roles in plant growth and development. Foxtail millet (Setaria italica) is one new C₄ model crop. However, the functions of foxtail millet LBD genes are unknown. In this study, a genome-wide identification of foxtail millet LBD genes and a systematical analysis were conducted. A total of 33 SiLBD genes were identified. They are unevenly distributed on nine chromosomes. Among these SiLBD genes, six segmental duplication pairs were detected. The thirty-three encoded SiLBD proteins could be classified into two classes and seven clades. Members in the same clade have similar gene structure and motif composition. Forty-seven kinds of cis-elements were found in the putative promoters, and they are related to development/growth, hormone, and abiotic stress response, respectively. Meanwhile, the expression pattern was investigated. Most SiLBD genes are expressed in different tissues, while several genes are mainly expressed in one or two kinds of tissues. In addition, most SiLBD genes respond to different abiotic stresses. Furthermore, the function of SiLBD21, which is mainly expressed in roots, was characterized by ectopic expression in Arabidopsis and rice. Compared to controls, transgenic plants generated shorter primary roots and more lateral roots, indicating the function of SiLBD21 in root development. Overall, our study laid the foundation for further functional elucidation of SiLBD genes.

Genome-Wide Association Study for Grain Protein, Thousand Kernel Weight, and Normalized Difference Vegetation Index in Bread Wheat (Triticum aestivum L.)

Genomic regions governing grain protein content (GPC), 1000 kernel weight (TKW), and normalized difference vegetation index (NDVI) were studied in a set of 280 bread wheat genotypes. The genome-wide association (GWAS) panel was genotyped using a 35K Axiom array and phenotyped in three environments. A total of 26 marker-trait associations (MTAs) were detected on 18 chromosomes covering the A, B, and D subgenomes of bread wheat. The GPC showed the maximum MTAs (16), followed by NDVI (6), and TKW (4). A maximum of 10 MTAs was located on the B subgenome, whereas, 8 MTAs each were mapped on the A and D subgenomes. In silico analysis suggest that the SNPs were located on important putative candidate genes such as NAC domain superfamily, zinc finger RING-H2-type, aspartic peptidase domain, folylpolyglutamate synthase, serine/threonine-protein kinase LRK10, pentatricopeptide repeat, protein kinase-like domain superfamily, cytochrome P450, and expansin. These candidate genes were found to have different roles including regulation of stress tolerance, nutrient remobilization, protein accumulation, nitrogen utilization, photosynthesis, grain filling, mitochondrial function, and kernel development. The effects of newly identified MTAs will be validated in different genetic backgrounds for further utilization in marker-aided breeding.

Genome-Wide Analysis of the LATERAL ORGAN BOUNDARIES Domain (LBD) Members in Alfalfa and the Involvement of MsLBD48 in Nitrogen Assimilation

The LATERAL ORGAN BOUNDARIES DOMAIN (LBD) proteins, a transcription factor family specific to the land plants, have been implicated in multiple biological processes including organ development, pathogen response and the uptake of inorganic nitrogen. The study focused on LBDs in legume forage Alfalfa. The genome-wide analysis revealed that in Alfalfa 178 loci across 31 allelic chromosomes encoded 48 unique LBDs (MsLBDs), and the genome of its diploid progenitor M. sativa spp. Caerulea encoded 46 LBDs. Synteny analysis indicated that the expansion of AlfalfaLBDs was attributed to the whole genome duplication event. The MsLBDs were divided into two major phylogenetic classes, and the LOB domain of the Class I members was highly conserved relative to that of the Class II. The transcriptomic data demonstrated that 87.5% of MsLBDs were expressed in at least one of the six test tissues, and Class II members were preferentially expressed in nodules. Moreover, the expression of Class II LBDs in roots was upregulated by the treatment of inorganic nitrogen such as KNO3 and NH4Cl (0.3 mM). The overexpression of MsLBD48, a Class II member, in Arabidopsis resulted in growth retardance with significantly declined biomass compared with the non-transgenic plants, and the transcription level of the genes involved in nitrogen uptake or assimilation, including NRT1.1, NRT2.1, NIA1 and NIA2 was repressed. Therefore, the LBDs in Alfalfa are highly conserved with their orthologs in embryophytes. Our observations that ectopic expression of MsLBD48 inhibited Arabidopsis growth by repressing nitrogen adaption suggest the negative role of the transcription factor in plant uptake of inorganic nitrogen. The findings imply the potential application of MsLBD48 in Alfalfa yield improvement via gene editing.

Mapping quantitative trait loci and predicting candidate genes for Striga resistance in maize using resistance donor line derived from Zea diploperennis

The parasitic weed, Striga is a major biological constraint to cereal production in sub-Saharan Africa (SSA) and threatens food and nutrition security. Two hundred and twenty-three (223) F_2:3 mapping population involving individuals derived from TZdEI 352 x TZEI 916 were phenotyped for four Striga-adaptive traits and genotyped using the Diversity Arrays Technology (DArT) to determine the genomic regions responsible for Striga resistance in maize. After removing distorted SNP markers, a genetic linkage map was constructed using 1,918 DArTseq markers which covered 2092.1 cM. Using the inclusive composite interval mapping method in IciMapping, twenty-three QTLs influencing Striga resistance traits were identified across four Striga-infested environments with five stable QTLs (qGY4, qSC2.1, qSC2.2, qSC5, and qSC6) detected in more than one environment. The variations explained by the QTLs ranged from 4.1% (qSD2.3) to 14.4% (qSC7.1). Six QTLs each with significant additive × environment interactions were also identified for grain yield and Striga damage. Gene annotation revealed candidate genes underlying the QTLs, including the gene models GRMZM2G077002 and GRMZM2G404973 which encode the GATA transcription factors, GRMZM2G178998 and GRMZM2G134073 encoding the NAC transcription factors, GRMZM2G053868 and GRMZM2G157068 which encode the nitrate transporter protein and GRMZM2G371033 encoding the SBP-transcription factor. These candidate genes play crucial roles in plant growth and developmental processes and defense functions. This study provides further insights into the genetic mechanisms of resistance to Striga parasitism in maize. The QTL detected in more than one environment would be useful for further fine-mapping and marker-assisted selection for the development of Striga resistant and high-yielding maize cultivars.

Microspore embryogenesis induction by mannitol and TSA results in a complex regulation of epigenetic dynamics and gene expression in bread wheat

Reprogramming of microspores development towards embryogenesis mediated by stress treatment constitutes the basis of doubled haploid production. Recently, compounds that alter histone post-translational modifications (PTMs) have been reported to enhance microspore embryogenesis (ME), by altering histones acetylation or methylation. However, epigenetic mechanisms underlying ME induction efficiency are poorly understood. In this study, the epigenetic dynamics and the expression of genes associated with histone PTMs and ME induction were studied in two bread wheat cultivars with different ME response. Microspores isolated at 0, 3 and 5 days, treated with 0.7M mannitol (MAN) and 0.7M mannitol plus 0.4µM trichostatin A (TSA), which induced ME more efficiently, were analyzed. An additional control of gametophytic development was included. Microspores epigenetic state at the onset of ME induction was distinctive between cultivars by the ratio of H3 variants and their acetylated forms, the localization and percentage of labeled microspores with H3K9ac, H4K5ac, H4K16ac, H3K9me2 and H3K27me3, and the expression of genes related to pollen development. These results indicated that microspores of the high responding cultivar could be at a less advanced stage in pollen development. MAN and TSA resulted in a hyperacetylation of H3.2, with a greater effect of TSA. Histone PTMs were differentially affected by both treatments, with acetylation being most concerned. The effect of TSA was observed in the H4K5ac localization pattern at 3dT in the mid-low responding cultivar. Three gene networks linked to ME response were identified. TaHDT1, TaHAG2, TaYAO, TaNFD6-A, TabZIPF1 and TaAGO802-B, associated with pollen development, were down-regulated. TaHDA15, TaHAG3, TaHAM, TaYUC11D, Ta-2B-LBD16 TaMS1 and TaDRM3 constituted a network implicated in morphological changes by auxin signaling and cell wall modification up-regulated at 3dT. The last network included TaHDA18, TaHAC1, TaHAC4, TaABI5, TaATG18fD, TaSDG1a-7A and was related to ABA and ethylene hormone signaling pathways, DNA methylation and autophagy processes, reaching the highest expression at 5dT. The results indicated that TSA mainly modified the regulation of genes related to pollen and auxin signaling. This study represents a breakthrough in identifying the epigenetic dynamics and the molecular mechanisms governing ME induction efficiency, with relevance to recalcitrant wheat genotypes and other crops.

ZmLBD2 a maize (Zea mays L.) lateral organ boundaries domain (LBD) transcription factor enhances drought tolerance in transgenic Arabidopsis thaliana

Maize (Zea mays L.) is an annual gramineous herb and is among the world's most important crop species. Drought is the main factor contributing to maize yield reduction. The lateral organ boundaries domain (LBD) proteins belong to a class of higher-plant-specific transcription factors. LBD proteins usually include the highly conserved lateral organ boundaries (LOB) domains that play essential roles in plant growth and response to biotic stresses. However, few studies have addressed the biological functions of LBD genes associated with maize response to drought. Here we cloned the ZmLBD2 gene from maize and described its role in combating drought. Investigating ZmLBD2 subcellular localization, we show that it localizes to the cell nucleus and can specifically bind with inverted repeats of "GCGGCG". Under drought stress, Arabidopsis thaliana overexpressing ZmLBD2 performed better than the wild-type plants in terms of seed germination rates, root length, relative water content, fresh weight, chlorophyll content, proline content, and antioxidant enzyme content. Arabidopsis overexpressing ZmLBD2 contained less MDA, H₂O₂, and O 2 - than the wild-type plants. Our protein-protein interaction results indicate an interaction between the ZmLBD2 and ZmIAA5 genes. In conclusion, the ZmLBD2 gene positively regulates H₂O₂ homeostasis in plants, strengthening drought resistance.

Genome-wide association study for grain yield and component traits in bread wheat (Triticum aestivum L.)

Genomic regions governing days to heading (DH), grain filling duration (GFD), grain number per spike (GNPS), grain weight per spike (GWPS), plant height (PH), and grain yield (GY) were investigated in a set of 280 diverse bread wheat genotypes. The genome-wide association studies (GWAS) panel was genotyped using a 35K Axiom Array and phenotyped in five environments. The GWAS analysis showed a total of 27 Bonferroni-corrected marker-trait associations (MTAs) on 15 chromosomes representing all three wheat subgenomes. The GFD showed the highest MTAs (8), followed by GWPS (7), GY (4), GNPS (3), PH (3), and DH (2). Furthermore, 20 MTAs were identified with more than 10% phenotypic variation. A total of five stable MTAs (AX-95024590, AX-94425015, AX-95210025 AX-94539354, and AX-94978133) were identified in more than one environment and associated with the expression of DH, GFD, GNPS, and GY. Similarly, two novel pleiotropic genomic regions with associated MTAs i.e. AX-94978133 (4D) and AX-94539354 (6A) harboring co-localized QTLs governing two or more traits were also identified. In silico analysis revealed that the SNPs were located on important putative candidate genes such as F-box-like domain superfamily, Lateral organ boundaries, LOB, Thioredoxin-like superfamily Glutathione S-transferase, RNA-binding domain superfamily, UDP-glycosyltransferase family, Serine/threonine-protein kinase, Expansin, Patatin, Exocyst complex component Exo70, DUF1618 domain, Protein kinase domain involved in the regulation of grain size, grain number, growth and development, grain filling duration, and abiotic stress tolerance. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection (MAS).

Genome-Wide Identification of the Ginkgo (Ginkgo biloba L.) LBD Transcription Factor Gene and Characterization of Its Expression

Lateral organ boundaries domain (LBD) proteins are plant-specific transcription factors involved in various transcriptional regulation processes. We identified a total of 37 GbLBD genes in ginkgo, and based on gene structure and phylogenetic analysis, the GbLBD gene family was classified into class I (33, with the largest number of Id genes (16)) and class II (4). The ginkgo LBD gene was also analyzed regarding its chromosomal distributions, gene duplications, promoters, and introns/exons. In addition, gene expression profiling and real-time quantitative PCR analysis showed that the expression of 14 GbLBD genes differed in six different tissues and three developmental stages. The GbLBD gene of class II were highly expressed relative to the class I gene in all tissues and developmental stages, while class Id gene were generally at low levels or were not expressed, especially in seed developmental stages. The expression pattern analysis of cold/drought treatment and IAA/ABA hormone treatment showed that abiotic stress treatment could significantly induce the expression of GbLBD gene, of which class II genes played a key role in stress treatment. Our study provides a solid foundation for further evolutionary and functional analysis of the ginkgo LBD gene family.