Next Article in Journal
Molecular Regulation of α3β4 Nicotinic Acetylcholine Receptors by Lupeol in Cardiovascular System
Next Article in Special Issue
High-Throughput Sequencing and Expression Analysis Suggest the Involvement of Pseudomonas putida RA-Responsive microRNAs in Growth and Development of Arabidopsis
Previous Article in Journal
Inherited Genetic Mutations and Polymorphisms in Malignant Mesothelioma: A Comprehensive Review
Previous Article in Special Issue
Small RNA and Transcriptome Sequencing Reveals miRNA Regulation of Floral Thermogenesis in Nelumbo nucifera
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 12
10.3390/ijms21124326
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

BrmiR828 Targets BrPAP1, BrMYB82, and BrTAS4 Involved in the Light Induced Anthocyanin Biosynthetic Pathway in Brassica rapa

by
Bo Zhou
1,2,*,
Jingtong Leng
2,
Yanyun Ma
2,
Pengzhen Fan
2,
Yuhua Li
1,2,
Haifang Yan
1,2 and
Qijiang Xu
1,2,*
1
Key Laboratory of Saline-alkali Vegetation Ecology Restoration (Northeast Forestry University), Ministry of Education, Harbin 150040, China
2
College of Life Science, Northeast Forestry University, Harbin 150040, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(12), 4326; https://doi.org/10.3390/ijms21124326
Submission received: 8 May 2020 / Revised: 12 June 2020 / Accepted: 13 June 2020 / Published: 17 June 2020
(This article belongs to the Special Issue Molecular Bases of Stress Adaptation in Plants: The Contribute of Regulation by Small RNAs in Plant Development and Stress Response)
Browse Figures
Versions Notes

Abstract

:
Comprehensive research in various plants shows that the metabolic pathway of anthocyanin biosynthesis is affected by environmental factors and regulated by microRNAs through post-transcriptional regulation. In seedlings of Brassica rapa Tsuda, the accumulation of anthocyanin is induced by light. However, the roles of BrmiR828 in the light-induced synthesis of anthocyanin in Brassica rapa remain to be explored. Here, a primary transcript of BrmiR828 was identified to be located on the chromosomes of the A03 sub-genome. Five candidate MYB family genes were predicted as targets of BrmiR828 in the database of Brassica rapa (BRAD, V1.1) by using psRNATarget. The transcript abundance of mature BrmiR828 was reduced in seedlings of Brassica rapa Tsuda under blue light irradiation comparing with dark treatment. However, Real-time PCR showed the transcript level of the five candidate targets, Bra004162, Bra022602, Bra001917, Bra029113, and Bra039763 was up-regulated when the seedlings exposed to blue or UV-A light. Trans-acting siRNA gene 4 (BrTAS4) was also identified to have a higher transcript level under blue and UV-A light irradiation than that in dark treatment. RNA ligase mediated 5′amplification of cDNA ends (RLM-5′ RACE) showed that BrmiR828 can splice the mRNA of Bra039763, Bra022602, and BrTAS4 on binding sites. Phylogenetic analysis of candidate BrMYBs targets along with MYBs from Arabidopsis thaliana showed that Bra039763, Bra004162, Bra001917, Bra029113, and Bra022602 are classified to the same group with AtMYB75, AtMYB114, AtMYB90, AtMYB113, and AtMYB82 which are involved in the anthocyanin biosynthetic pathway. As a result, light-induced down-regulation of BrmiR828 can target BrTAS4, BrPAP1 (Bra039763), MYB82 (Bra022602) to negatively regulate their transcript levels leading to the accumulation of MYB transcription factors that positively regulate anthocyanin biosynthesis in light-exposed seedlings of Brassica rapa.

1. Introduction

Anthocyanins are widely distributed in the plant kingdom and affect the color of fruits and vegetables [1,2]. They take important roles in protecting against UV radiation, attracting insects for pollination, pathogen resistance, and drought tolerance [3,4,5]. Also, anthocyanins have been demonstrated to possess antioxidative activity and have potential benefits for human health [6,7]. The pathway of anthocyanin biosynthesis has been clarified and most of the genes involved in anthocyanin biosynthesis and regulation have been identified in model plants [8]. However, the post-transcriptional regulation especially miRNAs involved in the molecular regulation of anthocyanin biosynthesis is still an important research objective.
Light is an important environmental signal to control metabolic processes during plant growth and development [9,10]. Light signaling components such as PIF3, CYR1, UVR8, and HY5 have been demonstrated to be key light signal transduction factors involved in anthocyanin biosynthesis under different light conditions [11,12,13,14]. Blue and UV-A light have also been reported to induce the expression of genes involved in anthocyanin biosynthesis in Brassica rapa [15,16], Litchi chinensis [17], Lactuca sativa [18], Arabidopsis [12], and Pyrus [19]. miRNAs negatively regulate the expression of genes involved in light-regulated processes [20,21]. Analysis of UV-B regulated miRNAs showed stress-responsive cis-elements located in promoters of MIR genes are speculated to be involved in responses to abiotic and biotic stresses [22,23,24]. Also, HY5 was shown to bind the promoters of MIR156D, MIR402, MIR408, and MIR858A and regulate their expression in Arabidopsis [25]. Analysis of the upstream regions of UV-B regulated MIRNA genes showed several light-related motifs, such as G box, GT-1 site, I-box, CCAAT-box, and GATA-box exist in regulatory regions and consistent with the effect of light on these genes [22,26]. Then light strongly affects the transcription level of MIRs and further research needs to be conducted to explore the roles of miRNAs and their target genes in light-induced anthocyanin biosynthesis process.
Transcription factors take important roles in regulating the transcription of genes involved in anthocyanin biosynthesis. R2R3 MYB, basic helix-loop-helix (bHLH) and WD repeat (WDR) transcription factors can form a trimeric MYB -bHLH- WD40 (MBW) complex to activate the expression of late pathway genes including dihydroflavonol reductase (DFR), anthocyanidin synthase (ANS), and flavonol-3-glucosyltransferase (3GT) [27,28,29]. In Arabidopsis, four genes PAP1/MYB75, PAP2/MYB90, MYB113, and MYB114 encoded R2R3-MYB proteins with high sequence similarities have been identified to positively regulate anthocyanin biosynthesis in vegetative tissues [30,31]. Also, GL3, EGL3, and TT8, the members of the bHLH transcription factor family have been proved to regulate anthocyanin biosynthesis [32,33,34]. Moreover, TTG1, the member of the WD40 protein determines the regulation of anthocyanin biosynthesis in Arabidopsis [35]. However, recent research shows that MYB-LIKE2 (MYBL2) and Squamosa promoter binding protein-like 9 (SPL9) negatively regulate anthocyanin biosynthesis through interfering with the formation of the MBW complex [36,37]. Then the quantitative competition between positive MBW regulators and negative regulators MYBL2/SPL9 can determine the activation/repression of the expression of anthocyanin biosynthetic pathway genes [38].
In plants, miRNAs play critical regulatory roles by negatively regulating the expression of transcription factors and metabolism-related proteins that modulate growth, development, and abiotic stress responses at the post-transcriptional level [36,39,40,41]. Recently, miRNAs were reported to be involved in anthocyanin biosynthesis in Arabidopsis [25,36,42,43]. MiR156 has been demonstrated to negatively regulate the transcription of SPL9, the protein product of which can compete with bHLH proteins to interact with PAP1 (Production of Anthocyanin Pigment 1) and positively regulate the accumulation of anthocyanins in the inflorescent stem [36]. MiR408, a positive regulator of photomorphogenesis, can lead to anthocyanin accumulation in seedlings [25,42]. MiR858a has also been demonstrated to enhance anthocyanin biosynthesis in seedlings by inhibiting the expression of MYBL2 through translational repression [43]. However, in Solanum lycopersicum, miR858 negatively regulates the anthocyanin pathway by modulating SlMYB48-like and SlMYB7-like transcripts [44]. Moreover, anthocyanin biosynthesis is also reported to be negatively regulated by miR828 or TAS D4(−) (previously TAS4-siRNA81(−)) in Arabidopsis [45,46,47]. Then miRNAs seem to have dual regulating roles in anthocyanin biosynthesis in different plants.
MiR828 was initially identified in Arabidopsis and functioned as the negative regulator of anthocyanin biosynthesis [45,47]. MiR828 mediates the cleavage of TAS4 and initiates the production of TAS D4(−) which is predicted to target MYB113, MYB75, and MYB90 genes that are involved in anthocyanin biosynthesis [47]. MiR828 has also been predicted to directly target MYB113 [46,48]. Over-expression of miR828 results in reduced expression of MYB113, MYB75, MYB82, and MYB90 and anthocyanin accumulation in Arabidopsis [45]. In apple (Malus × domestica), miR828, miR858, and TAS4 D4(−) has also been identified in high-throughput sequencing and was found to target 81 MYB genes [39]. Our previous research has characterized the miRNAs in response to blue/UV-A light in Brassica rapa [49], but the function of miR828 involved in light-induced anthocyanin biosynthesis in Brassica rapa remains to be elucidated.
Turnip (Brassica rapa subsp. Rapa; Brassicaceae) is an important and popular cruciferous root vegetable. The seedlings of turnip cv. Tsuda can accumulate anthocyanin in hypocotyls after exposure to sunlight, blue light, and UV-A light [16,49]. Among UV-B, UV-A, blue, red, and far-red lights, UV-A specific induces anthocyanin biosynthesis in the hypocotyls of turnip [50]. This study aimed to determine the regulation of BrmiR828 and its targets in blue light and UV-A induced anthocyanin biosynthesis of Brassica rapa. Based on the data, we demonstrate that BrmiR828 targets the transcripts of BrPAP1, BrMYB82, and BrTAS4 by RLM-5′ RACE. The transcript abundance of BrMIR828 decreased while the expression of its targets BrPAP1, BrPAP2, BrMYB82, and BrTAS4 greatly increased after blue light induction, indicating that BrmiR828 negatively controls light-induced anthocyanin biosynthetic pathway by repressing the expression of BrPAP1, BrPAP2, BrMYB82, and BrTAS4 genes in Brassica rapa.

