Faecal cortisol metabolites to assess stress in wildlife: evaluation of a field method in free‐ranging chamois

Spontaneous Lesions of Endangered Geriatric Julia Creek Dunnarts (Sminthopsis douglasi, Archer 1979) with Emphasis in Reproductive Pathology

Julia Creek dunnarts are an endangered species of carnivorous marsupials and the focus of multiple conservation strategies involving significant resources such as captive breeding programs. Despite the relevance for conservation, no study to date has focused on evaluating geriatric diseases in dunnarts. This study describes the pathology findings in a group of one wild and thirty-five captive-born, mostly geriatric Julia Creek dunnarts that failed to produce offspring over multiple breeding periods. A total of 20 females and 16 males were submitted for a postmortem examination, with ages ranging from 9 to 42 and 12 to 42 months for females and males, respectively. Of these, 10 had unremarkable findings. The most common condition in females was cystic glandular hyperplasia (n = 8), typical of hormonal dysregulation profiles in senescence, particularly hyperestrogenism. Rarely, cutaneous disease represented by unidentified dermal round cell infiltrates was observed in females (n = 2). Primary reproductive hormonal dysregulation was also suspected in males diagnosed with testicular degeneration, aspermatogenesis and/or atrophy (n = 3). Cutaneous round cell infiltrates, possibly compatible with epitheliotropic lymphomas, were seen in males (n = 3), and 2/3 affected males also had concurrent testicular degeneration or atrophy, indicating male sex could be a predictor for lymphoid neoplasia in aged dunnarts, especially in individuals with concurrent testosterone-luteinizing hormone dysregulation as it occurs in gonadectomized animals. The role of an underlying viral etiology is also explored. This study is the first to describe major spontaneous diseases in endangered aged Julia Creek dunnarts, providing an important understanding of senescence and geriatric diseases within a conservation context.

Predicting Sex in White Rhinoceroses: A Statistical Model for Conservation Management

Ensuring the effective management of every rhinoceros population is crucial for securing a future for the species, especially considering the escalating global threat of poaching and the challenges faced in captive breeding programs for this endangered species. Steroid hormones play pivotal roles in regulating diverse biological processes, making fecal hormonal determinations a valuable non-invasive tool for monitoring adrenal and gonadal endocrinologies and assessing reproductive status, particularly in endangered species. The purpose of this study was to develop a statistical model for predicting the sex of white rhinoceroses using hormonal determinations obtained from a single fecal sample. To achieve this, 562 fecal samples from 15 individuals of the Ceratotherium simum species were collected, and enzyme immunoassays were conducted to determine the concentrations of fecal cortisol, progesterone, estrone, and testosterone metabolites. The biological validation of the method provided an impressive accuracy rate of nearly 80% in predicting the sex of hypothetically unknown white rhinoceroses. Implementing this statistical model for sex identification in white rhinoceroses would yield significant benefits, including a better understanding of the structure and dynamics of wild populations. Additionally, it would enhance conservation management efforts aimed at protecting this endangered species. By utilizing this innovative approach, we can contribute to the preservation and long-term survival of white rhinoceros populations.

Non-invasive Assessment of Fecal Glucocorticoid and Androgen Metabolites in the Pygmy Hippopotamus (Choeropsis liberiensis)

