Genome-wide analysis of the R2R3-MYB transcription factor genes in Chinese cabbage (Brassica rapa ssp. pekinensis) reveals their stress and hormone responsive patterns

Plant Signaling Hormones and Transcription Factors: Key Regulators of Plant Responses to Growth, Development, and Stress

Plants constantly encounter a wide range of biotic and abiotic stresses that adversely affect their growth, development, and productivity. Phytohormones such as abscisic acid, jasmonic acid, salicylic acid, and ethylene serve as crucial regulators, integrating internal and external signals to mediate stress responses while also coordinating key developmental processes, including seed germination, root and shoot growth, flowering, and senescence. Transcription factors (TFs) such as WRKY, NAC, MYB, and AP2/ERF play complementary roles by orchestrating complex transcriptional reprogramming, modulating stress-responsive genes, and facilitating physiological adaptations. Recent advances have deepened our understanding of hormonal networks and transcription factor families, revealing their intricate crosstalk in shaping plant resilience and development. Additionally, the synthesis, transport, and signaling of these molecules, along with their interactions with stress-responsive pathways, have emerged as critical areas of study. The integration of cutting-edge biotechnological tools, such as CRISPR-mediated gene editing and omics approaches, provides new opportunities to fine-tune these regulatory networks for enhanced crop resilience. By leveraging insights into transcriptional regulation and hormone signaling, these advancements provide a foundation for developing stress-tolerant, high-yielding crop varieties tailored to the challenges of climate change.

Genome-Wide Analysis of the Common Fig (Ficus carica L.) R2R3-MYB Genes Reveals Their Structure, Evolution, and Roles in Fruit Color Variation

The R2R3-MYB transcription factor (TF) family is crucial for regulating plant growth, stress response, and fruit ripening. Although this TF family has been examined in a multitude of plants, the R2R3-MYB TFs in Ficus carica, a Mediterranean fruit species, have yet to be characterized. This study identified and classified 63 R2R3-MYB genes (FcMYB1 to FcMYB63) in the F. carica genome. We analyzed these genes for physicochemical properties, conserved motifs, phylogenetic relationships, gene architecture, selection pressure, and gene expression profiles and networks. The genes were classified into 29 clades, with members of the same clade showing similar exon-intron structures and motif compositions. Of the 54 orthologous gene pairs shared with mulberry (Morus notabilis), 52 evolved under negative selection, while two pairs (FcMYB55/MnMYB20 and FcMYB59/MnMYB31) experienced diversifying selection. RNA-Seq analysis showed that FcMYB26, FcMYB33, and FcMYB34 were significantly overexpressed in fig fruit peel during maturation phase III. Weighted gene co-expression network analysis (WGCNA) indicated that these genes are part of an expression module associated with the anthocyanin pathway. RT-qPCR validation confirmed these findings and revealed that the Tunisian cultivars 'Zidi' and 'Soltani' have cultivar-specific R2R3-FcMYB genes highly overexpressed during the final stage of fruit maturation and color acquisition. These genes likely influence cultivar-specific pigment synthesis. This study provides a comprehensive overview of the R2R3-MYB TF family in fig, offering a framework for selecting genes related to fruit peel color in breeding programs.

A R2R3-MYB transcription factor, FeR2R3-MYB, positively regulates anthocyanin biosynthesis and drought tolerance in common buckwheat (Fagopyrum esculentum)

The R2R3-MYB transcription factors (TFs) play a crucial role in regulating plant secondary metabolism and abiotic stress responses, yet they are still poorly understood in common buckwheat (Fagopyrum esculentum), a valuable minor grain crop resource. In this study, a candidate gene, FeR2R3-MYB, was cloned from the anthocyanin-rich common buckwheat variety 'QZZTQ'. FeR2R3-MYB was found to contain two MYB DNA-binding domains and be located at the nucleus with transcriptional activation activity. Molecular analysis indicated that FeR2R3-MYB is predominantly expressed in flowering tissue and is highly responsive to environmental factors such as light, drought, and cold. In addition, the promoter of FeR2R3-MYB showed a positive correlation with fragment length. Further functional analysis suggested that FeR2R3-MYB not only participates in the anthocyanin biosynthetic pathway by interacting with leucoanthocyanidin reductase (FeLAR), but also enhances drought tolerance in common buckwheat. To sum up, FeR2R3-MYB exhibits positive effects on both pigment production (e.g., anthocyanin) and abiotic stress resistance, providing valuable insights for future research in buckwheat molecular breeding and resource development.

Multi-omics landscape to decrypt the distinct flavonoid biosynthesis of Scutellaria baicalensis across multiple tissues

Scutellaria baicalensis Georgi, also known as huang-qin in traditional Chinese medicine, is a widely used herbal remedy due to its anticancer, antivirus, and hepatoprotective properties. The S. baicalensis genome was sequenced many years ago; by contrast, the proteome as the executer of most biological processes of S. baicalensis in the aerial parts, as well as the secondary structure of the roots (xylem, phloem, and periderm), is far less comprehensively characterized. Here we attempt to depict the molecular landscape of the non-model plant S. baicalensis through a multi-omics approach, with the goal of constructing a highly informative and valuable reference dataset. Furthermore, we provide an in-depth characterization dissection to explain the two distinct flavonoid biosynthesis pathways that exist in the aerial parts and root, at the protein and phosphorylated protein levels. Our study provides detailed spatial proteomic and phosphoproteomic information in the context of secondary structures, with implications for the molecular profiling of secondary metabolite biosynthesis in non-model medicinal plants.

Multi-Omic Analysis Reveals the Molecular Mechanism of UV-B Stress Resistance in Acetylated RcMYB44 in Rhododendron chrysanthum

Ultraviolet-B (UV-B) radiation is a significant environmental factor influencing the growth and development of plants. MYBs play an essential role in the processes of plant responses to abiotic stresses. In the last few years, the development of transcriptome and acetylated proteome technologies have resulted in further and more reliable data for understanding the UV-B response mechanism in plants. In this research, the transcriptome and acetylated proteome were used to analyze Rhododendron chrysanthum Pall. (R. chrysanthum) leaves under UV-B stress. In total, 2348 differentially expressed genes (DEGs) and 685 differentially expressed acetylated proteins (DAPs) were found. The transcriptome analysis revealed 232 MYB TFs; we analyzed the transcriptome together with the acetylated proteome, and screened 4 MYB TFs. Among them, only RcMYB44 had a complete MYB structural domain. To investigate the role of RcMYB44 under UV-B stress, a homology tree was constructed between RcMYB44 and Arabidopsis MYBs, and it was determined that RcMYB44 shares the same function with ATMYB44. We further constructed the hormone signaling pathway involved in RcMYB44, revealing the molecular mechanism of resistance to UV-B stress in R. chrysanthum. Finally, by comparing the transcriptome and the proteome, it was found that the expression levels of proteins and genes were inconsistent, which is related to post-translational modifications of proteins. In conclusion, RcMYB44 of R. chrysanthum is involved in mediating the growth hormone, salicylic acid, jasmonic acid, and abscisic acid signaling pathways to resist UV-B stress.

Identification of Vitis vinifera MYB transcription factors and their response against grapevine berry inner necrosis virus

BACKGROUND

The myeloblastosis (MYB) superfamily is the largest transcription factor family in plants that play diverse roles during stress responses. However, the biotic stress-responsive MYB transcription factors of the grapevine have not been systematically studied. In China, grapevine berries are often infected with the grapevine berry inner necrosis virus (GINV), which eventually reduces the nutritional quality and commodity value.

RESULTS

The present study identified and characterized 265 VvMYB or VvMYB-related genes of the "Crimson seedless" grapevine. Based on DNA-binding domain analysis, these VvMYB proteins were classified into four subfamilies, including MYB-related, 2R-MYB, 3R-MYB, and 4R-MYB. Phylogenetic analysis divided the MYB transcription factors into 26 subgroups. Overexpression of VvMYB58 suppressed GINV abundance in the grapevine. Further qPCR indicated that among 41 randomly selected VvMYB genes, 12 were induced during GINV infection, while 28 were downregulated. These findings suggest that VvMYB genes actively regulate defense response in the grapevine.

CONCLUSION

A deeper understanding of the MYB TFs engaged in GINV defense response will help devise better management strategies. The present study also provides a foundation for further research on the functions of the MYB transcription factors.

Genome-Wide Analysis of MYB Transcription Factor Gene Superfamily Reveals BjPHL2a Involved in Modulating the Expression of BjCHI1 in Brassica juncea

