The impact of invasive cane toads on native wildlife in southern Australia

Abstract Commonly, invaders have different impacts in different places. The spread of cane toads (Rhinella marina: Bufonidae) has been devastating for native fauna in tropical Australia, but the toads' impact remains unstudied in temperate‐zone Australia. We surveyed habitat characteristics and fauna in campgrounds along the central eastern coast of Australia, in eight sites that have been colonized by cane toads and another eight that have not. The presence of cane toads was associated with lower faunal abundance and species richness, and a difference in species composition. Populations of three species of large lizards (land mullets Bellatorias major, eastern water dragons Intellagama lesueurii, and lace monitors Varanus varius) and a snake (red‐bellied blacksnake Pseudechis porphyriacus) were lower (by 84 to 100%) in areas with toads. The scarcity of scavenging lace monitors in toad‐invaded areas translated into a 52% decrease in rates of carrion removal (based on camera traps at bait stations) and an increase (by 61%) in numbers of brush turkeys (Alectura lathami). The invasion of cane toads through temperate‐zone Australia appears to have reduced populations of at least four anurophagous predators, facilitated other taxa, and decreased rates of scavenging. Our data identify a paradox: The impacts of cane toads are at least as devastating in southern Australia as in the tropics, yet we know far more about toad invasion in the sparsely populated wilderness areas of tropical Australia than in the densely populated southeastern seaboard.


Introduction
Invasive species imperil native biodiversity (Mack et al. 2000;McGeoch et al. 2010), but invader impacts are highly heterogeneous (Melbourne et al. 2007). Some invaders have catastrophic impacts, whereas others may benefit native taxa; some native taxa are more vulnerable than others (Wonham et al. 2005;King et al. 2006;Brown et al. 2011;Simberloff 2011). Invader impact can vary even among populations of the same native species (e.g., Letnic et al. 2008;Somaweera and Shine 2012). The details of most invader impacts remain unclear because of complex direct and indirect pathways of effect (Shine 2012); a decline in one species may have positive or negative flow-on effects to many others (Doody et al. 2006;Brown et al. 2011Brown et al. , 2013aEstes et al. 2011).
Impact heterogeneity necessitates studies at multiple locations throughout an invader's range (Parker et al. 1999;Melbourne et al. 2007). To assess an invader's effects, we must measure changes in abundance of native species coincident with the presence of an invader. This task may be simple in some cases (e.g., for sessile organisms), but is more challenging if the impact falls on vagile, rare predators (Caughley 1977;Woinarski et al. 2001Woinarski et al. , 2004. To incorporate the influence of habitat and climate, studies of invader impacts on large, mobile vertebrates must employ multiple sites with a sampling design that captures a range of natural heterogeneity. Since its introduction in 1935, the cane toad (Rhinella marina) has spread rapidly through Australia (Lever 2001;Kolbe, Kearney & Shine 2010). Extensive research in tropical Australia has demonstrated that the arrival of cane toads is consistently followed by population-level declines of some species of large predators that are fatally poisoned when they eat the toxic toads (Letnic et al. 2008;Doody et al. 2009). However, impacts on the (virtually unstudied) southern edge of the toads' range expansion may be very different from those in the tropics. For example, the tropical invasion front is dominated by large adult toads; any predator that consumes a large toad will be fatally poisoned (Shine 2010). In contrast, the southern front contains small as well as large toads (McCann 2014). A small toad offers a nauseating but nonfatal meal that may allow aversion learning by predators, thereby ameliorating population-level impact (O'Donnell, Webb & Shine 2010).
The temperate-zone invasion front also facilitates spatial comparisons of affected versus unaffected sites. In the tropics, toads invade in a rapidly advancing, continuous front hundreds of kilometers wide, precluding comparisons between adjacent sites with versus without toads (Shine 2010). Thus, these studies yield before-after data on population-level impacts of toad invasion, a critical difficulty in the wet-dry tropics where even minor yearto-year variation in rainfall patterns can induce massive changes in faunal populations (Brown et al. 2002;Madsen et al. 2006;Shine and Brown 2008). Those temporal shifts make it difficult to attribute specific faunal shifts to cane toad invasion. For example, some declines in predator populations coincident with toad invasion in tropical Australia were caused by stochastic weather events, not toads (Brown et al. 2011).
The toad invasion into southern Australia is progressing slowly  and in a patchwork fashion (likely due to multiple translocations of toads by humans: White and Shine 2009). Adjacent areas with versus without cane toads create an opportunity for direct and concurrent comparisons. Such spatial comparisons do not overcome all of the problems associated with inferring toad impact-for example, invaded and uninvaded sites may differ in ways that affect faunal assemblages (Suarez et al. 1998;Bolger et al. 2000;Holway et al. 2002a,b). Nonetheless, confounding variables can be measured; and such a study is not weakened by the influence of annual variation.
Despite logistical advantages (a fragmented invasion front, a significantly higher human population, and infrastructure that facilitates research), the impact of cane toads in southern Australia has been largely ignored. We conducted surveys to quantify the characteristics of native faunal assemblages in adjacent sites that contained versus those that did not contain cane toads to test what impact toads have on fauna assemblages. We predicted that native fauna that prey on toads would be less abundant in toad-invaded areas than in areas that did not contain toads. We also surveyed habitat variables to test for potential confounding factors between invaded versus uninvaded sites, and measured rates of carrion removal to test the hypothesis that this critical ecological function would be affected by toad-induced declines of scavengers.

