Testing for multiple invasion routes and source populations for the invasive brown treesnake (Boiga irregularis) on Guam: implications for pest management

Mapping environmental crime to characterize human impacts on islands: an applied and methodological research in Canary Islands

Environmental crimes are a global issue due to the damage they cause to landscapes and ecosystems. This study focused on characterizing environmental crimes in the Canary Islands (Spain). Four categories of environmental crimes related to construction, mining and tilling, solid waste, and liquid waste) were defined and analysed. A total of 28 databases were generated, corresponding to each of the 7 major islands and each environmental crime typology. Each database was linked to information on land use and the socioeconomic and physical characteristics of the territory. For each database, firstly a descriptive statistical analysis was conducted, followed by the generation of a regularized Random Forest model with the aim of identifying characteristics that may be related to the location of environmental crimes. The results showed that, in most cases, proximity to residential accommodations, agricultural areas and industrial zones act as the main explanatory features of the distribution of environmental crimes. Furthermore, a marked pattern of concentration of environmental crimes in the coastal belt of the islands was observed, mainly associated with urban-tourist development since the 1960s and 1970s.

Evaluation of camera trap-based abundance estimators for unmarked populations

Estimates of species abundance are critical to understand population processes and to assess and select management actions. However, capturing and marking individuals for abundance estimation, while providing robust information, can be economically and logistically prohibitive, particularly for species with cryptic behavior. Camera traps can be used to collect data at temporal and spatial scales necessary for estimating abundance, but the use of camera traps comes with limitations when target species are not uniquely identifiable (i.e., "unmarked"). Abundance estimation is particularly useful in the management of invasive species, with herpetofauna being recognized as some of the most pervasive and detrimental invasive vertebrate species. However, the use of camera traps for these taxa presents additional challenges with relevancy across multiple taxa. It is often necessary to use lures to attract animals in order to obtain sufficient observations, yet lure attraction can influence species' landscape use and potentially induce bias in abundance estimators. We investigated these challenges and assessed the feasibility of obtaining reliable abundance estimates using camera-trapping data on a population of invasive brown treesnakes (Boiga irregularis) in Guam. Data were collected using camera traps in an enclosed area where snakes were subject to high-intensity capture-recapture effort, resulting in presumed abundance of 116 snakes (density = 23/ha). We then applied spatial count, random encounter and staying time, space to event, and instantaneous sampling estimators to photo-capture data to estimate abundance and compared estimates to our presumed abundance. We found that all estimators for unmarked populations performed poorly, with inaccurate or imprecise abundance estimates that limit their usefulness for management in this system. We further investigated the sensitivity of these estimators to the use of lures (i.e., violating the assumption that animal behavior is unchanged by sampling) and camera density in a simulation study. Increasing the effective distances of a lure (i.e., lure attraction) and camera density both resulted in biased abundance estimates. Each estimator rarely recovered truth or suffered from convergence issues. Our results indicate that, when limited to unmarked estimators and the use of lures, camera traps alone are unlikely to produce abundance estimates with utility for brown treesnake management.

Influence of niche breadth and position on the historical biogeography of seafaring scincid lizards

Niche breadth and position can influence diversification among closely related species or populations, yet limited empirical data exist concerning the predictability of the outcomes. We explored the effects of these factors on the evolution of the Emoia atrocostata species group, an insular radiation of lizards in the western Pacific Ocean and Indo-Australasia composed of both endemic and widespread species that differ in niche occupancy. We used molecular data and phylogeographical diffusion models to estimate the timing and patterns of range expansion, and ancestral reconstruction methods to infer shifts in ecology. We show evidence of multidirectional spread from a centre of origin in western Micronesia, and that the phyletic diversity of the group is derived from a putative habitat specialist that survives in the littoral zone. This species is composed of paraphyletic lineages that represent stages or possible endpoints in the continuum toward speciation. Several descendant species have transitioned to either strand or interior forest habitat, but only on remote islands with depauperate terrestrial faunas. Our results suggest that the atrocostata group might be in the early phases of a Wilsonian taxon cycle and that the capacity to tolerate salt stress has promoted dispersal and colonization of remote oceanic islands. Divergence itself, however, is largely driven by geographical isolation rather than shifts in ecology.

