Sub-zero temperatures and large-scale weather patterns induce tooth damage in Icelandic arctic foxes.

Tooth damage in carnivores can reflect shifts in both diet and feeding habits, and in large carnivores, it is associated with increased bone consumption. Variation in tooth condition in Icelandic arctic foxes, a mesocarnivore, was recorded from 854 individual foxes spanning 29 years. We hypothesized that annual climatic variations, which can influence food abundance and accessibility, will influence tooth condition by causing dietary shifts toward less edible prey. We examined tooth condition in relation to four climatic predictors: mean annual winter temperature, indices of both the El Niño anomaly and North Atlantic subpolar gyre (SPG), and the number of rain-on-snow days (ROS). We found unequivocal evidence for a strong effect of annual climate on tooth condition. Teeth of Icelandic foxes were in better condition when winter temperatures were higher, when the SPG was more positive, and when the number of ROS was low. We also found a substantial subregional effect with foxes from northeastern Iceland having lower tooth damage than those from two western sites. Contradicting our original hypothesis that foxes from northeastern Iceland, where foxes are known to scavenge on large mammal remains (e.g., sheep and horses), would show the highest tooth damage, we suggest that western coastal sites exhibited greater tooth damage because cold winter temperatures lowered the availability of seabirds, causing a shift in diet toward abrasive marine subsidies (e.g., bivalves) and frozen beach wrack. Our study shows that monitoring tooth breakage and wear can be a useful tool for evaluating the impact of climate on carnivore populations and that climate change may influence the condition and fitness of carnivores in complex and potentially conflicting ways.

Male rock hyraxes that maintain an isochronous song rhythm achieve higher reproductive success.

Rhythmic stability (nonrandom temporal structure) is required for many neural and physiological functions, whereas rhythmic irregularities can indicate genetic or developmental deficiencies. Therefore, rhythmic courtship or contest signals are widespread in nature as honest advertisement displays. Examination of bird songs revealed the pervasiveness of categorical rhythmic patterns that can be described as small integer ratios between sequential inter-call intervals. As similar rhythmic profiles are prevalent in human music, it was suggested that a shared functionality could drive both animal songs and human musical rhythms, facilitating synchrony between signallers and enabling easy identification of performance errors. Here we examined whether the rhythmic structure and the rhythmic stability of vocal displays are related to reproductive success in male rock hyraxes (Procavia capensis), which presents an unusual case of a terrestrial singing mammal. We combined long-term parentage analysis of 13 male hyraxes (22 male/years) with an analysis of an audio library of 105 hyrax songs. Male annual reproductive success was determined by the number of offspring that survived to the age of 1 year. The frequency of singing events was used to determine the seasonal singing effort for each male. Songs were analysed for rhythmic structure, focusing on the presence of categorical rhythms and the contribution of rhythmic stability to annual reproductive success. We found that male hyraxes that sing more frequently tend to have more surviving offspring and that the rhythmic profile of hyrax songs is predominantly isochronous with sequential vocal element pairs nearly equally spaced. The ratio of isochronous vocal element transitions (on-integer) to element transitions that deviate from an isochronous pattern (off-integer) in hyrax songs is positively correlated with male reproductive success. Our findings support the notion that isochronous rhythmic stability can serve as an indication of quality in sexually selected signals and is not necessarily driven by the need for multiple caller synchronization. The relative scarcity of nonisochronous rhythmic categories in individually performed hyrax songs raises the question of whether such rhythmic categories could be a product of collective, coordinated signalling, while being selected against in individual performance.

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates).

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.

Social context mediates testosterone's effect on snort acoustics in male hyrax songs.

Testosterone affects physical and motivational states, both of which may strongly influence vocalization structure and acoustics. The loud complex calls (i.e., songs) of male rock hyraxes (Procavia capensis) are used as honest signals for advertising physical and social states. The snort, a low frequency, noisy element of the song, encodes information on the singer's age and social rank via harshness, as measured by jitter (i.e., acoustic frequency stability) and duration; suggesting that the snort concomitantly advertises both vocal stability and aggression. Our past findings revealed that testosterone levels are related to both vocal elements and social status of male hyraxes, suggesting that hormonal mechanisms mediate the motivation for aggressive and courtship behaviors. Here we examined whether long-term androgen levels are related to snort acoustics and song structure by comparing levels of testosterone in hair with acoustic and structural parameters. We found that songs performed by individuals with higher testosterone levels include more singing bouts and longer, smoother snorts, but only in those songs induced by external triggers. It is possible that hyraxes with higher levels of testosterone possess the ability to perform higher-quality singing, but only invest in situations of high social arousal and potential benefit. Surprisingly, in spontaneous songs, hyraxes with high testosterone were found to snort more harshly than low-testosterone males. The context dependent effects of high testosterone on snort acoustics suggest that the aggressive emotional arousal associated with testosterone is naturally reflected in the jittery hyrax snort, but that it can be masked by high-quality performance.

