Genome-Wide Characterization, Expression Profile Analysis of WRKY Family Genes in Santalum album and Functional Identification of Their Role in Abiotic Stress

Genome-wide characterization and expression analysis of WRKY family genes in the biosynthesis of triptolide in Tripterygium wilfordii

BACKGROUND

WRKY transcription factors play a vital role in regulating plant growth, development, and secondary metabolism. Tripterygium wilfordii is a medicinal plant that has been widely utilized in rheumatoid arthritis therapy; it contains triptolide, a prominent bioactive constituent exhibiting potent anti-inflammatory and anti-tumor properties. However, the mechanism underlying the regulatory effects of WRKY on triptolide biosynthesis is poorly understood.

RESULTS

In this study, 95 TwWRKY genes were identified in the T. wilfordii genome, which were divided into three groups. Phylogenetic analysis indicated that the TwWRKY were conservative relative to other plants. Collinearity analysis revealed that gene duplications played a crucial role in the evolution of this gene family. Transcriptome data from various plant tissues were integrated by correlation analysis, and a gene-to-metabolite network was successfully mapped; consequently, 32 TwWRKY genes were selected as potential regulators of triptolide biosynthesis. Furthermore, the expression changes in the 32 TwWRKY genes were analyzed following methyl jasmonate (MeJA) induction, and the key candidates likely to regulate the biosynthesis of triptolide were screened. Finally, we performed subcellular localization on the key candidate gene TW23G00056.1 and found that it plays its biological role in the nucleus.

CONCLUSION

Our study provides a valuable resource for further research on TwWRKY in T. wilfordii. The candidate genes reported here lay the foundation for elucidating the regulatory mechanism of triptolide.

Genome-Wide Identification and Characterization of the WRKY Gene Family and Their Associated Regulatory Elements in Fortunella hindsii

Background

The WRKY gene family is identified as one of the most prominent transcription factor families in plants and is involved in various biological processes such as metabolism, growth and development, and response to biotic and abiotic stresses. In many plant species, the WRKY gene family was widely studied and analyzed but little to no information for Fortunella hindsii. However, the completion of the whole genome sequencing of Fortunella hindsii allowed us to investigate the genome-wide analysis of WRKY proteins.

Objective

The main objective of this study was to analyze and identify the WRKY gene family in Fortunella hindsii genome.

Methodology

Various bioinformatics approaches have been used to conduct this study.

Results

We constituted 46 members of the Fortunella hindsii WRKY gene family, which were unevenly distributed on all nine chromosomes. The phylogenetic relationship of predicted WRKY proteins of Fortunella hindsii with the WRKY proteins of Arabidopsis showed that 46 FhWRKY genes were divided into three main groups (G1, G2, G3) with five subgroups (2A, 2B, 2C, 2D, and 2E) of G2 group. Domain, conserved motif identification, and gene structure were conducted and the results found that these FhWRKY proteins have conserved identical characteristics within groups and maintain differences between groups. In silico subcellular localization, results showed that FhWRKY genes are located in the nucleus. The cis-regulatory element analysis identified several key CREs that are significantly associated with light, hormone responses, and stress. The gene ontology analysis of these predicted FhWRKY genes showed that these genes are significantly enriched in sequence-specific DNA binding, transcriptional activity, cellular biosynthesis, and metabolic processes.

Conclusion

Therefore, overall, our results provided an excellent foundation for further functional characterization of WRKY genes with an aim of Fortunella hindsii citrus crop improvement.

Comparative transcriptomics reveals potential regulators of climate adaptation in Santalum album L. (Indian Sandalwood)

Santalum album L. (Indian Sandalwood), a valued tree species known for its fragrant heartwood and essential oil is facing increasing threat due to severe anthropogenic pressures compounded by climate change which has resulted in depletion of its adaptive gene pool. The present study investigates the transcriptome-level responses of nine sandalwood genotypes sourced from diverse climatic zones to identify adaptive genes in the species. Comparative transcriptomics predicted 727, 1141 and 479 differentially expressed transcripts (DETs) across wet vs. dry; monsoon vs. dry and wet vs. monsoon conditions, respectively, and majority of DETs were up-regulated in samples sourced from high rainfall areas. Transcripts including heat shock proteins, Zinc finger binding protein, ribosomal proteins, transcription factors and protein kinase were identified as probable regulators of climate adaptation in S. album. The expression changes of eight selected transcripts were further validated by real-time quantitative PCR. Protein-protein interaction analysis revealed key hub transcripts involved in climate response, while alternative splicing events in transcripts such as SURP and G-patch domain-containing protein 1-like protein, G-type lectin S-receptor-like serine/threonine protein kinase B120, Tetraspanin-3 and ARM repeat superfamily protein indicated the probable role of alternate splicing in increasing the transcript diversity during adaptation. This study presents the first insight into the molecular mechanisms of climate adaptation in the species and can form the basis for specific interventions such as selective breeding, genetic manipulation, and habitat management for conservation and long-term survival of sandalwood.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-025-04218-4.

Genome-wide identification and expression profiling of WRKY gene family in grain Amaranth (Amaranthus hypochondriacus L.) under salinity and drought stresses

BACKGROUND

The WRKY gene family plays a significant role in plant growth, development, and responses to biotic and abiotic stresses. However, the role of the WRKY gene family has not been reported in Amaranthus hypochondriacus. This study presents a comprehensive genome-wide analysis of the WRKY gene family in grain amaranth (A. hypochondriacus L.), a resilient crop known for its high nutritional value and adaptability to challenging environments.

RESULTS

