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

Maladaptive learning and memory in hybrids as a reproductive isolating barrier

Selection against hybrid offspring, or postzygotic reproductive isolation, maintains species boundaries in the face of gene flow from hybridization. In this review, we propose that maladaptive learning and memory in hybrids is an important, but overlooked form of postzygotic reproductive isolation. Although a role for learning in premating isolation has been supported, whether learning deficiencies can contribute to postzygotic isolation has rarely been tested. We argue that the novel genetic combinations created by hybridization have the potential to impact learning and memory abilities through multiple possible mechanisms, and that any displacement from optima in these traits is likely to have fitness consequences. We review evidence supporting the potential for hybridization to affect learning and memory, and evidence of links between learning abilities and fitness. Finally, we suggest several avenues for future research. Given the importance of learning for fitness, especially in novel and unpredictable environments, maladaptive learning and memory in hybrids may be an increasingly important source of postzygotic reproductive isolation.

A high-density SNP chip for genotyping great tit (Parus major) populations and its application to studying the genetic architecture of exploration behaviour

High-density SNP microarrays ("SNP chips") are a rapid, accurate and efficient method for genotyping several hundred thousand polymorphisms in large numbers of individuals. While SNP chips are routinely used in human genetics and in animal and plant breeding, they are less widely used in evolutionary and ecological research. In this article, we describe the development and application of a high-density Affymetrix Axiom chip with around 500,000 SNPs, designed to perform genomics studies of great tit (Parus major) populations. We demonstrate that the per-SNP genotype error rate is well below 1% and that the chip can also be used to identify structural or copy number variation. The chip is used to explore the genetic architecture of exploration behaviour (EB), a personality trait that has been widely studied in great tits and other species. No SNPs reached genomewide significance, including at DRD4, a candidate gene. However, EB is heritable and appears to have a polygenic architecture. Researchers developing similar SNP chips may note: (i) SNPs previously typed on alternative platforms are more likely to be converted to working assays; (ii) detecting SNPs by more than one pipeline, and in independent data sets, ensures a high proportion of working assays; (iii) allele frequency ascertainment bias is minimized by performing SNP discovery in individuals from multiple populations; and (iv) samples with the lowest call rates tend to also have the greatest genotyping error rates.

Transgenerational epigenetic inheritance in birds

While it has been shown that epigenetics accounts for a portion of the variability of complex traits linked to interactions with the environment, the real contribution of epigenetics to phenotypic variation remains to be assessed. In recent years, a growing number of studies have revealed that epigenetic modifications can be transmitted across generations in several animal species. Numerous studies have demonstrated inter- or multi-generational effects of changing environment in birds, but very few studies have been published showing epigenetic transgenerational inheritance in these species. In this review, we mention work conducted in parent-to-offspring transmission analyses in bird species, with a focus on the impact of early stressors on behaviour. We then present recent advances in transgenerational epigenetics in birds, which involve germline linked non-Mendelian inheritance, underline the advantages and drawbacks of working on birds in this field and comment on future directions of transgenerational studies in bird species.

CNVs are associated with genomic architecture in a songbird

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

SERT gene polymorphisms are associated with risk-taking behaviour and breeding parameters in wild great tits

Individual differences in coping with potentially dangerous situations are affected by a combination of genetic and environmental factors. How genetic polymorphisms and behavioural variations are related to fitness is unknown. One of the candidate genes affecting a variety of behavioural processes, including impulsivity, anxiety and mood fluctuations in both humans and other vertebrates, is the serotonin transporter gene (SERT/SLC6A). The aim of this study was to assess an association between SERT genotypes and novelty-seeking and risk-taking behaviours as well as breeding parameters of great tits (Parus major) in a natural environment. We associated polymorphisms in the promoter exonic regions of the SERT gene with parental risk-taking-related behaviour and fitness traits. Our results show that: (1) risk-taking behaviour in our great tit population is linked to single nucleotide polymorphisms in the SERT gene exon 3 and exon 8; (2) the genotype-behaviour associations are consistent with the presence of different stressors; and (3) polymorphisms in exon 8 could be associated with fitness-related traits, such as the start of egg-laying and hatching success. We showed for the first time that genetic variability of SERT plays an important role in shaping individual decision-making that affects fitness in a wild population. However, the results are based on one population and on the polymorphisms that are in a single gene. Therefore, replication studies are needed in order to confirm these preliminary results.

