Interplay between metabolic rate and diet quality in the South American fox, Pseudalopex culpaeus

We studied the metabolic costs associated with the ingestion of peppertree fruits (Schinus molle) in the culpeo fox, Pseudalopex culpaeus, the second largest canid in South America. Throughout its range of distribution, this fox feeds on rodents and other small vertebrates, and also on peppertree fruits, which represent 98% of total fruits consumed in semiarid Chile. Peppertree contains a high diversity of phytochemicals. Foxes feeding on diets containing rats and peppertree fruits (mixed diets) exhibited a 98.9% increase in basal rate of metabolism when compared to rat-acclimated foxes. Thus, acute ingestion of chemically defended fruits has an energetic cost for the fox, reflected in higher values of basal metabolism. Increased metabolic rates may be associated with increased protein synthesis for detoxification and for tissue repair, including the production of biotransformation enzymes.     

Contrasting genetic trajectories of endangered and expanding red fox populations in the western U.S

As anthropogenic disturbances continue to drive habitat loss and range contractions, the maintenance of evolutionary processes will increasingly require targeting measures to the population level, even for common and widespread species. Doing so requires detailed knowledge of population genetic structure, both to identify populations of conservation need and value, as well as to evaluate suitability of potential donor populations. We conducted a range-wide analysis of the genetic structure of red foxes in the contiguous western U.S., including a federally endangered distinct population segment of the Sierra Nevada subspecies, with the objectives of contextualizing field observations of relative scarcity in the Pacific mountains and increasing abundance in the cold desert basins of the Intermountain West. Using 31 autosomal microsatellites, along with mitochondrial and Y-chromosome markers, we found that populations of the Pacific mountains were isolated from one another and genetically depauperate (e.g., estimated Ne range = 3–9). In contrast, red foxes in the Intermountain regions showed relatively high connectivity and genetic diversity. Although most Intermountain red foxes carried indigenous western matrilines (78%) and patrilines (85%), the presence of nonindigenous haplotypes at lower elevations indicated admixture with fur-farm foxes and possibly expanding midcontinent populations as well. Our findings suggest that some Pacific mountain populations could likely benefit from increased connectivity (i.e., genetic rescue) but that nonnative admixture makes expanding populations in the Intermountain basins a non-ideal source. However, our results also suggest contact between Pacific mountain and Intermountain basin populations is likely to increase regardless, warranting consideration of risks and benefits of proactive measures to mitigate against unwanted effects of Intermountain gene flow.Subject terms: Conservation biology, Population genetics     

Population persistence in a landscape context: the case of endangered arctic fox populations in Fennoscandia

Anthropogenic fragmentation of habitat and populations is recognized as one of the most important factors influencing loss of biodiversity. Since it is difficult to quantify demographic parameters in small populations, we need alternative methods to elucidate important factors affecting the viability of local populations. The Fennoscandian arctic fox inhabits a naturally fragmented alpine tundra environment, but historic anthropogenic impacts have further fragmented its distribution. After almost 80 yr of protection, the population remains critically endangered. Both intrinsic factors (related to the isolation and size of sub‐populations) and extrinsic factors (related to environmental conditions influencing patch quality and interspecific competition) have been proposed as explanations for the lack of population growth. To distinguish between these hypotheses, we conducted a spatially explicit analysis that compares areas where the species has persisted with areas where it has become locally extinct. We used characteristics of the fragments of alpine tundra habitat and individual arctic fox breeding dens (including both currently active dens and historically active dens) within the fragments to evaluate the importance of habitat characteristics and connectivity in explaining variation in persistence within a fragment. The number of reproductive events in a fragment was related to the size of the fragment, but not more than expected following a 1:1 relationship, suggesting little effect of fragment size on the relative number of reproductions. The likelihood of a den being used for breeding was positively associated with factors minimising interspecific competition as well as increasing within‐fragment connectivity. These results support the idea that the failure of Fennoscandian arctic fox to recover is caused by demographic factors that can be related to fine‐scale Allee or Allee‐like effects, as well as environmental influences related to increased competition and exclusion by red foxes.     

Genetic variation in natural and translocated populations of the endangered Delmarva fox squirrel (Sciurus niger cinereus)

