Journal of Threatened Taxa | www.threatenedtaxa.org | 17 February 2020 | 12(2): 15229–15237

 

 

ISSN 0974-7907 (Online) | ISSN 0974-7893 (Print) 

doi: https://doi.org/10.11609/jott.5506.12.2.15229-15237

#5506 | Received 28 October 2019 | Final received 25 January 2020 | Finally accepted 10 February 2020

 

 

Iberian Lynx Lynx pardinus Temminck, 1827 (Mammalia: Carnivora: Felidae) in central Spain: trophic niche of an isolated population

 

Pedro Alfaya 1, Ariadna Invernón 2 & Germán Alonso 3

 

1,2,3 Department of Biodiversity, Ecology and Evolution, Complutense University of Madrid. Avda. José Antonio Nováis 12, 28040 Madrid, Spain.

1 palfaya@ucm.es (correspondence author), 2 arinvern@ucm.es, 3 galonso@ucm.es

 

 

 

Editor: Angie Appel, Wild Cat Network, Bad Marienberg, Germany.        Date of publication: 17 January 2020 (online & print)

 

Citation: Alfaya, P., A. Invernon & G. Alonso (2020). Iberian Lynx Lynx pardinus Temminck 1827 (Carnivora: Felidae) in central Spain: trophic niche of an isolated population.  Journal of Threatened Taxa 12(2): 15229–15237. https://doi.org/10.11609/jott.5506.12.2.15229-15237

 

Copyright: © Alfaya et al. 2020. Creative Commons Attribution 4.0 International License.  JoTT allows unrestricted use, reproduction, and distribution of this article in any medium by providing adequate credit to the author(s) and the source of publication.

 

Funding: Spanish Ministry of Agriculture, Food and Environment.

 

Competing interests: The authors declare no competing interests.

 

Author details: Pedro Alfaya is a honorary research fellow in the Ecology Department of the Universidad Complutense de Madrid. He is interested in species conservation, landscape ecology and statistical modelling, but currently he is working on environmental impact consulting. Ariadna Invernón: she is a Msc in ecology and collaborator of the Ecology Department of the Universidad Complutense de Madrid. Germán Alonso is a lecturer of Ecology in the Universidad Complutense Madrid. He has alternated his research activity with public positions, being always linked to environmental conservation, landscape ecology and land management.

 

Author contribution: PA made the bibliographic review, statistical analysis, wrote the manuscript and participated in scat collection.  AI designed and carried out the identification of scat samples content.  GA provides support and experience in field work, identification of scat content and also reviewed the manuscript.

 

Acknowledgements: We thank the Spanish Ministry of Agriculture, Food and Environment (16MNSV002) for funding our study.  We also appreciate the effort of Pilar, Beatriz, and Alejandra for their hard work in the laboratory, and also to Javier and Jorge for their assistance in the fieldwork.

 

 

 

Abstract: Understanding predator-prey relationships is fundamental to develop effective conservation plans.  Between 2015 and 2018, we combed 21 transects, each 7km long, searching for Iberian Lynx Lynx pardinus scat within the province of Madrid in central Spain.  In order to minimise inherent subjectivity of visual identification as much as possible, we performed a double specific nested polymerase chain reaction (PCR) followed by a primer extension assay addressed to two Iberian Lynx diagnostic single nucleotide polymorphisms.  Forty-six scat samples were positively identified as belonging to Iberian Lynx through genetic analysis.  From these, we extracted remains of consumed prey, which we determined to the lowest possible taxonomic level, mainly through hair identification.  Identified prey was divided into four types: lagomorphs, small mammals, birds, and ungulates.  The species’ diet composition was described based on the frequency of occurrence (FO) of each prey and niche breadth, and also compared with prior knowledge of the species using four prior studies as a comparative reference through the calculation of the niche overlap value.  The FO of lagomorphs (39%) was the lowest, while the FO of small mammals (54%) was the highest recorded to date.  The niche breath (0.36) was higher than recorded in prior studies, but still showing the specialist character of the Iberian Lynx.  Niche overlap was low (C = 0.49), showing differences in trophic niche between the population in our study area and the one studied in southern Spain.  This indicates that the Iberian Lynx is adept at switching its main prey, an ability that has previously been firmly rejected.  It is, however, capable of adapting to alternative prey more often than recorded to date, which could be a behavioural response to the patchy distribution of European Rabbit Oryctolagus cuniculus in the study area.

 

Keywords: Diet, lagomorphs, niche breadth, niche overlap, single nucleotide polymorphism.

 

 

 

 

Introduction

 

The Iberian Lynx Lynx pardinus is endemic to the Iberian peninsula (Rodríguez & Delibes 1992), and is regarded as a trophic super specialist (Ferrer & Negro 2004).  Since the 1950s, the Iberian Lynx population has declined continuously (Valverde 1963; Cabezas-Díaz et al. 2009).  Only 93 individuals were recorded in 2002 (Guzmán et al. 2004).  Following conservation measures such as reintroductions of captive-bred Iberian Lynxes in southern Spain, this population experienced a constant growth (Simón et al. 2011; Rodríguez & Calzada 2015), reaching 589 individuals in 2017 (Simón 2018).  Furthermore, Cruz et al. (2019) confirmed the presence of Iberian Lynxes outside the currently known range of the species in the southern Iberian peninsula, suggesting the continued existence of a stable population in central Spain within the province of Madrid.

