Growth regulatory factors (GRF) are plant-specific transcription factors that play pivotal roles in growth and various abiotic stresses regulation. However, adaptive evolution of GRF gene family in land plants are still being elucidated. Here, we performed the evolutionary and expression analysis of GRF gene family from seven representative species. Extensive phylogenetic analyses and gene structure analysis revealed that the number of genes, QLQ domain and WRC domain identified in higher plants was significantly greater than those identified in lower plants. Besides, dispersed duplication and WGD/segmental duplication effectively promoted expansion of the GRF gene family. The expression patterns of GRF gene family and target genes were found in multiple floral organs and abundant in actively growing tissues. They were also found to be particularly expressed in response to various abiotic stresses, with stress-related elements in promoters, implying potential roles in floral development and abiotic stress. Our analysis in GRF gene family interaction network indicated the similar results that GRFs resist to abiotic stresses with the cooperation of other transcription factors like GIFs. This study provides insights into evolution in the GRF gene family, together with expression patterns valuable for future functional researches of plant abiotic stress biology.

(Short Internodes/Stylish/SHI-Related Sequence) genes play a crucial role in plant developments, encompassing organ morphogenesis. However, at present, the biological significance of potato StSRS genes remains unknown. This study comprehensively identified and analysed the potential functions of StSRS genes in potato tuberization. Through phylogenetic and conserved motif analyses, eight members of the StSRS family were successfully identified and classified into five distinct subfamilies. The intraspecies and interspecies collinearity analyses offered insights into the evolutionary clues of the StSRS family, showing closer homology to tomato and tobacco. High-quality RNA-sequencing (Q30 > 97.16%) revealed dynamic gene expression profiles, with 146 to 997 differentially expressed transcription factors (DETFs), particularly from the MYB-related, GATA, and bHLH families. Among them, the expression of two StSRS genes, StSRS1 and StSRS8, were dramatically changed during potato tuber formation. Weighted Gene Co-expression Network Analysis (WGCNA) uncovered stage-specific gene modules and spotlighted hub genes. Specifically, StSRS8 was strongly associated with the Stolon stage, while StSRS1 was evidently linked to the Tuber stage. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis revealed that the transcripts of both StSRS1 and StSRS8 were predominantly expressed during potato tuberization, as compared to those in other tissues. However, it is scarcely detectable for other six StSRS genes. Furthermore, the transcripts of StSRS8 were repressed under short-day (SD) conditions. These findings yield critical insights into the potato StSRS gene family and furnish essential information for further function investigation of StSRS genes, particularly StSRS1 and StSRS8, in the context of potato tuberization.

Background: As the characteristic functional component in ginger, gingerols possess several health-promoting properties. Long non-coding RNAs (lncRNAs) act as crucial regulators of diverse biological processes. However, lncRNAs in ginger are not yet identified so far, and their potential roles in gingerol biosynthesis are still unknown. In this study, metabolomic and transcriptomic analyses were performed in three main ginger cultivars (leshanhuangjiang, tonglingbaijiang, and yujiang 1 hao) in China to understand the potential roles of the specific lncRNAs in gingerol accumulation. Results: A total of 744 metabolites were monitored by metabolomics analysis, which were divided into eleven categories. Among them, the largest group phenolic acid category contained 143 metabolites, including 21 gingerol derivatives. Of which, three gingerol analogs, [8]-shogaol, [10]-gingerol, and [12]-shogaol, accumulated significantly. Moreover, 16,346 lncRNAs, including 2,513, 1,225, and 2,884 differentially expressed (DE) lncRNA genes (DELs), were identified in all three comparisons by transcriptomic analysis. Gene ontology enrichment (GO) analysis showed that the DELs mainly enriched in the secondary metabolite biosynthetic process, response to plant hormones, and phenol-containing compound metabolic process. Correlation analysis revealed that the expression levels of 11 DE gingerol biosynthesis enzyme genes (GBEGs) and 190 transcription factor genes (TF genes), such as MYB1, ERF100, WRKY40, etc. were strongly correlation coefficient with the contents of the three gingerol analogs. Furthermore, 7 and 111 upstream cis-acting lncRNAs, 1,200 and 2,225 upstream trans-acting lncRNAs corresponding to the GBEGs and TF genes were identified, respectively. Interestingly, 1,184 DELs might function as common upstream regulators to these GBEGs and TFs genes, such as LNC_008452, LNC_006109, LNC_004340, etc. Furthermore, protein–protein interaction networks (PPI) analysis indicated that three TF proteins, MYB4, MYB43, and WRKY70 might interact with four GBEG proteins (PAL1, PAL2, PAL3, and 4CL-4). Conclusion: Based on these findings, we for the first time worldwide proposed a putative regulatory cascade of lncRNAs, TFs genes, and GBEGs involved in controlling of gingerol biosynthesis. These results not only provide novel insights into the lncRNAs involved in gingerol metabolism, but also lay a foundation for future in-depth studies of the related molecular mechanism.

Herba Epimedii’s leaves are highly valued in traditional Chinese medicine for their substantial concentration of flavonoids, which play a crucial role in manifesting the plant’s therapeutic properties. This study investigated the metabolomic, transcriptomic and proteomic profiles of leaves from two Herba Epimedii cultivars, Epipremnum sagittatum (J) and Epipremnum pubescens (R), at three different developmental stages. Metabolite identification and analysis revealed a total of 1,412 and 1,421 metabolites with known structures were found. Flavonoids made up of 33%, including 10 significant accumulated icariin analogues. Transcriptomic analysis unveiled totally 41,644 differentially expressed genes (DEGs) containing five encoded genes participated in icariin biosynthesis pathways. Totally, 9,745 differentially expressed proteins (DEPs) were found, including Cluster-47248.2.p1 (UDP-glucuronosy/UDP-glucosyltransferase), Cluster-30441.2.p1 (O-glucosyltransferase), and Cluster-28344.9.p1 (anthocyanidin 3-O-glucoside 2 “-O-glucosyltransferase-like) through proteomics analysis which are involved to icariin biosynthesis. Protein-protein interaction (PPI) assay exhibited, totally 12 proteins showing a strong relationship of false discovery rate (FDR) <0.05 with these three proteins containing 2 leucine-rich repeat receptor kinase-like protein SRF7, and 5 methyl jasmonate esterase 1. Multi-omics connection networks uncovered 237 DEGs and 72 DEPs exhibited significant associations with the 10 icariin analogues. Overall, our integrated omics approach provides comprehensive insights into the regulatory network underlying icariin synthesis in Herba Epimedii, offering valuable resources for further research and development in medicinal plant cultivation and pharmaceutical applications.