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