The diet of a species is a fundamental aspect of its ecology that depends mainly on the abundance and availability of prey types (Terraube & Arroyo 2011), but also on learning and experience of individuals (Shipley et al. 2009).  A widespread phenomenon in many vertebrate and invertebrate taxa (Bolnick et al. 2003) is the so- called ‘niche variation hypothesis’.  This occurs when some co-occurring individuals of a species actively select different prey types in their shared environment (Araujo et al. 2011).  The niche variation could be a response to two main factors: (i) change in environmental conditions that affects prey availability and prompts all individuals of a population to use a larger spectrum of resources, or (ii) each individual continues to use a narrow range of resources that diverges from conspecifics, thus minimizing the interspecific competition (Costa et al. 2008).

Understanding predator-prey relationships is fundamental to identify conservation priorities, prior to the design of conservation programmes for vulnerable or endangered species (Popp et al. 2018).  Lacking information on these relationships could result in ill-informed conservation strategies that lead to a failure of reaching conservation goals and at the same time to a gross waste of resources, as occurred in the Doñana National Park with the restocking of European Rabbits Oryctolagus cuniculus (Carro et al. 2019).

The knowledge about diet plasticity of a species is a keystone to assess the transferability of results obtained in a certain area to another, and for assessing whether an alternative management will provide similar results (Terraube & Arroyo 2011).  A relevant descriptor of niche is breadth, which is a function of the proportion of each resource used with regard to total consumed resources (Smith 1982).  Therefore, a species that uses a wide range of trophic resources in a similar proportion will show a high niche breadth and, consequently, will be regarded as a generalist for studied resources (Symondson et al. 2002).  On the contrary, a species that uses a high proportion of a narrow range of resources will be regarded as a specialist (Shipley et al. 2009).

Rodríguez & Delibes (1992) were the last authors who reported an Iberian Lynx population in the province of Madrid before Cruz et al. (2019).  The territorial and solitary behaviour of the Iberian Lynx (San Miguel 2006; Calzada et al. 2007; Martín et al. 2007) results in a low-density spatial organization that makes it extremely difficult to find and track (Alfaya et al. 2019).  The central Spanish population was already small in the early 1990s (Rodríguez & Delibes 1992), remained elusive and was not considered in conservation programs initiated in 2002 (Rodríguez & Calzada 2015) that lead to the recovery of the population in southern Spain (Simón 2018).

Different studies carried out in Doñana National Park and in Sierra Morena showed that the European Rabbit is its main prey, being present in 70–99% of analysed samples (Delibes 1980; Beltrán & Delibes 1991; Palomares et al. 2001; Gil-Sánchez et al. 2006).  The Iberian Lynx, however, also consumes other prey species in lower proportions, but their relative importance increases when the availability of Rabbits decreases (Beltrán et al. 1985; Beltrán & Delibes 1991).  Therefore, Iberian Lynxes also prey on small mammals, e.g., Wood Mouse Apodemus sylvaticus and Garden Dormouse Eliomys quercinus (Aymerich 1982; Gil-Sánchez et al. 2006), wild ungulates like Red Deer Cervus elaphus, Fallow Deer Dama dama (Delibes 1980; Beltrán et al. 1985; Beltrán & Delibes 1991; Gil-Sánchez et al. 2006), European Mouflon Ovis orientalis musimon (Gil-Sánchez et al. 2006), and birds such as Red-legged Partridge Alectoris rufa, Mallard Anas platyrhynchos and Eurasian Magpie Pica pica (Delibes 1980; Aymerich 1982; Beltrán & Delibes 1991; Gil-Sánchez et al. 2006).

In this article, we report the diet composition of Iberian Lynx, based on analysis of scat collected in a study area in central Spain.  We discuss the trophic niche breadth of this population in the light of research conducted on the species’ diet in southern Spain.

 

Study area

 

The research was performed in the western region of the province of Madrid (Figure 1), which is delimited by boundaries with the community of Castilla–León in the north and northwest, the community of Castilla–La Mancha in the south, and the Manzanares River basin in the east.  The study area ranges in elevation from 440 to 2,320 m.  It represents three main landscape regions: (i) the Guadarrama Mountains, a mountainous granitic zone, (ii) the foothills with a gradient of siliceous sand and soft slopes, and (iii) the depression, a terrain characterised by interfluvial hills (Rivas-Martínez 1982; Zabía & del Olmo 2007).

The meso-Mediterranean zone is the dominant bioclimatic belt within the study area, but the oro- and supra-Mediterranean zones are also present in the Guadarrama Mountains (Rivas-Martínez 1982).  The main climatic features of the study area are the seasonal variation in temperature between -8°C and 44°C, summer drought and irregular precipitation ranging from 400 to 2,000 mm per year (Zabía & del Olmo 2007).

The landscape in the study area is a mosaic of pastures with scrub and Holm Oak Quercus ilex groves interspersed with villages and patches of agricultural land (Schmitz et al. 2007; Image 1).  Local people use pastures traditionally during the summer for grazing transhumant cattle, periodically perform selective logging in the forests and clean the understorey (Arnaiz-Schmitz et al. 2018).

 

Material and Methods

 

Sample collection

Evidence of the presence of Iberian Lynxes within the study area was collected between January 2015 and May 2018.  We designed 21 transects of 7km each that were combed by at least two researchers.  We searched for scat on foot along pathways and firebreaks, since both Iberian Lynx and European Wildcat Felis silvestris usually move along such linear structures (Lozano et al. 2013; Garrote et al. 2014).  Sometimes, we also combed other less regular landscape features such as the bases of large rocks, around Rabbit holes and near rivers, where scat was more likely to be found (Martín et al. 2007).  Along these transects, we searched for scat that is morphologically compatible with scat of the Iberian Lynx.  This ranges in length from 5cm to 9cm and in width from 1.5cm to 2cm, and is divided into several fragments (Rodríguez 1993).  It ranges in colour from ash-grey to dark-brown, and is entirely covered by a mucous patina when fresh (Iglesias & España 2010).

