<<

Biological Control 59 (2011) 354–360

Contents lists available at SciVerse ScienceDirect

Biological Control

journal homepage: www.elsevier.com/locate/ybcon

Evaluation of mortality factors using life table analysis of boliviana, a biological control agent of tropical soda apple in ⇑ Veronica Manrique a, , Rodrigo Diaz a, Stephen D. Hight b, William A. Overholt a

a Indian River Research & Education Center, University of Florida, 2199 South Rock Road, Fort Pierce, FL 34945, United States b USDA-ARS-CMAVE, Behavior and Biocontrol Research Unit, 6383 Mahan Drive, Tallahassee, FL 32308, United States

highlights graphical abstract

" is a biological control agent of in Florida. " Survival to adulthood in open cages (15%) was similar in North and Central Florida. " Survival was higher in closed cages (51%) than open cages (15%). " survival in a natural population determined from vertical life tables was 5%.

" Positive growth rate (rm = 0.3) results in reduction of S. viarum in Central Florida.

article info abstract

Article history: Tropical soda apple (TSA), Solanum viarum Dunal (), has invaded many pastures and natural Received 1 April 2011 areas in Florida. The biological control agent Gratiana boliviana Spaeth (Coleoptera:Chrysomelidae) is pro- Accepted 21 August 2011 viding adequate control of TSA stands in South and Central Florida. However, poor or no establishment of Available online 27 August 2011 this agent has occurred in northern Florida. The goal of this study was to examine the mortality factors that influence the population dynamics of G. boliviana in Florida. Horizontal life tables were constructed Keywords: by following cohorts of individuals in the laboratory and inside closed and open cages at field sites in Cen- Solanum viarum tral and North Florida. Fertility life table parameters were estimated using laboratory and field data. In addition, as part of a vertical life table analysis, TSA plants were sampled every two weeks in pastures Biological control Horizontal life-table in Central Florida, and counts of all G. boliviana and other herbivores and predators were recorded. Sur- Vertical life-table vival to adulthood was similar between Central and North Florida (open cages: 12–19%). Intrinsic mor- Fertility life-table tality (laboratory data) and biotic factors (predation) together accounted for 75% of the mortality of immature stages. Survival of in a natural population determined from vertical life tables was 5%. A complex of three mirid species (Engytatus modesta Distant, notatus Distant, and Macrol- ophus sp.) were the most abundant predators found in the field, and are known to feed on G. boliviana

and larvae. Positive growth rates (rm = 0.3) during the summer and early fall allow the beetle population to increase and provide suppression of TSA in Central Florida. Ó 2011 Elsevier Inc. All rights reserved.

