Comparative Transcriptome Analysis Revealing the Effect of Light on

Mar 9, 2018 - To study the mechanism of anthocyanin-biosynthesis regulation, we examined light-regulated gene expression involved in anthocyanin ...
0 downloads 0 Views 2MB Size
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