Phylogenetic relationships and biogeographic range evolution in cat-eyed snakes, Boiga (Serpentes: Colubridae)

The genus Boiga includes 35, primarily arboreal snake species distributed from the Middle East to Australia and many islands in the western Pacific, with particularly high species diversity in South-East Asia. Despite including the iconic mangrove snakes (Boiga dendrophila complex) and the brown tree snake (Boiga irregularis; infamous for avian extinctions on small islands of the Pacific), species-level phylogenetic relationships and the biogeographic history of this ecologically and morphologically distinct clade are poorly understood. In this study, we sequenced mitochondrial and nuclear DNA for 24 Boiga species and used these data to estimate a robust phylogenetic inference, in order to (1) test the hypothesis that Boiga is monophyletic, (2) evaluate the validity of current species-level taxonomy and (3) examine whether geographic range evolution in Boiga is consistent with expectations concerning dispersal and colonization of vertebrates between continents and islands. Our results support the prevailing view that most dispersal events are downstream – from continents to oceanic islands – but we also identify a role for upstream dispersal from oceanic islands to continents. Additionally, the novel phylogeny of Boiga presented here is informative for updating species-level taxonomy within the genus.

Genomic pedigree reconstruction identifies predictors of mating and reproductive success in an invasive vertebrate

The persistence of an invasive species is influenced by its reproductive ecology, and a successful control program must operate on this premise. However, the reproductive ecology of invasive species may be enigmatic due to factors that also limit their management, such as cryptic coloration and behavior. We explored the mating and reproductive ecology of the invasive Brown Treesnake (BTS: Boiga irregularis) by reconstructing a multigenerational genomic pedigree based on 654 single nucleotide polymorphisms for a geographically closed population established in 2004 on Guam (N = 426). The pedigree allowed annual estimates of individual mating and reproductive success to be inferred for snakes in the study population over a 14-year period. We then employed generalized linear mixed models to gauge how well phenotypic and genomic data could predict sex-specific annual mating and reproductive success. Average snout-vent length (SVL), average body condition index (BCI), and trappability were significantly related to annual mating success for males, with average SVL also related to annual mating success for females. Male and female annual reproductive success was positively affected by SVL, BCI, and trappability. Surprisingly, the degree to which individuals were inbred had no effect on annual mating or reproductive success. When juxtaposed with current control methods, these results indicate that baited traps, a common interdiction tool, may target fecund BTS in some regards but not others. Our study emphasizes the importance of reproductive ecology as a focus for improving BTS control and promotes genomic pedigree reconstruction for such an endeavor in this invasive species and others.

Behavioral differences following ingestion of large meals and consequences for management of a harmful invasive snake: A field experiment