2. Results

2.1. Light-Induced Anthocyanin Biosynthesis in Outer Epidermal Cells of Hypocotyls

When the seedlings of Brassica rapa were exposed to sunlight, pigments could accumulate in the hypocotyls of seedlings and showed a light-induced photomorphogenic phenotype with de-etiolation, phototropism, cotyledon expansion, chlorophyll synthesis, and short hypocotyls. Also, the content of anthocyanin in light-induced seedlings is much higher than that in dark treated seedlings (Figure 1). To clarify the distribution of anthocyanin in epidermal cells of hypocotyls, a hypocotyl cross-section was cut and examined using a confocal microscope. The result showed that anthocyanins synthesize in the outermost layer cells of the light treated hypocotyls and no anthocyanins accumulate in outer epidermal cells of dark treated hypocotyls (Figure 2).

2.2. Sequence Analysis of Primary miR828 Transcripts and Targets Identification of BrmiR828 in Brassica rapa

Through high throughput sequencing, the sequence of BrmiR828 was obtained in Brassica rapa [49]. To further study the regulatory roles of BrmiR828 in the anthocyanin biosynthesis pathway, primary candidate miR828 transcripts were identified and the target genes of miR828 were also predicted for analysis. Three candidate primary miR828 transcripts were identified by homology cloning and sequences alignment showed they are respectively located in B. rapa chromosome A03, A08, and A01. The sequences of these primary miR828 transcripts are highly similar except that primiR828-A01 has a 13–14 bp sequence deletion compared with primiR828-A01 and primiR828-A08 (Figure 3) and they are all predicted to have a stem-loop structure by Mfold [51] (Figure S1). According to the predicted stem-loop structure, there should be two kinds of mature miR828 sequences, BrmiR828a and BrmiR828b. BrmiR828a (from A01 and A08) is identical with AtmiR828 from Arabidopsis, and BrmiR828b (from A03) has only one base different with either AtmiR828 from Arabidopsis or SlmiR828 from Solanum lycopersicum (Figure S2). Only BrmiR828b was identified by sequencing and the obtained mature BrmiR828 might derive from the strand of the pre-miRNA (from A03) stem which was processed in the nucleus by RNase complex. Then the mature miR828 from Brassica rapa was deduced to have the same function as miR828 from Arabidopsis or Solanum lycopersicum. The prediction of miR828 targets showed five candidate MYB transcripts and one protein kinase family gene are directly bound by miR828 and regulated by cleavage using psRNATarget [52] (Table 1). Moreover, the transcript of a non-coding RNA, Trans-acting siRNA gene 4 (TAS4) cloned from Brassica rapa was also predicted to be cleaved by miR828.

2.3. Phylogenetic Analysis of the Candidate Target MYB Genes in Brassica rapa with MYB Gene Family in Arabidopsis

To investigate the phylogenic relationship between the target MYB genes of miR828 in Brassica rapa and MYB family genes in Arabidopsis, an unrooted phylogenetic tree was constructed using Neighbor-Joining method, based on the alignment amino acid sequences of candidate BrMYBs with 132 AtMYBs (Table S3). The phylogenic tree showed that all the MYB family proteins were divided into seven subgroups and the target MYB genes of miR828 were all clustered into subgroup VI (Figure 4A). Moreover, in subgroup VI, five Arabidopsis MYB genes encoding AtMYB75, AtMYB113, AtMYB114, AtMYB90, and AtMYB82 have been identified to be involved in anthocyanin biosynthesis [38,45,53]. Phylogenetic analysis indicates that the BrMYBs clustered with already characterized AtMYBs have similar functions. Multiple sequence alignment of MYBs in subgroup VI showed two DNA-binding domains, R2 and R3 are conserved and each domain has three α-helices which are reported to be critical for contacting with DNA (Figure 4B) [54,55]. Furthermore, a conserved motif DLX2RX3LX6LX3R in the first two α-helices of the R3 domain was also reported to be involved in MYB-bHLH interactions (Figure 4B) [54,55,56,57]. Three of five BrMYBs (BrPAP1 (Bra039763), BrPAP2 (Bra004162 and Bra001917)) showed similarity to AtMYB75, AtMYB113, AtMYB114, and AtMYB90. Bra039763 had a substitution at Arg63 (R) by Lys (K). The other two BrMYBs (Bra022602, Bra029113) were similar to AtMYB82 and putative BrMYB82 (Bra022602) had a substitution at Ser71 (S) by Thr (T) in the region between R2 and R3 domains (Figure 4C). The frequency analysis of the most prevalent amino acids at each position within each repeat of the R2R3-MYB domains showed highly conserved triplet tryptophan (Trp, W) residues are located at positions 6, 26, and 46 of the R2 domain and 59, 78, and 97 of the R3 domain in MYBs of Brassica rapa. In the R3 repeat, the first tryptophan (Trp59) residue was also usually replaced by isoleucine or phenylalanine (Ile, I or Phe, F) (Figure 4C).

2.4. Light-Dependent Expression Analysis of BrMIR828 and Its Candidate Target Genes under Blue Light and UV-A

To explore the expression character of BrMIR828 and its candidate target genes responsive to dark, blue, and UV-A light, their transcription level was detected by real-time PCR. The results showed that a significant change in abundance of BrmiR828 after blue light treatment compared with the dark treatment (p < 0.05) (two-tailed t-test, the p-value is 0.0267, t = 4.303, df = 3), but under UV-A light treatment compared with the dark treatment the abundance did not differ significantly (p < 0.05) (p-value is 0.0602, t = 4.303, df = 3) (Figure 5). In the candidate targets of BrmiR828, the transcript level of Bra039763 had no obvious difference (p < 0.05) (p-value is 0.0660, t = 4.303, df = 3) under UV-A light compared with dark treatment, and the others showed both blue and UV-A light-induced expression. The mRNA level of predicted targets Bra001917, Bra029113, and Bra004162 induced by blue light and UV-A was over five times as much as that in the dark treatment. However, the transcript level of Bra039763 was only about a 1.5-fold change in blue and UV-A light compared with that in dark treatment. The transcription level of BrTAS4 was also greatly increased after blue and UV-A light induction (Figure 5). The reversed expression patterns of target genes and BrMIR828 suggest post-transcription level regulation occurs between them.

2.5. BrmiR828 Guided the Cleavage of BrPAP1, BrMYB82 and BrTAS4 to Inhibit Their Expression

To confirm the regulation between BrmiR828 and their candidate targets, RLM-5′ RACE was used to detect the binding site and cleavage site. Sequencing of the 5′RACE PCR clones revealed that the miR828-guided cleavage of target BrPAP1 mRNA (Bra039763), BrMYB82 mRNA (Bra022602), and BrTAS4 mRNA occur as expected. Also, the cleavage sites were identified between 10 and 11 bases from the 5′ ends of the binding region in the BrmiR828 sequence (Figure 6A) and the motif site of MYB was targeted by miR828 coding for helix 3 of specific R3 domain (Figure 6B). Then, a conserved binding site targeted by BrmiR828 verified through RLM 5′-RACE can be concluded (Figure 6C). Unfortunately, we have not detected the cleavage sites of other candidate targets such as Bra001917, Bra004162, and Bra029113. Further analysis showed BrTAS4 is located on chromosome 9 of Brassica rapa (Figure 6D). After cleavage by BrmiR828, small RNAs (TAS4 siRNAs) were identified (Table 2) by high throughput sequencing [49] and one of the TAS4 siRNAs or siRNA TAS D4(−) was also predicted to target Bra039763 (BrPAP1), Bra004162 (BrPAP2), and Bra001917 (BrPAP2) (Figure S3).

3. Discussion