Brassica juncea is an economically important vegetable and oilseed crop. The MYB transcription factor superfamily is one of the largest transcription factor families in plants, and plays crucial roles in regulating the expression of key genes involved in a variety of physiological processes. However, a systematic analysis of the MYB transcription factor genes in Brassica juncea (BjMYB) has not been performed. In this study, a total of 502 BjMYB superfamily transcription factor genes were identified, including 23 1R-MYBs, 388 R2R3-MYBs, 16 3R-MYBs, 4 4R-MYBs, 7 atypical MYBs, and 64 MYB-CCs, which is approximately 2.4-fold larger than that of AtMYBs. Phylogenetic relationship analysis revealed that the MYB-CC subfamily consists of 64 BjMYB-CC genes. The expression pattern of members of PHL2 subclade homologous genes in Brassica juncea (BjPHL2) after Botrytis cinerea infection were determined, and BjPHL2a was isolated from a yeast one-hybrid screen with the promoter of BjCHI1 as bait. BjPHL2a was found to localize mainly in the nucleus of plant cells. An EMSA assay confirmed that BjPHL2a binds to the Wbl-4 element of BjCHI1. Transiently expressed BjPHL2a activates expression of the GUS reporter system driven by a BjCHI1 mini-promoter in tobacco (Nicotiana benthamiana) leaves. Taken together, our data provide a comprehensive evaluation of BjMYBs and show that BjPHL2a, one of the members of BjMYB-CCs, functions as a transcription activator by interacting with the Wbl-4 element in the promoter of BjCHI1 for targeted gene-inducible expression.

Identification and Phylogenetic Analysis of the R2R3-MYB Subfamily in Brassica napus

The R2R3-MYB sub-family proteins are composed of most members of MYB (v-Myb avian myeloblastosis viral oncogene homolog) protein, a plant-specific transcription factor (TF) that is classified into four classes depending on the number of MYB repeats. R2R3-MYB TFs are involved in physiological and biochemical processes. However, the functions of the Brassica napus R2R3-MYB genes are still mainly unknown. In this study, 35 Brassica napus MYB (BnaMYB) genes were screened in the genome of Brassica napus, and details about their physical and chemical characteristics, evolutionary relationships, chromosome locations, gene structures, three-dimensional protein structures, cis-acting promoter elements, and gene duplications were uncovered. The BnaMYB genes have undergone segmental duplications and positive selection pressure, according to evolutionary studies. The same subfamilies have similar intron-exon patterns and motifs, according to the genes' structure and conserved motifs. Additionally, through cis-element analysis, many drought-responsive and other stress-responsive cis-elements have been found in the promoter regions of the BnaMYB genes. The expression of the BnaMYB gene displays a variety of tissue-specific patterns. Ten lignin-related genes were chosen for drought treatment. Our research screened four genes that showed significant upregulation under drought stress, and thus may be important drought-responsive genes. The findings lay a new foundation for understanding the complex mechanisms of BnaMYB in multiple developmental stages and pathways related to drought stress in rapeseed.

Genome-Wide Characterization and Analysis of R2R3-MYB Genes Related to Fruit Ripening and Stress Response in Banana (Musa acuminata L. AAA Group, cv. 'Cavendish')

MYB is an important type of transcription factor in eukaryotes. It is widely involved in a variety of biological processes and plays a role in plant morphogenesis, growth and development, primary and secondary metabolite synthesis, and other life processes. In this study, bioinformatics methods were used to identify the R2R3-MYB transcription factor family members in the whole Musa acuminata (DH-Pahang) genome, one of the wild ancestors of banana. A total of 280 MaMYBs were obtained, and phylogenetic analysis indicated that these MaMYBs could be classified into 33 clades with MYBs from Arabidopsis thaliana. The amino acid sequences of the R2 and R3 Myb-DNA binding in all MaMYB protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 277 MaMYBs were localized on the 11 chromosomes in the Musa acuminata genome. The MaMYBs were distributed unevenly across the 11 chromosomes. More than 40.0% of the MaMYBs were located in collinear fragments, and segmental duplications likely played a key role in the expansion of the MaMYBs. Moreover, the expression profiles of MaMYBs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Weighted gene co-expression network analysis (WGCNA) was used to analyze transcriptome data of banana from the above 11 samples. We found MaMYBs participating in important metabolic biosynthesis pathways in banana. Collectively, our results represent a comprehensive genome-wide study of the MaMYB gene family, which should be helpful in further detailed studies on MaMYBs functions related to fruit development, postharvest ripening, and the seedling response to stress in an important banana cultivar.

A Novel R2R3-MYB Transcription Factor SbMYB12 Positively Regulates Baicalin Biosynthesis in Scutellaria baicalensis Georgi

Scutellaria baicalensis Georgi is an annual herb from the Scutellaria genus that has been extensively used as a traditional medicine for over 2000 years in China. Baicalin and other flavonoids have been identified as the principal bioactive ingredients. The biosynthetic pathway of baicalin in S. baicalensis has been elucidated; however, the specific functions of R2R3-MYB TF, which regulates baicalin synthesis, has not been well characterized in S. baicalensis to date. Here, a S20 R2R3-MYB TF (SbMYB12), which encodes 263 amino acids with a length of 792 bp, was expressed in all tested tissues (mainly in leaves) and responded to exogenous hormone methyl jasmonate (MeJA) treatment. The overexpression of SbMYB12 significantly promoted the accumulation of flavonoids such as baicalin and wogonoside in S. baicalensis hairy roots. Furthermore, biochemical experiments revealed that SbMYB12 is a nuclear-localized transcription activator that binds to the SbCCL7-4, SbCHI-2, and SbF6H-1 promoters to activate their expression. These results illustrate that SbMYB12 positively regulates the generation of baicalin and wogonoside. In summary, this work revealed a novel S20 R2R3-MYB regulator and enhances our understanding of the transcriptional and regulatory mechanisms of baicalin biosynthesis, as well as sheds new light on metabolic engineering in S. baicalensis.

Comprehensive Analysis of Betula platyphylla Suk. PIF Gene Family and Their Potential Functions in Growth and Development

Phytochrome-interacting factors (PIFs) are transcription factors with the basic helix-loop-helix (bHLH) domain. As integration factors between different signal pathways, members of the PIF protein family regulate many aspects of plant growth and development, such as seed germination, photomorphogenesis, thermomorphogenesis, rhythm regulation, flowering response, stomatal development, and stress responses. Our previous studies have shown that the BpSPL2 gene may regulate plants' adventitious root development through PIF genes. Within the Betula platyphylla genome, we identified eight PIF (BpPIFs) genes. We analysed and named them based on a phylogenetic tree, gene structures, and conserved motifs. Synteny analysis indicated that transposition or segmental duplication events played a minor role in the expansion of BpPIFs. The comparative syntenic analysis combined with phylogenetic analysis provided a deep insight into the phylogenetic relationships of BpPIF genes, suggesting that BpPIF proteins are closer to PtPIF than to AtPIF. The analysis of cis-acting elements in promoter regions of BpPIF genes indicated that various elements were related to light, abiotic stress, and plant hormone responsiveness. In addition, we found that these promoters have the transcription factor of B. platyphylla SPL2 (BpSPL2) binding motif GTAC. Expression analysis demonstrated that BpPIF genes, especially BpPIF4, BpPIF9b, and BpPIF10, might be the potential target genes of BpSPL2 in the process of adventitious root formation. Besides providing a comprehensive understanding of the BpPIF family, we propose a hypothetical gene network regulatory model for adventitious root formation.

Genomic Characteristics of Elite Maize Inbred Line 18-599 and Its Transcriptional Response to Drought and Low-Temperature Stresses

Elite inbred line 18-599 was developed via triple test cross from introduced hybrid P78599 and used as parents of dozens of maize hybrids adapting to the diverse ecological conditions of the maize ecological region in Southwest China. In this study, its genomic DNA was resequenced and aligned with the B73 genome sequence to identify single nucleotide polymorphism (SNP), and insertion (In) and deletion (Del) loci. These loci were aligned with those between B73 and 1020 inbred lines in the HapMap database to identify specific variation loci of 18-599. The results showed that there were 930,439 specific SNPs and 358,750 InDels between 18-599 and the 1020 lines. In total, 21,961 of them showed significant impacts on the functions of 12,297 genes, such as frameshift, change of splicing site, stop gain, change of start site, and stop loss. Phylogenetic analysis showed that 18-599 was closely related to inbred lines ZEAxujRAUDIAAPE and 2005-4, but far from some inbred lines directly isolated from P78599. This result indicated that 18-599 not only pyramided the elite genes of P78599, but also acquired genetic divergence during the repetitive backcrosses of triple test cross to confer its elite agronomic characteristics. Subsequently, the RNA of 18-599 was sequenced. The aligned 9713 and 37,528 of the 165,098 unigenes were screened and aligned with annotated transcripts of the B73 genome differentially expressed under drought and low-temperature stress, respectively, and their functions were involved in the responses to these stresses. The quantitative PCR results of fourteen random genes verified the RNA sequencing results. These findings suggest that the transcriptional responses of many resistance-related genes were an important mechanism for 18-599 to adapt to diverse ecological conditions.

Transcriptome Analysis of Chinese Cabbage Provides Insights into the Basis of Understanding the Lignin Affected by Low Temperature

