Yellow-green leaf mutants are common in higher plants, and these non-lethal chlorophyll-deficient mutants are ideal materials for research on photosynthesis and plant development. A novel xantha mutant of Ginkgo biloba displaying yellow-colour leaves (YL) and green-colour leaves (GL) was identified in this study. The chlorophyll content of YL was remarkably lower than that in GL. The chloroplast ultrastructure revealed that YL had less dense thylakoid lamellae, a looser structure and fewer starch grains than GL. Analysis of the photosynthetic characteristics revealed that YL had decreased photosynthetic activity with significantly high nonphotochemical quenching. To explain these phenomena, we analysed the proteomic differences in leaves and chloroplasts between YL and GL of ginkgo using two-dimensional gel electrophoresis (2-DE) coupled with MALDI-TOF/TOF MS. In total, 89 differential proteins were successfully identified, 82 of which were assigned functions in nine metabolic pathways and cellular processes. Among them, proteins involved in photosynthesis, carbon fixation in photosynthetic organisms, carbohydrate/energy metabolism, amino acid metabolism, and protein metabolism were greatly enriched, indicating a good correlation between differentially accumulated proteins and physiological changes in leaves. The identifications of these differentially accumulated proteins indicates the presence of a specific different metabolic network in YL and suggests that YL possess slower chloroplast development, weaker photosynthesis, and a less abundant energy supply than GL. These studies provide insights into the mechanism of molecular regulation of leaf colour variation in YL mutants.

3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), as the rate-limiting enzyme in the mevalonate pathway, is essential for the biosynthesis of shikonin in Lithospermum erythrorhizon. However, in the absence of sufficient data, the principles of a genome-wide in-depth evolutionary exploration of HMGR family members in plants, as well as key members related to shikonin biosynthesis, remain unidentified. In this study, 124 HMGRs were identified and characterized from 36 representative plants, including L. erythrorhizon. Vascular plants were found to have more HMGR family genes than nonvascular plants. The phylogenetic tree revealed that during lineage and species diversification, the HMGRs evolved independently and intronless LerHMGRs emerged from multi-intron HMGR in land plants. Among them, Pinus tabuliformis and L. erythrorhizon had the most HMGR gene duplications, with 11 LerHMGRs most likely expanded through WGD/segmental and tandem duplications. In seedling roots and M9 cultured cells/hairy roots, where shikonin biosynthesis occurs, LerHMGR1 and LerHMGR2 were expressed significantly more than other genes. The enzymatic activities of LerHMGR1 and LerHMGR2 further supported their roles in catalyzing the conversion of HMG-CoA to mevalonate. Our findings provide insight into the molecular evolutionary properties and function of the HMGR family in plants and a basis for the genetic improvement of efficiently produced secondary metabolites in L. erythrorhizon.

To elucidate the transcriptomic regulation mechanisms that underlie the response of Ginkgo biloba to dehydration and rehydration, we used ginkgo saplings exposed to osmotically driven water stress and subsequent rewatering. When compared with a control group, 137, 1453, 1148, and 679 genes were differentially expressed in ginkgo leaves responding to 2, 6, 12, and 24 h of water deficit, and 796 and 1530 genes were differentially expressed responding to 24 and 48 h of rewatering. Upregulated genes participated in the biosynthesis of abscisic acid, eliminating reactive oxygen species (ROS), and biosynthesis of flavonoids and bilobalide, and downregulated genes were involved in water transport and cell wall enlargement in water stress-treated ginkgo leaves. Under rehydration conditions, the genes associated with water transport and cell wall enlargement were upregulated, and the genes that participated in eliminating ROS and the biosynthesis of flavonoids and bilobalide were downregulated in the leaves of G. biloba. Furthermore, the weighted gene coexpression networks were established and correlated with distinct water stress and rewatering time-point samples. Hub genes that act as key players in the networks were identified. Overall, these results indicate that the gene coexpression networks play essential roles in the transcriptional reconfiguration of ginkgo leaves in response to water stress and rewatering.

