The R2R3-MYB gene family in banana (Musa acuminata): Genome-wide identification, classification and expression patterns

Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.

Genome-wide identification and expression analysis of the MADS gene family in sweet orange (Citrus sinensis) infested with pathogenic bacteria

The risk of pathogenic bacterial invasion in plantations has increased dramatically due to high environmental climate change and has seriously affected sweet orange fruit quality. MADS genes allow plants to develop increased resistance, but functional genes for resistance associated with pathogen invasion have rarely been reported. MADS gene expression profiles were analyzed in sweet orange leaves and fruits infested with Lecanicillium psalliotae and Penicillium digitatum, respectively. Eighty-two MADS genes were identified from the sweet orange genome, and they were classified into five prime subfamilies concerning the Arabidopsis MADS gene family, of which the MIKC subfamily could be subdivided into 13 minor subfamilies. Protein structure analysis showed that more than 93% of the MADS protein sequences of the same subfamily between sweet orange and Arabidopsis were very similar in tertiary structure, with only CsMADS8 and AG showing significant differences. The variability of MADS genes protein structures between sweet orange and Arabidopsis subgroups was less than the variabilities of protein structures within species. Chromosomal localization and covariance analysis showed that these genes were unevenly distributed on nine chromosomes, with the most genes on chromosome 9 and the least on chromosome 2, with 36 and two, respectively. Four pairs of tandem and 28 fragmented duplicated genes in the 82 MADS gene sequences were found in sweet oranges. GO (Gene Ontology) functional enrichment and expression pattern analysis showed that the functional gene CsMADS46 was strongly downregulated of sweet orange in response to biotic stress adversity. It is also the first report that plants' MADS genes are involved in the biotic stress responses of sweet oranges. For the first time, L. psalliotae was experimentally confirmed to be the causal agent of sweet orange leaf spot disease, which provides a reference for the research and control of pathogenic L. psalliotae.

Genome-wide identification, phylogeny and expression analysis of the R2R3-MYB gene family in quinoa (Chenopodium quinoa) under abiotic stress

The MYB transcription factor (TF) are among the largest gene families of plants being responsible for several biological processes. The R2R3-MYB gene family are integral player regulating plant primary and secondary metabolism, growth and development, and responses to hormones and stresses. The phylogenetic analysis combined with gene structure analysis and motif determination resulted in division of R2R3-MYB gene family into 27 subgroups. Evidence generated from synteny analyses indicated that CqR2R3-MYBs gene family is featured by tandem and segmental duplication events. On the basis of RNA-Seq data, the expression patterns of different tissues under salt treatment were investigated resulting CqR2R3-MYB genes high expression both in roots and stem of quinoa (Chenopodium quinoa ) plants. More than half of CqR2R3-MYB genes showed expression under salt stress. Based on this result, CqR2R3-MYB s may regulate quinoa plant growth development and resistance to abiotic stresses. These findings provided comprehensive insights on role of CqR2R3-MYBs gene family members in quinoa and candidate MYB gene family members can be further studies on their role for abiotic stress tolerance in crop plants.

Petunia as a model for MYB transcription factor action under salt stress

Salinity is a current and growing problem, affecting crops worldwide by reducing yields and product quality. Plants have different mechanisms to adapt to salinity; some crops are highly studied, and their salinity tolerance mechanisms are widely known. However, there are other crops with commercial importance that still need characterization of their molecular mechanisms. Usually, transcription factors are in charge of the regulation of complex processes such as the response to salinity. MYB-TFs are a family of transcription factors that regulate various processes in plant development, and both central and specialized metabolism. MYB-TFs have been studied extensively as mediators of specialized metabolism, and some are master regulators. The influence of MYB-TFs on highly orchestrated mechanisms, such as salinity tolerance, is an attractive research target. The versatility of petunia as a model species has allowed for advances to be made in multiple fields: metabolomic pathways, quality traits, stress resistance, and signal transduction. It has the potential to be the link between horticultural crops and lab models, making it useful in translating discoveries related to the MYB-TF pathways into other crops. We present a phylogenetic tree made with Petunia axillaris and Petunia inflata R2R3-MYB subfamily sequences, which could be used to find functional conservation between different species. This work could set the foundations to improve salinity resistance in other commercial crops in later studies.

Mutation of MYB36 affects isoflavonoid metabolism, growth, and stress responses in Lotus japonicus

