Next Article in Journal
Genome-Wide Analysis of the FABP Gene Family in Liver of Chicken (Gallus gallus): Identification, Dynamic Expression Profile, and Regulatory Mechanism
Previous Article in Journal
Clinical Implications of Extracellular HMGA1 in Breast Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Flower-Specific Promoters through Comparative Transcriptome Analysis in Brassica napus

Key Laboratory of Biology and Genetics Improvement of Oil Crops, Oil Crops Research Institute of Chinese Academy of Agricultural Sciences, Ministry of Agriculture and Rural Affairs, Wuhan 430062, Hubei, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(23), 5949; https://doi.org/10.3390/ijms20235949
Submission received: 31 October 2019 / Revised: 20 November 2019 / Accepted: 25 November 2019 / Published: 26 November 2019
(This article belongs to the Section Molecular Plant Sciences)

Abstract

:
Brassica napus (oilseed rape) is an economically important oil crop worldwide. Sclerotinia stem rot (SSR) caused by Sclerotinia sclerotiorum is a threat to oilseed rape production. Because the flower petals play pivotal roles in the SSR disease cycle, it is useful to express the resistance-related genes specifically in flowers to hinder further infection with S. sclerotiorum. To screen flower-specific promoters, we first analyzed the transcriptome data from 12 different tissues of the B. napus line ZS11. In total, 249 flower-specific candidate genes with high expression in petals were identified, and the expression patterns of 30 candidate genes were verified by quantitative real-time transcription-PCR (qRT-PCR) analysis. Furthermore, two novel flower-specific promoters (FSP046 and FSP061 promoter) were identified, and the tissue specificity and continuous expression in petals were determined in transgenic Arabidopsis thaliana with fusing the promoters to β-glucuronidase (GUS)-reporter gene. GUS staining, transcript expression pattern, and GUS activity analysis indicated that FSP046 and FSP061 promoter were strictly flower-specific promoters, and FSP046 promoter had a stronger activity. The two promoters were further confirmed to be able to direct GUS expression in B. napus flowers using transient expression system. The transcriptome data and the flower-specific promoters screened in the present study will benefit fundamental research for improving the agronomic traits as well as disease and pest control in a tissue-specific manner.

1. Introduction

Brassica napus L. is an evolutionarily young allotetraploid species formed about 7500 years ago by natural hybridization between Brassica oleracea and Brassica rapa, followed by chromosome doubling [1,2]. B. napus is an economically important oil crop all over the world. It not only provides vegetable oil for people, but also provides high quality fodder for animals. In addition, owing to its favorable agronomic properties, such as cultivation under different seasons and rotation with cereals, B. napus is preferred by farmers worldwide.
However, Sclerotinia stem rot (SSR) is destructive to oilseed rape production. SSR is caused by Sclerotinia sclerotiorum (Lib.) de Bary, a cosmopolitan pathogen of many economically important crops. As a necrotrophic pathogen, it infects more than 600 plant species, including important oil crops such as oilseed rape, soybean, and sunflower [3,4,5,6]. SSR not only deteriorates the quality of the seed, but also significantly reduces the oil content [7].
Numerous studies have shown that flower petals play vital roles across the whole infection cycle of SSR, by supplying nutrition for ascospore germination and hyphae development [8,9,10]. As senescent petals fall onto leaves, petioles, and stems, secondary infection of SSR disease is initiated [8,9]. Therefore, the disease incidence is significantly positively correlated with petal infestation [9]. When the conditions for infection are favorable, S. sclerotirum causes severe yield losses in oilseed rape in the field [10].
Genetic engineering of crops, where resistant genes are introduced to improve crops resistance against disease, is a rapid breeding approach [11]. The cauliflower mosaic virus 35S promoter, actin, and maize ubiquitin promoters are the most commonly used constitutive promoters. However, these constitutive promoters trigger gene overexpression in all tissues, leading to some negative or undesirable pleiotropic effects, which interfere with normal gene expression and agronomic performance [12,13,14]. Tissue-specific promoters driving transgene expression within a given tissue are thus preferred. With control by these promoters, transgenes could result in spatio-temporal expression, which could eliminate unnecessary energy waste during normal plant growth and would make it easier to achieve precise and reliable results [15,16,17,18]. Consequently, specifically expressing resistant genes in oilseed rape flowers to control S. sclerotirum infection and extension will likely be effective.
To date, several flower-specific promoters have been reported and applied to research in the cut-flower industry. Most of these promoters are flavonoid synthesis-related gene promoters, and are used for altering corolla colors to create novel varieties in ornamental flowers, such as lily, petunia, and rose [19,20,21,22]. In Brassicacea, Arabidopsis thaliana PISTILLATA (PI) and APETALA3 (AP3) are B-class organ identity genes in the ABC model, which are required for petal and stamen development [23,24,25,26]. AtPI promoter driven GUS expression only occurs in the petals and stamens in A. thaliana [27]. However, the promoter of B. napus homolog gene PISTILLATA-1 showed low flower specificity, under which GUS could also express in leaf and silique [28]. BnQRT3 promoter drives GUS expressing in branch connective tissue between the pedicels and stem, besides in flowers [29]. The expression patterns of homologous genes between A. thaliana and B. napus are different, even though they all belong to Brassicacea. Consequently, it is necessary to identify new flower-specific promoters in B. napus per se. In this study, flower-specific promoters in B. napus were identified based on transcriptome data from 12 different tissues including root, stem, leaf, flower bud, and six dissected flower parts. Two novel strictly flower-specific promoters were cloned and experimentally confirmed by a GUS reporter gene expression system. The transcriptome data and the flower-specific promoters screened in the present study will benefit fundamental research for improving the agronomic traits as well as disease and pest control in a tissue-specific manner.

2. Results

2.1. Transcriptome Analysis and Identification of Tissue-Specific Genes

