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

Genome-Wide Analysis of Transcription Factor R2R3-MYB Gene Family and Gene Expression Profiles during Anthocyanin Synthesis in Common Walnut (Juglans regia L.)

The R2R3-MYB gene family, encoding plant transcriptional regulators, participates in many metabolic pathways of plant physiology and development, including flavonoid metabolism and anthocyanin synthesis. This study proceeded as follows: the JrR2R3-MYB gene family was analyzed genome-wide, and the family members were identified and characterized using the high-quality walnut reference genome "Chandler 2.0". All 204 JrR2R3-MYBs were established and categorized into 30 subgroups via phylogenetic analysis. JrR2R3-MYBs were unevenly distributed over 16 chromosomes. Most JrR2R3-MYBs had similar structures and conservative motifs. The cis-acting elements exhibit multiple functions of JrR2R3-MYBs such as light response, metabolite response, and stress response. We found that the expansion of JrR2R3-MYBs was mainly caused by WGD or segmental duplication events. Ka/Ks analysis indicated that these genes were in a state of negative purifying selection. Transcriptome results suggested that JrR2R3-MYBs were widely entangled in the process of walnut organ development and differentially expressed in different colored varieties of walnuts. Subsequently, we identified 17 differentially expressed JrR2R3-MYBs, 9 of which may regulate anthocyanin biosynthesis based on the results of a phylogenetic analysis. These genes were present in greater expression levels in 'Zijing' leaves than in 'Lvling' leaves, as revealed by the results of qRT-PCR experiments. These results contributed to the elucidation of the functions of JrR2R3-MYBs in walnut coloration. Collectively, this work provides a foundation for exploring the functional characteristics of the JrR2R3-MYBs in walnuts and improving the nutritional value and appearance quality of walnuts.

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.

Genome-wide identification of PEBP gene family in pineapple reveal its potential functions in flowering

Phosphatidylethanolamine binding protein (PEBP) plays an important role in regulating flowering time and morphogenesis of plants. However, the identification and functional analysis of PEBP gene in pineapple (AcPEBP) have not been systematically studied. The pineapple genome contained 11 PEBP family members, which were subsequently classified into three subfamilies (FT-like, TFL-like and MFT-like) based on phylogenetic relationships. The arrangement of these 11 shows an unequal pattern across the six chromosomes of pineapple the pineapple genome. The anticipated outcomes of the promoter cis-acting elements indicate that the PEBP gene is subject to regulation by diverse light signals and endogenous hormones such as ethylene. The findings from transcriptome examination and quantitative real-time polymerase chain reaction (qRT-PCR) indicate that FT-like members AcFT3 and AcFT4 display a heightened expression level, specifically within the floral structures. The expression of AcFT3 and AcFT4 increases sharply and remains at a high level after 4 days of ethylene induction, while the expression of AcFT7 and AcMFT1 decreases gradually during the flowering process. Additionally, AcFT3, AcFT4 and AcFT7 show specific expression in different floral organs of pineapple. These outcomes imply that members belonging to the FT-like subfamily may have a significant impact on the process of bud differentiation and flower development. Through transcriptional activation analysis, it was determined that AcFT4 possesses transcriptional activation capability and is situated in the nucleus and peripheral cytoplasm. Overexpression of AcFT4 in Arabidopsis resulted in the promotion of early flowering by 6-7 days. The protein interaction prediction network identified potential flower regulators, including CO, AP1, LFY and SOC1, that may interact with PEBP proteins. This study explores flower development in pineapple, thereby serving as a valuable reference for future research endeavors in this domain.

A New Function of MbIAA19 Identified to Modulate Malus Plants Dwarfing Growth

