Secondary Metabolite Accumulation Associates with Ecological

6 days ago - Cynomorium songaricum Rupr. is a rare root-parasitic plant distributed in the desert ecosystem. Little is known about the role of endophy...
0 downloads 3 Views 1MB Size
Subscriber access provided by University | of Minnesota Libraries

Agricultural and Environmental Chemistry

Secondary Metabolite Accumulation Associates with Ecological Succession of Endophytic Fungi in Cynomorium songaricum Rupr. Jin-Long Cui, Yan-Yan Zhang, Vinod Vijayakumar, Gang Zhang, Meng-Liang Wang, and Jun-Hong Wang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01737 • Publication Date (Web): 11 May 2018 Downloaded from http://pubs.acs.org on May 12, 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 43

Journal of Agricultural and Food Chemistry

1

Secondary Metabolite Accumulation Associates with Ecological Succession of Endophytic

2

Fungi in Cynomorium songaricum Rupr.

3

Jin-Long Cui,*,† Yan-Yan Zhang,† Vinod Vijayakumar,§ Gang Zhang,*,‡ Meng-Liang Wang,†

4

Jun-Hong Wang†

5 6





Institute of Applied Chemistry, Shanxi University, Taiyuan, Shanxi 030006 P. R. China. College of Pharmacy, Shaanxi University of Chinese Medicine, Xianyang, Shaanxi

7

712046, P. R. China.

8

§

9

Technology, The Ohio State University, Columbus, OH, 43210, USA.

College of Food, Agricultural and Environmental Sciences, Department of Food Science and

10

*

11

*

12

E-mail: [email protected]

13

*

14

E-mail: [email protected]

15

Keywords: Root-parasitic interactions, SM active ingredients, fungal endophytes, medicinal

16

plants

Correspondence: Jin-Long Cui Phone: +86-351-7016101

Gang Zhang

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

17

INTRODUCTION

18

Cynomorium songaricum Rupr. belongs to the Cynomoriaceae family, has no chlorophyll

19

but tends to parasitize the roots of Nitrariaceae plant.1 It is distributed in the peripheral

20

regions of desert with extreme environments of dry, rocky, sandy and large temperature

21

differences in Asia, Africa and Europe.2 As a rare folk medicine and healthy food, it is

22

traditionally used to treat the conditions of impotence, premature ejaculation, low sexual

23

function, spermatorrhea colic, and stomach ulcers.3 While modern phytochemical and

24

ethnopharmacological studies have indicated that its effective properties are based on diverse

25

chemical components especially active ingredients such as ursolic acid, protocatechuic acid,

26

catechins, gallic acid, flavanoids, tannins and polysaccharides (Figure 1),4 which are

27

secondary metabolites (SMs) produced during the root-parasitic plants developmental

28

process.

29

To obtain high content of active ingredients, a medicinal plant needs to be collected at the

30

appropriate time of the season which requires experience.5 Further, depending on plant age

31

and phenology the content of the active ingredient varies largely, mainly affected by

32

extracorporeal environment factors, such as soil pH, temperature and humidity.6 Besides these

33

factors, it is noteworthy that the accumulation of SMs are also regulated by plant internal

34

factors, such as endophytes.7 Which not only serve as reservoirs of novel bioactive SMs

35

themselves, but also as potential candidates for drug discovery.8 However, lots of knowledge

36

in this field is largely undetermined including mechanisms of interaction, signal transmission

37

and transformation, co-evolution and specificity between endophyte and host, and so on.9

38

Endophytes can colonize host but do not cause strong hypersensitive reactions to plant.9 2

ACS Paragon Plus Environment

Page 2 of 43

Page 3 of 43

Journal of Agricultural and Food Chemistry

39

The enormous potential of endophytes especially endophytic fungi affecting the accumulation

40

of SMs has received much attention as an integral part of their host’s biology, have been

41

proposed to be akin to host organs in and of themselves.10 Though limited data is available to

42

date, at least three possible ways are broadly accepted that highlight the effects of endophytes

43

on host physiology and metabolism: (ⅰ) secretion of versatile chemical entities including plant

44

