Genome-wide identification and characterization of R2R3-MYB family in Hypericum perforatum under diverse abiotic stresses

Genome-wide identification of the Pyrus R2R3-MYB gene family and PhMYB62 regulation analysis in Pyrus hopeiensis flowers at low temperature

The R2R3-MYB gene family play an important role in plant growth, development and stress responses. In this study, a total of 122 PcoR2R3-MYB genes were identified and grouped into 26 clades in pear. And these PcoMYBs were unevenly distributed among 17 chromosomes. The sequence characteristics, conversed motifs, exon/intron structures, classification, duplication events and cis-acting elements were also investigated. The gene duplication events showed that segmental duplication may play key roles in expansion of the PcoMYB gene family. Pyrus hopeiensis, which is a valuable wild resource, has strong cold resistance. An integrative analyses of miRNA and mRNA showed that PhMYB62 was involved in regulating low-temperature stress in P. hopeiensis flower organs. Subcellular localization analysis showed that PhMYB62 protein was specifically localized to the nucleus. The result of DAP-seq showed that PhMYB62 responded to low-temperature stress in P. hopeiensis by regulating TFs, which were associated with plant stress resistance, and POD, GAUT12, AUX28 and CHS genes. Subsequently, yeast one-hybrid verified that PhMYB62 could bind and activate the promoter of POD gene. The current study would provide a comprehensive information for further functional research on the stress-responsive R2R3-MYB gene candidates in pear, and may help to identify the genes associated with cold resistance for the cultivation of cold-resistant pear varieties.

Systematic Identification and Functional Analysis of the Hypericum perforatum L. bZIP Gene Family Indicating That Overexpressed HpbZIP69 Enhances Drought Resistance

Basic leucine zipper (bZIP) transcription factors play significant roles in plants' growth and development processes, as well as in response to biological and abiotic stresses. Hypericum perforatum is one of the world's top three best-selling herbal medicines, mainly used to treat depression. However, there has been no systematic identification or functional analysis of the bZIP gene family in H. perforatum. In this study, 79 HpbZIP genes were identified. Based on phylogenetic analysis, the HpbZIP gene family was divided into ten groups, designated A-I and S. The physicochemical properties, gene structures, protein conserved motifs, and Gene Ontology enrichments of all HpbZIPs were systematically analyzed. The expression patterns of all genes in different tissues of H. perforatum (i.e., root, stem, leaf, and flower) were analyzed by qRT-PCR, revealing the different expression patterns of HpbZIP under abiotic stresses. The HpbZIP69 protein is localized in the nucleus. According to the results of the yeast one-hybrid (Y1H) assays, HpbZIP69 can bind to the HpASMT2 (N-acetylserotonin O-methyltransferase) gene promoter (G-box cis-element) to activate its activity. Overexpressing HpbZIP69 in Arabidopsis wild-type lines enhanced their tolerance to drought. The MDA and H2O2 contents were significantly decreased, and the activity of superoxide dismutase (SOD) was considerably increased under the drought stress. These results may aid in additional functional studies of HpbZIP transcription factors, and in cultivating drought-resistant medicinal plants.

Genome-wide identification of Apetala2 gene family in Hypericum perforatum L and expression profiles in response to different abiotic and hormonal treatments

The Apetala2 (AP2) gene family of transcription factors (TFs) play important functions in plant development, hormonal response, and abiotic stress. To reveal the biological functions and the expression profiles of AP2 genes in Hypericum perforatum, genome-wide identification of HpAP2 family members was conducted.

Methods

We identified 21 AP2 TFs in H. perforatum using bioinformatic methods; their physical and chemical properties, gene structures, conserved motifs, evolutionary relationships, cis-acting elements, and expression patterns were investigated.

Results

We found that based on the structural characteristics and evolutionary relationships, the HpAP2 gene family can be divided into three subclasses: euANT, baselANT, and euAP2. A canonical HpAP2 TF shared a conserved protein structure, while a unique motif 6 was found in HpAP2_1, HpAP2_4, and HpAP2_5 from the euANT subgroup, indicating potential biological and regulatory functions of these genes. Furthermore, a total of 59 cis-acting elements were identified, most of which were associated with growth, development, and resistance to stress in plants. Transcriptomics data showed that 57.14% of the genes in the AP2 family were differentially expressed in four organs. For example, HpAP2_18 was specifically expressed in roots and stems, whereas HpAP2_17 and HpAP2_11 were specifically expressed in leaves and flowers, respectively. HpAP2_5, HpAP2_11, and HpAP2_18 showed tissue-specific expression patterns and responded positively to hormones and abiotic stresses.

