Systematic Analysis of MYB Family Genes in Potato and Their Multiple Roles in Development and Stress Responses

Transcriptome Analysis Reveals Key Genes Involved in Trichome Formation in Pepper (Capsicum annuum L.)

Trichomes are specialized organs located in the plant epidermis that play important defense roles against biotic and abiotic stresses. However, the mechanisms regulating the development of pepper epidermal trichomes and the related regulatory genes at the molecular level are not clear. Therefore, we performed transcriptome analyses of A114 (less trichome) and A115 (more trichome) to dig deeper into the genes involved in the regulatory mechanisms of epidermal trichome development in peppers. In this study, the epidermal trichome density of A115 was found to be higher by phenotypic observation and was highest in the leaves at the flowering stage. A total of 39,261 genes were quantified by RNA-Seq, including 11,939 genes not annotated in the previous genome analysis and 18,833 differentially expressed genes. Based on KEGG functional enrichment, it was found that DEGs were mainly concentrated in three pathways: plant-pathogen interaction, MAPK signaling pathway-plant, and plant hormone signal transduction. We further screened the DEGs associated with the development of epidermal trichomes in peppers, and the expression of the plant signaling genes GID1B-like (Capana03g003488) and PR-6 (Capana09g001847), the transcription factors MYB108 (Capana05g002225) and ABR1-like (Capana04g001261), and the plant resistance genes PGIP-like (Capana09g002077) and At5g49770 (Capana08g001721) in the DEGs were higher at A115 compared to A114, and were highly expressed in leaves at the flowering stage. In addition, based on the WGCNA results and the establishment of co-expression networks showed that the above genes were highly positively correlated with each other. The transcriptomic data and analysis of this study provide a basis for the study of the regulatory mechanisms of pepper epidermal trichomes.

Genome-Wide Identification and Analysis of Catharanthus roseus Receptor-like Kinase 1-like Proteins in Eggplant

As an important member of the plant receptor-like kinases, Catharanthus roseus receptor-like kinase 1-like (CrRLK1L) plays vital roles in plant growth and development, as well as biotic and abiotic stress response. Numerous CrRLK1Ls have been identified and analyzed in various plant species, while our knowledge about eggplant (Solanum melongena L.) CrRLK1Ls is still scarce. Utilizing state-of-the-art genomic data, we conducted the first genome-wide identification and analysis of CrRLK1L proteins in eggplant. In this study, 32 CrRLK1L proteins were identified and analyzed in eggplant. A subsequent gene structure and protein domain analysis showed that the identified eggplant CrRLK1Ls possessed typical features of CrRLK1Ls. A subcellular localization prediction demonstrated that these proteins mostly localized on the plasma membrane. A collinearity analysis showed that some eggplant CrRLK1L genes had predicted intraspecies or interspecies evolutionary duplication events. Promoter analysis suggests that eggplant CrRLK1Ls may be involved in plant hormone signaling, host-pathogen interactions, and environmental responses. Based on transcriptomic gene expression analysis, it is indicated that eggplant CrRLK1Ls may be involved in the resistance response of eggplant to Botrytis cinerea. Together, these results will give us a theoretical foundation and guidance for elaborating the biological functions of CrRLK1Ls in eggplant growth, development, and resistance response.

In-silico mining and characterization of MYB family genes in wilt-resistant hybrid guava (Psidium guajava × Psidium molle)

BACKGROUND

The MYB family is one of the most significant groups of transcription factors in plants. However, several MYBs have been linked to secondary metabolism and are important for determining the color of fruit's peel and pulp. Despite being a substantial fruit crop in tropical and subtropical areas of the world, wilt-resistant hybrid guava (Psidium guajava × Psidium molle; PGPM) has not yet been the subject of a thorough examination. This study's goal was to assess the expression of MYB in guava fruit pulp, roots, and seeds to predict its function by in silico analysis of the guava root transcriptome data.

RESULTS

In the current study, we have mined the MYBs family of MYB genes from the transcriptome of the PGPM guava root. We have mined 15 distinct MYB transcription factor genes/transcripts viz MYB3, MYB4, MYB23, MYB86, MYB90, MYB308, MYB5, MYB82, MYB114, MYB6, MYB305, MYB44, MYB51, MYB46, and MYB330. From the analyses, it was found that R2-MYB and R3-MYB domains are conserved in all known guava MYB proteins. The expression of six different MYB TFs was examined using semi-quantitative RT-PCR in "Shweta" pulp (white colour pulp), "Lalit" pulp (red color pulp), "Lalit" root, and "Lalit" seed.

CONCLUSION

There were 15 MYB family members observed in guava. They were unequally distributed across the chromosomes, most likely as a result of gene duplication. Additionally, the expression patterns of the particular MYBs showed that MYB may be involved in the control of wilt, fruit ripening, seed development, and root development. Our results allow for a more thorough functional characterization of the guava MYB family genes and open the door to additional research into one essential MYB transcription factor family of genes and its involvement in the growth and ripening of guava fruit.

Over-expression of CcMYB24, encoding a R2R3-MYB transcription factor from a high-leaf-number mutant of Cymbidium, increases the number of leaves in Arabidopsis

Ornamental foliage plants have long been cultivated for their attractive leaves. Variation in leaf traits of ornamental foliage plants is one of the goals in breeding. MYB transcription factors regulate many aspects of leaf development, and thus influence morphological traits of leaves. However, little is known about the function of MYB transcription factors in leaf development of Cymbidium, one of the most economically important ornamental plants in the world. In the present study, a MYB transcription factor, CcMYB24, was identified and the corresponding gene cloned from a new orchid mutant, TRIR-2, which produces more leaves than control plants. The CcMYB24 showed a higher expression level in 'TRIR-2' than in control plants, and the protein was located in the nucleus. The sequence of CcMYB24 showed a high similarity with RAX2-like genes which belong to the R2R3-MYB gene family in other Cymbidium plants. Overexpression of CcMYB24 resulted in a phenotype with an increased number of leaves, elevated chlorophyll content, and decreased contents of carotenoids and flavonoids in Arabidopsis. These results provide functional evidence for the role of CcMYB24 in promoting the production of leaves in 'TRIR-2'. Understanding the role of CcMYB24 in Cymbidium will be beneficial for the molecular breeding of ornamental foliage plants.

SlMYB1 regulates the accumulation of lycopene, fruit shape, and resistance to Botrytis cinerea in tomato