hormones such as GA, IAA and ethylene,11 antibiotics or other products,12 which promote

45

plant metabolite accumulation, alter plant metabolic status or help plants respond to other

46

biotic and abiotic stresses. (ⅱ) by “induced effect” or “lateral gene transfer” altering host gene

47

expression and metabolism. Upon recognition of elicitors from endophyte, plants produce

48

reactive oxygen species (ROS) including H2O2 and O2-; resulting in ion flow of Cl−/Ca2+ and

49

K+/H+, signal net cross-talk such as cAMP, NO and SA; transcription factor activation and

50

regulation of gene expression culminating in the production of defense products.13 (ⅲ)

51

bio-conversion affecting plant chemical profiles. As shown for the endophytic

52

Paraconiothyrium variabile which specifically bio-transforms glycosylated flavonoids to

53

aglycone or other chemicals, which in turn alters the host Cephalotaxus harringtonia

54

metabolomic composition.14

55

Host plants on the other hand present selective pressure on endophytic fungal community

56

and shape the assembly of fungal endophytes on any given host spp.15 A plant of certain

57

genotype and age located in an environment/ecosystem harbors its own special endophyte

58

composition that may govern aspects of bi-directional SM accumulation and cross-talk.16 For

59

example, it is reported that comparison of Fusarium-infected spring wheat cultivars and

60

winter wheat cultivars showed altered accumulation of SMs (benzoxazinoids, phenolic acids, 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

61

carotenoids, and flavonoids), whose levels exhibited varying selective pressures on the fungal

62

pathogen Fusarium at five time points. Furthermore, in addition to certain other wheat SMs

63

the flavonoids, homoorientin and orientin were identified as potential inhibitors of the

64

trichothecene mycotoxin deoxynivalenol (DON) produced under Fusarium infection.17

65

The interaction between endophytic fungi and host plant is complex. The many challenges

66

arise from: (ⅰ) the enormous number of endophytic fungi in association with the plant

67

ecosystem; (ⅱ) complex interactions of host plants with numerous individual fungal species at

68

the same time; (ⅲ) the interactions occurring among endophytes themselves; and (ⅳ) the

69

elucidation of interactions that are authentic i.e., reproducibility under standard laboratory in

70

vitro conditions.9 However, objectively elucidating the correlation between the special

71

chemical composition and the corresponding endophytic fungus is a key first step in scientific

72

exploration and application of fungal endophytes in providing the much needed and sort after

73

resources of drugs in medicinal use against various forms of diseases.8,17,18 Recently, it is

74

reported that Humicola fuscoatra and Aspergillus niger significantly (p≤0.05) relate to host

75

development under salinity stress.19 The other notable example is that of the health hazard

76

