Molecular evolution and expression assessment of DFRs in apple

Anthocyanins are the secondary metabolites of flavonoids in plants. As a key enzyme in the biosynthetic pathway of anthocyanin, dihydroflavonol 4-reductase (DFR) act as an important regulatory point, but DFR family genes has not been systematically characterized in apple (Malus domestica Borkh.). The members of DFR genes in apple were identified and their gene structure, chromosome distribution, evolutionary relationships, collinearity, cis-component and protein interaction relationships were predicted using bioinformatics methods. The expression patterns of MdDFRs in various organs, such as leaves, fruit flushes, fruits, ripe fruit peels, flowers and stems were analyzed using GeneChip expression array analysis. qRT-PCR was employed to analyze the expression levels of MdDFRs in different apple varieties with varying levels of fruit skin at maturity. The apple database revealed 96 DFR genes, which are distributed on 17 chromosomes and can be divided into 3 subfamilies. These 96 DFR genes were mostly composed of α-helix and random coil according to secondary structure prediction, and were mainly expressed in chloroplasts and cytoplasm. MYB binding site involved in flavonoid biosynthetic genes regulation element (MBSI) was identified in the promoter of MdDFR15/76/81/89/90/91/93/94. Lignin/flavonoid synthesis-related elements of MYB recognition site and MYB-binding site were identified in the promoters of MdDFR05/09/13/19/22/24/26/30/31/33/34/46/50/52/54/64/65/69/75/76/79/86. The internal collinearity analysis of the apple MdDFR genome revealed a total of 34 pairs of duplicated gene pairs. Interspecific collinearity analysis showed that there were 66 and 57 homologous gene pairs in apple/tomato and apple/grape, respectively. GeneChip expression array analysis showed that MdDFR72 and MdDFR96 were higher expressed in ripe fruit fleshes and peel, MdDFR01/06/67/49/54/91 were higher expressed in flowers, MdDFR64 was higher expressed in ripe fruit peels and flowers than those of other tissues. Besides, 75 MdDFR proteins interacted directly or indirectly with anthocyanidin synthesis related proteins MdANS, MdF3H, MdMYB1, MdMYBPA1 to form a protein interaction network. Interestingly, MdDFR69 and MdDFR87 had direct interactions with these four proteins, MdDFR64 had direct interactions with MdANS and MdF3H. qRT-PCR analysis showed that the expression levels of MdDFR01/05/31/53/64/69/73/84/87/94/96 were up-regulated with the accumulation of anthocyanins. This study lays a foundation for further research on the function of DFR genes in apple.


Background
Apple (Malus domestica Borkh.) is a perennial woody plant of the Rosaceae family and is one of the most popular fruits in China, with a wide cultivation area and high yield throughout the country.Apple fruit is rich in minerals and vitamins, and its nutrient components are soluble, easy to be absorbed by the human body, with the role of calming the nerves, helping sleep, and reducing cholesterol [1].Apple fruits also contain a certain amount of anthocyanins, which has the function of anti-aging and preventing cardiovascular and cerebrovascular diseases [2].Anthocyanidin and proanthocyanidin are two end products of flavonoid biosynthesis pathway, which are believed to be synthesized and modificated in cytoplasm and endoplasmic reticulum before being sequestered into vacuoles [3].Anthocyanins are synthesized by phenylalanine via phenylalanine lyase (PAL), cinnamate-4-hydroxylase (C4H), 4-coumaryl CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3′-hydroxylase (F3′H), and flavonoid 3′5′-hydroxylase (F3′5′H), dihydroflavonol 4-reductase (DFR), and anthocyanin synthase (ANS).DFR catalyzed stereospecific reduction reaction of dihydroflavonol at the C4 position, which enables dihydroflavonol to generate colorless proanthocyanidins, and is one of the most critical enzymes in anthocyanin synthesis [4].Notably, DFR enzymes exhibit striking substrate specificities, resulting in accumulation of distinct patterns of anthocyanins, catechins and proanthocyanidins [5].
F3H, as the upstream of DFR, can catalyze the synthesis of dihydroflavonols from flavonoids [6].DFR further catalyzes catalyze three different colorless dihydroflavonols, including dihydro kaempferol (DHK), dihydro quercetin (DHQ) and dihydromyricetin (DHM), to generate different anthocyanidin precursors, and then synthesize various anthocyanin under the action of ANS [7].It can be seen that DFR is the first enzyme to convert dihydroflavonol into anthocyanin.DFR can use DHK, DHQ, and DHM as substrates, respectively.The functional expression of Saussurea medusa gene SmDFR in yeast (Saccharomyces cerevisiae) confirms that SmDFR has a reducing effect on DHQ and DHK, but not on DHM [8].The heterologous expression of grape hyacinth (Muscari spp.) gene MaDFR in tobacco (Nicotiana tabacum) resulted in increasing anthocyanin accumulation, leading to a darker flower color [9].However, due to the different substrate selectivity of different species of DFR, different anthocyanins are synthesized, presenting different plant colors.Structural genes encoding enzymes of the anthocyanin pathways are regulated at the transcriptional level by a MYB-bHLH-WD40 (MBW) ternary complex of R2R3-MYB, basic helix-loop-helix (bHLH), and tryptophan aspartic (WD40)-repeat proteins [10,11].Among these

