Structural variation, functional differentiation and expression characteristics of the AP2/ERF gene family and its response to cold stress and methyl jasmonate in Panax ginseng C.A. Meyer

Transcription factor PgNAC72 activates DAMMARENEDIOL SYNTHASE expression to promote ginseng saponin biosynthesis

Ginsenosides, the primary bioactive constituents in ginseng (Panax ginseng), possess substantial pharmacological potential and are in high demand in the market. The plant hormone methyl jasmonate (MeJA) effectively elicits ginsenoside biosynthesis in P. ginseng, though the regulatory mechanism remains largely unexplored. NAC transcription factors are critical in intricate plant regulatory networks and participate in numerous plant physiological activities. In this study, we identified a MeJA-responsive NAC transcription factor gene, PgNAC72, from a transcriptome library produced from MeJA-treated P. ginseng callus. Predominantly expressed in P. ginseng flowers, PgNAC72 localizes to the nucleus. Overexpressing PgNAC72 (OE-PgNAC72) in P. ginseng callus notably elevated total saponin levels, particularly dammarane-type ginsenosides, by upregulating dammarenediol synthase (PgDDS), encoding a key enzyme in the ginsenoside biosynthesis pathway. Electrophoretic mobility shift assays and dual-luciferase assays confirmed that PgNAC72 binds to the NAC-binding elements in the PgDDS promoter, thereby activating its transcription. Further RNA-seq and terpenoid metabolomic data in the OE-PgNAC72 line confirmed that PgNAC72 enhances ginsenoside biosynthesis. These findings uncover a regulatory role of PgNAC72 in MeJA-mediated ginsenoside biosynthesis, providing insights into the ginsenoside regulatory network and presenting a valuable target gene for metabolic engineering.

Beyond genome: Advanced omics progress of Panax ginseng

Ginseng is a perennial herb of the genus Panax in the family Araliaceae as one of the most important traditional medicine. Genomic studies of ginseng assist in the systematic discovery of genes related to bioactive ginsenosides biosynthesis and resistance to stress, which are of great significance in the conservation of genetic resources and variety improvement. The transcriptome reflects the difference and consistency of gene expression, and transcriptomics studies of ginseng assist in screening ginseng differentially expressed genes to further explore the powerful gene source of ginseng. Protein is the ultimate bearer of ginseng life activities, and proteomic studies of ginseng assist in exploring the biosynthesis and regulation of secondary metabolites like ginsenosides and the molecular mechanism of ginseng adversity adaptation at the overall level. In this review, we summarize the current status of ginseng research in genomics, transcriptomics and proteomics, respectively. We also discuss and look forward to the development of ginseng genome allele mapping, ginseng spatiotemporal, single-cell transcriptome, as well as ginseng post-translational modification proteome. We hope that this review will contribute to the in-depth study of ginseng and provide a reference for future analysis of ginseng from a systems biology perspective.

BnAP2-12 overexpression delays ramie flowering: evidence from AP2/ERF gene expression

Introduction

The APETALA2/ethylene response factor (AP2/ERF) superfamily plays a significant role in regulating plant gene expression in response to growth and development. To date, there have been no studies into whether the ramie AP2/ERF genes are involved in the regulation of flower development.

Methods

Here, 84 BnAP2/ERF members were identified from the ramie genome database, and various bioinformatics data on the AP2/ERF gene family, structure, replication, promoters and regulatory networks were analysed. BnAP2-12 was transferred into Arabidopsis through the flower-dipping method.

Results

Phylogenetic analysis classified the 84 BnAP2/ERF members into four subfamilies: AP2 (18), RAV (3), ERF (42), and DREB (21). The functional domain analysis of genes revealed 10 conserved motifs. Genetic mapping localised the 84 members on 14 chromosomes, among which chromosomes 1, 3, 5, and 8 had more members. Collinearity analysis revealed that 43.37% possibly resulted from replication events during the evolution of the ramie genome. Promoter sequence analysis identified classified cis-acting elements associated with plant growth and development, and responses to stress, hormones, and light. Transcriptomic comparison identified 3,635 differentially expressed genes (DEGs) between male and female flowers (1,803 and 1,832 upregulated and downregulated genes, respectively). Kyoto Encyclopaedia of Genes and Genomes pathway analysis categorised DEGs involved in metabolic pathways and biosynthesis of secondary metabolites. Gene Ontology enrichment analysis further identified enriched genes associated with pollen and female gamete formations. Of the 84 BnAP2/ERFs genes, 22 and 8 upregulated and downregulated genes, respectively, were present in female flowers. Co-expression network analysis identified AP2/ERF members associated with flower development, including BnAP2-12. Subcellular localisation analysis showed that the BnAP2-12 protein is localised in the nucleus and cell membrane. Overexpression BnAP2-12 delayed the flowering time of Arabidopsis thaliana.

