1. Introduction
Incubation behavior is an important and specifically evolved reproductive behavior in birds, which is critical for the survival of offspring and thus adult reproductive success. Incubation behaviors lead to ovarian atrophy and the cessation of egg production, which seriously affects the egg production performance of domestic birds. Incubation behavior has been effectively genetically eliminated underlying the high intensity of artificial selection pressures [1]. However, geese still stubbornly maintain incubation behaviors.
Goose reproductive behavior is driven by light stimulation and is tightly regulated by the hypothalamus-pituitary-gonadal axis (HPG) [2,3,4]. The precise coordination of hormonal regulation at all levels of the HPG is essential for the formation of a complete reproductive cycle in chickens. At the neuroendocrine level, proper expression of key reproductive hormones from different pituitary cell types is required for incubation behavior onset and maintenance. The fluctuation of GnRH-FSH/LH and VIP-PRL in the neuroendocrine system determines the formation of the laying and incubation cycle in chickens [5]. High levels of PRL and low concentrations of ovarian steroids are a hallmark of incubation behavior [6,7].
Incubation behavior can be induced via the combined action of estradiol, progesterone, and PRL in ovariectomized female avian [8]. It has been suggested that PRL, estradiol, and progesterone are essential for modulating the reproductive axis for the onset of incubation behavior. However, the molecular role of hormonal signaling is still unclear. The regulatory architecture of incubation behavior and connections between incubation behavior-related gene regulatory networks and behavior maintenance are poorly understood.
Epigenome changes were thought to affect the domesticated phenotype by influencing gene expression in birds [9,10]. The primary regulators of gene expression are transcription factors. Progressively, more and more evidence shows that transcription factors are involved in endocrinology regulation [11]. Transcription factors are the key drivers of cell function and phenotype via their binding to a specific regulatory sequence (cis-regulatory element, CRE). Various elements regulate the actions of genes and alter phenotypic variation [12]. Epigenome alterations lead to heritable phenotypic and transcriptomic changes by regulating transcription factor binding [13,14]. Transcriptome comparisons between laying and brooding chickens have provided potential mechanisms that contribute to incubation behavior. Many transcription factors are activated or inactivated to regulated various genes involved in maintaining low concentrations of gonadotropin and low steroid biosynthesis in the brooding stage, at the pituitary level [15,16]. These reports highlight the key role of transcription factors in the regulation of incubation behavior at transcriptomic and epigenomic levels.
The pituitary gland plays critical roles in regulating many key physiological functions related to growth and reproduction [17]. Previously, our results showed the highly divergent transcriptome in the pituitary between laying and broody stages. Similar results were obtained by others [16,18,19]. The pituitary is a key component of the endocrine system in the modulation of avian incubation behavior. Incubation behavior is related to the cell types transiting on the pituitary [20,21,22]. The differentiation of both lactotrophs and somatotrophs in the pituitary is strictly modulated by transcription factors (TF) [23]. Transcription factors influence the epigenetic priming of cells toward different cell fates [24]. Transcription factors are an important type of gene expression regulator that plays a central role in regulating various physiological processes such as pituitary cell fate and the secretion of cell-type-specific pituitary hormones [23,25].
Gene expression is controlled by accessible chromatin. Genes with chromatin accessibility are more likely to be regulated by TF and are differentially expressed at the mRNA level. We suspect that the potential mechanism of hormonal signaling in the formation of incubation behavior is related to chromatin accessibility. To characterize the chromatin state landscape of the pituitary of brooding geese, integrative transcriptome and chromatin landscape analysis was performed to identify cis-regulatory elements and their potential transcription factors involved in regulating goose incubation behavior. These results may provide new evidence concerning the correlations among chromatin accessibility, gene expression, and hormone signaling for goose incubation behavior.
2. Materials and Methods
2.1. Ethics Statement
This study was carried out in strict accordance with the Experimental Animal Committee of Zhongkai University of Agriculture and Engineering (NO. 20201226012). All efforts have been made to minimize animal suffering.
2.2. Animals and Sample Collection
The Chinese natively bred Magang goose was obtained from the Wen Jianmin goose farm (Panyu, Guangzhou, Guangdong, China). We used twelve female geese in the present study; six were brooding (broodiness group) and the other six were continuously laying (laying group). Pituitary samples were collected, and transcriptome sequencing (RNA-seq) was performed (n = 3 for each stage). Assays were taken for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) (n = 2 for each stage). Serum samples were collected via venipuncture (n = 6 for each stage) and then stored at −80 °C before the next step. A section of the ovary was obtained and stained for observation (n = 6 for each stage).
2.3. ATAC-seq Library for Frozen Pituitary Tissues
Nuclei were isolated from frozen pituitary according to the Ryan Corces protocol [26]. The frozen pituitary was placed into a pre-chilled 2 mL Dounce tissue grinder set (Cat. No. D8938-1SET) (Sigma-Aldrich, Darmstadt, Germany) containing 2 mL cold 1× HB and then thawed for 5 min. The tissue was filtered during transfer using a 70 μm cell strainer (Cat. No. 431751) (Corning, New York, NY, USA), and the homogenate was transferred to a pre-chilled Eppendorf 2 mL Lo-Bind tube (Cat. No. Z666556-250EA) (Sigma-Aldrich, Darmstadt, Germany). Pellet nuclei were generated by spinning the homogenate for 5 min at 4 °C at 350 RCF in a fixed-angle centrifuge. Iodixanol (Cat. No. D1556-250 mL) (Sigma-Aldrich, Darmstadt, Germany) was used for collecting the nuclei. Isolated nuclei were counted with trypan blue staining in a cell counter, and 30 K high-quality nuclei were used to prepare an ATAC-seq library.
This ATAC-seq library was prepared following the instructions of the TruePrep DNA Library Prep Kit V2 for Illumina (Vazyme, Nanjing, China). Nuclei pellets were incubated in a 50 µL transposition mix (10 µL 5 × TTBL buffer, 5 µL TTE Mix, 35 µL ddH2O) for 30 min at 37 °C. Transposed DNA was then purified with VAHTS DNA cleaning beads (Vazyme, Nanjing, China) and then amplified for 12~15 cycles with PCR. Libraries were purified and assessed for fragment length distribution with a Bioanalyzer Qseq 100 Bio-Fragment Analyzer (Bioptic Inc., New Taipei City, Taiwan) and submitted for paired-end 150 bp sequencing on the NextSeq platform.
2.4. ATAC-seq Data Analysis
The Fastp [27] (version 0.23.1) software was used to remove adapter sequences, poly-X in 3′ ends, and reads that had a phred-scaled quality score of less than 15 for more than 15% of the bases. The trimmed fastq files were mapped to the Anser cygnoides genome obtained from the NCBI database (AnsCyg_PRJNA183603_v1.0) using Bowtie2 [28] (version 2.3.4 with named parameters ‘--very-sensitive -X 2000′. The sequence alignment map (SAM) files were compressed to the binary alignment map (BAM) version on which Samtools [29] (version 1.7) was used to filter reads that were unmapped, mate unmapped, not primary aligned, or failed platform quality checks. Read pairs mapped to mitochondria DNA were discarded using Samtools (version 1.7). Redundancy read pairs from PCR amplification were also removed afterward using Picard [30] MarkDuplicates (version 2.18.29). The Deeptools [31] alignmentSieve (version 3.1.3) software was used to move the forward chain of the BAM file forward by 4 bp and the reverse chain backward by 5 bp with parameter ‘-ATACshift’. Then, open accessible-regions for each shifted bam file were defined by the peaks called by MACS2 [32] (version 2.1.1) with the parameters “-g ‘genome_size’ -f BAMPE –q 0.05 --keep-dup all”. Motif analysis on peak regions was performed by Hypergeometric Optimization of Motif Enrichment [33] (HOMER version 4.9.1), using the function findMotifsGenome.pl.
2.5. The Identification of DARs and Peak Annotation
Analysis of ATAC-seq for differential accessibility was carried out in the DiffBind [34] package (version 3.6.5) using the default threshold parameters of |Fold| ≥ 0.6 and p-adj ≤ 0.01 by DESeq2. Annotations of the peaks were achieved using the HORMER annotate function.
2.6. Computational Footprint Analysi
The non-redundant vertebrate TF motifs were downloaded from the JASPAR database (release 2020) [35]. Computational footprint analysis was conducted across each merged bam file using TOBIAS [36]. Briefly, the Tn5 transposase sequence preference of cutting sites was estimated and corrected using the parameter of ‘ATACorrect’. The deletion of ATAC-seq signals given rise from protein binding and the neighboring signals around binding sites were calculated using the parameter of ‘FootprintScores’. The differential binding TFs were detected using the parameter of ‘BINDetect’. All transcription factors with −log10 (p-value) above the 95% quantile or differential binding scores smaller/larger than the 5 and 95% quantiles were considered differential binding TFs.
2.7. RNA-seq and Differential Gene Expression Analysis
Total RNA extracted as described above was used for transcriptome sequencing on Illumina HiSeq 4000 to generate 150 bp paired-end reads. Clean reads filtered by Fastp software were mapped to the goose genome (AnsCyg_PRJNA183603_v1.0) using STAR 2.0.4 software with default parameters, DESeq2 [37] was performed to identifyDEGs. The threshold level was set with a standard: FDR ≤ 0.05, |log2 fold change| ≥ 1.0. Each pairwise combination of the two reproductive stages was investigated. Real-time PCR was used to validate the differential expression genes.
2.8. GO and KEGG Pathway Enrichment Analyses
The gene ontology (GO) enrichment analysis of genes was performed using Metascape (
2.9. Ovarian Histology and Serum Hormone
The geese ovaries were collected immediately after slaughter. The dissected ovarian tissues were immediately fixed in paraformaldehyde. The slides were stained by Hematoxylin and Eosin (H&E) to investigate the structure, morphology, and histological changes. Serum samples were isolated from blood samples. The concentrations of PRL and LH in serum were measured by Enzyme-Linked Immunosorbent Assay (ELISA), according to the instructions of the Prolactin ELISA kit (Cat. No. EK18156) (SAB, Baltimore, MD, USA) and LH ELISA kit (Cat. No. EK11984) (SAB, Baltimore, MD, USA). Statistical analysis was performed by the Analysis of Variance (ANOVA). A probability level of p < 0.05 (*) or p < 0.01 (**) was used to indicate significance.
3. Results
3.1. Ovarian Histology Structure and Hormone Concentration
The ovarian cortex of laying geese contains a hierarchy of all stages of developing follicles, five to eight of the largest hierarchy of preovulatory follicles are present, pre-hierarchical follicles including small yellow follicles (SYF), large white follicles (LWFs), and small white follicles (SWFs), as well as the post-ovulatory follicles (POF) are visible, whereas the ovaries manifested as regression and atresia in the brooding stage geese. In addition, ovarian histology of the laying geese and broody geese were investigated. Consistent with previous results, many primary follicles and secondary follicles can be observed in the ovaries of the laying stage; however, the brooding geese had numerous primary follicles (Supplementary Figure S1a). We measured serum PRL and LH concentrations of individuals in both stages. Compared to the laying stage, the concentration of PRL in the serum was significantly higher in the brooding stage (p < 0.01) (Supplementary Figure S1b), while a significantly lower concentration of LH was observed in the serum (p < 0.05) (Supplementary Figure S1c).
3.2. Characteristics of Chromatin Accessibility in Goose Pituitary Underlying Two Reproductive Stages
We isolated the nuclei from the geese pituitary and investigated the chromatin accessibility involved in incubation behavior. We profiled the chromatin accessibility regions in both the laying stages and the brooding stages (Figure 1). We generated more than 48 million paired-end reads for each sample. The Q20 ratio, Q30 ratio, and GC content in the clean data were calculated, which ranged from 95.2% to 98.1%, 92.2% to 94.6%, and 46.2% to 47.8%, respectively, 94.6 to 96.2% of reads from all libraries were mapped to the goose genomes from NBCI (Supplementary Table S1).
Pearson correlation analysis and PCA analysis of ATAC-seq showed the profiles of chromatin accessibility in the pituitary were strikingly divergent between laying stage and brooding stage geese (Figure 2a and Figure S2). A distribution plot of insert fragment length clearly shows that nucleosome packing and chromatin sequence insertion exhibit a distinct periodicity of approximately 200 bp in the chromatin landscape of the pituitary gland, which is the normal insertion fragment pattern for the ATAC-seq library (Figure 2b,c and Figure S3). The majority of reads were located in regions less than 100 bp, indicating that the nucleosome-free open regions were cut and sequenced. Mononucleosome, dinucleosome, and trinucleosome patterns were observed. We also assessed the quality of all libraries based on the peak signal distributions and the transcription start site enrichment (TSS) score (Supplementary Figure S4).
Open chromatin regions can generally combine specific transcription factors to affect the transcription process. We investigated the differences between transcription factors in the pituitary during the laying and brooding stages at a genome-wide level. For comparison, we first calculated the proportion of genomic positions of peaks in each library and identified 9174 and 7155 differentially accessible regions of chromatin by merging results from the laying stage and the brooding stage, respectively. (Figure 2d and Figure S5). Our data revealed that the profile of chromatin accessibility was strikingly different between the laying stage and the brooding stage. Then, we conducted motif enrichment analysis on the consensus peaks under the two states.
The motif analysis suggested that the ranking of motif usage is a little different in the pituitary between the laying stage and the brooding stage. As expected, a significant enrichment of the binding motif of CTCF, SIX2, SIX1, BORIS, RORIN, and helix-turn-helix (HTH) superclass (RFX, RFX2, RFX1, and X-box) was identified in the laying stage; likely reflecting that these transcription factors are the master regulator responsible for maintaining the physiological state of laying in pituitary (Figure 2e). However, transcriptional regulatory regions with motif binding of CTCF, BORIS, Homeobox families (SIX2, SIX1), NR family (GRE, PGR, ARE), and HTH families (RFX, RFX2) showed open chromatin accessibility in the brooding stage (Figure 2); the detailed information on all transcription factors identified is presented (Supplementary Tables S2 and S3). The observations for enrichment motif showed the epigenetic landscape of the pituitary is stable in both reproductive stages, suggesting that open chromatin regions are required for essential pathways to maintain a multitude of the endocrine function of the pituitary in geese. Thus, further analysis of the data is required to filter and identify key chromatin regions that regulate incubation behavior. In addition, we showed the trace near the TRHR gene, indicating the authenticity and repeatability of our data (Figure 2).
3.3. Differential Chromatin Accessibility Regions between Two Reproductive Stages
To identify the potential cis-regulatory elements in the pituitary at the chromatin level related to incubation behavior, we compared chromatin accessibility in the pituitary between laying and brooding. A total of 920 significant differential accessible regions (DARs) were identified between the laying and brooding pituitary using DiffBind (abs (Fold) ≥ 0.6, FDR ≤ 0.01). These DARs were annotated into 847 genes based on their distance to the nearest TSSs. The DARs were clearly classified into two distinct clusters based on the accessibility in each fraction based on the K-means clustering method (Figure 3a). The two clusters revealed strikingly divergent chromatin accessibility between the two stages. Cluster 1 (114 peaks) showed a higher average chromatin accessibility in the laying samples, while cluster 2 (n = 806) showed higher accessibility in the broody samples (Figure 3a).
Based on the genome-wide functional regions, the majority of up-regulated peaks were mapped to intronic regions, followed by intergenic regions, and promoter regions; while the down-regulated peaks mostly come from intergenic regions, intronic, and exon regions of genomes. There was no significant proportional change between up-regulated and down-regulated DARs in the intergenic and promoter-TSS regions. However, the proportion of down-regulated DARs, in both exon and TTS regions, was significantly increased, whereas the proportion of up-regulated DARs was relatively decreased. On the contrary, a higher proportion of peak up-regulation occurred in intron regions, which may suggest that enhancers are an important regulatory layer in the broody stage (Figure 3b).
We then compared the frequencies of the transcription factor binding motifs at differentially accessible regions between the two stages by using homer software and calculated the p-values and enrichment score of the motifs. Here we refer to the peaks in cluster1 as closed DARs, and the peaks in cluster2 as open DARs. Motif analysis of open DARs at the incubation behavior stage showed that the most significant transcription factor occupies sites enriched in motifs for the HTH superclass (RFX5, RFX2, RFX1, RFX and X-box), homeobox (six2, six1, Lhx3, LHX9, and EN1), NF1, and NeuroD1 in open DARs, while an enrichment of motifs for NR (ARE, GRE, and PGR), homeobox (Bapx1), EFL-1,AT3G10030,ERF115, LIN-15B, and MYB (MYB51, MYB40, and MYB39) were identified in the closed DARs at the incubation behavior stage (Figure 3c). The motif enrichment analysis in DARs reveals a set of potential TFs highly related to incubation behavior in the pituitary.
To explore the potential function of DARs, GO enrichment analysis was performed for DARs related to genes. GO annotation of biological processes related to open DARs showed highly enriched neuron projection development and neuronal differentiation in the brooding stage (Figure 3e); however, closed DARs related to genes were strongly enriched in negative regulation of transport and positive regulation of protein phosphorylation in biological processes (Figure 3).
3.4. Transcriptome Change in Pituitary Underlying Two Reproductive Stages
To assess the potential impact of chromatin changes on gene expression, transcriptome sequencing was conducted to investigate the differentially expressed genes (DEGs) in the pituitary underlying the two reproductive stages (Figure 1 and Supplementary Table S4). Principal component analysis (PCA) of the normalized count data separated the sample of the two stages (Figure 4a). We identified 225 DEGs (Fold Change ≥ 1.5, FDR ≤ 0.05) between the laying stage and the brooding stage, and 102 genes were upregulated in the brooding stage, while 123 genes were down-regulated (Supplementary Table S5). The functional genes involved in incubation behavior, such as TSHR and SREBF2 were markedly down-regulated in the brooding pituitary [15]. A bar plot showed the top 10 DEGs ranked by the FPKM value of the fold change (Figure 4b). Remarkably, the prolactin gene PRL was the most highly expressed gene in the pituitary gland and showed a clear and consistent difference between the two states (log2FoldChange = −1.53, FDR = 0.009). qPCR was performed to validate the gene expression difference (Supplementary Figure S6).
GO enrichment analysis revealed that up-regulated DEGs were associated with lipid transfer activity and transporter activity in molecular function (Figure 4c), and down-regulated DEGs were associated with monooxygenase activity, glycosyltransferase activity, and oxidoreductase activity (Figure 4d). Interestingly, the GO analysis result indicated that the up-regulation genes and down-regulation genes are linked to lipid metabolism in the biological process and play different roles in different processes of lipid metabolism. KEGG analysis revealed that the highlighted pathways are significantly overrepresented in steroid biosynthesis, metabolic pathways, and glycerophospholipid metabolism (Figure 4e). Our results indicated that the alteration of lipid metabolism in the pituitary participated in the regulation of incubation behavior.
3.5. Expression Change in DAR-Related Genes to Goose Incubation Behavior
In our study, the pituitary gland showed comparatively less variation at the level of gene expression and more change involving chromatin accessibility in goose incubation behavior. To determine the effects of chromatin accessibility on gene expression, no significant correlation was observed. We performed overlapping analysis of DARs and DEGs. Change in chromatin accessibility regulated a small set of DEGs. Five genes were both differently expressed and showed differential chromatin accessibility between the laying stage and the brooding stage (Figure 5, Supplementary Table S6). LOC106036132, OTOGL and PCDH18 showed increased expression and altered chro-matin accessibility during incubation. The expression levels of PRL and GPX3 genes were higher during the laying period, but chromatin accessibility was decreased. The ex-pression and chromatin landscape of the above five genes are shown in Supplementary Figures S7–S11.We next investigated the transcription factors within open chromatin regions of DAR-related DEGs in the goose pituitary and found SREBF2 was the unique differentially expressed transcription factor bound to DARs in the brooding stage (Supplementary Table S7). SREBF2 were significantly down-regulated and enriched in hyper-accessible regions of PRL in the broody stage (Supplementary Tables S8 and S9).
3.6. Key Transcription Factors Were Explored by Footprinting Analysis
We performed footprinting analysis to estimate the transcription factor (TF) binding on different open chromatin regions. We collected the vertebrate TF binding profiles from the JAPSAR database [35]. After correction for deviation, factors with −log10 (p-value) above the 95% quantile or differential binding scores smaller/larger than the 5 and 95% quantiles were defined as differentially bound TFs between the laying stage and the brooding stage (Figure 6). The differentially bound TFs are labeled in the Figure 6. We observed that TFs exhibit similar binding dynamics, as TFs from the same family often have highly similar binding motifs. These transcription factors may play an important regulatory role in the maintenance of incubation behavior. Interestingly, we compared the TF binding between the laying stage and the broody stage in the pituitary. Footprints of a set of transcription factors (RFX1, RFX2, RFX3, RFX5, BHLHA15, SIX1, and DUX) displayed the highest activity at the brooding stage, suggesting that factors are likely involved in incubation behavior (Figure 6).
4. Discussion
Magang geese are typical short-day breeders, which have strong incubation behavior [5]. Ninety percent of female geese present incubation behavior after laying one clutch of eggs with a high incubation constancy and a long incubation duration time. The seasonal breeding characteristics and the poor laying performance caused by incubation behavior have severely restricted the development of the waterflow industry. Over the past 20 years, avian incubation behavior has been studied extensively at the molecular level, with many previous reports focusing on identifying functional genes and genetic mutations associated with this behavior [4,40,41,42,43]. In spite of these major advances, the underlying mechanisms remain unknown as incubation behavior is a low heritability trait. Incubation behavior had been demonstrated to be induced by the combined action of three hormones [8], suggesting that changes involved in incubation behavior many occur at different levels, such as transcription and post-transcriptional regulation. The transcriptome is a dynamic component that plays a critical role in determining phenotypes. Transcriptome profiling of the HPG axis has revealed distinct expression patterns of gene and non-coding RNA associated with incubation behavior [15,16,44]. However, the genetic architecture and epigenetic regulatory mechanisms of incubation behavior is still unclear.
In the present study, we comprehensively profiled the chromatin accessibility and transcriptome in the pituitary tissue of both laying and brooding geese. Accumulating evidence has emerged showing that changes in chromatin accessibility are associated with their regulatory role in the endocrinology function of the pituitary gland [17,23,45]. Although alterations in chromatin accessibility might be related to different expression patterns between the laying stage and the brooding stage, how these alterations relate to phenotype remains unclear. Transcription factors influencing gene expression can be modulated by chromatin accessibility at transcription factor binding sites. We previously reported that the transcription factors SREBF2, PGR, and SF1 act as central signal modulators during the transition from laying to brooding at the molecular level [15]. We investigated the characteristics of open chromatin regions in the pituitary gland underlying both the laying stage and the brooding stage in terms of chromatin level; the open chromatin region showed alterations in chromatin accessibility preferentially occurring during the brooding stage. As expected, CTCF was enriched at the two reproductive stages. Our results showed that the nuclear receptor (NR) superfamily (GRE, PGR, and ARE) was significantly enriched at the brooding stage. Transcription factors of the NR superfamily were involved in a numbers of physiological processes such as reproduction, circadian rhythms, and metabolism [46]. Most members of the NR superfamily are activated by hormones, such as thyroid hormones, steroids, vitamin D, or retinoic acid [47]. Three transcription factors, glucocorticoid response elements (GRE), progesterone receptors (PGR), and androgen response element (ARE), were strongly enriched in the brooding stage. Glucocorticoid release is a classic endocrine response to stress, which can directly modulate reproductive function in the HPG axis [48]. In quail, GRE corticosterone treatments increased the expression level of GnIH through GRE, which is located in the GnIH promotor [49]. ARE is the target binding element present in promoters or enhancers of genes targeted by the androgen receptor (AR) [50]. AREs can function as glucocorticoid or progesterone response elements and vice versa as steroid hormone receptors are all similar in their receptor structures [51]. PGR belongs to steroid receptors and plays a key role in avian reproduction. Progesterone regulated chicken ovulation and incubation behavior via the pituitary [8]. PGR is a nuclear receptor transcription factor which is essential for female reproduction and plays a pleiotropic role in different tissues of the female reproductive systems [52]. In the female reproductive systems, PGR transcriptional regulation is highly diverse and tissue-specific and is related to different physiological roles in different target tissues [53]. The expression level of PGR significantly declined at the brooding stage. We suspected that PGR acted as a central signal modulator during the transition from laying to brooding via the pituitary [15].
Phenotypic changes mostly result in the activity of cis-regulatory DNA elements, driven by changes in transcription factor (TF) binding [54]. At the chromatin level, we identified 920 DARs related to the reproductive stage. The differentially open accessible regions in the brooding stage are enriched for non-promoter regions, with binding motifs of 153 transcription factors enriched. The RFX HTH family (RFX2, RFX, Rfx1, and Rfx5), the sine oculis (six) homeobox family (six2 and six1), the NeuroD bHLH family (NeuroD1 and NeuroD2), the Lhx family (LHX3 and LHX9), NF1, and EN1 were the significantly enriched in open DARs at the incubation behavior stage. The RFX HTH family are highly conserved in animals and act as master regulators of central nervous system development and ciliogenesis [55,56,57]. RFX transcription factors regulate their target genes through X-box in the promoter and are highly expressive in four organ systems: the immune system, gastrointestinal tract, reproductive system, and nervous system. RFX2 was reported as a key regulator involved in mouse spermiogenesis [58,59,60]. Through footprint analysis, similar to six1 and six2, RFX1 and RFX2 had a higher footprint score in the brooding stage. RFX5 was an essential and highly specific regulator of major histocompatibility class II (MHCII) gene expression in the immune system. Members of the six homeobox transcription factor family play an important role in organogenesis and differentiation in a wide range of animal species [61]. NeuroD transcription factors govern photoreceptor genesis and regeneration, while NeuroD knockdown prevents cell cycle exit and photoreceptor regeneration in the retina of mice [62]. NeuroD1 regulates the thyroid hormone receptor β2 and cone opsins in the developing mouse retina [63]. NeuroD1 is strongly expressed in the gonadotrope progenitor, and over-expression of NeuroD1 increases the mRNA expression level of GnRH receptor gene in mouse gonadotrope cells [64]. LIM homeodomain transcription factors are required for pituitary gland and nervous system development. In the pituitary, LHX3 is reportedly involved in the activation of the FSH β-subunit gene gonadotrope cell [65]. Lhx9 plays an important role in the regulation of cell proliferation and migration, which is essential for mouse gonad formation [66]. The role of NF-1 in the pituitary seems to be as a repressor related to human growth hormone [67].
However, the NR superfamilies (ARE, GRE, and PGR) were significantly enriched in closed DARs at the incubation behavior stage. The transcriptional activity of NRs is modulated by various ligands, including hormones and lipids [68]. By sensing changes in lipid metabolite levels, NRs drive differential gene expression producing different physiologic effects [46]. Androgens start puberty and play a role in reproduction. Pituitary androgen receptor signaling is previously reported to regulate prolactin release in males [69]. ARE motifs were present in the promoters or the enhancers of genes targeted by the androgen receptor. Previous studies showed that androgen receptors positively regulated gonadotropin-releasing hormone receptors in pituitary gonadotropes [70]. The glucocorticoid receptor is a steroid hormone-activated transcription factor that regulates the transcription of thousands of target genes by binding to GRE upstream of target genes, thus producing different physiological functions [71]. Progesterone plays critical roles in reproduction across vertebrates. Nuclear PGR is a ligand-dependent transcription factor responsible for mediating progesterone action related to reproduction [72]. Our results indicated that closed DARs related to the transcriptional activity of NRs, which sense change among androgen, glucocorticoid, and progesterone in the pituitary.
Conversely, the majority of differential gene expression between the laying and brooding stages were largely similar to those reported in previous bird studies in transcriptome level [15,16]. Cholesterol metabolic processes were significantly down-regulated in the laying stage. To the best of our knowledge, high levels of prolactin are one of the recognized characteristics of incubation behavior. However, in our data the expression pattern of PRL significantly declined in the brooding stage; notably, the serum concentration of the PRL hormone was higher in the stage of broodiness than in the laying stage, which is consistent with the observations in other breeds of geese [73]. We suspect that the brooding geese, deprived of their eggs and nests, exhibited a significant reduction in the mRNA level of PRL [45,74,75]. ATAC-seq showed that chromatin changes in intergenic regions of the PRL gene. Compared to the laying stage, steroid biosynthesis, metabolic pathways, and glycerophospholipid metabolism pathways were all significantly changed, which is consistent with previous results for birds [15]. SREBF2 is the central player of lipogenesis, as it is the master regulator of cholesterol synthesis [76]. In our study, SREBF was the differentially expressed transcription factor whose mRNA expression level was significantly down-regulated and enriched in hyper-accessible regions of PRL during the broody stage. Moreover, the lipoprotein metabolism gene is rarely upregulated during the brooding stage. Cholesteryl ester transfer protein (CETP) was the most expressed gene in the brooding stage, indicating that lipoprotein metabolism was the most dynamic change between brooding and laying.
The integration of ATAC-seq and RNA-seq results showed that chromatin changes are uncoupled from differential gene expression. We observed five gene changes in both chromatin accessibility and expression level. GPX3, which is known to be epigenetically regulated, showed significant down-regulation of mRNA expression levels, and the chromatin accessibility was altered. Loss of Gpx3 induces oxidative stress and increases prostatic intraepithelial neoplasia proliferation in prostatic cancer [77], and silencing GPX3 expression has been reported to enhance metastasis of human thyroid cancer. DAR-related DEGs, which have low expression in the brooding stage, demonstrated their function on the reproductive system. PCDH18 is thought to play a role in cell–cell connections in the brain, and haploid defects in this gene may lead to alterations in brain development and associated malformations [78]. OTOGL mutations have been found to be associated with hearing loss [79]. In addition, correlation analysis of our nine-quadrant plot shows no significant correlation between chromatin accessibility and gene expression level.
We downloaded the motif probability matrix corresponding to transcription factors from the Jaspar database for evaluating the activity of transcription factors at different stages. Four RFX transcription factors (RFX1, RFX2, RFX3, and RFX5), BHLHA15, SIX1, and DUX4 were highly active in the brooding phase, suggesting that these factors are likely involved in goose incubation behavior. The RFX transcription factor family was identified as the highly active transcription factor in the brooding stage in our footprint analysis with a significantly up-regulated chromatin accessibility. This finding is very interesting because several previous studies have shown that the RFX transcription factors play a crucial role in the ciliogenesis involved in the spermatogenesis [58,59,60].TFs binding footprints have been observed in some promoter regions. Further footprint analysis will help reveal the dynamics of open chromatin regions and chromatin accessibility of these genes and identify cis-regulatory elements associated with their specific transcription during the transition phase of incubation behavior. These potentially cis-regulatory elements provide candidates for experimental verification of protein–DNA interactions for research into the genetic breeding of poultry.
5. Conclusions
In this study, we applied ATAC-seq and RNA-seq to identify the key DARs and DEGs in the pituitary that affect the incubation behavior of geese. We first identified 920 DARs and important transcription factors that might regulate the emergence of incubation behaviors by ATAC-seq from pituitary tissues of laying geese and broody geese. We obtained genes significantly related to steroid biosynthesis in the transcriptome. In addition, the combined analysis of ATAC-seq and RNA-seq identified five genes closely related to the incubation behavior of geese; SREBF2 may bind to the hyper-accessible region of the PRL gene and affect its expression changes. Finally, footprinting analysis revealed a set of transcription factors that displayed the highest activity at the brooding stage. Taken together, we provide a new insight into the regulatory elements of goose incubation behavior from a pituitary perspective, help reveal epigenetic mechanisms regulating bird incubation behavior, and provide a new clue for genome-assisted breeding in the poultry industry.
Formal analysis, J.C.; funding acquisition, J.R., X.S. and Y.H.; investigation, D.F. and Y.X.; methodology, J.C., D.F., J.L., Y.X., X.H. and J.R.; project administration, J.X. and X.S.; resources, X.S. and Y.H.; software, J.C., J.X. and J.R.; supervision, Y.T., J.X., J.R., X.S. and Y.H.; validation, J.L.; writing—original draft, J.C., J.X. and X.S.; writing—review and editing, J.R., X.S. and Y.H. All authors have read and agreed to the published version of the manuscript.
This study was carried out in strict accordance with the Experimental Animal Committee of Zhongkai University of Agriculture and Engineering (NO. 20201226012). All efforts have been made to minimize animal suffering.
No applicable.
The data sets supporting the results of this article were included within the article and additional files. Raw sequencing reads are available at Genome Sequence Archive (GSA) database of the National Geoscience Data Centre (NGDC) (bio-project accession: PRJCA012997, ATAC-seq accession CRA008881, and RNA-seq accessions CRA008882).
All support from Xu Lab at Sun Yat-Sen University, such as technical guidance and material supply, are appreciated.
The authors declare no conflict of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Scanes, C.G.; Butler, L.D.; Kidd, M.T. Reproductive management of poultry. Animal Agriculture; Elsevier: Amsterdam, The Netherlands, 2020; pp. 349-366. [DOI: https://dx.doi.org/10.1016/b978-0-12-817052-6.00020-3]
2. Mauro, L.J.; Elde, R.P.; Youngren, O.M.; Phillips, R.E.; El Halawani, M.E. Alterations in hypothalamic vasoactive intestinal peptide-like immunoreactivity are associated with reproduction and prolactin release in the female turkey. Endocrinology; 1989; 125, pp. 1795-1804. [DOI: https://dx.doi.org/10.1210/endo-125-4-1795] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/2676472]
3. El Halawani, M.E.; Rozenboim, I. The Ontogeny and Control of Incubation Behavior in Turkeys. Poult. Sci.; 1993; 72, pp. 906-911. [DOI: https://dx.doi.org/10.3382/ps.0720906]
4. Dunn, I.C.; Miao, Y.W.; Morris, A.; Romanov, M.N.; Wilson, P.W.; Waddington, D. A study of association between genetic markers in candidate genes and reproductive traits in one generation of a commercial broiler breeder hen population. Heredity; 2004; 92, pp. 128-134. [DOI: https://dx.doi.org/10.1038/sj.hdy.6800396] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14679392]
5. Huang, Y.M.; Shi, Z.D.; Liu, Z.; Liu, Y.; Li, X.W. Endocrine regulations of reproductive seasonality, follicular development and incubation in Magang geese. Anim. Reprod. Sci.; 2008; 104, pp. 344-358. [DOI: https://dx.doi.org/10.1016/j.anireprosci.2007.02.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17368765]
6. Donham, R.S. Annual cycle of plasma luteinizing hormone and sex hormones in male and female mallards (Anas platyrhynchos). Biol Reprod.; 1979; 21, pp. 1273-1285. [DOI: https://dx.doi.org/10.1095/biolreprod21.5.1273]
7. Ramsey, S.M.; Goldsmith, A.; Silver, R. Stimulus requirements for prolactin and LH secretion in incubating ring doves. Gen. Comp. Endocrinol.; 1985; 59, pp. 246-256. [DOI: https://dx.doi.org/10.1016/0016-6480(85)90376-4]
8. El Halawani, M.E.; Silsby, J.L.; Behnke, E.J.; Fehrer, S.C. Hormonal induction of incubation behavior in ovariectomized female turkeys (Meleagris gallopavo). Biol. Reprod.; 1986; 35, pp. 59-67. [DOI: https://dx.doi.org/10.1095/biolreprod35.1.59]
9. Hoglund, A.; Henriksen, R.; Fogelholm, J.; Churcher, A.M.; Guerrero-Bosagna, C.M.; Martinez-Barrio, A.; Johnsson, M.; Jensen, P.; Wright, D. The methylation landscape and its role in domestication and gene regulation in the chicken. Nat. Ecol. Evol; 2020; 4, pp. 1713-1724. [DOI: https://dx.doi.org/10.1038/s41559-020-01310-1]
10. Bélteky, J.; Agnvall, B.; Bektic, L.; Höglund, A.; Jensen, P.; Guerrero-Bosagna, C. Epigenetics and early domestication: Differences in hypothalamic DNA methylation between red junglefowl divergently selected for high or low fear of humans. Genet. Sel. Evol.; 2018; 50, 13. [DOI: https://dx.doi.org/10.1186/s12711-018-0384-z]
11. Frankl-Vilches, C.; Gahr, M. Androgen and estrogen sensitivity of bird song: A comparative view on gene regulatory levels. J. Comp. Physiol. A Neuroethol. Sens. Neural. Behav. Physiol.; 2018; 204, pp. 113-126. [DOI: https://dx.doi.org/10.1007/s00359-017-1236-y]
12. Deplancke, B.; Alpern, D.; Gardeux, V. The Genetics of Transcription Factor DNA Binding Variation. Cell; 2016; 166, pp. 538-554. [DOI: https://dx.doi.org/10.1016/j.cell.2016.07.012]
13. Li, J.; Lee, M.; Davis, B.W.; Lamichhaney, S.; Dorshorst, B.J.; Siegel, P.B.; Andersson, L. Mutations Upstream of the TBX5 and PITX1 Transcription Factor Genes Are Associated with Feathered Legs in the Domestic Chicken. Mol. Biol. Evol.; 2020; 37, pp. 2477-2486. [DOI: https://dx.doi.org/10.1093/molbev/msaa093]
14. Castiglione, G.M.; Xu, Z.; Zhou, L.; Duh, E.J. Adaptation of the master antioxidant response connects metabolism, lifespan and feather development pathways in birds. Nat. Commun.; 2020; 11, 2476. [DOI: https://dx.doi.org/10.1038/s41467-020-16129-4]
15. Shen, X.; Bai, X.; Xu, J.; Zhou, M.; Xu, H.; Nie, Q.; Lu, X.; Zhang, X. Transcriptome sequencing reveals genetic mechanisms underlying the transition between the laying and brooding phases and gene expression changes associated with divergent reproductive phenotypes in chickens. Mol. Biol. Rep.; 2016; 43, pp. 977-989. [DOI: https://dx.doi.org/10.1007/s11033-016-4033-8]
16. Ye, P.; Ge, K.; Li, M.; Yang, L.; Jin, S.; Zhang, C.; Chen, X.; Geng, Z. Egg-laying and brooding stage-specific hormonal response and transcriptional regulation in pituitary of Muscovy duck (Cairina moschata). Poult. Sci.; 2019; 98, pp. 5287-5296. [DOI: https://dx.doi.org/10.3382/ps/pez433]
17. Ruf-Zamojski, F.; Zhang, Z.; Zamojski, M.; Smith, G.R.; Mendelev, N.; Liu, H.; Nudelman, G.; Moriwaki, M.; Pincas, H.; Castanon, R.G. et al. Single nucleus multi-omics regulatory landscape of the murine pituitary. Nat. Commun.; 2021; 12, 2677. [DOI: https://dx.doi.org/10.1038/s41467-021-22859-w]
18. Fallahshahroudi, A.; Løtvedt, P.; Bélteky, J.; Altimiras, J.; Jensen, P. Changes in pituitary gene expression may underlie multiple domesticated traits in chickens. Heredity; 2019; 122, pp. 195-204. [DOI: https://dx.doi.org/10.1038/s41437-018-0092-z]
19. Liu, H.; Wang, J.; Li, L.; Han, C.; He, H.; Xu, H. Transcriptome analysis revealed the possible regulatory pathways initiating female geese broodiness within the hypothalamic-pituitary-gonadal axis. PLoS ONE; 2018; 13, e0191213. [DOI: https://dx.doi.org/10.1371/journal.pone.0191213]
20. Lopez, M.E.; Gregory, C.C.; Porter, T. Cellular basis for elevated prolactin secretion during incubation behavior in Bantam chickens: Analysis by reverse hemolytic plaque assay. Biol. Reprod.; 1996; 54, pp. 826-833. [DOI: https://dx.doi.org/10.1095/biolreprod54.4.826]
21. Ramesh, R.; Proudman, J.A.; Kuenzel, W.J. Changes in pituitary somatotroph and lactotroph distribution in laying and incubating turkey hens. Gen. Comp. Endocrinol.; 1996; 104, pp. 67-75. [DOI: https://dx.doi.org/10.1006/gcen.1996.0142]
22. Ramesh, R.; Solow, R.; Proudman, J.A.; Kuenzel, W.J. Identification of mammosomatotrophs in the turkey hen pituitary: Increased abundance during hyperprolactinemia. Endocrinology; 1998; 139, pp. 781-786. [DOI: https://dx.doi.org/10.1210/endo.139.2.5725] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9449653]
23. Zhang, Z.; Zamojski, M.; Smith, G.R.; Willis, T.L.; Yianni, V.; Mendelev, N.; Pincas, H.; Seenarine, N.; Amper, M.A.S.; Vasoya, M. et al. Single nucleus transcriptome and chromatin accessibility of postmortem human pituitaries reveal diverse stem cell regulatory mechanisms. Cell Rep.; 2022; 38, 110467. [DOI: https://dx.doi.org/10.1016/j.celrep.2022.110467] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35263594]
24. Baskar, R.; Chen, A.F.; Favaro, P.; Reynolds, W.; Mueller, F.; Borges, L.; Jiang, S.; Park, H.S.; Kool, E.T.; Greenleaf, W.J. et al. Integrating transcription-factor abundance with chromatin accessibility in human erythroid lineage commitment. Cell Rep. Methods; 2022; 2, 100188. [DOI: https://dx.doi.org/10.1016/j.crmeth.2022.100188] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35463156]
25. Lee, E.J.; Russell, T.; Hurley, L.; Jameson, J.L. Pituitary transcription factor-1 induces transient differentiation of adult hepatic stem cells into prolactin-producing cells in vivo. Mol. Endocrinol; 2005; 19, pp. 964-971. [DOI: https://dx.doi.org/10.1210/me.2004-0034] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15637144]
26. Corces, M.R.; Trevino, A.E.; Hamilton, E.G.; Greenside, P.G.; Sinnott-Armstrong, N.A.; Vesuna, S.; Satpathy, A.T.; Rubin, A.J.; Montine, K.S.; Wu, B. et al. An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues. Nat. Methods.; 2017; 14, pp. 959-962. [DOI: https://dx.doi.org/10.1038/nmeth.4396] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28846090]
27. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics; 2018; 34, pp. i884-i890. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty560] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30423086]
28. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods; 2012; 9, pp. 357-359. [DOI: https://dx.doi.org/10.1038/nmeth.1923]
29. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Genome Project Data Processing, S. The Sequence Alignment/Map format and SAMtools. Bioinformatics; 2009; 25, pp. 2078-2079. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp352]
30. Broad Institute. Picard Toolkit; Broad Institute: Cambridge, MA, USA, 2019; Available online: https://github.com/broadinstitute/picard (accessed on 22 March 2023).
31. Ramirez, F.; Ryan, D.P.; Gruning, B.; Bhardwaj, V.; Kilpert, F.; Richter, A.S.; Heyne, S.; Dundar, F.; Manke, T. deepTools2: A next generation web server for deep-sequencing data analysis. Nucleic Acids Res.; 2016; 44, pp. W160-W165. [DOI: https://dx.doi.org/10.1093/nar/gkw257]
32. Liu, T. Use model-based Analysis of ChIP-Seq (MACS) to analyze short reads generated by sequencing protein-DNA interactions in embryonic stem cells. Methods Mol. Biol.; 2014; 1150, pp. 81-95. [DOI: https://dx.doi.org/10.1007/978-1-4939-0512-6_4]
33. Heinz, S.; Benner, C.; Spann, N.; Bertolino, E.; Lin, Y.C.; Laslo, P.; Cheng, J.X.; Murre, C.; Singh, H.; Glass, C.K. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell; 2010; 38, pp. 576-589. [DOI: https://dx.doi.org/10.1016/j.molcel.2010.05.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20513432]
34. Stark, R.; Brown, G. DiffBind: Differential Binding Analysis of ChIP-Seq Peak Data. 2012; Available online: https://bioconductor.org/packages/release/bioc/html/DiffBind.html (accessed on 22 March 2023).
35. Castro-Mondragon, J.A.; Riudavets-Puig, R.; Rauluseviciute, I.; Lemma, R.B.; Turchi, L.; Blanc-Mathieu, R.; Lucas, J.; Boddie, P.; Khan, A.; Manosalva Pérez, N. et al. JASPAR 2022: The 9th release of the open-access database of transcription factor binding profiles. Nucleic Acids Res.; 2022; 50, pp. D165-D173. [DOI: https://dx.doi.org/10.1093/nar/gkab1113] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34850907]
36. Bentsen, M.; Goymann, P.; Schultheis, H.; Klee, K.; Petrova, A.; Wiegandt, R.; Fust, A.; Preussner, J.; Kuenne, C.; Braun, T. et al. ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation. Nat. Commun.; 2020; 11, 4267. [DOI: https://dx.doi.org/10.1038/s41467-020-18035-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32848148]
37. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.; 2014; 15, 550. [DOI: https://dx.doi.org/10.1186/s13059-014-0550-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25516281]
38. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun.; 2019; 10, 1523. [DOI: https://dx.doi.org/10.1038/s41467-019-09234-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30944313]
39. Bu, D.; Luo, H.; Huo, P.; Wang, Z.; Zhang, S.; He, Z.; Wu, Y.; Zhao, L.; Liu, J.; Guo, J. et al. KOBAS-i: Intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res.; 2021; 49, pp. W317-W325. [DOI: https://dx.doi.org/10.1093/nar/gkab447]
40. Xu, H.; Zeng, H.; Luo, C.; Zhang, D.; Wang, Q.; Sun, L.; Yang, L.; Zhou, M.; Nie, Q.; Zhang, X. Genetic effects of polymorphisms in candidate genes and the QTL region on chicken age at first egg. BMC Genet.; 2011; 12, 33. [DOI: https://dx.doi.org/10.1186/1471-2156-12-33]
41. Zhou, M.; Lei, M.; Rao, Y.; Nie, Q.; Zeng, H.; Xia, M.; Liang, F.; Zhang, D.; Zhang, X. Polymorphisms of vasoactive intestinal peptide receptor-1 gene and their genetic effects on broodiness in chickens. Poult. Sci.; 2008; 87, pp. 893-903. [DOI: https://dx.doi.org/10.3382/ps.2007-00495]
42. Namken, S.; Sinpru, P.; Kamkrathok, B.; Sartsoongnoen, N.; Chaiseha, Y. Role of vasoactive intestinal peptide during the transition from incubation behavior to rearing behavior in the female native Thai chicken. Poult. Sci.; 2017; 96, pp. 3768-3774. [DOI: https://dx.doi.org/10.3382/ps/pex180]
43. Romanov, M.N.; Talbot, R.T.; Wilson, P.; Sharp, P.J. Genetic Control of Incubation Behavior in the Domestic Hen. Poult. Sci.; 2002; 81, pp. 928-931. [DOI: https://dx.doi.org/10.1093/ps/81.7.928]
44. Liu, L.; Xiao, Q.; Gilbert, E.R.; Cui, Z.; Zhao, X.; Wang, Y.; Yin, H.; Li, D.; Zhang, H.; Zhu, Q. Whole-transcriptome analysis of atrophic ovaries in broody chickens reveals regulatory pathways associated with proliferation and apoptosis. Sci. Rep.; 2018; 8, 7231. [DOI: https://dx.doi.org/10.1038/s41598-018-25103-6]
45. Liu, T.; Mi, L.; Xiong, J.; Orchard, P.; Yu, Q.; Yu, L.; Zhao, X.Y.; Meng, Z.X.; Parker, S.C.J.; Lin, J.D. et al. BAF60a deficiency uncouples chromatin accessibility and cold sensitivity from white fat browning. Nat. Commun.; 2020; 11, 2379. [DOI: https://dx.doi.org/10.1038/s41467-020-16148-1]
46. Weikum, E.R.; Liu, X.; Ortlund, E.A. The nuclear receptor superfamily: A structural perspective. Protein Sci.; 2018; 27, pp. 1876-1892. [DOI: https://dx.doi.org/10.1002/pro.3496]
47. Schubert, M.; Brunet, F.; Paris, M.; Bertrand, S.; Benoit, G.; Laudet, V. Nuclear hormone receptor signaling in amphioxus. Dev. Genes Evol.; 2008; 218, pp. 651-665. [DOI: https://dx.doi.org/10.1007/s00427-008-0251-y]
48. Whirledge, S.; Cidlowski, J.A. Glucocorticoids and Reproduction: Traffic Control on the Road to Reproduction. Trends Endocrinol. Metab.; 2017; 28, pp. 399-415. [DOI: https://dx.doi.org/10.1016/j.tem.2017.02.005]
49. Son, Y.L.; Ubuka, T.; Narihiro, M.; Fukuda, Y.; Hasunuma, I.; Yamamoto, K.; Belsham, D.D.; Tsutsui, K. Molecular basis for the activation of gonadotropin-inhibitory hormone gene transcription by corticosterone. Endocrinology; 2014; 155, pp. 1817-1826. [DOI: https://dx.doi.org/10.1210/en.2013-2076]
50. Wilson, S.; Qi, J.; Filipp, F.V. Refinement of the androgen response element based on ChIP-Seq in androgen-insensitive and androgen-responsive prostate cancer cell lines. Sci. Rep.; 2016; 6, 32611. [DOI: https://dx.doi.org/10.1038/srep32611]
51. Schauwaers, K.; De Gendt, K.; Saunders, P.T.; Atanassova, N.; Haelens, A.; Callewaert, L.; Moehren, U.; Swinnen, J.V.; Verhoeven, G.; Verrijdt, G. et al. Loss of androgen receptor binding to selective androgen response elements causes a reproductive phenotype in a knockin mouse model. Proc. Natl. Acad. Sci. USA; 2007; 104, pp. 4961-4966. [DOI: https://dx.doi.org/10.1073/pnas.0610814104]
52. Akison, L.K.; Robker, R.L. The critical roles of progesterone receptor (PGR) in ovulation, oocyte developmental competence and oviductal transport in mammalian reproduction. Reprod. Domest. Anim.; 2012; 47, (Suppl. S4), pp. 288-296. [DOI: https://dx.doi.org/10.1111/j.1439-0531.2012.02088.x]
53. Dinh, D.T.; Breen, J.; Akison, L.K.; DeMayo, F.J.; Brown, H.M.; Robker, R.L.; Russell, D.L. Tissue-specific progesterone receptor-chromatin binding and the regulation of progesterone-dependent gene expression. Sci. Rep.; 2019; 9, 11966. [DOI: https://dx.doi.org/10.1038/s41598-019-48333-8]
54. Huang, J.H.; Kwan, R.S.; Tsai, Z.T.; Lin, T.C.; Tsai, H.K. Borders of Cis-Regulatory DNA Sequences Preferentially Harbor the Divergent Transcription Factor Binding Motifs in the Human Genome. Front. Genet.; 2018; 9, 571. [DOI: https://dx.doi.org/10.3389/fgene.2018.00571] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30524473]
55. Choksi, S.P.; Lauter, G.; Swoboda, P.; Roy, S. Switching on cilia: Transcriptional networks regulating ciliogenesis. Development; 2014; 141, pp. 1427-1441. [DOI: https://dx.doi.org/10.1242/dev.074666] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24644260]
56. Sugiaman-Trapman, D.; Vitezic, M.; Jouhilahti, E.M.; Mathelier, A.; Lauter, G.; Misra, S.; Daub, C.O.; Kere, J.; Swoboda, P. Characterization of the human RFX transcription factor family by regulatory and target gene analysis. BMC Genom.; 2018; 19, 181. [DOI: https://dx.doi.org/10.1186/s12864-018-4564-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29510665]
57. Harris, H.K.; Nakayama, T.; Lai, J.; Zhao, B.; Argyrou, N.; Gubbels, C.S.; Soucy, A.; Genetti, C.A.; Suslovitch, V.; Rodan, L.H. et al. Disruption of RFX family transcription factors causes autism, attention-deficit/hyperactivity disorder, intellectual disability, and dysregulated behavior. Genet. Med.; 2021; 23, pp. 1028-1040. [DOI: https://dx.doi.org/10.1038/s41436-021-01114-z]
58. Her, Y.R.; Wang, L.; Chepelev, I.; Manterola, M.; Berkovits, B.; Cui, K.; Zhao, K.; Wolgemuth, D.J. Genome-wide chromatin occupancy of BRDT and gene expression analysis suggest transcriptional partners and specific epigenetic landscapes that regulate gene expression during spermatogenesis. Mol. Reprod. Dev.; 2021; 88, pp. 141-157. [DOI: https://dx.doi.org/10.1002/mrd.23449]
59. Wu, Y.; Hu, X.; Li, Z.; Wang, M.; Li, S.; Wang, X.; Lin, X.; Liao, S.; Zhang, Z.; Feng, X. et al. Transcription Factor RFX2 Is a Key Regulator of Mouse Spermiogenesis. Sci. Rep.; 2016; 6, 20435. [DOI: https://dx.doi.org/10.1038/srep20435]
60. Kistler, W.S.; Baas, D.; Lemeille, S.; Paschaki, M.; Seguin-Estevez, Q.; Barras, E.; Ma, W.; Duteyrat, J.L.; Morle, L.; Durand, B. et al. RFX2 Is a Major Transcriptional Regulator of Spermiogenesis. PLoS Genet.; 2015; 11, e1005368. [DOI: https://dx.doi.org/10.1371/journal.pgen.1005368]
61. Meurer, L.; Ferdman, L.; Belcher, B.; Camarata, T. The SIX Family of Transcription Factors: Common Themes Integrating Developmental and Cancer Biology. Front. Cell Dev. Biol.; 2021; 9, 707854. [DOI: https://dx.doi.org/10.3389/fcell.2021.707854]
62. Taylor, S.M.; Alvarez-Delfin, K.; Saade, C.J.; Thomas, J.L.; Thummel, R.; Fadool, J.M.; Hitchcock, P.F. The bHLH Transcription Factor NeuroD Governs Photoreceptor Genesis and Regeneration Through Delta-Notch Signaling. Investig. Ophthalmol. Vis. Sci.; 2015; 56, pp. 7496-7515. [DOI: https://dx.doi.org/10.1167/iovs.15-17616]
63. Liu, H.; Etter, P.; Hayes, S.; Jones, I.; Nelson, B.; Hartman, B.; Forrest, D.; Reh, T.A. NeuroD1 regulates expression of thyroid hormone receptor 2 and cone opsins in the developing mouse retina. J. Neurosci.; 2008; 28, pp. 749-756. [DOI: https://dx.doi.org/10.1523/JNEUROSCI.4832-07.2008]
64. Cherrington, B.D.; Bailey, J.S.; Diaz, A.L.; Mellon, P.L. NeuroD1 and Mash1 temporally regulate GnRH receptor gene expression in immortalized mouse gonadotrope cells. Mol. Cell. Endocrinol.; 2008; 295, pp. 106-114. [DOI: https://dx.doi.org/10.1016/j.mce.2008.07.017]
65. West, B.E.; Parker, G.E.; Savage, J.J.; Kiratipranon, P.; Toomey, K.S.; Beach, L.R.; Colvin, S.C.; Sloop, K.W.; Rhodes, S.J. Regulation of the follicle-stimulating hormone beta gene by the LHX3 LIM-homeodomain transcription factor. Endocrinology; 2004; 145, pp. 4866-4879. [DOI: https://dx.doi.org/10.1210/en.2004-0598]
66. Birk, O.S.; Casiano, D.E.; Wassif, C.A.; Cogliati, T.; Zhao, L.; Zhao, Y.; Grinberg, A.; Huang, S.; Kreidberg, J.A.; Parker, K.L. et al. The LIM homeobox gene Lhx9 is essential for mouse gonad formation. Nature; 2000; 403, pp. 909-913. [DOI: https://dx.doi.org/10.1038/35002622]
67. Norquay, L.D.; Jin, Y.; Surabhi, R.M.; Gietz, R.D.; Tanese, N.; Cattini, P.A. A member of the nuclear factor-1 family is involved in the pituitary repression of the human placental growth hormone genes. Biochem. J.; 2001; 354, pp. 387-395. [DOI: https://dx.doi.org/10.1042/bj3540387]
68. Xiao, Y.; Kim, M.; Lazar, M.A. Nuclear receptors and transcriptional regulation in non-alcoholic fatty liver disease. Mol. Metab.; 2021; 50, 101119. [DOI: https://dx.doi.org/10.1016/j.molmet.2020.101119]
69. O’Hara, L.; Curley, M.; Tedim Ferreira, M.; Cruickshanks, L.; Milne, L.; Smith, L.B. Pituitary androgen receptor signalling regulates prolactin but not gonadotrophins in the male mouse. PLoS ONE; 2015; 10, e0121657. [DOI: https://dx.doi.org/10.1371/journal.pone.0121657]
70. Ryan, G.E.; Bohaczuk, S.C.; Cassin, J.; Witham, E.A.; Shojaei, S.; Ho, E.V.; Thackray, V.G.; Mellon, P.L. Androgen receptor positively regulates gonadotropin-releasing hormone receptor in pituitary gonadotropes. Mol. Cell. Endocrinol.; 2021; 530, 111286. [DOI: https://dx.doi.org/10.1016/j.mce.2021.111286]
71. Frank, F.; Ortlund, E.A.; Liu, X. Structural insights into glucocorticoid receptor function. Biochem. Soc. Trans.; 2021; 49, pp. 2333-2343. [DOI: https://dx.doi.org/10.1042/BST20210419]
72. Huang, J.; Zhang, T.T.; Jiang, K.; Hong, W.S.; Chen, S.X. GFP expression pattern in pituitary and gonads under the control of nuclear progesterone receptor promoter in transgenic zebrafish. Dev. Dyn.; 2020; 249, pp. 1365-1376. [DOI: https://dx.doi.org/10.1002/dvdy.213]
73. Yao, Y.; Yang, Y.Z.; Gu, T.T.; Cao, Z.F.; Zhao, W.M.; Qin, H.R.; Xu, Q.; Chen, G.H. Comparison of the broody behavior characteristics of different breeds of geese. Poult. Sci.; 2019; 98, pp. 5226-5233. [DOI: https://dx.doi.org/10.3382/ps/pez366]
74. El Halawani, M.E.; Burke, W.H.; Dennison, P.T. Effect of nest-deprivation on serum prolactin level in nesting female turkeys. Biol. Reprod.; 1980; 23, pp. 118-123. [DOI: https://dx.doi.org/10.1095/biolreprod23.1.118] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/7417659]
75. Tong, Z.; Pitts, G.R.; Foster, D.N.; El Halawani, M.E. Transcriptional and post-transcriptional regulation of prolactin during the turkey reproductive cycle. J. Mol. Endocrinol.; 1997; 18, pp. 223-231. [DOI: https://dx.doi.org/10.1677/jme.0.0180223] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9195476]
76. Hong, X.; Roh, W.; Sullivan, R.J.; Wong, K.H.K.; Wittner, B.S.; Guo, H.; Dubash, T.D.; Sade-Feldman, M.; Wesley, B.; Horwitz, E. et al. The Lipogenic Regulator SREBP2 Induces Transferrin in Circulating Melanoma Cells and Suppresses Ferroptosis. Cancer Discov.; 2021; 11, pp. 678-695. [DOI: https://dx.doi.org/10.1158/2159-8290.CD-19-1500] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33203734]
77. Kim, U.; Kim, C.Y.; Lee, J.M.; Ryu, B.; Kim, J.; Bang, J.; Ahn, N.; Park, J.H. Loss of glutathione peroxidase 3 induces ROS and contributes to prostatic hyperplasia in Nkx3.1 knockout mice. Andrology; 2020; 8, pp. 1486-1493. [DOI: https://dx.doi.org/10.1111/andr.12828]
78. Kasnauskiene, J.; Ciuladaite, Z.; Preiksaitiene, E.; Matulevičienė, A.; Alexandrou, A.; Koumbaris, G.; Sismani, C.; Pepalytė, I.; Patsalis, P.C.; Kučinskas, V. A single gene deletion on 4q28.3: PCDH18-A new candidate gene for intellectual disability?. Eur. J. Med. Genet.; 2012; 55, pp. 274-277. [DOI: https://dx.doi.org/10.1016/j.ejmg.2012.02.010]
79. Yariz, K.O.; Duman, D.; Zazo Seco, C.; Dallman, J.; Huang, M.; Peters, T.A.; Sirmaci, A.; Lu, N.; Schraders, M.; Skromne, I. et al. Mutations in OTOGL, encoding the inner ear protein otogelin-like, cause moderate sensorineural hearing loss. Am. J. Hum. Genet.; 2012; 91, pp. 872-882. [DOI: https://dx.doi.org/10.1016/j.ajhg.2012.09.011]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
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.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
1 Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan 430070, China; Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
2 College of Animal Science & Technology, Zhongkai University of Agriculture and Engineering, Guangzhou 510225, China; Guangdong Province Key Laboratory of Waterfowl Healthy Breeding, Guangzhou 510225, China
3 College of Life Science, Sun Yat-Sen University, Guangzhou 510642, China