Melia azedarach L. is a native tree species that can be used in a comprehensive way and is widely distributed in all provinces south of the Yellow River in China. Genetic diversity analysis of different M. azedarach germplasm sources is an important basic work for the selection, evaluation, and genetic improvement of M. azedarach germplasm resources. In this study, 100 pairs of SSR primers were designed and synthesized based on M. azedarach transcriptome data, and 16 pairs of reliable SSR primers were finally selected. The developed primers were used to analyze the genetic diversity of M. azedarach from 15 sources in 10 provinces in East, Central, and South China. The results showed that the frequency of the M. azedarach transcriptome SSR loci was high, and the distribution density was high. There were 15 sources of M. azedarach genetic diversity at a moderate level, and genetic variation was mainly present within the sources. The present study further enriches the existing SSR marker database of the M. azedarach family and can provide a reference for genetic diversity analysis and molecularly assisted breeding of M. azedarach plants at the genomic level.

DNA methylation plays a vital role in diverse biological processes. DNA methyltransferases (DNMTs) genes and RNA-directed DNA methylation (RdDM)-related genes are key genes responsible for establishing and maintaining genome DNA methylation in plants. In the present study, we systematically identified nine GbDNMTs in Ginkgo biloba, including the three common families of GbMET1a/1b, GbCMT2, and GbDRMa/b/2a/2b/2c, and a fourth family—GbDNMT3—which is absent in most angiosperms. We also identified twenty RdDM-related genes, including four GbDCLs, six GbAGOs, and ten GbRDRs. Expression analysis of the genes showed the different patterns of individual genes, and 15 of 29 genes displayed expression change under five types of abiotic stress. Gene coexpression analysis and weighted gene co-expression network analysis (WGCNA) using 126 public transcriptomic datasets revealed that these genes were clustered into two groups. In group I, genes covered members from all six families which were preferentially expressed in the ovulate strobile and fruit. A gene ontology (GO) enrichment analysis of WGCNA modules indicated that group I genes were most correlated with the biological process of cell proliferation. Group II only consisted of RdDM-related genes, including GbDRMs, GbAGOs, and GbRDRs, but no GbDCLs, and these genes were specifically expressed in the cambium, suggesting that they may function in a dicer-like (DCL)-independent RdDM pathway in specific tissues. The gene module related to group II was most enriched in signal transduction, cell communication, and the response to the stimulus. These results demonstrate that gene family members could be conserved or diverged across species, and multi-member families in the same pathway may cluster into different modules to function differentially. The study provides insight into the DNA methylation genes and their possible functions in G. biloba, laying a foundation for the further study of DNA methylation in gymnosperms.

The BAHD acyltransferase family is a unique class of plant proteins that acylates plant metabolites and participates in plant secondary metabolic processes. However, the BAHD members in Lithospermum erythrorhizon remain unknown and uncharacterized. Although the heterologously expressed L. erythrorhizon BAHD family member LeSAT1 in Escherichia coli has been shown to catalyze the conversion of shikonin to acetylshikonin in vitro, its in vivo role remains unknown. In this study, the characterization, evolution, expression patterns, and gene function of LeBAHDs in L. erythrorhizon were explored by bioinformatics and transgenic analysis. We totally identified 73 LeBAHDs in the reference genome of L. erythrorhizon. All LeBAHDs were phylogenetically classified into five clades likely to perform different functions, and were mainly expanded by dispersed and WGD/segmental duplication. The in vivo functional investigation of the key member LeBAHD1/LeSAT1 revealed that overexpression of LeBAHD1 in hairy roots significantly increased the content of acetylshikonin as well as the conversion rate of shikonin to acetylshikonin, whereas the CRISPR/Cas9-based knockout of LeBAHD1 in hairy roots displayed the opposite trend. Our results not only confirm the in vivo function of LeBAHD1/LeSAT1 in the biosynthesis of acetylshikonin, but also provide new insights for the biosynthetic pathway of shikonin and its derivatives.