Transcriptomic and Proteomic Research To ... - ACS Publications

Aug 10, 2016 - ABSTRACT: Mungbean (Vigna radiata (L.) Wilczek) is an important rotation legume crop for human nutrition in Asia. Bruchids (Callosobruc...
0 downloads 0 Views 1MB Size
Subscriber access provided by SUNY DOWNSTATE

Article

Transcriptomic and proteomic research to explore bruchid-resistant genes in mungbean isogenic lines Wu-Jui Lin, Chia-Yun Ko, Mao-Sen Liu, Chien-Yen Kuo, Dung-Chi Wu, ChienYu Chen, Roland Schafleitner, Long-Fang Oliver Chen, and Hsiao-Feng Lo J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.6b03015 • Publication Date (Web): 10 Aug 2016 Downloaded from http://pubs.acs.org on August 16, 2016

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 free 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 accessible to all readers and 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.

Journal of Agricultural and Food Chemistry 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 44

Journal of Agricultural and Food Chemistry

Transcriptomic and proteomic research to explore bruchid-resistant genes in mungbean isogenic lines

Wu-Jui Lin,†,‡,# Chia-Yun Ko,‡,# Mao-Sen Liu,‡ Chien-Yen Kuo,§ Dung-Chi Wu,§ Chien-Yu Chen,§ Roland Schafleitner,|| Long-Fang O. Chen,‡ and Hsiao-Feng Lo*,†



Department of Horticulture and Landscape Architecture, National Taiwan University, Taipei,

10617, Taiwan ‡

Institute of Plant and Microbial Biology, Academia Sinica, Taipei, 11529, Taiwan

§

Department of Bio-industrial Mechatronics Engineering, National Taiwan University, Taipei,

10617, Taiwan ||

AVRDC- World Vegetable Center, Shanhua, Tainan, 74151, Taiwan

#These authors contributed equally to this work *CORRESPONDING AUTHORS

Dr. Hsiao-Feng Lo Department of Horticulture and Landscape Architecture, National Taiwan University, Taipei, 10617, Taiwan (Tel: +886-2-3366-4871; Email: [email protected])

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

1

ABSTRACT

2

Mungbean (Vigna radiata (L.) Wilczek) is an important rotation legume crop for

3

human nutrition in Asia. Bruchids (Callosobruchus spp.) currently cause heavy damage as

4

pests of grain legumes during storage. We used omics-related technologies to study the

5

mechanisms of bruchid-resistance in seeds of the nearly isogenic lines VC1973A

6

(bruchid-susceptible) and VC6089A (bruchid-resistant). 399 differentially expressed genes

7

(DEGs) were identified between the two lines by transcriptome sequencing. Among these

8

DEGs, 251 DEGs exhibited high expression levels and 148 DEGs expressed low expression

9

levels in seeds of VC6089A. 45 differential proteins (DPs) were identified by isobaric tags

10

for relative and absolute quantification (iTRAQ). The 21 DPs had higher abundances in

11

VC6089A, and 24 DPs had higher abundances in VC1973A. According to transcriptome

12

and proteome data, only three DEGs/DPs, including resistant-specific protein (g39185),

13

gag/pol polyprotein (g34458), and aspartic proteinase (g5551) were identified and located

14

on chromosome 5, 1 and 7, respectively. Both g39185 and g34458 genes encode a protein

15

containing a BURP domain. In previous research on bruchid molecular markers, the g39185

16

gene located closely to the molecular markers of major bruchid-resistant locus may be a

17

bruchid-resistant gene.

18 19

KEYWORDS: mungbean, bruchid, transcriptome, proteome, iTRAQ 2

ACS Paragon Plus Environment

Page 2 of 44

Page 3 of 44

20

Journal of Agricultural and Food Chemistry

INTRODUCTION

21

Mungbean, Vigna radiata (L.) Wilczek, is an important rotation leguminous crop in

22

south and southeast Asia.1-4 For human nutrition, mungbean seeds constitute valuable protein

23

and carbohydrate sources, and sprouts are important vitamin and mineral sources.2 However,

24

two bruchids, the azuki bean weevil (Callosobruchus chinensis) and cowpea weevil (C.

25

maculatus), attack mungbean and cause severe losses during storage. To overcome this

26

challenge, bruchid-resistance has constituted an important breeding goal of mungbean.5 The wild mungbean accession TC1966 (Vigna radiata var. sublobata), which is

27 28

completely resistant to C. chinensis, C. maculates, C. phaseoli and Z. subfasciatus, was used

29

for bruchid-resistance breeding.6 In the first report on bruchid-resistance, a wild mungbean

30

was controlled with bruchid-resistance by a single dominant locus.7 Based on segregation

31

populations of TC1966 and a mungbean cultivar (NM92) with different levels of

32

bruchid-resistance and -susceptibility, Chen et al. first reported that the bruchid-resistant

33

genes (Br genes) were controlled by one major locus and two minor loci in wild mungbean.5,

34

8

35

random amplified polymorphic DNA (RAPD), and DMB-SSR 158 was mapped on linkage

36

group 9 by simple sequence repeat (SSR).5, 8 The two minor loci that were found to be tightly

37

linked with molecular markers mg7pgc325 and ma3pat361 through analysis by amplified

38

fragment length polymorphism (AFLP) were located at 114 and 132 cM on linkage group 7,

The major locus was tightly linked with molecular markers W02a4 through analysis by

3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

39 40

respectively.8 Resistance mechanisms of plants against insect pests could be associated with antixenosis

41

and antibiosis.9 Hence, bruchid-resistance in legumes relies on anti-nutritional compounds

42

and secondary metabolites which are toxic to bruchids.10 A 4-week feeding study, comparing

43

a commercial mungbean with a bruchid-resistant isogenic line, showed no negative effects on

44

growth or any pathological effects on mice.11 Mungbean seeds contain 58.2~61.8%

45

carbohydrates, 22.9~23.6% protein, and 1.2% oil.12 Consequently, mungbean seeds can

46

provide the major nutrients, starch and protein, for bruchid larvae. The specific alleles of the

47

arcelin-phytohaemagglutinin-α-amylase (APA) locus from wild bean (Pbaseolus vugaris)

48

provide some resistance to bruchids.13 An arcelin of common bean is a lectin-like protein

49

bound to carbohydrates in the intestinal epithelium of insects and causes an alteration of

50

insect gut structure.14, 15 The VrD1/VrCRP protein detected in the mungbean seed coat

51

inhibited the development of C. maculates to adults in artificial seeds.16-18 VrD1 protein is a

52

specific α-amylase inhibitor that inhibited α-amylase of insects, but not of animals.16, 19

53

Moreover, although two novel cyclopeptide alkaloids, vignatic acid A and B, were isolated

54

from mungbean,20 they were not the principal factors responsible for the bruchid-resistance.21

55

Thus, the mechanisms of bruchid-resistance in TC1966 are still not clearly understood.

56 57

The mungbean genome size (2n = 22 chromosomes) was estimated as 494-579 Mb by flow cytometry.22, 23 Recent mungbean genome-sequencing and gene annotation of 4

ACS Paragon Plus Environment

Page 4 of 44

Page 5 of 44

Journal of Agricultural and Food Chemistry

58

bruchid-resistant mungbean, RIL59,22 and bruchid-susceptible mungbean, VC1973A,24

59

provide an excellent foundation for omic-related analyses. To elucidate the mechanisms of

60

bruchid-resistance, we used Illumina RNA-seq technology and the iTRAQ method to

61

investigate the differentially expressed transcripts and proteins in two isogenic mungbean

62

lines, VC6089A and VC1973A. By combining these omic-related technologies, our results

63

reveal that three Br candidate genes/proteins are involved in bruchid-resistance mechanisms.

64

The results will be applicable in mungbean breeding and in insect-resistance research of other

65

legume crops.

66 67

MATERIALS AND METHODS

68

Plant Materials

69

Mungbean bruchid-susceptible line (VC1973A, NM92) and bruchid-resistant line

70

(VC6089A, TC1966, and recombinant inbred line 59 (RIL59)) were obtained from the World

71

Vegetable Center - Asian Vegetable Research and Development Center (AVRDC). VC6089A

72

is a nearly isogenic line (BC6F2) of VC1973A. It was derived from cross of VC1973A and a

73

wild bruchid-resistant mungbean TC1966 (Vigna radiate var. sublobata), and then

74

bruchid-resistant hybrids were back-crossed for six times to VC1973A. The genome of

75

VC6089A is only 1% different from VC1973A.17, 25 The contents of protein and soluble sugar

76

are shown in Table S1. RIL59, one of 200 F12 RIL, was generated from an inter-subspecific 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

77

cross between bruchid-susceptible mungbean variety NM92 and a bruchid-resistant accession

78

TC1966 (Vigna radiata var. sublobata).5,8,22 The bruchid-resistance analysis with 40 seeds

79

was performed following the Chen et al. method.5 Mungbean seeds with 0% damage were

80

defined as bruchid-resistant and those with more than 80% damage as bruchid-susceptible.

81

RNA Extraction

82

Total RNA was extracted following the Pine Tree method26 with slight modification.

83

Seeds were ground into powder with a mortar in a pestle filled with liquid nitrogen. The

84

extraction buffer containing 2% CTAB, 2% PVP K30, 100 mM Tris-HCl pH 8.0, 25 mM

85

EDTA, 2.0 M NaCl, 0.5 gL-1 spermidine, and 2% beta-mercaptoethanol (added prior to use)

86

was pre-heated at 65 °C in a water bath. One gram of seed powder was extracted with a 10

87

mL extraction buffer by vortexing in a 50 mL Falcon tube. Then, an equal volume of

88

chloroform : isoamyl alcohol (24:1) was added and mixed well by vortexing. After

89

centrifuging under 12 000 g at room temperature for 10 min, the supernatant was transferred

90

to a new 50 mL tube. A half volume of phenol and 1 volume of chloroform:isoamyl alcohol

91

(24:1) were added to the supernatant and mixed by vortexing. After centrifugation, the

92

supernatant was transferred to a new tube and mixed with 1/3 volume of 8 M LiCl. RNA was

93

precipitated at 4 °C overnight. RNA was harvested by centrifuging under 12 000 g at 4 °C for

94

30 min. This was followed by discarding the supernatant, air-drying for 5 to 10 min, and then

95

dissolution in 200 µL DEPC-H2O. Contaminated DNA was removed by the TURBO 6

ACS Paragon Plus Environment

Page 6 of 44

Page 7 of 44

Journal of Agricultural and Food Chemistry

96

DNA-free kit (Ambion) following the manufacturer’s instructions. The DNA-free RNA

97

sample was precipitated by adding 1/10 volume of sodium acetate, 2 volume of 100%

98

ethanol, and 15 µg/mL linear acrylamide (Ambion) at -70 °C overnight. After centrifugation

99

and discarding the supernatant, RNA was pelleted by 1 mL 75% ethanol, air-dried for a few

100

minutes, and then dissolved in 32 µL DEPC-H2O. The RNA quality was confirmed by a

101

Bioanalyzer RNA 6000 NanoChip (Agilent Technologies, Santa Clara, CA, U.S.A.) coupled

102

with an Agilent 2100 Bioanalyzer (Agilent Technologies) at the DNA Microarray Core

103

Laboratory, Institute of Plant and Microbial Biology (IPMB), Academia Sinica, Taiwan.

104

RNA Sequencing

105

For transcriptomic analysis, two biological repeats of total RNA from VC1973A and

106

VC6089A were performed. Four paired-end RNA libraries were constructed for sequencing

107

by the Illumina Hiseq 2000 platform. The RNA-seq data of VC1973A and VC6089A are

108

available on the NCBI Sequence Read Archive under accession of SRP070726. The RNA-seq

109

reads (1973A-seed-1: 19052484 reads, 2.88 Gb; 1973A-seed-2: 17591440 reads, 2.66 Gb;

110

6089A-seed-1: 24883022 reads, 2.76 Gb; 6089A-seed-2: 19099398 reads, 2.88Gb) of all

111

samples were trimmed for low quality bases and then individually aligned to the set of

112

annotated transcripts using BWA MEM.27 For each data set, a quantification of transcript

113

expression was performed by using eXpress28 to calculate the transcripts per million (TPM)

114

for each transcript. 7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

115

Computing Differentially Expressed Genes

116

DESeq was used for differential expression analysis by calculating the total read counts

117

of a gene in each sample. A transcript was denoted as differential expression genes (DEGs) if

118

Padj> 0.1, P> 0.0529 and the fold change (FC) (resistant/susceptible) was greater than 2 or less

119

than 0.5, respectively. Otherwise, a transcript was denoted as non-differentially expressed.

120

Protein Extraction

121

Total protein was extracted by TRIzol Reagent (Ambion) following the manufacturer’s

122

instructions. Seeds were ground into powder with a mortar in a pestle filled with liquid

123

nitrogen. One hundred milligrams of seed powder was extracted with 1 mL TRIzol Reagent

124

and 0.2 mL chloroform. After centrifugation with 1 200 g at 4 °C for 5 min, the proteins in

125

the organic phase were transferred to a new tube, 0.3 mL 100% ethanol was added, and then

126

centrifuged with 2 000 g at 4 °C for 5 min. The supernatant was precipitated with 1.5 ml

127

isopropanol. After centrifugation with 1 200 g at 4 °C for 10 min, the protein pellets were

128

washed for three times with 2 mL 0.3 M guanidine hydrochloride in 95% ethanol. The

129

protein pellets were harvested, air-dried, and dissolved in100~200 µL 9 M urea for usage.

130

Protease Digestion

131

Protein treatment, protease digestion, and labeling prior to LC-MS analysis were

132

performed as described previously30 with minor modifications. Protein concentration was

133

measured with a Pierce 660 nm Protein Assay Kit (Thermo Scientific). Proteins of 100 µg 8

ACS Paragon Plus Environment

Page 8 of 44

Page 9 of 44

Journal of Agricultural and Food Chemistry

134

extracted from each sample were reduced in 10 mM dithiothreitol (DTT) for 1 h at 37 °C.

135

Cysteines (Cys) were alkylated in 50 mM iodoacetamide at room temperature for 30 min in

136

the dark. Protein solutions were then diluted to 4 M urea with 50 mM Tris-Cl pH 8.5 and

137

digested with 250 units/ml benzonase (Sigma-Aldrich) at room temperature for 2 h, followed

138

by Lys-C (Wako, Japan) digestion [1:200 (w/w)] at room temperature for 4 h. Protein

139

solutions were further diluted to < 2 M urea with 50 mM Tris-Cl pH 8.0, and incubated with

140

2 µg of modified trypsin (w/w, 1:50, Promega) at 37 °C overnight. These protease-digested

141

solutions were acidified with 10% trifluoroacetic acid, desalted using an Oasis HLB cartridge

142

(Waters, U.S.A.), and then dried with SpeedVac.

143

Peptide Labeling with Isobaric Tags and SCX Fractionation

144

The dissolution of dried peptides in dissolution buffer and labeling with iTRAQ reagents

145

were performed according to the manufacturer’s instructions (Applied Biosystems). Digested

146

peptides from mungbean VC6089A were labeled with iTRAQ 113, and those from VC1973A

147

were labeled with iTRAQ 117 reagents. The labeling reactions with iTRAQ reagents were

148

incubated for 1 h at room temperature. Following the reactions, solutions from all of the

149

different iTRAQ labels were combined and further fractionated on a strong cation-exchange

150

(SCX, PolySulfoethyl A, 4.6×200 mm, 5 µm, 200 Å, PolyLC) high-performance liquid

151

chromatography (HPLC). SCX chromatography was performed with initial equilibrium

152

buffer A containing 10 mM KH2PO4, 25% acetonitrile (ACN), pH 2.65, and followed by a 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

153

0~15% buffer B (1 M KCl in buffer B, pH 2.65) gradient for 20 min, 15~30% buffer B

154

gradient for 10 min, 30~50% buffer B gradient for 5 min, 50~100% buffer B gradient for 1

155

min, and 100% buffer B for 5 min. The flow rate was 1 mL/min. Chromatography was

156

recorded with Abs at 214 nm UV light. Fractions (0.5 min/fraction) were collected, pooled,

157

and separated into 25 final fractions. Samples were desalted using an Oasis HLB cartridge

158

(Waters, U.S.A.) prior to LC-MS/MS analysis.

159

LC-MS/MS Analysis

160

Pooled and desalted fractions were re-dissolved in 0.1% formic acid and analyzed using

161

LC-MS/MS (a DionexUltiMate 3000 RSLCnano LC system coupled to a Q Exactive hybrid

162

quadrupole-Orbitrap mass spectrometer equipped with a nanospray Flex ion source, Thermo

163

Scientific). A C18 capillary column (Acclaim PepMap RSLC, 75 µm × 250 mm, Thermo

164

Scientific) was utilized to separate peptides with a 120 min linear gradient from 3% to 30%

165

solvent B (0.1% formic acid in acetonitrile, ACN) at a flow rate of 300 ml/min. The Q

166

Exactive MS was operated in the data-dependent mode, with the top 10 ions (charge states ≥2)

167

for MS/MS analysis following the MS survey scan for each acquisition cycle. The selected

168

ions were isolated in the quadrupole, and subsequently activated using higher-energy

169

collisional dissociation (HCD) and analyzed in an Orbitrap cell. The dynamic exclusion

170

duration of ion selection was 15 s. The MS was set as follows: m/z 350~1,600 range,

171

resolving power of 70,000, automatic gain control (AGC) target of 3×106, and maximum IT 10

ACS Paragon Plus Environment

Page 10 of 44

Page 11 of 44

Journal of Agricultural and Food Chemistry

172

of 30 ms. MS/MS was set as follows: resolving power of 17,500, AGC target of 1×105, and

173

maximum IT of 150 ms. HCD was set at a collision energy of 30% normalized collision

174

energy (NCE).

175

Protein Identification and Quantification

176

Peptide identification was performed using Proteome Discoverer software (v1.4,

177

Thermo Fisher Scientific) with SEQUEST and Mascot (v2.5, Matrix Sciences) search

178

engines. MS data were searched against the mungbean protein sequence database with 49,951

179

entries. Search conditions were set as full trypsin digestion, two maximum missed cleavage

180

sites, precursor mass tolerance of 10 ppm, fragment mass tolerance of 50 mmu, dynamic

181

modifications

182

carbamidomethyl (C), and iTRAQ8plex (N-terminus and K). Peptide spectrum matches

183

(PSM) were validated using the target decoy PSM validator algorithm, which automatically

184

conducted a decoy database search, and rescored PSMs using q-values and posterior error

185

probabilities (PEP). All PSMs were filtered with a q-value threshold of 0.01 (1% false

186

discovery rate), and proteins were filtered with a minimum of two distinct peptides identified

187

per protein. For comparative protein quantification, the ratios (113/117) of iTRAQ reporter

188

ion intensities in the MS/MS spectra of PSMs were utilized to calculate fold changes between

189

samples. Only unique peptides were used for protein quantification. All peptide ratios were

190

normalized by the median protein ratio.

of

oxidation

(M)

and

iTRAQ8plex

(Y),

11

ACS Paragon Plus Environment

static

modifications

of

Journal of Agricultural and Food Chemistry

191

Statistical Analysis

192

Proteins with significant changes in abundance between the two mungbean isogenic

193

lines were selected using a method described previously.30 The mean and S.D. from the log 2

194

ratios of the 1,271 proteins overlapping in two biological repeats were calculated. Next, 95%

195

confidence (Z score = 1.96) was used to select those proteins with distribution removed from

196

the main distribution. For down-regulated proteins, the confidence interval was 0.00318

197

(mean ratio of the 1,271 proteins) – 1.96 × 0.189418 (S.D.), corresponding to a protein ratio

198

of 0.774814. Similarly, for up-regulated proteins, the mean confidence interval was

199

calculated (mean ratio + 1.96 × S.D.) to a protein ratio of 1.296338. Protein ratios outside of

200

this range were defined as being significantly different at P = 0.05. The cutoff value for

201

down-regulated proteins was 0.77-fold, and for up-regulated proteins was 1.30-fold.

202

Combining Transcriptomic and Proteomic Data

203

The lists of 399 DEGs and 45 DPs were saved as DEG.csv and DP.csv files, respectively. Two

204

files were analyzed, and then a Venn diagram was created by Partek Genomics Suite software

205

(Partek, Inc.). The overlap area of the Venn diagram represents the DEGs/DPs that belong

206

equally to both transcriptomic and proteomic data.

207

Chromosome Location Analysis

208 209

The sequences of genes that encode proteins with differential abundances were blasted on the mungbean genome database 12

ACS Paragon Plus Environment

Page 12 of 44

Page 13 of 44

Journal of Agricultural and Food Chemistry

210

(http://plantgenomics.snu.ac.kr/mediawiki-1.21.3/index.php/Main_Page)24 to identify their

211

locations on the chromosome.

212

Predicted protein domain and glycosylation sites

213

The predicted protein domains were analyzed by NCBI BLAST

214

(http://blast.ncbi.nlm.nih.gov/Blast.cgi). The predicted N-glycosylation and O-glycosylation

215

sites were analyzed by NetNGlyc and NetOGlyc software in the ExPASy website

216

(http://www.expasy.org/). Prediction of transmembrane helices in proteins was analyzed by

217

TMHMM server (http://www.cbs.dtu.dk/services/TMHMM/).

218 219

RESULTS AND DISCUSSION

220

Transcriptomic Analyses of Bruchid-Resistant and -Susceptible Mungbeans

221

Total RNA in VC1973A (bruchid-susceptible line) and VC6089A (bruchid-resistant line)

222

were extracted for RNA sequencing (RNA-seq) by Illumina Hiseq 2000 to search for

223

potential Br genes. Approximately 23,000 transcripts (54.8% of all transcripts in the

224

mungbean genome) were detected and annotated by BWA MEM in two biological repeats of

225

mungbean lines. The transcripts between VC1973A and VC6089A were analyzed by

226

DESeq.29 The transcripts with Padj lower than 0.1, P value lower than 0.05, and fold change

227

(FC) higher than 2 or lower than 0.5 were selected and defined as differential expression

228

genes (DEGs). A total of 399 DEGs were identified, among which 251 DEGs showed 13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

229

up-expression and 148 DEGs showed down-expression in the seeds of bruchid-resistant line

230

(Table S2). The numbers of DEG exhibited a 1.73% difference in seed transcriptome

231

between VC1973A and VC6089A. Six DEGs (g11456, g4739, g34480, g21777, g27780, and

232

g39181) were only expressed in the resistant line and one DEG (g13480) in the susceptible

233

line. However, these seven DEGs showed lower expression levels (TPM 5.2 or < 0.35 in Table S2. Inf represented the genes expressed in VC6089A, but not in VC1973A.

b

35

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 36 of 44

Table 2. Differential proteins between seeds of mungbean isogenic lines VC6089A and VC1973A ID

Description

Accession

Coverage

E-value

Ident

Fold change

Location

Functional category

g34458.t1 Gag/pol polyprotein

AAQ82037.1

76%

0

59%

4.460

Vr01

Other

g5551.t1

Aspartic proteinase

AAB03843.2

99%

0

94%

4.324

Vr07

Enzyme

g39185.t1 Resistant-specific protein-1(4)

BAC22498.1

99%

0

83%

3.772

Vr05

Other

BAA77691.1

99%

8.00E-96

94%

2.388

Vr09

Defense-related protein

g39724.t1 Small ubiquitin-related modifier 1

XP_008456407.1

98%

7.00E-59

91%

2.051

scaffold_335

g8174.t1

XP_003553937.1

99%

0

89%

1.584

Vr01

Enzyme

g21583.t1 40S ribosomal protein S4-1

KHN38535.1

99%

5.00E-161

98%

1.557

Vr06

Protein synthesis and digestion

g8171.t1

Alpha-dioxygenase 1

XP_002279884.1

98%

0

70%

1.542

Vr01

Enzyme

g1322.t1

Stem-specific protein TSJT1-like

XP_003545925.1

99%

4.00E-160

86%

1.537

Vr10

Development

g30519.t1 40S ribosomal protein S7-2

KHN47668.1

99%

1.00E-123

92%

1.512

Vr01

Protein synthesis and digestion

g3584.t1

KHN16793.1

53%

6.00E-64

99%

1.415

Vr07

Protein synthesis and digestion

g10115.t1 Early nodulin-like protein 1

KHN37083.1

96%

1.00E-25

66%

1.398

Vr06

Other

g6021.t1

Signal peptide peptidase-like isoform X1

XP_003526473.1

99%

0

96%

1.392

Vr05

Enzyme

g8483.t1

PITH domain-containing protein

KHN16193.1

99%

1.00E-118

95%

1.368

Vr08

Other

g19167.t1 TPR superfamily protein

KEH24003.1

75%

2.00E-139

74%

1.357

Vr03

Other

g38165.t1 Isoflavone reductase

KHN41959.1

99%

4.00E-178

82%

1.345

Vr11

Enzyme

g17708.t1 Cu/Zn superoxide dismutase

ADZ72850.1

99%

8.00E-92

100%

1.338

Vr10

Enzyme

g42017.t1 Lipid transfer protein I

AAQ74627.1

99%

2.00E-63

100%

1.331

scaffold_91

g17385.t1 tRNA (mo5U34)-methyltransferase

KHG13272.1

94%

5.00E-133

68%

1.329

Vr01

Enzyme

XP_003548805.1

99%

0

97%

1.312

scaffold_7

Enzyme

XP_004508989.1

99%

3.00E-144

96%

1.310

Vr11

Other

KHN12858.1

91%

3.00E-118

94%

0.767

Vr10

Chaperone

g27390.t1 Cowpea pathogenesis-related protein 3 (CpPR3)

Alpha-dioxygenase 1-like

40S ribosomal protein S10

DNA repair

Defense-related protein

g21722.t1 Glyceraldehyde-3-phosphate dehydrogenase A, chloroplastic-like g1820.t1

Ras-related protein RABH1b-like

g21463.t1 Putative prefoldin subunit 3

36

ACS Paragon Plus Environment

Page 37 of 44

g2325.t1

Journal of Agricultural and Food Chemistry

Sulfite reductase [ferredoxin], chloroplastic

XP_003537728.1

91%

0

93%

0.764

Vr02

Enzyme

g32355.t1 Actin-related protein 4-like

XP_003521054.1

99%

0

96%

0.762

Vr09

Development

g18315.t1 Ubiquitin receptor RAD23d

XP_008230399.1

99%

3.00E-177

71%

0.761

Vr08

Protein synthesis and digestion

g2293.t1

XP_003607456.1

54%

0

89%

0.759

Vr02

Other

g17312.t1 60S ribosomal protein L7-4

KHN24635.1

58%

9.00E-120

85%

0.754

Vr05

Protein synthesis and digestion

g9314.t1

XP_003546884.1

99%

3.00E-132

79%

0.754

Vr05

Enzyme

0.742

Vr06

Other

Prohibitin 1-like protein

ATP synthase subunit O, mitochondrial-like

g10116.t1 Uncharacterized protein g31165.t1 Glutamate 5-kinase

XP_003613945.1

97%

0

89%

0.738

scaffold_100

Enzyme

g36464.t1 Nucleic acid-binding, OB-fold-like protein

XP_007044641.1

88%

2.00E-73

90%

0.731

scaffold_293

DNA repair

g176.t1

KHN43543.1

99%

2.00E-175

93%

0.726

Vr03

Protein synthesis and digestion

g19307.t1 Non-specific lipid-transfer protein A

KHN12008.1

90%

5.00E-23

47%

0.725

Vr06

Defense-related protein

g8144.t1

Formate dehydrogenase

ACZ74696.1

99%

0

96%

0.724

Vr01

Enzyme

g5544.t1

L-ascorbate oxidase homolog

XP_003549771.1

99%

0

91%

0.722

Vr07

Enzyme

XP_003547730.1

94%

1.00E-50

46%

0.720

scaffold_2330

KF033519.1

100%

2.00E-171

89%

0.702

Vr10

Other

g26926.t1 Nuclear cap-binding protein subunit 1-like

XP_006599000.1

99%

0

91%

0.695

Vr04

Gene expression and regulation

g27806.t1 Pantothenate kinase

AES75860.2

99%

0

91%

0.684

Vr01

Enzyme

g32713.t1 Probable aldo-keto reductase 1-like

XP_006577284.1

99%

0

88%

0.682

Vr10

Enzyme

g5877.t1

AES88942.2

98%

0

77%

0.666

Vr02

Other

XP_003531365.1

99%

0

92%

0.657

Vr01

Enzyme

OTU domain-containing protein 6B

g18113.t1 Germin-like protein-like g30742.t1 Phaseolus vulgaris clone BE5d2113

Development

glutamate-rich protein mRNA

DUF1680 domain protein

g33674.t1 Bifunctional aspartate aminotransferase and glutamate/aspartate-prephenate aminotransferase-like g6005.t1

Glycine-rich RNA-binding protein GRP1A

XP_008236997.1

48%

7.00E-45

84%

0.636

Vr02

Gene expression and regulation

g851.t1

Pathogenesis-related protein 10

AAX19889.1

99%

4.00E-104

100%

0.557

Vr07

Defense-related protein

XP_006575386.1

99%

6.00E-159

98%

0.556

Vr06

Protein synthesis and digestion

g19819.t1 60S ribosomal protein L10

37

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

521

TOC Graphic: Exploring bruchid-resistant candidate genes in mungbean

522

isogenic lines by omics technologies

523

38

ACS Paragon Plus Environment

Page 38 of 44

Page 39 of 44

Journal of Agricultural and Food Chemistry

Figure 1. The functional categories (A) and chromosome location (B) of 399 differential expressed genes (DEGs) analyzed in two isogenic lines. The up-expressed DEGs (dark gray) and down-expressed DEGs (light gray) in VC6089A were classified into 18 functional categories and chromosome locations. 122x177mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 2. The functional categories (A) and chromosome location (B) of 45 differential proteins (DPs) analyzed in two isogenic lines. The DPs with higher abundance (dark gray) and lower abundance (light gray) in VC6089A were classified into eight functional categories and chromosome locations. 122x177mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 40 of 44

Page 41 of 44

Journal of Agricultural and Food Chemistry

Figure 3. The bruchid-resistant candidate genes identified with transcriptomic and proteomic analyses. The 399 DEGs and 45 DPs were selected from RNA-seq and iTRAQ data between two isogenic lines. The three DEGs/DPs were identified from overlap of a Venn diagram and defined bruchid-resistant candidate genes. 62x46mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 4. Expression levels of bruchid-resistant candidate genes in resistant and susceptible lines. Gene expressions of g39185 (A), g34458 (B), and g5551 (C) were measured by qRT-PCR in resistant lines (RIL59, TC1966, and VC6089A) and susceptible lines (NM92 and VC1973A). Data represent the mean ± S.D. in three independent experiments. 60x21mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 42 of 44

Page 43 of 44

Journal of Agricultural and Food Chemistry

Figure 5. Map of the differential expressed genes (DEGs) and differential proteins (DPs) on chromosome 5. The DEGs located on chromosome 5 were represented on the left side. The DPs located on chromosome 5 were represented on the right side. The bruchid-resistant markers, DMB-SSR 158 and W02a4, were indicated by the letters with black color. The potential major Br gene was indicated in black color and with an asterisk. 122x177mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Exploring bruchid-resistant candidate genes in mungbean isogenic lines by omics technologies 47x26mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 44 of 44