Comparison of whole mitochondrial genome sequences of northern and southern white rhinoceroses (<Emphasis Type="Italic">Ceratotherium simum</Emphasis>): the conservation consequences of species definitions

The probable extinction of the last confirmed population of northern white rhinoceros (Ceratotherium simum cottoni) in the world has ignited debate regarding its species status compared to the southern white rhinoceros (Ceratotherium simum simum). Previous studies, based on partial mitochondrial sequences, have reported conflicting results regarding the species status of the northern white rhinoceros. We use whole mitochondrial genomes obtained using Next Generation Sequencing of four northern white rhinoceros and three southern white rhinoceros using novel primers in three overlapping fragments. Phylogenetic relationships were constructed, using Maximum Likelihood, and recovered monophyletic clades for northern white rhinoceros and southern white rhinoceros. The divergence time between the two mitochondrial DNA lineages was estimated to be between 0.46 and 0.97 million years ago using Bayesian inferences. Since there are currently only three surviving northern white rhinoceros individuals these results put into sharp focus the ongoing debate regarding the methods of species definition, and in particular the consequences of such definitions in conservation management of endangered species and subspecies. We conclude that the designation of sub-species status is more applicable to northern and southern white rhinoceros.     

Polymorphic microsatellites in white rhinoceros

The southern white rhinoceros (Ceratotherium simum simum) has suffered severe reductions in population size over the last 150 years as a result of overhunting. We optimized 10 southern white rhinoceros microsatellite loci and tested them on 30 individuals from the largest remaining population of this species. Five of the 10 loci were polymorphic with mean expected heterozygosity of 0.578, mean polymorphic information content of 0.481 and mean number of alleles 2.8. Although these data suggest low genetic variability in C. s. simum, an accurate comparison of variability awaits results of ongoing broad‐scale microsatellite surveys in this family.     

The response of forest plant regeneration to temperature variation along a latitudinal gradient

Background and Aims

The response of forest herb regeneration from seed to temperature variations across latitudes was experimentally assessed in order to forecast the likely response of understorey community dynamics to climate warming.

Methods

Seeds of two characteristic forest plants (Anemone nemorosa and Milium effusum) were collected in natural populations along a latitudinal gradient from northern France to northern Sweden and exposed to three temperature regimes in growth chambers (first experiment). To test the importance of local adaptation, reciprocal transplants were also made of adult individuals that originated from the same populations in three common gardens located in southern, central and northern sites along the same gradient, and the resulting seeds were germinated (second experiment). Seedling establishment was quantified by measuring the timing and percentage of seedling emergence, and seedling biomass in both experiments.

Key Results

Spring warming increased emergence rates and seedling growth in the early-flowering forb A. nemorosa. Seedlings of the summer-flowering grass M. effusum originating from northern populations responded more strongly in terms of biomass growth to temperature than southern populations. The above-ground biomass of the seedlings of both species decreased with increasing latitude of origin, irrespective of whether seeds were collected from natural populations or from the common gardens. The emergence percentage decreased with increasing home-away distance in seeds from the transplant experiment, suggesting that the maternal plants were locally adapted.

Conclusions

Decreasing seedling emergence and growth were found from the centre to the northern edge of the distribution range for both species. Stronger responses to temperature variation in seedling growth of the grass M. effusum in the north may offer a way to cope with environmental change. The results further suggest that climate warming might differentially affect seedling establishment of understorey plants across their distribution range and thus alter future understorey plant dynamics.     

Identifying shared genetic structure patterns among Pacific Northwest forest taxa: insights from use of visualization tools and computer simulations

Background

Identifying causal relationships in phylogeographic and landscape genetic investigations is notoriously difficult, but can be facilitated by use of multispecies comparisons.

Methodology/Principal Findings