In this study, the role of miR828 involved in light-induced anthocyanin biosynthesis was investigated. Numerous studies indicate that various abiotic factors such as light [66], temperature [67], nitrogen and phosphate deficiencies [68], and sucrose [69] induce anthocyanin biosynthesis. Light is one of the most important environmental factors regulating the biosynthesis of anthocyanins. For example, blue and UV-B light have also been reported to mediate miRNAs that regulate genes in maize and Arabidopsis [23,70]. Anthocyanins are often presented in the epidermis of flowers and fruits to attract pollinators and seed dispersers [71]. In the hypocotyl of Brassica rapa seedlings, the pigment also accumulated in the outer layer cells of the epidermis after light induction (Figure 2). The anthocyanin accumulation is related to light-mediated photomorphogenic responses through photoreceptors and they might take resistance roles to UV damage.
miRNAs are known to be involved in various multiple biological processes such as development, primary and secondary metabolism, and stress responses [36,39]. In Arabidopsis, miR828 acts as a negative regulator and inhibits anthocyanin biosynthesis at different developmental stages [45]. In tomato, anthocyanin content was also greatly reduced in overexpressed AtmiR828 transgenic plants [72]. AtmiR828 encoded by a single locus lies in chromosome 4 (AT4G27765) which can form a stem-loop hairpin. By high throughput sequencing [49], a 22 nt long mature miR828 was also identified with almost the same sequences to AtmiR828 in Brassica rapa. MIR genes encoding identical mature miRNAs can differ in gene structure and regulatory sequences and express in response to various environmental stimuli [73]. It is necessary to identify the primary transcripts and the locus of miR828 in Brassica rapa. Although, three candidate MIR loci in chromosome A03, A08, A01 were obtained (Figure 3B) and predicted to form a stem-loop structure (Figure S1), only the mature miR828 originated from chromosome A03 was detected by sequencing [49]. The candidate mature miRNA from chromosome A08 and A01 still need experimental verification. Also, the candidate loci in chromosome A01 and A08 might be pseudo MIR gene. This might be due to several chromosome rearrangements events, leading to the evolution of Brassica rapa (AA, 2n = 20) [74]. Then the mature miR828 originated from chromosome A03 was used to predict the targets of miR828.
Among the predicted target genes of miR828, five genes were annotated to encode MYB protein. Phylogenetic analysis showed that candidate target MYB genes in Brassica rapa are clustered with AtMYB75, AtMYB113, AtMYB114, AtMYB90, and AtMYB82 in Arabidopsis (Figure 4A). AtMYB75 also known as Production of Anthocyanin Pigment 1 (PAP1) have been identified to control anthocyanin biosynthesis in vegetative tissues [75]. Overexpression of AtMYB90 (PAP2), AtMYB113, and AtMYB114 can also lead to anthocyanin accumulation in leaves and seedlings of Arabidopsis [33,75]. In addition, MYB82 was demonstrated to function in the regulation of trichome development [76] and might also take a potential role in the regulation of anthocyanin biosynthesis [45]. The results of the phylogenetic analysis indicate that the candidate BrPAP1 (Bra039763), BrPAP2 (Bra001917 and Bra004162), and BrMYB82 (Bra029113) might also function as a regulator of anthocyanin biosynthesis in Brassica rapa. Among these BrMYBs, each of the R2 and R3 domains contains three helices (Figure 4B), and helices 1 and 2 of the R3 domain can form a bHLH binding domain to interact with bHLH proteins [57]. The third helix is relatively conserved and important for DNA binding activity and the residues within the region between R2 and R3 domains are related to protein-DNA binding ability [77,78]. The results indicate that these conserved MYB domains can interact with bHLH protein together with WD repeat transcription factors to form MBW complex to activate the transcription of anthocyanin biosynthetic pathway genes including DFR, ANS, and 3-GT just like AtMYB75 in Arabidopsis [27,38,79].
Co-expression analysis of miRNAs and target genes is essential to understand the regulation relationship between them. A negative correlation between BrmiR828 and its predicted targets Bra022602 (MYB82), Bra029113 (MYB82), Bra001917 (PAP2), Bra004162 (PAP2), and Bra039763 (PAP1) was detected in dark and blue light treated seedlings by qRT-PCR analysis (Figure 5) and suggested that BrmiR828 directly mediate their cleavages in mRNA level. The mechanism is different from miR858a regulating MYBL2 through repressing translation [43]. Although the expression of BrmiR828 and Bra039763 did not show obvious differences under UV-A light and dark treatment, their regulatory roles in light-induced anthocyanin biosynthesis could not be ignored. The final expression level of the miRNA and its target gene which may not be correlated has been reported to be related to spatial restriction, mutual exclusion, dampening, and temporal regulation [80]. We also found that the abundance of BrmiR828 was also lower than that of other miRNAs such as miR156, miR159 detected by high throughput sequencing, and qRT-PCR [49]. Additionally, miR828 was also detected with low abundance in Malus × domestica [39], Arabidopsis [45], Prunus persica [81], and Vitis [82], and the functional role of miR828 has been proved to regulate the development of Arabidopsis trichome and cotton fiber [83], the anthocyanin biosynthesis of Vitis [82] and Arabidopsis [45]. A Brassica rapa TAS4 ortholog, defined as BrTAS4, was also identified to be a candidate target of miR828 and the sequence of miR828 binding site is similar to that of TAS4 in Arabidopsis except for a substitution of C by T base (Figure 6D). Only 24 bp length siRNAs from BrTAS4 were identified in Brassica rapa [49] and in Arabidopsis, the siRNAs were also mostly 24 mers [48]. Moreover, the expression of BrTAS4 was increased with blue and UV-A light treatment and showed a negative correlation with the expression of BrmiR828 under dark and blue light treatment. TAS4-siR81(−), one of the small RNAs originated from TAS4, was shown to target PAP1, PAP2, and MYB113 [46]. Furthermore, miR828 and TAS4 are also positively regulated by PAP1, and a negative feedback loop regulates the expression of PAP1 and its homologs [47]. These results have shown that a potential autoregulatory mechanism of the expression of anthocyanin pathway genes through sRNAs including miR828, TAS4-siR81(−), and miR156/157 [38,49].
Through RLM-5′ RACE, the cleavage site between miRNA and its target can be identified experimentally [84]. In the present study, we have obtained the binding cleavage site of BrmiR828 in BrPAP1 (Bra039763), BrMYB82 (Bra022602), and BrTAS4 and confirmed the direct cleavage regulation of miR828 to them in vivo. However, the cleavage site of other candidate targets (Bra001917, Bra004162, and Bra029113) have not been obtained in our research. Extensive base-pairing between the miRNA and the mRNA was not always sufficient to induce cleavage, and there can be additional requirements for a RNA-induced silencing complex (RISC) complex to catalyze endonucleolytic cleavage [85]. Also, these candidate MYBs might be targeted by other miRNAs/siRNAs, such as miR159, miR858, and tasiRNAs which have been reported to regulate MYBs in Malus × domestica [39]. To elucidate the regulation of miR828 and its targets, comparing the phenotypic impact of the expression of the miRNA-resistant target to its wild-type counterpart can be adopted [86]. In Arabidopsis, miR828 targeting AtMYB75, AtMYB113, AtMYB82, and TAS4 regulates the expression of genes in the anthocyanin biosynthetic pathway including DFR and LDOX [45,47]. MiR828 and miR858 have also been proved to regulate MYB114 to accumulate anthocyanin in Vitis [82]. Our present results are consistent with these findings and demonstrate that the mechanism of miR828 involved in anthocyanin biosynthesis is also conserved in Brassica rapa. Here, we also confirmed that BrmiR828 directly mediates the cleavage of BrMYB82 in Brassica rapa by RLM-5′ RACE. In Arabidopsis, the elevated level of miR828 inhibits the transcription of MYB82, indicating a potential role of MYB82 in the regulation of anthocyanin biosynthesis [45]. However, direct involvement of MYB82 in anthocyanin biosynthesis has not been demonstrated in any plant species except that MYB82 functions in the regulation of trichome development in Arabidopsis [76]. The regulated role of BrMYB82 on anthocyanin biosynthesis might be indirect in the pathway and our previous research showed miR156/157-guided cleavage of target SPL9 and SPL15 [49] also indirectly control anthocyanin accumulation by disrupting the stabilization of the MYB-bHLH-WD40 MBW complex [36]. Then the function for target MYB82 of miR828 in the pathway of anthocyanin biosynthesis needs to be done in the future.
The function of BrmiR828 and its transcription factor targets involved in anthocyanin biosynthesis should be detected through over or knockdown expression of these genes in Brassica rapa (turnip) plants. However, the generation of transgenic Brassica rapa plants has been reported to be difficult [87]. Although different transformation efficiencies vary among different genotypes and species of Brassica, compared with other Brassica species, Brassica rapa is the most recalcitrant species to transform [88]. Nowadays, the genetic transformation of turnip (Brassica rapa subsp. rapa) which form tubers has not been overcome. Overexpressing SlMYB75 (ortholog of BrPAP1) has been reported to effectively induce the accumulation of anthocyanin in both vegetative and reproductive organs of tomato [89]. Then the agroinfiltration-mediated transient expression system in tomato can be used to investigate the function of BrPAP1. When BrPAP1 overexpressed in tomato fruits through transient expression, pigment accumulated and the expression levels of SlCHS, SlDFR, and SlANS genes which involved in anthocyanin biosynthesis also increased [90]. The research suggested that BrPAP1, the target of BrmiR828, is involved in anthocyanin biosynthesis in Brassica rapa. Also, the function of BrMYB82 and BrTAS4 needs to be further explored by overexpression or knockdown analysis in the pathway of anthocyanin synthesis and metabolism in the future.

4. Materials and Methods

4.1. Plant Materials

