MYB transcription factor family in sweet cherry (Prunus avium L.): genome-wide investigation, evolution, structure, characterization and expression patterns

Identifying potential anthocyanin biosynthesis regulator in Chinese cherry by comprehensive genome-wide characterization of the R2R3-MYB transcription factor gene family

BACKGROUND

Chinese cherry [Cerasus pseudocerasus (Lindl.) G.Don] (syn. Prunus pseudocerasus Lindl.) is an economically important fruiting cherry species with a diverse range of attractive colors, spanning from the lightest yellow to the darkest black purple. However, the MYB transcription factors involved in anthocyanin biosynthesis underlying fruit color variation in Chinese cherry remain unknown.

RESULTS

In this study, we characterized the R2R3-MYB gene family of Chinese cherry by genome-wide identification and compared it with those of 10 Rosaceae relatives and Arabidopsis thaliana. A total of 1490 R2R3-MYBs were classified into 43 subfamilies, which included 29 subfamilies containing both Rosaceae MYBs and AtMYBs. One subfamily (S45) contained only Rosaceae MYBs, while three subfamilies (S12, S75, and S77) contained only AtMYBs. The variation in gene numbers within identical subfamilies among different species and the absence of certain subfamilies in some species indicated the species-specific expansion within MYB gene family in Chinese cherry and its relatives. Segmental and tandem duplication events primarily contributed to the expansion of Chinese cherry R2R3-CpMYBs. The duplicated gene pairs underwent purifying selection during evolution after duplication events. Phylogenetic relationships and transcript profiling revealed that CpMYB10 and CpMYB4 are involved in the regulation of anthocyanin biosynthesis in Chinese cherry fruits. Expression patterns, transient overexpression and VIGS results confirmed that CpMYB10 promotes anthocyanin accumulation in the fruit skin, while CpMYB4 acts as a repressor, inhibiting anthocyanin biosynthesis of Chinese cherry.

CONCLUSIONS

This study provides a comprehensive and systematic analysis of R2R3-MYB gene family in Chinese cherry and Rosaceae relatives, and identifies two regulators, CpMYB10 and CpMYB4, involved in anthocyanin biosynthesis in Chinese cherry. These results help to develop and utilize the potential functions of anthocyanins in Chinese cherry.

Transcriptomic Analysis Reveals the Flavonoid Biosynthesis Pathway Involved in Rhizome Development in Polygonatum cyrtonema Hua

Polygonatum cyrtonema Hua (P. cyrtonema) rhizomes are rich in flavonoids and other secondary metabolites, exhibiting remarkable antioxidant, anti-tumor, and immunomodulatory effects. Polygonatum flavonoid-biosynthesis-related genes have been characterized already. However, a comprehensive overview of Polygonatum flavonoid biosynthesis pathways is still absent. To articulate the accumulation of the flavonoid biosynthesis pathways, we examined transcriptome changes using Illumina HiSeq from five different tissues and the RNA-seq of 15 samples had over 105 Gb of a clean base, generating a total of 277,955 unigenes. The cDNA libraries of the fruits (F), leaves (L), roots (R), stems (S), and rhizomes (T) of three-year-old P. cyrtonema plants generated 57,591, 53,578, 60,321, 51,530, and 54,935 unigenes. Comparative transcriptome analysis revealed that 379 differentially expressed genes (DEGs) were in the group of F _vs_ T, L _vs_ T, R _vs_ T, and S _vs_ T, and the transcripts of flavonoid-biosynthesis-related DEGs were principally enriched in rhizomes. In addition, combined with WGCNA and the FPKM of five tissues' transcription, nine differentially expressed transcription factor families (MYB, WRKY, AP2/ERF, etc.) were characterized in the red module, the red module positively correlated with rhizome flavonoid accumulation. Quantitative real-time PCR (qRT-PCR) further indicated that BZIP1, C3H31, ERF114, and DREB21 are differentially expressed in rhizomes, accompanied in rhizome development in P. cyrtonema. Therefore, this study provides a foundation for further research into uncovering the accumulation of flavonoid biosynthesis in the rhizomes of P. cyrtonema.

Abscisic acid-induced transcription factor PsMYB306 negatively regulates tree peony bud dormancy release

