Lifelong leukocyte telomere dynamics and survival in a free‐living mammal

Telomeres play a fundamental role in the maintenance of genomic integrity at a cellular level, and average leukocyte telomere length (LTL) has been proposed as a biomarker of organismal aging. However, studies tracking LTL across the entire life course of individuals are lacking. Here, we examined lifelong patterns of variation in LTL among four birth cohorts of female Soay sheep (Ovis aries) that were longitudinally monitored and sampled from birth to death. Over the first 4 months of life, there was within‐individual loss of LTL, consistent with findings in the human and primate literature, but there was little evidence of consistent LTL loss associated with age after this point. Overall, we observed only weak evidence of individual consistency in LTL across years and over the entire lifespan: Within‐individual variation was considerable, and birth cohorts differed markedly in their telomere dynamics. Despite the high levels of LTL variation within the lifetimes of individuals, there remained significant associations between LTL and longevity. Detailed analysis of the longitudinal data set showed that this association was driven by improved survival of individuals with longer LTL over the first 2 years of life. There was no evidence that LTL predicted survival in later adulthood. Our data provide the first evidence from a mammal that LTL can predict mortality and lifespan under natural conditions, and also highlight the potentially dynamic nature of LTL within the lifetimes of individuals experiencing a complex and highly variable environment.     

Fecal antibody levels as a noninvasive method for measuring immunity to gastrointestinal nematodes in ecological studies

Among‐individual variation in antibody‐associated immunity to gastrointestinal nematode parasites (GIN) is known be associated with life‐history traits and vital rates in wild vertebrate systems. To date, measurement of levels of antibodies against GIN antigens in natural populations has exclusively been based on invasive blood sampling techniques. Previous work in laboratory rodents and ruminant livestock suggests that antibody measures from feces may provide a viable noninvasive approach. We measured total and anti‐GIN antibodies of different isotypes (immunoglobulin (Ig) G, IgA and IgE) from paired samples of plasma and feces from free‐living Soay sheep of different ages and sexes. We tested the correlations among these measures as well as their associations with body mass and Strongyle nematode fecal egg counts (FEC). Significant positive correlations were present among plasma and fecal anti‐GIN antibody levels for IgG and IgA. Generally, correlations between total antibody levels in plasma and feces were weaker and not significant. No significant relationships were found between any antibody measures and body mass; however, fecal anti‐GIN antibody levels were significantly negatively correlated with FEC. Our data clearly demonstrate the feasibility of measuring anti‐GIN antibodies from fecal samples collected in natural populations. Although associations of fecal antibody levels with their plasma counterparts and FEC were relatively weak, the presence of significant correlations in the predicted direction in a relatively small and heterogeneous sample suggests fecal antibody measures could be a useful, noninvasive addition to current eco‐immunological studies.     

Sex‐ And tissue‐specific differences in telomere length in a reptile

The usage of telomere length (TL) in blood as a proxy for the TL of other tissues relies on the assumption that telomere dynamics across all tissues are similar. However, telomere attrition can be caused by reactive oxygen species (ROS) which may vary with metabolic rate, which itself varies across organs depending upon the life history strategy of an organism. Thus, we chose to measure the telomeres of various cell types in juvenile painted dragon lizards, Ctenophorus pictus, given their unusual life history strategy. Individuals typically only experience a single mating season. We measured the TL of male and female dragons using qPCR and observed that TL varied with tissue type and sex. Telomeres of blood cells were longer than those of liver, heart, brain, and spleen, and females had longer telomeres than males. Brain telomeres in males were approximately half the length of those in females. Telomeric attrition in the male brain may be due to the need for rapid learning of reproductive tactics (territory patrol and defense, mate‐finding). Significant correlations between the TL of tissue types suggest that blood TL may be a useful proxy for the TL of other tissues. Our comparison of organ‐specific telomere dynamics, the first in a reptile, suggests that the usage of blood TL as a proxy requires careful consideration of the life history strategy of the organism.     

Mitochondrial DNA and Y‐chromosome diversity in East Adriatic sheep

