The nearly complete genome of Ginkgo biloba illuminates gymnosperm evolution

The evolution and development of sporangia-The fundamental reproductive organ of land plant sporophytes

A key innovation of land plants is the origin and evolution of the sporangium, the fundamental reproductive structure of the diploid sporophyte. In vascular plants, whether the structure is a cone, fertile leaf, or flower-all are clusters of sporangia. The evolution of morphologically distinct sporangia (heterospory) and retention of the gametophyte evolved three times independently as a prerequisite for the evolution of seeds. This review summarizes the development of vascular plant sporangia, molecular genetics of angiosperm sporangia, and provides a framework to investigate evolution and development in vascular plant sporangia.

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

Ginkgo Biloba Bioactive Phytochemicals against Age-Related Diseases: Evidence from a Stepwise, High-Throughput Research Platform

The seeds of ginkgo biloba L (GB) have been widely used worldwide. This study investigated the bioefficacies of whole GB seed powder (WGP) retaining the full nutrients of ginkgo against aging, atherosclerosis, and fatigue. The experimental results indicated that WGP lowered brain monoamine oxidase and serum malondialdehyde levels, enhanced thymus/spleen indexes, and improved learning ability, and delayed aging in senescent mice. WGP regulated lipid levels and prevented atherosclerosis by reducing triglycerides, lowering low-density lipoprotein cholesterol, increasing high-density lipoprotein cholesterol, and decreasing the atherosclerosis index. WGP improved exercise performance by reducing blood lactate accumulation and extending exhaustive swimming and climbing times, improved energy storage by increasing muscle/liver glycogen levels, and relieved physical fatigue. Network pharmacology analysis revealed 270 potential targets of WGP that play roles in cellular pathways related to inflammation inhibition, metabolism regulation, and anti-cellular senescence, etc. Protein-protein interaction analysis identified 10 hub genes, including FOS, ESR1, MAPK8, and SP1 targets. Molecular docking and molecular dynamics simulations showed that the bioactive compounds of WGP bound well to the targets. This study suggests that WGP exerts prominent health-promoting effects through multiple components, targets, and pathways.

Tree Longevity: Multifaceted Genetic Strategies and Beyond

Old trees are remarkable for their ability to endure for centuries or even millennia, acting as recordkeepers of historical climate and custodians of genetic diversity. The secret to their longevity has long been a subject of fascination. Despite the challenges associated with studying old trees, such as massive size, slow growth rate, long lifespan and often remote habitat, accumulating studies have investigated the mechanisms underlying tree aging and longevity over the past decade. The recent publication of high-quality genomes of long-lived tree species, coupled with research on stem cell function and secondary metabolites in longevity, has brought us closer to unlocking the secrets of arboreal longevity. This review provides an overview of the global distribution of old trees and examines the environmental and anthropogenic factors that shape their presence. We summarize the contributions of physiological characteristics, stem cell activity, and immune system responses to their extraordinary longevity. We also explore the genetic and epigenetic 'longevity code', which consists of resistance and defense genes, DNA repair genes and patterns of DNA methylation modification. Further, we highlight key areas for future research that could enhance our understanding of the mechanisms underlying tree longevity.

Hyperspectral Imaging Combined with Deep Transfer Learning to Evaluate Flavonoids Content in Ginkgo biloba Leaves

Ginkgo biloba is a famous economic tree. Ginkgo leaves have been utilized as raw materials for medicines and health products due to their rich active ingredient composition, especially flavonoids. Since the routine measurement of total flavones is time-consuming and destructive, rapid, non-destructive detection of total flavones in ginkgo leaves is of significant importance to producers and consumers. Hyperspectral imaging technology is a rapid and non-destructive technique for determining the total flavonoid content. In this study, we discuss five modeling methods, and three spectral preprocessing methods are discussed. Bayesian Ridge (BR) and multiplicative scatter correction (MCS) were selected as the best model and the best pretreatment method, respectively. The spectral prediction results based on the BR + MCS treatment were very accurate (RTest2 = 0.87; RMSETest = 1.03 mg/g), showing a high correlation with the analytical measurements. In addition, we also found that the more and deeper the leaf cracks, the higher the flavonoid content, which helps to evaluate leaf quality more quickly and easily. In short, hyperspectral imaging is an effective technique for rapid and accurate determination of total flavonoids in ginkgo leaves and has great potential for developing an online quality detection system for ginkgo leaves.

De novo Whole-Genome Assembly of the 10-Gigabase Fokienia Hodginsii Genome to Reveal Differential Epigenetic Events Between Callus and Xylem

Fokienia hodginsii (F. hodginsii), belonging to the genus Fokienia of the Cupressaceae. F. hodginsii has significant application value due to its wood properties and great research value in evolutionary studies as a gymnosperm. However, the genome of F. hodginsii remains unknown due to the large size of gymnosperms genome. Pacific Bioscience sequencing, Hi-C mapping, whole-genome Bisulfite Sequencing (BS-Seq), long-read isoform sequencing (Iso-Seq), direct RNA sequencing (DRS), quantitative proteomics, and metabonomics analysis are employed to facilitate genome assembly, gene annotation, and investigation into epigenetic mechanisms. In this study, the 10G F. hodginsii genome is assembled into 11 chromosomes. Furthermore, 50 521 protein-coding genes are annotated and determined that 65% of F. hodginsii genome comprises repetitive sequences. It is discovered that transposable element (TE)-including introns is associated with higher expression. The DNA methylome of F. hodginsii reveals that xylem has a higher DNA methylation level compared to callus. Moreover, DRS reveals the significant alterations in RNA full-length ratio, which potentially associated with poly(A) length (PAL) and alternative polyadenylation (APA). Finally, the morphology measurement and metabonomics analysis revealed the difference of 14 cultivars. In summary, the genomes and epigenetics datasets provide a molecular basis for callus formation in the gymnosperm family.

