Evolutionary signals of selection on cognition from the great tit genome and methylome

The great tit HapMap project: A continental-scale analysis of genomic variation in a songbird

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude - almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics.

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.

Early developmental carry-over effects on exploratory behaviour and DNA methylation in wild great tits (Parus major)

Adverse, postnatal conditions experienced during development are known to induce lingering effects on morphology, behaviour, reproduction and survival. Despite the importance of early developmental stress for shaping the adult phenotype, it is largely unknown which molecular mechanisms allow for the induction and maintenance of such phenotypic effects once the early environmental conditions are released. Here we aimed to investigate whether lasting early developmental phenotypic changes are associated with post-developmental DNA methylation changes. We used a cross-foster and brood size experiment in great tit (Parus major) nestlings, which induced post-fledging effects on biometric measures and exploratory behaviour, a validated personality trait. We investigated whether these post-fledging effects are associated with DNA methylation levels of CpG sites in erythrocyte DNA. Individuals raised in enlarged broods caught up on their developmental delay after reaching independence and became more explorative as days since fledging passed, while the exploratory scores of individuals that were raised in reduced broods remained stable. Although we previously found that brood enlargement hardly affected the pre-fledging methylation levels, we found 420 CpG sites that were differentially methylated between fledged individuals that were raised in small versus large sized broods. A considerable number of the affected CpG sites were located in or near genes involved in metabolism, growth, behaviour and cognition. Since the biological functions of these genes line up with the observed post-fledging phenotypic effects of brood size, our results suggest that DNA methylation provides organisms the opportunity to modulate their condition once the environmental conditions allow it. In conclusion, this study shows that nutritional stress imposed by enlarged brood size during early development associates with variation in DNA methylation later in life. We propose that treatment-associated DNA methylation differences may arise in relation to pre- or post-fledging phenotypic changes, rather than that they are directly induced by the environment during early development.

The regulation of methylation on the Z chromosome and the identification of multiple novel Male Hyper-Methylated regions in the chicken

DNA methylation is a key regulator of eukaryote genomes, and is of particular relevance in the regulation of gene expression on the sex chromosomes, with a key role in dosage compensation in mammalian XY systems. In the case of birds, dosage compensation is largely absent, with it being restricted to two small Male Hyper-Methylated (MHM) regions on the Z chromosome. To investigate how variation in DNA methylation is regulated on the Z chromosome we utilised a wild x domestic advanced intercross in the chicken, with both hypothalamic methylomes and transcriptomes assayed in 124 individuals. The relatively large numbers of individuals allowed us to identify additional genomic MHM regions on the Z chromosome that were significantly differentially methylated between the sexes. These regions appear to down-regulate local gene expression in males, but not remove it entirely (unlike the lncRNAs identified in the initial MHM regions). These MHM regions were further tested and the most balanced genes appear to show decreased expression in males, whilst methylation appeared to be far more correlated with gene expression in the less balanced, as compared to the most balanced genes. In addition, quantitative trait loci (QTL) that regulate variation in methylation on the Z chromosome, and those loci that regulate methylation on the autosomes that derive from the Z chromosome were mapped. Trans-effect hotspots were also identified that were based on the autosomes but affected the Z, and also one that was based on the Z chromosome but that affected both autosomal and sex chromosome DNA methylation regulation. We show that both cis and trans loci that originate from the Z chromosome never exhibit an interaction with sex, whereas trans loci originating from the autosomes but affecting the Z chromosome always display such an interaction. Our results highlight how additional MHM regions are actually present on the Z chromosome, and they appear to have smaller-scale effects on gene expression in males. Quantitative variation in methylation is also regulated both from the autosomes to the Z chromosome, and from the Z chromosome to the autosomes.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

Genetic and context-specific effects on individual inhibitory control performance in the guppy (Poecilia reticulata)

Among-individual variation in cognitive traits, widely assumed to have evolved under adaptive processes, is increasingly being demonstrated across animal taxa. As variation among individuals is required for natural selection, characterizing individual differences and their heritability is important to understand how cognitive traits evolve. Here, we use a quantitative genetic study of wild-type guppies repeatedly exposed to a 'detour task' to test for genetic variance in the cognitive trait of inhibitory control. We also test for genotype-by-environment interactions (GxE) by testing related fish under alternative experimental treatments (transparent vs. semi-transparent barrier in the detour-task). We find among-individual variation in detour task performance, consistent with differences in inhibitory control. However, analysis of GxE reveals that heritable factors only contribute to performance variation in one treatment. This suggests that the adaptive evolutionary potential of inhibitory control (and/or other latent variables contributing to task performance) may be highly sensitive to environmental conditions. The presence of GxE also implies that the plastic response of detour task performance to treatment environment is genetically variable. Our results are consistent with a scenario where variation in individual inhibitory control stems from complex interactions between heritable and plastic components.

Dietary nucleotides can prevent glucocorticoid-induced telomere attrition in a fast-growing wild vertebrate

Telomeres are chromosome protectors that shorten during eukaryotic cell replication and in stressful conditions. Developing individuals are susceptible to telomere erosion when their growth is fast and resources are limited. This is critical because the rate of telomere attrition in early life is linked to health and life span of adults. The metabolic telomere attrition hypothesis (MeTA) suggests that telomere dynamics can respond to biochemical signals conveying information about the organism's energetic state. Among these signals are glucocorticoids, hormones that promote catabolic processes, potentially impairing costly telomere maintenance, and nucleotides, which activate anabolic pathways through the cellular enzyme target of rapamycin (TOR), thus preventing telomere attrition. During the energetically demanding growth phase, the regulation of telomeres in response to two contrasting signals - one promoting telomere maintenance and the other attrition - provides an ideal experimental setting to test the MeTA. We studied nestlings of a rapidly developing free-living passerine, the great tit (Parus major), that either received glucocorticoids (Cort-chicks), nucleotides (Nuc-chicks) or a combination of both (NucCort-chicks), comparing these with controls (Cnt-chicks). As expected, Cort-chicks showed telomere attrition, while NucCort- and Nuc-chicks did not. NucCort-chicks was the only group showing increased expression of a proxy for TOR activation (the gene TELO2), of mitochondrial enzymes linked to ATP production (cytochrome oxidase and ATP-synthase) and a higher efficiency in aerobically producing ATP. NucCort-chicks had also a higher expression of telomere maintenance genes (shelterin protein TERF2 and telomerase TERT) and of enzymatic antioxidant genes (glutathione peroxidase and superoxide dismutase). The findings show that nucleotide availability is crucial for preventing telomere erosion during fast growth in stressful environments.

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.

From maternal glucocorticoid and thyroid hormones to epigenetic regulation of offspring gene expression: An experimental study in a wild bird species

Offspring phenotype at birth is determined by its genotype and the prenatal environment including exposure to maternal hormones. Variation in both maternal glucocorticoids and thyroid hormones can affect offspring phenotype, but the underlying molecular mechanisms, especially those contributing to long-lasting effects, remain unclear. Epigenetic changes (such as DNA methylation) have been postulated as mediators of long-lasting effects of early-life environment. In this study, we determined the effects of elevated prenatal glucocorticoid and thyroid hormones on handling stress response (breath rate) as well as DNA methylation and gene expression of glucocorticoid receptor (GR) and thyroid hormone receptor (THR) in great tits (Parus major). Eggs were injected before incubation onset with corticosterone (the main avian glucocorticoid) and/or thyroid hormones (thyroxine and triiodothyronine) to simulate variation in maternal hormone deposition. Breath rate during handling and gene expression of GR and THR were evaluated 14 days after hatching. Methylation status of GR and THR genes was analyzed from the longitudinal blood cells sampled 7 and 14 days after hatching, as well as the following autumn. Elevated prenatal corticosterone level significantly increased the breath rate during handling, indicating an enhanced metabolic stress response. Prenatal corticosterone manipulation had CpG-site-specific effects on DNA methylation at the GR putative promoter region, while it did not significantly affect GR gene expression. GR expression was negatively associated with earlier hatching date and chick size. THR methylation or expression did not exhibit any significant relationship with the hormonal treatments or the examined covariates, suggesting that TH signaling may be more robust due to its crucial role in development. This study provides some support to the hypothesis suggesting that maternal corticosterone may influence offspring metabolic stress response via epigenetic alterations, yet their possible adaptive role in optimizing offspring phenotype to the prevailing conditions, context-dependency, and the underlying molecular interplay needs further research.

Epigenetic modification of the hypothalamic-pituitary-adrenal (HPA) axis during development in the house sparrow (Passer domesticus)

Epigenetic modifications such as DNA methylation are important mechanisms for mediating developmental plasticity, where ontogenetic processes and their phenotypic outcomes are shaped by early environments. In particular, changes in DNA methylation of genes within the hypothalamic-pituitary-adrenal (HPA) axis can impact offspring growth and development. This relationship has been well documented in mammals but is less understood in other taxa. Here, we use target-enriched enzymatic methyl sequencing (TEEM-seq) to assess how DNA methylation in a suite of 25 genes changes over development, how these modifications relate to the early environment, and how they predict differential growth trajectories in the house sparrow (Passer domesticus). We found that DNA methylation changes dynamically over the postnatal developmental period: genes with initially low DNA methylation tended to decline in methylation over development, whereas genes with initially high DNA methylation tended to increase in methylation. However, sex-specific differentially methylated regions (DMRs) were maintained across the developmental period. We also found significant differences in post-hatching DNA methylation in relation to hatch date, with higher levels of DNA methylation in nestlings hatched earlier in the season. Although these differences were largely absent by the end of development, a number of DMRs in HPA-related genes (CRH, MC2R, NR3C1, NR3C2, POMC)-and to a lesser degree HPG-related genes (GNRHR2)-predicted nestling growth trajectories over development. These findings provide insight into the mechanisms by which the early environment shapes DNA methylation in the HPA axis, and how these changes subsequently influence growth and potentially mediate developmental plasticity.

Genomic Basis of Adaptive Divergence in Leg Length between Ground- and Tree-Dwelling Species within a Bird Family