Conclusion

These results demonstrated that the HpAP2 family genes are involved in diverse developmental processes and generate responses to abiotic stress conditions in H. perforatum. This article, for the first time, reports the identification and expression profiles of the AP2 family genes in H. perforatum, laying the foundation for future functional studies with these genes.

Genomic-wide identification and expression analysis of R2R3-MYB transcription factors related to flavonol biosynthesis in Morinda officinalis

BACKGROUND

The R2R3-MYB transcription factors are a crucial and extensive gene family in plants, which participate in diverse processes, including development, metabolism, defense, differentiation, and stress response. In the Lingnan region of China, Morinda officinalis is extensively grown and is renowned for its use as both a medicinal herb and food source. However, there are relatively few reports on the R2R3-MYB transcription factor family in M.officinalis.

RESULTS

In this study, we identified 97 R2R3-MYB genes in the genome of Morinda officinalis and classified them into 32 subgroups based on phylogenetic comparison with Arabidopsis thaliana. The lack of recent whole-genome duplication events in M.officinalis may be the reason for the relatively few members of the R2R3-MYB family. We also further analyzed the physical and chemical characteristics, conserved motifs, gene structure, and chromosomal location. Gene duplication events found 21 fragment duplication pairs and five tandem duplication event R2R3-MYB genes in M.officinalis may also affect gene family expansion. Based on phylogenetic analysis, cis-element analysis, co-expression analysis and RT-qPCR, we concluded that MoMYB33 might modulate flavonol levels by regulating the expression of 4-coumarate-CoA ligase Mo4CL2, chalcone isomerase MoCHI3, and flavonol synthase MoFLS4/11/12. MoMYB33 and AtMYB111 showed the highest similarity of 79% and may be involved in flavonol synthase networks by the STRING database. Moreover, we also identified MoMYB genes that respond to methyl Jasmonate (MeJA) and abscisic acid (ABA) stress by RT-qPCR.

CONCLUSIONS

This study offers a thorough comprehension of R2R3-MYB in M.officinalis, which lays the foundation for the regulation of flavonol synthesis and the response of MoMYB genes to phytohormones in M.officinalis.

Two R2R3-MYB transcription factors from Chinese cedar (Cryptomeria fortunei Hooibrenk) are involved in the regulation of secondary cell wall formation

As the most abundant renewable energy source, wood comprises the secondary cell wall (SCW). SCW biosynthesis involves lignin and cellulose deposition. Increasing studies have illustrated that R2R3-MYB transcription factors (TFs) play pivotal roles in affecting lignin accumulation and SCW formation. Nevertheless, the regulatory roles of R2R3-MYBs are still unresolved in Cryptomeria fortunei Hooibrenk cambium and wood formation. To dissect the potentials of CfMYBs, we successfully cloned and intensively studied the functions of CfMYB4 and CfMYB5 in SCW formation and abiotic stress response. They both contained the conserved MYB domain capable of forming a special structure that could bind to the core motifs of downstream genes. The phylogenetic tree implied that two CfMYBs clustered into different evolutionary branches. They were predominantly expressed in the stem and were localized to the nucleus. Furthermore, CfMYB4 functioned as an activator to enhance lignin and cellulose accumulation, and increase the SCW thickness by elevating the expression levels of SCW-related genes. By contrast, CfMYB5 negatively regulated lignin and cellulose biosynthesis, and decreased SCW formation by reducing the expression of SCW biosynthetic genes. Our data not only highlight the regulatory functions of CfMYBs in lignin deposition but also provide critical insights into the development of strategies for the genetic improvement of Cryptomeria fortunei wood biomass.

Modulation of Anthocyanin Biosynthesis-Related Genes during the Ripening of Olea europaea L. cvs Carolea and Tondina Drupes in Relation to Environmental Factors

