# Transcriptional profiling and pathway analysis reveal differences in pituitary gland function, morphology, and vascularization in chickens genetically selected for high or low body weight.

Laura E Ellestad, Larry A Cogburn, Jean Simon, Elisabeth Le Bihan-Duval, Samuel E Aggrey, Mardi S Byerly, Michel Jacques M.J. Duclos, Tom E Porter

► To cite this version:

Laura E Ellestad, Larry A Cogburn, Jean Simon, Elisabeth Le Bihan-Duval, Samuel E Aggrey, et al. Transcriptional profiling and pathway analysis reveal differences in pituitary gland function, morphology, and vascularization in chickens genetically selected for high or low body weight. BMC Genomics, 2019, 20, pp.316. 10.1186/s12864-019-5670-9. hal-02627980

HAL Id: hal-02627980

https://hal.inrae.fr/hal-02627980

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Transcriptional profiling and pathway analysis reveal differences in pituitary gland function, morphology, and vascularization in chickens genetically selected for high or low body weight

Laura E. Ellestad¹,², Larry A. Cogburn³, Jean Simon⁴, Elisabeth Le Bihan-Duval⁴, Samuel E. Aggrey¹, Mardi S. Byerly², Michel J. Duclos⁴ and Tom E. Porter²*

## Abstract

**Background:** Though intensive genetic selection has led to extraordinary advances in growth rate and feed efficiency in production of meat-type chickens, endocrine processes controlling these traits are still poorly understood. The anterior pituitary gland is a central component of the neuroendocrine system and plays a key role in regulating important physiological processes that directly impact broiler production efficiency, though how differences in pituitary gland function contribute to various growth and body composition phenotypes is not fully understood.

**Results:** Global anterior pituitary gene expression was evaluated on post-hatch weeks 1, 3, 5, and 7 in male broiler chickens selected for high (HG) or low (LG) growth. Differentially expressed genes (DEGs) were analyzed with gene ontology categorization, self-organizing maps, gene interaction network determination, and upstream regulator identification to uncover novel pituitary genes and pathways contributing to differences in growth and body composition. A total of 263 genes were differentially expressed between HG and LG anterior pituitary glands ($P \le 0.05$ for genetic line-by-age interaction or main effect of line; $\ge 1.6$-fold difference between lines), including genes encoding four anterior pituitary hormones. Genes involved in signal transduction, transcriptional regulation, and vesicle-mediated transport were differentially expressed and are predicted to influence expression and secretion of pituitary hormones. DEGs involved in immune regulation provide evidence that inflammation and response to cellular stressors may compromise pituitary function in LG birds, affecting their ability to adequately produce pituitary hormones. Many DEGs were also predicted to function in processes that regulate organ morphology and angiogenesis, suggesting pituitary gland structure differs between the divergently selected lines.

**Conclusions:** The large number of DEGs within the anterior pituitary gland of birds selected for high or low body weight highlights the importance of this gland in regulating economically important traits such as growth and body composition in broiler chickens. Intracellular signaling, transcriptional regulation, and membrane trafficking are important cellular processes contributing to proper hormone production and secretion. The data also suggest that pituitary function is intimately tied to structure, and organization of the gland could influence hypothalamic and systemic metabolic inputs and delivery of hormones regulating growth and metabolism into peripheral circulation.

**Keywords:** Adenohypophysis, Growth, Metabolism, Body composition, Abdominal fatness, Pituitary hormones, Secretion, Angiogenesis, Cellular compromise, Microarray analysis

<footnote>* Correspondence: teporter@umd.edu
²Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA
Full list of author information is available at the end of the article</footnote>
## Background

The hypothalamus and pituitary gland are central components of the neuroendocrine system, which integrates internal and external cues to regulate important physiological processes including growth, metabolism, response to stress, and reproduction. At the central level, the hypothalamus transmits central nervous system signals in the form of releasing factors and release-inhibiting factors to the anterior pituitary gland, which relays information to target endocrine organs through secretion of trophic hormones [1]. Three of the six anterior pituitary hormones, adrenocorticotropic hormone [ACTH; derived from proteolytic cleavage of pro-opiomelanocortin (POMC)], thyroid-stimulating hormone (TSH), and growth hormone (GH), are important regulators of growth, metabolism, and body composition in all vertebrates.

Primarily through genetic improvements aimed at enhancing growth and nutrient utilization, broiler (meat-type) chickens have emerged as a fast-growing and efficient high-quality protein source for human consumption. Despite the importance of understanding endocrine control of growth and metabolism in chickens as it relates to optimizing efficiency of feed utilization, these processes are still not clearly understood [2, 3]. Two valuable experimental genetic model systems of selection for high and low body weight are the high weight selected (HWS) and low weight (LWS) selected lines generated and maintained at Virginia Tech [4, 5], and the high growth (HG) and low growth (LG) lines of chickens developed and maintained at the Institut National de la Recherche Agronomique (INRA) [6]. The HWS and LWS lines have been continually selected for over 50 generations and now differ almost 10-fold in body weight at the age of selection (8 weeks) [4, 5, 7]. They have been widely used to investigate growth, metabolism, and feed intake in broiler chickens. Alterations in expression of components of the somatotropic axis, including minor differences in pituitary GH expression [8], and muscle regulatory genes [9] have been observed between HWS and LWS birds. Metabolic differences between these genetic groups are highlighted by alterations in pancreatic function, glucose homeostasis, and metabolic flux in skeletal muscle and adipose tissue [10–12]. The HG and LG lines of broiler chickens have been divergently selected for a large difference in body weight at two developmental stages (juveniles at 8 weeks and adults at 36 weeks) [6], and they differ markedly in their growth curves [13] and body composition, exhibiting more than a 3-fold difference in body weight [14, 15] and an almost 20-fold difference in abdominal fatness [15] by 6 to 7 weeks post-hatch. Heavier HG birds have increased muscle fiber number and size relative to LG birds [16, 17], and their muscle cells are more responsive to in vitro stimulation with insulin-like growth factor

(IGF) 1 [18]. Despite this, muscle levels of IGF receptor do not appear to differ between the two lines [19]. Hepatic IGF1 and IGF2 mRNA, as well as circulating IGF1, IGF2, and insulin, are higher in HG birds between 1 and 6 weeks post-hatch [14]. Transcriptional profiling of abdominal fat in the two lines suggests that HG birds express higher levels of transcription factors linked to adipogenesis, while LG birds exhibit increased expression of transcripts for genes that promote energy expenditure and are involved in hemostasis [15]. In addition to these molecular and endocrine differences, several quantitative trait loci for growth [20], metabolism and body composition [21], and breast meat quality [22] have been identified through use of the HG and LG lines. The large difference in body weight and ancillary changes in body composition due to divergent selection make the HG and LG chickens an ideal model to investigate the genetic basis for extreme differences in growth and metabolism, two significant economic traits in this important agricultural species. Presently, we have used this population for the discovery of additional neuroendocrine factors that regulate growth and body composition.

Consumer demand for food produced in the absence of antibiotics continues to increase, providing particular challenges for poultry production systems in terms of maintaining bird health and growth efficiency [23]. Therefore, it is critical to understand fundamental biological processes that regulate growth and body composition of broiler chickens for successful development of novel alternatives to the use of antibiotic growth promoters. Uncovering genetic mechanisms that govern these processes in the anterior pituitary gland will provide much-needed information on factors regulating growth, feed efficiency, and nutrient utilization. In the present study, we used our Del-Mar 14 K Chicken Integrated Systems Microarray [24, 25] for transcriptional analysis of the anterior pituitary gland in HG and LG broiler chicken lines at 1, 3, 5, and 7 weeks post-hatch. Novel genes and gene interaction networks were identified that regulate expression and secretion of pituitary hormones and the structure and function of the anterior pituitary gland.

## Results

### Phenotypic characterization of HG and LG chickens

In order to assess phenotypic characteristics of birds used in this study (previously reported in [15]), body weights were determined and blood was collected for hormone analysis at time of tissue collection on post-hatch weeks 1, 3, 5, and 7. As shown in Fig. 1, body weight and percentage abdominal fat were significantly different between lines at most ages ($P \le 0.05$; $n = 8$). Divergence in body weight and abdominal fatness was most evident at week 7, where a 3.2-fold difference in body weight and a 19.6-fold difference in abdominal fat
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Phenotypic characterization of juvenile male HG and LG chickens. **a** Body weight (BW, kg; top) and **b** abdominal fat (% BW; bottom) of birds from each line at post-hatch weeks 1, 3, 5, and 7 (n = 8 birds per line). Values are expressed as mean ± SEM. When the two-way ANOVA demonstrated a significant line-by-age interaction (P ≤ 0.05), the presence of an asterisk (*) denotes a significant difference between the lines at the indicated age (P ≤ 0.05). Data included in this figure have been previously published in (21)</figcaption></figure>

were observed. Thyroid hormones are primary regulators of longitudinal growth in birds. Plasma levels of the pro-hormone thyroxine (T₄) exhibited a line-by-age interaction and were higher in HG birds than LG birds on week 7 (Table 1). Plasma levels of metabolically active 3,5,3'-triodothyronine (T₃) showed significant main effects of genetic line and age. Levels were elevated in the HG birds and decreased between weeks 1 and 7 in both lines (Table 1). Thus, LG chickens could be considered to be slightly hypothyroid when compared to the HG chickens.

### Expression profiles for pituitary hormones in HG and LG birds

The primary function of the anterior pituitary gland is synthesis and secretion of five major hormones: ACTH,

TSH, GH, prolactin (PRL), follicle-stimulating hormone (FSH), and luteinizing hormone (LH). TSH, FSH, and LH are heterodimers consisting of a common α-glycoprotein hormone subunit (CGA) and distinct hormone-specific β-subunits. Differences in pituitary hormone expression between HG and LG birds were determined by reverse transcription-quantitative PCR (RT-qPCR) analysis (Fig. 2). There was a significant increase in expression levels of POMC, TSHβ and LHβ mRNA between post-hatch weeks 1 and 7, and a main effect of line where HG birds had higher levels of TSHβ, PRL, FSHβ, LHβ, and CGA mRNA than LG birds (P ≤ 0.05; n = 4). There were no significant differences in pituitary GH expression over time or between lines as determined by RT-qPCR (P > 0.05; n = 4),