The plant component of an Acheulian diet at Gesher Benot Ya'aqov, Israel.

Diet is central for understanding hominin evolution, adaptation, and environmental exploitation, but Paleolithic plant remains are scarce. A unique macrobotanical assemblage of 55 food plant taxa from the Acheulian site of Gesher Benot Ya'aqov, Israel includes seeds, fruits, nuts, vegetables, and plants producing underground storage organs. The food plant remains were part of a diet that also included aquatic and terrestrial fauna. This diverse assemblage, 780,000 y old, reflects a varied plant diet, staple plant foods, environmental knowledge, seasonality, and the use of fire in food processing. It provides insight into the wide spectrum of the diet of mid-Pleistocene hominins, enhancing our understanding of their adaptation from the perspective of subsistence. Our results shed light on hominin abilities to adjust to new environments, facilitating population diffusion and colonization beyond Africa. We reconstruct the major vegetal foodstuffs, while considering the possibility of some detoxification by fire. The site, located in the Levantine Corridor through which several hominin waves dispersed out of Africa, provides a unique opportunity to study mid-Pleistocene vegetal diet and is crucial for understanding subsistence aspects of hominin dispersal and the transition from an African-based to a Eurasian diet.

Demographically-Based Evaluation of Genomic Regions under Selection in Domestic Dogs.

Controlling for background demographic effects is important for accurately identifying loci that have recently undergone positive selection. To date, the effects of demography have not yet been explicitly considered when identifying loci under selection during dog domestication. To investigate positive selection on the dog lineage early in the domestication, we examined patterns of polymorphism in six canid genomes that were previously used to infer a demographic model of dog domestication. Using an inferred demographic model, we computed false discovery rates (FDR) and identified 349 outlier regions consistent with positive selection at a low FDR. The signals in the top 100 regions were frequently centered on candidate genes related to brain function and behavior, including LHFPL3, CADM2, GRIK3, SH3GL2, MBP, PDE7B, NTAN1, and GLRA1. These regions contained significant enrichments in behavioral ontology categories. The 3rd top hit, CCRN4L, plays a major role in lipid metabolism, that is supported by additional metabolism related candidates revealed in our scan, including SCP2D1 and PDXC1. Comparing our method to an empirical outlier approach that does not directly account for demography, we found only modest overlaps between the two methods, with 60% of empirical outliers having no overlap with our demography-based outlier detection approach. Demography-aware approaches have lower-rates of false discovery. Our top candidates for selection, in addition to expanding the set of neurobehavioral candidate genes, include genes related to lipid metabolism, suggesting a dietary target of selection that was important during the period when proto-dogs hunted and fed alongside hunter-gatherers.

Living quarters of a living fossil-Uncovering the current distribution pattern of the rediscovered Hula painted frog (Latonia nigriventer) using environmental DNA.

One of the greatest challenges of effective conservation measures is the correct identification of sites where rare and elusive organisms reside. The recently rediscovered Hula painted frog (Latonia nigriventer) has not been seen for many decades and was therefore categorized extinct. Since its rediscovery in 2011, individuals from the critically endangered species have been found, with great effort, only in four restricted sites. We applied the environmental DNA (eDNA) approach to search for new populations of the Hula painted frog in suitable aquatic habitats. We further used the eDNA data to classify the landscape factors associated with the species distribution and to predict its suitable habitats. We sampled 52 aquatic sites in the Hula Valley during the spring of 2015 and 2016 and amplified the samples with a species-specific qPCR assay. DNA of the Hula painted frog was detected in 22 of the sites, all of which clustered within three main areas. A boosting classification model showed that soil type, vegetation cover and the current and former habitats are all key predictors of the frog's current distribution. Intriguingly, the habitat suitability models reveal a high affinity of the species to its long-lost habitat of the historical wetlands. Our findings encourage a series of informed searches for new populations of this threatened frog and provide guidance for future conservation management programmes. In the era of global conservation crisis of amphibians, developing the eDNA approach, a reliable detection method for many critically endangered and elusive amphibians, is of particular importance.

Genome sequencing highlights the dynamic early history of dogs.