Isoflavonoids are mostly produced by legumes although little is known about why and how legumes are able to regulate the biosynthesis of these particular compounds. Understanding the role of potential regulatory genes of the isoflavonoid biosynthetic pathway constitutes an important topic of research. The LORE1 mutation of the gene encoding the transcription factor MYB36 allowed the identification of this gene as a regulator of isoflavonoid biosynthesis in Lotus japonicus plants. The levels of several isoflavonoid compounds were considerably lower in two lines of Ljmyb36 mutant plants compared to the WT. In addition, we found that Ljmyb36 mutant plants were significantly smaller and showed a substantial decrease in the chlorophyll levels under normal growth conditions. The analysis of plants subjected to different types of abiotic stress conditions further revealed that mutant plants presented a higher sensitivity than WT plants, indicating that the MYB36 transcription factor is also involved in the stress response in L. japonicus plants.

Genome-wide analysis of the R2R3-MYB gene family in Spatholobus suberectus and identification of its function in flavonoid biosynthesis

Spatholobus suberectus Dunn (S. suberectus), a plant species within the Leguminosae family, has a long history of use in traditional medicines. The dried stem of S. suberectus exhibits various pharmacological activities because it contains various flavonoids. Diverse functions in plants are associated with the R2R3-MYB gene family, including the biosynthesis of flavonoids. Nonetheless, its role remains unelucidated in S. suberectus. Therefore, the newly sequenced S. suberectus genome was utilized to conduct a systematic genome-wide analysis of the R2R3-MYB gene family. The resulting data identified 181 R2R3-SsMYB genes in total, which were then categorized by phylogenetic analysis into 35 subgroups. Among the R2R3-SsMYB genes, 174 were mapped to 9 different chromosomes, and 7 genes were not located on any chromosome. Moreover, similarity in terms of exon-intron structures and motifs was exhibited by most genes in the same subgroup. The expansion of the gene family was primarily driven by segmental duplication events, as demonstrated by collinearity analysis. Notably, most of the duplicated genes underwent purifying selection, which was depicted through the Ka/Ks analysis. In this study, 22 R2R3-SsMYB genes were shown to strongly influence the level of flavonoids. The elevated expression level of these genes was depicted in the tissues with flavonoid accumulation in contrast with other tissues through qRT-PCR data. The resulting data elucidate the structural and functional elements of R2R3-SsMYB genes and present genes that could potentially be utilized for enhancing flavonoid biosynthesis in S. suberectus.

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.

Three R2R3-MYB transcription factors from banana (Musa acuminata) activate structural anthocyanin biosynthesis genes as part of an MBW complex

OBJECTIVE

Bananas are one of the most popular fruits in the world, providing food security and employment opportunities in several developing countries. Increasing the anthocyanin content of banana fruit could improve the health-promoting properties. Anthocyanin biosynthesis is largely regulated at the transcriptional level. However, relatively little is known about the transcriptional activation of anthocyanin biosynthesis in banana.

RESULTS

We analysed the regulatory activity of three Musa acuminata MYBs that were predicted by bioinformatic analysis to transcriptionally regulate anthocyanin biosynthesis in banana. MaMYBA1, MaMYBA2 and MaMYBPA2 did not complement the anthocyanin-deficient phenotype of the Arabidopsis thaliana pap1/pap2 mutant. However, co-transfection experiments in A. thaliana protoplasts showed that MaMYBA1, MaMYBA2 and MaMYBPA2 function as components of a transcription factor complex with a bHLH and WD40 protein, the so called MBW complex, resulting in the activation of the A. thaliana ANTHOCYANIDIN SYNTHASE and DIHYDROFLAVONOL 4-REDUCTASE promoters. The activation potential of MaMYBA1, MaMYBA2 and MaMYBPA2 was increased when combined with the monocot Zea mays bHLH ZmR instead of the dicot AtEGL3. This work paves the path towards decoding the MBW complex-mediated transcriptional activation of anthocyanin biosynthesis in banana. It will also facilitate research towards increased anthocyanin content in banana and other monocot crops.

Genome-wide analysis of R2R3-MYB transcription factors reveals their differential responses to drought stress and ABA treatment in desert poplar (Populus euphratica)