In this study, 55 WRKY genes (AhyWRKY1-55) were identified in A. hypochondriacus and distributed unevenly across 16 scaffolds. Of these, 50 contained conserved WRKY domains and were classified into three main groups. Group II was further divided into five subgroups (IIa-IIe) based on phylogenetic analysis, with each clade being well supported by conserved motifs. Additionally, the gene structure analysis revealed variations in exon-intron organization. In contrast, motif analysis showed the presence of conserved domains that were similar within the group but differed between groups, suggesting their functional diversity. Cis-acting elements related to plant growth and development and light, hormones, and stress responses were identified. Synteny analysis revealed that 34 (61.8%) of the genes originated from tandem duplication, indicating the role of tandem duplication in the expansion of the A. hypochondriacus WRKY gene family. Protein-protein interaction analysis suggested that AhyWRKY3, AhyWRKY27, AhyWRKY28, AhyWRKY36, and AhyWRKY52 were hub genes involved in the complex protein interaction network. Using in silico and real-time quantitative PCR, expression analysis revealed tissue- and condition-specific expression patterns of AhyWRKY genes. Notably, under drought stress, AhyWRKY39, AhyWRKY40, AhyWRKY54, and AhyWRKY01 showed increased expression, while under salt stress, AhyWRKY40, AhyWRKY54, AhyWRKY39, AhyWRKY49, and AhyWRKY8 were upregulated at 30 days, suggesting that these genes may play key role in response to salinity stress.

CONCLUSIONS

The present study provides valuable insights into the organization and evolutionary patterns of the WRKY gene family in amaranth. It also identifies putative candidate WRKY genes that may play a role in conferring drought and salt tolerance. Overall, this study lays a foundation for further functional validation of these WRKY candidate genes, facilitating their exploitation in the amaranth genetic improvement programs to develop stress-resilient varieties.

Comprehensive Expression Analysis of the WRKY Gene Family in Phoebe bournei under Drought and Waterlogging Stresses

In response to biotic and abiotic stresses, the WRKY gene family plays a crucial role in plant growth and development. This study focused on Phoebe bournei and involved genome-wide identification of WRKY gene family members, clarification of their molecular evolutionary characteristics, and comprehensive mapping of their expression profiles under diverse abiotic stress conditions. A total of 60 WRKY gene family members were identified, and their phylogenetic classification revealed three distinct groups. A conserved motif analysis underscored the significant conservation of motif 1 and motif 2 among the majority of PbWRKY proteins, with proteins within the same class sharing analogous gene structures. Furthermore, an examination of cis-acting elements and protein interaction networks revealed several genes implicated in abiotic stress responses in P. bournei. Transcriptomic data were utilized to analyze the expression patterns of WRKY family members under drought and waterlogged conditions, with subsequent validation by quantitative real-time PCR (RT-qPCR) experiments. Notably, PbWRKY55 exhibited significant expression modulation under drought stress; PbWRKY36 responded prominently to waterlogging stress; and PbWRKY18, PbWRKY38, and PbWRKY57 demonstrated altered expression under both drought and waterlogging stresses. This study revealed the PbWRKY candidate genes that potentially play a pivotal role in enhancing abiotic stress resilience in P. bournei. The findings have provided valuable insights and knowledge that can guide further research aimed at understanding and addressing the impacts of abiotic stress within this species.

Genome-wide identification and expression analysis of WRKY gene family members in red clover (Trifolium pratense L.)

Trifolium pratense is an important legume forage grass and a key component of sustainable livestock development. Serving as an essential component, the WRKY gene family, a crucial group of regulatory transcription factors in plants, holds significant importance in their response to abiotic stresses. However, there has been no systematic analysis conducted on the WRKY gene family in Trifolium pratense. This study conducted a comprehensive genomic characterization of the WRKY gene family in Trifolium pratense, utilizing the latest genomic data, resulting in the identification of 59 TpWRKY genes. Based on their structural features, phylogenetic characteristics, and conserved motif composition, the WRKY proteins were classified into three groups, with group II further subdivided into five subgroups (II-a, II-b, II-c, II-d, and II-e). The majority of the TpWRKYs in a group share a similar structure and motif composition. Intra-group syntenic analysis revealed eight pairs of duplicate segments. The expression patterns of 59 TpWRKY genes in roots, stems, leaves, and flowers were examined by analyzing RNA-seq data. The expression of 12 TpWRKY genes under drought, low-temperature (4°C), methyl jasmonate (MeJA) and abscisic acid (ABA) stresses was analyzed by RT-qPCR. The findings indicated that TpWRKY46 was highly induced by drought stress, and TpWRKY26 and TpWRKY41 were significantly induced by low temperature stress. In addition, TpWRKY29 and TpWRKY36 were greatly induced by MeJA stress treatment, and TpWRKY17 was significantly upregulated by ABA stress treatment. In this research, we identified and comprehensively analyzed the structural features of the WRKY gene family in T.pratense, along with determined the possible roles of WRKY candidate genes in abiotic stress. These discoveries deepen our understandings of how WRKY transcription factors contribute to species evolution and functional divergence, laying a solid molecular foundation for future exploration and study of stress resistance mechanisms in T.pratense.

Chromosome-level genome assemblies from two sandalwood species provide insights into the evolution of the Santalales

Sandalwood is one of the most expensive woods in the world and is well known for its long-lasting and distinctive aroma. In our study, chromosome-level genome assemblies for two sandalwood species (Santalum album and Santalum yasi) were constructed by integrating NGS short reads, RNA-seq, and Hi-C libraries with PacBio HiFi long reads. The S. album and S. yasi genomes were both assembled into 10 pseudochromosomes with a length of 229.59 Mb and 232.64 Mb, containing 21,673 and 22,816 predicted genes and a repeat content of 28.93% and 29.54% of the total genomes, respectively. Further analyses resolved a Santalum-specific whole-genome triplication event after divergence from ancestors of the Santalales lineage Malania, yet due to dramatic differences in transposon content, the Santalum genomes were only one-sixth the size of the Malania oleifera genome. Examination of RNA-seq data revealed a suite of genes that are differentially expressed in haustoria and might be involved in host hemiparasite interactions. The two genomes presented here not only provide an important comparative dataset for studying genome evolution in early diverging eudicots and hemiparasitic plants but will also hasten the application of conservation genomics for a lineage of trees recovering from decades of overexploitation.