The primary determinants of apple (Malus) tree architecture include plant height and internode length, which are the significant criteria for evaluating apple dwarf rootstocks. Plant height and internode length are predominantly governed by phytohormones. In this study, we aimed to assess the mechanisms underlying dwarfism in a mutant of Malus baccata. M. baccata dwarf mutant (Dwf) was previously obtained through natural mutation. It has considerably reduced plant height and internode length. A comparative transcriptome analysis of wild-type (WT) and Dwf mutant was performed to identify and annotate the differentially expressed genes responsible for the Dwf phenotype using RNA-seq and GO and KEGG pathway enrichment analyses. Multiple DEGs involved in hormone signaling pathways, particularly auxin signaling pathways, were identified. Moreover, the levels of endogenous indole-3-acetic acid (IAA) were lower in Dwf mutant than in WT. The Aux/IAA transcription factor gene MbIAA19 was downregulated in Dwf mutant due to a single nucleotide sequence change in its promoter. Genetic transformation assay demonstrated strong association between MbIAA19 and the dwarf phenotype. RNAi-IAA19 lines clearly exhibited reduced plant height, internode length, and endogenous IAA levels. Our study revealed that MbIAA19 plays a role in the regulation of dwarfism and endogenous IAA levels in M. baccata.

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.

Transcriptomic analysis reveals the mechanism underlying the anthocyanin changes in Fragaria nilgerrensis Schlecht. and its interspecific hybrids

BACKGROUND

Fragaria nilgerrensis (FN) provides a rich source of genetic variations for strawberry germplasm innovation. The color of strawberry fruits is a key factor affecting consumer preferences. However, the genetic basis of the fruit color formation in F. nilgerrensis and its interspecific hybrids has rarely been researched.

RESULTS

In this study, the fruit transcriptomes and flavonoid contents of FN (white skin; control) and its interspecific hybrids BF1 and BF2 (pale red skin) were compared. A total of 31 flavonoids were identified. Notably, two pelargonidin derivatives (pelargonidin-3-O-glucoside and pelargonidin-3-O-rutinoside) were revealed as potential key pigments for the coloration of BF1 and BF2 fruits. Additionally, dihydroflavonol 4-reductase (DFR) (LOC101293459 and LOC101293749) and anthocyanidin 3-O-glucosyltransferase (BZ1) (LOC101300000), which are crucial structural genes in the anthocyanidin biosynthetic pathway, had significantly up-regulated expression levels in the two FN interspecific hybrids. Moreover, most of the genes encoding transcription factors (e.g., MYB, WRKY, TCP, bHLH, AP2, and WD40) related to anthocyanin accumulation were differentially expressed. We also identified two DFR genes (LOC101293749 and LOC101293459) that were significantly correlated with members in bHLH, MYB, WD40, AP2, and bZIP families. Two chalcone synthase (CHS) (LOC101298162 and LOC101298456) and a BZ1 gene (LOC101300000) were highly correlated with members in bHLH, WD40 and AP2 families.

CONCLUSIONS

Pelargonidin-3-O-glucoside and pelargonidin-3-O-rutinoside may be the key pigments contributing to the formation of pale red fruit skin. DFR and BZ1 structural genes and some bHLH, MYB, WD40, AP2, and bZIP TF family members enhance the accumulation of two pelargonidin derivatives. This study provides important insights into the regulation of anthocyanidin biosynthesis in FN and its interspecific hybrids. The presented data may be relevant for improving strawberry fruit coloration via genetic engineering.

Identification and characterization of the gene BraANS.A03 associated with purple leaf color in pak choi (Brassica rapa L. ssp. chinensis)

MAIN CONCLUSION

BraANS.A3 was the key gene controlling purple leaf color in pak choi, and two short fragments of promoter region in green pak choi might be interfering its normal expression. Pak choi (B. rapa L. ssp. chinensis) is an influential and important vegetable with green, yellow, or purple leaves that is cultivated worldwide. The purple leaves are rich in anthocyanins, but the underlying genetics and evolution have yet to be extensively studied. Free-hand sections of the purple leaves indicated that anthocyanins mainly accumulate throughout the adaxial and abaxial epidermal leaf cells. Segregation analyses of an F2 population of a B. rapa ssp. chinensis L. purple leaf mutant ZBC indicated that the purple trait is controlled by an incompletely dominant nuclear gene. Bulked segregant analysis (BSA) showed that the key genes controlling the trait were between 24.25 and 38.10 Mb on chromosome A03 of B. rapa. From the annotated genes, only BraA03g050560.3C, homologous to Arabidopsis AtANS, was related to the anthocyanin synthesis pathway. Genome annotation results and transcriptional sequencing analyses revealed that the BraANS.A3 gene was involved in the purple leaf trait. qRT-PCR analyses showed that BraANS.A3 was highly upregulated in ZBC but hardly expressed in the leaves of an inbred homozygous line of B. campestris ssp. chinensis L. green leaf mutant WTC, indicating that BraANS.A3 played a key role catalyzing anthocyanin synthesis in ZBC. Full-length sequence alignment of BraANS.A3 in WTC and ZBC showed that it was highly conserved in the gene region, with significant variation in the promoter region. In particular, the insertion of two short fragments of the promoter region in WTC may interfere with its normal expression. The promoter regions of ANS in six Brassica species all had multiple cis-elements involved in responses to abscisic acid, light, and stress, suggesting that ANS may be involved in multiple metabolic pathways or biological processes. Protein-protein interactions predicted that BraANS.A3 interacts with virtually all catalytic proteins in the anthocyanin synthesis pathway and has a strong relationship with Transparent Testa 8 (TT8). These results suggest that BraANS.A3 promotes anthocyanin accumulation in purple pak choi and provide new insights into the functional analysis of anthocyanin-related genes in Chinese cabbage and transcriptional regulatory networks.