Chinese cabbage, which is a cold season crop, can still be damaged at an overly low temperature. It is crucial to study the mechanism of the resistance to low temperature of Chinese cabbage. In this study, the Chinese cabbage 'XBJ' was used as the material, and nine different low temperatures and control samples were treated. Using RNA-seq and lignin content determination, we analyzed 27 samples, and the stained sections of them were observed. A total of 8845 genes were screened for the WGCNA analysis, yielding 17 modules. The GO and KEGG analyses of the modules was highly associated with a low-temperature treatment. The pathways such as 'starch and sucrose metabolism' and 'plant hormone signal transduction' were enriched in modules related to low temperature. Interestingly, L-15DAT-associated MEcoral2 was found to have 14 genes related to the 'lignin biosynthetic process' in the GO annotation. The combination of the determination of the lignin content and the treatment of the stained sections showed that the lignin content of the low-temperatures samples were indeed higher than that of the control. We further explored the expression changes of the lignin synthesis pathway and various genes and found that low temperature affects the expression changes of most genes in the lignin synthesis pathway, leading to the speculation that the lignin changes at low temperature are a defense mechanism against low temperatures. The 29 BrCOMT gene sequence derived from the RNA-seq was non-conserved, and eight BrCOMT genes were differentially expressed. This study provides a new insight into how lignin is affected by low temperature.

IbMYB308, a Sweet Potato R2R3-MYB Gene, Improves Salt Stress Tolerance in Transgenic Tobacco

The MYB (v-myb avian myeloblastosis viral oncogene homolog) transcription factor family plays an important role in plant growth, development, and response to biotic and abiotic stresses. However, the gene functions of MYB transcription factors in sweet potato (Ipomoea batatas (L.) Lam) have not been elucidated. In this study, an MYB transcription factor gene, IbMYB308, was identified and isolated from sweet potato. Multiple sequence alignment showed that IbMYB308 is a typical R2R3-MYB transcription factor. Further, quantitative real-time PCR (qRT-PCR) analysis revealed that IbMYB308 was expressed in root, stem, and, especially, leaf tissues. Moreover, it showed that IbMYB308 had a tissue-specific profile. The experiment also showed that the expression of IbMYB308 was induced by different abiotic stresses (20% PEG-6000, 200 mM NaCl, and 20% H₂O₂). After a 200 mM NaCl treatment, the expression of several stress-related genes (SOD, POD, APX, and P5CS) was upregulation in transgenic plants, and the CAT activity, POD activity, proline content, and protein content in transgenic tobacco had increased, while MDA content had decreased. In conclusion, this study demonstrated that IbMYB308 could improve salt stress tolerance in transgenic tobacco. These findings lay a foundation for future studies on the R2R3-MYB gene family of sweet potato and suggest that IbMYB308 could potentially be used as an important positive factor in transgenic plant breeding to improve salt stress tolerance in sweet potato plants.

The Pea R2R3-MYB Gene Family and Its Role in Anthocyanin Biosynthesis in Flowers

Pea (Pisum sativum L.) is one of the most important legume crops in the world, and it has attracted great attention for its high nutritive values. Recently, the crop breeding program has been focused on the crop metabolic engineering (i.e., color, flavor, nutrition) to improve the quality of crop. As a major group of transcription factors forming the ternary MYB–bHLH–WD repeat protein (MBW) complex to regulate the anthocyanin biosynthesis pathway, members of R2R3-MYB gene family have always been the focus of research targets to improve the valuable metabolic product of crops. Until now, few report about the R2R3-MYB gene family of pea has been released. In this study, we identified 119 R2R3-MYB genes in the assembled pea genome (Version 1a), of which 111 were distributed across 14 chromosomes. Combining with the 126 R2R3-MYB protein sequences of Arabidopsis, we categorized 245 R2R3-MYB proteins into 36 subgroups according to sequence similarity and phylogenetic relationships. There was no member from subgroup 12, 15 and 29 existing in pea genome, whereas three novel subgroups were found in pea and named as N1-N3. Further analyses of conserved domains and Motifs, gene structures, and chromosomal locations showed that the typical R2 and R3 domains were present across all R2R3-MYB proteins, and Motif 1, 2, and 3 were identified in most members. Most of them had no more than two introns. Additionally, 119 pea R2R3-MYB genes did not experience large-scale duplication events. Finally, we concluded that several candidate genes may be responsible for the spatiotemporal accumulation of anthocyanins in pea petals. PsMYB116 was predominantly expressed in the dorsal petals to presumably activate the anthocyanin biosynthesis pathway, while PsMYB37 and PsMYB32 may positively regulates the anthocyanin accumulation in the lateral petals. This study not only provides a good reference to further characterize the diverse functions of R2R3-MYB genes but also helps researchers to understand the color formation of pea flowers.

Genome-wide identification of R2R3-MYB gene family and association with anthocyanin biosynthesis in Brassica species

Brassica species include important oil crops and vegetables in the world. The R2R3-MYB gene participates in a variety of plant functions, including the activation or inhibition of anthocyanin biosynthesis. Although previous studies have reported its phylogenetic relationships, gene structures, and expression patterns in Arabidopsis, the number and sequence variation of this gene family in Brassica crops and its involvement in the natural quantitative variation in anthocyanin biosynthesis regulation are still largely unknown. In this study, by using whole genome sequences and comprehensive genome-wide comparative analysis among the six cultivated Brassica species, 2120 R2R3-MYB genes were identified in six Brassica species, in total These R2R3-MYB genes were phylogenetically clustered into 12 groups. The R2R3-MYB family between A and C subgenomes showed better collinearity than between B and C and between A and B. From comparing transcriptional changes of five Brassica species with the purple and green leaves for the detection of the R2R3-MYB genes associated with anthocyanin biosynthesis, 7 R2R3-MYB genes were co-differentially expressed. The promoter and structure analysis of these genes showed that some variations between non-coding region, but they were highly conserved at the protein level and spatial structure. Co-expression analysis of anthocyanin-related genes and R2R3-MYBs indicated that MYB90 was strongly co-expressed with TT8, and they were co-expressed with structural genes F3H, LDOX, ANS and UF3GT at the same time. These results further clarified the roles of the R2R3-MYBs for leaf coloration in Brasica species, which provided new insights into the functions of the R2R3-MYB gene family in Brasica species.

Genome wide characterization of R2R3 MYB transcription factor from Apocynum venetum revealed potential stress tolerance and flavonoid biosynthesis genes

MYB transcription factors are crucial in regulating stress tolerance and expression of major genes involved in flavonoid biosynthesis. The functions of MYBs is well explored in a number of plants, yet no studies is reported in Apocynum venetum. We identified a total of 163 MYB candidates, that comprised of 101 (61.96%) R2R3, 6 3R, 1 4R and 55 1R. Syntenic analysis of A. venetum R2R3 (AvMYB) showed highest orthologous pairs with Vitis vinifera MYBs followed by Arabidopsis thaliana among the four species evaluated. Thirty segmental duplications and 6 tandem duplications were obtained among AvMYB gene pairs signifying their role in the MYB gene family expansion. Nucleotide substitution analysis (Ka/Ks) showed the AvMYBs to be under the influence of strong purifying selection. Expression analysis of selected AvMYB under low temperature and cadmium stresses resulted in the identification of AvMYB48, AvMYB97, AvMYB8,AvMYB4 as potential stress responsive genes and AvMYB10 and AvMYB11 in addition, proanthocyanidin biosynthesis regulatory genes which is consistent with their annotated homologues in Arabidopsis. Tissue specific expression profile analysis of AvMYBs further supported the qPCR analysis result. MYBs with higher transcript levels in root, stem and leaf like AvMYB4 forexample, was downregulated under the stresses and such with low transcript level such as AvMYB48 which had low transcript in the leaf was upregulated under both stresses. Transcriptome and phylogenetic analysis suggested AvMYB42 as a potential regulator of anthocyanin biosynthesis. Thus, this study provided valuable information on AvR2R3-MYB gene family with respect to stress tolerance and flavonoid biosynthesis.

The BrGI Circadian Clock Gene Is Involved in the Regulation of Glucosinolates in Chinese Cabbage

Circadian clocks integrate environmental cues with endogenous signals to coordinate physiological outputs. Clock genes in plants are involved in many physiological and developmental processes, such as photosynthesis, stomata opening, stem elongation, light signaling, and floral induction. Many Brassicaceae family plants, including Chinese cabbage (Brassica rapa ssp. pekinensis), produce a unique glucosinolate (GSL) secondary metabolite, which enhances plant protection, facilitates the design of functional foods, and has potential medical applications (e.g., as antidiabetic and anticancer agents). The levels of GSLs change diurnally, suggesting a connection to the circadian clock system. We investigated whether circadian clock genes affect the biosynthesis of GSLs in Brassica rapa using RNAi-mediated suppressed transgenic Brassica rapa GIGENTEA homolog (BrGI knockdown; hereafter GK1) Chinese cabbage. GIGANTEA plays an important role in the plant circadian clock system and is related to various developmental and metabolic processes. Using a validated GK1 transgenic line, we performed RNA sequencing and high-performance liquid chromatography analyses. The transcript levels of many GSL pathway genes were significantly altered in GK1 transgenic plants. In addition, GSL contents were substantially reduced in GK1 transgenic plants. We report that the BrGI circadian clock gene is required for the biosynthesis of GSLs in Chinese cabbage plants.

Crocus transcription factors CstMYB1 and CstMYB1R2 modulate apocarotenoid metabolism by regulating carotenogenic genes