Based on the previous research, we have known that blue light and UV-A light spectra were effect to accumulate anthocyanin in sunlight [16]. Then seeds of the turnip Brassica rapa L. subsp. rapa cv. Tsuda were sown in a row on wet filter paper and grown in the dark at 25 °C for 3 days, then irradiated with blue light (470 nm light-emitting diode LED, NSPB5205, Nichia, Tokushima, Japan) at 10 W m-2 or UV-A light (UV-A fluorescent lamp, FL10BLB, Toshiba, Japan, filtered through soda-lime glass plates, peak at 350 nm) at 3 W m-2 for 24 h. For the control group, the seedlings were still grown in the dark. Then light treated and dark-grown seedlings were used as materials for RNA isolation and real-time PCR. For anthocyanin measurements and cross-section analysis, the seedlings grew in sunlight or darkness for four days separately.

4.2. Anthocyanin Measurements

Approximate 0.1 g of four-day-old seedlings was weighed and ground in liquid nitrogen, and the pigments were extracted in 1 mL of 1% HCl in methanol (v/v) overnight at 4 °C with continuous shaking. For removing chlorophylls, 1 mL chloroform and 0.5 mL of water were added into the extracts and mixed to centrifuge at 15,000× g for 2 min. The top layer was collected and the quantity of anthocyanins was determined by spectrophotometric measurement of the aqueous phase (A530-A657) (Varian, Cary 50 Bio, Melbourne, Australia), and normalized to the total fresh weight of tissue used in each sample. (Modified after Martin et al. 2002). Eight biological replicates were prepared for light and dark treatment in anthocyanin analysis. Anthocyanin was determined by calculating the ratio of Abs530 to gram fresh mass [49].

4.3. Fluorescent Microscopic Observation

The hypocotyls of seedlings grown in dark and sunlight were used to obtain tissue slices. The tissue sections were made by a hand-sliced method, and put on a glass slide covered with a cover slip. Then the crosscut slices of hypocotyls were observed with a confocal laser scanning microscope (510 META, from Carl Zeiss, Jena, Germany). The argon laser with output at 488 nm unit was set for anthocyanin excitation and an emission filter 530 nm was used to detect the emission spectrum of fluorescence.

4.4. Prediction Potential Target Genes of BrmiR828

The potential targets of BrmiR828 were predicted by the psRNA Target program [52] using default parameters except Expectation set 3.0. Brassica rapa, de novo scaffolds assembly v1.1 (2011–08–30) cds were used as the cDNA library for the target search. The functional annotations of these target genes were then analyzed using the BLASTn search in the NCBI database.

4.5. WebLogo and Phylogenetic Analysis

The distribution of conserved R2R3-type MYB amino acid residues at the corresponding positions in domain profiles were generated using the Weblogo program with default parameters (http://weblogo.berkeley.edu/logo.cgi, v2.8.2) [65]. The alignments of full-length amino acid sequences of MYBs from Arabidopsis and Brassica rapa (Table S3) were aligned using ClustalX 2.1 [62] with default settings. Then the phylogenetic tree was constructed based on the alignments using MEGA 7.0 [61] with the Neighbor-joining method [58]. The parameters used in the tree construction were the Poisson model [60] plus uniform rates and 1000 bootstraps [59].

4.6. RNA Extraction and Quality Detection

Six seedlings in each sample with UV-A, blue light, and dark treatment were collected and isolated total RNA using TRNzol-A+ reagent according to the manufacturer’s instructions (TIANGEN Biotech, Beijing, China). The concentration and quality of total RNA were detected using 1% agar gel electrophoresis and NanoDrop 2000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, NC, USA).

4.7. Quantitative Real-Time PCR

Total RNA of each sample was extracted from six seedlings of Brassica rapa under dark, blue light or UV-A light treatment. For cDNA synthesis, 2 µg of total RNA of each sample was reverse transcribed with stem-loop primers for miRNAs and specific reverse primers for the housekeeping and target genes using Reagent Kit with gDNA Eraser (Takara, Dalian, China). All primer sequences are listed in Table S1. The qRT-PCR was performed using the POWER SYBR GREEN PCR Master Mix (Applied Biosystems, Foster City, CA, USA) and ABI 7500 real-time system (Applied Biosystems). Each reaction with a total volume of 20 µL contained 10 µL of PCR Master Mix, 1 µL of the first-strand cDNA as a template and 0.5 µM of each primer. The amplification program consisted of 95 °C for 30 s, followed by 40 cycles of 95 °C for 15 s, 62 °C for 35 s, and 72 °C for 35 s, then added melt gradient-dissociation curves. Three light and dark treatment replicates and three qRT-PCR experimental replicates were performed for each of the detected genes and samples. The fragment of the UBQ gene was amplified as an internal reference control for the normalization of data. A comparative CT (ΔΔCT) method was used for the analysis of the relative amounts of the transcripts [91]. Student’s t-tests were used to analyze the data, and a difference was considered to be statistically significant when p < 0.05, the data of three biological replicates was used for statistical analysis.

4.8. 5′ RNA Ligase-Mediated RACE (RLM-5′ RACE) PCR

Total RNA (3 µg) from different light-treated seedlings was mixed equally and used for RLM-5′ RACE analysis. The 5′ RACE adapter was directly ligated to mRNA according to the manufacturer’s instructions of the 5′-Full RACE Kit (TaKaRa, Dalian, China) without the alkaline phosphatase and tobacco acid pyrophosphatase treatment. The nested 5′ RACE Outer/Inner primer and two or three gene-specific nested primers were used for the rapid amplification of cDNA ends (Table S2). Then the PCR products were purified and cloned into the pEASY-T5 zero cloning vector (TransGen Biotech, Beijing, China) for sequencing.

5. Conclusions

We show here that during anthocyanin accumulating in light-induced seedlings of Brassica rapa, the transcript level of BrPAP1, BrMYB82, and BrTAS4 is positively related to the content of pigment and BrPAP1 has been proved to be involved in anthocyanin biosynthetic pathway. BrmiR828 directly splices mRNA of BrPAP1, BrMYB82, and BrTAS4 and negatively regulates their expression in seedlings of light treatment. Therefore, BrmiR828 targets BrPAP1, BrMYB82, and BrTAS4 involved in light-induced anthocyanin biosynthetic pathway in Brassica rapa. The current findings supply a clue to facilitate the production of rich anthocyanin in plants for food with potential benefits to human health.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/12/4326/s1, Figure S1: The prediction stem-loop structures of premiR828 homologs isolated in Brassica rapa by mfold, Figure S2: Comparison of miR828 sequences in turnip, Arabidopsis and tomato, Figure S3: The complementarity binding profiles for BrmiR828 with BrTAS4 and for TAS4-D4 (−) with Bra001917 (BrPAP2), Bra039763 (BrPAP1) and Bra004162 (BrPAP2), Table S1: Primers for premiR828, TAS4 and Real-time PCR primers for detected miR828 and targeted genes, Table S2: Nest primers for RLM 5′-RACE PCR used in this topic. Table S3: Full-length amino acid sequences of MYBs from Arabidopsis and Brassica rapa.

Author Contributions

Conceptualization, B.Z. and Q.X.; methodology, B.Z.; validation, B.Z., H.Y., and Q.X.; formal analysis, B.Z.; investigation, J.L., Y.M., P.F., H.Y., and B.Z.; resources, Y.L.; data curation, X.X.; writing—original draft preparation, B.Z.; writing—review and editing, B.Z. and Q.X.; visualization, B.Z.; supervision, B.Z.; project administration, B.Z.; funding acquisition, Q.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities (2572017CA23, 2572014EA03-02) and Heilongjiang Provincial Natural Science Foundation of China(C2018002).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

3GTflavonol-3-glucosyltransferase
ANSanthocyanidin synthase
ArgArginine
AsnAsparagine
bHLHbasic helix-loop-helix
CYR1cryptochrome 1
DFRdihydroflavonol reductase
EGL3ENHANCER OF GLABRA 3
GL3GLABRA3
HY5Elongated Hypocotyl 5
LDOXleucoanthocyandin dioxygenase
LysLysine
MBWMYB-bHLH-WD40
miRNAsmicroRNAs
MYBL2MYB-LIKE2
PAP1Production of Anthocyanin Pigment 1
PIF3phytochrome-interacting factor 3
qRT-PCRquantitative reverse transcription PCR
RISCRNA-induced silencing complex
RLM-5′ RACERNA ligasemediated 5′amplification of cDNA ends
SerSerine
siRNASmall interfering RNA
SPLSQUAMOSA Promoter-Binding Protein-Like
SPL9SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9
sRNAssmall RNAs
TAS4Trans-acting siRNA gene 4
ThrThreonine
TrpTryptophan
TT8TRANSPARENT TESTA8
TTG1TRANSPARENT TESTA GLABRA1
UBQUbiquitin
UV-Aultraviolet A
UVR8ultraviolet-B receptor UVR8
WDRWD repeat