Twelve different tissues from B. napus line ZS11 were collected and used to construct libraries for RNA sequencing. After quality filtering, we obtained 18 million to 51 million clean reads for each tissue library (SRA accession: PRJNA474576). Through mapping to the B. napus genome database, 67.85 to 87.88% of the clean reads could be mapped. Over 12 million to 26 million unique reads were obtained, which occupied 51.02 to 75.54% of the mapped data (Table 1).
The number of the tissue-specific genes were determined and summarized in Table S1. There are the most specific genes: 1074 genes in Flower bud, while the specific genes in the stem were the least: 38 genes. After filtering with high stringency based on the FPKM (fragments per kilobase of transcript per million mapped reads) value (Figure 1a), 1317 genes were detected as specifically expressed in floral tissues (including bud, pistil, stamen, sepal, lossomy petal (BP), and wilting petal (WP)). We further narrowed this to a list of 249 candidate genes that exhibited high expression in two kinds of petals (Figure 1a and Table S2). The 249 genes were regarded as flower-specific candidate genes. After Gene Ontology (GO) classification and enrichment using PlantTFDB 4.0 online tool, 143 genes had GO annotations (Figure 1b). Many genes were involved in the metabolic process of biological processes, and in the binding of molecular function. These results suggested that flower-specific candidate genes play different roles in flowers. Thirty-one GO terms were enriched, many of which were related with flower organ formation, development or specification. As shown in Figure 1c, five genes were related with flower development and reproductive shoot system development; four genes participated in pattern specification process, floral organ development and organ morphogenesis. These results indicated that the flower-specific candidates we screened were reliable.

2.2. qRT-PCR Analysis of Flower-Specific Candidate Genes

Thirty of 249 flower-specific candidates, most of which specifically expressed in BP or WP, were selected for quantitative real-time transcription-PCR (qRT-PCR) analysis. In the first round of detection, we defined flower-specific genes as those with expression values in the flower buds that were at least three times higher than that in other tissues (Figure S1). Six genes (FSP046—BnaA03g02790D, FSP061—BnaA03g18510D, FSP102—BnaC07g09080D, FSP089—BnaA07g11640D, FSP187—BnaUnng01380D, FSP098—BnaC01g10920D) exhibited specific and high expression in flower buds (Figure 2a). Then, the expression patterns of these six genes were detected in dissected flower parts (Figure 2b). To screen the continuous expression during the flowering stage, we detected the expression levels of the six genes in two kinds of petals (BP and WP) in detail. As illustrated in Figure 2b, five of the six candidates (excepting FSP102) were highly expressed in BP or WP, while only FSP046 and FSP061 were highly expressed in BP and WP at the same time. In addition, the expression levels of the four candidates (FSP046, FSP061, FSP187, and FSP098) in the stamen were relatively high.

2.3. FSP046 and FSP061 Promoters are Flower-Specific in Transgenic A. thaliana

To detect the flower-specificity, about 3.5 Kb of upstream sequences of two flower-specific candidate genes (FSP046 and FSP061) from the ZS11 cultivar were cloned and sequenced. The cis elements in these two promoter regions were analyzed using the PlantCARE online tool; and the results are illustrated in Figure 3. One P-box, which contributes flower-specific gene expression, is found in the FSP061 promoter sequence, but none are identified in the FSP046 promoter sequence. There are several MYB, and MYC elements located in both promoter regions. Many cis elements related to light responsiveness, such as Box 4, I-box, TCT-motif, and GT1-motif, are found in both FSP046 and FSP061 promoters.
To detect the tissue-specific expression, the whole upstream sequences of FSP046 and FSP061 gene were fused to the GUS reporter gene (FSP046::GUS and FSP061::GUS) and transformed into A. thaliana. Histochemical staining was detected during the two-leaf, six-leaf, and flowering stages of T2 generation plants by histochemical assay (Figure 4 and Figure S2). The GUS protein within FSP046::GUS transgenic plants was specifically expressed in flowers (Figure 4e–h). Moreover, there is the highest GUS activity in stamen (Figure 4f) and less activity in the petals (Figure 4f–h), which is in good accordance with FSP046 gene expression pattern in B. napus flowers (Figure 2b). However, no GUS staining was detected in the roots (Figure 4a,b), leaves (Figure 4a–c), siliques (Figure 4i), or seeds (Figure 4j). The FSP061::GUS transgenic plants showed also flower-specific and similar GUS expression pattern (Figure S2).
To accurately define the tissue-specific expression pattern, we performed qRT-PCR analysis and GUS activity assays in transgenic A. thaliana (Figure 5). The GUS gene expression was exceedingly high in flower comparing with that in root, stem, and leaf. Moreover, the expression levels of the GUS gene in the flower were about seven and 25 fold higher than that in silique under FSP046 and FSP061 promoter, respectively (Figure 5a). The GUS activity results also showed the same expression pattern among different tissues (Figure 5b). Furthermore, we also detected GUS activity in transgenic A. thaliana seeds. Of which the activity was obviously lower than that in silique. In summary, FSP046 and FSP061 promoter could drive GUS specifically expressed in flowers, and the activity of FSP046 promoter was higher than that of FSP061 promoter.

2.4. FSP046 and FSP061 promoter can Drive GUS Expressing in Flowers of B. napus

Based on the flower-specific expression in transgenic A. thaliana, we further investigated whether the GUS gene could be expressed under FSP046 and FSP061 promoter driving in B. napus flowers. Using the transient expression system, we detected that both of the flowers can be stained by GUS reaction buffer. As shown in Figure 6, comparing with inoculating empty Agrobacterium, the sepals, petals, stamen, and stigmas of FSP046::GUS and FSP061::GUS could be clearly stained, which was similarly found for the 35S::GUS positive control transiently transformed flowers. These results confirmed that both promoters can drive GUS expressing in flowers of B. napus.

3. Discussion

