RNA-Sequencing Analysis Reveals L-theanine Regulating

5 days ago - L-theanine, a unique amino acid in tea leaves, is known to have beneficial effects on stress relief, tumor suppression, prevention of ...
0 downloads 0 Views 1MB Size
Subscriber access provided by Iowa State University | Library

Bioactive Constituents, Metabolites, and Functions

RNA-Sequencing Analysis Reveals L-theanine Regulating Transcriptional Rhythm Alteration in Vascular Smooth Muscle Cells Induced by Dexamethasone zhongwen Xie J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b05057 • Publication Date (Web): 27 Jan 2019 Downloaded from http://pubs.acs.org on January 27, 2019

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 36

Journal of Agricultural and Food Chemistry

1

TITLE: RNA-Sequencing Analysis Reveals L-theanine Regulating Transcriptional

2

Rhythm Alteration in Vascular Smooth Muscle Cells Induced by Dexamethasone

3

AUTHORSHIP:

4

Ruru Wang † §, Menzhao Xiao † §, Yujing Zhang † §, Chi-Tang Ho§ ‡ , Xiaochun Wan † §

5

Daxiang Li*,†,§,and Zhongwen Xie *,†,§

6

†State Key Laboratory of Tea Plant Biology and Utilization, School of Tea and Food

7

Sciences and Technology and §International Joint Laboratory on Tea Chemistry and

8

Health Effects of Ministry of Education, Anhui Agricultural University, Hefei, Anhui

9

230036, PR China

10

‡Department of Food Science, Rutgers University, 65 Dudley Road, New Brunswick,

11

New Jersey 08901-8520, United States

12

*

13

+86-551-65786153

Correspondences:

[email protected]

or

[email protected];

ACS Paragon Plus Environment

Tel.:

Journal of Agricultural and Food Chemistry

14

ABSTRACT: L-theanine, a unique amino acid in tea leaves, is known to have

15

beneficial effects on stress relief, tumor suppression, prevention of hypertension

16

and cardiovascular diseases (CADs). The disruption of the circadian rhythm has been

17

implied in the pathogenesis of CADs. However, it is unknown whether the

18

L-theanine has the modulatory effect on the vascular circadian rhythm. In this

19

research, we have established a circadian gene expression model in rat vascular

20

smooth muscle cells (VSMCs) by dexamethasone induction. L-Theanine treatment

21

enhanced the expression amplitude of clock genes including Bmal1, Cry1, Rev-erbα

22

and Per2. Moreover, pairwise comparisons of the RNA-seq data showed that

23

L-theanine is able to up-regulate a ray of the rhythm genes, and differentially

24

expressed genes (DEGs) that are involved in vasoconstriction and actin cytoskeleton

25

regulation pathways. Our data suggest that L-theanine changes the circadian gene

26

rhythm involving in the process of vascular smooth muscle restructure.

27 28

KEYWORDS: L-theanine, Biological rhythms, RNA-seq, Circadian genes, VSMCs

29

1

ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

Journal of Agricultural and Food Chemistry

31

INTRODUCTION

32

L-Theanine is a unique non-protein amino acid naturally found in tea plants. It

33

supplies the umami taste of green tea infusion.1,

34

reported that L-theanine has beneficial effects on attenuating the brain and liver

35

damage, improving the immune ability of immune cells, reducing the apoptosis rate

36

of nerve cells, inhibiting the tumor formation and cancer cell proliferation, and

37

reducing blood pressure.3,4 However, little is known whether L-theanine regulates

38

biological clock of smooth muscle cells.

39

Biological clock, which is also known as circadian rhythm, is the cyclical fluctuation

40

of biological behavior and physiological phenomena around 24 hours in organism. In

41

mammals, the central clock is located in the suprachiasmatic nucleus (SCN) of the

42

brain hypothalamus, and peripheral clocks are located in most mammalian

43

peripheral cells. The central clock is entrained by light/dark cycles, whereas

44

peripheral clocks are entrained by feeding cycles. Peripheral organs have their

45

independent autonomous circadian clocks but are synchronized by the SCN via

46

neural and endocrine signaling pathways. The Clock and Bmal1 are core circadian

47

genes, and also are transcriptional factors, which hetero-dimerize together to

48

activate expression of Period (Per1/2/3) and Cryptochrome (Cry1/2) via the E/E’box

49

sequence. PER/CRY protein complexes then inhibit the transcription-activating role

50

of ClOCK/BMAL1, decreasing Per and Cry expression, thereby forming the negative

51

feedback loop.5 Additional circadian clock components such as Rev-erbα and Rorα

