Role of the Tomato Non-Ripening Mutation in Regulating Fruit Quality Elucidated Using iTRAQ Protein Profile Analysis

Natural mutants of the Non-ripening (Nor) gene repress the normal ripening of tomato fruit. The molecular mechanism of fruit ripening regulation by the Nor gene is unclear. To elucidate how the Nor gene can affect ripening and fruit quality at the protein level, we used the fruits of Nor mutants and wild-type Ailsa Craig (AC) to perform iTRAQ (isobaric tags for relative and absolute quantitation) analysis. The Nor mutation altered tomato fruit ripening and affected quality in various respects, including ethylene biosynthesis by down-regulating the abundance of 1-aminocyclopropane-1-carboxylic acid oxidase (ACO), pigment biosynthesis by repressing phytoene synthase 1 (PSY1), ζ-carotene isomerase (Z-ISO), chalcone synthase 1 (CHS1) and other proteins, enhancing fruit firmness by increasing the abundance of cellulose synthase protein, while reducing those of polygalacturonase 2 (PG2) and pectate lyase (PL), altering biosynthesis of nutrients such as carbohydrates, amino acids, and anthocyanins. Conversely, Nor mutation also enhanced the fruit’s resistance to some pathogens by up-regulating the expression of several genes associated with stress and defense. Therefore, the Nor gene is involved in the regulation of fruit ripening and quality. It is useful in the future as a means to improve fruit quality in tomato.


Introduction
The tomato (Solanum lycopersicum) is an important model to study climacteric fruit ripening because it has advantages, such as a simple genetic background, efficient genetic transformation method, fully sequenced genome, and large agricultural acreage [1,2]. Its fruit constitutes an important part of a healthy human diet because it contains high levels of nutrients, such as lycopene, vitamins, minerals, etc. Quality enhancement and nutrient synthesis in tomato fruits are closely related to its ripening process, which is regulated by various genes and environmental conditions [3,4]. For example, lycopene is present at extremely low levels in green tomato

Plant materials
The seeds of wild type (WT) tomato (Solanum lycopersicum 'Ailsa Craig') and its late-ripening mutant Nor were germinated and grown in the same green house at 25°C with 75% relative humidity under a16 h light/8 h dark regime. Flowers were gently shaken with an electronic toothbrush at full-boom to obtain pollen and were tagged one day post anthesis (DPA). Fruit samples were harvested at approximately 44 DPA when the WT fruit were at the breaker stage (BK). Immediately after harvesting, the pericarps were collected, snap-frozen in liquid nitrogen, and stored at-80°C until use.

Plasmid construction and virus-induced gene silencing (VIGS) assay
The tobacco rattle virus (TRV)-based vectors pTRV1 and pTRV2 were used for a VIGS of the tomato genes. A 540-bp fragment of the Nor gene corresponding to bases 617-1156 of the Nor gene sequence (NM_001247723.2) and a 421-bp fragment of the phytoene desaturase (PDS) gene corresponding to bases 1324-1744 of the PDS gene sequence (XM_010320112.1) were PCR-amplified from tomato cDNA using the primers listed in S1 Table. The resulting products were cloned into pTRV2 to form pTRV2-Nor or pTRV2-PDS.
The VIGS assay was carried out as previously described [9] with minor modifications. The Agrobacterium tumefaciens strain GV3101 containing the TRV-VIGS vectors was grown at 28°C in Luria-Bertani medium containing 10 mM MES and 20 mM acetosyringone with appropriate antibiotics. After incubation with shaking for 24 h, the cells were harvested and resuspended in infiltration buffer (10 mM MgCl 2 , 10 mM MES, pH 5.6, 200 mM acetosyringone) to a final OD 600 of 6.0 (for both pTRV1 or pTRV2 and its derivatives) and left to stand for 3-4 h at room temperature (20-25°C) before infiltration. Resuspensions of pTRV1 and pTRV2 or their derived vectors were mixed at a ratio of 1:1 and infiltrated into the carpopodium of the tomato fruits that were still attached to the plant at approximately 7-10 DPA with a 1-ml syringe. Tomato fruits infiltrated with pTRV2 or pTRV2-PDS alone were used as control.

