Sex-Specific Methylomic and Transcriptomic Responses of the Avian Pineal Gland to Unpredictable Illumination Patterns

In the production environment of chickens, exposure to unpredictable light patterns is a common painless stressor, widely used to influence growth rate and egg production efficiency. The pineal gland, a key regulator of circadian rhythms through melatonin secretion, responds to environmental light cues, and its function is modulated by epigenetic mechanisms. In this study, we investigated how the pineal gland methylome and transcriptome (including micro-RNAs) interact to respond to a rearing exposure to unpredictable illumination patterns, with a particular focus on sex differences. We conducted an integrative multi-omic analysis-including methylomic (MeDIP-seq), transcriptomic (RNA-seq), and miRNA expression profiling-on the pineal gland of Hy-Line White chickens (n = 34, 18 females, 16 males) exposed to either a standard 12:12 light-dark cycle (control) or a randomized, unpredictable light schedule from Days 3 to 24 post-hatch. Our findings reveal that unpredictable light exposure alters the pineal gland methylome and transcriptome in a sex-specific manner. However, while transcriptomic differences between sexes increased due to the stress, methylomic differences decreased, particularly on the Z chromosome. These changes were driven by females (the heterogametic sex in birds), which became more male-like in their pineal methylome after exposure to the illumination stress, leading to reduced epigenetic sexual dimorphism while maintaining differences at the gene expression level. Further, we implemented a fixed sex effect model as a biological proof of concept, identifying a network of 12 key core genes interacting with 102 other genes, all linked to circadian regulation and stress adaptation. This network of genes comprises a core regulatory framework for circadian response. Additionally, tissue-specific expression analysis and cell-type specific expression analysis revealed enrichment in brain regions critical for circadian function, including neuronal populations involved in circadian regulation and the hypothalamic-pituitary-thyroid axis. Together, these findings provide strong evidence of sex-specific epigenetic transcriptomic responses of the pineal gland upon illumination stress and offer valuable insights into the interplay of different omic levels in relation to circadian response.