ORIGINAL RESEARCH article
Sec. Plant Pathogen Interactions
Volume 12 - 2021 | https://doi.org/10.3389/fpls.2021.613568
Oryza sativa cv. Nipponbare and Oryza barthii as Unexpected Tolerance and Susceptibility Sources Against Schizotetranychus oryzae (Acari: Tetranychidae) Mite Infestation
- 1Graduate Program in Biotechnology, University of Taquari Valley-Univates, Lajeado, Brazil
- 2Biological Sciences and Health Center, University of Taquari Valley-Univates, Lajeado, Brazil
- 3Graduate Program in Botany, Federal University of Rio Grande do Sul, Porto Alegre, Brazil
- 4Graduate Program in Biotechnology, University of Caxias do Sul, Caxias do Sul, Brazil
- 5Graduate Program in Molecular and Cellular Biology, Federal University of Rio Grande do Sul, Porto Alegre, Brazil
- 6Laboratory of Biotechnology, Bioscience and Biotechnology Center, State University of Northern Rio de Janeiro Darcy Ribeiro, Campos dos Goytacazes, Brazil
- 7Rice Research Institute (Instituto Rio-Grandense do Arroz), Cachoeirinha, Brazil
Cultivated rice (Oryza sativa L.) is frequently exposed to multiple stresses, including Schizotetranychus oryzae mite infestation. Rice domestication has narrowed the genetic diversity of the species, leading to a wide susceptibility. This work aimed to analyze the response of two African rice species (Oryza barthii and Oryza glaberrima), weedy rice (O. sativa f. spontanea), and O. sativa cv. Nipponbare to S. oryzae infestation. Surprisingly, leaf damage, histochemistry, and chlorophyll concentration/fluorescence indicated that the African species present a higher level of leaf damage, increased accumulation of H2O2, and lower photosynthetic capacity when compared to O. sativa plants under infested conditions. Infestation decreased tiller number, except in Nipponbare, and caused the death of O. barthii and O. glaberrima plants during the reproductive stage. While infestation did not affect the weight of 1,000 grains in both O. sativa, the number of panicles per plant was affected only in O. sativa f. spontanea, and the percentage of full seeds per panicle and seed length were increased only in Nipponbare. Using proteomic analysis, we identified 195 differentially abundant proteins when comparing susceptible (O. barthii) and tolerant (Nipponbare) plants under control and infested conditions. O. barthii presents a less abundant antioxidant arsenal and is unable to modulate proteins involved in general metabolism and energy production under infested condition. Nipponbare presents high abundance of detoxification-related proteins, general metabolic processes, and energy production, suggesting that the primary metabolism is maintained more active compared to O. barthii under infested condition. Also, under infested conditions, Nipponbare presents higher levels of proline and a greater abundance of defense-related proteins, such as osmotin, ricin B-like lectin, and protease inhibitors (PIs). These differentially abundant proteins can be used as biotechnological tools in breeding programs aiming at increased tolerance to mite infestation.
Rice (Oryza sativa L.) is extremely important for human nutrition, representing around 20% of the daily intake calories in the world. It is considered a staple food for over half of the world’s population (Muthayya et al., 2014). To meet this nutritional demand, natural resources must be used efficiently, and crop production needs to be increased (Godfray et al., 2010). However, the search for new agricultural frontiers and the availability of natural resources is limited, exposing rice cultivation to several disturbances caused by biotic and abiotic stresses, which negatively impact grain yield (Palmgren et al., 2015). Among these important stresses that influence plant yield, the presence of arthropods (Oerke, 2006) impacts about 18–26% of annual crop production worldwide, causing losses of over US$ 470 billion (Culliney, 2014). Most of these losses (13–16%) occur in the field before harvesting, and losses are most often described in developing countries (Culliney, 2014). One of the arthropods that attack rice crops in Brazil, and has been reported in other South American countries, is Schizotetranychus oryzae Rossi de Simons mite species, which can cause more than 60% loss in rice grain yield (Buffon et al., 2018). An interesting strategy for reducing these losses is the search for mite-tolerant cultivars, as these plants tolerate mite infestation and maintain grain yield without the need to apply acaricides (Sperotto et al., 2018a).
Plant defense to herbivory includes both resistance (the ability of plants to escape attacking enemies) and tolerance (the ability of plants to withstand herbivory without any decline in yield) (Stenberg and Muola, 2017). Both defense strategies can be plant focused and can be useful to measure the plant’s ability to defend from herbivores. Several molecular and physiological characteristics can be used to detect resistance/tolerance plant capacity, including expression of defense-related genes, levels of soluble protein/sugar, reactive oxygen species, lipid peroxidation, carotenoids/chlorophyll, photosynthetic efficiency, cell membrane damage, and others (Blasi et al., 2015).
The search for genetic variability in the Oryza genus is a shortcut to obtain resistant/tolerant rice genotypes and to develop cultivars that adapt to natural adversities (Menguer et al., 2017; Szareski et al., 2018). Several authors reported wild rice plants resistant/tolerant to biotic and abiotic factors (Brar and Khush, 2003). The use of wild rice species in breeding programs can facilitate adaptation to biotic and abiotic stresses, as well as meeting the demand for food security in the current scenario of rapidly growing world population (Henry, 2014). The Oryza genus is composed of 24 species, two cultivated (O. sativa and Oryza glaberrima) and other 22 wild species. These wild species represent between 15 and 25 million years of evolutionary diversification (Vaughan, 1994; Menguer et al., 2017).
Oryza barthii is an annual wild African species that is commonly known as the progenitor of O. glaberrima, the rice cultivated species grown in Africa. O. barthii is tolerant to several biotic and abiotic stresses (Khush, 1997), being adapted to more adverse ecological conditions, and resistant to multiple environmental restrictions (Sarla and Swamy, 2005). Weedy rice (O. sativa f. spontanea), which is the result of rice de-domestication process, also shows tolerance to environmental stresses. Therefore, the genetic background, morphology, and growth behavior are similar to cultivated rice (Qiu et al., 2017). However, differences in tolerance to stressful conditions may occur during de-domestication, e.g., weedy rice is tolerant to cold, high salinity and drought, and resistant to blast disease. Due to its greater stress tolerance/resistance, weedy rice has become one of the most feared and harmful weeds in rice producing regions worldwide (Dai et al., 2014). Nonetheless, interesting traits found in weedy rice can be eventually transferred to cultivated genotypes.
Aiming to identify novel resistance/tolerance mechanisms in rice plants exposed to S. oryzae mite infestation, we analyzed two African species (wild O. barthii and cultivated O. glaberrima), weedy rice (O. sativa f. spontanea) and cv. Nipponbare (O. sativa ssp. japonica). Although we expected wild rice species to show higher tolerance to mite infestation, we found unexpected results. According to our data, O. barthii and O. glaberrima are extremely sensitive to mite infestation, while cv. Nipponbare seems to be more tolerant. Therefore, we tried to understand the molecular and physiological mechanisms behind Nipponbare tolerance and African species susceptibility to this mite. Our results may be useful for future breeding programs aiming at resistance/tolerance to S. oryzae mite infestation.
Materials and Methods
Plant Growth Conditions and Mite Infestation
Seeds of O. sativa cv. Nipponbare (hereafter “Nipponbare”) were obtained from the seed bank Rice Genome Resource Center1. Seeds of weedy and African rice species are part of the International Rice Research Institute (IRRI) germplasm collection: O. sativa f. spontanea (IRGC 80590, hereafter “weedy rice”), O. barthii (IRGC 86524), and O. glaberrima (IRGC 103959). Seeds were surface sterilized and germinated for 4 days in an incubator (28°C) on paper soaked with distilled water. After germination, plantlets were transferred to vermiculite/soil mixture (1:3) for additional 14 days in greenhouse conditions, and then transferred to plastic buckets containing soil and water. Plastic buckets containing rice plants highly infested by S. oryzae were kindly provided by Instituto Rio-Grandense do Arroz (IRGA, Cachoeirinha, RS), and were used to infest rice plants in our experiment. Fifty plants (V10-13 stage, according to Counce et al., 2000) of each cultivar/species (five plants per bucket) were infested by proximity with the bucket containing the highly infested plants placed in the center of the other buckets. For greater homogeneity of infestation and contact, buckets of each cultivar/species were rotated at a 90° angle counterclockwise every 2 days. Ten plants of each cultivar were cultivated without infestation (control condition).
The level of damage caused by S. oryzae was analyzed from V10-13 stage until the plants reached the final stage of reproductive development (panicle maturity, R9 stage; Counce et al., 2000). Evaluation of damage in the abaxial and adaxial faces of leaves was based on a classification of four infestation levels: Level 1: control condition, without any sign of infestation; Level 2: early infested (EI) leaves, 10–20% of damaged leaf area, average of 168 h of exposure to the mite; Level 3: intermediate infested (II) leaves, 40–50% of damaged leaf area, average of 360 h; and Level 4: late infested (LI) leaves, more than 80% of damaged leaf area, average of 720 h, according to Buffon et al. (2018).
Plant Height and Tiller Number
Plant height and tiller number were evaluated (n = 50) during the vegetative stage (V10-13, before infestation) and during the last reproductive stage (R9, control and infested plants).
Chlorophyll a Fluorescence Transients
The chlorophyll a fluorescence transient was measured (n = 10) on the third upper leaves of control and infested plants in two different exposure times, 1 week (which correspond to EI) and 4 weeks (which correspond to LI), using a portable fluorometer (OS30p, Optisciences, United Kingdom). Before the measurements, plants were dark adapted for 20 min, and the fluorescence intensity was measured by applying a saturating pulse of 3,000 μmol photons m–2 s–1. The resulting fluorescence of the chlorophyll a was measured from 0 to 1 s. The chlorophyll fluorescence intensity rises from a minimum level (the O level, FO = 160 μs) to a maximum level (the P level, FP = 300 ms), via two intermediate steps labeled J (FJ = 2 ms) and I (FI = 30 ms) phases (Stirbet and Govindjee, 2011), also known as OJIP curve (Strasser et al., 2000).
Total Chlorophyll Concentration
Samples (n = 3) of leaves from rice plants submitted to control and LI conditions were collected, and total chlorophyll concentration was quantified. In brief, 100 mg of sample was sprayed in liquid nitrogen, homogenized in 1.5 mL of 85% (v/v) acetone, maintained in the dark for 30 min, and centrifuged at 12,000 g for 3 min. The supernatant was transferred to a new tube, and the process was repeated until complete discoloration of the plant tissue. Chlorophyll a and b were quantified by measuring absorbance at 663 and 645 nm, and the concentrations were calculated according to Ross (1974).
In situ Histochemical Localization of H2O2
In situ accumulation of H2O2 in control and LI leaves was detected by histochemical staining with DAB, according to Shi et al. (2010), with minor modifications. Leaves were excised and immersed in DAB solution (1 mg mL–1, pH 3.8) in 10 mM phosphate buffer (pH 7.8) and incubated at room temperature for 8 h in the light until brown spots were visible, which are derived from the reaction of DAB with H2O2. Leaves were bleached in boiling concentrated ethanol to visualize the brown spots, which were kept in 70% ethanol for taking pictures with a digital camera coupled to a stereomicroscope.
Seeds of Nipponbare and weedy rice were collected in R9 stage, and the following agronomical parameters were evaluated: number of panicles per plant, number of seeds (empty + full) per panicle, percentage of full seeds per panicle, grain length, and weight of 1,000 full grains. Yield reduction caused by S. oryzae infestation in Nipponbare was calculated using the following equation, for each condition (control and infested): number of seeds (empty + full) per plant × percentage of full seeds × weight of one seed = seed weight per plant. Seeds from African rice species were not analyzed as these plants did not reach the reproductive stage.
Control and LI leaves of Nipponbare and weedy rice were collected, and proline accumulation was quantified (n = 3) according to Shukla et al. (2012). In brief, 250 mg of sample was sprayed in liquid nitrogen, homogenized in 2 mL of sulfosalicylic acid (3%), and centrifuged at 10,000 g for 30 min. Then, 1 mL of the supernatant was transferred to a new tube, and 1 mL of glacial acetic acid and 1 mL of ninhydrin (1.25 g of ninhydrin + 30 mL of glacial acetic acid + 20 mL of 6 M phosphoric acid) were added. The homogenized mixture was boiled in a water bath at 100°C for 30 min. Subsequently, the reaction was cooled on ice and 2 mL of toluene was added, mixing in the vortex for 30 s. The upper phase containing proline was measured on the spectrophotometer at 520 nm. The proline level (mmol g–1 fresh weight) was quantified using L-proline as standard. The standard curve was prepared from a stock solution of proline (5 mM) in sulfosalicylic acid for a concentration range between 10 and 320 mmol mL–1. The free proline concentration of each sample was calculated from the linear regression equation obtained from the standard curve. The procedure for constructing the standard curve was the same as for free proline determination of the samples.
Plant Protein Extraction and Quantification
Three biological samples (250 mg of fresh weight) of control and EI leaves from O. barthii (mite-sensitive) and Nipponbare (mite-tolerant), each containing three leaves from three different plants, were subjected to protein extraction using Plant Total Protein Extraction Kit (Sigma–Aldrich). The protein concentration was measured using 2-D Quant Kit (GE Healthcare, Piscataway, NJ, United States).
For protein digestion, three biological replicates of 100 μg of proteins from O. barthii and Nipponbare leaves were used. Before the trypsin digestion step, protein samples were precipitated using the methanol/chloroform methodology to remove any detergent from samples (Nanjo et al., 2012), followed by resuspension in urea 7 M and thiourea 2 M buffer in 50 mM ammonium bicarbonate. Right after, protein digestion was performed using the filter-aided sample preparation (FASP) methodology (Wiśniewski et al., 2009), with modifications described by Burrieza et al. (2019). In brief, after checking the integrity of the Microcon-30 kDa (Merck Millipore, Germany) filter units, protein aliquots were added to the filter units, washed with 200 μL of 50 mM ammonium bicarbonate (solution A), and centrifuged at 10,000 g for 15 min at 25°C (unless otherwise stated, all centrifugation steps were performed under this condition). This step was repeated once for complete removal of urea before reduction of proteins. Next, 100 μL of 50 mM DTT (GE Healthcare), freshly made in solution A, was added, gently vortexed, and incubated for 20 min at 60°C (1 min agitation and 4 min resting, at 650 r/min). Then, 200 μL of 8 M urea in 50 mM ammonium bicarbonate (solution B) was added, and centrifuged for 15 min. For protein alkylation, 100 μL of 50 mM iodoacetamide (GE Healthcare), freshly prepared in solution B, was added, gently vortexed, and incubated for 20 min at 25°C in the dark (1 min agitation and 19 min resting, at 650 r/min). Next, 200 μL of solution B was added and centrifuged for 15 min. This step was repeated once. Then, 200 μL of solution A was added and centrifuged for 15 min. This step was repeated twice. Approximately 50 μL of the sample should remain in the last washing. For protein digestion, 25 μL of 0.2% (v/v) RapiGest (Waters, Milford, CT, United States) and 25 μL of trypsin solution (1:100 enzyme:protein, V5111, Promega, Madison, WI, United States) were added, gently vortexed, and incubated for 18 h at 37°C (1 min agitation and 4 min resting, at 650 r/min). For peptide elution, the filter units were transferred to new microtubes and centrifuged for 10 min. Then, 50 μL of solution A was added and centrifuged for 15 min. This step was repeated once. For RapiGest precipitation and trypsin inhibition, 5 μL of 15% trifluoroacetic acid (Sigma–Aldrich) was added, gently vortexed, and incubated for 30 min at 37°C. Then, samples were centrifuged for 15 min, and the supernatants collected and vacuum dried. Peptides were resuspended in 100 μL solution of 95% 50 mM ammonium bicarbonate, 5% acetonitrile, and 0.1% formic acid. The resulting peptides were quantified by the A205 nm protein and peptide methodology, using a NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific).
Mass Spectrometry Analysis
A nanoAcquity UPLC connected to a Synapt G2-Si HDMS mass spectrometer (Waters, Manchester, United Kingdom) was used for ESI-LC-MS/MS analysis. Runs consisted of injecting 1 μg of digested peptides from each biological replicate. During separation, samples were loaded onto the nanoAcquity UPLC 5 μm C18 trap column (180 μm × 20 mm) at 5 μL/min during 3 min and then onto the nanoAcquity HSS T3 1.8 μm analytical reversed phase column (75 μm × 150 mm) at 400 nL/min, with a column temperature of 45°C. For peptide elution, a binary gradient was used, with mobile phase A consisting of water (Tedia, Fairfield, OH, United States) and 0.1% formic acid (Sigma–Aldrich), and mobile phase B consisting of acetonitrile (Sigma–Aldrich) and 0.1% formic acid. Gradient elution started at 7% B, then ramped from 7% B to 40% B up to 91.12 min, and from 40% B to 99.9% B until 92.72 min, being maintained at 99.9% until 106.00 min, then decreasing to 7% B until 106.1 min and kept 7% B until the end of the experiment at 120.00 min. Mass spectrometry was performed in positive and resolution mode (V mode), 35,000 FWHM, with ion mobility, and in data-independent acquisition (DIA) mode; ion mobility separation (HDMSE) using IMS wave velocity of 600 m/s, and helium and IMS gas flow of 180 and 90 mL/min, respectively; the transfer collision energy ramped from 19 V to 55 V in high-energy mode; cone and capillary voltages of 30 and 2750 V, respectively; and a source temperature of 70°C. In TOF parameters, the scan time was set to 0.5 s in continuum mode with a mass range of 50–2,000 Da. The human [Glu1]-fibrinopeptide B (Sigma–Aldrich) at 100 fmol/μL was used as an external calibrant and lock mass acquisition was performed every 30 s. Mass spectra acquisition was performed by MassLynx v4.0 software.
Proteomics Data Analysis
Spectral processing and database searching conditions were performed by ProteinLynx global server (PLGS; version 3.0.2) (Waters) and ISOQuant workflow software (Distler et al., 2014). The analysis used the following parameters: Apex3D of 150 counts for low energy threshold, 50 counts for elevated energy threshold, and 750 counts for intensity threshold; one missed cleavage, minimum fragment ion per peptide equal to three, minimum fragment ion per protein equal to seven, minimum peptide per protein equal to two, fixed modifications of carbamidomethyl (C), and variable modifications of oxidation (M) and phosphoryl (STY). The false discovery rate (FDR) for peptide and protein was set to a maximum of 1%, with a minimum peptide length of six amino acids. The analysis used the O. sativa ssp. japonica reference proteome from Uniprot2 (48,903 protein sequences downloaded on 18 August 2018).
Comparative label-free quantification analysis was performed using ISOQuant software (Distler et al., 2014). Briefly, the analysis included retention time alignment, exact mass retention time (EMRT) and IMS clustering, as well as data normalization and protein homology filtering. ISOQuant annotates the resulting feature clusters by evaluating consensus peptide identifications and identification probabilities. Protein identification parameters in ISOQuant were set to FDR 1%, peptide score greater than six, minimum peptide length of six amino acids, and at least two peptides per protein. Label-free quantification was estimated using the TOP3 quantification approach (Silva et al., 2006), followed by the multidimensional normalization process implemented within ISOQuant (Distler et al., 2014). After data processing and to ensure the quality of results, only proteins present or absent (for unique proteins) in three out of three runs were accepted and submitted to differential abundance analysis. Data were analyzed using Student’s t-test (two-tailed). Proteins with P < 0.05 were considered more abundant if the Log2 of fold change (FC) was greater than 0.5 and less abundant if the Log2 of FC was less than −0.5. All proteomics data were deposited to the ProteomeXchange Consortium (Deutsch et al., 2017) via the PRIDE repository3 (Perez-Riverol et al., 2019) and can be accessed using the code PXD020940.
Except for proteomics, data from other experiments were analyzed using the Student’s t-test (P ≤ 0.05 and 0.01) or one-way ANOVA followed by Tukey test (P ≤ 0.05), using SPSS Base 23.0 for Windows (SPSS Inc., United States).
Different Responses of Rice Cultivars/Species to S. oryzae Infestation
Unexpectedly, we verified that African rice species (O. barthii and O. glaberrima) and weedy rice (O. sativa f. spontanea) presented high susceptibility to S. oryzae infestation, while O. sativa cv. Nipponbare showed low levels of leaf damage (Figure 1). After 4 weeks of infestation, leaves of O. barthii were dry and the leaf area was totally injured. The species O. glaberrima presented infestation level 4 (above 80% of the damaged leaf area), while Nipponbare presented level 2 (10–20% of damaged leaf area), and weedy rice revealed level 3 (40–50% of the damaged leaf area), as shown in Figure 1. Both African species died before reaching the reproductive stage, while both O. sativa were able to complete their life cycle and set seeds. Although the plant height was not negatively affected by mite infestation in any of the tested cultivars/species, tiller number decreased under infested conditions in all plants, except in Nipponbare (Figure 2).
Figure 1. Visual characteristics of leaves from control and late infested (LI) plants of Oryza barthii, Oryza glaberrima, Oryza sativa cv. Nipponbare, and weedy rice (Oryza sativa f. spontanea). Bars indicate 0.5 cm.
Figure 2. Plant height (m) and tiller number of Oryza barthii (A), Oryza glaberrima (B), Oryza sativa cv. Nipponbare (C), and weedy rice (Oryza sativa f. spontanea) (D) at the vegetative stage (no infestation, V10-13) and full maturity stage (control or infested conditions, R9). Represented values are the averages of 50 samples ± SE. Different letters indicate that the means are different by the Tukey HSD test (P ≤ 0.05).
Chlorophyll a fluorescence analysis showed that African rice species (O. barthii and O. glaberrima) and weedy rice presented a decrease in at least one of the OJIP curve-times. O. barthii seems to be the most affected species, with a decrease in all curve-times after four weeks of infestation (LI leaves), while O. glaberrima presented a decrease in J-I-P stages, and weedy rice in I-P stages. Nipponbare was the only one with no decrease during the entire OJIP curve under infested condition (Figure 3), suggesting that mite infestation triggers lower damage in Nipponbare photosynthetic apparatus than in the other tested cultivars/species. These results agree with total chlorophyll concentration analysis, where Nipponbare was the only tested cultivars/species with no decrease under infested conditions, while O. barthii presented the highest decrease comparing infested and control conditions (Figure 4). Also, histochemical staining assay using DAB indicated that Nipponbare leaves accumulate lower levels of H2O2 than the other tested cultivars/species, especially when compared to O. barthii (Figure 5). Therefore, S. oryzae infestation differentially affects the generation of oxidative stress in LI leaves of the tested plants. Altogether, we suggest that O. barthii is extremely susceptible to S. oryzae infestation, while Nipponbare can be considered tolerant.
Figure 3. OJIP-test parameters calculated from the chlorophyll a fluorescence transient in control (C), early infested (EI), and late infested (LI) leaves of Oryza barthii, Oryza glaberrima, Oryza sativa cv. Nipponbare, and weedy rice (Oryza sativa f. spontanea). (A) O; (B) J; (C) I; (D) P. Represented values are the averages of 10 samples ± SE. Mean values (from each species/cultivar and each exposure time: C × EI; C × LI) with one, two, or three asterisks are significantly different as determined by a Student’s t-test (P ≤ 0.05, 0.01, and 0.001, respectively). r.u., Raman units.
Figure 4. Total chlorophyll concentration in control and late infested (LI) leaves of Oryza barthii, Oryza glaberrima, Oryza sativa cv. Nipponbare, and weedy rice (Oryza sativa f. spontanea) plants. Represented values are the averages of three samples ± SE. Mean values (from each species/cultivar: C × LI) with one asterisk are significantly different as determined by a Student’s t-test (p-value ≤ 0.05). DW, dry weight.
Figure 5. Histochemical staining assay of H2O2 by diaminobenzidine (DAB) in control and late infested (LI) leaves of Oryza barthii, Oryza glaberrima, Oryza sativa cv. Nipponbare, and weedy rice (Oryza sativa f. spontanea). The positive staining (detected in higher levels on infested leaves) in the photomicrographs shows as bright images (brown-color). Bars indicate 0.5 cm.
As African rice species were not able to reach the reproductive stage, seeds from Nipponbare and weedy rice were evaluated to verify whether S. oryzae infestation could impact seed production. Weedy rice presents higher number of panicles per plant than Nipponbare under both control and infested conditions (Figure 6A), even though the number of panicles per plant was affected by S. oryzae infestation only in weedy rice (Figure 6A). In Nipponbare, none of the analyzed parameters (number of panicles per plant— Figure 6A, number of seeds per panicle—Figure 6B, percentage of full seeds per panicle—Figure 6C, grain length—Figures 6E,F, and weight of 1,000 full grains—Figure 6G) was negatively affected by mite infestation. As seen in Figure 6C, 93% of the seeds produced by Nipponbare plants under infested condition were full, while in the same condition, weedy rice produced 87% of full seeds. Such difference, together with the higher number of panicles per plant in weedy rice than Nipponbare under infested condition (Figure 6A), probably explains the high number of empty seeds in weedy rice plants infested by S. oryzae (Figure 6D). Surprisingly, two parameters (percentage of full seeds per panicle—Figure 6C and grain length—Figures 6E,F) presented higher values in Nipponbare plants under infested than control condition, resulting in the maintenance of seed weight per plant (an estimate of yield) under infested condition (Figure 7). Therefore, Nipponbare reveals physiological characteristics of tolerance, since even during infestation, plants can maintain development, photosynthetic/antioxidant activities, and grain yield.
Figure 6. Seeds analysis of Oryza sativa cv. Nipponbare and weedy rice (Oryza sativa f. spontanea). (A) Number of panicles per plant. (B) Number of seeds (empty + full) per panicle. (C) Percentage of full seeds per panicle. (D) Number of empty seeds (10 plants) under infested condition. (E) Grain length (mm). (F) Grain visual characteristics. (G) Weight of 1,000 full grains. Represented values are the averages of 50 samples ± SE. Mean values (from each cultivar: Control × Infested) with one or three asterisks are significantly different as determined by a Student’s t-test (P ≤ 0.05 and 0.001).
Figure 7. Seed weight per plant (estimative of yield) of Oryza sativa cv. Nipponbare under control and infested conditions. Represented values are the averages of 10 samples ± SE.
Overview of Proteomic Analysis
A total of 1,041 proteins were identified comparing control and infested conditions in O. barthii and Nipponbare, with 195 (18.7%) proteins unique to or differentially abundant between treatments. As seen in Figure 8, comparing control and infested leaves of O. barthii, we detected 48 differentially abundant proteins, being two more abundant (and one unique) in control condition, and 44 more abundant (and one unique) in infested condition. We identified 147 differentially abundant proteins in control and infested conditions of Nipponbare, being 35 more abundant (and two unique) in control condition, and 104 more abundant (and six unique) in infested condition. Due to the different genetic background between O. barthii and O. sativa, we prefer not to compare the same conditions between species, in order to avoid the emergence of species-specific sequences that are not mite related.
Figure 8. Venn diagram showing the overlap of rice proteins identified in control and early infested (EI) leaves of mite susceptible Oryza barthii and mite-tolerant Oryza sativa cv. Nipponbare. Dark green circles: control leaves; yellow circles: EI leaves. Light green means overlap.
The corresponding sequence of each identified protein was compared to NCBI BLASTp to identify specific domains, molecular functions, and protein annotations. Afterward, proteins were arbitrarily categorized in functional categories, according to the available literature and their putative molecular function. The lists of all unique or differentially abundant proteins identified in this work are presented in Supplementary Tables 1, 2. In order to facilitate the overall understanding of the differential protein abundance in each functional category presented in Supplementary Tables 1, 2, we summarize these data in Table 1. We detected four different patterns, based on the number of differentially abundant proteins: (1) most of the functional categories (antioxidant system, carbohydrate metabolism/energy production, general metabolic processes, hormone-related, lipid metabolism, protein modification/degradation, stress response, and others) were more represented in infested conditions, regardless the analyzed species; (2) four categories were differentially represented only in one species (amino acid metabolism in O. barthii under infested condition; cytoskeleton and transport in Nipponbare under control condition; and protease inhibitor (PI) also in Nipponbare, but under infested condition); (3) photosynthesis (over represented in O. barthii under infested condition, and in Nipponbare under control condition); and (4) translation, with no difference between control and infested conditions, regardless of the analyzed species.
Table 1. Schematic representation of the functional categories differently represented in leaves of mite susceptible Oryza barthii and mite-tolerant Oryza sativa cv. Nipponbare plants under control and infested conditions.
Considering the tolerant Nipponbare, S. oryzae infestation seems to be less damaging and to generate a more complex defense response. It is interesting to highlight the higher diversity and expression level of antioxidant proteins in Nipponbare than O. barthii under infested conditions, which agrees with the lower oxidative stress seen in Nipponbare leaves (Figure 5), and also the higher diversity of proteins involved in general metabolic processes and carbohydrate metabolism/energy production, suggesting that Nipponbare can maintain the basal/primary metabolism more active than O. barthii under infested condition. This finding also agrees with our observation that photosynthesis is maintained in Nipponbare, whereas it is drastically affected in O. barthii. At the same time, Nipponbare seems more able to fight S. oryzae infestation than O. barthii, presenting higher diversity and/or representation level of PIs and proteins involved in stress response. As seen in Figure 9, we also detected enhanced proline levels only in Nipponbare leaves under infested condition, suggesting an active defense mechanism. Several proteins belonging to these functional categories are discussed in more detail below.
Figure 9. Proline levels in control and late infested (LI) leaves of mite susceptible Oryza barthii and mite-tolerant Oryza sativa cv. Nipponbare. Represented values are the averages of three samples ± SE. Mean values with one asterisk are significantly different as determined by a Student’s t-test (P ≤ 0.05). FW, fresh weight.
Our first screening of Oryza cultivars/species responses to S. oryzae infestation showed that none presents a classical resistance response involving antibiosis and antixenosis (Stenberg and Muola, 2017), since mite population increased over time for all cultivars/species, even though with different kinetics (Figure 1). Physiological analysis and agronomical parameters showed that O. glaberrima and especially O. barthii are extremely susceptible to S. oryzae infestation (being unable to reach the reproductive stage), while weedy rice and especially Nipponbare are less affected by mite infestation (Figures 1–7). Recently, we detected an indica rice cultivar able to maintain seed yield under S. oryzae infestation condition, which we characterized as mite tolerant (Buffon et al., 2018). Even though wild plant species have been widely recognized as valuable source of resistance genes for developing herbivore-resistant cultivars, Veasey et al. (2008) and Chandrasena et al. (2016) found no signs of mite resistance when different wild rice species were, respectively, tested against the infestation of S. oryzae and Steneotarsonemus spinki. These data differ from previously described ones, where O. barthii and O. glaberrima are recognized as resistant to diseases and pests (Khush, 1997; Linares, 2002). Still, our report is one of the first to test mite tolerance in wild and cultivated rice species.
Grain size is an important component trait of grain yield, which is frequently threatened by stress conditions (Lv et al., 2019; Dong et al., 2020). Interestingly, in our analysis, Nipponbare plants under infested condition presented higher seed length when compared to control condition (Figures 6E,F). Several studies have identified key grain size regulator genes, involved in different mechanisms (reviewed by Li et al., 2018). Even though none of the 43 genes cited by Li et al. (2018) as involved in rice grain size control code proteins differentially abundant in our proteomic analysis, it would be interesting to go deeper in such question, since we found similar results in a previous study with another mite-tolerant indica rice cultivar (IRGA 423—Buffon et al., 2018). Recent studies have identified molecular links between regulation of grain size and abiotic stress tolerance (Chen et al., 2016, 2019; Lv et al., 2019; Dong et al., 2020). However, as far as we know, any link between grain size regulation and biotic stress tolerance is unknown.
Recently, Sperotto et al. (2018b) hypothesized that mite-sensitivity presented by wild rice species could be explained, at least partially, by a presumable high GA:JA ratio in these plants, since most of the wild rice species are tall plants and probably present high GA levels. The roles of GA-JA in growth-defense conflicts during herbivory are not yet fully characterized. However, it is already known that plants need to prioritize GA- or JA-induced responses (Yang et al., 2019). Even though proteins related to GA biosynthesis have not been detected as differentially abundant in O. barthii or Nipponbare under infested conditions (Supplementary Tables 1, 2), the height of these plants at full maturity stage is quite different (55% higher in O. barthii—Figure 2). It is important to highlight that O. barthii and Nipponbare plant heights are not affected by S. oryzae infestation (Figure 2), and this could explain the lack of differentially abundant proteins involved in GA biosynthesis when we compare control and infested conditions. Sperotto et al. (2018b) suggest short Oryza species (O. minuta, O. meyeriana, O. neocaledonica, and O. schlechteri) as primary sources of herbivory tolerance, including mites. According to these authors, plants with low GA level/sensitivity would probably amplify JA responses and drive plant resources toward defense instead of growth. Based on our data, the GA-based mechanism is not involved in the different tolerance level found in the tested species, although JA might play a role.
One of the most abundant proteins found in Nipponbare under infested condition is lipoxygenase (log2 FC 4.7), which is well known as a precursor of plant defense mechanisms. Once the herbivore feeds on plant tissues, polyunsaturated fatty acids are released from cell membranes, and accumulate at the wound site (Roach et al., 2015). Also, polyunsaturated fatty acids act as a substrate to lipoxygenase, leading to the production of several defense-related compounds, including different structures of JAs such as methyl jasmonate (MeJA) (Mochizuki et al., 2016; Yang et al., 2019). We detected a key enzyme related to JA biosynthesis, allene oxide cyclase (AOC), which converts allene oxide to 12-oxophytodienoic acid (Nguyen et al., 2020), more abundant in both species under infested conditions. However, the log2 FC is higher in Nipponbare than O. barthii (Supplementary Tables 1, 2). Therefore, it is possible that Nipponbare presents a more active JA biosynthesis/signaling pathway, which can probably contribute to mite tolerance. Interestingly, MYB and WRKY transcription factors, which are involved in the activation of defense responses in plants, including JA signaling (Wang et al., 2014; Zhou and Memelink, 2016), were not detected as differentially abundant in our work, raising the possibility of new signaling/defense pathways in O. barthii and Nipponbare infested by S. oryzae.
The functional category most affected by S. oryzae infestation was antioxidant system, with 10 proteins more abundant in the mite susceptible O. barthii (Supplementary Table 1), and 22 in the mite-tolerant Nipponbare (Supplementary Table 2) under infested condition. This indicates that both species try to defend and maintain its redox homeostasis, since it is governed by the presence of large antioxidant pools that absorb and buffer reductants and oxidants (Foyer and Noctor, 2005). Considering the leaf damage (Figure 1), H2O2 accumulation (Figure 5), and the high diversity of differentially abundant proteins (Supplementary Tables 1, 2), such strategy seems to be more efficient in the mite-tolerant Nipponbare than in O. barthii. Six proteins identified as more abundant under infested condition only in Nipponbare are ferredoxin, ferredoxin-NADP+ reductase, thioredoxin H-type, peroxiredoxin, glutathione reductase, and superoxide dismutase (Supplementary Table 2). The stromal ferredoxin–thioredoxin system, a regulatory mechanism linking light to the activity of associated enzymes, is the best-characterized redox signal transduction system in plants, functioning in the regulation of photosynthetic carbon metabolism (Foyer and Noctor, 2005; Schürmann and Buchanan, 2008). Ferredoxin-NADP+ reductase catalyzes the reduction of NADP+ to NADPH, using the electrons provided by reduced ferredoxin (Cassan et al., 2005). Both NADPH and reduced ferredoxin are required for reductive assimilation and light/dark activation/deactivation of enzymes. Ferredoxin-NADP+ reductase is therefore a hub, connecting photosynthetic electron transport to chloroplast redox metabolism (Kozuleva et al., 2016). Even though the precise mechanism remains unknown, a clear correlation between ferredoxin-NADP+ reductase content and tolerance to oxidative stress is well established (Giró et al., 2006; Kozuleva et al., 2016).
Thioredoxins H-type are involved in the cellular protection against oxidative stress (Gelhaye et al., 2004), and function as electron donors to several antioxidant enzymes, such as peroxiredoxin and glutathione reductase (Rouhier et al., 2001). Plants over-expressing thioredoxin H-type genes show oxidative stress tolerance and lower levels of H2O2 accumulation (Sun et al., 2010; Zhang et al., 2011), similar to our findings with the mite-tolerant Nipponbare (Figure 5). Peroxiredoxins are abundant low-efficiency peroxidases that act as antioxidant and are involved in modulating redox-dependent signaling cascades (Dietz, 2003). Also, after being reduced by thioredoxin, peroxiredoxins reductively convert H2O2 to H2O at the chloroplast (König et al., 2002), acting together with other peroxidase enzymes that use ascorbate as the electron donor. Ascorbate can be directly regenerated using photosynthetic electrons via reduced ferredoxin or using reduced glutathione (GSH). In the chloroplast, the oxidized glutathione (GSSG) is reduced by the glutathione reductase enzyme using NADPH (Kozuleva et al., 2016). It is important to highlight that very efficient scavenging mechanisms exist in chloroplasts to prevent oxidative damage, including the rapid conversion of superoxide radical (O2•–, one of the dominant species produced during light excitation) to H2O2 by superoxide dismutase enzymes (Mubarakshina and Ivanov, 2010). This is the most likely explanation for the lack of differences in O2•– accumulation in the leaves of plants subjected to control and infested conditions (data not shown). Altogether, these data show that the mite-tolerant Nipponbare under infested condition can maintain the correct balance between rapid removal of damaging oxidative species and appropriate concentrations necessary to initiate signaling cascades.
We verified that all tested cultivars/species revealed a decrease in chlorophyll fluorescence parameters and total chlorophyll concentration, except the mite-tolerant Nipponbare (Figures 3, 4). Increased photosynthesis is an important mechanism behind tolerance to herbivory (Thomson et al., 2003), and parameters of chlorophyll fluorescence (excess energy absorbed by chlorophyll that is reemitted as light) have been used to obtain quantitative information on the photosynthetic performance of plants (Maxwell and Johnson, 2000), since chlorophyll fluorescence curves allow to evaluate the physiological condition of photosystem II (PSII) components and photosynthetic electron transport chain (Kalaji et al., 2016), where ferredoxin protein acts. Recently, Madriaza et al. (2019) indicated that chlorophyll fluorescence is an accurate predictor of tolerance to herbivore damage, and claimed that this rapid methodology could be used to estimate plant tolerance to herbivory in natural populations.
Even though plant height was not affected by S. oryzae infestation in O. barthii and Nipponbare, the mite susceptible O. barthii presented a decrease in tiller number under infested condition, while Nipponbare maintained the same level found in control condition (Figure 2), resulting in no yield losses (Figure 7). Reduction in tillering (and consequently in grain weight) was already seen in rice plants infested by the rice water weevil, Lissorhoptrus oryzophilus Kuschel (Zou et al., 2004), and also in IR22 (a rice variety susceptible to the brown planthopper, Nilaparvata lugens) (Horgan et al., 2018). It is important to highlight that planthopper damage to IR22 in field cages was severe, and plant death occurred in most of IR22 plants, similar to our findings with O. barthii and O. glaberrima under infested conditions. Both tillering and seed filling processes are expensive for the crop energy balance. As seen in Supplementary Tables 1, 2, we verified a larger diversity of proteins involved in general metabolic processes and carbohydrate metabolism/energy production under infested condition in tolerant Nipponbare than in susceptible O. barthii, suggesting that general metabolism and energy production in Nipponbare are less affected by S. oryzae infestation.
Two proteins related to general metabolic processes were identified as more abundant under infested condition only in Nipponbare: CBS domain-containing protein (CBSX1) and cysteine synthase (Supplementary Table 2). In Arabidopsis, AtCBSX1 and AtCBX2 positively regulate thioredoxins, which then reduce target proteins, such as peroxiredoxin, which in turn reduce cellular H2O2 and regulate its level (Yoo et al., 2011). Recently, Kumar et al. (2018) overexpressed a rice Two Cystathionine-β-Synthase Domain-containing Protein (OsCBSCBSPB4) in tobacco. Transgenic seedlings exhibit delayed leaf senescence, profuse root growth, and increased biomass in contrast to the wild-type seedlings when subjected to salinity, dehydration, oxidative, and extreme temperature treatments, suggesting that OsCBSCBSPB4 is involved in abiotic stress response. As far as we know, this is the first time that a CBS protein is related to biotic stress tolerance in plants. The biosynthesis of cysteine is a limiting step in the production of glutathione, a thiol implicated in various cellular functions, including scavenging of reactive oxygen species and resistance to biotic stresses (Youssefian et al., 2001). At the same time, several studies have pointed out the function of cysteine-thiols and sulfur-based primary and secondary metabolites in basal defense responses to abiotic and biotic stresses (for a comprehensive review, see Tahir and Dijkwel, 2016). Even though herbivory resistance can employ metabolites containing cysteine such as glucosinolates (Guo et al., 2013), to the best of our knowledge, this is the first work that links the expression of cysteine synthase protein with mite tolerance.
Two proteins related to stress response were also identified as more abundant under infested condition only in Nipponbare: osmotin-like protein and ricin B-like lectin. Osmotin protein belongs to the PR-5 family of pathogenesis-related (PR) proteins, which are produced in response to various biotic and abiotic stresses. Osmotin can increase its own cytotoxic efficiency promoting changes in the cell wall structure and inhibiting defensive cell wall barriers, which in turn enables osmotin penetration into the plasma membrane (Hakim et al., 2018). Overexpression of genes encoding osmotins confers tolerance against biotic and abiotic stresses in different plant species (Barthakur et al., 2001; Xue et al., 2016; Chowdhury et al., 2017), besides stimulating the accumulation of proline, which functions as a compatible osmolyte useful in responding to abiotic stress, especially drought and salinity (Hakim et al., 2018). Proline also functions suppressing free radicals and preventing chlorophyll degradation in vegetative cells, in order to maintain the photosynthetic processes (Barthakur et al., 2001; Chowdhury et al., 2017). In the mite-tolerant Nipponbare, we detected higher levels of proline accumulation under infested condition when compared to control (Figure 9), along with maintenance of chlorophyll concentration (Figure 4) after mite infestation. This is the first time that expression of an osmotin protein is related to mite tolerance. Since the generation of transgenic plants overexpressing osmotin genes is being recommended as a biotechnological tool to generate stress tolerance, we are currently overexpressing and knocking out OsOSM1 gene (LOC_Os12g38170) in rice plants.
Ricin B-like lectins, one of the most widespread families of carbohydrate-binding proteins (Lannoo and Van Damme, 2014), have been suggested to play a role in plant defense against pathogens (Vandenbussche et al., 2004) and insects (Wei et al., 2004; Shahidi-Noghabi et al., 2009). Since mutation of the carbohydrate-binding site can abolish or reduce the toxic effect, the entomotoxic properties of the proteins can be linked to their carbohydrate-binding activity (Shahidi-Noghabi et al., 2008). Our work is the first one that correlates a ricin B-like lectin expression with mite response/tolerance in plants. Curiously, it was shown that β-trefoil structure enables interactions between fungal ricin B-like lectins and PIs, and that these interactions modulate their biological activity (Žurga et al., 2015). It would be interesting to test whether the two PIs (putative Bowman Birk trypsin inhibitor) we found more abundant only in the mite-tolerant Nipponbare under infested condition (Supplementary Table 2) are able to interact with the ricin B-like lectins.
Among plant defenses, PIs have direct effects on herbivores, interfering with their physiology (Diaz and Santamaria, 2012; Martinez et al., 2016; Arnaiz et al., 2018). Once ingested by the arthropod, PIs inhibit proteolytic activities, preventing or hindering protein degradation (Terra and Ferreira, 2005). The first successfully overexpressed Bowman-Birk family PI gene in tobacco was the cowpea trypsin inhibitor (CpTI—Hilder et al., 1987). After, the CpTI gene was inserted into the genome of several plants, increasing the resistance to different herbivore species (Santamaria et al., 2012). Since then, different PIs have been obtained from cultivated plants such as rice, barley, soybean, sweet potato, and corn, and have been overexpressed in different plant species, conferring resistance to various pest insects (Clemente et al., 2019). Rice has seven Bowman Birk trypsin inhibitor genes, and overexpression of RBBI2-3 in transgenic rice plants resulted in resistance to the fungal pathogen Pyricularia oryzae, indicating that PIs confer in vivo resistance against the fungal pathogen, and might play a role in the defense system of the rice plant (Qu et al., 2003). Recently, Arnaiz et al. (2018) verified that Kunitz trypsin inhibitors (KTI) participate in the plant defense responses against the spider mite Tetranychus urticae. Thus, our work is the first one to suggest that Bowman Birk trypsin inhibitors may be involved in the tolerance of rice plants to S. oryzae mite infestation.
The model in Figure 10 summarizes the rice tolerance mechanisms putatively employed by Nipponbare plants.
Figure 10. Rice mechanisms employed by Oryza sativa cv. Nipponbare to tolerate Schizotetranychus oryzae mite infestation.
Tolerance is a sustainable pest management strategy as it only involves plant response, and therefore does not cause evolution of resistance in target pest populations (Peterson et al., 2017; Sperotto et al., 2018a). According to Gerszberg and Hnatuszko-Konka (2017), elucidation of the tolerance mechanisms at the biochemical, physiological, and morphological level remains one of the greatest challenges of contemporary plant physiology. Our physiological results showed that African rice species (wild O. barthii and cultivated O. glaberrima) and weedy rice (O. sativa f. spontanea) do not present tolerance to S. oryzae mite infestation. Still, we characterized O. sativa cv. Nipponbare as tolerant to mite infestation. Proteomic data showed that O. barthii presents a less diverse/efficient antioxidant apparatus under infested condition, not being able to modulate proteins involved in general metabolic processes and energy production, resulting in damage to its photosynthetic apparatus and development. On the other hand, Nipponbare presents more efficient antioxidant mechanisms that maintain cellular homeostasis and its photosynthetic capacity. Tolerant plants are able to modulate proteins involved in general metabolic processes and energy production, maintaining their development. In addition, these plants can modulate the expression of defense proteins, such as osmotin, ricin B-like lectin, and PIs, which may be used in future breeding programs for increasing rice tolerance to mite infestation.
Data Availability Statement
All proteomics data were deposited to the ProteomeXchange Consortium (Deutsch et al., 2017) via the PRIDE repository (https://www.ebi.ac.uk/pride/; Perez-Riverol et al., 2019) and can be accessed using the code PXD020940.
JS, FR, and RS conceived and designed the research. GB, ÉB, TL, JA, and AB conducted the experiments. VS and ML contributed with analytical tools. GB, ÉB, JA, AB, VS, and RS analyzed the data. GB and RS wrote the manuscript. All authors read and approved the manuscript.
This research was supported by University of Taquari Valley–Univates, Fundação de Amparo à Pesquisa do Estado do RS (FAPERGS) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Instituto Rio-Grandense do Arroz (IRGA) for technical support. Pre-print of this manuscript at bioRxiv: https://doi.org/10.1101/2020.01.22.914184.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.613568/full#supplementary-material
Supplementary Table 1 | Differentially abundant proteins in mite susceptible Oryza barthii (control × early infested condition).
Supplementary Table 2 | Differentially abundant proteins in mite-tolerant Oryza sativa cv. Nipponbare (control × early infested condition).
DAB, diaminobenzidine; EI, early infested; GA, gibberellin; JA, jasmonic acid; LI, late infested; PI, protease inhibitor.
- ^ https://www.rgrc.dna.affrc.go.jp/
- ^ https://www.uniprot.org/proteomes/UP000059680
- ^ https://www.ebi.ac.uk/pride/
Arnaiz, A., Talavera-Mateo, L., Gonzalez-Melendi, P., Martinez, M., Diaz, I., and Santamaria, M. E. (2018). Arabidopsis Kunitz trypsin inhibitors in defense against spider mites. Front. Plant Sci. 9:986. doi: 10.3389/fpls.2018.00986
Barthakur, S., Babu, V., and Bansa, K. C. (2001). Over-expression of osmotin induces proline accumulation and confers tolerance to osmotic stress in transgenic tobacco. J. Plant Biochem. Biotechnol. 10, 31–37. doi: 10.1007/bf03263103
Blasi, E. A. R., Buffon, G., da Silva, R. Z., Stein, C., Dametto, A., Ferla, N. J., et al. (2015). Alterations in rice, corn and wheat plants infested by phytophagous mite. Int. J. Acarol. 41, 10–18. doi: 10.1080/01647954.2014.988643
Brar, D. S., and Khush, G. S. (2003). “Utilization of wild species of genus Oryza in rice improvement,” in Monograph on Genus Oryza, eds J. S. Nanda and S. D. Sharma (Enfield, NH: Science Publishers.), 283–309.
Buffon, G., Blasi, ÉA. R., Rativa, A. G. S., Lamb, T. I., Gastmann, R., Adamski, J. M., et al. (2018). Unraveling rice tolerance mechanisms against Schizotetranychus oryzae mite infestation. Front. Plant Sci. 9:1341. doi: 10.3389/fpls.2018.01341
Burrieza, H. P., Rizzo, A. J., Moura Vale, E., Silveira, V., and Maldonado, S. (2019). Shotgun proteomic analysis of quinoa seeds reveals novel lysine-rich seed storage globulins. Food Chem. 293, 299–306. doi: 10.1016/j.foodchem.2019.04.098
Cassan, N., Lagoutte, B., and Sétif, P. (2005). Ferredoxin-NADP+ reductase - Kinetics of electron transfer, transient intermediates, and catalytic activities studied by flash-absorption spectroscopy with isolated photosystem I and ferredoxin. J. Biol. Chem. 280, 25960–25972.
Chandrasena, G. D. S. N., Jayawardane, J. D. K. M., Umange, S. D., and Gunawardana, A. D. B. U. (2016). Host range of panicle rice mite Steneotarsonemus spinki smiley (Acari: Tarsonemidae) in Sri Lanka. Univ. J. Agric. Res. 4, 21–24. doi: 10.13189/ujar.2016.040104
Chen, C., Begcy, K., Liu, K., Folsom, J. J., Wang, Z., Zhang, C., et al. (2016). Heat stress yields a unique MADS box transcription factor in determining seed size and thermal sensitivity. Plant Physiol. 171, 606–622. doi: 10.1104/pp.15.01992
Chen, S., Zhang, N., Zhang, Q., Zhou, G., Tian, H., Hussain, S., et al. (2019). Genome editing to integrate seed size and abiotic stress tolerance traits in Arabidopsis reveals a role for DPA4 and SOD7 in the regulation of inflorescence architecture. Int. J. Mol. Sci. 20:2695. doi: 10.3390/ijms20112695
Chowdhury, S., Basu, A., and Kundu, S. (2017). Overexpression of a new osmotin-like protein gene (SindOLP) confers tolerance against biotic and abiotic stresses in sesame. Front. Plant Sci. 8:410. doi: 10.3389/fpls.2017.00410
Clemente, M., Corigliano, M. G., Pariani, S., Sánchez-López, E. F., Sander, V. A., and Ramos-Duarte, V. A. (2019). Plant serine protease inhibitors: biotechnology application in agriculture and molecular farming. Int. J. Mol. Sci. 20:E1345.
Counce, P. A., Keisling, T. C., and Mitchell, A. J. (2000). A uniform, objective, and adaptive system for expressing rice development. Crop Sci. 40, 436–443. doi: 10.2135/cropsci2000.402436x
Culliney, T. (2014). “Crop losses to arthropods,” in Integrated Pest Management, eds D. Pimentel and R. Peshin (Dordrecht: Springer), 201–225. doi: 10.1007/978-94-007-7796-5_8
Dai, L., Dai, W., Song, X., Lu, B., and Qiang, S. (2014). A comparative study of competitiveness between different genotypes of weedy rice (Oryza sativa) and cultivated rice. Pest Manag. Sci. 70, 113–122. doi: 10.1002/ps.3534
Deutsch, E. W., Csordas, A., Sun, Z., Jarnuczak, A., Perez-Riverol, Y., Ternent, T., et al. (2017). The ProteomeXchange consortium in 2017: supporting the cultural change in proteomics public data deposition. Nucleic Acids Res. 45, D1100–D1106.
Diaz, I., and Santamaria, M. E. (2012). “Biotechnological approaches to combat phytophagous arthropods,” in Arthropod-Plant Interactions. Progress in Biological Control, eds G. Smagghe and I. Diaz (Dordrecht: Springer), 159–176. doi: 10.1007/978-94-007-3873-7_6
Dietz, K. J. (2003). Plant peroxiredoxins. Annu. Rev. Plant Biol. 54, 93–107.
Distler, U., Kuharev, J., Navarro, P., Levin, Y., Schild, H., and Tenzer, S. (2014). Drift time-specific collision energies enable deep-coverage data-independent acquisition proteomics. Nat. Methods 11, 167–170. doi: 10.1038/nmeth.2767
Dong, N. Q., Sun, Y., Guo, T., Shi, C. L., Zhang, Y. M., Kan, Y., et al. (2020). UDP-glucosyltransferase regulates grain size and abiotic stress tolerance associated with metabolic flux redirection in rice. Nat. Commun. 11:2629.
Foyer, C. H., and Noctor, G. (2005). Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell 17, 1866–1875. doi: 10.1105/tpc.105.033589
Gelhaye, E., Rouhier, N., and Jacquot, J. P. (2004). The thioredoxin h system of higher plants. Plant Physiol. Biochem. 42, 265–271. doi: 10.1016/j.plaphy.2004.03.002
Gerszberg, A., and Hnatuszko-Konka, K. (2017). Tomato tolerance to abiotic stress: a review of most often engineered target sequences. Plant Growth Regul. 83, 175–198. doi: 10.1007/s10725-017-0251-x
Giró, M., Carrillo, N., and Krapp, A. R. (2006). Glucose-6-phosphate dehydrogenase and ferredoxin-NADP(H) reductase contribute to damage repair during the soxRS response of Escherichia coli. Microbiology 152, 1119–1128. doi: 10.1099/mic.0.28612-0
Godfray, H. C., Beddington, J. R., Crute, I. R., Haddad, L., Lawrence, D., Muir, J. F., et al. (2010). Food security: the challenge of feeding 9 billion people. Science 327, 812–818.
Guo, R., Shen, W., Qian, H., Zhang, M., Liu, L., and Wang, Q. (2013). Jasmonic acid and glucose synergistically modulate the accumulation of glucosinolates in Arabidopsis thaliana. J. Exp. Bot. 64, 5707–5719. doi: 10.1093/jxb/ert348
Hakim, Ullah, A., Hussain, A., Shaban, M., Khan, A. H., Alariqi, M., et al. (2018). Osmotin: a plant defense tool against biotic and abiotic stresses. Plant Physiol. Biochem. 123, 149–159. doi: 10.1016/j.plaphy.2017.12.012
Henry, R. J. (2014). Genomics strategies for germplasm characterization and the development of climate resilient crops. Front. Plant Sci. 5:68. doi: 10.3389/fpls.2014.00068
Hilder, V. A., Gatehouse, A. M. R., Sheerman, S. E., Barker, R. F., and Boulter, D. (1987). A novel mechanism of insect resistance engineered into tobacco. Nature 330, 160–163. doi: 10.1038/330160a0
Horgan, F. G., Peñalver Cruz, A., Bernal, C. C., Ramal, A. F., Almazan, M. L. P., and Wilby, A. (2018). Resistance and tolerance to the brown planthopper, Nilaparvata lugens (Stål), in rice infested at different growth stages across a gradient of nitrogen applications. Field Crops Res. 217, 53–65. doi: 10.1016/j.fcr.2017.12.008
Kalaji, H. M., Jajoo, A., Oukarroum, A., Brestic, M., Zivcak, M., Samborska, I. A., et al. (2016). Chlorophyll a fluorescence as a tool to monitor physiological status of plants under abiotic stress conditions. Acta Physiol. Plant. 38:102.
Khush, G. S. (1997). Origin, dispersal, cultivation and variation of rice. Plant Mol. Biol. 35, 25–34. doi: 10.1007/978-94-011-5794-0_3
König, J., Baier, M., Horling, F., Kahmann, U., Harris, G., Schürmann, P., et al. (2002). The plant-specific function of 2-Cys peroxiredoxin-mediated detoxification of peroxides in the redox-hierarchy of photosynthetic electron flux. Proc. Natl. Acad. Sci. U.S.A. 99, 5738–5743. doi: 10.1073/pnas.072644999
Kozuleva, M., Goss, T., Twachtmann, M., Rudi, K., Trapka, J., Selinski, J., et al. (2016). Ferredoxin: NADP (H) oxidoreductase abundance and location influences redox poise and stress tolerance. Plant Physiol. 172, 1480–1493. doi: 10.1104/pp.16.01084
Kumar, R., Subba, A., Kaur, C., Ariyadasa, T. U., Sharan, A., Pareek, A., et al. (2018). OsCBSCBSPB4 is a two cystathionine-β-synthase domain-containing protein from rice that functions in abiotic stress tolerance. Curr. Genom. 19, 50–59.
Lannoo, N., and Van Damme, E. J. M. (2014). Lectin domains at the frontiers of plant defense. Front. Plant Sci. 5:397. doi: 10.3389/fpls.2014.00397
Li, N., Xu, R., Duan, P., and Li, Y. (2018). Control of grain size in rice. Plant Reprod. 31, 237–251. doi: 10.1007/s00497-018-0333-6
Linares, O. F. (2002). African rice (Oryza glaberrima): history and future potential. Proc. Natl. Acad. Sci. U.S.A. 99, 16360–16365. doi: 10.1073/pnas.252604599
Lv, M. J., Wan, W., Yu, F., and Meng, L. S. (2019). New insights into the molecular mechanism underlying seed size control under drought stress. J. Agric. Food Chem. 67, 9697–9704. doi: 10.1021/acs.jafc.9b02497
Madriaza, K., Saldaña, A., Salgado-Luarte, C., Escobedo, V. M., and Gianoli, E. (2019). Chlorophyll fluorescence may predict tolerance to herbivory. Int. J. Plant Sci. 180, 81–85. doi: 10.1086/700583
Martinez, M., Santamaria, M. E., Diaz-Mendoza, M., Arnaiz, A., Carrillo, L., Ortego, F., et al. (2016). Phytocystatins: defense proteins against phytophagous insects and Acari. Int. J. Mol. Sci. 17, E1747.
Maxwell, K., and Johnson, G. N. (2000). Chlorophyll fluorescence: a practical guide. J. Exp. Bot. 51, 659–668. doi: 10.1093/jxb/51.345.659
Menguer, P. K., Sperotto, R. A., and Ricachenevsky, F. K. (2017). A walk on the wild side: Oryza species as source for rice abiotic stress tolerance. Genet. Mol. Biol. 40, 238–252. doi: 10.1590/1678-4685-gmb-2016-0093
Mochizuki, S., Sugimoto, K., Koeduka, T., and Matsui, K. (2016). Arabidopsis lipoxygenase 2 is essential for formation of green leaf volatiles and five−carbon volatiles. FEBS Lett. 590, 1017–1027. doi: 10.1002/1873-3468.12133
Mubarakshina, M. M., and Ivanov, B. N. (2010). The production and scavenging of reactive oxygen species in the plastoquinone pool of chloroplast thylakoid membranes. Physiol. Plant. 140, 103–110. doi: 10.1111/j.1399-3054.2010.01391.x
Muthayya, S., Sugimoto, J. D., Montgomery, S., and Maberly, G. F. (2014). An overview of global rice production, supply, trade, and consumption. Ann. N. Y. Acad. Sci. 1324, 7–14. doi: 10.1111/nyas.12540
Nanjo, Y., Skultety, L., Uvacìkova, L., Klubicova, K., Hajduch, M., and Komatsu, S. (2012). Mass spectrometry-based analysis of proteomic changes in the root tips of flooded soybean seedlings. J. Prot. Res. 11, 372–385. doi: 10.1021/pr200701y
Nguyen, T. H., Mai, H. T. T., Moukouanga, D., Lebrun, M., Bellafiore, S., and Champion, A. (2020). “CRISPR/Cas9-mediated gene editing of the jasmonate biosynthesis OsAOC gene in rice,” in Jasmonate in Plant Biology. Methods in Molecular Biology, eds A. Champion and L. Laplaze (New York, NY: Humana Press), 199–209. doi: 10.1007/978-1-0716-0142-6_15
Oerke, E. C. (2006). Crop losses to pests. J. Agric. Sci. 144, 31–43. doi: 10.1017/s0021859605005708
Palmgren, M. G., Edenbrandt, A. K., Vedel, S. E., Andersen, M. M., Landes, X., Østerberg, J. T., et al. (2015). Are we ready for back-to-nature crop breeding? Trends Plant Sci. 20, 155–164. doi: 10.1016/j.tplants.2014.11.003
Perez-Riverol, Y., Csordas, A., Bai, J., Bernal-Llinares, M., Hewapathirana, S., Kundu, D. J., et al. (2019). The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 47, D442–D450.
Peterson, R. K. D., Varella, A. C., and Higley, L. G. (2017). Tolerance: the forgotten child of plant resistance. PeerJ 5:e3934. doi: 10.7717/peerj.3934
Qiu, J., Zhou, Y., Mao, L., Ye, C., Wang, W., Zhang, J., et al. (2017). Genomic variation associated with local adaptation of weedy rice during de-domestication. Nat. Commun. 8:15323.
Qu, L. J., Chen, J., Liu, M., Pan, N., Okamoto, H., Lin, Z., et al. (2003). Molecular cloning and functional analysis of a novel type of Bowman-Birk inhibitor gene family in rice. Plant Physiol. 133, 560–570. doi: 10.1104/pp.103.024810
Roach, T., Colville, L., Beckett, R. P., Minibayeva, F. V., Havaux, M., and Kranner, I. (2015). A proposed interplay between peroxidase, amine oxidase and lipoxygenase in the wounding-induced oxidative burst in Pisum sativum seedlings. Phytochemistry 112, 130–138. doi: 10.1016/j.phytochem.2014.06.003
Ross, C. W. (1974). Plant Physiology Laboratory Manual. Belmont, CA: Wadsworth Publ. Comp. Inc.
Rouhier, N., Gelhaye, E., Sautiere, P. E., Brun, A., Laurent, P., Tagu, D., et al. (2001). Isolation and characterization of a new peroxiredoxin from poplar sieve tubes that uses either glutaredoxin or thioredoxin as a proton donor. Plant Physiol. 127, 1299–1309. doi: 10.1104/pp.010586
Santamaria, M. E., Cambra, I., Martinez, M., Pozancos, C., González-Melendi, P., Grbic, V., et al. (2012). Gene pyramiding of peptidase inhibitors enhances plant resistance to the spider mite Tetranychus urticae. PLoS One 7:e43011. doi: 10.1371/journal.pone.0043011
Sarla, N., and Swamy, B. P. M. (2005). Oryza glaberrima: a source for the improvement of Oryza sativa. Curr. Sci. 89, 955–963.
Schürmann, P., and Buchanan, B. B. (2008). The ferredoxin/thioredoxin system of oxygenic photosynthesis. Antioxid. Redox Signal. 10, 1235–1274. doi: 10.1089/ars.2007.1931
Shahidi-Noghabi, S., Van Damme, E. J. M., and Smagghe, G. (2008). Carbohydrate-binding activity of the type-2 ribosome-inactivating protein SNA-I from elderberry (Sambucus nigra) is a determining factor for its insecticidal activity. Phytochemistry 69, 2972–2978. doi: 10.1016/j.phytochem.2008.09.012
Shahidi-Noghabi, S., Van Damme, E. J. M., and Smagghe, G. (2009). Expression of Sambucus nigra agglutinin (SNA-I’) from elderberry bark in transgenic tobacco plants results in enhanced resistance to different insect species. Transg. Res. 18, 249–259. doi: 10.1007/s11248-008-9215-2
Shi, J., Fu, X. Z., Peng, T., Huang, X. S., Fan, Q. J., and Liu, J. H. (2010). Spermine pretreatment confers dehydration tolerance of citrus in vitro plants via modulation of antioxidative capacity and stomatal response. Tree Physiol. 30, 914–922. doi: 10.1093/treephys/tpq030
Shukla, N., Awasthi, R. P., Rawat, L., and Kumar, J. (2012). Biochemical and physiological responses of rice (Oryza sativa L.) as influenced by Trichoderma harzianum under drought stress. Plant Physiol. Biochem. 54, 78–88. doi: 10.1016/j.plaphy.2012.02.001
Silva, J. C., Gorenstein, M. V., Li, G. Z., Vissers, J. P. C., and Geromanos, S. J. (2006). Absolute quantification of proteins by LCMSE: a virtue of parallel ms acquisition. Mol. Cell. Proteomics 5, 144–156. doi: 10.1074/mcp.m500230-mcp200
Sperotto, R. A., Buffon, G., Schwambach, J., and Ricachenevsky, F. K. (2018a). Crops responses to mite infestation: it’s time to look at plant tolerance to meet the farmers’ needs. Front. Plant Sci. 9:556. doi: 10.3389/fpls.2018.00556
Sperotto, R. A., Buffon, G., Schwambach, J., and Ricachenevsky, F. K. (2018b). Checkmite!? Is the resistance to phytophagous mites on short and stocky wild Oryza species? Front. Plant Sci. 9:321. doi: 10.3389/fpls.2018.00321
Stenberg, J. A., and Muola, A. (2017). How should plant resistance to herbivores be measured? Front. Plant Sci. 8:663. doi: 10.3389/fpls.2017.00663
Stirbet, A., and Govindjee. (2011). On the relation between the Kautsky effect (chlorophyll a fluorescence induction) and Photosystem II: basics and applications of the OJIP fluorescence transient. J. Photochem. Photobiol. B 104, 236–257. doi: 10.1016/j.jphotobiol.2010.12.010
Strasser, R. J., Srivastava, A., and Tsimilli-Michael, M. (2000). “The fluorescence transient as a tool to characterize and screen photosynthetic samples,” in Probing Photosynthesis: Mechanism, Regulation and Adaptation, eds M. Yunus, U. Pathre, and P. Mohanty (London: Taylor & Francis), 443–480.
Sun, L., Ren, H., Liu, R., Li, B., Wu, T., Sun, F., et al. (2010). An h-type thioredoxin functions in tobacco defense responses to two species of viruses and an abiotic oxidative stress. Mol. Plant Microbe Interact. 23, 1470–1485. doi: 10.1094/mpmi-01-10-0029
Szareski, V. J., Carvalho, I. R., da Rosa, T. C., Dellagostin, S. M., de Pelegrin, A. J., Barbosa, M. H., et al. (2018). Oryza wild species: an alternative for rice breeding under abiotic stress conditions. Am. J. Plant Sci. 9:1093. doi: 10.4236/ajps.2018.96083
Tahir, J., and Dijkwel, P. (2016). β-Substituting alanine synthases: roles in cysteine metabolism and abiotic and biotic stress signalling in plants. Funct. Plant Biol. 43, 307–323. doi: 10.1071/fp15272
Terra, W. R., and Ferreira, C. (2005). “Biochemistry of digestion,” in Comprehensive Molecular Insect Science. Pergamon, ed. L. I. Gilbert (Oxford: Elsevier), 171–204.
Thomson, V. P., Cunningham, S. A., Ball, M. C., and Nicotra, A. B. (2003). Compensation for herbivory by Cucumis sativus through increased photosynthetic capacity and efficiency. Oecologia 134, 167–175. doi: 10.1007/s00442-002-1102-6
Vandenbussche, F., Desmyter, S., Ciani, M., Proost, P., Peumans, W. J., and van Damme, E. J. M. (2004). Analysis of the in planta antiviral activity of elderberry ribosome-inactivating proteins. Eur. J. Biochem. 271, 1508–1515. doi: 10.1111/j.1432-1033.2004.04059.x
Vaughan, D. A. (1994). The Wild Relatives of Rice: A Genetic Resources Handbook. Los Baños: International Rice Research Institute, 137.
Veasey, E. A., da Silva, E. F., Schammass, E. A., Oliveira, G. C. X., and Ando, A. (2008). Morphoagronomic genetic diversity in american wild rice species. Braz. Arch. Biol. Technol. 51, 95–104.
Wang, Y., Kwon, S. J., Wu, J., Choi, J., Lee, Y. H., Agrawal, G. K., et al. (2014). Transcriptome analysis of early responsive genes in rice during Magnaporthe oryzae infection. Plant Pathol. J. 30, 343–354. doi: 10.5423/ppj.oa.06.2014.0055
Wei, G. Q., Liu, R. S., Wang, Q., and Liu, W. Y. (2004). Toxicity of two type II ribosome-inactivating proteins (cinnamomin and ricin) to domestic silkworm larvae. Arch. Insect Biochem. Physiol. 57, 160–165. doi: 10.1002/arch.20024
Wiśniewski, J. R., Zougman, A., Nagaraj, N., and Mann, M. (2009). Universal sample preparation method for proteome analysis. Nature Met. 6, 359–362. doi: 10.1038/nmeth.1322
Xue, X., Cao, Z. X., Zhang, X. T., Wang, Y., Zhang, Y. F., Chen, Z. X., et al. (2016). Overexpression of OsOSM1 enhances resistance to rice sheath blight. Plant Dis. 100, 1634–1642. doi: 10.1094/pdis-11-15-1372-re
Yang, J., Duan, G., Li, C., Liu, L., Han, G., Zhang, Y., et al. (2019). The crosstalks between Jasmonic Acid and other plant hormone signaling highlight the involvement of Jasmonic Acid as a core component in plant response to biotic and abiotic stresses. Front. Plant Sci. 10:1349. doi: 10.3389/fpls.2019.01349
Yoo, K. S., Ok, S. H., Jeong, B. C., Jung, K. W., Cui, M. H., Hyoung, S., et al. (2011). Single cystathionine β-synthase domain-containing proteins modulate development by regulating the thioredoxin system in Arabidopsis. Plant Cell 23, 3577–3594. doi: 10.1105/tpc.111.089847
Youssefian, S., Nakamura, M., Orudgev, E., and Kondo, N. (2001). Increased cysteine biosynthesis capacity of transgenic tobacco overexpressing an O-acetylserine(thiol) lyase modifies plant responses to oxidative stress. Plant Physiol. 126, 1001–1011. doi: 10.1104/pp.126.3.1001
Zhang, C. J., Zhao, B. C., Ge, W. N., Zhang, Y. F., Song, Y., Sun, D. Y., et al. (2011). An apoplastic h-type thioredoxin is involved in the stress response through regulation of the apoplastic reactive oxygen species in rice. Plant Physiol. 157, 1884–1899. doi: 10.1104/pp.111.182808
Zhou, M., and Memelink, J. (2016). Jasmonate-responsive transcription factors regulating plant secondary metabolism. Biotechnol. Adv. 34, 441–449. doi: 10.1016/j.biotechadv.2016.02.004
Zou, L., Stout, M. J., and Dunand, R. T. (2004). The effects of feeding by the rice water weevil, Lissorhoptrus oryzophilus Kuschel, on the growth and yield components of rice, Oryza sativa. Agr. Forest Entomol. 6, 47–54. doi: 10.1111/j.1461-9555.2004.00203.x
Žurga, S., Pohleven, J., Kos, J., and Sabotič, J. (2015). β-Trefoil structure enables interactions between lectins and protease inhibitors that regulate their biological functions. J. Biochem. 158, 83–90. doi: 10.1093/jb/mvv025
Keywords: osmotin, proline, protease inhibitors, proteome, resistance, wild species
Citation: Buffon G, Blasi ÉAR, Lamb TI, Adamski JM, Schwambach J, Ricachenevsky FK, Bertolazi A, Silveira V, Lopes MCB and Sperotto RA (2021) Oryza sativa cv. Nipponbare and Oryza barthii as Unexpected Tolerance and Susceptibility Sources Against Schizotetranychus oryzae (Acari: Tetranychidae) Mite Infestation. Front. Plant Sci. 12:613568. doi: 10.3389/fpls.2021.613568
Received: 02 October 2020; Accepted: 21 January 2021;
Published: 10 February 2021.
Edited by:Jeremy Astier, INRA UMR1347 Agroécologie, France
Reviewed by:Sang-Wook Han, Chung-Ang University, South Korea
Kanako Bessho-Uehara, Tohoku University, Japan
Qian Qian, Chinese Academy of Agricultural Sciences, China
Copyright © 2021 Buffon, Blasi, Lamb, Adamski, Schwambach, Ricachenevsky, Bertolazi, Silveira, Lopes and Sperotto. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Raul Antonio Sperotto, email@example.com; firstname.lastname@example.org