Progress in the understanding of WRKY transcription factors in woody plants

The WRKY transcription factor (TF) family, named for its iconic WRKY domain, is among the largest and most functionally diverse TF families in higher plants. WRKY TFs typically interact with the W-box of the target gene promoter to activate or inhibit the expression of downstream genes; these TFs are involved in the regulation of various physiological responses. Analyses of WRKY TFs in numerous woody plant species have revealed that WRKY family members are broadly involved in plant growth and development, as well as responses to biotic and abiotic stresses. Here, we review the origin, distribution, structure, and classification of WRKY TFs, along with their mechanisms of action, the regulatory networks in which they are involved, and their biological functions in woody plants. We consider methods currently used to investigate WRKY TFs in woody plants, discuss outstanding problems, and propose several new research directions. Our objective is to understand the current progress in this field and provide new perspectives to accelerate the pace of research that enable greater exploration of the biological functions of WRKY TFs.

Genome-wide characterization of R2R3-MYB gene family in Santalum album and their expression analysis under cold stress

Sandalwood (Santalum album) is a high-value multifunctional tree species that is rich in aromatic substances and is used in medicine and global cosmetics. Due to the scarcity of land resources in tropical and subtropical regions, land in temperate regions is a potential resource for the development of S. album plantations in order to meet the needs of S. album production and medicine. The R2R3-MYB transcription factor family is one of the largest in plants and plays an important role in the response to various abiotic stresses. However, the R2R3-MYB gene family of S. album has not been studied. In this study, 144 R2R3-MYB genes were successfully identified in the assembly genome sequence, and their characteristics and expression patterns were investigated under various durations of low temperature stress. According to the findings, 31 of the 114 R2R3-MYB genes showed significant differences in expression after cold treatment. Combining transcriptome and weighted gene co-expression network analysis (WGCNA) revealed three key candidate genes (SaMYB098, SaMYB015, and SaMYB068) to be significantly involved in the regulation of cold resistance in S. album. The structural characteristics, evolution, and expression pattern of the R2R3-MYB gene in S. album were systematically examined at the whole genome level for the first time in this study. It will provide important information for future research into the function of the R2R3-MYB genes and the mechanism of cold stress response in S. album.

WRKY transcription factors: evolution, regulation, and functional diversity in plants

The recent advancements in sequencing technologies and informatic tools promoted a paradigm shift to decipher the hidden biological mysteries and transformed the biological issues into digital data to express both qualitative and quantitative forms. The transcriptomic approach, in particular, has added new dimensions to the versatile essence of plant genomics through the large and deep transcripts generated in the process. This has enabled the mining of super families from the sequenced plants, both model and non-model, understanding their ancestry, diversity, and evolution. The elucidation of the crystal structure of the WRKY proteins and recent advancement in computational prediction through homology modeling and molecular dynamic simulation has provided an insight into the DNA-protein complex formation, stability, and interaction, thereby giving a new dimension in understanding the WRKY regulation. The present review summarizes the functional aspects of the high volume of sequence data of WRKY transcription factors studied from different species, till date. The review focuses on the dynamics of structural classification and lineage in light of the recent information. Additionally, a comparative analysis approach was incorporated to understand the functions of the identified WRKY transcription factors subjected to abiotic (heat, cold, salinity, senescence, dark, wounding, UV, and carbon starvation) stresses as revealed through various sets of studies on different plant species. The review will be instrumental in understanding the events of evolution and the importance of WRKY TFs under the threat of climate change, considering the new scientific evidences to propose a fresh perspective.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Akebiatrifoliata: A Bioinformatics Study

WRKY transcription factors have been found in most plants and play an important role in regulating organ growth and disease response. Outlining the profile of WRKY genes is a very useful project for studying morphogenesis and resistance formation. In the present study, a total of 63 WRKY genes consisting of 13 class I, 41 class II, and 9 class III genes were identified from the newly published A. trifoliata genome, of which 62 were physically distributed on all 16 chromosomes. Structurally, two AkWRKY genes (AkWRKY6 and AkWRKY52) contained four domains, and AkWRKY17 lacked the typical heptapeptide structure. Evolutionarily, 42, 16, and 5 AkWRKY genes experienced whole genome duplication (WGD) or fragmentation, dispersed duplication, and tandem duplication, respectively; 28 Ka/Ks values of 30 pairs of homologous genes were far lower than 1, while those of orthologous gene pairs between AkWRKY41 and AkWRKY52 reached up to 2.07. Transcriptome analysis showed that many of the genes were generally expressed at a low level in 12 fruit samples consisting of three tissues, including rind, flesh, and seeds, at four developmental stages, and interaction analysis between AkWRKY and AkNBS genes containing W-boxes suggested that AkWRKY24 could play a role in plant disease resistance by positively regulating AkNBS18. In summary, the WRKY gene family of A. trifoliata was systemically characterized for the first time, and the data and information obtained regarding AkWRKY could be very useful in further theoretically elucidating the molecular mechanisms of plant development and response to pathogens and practically improving favorable traits such as disease resistance.

Identification and Functional Analysis of SabHLHs in Santalum album L

