Genome-wide identification and expression analysis of the MYB transcription factor in moso bamboo (Phyllostachys edulis)

Identification of WRKY gene family in Dioscorea opposita Thunb. reveals that DoWRKY71 enhanced the tolerance to cold and ABA stress

WRKY transcription factors constitute one of the largest plant-specific gene families, regulating various aspects of plant growth, development, physiological processes, and responses to abiotic stresses. This study aimed to comprehensively analyze the WRKY gene family of yam (Dioscorea opposita Thunb.), to understand their expression patterns during the growth and development process and their response to different treatments of yam and analyze the function of DoWRKY71 in detail. A total of 25 DoWRKY genes were identified from the transcriptome of yam, which were divided into six clades (I, IIa, IIc, IId, IIe, III) based on phylogenetic analysis. The analysis of conserved motifs revealed 10 motifs, varying in length from 16 to 50 amino acids. Based on real-time quantitative PCR (qRT-PCR) analysis, DoWRKY genes were expressed at different stages of growth and development and responded differentially to various abiotic stresses. The expression level of DoWRKY71 genes was up-regulated in the early stage and then down-regulated in tuber enlargement. This gene showed responsiveness to cold and abiotic stresses, such as abscisic acid (ABA) and methyl jasmonate (MeJA). Therefore, further study was conducted on this gene. Subcellular localization analysis revealed that the DoWRKY71 protein was localized in the nucleus. Moreover, the overexpression of DoWRKY71 enhanced the cold tolerance of transgenic tobacco and promoted ABA mediated stomatal closure. This study presents the first systematic analysis of the WRKY gene family in yam, offering new insights for studying WRKY transcription factors in yam. The functional study of DoWRKY71 lays theoretical foundation for further exploring the regulatory function of the DoWRKY71 gene in the growth and development related signaling pathway of yam.

Integrated Transcriptome and Metabolome Analyses Reveal Bamboo Culm Color Formation Mechanisms Involved in Anthocyanin Biosynthetic in Phyllostachys nigra

Phyllostachys nigra has green young culms (S1) and purple black mature culms (S4). Anthocyanins are the principal pigment responsible for color presentation in ornamental plants. We employ a multi-omics approach to investigate the regulatory mechanisms of anthocyanins in Ph. nigra. Firstly, we found that the pigments of the culm of Ph. nigra accumulated only in one to four layers of cells below the epidermis. The levels of total anthocyanins and total flavonoids gradually increased during the process of bamboo culm color formation. Metabolomics analysis indicated that the predominant pigment metabolites observed were petunidin 3-O-glucoside and malvidin O-hexoside, exhibiting a significant increase of up to 9.36-fold and 13.23-fold, respectively, during pigmentation of Ph. nigra culm. Transcriptomics sequencing has revealed that genes involved in flavonoid biosynthesis, phenylpropanoid biosynthesis, and starch and sucrose metabolism pathways were significantly enriched, leading to color formation. A total of 62 differentially expressed structural genes associated with anthocyanin synthesis were identified. Notably, PnANS2, PnUFGT2, PnCHI2, and PnCHS1 showed significant correlations with anthocyanin metabolites. Additionally, certain transcription factors such as PnMYB6 and PnMYB1 showed significant positive or negative correlations with anthocyanins. With the accumulation of sucrose, the expression of PnMYB6 is enhanced, which in turn triggers the expression of anthocyanin biosynthesis genes. Based on these findings, we propose that these key genes primarily regulate the anthocyanin synthesis pathway in the culm and contribute to the accumulation of anthocyanin, ultimately resulting in the purple-black coloration of Ph. nigra.

Identification of MYB gene family in medicinal tea tree Melaleuca alternifolia (Maiden and Betche) cheel and analysis of members regulating terpene biosynthesis

BACKGROUND

The tea tree (Melaleuca alternifolia) is renowned for its production of tea tree oil, an essential oil primarily composed of terpenes extracted from its shoot. MYB transcription factors, which are one of the largest TF families, play a crucial role in regulating primary and secondary metabolite synthesis. However, knowledge of the MYB gene family in M. alternifolia is limited.

METHODS AND RESULTS

Here, we conducted a comprehensive genome-wide analysis of MYB genes in M. alternifolia, referred to as MaMYBs, including phylogenetic relationships, structures, promoter regions, and GO annotations. Our findings classified 219 MaMYBs into four subfamilies: one 5R-MYB, four 3R-MYBs, sixty-one MYB-related, and the remaining 153 are all 2R-MYBs. Seven genes (MYB189, MYB146, MYB44, MYB29, MYB175, MYB162, and MYB160) were linked to terpenoid synthesis based on GO annotation. Phylogenetic analysis with Arabidopsis homologous MYB genes suggested that MYB193 and MYB163 may also be involved in terpenoid synthesis. Additionally, through correlation analysis of gene expression and metabolite content, we identified 42 MYB genes associated with metabolite content.

CONCLUSION

The results provide valuable insights into the importance of MYB transcription factors in essential oil production in M. alternifolia. These findings lay the groundwork for a better understanding of the MYB regulatory network and the development of novel strategies to enhance essential oil synthesis in M. alternifolia.

Identification of the MYB gene family in Sorghum bicolor and functional analysis of SbMYBAS1 in response to salt stress