References

  1. Jaakola, L. New insights into the regulation of anthocyanin biosynthesis in fruits. Trends Plant Sci. 2013, 18, 477–483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Leong, S.Y.; Oey, I. Effects of processing on anthocyanins, carotenoids and vitamin C in summer fruits and vegetables. Food Chem. 2012, 133, 1577–1587. [Google Scholar] [CrossRef]
  3. Gould, K.S. Nature’s Swiss army knife: The diverse protective roles of anthocyanins in leaves. J. Biomed. Biotechnol. 2004, 2004, 314–320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tanaka, Y.; Sasaki, N.; Ohmiya, A. Biosynthesis of plant pigments: Anthocyanins, betalains and carotenoids. Plant J. 2008, 54, 733–749. [Google Scholar] [CrossRef] [PubMed]
  5. Buer, C.S.; Imin, N.; Djordjevic, M.A. Flavonoids: New Roles for Old Molecules. J. Integr. Plant Biol. 2010, 52, 98–111. [Google Scholar] [CrossRef]
  6. He, J.; Giusti, M.M. Anthocyanins: Natural colorants with health-promoting properties. Annu. Rev. Food Sci. Technol. 2010, 1, 163–187. [Google Scholar] [CrossRef]
  7. de Pascual-Teresa, S.; Moreno, D.A.; Garcia-Viguera, C. Flavanols and anthocyanins in cardiovascular health: A review of current evidence. Int. J. Mol. Sci. 2010, 11, 1679–1703. [Google Scholar] [CrossRef] [Green Version]
  8. Quattrocchio, F.; Wing, J.F.; Leppen, H.; Mol, J.; Koes, R.E. Regulatory genes controlling anthocyanin pigmentation are functionally conserved among plant species and have distinct sets of target genes. Plant Cell 1993, 5, 1497–1512. [Google Scholar] [CrossRef]
  9. Stracke, R.; Favory, J.J.; Gruber, H.; Bartelniewoehner, L.; Bartels, S.; Binkert, M.; Funk, M.; Weisshaar, B.; Ulm, R. The Arabidopsis bZIP transcription factor HY5 regulates expression of the PFG1/MYB12 gene in response to light and ultraviolet-B radiation. Plant. Cell Environ. 2010, 33, 88–103. [Google Scholar] [CrossRef]
  10. Krahmer, J.; Ganpudi, A.; Abbas, A.; Romanowski, A.; Halliday, K.J. Phytochrome, carbon sensing, metabolism, and plant growth plasticity. Plant Physiol. 2018, 176, 1039–1048. [Google Scholar] [CrossRef] [Green Version]
  11. Shin, J.; Park, E.; Choi, G. PIF3 regulates anthocyanin biosynthesis in an HY5-dependent manner with both factors directly binding anthocyanin biosynthetic gene promoters in Arabidopsis. Plant J. 2007, 49, 981–994. [Google Scholar] [CrossRef]
  12. Chen, D.Q.; Li, Z.Y.; Pan, R.C.; Wang, X.J. Anthocyanin accumulation mediated by blue light and cytokinin in Arabidopsis seedlings. J. Integr. Plant Biol. 2006, 48, 420–425. [Google Scholar] [CrossRef]
  13. Gu, K.D.; Wang, C.K.; Hu, D.G.; Hao, Y.J. How do anthocyanins paint our horticultural products? Sci. Hortic. Amst. 2019, 249, 257–262. [Google Scholar] [CrossRef]
  14. Shin, D.H.; Choi, M.; Kim, K.; Bang, G.; Cho, M.; Choi, S.B.; Choi, G.; Park, Y.I. HY5 regulates anthocyanin biosynthesis by inducing the transcriptional activation of the MYB75/PAP1 transcription factor in Arabidopsis. FEBS Lett. 2013, 587, 1543–1547. [Google Scholar] [CrossRef] [Green Version]
  15. Zhou, B.; Wang, Y.; Zhan, Y.; Li, Y.; Kawabata, S. Chalcone synthase family genes have redundant roles in anthocyanin biosynthesis and in response to blue/UV-A light in turnip (Brassica rapa; Brassicaceae). Am. J. Bot. 2013, 100, 2458–2467. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Y.; Zhou, B.; Sun, M.; Li, Y.; Kawabata, S. UV-A light induces anthocyanin biosynthesis in a manner distinct from synergistic blue + UV-B light and UV-A/blue light responses in different parts of the hypocotyls in turnip seedlings. Plant Cell Physiol. 2012, 53, 1470–1480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zhang, H.N.; Li, W.C.; Wang, H.C.; Shi, S.Y.; Shu, B.; Liu, L.Q.; Wei, Y.Z.; Xie, J.H. Transcriptome Profiling of Light-Regulated Anthocyanin Biosynthesis in the Pericarp of Litchi. Front. Plant Sci. 2016, 7, 963. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Zhang, Y.; Xu, S.; Cheng, Y.; Peng, Z.; Han, J. Transcriptome profiling of anthocyanin-related genes reveals effects of light intensity on anthocyanin biosynthesis in red leaf lettuce. Peerj 2018, 6, e4607. [Google Scholar] [CrossRef]
  19. Tao, R.Y.; Bai, S.L.; Ni, J.B.; Yang, Q.S.; Zhao, Y.; Teng, Y.W. The blue light signal transduction pathway is involved in anthocyanin accumulation in ‘Red Zaosu’ pear. Planta 2018, 248, 37–48. [Google Scholar] [CrossRef]
  20. Sun, W.; Xu, X.H.; Wu, X.; Wang, Y.; Lu, X.B.; Sun, H.W.; Xie, X.Z. Genome-wide identification of microRNAs and their targets in wild type and phyB mutant provides a key link between microRNAs and the phyB-mediated light signaling pathway in rice. Front. Plant Sci. 2015, 6, 372. [Google Scholar] [CrossRef] [Green Version]
  21. Sun, Z.F.; Li, M.; Zhou, Y.; Guo, T.T.; Liu, Y.; Zhang, H.; Fang, Y.D. Coordinated regulation of Arabidopsis microRNA biogenesis and red light signaling through Dicer-like 1 and phytochrome-interacting factor 4. PLoS Genet. 2018, 14, e1007247. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, B.; Sun, Y.F.; Song, N.; Wang, X.J.; Feng, H.; Huang, L.L.; Kang, Z.S. Identification of UV-B-induced microRNAs in wheat. Genet. Mol. Res. 2013, 12, 4213–4221. [Google Scholar] [CrossRef] [PubMed]
  23. Casati, P. Analysis of UV-B regulated miRNAs and their targets in maize leaves. Plant Signal. Behav. 2013, 8, e26758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Li, S.; Castillo-Gonzalez, C.; Yu, B.; Zhang, X. The functions of plant small RNAs in development and in stress responses. Plant J. 2017, 90, 654–670. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, H.Y.; He, H.; Wang, X.C.; Wang, X.F.; Yang, X.Z.; Li, L.; Deng, X.W. Genome-wide mapping of the HY5-mediated genenetworks in Arabidopsis that involve both transcriptional and post-transcriptional regulation. Plant J. 2011, 65, 346–358. [Google Scholar] [CrossRef]
  26. Zhou, X.F.; Wang, G.D.; Zhang, W.X. UV-B responsive microRNA genes in Arabidopsis thaliana. Mol. Syst. Biol. 2007, 3, 103. [Google Scholar] [CrossRef]
  27. Xu, W.J.; Dubos, C.; Lepiniec, L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 2015, 20, 176–185. [Google Scholar] [CrossRef]
  28. Broun, P. Transcriptional control of flavonoid biosynthesis: A complex network of conserved regulators involved in multiple aspects of differentiation in Arabidopsis. Curr. Opin. Plant Biol. 2005, 8, 272–279. [Google Scholar] [CrossRef]
  29. Stracke, R.; Ishihara, H.; Barsch, G.H.A.; Mehrtens, F.; Niehaus, K.; Weisshaar, B. Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 2007, 50, 660–677. [Google Scholar] [CrossRef] [Green Version]
  30. Stracke, R.; Werber, M.; Weisshaar, B. The R2R3-MYB gene family in Arabidopsis thaliana. Curr. Opin. Plant Biol. 2001, 4, 447–456. [Google Scholar] [CrossRef]
  31. Dubos, C.; Stracke, R.; Grotewold, E.; Weisshaar, B.; Martin, C.; Lepiniec, L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15, 573–581. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, F.; Gonzalez, A.; Zhao, M.; Payne, C.T.; Lloyd, A. A network of redundant bHLH proteins functions in all TTG1-dependent pathways of Arabidopsis. Development 2003, 130, 4859–4869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Gonzalez, A.; Zhao, M.; Leavitt, J.M.; Lloyd, A.M. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J. 2008, 53, 814–827. [Google Scholar] [CrossRef] [PubMed]
  34. Baudry, A.; Caboche, M.; Lepiniec, L. TT8 controls its own expression in a feedback regulation involving TTG1 and homologous MYB and bHLH factors, allowing a strong and cell-specific accumulation of flavonoids in Arabidopsis thaliana. Plant J. 2006, 46, 768–779. [Google Scholar] [CrossRef] [PubMed]
  35. Walker, A.R.; Davison, P.A.; Bolognesi-Winfield, A.C.; James, C.M.; Srinivasan, N.; Blundell, T.L.; Esch, J.J.; Marks, M.D.; Gray, J.C. The TRANSPARENT TESTA GLABRA1 locus, which regulates trichome differentiation and anthocyanin biosynthesis in Arabidopsis, encodes a WD40 repeat protein. Plant Cell 1999, 11, 1337–1350. [Google Scholar] [CrossRef] [Green Version]
  36. Gou, J.Y.; Felippes, F.F.; Liu, C.J.; Weigel, D.; Wang, J.W. Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor. Plant Cell 2011, 23, 1512–1522. [Google Scholar] [CrossRef] [Green Version]
  37. Matsui, K.; Umemura, Y.; Ohme-Takagi, M. AtMYBL2, a protein with a single MYB domain, acts as a negative regulator of anthocyanin biosynthesis in Arabidopsis. Plant J. 2008, 55, 954–967. [Google Scholar] [CrossRef]
  38. Shi, M.Z.; Xie, D.Y. Biosynthesis and metabolic engineering of anthocyanins in Arabidopsis thaliana. Recent Pat. Biotechnol. 2014, 8, 47–60. [Google Scholar] [CrossRef] [Green Version]
  39. Xia, R.; Zhu, H.; An, Y.Q.; Beers, E.P.; Liu, Z.R. Apple miRNAs and tasiRNAs with novel regulatory networks. Genome Biol. 2012, 13, R47. [Google Scholar] [CrossRef] [Green Version]
  40. Sun, C.; Zhao, Q.; Liu, D.D.; You, C.X.; Hao, Y.J. Ectopic expression of the apple Md-miRNA156h gene regulates flower and fruit development in Arabidopsis. Plant Cell Tiss. Org. 2013, 112, 343–351. [Google Scholar] [CrossRef]
  41. Jones-Rhoades, M.W.; Bartel, D.P.; Bartel, B. MicroRNAs and their regulatory roles in plants. Annu. Rev. Plant. Biol. 2006, 57, 19–53. [Google Scholar] [CrossRef] [PubMed]
  42. Zhang, H.Y.; Zhao, X.; Li, J.G.; Cai, H.Q.; Deng, X.W.; Li, L. MicroRNA408 Is Critical for the HY5-SPL7 Gene Network That Mediates the Coordinated Response to Light and Copper. Plant Cell 2014, 26, 4933–4953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wang, Y.L.; Wang, Y.Q.; Song, Z.Q.; Zhang, H.Y. Repression of MYBL2 by Both microRNA858a and HY5 Leads to the Activation of Anthocyanin Biosynthetic Pathway in Arabidopsis. Mol. Plant 2016, 9, 1395–1405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Jia, X.Y.; Shen, J.; Liu, H.; Li, F.; Ding, N.; Gao, C.Y.; Pattanaik, S.; Patra, B.; Li, R.Z.; Yuan, L. Small tandem target mimic-mediated blockage of microRNA858 induces anthocyanin accumulation in tomato. Planta 2015, 242, 283–293. [Google Scholar] [CrossRef]
  45. Yang, F.X.; Cai, J.; Yang, Y.; Liu, Z.B. Overexpression of microRNA828 reduces anthocyanin accumulation in Arabidopsis. Plant. Cell Tiss. Org. 2013, 115, 159–167. [Google Scholar] [CrossRef]
  46. Luo, Q.J.; Mittal, A.; Jia, F.; Rock, C.D. An autoregulatory feedback loop involving PAP1 and TAS4 in response to sugars in Arabidopsis. Plant Mol. Biol. 2012, 80, 117–129. [Google Scholar] [CrossRef] [Green Version]
  47. Hsieh, L.C.; Lin, S.I.; Shih, A.C.C.; Chen, J.W.; Lin, W.Y.; Tseng, C.Y.; Li, W.H.; Chiou, T.J. Uncovering Small RNA-Mediated Responses to Phosphate Deficiency in Arabidopsis by Deep Sequencing. Plant Physiol. 2009, 151, 2120–2132. [Google Scholar] [CrossRef] [Green Version]
  48. Rajagopalan, R.; Vaucheret, H.; Trejo, J.; Bartel, D.P. A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Gene Dev. 2006, 20, 3407–3425. [Google Scholar] [CrossRef] [Green Version]
  49. Zhou, B.; Fan, P.Z.; Li, Y.H.; Yan, H.F.; Xu, Q.J. Exploring miRNAs involved in blue/UV-A light response in Brassica rapa reveals special regulatory mode during seedling development. BMC Plant Biol. 2016, 16, 111. [Google Scholar] [CrossRef] [Green Version]
  50. Zhou, B.; Li, Y.H.; Xu, Z.R.; Yan, H.F.; Homma, S.; Kawabata, S. Ultraviolet A-specific induction of anthocyanin blosynthesis in the swollen hypocotyls of turnip (Brassica rapa). J. Exp. Bot. 2007, 58, 1771–1781. [Google Scholar] [CrossRef]
  51. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef] [PubMed]
  52. Dai, X.B.; Zhuang, Z.H.; Zhao, P.X.C. psRNATarget: A plant small RNA target analysis server (2017 release). Nucleic Acids Res. 2018, 46, W49–W54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Liu, J.Y.; Osbourn, A.; Ma, P.D. MYB Transcription Factors as Regulators of Phenylpropanoid Metabolism in Plants. Mol. Plant. 2015, 8, 689–708. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Gabrielsen, O.S.; Sentenac, A.; Fromageot, P. Specific DNA binding by c-Myb: Evidence for a double helix-turn-helix-related motif. Science 1991, 253, 1140–1143. [Google Scholar] [CrossRef]
  55. Jia, L.; Clegg, M.T.; Jiang, T. Evolutionary dynamics of the DNA-binding domains in putative R2R3-MYB genes identified from rice subspecies indica and japonica genomes. Plant Physiol. 2004, 134, 575–585. [Google Scholar] [CrossRef] [Green Version]
  56. Tominaga, R.; Iwata, M.; Okada, K.; Wada, T. Functional analysis of the epidermal-specific MYB genes CAPRICE and WEREWOLF in Arabidopsis. Plant Cell 2007, 19, 2264–2277. [Google Scholar] [CrossRef] [Green Version]
  57. Zhu, H.F.; Fitzsimmons, K.; Khandelwal, A.; Kranz, R.G. CPC, a single-repeat R3 MYB, is a negative regulator of anthocyanin biosynthesis in Arabidopsis. Mol. Plant 2009, 2, 790–802. [Google Scholar] [CrossRef]
  58. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef]
  59. Felsenstein, J. Confidence Limits on Phylogenies: An Approach Using the Bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef]
  60. Zuckerkandl, E.; Pauling, L. Evolutionary Divergence and Convergence in Proteins; Academic Press: New York, NY, USA, 1965. [Google Scholar]
  61. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  62. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [Green Version]
  63. Nicholas, K.B.; Nicholas, H.B.J. GeneDoc: Analysis and visualization of genetic variation. Embnew 1997, 4, 14. [Google Scholar]
  64. Li, C.L.; Lu, S.F. Genome-wide characterization and comparative analysis of R2R3-MYB transcription factors shows the complexity of MYB-associated regulatory networks in Salvia miltiorrhiza. BMC Genom. 2014, 15, 277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Crooks, G.E.; Hon, G.; Chandonia, J.M.; Brenner, S.E. WebLogo: A sequence logo generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [Green Version]
  66. Dubos, C.; Le Gourrierec, J.; Baudry, A.; Huep, G.; Lanet, E.; Debeaujon, I.; Routaboul, J.M.; Alboresi, A.; Weisshaar, B.; Lepiniec, L. MYBL2 is a new regulator of flavonoid biosynthesis in Arabidopsis thaliana. Plant J. 2008, 55, 940–953. [Google Scholar] [CrossRef] [PubMed]
  67. Rowan, D.D.; Cao, M.S.; Lin-Wang, K.; Cooney, J.M.; Jensen, D.J.; Austin, P.T.; Hunt, M.B.; Norling, C.; Hellens, R.P.; Schaffer, R.J.; et al. Environmental regulation of leaf colour in red 35S:PAP1 Arabidopsis thaliana. New Phytol. 2009, 182, 102–115. [Google Scholar] [CrossRef] [PubMed]
  68. Lea, U.S.; Slimestad, R.; Smedvig, P.; Lillo, C. Nitrogen deficiency enhances expression of specific MYB and bHLH transcription factors and accumulation of end products in the flavonoid pathway. Planta 2007, 225, 1245–1253. [Google Scholar] [CrossRef]
  69. Teng, S.; Keurentjes, J.; Bentsink, L.; Koornneef, M.; Smeekens, S. Sucrose-specific induction of anthocyanin biosynthesis in Arabidopsis requires the MYB75/PAP1 gene. Plant Physiol. 2005, 139, 1840–1852. [Google Scholar] [CrossRef] [Green Version]
  70. Jung, J.H.; Seo, Y.H.; Seo, P.J.; Reyes, J.L.; Yun, J.; Chua, N.H.; Park, C.M. The GIGANTEA-regulated MicroRNA172 mediates photoperiodic flowering independent of CONSTANS in Arabidopsis. Plant Cell 2007, 19, 2736–2748. [Google Scholar] [CrossRef] [Green Version]
  71. Rieseberg, L.H.; Blackman, B.K. Speciation genes in plants. Ann. Bot. Lond. 2010, 106, 439–455. [Google Scholar] [CrossRef] [Green Version]
  72. Jia, X.Y.; Yu, Z.Q.; Liang, J.P.; Tang, G.L.; Jin, L.H.; Zhang, L.; He, L.H.; Li, R.Z. Cloning of Arabidopsis At-pri-miR828 Gene and Its Genetic Transformation into Tomato. Acta Hortic. Sin. 2013, 40, 2419–2428. [Google Scholar]
  73. Zhao, X.; Zhang, H.; Li, L. Identification and analysis of the proximal promoters of microRNA genes in Arabidopsis. Genomics 2013, 101, 187–194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Waminal, N.E.; Perumal, S.; Lee, J.; Kim, H.H.; Yang, T.J. Repeat Evolution in Brassica rapa (AA), B. oleracea (CC), and B. napus (AACC) Genomes. Plant Breed. Biotech. 2016, 4, 15. [Google Scholar] [CrossRef] [Green Version]
  75. Borevitz, J.O.; Xia, Y.; Blount, J.; Dixon, R.A.; Lamb, C. Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 2000, 12, 2383–2394. [Google Scholar] [CrossRef] [Green Version]
  76. Liang, G.; He, H.; Li, Y.; Ai, Q.; Yu, D. MYB82 functions in regulation of trichome development in Arabidopsis. J. Exp. Bot. 2014, 65, 3215–3223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Dias, A.P.; Braun, E.L.; McMullen, M.D.; Grotewold, E. Recently duplicated maize R2R3 Myb genes provide evidence for distinct mechanisms of evolutionary divergence after duplication. Plant Physiol. 2003, 131, 610–620. [Google Scholar] [CrossRef] [Green Version]
  78. Heine, G.F.; Hernandez, J.M.; Grotewold, E. Two cysteines in plant R2R3 MYB domains participate in REDOX-dependent DNA binding. J. Biol. Chem. 2004, 279, 37878–37885. [Google Scholar] [CrossRef] [Green Version]
  79. Li, Z.; Peng, R.; Tian, Y.; Han, H.; Xu, J.; Yao, Q. Genome-Wide Identification and Analysis of the MYB Transcription Factor Superfamily in Solanum lycopersicum. Plant Cell Physiol. 2016, 57, 1657–1677. [Google Scholar] [CrossRef] [Green Version]
  80. Voinnet, O. Origin, biogenesis, and activity of plant microRNAs. Cell 2009, 136, 669–687. [Google Scholar] [CrossRef] [Green Version]
  81. Zhu, H.; Xia, R.; Zhao, B.; An, Y.Q.; Dardick, C.D.; Callahan, A.M.; Liu, Z. Unique expression, processing regulation, and regulatory network of peach (Prunus persica) miRNAs. BMC Plant Biol. 2012, 12, 149. [Google Scholar] [CrossRef] [Green Version]
  82. Tirumalai, V.; Swetha, C.; Nair, A.; Pandit, A.; Shivaprasad, P.V. miR828 and miR858 regulate VvMYB114 to promote anthocyanin and flavonol accumulation in grapes. J. Exp. Bot. 2019, 70, 4775–4792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Guan, X.; Pang, M.; Nah, G.; Shi, X.; Ye, W.; Stelly, D.M.; Chen, Z.J. miR828 and miR858 regulate homoeologous MYB2 gene functions in Arabidopsis trichome and cotton fibre development. Nat. Commun. 2014, 5, 3050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Zhao, T.; Tao, X.; Li, M.; Gao, M.; Chen, J.; Zhou, N.; Mei, G.; Fang, L.; Ding, L.; Zhou, B.; et al. Role of phasiRNAs from two distinct phasing frames of GhMYB2 loci in cis- gene regulation in the cotton genome. BMC Plant Biol. 2020, 20, 219. [Google Scholar] [CrossRef] [PubMed]
  85. Chen, X.M. A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science 2004, 303, 2022–2025. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Li, J.Y.; Millar, A.A. Expression of a microRNA-Resistant Target Transgene Misrepresents the Functional Significance of the Endogenous microRNA: Target Gene Relationship. Mol. Plant. 2013, 6, 577–580. [Google Scholar] [CrossRef] [Green Version]
  87. Baskar, V.; Gangadhar, B.H.; Park, S.W.; Nile, S.H. A simple and efficient Agrobacterium tumefaciens-mediated plant transformation of Brassica rapa ssp. pekinensis. 3 Biotech 2016, 6, 88. [Google Scholar] [CrossRef] [Green Version]
  88. Li, G.L.; Yue, L.X.; Li, F.; Zhang, S.F.; Zhang, H.; Qian, W.; Fang, Z.Y.; Wu, J.; Wang, X.W.; Zhang, S.J.; et al. Research progress on Agrobacterium tumefaciens-based transgenic technology in Brassica rapa. Hortic. Plant J. 2018, 4, 126–132. [Google Scholar] [CrossRef]
  89. Jian, W.; Cao, H.; Yuan, S.; Liu, Y.; Lu, J.; Lu, W.; Li, N.; Wang, J.; Zou, J.; Tang, N.; et al. SlMYB75, an MYB-type transcription factor, promotes anthocyanin accumulation and enhances volatile aroma production in tomato fruits. Hortic Res. 2019, 6, 22. [Google Scholar] [CrossRef] [Green Version]
  90. Zhang, Y.Y.; Wang, G.J.; Li, L.X.; Li, Y.H.; Zhou, B.; Yan, H.F. Identification and expression analysis of BrTT8 during anthocyanin biosynthesis and exposure to abiotic stress in turnip (Brassica rapa subsp. rapa ‘Tsuda’). Sci. Hortic. 2020, 268, 109332. [Google Scholar] [CrossRef]
  91. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Ijms 21 04326 g001 550
