Genome-Wide Identification and Expression Analysis of MYB Transcription Factor Superfamily in Dendrobium catenatum

Regulatory Dynamics of Plant Hormones and Transcription Factors under Salt Stress

The negative impacts of soil salinization on ion homeostasis provide a significant global barrier to agricultural production and development. Plant physiology and biochemistry are severely affected by primary and secondary NaCl stress impacts, which damage cellular integrity, impair water uptake, and trigger physiological drought. Determining how transcriptional factors (TFs) and hormone networks are regulated in plants in response to salt stress is necessary for developing crops that tolerate salt. This study investigates the complex mechanisms of several significant TF families that influence plant responses to salt stress, involving AP2/ERF, bZIP, NAC, MYB, and WRKY. It demonstrates how these transcription factors (TFs) help plants respond to the detrimental effects of salinity by modulating gene expression through mechanisms including hormone signaling, osmotic stress pathway activation, and ion homeostasis. Additionally, it explores the hormonal imbalances triggered by salt stress, which entail complex interactions among phytohormones like jasmonic acid (JA), salicylic acid (SA), and abscisic acid (ABA) within the hormonal regulatory networks. This review highlights the regulatory role of key transcription factors in salt-stress response, and their interaction with plant hormones is crucial for developing genome-edited crops that can enhance agricultural sustainability and address global food security challenges.

A Sweet Potato MYB Transcription Factor IbMYB330 Enhances Tolerance to Drought and Salt Stress in Transgenic Tobacco

MYB transcription factors (TFs) play vital roles in plant growth, development, and response to adversity. Although the MYB gene family has been studied in many plant species, there is still little known about the function of R2R3 MYB TFs in sweet potato in response to abiotic stresses. In this study, an R2R3 MYB gene, IbMYB330 was isolated from sweet potato (Ipomoea batatas). IbMYB330 was ectopically expressed in tobacco and the functional characterization was performed by overexpression in transgenic plants. The IbMYB330 protein has a 268 amino acid sequence and contains two highly conserved MYB domains. The molecular weight and isoelectric point of IbMYB330 are 29.24 kD and 9.12, respectively. The expression of IbMYB330 in sweet potato is tissue-specific, and levels in the root were significantly higher than that in the leaf and stem. It showed that the expression of IbMYB330 was strongly induced by PEG-6000, NaCl, and H2O2. Ectopic expression of IbMYB330 led to increased transcript levels of stress-related genes such as SOD, POD, APX, and P5CS. Moreover, compared to the wild-type (WT), transgenic tobacco overexpression of IbMYB330 enhanced the tolerance to drought and salt stress treatment as CAT activity, POD activity, proline content, and protein content in transgenic tobacco had increased, while MDA content had decreased. Taken together, our study demonstrated that IbMYB330 plays a role in enhancing the resistance of sweet potato to stresses. These findings lay the groundwork for future research on the R2R3-MYB genes of sweet potato and indicates that IbMYB330 may be a candidate gene for improving abiotic stress tolerance in crops.

Comparative transcriptional analysis of Persea americana MYB, WRKY and AP2/ERF transcription factors following Phytophthora cinnamomi infection

Plant cells undergo extensive transcriptional reprogramming following pathogen infection, with these reprogramming patterns becoming more complex when pathogens, such as hemibiotrophs, exhibit different lifestyles. These transcriptional changes are often orchestrated by MYB, WRKY and AP2/ERF transcription factors (TFs), which modulate both growth and defence-related gene expression. Transcriptional analysis of defence-related genes in avocado (Persea americana) infected with Phytophthora cinnamomi indicated differential immune response activation when comparing a partially resistant and susceptible rootstock. This study identified 226 MYB, 82 WRKY, and 174 AP2/ERF TF-encoding genes in avocado, using a genome-wide approach. Phylogenetic analysis revealed substantial sequence conservation within TF groups underscoring their functional significance. RNA-sequencing analysis in a partially resistant and susceptible avocado rootstock infected with P. cinnamomi was indicative of an immune response switch occurring in either rootstock after 24 and 6 h post-inoculation, respectively. Different clusters of co-expressed TF genes were observed at these times, suggesting the activation of necrotroph-related immune responses at varying intervals between the two rootstocks. This study aids our understanding of avocado immune response activation following P. cinnamomi infection, and the role of the TFs therein, elucidating the transcriptional reprogramming disparities between partially resistant and susceptible rootstocks.

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.