RNA extraction and quantitative RT-PCR analysis
Total RNA extraction from the pericarp of the fruit was conducted using an RNeasy MiniKit (QIAGEN, Gemany) excluding polysaccharides and polyphenols from the sample according to the manufacturer's instructions. The extracted RNA was reverse-transcribed into cDNA using the one-step reverse-transcription system (TransGen Biotech, China) as described in the manufacturer's protocol.
Quantitative real-time PCR was performed using the SYBR Green PCR TOP Mix (Trans-Gen Biotech, China) along with a BIO-RAD real-time PCR system. The reaction program was as follows: 95°C for 10 min, 40 cycles at 95°C for 15 s and 60°C for 30 s. The fluorescence signal was measured automatically during each cycle. The relative quantification of specific mRNA levels were measured using the 2 (-ΔΔCt) analysis method, and expression values were normalized using the Actin gene and then normalized against the AC. The primers used in this study are listed in S1 Table. Three independent biological replicates were analyzed for each sample. homogenization for 3 min on ice and centrifugation at 15, 000 rpm for 10 min. The upper aqueous phase was removed and re-extracted using Buffer A. Proteins in the final phenol phase were precipitated overnight with four volumes of ice-cold Buffer C (saturated ammonium acetate in methanol) at-20°C. The proteins were pelleted by centrifugation and washed with ice-cold Buffer C three times and then twice in ice-cold acetone. Protein pellets were immediately solubilized in lysis buffer [8 M urea, 30 mM HEPES, 1 mM PMSF, 2 mM EDTA and 10 mM dithiothreitol (DTT)] using 5 min of sonication (pulse-on 2 s, pulse-off 3 s, power 180 W), followed by centrifugation at 20, 000 rpm for 30 min. The supernatant was collected, and DTT was added to a final concentration of 10 mM. Then, the mixture was incubated at 56°C for 1 h and immediately reduced by adding iodoacetamide (IAM) to a final concentration of 55 mM, this was followed by a 1 h incubation at room temperature in the dark. Proteins were precipitated overnight using four volumes of ice-cold acetone at-20°C and centrifuged at 20, 000 rpm for 30 min at 4°C. The precipitate was collected and diluted in 400 μl of dilution buffer to a final concentration of 50% triethyl ammonium bicarbonate (TEAB) and 0.1% SDS. The mixture was sonicated for 3 min as previously described and centrifuged at 20, 000 rpm for 30 min at 4°C. Then, the supernatant was collected and the protein concentration was determined using the Bradford assay.

Protein digestion and iTRAQ labeling
One hundred micrograms of protein from each sample were brought with 50% TEAB and 0.1% SDS to an equal volume and digested with 3.3 μg of trypsin at 37°C overnight. Then, 1 μg of trypsin was added and the proteins were digested again at 37°C for 12 h. The digested sample was lyophilized and resuspended in 30 μl of TEAB (water: TEAB = 1:1).
The tryptic peptides were labeled using the iTRAQ-8plex reagents according to the manufacturer's instructions. Samples of fruit tissue from the Nor mutant were labeled with iTRAQ tags 113, 114 and 115, and samples from the WT were labeled with iTRAQ tags 116, 117 and 118. Three independent biological replicates were performed. The incubation was processed at room temperature for 2 h, and the iTRAQ-labeled samples were subsequently vacuum dried.

Peptide fractionation with strong cation exchange (SCX) chromatography
For SCX chromatography, iTRAQ-labeled peptides were reconstituted with Buffer A (25% ACN, 10 mM KH 2 PO 4 , pH 3.0) to a final ten-fold volume and were centrifuged at 15, 000 rpm for 10 min. The supernatant was collected and loaded onto a Phenomenex Luna SCX column (100 A). The peptides were eluted at a flow rate of 1 ml/min using Buffer B (25% ACN, 2 M KCl, 10 mM KH 2 PO 4 , pH 3.0) as follows: a gradient of 0%-5% for 1 min, 5%-30% for 20 min, and 30%-50% for 5 min; maintained at 50% for 5 min; and then maintained at 50%-100% for 5 min. The system was maintained in 100% Buffer B for 10 min before equilibrating with Buffer A for approximately 10 min until 20 min before the next injection. The elution was monitored by measuring the absorbance at 214 nm, and fractions were collected every 1 min. The eluted peptides were finally combined into 12 pooled samples and desalted on a Strata XC18 column. The fractions were dried using low-temperature centrifugation and resuspended in 1% (v/v) formic acid.
column and separated at a flow rate of 400 ml/min using a mobile phase B gradient of 5% for 10 min, 5%-30% for 30 min, 30%-60% for 5 min, and 60%-80% for 3 min; maintained at 80% for 7 min; returned to 5% for 3 min; and maintained at 5% for 7 min. The eluted peptides were detected using Q-Exactive, and MS data were acquired using a data-dependent top20 method by choosing the abundant precursor ions from the survey scan (350-20, 000 Da) using higher energy collision dissociation (HCD). Determination of the target value was based on automatic gain control (AGC). Survey scans were acquired at a resolution of 70, 000 and the resolution for the HCD spectra was set to 17, 500.