Figure 1. The phenotype of Brassica rapa subsp. rapa cv. Tsuda seedling and anthocyanin content after dark and sunlight treatment. (A) Seedlings were exposed to sunlight or darkness. (B) Anthocyanin was extracted from fresh whole seedlings, and the concentration was determined using absorbance at 530 nm per gram of fresh mass. Vertical bars indicate ± SE of replicate assays (n = 8). The statistically significant difference between light treatment and dark is shown (t-test, * p < 0.05).
Figure 1. The phenotype of Brassica rapa subsp. rapa cv. Tsuda seedling and anthocyanin content after dark and sunlight treatment. (A) Seedlings were exposed to sunlight or darkness. (B) Anthocyanin was extracted from fresh whole seedlings, and the concentration was determined using absorbance at 530 nm per gram of fresh mass. Vertical bars indicate ± SE of replicate assays (n = 8). The statistically significant difference between light treatment and dark is shown (t-test, * p < 0.05).
Ijms 21 04326 g001
Ijms 21 04326 g002 550
Figure 2. Anthocyanin distribution in the cross-section of the hypocotyl. Samples were prepared through the hand-sliced method and observed using a confocal microscope. Anthocyanin detection using 488 nm excitation and 545 nm emission. Samples of dark (A) and light (B) treatment viewed with anthocyanin detection field, bright field, and merged images. The scale bars in the figure represent 50 µm. The green color is indicative of anthocyanin presence.
Figure 2. Anthocyanin distribution in the cross-section of the hypocotyl. Samples were prepared through the hand-sliced method and observed using a confocal microscope. Anthocyanin detection using 488 nm excitation and 545 nm emission. Samples of dark (A) and light (B) treatment viewed with anthocyanin detection field, bright field, and merged images. The scale bars in the figure represent 50 µm. The green color is indicative of anthocyanin presence.
Ijms 21 04326 g002
Ijms 21 04326 g003 550
Figure 3. Three candidate primary miR828 transcripts were identified by qRT-PCR and sequences alignment showed they are respectively located in B. rapa chromosome A03, A08, and A01. (A) Electrophoresis of PCR product on agarose gel. (B) The sequences of these primary miR828 transcripts are highly similar except that primiR828-A01 has a 14 bp sequence deletion. The sequences in blue box and purple box respectively represent the region of miR828a and miR828b.
Figure 3. Three candidate primary miR828 transcripts were identified by qRT-PCR and sequences alignment showed they are respectively located in B. rapa chromosome A03, A08, and A01. (A) Electrophoresis of PCR product on agarose gel. (B) The sequences of these primary miR828 transcripts are highly similar except that primiR828-A01 has a 14 bp sequence deletion. The sequences in blue box and purple box respectively represent the region of miR828a and miR828b.
Ijms 21 04326 g003
Ijms 21 04326 g004 550
Figure 4. Phylogenetic analysis of candidate BrMYBs and multiple sequence alignment of R2R3 MYB family in Brassica rapa and AtMYBs in Arabidopsis. (A) The phylogenetic tree was constructed by using the Neighbor-Joining method [58] with 1000 bootstrap replicates [59]. Branches corresponding to partitions reproduced in less than 40% bootstrap replicates are collapsed. The evolutionary distances were computed using the Poisson correction method [60] and were conducted in MEGA7 [61] I to VII represents the MYB proteins divided subgroup I to subgroup VII, separately. (B) Protein sequences of subgroup VI were aligned by using Clustal X ver. 2 [62] with default settings, and the conserved amino acids were shaded by using GeneDoc (2.6) [63]. The R2, R3/bHLH-binding domains are marked with black bars and three α-helices of both R2 and R3 domains are indicated by square frame [64]. (C) The R2R3-domain sequences from subgroup VI were aligned and the sequence conservation at a particular position expressed as a stack of letters using the WebLogo program [65]. The overall height of each stack indicates the sequence conservation at that position (measured in bits), whereas the height of symbols within the stack reflects the relative frequency of the corresponding amino acid. The conserved tryptophan residues (Trp, W) in the MYB domain are marked with red asterisks. The replaced residues in the R3 repeat are shown by blue asterisks.
Figure 4. Phylogenetic analysis of candidate BrMYBs and multiple sequence alignment of R2R3 MYB family in Brassica rapa and AtMYBs in Arabidopsis. (A) The phylogenetic tree was constructed by using the Neighbor-Joining method [58] with 1000 bootstrap replicates [59]. Branches corresponding to partitions reproduced in less than 40% bootstrap replicates are collapsed. The evolutionary distances were computed using the Poisson correction method [60] and were conducted in MEGA7 [61] I to VII represents the MYB proteins divided subgroup I to subgroup VII, separately. (B) Protein sequences of subgroup VI were aligned by using Clustal X ver. 2 [62] with default settings, and the conserved amino acids were shaded by using GeneDoc (2.6) [63]. The R2, R3/bHLH-binding domains are marked with black bars and three α-helices of both R2 and R3 domains are indicated by square frame [64]. (C) The R2R3-domain sequences from subgroup VI were aligned and the sequence conservation at a particular position expressed as a stack of letters using the WebLogo program [65]. The overall height of each stack indicates the sequence conservation at that position (measured in bits), whereas the height of symbols within the stack reflects the relative frequency of the corresponding amino acid. The conserved tryptophan residues (Trp, W) in the MYB domain are marked with red asterisks. The replaced residues in the R3 repeat are shown by blue asterisks.
Ijms 21 04326 g004
Ijms 21 04326 g005 550
Figure 5. Relative expression of BrmiR828 and its candidate targets in seedlings of Brassica rapa with dark, blue, or UV-A light treatment was detected by qRT-PCR and calculated with 2−ΔΔCT. The ubiquitin (UBQ) gene served as the internal control to correct for template quantity. Each bar indicates the mean ± SE of three biological replicates assays after correcting for template quantity relative to the UBQ gene. Statistically significant differences between blue light treatment and dark, UV-A treatment and dark are shown (t-test, * p < 0.05; ns, not significant difference, df = 3)
Figure 5. Relative expression of BrmiR828 and its candidate targets in seedlings of Brassica rapa with dark, blue, or UV-A light treatment was detected by qRT-PCR and calculated with 2−ΔΔCT. The ubiquitin (UBQ) gene served as the internal control to correct for template quantity. Each bar indicates the mean ± SE of three biological replicates assays after correcting for template quantity relative to the UBQ gene. Statistically significant differences between blue light treatment and dark, UV-A treatment and dark are shown (t-test, * p < 0.05; ns, not significant difference, df = 3)
Ijms 21 04326 g005
Ijms 21 04326 g006 550
Figure 6. Sites of BrmiR828-mediated cleavage on target BrPAP1 (Bra039763), BrMYB82 (Bra022602) mRNAs and BrTAS4 and chromosome location and sequence alignment of TAS4 orthologs in Brassica rapa and Arabidopsis thaliana. (A) Positions of dominant 5′-RACE products of BrPAP1, BrMYB82, and BrTAS4 are indicated by a vertical arrow in the binding region. (B) Schematic diagram of the motif site of BrMYB targeted by BrmiR828, H1-3 are helices. (C) RLM 5′-RACE verified conserved binding site targeted by BrmiR828. (D) The sequences designated by dark line shows the binding site of miR828 and the complementary site of TAS4-siR81(−) [46]. The first siRNA, TAS4 D1(+) originated from the splicing site of BrmiR828 and TAS4.
Figure 6. Sites of BrmiR828-mediated cleavage on target BrPAP1 (Bra039763), BrMYB82 (Bra022602) mRNAs and BrTAS4 and chromosome location and sequence alignment of TAS4 orthologs in Brassica rapa and Arabidopsis thaliana. (A) Positions of dominant 5′-RACE products of BrPAP1, BrMYB82, and BrTAS4 are indicated by a vertical arrow in the binding region. (B) Schematic diagram of the motif site of BrMYB targeted by BrmiR828, H1-3 are helices. (C) RLM 5′-RACE verified conserved binding site targeted by BrmiR828. (D) The sequences designated by dark line shows the binding site of miR828 and the complementary site of TAS4-siR81(−) [46]. The first siRNA, TAS4 D1(+) originated from the splicing site of BrmiR828 and TAS4.
Ijms 21 04326 g006
Table
Table 1. The predict targets of miR828 and annotation by psRNATarget using Brassica rapa cds library (version 1.1, released on 30 August 2011).
Table 1. The predict targets of miR828 and annotation by psRNATarget using Brassica rapa cds library (version 1.1, released on 30 August 2011).
Target Gene IDAlignmentTarget Gene AnnotationExpect 1UPE 2
Bra001917miR828    22 ACCUUAUGGGUAAAUUCGUUCU 1
Target    310 UGGAACACCCAUCUGAGCAAGA 331		ATMYB90/PAP22.513.531
Bra004162miR828    22 ACCUUAUGGGUAAAUUCGUUCU 1
Target    310 UGGAACACCCAUUUGAGUAAGA 331		MYB domain protein MYB90/ PAP22.016.874
Bra004167miR828    22 ACCUUAUGGGUAAAUUCGUUCU 1
Target    1071 UCGAAAGCCCAUUGAAGUAAGA 1092		protein kinase family protein3.019.602
Bra022602miR828    22 ACCUUAUGGGUAAAUUCGUUCU 1
Target    322 UGGAAUACCCAUUUGAACAAGA 343		MYB domain protein MYB822.016.804
Bra029113miR828    22 ACCUUAUGGGUAAAUUCGUUCU 1
Target    322 UGGAAUACACAUUUGAACAAGA 343		MYB domain protein MYB823.013.444
Bra039763miR828    22 ACCUUAUGGGUAAAUUCGUUCU 1
Target    310 UGGAACACCCAUUUGAGUAAGA 331		PAP12.010.13
BraTAS4miRNA    22 ACCUUAUGGGUAAAUUCGUUCU 1
Target    663 UGGAAUAUUCAUUUGAGCAAGA 684		Non-coding RNA (Trans-acting siRNA gene 4)1.5N/A
1 Expect value, maximum cutoff of score based on given scoring schema. 2 UPE value, Target accessibility, the energy required to unpair target and small RNA.
Table
Table 2. The phased tas4 siRNAs identified by high throughput sequencing in UV-A treatment sRNA library of Brassica rapa.
Table 2. The phased tas4 siRNAs identified by high throughput sequencing in UV-A treatment sRNA library of Brassica rapa.
siRNA IDSequence (5′-3′)Length (bp)Reads NumbersiRNA Name [48]Phase Site in Sequence of BrTAS4 (From 5′ Terminal)
T0111568uugagcaagauguuggcaugaaau246TAS D1(+)675 uugagcaagauguuggcaugaaau 698
T0063765aauugcagcggugaaggacgagcu2411TAS D2(+)696 aauugcagcggugaaggacgagcu 719
T0038241auaacaacgugaaggaucgagguc2418TAS D4(-)742 gaccucgauccuucacguuguuau 765
T0229627acgugaaggaucgaggucgaggua243TAS D4(-)736 uaccucgaccucgauccuucacgu 759
T0197894aauaacaacgugaaggaucgaggu244TAS D4(-)743 accucgauccuucacguuguuauu 766
T0099445uuaauuuguaugcaugcugauaac247TAS D10(+)885 uuaauuuguaugcaugcugauaac 908
T0159302uuuguaugcaugcugauaacugaa244TAS D10(+)889 uuuguaugcaugcugauaacugaa 912