Salt stress adversely affects plant growth and development. It is necessary to understand the underlying salt response mechanism to improve salt tolerance in plants. MYB transcription factors can regulate plant responses to salt stress. However, only a few studies have explored the role of MYB TFs in Sorghum bicolor (L.) Moench. So we decided to make a systematic analysis and research on the sorghum MYB family. A total of 210 MYB genes in sorghum were identified in this study. Furthermore, 210 MYB genes were distributed across ten chromosomes, named SbMYB1-SbMYB210. To study the phylogeny of the identified TFs, 210 MYB genes were divided into six subfamilies. We further demonstrated that SbMYB genes have evolved under strong purifying selection. SbMYBAS1 (SbMYB119) was chosen as the study object, which the expression decreased under salt stress conditions. Further study of the SbMYBAS1 showed that SbMYBAS1 is located in the nucleus. Under salt stress conditions, Arabidopsis plants overexpressed SbMYBAS1 showed significantly lower dry/fresh weight and chlorophyll content but significantly higher membrane permeability, MDA content, and Na+/K+ ratio than the wild-type Arabidopsis plants. Yeast two-hybrid screening result showed that SbMYBAS1 might interact with proteins encoded by SORBI_302G184600, SORBI_3009G247900 and SORBI_3004G59600. Results also showed that SbMYBAS1 could regulate the expression of AtGSTU17, AtGSTU16, AtP5CS2, AtUGT88A1, AtUGT85A2, AtOPR2 and AtPCR2 under salt stress conditions. This work laid a foundation for the study of the response mechanism of sorghum MYB gene family to salt stress.

Comparative analysis of the MYB gene family in seven Ipomoea species

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

Comprehensive analysis of metabolome and transcriptome reveals the mechanism of color formation in different leave of Loropetalum Chinense var. Rubrum

BACKGROUND

Loropetalum chinense var. rubrum (L. chinense var. rubrum) is a precious, coloured-leaf native ornamental plant in the Hunan Province. We found an L. chinense var. rubrum tree with three different leaf colours: GL (green leaf), ML (mosaic leaf), and PL (purple leaf). The mechanism of leaf coloration in this plant is still unclear. Therefore, this study aimed to identify the metabolites and genes involved in determining the colour composition of L. chinense var. rubrum leaves, using phenotypic/anatomic observations, pigment content detection, and comparative metabolomics and transcriptomics.

RESULTS

