www.nature.com/scientificreports

OPEN on endangered by -subsidized domestic on Tokunoshima Island Tamao Maeda1, Rumiko Nakashita2, Kazumi Shionosaki3,4, Fumio Yamada2 & Yuya Watari 2*

It is important to unravel how impact native in order to control them efectively. The presence of abundant exotic prey promotes of invasive predators, thereby enhancing the predation pressure on native prey (hyper-predation). Not only the exotic prey but also feeding by is likely to cause “hyper-predation”. However, the contribution of artifcial resources to this was underestimated in previous studies. Here, we combined fecal and stable isotope analyses to reveal short- and long-term food habits of free-ranging cats on Tokunoshima Island. Although 20.1% of the feces contained evidence of forest-living species, stable isotope analysis suggested that the cats were mostly dependent on artifcial resources. In addition, a general linear model analysis showed that their diet was strongly correlated with landscape variables. These results indicate that the invasive free-ranging cats are aided by anthropogenic feeding, and they move from the human habituated area to natural areas with high . These fndings suggest the possibility of human feeding indirectly accelerates the efect of cat predation, and call for a further study on their demography. Cat management mainly involves trapping, but our fndings show that educating local residents to stop feeding free-ranging cats and keeping cats indoors are also important.

Biological invasion is one of the major causes of biodiversity loss globally1,2, especially in insular ecosystems worldwide3,4. To efectively control invasive species, it is important to determine how they survive and increase their populations in the areas that they invade. Te establishment of populations is one of the essential processes for successful invasion5,6; in other words, the impact of the introduced species may not become evident unless they have succeeded in increasing their numbers7,8. Various factors, such as species traits, biotic and abiotic envi- ronments, and temporal events, can afect a result of such establishment, subsequent increases in , and the ecological impact of an invasive species8,9. Hyper-predation is one of the possible processes enhancing the impact of invasive predators10. Tis hypothesis predicts that the presence of abundant primary prey subsidizes the predator population, allowing it to grow and then more severely impact the relatively scarce native prey11,12. Feral cats (Felis catus) are among the most infuential introduced species13, as they have been responsible for the or decrease in numerous mammals, , and , particularly in insular ecosystems13–19. Many studies have reported that introduced prey (e.g., European , black , and house mice) are sus- pected of causing hyper-predation to feral cats on native organisms11,18,20. Te most well-known example is the case on Macquarie Island, where invasive feral cats caused extinction of the endemic parakeet, Cyanoramphus novaezelandiae erythrotis7. Feral cats and parakeets coexisted on the island for 60 years, but afer rabbits were introduced, feral cats rapidly increased their number and ate up the parakeets in the next 10 years. Not only the introduced prey but also direct or indirect feeding by humans support cat populations and enhance predation pressure on , thereby accelerating their extinction10,21–23. Some previous studies reported that cats fed by people can substantially impact on the local ecosystem21,24. However, most of the studies were limited to native predators living in urban or peri-urban areas23. Tokunoshima Island is located in southwestern Japan and is a biodiversity hotspot with unique biota that evolved in the absence of native mammalian predators25. On this island and adjacent Amami-Oshima and Okinawa Island, free-ranging cats prey on endangered endemic species, such as Amami (Pentalagus fur- nessi), Ryukyu long-haired (Diplothrix legata), and spiny rat (Tokudaia tokunoshimensis, T. osimensis and

1Wildlife Research Center, Kyoto University, 2-24 Tanaka-Sekiden-cho, Sakyo, Kyoto, 606-8203, Japan. 2Forestry and Forest Products Research Institute (FFPRI), 1 Matsunosato, Tsukuba, Ibaraki, 305-8687, Japan. 3Amami Wild Research Center, 2662 Ogachi, Tatsugo-cho, Kagoshima, 894-0105, Japan. 4Amami Research Center Co., Ltd, 10-11-2F Naze Suehiro-cho, Amami, Kagoshima, 894-0027, Japan. *email: [email protected]

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 1 www.nature.com/scientificreports/ www.nature.com/scientificreports

Figure 1. Feral cats killing endemic mammals taken by censor cameras on Tokunoshima Island. (a) Amami rabbit, (b) Ryukyu long-haired rat. (a) Photograph was taken in 2017 and provided with permission by the Naha Nature Conservation Ofce, Ministry of the Environment. (b) Photograph taken in 2018 by the author (Y. Watari, Forestry and Forest Products Research Institute).

T. muenninki)26,27, and are believed to be responsible for the population reduction of these species on the islands28,29 (Fig. 1). Indeed, a camera-trapping study by the Japanese Ministry of the Environment found a neg- ative correlation between the number of cat appearances and that of endangered Amami rabbits, indicating that the presence of cats can limit the rabbit distribution30. Te local Tokunoshima Island government and the Ministry of the Environment have been capturing free-ranging cats since 2014 to conserve the endemic species. Cats caught in forests are referred to here as “feral” cats, whereas those caught in residential areas or farmlands are called “stray” cats. “Feral” cats are kept in a shelter, and some can be adopted by new owners afer sterilization, whereas “stray” cats are returned to the wild afer sterilization (TNR) because it is implicitly assumed that they do not enter the forest and prey on native . Tis management program has been achieving some degree of success, as the monthly route census conducted by the Ministry of the Environment showed that the encounter rates of three endangered mammals, Amami rabbits, Tokunoshima spiny rats (T. tokunoshimensis), and Ryukyu long-haired rats, have been increasing while that of cats are decreasing since 201431. It has also been suggested that free-ranging cats are placing substantial predation pressure on native species. Scientifcally, the term “feral” means completely independent and rarely interacting with humans, whereas “stray” cats do not have an owner but still depend on human care32. Te division of “feral” and “stray” cats by the government suggests their assumption that cats rarely migrate between forests and residential areas, but it is still unclear whether this division is scientifcally appropriate due to the lack of studies on free-ranging cats on Tokunoshima Island. Tokunoshima Island is characterized by small forested areas33, so it is rather likely that “feral” cats and “stray” cats have access to both wild animals in the forest and artifcial food in the villages. Even the core area of the forest is only a few kilometers away from farmlands, which is close enough for free-ranging cats to access both34–37. Human garbage was found in 7.1–50% of cat feces on the northern part of Okinawa Island (Yambaru) and Amami-Oshima Islands26,27. However, fecal analysis cannot be used for accurate estimation of the actual dependence on artifcial resources, owing to the diferent digestibility among food items38. Instead, stable isotope analysis is a powerful method to clarify the dependence of subjects on a given food item39. If the so-called “feral” cats are fed with human food, in order to reduce the high predation on native species, it could be efective to keep domestic cats indoors or stop feeding cats without owners. Evidence of dependence of cats would provide strong support in promoting public awareness of this issue and in developing an efective strategy for conserving native species. Tis study evaluated the diet of free-ranging cats on Tokunoshima Island to verify our hypothesis that both “feral” and “stray” cats are accessible to the forest and the residential area and they are highly dependent on the human-provided resources, although they also predate endemic species in the forest. Tis dietary state of cats is one of the necessary conditions for the occurrence of hyper-predation induced by human-derived food resources. We specifcally addressed three questions: (1) Do “feral” cats eat more ofen than “stray” cats? (2) Does providing artifcial food substantially support diet of free-ranging cats? (3) Do the contributions of forest prey (including endangered species) and artifcial food of free-ranging cats difer among capture locations with diferent surrounding landscapes? To answer these questions, we collected fecal and hair samples from trapped free-ranging cats and conducted fecal and stable isotope analyses. Results Fecal analysis. In total, 208 “feral” cats (75 females, 123 males, and 10 unidentifed) and 54 “stray” cats (22 females, 30 males, and 2 unidentifed) were captured, and 198 fecal samples (from 174 “feral” cats and 24 “stray” cats) were obtained (Fig. 2). A total of 13.4% of the cats (35; 31 “feral” and 4 “stray”) were ear-tipped, which means that they had been captured as “stray” cats and sterilized (Table S1). Evidence of forest-living species and farm- land-living animals were found in at least 17.7% and 30.8% of the fecal samples, respectively (Table 1). Overall, 23.7% of fecal samples contained artifcial objects, such as plastic or paper (Table 1). In addition, 20.1% of the fecal samples from “feral” cats contained evidence of forest animals, which was signifcantly higher than the rate for “stray” cats (0.0%) (Fisher’s exact test, p < 0.01). No signifcant diference was observed in the occurrence frequency of farmland species (“feral”: 31.6%, “stray”: 20.8%) and artifcial objects (“feral”: 24.1%, “stray”: 20.8%) (p > 0.05). Six threatened species (at least 43 individuals) were detected in 13.5% of the fecal samples: Ryukyu robins (Erithacus komadori komadori), which are endemic to southern Japan40, and Amami rabbit, Ryukyu

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 2 www.nature.com/scientificreports/ www.nature.com/scientificreports