Graphical Abstract
regulatory factors, the R2R3-MYB transcription factors, which play a crucial role in determining the activation of specific downstream genes to control the spatiotemporal accumulation of anthocyanins, have been widely reported in apple, including MdMYBPA1 and MdMYB1 [12,13].The regulatory interaction between transcription factors related to anthocyanin synthesis and structural genes has been a focus of research [14][15][16].
Our precious study reported that the expression of DFR was upregulated with the increase of anthocyanin accumulation revealed by comprehensive transcriptome analyses [17].Due to the substrate specificity of DFR, different members of its gene family have different catalytic efficiency for each substrate and different expression levels in different tissues, which causes different tissues or organs of plants show various colors [18].Given the importance of DFRs in the regulation of flavonoid synthesis and the production and accumulation of anthocyanidin, the expression and functional characterization of DFR family members in anthocyanidin biosynthesis of apple have not been systematically investigated, especially in different apple organs and varieties, as well as the relationships with other identified anthocyanidin synthesis related genes or transcription factors.Consequently, in light of the discovery of the apple DFR gene members, this experiment performed a bioinformatics analysis of the MdDFR gene members and predicted the interaction network between the MdDFRs and the anthocyanidin synthesis related proteins MdANS, MdF3H, MdMYB1, MdMYBPA1.The role of DFR genes in apple anthocyanin synthesis was examined through the analysis of GeneChip expression array of different apple tissues and expression levels of apples with different coloring degrees, which provided theoretical reference for the study of the regulatory network of apple anthocyanin synthesis.

Plant materials
Apple varieties 'Golden Delicious' , 'Ruixue' , 'Chengji No.1' and ' Astar' with different coloring degrees at maturity were selected as research materials, and named G1, G2, G3 and G4, respectively.Fifteen fruits were collected from each variety, with every five fruits being a duplicate.Fruit skin samples were accurately weighed and quickly frozen in liquid nitrogen for 2 h, and then stored in a refrigerator at − 80 °C for RNA extraction.