The pygmy hippopotamus (Choeropsis liberiensis) is an endangered species endemic to the Upper Guinea Forest ecosystem in West Africa. We have limited information concerning the species' reproduction and well-being under managed care. We therefore developed non-invasive methods for characterizing gonadal androgen and adrenal hormone profiles in pygmy hippos using fecal samples collected from 12 males and 12 females housed in North American zoological institutions. We aimed to: 1) identify and validate enzyme immunoassays (EIAs) for measuring metabolites of corticosteroids and testosterone in feces; and 2) test whether gonadal activity is correlated with previous breeding history, season or type of housing. For glucocorticoids, several EIAs for measuring metabolites were investigated. A group-specific EIA exhibiting cross-reactivity with 11,17-dioxoandrostane (DOA) metabolites of cortisol most clearly reflected adrenocortical activity in response to ACTH¹ challenge in both males and females. However, day-to-day concentrations of this metabolite in the feces of pygmy hippos that did not undergo ACTH challenge were near the detection limits of the assay, making this EIA impractical for assessing glucocorticoid activity in this species. Another group-specific EIA, exhibiting cross-reactivity with 5α-pregnane-3β,11β,21-triol-20-one, produced biologically relevant data and evidence of an appropriate response to pharmacological challenge with exogenous ACTH. The testosterone metabolite assay C196 (Arbor Assays, Ann Arbor, Michigan, USA) also produced biologically coherent data: adult males exhibited the highest mean androgen metabolite concentrations (477 ng/g), followed by adult females (259 ng/g) and juvenile males (160 ng/g). Proven breeding males had higher, but not significantly different, mean concentrations (472 ng/g) to unproven males (352 ng/g; P = 0.400). Similarly, adult males housed outdoors year-round in subtropical climates exhibited higher, but not statistically different mean concentrations (554 ng/g) to males in temperate climates that were housed indoors at least part of the year (412 ng/g; P = 0.208). There were, however, significant differences in mean concentrations among seasons for adult males, with higher values in spring (546 ng/g) and summer (542 ng/g) than in autumn (426 ng/g) and winter (388 ng/g, P = 0.003). In conclusion, we identified EIAs for the measurement of fecal metabolites of androgens and glucocorticoids that can be used for further studies to monitor gonadal activity in male pygmy hippos and adrenocortical activity in both sexes. We also identified a seasonal trend in male gonadal activity in this species under managed care in North America. Finally, our findings highlight an important consideration when using non-invasive methods for evaluating fecal cortisol metabolites: ACTH used for pharmacological validation of an EIA does not necessarily equate to biological relevance.

Mothers of disabled infants had higher cortisol levels in a free-ranging group of Japanese macaques (Macaca fuscata)

Glucocorticoids (GCs) are hormones released in response to stressors and can provide insight into an organism's physiological well-being. Experiencing chronic challenges to homeostasis is associated with significant deviations from baseline fecal GCs (fGCs) in many species, providing a noninvasive biomarker for assessing stress. In the group of free-ranging Japanese macaques (Macaca fuscata) at the Awajishima Monkey Center in Japan, ~17% have congenital limb malformations. We collected 646 fecal samples from 27 females over three consecutive birth seasons (May-August) and analyzed them using enzyme immunoassay to extract fGCs. We explored the relationship between fGC levels and individual (physical impairment and reproductive status), social (dominance rank and availability of kin for social support), and ecological variables (exposure to potential predators, rainfall, and wild fruit availability). A disabled infant was associated significantly with higher fGC in the mother; however, physical impairment in adult females was not significantly related to fGC levels. Females with higher dominance rank had significantly lower fGC levels than lower ranking females. Other factors did not relate significantly to fGC. These results suggest that providing care that meets the support needs of disabled infants poses a physiological challenge for mothers and suggests that physically impaired adults are able to effectively compensate for their disabilities with behavioral plasticity. Once an individual with congenital limb malformations survives infancy through their mother's care, physical impairment does not appear to influence fGC values, while social variables like dominance rank significantly influenced cortisol values in free-ranging female Japanese macaques.

Skin Cortisol and Acoustic Activity: Potential Tools to Evaluate Stress and Welfare in Captive Cetaceans

As people's focus broadens from animals on farms to zoos and aquaria, the field of welfare science and the public's concern for animal welfare continue to grow. In captive animals, stress and its causes are topics of interest in welfare issues, and the identification of an objective method that can be used to assess animals' stress as a physiological state is essential. Both behavioral and physiological parameters can be used as indicators in order to assess animal stress quantitatively. To validate this approach, acoustic activity and the sloughed scrape skin cortisol concentration were used to evaluate the animal welfare of captive beluga whales (Delphinapterus leucas). The acoustic activity (5 min at 10:00 am) of three captive D. leucas was routinely recorded by a transducer and analyzed using audio editing software. The calls were separated into three main categories: whistles, pulses, and combo calls. The sloughed scrape skin samples were collected non-invasively once a week from all three animals' fluke and/or flipper. Cortisol was extracted using a modified skin steroid extraction technique, and detected via commercially available enzyme immunoassays. The results showed that the cortisol concentration increased by varying levels when the whales encountered the same event. In addition, the number and distribution of the calls changed along with the events. This indicated that the changes in the cortisol concentration and acoustic behavior may have reflected the fluctuations in the environment and body condition. Therefore, the scrape cortisol measurement and acoustic recordings could be used to monitor stress levels in captive beluga whales. We recommend that aquaria consider incorporating skin scrape cortisol and acoustic activity monitoring into their standards for animal welfare.