Two MYB genes have been identified which regulate apocarotenoid metabolism in Crocus sativus. Apocarotenoids like crocin, picrocrocin and safranal are restricted to genus Crocus and are synthesized by oxidative cleavage of zeaxanthin followed by glycosylation reactions. In Crocus sativus, these apocarotenoids are synthesized in stigma part of the flower in developmentally regulated manner. Most of the genes of apocarotenoid pathway are known, however, the mechanism that regulates its tissue and stage specific biosynthesis remains elusive. MYB family was identified as the largest transcription factor family from Crocus transciptome which indicated its possible role in apocarotenoid regulation besides regulating other metabolic pathways. Towards this, we started with identification of 150 MYB genes from Crocus transcriptome databases. The phylogenetic analysis of Crocus MYB genes divided them into 27 clusters. Domain analysis resulted in identification of four groups of MYBs depending upon the number of R repeats present. Expression profiling indicated that 12 MYBs are upregulated in stigma out of which expression of four genes CstMYB1, CstMYB14, CstMYB16 and CstMYB1R2 correlated with crocin accumulation. Transient overexpression of two nuclear localized MYB genes (CstMYB1 and CstMYB1R2) in Crocus confirmed their role in regulating carotenoid metabolism. Yeast-one-hybrid confirmed that CstMYB1 binds to carotenoid cleavage dioxygenase 2 (CCD2) promoter while CstMYB1R2 binds to phytoene synthase (PSY) and CCD2 promoters. Overall, our study established that CstMYB1 and CstMYB1R2 regulate apocarotenoid biosynthesis by directly binding to promoters of pathway genes.

Genome-Wide Analysis of the R2R3-MYB Gene Family in Fragaria × ananassa and Its Function Identification During Anthocyanins Biosynthesis in Pink-Flowered Strawberry

The strawberry (Fragaria × ananassa) is an economically important fruit throughout the world. The large R2R3-MYB gene family participates in a variety of plant functions, including anthocyanin biosynthesis. The present study is the first genome-wide analysis of the MYB gene family in the octoploid strawberry and describes the identification and characterization of the family members using the recently sequenced F. × ananassa genome. Specifically, we aimed to identify the key MYBs involved in petal coloration in the pink-flowered strawberry, which increases its ornamental value. A comprehensive, genome-wide analysis of F. × ananassa R2R3-FaMYBs was performed, investigating gene structures, phylogenic relationships, promoter regions, chromosomal locations, and collinearity. A total of 393 R2R3-FaMYB genes were identified in the F. × ananassa genome and divided into 36 subgroups based on phylogenetic analysis. Most genes with similar functions in the same subgroup exhibited similar exon-intron structures and motif compositions. These R2R3-FaMYBs were unevenly distributed over 28 chromosomes. The expansion of the R2R3-FaMYB gene family in the F. × ananassa genome was found to be caused mainly by segmental duplication. The Ka/Ks analysis indicated that duplicated R2R3-FaMYBs mostly experienced purifying selection and showed limited functional divergence after the duplication events. To elucidate which R2R3-FaMYB genes were associated with anthocyanin biosynthesis in the petals of the pink-flowered strawberry, we compared transcriptional changes in different flower developmental stages using RNA-seq. There were 131 differentially expressed R2R3-FaMYB genes identified in the petals, of which three genes, FaMYB28, FaMYB54, and FaMYB576, appeared likely, based on the phylogenetic analysis, to regulate anthocyanin biosynthesis. The qRT-PCR showed that these three genes were more highly expressed in petals than in other tissues (fruit, leaf, petiole and stolon) and their expressions were higher in red compared to pink and white petals. These results facilitate the clarification on the roles of the R2R3-FaMYB genes in petal coloration in the pink-flowered strawberry. This work provides useful information for further functional analysis on the R2R3-FaMYB gene family in F. × ananassa.

Genome-Wide Characterization of R2R3-MYB Transcription Factors in Pitaya Reveals a R2R3-MYB Repressor HuMYB1 Involved in Fruit Ripening through Regulation of Betalain Biosynthesis by Repressing Betalain Biosynthesis-Related Genes

The MYB (myeloblastosis) superfamily constitutes one of the most abundant transcription factors (TFs) regulating various biological processes in plants. However, the molecular characteristics and functions of MYB TFs in pitaya remain unclear. To date, no genome-wide characterization analysis of this gene family has been conducted in the Cactaceae species. In this study, 105 R2R3-MYB members were identified from the genome data of Hylocereus undatus and their conserved motifs, physiological and biochemical characteristics, chromosome locations, synteny relationship, gene structure and phylogeny were further analyzed. Expression analyses suggested that three up-regulated HuMYBs and twenty-two down-regulated HuMYBs were probably involved in fruit ripening of pitaya. Phylogenetic analyses of R2R3-MYB repressors showed that seven HuMYBs (HuMYB1, HuMYB21, HuMYB48, HuMYB49, HuMYB72, HuMYB78 and HuMYB101) were in clades containing R2R3-MYB repressors. HuMYB1 and HuMYB21 were significantly down-regulated with the betalain accumulation during fruit ripening of ‘Guanhuahong’ pitaya (H. monacanthus). However, only HuMYB1 had R2 and R3 repeats with C1, C2, C3 and C4 motifs. HuMYB1 was localized exclusively to the nucleus and exhibited transcriptional inhibition capacities. Dual luciferase reporter assay demonstrated that HuMYB1 inhibited the expression of betalain-related genes: HuADH1, HuCYP76AD1-1 and HuDODA1. These results suggested that HuMYB1 is a potential repressor of betalain biosynthesis during pitaya fruit ripening. Our results provide the first genome-wide analyses of the R2R3-MYB subfamily involved in pitaya betalain biosynthesis and will facilitate functional analysis of this gene family in the future.

Evolution of the Subgroup 6 R2R3-MYB Genes and Their Contribution to Floral Color in the Perianth-Bearing Piperales

Flavonoids, carotenoids, betalains, and chlorophylls are the plant pigments responsible for floral color. Anthocyanins, a class of flavonoids, are largely responsible for the red, purple, pink, and blue colors. R2R3-MYB genes belonging to subgroup 6 (SG6) are the upstream regulatory factors of the anthocyanin biosynthetic pathway. The canonical members of these genes in Arabidopsis include AtMYB75, AtMYB90, AtMYB113, and AtMYB114. The Aristolochiaceae is an angiosperm lineage with diverse floral groundplans and perianth colors. Saruma henryi exhibits a biseriate perianth with green sepals and yellow petals. All other genera have sepals only, with colors ranging from green (in Lactoris) to a plethora of yellow to red and purple mixtures. Here, we isolated and reconstructed the SG6 R2R3-MYB gene lineage evolution in angiosperms with sampling emphasis in Aristolochiaceae. We found numerous species-specific duplications of this gene lineage in core eudicots and local duplications in Aristolochiaceae for Saruma and Asarum. Expression of SG6 R2R3-MYB genes examined in different developmental stages and plant organs of four Aristolochiaceae species, largely overlaps with red and purple pigments, suggesting a role in anthocyanin and flavonoid synthesis and accumulation. A directed RNA-seq analysis corroborated our RT-PCR analyses, by showing that these structural enzymes activate during perianth development in Aristolochia fimbriata and that the regulatory genes are expressed in correlation with color phenotype. Finally, the reconstruction of the flavonoid and anthocyanin metabolic pathways using predicted peptides from transcriptomic data show that all pivotal enzymes are present in the analyzed species. We conclude that the regulatory genes as well as the biosynthetic pathway are largely conserved across angiosperms. In addition, the Aristolochiaceae emerges as a remarkable group to study the genetic regulatory network for floral color, as their members exhibit an outstanding floral diversity with elaborate color patterns and the genetic complement for SG6 R2R3-MYB genes is simpler than in core eudicot model species.

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.

Overexpression of ThMYB8 mediates salt stress tolerance by directly activating stress-responsive gene expression

MYB transcription factors are important in abiotic stress responses; however, the detailed mechanisms are unclear. Tamarix hispida contains multiple MYB genes. The present study characterized T. hispida MYB8 (ThMYB8) during salt stress using transgenic T. hispida and Arabidopsis assays. ThMYB8 overexpression and ThMYB8 RNAi analysis demonstrated that ThMYB8 enhanced the salt stress tolerance. Transgenic Arabidopsis ectopic expression of ThMYB8 significantly increased root growth, fresh weight, and seed germination rate compared with that of the wild-type under salt stress. Physiological parameters analysis in T. hispida and Arabidopsis showed that ThMYB8 overexpressing plants had the lowest levels of O², H₂O₂, cell death, malondialdehyde, and electrolyte leakage. Overexpression of ThMYB8 regulated Na⁺ and K⁺ concentrations in plant tissues while maintaining K⁺/Na⁺ homeostasis. Analysis using qRT-PCR and ChIP-PCR identified possible downstream ThMYB8-regulated genes. ThMYB8 regulated the expression of ThCYP450-2 (cytochrome p450-2), Thltk (leucine-rich repeat transmembrane protein kinase), and ThTIP (aquaporin TIP) by binding to the MBSI motif ('CAACTG') in their promoters. The results indicated that ThMYB8 enhanced salt stress tolerance in T. hispida by regulating gene expression related to the activation of stress-associated physiological changes, such as enhanced reactive oxygen species scavenging capability, maintaining K⁺/Na⁺ homeostasis, and decreasing the malondialdehyde content and lipid peroxidation cell membranes.