The R2R3-MYB transcription factors are widely involved in the regulation of plant growth, biotic and abiotic stress responses. Meanwhile, seed germination, which is stimulated by internal and external environments, is a critical stage in the plant life cycle. However, the identification, characterization, and expression profiling of the Populus euphratica R2R3-MYB family in drought response during seed germination have been unknown. Our study attempted to identify and characterize the R2R3-MYB genes in P. euphratica (PeR2R3-MYBs) and explore how R2R3-MYBs trigger the drought and abscisic acid (ABA) response mechanism in its seedlings. Based on the analysis of comparative genomics, 174 PeR2R3-MYBs were identified and expanded driven by whole genome duplication or segment duplication events. The analysis of Ka/Ks ratios showed that, in contrast to most PeR2R3-MYBs, the other PeR2R3-MYBs were subjected to positive selection in P. euphratica. Further, the expression data of PeR2R3-MYBs under drought stress and ABA treatment, together with available functional data for Arabidopsis thaliana MYB genes, supported the hypothesis that PeR2R3-MYBs involved in response to drought are dependent or independent on ABA signaling pathway during seed germination, especially PeR2R3-MYBs with MYB binding sites (MBS) cis-element and/or tandem duplication. This study is the first report on the genome-wide analysis of PeR2R3-MYBs, as well as the other two Salicaceae species. The duplication events and differential expressions of PeR2R3-MYBs play important roles in enhancing the adaptation to drought desert environment. Our results provide a reference for prospective functional studies of R2R3-MYBs of poplars and lay the foundation for new breeding strategies to improve the drought tolerance of P. euphratica.

Genome-wide identification, characterization, and expression profile ofNBS-LRRgene family in sweet orange (Citrussinensis)

BACKGROUND

The NBS-LRR (nucleotide-binding site-leucine-rich repeat gene) gene family, known as the plant R (resistance) gene family with the most members, plays a significant role in plant resistance to various external adversity stresses. The NBS-LRR gene family has been researched in many plant species. Citrus is one of the most vital global cash crops, the number one fruit group, and the third most traded agricultural product world wild. However, as one of the largest citrus species, a comprehensive study of the NBS-LRR gene family has not been reported on sweet oranges.

METHODS

In this study, NBS-LRR genes were identified from the Citrus sinensis genome (v3.0), with a comprehensive analysis of this gene family performed, including phylogenetic analysis, gene structure, cis-acting element of a promoter, and chromosomal localization, among others. The expression pattern of NBS-LRR genes was analyzed when sweet orange fruits were infected by Penicillium digitatum, employing experimental data from our research group. It first reported the expression patterns of NBS-LRR genes under abiotic stresses, using three transcript data from NCBI (National Center for Biotechnology Information).

RESULTS

In this study, 111 NBS-LRR genes were identified in the C. sinensis genome (v3.0) and classified into seven subfamilies according to their N-terminal and C-terminal domains. The phylogenetic tree results indicate that genes containing only the NBS structural domain are more ancient in the sweet orange NBS-LRR gene family. The chromosome localization results showed that 111 NBS-LRR genes were distributed unevenly on nine chromosomes, with the most genes distributed on chromosome 1. In addition, we identified a total of 18 tandem duplication gene pairs in the sweet orange NBS-LRR gene family, and based on the Ka/Ks ratio, all of the tandem duplication genes underwent purifying selection. Transcriptome data analysis showed a significant number of NBS-LRR genes expressed under biotic and abiotic stresses, and some reached significantly different levels of expression. It indicates that the NBS-LRR gene family is vital in resistance to biotic and abiotic stresses in sweet oranges.

CONCLUSION

Our study provides the first comprehensive framework on the NBS-LRR family of genes, which provides a basis for further in-depth studies on the biological functions of NBS-LRR in growth, development, and response to abiotic stresses in sweet orange.

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

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

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

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

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

Interplay between R2R3 MYB-type activators and repressors regulates proanthocyanidin biosynthesis in banana (Musa acuminata)

Proanthocyanidins are oligomeric flavonoids that promote plant disease resistance and benefit human health. Banana is one of the world's most extensively farmed crops and its fruit pulp contain proanthocyanidins. However, the transcriptional regulatory network that fine tunes proanthocyanidin biosynthesis in banana remains poorly understood. We characterised two proanthocyanidin-specific R2R3 MYB activators (MaMYBPA1-MaMYBPA2) and four repressors (MaMYBPR1-MaMYBPR4) to elucidate the mechanisms underlying the transcriptional regulation of proanthocyanidin biosynthesis in banana. Heterologous expression of MaMYBPA1 and MaMYBPA2 partially complemented the Arabidopsis thaliana proanthocyanidin-deficient transparent testa2 mutant. MaMYBPA1 and MaMYBPA2 interacted physically with MaMYCs to transactivate anthocyanin synthase, leucoanthocyanidin reductase, and anthocyanidin reductase genes in vitro and form functional MYB-bHLH-WD Repeat (MBW) complexes with MaTTG1 to transactivate these promoters in vivo. Overexpression of MaMYBPAs alone or with MaMYC in banana fruits induced proanthocyanidin accumulation and transcription of proanthocyanidin biosynthesis-related genes. MaMYBPR repressors are also shown to interact with MaMYCs forming repressing MBW complexes, and diminished proanthocyanidin accumulation. Interestingly overexpression of MaMYBPA induces the expression of MaMYBPR, indicating an agile regulation of proanthocyanidin biosynthesis through the formation of competitive MBW complexes. Our results reveal regulatory modules of R2R3 MYB- that fine tune proanthocyanidin biosynthesis and offer possible targets for genetic manipulation for nutritional improvement of banana.

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

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