Conclusion

These findings provide insights into the mechanism of ramie flower development.

Molecular regulation of the key specialized metabolism pathways in medicinal plants

The basis of modern pharmacology is the human ability to exploit the production of specialized metabolites from medical plants, for example, terpenoids, alkaloids, and phenolic acids. However, in most cases, the availability of these valuable compounds is limited by cellular or organelle barriers or spatio-temporal accumulation patterns within different plant tissues. Transcription factors (TFs) regulate biosynthesis of these specialized metabolites by tightly controlling the expression of biosynthetic genes. Cutting-edge technologies and/or combining multiple strategies and approaches have been applied to elucidate the role of TFs. In this review, we focus on recent progress in the transcription regulation mechanism of representative high-value products and describe the transcriptional regulatory network, and future perspectives are discussed, which will help develop high-yield plant resources.

Jasmonate activates secondary cell wall biosynthesis through MYC2-MYB46 module

Formation of secondary cell wall (SCW) is tightly regulated spatiotemporally by various developmental and environmental signals. Successful fine-tuning of the trade-off between SCW biosynthesis and stress responses requires a better understanding of how plant growth is regulated under environmental stress conditions. However, the current understanding of the interplay between environmental signaling and SCW formation is limited. The lipid-derived plant hormone jasmonate (JA) and its derivatives are important signaling components involved in various physiological processes including plant growth, development, and abiotic/biotic stress responses. Recent studies suggest that JA is involved in SCW formation but the signaling pathway has not been studied for how JA regulates SCW formation. We tested this hypothesis using the transcription factor MYB46, a master switch for SCW biosynthesis, and JA treatments. Both the transcript and protein levels of MYB46, a master switch for SCW formation, were significantly increased by JA treatment, resulting in the upregulation of SCW biosynthesis. We then show that this JA-induced upregulation of MYB46 is mediated by MYC2, a central regulator of JA signaling, which binds to the promoter of MYB46. We conclude that this MYC2-MYB46 module is a key component of the plant response to JA in SCW formation.

Bioinformatic Assessment and Expression Profiles of the AP2/ERF Superfamily in the Melastoma dodecandrum Genome

AP2/ERF transcription factors play crucial roles in various biological activities, including plant growth, development, and responses to biotic and abiotic stressors. However, limited research has been conducted on the AP2/ERF genes of Melastoma dodecandrum for breeding of this potential fruit crop. Leveraging the recently published whole genome sequence, we conducted a comprehensive assessment of this superfamily and explored the expression patterns of AP2/ERF genes at a genome-wide level. A significant number of genes, totaling 218, were discovered to possess the AP2 domain sequence and displayed notable structural variations among five subfamilies. An uneven distribution of these genes was observed on 12 pseudochromosomes as the result of gene expansion facilitated by segmental duplications. Analysis of cis-acting elements within promoter sites and 87.6% miRNA splicing genes predicted their involvement in multiple hormone responses and abiotic stresses through transcriptional and post-transcriptional regulations. Transcriptome analysis combined with qRT-PCR results indicated that certain candidate genes are involved in tissue formation and the response to developmental changes induced by IAA hormones. Overall, our study provides valuable insights into the evolution of ERF genes in angiosperms and lays a solid foundation for future breeding investigations aimed at improving fruit quality and enhancing adaptation to barren land environments.

Functional Characterization of AP2/ERF Transcription Factors during Flower Development and Anthocyanin Biosynthesis Related Candidate Genes in Lycoris