The probability of an erroneous identification of the scat of Iberian Lynx, however, is high (Boshoff & Kerley 2010, Molinari-Jobin et al. 2012; Garrote & de Ayala 2015).  It has been often misidentified due to its similarity with the scat of European Wildcat and Red Fox Vulpes vulpes (Palomares et al. 2002).  To reduce this probability as much as possible, we performed a specific genetic identification analysis designed by Cruz et al. (2019).  This genetic analysis consists of a double specific nested PCR followed by a primer extension assay addressed to two Iberian Lynx diagnostic single nucleotide polymorphisms (SNPs).  The product of the double nested PCR is already specific for the Iberian Lynx since we used the primer DL7F [5’-CTT AAT CGT GCA TTA TAC CTTGT-3’] developed by Palomares et al. (2002), which was aligned to sequences of orthologue carnivores including Eurasian Lynx Lynx lynx, Canada Lynx L. canadensis, European Wildcat, and Domestic Cat Felis catus in order to select diagnostic positions.  Then we identified two SNPs specific of the Iberian Lynx.  These SNPs were marked with fluorescence and detected through a capillary electrophoresis.  This method of analysis provides an increase of sensitivity and straightforward verification of the belonging species through the diagnostic SNPs, being strongly protected against false positive results.  For further details see Cruz et al. (2019).

 

Content analysis

Subsequent to positive genetic identification of scat samples as belonging to the Iberian Lynx, we analysed the contents of these samples.  We used a stereomicroscope to identify and remove remains of consumed prey like broken bones, teeth, feathers, and hair.  Teeth and bone remains were identified with a stereomicroscope, while feathers and hair required the use of a 40x microscope.  We washed hairs, first with distilled water and detergent, and then with 70% alcohol as described in Teerink (1991).  After drying hairs, we poured a thin layer of transparent nail varnish over a slide and let it dry for 30 seconds.  Then we put each hair on the slide for 30 minutes and covered it with a cover glass.  That way, we obtained a hair cuticle mould with a scale pattern showing a certain, although limited, taxonomic value (Short 1978).

Removed remains were identified up to the family level, except those belonging to Wild Boar Sus scrofa, because of their easy identification.  We identified hair using Teerink (1991) and Valente et al. (2015), teeth using Dueñas et al. (1985), and feathers using Dove & Koch (2011).

 

Diet composition

Hutchinson (1957) defined the niche as an n-dimensional hypervolume where distribution of environmental variables and/or factors would allow a certain species to exist indefinitely.  This approach provides a quantitative perspective of the niche concept and, therefore, established the conceptual basis for the performance of studies in many different fields of ecology (Smith 1982).  We defined the trophic niche as the n-dimensional hypervolume, n being the number of prey types consumed by the target species, constrained by used trophic resources that would allow the species to exist indefinitely.

For diet description we grouped consumed prey into four categories: birds, lagomorphs, small mammals and ungulates.  We calculated the frequency of occurrence (FO) for each category regarding total analysed scat samples, and also the niche breadth using Levins Index (Levins 1968):

B = 1/(∑ [pi2])

where pi is the proportion of occurrence of the prey category i, regarding the total consumed prey.  To compare this with other populations, we used the standardisation suggested by Colwell & Futuyma (1971),

Bstand= (B-1)/ (n-1)

where n is the number of prey categories consumed.  This index shows the degree of specialisation of a certain species; a value close to 0 is indicative of a specialist predator, while a generalist predator shows values close to 1 (Colwell & Futuyma 1971).

Both FO and Bstand calculated from analysed scat samples were compared with prior knowledge of the trophic ecology of the Iberian Lynx.  For that, we selected four relevant studies as a comparative reference (Table 1), and regrouped their results to our four prey types.  This was not possible in regard to the study by Fedriani et al. (1999) who used a broad classification of prey items, e.g., other vertebrates, referring to all non-lagomorph vertebrates.  Therefore, we calculated the FO of each prey category and Bstand for all four reference studies, and compared results with those obtained in our study area.

We compared the trophic niche of the Iberian Lynx population in the province of Madrid (M) with that described in prior studies (A).  For the latter, we calculated the average FO of each prey category in the reference studies. Then we used the index formulated by Schoener (1970) for calculating the niche overlap between both populations, M and A:

C = 1- ½ ∑|piMpiA|

where piM is the proportion of occurrence of the category i within population M, and piA is the same but within population A.  C takes a minimum value of 0 when there is no overlap, and a maximum of 1 when the proportions of consumed resources are the same in both populations.

Lastly, we compared FO and Bstand obtained in the province of Madrid between the two periods when samples were collected, i.e., spring–summer and autumn–winter.  We used Fisher’s exact test, which is suitable for small sampling sizes.

 

 

Results

 

Between January 2015 and May 2018, we collected 98 scat samples along 21 transects that were each combed twice, once in spring–summer from May to July and once in autumn–winter from October to February.  Through genetic analysis we identified 46 of these samples positively as belonging to Iberian Lynx, with 31 collected in the spring–summer season and 15 collected in the autumn–winter season.  As our genetic method allows only for identifying scat of Iberian Lynx, we did not attempt to identify the remaining 52 scat samples to other species.