Figure 2. Map of Tokunoshima Island, which is located on the Amami Islands in the Rukyu Chain and the capture location of cats.

long-haired rat, Tokunoshima spiny rat, Crocidura spp., and Amami tip-nosed (Odorrana amamiensis), which are endemic to the Ryukyu Islands41. Black rats and were the only non-native prey in this study42. Te average weight of captured cats was 3.3 ± 1.0 kg (range: 1.0–6.0 kg; male: 3.6 ± 0.9 kg, female: 2.7 ± 0.7 kg). Tus, the estimated average daily consumed (DCB) of cats was 379 ± 143 g (range: 146–629 g). Te body mass of Amami rabbits and chickens exceeded the maximum DCB of the captured cats, so for these species, we used the maximum DCB (629 g) as the weight. Te results showed that this method explained only 24.2% of the cats’ diet (forest animals: 15.5%, farmland animals: 8.7%). Te contribution from forest animals was mostly based on two endangered mammals, Ryukyu long-haired rats (7.7%) and Amami rabbits (6.7%), and that of farmland animals on black rats (6.9%).

Isotopic mixing model. We analyzed the hair of 189 “feral” cats, 52 “stray” cats, and 9 indoor cats. Te sta- ble isotope ratios of carbon were −17.4 ± 1.4‰, −17.2 ± 1.2‰, and −16.9 ± 1.7‰, and those of nitrogen were 7.0 ± 0.9‰, 7.1 ± 0.8‰, and 6.8 ± 0.8‰, respectively (Figs 3 and 4). Analysis of variance (ANOVA) detected no signifcant diference in the stable isotope ratio among “feral”, “stray”, and indoor cats [δ13C: F(2,247) = 1.15, p = 0.319; δ15N: F(2,247) = 0.43, p = 0.651]. We obtained diferent brands of dried cat food (n = 9) and hair samples of Amami rabbit (n = 7), Ryukyu rat (n = 7), and black rats (n = 7) as candidate cat dietary resources. Te stable isotope ratios of carbon in for- est animals, farmland animals, and artifcial resources were −24.8 ± 2.6‰, −20.9 ± 2.1‰, and −18.9 ± 2.5‰, and those of nitrogen were 1.6 ± 1.3‰, 6.4 ± 1.3‰, and 4.6 ± 1.2‰, respectively. ANOVA on the isotope ratio revealed signifcant variation among resources, δ13C: F(2,33) = 16.43, p < 0.001; δ15N: F(2,33) = 43.84, p < 0.001. Post hoc Tukey’s test showed that forest animals had signifcantly lower δ13C than the rest (p < 0.01). Farmland animals had the highest δ15N, artifcial resources had the second highest, and forest animals had the lowest (p < 0.01). Te changes in δ13C and δ15N of sheltered cats versus time are shown in Fig. S1. Te estimated asymptotes of the regression model (TEF) of δ13C and δ15N were 2.3 ± 0.3 and 2.8 ± 0.1, respectively. Te stable isotope analysis in R (SIAR) indicated that artifcial resources were the largest component in energy consumption of the “feral” (67.8%; 95% highest density region: 62.8–72.8%) and “stray” cats (69.0%; 59.3–78.8%), followed by farmland animals (“feral”: 17.9%; 13.4–22.3%, “stray”: 18.5%; 9.7–27.3%) and forest animals (“feral”: 14.3%; 11.6–17.1%, “stray”: 12.4%; 9.7–27.3%) (Fig. 5). In addition, this high dependence on artifcial resources remained when we conducted the SIAR only on cats that had evidence of forest animals in their feces (Fig. S3), even though they tended to be captured closer to the forest than the undetected individuals (Fig. S2, Appendix). Te estimated dependence on wild animals is generally consistent with the contribution calculated from the fecal analysis (forest animals: 15.5%, farmland animals: 8.7%).

Efect of landscape elements on the cat diet. In factor analysis, residential area coverage and building density both loaded positively on factor 1, whereas forest coverage and farmland coverage loaded negatively and positively on factor 2, respectively (Table S2). Te results of the general linear model (GLM) revealed that dependence on artifcial resources was positively correlated with residential area coverage (factor 1) and weight, dependence on farmland animals was positively correlated with farmland coverage (factor 2), and dependence on residential area was positively correlated with forest coverage (factor 2) but negatively correlated with body weight (Table 2).

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 3 www.nature.com/scientificreports/ www.nature.com/scientificreports