Transcriptome-wide expression analysis of MYB gene family leads to functional characterization of flavonoid biosynthesis in fruit coloration of Ziziphus Mill

The Ziziphus mauritiana Lam. and Z. jujuba Mill. are the two most economically important members of the genus Ziziphus. The fruit color of Z. mauritiana remains green throughout fruit development in the majority of commercial cultivars, whereas its close relative, Z. jujuba Mill. turns from green to red in all cultivars. However, the lack of transcriptomic and genomic information confines our understanding of the molecular mechanisms underlying fruit coloration in Z. mauritiana (Ber). In the present study, we performed the transcriptome-wide analysis of MYB transcription factors (TFs) genes in Z. mauritiana and Z. jujuba, and identified 56 ZmMYB and 60 ZjMYB TFs in Z. mauritiana and Z. jujuba, respectively. Through transcriptomic expression analysis, four similar MYB genes (ZmMYB/ZjMYB13, ZmMYB/ZjMYB44, ZmMYB/ZjMYB50, and ZmMYB/ZjMYB56) from Z. mauritiana and Z. jujuba were selected as candidate key genes regulating flavonoid biosynthesis. Among these genes, the ZjMYB44 gene was transiently highly expressed in fruit, and flavonoid content accumulation also increased, indicating that this gene can influence flavonoid content during the period of fruit coloration in Z. jujuba. The current study adds to our understanding of the classification of genes, motif structure, and predicted functions of the MYB TFs, as well as identifying MYBs that regulate flavonoid biosynthesis in Ziziphus (Z. mauritiana and Z. jujuba). Based on this information, we concluded that MYB44 is involved in the flavonoids biosynthesis pathway during the fruit coloring of Ziziphus. Our research results provide an important understanding of the molecular mechanism of flavonoid biosynthesis resulting in fruit coloration and laying a foundation for further genetic improvement of fruit color in Ziziphus.

Genome-wide characterization and expression analysis of MYB transcription factors in Chrysanthemum nankingense

BACKGROUND

Chrysanthemum is a popular ornamental plant worldwide. MYB (v-myb avian myeloblastosis viral oncogene homolog) transcription factors play an important role in everything from stress resistance to plant growth and development. However, the MYB family of chrysanthemums has not been the subject of a detailed bioinformatics and expression investigation.

RESULTS

In this study, we examined 324 CnMYB transcription factors from Chrysanthemum nankingense genome data, which contained 122 Cn1R-MYB, 183 CnR2R3-MYB, 12 Cn3R-MYB, 2 Cn4R-MYB, and 5 atypical CnMYB. The protein motifs and classification of CnMYB transcription factors were analyzed. Among them, motifs 1, 2, 3, and 4 were found to encode the MYB DNA-binding domain in R2R3-MYB proteins, while in other-MYB proteins, the motifs 1, 2, 3, 4, 5, 6, 7, and 8 encode the MYB DNA-binding domain. Among all CnMYBs, 44 genes were selected due to the presence of CpG islands, while methylation is detected in three genes, including CnMYB9, CnMYB152, and CnMYB219. We analyzed the expression levels of each CnMYB gene in ray floret, disc floret, flower bud, leaf, stem, and root tissues. Based on phylogenetic analysis and gene expression analysis, three genes appeared likely to control cellulose and lignin synthesis in stem tissue, and 16 genes appeared likely to regulate flowering time, anther, pollen development, and flower color. Fifty-one candidate genes that may be involved in stress response were identified through phylogenetic, stress-responseve motif of promoter, and qRT-PCR analyses. According to genes expression levels under stress conditions, six CnMYB genes (CnMYB9, CnMYB172, CnMYB186, CnMYB199, CnMYB219, and CnMYB152) were identified as key stress-responsive genes.