Rapid regulatory evolution of a nonrecombining autosome linked to divergent behavioral phenotypes

In the white-throated sparrow (Zonotrichia albicollis), the second chromosome bears a striking resemblance to sex chromosomes. First, within each breeding pair of birds, one bird is homozygous for the standard arrangement of the chromosome (ZAL2/ZAL2) and its mate is heterozygous for a different version (ZAL2/ZAL2^m). Second, recombination is profoundly suppressed between the two versions, leading to genetic differentiation between them. Third, the ZAL2^m version is linked with phenotypic traits, such as bright plumage, high aggression, and low parental behavior, which are usually associated with males. These similarities to sex chromosomes suggest that the evolutionary mechanisms that shape sex chromosomes, in particular genetic degeneration of the heterogametic version due to the suppression of recombination, are likely important in this system as well. Here, we investigated patterns of protein sequence evolution and gene expression evolution between the ZAL2 and ZAL2^m chromosomes by whole-genome sequencing and transcriptome analyses. Patterns of protein evolution exhibited only weak signals of genetic degeneration, and few genes harbored signatures of positive selection. We found substantial evidence of transcriptome evolution, such as significant expression divergence between ZAL2 and ZAL2^m alleles and signatures of dosage compensation for highly expressed genes. These results suggest that, early in the evolution of heteromorphic chromosomes, gene expression divergence and dosage compensation can prevail before large-scale genetic degeneration. Our results show further that suppression of recombination between heteromorphic chromosomes can lead to the evolution of alternative (sex-like) behavioral phenotypes before substantial genetic degeneration.

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

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

DNA methylation changes and evolution of RNA-based duplication in Sus scrofa: based on a two-step strategy

AIM

This study aims to couple DNA methylation changes and evolution of retrogenes.

MATERIALS & METHODS

A new two-step strategy was developed to screen retrogenes. Further, reduced representation bisulfite sequencing and RNA-seq data of eight tissues were used to analyze retrogenes.

RESULTS

A total of 964 retrocopies were identified and new retrocopies were available for the synthesis of glycans and lipids corresponding to pig phenotypic traits. Retrogenes were consistently hypermethylated. Hypomethylation of parental genes presented more susceptibility to retroposition. Promoter DNA methylation of retrogenes was negatively correlated with evolutionary time and played important roles in regulating retrogene tissue-specific expression pattern.

CONCLUSION

A two-step procedure is effective and necessary for identifying retrogenes. DNA methylation drives origination, survival, evolution and expression of retrogenes.

Recent natural selection causes adaptive evolution of an avian polygenic trait

We used extensive data from a long-term study of great tits (Parus major) in the United Kingdom and Netherlands to better understand how genetic signatures of selection translate into variation in fitness and phenotypes. We found that genomic regions under differential selection contained candidate genes for bill morphology and used genetic architecture analyses to confirm that these genes, especially the collagen gene COL4A5, explained variation in bill length. COL4A5 variation was associated with reproductive success, which, combined with spatiotemporal patterns of bill length, suggested ongoing selection for longer bills in the United Kingdom. Last, bill length and COL4A5 variation were associated with usage of feeders, suggesting that longer bills may have evolved in the United Kingdom as a response to supplementary feeding.

Evolution of beak morphology in the Ground Tit revealed by comparative transcriptomics

Background

Beak morphology exhibits considerable adaptive plasticity in birds, which results in highly varied or specialized forms in response to variations in ecology and life history. As the only parid species endemic to the Qinghai-Tibet Plateau, the Ground Tit (Parus humilis) has evolved a distinctly long and curved beak from other parids. An integration of morphometrics, phylogenetics, transcriptomics and embryology allows us to address the evolutionary and developmental mechanisms of the adaptive beak structure observed in the Ground Tit.

