Survival and population growth of swordtail varieties were examined across increasing community complexity in experimental outdoor ponds (each treatment, n = 4). Treatments included control (swordtails only), simple communities with swordtails and a poeciliid congener (eastern mosquitofish; Gambusia holbrooki; hereafter low) to increasing complexity with weak piscivores (swordtails + mosquitofish + sunfishes Lepomis spp.; medium) and a strong piscivore (swordtails + mosquitofish + sunfishes + largemouth bass Micropterus salmoides; high). This additive sequence of treatments is similar to what a non-native fish would experience when moving from an aquaculture effluent to more natural environments at increasing distance (Hill and Tuckett 2018). Thus, our goal was to mimic natural systems by including species commonly fund near aquaculture effluents and not simply increase diversity as would be expected from a study examining diversity-invasibility.
Pond preparation began in early August 2014 by applying an aquatic herbicide, to ensure ponds began with a similar vegetative community. A fish toxicant was applied on 08-Aug-2014 using recommended dosages to ensure all fish were removed prior to experimentation (Finlayson et al. 2010; Lazur et al. 2013). Ponds were then filled with aerated well water and also groundwater infiltration. Mosquitofish were stocked on 20-Aug-2014 (1000 individuals of various size per pond). A subsample of 100 individuals from the total stock were individually weighed and measured prior to stocking; mosquitofish length (total length; TL) averaged 29.57 mm (SD = 6.93) and weight averaged 0.35 g (SD = 0.29). Average pond size (wetted area) was 168 m2; thus, mosquitofish were stocked at an average density of 6.05 / m2, a lower density than local streams and ditches and lower than that used in a previous experimental study (Thompson et al. 2012). However, mosquitofish typically utilize vegetative habitats, which will be more common around the pond perimeter (Moyle and Nichols 1973; Pyke 2008). Pond size and the resulting mosquitofish beginning density was similar among treatment groups (pond size: F3,12 = 0.537; P = 0.666; density: F2,9 = 0.583; P = 0.638). Ponds lacking a top predator (medium treatments) were stocked with ten bluegills (Lepomis macrochirus), one redbreast sunfish (Lepomis auritus), and one spotted sunfish (Lepomis punctatus), approximating the density and proportion of sunfishes near sources of production (Hill and Tuckett 2018). Ponds with the top predator (high), received the same compliment of bluegill, redbreast sunfish, and spotted sunfish, but also received two largemouth bass. Initial sunfish size (TL) was similar among ponds (nested ANOVA, pond nested within treatment; F6,88 = 0.784; P = 0.533) and between medium (average = 127.46 mm) and high treatments (mean = 123.94 mm; F1,88 = 0.530; P = 0.784). Similarly, mean bass TL among ponds varied from 258 to 304, but was similar among the four ponds (F3,4 = 0.365; P = 0.783).
Periodic propagule pressure was manipulated using weekly stockings of swordtails beginning on 28-Aug-2014; stocking continued for an additional 13 weeks for a total of 14 stocking dates. Twenty swordtails were stocked weekly, and each pond received 280 total swordtails. The proportion of stocked swordtails mirrored trade data and presumably escapement from aquaculture facilities. The top six varieties sold by a local wholesaler were the velvet wag (33.7%), velvet (21.9), neon (18.6), assorted (17.2; not classified to variety), pineapple (14.3), and wild-type variety (11.5). Velvet wag and velvet varieties were unavailable from producers during the required timeframe and were replaced with similarly vibrant varieties, marigold and painted. Seven marigold, four painted, four neon, three pineapple, and two wild type were stocked weekly in each pond (Fig. 1). Thus, propagule pressure was varied among varieties but was held constant among all treatments. A subsample of ten fish from the five varieties were measured for TL. Fish varied in size from 35 to 60 mm (mean = 44.84 ± 6.57 SD) and size was similar among varieties (F4,45 = 1.002; P = 0.417). All swordtail varieties, including the wild-type variety, were sourced from the same local producer using similar rearing protocols. Thus, the wild-type experienced similar rearing conditions as the more colorful varieties but was not deliberately bred by the producer for color patterns. Water quality parameters and temperature were monitored weekly, but always stayed within acceptable levels (Online Resource).
Ponds were randomly harvested over a three-day period, from 01-Dec-2014 to 03-Dec-2014, by first reducing the pool volume ~ 75% by pumping water through a screened intake; pond pumping concentrates fish and aids in capture. Fishes were captured using a fine-mesh seine net (1.59 mm) with three passes. A 1.59 mm mesh seine was used because it captures all but the smallest mosquitofish and swordtails. All captured fish except bass and sunfish were euthanized following institutional protocols and subsequently frozen for later analysis. Adult bass and sunfish were measured (TL) and removed to an on-site detention pond. Relative growth in TL was calculated for bass (relative growth = TL final – TL initial / TL initial). Juvenile sunfish were retained, weighed, and measured. Because of the large numbers of mosquitofish, up to 100 individuals for each pond were subsampled for length (TL) and weight (g); the remaining individuals were batch weighed to estimate numbers at the conclusion of the trial. All swordtails were measured, identified to variety, and sex was determined based on the presence of a gonopodium.
Predation Experiments
Two separate predation experiments were performed to determine bass prey selection: 1) between the swordtail and mosquitofish and 2) among swordtail varieties. Experiment #1 was conducted in an indoor re-circulating tank system. The tanks were concrete (221 cm x 79 cm x 58 cm) with a water depth of 30 cm. Water parameters fell within acceptable levels (Online Resource). All tanks had artificial vegetation as a refuge (645 stems m− 2) which consisted of black plastic strips tied to a plastic screen that covered 50% of the tank bottom (Savino and Stein 1989). Twenty each of marigold swordtails (TL = 32 ± 6 mm) and mosquitofish (TL = 33 ± 4 mm) were stocked concurrently into the tanks. A low density of mosquitofish was used because previous research showed that effects of mosquitofish (i.e., fin nipping of ornamental fish) are density-dependent (Thompson et al. 2012). Prey were stocked concurrently to provide a similar acclimation period but also because the effects of mosquitofish can be delayed when stocked concurrently with other small-bodied fishes (Thompson et al. 2012). A single bass per tank (TL = 235 ± 28 mm) was stocked 1 hour after the prey. Trials (n = 12) lasted 24 hours. All fish were counted at the end of the trials.
Experiment #2 used the same experimental system exhibiting similar water quality parameters. Ten each of marigold, neon, painted, pineapple, and wild-type swordtails were stocked concurrently into the tanks. Swordtails for experiment #2 were taken from the same source as the pond study and thus exhibited similar size variation. Bass (TL = 317 ± 35 mm) were stocked as above and included eight trials, each lasting 24 hours. Bass in experiment #2 were not used in experiment #1. All swordtails were counted at the conclusion of the trials. For both experiments the “case 2” approximation of Chesson’s α was used as a selection index (Chesson 1983; see Hill et al. 2011), which accounts for the depletion of prey through consumption without replacement:
where m is the number of prey categories, ri0 is the number of prey type, i is the diet, and ni0 is the number of prey type i in the environment. Values of α can range between 0 indicating complete avoidance and 1 indicate complete preference. Random prey selection in this experiment (1/m) equaled 0.5 for experiment #1 and 0.2 for experiment #2.
Data Analysis
Our general approach for the analysis of count data was to first fit the data to the Poisson distribution. Alternant distributions were applied to overdispersed data (quasi-Poisson; negative binomial) and the model fit was evaluated using multiple criteria. For the pond study, we first analyzed counts of fish within the artificial community. The number of mosquitofish at the conclusion of the experiment with a GLM (Poisson distribution), but our data were overdispersed (McCullagh and Nelder 1989), whereupon we fit the data to the negative binomial distribution. Mosquitofish size (logarithm transformed TL; total n = 1296) was analyzed with nested GLM (normal distribution; identity link) with factors treatment (low, medium, and high) and pond nested within treatment. Counts of surviving sunfish (pooled among species) and also young of the year (reproduction in the ponds) were both analyzed with GLM (Poisson distribution). Similar to the mosquitofish counts, swordtail counts across levels of biotic resistance were analyzed with GLM (negative binomial distribution).
For the analysis of variation in success of the five varieties across pond treatments, some individuals could not be identified to a specific variety. Most of these individuals were found in control ponds (324 of 339) due to reproduction, which could have been due to hybridization between varieties. Varietal survival (square root transformed # remaining/# stocked) was analyzed with a mixed-model ANOVA with fixed effects treatment, variety, and interaction. Pond was included as a random factor nested within treatment. The control treatment was analyzed separate from the other three treatments because of dramatic differences in variance and to streamline the analysis for the fixed effect, variety. For the control treatment we analyzed both the relative success (square root transformed # remaining/# stocked) and the total counts to determine if propagule pressure increased numbers in the experimental ponds. In both analyses, pond was included as a random factor. Tukey HSD was used to analyze pairwise differences.
Selectivity values were compared using the paired t-test (experiment #1) or within-subjects ANOVA on transformed selectivity values (log (x + 1); experiment #2) because selection values within each tank were not independent. The within-subjects ANOVA included the fixed effect of variety and tank was included as the random effect. For experiment #2, two bass ate more than half of the swordtail prey; one of these bass consumed all 50 swordtails within a 24-hour period, significantly more than expected. Because of a potential saturation effect of predation where the predator could switch to a less desirable or visible prey as the preferred prey becomes scarce (Allen and Greenwood 1988), we excluded both of these trials. Thus, our final replicate number was just six.