We observed that the mesophyll cells in PL were purple, while those in GL were green and those in ML were a mix of purple-green. The contents of chlorophyll a, b, carotenoids, and total chlorophyll in PL and ML were significantly lower than those in GL. While the anthocyanin content in PL and ML was significantly higher than that in GL. The metabolomics results showed the differences in the content of cyanidin 3-O-glucoside, delphinidin 3-O-glucoside, cyanidin 3,5-O-diglucoside, pelargonidin, and petunidin 3,5-diglucoside in ML, GL, and PL were significant. Considering that the change trend of anthocyanin content change was consistent with the leaf colour difference, we speculated that these compounds might influence the colour of L. chinense var. rubrum leaves. Using transcriptomics, we finally identified nine differentially expressed structural genes (one ANR (ANR1217); four CYP75As (CYP75A1815, CYP75A2846, CYP75A2909, and CYP75A1716); four UFGTs (UFGT1876, UFGT1649, UFGT1839, and UFGT3273) and nine transcription factors (two MYBs (MYB1057 and MYB1211), one MADS-box (MADS1235), two AP2-likes (AP2-like1779 and AP2-like2234), one bZIP (bZIP3720), two WD40s (WD2173 and WD1867) and one bHLH (bHLH1631) that might be related to flavonoid biosynthesis and then impacted the appearance of colour in L. chinense var. rubrum leaves.

CONCLUSION

This study revealed potential molecular mechanisms associated with leaf coloration in L. chinense var. rubrum by analyzing differential metabolites and genes related to the anthocyanin biosynthesis pathway. It also provided a reference for research on leaf colour variation in other ornamental plants.

Genome-wide survey of HMA gene family and its characterization in wheat (Triticum aestivum)

Background

Abiotic stresses, particularly drought and heavy metal toxicity, have presented a significant risk to long-term agricultural output around the world. Although the heavy-metal-associated domain (HMA) gene family has been widely explored in Arabidopsis and other plants, it has not been thoroughly studied in wheat (Triticum aestivum). This study was proposed to investigate the HMA gene family in wheat.

Methods

To analyze the phylogenetic relationships, gene structure, gene ontology, and conserved motifs, a comparative study of wheat HMA genes with the Arabidopsis genome was performed.

Results

A total of 27 T. aestivum proteins belonging to the HMA gene family were identified in this study, with amino acid counts ranging from 262 to 1,071. HMA proteins were found to be grouped into three subgroups in a phylogenetic tree, and closely related proteins in the tree showed the same expression patterns as motifs found in distinct subgroups. Gene structural study elucidated that intron and exon arrangement differed by family.

Conclusion

As a result, the current work offered important information regarding HMA family genes in the T. aestivum genome, which will be valuable in understanding their putative functions in other wheat species.

Identification of KFB Family in Moso Bamboo Reveals the Potential Function of PeKFB9 Involved in Stress Response and Lignin Polymerization

The Kelch repeat F-box (KFB) protein is an important E3 ubiquitin ligase that has been demonstrated to perform an important post-translational regulatory role in plants by mediating multiple biological processes. Despite their importance, KFBs have not yet been identified and characterized in bamboo. In this study, 19 PeKFBs were identified with F-box and Kelch domains; genes encoding these PeKFBs were unevenly distributed across 12 chromosomes of moso bamboo. Phylogenetic analysis indicated that the PeKFBs were divided into eight subclades based on similar gene structures and highly conserved motifs. A tissue-specific gene expression analysis showed that the PeKFBs were differentially expressed in various tissues of moso bamboo. All the promoters of the PeKFBs contained stress-related cis-elements, which was supported by the differentially expression of PeKFBs of moso bamboo under drought and cold stresses. Sixteen proteins were screened from the moso bamboo shoots' cDNA library using PeKFB9 as a bait through a yeast two-hybrid (Y2H) assay. Moreover, PeKFB9 physically interacted with PeSKP1-like-1 and PePRX72-1, which mediated the activity of peroxidase in proteolytic turnover. Taken together, these findings improved our understanding of PeKFBs, especially in response to stresses, and laid a foundation for revealing the molecular mechanism of PeKFB9 in regulating lignin polymerization by degrading peroxidase.

Multifaceted analyses reveal carbohydrate metabolism mainly affecting the quality of postharvest bamboo shoots

Bamboo shoot is one of nutritious vegetables in China. However, the edible quality of fresh bamboo shoots deteriorates easily after harvest. Here, morphological, physiological, transcriptomic and microRNA sequencing analyses were conducted to investigate the postharvest characteristics of moso bamboo (Phyllostachys edulis) shoots. Rapid decreases of soluble sugars, structural polysaccharides and hydrolyzed tannins, and increases of lignin and condensed tannins were observed in the postharvest bamboo shoots. Differentially expressed genes (DEGs) and miRNAs with opposite trends were mainly enriched in structural polysaccharide metabolism, starch and sucrose metabolism and glycolysis pathways, which were consistent with the changes of carbohydrates. A co-expression network of carbohydrate metabolism was constructed, which was verified by qPCR and yeast one-hybrid (Y1H) assay. Furthermore, the function of one hub glycosyltransferase gene was validated in Arabidopsis, which confirmed that it was involved in xylan biosynthesis. These results are of great significance for revealing the carbohydrate metabolism mechanisms of postharvest bamboo shoots and provide a potential candidate gene for molecular breeding related to xylan in the future.

Genome-wide identification and expression characterization of the DoG gene family of moso bamboo (Phyllostachys edulis)

BACKGROUND

The DoG (Delay of Germination1) family plays a key regulatory role in seed dormancy and germination. However, to date, there is no complete genomic overview of the DoG gene family of any economically valuable crop, including moso bamboo (Phyllostachys edulis), and no studies have been conducted to characterize its expression profile. To identify the DoG gene members of moso bamboo (PeDoG) and to investigate their family structural features and tissue expression profile characteristics, a study was conducted. Based on the whole genome and differential transcriptome data, in this investigation, we have scrutinized the physicochemical properties, gene structure, cis-acting elements, phylogenetic relationships, conserved structural (CS) domains, CS motifs and expression patterns of the PeDoG1 family of moso bamboo.

RESULTS

The DoG family genes of moso bamboo were found distributed across 16 chromosomal scaffolds with 24 members. All members were found to carry DoG1 structural domains, while 23 members additionally possessed basic leucine zipper (bZIP) structural domains. We could divide the PeDoG genes into three subfamilies based on phylogenetic relationships. Covariance analysis revealed that tandem duplication was the main driver of amplification of the PeDoG genes. The upstream promoter of these genes containing several cis-acting elements indicates a plausible role in abiotic stress and hormone induction. Gene expression pattern according to transcriptome data revealed participation of the PeDoG genes in tissue and organ development. Analysis using Short Time-series Expression Miner (STEM) tool revealed that the PeDoG gene family is also associated with rapid early shoot growth. Gene ontology (GO) and KEGG analyses showed a dual role of the PeDoG genes. We found that PeDoGs has a possible role as bZIP transcription factors by regulating Polar like1 (PL1) gene expression, and thereby playing a disease response role in moso bamboo. Quantitative gene expression of the PeDoG genes revealed that they were abundantly expressed in roots and leaves, and could be induced in response to gibberellin (GA).

CONCLUSION

In this study, we found that the PeDoG genes are involved in a wide range of activities such as growth and development, stress response and transcription. This forms the first report of PeDoG genes and their potential roles in moso bamboo.

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.

Genome-wide identification of Cymbidium sinense WRKY gene family and the importance of its Group III members in response to abiotic stress

Transcription factors (TFs) of the WRKY family play pivotal roles in defense responses and secondary metabolism of plants. Although WRKY TFs are well documented in numerous plant species, no study has performed a genome-wide investigation of the WRKY gene family in Cymbidium sinense. In the present work, we found 64 C. sinense WRKY (CsWRKY) TFs, and they were further divided into eight subgroups. Chromosomal distribution of CsWRKYs revealed that the majority of these genes were localized on 16 chromosomes, especially on Chromosome 2. Syntenic analysis implied that 13 (20.31%) genes were derived from segmental duplication events, and 17 orthologous gene pairs were identified between Arabidopsis thaliana WRKY (AtWRKY) and CsWRKY genes. Moreover, 55 of the 64 CsWRKYs were detectable in different plant tissues in response to exposure to plant hormones. Among them, Group III members were strongly induced in response to various hormone treatments, indicating their potential essential roles in hormone signaling. We subsequently analyzed the function of CsWRKY18 in Group III. The CsWRKY18 was localized in the nucleus. The constitutive expression of CsWRKY18 in Arabidopsis led to enhanced sensitivity to ABA-mediated seed germination and root growth and elevated plant tolerance to abiotic stress within the ABA-dependent pathway. Overall, our study represented the first genome-wide characterization and functional analysis of WRKY TFs in C. sinense, which could provide useful clues about the evolution and functional description of CsWRKY genes.

Genome-Wide Identification and Expression Analysis of the Plant U-Box Protein Gene Family in Phyllostachys edulis

The U-box gene encodes a ubiquitin ligase that contains a U-box domain. The plant U-box (PUB) protein plays an important role in the plant stress response; however, very few studies have investigated the role of these proteins in Moso bamboo (Phyllostachys edulis). Thus, more research on PUB proteins is necessary to understand the mechanisms of stress tolerance in P. edulis. In this study, we identified 121 members of the PUB family in P. edulis (PePUB), using bioinformatics based on the P. edulis V2 genome build. The U-box genes of P. edulis showed an uneven distribution among the chromosomes. Phylogenetic analysis of the U-box genes between P. edulis and Arabidopsis thaliana suggested that these genes can be classified into eight subgroups (Groups I–VIII) based on their structural and phylogenetic features. All U-box genes and the structure of their encoded proteins were identified in P. edulis. We further investigated the expression pattern of PePUB genes in different tissues, including the leaves, panicles, rhizomes, roots, and shoots. The qRT-PCR results showed that expression of three genes, PePUB15, PePUB92, and PePUB120, was upregulated at low temperatures compared to that at 25°C. The expression levels of two PePUBs, PePUB60 and PePUB120, were upregulated under drought stress. These results suggest that the PePUB genes play an important role in resistance to low temperatures and drought in P. edulis. This research provides new insight into the function, diversity, and characterization of PUB genes in P. edulis and provides a basis for understanding their biological roles and molecular mechanisms.

Long terminal repeats (LTR) and transcription factors regulate PHRE1 and PHRE2 activity in Moso bamboo under heat stress

BACKGROUND

LTR retrotransposons play a significant role in plant growth, genome evolution, and environmental stress response, but their regulatory response to heat stress remains unclear. We have investigated the activities of two LTR retrotransposons, PHRE1 and PHRE2, of moso bamboo (Phyllostachys edulis) in response to heat stress.

RESULTS

The differential overexpression of PHRE1 and PHRE2 with or without CaMV35s promoter showed enhanced expression under heat stress in transgenic plants. The transcriptional activity studies showed an increase in transposition activity and copy number among moso bamboo wild type and Arabidopsis transgenic plants under heat stress. Comparison of promoter activity in transgenic plants indicated that 5'LTR promoter activity was higher than CaMV35s promoter. Additionally, yeast one-hybrid (Y1H) system and in planta biomolecular fluorescence complementation (BiFC) assay revealed interactions of heat-dependent transcription factors (TFs) with 5'LTR sequence and direct interactions of TFs with pol and gag.

CONCLUSIONS

Our results conclude that the 5'LTR acts as a promoter and could regulate the LTR retrotransposons in moso bamboo under heat stress.

Identification and expression analysis of the glycosyltransferase GT43 family members in bamboo reveal their potential function in xylan biosynthesis during rapid growth

Background

Xylan is one of the most abundant hemicelluloses and can crosslink cellulose and lignin to increase the stability of cell walls. A number of genes encoding glycosyltransferases play vital roles in xylan biosynthesis in plants, such as those of the GT43 family. However, little is known about glycosyltransferases in bamboo, especially woody bamboo which is a good substitute for timber.

Results

A total of 17 GT43 genes (PeGT43–1 ~ PeGT43–17) were identified in the genome of moso bamboo (Phyllostachys edulis), which belong to three subfamilies with specific motifs. The phylogenetic and collinearity analyses showed that PeGT43s may have undergone gene duplication, as a result of collinearity found in 12 pairs of PeGT43s, and between 17 PeGT43s and 10 OsGT43s. A set of cis-acting elements such as hormones, abiotic stress response and MYB binding elements were found in the promoter of PeGT43s. PeGT43s were expressed differently in 26 tissues, among which the highest expression level was found in the shoots, especially in the rapid elongation zone and nodes. The genes coexpressed with PeGT43s were annotated as associated with polysaccharide metabolism and cell wall biosynthesis. qRT–PCR results showed that the coexpressed genes had similar expression patterns with a significant increase in 4.0 m shoots and a peak in 6.0 m shoots during fast growth. In addition, the xylan content and structural polysaccharide staining intensity in bamboo shoots showed a strong positive correlation with the expression of PeGT43s. Yeast one-hybrid assays demonstrated that PeMYB35 could recognize the 5′ UTR/promoter of PeGT43–5 by binding to the SMRE cis-elements.

Conclusions

PeGT43s were found to be adapted to the requirement of xylan biosynthesis during rapid cell elongation and cell wall accumulation, as evidenced by the expression profile of PeGT43s and the rate of xylan accumulation in bamboo shoots. Yeast one-hybrid analysis suggested that PeMYB35 might be involved in xylan biosynthesis by regulating the expression of PeGT43–5 by binding to its 5′ UTR/promoter. Our study provides a comprehensive understanding of PeGT43s in moso bamboo and lays a foundation for further functional analysis of PeGT43s for xylan biosynthesis during rapid growth.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08192-y.

Genome-Wide Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Petal Pigment Synthesis in Liriodendron

The MYB transcription factor family is one of the largest families in plants, and its members have various biological functions. R2R3-MYB genes are involved in the synthesis of pigments that yield petal colors. Liriodendron plants are widely cultivated as ornamental trees owing to their peculiar leaves, tulip-like flowers, and colorful petals. However, the mechanism underlying petal coloring in this species is unknown, and minimal information about MYB genes in Liriodendron is available. Herein, this study aimed to discern gene(s) involved in petal coloration in Liriodendron via genome-wide identification, HPLC, and RT-qPCR assays. In total, 204 LcMYB superfamily genes were identified in the Liriodendron chinense genome, and 85 R2R3-MYB genes were mapped onto 19 chromosomes. Chromosome 4 contained the most (10) R2R3-MYB genes, and chromosomes 14 and 16 contained the fewest (only one). MEME analysis showed that R2R3-MYB proteins in L. chinense were highly conserved and that their exon-intron structures varied. The HPLC results showed that three major carotenoids were uniformly distributed in the petals of L. chinense, while lycopene and β-carotene were concentrated in the orange band region in the petals of Liriodendron tulipifera. Furthermore, the expression profiles via RT-qPCR assays revealed that four R2R3-MYB genes were expressed at the highest levels at the S3P/S4P stage in L. tulipifera. This result combined with the HPLC results showed that these four R2R3-MYB genes might participate in carotenoid synthesis in the petals of L. tulipifera. This work laid a cornerstone for further functional characterization of R2R3-MYB genes in Liriodendron plants.

A regulatory network driving shoot lignification in rapidly growing bamboo

Woody bamboo is environmentally friendly, abundant, and an alternative to conventional timber. Degree of lignification and lignin content and deposition affect timber properties. However, the lignification regulatory network in monocots is poorly understood. To elucidate the regulatory mechanism of lignification in moso bamboo (Phyllostachys edulis), we conducted integrated analyses using transcriptome, small RNA, and degradome sequencing followed by experimental verification. The lignification degree and lignin content increased with increased bamboo shoot height, whereas phenylalanine ammonia-lyase and Laccase activities first increased and then decreased with shoot growth. Moreover, we identified 11,504 differentially expressed genes (DEGs) in different portions of the 13th internodes of different height shoots; most DEGs associated with cell wall and lignin biosynthesis were upregulated, whereas some DEGs related to cell growth were downregulated. We identified a total of 1,502 miRNAs, of which 687 were differentially expressed. Additionally, in silico and degradome analyses indicated that 5,756 genes were targeted by 691 miRNAs. We constructed a regulatory network of lignification, including 11 miRNAs, 22 transcription factors, and 36 enzyme genes, in moso bamboo. Furthermore, PeLAC20 overexpression increased lignin content in transgenic Arabidopsis (Arabidopsis thaliana) plants. Finally, we proposed a reliable miRNA-mediated "MYB-PeLAC20" module for lignin monomer polymerization. Our findings provide definite insights into the genetic regulation of bamboo lignification. In addition to providing a platform for understanding related mechanisms in other monocots, these insights could be used to develop strategies to improve bamboo timber properties.

Comparative genomic analysis of MYB transcription factors for cuticular wax biosynthesis and drought stress tolerance in Helianthus annuus L

Sunflower is an important oil-seed crop in Pakistan, it is mainly cultivated in the spring season. It is severely affected by drought stress resulting in lower yield. Cuticular wax acts as the first defense line to protect plants from drought stress condition. It seals the aerial parts of plants and reduce the water loss from leaf surfaces. Various myeloblastosis (MYB) transcription factors (TFs) are involved in biosynthesis of epicuticular waxes under drought-stress. However, less information is available for MYB, TFs in drought stress and wax biosynthesis in sunflower. We used different computational tools to compare the Arabidopsis MYB, TFs involved in cuticular wax biosynthesis and drought stress tolerance with sunflower genome. We identified three putative MYB genes (MYB16, MYB94 and MYB96) in sunflower along with their seven homologs in Arabidopsis. Phylogenetic association of MYB TFs in Arabidopsis and sunflower indicated strong conservation of TFs in plant species. From gene structure analysis, it was observed that intron and exon organization was family-specific. MYB TFs were unevenly distributed on sunflower chromosomes. Evolutionary analysis indicated the segmental duplication of the MYB gene family in sunflower. Quantitative Real-Time PCR revealed the up-regulation of three MYB genes under drought stress. The gene expression of MYB16, MYB94 and MYB96 were found many folds higher in experimental plants than control. The present study provided the first insight into MYB TFs family's characterization in sunflower under drought stress conditions and wax biosynthesis TFs.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-wide identification and expression analysis of dirigent-jacalin genes from plant chimeric lectins in Moso bamboo (Phyllostachys edulis)

Dirigent-jacalin (D-J) genes belong to the plant chimeric lectin family, and play vital roles in plant growth and resistance to abiotic and biotic stresses. To explore the functions of the D-J family in the growth and development of Moso bamboo (Phyllostachys edulis), their physicochemical properties, phylogenetic relationships, gene and protein structures, and expression patterns were analyzed in detail. Four putative PeD-J genes were identified in the Moso bamboo genome, and microsynteny and phylogenetic analyses indicated that they represent a new branch in the evolution of plant lectins. PeD-J proteins were found to be composed of a dirigent domain and a jacalin-related lectin domain, each of which contained two different motifs. Multiple sequence alignment and homologous modeling analysis indicated that the three-dimensional structure of the PeD-J proteins was significantly different compared to other plant lectins, primarily due to the tandem dirigent and jacalin domains. We surveyed the upstream putative promoter regions of the PeD-Js and found that they mainly contained cis-acting elements related to hormone and abiotic stress response. An analysis of the expression patterns of root, leaf, rhizome and panicle revealed that four PeD-J genes were highly expressed in the panicle, indicating that they may be required during the formation and development of several different tissue types in Moso bamboo. Moreover, PeD-J genes were shown to be involved in the rapid growth and development of bamboo shoots. Quantitative Real-time PCR (qRT PCR) assays further verified that D-J family genes were responsive to hormones and stresses. The results of this study will help to elucidate the biological functions of PeD-Js during bamboo growth, development and stress response.

Genome-Wide Analysis Reveals the Potential Role of MYB Transcription Factors in Floral Scent Formation in Hedychium coronarium

The MYB gene family is one of the largest groups of transcription factors (TFs) playing diverse roles in several biological processes. Hedychium coronarium (white ginger lily) is a renowned ornamental plant both in tropical and subtropical regions due to its flower shape and strong floral scent mainly composed of terpenes and benzenoids. However, there is no information available regarding the role of the MYB gene family in H. coronarium. In the current study, the MYB gene family was identified and extensively analyzed. The identified 253 HcMYB genes were unevenly mapped on 17 chromosomes at a different density. Promoter sequence analysis showed numerous phytohormones related to cis-regulatory elements. The majority of HcMYB genes contain two to three introns and motif composition analysis showed their functional conservation. Phylogenetic analysis revealed that HcMYBs could be classified into 15 distinct clades, and the segmental duplication events played an essential role in the expansion of the HcMYB gene family. Tissue-specific expression patterns of HcMYB genes displayed spatial and temporal expression. Furthermore, seven HcMYB (HcMYB7/8/75/79/145/238/248) were selected for further investigation. Through RT-qPCR, the response of candidates HcMYB genes toward jasmonic acid methyl ester (MeJA), abscisic acid (ABA), ethylene, and auxin was examined. Yeast one-hybrid (Y1H) assays revealed that candidate genes directly bind to the promoter of bottom structural volatile synthesis genes (HcTPS1, HcTPS3, HcTPS10, and HcBSMT2). Moreover, yeast two-hybrid (Y2H) assay showed that HcMYB7/8/75/145/248 interact with HcJAZ1 protein. In HcMYB7/8/79/145/248-silenced flowers, the floral volatile contents were decreased and downregulated the expression of key structural genes, suggesting that these genes might play crucial roles in floral scent formation in H. coronarium by regulating the expression of floral scent biosynthesis genes. Collectively, these findings indicate that HcMYB genes might be involved in the regulatory mechanism of terpenoids and benzenoid biosynthesis in H. coronarium.

Identification of Isoflavonoid Biosynthesis-Related R2R3-MYB Transcription Factors in Callerya speciosa (Champ. ex Benth.) Schot Using Transcriptome-Based Gene Coexpression Analysis

The R2R3-MYB family is one of the largest plant transcription factor (TF) families playing vital roles in defense, plant growth, and secondary metabolism biosynthesis. Although this gene family has been studied in many species, isoflavonoid biosynthesis-related R2R3-MYB TFs in Callerya speciosa (Champ. ex Benth.) Schot, a traditional Chinese medicinal herb, are poorly understood. Here, a total of 101 R2R3-MYB TFs were identified from C. speciosa transcriptome dataset. 25 clades divided into five functional groups were clustered based on the sequence similarity and phylogenetic tree. Conserved motifs and domain distribution, expression patterns, and coexpression networks were also employed to identify the potential R2R3-MYB TFs in the regulation of isoflavonoid biosynthesis. In silico evaluation showed that the deduced R2R3-CsMYB proteins contain highly conserved R2R3 repeat domain at the N-terminal region, that is the signature motif of R2R3-type MYB TFs. Eight potential TFs (CsMYB17, CsMYB36, CsMYB41, CsMYB44, CsMYB45, CsMYB46, CsMYB72, and CsMYB81) had high degrees of coexpression with four key isoflavonoid biosynthetic genes (CsIFS, CsCHS7, CsHID-1, and CsCHI3), in which CsMYB36 as a potential regulator possessed the highest degree. HPLC analysis showed that formononetin and maackiain contents were significantly increased during the development of tuberous roots, which might be controlled by both related R2R3-CsMYBs and structural genes involved in the isoflavonoid biosynthesis pathway. The transcriptome data were further validated by reverse transcription real-time PCR (RT-qPCR) analysis, and similar expression profiles between TFs and key structural genes were identified. This study was the first step toward the understanding of the R2R3-MYB TFs regulating isoflavonoid biosynthesis in C. speciosa. The results will provide information for further functional analysis and quality improvement through genetic manipulation of these potential R2R3-CsMYB genes in C. speciosa.

Genome-wide identification of R2R3-MYB family in wheat and functional characteristics of the abiotic stress responsive gene TaMYB344

Background

MYB superfamily is one of the most abundant families in plants, which plays important roles in plant growth, development, and productivity. However, to date, researches on MYBs in wheat (Triticum aestivum L.) are scattered mostly, not comprehensive.

Results

In this study, a total of 393 R2R3-MYBs and 12 R1R2R3-MYBs were identified and analyzed including gene structure, chromosomal distribution, synteny relationship, and evolutionary relationship. Then, 29 clusters tandem duplication and 8 clusters segmental duplication genes were discovered. The expression profile of the identified genes under abiotic and biotic stress was analyzed using RNA-seq data. Based on expression patterns analysis, we screened many candidate genes involved in plant response to abiotic and biotic stress. Among them, the functional characteristics of TaMYB344 were further studied. TaMYB344 was localized in the nucleus and functioned as a weak transcriptional activator. We demonstrated that TaMYB344-overexpressing transgenic tobacco plants had enhanced tolerance to drought, heat, and high salt stress.

Conclusions

In this study, 393 R2R3-MYBs and 12 R1R2R3-MYBs in wheat were systemically identified and analyzed. Differential expression analysis indicated that many R2R3-MYBs were involved in abiotic and biotic stress response. We identified a potential candidate gene TaMYB344, overexpression of which in tobacco plants enhanced drought, heat, and salt stress tolerance. These results will provide abundant molecular data for breeding new varieties of wheat in the future.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07175-9.

Genome-wide identification and expression analysis of the MYB transcription factor in Japanese plum (Prunus salicina)

MYB proteins constitute one of the largest transcription factor families in plants, members of which are involved in various plant physiological and biochemical processes. Japanese plum (Prunus salicina) is one of the important stone fruit crops worldwide. To date, no comprehensive study of the MYB family in Japanese plum has been reported. In this study, we performed genome-wide analysis of MYB genes in Japanese plum including the phylogeny, gene structures, protein motifs, chromosomal locations, collinearity and expression patterns analysis. A total of 96 Japanese plum R2R3-MYB (PsMYB) genes were characterized and distributed on 8 chromosomes at various densities. Collinearity analysis indicated that the segmental duplication events played a crucial role in the expansion of PsMYB genes, and the interspecies synteny analysis revealed the orthologous gene pairs between Japanese plum and other four selected Rosaceae species. The 96 PsMYB genes could be classified into 27 subgroups based on phylogenetic topology, as supported by the conserved gene structures and motif compositions. Further comparative phylogenetic analysis revealed the functional divergence of MYB gene family during evolution, and three subgroups which included only Rasaceae MYB genes were identified. Expression analysis revealed the distinct expression profiles of the PsMYB genes, and further functional predictions found some of them might be associated with the plum fruit quality traits. Our researches provide a global insight into the organization, phylogeny, evolution and expression patterns of the PsMYB genes, and contribute to the greater understanding of their functional roles in Japanese plum.

Genome-Wide Characterization and Comparative Analysis of MYB Transcription Factors in Ganoderma Species

Numerous studies in plants have shown the vital roles of MYB transcription factors in signal transduction, developmental regulation, biotic/abiotic stress responses and secondary metabolism regulation. However, less is known about the functions of MYBs in Ganoderma. In this study, five medicinal macrofungi of genus Ganoderma were subjected to a genome-wide comparative analysis of MYB genes. A total of 75 MYB genes were identified and classified into four types: 1R-MYBs (52), 2R-MYBs (19), 3R-MYBs (2) and 4R-MYBs (2). Gene structure analysis revealed varying exon numbers (3-14) and intron lengths (7-1058 bp), and noncanonical GC-AG introns were detected in G. lucidum and G. sinense. In a phylogenetic analysis, 69 out of 75 MYB genes were clustered into 15 subgroups, and both single-copy orthologous genes and duplicated genes were identified. The promoters of the MYB genes harbored multiple cis-elements, and specific genes were co-expressed with the G. lucidum MYB genes, indicating the potential roles of these MYB genes in stress response, development and metabolism. This comprehensive and systematic study of MYB family members provides a reference and solid foundation for further functional analysis of MYB genes in Ganoderma species.

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.

Insights into the Diversification and Evolution of R2R3-MYB Transcription Factors in Plants

As one of the largest families of transcription factors (TFs) in plants, R2R3-MYB proteins play crucial roles in regulating a series of plant-specific biological processes. Although the diversity of plant R2R3-MYB TFs has been studied previously, the processes and mechanisms underlying the expansion of these proteins remain unclear. Here, we performed evolutionary analyses of plant R2R3-MYB TFs with dense coverage of streptophyte algae and embryophytes. Our analyses revealed that ancestral land plants exhibited 10 subfamilies of R2R3-MYB proteins, among which orthologs of seven subfamilies were present in chlorophytes and charophycean algae. We found that asymmetric gene duplication events in different subfamilies account for the expansion of R2R3-MYB proteins in embryophytes. We further discovered that the largest subfamily of R2R3-MYBs in land plants, subfamily VIII, emerged in the common ancestor of Zygnematophyceae and embryophytes. During plant terrestrialization, six duplication events gave rise to seven clades of subfamily VIII. Subsequently, this TF subfamily showed a tendency for expansion in bryophytes, lycophytes, and ferns and extensively diversified in ancestral gymnosperms and angiosperms in clades VIII-A-1, VIII-D, and VIII-E. In contrast to subfamily VIII, other subfamilies of R2R3-MYB TFs have remained less expanded across embryophytes. The findings regarding phylogenetic analyses, auxiliary motifs, and DNA-binding specificities provide insight into the evolutionary history of plant R2R3-MYB TFs and shed light on the mechanisms underlying the extensive expansion and subsequent sub- and neofunctionalization of these proteins.

Identification of WRKY gene family and characterization of cold stress-responsive WRKY genes in eggplant

BACKGROUND

WRKY proteins play a vital role in the plants response to different stresses, growth and development. Studies of WRKY proteins have been mainly focused on model plant Arabidopsis and a few other vegetable plants. However, the systematical study of eggplant WRKY transcription factor superfamily is scarce.

METHODS

Bioinformatics has been used to identify and characterize the eggplant WRKY gene family. For the exploration of the differentially expressed WRKY genes, two cultivars with different cold-tolerance were used. Finally, we performed a virus-induced gene silencing (VIGS) experiment to verify the functions of SmWRKY26 and SmWRKY32.

RESULTS

Fifty eight (58) genes encoding eggplant WRKY proteins were identified through searching the eggplant genome. Eggplant WRKY proteins could be classified into three groups or seven subgroups in accordance with other plants. WRKY variants were identified from the eggplant. Gene structure analysis showed that the number of intron in eggplant WRKY family was from 0 to 11, with an average of 4.4. Conserved motif analysis suggested that WRKY DNA-binding domain was conserved in eggplant WRKY proteins. Furthermore, RNA-seq data showed that WRKY genes were differentially expressed in eggplant response to cold stress. By using VIGS, the two differentially expressed genes-SmWRKY26 and SmWRKY32 were verified in response to cold stress.

DISCUSSIONS

This study provides a foundation for further exploring the functions of WRKY proteins in eggplant response to stresses and eggplant genetic improvement in stresses.

Identification of Homeobox Genes Associated with Lignification and Their Expression Patterns in Bamboo Shoots

: Homeobox (HB) genes play critical roles in regulating various aspects of plant growth and development. However, little is known about HB genes in bamboo. In this study, a total of 115 HB genes (PeHB001‒PeHB115) were identified from moso bamboo (Phyllostachys edulis) and grouped into 13 distinct classes (BEL, DDT, HD-ZIP I‒IV, KNOX, NDX, PHD, PINTOX, PLINC, SAWADEE, and WOX) based on the conserved domains and phylogenetic analysis. The number of members in the different classes ranged from 2 to 24, and they usually varied in terms of exon‒intron distribution pattern and length. There were 20 conserved motifs found in 115 PeHBs, with motif 1 being the most common. Gene ontology (GO) analysis showed that PeHBs had diverse molecular functions, with 19 PeHBs being annotated as having xylem development, xylem, and phloem pattern formation functions. Co-expression network analysis showed that 10 of the 19 PeHBs had co-expression correlations, and three members of the KNOX class were hub proteins that interacted with other transcription factors (TFs) such as MYB, bHLH, and OVATE, which were associated with lignin synthesis. Yeast two-hybridization results further proved that PeHB037 (BEL class) interacted with PeHB057 (KNOX class). Transcriptome expression profiling indicated that all PeHBs except PeHB017 were expressed in at least one of the seven tissues of moso bamboo, and 90 PeHBs were expressed in all the tissues. The qRT-PCR results of the 19 PeHBs showed that most of them were upregulated in shoots as the height increased. Moreover, a KNOX binding site was found in the promoters of the key genes involved in lignin synthesis such as Pe4CL, PeC3H, PeCCR, and PeCOMT, which had positive expression correlations with five KNOX genes. Similar results were found in winter bamboo shoots with prolonged storage time, which was consistent with the degree of lignification. These results provide basic data on PeHBs in moso bamboo, which will be helpful for future functional research on PeHBs with positive regulatory roles in the process of lignification.

Genome-Wide Investigation of the NAC Gene Family and Its Potential Association with the Secondary Cell Wall in Moso Bamboo

NAC (NAM, ATAF, and CUC) transcription factors (TFs) are implicated in the transcriptional regulation of diverse processes and have been characterized in a number of plant species. However, NAC TFs are still not well understood in bamboo, especially their potential association with the secondary cell wall (SCW). Here, 94 PeNACs were identified and characterized in moso bamboo (Phyllostachys edulis). Based on their gene structures and conserved motifs, the PeNACs were divided into 11 groups according to their homologs in Arabidopsis. PeNACs were expressed variously in different tissues of moso bamboo, suggesting their functional diversity. Fifteen PeNACs associated with the SCW were selected for co-expression analysis and validation. It was predicted that 396 genes were co-expressed with the 15 PeNACs, in which 16 and 55 genes were involved in the lignin catabolic process and cellulose biosynthetic process respectively. As the degree of lignification in the growing bamboo shoots increased, all 15 PeNACs were upregulated with a trend of rising first and then decreasing except PeNAC37, which increased continuously. These results indicated that these PeNACs might play important roles in SCW biosynthesis and lignification in bamboo shoots. Seven of 15 PeNACs had been found positively co-expressed with seven PeMYBs, and they had similar expression patterns with those of the PeMYBs in bamboo shoots. The targeted sites of miR164 were found in 16 PeNACs, of which three PeNACs associated with SCW were validated to have an opposite expression trend to that of miR164 in growing bamboo shoots. In addition, three PeNACs were selected and verified to have self-activation activities. These results provide comprehensive information of the NAC gene family in moso bamboo, which will be helpful for further functional studies of PeNACs to reveal the molecular regulatory mechanisms of bamboo wood property.

Genome-wide identification of growth-regulating factors in moso bamboo (Phyllostachys edulis): in silico and experimental analyses

Growth-regulating factor (GRF), a small plant-specific transcription factor (TF) family, is extensively involved in the regulation of growth and developmental processes. However, the GRF family has not been comprehensively studied in moso bamboo (Phyllostachys edulis), a typical non-timber forest member. Here, 18 GRF genes were identified and characterized from the moso bamboo genome, and they clustered into three subfamilies (A, B and C). PeGRF genes were analyzed to determine their gene structures, conserved motifs and promoter. The non-synonymous/synonymous substitution ratios of paralogous and orthologous were less than 1, indicating that the GRF family mainly experienced purifying selection during evolution. According to the analysis of tissue-specific expression patterns, the participation of moso bamboo GRFs might be required during the formation and development of these five tissues. Moreover, PeGRF proteins might be involved in the regulation of plant development in biological processes. The qRT-PCR analysis demonstrated that PeGRF genes played essential roles in combating hormonal stresses and they might be involved in hormone regulation. PeGRF11, a nuclear localized protein as assessed by a subcellular localization assay, could interact with PeGIF3 in yeast and in planta according to yeast two-hybridization and bimolecular fluorescence complementation assays (BiFC) assays. But PeGRF11, as a TF, had no transcriptional activity in yeast. These results provide useful information for future functional research on the GRF genes in moso bamboo.