We used data visualizations to identify common spatial patterns within single lineages of four taxa inhabiting Pacific Northwest forests (northern spotted owl: Strix occidentalis caurina; red tree vole: Arborimus longicaudus; southern torrent salamander: Rhyacotriton variegatus; and western white pine: Pinus monticola). Visualizations suggested that, despite occupying the same geographical region and habitats, species responded differently to prevailing historical processes. S. o. caurina and P. monticola demonstrated directional patterns of spatial genetic structure where genetic distances and diversity were greater in southern versus northern locales. A. longicaudus and R. variegatus displayed opposite patterns where genetic distances were greater in northern versus southern regions. Statistical analyses of directional patterns subsequently confirmed observations from visualizations. Based upon regional climatological history, we hypothesized that observed latitudinal patterns may have been produced by range expansions. Subsequent computer simulations confirmed that directional patterns can be produced by expansion events.

Conclusions/Significance

We discuss phylogeographic hypotheses regarding historical processes that may have produced observed patterns. Inferential methods used here may become increasingly powerful as detailed simulations of organisms and historical scenarios become plausible. We further suggest that inter-specific comparisons of historical patterns take place prior to drawing conclusions regarding effects of current anthropogenic change within landscapes.     

Rapid growth reduces cold resistance: evidence from latitudinal variation in growth rate, cold resistance and stress proteins

Background

Physiological costs of rapid growth may contribute to the observation that organisms typically grow at submaximal rates. Although, it has been hypothesized that faster growing individuals would do worse in dealing with suboptimal temperatures, this type of cost has never been explored empirically. Furthermore, the mechanistic basis of the physiological costs of rapid growth is largely unexplored.

Methodology/Principal Finding

Larvae of the damselfly Ischnura elegans from two univoltine northern and two multivoltine southern populations were reared at three temperatures and after emergence given a cold shock. Cold resistance, measured by chill coma recovery times in the adult stage, was lower in the southern populations. The faster larval growth rates in the southern populations contributed to this latitudinal pattern in cold resistance. In accordance with their assumed role in cold resistance, Hsp70 levels were lower in the southern populations, and faster growing larvae had lower Hsp70 levels. Yet, individual variation in Hsp70 levels did not explain variation in cold resistance.

Conclusions/Significance

We provide evidence for a novel cost of rapid growth: reduced cold resistance. Our results indicate that the reduced cold resistance in southern populations of animals that change voltinism along the latitudinal gradient may not entirely be explained by thermal selection per se but also by the costs of time constraint-induced higher growth rates. This also illustrates that stressors imposed in the larval stage may carry over and shape fitness in the adult stage and highlights the importance of physiological costs in the evolution of life-histories at macro-scales.     

Population structure of Cynara cardunculus complex and the origin of the conspecific crops artichoke and cardoon

Background and Aims

Globe artichoke and leafy cardoon, two crops within the same species Cynara cardunculus, are traditionally cultivated in the Mediterranean region and play a significant role in the agricultural economy of this area. The two cultigens have different reproductive systems: artichoke is generally vegetatively propagated, while leafy cardoon is seed propagated. The domestication events underlying the origin of both artichoke and cultivated cardoon from their wild relative and the area of occurrence are not yet fully understood. The aim of this study was to investigate population structure in wild cardoon, globe artichoke and leafy cardoon material and infer domestication events.

Methods

Thirty-five microsatellite (simple sequence repeat) markers, distributed in the C. cardunculus genome, and a large geographical and numerical sampling in southern Europe and North Africa were used to assess population structure and diversity.

Key Results

The results suggest the presence of two distinct domestication events for artichoke and leafy cardoon, and also suggest a new possible scenario, with western wild cardoon having originated from cultivated cardoon escaped from cultivation. Evidence was found for a demographic bottleneck in the past history of globe artichoke.

Conclusions

The results shed new light on the relationships between the three taxa of C. cardunculus and highlight relevant aspects on the evolution of domestication of two crops with a different reproductive system within the same species. It is proposed that the probable centre of origin of artichoke is located in southern Italy, probably Sicily.     

Projected Future Distributions of Vectors of Trypanosoma cruzi in North America under Climate Change Scenarios

Background

Chagas disease kills approximately 45 thousand people annually and affects 10 million people in Latin America and the southern United States. The parasite that causes the disease, Trypanosoma cruzi, can be transmitted by insects of the family Reduviidae, subfamily Triatominae. Any study that attempts to evaluate risk for Chagas disease must focus on the ecology and biogeography of these vectors. Expected distributional shifts of vector species due to climate change are likely to alter spatial patterns of risk of Chagas disease, presumably through northward expansion of high risk areas in North America.

Methodology/Principal Findings

We forecast the future (2050) distributions in North America of Triatoma gerstaeckeri and T. sanguisuga, two of the most common triatomine species and important vectors of Trypanosoma cruzi in the southern United States. Our aim was to analyze how climate change might affect the future shift of Chagas disease in North America using a maximum entropy algorithm to predict changes in suitable habitat based on vector occurrence points and predictive environmental variables. Projections based on three different general circulation models (CCCMA, CSIRO, and HADCM3) and two IPCC scenarios (A2 and B2) were analyzed. Twenty models were developed for each case and evaluated via cross-validation. The final model averages result from all twenty of these models. All models had AUC >0.90, which indicates that the models are robust. Our results predict a potential northern shift in the distribution of T. gerstaeckeri and a northern and southern distributional shift of T. sanguisuga from its current range due to climate change.

Conclusions/Significance

The results of this study provide baseline information for monitoring the northward shift of potential risk from Chagas disease in the face of climate change.     

Comparative genomics of methicillin-resistant Staphylococcus aureus ST239: distinct geographical variants in Beijing and Hong Kong

Background

The ST239 lineage is a globally disseminated, multiply drug-resistant hospital-associated methicillin-resistant Staphylococcus aureus (HA-MRSA). We performed whole-genome sequencing of representative HA-MRSA isolates of the ST239 lineage from bacteremic patients in hospitals in Hong Kong (HK) and Beijing (BJ) and compared them with three published complete genomes of ST239, namely T0131, TW20 and JKD6008. Orthologous gene group (OGG) analyses of the Hong Kong and Beijing cluster strains were also undertaken.

Results

Homology analysis, based on highest-percentage nucleotide identity, indicated that HK isolates were closely related to TW20, whereas BJ isolates were more closely related to T0131 from Tianjin. Phylogenetic analysis, incorporating a total of 30 isolates from different continents, revealed that strains from HK clustered with TW20 into the ‘Asian clade’, whereas BJ isolates and T0131 clustered closely with strains of the ‘Turkish clade’ from Eastern Europe. HK isolates contained the typical φSPβ-like prophage with the SasX gene similar to TW20. In contrast, BJ isolates contained a unique 15 kb PT1028-like prophage but lacked φSPβ-like and φSA1 prophages. Besides distinct mobile genetic elements (MGE) in the two clusters, OGG analyses and whole-genome alignment of these clusters highlighted differences in genes located in the core genome, including the identification of single nucleotide deletions in several genes, resulting in frameshift mutations and the subsequent predicted truncation of encoded proteins involved in metabolism and antimicrobial resistance.

Conclusions

Comparative genomics, based on de novo assembly and deep sequencing of HK and BJ strains, revealed different origins of the ST239 lineage in northern and southern China and identified differences between the two clades at single nucleotide polymorphism (SNP), core gene and MGE levels. The results suggest that ST239 strains isolated in Hong Kong since the 1990s belong to the Asian clade, present mainly in southern Asia, whereas those that emerged in northern China were of a distinct origin, reflecting the complexity of dissemination and the dynamic evolution of this ST239 lineage.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-529) contains supplementary material, which is available to authorized users.     

Limited genetic divergence among Australian alpine Poa tussock grasses coupled with regional structuring points to ongoing gene flow and taxonomic challenges

Background and Aims

While molecular approaches can often accurately reconstruct species relationships, taxa that are incompletely differentiated pose a challenge even with extensive data. Such taxa are functionally differentiated, but may be genetically differentiated only at small and/or patchy regions of the genome. This issue is considered here in Poa tussock grass species that dominate grassland and herbfields in the Australian alpine zone.

Methods

Previously reported tetraploidy was confirmed in all species by sequencing seven nuclear regions and five microsatellite markers. A Bayesian approach was used to co-estimate nuclear and chloroplast gene trees with an overall dated species tree. The resulting species tree was used to examine species structure and recent hybridization, and intertaxon fertility was tested by experimental crosses.

Key Results

Species tree estimation revealed Poa gunnii, a Tasmanian endemic species, as sister to the rest of the Australian alpine Poa. The taxa have radiated in the last 0·5–1·2 million years and the non-gunnii taxa are not supported as genetically distinct. Recent hybridization following past species divergence was also not supported. Ongoing gene flow is suggested, with some broad-scale geographic structure within the group.

Conclusions

The Australian alpine Poa species are not genetically distinct despite being distinguishable phenotypically, suggesting recent adaptive divergence with ongoing intertaxon gene flow. This highlights challenges in using conventional molecular taxonomy to infer species relationships in recent, rapid radiations.     

Evolutionary diversification of the lizard genus Bassiana (Scincidae) across Southern Australia

Background

Relatively recent (Plio-Pleistocene) climatic variations had strong impacts on the fauna and flora of temperate-zone North America and Europe; genetic analyses suggest that many lineages were restricted to unglaciated refuges during this time, and have expanded their ranges since then. Temperate-zone Australia experienced less severe glaciation, suggesting that patterns of genetic structure among species may reflect older (aridity-driven) divergence events rather than Plio-Pleistocene (thermally-mediated) divergences. The lizard genus Bassiana (Squamata, Scincidae) contains three species that occur across a wide area of southern Australia (including Tasmania), rendering them ideally-suited to studies on the impact of past climatic fluctuations.

Methodology/Principal Findings

We performed molecular phylogenetic and dating analyses using two partial mitochondrial genes (ND2 and ND4) of 97 samples of Bassiana spp. Our results reveal a pattern of diversification beginning in the Middle Miocene, with intraspecific diversification arising from 5.7 to 1.7 million years ago in the Upper Miocene-Lower Pleistocene.

Conclusions/Significance

In contrast to the temperate-zone Northern Hemisphere biota, patterns of evolutionary diversification within southern Australian taxa appear to reflect geologically ancient events, mostly relating to east-west discontinuities imposed by aridity rather than (as is the case in Europe and North America) relatively recent recolonisation of northern regions from unglaciated refugia to the south.     

Associations of HLA-DP variants with hepatitis B virus infection in southern and northern Han Chinese populations: a multicenter case-control study

Background

Human leukocyte antigen DP (HLA-DP) locus has been reported to be associated with hepatitis B virus (HBV) infection in populations of Japan and Thailand. We aimed to examine whether the association can be replicated in Han Chinese populations.

Methodology/Principal Findings

Two HLA-DP variants rs2395309 and rs9277535 (the most strongly associated SNPs from each HLA-DP locus) were genotyped in three independent Han cohorts consisting of 2 805 cases and 1 796 controls. By using logistic regression analysis, these two SNPs in the HLA-DPA1 and HLA-DPB1 genes were significantly associated with HBV infection in Han Chinese populations (P = 0.021∼3.36×10⁻⁸ at rs2395309; P = 8.37×10⁻³∼2.68×10⁻¹⁰ at rs9277535). In addition, the genotype distributions of both sites (rs2395309 and rs9277535) were clearly different between southern and northern Chinese population (P = 8.95×10⁻⁵ at rs2395309; P = 1.64×10⁻⁹ at rs9277535). By using asymptomatic HBV carrier as control group, our study showed that there were no associations of two HLA-DP variants with HBV progression (P = 0.305∼0.822 and 0.163∼0.881 in southern Chinese population, respectively; P = 0.097∼0.697 and 0.198∼0.615 in northern Chinese population, respectively).