The APETALA2/ethylene-responsive transcription factor (AP2/ERF) family has been extensively investigated because of its significant involvement in plant development, growth, fruit ripening, metabolism, and plant stress responses. To date, there has been little investigation into how the AP2/ERF genes influence flower formation and anthocyanin biosynthesis in Lycoris. Herein, 80 putative LrAP2/ERF transcription factors (TFs) with complete open reading frames (ORFs) were retrieved from the Lycoris transcriptome sequence data, which could be divided into five subfamilies dependent on their complete protein sequences. Furthermore, our findings demonstrated that genes belonging to the same subfamily had structural similarities and conserved motifs. LrAP2/ERF genes were analyzed for playing an important role in plant growth, water deprivation, and flower formation by means of gene ontology (GO) enrichment analysis. The expression pattern of the LrAP2/ERF genes differed across tissues and might be important for Lycoris growth and flower development. In response to methyl jasmonate (MeJA) exposure and drought stress, the expression of each LrAP2/ERF gene varied across tissues and time. Moreover, a total of 20 anthocyanin components were characterized using ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis, and pelargonidin-3-O-glucoside-5-O-arabinoside was identified as the major anthocyanin aglycone responsible for the coloration of the red petals in Lycoris. In addition, we mapped the relationships between genes and metabolites and found that LrAP2/ERF16 is strongly linked to pelargonidin accumulation in Lycoris petals. These findings provide the basic conceptual groundwork for future research into the molecular underpinnings and regulation mechanisms of AP2/ERF TFs in anthocyanin accumulation and Lycoris floral development.

Genome-wide characterization, evolutionary analysis, and expression pattern analysis of the trihelix transcription factor family and gene expression analysis under MeJA treatment in Panax ginseng

Panax ginseng is a well-known medicinal plant with several pharmacological uses in China. The trihelix family transcription factors, also known as GT factors, can be involved in the regulation of growth and developmental processes in plants. There have been no in-depth reports or systematic studies about the trihelix transcription factor in ginseng. In this study, the structure, chromosomal localization, gene duplication, phylogeny, functional differentiation, expression patterns and coexpression interactions of trihelix transcripts were analysed using bioinformatics methods based on the ginseng transcriptome database. Thirty-two trihelix transcription factor genes were identified in ginseng, and these genes were alternatively spliced to obtain 218 transcripts. These transcripts were unevenly distributed on different chromosomes of ginseng, and phylogenetic analysis classified the PgGT transcripts into five subgroups. Gene Ontology (GO) analysis classified PgGT transcripts into eight functional subclasses, indicating that they are functionally diverse. The expression pattern analysis of 218 PgGT transcripts revealed that their expression was tissue-specific and spatiotemporally-specific in 14 different tissues of 4-year-old ginseng, 4 different ages of ginseng roots, and 42 farmers' cultivars of 4-year-old ginseng roots. Despite the differences in the expression patterns of these transcripts, coexpression network analysis revealed that these transcripts could be expressed synergistically in ginseng. In addition, two randomly selected PgGT transcripts in each of the five different subfamilies were subjected to methyl jasmonate treatment at different times, and PgGT was able to respond to the regulation of methy1 jasmonate. These results provide a theoretical basis and gene resources for an in-depth study of the function of trihelix genes in other plants.

Identification of the AP2/ERF transcription factor family of Eleutherococcus senticosus and its expression correlation with drought stress

In this study, through analysis of the genome of Eleutherococcus senticosus (ES). 228 AP2/ERF genes were identified and classified into 5 groups AP2 (47 genes), ERF (108 genes), RAV (6 genes), DREB (64 genes), and soloist (3 genes). According to the AP2/ERF classification of Arabidopsis thaliana, the ES AP2/ERF proteins were subdivided into 15 groups. The gene structure and motifs of each group of AP2/ERF in ES were highly similar, which confirmed the conservation of AP2/ERF genes. The ES AP2/ERF genes were unevenly distributed on chromosomes, and a total of four pairs of tandem repeats, and 84 co-linear gene pairs were found, so the AP2/ERF genes expanded in a fragment replication manner, and dominated by pure selection during evolution. By analyzing the transcriptome data of ES under different drought stress conditions, 87 AP2/ERF genes with differential expression were obtained, of which 10 genes with highly significant differences were further analyzed and screened for qRT-PCR validation. To the best of our knowledge, this is the first report on the AP2/ERF gene of Eleutherococcus senticosus, and the bioinformatics analysis and experimental validation provided valuable information about them, which is of great significance for further research on the molecular mechanisms of ES in response to drought stress.

Multilayered regulation of secondary metabolism in medicinal plants

