1. Introduction
The whitefly,
Bemisia tabaci MEAM1 (Gennadius 1889) (Hemiptera: Aleyrodidae), is considered one of the most destructive invasive pests of many crops, including soybean (
Glycine max (L.) Merril) [
1]. Whiteflies were considered a secondary pest in soybean until the early 2000s, but today, it represents a species of economic importance [
2], responsible for causing direct damage through the removal of phloem sap and indirectly favoring formation and growth of the fungus (
Capnodium sp.) that causes sooty mold due to the excretion of honeydew.
One of the obstacles to the management of B. tabaci is the ability to quantify the damage caused by the pest, as these are indirectly inferred through losses in productivity and are invariably confounded by the impacts of other pests during the soybean production cycle. An alternative to directly assessing the damage caused by B. tabaci is to estimate the effects of its feeding on photosynthesis and the sugar and starch composition of the leaves.
The interaction between the chemical composition of leaves and the performance of whiteflies has been investigated in eggplant, and results show that varieties with higher glucose and sucrose concentrations provide better conditions for the development of
B. tabaci [
3]. A positive correlation between carbohydrate levels and whitefly populations also has been reported in eggplant and other solanaceous crops [
4], and higher amounts of glucose and fructose also contribute to higher densities of whiteflies in cotton [
5]. In soybean, the profile of sugars in the leaves change with the phenological stage, and the presence and intensity of whitefly feeding also may affect the plant’s physiology. Sugars, such as glucose and sucrose, are associated with photosynthesis and thus yield, and they also are important to
B. tabaci development [
6,
7]. These elements play a role in the development of
B. tabaci and might be useful for understanding the influence of
B. tabaci on soybean at different phenological stages under different pest population densities.
Another way to investigate the injury caused by
B. tabaci in soybeans is by assessing how infestations interfere with photosynthesis. Reductions in photosynthetic rate caused by
B. tabaci infestations has been reported for many crops, including tobacco, eggplant, zucchini, and cotton [
8,
9,
10,
11].
B. tabaci feeding has been associated with a reduction in the net photosynthetic rate in zucchini [
10,
11], and infestations has been associated with reductions in the photosynthetic rate related to decreases in chlorophyll content, photosynthetic capacity, and stomatal conductance [
12]. It is known that
B. tabaci causes damage to soybean [
13], resulting in yield losses and reduction in grain weight. The objective of this study was to identify the influence of
B. tabaci feeding on soybean by assessing effects on photosynthetic parameters, chlorophyll, and the sugar and starch content of soybean leaves, with a goal of establishing better parameters to quantify crop damage by this pest.
2. Materials and Methods
2.1. Bemisia tabaci Rearing
Whiteflies were reared in cages in a greenhouse using soybean (
Glycine max (L.) Merril) as the host. Plants were grown in 2-L pots in a greenhouse and replaced every 20 days. The initial whitefly population was collected from a cotton field at the United States Department of Agriculture—Arid Land Agricultural Research Center (USDA-ALARC) in Maricopa, Arizona, USA, and was previously identified as
B. tabaci biotype B (MEAM1). Molecular characterization of the insects was made periodically during the study for confirmation of the insect biotype [
14].
2.2. Establishment of Soybean Plants
Soybean (
var. 01057408) was sown in seed trays containing a general-purpose growing medium. When the seedlings reached the V1 phenological stage [
15], they were transplanted individually to a 5-L pot containing 9 parts growing medium/5 parts sand, both sterilized. This process was repeated after two weeks to obtain soybean in two different phenological stages for further infestation with
B. tabaci. Each plant was covered with an anti-aphid screen attached to the bottom of the pot that served as a cage for eventual whitefly infestation. The experiment was carried out under greenhouse conditions.
2.3. Soybean Infestation with Bemisia tabaci
For the infestation, adults of B. tabaci were collected from the rearing culture with an aspirator made with an 8-mm transparent plastic hose, “voil” fabric, and a 10-mL pipette tip and were placed in Eppendorf® tubes and released into the cages containing the soybean plants.
Soybean at V2 and R1 phenological stages were infested at three different levels on the same date. The control treatment contained no insects (
Table 1). The experimental unit was a single cage.
The experiment was carried out in a completely randomized design consisting of eight treatments with 12 replications. The experiment was repeated for a total of 24 replicates per treatment. Sampling of insects and plants ended simultaneously for both phenological stages.
2.4. Bemisia tabaci Population Density
A leaflet from the middle third of each plant was collected every two weeks for the duration of the experiment for a total of six samples and placed individually in a transparent plastic bag for transport to the laboratory. Samples were refrigerated until they were processed. Counts were made from the leaflet abaxial surface with a stereoscopic microscope (40× magnification). Then, the leaf area of the sample was measured with a LI-3100C (LI-COR, Inc., Lincoln, NE, USA).
2.5. Photosynthetic Parameters by Chlorophyll Fluorescence Induction Kinetics
For the Chlorophyll Fluoresce Induction Kinetics (OJIP) protocol, a leaf on the plant was dark-adapted for at least 20 min with detachable leaf clips. Measurements were then performed using the FluorPen FP100 (Photo System Instruments, Czech Republic).
A clip was attached to one of the leaflets in the upper middle third of the plant before the measurements. The device emits a saturation pulse through a beam of light, which when reflected, is read by the device. The OJIP protocol is a dark-adapted chlorophyll fluorescence technique used to measure plant stress. The parameters measured include initial fluorescence (F0), performance index on absorption basis (Pi_Abs), and turn-over number (N). These assays were performed in situ in the greenhouse weekly.
2.6. Sampling and Sample Preparation
The leaflets were collected every two weeks, stored individually in 50-mL centrifuge tubes, and immediately frozen in liquid nitrogen, after which they were stored in a freezer at −80 °C. One leaflet per plant was collected, totaling 12 leaflets per treatment in each evaluation. Samples were collected bi-weekly for six weeks.
The leaflets were individually homogenized in liquid nitrogen using a pestle and a mortar, weighted in two aliquots of 20 mg each sample, placed in 2-mL Eppendorf® tubes, and stored in a −80 °C freezer for later analyses of chemical composition (sugars and starch) and chlorophyll (a, b, and a + b).
2.7. Chemical Composition Assessment of Soybean Leaves
An assessment of the quantity and quality of sugars and starch present in the infested soybean leaves was carried out for all the leaf samples previously prepared. The evaluation of the sugar content was carried out by first preparing the samples for extraction of non-structural carbohydrates and then extracting the starch for subsequent chromatographic analysis.
2.7.1. Extraction of Non-Structural Carbohydrates
A total of 25 µL of maltose (20 mg/mL of 80% ethanol) and 1 mL of 80% ethanol were added to the tubes containing the samples; the tubes were placed in a water bath at 85 °C for 10 min and then centrifuged for 5 min at 10,000 rpm. The supernatant was transferred to a new 2-mL tube. Then, 0.5 mL of 80% ethanol was added to the tubes containing the precipitate, submitting them again to the water bath and centrifugation as described above two more times, adding the supernatants each time. The tubes containing the total supernatant were adjusted to 2 mL with 80% ethanol and centrifuged once more for 5 min at 10,000 rpm. The samples were then transferred to HPLC vials using a syringe and a 0.45-µm filter. The sugars assessed were fructose, glucose, and sucrose. Maltose could not be analyzed because the test protocols for other sugars required the addition of maltose, which was used as a surrogate.
2.7.2. Starch Extraction
Then, 25 µL of maltose (20 mg/mL in 80% ethanol) and 500 µL of 2M potassium hydroxide were added to the tubes containing the centrifuged precipitate resulting from the extraction of non-structural carbohydrates, placed in a water bath at 100 °C for 10 min, and then allowing to cool to room temperature. Then, 100 µL of 1M acetic acid (PH = 4.5), 50 µL of Tris 1M buffer solution (PH = 7.2), and 100 µL of alpha-amylase were added. The tubes were placed in a water bath at 85 °C for 30 min and allowed to cool to room temperature. A total of 100 µL of 1M acetic acid (PH = 4.5) and 500 µL of amyl glucosidase were added, and again, tubes were placed in a water bath at 100 °C for 4 min. Finally, 75 µL of 2M sodium hydroxide was added and the volume of the tubes adjusted to 2 mL with 100% ethanol for centrifugation at 6000 rpm for 5 min. The sample solution was then transferred to HPLC vials using a syringe and a 0.45-µm filter.
2.7.3. HPLC Analyses
The analyses of content were determined in 2 mL of each sample following methods proposed for measurements of nonstructural carbohydrates in plant tissues [
16] using a high-performance liquid chromatography system (Shimadzu) composed of a Luna 5u NH2 100Å column with dimensions of 250 × 4.6 mm 5 µm column (Phenomenex) and a refractive index detector (Shimadzu).
2.8. Evaluation of Chlorophyll Content
For the extraction of chlorophyll, 1 mL of 100% methanol at 4 °C was added to the 2-mL Eppendorf® tubes containing the samples previously weighted. The tubes were completely covered with aluminum foil to avoid exposure to light and placed on a shaking plate in a refrigerator at 4 °C for 48 h, inverting the tubes twice a day to keep the samples homogeneous. After 48 h, the samples were centrifuged at 10,000 rpm for 10 min, then 200 µL of the supernatants were transferred to propylene microplates to perform the readings.
The chlorophyll contents were obtained by spectrophotometry at absorbances 652 and 665, performed with the Synergy
® HT microplate reader (BioTek, Winooski, VT, USA). The absorbances were then used in the calculations of chlorophyll
a,
b, and
a +
b by multiplying the absorbance values by the respective correction factor (A
665.2 = 0.478) or (A
652 = 0.466) divided by the sample weight (mg) to obtain the value in µg/mg using the following equations (Porra et al., 1989):
2.9. Statistical Analysis
Statistical analyses were performed using the R Core Team [
17] program. Data were checked for normality, homogeneity, and suitability by the Shapiro–Wilk [
18] and Bartlett [
19] tests for further analysis of variance (ANOVA). Multiple comparisons of means were assessed by Tukey’s test. We used the functions shapiro.test and bartlett.test from the package stats in R Core Team [
17] (2019) program to test normality and homogeneity, respectively.
In the univariate analyses, we observed no significant interactions between infestation treatment and plant growth stage for any of the variables. We also found no differences between treatments or stages (
p > 0.05,
Table S1). However, when we considered all variables in a multivariate analysis with the Permutation test using the package vegan in R Core Team [
17], we observed that there was a significant effect of plant growth stage (
Supplementary material).
We conducted a correlation network analysis with the package corrr in R Core Team [
17] to assess correlations among variables. The graphical displays show which variables appear closer together and are joined by stronger paths. The paths were colored by their sign (red for positive and blue for negative). The proximity of the points was determined using multidimensional clustering.
Correlation networks were created from data on sugar content, starch, and photosynthetic parameters, connecting these factors to the number of nymphs at each infestation level (low, medium, and high) in two phenological stages, vegetative and reproductive. The more solid the color of the line, the stronger the correlation between the factors. Correlation coefficients and significance levels are provided in
Tables S2–S4.
4. Discussion
Despite its size,
B. tabaci may affect the plant’s physiology in a way that ultimately translates into yield loss [
20]. Unlike defoliating pests, like caterpillars, which reduce the photosynthetic area of the leaves, or stinkbugs, which cause direct damage to grains, the whitefly affects plant performance by feeding directly from its energy sources.
When soybean was infested in the vegetative stage, nymphal density increased faster than with infestation in the reproductive stage. It is possible that during the vegetative stage, more nutrients are available in the leaves, as the plant is allocating photosynthate to the promotion of plant growth, consequently providing better conditions for the development of
B. tabaci. After reaching the reproductive stage, these resources are reallocated to the production of flowers and grains [
21], with fewer nutrients available in the leaves, and thus perhaps fewer nutrients available to promote
B. tabaci growth. Although infestations occurred at the same time, the number of nymphs in the vegetative stage was more than four times higher than in the reproductive stage. This high density of nymphs (208.89 per cm²) likely affected leaf growth, as was reflected in reduced leaf area (
Table 2).
Ideally, in a healthy plant, the initial fluorescence (F
0) from a photosynthetic organ like a leaf should be lower compared with treatments that subject the plant to some type of stress [
22]. Feeding by
B. tabaci induced stress that was measurable through changes in chlorophyll fluorescence. While there was a positive correlation between nymphal density and F
0 at both low- and high-infestation rates, there was a negative correlation for medium infestations. The same pattern of correlations was observed for starch. This is consistent with a reduction in chlorophyll
b content in the medium infestation during the reproductive stage. Previous studies have shown that increases in chlorophyll
b content results in increased amount of starch content, both indicators of growth vigor in plants [
23,
24]. Like starch, a similar pattern was observed for F
0, which is one most representative chlorophyll fluorescence parameters and widely used in studies to measure stress on plants [
25].
F
0 also showed a correlation with starch. The correlation between these parameters was positive in the reproductive stage and negative in the vegetative stage. This pattern was observed for all infestation treatments but was strongest under low infestation. Starch is an energy reservoir, synthesized in photosynthetic cells and degraded to provide substrates for photosynthesis and sucrose for the plant [
26]. We also observed a positive correlation between starch and F
0 at higher nymph density during the vegetative stage. At the same time, the sucrose content was inversely proportional to starch. Higher concentration of starch in the leaves indicates a healthier plant, capable of providing higher levels of sucrose for plant growth, but also better conditions for
B. tabaci development [
3].
F0 and starch were positively correlated in all the infestation levels during the reproductive stage. Even when nymphs were negatively correlated with F0 and with starch at medium infestations, the positive correlation between F0 and starch was still observed, indicating that the initial fluorescence is influenced not only by stress but also by starch content. During the vegetative stage, a negative correlation between nymphs and F0 was observed at all infestation levels. However, during this stage, the correlation between F0 and starch was negative in all the treatments, supporting the relationship between starch and initial fluorescence.
For the absorption flow performance index (Pi_Abs), the expected response is the inverse of that expected for F
0. When the plants are healthy and free from stress, they are expected to perform better in the absorption of light by chlorophyll pigments in the leaves. Non-stressed plants should present higher indexes of absorption, while plants subject to stress through higher insect infestations should show a negative correlation, corresponding to a lower vitality of the plant [
22]. Treatments in the reproductive stage showed a negative correlation as expected. No significant correlation was observed for medium and high infestation in the vegetative stage, probably due to the high density of nymphs interfering with the absorption of light by the leaves. However, a positive correlation was observed in the low infestation at the vegetative stage. This could be explained due to the negative correlation of Pi_Abs with N. The Pi_Abs indicates light absorbance performance, while the turn-over (N) represents the conversion rate of this light energy to chemical energy. In this case, the plant would be showing a higher absorbance performance in relation to energy conversion that is occurring at a lower frequency because of nymphal feeding. Thus, the higher concentration of chlorophyll
a in the vegetative stage results in better absorption, allowing the plant to store more energy than it is consumed (i.e., a lower conversion rate (N)).
The N, or turn-over number, represents how many times the primary quinone acceptor (Q
A) was reduced to Q
A- in the interval from time zero (t
0) to reach maximum fluorescence (F
m) after receiving electrons from water using light energy. This direct electron transfer is essential in the process of converting light energy into chemical energy [
22,
27]. An increase in the absorption index would be reflected in a lower reduction rate of the quinone acceptor (Q
A) such that the plant would be absorbing more light but not converting this light energy to photoassimilates. Previous study has shown that increases in F
0, which is inverse to Pi_Abs, represents more Q
A- formation [
27], i.e., the apparently higher light absorption does not mean the plant has better photosynthetic performance. In the medium-infestation treatment, N was negatively correlated with all the parameters studied except F
0 and nymphs. Higher levels of chlorophyll (
a +
b) were observed in the vegetative treatments, and during the reproductive treatments, the medium infestation showed the lowest concentration. We also observed the inversion in the correlations between nymphs and starch and between nymphs and F
0 in the medium-infestation-level treatments. In other words, besides a higher light absorption, the plant also depends on good photosynthetic performance to provide resources such as starch for plant health, making it a more suitable host to
B. tabaci.
Whiteflies are phloem-feeders, and phloem is responsible for transporting photoassimilates, with generally high concentrations of sugars [
28]. The primary sugars associated with the high occurrences of whiteflies are sucrose and glucose, responsible for improving survival and adult longevity and fecundity, and promoting feeding, oviposition preference, and immature development [
3]. As these sugars are also related to photosynthesis, a negative relationship between
B. tabaci infestation and photosystem activity is expected. Here, that situation resulted in poorer plant development, as indicated by smaller leaves in soybean under high insect infestation.
The vegetative stage appears to provide better conditions for the insect’s development, as evidenced by the faster growth of nymphs during this stage. Differing from that observed for the other sugars, the correlations between glucose with nymphs were non-significant, indicating that during the vegetative stage other components, such as fructose and starch, may provide better conditions for
B. tabaci development. Glucose is second only to maltose as a product of starch breakdown and a source of carbon exported from chloroplasts during its degradation [
29].
Even though glucose also is a product of starch breakdown, it did not present negative correlation with starch. First, it could be that maltose provided an energetic advantage over the export of glucose, representing two-thirds of the carbon export, while glucose represents the rest [
29]. Secondly, because it represented a smaller part in the leaf composition, it could have been affected by a reduction in the leaflet size, especially in the high-infestation treatment. In contrast, in the high infestation, there was no significant correlation between glucose and nymphs, and that might be due to compensation by the plant that reduced the concentration of resources (glucose), therefore mitigating stress on the leaves by reducing leaf size [
30].
One of the main components of honeydew excreted by whiteflies is fructose [
31,
32], which has also been reported for other hemipterans [
33,
34]. A strong positive correlation between fructose and nymphs was observed for all treatments at the vegetative stage, when the greatest number of
B. tabaci nymphs occurred. A similar pattern was observed for sucrose. Varieties with higher content of fructose and sucrose in their leaves were reported to promote higher densities and better conditions for the development of
B. tabaci in eggplants and cotton, respectively [
3,
5]. On the other hand, a positive correlation was observed between nymphs and glucose instead of fructose during the reproductive stage. This could be due to a variation in sugar availability in soybean leaves from the vegetative to reproductive stage [
6]. During the vegetative stage, sucrose, which is also a precursor of fructose, is abundant in the leaves, and levels start to decrease during the reproductive stage as it is being relocated to reproductive structures [
6,
7]. This pattern was observed for all except the high-infestation level, where a positive correlation occurred between nymphs and starch. Nonetheless, glucose is still available via starch breakdown, supplying
B. tabaci needs during the reproductive stage.
Among the photosynthetic parameters, the turn-over number N was positively correlated with nymphs at lower densities of nymphs and negatively correlated with nymphs at higher densities. These correlations also were affected by the phenological stage, whereby correlations were negative during the vegetative and positive during the reproductive stage. Nevertheless, this pattern is modified somewhat by the density of infestation. During the reproductive stage, the correlation is strongest at the lowest density of infestation and declines in strength with increasing density. During the vegetative stage, the negative correlation between nymphs and N is weakest at the lowest infestation and grows in strength as density increases at the medium infestation. Although the maximum correlation was not observed at the high infestation compared with the medium infestation, a strong negative correlation was still observed. This may have been influenced by the significant drop in chlorophyll content observed at the high infestation affecting N. The assessment of the correlation between nymphs and N may allow for development of a ranking scale of N values to predict density of nymphs and might represent a useful tool for evaluating the impact of B. tabaci in the crop, especially in large areas, where nymph monitoring can be time consuming.
B. tabaci feeding affects the plant’s physiology reflected in part by the relationships among photosynthetic parameters as well as the levels of sugars and starch. From our results, the greatest impacts of insect infestation occur during the vegetative phenological stage, when plants are more suitable for B. tabaci development, thus leading to plant stress that will likely have consequences on productivity.