The R2R3-MYB gene family in Cicer arietinum: genome-wide identification and expression analysis leads to functional characterization of proanthocyanidin biosynthesis regulators in the seed coat

We identified 119 typical CaMYB encoding genes and reveal the major components of the proanthocyanidin regulatory network. CaPARs emerged as promising targets for genetic engineering toward improved agronomic traits in C. arietinum. Chickpea (Cicer arietinum) is among the eight oldest crops and has two main types, i.e., desi and kabuli, whose most obvious difference is the color of their seeds. We show that this color difference is due to differences in proanthocyanidin content of seed coats. Using a targeted approach, we performed in silico analysis, metabolite profiling, molecular, genetic, and biochemical studies to decipher the transcriptional regulatory network involved in proanthocyanidin biosynthesis in the seed coat of C. arietinum. Based on the annotated C. arietinum reference genome sequence, we identified 119 typical CaMYB encoding genes, grouped in 32 distinct clades. Two CaR2R3-MYB transcription factors, named CaPAR1 and CaPAR2, clustering with known proanthocyanidin regulators (PARs) were identified and further analyzed. The expression of CaPAR genes correlated well with the expression of the key structural proanthocyanidin biosynthesis genes CaANR and CaLAR and with proanthocyanidin levels. Protein-protein interaction studies suggest the in vivo interaction of CaPAR1 and CaPAR2 with the bHLH-type transcription factor CaTT8. Co-transfection analyses using Arabidopsis thaliana protoplasts showed that the CaPAR proteins form a MBW complex with CaTT8 and CaTTG1, able to activate the promoters of CaANR and CaLAR in planta. Finally, transgenic expression of CaPARs in the proanthocyanidin-deficient A. thaliana mutant tt2-1 leads to complementation of the transparent testa phenotype. Taken together, our results reveal main components of the proanthocyanidin regulatory network in C. arietinum and suggest that CaPARs are relevant targets of genetic engineering toward improved agronomic traits.

Integrative Analysis of Metabolome and Transcriptome Identifies Potential Genes Involved in the Flavonoid Biosynthesis in Entada phaseoloides Stem

Entada phaseoloides stem is known for its high medicinal benefits and ornamental value. Flavonoids are one of the main active constituents in E. phaseoloides stem. However, the regulatory mechanism of flavonoids accumulation in E. phaseoloides is lacking. Here, phytochemical compounds and transcripts from stems at different developmental stages in E. phaseoloides were investigated by metabolome and transcriptome analysis. The metabolite profiling of the oldest stem was obviously different from young and older stem tissues. A total of 198 flavonoids were detected, and flavones, flavonols, anthocyanins, isoflavones, and flavanones were the main subclasses. The metabolome data showed that the content of acacetin was significantly higher in the young stem and older stem than the oldest stem. Rutin and myricitrin showed significantly higher levels in the oldest stem. A total of 143 MYBs and 143 bHLHs were identified and classified in the RNA-seq data. Meanwhile, 34 flavonoid biosynthesis structural genes were identified. Based on the expression pattern of structural genes involved in flavonoid biosynthesis, it indicated that flavonol, anthocyanin, and proanthocyanin biosynthesis were first active during the development of E. phaseoloides stem, and the anthocyanin or proanthocyanin biosynthesis branch was dominant; the flavone biosynthesis branch was active at the late developmental stage of the stem. Through the correlation analysis of transcriptome and metabolome data, the potential candidate genes related to regulating flavonoid synthesis and transport were identified. Among them, the MYBs, bHLH, and TTG1 are coregulated biosynthesis of flavonols and structural genes, bHLH and transporter genes are coregulated biosynthesis of anthocyanins. In addition, the WDR gene TTG1-like (AN11) may regulate dihydrochalcones and flavonol biosynthesis in specific combinations with IIIb bHLH and R2R3-MYB proteins. Furthermore, the transport gene protein TRANSPARENT TESTA 12-like gene is positively regulated the accumulation of rutin, and the homolog of ABC transporter B family member gene is positively correlated with the content of flavone acacetin. This study offered candidate genes involved in flavonoid biosynthesis, information of flavonoid composition and characteristics of flavonoids accumulation, improved our understanding of the MYBs and bHLHs-related regulation networks of flavonoid biosynthesis in E. phaseoloides stem, and provided references for the metabolic engineering of flavonoid biosynthesis in E. phaseoloides stem.

