Silencing of TT8 and HB12 Affects Nutritional ... - ACS Publications

the alias() function in STATS package was used to identify predictors that were. 178 ..... profiles were negatively related to SCP and EE (r < -0.48, ...
1 downloads 5 Views 624KB Size
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