Transcriptome analysis revealed the possible regulatory pathways initiating female geese broodiness within the hypothalamic-pituitary-gonadal axis

<div><p>Geese have the strongest tendency toward broodiness among all poultry. The mechanisms initiating broodiness within the goose hypothalamic-pituitary-gonadal axis (HPGA) are still unclear. Here, we reported the transcriptome differences between laying and initial nesting within the HPGA tissues of geese. We constructed a unigene database based on HPGA tissues and identified 128,148 unigenes, 100% of which have been annotated. By using Digital Gene Expression (DGE) sequencing, we screened 19, 110, 289, and 211 differentially expressed genes (DEGs) in the hypothalamus, pituitary gland, stroma ovarii, and follicles, respectively, between laying and nesting geese. Expression changes of hypocretin (HCRT) and pro-opiomelanocortin (POMC) in the hypothalamus of nesting geese may cause appetite reduction, which is possibly the first step and a prerequisite to initiate broodiness. In addition to prolactin (PRL), follicle-stimulating hormone (FSH) and luteinizing hormone (LH), genes including oxytocin-neurophysin (OXT), chordin-like protein 1 (CHRDL1) and growth hormone (GH), expressed in the pituitary gland, are new candidate molecules that may be involved in broodiness in geese. Heme oxygenase 1 (HMOX1) in the pituitary gland, the proto-oncogene c-Fos (FOS), heat shock protein 90-alpha (HSP90AA), and cyclin-dependent kinase 1 (CDK1) in the ovary that may consolidate and transduce signals regulating the HPGA during broodiness in geese.</p></div

The DEGs of interest that enriched in HPGA-related tissues.

<p>The DEGs of interest that enriched in HPGA-related tissues.</p

Volcano map of differential gene expression.

<p>The horizontal coordinates represent the fold changes of genes between the laying and initial nesting groups. The longitudinal coordinates represent the statistical significance of the changes in gene expression. The smaller the p value, the greater the value of -log₁₀p. Each dot in the image represents one gene; the blue ones indicate no significant difference in the gene, while the red ones indicate a significant difference. All the DEGs enriched in the hypothalamus, pituitary gland, stroma ovarii and walls of follicles with diameters of 8–10 mm between laying geese and early nesting geese are provided in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s008" target="_blank">S3</a>–<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s011" target="_blank">S6</a> Tables.</p

Hierarchical clusters of DEGs.

<p>A, Clustering of the DEGs listed in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.t004" target="_blank">Table 4</a>. B, Clustering of groups based on gene expression patterns among groups.</p

Quality analyses of DEG Library.

<p>Quality analyses of DEG Library.</p

Protein-protein interaction relationships based on the DEGs.

<p>A, <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.t004" target="_blank">Table 4</a> lists the DEGs that were analyzed. B, Networks are based on the DEGs that code for secreted proteins in the hypothalamus (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.t004" target="_blank">Table 4</a>) and all DEGs from the pituitary gland (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s009" target="_blank">S4 Table</a>). C, Networks based on the DEGs that code for secreted proteins in the hypothalamus and pituitary gland (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.t004" target="_blank">Table 4</a>) and all DEGs in the stroma ovarii and follicle wall (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s009" target="_blank">S4</a> and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s011" target="_blank">S6</a> Tables). The information on the interactions between the proteins was downloaded from the STRING database (<a href="http://string-db.org/" target="_blank">http://string-db.org/</a>).</p

Size distribution of valid unigenes compared with the public databases.

<p>Size distribution of valid unigenes compared with the public databases.</p

The entire experimental design and inclusion criteria for the groups of geese.

<p>A, The entire experimental design and sequencing procedure. L_Hy, L_Pitu, L_Ov and L_ThB, respectively, mean hypothalamus, pituitary gland, stroma ovarii and walls of follicles with a diameter range of 8–10 mm from laying geese. N_Hy, N_Pitu, N_Ov and N_ThB refer to the corresponding tissues from geese at the beginning of nesting. B, The standards for dividing geese into initial laying and nesting basing on the status of hierarchical follicles status (before the anatomical observation, the geese were distinguished first by their nesting behaviors, such as nesting, foraging times and fluffy feathers).</p

Regulatory processes within the HPGA that may take part in initiating broodiness in geese.

<p>Regulatory processes within the HPGA that may take part in initiating broodiness in geese.</p

Validation of the sequencing data by qRT-PCR.

<p>A-D, Validation of results for the hypothalamus, the pituitary gland, the stroma ovarii and the walls of follicles (diameter range 8–10 mm), respectively. E, The correlation coefficient between qRT-PCR and DEG data was analyzed based on all validated genes. The genes were selected randomly from <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s008" target="_blank">S3</a>–<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s011" target="_blank">S6</a> Tables, and the primers designed for qRT-PCR are listed in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0191213#pone.0191213.s006" target="_blank">S1 Table</a>. For both the qRT-PCR and DEG results, the fold changes were calculated based on the expression levels in laying and nesting group and were then converted to log₂(Fold change).</p