A haplotype-resolved reference genome of a long-distance migratory bat, Pipistrellus nathusii (Keyserling & Blasius, 1839).

We present a complete, chromosome-scale reference genome for the long-distance migratory bat Pipistrellus nathusii. The genome encompasses both haplotypic sets of autosomes and the separation of both sex chromosomes by utilizing highly accurate long-reads and preserving long-range phasing information through the use of three-dimensional chromatin conformation capture sequencing (Hi-C). This genome, accompanied by a comprehensive protein-coding sequence annotation, provides a valuable genomic resource for future investigations into the genomic bases of long-distance migratory flight in bats as well as uncovering the genetic architecture, population structure and evolutionary history of Pipistrellus nathusii. The reference-quality genome presented here gives a fundamental resource to further our understanding of bat genetics and evolution, adding to the growing number of high-quality genetic resources in this field. Here, we demonstrate its use in the phylogenetic reconstruction of the order Chiroptera, and in particular, we present the resources to allow detailed investigations into the genetic drivers and adaptations related to long-distance migration.

Epigenetic signatures of social status in wild female spotted hyenas (Crocuta crocuta).

In mammalian societies, dominance hierarchies translate into inequalities in health, reproductive performance and survival. DNA methylation is thought to mediate the effects of social status on gene expression and phenotypic outcomes, yet a study of social status-specific DNA methylation profiles in different age classes in a wild social mammal is missing. We tested for social status signatures in DNA methylation profiles in wild female spotted hyenas (Crocuta crocuta), cubs and adults, using non-invasively collected gut epithelium samples. In spotted hyena clans, female social status influences access to resources, foraging behavior, health, reproductive performance and survival. We identified 149 differentially methylated regions between 42 high- and low-ranking female spotted hyenas (cubs and adults). Differentially methylated genes were associated with energy conversion, immune function, glutamate receptor signalling and ion transport. Our results provide evidence that socio-environmental inequalities are reflected at the molecular level in cubs and adults in a wild social mammal.

The Addis Ababa Lions: Whole-Genome Sequencing of a Rare and Precious Population.

Lions are widely known as charismatic predators that once roamed across the globe, but their populations have been greatly affected by environmental factors and human activities over the last 150 yr. Of particular interest is the Addis Ababa lion population, which has been maintained in captivity at around 20 individuals for over 75 yr, while many wild African lion populations have become extinct. In order to understand the molecular features of this unique population, we conducted a whole-genome sequencing study on 15 Addis Ababa lions and detected 4.5 million distinct genomic variants compared with the reference African lion genome. Using functional annotation, we identified several genes with mutations that potentially impact various traits such as mane color, body size, reproduction, gastrointestinal functions, cardiovascular processes, and sensory perception. These findings offer valuable insights into the genetics of this threatened lion population.

First mitogenome phylogeny of the sun bear Helarctos malayanus reveals a deep split between Indochinese and Sundaic lineages.

The sun bear Helarctos malayanus is one of the most endangered ursids, and to date classification of sun bear populations has been based almost exclusively on geographic distribution and morphology. The very few molecular studies focussing on this species were limited in geographic scope. Using archival and non-invasively collected sample material, we have added a substantial number of complete or near-complete mitochondrial genome sequences from sun bears of several range countries of the species' distribution. We here report 32 new mitogenome sequences representing sun bears from Cambodia, Thailand, Peninsular Malaysia, Sumatra, and Borneo. Reconstruction of phylogenetic relationships revealed two matrilines that diverged ~295 thousand years ago: one restricted to portions of mainland Indochina (China, Cambodia, Thailand; "Mainland clade"), and one comprising bears from Borneo, Sumatra, Peninsular Malaysia but also Thailand ("Sunda clade"). Generally recent coalescence times in the mitochondrial phylogeny suggest that recent or historical demographic processes have resulted in a loss of mtDNA variation. Additionally, analysis of our data in conjunction with shorter mtDNA sequences revealed that the Bornean sun bear, classified as a distinct subspecies (H. m. euryspilus), does not harbor a distinctive matriline. Further molecular studies of H. malayanus are needed, which should ideally include data from nuclear loci.

