Genome-wide analyses of the GbAP2 subfamily reveal the function of GbTOE1a in salt and drought stress tolerance in Ginkgo biloba.

The APETALA2 (AP2) transcription factors play crucial roles in plant growth and stage transition. Ginkgo biloba is an important medicinal plant renowned for the rich flavonoid content in its leaves. In this study, 18 GbAP2s were identified from the G. biloba genome and classified into three clusters. We found that the members of the euAP2 cluster, including four TOEs (GbTOE1a/1b/1c/3), exhibited a higher expression level in most samples compared to other members. Specifically, GbTOE1a may have a positive regulatory role in salt and drought stress responses. The overexpression of GbTOE1a in G. biloba calli resulted in a significant increase in the flavonoid content and upregulation of flavonoid biosynthesis genes, including PAL, 4CL, CHS, F3H, FLSs, F3'Hs, OMT, and DFRs. By contrast, the silencing of GbTOE1a in seedlings decreased the flavonoid content and the expression of flavonoid synthesizing genes. In addition, the silenced seedlings exhibited decreased antioxidant levels and a higher sensitivity to salt and drought treatments, suggesting a crucial role of GbTOE1a in G. biloba salt and drought tolerance. To the best of our knowledge, this was the first investigation into the identification and characterization of GbAP2s in G. biloba. Our results lay a foundation for further research on the regulatory role of the AP2 family in flavonoid synthesis and stress responses.

Comprehensive analysis of JAZ family members in Ginkgo biloba reveals the regulatory role of the GbCOI1/GbJAZs/GbMYC2 module in ginkgolide biosynthesis.

Ginkgo biloba L., an ancient relict plant known as a 'living fossil', has a high medicinal and nutritional value in its kernels and leaves. Ginkgolides are unique diterpene lactone compounds in G. biloba, with favorable therapeutic effects on cardiovascular and cerebrovascular diseases. Thus, it is essential to study the biosynthesis and regulatory mechanism of ginkgolide, which will contribute to quality improvement and medication requirements. In this study, the regulatory roles of the JAZ gene family and GbCOI1/GbJAZs/GbMYC2 module in ginkgolide biosynthesis were explored based on genome and methyl jasmonate-induced transcriptome. Firstly, 18 JAZ proteins were identified from G. biloba, and the gene characteristics and expansion patterns along with evolutionary relationships of these GbJAZs were analyzed systematically. Expression patterns analysis indicated that most GbJAZs expressed highly in the fibrous root and were induced significantly by methyl jasmonate. Mechanistically, yeast two-hybrid assays suggested that GbJAZ3/11 interacted with both GbMYC2 and GbCOI1, and several GbJAZ proteins could form homodimers or heterodimers between the GbJAZ family. Moreover, GbMYC2 is directly bound to the G-box element in the promoter of GbLPS, to regulate the biosynthesis of ginkgolide. Collectively, these results systematically characterized the JAZ gene family in G. biloba and demonstrated that the GbCOI1/GbJAZs/GbMYC2 module could regulate ginkgolides biosynthesis, which provides a novel insight for studying the mechanism of JA regulating ginkgolide biosynthesis.

The long noncoding RNAs lnc10 and lnc11 regulating flavonoid biosynthesis in Ginkgo biloba.

Although long non-coding RNAs have been recognized to play important roles in plant, their possible functions and potential mechanism in Ginkgo biloba flavonoid biosynthesis are poorly understood. Flavonoids are important secondary metabolites and healthy components of Ginkgo biloba. They have been widely used in food, medicine, and natural health products. Most previous studies have focused on the molecular mechanisms of structural genes and transcription factors that regulate flavonoid biosynthesis. Few reports have examined the biological functions of flavonoid biosynthesis by long non-coding RNAs in G. biloba. Long noncoding RNAs associated with flavonoid biosynthesis in G. biloba have been identified through RNA sequencing, but the function of lncRNAs has not been reported. In this study, the expression levels of lnc10 and lnc11 were identified. Quantitative real-time polymerase chain reaction analysis revealed that lnc10 and lnc11 were expressed in all detected organs, and they showed significantly higher levels in immature and mature leaves than in other organs. In addition, to fully identify the function of lnc10 and lnc11 in flavonoid biosynthesis in G. biloba, lnc10 and lnc11 were cloned from G. biloba, and were transformed into Arabidopsis and overexpressed. Compared with the wild type, the flavonoid content was increased in transgenic plants. Moreover, the RNA-sequencing analysis of wild-type, lnc10-overexpression, and lnc11-overexpression plants screened out 2019 and 2552 differentially expressed genes, and the transcript levels of structural genes and transcription factors associated with flavonoid biosynthesis were higher in transgenic Arabidopsis than in the wild type, indicating that lnc10 and lnc11 activated flavonoid biosynthesis in the transgenic lines. Overall, these results suggest that lnc10 and lnc11 positively regulate flavonoid biosynthesis in G. biloba.