Temporal stability of fecal cortisol metabolites in mountain-dwelling ungulates

Fecal cortisol metabolites (FCMs) are widely used to track stress responses in wildlife and captive species. Rules of thumb suggest that samples should be collected as soon as possible after defecation, to avoid decay of FCMs. To date, however, only a few studies investigated the stability of defecated FCMs over time, and most of them were conducted in controlled laboratory conditions. Here, we investigated the stability of FCMs over seven consecutive days, in two mountain-dwelling ungulates, under natural environmental conditions using a semi-experimental approach. Fecal samples from Northern chamois Rupicapra rupicapra (n = 24) and red deer Cervus elaphus (n = 22) were collected in summer of 2020 within the Stelvio National Park, Italy, and placed in an open area above 2000 m a.s.l. For the next 7 days, we collected a portion of each sample, and all sub-samples were analyzed with an 11-oxoetiocholanolone enzyme immunoassay. Exposure, temperature, and precipitation were fitted as covariates in non-linear generalized mixed models to assess FCM variation over time, and competing models were selected using AICc. For chamois, the best model included only time as a predictor, while for red deer, it included time, precipitation, and exposure. For both species, FCM values decreased rapidly from the first days after deposition until the fourth day. For red deer, in northern-exposed samples, FCM values decreased slower than in south-exposed ones; furthermore, FCM values increased with increasing precipitation. Our results offer a solid methodological basis to wildlife researchers and practitioners interested in the investigation of the ecological factors affecting stress variation in wildlife and support the recommendation to collect samples as fresh as possible, to avoid misleading inference. Further studies are necessary to evaluate the stability of FCMs when other enzyme immunoassays are used.

The Application of NIRS to Determine Animal Physiological Traits for Wildlife Management and Conservation

The ability to measure and monitor wildlife populations is important for species management and conservation. The use of near-infrared spectroscopy (NIRS) to rapidly detect physiological traits from wildlife scat and other body materials could play an important role in the conservation of species. Previous research has demonstrated the potential for NIRS to detect diseases such as the novel COVID-19 from saliva, parasites from feces, and numerous other traits from animal skin, hair, and scat, such as cortisol metabolites, diet quality, sex, and reproductive status, that may be useful for population monitoring. Models developed from NIRS data use light reflected from a sample to relate the variation in the sample’s spectra to variation in a trait, which can then be used to predict that trait in unknown samples based on their spectra. The modelling process involves calibration, validation, and evaluation. Data sampling, pre-treatments, and the selection of training and testing datasets can impact model performance. We review the use of NIRS for measuring physiological traits in animals that may be useful for wildlife management and conservation and suggest future research to advance the application of NIRS for this purpose.

The benefits of being dominant: health correlates of male social rank and age in a marmot

The benefits of dominance may not come without costs, particularly for males. For example, the “immunocompetence handicap hypothesis” states that males with enhanced mating success allocate resources to enhance reproductive output at a cost to their current health, whereas the “resource quality hypothesis” predicts that high-ranking males may benefit from increased reproduction and good health. Whereas the predictions from each have been well tested in captive animals and in a variety of highly social primates, fewer studies have been carried out in free-living, facultatively social animals. Using adult male yellow-bellied marmots (Marmota flaviventer), we evaluated predictions of these hypotheses by examining the relationship between social rank and 2 health indicators—fecal glucocorticoid metabolite (FCM) levels, and neutrophil/lymphocyte (N/L) ratios—after accounting for variation explained by age, body mass, and seasonality. We found that higher-ranking males tended to have a lower N/L ratio (reflecting good health) than lower-ranking individuals, whereas FCM levels were not significantly related to rank. In addition, heavier male marmots had lower N/L ratios, whereas body mass was not associated with FCM levels. We also found that older adult males had lower FCM levels (reflecting less physiological stress) but higher N/L ratios than younger adults. Finally, we found that FCM levels decreased as the active season progressed and FCM levels were associated with the time of the day. Overall, our results suggest that socially-dominant male marmots enjoyed better, not worse health in terms of lower N/L ratios.

Monitoring corticoid metabolites in urine of white-tailed sea eagles: Negative effects of road proximity on breeding pairs

