Gene and transposable element methylation in great tit (Parus major) brain and blood

Manipulation of a social signal affects DNA methylation of a stress-related gene in a free-living bird

Social status directly affects the health of humans and other animals. Low status individuals receive more antagonistic encounters, have fewer supportive relationships and have worse health outcomes. However, the physiological and cellular processes that mediate the relationship between the social environment and health are incompletely known. Epigenetic regulation of the hypothalamic-pituitary-adrenal (HPA) axis, the neuroendocrine pathway that activates in response to stressors, may be one process that is sensitive to the social environment. Here, we experimentally manipulated plumage, a key social signal in female tree swallows (Tachycineta bicolor) and quantified methylation of four genes in the HPA axis before and after treatment. We found that dulling the white breast plumage affected methylation in one gene, CRHR1; however, the effect depended on the original brightness of the bird. Methylation in this gene was correlated with baseline corticosterone levels, suggesting that DNA methylation of CRHR1 helps regulate glucocorticoid production in this species. Methylation in two other genes, FKBP5 and GR, changed over the course of the experiment, independent of treatment. These results show that methylation of these genes is labile into adulthood and suggest that epigenetic regulation of the HPA axis could help birds respond to current environmental conditions.

Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

The genome of a globally invasive passerine, the common myna, Acridotheres tristis

In an era of global climate change, biodiversity conservation is receiving increased attention. Conservation efforts are greatly aided by genetic tools and approaches, which seek to understand patterns of genetic diversity and how they impact species health and their ability to persist under future climate regimes. Invasive species offer vital model systems in which to investigate questions regarding adaptive potential, with a particular focus on how changes in genetic diversity and effective population size interact with novel selection regimes. The common myna (Acridotheres tristis) is a globally invasive passerine and is an excellent model species for research both into the persistence of low-diversity populations and the mechanisms of biological invasion. To underpin research on the invasion genetics of this species, we present the genome assembly of the common myna. We describe the genomic landscape of this species, including genome wide allelic diversity, methylation, repeats, and recombination rate, as well as an examination of gene family evolution. Finally, we use demographic analysis to identify that some native regions underwent a dramatic population increase between the two most recent periods of glaciation, and reveal artefactual impacts of genetic bottlenecks on demographic analysis.

An ecologist's guide for studying DNA methylation variation in wild vertebrates

The field of molecular biology is advancing fast with new powerful technologies, sequencing methods and analysis software being developed constantly. Commonly used tools originally developed for research on humans and model species are now regularly used in ecological and evolutionary research. There is also a growing interest in the causes and consequences of epigenetic variation in natural populations. Studying ecological epigenetics is currently challenging, especially for vertebrate systems, because of the required technical expertise, complications with analyses and interpretation, and limitations in acquiring sufficiently high sample sizes. Importantly, neglecting the limitations of the experimental setup, technology and analyses may affect the reliability and reproducibility, and the extent to which unbiased conclusions can be drawn from these studies. Here, we provide a practical guide for researchers aiming to study DNA methylation variation in wild vertebrates. We review the technical aspects of epigenetic research, concentrating on DNA methylation using bisulfite sequencing, discuss the limitations and possible pitfalls, and how to overcome them through rigid and reproducible data analysis. This review provides a solid foundation for the proper design of epigenetic studies, a clear roadmap on the best practices for correct data analysis and a realistic view on the limitations for studying ecological epigenetics in vertebrates. This review will help researchers studying the ecological and evolutionary implications of epigenetic variation in wild populations.

The Epigenetics of Animal Personality

Animal personality, consistent individual differences in behaviour, is an important concept for understanding how individuals vary in how they cope with environmental challenges. In order to understand the evolutionary significance of animal personality, it is crucial to understand the underlying regulatory mechanisms. Epigenetic marks such as DNA methylation are hypothesised to play a major role in explaining variation in phenotypic changes in response to environmental alterations. Several characteristics of DNA methylation also align well with the concept of animal personality. In this review paper, we summarise the current literature on the role that molecular epigenetic mechanisms may have in explaining personality variation. We elaborate on the potential for epigenetic mechanisms to explain behavioural variation, behavioural development and temporal consistency in behaviour. We then suggest future routes for this emerging field and point to potential pitfalls that may be encountered. We conclude that a more inclusive approach is needed for studying the epigenetics of animal personality and that epigenetic mechanisms cannot be studied without considering the genetic background.

Detailed molecular and epigenetic characterization of the pig IPEC-J2 and chicken SL-29 cell lines

The pig IPEC-J2 and chicken SL-29 cell lines are of interest because of their untransformed nature and wide use in functional studies. Molecular characterization of these cell lines is important to gain insight into possible molecular aberrations. The aim of this paper is to provide a molecular and epigenetic characterization of the IPEC-J2 and SL-29 cell lines, a cell-line reference for the FAANG community, and future biomedical research. Whole genome sequencing, gene expression, DNA methylation, chromatin accessibility, and ChIP-seq of four histone marks (H3K4me1, H3K4me3, H3K27ac, H3K27me3) and an insulator (CTCF) are used to achieve these aims. Heteroploidy (aneuploidy) of various chromosomes was observed from whole genome sequencing analysis in both cell lines. Furthermore, higher gene expression for genes located on chromosomes with aneuploidy in comparison to diploid chromosomes was observed. Regulatory complexity of gene expression, DNA methylation, and chromatin accessibility was investigated through an integrative approach.