Many snakes are uniquely adapted to ingest large prey at infrequent intervals. Digestion of large prey is metabolically and aerobically costly, and large prey boluses can impair snake locomotion, increasing vulnerability to predation. Cessation of foraging and use of refugia with microclimates facilitating digestion are expected to be strategies employed by free-ranging snakes to cope with the demands of digestion while minimizing risk of predation. However, empirical observations of such submergent behavior from field experiments are limited. The brown treesnake (Serpentes: Colubridae: Boiga irregularis) is a nocturnal, arboreal, colubrid snake that was accidentally introduced to the island of Guam, with ecologically and economically costly consequences. Because tools for brown treesnake damage prevention generally rely on snakes being visible or responding to lures or baits while foraging, cessation of foraging activities after feeding would complicate management. We sought to characterize differences in brown treesnake activity, movement, habitat use, and detectability following feeding of large meals (rodents 33% of the snake's unfed body mass) via radio telemetry, trapping, and visual surveys. Compared to unfed snakes, snakes in the feeding treatment group showed drastic decreases in hourly and nightly activity rates, differences in refuge height and microhabitat type, and a marked decrease in detectability by trapping and visual surveys. Depression of activity lasted approximately 5-7 days, a period that corresponds to previous studies of brown treesnake digestion and cycles of detectability. Our results indicate that management strategies for invasive brown treesnakes need to account for cycles of unavailability and underscore the importance of preventing spread of brown treesnakes to new environments where large prey are abundant and periods of cryptic behavior are likely to be frequent. Characterization of postfeeding behavior changes provides a richer understanding of snake ecology and foraging models for species that consume large prey.

Inferring the absence of an incipient population during a rapid response for an invasive species

Successful eradication of invasives is facilitated by early detection and prompt onset of control. However, realizing or verifying that a colonization has occurred is difficult for cryptic species especially at low population densities. Responding to the capture or unconfirmed sighting of a cryptic invasive species, and the associated effort to determine if it indicates an incipient (small, localized) population or merely a lone colonizer, is costly and cannot continue indefinitely. However, insufficient detection effort risks erroneously concluding the species is not present, allowing the population to increase in size and expand its range. Evidence for an incipient population requires detection of ≥1 individual; its absence, on the other hand, must be inferred probabilistically. We use an actual rapid response incident and species-specific detection estimates tied to a known density to calculate the amount of effort (with non-sequential detections) necessary to assert, with a pre-defined confidence, that invasive brown treesnakes are absent from the search area under a wide range of hypothetical population densities. We illustrate that the amount of effort necessary to declare that a species is absent is substantial and increases with decreased individual detection probability, decreased density, and increased level of desired confidence about its absence. Such survey investment would be justified where the cost savings due to early detection are large. Our Poisson-based model application will allow managers to make informed decisions about how long to continue detection efforts, should no additional detections occur, and suggests that effort to do so is significantly higher than previously thought. While our model application informs how long to search to infer absence of an incipient population of brown treesnakes, the approach is sufficiently general to apply to other invasive species if density-dependent detection estimates are known or reliable surrogate estimates are available.

Objectives and integrated approaches for the control of brown tree snakes: An updated overview

After its inadvertent introduction to Guam, the brown tree snake (Boiga irregularis, BTS) extirpated most of the island's native terrestrial vertebrates, presented a health hazard to small children, and had considerable economic ramifications. Management of BTS is aimed at a number of objectives, the foremost of which has been to deter its dispersal from Guam to other locations. Further objectives include reclaiming areas on Guam as snake-free for reintroduction of native wildlife, to protect small sensitive sites on Guam from BTS intrusion (e.g, power stations, bird nesting sites), to contain and capture stowaway BTS incoming to vulnerable destinations, and to control incipient populations in other areas beyond their native range. A number of control tools have been developed, and the efficacy of each control method depends on the situation to which it is applied. Integration of control methods provides the most efficacious results for all objectives. Here, we outline the different objectives for managing BTS, and the tools and methods available for BTS management. We complete the picture by describing which tools and methods are best suited to accomplish each management objective.

Transcriptomics-guided bottom-up and top-down venomics of neonate and adult specimens of the arboreal rear-fanged Brown Treesnake, Boiga irregularis, from Guam

The Brown Treesnake (Boiga irregularis) is an arboreal, nocturnal, rear-fanged venomous snake native to northern and eastern regions of Australia, Papua New Guinea and the Solomon Islands. It was inadvertently introduced onto the island of Guam during the late 1940's to early 1950's, and it has caused massive declines and extirpations of the native bird, lizard, and mammal populations. In the current study, we report the characterization of the venom proteome of an adult and a neonate B. irregularis specimens from Guam by a combination of venom gland transcriptomic and venomic analyses. Venom gland transcriptomic analysis of an adult individual identified toxins belonging to 18 protein families, with three-finger toxin isoforms being the most abundantly expressed transcripts, comprising 94% of all venom protein transcript reads. Transcripts for PIII-metalloproteinases, C-type lectins, cysteine-rich secretory proteins, acetylcholinesterases, natriuretic peptides, ficolins, phospholipase A₂ (PLA₂) inhibitors, PLA₂s, vascular endothelial growth factors, Kunitz-type protease inhibitors, cystatins, phospholipase Bs, cobra venom factors, waprins, SVMP inhibitors, matrix metalloproteinases, and hyaluronidases were also identified, albeit, at very low abundances ranging from 0.05% to 1.7% of the transcriptome. The venom proteomes of neonate and adult B. irregularis were also both overwhelmingly (78 and 84%, respectively) dominated by monomeric and dimeric 3FTxs, followed by moderately abundant (21% (N) and 13% (A)) CRISPs, low abundance (1% (N), 3% (A)) PIII-SVMPs, and very low abundance (<0.01%) PLA₂ and SVMP inhibitors. The differences in relative toxin abundances identified between neonate and adult snakes likely correlates to shifts in prey preference between the two age classes, from nearly-exclusively lizards to lizards, birds and small mammals. Immunoaffinity antivenomics with experimentally designed rabbit anti-Brown Treesnake (anti-BTS) venom IgGs against homologous venom from adult snakes demonstrated that CRISPs, PIII-SVMPs, and 60-70% of 3FTxs were effectively immunocaptured. Western blot analysis showed that all venom proteins were recognized by anti-BTS IgGs, and cross-reactivity with other rear-fanged snake venoms was also observed. Incubation of anti-BTS venom IgGs with crude B. irregularis venom resulted in a significant decrease in proteolytic (SVMP) activity against azocasein. These results provide the first comparative venomic and anti-venomic analysis of neonate and adult B. irregularis from Guam, further highlighting evolutionary trends in venom composition among rear-fanged venomous snakes.

SIGNIFICANCE PARAGRAPH

The Brown Treesnake (Boiga irregularis) has caused extensive ecological and economic damage to the island of Guam where it has become a classic example of the negative impacts of invasive species. In the current study, we report the first combined transcriptomic and proteomic analysis of B. irregularis venom of Guam origin. The transcriptome of an adult snake contained toxin sequences belonging to 18 protein families, with three-finger toxin (3FTx) isoforms being the most abundant and representing 94% of all venom protein transcript reads. Our bottom-up and top-down venomic analyses confirmed that 3FTxs are the major components of B. irregularis venom, and a comparative analysis of neonate and adult venoms demonstrate a clear ontogenetic shift in toxin abundance, likely driven by dietary variation between the two age classes. Second-generation antivenomics and Western blot analysis using purified anti-Brown Treesnake rabbit serum IgGs (anti-BTS IgGs) showed strong immunoreactivity toward B. irregularis venom. Interestingly, our anti-BTS IgGs did not cross-react with 3FTxs found in several other rear-fanged snake venoms, or against 3FTxs in the venom of the elapid Ophiophagus hannah, indicating that epitopes in these 3FTx molecules are quite distinct.

Genetic reconstruction of a bullfrog invasion to elucidate vectors of introduction and secondary spread

Reconstructing historical colonization pathways of an invasive species is critical for uncovering factors that determine invasion success and for designing management strategies. The American bullfrog (Lithobates catesbeianus) is endemic to eastern North America, but now has a global distribution and is considered to be one of the worst invaders in the world. In Montana, several introduced populations have been reported, but little is known of their sources and vectors of introduction and secondary spread. We evaluated the genetic composition of introduced populations at local (Yellowstone River floodplain) and regional (Montana and Wyoming) scales in contrast to native range populations. Our objectives were to (1) estimate the number of introductions, (2) identify probable native sources, (3) evaluate genetic variation relative to sources, and (4) characterize properties of local‐ and regional‐scale spread. We sequenced 937 bp of the mitochondrial cytochrome b locus in 395 tadpoles collected along 100 km of the Yellowstone River, from three additional sites in MT and a proximate site in WY. Pairwise Φ_ST revealed high divergence among nonnative populations, suggesting at least four independent introductions into MT from diverse sources. Three cyt b haplotypes were identical to native haplotypes distributed across the Midwest and Great Lakes regions, and AMOVA confirmed the western native region as a likely source. While haplotype (H_d = 0.69) and nucleotide diversity (π = 0.005) were low in introduced bullfrogs, the levels of diversity did not differ significantly from source populations. In the Yellowstone, two identified haplotypes implied few introduction vectors and a significant relationship between genetic and river distance was found. Evidence for multiple invasions and lack of subsequent regional spread emphasizes the importance of enforcing legislation prohibiting bullfrog importation and the need for continuing public education to prevent transport of bullfrogs in MT. More broadly, this study demonstrates how genetic approaches can reveal key properties of a biological invasion to inform management strategies.

Invaded Invaders: Infection of Invasive Brown Treesnakes on Guam by an Exotic Larval Cestode with a Life Cycle Comprised of Non-Native Hosts

Background

Multiple host introductions to the same non-native environment have the potential to complete life cycles of parasites incidentally transported with them. Our goal was to identify a recently detected parasitic flatworm in the invasive Brown Treesnake (Boiga irregularis) on the remote Pacific island of Guam. We considered possible factors influencing parasite transmission, and tested for correlations between infection status and potential indicators of host fitness. We used genetic data from the parasite and information about the native ranges of other possible non-native hosts to hypothesize how it arrived on Guam and how its life cycle may be currently supported.

Methods

We identified the parasite by comparing larval morphology and mtDNA sequences with other Pseudophyllid tapeworms. We assessed probability of infection in individual snakes using logistic regression and examined different factors influencing presence of parasites in hosts.

Results

We identified the parasite as the pseudophyllid cestode Spirometra erinaceieuropaei, with all sampled worms from multiple snakes sharing a single mtDNA haplotype. Infection appears to be limited to the only freshwater watershed on the island, where infection prevalence was high (77.5%). Larger snakes had a higher probability of being infected, consistent with the chronic nature of such infections. While infection status was positively correlated with body condition, infected snakes tended to have lower intra-peritoneal fat body mass, potentially indicating a negative effect on energy stores.

Conclusions

We discovered that B. irregularis inhabiting a small area of forested habitat in a freshwater watershed on Guam are often infected by a novel parasite of Asian origin. While further work is needed, this species of Spirometra, itself a non-native species, likely depends on a suite of recently introduced hosts from different parts of the world to complete the life cycle. This baseline study provides little evidence of any effects on host fitness, but additional data are needed to more thoroughly explore the consequences of infection in this invasive snake population.

Development and characterization of 33 novel polymorphic microsatellite markers for the brown tree snake Boiga irregularis

BACKGROUND

Boiga irregularis is a widespread invasive species on Guam and has led to extirpation of most of the island's native avifauna. There are presently no microsatellite markers for this invasive species, hence we developed highly polymorphic microsatellite markers to allow for robust population genetic studies on Guam.

FINDINGS

We isolated and characterized 33 microsatellite loci for the brown tree snake, B. irregularis. The loci were screened across 32 individuals from Guam. The number of alleles per locus ranged from three to ten, with an average of 4.62. The expected (He) and observed heterozygosity (Ho) ranged from 0.294 to 0.856 and from 0.031 to 0.813, with an average of 0.648 and 0.524, respectively. Significant deviations from Hardy-Weinberg equilibrium were detected at seven loci after Bonferoni correction. Probability of identity values ranged from 0.043 to 0.539.

CONCLUSIONS

These genetic markers are useful for understanding a suite of post-invasion population genetic parameters, sources of invasions, and effectiveness of management strategies for this invasive species.