Results

A morphometric approach quantified that the Ground Tit has a comparatively longer and more decurved upper beaks than other parids. We estimated that the ancestor of the Ground Tit likely had a short straight upper beak similar to most current recognized parid species using an ancestral state reconstruction. This morphological specialization is considered an adaptation to its ground-oriented behavior on the high plateau. To identify genetic mechanisms behind this adaptive change, a comparative transcriptomic analysis was applied between the Ground Tit and its closely related species, the Great Tit (Parus major). We detected that 623 genes were significantly differentially expressed in embryonic upper beaks between the two species, 17 of which were functionally annotated to correlate with bone development and morphogenesis, although genes related to bone development were not found to undergo accelerated evolution in the Ground Tit. RT-qPCR validation confirmed differential expression of five out of eight genes that were selected from the 17 genes. Subsequent functional assays in chicken embryos demonstrated that two of these genes, FGF13 and ITGB3, may affect beak morphology by modulating levels of osteoblasts and osteoclasts.

Conclusions

Our results provide preliminary evidence that development of the long decurved beak of the Ground Tit is likely regulated by transcriptional activities of multiple genes coordinating osteoblasts and osteoclasts. The integration of multiple approaches employed here sheds light on ecological and genetic mechanisms in the evolution of avian morphology.

Characterization of a de novo assembled transcriptome of the Common Blackbird (Turdus merula)

BACKGROUND

In recent years, next generation high throughput sequencing technologies have proven to be useful tools for investigations concerning the genomics or transcriptomics also of non-model species. Consequently, ornithologists have adopted these technologies and the respective bioinformatics tools to survey the genomes and transcriptomes of a few avian non-model species. The Common Blackbird is one of the most common bird species living in European cities, which has successfully colonized urban areas and for which no reference genome or transcriptome is publicly available. However, to target questions like genome wide gene expression analysis, a reference genome or transcriptome is needed.

METHODS

Therefore, in this study two Common Blackbirds were sacrificed, their mRNA was isolated and analyzed by RNA-Seq to de novo assemble a transcriptome and characterize it. Illumina reads (125 bp paired-end) and a Velvet/Oases pipeline led to 162,158 transcripts. For the annotation (using Blast+), an unfiltered protein database was used. SNPs were identified using SAMtools and BCFtools. Furthermore, mRNA from three single tissues (brain, heart and liver) of the same two Common Blackbirds were sequenced by Illumina (75 bp single-end reads). The draft transcriptome and the three single tissues were compared by their BLAST hits with the package VennDiagram in R.

RESULTS

Following the annotation against protein databases, we found evidence for 15,580 genes in the transcriptome (all well characterized hits after annotation). On 18% of the assembled transcripts, 144,742 SNPs were identified which are, consequently, 0.09% of all nucleotides in the assembled transcriptome. In the transcriptome and in the single tissues (brain, heart and liver), 10,182 shared genes were found.

DISCUSSION

Using a next-generation technology and bioinformatics tools, we made a first step towards the genomic investigation of the Common Blackbird. The de novo assembled transcriptome is usable for downstream analyses such as differential gene expression analysis and SNP identification. This study shows the importance of the approach to sequence single tissues to understand functions of tissues, proteins and the phenotype.

Genomic differentiation and patterns of gene flow between two long-tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions ("genomic differentiation islands") spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene-flow patterns for autosomes using RAD-seq data between two closely related species of long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.

Environment-Dependent Genotype-Phenotype Associations in Avian Breeding Time