Changes in the rearing environment cause reorganization of molecular networks associated with DNA methylation

Disentangling the interaction between the genetic basis and environmental context underlying phenotypic variation is critical for understanding organismal evolution. Environmental change, such as increased rates of urbanization, can induce shifts in phenotypic plasticity with some individuals adapting to city life while others are displaced. A key trait that can facilitate adaptation is the degree at which animals respond to stressors. This stress response, which includes elevation of baseline circulating concentrations of glucocorticoids, has a heritable component and exhibits intra- and inter-individual variation. However, the mechanisms behind this variability and whether they might be responsible for adaptation to different environments are not known. Variation in DNA methylation can be a potential mechanism that mediates environmental effects on the stress response, as early-life stressors increase glucocorticoid concentrations and change adult phenotype. We used an inter- and intra-environmental cross-foster experiment to analyse the contribution of DNA methylation to early-life phenotypic variation. We found that at hatching, urban house wren (Troglodytes aedon) offspring had higher methylation frequencies compared with their rural counterparts. We also observed age-related patterns in offspring methylation, indicating the developmental effects of the rearing environment on methylation. At fledgling, differential methylation analyses showed that cellular respiration genes were differentially methylated in broods of different origins and behavioural and metabolism genes were differentially methylated in broods of different rearing environments. Lastly, hyper-methylation of a single gene (CNTNAP2) is associated with decreased glucocorticoid levels and the rearing environment. These differential methylation patterns linked to a specific physiological phenotype suggest that DNA methylation may be a mechanism by which individuals adjust to novel environments during their lifespan. Characterizing genetic and environmental influences on methylation is critical for understanding the role of epigenetic mechanisms in evolutionary adaptation.

Epigenetics and the city: Non-parallel DNA methylation modifications across pairs of urban-forest Great tit populations

Identifying the molecular mechanisms involved in rapid adaptation to novel environments and determining their predictability are central questions in evolutionary biology and pressing issues due to rapid global changes. Complementary to genetic responses to selection, faster epigenetic variations such as modifications of DNA methylation may play a substantial role in rapid adaptation. In the context of rampant urbanization, joint examinations of genomic and epigenomic mechanisms are still lacking. Here, we investigated genomic (SNP) and epigenomic (CpG methylation) responses to urban life in a passerine bird, the Great tit (Parus major). To test whether urban evolution is predictable (i.e. parallel) or involves mostly nonparallel molecular processes among cities, we analysed both SNP and CpG methylation variations across three distinct pairs of city and forest Great tit populations in Europe. Our analyses reveal a polygenic response to urban life, with both many genes putatively under weak divergent selection and multiple differentially methylated regions (DMRs) between forest and city great tits. DMRs mainly overlapped transcription start sites and promotor regions, suggesting their importance in modulating gene expression. Both genomic and epigenomic outliers were found in genomic regions enriched for genes with biological functions related to the nervous system, immunity, or behavioural, hormonal and stress responses. Interestingly, comparisons across the three pairs of city-forest populations suggested little parallelism in both genetic and epigenetic responses. Our results confirm, at both the genetic and epigenetic levels, hypotheses of polygenic and largely nonparallel mechanisms of rapid adaptation in novel environments such as urbanized areas.

Performance of methods to detect genetic variants from bisulphite sequencing data in a non-model species

The profiling of epigenetic marks like DNA methylation has become a central aspect of studies in evolution and ecology. Bisulphite sequencing is commonly used for assessing genome‐wide DNA methylation at single nucleotide resolution but these data can also provide information on genetic variants like single nucleotide polymorphisms (SNPs). However, bisulphite conversion causes unmethylated cytosines to appear as thymines, complicating the alignment and subsequent SNP calling. Several tools have been developed to overcome this challenge, but there is no independent evaluation of such tools for non‐model species, which often lack genomic references. Here, we used whole‐genome bisulphite sequencing (WGBS) data from four female great tits (Parus major) to evaluate the performance of seven tools for SNP calling from bisulphite sequencing data. We used SNPs from whole‐genome resequencing data of the same samples as baseline SNPs to assess common performance metrics like sensitivity, precision, and the number of true positive, false positive, and false negative SNPs for the full range of variant and genotype quality values. We found clear differences between the tools in either optimizing precision (bis‐snp), sensitivity (biscuit), or a compromise between both (all other tools). Overall, the choice of SNP caller strongly depends on which performance parameter should be maximized and whether ascertainment bias should be minimized to optimize downstream analysis, highlighting the need for studies that assess such differences.

The effect of experimental lead pollution on DNA methylation in a wild bird population

Anthropogenic pollution is known to negatively influence an organism's physiology, behaviour, and fitness. Epigenetic regulation, such as DNA methylation, has been hypothesized as a potential mechanism to mediate such effects, yet studies in wild species are lacking. We first investigated the effects of early-life exposure to the heavy metal lead (Pb) on DNA methylation levels in a wild population of great tits (Parus major), by experimentally exposing nestlings to Pb at environmentally relevant levels. Secondly, we compared nestling DNA methylation from a population exposed to long-term heavy metal pollution (close to a copper smelter), where birds suffer from pollution-related decrease in food quality, and a control population. For both comparisons, the analysis of about one million CpGs covering most of the annotated genes revealed that pollution-related changes in DNA methylation were not genome wide, but enriched for genes underlying developmental processes. However, the results were not consistent when using binomial or beta binomial regression highlighting the difficulty of modelling variance in CpGs. Our study indicates that post-natal anthropogenic heavy metal exposure can affect methylation levels of development related genes in a wild bird population.