Study region
The Northern Rivers region encompasses the northeastern corner (50,266 km 2 ) of New South Wales (NSW), Australia (DECCW 2010; Fig. 1). Although it occupies less than 10% of the state, this area supports more than 40% of the threatened species and 20% of the threatened communities of NSW (Goldingay et al. 1999;(DECCW 2010;Newell 2011). The region experiences a warm temperate to subtropical climate with warm humid summers (mean average maximum 26-30°C) and moderate winters (mean average minimum 6-10°C: Bureau of Meteorology 2014). Most rainfall occurs during summer and early autumn (Taylor and Goldingay 2003). Many native taxa (both ectotherms and endotherms) of this region are most active in the warmer, wetter months (September to February : Shine 1979;York et al. 1991;Christian and Weavers 1996;Kavanagh and Stanton 2005;Daly and Lemckert 2011), so we confined our surveys to the spring and summer months of 2013-2014.

Study species
Cane toads (Rhinella marina) are large (to 230-mm snout-urostyle length [SUL] and 500 g mass: Zug and Zug 1979;Brown et al. 2013a,b) toxic bufonid anurans (Tyler 1975;Pramuk 2006), native to Central and South America (Lever 2001). Introduced from Hawaii to northeastern QLD in 1935, this highly toxic anuran has since invaded more than 1.2 million km 2 of tropical and subtropical Australia (Lever 2001;Urban et al. 2007Urban et al. , 2008 and has also spread into southeastern QLD and northeastern NSW (Seabrook 1993;Lever 2003). Cane toads were introduced to NSW between 1964 and 1966, when a satellite population was established in Byron Bay (van Beurden and Grigg 1980). By 1989, this satellite population had merged to form a continuous population from the Northern Rivers region of NSW into QLD (Seabrook 1991(Seabrook , 1993. Presumably due to constraints of climatic conditions on toad breeding, feeding, and locomotion (Semeniuk et al. 2007;Kearney et al. 2008), cane toads are expanding their range far less rapidly in this region than at the tropical invasion front (1-3 km/year vs. 55 km/year : Seabrook 1993;Phillips et al. 2006Phillips et al. , 2007Urban et al. 2007). Although toads have been in NSW for almost 50 years (van Beurden and Grigg 1980), they remain patchily distributed (Seabrook 1993). That patchy distribution likely reflects anthropogenic habitat fragmentation, human-assisted dispersal, and heterogeneity of densely vegetated habitats (Seabrook 1993;Semeniuk et al. 2007). Hence, large tracts of suitable, but as yet  uncolonized, native bushland occur near long-term toadcolonized areas.

Study sites
We surveyed habitat characteristics and fauna at 16 campgrounds and picnic areas surrounded by bushland (from 28°22 0 S, 153°14 0 E to 29°57 0 S, 153°15 0 E) between October 2013 and February 2014. We selected eight sites in areas where toads are currently present (from 2 to 25 years postinvasion) and eight sites where toads are predicted to invade (Urban et al. 2007;Kolbe et al. 2010) but have not yet done so. Due to the patchy distribution of cane toads at the southern edge of their invasion, we were able to select interspersed study sites to minimize the confounding effects of latitude, longitude, elevation, climate, or vegetation. All campgrounds and picnic areas were adjacent to or within state forests or national parks ( Fig. 1, see Supplementary Information for Table S1). Each site consisted of a cleared campground or picnic area plus an adjacent 5-km section of access road through native bushland. Human-supplied food and water subsidies attract both native and invasive fauna to disturbed patches such as campgrounds, picnic areas, and landfill dumps located within reserves (Warnken et al. 2004;Piper and Catterall 2006;Jessop et al. 2012). Compared to the surrounding bushland, campgrounds and access roads also offer reduced canopy cover (and hence, greater sun exposure) at ground level, potentially facilitating reptile thermoregulation (Sartorius et al. 1999;Heard et al. 2004;McDonald 2012). Additionally, access roads provide corridors for animals to move between vegetation remnants (Taylor and Goldingay 2003), rendering the animals more easily observable during surveys. Cane toads regularly use roadways as dispersal corridors (Seabrook and Dettmann 1996;Brown et al. 2006) and prefer the sparse vegetation associated with human-modified areas within this region (e.g., campgrounds, golf courses, paddocks, parkland : Seabrook 1993;Semeniuk et al. 2007).

Fauna surveys
We examined the abundance and diversity of native predators using standardized surveys for a five-month period spanning the spring and summer peak period of faunal activity in this region. We recorded all reptile, bird (predatory and scavenger species), and mammal taxa that were encountered in these surveys (see Supporting Information for Table S2). Because populations of large elapid snakes, varanid and scincid lizards, and dasyurid marsupials have suffered severe declines as a result of the toad's tropical invasion, our surveys were designed to detect any toad-imposed differences in the abundance of these taxa.
We actively sampled each site to count reptiles, birds, and mammals. As sites differed in extent (see Supporting  Information for Table S1), survey effort was standardized (1 h/survey). For the first 15 min, we actively searched for fauna in and around campground areas on foot. The remaining 45 min was conducted from a slow-moving car (20-40 km/h along 5-km road transect), from which we scanned the road for crossing animals, and roadside areas (e.g., trees, branches, fallen logs, termite mounds) for sequestered animals. Data for each campground and road transect within a site were combined to give total number of species and total number of individuals of each species for each site (standardized per 15 h of survey time). All animals observed were identified to species, using field guides and keys (taxonomy follows Pizzey and Knight 2012;Wilson and Swan 2013;Van Dyck et al. 2013).
During the five-month survey period, each site was sampled in 10-h-long diurnal surveys (on sunny days, >23°C), conducted at 0900-1200 h and 1300-1800 h, and in 5-h-long nocturnal surveys (on dry nights, >16°C), conducted at 1930-2330 h. To avoid bias, we rotated the time of day that each site was surveyed from day to day. For practical reasons, groups of three adjacent sites were surveyed for at least three concurrent days, in random order with respect to each group of sites and sites within each group. For statistical analysis, data for rarely encountered species (<8 individuals across all sites) were excluded.

Habitat variables
At each site, we recorded habitat variables to determine whether any spatial differences in fauna composition might be attributable to such factors rather than to the impacts of cane toads. Twenty-eight structural habitat variables (see Supporting Information for Table S3) were measured following a protocol similar to that of Brown et al. (2008). Beginning at the center of each campground or picnic area, we marked out two transects (north-south and east-west) into the surrounding bushland. From the origin point, we marked out four 10 9 10 m square plots centered at 20-m intervals along each transect line (within the cleared area of the site) and a further four plots at 10-m intervals in the surrounding bushland. The distances between adjacent plots were smaller in the bushland because our primary interest was the transitional zone and because of the greater apparent homogeneity of habitat conditions within the open areas. Within each plot, we estimated the proportions of the substrate covered by bare ground, lawn, leaf litter, rock, log, grass/ herb, small shrub, tree trunk, and human debris (e.g., graveled roads, corrugated iron). To quantify structural aspects of the vegetation, we estimated percentage cover of the understory, midstory, and canopy. Percentage of the ground within each plot that was exposed to direct sunlight at midday was also recorded. We noted the numbers of small (1-10 cm diameter at breast height [dbh]), medium (11-30 cm dbh), and large (>31 cm dbh) trees per plot, as well as mean tree height . In each plot, the following structural variables were also measured: leaf litter depth (cm); number of tree stumps (<5 m in height); number of fallen logs; number of hollows in tree trunks and limbs; and number of termite mounds in trees and on the ground. Additionally, we recorded the distance from the center of each plot to the closest building, tent, or car; water source; road, track, or trail; and tree.
Additional whole-campground variables measured at the commencement of surveys comprised: area cleared (m 2 ); number of car spaces and campsites; number of campsites occupied; number of people and cars present during surveys; number of rubbish bins; number of barbeques; number of outside taps; number of access roads and walking trails; and distance (m) to nearest permanent freshwater. Landscape-scale parameters were assessed using satellite and aerial maps of the region (Google Earth). Connectivity, disturbance, and isolation were quantified by measuring cover of wooded vegetation (km 2 ), area cleared for agricultural or residential purposes (km 2 ), length of roadways (km), and average speed limit (km/h) of roadway, all within a 5 km radius of each site. Climatic variables were extracted from BIOCLIM climate layers to predict mean annual precipitation (mm) and mean annual maximum and minimum temperatures (°C) of each site as a function of latitude, longitude, and elevation (Busby 1991). Climatic variables were checked against long-term data sets of local weather stations (Bureau of Meteorology 2014).

Rates of carrion removal
To record rates of carrion removal in sites with and without cane toads, as well as the identity of the species responsible for carrion removal, we set up eight bait stations with remotely triggered cameras at each of 12 sites. Four bait stations were located in the bushland surrounding each campground and picnic area (in bushland vegetation plots; see above), and four bait stations were located 2 km away, along the ingress road (two on either side, and within 10 m of the road). All bait stations consisted of 10 chicken necks, lightly concealed with substrate, that were monitored by motion-sensitive infrared-triggered digital video cameras (Moultrie M-990i: EBSCO Industries, Birmingham, Alabama) for 48 h. Once this period had elapsed, we counted the remaining baits and reviewed the camera footage to determine the iden-tity of the scavenging species. Only data from baits that were deployed during a 48-h period of suitable weather (clear, sunny days >23°C, and fine, warm nights >16°C) were used in comparisons between sites.

Habitat variables
Using the software package SPSS (SAS Institute, Cary, NC), we carried out a principal components analysis (PCA) to reduce the number of intercorrelated variables. Variances were homogeneous (Levene's test; P > 0.05). Bartlett's test of sphericity rejected the null hypothesis that variables were not correlated with habitats (P < 0.05), and identified 12 principal components (PCs; eigenvalue >1) in the rotated component matrix (using varimax with Kaiser normalization, accounting for >95% of variance in the habitat data). However, only five of these PCs were deemed to be biologically relevant. We used a multivariate analysis of variance (MANOVA) with toad exposure as the factor (two levels; cane toad present vs. absent), to explore determinants of variation in each PC axis (dependent variables). To test the roles of toads versus preexisting habitat conditions as determinants of wildlife assemblages, we constructed a series of multivariate regression models to explain (1) species richness and (2) faunal abundance at the 16 study sites. In these models, the predictor variables were toad invasion status (coded as either presence/absence or number of years since invasion [uninvaded sites were scored as 0]), the five PC habitat axis scores, and all interactions between these factors. We assessed model fit using the corrected Akaike's information criterion (AIC). We ran separate models with toad status treated as either a nominal variable (present/absent) or a continuous variable (years since toad arrival) because plausibly, the impact of toads might be immediate; or might change through time (either increasing or decreasing) if native fauna are affected through indirect processes (which might thus take significant time to eventuate), or are vulnerable only under specific environmental conditions.

Fauna surveys
Univariate analyses were conducted using general linear models in the statistical software JMP Pro 9.0 (SAS Institute). Data were visually assessed for normality, and variances were homogeneous (Levene's test, P > 0.05). One-factor (toad exposure) analyses of variance (ANOVA) were used to test for effects of the independent variable (two levels; cane toad present vs. absent) on total abundance and species richness of native vertebrates. We used the program PRIMER v5 (PRIMER-E, Plymouth, UK) to assess differences in diversity and composition of native vertebrate species between sites with and without cane toads. Using the Bray-Curtis similarity coefficient (Bray and Curtis 1957), we calculated a similarity matrix with data from fauna surveys (abundance of each native vertebrate species at each site). Abundance data were fourth-root-transformed to reduce the influence of common species relative to rarer species (e.g., Clarke 1993;Quinn and Keough 2002;Lassau and Hochuli 2005). Standardization of the data was unnecessary due to prior standardization of sampling design. We employed nonmetric multidimensional scaling (nMDS: Clarke and Warwick 2001), using Bray-Curtis similarity measures, to plot a two-dimensional ordination and determine whether native vertebrate assemblages differed between sites with versus without cane toads. To compare between toad-present versus toad-absent sites in a multivariate similarity matrix, we performed analysis of similarities (ANOSIM) with 999 permutations (Sokal and Rohlf 1995;Anderson 2001;Bond and Lake 2003). As the presence of cane toads contributed to differences in species composition (see Results), we used a similarity percentage analysis (SIM-PER: Clarke 1993) to identify individual species that contributed most to those differences (Clarke and Warwick 2001).
Based on those analyses, we then proceeded to compare the abundances of three species of lizard (land mullets Bellatorias major, eastern water dragons Intellagama lesueurii, lace monitors Varanus varius) and one species of snake (red-bellied black snake Pseudechis porphyriacus) between sites where toads were present versus absent. Data on the three lizard taxa were assessed independently because they were identified as having declined by SIMPER. The abundance of P. porphyriacus was assessed because this species has been anecdotally reported to experience severe toad-imposed population declines (Rayward 1974;Covacevich and Archer 1975;Fearn 2003;Phillips et al. 2003;Phillips & Fitzgerald 2004). Using JMP Pro 9.0, we compared the abundance of each species between sites where toads were present versus absent (independent variable). Because the abundance data of these individual species could not be normalized via transformation, data were analyzed using nonparametric Kruskal-Wallis analysis of variance tests (Crossland 1998;Ujvari et al. 2011;Crossland and Shine 2012).

Rates of carrion removal
Data conformed to the assumptions of normality and variance homogeneity. Using JMP Pro 9.0, we compared the number of chicken necks removed from bait stations (both in campgrounds and in surrounding bushland) in both toad-present and toad-absent sites. The number of baits removed per bait station was the dependent variable, and toad exposure (two levels; toad present vs. absent) and location (two levels; campground vs. bushland) were the independent variables in a two-way analysis of variance (ANOVA). Site was included as a random factor. Differences in numbers of goanna vs. nongoanna scavengers between toad-present and toad-absent sites were compared using Fisher's exact test.

Habitat variables
From the 49 input habitat and climatic variables, we retained five PC axes (eigenvalue >4.3; each accounting for >8.9% of variance). The first axis (PC1, 19.9% of variance) loaded on campground size. PC2 (15.2% of variance) was linked to forested sites with high human resource subsidies. PC3 (11.9% of variance) was linked to leaf litter and dense vegetation. PC4 (9.2% of variance) was associated with naturally sparse vegetation, and PC5 (8.9% of variance) was associated with warm, woody sites. A MANOVA with toad presence/absence as the factor, and scores on these PC axes as dependent variables, showed no overall habitat differences between the two groups of sites (MANOVA: Wilks' Lambda = 1.63, P > 0.05).
Of the models we constructed to explain species richness and faunal abundance at the 16 study sites, the bestfitting models all included the impact of cane toads rather than (or as well as) habitat variation. If we included toad presence/absence as a factor, the best-fitting model for species richness included toad presence/absence and PC axis 4; the second-best (not significantly different) model included only toad presence/absence (Table 1). The results for faunal abundance were the same, except that PC1 took the role of PC4 (Table 2). If we included the impact of toads as a continuous variable (years since arrival) rather than a dichotomy, faunal species richness was best explained by a model that included only time since toad arrival and PC3 (Table 3). For faunal abundance, the top three models were equally well supported (DAICc < 2), and all included time since toad arrival (one also contained PC1, one contained PC2, and one contained both PC1 and PC2; see Table 4). Thus, all of the models to explain variation among sites in faunal richness and abundance included the effect of cane toads, even after the effects of habitat variation were taken into account.

Discussion
Invasive cane toads have massively affected the abundance and species richness of native fauna in temperate-zone Australia. In sites where toads were present, we recorded 40% fewer species and 31% fewer individuals. A causal role for toad invasion in faunal decline was supported by (1) the similarity in habitat structure and climate between toad-present versus toad-absent areas; and (2) the nature of faunal differences: The species that were less common in toad-invaded areas were reptile taxa that eat anurans. In contrast, the ground-dwelling brush turkey increased in abundance, consistent with reduced predation by goannas (Jones 1988;Goth and Vogel 2002). Our data thus suggest both direct and indirect impacts of cane toad invasion. Below, we consider native species that (1) were not affected; (2) decreased coincident with toads; and (3) increased in toad-invaded areas.
Most species were equally common in sites with and without cane toads with no negative effect on any surveyed native mammal or bird taxa. Studies on toad impact in tropical Australia have reported the same general result. Most Australian birds that prey on anurans can detect and tolerate the toad's toxin Shine 2009, 2011;. Although one of the bird species we counted can be killed by ingesting toads (laughing kookaburra Dacelo novaeguineae: Covacevich and Archer 1975; but see Ringma 2013), this vulnerability has not translated to population-level impacts. Broadly, then, direct impacts of cane toads on predatory birds appear to be minimal, as in tropical Australia.
Similarly, we found no significant changes in mammal abundance associated with toad presence. Only three of the ten species we surveyed are likely to consume toads (fawn-footed melomys Melomys cervinipes; bush rat Rattus fuscipes; feral cat Felis catus). Rodents can kill and consume cane toads without ill effects (Cassels 1966;Cabrera-Guzm an et al. 2015), reflecting ancestral exposure to Asian bufonids (Fitzgerald 1990;Shine 2010). The remaining seven mammal species we surveyed (see Supporting Information for Table S2) do not eat anurans; thus, any impact to these species would have been indirect (e.g., through a reduction in goanna predation). Small sample sizes (number of individuals encountered) weakened our ability to detect any changes, and further research would be valuable.
Our surveys failed to include some rare mammal species that might be affected by toads. In tropical Australia, the northern quoll (Dasyurus hallucatus) has experienced severe declines from toad invasion (Burnett 1997;O'Donnell et al. 2010;Shine 2010). Two related dasyurid species (spotted-tailed quoll Dasyurus maculatus and brush-tailed phascogale Phascogale tapoatafa) are found in northeastern (NSW) (Van Dyck et al. 2013). However, both are rare and neither was encountered during our surveys. We saw a single P. tapoatafa at one toad-free site, but not during standardized surveys.
Our data do not reveal any overall declines in snake abundance. Phillips et al. (2003) predicted that cane toads  would cause population declines in 30% of Australia's snake species (including eight of nine species recorded in the present study), but our results suggest a more encouraging scenario. Most snake taxa, even frog specialist species (e.g., Dendrelaphis punctulata, Tropidechis carinatus) predicted to suffer toad-imposed impacts (Phillips et al. 2003), were unaffected. That lack of effect may be due to small sample sizes, to morphological or behavioral traits that render the snakes invulnerable to toads (Llewelyn et al. 2012), or to indirect positive effects (e.g., of goanna mortality) that outweigh any direct negative effects (Brown et al. 2011(Brown et al. , 2013bDoody et al. 2013).
Although most taxa were unaffected, toads appear to have caused catastrophic declines in populations of four reptile taxa in temperate Australia. The three lizard species affected have not been studied previously in this respect, but toad-induced declines in red-bellied blacksnakes have been reported (Pockley 1965;Rayward 1974), based on anecdotal evidence (but see Seabrook 1993). This species was absent from our toad-invaded sites. Both lace monitors and water dragons are widely sympatric with cane toads in Australia (Lever 2001;Wilson and Swan 2013), so the lack of prior reports of toad impact on these taxa is surprising.
Although these reptile species span four phylogenetically diverse families (Varanidae, Agamidae, Scincidae, and Elapidae), they share three attributes common to predators that are vulnerable to toads (Shine 2010): 1 Low resistance to toad toxins. Lace monitors, land mullets, and red-bellied blacksnakes can be fatally poisoned by ingesting toads (Covacevich and Archer 1975;Fearn 2003;Phillips and Shine 2006a,b;Ujvari et al. 2013). Toxin resistance is unknown for water dragons, but is low in related species (e.g., Chlamydosaurus kingii: Pearson et al. 2014). 2 Anurophagy. All of these species eat anurans (Shine 1977;Phillips and Shine 2006a,b) or are generalists with broad diets (Shea 1999;Wilson 2012). A close relative of the land mullet, the blue-tongued lizard (Tiliqua scincoides intermedia), has shown dramatic population declines due to toad invasion in tropical Australia size (Phillips and Shine 2006a,b), larger predators (which attack larger toads) are most at risk (Shine 2010).
In tropical Australia, the decline of large predatory lizards (e.g., Varanus panoptes) means that some native species indirectly benefit from the toad's arrival (Brown  , 2013aDoody et al. 2013). Similarly, a reduction in the abundance of varanid lizards may explain the increase in the number of brush turkeys (Alectura lathami) in our toad-present sites. The ground-nesting habits of brush turkeys may render them (and their eggs) highly vulnerable to predation by monitors (Jones 1988;Goth and Vogel 2002). In addition to increasing brush turkey abundance, the toad-induced decline in lace monitors results in an indirect reduction to scavenging rates by 74% in toad-invaded sites. That situation also may allow increases in other scavengers (such as insects) and might create health problems for human users of the campgrounds if discarded food is left to rot. Scavengers can strongly influence the structure and function of faunal assemblages (Wilson and Wolkovich 2011).
What can be concluded, overall, about the impacts of the invasive cane toad (Rhinella marina) on the abundance, diversity, and composition of Australian wildlife? Although brief, our study incorporated more spatial replication and "control" sites than have been possible in the wet-dry tropics (Doody et al. 2006(Doody et al. , 2009(Doody et al. , 2013Letnic et al. 2008;Phillips et al. 2010;Brown et al. 2011Brown et al. , 2013bUjvari et al. 2011). The ability to compare faunal assemblages between areas that differ mostly in exposure to toads (rather than in environmental factors) provides convincing evidence of population-level declines in native taxa due to toad arrival. Our results support the conclusion that toad impacts on native taxa are complex and proceed via both direct and indirect pathways (Shine 2010;Brown et al. 2013a). Additionally, our study provides evidence of invader impact in a region that heretofore has been largely ignored.
Cane toads are publicly vilified Australia-wide, and the public views this invasive anuran with abhorrence (Clarke et al. 2009;Shine 2010). That passion has encouraged extensive research, as well as community-based "toad busting" (Shine and Doody 2011). Given this high public profile and substantial governmental investment into cane toad issues , the lack of previous research on toad impacts in southern Australia is remarkable. Our study paints a bleak cautionary tale; even in an intensively studied invasive species system, major impacts on iconic native species have been overlooked. Remarkably, even anecdotal reports of cane toad impacts are rare in southern Australia, despite the dense human population. Nonchalance among the general public appears to have had flow-on effects for political and research priorities.
The scientific neglect of the toad's impacts at the southern invasion front stands in contrast to the situation in tropical Australia. Many topics related to toad biology  and ecological impacts have been investigated in remote tropical regions (review by Shine 2010). Between 1963 and 2014, 102 scientific publications assessed the ecological impacts of cane toads on the native fauna of northern Australia (i.e., north of Brisbane) compared to 14 articles on cane toad impacts in southern Australian (Shine 2010(Shine , 2014. Of those 14 articles, the most recent was completed more than 20 years ago (Seabrook 1993). This disparity in research effort is perplexing in light of the logistical difficulties that have compromised experimental designs in these tropical studies.
In summary, we found that toad invasion has caused population declines in some large anurophagous predators in temperate Australia similar to those documented in the more intensively studied tropics (Shine 2010). Ultimately, the neglect of toad impact in temperate-zone Australia reflects the fact that debates about invasive species occur within a sociopolitical framework. Geographic differences in the priority given to competing issues can influence funding incentives and research effort. The forests of northern NSW have been the focus of vigorous "environmental" battles over many years relating to forestry practices (Lemckert 1999;Kavanagh and Stanton 2005). That focus may have drawn public attention away from feral species impacts, but it is difficult to understand why scientific efforts were equally scarce. Whatever the reasons behind the disproportionate allocation of research effort toward the cane toad problem in tropical versus southern Australia, the result is clear. We have failed to recognize a major ecological problem unfolding in a place close to major cities where logistics are straightforward, and robust experimental designs are possible. If that can happen with cane toad impact, it may well happen in other ecological issues as well.

Supporting Information
Additional Supporting Information may be found in the online version of this article: Table S1. Study sites located in the Northern Rivers region, New South Wales, Australia. Table S2. Vertebrate taxa recorded during surveys of campgrounds in northeastern New South Wales, Australia. Table S3. Habitat attributes of campgrounds and bushland in toad-present versus toad-absent sites.