Hind limbs of tetrapods vary greatly in length and the variability can be associated with locomotor adaptation. Although the phenotypic evolution has been well documented, the underlying genetic basis remains poorly understood. We address this issue by integrating comparative genomics and functional prediction with a study system consisting of ground-dwelling, long-legged and tree-dwelling, short-legged species within the avian family Paridae. Genome-wide divergence and phenotypic correlation analyses jointly identified five highly divergent genomic regions that are significantly related with the difference in leg length between these two groups. Gene annotation for these regions detected three genes involved in skeletal development, that is, PTPA, BRINP1, and MIGA2, with the first one being under the strongest selection. Furthermore, four single nucleotide polymorphisms (SNPs) in the coding region of PTPA can well distinguish the two groups with distinct leg length. Among the four SNPs, one is non-synonymous mutation, and according to the prediction for protein structure and function, it can affect the 3D structure of the encoded protein by altering the corresponding amino acid's position. The alleles of PTPA were found in all sequenced species of the orders Palaeognathae and Psittaciformes, which typically take a ground locomotion style. A whole-genome scanning across bird species uncovered that the four SNPs are more likely to be present in resident passerines with increased leg length/wing length ratios (a proxy of leg-dependent locomotion efficiency). Our findings provide insight into the molecular evolution of locomotion performance based on leg morphology in birds.

Chromatin accessibility, not 5mC methylation covaries with partial dosage compensation in crows

The evolution of genetic sex determination is often accompanied by degradation of the sex-limited chromosome. Male heterogametic systems have evolved convergent, epigenetic mechanisms restoring the resulting imbalance in gene dosage between diploid autosomes (AA) and the hemizygous sex chromosome (X). Female heterogametic systems (AAf Zf, AAm ZZm) tend to only show partial dosage compensation (0.5 < Zf:AAf < 1) and dosage balance (0.5 Collapse

Key Words Collapse

MESH Headings

• Animals
• Female
• Male
• Chromatin/genetics
• Crows/genetics
• Epigenesis, Genetic
• Methylation
• Dosage Compensation, Genetic
• Sex Chromosomes

Collapse

Grants

• HORIZON EUROPE European Research Council
• Swedish Research Council
• Knut och Alice Wallenbergs Stiftelse

Collapse

Affiliation(s)

Ana Catalán Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Justin Merondun Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Ulrich Knief Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany Evolutionary Biology & Ecology,Faculty of Biology, University of Freiburg, Freiburg, Germany
Jochen B. W. Wolf Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany

Collapse

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.

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.

Comparative analysis of genome-scale, base-resolution DNA methylation profiles across 580 animal species

Methylation of cytosines is a prototypic epigenetic modification of the DNA. It has been implicated in various regulatory mechanisms across the animal kingdom and particularly in vertebrates. We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale DNA methylation profiles of multiple organs. Bioinformatic analysis of this large dataset quantified the association of DNA methylation with the underlying genomic DNA sequence throughout vertebrate evolution. We observed a broadly conserved link with two major transitions-once in the first vertebrates and again with the emergence of reptiles. Cross-species comparisons focusing on individual organs supported a deeply conserved association of DNA methylation with tissue type, and cross-mapping analysis of DNA methylation at gene promoters revealed evolutionary changes for orthologous genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution.

The combination of genomic offset and niche modelling provides insights into climate change-driven vulnerability

Global warming is increasingly exacerbating biodiversity loss. Populations locally adapted to spatially heterogeneous environments may respond differentially to climate change, but this intraspecific variation has only recently been considered when modelling vulnerability under climate change. Here, we incorporate intraspecific variation in genomic offset and ecological niche modelling to estimate climate change-driven vulnerability in two bird species in the Sino-Himalayan Mountains. We found that the cold-tolerant populations show higher genomic offset but risk less challenge for niche suitability decline under future climate than the warm-tolerant populations. Based on a genome-niche index estimated by combining genomic offset and niche suitability change, we identified the populations with the least genome-niche interruption as potential donors for evolutionary rescue, i.e., the populations tolerant to climate change. We evaluated potential rescue routes via a landscape genetic analysis. Overall, we demonstrate that the integration of genomic offset, niche suitability modelling, and landscape connectivity can improve climate change-driven vulnerability assessments and facilitate effective conservation management.

Integrated Analyses of DNA Methylation and Gene Expression of Rainbow Trout Muscle under Variable Ploidy and Muscle Atrophy Conditions

Rainbow trout, Oncorhynchus mykiss, is an important cool, freshwater aquaculture species used as a model for biological research. However, its genome reference has not been annotated for epigenetic markers affecting various biological processes, including muscle growth/atrophy. Increased energetic demands during gonadogenesis/reproduction provoke muscle atrophy in rainbow trout. We described DNA methylation and its associated gene expression in atrophying muscle by comparing gravid, diploid females to sterile, triploid females. Methyl Mini-seq and RNA-Seq were simultaneously used to characterize genome-wide DNA methylation and its association with gene expression in rainbow trout muscle. Genome-wide enrichment in the number of CpGs, accompanied by depleted methylation levels, was noticed around the gene transcription start site (TSS). Hypermethylation of CpG sites within ±1 kb on both sides of TSS (promoter and gene body) was weakly/moderately associated with reduced gene expression. Conversely, hypermethylation of the CpG sites in downstream regions of the gene body +2 to +10 kb was weakly associated with increased gene expression. Unlike mammalian genomes, rainbow trout gene promotors are poor in CpG islands, at <1% compared to 60%. No signs of genome-wide, differentially methylated (DM) CpGs were observed due to the polyploidy effect; only 1206 CpGs (0.03%) were differentially methylated, and these were primarily associated with muscle atrophy. Twenty-eight genes exhibited differential gene expression consistent with methylation levels of 31 DM CpGs. These 31 DM CpGs represent potential epigenetic markers of muscle atrophy in rainbow trout. The DM CpG-harboring genes are involved in apoptosis, epigenetic regulation, autophagy, collagen metabolism, cell membrane functions, and Homeobox proteins. Our study also identified genes explaining higher water content and modulated glycolysis previously shown as characteristic biochemical signs of rainbow trout muscle atrophy associated with sexual maturation. This study characterized DNA methylation in the rainbow trout genome and its correlation with gene expression. This work also identified novel epigenetic markers associated with muscle atrophy in fish/lower vertebrates.

Epigenetics in Ecology, Evolution, and Conservation

Epigenetic variation is often characterized by modifications to DNA that do not alter the underlying nucleotide sequence, but can influence behavior, morphology, and physiological phenotypes by affecting gene expression and protein synthesis. In this review, we consider how the emerging field of ecological epigenetics (eco-epi) aims to use epigenetic variation to explain ecologically relevant phenotypic variation and predict evolutionary trajectories that are important in conservation. Here, we focus on how epigenetic data have contributed to our understanding of wild populations, including plants, animals, and fungi. First, we identified published eco-epi literature and found that there was limited taxonomic and ecosystem coverage and that, by necessity of available technology, these studies have most often focused on the summarized epigenome rather than locus- or nucleotide-level epigenome characteristics. We also found that while many studies focused on adaptation and heritability of the epigenome, the field has thematically expanded into topics such as disease ecology and epigenome-based ageing of individuals. In the second part of our synthesis, we discuss key insights that have emerged from the epigenetic field broadly and use these to preview the path toward integration of epigenetics into ecology. Specifically, we suggest moving focus to nucleotide-level differences in the epigenome rather than whole-epigenome data and that we incorporate several facets of epigenome characterization (e.g., methylation, chromatin structure). Finally, we also suggest that incorporation of behavior and stress data will be critical to the process of fully integrating eco-epi data into ecology, conservation, and evolutionary biology.

Sympatric speciation of the spiny mouse from Evolution Canyon in Israel substantiated genomically and methylomically

SignificanceWhether sympatric speciation (SS) is rare or common is still debated. Two populations of the spiny mouse, Acomys cahirinus, from Evolution Canyon I (EC I) in Israel have been depicted earlier as speciating sympatrically by molecular markers and transcriptome. Here, we investigated SS both genomically and methylomically, demonstrating that the opposite populations of spiny mice are sister taxa and split from the common ancestor around 20,000 years ago without an allopatric history. Mate choice, olfactory receptors, and speciation genes contributed to prezygotic/postzygotic reproductive isolation. The two populations showed different methylation patterns, facilitating adaptation to their local environment. They cope with abiotic and biotic stresses, due to high solar interslope radiation differences. We conclude that our new genomic and methylomic data substantiated SS.

Risk-taking coping style correlates with SERT SNP290 polymorphisms in free-living great tits

The coping style of an individual in relation to potentially dangerous situations has been suggested to be inherited in a polygenic fashion, being SERT one of the candidate genes. In this paper, we assessed in free-living great tits Parus major the association between SNP290 in the SERT promoter and three standard fear-related behaviors, namely the response of the birds to a black and white flag fixed to the top of the nest-box, distress calling rate of the birds in the hand once captured and the hissing call of incubating females when approached by a predator. We found a strong association between SNP290 polymorphism and the three risk-taking behaviors, with birds with genotype CT entering faster to the nest box with the flag and displaying more distress calls and less hissing calls. CT birds could therefore be described as more proactive than CC individuals. These results also suggest that hissing behavior should be regarded as a fear-induced shy behavior, and confirm that SERT has an important function in relation to risk aversion behaviors and coping style.

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).

Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots

The karyotype of most birds has remained considerably stable during more than 100 million years’ evolution, except for some groups, such as parrots. The evolutionary processes and underlying genetic mechanism of chromosomal rearrangements in parrots, however, are poorly understood. Here, using chromosome-level assemblies of four parrot genomes, we uncover frequent chromosome fusions and fissions, with most of them occurring independently among lineages. The increased activities of chromosomal rearrangements in parrots are likely associated with parrot-specific loss of two genes, ALC1 and PARP3, that have known functions in the repair of double-strand breaks and maintenance of genome stability. We further find that the fusion of the ZW sex chromosomes and chromosome 11 has created a pair of neo-sex chromosomes in the ancestor of parrots, and the chromosome 25 has been further added to the sex chromosomes in monk parakeet. Together, the combination of our genomic and cytogenetic analyses characterizes the complex evolutionary history of chromosomal rearrangements and sex chromosomes in parrots.

Parrots have undergone substantial karyotype evolution compared to most other birds. Here, Huang et al. analyze chromosome-level genome assemblies for four parrot species and elucidate the complex evolutionary history of parrot chromosomes.

Genomic architecture of phenotypic extremes in a wild cervid

Identifying the genes underlying fitness-related traits such as body size and male ornamentation can provide tools for conservation and management and are often subject to various selective pressures. Here we performed high-depth whole genome re-sequencing of pools of individuals representing the phenotypic extremes for antler and body size in white-tailed deer (Odocoileus virginianus). Samples were selected from a tissue repository containing phenotypic data for 4,466 male white-tailed deer from Anticosti Island, Quebec, with four pools representing the extreme phenotypes for antler and body size after controlling for age. Our results revealed a largely homogenous population but detected highly divergent windows between pools for both traits, with the mean allele frequency difference of 14% for and 13% for antler and body SNPs in outlier windows, respectively. Genes in outlier antler windows were enriched for pathways associated with cell death and protein metabolism and some of the most differentiated windows included genes associated with oncogenic pathways and reproduction, processes consistent with antler evolution and growth. Genes associated with body size were more nuanced, suggestive of a highly complex trait. Overall, this study revealed the complex genomic make-up of both antler morphology and body size in free-ranging white-tailed deer and identified target loci for additional analyses.