Whole-Transcriptome Sequencing Reveals the Global Molecular Responses and NAC Transcription Factors Involved in Drought Stress in Dendrobium catenatum

Dendrobium catenatum is a highly drought-tolerant herb, which usually grows on cliffs or in the branches of trees, yet the underlying molecular mechanisms for its tolerance remain poorly understood. We conducted a comprehensive study utilizing whole-transcriptome sequencing approaches to investigate the molecular response to extreme drought stress in D. catenatum. A large number of differentially expressed mRNAs, lncRNAs, and circRNAs have been identified, and the NAC transcription factor family was highly enriched. Meanwhile, 46 genes were significantly up-regulated in the ABA-activated signaling pathway. In addition to the 89 NAC family members accurately identified in this study, 32 members were found to have different expressions between the CK and extreme drought treatment. They may regulate drought stress through both ABA-dependent and ABA-independent pathways. Moreover, the 32 analyzed differentially expressed DcNACs were found to be predominantly expressed in the floral organs and roots. The ceRNA regulatory network showed that DcNAC87 is at the core of the ceRNA network and is regulated by miR169, miR393, and four lncRNAs. These investigations provided valuable information on the role of NAC transcription factors in D. catenatum's response to drought stress.

Integrated metabolome, full-length sequencing, and transcriptome analyses unveil the molecular mechanisms of color formation of the canary yellow and red bracts of Bougainvillea × buttiana 'Chitra'

Bougainvillea is a typical tropical flower of great ornamental value due to its colorful bracts. The molecular mechanism behind color formation is not well-understood. Therefore, this research conducted metabolome analysis, transcriptome analysis, and multi-flux full-length sequencing in two color bracts of Bougainvillea × buttiana 'Chitra' to investigate the significantly different metabolites (SDMs) and differentially expressed genes (DEGs). Overall, 261 SDMs, including 62 flavonoids and 26 alkaloids, were detected, and flavonols and betalains were significantly differentially accumulated among the two bracts. Furthermore, the complete-length transcriptome of Bougainvillea × buttiana was also developed, which contained 512 493 non-redundant isoforms. Among them, 341 210 (66.58%) displayed multiple annotations in the KOG, GO, NR, KEGG, Pfam, Swissprot, and NT databases. RNA-seq findings revealed that 3610 DEGs were identified between two bracts. Co-expression analysis demonstrated that the DEGs and SDMs involved in flavonol metabolism (such as CHS, CHI, F3H, FLS, CYP75B1, kaempferol, and quercetin) and betacyanin metabolism (DODA, betanidin, and betacyanins) were the main contributors for the canary yellow and red bract formation, respectively. Further investigation revealed that several putative transcription factors (TFs) might interact with the promoters of the genes mentioned above. The expression profiles of the putative TFs displayed that they may positively and negatively regulate the structural genes' expression profiles. The data revealed a potential regulatory network between important genes, putative TFs, and metabolites in the flavonol and betacyanin biosynthesis of Bougainvillea × buttiana 'Chitra' bracts. These findings will serve as a rich genetic resource for future studies that could create new color bracts.

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.

Genome-Wide Characterization and Expression Analysis of Transcription Factor Families in Desert Moss Syntrichia caninervis under Abiotic Stresses

Transcription factor (TF) families play important roles in plant stress responses. S. caninervis is a new model moss for plant desiccation tolerance studies. Here, we report a high-confidence identification and characterization of 591 TFs representing 52 families that covered all chromosomes in S. caninervis. GO term and KEGG pathway analysis showed that TFs were involved in the regulation of transcription, DNA-templated, gene expression, binding activities, plant hormone signal transduction, and circadian rhythm. A number of TF promoter regions have a mixture of various hormones-related cis-regulatory elements. AP2/ERF, bHLH, MYB, and C2H2-zinc finger TFs were the overrepresented TF families in S. caninervis, and the detailed classification of each family is performed based on structural features. Transcriptome analysis revealed the transcript abundances of some ScAP2/ERF, bHLH, MYB, and C2H2 genes were accumulated in the treated S. caninervis under cold, dehydration, and rehydration stresses. The RT-qPCR results strongly agreed with RNA-seq analysis, indicating these TFs might play a key role in S. caninervis response to abiotic stress. Our comparative TF characterization and classification provide the foundations for functional investigations of the dominant TF genes involved in S. caninervis stress response, as well as excellent stress tolerance gene resources for plant stress resistance breeding.