Fruit lycopene, shape, and resistance are essential traits in vegetables whose final product is fruit, and they are also closely related to and strictly regulated by multiple transcription factors. Lycopene, which cannot be synthesized by the human body and can only be ingested from the outside, was important in maintaining human health. During fruit ripening and post-harvest, tomato plants face a variety of biotic or abiotic stresses, which might inflict great damage to fruit quality due to its flat shape and pointed tip during storage and transportation. Therefore, there is an urgent need for key molecular switches to simultaneously improve fruit lycopene and resistance to biotic stress during ripening. Here, we identified the MYB transcription factor SlMYB1 in tomato plants which could bind to the promoters of lycopene synthesis-related genes, SlLCY1, SlPSY2, and the pathogen-related gene SlPR5 directly, to regulate the fruit lycopene and resistance to Botrytis cinerea in tomato. In addition to regulating lycopene synthesis, SlMYB1 also regulates the content of soluble sugar, soluble protein and flavonoid in tomato. What's more, SlMYB1 could regulate the tomato fruit shape, making it smoother or flatter to prevent skin damage caused by vibration on fruits. RNA sequencing (RNA-seq) further showed that SlMYB1 fruit-specific expression lines had multiple differentially expressed genes compared with those from wild-type plants, suggesting that SlMYB1 might have multiple roles in fruit nutritional quality control and resistance to stresses, which is a rare occurrence in previous studies. In summary, our results revealed that SlMYB1 was an essential multi-functional transcription factor that could regulate the lycopene and resistance to Botrytis cinerea, and change the shape of fruit in tomato plants.

Identification and Analysis of the CBF Gene Family in Three Species of Acer under Cold Stress

The C-Repeat Binding Factor (CBF) gene family has been identified and characterized in multiple plant species, and it plays a crucial role in responding to low temperatures. Presently, only a few studies on tree species demonstrate the mechanisms and potential functions of CBFs associated with cold resistance, while our study is a novel report on the multi-aspect differences of CBFs among three tree species, compared to previous studies. In this study, genome-wide identification and analysis of the CBF gene family in Acer truncatum, Acer pseudosieboldianum, and Acer yangbiense were performed. The results revealed that 16 CBF genes (five ApseCBFs, four AcyanCBFs, and seven AtruCBFs) were unevenly distributed across the chromosomes, and most CBF genes were mapped on chromosome 2 (Chr2) and chromosome 11 (Chr11). The analysis of phylogenetic relationships, gene structure, and conserved motif showed that 16 CBF genes could be clustered into three subgroups; they all contained Motif 1 and Motif 5, and most of them only spanned one exon. The cis-acting elements analysis showed that some CBF genes might be involved in hormone and abiotic stress responsiveness. In addition, CBF genes exhibited tissue expression specificity. High expressions of ApseCBF1, ApseCBF3, AtruCBF1, AtruCBF4, AtruCBF6, AtruCBF7, and ApseCBF3, ApseCBF4, ApseCBF5 were detected on exposure to low temperature for 3 h and 24 h. Low expressions of AtruCBF2, AtruCBF6, AtruCBF7 were detected under cold stress for 24 h, and AtruCBF3 and AtruCBF5 were always down-regulated under cold conditions. Taken together, comprehensive analysis will enhance our understanding of the potential functions of the CBF genes on cold resistance, thereby providing a reference for the introduction of Acer species in our country.

Genome-wide characterization and functional identification of MYB genes in Malus sieversii infected by Valsa mali

Among the most important transcription factors in plants, the v-myb avian myeloblastosis viral oncogene homolog (MYB) regulates the expression network of response genes under stresses such as fungal infection. In China, the canker disease Valsa mali threatens the survival of Malus sieversii, an ancestor of cultivated apples. Using the M. sieversii genome, we identified 457 MsMYB and 128 R2R3-MsMYB genes that were randomly distributed across 17 chromosomes. Based on protein sequence and structure, the R2R3-MsMYB genes were phylogenetically divided into 29 categories, and 26 conserved motifs were identified. We further predicted cis-elements in the 2000-kb promoter region of R2R3-MsMYBs based on the genome. Transcriptome analysis of M. sieversii under V. mali infection showed that 27 R2R3-MsMYBs were significantly differentially expressed, indicating their key role in the response to V. mali infection. Using transient transformation, MsMYB14, MsMYB24, MsMYB39, MsMYB78, and MsMYB108, which were strongly induced by V. mali infection, were functionally identified. Among the five MsMYBs, MsMYB14 and MsMYB78 were both important in enhancing resistance to diseases, whereas MsMYB24 inhibited resistance. Based on the results of this study, we gained a better understanding of the MsMYB transcription factor family and laid the foundation for a future research program on disease prevention strategies in M. sieversii.

Genome-wide identification and expression analysis of MYB gene family under nitrogen stress in Panax notoginseng

The myeloblastosis (MYB) gene family, involved in regulating many important physiological and biochemical processes, is one of the largest transcript factor superfamilies in plants. Since the identification of genome sequencing of Panax notoginseng has been completed, there was little known about the whole genome of its specific MYB gene family and the response to abiotic stresses, in consideration of the excessive application of nitrogen fertilizers in P. notoginseng. In this study, 123 PnMYB genes (MYB genes of P. notoginseng) have been identified and divided into 3 subfamilies by the phylogenetic analysis. These PnMYB genes were unevenly located on 12 chromosomes. Meanwhile, the gene structure and protein conserved domain were established by MEME Suite. The analysis of collinear relationships reflected that there were 121 homologous genes between P. notoginseng and Arabidopsis and 30 between P. notoginseng and rice. Moreover, cis-acting elements of PnMYB gene promoters were predicted which indicated that PnMYBs are involved in biotic, abiotic stress, and hormone induction. The expressions of PnMYB transcription factors in its roots, flowers, and leaves were detected by qRT-PCR and they had tissue-specific expressions and related to the growth of different tissues. Under nitrogen stress, MYB transcription factors had great feedback. Ten R2R3-MYB subfamily genes were significantly induced and indicated the possible function of protecting P. notoginseng from excess nitrogen. With further knowledge on identification of PnMYB gene related to tissue selectivity and abiotic stresses, this study laid the foundation for the functional development of PnMYB gene family and improved the cultivation of P. notoginseng.

Deciphering the roles of tobacco MYB transcription factors in environmental stress tolerance

The MYB members play important roles in development, metabolism, and stress tolerance in plants. In the current study, a total of 246 tobacco R2R3-MYB transcription factors were identified and systemically analyzed from the latest genome annotation. The newly identified tobacco members were divided into 33 subgroups together with the Arabidopsis members. Furthermore, 44 NtMYB gene pairs were identified to arise from duplication events, which might lead to the expansion of tobacco MYB genes. The expression patterns were revealed by transcriptomic analysis. Notably, the results from phylogenetic analysis, synthetic analysis, and expression analysis were integrated to predict the potential functions of these members. Particularly, NtMYB102 was found to act as the homolog of AtMYB70 and significantly induced by drought and salt treatments. The further assays revealed that NtMYB102 had transcriptional activities, and the overexpression of the encoding gene enhanced the drought and salt stress tolerance in transgenic tobacco. The results of this study may be relevant for future functional analyses of the MYB genes in tobacco.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Cymbidium goeringii