Genome-Wide Identification and Analysis of SRO Gene Family in Chinese Cabbage (Brassica rapa L.)

Similar to radical-induced cell death 1 (SROs) is a family of small proteins unique to plants. SRO transcription factors play an important role in plants’ response to biotic and abiotic stresses. In this study, we identified 12 BrSRO genes in Chinese cabbage (Brassica rapa L.). Among them, a comprehensive overview of the SRO gene family is presented, including physical and chemical characteristics, chromosome locations, phylogenetic analysis, gene structures, motif analysis, and cis-element analyses. The number of amino acids of BrSRO genes is between 77–779 aa, isoelectric point changed from 6.02 to 9.6. Of the 12 BrSRO genes, 11 were randomly distributed along the 7 chromosomes, while BrSRO12 was located along unassigned scaffolds. Phylogenetic analysis indicated that the SRO proteins from six species, including Arabidopsis, banana, rice, Solanum lycopersicum, Zea mays, and Chinese cabbage were divided into eleven groups. The exon-rich BrSRO6 and BrSRO12 containing 15 exons were clustered to group K. All 12 genes have motif 2, which indicate that motif 2 is a relatively conservative motif. There are many hormone and stress response elements in BrSRO genes. The relative expression levels of 12 BrSRO genes under high temperature, drought, salt, and low temperature conditions were analyzed by real-time fluorescence quantitative PCR. The results indicated the relative expression level of BrSRO8 was significantly up-regulated when plants were exposed to high temperature. The relative expression levels of BrSRO1, 3, 7, 8, and 9 were higher under low temperature treatment. The up-regulated genes response to drought and salt stresses were BrSRO1, 5, 9 and BrSRO1, 8, respectively. These results indicated that these genes have certain responses to different abiotic stresses. This work has provided a foundation for further functional analyses of SRO genes in Chinese cabbage.

Phylogenomic Analysis of R2R3 MYB Transcription Factors in Sorghum and their Role in Conditioning Biofuel Syndrome

Background

Large scale cultivation of sorghum for food, feed, and biofuel requires concerted efforts for engineering multipurpose cultivars with optimised agronomic traits. Due to their vital role in regulating the biosynthesis of phenylpropanoid-derived compounds, biomass composition, biotic, and abiotic stress response, R2R3-MYB family transcription factors are ideal targets for improving environmental resilience and economic value of sorghum.

Methods

We used diverse computational biology tools to survey the sorghum genome to identify R2R3-MYB transcription factors followed by their structural and phylogenomic analysis. We used in-house generated as well as publicly available high throughput expression data to analyse the R2R3 expression patterns in various sorghum tissue types.

Results

We have identified a total of 134 R2R3-MYB genes from sorghum and developed a framework to predict gene functions. Collating information from the physical location, duplication, structural analysis, orthologous sequences, phylogeny, and expression patterns revealed the role of duplications in clade-wise expansion of the R2R3-MYB family as well as intra-clade functional diversification. Using publicly available and in-house generated RNA sequencing data, we provide MYB candidates for conditioning biofuel syndrome by engineering phenylpropanoid biosynthesis and sugar signalling pathways in sorghum.

Conclusion

The results presented here are pivotal to prioritize MYB genes for functional validation and optimize agronomic traits in sorghum.

MYB43 in Oilseed Rape (Brassica napus) Positively Regulates Vascular Lignification, Plant Morphology and Yield Potential but Negatively Affects Resistance to Sclerotinia sclerotiorum

Arabidopsis thaliana MYB43 (AtMYB43) is suggested to be involved in cell wall lignification. PtrMYB152, the Populus orthologue of AtMYB43, is a transcriptional activator of lignin biosynthesis and vessel wall deposition. In this research, MYB43 genes from Brassica napus (rapeseed) and its parental species B. rapa and B. oleracea were molecularly characterized, which were dominantly expressed in stem and other vascular organs and showed responsiveness to Sclerotinia sclerotiorum infection. The BnMYB43 family was silenced by RNAi, and the transgenic rapeseed lines showed retardation in growth and development with smaller organs, reduced lodging resistance, fewer silique number and lower yield potential. The thickness of the xylem layer decreased by 28%; the numbers of sclerenchymatous cells, vessels, interfascicular fibers, sieve tubes and pith cells in the whole cross section of the stem decreased by 28%, 59%, 48%, 34% and 21% in these lines, respectively. The contents of cellulose and lignin decreased by 17.49% and 16.21% respectively, while the pectin content increased by 71.92% in stems of RNAi lines. When inoculated with S. sclerotiorum, the lesion length was drastically decreased by 52.10% in the stems of transgenic plants compared with WT, implying great increase in disease resistance. Correspondingly, changes in the gene expression patterns of lignin biosynthesis, cellulose biosynthesis, pectin biosynthesis, cell cycle, SA- and JA-signals, and defensive pathways were in accordance with above phenotypic modifications. These results show that BnMYB43, being a growth-defense trade-off participant, positively regulates vascular lignification, plant morphology and yield potential, but negatively affects resistance to S. sclerotiorum. Moreover, this lignification activator influences cell biogenesis of both lignified and non-lignified tissues of the whole vascular organ.

Root Epidermal Cell Patterning Is Modulated by a Critical Residue in the WEREWOLF Transcription Factor

The Arabidopsis (Arabidopsis thaliana) root epidermis exhibits a position-dependent pattern of root-hair and nonhair cell types. A highly orchestrated network of gene regulatory interactions, including the R2R3-type MYB transcription factor WEREWOLF (WER), is responsible for generating this cell pattern during root development. In this study, we identified a novel wer mutant from a genetic enhancer screen, designated wer-4, that exhibits an abnormal pattern of root-hair and nonhair cells. We established that wer-4 bears a single-residue substitution (D105N) in the DNA-binding R3 MYB repeat of WER, which differentially affects the transcription of WER target genes, including GLABRA2, CAPRICE, TRIPTYCHON, and ENHANCER OF TRY AND CPC1 This modulation of the gene regulatory network leads to altered levels and distributions of cell fate regulators in the differentiating epidermal cells that ultimately generate the abnormal cell-type pattern. We also created several WER variants with substitutions at the Asp-105 position, and these exhibited a variety of gene expression and cell-type pattern alterations, further supporting the critical role of this residue. These findings provide insight into WER protein function and its importance in generating the proper balance of downstream transcriptional factors in the gene regulatory network that establishes root epidermal cell fate.

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.

Mining MYB transcription factors from the genomes of orchids (Phalaenopsis and Dendrobium) and characterization of an orchid R2R3-MYB gene involved in water-soluble polysaccharide biosynthesis

Members of the MYB superfamily act as regulators in a wide range of biological processes in plants. Despite this, the MYB superfamily from the Orchidaceae has not been identified, and MYB genes related to bioactive water-soluble polysaccharide (WSP) biosynthesis are relatively unknown. In this study, we identified 159 and 165 MYB genes from two orchids, Phalaenopsis equestris and Dendrobium officinale, respectively. The MYB proteins were classified into four MYB classes in both orchids: MYB-related (MYBR), R2R3-MYB, 3R-MYB and atypical MYB proteins. The MYBR proteins in both orchids were classified into five subfamilies and 12 genes were strongly up-regulated in response to cold stress in D. officinale. The R2R3-MYB proteins were both divided into 31 clades in P. equestris and D. officinale. Among these clades, nine contained MYB TFs related to secondary cell wall biosynthesis or testa mucilage biosynthesis in Arabidopsis thaliana. In D. officinale, 10 candidate genes showed an expression pattern corresponding to changes in the WSP content. Overexpression of one of these candidate genes (DoMYB75) in A. thaliana increased seed WSP content by about 14%. This study provides information about MYB genes in two orchids that will further help to understand the transcriptional regulation of WSP biosynthesis in these orchids as well as other plant species.

Mining MYB transcription factors from the genomes of orchids (Phalaenopsis and Dendrobium) and characterization of an orchid R2R3-MYB gene involved in water-soluble polysaccharide biosynthesis

Members of the MYB superfamily act as regulators in a wide range of biological processes in plants. Despite this, the MYB superfamily from the Orchidaceae has not been identified, and MYB genes related to bioactive water-soluble polysaccharide (WSP) biosynthesis are relatively unknown. In this study, we identified 159 and 165 MYB genes from two orchids, Phalaenopsis equestris and Dendrobium officinale, respectively. The MYB proteins were classified into four MYB classes in both orchids: MYB-related (MYBR), R2R3-MYB, 3R-MYB and atypical MYB proteins. The MYBR proteins in both orchids were classified into five subfamilies and 12 genes were strongly up-regulated in response to cold stress in D. officinale. The R2R3-MYB proteins were both divided into 31 clades in P. equestris and D. officinale. Among these clades, nine contained MYB TFs related to secondary cell wall biosynthesis or testa mucilage biosynthesis in Arabidopsis thaliana. In D. officinale, 10 candidate genes showed an expression pattern corresponding to changes in the WSP content. Overexpression of one of these candidate genes (DoMYB75) in A. thaliana increased seed WSP content by about 14%. This study provides information about MYB genes in two orchids that will further help to understand the transcriptional regulation of WSP biosynthesis in these orchids as well as other plant species.

