Developmental stress does not induce genome-wide DNA methylation changes in wild great tit (Parus major) nestlings.

The environment experienced during early life is a crucial factor in the life of many organisms. This early life environment has been shown to have profound effects on morphology, physiology and fitness. However, the molecular mechanisms that mediate these effects are largely unknown, even though they are essential for our understanding of the processes that induce phenotypic variation in natural populations. DNA methylation is an epigenetic mechanism that has been suggested to explain such environmentally induced phenotypic changes early in life. To investigate whether DNA methylation changes are associated with experimentally induced early developmental effects, we cross-fostered great tit (Parus major) nestlings and manipulated their brood sizes in a natural study population. We assessed experimental brood size effects on pre-fledging biometry and behaviour. We linked this to genome-wide DNA methylation levels of CpG sites in erythrocyte DNA, using 122 individuals and an improved epiGBS2 laboratory protocol. Brood enlargement caused developmental stress and negatively affected nestling condition, predominantly during the second half of the breeding season, when conditions are harsher. Brood enlargement, however, affected nestling DNA methylation in only one CpG site and only if the hatch date was taken into account. In conclusion, this study shows that nutritional stress in enlarged broods does not associate with direct effects on genome-wide DNA methylation. Future studies should assess whether genome-wide DNA methylation variation may arise later in life as a consequence of phenotypic changes during early development.

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.

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.

Evidence from pyrosequencing indicates that natural variation in animal personality is associated with DRD4 DNA methylation.

Personality traits are heritable and respond to natural selection, but are at the same time influenced by the ontogenetic environment. Epigenetic effects, such as DNA methylation, have been proposed as a key mechanism to control personality variation. However, to date little is known about the contribution of epigenetic effects to natural variation in behaviour. Here, we show that great tit (Parus major) lines artificially selected for divergent exploratory behaviour for four generations differ in their DNA methylation levels at the dopamine receptor D4 (DRD4) gene. This D4 receptor is statistically associated with personality traits in both humans and nonhuman animals, including the great tit. Previous work in this songbird failed to detect functional genetic polymorphisms within DRD4 that could account for the gene-trait association. However, our observation supports the idea that DRD4 is functionally involved in exploratory behaviour but that its effects are mediated by DNA methylation. While the exact mechanism underlying the transgenerational consistency of DRD4 methylation remains to be elucidated, this study shows that epigenetic mechanisms are involved in shaping natural variation in personality traits. We outline how this first finding provides a basis for investigating the epigenetic contribution to personality traits in natural systems and its subsequent role for understanding the ecology and evolution of behavioural consistency.

Genetic and epigenetic differentiation in response to genomic selection for avian lay date.

Anthropogenic climate change has led to globally increasing temperatures at an unprecedented pace and, to persist, wild species have to adapt to their changing world. We, however, often fail to derive reliable predictions of species' adaptive potential. Genomic selection represents a powerful tool to investigate the adaptive potential of a species, but constitutes a 'blind process' with regard to the underlying genomic architecture of the relevant phenotypes. Here, we used great tit (Parus major) females from a genomic selection experiment for avian lay date to zoom into this blind process. We aimed to identify the genetic variants that responded to genomic selection and epigenetic variants that accompanied this response and, this way, might reflect heritable genetic variation at the epigenetic level. We applied whole genome bisulfite sequencing to blood samples of individual great tit females from the third generation of bidirectional genomic selection lines for early and late lay date. Genomic selection resulted in differences at both the genetic and epigenetic level. Genetic variants that showed signatures of selection were located within genes mostly linked to brain development and functioning, including LOC107203824 (SOX3-like). SOX3 is a transcription factor that is required for normal hypothalamo-pituitary axis development and functioning, an essential part of the reproductive axis. As for epigenetic differentiation, the early selection line showed hypomethylation relative to the late selection line. Sites with differential DNA methylation were located in genes important for various biological processes, including gonadal functioning (e.g., MSTN and PIK3CB). Overall, genomic selection for avian lay date provided insights into where within the genome the heritable genetic variation for lay date, on which selection can operate, resides and indicates that some of this variation might be reflected by epigenetic variants.