Conclusions

Our results confirmed that two SNPs (rs2395309 and rs9277535) in the HLA-DP loci were strongly associated with HBV infection in southern and northern Han Chinese populations, but not with HBV progression.     

Genetic Diversity of Brazilian Aedes aegypti: Patterns following an Eradication Program

Background

Aedes aegypti is the most important vector of dengue fever in Brazil, where severe epidemics have recently taken place. Ae. aegypti in Brazil was the subject of an intense eradication program in the 1940s and 50s to control yellow fever. Brazil was the largest country declared free of this mosquito by the Pan-American Health Organization in 1958. Soon after relaxation of this program, Ae. aegypti reappeared in this country, and by the early 1980s dengue fever had been reported. The aim of this study is to analyze the present-day genetic patterns of Ae. aegypti populations in Brazil.

Methodology/Principal Findings

We studied the genetic variation in samples of 11 widely spread populations of Ae. aegypti in Brazil based on 12 well-established microsatellite loci. Our principal finding is that present-day Brazilian Ae. aegypti populations form two distinct groups, one in the northwest and one in the southeast of the country. These two groups have genetic affinities to northern South American countries and the Caribbean, respectively. This is consistent with what has been reported for other genetic markers such as mitochondrial DNA and allele frequencies at the insecticide resistance gene, kdr.

Conclusions/Significance

We conclude that the genetic patterns in present day populations of Ae. aegypti in Brazil are more consistent with a complete eradication of the species in the recent past followed by re-colonization, rather than the alternative possibility of expansion from residual pockets of refugia. At least two colonizations are likely to have taken place, one from northern South American countries (e.g., Venezuela) that founded the northwestern group, and one from the Caribbean that founded the southeastern group. The proposed source areas were never declared free of Ae. aegypti.     

Rewinding the process of mammalian extinction

With only three living individuals left on this planet, the northern white rhinoceros (Ceratotherium simum cottoni) could be considered doomed for extinction. It might still be possible, however, to rescue the (sub)species by combining novel stem cell and assisted reproductive technologies. To discuss the various practical options available to us, we convened a multidisciplinary meeting under the name “Conservation by Cellular Technologies.” The outcome of this meeting and the proposed road map that, if successfully implemented, would ultimately lead to a self‐sustaining population of an extremely endangered species are outlined here. The ideas discussed here, while centered on the northern white rhinoceros, are equally applicable, after proper adjustments, to other mammals on the brink of extinction. Through implementation of these ideas we hope to establish the foundation for reversal of some of the effects of what has been termed the sixth mass extinction event in the history of Earth, and the first anthropogenic one. Zoo Biol. 35:280–292, 2016. © 2016 The Authors. Zoo Biology published by Wiley Periodicals, Inc.     

Phylogeography and Population Structure of Glossina fuscipes fuscipes in Uganda: Implications for Control of Tsetse

Background

Glossina fuscipes fuscipes, a riverine species of tsetse, is the main vector of both human and animal trypanosomiasis in Uganda. Successful implementation of vector control will require establishing an appropriate geographical scale for these activities. Population genetics can help to resolve this issue by characterizing the extent of linkage among apparently isolated groups of tsetse.

Methodology/Principal Findings