Does Arsenic Contamination Affect DNA Methylation Patterns in a Wild Bird Population? An Experimental Approach

Pollutants, such as toxic metals, negatively influence organismal health and performance, even leading to population collapses. Studies in model organisms have shown that epigenetic marks, such as DNA methylation, can be modulated by various environmental factors, including pollutants, influencing gene expression, and various organismal traits. Yet experimental data on the effects of pollution on DNA methylation from wild animal populations are largely lacking. We here experimentally investigated for the first time the effects of early-life exposure to environmentally relevant levels of a key pollutant, arsenic (As), on genome-wide DNA methylation in a wild bird population. We experimentally exposed nestlings of great tits (Parus major) to arsenic during their postnatal developmental period (3 to 14 days post-hatching) and compared their erythrocyte DNA methylation levels to those of respective controls. In contrast to predictions, we found no overall hypomethylation in the arsenic group. We found evidence for loci to be differentially methylated between the treatment groups, but for five CpG sites only. Three of the sites were located in gene bodies of zinc finger and BTB domain containing 47 (ZBTB47), HIVEP zinc finger 3 (HIVEP3), and insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1). Further studies are needed to evaluate whether epigenetic dysregulation is a commonly observed phenomenon in polluted populations and what are the consequences for organism functioning and for population dynamics.

Does Arsenic Contamination Affect DNA Methylation Patterns in a Wild Bird Population? An Experimental Approach

Pollutants, such as toxic metals, negatively influence organismal health and performance, even leading to population collapses. Studies in model organisms have shown that epigenetic marks, such as DNA methylation, can be modulated by various environmental factors, including pollutants, influencing gene expression, and various organismal traits. Yet experimental data on the effects of pollution on DNA methylation from wild animal populations are largely lacking. We here experimentally investigated for the first time the effects of early-life exposure to environmentally relevant levels of a key pollutant, arsenic (As), on genome-wide DNA methylation in a wild bird population. We experimentally exposed nestlings of great tits (Parus major) to arsenic during their postnatal developmental period (3 to 14 days post-hatching) and compared their erythrocyte DNA methylation levels to those of respective controls. In contrast to predictions, we found no overall hypomethylation in the arsenic group. We found evidence for loci to be differentially methylated between the treatment groups, but for five CpG sites only. Three of the sites were located in gene bodies of zinc finger and BTB domain containing 47 (ZBTB47), HIVEP zinc finger 3 (HIVEP3), and insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1). Further studies are needed to evaluate whether epigenetic dysregulation is a commonly observed phenomenon in polluted populations and what are the consequences for organism functioning and for population dynamics.

Epigenetics and Early Life Stress: Experimental Brood Size Affects DNA Methylation in Great Tits (Parus major)

Early developmental conditions are known to have life-long effects on an individual’s behavior, physiology and fitness. In altricial birds, a majority of these conditions, such as the number of siblings and the amount of food provisioned, are controlled by the parents. This opens up the potential for parents to adjust the behavior and physiology of their offspring according to local post-natal circumstances. However, the mechanisms underlying such intergenerational regulation remain largely unknown. A mechanism often proposed to possibly explain how parental effects mediate consistent phenotypic change is DNA methylation. To investigate whether early life effects on offspring phenotypes are mediated by DNA methylation, we cross-fostered great tit (Parus major) nestlings and manipulated their brood size in a natural study population. We assessed genome-wide DNA methylation levels of CpG sites in erythrocyte DNA, using Reduced Representation Bisulfite Sequencing (RRBS). By comparing DNA methylation levels between biological siblings raised in enlarged and reduced broods and between biological siblings of control broods, we assessed which CpG sites were differentially methylated due to brood size. We found 32 differentially methylated sites (DMS) between siblings from enlarged and reduced broods, a larger number than in the comparison between siblings from control broods. A considerable number of these DMS were located in or near genes involved in development, growth, metabolism, behavior and cognition. Since the biological functions of these genes line up with previously found effects of brood size and food availability, it is likely that the nestlings in the enlarged broods suffered from nutritional stress. We therefore conclude that early life stress might directly affect epigenetic regulation of genes related to early life conditions. Future studies should link such experimentally induced DNA methylation changes to expression of phenotypic traits and assess whether these effects affect parental fitness to determine if such changes are also adaptive.

The emergence of the brain non-CpG methylation system in vertebrates

Mammalian brains feature exceptionally high levels of non-CpG DNA methylation alongside the canonical form of CpG methylation. Non-CpG methylation plays a critical regulatory role in cognitive function, which is mediated by the binding of MeCP2, the transcriptional regulator that when mutated causes Rett syndrome. However, it is unclear whether the non-CpG neural methylation system is restricted to mammalian species with complex cognitive abilities or has deeper evolutionary origins. To test this, we investigated brain DNA methylation across 12 distantly related animal lineages, revealing that non-CpG methylation is restricted to vertebrates. We discovered that in vertebrates, non-CpG methylation is enriched within a highly conserved set of developmental genes transcriptionally repressed in adult brains, indicating that it demarcates a deeply conserved regulatory program. We also found that the writer of non-CpG methylation, DNMT3A, and the reader, MeCP2, originated at the onset of vertebrates as a result of the ancestral vertebrate whole-genome duplication. Together, we demonstrate how this novel layer of epigenetic information assembled at the root of vertebrates and gained new regulatory roles independent of the ancestral form of the canonical CpG methylation. This suggests that the emergence of non-CpG methylation may have fostered the evolution of sophisticated cognitive abilities found in the vertebrate lineage.