CONCLUSIONS

This research provides useful information for further functional analysis of the CnMYB gene family in chrysanthemums, as well as offers candidate genes for further study of cellulose and lignin synthesis, flowering traits, salt and drought stress mechanism.

Genome-wide identification and expression analysis of the anthocyanin-related genes during seed coat development in six Brassica species

Yellow seed is one favorite trait for the breeding of Brassica oilseed crops, but the performance of seed coat color is very complicated due to the involvement of various pigments. The change of seed coat color of Brassica crops is related to the specific synthesis and accumulation of anthocyanin, and the expression level of structural genes in anthocyanin synthesis pathway is specifically regulated by transcription factors. Despite some previous reports on the regulations of seed coat color from linkage marker development, gene fine-mapping and multi-omics association analysis, the trait of Brassica crops is affected by the evolutionary events such as genome triploidization, the regulatory mechanism is still largely unknown. In this study, we identified genes related to anthocyanin synthesis in six Brassica crops in U-triangle at the genome-wide level and performed collinearity analysis. A total of 1119 anthocyanin-related genes were identified, the collinear relationship of anthocyanin-related genes on subgenomic chromosomes was the best in B. napus (AACC) and the worst in B. carinata (BBCC). The comparisons of gene expressions for anthocyanin metabolic pathways in seed coats during seed development revealed differences in its metabolism among these species. Interestingly, the R2R3-MYB transcription factors MYB5 and TT2 were differentially expressed at all eight stages of seed coat development, indicating that they might be the key genes that caused the variation of the seed coat color. The expression curve and trend analyses of the seed coat development period showed that the main reason for the unexpressed copies of MYB5 and TT2 was likely gene silencing caused by gene structural variation. These results were valuable for the genetic improvement of Brassica seed coat color, and also provided new insights into gene multicopy evolution in Brassica polyploids.

Characterization of the REVEILLE family in Rosaceae and role of PbLHY in flowering time regulation

BACKGROUND

The circadian clock integrates endogenous and exogenous signals and regulates various physiological processes in plants. REVEILLE (RVE) proteins play critical roles in circadian clock system, especially CCA1 (CIRCADIAN CLOCK ASSOCIATED 1) and LHY (LATE ELONGATED HYPOCOTYL), which also participate in flowering regulation. However, little is known about the evolution and function of the RVE family in Rosaceae species, especially in Pyrus bretschneideri.

RESULTS

In this study, we performed a genome-wide analysis and identified 51 RVE genes in seven Rosaceae species. The RVE family members were classified into two groups based on phylogenetic analysis. Dispersed duplication events and purifying selection were the main drivers of evolution in the RVE family. Moreover, the expression patterns of ten PbRVE genes were diverse in P. bretschneideri tissues. All PbRVE genes showed diurnal rhythms under light/dark cycles in P. bretschneideri leaves. Four PbRVE genes also displayed robust rhythms under constant light conditions. PbLHY, the gene with the highest homology to AtCCA1 and AtLHY in P. bretschneideri, is localized in the nucleus. Ectopic overexpression of PbLHY in Arabidopsis delayed flowering time and repressed the expression of flowering time-related genes.

CONCLUSION

These results contribute to improving the understanding and functional research of RVE genes in P. bretschneideri.

Identification of MYB Transcription Factors Involving in Fruit Quality Regulation of Fragaria × ananassa Duch