Inferring the Regulatory Network of miRNAs on Terpene Trilactone Biosynthesis Affected by Environmental Conditions.

As a medicinal tree species, ginkgo (Ginkgo biloba L.) and terpene trilactones (TTLs) extracted from its leaves are the main pharmacologic activity constituents and important economic indicators of its value. The accumulation of TTLs is known to be affected by environmental stress, while the regulatory mechanism of environmental response mediated by microRNAs (miRNAs) at the post-transcriptional levels remains unclear. Here, we focused on grafted ginkgo grown in northwestern, southwestern, and eastern-central China and integrally analyzed RNA-seq and small RNA-seq high-throughput sequencing data as well as metabolomics data from leaf samples of ginkgo clones grown in natural environments. The content of bilobalide was highest among detected TTLs, and there was more than a twofold variation in the accumulation of bilobalide between growth conditions. Meanwhile, transcriptome analysis found significant differences in the expression of 19 TTL-related genes among ginkgo leaves from different environments. Small RNA sequencing and analysis showed that 62 of the 521 miRNAs identified were differentially expressed among different samples, especially the expression of miRN50, miR169h/i, and miR169e was susceptible to environmental changes. Further, we found that transcription factors (ERF, MYB, C3H, HD-ZIP, HSF, and NAC) and miRNAs (miR319e/f, miRN2, miRN54, miR157, miR185, and miRN188) could activate or inhibit the expression of TTL-related genes to participate in the regulation of terpene trilactones biosynthesis in ginkgo leaves by weighted gene co-regulatory network analysis. Our findings provide new insights into the understanding of the regulatory mechanism of TTL biosynthesis but also lay the foundation for ginkgo leaves' medicinal value improvement under global change.

Overview and Recent Progress on the Biosynthesis and Regulation of Flavonoids in Ginkgo biloba L.

Flavonoids and their derivatives play important roles in plants, such as exerting protective activity against biotic and abiotic stresses, functioning in visual signaling to attract pollinators, and regulating phytohormone activity. They are also important secondary metabolites that are beneficial to humans. Ginkgo biloba L. is a well-known relict plant considered to be a "living fossil". Flavonoids present in ginkgo leaves have antioxidant and anti-aging capacities and show good therapeutic effects on a variety of neurological diseases. To date, studies on flavonoids have mainly focused on their extraction, pharmacological effects, and component analysis and on the expression levels of the key genes involved. However, a systematic review summarizing the biosynthesis and regulatory mechanisms of ginkgo flavonoids is still lacking. Thus, this review was conducted to comprehensively introduce the biological characteristics, value, and utilization status of ginkgo; summarize the effects, biosynthetic pathways, and transcriptional regulation of flavonoids; and finally, discuss the factors (ecological factors, hormones, etc.) that regulate the biosynthesis of flavonoids in ginkgo. This review will provide a reference basis for future research on the biosynthesis and efficient utilization of flavonoids in ginkgo.

Genome-wide identification and expression analysis of the NRT genes in Ginkgo biloba under nitrate treatment reveal the potential roles during calluses browning.

Nitrate is a primary nitrogen source for plant growth, and previous studies have indicated a correlation between nitrogen and browning. Nitrate transporters (NRTs) are crucial in nitrate allocation. Here, we utilized a genome-wide approach to identify and analyze the expression pattern of 74 potential GbNRTs under nitrate treatments during calluses browning in Ginkgo, including 68 NITRATE TRANSPORTER 1 (NRT1)/PEPTIDE TRANSPORTER (PTR) (NPF), 4 NRT2 and 2 NRT3. Conserved domains, motifs, phylogeny, and cis-acting elements (CREs) were analyzed to demonstrate the evolutionary conservation and functional diversity of GbNRTs. Our analysis showed that the NPF family was divided into eight branches, with the GbNPF2 and GbNPF6 subfamilies split into three groups. Each GbNRT contained 108-214 CREs of 19-36 types, especially with binding sites of auxin and transcription factors v-myb avian myeloblastosis viral oncogene homolog (MYB) and basic helix-loop-helix (bHLH). The E1X1X2E2R motif had significant variations in GbNPFs, indicating changes in the potential dynamic proton transporting ability. The expression profiles of GbNRTs indicated that they may function in regulating nitrate uptake and modulating the signaling of auxin and polyphenols biosynthesis, thereby affecting browning in Ginkgo callus induction. These findings provide a better understanding of the role of NRTs during NO3- uptake and utilization in vitro culture, which is crucial to prevent browning and develop an efficient regeneration and suspension production system in Ginkgo.

