Differential population responses of native and alien rodents to an invasive predator, habitat alteration and plant masting

Adaptation to bioinvasions: When does it occur?

The presence of alien species represents a major cause of habitat degradation and biodiversity loss worldwide, constituting a critical environmental challenge of our time. Despite sometimes experiencing reduced propagule pressure, leading to a reduced genetic diversity and an increased chance of inbreeding depression, alien invaders are often able to thrive in the habitats of introduction, giving rise to the so-called "genetic paradox" of biological invasions. The adaptation of alien species to the new habitats is therefore a complex aspect of biological invasions, encompassing genetic, epigenetic, and ecological processes. Albeit numerous studies and reviews investigated the mechanistic foundation of the invaders' success, and aimed to solve the genetic paradox, still remains a crucial oversight regarding the temporal context in which adaptation takes place. Given the profound knowledge and management implications, this neglected aspect of invasion biology should receive more attention when examining invaders' ability to thrive in the habitats of introduction. Here, we discuss the adaptation mechanisms exhibited by alien species with the purpose of highlighting the timing of their occurrence during the invasion process. We analyze each stage of the invasion separately, providing evidence that adaptation mechanisms play a role in all of them. However, these mechanisms vary across the different stages of invasion, and are also influenced by other factors, such as the transport speed, the reproduction type of the invader, and the presence of human interventions. Finally, we provide insights into the implications for management, and identify knowledge gaps, suggesting avenues for future research that can shed light on species adaptability. This, in turn, will contribute to a more comprehensive understanding of biological invasions.

Climate match is key to predict range expansion of the world's worst invasive terrestrial vertebrates

Understanding the determinants of the range expansion of invasive alien species is crucial for developing effective prevention and control strategies. Nevertheless, we still lack a global picture of the potential factors influencing the invaded range expansion across taxonomic groups, especially for the world's worst invaders with high ecological and economic impacts. Here, by extensively collecting data on 363 distributional ranges of 19 of world's worst invasive terrestrial vertebrates across 135 invaded administrative jurisdictions, we observed remarkable variations in the range expansion across species and taxonomic groups. After controlling for taxonomic and geographic pseudoreplicates, model averaging analyses based on generalized additive mixed-effect models showed that species in invaded regions having climates more similar to those of their native ranges tended to undergo a larger range expansion. In addition, as proxies of propagule pressure and human-assisted transportation, the number of introduction events and the road network density were also important predictors facilitating the range expansion. Further variance partitioning analyses validated the predominant role of climate match in explaining the range expansion. Our study demonstrated that regions with similar climates to their native ranges could still be prioritized to prevent the spread of invasive species under the sustained global change.

Genetic relationships among populations of the small Indian mongoose (Urva auropunctata) introduced in Japan

The small Indian mongoose, Urva auropunctata, is one of the most notorious invasive species in the world. In Japan, the mongoose was introduced to Amami, Kyushu (Kiire region), and Okinawa Islands. Currently, an eradication program for the Amami population is proving successful. To prevent reinvasion and conserve biodiversity, it would be advantageous to investigate the origin and dispersal history of introduced mongoose populations in Japan. We determined the mitochondrial genome of a mongoose from the Okinawa population and confirmed that the systematic status of the mongoose is Urva auropunctata. We also found two haplotypes of cytochrome b in the Amami population, one of which is common among the Okinawa, Amami, and Kiire populations. Based on our microsatellite analysis, the three populations were genetically different from each other. Furthermore, with a low pairwise FST value, the Kiire population was closer to the Okinawa than the Amami population. Our results suggest that the Okinawa population could be the source of the Amami and Kiire populations. These results align with historical records of mongoose in Japan and clarify its genetic status. Our results should aid the prediction of the pathway through which a new invasion will likely occur.

Population fluctuations and synanthropy explain transmission risk in rodent-borne zoonoses

Population fluctuations are widespread across the animal kingdom, especially in the order Rodentia, which includes many globally important reservoir species for zoonotic pathogens. The implications of these fluctuations for zoonotic spillover remain poorly understood. Here, we report a global empirical analysis of data describing the linkages between habitat use, population fluctuations and zoonotic reservoir status in rodents. Our quantitative synthesis is based on data collated from papers and databases. We show that the magnitude of population fluctuations combined with species' synanthropy and degree of human exploitation together distinguish most rodent reservoirs at a global scale, a result that was consistent across all pathogen types and pathogen transmission modes. Our spatial analyses identified hotspots of high transmission risk, including regions where reservoir species dominate the rodent community. Beyond rodents, these generalities inform our understanding of how natural and anthropogenic factors interact to increase the risk of zoonotic spillover in a rapidly changing world.