Anthocyanins protect plants against various biotic and abiotic stresses, and anthocyanin-rich foods exert benefits on human health due to their antioxidant activity. Nevertheless, little information is available on the influence of genetic and environmental factors on the anthocyanin content in olive fruits. Based on this consideration, the total anthocyanin content, the genes involved in anthocyanin biosynthesis, and three putative R2R3-MYB transcription factors were evaluated at different ripening stages in the drupes of the Carolea and Tondina cultivars, sampled at different altitudes in the Calabria region, Italy. During drupe ripening, the total anthocyanin content and the transcript levels of analyzed genes gradually increased. In line with the anthocyanin content, a different level of expression of anthocyanin structural genes was observed in 'Carolea' compared to 'Tondina', and in relation to the cultivation area. Furthermore, we identified Oeu050989.1 as a putative R2R3-MYB involved in the regulation of anthocyanin structural genes correlated with the environmental temperature change response. We conclude that anthocyanin accumulation is strongly regulated by development, genotype, and also by environmental factors such as temperature, associated with the altitude gradient. The obtained results contribute to reducing the current information gap regarding the molecular mechanisms on anthocyanin biosynthesis regulation related to the environmental conditions in Olea europaea.

Genome-Wide Characterization of the Mitogen-Activated Protein Kinase Gene Family and Their Expression Patterns in Response to Drought and Colletotrichum Gloeosporioides in Walnut (Juglans regia)

Mitogen-activated protein kinases (MAPKs) are a family of Ser/Thr (serine/threonine) protein kinases that play very important roles in plant responses to biotic and abiotic stressors. However, the MAPK gene family in the important crop walnut (Juglans regia L.) has been less well studied compared with other species. We discovered 25 JrMAPK members in the Juglans genome in this study. The JrMAPK gene family was separated into four subfamilies based on phylogenetic analysis, and members of the same subgroup had similar motifs and exons/introns. A variety of cis-acting elements, mainly related to the light response, growth and development, stress response, and hormone responses, were detected in the JrMAPK gene promoters. Collinearity analysis showed that purification selection was the main driving force in JrMAPK gene evolution, and segmental and tandem duplications played key roles in the expansion of the JrMAPK gene family. The RNA-Seq (RNA Sequencing) results indicated that many of the JrMAPK genes were expressed in response to different levels of Colletotrichum gloeosporioides infection. JrMAPK1, JrMAPK3, JrMAPK4, JrMAPK5, JrMAPK6, JrMAPK7, JrMAPK9, JrMAPK11, JrMAPK12, JrMAPK13, JrMAPK17, JrMAPK19, JrMAPK20, and JrMAPK21 were upregulated at the transcriptional level in response to the drought stress treatment. The results of this study will help in further investigations of the evolutionary history and biological functions of the MAPK gene family in walnut.

Genome-Wide Identification of the Hypericum perforatum WRKY Gene Family Implicates HpWRKY85 in Drought Resistance

WRKY, named for its special heptapeptide conserved sequence WRKYGOK, is one of the largest transcription factor families in plants and is widely involved in plant responses to biotic, abiotic, and hormonal stresses, especially the important regulatory function in response to drought stress. However, there is no complete comprehensive analysis of this family in H. perforatum, which is one of the most extensively studied plants and is probably the best-known herbal medicine on the market today, serving as an antidepressant, neuroprotective, an antineuralgic, and an antiviral. Here, we identified 86 HpWRKY genes according to the whole genome database of H. perforatum, and classified them into three groups through phylogenetic analysis. Gene structure, conserved domain, motif, cis-elements, gene ontology, and expression profiling were performed. Furthermore, it was found that HpWRKY85, a homologous gene of AtWRKY75, showed obvious responses to drought treatment. Subcellular localization analysis indicated that this protein was localized in the nucleus by the Arabidopsis protoplasts transient transfection. Meanwhile, HpWRKY85-overexpressing Arabidopsis plants showed a stronger ability of root growth and scavenging endogenous reactive oxygen species. The results provide a reference for further understanding the role of HpWRKY85 in the molecular mechanism of drought resistance of H. perforatum.

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.