First Steps towards the Development of Epigenetic Biomarkers in Female Cheetahs (Acinonyx jubatus).

Free-ranging cheetahs (Acinonyx jubatus) are generally healthy, whereas cheetahs under human care, such as those in zoological gardens, suffer from ill-defined infectious and degenerative pathologies. These differences are only partially explained by husbandry management programs because both groups share low genetic diversity. However, mounting evidence suggests that physiological differences between populations in different environments can be tracked down to differences in epigenetic signatures. Here, we identified differentially methylated regions (DMRs) between free-ranging cheetahs and conspecifics in zoological gardens and prospect putative links to pathways relevant to immunity, energy balance and homeostasis. Comparing epigenomic DNA methylation profiles obtained from peripheral blood mononuclear cells (PBMCs) from eight free-ranging female cheetahs from Namibia and seven female cheetahs living in zoological gardens within Europe, we identified DMRs of which 22 were hypermethylated and 23 hypomethylated. Hypermethylated regions in cheetahs under human care were located in the promoter region of a gene involved in host-pathogen interactions (KLC1) and in an intron of a transcription factor relevant for the development of pancreatic ß-cells, liver, and kidney (GLIS3). The most canonical mechanism of DNA methylation in promoter regions is assumed to repress gene transcription. Taken together, this could indicate that hypermethylation at the promoter region of KLC1 is involved in the reduced immunity in cheetahs under human care. This approach can be generalized to characterize DNA methylation profiles in larger cheetah populations under human care with a more granular longitudinal data collection, which, in the future, could be used to monitor the early onset of pathologies, and ultimately translate into the development of biomarkers with prophylactic and/or therapeutic potential.

Genetic population structure defines wild boar as an urban exploiter species in Barcelona, Spain.

Urban wildlife ecology is gaining relevance as metropolitan areas grow throughout the world, reducing natural habitats and creating new ecological niches. However, knowledge is still scarce about the colonisation processes of such urban niches, the establishment of new communities, populations and/or species, and the related changes in behaviour and life histories of urban wildlife. Wild boar (Sus scrofa) has successfully colonised urban niches throughout Europe. The aim of this study is to unveil the processes driving the establishment and maintenance of an urban wild boar population by analysing its genetic structure. A set of 19 microsatellite loci was used to test whether urban wild boars in Barcelona, Spain, are an isolated population or if gene flow prevents genetic differentiation between rural and urban wild boars. This knowledge will contribute to the understanding of the effects of synurbisation and the associated management measures on the genetic change of large mammals in urban ecosystems. Despite the unidirectional gene flow from rural to urban areas, the urban wild boars in Barcelona form an island population genotypically differentiated from the surrounding rural ones. The comparison with previous genetic studies of urban wild boar populations suggests that forest patches act as suitable islands for wild boar genetic differentiation. Previous results and the genetic structure of the urban wild boar population in Barcelona classify wild boar as an urban exploiter species. These wild boar peri-urban island populations are responsible for conflict with humans and thus should be managed by reducing the attractiveness of urban areas. The management of peri-urban wild boar populations should aim at reducing migration into urban areas and preventing phenotypic changes (either genetic or plastic) causing habituation of wild boars to humans and urban environments.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae.

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes and searching for species-specific structural variants. Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observed species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This argues strongly for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

Assessing Coagulation Parameters in Healthy Asian Elephants (Elephas maximus) from European and Thai Populations.