In genetic engineering, using tissue-specific or temporal-specific promoters to drive gene expression can avoid unwanted influences on plant phenotypes caused by constitutive promoters [17,18,30]. To find tissue-specific promoters, researchers usually identify the genes expressed in a tissue-specific manner. Previously, gene expression was typically analyzed by subtractive hybridization [31] and gene expression microarray [32], which were time consuming or had high background noise. In contrast, RNA-seq is a low-cost and high-throughput sequencing technology for genome-scale gene expression analysis [33]. Analyzing transcriptome data from different tissues to obtain tissue-specific genes is a high throughput, practical, and rapid method, which has been applied in peanut and poplar [34,35]. In the present study, we sequenced the transcriptomes of 12 different tissues from B. napus. After analyzing the data, 249 flower-specific genes highly expressed in petals were identified, and two novel flower-specific promoters were selected for cloning and confirmed to be flower specific through driving GUS reporter gene expressing in transgenic A. thaliana plants.
In present study, both FSP046 and FSP061 promoters were proved to be flower-specific. It was reported that the P-box contributes to flower-specific gene expression of gPAL2 of Phaseolus vulgaris [36], chsA of Petunia hybrid [37] and chs of Phalaenopsis hybrid [38]. However, only the FSP061 promoter sequence contained the P-box (sequence AACCAAAC) according to the P-box consensus sequence [37]. Therefore, we presumed that there may be other unknown motifs contributing to flower-specific expression in B. napus. In future studies, GUS activity of different deletions of the two promoters will be implemented to find novel flower-specific motifs, and this method will also be used to identify the core regions of FSP046 and FSP061 promoters.
From Figure 4 and Figure S2, the stamen was deeply stained by GUS reaction buffer, and exhibit higher GUS activity than petal. AtPI and AtAP3 promoters also drive GUS expressing in petal and stamen [26,27]. Oilseed rape XY355 was reported as a petal-specific promoter [39], while XY355::OvPAP2 transgenic B. napus plants exhibited red petals and stamen [40]. It was worth noticing that genes that were expressed in petals were usually accompanied by expression in stamen. Up to now, there is only one promoter (InMYB1 promoter) driving GUS protein expression in petals distinctly and uniquely in A. thaliana [41]. Young and Werner reported that stamens were probably an important infection route of S. sclerotiorum for apetalous winter oilseed rape [42]. Consequently, FSP046 and FSP061 promoters are suitable to apply in SSR disease control in a flower-specific manner.
In transgenic A. thaliana, both of the GUS gene expression level and protein activity in FSP046::GUS transgenic plants were higher than that of FSP061::GUS in flowers (Figure 5). These results were consistent with the qRT-PCR results in B. napus (Figure 2b), which indicated that the FSP046 promoters have stronger activity to drive downstream gene expression. In addition, from Figure 5a, the GUS gene could also be detected in silique with a relative higher expression level. To further detect whether GUS was expressed in transgenic seeds, we directly assayed the GUS activity, and found that the activity were decreased by about 98% and 54% in FSP046::GUS and FSP061::GUS transgenic seeds, respectively (Figure 5b). These results indicated that FSP046 and FSP061 promoters have a strict flower specificity, which were more safety to environment for application as the flower is senescent tissue.
In the present study, we use the syringe-press transient expression system to over-express FSP046::GUS and FSP061::GUS in B. napus flowers. Using the transient expression system, all flowers (especially sepals and petals) could be stained by GUS reaction buffer (Figure 6), which indicated that both promoters can drive GUS expression in B. napus flowers. Compared with hypocotyl transformation, the transient expression system in flowers was fast and efficient, which could also be applied in subcellular localization, BiFC assay, and gene function analysis of B. napus genes.
The multiple functions of oilseed rape (including seed oil, vegetable bolt, flower sightseeing, and fodder for animals) are heavily researched. Many regions in China have constructed large-scale oilseed flower sightseeing tourism, which has become an important part of the local economy [43]. To add aesthetic value, genetic engineering under the direction of flower specific promoters is a timesaving way to produce rich flower colors of oilseed rape. Furthermore, the flower-specific promoters isolated in this study could also be used to modify colors for ornamental flowers such as lily and chrysanthemum. SSR is the major threat for oilseed rape production, and it was reported that pollen beetle (Meligethes aeneus) was a major pest at the inflorescence stage [44]. In recent years, cross-kingdom RNAi has been discovered between host plant and pathogens, pests [45,46,47]. Host-induced gene silencing (HIGS) by transgenic expression of pathogen/pest gene-targeting double-stranded (ds)RNA is a promising alternative to control disease/pest in plant protection [48]. Our flower-specific promoter has the potential to be combined with the HIGS technique to control SSR and pests in a tissue-specific manner.

4. Materials and Methods

4.1. Plant Material Collection

Seeds of B. napus line ZS11 were obtained from Key Laboratory of Biology and Genetic Improvement of Oil Crops at OCRI and sown in October in the field located in Wuhan (average annual temperature 13–22 °C and 70–80% humidity), Hubei province, China. Hubei province is the major production area for winter rapeseed planting in the middle and lower Yangtze River. The trial management followed standard breeding field protocols. Tissue samples were collected at the full-bloom stage in the following spring. Root, stem, leaf, flower bud, little siliques (about 4 cm long), and dissected flower parts including blossomy petal (petal at the blooming day, BP), wilting petal (one day after blooming, WP), pistil, stamen, sepal, ovule, and pericarp were collected from 15 plants in the same developmental stage in the same test plot.

4.2. RNA Preparation and RNA-seq Analysis

Total RNA was extracted from different tissues using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). Poly-A-containing mRNA was isolated from the total RNA using Oligotex mRNA Mini Kit (QIAGEN, San Francisco/Bay area, CA, USA) according to the manufacturers’ instructions. RNA libraries with insert sizes of 250 bp were constructed for each tissue sample and sequenced on an Illumina HiSeq 2000 platform (San Diego, CA, USA) at BGI Co. Ltd. The libraries were sequenced for paired-end reads of 150 bp.

4.3. Functional Annotation