Development and characterization of microsatellite markers in the small Indian mongoose (Urva auropunctata)

BACKGROUND

The small Indian mongoose (Urva auropunctata) is one of the world's worst invasive alien species and eradication programs are ongoing worldwide. The development of individual and sex identification markers will improve their management.

METHODS AND RESULTS

We searched for novel mongoose microsatellite markers using genome-wide screening and identified 115,265 tetra-nucleotide repeat loci. Of 96 loci tested, 17 were genotyped in 28 mongooses from the Okinawa population. The genetic diversity analysis showed that the average expected and observed heterozygosity and number of alleles were 0.55, 0.56, and 2.94, respectively. Of 17 loci, one deviated from Hardy-Weinberg equilibrium and six loci pairs were likely linked to each other. However, we succeed in identifying all individuals using all of the microsatellite loci. The novel sex identification markers worked successfully in a test using sex known samples.

CONCLUSION

Our novel microsatellite and sex identification markers should be useful in studies of individual identification and population genetics of the mongoose.

First synthesis of the economic costs of biological invasions in Japan

Despite the large body of knowledge recognising the impact of biological invasions on biodiversity, their economic impact has been less evaluated. However, the associated economic costs ought to provide useful information on many different aspects to prevent and manage invasions. Here, we describe the economic costs of biological invasions in Japan using InvaCost, a recently-published global database on monetary costs extracted from English and non-English sources, as well as a complementary search, thereby filling a gap in regional knowledge. We focused on the following four dimensions when analysing the economic costs of biological invasions: damage to biodiversity, damage to human livelihood, management for biodiversity and management for human livelihood. Interestingly, there was no information about biological invasion costs for Japan in English, but the Japanese search and our additional survey provided a total of 630 cost entries, with a total economic cost of 728 million USD (2017 value, equivalent to 62 billion JPY). These entries appeared in 33 documents and corresponded to a total of 54 species. We showed that: 1) damage costs from biological invasions tend not to be assessed as frequently as management costs and are more underestimated; 2) despite the numerous entries, an overwhelmingly limited amount of the management budget was allocated to biodiversity conservation compared to protecting human livelihood; 3) budgets have been intensively invested in invasive species management on small islands, which reflects the vulnerability of small island ecosystems and economies to biological invasions; 4) the recorded costs still seem to be greatly underestimated, mainly due to the lack of recording (and potentially limited access to recorded cost information). These findings are not only specific to Japan, but may also be widely applicable to most other countries. The future recording of economic costs will help to close the gap between actual and recorded costs, leading to more realistic guidelines for tackling biological invasions.

Contextual behavioural plasticity in Italian agile frog (Rana latastei) tadpoles exposed to native and alien predator cues

Predation is a strong driver for the evolution of prey behaviour. To properly assess the actual risk of predation, anuran tadpoles mostly rely on water-borne chemical cues, and their ability to evaluate environmental information is even more crucial when potential predators consist of unknown alien species. Behavioural plasticity - that is, the capacity to express changes in behaviour in response to different environmental stimuli - is crucial to cope with predation risk. We explored the defensive behaviour of Italian agile frog (Rana latastei) tadpoles when exposed to the chemical cues of two predator species, one native (dragonfly larvae) and one alien (red swamp crayfish). Firstly, we observed whether a plastic life history trait (i.e. hatching time) might be affected by native predatory cues. Secondly, we recorded a suite of behavioural responses (activity level, lateralization and sinuosity) to each cue. For assessing lateralization and sinuosity, we developed a C++ code for the automatic analysis of digitally recorded tadpole tracks. Hatching time seemed not to be affected by the potential risk of predation, while both predator species and diet affected tadpoles' defensive behaviour. Tadpoles responded to a predator threat by two main defensive strategies: freezing and 'zig-zagging'. While the first behaviour had previously been reported, the analysis of individual trajectories indicated that tadpoles can also increase path complexity, probably to prevent predators from anticipating their location. We also recorded a decrease in lateralization intensity, which suggests that under predation risk, tadpoles tend to scrutinize the surrounding environment equally on both sides.

Spatial optimization of invasive species control informed by management practices