The Delmarva fox squirrel, Sciurus niger cinereus, is a federally listed endangered subspecies whose range has been reduced by 90%. In an attempt to increase both population size and range, translocation sites were established beginning in the 1960's by moving squirrels from the natural range to sites outside the current range. Although translocations have served as the primary component of the DFS recovery program, there has been very little post-release examination of the genetics of the translocation sites. In this study, we developed ten microsatellite loci, screened the three polymorphic loci, and sequenced a 330 bp fragment of the mitochondrial control region in order to assess levels of genetic variation in natural and translocated regions of Delmarva fox squirrels and to compare them to Southeastern fox squirrels (S. n. niger). Although we found low levels of microsatellite polymorphism, there were no differences in heterozygosity between natural and translocated regions, or between Delmarva and Southeastern fox squirrels. We found high levels of polymorphism in the mitochondrial control region. Our patterns of haplotype diversity suggest incomplete lineage sorting of the two subspecies. In general, our data suggest that the current levels of genetic variation in the translocated sites are representative of those found in the natural population, and we encourage the continued use of translocations as a major component of Delmarva fox squirrel recovery. This revised version was published online in July 2006 with corrections to the Cover Date.     

Hematologic values of the endangered San Joaquin kit fox, Vulpes macrotis mutica

Between 1981 and 1982 blood samples were collected from 64 adult San Joaquin kit foxes, Vulpes macrotis mutica, in western Kern County, California. The goal of the study was to establish normal blood values for this endangered species, and to determine whether changes in them could be used to assess the possible effects of petroleum developments on these foxes. None of the values differed significantly between the sexes, or between foxes sampled in developed habitats compared with foxes sampled in undisturbed habitats. Mean values of Hb, MCH, MCHC, and WBC counts differed significantly between summer and winter. Average hematological characteristics were: RBC, 8.4 X 10(6)/microliter; Hb, 14.5 g/dl (summer), 15.6 g/dl (winter); PCV, 46.9%; MCV, 56.3 fl; MCH, 17.8 pg (summer), 18.4 pg (winter); MCHC, 31.2 g/dl (summer), 33.2 g/dl (winter); and WBC, 6,200/microliter (summer), 7,500/microliter (winter). Comparisons of hematological data for kit foxes, coyotes (Canis latrans), and wolves (Canis lupus) confirmed a previously published observation that within mammalian families RBC counts are correlated inversely with body weight, and that MCV is correlated directly with body weight.     

Possible contemporary evolution in an endangered species, the Santa Cruz Island fox

An ability to mount rapid evolutionary responses to environmental change may be necessary for species persistence in a human-dominated world. We present evidence of the possibility of such contemporary evolution in the anti-predator behaviour of the critically endangered Santa Cruz Island fox Urocyon littoralis . In 1994, golden eagles colonized Santa Cruz Island, CA and devastated the predator-naïve, endemic island fox population by 95% within 10 years. In 1992, just before the arrival of golden eagles, foxes showed substantial diurnal activity, but diurnal activity was 37.0% lower in 2003–2007, after golden eagle colonization; concurrently, overall activity declined and nocturnal activity increased. Moreover, on nearby Santa Catalina Island, where golden eagles were absent but where the fox population recently crashed due to a disease epidemic, remaining foxes were significantly more diurnally active than were those on Santa Cruz Island. The weight of evidence suggests that the change in activity pattern was a response to predation, not to low population density, and that this was probably a heritable rather than a learned behavioural trait. This behavioural change may allow for prolonged island fox persistence, but also potentially represents a loss of behavioural diversity in fox populations.     

Fission and fusion of Darwin's finches populations

This study addresses the causes and evolutionary consequences of introgressive hybridization in the sympatric species of Darwin's ground finches (Geospiza) on the small island of Daphne Major in the Galápagos archipelago. Hybridization occurs rarely (less than 2% of breeding pairs) but persistently across years, usually as a result of imprinting on the song of another species. Hybrids survive well under some ecological conditions, but not others. Hybrids mate according to song type. The resulting introgression increases phenotypic and genetic variation in the backcrossed populations. Effects of introgression on beak shape are determined by the underlying developmental genetic pathways. Introgressive hybridization has been widespread throughout the archipelago in the recent past, and may have been a persistent feature throughout the early history of the radiation, episodically affecting both the speed and direction of evolution. We discuss how fission through selection and fusion through introgression in contemporary Darwin's finch populations may be a reflection of processes occurring in other young radiations. We propose that introgression has the largest effect on the evolution of interbreeding species after they have diverged in morphology, but before the point is reached when genetic incompatibilities incur a severe fitness cost.     