Protein identification and quantification
Mascot version 2.3.0 and Thermo fisher Proteome Discoverer software version 1.3 were used to identify and quantify proteins. To enhance the precision of identification and quantification, only unique peptides used for protein quantification were selected. The raw MS data were first filtered using the Proteome Discoverer software with the following parameters: molecule weight of precursor ions set between 350-6, 000 Da; minimum peaks in MS/MS set to 10; and signal-to-noise threshold set to 1.5. Searches using Mascot were made against the database tomato uni4081. The following parameters were used in the database searches: tolerance for peptides and MS/MS: 15 ppm and 20 mmu, respectively; fixed modification: carbamidomethyl (C); variable modification: oxidation (M), Gln!Pyro-Glu (N-term Q), iTRAQ 8 plex(K), iTRAQ 8 plex (Y), iTRAQ 8 plex (N-term); maximum missed cleavages 1; and enzyme: Trypsin. The FDR was less than 1%. The following parameters were used for the quantification analysis using Proteome Discoverer: protein ratio type: median; minimum peptides: 1; normalization method: median; P-value<0.05; and fold change>1.4.

Proteomic data analysis
Proteins with a 1.4-fold change and a P-value<0.05 were identified as differentially abundant proteins. Gene Ontology (GO) analysis and a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the differentially abundant proteins were conducted using the Gene Ontology database (http://www.geneontology.org/) and the KEGG pathway database (http://www.kegg.jp/kegg/pathway.html), respectively. Based on the acquired GO terms, the WEGO (http://wego.genomics.org.cn/) was used for GO enrichment analysis and the proteins were filtered corresponding to cellular components, molecular functions and biological processes.

Results and Discussion
Strong inhibition of fruit ripening as a result of silencing the Nor gene With the completion of tomato genome sequencing, researchers soon discovered that the Nor mutant resulted from a mutation in a NAC family gene; this gene was cloned and its sequence was published in GenBank. To verify whether the fruit phenotype of the plants subjected to Nor gene silencing is consistent with that of the Nor mutant, VIGS was used [9] to suppress the Nor gene in the tomato plant. The 540-bp cDNA fragment of the Nor gene conserved region was PCR-amplified and cloned into the pTRV2 vector to generate pTRV2-Nor constructs. A mixture of A. tumefaciens (GV3101) cultures containing pTRV1 and pTRV2-Nor constructs at a 1:1 ratio were injected using a needle into the carpopodium of wild-type Ailsa Craig tomato fruits attached to the plant 7-10 days after pollination. Plants infiltrated with A. tumefaciens cultures, including pTRV1 and pTRV2-00, were used as controls ( Fig 1B). The phenotype of the fruit was observed 25-35 days later. Twenty fruits (approximately 70%) injected with pTRV1 and pTRV2-Nor produced regions that failed to ripen normally and displayed a distinct green color, which is indicative of delayed ripening, or produced yellow (partially ripened) fruits at a time when parts of the fruits not subjected to silencing turned red. The green region, indicating inhibited ripening on the Nor fruit subjected to VIGS, eventually turned orange and showed signs of slow ripening. All control fruits injected with the pTRV1 and pTRV2-00 vector ripened normally (Fig 1A).
To confirm Nor gene suppression at the molecular level, primers specific to the Nor genes outside the region targeted for silencing were designed for real-time PCR. The reduction in the Nor transcripts in the silenced green part was 86% compared with the control fruit (in the parts that had the normal red color). Levels of the ACT1 transcript were similar in tissues infected with pTRV2-Nor and pTRV2-00 vector constructs ( Fig 1C). Tomato fruit infected with pTRV1 and pTRV2-PDS turned yellow, as observed for the positive control. Tomato fruits infected with pTRV1 and pTRV2-00 turned red as observed for the empty control ( Fig 1A). These results indicate that the Nor gene plays an important role in controlling the ripening of tomato fruits, but its mechanism of regulating fruit ripening was not clear. This was further studied in the next set of experiments.

Primary data analysis and protein identification in the Nor mutant
To further unravel the mechanism of Nor regulation of fruit ripening at the level of protein, the iTRAQ technique will be performed using WT and Nor mutant tomato fruit. First, we need to make sure that Nor expression cannot be detected in Nor mutant. Total RNA from the break stage (BK) of the WT and Nor mutant tomato fruit was prepared and used to reverse transcript to cDNA for Nor gene RT-PCR. The results showed that the expression of the Nor gene was barely detected in the Nor mutant (S1 Fig). This indicates that the material can be used for the iTRAQ analysis of the Nor mutant. Next, total proteins were extracted from the pericarp tissue of Nor mutant and WT fruits and were analyzed using the iTRAQ technique at BK (Fig 2A).
After pooling data from three biological replicates, a total of 333, 859 spectra were generated from the iTRAQ analysis using the control and Nor mutant fruits at BK. Mascot identified a total of 81, 183 spectra matched to known spectra, which were further matched to 22, 661 peptides and finally matched to 3, 734 proteins with a FDR of less than 1% (Fig 2B). The distribution of mass coverage of proteins is provided in Fig 2C and 2D.
When screening for differentially abundant proteins between WT and the Nor mutant, a protein species was considered differentially accumulated if it exhibited a fold change>1.40 and a P-value<0.05. Using these two criteria, 218 differential abundant protein species including 131 uncharacterized proteins were screened from the Nor mutant sample compared with the wild type, of which 123 proteins were up regulated and 95 were down regulated (Fig 2D). Detailed information is provided in Table 1 and in S2 Table. Bioinformatic analysis of differentially abundant protein species (DAPS) identified using iTRAQ To identify significantly enriched GO functional groups of DAPS, GO annotation was carried out. A total of 87 DAPS (40%) between the Nor mutant and control were classified into 24 functional groups (Fig 3), of which biological processes accounted for 10 GO terms (the most representative were "metabolic process"), cellular components accounted for 7 GO terms (the most representative were "cell and cell part"), and molecular functions accounted for 7 GO terms (the most representative was "catalytic").
To further investigate the biological functions of these proteins, 87 (40%) were mapped to 20 pathways in the KEGG database (Table 1). "carbohydrate metabolic process" was the most represented pathway followed by "defense response" and "amino acid metabolic process" (Fig 4).
To verify whether the differentially abundant proteins have the same tendency as the expression of their target genes at the transcript level, real-time PCR was used to detect the expression of target genes in the Nor mutant and wild-type tomato fruit. Because we wanted to examine the pathways related to quality formation and ripening of Nor mutant samples (which we discuss in the next section), the target genes selected for analysis included genes related to ethylene biosynthesis (4), fruit softening (4), pigment synthesis (6), nutrient synthesis (11), and stress and defense pathways (2). iTRAQ data for the proteins encoded by these genes indicated that some were up regulated, such as the cellulose synthase gene, whereas some were down regulated, such as PSY1 and CHS1. The results of real-time PCR indicate that the trend of most of these gene expressions was consistent with that of the differentially abundant proteins (Fig 5), which indicates the reliability of our iTRAQ data and that these data can be used for further analysis.

Ethylene biosynthesis
Ethylene is a key factor promoting the ripening of climacteric fruits, such as tomatoes [4]. Methionine is the direct precursor of the ethylene biosynthesis pathway; ACC synthase and ACC oxidase catalyze two key steps that are coded by the gene family [6]. The expression of the SlACS2 and SlACO4 genes increases in response to considerable ethylene synthesis during fruit ripening [6]. In addition, SAM synthase (SAMS) catalyzes ethylene synthesis from methionine in the first step of the reaction; silencing this gene can effectively reduce ethylene biosynthesis and considerably alter the ripening of tomato fruit compared with the control [6]. Our protein data indicated that the relative abundances of the four ACO genes in Nor tomato fruit were 66%, 53%, 30%, and 51%, respectively, compared with the control fruits, whereas the abundance of two SAMS increased by 153% and 147%, respectively. Real-time PCR revealed that the expression of the four ACO genes in the Nor mutant was significantly less than that in the control fruit. Ethylene synthesis barely occurred in the Nor tomato fruit (data not shown). These results indicate that the Nor gene regulates the expression of the SAM1 and ACO genes, which are involved in the regulation of fruit ripening (Fig 6).

Fruit softening
With the ripening of tomato fruit, the firmness of the fruit also changes significantly, primarily due to the reduction in cellulose synthesis and pectin degradation. Firmness is an important factor affecting the storage period and quality of fruit after fruit harvest [4]. The process of fruit  Cellulose is the main structural component of the plant cell wall and is typically cross-linked with hemicellulose, pectin, and lignin, which all enhance fruit firmness [18]. Plant cellulose synthases belong to the family of glycosyltransferases, which includes proteins involved in cellulose biosynthesis [19]. The two cellulose synthase proteins in Nor mutant are more abundant than that of wild type: 1.96-and 1.64-fold, respectively, respectively. This indicates that the Nor mutant can synthesize considerably more cellulose than the wild type, which has positive implications for fruit firmness.
Polygalacturonase, which is also known as PG, is an enzyme that hydrolyzes the alpha-1, 4 glycosidic bonds between galacturonic acid residues [20]. The activity of the endogenous plant PGs results in the softening and sweetening of the fruit during the ripening process [21]. During the ripening process, PG gene expression significantly increases, promoting the softening of the fruit. Pectate lyase is an enzyme that is involved in the maceration and soft rotting of plant tissue and is responsible for the eliminative cleavage of pectate, yielding oligosaccharides with 4-deoxy-α-D-mann-4-enuronosyl groups at their non-reducing ends. The protein is maximally expressed in the late stages of fruit ripening [22]. The expression of pectate lyase (PL) genes has been suggested to be related to a requirement for fruit softening. Our iTRAQ data show that the levels of PG and PL protein in the Nor mutant were significantly less than in the control fruit: 0.31-and 0.37-fold less, respectively.
Based on these data, we suggest that the Nor mutation significantly enhances tomato fruit firmness by promoting the expression of the cellulose synthase gene and inhibiting the expression of the PG and PL genes (Fig 6).

Fruit nutrition
The process of fruit ripening is accompanied by the synthesis and degradation of specific nutrients that provide good quality food for humans. The changes that occur during fruit ripening include lycopene synthesis, chlorophyll degradation, synthesis of sugars, starch degradation, flavor compound synthesis, and organic acid degradation, as well as changes in amino acids and synthesis of secondary metabolites [3]. As a result of the silencing or mutation of key ripening-associated genes, such as Rin, Cnr, apetala 2a (AP2a), tomato agamous-like1 (TAGL1) and fruitfrull1 (FUL1/TDR4), the fruit ripening process was significantly inhibited [2]. The metabolic pathways were constructed based on the KEGG pathway analysis. Dashed arrows represent the multiple steps between two compounds, the red colored text indicates the accumulated proteins, the green colored text indicates decreased expressed proteins in the Nor mutant tomato fruits, and the boldfaced text represents the proteins listed in Table 1. Synthesis of the nutrients described above was also considerably affected. Silencing of the TAGL1 gene resulted in a decrease in the lycopene, lutein and chlorophyll contents of tomato fruit [23]. Tomato fruit with antisense TDR4 reduced the lycopene levels by down regulating the expression of 1-deoxy-D-xylulose 5-phosphate synthase 1 (DXS1), which affected lipid and cuticle metabolism, resulting in fruit wilting and a significant decrease in free glutamic acid by up regulating the expression of the glutamic acid decarboxylase gene [24]. Rin regulates the tomato lipoxygenase (TomloxC), polyneuridine aldehyde esterase (PNAE), phosphoglycerate kinase (PGK) and alcohol dehydrogenase 2 (ADH2) genes in the plant lipoxygenase (LOX) pathway, which affect flavor compound synthesis in tomato fruit [16]. The golden-like2 (GLK2) gene regulates the synthesis of fruit nutrients, such as lycopene, β-carotene, lutein, ascorbic acid, fructose, glucose, sucrose, and starch, as well as the synthesis of other nutrients. Together, the data indicate that ripening-associated transcription factors regulate the synthesis of important nutrients during the tomato ripening process [25]. Mutation in the Nor gene hampers the ripening of the fruit. To determine whether this mutation also affects the quality of the fruit with respect to the levels of nutrients, we analyzed the levels of differentially abundant proteins in Nor mutant and wild-type fruit. The mutation in the Nor gene was found to affect the synthesis of several nutrient biomolecules and enzymes involved in various metabolic processes ( Table 1, Fig 6), including carbohydrates (13), lipids (6), amino acids (8), secondary metabolites of brass (3), etc. The real-time PCR results of these genes encoding target proteins were consistent with the observed trend of protein expression. The results showed that the Nor gene was involved in the synthesis of nutrients during fruit quality enhancement and played an important role in the ripening and quality enhancement of tomato fruit.

Pigment synthesis
The pigmentation of ripe tomato fruit is gradually generated by the accumulation of pigments, such as carotenoids and flavonoids, along with the degradation of chlorophyll [26]. In our study, the protein level of eight pigmentation related proteins was significantly less in Nor mutant fruits compared with the controls; two of these proteins, PSY1 and Z-ISO, take part in carotenoid biosynthesis: PSY1 is an enzyme catalyzing the synthesis of phytoene from geranylgeranyl pyrophosphate (GGPP). Down-regulation of the PSY1 gene by transgenic technology or VIGS eliminates normal carotenoid formation [27]. Z-ISO gene is highly expressed in tomato fruit, and silencing this gene using VIGS leads to a dramatic reduction in lycopene content [27]. These results indicate that carotenoid biosynthesis was affected by Nor gene (Fig 6).
In addition, six proteins related to flavonoid biosynthesis were repressed in Nor mutant fruit: CHS1, chalcone-flavonone isomerase family protein (CHI), trans-cinnamate-4-monooxygenase, and three glycosyltransferase proteins. The repression of these proteins essentially abolished the flavonoid biosynthesis pathway (Fig 6). Moreover, the biosynthesis of phenylpropanoid, which is the precursor of flavonoid, was also suppressed, as revealed using KEGG pathway analysis. This evidence suggests that Nor functions is a positive regulator in these processes. In a previous study, a reduction in the level of phenolics was detected in the Nor fruit; these observations are consistent with our analysis [28]. CHS1 protein was significantly repressed (0.55-fold as revealed by iTRAQ and 0.005 as revealed by RT-PCR) in the Nor mutant. CHS1 down regulation reportedly reduces the level of yellow-pigmented naringenin chalcone, which is generally present at high levels in the epidermal cells of ripening fruit [29]. Notably, some of the ripening-associated transcription factors are suggested to function in the activation of the flavonoid pathway, such as TDR4 and TAGL1 [30,31]. Both TDR4 and TAGL1 RNAi tomato fruit exhibited a similar orange-yellow fruit phenotype with reduced the levels of CHS1 expression and naringenin chalcone accumulation in the cuticle. There are a few reports indicating the relationship between flavonoid biosynthesis and fruit ripening, but the mechanism is still unclear.

Stress and Defense Response
The stress and defense responses are plant defense mechanisms that act against invasion by bacteria and fungi (biotic stress) or as adaptations to changes in environmental conditions (abiotic stress). Pathogenesis-related proteins (PRs) are defined as a family of plant proteins that are induced and accumulate in plant tissues under pathological or related situations [32]. PRs possess antimicrobial activities in vitro through hydrolytic attacks on cell walls, contact toxicity, and their putative involvement in defense signaling; PRs have been classified into 17 families [33]. Our iTRAQ data revealed that seven PRs were differently displayed in Nor mutant fruits, of which five proteins were depressed to various extents: glucan endo-1, 3-β-glucosidase B (0.47-fold), an osmotin-like protein (NP24) (0.47-fold), a pathogenesis-related protein (PR1) (0.66-fold), a class I chitinase (CHI9) (0.53-fold), and peroxidase (0.66-fold). Glucan endo-1,3-β-glucosidase B and CHI9 can hydrolyze glucan or chitin, respectively, which degrade bacterial and fungal cell wall and inhabit pathogen growth [34]. Moreover, the protein level of glucan endo-1,3-β-glucosidase is reported to be regulated by Bradyrhizobium japonicum inoculation in seedlings under flooding [35]. NP24, which is a thaumatin-like protein in tomato fruit, has been tested as a growth inhibitor against six different fungi [36] and is up regulated in response to salt or osmotic stress [37,38]. Peroxidase belongs to the family of class III plant peroxidases, which partake in a wide-range of physiological processes, for instance, lignin and suberin formation, cell wall component cross linking, and phytoalexins synthesis, or participate in reaction oxygen species (ROS) and reactive nitrogen species (RNS) metabolism, which "turn on" the hypersensitive response [39]. PR1 is part of the PR-1 protein family, which is in addition to the known PR-1a, PR-1b, and PR-1c proteins found in tomato. The main PR-1 proteins induced by salicylic acid (SA) or pathogens are commonly used as markers of the state of systemic acquired resistance (SAR), but their biological function is still unknown [33,40]. In contrast, two non-specific lipid-transfer proteins (LTPs) are considerably up regulated by 3.39-fold and 1.65-fold. These proteins belong to the PR-14 family, whose members are approximately 70 and 90 amino-acid-(aa)-long cationic proteins with eight Cys residues. According to a previous study, LTPs are suggested to be extracellular proteins and promote the formation of pores in the membranes of pathogens rather than hosts cells, which inhibits the growth of pathogens, and are involved in the plant defense system [41].
Small heat shock proteins (sHSPs) weighing between 12 and 42 kDa are the most ubiquitous subgroup of heat shock protein (HSPs). They are encoded in response to heat stress and are commonly not found in normal plant tissues, whereas the accumulation of some sHSPs has been detected during plant developmental processes, for example, pollen and embryo development, seed germination, and fruit ripening [42]. Our iTRAQ data suggest that all five sHSPs protein significantly accumulated in the Nor mutant. HSP17.4, which is regulated by seed development and stress, is associated with seed desiccation tolerance in Arabidopsis [43]. HSP17.6 can be induced by Botrytis cinerea in red ripe tomato fruit instead of in mature green fruit, which indicates the ripening-regulated susceptibility of tomato fruit to this fungus [44]. HSP17.7 is related to the UVA, UVB, Cu 2+ , and Zn 2+ stress responses in sunflowers [45]. Moreover, MTSHP, which is a mitochondrial small heat shock protein in tomato, has been reported to possess molecular chaperone activity in vitro [46]. Additionally, viscosity 1 (VIS1) contributes to pectin depolymerization and juice viscosity in tomato [42], leading to speculation that other sHSPs could play a role in fruit ripening.
Snakin-2 is a member of the snakin peptide family, which is characterized by maintaining 12 cysteines in the highly conserved C-terminus, forming six disulfide bonds [47].
Overexpression of Snakin-2 restricts pathogen invasiveness and enhances tolerance to Clavibacter michiganensis in transgenic tomato plants [48]. The Snakin-2 gene is induced by wounding and responds to pathogen invasion in potatoes [49]. Snakin-2 protein harbors antimicrobial activity against bacteria and fungi, and this activity is mediated through the nonspecific perforation of the microbe's biomembrane [47,50]. Here, Snakin-2 shows a significantly increased level (1.86-fold) in the Nor mutant, indicating a changed tolerance to pathogen infection.
The methionine sulfoxide reductase (Msr) family contains ubiquitous enzymes that function to protect cells against oxidative damage [51]. The two key enzymes in the Msr family are MSRA and MSRB, which can reduce the S and R epimers of methionine sulfoxide to methionine, respectively [52]. Both enzymes have been shown to defend bacteria, yeast, and plant cells against the cytotoxic effects of ROS, which prevents undue amassing of oxidized proteins and aging-associated death [53]. In tomato, the fruit ripening gene E4, which is strongly induced by ethylene, encodes a MSRA protein [54]. This MSRA protein displayed significantly depressed accumulation (0.38-fold), which is likely attributed to the elimination of ethylene in the Nor mutant fruit.

Conclusions
Ripening is the crucial stage of fruit development that prepares the seed for reproduction and confers fruit quality traits, such as color, aroma, taste, and texture, to ultimately provide food fit for human consumption. Therefore, some important genes controlling fruit ripening will also affect fruit quality, including nutrient synthesis. The tomato fruit has several well-known ripening-repressed natural mutations, such as Rin, Cnr, Nr, and Nor. The involvement of these three genes in the regulation of fruit ripening and formation and in the development of a goodquality fruit has been extensively studied, but there is limited knowledge regarding the involvement of the Nor gene. Here, we analyzed differences in the abundances of proteins in Nor mutant and wild-type fruit using the iTRAQ proteomic technique. The results indicate that the Nor mutation has a significant effect on ethylene synthesis, pigment biosynthesis, softening, and nutrient biosynthesis, as well as impacts disease resistance-associated genes. Moreover, our study data show that Nor plays an important role in the development of fruit quality. We should, however, further verify whether the Nor gene, as a transcription factor, directly regulates the transcription of the genes encoding the differentially abundant proteins in the Nor mutant of tomato fruit.