The content analysis of the 46 scat samples revealed an overall niche breadth Bstand of 0.36, with small mammals constituting the majority of prey categories.  Fisher’s exact test shows the existence of marginally non-significant (p = 0.07) differences in diet composition between both seasons of the year considered (Figure 2).  Details are provided in Table 2.

 

 

Discussion

 

The Bstand (0.36) calculated for the population in our study area shows the specialist character of the Iberian Lynx. However, this value is higher than those obtained for comparative reference studies in Table 1.  Furthermore, the obtained C value of 0.49 shows the trophic niche shift of this  population, with regard to that known so far.  Figure 3 shows that the trophic niche of the Iberian Lynxes within the study area is directed towards predation on small mammals.

The Iberian Lynx is regarded as a trophic specialist, strictly dependent of the European Rabbit (Delibes 1980; Aymerich 1982; Beltrán et al. 1985; Beltrán & Delibes 1991; Calzada & Palomares 1996; Palomares et al. 2001; Ferrer & Negro 2004; Gil-Sánchez et al. 2006).  There are such strong links between these two species that the collapse of Rabbit populations can even inhibit the reproductive capability of the Iberian Lynx which has been interpreted as its ‘inability’ to switch its main prey (Ferreras et al. 2011).  In this research, we compared niche breadth and overlap between a potential population in the central Iberian peninsula and prior knowledge obtained from southern populations.  Our results show differences in comparison with those obtained from the four studies used as a comparative reference (Fig. 3).  Therefore, our study is the first record in which lagomorphs are not the main prey, showing a 30% lower FO than in the lowest record so far (70%, Beltrán & Delibes 1991).  On the contrary, the FO of small mammals is clearly over-represented (47.5% higher) in comparison with prior studies.

A similar pattern than the one observed here was already recorded in Delibes et al. (1975).  In this study carried out in the provinces of Cáceres and Salamanca (closer to our study area than to southern populations), the recorded FO for Rabbits was 56.5% while the small mammals and birds occur in the 27% and 12% of samples, respectively.  These results, although still different to ours, show a pattern of Iberian Lynxes farther inland feeding on alternative prey other than Rabbits more frequently.

Rabbit distribution within our study area shows clear differences between main landscape regions.  The population in the north is naturally fragmented, most likely because of the patchy distribution of suitable habitat (Virgós et al. 2003).  In the south, Rabbits are widespread (Blanco & Villafuerte 1993) due to the existence of a high density of boundaries between croplands and scrublands (Calvete et al. 2004) where they find a suitable combination of trophic resources and shelter (Tapia et al. 2014).  As far as we know there is no more actualized information about new population trends, but the described spatial arrangement coincides with our field observations throughout the sampling period.

The observed pattern in our study area could be a response to: (i) Iberian Lynx adaptation that shows a different trophic behaviour in different environments.  Note that 65% of the Iberian Lynx scat samples analysed were collected in the landscape region of the Guadarrama Mountains, where Rabbit distribution is patchy.  This could lead to the exploration of different trophic niches in areas where Rabbit abundance is lower.  A similar pattern was obtained by Sáez-Gómez et al. (2018) and Nájera et al. (2019), who recorded Iberian Lynxes preying on Red-necked Nightjar Caprimulgus ruficollis eggs and Domestic Cats, respectively, as a response to the decline of Rabbit abundance; (ii) an uncertain proportion of our Iberian Lynx scat samples could come from juvenile individuals, whose habitat requirements are less restrictive than those of resident individuals (Gastón et al. 2016).  Therefore, trophic plasticity could be wider too, which would add some noise to our results; and (iii) overestimations of the FO of small species might have been obtained (Torres et al. 2015).  These have more hair and other indigestible matter per unit of body mass, which can cause their occurrence in a higher number of scat samples per unit of consumed mass (Floyd et al. 1978).  Despite this, earlier studies on the trophic ecology of the Iberian Lynx did not suggest evidence of overrepresentation of small prey (Delibes 1980; Aymerich 1982; Beltrán et al. 1985; Beltrán & Delibes 1991; Calzada & Palomares 1996; Palomares et al. 2001; Ferrer & Negro 2004; Gil-Sánchez et al. 2006).  Therefore results are still comparable.

The observed seasonal variation in the diet of the Iberian Lynx in our study area corroborates results of previous studies on the species (Delibes 1977; Beltrán & Delibes 1991; Gil-Sánchez et al. 2006) as well as on the Eurasian Lynx (Krofel et al. 2011).  Lagomorph predation resulted in a 27% lower value during the spring–summer period, while small mammals consumed showed a 10% increase in comparison with the autumn–winter period.  Bird predation was only recorded in spring–summer (FO = 21%).  Bstand also shows differences between both seasons, being higher in spring–summer (0.58) than in autumn–winter (0.41).  Therefore, during the cold season of the year, the Iberian Lynx consumes a lower variety of trophic resources, whilst this pattern changes in the warm season.

This could be motivated by two facts that are likely to produce a synergic effect: (i) during autumn–winter, when high precipitation and low temperature occur, the daily activity of prey is reduced, being less available for Iberian Lynxes (Beltrán & Delibes 1994).  On the other hand, during the spring–summer season, climatic conditions are less adverse, which allows for an increase in daily activity and, therefore, higher availability of different prey species; (ii) the Rabbit reproduction period begins in October–November and can last until June–July, depending on environmental conditions.  This produces a maximum peak of abundance just before summer.  Then Rabbits become the most abundant prey and, as a consequence, predators apply the highest pressure to a single trophic resource.  Moreover, Rabbits do not reproduce during summer (Soriguer & Palacios 1994).  Therefore, a quick and deep decrease of Rabbits occurs, forcing Iberian Lynxes to prey on alternative trophic resources (Delibes 1980) for the rest of the summer.

Our results reinforce the key role that lagomorphs play in the diet of the Iberian Lynx.  This category is the most frequent prey when diversity of available prey is lower.  Here, however, we provide evidence for a lower trophic dependence of the Iberian Lynx on lagomorphs than in the areas of Doñana-Aljarafe and Andújar-Cardeña.  In our study area, the Iberian Lynx shows its adaptive capacity, adopting a relatively generalist strategy when trophic diversity is high, and a more specialist strategy when diversity is low.  Despite this, the low number of samples collected in autumn–winter season (n=15) must be taken into account and, therefore, the pattern showed here may change with a larger dataset.

Knowledge of predator-prey relationships is fundamental for the adequate design and implementation of species conservation plans (Popp et al. 2018).  Therefore, the results of our research provides base line information for designing conservation actions for the Iberian Lynx in central Spain.  We show that the Iberian Lynx is capable to adapt to a wider prey spectrum than previously assumed by Ferrer & Negro (2004) and Ferreras et al. (2011).  Based on the described pattern, we think that Iberian Lynxes can profit from an increase in prey diversity provided in enrichment programmes carried out at captive breeding centres (Rivas et al. 2016).  Familiarising them with a broader prey diversity may enhance the ability of reintroduced individuals to colonise and survive in new territories.

Future research efforts on the trophic ecology of the Iberian Lynx should focus on increasing the number of scat samples for analysis of diet composition, but also on prey availability and the estimation of ‘real’ proportion each prey species contributes to the diet by means of correction factors, as suggested by Wachter et al. (2012) and Klare et al. (2011).  This will provide more reliable information about the trophic needs of the Iberian Lynx.

 

 

Table 1. Previous studies used as comparative references, with authors, sample size (n), sample periods (years), studied population and niche breadth (Bstand) of respective populations.

Authors

n

Years

Population

Bstand

Delibes 1980

1573

19731976

Doñana-Aljarafe

0.060

Beltrán & Delibes 1991

209

19831984

Doñana-Aljarafe

0.279

Palomares 2001

1171

19931996

Doñana-Aljarafe

0.005

Gil-Sánchez et al. 2006

360

20012002

Andújar-Cardeña

0.035

 

Table 2. Frequency of occurrence (FO) of prey categories in scat samples of Iberian Lynx (n=46) collected in the study area between January 2015 and May 2018.

Prey category

FO in spring–summer

FO in autumn–winter

Total FO

Small mammals

50%

40%

54%

Lagomorphs

26.31%

53.33%

39%

Birds

21.05%

0%

17%

Ungulates

2.63%

6.67%

4%

 

 

For figures  - - click here

 

 

References

 

Alfaya, P., J.G. Casanovas, J. Lobón-Rovira, B. Matallanas, A. Cruz, P. Arana & G. Alonso (2019). Using MaxEnt algorithm to assess habitat suitability of a potential Iberian Lynx population in central Iberian Peninsula. Community Ecology 20(3): 266–267. https://doi.org/10.1556/168.2019.20.3.7

Araujo, M.S., D.I. Bolnick & C.A. Layman (2011). The ecological causes of individual specialisation. Ecology Letters 14: 948–958. https://doi.org/10.1111/j.1461-0248.2011.01662.x

Arnaiz-Schmitz, C., P. Díaz,  D. Ruiz-Labourdette, C. Herrero-Jáuregui, M. Molina, C. Montes, F.D. Pineda & M.F. Schmitz (2018). Modelling of socio-ecological connectivity. The rural-urban network in the surroundings of Madrid (Central Spain). Urban Ecosystems 21: 1199–1212. https://doi.org/10.1007/s11252-018-0797-z

Aymerich, M. (1982). Étude comparative des régimes alimentaires du Lynx Pardelle (Lynx pardina) et du Chat Sauvage (Felis sylvestris). Mammalia 46(4): 515–521. https://doi.org/10.1515/mamm.1982.46.4.515

Beltrán, J.F. & M. Delibes (1991). Ecología trófica del Lince Ibérico en Doñana durante un periodo seco. Doñana Acta Vertebrata 18(1): 113–122.

Beltrán, J.F. & M. Delibes (1994). Environmental determinants of circadian activity of free-ranging Iberian Lynxes. Journal of Mammalogy 75(2): 382–393.

Beltrán, J.F., C.S. José, M. Delibes & F. Braza (1985). An analysis of the Iberian Lynx predation upon Fallow Deer in the Coto Donana, SW Spain, pp. 961–967. In: Crombrugghe, S. (ed.). Transactions of XVIIth Congress of the International Union of Game Biologists. Ministry of Agriculture, Agricultural Research Administration Brussels, 1155pp.

Blanco, J.C. & R. Villafuerte (1993). Factores ecológicos que influyen sobre las poblaciones de Conejos. Incidencia de la enfermedad hemorrágica. Internal Technical Report Empresa de Transformación Agraria S.A.: 1–122.

Bolnick, D.I., R. Svanbäck, J.A. Fordyce, L.H. Yang, J.M. Davis, C.D. Hulsey & M.L. Forister (2003). The ecology of individuals: incidence and implications of individual specialisation. The American Naturalist 161(1): 1–28.