Share and Cite

MDPI and ACS Style

Zhou, B.; Leng, J.; Ma, Y.; Fan, P.; Li, Y.; Yan, H.; Xu, Q. BrmiR828 Targets BrPAP1, BrMYB82, and BrTAS4 Involved in the Light Induced Anthocyanin Biosynthetic Pathway in Brassica rapa. Int. J. Mol. Sci. 2020, 21, 4326. https://doi.org/10.3390/ijms21124326

AMA Style

Zhou B, Leng J, Ma Y, Fan P, Li Y, Yan H, Xu Q. BrmiR828 Targets BrPAP1, BrMYB82, and BrTAS4 Involved in the Light Induced Anthocyanin Biosynthetic Pathway in Brassica rapa. International Journal of Molecular Sciences. 2020; 21(12):4326. https://doi.org/10.3390/ijms21124326

Chicago/Turabian Style

Zhou, Bo, Jingtong Leng, Yanyun Ma, Pengzhen Fan, Yuhua Li, Haifang Yan, and Qijiang Xu. 2020. "BrmiR828 Targets BrPAP1, BrMYB82, and BrTAS4 Involved in the Light Induced Anthocyanin Biosynthetic Pathway in Brassica rapa" International Journal of Molecular Sciences 21, no. 12: 4326. https://doi.org/10.3390/ijms21124326

APA Style

Zhou, B., Leng, J., Ma, Y., Fan, P., Li, Y., Yan, H., & Xu, Q. (2020). BrmiR828 Targets BrPAP1, BrMYB82, and BrTAS4 Involved in the Light Induced Anthocyanin Biosynthetic Pathway in Brassica rapa. International Journal of Molecular Sciences, 21(12), 4326. https://doi.org/10.3390/ijms21124326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop