Transcriptomic and Chromatin Landscape Analysis Reveals That Involvement of Pituitary Level Transcription Factors Modulate Incubation Behaviors of Magang Geese

The incubation behavior of geese seriously affects their egg production performance. Studies on incubation behavior have identified functional genes, but the regulatory architecture relationship between functional genes and chromatin accessibility remains poorly understood. Here, we present an integrated analysis of open chromatin profiles and transcriptome to identify the cis-regulatory element and their potential transcription factors involved in regulating incubation behavior in goose pituitary. Assay for transposase-accessible chromatin sequencing (ATAC-seq) revealed that open chromatin regions increased in the pituitary during the transition from incubation behavior to laying. We identified 920 significant differential accessible regions (DARs) in the pituitary. Compared to the laying stage, most DARs had higher chromatin accessibility in the brooding stage. Motif analysis of open DARs showed that the most significant transcription factor (TF) occupied sites predominantly enriched in motifs binding to the RFX family (RFX5, RFX2, and RFX1). While the majority of TF motifs enriched under sites of the nuclear receptor (NR) family (ARE, GRE, and PGR) in closed DARs at the incubation behavior stage. Footprint analysis indicated that the transcription factor RFX family exhibited higher binding on chromatin at the brooding stage. To further elucidate the effect of changes in chromatin accessibility on gene expression levels, a comparison of the transcriptome revealed 279 differentially expressed genes (DEGs). The transcriptome changes were associated with processes of steroid biosynthesis. By integrating ATAC-seq and RNA-seq, few DARs directly affect incubation behavior by regulating the transcription levels of genes. Five DAR-related DEGs were found to be closely related to maintaining the incubation behavior in geese. Footprinting analysis revealed a set of transcription factors (RFX1, RFX2, RFX3, RFX5, BHLHA15, SIX1, and DUX) which displayed the highest activity at the brooding stage. SREBF2 was predicted to be the unique differentially expressed transcription factor whose mRNA level was down-regulated and enriched in hyper-accessible regions of PRL in the broody stage. In the present study, we comprehensively profiled the transcriptome and chromatin accessibility in the pituitary related to incubation behavior. Our findings provided insight into the identification and analysis of regulatory elements in goose incubation behavior. The epigenetic alterations profiled here can help decipher the epigenetic mechanisms that contribute to the regulation of incubation behavior in birds.