
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

61761
10.1038/s41598-024-61761-5
Article
Quantifying the importance of ontogeny and prey type in modeling top-down and bottom-up effects of an ectothermic predator
Nafus Melia G. mnafus@usgs.gov

1
Gray Levi N. levingray82@gmail.com

2
1 https://ror.org/03eg83y24 U.S. Geological Survey, Pacific Island Ecosystems Research Center, Hawai’i Volcanoes National Park, Hawai’i, USA
2 https://ror.org/01y64my43 grid.273335.3 0000 0004 1936 9887 Department of Biological Sciences, University at Buffalo, Buffalo, NY USA
16 9 2024
16 9 2024
2024
14 2160115 5 2023
9 5 2024
© Works prepared by Federal government employees as part of their official duties, cannot be copyrighted. Under the provisions of the Copyright Act Title 17 of the United States Code (U.S.C.), “[c]opyright protection under this title is not available for any work of the United States Government….” (17 U.S.C. § 105). 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Dietary decisions by predators can affect prey abundance and overall food web dynamics. Many predators do not forage on the same prey at the same frequency throughout their lives. Ontogenetic shifts in prey preference are not, however, often accounted for when modeling food web relationships, despite growing literature that suggests that stage specific dietary relationships may be an important consideration when modeling trophic interactions. We investigated the importance of considering size-structure of a predator population with ontogenetic diet shifts in evaluating relationships with prey response using a manipulative experiment with the brown treesnake (Boiga irregularis) in Guam. After removing ~ 40% of the snake population via toxic mammal carrion, we measured the strength of the relationship between snake density and the response of two types of prey (lizards and mammals). We evaluated these relationships based on total population size or division of the population into stage specific size categories based on diet preference predictions. We hypothesized that the density of juvenile snakes would correlate more strongly with lizard detections, while adult snakes would better correlate to rodent detections. We also measured reproductive output following changes in rates of prey detection. As expected by known ontogenetic shifts in dietary preference, explicit stage-based models better predicted shifts in rates of observed prey items than did total predator density for both lizards and mammals. Additionally, rodent detections were predictive of one reproductive pulse from snakes, while lizard detections were not predictive or correlated. Our findings support that consideration of predatory species stage-based dietary preference can be meaningful for understanding food web dynamics, particularly when a predator has a broad diet or one that changes through time.

Keywords

Food web
Ontogenetic diet shifts
Size-structured population
Trophic linkages
Subject terms

Community ecology
Ecosystem ecology
Population dynamics
http://dx.doi.org/10.13039/100020697 Office of Insular Affairs http://dx.doi.org/10.13039/100000005 U.S. Department of Defense issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Prey abundance is a major driver of predator vital rates, such that predators are predicted to have increased per capita birth rates as a function of prey availability1, while prey are anticipated to have reduced fecundity or survival in environments with greater predator abundances2. In traditional Lotka-Volterra dynamics, predators are predicted to overshoot their prey’s carrying capacity, causing prey declines and subsequently predator declines after a lag period3. However, the overall stability of food webs is generally governed by the strength of the relationships between community members, with weaker links contributing to greater stability4. Stage-structured populations, in which juveniles and adults use the same resources, but allocate their time differently, can thus have stabilizing effects on food webs by decreasing linkage strength5.

Based on the effects that population stage structure can exert on food webs, the importance of ontogenetic diet shifts (ODS)—transition of diet over the lifetime of an individual—as a theoretical consideration in food web stability is of growing interest and can be a key component of how trophic communities are structured6–9. Thus, the relative importance of top down control on food webs may be dependent on the stage structure of the predator population10. Moreover, failure to account for shifts in prey resource use based on ontogeny can affect the accuracy of consumption rates and downstream trophic effects11. In a system in which the apex predator undergoes ODS, the stage structure of the population, and thus the relative consumption rates, may determine the relative importance and strength of top-down or bottom-up control of trophic linkages12. Failure to account for the stage structure of the population when measuring or modeling food web dynamics may therefore dampen or mask the detection of direct trophic linkages.

Despite growing attention towards the effects of ODS in regulating trophic interactions, the majority of field studies investigating them have been conducted on aquatic organisms13,14. In reptiles, ontogenetic shifts often co-occur with habitat, physiological, metabolic, or morphometric changes that occur with age or size15–18. Therefore, understanding how predator ontogeny affects trophic dynamics may be particularly important for assessing the community linkages or their strength when evaluating reptiles. Reptiles are recognized to be important ecological effectors19–21. Reptiles can perturb relationships sufficiently to cause disruptive events in food webs, including causing ecological cascades20,22. However, additional studies that evaluate the effect of ontogenetic relationships may refine understanding of trophic linkages involving reptiles. Better understanding of the strength of trophic interactions based on the stage structure of populations will likely improve overall understanding as to the influence of terrestrial ectotherms in their communities.

In the island of Guam (540 km2; also known as Guåhan in CHamoru), introduced brown treesnakes (Boiga irregularis) are documented in their ability to influence native and non-native prey populations17,23–25. There effects on vertebrate loss are linked to trophic cascades in subsequent decades20. Although the potential for this species and other invasive reptiles to modify or destabilize prey communities has been documented, the stage dependencies of relationships that may result from ODS have not been examined21,26,27. In Guam, brown treesnakes have a diet composition that is primarily composed of ectothermic lizards until ~ 1250 mm snout-vent length (SVL), although they begin to incorporate endothermic prey (small mammals and birds) at ~ 650–700 mm SVL17. They are also documented to be reliably attracted to rodent lures at around 900 mm SVL28,29. Based on diets in the wild and subsequent laboratory trials, the attractivity of rodents appears to decline around 1550 mm SVL, with birds becoming an increasingly dominant or preferred food source for individuals when available17,30.

We explored the relationship between demographic stage structure of brown treesnakes (or snakes) on two prey types through a manipulative experiment involving snake removal in a large enclosure. After removing snakes that ate toxic rodent baits, we evaluated whether total density of snakes or ontogenetically segregated densities were better predictors of rat (endothermic) or lizard (ectothermic) prey detection rates, as estimated by sightings-per-unit-effort (SPUE) metrics for prey. Subsequently, we evaluated whether the abundance of rodents or lizards was driven by the abundance of particular snake stage classes to determine whether predator mediated effects were correlated to ontogenetic foraging shifts. We hypothesized that larger snakes would demonstrate tighter linkages to rats and smaller snakes to lizards, while rat detections would be the primary driver of reproductive events.

Methods

Study site

Guam is the largest and southernmost island in the Mariana Island archipelago, with temperatures between 25 and 32 °C. The island undergoes two rainfall seasons: dry and wet, with peak rainfall occurring from July–October and cumulative annual rainfall averaging 2133–2946 mm. Seasonal differences in rainfall can have substantive effects on plant and insect productivity31. We did not, however, monitor plants or insects during our experiment.

We completed this work in a 5-ha experimental snake enclosure located on Andersen Air Force Base, in northern Guam (13.639 N, 144.862 E). The enclosure was designed to prevent the ingress and egress of snakes. Therefore, the population could be accurately counted through mark-recapture, such that births represent the only source of new individuals and deaths the cause of disappearance29. Habitat within the enclosure was mixed native limestone forest and non-native Leucaena leucocephala forest, with a limestone karst substrate32. Within the enclosure, we maintained 27 forest-interior and 2 forest-edge transects that were 220 m long and spaced at ~ 8 m intervals (2-m width). In addition to snakes, the experimental enclosure housed several other vertebrate species, including two nocturnal geckos (Lepidodactylus lugubris [mourning gecko], Hemidactylus frenatus [common house gecko]), two diurnal skinks (Emoia caeruleocauda [Pacific blue-tail skink], Carlia ailanpalai [curious skink], one diurnal mangrove monitor [Varanus tsukamotoi]), and various non-native rodents24,33,34. The specific species of rodent present has not been identified. The enclosure was not closed to lizards or rats and prey immigration and emigration remained possible throughout the course of study.

Overall, rodents have occurred at substantially lower densities in Guam compared to other Mariana Islands, which is a difference thought to be caused by the snakes24. Additionally, there were no native forest birds observed at the study site during the period of this study. The only other potential predators of rats, lizards, and snakes were thus monitor lizards and cats. Although we were not directly monitoring either species, no change was detected in their abundance during nocturnal visual surveys over the period covered herein35.

Snake population manipulation

We defined discrete snake stage classes based on the reliability with which they are attracted to ectothermic (< 700 mm SVL) versus endothermic rodent (≥ 900 mm SVL) prey17,28,29,36. To quantify stage-based effects, we artificially manipulated the population by removing snakes that we expected to consume toxic rodent carrion laced with 80-mg of acetaminophen37. Specific baiting protocols and intervals are described more thoroughly in Nafus et al.38. Briefly, beginning on 31 March 2017, we distributed 594 toxic mouse baits in an 8 × 10 m array at 7 treatment periods. We applied baits in accordance with the Environmental Protection Agency (EPA) pesticide registration label for “Acetaminophen for Brown Treesnake Control” (80-mg uncoated tablets; Reg. No. 56228-34). Outside of exposure to toxicants we did not euthanatize snakes and thus removal was anticipated to affect snakes ≥ 850 mm SVL or those eating rodent carrion. Dead mice are not consumed by small lizards and mammals have much higher tolerances to acetaminophen, limiting the expected effect of baits on other species37. Monitors may consume the baits, which may be lethal for smaller individuals39, but as a large diurnal species these non-target effects should not have overly biased prey counts. Therefore, we did not expect that application of baits to have adverse effects on the other vertebrate species we monitored.

Visual surveys for snakes, lizards, and rats

To quantify the relationship between snake abundance and prey response, we divided surveys into two discrete periods: pre-snake removal (24 October 2016–30 March 2017, N = 67 surveys, mean = 11 surveys/month) and post-snake removal (31 March–31 January 2018, N = 128 surveys, mean = 13 surveys/month). We conducted surveys continuously during the period referenced, with survey effort typically ranging from 2 to 4 per week. Each survey was completed by two teams of two trained observers and we surveyed all 29 transects on a given survey night. Each observer scanned their respective side of a given transect using powerful headlamps (Wilma 3200-lumens, Lupine Lighting Systems, Pennsylvania, USA) per methods described in Christy et al.40. During surveys, observers documented all lizards, small mammals (rats), and snakes. We counted but did not mark rats and lizards along each transect. When we detected a snake, we attempted to capture and mark the individual each time creating a known population.

When we captured unmarked snakes, we clipped ventral scales in a unique pattern using scissors and team members trained by veterinary professionals implanted a passive integrated transponder (PIT) tag into the abdominal cavity. For each captured snake, we scanned for the presence of identification marks (numerical caudal scale clips) and a PIT tag to census the population during each survey. We collected morphometric data, including SVL by gently stretching the snakes along a tape measure each night a snake was captured. We also completed a series of trapping events using a 13 × 13 grid of traps (169 traps), with roughly 16 × 16 m spacing (as per methods described in41) during 18 January 2017 (51) days and 04 April 2017 (9 days). Methods of estimating total abundance of individuals are described in detail in Nafus et al.38,42. Briefly, however, for new snake detections, we estimated birth date by projecting growth using size and sex specific growth rates created for this population during this period to identify the first month each individual was predicted to be < 400 mm SVL, as well as size in a given month when the individual was not detected. We used the last month of detection as the estimate of when a snake was no longer present in the population, which may have underestimated the true duration of time individuals were present. At the close of the study, we thus had monthly estimates of absolute brown treesnake abundance by size and per transect SPUE (individuals sighted per km) for lizards and rats. For the purposes of this manuscript, we did not attempt to differentiate prey by species under the assumption that snakes eat all the species present in the enclosure.

Ethics statement and data availability

We completed all work in accordance with Institutional Animal Care and Use Committee protocols U.S. Geological Survey Fort Collins Science Center #2017-03 and Colorado State University protocol #15‐5892A, such that we performed all methods in accordance with the relevant guidelines and regulations. Methods are reported in accordance with ARRIVE guidelines. Experimental protocols were approved by U.S. Geological Survey. All data used are publicly available for download35,43.

Statistical analysis

The primary purpose of the analyses described herein was to test the relationship between stage-based snake abundance and rodent or lizard SPUE. We, therefore, selected two discrete snake size classes based on evidence for preference of consuming ectotherms (< 700 mm SVL) or reliable attraction to endothermic rodent prey (≥ 900 mm SVL)17,29,36. Based on their highly transitory state for diet we discounted present individuals from 700 to 899 mm SVL by removing individuals within those sizes from all analyses. For each prey dataset, we calculated monthly SPUE averages by averaging the SPUE across all transects surveyed by all observers in any given month (for full description of methods, refer to Nafus43). We opted to combine data across all transects and within months to help minimize bias from spatial, temporal, or observer effects on SPUE33,44.

We used two statistical approaches to evaluate the data: negative binomial regression and time series regression. We completed the negative binomial regression using the ‘MASS’ package in R45. We completed time series regression using the ‘forecast’ package in R46. We used two approaches based on imperfect ability to account for violating specific assumptions by each modelling approach. For example, data were not normally distributed and demonstrated evidence of overdispersion, which violates linear regression assumptions for time series models46. However, given the sequential nature of our study, individual data points were temporally correlated and thus not independent, as would be ideal for a negative binomial regression. We tested model fit using a chi-square goodness of fit test with an acceptable fit determined by P > 0.05. The time series regression demonstrated generally poor predictive power and both models identified the same predictors as significant.

In our first analysis, we tested whether snake density (snakes/ha) predicted rat or lizard SPUE from 24 October 2016 to 30 September 2017 using negative binomial regression with prey as the response and snake density as a fixed effect. Because rat SPUE was very low (< 0.1 rats/km), we multiplied monthly averages by 1000 to create whole numbers. We completed two additional negative binomial regressions with snake density evaluated as two independent fixed effects: < 700 mm SVL or N< 700 or ≥ 900 mm SVL or N≥ 900. One model had rat SPUE as the response and the other lizard SPUE. We also completed time series regressions for these models, using a univariate, autoregressive integrated moving average (ARIMA) approach. We used a 1-year interval subdivided into monthly estimates to forecast time series models, maintaining N< 700 or N≥ 900 as predictors. In order to test whether lizard or rats correlated with snake reproduction, we used lizard and rat SPUE as predictors in a negative binomial regression and time series regression (univariate, ARIMA, with monthly intervals over a 12 month period), for the response of density of juvenile snakes (< 700 mm SVL) after a four month lag period. We selected four months as the lag period to account for a three-month egg gestation period, as well as a one-month lag in detection of juveniles or follicle production. We report the regression estimate (β^), 95% confidence intervals (CI), and significance for negative binomial regressions, and β^ and 95% CI for time series regression by parameter. For all models, we accepted significance at α < 0.05 and we report descriptive statistics as mean ± standard error.

Results

The estimated number of snakes present prior to removal was 117 individuals (23.4 snakes/ha), with N< 700 = 16 and N≥ 900 = 87 (Fig. 1A). From 31 March to 30 September 2017, the population declined to 72 individuals (14.4 snakes/ha), with N< 700 = 22 and N≥ 900 = 37 mm SVL. By 31 January 2018, the number of juvenile snakes (< 700 mm SVL) had roughly doubled from the starting population to N< 700 = 33. The mean SVL of the population decreased from 967 ± 18 in March 2017 to 816 ± 31 mm in January 2018, supporting a shift from predominately endothermically interested snake sizes to ectothermically interested or transitioning into endothermic prey (Fig. 1B).Figure 1 The change in the number (A) and mean population snout-vent length (SVL, B) of brown treesnakes (Boiga irregularis, snakes) in a 5-ha enclosure in Guam proceeding and during experimental removal of snakes willing to eat mouse baits. The number (No.) of snakes is the total individuals in the population, as well as individuals < 700 or ≥ 900 mm SVL. The mean size is the average SVL of snakes present in the population ± the standard error (SE). The red line indicates when removal of snakes susceptible to toxic mouse baits (≥ 850 mm SVL) began.

During removal of snakes susceptible to toxic mouse baits, rat SPUE increased by approximately a factor of 3, from 0.02 ± 0.01 pre-snake removal to 0.06 ± 0.02 rodents/km (Fig. 2A). There was a significant relationship between total snake density and rat SPUE (β^ = − 0.11, 95% CI = − 0.17 to − 0.06, t = − 3.6, P = 0.005). Dividing the total population into stage-specific predictors supported that the relationship between rat SPUE and snake density was, however, stage specific. Based on negative binomial regression, the density of N≥ 900 was negatively correlated with rat SPUE (β^ = − 0.16, 95% CI = − 0.23 to − 0.09, t = − 3.6, P = 0.006); for each decrease of N≥ 900/ha corresponded to an increase in rat SPUE by a factor of 0.9 (Fig. 2B). In contrast, there was no relationship between snakes N< 700/ha and rat SPUE (β^ = − 0.26, 95% CI − 0.10 to 0.02, t = − 1.14, P = 0.28). Time series regression supported that N≥ 900/ha was negatively predictive of rat SPUE, with a CI that did not overlap zero (β^ = − 0.58, 95% CI = − 0.80 to − 0.37). Time series regression indicates that N<700/ha was also negatively correlated with rat SPUE (β^ = − 1.44, 95% CI = − 0.13 to − 2.76), with a wide 95% CI that was approaching but did not cross zero. The combined inclusion of both size classes forecasted a change in rat SPUE that encompassed actual values (Fig. 2C).Figure 2 The density of brown treesnakes (Boiga irregularis) in a 5-ha enclosure in Guam proceeding and during experimental removal, and the monthly mean sightings of rats (prey sightings per unit effort [SPUE]) based on visual surveys. Snakes are divided into two size classes based on ontogenetic interest in rodent (≥ 900 mm snout-vent length [SVL]) or lizard (≤ 700 mm SVL) prey (A). The red line indicates when snake removal began. The negative binomial regression and 95% confidence interval (CI) between rat SPUE relative and snakes interested in rodent prey (≥ 900 mm) was negative (B). The predicted sightings and 95% CI of rats during each time step based on the density of snakes in each size class using time series regression forecasted an increase as snakes ≥ 900 mm declined (C). Actual rat SPUE is represented by black circles.

In contrast to rodents, there was no change in average monthly lizard SPUE in the six months after snake removal (6 ± 1 lizard per km), as compared to the six months during toxic mouse baiting (5 ± 1 lizard per km), although 3 months after snake removal began there was a temporary increase (Fig. 3A). There was no relationship between total snake density and lizard SPUE as estimated by binomial regression (β^ = − 0.01, 95% CI = − 0.03 to 0.01, t = − 0.65, P = 0.53). Dividing the total population into stage-specific predictors, illuminated a significant relationship between lizard SPUE and snake density that was undetectable without stage specific predictors. Based on negative binomial regression density of N<700 was negatively correlated with lizard SPUE (β^ = − 0.23, 95% CI = − 0.34 to − 0.13, t = − 2.80, P = 0.02), such that for each decrease of N< 700/ha corresponded to an increase in lizard SPUE by a factor of 0.8 (Fig. 3B). In contrast, there was no relationship between snakes N≥ 900/ha and lizard SPUE (β^ = − 0.01, 95% CI − 0.04 to 0.01, t = − 1.01, P = 0.33). Time series regression suggest that N<700/ha was also negatively correlated with lizard SPUE (β^ = − 1.47, 95% CI = − 2.38 to − 0.56), with a wide 95% CI that was approaching but did not cross zero. Time series regression supported that N≥ 900/ha was not correlated to lizard CPUE, with a CI that overlapped zero (β^ = − 0.13, 95% CI = − 0.29 to 0.02). The time series regression forecasted a change in lizard SPUE that encompassed actual values (Fig. 2C).Figure 3 The number of brown treesnakes (Boiga irregularis) in a 5 ha enclosure in Guam proceeding and during experimental removal, and the mean monthly sightings of lizards (prey sightings per unit effort [SPUE]), based on visual surveys. Snakes are divided into two size classes based on ontogenetic preferences for rodents (≥ 900 mm snout-vent length [SVL]) or lizards (≤ 700 mm SVL) (A). The negative binomial regression and 95% confidence interval (CI) for the relationship between lizard SPUE relative to snakes interested in only lizard prey (< 700 mm) was negative (B). The predicted and 95% CI of lizard sightings during each time step based on the density of snakes in each size class using time series regression forecasted a decrease as snakes < 700 mm increased (C). Actual lizard SPUE is represented by black circles.

