# Effect of substrate voltage on EPG recordings of ingestion and probing behavior in Diaphorina citri (Hemiptera: Liviidae).

The Asian citrus psyllid, Diaphorina citri Kuyawama (Hemiptera:
Liviidae), transmits a phloem-limited alpha proteobacterium (Candidatus
Liberibacter asiaticus) that is the putative causal agent of Asiatic
huanglongbing. The psyllid was first reported in Florida in 1998, and
the disease was first detected in Florida in 2005 (Halbert &
Manjunath 2004; Grafton-Cardwell et al. 2013). The disease kills citrus
trees, resulting in extensive economic losses (Farnsworth et al. 2014).
In the USA, the Asian citrus psyllid was first reported in Texas in 2001
(French et al. 2001) and California in 2008 (Grafton-Cardwell 2010), and
the disease was subsequently detected in both Texas (Kunta et al. 2014)
and California in 2012 (Kumagai et al. 2013). Because Candidatus
Liberibacter asiaticus is phloem limited, it is important to understand
the phloem-related feeding behaviors of D. citri. However, in this
study, we examined how experimental methods might influence results and
impact the ability to compare those results among laboratories that have
used different methods.

The most rigorous experimental method for studying hemipteran probing and ingesting behavior is through electropenetrography (EPG), a technique first developed over 50 yr ago (McLean & Kinsey 1964). The original equipment design has been improved greatly in terms of both sensitivity and usability (Tjallingii 1978, 1985; Backus & Bennett 2009). The output is a temporal pattern of voltage shifts reminiscent of an electrocardiogram or electroencephalogram. Stereotypical repetitive patterns within these recordings are called waveforms. This method has been used to study a number of psyllid species: D. citri (Bonani et al. 2010; Zhu et al. 2010; Serikawa et al. 2013; Luo et al. 2015), Bactericera cockerelli (Sulc) (Hemiptera: Triozidae) (Butler et al. 2012; Pearson et al. 2014; Sandanayaka et al. 2014), and Cacopsylla pyri L. (Hemiptera: Psyllidae) (Civolani et al. 2011, 2013).

Waveforms in D. citri were originally characterized and correlated as non-probing (NP), pathway (C), xylem ingestion (G), phloem contact (D), phloem salivation (E1), and phloem ingestion (E2) (Bonani et al. 2010). Subsequently, the non-probing waveform was subdivided into walking (NP), and standing still (Z) (Youn et al. 2011). Furthermore, pathway has been subdivided into initial stylet penetration and salivation (A), epidermis mesophyll sheath salivation (B), and pathway (C) (Yang et al. 2011; Cen et al. 2012).

The objective of this study was to document behavioral changes in D. citri as the applied voltage changed from 20 mV to 150 mV to 600 mV DC. Recently, it was suggested that an upper limit of 150 mV be used for the applied voltage due to distortions in the appearance of the recorded waveforms at higher voltages (Pearson et al. 2014), but this was for another psyllid species with no examinations made of differing voltages on insect feeding behavior. Previous studies in M. E. Roger's laboratory used a substrate voltage of 150 mV DC. However, this voltage is higher than typically reported in the literature (Table 1). In the present study, we examined a range of applied voltages to determine if the applied voltage might influence D. citri's feeding behavior.

Changing the applied voltage in EPG experiments can be used to improve the quality of the recorded signal. Although we know that the applied voltage will affect the insect at some point, we do not know exactly where that critical point occurs. As apparent in Table 1, many researchers do not report the applied voltage. This may be because of a doubt expressed by McLean that the current flowing through an aphid was sufficient to be detectable by the insect (McLean 1977). In turn, this doubt was influenced by earlier work in mosquitoes suggesting that the current flowing through the insect was less than the conjectured static charge generated during flight (Kashin & Wakeley 1965). These assumptions were not tested. Recent data indicate that some insects will alter their behavior in response to the applied voltage (Cervantes & Backus 2015). Because it has been shown that some insects respond to current within the ranges achievable by modern equipment, it is important to determine the effect this might have on interpreting results or comparing results among researchers.

Materials and Methods

PLANTS AND PSYLLIDS

We used Citrus sinensis (L. Osbeck.) Valencia orange scion on Kuharske citrange (C. sinensis x Poncirus trifoliate (L.) Raf.) rootstock (Rutaceae) planted in 3.92 L black plastic pots measuring 18 cm at the rim and 18 cm deep filled with Fafard Professional Custom Mix (Agawam, Massachusetts). Plants were pruned to a height of 51 cm from soil surface. As needed, artificial light was provided by high-pressure sodium lamps (16:8 h L:D photoperiod). Additional data were used in this manuscript (described later), and some of those data come from experiments that used C. sinensis Midsweet orange scion on Kuharske rootstock grown under the aforementioned conditions.

Psyllids used for EPG recordings were obtained from a Candidatus Liberibacter asiaticus -negative colony maintained in the laboratory at the Citrus Research and Education Center, Lake Alfred, Florida. The colony was tested periodically for Candidatus Liberibacter asiaticus by using polymerase chain reaction methods (Li et al. 2006) but has never tested positive (results not shown). We did not sex or age the psyllids, although sex-based behavioral differences had been reported before (Serikawa 2011). The plants used for colony maintenance were the same cultivar and from the same source as the plants used in the experiment.

EPG METHODS

All insects were recorded feeding on the abaxial surface of young leaves. Young leaves were light green, not fully expanded, and of unknown age. The main criterion for defining a young leaf was stickiness. If you rub your finger over a mature leaf, the waxy coating makes the leaf feel slick. You can apply pressure, and your finger will still slide across the surface. Young leaves lacked this waxy coating, and the leaf would get damaged if you forced your finger across the surface. Such leaves are found only within 3 to 5 leaves of the apical meristem.

We used two 4-channel AC-DC monitors (Backus & Bennett 2009) custom-built by William H. Bennett (EPG Equipment Co., Otterville, Missouri) in DC mode with 20, 150, or 600 mV substrate voltage and 160x adjusted amplification at the control box, and 100x fixed amplification in the head-stage amplifier. Data were acquired through a DI710 AD converter (Akron, Ohio) using Windaq software at a sampling rate of 100 Hz per channel. Psyllids were tethered using a 2 cm long and 25.4 [micro]m diameter gold wire (Sigma Cohn Corp., Vernon, New York) attached to thoracic tergites by using silver glue (1:1:1 w/w/w, white glue:water:silver flake [8-10 [micro]m; Inframat Advanced Materials, Manchester, Connecticut]). The other end of the gold wire was attached to a 23 mm long (0.48 mm diameter) copper wire by using the silver glue. This wire was soldered to a 20 mm long and 1.14 mm diameter brass nail that was inserted into the unit's head amp that was set to an impedance of 109 ohms. To complete the circuit, a 10 cm long (2 mm diameter) copper wire was inserted into the water-saturated soil of the pot with the plant. There was a 30 min starvation period from the time the insects were removed from the colony until they were placed on the plant. All insects were wired during this period. Psyllids were not chilled or anesthetized with C[O.sub.2]. Recording was started before psyllids were placed on the plant to ensure that all recordings started in the NP behavior. Recordings were 23 h in duration. We chose to measure only the original 6 waveforms (NP, C, D, E1, E2, G) (Bonani et al. 2010).

To reduce electronic noise, the plants and insects were placed in a Faraday cage built as an aluminum frame covered with pure copper screen (16 x 16 mesh = 0.15 mm wire spaced 1.58 mm apart). A ring stand was placed in the Faraday cage to support the electronics and facilitate placement of the insect on a specific leaf. In all cases, the leaf was held with double-sided tape to the outside flat surface of at 35 mm diameter plastic Petri dish. This helped ensure that the insect was feeding at a specific location, and it prevented plant movement from breaking the electrical circuit. Light was provided by overhead fluorescent lights (24:0 h L:D photoperiod). Room temperature was maintained at 26.6 [degrees]C.

PROCESSING RAW DATA