Rapid changes in DNA methylation associated with the initiation of reproduction in a small songbird

Species with a circannual life cycle need to match the timing of their life history events to the environment to maximize fitness. However, our understanding of how circannual traits such as timing of reproduction are regulated on a molecular level remains limited. Recent studies have implicated that epigenetic mechanisms can be an important part in the processes that regulate circannual traits. Here, we explore the role of DNA methylation in mediating reproductive timing in a seasonally breeding bird species, the great tit (Parus major), using genome‐wide DNA methylation data from individual females that were blood sampled repeatedly throughout the breeding season. We demonstrate rapid and directional changes in DNA methylation within the promoter region of several genes, including a key transcription factor (NR5A1) known from earlier studies to be involved in the initiation of timing of reproduction. Interestingly, the observed changes in DNA methylation at NR5A1 identified here are in line with earlier gene expression studies of reproduction in chicken, indicating that the observed shifts in DNA methylation at this gene can have a regulatory role. Our findings provide an important step towards elucidating the genomic mechanism that mediates seasonal timing of a key life history traits and provide support for the idea that epigenetic mechanisms may play an important role in circannual traits.

see also the Perspective by Melanie J. Heckwolf and Britta S. Meyer

Temporal changes in DNA methylation and RNA expression in a small song bird: within- and between-tissue comparisons

Background

DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n = 6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled.

Results

We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10 kb up- and downstream regions adjacent to the gene body.

Conclusion

Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case. However, associations between temporal changes in DNA methylation with changes in gene expression are mostly tissue- and genomic location-dependent. The observation that temporal changes in DNA methylation within RBCs can relate to changes in gene expression in less accessible tissues is important for a better understanding of how environmental conditions shape traits that temporally change in expression in wild populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07329-9.

Urbanization is associated with modifications in DNA methylation in a small passerine bird

Urbanization represents a fierce driver of phenotypic change, yet the molecular mechanisms underlying observed phenotypic patterns are poorly understood. Epigenetic changes are expected to facilitate more rapid adaption to changing or novel environments, such as our towns and cities, compared with slow changes in gene sequence. A comparison of liver and blood tissue from great tits Parus major originating from an urban and a forest site demonstrated that urbanization is associated with variation in genome-wide patterns of DNA methylation. Combining reduced representation bisulphite sequencing with transcriptome data, we revealed habitat differences in DNA methylation patterns that suggest a regulated and coordinated response to the urban environment. In the liver, genomic sites that were differentially methylated between urban- and forest-dwelling birds were over-represented in regulatory regions of the genome and more likely to occur in expressed genes. DNA methylation levels were also inversely correlated with gene expression at transcription start sites. Furthermore, differentially methylated CpG sites, in liver, were over-represented in pathways involved in (i) steroid biosynthesis, (ii) superoxide metabolism, (iii) secondary alcohol metabolism, (iv) chylomicron remodelling, (v) cholesterol transport, (vi) reactive oxygen species (ROS) metabolic process and (vii) epithelial cell proliferation. This corresponds with earlier studies identifying diet and exposure to ROS as two of the main drivers of divergence between organisms in urban and nonurban environments. Conversely, in blood, sites that were differentially methylated between urban- and forest-dwelling birds were under-represented in regulatory regions, more likely to occur in nonexpressed genes and not over-represented in specific biological pathways. It remains to be determined whether diverging patterns of DNA methylation represent adaptive evolutionary responses and whether the conclusions can be more widely attributed to urbanization.

On the Use of Blood Samples for Measuring DNA Methylation in Ecological Epigenetic Studies

There is increasing interest in understanding the potential for epigenetic factors to contribute to phenotypic diversity in evolutionary biology. One well studied epigenetic mechanism is DNA methylation, the addition of a methyl group to cytosines, which have the potential to alter gene expression depending on the genomic region in which it takes place. Obtaining information about DNA methylation at genome-wide scale has become straightforward with the use of bisulfite treatment in combination with reduced representation or whole-genome sequencing. While it is well recognized that methylation is tissue specific, a frequent limitation for many studies is that sampling-specific tissues may require sacrificing individuals, something which is generally undesirable and sometimes impossible. Instead, information about DNA methylation patterns in the blood is frequently used as a proxy tissue. This can obviously be problematic if methylation patterns in the blood do not reflect that in the relevant tissue. Understanding how, or if, DNA methylation in blood reflect DNA methylation patterns in other tissues is therefore of utmost importance if we are to make inferences about how observed differences in methylation or temporal changes in methylation can contribute to phenotypic variation. The aim of this review is to examine what we know about the potential for using blood samples in ecological epigenetic studies. I briefly outline some methods by which we can measure DNA methylation before I examine studies that have compared DNA methylation patterns across different tissues and, finally, examine how useful blood samples may be for ecological studies of DNA methylation. Ecological epigenetic studies are in their infancy, but it is paramount for the field to move forward to have detailed information about tissue and time dependence relationships in methylation to gain insights into if blood DNA methylation patterns can be a reliable bioindicator for changes in methylation that generate phenotypic variation in ecologically important traits.