Variation in mitochondrial DNA (mtDNA) and Y‐chromosome haplotypes was analysed in nine domestic sheep breeds (159 rams) and 21 mouflon ( Ovis musimon) sampled in the East Adriatic. Mitochondrial DNA analyses revealed a high frequency of type B haplotypes, predominantly in European breeds, and a very low frequency of type A haplotypes, which are more frequent in some Asian breeds. Mitochondrial haplotype Hmt‐3 was the most frequent (26.4%), and 37.1%, 20.8% and 7.6% of rams had haplotypes one, two and three mutations remote from Hmt‐3 respectively. In contrast, Y‐chromosome analyses revealed extraordinary paternal allelic richness: HY‐6, 89.3%; HY‐8, 5.0%; HY‐18, 3.1%; HY‐7, 1.3%; and HY‐5, 1.3%. In fact, the number of haplotypes observed is comparable to the number found in Turkish breeds and greater than the number found in European breeds so far. Haplotype HY‐18 (A‐oY1/135‐SRYM18), identified here for the first time, provides a link between the haplotype HY‐12 (A‐oY1/139‐SRYM18) found in a few rams in Turkey and haplotype HY‐9 (A‐oY1/131‐SRYM18) found in one ram in Ethiopia. All mouflons had type B mtDNA haplotypes, including the private haplotype (Hmt‐55), and all were paternally monomorphic for haplotype HY‐6. Our data support a quite homogeneous maternal origin of East Adriatic sheep, which is a characteristic of European breeds. At the same time, the high number of haplotypes found was surprising and intriguing, and it begs for further analysis. Simultaneous analysis of mtDNA and Y‐chromosome information allowed us to detect a large discrepancy between maternal and paternal lineages in some populations. This is most likely the result of breeder efforts to ‘upgrade’ local populations using rams with different paternal origins.     

Joint associations of blood plasma proteins with overwinter survival of a large mammal

In many wild animal populations, hosts are at risk of parasites and malnutrition and resource costs of defence may be difficult to afford. We postulate that proteins, important in homeostasis and immunity, play a complex but central role in condition dependence and resource costs of mammalian immune defence. To test this, we measured plasma concentrations of albumin, total proteins. Self‐reactive antibodies and parasite‐specific IgG in female Soay sheep. Using a principal component analysis, we found a new metric of condition reflecting individual variation in acquisition, assimilation and/or recycling of plasma proteins that predicted overwinter survival. Controlling for this metric, an age‐dependent trade‐off between antibody titres and protein reserves emerged, indicating costs of mounting an antibody response: younger individuals survived best when prioritising immunity while older individuals fared better when maintaining high‐protein nutritional plane. These findings suggest fascinating roles for protein acquisition and allocation in influencing survival in wild animal populations.     

Cellular and humoral immunity in a wild mammal: Variation with age & sex and association with overwinter survival

Immune defenses are expected to be crucial for survival under the considerable parasite pressures experienced by wild animals. However, our understanding of the association between immunity and fitness in nature remains limited due to both the complexity of the vertebrate immune system and the often‐limited availability of immune reagents in nonmodel organisms. Here, we use methods and reagents developed by veterinary researchers for domestic ungulates on blood samples collected from a wild Soay sheep population, to evaluate an unusually broad panel of immune parameters. Our evaluation included different innate and acquired immune cell types as well as nematode parasite‐specific antibodies of different isotypes. We test how these markers correlate with one another, how they vary with age‐group and sex, and, crucially, whether they predict overwinter survival either within or among demographic groups. We found anticipated patterns of variation in markers with age, associated with immune development, and once these age trends were accounted for, correlations among our 11 immune markers were generally weak. We found that females had higher proportions of naïve T cells and gamma–delta T cells than males, independent of age, while our other markers did not differ between sexes. Only one of our 11 markers predicted overwinter survival: sheep with higher plasma levels of anti‐nematode IgG antibodies were significantly more likely to survive the subsequent high mortality winter, independent of age, sex, or weight. This supports a previous finding from this study system using a different set of samples and shows that circulating antibody levels against ecologically relevant parasites in natural systems represent an important parameter of immune function and may be under strong natural selection. Our data provide rare insights into patterns of variation among age‐ and sex groups in different T‐cell subsets and antibody levels in the wild, and suggest that certain types of immune response—notably those likely to be repeatable within individuals and linked to resistance to ecologically relevant parasites—may be most informative for research into the links between immunity and fitness under natural conditions.     