The raw reads were filtered to remove adapter sequences, low-quality reads, and reads containing poly-N using the tools from the NGS QC Toolkit (v2.3) [49]. The derived clean reads were mapped to the reference genome of B. napus line Darmor-bzh (Brassica napus Annotation v5, http://www.genoscope.cns.fr/brassicanapus/) [1] using TopHat software (v2.0.9) [50]. The mapped reads were quantified by fragments per kilobase of transcript per million mapped reads (FPKM) for each gene using Cufflinks (v2.1.0) software [50].

4.4. Screening of Tissue-Specific Candidate Genes

Tissue-specific genes were screened and calculated using the threshold of FPKM value >2. We defined flower-specific genes as genes with threshold FPKM values ≥15 in either flower tissues (including bud, pistil, stamen, sepal, BP, and WP), while the FPKM values were <1 in all of the other tissues (root, stem, leaf, ovule, and pericarp). Furthermore, according to the trait of S. sclerotiorum tending to infect the plant petals primarily in the SSR disease cycle, the genes with FPKM values ≥15 in BP or WP were finally selected as flower-specific candidates. GO functional annotation and classification of the flower-specific candidate genes were conducted in the PlantTFDB 4.0 online tool (http://planttfdb.cbi.pku.edu.cn/) with the threshold of p-value ≤ 0.01 [51,52]. The GO enrichment result was displayed using the ImageGP online tool (http://www.ehbio.com/ImageGP/).

4.5. Quantitative Real-Time PCR Analysis

To verify the tissue specificity of flower-specific candidate genes, qRT-PCR primers of flower-specific candidate genes were designed (Table S3). Primer specificity was examined by PCR and agarose gel analysis. Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). About 1 μg of total RNA was reverse transcribed using the PrimeScript™ RT reagent Kit with gDNA Eraser (TaKaRa Co., LTD, Beijing, China). To screen the genes continuously expressed during the flowering stage (from flower bud to WP), the expression patterns of flower-specific candidate genes were first detected in five tissues including root, stem, leaf, flower bud, and silique by qRT-PCR analysis. The genes with specific expression in flower bud were used for subsequent qRT-PCR in dissected flower parts including sepal, pistil, stamen, BP, and WP. Three biological repeats were conducted. The B. napus β-actin gene (AF111812) was used as a reference standard. The relative expression was calculated using the 2−∆∆Ct method [53].

4.6. Cloning and Characterization of Flower-Specific Promoters

About 3.5 Kb sequences upstream of the start codon ATG were selected, and the primers (Table S3) were designed according to B. napus line ZS11 genome [54]. To analyze the significant cis elements of the flower-specific candidate genes, accurate promoter sequences were analyzed using the PlantCARE online tool (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [55].

4.7. Binary Vector Construction and Plant Transformation

To construct the promoter::GUS reporter system, flower-specific gene promoters were cloned into pBI121 at HindIII and BamHI sites to replace the CaMV 35S promoter using the ClonExpress II One Step Cloning Kit (Vazyme Biotech Co., LTD, Nanjing, China). The primer pairs are listed in Table S3. Recombinant plasmids were introduced into Agraobecterium tumefaciens GV3101 competent cell (AngYu Biotech Co., LTD, Shanghai, China). A. thaliana transformation was performed using the floral dipping method [56]. For transient expression in flowers of B. napus, a syringe-press method was used for Agro-infiltration [57]. Briefly, the Agrobacterium cells were harvested and adjusted to OD600 2.0. The resuspension solutions were put into a 50 mL syringe, and then, 10–15 flowers from B. napus line Westar were added into the solution. Positive pressure was used for 30 s, and the inoculated flowers were placed into wet box for preserving moisture. Forty-eight hours later, the flower tissues were collected for GUS staining. Three independent biological repeats were performed.

4.8. Histochemical GUS Assays and Evaluation of GUS Activity

Histochemical detection of GUS staining was performed as described by Jefferson et al. [44]. Transgenic positive A. thaliana seedlings, organs, and tissues from different T2 generation lines and flowers of B. napus were infiltrated in GUS reaction solution (Coolaber Science & Technology Co., LTD, Beijing, China) for 12 h at 37 °C in the dark. Several washes of 70% ethanol were performed to stop the GUS reaction and remove chlorophyll at room temperature. Then the stained tissues were observed under an anatomical lens (Olympus SZX16, Ac Adapter, Tokyo, Japan).
GUS activity from different transgenic positive A. thaliana tissues was measured according to Jefferson et al. [58] by measuring 4-methylumbelliferone (4-MU) fluorometric quantity. Different tissues were grounded to a finely pulverized powder in liquid nitrogen, and extracted in GUS extraction buffer. The protein concentration was measured with the Bradford method [59]. Fluorescence values were recorded with Glomax® Explorer Multimode Microplate Reader (Promega, Madison, WI, USA). The GUS enzyme activity was calculated as picomoles of 4-MU produced per milligrams of protein per minute. The GUS gene expression in different tissues was examined by qRT-PCR analysis. The A. thaliana actin gene (AT3G18780) was used as a reference gene. The used primers were listed in Table S3.

5. Conclusions

In the present study, we sequenced the transcriptomes of 12 different tissues from B. napus, which provide a foundation for the expression patterns of gene family studies. Moreover, based on the RNA sequencing data, other tissue-specific promoters (such as in the root and stem) could also be screened using similar procedures, which could be used for disease and pest management in a tissue-specific manner, for example with clubroot disease and aphids. Therefore, the transcriptome data present in this study will provide practical guidance for other studies. In conclusion, the transcriptome data and the flower-specific promoters screened in the present study will benefit fundamental research and disease and pest control in a tissue-specific manner.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/20/23/5949/s1. Figure S1. Expression profiles of 30 flower-specific candidate genes from five different tissues by qRT-PCR. Tissues from five plants were collected together, three technical repeats were performed and three biological repeats were conducted. The B. napus β-actin gene (AF111812) was used as a reference standard. R (root), S (stem), L (leaf), F (flower bud), SQ (silique). The asterisk indicates that the expression value in the flower bud was at least three times higher than that in other tissues. Figure S2. Histochemical staining of FSP061::GUS transgenic A. thaliana plants. There is no GUS staining in the two-leaf stage (a), six-leaf stage (b), rosette leaf (c), silique (i), or seed (j). GUS activity driven by FSP061 promoter was observed in flower bud (d), flowers one day before blooming (e), flowers at the blooming day (f), petals one day after blooming (g), and petals two days after blooming (h). Moreover, there is the highest GUS activity in the stamen/stigma (f) and less activity in the petals (f, g, and h). Bar = 1 mm. Table S1 Tissue-specific gene data of 12 different tissues from Brassica napus. Table S2 Transcriptome data of 249 flower-specific candidate genes from 12 different tissues. Table S3 Primers used in this study.

Author Contributions

Data curation, M.H., C.T., and M.C.; formal analysis, Y.L. (Yan Li), C.D., C.T., and R.Z.; funding acquisition, S.L.; methodology, C.D. and Z.B.; resources, Y.L. (Yueying Liu) and X.C.; software, M.H.; supervision, J.H. and S.L; writing―original draft, Y.L. (Yan Li); writing―review and editing, Y.L. (Yan Li), Z.B., J.H., and S.L.

Funding

This research was funded by National Key Research and Development Program of China (2016YFD0101001, 2018YFD0200904, and 2018YFE0108000), National Natural Science Foundation of China (NO. 31901871, 31770250), Postdoctoral Science Foundation of China (2016M601190), Earmarked Fund for China Agriculture Research System (CARS-12), and Agricultural Science and Technology Innovation Program (ASTIP) of Chinese Academy of Agricultural Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

SSRSclerotinia stem rot
BPBlossomy petal
WPWilting petal
FPKMFragments per kilobase of transcript per million mapped reads
GOGene Ontology
Qrt-PCRQuantitative real-time PCR
4-MU4-methylumbelliferone

References

  1. Chalhoub, B.; Denoeud, F.; Liu, S.Y.; Parkin, I.A.P.; Tang, H.B.; Wang, X.Y.; Chiquet, J.; Belcram, H.; Tong, C.B.; Samans, B.; et al. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 2014, 345, 950–953. [Google Scholar] [CrossRef]
  2. Wang, J.; Lydiate, D.J.; Parkin, I.A.P.; Falentin, C.; Delourme, R.; Carion, P.W.C.; King, G.J. Integration of linkage maps for the Amphidiploid Brassica napus and comparative mapping with Arabidopsis and Brassica rapa. BMC Genom. 2011, 12. [Google Scholar] [CrossRef] [PubMed]
  3. Boland, G.J.; Hall, R. Index of plant hosts of Sclerotinia sclerotiorum. Can. J. Plant. Pathol. 1994, 16, 93–108. [Google Scholar] [CrossRef]
  4. Bolton, M.D.; Thomma, B.P.H.J.; Nelson, B.D. Sclerotinia sclerotiorum (Lib.) de Bary: Biology and molecular traits of a cosmopolitan pathogen. Mol. Plant Pathol. 2006, 7, 1–16. [Google Scholar] [CrossRef] [PubMed]
  5. Kamal, M.M.; Savocchia, S.; Lindbeck, K.D.; Ash, G.J. Biology and biocontrol of Sclerotinia sclerotiorum (Lib.) de Bary in oilseed Brassicas. Australa. Plant Path. 2016, 45, 1–14. [Google Scholar] [CrossRef]
  6. Liang, X.; Rollins, J.A. Mechanisms of broad host range necrotrophic pathogenesis in Sclerotinia sclerotiorum. Phytopathology 2018, 108, 1128–1140. [Google Scholar] [CrossRef]
  7. Aggarwal, R.A.K.; Kumar, A.; Thakur, H.L. Effect of Sclerotinia rot on oil quality in low erucic acid cultivars of rapeseed. Cruciferae Newsl. 1997, 19, 103–104. [Google Scholar]
  8. Turkington, T.K.; Raa, M. Use of petal infestation to forecast Sclerotinia stem rot of canola: The influence of inoculum variation over the flowering period and canopy density. Phytopathology 1993, 83, 561–565. [Google Scholar] [CrossRef]
  9. Inglis, G.D.; Boland, G.J. The microflora of bean and rapeseed petals and the influence of the microflora of bean petals on white mold. Can. J. Plant Pathol. 1990, 2, 129–134. [Google Scholar] [CrossRef]
  10. Saharan, G.S.; Mehta, N. Sclerotinia Diseases of Crop Plants: Biology, Ecology and Disease Management, 1st ed.; Springer: Dordrecht, The Netherlands, 2008; pp. 71–75. [Google Scholar]
  11. Van der Biezen, E.A. Quest for antimicrobial genes to engineer disease-resistant crops. Trends Plant Sci. 2001, 6, 89–91. [Google Scholar] [CrossRef]
  12. Shikata, M.; Narumi, T.; Yamaguchi, H.; Sasaki, K.; Aida, R.; Oshima, Y.; Takiguchi, Y.; Ohme-Takagi, M.; Mitsuda, N.; Ohtsubo, N. Efficient production of novel floral traits in torenia by collective transformation with chimeric repressors of Arabidopsis transcription factors. Plant Biotechnol.-Nar 2011, 28, 189–199. [Google Scholar] [CrossRef]
  13. Pino, M.T.; Skinner, J.S.; Park, E.J.; Jeknić, Z.; Hayes, P.M.; Thomashow, M.F.; Chen, T.H. Use of a stress inducible promoter to drive ectopic AtCBF expression improves potato freezing tolerance while minimizing negative effects on tuber yield. Plant Biotechnol. J. 2010, 5, 591–604. [Google Scholar] [CrossRef] [PubMed]
  14. Yoshimi, O.; Masahito, S.; Tomotsugu, K.; Norihiro, O.; Nobutaka, M.; Masaru, O.T. MIXTA-like transcription factors and WAX INDUCER1/SHINE1 coordinately regulate cuticle development in Arabidopsis and Torenia fournieri. Plant Cell 2013, 25, 1609–1624. [Google Scholar]
  15. Hsieh, T.H.; Lee, J.T.; Charng, Y.Y.; Chan, M.T. Tomato plants ectopically expressing Arabidopsis CBF1 show enhanced resistance to water deficit stress. Plant Physiol. 2002, 130, 618–626. [Google Scholar] [CrossRef] [PubMed]
  16. Zavallo, D.; Bilbao, M.L.; Hopp, H.E.; Heinz, R. Isolation and functional characterization of two novel seed-specific promoters from sunflower (Helianthus annuus L.). Plant Cell Rep. 2010, 29, 239–248. [Google Scholar] [CrossRef]
  17. Hou, J.; Jiang, P.; Qi, S.; Ke, Z.; He, Q.; Xu, C.; Ding, Z.; Zhang, K.; Li, K. Isolation and Functional Validation of Salinity and Osmotic Stress Inducible Promoter from the Maize Type-II H+-Pyrophosphatase Gene by Deletion Analysis in Transgenic Tobacco Plants. Plos ONE 2016, 11, e0154041. [Google Scholar] [CrossRef]
  18. Zhang, H.; Jing, R.; Mao, X. Functional Characterization of TaSnRK2.8 Promoter in Response to Abiotic Stresses by Deletion Analysis in Transgenic Arabidopsis. Front. Plant Sci. 2017, 8, 1198. [Google Scholar] [CrossRef]
  19. Liu, Y.L.; Lou, Q.; Xu, W.R.; Xin, Y.; Bassett, C.; Wang, Y.J. Characterization of a chalcone synthase (CHS) flower-specific promoter from Lilium orential ‘Sorbonne’. Plant Cell Rep. 2011, 30, 2187–2194. [Google Scholar] [CrossRef]
  20. Sasaki, K.; Yamaguchi, H.; Narumi, T.; Shikata, M.; Oshima, Y.; Nakata, M.; Mitsuda, N.; Ohme-Takagi, M.; Ohtsubo, N. Utilization of a floral organ-expressing AP1 promoter for generation of new floral traits in Torenia fournieri Lind. Plant. Biotech.-Nar. 2011, 28, 181–188. [Google Scholar] [CrossRef]
  21. An, N.U.; Hu, Y.; Li, P.; Yuchi, Z.; Chen, Y.; Zhang, Y. Cloning and expression of a non-ribosomal peptide synthetase to generate blue rose. ACS Synth Biol. 2019, 8, 1698–1704. [Google Scholar] [CrossRef]
  22. Ma, N.; Chen, W.; Fan, T.; Tian, Y.; Zhang, S.; Zeng, D.; Li, Y. Low temperature-induced DNA hypermethylation attenuates expression of RhAG, an AGAMOUS homolog, and increases petal number in rose (Rosa hybrida). BMC Plant Biol. 2015, 15, 1–13. [Google Scholar] [CrossRef] [PubMed]
  23. Bowman, J.L.; Smyth, D.R.; Meyerowitz, E.M. Genes directing flower development in Arabidopsis. Plant Cell 1989, 1, 37–52. [Google Scholar] [PubMed]
  24. Hill, J.P.; Lord, E.M. Floral development in Arabidopsis thaliana: Comparison of the wildtype and the homeotic pistillata mutant. Can. J. Bot. 1989, 67, 2922–2936. [Google Scholar] [CrossRef]
  25. Jack, T.; Brockman, L.L.; Meyerowitz, E.M. The homeotic gene APETALA3 of Arabidopsis thaliana encodes a MADS box and is expressed in petals and stamens. Cell 1992, 68, 683–697. [Google Scholar] [CrossRef]
  26. Hill, T.A.; Day, C.D.; Zondlo, S.C.; Thackeray, A.G.; Irish, V.F. Discrete spatial and temporal cis-acting elements regulate transcription of the Arabidopsis floral homeotic gene APETALA3. Development 1998, 125, 1711–1721. [Google Scholar]
  27. Honma, T.; Goto, K. The Arabidopsis floral homeotic gene PISTILLATA is regulated by discrete cis-elements responsive to induction and maintenance signals. Development 2000, 127, 2021–2030. [Google Scholar]
  28. Roh, K.H.; Choi, S.B.; Kang, H.C.; Kim, J.B.; Kim, H.U.; Lee, K.R.; Sun, H.K. Isolation and functional characterization of a PISTILLATA-1 gene promoter from Brassica napus. J. Korean Soc. Appl. Biol. 2014, 57, 759–768. [Google Scholar] [CrossRef]
  29. Wan, L.; Xia, X.; Hong, D.; Yang, G. Molecular analysis and expression of a floral organ-specific polygalacturonase gene isolated from rapeseed (Brassica napus L.). Mol. Biol. Rep. 2010, 37, 3851–3862. [Google Scholar] [CrossRef]
  30. Chiou, C.Y.; Wu, K.Q.; Yeh, K.W. Characterization and promoter activity of chromoplast specific carotenoid associated gene (CHRC) from Oncidium Gower Ramsey. Biotechnol. Lett. 2008, 30, 1861–1866. [Google Scholar] [CrossRef]
  31. Diatchenko, L.; Lau, Y.F.; Campbell, A.P.; Chenchik, A.; Moqadam, F.; Huang, B.; Lukyanov, S.; Lukyanov, K.; Gurskaya, N.; Sverdlov, E.D. Suppression subtractive hybridization: A method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc. Natl. Acad. Sci. USA 1996, 93, 6025–6030. [Google Scholar] [CrossRef]
  32. Kononen, J.; Bubendorf, L.; Kallioniemi, A.; Barlund, M.; Schraml, P.; Leighton, S.; Torhorst, J.; Mihatsch, M.J.; Sauter, G.; Kallioniemi, O.P. Tissue microarrays for high-throughput molecular profiling of tumor specimens. Nat. Med. 1998, 4, 844–847. [Google Scholar] [CrossRef] [PubMed]
  33. ‘t Hoen, P.A.C.; Ariyurek, Y.; Thygesen, H.H.; Vreugdenhil, E.; Vossen, R.H.A.M.; de Menezes, R.X.; Boer, J.M.; van Ommen, G.J.B.; den Dunnen, J.T. Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36, e141. [Google Scholar] [CrossRef] [PubMed]
  34. Ko, J.H.; Kim, H.T.; Hwang, I.; Han, K.H. Tissue-type-specific transcriptome analysis identifies developing xylem-specific promoters in poplar. Plant Biotechnol. J. 2012, 10, 587–596. [Google Scholar] [CrossRef] [PubMed]
  35. Geng, L.; Duan, X.; Liang, C.; Shu, C.; Song, F.; Zhang, J. Mining tissue-specific contigs from peanut (Arachis hypogaea L.) for promoter cloning by deep transcriptome sequencing. Plant Cell Physiol. 2014, 55, 1793–1801. [Google Scholar] [CrossRef] [PubMed]
  36. Meer, I.M.V.D.; Spelt, C.E.; Mol, J.N.M.; Stuitje, A.R. Promoter analysis of the chalcone synthase (chs A) gene of Petunia hybrida: A 67 bp promoter region directs flower-specific expression. Plant Mol. Biol. 1990, 15, 95–109. [Google Scholar] [CrossRef]
  37. Sablowski, R.W.; Moyano, E.; Culianez-Macia, F.A.; Schuch, W.; Martin, C.; Bevan, M. A flower-specific Myb protein activates transcription of phenylpropanoid biosynthetic genes. EMBO J. 1994, 13, 128–137. [Google Scholar] [CrossRef]
  38. Han, Y.Y.; Ming, F.; Wang, J.W.; Ye, M.M.; Shen, D.L. A novel chalcone synthase gene from Phalaenopsis orchid that alters floral morphology in transgenic tobacco plants. Plant Mol. Bio. Rep. 2005, 23, 193–194. [Google Scholar] [CrossRef]
  39. Brocard, I.; Charlot, F.; Teoule, E.; Guerche, P. Petal Specific Promoter and Method for Obtaining Plants Having (P). European: Institut National De La Recherche Agronomique. JP 2001517450-A 309-0CT-2001, 2001. [Google Scholar]
  40. Fu, W.; Chen, D.; Pan, Q.; Li, F.; Zhao, Z.; Ge, X.; Li, Z. Production of red-flowered oilseed rape via the ectopic expression of Orychophragmus violaceus OvPAP2. Plant Biotechnol. J. 2018, 16, 367–380. [Google Scholar] [CrossRef] [Green Version]
  41. Azuma, M.; Morimoto, R.; Hirose, M.; Morita, Y.; Hoshino, A.; Iida, S.; Oshima, Y.; Mitsuda, N.; Ohme-Takagi, M.; Shiratake, K. A petal-specific InMYB1 promoter from Japanese morning glory: A useful tool for molecular breeding of floricultural crops. Plant Biotechnol. J. 2016, 14, 354–363. [Google Scholar] [CrossRef]
  42. Young, C.S.; Werner, C.P. Infection routes for Sclerotinia sclerotiorum in apetalous and fully petalled winter oilseed rape. Plant Pathol. 2012, 61, 730–738. [Google Scholar] [CrossRef]
  43. Dong-Hui, F.U.; Jiang, L.Y.; Mason, A.S.; Xiao, M.L.; Zhu, L.R.; Li-Zhi, L.I.; Zhou, Q.H.; Shen, C.J.; Huang, C.H. Research progress and strategies for multifunctional rapeseed: A case study of China. J. Integr. Agric. 2016, 15, 1673–1684. [Google Scholar]
  44. Doring, T.F.; Skellern, M.; Watts, N.; Cook, S.M. Colour choice behaviour in the pollen beetle Meligethes aeneus (Coleoptera: Nitidulidae). Physiol. Entomol. 2012, 37, 360–378. [Google Scholar] [CrossRef]
  45. Zhang, H.; Li, H.C.; Miao, X.X. Feasibility, limitation and possible solutions of RNAi-based technology for insect pest control. Insect Sci. 2013, 20, 15–30. [Google Scholar] [CrossRef] [PubMed]
  46. Lilley, C.J.; Davies, L.J.; Urwin, P.E. RNA interference in plant parasitic nematodes: A summary of the current status. Parasitology 2012, 139, 630–640. [Google Scholar] [CrossRef]
  47. Cai, Q.; Qiao, L.; Wang, M.; He, B.; Lin, F.M.; Palmquist, J.; Huang, H.; Jin, H. Plants send small RNAs in extracellular vesicles to fungal pathogen to silence virulence genes. Science 2018, 360, 1126–1129. [Google Scholar] [CrossRef] [Green Version]
  48. Cai, Q.; He, B.; Kogel, K.H.; Jin, H. Cross-kingdom RNA trafficking and environmental RNAi – nature’s blueprint for modern crop protection strategies. Curr. Opin. Microbiol. 2018, 46, 58–64. [Google Scholar] [CrossRef]
  49. Patel, R.K.; Jain, M. NGS QC Toolkit: A Toolkit for Quality Control of Next Generation Sequencing Data. Plos ONE 2012, 7. [Google Scholar] [CrossRef]
  50. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [Green Version]
  51. Jin, J.P.; Tian, F.; Yang, D.C.; Meng, Y.Q.; Kong, L.; Luo, J.C.; Gao, G. PlantTFDB 4.0: Toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2017, 45, D1040–D1045. [Google Scholar] [CrossRef] [Green Version]
  52. Jin, J.P.; He, K.; Tang, X.; Li, Z.; Lv, L.; Zhao, Y.; Luo, J.C.; Gao, G. An Arabidopsis Transcriptional Regulatory Map Reveals Distinct Functional and Evolutionary Features of Novel Transcription Factors. Mol. Biol. Evol. 2015, 32, 1767–1773. [Google Scholar] [CrossRef] [Green Version]
  53. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  54. Sun, F.M.; Fan, G.Y.; Hu, Q.; Zhou, Y.M.; Guan, M.; Tong, C.B.; Li, J.N.; Du, D.Z.; Qi, C.K.; Jiang, L.C.; et al. The high-quality genome of Brassica napus cultivar ‘ZS11′ reveals the introgression history in semi-winter morphotype. Plant J. 2017, 92, 452–468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Higo, K.; Ugawa, Y.; Iwamoto, M.; Korenaga, T. Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27, 297–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Clough, S.J.; Bent, A.F. Floral dip: A simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998, 16, 735–743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Sung, Y.C.; Lin, C.P.; Chen, J.C. Optimization of virus-induced gene silencing in Catharanthus roseus. Plant Pathol. 2015, 63, 1159–1167. [Google Scholar] [CrossRef]
  58. Jefferson, R.A.; Kavanagh, T.A.; Bevan, M.W. GUS fusions: Beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 1987, 6, 3901–3907. [Google Scholar] [CrossRef]
  59. Bradford, M.M. A Rapid Method for the Quantitation of Microgram Quantities of Protein Utilizing the Principle of Protein-Dye Binding. Anal. Biochem. 1976, 72, 248–254. [Google Scholar]
Figure 1. Flower-specific gene screening procedure and Gene Ontology (GO) classification and enrichment of flower-specific candidates: (a) Screening procedures for flower-specific candidate genes. (b) GO classification of flower-specific candidate genes. A total of 143 genes of 249 gene candidates have GO annotations. (c) GO enrichment related with organ formation, development or specification of flower-specific candidate genes. 31 GO terms were enriched, 15 of which were related with organ formation, development or specification.
Figure 1. Flower-specific gene screening procedure and Gene Ontology (GO) classification and enrichment of flower-specific candidates: (a) Screening procedures for flower-specific candidate genes. (b) GO classification of flower-specific candidate genes. A total of 143 genes of 249 gene candidates have GO annotations. (c) GO enrichment related with organ formation, development or specification of flower-specific candidate genes. 31 GO terms were enriched, 15 of which were related with organ formation, development or specification.
Ijms 20 05949 g001
Figure 2. Expression profiles of six flower-specific candidate genes in different tissues by quantitative real-time transcription-PCR (qRT-PCR) analysis. (a) Expression patterns of six candidate genes in five different tissues. R (root), S (stem), L (leaf), F (flower bud), SQ (silique). The asterisk indicates that the expression value in the flower bud was at least three times higher than that in other tissues. (b) Expression profiles of six flower-specific candidate genes in dissected flower parts. The relative expression values were log10 transformed. Se (sepal), P (pistil), St (stamen), BP (blossomy petals), WP (wilting petals). Tissues from five plants were collected together, three technical repeats were performed and three biological repeats were conducted. The B. napus β-actin gene (AF111812) was used as a reference standard.
Figure 2. Expression profiles of six flower-specific candidate genes in different tissues by quantitative real-time transcription-PCR (qRT-PCR) analysis. (a) Expression patterns of six candidate genes in five different tissues. R (root), S (stem), L (leaf), F (flower bud), SQ (silique). The asterisk indicates that the expression value in the flower bud was at least three times higher than that in other tissues. (b) Expression profiles of six flower-specific candidate genes in dissected flower parts. The relative expression values were log10 transformed. Se (sepal), P (pistil), St (stamen), BP (blossomy petals), WP (wilting petals). Tissues from five plants were collected together, three technical repeats were performed and three biological repeats were conducted. The B. napus β-actin gene (AF111812) was used as a reference standard.
Ijms 20 05949 g002
Figure 3. Conserved cis-elements in FSP046 and FSP061 promoters. Motif positions are indicated relative to the start codon.
Figure 3. Conserved cis-elements in FSP046 and FSP061 promoters. Motif positions are indicated relative to the start codon.
Ijms 20 05949 g003
Figure 4. Histochemical staining of FSP046::GUS transgenic A. thaliana plants. There is no GUS staining in the two-leaf stage (a), six-leaf stage (b), rosette leaf (c), flower bud (d), silique (i), or seed (j). GUS activity driven by the FSP046 promoter was observed in flowers one day before blooming (e), flowers at the blooming day (f), petals one day after blooming (g), and petals two days after blooming (h). Moreover, there is the highest GUS activity in the stamen/stigma (f), and less activity in the petals (fh). Bar = 1 mm
Figure 4. Histochemical staining of FSP046::GUS transgenic A. thaliana plants. There is no GUS staining in the two-leaf stage (a), six-leaf stage (b), rosette leaf (c), flower bud (d), silique (i), or seed (j). GUS activity driven by the FSP046 promoter was observed in flowers one day before blooming (e), flowers at the blooming day (f), petals one day after blooming (g), and petals two days after blooming (h). Moreover, there is the highest GUS activity in the stamen/stigma (f), and less activity in the petals (fh). Bar = 1 mm
Ijms 20 05949 g004
Figure 5. GUS gene expression in transgenic A. thaliana plants. (a) GUS gene expression patterns in five different tissues (root, stem, leaf, flower, silique) under FSP046 and FSP061 promoter direction. The A. thaliana Actin gene (AT3G18780) was used as reference gene. (b) GUS activity in six different tissues (root, stem, leaf, flower, silique, and seed) under FSP046 and FSP061 promoter direction. GUS activity was calculated as pmol 4-MU per min per mg protein. 4-MU, 4-methylumbelliferone.
Figure 5. GUS gene expression in transgenic A. thaliana plants. (a) GUS gene expression patterns in five different tissues (root, stem, leaf, flower, silique) under FSP046 and FSP061 promoter direction. The A. thaliana Actin gene (AT3G18780) was used as reference gene. (b) GUS activity in six different tissues (root, stem, leaf, flower, silique, and seed) under FSP046 and FSP061 promoter direction. GUS activity was calculated as pmol 4-MU per min per mg protein. 4-MU, 4-methylumbelliferone.
Ijms 20 05949 g005
Figure 6. Histochemical staining of B. napus flowers under transient expression. GUS staining can be observed in flowers of the positive control 35S::GUS, FSP046::GUS, and FSP061::GUS, respectively (especially in their sepals, petals, and stamens). While no GUS stained could be detected in negative control (CK-, the empty Agrobacterium). Bar = 1 mm.
Figure 6. Histochemical staining of B. napus flowers under transient expression. GUS staining can be observed in flowers of the positive control 35S::GUS, FSP046::GUS, and FSP061::GUS, respectively (especially in their sepals, petals, and stamens). While no GUS stained could be detected in negative control (CK-, the empty Agrobacterium). Bar = 1 mm.
Ijms 20 05949 g006
Table 1. Transcriptome statistics of twelve different tissues from Brassica napus.
Table 1. Transcriptome statistics of twelve different tissues from Brassica napus.
RootLeafBudSiliqueStamenPistilBlossomy Petal *Wilting Petal **StemSepalOvulePericarp
Clean data47,855,96751,669,76550,250,95047,583,78920,822,44323,301,85324,872,04423,779,16918,987,58525,334,03823,375,36025,707,408
All data mapping to genome34,126,090 35,693,473 34,969,63632,285,60118,163,41720,477,66820,852,72217,021,12916,418,56521,868,34220,091,12218,195,703
The percent of all data mapping to genome71.31%69.08%69.59%67.85%87.23%87.88%83.84%71.58%86.47%86.32%85.95%70.78%
Unique mapping data25,997,36226,362,59725,777,94023,914,69215,536,38617,602,36417,884,29512,917,63514,049,82918,288,12317,027,68913,331,423
The percent of unique mapped data54.32%51.02%51.30%50.26%74.61%75.54%71.91%54.32%73.99%72.19%72.84%51.86%
* The abbreviation of blossomy petal is BP; ** the abbreviation of wilting petal is WP.

Share and Cite

MDPI and ACS Style

Li, Y.; Dong, C.; Hu, M.; Bai, Z.; Tong, C.; Zuo, R.; Liu, Y.; Cheng, X.; Cheng, M.; Huang, J.; et al. Identification of Flower-Specific Promoters through Comparative Transcriptome Analysis in Brassica napus. Int. J. Mol. Sci. 2019, 20, 5949. https://doi.org/10.3390/ijms20235949

AMA Style

Li Y, Dong C, Hu M, Bai Z, Tong C, Zuo R, Liu Y, Cheng X, Cheng M, Huang J, et al. Identification of Flower-Specific Promoters through Comparative Transcriptome Analysis in Brassica napus. International Journal of Molecular Sciences. 2019; 20(23):5949. https://doi.org/10.3390/ijms20235949

Chicago/Turabian Style

Li, Yan, Caihua Dong, Ming Hu, Zetao Bai, Chaobo Tong, Rong Zuo, Yueying Liu, Xiaohui Cheng, Mingxing Cheng, Junyan Huang, and et al. 2019. "Identification of Flower-Specific Promoters through Comparative Transcriptome Analysis in Brassica napus" International Journal of Molecular Sciences 20, no. 23: 5949. https://doi.org/10.3390/ijms20235949

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