Expression profiling of the phenylalanine ammonia-lyase (PAL) gene family in ginkgo biloba L.

The PAL gene family plays an important role in plant growth, development, and response to abiotic stresses and has been identified in a variety of plants. However, a systematic characterization is still lacking in Ginkgo biloba. Using a bioinformatics approach, 11 GbPAL members of the PAL gene family identified in ginkgo were identified in this study. The protein structure and physicochemical properties indicated that the GbPAL genes were highly similar. Based on their exon-intron structures, they can be classified into three groups. A total of 62 cis-elements for hormone, light, and abiotic stress responses were identified in the promoters of GbPAL genes, indicating that PAL is a multifunctional gene family. GbPAL genes were specifically expressed in different tissues and ploidy of ginkgo. These results provide a theoretical basis for further studies on the functional expression of the GbPAL genes.

Ginkgo biloba GbbZIP08 transcription factor is involved in the regulation of flavonoid biosynthesis.

Ginkgo biloba is the oldest relict plant on Earth and an economic plant resource derived from China. Flavonoids extracted from G. biloba are beneficial to the prevention and treatment of cardiovascular and cerebrovascular diseases. Basic leucine zipper (bZIP) transcription factors (TFs) have been recognized to play important roles in plant secondary metabolism. In this study, GbbZIP08 was isolated and characterized. It encodes a protein containing 154 amino acids, which belongs to hypocotyl 5 in group H of the bZIP family. Tobacco transient expression assay indicated that GbbZIP08 was localized in the plant nucleus. GbbZIP08 overexpression showed that the contents of total flavonoids, kaempferol, and anthocyanin in transgenic tobacco were significantly higher than those in the wild type. Transcriptome sequencing analysis revealed significant upregulation of structural genes in the flavonoid biosynthesis pathway. In addition, phytohormone signal transduction pathways, such as the abscisic acid, salicylic acid, auxin, and jasmonic acid pathways, were enriched with a large number of differentially expressed genes. TFs such as MYB, AP2, WRKY, NAC, bZIP, and bHLH, were also differentially expressed. The above results indicated that GbbZIP08 overexpression promoted flavonoid accumulation and increased the transcription levels of flavonoid-synthesis-related genes in plants.

Genome-Wide Identification of Expansin Gene Family and Their Response under Hormone Exposure in Ginkgo biloba L.

Expansins are pH-dependent enzymatic proteins that irreversibly and continuously facilitate cell-wall loosening and extension. The identification and comprehensive analysis of Ginkgo biloba expansins (GbEXPs) are still lacking. Here, we identified and investigated 46 GbEXPs in Ginkgo biloba. All GbEXPs were grouped into four subgroups based on phylogeny. GbEXPA31 was cloned and subjected to a subcellular localization assay to verify our identification. The conserved motifs, gene organization, cis-elements, and Gene Ontology (GO) annotation were predicted to better understand the functional characteristics of GbEXPs. The collinearity test indicated segmental duplication dominated the expansion of the GbEXPA subgroup, and seven paralogous pairs underwent strong positive selection during expansion. A majority of GbEXPAs were mainly expressed in developing Ginkgo kernels or fruits in transcriptome and real-time quantitative PCR (qRT-PCR). Furthermore, GbEXLA4, GbEXLA5, GbEXPA5, GbEXPA6, GbEXPA8, and GbEXPA24 were inhibited under the exposure of abiotic stresses (UV-B and drought) and plant hormones (ABA, SA, and BR). In general, this study expanded our understanding for expansins in Ginkgo tissues' growth and development and provided a new basis for studying GbEXPs in response to exogenous phytohormones.

Genome-Wide Identification, Evolutionary and Functional Analyses of WRKY Family Members in Ginkgo biloba.