GsMYB7 encoding a R2R3-type MYB transcription factor enhances the tolerance to aluminum stress in soybean (Glycine max L.)

Background

MYB transcription factor (TF) is one of the largest families of TFs in plants and play essential roles in plant growth and development, and is involved in responses to biological and abiotic stress. However, there are few reports on GsMYB7 gene in soybean under aluminum acid stress, and its regulatory mechanism remains unclear.

Results

The GsMYB7 protein is localized in the nucleus and has transcriptional activation ability. Quantitative real-time PCR (qRT-PCR) results showed that GsMYB7 held a constitutive expression pattern rich in roots. When AlCl₃ concentration was 25 µM, the total root surface area (SA) of GsMYB7 transgenic lines were 34.97% higher than that of wild-type Huachun 6 (HC6). While the accumulation of Al³⁺ in root tip of transgenic plants after aluminum treatment was 17.39% lower than that of wild-type. RNA-sequencing analysis indicated that over 1181 genes were regulated by GsMYB7 and aluminum stress. Among all the regulated genes, the expression levels of glutathione peroxidase, protein kinase, cytochrome and other genes in the transgenic lines were significantly higher than those in wild type by acidic aluminum stress. The bioinformatics and qRT-PCR results showed that 9 candidate genes were induced under the treatments of acidic aluminum stress which were indirectly and/or directly regulated by GsMYB7. After AlCl₃ treatments, the transcripts of these genes in GsMYB7 transgenic seedlings were significantly higher than those of wide-type HC6.

Conclusions

The results suggested that GsMYB7 may enhance soybean tolerance to acidic aluminum stress by regulating the downstream genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08744-w.

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

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

Anthocyanin Biosynthesis Associated with Natural Variation in Autumn Leaf Coloration in Quercus aliena Accessions

Quercus aliena is an economically important tree species and one of the dominant native oak species in China. Although its leaves typically turn yellow in autumn, we observed natural variants with red leaves. It is important to understand the mechanisms involved in leaf color variation in this species. Therefore, we compared a Q. aliena tree with yellow leaves and three variants with red leaves at different stages of senescence in order to determine the causes of natural variation. We found that the accumulation of anthocyanins such as cyanidin 3-O-glucoside and cyanidin 3-O-sambubiglycoside had a significant effect on leaf coloration. Gene expression analysis showed upregulation of almost all genes encoding enzymes involved in anthocyanin synthesis in the red-leaved variants during the early and main discoloration stages of senescence. These findings are consistent with the accumulation of anthocyanin in red variants. Furthermore, the variants showed significantly higher expression of transcription factors associated with anthocyanin synthesis, such as those encoded by genes QaMYB1 and QaMYB3. Our findings provide new insights into the physiological and molecular mechanisms involved in autumn leaf coloration in Q. aliena, as well as provide genetic resources for further development and cultivation of valuable ornamental variants of this species.

Identification and Characterization of Jasmonic Acid Biosynthetic Genes in Salvia miltiorrhiza Bunge

Jasmonic acid (JA) is a vital plant hormone that performs a variety of critical functions for plants. Salvia miltiorrhiza Bunge (S. miltiorrhiza), also known as Danshen, is a renowned traditional Chinese medicinal herb. However, no thorough and systematic analysis of JA biosynthesis genes in S. miltiorrhiza exists. Through genome-wide prediction and molecular cloning, 23 candidate genes related to JA biosynthesis were identified in S. miltiorrhiza. These genes belong to four families that encode lipoxygenase (LOX), allene oxide synthase (AOS), allene oxide cyclase (AOC), and 12-OPDA reductase3 (OPR3). It was discovered that the candidate genes for JA synthesis of S. miltiorrhiza were distinct and conserved, in contrast to related genes in other plants, by evaluating their genetic structures, protein characteristics, and phylogenetic trees. These genes displayed tissue-specific expression patterns concerning to methyl jasmonate (MeJA) and wound tests. Overall, the results of this study provide valuable information for elucidating the JA biosynthesis pathway in S. miltiorrhiza by comprehensive and methodical examination.

Systematic Analysis and Functional Characterization of R2R3-MYB Genes in Scutellaria baicalensis Georgi