The MYB gene family plays a vital regulatory role in plant metabolism, stress response, and floral color. The R2R3-MYB gene family of C. goeringii was identified, and its expression was analyzed using bioinformatics in this article. The R2R3-MYB genes of Arabidopsis thaliana were used as a reference to determine 104 CgMYB genes and categorize them into 22 subfamilies. Exon/intron organizations and conserved motif analysis revealed that the majority of CgMYB genes were highly conserved, and chromosome localization and collinearity analysis provided evidence of tandem duplication and segmental duplication events, indicating the phenomenon of gene family expansion and contraction. The function of CgMYB genes was analyzed by cis-acting element and gene ontology (GO) enrichment. In addition, we selected CgMYB91 and CgMYB32 for RT-qPCR, suggesting that CgMYB91 and CgMYB32 are associated with anthocyanin formation. In short, this study provides a comprehensive and specific function of the R2R3-MYB transcription factors (TFs) in orchids.

Systematic analysis of the pectin methylesterase gene family in Nicotiana tabacum and reveal their multiple roles in plant development and abiotic stresses

The pectin methylesterases (PMEs) play multiple roles in regulating plant development and responses to various stresses. In our study, a total of 121 PME genes were identified in the tobacco genome, which were clustered into two groups based on phylogenetic analysis together with Arabidopsis members. The investigations of gene structure and conserved motif indicated that exon/intron and motif organizations were relatively conserved in each group. Additionally, several stress-related elements were identified in the promoter region of these genes. The survey of duplication events revealed that segmental duplications were critical to the expansion of the PME gene family in tobacco. The expression profiles analysis revealed that these genes were expressed in various tissues and could be induced by diverse abiotic stresses. Notably, NtPME029 and NtPME043, were identified as homologues with AtPME3 and AtPME31, respectively. Furthermore, NtPME029 was highly expressed in roots and the over-expression of the NtPME029 gene could promote the development of roots. While NtPME043 could be induced by salt and ABA treatments, and the over-expression of the NtPME043 gene could significantly enhance the salt-stress tolerance in tobacco. Overall, these findings may shed light on the biological and functional characterization of NtPME genes in tobacco.

Genome-wide analysis of MYB family genes in Tripterygium wilfordii and their potential roles in terpenoid biosynthesis

Terpenoids are a class of significant bioactive components in the woody vine of Tripterygium wilfordii. Previous studies have shown that MYB transcription factors play important roles in plant secondary metabolism, growth, and developmental processes. However, the MYB involved in terpenoid biosynthesis in Tripterygium wilfordii are unknown. To identify Tripterygium wilfordii MYB (TwMYB) genes that are involved in terpenoid biosynthesis, we conducted the genome-wide analysis of the TwMYB gene family. A total of 207 TwMYBs were identified including 84 1R-TwMYB, 117 R2R3-TwMYB, four 3R-TwMYB, and two 4R-TwMYB genes. The most abundant R2R3-TwMYBs together with their Arabidopsis homologs were categorized into 26 subgroups. Intraspecific collinearity analysis found that the 74.9% of the TwMYBs may be generated by segmental duplication events, and 36.7% of duplicated gene pairs were derived from the specific whole genome duplication (WGD) event in Tripterygium wilfordii. In addition, interspecies collinearity analysis found that 16 TwMYB genes formed homologous gene pairs with MYB genes in seven representative species, which indicated they may have a key role in evolution. Notably, we found that the TwMYB genes were differentially expressed in various tissues by expression pattern analysis. In order to further select the candidate genes related to terpenoid biosynthesis, the assay of Methyl jasmonate (MeJA) induction and analysis of phylogenetic tree was conducted. It was speculated that six candidate TwMYB genes (TwMYB33, TwMYB34, TwMYB45, TwMYB67, TwMYB102, and TwMYB103) are involved in regulating terpenoid biosynthesis. This study is the first systematic analysis of the TwMYB gene family and will lay a foundation for the functional characterization of TwMYB genes in the regulation of terpenoid biosynthesis.

Systematic Analysis of Tobacco CrRLK1L Family Genes and Functional Identification of NtCrRLK1L47 in Environmental Stresses

The Catharanthus roseus RLK1-like (CrRLK1L) family is involved in the regulation of plant reproduction, growth and development, cell wall integrity sensing, as well as responses to both biotic and abiotic stress conditions. Extraordinary progress has been made in elucidating the CrRLK1L family receptor kinases-mediated signaling pathway, while limited research addressed the functions of CrRLK1L proteins in tobacco. In this study, we identified and analyzed 48 NtCrRLK1L members from the tobacco genome. The newly identified NtCrRLK1L members were divided into seven groups together with the Arabidopsis CrRLK1L members. The syntenic analysis revealed that four pairs of NtCrRLK1L genes were predicted to have arisen from segmental duplication events. Expression profiling showed that the NtCrRLK1L genes were expressed in various tissues, and most NtCrRLK1L genes were induced by salt and drought stress conditions. Notably, NtCrRLK1L47 was upregulated under drought and salinity stresses, and the NtCrRLK1L47-GFP fusion protein was located in the cell membrane. Furthermore, overexpression of the NtCrRLK1L47 gene enhanced the salt tolerance in tobacco seedlings.

Identification and Analysis of MYB Gene Family for Discovering Potential Regulators Responding to Abiotic Stresses in Curcuma wenyujin

MYB superfamily is one of the most abundant families in plants, and plays critical role in plant growth, development, metabolism regulation, and stress response. Curcuma wenyujin is the main source plant of three traditional Chinese medicines, which are widely used in clinical treatment due to its diverse pharmacological activities. In present study, 88 CwMYBs were identified and analyzed in C. wenyujin, including 43 MYB-related genes, 42 R2R3-MYB genes, two 3R-MYB genes, and one 4R-MYB gene. Forty-three MYB-related proteins were classified into several types based on conserved domains and specific motifs, including CCA1-like type, R-R type, Myb-CC type, GARP-like type, and TBR-like type. The analysis of motifs in MYB DBD and no-MYB regions revealed the relevance of protein structure and function. Comparative phylogeny analysis divided 42 R2R3-MYB proteins into 19 subgroups and provided a reference for understanding the functions of some CwMYBs based on orthologs of previously characterized MYBs. Expression profile analysis of CwMYB genes revealed the differentially expressed genes responding to various abiotic stresses. Four candidate MYB genes were identified by combining the results of phylogeny analysis and expression analysis. CwMYB10, CwMYB18, CwMYB39, and CwMYB41 were significantly induced by cold, NaCl, and MeJA stress treatments. CwMYB18 and CwMYB41 were proved as regulators with activity of transcriptional activation, whereas CwMYB39 and CwMYB10 were not. They may participate in the response to abiotic stresses through different mechanisms in C. wenyujin. This study was the first step toward understanding the CwMYB family and the response to abiotic stresses in C. wenyujin.

CstMYB14 links ROS signaling, apocarotenoid metabolism, and stress response in Crocus sativus L

