#These authors contribute equally.
Rosa rugosa, as a deciduous shrub in the Rosacea family, has high ornamental value and fragrant aroma, and plays an important role in landscape architecture. But the color of R. rugosa is simple, most is red, pink and white, and is fresh in other color, which have seriously restricted its application in gardens for a long time. In higher plants, the color of petals, tissues and fruits is determined by anthocyanin, which is a secondary metabolite belonging to flavonoids. Anthocyanins are synthesized through three steps with a series of enzymes  . The anthocyanins pathway has been explored in depth in Petunia hybrid  , Zea mays  , Antirrhinum majus  and so on. In higher plants, most anthocyanin biosynthesis is regulated by the combination of many transcription factors in different time and space. TFs like R2R3-MYB, bHLH, WD-40 can regulate anthocyanin biosynthesis by forming a MBW ternary complex   . As a key transcription factor, the up-regulation of R2R3-MYB can activate the expression of one or more structural genes involved in anthocyanin synthesis, so leads to anthocyanin accumulating. The related genes have been cloned and analyzed in its function from species such as Epimedium sagittatum  , Pyrus pyrifolia  and Prunus avium L.  . At present, there are few reports on the anthocyanin biosynthesis mechanism of R. rugose, so we don’t know exactly how it works. Based on the rose transcriptome data, we cloned one R2R3-MYB gene related with anthocyanin synthesis and figured out its function by bioinformatics analysis and real-time PCR to provide foundation for the regulation mechanism of anthocyanin synthesis in the petals of R. rugosa.
2. Materials and Methods
The experiment was conducted from April 2016 to January 2017 at the flower germplasm resource nursery of Shandong Agricultural University and the Flower Institute of Forestry College.
2.1. Plant Materials
The plant material, Chinese representative R. rugosa “Zi zhi”, “Bai zizhi”, “Fen zizhi” were from the rose germplasm resources garden at Shandong Agricultural College. R. rugosa “Zi zhi”, “Bai zizhi”, “Fen zizhi” are the most representative traditional rose in China. And the colors of three are purple, white and pink, respectively. From April to May, 2016, the petals in initial opening stage of a single plant with strong growth potential and stable hues in above three cultivars were collected for the experiment of gene cloning and difference expression; Sepals, stems, leaves, stamens, pistils and petals samples for gene spatial expression analysis were collected from R. rugosa “Zi zhi” at full opening stage. All samples were collected directly frozen with liquid nitrogen, and finally stored at −80˚C until used.
2.2.1. RNA Extraction and Gene Cloning
The extraction of RNA from various tissues of plants was operated according to the instructions of EASYspin plant RNA rapid extraction kit, and its concentration and purity were determined by ultraviolet spectrophotometer. At the same time, the integrity of RNA was detected by 1% agarose gel electrophoresis. cDNA was synthesized by RNA reverse transcription kit according to the description of abm reverse transcription kit. Fragment amplification: according to the data of rose transcriptome, the relevant MYB gene fragments were selected, and specific primers were designed by Oligo 7.0 software (Table 1). PCR amplification was carried out with the synthesized cDNA as the template. The reaction system of PCR was as follows: Mix 12.5 ul, the target gene upstream and downstream primers each 1 ul, template cDNA 1 ul, sterilization ddH2O supplement to 25 ul; The PCR reaction conditions are as follows: 94˚C for 5 min; 94˚C for 30 s, 53˚C for 30 s, and 72˚C for 1 min for a total of 35 cycles; and then extension at 72˚C for 10 min. Next, 1% agarose gel electrophoresis was used to detect the PCR products. The full-length amplification was carried out at 72˚C for 10 min. Specific primers (Table 1) were used to amplify the three races and the complete open reading frame (ORF). The reaction system and PCR reaction conditions were described above. The PCR product was reclaimed according to the description of Hipure Gel Pure DNA Mini Kit(Magen), then ligated with the carrier pMD18-T to transform Escherichia coli DH5α. Positive clones were selected and sequenced at BGI.
2.2.2. Bioinformatics Analysis of RrMYB6
BLASTX (NCBI) was used to study the homology of the nucleotide sequence and the deduced amino acid sequence. DNAMAN5.2.2 was used to conduct multiple sequence alignment. The ORF finder (NCBI) was used to search for an
Table 1. Primers used to clone and expression analysis of RrMYB6 in R. rugosa.
open reading frame, and the Conserved Domains database (NCBI) was used to analyze the conserved domains. ExPaSy-SOPMA was used to predict protein secondary structure. The ProtParam Tool was used to analyze protein physical and chemical properties. Furthermore, the ProtScale was used to predict hydrophilic or hydrophobic protein proper-ties. The NetPhos 3.1 Server was used to predict potential protein phosphorylation sites, and the NetOGlyc 4.0 Server was used to predict potential protein glycosylation sites. In order to study the evolutionary relationship between RrMYB6 and MYB6 proteins in other species, we chose 13 MYB6 genes in different spices with high level of homology to RrMYB6 by BLAST, and constructed evolutionary tree with them by MEGA5.0 software. The phylogenetic tree was constructed according to the neighbor-joining method, and tested by bootstrap, which was repeated 1000 times.
2.2.3. Expression Analysis of RrMYB6 in Different Tissues and Different Cultivars
Using CFX96TM Real-Time System RT-qPCR instrument and SYBR®Premix Ex TaqTM kit as template for real-time fluorescence quantitative analysis. The primers and internal reference primers (Actinin) are shown in Table 1. The reaction system is: Mix10ul, cDNA 1 ul, 0.4 ul each of upstream and downstream primer, adding sterilizing ddH2O up to 20 ul. The reaction procedure was: predenaturation at 95˚C for 30 s; 95˚C for 30 s, 60˚C for 30 s for a total of 39 cycles. The cycle threshold (Ct) value for each PCR reaction was calculated. After completion of the amplification steps, the melting curve was determined for each analysis. Gene transcripts were quantified using the comparative Ct method, which compares the transcript level of the target gene with that of the reference gene. The gene was set up with three repeats and the experimental data were processed by 2−ΔΔCT. And the data analysis was graphed through software Prism 5.
3. Results and Analysis
3.1. Cloning and Sequence Analysis of RrMYB6
One R2R3-MYB transcription factor, RrMYB6, was cloned from the petals of Rosa rugosa “Zi zhi”. Its Genbank accession Nos: MG745778. The RrMYB6 gene has a full length of 1107 bp, an open reading frames (ORFs) with length of 921 bp (Figure 1), and encoding 306 amino acids.
Multiple sequence alignment analysis showed that the amino acid sequence encoded by RrMYB6 had five conserved regions, including R2R3-MYB domain, C 1 and C 2 suppressor motifs domain, PCCEK (N1) motif and the domain interacting with bHLH protein. It is inferred that RrMYB6 is a typical R2R3-MYB transcription factor (Figure 2). The evolutionary tree results revealed that RrMYB6 exhibited the highest homology to FvMYB6 in Fragaria vesca. And six MYB6 proteins from other plants in Rosaceae family, such as MdMYB6, PmMYB6 and so on, were shown to be closely related with RrMYB6, and
Figure 1. PCR amplification of RrMYB6. M:Marker; A: Intermediate ragment; B: 3’-RACE; C: Full-length fragment.
Figure 2. Multiple alignment of theRrMYB6 with other MYB TFs Notes: The red linear indicate the conserved R2-domain and R3-domain, the black linear indicate the conserved residuals interacting with bHLH proteins.
clustering in a clade. While RrMYB6 was relatively distant from MaMYB6, JrMYB6, and other MYB6 proteins in other families, which were clustered into another branch (Figure 3).
3.2. Bioinformatics Analysis of RrMYB6 Gene
The RrMYB6 gene encoded 306 amino acids. The protein had a molecular weight of 34690.97 Da, an isoelectric point of 8.74, and its formula was C1491H2368N452O470S17. There were 42 basic amino acids (Arg + Lys), 35 acidic amino acids, and 229 neutral amino acids. It belonged to unstable protein with an unstable index at 50.59, and it was also a hydrophilic protein with the total average hydrophobic index at −0.847. The secondary structure prediction results demonstrated that there was 99 α-helix, 145 random coil, 34 extended peptide
Figure 3. The phylogenetic tree derived from the alignment of amino acid sequences of RrMYB6 and other MYB6 TFs.
chain, and 28 β-turn. The phosphorylation site prediction results demonstrated that there were 19 Ser phosphorylation sites, 22 Thr phosphorylation sites, and 6 Tyr phosphorylation sites, thus speculating that it was related to phosphorylation regulation.
3.3. Expression Patterns of RrMYB6 in Different Tissues and Different Varieties
Relative expression of RrMYB6 was profiled using q-PCR in R. rugose “Zi zhi” leaves, sepals, stems, stamens, pistils, and petals to analyze its tissue specificity, while it was also detected in the petals of R. rugose “Zi zhi”, R. rugosa “Fen zizhi”, R. rugosa “Bai zizhi” to analyze its variety specificity. The expression analysis of RrMYB6 in different tissues showed that RrMYB6 expressed differentially among stems, leaves, stamens, sepals, pistils and petals. RrMYB6 was more abundant in stamens than stems, leaves, petals, sepals and pistils. The highest expression level of RrMYB6 was observed in stamens, while it expressed slightly in sepals, pistils and petals, and almost didn’t express in stems and leaves. In addition, the expression analysis of RrMYB6 in different varieties showed that the RrMYB6 decreased with the color deepening among the three cultivars, lowest in the petals of R. rugose “Zi zhi”, followed by R. rugosa “Fen zizhi”, and highest in R. rugosa “Bai zizhi” (Figure 4).
One MYB transcription factor, RrMYB6, was cloned, and the sequence analysis of its protein showed that it contained a R2R3 domain, C1 and C2 suppressor motifs domain, PCCEK (N1) motif and the domain interacting with bHLH protein. Furthermore, its amino acid sequence is highly conserved at N1 motif and R2R3 domains, thus belonging to R2R3-MYB protein in Sg4 subfamily       . Studies show that the C1 repression motif of Sg4 subfamily affect the secondary metabolism of plants by inhibiting the expression of cinnamic acid 4-hydroxylase gene   . The R2R3-MYB proteins in different
Figure 4. Relative expression levels of RrMYB6.
spices have different effects on anthocyanin synthesis. For example, the overexpression of FaMYB1 in tobacco made its flower color shallower with the anthocyanins decreased significantly  , while the MdMYB3 isolated from Malus domstica exhibited oppositely   . Evolutionary analysis showed that RrMYB6 shared the closet homology to FvMYB6, so it was speculated that RrMYB6 could regulate the anthocyanin synthesis in R. rugosa. Specifically, whether it exhibited in a positive way or negative way needs further discussion.
The bioinformatics analysis showed that in the secondary structure of RrMYB6 protein, α-helix and random coil both accounted for considerable proportion. However, the proportion of extended strand and β-turn was smaller. Studies have proved that α-helix plays an important role in the R motif of the MYB domain. Each R motif is generally composed of three α-helix, in which the α-helix at the secondary and third position will form the HTH structure. And then combine with the first R motif, further forming a HTH domain with a hydrophobic core. What’s more, the third alpha helices in R motif have a role of identifying DNA, which makes the MYB protein more specific. Therefore, it can be further determined that RrMYB6 gene belongs to R2R3-MYB type   . The irregular curl structure is very beneficial to the binding of cell and water molecule, and both belong to hydrophilic protein. It can be inferred that both of them may protect the plant from osmotic stress. In addition, RrMYB6 contains a certain phosphorylation site, which may be related to the regulation of phosphorylation.
The qRT-PCR results showed that RrMYB6 expressed slightly in the petals of R. rugose “Zi zhi”, followed by R. rugosa “Fen zizhi”, and highly expressed in R. rugosa “Bai zizhi”. Furthermore, RrMYB6 exhibited a characteristic low expression in petals while expressed at a high level in stamen in the different tissues of R. rugosa “Zi zhi”. Studies show that there are positive and negative mechanisms for anthocyanin biosynthesis regulation by MYB protein in plants  . The MdMYB isolated from Malus domestica regulates the synthesis of anthocyanin. Moreover, lightness is the essential factor in the regulation of biosynthetic gene by MdMYB. Overexpression of MdMYB10 isolated from leaves and pulp can increase anthocyanin accumulation in apple seedlings. What’s more, the overexpression of MdMYB16, MdMYB17 and MdMYB111 in tobacco could inhibit the activity of DFR promoter and thus inhibiting anthocyanin synthesis   . Li also find that MYB protein has positive and negative mechanisms for anthocyanin synthesis in Morus alba  . The expression of RrMYB6 gene decreased with the deepening of petal color. And in the expression analysis of R. rugosa “Zi zhi” tissues, the expression of petals was very low, which was consistent with the analysis of expression among varieties. The gene expressed slightly in the petals, which may caused by different flowering stages. These results suggested that the regulation of RrMYB6 gene on anthocyanin biosynthesis in R. rugosa was a negative regulation mechanism. The regulation of genes involved in this mechanism needs further verification.
In conclusion, one R2R3-MYB TF, RrMYB6, was isolated from R. rugosa and was found to be involved in regulating anthocyanin biosynthetic pathway. The results of this study provided important information on the anthocyanin synthesis of R. rugosa. In future work, we will test whether the overexpression of RrMYB6 leads to anthocyanin accumulation in Arabidopsis thaliana and Nicotiana tabacum.
This work was funded by Shandong Province Agricultural Engineering project of breeding ( No. 96).
 Schwinn, K.E., Boase, M.R., Bradley, J.M., et al. (2014) MYB and bHLH transcription Factor Transgenes Increase Anthocyanin Pigmentation in Petunia and Lisianthus Plants, and the Petunia Phenotypes Are Strongly Enhanced under Field Conditions. Frontiers in Plant Science, 5, 603.
 Ibraheem, F., Gaffoor, I., Tan, Q., et al. (2015) A Sorghum MYB Transcription Factor Induces 3-deoxyanthocyanidins and Enhances Resistance against Leaf Blights in Maize. Molecules, 20, 2388-2404.
 Schwinn, K., Venail, J., Shang, Y., et al. (2006) A Small Family of MYB-Regulatory Genes Controls Floral Pigmentation Intensity and Patterning in the Genus Antirrhinum. The Plant Cell, 18, 831-851.
 Ravaglia, D., Espley, R., Henry, K.R., et al. (2013) Transcriptional Regulation of Flavonoid Biosynthesis in Nectarine (Prunus persica) by a Set of R2R3-MYB Transcription Factors. BMC Plant Biology, 13, 68.
 Li, L., Ban, Z.J., Li, X.H., et al. (2012) Differential Expression of Anthocyanin Biosynthetic Genes and Transcription Factor PcMYB10 in Pears (Pyrus communis L.). Plos One, 7, e46070.
 Huang, W., Sun, W., Lv, H., et al. (2013) A R2R3-MYB Transcription Factor from Epimedium sagittatum Regulates the Flavonoid Biosynthetic Pathway. PLoS ONE, 8, e70778.
 Feng, S., Wang, Y., Yang, S., et al. (2010) Anthocyanin Biosynthesis in Pears Is Regulated by a R2R3-MYB Transcription Factor PyMYB10. Planta, 232, 245-255.
 Starkevic, P., Paukstyte, J., Kazanaviciute, V., et al. (2015) Expression and Anthocyanin Biosynthesis-Modulating Potential of Sweet Cherry (Prunus avium L.) MYB10 and bHLH Genes. PLoS ONE, 10, e0126991.
 Xu, Z.S., Feng, K., Que, F., et al. (2017) A MYB Transcription Factor, DcMYB6, Is Involved in Regulating Anthocyanin Biosynthesis in Purple Carrot Taproots. Scientific Reports 7, Article ID: 45324.
 Butelli, E., Licciardello, C., Zhang, Y., Liu, J., et al. (2012) Retrotransposons Control Fruit-Specific, Cold-Dependent Accumulation of Anthocyanins in Blood Oranges. The Plant Cell, 24, 1242-1255.
 Yamagishi, M., Shimoyamada, Y., Nakatsuka, T., et al. (2010) Two R2R3-MYB Genes, Homologs of Petunia AN2, Regulate Anthocyanin Biosyntheses in Flower Tepals, Tepal Spots and Leaves of Asiatic Hybrid Lily. Plant and Cell Physiology, 51, 463-474.
 Xu, Z.S., et al. (2014) Transcript Profiling of Structural Genes Involved in Cyanidin-Based Anthocyanin Biosynthesis between Purple and Non-Purple Carrot (Daucus carota L.) Cultivars Reveals Distinct Patterns. BMC Plant Biology, 14, 1583-1588.
 Wang, H.Z., Qu, H.Y., Zhou, T.T., et al. (2017) Cloning and Expression Analysis of Anthocyanin Biosynthesis-Associated DFR and MYB Genes in Calyx of Eggplant (Solanum melongena L.). Scientia Agricultura Sinica, 50, 2781-2792.
 Jin, H., Cominelli, E., Bailey, P., et al. (2000) Transcriptional Repression by AtMYB4 Controls Production of UV-Protecting Sunscreens in Arabidopsis. The EMBO Journal, 19, 6150-6161.
 Lin, W.K., Bolitho, K., Grafton, K., et al. (2010) An R2R3-MYB Transcription Factor Associated with Regulation of the Anthocyanin Biosynthetic Pathway in Rosaceae. BMC Plant Biology, 10, 50.
 Vimolmangkang, S., Han, Y., Wei, G., et al. (2013) An Apple MYB Transcription Factor, MdMYB3, Is Involved in Regulation of Anthocyanin Biosynthesis and Flower Development. BMC Plant Biology, 13, 176.
 Aharoni, A., De, V.C.H., Wein, M., et al. (2001) The Strawberry FaMYB1 Transcription Factor Suppresses Anthocyanin and Flavonol Accumulation in Transgenic Tobacco. The Plant Journal, 28, 319-332.
 Wang, L., Han, W., Xie, C., et al. (2013) Comparing the Acceptor Promiscuity of a Rosa hybrida Glucosyltransferase RhGT1 and an Engineered Microbial Glucosyltransferase OleD (PSA) toward a Small Flavonoid Library. Carbohydrate Research, 368, 73-77.
 Liu, X.F., Li, F., Yin, X.R., et al. (2013) Recent Advances in the Transcriptional Regulation of Anthocyanin Biosynthesis. Acta Horticulturae Sinica, 40, 2295-2306.
 Bendahmane, M., Dubois, A., Raymond, O., et al. (2013) Genetics and Genomics of Flower Initiation and Development in Roses. Journal of Experimental Botany, 64, 847-857.
 Albert, N.W., Davies, K.M., Lewis, D.H., et al. (2014) A Conserved Network of Transcriptional Activators and Repressors Regulates Anthocyanin Pigmentation in Eudicots. The Plant Cell, 26, 962-980.
 Takos, A.M., Jaffe, F.W., Jacob, S.R., et al. (2006) Light-Induced Expression of a MYB Gene Regulates Anthocyanin Biosynthesis in Red Apples. Plant Physiology, 142, 1216-1232.
 Lin, W.K., Micheletti, D., Palmer, J., et al. (2011) High Temperature Reduces Apple Fruit Colour via Modulation of the Anthocyanin Regulatory Complex. Plant, Cell & Environment, 34, 1176-1190.