WRKY transcription factors (TFs) are one of the largest families in plants which play essential roles in plant growth and stress response. Ginkgo biloba is a living fossil that has remained essentially unchanged for more than 200 million years, and now has become widespread worldwide due to the medicinal active ingredients in its leaves. Here, 37 WRKY genes were identified, which were distributed randomly in nine chromosomes of G. biloba. Results of the phylogenetic analysis indicated that the GbWRKY could be divided into three groups. Furthermore, the expression patterns of GbWRKY genes were analyzed. Gene expression profiling and qRT-PCR revealed that different members of GbWRKY have different spatiotemporal expression patterns in different abiotic stresses. Most of the GbWRKY genes can respond to UV-B radiation, drought, high temperature and salt treatment. Meanwhile, all GbWRKY members performed phylogenetic tree analyses with the WRKY proteins of other species which were known to be associated with abiotic stress. The result suggested that GbWRKY may play a crucial role in regulating multiple stress tolerances. Additionally, GbWRKY13 and GbWRKY37 were all located in the nucleus, while GbWRKY15 was located in the nucleus and cytomembrane.

Overexpression of the Ginkgo biloba dihydroflavonol 4-reductase gene GbDFR6 results in the self-incompatibility-like phenotypes in transgenic tobacco.

Although flavonoids play multiple roles in plant growth and development, the involvement in plant self-incompatibility (SI) have not been reported. In this research, the fertility of transgenic tobacco plants overexpressing the Ginkgo biloba dihydroflavonol 4-reductase gene, GbDFR6, were investigated. To explore the possible physiological defects leading to the failure of embryo development in transgenic tobacco plants, functions of pistils and pollen grains were examined. Transgenic pistils pollinated with pollen grains from another tobacco plants (either transgenic or wild-type), developed full of well-developed seeds. In contrast, in self-pollinated transgenic tobacco plants, pollen-tube growth was arrested in the upper part of the style, and small abnormal seeds developed without fertilization. Although the mechanism remains unclear, our research may provide a valuable method to create SI tobacco plants for breeding.

Transcriptome and proteome associated analysis of flavonoid metabolism in haploid Ginkgo biloba.

Having different number if genome copies affect transcription and metabolite production of plants. This may be due to different gene transcription and protein expression, but the reasons for this remains poorly known. Here we measured flavonoid content in leaves of three haploid and diploid grafted plants of Ginkgo biloba, a model gymnosperm important economically for its flavonoid content. We reported the first combined transcriptomic and proteomic analysis of the difference in flavonoid content in three haploid ginkgos to investigate the effect of haploidy. Haploids had always smaller leaves and flavonoid content than the diploids. The selected haploid had also generally lower gene dosage than the selected diploid, with 1149 up-regulated (46.8 %) and 1309 down-regulated (53.2 %) among 2452 differentially expressed genes (DEGs). Of 686 differentially expressed proteins (DEPs) detected, 289 proteins (42.1 %) were upregulated, and 397 proteins (57.9 %) were downregulated in haploids. A particular attention deserves the downregulation of PAL, PAM, FLS, OMT1 hub genes involved in flavonoid biosynthesis regulation. Our study confirms the trend of haploids to have lower metabolic contents and points that lower flavonoid content in ginkgo monoploids could be due to reduced dosage of the corresponding regulatory genes and downregulation of genes involved in flavonoid synthesis.

Genome-wide identification and expression analysis of NAC family genes in Ginkgo biloba L.

NAC (NAM, ATAF, CUC2) transcription factors constitute one of the largest families of plant-specific transcription factors with important roles in plant growth and development and in biotic and abiotic stresses. The physicochemical properties, gene structure, cis-acting elements and expression patterns of NAC transcription factors in Ginkgo biloba were analysed using bioinformatics, and expression of this gene family was analysed via quantitative reverse transcription PCR. The family of G. biloba NAC transcription factors had 50 members, distributed on 12 chromosomes and divided into 11 groups. Members in the same group share a similar gene structure and motif distribution. Transcriptome data analysis of G. biloba showed that 35 genes were expressed in eight tissues. Correlation analysis suggested that GbNAC007 and GNAC008 might be involved in flavonoid biosynthesis. Expression levels of 12 GbNACs under cold, het, and salt stresses were analysed. Results indicate that NAC transcription factors play an important role in response to abiotic stresses. This study provides a reference for the functional analysis of the G. biloba family of NAC transcription factors, as well as a resource for studies on the involvement of this family in responses to abiotic stresses and flavonoid biosynthesis.