The Asian elephant population is continuously declining due to several extrinsic reasons in their range countries, but also due to diseases in captive populations worldwide. One of these diseases, the elephant endotheliotropic herpesvirus (EEHV) hemorrhagic disease, is very impactful because it particularly affects Asian elephant calves. It is commonly fatal and presents as an acute and generalized hemorrhagic syndrome. Therefore, having reference values of coagulation parameters, and obtaining such values for diseased animals in a very short time, is of great importance. We analyzed prothrombin time (PT), activated partial thromboplastin time (aPTT), and fibrinogen concentrations using a portable and fast point-of-care analyzer (VetScan Pro) in 127 Asian elephants from Thai camps and European captive herds. We found significantly different PT and aPTT coagulation times between elephants from the two regions, as well as clear differences in fibrinogen concentration. Nevertheless, these alterations were not expected to have biological or clinical implications. We have also sequenced the coagulation factor VII gene of 141 animals to assess the presence of a previously reported hereditary coagulation disorder in Asian elephants and to investigate the presence of other mutations. We did not find the previously reported mutation in our study population. Instead, we discovered the presence of several new single nucleotide polymorphisms, two of them being considered as deleterious by effect prediction software.

Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx.

Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species' evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies' current range. L. l. dinniki harbors the second most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species.

Determinants of defence strategies of a hibernating European bat species towards the fungal pathogen Pseudogymnoascus destructans.

Pseudogymnoascus destructans (Pd), the causative agent of white-nose syndrome in North America, has decimated bat populations within a decade. The fungus impacts bats during hibernation when physiological functions, including immune responses, are down-regulated. Studies have shown that Pd is native to Europe, where it is not associated with mass mortalities. Moreover, genomic and proteomic studies indicated that European bats may have evolved an effective immune defence, which is lacking in North American bats. However, it is still unclear which defence strategy enables European bats to cope with the pathogen. Here, we analyzed selected physiological and immunological parameters in torpid, Pd infected European greater mouse-eared bats (Myotis myotis) showing three different levels of infection (asymptomatic, mild and severe symptoms). From a subset of the studied bats we tracked skin temperatures during one month of hibernation. Contrasting North American bats, arousal patterns remained unaffected by Pd infections in M. myotis. In general, heavier M. myotis aroused more often from hibernation and showed less severe disease symptoms than lean individuals; most likely because heavy bats were capable of reducing the Pd load more effectively than lean individuals. In the blood of severely infected bats, we found higher gene expression levels of an inflammatory cytokine (IL-1ß), but lower levels of an acute phase protein (haptoglobin), reactive oxygen metabolites (ROMs) and plasma non-enzymatic antioxidant capacity (OXY) compared to conspecifics with lower levels of infection. We conclude that M. myotis, and possibly also other European bat species, tolerate Pd infections during torpor by using selected acute phase response parameters at baseline levels, yet without arousing from torpor and without synthesizing additional immune molecules.

Intravital Dynamic and Correlative Imaging of Mouse Livers Reveals Diffusion-Dominated Canalicular and Flow-Augmented Ductular Bile Flux.

BACKGROUND AND AIMS: Small-molecule flux in tissue microdomains is essential for organ function, but knowledge of this process is scant due to the lack of suitable methods. We developed two independent techniques that allow the quantification of advection (flow) and diffusion in individual bile canaliculi and in interlobular bile ducts of intact livers in living mice, namely fluorescence loss after photoactivation and intravital arbitrary region image correlation spectroscopy. APPROACH AND RESULTS: The results challenge the prevailing "mechano-osmotic" theory of canalicular bile flow. After active transport across hepatocyte membranes, bile acids are transported in the canaliculi primarily by diffusion. Only in the interlobular ducts is diffusion augmented by regulatable advection. Photoactivation of fluorescein bis-(5-carboxymethoxy-2-nitrobenzyl)-ether in entire lobules demonstrated the establishment of diffusive gradients in the bile canalicular network and the sink function of interlobular ducts. In contrast to the bile canalicular network, vectorial transport was detected and quantified in the mesh of interlobular bile ducts. CONCLUSIONS: The liver consists of a diffusion-dominated canalicular domain, where hepatocytes secrete small molecules and generate a concentration gradient and a flow-augmented ductular domain, where regulated water influx creates unidirectional advection that augments the diffusive flux.

The complete mitochondrial genome of the meerkat (Suricata suricatta) and its phylogenetic relationship with other feliform species.