FO (%) NI body TOTAL feral stray TOTAL feral stray Contribution weight IUCN Red List Reference for Items (n = 198) (n = 174) (n = 24) (n = 198) (n = 174) (n = 24) to DCB (%) (g) 2017 body weight Forest species 17.7 20.1 0.0 40 40 0 15.5 Diplothrix legata 6.1 6.9 0.0 12 12 0 7.7 483 EN 27 Pentalagus furnessi 4.0 4.6 0.0 8 8 0 6.7 2880* EN 75 Tokudaia tokunoshimensis 3.0 3.4 0.0 6 6 0 1.3 162.4 EN Yamada (unpublished) Erithacus komadori komadori 0.5 0.6 0.0 1 1 0 0.0 22.4 VU 55 Turdus pallidus 0.5 0.6 0.0 1 1 0 0.1 78 27 Odorrana amamiensis 0.5 0.6 0.0 1 1 0 0.1 60 VU 76 Diestrammena gigas 5.1 5.7 0.0 10 10 0 0.0 3 76 Tereuopoda clunifera 0.5 0.6 0.0 1 1 0 0.0 3.5 27 Farmland species 30.8 31.6 20.8 70 65 5 8.7 Rattus rattus 24.2 26.4 20.8 53 48 5 6.9 98 ** (C. orii) EN Crocidura spp. 6.6 7.5 0.0 15 15 0 0.1 7 77 (C. watase) NT Gallus gallus domesticus 1.0 1.1 0.0 2 2 0 1.7 1500-* 78 Unknown — — — Horornis diphone 1.0 1.1 0.0 2 2 0 0.0 15.8 79 Unidentifed birds 15.7 17.8 0.0 31 31 0 — — — /Reptiles 3.5 4.0 0.0 7 7 0 — — — Orthoptera 8.6 6.3 25.0 17 11 6 — — — Mantodea 1.0 1.1 0.0 2 2 0 — — — Coleoptera 3.0 2.9 4.2 6 5 1 — — — Unidentifed insects 41.4 46.6 4.2 82 81 1 — — — Crustacea 1.5 1.7 0.0 3 3 0 — — — Gastropods 0.5 0.0 4.2 1 0 1 — — — Artifcial objects 23.7 24.1 20.8 47 42 5 — — — Plants 42.4 32.2 50.0 84 72 12 — — —

Table 1. Result of fecal analysis. FO, frequency of occurrence of items and NI, number of prey items. *As these species (Pentalagus furnessi and Gallus domesticus) were heavier than the max DCB, we used max DCB to calculate the contribution to DCB. **We calculated the mean body weight of the captured Rattus rattus (n = 12).

Figure 3. Stable isotope ratios of (a) carbon and (b) nitorogen in feral, stray, and indoor cats. Error bars represent standard deviations and the bold bars in the box represent medians.

Discussion Our study revealed the dietary habits of invasive cats on Tokunoshima Island by combining fecal analysis with stable isotope analysis. A dietary diference between “feral” and “stray” cats was only detected in the fecal analysis, which refects the diet of the past few days in each where the cats were captured. According to the isotopic mixing model, the cats’ long-term diet did not difer signifcantly, and even the cats that had evidence of forest ani- mals in their feces largely depended on artifcial resources. It is likely that the cats visited the forests for a few days, where they hunted native endangered animals, and then traveled back to the villages to eat cat food, which was

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 4 www.nature.com/scientificreports/ www.nature.com/scientificreports

Figure 4. Stable isotope ratios of the cats and their potential resources. Error bars represent standard deviations. Farmland animals and artifcial resources are represented by black rats and pet food, respectively. Te stable isotope ratio of the forest animals was the average of Amami rabbits and Ryukyu long-haired rats.

Figure 5. Dependency of captured cats on three resource types, including artifcial resources, farmland animals, and forest animals. Error bars represent the 95% high density region. Values were derived from stable isotope analysis in R (SIAR)80.

Artifcial Resources Farmland Animals Forest Animals Estimate Estimate Estimate Explanatory variables (×10−3) z p (×10−3) z p (×10−3) z p Residential Area (Factor1) 10.65 2.59 ** −0.11 0.06 −5.83 0.83 +Farmland/−Forest (Factor2) −0.10 0.08 16.08 2.84 ** −2.09 2.85 ** weight 16.62 3.54 *** 0.72 0.27 −22.15 3.33 *** sex (male) −4.43 0.54 1.14 0.23 0.56 0.13 MEM1 −9.43 2.24 * −11.83 2.29 * 22.82 3.53 *** MEM2 −0.38 0.21 11.16 2.13 * −6.36 0.89 MEM3 0.00 0.00 −10.79 2.09 * 5.85 0.84 MEM4 −0.03 0.03 0.02 0.01 −0.03 0.02 MEM5 −1.31 0.45 −16.08 3.24 ** 24.45 3.93 *** MEM6 −0.04 0.04 14.00 2.62 ** −14.62 2.18 *

Table 2. Te results of model averaging of the selected general linear model (GLM) with corrected version of Akaike information criterion Δ (AIC) < 2 (family = Gaussian). Te heading row represents the response variables and dependency on each of the three resources; *p < 0.05, **p < 0.01, ***p < 0.001.

their main source of food. In addition, many of the captured “feral” cats were ear-tipped, i.e. they had been captured as “stray” cats in the residential area, and this ratio was not signifcantly diferent from that of “stray” cats. It also suggests movements of cats between the forest and the villages. Te endemic mammal population has been substan- tially impacted by cats28,31, whereas our study shows that cats themselves depend on human-derived resources. In addition, both stable isotope analysis and fecal analysis suggested relatively low dependence on farmland animals, namely, black rats, unlike the case on many other islands where introduced prey is available16. Tis may be due to the low density of black rats in the forests of Tokunoshima Island30. In fact, a rat population survey by Jogahara showed that spiny rats dominated in the forests, whereas black rats were rarely captured (unpublished data). Free-ranging unowned cats are usually divided into two categories: feral cats, which depend on native resources, and stray cats, which depend on artifcial resources32. Te GLM results show that dependence on each resource had a positive association with the land use where it was assumed to have been obtained. Te

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 5 www.nature.com/scientificreports/ www.nature.com/scientificreports

