Genome-wide analysis of MYB transcription factors and their responses to salt stress in Casuarina equisetifolia

The MYB family and their response to abiotic stress in ginger (Zingiber officinale Roscoe)

BACKGROUND

Zingiber officinale Roscoe, colloquially known as ginger, is a crop of significant medicinal and culinary value that frequently encounters adversity stemming from inhospitable environmental conditions. The MYB transcription factors have garnered recognition for their pivotal role in orchestrating a multitude of plant biological pathways. Nevertheless, the enumeration and characterization of the MYBs within Z. officinale Roscoe remains unknown. This study embarks on a genome-wide scrutiny of the MYB gene lineage in ginger, with the aim of cataloging all ZoMYB genes implicated in the biosynthesis of gingerols and curcuminoids, and elucidating their potential regulatory mechanisms in counteracting abiotic stress, thereby influencing ginger growth and development.

RESULTS

In this study, we identified an MYB gene family comprising 231 members in ginger genome. This ensemble comprises 74 singular-repeat MYBs (1R-MYB), 156 double-repeat MYBs (R2R3-MYB), and a solitary triple-repeat MYB (R1R2R3-MYB). Moreover, a comprehensive analysis encompassing the sequence features, conserved protein motifs, phylogenetic relationships, chromosome location, and gene duplication events of the ZoMYBs was conducted. We classified ZoMYBs into 37 groups, congruent with the number of conserved domains and gene structure analysis. Additionally, the expression profiles of ZoMYBs during development and under various stresses, including ABA, cold, drought, heat, and salt, were investigated in ginger utilizing both RNA-seq data and qRT-PCR analysis.

CONCLUSION

This work provides a comprehensive understanding of the MYB family in ginger and lays the foundation for the future investigation of the potential functions of ZoMYB genes in ginger growth, development and abiotic stress tolerance of ginger.

Genome-wide identification of WRKY transcription factors in Casuarina equisetifolia and the function analysis of CeqWRKY11 in response to NaCl/NaHCO3 stresses

BACKGROUND

Casuarina equisetifolia (C. equisetifolia) is a woody species with many excellent features. It has natural resistance against drought, salt and saline-alkali stresses. WRKY transcription factors (TFs) play significant roles in plant response to abiotic stresses, therefore, molecular characterization of WRKY gene family under abiotic stresses holds great significance for improvement of forest trees through molecular biological tools. At present, WRKY TFs from C. equisetifolia have not been thoroughly studied with respect to their role in salt and saline-alkali stresses response. The current study was conducted to bridge the same knowledge gap.

RESULTS

A total of 64 WRKYs were identified in C. equisetifolia and divided into three major groups i.e. group I, II and III, consisting of 10, 42 and 12 WRKY members, respectively. The WRKY members in group II were further divided into 5 subgroups according to their homology with Arabidopsis counterparts. WRKYs belonging to the same group exhibited higher similarities in gene structure and the presence of conserved motifs. Promoter analysis data showed the presence of various response elements, especially those related to hormone signaling and abiotic stresses, such as ABRE (ABA), TGACG (MeJA), W-box ((C/T) TGAC (T/C)) and TC-rich motif. Tissue specific expression data showed that CeqWRKYs were mainly expressed in root under normal growth conditions. Furthermore, most of the CeqWRKYs were up-regulated by NaCl and NaHCO3 stresses with few of WRKYs showing early responsiveness to both stresses while few others exhibiting late response. Although the expressions of CeqWRKYs were also induced by cold stress, the response was delayed compared with other stresses. Transgenic C. equisetifolia plants overexpressing CeqWRKY11 displayed lower electrolyte leakage, higher chlorophyll content, and enhanced tolerance to both stresses. The higher expression of abiotic stress related genes, especially CeqHKT1 and CeqPOD7, in overexpression lines points to the maintenance of optimum Na+/K+ ratio, and ROS scavenging as possible key molecular mechanisms underlying salt stress tolerance.

CONCLUSIONS

Our results show that CeqWRKYs might be key regulators of NaCl and NaHCO3 stresses response in C. equisetifolia. In addition, positive correlation of CeqWRKY11 expression with increased stress tolerance in C. equisetifolia encourages further research on other WRKY family members through functional genomic tools. The best candidates could be incorporated in other woody plant species for improving stress tolerance.

Comparative transcriptome analysis and identification of candidate R2R3-MYB genes involved in anthraquinone biosynthesis in Rheum palmatum L

BACKGROUND

Rheum palmatum L. has important medicinal value because it contains biologically active anthraquinones. However, the key genes and TFs involved in anthraquinone biosynthesis and regulation in R. palmatum remain unclear.

METHODS

Based on full length transcriptome data, in this study, we screened the differentially expressed genes in the anthraquinone biosynthesis pathway. The R2R3-MYB family genes of R. palmatum were systematically identified based on full-length transcriptome sequencing followed by bioinformatics analyses. The correlation analysis was carried out by using co-expression analysis, protein interaction analysis, and real-time fluorescence quantitative analysis after MeJA treatment. The RpMYB81 and RpMYB98 genes were amplified by RT-PCR, and their subcellular localization and self-activation characteristics were analyzed.

RESULTS

Comparative transcriptome analysis results revealed a total of 3525 upregulated and 6043 downregulated DEGs in the CK versus MeJA group; 28 DEGs were involved in the anthraquinone pathway. Eleven CHS genes that belonged to the PKS family were differentially expressed and involved in anthraquinone biosynthesis. Twelve differentially expressed MYBs genes were found to be co-expressed and interact with CHS genes. Furthermore, 52 MYB genes were identified as positive regulators of anthraquinone biosynthesis and were further characterized. Three MYB genes including RpMYB81, RpMYB98, and RpMYB100 responded to MeJA treatment in R. palmatum, and the levels of these genes were verified by qRT-PCR. RpMYB81 was related to anthraquinone biosynthesis. RpMYB98 had an interaction with genes in the anthraquinone biosynthesis pathway. RpMYB81 and RpMYB98 were mainly localized in the nucleus. RpMYB81 had self-activation activity, while RpMYB98 had no self-activation activity.