R2R3-MYB transcription factors participate in multiple critical biological processes, particularly as relates to the regulation of secondary metabolites. The dried root of Scutellaria baicalensis Georgi is a traditional Chinese medicine and possesses various bioactive attributes including anti-inflammation, anti-HIV, and anti-COVID-19 properties due to its flavonoids. In the current study, a total of 95 R2R3-MYB genes were identified in S. baicalensis and classified into 34 subgroups, as supported by similar exon–intron structures and conserved motifs. Among them, 93 R2R3-SbMYBs were mapped onto nine chromosomes. Collinear analysis revealed that segmental duplications were primarily responsible for driving the evolution and expansion of the R2R3-SbMYB gene family. Synteny analyses showed that the ortholog numbers of the R2R3-MYB genes between S. baicalensis and other dicotyledons had a higher proportion compared to that which is found from the monocotyledons. RNA-seq data indicated that the expression patterns of R2R3-SbMYBs in different tissues were different. Quantitative reverse transcriptase-PCR (qRT-PCR) analysis showed that 36 R2R3-SbMYBs from different subgroups exhibited specific expression profiles under various conditions, including hormone stimuli treatments (methyl jasmonate and abscisic acid) and abiotic stresses (drought and cold shock treatments). Further investigation revealed that SbMYB18/32/46/60/70/74 localized in the nucleus, and SbMYB18/32/60/70 possessed transcriptional activation activity, implying their potential roles in the regulatory mechanisms of various biological processes. This study provides a comprehensive understanding of the R2R3-SbMYBs gene family and lays the foundation for further investigation of their biological function.

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.

Structural and interactions analysis of a transcription factor PnMYB2 in Panax notoginseng

The main active ingredients of the traditional Chinese medicinal plant, Panax notoginseng, are the Panax notoginseng saponins (PNS). They can be synthesized via the mevalonate pathway; PnSS and PnSE1 are the key rate-limiting enzymes in this pathway. In this study, an interaction between PnMYB2 and the key enzymes was identified and characterized from the P. notoginseng cDNA library using the Y1H technique. Subsequently, X-α-gal color reaction confirmed the interaction between PnMYB2 and the upstream sequences of PnSS and PnSE1 promoters. Full-length cDNA sequence of PnMYB2 was isolated and characterized. PnMYB2 has an open reading frame of 864 bp, encoding 287 amino acids. 3D structural analysis of PnMYB2 indicated that its structure was similar to that of the template. Phylogenetic analysis revealed that PnMYB2 and PgMYB2 are highly homologous and belong to the R2R3 MYB transcription factor (TF). Subcellular localization analysis showed that PnMYB2 was localized in the nucleus. The recombinant protein PnMYB2 was successfully obtained through prokaryotic expression and was confirmed to be an inclusion body protein. Furthermore, electrophoretic mobility shift assay (EMSA) experiments demonstrated that PnMYB2 specifically binds to MYB core and AC-rich elements. This study provides a theoretical basis for transcriptional regulation of saponin biosynthesis in P. notoginseng.

Convergent gene clusters underpin hyperforin biosynthesis in St John's wort

The meroterpenoid hyperforin is responsible for the antidepressant activity of St John's wort extracts, but the genes controlling its biosynthesis are unknown. Using genome mining and biochemical work, we characterize two biosynthetic gene clusters (BGCs) that encode the first three steps in the biosynthesis of hyperforin precursors. The findings of syntenic and phylogenetic analyses reveal the parallel assembly of the two BGCs. The syntenous BGC in Mesua ferrea indicates that the first cluster was assembled before the divergence of the Hypericaceae and Calophyllaceae families. The assembly of the second cluster is the result of a coalescence of genomic fragments after a major duplication event. The differences between the two BGCs - in terms of gene expression, response to methyl jasmonate, substrate specificity and subcellular localization of key enzymes - suggest that the presence of the two clusters could serve to generate separate pools of precursors. The parallel assembly of two BGCs with similar compositions in a single plant species is uncommon, and our work provides insights into how and when these gene clusters form. Our discovery helps to advance our understanding of the evolution of plant specialized metabolism and its genomic organization. Additionally, our results offer a foundation from which hyperforin biosynthesis can be more fully understood, and which can be used in future metabolic engineering applications.