dependence on forest-living animals increased at capture locations close to the forest, suggesting that feral and stray cats cannot be clearly separated and that free-ranging cats can be a threat to native species, particularly in human residential areas adjacent to the . If the forest was large enough, their dependence on artifcial resources would be minimized, which would mean that pure feral cats would breed. In other words, small , such as the forest on Tokunoshima Island, may be more susceptible to the efect of human-derived resource subsidization as well as other kinds of efects27. In addition, roads may also increase the accessibility of forests from residential areas, as ofen prefer to move on tracks or roads43,44, and people can use them to abandon their in natural areas more easily. According to the local Tokunoshima Island government, 2,797 “stray” cats were captured and sterilized from April 2014 to March 2018. However, only 13% of the captured cats were ear-tipped, and this proportion is not increasing. The results imply the huge number of cats on the island and their successful reproduc- tion. We assume that stable and inexhaustible human-derived resources enable cats to sustain this large population, but further investigation would be needed to evaluate the effect of artificial resources on cat demographics. Although the management may have been successful in achieving some recovery of endemic mammals, it might be difficult to eliminate free-ranging cats unless the resource subsidization by humans is controlled. Overall, our study indicates that invasive free-ranging cats depend on anthropogenic feeding, the efect of which may reach far from the habituated area to natural areas with high biodiversity. Tis fnding provides new insight into how to best manage invasive cat populations. Te main predator management options are trapping, which has occurred on Tokunoshima Island, and lethal control10,45. However, owing to the necessity for contin- uous intervention, such management is ofen very expensive, which sometime leads to failure of the whole pro- ject45. In the case of human-driven hyper-predation, preventing the access of cats to artifcial resources is a more cost-efective way of reducing the predator population in the long term. Studies on this topic will be important to better plan predator control programs. If the local people provide additional resources to predators without being aware of their impact on the , introducing scientifc evidence for human-driven hyper-predation may improve their awareness of how to treat their pets and neighboring wildlife appropriately. Although this method would efectively reduce the predator population in the long term, a sudden decrease of resource subsidization usually causes a temporary increase in predation pressure on the native prey29,46. Several methods should thus be combined, including lethal control and resource subsidization control, to develop an efective conservation strategy. Tokunoshima Island has regulations about keeping pet cats indoors and prohibits the feeding of unowned cats. However, as the mixing model showed high dependence on artifcial resources, it is likely that many people are not following these regulations. Our study provides important scientifc evidence to prove the need to edu- cate people and support such regulations. In addition, this study suggests the potential impact of TNR cats on endangered species in the short term. Te efectiveness and validity of this method should thus be reconsidered. Our study also points out the limitations of fecal analysis in detecting the efect of anthropogenic subsidization because only 24% of the feces contained artifcial objects, despite the high dependence suggested by the isotopic mixing model. Previous studies detected artifcial materials in cat feces or stomach/gut contents, but they were usually excluded from the subsequent dependence calculation due to the low frequency of occurrence of artifcial items and the difculty in estimating their mass27,47,48. It is possible that these studies largely underestimated the efect of feeding by humans; thus, it would be better to combine multiple methods, such as stable isotope analysis, to obtain accurate estimates of the diet49. Although our study suggested the occurrence of resource subsidization by humans on the cat population, we still lack the demographic and ethological studies on free-ranging cats and endemic species. Te high dependency of diet does not necessarily mean that the resource is essential for supporting the population because it might be just a consequences of resource selectivity. It is required to research whether human feeding actually gives posi- tive impacts on cat population, and predation by cats gives negative impacts on the endemic mammals in order to fully verify our assumption of hyper-predation by human feeding and to understand its precise process and consequences. For example, if cat population increases due to anthropogenic resource subsidization, the density of cats may positively correlate with that of residents. In conclusion, our study provides strong circumstantial evidence of anthropogenic resource subsidization on free-ranging cats. It points out the possibility of human-driven hyper-predation and provides important support for promoting local and global invasive predator control management. Methods Study area. Tokunoshima Island (N 27°45′, E 128°58′) is located in the Ryukyu Archipelago, southwestern Japan (Fig. 2). It has an area of 247.85 km2 and a population of ~25,000 inhabitants33. Tokunoshima Island is in a subtropical region with high precipitation (mean temperature, 21.6 °C; mean annual rainfall, 1912 mm). It has mountains (highest peak, 645 m) that run north to south, surrounded by a plateau of Ryukyu limestone. Tis island consists primarily of crop felds and broad-leaved evergreen forests, covering 28% and 43% of the area, respectively. Sugarcane predominates as a crop feld, and evergreen oak species, such as Castanopsis sieboldii and Quercus miyagii, are the major components of the forest. Te Ryukyu Archipelago, especially the Central Ryukyus, including Tokunoshima Island, had separated from the Eurasian continent at least by the late Miocene (11.63–5.33 million years ago)50. Te native top predators are habu vipers (Protobothrops favoviridis and okinavensis)51, and a great number of endemic species and subspecies have evolved in the absence of native mammalian predators40. Many of them are highly endemic, e.g., Amami rabbit is endemic to Tokunoshima and adjacent Amami-Oshima Island; Ryukyu long-haired rat is endemic to Tokunoshima, Amami-Oshima, and the northern part of Okinawa Island; and Tokunoshima spiny rat

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 6 www.nature.com/scientificreports/ www.nature.com/scientificreports

only lives on Tokunoshima Island41. Most of the endemic species, including these three mammals, are threatened and listed on the Red List of the International Union for Conservation of Nature and Natural Resources (IUCN 2017) and the Japanese Ministry of the Environment (2017).

Sample collection and dietary analysis. Cat diets are usually evaluated by stomach content or fecal analysis16,52. Tese methods provide a short-term picture of the diet but underestimate the contribution of highly digestible material, such as pet food, and immeasurable objects, such as garbage. Stable isotope analysis is an alternative way to identify major food items, as it provides information on the long-term contributions of major foods and is less afected by diferences in digestibility. However, as potential prey species may have similar iso- topic values, taxonomic resolution is occasionally low, particularly for generalist and opportunistic predators such as cats. A combination of fecal and stable isotope analyses can compensate for the other’s disadvantages and reveal a more precise, long-term dietary history53,54. In this study, hair and fecal samples were obtained from feral and stray cats captured in the population con- trol program conducted on Tokunoshima Island. Feral cats were trapped in forest areas more than 500 m away from villages, whereas stray cats were trapped in or around villages. Both feral and stray cats were captured using metal box traps with cat food or fried inside and then brought to an animal hospital for sterilization. Hair samples were collected by a vet during surgery. Afer sterilization, the cats were kept in separate cages for a few days for the collection of feces. samples were collected from December 2014 to January 2018, and stray cat samples were collected in November 2017. Te capture location, capture date, sex, and body weight were recorded for each cat. When ear-tipped cats, individuals which had experienced TNR in the past, were captured afer November 2017, we compared them with the photos of stray cats captured previously to make sure that they had not yet been sampled.