Three-dimensional reconstruction and multiomics analysis reveal a unique pattern of embryogenesis in Ginkgo biloba

Ginkgo (Ginkgo biloba L.) is one of the earliest extant species in seed plant phylogeny. Embryo development patterns can provide fundamental evidence for the origin, evolution, and adaptation of seeds. However, the architectural and morphological dynamics during embryogenesis in G. biloba remain elusive. Herein, we obtained over 2,200 visual slices from 3 stages of embryo development using micro-computed tomography imaging with improved staining methods. Based on 3-dimensional (3D) spatiotemporal pattern analysis, we found that a shoot apical meristem with 7 highly differentiated leaf primordia, including apical and axillary leaf buds, is present in mature Ginkgo embryos. 3D rendering from the front, top, and side views showed 2 separate transport systems of tracheids located in the hypocotyl and cotyledon, representing a unique pattern of embryogenesis. Furthermore, the morphological dynamic analysis of secretory cavities indicated their strong association with cotyledons during development. In addition, we identified genes GbLBD25a (lateral organ boundaries domain 25a), GbCESA2a (cellulose synthase 2a), GbMYB74c (myeloblastosis 74c), GbPIN2 (PIN-FORMED 2) associated with vascular development regulation, and GbWRKY1 (WRKYGOK 1), GbbHLH12a (basic helix-loop-helix 12a), and GbJAZ4 (jasmonate zim-domain 4) potentially involved in the formation of secretory cavities. Moreover, we found that flavonoid accumulation in mature embryos could enhance postgerminative growth and seedling establishment in harsh environments. Our 3D spatial reconstruction technique combined with multiomics analysis opens avenues for investigating developmental architecture and molecular mechanisms during embryogenesis and lays the foundation for evolutionary studies of embryo development and maturation.

Genome-wide characterization of RR gene family members in Zanthoxylum armatum and the subsequent functional characterization of the C-type RR

Response Regulators (RRs) are crucial regulators in plant development and stress responses, comprising A-type, B-type, C-type, and pseudo-RR subfamilies. However, previous studies have often focused on specific subfamilies, which restricts our understanding of the complete RR gene family. In this study, we conducted a comprehensive analysis of 63 RR members from Zanthoxylum armatum, using phylogenetic relationships, motif composition, cis-acting elements, gene duplication and collinearity analyses. Segmental repeats among ZaRR genes enhanced the various environmental adaptabilities of Z. armatum, and the B-type ZaRR exhibited significant collinearity with the RRs in P. trichocarpa and C. sinensis. Cis-element analysis indicated ZaRRs play a significant role in abiotic stress and phytohormone pathways, particularly in light, drought, cold, abscisic acid (ABA) and salicylic acid (SA) responses. Abundant Ethylene Response Factor (ERF) and reproduction-associated binding sites in ZaRR promoters suggested their roles in stress and reproductive processes. A-type ZaRRs were implicated in plant vegetative and reproductive growth, whereas B-type ZaRRs contributed to both growth and stress responses. C-type ZaRRs were associated with plant reproductive growth, whereas pseudo-RRs may function in plant stress responses, such as water logging, cold, and response to ethylene (ETH), SA, and jasmonic acid (JA). Ectopic expression of ZaRR24, a C-type RR, inhibits growth, induces early flowering, and shortens fruit length in Arabidopsis. ZaRR24 overexpression also affected the expression of A- and B-type RRs, as well as floral meristem and organ identity genes. These findings establish a solid and comprehensive foundation for RR gene research in Z. armatum, and provide a platform for investigating signal transduction in other woody plants.

Functional characterization of AcWRKY94 in response to Pseudomonas syringae pv. actinidiae in kiwifruit

WRKY transcription factors are essential for coping with various biotic stresses. Pseudomonas syringae pv. actinidiae (Psa)-induced kiwifruit canker is a major problem restricting kiwifruit yield. Nevertheless, it's unclear how the kiwifruit WRKY genes respond to Psa. Through genome-wide identification, 112 WRKY members were found in 'Hongyang' genome in this work. Promoter analysis revealed that there were many cis-acting elements associated with stress responses in the AcWRKY gene's promoter region. According to transcriptomic analysis, 90 of the AcWRKY genes were differently expressed following Psa, salicylic acid (SA), or methyl jasmonate (MeJA) treatments. Almost all group III WRKYs were responsive to at least one of these treatments, with tissue-specific expression patterns. Quantitative RT-PCR study provided more evidence that Psa and SA treatments significantly induced the expression of the group III WRKY gene AcWRKY94, whereas MeJA treatment repressed it. AcWRKY94 was a transcriptionally active protein localized in the nucleus. Transient overexpression of AcWRKY94 in the leaves of 'Hongyang' enhanced the resistance of kiwifruit to Psa. Overexpression of AcWRKY94 in kiwifruit callus remarkably promoted the expression of PR and JAZ genes associated with SA and JA signals, respectively. These data imply that AcWRKY94 controls the signaling pathway dependent on SA and JA, thereby enhancing resistance to Psa. Taken together, this study establishes the basis for functional research on WRKY genes and provides important information for elucidating the resistance mechanism of kiwifruit canker disease.

Beyond co-expression: pathway discovery for plant pharmaceuticals

Plant natural products have been an important source of medicinal molecules since ancient times. To gain access to the whole diversity of these molecules for pharmaceutical applications, it is important to understand their biosynthetic origins. Whilst co-expression is a reliable tool for identifying gene candidates, a variety of complementary methods can aid in screening or refining candidate selection. Here, we review recently employed plant biosynthetic pathway discovery approaches, and highlight future directions in the field.