Reactive oxygen species (ROS) behave as signaling molecules and induce biosynthesis of many secondary metabolites, including apocarotenoids, which play critical roles in stress tolerance through radical scavenging. However, the mechanism that regulates ROS responsive apocarotenoid metabolism and subsequent stress response is unknown. In this study, an R2R3-MYB transcription factor (CstMYB14) was identified from Crocus sativus L., which acts as a regulator of apocarotenoid biosynthesis. CstMYB14 expression increases in response to H₂ O₂ in a concentration and time-dependent manner. CstMYB14 localizes to the nucleus and acts as a transcriptional activator. Over-expression of CstMYB14 in Crocus stigmas enhanced apocarotenoid biosynthesis. Yeast-one-hybrid demonstrated binding of CstMYB14 to promoters of two apocarotenoid pathway genes (phytoene synthase and carotenoid cleavage dioxygenase 2). Nicotiana benthamiana plants overexpressing CstMYB14 showed better growth and higher stress tolerance than wild type plants. Higher antioxidant activity in CstMYB14-Ox plants indicated that stress tolerance might be due to ROS scavenging. These results establish a molecular link between ROS signaling, apocarotenoid metabolism and stress tolerance. Further, CstMYB14 is shown to act as a key regulator which modulates ROS responsive biosynthesis of apocarotenoids which in turn impart stress tolerance through ROS scavenging.

Calcium-dependent protein kinase 2 plays a positive role in the salt stress response in potato

StCDPK2 is an early player in the salt stress response in potato plants; its overexpression promoted ROS scavenging, chlorophyll stability, and the induction of stress-responsive genes conferring tolerance to salinity. The salinity of soils affects plant development and is responsible for great losses in crop yields. Calcium-dependent protein kinases (CDPKs) are sensor-transducers that decode Ca²⁺ signatures triggered by abiotic stimuli and translate them into physiological responses. Histochemical analyses of potato plants harboring StCDPK2 promoter fused to the reporter gene β-glucuronidase (Pro_StCDPK2:GUS) revealed that GUS activity was high in the leaf blade and veins, it was restricted to root tips and lateral root primordia, and was observed upon stolon swelling. Comparison with Pro_StCDPK1:GUS and Pro_StCDPK3:GUS plants revealed their differential activities in the plant tissues. Pro_StCDPK2:GUS plants exposed to high salt presented enhanced GUS activity in roots which correlated with the numerous stress-responsive sites predicted in its promoter sequence. Moreover, StCDPK2 expression increased in in vitro potato plants after 2 h of high salt exposure and in greenhouse plants exposed to a dynamic stress condition. As inferred from biometric data and chlorophyll content, plants that overexpress StCDPK2 were more tolerant than wild-type plants when exposed to high salt. Overexpressing plants have a more efficient antioxidant system; they showed reduced accumulation of peroxide and higher catalase activity under salt conditions, and enhanced expression of WRKY6 and ERF5 transcription factors under control conditions. Our results indicate that StCDPK2 is an early player in the salt stress response and support a positive correlation between StCDPK2 overexpression and tolerance towards salt stress.

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 genetic characterization of the CaMYB family and its response to five types of heavy metal stress in hot pepper (Capsicum annuum cv. CM334)

MYB proteins play a crucial role in plant growth and development and stress responses. In this study, 160 members of the MYB gene family from the pepper genome database were used to analyze gene structures, chromosome localization, collinearity, genetic affinity and expression in response to heavy metals. The results identified R2R3-MYB members and further phylogenetically classified them into 35 subgroups based on highly conserved gene structures and motifs. Collinearity analysis showed that segmental duplication events played a crucial role in the functional expansion of the CaMYB gene family by intraspecific collinearity, and at least 12 pairs of CaMYB genes existed between species prior to the differentiation between monocots and dicots. Moreover, the upstream CaMYB genes were mainly localized to the phytohormone elements ABRE and transcription factor elements MYB and MYC. Further analysis revealed that MYB transcription factors were closely associated with a variety of abiotic stress-related proteins (e.g., MAC-complex and SKIP). Under the stress of five metal ions, Cd²⁺, Cu²⁺, Pb²⁺, Zn²⁺, and Fe³⁺, the expression levels of some CaMYB family genes were upregulated. Of these genes, pairing homologous 1 (PH-1), PH-13, and PH-15 in the roots of Capsicum annuum were upregulated to the greatest extent, indicating that these three MYB family members are particularly sensitive to these five metals. This study provides a theoretical reference for the analysis of the molecular regulatory mechanism of MYB family genes in mediating the response to heavy metals in plants. This study reveals the mode of interaction between MYB and a variety of abiotic stress proteins and clarifies the biological functions of CaMYB family members in the regulation of heavy metal stress.

Genome-wide screening and identification of nuclear Factor-Y family genes and exploration their function on regulating abiotic and biotic stress in potato (Solanum tuberosum L.)

The Nuclear Factor-Y (NF-Y) transcription factor (TF), which includes three distinct subunits (NF-YA, NF-YB and NF-YC), is known to manipulate various aspects of plant growth, development, and stress responses. Although the NF-Y gene family was well studied in many species, little is known about their functions in potato. In this study, a total of 37 potato NF-Y genes were identified, including 11 StNF-YAs, 20 StNF-YBs, and 6 StNF-YCs. The genetic features of these StNF-Y genes were investigated by comparing their evolutionary relationship, intron/exon organization and motif distribution pattern. Multiple alignments showed that all StNF-Y proteins possessed clearly conserved core regions that were flanked by non-conserved sequences. Gene duplication analysis indicated that nine StNF-Y genes were subjected to tandem duplication and eight StNF-Ys arose from segmental duplication events. Synteny analysis suggested that most StNF-Y genes (33 of 37) were orthologous to potato's close relative tomato (Solanum lycopersicum L.). Tissue-specific expression of the StNF-Y genes suggested their potential roles in controlling potato growth and development. The role of StNF-Ys in regulating potato responses to abiotic stress (ABA, drought and salinity) was also confirmed: twelve StNF-Y genes were up-regulated and another two were down-regulated under different abiotic treatments. In addition, genes responded differently to pathogen challenges, suggesting that StNF-Y genes may play distinct roles under certain biotic stress. In summary, insights into the evolution of NF-Y family members and their functions in potato development and stress responses are provided.

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.

Transcriptome Profiles of Leaves and Roots of Goldenrain Tree (Koelreuteria paniculata Laxm.) in Response to Cadmium Stress

Cadmium (Cd) pollution is a widespread environmental problem. In this study, we explored the transcriptome and biochemical responses of goldenrain tree (Koelreuteria paniculata Laxm.) leaves and roots to Cd stress. Leaf and root growth decreased substantially under Cd stress (50 mg/L CdCl₂), but leaf and root antioxidant mechanisms were significantly activated. In RNA-seq analysis, roots treated with 25 mg/L CdCl₂ featured enriched GO terms in cellular components related to intracellular ribonucleoprotein complex, ribonucleoprotein complex, and macromolecular complex. In leaves under Cd stress, most differentially expressed genes were enriched in the cellular component terms intrinsic component of membrane and membrane part. Weighted gene co-expression network analysis and analysis of module–trait relations revealed candidate genes associated with superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities and malondialdehyde (MDA). Ten transcription factors responded to Cd stress expression, including those in C2H2, MYB, WRKY, and bZIP families. Transcriptomic analysis of goldenrain tree revealed that Cd stress rapidly induced the intracellular ribonucleoprotein complex in the roots and the intrinsic component of membrane in the leaves. The results also indicate directions for further analyses of molecular mechanisms of Cd tolerance and accumulation in goldenrain tree.