Optimization of spatial resource allocation is crucial for the successful control of invasive species under a limited budget but requires labor-intensive surveys to estimate population parameters. In this study, we devised a novel framework for the spatially explicit optimization of capture effort allocation using state-space population models from past capture records. We applied it to a control program for invasive snapping turtles to determine effort allocation strategies that minimize the population density over the whole area. We found that spatially heterogeneous density dependence and capture pressure limit the abundance of snapping turtles. Optimal effort allocation effectively improved the control effect, but the degree of improvement varied substantially depending on the total effort. The degree of improvement by the spatial optimization of allocation effort was only 3.21% when the total effort was maintained at the 2016 level. However, when the total effort was increased by two, four, and eight times, spatial optimization resulted in improvements of 4.65%, 8.33%, and 20.35%, respectively. To achieve the management goal for snapping turtles in our study area, increasing the current total effort by more than four times was necessary, in addition to optimizing the spatial effort. The snapping turtle population is expected to reach the target density one year after the optimal management strategy is implemented, and this rapid response can be explained by high population growth rate coupled with density-dependent feedback regulation. Our results demonstrated that combining a state-space model with optimization makes it possible to adaptively improve the management of invasive species and decision-making. The method used in this study, based on removal records from an invasive management program, can be easily applied to monitoring data for wildlife and pest control management using traps in a variety of ecosystems.

The necessity of tailored control of irrupting pest populations driven by pulsed resources

Resource pulses are widespread phenomena in diverse ecosystems. Irruptions of generalist consumers and corresponding generalist predators often follow such resource pulses. This can have severe implications on the ecosystem and also on the spread of diseases or on regional famines. Suitable management strategies are necessary to deal with these systems. In this study, we develop a general model to investigate optimal control for such a system and apply this to a case study from New Zealand. In particular, we consider the dynamics of beech masting (episodic synchronous seed production) leading to rodent outbreaks and subsequent stoat (Mustela erminea) irruptions. Here, stoat control happens via secondary poisoning. The results show that the main driver of the optimal control timing (June) is the population density of the control vector. Intermediate control levels are superior to higher levels if the generalist consumer is necessary as a control vector. Finally, we extend the model to a two-patch metapopulation model, which indicates that, as a consequence of the strong vector dependence, a strategy of alternating control patches yields better results than static control. This highlights that besides control level, also the design impacts the control success. The results presented in this study reveal important insights for proper pest management in the New Zealand case study. However, they also generally indicate the necessity of tailored control in such systems.

Long-term impact of agricultural practices on the diversity of small mammal communities: a case study based on owl pellets

Small mammals have been seldom used as indicators of biodiversity responses to environmental changes, probably because their long-term population trend in a given area is not easy to monitor. To assess the impact of agricultural intensification in a protected area of northern Italy, we compared the composition of its small mammal communities, as assessed in 1994-1995 and 2015-2016 by the analysis of owl pellets (N = 265 and 302, respectively), which provides an effective and affordable method for assessing changes in the diversity and structure of small mammal assemblages over time. We recorded a sharp reduction in the frequency of occurrence of shrews (Sorex spp. and Crocidura spp.), which were replaced by generalist/anthropophilic rats (Rattus norvegicus) and house mice (Mus domesticus). Overall richness and diversity of the community varied only slightly, while trophic level and functional diversity indices clearly reflected the decline of the predator-level fraction of the community. We could reliably exclude both broad-scale land use- and climate changes as drivers of variation in the composition of small mammal communities and ascribe the decline of insectivores to changes in agricultural practices, namely the increase in cover of maize fields and spread of both herbicides and insecticides. Our results are consistent with the general opinion that crop specialization and increasing chemical inputs reduce the diversity and abundance of invertebrate prey, with bottom-up effects on higher trophic levels.

Strengths and potential pitfalls of hay transfer for ecological restoration revealed by RAD-seq analysis in floodplain Arabis species

Achieving high intraspecific genetic diversity is a critical goal in ecological restoration as it increases the adaptive potential and long-term resilience of populations. Thus, we investigated genetic diversity within and between pristine sites in a fossil floodplain and compared it to sites restored by hay transfer between 1997 and 2014. RAD-seq genotyping revealed that the stenoecious floodplain species Arabis nemorensis is co-occurring with individuals that, based on ploidy, ITS-sequencing and morphology, probably belong to the close relative Arabis sagittata, which has a documented preference for dry calcareous grasslands but has not been reported in floodplain meadows. We show that hay transfer maintains genetic diversity for both species. Additionally, in A. sagittata, transfer from multiple genetically isolated pristine sites resulted in restored sites with increased diversity and admixed local genotypes. In A. nemorensis, transfer did not create novel admixture dynamics because genetic diversity between pristine sites was less differentiated. Thus, the effects of hay transfer on genetic diversity also depend on the genetic make-up of the donor communities of each species, especially when local material is mixed. Our results demonstrate the efficiency of hay transfer for habitat restoration and emphasize the importance of prerestoration characterization of microgeographic patterns of intraspecific diversity of the community to guarantee that restoration practices reach their goal, that is maximize the adaptive potential of the entire restored plant community. Overlooking these patterns may alter the balance between species in the community. Additionally, our comparison of summary statistics obtained from de novo- and reference-based RAD-seq pipelines shows that the genomic impact of restoration can be reliably monitored in species lacking prior genomic knowledge.