Transcriptomic de novo analysis of pitaya (Hylocereus polyrhizus) canker disease caused by Neoscytalidium dimidiatum

BACKGROUND

Canker disease caused by Neoscytalidium dimidiatum is the most serious disease that attacks the pitaya industry. One pathogenic fungus, referred to as ND8, was isolated from the wild-type red-fleshed pitaya (Hylocereus polyrhizus) of Hainan Province. In the early stages of this disease, stems show little spots and a loss of green color. These spots then gradually spread until the stems became rotten due to infection by various strains. Canker disease caused by Neoscytalidium dimidiatum poses a significant threat to pitaya commercial plantations with the growth of stems and the yields, quality of pitaya fruits. However, a lack of transcriptomic and genomic information hinders our understanding of the molecular mechanisms underlying the pitaya defense response.

RESULTS

We investigated the host responses of red-fleshed pitaya (H. polyrhizus) cultivars against N. dimidiatum using Illumina RNA-Seq technology. Significant expression profiles of 23 defense-related genes were further analyzed by qRT-PCR. The total read length based on RNA-Seq was 25,010,007; mean length was 744, the N50 was 1206, and the guanine-cytosine content was 44.48%. Our investigation evaluated 33,584 unigenes, of which 6209 (18.49%) and 27,375 (81.51%) were contigs and singlets, respectively. These unigenes shared a similarity of 16.62% with Vitis vinifera, 7.48% with Theobroma cacao, 6.6% with Nelumbo nucifera and 5.35% with Jatropha curcas. The assembled unigenes were annotated into non-redundant (NR, 25161 unigenes), Kyoto Encyclopedia of Genes and Genomes (KEGG, 17895 unigenes), Clusters of Orthologous Groups (COG, 10475 unigenes), InterPro (19,045 unigenes), and Swiss-Prot public protein databases (16,458 unigenes). In addition, 24 differentially expressed genes, which were mainly associated with plant pathology pathways, were analyzed in-depth.

CONCLUSIONS

This study provides a basis for further in-depth research on the protein function of the annotated unigene assembly with cDNA sequences.

Tomato MYB49 enhances resistance to Phytophthora infestans and tolerance to water deficit and salt stress

MYB49-overexpressing tomato plants showed significant resistance to Phytophthora infestans and tolerance to drought and salt stresses. This finding reveals the potential application of tomato MYB49 in future molecular breeding. Biotic and abiotic stresses severely reduce the productivity of tomato worldwide. Therefore, it is necessary to find key genes to simultaneously improve plant resistance to pathogens and tolerance to various abiotic stresses. In this study, based on homologous relationships with Arabidopsis R2R3-MYBs (AtMYBs) involved in responses to biotic and abiotic stresses, we identified a total of 24 R2R3-MYB transcription factors in the tomato genome. Among these tomato R2R3-MYBs, MYB49 (Solyc10g008700.1) was clustered into subgroup 11 by phylogenetic analysis, and its expression level was significantly induced after treatment with P. infestans, NaCl and PEG6000. Overexpression of MYB49 in tomato significantly enhanced the resistance of tomato to P. infestans, as evidenced by decreases in the number of necrotic cells, sizes of lesion, abundance of P. infestans, and disease index. Likewise, MYB49-overexpressing transgenic tomato plants also displayed increased tolerance to drought and salt stresses. Compared to WT plants, the accumulation of reactive oxygen species (ROS), malonaldehyde content, and relative electrolyte leakage was decreased, and peroxidase activity, superoxide dismutase activity, chlorophyll content, and photosynthetic rate were increased in MYB49-overexpressing tomato plants under P. infestans, salt or drought stress. These results suggested that tomato MYB49, as a positive regulator, could enhance the capacity to scavenge ROS, inhibit cell membrane damage and cell death, and protect chloroplasts, resulting in an improvement in resistance to P. infestans and tolerance to salt and drought stresses, and they provide a candidate gene for tomato breeding to enhance biotic stress resistance and abiotic stress tolerance.

Genome-Wide Characterization and Expression Analyses of Pleurotus ostreatus MYB Transcription Factors during Developmental Stages and under Heat Stress Based on de novo Sequenced Genome

Pleurotus ostreatus is a commercially grown mushroom species in China. However, studies on the mechanisms of the fruiting body development and stress response of P. ostreatus are still at a primary stage. In this study, we report the entire genome sequence of P. ostreatus CCMSSC03989. Then, we performed comprehensive genome-wide characterization and expression analysis of the MYB transcription factor family during a series of developmental stages and under the condition of heat stress. A 34.76 Mb genome was obtained through next-generation sequencing (NGS) and Bionano optical mapping approaches. The genome has a scaffold N50 of 1.1 Mb and contains 10.11% repeats, and 10,936 gene models were predicted. A total of 20 MYB genes (PoMYB) were identified across the genome, and the full-length open reading frames were isolated. The PoMYBs were classified into 1 repeat (1R), 2R, and 3R-MYB groups according to their MYB domain repeat numbers, and 3R-MYBs possessed relatively more introns than 1R and 2R-MYBs. Based on phylogenetic analysis, the PoMYBs were divided into four groups and showed close relationships with the MYB genes of plants and fungi. RNA-sequencing (RNA-Seq) and quantitative PCR (qPCR) analyses revealed that PoMYB expression showed stage-specific patterns in reproductive stages and could be induced by heat stress. The P. ostreatus draft genome will promote genome-wide analysis, and our study of PoMYBs will promote further functional analysis of MYB genes in mushrooms.

The R2R3MYB Gene Family in Phyllostachys edulis: Genome-Wide Analysis and Identification of Stress or Development-Related R2R3MYBs

The MYB transcription factor (TF) is one of the largest gene families in plants and involved to multiple biological processes. However, little is known about the MYB family and its functional role in the genome of moso bamboo. In the present study, a total of 114 R2R3MYB genes were first identified from moso bamboo genome and full-length non-chimeric (FLNC) reads. Phylogenetic analysis coupled with gene structure analysis and motif determination resulted in the division of these PheR2R3MYBs into 17 subgroups. The position of eight proteins along an external branch in the phylogenetic tree suggested their relatively ancient origin. The genes in this group were all substituted by (Met, M)/(Arg, R) at conservative W residues in both R2 and R3 repeats, and half were found to possess no transcriptional activation activity. The analysis of evolutionary patterns and divergence suggests that the expansion of PheMYBs was mainly attributable to whole genome duplication (WGD) under different selection pressures. Expressional analysis based on microarray and qRT-PCR data performed diverse expression patterns of R2R3MYBs in response to both various abiotic stimuli and flower development. Furthermore, the co-expression analysis of R2R3MYBs suggested an intricate interplay of growth- and stress-related responses. Finally, we found a hub gene, PheMYB4, was involved in a complex proteins interaction network. Further functional analysis indicated that ectopic overexpression of its homologous gene, PheMYB4-1, could increase tolerance to cold treatment and sensitivity to drought and salt treatment of transgenic Arabidopsis seedlings. These findings provide comprehensive insights into the MYB family members in moso bamboo and offer candidate MYB genes for further studies on their roles in stress resistance.

SPL13 regulates shoot branching and flowering time in Medicago sativa

Our results show SPL13 plays a crucial role in regulating vegetative and reproductive development in Medicago sativa L. (alfalfa), and that MYB112 is targeted and downregulated by SPL13 in alfalfa. We previously showed that transgenic Medicago sativa (alfalfa) plants overexpressing microRNA156 (miR156) show a bushy phenotype, reduced internodal length, delayed flowering time, and enhanced biomass yield. In alfalfa, transcripts of seven SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors, including SPL13, are targeted for cleavage by miR156. Thus, association of each target SPL gene to a trait or set of traits is essential for developing molecular markers for alfalfa breeding. In this study, we investigated SPL13 function using SPL13 overexpression and silenced alfalfa plants. Severe growth retardation, distorted branches and up-curled leaves were observed in miR156-impervious 35S::SPL13m over-expression plants. In contrast, more lateral branches and delayed flowering time were observed in SPL13 silenced plants. SPL13 transcripts were predominantly present in the plant meristems, indicating that SPL13 is involved in regulating shoot branch development. Accordingly, the shoot branching-related CAROTENOID CLEAVAGE DIOXYGENASE 8 gene was found to be significantly downregulated in SPL13 RNAi silencing plants. A R2R3-MYB gene MYB112 was also identified as being directly silenced by SPL13 based on Next Generation Sequencing-mediated transcriptome analysis and chromatin immunoprecipitation assays, suggesting that MYB112 may be involved in regulating alfalfa vegetative growth.

Genome-wide characterization and expression profiling of NAC transcription factor genes under abiotic stresses in radish (Raphanus sativus L.)

