Search

Next-Generation sequencing technology has enabled accurate insights into the diet of wildlife species. The protocols for faecal sample collection and extraction of DNA for diet analysis have been different from those focusing on focal species DNA, even in studies combining questions on a focal species and its diet. We designed an experiment to evaluate two protocols and select a single one that will generate both focal species (Asiatic wild ass, Equus hemionus, in Israel) and diet DNA, as an effective strategy to minimize time, efforts, and cost without hampering efficiency. In Protocol A, we swabbed the outer surface of faecal boluses and extracted DNA using a Stool Kit, while for Protocol B, we homogenised faecal matter from inside the bolus followed by extraction using a Powersoil Kit. Protocol A performed significantly better for most of the measures, which included, for the focal species, microsatellite amplification success and the quantity of the GAPDH gene; and for its diet, amplification success of a chloroplast barcoding marker, trnL, number of reads, and plant genus richness and diversity. However, there was no significant difference for the amplification success of a sex-linked marker and plant genus composition. Both protocols yielded results for the focal species and its diet, demonstrating that one single protocol can be used for both purposes; although, a pilot study is recommended to optimise the protocol for specific systems. This strategy may also be useful for studies combining focal species and their gut microbiome and parasitic load.

Translocated animals typically find themselves in a novel environment in which they must establish a home range in a manner that will maximize their fitness. We hypothesized that the initial establishment of a home range is followed by adjustments expressed as home range shifting, and occurs as familiarity with the landscape increases, until the home range is stabilized. We studied the process of home range shifting in 42 female Persian fallow deer (Dama mesopotamica) reintroduced into the Galilee, Israel over a period of 2–5 years. We used changes in the degree of home range overlap between consecutive years as an indicator of stabilization. We then compared how the mean percent cover of the key vegetation types (woodland, scrubland and open pastures) differed between the areas abandoned in the first year's home range and the areas added to the last year's home range relative to the first (using a weighted paired t-test). We also compared the distribution (using χ2 test of independence and Levene's test for homogeneity of variance) of %cover of the 3 vegetation types between the first and last year's home range. The average home range overlap increased over the 5 years following the first release. During the first-year post release, deer avoided open pastures and preferred woodland. In later years deer increase in the % open pastures (weighted t-test: p < 0.001) and decreased the % woodland cover (weighted t-test: p = 0.07) by abandoning areas with little open pasture and steeper terrain and moving into areas with more open pasture and moderate terrain. Variance of the cover types across individuals increased with time. We conclude that the home ranges of the reintroduced deer stabilized with time. The changes in vegetation and slope are driven by time-dependent changing needs reflecting a tradeoff between safety (refuge) and foraging. Our findings suggest that using the initially established home range to determine species preferences can create a misleading picture of what the optimal home range of the species really is. Individual variation in term of preferences can take a few years to be expressed due to the initial high-risk perceived by individuals in a novel environment.

SummaryPlant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts.We generated an A. hierochuntica reference transcriptome and pinpointed extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat and identifying positively selected genes in A. hierochuntica.The two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat “stress-ready” state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype.We suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle.

Reintroductions are a powerful tool for the recovery of endangered species. However, their long‐term success is strongly influenced by the genetic diversity of the reintroduced population. The chances of population persistence can be improved by enhancing the population's adaptive ability through the mixing of individuals from different sources. However, where source populations are too diverse the reintroduced population could also suffer from outbreeding depression or unsuccessful admixture due to behavioural or genetic barriers. For the reintroduction of Asiatic wild ass Equus hemionus ssp. in Israel, a breeding core was created from individuals of two different subspecies (E. h. onager & E. h. kulan). Today the population comprises approximately 300 individuals and displays no signs of outbreeding depression. The aim of this study was a population genomic evaluation of this conservation reintroduction protocol. We used maximum likelihood methods and genetic clustering analyses to investigate subspecies admixture and test for spatial autocorrelation based on subspecies ancestry. Further, we analysed heterozygosity and effective population sizes in the breeding core prior to release and the current wild population. We discovered high levels of subspecies admixture in the breeding core and wild population, consistent with a significant heterozygote excess in the breeding core. Furthermore, we found no signs of spatial autocorrelation associated with subspecies ancestry in the wild population. Inbreeding and variance effective population size estimates were low. Our results indicate no genetic or behavioural barriers to admixture between the subspecies and suggest that their hybridization has led to greater genetic diversity in the reintroduced population. The study provides rare empirical evidence of the successful application of subspecies hybridization in a reintroduction. It supports use of intraspecific hybridization as a tool to increase genetic diversity in conservation translocations.

The long-term success of species reintroductions is strongly dependent on the availability of large areas of suitable habitat and the genetic make-up of the population. If available habitat is poorly connected this can hinder gene flow and lead to genetic fragmentation of the population, potentially increasing its extinction risk. We employed a conservation genomics approach in which we combined analyses of genetic structure with testing for potential landscape effects on habitat selection and gene flow in reintroduced Asiatic wild ass Equus hemionus ssp. in the Israeli Negev desert. Genetic structure and pairwise relatedness were first investigated followed by examination of landscape effects on individual habitat selection using records of GPS collared individuals. We then built habitat resistance surfaces and used electrical circuit theory to test for landscape effects on genetic relatedness. We detected weak genetic structuring, yet low spatial coherence among individuals from the same genetic cluster. Landscape variables had a significant impact on individual habitat selection, with wild ass avoiding steep slopes and habitats of low suitability as predicted by a species distribution model. However, the landscape genetic analysis revealed no effect of habitat resistance on genetic relatedness. These results suggest that gene flow in the reintroduced population is not impacted by landscape resistance. Indeed, the high mobility of the species may increase its resistance to the genetic effects of habitat fragmentation, at least over a small number of generations. We discuss other potential causes for the observed genetic structure including a behavioural effect. Our study highlights the importance of understanding species-habitat interactions for the long-term success of reintroductions.

In a reintroduced population, the social behavior of the species can strongly affect the long-term viability of the population through its effects on movement, information flow, disease spread and the population's genetic variability. Therefore, information on the social behavior of a reintroduced population can contribute to conservation practices; however, its importance is often underestimated. The initial phase of the Asiatic wild ass's (Equus hemionus) reintroduction in Israel has been considered a success, and the population is currently estimated at more than 250 individuals. However, the current social structure of the population remained unknown. We aimed to study this important population trait and to provide helpful information for efficient conservation and management protocols. The study was based on direct observations that were conducted over four consecutive years, and on the analyses of groups' composition and female groups' stability. Female groups accompanied by males constituted only 5% of the total 659 observations, males were observed to be mainly solitary or in groups of various sizes, and females were organized in non-stable groups, indicating that the reintroduced population exhibits a fission-fusion social structure. Identifying the social structure for the species in the expanding Negev population of the Asiatic wild ass can assist in implementing future reintroductions and can contribute to effective management decisions aimed at protecting the species.

Effective population size, a central concept in conservation biology, is now routinely estimated from genetic surveys and can also be theoretically predicted from demographic, life-history, and mating-system data. By evaluating the consistency of theoretical predictions with empirically estimated effective size, insights can be gained regarding life-history characteristics and the relative impact of different life-history traits on genetic drift. These insights can be used to design and inform management strategies aimed at increasing effective population size. We demonstrated this approach by addressing the conservation of a reintroduced population of Asiatic wild ass (Equus hemionus). We estimated the variance effective size (Nev ) from genetic data (N ev =24.3) and formulated predictions for the impacts on Nev of demography, polygyny, female variance in lifetime reproductive success (RS), and heritability of female RS. By contrasting the genetic estimation with theoretical predictions, we found that polygyny was the strongest factor affecting genetic drift because only when accounting for polygyny were predictions consistent with the genetically measured Nev . The comparison of effective-size estimation and predictions indicated that 10.6% of the males mated per generation when heritability of female RS was unaccounted for (polygyny responsible for 81% decrease in Nev ) and 19.5% mated when female RS was accounted for (polygyny responsible for 67% decrease in Nev ). Heritability of female RS also affected Nev ; hf2=0.91 (heritability responsible for 41% decrease in Nev ). The low effective size is of concern, and we suggest that management actions focus on factors identified as strongly affecting Nev, namely, increasing the availability of artificial water sources to increase number of dominant males contributing to the gene pool. This approach, evaluating life-history hypotheses in light of their impact on effective population size, and contrasting predictions with genetic measurements, is a general, applicable strategy that can be used to inform conservation practice.

The Near Eastern fire salamander, Salamandra infraimmaculata , is considered an endangered species in Israel and is near-threatened regionally. For 25 years, our laboratory has sought ethical sampling methods to protect individuals and populations of Salamandra . To “mark” individuals for estimating dispersal and population size, we use non-invasive individual-specific markings from photographs of larvae and adults. We demonstrated through mesocosm experiments (which are less mortality-driven than in nature) that exotic Gambusia affinis have extreme negative mortality effects on Salamandra larvae. From a compassionate conservation aspect, G. affinis should not be killed and placed in habitats where amphibians are not in danger and mosquitoes can be controlled. We identified breeding-site characteristics demonstrating that permanent breeding sites support larger adult populations than temporary breeding sites. For population genetics studies, we take minimal sized tail tips from adults (which have no adverse effects) for microsatellite data. For gene expression studies, rather than sacrifice entire bodies, we demonstrated that by taking only small larval tail tips, we could follow gene expression. We additionally demonstrated that tail tip removal does not affect survival, time to or size at metamorphosis. We documented high road kill rates at a specific breeding site. To prevent potential disease spread, we sterilize boots and sampling gear. We use results for implementing or recommending conservation of individuals and populations – e.g., identifying: movement corridors for breeding site dispersal; roadkill hotspots for under-road tunnels; suitable habitat for pool construction for more effective conservation; utilizing population genetics for recommending management units; information on demography and genetic diversity to identify hotspots for conservation; removal of Gambusia for amphibian protection.

Conservation biological control (CBC) seeks to minimize the deleterious effects of agricultural pests by enhancing the efficiency of natural enemies. Despite the documented potential of insectivorous bats to consume pests, many synanthropic bat species are still underappreciated as beneficial species. We investigated the diet of Kuhl's pipistrelle (Pipistrellus kuhlii), a common synanthropic insectivorous bat that forages in urban and agricultural areas, to determine whether it may function as a natural enemy in CBC. Faecal samples of P. kuhlii were collected throughout the cotton-growing season from five roost sites near cotton fields located in a Mediterranean agroecosystem, Israel, and analyzed using DNA metabarcoding. Additionally, data on estimated abundance of major cotton pests were collected. We found that the diet of P. kuhlii significantly varied according to sites and dates and comprised 27 species of agricultural pests that were found in 77.2% of the samples, including pests of key economic concern. The dominant prey was the widespread cotton pest, the pink bollworm, Pectinophora gossypiella, found in 31% of the samples and in all the roosts. Pink bollworm abundance was positively correlated with its occurrence in the bat diet. Furthermore, the bats' dietary breadth narrowed, while temporal dietary overlap increased, in relation to increasing frequencies of pink bollworms in the diet. This suggests that P. kuhlii exploits pink bollworm irruptions by opportunistic feeding. We suggest that synanthropic bats provide important pest suppression services, may function as CBC agents of cotton pests and potentially contribute to suppress additional deleterious arthropods found in their diet in high frequencies.

We evaluated a 20-yr-old spatially explicit model (SEM) that predicted the spatial expansion of reintroduced Persian fallow deer in northern Israel. Using the current distribution of the deer and based on multi-model inference we assessed the accuracy of the SEM's prediction and what other factors affected the population's current distribution. If the SEM's projection was still valid, the leading model in the multi-model inference would include only the SEM's projection as an explanatory variable with a good fit. Different leading models would reveal key variables overlooked when the SEM was constructed or changes in the landscape unforeseen at the time, thus assisting adaptive management and decision-making. We assessed deer presence from camera trap encounter counts analyzed using N-mixture models. Models included various combinations of seven predictors: the 20-yr predictions of an SEM developed during the initial phases of the reintroduction, three key landscape characteristics on which the SEM was originally based but updated to reflect current conditions, distance from the release site, elevation, and the distribution of gray wolves (a predator that was absent from the area when the SEM was developed). Competing models were ranked by Akaike information criterion (AIC). Wolf distribution was the key predictor explaining the current deer distribution, appearing in all three leading models (∆AIC

Species distribution models (SDMs) are often limited by the use of coarse-resolution environmental variables and by the number of observations required for their calibration. This is particularly true in the case of elusive animals. Here, we developed a SDM by combining three elements: a database of explanatory variables, mapped at a fine resolution; a systematic sampling scheme; and an intensive survey of indirect observations. Using MaxEnt, we developed the SDM for the population of the Asiatic wild ass (Equus hemionus), a rare and elusive species, at three spatial scales: 10, 100, and 1000 m per pixel. We used indirect observations of feces mounds. We constructed 14 layers of explanatory variables, in five categories: water, topography, biotic conditions, climatic variables and anthropogenic variables. Woody vegetation cover and slopes were found to have the strongest effect on the distribution of wild ass and were included as the main predictors in the SDM. Model validation revealed that an intensive survey of feces mounds and high-resolution predictor layers resulted in a highly accurate and informative SDM. Fine-grain (10 and 100 m) SDMs can be utilized to: (1) characterize the variables influencing species distribution at high resolution and local scale, including anthropogenic effects and geomorphologic features; (2) detect potential population activity centers; (3) locate potential corridors of movement and possible isolated habitat patches. Such information may be useful for the conservation efforts of the Asiatic wild ass. This approach could be applied to other elusive species, particularly large mammals.

Genetic diversity and structure of populations at the edge of the species' spatial distribution are important for potential adaptation to environmental changes and consequently, for the long-term survival of the species. Here, we combined classical population genetic methods with newly developed network analyses to gain complementary insights into the genetic structure and diversity of Acacia tortilis, a keystone desert tree, at the northern edge of its global distribution, where the population is under threat from climatic, ecological, and anthropogenic changes. We sampled A. tortilis from 14 sites along the Dead Sea region and the Arava Valley in Israel and in Jordan. In addition, we obtained samples from Egypt and Sudan, the hypothesized origin of the species. Samples from all sites were genotyped using six polymorphic microsatellite loci.Our results indicate a significant genetic structure in A. tortilis along the Arava Valley. This was detected at different hierarchical levels-from the basic unit of the subpopulation, corresponding to groups of trees within ephemeral rivers (wadis), to groups of subpopulations (communities) that are genetically more connected relative to others. The latter structure mostly corresponds to the partition of the major drainage basins in the area. Network analyses, combined with classical methods, allowed for the identification of key A. tortilis subpopulations in this region, characterized by their relatively high level of genetic diversity and centrality in maintaining gene flow in the population. Characterizing such key subpopulations may enable conservation managers to focus their efforts on certain subpopulations that might be particularly important for the population's long-term persistence, thus contributing to species conservation within its peripheral range.

Noninvasive genetic methods enable the sampling of natural populations while minimizing detrimental effects on them. However, noninvasive methods are marred by DNA extraction and amplification difficulties that can be mitigated by improved methodology. Past studies have shown that noninvasive genetic protocols are system specific and highlighted the importance of pilot studies in the establishment of genetic monitoring programs. We tested, using a factorial design experiment, the effect of different methods for the collection-preservation and extraction on the amplification of DNA from Nubian ibex (Capra nubiana) fecal samples. We found collection-preservation using paper bags and extraction with QIAamp® fast DNA stool mini kit to significantly enhance success rate compared to collection-preservation in ethanol and extraction with QIAamp® DNA mini kit. Our results will contribute to the studies of natural populations of the vulnerable C. nubiana.

Effective population size, a central concept in conservation biology, is now routinely estimated from genetic surveys and can also be theoretically predicted from demographic, life-history, and mating-system data. By evaluating the consistency of theoretical predictions with empirically estimated effective size, insights can be gained regarding life-history characteristics and the relative impact of different life-history traits on genetic drift. These insights can be used to design and inform management strategies aimed at increasing effective population size. We demonstrated this approach by addressing the conservation of a reintroduced population of Asiatic wild ass (Equus hemionus). We estimated the variance effective size (N(ev)) from genetic data (N(ev) = 24.3) and formulated predictions for the impacts on N(ev) of demography, polygyny, female variance in lifetime reproductive success (RS), and heritability of female RS. By contrasting the genetic estimation with theoretical predictions, we found that polygyny was the strongest factor affecting genetic drift because only when accounting for polygyny were predictions consistent with the genetically measured N(ev). The comparison of effective-size estimation and predictions indicated that 10.6% of the males mated per generation when heritability of female RS was unaccounted for (polygyny responsible for 81% decrease in N(ev)) and 19.5% mated when female RS was accounted for (polygyny responsible for 67% decrease in N(ev)). Heritability of female RS also affected N(ev); [Formula: see text] (heritability responsible for 41% decrease in N(ev)). The low effective size is of concern, and we suggest that management actions focus on factors identified as strongly affecting N(ev), namely, increasing the availability of artificial water sources to increase number of dominant males contributing to the gene pool. This approach, evaluating life-history hypotheses in light of their impact on effective population size, and contrasting predictions with genetic measurements, is a general, applicable strategy that can be used to inform conservation practice.

Taxonomic over-splitting of extinct or endangered taxa, due to an incomplete knowledge of both skeletal morphological variability and the geographical ranges of past populations, continues to confuse the link between isolated extant populations and their ancestors. This is particularly problematic with the genusEquus.To more reliably determine the evolution and phylogeographic history of the endangered Asiatic wild ass, we studied the genetic diversity and inter-relationships of both extinct and extant populations over the last 100,000 years, including samples throughout its previous range from Western Europe to Southwest and East Asia. Using 229 bp of the mitochondrial hypervariable region, an approach which allowed the inclusion of information from extremely poorly preserved ancient samples, we classify all non-African wild asses into nine clades that show a clear phylogeographic structure revealing their phylogenetic history. This study places the extinct European wild ass,E. hydruntinus, the phylogeny of which has been debated since the end of the 19thcentury, into its phylogenetic context within the Asiatic wild asses and reveals recent gene flow events between populations currently regarded as separate species. The phylogeographic organization of clades resulting from these efforts can be used not only to improve future taxonomic determination of a poorly characterized group of equids, but also to identify historic ranges, interbreeding events between various populations, and the impact of ancient climatic changes. In addition, appropriately placing extant relict populations into a broader phylogeographic and genetic context can better inform ongoing conservation strategies for this highly endangered species.

Movement of animals is a key process affecting population dynamics. Information on factors that affect pathway use is essential for identifying and protecting pathways, and important for maintaining connectivity among populations. We present an innovative, non-invasive, approach for predicting pathways of reintroduced Asiatic wild ass (Equus hemionus) in Israel, which is based on understanding the effects of landscape factors on pathways use. The approach includes: Predicting pathways, by employing a least cost pathway (LCP) GIS models based on several landscape factors, so as to efficiently direct a field survey and explore the wild ass’s general preferences of pathway types; Collecting empirical data by surveying the dung density of wild ass along each of the predicted pathways and using the data as an index of pathway use; Evaluating the predicted pathways against the empirical data collected, to estimate the general pathway preferences of the wild ass; and Developing and evaluating alternative generalized linear models, according to a priori hypotheses based on empirical data so as to quantify the effect of different landscape factors on pathway use. The analyses were conducted for the entire landscape, and then for two distinct landscape types, open landscape and landscape-barriers (mountain ridges), as subsets of the entire landscape. There were clear differences in the mean number of faeces counts between the LCPs, indicating that the wild ass prefers certain pathway types as a function of landscape features. We further found that the factors affecting E. hemionus pathway usage—vegetation; slopes; canyons; and 4-wheel drive routes—varied largely between the two major landscape types studied, demonstrating the importance of studying space use patterns at different landscape terrains. This information can be applicable to landscape planning measures that aim to enhance protection of the species. This approach provides a framework for studying animal space-use patterns of a variety of species, including elusive species, in a heterogeneous landscape.

The negative effects of habitat fragmentation and population isolation on population viability, genetic variability and structuring are well documented, and conservation plans failing to take into account spatial population structure and connectivity can be ineffectual. Of special concern are populations at the periphery of the species range that might show reduced genetic diversity, thus affecting their adaptive potential at environmental margins. We investigated genetic variability and differentiation of the globally near threatened and locally endangered fire salamander Salamandra infraimmaculata in northern Israel, an area that represents the periphery of this species' distribution range. Analyses of variability in 15 microsatellite loci from 20 sites revealed substantial population structuring, most of which was due to a strong subdivision between two regions separated by a heavily urbanized valley. In addition, levels of genetic variability within populations were lowest in the peripheral, southernmost populations. These results suggest that the conservation plans for this species should recognize the lower diversity and increased divergence in the peripheral regions, and take into account the observed spatial population structure when devising strategies and measures to ensure the species persistence.

Noninvasive genetic sampling has increasingly been used in ecological and conservation studies during the last decade. A major part of the noninvasive genetic literature is dedicated to the search for optimal protocols, by comparing different methods of collection, preservation and extraction of DNA from noninvasive materials. However, the lack of quantitative comparisons among these studies and the possibility that different methods are optimal for different systems make it difficult to decide which protocol to use. Moreover, most studies that have compared different methods focused on a single factor - collection, preservation or extraction - while there could be interactions between these factors. We designed a factorial experiment, as a pilot study, aimed at exploring the effect of several collection, preservation and extraction methods, and the interactions between them, on the quality and amplification success of DNA obtained from Asiatic wild ass (Equus hemionus) faeces in Israel. The amplification success rates of one mitochondrial DNA and four microsatellite markers differed substantially as a function of collection, preservation and extraction methods and their interactions. The most efficient combination for our system integrated the use of swabs as a collection method with preservation at -20 °C and with the Qiagen DNA Stool Kit with modifications as the DNA extraction method. The significant interaction found between the collection, preservation methods and the extraction methods reinforces the importance of conducting a factorial design experiment, rather than examining each factor separately, as a pilot study before initiating a full-scale noninvasive research project.

Releasing animals in more than one location may increase or decrease the probability of success of a reintroduction project, yet the question of how many release sites to use has received little attention. We used empirical data from the reintroduction program of the Persian fallow deer (Dama mesopotamica) (Galilee region in northern Israel) in an individual-based spatially explicit simulation model to assess the effects of releasing deer from multiple sites. We examined whether multiple release sites increase reintroduction success, and if so, whether the optimal number of sites for a given scenario can be determined and whether the outcome differs if animals are released alternately (i.e., the location of the release alternates yearly between sites) or consecutively (i.e., one release site is used for several years, then another is used, and so forth). We selected 8 potential release sites in addition to the original site and simulated the release of 180 individuals at a rate of 10 individuals per year in different combinations of the original site and 1-4 additional sites. In our model, releasing animals into the wild at multiple sites produced higher population growth and greater spatial expansion than releasing animals at only one site and a consecutive-release approach was superior to an alternate-release approach. We suggest that through the use of simulation modeling that is based on empirical data from previous releases, managers can make better-informed decisions regarding the use of multiple release sites and greatly improve the probability of reintroduction success.

Recent developments of automated methods for monitoring animal movement, e.g., global positioning systems (GPS) technology, yield high-resolution spatiotemporal data. To gain insights into the processes creating movement patterns, we present two new techniques for extracting information from these data on repeated visits to a particular site or patch (“recursions”). Identification of such patches and quantification of recursion pathways, when combined with patch-related ecological data, should contribute to our understanding of the habitat requirements of large herbivores, of factors governing their space-use patterns, and their interactions with the ecosystem. We begin by presenting output from a simple spatial model that simulates movements of large-herbivore groups based on minimal parameters: resource availability and rates of resource recovery after a local depletion. We then present the details of our new techniques of analyses (recursion analysis and circle analysis) and apply them to data generated by our model, as well as two sets of empirical data on movements of African buffalo (Syncerus caffer): the first collected in Klaserie Private Nature Reserve and the second in Kruger National Park, South Africa. Our recursion analyses of model outputs provide us with a basis for inferring aspects of the processes governing the production of buffalo recursion patterns, particularly the potential influence of resource recovery rate. Although the focus of our simulations was a comparison of movement patterns produced by different resource recovery rates, we conclude our paper with a comprehensive discussion of how recursion analyses can be used when appropriate ecological data are available to elucidate various factors influencing movement. Inter alia, these include the various limiting and preferred resources, parasites, and topographical and landscape factors.

Summary Successfully reintroduced species offer an exceptional opportunity to study the rules that govern population range expansion. We suggest that models that are based on knowledge of space-use patterns of reintroduced species can be used as projection models in the context of landscape planning. We illustrate this idea through a case study in which we used the population range-expansion patterns of Persian fallow deer ( Dama mesopotamica ) reintroduced in northern Israel, to explore potential impact of landscape development on species space-use patterns. In this approach we used an individual-based spatially explicit model that simulates species range expansion on a realistic landscape. We previously validated the model over the short-term (5 years) and apply it here for long-term, 50-year, projections using two scenarios: (1) current landscape conditions; and (2) the government's development plans for the region. Based on the long-term projections using the current landscape, we identified and mapped preferred deer habitat in which population activity centres were established. By applying the model to the government's development scenario for the region we examined the potential effect of landscape development on habitat availability and the future spatial distribution of the deer. By comparing the outputs from both scenarios threats to landscape connectivity could be identified. A map presenting key areas for landscape connectivity, important for the reintroduced species, can be used for long-term planning of local and regional development. Under the umbrella paradigm, protecting these areas may help to also conserve other species in the same ecosystem.

Securing natural habitats into which reintroduced species are released and those into which they will radiate is essential for the long-term establishment and viability of reintroduced populations. Consequently, the ability to project future range expansion is critical in reintroduction programs. We evaluated whether numerical growth and spatial expansion of a reintroduced population can be projected based on short-term monitoring of a newly reintroduced population, by studying a Persian fallow deer (Dama mesopotamica) population reintroduced to Israel. Over a five-year period, five to six adult (>1 yr) females were released semi-annually, for a total of 53 over the research period. We developed an individual-based spatially explicit model that projected spatial expansion of the female population in the wild. Model parameters were derived from empirical studies on movement patterns, home range establishment, behavioral adjustments, and demographic dynamics of reintroduced females during the first three years of the reintroduction project. Significant factors that were found to affect home-range establishment and that were incorporated into the model were: occurrence of Mediterranean woodland; occurrence of moderate terrain; roads; built-up areas; and presence of conspecifics. We tested the model by comparing its projections with field observations at the end of the five-year period since the project's onset, i.e., two years beyond the time on which parameters were based. The model was able to predict the direction and rate of expansion of the population and the main activity areas in the wild during the last two years. By using the model for long-term projection, future potential activity centers and landscape connectivity, important for securing the lands required for the Persian fallow deer future radiation, can be identified.

Because most reintroduced species are rare, data on their dynamics are scarce. Consequently, reintroduction programs often rely on data from other species or captive populations to project the performance of the reintroduced population in the wild. We compared the reproductive success and survival of a Persian fallow deer ( Dama mesopotamica) population reintroduced in Israel over the first 5 years of the project with the survival and reproduction parameters estimated while planning the reintroduction. In addition, we compared the actual growth of the wild population with the growth originally projected by a computer model in the original reintroduction program. We monitored 74 radio-collared individuals (57 females and 17 males) released semiannually 1996-2001. Survival during the first year after release was lower than later years (0.90 and 0.82 versus 0.95 and 0.88, for females and males, respectively). Such an impact was not anticipated in the original plan, but overall survival was higher than originally projected. As assumed in the reintroduction program, reproductive success improved significantly with time since release and overall, was higher than expected. The mean number of animals released annually was lower than planned. Overall, the growth of the reintroduced population was slower than projected, but the deviation was close to confidence limits and the pattern similar. After 5 years it appears that the original time frame of 8-10 years for project completion can be met or at worst will cause a 1-year delay. Over the short term of 5 years, projection models in reintroduction programs are useful tools for assessing the sustained use of the breeding core, depicting the dynamics of the population in the wild, providing a relatively accurate time frame for the successful completion of the project, and assessing project success.

Small populations are prone to loss of genetic variation and hence to a reduction in their evolutionary potential. Therefore, studying the mating system of small populations and its potential effects on genetic drift and genetic diversity is of high importance for their viability assessments. The traditional method for studying genetic mating systems is paternity analysis. Yet, as small populations are often rare and elusive, the genetic data required for paternity analysis are frequently unavailable. The endangered Asiatic wild ass (Equus hemionus), like all equids, displays a behaviourally polygynous mating system; however, the level of polygyny has never been measured genetically in wild equids. Combining noninvasive genetic data with stochastic modelling of shifts in allele frequencies, we developed an alternative approach to paternity analysis for studying the genetic mating system of the re-introduced Asiatic wild ass in the Negev Desert, Israel. We compared the shifts in allele frequencies (as a measure of genetic drift) that have occurred in the wild ass population since re-introduction onset to simulated scenarios under different proportions of mating males. We revealed a strongly polygynous mating system in which less than 25% of all males participate in the mating process each generation. This strongly polygynous mating system and its potential effect on the re-introduced population's genetic diversity could have significant consequences for the long-term persistence of the population in the Negev. The stochastic modelling approach and the use of allele-frequency shifts can be further applied to systems that are affected by genetic drift and for which genetic data are limited.

Releasing animals in more than one location may increase or decrease the probability of success of a reintroduction project, yet the question of how many release sites to use has received little attention. We used empirical data from the reintroduction program of the Persian fallow deer (Dama mesopotamica) (Galilee region in northern Israel) in an individual-based spatially explicit simulation model to assess the effects of releasing deer from multiple sites. We examined whether multiple release sites increase reintroduction success, and if so, whether the optimal number of sites for a given scenario can be determined and whether the outcome differs if animals are released alternately (i.e., the location of the release alternates yearly between sites) or consecutively (i.e., one release site is used for several years, then another is used, and so forth). We selected 8 potential release sites in addition to the original site and simulated the release of 180 individuals at a rate of 10 individuals per year in different combinations of the original site and 1-4 additional sites. In our model, releasing animals into the wild at multiple sites produced higher population growth and greater spatial expansion than releasing animals at only one site and a consecutive-release approach was superior to an alternate-release approach. We suggest that through the use of simulation modeling that is based on empirical data from previous releases, managers can make better-informed decisions regarding the use of multiple release sites and greatly improve the probability of reintroduction success.

Recursive movement—returns to previously visited areas—is a widespread phenomenon exhibited by a large range of species from bees and birds to primates and large felines, at different spatial scales. Nevertheless, the wide scope and generality of this phenomenon remain underestimated by the scientific community. This limited appreciation for the pervasiveness of recursive movement can be attributed to its study by parallel lines of research, with different methodologies and nomenclature, and almost no cross referencing among them. Among these lines of studies are traplining behavior in foraging ecology, path recursions in movement ecology and the ecology of fear in predator–prey studies. We synthesize these three lines of research, to underline the mechanisms driving these patterns and create a conceptual model for recursive movement behavior across species and spatio-temporal scales. The emergence and complexity of recursive movement patterns are determined by the rate of resource recovery, environmental...

The introduction of chronic, infectious diseases by colonizing populations (invasive or reintroduced) is a serious hazard in conservation biology, threatening the original host and other spillover species. Most research on spatial invasion of diseases has pertained to established host populations, either at steady state or fluctuating through time. Within a colonizing population, however, the spread of disease may be influenced by the expansion process of the population itself. Here we explore the simultaneous expansion of a colonizing population and a chronic, nonlethal disease introduced with it, describing basic patterns in homogeneous and structured landscapes and discussing implications for disease management. We describe expected outcomes of such introductions for three qualitatively distinct cases, depending on the relative velocities at which the population and epidemic expand. (1) If transmissibility is low the disease cannot be sustained, although it may first expand its range somewhat around the point of introduction. (2) If transmissibility is moderate but the wave- front velocity for the population, vp, is higher than that for the disease, vd, the disease wave front lags behind that of the population. (3) A highly transmissible disease, with vd . vp, will invade sufficiently rapidly to track the spread of the host. To test these elementary theoretical predictions, we simulated disease outbreaks in a spatially structured host population occupying a real landscape. We used a spatially explicit, individual-based model of Persian fallow deer (Dama mesopotamica) reintroduced in northern Israel, considering a hypothetical introduction of bovine tuberculosis. Basic patterns of disease expansion in this realistic setting were similar to our conceptual predictions for homogeneous landscapes. Landscape heterogeneity, however, induced the establishment of population activity centers and disease foci within them, leading to jagged wave fronts and causing local variation in the relative velocities at which the population and epidemic expanded. Based on predictions from simple theory and simulations of managed outbreaks, we suggest that the relative velocities at which the population and epidemic expand have important implications for the impact of different management strategies. Recognizing which of our three general cases best describes a particular outbreak will aid in planning an efficient strategy to contain the disease.

Allelic richness (number of alleles) is a measure of genetic diversity indicative of a population's long-term potential for adaptability and persistence. It is used less commonly than heterozygosity as a genetic diversity measure, partially because it is more mathematically difficult to take into account the stochastic process of genetic drift for allelic richness. This paper presents a stochastic model for the allelic richness of a newly founded population experiencing genetic drift and gene flow. The model follows the dynamics of alleles lost during the founder event and simulates the effect of gene flow on maintenance and recovery of allelic richness. The probability of an allele's presence in the population was identified as the relevant statistical property for a meaningful interpretation of allelic richness. A method is discussed that combines the probability of allele presence with a population's allele frequency spectrum to provide predictions for allele recovery. The model's analysis provides insights into the dynamics of allelic richness following a founder event, taking into account gene flow and the allele frequency spectrum. Furthermore, the model indicates that the “One Migrant per Generation” rule, a commonly used conservation guideline related to heterozygosity, may be inadequate for addressing preservation of diversity at the allelic level. This highlights the importance of distinguishing between heterozygosity and allelic richness as measures of genetic diversity, since focusing merely on the preservation of heterozygosity might not be enough to adequately preserve allelic richness, which is crucial for species persistence and evolution.

The space-use patterns of reintroduced animals are not well studied. We studied the space-use patterns of reintroduced Persian fallow deer (Dama mesopotamica) to assess how repeated releases from a single habituation enclosure and the presence of conspecifics influence the establishment of home range. We released 52 animals (33 were radiocollared). Movements of 27 radiocollared animals (5 M, 22 F) and 20 uncollared subadult males were studied by radiotracking and direct observation for ≤2 years from the time of reintroduction. Deer were released from the same habituation enclosure on 4 occasions over 20 months. During the first 3-6 months after release, deer usually occupied a limited range within 500 m from the enclosure. Later, animals moved to 1 of 3 areas 0.5-1.5 km from the enclosure, forming 3 subpopulations. The subpopulations consisted of animals from all releases and both sexes. Home ranges of animals from the first and second releases stabilized 8-10 months after release. Animals from the third release appeared to establish a stable home range almost immediately. Home-range size between September 1997 and April 1998 (7-8 mo) ranged from 86 to 365 ha, with home ranges of females not differing significantly from those of males. Later releases had little effect on the home-range stability of animals from previous releases. Our results suggest that over the short term, repeated releases from a single enclosure have no detrimental effects and actually may enhance the establishment of females from later releases. Movement patterns of released Persian fallow deer in this reintroduction indicate a slow, gradual movement away from the release site and establishment of a home range within a year.

The Persian fallow deer (Dama mesopotamica)-among the rarest deer species in the world-has been gradually reintroduced, using individuals from a captive-bred population, in northern Israel since September 1996. As of October 2000, >80 animals were in the wild population. We studied seasonal and circadian attributes of deer home ranges to assess the success of the reintroduction in terms of behavioral adjustments to the wild. We used radiotelemetry to determine locations and analyzed home ranges with the adaptive kernel method. We defined 3 seasons: fawning (Mar-Jun), rut (Jul-Oct), and winter (Nov-Feb). For females (n = 16), rut home ranges were significantly larger than winter home ranges (449 ± 45 ha [mean ± SE] vs. 384 ± 36 ha, P 1,15 = 0.013). During fawning, female home ranges were intermediate (424 ± 51 ha). Males (n = 5) increased their home ranges in rut season (820 ± 162 ha [mean ± SE], P < 0.012) and shifted their locations toward the release point. In winter, males significantly decreased their home ranges (584 ± 158 ha, P < 0.012), shifted their locations away from the release point, and almost no overlap of core areas was noticeable (1.8% overlap). In fawning (the time of antler casting and regrowth), males continued to shift away from the release point and decreased home ranges (358 ± 66 ha, P= 0.049) with almost no overlapping of core areas (0.06% overlap). No statistically significant differences were found between day home ranges (males [n = 5]: 621 ± 220 ha [mean ± SD], females [n = 16]: 402 ± 164) and night home ranges (males: 482 ± 145, females: 389 ± 183), although day core areas tended to be larger (in all males and 12 of 16 females). All documented aspects of seasonally, in female and male home ranges are in accordance with the annual reproduction cycle, and are related to seasonal food availability. These results, combined with previous works, suggest that so far, the reintroduced Persian fallow deer have adjusted well to living in the wild and that the chances of achieving A self-wild population are good. However, further research for extended period should verify these conclusions.