Medicinal plants represent a huge reservoir of secondary metabolites (SMs), substances with significant pharmaceutical and industrial potential. However, obtaining secondary metabolites remains a challenge due to their low-yield accumulation in medicinal plants; moreover, these secondary metabolites are produced through tightly coordinated pathways involving many spatiotemporally and environmentally regulated steps. The first regulatory layer involves a complex network of transcription factors; a second, more recently discovered layer of complexity in the regulation of SMs is epigenetic modification, such as DNA methylation, histone modification and small RNA-based mechanisms, which can jointly or separately influence secondary metabolites by regulating gene expression. Here, we summarize the findings in the fields of genetic and epigenetic regulation with a special emphasis on SMs in medicinal plants, providing a new perspective on the multiple layers of regulation of gene expression.

Transcriptome-Wide Integrated Analysis of the PgGT25-04 Gene in Controlling Ginsenoside Biosynthesis in Panax ginseng

Panax ginseng is a valuable medicinal herb of the Araliaceae family with various pharmacological activities. The Trihelix transcription factors family is involved in growth and secondary metabolic processes in plants, but no studies have been reported on the involvement of Trihelix genes in secondary metabolic processes in ginseng. In this study, weighted co-expression network analysis, correlation analysis between PgGTs and ginsenosides and key enzyme genes, and interaction network analysis between PgGTs and key enzyme genes were used to screen out the PgGT25-04 gene, which was negatively correlated with ginsenoside synthesis. Using ABA treatment of ginseng hair roots, PgGT genes were found to respond to ABA signals. Analysis of the sequence characteristics and expression pattern of the PgGT25-04 gene in ginseng revealed that its expression is spatiotemporally specific. The interfering vector pBI121-PgGT25-04 containing the PgGT25-04 gene was constructed, and the ginseng adventitious roots were transformed using the Agrobacterium-mediated method to obtain the pBI121-PgGT25-04 positive hairy root monocot line. The saponin contents of positive ginseng hair roots were measured by HPLC, and the changes in PgGT25-04 and key enzyme genes in positive ginseng hair roots were detected via fluorescence quantitative RT-PCR. These results preliminarily identified the role of the PgGT25-04 gene in the secondary metabolism of ginseng in Jilin to provide a theoretical basis for the study of Trihelix transcription factors in Panax ginseng.

Genome-Wide Identification, Evolution, and Expression Analyses of AP2/ERF Family Transcription Factors in Erianthus fulvus

The AP2/ERF transcription factor family is one of the most important gene families in plants and plays a vital role in plant abiotic stress responses. Although Erianthus fulvus is very important in the genetic improvement of sugarcane, there are few studies concerning AP2/ERF genes in E. fulvus. Here, we identified 145 AP2/ERF genes in the E. fulvus genome. Phylogenetic analysis classified them into five subfamilies. Evolutionary analysis showed that tandem and segmental duplication contributed to the expansion of the EfAP2/ERF family. Protein interaction analysis showed that twenty-eight EfAP2/ERF proteins and five other proteins had potential interaction relationships. Multiple cis-acting elements present in the EfAP2/ERF promoter were related to abiotic stress response, suggesting that EfAP2/ERF may contribute to adaptation to environmental changes. Transcriptomic and RT-qPCR analyses revealed that EfDREB10, EfDREB11, EfDREB39, EfDREB42, EfDREB44, EfERF43, and EfAP2-13 responded to cold stress, EfDREB5 and EfDREB42 responded to drought stress, and EfDREB5, EfDREB11, EfDREB39, EfERF43, and EfAP2-13 responded to ABA treatment. These results will be helpful for better understanding the molecular features and biological role of the E. fulvus AP2/ERF genes and lay a foundation for further research on the function of EfAP2/ERF genes and the regulatory mechanism of the abiotic stress response.

Analysis of metabolites and metabolism-mediated biological activity assessment of ginsenosides on microfluidic co-culture system

In vivo, the complex process of drugs metabolism alters the change in drug composition and determines the final pharmacological properties of oral drugs. Ginsenosides are primary constituents of ginseng, whose pharmacological activities are greatly affected by liver metabolism. However, the predictive power of existing in vitro models is poor due to their inability to mimic the complexity of drug metabolism in vivo. The advance of organs-on-chip-based microfluidics system could provide a new in vitro drug screening platform by recapitulating the metabolic process and pharmacological activity of natural product. In this study, an improved microfluidic device was employed to establish an in vitro co-culture model by culturing multiple cell types in compartmentalized microchambers. Different cell lines were seeded on the device to examine the metabolites of ginsenosides from the hepatocytes in top layer and its resulting efficacy on the tumors in bottom layer. Metabolism dependent drug efficacy of Capecitabine in this system demonstrated the model is validated and controllable. High concentrations of CK, Rh2 (S), and Rg3 (S) ginsenosides showed significant inhibitory effects on two types of tumor cells. In addition, apoptosis detection showed that Rg3 (S) through liver metabolism promoted early apoptosis of tumor cells and displayed better anticancer activity than prodrug. The detected ginsenoside metabolites indicated that some protopanaxadiol saponins were converted into other anticancer aglycones in varying degrees due to orderly de-sugar and oxidation. Ginsenosides exhibited different efficacy on target cells by impacting their viabilities, indicating hepatic metabolism plays an important role in determining ginsenosides efficacy. In conclusion, this microfluidic co-culture system is simple, scalable, and possibly widely applicable in evaluating anticancer activity and metabolism of drug during the early developmental phases of natural product.