The meerkat, Suricata suricatta, is a highly social member of the mongoose family (Herpestidae) and the only extant species of the genus Suricata. We present the first complete mitochondrial genome of the meerkat, assembled with a seed-and-extend algorithm using three closely related species as references. Phylogenetic analyses using 22 mitochondrial genome sequences confirm the position of meerkat within the Herpestidae family and the Feliformia, a suborder of Carnivora, with high support values. This position is in good agreement with formerly conducted studies based on a small number of mitochondrial and nuclear gene fragments. Our complete mitochondrial genome represents a valuable resource for further phylogenetic studies, especially of the underrepresented members of the Herpestidae family.

Unexpected Gene-Flow in Urban Environments: The Example of the European Hedgehog.

We use the European hedgehog (Erinaceus europaeus), a mammal with limited mobility, as a model species to study whether the structural matrix of the urban environment has an influence on population genetic structure of such species in the city of Berlin (Germany). Using ten established microsatellite loci we genotyped 143 hedgehogs from numerous sites throughout Berlin. Inclusion of all individuals in the cluster analysis yielded three genetic clusters, likely reflecting spatial associations of kin (larger family groups, known as gamodemes). To examine the potential bias in the cluster analysis caused by closely related individuals, we determined all pairwise relationships and excluded close relatives before repeating the cluster analysis. For this data subset (N = 65) both clustering algorithms applied (Structure, Baps) indicated the presence of a single genetic cluster. These results suggest that the high proportion of green patches in the city of Berlin provides numerous steppingstone habitats potentially linking local subpopulations. Alternatively, translocation of individuals across the city by hedgehog rescue facilities may also explain the existence of only a single cluster. We therefore propose that information about management activities such as releases by animal rescue centres should include location data (as exactly as possible) regarding both the collection and the release site, which can then be used in population genetic studies.

Genomic Adaptations and Evolutionary History of the Extinct Scimitar-Toothed Cat, Homotherium latidens.

Homotherium was a genus of large-bodied scimitar-toothed cats, morphologically distinct from any extant felid species, that went extinct at the end of the Pleistocene [1-4]. They possessed large, saber-form serrated canine teeth, powerful forelimbs, a sloping back, and an enlarged optic bulb, all of which were key characteristics for predation on Pleistocene megafauna [5]. Previous mitochondrial DNA phylogenies suggested that it was a highly divergent sister lineage to all extant cat species [6-8]. However, mitochondrial phylogenies can be misled by hybridization [9], incomplete lineage sorting (ILS), or sex-biased dispersal patterns [10], which might be especially relevant for Homotherium since widespread mito-nuclear discrepancies have been uncovered in modern cats [10]. To examine the evolutionary history of Homotherium, we generated a ∼7x nuclear genome and a ∼38x exome from H. latidens using shotgun and target-capture sequencing approaches. Phylogenetic analyses reveal Homotherium as highly divergent (∼22.5 Ma) from living cat species, with no detectable signs of gene flow. Comparative genomic analyses found signatures of positive selection in several genes, including those involved in vision, cognitive function, and energy consumption, putatively consistent with diurnal activity, well-developed social behavior, and cursorial hunting [5]. Finally, we uncover relatively high levels of genetic diversity, suggesting that Homotherium may have been more abundant than the limited fossil record suggests [3, 4, 11-14]. Our findings complement and extend previous inferences from both the fossil record and initial molecular studies, enhancing our understanding of the evolution and ecology of this remarkable lineage.

Epigenomics and gene regulation in mammalian social systems.

Social epigenomics is a new field of research that studies how the social environment shapes the epigenome and how in turn the epigenome modulates behavior. We focus on describing known gene-environment interactions (GEIs) and epigenetic mechanisms in different mammalian social systems. To illustrate how epigenetic mechanisms integrate GEIs, we highlight examples where epigenetic mechanisms are associated with social behaviors and with their maintenance through neuroendocrine, locomotor, and metabolic responses. We discuss future research trajectories and open questions for the emerging field of social epigenomics in nonmodel and naturally occurring social systems. Finally, we outline the technological advances that aid the study of epigenetic mechanisms in the establishment of GEIs and vice versa.