R2R3-MYB gene family: Genome-wide identification provides insight to improve the content of proanthocyanidins in Trifolium repens

The R2R3-MYB family is one of largest transcription factor families in plants playing significant roles in regulating anthocyanin and proanthocyanidin biosynthesis. Proanthocyanidins are one of major objectives to improve the quality of white clover (Trifolium repens L.), which have a beneficial effect on ruminant to prevent the lethal pasture bloat. A total of 133 TrR2R3-MYB genes were identified and distributed on all 16 chromosomes based on the whole genome information of white clover. Also, by exploring the gene structure, motifs and duplication events of TrR2R3-MYBs, as well as the evolutionary relationship with TrR2R3-MYB genes of other species, 10 TrR2R3-MYB genes with the potential to regulate the anthocyanins and proanthocyanidins biosynthesis were screened. These TrR2R3-MYB genes responded significantly to low temperature in white clover. In addition, they have different expression patterns in leaves, petioles and inflorescences of white clover. Importantly, TrMYB116 and TrMYB118 may positively regulate anthocyanin accumulation and low temperature response in white clover. TrMYB118 may also be associated with anthocyanin pigmentation pattern in Purple leaves. This study provides a basis for verifying the function of TrR2R3-MYB and breeding white clover cultivars with high proanthocyanidins.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

The XTH Gene Family in Schima superba: Genome-Wide Identification, Expression Profiles, and Functional Interaction Network Analysis

Xyloglucan endotransglucosylase/hydrolase (XTH), belonging to glycoside hydrolase family 16, is one of the key enzymes in plant cell wall remodeling. Schima superba is an important timber and fireproof tree species in southern China. However, little is known about XTHs in S. superba. In the present study, a total of 34 SsuXTHs were obtained, which were classified into three subfamilies based on the phylogenetic relationship and unevenly distributed on 18 chromosomes. Furthermore, the intron-exon structure and conserved motif composition of them supported the classification and the members belonging to the same subfamily shared similar gene structures. Segmental and tandem duplication events did not lead to SsuXTH gene family expansion, and strong purifying selection pressures during evolution led to similar structure and function of SsuXTH gene family. The interaction network and cis-acting regulatory elements analysis revealed the SsuXTH expression might be regulated by multiple hormones, abiotic stresses and transcription factors. Finally, expression profiles and GO enrichment analysis showed most of the tandem repeat genes were mainly expressed in the phloem and xylem and they mainly participated in glycoside metabolic processes through the transfer and hydrolysis of xyloglucan in the cell wall and then regulated fiber elongation.

Genome-wide identification of Pistacia R2R3-MYB gene family and function characterization of PcMYB113 during autumn leaf coloration in Pistacia chinensis

Pistacia chinensis is known for its biodiesel production. Several varieties of this plant have leaves that produce anthocyanin, which is responsible for their reddish coloration in autumn. This reddish hue is what makes them useful as ornamental plants. However, the mechanism of anthocyanin accumulation during autumn leaf coloration remains unclear. R2R3-MYB proteins reportedly regulated anthocyanin biosynthesis in many plant species. Here, we performed a genome-wide analysis and expression profiles of R2R3-MYB transcription factor in Pistacia. A total of 158 R2R3-MYB proteins were identified and grouped into 32 clades. Combining the data from RNA-seq and qRT-PCR, one key gene, EVM0016534, was screened and identified to have the highest correlation with anthocyanin accumulation. It was named PcMYB113 due to its sequence similarity to AtMYB113 and it could bind to the promoter of PcF3H. Furthermore, ectopic expression of PcMYB113 in Arabidopsis promoted the accumulation of anthocyanin in the seed coat, cotyledon, and mature leaves, thus confirming the function of PcMYB113 in anthocyanin biosynthesis. In addition, PcMYB113 had a specifically higher expression in senesced red leaves than in mature green leaves and young red leaves in P. chinensis, thereby suggesting the potential role of PcMYB113 in promoting anthocyanin biosynthesis during autumn leaf coloration. These findings enrich our understanding of the function of R2R3-MYB genes in anthocyanin biosynthesis and autumn leaf coloration.