Bud dormancy is a crucial strategy for perennial plants to withstand adverse winter conditions. However, the regulatory mechanism of bud dormancy in tree peony (Paeonia suffruticosa) remains largely unknown. Here, we observed dramatically reduced and increased accumulation of abscisic acid (ABA) and bioactive gibberellins (GAs) GA1 and GA3, respectively, during bud endodormancy release of tree peony under prolonged chilling treatment. An Illumina RNA sequencing study was performed to identify potential genes involved in the bud endodormancy regulation in tree peony. Correlation matrix, principal component, and interaction network analyses identified a downregulated MYB transcription factor gene, PsMYB306, the expression of which positively correlated with 9-CIS-EPOXYCAROTENOID DIOXYGENASE 3 (PsNCED3) expression. Protein modeling analysis revealed 4 residues within the R2R3 domain of PsMYB306 to possess DNA binding capability. Transcription of PsMYB306 was increased by ABA treatment. Overexpression of PsMYB306 in petunia (Petunia hybrida) inhibited seed germination and plant growth, concomitant with elevated ABA and decreased GA contents. Silencing of PsMYB306 accelerated cold-triggered tree peony bud burst and influenced the production of ABA and GAs and the expression of their biosynthetic genes. ABA application reduced bud dormancy release and transcription of ENT-KAURENOIC ACID OXIDASE 1 (PsKAO1), GA20-OXIDASE 1 (PsGA20ox1), and GA3-OXIDASE 1 (PsGA3ox1) associated with GA biosynthesis in PsMYB306-silenced buds. In vivo and in vitro binding assays confirmed that PsMYB306 specifically transactivated the promoter of PsNCED3. Silencing of PsNCED3 also promoted bud break and growth. Altogether, our findings suggest that PsMYB306 negatively modulates cold-induced bud endodormancy release by regulating ABA production.

Comparative genomics of N-acetyl-5-methoxytryptamine members in four Prunus species with insights into bud dormancy and abiotic stress responses in Prunus avium

KEY MESSAGE

This study provides novel insights into the evolution, diversification, and functions of melatonin biosynthesis genes in Prunus species, highlighting their potential role in regulating bud dormancy and abiotic stresses. The biosynthesis of melatonin (MEL) in plants is primarily governed by enzymatic reactions involving key enzymes such as serotonin N-acetyltransferase (SNAT), tryptamine 5-hydroxylase (T5H), N-acetylserotonin methyltransferase (ASMT) and tryptophan decarboxylase (TDC). In this study, we analyzed Melatonin genes in four Prunus species such as Prunus avium (Pavi), Prunus pusilliflora (Ppus), Prunus serulata (Pser), and Prunus persica (Pper) based on comparative genomics approach. Among the four Prunus species, a total of 29 TDCs, 998 T5Hs, 16 SNATs, and 115 ASMTs within the genome of four Prunus genomes. A thorough investigation of melatonin-related genes was carried out using systematic biological methods and comparative genomics. Through phylogenetic analysis, orthologous clusters, Go enrichment, syntenic relationship, and gene duplication analysis, we discovered both similarities and variations in Melatonin genes among these Prunus species. Additionally, our study revealed the existence of unique subgroup members in the Melatonin genes of these species, which were distinct from those found in Arabidopsis genes. Furthermore, the transcriptomic expression analysis revealed the potential significance of melatonin genes in bud dormancy regulation and abiotic stresses. Our extensive results offer valuable perspectives on the evolutionary patterns, intricate expansion, and functions of PavMEL genes. Given their promising attributes, PavTDCs, PavT5H, PavNAT, and three PavASMT genes warrant in-depth exploration as prime candidates for manipulating dormancy in sweet cherry. This was done to lay the foundation for future explorations into the structural and functional aspects of these factors in Prunus species. This study offers significant insights into the functions of ASMT, SNAT, T5H, and TDC genes and sheds light on their roles in Prunus avium. Moreover, it established a robust foundation for further exploration functional characterization of melatonin genes in fruit species.

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.

Genome-wide identification, evolution, and role of SPL gene family in beet (Beta vulgaris L.) under cold stress

BACKGROUND

SPL transcription factors play vital roles in regulating plant growth, development, and abiotic stress responses. Sugar beet (Beta vulgaris L.), one of the world's main sugar-producing crops, is a major source of edible and industrial sugars for humans. Although the SPL gene family has been extensively identified in other species, no reports on the SPL gene family in sugar beet are available.

RESULTS