Epigenetics of Animal Personality: DNA Methylation Cannot Explain the Heritability of Exploratory Behavior in a Songbird

The search for the hereditary mechanisms underlying quantitative traits traditionally focused on the identification of underlying genomic polymorphisms such as single-nucleotide polymorphisms. It has now become clear that epigenetic mechanisms, such as DNA methylation, can consistently alter gene expression over multiple generations. It is unclear, however, if and how DNA methylation can stably be transferred from one generation to the next and can thereby be a component of the heritable variation of a trait. In this study, we explore whether DNA methylation responds to phenotypic selection using whole-genome and genome-wide bisulfite approaches. We assessed differential erythrocyte DNA methylation patterns between extreme personality types in the Great Tit (Parus major). For this, we used individuals from a four-generation artificial bi-directional selection experiment and siblings from eight F2 inter-cross families. We find no differentially methylated sites when comparing the selected personality lines, providing no evidence for the so-called epialleles associated with exploratory behavior. Using a pair-wise sibling design in the F2 intercrosses, we show that the genome-wide DNA methylation profiles of individuals are mainly explained by family structure, indicating that the majority of variation in DNA methylation in CpG sites between individuals can be explained by genetic differences. Although we found some candidates explaining behavioral differences between F2 siblings, we could not confirm this with a whole-genome approach, thereby confirming the absence of epialleles in these F2 intercrosses. We conclude that while epigenetic variation may underlie phenotypic variation in behavioral traits, we were not able to find evidence that DNA methylation can explain heritable variation in personality traits in Great Tits.

Genome-Wide Assessment of DNA Methylation in Chicken Cardiac Tissue Exposed to Different Incubation Temperatures and CO2 Levels

Temperature and CO₂ concentration during incubation have profound effects on broiler chick development, and numerous studies have identified significant effects on hatch heart weight (HW) as a result of differences in these parameters. Early life environment has also been shown to affect broiler performance later in life; it has thus been suggested that epigenetic mechanisms may mediate long-term physiological changes induced by environmental stimuli. DNA methylation is an epigenetic modification that can confer heritable changes in gene expression. Using reduced-representation bisulfite sequencing (RRBS), we assessed DNA methylation patterns in cardiac tissue of 84 broiler hatchlings incubated at two egg shell temperatures (EST; 37.8°C and 38.9°C) and three CO₂ concentrations (0.1%, 0.4%, and 0.8%) from day 8 of incubation onward. We assessed differential methylation between EST treatments and identified 2,175 differentially methylated (DM) CpGs (1,121 hypermethylated, 1,054 hypomethylated at 38.9° vs. 37.8°) in 269 gene promoters and 949 intragenic regions. DM genes (DMGs) were associated with heart developmental processes, including cardiomyocyte proliferation and differentiation. We identified enriched binding motifs among DM loci, including those for transcription factors associated with cell proliferation and heart development among hypomethylated CpGs that suggest increased binding ability at higher EST. We identified 9,823 DM CpGs between at least two CO₂ treatments, with the greatest difference observed between 0.8 and 0.1% CO₂ that disproportionately impacted genes involved in cardiac muscle development and response to low oxygen levels. Using HW measurements from the same chicks, we performed an epigenome-wide association study (EWAS) for HW, and identified 23 significantly associated CpGs, nine of which were also DM between ESTs. We found corresponding differences in transcript abundance between ESTs in three DMGs (ABLIM2, PITX2, and THRSP). Hypomethylation of an exonic CpG in PITX2 at 38.9°C was associated with increased expression, and suggests increased cell proliferation in broiler hatchlings incubated at higher temperatures. Overall, these results identified numerous epigenetic associations between chick incubation factors and heart development that may manifest in long-term differences in animal performance.

DNA methylation reprogramming, TE derepression, and postzygotic isolation of nascent animal species

The genomic shock hypothesis stipulates that the stress associated with divergent genome admixture can cause transposable element (TE) derepression, which could act as a postzygotic isolation mechanism. TEs affect gene structure, expression patterns, and chromosome organization and may have deleterious consequences when released. For these reasons, they are silenced by heterochromatin formation, which includes DNA methylation. Here, we show that a significant proportion of TEs are differentially methylated between the "dwarf" (limnetic) and the "normal" (benthic) whitefish, two nascent species that diverged some 15,000 generations ago within the Coregonus clupeaformis species complex. Moreover, TEs are overrepresented among loci that were demethylated in hybrids, indicative of their transcriptional derepression. These results are consistent with earlier studies in this system that revealed TE transcriptional derepression causes abnormal embryonic development and death of hybrids. Hence, this supports a role of DNA methylation reprogramming and TE derepression in postzygotic isolation of nascent animal species.

Fine-tuning of seasonal timing of breeding is regulated downstream in the underlying neuro-endocrine system in a small songbird