After accounting for a four-month lag, the number of juvenile snakes (< 700 mm SVL) in the population increased to 33 individuals from the 16 originally present (Fig. 4A). Based on negative binomial regression, rat SPUE was a significant positive predictor of juvenile abundance (β^= 0.10, 95% CI = 0.06–0.14, t = 2.73, P = 0.02, Fig. 4B), while lizard SPUE was not correlated (β^= − 0.09, 95% CI = − 0.16 to 0.00, t = − 1.27, P = 0.26). The time series regression poorly fit the data and neither rat (β^= − 0.19, 95% CI = − 0.93 to 0.54) nor lizard SPUE (β^= − 0.17, 95% CI = − 1.39 to 1.05) was correlated with juvenile abundance 4 months later (Fig. 4C).Figure 4 The number of immature brown treesnakes (Boiga irregularis, < 700 mm snout-vent length [SVL]) in a 5 ha enclosure on Guam proceeding and during experimental removal of adults. The monthly mean sightings of lizards and rats (prey sightings per km surveyed [SPUE]) observed during nocturnal visual surveys are depicted relative to the abundance of immature snakes (A). After a 4-month lag, the negative binomial regression relationship for the number of juvenile brown treesnakes (< 700 mm snout-vent length [SVL]) was positively correlated to rodents sighted per km surveyed (rat SPUE) but not lizards (B). The predicted and 95% CI of juvenile snakes during each time step was not related to either prey type (C). Actual juvenile snake count is represented by black circles.

Discussion

In this study, models that accounted for ODS in a top predator illuminated stage specific predator–prey relationships within food web dynamics involving reptilian predators. Our findings, therefore, provide limited field support for the hypothesis that stage structure of predator populations with ODS is an important consideration in accurately assessing terrestrial trophic linkages. Of note, however, are quantitative factors limiting the interpretation of these data, including how we measured prey response, the lack of spatial replication, and temporal dependence between the time series data.

Based on the scale of the study and the resources required to effectively study cryptic species at the intensity required, replication of our study was not possible. These data thus represent a single event, in a single location confounded by environmental or temporal factors for which we did not control. Temporal autocorrelation is documented to have substantive effects on populations when modeling stochastic population growth47. Although the effects of temporal autocorrelation should be lower in our study, they are likely not entirely negligible based on the slightly different statistical relationships between the snake density and prey SPUE within the time series and negative binomial regressions. Thus, although quantitative limitations in interpretation should be considered as limiting the strength of support, the data add to broader theoretical ideas on ODS as a component of trophic relationships.

Additionally, because the study site described was closed only to snakes, rats and lizards could immigrate or emigrate freely. Thus, prey SPUE may have been a function of movement, activity levels, and reproduction. Linkages between predator and prey may be tighter in closed populations where recruitment or mortality is the only mechanism for population changes for both predator and prey. Likewise, our results also rely on the assumption that SPUE is a relatively accurate index of true abundance. In general, capture per unit effort (CPUE), as a measure of abundance, is coarse and may be problematic without prior species-specific validation48. There are many known and unknown sources of error in using CPUE to index population abundance for lizards in Guam forests33,34. The relationship between nocturnal visual sightings of rodents and their abundance in tropical forests is untested. Known concerns with visual mammal survey effectiveness include having well-worn paths and rested surveyors, which should not be hindrances in our study49. Regardless, visual detections may be driven by both direct effects (true abundance) and indirect effects, such as reduced activity or area avoidance or attraction, which hinder our ability to determine whether decreases in SPUE reflected fewer individuals or sightings of individuals that were present.

The mechanistic explanations for changes in relative prey SPUE in response to changing predator populations may, however, be less relevant when a biological response (i.e., changes in activity or abundance) was still observed. Specifically, the relationship between rat SPUE and the production of juvenile snakes after a suitable lag period supports the idea that the mechanism that caused changes in rat SPUE also affected their availability as prey for snakes. As measured by mass, rodents are a larger prey reward than most lizards that occur in Guam. Reproduction tends to metabolically costly, with a 30% increase in metabolic activity to support egg yolk deposition in viviparous snakes and measured rates of energy assimilation from meals at upwards of 90% in pythons50,51. When each brown treesnake egg is roughly 8 g52, a clutch of 3 or more eggs would require consuming quite a few < 5 g lizards to support normal metabolic costs as well as reproduction, while one 300 g rat may support both functions easily. Lizards may, therefore, primarily be intermediate meals that support survival, while larger prey items are needed to reproduce, but additional data would be necessary to corroborate this hypothesis. Our data support that a measurable increase in rat SPUE potentially predicted a reproductive pulse for BTS, but that changes in lizard SPUE did not.

Although brown treesnakes are reliably interested in endothermic prey around 900 mm SVL, lizards remain a major part of their diet into sizes well over 1100 mm SVL17,28,29. The maximum size of a snake in this population was 1299 mm SVL, indicating almost all individuals in the studied population likely consumed lizards. The long-term persistence of smaller prey items in the diet of species with ontogenetic prey changes is documented in several fish and reptile species53–56. In the forests of Guam, very few endothermic prey remain, with the few remnant native birds mostly isolated to small areas in urban landscapes23,57. Thus, daily energetic and metabolic needs may require that snakes maintain lizards in their diet. Snakes collected from forested areas on Guam, where the primary prey are anticipated to be various lizard species, maintained lizards as a higher proportion of their diet into larger size classes than snakes in urban landscapes17. Therefore, lizard populations are likely also affected by larger size classes and prior work has demonstrated lizard counts increase with brown treesnake eradication58. Trophic linkages between larger snakes and lizards may, however, be weaker as compared to those between juvenile snakes that are dependent on lizards as their primary food source. Thus, the ability for brown treesnakes to exert top-down control on ectothermic prey species is reasonable and documented. The importance, however, of considering the ontogeny in populations brown treesnake or other reptile populations in exerting top-down effects has not, to our knowledge, been previously reported.

Although the total population was reduced by ~ 40% after treatment, we did not detect a direct relationship between snake density and lizard response based on the total population density of snakes. The stronger response of rodents to the reduction of snakes, as compared to lizards, may have been biased by the use of rodents as the snake removal tool. This approach should favor the removal of snakes foraging on rodents. The separation of the population into different stages was, however, effective in unmasking a relationship between snakes and lizards that was otherwise not detectable based on the methods we used for removal. Thus, failing to consider ODS in trophic interactions and population stage structure may miss trophic linkages in community food web patterns13,59,60. For herpetofauna, this may contribute bias in understanding their overall ecological impacts or relative importance, as diet shifts in prey are relatively common in predatory herpetofauna.

Management implications

The stage-specific and ontogenetic effects of predators on prey may be especially important to consider for invasive snake management. Native prey species may have altered risk of depredation depending on the stage structure of the invasive predator population. Regarding brown treesnakes, the primary removal tools are based on food attractants that target a subset of the population29,61,62. This removal approach may have implications for broader management goals, if conservation or recovery goals are targeting different species than those used for control42. In general, reptiles have important roles in communities and injurious invasive species, with several examples of major perturbations following the introduction of snakes21,63. ODS are common among reptiles and may also co-occur with other types of niche shifts64. In addition to many-fold ontogenetic changes in size, which can affect decisions such as toxicant dosing, individuals may display stage-based preferences for prey states36,61,65. In systems in which an invasive predator affects many trophic levels and ecological recovery is desired, consideration of population structure may be meaningful for both planning species recovery and to achieve removal targets.