A genomic investigation of ecological differentiation between free‐living and Drosophila‐associated bacteria

Various bacterial taxa have been identified both in association with animals and in the external environment, but the extent to which related bacteria from the two habitat types are ecologically and evolutionarily distinct is largely unknown. This study investigated the scale and pattern of genetic differentiation between bacteria of the family Acetobacteraceae isolated from the guts of Drosophila fruit flies, plant material and industrial fermentations. Genome‐scale analysis of the phylogenetic relationships and predicted functions was conducted on 44 Acetobacteraceae isolates, including newly sequenced genomes from 18 isolates from wild and laboratory Drosophila. Isolates from the external environment and Drosophila could not be assigned to distinct phylogenetic groups, nor are their genomes enriched for any different sets of genes or category of predicted gene functions. In contrast, analysis of bacteria from laboratory Drosophila showed they were genetically distinct in their universal capacity to degrade uric acid (a major nitrogenous waste product of Drosophila) and absence of flagellar motility, while these traits vary among wild Drosophila isolates. Analysis of the competitive fitness of Acetobacter discordant for these traits revealed a significant fitness deficit for bacteria that cannot degrade uric acid in culture with Drosophila. We propose that, for wild populations, frequent cycling of Acetobacter between Drosophila and the external environment prevents genetic differentiation by maintaining selection for traits adaptive in both the gut and external habitats. However, laboratory isolates bear the signs of adaptation to persistent association with the Drosophila host under tightly defined environmental conditions.     

Heritable variation in circulating glucocorticoids and endocrine flexibility in a free‐living songbird

Phenotypic flexibility is a central way that organisms cope with challenging and changing environments. As endocrine signals mediate many phenotypic traits, heritable variation in hormone levels, or their context‐dependent flexibility, could present an important target for selection. Several studies have estimated the heritability of circulating glucocorticoid levels under acute stress conditions, but little is known about the potential for either baseline hormone levels or rapid endocrine flexibility to evolve. Here, we assessed the potential for selection to operate on the elevation (circulating hormone levels) and flexibility of glucocorticoid reaction norms to acute restraint stress. Multivariate animal models revealed low but significant heritability in baseline (h² = 0.13–0.14) and stress‐induced glucocorticoids (h² = 0.18), and moderate heritability in glucocorticoid flexibility in response to acute stress (h² = 0.38) in free‐living juvenile tree swallows (Tachycineta bicolor; n = 408). Baseline glucocorticoids were not genetically correlated with either stress‐induced glucocorticoids or glucocorticoid flexibility. These findings indicate that baseline glucocorticoids and the acute stress response are distinct traits that can be independently shaped by selection. Microevolutionary changes that influence the expression or flexibility of these endocrine mediators of phenotype may be an important way that populations adapt to changing environments and novel threats.     

Histories of host shifts and cospeciation among free‐living parasitoids of Rhagoletis flies

Host shifts by specialist insects can lead to reproductive isolation between insect populations that use different hosts, promoting diversification. When both a phytophagous insect and its ancestrally associated parasitoid shift to the same novel host plant, they may cospeciate. However, because adult parasitoids are free living, they can also colonize novel host insects and diversify independent of their ancestral host insect. Although shifts of parasitoids to new insect hosts have been documented in ecological time, the long‐term importance of such shifts to parasitoid diversity has not been evaluated. We used a genus of flies with a history of speciation via host shifting (Rhagoletis [Diptera: Tephritidae]) and three associated hymenopteran parasitoid genera (Diachasma, Coptera and Utetes) to examine cophylogenetic relationships between parasitoids and their host insects. We inferred phylogenies of Rhagoletis, Diachasma, Coptera and Utetes and used distance‐based cophylogenetic methods (ParaFit and PACo) to assess congruence between fly and parasitoid trees. We used an event‐based method with a free‐living parasitoid cost model to reconstruct cophylogenetic histories of each parasitoid genus and Rhagoletis. We found that the current species diversity and host–parasitoid associations between the Rhagoletis flies and parasitoids are the primary result of ancient cospeciation events. Parasitoid shifts to ancestrally unrelated hosts primarily occur near the branch tips, suggesting that host shifts contribute to recent parasitoid species diversity but that these lineages may not persist over longer time periods. Our analyses also stress the importance of biologically informed cost models when investigating the coevolutionary histories of hosts and free‐living parasitoids.     