caused by Lolium perenne an important pasture and forage plant to livestock`s which was

77

linked to the accumulation of alkaloids, lolitrem B and peramine which were found to

78

significantly positively correlate with endophytic wild type E (WT) strain Neotyphodium lolii

79

infection.20 Based on this find, N. lolii strains like AR1 (EAR1) and AR37 (EAR37) were

80

modified and inoculated on L. perenne, leading to the reduction of hazardous effects and

81

improved grass productivity.21

82

As an annual, herbaceous and succulent plant, C. songaricum must be collected before or 4

ACS Paragon Plus Environment

Page 4 of 43

Page 5 of 43

Journal of Agricultural and Food Chemistry

83

just at the developmental stage of sprouting to above ground emergence, which possesses the

84

highest pharmacological activity or else the plant would become hollow with low

85

pharmacological activity according to traditional experience and Chinese pharmacopoeia.22

86

We aimed at addressing the key question of, do endophytic fungi participate in the

87

physiological and developmental process of C. songaricum? To answer this question, we

88

asked: (1) what is the composition of endophytic fungi in C. songaricum? (2) how do

89

endophytic fungal communities evolve across its different developmental stages? (3) how do

90

the accumulations of main active ingredients change under endophytic associations? (4) are

91

there any correlations between active ingredients and the endophytic fungi? (5) If so, what are

92

these correlations and what do they imply? In part all or some of these cruces form the main

93

objective of this study. Our study shows a clear correlation between active ingredient

94

accumulation and endophyte assembly and presents new opportunities to better understand

95

plant fungal symbiosis illustrating the influential factors for improving medicinal quality of C.

96

songaricum.

97

MATERIALS AND METHODS

98

Plant Material Collection and Sample Preparation. C. songaricum developmental growth

99

was divided into five stages: Tubercle (T), Sprouting (S), Unearthing (U), Maturing (M) and

100

Atrophy (A) (Cui et al., 2018). The plants located around area of Xilin-gaole town (39°05′ 45′′

101

N, 105°23′27′′ E, ≈1133.98 m), Alashan League, Inner Mongolia of China, were collected in

102

November, March, May, June and September, respectively during the calendar year 2015 to

103

2016. The fresh rhizomes were collected and stored at 4℃, and transported to the laboratory

104

in < 8 hours. Collected sample materials were cleaned and surface-sterilized following 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

105

methods described earlier,18 then they were divided into two parts: one part of the collected

106

material was dried in the shade, and used for chemical analyses; while the other part was used

107

for DNA extraction and high-throughput sequencing of endophytic fungi.

108

Extraction and Analysis of Active Ingredients. The dried samples of C. songaricum were

109

powdered and passed through a 100 micron mesh sieve. A total of seven active ingredients in

110

C. songaricum were evaluated, and all standards and/or reference chemicals including

111

D-(+)-glucose (No.110833), rutin (No.100080), gallic acid (No. 110831), protocatechuic acid

112

(No.110809), catechins (No.110877) and ursolic acid (No.110742) were purchased from

113

National Institutes for Food and Drug Control, China. Total sugars, flavonoid and tannin

114

content were analyzed using UV spectrophotometry UV2450 (Shimadu, Japan), and the other

115

four ingredients (protocatechuic acid, catechins, gallic acid and ursolic acid) were determined

116

by HPLC 1200 LC system (Agilent Technologies, USA) (Supplement Figure S1). The

117

equations for linear regression, calibration coefficient, linear range, wavelength and retention

118

time (tR, in HPLC) are shown in Table 1. All the experimental treatments were done in

119

triplicate measurements.

120

Total sugars including oligosaccharides and polysaccharides were determined by modified

121

phenol-sulfate method.23 Firstly, lipid and protein content were removed with Soxhlet

122

circulation at 85.5℃ for 50 min and later the residue was dried and extracted with 80%

123

ethanol using circulation reflux at 80℃ for 2 h. The residue was dissolved with 25 mL

124

ultra-pure water by ultrasonication (350 w, 35 Hz) for 30 minutes. The mixture was

125

centrifuged at 5000 rpm for 5 min, to the supernatant 100 mL of 95% ethanol was added,

126

mixed well and keep static at 4 ℃ overnight, centrifuged at 5000 rpm for 5 min. The 6

ACS Paragon Plus Environment

Page 6 of 43

Page 7 of 43

Journal of Agricultural and Food Chemistry

127

pellet/residue was washed with 5 mL of 80% ethanol, followed by pure diethyl ether and

128

acetone, respectively dried at 50 ℃ under low pressure (0.09 MPa) using vacuum drying oven

129

DZF-6020 (Ruijia Scientific co. LTD, China). Lastly, the residue was weighed and dissolved

130

in 50 mL of ultra-pure water.

131

Flavonoids were determined by colorimetric method.24 Briefly, 1 g powder was placed in a

132

25 mL flask, sonicated (350 w, 35 Hz) with 25 mL of 50% ethanol for 35 min three times and

133

the extraction combined into a 100 mL volumetric flask, then adjusted to final volume of 100

134

mL with 50% ethanol and used as sample for determination of flavonoid content. The

135

calibration curve was generated based on ten concentrations of reference rutin from 12.50

136

μg/mL to 125.00 μg/mL in 50% ethanol, with increments of 12.50 μg/mL.

137

Analysis of tannins was done by the method described in Chinese Pharmacopoeia.22 A

138

powdered sample (1 g) of C. songaricum was placed in brown flask, ultra-pure water to a 25

139

mL final volume was added and incubated overnight. Next day, ultrasonic extraction was

140

done for 30 min, and then filtered through Whatman® qualitative filter paper, Grade 1. The

141

above process was repeated three times, the extracts combined, and the volume made up to

142

100 mL with ultra-pure water. Two mL of sample was mixed with 1 mL

143

phosphomolybdotungstic acid and 10 mL H2O, and the total phenol content was determined

144

by standard curve using UV spectrophotometry UV2450. For determination of unreacted

145

phenol, casein (0.6 g) was dissolved in 25 mL of sample in a 100 mL conical flask, sealed and

146

shaken at 30 ℃ for 1 h, then 2 mL of cold sample filtrate mixed with 1 mL

147

phosphomolybdotungstic acid and 7 mL of H2O. Finally, 29% Na2CO3 was added and the

148

final volume made up to 25 mL with ultra-pure water. The total tannin content (TC) was 7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

149

calculated using the formula: TC (mg) = total phenol (mg) – unreacted phenol (mg).

150

The contents of protocatechuic acid, catechins and gallic acid were simultaneously

151

determined by UV measurements at 280 nm by RP-HPLC-DAD.25 One gram of the sample

152

was weighed accurately and placed in a 50 mL flask and ultrasonicated with 40 mL of 50%

153

aqueous methanol (v/v) for 40 min. This was repeated twice, the combined extract volumes

154

were adjusted to a final volume of 50 mL, the solution was then passed through 0.45 μm

155

nylon membrane filter (Merck Millipore, USA). The samples were then analyzed using

156

Agilent 1200 LC system equipped with diode array detector (DAD), a quaternary pump

157

equipped with Eclipse XDB-C18 column (5.0 μm, 250×4.6 mm, Agilent, USA). The column

158

oven temperature was fixed at 25℃. A mobile phase consisting of acetonitrile (A)

159

(ThermoFisher Scientific, USA) and 0.1% glacial acetic acid:water buffer (B), was applied for

160

the gradient elution: 8A/92B to 100A/0B in 25 min with a flow rate of 1.0 mL /min with an

161

injection volume of 10 μL.

162

Ursolic acid content was determined using HPLC method.26 One gram of C. songaricum

163

sample powder was placed in a 50 mL volumetric flask, sonicated with 50 mL chloroform for

164

40 min, filtered, and impregnated twice with petroleum ether after evaporation to complete

165

dryness. After decanting the petroleum ether, the residue was taken up into a 10 mL

166

volumetric flask and the final volume adjusted with chloroform. The dissolved sample was

167

then filtered with 0.45 μm nylon membrane and used directly for HPLC analysis as described

168

above. The mobile phase consisted of acetonitrile (A) and 1% glacial acetic acid:water buffer

169

(B): 88A/12B in 15 min.

170

DNA Extraction and Amplicon Preparation. To characterize the endophytic fungal 8

ACS Paragon Plus Environment

Page 8 of 43

Page 9 of 43

Journal of Agricultural and Food Chemistry

171

community in C. songaricum in different developmental stages, genomic DNA was extracted

172

from surface-sterilized rhizome of C. songaricum using E.Z.N.A. Plant DNA Mini Kit

173

(Omega Bio-Tek, Doraville, GA, USA). PCR amplicons were then purified using PowerClean

174

Pro Clean-Up DNA Kit (MO BIO Laboratories, Inc., USA). The internal transcribed spacer 2

175

(ITS2) region was selected for PCR amplification with primer combinations of ITS3f

176

(GCATCGATGAAGAACGCAGC) and ITS4 (TCCTCCGCTTATTGATATGC) complete

177

with Illumina adapters and unique 7-bp barcode sequences that allowed for identification of

178

individual samples in multiplexed runs.18 PCR was performed in a 25 μL reaction system

179

containing of 1 μL of each primer (10 μM), 2 μL of template (50-100 ng/μL), 0.2 μL of

180

TaqDNA Polymerase (5 U/μL), 1 μL of dNTP (10 mM), 2.5 μL of 10×buffer, and 17.3 μL of

181

ddH2O. PCR amplifications were performed using: 3 min of initial denaturation at 94℃,

182

followed by 40 cycles of: 94℃ for 30 s, 47℃ for 30 s, and 72℃ for 45 s; and a final

183

extension of 7 min at 72℃. PCR products from triplicate runs were pooled to make one

184

sample type to avoid PCR bias. After purification using PureLink Quick gel extraction Kit

185

(Thermo Fisher Scientific, USA), DNA concentration was determined with Nanodrop

186

ND-1000 with software ver.3.3 (Thermo Scientific, USA). Finally, all bar-coded samples

187

were then pooled together in an equimolar ratio and were sent to high-throughput amplicon

188

sequencing using the Illumina Miseq platform (USA) at Sangon Biotech (Shanghai) Co., Ltd.

189

(Shanghai, China).

190

Next-generation Sequencing and Bioinformatic Analysis. Following sequencing, forward

191

sequence reads were de-multiplexed and sorted according to their barcodes. Similarly, the

192

reverse primer barcode sequences and low-quality regions of each sequenced sample were 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

193

removed through CUTADAPT and TRIMMOMATIC. Bioinformatics and statistical analyses

194

of pair-end sequences were done as described recently in Cui et al.18 Sequences obtained from

195

next-generation sequencing were deposited in the Sequence Read Archive (SRA) database in

196

NCBI (https://submit.ncbi.nlm.nih.gov) (SRA Accession: SRP128486).

197

The fungal OTU was picked with open picking method through the centroid based

198

algorithm.27 The similarity between samples was tested to determine the OTU occurrences

199

and statistical analysis on the diversity and evenness of the sample species was performed.

200

The alpha-diversity indices including ACE and Chao-1 richness index predicting the total

201

numbers of OTUs in each sample, Shannon and Simpson index were used to investigate the

202

biodiversity based on the OTU richness. The similarity criteria that were used was: species

203

98%, genus 94%, family 90%, order 85%, class 80%, and phylum 75%.28 The difference of

204

endophytic fungal communities based on OTUs in different developmental stages of C.

205

songaricum was evaluated by ANOVA statistics, multi-dimension and Venn Diagram

206

package in R software. Furthermore, the difference among all developmental groups was

207

analyzed through beta-diversity including non_multi_dimension_analysis (NMDS), Principal

208

Co_ordinates (PCoA), Principal Component Analysis (PCA), Unweighted pair group method

209

with arithmetic mean (UPGMA) and Unweighted unifrac. The phylogenetic relationship in

210

representative OTUs was analyzed by Hierarchical Clustering dendrogram, Bray_Crutis_tree

211

and “OTU_co_network”, based on R software. Test of the null hypotheses of no differences

212

among the microbial community compositions were examined using permutation multivariate

213

analysis of variance (PERMANOVA)18 and the p-value was recalculated based on the

214

Bonferroni significance. 10

ACS Paragon Plus Environment

Page 10 of 43

Page 11 of 43

Journal of Agricultural and Food Chemistry

215

Dynamic Analysis of Endophytic Fungi. The assembly and dynamics of endophytic fungal

216

community was characterized in C. songaricum across T, S, U, M and A developmental

217

stages corresponding to their individual sampling dates in November, March, May, June and

218

September, respectively encompassing one developmental cycle. All samples analyzed in the

219

five developmental stages were replicated three times.

220

Analysis of Correlation. Unitary and multiple linear regression analyses was performed to

221

evaluate the correlation between endophytic fungal richness and the accumulation of active

222

components in C. songaricum across T, S, U, M and A developmental stages, which included

223

simple correlation (Pearson correlation coefficient) and multi-correlation coefficient. As for

224

samples not conforming to a normal distribution, non-parametric test of Wilcoxon Signed

225

Rank was used to evaluate the variation in fungal richness and the contents of seven active

226

components using IBM SPSS Statistics 19.0 (Chigago, USA) with the parameters of

227

zero-order correlation coefficient (Pearson) and two-tailed test. In the matrices, correlation

228

coefficients are marked in red and green for positive and negative correlations, respectively.

229

While, light blue and dark blue colors are used to mark significant difference (p