Boshoff, A.F. & G.I. Kerley (2010). Historical mammal distribution data: How reliable are written records? South African Journal of Science 106(1–2): 26–33.

Cabezas-Díaz, S., J. Lozano & E. Virgós (2009). The decline of the wild Rabbit (Oryctolagus cuniculus) and the Iberian Lynx (Lynx pardinus) in Spain: redirecting conservation efforts, pp. 283–310. In: Aranoff, J.B. (ed.). Handbook of Nature Conservation: Global, Environmental and Economic Issues. Nova Science Publishers, New York, 479pp.

Calvete, C., R. Estrada, E. Angulo & S. Cabezas-Ruiz (2004). Habitat factors related to wild Rabbit conservation in an agricultural landscape. Landscape Ecology 19: 531–542. https://doi.org/10.1023/B:LAND.0000036139.04466.06

Calzada, J. & F. Palomares (1996). Frecuencia de aparición de diferentes restos de conejo en excrementos de Lince y Zorro. Doñana Acta Vertebrata 23(2): 243–252.

Calzada, J., P. Muñoz, A. Sánchez, M.J. Palacios, C. Dávila, M.A. Simón, R. Cadenas, A. Leiva, S. Saldaña, F. Jiménez, R. Carbonell, J.I. Molina, A. Aranda, I. Mosqueda, L. Prada, L.M. González, B. Heredia, M. Aymerich, F.G. Domínguez, M. Delibes, F. Palomares, J.A. Godoy, A. Píriz, J. Calzada, A. Vargas, F. Martínez, J.N. Guzmán, S. Centenera, P. Garzón & J.L. González (2007). Estrategia para la Conservación del Lince Ibérico (Lynx pardinus). Comisión Nacional de Protección de la Naturaleza, Madrid, 44pp.

Carro, F., M. Ortega & C. Soriguer (2019). Is restocking a useful tool for increasing Rabbit densities? Global Ecology and Conservation 17: e00560. https://doi.org/10.1016/j.gecco.2019.e00560

Colwell, R.K. & D.J. Futuyma (1971). On the measurement of niche breadth and overlap. Ecology 52(4): 567–576.

Costa, G.C., D.O. Mesquita, G.R. Colli & L.J. Vitt (2008). Niche expansion and niche variation hypothesis: does the degree of individual variation increase in depauperate assemblages? The American Naturalist 172(6): 868–877.

Cruz, A., B. Matallanas, J. Lobón-Rovira,  J.G. Casanovas, G. Alonso & P. Arana (2019). Double specific nested PCR and diagnostic SNP assay for species identification in lynx fecal critical samples. Conservation Genetics Resources 11(2): 173–175. https://doi.org/10.1007/s12686-018-0993-4

Delibes, M. (1977). Ecología y comportamiento alimenticio de Lynx pardina (Temminck 1824) en el coto de Doñana. PhD Thesis. Universidad Autónoma de Madrid, Madrid.

Delibes, M. (1980). Feeding Ecology of the Spanish Lynx in the Coto Doñana. Acta Theriologica 25(24): 309–324.

Delibes, M., F. Palacios, J. Garzónz & J. Castroviejo (1975). Notes sur l’alimentation et la biologie du Lynx Pardelle (Lynx pardina Temminck, 1824) en Espagne. Mammalia 39: 387–393.

Dove, C.J. & S.L. Koch (2011). Microscopy of feathers: a practical guide for forensic feather identification. The Microscope 59(2): 51–71.

Fedriani, J.M, F. Palomares & M. Delibes (1999). Niche relations among three sympatric Mediterranean carnivores. Oecologia 121(1): 138–148. https://doi.org/10.1007/s004420050915

Ferrer, M. & J.J. Negro (2004). The enar extinction of two large European predators: super specialists pay a price. Conservation Biology 18(2): 344–349. https://doi.org/10.1111/j.1523-1739.2004.00096.x

Ferreras, P., A. Travaini, S.C. Zapata & M. Delibes (2011). Short-term responses of mammalian carnivores to a sudden collapse of Rabbits in Mediterranean Spain. Basic and Applied Ecology 12(2): 116–124. https://doi.org/10.1016/j.baae.2011.01.005

Garrote, G. & R.P. de Ayala (2015). Assessing unverified observation data used for estimating Iberian Lynx distribution. European Journal of Wildlife Research 61(5): 801–806. https://doi.org/10.1007/s10344-015-0941-5

Garrote, G., R.P. de Ayala & J.L. Tellería (2014). A comparison of scat counts and camera-trapping as means of assessing Iberian Lynx abundance. European Journal of Wildlife Research 60(6): 885–889. https://doi.org/10.1007/s10344-014-0855-7

Gastón, A., S. Blázque-Cabrera, G. Garrote, M.C. Mateo-Sánchez, P. Beier, M.A. Simón & S. Saura (2016). Response to agriculture by a woodland species depends on cover type and behavioural state: insights from resident and dispersal Iberian Lynx. Journal of Applied Ecology 53 (3): 814–824.  https://doi.org/10.1111/1365-2664.12629

Gil-Sánchez, J.M., E. Ballesteros-Duperón & J.F. Bueno-Segura (2006). Feeding ecology of the Iberian Lynx Lynx pardinus in eastern Sierra Morena (southern Spain). Acta Theriologica 51(1): 85–90. https://doi.org/10.1007/BF03192659

Guzmán, J.N., F.J. García, G. Garrote, R. Pérez & C. Iglesias (2004). Capítulo VI. Sínteis y Discusión, pp. 133–138. In: El Lince Ibérico (Lynx pardinus) en España y Portugal. Censo-diagnóstico de sus poblaciones. Dirección General para la Biodiversidad, Madrid, 174pp.