Beta‐globin gene evolution in the ruminants: evidence for an ancient origin of sheep haplotype B

Domestic sheep (Ovis aries) can be divided into two groups with significantly different responses to hypoxic environments, determined by two allelic beta‐globin haplotypes. Haplotype A is very similar to the goat beta‐globin locus, whereas haplotype B has a deletion spanning four globin genes, including beta‐C globin, which encodes a globin with high oxygen affinity. We surveyed the beta‐globin locus using resequencing data from 70 domestic sheep from 42 worldwide breeds and three Ovis canadensis and two Ovis dalli individuals. Haplotype B has an allele frequency of 71.4% in O. aries and was homozygous (BB) in all five wild sheep. This shared ancestry indicates haplotype B is at least 2–3 million years old. Approximately 40 kb of the sequence flanking the ~37‐kb haplotype B deletion had unexpectedly low identity between haplotypes A and B. Phylogenetic analysis showed that the divergent region of sheep haplotype B is remarkably distinct from the beta‐globin loci in goat and cattle but still groups with the Ruminantia. We hypothesize that this divergent ~40‐kb region in haplotype B may be from an unknown ancestral ruminant and was maintained in the lineage to O. aries, but not other Bovidae, evolving independently of haplotype A. Alternatively, the ~40‐kb sequence in haplotype B was more recently acquired by an ancestor of sheep from an unknown non‐Bovidae ruminant, replacing part of haplotype A. Haplotype B has a lower nucleotide diversity than does haplotype A, suggesting a recent bottleneck, whereas the higher frequency of haplotype B suggests a subsequent spread through the global population of O. aries.     

Maternal and paternal genetic diversity of ancient sheep in Estonia from the Late Bronze Age to the post‐medieval period and comparison with other regions in Eurasia

Sheep were among the first domesticated animals to appear in Estonia in the late Neolithic and became one of the most widespread livestock species in the region from the Late Bronze Age onwards. However, the origin and historical expansion of local sheep populations in Estonia remain poorly understood. Here, we analysed fragments of the hypervariable D‐loop of mitochondrial DNA (mtDNA; 213 bp) and the Y‐chromosome SRY gene (130 bp) extracted from 31 archaeological sheep bones dated from approximately 800 BC to 1700 AD. The ancient DNA data of sheep from Estonia were compared with ancient sheep from Finland as well as a set of contemporary sheep breeds from across Eurasia in order to place them in a wider phylogeographical context. The analysis shows that: (i) 24 successfully amplified and analysed mtDNA sequences of ancient sheep cluster into two haplogroups, A and B, of which B is predominant; (ii) four of the ancient mtDNA haplotypes are novel; (iii) higher mtDNA haplotype diversity occurred during the Middle Ages as compared to other periods, a fact concordant with the historical context of expanding international trade during the Middle Ages; (iv) the proportion of rarer haplotypes declined during the expansion of sheep from the Near Eastern domestication centre to the northern European region; (v) three male samples showed the presence of the characteristic northern European haplotype, SNP G‐oY1 of the Y‐chromosome, and represent the earliest occurrence of this haplotype. Our results provide the first insight into the genetic diversity and phylogeographical background of ancient sheep in Estonia and provide basis for further studies on the temporal fluctuations of ancient sheep populations.     

Male reproductive success in a promiscuous mammal: behavioural estimates compared with genetic paternity