Eight BvSPL genes were identified at the whole-genome level and were renamed based on their positions on the chromosome. The gene structure, SBP domain sequences, and phylogenetic relationship with Arabidopsis were analyzed for the sugar beet SPL gene family. The eight BvSPL genes were divided into six groups (II, IV, V, VI, VII, and VIII). Of the BvSPL genes, no tandem duplication events were found, but one pair of segmental duplications was present. Multiple cis-regulatory elements related to growth and development were identified in the 2000-bp region upstream of the BvSPL gene start codon (ATG). Using quantitative real-time polymerase chain reaction (qRT-PCR), the expression profiles of the eight BvSPL genes were examined under eight types of abiotic stress and during the maturation stage. BvSPL transcription factors played a vital role in abiotic stress, with BvSPL3 and BvSPL6 being particularly noteworthy.

CONCLUSION

Eight sugar beet SPL genes were identified at the whole-genome level. Phylogenetic trees, gene structures, gene duplication events, and expression profiles were investigated. The qRT-PCR analysis indicated that BvSPLs play a substantial role in the growth and development of sugar beet, potentially participating in the regulation of root expansion and sugar accumulation.

Genome-Wide Analysis of Cotton MYB Transcription Factors and the Functional Validation of GhMYB in Response to Drought Stress

MYB transcription factors play important roles during abiotic stress responses in plants. However, little is known about the accurate systematic analysis of MYB genes in the four cotton species, Gossypium hirsutum, G. barbadense, G. arboreum and G. raimondii. Herein, we performed phylogenetic analysis and showed that cotton MYBs and Arabidopsis MYBs were clustered in the same subfamilies for each species. The identified cotton MYBs were distributed unevenly on chromosomes in various densities for each species, wherein genome-wide tandem and segment duplications were the main driving force of MYB family expansion. Synteny analysis suggested that the abundant collinearity pairs of MYBs were identified between G. hirsutum and the other three species, and that they might have undergone strong purification selection. Characteristics of conserved motifs, along with their consensus sequence, promoter cis elements and gene structure, revealed that MYB proteins might be highly conserved in the same subgroups for each species. Subsequent analysis of differentially expressed genes and expression patterns indicated that most GhMYBs might be involved in response to drought (especially) and salt stress, which was supported by the expression levels of nine GhMYBs using real-time quantitative PCR. Finally, we performed a workflow that combined virus-induced gene silencing and the heterologous transformation of Arabidopsis, which confirmed the positive roles of GhMYBs under drought conditions, as validated by determining the drought-tolerant phenotypes, damage index and/or water loss rate. Collectively, our findings not only expand our understanding of the relationships between evolution and function of MYB genes, but they also provide candidate genes for cotton breeding.

Genome-wide identification, classification and expression analysis of MYB gene family in coconut (Cocos nucifera L.)

MYB transcription factors regulate the growth, development, and secondary metabolism of plant species. To investigate the origin of color variations in coconut pericarp, we identified and analyzed the MYB gene family present in coconut. According to the sequence of MYB genes in Arabidopsis thaliana, homologous MYB gene sequences were found in the whole genome database of coconut, the conserved sequence motifs within MYB proteins were analyzed by Motif Elicitation (MEME) tool, and the sequences without conservative structure were eliminated. Additionally, we employed RNA-seq technology to generate gene expression signatures of the R2R3-MYB genes across distinctive coconut parts exhibiting diverse colors. To validate these profiles, we conducted quantitative PCR (qPCR). Through comprehensive genome-wide screening, we successfully identified a collection of 179 MYB genes in coconut. Subsequent phylogenetic analysis categorized these 179 coconut MYB genes into 4-subfamilies: 124 R2R3-MYB, 4 3R-MYB types, 4 4R-MYB type, and 47 unknown types. Furthermore, these genes were further divided into 34 subgroups, with 28 of these subgroups successfully classified into known subfamilies found in Arabidopsis thaliana. By mapping the CnMYB genes onto the 16 chromosomes of the coconut genome, we unveiled a collinearity association between them. Moreover, a preservation of gene structure and motif distribution was observed across the CnMYB genes. Our research encompassed a thorough investigation of the R2R3-MYB genes present in the coconut genome, including the chromosomal localization, gene assembly, conserved regions, phylogenetic associations, and promoter cis-acting elements of the studied genes. Our findings revealed a collection of 12 R2R3-MYB candidate genes, namely CnMYB8, CnMYB15, CnMYB27, CnMYB28, CnMYB61, CnMYB63, CnMYB68, CnMYB94, CnMYB101, CnMYB150, CnMYB153, and CnMYB164. These genes showed differential expressions in diverse tissues and developmental stages of four coconut species, such as CnMYB68, CnMYB101, and CnMYB28 exhibited high expression in majority of tissues and coconut species, while CnMYB94 and CnMYB164 showed lower expression. These findings shed light on the crucial functional divergence of CnMYB genes across various coconut tissues, suggesting these genes as promising candidate genes for facilitating color development in this important crop.