Genome-wide identification, classification and expression profile analysis of the HSF gene family in Hypericum perforatum

Heat shock transcription factors (HSFs) are critical regulators of plant responses to various abiotic and biotic stresses, including high temperature stress. HSFs are involved in regulating the expression of heat shock proteins (HSPs) by binding with heat stress elements (HSEs) to defend against high-temperature stress. The H. perforatum genome was recently fully sequenced; this provides a valuable resource for genetic and functional analysis. In this study, 23 putative HpHSF genes were identified and divided into three groups (A, B, and C) based on phylogeny and structural features. Gene structure and conserved motif analyses were performed on HpHSFs members; the DNA-binding domain (DBD), hydrophobic heptad repeat (HR-A/B), and exon-intron boundaries exhibited specific phylogenetic relationships. In addition, the presence of various cis-acting elements in the promoter regions of HpHSFs underscored their regulatory function in abiotic stress responses. RT-qPCR analyses showed that most HpHSF genes were expressed in response to heat conditions, suggesting that HpHSFs play potential roles in the heat stress resistance pathway. Our findings are advantageous for the analysis and research of the function of HpHSFs in high temperature stress tolerance in H. perforatum.

R2R3-MYB transcription factor family in tea plant (Camellia sinensis): Genome-wide characterization, phylogeny, chromosome location, structure and expression patterns

MYB transcription factors play essential roles in many biological processes and environmental stimuli. However, the functions of the MYB transcription factor family in tea plants have not been elucidated. Here, a total of 122 CsR2R3-MYB genes were identified from the chromosome level genome of tea plant (Camellia sinensis). The CsR2R3-MYB genes were phylogenetically classified into 25 groups. Results from the structure analysis of the gene, conserved motifs, and chromosomal distribution supported the relative conservation of the R2R3-MYB genes family in the tea plant. Synteny analysis indicated that 122, 34, and 112 CsR2R3-MYB genes were orthologous to Arabidopsis thaliana, Oryza sativa and C. sinensis var. 'huangdan' (HD), respectively. Tissue-specific expression showed that all CsR2R3-MYB genes had different expression patterns in the tea plant tissues, indicating that these genes may perform diverse functions. The expression patterns of representative R2R3-MYB genes and the regulatory network of the main anthocyanin components were analyzed, which suggested that CsMYB17 may played a key role in the regulation of cya-3-O-gal, del-3-O-gal, cya-3-O-glu and pel-3-O-glu. Results from the qRT-PCR validation of selected genes suggested that CsR2R3-MYB genes were induced in response to drought, cold, GA, and ABA treatments. Overall, this study provides comprehensive and systematic information for research on the function of R2R3-MYB genes in tea plants.

Whole-genome sequence data of Hypericum perforatum and functional characterization of melatonin biosynthesis by N-acetylserotonin O-methyltransferase

Hypericum perforatum is among the most commonly used herbal remedies and supplements. The aerial plant parts are often used to treat depression. Due to the lack of genomic information of H. perforatum, the gene networks regulating secondary metabolite synthesis remain unclear. Here, we present a high-quality genome for H. perforatum with a 2.3-Mb scaffold N50. The draft assembly covers 91.9% of the predicted genome and represents the fourth sequenced genus in the order Malpighiales. Comparing this sequence with model or related species revealed that Populus trichocarpa and Hevea brasiliensis could be grouped into one branch, while H. perforatum and Linum usitatissimum are grouped in another branch. Combined with transcriptome data, 40 key genes related to melatonin, hyperforin, and hypericin synthesis were screened and analyzed. Five N-acetylserotonin O-methyltransferases (HpASMT1-HpASMT5) were cloned and functionally characterized. Purified HpASMT3 protein converted N-acetylserotonin into melatonin with a V_max of about 1.35 pkat/mg protein. HpASMT1 and HpASMT3 overexpression in Arabidopsis mutants caused 1.5-2-fold higher melatonin content than in mutant and wild-type plants. The endogenous reactive oxygen species (ROS) in transgenic plants was significantly lower than ROS in mutant and wild-type plants, suggesting higher drought tolerance. The obtained genomic data offer new resources for further study on the evolution of Hypericaceae family, but also provide a basis for further study of melatonin biosynthetic pathways in other plants.