1. Introduction in the introduced areas (e.g., Sutherst and Maywald, 1985; Byrne et al., 2003; Ghosheh, 2005). Climate matching is commonly used Several studies suggest that the establishment of biological to select biocontrol agents that are well adapted to the climate in control agents is greatly affected by biotic and abiotic factors found the introduced range (Sutherst and Maywald, 1985; Byrne et al., 2003; Senaratne et al., 2006; Mukherjee et al., 2011). In addition, predation and have been shown to limit the success of ⇑ Corresponding author. Fax: +1 772 460 3673. many biological control agents (e.g., Goeden and Louda, 1976; E-mail address: vero72@ufl.edu (V. Manrique).

1049-9644/$ - see front matter Ó 2011 Elsevier Inc. All rights reserved. doi:10.1016/j.biocontrol.2011.08.007 V. Manrique et al. / Biological Control 59 (2011) 354–360 355

Briese, 1986; Semple and Forno, 1987; Ghosheh, 2005). Life table pare fertility life tables constructed from laboratory and field data, studies are often employed to examine the mortality factors affect- (3) compare survivorship derived from a vertical life table with ing insect populations in the field (e.g., Royama, 1981; Carey, 1993, that calculated from horizontal life tables, and (4) quantify preda- 2001; Kuhar et al., 2002; Zanuncio et al., 2006; Yang et al., 2010). tor densities found in the field. Moreover, life tables provide estimations of demographic parame- ters including age structure, reproductive and intrinsic rates of 2. Materials and methods increase, and generation times. These parameters are useful when comparing potential population growth of different agents and for 2.1. and plants the development of mass rearing methods. Horizontal life tables are constructed by following the fate of a cohort of individuals G. boliviana used for the experiments were obtained from a over the course of their lives, while vertical life tables examine the mass rearing facility operated by the Florida Department of Agri- age structure of a population at a single moment in time culture & Consumer Services, Division of Plant Industry in Fort (Bellows et al., 1992). The latter approach is commonly used for Pierce, FL. Insects were originally collected in Southeast Brazil in insects with overlapping generations. the early 2000s. Tropical soda apple (hereafter TSA), Solanum viarum Dunal TSA plants were grown from seeds, planted in nursery pots (Solanaceae), native to South America, is considered a problematic (20 cm tall, 21 cm diameter) and fertilized with 15 g of a slow re- weed of pastures and natural areas in Florida (Mullahey, 1996). lease fertilizer (19-6-12N-P-K OsmocoteÒ, Scotts-Sierra Horticul- TSA has spread rapidly throughout Florida since it was first tural Products Co., Marysville, OH, USA). All plants were kept reported in 1988 (Mullahey et al., 1993), and the plant is now pres- inside screen cages at the University of Florida, Indian River Re- ent in several other states (e.g., Arizona, Arkansas, Georgia, Tennes- search & Education Center (IRREC) located in Fort Pierce, FL. Plants see, Texas) (National Plant Data Center, 2010). It is a perennial were watered daily using an automatic drip irrigation system. weed characterized by having thorn-like prickles on leaves and stems, which are unpalatable to , but has yellow fruits which 2.2. Data collection in the laboratory are consumed by livestock, aiding in the dispersal of seed (Mullahey et al., 1993). Management practices against this weed G. boliviana pupae were held individually in Petri dishes at 24 °C are typically based on expensive herbicide treatments and mow- with a fresh TSA leaf until adult emergence. After one week, adults ing, costing Florida ranchers an estimated $6.5 to 16 million annu- were sexed by examining the ventral side of the abdomen. Males ally (Thomas, 2007). In addition, TSA harbors several viruses that were distinguished from females by the presence of two rounded can infect solanaceous crops (McGovern et al., 1994a, 1994b, orange testes which were visible through the ventral side of the 1996) and serves as an alternative for key insect pests of cul- abdominal sternites (Medal et al., 2002). Single pairs of adults were tivated crops (e.g., Leptinotarsa decemlineata Say and sexta placed inside cages (60 Â 60 Â 60 cm BugDorms, Bioquip, Rancho L.) (Habeck et al., 1996; Sudbrink et al., 1999). Dominguez, CA, USA) with one TSA plant and maintained at Foreign exploration for natural enemies of TSA was initiated in 24 ± 2 °C and a 14:10 L:D photoperiod (9 replicates total). Plants South America in 1994, and several potential biocontrol agents were inspected daily and all eggs were counted, removed and were identified (Medal et al., 1996). One of the most promising placed inside Petri dishes. If a male died before the female, the agents was a leaf-feeding beetle, Gratiana boliviana Spaeth (Coleop- male was replaced. Lifetime fecundity was recorded for each fe- tera: Chrysomelidae), from Argentina and Paraguay. Host range male. Petri dishes containing fresh eggs (<24 h old) and moist filter studies demonstrated that this beetle is highly host specific, feed- paper were kept inside environmental growth chambers (24 ± 2 °C, ing and completing development only on TSA (Medal et al., 2002). 14:10 L:D photoperiod). TSA leaves were replaced as needed, and G. boliviana was introduced into Florida in 2003, and since then it survival and date of molting were recorded daily. A total of 1247 has successfully established throughout South and Central Florida eggs (<24 h old) were followed until adult emergence or death. (Medal et al., 2006; Overholt et al., 2009). Both larvae and adults feed on TSA foliage by making shot-gun holes in the leaves and eventually defoliating the whole plant (Medal et al., 2002, 2006). 2.3. Data collection in the field Adults of G. boliviana lay single eggs on the leaves and petioles, and the larvae complete five and pupate on the plant foliage Life-history parameters of G. boliviana were studied in Central (Diaz et al., 2008). Beetle populations increase during the spring– and North Florida during the summer and fall of 2010. Experiments summer (May–Oct) and adults enter reproductive during were conducted at two field sites in Saint Lucie County (Central Flor- the coldest months (Nov–April) in Florida (Overholt et al., 2010; ida), and a third site in Madison County (North Florida) (Table 1). Diaz et al., unpublished results). Overwintering adults become Ten potted TSA plants (25 cm tall) were caged with reproductive reproductive between the second and third week of March and adults (15 females and 5 males per cage) of G. boliviana for 24– continue to lay eggs for several weeks resulting in overlapping gen- 48 h to allow oviposition. All plants were inspected for eggs, and erations. Substantial suppression of TSA has occurred in some 30 eggs per plant were marked with a red indelible marker pen areas with density declining as much as 90% within three years (red line next to, but not touching the ) while extra eggs were re- in South and Central Florida, while poor establishment of G. boliv- moved from the plants. TSA plants with eggs were then placed in the iana has been reported in North Florida (Overholt et al., 2009, field where they were enclosed inside a fine mesh cloth cage (46 cm 2010). Several factors may explain the differential impact in North diameter, 70 cm tall) placed over a wire plant support. A 60- and South/Central Florida including abiotic factors and biotic fac- cm vertical zipper in the cage allowed access to the plant. Five TSA tors. For example, three mirid species (Engytatus modesta Distant, plants were randomly assigned to each treatment (open or closed Distant, and Macrolophus sp.) are abundant on cages). For the open cages, the bottom of the mesh was folded such TSA plants, and are known to feed on G. boliviana eggs and larvae that only the top half of the plant was covered. Pots were buried in in the laboratory and field (R. Diaz, personal observation). the ground, straw mulch was placed around the pot, and pots were The goal of this study was to examine the mortality factors that separated by 1.5 m and arranged in two rows (five plants per row). influence the population dynamics of G. boliviana in Central and Data loggers measured temperature and humidity inside open and North Florida. The specific objectives were to: (1) compare biotic closed cages during the experiment. The experiment was replicated and abiotic mortality between the laboratory and field, (2) com- twice at the Central Florida sites and once at the North Florida site. 356 V. Manrique et al. / Biological Control 59 (2011) 354–360

Table 1 Description of three Florida field sites used for Gratiana boliviana horizontal life table studies. CF refers to Central Florida, NF refers to North Florida.

Location GPS coordinates description Site 1 CF Privately owned ranch, Fort Pierce, N 27.3455 Mosaic of open pastures and wooded areas. TSA density has been reduced by G. boliviana St. Lucie Co. W 80.4966 Site 2 CF IRREC, Fort Pierce, St. Lucie Co. N 27.4262 Open field at research station next to field crops. Previously used for TSA-G. boliviana studies W 80.4056 Site 3 NF Privately owned ranch Madison, N 30.4025 Mosaic of open pastures and wooded areas. G. boliviana has not reduced the density of TSA Madison Co. W 83.4566