CONCLUSION

RpMYB81, RpMYB98, and RpMYB100 were significantly induced by MeJA treatment. RpMYB81 and RpMYB98 are located in the nucleus, and RpMYB81 has transcriptional activity, suggesting that it might be involved in the transcriptional regulation of anthraquinone biosynthesis in R. palmatum.

Advances in regulatory mechanism(s) and biotechnological approaches to modulate nicotine content in tobacco

More than 8 million deaths are caused by tobacco-related diseases every year. A staggering 1.2 million of those fatalities occur due to second-hand smoke exposure among non-smokers, but more than 7 million are due to direct tobacco use among smokers. Nicotine acts as the key ingredient triggering the addiction. The United States Food and Drug Administration (FDA) has classified more than 90 chemical components of tobacco and related smoke as hazardous or potentially hazardous leading to cancer, cardiovascular, respiratory, and reproductive disorders. Hence, reducing nicotine content has been the foremost objective to reduce health and death risks. Therefore, various biotechnological approaches for developing tobacco varieties with low nicotine concentrations are urgently required for the welfare of humankind. In recent years, numerous advancements have been made in nicotine-based tobacco research, suggesting regulatory components involved in nicotine biosynthesis and developing nicotine-less tobacco varieties through biotechnological approaches. This review highlights the various regulatory components and major approaches used to modulate nicotine content in tobacco cultivars.

Genome-wide identification, classification and expression analysis of MYB gene family in coconut (Cocos nucifera L.)

MYB transcription factors regulate the growth, development, and secondary metabolism of plant species. To investigate the origin of color variations in coconut pericarp, we identified and analyzed the MYB gene family present in coconut. According to the sequence of MYB genes in Arabidopsis thaliana, homologous MYB gene sequences were found in the whole genome database of coconut, the conserved sequence motifs within MYB proteins were analyzed by Motif Elicitation (MEME) tool, and the sequences without conservative structure were eliminated. Additionally, we employed RNA-seq technology to generate gene expression signatures of the R2R3-MYB genes across distinctive coconut parts exhibiting diverse colors. To validate these profiles, we conducted quantitative PCR (qPCR). Through comprehensive genome-wide screening, we successfully identified a collection of 179 MYB genes in coconut. Subsequent phylogenetic analysis categorized these 179 coconut MYB genes into 4-subfamilies: 124 R2R3-MYB, 4 3R-MYB types, 4 4R-MYB type, and 47 unknown types. Furthermore, these genes were further divided into 34 subgroups, with 28 of these subgroups successfully classified into known subfamilies found in Arabidopsis thaliana. By mapping the CnMYB genes onto the 16 chromosomes of the coconut genome, we unveiled a collinearity association between them. Moreover, a preservation of gene structure and motif distribution was observed across the CnMYB genes. Our research encompassed a thorough investigation of the R2R3-MYB genes present in the coconut genome, including the chromosomal localization, gene assembly, conserved regions, phylogenetic associations, and promoter cis-acting elements of the studied genes. Our findings revealed a collection of 12 R2R3-MYB candidate genes, namely CnMYB8, CnMYB15, CnMYB27, CnMYB28, CnMYB61, CnMYB63, CnMYB68, CnMYB94, CnMYB101, CnMYB150, CnMYB153, and CnMYB164. These genes showed differential expressions in diverse tissues and developmental stages of four coconut species, such as CnMYB68, CnMYB101, and CnMYB28 exhibited high expression in majority of tissues and coconut species, while CnMYB94 and CnMYB164 showed lower expression. These findings shed light on the crucial functional divergence of CnMYB genes across various coconut tissues, suggesting these genes as promising candidate genes for facilitating color development in this important crop.

Function of the R2R3-MYB Transcription Factors in Dalbergia odorifera and Their Relationship with Heartwood Formation

R2R3-MYB transcription factors (TFs) form one of the most important TF families involved in regulating various physiological functions in plants. The heartwood of Dalbergia odorifera is a kind of high-grade mahogany and valuable herbal medicine with wide application. However, the role of R2R3-MYB genes in the growth and development of D. odorifera, especially their relevance to heartwood formation, has not been revealed. A total of 126 R2R3-MYBs were screened from the D. odorifera genome and named DodMYB1-126 based on their location on 10 chromosomes. The collinearity results showed that purification selection was the main driving force for the evolution of the R2R3-MYB TFs family, and whole genome/fragment replication event was the main form for expanding the R2R3-MYB family, generating a divergence of gene structure and function. Comparative phylogenetic analysis classified the R2R3-MYB TFs into 33 subfamilies. S3-7,10,12-13,21 and N4-7 were extensively involved in the metabolic process; S9,13,16-19,24-25 and N1-3,8 were associated with the growth and development of D. odorifera. Based on the differential transcriptional expression levels of R2R3-MYBs in different tissues, DodMYB32, DodMYB55, and DodMYB89 were tentatively screened for involvement in the regulatory process of heartwood. Further studies have shown that the DodMYB89, localized in the nucleus, has transcriptional activation activity and is involved in regulating the biosynthesis of the secondary metabolites of heartwood by activating the promoters of the structural genes DodI2'H and DodCOMT. This study aimed to comprehensively analyze the functions of the R2R3-MYB TFs and screen for candidate genes that might be involved in heartwood formation of D. odorifera.

Genomic survey of MYB gene family in six pearl millet (Pennisetum glaucum) varieties and their response to abiotic stresses