Evolutionary history and genetic connectivity across highly fragmented populations of an endangered daisy

Conservation management can be aided by knowledge of genetic diversity and evolutionary history, so that ecological and evolutionary processes can be preserved. The Button Wrinklewort daisy (Rutidosis leptorrhynchoides) was a common component of grassy ecosystems in south-eastern Australia. It is now endangered due to extensive habitat loss and the impacts of livestock grazing, and is currently restricted to a few small populations in two regions >500 km apart, one in Victoria, the other in the Australian Capital Territory and nearby New South Wales (ACT/NSW). Using a genome-wide SNP dataset, we assessed patterns of genetic structure and genetic differentiation of 12 natural diploid populations. We estimated intrapopulation genetic diversity to scope sources for genetic management. Bayesian clustering and principal coordinate analyses showed strong population genetic differentiation between the two regions, and substantial substructure within ACT/NSW. A coalescent tree-building approach implemented in SNAPP indicated evolutionary divergence between the two distant regions. Among the populations screened, the last two known remaining Victorian populations had the highest genetic diversity, despite having among the lowest recent census sizes. A maximum likelihood population tree method implemented in TreeMix suggested little or no recent gene flow except potentially between very close neighbours. Populations that were more genetically distinctive had lower genetic diversity, suggesting that drift in isolation is likely driving population differentiation though loss of diversity, hence re-establishing gene flow among them is desirable. These results provide background knowledge for evidence-based conservation and support genetic rescue within and between regions to elevate genetic diversity and alleviate inbreeding.Subject terms: Ecological genetics, Population genetics     

Similarity of contemporary and historical gene flow among highly fragmented populations of an endangered rattlesnake

Populations of endangered taxa in recently fragmented habitats often show high levels of genetic structure, but the role that contemporary versus historical processes play in generating this pattern is unclear. The eastern massasauga rattlesnake (Sistrurus c. catenatus) is an endangered snake that presently occurs throughout central and eastern North America in a series of populations that are isolated because of habitat fragmentation and destruction. Here, we use data from 19 species-specific microsatellite DNA loci to assess the levels of genetic differentiation, genetic effective population size, and contemporary and historical levels of gene flow for 19 populations sampled across the range of this snake. Eastern massasaugas display high levels of genetic differentiation (overall θ(Fst) = 0.21) and a Bayesian clustering method indicates that each population represents a unique genetic cluster even at regional spatial scales. There is a twofold variation in genetically effective population sizes but little genetic evidence that populations have undergone recent or historical declines in size. Finally, both contemporary and historical migration rates among populations were low and similar in magnitude even for populations located <7 km apart. A test of alternate models of population history strongly favours a model of long-term drift-migration equilibrium over a recent isolation drift-only model. These results suggest that recent habitat fragmentation has had little effect on the genetic characteristics of these snakes, but rather that this species has historically existed in small isolated populations that may be resistant to the long-term negative effects of inbreeding.     

Ecological niche modeling re‐examined: A case study with the Darwin's fox

Many previous studies have attempted to assess ecological niche modeling performance using receiver operating characteristic (ROC) approaches, even though diverse problems with this metric have been pointed out in the literature. We explored different evaluation metrics based on independent testing data using the Darwin's Fox (Lycalopex fulvipes) as a detailed case in point. Six ecological niche models (ENMs; generalized linear models, boosted regression trees, Maxent, GARP, multivariable kernel density estimation, and NicheA) were explored and tested using six evaluation metrics (partial ROC, Akaike information criterion, omission rate, cumulative binomial probability), including two novel metrics to quantify model extrapolation versus interpolation (E‐space index I) and extent of extrapolation versus Jaccard similarity (E‐space index II). Different ENMs showed diverse and mixed performance, depending on the evaluation metric used. Because ENMs performed differently according to the evaluation metric employed, model selection should be based on the data available, assumptions necessary, and the particular research question. The typical ROC AUC evaluation approach should be discontinued when only presence data are available, and evaluations in environmental dimensions should be adopted as part of the toolkit of ENM researchers. Our results suggest that selecting Maxent ENM based solely on previous reports of its performance is a questionable practice. Instead, model comparisons, including diverse algorithms and parameterizations, should be the sine qua non for every study using ecological niche modeling. ENM evaluations should be developed using metrics that assess desired model characteristics instead of single measurement of fit between model and data. The metrics proposed herein that assess model performance in environmental space (i.e., E‐space indices I and II) may complement current methods for ENM evaluation.     