Sweet cherry TCP gene family analysis reveals potential functions of PavTCP1, PavTCP2 and PavTCP3 in fruit light responses

BACKGROUND

TCP proteins are plant specific transcription factors that play important roles in plant growth and development. Despite the known significance of these transcription factors in general plant development, their specific role in fruit growth remains largely uncharted. Therefore, this study explores the potential role of TCP transcription factors in the growth and development of sweet cherry fruits.

RESULTS

Thirteen members of the PavTCP family were identified within the sweet cherry plant, with two, PavTCP1 and PavTCP4, found to contain potential target sites for Pav-miR159, Pav-miR139a, and Pav-miR139b-3p. Analyses of cis-acting elements and Arabidopsis homology prediction analyses that the PavTCP family comprises many light-responsive elements. Homologs of PavTCP1 and PavTCP3 in Arabidopsis TCP proteins were found to be crucial to light responses. Shading experiments showed distinct correlation patterns between PavTCP1, 2, and 3 and total anthocyanins, soluble sugars, and soluble solids in sweet cherry fruits. These observations suggest that these genes may contribute significantly to sweet cherry light responses. In particular, PavTCP1 could play a key role, potentially mediated through Pav-miR159, Pav-miR139a, and Pav-miR139b-3p.

CONCLUSION

This study is the first to unveil the potential function of TCP transcription factors in the light responses of sweet cherry fruits, paving the way for future investigations into the role of this transcription factor family in plant fruit development.

Unveiling the effect of gibberellin-induced iron oxide nanoparticles on bud dormancy release in sweet cherry (Prunus avium L.)

Hydrogen cyanide has been extensively used worldwide for bud dormancy break in fruit trees, consequently enhancing fruit production via expedited cultivation, especially in areas with controlled environments or warmer regions. A novel and safety nanotechnology was developed since the hazard of hydrogen cyanide for the operators and environments, there is an urgent need for the development of novel and safety approaches to replace it to break bud dormancy for fruit trees. In current study, we have systematically explored the potential of iron oxide nanoparticles, specifically α-Fe2O3, to modulate bud dormancy in sweet cherry (Prunus avium). The synthesized iron oxide nanoparticles underwent meticulous characterization and assessment using various techniques, including Fourier-transform infrared spectroscopy (FTIR), X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), transmission electron microscopy (TEM), and ultraviolet-visible infrared (UV-Vis) spectroscopy. Remarkably, when applied at a concentration of 10 mg L-1 alongside gibberellin (GA4+7), these iron oxide nanoparticles exhibited a substantial 57% enhancement in bud dormancy release compared to control groups, all achieved within a remarkably short time span of 4 days. Our RNA-seq analyses further unveiled that 2757 genes within the sweet cherry buds were significantly up-regulated when treated with 10 mg L-1 α-Fe2O3 nanoparticles in combination with GA, while 4748 genes related to dormancy regulation were downregulated in comparison to the control. Moreover, we discovered an array of 58 transcription factor families among the crucial differentially expressed genes (DEGs). Through hormonal quantification, we established that the increased bud burst was accompanied by a reduced concentration of abscisic acid (ABA) at 761.3 ng/g fresh weight in the iron oxide treatment group, coupled with higher levels of gibberellins (GAs) in comparison to the control. Comprehensive transcriptomic and metabolomic analyses unveiled significant alterations in hormone contents and gene expression during the bud dormancy-breaking process when α-Fe2O3 nanoparticles were combined with GA. In conclusion, our findings provide valuable insights into the intricate molecular mechanisms underlying the impact of iron oxide nanoparticles on achieving uniform bud dormancy break in sweet cherry trees.

Genome-wide identification, expression analysis, and potential roles under low-temperature stress of bHLH gene family in Prunus sibirica