In addition to their roles in developmental and metabolic processes, MYB transcription factors play crucial roles in plant defense mechanisms and stress responses. A comprehensive analysis of six pearl millet genomes revealed the presence of 1133 MYB genes, which can be classified into four phylogenetically distinct subgroups. The duplication pattern of MYB genes across the pearl millet genomes demonstrates their conserved and similar evolutionary history. Overall, MYB genes were observed to be involved in drought and heat stress responses, with stronger differential expressed observed in root tissues. Multiple analyses indicated that MYB genes mediate abiotic stress responses by modulating abscisic acid-related pathways, circadian rhythms, and histone modification processes. A substantial number of duplicated genes were determined to exhibit differential expression under abiotic stress. The consistent positive expression trend observed in duplicated gene pairs, such as PMA5G04432.1 and PMA2G00728.1, across various abiotic stresses suggests that duplicated MYB genes plays a key role in the evolution of adaptive responses of pearl millet to abiotic stresses.

Transcriptome-wide expression analysis of MYB gene family leads to functional characterization of flavonoid biosynthesis in fruit coloration of Ziziphus Mill

The Ziziphus mauritiana Lam. and Z. jujuba Mill. are the two most economically important members of the genus Ziziphus. The fruit color of Z. mauritiana remains green throughout fruit development in the majority of commercial cultivars, whereas its close relative, Z. jujuba Mill. turns from green to red in all cultivars. However, the lack of transcriptomic and genomic information confines our understanding of the molecular mechanisms underlying fruit coloration in Z. mauritiana (Ber). In the present study, we performed the transcriptome-wide analysis of MYB transcription factors (TFs) genes in Z. mauritiana and Z. jujuba, and identified 56 ZmMYB and 60 ZjMYB TFs in Z. mauritiana and Z. jujuba, respectively. Through transcriptomic expression analysis, four similar MYB genes (ZmMYB/ZjMYB13, ZmMYB/ZjMYB44, ZmMYB/ZjMYB50, and ZmMYB/ZjMYB56) from Z. mauritiana and Z. jujuba were selected as candidate key genes regulating flavonoid biosynthesis. Among these genes, the ZjMYB44 gene was transiently highly expressed in fruit, and flavonoid content accumulation also increased, indicating that this gene can influence flavonoid content during the period of fruit coloration in Z. jujuba. The current study adds to our understanding of the classification of genes, motif structure, and predicted functions of the MYB TFs, as well as identifying MYBs that regulate flavonoid biosynthesis in Ziziphus (Z. mauritiana and Z. jujuba). Based on this information, we concluded that MYB44 is involved in the flavonoids biosynthesis pathway during the fruit coloring of Ziziphus. Our research results provide an important understanding of the molecular mechanism of flavonoid biosynthesis resulting in fruit coloration and laying a foundation for further genetic improvement of fruit color in Ziziphus.

BcMYB111 Responds to BcCBF2 and Induces Flavonol Biosynthesis to Enhance Tolerance under Cold Stress in Non-Heading Chinese Cabbage

Flavonols have been shown to respond to a variety of abiotic stresses in plants, including cold stress. Higher total flavonoid content was found in non-heading Chinese cabbage (NHCC, Brassica campestris (syn. Brassica rapa) ssp. chinensis) after cold stress. A non-targeted metabolome analysis showed a significant increase in flavonol content, including that of quercetin and kaempferol. Here, we found that an R2R3-MYB transcription factor, BcMYB111, may play a role in this process. BcMYB111 was up-regulated in response to cold treatment, with an accompanying accumulation of flavonols. Then, it was found that BcMYB111 could regulate the synthesis of flavonols by directly binding to the promoters of BcF3H and BcFLS1. In the transgenic hairy roots of NHCC or stable transgenic Arabidopsis, overexpression of BcMYB111 increased flavonol synthesis and accumulation, while these were reduced in virus-induced gene silencing lines in NHCC. After cold stress, the higher proline content and lower malondialdehyde (MDA) content showed that there was less damage in transgenic Arabidopsis than in the wild-type (WT). The BcMYB111 transgenic lines performed better in terms of antioxidant capacity because of their lower H₂O₂ content and higher superoxide dismutase (SOD) and peroxidase (POD) enzyme activities. In addition, a key cold signaling gene, BcCBF2, could specifically bind to the DRE element and activate the expression of BcMYB111 in vitro and in vivo. The results suggested that BcMYB111 played a positive role in enhancing the flavonol synthesis and cold tolerance of NHCC. Taken together, these findings reveal that cold stress induces the accumulation of flavonols to increase tolerance via the pathway of BcCBF2-BcMYB111-BcF3H/BcFLS1 in NHCC.

Genome-Wide Identification and Analysis of MYB Transcription Factor Family in Hibiscus hamabo

MYB transcription factors constitute one of the largest gene families in plants and play essential roles in the regulation of plant growth, responses to stress, and a wide variety of physiological and biochemical processes. In this study, 204 MYB proteins (HhMYBs) were identified in the Hibiscus hamabo Sieb. et Zucc (H. hamabo) genome and systematically analyzed based on their genomic sequence and transcriptomic data. The candidate HhMYB proteins and MYBs of Arabidopsis thaliana were divided into 28 subfamilies based on the analysis of their phylogenetic relationships and their motif patterns. Expression analysis using RNA-seq and quantitative real-time PCR (qRT-PCR) indicated that most HhMYBs are differentially regulated under drought and salt stresses. qRT-PCR analysis of seven selected HhMYBs suggested that the HhMYB family may have regulatory roles in the responses to stress and hormones. This study provides a framework for a more comprehensive analysis of the role of MYBs in the response to abiotic stress in H. hamabo.

Comparative analysis of the MYB gene family in seven Ipomoea species