We conducted genetic analyses on mitochondrial and microsatellite data accumulated from approximately 1000 individual tsetse captured in Uganda and neighboring regions of Kenya and Sudan. Phylogeographic analyses suggested that the largest scale genetic structure in G. f. fuscipes arose from an historical event that divided two divergent mitochondrial lineages. These lineages are currently partitioned to northern and southern Uganda and co-occur only in a narrow zone of contact extending across central Uganda. Bayesian assignment tests, which provided evidence for admixture between northern and southern flies at the zone of contact and evidence for northerly gene flow across the zone of contact, indicated that this structure may be impermanent. On the other hand, microsatellite structure within the southern lineage indicated that gene flow is currently limited between populations in western and southeastern Uganda. Within regions, the average F_ST between populations separated by less than 100 km was less than ∼0.1. Significant tests of isolation by distance suggested that gene flow is ongoing between neighboring populations and that island populations are not uniformly more isolated than mainland populations.

Conclusions/Significance

Despite the presence of population structure arising from historical colonization events, our results have revealed strong signals of current gene flow within regions that should be accounted for when planning tsetse control in Uganda. Populations in southeastern Uganda appeared to receive little gene flow from populations in western or northern Uganda, supporting the feasibility of area wide control in the Lake Victoria region by the Pan African Tsetse and Trypanosomiasis Eradication Campaign.     

The role of short-tongued insects in floral variation across the range of a style-dimorphic plant

Background and Aims

Heterostyly and related style polymorphisms are suitable model systems to evaluate the importance of functional pollinators in the maintenance of population variability. In Narcissus papyraceus different functional pollinators, incompatibility system and flower morphology have been proposed to influence the maintenance of polymorphism through their effect on disassortative mating. Here a test is done to find out if the visitation rate of long- versus short-tongued pollinators correlates with the morph ratio and if the latter is related to other flower traits of the species across its main geographic range.

Methods

Floral traits from 34 populations in the south-west of the Iberian Peninsula and in north-west Africa were measured, perianth variation was described and a comparison was made of allometric relationships between sex organs and floral tube. Correlations between pollinator guilds, stigma–anther separation of reciprocal morphs (our proxy for disassortative mating) and morph-ratio variation were analysed. Finally, the incompatibility system of the species in the northern and southern borders of its distribution are described.

Key Results

Flowers from southern populations were significantly larger than flowers from centre and northern populations. The abundance of short-styled plants decreased gradually with increasing distance from the core region (the Strait of Gibraltar), with these disappearing only in the northern range. Although there was a significant difference in stigma–anther separation among populations, morph ratio was not associated with reciprocity or floral tube length. Long-style morph frequency increased with short-tongued pollinator visitation rate. Populations from both edges of the distribution range were self-incompatible and within- and between-morph compatible.

Conclusions

The style morph ratio changed gradually, whereas perianth trait variation showed abrupt changes with two morphotypes across the range. The positive relationship between the visitation rate of short-tongued pollinators and the decrease of the short-style morph supports our initial hypothesis. The results highlight the importance of different pollinators in determining the presence of style polymorphism.     

Interpopulation variation in contour feather structure is environmentally determined in great tits

Background

The plumage of birds is important for flying, insulation and social communication. Contour feathers cover most of the avian body and among other functions they provide a critical insulation layer against heat loss. Feather structure and composition are known to vary among individuals, which in turn determines variation in the insulation properties of the feather. However, the extent and the proximate mechanisms underlying this variation remain unexplored.

Methodology/Principal Findings

We analyzed contour feather structure from two different great tit populations adapted to different winter regimes, one northern population in Oulu (Finland) and one southern population in Lund (Sweden). Great tits from the two populations differed significantly in feather structure. Birds from the northern population had a denser plumage but consisting of shorter feathers with a smaller proportion containing plumulaceous barbs, compared with conspecifics from the southern population. However, differences disappeared when birds originating from the two populations were raised and moulted in identical conditions in a common-garden experiment located in Oulu, under ad libitum nutritional conditions. All birds raised in the aviaries, including adult foster parents moulting in the same captive conditions, developed a similar feather structure. These feathers were different from that of wild birds in Oulu but similar to wild birds in Lund, the latter moulting in more benign conditions than those of Oulu.