Formation mechanism and regulation analysis of trumpet leaf in Ginkgo biloba L

Introduction

The research on plant leaf morphology is of great significance for understanding the development and evolution of plant organ morphology. As a relict plant, the G. biloba leaf morphology typically exhibits bifoliate and peltate forms. However, throughout its long evolutionary history, Ginkgo leaves have undergone diverse changes.

Methods

This study focuses on the distinct "trumpet" leaves and normal fan-shaped leaves of G. biloba for analysis of their phenotypes, photosynthetic activity, anatomical observations, as well as transcriptomic and metabolomic analyses.

Results

The results showed that trumpet-shaped G. biloba leaves have fewer cells, significant morphological differences between dorsal and abaxial epidermal cells, leading to a significantly lower net photosynthetic rate. Additionally, this study found that endogenous plant hormones such as GA, auxin, and JA as well as metabolites such as flavonoids and phenolic acids play roles in the formation of trumpet-shaped G. biloba leaves. Moreover, the experiments revealed the regulatory mechanisms of various key biological processes and gene expressions in the trumpet-shaped leaves of G. biloba.

Discussion

Differences in the dorsal and abdominal cells of G. biloba leaves can cause the leaf to curl, thus reducing the overall photosynthetic efficiency of the leaves. However, the morphology of plant leaves is determined during the primordia leaf stage. In the early stages of leaf development, the shoot apical meristem (SAM) determines the developmental morphology of dicotyledonous plant leaves. This process involves the activity of multiple gene families and small RNAs. The establishment of leaf morphology is complexly regulated by various endogenous hormones, including the effect of auxin on cell walls. Additionally, changes in intracellular ion concentrations, such as fluctuations in Ca2+ concentration, also affect cell wall rigidity, thereby influencing leaf growth morphology.

Metabolomic and Transcriptomic Profiles in Diverse Brassica oleracea Crops Provide Insights into the Genetic Regulation of Glucosinolate Profiles

Glucosinolates (GSLs) are plant secondary metabolites commonly found in the cruciferous vegetables of the Brassicaceae family, offering health benefits to humans and defense against pathogens and pests to plants. In this study, we investigated 23 GSL compounds' relative abundance in four tissues of five different Brassica oleracea morphotypes. Using the five corresponding high-quality B. oleracea genome assemblies, we identified 183 GSL-related genes and analyzed their expression with mRNA-Seq data. GSL abundance and composition varied strongly, among both tissues and morphotypes, accompanied by different gene expression patterns. Interestingly, broccoli exhibited a nonfunctional AOP2 gene due to a conserved 2OG-FeII_Oxy domain loss, explaining the unique accumulation of two health-promoting GSLs. Additionally, transposable element (TE) insertions were found to affect the gene structure of MAM3 genes. Our findings deepen the understanding of GSL variation and genetic regulation in B. oleracea morphotypes, providing valuable insights for breeding with tailored GSL profiles in these crops.

A reference-grade genome of the xerophyte Ammopiptanthus mongolicus sheds light on its evolution history in legumes and drought-tolerance mechanisms

Plants that grow in extreme environments represent unique sources of stress-resistance genes and mechanisms. Ammopiptanthus mongolicus (Leguminosae) is a xerophytic evergreen broadleaf shrub native to semi-arid and desert regions; however, its drought-tolerance mechanisms remain poorly understood. Here, we report the assembly of a reference-grade genome for A. mongolicus, describe its evolutionary history within the legume family, and examine its drought-tolerance mechanisms. The assembled genome is 843.07 Mb in length, with 98.7% of the sequences successfully anchored to the nine chromosomes of A. mongolicus. The genome is predicted to contain 47 611 protein-coding genes, and 70.71% of the genome is composed of repetitive sequences; these are dominated by transposable elements, particularly long-terminal-repeat retrotransposons. Evolutionary analyses revealed two whole-genome duplication (WGD) events at 130 and 58 million years ago (mya) that are shared by the genus Ammopiptanthus and other legumes, but no species-specific WGDs were found within this genus. Ancestral genome reconstruction revealed that the A. mongolicus genome has undergone fewer rearrangements than other genomes in the legume family, confirming its status as a "relict plant". Transcriptomic analyses demonstrated that genes involved in cuticular wax biosynthesis and transport are highly expressed, both under normal conditions and in response to polyethylene glycol-induced dehydration. Significant induction of genes related to ethylene biosynthesis and signaling was also observed in leaves under dehydration stress, suggesting that enhanced ethylene response and formation of thick waxy cuticles are two major mechanisms of drought tolerance in A. mongolicus. Ectopic expression of AmERF2, an ethylene response factor unique to A. mongolicus, can markedly increase the drought tolerance of transgenic Arabidopsis thaliana plants, demonstrating the potential for application of A. mongolicus genes in crop improvement.

Evolutionary divergence of CXE gene family in green plants unveils that PtoCXEs overexpression reduces fungal colonization in transgenic Populus