The timing of breeding is under selection in wild populations as a result of climate change, and understanding the underlying physiological processes mediating this timing provides insight into the potential rate of adaptation. Current knowledge on this variation in physiology is, however, mostly limited to males. We assessed whether individual differences in the timing of breeding in females are reflected in differences in candidate gene expression and, if so, whether these differences occur in the upstream (hypothalamus) or downstream (ovary and liver) parts of the neuroendocrine system. We used 72 female great tits from two generations of lines artificially selected for early and late egg laying, which were housed in climate-controlled aviaries and went through two breeding cycles within 1 year. In the first breeding season we obtained individual egg-laying dates, while in the second breeding season, using the same individuals, we sampled several tissues at three time points based on the timing of the first breeding attempt. For each tissue, mRNA expression levels were measured using qPCR for a set of candidate genes associated with the timing of reproduction and subsequently analysed for differences between generations, time points and individual timing of breeding. We found differences in gene expression between generations in all tissues, with the most pronounced differences in the hypothalamus. Differences between time points, and early- and late-laying females, were found exclusively in the ovary and liver. Altogether, we show that fine-tuning of the seasonal timing of breeding, and thereby the opportunity for adaptation in the neuroendocrine system, is regulated mostly downstream in the neuro-endocrine system.

Exploration of tissue-specific gene expression patterns underlying timing of breeding in contrasting temperature environments in a song bird

BACKGROUND

Seasonal timing of breeding is a life history trait with major fitness consequences but the genetic basis of the physiological mechanism underlying it, and how gene expression is affected by date and temperature, is not well known. In order to study this, we measured patterns of gene expression over different time points in three different tissues of the hypothalamic-pituitary-gonadal-liver axis, and investigated specifically how temperature affects this axis during breeding. We studied female great tits (Parus major) from lines artificially selected for early and late timing of breeding that were housed in two contrasting temperature environments in climate-controlled aviaries. We collected hypothalamus, liver and ovary samples at three different time points (before and after onset of egg-laying). For each tissue, we sequenced whole transcriptomes of 12 pools (n = 3 females) to analyse gene expression.

RESULTS

Birds from the selection lines differed in expression especially for one gene with clear reproductive functions, zona pellucida glycoprotein 4 (ZP4), which has also been shown to be under selection in these lines. Genes were differentially expressed at different time points in all tissues and most of the differentially expressed genes between the two temperature treatments were found in the liver. We identified a set of hub genes from all the tissues which showed high association to hormonal functions, suggesting that they have a core function in timing of breeding. We also found ample differentially expressed genes with largely unknown functions in birds.

CONCLUSIONS

We found differentially expressed genes associated with selection line and temperature treatment. Interestingly, the latter mainly in the liver suggesting that temperature effects on egg-laying date may happen down-stream in the physiological pathway. These findings, as well as our datasets, will further the knowledge of the mechanisms of tissue-specific avian seasonality in the future.

Temporally replicated DNA methylation patterns in great tit using reduced representation bisulfite sequencing

Seasonal timing of reproduction is an important fitness trait in many plants and animals but the underlying molecular mechanism for this trait is poorly known. DNA methylation is known to affect timing of reproduction in various organisms and is therefore a potential mechanism also in birds. Here we describe genome wide data aiming to detect temporal changes in methylation in relation to timing of breeding using artificial selection lines of great tits (Parus major) exposed to contrasting temperature treatments. Methylation levels of DNA extracted from erythrocytes were examined using reduced representation bisulfite sequencing (RRBS). In total, we obtained sequencing data from 63 libraries over four different time points from 16 birds with on average 20 million quality filtered reads per library. These data describe individual level temporal variation in DNA methylation throughout the breeding season under experimental temperature regimes and provides a resource for future studies investigating the role of temporal changes in DNA methylation in timing of reproduction.

Association of PON1 gene promoter DNA methylation with the risk of Clopidogrel resistance in patients with coronary artery disease

BACKGROUND AND AIMS

The failure of therapeutic response to clopidogrel in platelet inhibition, which is called clopidogrel resistance (CR), is more likely to cause cardiovascular events. We aimed to study the contribution of promoter DNA methylation of paraoxonase 1 (PON1) to the risk of clopidogrel poor response.

METHODS

Through VerifyNow P2Y12 assay, patient' platelet functions were measured. Among 57 non-CR and 49 CR patients, the levels of DNA methylation in four CpG dinucleotides on the PON1 promoter were tested using bisulfite pyrosequencing technology. Besides, the relative expression of PON1 mRNA was analyzed by quantitative real-time PCR. Logistic regression was applied to investigate the interaction of PON1 methylation and clinical factors in CR.

RESULTS

In the subgroup with dyslipidemia, we discovered that higher CpG4 levels of the PON1 promoter indicated a poorer clopidogrel response (cases versus controls (%): 51.500 ± 14.742 vs 43.308 ± 10.891, P = 0.036), and the PON1 mRNA expression was reduced in CR patients. Additionally, the logistic regression indicated that higher level of albumin and the index of ALT were related to a lower risk of CR, and the index of AST as well as the quantity of stent may be positively associated with CR.

CONCLUSIONS

The DNA methylation of CpG4 in the PON1 promoter would lead to a low expression of PON1 mRNA, which might induce clopidogrel resistance in the patients with dyslipidemia, and the number of stents might be a risk for CR.

Seasonal Variation in Genome-Wide DNA Methylation Patterns and the Onset of Seasonal Timing of Reproduction in Great Tits