2

2

ACS Paragon Plus Environment

Previous researches have

Journal of Agricultural and Food Chemistry

52

also contribute to maintaining proper circadian rhythms and are thought to have a

53

role in establishing the night-time transcriptome.6

54

Many cardiovascular physiological functions have obvious circadian rhythm, and

55

many cardiovascular diseases, such as myocardial ischemia, myocardial infarction,

56

also show rhythmically and these diseases usually occur with high frequency in the

57

morning.7 This suggests that cardiovascular physiology and pathology are closely

58

related to circadian clock.8, 9 Indeed, clock rhythm is found in mainly all tissues and

59

various types of cells including endothelial cells, smooth muscle cells and vascular

60

stem cells in the vascular system. This vascular clock rhythm physiologically

61

contributes to regulate the daily vascular function, and plays a pivotal role in the

62

pathophysiology of vascular diseases.8 Moreover, vascular circadian rhythm inserts

63

regulatory effects on vessel contraction, haemodynamics, inflammatory response

64

and endothelium-derived NO synthesis and release, resulting in plaque formation

65

and instability which participated in CADs. 10,11

66

Recent findings showed that nutrients reset peripheral circadian oscillation and the

67

local clock genes govern downstream metabolic pathways.12, 13 A number of studies

68

also showed that consumption of beneficial components at right times would have

69

similar health effects as medication administered at specific times in

70

chronopharmacology. 14,15 Nutrients can modulate circadian rhythmicity, and some

71

studies reported that natural components administration changes biological rhythm,

72

and increasing the amplitude of rhythm leads to resist obesity, strengthen immunity

73

and make the body better adapt to the environment.16,17 These natural components 3

ACS Paragon Plus Environment

Page 4 of 36

Page 5 of 36

Journal of Agricultural and Food Chemistry

74

are able to alter signaling pathways that impact the molecular oscillator in

75

peripheral cells. The molecular clock inside the cells of the cardiovascular system

76

allows it to respond appropriately to the stimulation of the external nutrition.18 The

77

attenuation of this molecular clock's oscillation can affect the timely and effective

78

reaction of the cardiovascular to the external environment, thus affecting the

79

occurrence and development of some CADs.19 Previous studies majorly focused on

80

the effects of drugs and bioactive substances on the binding of Clock/Bmal1

81

heterodimer to E-box of target genes promoter region and insert the effects of

82

nutrients on the biological clock. 20,21

83

Changes in lifestyle, dietary habits and increased life expectancy have led to

84

increase

85

hypercholesterolemia, obesity, diabetes and hypertension over recent decades.

86

According to the World Health Organization (WHO), CADs are the leading cause of

87

death worldwide. Growing evidence links circadian alterations to metabolic

88

disorders. And clock mutation reduces circadian pacemaker amplitude in mice.

89

And alteration of cardiovascular circadian genes results in CADs. Disruption of

90

smooth muscle Bmal1 leads to alteration in rhythmic blood pressure and protection

91

of abdominal aortic aneurysm,23,24 endothelium specific disruption of the circadian

92

clock exacerbates the thrombogenic response and changes blood pressure, and

93

vascular transplantation of Bmal1-KO mice to WT mice leads to atherosclerosis.25

94

These results indicate that modulation of circadian rhythms in vascular system

widespread

presence

of

metabolic

4

ACS Paragon Plus Environment

syndrome,

such

as

22

Journal of Agricultural and Food Chemistry

95

might be a new strategy to prevention of CADs. In addition, L-theanine treatment

96

decreases blood pressure in human study.4

97

Furthermore, VSMCs are the major components in a vessel wall, and play an

98

important role in cardiovascular pathophysiology. However, little is known about

99

whether L-theanine modulates circadian rhythm in vasculature and if so, how the

100

L-theanine alters circadian genes expression.

101

In the present study, a circadian gene expression model using rat VSMCs by

102

dexamethasone induction was established. We then investigated how L-theanine

103

addition changed the rhythmic expression of major clock genes. A RNA-seq analysis

104

was employed on VSMCs to dissect the transcriptomic profiles in response to the

105

L-theanine treatments at different Zeitgeber times. Gene ontology enrichment and

106

KEGG pathway analysis of RNA-seq data were employed in deciphering the network

107

of metabolic and signaling pathways involved in the regulation of vascular smooth

108

function to L-theanine induction, with specifically focused on vascular clock gene

109

expression profiles. To the best of our knowledge, this is the first RNA-seq-based

110

study to assess the effect of L-theanine on the vascular circadian rhythms and

111