The MYB transcription factors regulate plant growth, development, and defense responses. However, information about the MYB gene family in Ipomoea species is rare. Herein, we performed a comprehensive genome-wide comparative analysis of this gene family among seven Ipomoea species, sweet potato (I. batatas), I. trifida, I. triloba, I. nil, I. purpurea, I. cairica, and I. aquatic, and identified 296, 430, 411, 291, 226, 281, and 277 MYB genes, respectively. The identified MYB genes were classified into five types: 1R-MYB (MYB-related), 2R-MYB (R2R3-MYB), 3R-MYB (R1R2R3-MYB), 4R-MYB, and 5R-MYB, and the MYB-related or R2R3-MYB type was the most abundant MYB genes in the seven species. The Ipomoea MYB genes were classed into distinct subgroups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis. Analysis of gene structure and protein motifs revealed that members within the same phylogenetic group presented similar exon/intron and motif organization. The identified MYB genes were unevenly mapped on the chromosomes of each Ipomoea species. Duplication analysis indicated that segmental and tandem duplications contribute to expanding the Ipomoea MYB genes. Non-synonymous substitution (Ka) to synonymous substitution (Ks) [Ka/Ks] analysis showed that the duplicated Ipomoea MYB genes are mainly under purifying selection. Numerous cis-regulatory elements related to stress responses were detected in the MYB promoters. Six sweet potato transcriptome datasets referring to abiotic and biotic stresses were analyzed, and MYB different expression genes' (DEGs') responses to stress treatments were detected. Moreover, 10 sweet potato MYB DEGs were selected for qRT-PCR analysis. The results revealed that four responded to biotic stress (stem nematodes and Ceratocystis fimbriata pathogen infection) and six responded to the biotic stress (cold, drought, and salt). The results may provide new insights into the evolution of MYB genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Genome-wide characterization and expression analysis of MYB transcription factors in Chrysanthemum nankingense

BACKGROUND

Chrysanthemum is a popular ornamental plant worldwide. MYB (v-myb avian myeloblastosis viral oncogene homolog) transcription factors play an important role in everything from stress resistance to plant growth and development. However, the MYB family of chrysanthemums has not been the subject of a detailed bioinformatics and expression investigation.

RESULTS

In this study, we examined 324 CnMYB transcription factors from Chrysanthemum nankingense genome data, which contained 122 Cn1R-MYB, 183 CnR2R3-MYB, 12 Cn3R-MYB, 2 Cn4R-MYB, and 5 atypical CnMYB. The protein motifs and classification of CnMYB transcription factors were analyzed. Among them, motifs 1, 2, 3, and 4 were found to encode the MYB DNA-binding domain in R2R3-MYB proteins, while in other-MYB proteins, the motifs 1, 2, 3, 4, 5, 6, 7, and 8 encode the MYB DNA-binding domain. Among all CnMYBs, 44 genes were selected due to the presence of CpG islands, while methylation is detected in three genes, including CnMYB9, CnMYB152, and CnMYB219. We analyzed the expression levels of each CnMYB gene in ray floret, disc floret, flower bud, leaf, stem, and root tissues. Based on phylogenetic analysis and gene expression analysis, three genes appeared likely to control cellulose and lignin synthesis in stem tissue, and 16 genes appeared likely to regulate flowering time, anther, pollen development, and flower color. Fifty-one candidate genes that may be involved in stress response were identified through phylogenetic, stress-responseve motif of promoter, and qRT-PCR analyses. According to genes expression levels under stress conditions, six CnMYB genes (CnMYB9, CnMYB172, CnMYB186, CnMYB199, CnMYB219, and CnMYB152) were identified as key stress-responsive genes.

CONCLUSIONS

This research provides useful information for further functional analysis of the CnMYB gene family in chrysanthemums, as well as offers candidate genes for further study of cellulose and lignin synthesis, flowering traits, salt and drought stress mechanism.

Full-length transcriptome combined with RNA sequence analysis of Fraxinus chinensis

BACKGROUND

The dry root or stem bark of Fraxinus chinensis is a famous herb Qin Pi which is known for its anti-inflammatory, analgesic, anti-tumor, liver protective and diuretic pharmacological effects, the fundamental chemical components are coumarin, phenylethanol glycosides and flavonoids. However, it is difficult to clarify the secondary metabolite synthesis pathway and key genes involved in the pathway because of lack genome information of Fraxinus chinensis.

OBJECTIVE

To generate a complete transcriptome of Fraxinus chinensis and to clarify the differentially expressed genes (DEGs) in leaves and stem barks.

METHODS

In this study, full-length transcriptome analysis and RNA-Seq were combined to characterize Fraxinus chinensis transcriptome.

RESULTS

A total of 69,145 transcripts were acquired and regarded as reference transcriptome, 67,441 transcripts (97.47%) were annotated to NCBI non-redundant protein (Nr), SwissProt, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and eukaryotic orthologous groups (KOG) databases. A total of 18,917 isoforms were annotated to KEGG database and classified to 138 biological pathways. In total, 10,822 simple sequence repeat (SSRs) and 11,319 resistance (R) gene were classified to 18 types, and 3947 transcription factors (TFs) were identified in full-length transcriptome analysis. Additionally, 15,095 DEGs were detected by RNA-seq in leaves and barks, including 4696 significantly up-regulated and 10,399 significantly down-regulated genes. And 254 transcripts were annotated into phenylpropane metabolism pathway containing 86 DEGs and ten of these enzyme genes were verified by qRT-PCR.

CONCLUSION

It laid the foundation for further exploration of the biosynthetic pathway of phenylpropanoids and related key enzyme genes.

Genome-wide identification of bHLH transcription factors and their response to salt stress in Cyclocarya paliurus

As a highly valued and multiple function tree species, the leaves of Cyclocarya paliurus are enriched in diverse bioactive substances with healthy function. To meet the requirement for its leaf production and medical use, the land with salt stress would be a potential resource for developing C. paliurus plantations due to the limitation of land resources in China. The basic helix-loop-helix (bHLH) transcription factor protein family, the second largest protein family in plants, has been found to play essential roles in the response to multiple abiotic stresses, especially salt stress. However, the bHLH gene family in C.paliurus has not been investigated. In this study, 159 CpbHLH genes were successfully identified from the whole-genome sequence data, and were classified into 26 subfamilies. Meanwhile, the 159 members were also analyzed from the aspects of protein sequences alignment, evolution, motif prediction, promoter cis-acting elements analysis and DNA binding ability. Based on transcriptome profiling under a hydroponic experiment with four salt concentrations (0%, 0.15%, 0.3%, and 0.45% NaCl), 9 significantly up- or down-regulated genes were screened, while 3 genes associated with salt response were selected in term of the GO annotation results. Totally 12 candidate genes were selected in response to salt stress. Moreover, based on expression analysis of the 12 candidate genes sampled from a pot experiment with three salt concentrations (0%, 0.2% and 0.4% NaCl), CpbHLH36/68/146 were further verified to be involved in the regulation of salt tolerance genes, which is also confirmed by protein interaction network analysis. This study was the first analysis of the transcription factor family at the genome-wide level of C. paliurus, and our findings would not only provide insight into the function of the CpbHLH gene family members involved in salt stress but also drive progress in genetic improvement for the salt tolerance of C. paliurus.