The basic helix-loop-helix (bHLH) family is one of the most well-known transcription factor families in plants, and it regulates growth, development, and abiotic stress responses. However, systematic analyses of the bHLH gene family in Prunus sibirica have not been reported to date. In this study, 104 PsbHLHs were identified and classified into 23 subfamilies that were unevenly distributed on eight chromosomes. Nineteen pairs of segmental replication genes and ten pairs of tandem replication genes were identified, and all duplicated gene pairs were under purifying selection. PsbHLHs of the same subfamily usually share similar motif compositions and exon-intron structures. PsbHLHs contain multiple stress-responsive elements. PsbHLHs exhibit functional diversity by interacting and coordinating with other members. Twenty PsbHLHs showed varying degrees of expression. Eleven genes up-regulated and nine genes down-regulated in -4°C. The majority of PsbHLHs were highly expressed in the roots and pistils. Transient transfection experiments demonstrated that transgenic plants with overexpressed PsbHLH42 have better cold tolerance. In conclusion, the results of this study have significant implications for future research on the involvement of bHLH genes in the development and stress responses of Prunus sibirica.

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.

Genome-wide identification, bioinformatics and expression analysis of HD-Zip gene family in peach

BACKGROUND

HD-Zips (Homeodomain-Leucine Zippers) are a class of plant-specific transcription factors that play multiple roles in plant growth and development. Although some functions of HD-Zip transcription factor have been reported in several plants, it has not been comprehensively studied in peach, especially during adventitious root formation of peach cuttings.

RESULTS

In this study, 23 HD-Zip genes distributed on 6 chromosomes were identified from the peach (Prunus persica) genome, and named PpHDZ01-23 according to their positions on the chromosomes. These 23 PpHDZ transcription factors all contained a homeomorphism box domain and a leucine zipper domain, were divided into 4 subfamilies(I-IV) according to the evolutionary analysis, and their promoters contained many different cis-acting elements. Spatio-temporal expression pattern showed that these genes were expressed in many tissues with different levels, and they had distinct expression pattern during adventitious root formation and development.

CONCLUSION

Our results showed the roles of PpHDZs on root formation, which is helpful to better understand the classification and function of peach HD-Zip genes.

Genome-Wide Analysis of MYB Transcription Factors and Screening of MYBs Involved in the Red Color Formation in Rhododendron delavayi

Flower color is one of the crucial traits of ornamental plants. Rhododendron delavayi Franch. is a famous ornamental plant species distributed in the mountain areas of Southwest China. This plant has red inflorescence and young branchlets. However, the molecular basis of the color formation of R. delavayi is unclear. In this study, 184 MYB genes were identified based on the released genome of R. delavayi. These genes included 78 1R-MYB, 101 R2R3-MYB, 4 3R-MYB, and 1 4R-MYB. The MYBs were divided into 35 subgroups using phylogenetic analysis of the MYBs of Arabidopsis thaliana. The members of the same subgroup in R. delavayi had similar conserved domains and motifs, gene structures, and promoter cis-acting elements, which indicate their relatively conserved function. In addition, transcriptome based on unique molecular identifier strategy and color difference of the spotted petals, unspotted petals, spotted throat, unspotted throat, and branchlet cortex were detected. Results showed significant differences in the expression levels of R2R3-MYB genes. Weighted co-expression network analysis between transcriptome and chromatic aberration values of five types of red samples showed that the MYBs were the most important TFs involved in the color formation, of which seven were R2R3-MYB, and three were 1R-MYB. Two R2R3-MYB (DUH019226.1 and DUH019400.1) had the highest connectivity in the whole regulation network, and they were identified as hub genes for red color formation. These two MYB hub genes provide references for the study of transcriptional regulation of the red color formation of R. delavayi.

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

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

Integrated Transcriptome and Metabolome Analyses Provide Insights into the Coloring Mechanism of Dark-red and Yellow Fruits in Chinese Cherry [Cerasus pseudocerasus (Lindl.) G. Don]