The white-tailed sea eagle (Haliaeetus albicilla) is known to be sensitive to disturbance. To better understand potential stressors, we measured corticosterone metabolite levels in H. albicilla excreta and recorded the nest success of breeding pairs. We tested the ability of four enzyme immunoassays (EIA) to measure urinary glucocorticoid metabolites (uGM) in the excreta of one adult female eagle subjected to a controlled physiological stress treatment. We identified corticosterone-21-HS to be the most sensitive EIA to changes in uGM concentration. To exclude a sex bias, we confirmed the assay's applicability with samples collected from similar stress treatments in two juvenile males. We used the identified EIA to measure uGM in wild breeding pairs and tested effects of disturbance. Breeding pairs nesting closer to roads and paths had higher uGM concentrations (p = 0.02), which is likely an effect of human recreational activity and disturbance. There was no difference in uGM concentrations between failed and successful nests. Our results highlight the potential impact of road and path proximity on white-tailed sea eagles, with potential importance for species management and conservation, particularly with respect to nest protection zone legislation.

Faecal metabolites and hair cortisol as biological markers of HPA-axis activity in the Rocky mountain goat

Monitoring glucocorticoids in faeces and hair is increasingly used in ecological studies and provides a powerful and minimally intrusive mean to identify physiological challenges faced by wild animals. Using a cortisol and a corticosterone immunoassays, we conducted an adrenocorticotropic (ACTH) challenge with five weekly repeated injections to validate the use of faecal glucocorticoid metabolites and hair cortisol concentration as biological markers of the HPA-axis activity in captive mountain goats (Oreamnos americanus). We also investigated the effect of endogenous (age, sex, reproductive status) and methodological (faecal sample collection date, freezing delay and hair type) variables on cortisol values using faecal and hair samples collected from marked wild mountain goats during a long-term study. The cortisol enzyme immunoassay was reliable for mountain goat faeces and hair, and was sensitive enough to detect a clear rise in glucocorticoid concentration following ACTH injections for both matrices. Age and sex had no detectable effect on faecal glucocorticoid metabolites, but hair cortisol concentration was higher in kids and yearlings than in older goats, and lower in adult males compared to adult females. Reproductive status had no detectable effect on both faecal and hair measurements. Faecal metabolite concentrations increased with sample collection date in late spring until mid-summer and decreased afterward until early fall. Guard hair had nearly twice as much cortisol per gram as undercoat hair. Prolonged delay to freezing reduced the concentration of faecal glucocorticoid metabolites, but degradation seemed limited when samples were exposed to wind and sun or when ambient temperature was low. We conclude that faeces and hair can be used as valid biomarkers of the HPA-axis activity in mountain goat provided that confounding variables are taken into account when interpreting measurements.

Rapid identification of species, sex and maturity by mass spectrometric analysis of animal faeces

BACKGROUND

We describe a new approach to the recovery of information from faecal samples, based on the analysis of the molecular signature generated by rapid evaporative ionisation mass spectrometry (REIMS).

RESULTS

Faecal pellets from five different rodent species were analysed by REIMS, and complex mass spectra were acquired rapidly (typically a few seconds per sample). The uninterpreted mass spectra (signatures) were then used to seed linear discriminant analysis and classification models based on random forests. It was possible to classify each species of origin with a high rate of accuracy, whether faeces were from animals maintained under standard laboratory conditions or wild-caught. REIMS signatures were stable to prior storage of the faecal material under a range of different conditions and were not altered rapidly or radically by changes in diet. Further, within species, REIMS signatures could be used to discriminate faeces from adult versus juvenile mice, male versus female mice and those from three different laboratory strains.

CONCLUSIONS

REIMS offers a completely novel method for the rapid analysis of faecal samples, extending faecal analysis (previously focused on DNA) to an assessment of phenotype, and has considerable potential as a new tool in the armamentarium of the field biologist.

Physiological stress responses of tigers due to anthropogenic disturbance especially tourism in two central Indian tiger reserves