Fecal analysis. Fecal samples were kept in plastic bags frozen at −20 °C. Te feces were washed over a 1-mm mesh sieve under a stream of water and dried in an oven at 65 °C for more than 12 h. Each food item was identi- fed to the species level and assigned to one of the four main habitat types: forest-living species (forest animals), farmland- and residential area-living species (farmland animals), artifcial resources, and unidentifed animal/ plant materials. Most of the species exclusively live in either forest or non-forest areas30,41,55, except Horornis diphone; thus, we categorized this species as “unidentifed.” We categorized black rats as “farmland animals” because black rats rarely occurred in the forest where endangered species inhabited. Unidentifed animal/plant materials were excluded from the following dietary analysis. Te number of individual prey in each scat was counted based on distinctive bones, such as jaws and incisors. We estimated the frequency of occurrence and the minimum number of individuals for each prey species. To narrow down the candidate prey species for the stable isotope analysis, we estimated the contribution of each prey species to DCB of cats following the methods of Bonnaud et al.56 and Shionosaki et al.27. As cats usually defecate once per day57,58, the formula can be written as follows: mean body weight of prey × NnI/ Contribution = × 100(%) mean DCBofcats (1) where n is the total number of scat samples and NI is the minimum total of individual prey found in the scats. DCB of free-living, eutherian predators can be estimated using the allometric equation: DCB = 3.358 × (body weight of predator)0.813 × 2.86/18 (g)59,60. Here, 2.86 is included to account for the 65% water content of prey (100/(100 − 65) = 2.86), and 18 represents the mean energy content in kJ of metabolizable energy per gram of dry prey59,60. We set the upper limitation of body weight of prey as the maximum DCB of the cats because, when cats catch large prey, they are likely to eat some and leave the rest61. We previously ran SIAR with prey species whose con- tribution was >1% and obtained the result which reveals that the smallest 95% HDR interval was 3.0%. Tus, we defned important prey species for cats as those whose contribution was >3%. Te results showed that two forest-living species (Amami rabbits and Ryukyu long-haired rats) and one farmland- and residential area-living species (black rats) satisfed this threshold and were used for the following stable isotope analysis. We compared the frequency of occurrence of each prey category between feral cat feces and stray cat feces using Fisher’s exact test.

Stable isotope analysis. Te entire hair of cats and prey species, including the root, was plucked and kept in a plastic bag. Hair samples of indoor pet cats, which had been fed only pet food, were also collected for com- parison with the feral and stray cats. In addition, hairs of “sheltered cats” (cats captured as feral and thereafer kept in a shelter for 23–536 days and supplied with pet food) were taken to estimate the trophic enrichment factor (TEF). All hair samples were provided by the cat population control programs by the local government on Tokunoshima Island and the Japanese Ministry of the Environment, which were carried out in accordance with the Act on Welfare and Management of Animals and Protection and Control of Wild Birds and Mammals and Management Law, respectively. We assumed three types of dietary resources for the cats according to the results of the fecal analysis: for- est animals (Amami rabbit and Ryukyu rats), farmland animals (black rats), and artifcial resources (pet food). Black rats were captured using metal box traps in villages and farmlands, and hairs were plucked from their necks. Samples of endangered Amami rabbits and Ryukyu rats were obtained, with permission, from frozen car- casses (mostly killed in trafc accidents) stored by the Ministry of the Environment. As a representative artifcial resource, pet food was analyzed because it is the major food people feed to cats on Tokunoshima Island, and it is likely to have an isotope ratio similar to that of other possible artifcial resources, such as lefover meals and garbage, as the ingredients of pet food resemble the human diet, namely, grains, fsh, meat, and soy62.

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 7 www.nature.com/scientificreports/ www.nature.com/scientificreports

We analyzed the carbon and nitrogen stable isotope ratios in the samples using a method similar to that of Mizukami et al. (2005a, 2005b)63,64. Te hair was rinsed with a 2:1 chloroform–methanol solution to remove lipids and was air-dried. It is recommended that lipids be removed because they are depleted in 13C relative to carbo- hydrates and proteins65, and their amount can vary greatly among individuals66. Pet food was dried in an oven at 65 °C for >12 h and pulverized with a food mill. Samples were enclosed in a tin cup and combusted in a FlashEA 1112 elemental analyzer (Termo Fisher Scientifc, Bremen, Germany) interfaced to a Delta V isotope ratio mass spectrometer (Termo Fisher Scientifc). Te analytical errors for the isotope analysis were within 0.1‰ for δ13C and 0.2‰ for δ15N.

Isotopic mixing model. Te Bayesian mixing model SIAR was applied using the R package “siar” to estimate the cats’ dependence on each resource67. Te SIAR model is suitable for Markov chain Monte Carlo (MCMC) methods in fnding a plausible dietary composition using Dirichlet prior distribution67. It is usually assumed that δ13C increases by 0‰–1‰ and δ15N by 3.4‰ from one to the next65,68. However, TEFs can difer among environments, trophic levels, tissues, species, and sample treatment proce- dures39,53,54,69,70. To estimate the appropriate TEF for this study, we analyzed the isotope ratio of sheltered cats. As sheltered cats were fed the same pet food afer being captured, their isotope ratio would converge with the isotope ratio of pet food +TEF. We used the asymptotic exponential model y = AeBx + C with ΔδX (δX of a shelter cat – δX of the pet food; X = 13C or 15N) as a response variable and time in days that a cat spent in the shelter as an explanatory variable. We defned TEF as the estimated asymptote (parameter C). We calculated the isotope ratio of forest animals (Amami rabbits and Ryukyu long-haired rats), farmland ani- mals (black rats), and artifcial resources (pet food). Te isotope ratio of forest animals was defned as the mean isotope ratio of the two species. Te MCMC was run 50,000 times, discarding the frst 5,000 samples and thinning by 10 {i.e., [number of groups × (number of sources + number of isotopes)] = 2 × (3 + 2)} to avoid sample autocorrelation.