Chinese cherry [Cerasus pseudocerasus (Lindl.) G. Don] is an important fruit tree from China that has excellent ornamental, economic, and nutritional values with various colors. The dark-red or red coloration of fruit, an attractive trait for consumers, is determined by anthocyanin pigmentation. In this study, the coloring patterns during fruit development in dark-red and yellow Chinese cherry fruits were firstly illustrated by integrated transcriptome and widely-targeted metabolome analyses. Anthocyanin accumulation in dark-red fruits was significantly higher compared with yellow fruits from the color conversion period, being positively correlated to the color ratio. Based on transcriptome analysis, eight structural genes (CpCHS, CpCHI, CpF3H, CpF3'H, CpDFR, CpANS, CpUFGT, and CpGST) were significantly upregulated in dark-red fruits from the color conversion period, especially CpANS, CpUFGT, and CpGST. On contrary, the expression level of CpLAR were considerably higher in yellow fruits than in dark-red fruits, especially at the early stage. Eight regulatory genes (CpMYB4, CpMYB10, CpMYB20, CpMYB306, bHLH1, CpNAC10, CpERF106, and CpbZIP4) were also identified as determinants of fruit color in Chinese cherry. Liquid chromatography-tandem mass spectrometry identified 33 and 3 differential expressed metabolites related to anthocyanins and procyanidins between mature dark-red and yellow fruits. Cyanidin-3-O-rutinoside was the predominant anthocyanin compound in both fruits, while it was 6.23-fold higher in dark-red than in yellow fruits. More accumulated flavanol and procyanidin contents resulted in less anthocyanin content in flavonoid pathway in yellow fruits due to the higher expression level of CpLAR. These findings can help understand the coloring mechanism of dark-red and yellow fruits in Chinese cherry, and provide genetic basis for breeding new cultivars.

Investigation and Expression Analysis of R2R3-MYBs and Anthocyanin Biosynthesis-Related Genes during Seed Color Development of Common Bean (Phaseolus vulgaris)

Anthocyanins are responsible for the coloration of common bean seeds, and their accumulation is positively correlated with the expression level of anthocyanin biosynthetic genes. The MBW (MYB-bHLH-WD40) complex is thought to regulate the expression of these genes, and MYB proteins, which are a key factor in activating anthocyanin pathway genes, have been identified in several plants. This study demonstrated gene structures, chromosomal placements, gene duplications of R2R3-MYBs, miRNAs associated with R2R3-MYBs, and the interaction of these genes with other flavonoid regulatory genes. qRT-PCR was used to investigate the role of specific R2R3-MYBs and flavonoid genes in common bean seed color development. As a result of a comprehensive analysis with the help of in silico tools, we identified 160 R2R3-MYB genes in the common bean genome. We divided these genes into 16 classes on the basis of their intron-exon and motif structures. Except for three, the rest of the common bean R2R3-MYB members were distributed to all chromosomes with different densities, primarily located on chromosomes 3 and 8. We identified a total of 44 duplicated gene pairs dispersed across 11 chromosomes and evolved under purifying selection (Ka/Ks  <  1), 19 of which were derived from a whole-genome duplication. Our research uncovered 25 putative repressor PvMYB proteins that contain the EAR motif. Additionally, fifty different cis-regulatory elements regulated by light, stress, and hormone were identified. Within the genome of the common bean, we discovered a total of 36 microRNAs that target a total of 72 R2R3-MYB transcripts. The effect of 16 R2R3-MYB genes and 16 phenylpropanoid pathway genes, selected on the basis of their interaction in the protein-protein interaction map, playing role in the regulation of seed coat color development was evaluated using qRT-PCR in 5 different tissues at different developmental stages. The results revealed that these specific genes have different expression levels during different developmental periods, with higher levels in the pod filling and early pod stages than in the rest of the developmental periods. Furthermore, it was shown that PvTT8 (bHLH), PvTT2 (PvMYB42), PvMYB113, PvTTG1, and PvWD68 genes have effects on the regulation of seed coat color. The findings of this study, which is the first to use whole-genome analysis to identify and characterize the R2R3-MYB genes in common bean, may serve as a reference for future functional research in the legume.

High-resolution genome-wide association study of a large Czech collection of sweet cherry (Prunus avium L.) on fruit maturity and quality traits

In sweet cherry (Prunus avium L.), quantitative trait loci have been identified for fruit maturity, colour, firmness, and size to develop markers for marker-assisted selection. However, resolution is usually too low in those analyses to directly target candidate genes, and some associations are missed. In contrast, genome-wide association studies are performed on broad collections of accessions, and assemblies of reference sequences from Tieton and Satonishiki cultivars enable identification of single nucleotide polymorphisms after whole-genome sequencing, providing high marker density. Two hundred and thirty-five sweet cherry accessions were sequenced and phenotyped for harvest time and fruit colour, firmness, and size. Genome-wide association studies were used to identify single nucleotide polymorphisms associated with each trait, which were verified in breeding material consisting of 64 additional accessions. A total of 1 767 106 single nucleotide polymorphisms were identified. At that density, significant single nucleotide polymorphisms could be linked to co-inherited haplotype blocks (median size ~10 kb). Thus, markers were tightly associated with respective phenotypes, and individual allelic combinations of particular single nucleotide polymorphisms provided links to distinct phenotypes. In addition, yellow-fruit accessions were sequenced, and a ~ 90-kb-deletion on chromosome 3 that included five MYB10 transcription factors was associated with the phenotype. Overall, the study confirmed numerous quantitative trait loci from bi-parental populations using high-diversity accession populations, identified novel associations, and genome-wide association studies reduced the size of trait-associated loci from megabases to kilobases and to a few candidate genes per locus. Thus, a framework is provided to develop molecular markers and evaluate and characterize genes underlying important agronomic traits.

Transcriptome-wide identification of walnut PP2C family genes in response to external stimulus

Walnut is an important economic tree species while confronting with global environmental stress, resulting in decline in quality and yield. Therefore, it is urgent to elucidate the molecular mechanism for the regulation of walnut response to adversity. The protein phosphatase 2C (PP2C) gene family participates in cellular processes in eukaryotes through reversible phosphorylation of proteins and signal transduction regulation. However, the stress response function of PP2C genes was far to be clarified. Therefore, to understand the stress response mechanism of walnut tree, in this study, a total of 41 PP2C genes with complete ORFs were identified from Juglans regia, whose basic bio-information and expression patterns in response to multiple stresses and ABA were confirmed. The results showed that the ORFs of JrPP2Cs were 495 ~ 3231 bp in length, the predicted JrPP2C proteins contained 164 to 1076 amino acids and the molecular weights were 18,581.96 ~ 118,853.34 Da, the pI was 4.55 ~ 9.58. These JrPP2C genes were unevenly distributed on 14 chromosomes, among which Chr11 and Chr13 contained the most genes. Phylogenetic analysis found that these JrPP2C proteins were classed into 9 subfamilies, among which group F covered most JrPP2Cs. The JrPP2Cs in the same subfamily exhibited similarities in the composition of conserved domains, amino acid sequences of motifs and exon/intron organization in DNA sequences. Each JrPP2C includes 4 ~ 10 motifs and each motif contained 15 ~ 37 amino acids. Among the motifs, motif1, motif2, motif3 and motif8 were most abundant. Most of the JrPP2C genes diversely response to osmotic, cadmium, and Colletotrichum gloeosporioide stress as well as ABA treatments, among which JrPP2C28, JrPP2C17, JrPP2C09, JrPP2C36 were more obvious and deserves further attention. All these results indicated that JrPP2C genes play potential vital roles in plant response to multiple stimulus, and are possibly involved in ABA-dependent signaling pathway.

Evolutionary and Integrative Analysis of Gibberellin-Dioxygenase Gene Family and Their Expression Profile in Three Rosaceae Genomes (F. vesca, P. mume, and P. avium) Under Phytohormone Stress

The gibberellin-dioxygenase (GAox) gene family plays a crucial role in regulating plant growth and development. GAoxs, which are encoded by many gene subfamilies, are extremely critical in regulating bioactive GA levels by catalyzing the subsequent stages in the biosynthesis process. Moreover, GAoxs are important enzymes in the GA synthesis pathway, and the GAox gene family has not yet been identified in Rosaceae species (Prunus avium L., F. vesca, and P. mume), especially in response to gibberellin and PCa (prohexadione calcium; reduce biologically active GAs). In the current investigation, 399 GAox members were identified in sweet cherry, Japanese apricot, and strawberry. Moreover, they were further classified into six (A-F) subgroups based on phylogeny. According to motif analysis and gene structure, the majority of the PavGAox genes have a remarkably well-maintained exon-intron and motif arrangement within the same subgroup, which may lead to functional divergence. In the systematic investigation, PavGAox genes have several duplication events, but segmental duplication occurs frequently. A calculative analysis of orthologous gene pairs in Prunus avium L., F. vesca, and P. mume revealed that GAox genes are subjected to purifying selection during the evolutionary process, resulting in functional divergence. The analysis of cis-regulatory elements in the upstream region of the 140 PavGAox members suggests a possible relationship between genes and specific functions of hormone response-related elements. Moreover, the PavGAox genes display a variety of tissue expression patterns in diverse tissues, with most of the PavGAox genes displaying tissue-specific expression patterns. Furthermore, most of the PavGAox genes express significant expression in buds under phytohormonal stresses. Phytohormones stress analysis demonstrated that some of PavGAox genes are responsible for maintaining the GA level in plant-like Pav co4017001.1 g010.1.br, Pav sc0000024.1 g340.1.br, and Pav sc0000024.1 g270.1.mk. The subcellular localization of PavGAox protein utilizing a tobacco transient transformation system into the tobacco epidermal cells predicted that GFP signals were mostly found in the cytoplasm. These findings will contribute to a better understanding of the GAox gene family's interaction with prohexadione calcium and GA, as well as provide a strong framework for future functional characterization of GAox genes in sweet cherry.

Metabolites in Cherry Buds to Detect Winter Dormancy

Winter dormancy is still a “black box” in phenological models, because it evades simple observation. This study presents the first step in the identification of suitable metabolites which could indicate the timing and length of dormancy phases for the sweet cherry cultivar ‘Summit’. Global metabolite profiling detected 445 named metabolites in flower buds, which can be assigned to different substance groups such as amino acids, carbohydrates, phytohormones, lipids, nucleotides, peptides and some secondary metabolites. During the phases of endo- and ecodormancy, the energy metabolism in the form of glycolysis and the tricarboxylic acid (TCA) cycle was shut down to a minimum. However, the beginning of ontogenetic development was closely related to the up-regulation of the carbohydrate metabolism and thus to the generation of energy for the growth and development of the sweet cherry buds. From the 445 metabolites found in cherry buds, seven were selected which could be suitable markers for the ecodormancy phase, whose duration is limited by the date of endodormancy release (t₁) and the beginning of ontogenetic development (t₁*). With the exception of abscisic acid (ABA), which has been proven to control bud dormancy, all of these metabolites show nearly constant intensity during this phase.

Identification and Comprehensive Genome-Wide Analysis of Glutathione S-Transferase Gene Family in Sweet Cherry (Prunus avium) and Their Expression Profiling Reveals a Likely Role in Anthocyanin Accumulation

Glutathione S-transferases (GSTs) in plants are multipurpose enzymes that are involved in growth and development and anthocyanins transportation. However, members of the GST gene family were not identified in sweet cherry (Prunus avium). To identify the GST genes in sweet cherry, a genome-wide analysis was conducted. In this study, we identified 67 GST genes in P. avium genome and nomenclature according to chromosomal distribution. Phylogenetic tree analysis revealed that PavGST genes were classified into seven chief subfamily: TCHQD, Theta, Phi, Zeta, Lambda, DHAR, and Tau. The majority of the PavGST genes had a relatively well-maintained exon-intron and motif arrangement within the same group, according to gene structure and motif analyses. Gene structure (introns-exons) and conserved motif analysis revealed that the majority of the PavGST genes showed a relatively well-maintained motif and exons-introns configuration within the same group. The chromosomal localization, GO enrichment annotation, subcellular localization, syntenic relationship, Ka/Ks analysis, and molecular characteristics were accomplished using various bioinformatics tools. Mode of gene duplication showed that dispersed duplication might play a key role in the expansion of PavGST gene family. Promoter regions of PavGST genes contain numerous cis-regulatory components, which are involved in multiple stress responses, such as abiotic stress and phytohormones responsive factors. Furthermore, the expression profile of sweet cherry PavGSTs showed significant results under LED treatment. Our findings provide the groundwork for future research into induced LED anthocyanin and antioxidants deposition in sweet cherries.

Genome-wide identification of myeloblastosis gene family and its response to cadmium stress in Ipomoea aquatica

The myeloblastosis (MYB) proteins perform key functions in mediating cadmium (Cd) tolerance of plants. Ipomoea aquatica has strong adaptability to Cd Stress, while the roles of the I. aquatica MYB gene family with respect to Cd stress are still unclear. Here, we identified a total of 183 MYB genes in the I. aquatica genome (laMYB), which were classified into 66 1R-type IaMYB, 112 2R-type IaMYB, four 3R-type IaMYB, and one 4R-type IaMYB based on the number of the MYB repeat in each gene. The analysis of phylogenetic tree indicated that most of IaMYB genes are associated with the diverse biological processes including defense, development and metabolism. Analysis of sequence features showed that the IaMYB genes within identical subfamily have the similar patterns of the motif distributions and gene structures. Analysis of gene duplication events revealed that the dispersed duplication (DSD) and whole-genome duplication (WGD) modes play vital roles in the expansion of the IaMYB gene family. Expression profiling manifests that approximately 20% of IaMYB genes had significant role in the roots of I. aquatica under Cd stress. Promoter profiling implied that the differentially expressed genes might be induced by environmental factors or inherent hormones and thereby execute their function in Cd response. Remarkably, the 2R-type IaMYB157 with abundant light-responsive element G-box and ABA-responsive element ABRE in its promoter region exhibited very strong response to Cd stress. Taken together, our findings provide an important candidate IaMYB gene for further deciphering the molecular regulatory mechanism in plant with respect to Cd stress.