Santalum album L., a semi-parasitic evergreen tree, contains economically important essential oil, rich in sesquiterpenoids, such as (Z) α- and (Z) β-santalol. However, their transcriptional regulations are not clear. Several studies of other plants have shown that basic-helix-loop-helix (bHLH) transcription factors (TFs) were involved in participating in the biosynthesis of sesquiterpene synthase genes. Herein, bHLH TF genes with similar expression patterns and high expression levels were screened by co-expression analysis, and their full-length ORFs were obtained. These bHLH TFs were named SaMYC1, SaMYC3, SaMYC4, SaMYC5, SabHLH1, SabHLH2, SabHLH3, and SabHLH4. All eight TFs had highly conserved bHLH domains and SaMYC1, SaMYC3, SaMYC4, and SaMYC5, also had highly conserved MYC domains. It was indicated that the eight genes belonged to six subfamilies of the bHLH TF family. Among them, SaMYC1 was found in both the nucleus and the cytoplasm, while SaMYC4 was only localized in the cytoplasm and the remaining six TFs were localized in nucleus. In a yeast one-hybrid experiment, we constructed decoy vectors pAbAi-SSy1G-box, pAbAi-CYP2G-box, pAbAi-CYP3G-box, and pAbAi-CYP4G-box, which had been transformed into yeast. We also constructed pGADT7-SaMYC1 and pGADT7-SabHLH1 capture vectors and transformed them into bait strains. Our results showed that SaMYC1 could bind to the G-box of SaSSy, and the SaCYP736A167 promoter, which SaSSy proved has acted as a key enzyme in the synthesis of santalol sesquiterpenes and SaCYP450 catalyzed the ligation of santalol sesquiterpenes into terpene. We have also constructed pGreenII 62-SK-SaMYC1, pGreenII 0800-LUC-SaSSy and pGreenII 0800-LUC-SaCYP736A167 via dual-luciferase fusion expression vectors and transformed them into Nicotiana benthamiana using an Agrobacterium-mediated method. The results showed that SaMYC1 was successfully combined with SaSSy or SaCYP736A167 promoter and the LUC/REN value was 1.85- or 1.55-fold higher, respectively, than that of the control group. Therefore, we inferred that SaMYC1 could activate both SaSSy and SaCYP736A167 promoters.

Characterization of the WRKY gene family in Akebia trifoliata and their response to Colletotrichum acutatum

BACKGROUND

Akebia trifoliata, belonging to the Lardizabalaceae family, is a well-known Chinese traditional medicinal plant, susceptible to many diseases, such as anthracnose and powdery mildew. WRKY is one of the largest plant-specific transcription factor families and plays important roles in plant growth, development and stress response, especially in disease resistance. However, little was known about the numbers, characters, evolutionary relationship and expression of WRKY genes in A. trifoliata in response to plant disease due to lacking of A. trifoliata genome.

RESULTS

A total of 42 putative AktWRKY genes were identified based on the full-length transcriptome-sequencing data of A. trifoliata. Then 42 AktWRKY genes were divided into three major groups (Group I-III) based on the WRKY domains. Motif analysis showed members within same group shared a similar motif composition, implying a functional conservation. Tissue-specific expression analysis showed that AktWRKY genes could be detected in all tissues, while few AktWRKY genes were tissue specific. We further evaluated the expression of AktWRKY genes in three varieties in response to Colletotrichum acutatum by qRT-PCR. The expression patterns of AktWRKY genes were similar between C01 and susceptible variety I02, but distinctly different in resistant variety H05. In addition, it showed that more than 64 percentages of AktWRKY genes were differentially expressed during fungal infection in I02 and H05. Furthermore, Gene ontology (GO) analysis showed that AktWRKY genes were categorized into 26 functional groups under cellular components, molecular functions and biological processes, and a predicted protein interaction network was also constructed.

CONCLUSIONS

Results of bioinformation analysis and expression patterns implied that AktWRKYs might play multiple function in response to biotic stresses. Our study could facilitate to further investigate the function and regulatory mechanism of the WRKY in A. trifoliata during pathogen response.

Genome-wide identification of the Liriodendron chinense WRKY gene family and its diverse roles in response to multiple abiotic stress

BACKGROUND

Liriodendron chinense (Lchi) is a tree species within the Magnoliaceae family and is considered a basal angiosperm. The too low or high temperature or soil drought will restrict its growth as the adverse environmental conditions, thus improving L. chinense abiotic tolerance was the key issues to study. WRKYs are a major family of plant transcription factors known to often be involved in biotic and abiotic stress responses. So far, it is still largely unknown if and how the LchiWRKY gene family is tied to regulating L. chinense stress responses. Therefore, studying the involvement of the WRKY gene family in abiotic stress regulation in L. chinense could be very informative in showing how this tree deals with such stressful conditions.

RESULTS

In this research, we performed a genome-wide analysis of the Liriodendron chinense (Lchi) WRKY gene family, studying their classification relationships, gene structure, chromosomal locations, gene duplication, cis-element, and response to abiotic stress. The 44 members of the LchiWRKY gene family contain a significant amount of sequence diversity, with their lengths ranging from 525 bp to 40,981 bp. Using classification analysis, we divided the 44 LchiWRKY genes into three phylogenetic groups (I, II, II), with group II then being further divided into five subgroups (IIa, IIb, IIc, IId, IIe). Comparative phylogenetic analysis including the WRKY families from 17 plant species suggested that LchiWRKYs are closely related to the Magnolia Cinnamomum kanehirae WRKY family, and has fewer family members than higher plants. We found the LchiWRKYs to be evenly distributed across 15 chromosomes, with their duplication events suggesting that tandem duplication may have played a major role in LchiWRKY gene expansion model. A Ka/Ks analysis indicated that they mainly underwent purifying selection and distributed in the group IId. Motif analysis showed that LchiWRKYs contained 20 motifs, and different phylogenetic groups contained conserved motif. Gene ontology (GO) analysis showed that LchiWRKYs were mainly enriched in two categories, i.e., biological process and molecular function. Two group IIc members (LchiWRKY10 and LchiWRKY37) contain unique WRKY element sequence variants (WRKYGKK and WRKYGKS). Gene structure analysis showed that most LchiWRKYs possess 3 exons and two different types of introns: the R- and V-type which are both contained within the WRKY domain (WD). Additional promoter cis-element analysis indicated that 12 cis-elements that play different functions in environmental adaptability occur across all LchiWRKY groups. Heat, cold, and drought stress mainly induced the expression of group II and I LchiWRKYs, some of which had undergone gene duplication during evolution, and more than half of which had three exons. LchiWRKY33 mainly responded to cold stress and LchiWRKY25 mainly responded to heat stress, and LchiWRKY18 mainly responded to drought stress, which was almost 4-fold highly expressed, while 5 LchiWRKYs (LchiWRKY5, LchiWRKY23, LchiWRKY14, LchiWRKY27, and LchiWRKY36) responded equally three stresses with more than 6-fold expression. Subcellular localization analysis showed that all LchiWRKYs were localized in the nucleus, and subcellular localization experiments of LchiWRKY18 and 36 also showed that these two transcription factors were expressed in the nucleus.