Genome-wide identification and systematic analysis of the HD-Zip gene family and its roles in response to pH in Panax ginseng Meyer

BACKGROUND

Ginseng, Panax ginseng Meyer, is a traditional herb that is immensely valuable both for human health and medicine and for medicinal plant research. The homeodomain leucine zipper (HD-Zip) gene family is a plant-specific transcription factor gene family indispensable in the regulation of plant growth and development and plant response to environmental stresses.

RESULTS

We identified 117 HD-Zip transcripts from the transcriptome of ginseng cv. Damaya that is widely grown in Jilin, China where approximately 60% of the world's ginseng is produced. These transcripts were positioned to 64 loci in the ginseng genome and the ginseng HD-Zip genes were designated as PgHDZ genes. Identification of 82 and 83 PgHDZ genes from the ginseng acc. IR826 and cv. ChP genomes, respectively, indicated that the PgHDZ gene family consists of approximately 80 PgHDZ genes. Phylogenetic analysis showed that the gene family originated after Angiosperm split from Gymnosperm and before Dicots split from Monocots. The gene family was classified into four subfamilies and has dramatically diverged not only in gene structure and functionality but also in expression characteristics. Nevertheless, co-expression network analysis showed that the activities of the genes in the family remain significantly correlated, suggesting their functional correlation. Five hub PgHDZ genes were identified that might have central functions in ginseng biological processes and four of them were shown to be actively involved in plant response to environmental pH stress in ginseng.

CONCLUSIONS

The PgHDZ gene family was identified from ginseng and analyzed systematically. Five potential hub genes were identified and four of them were shown to be involved in ginseng response to environmental pH stress. The results provide new insights into the characteristics, diversity, evolution, and functionality of the PgHDZ gene family in ginseng and lay a foundation for comprehensive research of the gene family in plants.

Genome-Wide Identification of AP2/ERF Transcription Factor Family and Functional Analysis of DcAP2/ERF#96 Associated with Abiotic Stress in Dendrobium catenatum

APETALA2/Ethylene Responsive Factor (AP2/ERF) family plays important roles in reproductive development, stress responses and hormone responses in plants. However, AP2/ERF family has not been systematically studied in Dendrobium catenatum. In this study, 120 AP2/ERF family members were identified for the first time in D. catenatum, which were divided into four groups (AP2, RAV, ERF and DREB subfamily) according to phylogenetic analysis. Gene structures and conserved motif analysis showed that each DcAP2/ERF family gene contained at least one AP2 domain, and the distribution of motifs varied among subfamilies. Cis-element analysis indicated that DcAP2/ERF genes contained abundant cis-elements related to hormone signaling and stress response. To further identify potential genes involved in drought stress, 12 genes were selected to detect their expression under drought treatment through qRT-PCR analysis and DcAP2/ERF#96, a nuclear localized ethylene-responsive transcription factor, showed a strong response to PEG treatment. Overexpression of DcAP2/ERF#96 in Arabidopsis showed sensitivity to ABA. Molecular, biochemical and genetic assays indicated that DcAP2ERF#96 interacts with DREB2A and directly inhibits the expression of P5CS1 in response to the ABA signal. Taken together, our study provided a molecular basis for the intensive study of DcAP2/ERF genes and revealed the biological function of DcAP2ERF#96 involved in the ABA signal.

Regulatory network of ginsenoside biosynthesis under Ro stress in the hairy roots of Panax ginseng revealed by RNA sequencing

