RESEARCH ARTICLE
Reduced Fitness of Virulent Aphis glycines (Hemiptera: Aphididae) Biotypes May Influence the Longevity of Resistance Genes in Soybean Adam J. Varenhorst, Michael T. McCarville, Matthew E. O’Neal* Department of Entomology, Iowa State University, Ames, IA 50011, United States of America * [email protected]
Abstract
OPEN ACCESS Citation: Varenhorst AJ, McCarville MT, O’Neal ME (2015) Reduced Fitness of Virulent Aphis glycines (Hemiptera: Aphididae) Biotypes May Influence the Longevity of Resistance Genes in Soybean. PLoS ONE 10(9): e0138252. doi:10.1371/journal. pone.0138252 Editor: Mark Gijzen, Agriculture and Agri-Food Canada, CANADA Received: May 19, 2015 Accepted: August 27, 2015 Published: September 15, 2015 Copyright: © 2015 Varenhorst et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Sustainable use of insect resistance in crops require insect resistance management plans that may include a refuge to limit the spread of virulence to this resistance. However, without a loss of fitness associated with virulence, a refuge may not prevent virulence from becoming fixed within a population of parthenogenetically reproducing insects like aphids. Aphid-resistance in soybeans (i.e., Rag genes) prevent outbreaks of soybean aphid (Aphis glycines), yet four biotypes defined by their capacity to survive on aphid-resistant soybeans (e.g., biotype-2 survives on Rag1 soybean) are found in North America. Although fitness costs are reported for biotype-3 on aphid susceptible and Rag1 soybean, it is not clear if virulence to aphid resistance in general is associated with a decrease in fitness on aphid susceptible soybeans. In laboratory assays, we measured fitness costs for biotype 2, 3 and 4 on an aphid-susceptible soybean cultivar. In addition, we also observed negative crossresistance for biotype-2 on Rag3, and biotype-3 on Rag1 soybean. We utilized a simple deterministic, single-locus, four compartment genetic model to account for the impact of these findings on the frequency of virulence alleles. When a refuge of aphid susceptible was included within this model, fitness costs and negative cross-resistance delayed the increase of virulence alleles when virulence was inherited recessively or additively. If virulence were inherited additively, fitness costs decreased the frequency of virulence. Combined, these results suggest that a refuge may prevent virulent A. glycines biotypes from overcoming Rag genes if this aphid-resistance were used commercially in North America.
Data Availability Statement: All relevant data are within the paper. Funding: North Central Soybean Research Program USDA National Institute of Food and Agriculture, Hatch Project #233507. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist.
Introduction In 2000, Aphis glycines Matsumura was first observed in the US. Prior to 2000, insecticide use in north central US soybean was infrequent [1]; however, after the establishment of A. glycines insecticide use on soybean in north central US dramatically increased [2]. The reason for the increase in insecticide use is due to soybean yield reductions of up to 40% caused by A. glycines
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
1 / 14
Fitness Costs Associated with A. glycines Biotypes
feeding [3]. Insecticides are effective at reducing A. glycines populations and preventing associated yield loss while also being cost effective [3, 4]. The insecticides commonly used to manage A. glycines populations are broad-spectrum and reduce populations of natural enemies present in soybean during application [5, 6, 7]. However, concerns for the future development of insecticide-resistant A. glycines populations if insecticides are consistently used suggest a need for additional management tools. An alternative strategy that is potentially more cost effective with negligible effects on natural enemies is soybean that contain one or more resistant to A. glycines genes (or Rag genes) [3, 8, 9]. Although there is evidence that Rag genes are effective, a limited number of varieties containing Rag1, Rag2, or Rag1+Rag2 are commercially available and their adoption is limited [8, 9, 10]. Initially the low adoption of Rag soybean was attributed to the concern for reduced yields associated with soybean varieties containing Rag genes. However, no yield drag is associated with the presence of Rag1, Rag2, or both genes [11, 12, 13, 14]. A second factor limiting the production and adoption of A. glycines-resistant soybean is the discovery of virulent (i.e., able to feed on soybean containing aphid-resistance genes) A. glycines biotypes [15, 16, 17]. To date four biotypes have been confirmed in North America [15, 16, 17]. Biotype-1 is avirulent to all of the Rag genes currently known. Biotype-2 is virulent towards the Rag1 gene [15]. Biotype-3 is virulent towards the Rag2 gene [16], and biotype-4 is virulent towards both the Rag1 and Rag2 genes, as well as pyramids constructed with both Rag1 and Rag2 [17]. The occurrence of virulence in North America is disconcerting, as their presence may limit the durability of aphid-resistance as a management tool if the use of Rag genes becomes more common. However, insect resistance management (IRM) programs that incorporate a refuge of susceptible plants can limit the frequency of virulent biotypes, but their utility is limited based on the biology and life-history traits of the target pest [18]. For example, a refuge has limited value in reducing the increase of a virulent insect that reproduces parthenogenetically, like A. glycines. However, when virulence is associated with a reduction in fitness on a susceptible host plant then a refuge can delay an increase in virulence. Gassmann et al. [19] define fitness costs as trade-offs in which alleles that confer higher fitness in one environment (e.g., on Rag soybean) reduce fitness in an alternative environment (e.g., on an aphid-susceptible soybean). Fitness costs in a virulent population result in a reduction of the frequency of virulent alleles when refuges of susceptible plants are present [20]. In addition to fitness costs, negative cross-resistance can also reduce the frequency of virulence in a population. Negative crossresistance occurs when the allele(s) that confer virulence to one resistance gene also confer hypersensitivity to another resistance gene. An alternative to negative cross-resistance is negatively correlated resistance, which occurs when the gene responsible for virulence to one source of resistance is not the same gene that is responsible for hypersensitivity to another source of resistance [20]. Crowder and Carrière [18] determined that a resistant crop for parthenogenic insects would only be effective for a short period of time unless fitness costs were associated with virulence. To date, only two studies have evaluated the fitness of a single virulent biotype, biotype-3. Wenger et al. [21] evaluated the fitness costs associated with virulence to Rag2 (i.e., biotype-3’s fitness on susceptible soybean), while Enders et al. [22] examined biotype-3 for negative crossresistance to Rag1. Wenger et al. [21] observed fitness costs for biotype-3 on susceptible soybean, and also concluded that virulence is not complete. Enders et al. [22] observed negative cross-resistance for biotype-3 on Rag1 soybean. Therefore, it is possible that the presence of these factors in virulent A. glycines biotypes could affect the rate at which virulence alleles increase in the environment. If either a fitness cost or negative cross-resistance occur within A. glycines, the impact could reduce the perceived importance of virulence as a hindrance to the successful and sustainable adoption of A. glycines-resistant soybean.
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
2 / 14
Fitness Costs Associated with A. glycines Biotypes
The objective of this study was to determine if fitness costs or negative cross-resistance are associated with virulence to Rag1 or Rag2 for A. glycines biotype-1, biotype-2, biotype-3, and biotype-4 populations on near-isogenic resistant and susceptible soybean cultivars. In addition to this evaluation, we sought to determine if A. glycines biotype-1 could obviate a fitness cost associated with virulence to Rag1 or Rag2 [23]. Finally we used a deterministic genetic model to predict the relative frequency of virulent A. glycines in light of fitness costs and negative cross-resistance observed herein.
Materials and Methods Aphid colonies and soybean cultivars Aphis glycines populations used for this experiment were obtained from The Ohio State University and the University of Wisconsin. Four populations that were defined by their response to Rag1 and Rag2 genes were utilized. A biotype avirulent to Rag1 and Rag2 (biotype-1; The Ohio State University), a biotype virulent to Rag1 but not Rag2 (biotype-2; The Ohio State University), a biotype virulent towards Rag2 but not Rag1 (biotype-3; The Ohio State University), and a biotype virulent towards Rag1 and Rag2 (biotype-4; University of Wisconsin) [11, 16, and 17]. These populations were initially collected and identified in Ohio and Wisconsin using detached leaf assays (described in [24]). Biotype-1 A. glycines were reared and maintained on susceptible soybean, biotype-2 A. glycines reared and maintained on Rag1 soybean, biotype-3 A. glycines reared and maintained on Rag2 soybean, and biotype-4 A. glycines reared and maintained on Rag1+Rag2 soybean. Cultivars used for rearing and maintaining the aphids are nearisogenic ( 75% of genes from the recurrent parent IA3027). Soybean plants used contained either no Rag genes (IA3027), Rag1 (IA3027RA1), or Rag1+Rag2 (IA3027RA12) and are nearisolines for the resistance genes Rag1 and Rag2 (approximately 93.75% genetically identical) [12, 13]. The near-isogenic line containing only the Rag2 gene is an experimental soybean line with 75% of its genes derived from the recurrent parent line IA3027 [25].
Fitness costs associated with virulence of biotype-2 and biotype-3 on susceptible soybean We hypothesized that the fitness of biotype-2 and biotype-3 A. glycines would be lower on susceptible soybean, when compared to biotype-1 (i.e., fitness costs would be associated with virulence to Rag1 or Rag2). We estimated the fitness of each biotype on Rag1, Rag2, and susceptible soybean by infesting plants with five mixed age apterous A. glycines of either biotype-1, biotype-2, or biotype-3, and measuring the population density on the plant 11 d after infestation. Fitness costs were identified if the population densities of the virulent biotype-2 and biotype-3 A. glycines were significantly lower than that of biotype-1 on the susceptible soybean. In addition, negative cross resistance was identified if the population density of a virulent biotype was significantly lower than that of avirulent biotype on soybean containing a resistance gene that it is not virulent towards. For example, if biotype-3 (virulent to Rag2, but avirulent to Rag1) has fewer aphids after 11 days on a Rag1 variety than biotype-1 (avirulent) then we would conclude that biotype-3 has negative cross-resistance. In total, we used nine treatments; each treatment was a combination of two factors, soybean cultivar (3 levels) and A. glycines biotype (3 levels). Aphis glycines individuals were transferred from colonies maintained at Iowa State University to the first full trifoliate of individually potted plants at the second trifoliate growth stage (V2 according to [26]). Plants were enclosed within mesh nets to prevent plant-to-plant movement of aphids. After 24 h we examined A. glycines populations to confirm their successful
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
3 / 14
Fitness Costs Associated with A. glycines Biotypes
establishment. The total number of A. glycines (both nymphs and adults) present on each plant was counted 11 d after initial infestation. This experiment was repeated twice using a randomized complete block design (RCBD) with five blocks per repetition (10 experimental units per treatment). Individually potted soybean plants were grown in 16-cm diameter pots in a growth chamber (E41L2C9, Percival Scientific, Incorporated, Perry, IA) using a 14:10 light:dark cycle and a constant temperature of 27°C with a relative humidity of 60%.
Fitness cost associated with virulence of biotype-4 on susceptible soybean Biotype-4 was not included in the previous experiment because a colony had not been established at the onset of that experiment. We next hypothesized that A. glycines biotype-4 fitness would be lower on susceptible soybean when compared to biotype-1 (i.e., fitness costs would be present). We tested for fitness costs using the same experimental procedures as described previously. Plants were infested with biotype-1 or biotype-4, and the population density was measured 11 d after initial infestation. We utilized a two factor experimental design, with eight total treatments. Each treatment was a combination of four soybean cultivars and two A. glycines biotypes. Three of the soybean cultivars used were the same as described in the previous experiment, with the addition of a resistant cultivar containing Rag1+Rag2 (IA3027RA12). Each cultivar was infested with either five biotype-1 or five biotype-4. The fitness of these two biotypes on each soybean cultivar was measured in population density and compared to determine the presence of fitness costs. The method for infesting A. glycines populations from the previous experiment was used. The same planting procedure and growth chamber specifications as the previous experiment were used. This experiment was repeated twice using a RCBD with five blocks within each repetition within each repetition (10 experimental units per treatment).
Impact of induced susceptibility on fitness costs of biotype-2 and biotype-3 on susceptible soybean Varenhorst et al. [23] observed obviation of Rag resistance by an improvement in fitness of avirulent A. glycines when the plant is co-infested by virulent biotypes. Our third hypothesis was that the herbivory by biotype-1 would improve the quality of susceptible soybean for biotype-2 and biotype-3. We tested this hypothesis using the same experimental design as outlined by Varenhorst et al. [23]. This design involved infesting soybean plants with an initial population (referred to as the inducer population) for 24 h prior to infesting plants with a second population (referred to as the response population). The population density of the response population 11 d after infestation is a measurement of the effect of the inducer population. The effect of the inducer population on the performance of the response population can then be assessed in comparison to plants receiving a response population but not an initial inducer population. If response populations were greater in the presence of a biotype-1 inducer population than in their absence, then we infer that induced susceptibility alleviates fitness costs. For this experiment, we used a three-factor design with ten total treatments. The three factors included plant cultivar, inducer population biotype, and response population biotype. We utilized three soybean cultivars: susceptible, Rag1, and Rag2. Four inducer populations were used: no inducer (none), 50 biotype-1 (B1), 50 biotype-2 (B2), or 50 biotype-3 (B3). Three response populations were used: five biotype-1, five biotype-2, or five biotype-3. Inducer populations were applied to the first full trifoliate when the plants reached the second trifoliate growth stage. The entire first trifoliate was then enclosed within a mesh net. The inducer
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
4 / 14
Fitness Costs Associated with A. glycines Biotypes
populations remained on the plans for the duration of the experiments. For plants that did not receive an inducer population, the first trifoliate was caged but no aphids were added. After 24 h, the response populations were added to the second full trifoliate and allowed to move freely about the plant, with the exception of the first trifoliate, which was enclosed in a net that contained the inducer population. The mesh nets are effective at separating the inducer and response populations [23]. Each potted plants was enclosed within mesh nets to prevent plantto-plant movement. Response populations were examined after 24 h to confirm successful establishment. The total number of A. glycines present in the response population on each plant was counted 11 d after the initial infestation. The experiment was repeated twice using a RCBD with three blocks per repetition (six total experimental units per treatment). The same planting procedure and growth chamber specifications as the previous experiments were used.
Statistical analysis To address each hypothesis, we analyzed the number of A. glycines per plant after 11 d. The A. glycines per plant data were log transformed to reduce heteroscadacity. All data for the first two experiments were analyzed using the PROC MIXED procedure with SAS statistical software version 9.3 (SAS Institute, Cary, NC). For both experiments, data were analyzed using an analysis of variance (ANOVA). Significant treatment effects were then separated using F-protected least-squares means with a significance level of P < 0.05. The statistical model used to analyze data from the first two experiments (i.e., fitness costs of biotype-2 and biotype-3, and fitness costs of biotype-4) included the main effects of repetition, block, soybean cultivar, and A. glycines biotype. All two and three-way interactions of the main effects were included in the model. All data from the third experiment (i.e., biotype-1 obviation of fitness costs on susceptible soybean) were analyzed using the PROC GLM procedure. Data were analyzed using an analysis of variance (ANOVA) with significant treatment effects separated using Student-NewmanKeuls (SNK) grouping with a significance level of P < 0.05. The statistical model included the main effects of repetition, block, inducer population biotype, and response population biotype. All two- and three-way interaction terms among the main effects were included in the model.
Modeling the consequences of fitness costs and negative crossresistance We hypothesized that the fitness cost and negative cross-resistance we observed would affect the rate at which virulence alleles increased within a population of A. glycines. We utilized a simple deterministic, single-locus, four compartment genetic model developed for A. glycines to track changes in the frequency of virulence alleles [23]. The model was adapted from one created for parthenogenic reproducing insects similar to A. glycines [18]. We tracked the change in the frequency of virulence alleles across 25 years with 14 generations of asexual reproduction and one generation of sexual reproduction occurring within each year. We assumed virulence to the Rag1 and Rag2 genes to be conferred by two independently segregating genes. Each virulence gene was assumed to have two alleles, with one allele conferring virulence and another conferring avirulence. Mating was assumed to be completely random with alleles returning to Hardy-Weinberg equilibrium after each year’s generation of sexual reproduction. We tracked the frequency of virulence alleles to Rag1 for 25 years beginning at the initial deployment of resistant cultivars. We report the number of years for the frequency of the Rag1 virulence allele to surpass 50% in the population. We report the allele frequency after 25 years for the Rag1 virulence allele in cases where the frequency fails to surpass 50% in 25 years.
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
5 / 14
Fitness Costs Associated with A. glycines Biotypes
Our goal was to evaluate the relative potential importance of fitness costs and negative cross-resistance for the development of virulence, not to evaluate all possible scenarios for the development of virulence to Rag genes. Therefore, we assessed a small proportion of possible scenarios for the development of virulence in A. glycines. We used values from empirical data for specific parameters in the A. glycines-soybean system whenever possible, including the field-to-field movement rate of A. glycines (5% static rate, [27]), the efficacy of the Rag1 and Rag2 genes (41% reduction in A. glycines fitness, [28]), the efficacy of a Rag1+Rag2 pyramid (59% reduction in A. glycines fitness, [28]). Fitness costs and negative cross-resistance were included as reductions in the fitness of virulent individuals on susceptible and resistant plants, respectively. The exact values used for the reduction in fitness due to fitness costs and negative cross-resistance were based on the results of experiments one and two. We included induced susceptibility effects in all of our models [23] as these can affect the rate at which virulence alleles increase in the population. The inclusion of a density-dependent increase in fitness was used to model induced susceptibility. A starting population density of 1.0×102 individuals per 1% of the landscape was used for each year of the model. Fitness of individuals in a compartment increased to 1.65 if the population density surpassed 1.0×1013 individuals per 1% of landscape. Obviation of resistance was modeled by setting the fitness of all individuals in a resistant compartment to 1.0 when the population density of homozygous virulent individuals surpassed 1.0×1013 per 1% of landscape. Obviation of fitness costs was modeled using the same method as obviation of resistance, except the obviation of fitness costs was based on the population density of homozygous avirulent individuals in the susceptible compartment of the landscape. We ran the model for a range of values for unknown parameters, including the dominance of virulence (recessive, additive, and dominant), and initial virulence allele frequency (0.02 and 0.2, [29]). We assumed a single static 25% refuge size of susceptible plants in the landscape for each run of the model. We investigated two resistance gene deployment strategies by partitioning the remaining 75% of the landscape to either plants containing only the Rag1 gene or plants containing the Rag1+Rag2 genes. In total, we ran the model 24 times once with fitness costs and negative cross-resistance included and once without for each combination of dominance of virulence (three levels), initial allele frequency (two levels), and gene deployment strategy (two levels).
Results Fitness costs associated with virulence of biotype-2 and biotype-3 on susceptible soybean We confirmed our hypothesis that the population densities of A. glycines biotype-2 and biotype-3 would be lower on susceptible soybean when compared to biotype-1 (i.e., a fitness costs associated with virulence). This was observed by analyzing data for the significance of the main effects of soybean cultivar, A. glycines biotype, and the interaction of soybean cultivar by A. glycines biotype. The population densities of the biotypes of A. glycines varied significantly by the main effects of soybean cultivar (F = 14.41; df = 2, 76; P < 0.0001), A. glycines biotype (F = 5.83; df = 2, 76; P < 0.0044), and their interaction (F = 58.00; df = 4, 76; P < 0.0001). We observed unequal population densities of the biotypes among the soybean cultivars, therefore we analyzed these data by soybean cultivar. On susceptible soybean, the population density of A. glycines biotype-1 was significantly greater than those of biotype-2 (t = 4.13; df = 2, 22; P < 0.0004) and biotype-3 (t = 6.80; df = 2, 22; P < 0.001) (Fig 1). The population density of biotype-2 was also significantly greater than that of biotype-3 (t = 2.68; df = 2, 22; P < 0.0138) on susceptible soybean (Fig 1). For the Rag1
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
6 / 14
Fitness Costs Associated with A. glycines Biotypes
Fig 1. Biotype-2 and biotype-3 A. glycines reproduction reduced on a susceptible soybean compared to aphid resistant cultivars. Note that biotype 1 is avirulent to both aphid-resistant cultivars, but virulent to the susceptible cultivar. The susceptible soybean cultivar was IA3027, the Rag1 cultivar was IA3027RA1, and the Rag2 cultivar was IA3027RA2. Data were analyzed by soybean cultivar, and capital letters indicate significance differences across biotypes (P < 0.05). doi:10.1371/journal.pone.0138252.g001
soybean cultivar, the population density of biotype-2 was significantly greater than those for biotype-1 (t = 6.72; df = 2, 22; P < 0.0001) or biotype-3 (t = 11.76; df = 2, 22; P < 0.0001) (Fig 1). The population density of biotype-1 was also significantly greater than that of biotype-3 (t = 5.04; df = 2, 22; P < 0.0001) on Rag1 soybean. On Rag2 soybean, the population density of biotype-3 was significantly greater than those of biotype-1 (t = 5.15; df = 2, 22; P < 0.0001) or biotype-2 (t = 8.68; df = 2, 22; P < 0.0001). The population density of biotype-1 was also significantly greater than that of biotype-2 (t = 3.53; df = 2, 22; P < 0.0019) on Rag2 soybean (Fig 1). In summary, these results indicate that for the susceptible soybean cultivar tested, the population densities of both biotype-2 and biotype-3 were lower than those of biotype-1. In addition, these results indicate that population densities of biotype-2 and biotype-3 were lower on Rag2 and Rag1 respectively when compared to biotype-1.
Fitness cost associated with virulence of biotype-4 on susceptible soybean In our second experiment, we confirmed our hypothesis that the population density of A. glycines biotype-4 would be lower on susceptible soybean when compared to biotype-1 (i.e., fitness costs for this biotype exist on susceptible soybean). As in the first experiment, data were analyzed for the significance of the main effects of soybean cultivar, A. glycines biotype, and the interaction of soybean cultivar by A. glycines biotype. The population densities of the biotypes of A. glycines varied significantly by the main effects of soybean cultivar (F = 50.31; df = 3, 67; P < 0.0001) and A. glycines biotype (F = 138.00; df = 1, 67; P < 0.0001). The interaction of soybean cultivar by A. glycines biotype was significant (F = 98.36; df = 3, 67; P < 0.0001), and indicated that the two A. glycines biotypes did not respond similarly to the soybean cultivars. Similar to the previous experiment we observed unequal reduction in the populations of the A. glycines biotypes among the soybean cultivars, therefore data were analyzed by soybean cultivar. On susceptible soybean, the population density of A. glycines biotype-1 was significantly greater than that of biotype-4 (t = 6.91; df = 1, 13; P < 0.0001) (Fig 2). For the Rag1
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
7 / 14
Fitness Costs Associated with A. glycines Biotypes
Fig 2. Biotype-4 A. glycines reproduction reduced on a susceptible soybean compared to an aphid resistant cultivar. Note that biotype 1 is avirulent to both aphid-resistant cultivars, but virulent to the susceptible cultivar. For this experiment the susceptible soybean cultivar was IA3027, the Rag1 cultivar was IA3027RA1, the Rag2 cultivar was IA3027RA2, and the Rag1+Rag2 cultivar was IA3027RA12. Data were analyzed by soybean cultivar, and capital letters indicate significance differences between biotype populations (P < 0.05). doi:10.1371/journal.pone.0138252.g002
(t = 4.19; df = 1, 13; P < 0.0011), Rag2 (t = 9.29; df = 1, 13; P < 0.0001), and the Rag1+Rag2 (t = 18.74; df = 1, 13; P < 0.0001) soybean cultivars, the population density of biotype-4 was significantly greater than that of biotype-1 (Fig 2). These results indicate that for the susceptible soybean cultivar that was tested, the population density of biotype-4 was lower on susceptible soybean when compared to biotype-1.
Impact of induced susceptibility on fitness costs of biotype-2 and biotype-3 on susceptible soybean For our third experiment, we confirmed that A. glycines biotype-1 is capable of alleviating the fitness costs observed for biotype-2 and biotype-3 on the susceptible soybean tested. To test this hypothesis, we analyzed these data for the significance of the main effects of inducer population biotype, response population biotype, and the interaction of inducer population biotype by response population biotype. The main effects of inducer population biotype (F = 31.05; df = 1, 41; P < 0.0001) and response population biotype (F = 28.83; df = 2, 41; P < 0.0001) both significantly affected the density of the response population. The interaction of inducer population biotype by response population biotype was also significant (F = 5.06; df = 2, 41; P < 0.0209). In agreement with results from Varenhorst et al. [23], both the presence and herbivory of biotype-1 positively affected the density of the biotype-2 and biotype-3 response populations on susceptible soybean (Fig 3). When inducer populations were absent the biotype-1 response population was 184% greater than biotype-2 and 196% greater than biotype-3. The biotype-1 response population without an inducer population was 21% greater than biotype-2 and 19% greater than biotype-3 when both had their respective biotype inducer populations present. In the presence of the biotype-1 inducer population, there were no significant differences among the population densities of biotype-1, biotype-2, or biotype-3 response populations (Fig 3). In summary, the differences between the population densities of avirulent and virulent biotypes on susceptible soybean were diminished with the addition of an inducer population with the same biotype.
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
8 / 14
Fitness Costs Associated with A. glycines Biotypes
Fig 3. Addition of an inducer population affects subsequent fitness of virulent biotypes on a susceptible cultivar. This effect was observed by varying inducer population biotypes on varying response population biotypes on susceptible soybean. For this experiment the susceptible cultivar IA3027, Rag1 cultivar IA3027RA1, and Rag2 cultivar IA3027RA2 were used. Capital letters indicate significance among treatments (P < 0.05). doi:10.1371/journal.pone.0138252.g003
Modeling the consequences of virulent biotype fitness costs and induced susceptibility In general, fitness costs and negative cross-resistance delayed the development of virulence. This delay was observed when virulence was rare (i.e., 0.02), regardless of the mode of inheritance. When virulence was common (i.e., 0.2), the delay was only observed with recessive or additive modes of inheritance (Table 1). Fitness costs resulted in the frequency of additively
Table 1. Simulated effect of induced susceptibility and fitness costs on virulence development. Defaulta
Model Factors Included Initial Vir. Allele Frequency
0.02
Default Including Fitness Costs 0.2
0.02
0.2
>25c (0.02)d
4
>25 (0.02)
18
>25 (0.02)
25
>25 (0.02)
>25 (0.03)
Recessiveb Rag1 Alone Rag1+Rag2 Pyramid Additive Rag1 Alone
>25 (0.03)
4
>25 (25 (25 (0.02)
14
>25 (25 (25 (0.06)
3
>25 (0.03)
3
Rag1+Rag2 Pyramid
>25 (0.02)
4
>25 (
Reduced Fitness of Virulent Aphis glycines (Hemiptera: Aphididae) Biotypes May Influence the Longevity of Resistance Genes in Soybean Adam J. Varenhorst, Michael T. McCarville, Matthew E. O’Neal* Department of Entomology, Iowa State University, Ames, IA 50011, United States of America * [email protected]
Abstract
OPEN ACCESS Citation: Varenhorst AJ, McCarville MT, O’Neal ME (2015) Reduced Fitness of Virulent Aphis glycines (Hemiptera: Aphididae) Biotypes May Influence the Longevity of Resistance Genes in Soybean. PLoS ONE 10(9): e0138252. doi:10.1371/journal. pone.0138252 Editor: Mark Gijzen, Agriculture and Agri-Food Canada, CANADA Received: May 19, 2015 Accepted: August 27, 2015 Published: September 15, 2015 Copyright: © 2015 Varenhorst et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Sustainable use of insect resistance in crops require insect resistance management plans that may include a refuge to limit the spread of virulence to this resistance. However, without a loss of fitness associated with virulence, a refuge may not prevent virulence from becoming fixed within a population of parthenogenetically reproducing insects like aphids. Aphid-resistance in soybeans (i.e., Rag genes) prevent outbreaks of soybean aphid (Aphis glycines), yet four biotypes defined by their capacity to survive on aphid-resistant soybeans (e.g., biotype-2 survives on Rag1 soybean) are found in North America. Although fitness costs are reported for biotype-3 on aphid susceptible and Rag1 soybean, it is not clear if virulence to aphid resistance in general is associated with a decrease in fitness on aphid susceptible soybeans. In laboratory assays, we measured fitness costs for biotype 2, 3 and 4 on an aphid-susceptible soybean cultivar. In addition, we also observed negative crossresistance for biotype-2 on Rag3, and biotype-3 on Rag1 soybean. We utilized a simple deterministic, single-locus, four compartment genetic model to account for the impact of these findings on the frequency of virulence alleles. When a refuge of aphid susceptible was included within this model, fitness costs and negative cross-resistance delayed the increase of virulence alleles when virulence was inherited recessively or additively. If virulence were inherited additively, fitness costs decreased the frequency of virulence. Combined, these results suggest that a refuge may prevent virulent A. glycines biotypes from overcoming Rag genes if this aphid-resistance were used commercially in North America.
Data Availability Statement: All relevant data are within the paper. Funding: North Central Soybean Research Program USDA National Institute of Food and Agriculture, Hatch Project #233507. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist.
Introduction In 2000, Aphis glycines Matsumura was first observed in the US. Prior to 2000, insecticide use in north central US soybean was infrequent [1]; however, after the establishment of A. glycines insecticide use on soybean in north central US dramatically increased [2]. The reason for the increase in insecticide use is due to soybean yield reductions of up to 40% caused by A. glycines
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
1 / 14
Fitness Costs Associated with A. glycines Biotypes
feeding [3]. Insecticides are effective at reducing A. glycines populations and preventing associated yield loss while also being cost effective [3, 4]. The insecticides commonly used to manage A. glycines populations are broad-spectrum and reduce populations of natural enemies present in soybean during application [5, 6, 7]. However, concerns for the future development of insecticide-resistant A. glycines populations if insecticides are consistently used suggest a need for additional management tools. An alternative strategy that is potentially more cost effective with negligible effects on natural enemies is soybean that contain one or more resistant to A. glycines genes (or Rag genes) [3, 8, 9]. Although there is evidence that Rag genes are effective, a limited number of varieties containing Rag1, Rag2, or Rag1+Rag2 are commercially available and their adoption is limited [8, 9, 10]. Initially the low adoption of Rag soybean was attributed to the concern for reduced yields associated with soybean varieties containing Rag genes. However, no yield drag is associated with the presence of Rag1, Rag2, or both genes [11, 12, 13, 14]. A second factor limiting the production and adoption of A. glycines-resistant soybean is the discovery of virulent (i.e., able to feed on soybean containing aphid-resistance genes) A. glycines biotypes [15, 16, 17]. To date four biotypes have been confirmed in North America [15, 16, 17]. Biotype-1 is avirulent to all of the Rag genes currently known. Biotype-2 is virulent towards the Rag1 gene [15]. Biotype-3 is virulent towards the Rag2 gene [16], and biotype-4 is virulent towards both the Rag1 and Rag2 genes, as well as pyramids constructed with both Rag1 and Rag2 [17]. The occurrence of virulence in North America is disconcerting, as their presence may limit the durability of aphid-resistance as a management tool if the use of Rag genes becomes more common. However, insect resistance management (IRM) programs that incorporate a refuge of susceptible plants can limit the frequency of virulent biotypes, but their utility is limited based on the biology and life-history traits of the target pest [18]. For example, a refuge has limited value in reducing the increase of a virulent insect that reproduces parthenogenetically, like A. glycines. However, when virulence is associated with a reduction in fitness on a susceptible host plant then a refuge can delay an increase in virulence. Gassmann et al. [19] define fitness costs as trade-offs in which alleles that confer higher fitness in one environment (e.g., on Rag soybean) reduce fitness in an alternative environment (e.g., on an aphid-susceptible soybean). Fitness costs in a virulent population result in a reduction of the frequency of virulent alleles when refuges of susceptible plants are present [20]. In addition to fitness costs, negative cross-resistance can also reduce the frequency of virulence in a population. Negative crossresistance occurs when the allele(s) that confer virulence to one resistance gene also confer hypersensitivity to another resistance gene. An alternative to negative cross-resistance is negatively correlated resistance, which occurs when the gene responsible for virulence to one source of resistance is not the same gene that is responsible for hypersensitivity to another source of resistance [20]. Crowder and Carrière [18] determined that a resistant crop for parthenogenic insects would only be effective for a short period of time unless fitness costs were associated with virulence. To date, only two studies have evaluated the fitness of a single virulent biotype, biotype-3. Wenger et al. [21] evaluated the fitness costs associated with virulence to Rag2 (i.e., biotype-3’s fitness on susceptible soybean), while Enders et al. [22] examined biotype-3 for negative crossresistance to Rag1. Wenger et al. [21] observed fitness costs for biotype-3 on susceptible soybean, and also concluded that virulence is not complete. Enders et al. [22] observed negative cross-resistance for biotype-3 on Rag1 soybean. Therefore, it is possible that the presence of these factors in virulent A. glycines biotypes could affect the rate at which virulence alleles increase in the environment. If either a fitness cost or negative cross-resistance occur within A. glycines, the impact could reduce the perceived importance of virulence as a hindrance to the successful and sustainable adoption of A. glycines-resistant soybean.
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
2 / 14
Fitness Costs Associated with A. glycines Biotypes
The objective of this study was to determine if fitness costs or negative cross-resistance are associated with virulence to Rag1 or Rag2 for A. glycines biotype-1, biotype-2, biotype-3, and biotype-4 populations on near-isogenic resistant and susceptible soybean cultivars. In addition to this evaluation, we sought to determine if A. glycines biotype-1 could obviate a fitness cost associated with virulence to Rag1 or Rag2 [23]. Finally we used a deterministic genetic model to predict the relative frequency of virulent A. glycines in light of fitness costs and negative cross-resistance observed herein.
Materials and Methods Aphid colonies and soybean cultivars Aphis glycines populations used for this experiment were obtained from The Ohio State University and the University of Wisconsin. Four populations that were defined by their response to Rag1 and Rag2 genes were utilized. A biotype avirulent to Rag1 and Rag2 (biotype-1; The Ohio State University), a biotype virulent to Rag1 but not Rag2 (biotype-2; The Ohio State University), a biotype virulent towards Rag2 but not Rag1 (biotype-3; The Ohio State University), and a biotype virulent towards Rag1 and Rag2 (biotype-4; University of Wisconsin) [11, 16, and 17]. These populations were initially collected and identified in Ohio and Wisconsin using detached leaf assays (described in [24]). Biotype-1 A. glycines were reared and maintained on susceptible soybean, biotype-2 A. glycines reared and maintained on Rag1 soybean, biotype-3 A. glycines reared and maintained on Rag2 soybean, and biotype-4 A. glycines reared and maintained on Rag1+Rag2 soybean. Cultivars used for rearing and maintaining the aphids are nearisogenic ( 75% of genes from the recurrent parent IA3027). Soybean plants used contained either no Rag genes (IA3027), Rag1 (IA3027RA1), or Rag1+Rag2 (IA3027RA12) and are nearisolines for the resistance genes Rag1 and Rag2 (approximately 93.75% genetically identical) [12, 13]. The near-isogenic line containing only the Rag2 gene is an experimental soybean line with 75% of its genes derived from the recurrent parent line IA3027 [25].
Fitness costs associated with virulence of biotype-2 and biotype-3 on susceptible soybean We hypothesized that the fitness of biotype-2 and biotype-3 A. glycines would be lower on susceptible soybean, when compared to biotype-1 (i.e., fitness costs would be associated with virulence to Rag1 or Rag2). We estimated the fitness of each biotype on Rag1, Rag2, and susceptible soybean by infesting plants with five mixed age apterous A. glycines of either biotype-1, biotype-2, or biotype-3, and measuring the population density on the plant 11 d after infestation. Fitness costs were identified if the population densities of the virulent biotype-2 and biotype-3 A. glycines were significantly lower than that of biotype-1 on the susceptible soybean. In addition, negative cross resistance was identified if the population density of a virulent biotype was significantly lower than that of avirulent biotype on soybean containing a resistance gene that it is not virulent towards. For example, if biotype-3 (virulent to Rag2, but avirulent to Rag1) has fewer aphids after 11 days on a Rag1 variety than biotype-1 (avirulent) then we would conclude that biotype-3 has negative cross-resistance. In total, we used nine treatments; each treatment was a combination of two factors, soybean cultivar (3 levels) and A. glycines biotype (3 levels). Aphis glycines individuals were transferred from colonies maintained at Iowa State University to the first full trifoliate of individually potted plants at the second trifoliate growth stage (V2 according to [26]). Plants were enclosed within mesh nets to prevent plant-to-plant movement of aphids. After 24 h we examined A. glycines populations to confirm their successful
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
3 / 14
Fitness Costs Associated with A. glycines Biotypes
establishment. The total number of A. glycines (both nymphs and adults) present on each plant was counted 11 d after initial infestation. This experiment was repeated twice using a randomized complete block design (RCBD) with five blocks per repetition (10 experimental units per treatment). Individually potted soybean plants were grown in 16-cm diameter pots in a growth chamber (E41L2C9, Percival Scientific, Incorporated, Perry, IA) using a 14:10 light:dark cycle and a constant temperature of 27°C with a relative humidity of 60%.
Fitness cost associated with virulence of biotype-4 on susceptible soybean Biotype-4 was not included in the previous experiment because a colony had not been established at the onset of that experiment. We next hypothesized that A. glycines biotype-4 fitness would be lower on susceptible soybean when compared to biotype-1 (i.e., fitness costs would be present). We tested for fitness costs using the same experimental procedures as described previously. Plants were infested with biotype-1 or biotype-4, and the population density was measured 11 d after initial infestation. We utilized a two factor experimental design, with eight total treatments. Each treatment was a combination of four soybean cultivars and two A. glycines biotypes. Three of the soybean cultivars used were the same as described in the previous experiment, with the addition of a resistant cultivar containing Rag1+Rag2 (IA3027RA12). Each cultivar was infested with either five biotype-1 or five biotype-4. The fitness of these two biotypes on each soybean cultivar was measured in population density and compared to determine the presence of fitness costs. The method for infesting A. glycines populations from the previous experiment was used. The same planting procedure and growth chamber specifications as the previous experiment were used. This experiment was repeated twice using a RCBD with five blocks within each repetition within each repetition (10 experimental units per treatment).
Impact of induced susceptibility on fitness costs of biotype-2 and biotype-3 on susceptible soybean Varenhorst et al. [23] observed obviation of Rag resistance by an improvement in fitness of avirulent A. glycines when the plant is co-infested by virulent biotypes. Our third hypothesis was that the herbivory by biotype-1 would improve the quality of susceptible soybean for biotype-2 and biotype-3. We tested this hypothesis using the same experimental design as outlined by Varenhorst et al. [23]. This design involved infesting soybean plants with an initial population (referred to as the inducer population) for 24 h prior to infesting plants with a second population (referred to as the response population). The population density of the response population 11 d after infestation is a measurement of the effect of the inducer population. The effect of the inducer population on the performance of the response population can then be assessed in comparison to plants receiving a response population but not an initial inducer population. If response populations were greater in the presence of a biotype-1 inducer population than in their absence, then we infer that induced susceptibility alleviates fitness costs. For this experiment, we used a three-factor design with ten total treatments. The three factors included plant cultivar, inducer population biotype, and response population biotype. We utilized three soybean cultivars: susceptible, Rag1, and Rag2. Four inducer populations were used: no inducer (none), 50 biotype-1 (B1), 50 biotype-2 (B2), or 50 biotype-3 (B3). Three response populations were used: five biotype-1, five biotype-2, or five biotype-3. Inducer populations were applied to the first full trifoliate when the plants reached the second trifoliate growth stage. The entire first trifoliate was then enclosed within a mesh net. The inducer
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
4 / 14
Fitness Costs Associated with A. glycines Biotypes
populations remained on the plans for the duration of the experiments. For plants that did not receive an inducer population, the first trifoliate was caged but no aphids were added. After 24 h, the response populations were added to the second full trifoliate and allowed to move freely about the plant, with the exception of the first trifoliate, which was enclosed in a net that contained the inducer population. The mesh nets are effective at separating the inducer and response populations [23]. Each potted plants was enclosed within mesh nets to prevent plantto-plant movement. Response populations were examined after 24 h to confirm successful establishment. The total number of A. glycines present in the response population on each plant was counted 11 d after the initial infestation. The experiment was repeated twice using a RCBD with three blocks per repetition (six total experimental units per treatment). The same planting procedure and growth chamber specifications as the previous experiments were used.
Statistical analysis To address each hypothesis, we analyzed the number of A. glycines per plant after 11 d. The A. glycines per plant data were log transformed to reduce heteroscadacity. All data for the first two experiments were analyzed using the PROC MIXED procedure with SAS statistical software version 9.3 (SAS Institute, Cary, NC). For both experiments, data were analyzed using an analysis of variance (ANOVA). Significant treatment effects were then separated using F-protected least-squares means with a significance level of P < 0.05. The statistical model used to analyze data from the first two experiments (i.e., fitness costs of biotype-2 and biotype-3, and fitness costs of biotype-4) included the main effects of repetition, block, soybean cultivar, and A. glycines biotype. All two and three-way interactions of the main effects were included in the model. All data from the third experiment (i.e., biotype-1 obviation of fitness costs on susceptible soybean) were analyzed using the PROC GLM procedure. Data were analyzed using an analysis of variance (ANOVA) with significant treatment effects separated using Student-NewmanKeuls (SNK) grouping with a significance level of P < 0.05. The statistical model included the main effects of repetition, block, inducer population biotype, and response population biotype. All two- and three-way interaction terms among the main effects were included in the model.
Modeling the consequences of fitness costs and negative crossresistance We hypothesized that the fitness cost and negative cross-resistance we observed would affect the rate at which virulence alleles increased within a population of A. glycines. We utilized a simple deterministic, single-locus, four compartment genetic model developed for A. glycines to track changes in the frequency of virulence alleles [23]. The model was adapted from one created for parthenogenic reproducing insects similar to A. glycines [18]. We tracked the change in the frequency of virulence alleles across 25 years with 14 generations of asexual reproduction and one generation of sexual reproduction occurring within each year. We assumed virulence to the Rag1 and Rag2 genes to be conferred by two independently segregating genes. Each virulence gene was assumed to have two alleles, with one allele conferring virulence and another conferring avirulence. Mating was assumed to be completely random with alleles returning to Hardy-Weinberg equilibrium after each year’s generation of sexual reproduction. We tracked the frequency of virulence alleles to Rag1 for 25 years beginning at the initial deployment of resistant cultivars. We report the number of years for the frequency of the Rag1 virulence allele to surpass 50% in the population. We report the allele frequency after 25 years for the Rag1 virulence allele in cases where the frequency fails to surpass 50% in 25 years.
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
5 / 14
Fitness Costs Associated with A. glycines Biotypes
Our goal was to evaluate the relative potential importance of fitness costs and negative cross-resistance for the development of virulence, not to evaluate all possible scenarios for the development of virulence to Rag genes. Therefore, we assessed a small proportion of possible scenarios for the development of virulence in A. glycines. We used values from empirical data for specific parameters in the A. glycines-soybean system whenever possible, including the field-to-field movement rate of A. glycines (5% static rate, [27]), the efficacy of the Rag1 and Rag2 genes (41% reduction in A. glycines fitness, [28]), the efficacy of a Rag1+Rag2 pyramid (59% reduction in A. glycines fitness, [28]). Fitness costs and negative cross-resistance were included as reductions in the fitness of virulent individuals on susceptible and resistant plants, respectively. The exact values used for the reduction in fitness due to fitness costs and negative cross-resistance were based on the results of experiments one and two. We included induced susceptibility effects in all of our models [23] as these can affect the rate at which virulence alleles increase in the population. The inclusion of a density-dependent increase in fitness was used to model induced susceptibility. A starting population density of 1.0×102 individuals per 1% of the landscape was used for each year of the model. Fitness of individuals in a compartment increased to 1.65 if the population density surpassed 1.0×1013 individuals per 1% of landscape. Obviation of resistance was modeled by setting the fitness of all individuals in a resistant compartment to 1.0 when the population density of homozygous virulent individuals surpassed 1.0×1013 per 1% of landscape. Obviation of fitness costs was modeled using the same method as obviation of resistance, except the obviation of fitness costs was based on the population density of homozygous avirulent individuals in the susceptible compartment of the landscape. We ran the model for a range of values for unknown parameters, including the dominance of virulence (recessive, additive, and dominant), and initial virulence allele frequency (0.02 and 0.2, [29]). We assumed a single static 25% refuge size of susceptible plants in the landscape for each run of the model. We investigated two resistance gene deployment strategies by partitioning the remaining 75% of the landscape to either plants containing only the Rag1 gene or plants containing the Rag1+Rag2 genes. In total, we ran the model 24 times once with fitness costs and negative cross-resistance included and once without for each combination of dominance of virulence (three levels), initial allele frequency (two levels), and gene deployment strategy (two levels).
Results Fitness costs associated with virulence of biotype-2 and biotype-3 on susceptible soybean We confirmed our hypothesis that the population densities of A. glycines biotype-2 and biotype-3 would be lower on susceptible soybean when compared to biotype-1 (i.e., a fitness costs associated with virulence). This was observed by analyzing data for the significance of the main effects of soybean cultivar, A. glycines biotype, and the interaction of soybean cultivar by A. glycines biotype. The population densities of the biotypes of A. glycines varied significantly by the main effects of soybean cultivar (F = 14.41; df = 2, 76; P < 0.0001), A. glycines biotype (F = 5.83; df = 2, 76; P < 0.0044), and their interaction (F = 58.00; df = 4, 76; P < 0.0001). We observed unequal population densities of the biotypes among the soybean cultivars, therefore we analyzed these data by soybean cultivar. On susceptible soybean, the population density of A. glycines biotype-1 was significantly greater than those of biotype-2 (t = 4.13; df = 2, 22; P < 0.0004) and biotype-3 (t = 6.80; df = 2, 22; P < 0.001) (Fig 1). The population density of biotype-2 was also significantly greater than that of biotype-3 (t = 2.68; df = 2, 22; P < 0.0138) on susceptible soybean (Fig 1). For the Rag1
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
6 / 14
Fitness Costs Associated with A. glycines Biotypes
Fig 1. Biotype-2 and biotype-3 A. glycines reproduction reduced on a susceptible soybean compared to aphid resistant cultivars. Note that biotype 1 is avirulent to both aphid-resistant cultivars, but virulent to the susceptible cultivar. The susceptible soybean cultivar was IA3027, the Rag1 cultivar was IA3027RA1, and the Rag2 cultivar was IA3027RA2. Data were analyzed by soybean cultivar, and capital letters indicate significance differences across biotypes (P < 0.05). doi:10.1371/journal.pone.0138252.g001
soybean cultivar, the population density of biotype-2 was significantly greater than those for biotype-1 (t = 6.72; df = 2, 22; P < 0.0001) or biotype-3 (t = 11.76; df = 2, 22; P < 0.0001) (Fig 1). The population density of biotype-1 was also significantly greater than that of biotype-3 (t = 5.04; df = 2, 22; P < 0.0001) on Rag1 soybean. On Rag2 soybean, the population density of biotype-3 was significantly greater than those of biotype-1 (t = 5.15; df = 2, 22; P < 0.0001) or biotype-2 (t = 8.68; df = 2, 22; P < 0.0001). The population density of biotype-1 was also significantly greater than that of biotype-2 (t = 3.53; df = 2, 22; P < 0.0019) on Rag2 soybean (Fig 1). In summary, these results indicate that for the susceptible soybean cultivar tested, the population densities of both biotype-2 and biotype-3 were lower than those of biotype-1. In addition, these results indicate that population densities of biotype-2 and biotype-3 were lower on Rag2 and Rag1 respectively when compared to biotype-1.
Fitness cost associated with virulence of biotype-4 on susceptible soybean In our second experiment, we confirmed our hypothesis that the population density of A. glycines biotype-4 would be lower on susceptible soybean when compared to biotype-1 (i.e., fitness costs for this biotype exist on susceptible soybean). As in the first experiment, data were analyzed for the significance of the main effects of soybean cultivar, A. glycines biotype, and the interaction of soybean cultivar by A. glycines biotype. The population densities of the biotypes of A. glycines varied significantly by the main effects of soybean cultivar (F = 50.31; df = 3, 67; P < 0.0001) and A. glycines biotype (F = 138.00; df = 1, 67; P < 0.0001). The interaction of soybean cultivar by A. glycines biotype was significant (F = 98.36; df = 3, 67; P < 0.0001), and indicated that the two A. glycines biotypes did not respond similarly to the soybean cultivars. Similar to the previous experiment we observed unequal reduction in the populations of the A. glycines biotypes among the soybean cultivars, therefore data were analyzed by soybean cultivar. On susceptible soybean, the population density of A. glycines biotype-1 was significantly greater than that of biotype-4 (t = 6.91; df = 1, 13; P < 0.0001) (Fig 2). For the Rag1
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
7 / 14
Fitness Costs Associated with A. glycines Biotypes
Fig 2. Biotype-4 A. glycines reproduction reduced on a susceptible soybean compared to an aphid resistant cultivar. Note that biotype 1 is avirulent to both aphid-resistant cultivars, but virulent to the susceptible cultivar. For this experiment the susceptible soybean cultivar was IA3027, the Rag1 cultivar was IA3027RA1, the Rag2 cultivar was IA3027RA2, and the Rag1+Rag2 cultivar was IA3027RA12. Data were analyzed by soybean cultivar, and capital letters indicate significance differences between biotype populations (P < 0.05). doi:10.1371/journal.pone.0138252.g002
(t = 4.19; df = 1, 13; P < 0.0011), Rag2 (t = 9.29; df = 1, 13; P < 0.0001), and the Rag1+Rag2 (t = 18.74; df = 1, 13; P < 0.0001) soybean cultivars, the population density of biotype-4 was significantly greater than that of biotype-1 (Fig 2). These results indicate that for the susceptible soybean cultivar that was tested, the population density of biotype-4 was lower on susceptible soybean when compared to biotype-1.
Impact of induced susceptibility on fitness costs of biotype-2 and biotype-3 on susceptible soybean For our third experiment, we confirmed that A. glycines biotype-1 is capable of alleviating the fitness costs observed for biotype-2 and biotype-3 on the susceptible soybean tested. To test this hypothesis, we analyzed these data for the significance of the main effects of inducer population biotype, response population biotype, and the interaction of inducer population biotype by response population biotype. The main effects of inducer population biotype (F = 31.05; df = 1, 41; P < 0.0001) and response population biotype (F = 28.83; df = 2, 41; P < 0.0001) both significantly affected the density of the response population. The interaction of inducer population biotype by response population biotype was also significant (F = 5.06; df = 2, 41; P < 0.0209). In agreement with results from Varenhorst et al. [23], both the presence and herbivory of biotype-1 positively affected the density of the biotype-2 and biotype-3 response populations on susceptible soybean (Fig 3). When inducer populations were absent the biotype-1 response population was 184% greater than biotype-2 and 196% greater than biotype-3. The biotype-1 response population without an inducer population was 21% greater than biotype-2 and 19% greater than biotype-3 when both had their respective biotype inducer populations present. In the presence of the biotype-1 inducer population, there were no significant differences among the population densities of biotype-1, biotype-2, or biotype-3 response populations (Fig 3). In summary, the differences between the population densities of avirulent and virulent biotypes on susceptible soybean were diminished with the addition of an inducer population with the same biotype.
PLOS ONE | DOI:10.1371/journal.pone.0138252 September 15, 2015
8 / 14
Fitness Costs Associated with A. glycines Biotypes
Fig 3. Addition of an inducer population affects subsequent fitness of virulent biotypes on a susceptible cultivar. This effect was observed by varying inducer population biotypes on varying response population biotypes on susceptible soybean. For this experiment the susceptible cultivar IA3027, Rag1 cultivar IA3027RA1, and Rag2 cultivar IA3027RA2 were used. Capital letters indicate significance among treatments (P < 0.05). doi:10.1371/journal.pone.0138252.g003
Modeling the consequences of virulent biotype fitness costs and induced susceptibility In general, fitness costs and negative cross-resistance delayed the development of virulence. This delay was observed when virulence was rare (i.e., 0.02), regardless of the mode of inheritance. When virulence was common (i.e., 0.2), the delay was only observed with recessive or additive modes of inheritance (Table 1). Fitness costs resulted in the frequency of additively
Table 1. Simulated effect of induced susceptibility and fitness costs on virulence development. Defaulta
Model Factors Included Initial Vir. Allele Frequency
0.02
Default Including Fitness Costs 0.2
0.02
0.2
>25c (0.02)d
4
>25 (0.02)
18
>25 (0.02)
25
>25 (0.02)
>25 (0.03)
Recessiveb Rag1 Alone Rag1+Rag2 Pyramid Additive Rag1 Alone
>25 (0.03)
4
>25 (25 (25 (0.02)
14
>25 (25 (25 (0.06)
3
>25 (0.03)
3
Rag1+Rag2 Pyramid
>25 (0.02)
4
>25 (