General linear model. We used a GLM with a Gaussian structure (link = “identity”) under the R environ- ment to analyze the efect of landscape elements on the diet of free-ranging cats. Te response variable was the dependence of an individual cat on each of the three resources estimated by stable isotope analysis, which was taken as the arcsine square root transformation. Te explanatory variables were land-use variables surrounding capture locations, sex, body weight, and spatial autocorrelations. Land-use variables included forest, residential area, and farmland coverage and density of buildings. Te land coverage data were obtained from the National Land Numerical Information download service, and building locations were taken from the Geospatial Information Authority of Japan website. Cats fed by humans usually do not travel more than 300–800 m from the feeder’s house71,72. In addition, feral cats in our study were defned as those captured at least 500 m away from villages. Tus, we created radius bufers of 100, 200, and 500 m from the capture location. As land-use variables were strongly correlated with each other, we summarized them using an explanatory factor analysis. An exploratory factor analysis was conducted using maximum likelihood factor extraction to determine the factor structure of the landscape elements around the 165 capture locations of the 237 cats. Te reason for using factor analysis rather than a principal component analysis is that the axes of a factor analysis are easier to interpret in terms of land-use patterns. A parallel analysis recommended a two-factor solu- tion. We employed Promax (oblique) rotation to interpret the two factors. Spatial autocorrelation variables were added to consider the efect of spatial proximity on the cat diets. We constructed Moran’s eigenvector maps (MEM) using the Delaunay triangulation method and calculated the scores for each capture location using the R package “adespatial”73. As the larger MEM values representing a fner spatial structure may overlap with land use within the 100–500 m radius bufer, we considered MEM1-10 frst and then selected the model. Te largest signifcant MEM value was MEM6, so we used only MEM1–6 in subsequent analyses. Te model was selected using a multi-model inference approach. We used “MuMIn” package74 to produce all subsets of models based on the global model and ranked them based on the corrected version of Akaike informa- tion criterion (AIC). We used model averaging to produce the averaged parameter estimates of all models with ΔAIC < 2.

Received: 22 October 2018; Accepted: 14 October 2019; Published: xx xx xxxx

References 1. Bellard, C., Cassey, P. & Blackburn, T. M. Alien species as a driver of recent . Biol. Lett. 12 (2016). 2. Clavero, M. & García-Berthou, E. Invasive species are a leading cause of animal extinctions. Trends in and Evolution 20, 110 (2005). 3. Roemer, G. W., Donlan, C. J. & Courchamp, F. Golden eagles, feral , and insular carnivores: How exotic species turn native predators into prey. Proc. Natl. Acad. Sci. 99, 791–796 (2002). 4. Blackburn, T. M., Cassey, P., Duncan, R. P., Evans, K. L. & Gaston, K. J. Avian extinction and mammalian introductions on oceanic islands. Science (80-). 305, 1955–1958 (2004). 5. Sakai, A. K. et al. Te population biology of invasive species. Annu. Rev. Ecol. Syst. 32, 305–332 (2001). 6. Blackburn, T. M. et al. A proposed unifed framework for biological invasions. Trends in Ecology and Evolution 26, 333–339 (2011). 7. Taylor, H. R. How the Macquarie Island Parakeet became extinct. N. Z. J. Ecol. 2, 42–45 (1979). 8. Zenni, R. D. & Nuñez, M. A. Te elephant in the room: Te role of failed invasions in understanding invasion biology. Oikos 122, 801–815 (2013). 9. Grarock, K., Lindenmayer, D. B., Wood, J. T. & Tidemann, C. R. Using invasion process theory to enhance the understanding and management of introduced species: A case study reconstructing the invasion sequence of the common myna (Acridotheres tristis). J. Environ. Manage. 129, 398–409 (2013).

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 8 www.nature.com/scientificreports/ www.nature.com/scientificreports