NAC (NAM, no apical meristem; ATAF, Arabidopsis transcription activation factor and CUC, cup-shaped cotyledon) proteins are among the largest transcription factor (TF) families playing fundamental biological processes, including cell expansion and differentiation, and hormone signaling in response to biotic and abiotic stresses. In this study, 172 RsNACs comprising 17 membrane-bound members were identified from the whole radish genome. In total, 98 RsNAC genes were non-uniformly distributed across the nine radish chromosomes. In silico analysis revealed that expression patterns of several NAC genes were tissue-specific such as a preferential expression in roots and leaves. In addition, 21 representative NAC genes were selected to investigate their responses to heavy metals (HMs), salt, heat, drought and abscisic acid (ABA) stresses using real-time polymerase chain reaction (RT-qPCR). As a result, differential expressions among these genes were identified where RsNAC023 and RsNAC080 genes responded positively to all stresses except ABA, while RsNAC145 responded more actively to salt, heat and drought stresses compared with other genes. The results provides more valuable information and robust candidate genes for future functional analysis for improving abiotic stress tolerances in radish.

Functional Characterization of the Versatile MYB Gene Family Uncovered Their Important Roles in Plant Development and Responses to Drought and Waterlogging in Sesame

The MYB gene family constitutes one of the largest transcription factors (TFs) modulating various biological processes in plants. Although genome-wide analysis of this gene family has been carried out in some species, only three MYB members have been functionally characterized heretofore in sesame (Sesamum indicum L.). Here, we identified a relatively high number (287) of sesame MYB genes (SIMYBs) with an uncommon overrepresentation of the 1R-subfamily. A total of 95% of SIMYBs was mapped unevenly onto the 16 linkage groups of the sesame genome with 55 SIMYBs tandemly duplicated. In addition, molecular characterization, gene structure, and evolutionary relationships of SIMYBs were established. Based on the close relationship between sesame and Arabidopsis thaliana, we uncovered that the functions of SIMYBs are highly diverse. A total of 65% of SIMYBs were commonly detected in five tissues, suggesting that they represent key TFs modulating sesame growth and development. Moreover, we found that SIMYBs regulate sesame responses to drought and waterlogging, which highlights the potential of SIMYBs towards improving stress tolerance in sesame. This work presents a comprehensive picture of the MYB gene family in sesame and paves the way for further functional validation of the members of this versatile gene family.

Transcriptome-Wide Identification and Characterization of MYB Transcription Factor Genes in the Laticifer Cells of Hevea brasiliensis

MYB transcription factors hold vital roles in the regulation of plant secondary metabolic pathways. Laticifers in rubber trees (Hevea brasiliensis) are of primary importance in natural rubber production because natural rubber is formed and stored within these structures. To understand the role of MYB transcription factors in the specialized cells, we identified 44 MYB genes (named HblMYB1 to HblMYB44) by using our previously obtained transcriptome database of rubber tree laticifer cells and the public rubber tree genome database. Expression profiles showed that five MYB genes were highly expressed in the laticifers. HblMYB19 and HblMYB44 were selected for further study. HblMYB19 and HblMYB44 bound the promoters of HbFDPS1, HbSRPP, and HRT1 in yeast. Furthermore, the transient overexpression of HblMYB19 and HblMYB44 in tobacco plants significantly increased the activity of the promoters of HbFDPS1, HbSRPP, and HRT1. Basing on this information, we proposed that HblMYB19 and HblMYB44 are the regulators of HbFDPS1, HbSRPP, and HRT1, which are involved in the biosynthesis pathway of natural rubber.

Genome-wide characterization of the WRKY gene family in radish (Raphanus sativus L.) reveals its critical functions under different abiotic stresses

The radish WRKY gene family was genome-widely identified and played critical roles in response to multiple abiotic stresses. The WRKY is among the largest transcription factors (TFs) associated with multiple biological activities for plant survival, including control response mechanisms against abiotic stresses such as heat, salinity, and heavy metals. Radish is an important root vegetable crop and therefore characterization and expression pattern investigation of WRKY transcription factors in radish is imperative. In the present study, 126 putative WRKY genes were retrieved from radish genome database. Protein sequence and annotation scrutiny confirmed that RsWRKY proteins possessed highly conserved domains and zinc finger motif. Based on phylogenetic analysis results, RsWRKYs candidate genes were divided into three groups (Group I, II and III) with the number 31, 74, and 20, respectively. Additionally, gene structure analysis revealed that intron-exon patterns of the WRKY genes are highly conserved in radish. Linkage map analysis indicated that RsWRKY genes were distributed with varying densities over nine linkage groups. Further, RT-qPCR analysis illustrated the significant variation of 36 RsWRKY genes under one or more abiotic stress treatments, implicating that they might be stress-responsive genes. In total, 126 WRKY TFs were identified from the R. sativus genome wherein, 35 of them showed abiotic stress-induced expression patterns. These results provide a genome-wide characterization of RsWRKY TFs and baseline for further functional dissection and molecular evolution investigation, specifically for improving abiotic stress resistances with an ultimate goal of increasing yield and quality of radish.

Genome-wide identification, functional prediction, and evolutionary analysis of the R2R3-MYB superfamily in Brassica napus

R2R3-MYB transcription factors (TFs) have been shown to play important roles in plants, including in development and in various stress conditions. Phylogenetic analysis showed the presence of 249 R2R3-MYB TFs in Brassica napus, called BnaR2R3-MYB TFs, clustered into 38 clades. BnaR2R3-MYB TFs were distributed on 19 chromosomes of B. napus. Sixteen gene clusters were identified. BnaR2R3-MYB TFs were characterized by motif prediction, gene structure analysis, and gene ontology. Evolutionary analysis revealed that BnaR2R3-MYB TFs are mainly formed as a result of whole-genome duplication. Orthologs and paralogs of BnaR2R3-MYB TFs were identified in B. napus, B. rapa, B. oleracea, and Arabidopsis thaliana using synteny-based methods. Purifying selection was pervasive within R2R3-MYB TFs. K_n/K_s values lower than 0.3 indicated that BnaR2R3-MYB TFs are being functionally converged. The role of gene conversion in the formation of BnaR2R3-MYB TFs was significant. Cis-regulatory elements in the upstream regions of BnaR2R3-MYB genes, miRNA targeting BnaR2R3MYB TFs, and post translational modifications were identified. Digital expression data revealed that BnaR2R3-MYB genes were highly expressed in the roots and under high salinity treatment after 24 h. BnaMYB21, BnaMYB141, and BnaMYB148 have been suggested for improving salt-tolerant B. napus. BnaR2R3-MYB genes were mostly up regulated on the 14th day post inoculation with Leptosphaeria biglobosa and L. maculan. BnaMYB150 is a candidate for increased tolerance to Leptospheria in B. napus.

Comparative Transcriptome Analysis Reveals Adaptive Evolution of Notopterygium incisum and Notopterygium franchetii, Two High-Alpine Herbal Species Endemic to China

The extreme conditions (e.g., cold, low oxygen, and strong ultraviolet radiation) of the high mountains provide an ideal natural laboratory for studies on speciation and the adaptive evolution of organisms. Up to now, few genome/transcriptome-based studies have been carried out on how plants adapt to conditions at extremely high altitudes. Notopterygium incisum and Notopterygiumfranchetii (Notopterygium, Apiaceae) are two endangered high-alpine herbal plants endemic to China. To explore the molecular genetic mechanisms of adaptation to high altitudes, we performed high-throughput RNA sequencing (RNA-seq) to characterize the transcriptomes of the two species. In total, more than 130 million sequence reads, 81,446 and 63,153 unigenes with total lengths of 86,924,837 and 62,615,693 bp, were generated for the two herbal species, respectively. OrthoMCL analysis identified 6375 single-copy orthologous genes between N. incisum and N. franchetii. In total, 381 positively-selected candidate genes were identified for both plants by using estimations of the non-synonymous to synonymous substitution rate. At least 18 of these genes potentially participate in RNA splicing, DNA repair, glutathione metabolism and the plant–pathogen interaction pathway, which were further enriched in various functional gene categories possibly responsible for environment adaptation in high mountains. Meanwhile, we detected various transcription factors that regulated the material and energy metabolism in N. incisum and N. franchetii, which probably play vital roles in the tolerance to stress in surroundings. In addition, 60 primer pairs based on orthologous microsatellite-containing sequences between the both Notopterygium species were determined. Finally, 17 polymorphic microsatellite markers (SSR) were successfully characterized for the two endangered species. Based on these candidate orthologous and SSR markers, we detected that the adaptive evolution and species divergence of N. incisum and N. franchetii were significantly associated with the extremely heterogeneous environments and climatic oscillations in high-altitude areas. This work provides important insights into the molecular mechanisms of adaptation to high-altitudes in alpine herbal plants.

Identification and expression analysis under abiotic stress of the R2R3-MYB genes in Ginkgo biloba L