CONCLUSIONS

This study shows that in Liriodendron chinense, several WRKY genes like LchiWRKY33, LchiWRKY25, and LchiWRKY18, respond to cold or heat or drought stress, suggesting that they may indeed play a role in regulating the tree's response to such conditions. This information will prove a pivotal role in directing further studies on the function of the LchiWRKY gene family in abiotic stress response and provides a theoretical basis for popularizing afforestation in different regions of China.

Identification of WRKY transcription factors involved in regulating the biosynthesis of the anti-cancer drug camptothecin in Ophiorrhiza pumila

Camptothecin is a chemotherapeutic drug widely used to treat various cancers. Ophiorrhiza pumila is an ideal plant model for the study of camptothecin production, with various advantages for studying camptothecin biosynthesis and regulation. The DNA-binding WRKY transcription factors have a key regulatory role in secondary metabolite biosynthesis in plants. However, little is currently known about their involvement in camptothecin biosynthesis in O. pumila. We identified 46 OpWRKY genes unevenly distributed on the 11 chromosomes of O. pumila. Phylogenetic and multiple sequence alignment analyses divided the OpWRKY proteins into three subfamilies. Based on spatial expression and co-expression, we targeted the candidate gene OpWRKY6. Overexpression of OpWRKY6 significantly reduced the accumulation of camptothecin compared with the control. Conversely, camptothecin accumulation increased in OpWRKY6 knockout lines. Further biochemical assays showed that OpWRKY6 negatively regulates camptothecin biosynthesis from both the iridoid and shikimate pathways by directly downregulating the gene expression of OpGES, Op10HGO, Op7DLH, and OpTDC. Our data provide direct evidence for the involvement of WRKYs in the regulation of camptothecin biosynthesis and offer valuable information for enriching the production of camptothecin in plant systems.

Genome-Wide Identification and Analysis of NAC Transcription Factor Family in Two Diploid Wild Relatives of Cultivated Sweet Potato Uncovers Potential NAC Genes Related to Drought Tolerance

NAC (NAM, ATAF1/2, and CUC2) proteins play a pivotal role in modulating plant development and offer protection against biotic and abiotic stresses. Until now, no systematic knowledge of NAC family genes is available for the food security crop, sweet potato. Here, a comprehensive genome-wide survey of NAC domain-containing proteins identified 130 ItbNAC and 144 ItfNAC genes with full length sequences in the genomes of two diploid wild relatives of cultivated sweet potato, Ipomoea triloba and Ipomoea trifida, respectively. These genes were physically mapped onto 15 I. triloba and 16 I. trifida chromosomes, respectively. Phylogenetic analysis divided all 274 NAC proteins into 20 subgroups together with NAC transcription factors (TFs) from Arabidopsis. There were 9 and 15 tandem duplication events in the I. triloba and I. trifida genomes, respectively, indicating an important role of tandem duplication in sweet potato gene expansion and evolution. Moreover, synteny analysis suggested that most NAC genes in the two diploid sweet potato species had a similar origin and evolutionary process. Gene expression patterns based on RNA-Seq data in different tissues and in response to various hormone, biotic or abiotic treatments revealed their possible involvement in organ development and response to various biotic/abiotic stresses. The expression of 36 NAC TFs, which were upregulated in the five tissues and in response to mannitol treatment, was also determined by real-time quantitative polymerase chain reaction (RT-qPCR) in hexaploid cultivated sweet potato exposed to drought stress. Those results largely corroborated the expression profile of mannitol treatment uncovered by the RNA-Seq data. Some significantly up-regulated genes related to drought stress, such as ItbNAC110, ItbNAC114, ItfNAC15, ItfNAC28, and especially ItfNAC62, which had a conservative spatial conformation with a closely related paralogous gene, ANAC019, may be potential candidate genes for a sweet potato drought tolerance breeding program. This analysis provides comprehensive and systematic information about NAC family genes in two diploid wild relatives of cultivated sweet potato, and will provide a blueprint for their functional characterization and exploitation to improve the tolerance of sweet potato to abiotic stresses.

Genome-Wide Analysis of WRKY Gene Family and the Dynamic Responses of Key WRKY Genes Involved in Ostrinia furnacalis Attack in Zea mays

WRKY transcription factors comprise one of the largest gene families and serve as key regulators of plant defenses against herbivore attack. However, studies related to the roles of WRKY genes in response to herbivory are limited in maize. In this study, a total of 128 putative maize WRKY genes (ZmWRKYs) were identified from the new maize genome (v4). These genes were divided into seven subgroups (groups I, IIa–e, and III) based on phylogenomic analysis, with distinct motif compositions in each subgroup. Syntenic analysis revealed that 72 (56.3%) of the genes were derived from either segmental or tandem duplication events (69 and 3, respectively), suggesting a pivotal role of segmental duplication in the expansion of the ZmWRKY family. Importantly, transcriptional regulation prediction showed that six key WRKY genes contribute to four major defense-related pathways: L-phenylalanine biosynthesis II and flavonoid, benzoxazinoid, and jasmonic acid (JA) biosynthesis. These key WRKY genes were strongly induced in commercial maize (Jingke968) infested with the Asian corn borer, Ostrinia furnacalis, for 0, 2, 4, 12 and 24 h in the field, and their expression levels were highly correlated with predicted target genes, suggesting that these genes have important functions in the response to O. furnacalis. Our results provide a comprehensive understanding of the WRKY gene family based on the new assembly of the maize genome and lay the foundation for further studies into functional characteristics of ZmWRKY genes in commercial maize defenses against O. furnacalis in the field.