10. Doherty, T. S., Dickman, C. R., Nimmo, D. G. & Ritchie, E. G. Multiple threats, or multiplying the threats? Interactions between invasive predators and other ecological disturbances. Biol. Conserv. 190, 60–68 (2015). 11. Courchamp, F., Langlais, M. & Sugihara, G. Rabbits killing birds: Modelling the hyperpredation process. J. Anim. Ecol. 69, 154–164 (2000). 12. DeCesare, N. J., Hebblewhite, M., Robinson, H. S. & Musiani, M. Endangered, apparently: Te role of apparent in endangered species conservation. Animal Conservation 13, 353–362 (2010). 13. Doherty, T. S., Glen, A. S., Nimmo, D. G., Ritchie, E. G. & Dickman, C. R. Invasive predators and global biodiversity loss. Proc. Natl. Acad. Sci. 113, 11261–11265 (2016). 14. Fitzgerald, B. M., Karl, B. J. & Veitch, C. R. Te diet of feral cats (Felis catus L.) on Raoul Island, Kermadec Group. New Zeal. J. Zool. 15, 123–129 (1991). 15. Courchamp, F., Chapuis, J. L. & Pascal, M. Mammal invaders on islands: Impact, control and control impact. Biol. Rev. Camb. Philos. Soc. 78, 347–383 (2003). 16. Bonnaud, E. et al. Te diet of feral cats on islands: A review and a call for more studies. Biological Invasions 13, 581–603 (2011). 17. Campbell, K. J., Harper, G. A. & Algar, L. Review of feral cat eradications on islands. In Island invasives: eradication and management (eds Veitch, C. R., Clout, M. N. & Towns, D. R.) 37–46 (IUCN, 2011). 18. Medina, F. M. et al. A global review of the impacts of invasive cats on island endangered . Glob. Chang. Biol. 17, 3503–3510 (2011). 19. Palmas, P. et al. Feral cats threaten the outstanding endemic fauna of the New Caledonia biodiversity hotspot. Biol. Conserv. 214, 250–259 (2017). 20. Donlan, C. J. & Wilcox, C. Diversity, invasive species and extinctions in insular ecosystems. J. Appl. Ecol. 45, 1114–1123 (2008). 21. Grant, W. E. & Longnecker, M. T. Efect of subsidized house cats on Calfornia birds and . Trans. West. Sect. Wildl. Soc. 35, 29–33 (1999). 22. Bateman, P. W. & Fleming, P. A. Big life: Carnivores in urban environments. J. Zool. 287, 1–23 (2012). 23. Newsome, T. M. et al. Te ecological efects of providing resource subsidies to predators. Glob. Ecol. Biogeogr. 24, 1–11 (2015). 24. Loss, S. R., Will, T. & Marra, P. P. Te impact of free-ranging domestic cats on wildlife of the United States. Nat. Commun. 4, https:// doi.org/10.1038/ncomms2380 (2013). 25. Okano, T. Biodiversity in the islands of Kagoshima. In Te Islands of Kagoshima (eds Kawai, K., Terada, R. & Kuwahara, S.) 136–145 (Kagoshima University Research Center for the Pacifc Islands, 2013). 26. Jogahara, T., Ogura, G., Sasaki, T., Takehara, K. & Kawashima, Y. Food habits of cats (Felis catus) in forests and villages and their impacts on native animals in the Yambaru area, northern part of Okinawa Island, Japan. Mamm. Sci. 43, 29–37 (2003). 27. Shionosaki, K., Yamada, F., Ishikawa, T. & Shibata, S. Feral cat diet and predation on endangered endemic mammals on a biodiversity hot spot (Amami-Ohshima Island, Japan). Wildl. Res. 42, 343–352 (2015). 28. Yamada, F. et al. Rediscovery afer thirty years since the last capture of the critically endangered Okinawa spiny rat Tokudaia muenninki in the northern part of Okinawa Island. Mammal Study 35, 243–255 (2010). 29. Shionosaki, K., Sakaki, S., Yamada, F. & Shibata, S. Changes in free-roaming cat activity following a regulation prohibiting feeding: A case study at a mountain forest near residential area in Amami City on Amami-Ohshima Island, Japan. Wildl. Hum. Soc. 3, 1–13 (2016). 30. Japan Wildlife Research Center. Annual report 2014: A survey on invasive mammal inhabitation for the ecosystem restoration of Tokunoshima Island. Ministry of the Environment & Japan Wildlife Research Center 1–96 (2015). 31. Sawanobori, R. Feral cat issue in Tokunoshima Island; II. Management by the Ministry of Environment. In Feral cats in Amami Islands (ed. Kagoshima Environmental Study Group) 120–126 (Nanpo-shinsya, 2019). 32. Levy, J. K. & Crawford, P. C. Humane strategies for controlling feral cat populations. J. Am. Vet. Med. Assoc. 225, 1354–1360 (2005). 33. Kagoshima Prefecture. The overview of Amami Islands. Available at, https://www.pref.kagoshima.jp/aq01/20170403.html (Accessed: 11th January 2018) (2017). 34. Jones, E. & Coman, B. J. Ecology of the Feral Cat, Felis catus (L.), in South-Eastern III. Home Ranges and in Semiarid North-West . Austrarian Wildl. Res. 9, 409–420 (1982). 35. Moseby, K. E., Stott, J. & Crisp, H. Movement patterns of feral predators in an arid environment – implications for control through poison baiting. Wildl. Res. 36, 422 (2009). 36. Horn, J. A., Mateus-Pinilla, N., Warner, R. E. & Heske, E. J. Home range, habitat use, and activity patterns of free-roaming domestic cats. J. Wildl. Manage. 75, 1177–1185 (2011). 37. Burrows, N. D. et al. Controlling introduced predators in the Gibson Desert of Western Australia. J. Arid Environ. 55, 691–713 (2003). 38. Barrett, R. T. et al. Diet studies of : A review and recommendations. ICES Journal of Marine Science 64, 1675–1691 (2007). 39. McCutchan, J. H., Lewis, W. M., Kendall, C. & McGrath, C. C. Variation in trophic shif for stable isotope ratios of carbon, nitrogen, and sulfur. Oikos 102, 378–390 (2003). 40. Sugimura, K., Yamada, F. & Miyamoto, A. Population trend, habitat change and conservation of the unique wildlife species on Amami Island, Japan. Glob. Environ. Res. 7, 79–89 (2003). 41. Ohdachi, S. D., Ishibashi, Y., Iwasa, M. A., Fukui, D. & Saitoh, T. Te wild mammals of Japan. (SHOUKADOH Book Sellers,Te Mammalogical Society of Japan, 2015). 42. Kambe, Y. et al. Introgressive hybridization of two major lineages of invasive black rats, Rattus rattus and R. tanezumi on the Japanese Islands inferred from Mc1r sequences. Mamm. Sci. 53, 289–299 (2013). 43. Suárez-Esteban, A., Delibes, M. & Fedriani, J. M. Barriers or corridors? Te overlooked role of unpaved roads in endozoochorous seed dispersal. J. Appl. Ecol. 50, 767–774 (2013). 44. Recio, M. R., Seddon, P. J. & Moore, A. B. Niche and movement models identify corridors of introduced feral cats infringing ecologically sensitive areas in . Biol. Conserv. 192, 48–56 (2015). 45. Nogales, M. et al. A review of feral cat eradication on islands. Conserv. Biol. 18, 310–319 (2004). 46. Yirga, G. et al. Adaptability of large carnivores to changing anthropogenic food sources: Diet change of spotted hyena (Crocuta crocuta) during Christian fasting period in northern Ethiopia. J. Anim. Ecol. 81, 1052–1055 (2012). 47. Meckstroth, A. M. & Miles, A. K. Diets of introduced predators using stable isotopes and stomach contents. J. Wildl. Manage. 71, 2387–2392 (2007). 48. Lanszki, J. et al. Feeding habits of house and feral cats (Felis catus) on small Adriatic islands (Croatia). North. West. J. Zool. 12, 336–348 (2016). 49. Polito, M. J. et al. Integrating stomach content and stable isotope analyses to quantify the diets of pygoscelid penguins. PLoS One 6, https://doi.org/10.1371/journal.pone.0026642 (2011). 50. Kimura, M. Paleogeography of the Ryukyu Islands. Tropics 10, 5–24 (2000). 51. Shibata, H. et al. The taxonomic position and the unexpected divergence of the Habu viper, Protobothrops among Japanese subtropical islands. Mol. Phylogenet. Evol. 101, 91–100 (2016). 52. Plantinga, E. A., Bosch, G. & Hendriks, W. H. Estimation of the dietary nutrient profile of free-roaming feral cats: possible implications for nutrition of domestic cats. Br. J. Nutr. 106, 35–48 (2011). 53. Tieszen, L. L., Boutton, T. W., Tesdahl, K. G. & Slade, N. A. Fractionation and turnover of stable carbon isotopes in animal tissues: Implications for δ13C analysis of diet. Oecologia 57, 32–37 (1983).

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 9 www.nature.com/scientificreports/ www.nature.com/scientificreports