Rosa1, a Transposable Element-Like Insertion, Produces Red Petal Coloration in Rose Through Altering RcMYB114 Transcription

Rose (Rosa sp.) flowers have a rich diversity of colors resulting from the differential accumulation of anthocyanins, flavonols, and carotenoids. However, the genetic and molecular determinants of the red-petal trait in roses remains poorly understood. Here we report that a transposable element-like insertion (Rosa1) into RcMYB114, a R2R3-MYB transcription factor's promoter region causes its transcription, resulting in red petals. In red-petal varieties, RcMYB114 is expressed specifically in flower organs, but is absent from non-red varieties. Sequencing, yeast two-hybrid, transient transformation, and promoter activity assays of RcMYB114 independently confirmed the role of Rosa1 in altering RcMYB114's transcription and downstream effects on flower color. Genetic and molecular evidence confirmed that the Rosa1 transposable element-like insertion, which is a previously unknown DNA transposable element, is different from those in other plants and is a reliable molecular marker to screen red-petal roses.

Automatic identification and annotation of MYB gene family members in plants

BACKGROUND

MYBs are among the largest transcription factor families in plants. Consequently, members of this family are involved in a plethora of processes including development and specialized metabolism. The MYB families of many plant species were investigated in the last two decades since the first investigation looked at Arabidopsis thaliana. This body of knowledge and characterized sequences provide the basis for the identification, classification, and functional annotation of candidate sequences in new genome and transcriptome assemblies.

RESULTS

A pipeline for the automatic identification and functional annotation of MYBs in a given sequence data set was implemented in Python. MYB candidates are identified, screened for the presence of a MYB domain and other motifs, and finally placed in a phylogenetic context with well characterized sequences. In addition to technical benchmarking based on existing annotation, the transcriptome assembly of Croton tiglium and the annotated genome sequence of Castanea crenata were screened for MYBs. Results of both analyses are presented in this study to illustrate the potential of this application. The analysis of one species takes only a few minutes depending on the number of predicted sequences and the size of the MYB gene family. This pipeline, the required bait sequences, and reference sequences for a classification are freely available on github: https://github.com/bpucker/MYB_annotator .

CONCLUSIONS

This automatic annotation of the MYB gene family in novel assemblies makes genome-wide investigations consistent and paves the way for comparative studies in the future. Candidate genes for in-depth analyses are presented based on their orthology to previously characterized sequences which allows the functional annotation of the newly identified MYBs with high confidence. The identification of orthologs can also be harnessed to detect duplication and deletion events.

Genome-Wide Identification of R2R3-MYB Transcription Factors: Discovery of a "Dual-Function" Regulator of Gypenoside and Flavonol Biosynthesis in Gynostemma pentaphyllum

The R2R3-MYB gene family participates in several plant physiological processes, especially the regulation of the biosynthesis of secondary metabolites. However, little is known about the functions of R2R3-MYB genes in Gynostemma pentaphyllum (G. pentaphyllum), a traditional Chinese medicinal herb that is an excellent source of gypenosides (a class of triterpenoid saponins) and flavonoids. In this study, a systematic genome-wide analysis of the R2R3-MYB gene family was performed using the recently sequenced G. pentaphyllum genome. In total, 87 R2R3-GpMYB genes were identified and subsequently divided into 32 subgroups based on phylogenetic analysis. The analysis was based on conserved exon-intron structures and motif compositions within the same subgroup. Collinearity analysis demonstrated that segmental duplication events were majorly responsible for the expansion of the R2R3-GpMYB gene family, and Ka/Ks analysis indicated that the majority of the duplicated R2R3-GpMYB genes underwent purifying selection. A combination of transcriptome analysis and quantitative reverse transcriptase-PCR (qRT-PCR) confirmed that Gynostemma pentaphyllum myeloblastosis 81 (GpMYB81) along with genes encoding gypenoside and flavonol biosynthetic enzymes exhibited similar expression patterns in different tissues and responses to methyl jasmonate (MeJA). Moreover, GpMYB81 could bind to the promoters of Gynostemma pentaphyllum farnesyl pyrophosphate synthase 1 (GpFPS1) and Gynostemma pentaphyllum chalcone synthase (GpCHS), the key structural genes of gypenoside and flavonol biosynthesis, respectively, and activate their expression. Altogether, this study highlights a novel transcriptional regulatory mechanism that suggests that GpMYB81 acts as a "dual-function" regulator of gypenoside and flavonol biosynthesis in G. pentaphyllum.