circadian gene expression profiles in primary cultured VSMCs. Our next-generation

112

sequencing data contributed to the molecular effects exerted by L-theanine on

113

vascular circadian rhythms and circadian gene expression profiles .

5

ACS Paragon Plus Environment

Page 6 of 36

Page 7 of 36

Journal of Agricultural and Food Chemistry

115

■ MATERIALS AND METHODS

116

Cell Culture and Samples Collection

117

VSMCs were isolated from the rat thoracic aorta by enzymatic dissociation, as

118

described previously.26,27 Briefly, the VSMCs were cultured in Dulbecco’s Modified

119

Eagle’s Medium (DMEM, Thermo Fisher Scientific, Waltham, MA, USA) supplied

120

with 10% fetal bovine serum (FBS, Thermo Fisher Scientific, Waltham, MA, USA),

121

100 U/mL penicillin and 100 μg/mL streptomycin (Thermo Fisher Scientific, Waltham,

122

MA, USA) at 37 °C with 5% CO2. When reaching 70-80% confluence, the cells were

123

trypsinized with Trypsin-EDTA Solution (Thermo Fisher Scientific, Waltham, MA, USA)

124

for 2 min, then washed in PBS and resuspended in DMEM growth medium. Finally,

125

the cells were passages with one to three ratio for subculture. The cells between the

126

fifth and tenth passages were used in all experiments. The cells were plated in 35

127

mm dishes for growth. The day before the experiment, cells at approximately 90%

128

of confluence were incubated with serum free starvation DMEM for 24 hours.

129

Subsequently, the cells were treated with medium containing dexamethasone (100

130

nM, Sigma, St. Louis, MO, USA) for 15 min and then changed back to a serum free

131

DMEM with or without L-theanine (40 µM, Sigma, St. Louis, MO, USA) for the

132

various time points when the cells samples were collected. The timing of the

133

beginning dexamethasone induction was defined as Zeitgeber time 0 (ZT0), and cells

134

were harvested for RNA extraction with 4 hour interval for two days at ZT0, ZT4, ZT8,

135

ZT12, ZT16, ZT20, ZT24, ZT28, ZT32, ZT36, ZT40, ZT44 and ZT48. And the cells were

136

immediately fixed with TRIzol™ (Thermo Fisher Scientific, Waltham, MA, USA). 6

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

137

RNA Isolation, cDNA Synthesis, and Real Time PCR

138

Total RNA was extracted using RNA isolator according to the protocols of

139

manufacturers (Vazyme Biotech, Nanjing, China). Reverse transcription was

140

conducted using first strand cDNA synthesis kit (Vazyme Biotech), and the Real-time

141

PCR was performed using qPCR SYBR Green Master Mix kit (Vazyme Biotech)

142

following the method described previously.27 Primer sequences were designed for

143

rat and listed in Supplementary Table S1.

144

RNA Sequencing

145

Two groups, control (C) and L-theanine treatment (T) were included for RNA-seq

146

analysis, Four time points (ZT0, ZT12, ZT24, ZT36,) were chosen for each of group

147

and three biological replicates (RI, RII, and RIII) were performed. We used 2 µg of

148

total RNA from each of samples to prepare the TruSeq library. The NEB Next TruSeq

149

RNA Sample Preparation Kit was employed to prepare barcoded cDNA (E7530S;

150

E7490S; E7335S; E7500S, New England Biolabs, Ipswich, MA, USA).28 The fragments

151

per Kilobase of transcript per Million fragments mapped (FPKM) was calculated by

152

dividing the read count of each transcriptional model with its length and scaling the

153

total per sample to one million, and was used to indicate the expression levels of

154

each sample. A gene expression value in FPKM equal or greater than 1 was

155

considered to be expressed in the sample.

156

Transcriptome Data Analysis

157

The aligning of the raw sequences data was analyzed by proper quality control.

158

FastQC and Trimmomatic v. 0.33. were used for Quality control (QC) and trimming. 7

ACS Paragon Plus Environment

Page 8 of 36

Page 9 of 36

Journal of Agricultural and Food Chemistry

159

TopHat/Bowtie2 was employed to map quality checked reads after QC . DESeq R

160

package (1.10.1) was applied to analyze differential expression of two groups.

161

DESeq provides statistical tools for evaluating differential expression in digital gene

162

expression data using a model based on the negative binomial distribution. The

163

Benjamini and Hochberg’s approach was used to adjust the resulting P values for

164

controlling the false discovery rate. Genes with a DESeq adjusted P-value < 0.05

165

were assigned as differentially expressed.29,30