Opposite differential allocation by males and females of the same species.

Differential allocation (DA)-the adjustment of an individual's parental investment in relation to its mate's attractiveness-is increasingly recognized as an important component of sexual selection. However, although DA is expected by both sexes of parents in species with biparental care, DA by males has rarely been investigated. We have previously demonstrated a decrease in the feeding rates of female blue tits Cyanistes caeruleus when their mate's UV coloration was experimentally reduced (i.e. positive DA). In this study, we used the same experimental protocol in the same population to investigate DA by male blue tits in relation to their female's UV coloration. Males mated to UV-reduced females had higher feeding rates than those mated to control females (i.e. negative DA). Thus, male and female blue tits display opposite DA for the same component of parental effort (chick provisioning), the first time that this has been reported for any species.

Genotypes selected for early and late avian lay date differ in their phenotype, but not fitness, in the wild.

Global warming has shifted phenological traits in many species, but whether species are able to track further increasing temperatures depends on the fitness consequences of additional shifts in phenological traits. To test this, we measured phenology and fitness of great tits (Parus major) with genotypes for extremely early and late egg lay dates, obtained from a genomic selection experiment. Females with early genotypes advanced lay dates relative to females with late genotypes, but not relative to nonselected females. Females with early and late genotypes did not differ in the number of fledglings produced, in line with the weak effect of lay date on the number of fledglings produced by nonselected females in the years of the experiment. Our study is the first application of genomic selection in the wild and led to an asymmetric phenotypic response that indicates the presence of constraints toward early, but not late, lay dates.

Integrated molecular and behavioural data reveal deep circadian disruption in response to artificial light at night in male Great tits (Parus major).

Globally increasing levels of artificial light at night (ALAN) are associated with shifting rhythms of behaviour in many wild species. However, it is unclear whether changes in behavioural timing are paralleled by consistent shifts in the molecular clock and its associated physiological pathways. Inconsistent shifts between behavioural and molecular rhythms, and between different tissues and physiological systems, disrupt the circadian system, which coordinates all major body functions. We therefore compared behavioural, transcriptional and metabolomic responses of captive great tits (Parus major) to three ALAN intensities or to dark nights, recording activity and sampling brain, liver, spleen and blood at mid-day and midnight. ALAN advanced wake-up time, and this shift was paralleled by advanced expression of the clock gene BMAL1 in all tissues, suggesting close links between behaviour and clock gene expression across tissues. However, further analysis of gene expression and metabolites revealed that clock shifts were inconsistent across physiological systems. Untargeted metabolomic profiling showed that only 9.7% of the 755 analysed metabolites followed the behavioural shift. This high level of desynchronization indicates that ALAN disrupted the circadian system on a deep, easily overlooked level. Thus, circadian disruption could be a key mediator of health impacts of ALAN on wild animals.

epiGBS2: Improvements and evaluation of highly multiplexed, epiGBS-based reduced representation bisulfite sequencing.

