Subscriber access provided by UNIVERSITY OF MICHIGAN LIBRARY
Biotechnology and Biological Transformations
Silencing of TT8 and HB12 Affects Nutritional Profiles and in vitro gas production Relating to Molecular Structures of Alfalfa (Medicago sativa) Plants Yaogeng Lei, Abdelali Hannoufa, Luciana Louzada Prates, Haitao Shi, Yuxi Wang, Bill Biligetu, David Christensen, and Peiqiang Yu J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01573 • 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 32
Journal of Agricultural and Food Chemistry
1
Silencing of TT8 and HB12 Affects Nutritional Profiles and
2
in vitro gas production Relating to Molecular Structures of
3
Alfalfa (Medicago sativa) Plants
4 5
Yaogeng Lei1, Abdelali Hannoufa2, Luciana Louzada Prates1, Haitao Shi1, Yuxi
6
Wang3, Bill Biligetu3, David Christensen1, Peiqiang Yu1*
7 8
1
Department of Animal and Poultry Science, University of Saskatchewan, 51 Campus
9
Drive, Saskatoon, SK S7N5A8, Canada
10
2
London Research and Development Centre, Agriculture and Argi-Food Canada,
11
1391 Sandford Street, London, ON N5V 4T3, Canada
12
3
13
Lethbridge Research and Development Centre, Agriculture and Argi-Food Canada, 5403 1st Avenue South, Lethbridge, AB T1J 4B1, Canada
14
4
Department of Plant Sciences, University of Saskatchewan, 51 Campus Drive,
15
Saskatoon, SKS7N5A8, Canada
16 17
Running Head: Gene Transformations of TT8 and HB12 RNAi Affected Nutrient
18
Profiles and gas production correlated with molecular structure of alfalfa
19 20 21
Corresponding author:
22
Peiqiang Yu
23
College of Agriculture and Bioresources, University of Saskatchewan
24
51 Campus Drive, Saskatoon, SK, S7N 5A8 Canada.
25
Tel:
1-306-966-4128
26
Fax:
1-306-966-4151
27
Email:
[email protected] 1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
28
Abstract
29
The objective of this study was to investigate the effects of silencing TT8 and HB12
30
genes on the nutritive profiles and in vitro gas production of alfalfa relating to spectral
31
molecular structures of alfalfa. TT8-silenced (TT8i, n=5) and HB12-silenced (HB12i,
32
n=11) alfalfa were generated by RNA interference (RNAi), and were grown with non-
33
transgenic control (NT, n=4) in a greenhouse. Alfalfa plants were harvested at early-
34
to-mid vegetative stage. Samples were analyzed for chemical compositions, CNCPS
35
fractions and in vitro gas production. Correlations and regressions of chemical and
36
CNCPS profiles with molecular spectral structures were also determined. Results
37
showed that transformed alfalfa had higher digestible fiber and lower crude protein
38
with higher proportion of PC. HB12 RNAi had lower gas production compared to
39
others. Some chemical, CNCPS and gas production profiles were closely correlated
40
with spectral structures and could be well predicted from spectral parameters. In
41
conclusion, RNAi silencing of TT8 and HB12 in alfalfa altered the chemical profiles
42
and CNCPS fractions of alfalfa, and such alterations were closely correlated with
43
inherent spectral structures of alfalfa.
44 45
Keywords: Genetic Transformation, Alfalfa (Medicago sativa), CNCPS Fractions,
46
Spectral Structures, Correlation and Regression, ATR-FTIR, HB12, TT8
47 48
2
ACS Paragon Plus Environment
Page 2 of 32
Page 3 of 32
Journal of Agricultural and Food Chemistry
49
Introduction
50
Significant advances have been achieved in genetic engineering-based breeding of
51
agricultural crops in recent years1. Compared to conventional breeding, genetic
52
engineering is more efficient, more accurate, and less time-consuming. Taking
53
advantage of the advanced molecular techniques and up-to-date bioinformatics tools,
54
genetic engineering is extremely powerful in manipulating important agronomic traits
55
of plants to meet human needs2–6. Apart from developing new genetic technologies,
56
research efforts have also focused on identification and characterization of genes that
57
have the potential to improve forage quality. Transparent Testa 8 (TT8) and Homo
58
box 12 (HB12) are two transcription factors in plants. TT8 encodes a bHLH (basic
59
helix-loop-helix) protein which interacts with myeloblastosis oncogene (MYB) and
60
WD40-repeat proteins to form a ternary complex7,8. This complex controls
61
transcription of late biosynthesis genes in the phenylpropanoid pathway, thereby
62
regulating biosynthesis of anthocyanins and proanthocyanidins9,10. The HB12, on the
63
other hand, belongs to the homeodomain-leucine zipper class I (HD-Zip I) family and
64
is more related to drought resistance in plants. The expression of HB12 increased
65
when plants were treated with salt or abscisic acid11,12. Observations from Brassica
66
napus showed that expression levels of TT8 and HB12 are positively associated with
67
lignin content13, indicating their potential use as targets towards lignin reduction in
68
forage crops.
69
Although alfalfa (Medicago sativa L.) contains favorable nutrient content and has
70
good palatability, the relatively high lignin content and rapid degradation of protein
71
are two drawbacks that limit the utilization of alfalfa forage1. The rapid degradation of
72
protein in the rumen can lead to rumen bloat in grazing animals, causing huge
73
economic losses to producers14. In contrast, high lignin content of alfalfa hinders the 3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
74
degradation of carbohydrate and other nutrients1. Altogether, these two drawbacks
75
make alfalfa more imbalanced in terms of ruminal nutrient synchronization.
76
Therefore, we transformed alfalfa with TT8 and HB12 RNAi constructs to explore the
77
effect of silencing these two genes on chemical profiles, CNCPS fractions and in vitro
78
gas productions of alfalfa. Meanwhile, the Attenuated Total Reflectance Fourier
79
Transform Infrared (ATR-FTIR) spectroscopy was used to explore relationships
80
between nutritional and fermentation profiles with structural spectral structures. We
81
hypothesized that genetic transformations will alter the chemical and nutritional
82
profiles of alfalfa, and that such alterations will affect the in vitro fermentation of
83
alfalfa and will be closely related to spectral structures.
84 85
Materials and Methods
86
Alfalfa Samples
87
Both transformed and non-transformed (NT) alfalfa samples were obtained from
88
Agriculture Agri-Food Canada (AAFC) Saskatoon Research Center. Details about
89
designing and making RNAi constructs of TT8 and HB12, transformation of alfalfa,
90
growth conditions of alfalfa plants, and harvests method were described in our
91
previous publications13,15. Briefly, total RNA was extracted from alfalfa for cDNA
92
synthesis, and making RNAi constructs for TT8 and HB12 by using Gateway systems
93
and pHellsgate12 vector. Afterwards, RNAi constructs were used to transform alfalfa
94
explants via Agrobacterium tumefaciens according to Aung et al.16. All plants were
95
grown in a greenhouse under normal conditions and harvested at early-to-mid
96
vegetative stage. Alfalfa plants were then freeze-dried and stored in individual bags
97
with each bag containing samples from one pot in the greenhouse. There were 11
98
HB12 RNAi (HB12i), 5 TT8 RNAi (TT8i) and 4 NT plants. Samples were ground 4
ACS Paragon Plus Environment
Page 4 of 32
Page 5 of 32
Journal of Agricultural and Food Chemistry
99
through a 1 mm sieve (Retsch SM-3000, Brinkmann Instruments, ON, Canada) for
100
chemical analysis, and through a 0.02 mm sieve for spectra collection, respectively.
101
Chemical Analysis
102
Contents of dry matter (DM), ash, crude protein (CP) and ether extract (EE) in
103
alfalfa samples were analyzed according to the AOAC methods17. Neutral detergent
104
fiber (NDF), acid detergent fiber (ADF) and acid detergent lignin (ADL) were
105
analyzed using ANKOM A200 Filter Bag Technique following the procedure
106
provided by ANKOM Technology18. Neutral detergent insoluble crude protein
107
(NDICP), acid detergent insoluble crude protein (ADICP) and non-protein nitrogen
108
(NPN) were determined according to the procedure described by Licitra et al.19.
109
Soluble crude protein (SCP) was analyzed according to the method of Roe et al.20.
110
Total starch was determined using the Megazyme method21, and sugar content was
111
analyzed according to Dubois et al22. Moreover, total carbohydrate content (CHO)
112
and non-fiber carbohydrate (NFC) were calculated according to Higgs et al.23.
113
CNCPS Fractions
114
The updated version (v6.5) of Cornell Net Carbohydrate and Protein System
115
(CNCPS) was used to calculate CNCPS carbohydrate and protein fractions23,24. The
116
CNCPS divides carbohydrate into eight fractions, CA1, CA2, CA3, CA4, CB1, CB2,
117
CB3, and CC; divides protein into five fractions, PA1, PA2, PB1, PB2, and PC. As
118
alfalfa hay rarely contains organic acids13 and ammonia, we calculated CA4 (sugar),
119
CB1 (starch), CB2 (soluble fiber), CB3 (digestible fiber), CC (indigestible fiber) in
120
carbohydrate fractions and PA2 (soluble true protein), PB1 (insoluble true protein),
121
PB2 (fiber-bound protein), PC (indigestible protein) in protein fractions in the current
122
study. The CNCPS carbohydrate and protein fractions were calculated on total
123
carbohydrate (CHO) and CP basis, respectively. 5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
124
ATR-FTIR spectrum and univariate analysis
125
Alfalfa samples that were finely-ground (0.02 mm) were used for FTIR spectra
126
collection with JASCO FT/IR-4200 with ATR (JASCO Corp., Tokyo, Japan) at
127
University of Saskatchewan, Canada. Five spectra of each sample were collected at
128
the mid-IR region (ca. 4000-700 cm-1) with 128 scans at a resolution of 4 cm-1.
129
Afterwards, FTIR spectra were normalized and saved as “csv” file along with their
130
second derivatives by using OMNIC 7.3 software (Spectra Tech, Madison, WI, USA).
131
Spectral parameters of peak heights and peak areas were calculated with Excel®
132
macro in carbohydrate region (ca. 1484-941 cm-1), amide (ca. 1710-1484 cm-1) and
133
lipid-related region (ca. 3000-2761 cm-1). In carbohydrate region, there were four
134
peaks (TC1, ca. 1025 cm-1; TC2, ca. 1074 cm-1; TC3, ca. 1104 cm-1; TC4, ca. 1149
135
cm-1) in total carbohydrate region (TC ca. 1178-941 cm-1); four peaks (STC1, ca.
136
1317 cm-1; STC2, ca. 1370 cm-1; STC3 ca. 1397 cm-1; STC4, ca. 1453 cm-1) in
137
structural carbohydrate region (STC, ca. 1484-1178 cm-1); and one peak (CEC, ca.
138
1237 cm-1) in cellulosic compounds region (CEC, ca. 1283-1178 cm-1). In amide
139
region, peaks of amide I (ca. 1649 cm-1) and amide II (ca. 1540 cm-1), and protein
140
secondary structure of α-helix (ca. 1653 cm-1) and β-sheet (ca. 1629 cm-1) were
141
measured. In lipid-related region, asymmetric CH3 (AsCH3, ca. 2955 cm-1),
142
asymmetric CH2 (AsCH2, ca. 2920 cm-1), symmetric CH3 (SyCH3, ca. 2872 cm-1)
143
and symmetric CH2 (SyCH2, ca. 2850 cm-1) were measured. Peak areas of total
144
carbohydrate (TCA), structural carbohydrate (STCA), cellulosic compounds (CECA),
145
total amide (AA), amide I (AIA), amide II (AIIA) and (a)symmetric CH2/CH3
146
(ASCCA) were also determined.
147
In vitro gas production fermentation
148
In vitro gas production fermentation was conducted at Lethbridge Research and 6
ACS Paragon Plus Environment
Page 6 of 32
Page 7 of 32
Journal of Agricultural and Food Chemistry
149
Development Center of AAFC (Alberta, Canada), according to the description of
150
Wang et al25 with minor modifications. About 0.3 g of ground samples was placed in
151
ANKOM F57 bags and then incubated with rumen fluid and buffer solution in a 125
152
mL serum vial. Rumen fluid was collected from three angus heifers at the Research
153
Center farm after morning feeding. Samples were incubated for 48h in the incubator
154
at 39 Celsius degrees with two experimental runs. Two replicates of each alfalfa
155
sample and two blanks were withdrawn from the incubator after 2, 4, 8, 12, 24 and 48
156
h of incubation. Gas production was measured at 2, 4, 8, 12, 24, 48 h for all samples
157
available by using a water replacement equipment, according to Wang et al.25 .
158
Statistical Analyses
159
Gas production kinetics model
160
Gas production kinetics were calculated with the non-linear model of ܲ = ܽ(1 −
161
݁ ି(௧ିሻ ሻ according to the description of Jonker
162
production at the time of t, a is asymptotic gas production, c is the fractional rate
163
(%/h), t is the incubation hours, and lag is the initial delay of gas production onset.
164
Average production (AP) was calculated as ( ÷ ܿ × ܽ = ܲܣ2 × (݈݊2 + ܿ × ݈ܽ݃ሻሻ ,
165
according to Jonker 26.
166
Correlations between nutritional and gas profiles with spectral parameters
167
26
. In this model, P is gas
Correlations and regressions between nutritional profiles and spectral parameters
168
were performed with R software
169
were obtained with the rcorr() function in HMISC package. Prior to correlations
170
analysis, normality of each variable was tested with the shapiro-test() function in
171
STATS package. Correlations involved non-normally distributed variables were
172
performed with “Spearman method”, while others were performed with “Pearson”
173
method.
27
. Correlation coefficients and their significances
7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
174
Multilinear regression of predicting nutritional and gas production profiles from
175
spectral parameters
176
Multiple linear regressions of predicting chemical and CNCPS profiles from
177
spectral parameters were conducted with the lm() function in STATS package. Then,
178
the alias() function in STATS package was used to identify predictors that were
179
completely linearly dependent on other predictors. The complete aliased predictors
180
were then deleted, and the vif() function from CAR package was used to check the
181
variance inflation factor (VIF) of each predictor. The predictor that had the highest
182
VIF, which was greater than 10, was removed from the model. The linear modeling
183
was re-conducted until the VIF values of all predictors were lower than 10.
184
Afterwards, the step() function with “direction = both” was used to select the linear
185
model that had the lowest AIC value. Then, predictors that were insignificant at
186
α=0.05 were removed from the model.
187
Statistical analysis
188
PROC MIXED procedure in SAS 9.4 (SAS Institute, Inc., Cary, NC, USA) was
189
used for analyzing chemical composition and gas production and kinetics data in this
190
study. The model for chemical data was ܻ = μ + trt + ߝ , where Yij is the
191
dependent variable, µ is the population mean, trti is the treatment effect, and εij is the
192
random error. The model for gas production and kinetics was ܻ = μ + trt +
193
ܾ݈݇ܿ + ߝ , where blockj is the fermentation run and other parameters are the same
194
as the previous model. Prior to variance analysis, observations with Studentized
195
residual greater than 2.5 were considered as outliers and removed from the dataset.
196
Differences between NT and transformed alfalfa were determined with contract
197
statement. The Tukey-Kramer method was used in multi-comparison after variance
198
analysis and a SAS macro called “pdmix800”28 was used to denote the letter for each 8
ACS Paragon Plus Environment
Page 8 of 32
Page 9 of 32
Journal of Agricultural and Food Chemistry
199
treatment mean at the significance level of 0.05. Normality of residual of each
200
variable was tested by using PROC UNIVARIATE in SAS 9.4 with normal and plot
201
options.
202 203
Results and Discussion
204
Effect of TT8 and HB12 silencing on Chemical Profiles
205
Chemical compositions of transformed and NT alfalfa are shown in Table 1.
206
Compared with NT control, transformed alfalfa had lower DM, ash, CP, EE, and
207
starch, but higher OM, CHO, NDF, ADF, ADL, NDICP, ADICP, and sugar (P 0.59, P < 0.01). Moreover, TC2, TC3 and TCA were negatively correlated with
358
CB1 (r < -0.45, P < 0.05). The CEC height was found to be negatively correlated with
359
PC fraction (r = -0.61, P < 0.01). In STC profiles, all STC peaks and STCA were
360
negatively correlated with CB1 and PB1 (r < -0.56, P < 0.01), while positively
361
correlated PC (r > 0.47, P < 0.05). Except for STC2, all other STC profiles were
362
positively correlated with CC fraction (r > 0.58, P < 0.01). In amide region, β-sheet,
363
AA and AIA were negatively correlated with CB1 and PB1 (r < -0.57, P < 0.01), but
364
positively correlated with CC and PC fraction (r > 0.57, P < 0.01). In lipid-related
365
ASCC region, SyCH2 and AsCH2 were negatively correlated with CB1 and PB1 (r
0.63, P < 0.01).
367
In our pilot study15, STC3 was positively correlated with CA4 and CB1 (r = 0.83),
368
while STC2 was negatively correlated with CB3 (r = -0.81). Moreover, TC1, TC3,
369
TC4 and TCA tended to be positively correlated with CA4 and CB1 in our pilot
370
study15 (r > 0.75, P < 0.1). Except for TC1 and TC4, which also tended to have
371
weakly positive correlations with CB1 in the present study, all other correlations
372
found in the pilot study were inconsistent with the present study. In CNCPS fractions,
373
CB1, CC, PB1, PA2 and PC are starch, indigestible fiber, insoluble true protein, fiber15
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
374
bound protein and indigestible protein, respectively23. In the present study, structural
375
carbohydrate spectral profiles had negative correlations with CB1, and positive
376
correlations with CC and PC, which were consistent with the CNCPS descriptions.
377
This implies that samples with higher STC spectral heights will be more resistant to
378
rumen degradation. As for amide profiles, β-sheet was found to be negatively related
379
to CB1 and PB1, which are easily degraded in the rumen. On the other hand, β-sheet
380
had positively relationship with PC, the indigestible protein fraction. It was reported
381
that a higher intensity of β-sheet in samples would lead to lower protein digestibility
382
because of the resistance of
383
digestion44, which is in line with our results.
384
Correlations between gas production profiles and spectral parameters
β-sheet to microbial attachment and enzymatic
385
Correlations between gas production profiles and spectral parameters were mostly
386
negative, except that CEC was positively correlated with gas production after 4h and
387
12h of fermentation, as well as asymptotic gas production (Table 6). The TC2 and
388
TC3 in TC region, all STC profiles (except for STC2), β-sheet, AA, and AIA in amide
389
region, SyCH2, AsCH2 and ASCCA in lipid-related region, were all negatively
390
correlated with gas production during fermentation, as well as asymptotic and average
391
gas production (P