Genome-Wide Identification and Expression Analysis of the R2R3-MYB Transcription Factor Family Revealed Their Potential Roles in the Flowering Process in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is a productive fruit crop with high nutritional and medical value in tropical and subtropical regions. The MYB gene family is one of the most widespread plant transcription factor (TF) families participating in the flowering regulation. However, little is known about the MYB TFs involved in the flowering process in longan and its regulatory network. In this study, a total of 119 DlR2R3-MYB genes were identified in the longan genome and were phylogenetically grouped into 28 subgroups. The groupings were supported by highly conserved gene structures and motif composition of DlR2R3-MYB genes in each subgroup. Collinearity analysis demonstrated that segmental replications played a more crucial role in the expansion of the DlR2R3-MYB gene family compared to tandem duplications, and all tandem/segmental duplication gene pairs have evolved under purifying selection. Interspecies synteny analysis among longan and five representative species implied the occurrence of gene duplication events was one of the reasons contributing to functional differentiation among species. RNA-seq data from various tissues showed DlR2R3-MYB genes displayed tissue-preferential expression patterns. The pathway of flower development was enriched with six DlR2R3-MYB genes. Cis-acting element prediction revealed the putative functions of DlR2R3-MYB genes were related to the plant development, phytohormones, and environmental stresses. Notably, the orthologous counterparts between Arabidopsis and longan R2R3-MYB members tended to play conserved roles in the flowering regulation and stress responses. Transcriptome profiling on off-season flower induction (FI) by KClO₃ indicated two up-regulated and four down-regulated DlR2R3-MYB genes involved in the response to KClO₃ treatment compared with control groups. Additionally, qRT-PCR confirmed certain genes exhibited high expression in flowers/flower buds. Subcellular localization experiments revealed that three predicted flowering-associated MYB proteins were localized in the nucleus. Future functional studies on these potential candidate genes involved in the flowering development could further the understanding of the flowering regulation mechanism.

Genomic Approaches for Improvement of Tropical Fruits: Fruit Quality, Shelf Life and Nutrient Content

The breeding of tropical fruit trees for improving fruit traits is complicated, due to the long juvenile phase, generation cycle, parthenocarpy, polyploidy, polyembryony, heterozygosity and biotic and abiotic factors, as well as a lack of good genomic resources. Many molecular techniques have recently evolved to assist and hasten conventional breeding efforts. Molecular markers linked to fruit development and fruit quality traits such as fruit shape, size, texture, aroma, peel and pulp colour were identified in tropical fruit crops, facilitating Marker-assisted breeding (MAB). An increase in the availability of genome sequences of tropical fruits further aided in the discovery of SNP variants/Indels, QTLs and genes that can ascertain the genetic determinants of fruit characters. Through multi-omics approaches such as genomics, transcriptomics, metabolomics and proteomics, the identification and quantification of transcripts, including non-coding RNAs, involved in sugar metabolism, fruit development and ripening, shelf life, and the biotic and abiotic stress that impacts fruit quality were made possible. Utilizing genomic assisted breeding methods such as genome wide association (GWAS), genomic selection (GS) and genetic modifications using CRISPR/Cas9 and transgenics has paved the way to studying gene function and developing cultivars with desirable fruit traits by overcoming long breeding cycles. Such comprehensive multi-omics approaches related to fruit characters in tropical fruits and their applications in breeding strategies and crop improvement are reviewed, discussed and presented here.

Transcriptome analysis of low-temperature-affected ripening revealed MYB transcription factors-mediated regulatory network in banana fruit