Characterization of HAK protein family in Casuarina equisetifolia and the positive regulatory role of CeqHAK6 and CeqHAK11 genes in response to salt tolerance

The potassium transporter group of the HAK/KUP/KT (high-affinity K⁺)/KUP (K⁺ uptake)/KT (K⁺ transporter) family plays a crucial role in plant growth and development as well as in environmental adaptation such as tolerance to salt stress. HAK/KUP/KT genes and their functions have been characterized for a number of plant species, but they remain unknown for Casuarina equisetifolia, an important tree species for coastal protection in southern China and many other countries. In this study, 25 HAK genes were identified in the C. equisetifolia genome. Their gene structure, conserved motif, phylogeny, and expression were comprehensively and systematically analyzed to understand their functions. All HAK genes were relatively conserved and could be divided into four clusters. The expression level of two particular genes, CeqHAK11 and CeqHAK6, increased significantly with the duration of salt treatment. To further elucidated their function in response to salt stress, subcellular localization, and their functional analysis were developed. Results revealed that CeqHAK11 and CeqHAK6 were localized on the plasma membrane, which mainly mediated high-affinity K⁺ uptake. Overexpression of CeqHAK6 or CeqHAK11 in Arabidopsis showed higher germination and survival rates and longer root length than wild-type (WT) under salt stress, suggesting that both genes improve tolerance to salt stress. Moreover, CeqHAK6 and CeqHAK11 improved their ability to tolerate salt stress by increasing the K⁺/Na⁺ ratio and antioxidant enzyme activities (CAT, POD, and SOD), and decreasing reactive oxygen species (ROS) accumulation. Consequently, CeqHAK6 and CeqHAK11 were verified as potassium transport proteins and could be applied for further molecular breeding for salt tolerance in C. equisetifolia or other crops to increasing salt tolerance.

Comprehensive identification and expression analyses of the SnRK gene family in Casuarina equisetifolia in response to salt stress

BACKGROUND

Sucrose nonfermenting-1 (SNF1)-related protein kinases (SnRKs) play crucial roles in plant signaling pathways and stress adaptive responses by activating protein phosphorylation pathways. However, there have been no comprehensive studies of the SnRK gene family in the widely planted salt-tolerant tree species Casuarina equisetifolia. Here, we comprehensively analyze this gene family in C. equisetifolia using genome-wide identification, characterization, and profiling of expression changes in response to salt stress.

RESULTS

A total of 26 CeqSnRK genes were identified, which were divided into three subfamilies (SnRK1, SnRK2, and SnRK3). The intron-exon structures and protein‑motif compositions were similar within each subgroup but differed among groups. Ka/Ks ratio analysis indicated that the CeqSnRK family has undergone purifying selection, and cis-regulatory element analysis suggested that these genes may be involved in plant development and responses to various environmental stresses. A heat map was generated using quantitative real‑time PCR (RT-qPCR) data from 26 CeqSnRK genes, suggesting that they were expressed in different tissues. We also examined the expression of all CeqSnRK genes under exposure to different salt concentrations using RT-qPCR, finding that most CeqSnRK genes were regulated by different salt treatments. Moreover, co-expression network analysis revealed synergistic effects among CeqSnRK genes.

CONCLUSIONS

Several CeqSnRK genes (CeqSnRK3.7, CeqSnRK3.16, CeqSnRK3.17) were up-regulated following salt treatment. Among them, CeqSnRK3.16 expression was significantly up-regulated under various salt treatments, identifying this as a candidate gene salt stress tolerance gene. In addition, CeqSnRK3.16 showed significant expression change correlations with multiple genes under salt stress, indicating that it might exhibit synergistic effects with other genes in response to salt stress. This comprehensive analysis will provide a theoretical reference for CeqSnRK gene functional verification and the role of these genes in salt tolerance.

Integrated metabolomic and transcriptomic analysis reveals the role of phenylpropanoid biosynthesis pathway in tomato roots during salt stress

As global soil salinization continues to intensify, there is a need to enhance salt tolerance in crops. Understanding the molecular mechanisms of tomato (Solanum lycopersicum) roots' adaptation to salt stress is of great significance to enhance its salt tolerance and promote its planting in saline soils. A combined analysis of the metabolome and transcriptome of S. lycopersicum roots under different periods of salt stress according to changes in phenotypic and root physiological indices revealed that different accumulated metabolites and differentially expressed genes (DEGs) associated with phenylpropanoid biosynthesis were significantly altered. The levels of phenylpropanoids increased and showed a dynamic trend with the duration of salt stress. Ferulic acid (FA) and spermidine (Spd) levels were substantially up-regulated at the initial and mid-late stages of salt stress, respectively, and were significantly correlated with the expression of the corresponding synthetic genes. The results of canonical correlation analysis screening of highly correlated DEGs and construction of regulatory relationship networks with transcription factors (TFs) for FA and Spd, respectively, showed that the obtained target genes were regulated by most of the TFs, and TFs such as MYB, Dof, BPC, GRAS, and AP2/ERF might contribute to the regulation of FA and Spd content levels. Ultimately, FA and Spd attenuated the harm caused by salt stress in S. lycopersicum, and they may be key regulators of its salt tolerance. These findings uncover the dynamics and possible molecular mechanisms of phenylpropanoids during different salt stress periods, providing a basis for future studies and crop improvement.

Genome-Wide Identification and Expression Analysis of the Zinc Finger Protein Gene Subfamilies under Drought Stress in Triticum aestivum