Wild epigenetics: insights from epigenetic studies on natural populations

Epigenetic mechanisms such as DNA methylation, histone modifications and non-coding RNAs are increasingly targeted in studies of natural populations. Here, I review some of the insights gained from this research, examine some of the methods currently in use and discuss some of the challenges that researchers working on natural populations are likely to face when probing epigenetic mechanisms. While studies supporting the involvement of epigenetic mechanisms in generating phenotypic variation in natural populations are amassing, many of these studies are currently correlative in nature. Thus, while empirical data point to widespread contributions of epigenetic mechanisms in generating phenotypic variation, there are still concerns as to whether epigenetic variation is instead ultimately controlled by genetic variation. Disentangling these two sources of variation will be a key to resolving the debate about the importance of epigenetic mechanisms, and studies on natural populations that partition the relative contribution of genetic and epigenetic factors to phenotypic variation can play an important role in this debate.

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.

Local selection signals in the genome of Blue tits emphasize regulatory and neuronal evolution

Understanding the genomic landscape of adaptation is central to the understanding of microevolution in wild populations. Genomic targets of selection and the underlying genomic mechanisms of adaptation can be elucidated by genome-wide scans for past selective sweeps or by scans for direct fitness associations. We sequenced and assembled 150 haplotypes of 75 Blue tits (Cyanistes caeruleus) of a single central-European population by a linked-read technology. We used these genome data in combination with coalescent simulations (1) to estimate an historical effective population size of ~250,000, which recently declined to ~10,000, and (2) to identify genome-wide distributed selective sweeps of beneficial variants most likely originating from standing genetic variation (soft sweeps). The genes linked to these soft sweeps, but also the ones linked to hard sweeps based on new beneficial mutants, showed a significant enrichment for functions associated with gene expression and transcription regulation. This emphasizes the importance of regulatory evolution in the population's adaptive history. Soft sweeps were further enriched for genes related to axon and synapse development, indicating the significance of neuronal connectivity changes in the brain potentially linked to behavioural adaptations. A previous scan of heterozygosity-fitness correlations revealed a consistent negative effect on arrival date at the breeding site for a single microsatellite in the MDGA2 gene. Here, we used the haplotype structure around this microsatellite to explain the effect as a local and direct outbreeding effect of a gene involved in synapse development.

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.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

Avian Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and W Degeneration

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed "neo-sex chromosomes" consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Exploring potentialities of avian genomic research in Nepalese Himalayas

Nepal, a small landlocked country in South Asia, holds about 800 km of Himalayan Mountain range including the Earth's highest mountain. Within such a mountain range in the north and plain lowlands in the south, Nepal provides a habitat for about 9% of global avian fauna. However, this diversity is underrated because of the lack of enough studies, especially using molecular tools to quantify and understand the distribution patterns of diversity. In this study, we reviewed the studies in the last two decades (2000‒2019) that used molecular methods to study the biodiversity in Nepal to examine the ongoing research trend and focus. Although Nepalese Himalaya has many opportunities for cutting-edge molecular research, our results indicated that the rate of genetic/genomic studies is much slower compared to the regional trends. We found that genetic research in Nepal heavily relies on resources from international institutes and that too is mostly limited to research on species monitoring, distribution, and taxonomic validations. Local infrastructures to carry out cutting-edge genomic research in Nepal are still in their infancy and there is a strong need for support from national/international scientists, universities, and governmental agencies to expand such genomic infrastructures in Nepal. We particularly highlight avian fauna as a potential future study system in this region that can be an excellent resource to explore key biological questions such as understanding eco-physiology and molecular basis of organismal persistence to changing environment, evolutionary processes underlying divergence and speciation, or mechanisms of endemism and restrictive distribution of species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s40657-021-00290-5.

A Chromosome-Level Genome Assembly of the Reed Warbler (Acrocephalus scirpaceus)

The reed warbler (Acrocephalus scirpaceus) is a long-distance migrant passerine with a wide distribution across Eurasia. This species has fascinated researchers for decades, especially its role as host of a brood parasite, and its capacity for rapid phenotypic change in the face of climate change. Currently, it is expanding its range northwards in Europe, and is altering its migratory behavior in certain areas. Thus, there is great potential to discover signs of recent evolution and its impact on the genomic composition of the reed warbler. Here, we present a high-quality reference genome for the reed warbler, based on PacBio, 10×, and Hi-C sequencing. The genome has an assembly size of 1,075,083,815 bp with a scaffold N50 of 74,438,198 bp and a contig N50 of 12,742,779 bp. BUSCO analysis using aves_odb10 as a model showed that 95.7% of BUSCO genes were complete. We found unequivocal evidence of two separate macrochromosomal fusions in the reed warbler genome, in addition to the previously identified fusion between chromosome Z and a part of chromosome 4A in the Sylvioidea superfamily. We annotated 14,645 protein-coding genes, and a BUSCO analysis of the protein sequences indicated 97.5% completeness. This reference genome will serve as an important resource, and will provide new insights into the genomic effects of evolutionary drivers such as coevolution, range expansion, and adaptations to climate change, as well as chromosomal rearrangements in birds.