Molecular techniques have enabled behavioural ecologists to reassess mating systems from a genetic perspective. Studies of paternity frequently reveal that mating behaviour does not always reflect parentage, and may bring to light alternative mating tactics. Here we present a comparison of behavioural and genetic measures of male reproductive success in a mammalian mating system in which both sexes are highly promiscuous. Rather than having a stable harem social structure, Soay rams (Ovis aries) on the island of Hirta, St Kilda, UK usually consort with individual oestrous ewes sequentially. Not all matings occur between consort pairs, however, and ewes have been seen to mate with up to 10 different rams on the same day. Using locus-specific polymorphism at five protein and 10 microsatellite DNA loci, we determined paternity for 236 lambs born into three cohorts, and compared paternity with estimates of mating success derived from more than one census of rutting behaviour. The correlation between the number of ewes with which each ram was witnessed in consort and the number of paternities assigned was positive and statistically significant, and rams that were observed in consort with a ewe were 18 times more likely to have sired her offspring than other candidate rams. However, most lambs (73%) were not sired by a ram seen in consort with the oestrous mother. Many juveniles, yearlings and some adult rams were rarely seen in consort with ewes, yet were assigned a significant number of paternities. These results suggest that mating tactics differ between age groups, and that alternative mating strategies among adults that do not involve forming consorts with many females also confer mating success. For these reasons, census-based observations of consort associations between individuals cannot be used to accurately estimate individual male reproductive success in this population.     

Occurrence,costs and heritability of delayed selfing in a free‐living flatworm

Evolutionary theory predicts that in the absence of outcrossing opportunities, simultaneously hermaphroditic organisms should eventually switch to self‐fertilization as a form of reproductive assurance. Here, we report the existence of facultative self‐fertilization in the free‐living flatworm Macrostomum hystrix, a species in which outcrossing occurs via hypodermic insemination of sperm into the parenchyma of the mating partner. First, we show that isolated individuals significantly delay the onset of reproduction compared with individuals with outcrossing opportunities (‘delayed selfing’) as predicted by theory. Second, consistent with the idea of M. hystrix being a preferential outcrosser under natural conditions, we report likely costs of selfing manifested via reduced hatchling production and offspring survival. Third, we demonstrate that selfing propensity has a genetic basis in this species, with a heritability estimated at 0.43 ± 0.11. Variation in selfing propensity could arise due to differing costs of inbreeding among families; despite marked inter‐family variation in apparent costs of inbreeding, we found no evidence for such a link. Alternatively, selfing propensity might differ across families because of heritable variation in reproductive traits that determine the likelihood of selfing. We speculate that adaptations to hypodermic insemination under outcrossing, most notably a highly modified copulatory stylet (male copulatory organ) and reduced sperm complexity, could also facilitate facultative selfing in this species.     

Biallelic β‐carotene oxygenase 2 knockout results in yellow fat in sheep via CRISPR/Cas9

The recently emerged CRISPR/Cas9 approach represents an efficient and versatile genome editing tool for producing genetically modified animals. Β‐carotene oxygenase 2 (BCO2) is a key enzyme in the progress of β‐carotene metabolism and is associated with yellow adipose tissue color in sheep. We have recently demonstrated targeted multiplex mutagenesis in sheep and have generated a group of BCO2‐disrupted sheep by zygote injection of the CRISPR/Cas9 components. Here, we show that biallelic modification of BCO2 resulted in yellow fat, compared with the fat color in monoallelic individuals and wild types (snow‐flower white). We subsequently characterized the effects of gene modifications at genetic levels employing sequencing and Western blotting, highlighting the importance of the BCO2 gene for the determination of fat color in sheep. These results indicate that genetic modification via CRISPR/Cas9 holds great potential for validating gene functions as well as for generating desirable phenotypes for economically important traits in livestock.     

Diving behavior in a free‐living,semi‐aquatic herbivore,the Eurasian beaver Castor fiber

Semi‐aquatic mammals have secondarily returned to the aquatic environment, although they spend a major part of their life operating in air. Moving both on land, as well as in, and under water is challenging because such species are considered to be imperfectly adapted to both environments. We deployed accelerometers combined with a depth sensor to study the diving behavior of 12 free‐living Eurasian beavers Castor fiber in southeast Norway between 2009 and 2011 to examine the extent to which beavers conformed with mass‐dependent dive capacities, expecting them to be poorer than wholly aquatic species. Dives were generally shallow (<1 m) and of short duration (<30 s), suggesting that the majority of dives were aerobic. Dive parameters such as maximum diving depth, dive duration, and bottom phase duration were related to the effort during different dive phases and the maximum depth reached. During the descent, mean vectorial dynamic body acceleration (VeDBA—a proxy for movement power) was highest near the surface, probably due to increased upthrust linked to fur‐ and lung‐associated air. Inconsistently though, mean VeDBA underwater was highest during the ascent when this air would be expected to help drive the animals back to the surface. Higher movement costs during ascents may arise from transporting materials up, the air bubbling out of the fur, and/or the animals’ exhaling during the bottom phase of the dive. In a manner similar to other homeotherms, beavers extended both dive and bottom phase durations with diving depth. Deeper dives tended to have a longer bottom phase, although its duration was shortened with increased VeDBA during the bottom phase. Water temperature did not affect diving behavior. Overall, the beavers’ dive profile (depth, duration) was similar to other semi‐aquatic freshwater divers. However, beavers dived for only 2.8% of their active time, presumably because they do not rely on diving for food acquisition.     