Understanding how genes shape phenotypes is essential to assess the evolutionary potential of a trait. Identifying the genes underlying quantitative behavioral or life-history traits has, however, proven to be a major challenge. The majority of these traits are phenotypically plastic and different parts of the genome can be involved in shaping the trait under different environmental conditions. These variable genotype-phenotype associations could be one explanation for the limited success of genome-wide association studies in such traits. We here use avian seasonal timing of breeding, a trait that is highly plastic in response to spring temperature, to explore effects of such genotype-by-environment interactions in genome-wide association studies. We genotyped 2045 great tit females for 384081 single nucleotide polymorphisms (SNPs) and recorded their egg-laying dates in the wild. When testing for associations between SNPs and egg-laying dates, no SNP reached genome-wide significance. We then explored whether SNP effects were modified by annual spring temperature by formally testing for an interaction between SNP effect and temperature. The models including the SNP^∗temperature interaction performed consistently better although no SNP reached genome-wide significance. Our results suggest that the effects of genes shaping seasonal timing depended on annual spring temperature. Such environment-dependent effects are expected for any phenotypically plastic trait. Taking these effects into account will thus improve the success of detecting genes involved in phenotypically plastic traits, thereby leading to a better understanding of their evolutionary potential.

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

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

Sociomics: Using Omic Approaches to Understand Social Evolution

All of life is social, from genes cooperating to form organisms, to animals cooperating to form societies. Omic approaches offer exceptional opportunities to solve major outstanding problems in the study of how sociality evolves. First, omics can be used to clarify the extent and form of sociality in natural populations. This is especially useful in species where it is difficult to study social traits in natural populations, such as bacteria and other microbes. Second, omics can be used to examine the consequences of sociality for genome evolution and gene expression. This is especially useful in cases where there is clear variation in the level of sociality, such as the social insects. Major tasks for the future are to apply these approaches to a wider range of non-model organisms, and to move from exploratory analyses to the testing of evolutionary theory.

Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.

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

BACKGROUND

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

METHODS

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

CONCLUSIONS

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

Transcriptome analysis of a wild bird reveals physiological responses to the urban environment

Identifying the molecular basis of environmentally induced phenotypic variation presents exciting opportunities for furthering our understanding of how ecological processes and the environment can shape the phenotype. Urban and rural environments present free-living organisms with different challenges and opportunities, which have marked consequences for the phenotype, yet little is known about responses at the molecular level. We characterised transcriptomes from an urban and a rural population of great tits Parus major, demonstrating striking differences in gene expression profiles in both blood and liver tissues. Differentially expressed genes had functions related to immune and inflammatory responses, detoxification, protection against oxidative stress, lipid metabolism, and regulation of gene expression. Many genes linked to stress responses were expressed at higher levels in the urban birds, in accordance with our prediction that urban animals are exposed to greater environmental stress. This is one of the first studies to reveal transcriptional differences between urban- and rural-dwelling animals and suggests an important role for epigenetics in mediating environmentally induced physiological variation. The study provides valuable resources for developing further in-depth studies of the mechanisms driving phenotypic variation in the urban context at larger spatial and temporal scales.

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

Methylation of the circadian Clock gene in the offspring of a free-living passerine bird increases with maternal and individual exposure to PM10