Conclusions/Significance

Wild populations exposed to different conditions develop contour feather differences either due to plastic responses or constraints. Environmental conditions, such as nutrient availability during feather growth play a crucial role in determining such differences in plumage structure among populations.     

Tracing the Tiger: Population Genetics Provides Valuable Insights into the Aedes (Stegomyia) albopictus Invasion of the Australasian Region

Background

The range of the Asian tiger mosquito Aedes albopictus is expanding globally, raising the threat of emerging and re-emerging arbovirus transmission risks including dengue and chikungunya. Its detection in Papua New Guinea''s (PNG) southern Fly River coastal region in 1988 and 1992 placed it 150 km from mainland Australia. However, it was not until 12 years later that it appeared on the Torres Strait Islands. We hypothesized that the extant PNG population expanded into the Torres Straits as an indirect effect of drought-proofing the southern Fly River coastal villages in response to El Nino-driven climate variability in the region (via the rollout of rainwater tanks and water storage containers).

Methodology/Principal Findings

Examination of the mosquito''s mitochondrial DNA cytochrome oxidase I (COI) sequences and 13 novel nuclear microsatellites revealed evidence of substantial intermixing between PNG''s southern Fly region and Torres Strait Island populations essentially compromising any island eradication attempts due to potential of reintroduction. However, two genetically distinct populations were identified in this region comprising the historically extant PNG populations and the exotic introduced population. Both COI sequence data and microsatellites showed the introduced population to have genetic affinities to populations from Timor Leste and Jakarta in the Indonesian region.

Conclusions/Significance

The Ae. albopictus invasion into the Australian region was not a range expansion out of PNG as suspected, but founded by other, genetically distinct population(s), with strong genetic affinities to populations sampled from the Indonesian region. We now suspect that the introduction of Ae. albopictus into the Australian region was driven by widespread illegal fishing activity originating from the Indonesian region during this period. Human sea traffic is apparently shuttling this mosquito between islands in the Torres Strait and the southern PNG mainland and this extensive movement may well compromise Ae. albopictus eradication attempts in this region.     

Kestrel-Prey Dynamic in a Mediterranean Region: The Effect of Generalist Predation and Climatic Factors

Background

Most hypotheses on population limitation of small mammals and their predators come from studies carried out in northern latitudes, mainly in boreal ecosystems. In such regions, many predators specialize on voles and predator-prey systems are simpler compared to southern ecosystems where predator communities are made up mostly of generalists and predator-prey systems are more complex. Determining food limitation in generalist predators is difficult due to their capacity to switch to alternative prey when the basic prey becomes scarce.

Methodology

We monitored the population density of a generalist raptor, the Eurasian kestrel Falco tinnunculus over 15 years in a mountainous Mediterranean area. In addition, we have recorded over 11 years the inter-annual variation in the abundance of two main prey species of kestrels, the common vole Microtus arvalis and the eyed lizard Lacerta lepida and a third species scarcely represented in kestrel diet, the great white-toothed shrew Crocidura russula. We estimated the per capita growth rate (PCGR) to analyse population dynamics of kestrel and predator species.

Principal Findings

Multimodel inference determined that the PCGR of kestrels was better explained by a model containing the population density of only one prey species (the common vole) than a model using a combination of the densities of the three prey species. The PCGR of voles was explained by kestrel abundance in combination with annual rainfall and mean annual temperature. In the case of shrews, growth rate was also affected by kestrel abundance and temperature. Finally, we did not find any correlation between kestrel and lizard abundances.

Significance

Our study showed for the first time vertebrate predator-prey relationships at southern latitudes and determined that only one prey species has the capacity to modulate population dynamics of generalist predators and reveals the importance of climatic factors in the dynamics of micromammal species and lizards in the Mediterranean region.