Genome-wide analysis and characterization of GRAS family in switchgrass

Panicum virgatum, a model plant of cellulosic ethanol conversion, not only has high large biomass and strong adaptability to soil, but also grows well in marginal soil and has the advantage of improving saline-alkali soil. GRAS transcription factor gene family play important roles in individual environment adaption, and these vital functions has been proved in several plants, however, the research of GRAS in the development of switchgrass (Panicum virgatum) were limited. A comprehensive study was investigated to explore the relationship between GRAS gene family and resistance. According to the phylogenetic analysis, a total of 144 GRAS genes were identified and renamed which were classified into eight subfamilies. Chromosome distribution, tandem and segmental repeats analysis indicated that gene duplication events contributed a lot to the expansion of GRAS genes in the switchgrass genome. Sixty-six GRAS genes in switchgrass were identified as having orthologous genes with rice through gene duplication analysis. Most of these GRAS genes contained zero or one intron, and closely related genes in evolution shared similar motif composition. Interaction networks were analyzed including DELLA and ten interaction proteins that were primarily involved in gibberellin acid mediated signaling. Notably, online analysis indicated that the promoter regions of the identified PvGRAS genes contained many cis-elements including light responsive elements, suggesting that PvGRAS might involve in light signal cross-talking. This work provides key insights into resistance and bioavailability in switchgrass and would be helpful to further study the function of GRAS and GRAS-mediated signal transduction pathways.

Recent Advanced Metabolic and Genetic Engineering of Phenylpropanoid Biosynthetic Pathways

The MYB transcription factors (TFs) are evolving as critical role in the regulation of the phenylpropanoid and tanshinones biosynthetic pathway. MYB TFs relate to a very important gene family, which are involved in the regulation of primary and secondary metabolisms, terpenoids, bioactive compounds, plant defense against various stresses and cell morphology. R2R3 MYB TFs contained a conserved N-terminal domain, but the domain at C-terminal sorts them different regarding their structures and functions. MYB TFs suppressors generally possess particular repressive motifs, such as pdLNLD/ELxiG/S and TLLLFR, which contribute to their suppression role through a diversity of complex regulatory mechanisms. A novel flower specific "NF/YWSV/MEDF/LW" conserved motif has a great potential to understand the mechanisms of flower development. In the current review, we summarize recent advanced progress of MYB TFs on transcription regulation, posttranscriptional, microRNA, conserved motif and propose directions to future prospective research. We further suggest there should be more focus on the investigation for the role of MYB TFs in microalgae, which has great potential for heterologous protein expression system for future perspectives.

Genome-Wide Identification, Phylogenetic and Expression Pattern Analysis of GATA Family Genes in Cucumber (Cucumis sativus L.)

GATA transcription factors are a class of transcriptional regulatory proteins that contain a characteristic type-IV zinc finger DNA-binding domain, which play important roles in plant growth and development. The GATA gene family has been characterized in various plant species. However, GATA family genes have not been identified in cucumber. In this study, 26 GATA family genes were identified in cucumber genome, whose physicochemical characteristics, chromosomal distributions, phylogenetic tree, gene structures conserved motifs, cis-regulatory elements in promoters, homologous gene pairs, downstream target genes were analyzed. Tissue expression profiles of cucumber GATA family genes exhibited that 17 GATA genes showed constitutive expression, and some GATA genes showed tissue-specific expression patterns. RNA-seq analysis of green and virescent leaves revealed that seven GATA genes might be involved in the chloroplast development and chlorophyll biosynthesis. Importantly, expression patterns analysis of GATA genes in response to abiotic and biotic stresses indicated that some GATA genes respond to either abiotic stress or biotic stress, some GATA genes such as Csa2G162660, Csa3G017200, Csa3G165640, Csa4G646060, Csa5G622830 and Csa6G312540 were simultaneously functional in resistance to abiotic and biotic stresses. Overall, this study will provide useful information for further analysis of the biological functions of GATA factors in cucumber.

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

BACKGROUND

MYB transcription factors are a kind of DNA binding protein that can specifically interact with the promoter region. Members of MYB TFs are widely involved in plant growth and development, secondary metabolism, stress response, and hormone signal transduction. However, there is no report of comprehensive bioinformatics analysis on the MYB family of Casuarina equisetifolia.

RESULTS

In this study, bioinformatics methods were used to screen out 182 MYB transcription factors from the Casuarina equisetifolia genome database, including 69 1R-MYB, 107 R2R3-MYB, 4 R1R2R3-MYB, and 2 4R-MYB. The C. equisetifolia R2R3-MYB genes were divided into 29 groups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis thaliana, while the remaining MYB genes (1R-MYB, R1R2R3-MYB, and 4R-MYB) was divided into 19 groups. Moreover, the conserved motif and gene structure analysis shown that the members of the CeqMYBs were divided into the same subgroups with mostly similar gene structures. In addition, many conserved amino acids in the R2 and R3 domains of CeqMYBs by WebLogo analysis, especially tryptophan residues (W), with 3 conserved W in R2 repeat and 2 conserved W in R3 repeat. Combining promoter and GO annotation analysis, speculated on the various biological functions of CeqMYBs, thus 32 MYB genes were selected to further explore its response to salt stress by using qPCR analysis technique. Most CeqMYB genes were differentially regulated following multiple salt treatments.

CONCLUSIONS

Seven genes (CeqMYB164, CeqMYB4, CeqMYB53, CeqMYB32, CeqMYB114, CeqMYB71 and CeqMYB177) were assigned to the "response to salt stress" by GO annotation. Among them, the expression level of CeqMYB4 was up-regulated under various salt treatments, indicating CeqMYB4 might participated in the response to salt stress. Our results provide important information for the biological function of C. equisetifolia, as well as offer candidate genes for further study of salt stress mechanism.

Genome-Wide Analysis of Myeloblastosis-Related Genes in Brassica napus L. and Positive Modulation of Osmotic Tolerance by BnMRD107