Plant enzymes significantly contribute to the rapidly diversified metabolic repertoire since the colonization of land by plants. Carboxylesterase is just one of the ubiquitous, multifunctional and ancient enzymes that has particularly diversified during plant evolution. This study provided a status on the carboxylesterase landscape within Viridiplantae. A total of 784 carboxylesterases were identified from the genome of 31 plant species representing nine major lineages of sequenced Viridiplantae and divided into five clades based on phylogenetic analysis. Clade I carboxylesterase genes may be of bacterial origin and then expanded and diversified during plant evolution. Clade II was first gained in the ancestor of bryophytes after colonization of land by plants, Clade III and Clade IV in ferns which were considered the most advanced seedless vascular plants, while Clade V was gained in seed plants. To date, the functions of carboxylesterase genes in woody plants remain unclear. In this study, 51 carboxylesterase genes were identified from the genome of Populus trichocarpa and further divided into eight classes. Tandem and segmental duplication events both contributed to the expansion of carboxylesterase genes in Populus. Although carboxylesterase genes were proven to enhance resistance to pathogens in many herbaceous species, relevant researches on forest trees are still needed. In this study, pathogen incubation assays showed that overexpressing of six Class VI carboxylesterases in Populus tomentosa, to a greater or lesser degree, reduced colonization of detached leaves by fungus Cytospora chrysosperma. A significant difference was also found in functional divergence patterns for genes derived from different gene duplication events. Functional differentiation of duplicated carboxylesterase genes in Populus was proved for the first time by in vivo physiological analysis. The identification of the potentially anti-fungal PtoCXE06 gene also laid a theoretical foundation for promoting the genetic improvement of disease-resistance traits in forest trees.

Research progress on the types, functions, biosynthesis, and metabolic regulation of ginkgo terpenoids

Ginkgo biloba L. is a relict plant endemic to China that is commonly considered a "living fossil". It contains unique medicinal compounds that play important roles in its response to various stresses and help maintain human health. Ginkgo terpenoids are known to be important active ingredients but have received less attention than flavonoids. Hence, this review focuses on recent progress in research on the pharmacological effects of ginkgo terpenoid and the bioactivities of different terpenoid monomers. Many key structural genes, enzyme-encoding genes, transcription factors, and noncoding RNAs involved in the ginkgo terpenoid pathway were identified. Finally, many external factors (ecological factors, hormones, etc.) that regulate the biosynthesis and metabolism of terpenoids were proposed. All these findings improve the understanding of the biosynthesis, accumulation, and medicinal functions of terpenoids. Finally, this review includes an in-depth discussion regarding the limitations of terpenoid-related studies and potential future research directions.

Deciphering recursive polyploidization in Lamiales and reconstructing their chromosome evolutionary trajectories

Lamiales is an order of core eudicots with abundant diversity, and many Lamiales plants have important medicinal and ornamental values. Here, we comparatively reanalyzed 11 Lamiales species with well-assembled genome sequences and found evidence that Lamiales plants, in addition to a hexaploidization or whole-genome triplication (WGT) shared by core eudicots, experienced further polyploidization events, establishing new groups in the order. Notably, we identified a whole-genome duplication (WGD) occurred just before the split of Scrophulariaceae from the other Lamiales families, such as Acanthaceae, Bignoniaceae, and Lamiaceae, suggesting its likely being the causal reason for the establishment and fast divergence of these families. We also found that a WGT occurred ∼68 to 78 million years ago (Mya), near the split of Oleaceae from the other Lamiales families, implying that it may have caused their fast divergence and the establishment of the Oleaceae family. Then, by exploring and distinguishing intra- and intergenomic chromosomal homology due to recursive polyploidization and speciation, respectively, we inferred that the Lamiales ancestral cell karyotype had 11 proto-chromosomes. We reconstructed the evolutionary trajectories from these proto-chromosomes to form the extant chromosomes in each Lamiales plant under study. We must note that most of the inferred 11 proto-chromosomes, duplicated during a WGD thereafter, have been well preserved in jacaranda (Jacaranda mimosifolia) genome, showing the credibility of the present inference implementing a telomere-centric chromosome repatterning model. These efforts are important to understand genome repatterning after recursive polyploidization, especially shedding light on the origin of new plant groups and angiosperm cell karyotype evolution.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Piece and parcel of gymnosperm organellar genomes

MAIN CONCLUSION

Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity.

Efficient purging of deleterious mutations contributes to the survival of a rare conifer

Cupressaceae is a conifer family rich in plants of horticultural importance, including Cupressus, Chamaecyparis, Juniperus, and Thuja, yet genomic surveys are lacking for this family. Cupressus gigantea, one of the many rare conifers that are threatened by climate change and anthropogenic habitat fragmentation, plays an ever-increasing role in ecotourism in Tibet. To infer how past climate change has shaped the population evolution of this species, we generated a de novo chromosome-scale genome (10.92 Gb) and compared the species' population history and genetic load with that of a widespread close relative, C. duclouxiana. Our demographic analyses, based on 83 re-sequenced individuals from multiple populations of the two species, revealed a sharp decline of population sizes during the first part of the Quaternary. However, populations of C. duclouxiana then started to recover, while C. gigantea populations continued to decrease until recently. The total genomic diversity of C. gigantea is smaller than that of C. duclouxiana, but contrary to expectations, C. gigantea has fewer highly and mildly deleterious mutations than C. duclouxiana, and simulations and statistical tests support purifying selection during prolonged inbreeding as the explanation. Our results highlight the evolutionary consequences of decreased population size on the genetic burden of a long-lived endangered conifer with large genome size and suggest that genetic purging deserves more attention in conservation management.

The Ginkgo biloba microRNA160-ERF4 module participates in terpene trilactone biosynthesis