To identify genetic changes underlying dog domestication and reconstruct their early evolutionary history, we generated high-quality genome sequences from three gray wolves, one from each of the three putative centers of dog domestication, two basal dog lineages (Basenji and Dingo) and a golden jackal as an outgroup. Analysis of these sequences supports a demographic model in which dogs and wolves diverged through a dynamic process involving population bottlenecks in both lineages and post-divergence gene flow. In dogs, the domestication bottleneck involved at least a 16-fold reduction in population size, a much more severe bottleneck than estimated previously. A sharp bottleneck in wolves occurred soon after their divergence from dogs, implying that the pool of diversity from which dogs arose was substantially larger than represented by modern wolf populations. We narrow the plausible range for the date of initial dog domestication to an interval spanning 11-16 thousand years ago, predating the rise of agriculture. In light of this finding, we expand upon previous work regarding the increase in copy number of the amylase gene (AMY2B) in dogs, which is believed to have aided digestion of starch in agricultural refuse. We find standing variation for amylase copy number variation in wolves and little or no copy number increase in the Dingo and Husky lineages. In conjunction with the estimated timing of dog origins, these results provide additional support to archaeological finds, suggesting the earliest dogs arose alongside hunter-gathers rather than agriculturists. Regarding the geographic origin of dogs, we find that, surprisingly, none of the extant wolf lineages from putative domestication centers is more closely related to dogs, and, instead, the sampled wolves form a sister monophyletic clade. This result, in combination with dog-wolf admixture during the process of domestication, suggests that a re-evaluation of past hypotheses regarding dog origins is necessary.

Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication.

Advances in genome technology have facilitated a new understanding of the historical and genetic processes crucial to rapid phenotypic evolution under domestication. To understand the process of dog diversification better, we conducted an extensive genome-wide survey of more than 48,000 single nucleotide polymorphisms in dogs and their wild progenitor, the grey wolf. Here we show that dog breeds share a higher proportion of multi-locus haplotypes unique to grey wolves from the Middle East, indicating that they are a dominant source of genetic diversity for dogs rather than wolves from east Asia, as suggested by mitochondrial DNA sequence data. Furthermore, we find a surprising correspondence between genetic and phenotypic/functional breed groupings but there are exceptions that suggest phenotypic diversification depended in part on the repeated crossing of individuals with novel phenotypes. Our results show that Middle Eastern wolves were a critical source of genome diversity, although interbreeding with local wolf populations clearly occurred elsewhere in the early history of specific lineages. More recently, the evolution of modern dog breeds seems to have been an iterative process that drew on a limited genetic toolkit to create remarkable phenotypic diversity.

Effects of parasite pressure on parasite mortality and reproductive output in a rodent-flea system: inferring host defense trade-offs.

Evaluating host resistance via parasite fitness helps place host-parasite relationships within evolutionary and ecological contexts; however, few studies consider both these processes simultaneously. We investigated how different levels of parasite pressure affect parasite mortality and reproductive success in relationship to host defense efforts, using the rodent Gerbillus nanus and the flea Xenopsylla conformis as a host-parasite system. Fifteen immune-naïve male rodents were infested with 20, 50, or 100 fleas for four weeks. During this time number of new imagoes produced per adult flea (our flea reproductive output metric), flea mortality, and change in circulating anti-flea immunoglobulin G (our measure of adaptive immune defense) were monitored. Three hypotheses guided this work: (1) increasing parasite pressure would heighten host defenses; (2) parasite mortality would increase and parasite reproductive output would decrease with increasing investment in host defense; and (3) hosts under high parasite pressure could invest in behavioral and/or immune responses. We predicted that at high infestation levels (a) parasite mortality would increase; (b) flea reproductive output per individual would decrease; and (c) host circulating anti-flea antibody levels would increase. The hypotheses were partially supported. Flea mortality significantly increased and flea reproductive output significantly decreased as flea pressure increased. Host adaptive immune defense did not significantly change with increasing flea pressure. Therefore, we inferred that investment in host behavioral defense, either alone or in combination with density-dependent effects, may be more efficient at increasing flea mortality and decreasing flea reproductive output than antibody production during initial infestation in this system.

Molecular assessment of the phylogeny and biogeography of a recently diversified endemic group of South American canids (Mammalia: Carnivora: Canidae).