Hutchinson, G.E. (1957). Concluding remarks. Cold Spring Quantitative Biology 22: 415–427.

Iglesias, A. & A.J. España (eds.) (2010). Rastros y huellas de carnívoros ibéricos. Ediciones Jaguar, Spain, 350pp.

Klare, U., J.F. Kamler & D.W. Macdonald (2011). A comparison and critique of different scat-analysis methods for determining carnivore diet. Mammal Review 41(4): 294–312. https://doi.org/10.1111/j.1365-2907.2011.00183.x

Krofel, M., D. Huber & M. Kos (2011). Diet of the Eurasian Lynx Lynx lynx in the northern Dinaric Mountains (Slovenia and Croatia). Acta Theriologica 56: 315–322. https://doi.org/10.1007/s13364-011-0032-2

Levins, R. (1968). Evolution in Changing Environments. Princeton University Press, Princeton, 120pp.

Lozano, J., E. Virgós & S. Cabezas-Díaz (2013). Monitoring European Wildcat Felis silvestris populations using scat surveys in central Spain: are population trends related to wild Rabbit dynamics or to landscape features? Zoological Studies 52: 16. https://doi.org/10.1186/1810-522X-52-16

Martín, J., R. Grande, A. Hernando, J. Garzón, J. Eliseo, J. Echegaray, P. García-Buendía, J.G. González-Vázquez, A. Gutierrez, J.A. Manzanedo, F.M. de Lecea, J. Moro, J. Porras, J. Nicolás & J.C. Zamarreño (2007). Prospección del Lince Ibérico en los montes del centro de España. Fundación Abertis, Cantabria, 26pp.

Molinari-Jobin, A., M. Kéry, E. Marboutin, P. Molinari, I. Koren, C. Fuxjäger, C. Breitenmoser-Würsten, S. Wölfl, M. Fasel, I. Kos, M. Wölfl & U. Breitenmoser  (2012). Monitoring in the presence of species misidentification: the case of the Eurasian Lynx in the Alps. Animal Conservation 15(3): 266–273. https://doi.org/10.1111/j.1469-1795.2011.00511.x

Nájera, F., S. Sánchez-Cuerda, G. López, T. del Rey-Wamba, C. Rueda, N. Vallverdú-Coll, J. Panadero, M.J. Palacios, J.V. López-Vao & J. Jiménez (2019). Lynx eats cat: disease risk assessment during an Iberian Lynx intraguild predation. European Journal of Wildlife Research 64: 39. https://doi.org/10.1007/s10344-019-1275-5

Palomares, F., M. Delibes, E. Revilla, J. Calzada & J.M. Fedriani (2001). Spatial ecology of Iberian Lynx and abundance of European Rabbits in southwestern Spain. Wildlife Monographs 148: 1–36.

Palomares, F., J.A. Godoy, A. Píriz & S.J. O’Brien (2002). Faecal genetic analysis to determine the presence and distribution of elusive carnivores: design and feasibility for the Iberian Lynx. Molecular Ecology 11(10): 2171–2182. https://doi.org/10.1046/j.1365-294X.2002.01608.x

Popp, J.N., J. Hamr, J.L. Larkin & F.F. Mallory (2018). Black Bear (Ursus americanus) and Wolf (Canis spp.) summer diet composition and ungulate prey selectivity in Ontario, Canada. Mammal Research 63: 433–441. https://doi.org/10.1007/s13364-018-0368-y

Rivas, A., J. Boixader, A. Vargas, M.J. Pérez, R. Serra, V. Asensio, I., Sánchez, V. Lobato, A. Vázquez, A.J. Pardo, J.M. Chaparro, L. Penabad, J. Reeves, O.M. Jiménez, C. Ferreira, J. Yerga, N. Gonçalves, A. Azevedo, L. Ferreira, J.L. Mendoza, R. González, J. Martín, A. Mamán, E. Díaz. J. Gañán, R. Canales, C. Rojo, J.M. Aguilar, M.J. Coca, M. Cuadrado, L. Flores, M.A. Quevedo, J. Valkenburg, J. Pechém, M. Moreira, M. Fjaere, T. Costa, T. Lopes, N. Sánchez, V. Requeijão, V. Madeira, A. Grancho, S. Agudo, S. Casado, M.J. Fernández, F.J. Gil, L. González, M. González, J.M. Hebrero, A. López, J. Ramírez, L. Coloma, D. Martín, A.A. Esteban, A. Martínez-Cifuentes, A. Martínez-Martín, A. Flores, S. Ruiz, E. Parra, J. Sánchez, O. Muñoz, J. Bergara, F.  Martínez & L. Klink. (2016). Manual de manejo del Lince Ibérico en cautividad. Grupo Asesor de Manejo en Cautividad del Comité de Cría del Lince Ibérico. Programa de Conservación Ex-situ del Lince Ibérico, Jeréz de la Frontera, 248pp.

Rivas-Martínez, S. (1982). Memoria del mapa de las series de vegetación de la provincia de Madrid. Escala 1:200,000. Servicio Forestal, del Medio Ambiente y Contra Incendios, Diputación de Madrid, 47pp.

Rodríguez, A. & J. Calzada (2015). Lynx pardinus. The IUCN Red List of Threatened Species 2015: e.T12520A50655794. Downloaded on 08 February 2020. https://doi.org/10.2305/IUCN.UK.2015-2.RLTS.T12520A50655794.en   