Acknowledgements

We would like to thank L. Bonewell, A. Collins, and R.N. Reed for administrative and logistical assistance. We would also like to thank B. Lardner, J. Savidge, S. Goetz, A. Yackel Adams, and R.N. Reed for topically related intellectual discourse. We thank S. Goetz and five anonymous reviewers for providing comments on an earlier draft of this manuscript. This work was supported by funding provided by the Department of the Navy, Joint Region Marianas (IAA #N61128-20140529-0001) and Department of Interior, Office of Insular Affairs. We thank Andersen Air Force Base for access to the enclosure and M. Mazurek and H. McCoy for coordination of financial support. We thank the many field technicians that participated in the collection of these data. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Author contributions

MGN conceived and designed experiments and statistical analysis. MGN and LNG wrote the manuscript. Both authors contributed to revisions.

Data availability

The datasets generated and/or analyzed during the current study are available for download (Nafus et al.35, 10.5066/P9QRWKQB, Nafus43, 10.5066/P9JQ9HG0).

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Nancy AK Wise DH Prey availability limits fecundity and influences the movement pattern of female fishing spiders Oecologia 2001 127 417 424 10.1007/s004420000607 28547112
Nancy, A. K. & Wise, D. H. Prey availability limits fecundity and influences the movement pattern of female fishing spiders. Oecologia 127, 417–424 (2001).28547112 10.1007/s004420000607
2. Oro D Influence of density dependence on predator-prey seabird interactions at large spatio-temporal scales Proc. R. Soc. B: Biol. Sci. 2006 273 379 383 10.1098/rspb.2005.3287
Oro, D. et al. Influence of density dependence on predator-prey seabird interactions at large spatio-temporal scales. Proc. R. Soc. B: Biol. Sci. 273, 379–383 (2006).10.1098/rspb.2005.3287
3. Wangersky PJ Lotka-volterra population models Annu. Rev. Ecol. Evol. Syst. 1978 1978 9
Wangersky, P. J. Lotka-volterra population models. Annu. Rev. Ecol. Evol. Syst. 1978, 9 (1978).
4. McCann K Hastings A Huxel GR Weak trophic interactions and the balance of nature Nature. 1998 395 794 798 10.1038/27427
McCann, K., Hastings, A. & Huxel, G. R. Weak trophic interactions and the balance of nature. Nature. 395, 794–798 (1998).10.1038/27427
5. de Roos AM Dynamic population stage structure due to juvenile–adult asymmetry stabilizes complex ecological communities Proc. Natl. Acad. Sci. USA 2021 118 e2023709118 10.1073/pnas.2023709118 34021084
de Roos, A. M. Dynamic population stage structure due to juvenile–adult asymmetry stabilizes complex ecological communities. Proc. Natl. Acad. Sci. USA 118, e2023709118 (2021).34021084 10.1073/pnas.2023709118
6. de Roos AM Persson L Competition in size-structured populations: Mechanisms inducing cohort formation and population cycles Theor. Popul. Biol. 2003 63 1 16 10.1016/S0040-5809(02)00009-6 12464491
de Roos, A. M. & Persson, L. Competition in size-structured populations: Mechanisms inducing cohort formation and population cycles. Theor. Popul. Biol. 63, 1–16 (2003).12464491 10.1016/S0040-5809(02)00009-6
7. Nilsson KA McCann KS Caskenette AL Interaction strength and stability in stage-structured food web modules Oikos. 2018 127 1494 1505 10.1111/oik.05029
Nilsson, K. A., McCann, K. S. & Caskenette, A. L. Interaction strength and stability in stage-structured food web modules. Oikos. 127, 1494–1505 (2018).10.1111/oik.05029
8. van Leeuwen A Predators with multiple ontogenetic niche shifts have limited potential for population growth and top-down control of their prey Am. Nat. 2013 182 53 66 10.1086/670614 23778226
van Leeuwen, A. et al. Predators with multiple ontogenetic niche shifts have limited potential for population growth and top-down control of their prey. Am. Nat. 182, 53–66 (2013).23778226 10.1086/670614
9. van Leeuwen A Huss M Gårdmark A de Roos AM Ontogenetic specialism in predators with multiple niche shifts prevents predator population recovery and establishment Ecology. 2014 95 2409 2422 10.1890/13-0843.1
van Leeuwen, A., Huss, M., Gårdmark, A. & de Roos, A. M. Ontogenetic specialism in predators with multiple niche shifts prevents predator population recovery and establishment. Ecology. 95, 2409–2422 (2014).10.1890/13-0843.1
10. De Roos AM Persson L McCauley E The influence of size-dependent life-history traits on the structure and dynamics of populations and communities Ecol. Lett. 2003 6 473 487 10.1046/j.1461-0248.2003.00458.x
De Roos, A. M., Persson, L. & McCauley, E. The influence of size-dependent life-history traits on the structure and dynamics of populations and communities. Ecol. Lett. 6, 473–487 (2003).10.1046/j.1461-0248.2003.00458.x
11. Lawson CL The influence of ontogenetic diet variation on consumption rate estimates: A marine example Sci. Rep. 2018 8 10725 10.1038/s41598-018-28479-7 30013084
Lawson, C. L. et al. The influence of ontogenetic diet variation on consumption rate estimates: A marine example. Sci. Rep. 8, 10725 (2018).30013084 10.1038/s41598-018-28479-7
12. Rudolf HW Rasmussen RL Ontogenetic functional diversity: Size structure of a keystone predator drives functioning of a complex ecosystem Ecology 2013 94 1047 1056 10.1890/12-0378.1
Rudolf, H. W. & Rasmussen, R. L. Ontogenetic functional diversity: Size structure of a keystone predator drives functioning of a complex ecosystem. Ecology 94, 1047–1056 (2013).10.1890/12-0378.1
13. Sánchez-Hernández J Nunn AD Adams CE Amundsen P-A Causes and consequences of ontogenetic dietary shifts: A global synthesis using fish models Biol. Rev. 2019 94 539 554 10.1111/brv.12468 30251433
Sánchez-Hernández, J., Nunn, A. D., Adams, C. E. & Amundsen, P.-A. Causes and consequences of ontogenetic dietary shifts: A global synthesis using fish models. Biol. Rev. 94, 539–554 (2019).30251433 10.1111/brv.12468
14. Simonis JL Predator ontogeny determines trophic cascade strength in freshwater rock pools Ecosphere. 2013 4 62 10.1890/ES13-00019.1
Simonis, J. L. Predator ontogeny determines trophic cascade strength in freshwater rock pools. Ecosphere. 4, 62. 10.1890/ES13-00019.1 (2013).10.1890/ES13-00019.1
15. Bouchard SS Bjorndal KA Ontogenetic diet shifts and digestive constraints in the omnivorous freshwater turtle Trachemys scripta Physiol. Biochem. Zool. 2006 79 150 158 10.1086/498190 16380936
Bouchard, S. S. & Bjorndal, K. A. Ontogenetic diet shifts and digestive constraints in the omnivorous freshwater turtle Trachemys scripta. Physiol. Biochem. Zool. 79, 150–158 (2006).16380936 10.1086/498190
16. Mautz WJ Nagy KA Ontogenetic changes in diet, field metabolic rate, and water flux in the herbivorous lizard Dipsosaurus dorsalis Physiol. Zool. 1987 60 640 658 10.1086/physzool.60.6.30159980
Mautz, W. J. & Nagy, K. A. Ontogenetic changes in diet, field metabolic rate, and water flux in the herbivorous lizard Dipsosaurus dorsalis. Physiol. Zool. 60, 640–658 (1987).10.1086/physzool.60.6.30159980
17. Savidge JA Food habits of Boiga irregularis, an introduced predator on Guam J. Herpetol. 1988 22 275 282 10.2307/1564150
Savidge, J. A. Food habits of Boiga irregularis, an introduced predator on Guam. J. Herpetol. 22, 275–282 (1988).10.2307/1564150
18. Palci A Lee MSY Hutchinson MN Patterns of postnatal ontogeny of the skull and lower jaw of snakes as revealed by micro-CT scan data and three-dimensional geometric morphometrics J. Anat. 2016 229 723 754 10.1111/joa.12509 27329823
Palci, A., Lee, M. S. Y. & Hutchinson, M. N. Patterns of postnatal ontogeny of the skull and lower jaw of snakes as revealed by micro-CT scan data and three-dimensional geometric morphometrics. J. Anat. 229, 723–754 (2016).27329823 10.1111/joa.12509
19. de Miranda EBP The plight of reptiles as ecological actors in the tropics Front. Ecol. Evol. 2017 10.3389/fevo.2017.00159
de Miranda, E. B. P. The plight of reptiles as ecological actors in the tropics. Front. Ecol. Evol.10.3389/fevo.2017.00159 (2017).10.3389/fevo.2017.00159
20. Rogers HS Effects of an invasive predator cascade to plants via mutualism disruption Nat. Commun. 2017 8 14557 10.1038/ncomms14557 28270682
Rogers, H. S. et al. Effects of an invasive predator cascade to plants via mutualism disruption. Nat. Commun. 8, 14557 (2017).28270682 10.1038/ncomms14557
21. Dorcas ME Severe mammal declines coincide with proliferation of invasive Burmese pythons in Everglades National Park Proc. Natl. Acad. Sci. USA 2012 109 2418 2422 10.1073/pnas.1115226109 22308381
Dorcas, M. E. et al. Severe mammal declines coincide with proliferation of invasive Burmese pythons in Everglades National Park. Proc. Natl. Acad. Sci. USA 109, 2418–2422 (2012).22308381 10.1073/pnas.1115226109
22. Schoener, T. W. & Spiller, D. A. Trophic Cascades on Islands in Trophic Cascades: Predators, Prey, and the Changing Dynamics of Nature, (eds Terborgh, J. & Estes, J.) 179–202 (University of Chicago Press, 2010).
23. Wiles GJ Bart J Beck RE Aguon CF Impacts of the brown tree snake: Patterns of decline and species persistence in Guam's avifauna Conserv. Biol. 2003 17 1350 1360 10.1046/j.1523-1739.2003.01526.x
Wiles, G. J., Bart, J., Beck, R. E. & Aguon, C. F. Impacts of the brown tree snake: Patterns of decline and species persistence in Guam’s avifauna. Conserv. Biol. 17, 1350–1360 (2003).10.1046/j.1523-1739.2003.01526.x
24. Wiewel AS YackelAdams AA Rodda GH Distribution, density, and biomass of introduced small mammals in the southern Mariana Islands Pac. Sci. 2009 63 205 222 10.2984/049.063.0204
Wiewel, A. S., YackelAdams, A. A. & Rodda, G. H. Distribution, density, and biomass of introduced small mammals in the southern Mariana Islands. Pac. Sci. 63, 205–222 (2009).10.2984/049.063.0204
25. Rodda GH Savidge JA Biology and impacts of Pacific Island invasive species. 2. Boiga irregularis, the brown tree snake (Reptilia: Colubridae) Pac. Sci. 2007 61 307 324 10.2984/1534-6188(2007)61[307:BAIOPI]2.0.CO;2
Rodda, G. H. & Savidge, J. A. Biology and impacts of Pacific Island invasive species. 2. Boiga irregularis, the brown tree snake (Reptilia: Colubridae). Pac. Sci. 61, 307–324 (2007).10.2984/1534-6188(2007)61[307:BAIOPI]2.0.CO;2
26. McCleery RA Marsh rabbit mortalities tie pythons to the precipitous decline of mammals in the Everglades Proc. R. Soc. B: Biol. Sci. 2015 282 e20120120 10.1098/rspb.2015.0120
McCleery, R. A. et al. Marsh rabbit mortalities tie pythons to the precipitous decline of mammals in the Everglades. Proc. R. Soc. B: Biol. Sci. 282, e20120120 (2015).10.1098/rspb.2015.0120
27. Savidge JA Extinction of an island forest avifauna by an introduced snake Ecology. 1987 68 660 668 10.2307/1938471
Savidge, J. A. Extinction of an island forest avifauna by an introduced snake. Ecology. 68, 660–668 (1987).10.2307/1938471
28. Rodda GH Dean-Bradley K Size selectivity of brown treesnake traps Micronesica. 2004 37 180 181
Rodda, G. H. & Dean-Bradley, K. Size selectivity of brown treesnake traps. Micronesica. 37, 180–181 (2004).
29. Rodda GH Size bias in visual searches and trapping of brown treesnakes on Guam J. Wildl. Manag. 2007 71 656 661 10.2193/2005-742
Rodda, G. H. et al. Size bias in visual searches and trapping of brown treesnakes on Guam. J. Wildl. Manag. 71, 656–661 (2007).10.2193/2005-742
30. Nafus MG Foraging behavior in a generalist snake (brown treesnake, Boiga irregularis) with implications for avian reintroduction and recovery Appl. Anim. Behav. Sci. 2021 243 105450 10.1016/j.applanim.2021.105450
Nafus, M. G. et al. Foraging behavior in a generalist snake (brown treesnake, Boiga irregularis) with implications for avian reintroduction and recovery. Appl. Anim. Behav. Sci. 243, 105450 (2021).10.1016/j.applanim.2021.105450
31. Wolda H Seasonal fluctuations in rainfall, food and abundance of tropical insects J. Anim. Ecol. 1978 47 369 381 10.2307/3789
Wolda, H. Seasonal fluctuations in rainfall, food and abundance of tropical insects. J. Anim. Ecol. 47, 369–381 (1978).10.2307/3789
32. Nafus MG Yackel Adams AA Klug PE Rodda GH Habitat type and structure affect trap capture success of an invasive snake across variable densities Ecosphere. 2018 9 e02339 10.1002/ecs2.2339
Nafus, M. G., Yackel Adams, A. A., Klug, P. E. & Rodda, G. H. Habitat type and structure affect trap capture success of an invasive snake across variable densities. Ecosphere. 9, e02339 (2018).10.1002/ecs2.2339
33. Lardner B Detection rates of geckos in visual surveys: Turning confounding variables into useful knowledge J. Herpetol. 2015 49 522 532 10.1670/14-048
Lardner, B. et al. Detection rates of geckos in visual surveys: Turning confounding variables into useful knowledge. J. Herpetol. 49, 522–532 (2015).10.1670/14-048
34. Rodda GH Stability of detectability over 17 years at a single site and other lizard detection comparisons from Guam J. Herpetol. 2015 49 513 521 10.1670/14-085
Rodda, G. H. et al. Stability of detectability over 17 years at a single site and other lizard detection comparisons from Guam. J. Herpetol. 49, 513–521 (2015).10.1670/14-085
35. Nafus, M. G. et al. Guam, USGS Closed Population (NWFN), an experimental eradication of brown treesnakes in a 5-ha study site, 2016–2023, U.S. Geological Survey data release. 10.5066/P9QRWKQB (2023).
36. Lardner B Savidge JA Rodda GH Reed RN Prey preferences and prey acceptance in juvenile brown treesnakes (Boiga irregularis) Herpetol. Conserv. Biol. 2009 4 313 323
Lardner, B., Savidge, J. A., Rodda, G. H. & Reed, R. N. Prey preferences and prey acceptance in juvenile brown treesnakes (Boiga irregularis). Herpetol. Conserv. Biol. 4, 313–323 (2009).
37. Johnston JJ Risk assessment of an acetaminophen baiting program for chemical control of brown tree snakes on Guam: Evaluation of baits, snake residues, and potential primary and secondary hazards Environ. Sci. Technol. 2002 36 3827 3833 10.1021/es015873n 12322757
Johnston, J. J. et al. Risk assessment of an acetaminophen baiting program for chemical control of brown tree snakes on Guam: Evaluation of baits, snake residues, and potential primary and secondary hazards. Environ. Sci. Technol. 36, 3827–3833 (2002).12322757 10.1021/es015873n
38. Nafus MG Demographic response of brown treesnakes to extended population suppression J. Wildl. Manag. 2022 86 e22136 10.1002/jwmg.22136
Nafus, M. G. et al. Demographic response of brown treesnakes to extended population suppression. J. Wildl. Manag. 86, e22136 (2022).10.1002/jwmg.22136
39. Mauldin RE Savarie PJ Acetaminophen as an oral toxicant for Nile monitor lizards (Varanus niloticus) and Burmese pythons (Python molurus bivittatus) Wildl. Res. 2010 37 215 222 10.1071/WR08168
Mauldin, R. E. & Savarie, P. J. Acetaminophen as an oral toxicant for Nile monitor lizards (Varanus niloticus) and Burmese pythons (Python molurus bivittatus). Wildl. Res. 37, 215–222 (2010).10.1071/WR08168
40. Christy MT Modeling detection probabilities to evaluate management and control tools for an invasive species J. Appl. Ecol. 2010 47 106 113 10.1111/j.1365-2664.2009.01753.x
Christy, M. T. et al. Modeling detection probabilities to evaluate management and control tools for an invasive species. J. Appl. Ecol. 47, 106–113 (2010).10.1111/j.1365-2664.2009.01753.x
41. Tyrrell CL Evaluation of trap capture in a geographically closed population of brown treesnakes on Guam J. Appl. Ecol. 2009 46 128 135 10.1111/j.1365-2664.2008.01591.x
Tyrrell, C. L. et al. Evaluation of trap capture in a geographically closed population of brown treesnakes on Guam. J. Appl. Ecol. 46, 128–135 (2009).10.1111/j.1365-2664.2008.01591.x
42. Nafus MG Reyes A Fies T Goetz SM Adaptive resource management (ARM): Achieving functional eradication of invasive snakes to benefit avian conservation J. Appl. Ecol. 2024 61 733 745 10.1111/1365-2664.14597
Nafus, M. G., Reyes, A., Fies, T. & Goetz, S. M. Adaptive resource management (ARM): Achieving functional eradication of invasive snakes to benefit avian conservation. J. Appl. Ecol. 61, 733–745 (2024).10.1111/1365-2664.14597
43. Nafus, M. G. Guam, USGS Closed Population (NWFN) data relating to brown treesnake and prey interactions processed into monthly intervals from 10/2016–2/2023, U.S. Geological Survey data release. 10.5066/P9JQ9HG0 (2023).
44. Lardner B Do observer fatigue and taxon bias compromise visual encounter surveys for small vertebrates? Wildl. Res. 2019 46 127 135 10.1071/WR18016
Lardner, B. et al. Do observer fatigue and taxon bias compromise visual encounter surveys for small vertebrates?. Wildl. Res. 46, 127–135 (2019).10.1071/WR18016
45. Venables, W. N. & Ripley, B. D. Generalized linear models. In Modern Applied Statistics, 183–210 (Springer, 2002).
46. Hyndman, R. et al. Forecast: Forecasting functions for time series and linear models. R package version 8.21.1 (2023).
47. Tuljapurkar S Haridas CV Temporal autocorrelation and stochastic population growth Ecol. Lett. 2006 9 327 337 10.1111/j.1461-0248.2006.00881.x 16958899
Tuljapurkar, S. & Haridas, C. V. Temporal autocorrelation and stochastic population growth. Ecol. Lett. 9, 327–337 (2006).16958899 10.1111/j.1461-0248.2006.00881.x
48. Harley SJ Myers RA Dunn A Is catch-per-unit-effort proportional to abundance? Can. J. Fish Aquat. Sci. 2001 58 1760 1772 10.1139/f01-112
Harley, S. J., Myers, R. A. & Dunn, A. Is catch-per-unit-effort proportional to abundance?. Can. J. Fish Aquat. Sci. 58, 1760–1772 (2001).10.1139/f01-112
49. Munari DP Keller C Venticinque EM An evaluation of field techniques for monitoring terrestrial mammal populations in Amazonia Mamm. Biol. 2011 76 401 408 10.1016/j.mambio.2011.02.007
Munari, D. P., Keller, C. & Venticinque, E. M. An evaluation of field techniques for monitoring terrestrial mammal populations in Amazonia. Mamm. Biol. 76, 401–408 (2011).10.1016/j.mambio.2011.02.007
50. Van Dyke JU Beaupre SJ Bioenergetic components of reproductive effort in viviparous snakes: Costs of vitellogenesis exceed costs of pregnancy Comp. Biochem. Physiol. A Mol. Integr. Physiol. 2011 160 504 515 10.1016/j.cbpa.2011.08.011 21884815
Van Dyke, J. U. & Beaupre, S. J. Bioenergetic components of reproductive effort in viviparous snakes: Costs of vitellogenesis exceed costs of pregnancy. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 160, 504–515 (2011).21884815 10.1016/j.cbpa.2011.08.011
51. Cox CL Secor SM Effects of meal size, clutch, and metabolism on the energy efficiencies of juvenile Burmese pythons, Python molurus Comp. Biochem. Physiol. A Mol. Integr. Physiol. 2007 148 861 868 10.1016/j.cbpa.2007.08.029 17913527
Cox, C. L. & Secor, S. M. Effects of meal size, clutch, and metabolism on the energy efficiencies of juvenile Burmese pythons, Python molurus. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 148, 861–868 (2007).17913527 10.1016/j.cbpa.2007.08.029
52. Mathies T Miller LA Cool temperatures elicit reproduction in a biologically invasive predator, the brown treesnake (Boiga irregularis) Zoo. Biol. 2003 22 227 238 10.1002/zoo.10084
Mathies, T. & Miller, L. A. Cool temperatures elicit reproduction in a biologically invasive predator, the brown treesnake (Boiga irregularis). Zoo. Biol. 22, 227–238 (2003).10.1002/zoo.10084
53. Schoener TW Gorman GC Some niche differences in three Lesser Antillean lizards of the genus Anolis Ecology. 1968 49 819 830 10.2307/1936533
Schoener, T. W. & Gorman, G. C. Some niche differences in three Lesser Antillean lizards of the genus Anolis. Ecology. 49, 819–830 (1968).10.2307/1936533
54. Eggold BT Motta PJ Ontogenetic dietary shifts and morphological correlates in striped mullet, Mugil cephalus Environ. Biol. Fishes. 1992 34 139 158 10.1007/BF00002390
Eggold, B. T. & Motta, P. J. Ontogenetic dietary shifts and morphological correlates in striped mullet, Mugil cephalus. Environ. Biol. Fishes. 34, 139–158 (1992).10.1007/BF00002390
55. Muñoz AA Ojeda FP Guild structure of carnivorous intertidal fishes of the Chilean coast: Implications of ontogenetic dietary shifts Oecologia. 1998 114 563 573 10.1007/s004420050481 28307906
Muñoz, A. A. & Ojeda, F. P. Guild structure of carnivorous intertidal fishes of the Chilean coast: Implications of ontogenetic dietary shifts. Oecologia. 114, 563–573 (1998).28307906 10.1007/s004420050481
56. Lahti ME Beck DD Ecology and ontogenetic variation of diet in the pigmy short-horned lizard (Phrynosoma douglasii) Am. Midl. Nat. 2008 159 327 339 10.1674/0003-0031(2008)159[327:EAOVOD]2.0.CO;2
Lahti, M. E. & Beck, D. D. Ecology and ontogenetic variation of diet in the pigmy short-horned lizard (Phrynosoma douglasii). Am. Midl. Nat. 159, 327–339 (2008).10.1674/0003-0031(2008)159[327:EAOVOD]2.0.CO;2
57. Pollock HS Pervasive impacts of invasive brown treesnakes drive low fledgling survival in endangered Micronesian Starlings (Aplonis opaca) on Guam Condor. 2019 121 1 11 10.1093/condor/duz014
Pollock, H. S. et al. Pervasive impacts of invasive brown treesnakes drive low fledgling survival in endangered Micronesian Starlings (Aplonis opaca) on Guam. Condor. 121, 1–11 (2019).10.1093/condor/duz014
58. Campbell EW Do predators control prey species abundance? An experimental test with brown treesnakes on Guam Ecology. 2012 93 1194 1203 10.1890/11-1359.1 22764505
Campbell, E. W. et al. Do predators control prey species abundance? An experimental test with brown treesnakes on Guam. Ecology. 93, 1194–1203 (2012).22764505 10.1890/11-1359.1
59. Harvey CJ Kitchell JF A stable isotope evaluation of the structure and spatial heterogeneity of a Lake Superior food web Can. J. Fish Aquat. Sci. 2000 57 1395 1403 10.1139/f00-072
Harvey, C. J. & Kitchell, J. F. A stable isotope evaluation of the structure and spatial heterogeneity of a Lake Superior food web. Can. J. Fish Aquat. Sci. 57, 1395–1403 (2000).10.1139/f00-072
60. Cross WF Food-web dynamics in a large river discontinuum Ecol. Monogr. 2013 83 311 337 10.1890/12-1727.1
Cross, W. F. et al. Food-web dynamics in a large river discontinuum. Ecol. Monogr. 83, 311–337 (2013).10.1890/12-1727.1
61. Shivik JA Clark L Ontogenetic shifts in carion attractiveness to brown tree snakes (Boiga irregularis) J. Herpetol. 1999 33 334 336 10.2307/1565737
Shivik, J. A. & Clark, L. Ontogenetic shifts in carion attractiveness to brown tree snakes (Boiga irregularis). J. Herpetol. 33, 334–336 (1999).10.2307/1565737
62. Clark L Clark C Siers SR Pitt WC Beasley JC Witmer GW Brown tree snakes methods and approaches for control Ecology and Management of Terrestrial Vertebrate Invasive Species in the United States 2018 CRC Press 107 134
Clark, L., Clark, C. & Siers, S. R. Brown tree snakes methods and approaches for control. In Ecology and Management of Terrestrial Vertebrate Invasive Species in the United States (eds Pitt, W. C. et al.) 107–134 (CRC Press, 2018).
63. Piquet JC López-Darias M Invasive snake causes massive reduction of all endemic herpetofauna on Gran Canaria Proc. R. Soc. B: Biol. Sci. 2021 288 20211939 10.1098/rspb.2021.1939
Piquet, J. C. & López-Darias, M. Invasive snake causes massive reduction of all endemic herpetofauna on Gran Canaria. Proc. R. Soc. B: Biol. Sci. 288, 20211939 (2021).10.1098/rspb.2021.1939
64. Whitfield SM Donnelly MA Ontogenetic and seasonal variation in the diets of a Costa Rican leaf-litter herpetofauna J. Trop. Ecol. 2006 22 409 417 10.1017/S0266467406003245
Whitfield, S. M. & Donnelly, M. A. Ontogenetic and seasonal variation in the diets of a Costa Rican leaf-litter herpetofauna. J. Trop. Ecol. 22, 409–417 (2006).10.1017/S0266467406003245
65. Siers SR Goetz SM Volsteadt RM Nafus MG Evaluating lethal toxicant doses for the largest individuals of an invasive vertebrate predator with indeterminate growth Manag. Biol. Invasions. 2021 12 476 494 10.3391/mbi.2021.12.2.17
Siers, S. R., Goetz, S. M., Volsteadt, R. M. & Nafus, M. G. Evaluating lethal toxicant doses for the largest individuals of an invasive vertebrate predator with indeterminate growth. Manag. Biol. Invasions. 12, 476–494 (2021).10.3391/mbi.2021.12.2.17