To investigate the evolution and biogeography of an endemic group of South American foxes, we examined mitochondrial DNA control region sequences for 118 individuals belonging to all six extant species of the genus Lycalopex. Phylogenetic and molecular dating analyses supported the inference that this genus has undergone a very recent and rapid radiation, stemming from a common ancestor that lived ca. 1 million years ago. The Brazilian endemic L. vetulus was supported as the most basal species in this genus, whereas the most internal group is comprised by the recently diverged (ca. 350,000 years ago) Andean/Patagonian species L. griseus and L. culpaeus. We discuss the inferred phylogenetic relationships and divergence times in the context of the current geographic distributions of these species, and the likely effects of Pleistocene climatic changes on the biogeography of this group. Furthermore, a remarkable finding was the identification of multiple individuals classified as L. gymnocercus bearing mtDNA haplotypes clearly belonging to L. griseus, sampled in regions where the latter is not known to occur. At a minimum, this result implies the need to clarify the present-day geographic distribution of each of these fox species, while it may also indicate an ongoing hybridization process between them. Future testing of this hypothesis with in-depth analyses of these populations is thus a priority for understanding the history, evolutionary dynamics and present-day composition of this endemic Neotropical genus.

A genome-wide perspective on the evolutionary history of enigmatic wolf-like canids.

High-throughput genotyping technologies developed for model species can potentially increase the resolution of demographic history and ancestry in wild relatives. We use a SNP genotyping microarray developed for the domestic dog to assay variation in over 48K loci in wolf-like species worldwide. Despite the high mobility of these large carnivores, we find distinct hierarchical population units within gray wolves and coyotes that correspond with geographic and ecologic differences among populations. Further, we test controversial theories about the ancestry of the Great Lakes wolf and red wolf using an analysis of haplotype blocks across all 38 canid autosomes. We find that these enigmatic canids are highly admixed varieties derived from gray wolves and coyotes, respectively. This divergent genomic history suggests that they do not have a shared recent ancestry as proposed by previous researchers. Interspecific hybridization, as well as the process of evolutionary divergence, may be responsible for the observed phenotypic distinction of both forms. Such admixture complicates decisions regarding endangered species restoration and protection.

The phylogeography of Middle Eastern tree frogs in Israel.

Western Palearctic treefrogs of the genus Hyla provide an example of a morphologically and ecologically cryptic group. Up to three distinct Hyla species have been proposed as resident in Israel and this number has consistently been subject to taxonomical debates. Here, we analyzed 16S rRNA and COI gene fragments of 658 individuals sampled at 47 pools in nine regions across Israel and the West Bank in order to resolve the taxonomic status of Hyla frogs. We generated both Bayesian and Maximum Likelihood phylogenies, and constructed time-calibrated trees to provide an evolutionary and historical context of sequence variations. We further applied SAMOVA as well as Monmonier's maximum-difference algorithm to study the genetic structure among populations and to identify potential zones acting as barriers to gene flow across locations. Our results revealed two distinct haplogroups for each gene fragment, with 95% CI divergence times dated from 8.9-17.1 Mya (16S) and 7.1-23.6 Mya (COI), respectively. SAMOVA and barrier analyses partitioned the populations into three groups. Our results highlight that, while there are probably only two Hyla species in Israel, one population of one of the species might qualify as a separate evolutionarily significant unit. Our findings elucidate the taxonomic status of Hyla frogs in Israel and provide the basis for determining appropriate management and conservation priorities.

Genetic and behavioural factors affecting interpopulation colour pattern variation in two congeneric chameleon species.

We conducted a study on interpopulation variation of colour patterns in two congeneric chameleon species, which have an analogous life history. Both species are able to rapidly change colour pattern, and their context-dependent colour patterns often vary across a wide geographical range. Specifically, we tested four hypotheses that can explain the observed interpopulation variation of colour patterns by a series of behavioural field trials where the colour patterns of individuals were recorded and later analysed by a deep neural network algorithm. We used redundancy analysis to relate genetic, spectral and behavioural predictors to interpopulation colour pattern distance. Our results showed that both isolation by distance (IBD) and alternative mating tactics were significant predictors for interpopulation colour pattern variation in Chamaeleo chamaeleon males. By contrast, in Chamaeleo dilepis, the interpopulation colour pattern variation was largely explained by IBD, and evidence for alternative mating tactics was absent. In both chameleon species, the environmental colours showed no evidence of influencing chameleon interpopulation colour pattern variation, regardless of sex or behavioural context. This contrasting finding suggests that interpopulation context-dependent colour pattern variations in each species are maintained under a different set of selective pressures or circumstances.

Macroevolution of life-history traits in passerine birds: adaptation and phylogenetic inertia.

We used a recent passerine phylogeny and comparative method to evaluate the macroevolution of body and egg mass, incubation and fledging periods, time to independence and time with parents of the main passerine lineages. We hypothesised that passerine reproductive traits are affected by adaptation to both past and present environmental factors and phenotypic attributes such as body mass. Our results suggest that the evolution of body and egg mass, time to independence, incubation and fledging times are affected by strong phylogenetic inertia and that these breeding traits are all affected by body mass. Time with parents, where major lineages exhibit their own fixed optima and body mass does not have an effect, and clutch size which is affected by body mass and additionally by climate regimes, do not exhibit any phylogenetic inertia.