54. Phillips, D. L. et al. Best practices for use of stable isotope mixing models in food-web studies. Can. J. Zool. 92, 823–835 (2014). 55. Seki, S. Population trends and molecular phylogenetic position of the isolated Ryukyu Robin Erithacus komadori on the Danjo Islands. Japanese J. Ornithol. 58, 18–27 (2009). 56. Bonnaud, E. et al. Feeding ecology of a feral cat population on a small mediterranean island. J. Mammal. 88, 1074–1081 (2007). 57. Liberg, O. Correction factors for important prey categories of domestic cats. Acta Teriol. 27, 115–122 (1982). 58. Graham, L. H., Swanson, W. F. & Brown, J. L. Chorionic gonadotropin administration in domestic cats causes an abnormal endocrine environment that disrupts oviductal embryo transport. Teriogenology 54, 1117–1131 (2000). 59. Nagy, K. A. Field metabolic rate and food requirement scaling in mammals and birds. Ecol. Monogr. 57, 111–128 (1987). 60. Keitt, B. S., Wilcox, C., Tershy, B. R., Croll, D. A. & Donlan, C. J. Te efect of feral cats on the population viability of black-vented shearwaters (Pufnus opisthomelas) on Natividad Island, Mexico. Anim. Conserv. 5, 217–223 (2002). 61. Fitzgerald, B. M. & Karl, B. J. Food of the feral house cat, Felis catus in forest of the Orongorongo Valley, Wellington, New Zealand. New Zeal. J. Zool. 6, 107–126 (1979). 62. Oshima, S. Types of pet food (3): Classifcation of pet food by form (or proceeding method): 1. Dry Foods. J. Pet. Anim. Nutr. 18, 117–120 (2016). 63. Mizukami, N. R., Goto, M., Izumiyama, S., Hayashi, H. & Yoh, M. Estimation of feeding history by measuring carbon and nitrogen stable isotope ratios in hair of Asiatic black bears. Ursus 16, 93–101 (2005). 64. Mizukami, R. N. et al. Temporal diet changes recorded by stable isotopes in Asiatic black bear (Ursus thibetanus) hair. Isotopes Environ. Health Stud. 41, 87–94 (2005). 65. DeNiro, M. J. & Epstein, S. Infuence of diet on the distribution of carbon isotopes in animals. Geochim. Cosmochim. Acta 42, 495–506 (1978). 66. Bond, A. L. & Jones, I. L. A practical introduction to stable-isotope analysis for biologists: Approaches, cautions and caveats. Mar. Ornithol. 37, 183–188 (2009). 67. Parnell, A. & Jackson, A. Package ‘siar’: Stable isotope analysis in R. Available at, https://cran.r-project.org/web/packages/siar/index. html (Accessed: 10th July 2018) (2015). 68. Miyake, Y. & Wada, E. Te abundance ratio of 15N/14N in marine environments. Records of Oceanographic Works in Japan 9, 37–53 (1967). 69. Vanderklif, M. A. & Ponsard, S. Sources of variation in -diet δ15N enrichment: A meta-analysis. Oecologia 136, 169–182 (2003). 70. Bond, A. L. & Diamond, A. W. Recent Bayesian stable-isotope mixing models are highly sensitive to variation in discrimination factors. Ecol. Appl. 21, 1017–1023 (2011). 71. Hervías, S. et al. Assessing the impact of introduced cats on island biodiversity by combining dietary and movement analysis. J. Zool. 292, 39–47 (2014). 72. Lilith, M., Claver, M. & Garkaklis, M. Roaming habits of pet cats on the suburban fringe in Perth, Western Australia: what size bufer zone is needed to protect wildlife in reserve? In Too Close to Comfort: Contentious Issues in Human–Wildlife Encounters, Royal Zoological Society of New South Wales, (eds Lunney, D., Munn, A. & Melkle, W.) 65–72 (Royal Zoological Society of New South Wales, 2008). 73. Dray, S. et al. Package ‘adespatial’: Multivariate multiscale spatial analysis. Available at, https://cran.r-project.org/package=adespatial (Accessed: 10th July 2018) (2018). 74. Barton, K. Package ‘MuMIn’. Available at, https://cran.r-project.org/web/packages/MuMIn/index.html (Accessed: 22nd July 2018) (2018). 75. Yamada, F. Lagomorpholody: Biology of evasion and escaping strategy. 182–187 (University of Tokyo Press, 2017). 76. Watari, Y., Takatsuki, S. & Miyashita, T. Efects of exotic (Herpestes javanicus) on the native fauna of Amami-Oshima Island, southern Japan, estimated by distribution patterns along the historical gradient of mongoose invasion. Biol. Invasions 10, 7–17 (2008). 77. Abe, H. Classifcation and biology of Japanese Insectivora(Mammalia): I. Studies on variation and classifcation. J. Fac. Agric. Hokkaido Univ. 55, 191–265 (1967). 78. Fujimoto, T., Miyazaki, Y. & Sawa, N. Egg-laying performance test on diferent breeds of egg-laying hens. Bull. Tokushima. Pref. Lives. Res. Ins. 9, 40–50 (2007). 79. Asai, S. & Asai, S. Te bush warbler and other birds inhabiting the Tokyo University Forest in Chichibu in 2005. J. Yamashina Inst. Ornithol. 40, 144–153 (2009). 80. Parnell, A. C., Inger, R., Bearhop, S. & Jackson, A. L. Source partitioning using stable isotopes: Coping with too much variation. PLoS One 5, 9672 (2010). Acknowledgements Te authors are grateful to T. Miyashita for insightful discussion and review of the manuscript. We thank R. Sawanobori, K. Kazato and members of the NPO Tokunoshima Nijinokai, the Tokunoshima Cat Management Council, and Tokunoshima Animal Hospital for their hospitality during our visit and their kind help collecting the samples. We also thank A. Nakajima for assistance with the analysis and the datasets. Tis study was supported by the Environment Research and Technology Development Fund (4-1804) of the Environmental Restoration and Conservation Agency of Japan, and the Sumitomo Foundation Fiscal 2017 Grant for Environmental Research Projects. Author contributions F. Yamada and Y. Watari designed the concept and managed the project. K. Shionosaki, R. Nakashita and T. Maeda collected samples and carried out the experiments. T. Maeda conducted the analysis and R. Nakashita, T. Maeda and Y. Watari interpreted the results. T. Maeda wrote the manuscript with help from F. Yamada and Y. Watari. All authors have approved the fnal version of the manuscript and agree to be accountable for all aspects of the work related to the accuracy or integrity of any part of the work. Competing interests Te datasets generated and analyzed during the current study are available from the corresponding author upon a reasonable request. Additional information Supplementary information is available for this paper at https://doi.org/10.1038/s41598-019-52472-3. Correspondence and requests for materials should be addressed to Y.W.

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 10 www.nature.com/scientificreports/ www.nature.com/scientificreports

Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits 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 Cre- ative Commons license, and indicate if changes were made. Te images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© Te Author(s) 2019

Scientific Reports | (2019)9:16200 | https://doi.org/10.1038/s41598-019-52472-3 11