P. ginseng C.A. Meyer is a valuable Chinese herbal medicine that belongs to the Araliaceae family. Major obstacles to the continuous cropping of ginseng have severely restricted the sustainable development of the ginseng industry. The allelopathic effects of triterpenoid saponins play an important role in disorders related to continuous cropping; however, the mechanisms underlying the allelopathic autotoxicity of triterpenoid ginsenosides remain unknown. In this study, we performed mRNA and miRNA sequencing analyses to identify candidate genes and miRNAs that respond differentially to ginsenoside Ro stress in ginseng and their targets. The growth of the ginseng hairy roots was significantly inhibited under Ro stress (0.5 mg/L, Ro-0.5). The inhibition of root growth and injury to root-tip cells promoted the accumulation of the endogenous hormones indole-3-acetic acid and salicylic acid and inhibited the accumulation of abscisic acid and jasmonate acid. The accumulation of ginsenosides, except Rg3, was significantly inhibited under Ro-0.5 stress. An mRNA analysis of the Ro-0.5 and control groups showed that differentially expressed genes were mostly concentrated in the hormone signal transduction pathway. ARF7 and EFM were upregulated, whereas XTH23 and ZOX1 were downregulated. These genes represent important potential candidates for hormone-responsive continuous cropping diseases. In total, 74 differentially expressed miRNAs were identified based on the miRNA sequencing analysis, of which 22 were upregulated and 52 were downregulated. The target genes of ptc-miR156k_L + 1, mtr-miR156b-5p, gma-miR156a_R + 1, and mtr-miR156e all belonged to TRINITY_DN14567_c0_g4, which is a gene in the plant hormone signal transduction pathway. These four miRNAs were all negatively correlated with mRNA, indicating their likely involvement in the response of ginseng to continuous cropping disorders and the regulation of ginsenoside synthesis. Our findings provide useful insights for removing the barriers to continuous ginseng cropping and have important implications in the genetic engineering of plant stress responses.

Functional analysis of mitogen-activated protein kinases (MAPKs) in potato under biotic and abiotic stress

Biotic and abiotic stresses are the main constrain of potato (Solanum tuberosum L.) production all over the world. To overcome these hurdles, many techniques and mechanisms have been used for increasing food demand for increasing population. One of such mechanism is mitogen-activated protein kinase (MAPK) cascade, which is significance regulators of MAPK pathway under various biotic and abiotic stress conditions in plants. However, the acute role in potato for various biotic and abiotic resistance is not fully understood. In eukaryotes including plants, MAPK transfer information from sensors to responses. In potato, biotic and abiotic stresses, as well as a range of developmental responses including differentiation, proliferation, and cell death in plants, MAPK plays an essential role in transduction of diverse extracellular stimuli. Different biotic and abiotic stress stimuli such as pathogen (bacteria, virus, and fungi, etc.) infections, drought, high and low temperatures, high salinity, and high or low osmolarity are induced by several MAPK cascade and MAPK gene families in potato crop. The MAPK cascade is synchronized by numerous mechanisms, including not only transcriptional regulation but also through posttranscriptional regulation such as protein-protein interactions. In this review, we will discuss the recent detailed functional analysis of certain specific MAPK gene families which are involved in resistance to various biotic and abiotic stresses in potato. This study will also provide new insights into functional analysis of various MAPK gene families in biotic and abiotic stress response as well as its possible mechanism.

Genome-wide characterization and identification of candidate ERF genes involved in various abiotic stress responses in sesame (Sesamum indicum L.)

BACKGROUND

The adverse effects of climate change on crop production are constraining breeders to develop high-quality environmentally stable varieties. Hence, efforts are being made to identify key genes that could be targeted for enhancing crop tolerance to environmental stresses. ERF transcription factors play an important role in various abiotic stresses in plants. However, the roles of the ERF family in abiotic stresses tolerance are still largely unknown in sesame, the "queen" of oilseed crops.

RESULTS

In total, 114 sesame ERF genes (SiERFs) were identified and characterized. 96.49% of the SiERFs were distributed unevenly on the 16 linkage groups of the sesame genome. The phylogenetic analysis with the Arabidopsis ERFs (AtERFs) subdivided SiERF subfamily proteins into 11 subgroups (Groups I to X; and VI-L). Genes in the same subgroup exhibited similar structure and conserved motifs. Evolutionary analysis showed that the expansion of ERF genes in sesame was mainly induced by whole-genome duplication events. Moreover, cis-acting elements analysis showed that SiERFs are mostly involved in environmental responses. Gene expression profiles analysis revealed that 59 and 26 SiERFs are highly stimulated under drought and waterlogging stress, respectively. In addition, qRT-PCR analyses indicated that most of SiERFs are also significantly up-regulated under osmotic, submerge, ABA, and ACC stresses. Among them, SiERF23 and SiERF54 were the most induced by both the abiotic stresses, suggesting their potential for targeted improvement of sesame response to multiple abiotic stresses.