Genome-Wide Identification and Analysis of the WRKY Gene Family and Cold Stress Response in Acer truncatum

WRKY transcription factors constitute one of the largest gene families in plants and are involved in many biological processes, including growth and development, physiological metabolism, and the stress response. In earlier studies, the WRKY gene family of proteins has been extensively studied and analyzed in many plant species. However, information on WRKY transcription factors in Acer truncatum has not been reported. In this study, we conducted genome-wide identification and analysis of the WRKY gene family in A. truncatum, 54 WRKY genes were unevenly located on all 13 chromosomes of A. truncatum, the highest number was found in chromosomes 5. Phylogenetic relationships, gene structure, and conserved motif identification were constructed, and the results affirmed 54 AtruWRKY genes were divided into nine subgroup groups. Tissue species analysis of AtruWRKY genes revealed which were differently exhibited upregulation in flower, leaf, root, seed and stem, and the upregulation number were 23, 14, 34, 18, and 8, respectively. In addition, the WRKY genes expression in leaf under cold stress showed that more genes were significantly expressed under 0, 6 and 12 h cold stress. The results of this study provide a new insight the regulatory function of WRKY genes under abiotic and biotic stresses.

Genome-Wide Identification and Analysis of NAC Transcription Factor Family in Two Diploid Wild Relatives of Cultivated Sweet Potato Uncovers Potential NAC Genes Related to Drought Tolerance

NAC (NAM, ATAF1/2, and CUC2) proteins play a pivotal role in modulating plant development and offer protection against biotic and abiotic stresses. Until now, no systematic knowledge of NAC family genes is available for the food security crop, sweet potato. Here, a comprehensive genome-wide survey of NAC domain-containing proteins identified 130 ItbNAC and 144 ItfNAC genes with full length sequences in the genomes of two diploid wild relatives of cultivated sweet potato, Ipomoea triloba and Ipomoea trifida, respectively. These genes were physically mapped onto 15 I. triloba and 16 I. trifida chromosomes, respectively. Phylogenetic analysis divided all 274 NAC proteins into 20 subgroups together with NAC transcription factors (TFs) from Arabidopsis. There were 9 and 15 tandem duplication events in the I. triloba and I. trifida genomes, respectively, indicating an important role of tandem duplication in sweet potato gene expansion and evolution. Moreover, synteny analysis suggested that most NAC genes in the two diploid sweet potato species had a similar origin and evolutionary process. Gene expression patterns based on RNA-Seq data in different tissues and in response to various hormone, biotic or abiotic treatments revealed their possible involvement in organ development and response to various biotic/abiotic stresses. The expression of 36 NAC TFs, which were upregulated in the five tissues and in response to mannitol treatment, was also determined by real-time quantitative polymerase chain reaction (RT-qPCR) in hexaploid cultivated sweet potato exposed to drought stress. Those results largely corroborated the expression profile of mannitol treatment uncovered by the RNA-Seq data. Some significantly up-regulated genes related to drought stress, such as ItbNAC110, ItbNAC114, ItfNAC15, ItfNAC28, and especially ItfNAC62, which had a conservative spatial conformation with a closely related paralogous gene, ANAC019, may be potential candidate genes for a sweet potato drought tolerance breeding program. This analysis provides comprehensive and systematic information about NAC family genes in two diploid wild relatives of cultivated sweet potato, and will provide a blueprint for their functional characterization and exploitation to improve the tolerance of sweet potato to abiotic stresses.

Screening and identification of salt-tolerance genes in Sophora alopecuroides and functional verification of SaAQP

Overexpression of SaAQP can improve the salt tolerance of transgenic soybean hairy roots and A. thaliana. Salt stress severely affects crop yield and food security. There is a need to improve the salt tolerance of crops, but the discovery and utilization of salt-tolerance genes remains limited. Owing to its strong stress tolerance, Sophora alopecuroides is ideal for the identification of salt-tolerance genes. Therefore, we aimed to screen and identify the salt-tolerance genes in S. alopecuroides. With a yeast expression library of seedlings, salt-tolerant genes were screened using a salt-containing medium to simulate salt stress. By combining salt-treatment screening and transcriptome sequencing, 11 candidate genes related to salt tolerance were identified, including genes for peroxidase, inositol methyltransferase, aquaporin, cysteine synthase, pectinesterase, and WRKY. The expression dynamics of candidate genes were analyzed after salt treatment of S. alopecuroides, and salt tolerance was verified in yeast BY4743. The candidate genes participated in the salt-stress response in S. alopecuroides, and their overexpression significantly improved the salt tolerance of yeast. Salt tolerance mediated by SaAQP was further verified in soybean hairy roots and Arabidopsis thaliana, and it was found that SaAQP might enhance the salt tolerance of A. thaliana by participating in a reactive oxygen species scavenging mechanism. This result provides new genetic resources in plant breeding for salt resistance.

Genome-wide analysis of BpDof genes and the tolerance to drought stress in birch (Betula platyphylla)

Background

DNA binding with one finger (Dof) proteins are plant-specific transcription factors playing vital roles in developmental processes and stress responses in plants. Nevertheless, the characterizations, expression patterns, and functions of the Dof family under drought stress (a key determinant of plant physiology and metabolic homeostasis) in woody plants remain unclear.

Methods