A gene cluster in Ginkgo biloba encodes unique multifunctional cytochrome P450s that initiate ginkgolide biosynthesis.

The ginkgo tree (Ginkgo biloba) is considered a living fossil due to its 200 million year's history under morphological stasis. Its resilience is partly attributed to its unique set of specialized metabolites, in particular, ginkgolides and bilobalide, which are chemically complex terpene trilactones. Here, we use a gene cluster-guided mining approach in combination with co-expression analysis to reveal the primary steps in ginkgolide biosynthesis. We show that five multifunctional cytochrome P450s with atypical catalytic activities generate the tert-butyl group and one of the lactone rings, characteristic of all G. biloba trilactone terpenoids. The reactions include scarless C-C bond cleavage as well as carbon skeleton rearrangement (NIH shift) occurring on a previously unsuspected intermediate. The cytochrome P450s belong to CYP families that diversifies in pre-seed plants and gymnosperms, but are not preserved in angiosperms. Our work uncovers the early ginkgolide pathway and offers a glance into the biosynthesis of terpenoids of the Mesozoic Era.

Dynamic transcriptome and network-based analysis of yellow leaf mutant Ginkgo biloba.

BACKGROUND: Golden leaf in autumn is a prominent feature of deciduous tree species like Ginkgo biloba L., a landscape tree widely cultivated worldwide. However, little was known about the molecular mechanisms of leaf yellowing, especially its dynamic regulatory network. Here, we performed a suite of comparative physiological and dynamic transcriptional analyses on the golden-leaf cultivar and the wild type (WT) ginkgo to investigate the underlying mechanisms of leaf yellowing across different seasons. RESULTS: In the present study, we used the natural bud mutant cultivar with yellow leaves "Wannianjin" (YL) as materials. Physiological analysis revealed that higher ratios of chlorophyll a to chlorophyll b and carotenoid to chlorophyll b caused the leaf yellowing of YL. On the other hand, dynamic transcriptome analyses showed that genes related to chlorophyll metabolism played key a role in leaf coloration. Genes encoding non-yellow coloring 1 (NYC1), NYC1-like (NOL), and chlorophyllase (CLH) involved in the degradation of chlorophyll were up-regulated in spring. At the summer stage, down-regulated HEMA encoding glutamyl-tRNA reductase functioned in chlorophyll biosynthesis, while CLH involved in chlorophyll degradation was up-regulated, causing a lower chlorophyll accumulation. In carotenoid metabolism, genes encoding zeaxanthin epoxidase (ZEP) and 9-cis-epoxy carotenoid dioxygenase (NCED) showed significantly different expression levels in the WT and YL. Moreover, the weighted gene co-expression network analysis (WGCNA) suggested that the most associated transcriptional factor, which belongs to the AP2/ERF-ERF family, was engaged in regulating pigment metabolism. Furthermore, quantitative experiments validated the above results. CONCLUSIONS: By comparing the golden-leaf cultivar and the wide type of ginkgo across three seasons, this study not only confirm the vital role of chlorophyll in leaf coloration of YL but also provided new insights into the seasonal transcriptome landscape and co-expression network. Our novel results pinpoint candidate genes for further wet-bench experiments in tree species.

High-Depth Transcriptome Reveals Differences in Natural Haploid Ginkgo biloba L. Due to the Effect of Reduced Gene Dosage.

As a representative of gymnosperms, the discovery of natural haploids of Ginkgo biloba L. has opened a new door for its research. Haploid germplasm has always been a research material of interest to researchers because of its special characteristics. However, we do not yet know the special features and mechanisms of haploid ginkgo following this significant discovery. In this study, we conducted a homogenous garden experiment on haploid and diploid ginkgo to explore the differences in growth, physiology and biochemistry between the two. Additionally, a high-depth transcriptome database of both was established to reveal their transcriptional differences. The results showed that haploid ginkgo exhibited weaker growth potential, lower photosynthesis and flavonoid accumulation capacity. Although the up-regulated expression of DEGs in haploid ginkgo reached 46.7% of the total DEGs in the whole transcriptome data, the gene sets of photosynthesis metabolic, glycolysis/gluconeogenesis and flavonoid biosynthesis pathways, which were significantly related to these differences, were found to show a significant down-regulated expression trend by gene set enrichment analysis (GSEA). We further found that the major metabolic pathways in the haploid ginkgo transcriptional database were down-regulated in expression compared to the diploid. This study reveals for the first time the phenotypic, growth and physiological differences in haploid ginkgos, and demonstrates their transcriptional patterns based on high-depth transcriptomic data, laying the foundation for subsequent in-depth studies of haploid ginkgos.