All plants were inspected twice per week, and the total number method was used to estimate variance and bias of estimators in or- of eggs, small larvae (instars I, II, and III), large larvae (instars IV der to quantify uncertainty associated with demographic parame- and V), pupae, and new adults were counted. All newly laid eggs ter estimates (for more details see Maia et al., 2000). and adults were removed from the plants in the open cages. In Longevity of adults was analyzed using two factor analysis of addition, the number of predators found on each plant was re- variance (ANOVA) with set-up (laboratory vs. field) and gender (fe- corded during each sampling period (Central Florida only). Newly male, male) as factors. Significant interaction was obtained emerged adults were counted and placed together for 10 days on (F1,127 = 6.4, P = 0.01), therefore, single ANOVA’s were used to com- TSA plants in cages to obtain sexually mature adults (9–12 d pre- pare longevity between laboratory and field (P < 0.05). oviposition period, Medal et al., 2003). Two pairs of mature adults were placed in closed field cages (10 cages total, 40 adults in total) 2.6. Calculations of vertical life table parameters and the number of adults alive (females and males) was recorded weekly. All TSA plants were replaced every 3 weeks, and the total Vertical life tables of G. boliviana were constructed using natu- number of eggs per plant was recorded. The experiment was termi- rally occurring overlapping generations at two field sites in Central nated when all adults had died (2 replicates in Central Florida Florida. One TSA infested pasture located in Okeechobee Co. (N only). 27.283, W 81.042) and one in St. Lucie Co. (N 27.363, W 80.560) were sampled every two weeks from June to October 2010. Beetles 2.4. Calculations of horizontal life table parameters were originally released in 2006 at both locations and populations have been established since then. In each pasture, fifteen randomly Life tables were constructed from the laboratory and field data selected TSA plants separated by at least 5 m were sampled along a (open and closed cages). The following parameters were calcu- longitudinal transect. Each plant was cut at the ground level, placed lated: number of individuals alive at the beginning of each stage in a plastic bag and transported to the laboratory for data collection. (l ), the number individuals that died during each age interval x The numbers of eggs, larvae, pupae and adults of G. boliviana were (d ), and the percent of individuals that died during each age inter- x counted on each plant. In addition, all predators and herbivores val (100q ). The causes of mortality were separated into intrinsic x found on plants were recorded. Leaves from each plant were har- factors (laboratory mortality), abiotic factors (difference between vested, oven-dried at 70 °C for one week, and weighed to measure mortality in closed cages and in the laboratory), and biotic factors leaf biomass. The mean development times for eggs (6 d), small lar- (difference between mortality in closed and open cages). Compar- vae (13 d), large larvae (9 d), and pupae (7 d) were obtained from isons of mortality factors were made using the Wilson method for the laboratory horizontal life table, while adult longevity (70 d) two independent proportions (Newcombe, 1998). Significant dif- was obtained from the fertility life tables constructed from the field, ferences between mortality were obtained if the 95% confidence and used to calculate survivorship curves. The relationship of interval (CI) did not include zero, while differences between pro- survival to adult stage and sampling date was compared between portions were significant (P < 0.05) if there was no overlap be- sites using linear regression, and pooled after it was found that tween 95% CIs. the slopes of the regression lines were not different (slope Okeechobee = 0.00039, slope St. Lucie = 0.00025; t = 1.3, P = 0.224). 2.5. Calculations of fertility life table parameters A multiple linear regression with backward elimination of non- significant (P > 0.05) independent variables was performed to Fertility life tables were constructed using data from the labora- examine the relationship between the number of mirid adults tory and the field (open cages). The following demographic param- and the following variables: sampling date, leaf biomass, and num- eters were calculated: net reproductive rate (R0), intrinsic rate of bers of spiders, ants, G. boliviana adults and immature G. boliviana increase (rm), finite rate of increase (k), mean generation time (T) (sum of eggs, larvae and pupae). Sampling date was included in the (weeks) and doubling time (Dt) (weeks). Net reproductiveP rate is model to allow estimation of the effect of leaf biomass and number the total contribution of female offspring per female (Ro = lxmx , of beetles and predators on mirid counts independent of time. The where lx is the proportion of females surviving from birth to exact analysis was performed with PROC REG procedure of SAS (SAS age x and mx is the average number of offspring produced by a fe- Institute, 2008). male in theP interval x to x +1)(Carey, 1993). Intrinsic rate of in- crease [1 = lxmx exp(Àrm)] is the rate of increase of a population assuming the population is in stable age distribution and an unlim- 3. Results ited environment. Finite rate of increase k = exp(rm) is the value that the population must be multiplied by to obtain theP population 3.1. Horizontal life table parameters size in the next time unit. Mean generation time (T = vlxmx/Ro)is the period between births of subsequent generations. Doubling Data from the laboratory environment (constant favorable time (Dt = ln(2)/rm) is the time necessary for a population to conditions with no natural enemies) allowed the construction of double in size. Demographic parameters (R0, rm, k, T and Dt) and a horizontal life table in which the mortality obtained during each their confidence intervals were calculated using the SAS (SAS Insti- life stage was attributed to intrinsic mortality (Table 2). The use of tute, 2008) program provided by Maia et al. (2000). The jackknife open and closed cages in the field allowed mortality measurements V. Manrique et al. / Biological Control 59 (2011) 354–360 357 of abiotic factors (climatic conditions) and biotic factors (preda- Table 3 tion) (Tables 3 and 4). Mean temperature in the field varied Horizontal life table for Gratiana boliviana constructed from open and closed cage data from two sites in Central Florida (mean ± SE). between 26 and 30 °C, and no differences were detected between open and closed cages (F1,9181 = 2.1, P = 0.14). Two closed cages Parameter Egg L1–L3 L4–L5 Egg-adult were discarded because plants desiccated at site 1, and one closed Closed cage at site 2 (Central Florida) was discarded due to the presence of lx 127.5 ± 14.4 84.0 ± 10.3 73.2 ± 8.5 58.2 ± 3.4 127.5 ± 14.4 the red imported fire ant, Solenopsis invicta Buren. The egg stage dx 46.5 ± 9.6 13.2 ± 2.6 12.2 ± 6.3 3.0 ± 1.7 72.2 ± 13.3 100q 33.5 ± 4.8 12.3 ± 4.3 18.8 ± 4.9 4.9 ± 2.6 55.2 ± 4.5 experienced the highest mortality in both laboratory and in the x field (Tables 2–4). As expected, mortality was much higher in open Open l 150.0 ± 0 47.5 ± 7.3 33.0 ± 3.9 19.5 ± 3.2 150.0 ± 0 cages compared to closed cages (Tables 3 and 4). There was no dif- x dx 102.5 ± 7.3 14.5 ± 8.4 14.0 ± 2.3 1.7 ± 0.5 132.2 ± 2.8 ference in overall mortality between North and Central Florida 100qx 68.3 ± 4.9 26.4 ± 11.9 42.9 ± 7.6 8.5 ± 1.1 88.1 ± 1.8 (Wilson Independent Proportion Test, P > 0.05). However, egg mor- tality due to predation (open–closed cages) was greater in Central lx = number of individuals alive at the beginning of each stage, dx = number indi- viduals that died during each age interval, 100qx = percent of individuals that died Florida (35%) (Table 3), while larval predation was higher in North during each age interval. Florida (40%) (Table 4). Time of development from egg to adult was similar in the field (30.6 ± 0.9 d) and the laboratory at 24 °C