Testing the morphological constraint hypothesis of tail length in the sexually dimorphic Cerastes vipera and new perspectives.

The morphological constraint hypothesis (MCH) states that, in snakes, males typically have relatively longer tails than females to accommodate the hemipenes and retractor muscles. To date, most studies testing the MCH have been interspecific and results have been equivocal. We tested the MCH intraspecifically on Cerastes vipera, a species with a relatively short tail and suitable for testing the MCH. The relative tail length and length of the hemipenes pocket in Cerastes vipera were measured in preserved museum-maintained males (n = 35) and in free-ranging males and females (n = 277). Males exhibited relatively longer tails than females, which was explained fully by the length of the hemipenes pocket. The relatively short tail of C. vipera presents a constraint to the reproductive structures in males, as the length of the hemipenes pocket occupies a greater proportion in shorter- than longer-tailed individuals. This is the first report presenting these intraspecific findings in support of the MCH. Whether these relations are widespread among snake families, within Viperidae, or specifically within C. vipera warrants further studies.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal.

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

The interaction between cortisol and testosterone predicts leadership within rock hyrax social groups.

Group movement leadership is associated with higher risks for those in the front. Leaders are the first to explore new areas and may be exposed to predation. Individual differences in risk-taking behavior may be related to hormonal differences. In challenging circumstances, such as risk-taking leadership that may pose a cost to the leader, cortisol is secreted both to increase the likelihood of survival by restoring homeostasis, and to mediate cooperative behavior. Testosterone too has a well-established role in risk-taking behavior, and the dual-hormone hypothesis posits that the interaction of testosterone and cortisol can predict social behavior. Based on the dual-hormone hypothesis, we investigated here whether the interaction between testosterone and cortisol can predict risk-taking leadership behavior in wild rock hyraxes (Procavia capensis). We used proximity loggers, observations, and playback trials to characterize hyrax leaders in three different leadership contexts that varied in their risk levels. In support of the dual-hormone hypothesis, we found that cortisol and testosterone interactions predict leadership that involves risk. Across different circumstances that involved low or high levels of risk, testosterone was positively related to leadership, but only in individuals (both males and females) with low levels of cortisol. We also found an interaction between these hormone levels and age at the low-risk scenarios. We suggest that the close social interactions and affiliative behavior among hyrax females within small egalitarian groups may make female leadership less risky, and therefore less stressful, and allow female leaders to influence group activities.

Syntactic structure and geographical dialects in the songs of male rock hyraxes.

Few mammalian species produce vocalizations that are as richly structured as bird songs, and this greatly restricts the capacity for information transfer. Syntactically complex mammalian vocalizations have been previously studied only in primates, cetaceans and bats. We provide evidence of complex syntactic vocalizations in a small social mammal: the rock hyrax (Procavia capensis: Hyracoidea). We adopted three algorithms, commonly used in genetic sequence analysis and information theory, to examine the order of syllables in hyrax calls. Syntactic dialects exist, and the syntax of hyrax calls is significantly different between different regions in Israel. Call syntax difference is positively correlated to geographical distance over short distances. No correlation is found over long distances, which may reflect limited dispersal movement. These findings indicate that rich syntactic structure is more common in the vocalizations of mammalian taxa than previously thought and suggest the possibility of vocal production learning in the hyrax.

Non-breeding feather concentrations of testosterone, corticosterone and cortisol are associated with subsequent survival in wild house sparrows.

Potential mechanistic mediators of Darwinian fitness, such as stress hormones or sex hormones, have been the focus of many studies. An inverse relationship between fitness and stress or sex hormone concentrations has been widely assumed, although empirical evidence is scarce. Feathers gradually accumulate hormones during their growth and provide a novel way to measure hormone concentrations integrated over time. Using liquid chromatography-tandem mass spectrometry, we measured testosterone, corticosterone and cortisol in the feathers of house sparrows (Passer domesticus) in a wild population which is the subject of a long-term study. Although corticosterone is considered the dominant avian glucocorticoid, we unambiguously identified cortisol in feathers. In addition, we found that feathers grown during the post-nuptial moult in autumn contained testosterone, corticosterone and cortisol levels that were significantly higher in birds that subsequently died over the following winter than in birds that survived. Thus, feather steroids are candidate prospective biomarkers to predict the future survival of individuals in the wild.