<table><caption>**Table 1** Plasma thyroxine (T₄) and 3,5,3'-triodothyronine (T₃) levels in HG and LG chickens</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Week 1</th><th rowspan="2">Week 3</th><th rowspan="2">Week 5</th><th rowspan="2">Week 7</th><th colspan="3">P-values</th></tr><tr><th>Line x Age</th><th>Line</th><th>Age</th></tr></thead><tbody><tr><td><sup>a</sup>T<sub>4</sub> (ng/ml)</td><td></td><td></td><td></td><td></td><td>0.0268</td><td>0.1154</td><td>0.0151</td></tr><tr><td>HG</td><td>24.9 ± 0.9</td><td>28.7 ± 1.2</td><td>27.9 ± 1.9</td><td>32.1 ± 0.5*</td><td></td><td></td><td></td></tr><tr><td>LG</td><td>23.2 ± 1.4</td><td>29.9 ± 1.1</td><td>29.0 ± 2.2</td><td>24.2 ± 2.3*</td><td></td><td></td><td></td></tr><tr><td><sup>a</sup>T<sub>3</sub> (ng/ml)</td><td></td><td></td><td></td><td></td><td>0.4427</td><td>0.0035</td><td>0.0027</td></tr><tr><td>HG</td><td>3.0 ± 0.2</td><td>2.7 ± 0.3</td><td>2.1 ± 0.1</td><td>2.2 ± 0.1</td><td></td><td></td><td></td></tr><tr><td>LG</td><td>2.3 ± 0.2</td><td>2.2 ± 0.1</td><td>1.9 ± 0.2</td><td>1.7 ± 0.2</td><td></td><td></td><td></td></tr></tbody></table>

<footnote><sup>a</sup>Values are mean ± SEM for n = 4 birds per line</footnote>

<footnote><sup>b</sup>When the two-way ANOVA demonstrated a significant line-by-age interaction (P ≤ 0.05), the presence of an asterisk (*) denotes a significant difference between the lines at the indicated age (P ≤ 0.05)</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Gene expression profiles for pituitary hormones in HG and LG chickens. Levels of mRNA for POMC (a), TSHβ (b), GH (c), PRL (d), FSHβ (e), LHβ (f), and CGA (g) were determined on post-hatch weeks 1, 3, 5, and 7 and normalized to levels of GAPDH mRNA. Values (mean + SEM; n = 4 birds per line) are expressed relative to the age and line with the highest expression level (set to 100% prior to log₂-transformation for statistical analysis). The line-by-age interaction was not significant for any hormone (P ≤ 0.05). Significant main effects of line and/or age (P ≤ 0.05) are indicated at the top of each graph, when present</figcaption></figure>

although birds in the LG line tended to have higher levels of pituitary GH mRNA on weeks 3, 5, and 7. These results indicate that normal synthesis, and potentially secretion, of a majority of the pituitary hormones are substantially compromised in LG birds.

### Transcriptional profiling of pituitary mRNA expression

In order to investigate differences in expression which could be contributing to altered pituitary hormone production and, ultimately, divergent growth and body composition, global gene expression profiling using the Del-Mar 14 K Chicken Integrated Systems Microarray was conducted in anterior pituitary glands of HG and LG broilers at weeks 1, 3, 5, and 7 post-hatch. A summary of the microarray data is given in Table 2. The array contains 19,200 spots representing 14,053 unique genes, over half of the known chicken genome. There were a total of 2312 cDNA probes that were significantly different across all

<table><caption>**Table 2** Summary of microarray analysis</caption><thead><tr><th>Category</th><th>Number</th></tr></thead><tbody><tr><td>Spots printed on array</td><td>19,200</td></tr><tr><td>Spots with cDNA probes</td><td>17,834</td></tr><tr><td>Unique genes represented</td><td>14,053</td></tr><tr><td><sup>a</sup>cDNA probes within trimmed dataset</td><td>10,437</td></tr><tr><td><sup>b</sup>Significant cDNA probes</td><td>2310</td></tr><tr><td>Line-by-age interaction</td><td>353</td></tr><tr><td>Main effect of line</td><td>1100</td></tr><tr><td>Main effect of age</td><td>1058</td></tr><tr><td><sup>c</sup>Differentially expressed genes</td><td>291</td></tr><tr><td>Line-by-age Interaction</td><td>110</td></tr><tr><td>Main effect of line</td><td>181</td></tr><tr><td>Main effect of age</td><td>72</td></tr><tr><td><sup>d</sup>Genes annotated with GeneBase</td><td>260</td></tr></tbody></table>

<footnote><sup>a</sup>cDNA probes remaining after quality control check for removal of any that were undetectable, below background, malformed, had saturated pixel intensity, or were not present on at least half of the microarray slides</footnote>

<footnote><sup>b</sup>cDNA probes with detectable background-corrected signal on at least 16 microarray slides and significantly different at P ≤ 0.05</footnote>

<footnote><sup>c</sup>cDNA probes with detectable background-corrected signal on at least 16 microarray slides, P ≤ 0.05 for line-by-age interaction or main effect of line, and a difference between groups of 1.6-fold or greater (0.68-fold on a log₂-scale)</footnote>