166

Real Time -PCR

167

To validate the reliability of RNA-seq results, core genes (Bmal1, Npas2) and clock

168

target genes (Cry1, Per2, RORα, Dbp) and related to vascular contraction genes

169

(Rock2, CPI-17) were selected for real time-PCR validation.

170

Western Blot Analysis

171

Western blot was performed following the method described previously.27 In brief,

172

the frozen cells were homogenized using a 2×SDS homogenization buffer. Equal

173

amounts of denatured proteins were loaded and separated by SDS-PAGE gels, and

174

the proteins were transferred to nitrocellulose membranes with constant current

175

model. The membranes were blocked with 5% skimmed milk in PBS-T buffer for 1 h

176

and were further incubated with the first antibodies ROCK2, CPI-17 (Santa Cruz

177

Biotechnology, Santa Cruz, CA, USA), and β-Actin (Proteintech, Wuhan, China) at 4

178

℃ overnight, then incubated with appropriate secondary antibodies (Proteintech)

179

for 1 h at room temperature. Enhanced chemiluminescent (ECL) reagent (Vazyme 8

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

180

Biotech) was used to detect protein bands, which were captured and analyzed using

181

the ChemicDoc Imaging System (Bio-Rad Laboratories, Hercules, CA, USA) with

182

ImageLab system (Bio-Rad Laboratories), respectively.

183

Statistical Analysis

184

GraphPad PRISM 6 was used for statistical analyses and data visualization. Data are

185

presented as mean ± SEM. Multiple group statistical comparisons were determined

186

by one-way or two-way ANOVA with Turkey’s or Dunnett’ tests. P value less than

187

0.05 was considered to be statistically significant. Peak times of circadian gene

188

expression were calculated by fitting a sinusoidal function to the real time-PCR data.

9

ACS Paragon Plus Environment

Page 10 of 36

Page 11 of 36

Journal of Agricultural and Food Chemistry

190

RESULTS

191

Desamethasone Synchronizes Clock Genes mRNA Rhythmic Expression in VSMCs

192

In the present study, the rat primary cultured VSMCs were used to determine

193

whether several core circadian genes were induced rhythmic expression by

194

dexamethasone. Our real time PCR data revealed that brief treatment of VSMCs

195

with dexamethasone resulted in the robust rhythmic expression of Bmal1, Cry1,

196

Rev-erbα, and Per2 for at least 2 circadian cycles. Bmal1 mRNA expression level

197

initiated a short-term increase and then decrease showing an around 24 hour

198

rhythm, which consisted of peaks at ZT12 and ZT36. The expression Levels of Cry1

199

mRNA showed peaks at ZT8 and ZT32. Robust approximate cycling of Rev-erbα and

200

Per2 mRNA expression was also observed, with mRNA expression levels

201

accumulating antiphase to Bmal1 RNA cycles. However, another core clock gene

202

Clock mRNA did not show rhythmic expression (Figure 1).

203

Transcriptome Sequencing and Assembly

204

Using this VSMCs circadian model, we attempted to figure out genes that responded

205

to L-theanine induction. To obtain a comprehensive circadian gene transcriptional

206

profile of the induced VSMCs, collecting sample at the right time point is crucial. In

207

general, gene expressions and regulations precede the presence of enzymes and

208

their products. Based on this rational and the several circadian genes expression

209

pattern showing in Figure 1, VSMCs samples at ZT0, ZT12, ZT24, ZT36 with or

210

without L-theanine induction were chosen to explore the potential rhythmic genes

211

response to L-theanine by RNA-seq analysis. The cDNA libraries for each sample 10

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

212

were constructed and sequenced. The sequencing data showed that around 25

213

million reads were obtained for each sample. After data filtering and stringent

214

quality control, we mapped clean reads to reference genome using Taphat2. Among

215

all reads, 83~95% of reads are mapped to rat genome, and the uniformity of the

216

mapping result for each sample suggests that the samples are reliable and

217

comparable. The details of mapping data are shown in Table 1. Our data suggested

218

that the assembled results were qualified for further analyses. The levels of gene

219

expression were quantified and expressed as fragments per kilobase of transcript

220

per million mapped reads (FPKM). Based on the gene expression information, we

221

preformed histogram to show the distribution of the gene expression level of each

222

sample. In addition, we observed the dispersion of the gene distribution. In order to

223

show the gene expression level under different FPKM value, we calculated the gene

224

expression level under eight different FPKM ranges (0 < FPKM < 0.1, 0.1=< FPKM< 1,

225

1=< FPKM < 5, 5 =< FPKM