Tigers continue to face unprecedented threats to their existence due to poaching, habitat loss, habitat fragmentation and anthropogenic disturbances. The present study examines the physiological stress response of tigers due to anthropogenic activities including wildlife tourism in Bandhavgarh Tiger Reserve and Kanha Tiger Reserve using faecal glucocorticoid metabolite (fGCM) measurement. We collected a total of 341 faecal samples from both reserves during tourism and non-tourism periods. Data on various anthropogenic disturbances including tourism activities like number of vehicles and visitors were also collected. We ascertained the species identity and sex of all the samples collected using genetic markers. fGCMs were extracted using a previously reported procedure, and fGCM concentrations were subsequently determined using an established enzyme immunoassay. There was no significant difference in overall mean fGCM concentrations between the two tiger reserves, but within each reserve, concentrations were significantly higher in tigers during the tourism period as compared to the non-tourism period. We also found that the number of tourist vehicles and disturbance level significantly correlated with fGCM concentrations. This study further supports the assumption that unbridled tourism associated with high anthropogenic disturbance can be related to perceived stress and consequently may have an impact on the reproductive fitness of tigers and long-term survival of isolated populations.

HPLC-QTOF method for quantifying 11-ketoetiocholanolone, a cortisol metabolite, in ruminants' feces: Optimization and validation

Studies of animal ecology can benefit from a quantified understanding of eco-physiological processes and, in particular, of the physiological responses in free-ranging animals to potential stressors. The determination of fecal cortisol metabolites as a noninvasive method for monitoring stress has proved to be a powerful tool. High-performance liquid chromatography coupled with tandem mass spectrometry (HPLC-MS/MS) has emerged as the most accurate method for avoiding problems related to the nonspecificity of immunoassays. In this study, we optimize and validate a reliable method using HPLC-MS/MS for quantifying 11-ketoetiocholanolone (11-k), a representative fecal cortisol metabolite in ruminants. An appropriate extraction and purification procedure was developed taking into account the complex nature of feces. The final extract obtained was then analyzed with HPLC-MS/MS using a quadrupole-time-of-fly (QTOF) tandem mass spectrometer with an electrospray ionization interface operating in positive mode, which allowed an unequivocal determination of the metabolite due to its accurate mass capabilities. After rigorous optimization of both sample extraction and the HPLC-QTOF parameters, making use of feces from free-ranging Iberian ibex, ideal conditions were established. Matrix-matched standards were used to calibrate the method. The limit of detection and quantification was 13- and 40- ng/g, respectively. The validation of the method was performed with recoveries in the range of 85-110%, a figure much higher than the 60% obtained with the previous extraction methods used in our laboratory, and with relative standard deviations (RSDs) no higher than 15% for the complete analytical procedure, including extraction and analysis. The time required for the fecal 11-k analysis was greatly reduced in comparison with the previous work carried out in our laboratory. This is the first time that QTOF mass detection coupled with HPLC has been validated for 11-k quantification in feces from free-ranging ruminants such as Iberian ibex. Given the high selectivity and sensitivity attained, our method could become a useful tool for noninvasive stress quantification in ruminants.

Seasonal Variations in Heart Rate Variability as an Indicator of Stress in Free-Ranging Pregnant Przewalski's Horses (E. ferus przewalskii) within the Hortobágy National Park in Hungary

Background: Ecosystems with seasonal fluctuations in climate and food availability present physiological challenges to resident mammals and may cause “stress.” The two predominant physiological responses to stressors are (1) the activation of the hypothalamic-pituitary-adrenal axis and (2) the modulation of the autonomic nervous system. To date, the primary indicator for “stress” in wildlife- and zoo animal research are glucocorticoid levels. By measuring the autonomic regulation of cardiac activity, particularly the vagal tone, heart rate variability (HRV) is presently emerging as a suitable indicator of “stress” in farm- and domestic animal research.

Objective: The aim of this study was to use HRV, a novel method in wildlife research, to assess seasonal patterns of “stress” in a group of free-ranging Przewalski's horses (Equus ferus przewalskii).

Methods: Six pregnant Przewalski's horses from one harem within the Hortobágy National Park in Hungary were subjected to the study. We used a dedicated telemetry system consisting of a subcutaneously implanted transmitter and a receiver and storage unit in a collar to record HRV, heart rate (HR), subcutaneous body temperature, and activity throughout a one-year study period—climate data was also collected. We defined “stress” as a decrease in parasympathetic nervous system tone and calculated RMSSD (root mean square of successive differences) as a measure of HRV. Linear mixed effects models with random intercept per individual were used for statistical analysis.