(33.6 ± 0.7 d) (F1,56 = 1.61, P = 0.21). Overall mortalities from egg to adult were as follows: 85% in Table 4 Horizontal life table for Gratiana boliviana constructed from open and closed cage data open cages > 50% in closed cages > 40% in the laboratory. Indepen- from one site in North Florida. dent tests of proportions showed that mortalities due to biotic fac- tors and intrinsic factors were similar, and both were higher than Parameter Egg L1–L3 L4–L5 Pupa Egg-adult mortality due to abiotic factors (Table 5). A complex of three mirid Closed species (E. modesta, T. notatus, and Macrolopus sp.) comprised 95% lx 150 87 86 86 150 d 63 1 0 1 65 of the predators found during the field studies. Repeated measures x 100qx 42 1 0 1 43 analysis of mirid densities showed no differences between repli- Open cates in Central Florida (F3,83 = 1.82, P = 0.16), but differences were lx 150 81 52 31 150 detected among months (F9,83 = 6.9, df = 9, 83, P < 0.0001). Overall, dx 69 29 21 3 122 mirid densities tended to increase during the season and were 100qx 46 36 40 10 81 highest 21 days after initiating the experiment (Fig. 1). lx = number of individuals alive at the beginning of each stage, dx = number indi-

viduals that died during each age interval, 100qx = percent of individuals that died 3.2. Fertility life table parameters during each age interval.

Fertility life tables of G. boliviana were constructed using data from the laboratory and field (Central Florida only). Adults were followed from July to February in the field, but most adults died be- Table 5 Comparison of mortality factors of Gratiana boliviana reared under different condi- fore the winter started. Two adults turned brown and migrated to tions using the Wilson method for two independent proportions (Newcombe, 1998). the bottom of the cage, indicative of entering reproductive dia- pause (R. Diaz, unpublished results), and were excluded from data Comparisons Proportion difference 95% CI Closed vs. open cage 0.317 (biotic factors) 0.214–0.414 analysis. The net reproductive rate (R0), intrinsic rate of increase Lab vs. open cage 0.481 (abiotic + biotic) 0.413–0.531 (rm), and finite rate of increase (k) were significantly greater in Lab vs. closed cage 0.163 (abiotic factors) 0.072–0.250 the laboratory compared to the field, while no differences were de- Laboratory onlya 0.399 (intrinsic mortality) 0.372–0.427 tected in mean generation time (T)(Table 6). In contrast, doubling a time (Dt) was higher in the field than the laboratory (Table 6). Life- Simple proportion. time fecundity was similar in the laboratory (271.86 ± 43.74 eggs) and the field (232.86 ± 37.52 eggs) (F = 0.38, df = 1, 27, P = 0.54). Fe- males lived longer in the field (80.4 ± 7.5 d) compared to the labo- (Table 7). A negative exponential model fitted the relationship be- ratory (60.8 ± 5.2 d) (F = 4.4, df = 1, 75, P = 0.03), but no differences tween the number of individuals per life stage and cumulative were detected for male longevity (ca.70 d) (F1,52 = 2.6, P = 0.1). In developmental time (Fig. 2). Only 5.1% of eggs reached the adult the field, females lived longer than males (F1,77 = 4.5, P = 0.03). stage. G. boliviana was the most common herbivore found on TSA 3.3. Vertical life table parameters plants followed by the Leptinotarsa decem- lineata (Say) (Coleoptera: Chrysomelidae) (Fig. 3). Other herbivores A total of 1909 eggs, 2148 larvae, 182 pupae and 1109 adults included larvae, stinkbugs () and long were collected during the five month sampling period. Calculated horn beetles (Cerambycidae). Mirids (E. modesta, T. notatus, and mortality was particularly high for eggs and large larvae Macrolopus sp.) were the most common predators on TSA plants followed by spiders (Fig. 3). Other predators included assassin bugs (), red fire ants (Formicidae), big-eyed bugs (: Table 2 Geocorinae), ground beetles (Carabidae) and tree crickets Horizontal life table for Gratiana boliviana constructed from laboratory data. (Grillidae) (Fig. 3). Sampling date, number of ants and G. boliviana Parameter Egg L1 L2 L3 L4 L5 Pupa Egg-adult adults were not significantly related to mirid numbers (P > 0.05)