The boon and bane of boldness: movement syndrome as saviour and sink for population genetic diversity.

BACKGROUND: Many felid species are of high conservation concern, and with increasing human disturbance the situation is worsening. Small isolated populations are at risk of genetic impoverishment decreasing within-species biodiversity. Movement is known to be a key behavioural trait that shapes both demographic and genetic dynamics and affects population survival. However, we have limited knowledge on how different manifestations of movement behaviour translate to population processes. In this study, we aimed to 1) understand the potential effects of movement behaviour on the genetic diversity of small felid populations in heterogeneous landscapes, while 2) presenting a simulation tool that can help inform conservation practitioners following, or considering, population management actions targeting the risk of genetic impoverishment. METHODS: We developed a spatially explicit individual-based population model including neutral genetic markers for felids and applied this to the example of Eurasian lynx. Using a neutral landscape approach, we simulated reintroductions into a three-patch system, comprising two breeding patches separated by a larger patch of differing landscape heterogeneity, and tested for the effects of various behavioural movement syndromes and founder population sizes. We explored a range of movement syndromes by simulating populations with various movement model parametrisations that range from 'shy' to 'bold' movement behaviour. RESULTS: We find that movement syndromes can lead to a higher loss of genetic diversity and an increase in between population genetic structure for both "bold" and "shy" movement behaviours, depending on landscape conditions, with larger decreases in genetic diversity and larger increases in genetic differentiation associated with bold movement syndromes, where the first colonisers quickly reproduce and subsequently dominate the gene pool. In addition, we underline the fact that a larger founder population can offset the genetic losses associated with subpopulation isolation and gene pool dominance. CONCLUSIONS: We identified a movement syndrome trade-off for population genetic variation, whereby bold-explorers could be saviours - by connecting populations and promoting panmixia, or sinks - by increasing genetic losses via a 'founder takes all' effect, whereas shy-stayers maintain a more gradual genetic drift due to their more cautious behaviour. Simulations should incorporate movement behaviour to provide better projections of long-term population viability and within-species biodiversity, which includes genetic diversity. Simulations incorporating demographics and genetics have great potential for informing conservation management actions, such as population reintroductions or reinforcements. Here, we present such a simulation tool for solitary felids.

Tissue-specific epigenetic inheritance after paternal heat exposure in male wild guinea pigs.

External temperature change has been shown to modify epigenetic patterns, such as DNA methylation, which regulates gene expression. DNA methylation is heritable, and as such provides a mechanism to convey environmental information to subsequent generations. Studies on epigenetic response to temperature increase are still scarce in wild mammals, even more so studies that compare tissue-specific epigenetic responses. Here, we aim to address differential epigenetic responses on a gene and gene pathway level in two organs, liver and testis. We chose these organs, because the liver is the main metabolic and thermoregulation organ, and epigenetic modifications in testis are potentially transmitted to the F2 generation. We focused on the transmission of DNA methylation changes to naive male offspring after paternal exposure to an ambient temperature increase of 10 °C, and investigated differential methylated regions of sons sired before and after the paternal exposure using Reduced Representation Bisulfite Sequencing. We detected both a highly tissue-specific epigenetic response, reflected in genes involved in organ-specific metabolic pathways, and a more general regulation of single genes epigenetically modified in both organs. We conclude that genomes are context-specifically differentially epigenetically regulated in response to temperature increase. These findings emphasize the epigenetic relevance in cell differentiation, which is essential for the specific function(s) of complex organs, and is represented in a diverse molecular regulation of genes and gene pathways. The results also emphasize the paternal contribution to adaptive processes.

Immune response of hibernating European bats to a fungal challenge.