CONCLUSION

This study provides a comprehensive understanding of the structure, classification, evolution, and abiotic stresses response of ERF genes in sesame. Moreover, it offers valuable gene resources for functional characterization towards enhancing sesame tolerance to multiple abiotic stresses.

The first high-quality chromosomal genome assembly of a medicinal and edible plant Arctium lappa

Arctium lappa has a long medicinal and edible history with great economic importance. Here, the first high-quality chromosome-level draft genome of A. lappa was presented by the Illumina and PacBio sequencing data. The assembled genome was approximately 1.79 Gb with a N50 contig size of 6.88 Mb. Approximately 1.70 Gb (95.4%) of the contig sequences were anchored onto 18 chromosomes using Hi-C data; the scaffold N50 was improved to be 91.64 Mb. Furthermore, we obtained 1.12 Gb (68.46%) of repetitive sequences and 32,771 protein-coding genes; 616 positively selected candidate genes were identified. Among candidate genes related to lignan biosynthesis, the following were found to be highly correlated with the accumulation of arctiin: 4-coumarate-CoA ligase (4CL), dirigent protein (DIR), and hydroxycinnamoyl transferase (HCT). Additionally, we compared the transcriptomes of A. lappa roots at three different developmental stages and identified 8,943 differentially expressed genes (DEGs) in these tissues. These data can be utilized to identify genes related to A. lappa quality or provide a basis for molecular identification and comparative genomics among related species.

Mitogen-Activated Protein Kinase Expression Profiling Revealed Its Role in Regulating Stress Responses in Potato (Solanum tuberosum)

Mitogen-activated protein kinase (MAPK) cascades are the universal signal transduction networks that regulate cell growth and development, hormone signaling, and other environmental stresses. However, their essential contribution to plant tolerance is very little known in the potato (Solanum tuberosum) plant. The current study carried out a genome-wide study of StMAPK and provided a deep insight using bioinformatics tools. In addition, the relative expression of StMAPKs was also assessed in different plant tissues. The similarity search results identified a total of 22 StMAPK genes in the potato genome. The sequence alignment also showed conserved motif TEY/TDY in most StMAPKs with conserved docking LHDXXEP sites. The phylogenetic analysis divided all 22 StMAPK genes into five groups, i.e., A, B, C, D, and E, showing some common structural motifs. In addition, most of the StMAPKs were found in a cluster form at the terminal of chromosomes. The promoter analysis predicted several stress-responsive Cis-acting regulatory elements in StMAPK genes. Gene duplication under selection pressure also indicated several purifying and positive selections in StMAPK genes. In potato, StMAPK2, StMAPK6, and StMAPK19 showed a high expression in response to heat stress. Under ABA and IAA treatment, the expression of the total 20 StMAPK genes revealed that ABA and IAA played an essential role in this defense process. The expression profiling and real-time qPCR (RT-qPCR) exhibited their high expression in roots and stems compared to leaves. These results deliver primary data for functional analysis and provide reference data for other important crops.

Integrating transcriptome and metabolome analyses of the response to cold stress in pumpkin (Cucurbita maxima)

Cucurbita maxima belong to the genus Cucurbita and are of nutritional and economic importance. Physiological activity, transcriptome, and metabolome analyses of leaf samples from the C. maxima inbreding line IL7 treated at 5 °C and 25 °C were performed. Cold stress resulted in a significant increase in the malondialdehyde content, relative electrical conductivity, soluble protein, sugar content, and catalase activity. A total of 5,553 differentially expressed genes were identified, of which 2,871 were up-regulated and 2,682 down-regulated. In addition, the transcription of differentially expressed genes in the plant hormone signal transduction pathway and transcription factor families of AP2/ERF, bHLH, WRKY, MYB, and HSF was activated. Moreover, 114 differentially expressed metabolites were identified by gas chromatography time-of-flight mass spectrometry, particularly through the analysis of carboxylic acids and derivatives, and organooxygen compounds. The demonstration of a series of potential metabolites and corresponding genes highlighted a comprehensive regulatory mechanism. These findings will provide novel insights into the molecular mechanisms associated with the response to cold stress in C. maxima.