<footnote><sup>d</sup>Number of 291 differentially expressed genes successfully annotated with protein identifiers using the GeneBase tool at http://cogburn.dbi.udel.edu/</footnote>
comparisons. Of these, the genes of most interest were identified by the 1631 cDNA probes that were significantly different (P ≤ 0.05; n = 4) for the main effect of line (1100 probes) or the line-by-age interaction (353 probes). These effects indicate that expression is consistently higher or lower in one of the lines (main effect of line) or developmental expression profiles differ by line (line-by-age interaction). In our dataset, differentially expressed genes (DEGs) were defined as those that were detected in at least half the samples, had a statistically significant line or line-by-age effect (P ≤ 0.05), and a 1.6-fold difference between lines. Based on these criteria, there were 291 DEGs between the HG and LG pituitaries, and these represented 263 unique DEGs due to duplication of cDNA probes on the array. Of note, this dataset of DEGs included four pituitary hormones, namely GH, TSHβ, LHβ, and FSHβ (Additional file 1). To facilitate downstream functional analysis, these DEGs were annotated using the GeneBase tool on our website (http://cogburn.dbi.udel.edu/), which provides protein identifiers for cDNA microarray probes. Of the 291 DEGs submitted, 260 were successfully annotated with protein identifiers.

DEGs were organized into gene ontology (GO) categories in order to identify cellular biological processes and molecular functions that could contribute to differences in growth and body composition between the genetically divergent lines (Fig. 3). For biological process, the largest number of genes were included in categories involved in signal transduction (104 total DEGs), transport (148 total DEGs), anatomical structure and development (94 total DEGs), and response to stress/immune system processes (140 total DEGs). Consistent with this, one of the top molecular functions identified was signal transducer activity (27 total DEGs), and a large number of genes were also categorized as playing a role in DNA binding or transcription (60 total DEGs). These results indicate that processes which are related to hormone expression and secretion, gland structure and morphology, and the immune response could lead to differences in pituitary function between the HG and LG lines.

To identify alterations in pituitary gene expression between the lines that could affect growth and metabolism, DEGs were also subject to self-organizing maps (SOMs) analysis. For each of the 291 DEGs, relative expression data (log₂ratio<sub>HGmean</sub> – log₂ratio<sub>LGmean</sub>) for each age were organized into 16 clusters (c0 – c15) in a 4 × 4 configuration (Fig. 4, Additional file 2). Genes contributing to differences in phenotype should be contained in clusters where relative expression between the lines changed in a manner consistent with divergence in growth and adiposity (i.e. primarily after week 1), namely c1 (21 DEGs), c2 (22 DEGs), c4 (29 DEGs), c5 (26 DEGs), c6 (8 DEGs), c8 (22 DEGs), c9 (17 DEGs), c14 (23 DEGs), and c15 (15 DEGs). The relative expression patterns

contained within these clusters changed substantially between weeks 1 and 3, the time period when body weight diverged, and remained different for the remainder of the ages examined. Genes for pituitary hormones were contained within c4 (*GH*), c14 (*TSHβ* and *LHβ*), and c15 (*FSHβ*), reflecting an increase in *GH* expression in LG birds over time and the increase in *TSHβ*, *LHβ*, and *FSHβ* at later ages in HG birds. Also contained within

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 3** GO analysis of DEGs in the anterior pituitary gland of HG and LG chickens. Genes with a significant line-by-age interaction or a significant main effect of line (P < 0.05; n = 4 birds per line) as determined by two-way ANOVA and a difference in expression of at least 1.6-fold between HG and LG lines (291 total) were organized into GO categories based on biological process (a) and molecular function (b), using the GORetriever and GOSlimViewer tools on AgBase (http://agbase.arizona.edu/). Shown are results for the top 10 groups within each annotation category</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Self-organizing maps (SOMs) analysis of DEGs in the anterior pituitary of HG and LG chickens. Genes with a significant line-by-age interaction or a significant main effect of line ($P \le 0.05$; $n = 4$ birds per line) as determined by two-way ANOVA and a difference in expression of at least 1.6-fold between HG and LG lines (291 total) were organized into 16 clusters (c0-c15) in a 4 × 4 configuration. Data were organized into clusters as mean log₂-ratio of HG – mean log₂-ratio of LG. The number of genes contained within each cluster is indicated in parentheses at the top center of each cluster. Examples of important genes upregulated in LG chickens (green) or HG chickens (red) are identified in c4 (GH), c14 (TSHβ and LHβ) and c15 (FSHβ)</figcaption></figure>

these and closely-related clusters are genes such as the following: secreted factors semaphorin 3F (SEMA3F) (c9) and semaphorin 3C (SEMA3C) (c15); receptors serotonin receptor 1B (5HTR1B) and glutamate receptor, ionotropic, AMPA 1 (GRIA1) (c8), frizzled-4 (FZD4) and fms-related tyrosine kinase 1 (FLT1) (c14), and plexin D1 (PLXND1) (c15); intracellular adenylate cyclase 2 (ADCY2) and dexamethasone-induced ras-related protein 1 (RASD1) (c4), chromogranin B (CHGB) (c9), and mitogen-activated protein kinase kinase kinase 7 (MAP 3 K7) (c14); and transcription factors UNC homeobox (UNCX) (c4), nuclear receptor subfamily 0, group B, member 1 (NR0B1) (c8), and mothers against decapentaplegic homolog 1 (SMAD1) (c14).

To confirm microarray expression patterns and selected relevant SOMs clusters, RT-qPCR was performed to measure mRNA expression for genes in five different clusters, as well as one that showed only a significant age effect (Fig. 5). Expression patterns of mRNA for avian leucosis virus envelope (ALVE) (c4), 5HTR1B (c8), CHGB (c9), TSHβ (c14), SEMA3C (c15), and stathmin 1 (STMN1) (no cluster; age effect only) are comparable when determined by microarray (top graph in each panel) or RT-qPCR (bottom graph in each panel), confirming validity of the genome-wide transcriptional

profiling and organization into SOMs clusters. Differences in expression between secreted factors and receptors, intracellular signaling molecules and chaperones, and transcription factors could affect pituitary function between HG and LG birds by contributing to differences in expression and secretion of pituitary hormones.

## Identification of gene interaction networks

To identify important pathways and gene networks associated with differences in pituitary mRNA expression between birds with different genetic growth potential, the functionally annotated gene list was submitted to Ingenuity Pathway Analysis (IPA) software for further functional annotation and analysis. Of the initial 260 DEGs, 249 had identifiers that mapped to the IPA annotated database and were eligible for incorporation into gene interaction networks and use in biological function identification. The top five gene interaction networks identified by IPA are listed in Table 3. The top network was associated with “Molecular Transport” (Fig. 6), and differential expression patterns of several genes within each network were confirmed by RT-qPCR (Fig. 7). A second network was identified as being associated with “Organ Morphology and Cellular Compromise” (Fig. 8), and RT-qPCR was used to confirm expression patterns
<figure><img src="image_8.png" /><figcaption>**Fig. 5** RT-qPCR confirmation of expression profiles determined by microarray and SOMs analysis for select DEGs. Microarray (a-f, top graph) and RT-qPCR (a-f, bottom graph) data are shown for each gene on post-hatch weeks 1, 3, 5, and 7. For RT-qPCR data, levels of mRNA for each gene were normalized to levels of GAPDH mRNA. Values (mean + SEM; n = 4 birds per line) are expressed relative to the age and line with the highest expression level for each technique (set to 100% prior to log₂-transformation for statistical analysis). The RT-qPCR data for TSHβ mRNA (d, bottom) are also shown in Fig. 2b. When the two-way ANOVA demonstrated a significant line-by-age interaction (P ≤ 0.05), the presence of an asterisk (*) denotes a significant difference between the lines at the indicated age (P ≤ 0.05). Significant main effects of line and/or age (P ≤ 0.05) are indicated at the top of each graph when the line-by-age interaction was not significant (P > 0.05).</figcaption></figure>

of genes in this network (Fig. 9). Of particular interest are DEGs whose relative expression switches between weeks 1 and 3 or those whose magnitude increases over time, as they likely drive differences in phenotype.

The top network included genes involved in molecular transport (Fig. 6). This network contained several genes

that play a role in signal transduction, including secreted ligands [vascular endothelial growth factor (VEGF) and SEMA3F], receptors [FLT1 and protein tyrosine phosphatase, receptor type Z1 (PTPRZ1)], and intracellular signaling molecules [RASD1, ADCY2, hematopoietic cell kinase (HCK), MAP3K7, mitogen-activated protein kinase/
<table><caption>**Table 3** Top gene interaction networks containing DEGs within pituitary glands of HG and LG chickens</caption><thead><tr><th><sup>a</sup>Associated Network Functions</th><th><sup>b</sup> DEGs (#)</th><th><sup>c</sup>Total Genes (#)</th></tr></thead><tbody><tr><td>Molecular transport, inflammatory disease</td><td>30</td><td>34</td></tr><tr><td>Cancer, developmental disorder, lipid metabolism</td><td>29</td><td>33</td></tr><tr><td>Organ morphology, cellular compromise, cardiovascular system development and function</td><td>21</td><td>34</td></tr><tr><td>Endocrine system development and function, carbohydrate metabolism, small molecule biochemistry</td><td>16</td><td>34</td></tr><tr><td>Cell-to-cell signaling and interaction, cellular function and maintenance, tumor morphology</td><td>16</td><td>34</td></tr></tbody></table>

extracellular signal-regulated kinase kinase 1 (MEK1), and extracellular signal-regulated kinase 1/2 (ERK1/2) and regulatory subunits [protein phosphatase 3 regulatory subunit B alpha (PPP3R1) and A-kinase anchoring protein 7 (AKAP7)]. Also within this network and the dataset as a whole are genes for several calcium binding proteins involved in intracellular vesicular trafficking and protein synthesis quality control [annexin A2 (ANXA2), copine 1 (CPNE1), calreticulin (CALR), and heat shock protein family A, member 5 (HSPA5)], as well as membrane trafficking [ankyrin repeat domain-containing protein 27 (ANKRD27) and microtubule associated protein RP/EB family member 2 (MAPRE2). Three of the pituitary hormones (GH, LH, and FSH) in our dataset were also included, and differences in their expression levels between lines either switches or increases in intensity over time. This indicates that alterations in cellular signaling, membrane trafficking components, and vesicular secretion may contribute to differences in hormone expression and secretion between HG and LG birds. Interestingly, the network was also populated with genes with roles in angiogenesis [(FLT1, VEGF, HCK1, podocalyxin-like protein 1 (PODXL), and hemoglobin alpha 1/hemoglobin alpha 2 (HBA1/HBA2)], suggesting that alterations in microvasculature formation and support differ between HG and LG birds.

A second network of interest contained genes with a role in cellular compromise (Fig. 7). Within this network are genes associated with antigen presentation [beta-2-microglobulin (B2M), major histocompatibility complex (MHC), human leukocyte antigen A (HLA-A), and T-cell receptor (TCR)], antiviral response [2'-5'-oligoadenylate synthetase 3 (OAS3)], and cytokine signaling or transcriptional regulation [interleukin-1, interferon alpha, inhibitor of kappa B, p38 mitogen-activated protein kinase, and signal transducer and activator of transcription 1 (STAT1)]. Genes which may be responsive to other cellular stressors were also included, such as those involved in apoptosis [caspase 8 (CASP8), cyclin E, protein kinase B/Akt, and retinoblastoma protein] or cellular responses to DNA damage [DNA damage-inducible transcript 4 (DDIT4)]. Comparing this network on post-hatch week 1 with those

of later ages, it is clear by the increase in numbers and intensity of green-colored DEGs that after divergence in growth, LG birds have increased expression of factors associated with the immune response, inflammation, and other cellular stressors. Altogether, this indicates that gene expression patterns related to the immune response and other compromising events may contribute to differences in pituitary function.

The top biological functions associated with DEGs in the anterior pituitary of HG and LG birds are summarized in Table 4; DEGs contained within each category are listed in Additional file 3. In terms of molecular and cellular function, 65 DEGs were associated with cellular growth and proliferation or cell cycle, and 52 DEGs were associated with cell death or cellular compromise. A large number were also associated with cellular signaling. Physiological system development and function categories were heavily populated with genes involved in vasculature formation, including hematological and cardiac system development and function, and tissue and organ morphology. These biological functions further emphasize that differences in signaling, cellular proliferation, microvasculature formation, and morphology all affect pituitary gland function in HG and LG birds.

## Upstream regulator analysis

Using IPA, upstream regulator analysis was performed to identify transcription factors with direct actions on DEGs in the pituitary glands of HG and LG birds. Transcriptional regulators associated with pituitary development and function that potentially contribute to differential expression of pituitary hormones are shown in Fig. 10a. Predicted upregulation of paired box gene 8 (PAX8) and early growth response protein 1 (EGR1) in HG birds is likely responsible for differential expression of genes associated with cell turnover that could change pituitary structure and morphology [collagen, type I, alpha 2 (COL1A2), clusterin (CLU), cyclin dependent kinase inhibitor 1C (CDKN1C), CASP8], as well as microvascular formation (FLT1). Other transcription factors identified as upstream regulators in this network influence pituitary cell-type

<footnote><sup>a</sup>As determined by Ingenuity Pathway Analysis gene interaction network prediction</footnote>

<footnote><sup>b</sup>The number of DEGs between HG and LG birds from our dataset contained within each network</footnote>

<footnote><sup>c</sup>The total number of genes contained within each network</footnote>
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Functional gene interaction network associated with pituitary hormone production in HG and LG chickens. A gene network showing direct (solid lines) and indirect (dashed lines) relationships associated with "Molecular Transport" was identified using Ingenuity Pathway Analysis software. Genes which are colored red (upregulated in HG line) or green (upregulated in LG line) are contained within the dataset of 291 genes with a significant line-by-age interaction or a significant main effect of line (P ≤ 0.05; n = 4 birds per line) and a difference in expression of at least 1.6-fold between lines. The intensity of the red or green colored gene symbols indicates the magnitude of the difference between genetic lines. The legend indicates the functional type of each gene</figcaption></figure>

differentiation and hormone expression, indicating that these could also be differentially expressed between HG and LG birds. A second upstream regulator identified by IPA was nuclear receptor subfamily 3, group C, member 1 (NR3C1) or glucocorticoid receptor, which is predicted to be upregulated in HG birds on week 1 and in LG birds at later ages (Fig. 10b). LG birds exhibit a substantial increase in expression of genes associated with immunity, inflammation, and other cellular stressors after week 1, such as B2M, DDIT4, interferon regulatory factor 1

(IRF1), 2'-5'-oligoadenylate synthetase like (OASL), and PPP3R1. This may be a result of increased pituitary glucocorticoid receptor activity in these birds as compared to those in the HG line on weeks 3, 5, and 7.

## Discussion

Transcriptional profiling of anterior pituitary gene expression during the juvenile growth period of broilers selected for high or low body weight has confirmed an important role for this gland in regulation of growth and
<figure><img src="image_10.png" /><figcaption>**Fig. 7** RT-qPCR confirmation of DEGs contained within the gene interaction network associated with pituitary hormone production. Microarray (a - h, top graph) and RT-qPCR (a - h, bottom graph) data are shown for each gene on post-hatch weeks 1, 3, 5, and 7. For RT-qPCR data, levels of mRNA for each gene were normalized to levels of GAPDH mRNA. Values (mean ± SEM; n = 4 birds per line) are expressed relative to the age and line with the highest expression level for each technique (set to 100% prior to log2-transformation for statistical analysis). The RT-qPCR data for GH (a, bottom), FSHβ (b, bottom), and LHβ (c, bottom) are also shown in Fig. 2. When the two-way ANOVA demonstrated a significant line-by-age interaction (*P* ≤ 0.05), the presence of an asterisk (*) denotes a significant difference between the lines at the indicated age (*P* ≤ 0.05). Significant main effects of line and/or age (*P* ≤ 0.05) are indicated at the top of each graph when the line-by-age interaction was not significant (*P* > 0.05).</figcaption></figure>

metabolism in chickens. A total of 363 DEGs were identified as being developmentally regulated with age in both lines (72 DEGs) or potentially contributing to phenotypic differences between HG and LG birds (291 DEGs). Of most interest are those genes in the latter category that exhibit overall line or line-by-age interactive effects, and gene expression profiles for 19 of these DEGs were confirmed by RT-qPCR. Several approaches were taken to analyze transcriptional profiling data in order to determine how these DEGs may contribute to

alterations in pituitary function and, ultimately, the divergent growth and compositional phenotypes.

Differences in pituitary hormone expression between the lines indicate that LG birds exhibit compromised pituitary function that could be manifest in altered production and secretion of several important pituitary hormones. Reduced levels of pituitary *TSH* mRNA in LG birds are consistent with our observation that circulating thyroid hormones are lower in these birds and findings of others that T₃, and to a lesser extent T₄, positively regulate overall body growth in chickens [26, 27] and may contribute
<figure><img src="image_11.png" /><figcaption>**Fig. 8** Functional gene interaction network associated with pituitary gland structure in HG and LG chickens. A gene network showing direct (solid lines) and indirect (dashed lines) relationships associated with "Organ Morphology and Cellular Compromise" was identified using Ingenuity Pathway Analysis software. Genes which are colored red (upregulated in HG line) or green (upregulated in LG line) are contained within the dataset of 291 genes with a significant line-by-age interaction or a significant main effect of line (P ≤ 0.05; n = 4 birds per line) and a difference in expression of at least 1.6-fold between lines. The intensity of the red or green colored gene symbols indicates the magnitude of the difference between genetic lines. The legend indicates the functional type of each gene</figcaption></figure>

to increased abdominal fat in HG birds [28]. Differences in expression of transcripts known to be regulated by thyroid hormones (e.g. thyroid hormone responsive spot 14 alpha (*THRSPA*), as well as enzymes responsible for activation and inactivation of *T*₃ (e.g. deiodinase (*DIO*) 1 and *DIO3*), were observed in adipose tissue from these birds [15], demonstrating that alterations in pituitary TSH production are driving differential thyrotropic axis activity that directly contributes to divergence in growth and adiposity. Pituitary mRNA levels of *DIO2*, another enzyme

responsible for local conversion of *T*₄ to *T*₃, were different between the lines. This indicates that, in addition to differences in circulating levels of thyroid hormones, there is also a potential difference between the lines related to local thyroid hormone action at the pituitary gland that could influence production of hormones such as TSH and GH. Despite an apparent increase in the number of pituitary somatotrophs in fast-growing lines of birds during late embryonic development [29], it is well established that circulating GH is elevated in slower-growing lines of
<figure><img src="image_12.png" /><figcaption>**Fig. 9** RT-qPCR confirmation of DEGs within the gene interaction network associated with pituitary gland structure. Microarray (a - e, top graph) and RT-qPCR (a - e, bottom graph) data are shown for each gene on post-hatch weeks 1, 3, 5, and 7. For RT-qPCR data, levels of mRNA for each gene were normalized to levels of GAPDH mRNA. Values (mean ± SEM; n = 4 birds per line) are expressed relative to the age and line with the highest expression level for each technique (set to 100% prior to log₂-transformation for statistical analysis). When the two-way ANOVA demonstrated a significant line-by-age interaction (*P* ≤ 0.05), the presence of an asterisk (*) denotes a significant difference between the genetic lines at the indicated age (*P* ≤ 0.05). Significant main effects of genetic line and/or age (*P* ≤ 0.05) are indicated at the top of each graph when the line-by-age interaction was not significant (*P* > 0.05)</figcaption></figure>

chickens after hatch [30, 31], and differences observed in pituitary GH mRNA levels in the current study were consistent with this. Pituitary GH mRNA was also shown to be higher in lighter birds 4 weeks after hatch using a similar genetic model system in which HWS birds weigh approximately 10-fold more than LWS birds at 8 weeks of

age [8]. It has been demonstrated in different birds from these same lines that circulating GH levels were 2.5-fold higher in LG chickens than HG chickens (unpublished observation by LAC, TEP, JS, and MJD). Despite increased pituitary GH, LG birds exhibit reduced levels of both hepatic IGF1 and IGF2 mRNA as well circulating IGFs [14].
<table><caption>**Table 4** Top biological functions of DEGs within pituitary glands of HG and LG chickens</caption><thead><tr><th><sup>a</sup>Top Biological Function</th><th><sup>b</sup>DEGs (#)</th><th><sup>c</sup>P-Value</th></tr></thead><tbody><tr><td><strong>Molecular and cellular function</strong></td><td></td><td></td></tr><tr><td>Cell death</td><td>46</td><td>9.96E-05 – 2.25E-02</td></tr><tr><td>Cell cycle</td><td>19</td><td>1.27E-04 – 2.25E-02</td></tr><tr><td>Cell-to-cell signaling and interaction</td><td>26</td><td>6.86E-04 – 2.25E-02</td></tr><tr><td>Cellular compromise</td><td>6</td><td>3.87E-04 – 2.09E-02</td></tr><tr><td>Cellular growth and proliferation</td><td>46</td><td>1.24E-03 – 2.20E-02</td></tr><tr><td><strong>Physiological system development and function</strong></td><td></td><td></td></tr><tr><td>Cardiovascular system development and function</td><td>16</td><td>5.36E-05 – 2.25E-02</td></tr><tr><td>Organ morphology</td><td>12</td><td>5.36E-05 – 1.13E-02</td></tr><tr><td>Hematological system development and function</td><td>32</td><td>1.16E-04 – 2.05E-02</td></tr><tr><td>Tissue morphology</td><td>24</td><td>1.16E-04 – 2.05E-02</td></tr><tr><td>Endocrine system development and function</td><td>6</td><td>7.13E-04 – 2.25E-02</td></tr></tbody></table>

This suggests that the slower-growing LG birds may be deficient in the hepatic response to GH, which may be a result of decreased expression of GH receptor (GHR) in liver or disruption of intracellular GHR signaling. In fact, hepatic GH-binding in LG chickens was only one-seventh that of the HG birds between 5 and 11 weeks post-hatch (unpublished observation by LAC, JS, and MJD). These findings are similar to those reported earlier in sex-link dwarf (dw/dw) chickens, which have elevated plasma GH, no detectable hepatic GH-binding activity, and maintain two-thirds of the plasma IGF-1 levels despite lacking a functional GHR gene [32–35]. In contrast, hepatic GHR mRNA expression was observed to be higher early post-hatch in LWS chickens and only modestly reduced at 4 weeks of age when compared to HWS chickens, and it was reported that GHR mRNA levels in breast muscle are higher in LWS birds during both embryogenesis and post-hatch [8]. Others have found that breast muscle cells from LG birds exhibit reduced sensitivity to IGF1 stimulation [18], and it has been suggested that their visceral fat may have a reduction of IGF1 signaling when compared to HG birds [15]. Taken together, it is apparent that differences in pituitary hormone production and downstream actions of these hormones play a major role in altering metabolic phenotypes in these birds.

Differential expression of receptors and intracellular signaling molecules may contribute to differences in expression of pituitary hormones between HG and LG birds. Midkine (MDK) is a secreted protein produced by folliculostellate cells within the embryonic rat pituitary gland [36], and its receptor, *PTPRZ1*, was more highly expressed in HG birds after divergence in growth and body composition. *PTPRZ1* expression has recently been

detected in the adult rat anterior pituitary [37], specifically within ACTH- and GH-producing cells, and the authors speculate that it mediates paracrine MDK signaling within these cell types. It is possible that MDK-PTRPZ signaling may be a novel regulator of these hormones, and differences in *PTPRZ1* expression between HG and LG birds may result in altered pituitary GH production. Many pituitary hormone releasing and release-inhibiting factors secreted by the hypothalamus activate G protein coupled receptors, which signal through generation of second messengers such as cyclic adenosine monophosphate (cAMP) and calcium, or MAPK pathways [38–40]. RASD1 has been shown to interfere with cAMP-stimulated peptide hormone secretion in a corticotroph cell line [41], and elevated RASD1 mRNA levels in LG birds suggests that they have increased RASD1 activity that may contribute to the reduced expression of *TSHβ*, *FSHβ*, and *LHβ* mRNA in a similar manner. An increase in *ADCY2* expression in these birds may increase intracellular cAMP in an attempt to maintain hormone expression levels in the face of this interference. Further supporting differences in cAMP-mediated intracellular signaling between the lines is the observation that LG birds exhibit elevated levels of *AKAP7* mRNA, a scaffolding protein that binds to regulatory subunits of cAMP-activated protein kinase A and influences its activity [42]. Increased expression of mRNA for a regulatory subunit of calcineurin, *PPP3R1*, in LG birds indicates that calcium-mediated intracellular signaling may also be altered between the lines and influence pituitary hormone expression levels and/or secretion. Activin-mediated induction of FSH production is regulated by MAP3K7 in sheep pituitary cells [43], and elevated levels

<footnote><sup>a</sup>As determined by Ingenuity Pathway Analysis (IPA) biological function prediction</footnote>

<footnote><sup>b</sup>The number of differentially expressed genes between HG and LG birds from our dataset contained within each category</footnote>

<footnote><sup>c</sup>P-values were determined by IPA software based on the number of DEGs from the dataset within each biological category divided by the total number of known genes assigned to that category within the database</footnote>
a

Week 1

<figure><img src="image_13.png" /></figure>

Week 3

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

Week 5

<figure><img src="image_16.png" /></figure>

Week 7

<figure><img src="image_17.png" /></figure>

b

Week 1

<figure><img src="image_18.png" /></figure>

Week 3

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

Week 5

<figure><img src="image_21.png" /></figure>

Week 7

<figure><img src="image_22.png" /><figcaption>**Fig. 10** Transcriptional regulators predicted to regulate DEGs observed in the pituitary of HG and LG chickens. Upstream regulator analysis using Ingenuity Pathway Analysis was performed on the dataset of 291 genes with a significant line-by-age interaction or a significant main effect of line (P ≤ 0.05; n = 4 birds per line) and a difference in expression of at least 1.6-fold between lines. This analysis identified transcription factors with direct actions on differentially expressed target genes. **a** Transcriptional regulators associated with pituitary development and function that could contribute to differential expression of pituitary hormones were identified. **b** A large number of differentially expressed genes between the two lines are also predicted to be regulated by glucocorticoid receptor (NR3C1). Upstream regulators colored orange are predicted to have higher activity in HG birds, and upstream regulators colored blue are predicted to have higher activity in LG birds. Genes within the dataset that are colored red were more highly expressed in the HG line, and genes within the dataset that are colored green were more highly expressed in the LG line. The color of the line connecting transcriptional regulators to differentially expressed genes indicates IPA prediction of whether transcription factor activity should lead to upregulation (orange line) or downregulation (blue line) of the DEG targets</figcaption></figure>
of MAP3K7 mRNA in HG birds suggests that this pathway may play a similar role in birds and also be a positive regulator of pituitary hormones. In addition to differential cell signaling, HG birds exhibited higher levels of mRNA for genes associated with the basic cellular functions of transcription and translation, such as those involved in processing nascent mRNA transcripts [heterogeneous nuclear ribonucleoprotein A2/B1 (HNRNPA2B1)] [44], and chaperones involved in the misfolded protein response (CALR and HSPA5) [45]. Alterations in intracellular signaling within pituitary glands of HG and LG birds, as well as the ability to process newly transcribed and translated genes, likely leads to differences in pituitary hormone production between the lines.

In addition to differences in pituitary hormone expression, it is clear that hormone secretion may also be altered, as HG and LG birds exhibit differential expression of several genes with roles in molecular transport and calcium-dependent membrane trafficking. The CHGB gene is one of the most abundantly transcribed genes in the pituitary gland [46] and is a major component of dense-core secretory vesicles within endocrine, neuronal, and other secretory cell types [47]. The observation that CHGB expression was elevated in the anterior pituitary of HG birds indicates that the formation of secretory vesicles containing pituitary hormones may be compromised in birds from the LG line. Vesicular transport, fluctuation in intracellular calcium levels, and membrane fusion events are necessary for proper hormone secretion from these dense-core vesicles, and differences in expression levels for genes likely involved in these processes are evident between HG and LG birds. ANKRD27 is thought to interfere with vesicle-associated membrane protein 7-mediated vesicular fusion by trapping it in a conformation unable to interact with snap receptor proteins on the opposing membrane [48]. This interaction is crucial for proper membrane fusion and release of vesicular contents at the plasma membrane, and increased levels of ANKRD27 mRNA indicate that this process may be less efficient and interfere with hormone secretion in LG birds. In mammalian cells, CPNE1, MAPRE2, an ANXA2 are thought to play a role in proper calcium-induced trafficking to and docking with the plasma membrane [47, 49, 50]. Upregulation of mRNA expression for these genes was observed in LG pituitary glands and may be a compensation for reduced synthesis and/or functionality of hormone-secreting vesicles in these birds.

Genes with a role in antigen presentation (B2M, MHC, and HLA-A) [51], the antiviral response (OAS3) [52], cytokine signaling and transcriptional regulation (STAT1) [53], apoptosis (CASP8) [54], and response to cellular stressors such as hypoxia and DNA damage response (DDIT4) [55] were elevated at later ages in the pituitary of

LG birds. This may lead to complications similar to that seen in lymphocytic hypophysitis, which is characterized by immune cell infiltration, inflammation, and damage to pituitary cells that can result in differing levels of hypopituitarism [56–58]. Interestingly, in addition to its role in response to cellular stressors, DDIT4 has recently been implicated in a wide range of cellular processes which impact energy homeostasis and metabolic function [55]. Of note is that many of these genes are upregulated in LG birds after week 3, when the metabolic and growth phenotypes diverge. Similarly, after 3 weeks of age, pituitary glucocorticoid receptor activity is predicted to be upregulated in LG birds. This may be a response to alterations in expression of genes associated with cellular stressors and further interfere with normal pituitary function. Elevated pituitary glucocorticoid receptor activity in LG birds suggests that these birds may have higher levels of circulating corticosterone, leading to metabolic alterations that negatively impact growth.

Results from several studies have demonstrated that ALVE transcripts are overexpressed in several tissues from slower growing chickens at most ages examined [15, 59–62]. In the present study, ALVE mRNA expression was higher in LG chickens from week 3 onwards, as determined by both global transcriptional profiling and RT-qPCR (Fig. 8). In fact, it was one of the most highly DEGs between the genetic lines. Transcriptional analysis of important somatic metabolic tissues such as abdominal fat, liver, and breast muscle has also revealed that ALVE expression is highly upregulated in LG chickens [15, 59]. In a second, independently selected genetic model system, slower growing LWS chickens were observed to exhibit higher levels of ALVE mRNA in the hypothalamus, whole brain, liver, breast muscle, and adipose tissue [60–62]. These authors determined that more ALVE integration sites were detected in LWS than HWS chickens and contributed to differences in expression between the lines [62]. It is likely a similar phenomenon has occurred in our genetically divergent LG versus HG birds. While it is not known if this differential expression is a direct response to selection for high or low growth (i.e. whether the ALVE sequences are directly affecting growth or are linked to loci regulating growth), it is clear that presence and expression of the ALVE allele are highly correlated with reduced growth. In addition to the above possibilities, multiple integrations of ALVE may disrupt key functional genes that control growth and metabolism in chickens.

Several categories of differentially expressed genes indicate that it is likely there are differences in structure and organization within the neuroendocrine system of HG and LG broiler chickens. The homeobox transcription factor UNCX is involved in developing neuronal contacts between the hypothalamus and pituitary [63]
and was differentially expressed between the lines, providing evidence that this connectivity may differ in HG and LG birds. Semaphorins are a class of secreted and membrane-tethered molecules that were first identified for their role in axon guidance, and plexins act as their cell-surface receptors [64]. Semaphorin-plexin D1 signaling has been shown to mediate VEGF-induced angiogenesis in mouse retina [65], with loss of signaling resulting in compromised vasculature formation. Both SEMA3C and its cognate receptor, PLXND1, were expressed at higher levels in the anterior pituitary of HG birds. The VEGF receptor FLT1 was also expressed at higher levels in HG birds at all ages. The Wnt protein receptor FZD4 has been shown to be necessary for blood vessel morphogenesis [66] and was expressed at higher levels in HG birds at later ages. Elevated PTPRZ1 mRNA levels in the pituitary gland of HG birds may also contribute to alterations in MDK-PTPRZ1 signaling leading to differences in pituitary gland structure and organization in addition to putative regulation of hormone production, as MDK has been shown to positively influence cell proliferation, cell migration, and angiogenesis [67]. Together, these results indicate that there may be differences between the lines related to connectivity between hypothalamic neurons, the hypophyseal portal vasculature, and the anterior pituitary gland. Alterations in microvascular formation and support may contribute to differences in metabolic feedback to the central neuroendocrine system between the genetic lines, as well as affect delivery of hypothalamic inputs to the pituitary gland or the delivery of anterior pituitary hormones into peripheral circulation.

## Conclusions

While intensive genetic selection within the broiler industry has led to rapid advances in growth rate, body composition (improved muscle accretion accompanied by increases in abdominal fat), and feed efficiency, there is still a lack of understanding related to endocrine control of growth and metabolism in these animals. As animal production systems are forced to reduce the use of antibiotics due to consumer demand and regulatory requirements, understanding biological mechanisms governing growth and metabolism in broiler chickens is crucial to maintaining production efficiency. In this study, transcriptional profiling followed by several bioinformatics approaches were used to evaluate gene expression within anterior pituitary glands of juvenile male broiler chickens selected for high or low body weight in order to identify endocrine mechanisms regulating phenotypic differences between these animals. A total of 263 candidate genes with a potential role in altering growth and metabolism were identified, and the results clearly demonstrate that molecular events associated with all levels of pituitary gland function were impacted by divergent selection. Among the candidate genes

were five of the six pituitary hormones. Organization of the dataset using SOMs identified clusters of genes which could be driving phenotypic differences through altered regulation of pituitary hormone production. GO categorization and pathway analyses confirmed that processes and gene networks related to signal transduction, transcriptional regulation, and membrane trafficking and vesicle-mediated transport were different between lines and may be influencing hormone expression and secretion. Differential expression of genes involved in immune regulation was observed, suggesting that processes such as inflammation and response to cellular stressors may compromise pituitary function. Finally, genes playing a role in processes related to morphology and angiogenesis were highlighted, providing evidence that pituitary function is intimately tied to structure and that pituitary gland organization will influence hypothalamic and systemic metabolic inputs and delivery of hormones regulating growth and metabolism into peripheral circulation.

## Methods

### Animals and sample collection

Meat-type chickens maintained at INRA (Nouzilly, France) that have been divergently selected for high or low juvenile body weight at 8 weeks of age were used in this study [6]. Samples used in this study were collected from male birds that were reared as previously reported [15]. Briefly, birds were raised in floor pens (4.4 m × 3.9 m) according to standard broiler practice and given ad libitum access to water and conventional starter (weeks 0–3; 22% crude protein and 3050 kcal ME/kg) or grower (weeks 3–11; 20% crude protein, and 3100 kcal) rations. To reduce early mortality of LG birds and minimize environmental differences, the two lines were raised separately through 3 weeks of age and then placed together and raised as a mixed population for the remainder of the study. At post-hatch weeks 1, 3, 5, and 7, randomly selected birds from each group (n = 8 per line) were weighed, and blood was collected from the brachial vein using heparinized syringes for plasma isolation. Following cervical dislocation, abdominal fat pad weights were determined, and anterior pituitary glands were collected, immediately snap-frozen in liquid nitrogen, and stored at -80 °C until RNA extraction. Plasma levels of total T₃ and T₄ were determined using commercially available coated tube radioimmunoassay kits (MP Biomedicals, Solon, OH). All procedures were approved by the Institutional Animal Care and Use Committee at each of the four institutions (INRA, University of Maryland, University of Delaware, and University of Georgia).

### RNA isolation and amplification

Anterior pituitary total RNA was isolated from four birds of both lines at each of four ages (1, 3, 5, and 7
weeks post-hatch; *n* = 4) using the RNeasy Mini Kit (Qiagen, Valencia, CA) with on-column DNase digestion to eliminate potential genomic DNA contamination. Isolated total RNA was quantified with a spectrophotometer by measuring absorbance at 260 nm, and the quality was verified with a bioanalyzer (Agilent Technologies, Palo Alto, CA). As chicken anterior pituitary glands do not yield sufficient total RNA for direct analysis by microarray, a previously detailed and validated [68–70] modification of the Eberwine procedure [71] was used to amplify mRNA for hybridization. Briefly, total RNA (0.5 µg) was reverse transcribed with SuperScript II (Invitrogen, Carlsbad, CA) using an oligo (dT) primer containing a T7 promoter site (5’-GGCCAGTGAATTGT AATACGACTCACTATAGGGAGGCGGT<sub>24</sub>-3’; Affymetrix, Santa Clara, CA). After second-strand synthesis, the double-stranded cDNA was phenol-chloroform extracted, purified using a Microcon-30 column (Millipore, Billerica, MA), and used as a template for in vitro transcription with the T7 MEGAscript Kit (Ambion, Austin, TX) according to the manufacturer’s protocol. The resulting amplified RNA (aRNA) was phenol-chloroform extracted, purified with a Spin Column-30 (Sigma, St. Louis, MO), and quantified using the Ribogreen RNA Quantitation Kit (Invitrogen).

## Microarray hybridization, processing, and analysis

The Del-Mar 14 K Chicken Integrated Systems Microarray (Geo Platform accession no. GPL1731), an annotated cDNA array that has been previously described [24, 72] was used for transcriptional profiling of anterior pituitary mRNA expression in HG and LG chickens. A total of 32 individual animals were analyzed, with four replicate samples (*n* = 4) from each of the two lines (HG and LG) at each of the four ages (weeks 1, 3, 5, and 7). A reference hybridization design [73] was used for microarray analysis, where an internal reference pool was generated from equal amounts of aRNA from each sample. Individual aRNA samples from each animal were labeled with Cy3, an aliquot of the reference sample was labeled with Cy5, and equal amounts of each (1 µg) were hybridized together on a slide. This design necessitated the use of 32 microarray slides (one per individual pituitary gland). Labeling with Cy3 and Cy5, microarray hybridization, and image scanning were performed at the University of Maryland Biotechnology Institute’s Microarray Core Facility as previously described [68, 69].

Data from the microarray analysis were processed, normalized, and trimmed as described earlier [28, 68, 69] using software that is part of the TM4 suite of microarray data analysis applications [74] freely available from The Institute for Genomic Research (TIGR, Rockville, MD). During image processing, cDNA spots that were flagged due to lack of detection, detection below

background, pixel intensity saturation, or malformation were rejected from further processing and normalization and removed from the dataset. Spots whose pixel intensity for all slides was < 90% of the lowest median pixel intensity for the salmon DNA control spots (background; 8 spots/slide) among all slides were also eliminated from further consideration, as were any spots that were not detectable on at least half of the slides. Microarray data files for the study were deposited in the GEO data repository (accession no. GSE122519, sample accession nos. GSM3473153 – GSM3473183).

The trimmed dataset consisted of 10,437 cDNA spots that were submitted for statistical analysis by two-way analysis of variance (ANOVA) using the general linear model (GLM) procedure of Statistical Analysis System (SAS Institute, Cary, NC) v. 9.4. Data were analyzed as log₂(normalized Cy3/raw Cy5), or log₂-ratio, for each spot, and differences were considered statistically significant at *P* ≤ 0.05. Genes with significant differences for the line-by-age interaction, main effect of line, and main effect of age are listed in Additional file 2. To further filter out false positives, these 2310 spots were considered identified as DEGs if the fold-difference among experimental groups was ≥1.6 (0.68 on a log₂ scale). This resulted in 291 DEGs of interest for the line-by-age interaction or main effect of line, which were further subjected to downstream analysis for GO categorization, SOMs clustering, and IPA gene network and upstream regulator analysis.

For purposes of comparison between microarray expression data and RT-qPCR expression data (see below), data were transformed as 2^(log2ratio) and divided by the mean value of the age and line with the highest expression for that spot so that data are expressed relative to the age and line with the highest expression level (set to 100%). Relative expression data for individual values were log₂-transformed prior to analysis using two-way ANOVA, and differences between groups were determined with the test of least significant difference (SAS) when the overall *P*-value for a given effect (genetic line-by-age, genetic line, or age) was statistically significant (*P* ≤ 0.05).

## Gene ontology and self-organizing maps analyses

DEGs were analyzed for GO enrichment terms based on biological process and molecular function using AgBase (http://agbase.arizona.edu/), a resource for functional analysis of agricultural gene products [75]. Of the 260 functionally annotated DEGs, 249 and 236 were successfully placed into GO categories based on biological process and molecular function, respectively. GeneCluster 2.0 [76] was used to organize the 291 DEGs into 16 SOMs clusters in a 4 × 4 configuration. Data were entered into the program as log₂ratio<sub>HGmean</sub> – log₂ratio<sub>LGmean</sub> at each age, in order to identify changes in relative expression between the lines which may be correlated with divergence in growth and
cell kinase; HG: High growth; HLA-A: Human leukocyte antigen A; HNRNPA2B1: Heterogeneous nuclear ribonucleoprotein A2/B1; HSPA5: Heat shock protein family A, member 5; HWS: High weight selected; IGF: Insulin-like growth factor; INRA: Institut National de la Recherche Agronomique; IPA: Ingenuity Pathway Analysis; IRF1: Interferon regulatory factor 1; LG: Low growth; LH: Luteinizing hormone; LWS: Low weight selected; MAP3K7: Mitogen-activated protein kinase kinase kinase 7; MAPRE2: Microtubule associated protein RP/EB family member 2; MDK: Midkine; MEK1: Mitogen-activated protein kinase/extracellular signal-regulated kinase kinase 1; MHC: Major histocompatibility complex; NR0B1: Nuclear receptor subfamily 0, group B, member 1; NR3C1: Nuclear receptor subfamily 3, group C, member 1 (glucocorticoid receptor); OAS3: 2'-5'-oligoadenylate synthetase 3; OASL: 2'-5'-oligoadenylate synthetase like; PLXND1: Plexin D1; PODXL: Podocalyxin-like protein 1; POMC: Pro-opiomanocortin; PPP3R1: Protein phosphatase 3 regulatory subunit B alpha; PRL: Prolactin; PTPRZ1: Protein tyrosine phosphatase, receptor type Z1; RASD1: Dexamethasone-induced ras-related protein 1; RT-qPCR: Reverse transcription-quantitative PCR; SEMA3C: Semaphorin 3C; SEMA3F: Semaphorin 3F; SMAD1: Mothers against decapentaplegic homolog 1; SOMS: Self-organizing maps; STAT1: Signal transducer and activator of transcription 1; STMN1: Stathmin 1; T3: 3,5,3'-triiodothyronine; T4: Thyroxine; TCR: T-cell receptor; TSH: Thyroid-stimulating hormone; UNCX: UNC homeobox; VEGF: Vascular endothelial growth factor

## Acknowledgements

The longitudinal animal study and associated sample collection required the participation of several scientists and technicians from laboratories and breeding facilities at the Unite de Recherches Avicoles and Unite Experimentale Pole d'Experimentation Avicole de Tours, INRA, Nouzily, France. Production of microarrays used for this study also required the participation of several scientists, post-doctoral associates, technicians, and students from laboratories at University of Maryland, University of Delaware, and Institut National de la Recherche Agronomique. Although not individually named, each person's contribution towards this project was essential and is gratefully acknowledged.

## Funding

This project was supported by United States Department of Agriculture, Cooperative State Research, Education, and Extension Service, Initiative for Future Agriculture and Food Systems (USDA-CREES-IFAFS) Program Grant No. 00-52100-9614. The funding agency provided financial support for the research but was not involved in the design of the study nor collection, analysis, and interpretation of data or in writing the manuscript.

## Availability of data and materials

MIAME-compliant microarray data from the longitudinal analysis of gene expression in anterior pituitary glands of HG and LG male broiler chickens were deposited in the NCBI Gene Expression Omnibus database as Series GSE122519.

## Authors' contributions

LEE assisted with production of the microarrays, extracted pituitary total RNA, analyzed microarray data, conducted all bioinformatics analyses, performed RT-qPCR, and drafted the manuscript. TEP isolated the anterior pituitary glands. LAC and TEP designed, constructed, and printed the Del-Mar 14 K Chicken Integrated Systems Microarray. LAC, TEP, SEA, and JS secured funding and designed the experiment. JS, ELBD, and MJD carried out the animal experiment. MSB performed microarray image processing and data normalization. MJD was involved in the design of the study and interpretation of the results. All authors reviewed and approved the final version of the manuscript.

## Ethics approval and consent to participate

Husbandry and experimental procedures for birds was performed according to the European Guidelines for Animal Care under authorization #37-04361 issued to Jean Simon. All animal procedures were approved by Veterinary Services of the Indre-et-Loire Department. All procedures were also in compliance with Institutional Care and Use Committee Guidelines at University of Maryland, University of Delaware, and University of Georgia.

## Consent for publication

Not applicable.

## Competing interests

The authors declare that they have no competing interests.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Author details

¹Department of Poultry Science, University of Georgia, Athens, GA 30602, USA. ²Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA. ³Department of Animal and Food Sciences, University of Delaware, Newark, DE 19716, USA. ⁴Biologie des Oiseaux et Aviculture, Institut National de la Recherche Agronomique (INRA), Université de Tours, UR83 Recherches Avicoles, 37380 Nouzilly, France.

Received: 20 December 2018 Accepted: 8 April 2019

Published online: 25 April 2019

## References

1. Sam S, Frohman LA. Normal physiology of hypothalamic pituitary regulation. Endocrinol Metab Clin N Am. 2008;37:1-22.

2. Harvey S. Growth hormone and growth? Gen Comp Endocrinol. 2013;190:3-9.

3. Scanes CG. Perspectives on the endocrinology of poultry growth and metabolism. Gen Comp Endocrinol. 2009;163:24-32.

4. Dunnington EA, Siegel PB. Long-term divergent selection for eight-week body weight in white Plymouth rock chickens. Poult Sci. 1996;75:1168-79.

5. Liu G, Dunnington EA, Siegel PB. Responses to long-term divergent selection for eight-week body weight in chickens. Poult Sci. 1994;73:1642-50.

6. Ricard FH, Marche G, Bourdiol R, Poisson H. Essai de sélection sur la forme de la courbe de croissance chez le poulet. Dispositif expérimental et premiers résultats d'ensemble. Ann Genet Sel Anim. 1975;7:427-43.

7. Dunnington EA, Honaker CF, McGillard ML, Siegel PB. Phenotypic responses of chickens to long-term, bidirectional selection for juvenile body weight—historical perspective. Poult Sci. 2013;92:1724-34.

8. Wu G, Siegel PB, Gilbert ER, Yang N, Wong EA. Expression profiles of somatotropin axis genes in lines of chickens divergently selected for 56-day body weight. Anim Biotechnol. 2011;22:100-10.

9. Yin H, Zhang S, Gilbert ER, Siegel PB, Zhu Q, Wong EA. Expression profiles of muscle genes in postnatal skeletal muscle in lines of chickens divergently selected for high and low body weight. Poult Sci. 2014;93:147-54.

10. Sumners LH, Zhang W, Zhao X, Honaker CF, Zhang S, Cline MA, Siegel PB, Gilbert ER. Chickens from lines artificially selected for juvenile low and high body weight differ in glucose homeostasis and pancreas physiology. Comp Biochem Physiol A Mol Integr Physiol. 2014;172:57-65.

11. Zhang S, McMillan RP, Hulver MW, Siegel PB, Sumners LH, Zhang W, Cline MA, Gilbert ER. Chickens from lines selected for high and low body weight show differences in fatty acid oxidation efficiency and metabolic flexibility in skeletal muscle and white adipose tissue. Int J Obes. 2014;38:1374-82.

12. Parker GA, Sumners LH, Zhao X, Honaker CF, Siegel PB, Cline MA, Gilbert ER. Delayed access of low body weight-selected chicks to food at hatch is associated with up-regulated pancreatic glucagon and glucose transporter gene expression. Comp Biochem Physiol A Mol Integr Physiol. 2015;189:124-9.

13. Mignon-Grasteau S, Beaumont C, Ricard FH. Genetic analysis of a selection experiment on the growth curve of chickens. Poult Sci. 2001;80:849-54.

14. Beccavin C, Chevalier B, Cogburn LA, Simon J, Duclos MJ. Insulin-like growth factors and body growth in chickens divergently selected for high or low growth rate. J Endocrinol. 2001;168:297-306.

15. Resnyk CW, Carré W, Wang X, Porter TE, Simon J, Le Bihan-Duval E, Duclos MJ, Aggrey SE, Cogburn LA. Transcriptional analysis of abdominal fat in chickens divergently selected on bodyweight at two ages reveals novel mechanisms controlling adiposity: validating visceral adipose tissue as a dynamic endocrine and metabolic organ. BMC Genomics. 2017;18:626.

16. Remignon H, Gardahaut MF, Marche G, Ricard FH. Selection for rapid growth increases the number and the size of muscle fibres without changing their typing in chickens. J Muscle Res Cell Motil. 1995;16:95-102.

17. Remignon H, Lefaucheur L, Blum JC, Ricard FH. Effects of divergent selection for body weight on three skeletal muscles characteristics in the chicken. Br Poult Sci. 1994;35:65-76.

18. Duclos MJ, Chevalier B, Remignon H, Ricard FH, Goddard C, Simon J. Divergent selection for high or low growth rate modifies the response of
muscle cells to serum or insulin-like growth factor-I in vitro. Growth Regul. 1996;6:176-84.

19. Oudin A, Chevalier B, Simon J, Duclos MJ. Muscle insulin-like growth factor-I (IGF-I) receptors in chickens with high or low body weight: effects of age and muscle fibre type. Growth Hormone IGF Res. 1998;8:243-50.

20. Ankra-Badu GA, Le Bihan-Duval E, Mignon-Grasteau S, Pitel F, Beaumont C, Duclos MJ, Simon J, Carré W, Porter TE, Vignal A, et al. Mapping QTL for growth and shank traits in chickens divergently selected for high or low body weight. Anim Genet. 2010;41:400-5.

21. Nadaf J, Pitel F, Gilbert H, Duclos MJ, Vignoles F, Beaumont C, Vignal A, Porter TE, Cogburn LA, Aggrey SE, et al. QTL for several metabolic traits map to loci controlling growth and body composition in an F2 intercross between high- and low-growth chicken lines. Physiol Genomics. 2009;38:241-9.

22. Nadaf J, Gilbert H, Pitel F, Berri CM, Feve K, Beaumont C, Duclos MJ, Vignal A, Porter TE, Simon J, et al. Identification of QTL controlling meat quality traits in an F2 cross between two chicken lines selected for either low or high growth rate. BMC Genomics. 2007;8:155.

23. Cervantes HM. Antibiotic-free poultry production: is it sustainable? J Appl Poult Res. 2015;24:91-7.

24. Cogburn LA, Wang X, Carré W, Rejto L, Aggrey SE, Duclos MJ, Simon J, Porter TE. Functional genomics in chickens: development of integrated-systems microarrays for transcriptional profiling and discovery of regulatory pathways. Comp Funct Genomics. 2004;5:253-61.

25. Cogburn LA, Wang X, Carré W, Rejto L, Porter TE, Aggrey SE, Simon J. Systems-wide chicken DNA microarrays, gene expression profiling, and discovery of functional genes. Poult Sci. 2003;82:939-51.

26. Tixier-Boichard M, Decuypere E, Huybrechts L, Kühn E, Mérat P. Effects of dietary T3 on growth parameters and hormone levels in normal and sex-linked dwarf chickens. Domest Anim Endocrinol. 1990;7:573-85.

27. Xiao Y, Wu C, Li K, Gui G, Zhang G, Yang H. Association of growth rate with hormone levels and myogenic gene expression profile in broilers. J Anim Sci Biotechnol. 2017;8:43.

28. Byerly MS, Simon J, Cogburn LA, Le Bihan-Duval E, Duclos MJ, Aggrey SE, Porter TE. Transcriptional profiling of hypothalamus during development of adiposity in genetically selected fat and lean chickens. Physiol Genomics. 2010;42:157-67.

29. Porter TE. Differences in embryonic growth hormone secretion between slow and fast growing chicken strains. Growth Horm IGF Res. 1998;8:133-9.

30. Goddard C, Wilkie RS, Dunn IC. The relationship between insulin-like growth factor-1, growth hormone, thyroid hormones and insulin in chickens selected for growth. Domest Anim Endocrinol. 1988;5:165-76.

31. Mao JNC, Burnside J, Postel-Vinay MC, Pesek JD, Chambers JR, Cogburn LA. Ontogeny of growth hormone receptor gene expression in tissue of growth-selected strains of broiler chickens. J Endocrinol. 1998;156:67-75.

32. Agarwal SK, Cogburn LA, Burnside J. Dysfunctional growth hormone receptor in a strain of sex-linked dwarf chicken: evidence for a mutation in the intracellular domain. J Endocrinol. 1994;142:427-34.

33. Burnside J, Liou SS, Cogburn LA. Molecular cloning of the chicken growth hormone receptor complementary deoxyribonucleic acid: mutation of the gene in sex-linked dwarf chickens. Endocrinology. 1991;128:3183-92.

34. Burnside J, Liou SS, Zhong C, Cogburn LA. Abnormal growth hormone receptor gene expression in the sex-linked dwarf chicken. Gen Comp Endocrinol. 1992;88:20-8.

35. Huang N, Cogburn LA, Agarwal SK, Marks HL, Burnside J. Overexpression of a truncated growth hormone receptor in the sex-linked dwarf chicken: evidence for a splice mutation. Mol Endocrinol. 1993;7:1391-8.

36. Fujiwara K, Maliza R, Tofrizal A, Batchuluun K, Ramadhani D, Tsukada T, Azuma M, Horiguchi K, Kikuchi M, Yashiro T. In situ hybridization analysis of the temporospatial expression of the midkine/pleiotrophin family in rat embryonic pituitary gland. Cell Tissue Res. 2014;357:337-44.

37. Fujiwara K, Horiguchi K, Maliza R, Tofrizal A, Batchuluun K, Ramadhani D, Syaidah R, Tsukada T, Azuma M, Kikuchi M, et al. Expression of the heparin-binding growth factor midkine and its receptor, Pptrz1, in adult rat pituitary. Cell Tissue Res. 2015;359:909-14.

38. Naor Z, Benard O, Seger R. Activation of MAPK cascades by G-protein-coupled receptors: the case of gonadotropin-releasing hormone receptor. Trends Endocrinol Metab. 2000;11:91-9.

39. Naor Z, Huhtaniemi I. Interactions of the GnRH receptor with heterotrimeric G proteins. Front Neuroendocrinol. 2013;34:88-94.

40. Wang Y, Li J, Want CY, Kwok AY, Zhang X, Leung FC. Characterization of the receptors for chicken GHRH and GHRH-related peptides: identification of a

novel receptor for GHRH and the receptor for GHRH-LP (PRP). Domest Anim Endocrinol. 2010;38:13-31.

41. Graham TE, Key TA, Kilpatrick K, Dorin RL. Dexas1/AGS-1, a steroid hormone-induced guanosine triphosphate-binding protein, inhibits 3',5'-cyclic adenosine monophosphate-stimulated secretion in AtT-20 corticotroph cells. Endocrinology. 2001;142:2631-40.

42. Edwards AS, Scott JD. A-kinase anchoring proteins: protein kinase a and beyond. Curr Opin Cell Biol. 2000;12:217-222.

43. Safwat N, Ninomiya-Tsuji J, Gore AJ, Miller WL. Transforming growth factor beta-activated kinase 1 is a key mediator of ovine follicle-stimulating hormone beta-subunit expression. Endocrinology. 2005;146:4814-24.

44. He Y, Smith R. Nuclear functions of heterogeneous nuclear ribonucleoproteins a/B. Cell Mol Life Sci. 2009;66:1239-56.

45. Wang J, Lee J, Liem D, Ping P. HSPA5 gene encoding Hsp70 chaperone BiP in the endoplasmic reticulum. Gene. 2017;618:14-23.

46. Tanaka S, Tatsumi K, Okubo K, Itoh K, Kawamoto S, Matsubara K, Amino N. Expression profile of active genes in the human pituitary gland. J Mol Endocrinol. 2002;28:33-44.

47. Wacker I, Kaether C, Krömer A, Migala A, Almers W, Gerdes HH. Microtubule-dependent transport of secretory vesicles visualized in real time with a GFP-tagged secretory protein. J Cell Sci. 1997;110:1453-63.

48. Schäfer IB, Hesketh GG, Bright NA, Gray SR, Pryor PR, Evans PR, Luzio JP, Owen DJ. The binding of Varp to VAMP7 traps VAMP7 in a closed, fusogenically inactive conformation. Nat Struct Mol Biol. 2012;19:1300-9.

49. Gabel M, Delavoe F, Demais V, Royer C, Bailly Y, Vitale N, Bader MF, Chasserot-Golaz S. Annexin A2-dependent actin bundling promotes secretory granule docking to the plasma membrane and exocytosis. J Cell Biol. 2015;210:785-800.

50. Perestenko PV, Pooler AM, Noorbakhshnia M, Gray A, Bauccio C, McIlhinney RAJ. Copines-1, -2, -3, -6 and -7 show different calcium-dependent intracellular membrane translocation and targeting. FEBS J. 2010;277:5174-89.

51. Kalish RS, Askenase PW. Molecular mechanisms of CD8+ T cell-mediated delayed hypersensitivity: implications for allergies, asthma, and autoimmunity. J Allergy Clin Immunol Pract. 1999;103:192-9.

52. Choi UY, Kang JS, Hwang YS, Kim YJ. Oligoadenylate synthase-like (OASL) proteins: dual functions and associations with diseases. Exp Mol Med. 2015;47:e144.

53. Ivashkov LB. IFNy: signalling, epigenetics and roles in immunity, metabolism, disease and cancer immunotherapy. Nat Rev Immunol. 2018;18:545-58.

54. Keller N, Ozmadenci D, Ichim G, Stupack D. Caspase-8 function, and phosphorylation, in cell migration. Semin Cell Dev Biol. 2018;82:105-17.

55. Lipina C, Hundal HS. Is REDD1 a metabolic eminence grise? Trends Endocrinol Metab. 2016;27:868-80.

56. Ng WH, Gonzales M, Kaye AH. Lymphocytic hypophysitis. J Clin Neurosci. 2003;10:409-13.

57. Takahashi Y. Autoimmune hypothyroidism: new developments. Handb Clin Neuroan. 2014;124:417-22.

58. Carmichael JD. Update on the diagnosis and management of hypophysitis. Curr Opin Endocrinol Diabetes Obes. 2012;19:314-21.

59. Cogburn LA, Porter TE, Duclos MJ, Simon J, Burgess SC, Zhu JJ, Cheng HH, Dodgson JB, Burnside J. Functional genomics of the lnc-a model organism. Poult Sci. 2007;86:2059-94.

60. Ka S, Albert FW, Denbow DM, Pääbo S, Siegel PB, Andersson L, Hallböök F. Differentially expressed genes in hypothalamus in relation to genomic regions under selection in two chicken lines resulting from divergent selection for high or low body weight. Neurogenetics. 2011;12:211-21.

61. Ka S, Fitzsimmons C, Jacobsson L, Savolainen P, Lunderberg J, Siegel PB, Andersson L, Hallböök F. Expression analysis of growth and energy regulation-associated genes in two divergent chicken strains. Ann N Y Acad Sci. 2005;1040:357-9.

62. Ka S, Kerje S, Bornold L, Liljegren U, Siegel PB, Andersson L, Hallböök F. Proviral integrations and expression of endogenous avian leucosis virus during long term selection for high and low body weight in two chicken lines. Retrovirology. 2009;6:68.

63. Asbreuk CH, van Doorninck JH, Mansouri A, Smidt MP, Burbach JP. Neurohypophysial dysmorphogenesis in mice lacking the homeobox gene Uncx4.1. J Mol Endocrinol. 2006;36:65-71.

64. Raper JA. Semaphorins and their receptors in vertebrates and invertebrates. Curr Opin Neurobiol. 2000;10:88-94.

65. Kim J, Oh WJ, Gaiano N, Yoshida Y, Gu C. Semaphorin 3E-Plexin-D1 signaling regulates VEGF function in developmental angiogenesis via a feedback mechanism. Genes Dev. 2011;25:1399-411.

66. Marković R, Peltan J, Gosak M, Horvat D, Žalik B, Seguy B, Chauvel R, Malandain G, Couffinhal T, Duplåa C, et al. Planar cell polarity genes
frizzled4 and frizzled6 exert patterning influence on arterial vessel morphogenesis. PLoS One. 2017;12:e1371.

67. Weckbach LT, Preissner KT, Deindl E. The role of Midkine in Arteriogenesis, involving Mechanosensing, endothelial cell proliferation, and vasodilation. Int J Mol Sci. 2018;19:2559.

68. Ellestad LE, Carré W, Muchow M, Jenkins SA, Wang X, Cogburn LA, Porter TE. Gene expression profiling during cellular differentiation in the embryonic pituitary gland using cDNA microarrays. Physiol Genomics. 2006;25:414–25.

69. Jenkins SA, Ellestad LE, Mukherjee M, Narayana J, Cogburn LA, Porter TE. Glucocorticoid-induced changes in gene expression in embryonic anterior pituitary cells. Physiol Genomics. 2013;45:422–33.

70. Porter TE, Ellestad LE. Gene expression profiling in the developing neuroendocrine system of the chick. In: Dawson A, Sharp PJ, editors. Functional avian endocrinology. New Dehli: Narosa Publishing House; 2005. p. 45–56.

71. Phillips J, Eberwine JH. Antisense RNA amplification: a linear amplification method for analyzing the mRNA population from single living cells. Methods. 1996;10:283–8.

72. Carré W, Wang X, Porter TE, Nys Y, Tang J, Bernberg E, Morgan R, Burnside J, Aggrey SE, Simon J, et al. Chicken genomics resource: sequencing and annotation of 35,407 ESTs from single and multiple tissue cDNA libraries and CAP3 assembly of a chicken gene index. Physiol Genomics. 2006;25:514–24.

73. Simon R, Radmacher MD, Dobbin K. Design of studies using DNA microarrays. Genet Epidemiol. 2002;23:21–36.

74. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N. Braisted J, Klapa M, Currier T, Thiagarajan M, et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003;34:374–8.

75. McCarthy FM, Want N, Magee GB, Nanduri B, Lawrence ML, Camon EB, Barrell DG, Hill DP, Dolan ME, Williams WP, et al. AgBase: a functional genomics resource for agriculture. BMC Genomics. 2006;7:229.

76. Reich M, Ohm K, Angelo M, Tamayo P. Mesirov JP. GeneCluster 2.0: an advanced toolset for bioarray analysis. Bioinformatics. 2004;20:1797–8.

77. Ellestad LE, Saliba J, Porter TE. Ontogenic characterization of gene expression in the developing neuroendocrine system of the chick. Gen Comp Endocrinol. 2011;171:82–93.

78. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25:402–8.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_23.png" /></figure>