Immunological responses of hibernating mammals are suppressed at low body temperatures, a possible explanation for the devastating effect of the white-nose syndrome on hibernating North American bats. However, European bats seem to cope well with the fungal causative agent of the disease. To better understand the immune response of hibernating bats, especially against fungal pathogens, we challenged European greater mouse-eared bats (Myotis myotis) by inoculating the fungal antigen zymosan. We monitored torpor patterns, immune gene expressions, different aspects of the acute phase response and plasma oxidative status markers, and compared them with sham-injected control animals at 30 min, 48 h and 96 h after inoculation. Torpor patterns, body temperatures, body masses, white blood cell counts, expression of immune genes, reactive oxygen metabolites and non-enzymatic antioxidant capacity did not differ between groups during the experiment. However, zymosan injected bats had significantly higher levels of haptoglobin than the control animals. Our results indicate that hibernating greater mouse-eared bats mount an inflammatory response to a fungal challenge, with only mild to negligible consequences for the energy budget of hibernation. Our study gives a first hint that hibernating European bats may have evolved a hibernation-adjusted immune response in order to balance the trade-off between competent pathogen elimination and a prudent energy-saving regime.

Non-invasive faecal sampling reveals spatial organization and improves measures of genetic diversity for the conservation assessment of territorial species: Caucasian lynx as a case species.

The Caucasian lynx, Lynx lynx dinniki, has one of the southernmost distributions in the Eurasian lynx range, covering Anatolian Turkey, the Caucasus and Iran. Little is known about the biology and the genetic status of this subspecies. To collect baseline genetic, ecological and behavioural data and benefit future conservation of L. l. dinniki, we monitored 11 lynx territories (396 km2) in northwestern Anatolia. We assessed genetic diversity of this population by non-invasively collecting 171 faecal samples and trapped and sampled 12 lynx individuals using box traps. We observed high allelic variation at 11 nuclear microsatellite markers, and found no signs of inbreeding despite the potential isolation of this population. We obtained similar numbers of distinct genotypes from the two sampling sources. Our results indicated that first order female relatives occupy neighbouring territories (female philopatry) and that territorial male lynx were highly unrelated to each other and to female territorial lynx, suggesting long distance male dispersal. Particular male and female resident territorial lynx and their offspring (kittens and subadults) were more likely to be trapped than resident floaters or dispersing (unrelated) lynx. Conversely, we obtained more data for unrelated lynx and higher numbers of territorials using non-invasive sampling (faeces). When invasive and non-invasive samples were analysed separately, the spatial organisation of lynx (in terms of female philopatry and females and males occupying permanent ranges) affected measures of genetic diversity in such a way that estimates of genetic diversity were reduced if only invasive samples were considered. It appears that, at small spatial scales, invasive sampling using box traps may underestimate the genetic diversity in carnivores with permanent ranges and philopatry such as the Eurasian lynx. As non-invasive sampling can also provide additional data on diet and spatial organisation, we advocate the use of such samples for conservation genetic studies of vulnerable, endangered or data deficient territorial species.

Environmental Change-Dependent Inherited Epigenetic Response.

Epigenetic modifications are a mechanism conveying environmental information to subsequent generations via parental germ lines. Research on epigenetic responses to environmental changes in wild mammals has been widely neglected, as well as studies that compare responses to changes in different environmental factors. Here, we focused on the transmission of DNA methylation changes to naive male offspring after paternal exposure to either diet (~40% less protein) or temperature increase (10 °C increased temperature). Because both experiments focused on the liver as the main metabolic and thermoregulation organ, we were able to decipher if epigenetic changes differed in response to different environmental changes. Reduced representation bisulfite sequencing (RRBS) revealed differentially methylated regions (DMRs) in annotated genomic regions in sons sired before (control) and after the fathers' treatments. We detected both a highly specific epigenetic response dependent on the environmental factor that had changed that was reflected in genes involved in specific metabolic pathways, and a more general response to changes in outer stimuli reflected by epigenetic modifications in a small subset of genes shared between both responses. Our results indicated that fathers prepared their offspring for specific environmental changes by paternally inherited epigenetic modifications, suggesting a strong paternal contribution to adaptive processes.