The AP2/ERF Gene Family in Triticum durum: Genome-Wide Identification and Expression Analysis under Drought and Salinity Stresses

Members of the AP2/ERF transcription factor family play critical roles in plant development, biosynthesis of key metabolites, and stress response. A detailed study was performed to identify TtAP2s/ERFs in the durum wheat (Triticum turgidum ssp. durum) genome, which resulted in the identification of 271 genes distributed on chromosomes 1A-7B. By carrying 27 genes, chromosome 6A had the highest number of TtAP2s/ERFs. Furthermore, a duplication assay of TtAP2s/ERFs demonstrated that 70 duplicated gene pairs had undergone purifying selection. According to RNA-seq analysis, the highest expression levels in all tissues and in response to stimuli were associated with DRF and ERF subfamily genes. In addition, the results revealed that TtAP2/ERF genes have tissue-specific expression patterns, and most TtAP2/ERF genes were significantly induced in the root tissue. Additionally, 13 TtAP2/ERF genes (six ERFs, three DREBs, two DRFs, one AP2, and one RAV) were selected for further analysis via qRT-PCR of their potential in coping with drought and salinity stresses. The TtAP2/ERF genes belonging to the DREB subfamily were markedly induced under both drought-stress and salinity-stress conditions. Furthermore, docking simulations revealed several residues in the pocket sites of the proteins associated with the stress response, which may be useful in future site-directed mutagenesis studies to increase the stress tolerance of durum wheat. This study could provide valuable insights for further evolutionary and functional assays of this important gene family in durum wheat.

Identification, characterization and functional differentiation of the NAC gene family and its roles in response to cold stress in ginseng, Panax ginseng C.A. Meyer

The NAC gene family is one of the important plant-specific transcription factor families involved in variety of physiological processes. It has been found in several plant species; however, little is known about the gene family in ginseng, Panax ginseng C.A. Meyer. Here we report identification and systematic analysis of this gene family in ginseng. A total of 89 NAC genes, designated PgNAC01 to PgNAC89, are identified. These genes are alternatively spliced into 251 transcripts at fruiting stage of a four-year-old ginseng plant. The genes of this gene family have five conserved motifs and are clustered into 11 subfamilies, all of which are shared with the genes of the NAC gene families identified in the dicot and monocot model plant species, Arabidopsis and rice. This result indicates that the PgNAC gene family is an ancient and evolutionarily inactive gene family. Gene ontology (GO) analysis shows that the functions of the PgNAC gene family have been substantially differentiated; nevertheless, over 86% the PgNAC transcripts remain functionally correlated. Finally, five of the PgNAC genes, PgNAC05-2, PgNAC41-2, PgNAC48, PgNAC56-1, and PgNAC59, are identified to be involved in plant response to cold stress, suggesting that this gene family plays roles in response to cold stress in ginseng. These results, therefore, provide new insights into functional differentiation and evolution of a gene family in plants and gene resources necessary to comprehensively determine the functions of the PgNAC gene family in response to cold and other abiotic stresses in ginseng.

Genome-wide identification of mitogen-activated protein kinase (MAPK) cascade and expression profiling of CmMAPKs in melon (Cucumis melo L.)

Mitogen-activated protein kinase (MAPK) is a form of serine/threonine protein kinase that activated by extracellular stimulation acting through the MAPK cascade (MAPKKK-MAPKK-MAPK). The MAPK cascade gene family, an important family of protein kinases, plays a vital role in responding to various stresses and hormone signal transduction processes in plants. In this study, we identified 14 CmMAPKs, 6 CmMAPKKs and 64 CmMAPKKKs in melon genome. Based on structural characteristics and a comparison of phylogenetic relationships of MAPK gene families from Arabidopsis, cucumber and watermelon, CmMAPKs and CmMAPKKs were categorized into 4 groups, and CmMAPKKKs were categorized into 3 groups. Furthermore, chromosome location revealed an unevenly distribution on chromosomes of MAPK cascade genes in melon, respectively. Eventually, qRT-PCR analysis showed that all 14 CmMAPKs had different expression patterns under drought, salt, salicylic acid (SA), methyl jasmonate (MeJA), red light (RL), and Podosphaera xanthii (P. xanthii) treatments. Overall, the expression levels of CmMAPK3 and CmMAPK7 under different treatments were higher than those in control. Our study provides an important basis for future functional verification of MAPK genes in regulating responses to stress and signal substance in melon.