Low temperature leads to abnormal ripening and poor quality of the harvested banana fruit, which is an urgent problem limiting the development of industry in China. To comprehensively understand the mechanism underlying low-temperature-affected ripening, we performed comparative RNA-Seq analysis of ethylene-induced ripening of banana fruit after 3 days of pre-storage at 7 °C and 22 °C. A total of 986 differentially expressed genes (DEGs) were identified in both RT-0 d versus RT-3 d, LT-0 d versus LT-3 d, RT-0 d versus LT-0 d and RT-3 d versus LT-3 d, and the RNA-Seq outputs of 15 randomly selective DEGs were verified using qRT-PCR. Among the 986 DEGs obtained in the four groups, 9 MYB genes (MaMYB75/281/219/4/151/156/3/37 and MaMYB3R1) and 32 genes related to carotenoid biosynthesis (MaPSY1/2a), flavor formation (MaLOX6, MaADH7, MaAAT1), sucrose transport (MaSUS2/4), ethylene production (MaSAM1, MaACO9/10/12, MaACS1/12), starch degradation (MaAMY1A/1B, MaPHS1/2, MaMEX2, MapGlcT1) and cell wall degradation (MaPG3/X1, MaPME25/41, MaXTH5/7/22/23/25, MaEXP2/20/A1/A15) were characterized. Combining transcription factor binding site (TFBS) analysis as well as cis-acting element analysis, the regulatory network of low-temperature-affected ripening mediated by MYBs were constructed. The data generated in this study may unravel the transcriptional regulatory network of MYBs associated with low-temperature-affected ripening and provide a solid foundation for future studies.

Functional Characterisation of Banana (Musa spp.) 2-Oxoglutarate-Dependent Dioxygenases Involved in Flavonoid Biosynthesis

Bananas (Musa) are non-grass, monocotyledonous, perennial plants that are well known for their edible fruits. Their cultivation provides food security and employment opportunities in many countries. Banana fruits contain high levels of minerals and phytochemicals, including flavonoids, which are beneficial for human nutrition. To broaden the knowledge on flavonoid biosynthesis in this major crop plant, we aimed to identify and functionally characterise selected structural genes encoding 2-oxoglutarate-dependent dioxygenases, involved in the formation of the flavonoid aglycon. Musa candidates genes predicted to encode flavanone 3-hydroxylase (F3H), flavonol synthase (FLS) and anthocyanidin synthase (ANS) were assayed. Enzymatic functionalities of the recombinant proteins were confirmed in vivo using bioconversion assays. Moreover, transgenic analyses in corresponding Arabidopsis thaliana mutants showed that MusaF3H, MusaFLS and MusaANS were able to complement the respective loss-of-function phenotypes, thus verifying functionality of the enzymes in planta. Knowledge gained from this work provides a new aspect for further research towards genetic engineering of flavonoid biosynthesis in banana fruits to increase their antioxidant activity and nutritional value.

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

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

The R2R3-MYB transcription factor MtMYB134 orchestrates flavonol biosynthesis in Medicago truncatula

Our results provide insights into the flavonol biosynthesis regulation of M. truncatula. The R2R3-MYB transcription factor MtMYB134 emerged as tool to improve the flavonol biosynthesis. Flavonols are plant specialized metabolites with vital roles in plant development and defense and are known as diet compound beneficial to human health. In leguminous plants, the regulatory proteins involved in flavonol biosynthesis are not well characterized. Using a homology-based approach, three R2R3-MYB transcription factor encoding genes have been identified in the Medicago truncatula reference genome sequence. The gene encoding a protein with highest similarity to known flavonol regulators, MtMYB134, was chosen for further experiments and was characterized as a functional flavonol regulator from M. truncatula. MtMYB134 expression levels are correlated with the expression of MtFLS2, encoding a key enzyme of flavonol biosynthesis, and with flavonol metabolite content. MtMYB134 was shown to activate the promoters of the A. thaliana flavonol biosynthesis genes AtCHS and AtFLS1 in Arabidopsis protoplasts in a transactivation assay and to interact with the Medicago promoters of MtCHS2 and MtFLS2 in yeast 1-hybrid assays. To ascertain the functional aspect of the identified transcription factor, we developed a sextuple mutant, which is defective in anthocyanin and flavonol biosynthesis. Ectopic expression of MtMYB134 in a multiple myb A. thaliana mutant restored flavonol biosynthesis. Furthermore, overexpression of MtMYB134 in hairy roots of M. truncatula enhanced the biosynthesis of various flavonol derivatives. Taken together, our results provide insight into the understanding of flavonol biosynthesis regulation in M. truncatula and provides MtMYB134 as tool for genetic manipulation to improve flavonol synthesis.

Transcriptome Sequence Reveals Candidate Genes Involving in the Post-Harvest Hardening of Trifoliate Yam Dioscorea dumetorum