The birch (Betula platyphylla var. mandshuric) genome and plant TFDB database were used to identify Dof gene family members in birch plants. ClustalW2 of BioEdit v7.2.1, MEGA v7.0, ExPASy ProtParam tool, Subloc, TMHMM v2.0, GSDS v2.0, MEME, TBtools, KaKs Calculator v2.0, and PlantCARE were respectively used to align the BpDof sequences, build a phylogenetic tree, identify the physicochemical properties, analyze the chromosomal distribution and synteny, and identify the cis-elements in the promoter regions of the 26 BpDof genes. Additionally, the birch seedlings were exposed to PEG6000-simulated drought stress, and the expression patterns of the BpDof genes in different tissues were analyzed by qRT-PCR. The histochemical staining and the evaluation of physiological indexes were performed to assess the plant tolerance to drought with transient overexpression of BpDof4, BpDof11, and BpDof17 genes. SPSS software and ANOVA were used to conduct all statistical analyses and determine statistically significant differences between results.

Results

A total of 26 BpDof genes were identified in birch via whole-genome analysis. The conserved Dof domain with a C(x)2C(x)21C(x)2C zinc finger motif was present in all BpDof proteins. These birch BpDofs were classified into four groups (A to D) according to the phylogenetic analysis of Arabidopsis thaliana Dof genes. BpDof proteins within the same group mostly possessed similar motifs, as detected by conserved motif analysis. The exon–intron analysis revealed that the structures of BpDof genes differed, indicating probable gene gain and lose during the BpDof evolution. The chromosomal distribution and synteny analysis showed that the 26 BpDofs were unevenly distributed on 14 chromosomes, and seven duplication events among six chromosomes were found. Cis-acting elements were abundant in the promoter regions of the 26 BpDof genes. qRT-PCR revealed that the expression of the 26 BpDof genes was differentially regulated by drought stress among roots, stems, and leaves. Most BpDof genes responded to drought stress, and BpDof4, BpDof11, and BpDof17­ were significantly up-regulated. Therefore, plants overexpressing these three genes were generated to investigate drought stress tolerance. The BpDof4-, BpDof11-, and BpDof17­-overexpressing plants showed promoted reactive oxygen species (ROS) scavenging capabilities and less severe cell damage, suggesting that they conferred enhanced drought tolerance in birch. This study provided an in-depth insight into the structure, evolution, expression, and function of the Dof gene family in plants.

The whole-genome and expression profile analysis of WRKY and RGAs in Dactylis glomerata showed that DG6C02319.1 and DgWRKYs may cooperate in the immunity against rust

Orchardgrass (Dactylis glomerata) is one of the top four perennial forages worldwide and, despite its large economic advantages, often threatened by various environmental stresses. WRKY transcription factors (TFs) can regulate a variety of plant processes, widely participate in plant responses to biotic and abiotic stresses, and are one of the largest gene families in plants. WRKYs can usually bind W-box elements specifically. In this study, we identified a total of 93 DgWRKY genes and 281 RGAs, including 65, 169 and 47 nucleotide-binding site-leucine-rich repeats (NBS-LRRs), leucine-rich repeats receptor-like protein kinases (LRR-RLKs), and leucine-rich repeats receptor-like proteins (LRR-RLPs), respectively. Through analyzing the expression of DgWRKY genes in orchardgrass under different environmental stresses, it was found that many DgWRKY genes were differentially expressed under heat, drought, submergence, and rust stress. In particular, it was found that the greatest number of genes were differentially expressed under rust infection. Consistently, GO and KEGG enrichment analysis of all genes showed that 78 DgWRKY TFs were identified in the plant–pathogen interaction pathway, with 59 of them differentially expressed. Through cis-acting element prediction, 154 RGAs were found to contain W-box elements. Among them, DG6C02319.1 (a member of the LRR-RLK family) was identified as likely to interact with 14 DGWRKYs. Moreover, their expression levels in susceptible plants after rust inoculation were first up-regulated and then down-regulated, while those in the resistant plants were always up-regulated. In general, DgWRKYs responded to both biotic stress and abiotic stress. DgWRKYs and RGAs may synergistically respond to the response of orchardgrass to rust. This study provides meaningful insight into the molecular mechanisms of WRKY proteins in orchardgrass.

WRKY Transcription Factors Actively Respond to Fusarium oxysporum in Lilium regale

The WRKY transcription factors form a plant-specific superfamily important for regulating plant development, stress responses, and hormone signal transduction. In this study, many WRKY genes (LrWRKY1-35) were identified in Lilium regale, which is a wild lily species highly resistant to Fusarium wilt. These WRKY genes were divided into three classes (I to III) based on a phylogenetic analysis. The Class-II WRKY transcription factors were further divided into five subclasses (IIa, IIb, IIc, IId, and IIe). Moreover, the gene expression patterns based on a quantitative real-time PCR analysis revealed the WRKY genes were differentially expressed in the L. regale roots, stems, leaves, and flowers. Additionally, the expression of the WRKY genes was affected by an infection by Fusarium oxysporum as well as by salicylic acid, methyl jasmonate, ethephon, and hydrogen peroxide treatments. Moreover, the LrWRKY1 protein was localized to the nucleus of onion epidermal cells. The recombinant LrWRKY1 protein purified from Escherichia coli bound specifically to DNA fragments containing the W-box sequence, and a yeast one-hybrid assay indicated that LrWRKY1 can activate transcription. A co-expression assay in tobacco (Nicotiana tabacum) confirmed LrWRKY1 regulates the expression of LrPR10-5. Furthermore, the overexpression of LrWRKY1 in tobacco and the Oriental hybrid 'Siberia' (susceptible to F. oxysporum) increased the resistance of the transgenic plants to F. oxysporum. Overall, LrWRKY1 regulates the expression of the resistance gene LrPR10-5 and is involved in the defense response of L. regale to F. oxysporum. This study provides valuable information regarding the expression and functional characteristics of L. regale WRKY genes.