lx 1247 966 918 870 846 821 797 1247 and therefore removed from the multiple regression model. The dx 281 48 48 24 24 24 48 498 number of adult mirids per plant was positively related to plant 100qx 22.5 5.0 5.3 2.8 2.9 6.1 6.1 40 biomass (partial correlation coefficient = 0.45, P < 0.0001), the lx = number of individuals alive at the beginning of each stage, dx = number indi- number of immature G. boliviana (eggs, larvae and pupae) (partial viduals that died during each age interval, and 100qx = percent of individuals that correlation coefficient = 0.03, P = 0.004) and spiders (partial corre- died during each age interval. lation coefficient = 0.03, P = 0.0001). 358 V. Manrique et al. / Biological Control 59 (2011) 354–360

Fig. 1. Number of mirids (mean ± SE) in open cages during Gratiana boliviana life table studies in Central Florida.

Table 6 Jackknife estimates of fertility life tables for Gratiana boliviana under laboratory and field conditions. Different letters indicate significant differences between laboratory and field parameters (no overlapping 95% CI correspond to significant differences, P < 0.05).

Parameter True calculation Jacknife estimate 95% CI

R0 Laboratory 79.489 a 79.489 49.996–108.981 Field 14.089 b 14.089 11.477–16.703

rm Laboratory 0.430 a 0.423 0.352–0.493 Field 0.319 b 0.319 0.307–0.330 k Laboratory 1.537 a 1.525 1.416–1.635 Field 1.377 b 1.376 1.360–1.393 T Laboratory 10.175 a 10.310 8.038–12.582 Field 8.260 a 8.288 7.580–8.990 Dt Laboratory 1.611 b 1.631 1.377–1.885 Field 2.160 a 2.168 2.088–2.248

Fig. 2. Survivorship curve of Gratiana boliviana instars (eggs, larvae small and large, pupae and adults), in Central Florida pastures. Open circles are observed values and 4. Discussion the line is the fitted model.