Storage ability of trifoliate yam (Dioscorea dumetorum) is restricted by a severe post-harvest hardening (PHH) phenomenon, which starts within the first 24 h after harvest and renders tubers inedible. Previous work has only focused on the biochemical changes affecting PHH in D. dumetorum. To the best of our knowledge, the candidate genes responsible for the hardening of D. dumetorum have not been identified. Here, transcriptome analyses of D. dumetorum tubers were performed in yam tubers of four developmental stages: 4 months after emergence (4MAE), immediately after harvest (AH), 3 days after harvest (3DAH) and 14 days after harvest (14DAH) of four accessions (Bangou 1, Bayangam 2, Fonkouankem 1, and Ibo sweet 3) using RNA-Seq. In total, between AH and 3DAH, 165, 199, 128 and 61 differentially expressed genes (DEGs) were detected in Bayangam 2, Fonkouankem 1, Bangou 1 and Ibo sweet 3, respectively. Functional analysis of DEGs revealed that genes encoding for CELLULOSE SYNTHASE A (CESA), XYLAN O-ACETYLTRANSFERASE (XOAT), CHLOROPHYLL A/B BINDING PROTEIN1, 2, 3, 4 (LHCB1, LHCB2, LHCB3, and LCH4) and an MYB transcription factor were predominantly and significantly up-regulated 3DAH, implying that these genes were potentially involved in the PHH as confirmed by qRT-PCR. A hypothetical mechanism of this phenomenon and its regulation has been proposed. These findings provide the first comprehensive insights into gene expression in yam tubers after harvest and valuable information for molecular breeding against the PHH.

Characterization of the Liriodendron Chinense MYB Gene Family and Its Role in Abiotic Stress Response

Liriodendron chinense (Lchi) is a Magnoliaceae plant, which is a basic angiosperm left behind by the Pleistocene and mainly distributed in the south of the Yangtze River. Liriodendron hybrids has good wood properties and is widely used in furniture and in other fields. It is not clear if they can adapt to different environmental conditions, such as drought and high and low temperatures, and the molecular mechanisms for this adaptation are unknown. Among plant transcription factors (TFs), the MYB gene family is one of the largest and is often involved in stress or adversity response signaling, growth, and development. Therefore, studying the role of MYBTFs in regulating abiotic stress signaling, growth, and development in Lchi is helpful to promote afforestation in different environments. In our research, a genome-wide analysis of the LchiMYB gene family was performed, including the phylogenetic relationship tree, gene exon-intron structure, collinearity, and chromosomal position. According to the evolutionary tree, 190 LchiMYBs were divided into three main branches. LchiMYBs were evenly distributed across 19 chromosomes, with their collinearity, suggesting that segment duplication events may have contributed to LchiMYB gene expansion. Transcriptomes from eight tissues, 11 stages of somatic embryogenesis, and leaves after cold, heat, and drought stress were used to analyze the function of the MYB gene family. The results of tissue expression analysis showed that most LchiMYB genes regulated bark, leaf, bud, sepal, stigma, and stamen development, as well as the four important stages (ES3, ES4, ES9, and PL) of somatic embryogenesis. More than 60 LchiMYBs responded to heat, cold, and drought stress; some of which underwent gene duplication during evolution. LchiMYB3 was highly expressed under all three forms of stress, while LchiMYB121 was strongly induced by both cold and heat stress. Eight genes with different expression patterns were selected and verified by quantitative real-time PCR (qRT-PCR) experiments. The results suggested that these LchiMYBs may regulate Lchi growth development and resistance to abiotic stress. This study shows the cross-regulatory function of LchiMYBs in the growth and development, asexual reproduction, and abiotic resistance of Lchi. This information will prove pivotal to directing further studies on the biological function of Lchi MYBTFs in genetic improvement and abiotic stress response.

Automatic Identification of Players in the Flavonoid Biosynthesis with Application on the Biomedicinal Plant Croton tiglium

The flavonoid biosynthesis is a well-characterised model system for specialised metabolism and transcriptional regulation in plants. Flavonoids have numerous biological functions such as UV protection and pollinator attraction, but also biotechnological potential. Here, we present Knowledge-based Identification of Pathway Enzymes (KIPEs) as an automatic approach for the identification of players in the flavonoid biosynthesis. KIPEs combines comprehensive sequence similarity analyses with the inspection of functionally relevant amino acid residues and domains in subjected peptide sequences. Comprehensive sequence sets of flavonoid biosynthesis enzymes and knowledge about functionally relevant amino acids were collected. As a proof of concept, KIPEs was applied to investigate the flavonoid biosynthesis of the medicinal plant Croton tiglium on the basis of a transcriptome assembly. Enzyme candidates for all steps in the biosynthesis network were identified and matched to previous reports of corresponding metabolites in Croton species.