Ebert 2.0 (http://www.crec.ifas.ufl.edu/extension/epg/sas.shtml) was used to analyze the data. Ebert 2.0 is a version of Ebert 1.0 (Ebert et al. 2015), modified specifically for analysis of data from psyllids. Ebert 2.0 calculated all variables in the Sarria workbook (Sarria et al. 2009) except those involving pd, F, and E1e. While some variables were removed, others were added for a more detailed analysis of the D and G waveforms. Variables measured included: number of probes to first D, number of probes to first G, number of D, total duration of D, duration of Np before first D, duration of Np before first G, mean duration of D, average number of D per probe, average number of G per probe, time from first probe to first D, time from first probe to first G, time from start of probe with first D to first D, time from start of probe with first G to first G, number of probes after first D, number of probes after first G, number of probes <3 min after first D, number of probes <3 min after first G, number of sustained G (>10 min), time to first sustained G (>10 min), duration of longest D, duration of longest G, and percentage of probing in D. We also added standard deviation of the duration of waveforms Np, C, G, D, E1, and E2. There were in total 85 variables available in the following analyses (for a list, see supplementary file online at http://purl.fcla.edu/fcla/entomologist/browse).

DATA ANALYSES

Synopsis

All variables were transformed to improve model fit based on a quantile plot of the residuals. Treatment differences were assessed using the Tukey-Kramer multiple-comparison procedure as implemented in Proc Glimmix in SAS (SAS Institute 1988). Use of mixed-model ANOVA results in improved power of statistical tests compared with non-parametric ANOVA (Gbur et al. 2012). Given the results from this analysis, additional analyses were performed using data from previous experiments as a means to provide further context (Table 2). The equipment and setup were the same for all of these past experiments as for this experiment, except that the voltage used in these other trials was always 150 mV. These data provide small differences that can be compared with the observed effects of voltage.

These data were analyzed using Proc StepDisc to reduce the number of variables. The variables identified in this step were subsequently used in the Proc Discrim procedure in SAS to get Mahalanobis distances and to get a classification error rate. The Mahalanobis distance is a unitless measure of the difference between groups or treatments. It is useful for comparing multivariate treatment effects when the units of the independent variables are different. The classification error rate is where SAS builds a statistical model but leaves one observation out of the data set. It then uses the statistical model to try and predict where that observation comes from. SAS repeats the procedure for every insect. The goal is to have more correct classifications than would be predicted by chance alone.

Details of Data Analysis

There are two options when using Proc Glimmix for dealing with heteroscedastic variables, and models with residuals that are not normally distributed. Historically, these problems were handled by transforming the data, but in Proc Glimmix one can change the underlying distribution. Both options were tried, with counts either Poisson distributed or square root transformed, and time was either Gamma distributed or log transformed. Percentages were either Beta distributed or logit transformed. If logit transformation did not help, then arcsine square root transformation was attempted. The other continuous variables were either log transformed or distributed lognormal, exponential, or tcentral (t-distribution). A quantile plot was examined to determine the best approach. In general, counts were square root transformed, percentages were logit transformed, and everything else was log(variable + 0.1) transformed (for specific cases, see supplemental material online at http://purl.fcla.edu/fcla/entomologist/browse). These transformations were done with the sole purpose of improving the statistical models. We also attempted transforming all variables to a standard mean and variance using Proc Standardize, but no qualitative difference in outcome from multivariate procedures was observed. Therefore, the results we present are from non-standardized values.

We used a stepwise discriminant analysis (Proc StepDisc) to select a limited number of variables from the list of 85. Because these 85 variables were all correlated to some degree, it was possible that StepDisc could find a solution that was dependent on a specific variable entering the model and could find a better model if that variable was made unavailable. After the first run, one or more of the variables that entered this first model were removed and Proc StepDisc was run again. Stepwise discriminant analysis is similar in concept to stepwise regression analysis. There is a dependent variable and many independent variables. The method searches through the independent variables to find the one that is most significant (P-value to enter). It adds that variable to the model, and then recalculates all the P-values for all variables in the model. It searches these values and removes variables that have become non-significant (P-value to remove). Initially we ran the analysis using a P-value to enter of 0.08 and to remove of 0.06, but values of 0.3 and 0.2, respectively, were also tried. The best model was selected based on having the fewest number of independent variables and the ability of the model to distinguish between treatments. This model is not presented, but the variables from this model were used in Proc Discrim. By using a smaller number of variables, we eliminated most of the problems with correlation amongst the variables, and we were able to use nearly all the insects (287 out of 290 in the data set). We used Mahalanobis distances to understand how far apart the treatments were from each other. We used misclassifications to understand how accurate the model was in classifying the data.

We next focus on the problem of running an experiment with a large number of variables. It is seldom clear in the EPG literature how many variables are subjected to a statistical analysis. It is possible that researchers have reported all the variables that they tried, or that they tried all the variables and only reported the ones that showed a significant difference. In the latter case, this probably represents examination of 30 or more variables. Given this number, what is the probability of finding a few significant variables by chance alone?

In our case, we had 85 variables. Given that the null hypothesis was true and assuming that the outcome for each variable was independent of the outcome for the others, the probability of finding exactly 2 significant variables is [(0.05).sup.2*] [(0.95).sup.85-2] (the general formula is [p.sup.r][q.sup.n-r] where q = 1 - p) multiplied by the number of ways that one can withdraw 2 objects from a list of 85. The general formula for doing this is where r is the number of items to be selected from a list of n items. In our specific case, this is to get a probability of about 0.13. Do this for all outcomes from 0 to 85, as in . Take the sum of the outcomes for the first three values in this series to find the probability of observing two or fewer significant variables. The probability of finding more than three significant variables is 1 minus this value or about 0.8 (Ross 1984). This does not mean that one should dismiss these significant outcomes. However, with so few significant differences, it does mean one should view these outcomes with some suspicion as the significant outcomes could have happened by chance alone. The exact number of significant variables necessary to rise above blind chance is related to the a used to declare significance, and the number of variables examined. As stated earlier, the assumption of independence needed in these calculations is not valid. It is unclear if these equations overestimate or underestimate the problem because that would be determined by whether the correlations between variables with significant differences was greater or less than the correlations between variables without significant differences.

All variables discussed in this manuscript are by insect. We used the name of this variable even after statistical analysis. However, the value displayed no longer matches the name. Thus, the total duration of E1 for each insect becomes a mean when these totals are averaged across insects. An elegant naming system was proposed to solve this issue for non-sequential variables (Backus et al. 2007), but the naming conventions become somewhat less user friendly when applied to the already long names for the sequential variables.

Results

The Tukey-Kramer multiple-comparison procedure resulted in 2 variables that showed significant treatment differences: number of probes after the first G, and number of short probes after the first G. If significance is determined by a P-value less than 0.05, then there is a 0.2 chance of finding 2 or fewer significant results and a 0.8 chance of finding 3 or more significant results. So there were fewer significant differences than we would expect by chance alone.

We could assume symmetry where if 2 variables were withdrawn from a group of 85 variables, the probability that they have significantly different means (significance is declared if a < 0.05) equals the probability that the 2 samples have significantly similar means (by symmetry, the null-hypothesis is true if a > 0.95) (Ebert et al. 1998). This should be true if the null-hypothesis is true. In the current study, there were 6 variables with P-values above 0.95: DurScndZ, NumD, NumE2, TtlDurD, PrcntPrbD, and TmFrstSusDFrstPrb (see supplemental data). Thus, there were more variables showing similarity than we should expect given that there were only 2 variables with P-values less than 0.05.

In the next step, we took the control treatments from experiments that we had done over the last 4 yr using the same methods as reported herein, but always at 150 mV applied voltage. We then compared the results from this experiment against the control treatments from those other experiments. The goal was to understand how the differences between the voltage treatments compared to the differences among other controls (Table 2). We noted that the distances between the voltage treatments (square box) were small relative to the other differences in the table, and within the voltage treatments there were no significant differences at a = 0.05 (Table 3).

The distance between the different voltage treatments was less than the distance between any of the other Valencia trials. The average Mahalanobis distance between the voltage treatments was 0.66 (Table 4). The average within Valencia trial distance excluding the voltage treatments was greater than this distance (1.97 versus 0.66) and greater than the distance between any of the voltage treatments and the other Valencia trials (1.56). The same pattern holds if one examines the number of significant differences as a fraction of the total number of differences. None of the distances within the voltage treatment were significant, but many of the other distances were significant (Table 4). To eliminate the possibility that the results are an artifact of having a treatment that had pesticide residues, we also present the results having eliminated trial A from the analysis. While this changed the numbers, it had no effect on the overall conclusion that the voltage treatments were more similar to each other than the controls from other experiments conducted over the last 4 yr.

We looked at the ability of this model to predict the placement of a new insect. If the model was good, then all new insects would be correctly placed into their respective treatment/trial. Trial L was the most successful, with half of the insects from this trial being correctly placed (Table 5). Trial B was next at 29%. In general, the classification success rate was very low (median 10.5% correct). Although there may be some significant differences in these data, the best model that we could develop was ineffective at distinguishing differences between the voltage treatments. The misclassified insects for E and G would often be placed in trial A. Misclassified insects for F were frequently misclassified as coming from trial H. So the effect of voltage on the psyllid was generally less than the difference between two controls run a few months apart.

We next tested for differences between the psyllid feeding on 2 suitable hosts: Valencia and Midsweet. There were 23 variables that showed significant differences (Table 6). There are several features that are worth noting. Firstly, there were no significant variables associated with E2. Secondly, psyllids on Midsweet spent more time ingesting xylem than did psyllids on Valencia. Furthermore, the psyllid spent less time trying to find the xylem, and they made fewer repeated probes after finding the xylem. Thirdly, psyllids on Midsweet spent less time in C, thereby indicating that Midsweet had some feature that made it easier for the psyllid to reach phloem or xylem (Table 6). The Mahalanobis distance in psyllid ingesting behavior between Valencia and Midsweet was 5.4 with a P-value of 0.0004. It is therefore likely that there was a biologically significant difference in the behavior of the psyllid on these hosts.

Discussion

The original goal of this experiment was to document the effect of voltage on D. citri's probing and ingesting behavior. We looked at 85 variables representing accepted measures for assessing significant differences in EPG data (Sarria et al. 2009, http://www.epgsystems.eu/systems.htm) or simple extensions of variables therein (e.g., Sarria has "time to first E2" so we added "time to first D"). Only 2 variables showed a significant treatment difference, and we showed that by chance alone one should expect more than 2 significant variables. We also showed that there were an unusually high number of variables with P-values of 0.95 or higher. We showed through a discriminant analysis that all the voltage treatments were clustered closer together than one should expect from running a series of controls. Furthermore, the classification model was unable to correctly identify the treatment from which the voltage data originated. For these reasons, we conclude that D. citri is unable to detect voltages of 600 mV or less under these experimental conditions. It would strengthen this conclusion if we had been able to identify a voltage that was detectable by D. citri, but that value lies outside the range of methods currently reported in the literature. That said, it is possible that waveforms may change shape somewhat at higher voltages as was observed for B. cockerelli (Pearson et al. 2014). This should be expected because the increased voltage improves the signal-to-noise ratio from the resistance component while making it worse for the emf component. It is therefore unclear if a change in waveform in B. cockerelli was associated with a concomitant change in behavior, though this may be the case.

To better understand the effect of voltage, we included some metadata. As an unexpected outcome, we were able to show that there was a detectable shift in the behavior of D. citri feeding on Midsweet versus Valencia oranges. The difference was associated with xylem ingestion being more attractive for psyllids on Midsweet and for psyllids on Midsweet to have fewer probes. There are a few possibilities for how this occurs. It is possible that host plant volatiles are different, thereby influencing a psyllid's perception of host quality. It is possible that small changes in morphology provide some mechanical protection in Valencia, or that Valencia has more chemical defenses in terms of the number of oil glands per leaf. Finally, it is possible that Midsweet has a xylem composition that is more favored by the psyllid or that Valencia has more chemical defenses in the xylem that discourage psyllid ingestion. It has already been shown that small morphological differences can influence psyllid behavior (Ammar et al. 2013, 2014). If future research can identify specific causes, it might be possible to breed citrus to be a less favorable host for the psyllid.

We end with a few words about the data analysis. We cannot say that past research using only univariate methods is inappropriate. Partly this is because it is unclear how past analyses were performed. Do all researchers using the Sarria workbook (Sarria et al. 2009) perform a statistical analysis of all the variables calculated therein? If so, do they report all the variables that were significant or only a subset of the significant variables that made the most sense given their understanding of the biology of the organisms studied? How many variables were calculated in manuscripts that did not list the software used in the analysis? Our perception of the results should be modified by the number of significant variables compared to the total number of variables examined. One approach to get around the problems in having a large number of univariate tests is to use multivariate methods. These would include approaches like discriminant analysis (as we described here), factor analysis (which would include principle component analysis [Serrano et al. 2000]), and cluster analysis. They cannot replace more simple methods; rather, they provide a more holistic perspective as part of a comprehensive and balanced analysis of the data from experiments such as the one presented herein.

Acknowledgments

This research was conducted through funding received from the Citrus Research and Development Foundation.

References Cited

Ammar E-D, Hall DG, Shatters Jr RG. 2013. Stylet morphometrics and citrus leaf vein structure in relation to feeding behavior of the Asian citrus psyllid Diaphorina citri, vector of citrus huanglongbing bacterium. PLoS One 8: e59914.

Ammar E-D, Richardson ML, Abdo Z, Hall DG, Shatters Jr RG. 2014. Differences in stylet sheath occurrence and the fibrous ring (sclerenchyma) between xCitroncirus plants relatively resistant or susceptible to adults of the Asian citrus psyllid Diaphorina citri (Hemiptera: Liviidae). PLoS One 9: e110919.

Backus EA, Bennett WH. 2009. The AC-DC correlation monitor: new EPG design with flexible input resistors to detect both R and emf components for any piercingsucking hemipteran. Journal of Insect Physiology 55: 869-884.

Backus EA, Cline AR, Ellerseick MR, Serrano MS. 2007. Lygus herperus (Hemiptera: Miridae) feeding on cotton: new methods and parameters for analysis of nonsequential electrical penetration graph data. Annals of the Entomological Society of America 100: 296-310.

Bonani JP, Fereres A, Garzo E, Miranda MP, Appezzato-Da-Gloria B, Lopes JRS. 2010. Characterization of electrical penetration graphs of the Asian citrus psyllid, Diaphorina citri, in sweet orange seedlings. Entomologia Experimentalis et Applicata 134: 35-49.

Butler CD, Walker GP, Trumble JT. 2012. Feeding disruption of potato psyllid, Bactericera cockerelli, by imidacloprid as measured by electrical penetration graphs. Entomologia Experimentalis et Applicata 142: 247-257.

Cen Y, Yang C, Holford P, Beattie GAC, Spooner-Hart RN, Liang G, Deng X. 2012. Feeding behaviour of the Asiatic citrus psyllid, Diaphorina citri, on healthy and huanglongbing-infected citrus. Entomologia Experimentalis et Applicata 143: 13-22.

Cervantes FA, Backus EA. 2015. Amount of applied signal voltage affects behaviors in EPG recordings of Lygus lineolaris adults on pinhead cotton squares. Annual Meeting of the Entomological Society of America, Nov 2015. Entomological Society of America, Minneapolis, Minnesota.

Civolani S, Leis M, Grandi G, Garzo E, Pasqualini E, Musacchi S, Chicca M, Castaldelli G, Rossi R, Freddy Tjallingii W. 2011. Stylet penetration of Cacopsylla pyri; an electrical penetration graph (EPG) study. Journal of Insect Physiology 57: 14071419.

Civolani S, Grandi G, Chicca M, Pasqualini E, Fano EA, Musacchi S. 2013. Probing behaviour of Cacopsylla pyri on a resistant pear selection. Journal of Applied Entomology 137: 365-375.

Ebert TA, Fargo WS, Cartwright B, Hall FR. 1998. Randomization tests: an example using morphological differences in Aphis gossypii (Hemiptera: Aphididae). Annals of the Entomological Society of America 91: 761-770.

Ebert TA, Backus E, Cid M, Fereres A, Rogers M. 2015. A new SAS program for behavioral analysis of electrical penetration graph data. Computers and Electronics in Agriculture 116: 80-87.

Farnsworth D, Grogan KA, van Bruggen AHC, Moss CB. 2014. The potential economic cost and response to greening in Florida. Choices 29: Available at: http://choicesmagazine.org/choices-magazine/submitted-articles/the-potential-economic-cost-and-response-to-greening-in-florida-citrus#sthash.7raVXfTy.dpuf (last accessed 12 Jun 2016)

French JV, Kahlke CJ, Da Graca JV. 2001. First record of the Asian citrus psylla, Diaphorina citri Kuwayama (Homoptera: Psyllidae), in Texas. Subtropical Plant Science 53: 14-15.

Gbur EE, Stroup WW, McCarter KS, Durham S, Young LJ, Christman M, West M, Kramer M. 2012. Analysis of Generalized Linear Mixed Models in the Agricultural and Natural Resources Sciences. American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America, Madison, Wisconsin.

Grafton-Cardwell EE. 2010. How serious is the threat of Asian citrus psyllid and huanglongbing? Citrograph 1: 8-10.

Grafton-Cardwell EE, Stelinski LL, Stansly PA. 2013. Biology and management of Asian citrus psyllid, vector of the huanglongbing pathogens. Annual Review of Entomology 58: 413-432.

Halbert SE, Manjunath KL. 2004. Asian citrus psyllids (Sternorrhyncha: Psyllidae) and greening disease of citrus: a literature review and assessment of risk in Florida. Florida Entomologist 87: 330-353.

Kashin P, Wakeley HG. 1965. An insect "bitometer". Nature (London) 208: 462-464.

Kumagai LB, LeVasque CS, Bloomquist CL, Madishetty K, Guo YY, Woods PW, Rooney-Latham S, Rascoe J, Gallindo T, Schnabel D, Polek M. 2013. First report of Candidatus Liberibacter asiaticus associated with citrus huanglongbing in California. Plant Disease 97: 283.

Kunta M, Setamou M, Louzada ES, da Graca J, Bartels DW. 2014. Citrus huanglongbing incidence, spread, and current situation in Texas. Phytopathology 104 (Suppl. 3): S3.64.

Li WB, Hartung JS, Levy L. 2006. Quantitative real-time PCR for detection and identification of Candidatus Liberibacter species associated with citrus huanglongbing. Journal of Microbiological Methods 66: 104-115.

Luo X, Yen AL, Powell KS, Wu F, Wang Y, Zeng L, Yang Y, Cen Y. 2015. Feeding behavior of Diaphorina citri (Hemiptera: Liviidae) and its acquisition of 'Candidatus Liberibacter asiaticus', on huanglongbing-infected Citrus reticulata leaves of several maturity stages. Florida Entomologist 98: 186-192.

McLean DL. 1977. An electrical measurement system for studying aphid probing behavior, pp. 277-290 In Harris KF, Maramorosch K [eds.], Aphids as Virus Vectors. Academic Press, Inc., New York, New York.

McLean DL, Kinsey MG. 1964. A technique for electronically recording aphid feeding and salivation. Nature (London) 202: 1358-1359.

Miranda MP, Yamamoto PT, Garcia RB, Lopes JPA, Lopes JRS. 2016. Thiamethoxam and imidacloprid drench applications on sweet orange nursery trees disrupt the feeding and settling behaviour of Diaphorina citri (Hemiptera: Liviidae). Pest Management Science 10.1002/ps4213.

Mustafa T, Horton DR, Cooper WR, Swisher KD, Zack RS, Pappu HR, Munyaneza JE. 2015. Use of electrical penetration graph technology to examine transmission of 'Candidatus Liberibacter solanacearum' to potato by three haplotypes of potato psyllid (Bactericera cockerelli; Hemiptera: Triozidae). PLoS One 10: e0138946.

Pearson CC, Backus EA, Shugart HJ, Munyaneza JE. 2014. Characterization and correlation of EPG waveforms of Bactericera cockerelli (Hemiptera: Triozidae): variability in waveform appearance in relation to applied signal. Annals of the Entomological Society of America 107: 650-666.

Ross S. 1984. A First Course in Probability (2nd ed.). Macmillan Publishing Company, New York, New York.

Sandanayaka WRM, Moreno A, Tooman LK, Page-Weir NEM, Fereres A. 2014. Stylet penetration activities linked to the acquisition and inoculation of Candidatus Liberibacter solanacearum by its vector tomato potato psyllid. Entomologia Experimentalis et Applicata 151: 170-181.

Sarria EV, Cid M, Garzo E, Fereres A. 2009. Excel workbook for automatic parameter calculation of EPG data. Computers and Electronics in Agriculture 67: 35-42.

SAS Institute. 1988. SAS/STAT User's Guide, release 6.03 edition. SAS Institute Inc., Cary, North Carolina.

Serikawa RH. 2011. Electrical penetration graph investigations of Asian citrus psyllid (Diaphorina citri Kuwayama) feeding behavior: effects of insecticides on the potential transmission of Candidatus Liberibacter asiaticus. Ph.D. dissertation, University of Florida, Gainesville, Florida.

Serikawa RH, Backus EA, Rogers ME. 2013. Probing behaviors of adult Asian citrus psyllid (Hemiptera: Liviidae) are not appreciably affected by soil application of field rate aldicarb to citrus. Florida Entomologist 96: 1334-1342.

Serrano MS, Backus EA, Cardona C. 2000. Comparison of AC electronic monitoring and field data for estimating tolerance to Empoasca kraemeri (Homoptera: Cicadellidae) in common bean genotypes. Journal of Economic Entomology 93: 1796-1809.

Tjallingii WF. 1978. Electronic recording of penetration behaviour by aphids. Entomologia Experimentalis et Applicata 24: 721-730.

Tjallingii WF. 1985. Electrical nature of recorded signals during stylet penetration by aphids. Entomologia Experimentalis et Applicata 38: 177-186.

Yang C-l, Cen Y-j, Liang G-w, Chen H-y. 2011. Study on the electrical penetration graph of Diaphorina citri. Journal of South China Agricultural University 32: 49-56.

Youn Y, Backus EA, Serikawa RH, Stelinski LL. 2011. Correlation of an electrical penetration graph waveform with walking by Asian citrus psyllid, Diaphorina citri (Hemiptera: Psyllidae). Florida Entomologist 94: 1084-1087.

Zhu HM, Zeng XN, Zaka SM, Cen YJ. 2010. Influence of the essential oil of guava leaves on the feeding behavior of the Asian citrus psyllid. Huan Jing Kun Chong Xue Bao 32: 483-487.

Timothy A. Ebert * and Michael E. Rogers

University of Florida, Citrus Research and Education Center, Department of Entomology and Nematology, 700 Experiment Station Rd., Lake Alfred, Florida 33850, USA; E-mail: tebert@ufl.edu (T. A. E.), mrgrs@ufl.edu (M. E. R.)

* Corresponding author; E-mail: tebert@ufl.edu (T. A. E.)

Supplementary material in Florida Entomologist 99(3) (Sep 2016) is online at http://purl.fcla.edu/fcla/entomologist/browse

The most rigorous experimental method for studying hemipteran probing and ingesting behavior is through electropenetrography (EPG), a technique first developed over 50 yr ago (McLean & Kinsey 1964). The original equipment design has been improved greatly in terms of both sensitivity and usability (Tjallingii 1978, 1985; Backus & Bennett 2009). The output is a temporal pattern of voltage shifts reminiscent of an electrocardiogram or electroencephalogram. Stereotypical repetitive patterns within these recordings are called waveforms. This method has been used to study a number of psyllid species: D. citri (Bonani et al. 2010; Zhu et al. 2010; Serikawa et al. 2013; Luo et al. 2015), Bactericera cockerelli (Sulc) (Hemiptera: Triozidae) (Butler et al. 2012; Pearson et al. 2014; Sandanayaka et al. 2014), and Cacopsylla pyri L. (Hemiptera: Psyllidae) (Civolani et al. 2011, 2013).

Waveforms in D. citri were originally characterized and correlated as non-probing (NP), pathway (C), xylem ingestion (G), phloem contact (D), phloem salivation (E1), and phloem ingestion (E2) (Bonani et al. 2010). Subsequently, the non-probing waveform was subdivided into walking (NP), and standing still (Z) (Youn et al. 2011). Furthermore, pathway has been subdivided into initial stylet penetration and salivation (A), epidermis mesophyll sheath salivation (B), and pathway (C) (Yang et al. 2011; Cen et al. 2012).

The objective of this study was to document behavioral changes in D. citri as the applied voltage changed from 20 mV to 150 mV to 600 mV DC. Recently, it was suggested that an upper limit of 150 mV be used for the applied voltage due to distortions in the appearance of the recorded waveforms at higher voltages (Pearson et al. 2014), but this was for another psyllid species with no examinations made of differing voltages on insect feeding behavior. Previous studies in M. E. Roger's laboratory used a substrate voltage of 150 mV DC. However, this voltage is higher than typically reported in the literature (Table 1). In the present study, we examined a range of applied voltages to determine if the applied voltage might influence D. citri's feeding behavior.

Changing the applied voltage in EPG experiments can be used to improve the quality of the recorded signal. Although we know that the applied voltage will affect the insect at some point, we do not know exactly where that critical point occurs. As apparent in Table 1, many researchers do not report the applied voltage. This may be because of a doubt expressed by McLean that the current flowing through an aphid was sufficient to be detectable by the insect (McLean 1977). In turn, this doubt was influenced by earlier work in mosquitoes suggesting that the current flowing through the insect was less than the conjectured static charge generated during flight (Kashin & Wakeley 1965). These assumptions were not tested. Recent data indicate that some insects will alter their behavior in response to the applied voltage (Cervantes & Backus 2015). Because it has been shown that some insects respond to current within the ranges achievable by modern equipment, it is important to determine the effect this might have on interpreting results or comparing results among researchers.

Materials and Methods

PLANTS AND PSYLLIDS

We used Citrus sinensis (L. Osbeck.) Valencia orange scion on Kuharske citrange (C. sinensis x Poncirus trifoliate (L.) Raf.) rootstock (Rutaceae) planted in 3.92 L black plastic pots measuring 18 cm at the rim and 18 cm deep filled with Fafard Professional Custom Mix (Agawam, Massachusetts). Plants were pruned to a height of 51 cm from soil surface. As needed, artificial light was provided by high-pressure sodium lamps (16:8 h L:D photoperiod). Additional data were used in this manuscript (described later), and some of those data come from experiments that used C. sinensis Midsweet orange scion on Kuharske rootstock grown under the aforementioned conditions.

Psyllids used for EPG recordings were obtained from a Candidatus Liberibacter asiaticus -negative colony maintained in the laboratory at the Citrus Research and Education Center, Lake Alfred, Florida. The colony was tested periodically for Candidatus Liberibacter asiaticus by using polymerase chain reaction methods (Li et al. 2006) but has never tested positive (results not shown). We did not sex or age the psyllids, although sex-based behavioral differences had been reported before (Serikawa 2011). The plants used for colony maintenance were the same cultivar and from the same source as the plants used in the experiment.

EPG METHODS

All insects were recorded feeding on the abaxial surface of young leaves. Young leaves were light green, not fully expanded, and of unknown age. The main criterion for defining a young leaf was stickiness. If you rub your finger over a mature leaf, the waxy coating makes the leaf feel slick. You can apply pressure, and your finger will still slide across the surface. Young leaves lacked this waxy coating, and the leaf would get damaged if you forced your finger across the surface. Such leaves are found only within 3 to 5 leaves of the apical meristem.

We used two 4-channel AC-DC monitors (Backus & Bennett 2009) custom-built by William H. Bennett (EPG Equipment Co., Otterville, Missouri) in DC mode with 20, 150, or 600 mV substrate voltage and 160x adjusted amplification at the control box, and 100x fixed amplification in the head-stage amplifier. Data were acquired through a DI710 AD converter (Akron, Ohio) using Windaq software at a sampling rate of 100 Hz per channel. Psyllids were tethered using a 2 cm long and 25.4 [micro]m diameter gold wire (Sigma Cohn Corp., Vernon, New York) attached to thoracic tergites by using silver glue (1:1:1 w/w/w, white glue:water:silver flake [8-10 [micro]m; Inframat Advanced Materials, Manchester, Connecticut]). The other end of the gold wire was attached to a 23 mm long (0.48 mm diameter) copper wire by using the silver glue. This wire was soldered to a 20 mm long and 1.14 mm diameter brass nail that was inserted into the unit's head amp that was set to an impedance of 109 ohms. To complete the circuit, a 10 cm long (2 mm diameter) copper wire was inserted into the water-saturated soil of the pot with the plant. There was a 30 min starvation period from the time the insects were removed from the colony until they were placed on the plant. All insects were wired during this period. Psyllids were not chilled or anesthetized with C[O.sub.2]. Recording was started before psyllids were placed on the plant to ensure that all recordings started in the NP behavior. Recordings were 23 h in duration. We chose to measure only the original 6 waveforms (NP, C, D, E1, E2, G) (Bonani et al. 2010).

To reduce electronic noise, the plants and insects were placed in a Faraday cage built as an aluminum frame covered with pure copper screen (16 x 16 mesh = 0.15 mm wire spaced 1.58 mm apart). A ring stand was placed in the Faraday cage to support the electronics and facilitate placement of the insect on a specific leaf. In all cases, the leaf was held with double-sided tape to the outside flat surface of at 35 mm diameter plastic Petri dish. This helped ensure that the insect was feeding at a specific location, and it prevented plant movement from breaking the electrical circuit. Light was provided by overhead fluorescent lights (24:0 h L:D photoperiod). Room temperature was maintained at 26.6 [degrees]C.

PROCESSING RAW DATA

Ebert 2.0 (http://www.crec.ifas.ufl.edu/extension/epg/sas.shtml) was used to analyze the data. Ebert 2.0 is a version of Ebert 1.0 (Ebert et al. 2015), modified specifically for analysis of data from psyllids. Ebert 2.0 calculated all variables in the Sarria workbook (Sarria et al. 2009) except those involving pd, F, and E1e. While some variables were removed, others were added for a more detailed analysis of the D and G waveforms. Variables measured included: number of probes to first D, number of probes to first G, number of D, total duration of D, duration of Np before first D, duration of Np before first G, mean duration of D, average number of D per probe, average number of G per probe, time from first probe to first D, time from first probe to first G, time from start of probe with first D to first D, time from start of probe with first G to first G, number of probes after first D, number of probes after first G, number of probes <3 min after first D, number of probes <3 min after first G, number of sustained G (>10 min), time to first sustained G (>10 min), duration of longest D, duration of longest G, and percentage of probing in D. We also added standard deviation of the duration of waveforms Np, C, G, D, E1, and E2. There were in total 85 variables available in the following analyses (for a list, see supplementary file online at http://purl.fcla.edu/fcla/entomologist/browse).

DATA ANALYSES

Synopsis

All variables were transformed to improve model fit based on a quantile plot of the residuals. Treatment differences were assessed using the Tukey-Kramer multiple-comparison procedure as implemented in Proc Glimmix in SAS (SAS Institute 1988). Use of mixed-model ANOVA results in improved power of statistical tests compared with non-parametric ANOVA (Gbur et al. 2012). Given the results from this analysis, additional analyses were performed using data from previous experiments as a means to provide further context (Table 2). The equipment and setup were the same for all of these past experiments as for this experiment, except that the voltage used in these other trials was always 150 mV. These data provide small differences that can be compared with the observed effects of voltage.

These data were analyzed using Proc StepDisc to reduce the number of variables. The variables identified in this step were subsequently used in the Proc Discrim procedure in SAS to get Mahalanobis distances and to get a classification error rate. The Mahalanobis distance is a unitless measure of the difference between groups or treatments. It is useful for comparing multivariate treatment effects when the units of the independent variables are different. The classification error rate is where SAS builds a statistical model but leaves one observation out of the data set. It then uses the statistical model to try and predict where that observation comes from. SAS repeats the procedure for every insect. The goal is to have more correct classifications than would be predicted by chance alone.

Details of Data Analysis

There are two options when using Proc Glimmix for dealing with heteroscedastic variables, and models with residuals that are not normally distributed. Historically, these problems were handled by transforming the data, but in Proc Glimmix one can change the underlying distribution. Both options were tried, with counts either Poisson distributed or square root transformed, and time was either Gamma distributed or log transformed. Percentages were either Beta distributed or logit transformed. If logit transformation did not help, then arcsine square root transformation was attempted. The other continuous variables were either log transformed or distributed lognormal, exponential, or tcentral (t-distribution). A quantile plot was examined to determine the best approach. In general, counts were square root transformed, percentages were logit transformed, and everything else was log(variable + 0.1) transformed (for specific cases, see supplemental material online at http://purl.fcla.edu/fcla/entomologist/browse). These transformations were done with the sole purpose of improving the statistical models. We also attempted transforming all variables to a standard mean and variance using Proc Standardize, but no qualitative difference in outcome from multivariate procedures was observed. Therefore, the results we present are from non-standardized values.

We used a stepwise discriminant analysis (Proc StepDisc) to select a limited number of variables from the list of 85. Because these 85 variables were all correlated to some degree, it was possible that StepDisc could find a solution that was dependent on a specific variable entering the model and could find a better model if that variable was made unavailable. After the first run, one or more of the variables that entered this first model were removed and Proc StepDisc was run again. Stepwise discriminant analysis is similar in concept to stepwise regression analysis. There is a dependent variable and many independent variables. The method searches through the independent variables to find the one that is most significant (P-value to enter). It adds that variable to the model, and then recalculates all the P-values for all variables in the model. It searches these values and removes variables that have become non-significant (P-value to remove). Initially we ran the analysis using a P-value to enter of 0.08 and to remove of 0.06, but values of 0.3 and 0.2, respectively, were also tried. The best model was selected based on having the fewest number of independent variables and the ability of the model to distinguish between treatments. This model is not presented, but the variables from this model were used in Proc Discrim. By using a smaller number of variables, we eliminated most of the problems with correlation amongst the variables, and we were able to use nearly all the insects (287 out of 290 in the data set). We used Mahalanobis distances to understand how far apart the treatments were from each other. We used misclassifications to understand how accurate the model was in classifying the data.

We next focus on the problem of running an experiment with a large number of variables. It is seldom clear in the EPG literature how many variables are subjected to a statistical analysis. It is possible that researchers have reported all the variables that they tried, or that they tried all the variables and only reported the ones that showed a significant difference. In the latter case, this probably represents examination of 30 or more variables. Given this number, what is the probability of finding a few significant variables by chance alone?

In our case, we had 85 variables. Given that the null hypothesis was true and assuming that the outcome for each variable was independent of the outcome for the others, the probability of finding exactly 2 significant variables is [(0.05).sup.2*] [(0.95).sup.85-2] (the general formula is [p.sup.r][q.sup.n-r] where q = 1 - p) multiplied by the number of ways that one can withdraw 2 objects from a list of 85. The general formula for doing this is where r is the number of items to be selected from a list of n items. In our specific case, this is to get a probability of about 0.13. Do this for all outcomes from 0 to 85, as in . Take the sum of the outcomes for the first three values in this series to find the probability of observing two or fewer significant variables. The probability of finding more than three significant variables is 1 minus this value or about 0.8 (Ross 1984). This does not mean that one should dismiss these significant outcomes. However, with so few significant differences, it does mean one should view these outcomes with some suspicion as the significant outcomes could have happened by chance alone. The exact number of significant variables necessary to rise above blind chance is related to the a used to declare significance, and the number of variables examined. As stated earlier, the assumption of independence needed in these calculations is not valid. It is unclear if these equations overestimate or underestimate the problem because that would be determined by whether the correlations between variables with significant differences was greater or less than the correlations between variables without significant differences.

All variables discussed in this manuscript are by insect. We used the name of this variable even after statistical analysis. However, the value displayed no longer matches the name. Thus, the total duration of E1 for each insect becomes a mean when these totals are averaged across insects. An elegant naming system was proposed to solve this issue for non-sequential variables (Backus et al. 2007), but the naming conventions become somewhat less user friendly when applied to the already long names for the sequential variables.

Results

The Tukey-Kramer multiple-comparison procedure resulted in 2 variables that showed significant treatment differences: number of probes after the first G, and number of short probes after the first G. If significance is determined by a P-value less than 0.05, then there is a 0.2 chance of finding 2 or fewer significant results and a 0.8 chance of finding 3 or more significant results. So there were fewer significant differences than we would expect by chance alone.

We could assume symmetry where if 2 variables were withdrawn from a group of 85 variables, the probability that they have significantly different means (significance is declared if a < 0.05) equals the probability that the 2 samples have significantly similar means (by symmetry, the null-hypothesis is true if a > 0.95) (Ebert et al. 1998). This should be true if the null-hypothesis is true. In the current study, there were 6 variables with P-values above 0.95: DurScndZ, NumD, NumE2, TtlDurD, PrcntPrbD, and TmFrstSusDFrstPrb (see supplemental data). Thus, there were more variables showing similarity than we should expect given that there were only 2 variables with P-values less than 0.05.

In the next step, we took the control treatments from experiments that we had done over the last 4 yr using the same methods as reported herein, but always at 150 mV applied voltage. We then compared the results from this experiment against the control treatments from those other experiments. The goal was to understand how the differences between the voltage treatments compared to the differences among other controls (Table 2). We noted that the distances between the voltage treatments (square box) were small relative to the other differences in the table, and within the voltage treatments there were no significant differences at a = 0.05 (Table 3).

The distance between the different voltage treatments was less than the distance between any of the other Valencia trials. The average Mahalanobis distance between the voltage treatments was 0.66 (Table 4). The average within Valencia trial distance excluding the voltage treatments was greater than this distance (1.97 versus 0.66) and greater than the distance between any of the voltage treatments and the other Valencia trials (1.56). The same pattern holds if one examines the number of significant differences as a fraction of the total number of differences. None of the distances within the voltage treatment were significant, but many of the other distances were significant (Table 4). To eliminate the possibility that the results are an artifact of having a treatment that had pesticide residues, we also present the results having eliminated trial A from the analysis. While this changed the numbers, it had no effect on the overall conclusion that the voltage treatments were more similar to each other than the controls from other experiments conducted over the last 4 yr.

We looked at the ability of this model to predict the placement of a new insect. If the model was good, then all new insects would be correctly placed into their respective treatment/trial. Trial L was the most successful, with half of the insects from this trial being correctly placed (Table 5). Trial B was next at 29%. In general, the classification success rate was very low (median 10.5% correct). Although there may be some significant differences in these data, the best model that we could develop was ineffective at distinguishing differences between the voltage treatments. The misclassified insects for E and G would often be placed in trial A. Misclassified insects for F were frequently misclassified as coming from trial H. So the effect of voltage on the psyllid was generally less than the difference between two controls run a few months apart.

We next tested for differences between the psyllid feeding on 2 suitable hosts: Valencia and Midsweet. There were 23 variables that showed significant differences (Table 6). There are several features that are worth noting. Firstly, there were no significant variables associated with E2. Secondly, psyllids on Midsweet spent more time ingesting xylem than did psyllids on Valencia. Furthermore, the psyllid spent less time trying to find the xylem, and they made fewer repeated probes after finding the xylem. Thirdly, psyllids on Midsweet spent less time in C, thereby indicating that Midsweet had some feature that made it easier for the psyllid to reach phloem or xylem (Table 6). The Mahalanobis distance in psyllid ingesting behavior between Valencia and Midsweet was 5.4 with a P-value of 0.0004. It is therefore likely that there was a biologically significant difference in the behavior of the psyllid on these hosts.

Discussion

The original goal of this experiment was to document the effect of voltage on D. citri's probing and ingesting behavior. We looked at 85 variables representing accepted measures for assessing significant differences in EPG data (Sarria et al. 2009, http://www.epgsystems.eu/systems.htm) or simple extensions of variables therein (e.g., Sarria has "time to first E2" so we added "time to first D"). Only 2 variables showed a significant treatment difference, and we showed that by chance alone one should expect more than 2 significant variables. We also showed that there were an unusually high number of variables with P-values of 0.95 or higher. We showed through a discriminant analysis that all the voltage treatments were clustered closer together than one should expect from running a series of controls. Furthermore, the classification model was unable to correctly identify the treatment from which the voltage data originated. For these reasons, we conclude that D. citri is unable to detect voltages of 600 mV or less under these experimental conditions. It would strengthen this conclusion if we had been able to identify a voltage that was detectable by D. citri, but that value lies outside the range of methods currently reported in the literature. That said, it is possible that waveforms may change shape somewhat at higher voltages as was observed for B. cockerelli (Pearson et al. 2014). This should be expected because the increased voltage improves the signal-to-noise ratio from the resistance component while making it worse for the emf component. It is therefore unclear if a change in waveform in B. cockerelli was associated with a concomitant change in behavior, though this may be the case.

To better understand the effect of voltage, we included some metadata. As an unexpected outcome, we were able to show that there was a detectable shift in the behavior of D. citri feeding on Midsweet versus Valencia oranges. The difference was associated with xylem ingestion being more attractive for psyllids on Midsweet and for psyllids on Midsweet to have fewer probes. There are a few possibilities for how this occurs. It is possible that host plant volatiles are different, thereby influencing a psyllid's perception of host quality. It is possible that small changes in morphology provide some mechanical protection in Valencia, or that Valencia has more chemical defenses in terms of the number of oil glands per leaf. Finally, it is possible that Midsweet has a xylem composition that is more favored by the psyllid or that Valencia has more chemical defenses in the xylem that discourage psyllid ingestion. It has already been shown that small morphological differences can influence psyllid behavior (Ammar et al. 2013, 2014). If future research can identify specific causes, it might be possible to breed citrus to be a less favorable host for the psyllid.

We end with a few words about the data analysis. We cannot say that past research using only univariate methods is inappropriate. Partly this is because it is unclear how past analyses were performed. Do all researchers using the Sarria workbook (Sarria et al. 2009) perform a statistical analysis of all the variables calculated therein? If so, do they report all the variables that were significant or only a subset of the significant variables that made the most sense given their understanding of the biology of the organisms studied? How many variables were calculated in manuscripts that did not list the software used in the analysis? Our perception of the results should be modified by the number of significant variables compared to the total number of variables examined. One approach to get around the problems in having a large number of univariate tests is to use multivariate methods. These would include approaches like discriminant analysis (as we described here), factor analysis (which would include principle component analysis [Serrano et al. 2000]), and cluster analysis. They cannot replace more simple methods; rather, they provide a more holistic perspective as part of a comprehensive and balanced analysis of the data from experiments such as the one presented herein.

Acknowledgments

This research was conducted through funding received from the Citrus Research and Development Foundation.

References Cited

Ammar E-D, Hall DG, Shatters Jr RG. 2013. Stylet morphometrics and citrus leaf vein structure in relation to feeding behavior of the Asian citrus psyllid Diaphorina citri, vector of citrus huanglongbing bacterium. PLoS One 8: e59914.

Ammar E-D, Richardson ML, Abdo Z, Hall DG, Shatters Jr RG. 2014. Differences in stylet sheath occurrence and the fibrous ring (sclerenchyma) between xCitroncirus plants relatively resistant or susceptible to adults of the Asian citrus psyllid Diaphorina citri (Hemiptera: Liviidae). PLoS One 9: e110919.

Backus EA, Bennett WH. 2009. The AC-DC correlation monitor: new EPG design with flexible input resistors to detect both R and emf components for any piercingsucking hemipteran. Journal of Insect Physiology 55: 869-884.

Backus EA, Cline AR, Ellerseick MR, Serrano MS. 2007. Lygus herperus (Hemiptera: Miridae) feeding on cotton: new methods and parameters for analysis of nonsequential electrical penetration graph data. Annals of the Entomological Society of America 100: 296-310.

Bonani JP, Fereres A, Garzo E, Miranda MP, Appezzato-Da-Gloria B, Lopes JRS. 2010. Characterization of electrical penetration graphs of the Asian citrus psyllid, Diaphorina citri, in sweet orange seedlings. Entomologia Experimentalis et Applicata 134: 35-49.

Butler CD, Walker GP, Trumble JT. 2012. Feeding disruption of potato psyllid, Bactericera cockerelli, by imidacloprid as measured by electrical penetration graphs. Entomologia Experimentalis et Applicata 142: 247-257.

Cen Y, Yang C, Holford P, Beattie GAC, Spooner-Hart RN, Liang G, Deng X. 2012. Feeding behaviour of the Asiatic citrus psyllid, Diaphorina citri, on healthy and huanglongbing-infected citrus. Entomologia Experimentalis et Applicata 143: 13-22.

Cervantes FA, Backus EA. 2015. Amount of applied signal voltage affects behaviors in EPG recordings of Lygus lineolaris adults on pinhead cotton squares. Annual Meeting of the Entomological Society of America, Nov 2015. Entomological Society of America, Minneapolis, Minnesota.

Civolani S, Leis M, Grandi G, Garzo E, Pasqualini E, Musacchi S, Chicca M, Castaldelli G, Rossi R, Freddy Tjallingii W. 2011. Stylet penetration of Cacopsylla pyri; an electrical penetration graph (EPG) study. Journal of Insect Physiology 57: 14071419.

Civolani S, Grandi G, Chicca M, Pasqualini E, Fano EA, Musacchi S. 2013. Probing behaviour of Cacopsylla pyri on a resistant pear selection. Journal of Applied Entomology 137: 365-375.

Ebert TA, Fargo WS, Cartwright B, Hall FR. 1998. Randomization tests: an example using morphological differences in Aphis gossypii (Hemiptera: Aphididae). Annals of the Entomological Society of America 91: 761-770.

Ebert TA, Backus E, Cid M, Fereres A, Rogers M. 2015. A new SAS program for behavioral analysis of electrical penetration graph data. Computers and Electronics in Agriculture 116: 80-87.

Farnsworth D, Grogan KA, van Bruggen AHC, Moss CB. 2014. The potential economic cost and response to greening in Florida. Choices 29: Available at: http://choicesmagazine.org/choices-magazine/submitted-articles/the-potential-economic-cost-and-response-to-greening-in-florida-citrus#sthash.7raVXfTy.dpuf (last accessed 12 Jun 2016)

French JV, Kahlke CJ, Da Graca JV. 2001. First record of the Asian citrus psylla, Diaphorina citri Kuwayama (Homoptera: Psyllidae), in Texas. Subtropical Plant Science 53: 14-15.

Gbur EE, Stroup WW, McCarter KS, Durham S, Young LJ, Christman M, West M, Kramer M. 2012. Analysis of Generalized Linear Mixed Models in the Agricultural and Natural Resources Sciences. American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America, Madison, Wisconsin.

Grafton-Cardwell EE. 2010. How serious is the threat of Asian citrus psyllid and huanglongbing? Citrograph 1: 8-10.

Grafton-Cardwell EE, Stelinski LL, Stansly PA. 2013. Biology and management of Asian citrus psyllid, vector of the huanglongbing pathogens. Annual Review of Entomology 58: 413-432.

Halbert SE, Manjunath KL. 2004. Asian citrus psyllids (Sternorrhyncha: Psyllidae) and greening disease of citrus: a literature review and assessment of risk in Florida. Florida Entomologist 87: 330-353.

Kashin P, Wakeley HG. 1965. An insect "bitometer". Nature (London) 208: 462-464.

Kumagai LB, LeVasque CS, Bloomquist CL, Madishetty K, Guo YY, Woods PW, Rooney-Latham S, Rascoe J, Gallindo T, Schnabel D, Polek M. 2013. First report of Candidatus Liberibacter asiaticus associated with citrus huanglongbing in California. Plant Disease 97: 283.

Kunta M, Setamou M, Louzada ES, da Graca J, Bartels DW. 2014. Citrus huanglongbing incidence, spread, and current situation in Texas. Phytopathology 104 (Suppl. 3): S3.64.

Li WB, Hartung JS, Levy L. 2006. Quantitative real-time PCR for detection and identification of Candidatus Liberibacter species associated with citrus huanglongbing. Journal of Microbiological Methods 66: 104-115.

Luo X, Yen AL, Powell KS, Wu F, Wang Y, Zeng L, Yang Y, Cen Y. 2015. Feeding behavior of Diaphorina citri (Hemiptera: Liviidae) and its acquisition of 'Candidatus Liberibacter asiaticus', on huanglongbing-infected Citrus reticulata leaves of several maturity stages. Florida Entomologist 98: 186-192.

McLean DL. 1977. An electrical measurement system for studying aphid probing behavior, pp. 277-290 In Harris KF, Maramorosch K [eds.], Aphids as Virus Vectors. Academic Press, Inc., New York, New York.

McLean DL, Kinsey MG. 1964. A technique for electronically recording aphid feeding and salivation. Nature (London) 202: 1358-1359.

Miranda MP, Yamamoto PT, Garcia RB, Lopes JPA, Lopes JRS. 2016. Thiamethoxam and imidacloprid drench applications on sweet orange nursery trees disrupt the feeding and settling behaviour of Diaphorina citri (Hemiptera: Liviidae). Pest Management Science 10.1002/ps4213.

Mustafa T, Horton DR, Cooper WR, Swisher KD, Zack RS, Pappu HR, Munyaneza JE. 2015. Use of electrical penetration graph technology to examine transmission of 'Candidatus Liberibacter solanacearum' to potato by three haplotypes of potato psyllid (Bactericera cockerelli; Hemiptera: Triozidae). PLoS One 10: e0138946.

Pearson CC, Backus EA, Shugart HJ, Munyaneza JE. 2014. Characterization and correlation of EPG waveforms of Bactericera cockerelli (Hemiptera: Triozidae): variability in waveform appearance in relation to applied signal. Annals of the Entomological Society of America 107: 650-666.

Ross S. 1984. A First Course in Probability (2nd ed.). Macmillan Publishing Company, New York, New York.

Sandanayaka WRM, Moreno A, Tooman LK, Page-Weir NEM, Fereres A. 2014. Stylet penetration activities linked to the acquisition and inoculation of Candidatus Liberibacter solanacearum by its vector tomato potato psyllid. Entomologia Experimentalis et Applicata 151: 170-181.

Sarria EV, Cid M, Garzo E, Fereres A. 2009. Excel workbook for automatic parameter calculation of EPG data. Computers and Electronics in Agriculture 67: 35-42.

SAS Institute. 1988. SAS/STAT User's Guide, release 6.03 edition. SAS Institute Inc., Cary, North Carolina.

Serikawa RH. 2011. Electrical penetration graph investigations of Asian citrus psyllid (Diaphorina citri Kuwayama) feeding behavior: effects of insecticides on the potential transmission of Candidatus Liberibacter asiaticus. Ph.D. dissertation, University of Florida, Gainesville, Florida.

Serikawa RH, Backus EA, Rogers ME. 2013. Probing behaviors of adult Asian citrus psyllid (Hemiptera: Liviidae) are not appreciably affected by soil application of field rate aldicarb to citrus. Florida Entomologist 96: 1334-1342.

Serrano MS, Backus EA, Cardona C. 2000. Comparison of AC electronic monitoring and field data for estimating tolerance to Empoasca kraemeri (Homoptera: Cicadellidae) in common bean genotypes. Journal of Economic Entomology 93: 1796-1809.

Tjallingii WF. 1978. Electronic recording of penetration behaviour by aphids. Entomologia Experimentalis et Applicata 24: 721-730.

Tjallingii WF. 1985. Electrical nature of recorded signals during stylet penetration by aphids. Entomologia Experimentalis et Applicata 38: 177-186.

Yang C-l, Cen Y-j, Liang G-w, Chen H-y. 2011. Study on the electrical penetration graph of Diaphorina citri. Journal of South China Agricultural University 32: 49-56.

Youn Y, Backus EA, Serikawa RH, Stelinski LL. 2011. Correlation of an electrical penetration graph waveform with walking by Asian citrus psyllid, Diaphorina citri (Hemiptera: Psyllidae). Florida Entomologist 94: 1084-1087.

Zhu HM, Zeng XN, Zaka SM, Cen YJ. 2010. Influence of the essential oil of guava leaves on the feeding behavior of the Asian citrus psyllid. Huan Jing Kun Chong Xue Bao 32: 483-487.

Timothy A. Ebert * and Michael E. Rogers

University of Florida, Citrus Research and Education Center, Department of Entomology and Nematology, 700 Experiment Station Rd., Lake Alfred, Florida 33850, USA; E-mail: tebert@ufl.edu (T. A. E.), mrgrs@ufl.edu (M. E. R.)

* Corresponding author; E-mail: tebert@ufl.edu (T. A. E.)

Supplementary material in Florida Entomologist 99(3) (Sep 2016) is online at http://purl.fcla.edu/fcla/entomologist/browse

Table 1. Voltages used in EPG studies on members of the Psylloidea. Voltage AC Species Reference in mV or DC Unknown DC Diaphorina (Bonani et al. 2010; Zhu et al. citri 2010; Yang et al. 2011; Youn et al. 2011; Cen et al. 2012; Luo et al. 2015; Miranda et al. 2016) DC Bactericera (Sandanayaka et al. 2014; cockerelli Mustafa et al. 2015) DC Cacopsylla (Civolani et al. 2013) pyri 15 DC Diaphorina (Serikawa 2011) citri 30 DC Bactericera (Butler et al. 2012) cockerelli 25 to 75 AC Bactericera (Pearson et al. 2014) cockerelli 20 to 600 DC Bactericera (Pearson et al. 2014) cockerelli Table 2. A list of supplemental sources of data used in this research. All plants used the Kuharske rootstock with either Valencia or Midsweet scion. The number of insects in each trial is n. Trial n Date Plant Voltage Insecticide Aa 8 17 Mar 2015 Valencia 150 mV Dibrom (b) B 37 10 Jul 2013 Midsweet 150 mV C 27 11 Jul 2013 Midsweet 150 mV D 22 10 Sep 2014 Midsweet 150 mV E 22 10 Jul 2014 Valencia 20 mV F 21 10 Jul 2014 Valencia 150 mV G 20 10 Jul 2014 Valencia 600 mV H 21 14 Nov 2014 Valencia 150 mV I 18 16 May 2013 Midsweet 150 mV K 21 16 Nov 2012 Midsweet 150 mV L 14 3 Jul 2013 Midsweet 150 mV M 20 4 Sep 2013 Midsweet 150 mV N 18 24 Jan 2014 Valencia 150 mV Oc 20 2 Apr 2015 Valencia 150 mV Admire (d) (a) Trial A had 8 insects feeding on plants treated with dibrom 21 d earlier. Plants were held in the greenhouse and subject to standard growing practices. (b) 1,2-dibromo-2,2-dichloroethyl dimethyl phosphate, CAS 300-76-5. (c) Trial O used plants previously treated with imidacloprid with residues at 30 ppb including analytes. This value was from a combined sample of many flushes, but lacks replication. (d) N-(1-[(6-chloro-3-pyridyl)methyl]-4,5-dihydroimidazol-2-yl) nitramide, CAS 105827-78-9. Table 3. Mahalanobis distances between trials (lower triangular matrix) and the associated P-value (upper triangular matrix). The box highlights the voltage experiment results. The letters are for the trials (Trl) as described in Table 2. Trl A B C D E F G H A 0.00 0.02 0.00 0.01 0.01 0.03 0.08 B 3.55 0.00 0.38 0.16 0.11 0.00 0.01 C 2.19 2.77 0.00 0.00 0.00 0.00 0.02 D 3.00 0.40 2.29 0.75 0.88 0.08 0.49 E 2.90 0.60 2.88 0.25 0.96 0.14 0.22 F 2.57 0.71 2.84 0.17 0.10 0.06 0.42 G 2.31 1.51 1.61 0.96 0.81 1.06 0.19 H 1.75 1.25 1.15 0.42 0.67 0.48 0.74 I 2.56 1.27 0.63 0.94 1.74 1.52 1.65 0.58 K 3.15 1.32 0.89 1.09 1.80 1.62 2.12 0.79 L 3.01 0.38 1.85 1.10 1.48 1.60 1.81 1.56 M 3.54 0.49 3.56 0.77 1.38 1.16 2.30 1.72 N 3.71 1.35 2.21 0.51 0.91 0.66 1.91 0.49 O 2.97 2.28 0.60 2.41 2.46 2.84 1.47 1.70 Trl I K L M N O A 0.02 0.00 0.01 0.00 0.00 0.01 B 0.01 0.01 0.60 0.30 0.01 0.00 C 0.23 0.08 0.01 0.00 0.00 0.21 D 0.10 0.05 0.10 0.16 0.42 0.00 E 0.00 0.00 0.03 0.02 0.12 0.00 F 0.01 0.01 0.02 0.04 0.28 0.00 G 0.01 0.00 0.01 0.00 0.00 0.01 H 0.34 0.16 0.03 0.00 0.46 0.00 I 0.93 0.20 0.01 0.16 0.02 K 0.14 0.15 0.00 0.26 0.03 L 0.92 1.01 0.20 0.01 0.04 M 1.62 2.04 0.91 0.00 0.00 N 0.89 0.71 2.04 2.06 0.00 O 1.31 1.26 1.44 3.90 2.60 Table 4. Average Mahalanobis distances and the mean number of significant distances for differences in probing and ingesting behavior of D. citri. "All" includes all the data; "Not A" includes all treatments except Trial A. EFG is restricted to values from treatments E (20 mV), F (150 mV), and G (600 mV). "Val-EFG" are all the trials on Valencia but not E, F, or G. "EFG*Val" are only values between EFG and all other Valencia. Val+EFG consists of trials A, E, F, G, H, N, and O. Mid-sweet consists of trials B, C, D, I, K, L, and M. "Val*Mid" is the distances between Valencia and Midsweet. Mahalanobis Mean number distance of significant distances All Not A All Not A EFG 0.66 0.66 0 0 Val-EFG 1.97 1.45 0.67 0.67 EFG*Val 1.56 1.33 0.58 0.44 Val+EFG 1.69 1.26 0.54 0.42 Midsweet 1.28 1.28 0.47 0.47 Val*Mid 1.68 1.45 0.74 0.69 Table 5. Classification results for a "new" insect. Trial is the letter designation for the trials (Table 1). Correct was the number of insects from each trial correctly placed into that trial. Error was the number of times a mistake was made. "Most common" was the trial where most of the individuals from each treatment were placed. "Number misclassified" was the number of insects placed into each trial that did not belong there. Misclassified percentage is the number misclassified divided by the sum of all misclassified for that treatment. Misclassified source is the treatment most likely to be misclassified into the listed trial. Trial Correct Error Most common A 1 7 E B 10 24 B C 4 23 O,C,D,M D 1 21 B,M E 2 20 B F 2 19 B,C,H,I G 0 20 I H 2 19 O,F,I I 2 17 H K 4 16 K L 7 7 L M 0 20 B,D,N N 2 16 E O 2 20 C Trial Number Misclassified Misclassified misclassified percentage source A 5 2.0 B,C,E,G,H B 30 12.0 D C 22 8.8 O D 27 10.8 B E 20 8.0 B F 11 4.4 H G 12 4.8 B,C H 17 6.8 F,I I 24 9.6 G K 13 5.2 O L 11 4.4 I M 20 8.0 C,D N 19 7.6 B,M O 18 7.2 C Table 6. ANOVA for differences between probing behavior of D. citri on Midsweet (M) versus Valencia (V). M > V are the variables where the mean from Midsweet exceeded that for Valencia. Numerator and denominator degrees of freedom (num df and den df, respectively), and probability of a greater difference using F test (P > F). Untransformed means and standard error of mean are also provided. M > V Variable Num df Den df P > F Probes MnPrbs 1 288 0.0206 sdPrbs 1 285 0.0499 Np MnDurNP 1 288 0.0004 sdNP 1 286 0.0007 G DurG 1 288 0.0411 MeanG 1 288 0.0492 meanNGPrb 1 267 0.0187 NumLngG 1 267 0.0186 PrcntPrbG 1 288 0.0086 D DurNnprbBfrFrstD 1 288 0.0006 meanNDPrb 1 288 0.0355 PrcntPrbD 1 288 0.0264 E1 MnDurE1 1 288 0.0045 PrcntPrbE1 1 288 0.0291 V > M Probes NumPrbs 1 288 0.0007 NP NumNP 1 288 0.0004 C NmbrC 1 288 0.0020 TtlDurC 1 288 0.0155 NmbrShrtC 1 288 0.0085 PrcntPrbC 1 288 0.0261 G CtoFrstG 1 267 0.0002 TmFrmFrstPrbFrstG 1 267 0.0007 NumPrbsAftrFrstG 1 288 0.0537 Midsweet Valencia M > V Variable Mean SE Mean SE Probes MnPrbs 4,153.9 282.50 3,511.1 309.07 sdPrbs 5,442.2 344.44 4,594.2 375.87 Np MnDurNP 3,977.0 450.73 2,532.1 493.12 sdNP 4,806.6 437.15 3,523.8 475.23 G DurG 5,118.8 318.55 4,216.8 348.51 MeanG 1,764.5 109.05 1,478.5 119.31 meanNGPrb 0.34 0.031 0.25 0.034 NumLngG 2.68 0.124 2.25 0.136 PrcntPrbG 14.00 0.981 10.17 1.073 D DurNnprbBfrFrstD 16,789 867.5 10,261 867.5 meanNDPrb 0.22 0.023 0.16 0.025 PrcntPrbD 0.38 0.035 0.27 0.039 E1 MnDurE1 70.77 6.333 42.89 6.929 PrcntPrbE1 0.76 0.093 0.46 0.102 V > M Probes NumPrbs 16.89 1.167 22.52 1.276 NP NumNP 17.20 1.165 23.00 1.275 C NmbrC 22.29 1.249 27.92 1.366 TtlDurC 25,948 942.3 31,092 942.3 NmbrShrtC 2.86 0.475 4.56 0.520 PrcntPrbC 57.00 1.915 62.97 2.095 G CtoFrstG 3.56 0.380 5.70 0.418 TmFrmFrstPrbFrstG 8,736.6 826.26 12,822 826.3 NumPrbsAftrFrstG 12.24 1.032 15.92 1.129

Printer friendly Cite/link Email Feedback | |

Title Annotation: | electropenetrography |
---|---|

Author: | Ebert, Timothy A.; Rogers, Michael E. |

Publication: | Florida Entomologist |

Article Type: | Report |

Geographic Code: | 1U5FL |

Date: | Sep 1, 2016 |

Words: | 7288 |

Previous Article: | Polistes major major and Polistes apachus (Hymenoptera: Vespidae) in Georgia, USA. |

Next Article: | Evidence of selective mating and triploidy among two social forms of Solenopsis invicta (Hymenoptera: Formicidae). |

Topics: |