Myeloblastosis (MYB)-related transcription factors comprise a large subfamily of the MYB family. They play significant roles in plant development and in stress responses. However, MYB-related proteins have not been comprehensively investigated in rapeseed (Brassica napus L.). In the present study, a genome-wide analysis of MYB-related transcription factors was performed in rapeseed. We identified 251 Brassica napus MYB (BnMYB)-related members, which were divided phylogenetically into five clades. Evolutionary analysis suggested that whole genome duplication and segmental duplication events have played a significant role in the expansion of BnMYB-related gene family. Selective pressure of BnMYB-related genes was estimated using the Ka/Ks ratio, which indicated that BnMYB-related genes underwent strong purifying selection during evolution. In silico analysis showed that various development-associated, phytohormone-responsive, and stress-related cis-acting regulatory elements were enriched in the promoter regions of BnMYB-related genes. Furthermore, MYB-related genes with tissue or organ-specific, stress-responsive expression patterns were identified in B. napus based on temporospatial and abiotic stress expression profiles. Among the stress-responsive MYB-related genes, BnMRD107 was strongly induced by drought stress, and was therefore selected for functional study. Rapeseed seedlings overexpressing BnMRD107 showed improved resistance to osmotic stress. Our findings not only lay a foundation for further functional characterization of BnMYB-related genes, but also provide valuable clues to determine candidate genes for future genetic improvement of B. napus.

Wolfberry genomes and the evolution of Lycium (Solanaceae)

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

The TAZ domain-containing proteins play important role in the heavy metals stress biology in plants

TAZ (transcriptional coactivator with PDZ-binding) zinc finger domains, also known as transcription adaptor putative zinc finger domains, that control diverse function in plant growth and development. Here, in the present study, we evaluated the role of the TAZ domain-containing gene in response to various heavy metals. Initially, we found a total of 3, 7, 8, 9, 9, 9, 7, 14, 6, 10, and 6 proteins containing TAZ domain in stiff brome, millet, sorghum, potato, pepper, maize, rice, apple, peach, pear, and tomato genome that could trigger the plant resistance against various heavy metals, respectively. Various in-silico approaches were applied such as duplication, phylogenetic analysis, and gene structure, to understand the basic features of the TAZ domain-containing genes in plants. Gene expression analyses were also performed under heavy metals (Cr, Zn, Ni, Cd, Co, Fe, Mn, and Pb). The results of quantitative real-time PCR analysis indicated that the TAZ gene family members were differentially expressed under different heavy metals. We further characterized the functions of the TAZ domain-containing gene under the heavy metal stresses by overexpressing the OsTAZ4 gene in Arabidopsis. The TAZ genes could promote plant resistance against various heavy metals by interacting with OsMYB34 and OsFHA9 transcription factors. The results will contribute to elucidate the relationship of TAZ proteins with heavy metals stresses and also ascertain the biological function in plant growth and development.

Identification and characterization of 5 walnut MYB genes in response to drought stress involved in ABA signaling

Walnut is a popular nut tree species and usually suffers from drought stress. However, little information is available on the mechanism of walnut responding to drought stress, resulting in lack of basic understanding for its resistance. In order to excavate more functional genes that can respond to stressors, and enrich the theoretical basis for walnut resistance, in this study, 5 MYB genes with complete ORFs were identified from J. regia and the basic bio-information as well as expression patterns in different tissues and response to drought and ABA stresses were confirmed using qRT-PCR assay. The results showed that 2 JrMYB genes belong to R1-MYB subfamily and 3 JrMYBs belong to R2R3-MYB, encoding the proteins from 212 to 362 aa in length. The phylogenetic analysis categorized proteins of 5 JrMYBs and 40 Arabidopsis AtMYBs into 10 subgroups. JrMYBs in the same subgroup exhibited significant similarities in the composition of conserved domains and motifs in amino acid sequences and exon/intron organization in DNA sequences. The results of qRT-PCR analysis revealed that JrMYB genes diversely expressed in various tissues. Moreover, the expression values of JrMYBs were upregulated or downregulated significantly under drought and ABA stresses. Most attractively, in contrast with suffering from drought stress alone, the treatments with drought and additional ABA greatly enhanced the transcript levels of JrMYBs. All these results suggested that JrMYB genes play a vital role in plant biological processes and drought as well as ABA stress response, and possibly perform as ABA-dependent drought response transcription factors in plant.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01008-z.

Genome-Wide Identification, Classification and Expression Analysis of the MYB Transcription Factor Family in Petunia

A lot of researches have been focused on the evolution and function of MYB transcription factors (TFs). For revealing the formation of petunia flower color diversity, MYB gene family in petunia was identified and analyzed. In this study, a total of 155 MYB genes, including 40 1R-MYBs, 106 R2R3-MYBs, 7 R1R2R3-MYBs and 2 4R-MYBs, have been identified in the Petunia axillaris genome. Most R2R3 genes contain three exons and two introns, whereas the number of PaMYB introns varies from 0 to 12. The R2R3-MYB members could be divided into 28 subgroups. Analysis of gene structure and protein motifs revealed that members within the same subgroup presented similar exon/intron and motif organization, further supporting the results of phylogenetic analysis. Genes in subgroup 10, 11 and 21 were mainly expressed in petal, not in vegetative tissues. Genes in subgroup 9, 19, 25 and 27 expressed in all tissues, but the expression patterns of each gene were different. According to the promoter analysis, five R2R3-MYB and two MYB-related genes contained MBSI cis-element, which was involved in flavonoid biosynthetic regulation. PaMYB100/DPL has been reported to positively regulate to pigmentation. However, although PaMYB82, PaMYB68 and Pa1RMYB36 contained MBSI cis-element, their function in flavonoid biosynthesis has not been revealed. Consistent with existing knowledge, PaMYBs in subgroup 11 had similar function to AtMYBs in subgroup 6, genes in which played an important role in anthocyanin biosynthesis. In addition, PaMYB1 and PaMYB40 belonged to P9 (S7) and were potentially involved in regulation of flavonoid synthesis in petunia vegetative organs. This work provides a comprehensive understanding of the MYB gene family in petunia and lays a significant foundation for future studies on the function and evolution of MYB genes in petunia.

Bioinformatics and Expression Analysis of IDA-Like Genes Reveal Their Potential Functions in Flower Abscission and Stress Response in Tobacco (Nicotiana tabacum L.)

The inflorescence deficient in abscission-like (IDL) genes have been shown to play critical roles in floral organ abscission, lateral root formation and various stress responses in Arabidopsis. The IDL gene family has been characterized in a number of plant species, while limited information is available about IDL genes of tobacco. In the current study, 15 NtIDL members were identified in the tobacco genome, and were classified into six groups together with IDL members from other species. Evolution analysis suggested that the NtIDL members form group VI might have originated from duplication events. Notably, NtIDL06 shared high similarities with AtIDA in the EPIP sequence, and its encoding gene was highly expressed in the abscission zone of flowers at late developmental stages, implying that NtIDL06 might regulate tobacco flower abscission. In addition, the results from cis-elements analysis of promoters and expression after stress treatments suggested that NtIDL members might be involved in various stress responses of tobacco. The results from this study provide information for further functional analysis related to flower abscission and stress responses of NtIDL genes.

Potato NAC Transcription Factor StNAC053 Enhances Salt and Drought Tolerance in Transgenic Arabidopsis