Genome-Scale Identification, Classification, and Expression Profiling of MYB Transcription Factor Genes in Cinnamomum camphora

The camphor tree (Cinnamomum camphora (L.) Presl.) is the representative species of subtropical evergreen broadleaved forests in eastern Asia and an important raw material for essential oil production worldwide. Although MYBs have been comprehensively characterized and their functions have been partially resolved in many plants, it has not been explored in C. camphora. In this study, 121 CcMYBs were identified on 12 chromosomes in the whole genome of C. camphora and found that CcMYBs were mainly expanded by segmental duplication. They were divided into 28 subgroups based on phylogenetic analysis and gene structural characteristics. In the promoter regions, numerous cis-acting elements were related to biological processes. Analysis of RNA sequencing data from seven tissues showed that CcMYBs exhibited different expression profiles, suggesting that they have various roles in camphor tree development. In addition, combined with the correlation analysis of structural genes in the flavonoid synthesis pathway, we identified CcMYBs from three subgroups that might be related to the flavonoid biosynthesis pathway. This study systematically analyzed CcMYBs in C. camphora, which will set the stage for subsequent research on the functions of CcMYBs during their lifetime and provide valuable insights for the genetic improvement of camphor trees.

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

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

DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium hybrid

BACKGROUND

Flower pigment and shape are determined by the coordinated expression of a set of structural genes during flower development. R2R3-MYB transcription factors are known regulators of structural gene expression. The current study focused on two members of this large family of transcription factors that were predicted to have roles in pigment biosynthesis and organ shape development in orchids.

METHODS

Phylogenetic analysis was used to identify candidate Dendrobium catenatum R2R3-MYB (DcaMYB) sequences associated with pigment and cell shape development. Gene silencing of candidate DhMYBs in Dendrobium hybrid by direct application of dsRNA to developing flowers was followed by observation of gene expression level and flower phenotypes. Silencing of the structural gene chalcone synthase was used as a comparative control.

KEY RESULTS

