- Research Article
- Open Access
Epigenetics and early domestication: differences in hypothalamic DNA methylation between red junglefowl divergently selected for high or low fear of humans
© The Author(s) 2018
Received: 15 August 2017
Accepted: 9 March 2018
Published: 2 April 2018
Domestication of animals leads to large phenotypic alterations within a short evolutionary time-period. Such alterations are caused by genomic variations, yet the prevalence of modified traits is higher than expected if they were caused only by classical genetics and mutations. Epigenetic mechanisms may also be important in driving domesticated phenotypes such as behavior traits. Gene expression can be modulated epigenetically by mechanisms such as DNA methylation, resulting in modifications that are not only variable and susceptible to environmental stimuli, but also sometimes transgenerationally stable. To study such mechanisms in early domestication, we used as model two selected lines of red junglefowl (ancestors of modern chickens) that were bred for either high or low fear of humans over five generations, and investigated differences in hypothalamic DNA methylation between the two populations.
Twenty-two 1-kb windows were differentially methylated between the two selected lines at p < 0.05 after false discovery rate correction. The annotated functions of the genes within these windows indicated epigenetic regulation of metabolic and signaling pathways, which agrees with the changes in gene expression that were previously reported for the same tissue and animals.
Our results show that selection for an important domestication-related behavioral trait such as tameness can cause divergent epigenetic patterns within only five generations, and that these changes could have an important role in chicken domestication.
Animal domestication is a process that occurred in an evolutionary short timespan—less than 10,000 years —and substantially altered several traits such as size, behavior, coat color, and physiological and morphological features in all domesticated species. Similar phenotypic variations tend to reoccur even in widely unrelated species, which are referred to as the ‘domestication phenotype’ [1–3]. This domestication process has altered the genetic structure that underlies these traits and has created distinct patterns that differentiate domestic animals from their wild counterparts [4–7]. Similar phenotypic changes are also seen in recent experimental domestication processes such as in the silver fox . Selection for tameness, or decreased fear of humans, was undoubtedly strong during the early period of domestication and has been proposed to be the major factor that drove the domesticated phenotype, which may, accordingly, have developed partly as correlated selection responses [8–11].
Accumulating evidence shows that phenotypic variation in behavior and other traits within a species or breed can be influenced by epigenetic factors , in addition to genetic factors. For example, variation in personality of the great tit is statistically associated with DNA methylation levels at a dopamine receptor gene . Similarly, in rat pups, the DNA methylation status of the hippocampal glucocorticoid receptor gene is affected by maternal behavior . In domestic chickens, differences in DNA methylation are related to susceptibility to disease , immune responses , growth, and metabolism . Moreover, epigenetic differences (e.g. in red blood cells) can emerge simply after individuals are subjected to different rearing conditions .
DNA methylation, i.e. the addition of a methyl group to cytosine in CpG dinucleotides, has the potential to affect gene expression. Although DNA methylation patterns are generally maintained after cell division, they can sometimes be modified by the action of external stimuli [14, 19]. Environmentally-altered DNA methylation patterns can be transmitted through the germ line  and be stable in somatic tissues over generations [21–25]. Somatic epigenetic differences, whether shaped by the environment or intentionally or unintentionally selected, can in turn affect phenotypic traits. Thus, epigenetic mechanisms could be an important factor in the rapid phenotypic changes that occur during domestication. This is supported, for example, by evidence of significant hypermethylation in purebred dogs compared to wolves , in domestic compared to wild worms , and in domesticated White Leghorn chickens compared to the ancestral red junglefowl (RJF) . Interestingly, in the example of Darwin’s finches, DNA methylation in blood was shown to be a better indicator of evolutionary phenotypic diversification than overall genetic changes (copy number variations) .
The chicken, the world’s most produced species for food, was domesticated from the RJF (Gallus gallus) approximately 8000 years ago and is today one of the most phenotypically diverse species of vertebrates [30–32]. With both wild and domestic chickens being available, they have been used in many comparative studies that evaluated domestication effects [33–35]. These studies have, for example, shown that domesticated chickens have variable coloration patterns, grow faster, reproduce more efficiently, and differ in their social behavior and general activity compared to RJF. In the current experiment, we attempted to experimentally recreate the early domestication process of chickens by selecting RJF bi-directionally for low or high levels of fear of humans . Five generations of repeated selection generated significant phenotypic variation that affected size, behavior, and metabolism [11, 36, 37], as well as transcriptomic differences in both hypothalamic and frontal cortex tissues [38, 39]. Hence, domesticated phenotypes developed in a few generations (five) as side effects of continuous selection for high or low fear of humans.
In this study, DNA samples from the hypothalamus of the fifth generation of a lineage of RJF selected for high or low fear of humans were used to perform methylated DNA immunoprecipitation (MeDIP) sequencing, in order to identify differentially-methylated (DM) regions. Previous studies on the same populations found significant differences in brain gene expression as a consequence of selection . Hence, our hypothesis was that epigenetic differences, measured as differential DNA methylation, would be observed between individuals selected for high or low fear of humans over five generations. We analyzed whether these differences mainly targeted genes and genomic regions that are relevant for tameness and stress.
Animals and sampling
Two unrelated populations of RJF were used to generate a parental population (P0), from which two selection lines were bred for high or low fear of humans for five generations (S5). For a detailed description of the breeding scheme, selection, and housing conditions of animals, see [36, 37]. Briefly, after an outbreeding scheme that lasted for two generations, the P0 generation consisted of approximately 70 birds. The P0 generation was divided into two groups composed of animals with high or low levels of fearfulness towards humans based on scores recorded from a standardized test of fear of humans, described in detail in . Individuals were then bi-directionally selected based on scores of fear of humans, resulting in two lines, i.e. one with a high fear of humans (HFH) and one with a low fear of humans (LFH). Approximately 50 animals per selection line were hatched in each generation from 5 to 10 families, and animals from the two lines were housed and reared together upon hatching in order to standardize rearing conditions. All animals received food and water ad libitum.
Animals from generation S5 were sacrificed at the age of 48 weeks by rapid decapitation. Brains were dissected from the skull and snap-frozen in liquid nitrogen within 10 min, as reported in . To access the hypothalamus, which is located in the inferior region of the brain, the whole brain was dissected into smaller parts. The region between the two optic lobes was rotated upside down to reveal a white butterfly-like pattern, with the thalamus/hypothalamus centered as a small, red protruding structure. In total, 12 S5 individuals were used for this study: six HFH and six LFH birds. Each group consisted of three males and three females, which were chosen randomly from each experimental group. These same 12 individuals had previously been used in studies that examined brain transcriptomic changes [38, 39].
DNA was extracted from hypothalamus tissue using an AllPrep RNA/DNA kit (Qiagen), following the manufacturer’s instructions. In short, approximately 20 to 30 mg of tissue were homogenized with 600 µL of Buffer RTL Plus using FastPrep®-24 (MP Biomedicals). DNA was separated from RNA using AllPrep DNA spin columns, which were kept on ice while RNA was purified for gene expression analysis . 350 µL of Buffer AW1 were added to each spin column, which was then centrifuged for 15 s at 17,000 g. The supernatant was discarded and a mixture of 20 µL of proteinase K (20 mg/mL) and 60 µL of Buffer AW1 was added to each sample. The samples were then incubated for 5 min at room temperature, cleaned with 350 µL of Buffer AW1, centrifuged for 15 s at full speed (17,000 g), and centrifuged again with 500 µL of Buffer AW2 for 2 min at full speed. Columns were dried by an additional centrifugation at 17,000 g for 1 min and then placed in 1.5 mL microcentrifuge tubes. For elution of DNA, 50 µL of Buffer EB were added to the spin column, followed by a 10-min incubation at room temperature and centrifugation for 1 min at 8000 g. DNA concentration and purity were measured using a NanoDrop® ND-2000c (ThermoFisher Scientific).
DNA methylation analysis was performed through enrichment of the methylated fraction of the genome by immunoprecipitation with an anti-methyl-cytosine antibody (MeDIP) , followed by next-generation sequencing (MeDIP-seq) . MeDIP-seq is a validated method for determining DNA methylation that is widely used in epigenetic research, and shows excellent performance for genome coverage compared to equivalent methods . The MeDIP procedure was carried out according to a protocol that was previously optimized in our laboratory using chicken DNA . From each sample, 4 µg DNA were used for the MeDIP capture. DNA was diluted in H2O to a total volume of 8.0 µL, and sonicated at “high” setting for six 30-s intervals using a Bioruptor® Standard sonicator (Diagenode). Samples were then run on a 2% agarose gel for confirmation of fragment lengths. Samples were cleaned from excessively long fragments using a PCR purification kit (Qiagen), following the manufacturer’s instructions, and sample concentrations were measured using a NanoDrop® ND-2000c (ThermoFisher Scientific). The sonicated DNA was diluted with 450 µL of TE buffer (10 mM Tris HCl, pH 7.5; 1 mM EDTA) and heat-denatured at 95 °C for 10 min, followed by cooling on ice for 5 min. To each sample, 51 µL of 10× IP buffer (100 mM NaPhosphate, pH 7; 5 M NaCl, 250 µL Triton-X 100) were added, followed by the addition of 10 µg of antibody (monoclonal mouse anti 5-methylcytosine (2 µg/µL), 5-mC, Diagenode). Then, samples were incubated at 4 °C for 2 h on a rotating platform. Agarose beads (Protein A/G Plus-Agarose, Santa Cruz Biotechnology) were washed before use by centrifuging 50 µL of bead suspension in a 1.5 mL microcentrifuge tube for 2 min at 6000 rpm at 4 °C. The supernatant was removed and 1 mL of PBS-BSA 0.1% solution (1 mL 1× PBS + 2 µL 50 mg/ml BSA) was added. The mixture was incubated for 5 min at 4 °C on a rotating platform and then centrifuged for 2 min at 6000 rpm and 4 °C. The supernatant was subsequently removed, and the cleaning steps were repeated three more times. After cleaning, 50 µL of 1× IP buffer were added to the washed beads, and the DNA-antibody mixture was transferred to the bead mixture. The solution was incubated for 2 h at 4 °C on a rotating platform. The beads and captured DNA-antibodies were washed by centrifugation of the mixture for 2 min at 6000 rpm and 4 °C. The supernatant was removed and 1 mL of 1× IP buffer was added. The mixture was incubated for 5 min at 4 °C on a rotating platform, followed by centrifugation for 2 min at 6000 rpm and 4 °C, and removal of the supernatant. This procedure was repeated three times. To digest the beads, 210 µL of digestion buffer and 20 µL of Proteinase K (20 mg/mL) were added, followed by incubation for 2 h at 55 °C on a rotating platform. The DNA was cleaned from the bead fragments by filtering through Pierce™ Spin Cups—Paper Filter (ThermoFisher Scientific) while centrifuging at max speed for 30 s. The flow-through was collected and 3 µL of glycogen (5 mg/mL) were added. DNA was precipitated by adding 20 µL of 5 M NaCl and 750 µL ethanol, both ice cold. The solution was mixed and incubated for 30 min on ice before centrifugation at 14,000 rpm for 30 min at 4 °C. The supernatant was carefully removed and samples dried in a heating block at 50 °C for 5 min. Samples were resuspended in 30 µL of H2O and heated on a heating block at 50 °C for 5 min before measuring DNA concentrations on a Nanodrop. The DNA samples were then used for whole-genome amplification using a WGA2 kit (Sigma-Aldrich), following the manufacturer’s instructions. WGA samples were cleaned using a QIAquick PCR Purification Kit (Qiagen) following the manufacturer’s instructions, and used for next-generation sequencing.
The 12 samples of DNA that were extracted from hypothalamus tissue from six HRH and six LFH individuals and enriched for the methylated fraction were sequenced. Sequencing was performed on IonProton chips (Ion Torrent Systems, Inc) at the National Genomics Infrastructure (NGI), using the Ion Fragment Library kit according to the manufacturer’s protocol. Data generated from the sequencing was processed using the Torrent Suite (version 5.0.2) software (ThermoFisher Scientific), which is a complete analysis solution for Ion Torrent data that includes signal processing, base calling, trimming of low-quality reads, as well as alignment to a reference genome. After quality checks, the data were aligned to the Galgal4 chicken reference genome (International Chicken Genome Consortium). Two samples of DNA enriched for the methylated fraction were loaded per Ion PI chip. Thus, we used six chips for 12 samples. The raw sequencing data were uploaded to Array Express (http://www.ebi.ac.uk/arrayexpress/) under accession number E-MTAB-6407.
Pairwise sample comparisons of methylated regions between selection groups, sexes or the interaction between sex and selection were performed in R (version 3.3.1) using edgeR, integrated in the Bioconductor (release 3.3) package MEDIPS , by dividing the genome into 1-kb windows, starting at the first position on chromosome 1, and trimming away duplicated sequences. Benjamini–Hochberg correction via R’s p.adjust function was used to adjust for false discovery rate (FDR) in the multiple testing, with a cut-off of p < 0.1 used for significance. The cut-off was used to detect changes in a larger number of windows, allowing for a better overview of the effects of selection. In order to validate the results, random resampling was performed by randomizing samples into two random groups of six samples each, regardless of HFH or LFH selection, and running the entire MEDIPS analysis as above. The randomized groups were balanced regarding sex, to avoid sex differences in the results. Twenty replicated runs were carried out.
DM regions that covered an annotated gene region or its promoter, defined as starting 7.25 kb upstream of the transcriptional start site , were included in the set of regions for functional analyses. The software WebGestalt2017 (updated 1/27/2017) was used for gene ontology enrichment and KEGG pathway analysis on gene symbols obtained via Ensembl . WebGestalt uses a hypergeometric distribution for significance estimation and combines it with Bonferroni–Hochberg adjustment of p values. All known chicken genes in the chicken genome from Ensembl were used as a base for enrichment analysis. In order to identify gene network modules that were over-represented in the DM regions, the web-based ConsensusPathDB tool was used  with standard settings and allowing for intermediate nodes.
A total of 430 million sequencing reads were generated, of which more than 99% were aligned to the chicken reference genome, resulting in a 9× genome coverage per sample. Each sample generated approximately 8.5 Gb of sequence data, with an average read length of 160 bp.
Summary statistics of sequenced methylated regions in the selection lines
CpGs in windows
CpGs in significantly DM windows
Significant windows males
Significant windows females
Differentially-methylated windows between the selection lines
Window start (bp)
Adjusted p value
9.52e − 02
Adaptor related protein complex 1 sigma 2 subunit [Source:HGNC Symbol;Acc:HGNC:560]
5.87e − 02
Sodium leak channel, non-selective [Source:CGNC Symbol;Acc:12660]
9.52e − 02
5.14e − 02
Discs, large homolog 2 (Drosophila) [Source:CGNC Symbol;Acc:12965]
3.27e − 07
Calpain 5 [Source:CGNC Symbol;Acc:479]
8.12e − 02
9.52e − 02
3.22e − 02
EH domain binding protein 1 [Source:CGNC Symbol;Acc:6763]
2.73e − 02
ARFGEF family member 3 [Source:CGNC Symbol;Acc:10361]
4.83e − 03
Myotubularin related protein 9 [Source:CGNC Symbol;Acc:12476]
9.24e − 02
3.22e − 02
7.28e − 02
Gallus gallus chromosome 5 open reading frame, mRNA. [Source:RefSeq mRNA;Acc:NM_205369]
7.79e − 02
5.91e − 02
8.12e − 02
3.22e − 02
Ubiquitin specific peptidase 30 [Source:CGNC Symbol;Acc:50370]
8.12e − 02
Mitochondrial ribosomal protein L41 [Source:CGNC Symbol;Acc:6639]
8.12e − 02
Patatin-like phospholipase domain containing 7 [Source:CGNC Symbol;Acc:6632]
5.14e − 02
Armadillo repeat containing 7 [Source:CGNC Symbol;Acc:6023]
1.01e − 02
DEAD-box helicase 25 [Source:HGNC Symbol;Acc:HGNC:18698]
3.22e − 02
7.28e − 02
Guanine nucleotide binding protein-like 2 (nucleolar) [Source:CGNC Symbol;Acc:1398]
The random resampling did not result in significantly differentially-methylated windows in any of the 20 replicates [see Additional file 1: Table S1].
Sex-specific effects on DNA methylation, estimated when comparing combined HFH/LFH females against HFH/LFH males, were much larger than the effects of selection. Comparison of DNA methylation patterns between sexes, with line excluded as a factor, revealed many differences, with over 2500 DM windows on the autosomes and almost 2000 on the sex chromosomes W and Z [see Additional file 2: Figure S1]. A more detailed analysis of the DM regions indicated that 79% of the windows were hypermethylated in males.
Successive selection of HFH and LFH birds seemed to result in large sex-specific effects in DNA methylation differences, with 66 windows being DM between HFH and LFH females (FDR adjusted p < 0.1) [see Additional file 3: Table S2]. Moreover, 51 DM windows were identified between HFH and LFH males (FDR adjusted p < 0.1). Comparison of DM windows identified within males and within females showed overlaps at only four windows, three of which were located on chromosome 7 and covered the expressed sequence tag (EST) Gga.15462, which was previously detected in testis mRNA in RJF and White Leghorn . The methylation patterns for this EST region differed between selection lines and between sexes, i.e. it was highly methylated in HFH males and in LFH females. The fourth window was located on chromosome 1 in an intron of the CAPN5 gene and was hypermethylated in HFH animals of both sexes.
Enrichment results for gene ontology (GO) analysis of differentially-methylated regions between selection lines across and within sexes
Cation channel activity
Calcium ion transport
Regulation of membrane potential
Passive transmembrane transporter activity
Amino acid binding
Glutamate receptor activity
GABA receptor complex
Response to histamine
Chloride channel activity
DNA-dependent ATPase activity
We examined the overlap between the DM windows identified here and the genomic regions that included genes that were previously reported as differentially expressed (DE) in the hypothalamus between the same selection groups and generation . Genomic regions located 7.25 kb upstream from DE transcripts were evaluated in order to include promoter regions . No direct overlaps were detected between DM and DE regions, even when including DM regions with a fold change higher than 1.5 and p < 0.05 before FDR correction. ConsensusPathDB analysis of these regions revealed that the DM genes MRPL41 and ARMC7 (Table 2) were significantly connected with the DE genes MAEA and ANKRD1 that were detected by Bélteky et al.  through one intermediate node. DM windows within males and within females were also compared to previously reported DE genes, but neither sex showed any overlap between expression and methylation changes .
In this experiment, we studied the effects of bidirectional selection for tameness on the DNA methylome of the hypothalamus in RJF. After five generations of selection for high and low fear of humans, 22 genomic regions with differential DNA methylation in the hypothalamus were detected. The DM regions targeted several genes with functions that may be related to previously reported phenotypic selection responses [11, 36], such as metabolism and signaling. However, the number of significant windows before correcting for multiple-testing was much larger than 22 windows, with over 51,000 differently methylated windows observed at p > 0.05. This large reduction in significant windows after correcting for FDR may indicate that extensive variation exists within the groups for methylation patterns in regions that are not strongly affected by the selection process. Such naturally occurring inter-individual variation has already been reported in humans .
Among the 22 DM windows, none overlapped with DE genes reported for the same tissue of birds from the same populations, which agrees with other studies on RJF and domesticated White Leghorn in which no overlap between DM promoters and gene expression of the same genes was observed. This suggests that other epigenetic mechanisms are involved in regulating gene expression , or that DNA methylation changes in regions outside of promoter and gene regions, i.e. in intergenic regions, could affect transcription in a trans-acting fashion [50, 51], or that the DM regions affect genes that are not activated by any pathways at the life stage studied here. Furthermore, gene expression may be affected primarily by methylation changes of single CpG sites within a CpG-dense region, which may not be detected by MeDIP-seq .
Although we observed no direct overlap between DM and DE genes, GO analysis of genes included in DM windows revealed functions similar to those of the DE genes in the same group of animals. Interestingly, some DM regions were related to behavioral and metabolic pathways. Further research should confirm whether and how these pathways are affected by selection for tameness. However, based on the genes targeted by the DM regions between HFH and LFH animals, we can derive some hints. For example, calpain 5 (CAPN5), which encodes a cytosolic cysteine protease that acts on signaling-related molecules and is involved in cell differentiation and proliferation , was reported to interact with nuclear receptors and thus to impact metabolism . Mutations in the gene sodium leak channel non-selective (NALCN), which is involved in the control of neuronal excitability, can lead to speech impairment and intellectual disability in humans [55, 56]. The gene DLG2, which encodes a scaffold protein belonging to the membrane-associated guanylate kinase (MAGUK) family and is active at postsynaptic sites, is associated with neurodevelopmental disorders, schizophrenia and cognition in humans [57, 58]. Downregulation of the mitochondrial gene deubiquitinase USP30 is reported to enhance degradation of damage mitochondria in human neurons, which is beneficial during neurological disorders . The fact that several of the genes related to DM regions are associated with neuronal functions suggests that molecular mechanisms involved in the behavioral changes have emerged between HFH and LFH birds .
We compared the DM regions detected here with DM probes that were previously detected in the hypothalamus by Nätt et al.  between RJF and the domesticated White Leghorn, in order to test whether the same genes with a modified methylation status were found in the current experiment and during the domestication of White Leghorn. We did not identify any common genes between these two selection processes, which could be explained by the fact that the driver of domestication of White Leghorn was mainly selection for egg size, while in our selection lines the main driver was tameness. This suggests that different selection pressures generate distinctive sets of epigenetic changes, which in turn are related to specific phenotypic traits. This concept can also apply to the transcriptome, since we observed no concordance between previously published DE genes from the same birds as those used in this study and DE genes observed in the hypothalamus of RJF and White Leghorn [28, 38].
A possible explanation for the lack of overlap between DM regions in the current and previous studies of the chicken hypothalamus is that the observed differences may be related to differential genetic drift in HFH and LFH populations . However, previous studies on gene expression differences in the same populations, along with an unselected line of RJF, indicated that the expression changes that occurred between the parental generation and the fifth selected generations resulted from the imposed selection and not from genetic drift [38, 39]. These results, together with the random resampling performed in the current experiment, strengthen the argument that the epigenetic changes observed are a consequence of the artificial selection imposed within a rather short period of time (5 generations), and not of genetic drift.
In addition to the DM regions detected between the HFH and LFH lines in this experiment, we found numerous sex-specific changes, which are consistent with previously reported results . However, in our study, differences were larger than observed previously, because it covered the entire genome, in contrast to only the promoter regions in a previous study that analyzed DNA methylation using promoter arrays . The largest sex-specific DM regions were associated to the genes ZFR on chromosome 1 and MHM on chromosome Z, which were both highly methylated in males. This large number of sex-specific differences in DNA-methylation suggests that the two sexes respond differently to selection pressures related to behavioral traits, in agreement with previous findings , but raising some interesting questions about underlying mechanisms. An interesting finding was that many of these sex-specific DNA methylation changes are on autosomes, which suggests that DNA methylation of autosomes plays a more important role than expected in gender-specific characteristics in the brain of vertebrates.
We detected 22 DM regions by comparing hypothalamic DNA from RJF selected during five generations for high versus low fear of humans. Functional annotation of the genes associated with these DM regions showed that they are related to, for example, cellular metabolism and neural signaling, similar to what was previously reported in terms of gene expression differences for the same animals. Our results suggest that bidirectional selection for tameness involves epigenetic factors that can even differ in a sex-specific manner. Observation of divergent DNA methylation patterns in the hypothalamus after only five generations of artificial selection highlights the importance of epigenetic mechanisms, in addition to genetic composition, in evolutionary phenotypic variation that emerges in response to selection pressures. Future research should delve into the molecular mechanism involved in the emergence of somatic epigenetic differences during selection. Two options could explain this. One is that epigenetic differences are linked to genetic differences emerging during selection, as suggested by Verhulst et al. . Another possibility is that selection on phenotypes would concomitantly select specific germ line epigenomes. Different germ line epigenomes in divergent selection lines would then influence somatic epigenomes later on in the ontogeny of each individual. Interestingly, we detected several sex-specific epigenetic changes on the autosomes, which raises the question about whether epigenetic differences have a role in gender-specific behavioral responses that are unrelated to sex chromosomes.
All authors contributed extensively to the work presented in this paper. JB, BA and PJ designed the experiment. BA supervised breeding of the selection lines. JB, LB, AH and CGB analyzed and interpreted the sequencing data. JB administered the experiment and wrote the first version of the manuscript. JB and CGB supervised the final version. All authors read and approved the final manuscript.
The authors would like to acknowledge support of the National Genomics Infrastructure (NGI)/Uppsala Genome Center and UPPMAX for providing assistance in massive parallel sequencing and computational infrastructure. Work performed at NGI/Uppsala Genome Center has been funded by RFI/VR and Science for Life Laboratory, Sweden. We also extend thanks to Petros Batakis for sample preparations, and special thanks to Markus Jöngren for help in developing the breeding programme.
The authors declare that they have no competing interests.
Availability of data and materials
The raw sequencing data have been uploaded to Array Express (http://www.ebi.ac.uk/arrayexpress/) under accession number E-MTAB-6407.
Consent for publication
Ethics approval and consent to participate
The experiments reported in this paper were carried out in accordance with ethical guidelines, and were approved by the Linköping Animal Ethics Committee, license no 122-10.
The work was supported by grants from the research council Formas, Vetenskapsrådet, and the European Research Council, ERC (Grant No. 322206 GENEWELL).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Zeder MA. The domestication of animals. J Anthropol Res. 2012;68:161–90.View ArticleGoogle Scholar
- Price EO. Animal domestication and behavior. Wallingford: CABI Publishing; 2002.View ArticleGoogle Scholar
- Sánchez-Villagra MR, Geiger M, Schneider RA. The taming of the neural crest: a developmental perspective on the origins of morphological covariation in domesticated mammals. R Soc Open Sci. 2016;3:160107.View ArticlePubMedPubMed CentralGoogle Scholar
- Rubin CJ, Zody MC, Eriksson J, Meadows JR, Sherwood E, Webster MT, et al. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010;464:587–91.View ArticlePubMedGoogle Scholar
- Andersson L, Georges M. Domestic-animal genomics: deciphering the genetics of complex traits. Nat Rev Genet. 2004;5:202–12.View ArticlePubMedGoogle Scholar
- Heyne HO, Lautenschlager S, Nelson R, Besnier F, Rotival M, Cagan A, et al. Genetic influences on brain gene expression in rats selected for tameness and aggression. Genetics. 2014;198:1277–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Kukekova AV, Johnson JL, Teiling C, Li L, Oskina IN, Kharlamova AV, et al. Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes). BMC Genomics. 2011;12:482.View ArticlePubMedPubMed CentralGoogle Scholar
- Trut L, Oskina I, Kharlamova A. Animal evolution during domestication: the domesticated fox as a model. BioEssays. 2009;31:349–60.View ArticlePubMedPubMed CentralGoogle Scholar
- The Belyaev DK, Wilhelmine E. Key 1978 invitational lecture. Destabilizing selection as a factor in domestication. J Hered. 1979;70:301–8.View ArticleGoogle Scholar
- Albert FW, Shchepina O, Winter C, Rompler H, Teupser D, Palme R, et al. Phenotypic differences in behavior, physiology and neurochemistry between rats selected for tameness and for defensive aggression towards humans. Horm Behav. 2008;53:413–21.View ArticlePubMedGoogle Scholar
- Agnvall B, Katajamaa R, Altimiras J, Jensen P. Is domestication driven by reduced fear of humans? Boldness, metabolism and serotonin levels in divergently selected red junglefowl (Gallus gallus). Biol Lett. 2015;11:20150509.View ArticlePubMedPubMed CentralGoogle Scholar
- Jensen P. Adding ‘epi-’ to behaviour genetics: implications for animal domestication. J Exp Biol. 2015;218:32–40.View ArticlePubMedGoogle Scholar
- Verhulst EC, Mateman AC, Zwier MV, Caro SP, Verhoeven KJF, Van Oers K, et al. Evidence from pyrosequencing indicates that natural variation in animal personality is associated with DRD4 DNA methylation. Mol Ecol. 2016;25:1801–11.View ArticlePubMedGoogle Scholar
- Weaver IC, Cervoni N, Champagne FA, D’Alessio AC, Sharma S, Seckl JR, et al. Epigenetic programming by maternal behavior. Nat Neurosci. 2004;7:847–54.View ArticlePubMedGoogle Scholar
- Tian F, Zhan F, VanderKraats ND, Hiken JF, Edwards JR, Zhang H, et al. DNMT gene expression and methylome in Marek’s disease resistant and susceptible chickens prior to and following infection by MDV. Epigenetics. 2013;8:431–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Berghof TVL, Parmentier HK, Lammers A. Transgenerational epigenetic effects on innate immunity in broilers: an underestimated field to be explored? Poult Sci. 2013;92:2904–13.View ArticlePubMedGoogle Scholar
- Hu Y, Xu H, Li Z, Zheng X, Jia X, Nie Q, et al. Comparison of the genome-wide DNA methylation profiles between fast-growing and slow-growing broilers. PLoS One. 2013;8:e56411.View ArticlePubMedPubMed CentralGoogle Scholar
- Pertille F, Brantsaeter M, Nordgreen J, Coutinho LL, Janczak AM, Jensen P, et al. DNA methylation profiles in red blood cells of adult hens correlate to their rearing conditions. J Exp Biol. 2017;220:3579–87.View ArticlePubMedGoogle Scholar
- Raynal NJ, Si J, Taby RF, Gharibyan V, Ahmed S, Jelinek J, et al. DNA methylation does not stably lock gene expression but instead serves as a molecular mark for gene silencing memory. Cancer Res. 2012;72:1170–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Guerrero-Bosagna C, Settles M, Lucker B, Skinner MK. Epigenetic transgenerational actions of vinclozolin on promoter regions of the sperm epigenome. PLoS One. 2010;5:e13100.View ArticlePubMedPubMed CentralGoogle Scholar
- Lindqvist C, Janczak AM, Natt D, Baranowska I, Lindqvist N, Wichman A, et al. Transmission of stress-induced learning impairment and associated brain gene expression from parents to offspring in chickens. PLoS One. 2007;2:e364.View ArticlePubMedPubMed CentralGoogle Scholar
- Goerlich VC, Natt D, Elfwing M, Macdonald B, Jensen P. Transgenerational effects of early experience on behavioral, hormonal and gene expression responses to acute stress in the precocial chicken. Horm Behav. 2012;61:711–8.View ArticlePubMedGoogle Scholar
- Gapp K, Soldado-Magraner S, Alvarez-Sánchez M, Bohacek J, Vernaz G, Shu H, et al. Early life stress in fathers improves behavioural flexibility in their offspring. Nat Commun. 2014;5:5466.View ArticlePubMedGoogle Scholar
- Franklin TB, Russig H, Weiss IC, Gräff J, Linder N, Michalon A, et al. Epigenetic transmission of the impact of early stress across generations. Biol Psychiatry. 2010;68:408–15.View ArticlePubMedGoogle Scholar
- Nätt D, Lindqvist N, Stranneheim H, Lundeberg J, Torjesen PA, Jensen P. Inheritance of acquired behaviour adaptations and brain gene expression in chickens. PLoS One. 2009;4:e6405.View ArticlePubMedPubMed CentralGoogle Scholar
- Janowitz Koch I, Clark MM, Thompson MJ, Deere-Machemer KA, Wang J, Duarte L, et al. The concerted impact of domestication and transposon insertions on methylation patterns between dogs and grey wolves. Mol Ecol. 2016;25:1838–55.View ArticlePubMedGoogle Scholar
- Xiang H, Li X, Dai F, Xu X, Tan A, Chen L, et al. Comparative methylomics between domesticated and wild silkworms implies possible epigenetic influences on silkworm domestication. BMC Genomics. 2013;14:646.View ArticlePubMedPubMed CentralGoogle Scholar
- Nätt D, Rubin CJ, Wright D, Johnsson M, Bélteky J, Andersson L, et al. Heritable genome-wide variation of gene expression and promoter methylation between wild and domesticated chickens. BMC Genomics. 2012;13:59.View ArticlePubMedPubMed CentralGoogle Scholar
- Skinner MK, Gurerrero-Bosagna C, Haque MM, Nilsson EE, Koop JAH, Knutie SA, et al. Epigenetics and the evolution of Darwin’s finches. Genome Biol Evol. 2014;6:1972–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Niu D, Fu Y, Luo J, Ruan H, Yu XP, Chen G, et al. The origin and genetic diversity of Chinese native chicken breeds. Biochem Genet. 2002;40:163–74.View ArticlePubMedGoogle Scholar
- Zuidhof MJ, Schneider BL, Carney VL, Korver DR, Robinson FE. Growth, efficiency, and yield of commercial broilers from 1957, 1978, and 2005. Poult Sci. 2014;93:2970–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Nicol CJ. The behavioural biology of ckickens. Wallingford: CABI; 2015.View ArticleGoogle Scholar
- Schutz KE, Jensen P. Effects of resource allocation on behavioural strategies: a comparison of red junglefowl (Gallus gallus) and two domesticated breeds of poultry. Ethology. 2001;107:753–65.View ArticleGoogle Scholar
- Vaisanen J, Jensen P. Social versus exploration and foraging motivation in young red junglefowl (Gallus gallus) and White Leghorn layers. Appl Anim Behav Sci. 2003;84:139–58.View ArticleGoogle Scholar
- Lindgren I, Altimiras J. Sensitivity of organ growth to chronically low oxygen levels during incubation in red junglefowl and domesticated chicken breeds. Poult Sci. 2011;90:126–35.View ArticlePubMedGoogle Scholar
- Agnvall B, Jöngren M, Strandberg E, Jensen P. Heritability and genetic correlations of fear-related behaviour in red junglefowl–possible implications for early domestication. PLoS One. 2012;7:e35162.View ArticlePubMedPubMed CentralGoogle Scholar
- Agnvall B, Ali A, Olby S, Jensen P. Red Junglefowl (Gallus gallus) selected for low fear of humans are larger, more dominant and produce larger offspring. Animal. 2014;8:1498–505.View ArticlePubMedGoogle Scholar
- Bélteky J, Agnvall B, Johnsson M, Wright D, Jensen P. Domestication and tameness: brain gene expression in red junglefowl selected for less fear of humans suggests effects on reproduction and immunology. R Soc Open Sci. 2016;3:160033.View ArticlePubMedPubMed CentralGoogle Scholar
- Bélteky J, Agnvall B, Jensen P. Gene expression of behaviorally relevant genes in the cerebral hemisphere changes after selection for tameness in Red Junglefowl. PLoS One. 2017;12:e0177004.View ArticlePubMedPubMed CentralGoogle Scholar
- Weber M, Davies JJ, Wittig D, Oakeley EJ, Haase M, Lam WL, et al. Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. Nat Genet. 2005;37:853–62.View ArticlePubMedGoogle Scholar
- Taiwo O, Wilson GA, Morris T, Seisenberger S, Reik W, Pearce D, et al. Methylome analysis using MeDIP-seq with low DNA concentrations. Nat Protoc. 2012;7:617–36.View ArticlePubMedGoogle Scholar
- Bock C, Tomazou EM, Brinkman AB, Muller F, Simmer F, Gu H, et al. Quantitative comparison of genome-wide DNA methylation mapping technologies. Nat Biotechnol. 2010;28:1106–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Guerrero-Bosagna C, Jensen P. Optimized method for methylated DNA immuno-precipitation. MethodsX. 2015;2:432–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.View ArticlePubMedGoogle Scholar
- McLaren W, Pritchard B, Rios D, Chen Y, Flicek P, Cunningham F. Deriving the consequences of genomic variants with the Ensembl API and SNP effect predictor. Bioinformatics. 2010;26:2069–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Kamburov A, Pentchev K, Galicka H, Wierling C, Lehrach H, Herwig R. ConsensusPathDB: toward a more complete picture of cell biology. Nucleic Acids Res. 2011;39:D712–7.View ArticlePubMedGoogle Scholar
- Skinner MK, Guerrero-Bosagna C. Role of CpG deserts in the epigenetic transgenerational inheritance of differential DNA methylation regions. BMC Genomics. 2014;15:692.View ArticlePubMedPubMed CentralGoogle Scholar
- Savolainen P, Fitzsimmons C, Arvestad L, Andersson L, Lundeberg J. ESTs from brain and testis of White Leghorn and red junglefowl: annotation, bioinformatic classification of unknown transcripts and analysis of expression levels. Cytogenet Genome Res. 2005;111:79–87.View ArticlePubMedGoogle Scholar
- Rakyan VK, Hildmann T, Novik KL, Lewin J, Tost J, Cox AV, et al. DNA methylation profiling of the human major histocompatibility complex: a pilot study for the human epigenome project. PLoS Biol. 2004;2:e405.View ArticlePubMedPubMed CentralGoogle Scholar
- Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, Souza C, Fouse SD, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010;466:253–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith ZD, Chan MM, Mikkelsen TS, Gu H, Gnirke A, Regev A, et al. A unique regulatory phase of DNA methylation in the early mammalian embryo. Nature. 2012;484:339–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Jeong HM, Lee S, Chae H, Kim R, Kwon MJ, Oh E, et al. Efficiency of methylated DNA immunoprecipitation bisulphite sequencing for whole-genome DNA methylation analysis. Epigenomics. 2016;8:1061–77.View ArticlePubMedGoogle Scholar
- Sato K, Kawashima S. Calpain function in the modulation of signal transduction molecules. Biol Chem. 2001;382:743–51.View ArticlePubMedGoogle Scholar
- Sáez ME, Grilo A, Morón FJ, Manzano L, Martínez-Larrad MT, González-Pérez A, et al. Interaction between calpain 5, peroxisome proliferator-activated receptor-gamma and peroxisome proliferator-activated receptor-delta genes: a polygenic approach to obesity. Cardiovasc Diabetol. 2008;7:23.View ArticlePubMedPubMed CentralGoogle Scholar
- Al-Sayed Moeenaldeen D, Al-Zaidan H, Albakheet A, Hakami H, Kenana R, Al-Yafee Y, et al. Mutations in NALCN cause an autosomal-recessive syndrome with severe hypotonia, speech impairment, and cognitive delay. Am J Hum Genet. 2013;93:721–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Köroğlu Ç, Seven M, Tolun A. Recessive truncating NALCN mutation in infantile neuroaxonal dystrophy with facial dysmorphism. J Med Genet. 2013;50:515–20.View ArticlePubMedGoogle Scholar
- Ingason A, Giegling I, Hartmann AM, Genius J, Konte B, Friedl M, et al. Expression analysis in a rat psychosis model identifies novel candidate genes validated in a large case-control sample of schizophrenia. Transl Psychiatry. 2015;5:e656.View ArticlePubMedPubMed CentralGoogle Scholar
- Nithianantharajah J, Komiyama NH, McKechanie A, Johnstone M, Blackwood DH, Clair DS, et al. Synaptic scaffold evolution generated components of vertebrate cognitive complexity. Nat Neurosci. 2013;16:16–24.View ArticlePubMedGoogle Scholar
- Bingol B, Tea JS, Phu L, Reichelt M, Bakalarski CE, Song Q, et al. The mitochondrial deubiquitinase USP30 opposes parkin-mediated mitophagy. Nature. 2014;510:370–5.View ArticlePubMedGoogle Scholar
- Agnvall B, Jensen P. Effects of divergent selection for fear of humans on behaviour in Red Junglefowl. PLoS One. 2016;11:e0166075.View ArticlePubMedPubMed CentralGoogle Scholar
- Nätt D, Agnvall B, Jensen P. Large sex differences in chicken behavior and brain gene expression coincide with few differences in promoter DNA-methylation. PLoS One. 2014;9:e96376.View ArticlePubMedPubMed CentralGoogle Scholar