Terpene trilactones (TTLs) are important secondary metabolites in ginkgo (Ginkgo biloba); however, their biosynthesis gene regulatory network remains unclear. Here, we isolated a G. biloba ethylene response factor 4 (GbERF4) involved in TTL synthesis. Overexpression of GbERF4 in tobacco (Nicotiana tabacum) significantly increased terpenoid content and upregulated the expression of key enzyme genes (3-hydroxy-3-methylglutaryl-CoA reductase [HMGR], 3-hydroxy-3-methylglutaryl-CoA synthase [HMGS], 1-deoxy-D-xylulose-5-phosphate reductoisomerase [DXR], 1-deoxy-D-xylulose-5-phosphate synthase [DXS], acetyl-CoA C-acetyltransferase [AACT], and geranylgeranyl diphosphate synthase [GGPPS]) in the terpenoid pathway in tobacco, suggesting that GbERF4 functions in regulating the synthesis of terpenoids. The expression pattern analysis and previous microRNA (miRNA) sequencing showed that gb-miR160 negatively regulates the biosynthesis of TTLs. Transgenic experiments showed that overexpression of gb-miR160 could significantly inhibit the accumulation of terpenoids in tobacco. Targeted inhibition and dual-luciferase reporter assays confirmed that gb-miR160 targets and negatively regulates GbERF4. Transient overexpression of GbERF4 increased TTL content in G. biloba, and further transcriptome analysis revealed that DXS, HMGS, CYPs, and transcription factor genes were upregulated. In addition, yeast 1-hybrid and dual-luciferase reporter assays showed that GbERF4 could bind to the promoters of the HMGS1, AACT1, DXS1, levopimaradiene synthase (LPS2), and GGPPS2 genes in the TTL biosynthesis pathway and activate their expression. In summary, this study investigated the molecular mechanism of the gb-miR160-GbERF4 regulatory module in regulating the biosynthesis of TTLs. It provides information for enriching the understanding of the regulatory network of TTL biosynthesis and offers important gene resources for the genetic improvement of G. biloba with high contents of TTLs.

Brassinosteroid-dependent phosphorylation of PHOSPHATE STARVATION RESPONSE2 reduces its DNA-binding ability in rice

Brassinosteroids (BRs) are widely used as plant growth regulators in modern agriculture. Understanding how BRs regulate nutrient signaling is crucial for reducing fertilizer usage. Here we elucidate that the central BR signaling inhibitor GSK3/SHAGGY-LIKE KINASE2 (GSK2) interacts directly with and phosphorylates PHOSPHATE STARVATION RESPONSE2 (OsPHR2), the key regulator of phosphate (Pi) signaling, to suppress its transcription factor activity in rice (Oryza sativa). We identify a critical phosphorylation site at serine residue S269 of OsPHR2 and demonstrate that phosphorylation by GSK2 or phosphor-mimic mutation of S269 substantially impairs the DNA-binding activity of OsPHR2, and thus diminishes expression of OsPHR2-induced genes and reduces Pi levels. Like BRs, Pi starvation noticeably induces GSK2 instability. We further show that this site-specific phosphorylation event is conserved in Arabidopsis (Arabidopsis thaliana), but varies among the PHR-family members, being present only in most land plants. These results unveil a distinctive post-transcriptional regulatory mechanism in Pi signaling by which BRs promote Pi acquisition, with a potential contribution to the environmental adaptability of plants during their evolution.

Investigation of a Perspective Urban Tree Species, Ginkgo biloba L., by Scientific Analysis of Historical Old Specimens

In this study, we examined over 200-year-old Ginkgo biloba L. specimens under different environmental conditions. The overall aim was to explore which factors influence their vitality and general fitness in urban environments and thus their ability to tolerate stressful habitats. In order to determine this, we used a number of different methods, including histological examinations (stomatal density and size) and physiological measurements (peroxidase enzyme activity), as well as assessing the air pollution tolerance index (APTI). The investigation of the genetic relationships between individuals was performed using flow cytometry and miRNA marker methods. The genetic tests revealed that all individuals are diploid, whereas the lus-miR168 and lus-miR408 markers indicated a kinship relation between them. These results show that the effect of different habitat characteristics can be detected through morphological and physiological responses, thus indicating relatively higher stress values for all studied individuals. A significant correlation can be found between the level of adaptability and the relatedness of the examined individuals. These results suggest that Ginkgo biloba L. is well adapted to an environment with increased stress factors and therefore suitable for use in urban areas.

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 deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid-flavonoid pathway in Thymus mongolicus

BACKGROUND

Thymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.

RESULTS

Chromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid-flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.

CONCLUSION

This study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Decoding the leaf apical meristem of Guarea glabra Vahl (Meliaceae): insight into the evolution of indeterminate pinnate leaves

In seed plants, growth of shoots and roots is indeterminate, while leaves are typically determinate organs that cease to grow after a certain developmental stage. This is due to the characteristics of the leaf meristem, where cell proliferation activity is retained only for a limited period. However, several plants exhibit indeterminacy in their leaves, exemplified by the pinnate compound leaves of Guarea and Chisocheton genera in the Meliaceae family. In these plants, the leaf meristem at the tip of the leaf retains meristematic activity and produces leaflets over years, resulting in a single leaf that resembles a twig. The molecular mechanism underlying the indeterminate leaf meristem of these plants has not been examined. In this research, we used Guarea glabra as a model to investigate the development of indeterminate pinnate leaves. Transcriptome analyses revealed that the gene expression profile in leaf apex tissue differed from that in the shoot apex. However, a class 1 KNOTTED-LIKE HOMEOBOX (KNOX1) gene which is lost in Brassicaceae was highly expressed in both tissues. We established an in situ hybridisation system for this species using Technovit 9100 to analyse the spatial expression patterns of genes. We revealed that the leaf meristematic region of G. glabra expresses KNOX1, LEAFY and ANGUSTIFORIA3 simultaneously, suggesting the involvement of these genes in the indeterminacy of the leaf meristem.

Wood of trees: Cellular structure, molecular formation, and genetic engineering

