Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE
Biotechnology and Biological Transformations
Comparative Transcriptome Analysis Reveals Effects of Light on Anthocyanin Biosynthesis in Purple Grain of Wheat Fang Wang, Yu Xiu Dong, Xiao Zhen Tang, Tian Li Tu, Bin Zhao, Na Sui, Daolin Fu, and Xian Sheng Zhang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b05435 • Publication Date (Web): 09 Mar 2018 Downloaded from http://pubs.acs.org on March 9, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 46
Journal of Agricultural and Food Chemistry
1
Comparative Transcriptome Analysis Reveals Effects of
2
Light on Anthocyanin Biosynthesis in Purple Grain of
3
Wheat
4
Fang Wang1§, Yu Xiu Dong1§, Xiao Zhen Tang2§, Tian Li Tu1, Bin Zhao1, Na Sui3,
5
Dao Lin Fu1, Xian Sheng Zhang1*
6
1
7
Agricultural University, Tai’an, 271018, Shandong, China
8
2
9
271018 , Shandong, China
State Key Laboratory of Crop Biology, College of Life Sciences, Shandong
College of Food Science and Engineering,Shandong Agricultural University, Tai’an,
10
3
11
Shandong Normal University, Jinan, 250014, China
12
§
13
Corresponding Author
14
*E-mail:
[email protected] Shandong Provincial Key Laboratory of Plant Stress, College of Life Science,
These authors contributed equally to this work.
15
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
16
ABSTRACT
17
To study mechanism of anthocyanin biosynthesis regulation, we examined
18
light-regulated gene expression involved in anthocyanin biosynthesis in purple grain
19
of wheat. Ten kinds of anthocyanins were identified from a purple-grain wheat
20
cultivar by HPLC ESI-MS/MS analysis. The libraries constructed from total RNA of
21
purple grain under light (L) or dark (D) conditions for 15 and 20 days were sequenced.
22
In total, 1,874 differentially expressed genes (DEGs) in L20 vs. L15, 1,432 DEGs in
23
D20 vs. D15, 862 DEGs in D15 vs. L15, and 1,786 DEGs in D20 vs. L20 were
24
identified. DEG functional enrichments suggested that light signal transduction is
25
critical to anthocyanin biosynthesis. 911 DEGs were referred to as LDEGs
26
(light-regulated DEGs) involved in a number of genes in anthocyanin biosynthesis,
27
transcription regulation, sugar- and calcium-signaling pathways and hormone
28
metabolisms. These findings lay the foundation for future studies on regulatory
29
mechanism of anthocyanin biosynthesis in purple grain of wheat.
30
KEY WORDS: anthocyanin biosynthesis, light, RNA-seq, purple-grain wheat,
31
transcriptomic analysis
32 33 34 35 36 37
ACS Paragon Plus Environment
Page 2 of 46
Page 3 of 46
Journal of Agricultural and Food Chemistry
38
INTRODUCTION
39
Anthocyanins are a class of secondary metabolites that contribute to the red, blue, and
40
purple colors in plants. In flowers, these pigments attract pollinators, and in fruit skin
41
they attract animals to aid in seed dispersal. Anthocyanins are also important for
42
maintaining human health by preventing cardiovascular disease, carcinogenesis,
43
inflammation, and many other pathological states.
44
Most of the genes involved in anthocyanin biosynthesis and its regulation have
45
been isolated and characterized in a number of plant species. 1 The genes involved in
46
the early step of dihydroflavonol biosynthesis include phenylalanine ammonia lyase
47
(PAL), chalcone synthase (CHS), chalcone isomerase (CHI), and flavanone
48
3-hydroxylase (F3H), and those involved in successive reactions in anthocyanin
49
production include dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS),
50
and
51
biosynthesis is regulated by the combined action of R2R3-MYB and bHLH
52
transcription factors (TFs), as well as WD40-repeat proteins.3-5
UDP-glucose, flavonoid 3-O-glucosyltran-ferase
(UFGT).2
Anthocyanin
53
The purple colors are not natural for seeds of common hexaploid wheat.
54
Cyanidin-3-glucoside and peonidine-3-glucoside are the main anthocyanins in the
55
grain coat of purple-grain wheat.6 Most of the structural genes for anthocyanin
56
biosynthesis in wheat have been identified. In previous studies, two closely linked,
57
highly homologous genes encoding phenylalanine ammonia-lyase (PAL1 and PAL2)
58
were isolated from the wheat genomic library, 7 and the loci for PAL, CHS, and CHI in
59
wheat were determined using Southern blot hybridization analyses. 8 Analyses of four
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
60
early flavonoid biosynthesis genes (CHS, CHI, F3H, and DFR) in wheat showed that
61
they were expressed at higher levels in the grain coat tissue of red-grained lines than
62
in white-grained ones.9 Three homoeologous full-length CHI copies were isolated and
63
precisely mapped to the long arms of chromosomes 5A, 5B, and 5D.10 F3H-A1,
64
F3H-B1, and F3H-D1 were mapped to chromosomes 2A, 2B, and 2D of wheat, and
65
were shown to be highly expressed in red grains and coleoptiles.11 Three copies of
66
DFR were mapped to homoeologous group 3 chromosomes,12 and five copies of ANS
67
were mapped to chromosomes 6A, 6B, and 6D.13 In addition, a partial UFGT
68
sequence has been isolated from the wheat genome.14
69
Light signaling plays a pivotal role in controlling plant morphogenesis, metabolism,
70
and development. Light is sensed by plants through several classes of photoreceptors
71
such as the red- and far-red light-sensing phytochromes, the blue/ultraviolet
72
(UV)-A-perceiving cryptochromes and phototropins, as well as the UV-B-sensing
73
photoreceptor UVR8.15 Upon light activation, photoreceptors induce developmental
74
responses such as seedling deetiolation, phototropism, and flowering induction, as
75
well as affecting metabolism. In Arabidopsis, anthocyanins accumulate in light-grown
76
plants but not in dark-grown ones, implying that light is required for anthocyanin
77
biosynthesis. However, little is known about the regulatory genes that control the
78
biosynthesis of anthocyanins in response to light signals.16 Wheat (Triticum aestivum
79
L.) cultivars possessing purple grain are thought to be more nutritious to human
80
because of their high anthocyanin contents in their grain coats. In this study, the
81
components of anthocyanin in the purple-grain cultivar (Luozhen No.1) were
ACS Paragon Plus Environment
Page 4 of 46
Page 5 of 46
Journal of Agricultural and Food Chemistry
82
examined by HPLC ESI-MS/MS analysis. Because light is critical to regulate the
83
production of anthocyanins in grain coats, we made the transcriptome analysis on the
84
coat tissues of purple-grain wheat from plants grown under light conditions for 15 and
85
20 days and under dark conditions for 15 and 20 days, respectively. The results of this
86
study might provide important information for understanding mechanism of
87
light-mediated anthocyanin biosynthesis in purple-grain wheat.
88
MATERIALS AND METHODS
89
Plant Materials and Growth Conditions. The purple-grain wheat cv. Luozhen No.1
90
was used for RNA-seq analyses. The wheat plants were cultivated in a field at the
91
Experimental Station of Shandong Agricultural University. Whole ears were covered
92
with two layers of dark paper bags for dark treatment. The grains were harvested at 0,
93
10, 15, 20, 25, and 30 DAP. The grain coats were isolated and collected, frozen in
94
liquid nitrogen, and stored at −80℃ for RNA extraction.
95
Total anthocyanins analysis. Total anthocyanins were measured using a
96
spectrophotometric differential pH method. 1.0 g grains were crushed into powder
97
and extracted separately with 2 mL of pH 1.0 buffer containing 50 mM KCl and 150
98
mM HCl, as well as 2 mL of pH 4.5 buffer containing 400 mM sodium acetate and
99
240 mM HCl. The mixtures were centrifuged at 10,000g for 3 min. Supernatants were
100
collected and diluted for direct measurement of absorbance at 510 nm. Total
101
anthocyanin content was calculated using the following equation:
102
Amount (mg/g) = (A510 at pH 1.0−A510 at pH 4.5)×445.2/29,600×dilution factor.
103
The number 445.2 is the molecular mass of cyanidin-3-glucoside and 29,600 is its
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
104
molar absorptivity (ε) at 510 nm. Each sample was analyzed in triplicate and the
105
results were expressed as the average of the three measurements.
106
Anthocyanin Extraction and HPLC-ESI-MS/MS Analysis. Coats of purple-grain
107
wheat were used for the anthocyanin components analysis. About 1.0 g grain coats
108
were weighed, and then, 5 mL of acidified ethanol (ethanol (95%) and HCl (1.0 N),
109
85:15, v/v) was added. The solution was mixed and adjusted to pH 1 with 4 N HCl,
110
and shaken for 20 min. The solution was centrifuged at 4,800 rpm for 5 min. The
111
extraction was repeated four times. The supernatant was poured into a 50 mL
112
volumetric flask and vaporized under 40°C by a rotary evaporator. The pigments were
113
dissolved in 5 ml 0.1% methanoic acid and purified using SPE C18 column (Waters).
114
The samples were then analyzed by ACQUITY Ultra Performance Liquid
115
Chromatography (Waters, USA), equipped with mass spectrum system (Micromass
116
Quattro Ultima IMPT, Waters, USA). The chromatographic separation was performed
117
on Sun FireTM C18 column (4. 6 mm × 150 mm, 5 μm, Waters, USA). The injection
118
volume was 10.0 μL. Elution was performed using mobile phase A (aqueous 0.1%
119
formic acid solution) and mobile phase B (methanol). The column oven temperature
120
was set at 40 °C. The flow rate was 0.8 mL/min. The gradient program is described as
121
follows: 0−13 min, 10-52% B; 13−15 min, 52-90% B; 15−17.5 min, 90−10% B.
122
Quantification of the different anthocyanins was based on peak areas and calculated
123
as equivalents of the standard compounds (Standard substances, purity > 90%, i.e.
124
cyanidin-3-glucoside, peonidin-3-glucoside and delphinidin-3-glucoside). Scanning
125
wavelength range of diode array detector was 200-800 nm. Data were acquired by
ACS Paragon Plus Environment
Page 6 of 46
Page 7 of 46
Journal of Agricultural and Food Chemistry
126
Masslinx 4.0 software (Micromass, Beverly, MA) with Quan-Optimize option for the
127
fragmentation study. The experimental conditions were as follows: ionization mode is
128
atmospheric pressure electrospray particle source (ESI), positive ion mode; capillary
129
voltage 1.50kV, cone voltage 50V, Ion source temperature, 120 ℃, dissolvent gas
130
temperature, 450 ℃, cone-hole gas flow, 81L / h, desolvent gas flow, 618L / h. The
131
voltage of photoelectric multiplier is 650 V. MS/MS, scan from m/z 200 to 1500; ion
132
trap, scan from m/z 200 to 1500; source accumulation, 50 ms; ion accumulation time,
133
300 ms; flight time to acquisition cell, 1 ms; smart parameter setting (SPS),
134
compound stability, 50%; trap drive level, 60%.
135
RNA extraction and library construction. Total RNAs were extracted using an
136
RNeasy Plus Micro kit (Qiagen, Valencia, CA, USA). At least 500 ng total RNAs
137
were extracted from each material used for library construction. First, mRNAs
138
extracted from each material were enriched by using oligo(dT) magnetic beads
139
(Illumina, San Diego, CA, USA). The purity and quantity of total RNAs were
140
checked. Then, mRNAs were further enriched by removing rRNAs from the total
141
RNAs. Mixed with the fragmentation buffer, the mRNAs were fragmented into short
142
fragments. Then cDNA was synthesized using the mRNA fragments as templates.
143
Short fragments were purified and resolved with EB buffer for end reparation and
144
single nucleotide A (adenine) addition. After that, the short fragments were ligated
145
with adapters. After agarose gel electrophoresis, the suitable fragments were selected
146
for the PCR amplification as templates. During the quality control (QC) steps, Agilent
147
2100 Bioanaylzer and ABI StepOnePlus Real-Time PCR System were used for
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 46
148
quantification and qualification of the sample libraries.
149
Illumina sequencing and data analysis. Each library was sequenced using the
150
Illumina HiSeq™ 2000 system in the Beijing Genomics Institute (Shenzhen, China).
151
Raw reads were subjected to QC that determines if a resequencing step was needed.
152
After QC, raw reads were filtered into clean reads which will be aligned to the
153
reference sequences. QC of alignment was performed to determine if resequencing is
154
needed. The alignment data was utilized to calculate distribution of reads on reference
155
genes and mapping ratio. If alignment result passed QC, we could proceed with
156
downstream analysis. RNA-seq data were obtained in three biological replicates and
157
deposited in the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress)
158
under accession number E-MTAB-5975 and E-MTAB-6398.
159
The filtered clean reads were mapped to the wheat reference genome
160
(ftp://ftp.ensemblgenomes.org/pub/plants/release-26/fasta/triticum_aestivum/dna/Triti
161
cum_aestivum.IWGSC2.26.dna.toplevel.fa.gz) using Bowtie software, and mapped to
162
the
163
(ftp://ftp.ensemblgenomes.org/pub/plants/release-26/fasta/triticum_aestivum/cds/Triti
164
cum_aestivum.IWGSC2.26.cds.all.fa.gz) using BWA software. Mismatches of less
165
than two bases were allowed in this process. According to the alignment, clean reads
166
were divided into unmapped reads, multi-position matched reads, and unique matched
167
reads. For all mapped transcripts with unique matched reads, the original digital gene
168
expression levels were calculated using fragments Per Kilobase of transcript per
169
Million fragments mapped (FPKM) method. The DEGs were identified based on a
wheat
reference
ACS Paragon Plus Environment
gene
Page 9 of 46
Journal of Agricultural and Food Chemistry
170
-fold change of ≥2 and diverge probability ≥0.8. Gene annotations were obtained
171
from the wheat sequence (http://plants.ensembl.org/Triticum_aestivum/Info/Index).
172
Genes were classified using MapMan software.
173
The gene ontology (GO) analysis was conducted using the singular enrichment
174
analysis tool (http://bioinfo.cau.edu.cn/agriGO/analysis.php). GO terms with p-value
175