Rodríguez, A. & M. Delibes (1992). Current range and status of the Iberian Lynx Felis pardina Temminck, 1824 in Spain. Biological Conservation 61(3): 189–196.

Rodríguez, J.L. (1993). Guía de campo de los mamíferos terrestres de España. Ediciones Omega, Barcelona, 314pp.

Sáez-Gómez, P., S. Palacios & C. Camacho (2018). Landscape change promotes the emergence of a rare predator-prey interaction. Food Webs 15: e00078. https://doi.org/10.1016/j.fooweb.2018.e00078

San Miguel, A. (Coord.) (2006).  Manual  para  la  gestión del hábitat del Lince Ibérico (Lynx pardinus Temminck) y de su presa principal, el Conejo de Monte (Oryctolagus cuniculus L.). 1ª Reimpresión, Fundación CBD-Habitat, Madrid, 263pp.

Schmitz, M.F., I. de Aranzabal & F.D. Pineda (2007). Spatial analysis of visitor preferences in the outdoor recreational niche of Mediterranean cultural landscapes. Environmental Conservation 34(4): 300–312. https://doi.org/10.1017/S0376892907004249

Schoener, T.W. (1970). Nonsynchronous spatial overlap of lizards in patchy habitats. Ecology 51(3): 408–418.

Shipley, L.A., J.S. Forbey & B.D. Moore (2009). Revisiting the dietary niche: when is a mammal herbivore a specialist? Integrative and Comparative Biology 49(3): 274–290. https://doi.org/10.1093/icb/icp051

Short, L.H. (1978). Analysis of cuticular scales on hairs using the scanning electron microscope. Journal of Mammalogy 59(2): 261–268.

Simón, M.A. (2018). Censo de las poblaciones Andaluzas de Lince Ibérico año 2017. Life+ Iberlince, Junta de Andalucía, Sevilla, 11pp.

Simón, M.A., J.M. Gil-Sánchez, G. Ruiz, G. Garrote, E.B. McCain, L. Fernández, M. López-Parra, E. Rojas, R. Arenas-Rojas, T. del Rey, M. García-Tardío & G. López. (2011). Reverse of the decline of the endangered Iberian Lynx. Conservation Biology 26(4): 731–736. https://doi.org/10.1111/j.1523-1739.2012.01871.x

Smith, E.P. (1982). Niche breadth, resource availability and inference. Ecology 63(6): 1675–1681.

Soriguer, R.C. & F. Palacios (1994). Los lagomorfos ibéricos: liebres y conejos, pp. 63–82. In: García-Granados, A. & L. de Hierro (eds.). Curso de Gestión y Ordenación Cinegética. Colegio Oficial de Biólogos, Granada, 242pp.

Symondson, W.O.C., K.D. Sunderland & M.H. Greenstone (2002). Can generalist predators be effective biocontrol agents? Annual Review of Entomology 47(1): 561–594. https://doi.org/10.1146/annurev.ento.47.091201.145240

Tapia, L., J. Domínguez, A. Regos & M. Vidal (2014). Using remote sensing data to model European wild Rabbit (Oryctolagus cuniculus) occurrence in a highly fragmented landscape in northwestern Spain. Acta Theriologica 59: 289–298. https://doi.org/10.1007/s13364-013-0169-2

Teerink, B.J. (1991). Hair of west European mammals: Atlas and Identification Key. Cambridge University Press, Cambridge, 224pp.

Terraube, J. & B. Arroyo (2011). Factors influencing diet variation in generalist predator across its range distribution. Biodiversity and Conservation 20: 2111–2131. https://doi.org/10.1007/s10531-011-0077-1

Torres, R.T., N. Silva, G. Brotas & C. Fonseca (2015). To eat or not to eat? The diet of endangered Iberian Wolf (Canis lupus signatus) in a human-dominated landscape in central Portugal. PLoS ONE 10 (6): e0129379. https://doi.org/10.1371/journal.pone.0129379

Valente, A.M., R.G. Rocha, E. Lozano, J.P. Ferreira & C. Fonseca (2015). Atlas dos Pelos dos Mamíferos Terrestres Ibéricos. Ediçoes Afrontamento, Porto, 188pp.

Valverde, J.A. (1963). Información sobre el Lince Ibérico en España. Boletín técnico, serie cinegética. Spanish Agricultural Ministry, Madrid, 37pp.

Villafuerte, R. (2007). Oryctolagus cuniculus (Linnaeus, 1758). Ficha Libro Rojo, pp. 487–491. In: Palomo, L.J., J. Gisbert & J.C. Blanco (eds). Atlas y Libro Rojo de los Mamíferos terrestres de España. Dirección General para la Biodiversidad SECEM-SECEMU, Madrid, 588pp.

Virgós, E., S. Cabezas-Díaz, A. Malo, J. Lozano & D. López-Huertas (2003). Factors shaping European Rabbit abundance in continuous and fragmented populations of central Spain. Acta Theriologica 48(1): 113–122. https://doi.org/10.1007/BF03194271

Wachter, B., A.S. Blanc, J. Melzheimer, O.P. Höner, M. Jago & H. Hofer (2012). An advanced method to assess the diet of free-ranging large carnivores based on scats. PLoS One 7: e38066. https://doi.org/10.1371%2Fjournal.pone.0038066

Zabía, M. & L. del Olmo (2007). Atlas Medio Ambiente de la Comunidad de Madrid. Consejería de Medio Ambiente y Ordenación del Territorio, Madrid, 83pp.