Connecting the data landscape of long-term ecological studies: The SPI-Birds data hub

The integration and synthesis of the data in different areas of science is drastically slowed and hindered by a lack of standards and networking programmes. Long-term studies of individually marked animals are not an exception. These studies are especially important as instrumental for understanding evolutionary and ecological processes in the wild. Furthermore, their number and global distribution provides a unique opportunity to assess the generality of patterns and to address broad-scale global issues (e.g. climate change). To solve data integration issues and enable a new scale of ecological and evolutionary research based on long-term studies of birds, we have created the SPI-Birds Network and Database (www.spibirds.org)-a large-scale initiative that connects data from, and researchers working on, studies of wild populations of individually recognizable (usually ringed) birds. Within year and a half since the establishment, SPI-Birds has recruited over 120 members, and currently hosts data on almost 1.5 million individual birds collected in 80 populations over 2,000 cumulative years, and counting. SPI-Birds acts as a data hub and a catalogue of studied populations. It prevents data loss, secures easy data finding, use and integration and thus facilitates collaboration and synthesis. We provide community-derived data and meta-data standards and improve data integrity guided by the principles of Findable, Accessible, Interoperable and Reusable (FAIR), and aligned with the existing metadata languages (e.g. ecological meta-data language). The encouraging community involvement stems from SPI-Bird's decentralized approach: research groups retain full control over data use and their way of data management, while SPI-Birds creates tailored pipelines to convert each unique data format into a standard format. We outline the lessons learned, so that other communities (e.g. those working on other taxa) can adapt our successful model. Creating community-specific hubs (such as ours, COMADRE for animal demography, etc.) will aid much-needed large-scale ecological data integration.

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.

Most Genomic Loci Misrepresent the Phylogeny of an Avian Radiation Because of Ancient Gene Flow

Phylogenetic trees based on genome-wide sequence data may not always represent the true evolutionary history for a variety of reasons. One process that can lead to incorrect reconstruction of species phylogenies is gene flow, especially if interspecific gene flow has affected large parts of the genome. We investigated phylogenetic relationships within a clade comprising eight species of passerine birds (Phylloscopidae, Phylloscopus, leaf warblers) using one de novo genome assembly and 78 resequenced genomes. On the basis of hypothesis-exclusion trials based on D-statistics, phylogenetic network analysis, and demographic inference analysis, we identified ancient gene flow affecting large parts of the genome between one species and the ancestral lineage of a sister species pair. This ancient gene flow consistently caused erroneous reconstruction of the phylogeny when using large amounts of genome-wide sequence data. In contrast, the true relationships were captured when smaller parts of the genome were analyzed, showing that the "winner-takes-all democratic majority tree" is not necessarily the true species tree. Under this condition, smaller amounts of data may sometimes avoid the effects of gene flow due to stochastic sampling, as hidden reticulation histories are more likely to emerge from the use of larger data sets, especially whole-genome data sets. In addition, we also found that genomic regions affected by ancient gene flow generally exhibited higher genomic differentiation but a lower recombination rate and nucleotide diversity. Our study highlights the importance of considering reticulation in phylogenetic reconstructions in the genomic era.[Bifurcation; introgression; recombination; reticulation; Phylloscopus.].

The Genomic Signature of Allopatric Speciation in a Songbird Is Shaped by Genome Architecture (Aves: Certhia americana)

The genomic signature of speciation with gene flow is often attributed to the strength of divergent selection and recombination rate in regions harboring targets for selection. In contrast, allopatric speciation provides a different geographic context and evolutionary scenario, whereby introgression is limited by isolation rather than selection against gene flow. Lacking shared divergent selection or selection against hybridization, we would predict the genomic signature of allopatric speciation would largely be shaped by genomic architecture-the nonrandom distribution of functional elements and chromosomal characteristics-through its role in affecting the processes of selection and drift. Here, we built and annotated a chromosome-scale genome assembly for a songbird (Passeriformes: Certhia americana). We show that the genomic signature of allopatric speciation between its two primary lineages is largely shaped by genomic architecture. Regionally, gene density and recombination rate variation explain a large proportion of variance in genomic diversity, differentiation, and divergence. We identified a heterogeneous landscape of selection and neutrality, with a large portion of the genome under the effects of indirect selection. We found higher proportions of small chromosomes under the effects of indirect selection, likely because they have relatively higher gene density. At the chromosome scale, differential genomic architecture of macro- and microchromosomes shapes the genomic signatures of speciation: chromosome size has: 1) a positive relationship with genetic differentiation, genetic divergence, rate of lineage sorting in the contact zone, and proportion neutral evolution and 2) a negative relationship with genetic diversity and recombination rate.

The time is ripe for functional genomics: Can epigenetic changes mediate reproductive timing?