The zinc finger protein (ZFP) family is one of plants’ most diverse family of transcription factors. These proteins with finger-like structural domains have been shown to play a critical role in plant responses to abiotic stresses such as drought. This study aimed to systematically characterize Triticum aestivum ZFPs (TaZFPs) and understand their roles under drought stress. A total of 9 TaC2H2, 38 TaC3HC4, 79 TaCCCH, and 143 TaPHD were identified, which were divided into 4, 7, 12, and 14 distinct subgroups based on their phylogenetic relationships, respectively. Segmental duplication dominated the evolution of four subfamilies and made important contributions to the large-scale amplification of gene families. Syntenic relationships, gene duplications, and Ka/Ks result consistently indicate a potential strong purifying selection on TaZFPs. Additionally, TaZFPs have various abiotic stress-associated cis-acting regulatory elements and have tissue-specific expression patterns showing different responses to drought and heat stress. Therefore, these genes may play multiple functions in plant growth and stress resistance responses. This is the first comprehensive genome-wide analysis of ZFP gene families in T. aestivum to elucidate the basis of their function and resistance mechanisms, providing a reference for precise manipulation of genetic engineering for drought resistance in T. aestivum.

Overexpression of McHB7 Transcription Factor from Mesembryanthemum crystallinum Improves Plant Salt Tolerance

Mesembryanthemum crystallinum (common ice plant) is one of the facultative halophyte plants, and it serves as a model for investigating the molecular mechanisms underlying its salt stress response and tolerance. Here we cloned one of the homeobox transcription factor (TF) genes, McHB7, from the ice plant, which has 60% similarity with the Arabidopsis AtHB7. Overexpression of the McHB7 in Arabidopsis (OE) showed that the plants had significantly elevated relative water content (RWC), chlorophyll content, superoxide dismutase (SOD), and peroxidase (POD) activities after salt stress treatment. Our proteomic analysis identified 145 proteins to be significantly changed in abundance, and 66 were exclusively increased in the OE plants compared to the wild type (WT). After salt treatment, 979 and 959 metabolites were significantly increased and decreased, respectively, in the OE plants compared to the WT. The results demonstrate that the McHB7 can improve photosynthesis, increase the leaf chlorophyll content, and affect the TCA cycle by regulating metabolites (e.g., pyruvate) and proteins (e.g., citrate synthase). Moreover, McHB7 modulates the expression of stress-related proteins (e.g., superoxide dismutase, dehydroascorbate reductase, and pyrroline-5-carboxylate synthase B) to scavenge reactive oxygen species and enhance plant salt tolerance.

Transcriptome analysis of Tamarix ramosissima leaves in response to NaCl stress

Halophyte Tamarix ramosissima. Lcdcb (T. ramosissima) are known as the representative of Tamarix plants that are widely planted in salinized soil. However, molecular mechanisms towards salt tolerance and adaptation are largely rare. In this study, we carried out RNA-sequence and transcriptome analysis of T. ramosissima in response to NaCl stress, screened differentially expressed genes (DEGs) and further verified by qRT-PCR. Results showed that 105702 unigenes were spliced from the raw data of transcriptome sequencing, where 54238 unigenes were retrieved from KEGG, KOG, NR, and SwissProt. After 48 hours of NaCl treatment, the expression levels of 6374 genes were increased, and 5380 genes were decreased in leaves. After 168 hours, the expression levels of 3837 genes were up-regulated and 7808 genes were down-regulated. In particular, 8 transcription factors annotated to the KEGG Pathway were obtained, involving the WRKY and bZIP transcription family. In addition, KEGG pathway annotation showed that expression of 39 genes involved in ROS scavenging mechanisms were significantly changed, in which 21 genes were up-regulated and 18 genes were down-regulated after 48 hours as well as 15 genes were up-regulated and 24 genes were down-regulated after 168h. Simultaneously, the enzyme activities of SOD and POD were significantly enhanced under NaCl treatment, but the enzyme activity of CAT was not significantly enhanced. Moreover, WRKY, MYB and bZIP may participate in the process of salt resistance in T. ramosissima. This study provides gene resources and a theoretical basis for further molecular mechanisms of salt tolerance in T. ramosissima.

Genome-wide Identification and Expression Analysis of RcMYB Genes in Rhodiola crenulata

Modern research has proved that the main medicinal component of Rhodiola crenulata, which has a wide range of medicinal value, is its secondary metabolite salidroside. The MYB transcription factor family is widely involved in biosynthesis of second metabolism and other roles in the stress response in plants, so a genome-wide identification and analysis for this family in R. crenulata is worth conducting. In this research, genome-wide analysis identified 139 MYB genes based on conserved domains in the R. crenulata genome, and 137 genes were used to construct a phylogenetic tree and modified with expression files to reveal evolutionary characteristics. Physical and chemical characteristics, gene structure, and conserved motif analysis were also used to further analyze RcMYBs. Additionally, cis-acting elements related to transcription, hormone, and MYB binding were found in the promoter region of the selected RcMYBs. Four RcMYBs were cloned, sequenced, and their gene expression pattern was analyzed for further analysis of their functions. The research results lay the foundation for further research on the function of RcMYB and R. crenulata.

Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Salt Tolerance in Cyclocarya paliurus

R2R3-MYB transcription factors are most abundant in the MYB superfamily, while the R2R3-MYB genes play an important role in plant growth and development, especially in response to environmental stress. Cyclocarya paliurus is a multifunction tree species, and the existing resources cannot meet the requirement for its leaf production and medical use. Therefore, lands with some environmental stresses would be potential sites for developing C. paliurus plantations. However, the function of R2R3-MYB genes in C.paliurus in response to environmental stress remains unknown. In this study, to identify the roles of R2R3-MYB genes associated with salt stress response, 153 CpaMYB genes and their corresponding protein sequences were identified from the full-length transcriptome. Based on the comparison with MYB protein sequences of Arabidopsis thaliana, 69 R2R3-MYB proteins in C. paliurus were extracted for further screening combined with conserved functional domains. Furthermore, the MYB family members were analyzed from the aspects of protein sequences alignment, evolution, motif prediction, promoter cis-acting element analysis, and gene differential expression under different salt treatments using both a pot experiment and hydroponic experiment. The results showed that the R2R3-MYB genes of C.paliurus conserved functional domains, whereas four R2R3-MYB genes that might respond to salt stress via regulating plant hormone signals were identified in this study. This work provides a basis for further functional characterization of R2R3-MYB TFs in C. paliurus.