Systematic investigation and expression profiles of the GbR2R3-MYB transcription factor family in ginkgo (Ginkgo biloba L.)

As one of the largest families of transcription factors, the R2R3-MYB family plays a significant role in plant growth, development, and response to hormone and environmental stress. To explore its evolutionary mechanism and potential function in Ginkgo biloba, a gymnosperm of great economic and ecological value, we presented a comprehensive analysis of the R2R3-MYB genes in ginkgo. Sixty-nine GbR2R3-MYB genes were identified and these genes could be classified into 33 groups based on the characteristics of the amino acid sequence of the R2R3-MYB domain and gene structure. Syntenic analyses indicated that few tandem and segmental duplications possibly resulted in the contraction of the GbR2R3-MYB gene family. Based on the transcriptome data, expression profiles of eight different tissues and different developmental stages of leaf and kernel showed that GbR2R3-MYB genes had distinct temporal and spatial expression characteristics. Specific expression patterns of the sixteen GbR2R3-MYB genes were also identified in response to different abiotic stresses and hormonal exposures. Further investigation revealed that GbR2R3-MYB19 was located in the nucleus and possessed transcriptional activity, implying its potential roles in the regulation of multiple biological processes. Our findings provide a robust basis for future comprehensive evolutionary and functional analyses of GbR2R3-MYB genes in ginkgo.

Genome-Wide Analysis and the Expression Pattern of the ERF Gene Family in Hypericum perforatum

Hypericum perforatum is a well-known medicinal herb currently used as a remedy for depression as it contains many high levels of secondary metabolites. The ethylene response factor (ERF) family encodes transcriptional regulators with multiple functions that play a vital role in the diverse developmental and physiological processes of plants, which can protect plants from various stresses by regulating the expression of genes. Although the function of several ERF genes from other plants has been further confirmed, H. perforatum is the first sequenced species in Malpighiales, and no information regarding the ERFs has been reported thus far. In this study, a total of 101 ERF genes were identified from H. perforatum. A systematic and thorough bioinformatic analysis of the ERF family was performed using the genomic database of H. perforatum. According to the phylogenetic tree analysis, HpERFs were further classified into 11 subfamilies. Gene ontology (GO) analysis suggested that most of the HpERFs likely participate in the biological processes of plants. The cis-elements were mainly divided into five categories, associated with the regulation of gene transcription, response to various stresses, and plant development. Further analysis of the expression patterns showed that the stress-responsive HpERFs responded to different treatments. This work systematically analyzed HpERFs using the genome sequences of H. perforatum. Our results provide a theoretical basis for further investigation of the function of stress-related ERFs in H. perforatum.

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 Analysis of the Shi-Related Sequence Family and Functional Identification of GmSRS18 Involving in Drought and Salt Stresses in Soybean

The plant-special SHI-RELATED SEQUENCE (SRS) family plays vital roles in various biological processes. However, the genome-wide analysis and abiotic stress-related functions of this family were less reported in soybean. In this work, 21 members of soybean SRS family were identified, which were divided into three groups (Group I, II, and III). The chromosome location and gene structure were analyzed, which indicated that the members in the same group may have similar functions. The analysis of stress-related cis-elements showed that the SRS family may be involved in abiotic stress signaling pathway. The analysis of expression patterns in various tissues demonstrated that SRS family may play crucial roles in special tissue-dependent regulatory networks. The data based on soybean RNA sequencing (RNA-seq) and quantitative Real-Time PCR (qRT-PCR) proved that SRS genes were induced by drought, NaCl, and exogenous abscisic acid (ABA). GmSRS18 significantly induced by drought and NaCl was selected for further functional verification. GmSRS18, encoding a cell nuclear protein, could negatively regulate drought and salt resistance in transgenic Arabidopsis. It can affect stress-related physiological index, including chlorophyll, proline, and relative electrolyte leakage. Additionally, it inhibited the expression levels of stress-related marker genes. Taken together, these results provide valuable information for understanding the classification of soybean SRS transcription factors and indicates that SRS plays important roles in abiotic stress responses.