Systematic Identification and Expression Analysis of the Auxin Response Factor (ARF) Gene Family in Ginkgo biloba L.

Auxin participates in various physiological and molecular response-related developmental processes and is a pivotal hormone that regulates phenotypic formation in plants. Auxin response factors (ARFs) are vital transcription factors that mediate downstream auxin signaling by explicitly binding to auxin-responsive genes' promoters. Here, to investigate the possible developmental regulatory functions of ARFs in Ginkgo biloba, through employing comprehensive bioinformatics, we recognized 15 putative GbARF members. Conserved domains and motifs, gene and protein structure, gene duplication, GO enrichment, transcriptome expression profiles, and qRT-PCR all showed that Group I and III members were highly conserved. Among them, GbARF10b and GbARF10a were revealed as transcriptional activators in the auxin response for the development of Ginkgo male flowers through sequences alignment, cis-elements analysis and GO annotation; the results were corroborated for the treatment of exogenous SA. Moreover, the GbARFs expansion occurred predominantly by segmental duplication, and most GbARFs have undergone purifying selection. The Ka/Ks ratio test identified the functional consistence of GbARF2a and GbARF2c, GbARF10b, and GbARF10a in tissue expression profiles and male flower development. In summary, our study established a new research basis for exploring Ginkgo GbARF members' roles in floral organ development and hormone response.

GinkgoDB: an ecological genome database for the living fossil, Ginkgo biloba.

Although an increasing number of species-specific databases have been constructed to facilitate the global research community, comprehensive databases remain rare with multidimensional research resources rather than genomic data solely, in particular for non-model species. Here we introduced GinkgoDB, the ecological genome database for the world-renowned living fossil, Ginkgo biloba L., a tree species with extreme longevity, exceptional resistance, global landscape application and dioecy. The present version of GinkgoDB consists of four modules of genome, occurrence, quadrat and gallery. The genome module includes two versions of chromosome-level assembled genomes with high-quality annotation, expression profiles of each sex for different tissues and variants from 545 ginkgo genomes. A suite of tools were also deployed online for search, blast and further analyses. The occurrence module contains distribution information of over one million records across the world with an emphasis on 8926 individual-level tagged trees with phenotypic trait data in eastern Asia. The quadrat module provides monitoring data of both species and habitats from 27 plots of natural ginkgo forests and periodic data of functional traits measured for the entire plant communities. The gallery module collects field morphological photos and scanned phenotypic images of various ginkgo materials. GinkgoDB showcases a comprehensive and persistently growing database for a single non-model species. Database URL: https://ginkgo.zju.edu.cn/.

Genome-wide transcriptome analysis reveals the regulatory network governing terpene trilactones biosynthesis in Ginkgo biloba.

Ginkgo biloba L. is currently the only remaining gymnosperm of the Ginkgoaceae Ginkgo genus, and its history can be traced back to the Carboniferous 200 million years ago. Terpene trilactones (TTLs) are one of the main active ingredients in G. biloba, including ginkgolides and bilobalide. They have a good curative effect on cardiovascular and cerebrovascular diseases because of their special antagonistic effect on platelet-activating factors. Therefore, it is necessary to deeply mine genes related to TTLs and to analyze their transcriptional regulation mechanism, which will hold vitally important scientific and practical significance for quality improvement and regulation of G. biloba. In this study, we performed RNA-Seq on the root, stem, immature leaf, mature leaf, microstrobilus, ovulate strobilus, immature fruit and mature fruit of G. biloba. The TTL regulatory network of G. biloba in different organs was revealed by different transcriptomic analysis strategies. Weighted gene co-expression network analysis (WGCNA) revealed that the five modules were closely correlated with organs. The 12 transcription factors, 5 structural genes and 24 Cytochrome P450 (CYP450) were identified as candidate regulators for TTL accumulation by WGCNA and cytoscape visualization. Finally, 6 APETALA2/ethylene response factors, 2 CYP450s and bHLH were inferred to regulate the metabolism of TTLs by correlation analysis. This study is the comprehensive in authenticating transcription factors, structural genes and CYP450 involved in TTL biosynthesis, thereby providing molecular evidence for revealing the comprehensive regulatory network involved in TTL metabolism in G. biloba.

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.