Populations are under strong selection to match reproductive timing with favourable environmental conditions. This becomes particularly important and challenging with increasing interannual environmental variability. Adjusting reproductive timing requires the ability to sense and interpret relevant environmental cues, while responding flexibly to their interannual variation. For instance, in seasonal species, reproductive timing is often dependent on photoperiod and temperature. Although many genes influencing the timing of reproduction have been identified, far less attention has been paid to the gene-regulatory cascades orchestrating these complex gene-environment interactions. In a From the Cover article in this issue of Molecular Ecology, Lindner, Laine, et al. (2021) addressed this knowledge gap by investigating the role of DNA methylation in mediating reproductive timing in the seasonally breeding great tit (Parus major). Using a clever blood sampling design, they investigated genome-wide DNA methylation changes following individual female birds across multiple reproductive stages. This approach revealed 10 candidate genes with a strong correlation between promoter methylation and reproductive status. Some of these genes are known to be involved in reproductive timing (e.g., MYLK-like or NR5A1), yet for others this function was previously unknown (Figure 1). Interestingly, NR5A1 is a key transcription factor, which may affect other genes that are part of the same regulatory network. The findings of Lindner, Laine, et al. (2021) provide a strong case for studying DNA methylation to uncover how gene-environment interactions influence important life-history traits, such as reproductive timing.

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.

Continent-wide genomic signatures of adaptation to urbanisation in a songbird across Europe

Urbanisation is increasing worldwide, and there is now ample evidence of phenotypic changes in wild organisms in response to this novel environment. Yet, the genetic changes and genomic architecture underlying these adaptations are poorly understood. Here, we genotype 192 great tits (Parus major) from nine European cities, each paired with an adjacent rural site, to address this major knowledge gap in our understanding of wildlife urban adaptation. We find that a combination of polygenic allele frequency shifts and recurrent selective sweeps are associated with the adaptation of great tits to urban environments. While haplotypes under selection are rarely shared across urban populations, selective sweeps occur within the same genes, mostly linked to neural function and development. Collectively, we show that urban adaptation in a widespread songbird occurs through unique and shared selective sweeps in a core-set of behaviour-linked genes.

Evolution of the "world's only alpine parrot": Genomic adaptation or phenotypic plasticity, behaviour and ecology?

Climate warming, in particular in island environments, where opportunities for species to disperse are limited, may become a serious threat to cold adapted alpine species. In order to understand how alpine species may respond to a warming world, we need to understand the drivers that have shaped their habitat specialisation and the evolutionary adaptations that allow them to utilize alpine habitats. The endemic, endangered New Zealand kea (Nestor notabilis) is considered the only alpine parrot in the world. As a species commonly found in the alpine zone it may be highly susceptible to climate warming. But is it a true alpine specialist? Is its evolution driven by adaptation to the alpine zone, or is the kea an open habitat generalist that simply uses the alpine zone to, for example, avoid lower lying anthropogenic landscapes? We use whole genome data of the kea and its close, forest adapted sister species, the kākā (Nestor meridionalis) to reconstruct the evolutionary history of both species and identify the functional genomic differences that underlie their habitat specialisations. Our analyses do not identify major functional genomic differences between kea and kākā in pathways associated with high-altitude. Rather, we found evidence that selective pressures on adaptations commonly found in alpine species are present in both Nestor species, suggesting that selection for alpine adaptations has not driven their divergence. Strongly divergent demographic responses to past climate warming between the species nevertheless highlight potential future threats to kea survival in a warming world.

Early acoustic experience alters genome-wide methylation in the auditory forebrain of songbird embryos

Early exposure to salient cues can critically shape the development of social behaviors. For example, both oscine birds and humans can hear and learn to recognize familiar sounds in ovo and in utero and recognize them following hatching and birth, respectively. Here we demonstrate that different chronic acoustic playbacks alter genome-wide methylation of the auditory forebrain in late-stage zebra finch (Taeniopygia guttata) embryos. Within the same subjects, immediate early gene activation in response to acute con- or heterospecific song exposure is negatively correlated with methylation extent in response to repeated daily prior exposure to the same type of stimuli. Specifically, we report less relative global methylation following playbacks of conspecific songs and more methylation following playbacks of distantly-related heterospecific songs. These findings offer a neuroepigenomic mechanism for the ontogenetic impacts of early acoustic experiences in songbirds.

Prenatal testosterone triggers long-term behavioral changes in male zebra finches: unravelling the neurogenomic mechanisms

BACKGROUND

Maternal hormones, like testosterone, can strongly influence developing offspring, even generating long-term organizational effects on adult behavior; yet, the mechanisms facilitating these effects are still unclear. Here, we experimentally elevated prenatal testosterone in the eggs of zebra finches (Taeniopygia guttata) and measured male aggression in adulthood along with patterns of neural gene expression (RNA-seq) and DNA methylation (MethylC-Seq) in two socially relevant brain regions (hypothalamus and nucleus taenia of the amygdala). We used enrichment analyses and protein-protein interaction networks to find candidate processes and hub genes potentially affected by the treatment. We additionally identified differentially expressed genes that contained differentially methylated regions.

RESULTS

We found that males from testosterone-injected eggs displayed more aggressive behaviors compared to males from control eggs. Hundreds of genes were differentially expressed, particularly in the hypothalamus, including potential aggression-related hub genes (e.g., brain derived neurotrophic factor). There were also enriched processes with well-established links to aggressive phenotypes (e.g., somatostatin and glutamate signaling). Furthermore, several highly connected genes identified in protein-protein interaction networks also showed differential methylation, including adenylate cyclase 2 and proprotein convertase 2.

CONCLUSIONS