Non-random fitness variation in two populations of Darwin's finches

Darwinian fitness of an individual is measured by the number of recruits it contributes to the next generation. We studied variation in fitness among members of three cohorts of two species of Darwin's finches living on the Galipagos island of Daphne Major: the medium ground finch (Geospiza fortis) and cactus finch (Geospiza scandens). Individuals of both species live for up to 16 years. Variation in fitness was neither random nor heritable. Non-randomness arises as a result of a few individuals living for an exceptionally long time and breeding many times. For each cohort, the number of recruits per breeder is strongly predicted by the number of fledglings per breeder. In turn, the number of fledglings is strongly predicted by longevity of the breeder. These results suggest that the most important determinant of fitness is the ability of an individual to survive to breed in many years. Morphological traits affect this ability. Although morphological traits are heritable they do not change unidirectionally because they are selected in opposite directions, and in different combinations, under fluctuating environmental conditions. Non-random fitness variation in fluctuating populations implies much smaller genetically effective sizes than breeding population sizes.     

Viability analysis of endangered crayfish populations

The noble crayfish Astacus astacus L. is a threatened freshwater invertebrate. Many of the remaining populations are isolated and there is considerable concern that diseases and the increased frequency of flooding events may drive these remnant populations to extinction. We performed a population viability analysis for a typical isolated noble crayfish population. We quantified the extinction risk by the mean time to extinction within 1000 years for several scenarios (flooding events, restocking of adults). For a set of parameters derived from field estimates, we estimated the mean time to extinction to be 240 years. However, the median was only 80 years. Multiple sensitivity analysis by logistic regression revealed that spawning probability, juvenile and adult mortality were the important parameters for the survival of the population. The mean time to extinction decreased with increasing frequency of floodings. This is alarming, considering the magnitude of the effect and the expectation of an increasing number of floodings with global warming. Restocking, however, was found to have only a minor effect on the mean time to extinction. Overall, our simulations suggested that for the long-term and self-sustaining survival of the noble crayfish, particularly where they remain isolated, we have to improve the extent and quality of habitats. Nevertheless, additional measures are necessary, especially the removal of dispersal barriers to allow some exchange of individuals between populations. However, this also calls for a control of invasive crayfish species.     

Conspecific versus heterospecific gene exchange between populations of Darwin's finches

This study addresses the extent and consequences of gene exchange between populations of Darwin''s finches. Four species of ground finches (Geospiza) inhabit the small island of Daphne Major in the centre of the Galápagos archipelago. We undertook a study of microsatellite DNA variation at 16 loci in order to quantify gene flow within species owing to immigration and between species owing to hybridization. A combination of pedigrees of observed breeders and assignments of individuals to populations by the program Structure enabled us to determine the frequency of gene exchange and the island of origin of immigrants in some cases. The relatively large populations of Geospiza fortis and G. scandens receive conspecific immigrants at a rate of less than one per generation. They exchange genes more frequently by rare but repeated hybridization. Effects of heterospecific gene flow from hybridization are not counteracted by lower fitness of the offspring. As a result, the standing genetic variation of the two main resident populations on Daphne Major is enhanced to a greater extent by introgressive hybridization than through breeding with conspecific immigrants. Immigrant G. fuliginosa also breeds with G. fortis. Conspecific immigration was highest in the fourth species, G. magnirostris. This species is much larger than the other three and perhaps for this reason it has not bred with any of them. The source island of most immigrants is probably the neighbouring island of Santa Cruz. Evolutionary change may be inhibited in G. magnirostris by continuing gene flow, but enhanced in G. fortis and G. scandens by introgressive hybridization.     

Genetic evidence for the persistence of the critically endangered Sierra Nevada red fox in California