Identification of DFR gene members in apple
The FASTA sequence of apple DFR gene CDS was obtained by searching gene names in NCBI network database (https:// www.ncbi.nlm.nih.gov/).InterPro number (IPR001509) of this gene was obtained from the Pfam website (https:// pfam.xfam.org/) using the obtained FASTA sequence.Collect all the genes containing Inter-Pro number in Phytozome v13 (https:// phyto zome-next.jgi.doe.gov/).The mRNA and CDS sequences of genes were downloaded from the apple reference genome website (https:// iris.angers.inra.fr/ gddh13/) using the accession number, and the corresponding protein sequences were obtained from the translation website (https:// web.expasy.org/ trans late/).Chromosome localization of DFR gene members of apple was performed using TBtools (Version 1.108) software.

Bioinformatics analysis of DFR gene members in apple
ProtParam tool in Expasy database (https:// web.expasy.org/ protp aram/) was used to analyze the basic physicochemical properties of protein, such as the number, relative molecular weight and instability coefficient.ClustalX, MEGA7.0 and EvolView (https:// www.evolg enius.info/) were used to construct a phylogenetic tree using the Neighbor-Joining method (NJ) with the bootstrap value of 1000.The gene structure was analyzed by GSDS 2.0 (http:// gsds.cbi.pku.edu.cn/).The conserved motifs of MdDFRs were determined using the MEME online tool (http:// meme-suite.org/ tools/ meme) the number of motifs output was 10.The motif analysis and gene structure were visualized using the TBtools software.The subcellular localization prediction of DFR gene members was performed online using WoLF PSORT (http:// www.gensc ript.com/ wolf-psort.html).The secondary structure prediction of DFR proteins was carried out using PredictProtein (https:// predi ctpro tein.org/).The chromosomal location information of MdDFRs was obtained from the apple genome website (https:// iris.angers.inra.fr/ gddh13/) and was drawn a picture by MG2C (http:// mg2c.iask.in/ mg2c_ v2.0/).

Intraspecific and interspecific collinearity analysis of MdDFR genes
Download the genomes and annotated files of tomato and grape using Phytozome v13.The intraspecific collinearity of MdDFR gene pairs as well as the interspecific collinearity of DFR members from Malus domestica/Solanum lycopersicum, Malus domestica/Vitis vinifera L. was calculated and visualized by MCScanX method at TBtools (Version 1.108).

Expression pattern and protein interaction analysis of DFR genes in apple
The FPKM values of 96 MdDFRs in different organs, including leaves, young leaves, fruit flushes, fruits, ripe fruit peels, flowers and stems, were obtained through the apple MDO multidimensional omics database (http:// bioin forma tics.cau.edu.cn/ Apple MDO/ gene_ detail.php) and standardized with log 10 .The expression heat maps were drawn using TBtools.The protein interaction network was predicted by STRING Version 11 (https:// string-db.org/).

Quantitative RT-PCR (qRT-PCR) analysis
Total RNA was extracted from 100 mg of frozen tissue using an RNA plant (Real Times, Beijing, China) basing on the manufacturer's instructions.The first strand cDNA fragment was synthesized from 1.0 μg of DNase created RNA in a reaction solution (total volume of 20 μL) by using the Evo M-MLV RT kit (Accurate Biology, Changsha, China).Gene specific primers (Additional file 1: Table S1) were designed and synthesized by Sangon Biotech Co., Ltd.(China).Gene expression was detected using qRT-PCR using SYBR ® Green Premix Pro Taq HS qPCR Kit (Accurate Biology).Using apple GADPH (Gen-Bank access NO.CB973647) gene as internal reference.Light Cycler ® 96 RealTime PCR System (Roche, Basel, Switzerland) was used to determine the expression levels of MdDFRs in the peel of four apple cultivars.The reaction system is 20 μL, including 2 mL cDNA, 1 μL upstream and downstream primers, 10 μL SYBR, and 20 μL ddH 2 O complement.The PCR reaction was as follows: 95 ℃ for 30 s, 40 cycles at 95 ℃ for 5 s and 60 ℃ for 1 min, 95 ℃ for 5 s, 55 ℃ for 30 s and 72 ℃ for 1 min.The relative expression levels of genes were calculated using the 2 − △△CT method [19].

Physicochemical properties and subcellular localization analysis of DFR genes in apple
A total of 96 DFR genes, named MdDFR01-MdDFR96, were obtained through the homologous search on the apple reference genome website, and distributed on 17 apple chromosomes (Fig. 1, Table 1).Among them, chromosome 2 (chr2) and chr11 were the most distributed with 12 genes, followed by chr7 with 9 genes, chr6 and chr8 with 8 genes, chr0 and chr14 with 1 gene each.The number of amino acids encoded by this gene members was ranged from 75 aa (MdDFR10) to 1204 aa (MdDFR70).The relative molecular weight ranged from 8344.64 D (MdDFR10) to 136,293.11D (MdDFR70).The isoelectric point (PI) ranged from 5.02 (MdDFR78) to 9.96 (MdDFR47), including 28 basic proteins (PI > 7.5), 44 acidic proteins (PI < 6.5), and 24 neutral proteins (PI 6.5-7.5).The instability index ranged from 18.26 (MdDFR29) to 49.93 (MdDFR05).There were 23 hydrophilic proteins and 73 unstable hydrophobic proteins.The full length of the mRNA ranged from 353 bp (MdDFR40) to 10,524 bp (MdDFR01).The full length of the cDNA ranged from 228 bp (MdDFR10) to 3615 bp (MdDFR70).Among the 96 apple DFR genes, 47 were located in the chloroplast, 31 were located in the cytoplasm, and there were 4 genes each in cytoskeleton, extracellar matrix, and mitochondria.There were 3 genes each in nucleus and endoplasmic reticulum, and only 1 member in peroxisome (Table 1).

Secondary structure prediction analysis of DFR proteins in apple
The secondary structure of the protein encoded by the apple DFR genes was mainly α-helix and random coil, followed by extended strand, and β-corner structure was the least (Table 2).The proportion of α-helix structure was between 29.84% (MdDFR23) and 56.06% (MdDFR70).The proportion of irregular curly structures ranged from 26.67% (MdDFR10) to 48.17% (MdDFR01).Extended strand accounted for 10.71% (MdDFR40) to 23.00% (MdDFR77).β-corner accounted for the smallest proportion of 4.13% (MdDFR18) to 12.00% (MdDFR10).According to the above analysis, it can be concluded that α-helix and random coil were the main factors in the formation of spiral structure of apple DFRs.

Phylogenetic analysis of DFR genes in apple
To further understand the evolution and functional characteristics of apple DFR genes, the phylogenetic tree of apple DFR members was constructed and divided into three subfamilies (Clade I-III, Fig. 2).Among them, the Clade I subfamily includes 7 genes, including MdDFR19, MdDFR22, MdDFR54, MdDFR62, MdDFR79, MdDFR89, and MdDFR95.The Clade II subfamily includes 40 genes, and the Clade III subfamily includes 49 apple DFR genes (Fig. 2).In each branch, apple DFR genes have high homology and close evolutionary relationship, but there are some differences, which preliminarily indicates that their evolution is conservative, and it is speculated that the functions of the lineal homologous genes are similar.

Gene structure and motif sequence analysis of apple DFRs
To gain more insight into the evolutional and structural diversity of MdDFRs, the conserved motifs in MdDFR proteins were analyzed using the MEME online software.A total of 10 distinct and highly conserved motifs were captured (Fig. 3A).The motif distribution pattern in most MdDFR proteins is highly conserved.Clade I mainly contain motif 1, motif 2, motif 3, motif 4, motif 7 and motif 8. Clade II and III mainly contains motif 7 and motif 10.To further explore the structural diversity of MdDFRs, the intron-exon organization of each gene was analyzed.As shown in Fig. 3B, the exon number varied from 1 to 12, and most of the MdDFR genes have 1-7 exons.In addition, 6 genes contain 9 exons, 3 genes contain 8 exons, 4 genes contain 10 exons, and 4 genes contain 12 exons.

Prediction of cis-component analysis of MdDFRs
To explore the mechanism of MdDFR genes in plant growth and development, the cis-regulatory elements in the 2 kb upstream sequence of MdDFR genes were predicted by PlantCARE (Fig. 4).In the MdDFR genes' promoter, a variety of hormones (ABA, Auxin, SA, MeJA, GA, ERF)-, stress (heat shock, drought, low temperature, defense and stress)-, MYB binding site involved in flavonoid biosynthetic genes regulation-, light responsiveness-and lignin/flavonoid biosynthesis-related cis-acting elements were identified.Among them, light responsive related cis-acting elements were the most at 751

Intraspecific and interspecific colinearity analysis of DFRs
The internal collinearity analysis of the apple MdDFR genome revealed a total of 34 pairs of duplicated gene pairs (Fig. 5A), which were located on chromosomes Fig. 2 Phylogenetic analysis of the MdDFRs.Phylogenetic trees were constructed using the DFR protein sequences.NJ method was adopted, and the bootstrap value was set to be equal to 1000 chr1, chr2, chr3, chr4, chr5, chr6, chr7, ch8, chr9 and chr13.Most of these duplicated gene pairs were located on chromosomes 2 and 3 with 8 and 9 pairs, respectively.This indicates that some MdDFR genes may be generated through gene duplication, which may have similar functions.However, no linked genes were found on chromosomes 10, 11, 12, 14, 15 and 16.To further elucidate the phylogenetic mechanism and homology relationship of the DFR genes in apple, synteny analysis was conducted between apple and two representative species, grape (Vitis vinifera L.) and tomato (Solanum lycopersicum).
The results showed that there were 66 pairs of homologous genes in apple and tomato, mainly concentrated in apple chr3, chr7 and chr11, with 8, 10 and 8 pairs of homologous genes, respectively (Fig. 5B, Dataset1).There were 57 pairs of homologous genes in apple and grape, which were mainly concentrated in apple chr3, chr7 and chr11, with 7, 7 and 9 pairs of homologous genes, respectively.Therefore, it is speculated that the DFR gene is relatively conservative in evolution.In addition, some DFR genes were identified to be associated with 2-3 pairs of orthologs, whereas some had no homologous pairs, which may be due to differences in evolutionary speed and genetic specificity between species.

Expression analysis of DFR gene members in different apple cultivars
Based on the above results, 20 aforementioned MdDFRs were selected for expression level analysis in the fruit skin of four different colored apple varieties.Our previous research revealed that the degree of pigmentation and anthocyanin content in the fruit skin from G1 to G4 apple varieties were gradually increased [17].The results of qRT-PCR analysis in this research showed that the expression levels of MdDFR01, MdDFR05, MdDFR31, MdDFR53, MdDFR64, MdDFR69, MdDFR73, MdDFR84, MdDFR87, MdDFR94 and MdDFR96 were up-regulated with the increase of fruit skin pigmentation and the accumulation of anthocyanins (Fig. 8).Among these 20 genes, the expression level of MdDFR87 was the highest.Besides, the expression levels of MdDFR06, MdDFR15, MdDFR24, MdDFR46, MdDFR50 and MdDFR78 were significantly higher in G3 than those of other three apple varieties.However, the expression levels of MdDFR34, MdDFR65, and MdDFR81 were all lower, and the expression levels in G1 and G3 were higher than those in G2 and G4.

Discussion
After being extracted from maize and snapdragon using the transposon tagging approach [20], the DFR gene was later discovered in other plants [21][22][23].The 96 gene members of DFR were categorized into three subclades (Fig. 2, 3), and there were 66 and 57 pairs of homologous genes in apple/tomato and apple/grape (Fig. 5b), respectively, exhibiting significant sequence conservation in apple.DFR is a structural gene that regulates the phenylalanine metabolic pathway and has a significant impact on plant coloring [24].In transgenic tobacco lines, overexpression of Hosta ventricose gene HvDFR enhanced anthocyanin accumulation and elevated critical genes, such as F3H, F3'H, ANS, and UFGT [25].Heterologous expression of Ophiorrhiza japonica gene OjDFR1 in transgenic tobacco contributed to ).The gray lines in the background show collinearity between the apple and tomato/grape genomes.The blue lines show collinearity between the MdDFRs and tomato/grape genome darker flower color via up-regulating the expressions of endogenous NtANS and NtUFGT [26].Physicochemical property analysis showed that most MdDFR genes were acidic proteins, followed by basic proteins and neutral proteins (Table 1).Subcellular predictions found that 47 MdDFRs were located in the chloroplast,  31 MdDFRs were located in the cytoplasm (Table 1), which is consistent with the conclusion that anthocyanins are synthesized mainly in the cytoplasm [27].Therefore, we speculate that DFR genes are involved in the synthesis of apple anthocyanins.
Gene duplication plays an important role in the evolution of organisms, including tandem duplication, local duplication and whole genome duplication [28].Duplicate genes may be maintained through subfunctionalization and neofunctionalization at the level of expression or coding sequence [29].Although 96 DFR genes were distributed on 17 apple chromosomes (Fig. 1, Table 1), chr2 and chr11 each had 12 genes, followed by chr7 with 9 genes, chr6 and chr8 with 8 genes each, which was speculated to be formed by tandem repeats.Internal collinearity analysis revealed a total of 34 DFR duplicate gene pairs, mostly located on chromosomes 2 and 3, with 8 and 9 pairs, respectively (Fig. 5a), suggesting that these genes may have similar functions.This suggests that tandem duplication may be the main reason for the amplification of DFR family in apple.
The different colors of flowers and fruits of higher plants are closely related to plant pigments, which are anthocyanins and carotenoids.Anthocyanin transport and accumulation are regulated by many factors, mainly internal factors, but also by environmental factors.DFR is also regulated by environmental factors, and many plants also synthesize anthocyanins to protect the organism under environmental stress, which is accompanied by the increase of DFR expression in vivo [22].Previous studies have shown that in most plants, anthocyanin synthesis requires light induction [30][31][32].Under the effect of strong light, the expression of structural genes and regulatory genes will be enhanced [33].Natural variations of DFR were analyzed using an expanded population of 282 accessions belonging to the spiny Solanum group, and the single-nucleotide polymorphisms were found at the MYB recognition site in the promoter region, which causes differential expression of DFR and affects anthocyanin accumulation in fruits of the detected accessions [27].In the present research, lignin/flavonoid biosynthesis-related cis-acting elements were present in almost all MdDFR promoter sequences (Fig. 4).The past decades have seen remarkable inroads made into our understanding of MYB proteins and their binding sites more specifically.The MBSI (aaaAaaC(G/C)GTTA) and MBSII (aaaAGT TAG TTA) were the binding sites of Ph3 in Petunia hybrid and involved in flavonoid biosynthetic genes regulation [34] , the expression level of MdDFR54 and MdDFR91 was higher in apple flowers, and MdDFR64 was higher in ripe apple fruit peels and flowers, whereas lower in other tissues (Fig. 6).The expression levels of MdDFR05/31/64/69/94 were up-regulated with the increase of fruit skin pigmentation verified by qRT-PCR (Fig. 8).Therefore, the position and number of cis-acting elements related to anthocyanin synthesis in DFR promoter region are the causes for the changes in apple anthocyanin accumulation.
The biosynthetic pathway of anthocyanins has been defined, but its regulatory network still needs to be perfected.MYB plays an important role in promoting plant growth and development, secondary metabolism regulation, hormone regulation, and stress response [35][36][37][38].These transcription factors can recognize and bind to structural gene promoters to enhance or inhibit anthocyanin synthesis [39].In tree peony (Paeonia suffruticosa Andr.), PsMYB44 negatively regulates anthocyanin biosynthesis by directly binding to the PsDFR promoter and subsequently inhibiting blotch formation, which helps to elucidate the molecular regulatory network of anthocyanin-mediated blotch formation in plants [40].The ectopic combined expression of transcription factors MYB, bHLH, and WDR in Lotus japonicus indicated that the DFR2 promoter was activated [41].MdMYB1 appears to coordinately regulate structural genes in the anthocyanin pathway and segregates with red-skinned apples [42].MdbHLH33 can directly binds to the MdMYBPA1 promoter and interacts with it to promote anthocyanidin biosynthesis in response to low temperatures red-fleshed apples [12].Recent study shows that mdm-miR858 targets transcription factor genes MdMYB9 and MdMYBPA1 to participate in anthocyanin accumulation in apple [43].In addition, F3H and ANS, as upstream and downstream structural genes of DFR, respectively, are crucial in the phenylalanine metabolic pathway [6,7].In the present research, 75 MdDFR proteins interacted directly or indirectly with anthocyanidin synthesis related proteins MdANS, MdF3H, MdMYB1 and MdMYBPA1 to form a protein interaction network (Fig. 7).Besides, MdDFR69 and MdDFR87 had direct interactions with these four proteins, MdDFR64 had direct interactions with MdANS and MdF3H (Fig. 7), and all of these three genes were up-regulated with the increase of fruit skin pigmentation (Fig. 8), suggesting that the expression and regulation of these genes and transcription factors were coordinated with the production of anthocyanin.Therefore, further research on the transcription factors regulated by DFR upstream and downstream is of great significance to improve the regulatory network of anthocyanidin synthesis.

Conclusions
In this study, a total of 96 members of the DFR genes were identified from the apple database and were mainly expressed in chloroplasts and cytoplasm.Lignin/flavonoid synthesis-related elements of MBSI, MYB recognition site and MYB-binding site were identified in the promoters of MdDFR05/09/ 13

Fig. 3
Fig. 3 Motif and gene structure analysis of MdDFRs.A Analysis of conserved motif of DFRs in apple.B The exon-intron structure of MdDFRs

Fig. 4
Fig. 4 Cis-regulatory element analysis of the MdDFRs.Different numbers represent the number of cis-regulatory elements, and the color from blue to red represents the number of elements from fewer to more

Fig. 5
Fig. 5 Collinearity analysis of DFRs.A Collinearity analysis of MdDFRs.B Collinearity analysis of DFRs in apple and two representative plants tomato (Solanum lycopersicum) and grape (Vitis vinifera L.).The gray lines in the background show collinearity between the apple and tomato/grape genomes.The blue lines show collinearity between the MdDFRs and tomato/grape genome

Fig. 6
Fig. 6 Expression of MdDFRs in different tissues of apple.Yellow and blue shading represented the up-regulated or down-regulated expression level, respectively.The scale denoted the relative expression level

Fig. 7
Fig. 7 Analysis of protein interaction of 96 DFRs, ANS (DQ381771.1),F3H (AF117270.1),MYB1 (JN315783.1)and MYBPA1 (KJ909759.1) in apple.Nodes indicate proteins.Empty nodes indicate the protein of unknown 3D structures, and filled nodes indicate that 3D structures are known or predicted.The connection between nodes indicates the interaction between proteins, and different colors correspond to different types of interaction

Table 1
Analysis of physicochemical properties and subcellular localization of DFR gene in Apple