Wood is an invaluable asset to human society due to its renewable nature, making it suitable for both sustainable energy production and material manufacturing. Additionally, wood derived from forest trees plays a crucial role in sequestering a significant portion of the carbon dioxide fixed during photosynthesis by terrestrial plants. Nevertheless, with the expansion of the global population and ongoing industrialization, forest coverage has been substantially decreased, resulting in significant challenges for wood production and supply. Wood production practices have changed away from natural forests toward plantation forests. Thus, understanding the underlying genetic mechanisms of wood formation is the foundation for developing high-quality, fast-growing plantation trees. Breeding ideal forest trees for wood production using genetic technologies has attracted the interest of many. Tremendous studies have been carried out in recent years on the molecular, genetic, and cell-biological mechanisms of wood formation, and considerable progress and findings have been achieved. These studies and findings indicate enormous possibilities and prospects for tree improvement. This review will outline and assess the cellular and molecular mechanisms of wood formation, as well as studies on genetically improving forest trees, and address future development prospects.

Unveiling the complexity of strigolactones: exploring structural diversity, biosynthesis pathways, and signaling mechanisms

Strigolactone is the collective name for compounds containing a butenolide as a part of their structure, first discovered as compounds that induce seed germination of root parasitic plants. They were later found to be rhizosphere signaling molecules that induce hyphal branching of arbuscular mycorrhizal fungi, and, finally, they emerged as a class of plant hormones. Strigolactones are found in root exudates, where they display a great variability in their chemical structure. Their structure varies among plant species, and multiple strigolactones can exist in one species. Over 30 strigolactones have been identified, yet the chemical structure of the strigolactone that functions as an endogenous hormone and is found in the above-ground parts of plants remains unknown. We discuss our current knowledge of the synthetic pathways of diverse strigolactones and their regulation, as well as recent progress in identifying strigolactones as plant hormones. Strigolactone is perceived by the DWARF14 (D14), receptor, an α/β hydrolase which originated by gene duplication of KARRIKIN INSENSITIVE 2 (KAI2). D14 and KAI2 signaling pathways are partially overlapping paralogous pathways. Progress in understanding the signaling mechanisms mediated by two α/β hydrolase receptors as well as remaining challenges in the field of strigolactone research are reviewed.

Plant secondary metabolites involved in the stress tolerance of long-lived trees

Ancient trees are natural wonders because of their longevity, having lived for hundreds or thousands of years, and their ability to withstand changing environments and a variety of stresses. These long-lived trees have sophisticated defense mechanisms, such as the production of specialized plant metabolites (SPMs). In this review, we provide an overview of the major biotic and abiotic stresses that long-lived trees often face, as well as an analysis of renowned ancient tree species and their unique protective SPMs against environmental stressors. We also discuss the synthesis and accumulation of defensive SPMs induced by environmental factors and endophytes in these trees. Furthermore, we conducted a comparative genomic analysis of 17 long-lived tree species and discovered significant expansions of SPM biosynthesis gene families in these species. Our comprehensive review reveals the crucial role of SPMs in high resistance in long-lived trees, providing a novel natural resource for plant defense, crop improvement and even the pharmaceutical industry.

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.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

PPGR: a comprehensive perennial plant genomes and regulation database