Several reduced-representation bisulfite sequencing methods have been developed in recent years to determine cytosine methylation de novo in nonmodel species. Here, we present epiGBS2, a laboratory protocol based on epiGBS with a revised and user-friendly bioinformatics pipeline for a wide range of species with or without a reference genome. epiGBS2 is cost- and time-efficient and the computational workflow is designed in a user-friendly and reproducible manner. The library protocol allows a flexible choice of restriction enzymes and a double digest. The bioinformatics pipeline was integrated in the Snakemake workflow management system, which makes the pipeline easy to execute and modular, and parameter settings for important computational steps flexible. We implemented bismark for alignment and methylation analysis and we preprocessed alignment files by double masking to enable single nucleotide polymorphism calling with Freebayes (epiFreebayes). The performance of several critical steps in epiGBS2 was evaluated against baseline data sets from Arabidopsis thaliana and great tit (Parus major), which confirmed its overall good performance. We provide a detailed description of the laboratory protocol and an extensive manual of the bioinformatics pipeline, which is publicly accessible on github (https://github.com/nioo-knaw/epiGBS2) and zenodo (https://doi.org/10.5281/zenodo.4764652).

Effects of experimental light at night on extra-pair paternity in a songbird.

Light pollution is increasing worldwide and significantly affects animal behavior. In birds, these effects include advancement of morning activity and onset of dawn song, which may affect extra-pair paternity. Advanced dawn song of males may stimulate females to engage in extra-pair copulations, and the earlier activity onset may affect the males' mate guarding behavior. Earlier work showed an effect of light at night on extra-pair behavior, but this was in an area with other anthropogenic disturbances. Here, we present a two-year experimental study on effects of light at night on extra-pair paternity of great tits (Parus major). Previously dark natural areas were illuminated with white, red, and green LED lamps and compared to a dark control. In 2014, the proportion of extra-pair young in broods increased with distance to the red and white lamps (i.e., at lower light intensities), but decreased with distance to the poles in the dark control. In 2013, we found no effects on the proportion of extra-pair young. The total number of offspring sired by a male was unaffected by artificial light at night in both years, suggesting that potential changes in female fidelity in pairs breeding close to white and red light did not translate into fitness benefits for the males of these pairs. Artificial light at night might disrupt the natural patterns of extra-pair paternity, possibly negates potential benefits of extra-pair copulations and thus could alter sexual selection processes in wild birds.

Female blue tits adjust parental effort to manipulated male UV attractiveness.

The differential allocation hypothesis predicts that parents should adjust their current investment in relation to perceived mate attractiveness if this affects offspring fitness. It should be selectively advantageous to risk more of their future reproductive success by investing heavily in current offspring of high reproductive value but to decrease investment if offspring value is low. If the benefits of mate attractiveness are limited to a particular offspring sex we would instead expect relative investment in male versus female offspring to vary with mate attractiveness, referred to as 'differential sex allocation'. We present strong evidence for differential allocation of parental feeding effort in the wild and show an immediate effect on a component of offspring fitness. By experimentally reducing male UV crown coloration, a trait known to indicate attractiveness and viability in wild-breeding blue tits (Parus caeruleus), we show that females, but not males, reduce parental feeding rates and that this reduces the skeletal growth of offspring. However, differential sex allocation does not occur. We conclude that blue tit females use male UV coloration as an indicator of expected offspring fitness and adjust their investment accordingly.

Transfer of maternal antibodies against avian influenza virus in mallards (Anas platyrhynchos).

Maternal antibodies protect chicks from infection with pathogens early in life and may impact pathogen dynamics due to the alteration of the proportion of susceptible individuals in a population. We investigated the transfer of maternal antibodies against avian influenza virus (AIV) in a key AIV host species, the mallard (Anas platyrhynchos). Combining observations in both the field and in mallards kept in captivity, we connected maternal AIV antibody concentrations in eggs to (i) female body condition, (ii) female AIV antibody concentration, (iii) egg laying order, (iv) egg size and (v) embryo sex. We applied maternity analysis to the eggs collected in the field to account for intraspecific nest parasitism, which is reportedly high in Anseriformes, detecting parasitic eggs in one out of eight clutches. AIV antibody prevalence in free-living and captive females was respectively 48% and 56%, with 43% and 24% of the eggs receiving these antibodies maternally. In both field and captive study, maternal AIV antibody concentrations in egg yolk correlated positively with circulating AIV antibody concentrations in females. In the captive study, yolk AIV antibody concentrations correlated positively with egg laying order. Female body mass and egg size from the field and captive study, and embryos sex from the field study were not associated with maternal AIV antibody concentrations in eggs. Our study indicates that maternal AIV antibody transfer may potentially play an important role in shaping AIV infection dynamics in mallards.