The cultivated strawberry (Fragaria × ananassa Duch.) is an important horticultural crop. The economic values of strawberry cultivars are decided by their fruit qualities including taste, color and aroma. The important role of MYB transcription factors in fruit quality regulation is recognized increasingly with the identification of MYB genes involved in metabolism. A total of 407 MYB genes of F. × ananassa (FaMYBs) were identified in the genome-wide scale and named according to subgenome locations. The 407 FaMYBs were clustered into 36 groups based on phylogenetic analysis. According to synteny analysis, whole genome duplication and segmental duplication contributed over 90% of the expansion of the FaMYBs family. A total of 101 FaMYB loci with 1-6 alleles were identified by the homologous gene groups on homologous chromosomes. The differentially expressed FaMYB profiles of three cultivars with different fruit quality and fruit ripe processes provided the 8 candidate loci involved in fruit quality regulation. In this experiment, 7, 5, and 4 FaMYBs were screeded as candidate genes involved in the regulation of metabolism/transportation of anthocyanins, sugars or organic acids and 4-hydroxy-2, 5-dimethyl-3(2H)-furanone, respectively. These results pointed out the key FaMYBs for further functional analysis of gene regulation of strawberry fruit quality and would be helpful in the clarification on ofe roles of MYBs in the metabolism of fruit crops.

Insights into transcription factors controlling strawberry fruit development and ripening

Fruit ripening is a highly regulated and complex process involving a series of physiological and biochemical changes aiming to maximize fruit organoleptic traits to attract herbivores, maximizing therefore seed dispersal. Furthermore, this process is of key importance for fruit quality and therefore consumer acceptance. In fleshy fruits, ripening involves an alteration in color, in the content of sugars, organic acids and secondary metabolites, such as volatile compounds, which influence flavor and aroma, and the remodeling of cell walls, resulting in the softening of the fruit. The mechanisms underlying these processes rely on the action of phytohormones, transcription factors and epigenetic modifications. Strawberry fruit is considered a model of non-climacteric species, as its ripening is mainly controlled by abscisic acid. Besides the role of phytohormones in the regulation of strawberry fruit ripening, a number of transcription factors have been identified as important regulators of these processes to date. In this review, we present a comprehensive overview of the current knowledge on the role of transcription factors in the regulation of strawberry fruit ripening, as well as in compiling candidate regulators that might play an important role but that have not been functionally studied to date.

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

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

Genome-Wide Analysis of Ascorbic Acid Metabolism Related Genes in Fragaria × ananassa and Its Expression Pattern Analysis in Strawberry Fruits

Ascorbic acid (AsA) is an important antioxidant for scavenging reactive oxygen species and it is essential for human health. Strawberry (Fragaria × ananassa) fruits are rich in AsA. In recent years, strawberry has been regarded as a model for non-climacteric fruit ripening. However, in contrast to climacteric fruits, such as tomato, the regulatory mechanism of AsA accumulation in strawberry fruits remains largely unknown. In this study, we first identified 125 AsA metabolism-related genes from the cultivated strawberry "Camarosa" genome. The expression pattern analysis using an available RNA-seq data showed that the AsA biosynthetic-related genes in the D-mannose/L-galactose pathway were downregulated remarkably during fruit ripening which was opposite to the increasing AsA content in fruits. The D-galacturonate reductase gene (GalUR) in the D-Galacturonic acid pathway was extremely upregulated in strawberry receptacles during fruit ripening. The FaGalUR gene above belongs to the aldo-keto reductases (AKR) superfamily and has been proposed to participate in AsA biosynthesis in strawberry fruits. To explore whether there are other genes in the AKR superfamily involved in regulating AsA accumulation during strawberry fruit ripening, we further implemented a genome-wide analysis of the AKR superfamily using the octoploid strawberry genome. A total of 80 FaAKR genes were identified from the genome and divided into 20 subgroups based on phylogenetic analysis. These FaAKR genes were unevenly distributed on 23 chromosomes. Among them, nine genes showed increased expression in receptacles as the fruit ripened, and notably, FaAKR23 was the most dramatically upregulated FaAKR gene in receptacles. Compared with fruits at green stage, its expression level increased by 142-fold at red stage. The qRT-PCR results supported that the expression of FaAKR23 was increased significantly during fruit ripening. In particular, the FaAKR23 was the only FaAKR gene that was significantly upregulated by abscisic acid (ABA) and suppressed by nordihydroguaiaretic acid (NDGA, an ABA biosynthesis blocker), indicating FaAKR23 might play important roles in ABA-mediated strawberry fruit ripening. In a word, our study provides useful information on the AsA metabolism during strawberry fruit ripening and will help understand the mechanism of AsA accumulation in strawberry fruits.

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

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