The R2R3-MYB gene family is the largest MYB subfamily in plants and is involved in the regulation of plant secondary metabolism and specific morphogenesis, as well as the response to biotic and abiotic stress. However, a systematic identification and characterization of this gene family has not been carried out in Ginkgo biloba. In this study, we performed a transcriptome-wide survey from four tissues of G. biloba to determine the genetic variation and expression pattern of the R2R3-MYB genes. We analyzed 45 GbMYBs and identified 42 with a complete coding sequence via conserved motif searches. The MYB domain and other motifs in GbMYBs are highly conserved with Arabidopsis thaliana AtMYBs. Phylogenetic analysis of the GbMYBs and AtMYBs categorized the R2R3-MYBs into 26 subgroups, of which 11 subgroups included proteins from both G. biloba and Arabidopsis, and 1 subgroup was specific to G. biloba. Moreover, the GbMYBs expression patterns were analyzed in different tissues and abiotic stress conditions. The results revealed that GbMYBs were differentially expressed in various tissues and following abiotic stresses and phytohormone treatments, indicating their possible roles in biological processes and abiotic stress tolerance and adaptation. Our study demonstrated the functional diversity of the GbMYBs and will provide a foundation for future research into their biological and molecular functions.

Genome-wide characterization and expression analysis enables identification of abiotic stress-responsive MYB transcription factors in cassava (Manihot esculenta)

The myeloblastosis (MYB) transcription factor superfamily is the largest transcription factor family in plants, playing different roles during stress response. However, abiotic stress-responsive MYB transcription factors have not been systematically studied in cassava (Manihot esculenta), an important tropical tuber root crop. In this study, we used a genome-wide transcriptome analysis to predict 299 putative MeMYB genes in the cassava genome. Under drought and cold stresses, many MeMYB genes exhibited different expression patterns in cassava leaves, indicating that these genes might play a role in abiotic stress responses. We found that several stress-responsive MeMYB genes responded to abscisic acid (ABA) in cassava leaves. We characterize four MeMYBs, namely MeMYB1, MeMYB2, MeMYB4, and MeMYB9, as R2R3-MYB transcription factors. Furthermore, RNAi-driven repression of MeMYB2 resulted in drought and cold tolerance in transgenic cassava. Gene expression assays in wild-type and MeMYB2-RNAi cassava plants revealed that MeMYB2 may affect other MeMYBs as well as MeWRKYs under drought and cold stress, suggesting crosstalk between MYB and WRKY family genes under stress conditions in cassava.

Genome-wide characterization and expression analysis of MYB transcription factors in Gossypium hirsutum

Background

MYB family proteins are one of the most abundant transcription factors in the cotton plant and play diverse roles in cotton growth and evolution. Previously, few studies have been conducted in upland cotton, Gossypium hirsutum. The recent release of the G. hirsutum genome sequence provides a great opportunity to identify and characterize the entire upland cotton MYB protein family.

Results

In this study, we undertook a comprehensive genome-wide characterization and expression analysis of the MYB transcription factor family during cotton fiber development. A total of 524 non-redundant cotton MYB genes, among 1986 MYB and MYB-related putative proteins, were identified and classified into four subfamilies including 1R-MYB, 2R-MYB, 3R-MYB, and 4R-MYB. Based on phylogenetic tree analysis, MYB transcription factors were divided into 16 subgroups. The results showed that the majority (69.1 %) of GhMYBs genes belong to the 2R-MYB subfamily in upland cotton.

Conclusion

Our comparative genomics analysis has provided novel insights into the roles of MYB transcription factors in cotton fiber development. These results provide the basis for a greater understanding of MYB regulatory networks and to develop new approaches to improve cotton fiber development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0436-8) contains supplementary material, which is available to authorized users.

Quantitative analyses of the hepatic proteome of methylmercury-exposed Atlantic cod (Gadus morhua) suggest oxidative stress-mediated effects on cellular energy metabolism

BACKGROUND

Methylmecury (MeHg) is a widely distributed environmental pollutant with considerable risk to both human health and wildlife. To gain better insight into the underlying mechanisms of MeHg-mediated toxicity, we have used label-free quantitative mass spectrometry to analyze the liver proteome of Atlantic cod (Gadus morhua) exposed in vivo to MeHg (0, 0.5, 2 mg/kg body weight) for 2 weeks.

RESULTS

Out of a toltal of 1143 proteins quantified, 125 proteins were differentially regulated between MeHg-treated samples and controls. Using various bioinformatics tools, we performed gene ontology, pathway and network enrichment analysis, which indicated that proteins and pathways mainly related to energy metabolism, antioxidant defense, cytoskeleton remodeling, and protein synthesis were regulated in the hepatic proteome after MeHg exposure. Comparison with previous gene expression data strengthened these results, and further supported that MeHg predominantly affects many energy metabolism pathways, presumably through its strong induction of oxidative stress. Some enzymes known to have functionally important oxidation-sensitive cysteine residues in other animals are among the differentially regulated proteins, suggesting their modulations by MeHg-induced oxidative stress. Integrated analysis of the proteomics dataset combined with previous gene expression dataset showed a more pronounced effect of MeHg on amino acid, glucose and fatty acid metabolic pathways, and suggested possible interactions of the cellular energy metabolism and antioxidant defense pathways.

CONCLUSIONS

MeHg disrupts mainly redox homeostasis and energy generating metabolic pathways in cod liver. The energy pathways appear to be modulated through MeHg-induced oxidative stress, possibly mediated by oxidation sensitive enzymes.

Characterization and expression profiling of MYB transcription factors against stresses and during male organ development in Chinese cabbage (Brassica rapa ssp. pekinensis)

MYB proteins comprise a large family of plant transcription factors that play regulatory roles in different biological processes such as plant development, metabolism, and defense responses. To gain insight into this gene superfamily and to elucidate its roles in stress resistance, we performed a comprehensive genome-wide identification, characterization, and expression analysis of MYB genes in Chinese cabbage (Brassica rapa ssp. pekinensis). We identified 475 Chinese cabbage MYB genes, among which most were from R2R3-MYB (256 genes) and MYB-related (202) subfamilies. Analysis of sequence characteristics, phylogenetic classification, and protein motif structures confirmed the existence of several categories (1R, 2R, 3R, 4R, and 5R) of Chinese cabbage MYB genes, which is comparable with MYB genes of other crops. An extensive in silico functional analysis, based on established functional properties of MYB genes from different crop species, revealed 11 and four functional clades within the Chinese cabbage R2R3-MYB and MYB-related subfamilies, respectively. In this study, we reported a MYB-like group within the MYB-related subfamily contains 77 MYB genes. Expression analysis using low temperature-treated whole-genome microarray data revealed variable transcript abundance of 1R/2R/3R/4R/5R-MYB genes in 11 clusters between two inbred lines of Chinese cabbage, Chiifu and Kenshin, which differ in cold tolerance. In further validation tests, we used qRT-PCR to examine the cold-responsive expression patterns of 27 BrMYB genes; surprisingly, the MYB-related genes were induced more highly than the R2R3-MYB genes. In addition, we identified 10 genes with corresponsive expression patterns from a set of salt-, drought-, ABA-, JA-, and SA-induced R2R3-MYB genes. We identified 11 R2R3-MYBs functioning in resistance against biotic stress, including 10 against Fusarium oxysporum f.sp. conglutinans and one against Pectobacterium carotovoram subsp. caratovorum. Furthermore, based on organ-specific expression data, we identified nine R2R3-MYBs that were constitutively expressed in male reproductive tissue, which may provide an important key for studying male sterility in Chinese cabbage. The extensive annotation and transcriptome profiling reported in this study will be useful for understanding the involvement of MYB genes in stress resistance in addition to their growth regulatory functions, ultimately providing the basis for functional characterization and exploitation of the candidate MYB genes for genetic engineering of Chinese cabbage.

Identification of Gene Modules Associated with Low Temperatures Response in Bambara Groundnut by Network-Based Analysis

Bambara groundnut (Vigna subterranea (L.) Verdc.) is an African legume and is a promising underutilized crop with good seed nutritional values. Low temperature stress in a number of African countries at night, such as Botswana, can effect the growth and development of bambara groundnut, leading to losses in potential crop yield. Therefore, in this study we developed a computational pipeline to identify and analyze the genes and gene modules associated with low temperature stress responses in bambara groundnut using the cross-species microarray technique (as bambara groundnut has no microarray chip) coupled with network-based analysis. Analyses of the bambara groundnut transcriptome using cross-species gene expression data resulted in the identification of 375 and 659 differentially expressed genes (p<0.01) under the sub-optimal (23°C) and very sub-optimal (18°C) temperatures, respectively, of which 110 genes are commonly shared between the two stress conditions. The construction of a Highest Reciprocal Rank-based gene co-expression network, followed by its partition using a Heuristic Cluster Chiseling Algorithm resulted in 6 and 7 gene modules in sub-optimal and very sub-optimal temperature stresses being identified, respectively. Modules of sub-optimal temperature stress are principally enriched with carbohydrate and lipid metabolic processes, while most of the modules of very sub-optimal temperature stress are significantly enriched with responses to stimuli and various metabolic processes. Several transcription factors (from MYB, NAC, WRKY, WHIRLY & GATA classes) that may regulate the downstream genes involved in response to stimulus in order for the plant to withstand very sub-optimal temperature stress were highlighted. The identified gene modules could be useful in breeding for low-temperature stress tolerant bambara groundnut varieties.