Transcription Factor CsWRKY65 Participates in the Establishment of Disease Resistance of Citrus Fruits to Penicillium digitatum

Penicillium digitatum is the primary pathogen that causes serious yield losses worldwide. In our previous study, CsWRKY transcription factors (TFs) and some genes associated with immunity were identified in citrus fruits after P. digitatum infection, but little information is available in the literature on the mechanisms of TFs in citrus disease resistance. In this study, the possible mechanisms of CsWRKY65 participating in the establishment of disease resistance were investigated. Results show that CsWRKY65 was a transcriptional activator in the nucleus. The dual-luciferase transient assays and electrophoretic mobility shift assays showed that CsWRKY65 bound with CsRbohB, CsRbohD, CsCDPK33, and CsPR10 promoters to activate gene transcription. Besides, the transient overexpression of CsWRKY65 induced reactive oxygen species accumulation and increased PR gene expression in Nicotiana benthamiana leaves. The transient overexpression of CsWRKY65 in the citrus peel enhanced the disease resistance against P. digitatum. In conclusion, CsWRKY65 is likely to be involved in regulating the disease resistance to P. digitatum of citrus fruits by directly activating the expressions of CsRbohB, CsRbohD, CsCDPK33, and CsPR10. This study provides new information for the mechanism of citrus WRKY TFs participating in the establishment of disease resistance.

Identification, evolution and expression analysis of WRKY gene family in Eucommia ulmoides

The WRKY transcription factors is one of the largest families of transcription factors (TFs) in plants and involved in multiple biological processes. However, the role of the WRKY family had not been reported in Eucommia ulmoides. In this study, 45 WRKY genes (EuWRKY1-45) with conserved WRKY domain were identified in E. ulmoides and classified into three groups. The group II was further divided into five subgroups based on phylogenetic analysis, and each clade was well supported by the conserved motifs. All the genes were located on 34 different scaffolds respectively. A number of development-, light-, hormone-, and stress-related elements were randomly distributed in the promoter sequences of EuWRKYs. Expression profiles indicated that EuWRKY genes were involved in leaf development, and majority of EuWRKYs genes were highly expressed in leaf buds. Co-expression analysis of WRKYs suggested an intricate interplay of growth-related responses. EuWRKY4 was involved in a complex proteins interaction network. Collectively, our results provide extensive insights into the WRKY gene family, thereby contributing to the screening of additional candidate genes in E. ulmoides.

Genomic Characterization of WRKY Transcription Factors Related to Andrographolide Biosynthesis in Andrographis paniculata

Andrographolide, which is enriched in the leaves of Andrographis paniculata, has been known as “natural antibiotic” due to its pharmacological activities such as anti-inflammatory, antimicrobial and antioxidant effects. Several key enzymes in andrographolide biosynthetic pathway have been studied since the genome sequences were released, but its regulatory mechanism remains unknown. WRKY transcription factors proteins have been reported to regulate plant secondary metabolism, development as well as biotic and abiotic stresses. Here, WRKY transcription factors related to andrographolide biosynthesis were systematically identified, including sequences alignment, phylogenetic analysis, chromosomal distribution, gene structure, conserved motifs, synteny, alternative splicing event and Gene ontology (GO) annotation. A total of 58 WRKYs were identified in Chuanxinlian genome and phylogenetically classified into three groups. Moreover, nine WRKY genes underwent alternative splicing events. Furthermore, the combination of binding site prediction, gene-specific expression patterns, and phylogenetic analysis suggested that 7 WRKYs (ApWRKY01, ApWRKY08, ApWRKY12, ApWRKY14, ApWRKY19, ApWRKY20, and ApWRKY50) might regulate andrographolide biosynthesis. This study laid a foundation for understanding the regulatory mechanism of andrographolide biosynthesis and the improvement and breeding of Andrographis paniculata varieties.

Lilium regale Wilson WRKY2 Regulates Chitinase Gene Expression During the Response to the Root Rot Pathogen Fusarium oxysporum

Root rot, mainly caused by Fusarium oxysporum, is the most destructive disease affecting lily (Lilium spp.) production. The WRKY transcription factors (TFs) have important roles during plant immune responses. To clarify the effects of WRKY TFs on plant defense responses to pathogens, a WRKY gene (LrWRKY2) was isolated from Lilium regale Wilson, which is a wild lily species highly resistant to F. oxysporum. The expression of LrWRKY2, which encodes a nuclear protein, is induced by various hormones (methyl jasmonate, ethephon, salicylic acid, and hydrogen peroxide) and by F. oxysporum infection. In this study, LrWRKY2-overexpressing transgenic tobacco plants were more resistant to F. oxysporum than the wild-type plants. Moreover, the expression levels of jasmonic acid biosynthetic pathway-related genes (NtAOC, NtAOS, NtKAT, NtPACX, NtJMT, NtOPR, and NtLOX), pathogenesis-related genes (NtCHI, NtGlu2, and NtPR-1), and antioxidant stress-related superoxide dismutase genes (NtSOD, NtCu-ZnSOD, and MnSOD) were significantly up-regulated in LrWRKY2 transgenic tobacco lines. Additionally, the transient expression of a hairpin RNA targeting LrWRKY2 increased the susceptibility of L. regale scales to F. oxysporum. Furthermore, an F. oxysporum resistance gene (LrCHI2) encoding a chitinase was isolated from L. regale. An electrophoretic mobility shift assay demonstrated that LrWRKY2 can bind to the LrCHI2 promoter containing the W-box element. Yeast one-hybrid assay results suggested that LrWRKY2 can activate LrCHI2 transcription. An examination of transgenic tobacco transformed with LrWRKY2 and the LrCHI2 promoter revealed that LrWRKY2 activates the LrCHI2 promoter. Therefore, in L. regale, LrWRKY2 is an important positive regulator that contributes to plant defense responses to F. oxysporum by modulating LrCHI2 expression.