California is home to both the native state-threatened Sierra Nevada red fox (Vulpes vulpes necator), which historically inhabited high elevations of the Sierra Nevada and Cascade mountains, and to multiple low-elevation red fox populations thought to be of exotic origin. During the past few decades the lowland populations have dramatically expanded their distribution, and possibly moved into the historic range of the native high-elevation fox. To determine whether the native red fox persists in its historic range in California, we compared mitochondrial cytochrome-b haplotypes of the only currently-known high-elevation population (n = 9 individuals) to samples from 3 modern lowland populations (n = 35) and historic (1911–1941) high-elevation (n = 22) and lowland (n = 7) populations. We found no significant population differentiation among the modern and historic high-elevation populations (average pairwise F _ST = 0.06), but these populations differed substantially from all modern and historic lowland populations (average pairwise F _ST = 0.52). Among lowland populations, the historic and modern Sacramento Valley populations were not significantly differentiated from one another (F _ST = −0.06), but differed significantly from recently founded populations in the San Francisco Bay region and in southern California (average pairwise F _ST = 0.42). Analysis of molecular variance indicated that 3 population groupings (mountain, Sacramento Valley, and other lowland regions) explained 45% of molecular variance (F _CT = 0.45) whereas only 4.5% of the variance was partitioned among populations within these groupings (F _SC = 0.08). These findings provide strong evidence that the native Sierra Nevada red fox has persisted in northern California. However, all nine samples from this population had the same haplotype, suggesting that several historic haplotypes may have become lost. Unidentified barriers have apparently prevented gene flow from the Sacramento Valley population to other eastern or southern populations in California. Future studies involving nuclear markers are needed to assess the origin of the Sierra Nevada red fox and to quantify levels of nuclear gene flow.     

Elevational isolation of red fox populations in the Greater Yellowstone Ecosystem

Foxes in the Greater Yellowstone Ecosystem are reported to show high frequencies of blonde and gray coat colors. A survey of park sighting records showed that the frequency of the novel coat colors significantly increases at elevations greater than 2300 m, suggesting some form of elevational isolation. We evaluated the degree of genetic separation between the high-elevation foxes (>2300 m) and contiguous populations of foxes at mid-elevations (1600–2300m). Low-elevation (>1600 m) foxes from North Dakota, >1000 km straight line distance from our populations, were used as a control group. We genotyped 15 high-elevation, 15 mid-elevation, and 10 low-elevation foxes at 10 microsatellite loci each. Heterozygosity was significantly lower in both the high-elevation and mid-elevation populations compared to the low-elevation foxes. The genetic differentiation was significantly greater between the high-elevation and mid-elevation foxes than between the mid-elevation and low-elevation foxes. Similarly, estimates of R_ST and F_ST suggest less gene flow occurs between the contiguous high- and mid-elevation fox populations than between the mid- and low-elevation fox populations separated by > 1000 km. The assignment test further supports this hypothesis. Although further work is needed, we suggest that the high-elevation foxes are remnant populations from the Wisconsin glaciation and should be managed as a unique population.     

Patterns of differentiation among endangered pondberry populations

Pondberry, Lindera melissifolia, is an endangered and partially clonally reproducing shrub species found in isolated populations that inhabit seasonally wet depressions in forested areas of the lower Mississippi River alluvial valley and southeastern regions of the United States. With eleven microsatellite loci, we quantified population genetic differentiation and diversity among 450 genets in 10 locations distributed across pondberry’s range. We used estimates of F _st and Jost’s D _est to measure genetic differences between populations and between geographic regions. The largest pairwise regional difference was found between eastern and western regional population groups (F _st = 0.23, D _est = 0.67), with the northern-most population groups in each region exhibiting larger divergence from each other than the southern-most population groups. Genetic diversity was lowest in the Sand Pond Conservation Area (A _e = 1.9, H _e = 0.36), which was the northern-most pondberry population, and highest in the Francis Marion National Forest (A _e = 4.1, H _e = 0.69), although we identified only 17 genets in that admixed population. Following adjustments for estimated null allele frequencies, we identified heterozygote excess in four eastern populations and found no evidence for inbreeding in any population. The observed patterns of differentiation indicate a phylogeography that exhibits an Appalachian Mountain discontinuity coupled with northward migrations along the Southern Atlantic Coastal Plain and into the Mississippi Alluvial Plain. The genetic consequences of this proposed phylogeographical structure may affect selection of germplasm sources for population reestablishment programs across pondberry’s range.