In seasonal environments, timing of reproduction is a trait with important fitness consequences, but we know little about the molecular mechanisms that underlie the variation in this trait. Recently, several studies put forward DNA methylation as a mechanism regulating seasonal timing of reproduction in both plants and animals. To understand the involvement of DNA methylation in seasonal timing of reproduction, it is necessary to examine within-individual temporal changes in DNA methylation, but such studies are very rare. Here, we use a temporal sampling approach to examine changes in DNA methylation throughout the breeding season in female great tits (Parus major) that were artificially selected for early timing of breeding. These females were housed in climate-controlled aviaries and subjected to two contrasting temperature treatments. Reduced representation bisulfite sequencing on red blood cell derived DNA showed genome-wide temporal changes in more than 40,000 out of the 522,643 CpG sites examined. Although most of these changes were relatively small (mean within-individual change of 6%), the sites that showed a temporal and treatment-specific response in DNA methylation are candidate sites of interest for future studies trying to understand the link between DNA methylation patterns and timing of reproduction.

Exploring the unmapped DNA and RNA reads in a songbird genome

BACKGROUND

A widely used approach in next-generation sequencing projects is the alignment of reads to a reference genome. Despite methodological and hardware improvements which have enhanced the efficiency and accuracy of alignments, a significant percentage of reads frequently remain unmapped. Usually, unmapped reads are discarded from the analysis process, but significant biological information and insights can be uncovered from these data. We explored the unmapped DNA (normal and bisulfite treated) and RNA sequence reads of the great tit (Parus major) reference genome individual. From the unmapped reads we generated de novo assemblies, after which the generated sequence contigs were aligned to the NCBI non-redundant nucleotide database using BLAST, identifying the closest known matching sequence.

RESULTS

Many of the aligned contigs showed sequence similarity to different bird species and genes that were absent in the great tit reference assembly. Furthermore, there were also contigs that represented known P. major pathogenic species. Most interesting were several species of blood parasites such as Plasmodium and Trypanosoma.

CONCLUSIONS

Our analyses revealed that meaningful biological information can be found when further exploring unmapped reads. For instance, it is possible to discover sequences that are either absent or misassembled in the reference genome, and sequences that indicate infection or sample contamination. In this study we also propose strategies to aid the capture and interpretation of this information from unmapped reads.

Age- and quality-dependent DNA methylation correlate with melanin-based coloration in a wild bird

Secondary sexual trait expression can be influenced by fixed individual factors (such as genetic quality) as well as by dynamic factors (such as age and environmentally induced gene expression) that may be associated with variation in condition or quality. In particular, melanin-based traits are known to relate to condition and there is a well-characterized genetic pathway underpinning their expression. However, the mechanisms linking variable trait expression to genetic quality remain unclear. One plausible mechanism is that genetic quality could influence trait expression via differential methylation and differential gene expression. We therefore conducted a pilot study examining DNA methylation at a candidate gene (agouti-related neuropeptide: AgRP) in the black grouse Lyrurus tetrix. We specifically tested whether CpG methylation covaries with age and multilocus heterozygosity (a proxy of genetic quality) and from there whether the expression of a melanin-based ornament (ultraviolet-blue chroma) correlates with DNA methylation. Consistent with expectations, we found clear evidence for age- and heterozygosity-specific patterns of DNA methylation, with two CpG sites showing the greatest DNA methylation in highly heterozygous males at their peak age of reproduction. Furthermore, DNA methylation at three CpG sites was significantly positively correlated with ultraviolet-blue chroma. Ours is the first study to our knowledge to document age- and quality-dependent variation in DNA methylation and to show that dynamic sexual trait expression across the lifespan of an organism is associated with patterns of DNA methylation. Although we cannot demonstrate causality, our work provides empirical support for a mechanism that could potentially link key individual factors to variation in sexual trait expression in a wild vertebrate.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

CNVs are associated with genomic architecture in a songbird

Background

Understanding variation in genome structure is essential to understand phenotypic differences within populations and the evolutionary history of species. A promising form of this structural variation is copy number variation (CNV). CNVs can be generated by different recombination mechanisms, such as non-allelic homologous recombination, that rely on specific characteristics of the genome architecture. These structural variants can therefore be more abundant at particular genes ultimately leading to variation in phenotypes under selection. Detailed characterization of CNVs therefore can reveal evolutionary footprints of selection and provide insight in their contribution to phenotypic variation in wild populations.

Results

Here we use genotypic data from a long-term population of great tits (Parus major), a widely studied passerine bird in ecology and evolution, to detect CNVs and identify genomic features prevailing within these regions. We used allele intensities and frequencies from high-density SNP array data from 2,175 birds. We detected 41,029 CNVs concatenated into 8,008 distinct CNV regions (CNVRs). We successfully validated 93.75% of the CNVs tested by qPCR, which were sampled at different frequencies and sizes. A mother-daughter family structure allowed for the evaluation of the inheritance of a number of these CNVs. Thereby, only CNVs with 40 probes or more display segregation in accordance with Mendelian inheritance, suggesting a high rate of false negative calls for smaller CNVs. As CNVRs are a coarse-grained map of CNV loci, we also inferred the frequency of coincident CNV start and end breakpoints. We observed frequency-dependent enrichment of these breakpoints at homologous regions, CpG sites and AT-rich intervals. A gene ontology enrichment analyses showed that CNVs are enriched in genes underpinning neural, cardiac and ion transport pathways.