Maternal predator‐exposure affects offspring size at birth but not telomere length in a live‐bearing fish

The perception of predation risk could affect prey phenotype both within and between generations (via parental effects). The response to predation risk could involve modifications in physiology, morphology, and behavior and can ultimately affect long‐term fitness. Among the possible modifications mediated by the exposure to predation risk, telomere length could be a proxy for investigating the response to predation risk both within and between generations, as telomeres can be significantly affected by environmental stress. Maternal exposure to the perception of predation risk can affect a variety of offspring traits but the effect on offspring telomere length has never been experimentally tested. Using a live‐bearing fish, the guppy (Poecilia reticulata), we tested if the perceived risk of predation could affect the telomere length of adult females directly and that of their offspring with a balanced experimental setup that allowed us to control for both maternal and paternal contribution. We exposed female guppies to the perception of predation risk during gestation using a combination of both visual and chemical cues and we then measured female telomere length after the exposure period. Maternal effects mediated by the exposure to predation risk were measured on offspring telomere length and body size at birth. Contrary to our predictions, we did not find a significant effect of predation‐exposure neither on female nor on offspring telomere length, but females exposed to predation risk produced smaller offspring at birth. We discuss the possible explanations for our findings and advocate for further research on telomere dynamics in ectotherms.     

Intragroup competition predicts individual foraging specialisation in a group‐living mammal

Individual foraging specialisation has important ecological implications, but its causes in group‐living species are unclear. One of the major consequences of group living is increased intragroup competition for resources. Foraging theory predicts that with increased competition, individuals should add new prey items to their diet, widening their foraging niche (‘optimal foraging hypothesis’). However, classic competition theory suggests the opposite: that increased competition leads to niche partitioning and greater individual foraging specialisation (‘niche partitioning hypothesis’). We tested these opposing predictions in wild, group‐living banded mongooses (Mungos mungo), using stable isotope analysis of banded mongoose whiskers to quantify individual and group foraging niche. Individual foraging niche size declined with increasing group size, despite all groups having a similar overall niche size. Our findings support the prediction that competition promotes niche partitioning within social groups and suggest that individual foraging specialisation may play an important role in the formation of stable social groupings.     

Introgression and the fate of domesticated genes in a wild mammal population

When domesticated species are not reproductively isolated from their wild relatives, the opportunity arises for artificially selected variants to be re‐introduced into the wild. However, the evolutionary consequences of introgression of domesticated genes back into the wild are poorly understood. By combining high‐throughput genotyping with 25 years of long‐term ecological field data, we describe the occurrence and consequences of admixture between a primitive sheep breed, the free‐living Soay sheep of St Kilda, and more modern breeds. Utilizing data from a 50 K ovine SNP chip, together with forward simulations of demographic scenarios, we show that admixture occurred between Soay sheep and a more modern breed, consistent with historical accounts, approximately 150 years ago. Haplotype‐sharing analyses with other breeds revealed that polymorphisms in coat colour and pattern in Soay sheep arose as a result of introgression of genetic variants favoured by artificial selection. Because the haplotypes carrying the causative mutations are known to be under natural selection in free‐living Soay sheep, the admixture event created an opportunity to observe the outcome of a ‘natural laboratory’ experiment where ancestral and domesticated genes competed with each other. The haplotype carrying the domesticated light coat colour allele was favoured by natural selection, while the haplotype associated with the domesticated self coat pattern allele was associated with decreased survival. Therefore, we demonstrate that introgression of domesticated alleles into wild populations can provide a novel source of variation capable of generating rapid evolutionary changes.