Genome-Wide Identification and Comparative Analysis of WRKY Transcription Factors Related to Momilactone Biosynthesis in Calohypnum plumiforme

Momilactones are diterpenoid phytoalexins with allelopathic functions, which have been found in the widely distributed bryophyte Calohypnum plumiforme. Clustered genes containing CpDTC1/HpDTC1, CpCYP970A14, CpCYP964A1, and CpMAS are involved in momilactone biosynthesis. Besides, momilactone concentration in C. plumiforme is affected by heavy metal treatment such as CuCl2. However, transcription factors which might regulate momilactone biosynthesis are unclear. WRKY transcription factors (TFs) regulate phytoalexin biosynthesis in many plant species. In this study, a systematic analysis of the WRKY TFs was performed according to the C. plumiforme genome. A total of 19 CpWRKY genes were identified and categorized into five subgroups based on their phylogenetic relationship. Conserved domain and motif analysis suggested that the WRKY domain was highly conserved, but there were some variations. Cis-acting elements and binding sites analysis implied that CpWRKY genes might be induced by stress and further regulate the biosynthesis of momilactones. Our study lays a foundation for further functional characterization of the candidate CpWRKY genes involved in the regulation of momilactone biosynthesis, and provides new strategies for increasing momilactone production.

MYB transcription factor family in sweet cherry (Prunus avium L.): genome-wide investigation, evolution, structure, characterization and expression patterns

BACK GROUND

MYB Transcription factors (TFs) are most imperative and largest gene family in plants, which participate in development, metabolism, defense, differentiation and stress response. The MYB TFs has been studied in various plant species. However, comprehensive studies of MYB gene family in the sweet cherry (Prunus avium L.) are still unknown.

RESULTS

In the current study, a total of 69 MYB genes were investigated from sweet cherry genome and classified into 28 subfamilies (C1-C28 based on phylogenetic and structural analysis). Microcollinearity analysis revealed that dispersed duplication (DSD) events might play an important role in the MYB genes family expansion. Chromosomal localization, the synonymous (Ks) and nonsynonymous (Ka) analysis, molecular characteristics (pI, weight and length of amino acids) and subcellular localization were accomplished using several bioinformatics tools. Furthermore, the members of distinct subfamilies have diverse cis-acting regions, conserved motifs, and intron-exon architectures, indicating functional heterogeneity in the MYB family. Moreover, the transcriptomic data exposed that MYB genes might play vital role in bud dormancy. The quantitative real-time qRT-PCR was carried out and the expression pattern indicated that MYB genes significantly expressed in floral bud as compared to flower and fruit.

CONCLUSION

Our comprehensive findings provide supportive insights into the evolutions, expansion complexity and functionality of PavMYB genes. These PavMYB genes should be further investigated as they seem to be brilliant candidates for dormancy manipulation in sweet cherry.

Genome-Wide Identification and Expression Analysis of the R2R3-MYB Transcription Factor Family Revealed Their Potential Roles in the Flowering Process in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is a productive fruit crop with high nutritional and medical value in tropical and subtropical regions. The MYB gene family is one of the most widespread plant transcription factor (TF) families participating in the flowering regulation. However, little is known about the MYB TFs involved in the flowering process in longan and its regulatory network. In this study, a total of 119 DlR2R3-MYB genes were identified in the longan genome and were phylogenetically grouped into 28 subgroups. The groupings were supported by highly conserved gene structures and motif composition of DlR2R3-MYB genes in each subgroup. Collinearity analysis demonstrated that segmental replications played a more crucial role in the expansion of the DlR2R3-MYB gene family compared to tandem duplications, and all tandem/segmental duplication gene pairs have evolved under purifying selection. Interspecies synteny analysis among longan and five representative species implied the occurrence of gene duplication events was one of the reasons contributing to functional differentiation among species. RNA-seq data from various tissues showed DlR2R3-MYB genes displayed tissue-preferential expression patterns. The pathway of flower development was enriched with six DlR2R3-MYB genes. Cis-acting element prediction revealed the putative functions of DlR2R3-MYB genes were related to the plant development, phytohormones, and environmental stresses. Notably, the orthologous counterparts between Arabidopsis and longan R2R3-MYB members tended to play conserved roles in the flowering regulation and stress responses. Transcriptome profiling on off-season flower induction (FI) by KClO₃ indicated two up-regulated and four down-regulated DlR2R3-MYB genes involved in the response to KClO₃ treatment compared with control groups. Additionally, qRT-PCR confirmed certain genes exhibited high expression in flowers/flower buds. Subcellular localization experiments revealed that three predicted flowering-associated MYB proteins were localized in the nucleus. Future functional studies on these potential candidate genes involved in the flowering development could further the understanding of the flowering regulation mechanism.

Integrative Analysis of the GRAS Genes From Chinese White Pear (Pyrus bretschneideri): A Critical Role in Leaf Regeneration