The NAC (NAM, ATAF1/2, and CUC2) transcription factors comprise one of the largest transcription factor families in plants and play important roles in stress responses. However, little is known about the functions of potato NAC family members. Here we report the cloning of a potato NAC transcription factor gene StNAC053, which was significantly upregulated after salt, drought, and abscisic acid treatments. Furthermore, the StNAC053-GFP fusion protein was found to be located in the nucleus and had a C-terminal transactivation domain, implying that StNAC053 may function as a transcriptional activator in potato. Notably, Arabidopsis plants overexpressing StNAC053 displayed lower seed germination rates compared to wild-type under exogenous ABA treatment. In addition, the StNAC053 overexpression Arabidopsis lines displayed significantly increased tolerance to salt and drought stress treatments. Moreover, the StNAC053-OE lines were found to have higher activities of superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) under multiple stress treatments. Interestingly, the expression levels of several stress-related genes including COR15A,DREB1A, ERD11, RAB18, ERF5, and KAT2, were significantly upregulated in these StNAC053-overexpressing lines. Taken together, overexpression of the stress-inducible StNAC053 gene could enhance the tolerances to both salt and drought stress treatments in Arabidopsis, likely by upregulating stress-related genes.

Genome-wide identification and function characterization of GATA transcription factors during development and in response to abiotic stresses and hormone treatments in pepper

Pepper (Capsicum annuum L.) is an economically important vegetable crop whose production and quality are severely reduced under adverse environmental stress conditions. The GATA transcription factors belonging to type IV zinc-finger proteins, play a significant role in regulating light morphogenesis, nitrate assimilation, and organ development in plants. However, the functional characteristics of GATA gene family during development and in response to environmental stresses have not yet been investigated in pepper. In this study, a total of 28 pepper GATA (CaGATA) genes were identified. To gain an overview of the CaGATAs, we analyzed their chromosomal distribution, gene structure, conservative domains, cis-elements, phylogeny, and evolutionary relationship. We divided 28 CaGATAs into four groups distributed on 10 chromosomes, and identified 7 paralogs in CaGATA family of pepper and 35 orthologous gene pairs between CaGATAs and Arabidopsis GATAs (AtGATAs). The results of promoter cis-element analysis and the quantitative real-time PCR (qRT-PCR) analysis revealed that CaGATA genes were involved in regulating the plant growth and development and the responses to various abiotic stresses and hormone treatments in pepper. Tissue-specific expression analysis showed that most CaGATA genes were preferentially expressed in flower buds, flowers, and leaves. Several CaGATA genes, especially CaGATA14, were significantly regulated under multiple abiotic stresses, and CaGATA21 and CaGATA27 were highly responsive to phytohormone treatments. Taken together, our results lay a foundation for the biological function analysis of GATA gene family in pepper.

R2R3-MYB transcription factors, StmiR858 and sucrose mediate potato flavonol biosynthesis

Flavonols and other phenylpropanoids protect plants from biotic and abiotic stress and are dietarily desirable because of their health-promoting properties. The ability to develop new potatoes (Solanum tuberosum) with optimal types and amounts of phenylpropanoids is limited by lack of knowledge about the regulatory mechanisms. Exogenous sucrose increased flavonols, whereas overexpression of the MYB StAN1 induced sucrolytic gene expression. Heterologous StAN1 protein bound promoter fragments from sucrolytic genes (SUSY1 and INV1). Two additional MYBs and one microRNA were identified that regulated potato flavonols. Overexpression analysis showed MYB12A and C increased amounts of flavonols and other phenylpropanoids. Endogenous flavonol amounts in light-exposed organs were much higher those in the dark. Expression levels of StMYB12A and C were high in flowers but low in tubers. Transient overexpression of miR858 altered potato flavonol metabolism. Endogenous StmiR858 expression was much lower in flowers than leaves and correlated with flavonol amounts in these organs. Collectively, these findings support the hypothesis that sucrose, MYBs, and miRNA control potato phenylpropanoid metabolism in a finely tuned manner that includes a feedback loop between sucrose and StAN1. These findings will aid in the development of potatoes with phenylpropanoid profiles optimized for crop performance and human health.

Genome-Wide Comparative Analysis of R2R3 MYB Gene Family in Populus and Salix and Identification of Male Flower Bud Development-Related Genes

The MYB transcription factor (TF) family is one of the largest plant transcription factor gene family playing vital roles in plant growth and development, including defense, cell differentiation, secondary metabolism, and responses to biotic and abiotic stresses. As a model tree species of woody plants, in recent years, the identification and functional prediction of certain MYB family members in the poplar genome have been reported. However, to date, the characterization of the gene family in the genome of the poplar's sister species willow has not been done, nor are the differences and similarities between the poplar and willow genomes understood. In this study, we conducted the first genome-wide investigation of the R2R3 MYB subfamily in the willow, identifying 216 R2R3 MYB gene members, and combined with the poplar R2R3 MYB genes, performed the first comparative analysis of R2R3 MYB genes between the poplar and willow. We identified 81 and 86 pairs of R2R3 MYB paralogs in the poplar and willow, respectively. There were 17 pairs of tandem repeat genes in the willow, indicating active duplication of willow R2R3 MYB genes. A further 166 pairs of poplar and willow orthologs were identified by collinear and synonymous analysis. The findings support the duplication of R2R3 MYB genes in the ancestral species, with most of the R2R3 MYB genes being retained during the evolutionary process. The phylogenetic trees of the R2R3 MYB genes of 10 different species were drawn. The functions of the poplar and willow R2R3 MYB genes were predicted using reported functional groupings and clustering by OrthoFinder. Identified 5 subgroups in general expanded in woody species, three subgroups were predicted to be related to lignin synthesis, and we further speculate that the other two subgroups also play a role in wood formation. We analyzed the expression patterns of the GAMYB gene of subgroup 18 (S18) related to pollen development in the male flower buds of poplar and willow at different developmental stages by qRT-PCR. The results showed that the GAMYB gene was specifically expressed in the male flower bud from pollen formation to maturity, and that the expression first increased and then decreased. Both the specificity of tissue expression specificity and conservation indicated that GAMYB played an important role in pollen development in both poplar and willow and was an ideal candidate gene for the analysis of male flower development-related functions of the two species.

MYB Superfamily in Brassica napus: Evidence for Hormone-Mediated Expression Profiles, Large Expansion, and Functions in Root Hair Development