Ten candidate flower-associated DcaMYBs were identified. Flowers treated with dsRNA of DhMYB22 and DhMYB60 sequences were less pigmented and had relatively low expression of anthocyanin biosynthetic genes (F3'H and DFR), lower total anthocyanin concentration and markedly lower levels of cyanidin-3-glucoside and cyanidin-3-rutinoside. Petals of DhMYB22-treated flowers and sepals of DhMYB60-treated flowers showed the greatest colour difference relative to the same organs in untreated flowers. DhMYB22-treated flowers had relatively narrow and constricted lips, while DhMYB60-treated flowers had narrow and constricted sepals. No significant difference in shape was observed for DhCHS-treated or untreated flowers.

CONCLUSIONS

Our results demonstrate that DhMYB22 and DhMYB60 regulate pigment intensity and floral organ shape in Dendrobium. This is a first report of MYB regulation of floral organ shape in orchids.

Natural alleles of CIRCADIAN CLOCK ASSOCIATED1 contribute to rice cultivation by fine-tuning flowering time

The timing of flowering is a crucial factor for successful grain production at a wide range of latitudes. Domestication of rice (Oryza sativa) included selection for natural alleles of flowering-time genes that allow rice plants to adapt to broad geographic areas. Here, we describe the role of natural alleles of CIRCADIAN CLOCK ASSOCIATED1 (OsCCA1) in cultivated rice based on analysis of single-nucleotide polymorphisms deposited in the International Rice Genebank Collection Information System database. Rice varieties harboring japonica-type OsCCA1 alleles (OsCCA1a haplotype) flowered earlier than those harboring indica-type OsCCA1 alleles (OsCCA1d haplotype). In the japonica cultivar "Dongjin", a T-DNA insertion in OsCCA1a resulted in late flowering under long-day and short-day conditions, indicating that OsCCA1 is a floral inducer. Reverse transcription quantitative PCR analysis showed that the loss of OsCCA1a function induces the expression of the floral repressors PSEUDO-RESPONSE REGULATOR 37 (OsPRR37) and Days to Heading 8 (DTH8), followed by repression of the Early heading date 1 (Ehd1)-Heading date 3a (Hd3a)-RICE FLOWERING LOCUS T 1 (RFT1) pathway. Binding affinity assays indicated that OsCCA1 binds to the promoter regions of OsPRR37 and DTH8. Naturally occurring OsCCA1 alleles are evolutionarily conserved in cultivated rice (O. sativa). Oryza rufipogon-I (Or-I) and Or-III type accessions, representing the ancestors of O. sativa indica and japonica, harbored indica- and japonica-type OsCCA1 alleles, respectively. Taken together, our results demonstrate that OsCCA1 is a likely domestication locus that has contributed to the geographic adaptation and expansion of cultivated rice.

Genome-wide identification and characterization of active ingredients related β-Glucosidases in Dendrobium catenatum

BACKGROUND

Dendrobium catenatum/D. officinale (here after D. catenatum), a well-known economically important traditional medicinal herb, produces a variety of bioactive metabolites including polysaccharides, alkaloids, and flavonoids with excellent pharmacological and clinical values. Although many genes associated with the biosynthesis of medicinal components have been cloned and characterized, the biosynthetic pathway, especially the downstream and regulatory pathway of major medicinal components in the herb, is far from clear. β-glucosidases (BGLUs) comprise a diverse group of enzymes that widely exist in plants and play essential functions in cell wall modification, defense response, phytohormone signaling, secondary metabolism, herbivore resistance, and scent release by hydrolyzing β-D-glycosidic bond from a carbohydrate moiety. The recent release of the chromosome-level reference genome of D. catenatum enables the characterization of gene families. Although the genome-wide analysis of the BGLU gene family has been successfully conducted in various plants, no systematic analysis is available for the D. catenatum. We previously isolated DcBGLU2 in the BGLU family as a key regulator for polysaccharide biosynthesis in D. catenatum. Yet, the exact number of DcBGLUs in the D. catenatum genome and their possible roles in bioactive compound production deserve more attention.

RESULTS

To investigate the role of BGLUs in active metabolites production, 22 BGLUs (DcBGLU1-22) of the glycoside hydrolase family 1 (GH1) were identified from D. catenatum genome. Protein prediction showed that most of the DcBGLUs were acidic and phylogenetic analysis classified the family into four distinct clusters. The sequence alignments revealed several conserved motifs among the DcBGLU proteins and analyses of the putative signal peptides and N-glycosylation site revealed that the majority of DcBGLU members dually targeted to the vacuole and/or chloroplast. Organ-specific expression profiles and specific responses to MeJA and MF23 were also determined. Furthermore, four DcBGLUs were selected to test their involvement in metabolism regulation. Overexpression of DcBGLU2, 6, 8, and 13 significantly increased contents of flavonoid, reducing-polysaccharide, alkaloid and soluble-polysaccharide, respectively.

CONCLUSION

The genome-wide systematic analysis identified candidate DcBGLU genes with possible roles in medicinal metabolites production and laid a theoretical foundation for further functional characterization and molecular breeding of D. catenatum.

Genome-Wide Analysis of the WOX Transcription Factor Genes in Dendrobium catenatum Lindl

The WUSCHEL-related homeobox (WOX) proteins are a class of transcription factors exclusive to plants. They can promote cell division or inhibit stem cell differentiation to regulate plant growth and development. However, the WOX transcription factor genes in the monocotyledon Dendrobium catenatum Lindl. remain relatively uncharacterized. Specifically, the effects of phytohormones on their expression levels are unclear. In this study, we identified and analyzed 10 candidate DcaWOX transcription factor genes in D. catenatum. The DcaWOX family was divided into the modern/WUS, intermediate, and ancient clades. The subcellular localization analysis detected DcaWOX-GFP fusion proteins in the tobacco epidermal leaf cell nucleus. In DcaWOX, members of the WUS clade with the WUS-box motif can significantly activate the expression of TPL in vivo, while members of the intermediate and ancient clades cannot. The expression of the DcaWOX genes varied among the examined tissues. Moreover, the DcaWOX expression patterns were differentially affected by the phytohormone treatments, with differences detected even between homologs of the same gene. Furthermore, the gene expression patterns were consistent with the predicted cis-acting elements in the promoters. The above results suggest that DcaWOX may have an important role in its growth and development and resistance to stress. The results of this comprehensive investigation of the DcaWOX gene family provide the basis for future studies on the roles of WOX genes in D. catenatum.

Genome-Wide Analysis of miR159 Gene Family and Predicted Target Genes Associated with Environmental Stress in Dendrobium officinale: A Bioinformatics Study

Dendrobium officinale (D. officinale) is a widely used traditional Chinese medicine with high economic value. MicroR159 (miR159) is an ancient and conserved microRNA (miRNA) family in land plants, playing roles in the progress of growth and development, as well as the stress response. In order to find out functions of miR159 in D. officinale, multiple bioinformatic approaches were employed and 10 MIR159 genes were found, localizing on seven chromosomes and an unanchored segment of the D. officinale genome. All of the precursor sequences of Dof-miR159 could form a stable stem-loop structure. The phylogenetic analysis revealed that the MIR159 genes of D. officinale were divided into five clades. Furthermore, the conservation analysis suggested that the 2 to 20 nt region of miR159 mature sequences were highly conserved among family members. The promoter analysis of MIR159s showed that the majority of the predicted cis-elements were related to environmental stress or hormones. In total, five classes of genes were predicted to be the target genes of Dof-miR159s, including GAMYB transcription factors, which had been confirmed in many other land plants. The expression patterns of predicted target genes revealed their potential roles in the growth and development of D. officinale, as well as in cold and drought stress responses. In conclusion, our results illustrated the stress-related miR159-targeted genes in D. officinale, which could provide candidate genes for resistance breeding in the future.

Identification and Expression Analysis of WRKY Gene Family in Response to Abiotic Stress in Dendrobium catenatum

Dendrobium catenatum has become a rare and endangered medicinal plant due to habitat loss in China. As one of the most important and largest transcription factors, WRKY plays a critical role in response to abiotic stresses in plants. However, little is known regarding the functions of the WRKY family in D. catenatum. In this study, a total of 62 WRKY genes were identified from the D. catenatum genome. Phylogenetic analysis revealed that DcWRKY proteins could be divided into three groups, a division supported by the conserved motif compositions and intron/exon structures. DcWRKY gene expression and specific responses under drought, heat, cold and salt stresses were analyzed through RNA-seq data and RT-qPCR assay. The results showed that these genes had tissue-specificity and displayed different expression patterns in response to abiotic stresses. The expression levels of DcWRKY22, DcWRKY36 and DcWRKY45 were up-regulated by drought stress. Meanwhile, DcWRKY22 was highly induced by heat in roots, and DcWRKY45 was significantly induced by cold stress in leaves. Furthermore, DcWRKY27 in roots and DcWRKY58 in leaves were extremely induced under salt treatment. Finally, we found that all the five genes may function in ABA- and SA-dependent manners. This study identified candidate WRKY genes with possible roles in abiotic stress and these findings not only contribute to our understanding of WRKY family genes, but also provide valuable information for stress resistance development in D. catenatum.

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.

Genome-Wide Identification of the MYB Gene Family in Cymbidiumensifolium and Its Expression Analysis in Different Flower Colors

MYB transcription factors of plants play important roles in flavonoid synthesis, aroma regulation, floral organ morphogenesis, and responses to biotic and abiotic stresses. Cymbidium ensifolium is a perennial herbaceous plant belonging to Orchidaceae, with special flower colors and high ornamental value. In this study, a total of 136 CeMYB transcription factors were identified from the genome of C. ensifolium, including 27 1R-MYBs, 102 R2R3-MYBs, 2 3R-MYBs, 2 4R-MYBs, and 3 atypical MYBs. Through phylogenetic analysis in combination with MYB in Arabidopsis thaliana, 20 clusters were obtained, indicating that these CeMYBs may have a variety of biological functions. The 136 CeMYBs were distributed on 18 chromosomes, and the conserved domain analysis showed that they harbored typical amino acid sequence repeats. The motif prediction revealed that multiple conserved elements were mostly located in the N-terminal of CeMYBs, suggesting their functions to be relatively conserved. CeMYBs harbored introns ranging from 0 to 13 and contained a large number of stress- and hormone-responsive cis-acting elements in the promoter regions. The subcellular localization prediction demonstrated that most of CeMYBs were positioned in the nucleus. The analysis of the CeMYBs expression based on transcriptome data showed that CeMYB52, and CeMYB104 of the S6 subfamily may be the key genes leading to flower color variation. The results lay a foundation for the study of MYB transcription factors of C. ensifolium and provide valuable information for further investigations of the potential function of MYB genes in the process of anthocyanin biosynthesis.