Results: HRV and HR varied considerably throughout the year. Similar to temperate ruminants and hibernating mammals, Przewalski's horses experienced lower HR and HRV during winter, when resources are limited indicating decreased metabolic rates coupled with “stress.” In spring, we observed a drop of HRV along with a peak in HR indicating an increase of allostatic load that is most likely associated with increased energy demands during pregnancy and/or seasonal routines such as the adjustment of the gastrointestinal system to better quality diet.

Conclusion: Measuring telemetric HRV is a proven method to study undisturbed reactions of wild animals to their changing environment over the long term. Przewalski's horses experience a loss of complexity in cardiovascular dynamics over the winter and particularly during spring, indicating seasonal “stress.”

How genetic data improve the interpretation of results of faecal glucocorticoid metabolite measurements in a free-living population

Measurement of glucocorticoid metabolites (GCM) in faeces has become a widely used and effective tool for evaluating the amount of stress experienced by animals. However, the potential sampling bias resulting from an oversampling of individuals when collecting “anonymous” (unknown sex or individual) faeces has rarely been investigated. We used non-invasive genetic sampling (NIGS) to investigate potential interpretation errors of GCM measurements in a free-living population of mountain hares during the mating and post-reproductive periods. Genetic data improved the interpretation of results of faecal GCM measurements. In general GCM concentrations were influenced by season. However, genetic information revealed that it was sex-dependent. Within the mating period, females had higher GCM levels than males, but individual differences were more expressed in males. In the post-reproductive period, GCM concentrations were neither influenced by sex nor individual. We also identified potential pitfalls in the interpretation of anonymous faecal samples by individual differences in GCM concentrations and resampling rates. Our study showed that sex- and individual-dependent GCM levels led to a misinterpretation of GCM values when collecting “anonymous” faeces. To accurately evaluate the amount of stress experienced by free-living animals using faecal GCM measurements, we recommend documenting individuals and their sex of the sampled population. In stress-sensitive and elusive species, such documentation can be achieved by using NIGS and for diurnal animals with sexual and individual variation in appearance or marked individuals, it can be provided by a detailed field protocol.

Selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes in hares (Lepus capensis L., 1758) from a steep ecological gradient in North Africa

Background

Recent studies of selection on mitochondrial (mt) OXPHOS genes suggest adaptation due mainly to environmental variation. In this context, Tunisian hares that display several external phenotypes with phylogenetically rather homogenous gene pool and shallow population structure provide a good precondition to detect positive selection on mt genes related to environmental/climatic variation, specifically ambient temperature and precipitation.

Results

We used codon-based methods along with population genetic data to test for positive selection on ATPase synthase 6 (ATP6) and NADH dehydrogenase 2 (ND2) of cape hares (Lepus capensis) collected along a steep ecological gradient in Tunisia. We found significantly higher differentiation at the ATP6 locus across Tunisia, with sub-humid Mediterranean, semi-arid, and arid Sahara climate than for fourteen unlinked supposedly neutrally evolving nuclear microsatellites and mt control region sequences. This suggested positive selection on ATP6 sequences, which was confirmed by several codon-based tests for one sequence site that together with a second site translated into four different amino acids. Positive selection on ND2 sequences was also confirmed by several codon-based tests. The corresponding frequencies of the two most prevalent variants at each locus varied significantly across climate regions, and our logistic general linear models of occurrence of those proteins indicated significant effects of mean annual temperature for ATP6 and mean minimum temperature of the coldest month of the year for ND2, independent of geographical location, annual precipitation, and the respective co-occurring protein at the second locus. Moreover, presence of the ancestral ATP6 protein, as inferred from phylogenetic networks, was positively affected by the simultaneous presence of the derived ND2 protein and vice versa, independent of temperature, precipitation, or geographic location. Finally, we obtained a significant coevolution signal for the ancestral ATP6 and derived ND2 sequences and vice versa.

Conclusions

positive selection was strongly suggested by the population genetic approach and the codon-based tests in both mtDNA genes. Moreover, the two most prevalent proteins at the ATP6 locus were distributed at significantly varying frequencies across the study area with a significant effect of mean annual temperature on the occurrence of the ATP6 proteins independent of geographical coordinates and the co-occuring ND2 protein variant. For ND2, occurrence of the two most frequent protein variants was significantly influenced by the mean minimum temperature of the coldest month, independent of the co-occurring ATP6 protein variant and geographical coordinates. This strongly suggests direct involvement of ambient temperature in the adaptation of the studied mtOXPHOS genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0896-0) contains supplementary material, which is available to authorized users.