Conclusion

Great tit CNVs are present in almost half of the genes and prominent at repetitive-homologous and regulatory regions. Although overlapping genes under selection, the high number of false negatives make neutrality or association tests on CNVs detected here difficult. Therefore, CNVs should be further addressed in the light of their false negative rate and architecture to improve the comprehension of their association with phenotypes and evolutionary history.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4577-1) contains supplementary material, which is available to authorized users.

Insights into Epigenome Evolution from Animal and Plant Methylomes

Evolutionary studies of DNA methylation offer insights into the mechanisms governing the variation of genomic DNA methylation across different species. Comparisons of gross levels of DNA methylation between distantly related species indicate that the size of the genome and the level of genomic DNA methylation are positively correlated. In plant genomes, this can be reliably explained by the genomic contents of repetitive sequences. In animal genomes, the role of repetitive sequences on genomic DNA methylation is less clear. On a shorter timescale, population-level comparisons demonstrate that genetic variation can explain the observed variability of DNA methylation to some degree. The amount of DNA methylation variation that has been attributed to genetic variation in the human population studies so far is substantially lower than that from Arabidopsis population studies, but this disparity might reflect the differences in the computational and experimental techniques used. The effect of genetic variation on DNA methylation has been directly examined in mammalian systems, revealing several causative factors that govern DNA methylation. On the other hand, studies from Arabidopsis have furthered our understanding of spontaneous mutations of DNA methylation, termed “epimutations.” Arabidopsis has an extremely high rate of spontaneous epimutations, which may play a major role in shaping the global DNA methylation landscape in this genome. Key missing information includes the frequencies of spontaneous epimutations in other lineages, in particular animal genomes, and how population-level variation of DNA methylation leads to species-level differences.

Genomic tools for behavioural ecologists to understand repeatable individual differences in behaviour

Behaviour is a key interface between an animal's genome and its environment. Repeatable individual differences in behaviour have been extensively documented in animals, but the molecular underpinnings of behavioural variation among individuals within natural populations remain largely unknown. Here, we offer a critical review of when molecular techniques may yield new insights, and we provide specific guidance on how and whether the latest tools available are appropriate given different resources, system and organismal constraints, and experimental designs. Integrating molecular genetic techniques with other strategies to study the proximal causes of behaviour provides opportunities to expand rapidly into new avenues of exploration. Such endeavours will enable us to better understand how repeatable individual differences in behaviour have evolved, how they are expressed and how they can be maintained within natural populations of animals.

Genome-wide DNA methylation analysis of the porcine hypothalamus-pituitary-ovary axis

Previous studies have suggested that DNA methylation in both CpG and CpH (where H = C, T or A) contexts plays a critical role in biological functions of different tissues. However, the genome-wide DNA methylation patterns of porcine hypothalamus-pituitary-ovary (HPO) tissues remain virtually unexplored. In this study, methylomes of HPO tissues were profiled to investigate their differences and similarities. We found that HPO methylomes displayed tissue-specific methylation patterns in both CpG and CpH contexts. At gene locations, the methylation and density of CpGs was negatively linked at transcription start sites but positively linked at transcription end sites. The densities of CpGs and CpHs at CpG island (CGI) locations were negatively correlated with their methylation. Moreover, the methylation interactions between CGIs and genes showed similar pattern in the CpG context but tissue-specific pattern in the CpH context. CpGs located in CGIs, upstream regions and exons were protected from methylation dynamics, whereas CGI shores, CGI shelves and intergenic regions were more likely to be targets of methylation changes. The methylation dynamics enriching in a tissue-specific manner appeared to maintain and establish the biological functions of HPO tissues. Our analyses provided valuable insights into the tissue-specific methylomes of porcine HPO tissues.

The elephant shark methylome reveals conservation of epigenetic regulation across jawed vertebrates

BACKGROUND

Methylation of CG dinucleotides constitutes a critical system of epigenetic memory in bony vertebrates, where it modulates gene expression and suppresses transposon activity. The genomes of studied vertebrates are pervasively hypermethylated, with the exception of regulatory elements such as transcription start sites (TSSs), where the presence of methylation is associated with gene silencing. This system is not found in the sparsely methylated genomes of invertebrates, and establishing how it arose during early vertebrate evolution is impeded by a paucity of epigenetic data from basal vertebrates.

METHODS

 We perform whole-genome bisulfite sequencing to generate the first genome-wide methylation profiles of a cartilaginous fish, the elephant shark Callorhinchus milii. Employing these to determine the elephant shark methylome structure and its relationship with expression, we compare this with higher vertebrates and an invertebrate chordate using published methylation and transcriptome data.  Results: Like higher vertebrates, the majority of elephant shark CG sites are highly methylated, and methylation is abundant across the genome rather than patterned in the mosaic configuration of invertebrates. This global hypermethylation includes transposable elements and the bodies of genes at all expression levels. Significantly, we document an inverse relationship between TSS methylation and expression in the elephant shark, supporting the presence of the repressive regulatory architecture shared by higher vertebrates.

CONCLUSIONS

 Our demonstration that methylation patterns in a cartilaginous fish are characteristic of higher vertebrates imply the conservation of this epigenetic modification system across jawed vertebrates separated by 465 million years of evolution. In addition, these findings position the elephant shark as a valuable model to explore the evolutionary history and function of vertebrate methylation.