MYB proteins are involved in diverse important biological processes in plants. Herein, we obtained the MYB superfamily from the allotetraploid Brassica napus, which contains 227 MYB-related (BnMYBR/Bn1R-MYB), 429 R2R3-MYB (Bn2R-MYB), 22 R1R2R3-MYB (Bn3R-MYB), and two R1R2R2R1/2-MYB (Bn4R-MYB) genes. Phylogenetic analysis classified the Bn2R-MYBs into 43 subfamilies, and the BnMYBRs into five subfamilies. Sequence characteristics and exon/intron structures within each subfamily of the Bn2R-MYBs and BnMYBRs were highly conserved. The whole superfamily was unevenly distributed on 19 chromosomes and underwent unbalanced expansion in B. napus. Allopolyploidy between B. oleracea and B. rapa mainly contributed to the expansion in their descendent B. napus, in which B. rapa-derived genes were more retained. Comparative phylogenetic analysis of 2R-MYB proteins from nine Brassicaceae and seven non-Brassicaceae species identified five Brassicaceae-specific subfamilies and five subfamilies that are lacking from the examined Brassicaceae species, which provided an example for the adaptive evolution of the 2R-MYB gene family alongside angiosperm diversification. Ectopic expression of four Bn2R-MYBs under the control of the viral CaMV35S and/or native promoters could rescue the lesser root hair phenotype of the Arabidopsis thaliana wer mutant plants, proving the conserved negative roles of the 2R-MYBs of the S15 subfamily in root hair development. RNA-sequencing data revealed that the Bn2R-MYBs and BnMYBRs had diverse transcript profiles in roots in response to the treatments with various hormones. Our findings provide valuable information for further functional characterizations of B. napusMYB genes.

Abscisic Acid and Jasmonate Metabolisms Are Jointly Regulated During Senescence in Roots and Leaves of Populus trichocarpa

Plant senescence is a highly regulated process that allows nutrients to be mobilized from dying tissues to other organs. Despite that senescence has been extensively studied in leaves, the senescence of ephemeral organs located underground is still poorly understood, especially in the context of phytohormone engagement. The present study focused on filling this knowledge gap by examining the roles of abscisic acid (ABA) and jasmonate in the regulation of senescence of fine, absorptive roots and leaves of Populus trichocarpa. Immunohistochemical (IHC), chromatographic, and molecular methods were utilized to achieve this objective. A transcriptomic analysis identified significant changes in gene expression that were associated with the metabolism and signal transduction of phytohormones, especially ABA and jasmonate. The increased level of these phytohormones during senescence was detected in both organs and was confirmed by IHC. Based on the obtained data, we suggest that phytohormonal regulation of senescence in roots and leaves is organ-specific. We have shown that the regulation of ABA and JA metabolism is tightly regulated during senescence processes in both leaves and roots. The results were discussed with respect to the role of ABA in cold tolerance and the role of JA in resistance to pathogens.

Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

Genome-wide analysis and expression profiles of the StR2R3-MYB transcription factor superfamily in potato (Solanum tuberosum L.)

MYB transcription factors comprise one of the largest families in plant kingdom, which play a variety of functions in plant developmental processes and defence responses, the R2R3-MYB members are the predominant form found in higher plants. In the present study, a total of 111 StR2R3-MYB transcription factors were identified and further phylogenetically classified into 31 subfamilies, as supported by highly conserved gene structures and motifs. Collinearity analysis showed that the segmental duplication events played a crucial role in the expansion of StR2R3-MYB gene family. Synteny analysis indicated that 37 and 13 StR2R3-MYB genes were orthologous to Arabidopsis and wheat (Triticum aestivum), respectively, and these gene pairs have evolved under strong purifying selection. RNA-seq data from different tissues and abiotic stresses revealed tissue-preferential and abiotic stress-responsive StR2R3-MYB genes. We further analyzed StR2R3-MYB genes might be involved in anthocyanin biosynthesis and drought stress by using RNA-seq data of pigmented tetraploid potato cultivars and drought-sensitive and -tolerant tetraploid potato cultivars under drought stress, respectively. Moreover, EAR motifs were found in 21 StR2R3-MYB proteins and 446 pairs of proteins were predicted to interact with 21 EAR motif-containing StR2R3-MYB proteins by constructing the interaction network with medium confidence (0.4). Additionally, Gene Ontology (GO) analysis of the 21 EAR motif-containing StR2R3-MYB proteins was performed to further investigate their functions. This work will facilitate future biologically functional studies of potato StR2R3-MYB transcription factors and enrich the knowledge of MYB superfamily genes in plant species.

Changes in Ecophysiology, Osmolytes, and Secondary Metabolites of the Medicinal Plants of Mentha piperita and Catharanthus roseus Subjected to Drought and Heat Stress

Global warming contributes to higher temperatures and reduces rainfall for most areas worldwide. The concurrent incidence of extreme temperature and water shortage lead to temperature stress damage in plants. Seeking to imitate a more natural field situation and to figure out responses of specific stresses with regard to their combination, we investigated physiological, biochemical, and metabolomic variations following drought and heat stress imposition (alone and combined) and recovery, using Mentha piperita and Catharanthus roseus plants. Plants were exposed to drought and/or heat stress (35 °C) for seven and fourteen days. Plant height and weight (both fresh and dry weight) were significantly decreased by stress, and the effects more pronounced with a combined heat and drought treatment. Drought and/or heat stress triggered the accumulation of osmolytes (proline, sugars, glycine betaine, and sugar alcohols including inositol and mannitol), with maximum accumulation in response to the combined stress. Total phenol, flavonoid, and saponin contents decreased in response to drought and/or heat stress at seven and fourteen days; however, levels of other secondary metabolites, including tannins, terpenoids, and alkaloids, increased under stress in both plants, with maximal accumulation under the combined heat/drought stress. Extracts from leaves of both species significantly inhibited the growth of pathogenic fungi and bacteria, as well as two human cancer cell lines. Drought and heat stress significantly reduced the antimicrobial and anticancer activities of plants. The increased accumulation of secondary metabolites observed in response to drought and/or heat stress suggests that imposition of abiotic stress may be a strategy for increasing the content of the therapeutic secondary metabolites associated with these plants.

MYB Gene Family in Potato (Solanum tuberosum L.): Genome-Wide Identification of Hormone-Responsive Reveals Their Potential Functions in Growth and Development

As an important nongrain crop, the growth and yield of potato (Solanum tuberosum L.) is often affected by an unfavorable external environment in the process of cultivation. The MYB family is one of the largest and most important gene families, participating in the regulation of plant growth and development and response to abiotic stresses. Several MYB genes in potato that regulate anthocyanin synthesis and participate in abiotic stress responses have been identified. To identify all Solanum tuberosum L. MYB (StMYB) genes involved in hormone or stress responses to potentially regulate potato growth and development, we identified the MYB gene family at the genome-wide level. In this work, 158 StMYB genes were found in the potato genome. According to the amino acid sequence of the MYB domain and gene structure, the StMYB genes were divided into R2R3-MYB and R1R2R3-MYB families, and the R2R3-MYB family was divided into 20 subgroups (SGs). The expression of 21 StMYB genes from different SGs in roots, stems, leaves, flowers, shoots, stolons, young tubers, and mature tubers was determined by quantitative real-time polymerase chain reaction (qRT-PCR). The expression patterns of StMYB genes in potatoes treated with abscisic acid (ABA), indole-3-acetic acid (IAA), gibberellin acid 3 (GA3), NaCl, mannitol, and heat were also measured. We have identified several potential candidate genes that regulate the synthesis of potato flavonoids or participate in hormone or stress responses. This work provides a comprehensive understanding of the MYB family in potato and will lay a foundation for the future investigation of the potential functions of StMYB genes in the growth and development of potato.