These results highlight genes and processes that may play an important role in mediating the effects of prenatal testosterone on long-term phenotypic outcomes, thereby providing insights into the molecular mechanisms that facilitate hormone-mediated maternal effects.

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.

Effects of maternal exposure to a bacterial antigen and altered post-hatching rearing conditions on avian offspring behaviour

Abstract

The early-life environment plays a crucial role in shaping morphological, physiological, and behavioural traits, with potential long-term consequences for fitness. Indeed, a set of factors experienced by offspring during prenatal and early post-natal development has been recognised to affect behavioural trait expression in later life. Several studies have shown that in birds, nutritional and social rearing conditions and maternal and/or neonatal immunisations may profoundly determine the development and establishment of behaviour in offspring. To our knowledge, no research has examined whether and how the interaction between immune-mediated maternal effects and post-hatching rearing conditions affects offspring behaviour. Here, we studied the effects of maternal exposure to a bacterial antigen and altered brood size on docility, breathing rate, and aggression in the offspring of great tit, Parus major. We used a 2 × 2 design to investigate the interactive effects of maternal immunisation and brood size manipulation on offspring behavioural development. We found no such interactive effect on offspring behaviour, although we observed it regarding to offspring body mass and tarsus length. Maternal immunisation itself did not affect offspring behaviour. However, we demonstrated that the offspring breathing rate and level of aggression were affected by brood size manipulation. Both breathing rate and aggression in offspring reared in enlarged broods were lower than those in offspring reared in non-manipulated broods. Our study did not confirm earlier reports that immune-mediated maternal effects modulate offspring behavioural development, but we showed that brood size during rearing might indeed be a factor that affects offspring behaviour.

Significance statement

The early environment experienced by offspring constitutes a significant source of developmental plasticity, which may profoundly affect the establishment of their behavioural traits. Food availability, social conditions, and maternal or offspring infection are crucial factors shaping various behavioural traits in birds. However, there remains a lack of studies emphasising the potential interactive effects of early-life conditions on behavioural trait development in natural bird populations. Here, to our knowledge for the first time, we experimentally examined how maternal immunisation and altered post-hatching rearing conditions interact to determine the behaviour of fledged offspring. We found that maternal treatment and brood size manipulation interactively affected offspring body mass and tarsus length, but this interaction had no effect on offspring behaviour. Our findings suggest that different mechanisms may underlie the development of morphological and behavioural traits.

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.

A Study of Faster-Z Evolution in the Great Tit (Parus major)

Sex chromosomes contribute substantially to key evolutionary processes such as speciation and adaptation. Several theories suggest that evolution could occur more rapidly on sex chromosomes, but currently our understanding of whether and how this occurs is limited. Here, we present an analysis of the great tit (Parus major) genome, aiming to detect signals of faster-Z evolution. We find mixed evidence of faster divergence on the Z chromosome than autosomes, with significantly higher divergence being found in ancestral repeats, but not at 4- or 0-fold degenerate sites. Interestingly, some 4-fold sites appear to be selectively constrained, which may mislead analyses that use these sites as the neutral reference (e.g., d_N/d_S). Consistent with other studies in birds, the mutation rate is significantly higher in males than females, and the long-term Z-to-autosome effective population size ratio is only 0.5, significantly lower than the expected value of 0.75. These are indicative of male-driven evolution and high variance in male reproductive success, respectively. We find no evidence for an increased efficacy of positive selection on the Z chromosome. In contrast, the Z chromosome in great tits appears to be affected by increased genetic drift, which has led to detectable signals of weakened intensity of purifying selection. These results provide further evidence that the Z chromosome often has a low effective population size, and that this has important consequences for its evolution. They also highlight the importance of considering multiple factors that can affect the rate of evolution and effective population sizes of sex chromosomes.

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

Genome-wide variation in DNA methylation linked to developmental stage and chromosomal suppression of recombination in white-throated sparrows

Much of our knowledge on regulatory impacts of DNA methylation has come from laboratory‐bred model organisms, which may not exhibit the full extent of variation found in wild populations. Here, we investigated naturally‐occurring variation in DNA methylation in a wild avian species, the white‐throated sparrow (Zonotrichia albicollis). This species offers exceptional opportunities for studying the link between genetic differentiation and phenotypic traits because of a nonrecombining chromosome pair linked to both plumage and behavioural phenotypes. Using novel single‐nucleotide resolution methylation maps and gene expression data, we show that DNA methylation and the expression of DNA methyltransferases are significantly higher in adults than in nestlings. Genes for which DNA methylation varied between nestlings and adults were implicated in development and cell differentiation and were located throughout the genome. In contrast, differential methylation between plumage morphs was concentrated in the nonrecombining chromosome pair. Interestingly, a large number of CpGs on the nonrecombining chromosome, localized to transposable elements, have undergone dramatic loss of DNA methylation since the split of the ZAL2 and ZAL2^m chromosomes. Changes in methylation predicted changes in gene expression for both chromosomes. In summary, we demonstrate changes in genome‐wide DNA methylation that are associated with development and with specific functional categories of genes in white‐throated sparrows. Moreover, we observe substantial DNA methylation reprogramming associated with the suppression of recombination, with implications for genome integrity and gene expression divergence. These results offer an unprecedented view of ongoing epigenetic reprogramming in a wild population.

see also the Perspective by Jordan A. Anderson and Jenny Tung.