GRAS is a transcription regulator factor, which plays an important role in plant growth and development. Previous analyses found that several GRAS functions have been identified, such as axillary bud meristem formation, radial root elongation, gibberellin signaling, light signaling, and abiotic stress. The GRAS family has been comprehensively evaluated in several species. However, little finding is on the GRAS transcription factors (TFs) in Chinese white pear. In this study, 99 PbGRAS were systemically characterized and renamed PbGRAS1 to PbGRAS99 according to their chromosomal localizations. Phylogenetic analysis and structural features revealed that could be classified into eight subfamilies (LISCL, Ls, SHR, HAM, SCL, PAT, SCR, and DELLA). Further analysis of introns/exons and conserved motifs revealed that they are diverse and functionally differentiated in number and structure. Synteny analysis among Pyrus bretschenedri, Prunus mume, Prunus avium, Fragaria vesca, and Prunus persica showed that GRAS duplicated regions were more conserved. Dispersed duplication events are the most common mechanism and may play a crucial role in the expansion of the GRAS gene family. In addition, cis-acting elements of the PbGRAS gene were found in promoter regions associated with hormone and environmental stress responses. Notably, the expression pattern detected by qRT-PCR indicated that PbGRAS genes were differentially expressed under gibberellin (GA), abscisic acid (ABA), and auxin (IAA) conditions, which are responsive to abiotic stress. PbGRAS89 and PbGRAS99 were highly expressed at different stages of hormone treatment and may play important role in leaf development. Therefore, we selected PbGRAS89 and PbGRAS99 to clone and construct pCAMBIA1301-PbGRAS89, 99 and transferred them into Arabidopsis thaliana. Finally, we observed and compared the changes of overexpressed plants and wild-type plants during regeneration. This method was used to analyze their roles in leaf regeneration of Chinese white pear. In addition, we also constructed pCAMBIA1305-PbGRAS89, 99, and transferred them into onion cells to determine the subcellular localization. Subcellular localization experiments showed that PbGRAS89 and PbGRAS99 were localized in the nucleus. In summary, the results of this study indicate that PbGRAS89 and PbGRAS99 are mainly responsible for leaf regeneration of Chinese white pear, which plays a positive role in callus formation and provides rich resources for studying GRAS gene functions.

Genome-wide identification of myeloblastosis gene family and its response to cadmium stress in Ipomoea aquatica

The myeloblastosis (MYB) proteins perform key functions in mediating cadmium (Cd) tolerance of plants. Ipomoea aquatica has strong adaptability to Cd Stress, while the roles of the I. aquatica MYB gene family with respect to Cd stress are still unclear. Here, we identified a total of 183 MYB genes in the I. aquatica genome (laMYB), which were classified into 66 1R-type IaMYB, 112 2R-type IaMYB, four 3R-type IaMYB, and one 4R-type IaMYB based on the number of the MYB repeat in each gene. The analysis of phylogenetic tree indicated that most of IaMYB genes are associated with the diverse biological processes including defense, development and metabolism. Analysis of sequence features showed that the IaMYB genes within identical subfamily have the similar patterns of the motif distributions and gene structures. Analysis of gene duplication events revealed that the dispersed duplication (DSD) and whole-genome duplication (WGD) modes play vital roles in the expansion of the IaMYB gene family. Expression profiling manifests that approximately 20% of IaMYB genes had significant role in the roots of I. aquatica under Cd stress. Promoter profiling implied that the differentially expressed genes might be induced by environmental factors or inherent hormones and thereby execute their function in Cd response. Remarkably, the 2R-type IaMYB157 with abundant light-responsive element G-box and ABA-responsive element ABRE in its promoter region exhibited very strong response to Cd stress. Taken together, our findings provide an important candidate IaMYB gene for further deciphering the molecular regulatory mechanism in plant with respect to Cd stress.

Transcriptome Analysis Reveals Genes Respond to Chlorophyll Deficiency in Green and Yellow Leaves of Chrysanthemum morifolium Ramat

Chlorophyll is vital for photosynthesis to produce sugars and other useful biochemical products in green plants. However, the molecular effects of chlorophyll deficiency in Chrysanthemum are largely unknown. In this study, we identified a bud sport mutant chrysanthemum belonging to the variety ‘Nannong Binyun’, which has yellow branches. Plant physiological studies have shown that the yellow color is revealed due to chlorophyll loss. RNA extracts of yellow and green tissues were analyzed using high-throughput RNA-sequencing, and a total of 11,649 tissue enriched unigenes that respond to chlorophyll deficiency were identified, including 4803 unigenes upregulated in yellow tissues and 6846 unigenes in green tissues. GO analysis revealed that these tissue-enriched genes may involve in the physiological processes of chlorophyll accumulation and photosynthesis. In addition, many DEGs from the families of AP2-EREBP, bHLH, MYB, and FAR1 that are associated with plant development and stress response were detected. Our study found that most of the genes from the GRAS family were downregulated in yellow leaves, indicating their putative roles in stem cell maintenance and possible contribution to leaf size determination.

Genome-Wide Identification of NAC Transcription Factor Family in Juglans mandshurica and Their Expression Analysis during the Fruit Development and Ripening

The NAC (NAM, ATAF and CUC) gene family plays a crucial role in the transcriptional regulation of various biological processes and has been identified and characterized in multiple plant species. However, genome-wide identification of this gene family has not been implemented in Juglans mandshurica, and specific functions of these genes in the development of fruits remain unknown. In this study, we performed genome-wide identification and functional analysis of the NAC gene family during fruit development and identified a total of 114 JmNAC genes in the J. mandshurica genome. Chromosomal location analysis revealed that JmNAC genes were unevenly distributed in 16 chromosomes; the highest numbers were found in chromosomes 2 and 4. Furthermore, according to the homologues of JmNAC genes in Arabidopsis thaliana, a phylogenetic tree was constructed, and the results demonstrated 114 JmNAC genes, which were divided into eight subgroups. Four JmNAC gene pairs were identified as the result of tandem duplicates. Tissue-specific analysis of JmNAC genes during different developmental stages revealed that 39 and 25 JmNAC genes exhibited upregulation during the mature stage in walnut exocarp and embryos, indicating that they may serve key functions in fruit development. Furthermore, 12 upregulated JmNAC genes were common in fruit ripening stage in walnut exocarp and embryos, which demonstrated that these genes were positively correlated with fruit development in J. mandshurica. This study provides new insights into the regulatory functions of JmNAC genes during fruit development in J. mandshurica, thereby improving the understanding of characteristics and evolution of the JmNAC gene family.