Resilience of antagonistic networks with regard to the effects of initial failures and degree-degree correlations

In this study we investigate the resilience of duplex networked layers α and β coupled with antagonistic interlinks, each layer of which inhibits its counterpart at the microscopic level, changing the following factors: whether the influence of the initial failures in α remains [quenched (case Q)] or not [free (case F)]; the effect of intralayer degree-degree correlations in each layer and interlayer degree-degree correlations; and the type of the initial failures, such as random failures or targeted attacks (TAs). We illustrate that the percolation processes repeat in both cases Q and F, although only in case F are nodes that initially failed reactivated. To analytically evaluate the resilience of each layer, we develop a methodology based on the cavity method for deriving the size of a giant component (GC). Strong hysteresis, which is ignored in the standard cavity analysis, is observed in the repetition of the percolation processes particularly in case F. To handle this, we heuristically modify interlayer messages for macroscopic analysis, the utility of which is verified by numerical experiments. The percolation transition in each layer is continuous in both cases Q and F. We also analyze the influences of degree-degree correlations on the robustness of layer α, in particular for the case of TAs. The analysis indicates that the critical fraction of initial failures that makes the GC size in layer α vanish depends only on its intralayer degree-degree correlations. Although our model is defined in a somewhat abstract manner, it may have relevance to ecological systems that are composed of endangered species (layer α) and invaders (layer β), the former of which are damaged by the latter whereas the latter are exterminated in the areas where the former are active.

Feral cat diet and predation on endangered endemic mammals on a biodiversity hot spot (Amami–Ohshima Island, Japan)

Context There has been concern that feral cats have negative impacts on the endangered endemic mammals of Amami-Ohshima Island, Japan, including the Amami rabbit, Pentalagus furnessi, Ryukyu long-tailed giant rat, Diplothrix legata, and Amami spiny rat, Tokudaia osimensis. However, no diet study of feral cat has been conducted to support the necessity of an urgent feasible feral-cat management for the island. Aims The aims of the present study were to analyse feral-cat diet on Amami-Ohshima Island by using scat analysis and estimate the potential predation impact of feral cats on endangered mammals on the island. Methods The diet of feral cats was studied using scat analysis. We estimated the number of prey, percentage of prey, frequency of occurrence (the percentage of scats in a sample containing a particular prey item), percentage of biomass (biomass of the same prey item divided by the total consumed biomass ×100) and daily consumed biomass (DCB). Key results Three endangered endemic mammals were the main prey species of the feral cat diet (65% of total DCB). The percentage contributions of these species on DCB were long-tailed giant rat (34.7%), Amami spiny rat (21.9%) and Amami rabbit (12%). Conclusions Mammals, especially endangered endemic mammals, were main prey species of feral cat on Amami Island. In Amami Island, where native and invasive rodents coexisted, feral cats consumed more native (56.6%) than invasive (22.2% for Rattus rattus) species. Implications Feral cats are likely to be having a significant impact on endangered endemic mammals on the island. To ensure the long-term survival of these endemic species, active management of the feral-cat population should be considered.

Evaluating the "recovery level" of endangered species without prior information before alien invasion

For maintaining social and financial support for eradication programs of invasive species, quantitative assessment of recovery of native species or ecosystems is important because it provides a measurable parameter of success. However, setting a concrete goal for recovery is often difficult owing to lack of information prior to the introduction of invaders. Here, we present a novel approach to evaluate the achievement level of invasive predator management based on the carrying capacity of endangered species estimated using long-term monitoring data. In Amami-Oshima Island, Japan, where the eradication project of introduced small Indian mongoose is ongoing since 2000, we surveyed the population densities of four endangered species threatened by the mongoose (Amami rabbit, the Otton frog, Amami tip-nosed frog, and Amami Ishikawa's frog) at four time points ranging from 2003 to 2011. We estimated the carrying capacities of these species using the logistic growth model combined with the effects of mongoose predation and environmental heterogeneity. All species showed clear tendencies toward increasing their density in line with decreased mongoose density, and they exhibited density-dependent population growth. The estimated carrying capacities of three endangered species had small confidence intervals enough to measure recovery levels by the mongoose management. The population density of each endangered species has recovered to the level of the carrying capacity at about 20–40% of all sites, whereas no individuals were observed at more than 25% of all sites. We propose that the present approach involving appropriate monitoring data of native organism populations will be widely applicable to various eradication projects and provide unambiguous goals for management of invasive species.