Life table studies provide measurements of stage-specific mortality factors and population parameters that are critical to feeding biology of herbivores influence susceptibility to natural understanding insect population dynamics in the field (Royama, enemies (Cornell and Hawkins, 1995; Hawkins et al., 1997). Higher 1981; Carey, 1993, 2001; Kuhar et al., 2002; Zanuncio et al., mortality rates by predators were found on ectophagous herbi- 2006; Yang et al., 2010). Surprisingly, this study is one of very vores, such as G. boliviana, compared to endophagous (Hawkins few that has examined life table parameters of a weed biological et al., 1997). A complex of three mirid species (E. modesta, T. not- control agent in the field (Briese, 1986; Jayanth and Bali, 1994). atus, and Macrolopus sp.) was the most abundant group of preda- G. boliviana immature stages, in particular the egg and larval tors found in the field, and they are known to feed on G. stages, suffered significant mortality with only 5–19% survival to boliviana eggs and larvae (R. Diaz, personal observation). These adulthood. Horizontal life tables were constructed using laboratory mirids are sap-feeders as well as predators in several crop systems and field data, which permitted the separation of different mortal- (Alomar et al., 1991; Sampson and King, 1996; Letourneau and ity factors. The intrinsic mortality (laboratory data) and biotic fac- Goldstein, 2001; Wheeler, 2001). This study was conducted after tors (predation) together accounted for 75% of the mortality of G. 4–6 years of field releases in Florida and could serve as a reference boliviana immature stages, while abiotic factors were less impor- for comparison of mortality factors evaluated in the future. tant. According to Cornell and Hawkins (1995), attack by natural Several studies have shown that G. boliviana significantly reduces enemies is the major source of mortality inflicted on immature TSA growth and reproduction, and has resulted in a decline of TSA stages of insect herbivores. Moreover, developmental stage and populations in Central Florida (Overholt et al., 2009, 2010; Medal

Table 7 Stage-specific survivorship of Gratiana boliviana in Central Florida pastures.

Stage Field counts (n) Duration of each stage (d) Number of individuals (lx) Number dying (dx) Proportion dying (qx) Eggs 212.1 6 35.4 22.3 0.6 L1–L3 169.4 13 13.0 5.3 0.4 L4–L5 69.2 9 7.7 4.8 0.6 Pupae 20.2 7 2.9 1.1 0.4 Adults 123.2 70 1.8 Survival to adult = 5.1% V. Manrique et al. / Biological Control 59 (2011) 354–360 359

ratory compared to the field, and doubling time was longer in the field. These results were expected since laboratory data were ob- tained under putatively ideal conditions for population growth (favorable environmental conditions and absence of natural ene- mies). The intrinsic rate of increase for G. boliviana was converted to days for comparison with other studies on chrysomelid beetles, and was higher (0.06 in laboratory) than Octodonta nipae (Maulik)

(rm = 0.01–0.03 in the laboratory) (Hou and Weng, 2010), but lower than Gastrophysa viridula (DeGeer) (rm = 0.1 in the laboratory) (Hon- ek et al., 2003). Overholt et al. (2010) derived rm from sampling a nat- urally occurring population of G. boliviana over three years and found that it varied between À0.13 and 0.15 depending on the time of year, which agrees well with the current finding. It was quite sur- prising that female G. boliviana lived longer in the field than in the laboratory. One likely explanation is that some adults in the field en- tered diapause in the late fall (November–December), resulting in increased longevity in the field. Vertical life tables are often employed for organisms with over- lapping generations, and assume stable age distribution and con- stant mortality factors acting on the populations (Bellows et al., 1992; Carey, 1993, 2001; Munga et al., 2007). The construction of vertical life tables provided estimates of age-specific mortality of G. boliviana at a larger population scale than our open and closed cage study. High mortality rates were found for all immature stages, particularly eggs and larvae, with an overall survival to adulthood of 5%, which falls in the lower survival range of other introduced insect herbivores (see review by Cornell and Hawkins, 1995). A possible explanation for the higher survival (12–19%) Fig. 3. Relative abundance (percent) of herbivores and predators found found in open cages was cage design; the half cages may have neg- on tropical soda apple in Central Florida. Note: CPB refers to Colorado potato beetle. atively influenced predator searching behavior and provided some protection from abiotic factors (e.g., rain and wind). and Cuda, 2010). However, poor establishment of G. boliviana has As indicated previously, mirids (76%) were the most common been reported in North Florida (Overholt et al., 2009). According to predators found on TSA plants throughout this field study. More- this study, the overall mortality to adulthood (85%) was similar be- over, the number of adult mirids per plant was positively corre- tween Central and North Florida. Although predation rates of imma- lated with plant biomass, G. boliviana immature stages and ture stages varied between locations with higher egg predation in spiders. These mirid species are known to be both phytophagous Central Florida (35%) and higher larval predation in North Florida and predaceous (Letourneau and Goldstein, 2001; Wheeler, (40%), overall mortality was similar (88% and 81% in Central and 2001), which may explain the positive correlation with plant bio- North Florida, respectively). As mortality was approximately equal mass and beetle densities. Both mirids and spiders may share a in North and Central Florida, other factors must be responsible for common prey such as G. boliviana larvae, which resulted in their the poor performance of G. boliviana in North Florida. Differential increased abundance in the field. In addition, G. boliviana was the survival during winter is a likely explanation, either through direct most common herbivore found during this study, supporting evi- effects of cold on beetles, or indirect effects through the host plant. dence from previous studies (Overholt et al., 2009, 2010; Medal Adult G. boliviana enter reproductive diapause at low temperatures and Cuda, 2010) indicating that this biocontrol agent exerts a and short photoperiods during winter months, which coincides with strong top-down impact on TSA populations. decreased host availability. Above ground parts of TSA die back in re- In summary, this study provides a better understanding of the dif- sponse to freezing temperatures, but roots survive and resprout ferent mortality factors affecting the population dynamics of G. boliv- when temperatures increase (Patterson et al., 1997). Based on our iana, a biocontrol agent of TSA in Florida. Even though significant field observations in Central and North Florida, TSA is absent for mortality to adulthood (81–95%) was found in the field, positive longer periods in North Florida than further south, which may result growth rates result in a reduction of TSA densities in the field (Over- in a mismatch in the seasonal phenologies of the beetle and its host- holt et al., 2010). In addition, herbivore damage in combination with plant. Beetle diapause may terminate in North Florida before TSA has other factors such as plant diseases may further reduce plant perfor- had sufficient time to recuperate from repeated freezes, leading to mance (see Caesar, 2005). TSA is known to host numerous plant starvation of beetles. Additionally, host-plant quality may vary be- viruses including the TSA mosaic virus (TSAMV) that affects plant tween locations, and this has been shown to affect performance of growth and is widespread in Florida (Adkins et al., 2007). A complex G. boliviana (Overholt et al., 2009). Additional studies are needed of three mirid species was the most common predator group found on to elucidate the factors which negatively influence beetle perfor- TSA plants, but further studies are needed to investigate individual mance at locations above 29° latitude (Overholt et al., 2009). predation rates of each mirid species on immature stages of G. boliv- Fertility life table studies have frequently been used to estimate iana. Finally, life table studies provide useful measurements of pop- population growth of biological control agents of insect pests, as ulation performance that can provide insight into the establishment they provide insight into the agent’s efficacy and potential for suc- and population growth of introduced biological control agents. cess in controlling populations (Bellows et al., 1992; Medeiros et al., 2000; Kuhar et al., 2002; Zanuncio et al., 2006; Gitau et al., Acknowledgments 2007). Fertility life tables of G. boliviana were constructed using lab- oratory and field data. The net reproductive rate (R0), intrinsic rate of The authors are grateful to Karen Stratman, Ben Anuforum, Joey increase (rm), and finite rate of increase (k) were greater in the labo- Fancher, and Katherine Jordan (University of Florida) for laboratory 360 V. Manrique et al. / Biological Control 59 (2011) 354–360 and field assistance in Central Florida, and to John Mass and Chris McGovern, R.J., Polston, J.E., Mullahey, J.J., 1996. Tropical soda apple (Solanum Albanese (USDA-ARS) in North Florida. Special thanks is extended viarum Dunal): host of tomato, pepper, and tobacco viruses in Florida. In: Mullahey, J.J. (Ed.), Proceedings of the Tropical Soda Apple Symposium. to Mr. Buzz Eaves, Ms. Pat Modine, the Micco Landing Management University of Florida-IFAS, Bartow, Florida, pp. 31–34. Area, and Mr. John Douglas for providing access to their ranches for Medal, J.C., Charudattan, R., Mullahey, J.J., Pitelli, R.A., 1996. An exploratory insect the field studies, and their constant support and interest in biolog- survey of tropical soda apple, Solanum viarum in Brazil and Paraguay. Florida Entomology 79, 70–73. ical control. In addition, Guillermo Logarzo (USDA-ARS SABCL, Medal, J.C., Sudbrink, D., Gandolfo, D., Ohashi, S., Cuda, J.P., 2002. Gratiana boliviana, Argentina) for his advice in statistical analysis. We also thank the a potential biocontrol agent of Solanum viarum: quarantine host-specificity Florida Department of Agriculture and Consumer Services for sup- testing in Florida and field surveys in South America. BioControl 47, 445–461. Medal, J.C., Gandolfo, D., Overholt, W., Stansly, P., Roda, A., Osborne, L., Hibbard, K., plying beetles and financial support for the study. Gaskalla, R., Burns, E., Chong, J., Sellers, B., Hight, S., Cuda, J.P., 2003. Biology of Gratiana boliviana, the First Biocontrol Agent Released to Control Tropical Soda References Apple in the USA. Entomology and Nematology Department, Florida Cooperative Extension Service, Institute of Food and Agricultural Sciences, University of Florida, Florida, . Adkins, S., Kamenova, I., Rosskopf, E.N., Lewandowski, D.J., 2007. Identification and Medal, J., Overholt, W., Stansly, P., Osborne, L., Roda, A., Chong, A.J., Gaskalla, R., characterization of a novel tobamovirus from tropical soda apple in Florida. Burns, E., Hibbard, K., Sellers, B., Gioeli, K., Munyan, S., Gandolfo, D., Hight, S., Plant Disease 91, 287–293. Cuda, J., 2006. Classical Biological Control of Tropical Soda Apple in the USA. Alomar, O., Castañe´, C., Gabarra, R., Arno, J., Arino, J., Albajes, R., 1991. Conservation Entomology and Nematology Department, Florida Cooperative Extension of native mirid bugs for biological control in protected and outdoor tomato Service, Institute of Food and Agricultural Sciences, University of Florida, crops. IOBC/WPRS Bulletin 14, 33–42. Florida, . Bellows, T.S., Van Driesche, R.G., Elkinton, J.S., 1992. Life-table construction and Medal, J., Cuda, J.P., 2010. Establishment and initial impact of the leaf-beetle analysis in the evaluation of natural enemies. Annual Review of Entomology 37, Gratiana boliviana (Chrysomelidae), First biocontrol agent released against 587–614. tropical soda apple in Florida. Florida Entomology 93, 493–500. Briese, D.T., 1986. Factors affecting the establishment and survival of Anaitis Medeiros, R.S., Ramalho, F.S., Lemos, W.P., Zanuncio, J.C., 2000. Age-dependent efformata (Lepidoptera: Geometridae) introduced into Australia for the fecundity and lifefertility tables for Podisus nigrispinus (Dallas) (Het., biological control of St. John’s Wort, Hypericum perforatum. II. Field trials. Pentatomidae). Journal of Applied Entomology 124, 319–324. Journal of Applied Ecology 23, 821–839. Mukherjee, A., Christman, M.C., Overholt, W.A., Cuda, J.P. 2011. Prioritizing areas in Byrne, M. J., Coetzee, J., McConnachie, A.J., Parasram, W., Hill, M.P. 2003. Predicting the native range of Hygrophila for surveys to collect biological control agents. climate compatibility of biological control agents in their region of introduction. Biological control. doi:10.1016/j.biocontrol.2010.11.006. In: Cullen, J.M., Briese, D.T., Kriticos, D.J., Londsdale, W.M., Morin, L., Scott, J.K. Mullahey, J.J., Nee, M., Wunderlin, R.P., Delaney, K.R., 1993. Tropical soda apple (Solanum (Eds.), Proceedings of the XI International Symposium on Biological Control of viarum): a new weed threat in subtropical regions. Weed Technology 7, 783–786. Weeds, Australia, pp. 28–35. Mullahey, J.J., 1996. Tropical soda apple (Solanum viarum), a biological pollutant Caesar, A., 2005. Melding ecology, classical weed biocontrol, and plant microbial threatening Florida. Castanea 61, 255–260. ecology can inform improved practices in controlling invasive plant species. Munga, S., Minakawa, N., Zhou, G., Githeko, A.K., Yan, G., 2007. Survivorship of Biological Control 35, 240–246. immature stages of Anopheles gambiae s.l. (Diptera: Culicidae) in natural Carey, J.R., 1993. Applied Demography for Biologists with Special Emphasis on in western Kenya highlands. Journal of Medical Entomology 44, 758–764. Insects. Oxford University Press, New York, p. 206. National Plant Data Center, 2010. The Plants Database. United States Department of Carey, J.R., 2001. Insect biodemography. Annual Review of Entomology 46, 79–110. Agriculture, Natural Resources Conservation Services, . Cornell, H.V., Hawkins, B.A., 1995. Survival patterns and mortality sources of herbivorous Newcombe, R.G., 1998. Interval estimation for the difference between independent insects: some demographic trends. American Naturalist 145, 563–593. proportions: comparison of eleven methods. Statistics in Medicine 17, 873–890. Diaz, R., Overholt, W.A., Samayoa, A., Sosa, F., Cordeau, D., Medal, J., 2008. Overholt, W.A., Diaz, R., Hibbard, K., Roda, A., Amalin, D., Fox, A., Hight, S.D., Medal, Temperature-dependent development, cold tolerance and potential J., Stansly, P., Carlisle, B., Walter, J., Hogue, P., Gary, L., Wiggins, L., Kirby, C., distribution of Gratiana boliviana (Coleoptera: Chrysomelidae), a biological Crawford, S., 2009. Releases, distribution and abundance of Gratiana boliviana control agent of tropical soda apple, Solanum viarum (Solanaceae). Biocontrol (Coleoptera: Chrysomelidae), a biological control agent of tropical soda apple Science and Technology 18, 193–207. (Solanum viarum, Solanaceae), in Florida. Florida Entomology 92, 450–457. Ghosheh, H.Z., 2005. Constraints in implementing biological weed control: a review. Overholt, W.A., Diaz, R., Markle, L., Medal, J.C., 2010. The effect of Gratiana boliviana Weed Biology Management 5, 83–92. (Coleoptera: Chrysomelidae) herbivory on growth and population density of Gitau, C.W., Ngi-Song, A.J., Otieno, S.A., Overholt, W.A., 2007. Host preference of tropical soda apple (Solanum viarum) in Florida. Biocontrol Science and Xanthopimpla stemmator (Hymenoptera: Ichneumonidae) and its reproductive Technology 20, 791–807. performance on selected African lepidopteran stem borers. Biocontrol Science Patterson, D.T., McGowan, M., Mullahey, J.J., Westbrooks, R.G., 1997. Effects of and Technology 17, 499–511. temperature and photoperiod on tropical soda apple (Solanum viarum Dunal) Goeden, R.D., Louda, S.M., 1976. Biotic interference with insects imported for weed and its potential range in the US. Weed Science 45, 404–408. control. Annual Review of Entomology 21, 325–342. Royama, T., 1981. Evaluation of mortality factors in insect life table analysis. Habeck, D.H., Medal, J.C., Cuda, J.P., 1996. Biological control of tropical soda apple. Ecological Monographs 51, 495–505. In: Mullahey, J.J. (Ed.), Proceedings of the Tropical Soda Apple Symposium. Sampson, A., King, V., 1996. Macrolophus caliginosus, field establishment and pest University of Florida-IFAS, Bartow, Florida, pp. 73–78. control effect in protected tomatoes. IOBC/WPRS Bulletin 19, 143–146. Hawkins, B.A., Cornell, H.B., Hochberg, M.E., 1997. Predators, , and SAS Institute, 2008. SAS User’s Manual, Version 9.2. SAS Institute, Cary, NC. pathogens as mortality agents in phytophagous insect populations. Ecology 78, Semple, J.L., Forno, I.W., 1987. Native parasitoids and pathogens attacking Samea 2145–2152. multiplicalis Guenée (Lepidoptera: Pyralidae) in Queensland. Journal of the Honek, A., Jarosik, V., Martinkova, Z., 2003. Effect of temperature on development Australian Entomological Society 26, 365–366. and reproduction in Gastrophysa viridula (Coleoptera: Chrysomelidae). Senaratne, K.A.D.W., Palmer, W.A., Sutherst,, R.W., 2006. Use of CLIMEX modeling to European Journal of Entomology 100, 295–300. identify prospective areas for exploration to find new biological control agents Hou, Y., Weng, Z., 2010. Temperature-dependent development and life table for prickly acacia. Australian Journal of Entomology 45, 298–302. parameters of Octodonta nipae (Coleoptera: Chrysomelidae). Environmental Sudbrink, J.R.D.L., Snodgrass, G.L., Bryson, C.T., Medal, J.C., Cuda, J.P., Gandolfo, D., Entomology 39, 1676–1684. 1999. associated with tropical soda apple, Solanum viarum in the Jayanth, K.P., Bali, G., 1994. Life table of the parthenium beetle, Zygogramma Southeastern USA. In: Spencer, N.R. (Ed.), Proceedings X International bicolorata Pallister (Coleoptera: Chrysomelidae) in Bangalore, India. Insect Symposium on Biological Control of Weeds, 4–9 July 1999. USDA-ARS/ Science and its Application 15, 19–23. Montana State University, Bozeman, MT, pp. 247–248. Kuhar, T.P., Wright, M.G., Hoffmann, M.P., Chenus, S.A., 2002. Life table studies of Sutherst, R.W., Maywald, G.E., 1985. A computerized system for matching climates European corn borer (Lepidoptera: Crambidae) with and without inoculative in ecology. Agriculture, Ecosystems and Environment 13, 281–299. releases of Trichogramma ostriniae (Hymenoptera: Trichogrammatidae). Thomas, M., 2007. Impact of tropical soda apple on Florida’s grazing land. The Environmental Entomology 31, 482–489. Florida Cattleman’s and Livestock Journal 71, 33. Letourneau, D.K., Goldstein, B., 2001. Pest damage and arthropod community Wheeler Jr., A.G., 2001. Biology of the plant bugs (: ): pests structure in organic vs conventional tomato production in California. Journal of predators opportunists. Cornell University Press, p. 507. Applied Ecology 38, 557–570. Yang, X.B., Zhang, Y.M., Hua, L., Liu, T.X., 2010. Life history and life tables of Maia, A.H.N., Luiz, A.J.B., Campanhola, C., 2000. Statistical inference on associated Bactericera cockerelli (Hemiptera: Psyllidae) on potato under laboratory and fertility life table parameters using Jackknife technique: computational aspects. field conditions in the Lower River Grande Valley of Texas. Journal of Economic Journal of Economic Entomology 93, 511–518. Entomology 103, 1729–1734. McGovern, R.J., Polston, J.E., Danyluk, G.M., Heibert, E., Abouzid, A.M., Stansly, P.A., Zanuncio, J.C., Lemos, W.P., Lacerda, M.C., Zanuncio, T.W., Serrao, J.E., Bauce, E., 1994a. Identification of a natural weed host of tomato mottle geminivirus in 2006. Age-dependent fecundity and fertility life tables of the predator Florida. Plant Disease 78, 1102–1106. Brontocoris tabidus (: Pentatomidae) under field conditions. McGovern, R.J., Polston, J.E., Mullahey, J.J., 1994b. Solanum viarum: weed reservoir of Journal of Economic Entomology 99, 401–407. plant viruses in Florida. International Journal of Pest Management 40, 270–273.