The consequences of exposure to particulate matter (PM) have been thoroughly investigated in humans and other model species, but there is a dearth of studies of the effects of PM on physiology and life-history traits of non-human organisms living in natural or semi-natural environments. Besides toxicological relevance, PM has been recently suggested to exert epigenetic effects by altering DNA methylation patterns. Here, we investigated for the first time the association between the exposure to free-air PM₁₀ and DNA methylation at two loci ('poly-Q exon' and '5'-UTR') of the Clock gene in blood cells of the nestlings of a synanthropic passerine bird, the barn swallow (Hirundo rustica). The Clock gene is a phylogenetically highly conserved gene playing a major role in governing circadian rhythms and circannual life cycles of animals, implying that change in its level of methylation can impact on important fitness traits. We found that methylation at both loci significantly increased with PM₁₀ levels recorded few days before blood sampling, and also with PM₁₀ exposure experienced by the mother during or shortly before egg laying. This study is the first where methylation at a functionally important gene has been shown to vary according to the concentration of anthropogenic pollutants in any animal species in the wild. Since early-life environmental conditions produce epigenetic effects that can transgenerationally be transmitted, DNA methylation of genes controlling photoperiodic response can have far reaching consequences for the ecology and the evolution of wild animal populations.

Impact of neonatal iron deficiency on hippocampal DNA methylation and gene transcription in a porcine biomedical model of cognitive development

BACKGROUND

Iron deficiency is a common childhood micronutrient deficiency that results in altered hippocampal function and cognitive disorders. However, little is known about the mechanisms through which neonatal iron deficiency results in long lasting alterations in hippocampal gene expression and function. DNA methylation is an epigenetic mark involved in gene regulation and altered by environmental factors. In this study, hippocampal DNA methylation and gene expression were assessed via reduced representation bisulfite sequencing and RNA-seq on samples from a previous study reporting reduced hippocampal-based learning and memory in a porcine biomedical model of neonatal iron deficiency.

RESULTS

In total 192 differentially expressed genes (DEGs) were identified between the iron deficient and control groups. GO term and pathway enrichment analysis identified DEGs associated with hypoxia, angiogenesis, increased blood brain barrier (BBB) permeability, and altered neurodevelopment and function. Of particular interest are genes previously implicated in cognitive deficits and behavioral disorders in humans and mice, including HTR2A, HTR2C, PAK3, PRSS12, and NETO1. Altered genome-wide DNA methylation was observed across 0.5 million CpG and 2.4 million non-CpG sites. In total 853 differentially methylated (DM) CpG and 99 DM non-CpG sites were identified between groups. Samples clustered by group when comparing DM non-CpG sites, suggesting high conservation of non-CpG methylation in response to neonatal environment. In total 12 DM sites were associated with 9 DEGs, including genes involved in angiogenesis, neurodevelopment, and neuronal function.

CONCLUSIONS

Neonatal iron deficiency leads to altered hippocampal DNA methylation and gene regulation involved in hypoxia, angiogenesis, increased BBB permeability, and altered neurodevelopment and function. Together, these results provide new insights into the mechanisms through which neonatal iron deficiency results in long lasting reductions in cognitive development in humans.

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

BACKGROUND

Studies on vertebrate DNA methylomes have revealed a regulatory role of tissue specific DNA methylation in relation to gene expression. However, it is not well known how tissue-specific methylation varies between different functional and structural components of genes and genomes. Using whole-genome bisulfite sequencing data we here describe both CpG and non-CpG methylation profiles of whole blood and brain tissue in relation to gene features, CpG-islands (CGIs), transposable elements (TE), and their functional roles in an ecological model species, the great tit (Parus major).

RESULTS

We show that hypomethylation at the transcription start site (TSS) is enriched in genes with functional classes that relate directly to processes specific to each tissue type. We find that 6877 (~21 %) of the CGIs are differentially methylated between blood and brain, of which 1186 and 2055 are annotated to promoter and intragenic regions, respectively. We observe that CGI methylation in promoter regions is more conserved between tissues compared to CGI methylation in intra and inter-genic regions. Differentially methylated CGIs in promoter and intragenic regions are overrepresented in genomic loci linked to development, suggesting a distinct role for CGI methylation in regulating expression during development. Additionally, we find significant non-CpG methylation in brain but not in blood with a strong preference for methylation at CpA dinucleotide sites. Finally, CpG hypermethylation of TEs is significantly stronger in brain compared to blood, but does not correlate with TE activity. Surprisingly, TEs showed significant hypomethylation in non-CpG contexts which was negatively correlated with TE expression.

CONCLUSION

The discovery that TSS methylation levels are directly linked to functional classes related to each tissue provides new insights in the regulatory role of DNA-methylation patterns. The dominant sequence motifs for brain non-CpG methylation, similar to those found in mammals, suggests that a conserved non-CpG regulatory mechanism was already present in the amniote ancestor. The negative correlation between brain non-CpG methylation and TE activity (not found for CpG methylation) suggests that non-CpG is the dominant regulatory form of methylation in TE silencing.

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.