Perennial woody plants hold vital ecological significance, distinguished by their unique traits. While significant progress has been made in their genomic and functional studies, a major challenge persists: the absence of a comprehensive reference platform for collection, integration and in-depth analysis of the vast amount of data. Here, we present PPGR (Resource for Perennial Plant Genomes and Regulation; https://ngdc.cncb.ac.cn/ppgr/) to address this critical gap, by collecting, integrating, analyzing and visualizing genomic, gene regulation and functional data of perennial plants. PPGR currently includes 60 species, 847 million protein-protein/TF (transcription factor)-target interactions, 9016 transcriptome samples under various environmental conditions and genetic backgrounds. Noteworthy is the focus on genes that regulate wood production, seasonal dormancy, terpene biosynthesis and leaf senescence representing a wealth of information derived from experimental data, literature mining, public databases and genomic predictions. Furthermore, PPGR incorporates a range of multi-omics search and analysis tools to facilitate browsing and application of these extensive datasets. PPGR represents a comprehensive and high-quality resource for perennial plants, substantiated by an illustrative case study that demonstrates its capacity in unraveling gene functions and shedding light on potential regulatory processes.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

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.

Approaches to increase the validity of gene family identification using manual homology search tools

Identifying homologs is an important process in the analysis of genetic patterns underlying traits and evolutionary relationships among species. Analysis of gene families is often used to form and support hypotheses on genetic patterns such as gene presence, absence, or functional divergence which underlie traits examined in functional studies. These analyses often require precise identification of all members in a targeted gene family. Manual pipelines where homology search and orthology assignment tools are used separately are the most common approach for identifying small gene families where accurate identification of all members is important. The ability to curate sequences between steps in manual pipelines allows for simple and precise identification of all possible gene family members. However, the validity of such manual pipeline analyses is often decreased by inappropriate approaches to homology searches including too relaxed or stringent statistical thresholds, inappropriate query sequences, homology classification based on sequence similarity alone, and low-quality proteome or genome sequences. In this article, we propose several approaches to mitigate these issues and allow for precise identification of gene family members and support for hypotheses linking genetic patterns to functional traits.

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.

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.

Morphological Characteristics, Ultrastructure, and Chemical Constituents of the Endotesta in Ginkgo (Ginkgo biloba L.)

Ginkgo biloba L. is a tree species of significant economic and ecological importance. Prior studies of the Ginkgo biloba seed coat have predominantly focused on the sarcotesta and sclerotesta, with less attention paid to the endotesta. In this study, the development and formation of Ginkgo endotesta were examined using light microscopy and transmission electron microscopy. The structural properties of the mature endotesta were analyzed using micro-CT imaging and scanning electron microscopy. The results indicate that the endotesta possess a membranous structure primarily originating from the inner bead peridium, a segment of bead core tissue, and the macrospore membrane. The endotesta from the middle constriction line to the chalazal end comprises a single layer with a greyish-white papery structure. In contrast, the endotesta was divided into two inner and two outer layers, from the middle constriction line to the micropylar end. The outer endosperm adheres closely to the sclerotesta, while the inner endosperm adheres to the seed kernel. The surface of the endotesta was irregularly raised, with thicker wax at the chalazal end, whereas the micropylar end demonstrated similar characteristics with thinner wax and tumor layers. The endotesta contained 17 amino acids, 18 fatty acids, 10 trace elements, and 7 vitamins. Overall, its nutritional value was relatively well balanced.

Genome and transcriptome of Selaginella kraussiana reveal evolution of root apical meristems in vascular plants

The evolution of roots allowed vascular plants to adapt to land environments. Fossil evidence indicates that roots evolved independently in euphyllophytes (ferns and seed plants) and lycophytes, the two lineages of extant vascular plants. Based on a high-quality genome assembly, mRNA sequencing (mRNA-seq) data, and single-cell RNA-seq data for the lycophyte Selaginella kraussiana, we show that the two root origin events in lycophytes and euphyllophytes adopted partially similar molecular modules in the regulation of root apical meristem (RAM) development. In S. kraussiana, the RAM initiates from the rhizophore primordium guided by auxin and duplicates itself by dichotomous branching. The auxin signaling pathway directly upregulates euAINTEGUMENTAb (SkeuANTb), and then SkeuANTb directly promotes the expression of SkeuANTa and the WUSCHEL-RELATED HOMEOBOX13b (SkWOX13b) for RAM maintenance, partially similar to the molecular pathway involving the euANT-branch PLETHORA (AtPLT) genes and AtWOX5 in root initiation in the seed plant Arabidopsis thaliana. Other molecular modules, e.g., SHORT-ROOT and SCARECROW, also have partially similar expression patterns in the RAMs of S. kraussiana and A. thaliana. Overall, our study not only provides genome and transcriptome tools of S. kraussiana but also indicates the employment of some common molecular modules in RAMs during root origins in lycophytes and euphyllophytes.

Large-Scale Identification and Characterization Analysis of VQ Family Genes in Plants, Especially Gymnosperms

VQ motif-containing (VQ) proteins are a class of transcription regulatory cofactors widely present in plants, playing crucial roles in growth and development, stress response, and defense. Although there have been some reports on the member identification and functional research of VQ genes in some plants, there is still a lack of large-scale identification and clear graphical presentation of their basic characterization information to help us to better understand this family. Especially in gymnosperms, the VQ family genes and their evolutionary relationships have not yet been reported. In this study, we systematically identified 2469 VQ genes from 56 plant species, including bryophytes, gymnosperms, and angiosperms, and analyzed their molecular and evolutionary features. We found that amino acids are only highly conserved in the VQ domain, while other positions are relatively variable; most VQ genes encode relatively small proteins and do not have introns. The GC content in Poaceae plants is the highest (up to 70%); these VQ proteins can be divided into nine subgroups. In particular, we analyzed the molecular characteristics, chromosome distribution, duplication events, and expression levels of VQ genes in three gymnosperms: Ginkgo biloba, Taxus chinensis, and Pinus tabuliformis. In gymnosperms, VQ genes are classified into 11 groups, with highly similar motifs in each group; most VQ proteins have less than 300 amino acids and are predicted to be located in nucleus. Tandem duplication is an important driving force for the expansion of the VQ gene family, and the evolutionary processes of most VQ genes and duplication events are relatively independent; some candidate VQ genes are preliminarily screened, and they are likely to be involved in plant growth and stress and defense responses. These results provide detailed information and powerful references for further understanding and utilizing the VQ family genes in various plants.

GAEP: a comprehensive genome assembly evaluating pipeline

With the rapid development of sequencing technologies, especially the maturity of third-generation sequencing technologies, there has been a significant increase in the number and quality of published genome assemblies. The emergence of these high-quality genomes has raised higher requirements for genome evaluation. Although numerous computational methods have been developed to evaluate assembly quality from various perspectives, the selective use of these evaluation methods can be arbitrary and inconvenient for fairly comparing the assembly quality. To address this issue, we have developed the Genome Assembly Evaluating Pipeline (GAEP), which provides a comprehensive assessment pipeline for evaluating genome quality from multiple perspectives, including continuity, completeness, and correctness. Additionally, GAEP includes new functions for detecting misassemblies and evaluating the assembly redundancy, which performs well in our testing. GAEP is publicly available at https://github.com/zy-optimistic/GAEP under the GPL3.0 License. With GAEP, users can quickly obtain accurate and reliable evaluation results, facilitating the comparison and selection of high-quality genome assemblies.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

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.

Screening of Ginkgo Individuals with Superior Growth Structural Characteristics in Different Genetic Groups Using Terrestrial Laser Scanning (TLS) Data

With the concept of sustainable management of plantations, individual trees with excellent characteristics in plantations have received attention from breeders. To improve and maintain long-term productivity, accurate and high-throughput access to phenotypic characteristics is essential when establishing breeding strategies. Meanwhile, genetic diversity is also an important issue that must be considered, especially for plantations without seed source information. This study was carried out in a ginkgo timber plantation. We used simple sequence repeat (SSR) markers for genetic background analysis and high-density terrestrial laser scanning for growth structural characteristic extraction, aiming to provide a possibility of applying remote sensing approaches for forest breeding. First, we analyzed the genetic diversity and population structure, and grouped individual trees according to the genetic distance. Then, the growth structural characteristics (height, diameter at breast height, crown width, crown area, crown volume, height to living crown, trunk volume, biomass of all components) were extracted. Finally, individual trees in each group were comprehensively evaluated and the best-performing ones were selected. Results illustrate that terrestrial laser scanning (TLS) point cloud data can provide nondestructive estimates of the growth structural characteristics at fine scale. From the ginkgo plantation containing high genetic diversity (average polymorphism information content index was 0.719) and high variation in growth structural characteristics (coefficient of variation ranged from 21.822% to 85.477%), 11 excellent individual trees with superior growth were determined. Our study guides the scientific management of plantations and also provides a potential for applying remote sensing technologies to accelerate forest breeding.

Identification and Expression of the MADS-box Gene Family in Different Versions of the Ginkgo biloba Genome

MADS-box transcription factors play important roles in many organisms. These transcription factors are involved in processes such as the formation of the flower organ structure and the seed development of plants. Ginkgo biloba has two genome versions (version 2019 and version 2021), and there is no analysis or comparison of the MADS-box gene family in these two genomes. In this study, 26 and 20 MADS-box genes were identified from the two genomes of Ginkgo, of which 12 pairs of genes reached more than 80% similarity. According to our phylogenetic analysis results, we divided these genes into type I (Mα and Mγ subfamilies) and type II (MIKC and Mδ subfamilies) members. We found that both sets of genomes lacked the Mβ gene, while the MIKC gene was the most numerous. Further analysis of the gene structure showed that the MIKC genes in the two genomes had extralong introns (≥20 kb); these introns had different splicing patterns, and their expression might be more abundant. The gene expression analysis proved that GbMADS genes were expressed to varying degrees in eight Ginkgo biological tissues. Type II GbMADS genes not only were found to be related to female flower bud differentiation and development but also are important in seed development. Therefore, MADS-box genes may play important roles in the development of Ginkgo reproductive organs, which may suggest a genetic role in sexual differentiation. This study further contributes to the research on MADS-box genes and provides new insights into sex determination in Ginkgo.

Genomic and Transcriptomic Insights into the Evolution and Divergence of MIKC-Type MADS-Box Genes in Carica papaya

MIKC-type MADS-box genes, also known as type II genes, play a crucial role in regulating the formation of floral organs and reproductive development in plants. However, the genome-wide identification and characterization of type II genes as well as a transcriptomic survey of their potential roles in Carica papaya remain unresolved. Here, we identified and characterized 24 type II genes in the C. papaya genome, and investigated their evolutional scenario and potential roles with a widespread expression profile. The type II genes were divided into thirteen subclades, and gene loss events likely occurred in papaya, as evidenced by the contracted member size of most subclades. Gene duplication mainly contributed to MIKC-type gene formation in papaya, and the duplicated gene pairs displayed prevalent expression divergence, implying the evolutionary significance of gene duplication in shaping the diversity of type II genes in papaya. A large-scale transcriptome analysis of 152 samples indicated that different subclasses of these genes showed distinct expression patterns in various tissues, biotic stress response, and abiotic stress response, reflecting their divergent functions. The hub-network of male and female flowers and qRT-PCR suggested that TT16-3 and AGL8 participated in male flower development and seed germination. Overall, this study provides valuable insights into the evolution and functions of MIKC-type genes in C. papaya.

Genome-wide identification and bioinformatics analysis of the WD40 transcription factor family and candidate gene screening for anthocyanin biosynthesis in Rhododendron simsii

WD40 transcription factors (TFs) constitute a large gene family in eukaryotes, playing diverse roles in cellular processes. However, their functions in the major ornamental plant, Rhododendron simsii, remain poorly understood. In this study, we identified 258 WD40 proteins in the R. simsii genome, which exhibited an uneven distribution across chromosomes. Based on domain compositions and phylogenetic analysis, we classified these 258 RsWD40 proteins into 42 subfamilies and 47 clusters. Comparative genomic analysis suggested that the expansion of the WD40 gene family predates the divergence of green algae and higher plants, indicating an ancient origin. Furthermore, by analyzing the duplication patterns of RsWD40 genes, we found that transposed duplication played a major role in their expansion. Notably, the majority of RsWD40 gene duplication pairs underwent purifying selection during evolution. Synteny analysis identified significant orthologous gene pairs between R. simsii and Arabidopsis thaliana, Oryza sativa, Vitis vinifera, and Malus domestica. We also investigated potential candidate genes involved in anthocyanin biosynthesis during different flower development stages in R. simsii using RNA-seq data. Specifically, we identified 10 candidate genes during the bud stage and 7 candidate genes during the full bloom stage. GO enrichment analysis of these candidate genes revealed the potential involvement of the ubiquitination process in anthocyanin biosynthesis. Overall, our findings provide a valuable foundation for further investigation and functional analysis of WD40 genes, as well as research on the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species.

The genome of the Wollemi pine, a critically endangered "living fossil" unchanged since the Cretaceous, reveals extensive ancient transposon activity

We present the genome of the living fossil, Wollemia nobilis, a southern hemisphere conifer morphologically unchanged since the Cretaceous. Presumed extinct until rediscovery in 1994, the Wollemi pine is critically endangered with less than 60 wild adults threatened by intensifying bushfires in the Blue Mountains of Australia. The 12 Gb genome is among the most contiguous large plant genomes assembled, with extremely low heterozygosity and unusual abundance of DNA transposons. Reduced representation and genome re-sequencing of individuals confirms a relictual population since the last major glacial/drying period in Australia, 120 ky BP. Small RNA and methylome sequencing reveal conservation of ancient silencing mechanisms despite the presence of thousands of active and abundant transposons, including some transferred horizontally to conifers from arthropods in the Jurassic. A retrotransposon burst 8-6 my BP coincided with population decline, possibly as an adaptation enhancing epigenetic diversity. Wollemia, like other conifers, is susceptible to Phytophthora, and a suite of defense genes, similar to those in loblolly pine, are targeted for silencing by sRNAs in leaves. The genome provides insight into the earliest seed plants, while enabling conservation efforts.