Skip to main content

Identification of candidate genes and regulatory factors related to growth rate through hypothalamus transcriptome analyses in broiler chickens

Abstract

Background

Intensive selection for growth rate (GR) in broiler chickens carries negative after-effects, such as aberrations in skeletal development and the immune system, heart failure, and deterioration of meat quality. In Poland, fast-growing chicken populations are highly non-uniform in term of growth rate, which is highly unprofitable for poultry producers. Therefore, the identification of genetic markers for boiler GR that could support the selection process is needed. The hypothalamus is strongly associated with growth regulation by inducing important pituitary hormones. Therefore, the present study used this tissue to pinpoint genes involved in chicken growth control.

Results

The experiment included male broilers of Ross 308 strain in two developmental stages, after 3rd and 6th week of age, which were maintained in the same housing and feeding conditions. The obtained results show for the overexpression of genes related to orexigenic molecules, such as neuropeptide Y (NPY), aldehyde dehydrogenase 1 family, member A1 (ALDH1A1), galanin (GAL), and pro-melanin concentrating hormone (PMCH) in low GR cockerels.

Conclusion

The results reveal strong associations between satiety centre and the growth process. The present study delivers new insights into hypothalamic regulation in broiler chickens and narrows the area for the searching of genetic markers for GR.

Graphical abstract

Background

The broiler chicken industry has grown over its 60-year history. The Global Information and Early Warning System (GIEWS) reported that the world poultry meat output was estimated at 121.6 million tons in 2018 with an upward tendency [1]. The increasing demand for poultry was an engine to bred highly efficient chicken populations with significantly improved productive values [2]. The conducted selection aimed to improve the growth rate, feed efficiency and carcass meat content [3]. The market broiler mass after sixth weeks from hatching has increased 4-fold, and pectoralis muscle nearly 2-fold over 25 years [4]. Unfortunately, the rate of genetic changes resulted in the alternations in bird physiology [5]. In turn, intensive selection led to negative after-effects, such as aberrations in skeletal development [6] and the immune system [7], heart failure [8] and deterioration of meat quality [9] what reflects elevated attendance of meat defects in fast-growing populations. The most frequent poultry meat defects are deep pectoralis myopathies (DPM) [10]; white striping; Pale, Soft and Exudative (PSE) meat [11], wooden breast [12]; and intramuscular connective tissue, termed spaghetti meat [13].

Ross 308 is the world’s most-popular commercial broiler hybrid. Ross 308 chickens are characterised by rapid growth and high muscle mass developed at an early age, providing the possibility of early slaughtering. Usually, they are slaughtered at sixth weeks from hatching, which is associated with the attainment of an optimum balance between growth rate, feed efficiency, high muscle mass, high performance and feather maturity. Regardless of the flock age, most of Ross 308 populations show low uniformity of growth rate (GR), where GR is estimated as average gain in body weight (BW) for a defined period of rearing. The differences are already visible in the first weeks of life, and after six-weeks, they can even reach 1.5 kg in the BW in one broiler population, which is unprofitable for poultry producers, butchers, and the future processing. A similar observation was also described in different broiler strains [14, 15].

On the other hand, the hypothalamus plays an overarching role in the regulation of primary metabolic processes also in bird organisms. It controls body temperature, hunger, sleep, growth and circadian rhythms, by stimulating the pituitary which secretes essential hormones controlling other glands/organs. As an example, the hypothalamus releases somatoliberin (GHRH) that induces pituitary somatotropin (GH) secretion - most crucial growth factor. Therefore, the hypothalamic molecular activity is regulated by numerous exogenous stimuli [16]. Thus, the newest high-throughput molecular genetic tools are used for the determination of these processes. Kamineni [17] in his doctoral dissertation used RNA-sequencing based on next-generation sequencing (NGS) technology to show that heat stress in broiler chickens increases hypothalamic gene expression of rate-limiting enzymes, such as sterol regulatory element-binding transcription factor 1 (SREBF1) and lipoprotein lipase (LPL). In turn, Han et al. [18] using analysis of miRNAom suggested that miRNAs play an essential role in the hypothalamus for timing the rapid development of chicken gonads. While Piórkowska et al. [19] used the RNA-seq method to pinpoint the hypothalamic gene expression changes during postnatal development in broiler chickens.

The present study aimed to identify whether different growth rate in broiler chickens found the reflection in hypothalamic genes’ expression and the activation of the essential molecular processes. This approach enables to narrow the search area of potential genetic markers associated with the determination of growth in broilers. The present experiment used transcriptome analysis based on next-generation sequencing technology.

Results

Cockerel characteristics

During the experiment, fattening and slaughter traits for all chickens were recorded. The high growth rate (HGR) cockerels were 176 g and 685 g heavier than low growth rate (LGR) birds, after the 3rd and 6th week from hatching, respectively (Table 1). The growth rate did not affect fat level but influenced leg and breast muscle mass. Moreover, low growth rate (LGR) cockerels slaughtered in 3rd and 6th week of age took 55 g and 461.2 g less feed respectively and showed a higher feed conversion ratios (FCR), which were measured from 2nd week of age to the slaughter.

Table 1 Traits of investigated cockerels used in RNA sequencing and qPCR analyses

Analysis of differentially expressed genes

Over 90% of the processed reads were mapped to the chicken reference genome (GCA_000002315.5) included 70% to the exonic regions and 9% to the introns (see supplementary File S1). The clustering analysis excluded three samples from further processing, one of 3-week-old and two of the 6-week-old group (Fig. 1 b, d). The statistic clustering analysis showed that these samples were significantly distinct than the rest samples in the particular groups, which could be a consequence of unknown factor. The GR-dependent DEG analysis was performed using the DESeq2 method, which pinpointed for 43 shared DEGs between both age groups (Fig. 1 a, c, e). They encode proteins involved in luteinizing and follicle-stimulating hormone secretion (TBX3 and CGA), feeding behaviour (OXT, NPY, PMCH, and BSX), and defence response to the virus (OASL, MX1, IFITM3, IFIT5, and IFI6).

Fig. 1
figure1

The plots showed hypothalamic DEGs of Ross308 identified in the present analysis. The volcano plots showing the number of significant (red dots) and not-significant (black dots) (a) for 3-week-old-chickens and (c) 6-week old chickens. The heatmap Principal Component Analysis (PCA) showing sample clustering based on normalised reads for Ross308 broiler chickens characterised high (HGR) and low growth rate (LGR) (b) for 3-week old chickens and (d) for 6-week-old chickens. On heatmap were indicated seven samples per age groups, because two of them were pinpointed during clustering analysis that are distinct from particular groups, similarly as 6LGR sample that was also removed from further processing. Venn diagram presented DEGs (e) identified for 3-week-old, and 6-week-old chickens and 43 shared DEGs

The hypothalamic response related to changes in growth rate indicates 195 and 339 DEGs (fold-change ≥1.5, by adjusted P-value ≤0.05) in 3- and 6-week-old cockerels, respectively (goo.gl/4XX5mt). In younger LGR cockerels, increased expression of genes encoding proteins mainly involved in the organ development (BSX, MGP, MMP2, MYH11, EGR1, PCDH15, FST, LECT2, TBX3, and NPPC), thyroid hormone regulation process (CRYM, CGA, and CPQ), regulation of hormone levels (POMC, SLCO1C1, EGR1, GHR, CGA, TBX3, BLK, GAL, SLC30A8, BTK, ALDH1A1, DIO2, TTR, and NR5A1) and feeding behaviour (BSX, PMCH, NPW, NMU, and NKX2–1) (Table 2) were observed. In turn, DEGs of older LGR birds are involved in the activity of inflammatory response processes included a response to interferon-gamma (IRF9, CCL19, IFITM1, GBP1, IRF8, IRF7, and IRF1), T and B cell activation (RAC2, DMB2, CD3E, B2M, CD3D, LCK, CD3E, PTPRC, FOS, BLNK, BTK, and PTPN6), and antigen presentation: folding, assembly and peptide loading of class I (ERAP1, BFIV21, B2M, BF1, and TAP2), as well as, different signalling pathways regulating hormone levels (BLK, TTR, CGA, POMC, SLCO1C1, DIO2, EGR1, ALDH1A1, SLC30A8, NR5A1, TBX3, GHR, and GAL) and feeding behaviour (GAL, FOS, PMCH, BSX, OXT, and STRA6). Three genes showing the highest FC were ALB, MGAM, FGB and AVD, IFI6, CCL19 in 3- and 6-week-old LGR cockerel (goo.gl/4XX5mt), respectively.

Table 2 Functional annotation of hypothalamic DEGs in response to a variable growth rate of cockerels after 3rd and 6th weeks from hatching that was performed based on STRING and PANTHER tools

Several genes in the hypothalamus of HGR cockerels showed over 3-fold higher expression (AVP, ARAP1, and OXT for younger and GBX2, TCF7L2, GRM2, and MGAM for older cockerels).

Validation of RNA-seq results and qPCR analysis

Eight DEGs have been validated by the qPCR method and then compared with RNA-seq results using Pearson correlation. The obtained results are accessed by following the link goo.gl/4XX5mt. The lowest R coefficient was 0.82 for ALDH1A1 (P-value = 2.7E-09). Therefore, the qPCR analysis confirmed the RNA-seq results. The comparison between LGR and HGR cockerels groups indicates significantly increased expression of POMC, ALDH1A1, BSX, PMCH, and NPY genes in LGR birds. In turn, OXT gene expression shows the opposite tendency (Fig. 2).

Fig. 2
figure2

DEGs involved in the regulation of hormonal level. Relative transcript abundance of genes evaluated in the hypothalamus of broilers characterised high (high GR) and low growth rate (low GR) (a) and the relationship of genes coding protein involved in the hormonal regulation up-regulated in 3-week-old low growth rate (LGR) cockerels (b) and up-regulated in 6-week-old LGR cockerels (c). The efficiency of PCR reactions was estimated based on the standard curve method. The gene expression levels were calculated using the delta-delta CT method [20], and the significant differences in gene expression levels between HGR and LGR individuals within each age group were determined by ANOVA (Duncan’s post hoc test; SAS Enterprise v. 7.1 with default settings; SAS Institute, Cary, USA)

Discussion

The artificial selection in broiler chickens carried out over the recent decades focused on the increase of bird efficiency. The growth rate was the prime selection trait since the 1950s, and significant progress achieved due to the rise in pectoralis muscle yield and feed efficiency. Currently, broilers reach the slaughter weight faster due to more effective digestion and better energy utilisation during growth [21]. However, the growth rate of broiler chickens is not uniform since the slaughter weight can differ by up to 1.5 kg in one population. In turn, the flock uniformity in ‘live weight’ is a crucial measure of performance when optimising feeding programs of broiler population, as it strongly relates to the yield of processed meat. In many countries, including the USA and Australia, the significant financial losses related to numerous meat defects have been observed [14, 15]. The fast growth in broilers is associated with meat defects that reduce production profits. Therefore in the present study, the idea was born to indicate metabolic and signalling processes that are regulated in the hypothalamus under the influence of variable growth rate in broilers. The goal of the study was to pinpoint the important agents, which may enable control of this economically important chicken trait.

The advantage of the study was that the examined birds were selected from the same chicken population (eggs obtained from one source), maintained under the same feeding and housing conditions to minimalize other than genetic factor. Nonetheless, the differences in BW between the groups were high and constituted 176 and 685 g in the 3rd and 6th week from hatching, respectively. The DEG was performed, and results reveal interesting observation related to neurotransmitters regulation of the feeding depends on ontogenesis stages.

POMC gene as a precursor for adrenocorticotropic hormones (ACTH), β-endorphin, αmelanocyte-stimulating hormone (αMSH) and βmelanocyte-stimulating hormone (βMSH). It is released by the hypothalamic arcuate nucleus (ARC) and belongs to the anorexigenic molecules [22]. In avians, the regulation by this neurotransmitter is somewhat distinct from that in mammals, since Tachibana et al. [23] observed that broiler chickens artificially selected for their high growth rate were sensitive only to high αMSH doses, which contributed to the suppression of feed intake. In turn, Rice et al. [24] showed an increased POMC expression in the hypothalamus of selected high growth rate broilers only after postprandial insulin injection. Honda et al. [25] also suggested that increased feed consumption in broiler chickens could be associated with a lack of β-MSH anorexigenic effect, and an increased level of hypothalamic POMC is not associated with suppressing of food intake. Concerning an ambiguous role of POMC in birds, Boswell and Dunn [26] proposed that avian hypothalamic POMC plays a more significant role in the production of ACTH and β-endorphin than of α-MSH; thus, increased POMC expression promotes body mass loss. However, our results show elevated POMC mRNA level in the LGR cockerels that characterised by lowered feed intake. This observation supports the thesis of the anorexigenic function of POMC; suppressing hunger and promoting body mass loss. Moreover, in the younger cockerels, positive POMC neuron regulatory element— homeobox protein NK-2 homolog A (Nkx2–1) was found to be overexpressed, and Nkx2–1 ablation from POMC neurons leads to decreased POMC expression in adult males and mildly increased their body weight and adiposity [27, 28] (Fig. 3).

Fig. 3
figure3

The connection of GR-regulated genes that were associated with feeding behaviour (a) in 3-week-old and (b) in 6-week-old chickens. The colour saturation is changed along with the level of fold-change, down-regulated genes (orbs) in blue colour, up-regulated in red colour, and in grey colour that were not regulated in response to the growth rate. → inducing, activation; •-• interaction; inhibition

Another hypothalamic hormone-neuropeptide Y (NPY) that is associated with the nutritional state was identified as DEG depending on GR. The increased NPY (+ 50%) expression was noticed in the examined LGR cockerels, which was in line with previous studies [29, 30]. The present results show that increased NPY expression is associated with lowered feed intake, particularly in younger 3-week-old cockerels. In the same analysed cockerel group, NPY neuron regulatory element— brain-specific homeobox protein (BSX) was found to be overexpressed. The role of BSX was also examined in the loss-of-function study included ob/ob mice [31], where authors confirmed the requirement of BSX for physiological expression of NPY/AgRP and stimuli of hyperphagic response. Therefore the increased BSX expression level could be associated with its NPY regulatory function in avians, as well. Moreover, in 6LGR broiler increased FOS expression was observed that encodes c-Fos proto-oncogene. FOS is induced by NPY in insulin companion that action promotes food intake [32]. Further, it is highly probable that TRH mediates the NPY hyperphagia effect via c-FOS [33] (Fig. 3).

ALDH1A1, PMCH, and GAL genes were also differentially expressed in response to different GR and showed increased expression in LGR birds. In the mammal hypothalamus, ALDH1A1 participates in the metabolism of dopamine [34] and decreases its extracellular level. Mebel et al. [35] indicated that mesolimbic dopamine (DA) released in the ventral tegmental area of the brain of mammals implicated in the incentive, reinforcing and motivational aspects of food intake. Moreover, this signalling is sensitive to palatable food with high fat and sugar content that activate DA reward circuitry [36]. However, in chickens, intracerebroventricular injection of dopamine decreases food intake and induces hypophagia [37]. The PMCH encodes pre-melanin-concentrating hormone, which is associated with the regulation of feeding behaviour by decreasing energy expenditure and increasing food intake [38]. Sun et al. [39] found in chickens the relationship between a missense PMCH variant and shear force measured in the breast and leg muscles, but without effects on growth performance. In the present study, increased ALDH1A1 and PMCH expressions in the LGR chickens at both ages were observed, and between 3rd and 6th week, these differences showed a downward trend. Moreover, the LGR birds showed lowered growth rate and feed intake, what contradicts of the orexigenic function of both molecules in the avian hypothalamus. GAL encodes galanin that is believed to be an acute feeding behaviour stimulator in the mammals [40]. Fang et al. [41] described that endogenous galanin contributes to regulating glucose uptake because it decreases insulin resistance and improves its sensitivity. In birds, the role of galanin is not precisely clear since contradictory literature evidence can be found [23]. While, the present study found increased hypothalamic GAL expression in 6LGR broilers, in which low body mass was likely the consequence of low feed intake. The present results deny of orexigenic galanin action at least at this developmental stage in broilers, although Tachibana et al. [40] suggested that galanin function in the vertebrates is constant.

In 3-week-old HGR chickens, increased expression of AVP and OXT was observed. Both neurotransmitters are known as regulators of nutritional status in the mammals [42, 43]. The paraventricular nucleus of the mammalian hypothalamus (PVH) is the central localisation of melanocortin action. The melanocortin system components are crucial for energy homeostasis. In turn, arginine vasopressin encoded by AVP gene is expressed in the mammalian PVH, and is induced by anorectic melanocortin agonist melanotan II (MTII). Pei et al. [42] examined the effect of AVP on feed intake by using a specific knock-out (KO) mice in combination with a viral vector to manipulate PVH-AVP, and stated that PVH-AVP neurons are critical components for transducing the anorectic effect of melanocortin agonists. The increased AVP expression in the hypothalamus was observed only in 3-week-old HGR cockerels with increased feed intake. Thus, a doubtful role of AVP as a mediator during the transduction of anorectic MTII action at this developmental stage in broilers, or different signalling pipeline of MTII in this context. The OXT gene encodes oxytocin that is expressed in the PVH, which was confirmed in mammals. Numerous studies reported an anorexigenic effect of oxytocin in rats [44, 45]. Early studies suspected that OXT action during appetite suppression displays in homeostatic modulation of sodium balance [46]. In turn, the recent research suggests that oxytocin limits food intake rich in sucrose by suppressing the reward pathway [43]. The increased OXT expression in the HGR group that was characterised by increased feed intake suggests for a different function of OXT during nutritional regulation. Moreover, in older 6-week-old broilers, 2-fold increased SIM1 expression was observed. SIM1 encodes transcription factor single-minded 1 that was found to interact with oxytocin [47]. The authors suggested that SIM1-expressing PVN neurons regulate feeding in response to melanocortin signalling but not energy expenditure. Therefore in the mammals, SIM1 is believed to be an anorexigenic molecule. In addition, Michaud et al. [48] indicated that SIM1 upstream regulation of POU3F2 expression is related to differentiation of the specific neuroendocrine lineage of the PVH and secretion of OXT. The present study, as a first report, shows that increased SIM1 expression was positively correlated with body mass and feed intake in broilers what suggest the orexigenic effect of this molecule in avians.

The present study shows that numerous orexigenic and anorexigenic molecules examined previously in the mammals seems to play a different role in the regulation of birds feeding behaviour. Nevertheless, the activity of neuromedin U (NMU) and neuropeptide W (NPW) in cockerels seems to be in line with previous studies in the mammals, because in the hypothalamus of 3LGR birds increased expression of both neurotransmitters was observed. The anorexigenic action of NMU was examined previously by Wren et al. [49], who found that NMU injected into the PVN or acute nucleus immediately decreased food intake in rats. Moreover, NMU released in the hypothalamus is stimulated by leptin, which is the major anorexigenic molecule. It has been confirmed that in pigs, neuropeptide W increases the feed intake in the bright phase and inhibits it in the dark period [50]. Its function was not fully elucidated, and the authors suggested that NPW regulates feeding behaviour in the hypothalamus by accompanying other feeding-regulating neurons. The present results suggest that NMU and NPW actions in feeding behaviour suppression in birds are dependent on age.

Conclusion

In the present study, low growth rate chickens with lowered feed intake showed increased expression of genes that are known as orexigenic molecules in the mammals. Our observation indicates for a little distinct feeding behaviour regulation at hormone gene expression levels in avian than it was described by numerous authors in the mammals. Moreover, we presented several candidate genes associated with the growth rate in broiler chickens such as POMC, NMU, NPW, PMCH, GAL, and FOS that are differentially expressed in the hypothalamus in response to variable growth. The findings narrow the search area for genetic markers, which could be used in the selection procedure, and to promote population homogeneity that determines poultry producers’ profits. The present findings enrich our knowledge in terms of the regulation of nutrition processes in the broiler hypothalamus showing interesting avian hormonal controls of growth and feed intake.

Methods

Birds used in the experiment

All conducted chicken trial was approved by the Approving Experiment Committee of National Research Institute of Animal Production (Kraków, Poland), and with the national authority according to the Polish Act on the Protection of Animals Used for Scientific or Educational Purposes of 15 January 2015 (which implements Directive 2010/63/EU of the European Parliament and the Council of 22 September 2010 on the protection of animals used for scientific purposes). Moreover, all procedures were following the guidelines and regulations of the Local Krakow Ethics Committee for Experiments with Animals.

The study included 36 Ross 308 cockerels (300 in total) that were hatched in a commercial poultry hatchery. Eggs were obtained from the Ross 308 parent stock farm located in Wolbrom (Poland). Chicks were delivered to the experimental farm of the National Research Institute of Animal Production located in Aleksandrowice (Poland). Birds were weighed on day 1st and day 14th. Based on broiler mass after 2nd week, they were classified into lighter and heavier groups (n = 36) that were kept in separate pens (n = 9) on the deep litter with electronically-controlled environmental conditions (temperature, air humidity, lighting regime). During test fattening, traits were measured, such as body weight gain, feed intake and feed conversion ratio (FCR). The cockerels were standardly fed complete starter, grower and finisher diets in the period between days 1–21, 22–35 and 36–45, respectively. In the particular nutrition periods, the diet contained 22, 20.5 and 20.5% crude protein, respectively; and 2990, 3130 and 3130 kcal/kg metabolizable energy, respectively. Feed and water were available ad libitum. Ross 308 chickens show intensive growth on the day 21 from hatching, and in two following weeks, they double their BW and then GR gradually decreases. The present study confirms the differences in the hypothalamic activity at the transcriptome level in HGR and LGR Ross 308 cockerels at these two crucial developmental stages; after week 3 and 6. Eighteen birds were kept until day 22 from hatching (after 3rd week), with 9 in the low growth rate (3LGR) and 9 in the high growth rate (3HGR) groups, and they reached a mean BW of 1.2 kg. The next 18 birds were kept until day 45 (6th week) in the 6LGR and 6HGR groups, and they reached a mean BW of 3.5 kg. Before the slaughter, birds were subjected to a 6-h feed withdrawal but had constant access to water. The birds were euthanized by decapitation after electrical stunning (150 mA, frequency of 200 Hz for 4 s). Then, the whole chicken hypothalamus was collected and stored as was described by Piórkowska et al. [19]. Carcass and growth traits were measured following the method described by Połtowicz et al. [51].

The statistical analysis between high and low growth rate chickens was estimated separately for 3rd and 6th-week old chickens. The data was analysed by the Student’s t-test using Statistica 13.1 (StatSoft Inc.). A significance level of P ≤ 0.05 was considered statistically significant.

cDNA library construction

RNA from hypothalamus was isolated, as was described by Piórkowska et al. [19]. The quality and concentration of RNA were measured using a TapeStation2200 (Agilent, Palo Alto, USA), and obtained RIN values were in the range between 7.3 and 8.4. RNA sequencing was performed for 16 randomly selected hypothalamus samples presenting HGR (n = 4) and LGR (n = 4) for both ages (Fig. 4). The cDNA libraries were prepared by TruSeq RNA Sample Preparation Kit v2 (Illumina, San Diego, CA, USA) according to the protocol and indexed with individual adaptors (File S1). The sequencing was performed on a HiScanSQ System (Illumina, San Diego, CA, USA) in 101 bp single end cycles using a TruSeq Kit v3-HS chemistry described by Piórkowska et al. [52].

Fig. 4
figure4

Bird selecting pathway for detecting the differentially expressed genes

Alignment

Raw reads were processed using the FastQC, Flexbar and CASAVA version 1.8.2 programs (default parameters). The processed reads were aligned to the GCA_000002315.5 genome assembly (Ensembl Genome Browser). The alignment was done by RSEM, v. 1.3.0, and SAMStat software, and to generate statistics RNA-SeQC v1.1.8 (default settings) was used. The RNA sequencing data have been submitted to the Gene Expression Omnibus NCBI (accession number GSE104255).

Differentially expressed gene analysis

Differentially expressed gene analysis was carried out using DESeq2 [53]. False Discovery Rate (FDR) was used as corrected p-values (≤ 0.05) due to multiple analysis. The threshold for gene expression fold change (FC) (comparison between HGR and LGR cockerels) was set at ≥1.5 level. The functional gene analysis was performed based on G. gallus reference using STRING v.10.5 and PANTHER™ Classification System v.14.0 tools (with default settings and established FDR at < 0.05 level).

RNA-seq result validation

Eighth genes were validated by qPCR using RPL4 and SDHA as endogenous controls [54]. qPCR reagents are presented in the supplementary material (File S2). Pearson’s correlation (SAS Enterprise v. 7.1) was used to compare RNA-seq and qPCR results. The genes for validation were chosen based on their relationship to the hormonal regulation.

The estimation of transcript abundance was performed for 36 cockerels at both ages (3- and 6-week-old chickens) and two GR groups (Fig. 4). The standard curve method was used to estimated qPCR. Relative transcript quantity (RQ) was calculated using delta-delta CT method [20], and significant results were determined by ANOVA (Duncan’s post hoc test; SAS Enterprise v. 7.1 with default settings; SAS Institute, Cary, USA).

The scheme of the whole experiment was shown in graphical abstract.

Availability of data and materials

The data generated during the study have been deposited in the GEO (NCBI) under accession number GSE104255 available by link.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104255.

Moreover, RNA-sequencing cleaned reads were aligned to the reference chicken genome GCA_000002315.5 available by link https://www.ensembl.org/Gallus_gallus/Info/Index. Functional analysis was perofroed using STRING: functional protein association networks v. 10.5 (https://string-db.org/) and PANTHER™ 14.0 (http://www.pantherdb.org/).

Abbreviations

GR:

Growth rate

NPY:

Neuropeptide Y

ALDH1A1:

Aldehyde dehydrogenase 1 family, member A1

GAL:

Galanin

PMCH:

Pro-Melanin concentrating hormone

GIEWS:

Global Information and early warning system

DPM:

Deep pectoralis myopathies

PSE:

Pale, soft and exudative meat

SREBF1:

Sterol regulatory element-binding transcription factor 1

LPL:

Lipoprotein lipase

HGR:

High growth rate

LGR:

Low growth rate

DEGs:

Differentially expressed genes

FCR:

Feed conversion ratio

ACTH:

Adrenocorticotropic hormones

αMSH:

αmelanocyte-stimulating hormone

βMSH:

βmelanocyte-stimulating hormone

ARC:

Arcuate nucleus

BSX:

Homeobox protein

POMC:

Pro-opiomelanocortin

PVH:

Mammalian hypothalamus

AVP:

Arginine vasopressin

MTII:

Melanocortin agonist melanotan II

OXT:

Oxytocin

SIM1:

Transcription factor single-minded 1

VP:

Vasopressin

TRH:

Thyrotropin-releasing hormone

CRH:

Corticotropin-releasing hormone

GHIH:

Somatostatin

NMU:

Neuromedin U

NPW:

Neuropeptide W

FDR:

False discovery rate

References

  1. 1.

    Meat Market Review. Overview of global meat market developments in 2018. 2018.

  2. 2.

    Buzala M, Janicki B. Review Review: effects of different growth rates in broiler breeder and layer hens on some productive traits. Poult Sci. 2016;95:2151–9.

    CAS  Article  Google Scholar 

  3. 3.

    Havenstein G, Ferket P, Qureshi M. Carcass composition and yield of 1957 versus 2001 broilers when fed representative 1957 and 2001 broiler diets. Poult Sci. 2003;82:1509–18.

    CAS  Article  Google Scholar 

  4. 4.

    Zuidhof MJ, Schneider BL, Carney VL, Korver DR, Robinson FE. Growth, efficiency, and yield of commercial broilers from 1957, 1978, and 2005. Poult Sci. 2014;93:2970–82.

    CAS  Article  Google Scholar 

  5. 5.

    McKay JC, Barton NF, Koerhuis ANM, McAdam J. The challenge of genetic change in the broiler chicken. BSAP Occas Publ. 2000;27:1–7.

    Article  Google Scholar 

  6. 6.

    Julian RJ. Production and growth related disorders and other metabolic diseases of poultry - a review. Vet J. 2005;169:350–69.

  7. 7.

    Koenen ME, Boonstra-Blom AG, SHM J. Immunological differences between layer- and broiler-type chickens. Vet Immunol Immunopathol. 2002;89:47–56.

  8. 8.

    Olkowski AA. Pathophysiology of heart failure in broiler chickens: structural, biochemical, and molecular characteristics. Poult Sci. 2007;86:999–1005.

    CAS  Article  Google Scholar 

  9. 9.

    Bailey RA, Watson KA, Bilgili SF, Avendano S. The genetic basis of pectoralis major myopathies in modern broiler chicken lines. Poult Sci. 2015;94:2870–9.

    CAS  Article  Google Scholar 

  10. 10.

    Bianchi M, Petracci M, Franchini A, Cavani C. The occurrence of deep pectoral myopathy in roaster chickens. Poult Sci. 2006;85:1843–6.

    CAS  Article  Google Scholar 

  11. 11.

    Owens CM, Alvarado CZ, Sams AR. Research developments in pale, soft, and exudative Turkey meat in North America. Poult Sci. 2009;88:1513–7.

    CAS  Article  Google Scholar 

  12. 12.

    Abasht B, Mutryn MF, Michalek RD, Lee WR. Oxidative stress and metabolic perturbations in wooden breast disorder in chickens. PLoS One. 2016;11:e0153750.

  13. 13.

    Purslow PP. Intramuscular connective tissue and its role in meat quality. Meat Sci. 2005;70:435–47.

    CAS  Article  Google Scholar 

  14. 14.

    Madsen T, Feedstuffs JP. Broiler flock uniformity. Feedstuffs. 2010;82:12–3.

    Google Scholar 

  15. 15.

    Hughes R, Heberle N, Barekatin R, Edwards N, Hynd P. Flock uniformity-is it important and how is it assessed. Sydney, Australia: Conference: Australian Poultry Science SymposiumAt; 2017. p. 28.

    Google Scholar 

  16. 16.

    Pearson CA, Placzek M. Development of the medial hypothalamus: forming a functional hypothalamic-Neurohypophyseal Interface. Curr Top Dev Biol. 2013;106:49–88.

  17. 17.

    Kamineni LP. Transcriptomic analysis of hypothalamic responses to heat stress in modern and legacy chicken lines. 2015.

  18. 18.

    Han W, Zou J, Wang K, Su Y, Zhu Y, Song C, et al. High-throughput sequencing reveals hypothalamic MicroRNAs as novel partners involved in timing the rapid development of chicken (Gallus gallus) gonads. PLoS One. 2015;10.

  19. 19.

    Piórkowska K, Żukowski K, Połtowicz K, Nowak J, Wojtysiak D, Derebecka N, et al. Transcriptomic changes in broiler chicken hypothalamus during growth and development. Int J Genomics. 2018;6049469:1–10.

  20. 20.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:45.

    Article  Google Scholar 

  21. 21.

    Tallentire CW, Leinonen I, Kyriazakis I. Breeding for efficiency in the broiler chicken: a review. Agron Sustain Dev. 2016;36:66.

    Article  Google Scholar 

  22. 22.

    Zhan C, Zhou J, Feng Q, Zhang J-E, Lin S, Bao J, et al. Behavioral/cognitive acute and long-term suppression of feeding behavior by POMC neurons in the brainstem and hypothalamus, Respectively. J Neurosci. 2013;33:3624–32.

    CAS  Article  Google Scholar 

  23. 23.

    Tachibana T, Mori M, Khan MSI, Ueda H, Sugahara K, Hiramatsu K. Central administration of galanin stimulates feeding behavior in chicks. Comp Biochem Physiol A Mol Integr Physiol. 2008;151:637–40.

    Article  Google Scholar 

  24. 24.

    Rice BB, Zhang W, Bai S, Siegel PB, Cline MA, Gilbert ER. Insulin-induced hypoglycemia associations with gene expression changes in liver and hypothalamus of chickens from lines selected for low or high body weight. Gen Comp Endocrinol. 2014;208:1–4.

    CAS  Article  Google Scholar 

  25. 25.

    Honda K, Saneyasu T, Hasegawa S, Peptides HK-, 2012 undefined. A comparative study of the central effects of melanocortin peptides on food intake in broiler and layer chicks. Peptides. 2012;37:13–7.

  26. 26.

    Boswell T, Dunn IC. Regulation of agouti-related protein and pro-opiomelanocortin gene expression in the avian arcuate nucleus. Front Endocrinol (Lausanne). 2017;8:75.

  27. 27.

    Yee CL, Wang Y, Anderson S, Ekker M, Rubenstein JLR. Arcuate nucleus expression of NKX2.1 and DLX and lineages expressing these transcription factors in neuropeptide Y+, proopiomelanocortin+, and tyrosine hydroxylase+ neurons in neonatal and adult mice. J Comp Neurol. 2009;517:37–50.

    CAS  Article  Google Scholar 

  28. 28.

    Orquera DP, Belén Tavella M, De Souza FSJ, Nasif S, Low XJ, Rubinstein XM. The Homeodomain Transcription Factor NKX2.1 Is Essential for the Early Specification of Melanocortin Neuron Identity and Activates Pomc Expression in the Developing Hypothalamus. J Neurosci. 2019;39:4023–35.

  29. 29.

    Ka S, Lindberg J, Strömstedt L, Fitzsimmons C, Lindqvist N, Lundeberg J, et al. Extremely different behaviours in high and low body weight lines of chicken are associated with differential expression of genes involved in neuronal plasticity. J Neuroendocrinol. 2009;21:208–16.

    CAS  Article  Google Scholar 

  30. 30.

    Yuan L, Ni Y, Barth S, Wang Y, Grossmann R, Zhao R. Layer and broiler chicks exhibit similar hypothalamic expression of orexigenic neuropeptides but distinct expression of genes related to energy homeostasis and obesity. Brain Res. 2009;1273:18–28.

    CAS  Article  Google Scholar 

  31. 31.

    Sakkou M, Wiedmer P, Anlag K, Hamm A, Seuntjens E, Ettwiller L, et al. A role for brain-specific homeobox factor Bsx in the control of hyperphagia and locomotory behavior. Cell Metab. 2007;5:450–63.

    CAS  Article  Google Scholar 

  32. 32.

    Boelen A, Wiersinga WM, Fliers E. Fasting-induced changes in the hypothalamus-pituitary-thyroid axis. Thyroid. 2008;18:123–9.

    CAS  Article  Google Scholar 

  33. 33.

    Fan S, Dakshinamoorthy J, Kim ER, Xu Y, Huang C, Tong Q. An indirect action contributes to C-Fos induction in Paraventricular hypothalamic nucleus by neuropeptide y. Sci Rep. 2016;6:19980.

    CAS  Article  Google Scholar 

  34. 34.

    Anderson DW, Schray RC, Duester G, Schneider JS. Functional significance of aldehyde dehydrogenase ALDH1A1 to the nigrostriatal dopamine system. Brain Res. 2011;1408:81–7.

    CAS  Article  Google Scholar 

  35. 35.

    Mebel DM, Wong JCY, Dong YJ, Borgland SL. Insulin in the ventral tegmental area reduces hedonic feeding and suppresses dopamine concentration via increased reuptake. Eur J Neurosci. 2012;36:2336–46.

    Article  Google Scholar 

  36. 36.

    Baik JH. Dopamine signaling in food addiction: role of dopamine D2 receptors. BMB Rep. 2013;46:519–26.

    CAS  Article  Google Scholar 

  37. 37.

    Zendehdel M, Hassanpour S. Central regulation of food intake in mammals and birds: a review. Neurotransmitter. 2014;1:1–7.

  38. 38.

    Date Y, Nakazato M, Matsukura S. A role for orexins and melanin-concentrating hormone in the central regulation of feeding behavior. Nihon rinsho Jpn J Clin Med. 2001;59:427–30.

  39. 39.

    Sun G, Li M, Li H, Tian Y, Chen Q, Bai Y, et al. Molecular cloning and SNP association analysis of chicken PMCH gene. Mol Biol Rep. 2013;40:5049–55.

    CAS  Article  Google Scholar 

  40. 40.

    Tachibana T, Tsutsui K. Neuropeptide control of feeding behavior in birds and its difference with mammals. Front Neurosci. 2016;10:485.

  41. 41.

    Fang PH, Yu M, Ma YP, Li J, Sui YM, Shi MY. Central nervous system regulation of food intake and energy expenditure: role of galanin-mediated feeding behavior. Neurosci Bull. 2011;27:407–12.

    CAS  Article  Google Scholar 

  42. 42.

    Pei H, Sutton AK, Burnett KH, Fuller PM, Olson DP. AVP neurons in the paraventricular nucleus of the hypothalamus regulate feeding. Mol Metab. 2014;3:209–15.

    CAS  Article  Google Scholar 

  43. 43.

    Sabatier N, Leng G, Menzies J. Oxytocin, feeding, and satiety. Front Endocrinol (Lausanne). 2013;4:35.

  44. 44.

    Arletti R, Benelli A, Bertolini A. Oxytocin inhibits food and fluid intake in rats. Physiol Behav. 1990;48:825–30.

    CAS  Article  Google Scholar 

  45. 45.

    Olson BR, Drutarosky MD, Chow MS, Hruby VJ, Stricker EM, Verbalis JG. Oxytocin and an oxytocin agonist administered centrally decrease food intake in rats. Peptides. 1991;12:113–8.

    CAS  Article  Google Scholar 

  46. 46.

    Flanagan LM, Blackburn RE, Verbalis JG, Stricker EM. Hypertonic NaCl inhibits gastric motility and food intake in rats with lesions in the rostral AV3V region. Am J Physiol Integr Comp Physiol. 1992;263:R9–14.

    CAS  Article  Google Scholar 

  47. 47.

    Kublaoui BM, Holder JL, Gemelli T, Zinn AR. Sim1 haploinsufficiency impairs melanocortin-mediated anorexia and activation of paraventricular nucleus neurons. Mol Endocrinol. 2006;20:2483–92.

    CAS  Article  Google Scholar 

  48. 48.

    Michaud JL, Rosenquist T, May NR, Fan CM. Development of neuroendocrine lineages requires the bHLH-PAS transcription factor SIM1. Genes Dev. 1998;12:3264–75.

    CAS  Article  Google Scholar 

  49. 49.

    Wren AM, Small CJ, Abbott CR, Jethwa PH, Kennedy AR, Murphy KG, et al. Hypothalamic actions of neuromedin U. Endocrinology. 2002;143:4227–34.

    CAS  Article  Google Scholar 

  50. 50.

    Takenoya F, Haruaki K, Date Y, Nakazato M, Shioda S. Feeding regulation in the brain: involvement of neuropeptide W. In: Transmitters and Modulators in Health and Disease. Japan: Springer; 2009. p. 31–9.

    Google Scholar 

  51. 51.

    Połtowicz K, Nowak J, Wojtysiak D. Effect of feed restriction on performance, carcass composition and physicochemical properties of the m. pectoralis superficialis of broiler chickens. Ann Anim Sci. 2015;15:1019–29.

    Article  Google Scholar 

  52. 52.

    Piórkowska K, Żukowski K, Nowak J, Połtowicz K, Ropka-Molik K, Gurgul A. Genome-wide RNA-Seq analysis of breast muscles of two broiler chicken groups differing in shear force. Anim Genet. 2016;47:68–80.

  53. 53.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15.

  54. 54.

    Piórkowska K, Nowak J, Połtowicz K. The normalisation of CAPN gene expression in M. pectoralis superficialis in broiler lines differing in growth rate and their relationship to breast muscle tenderness. Br Poult Sci. 2015;56:452–8.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The study was supported by the statutory activity of National Research Institute of Animal Production (no. 01–013.1) and the National Multidisciplinary Laboratory of Functional Nanomaterials NanoFun no. POIG.02.02.00–00-025/09 (Innovative Economy Operational Programme, Priority Axis 2: R&D Infrastructure, Action 2.2: Support of Formation of Common Research Infrastructure of Scientific Units). Adam Mickiewicz University, Faculty of Biology, is a member of KNOW Poznan RNA Centre (no. 01/KNOW2/2014).

Author information

Affiliations

Authors

Contributions

Conceptualization, K.P. 1; investigation, K.P. 1, J.N., N.D., D.W. methodology, K.P.1, J.N., D.W., K.P.2; resources, J.N., K.P.2; validation, K.P.1; visualization, K.P.1; K.Ż.; software, K.Ż., project administration, J.N., writing—original draft preparation, K.P.1, J.N., K.P.2; writing—review and editing, K.P.1, J.N., K.P.2, N.D., J.W., K.R.-M. supervision, K.R.-M., K.P.2; funding acquisition, K.P.2, J.W. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Katarzyna Połtowicz.

Ethics declarations

Ethics approval and consent to participate

All conducted chicken trial was approved by the Approving Experiment Committee of National Research Institute of Animal Production (Kraków, Poland) and the national authority according to the Polish Act on the Protection of Animals Used for Scientific or Educational Purposes of 15 January 2015 (which implements Directive 2010/63/EU of the European Parliament and the Council of 22 September 2010 on the protection of animals used for scientific purposes). Moreover, all procedures were following the guidelines and regulations of the Local Krakow Ethics Committee for Experiments with Animals.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: File S1.

Overall statistics and read annotations obtained for each library.

Additional file 2: File S2.

Primers and probes used in the qPCR analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Piórkowska, K., Żukowski, K., Połtowicz, K. et al. Identification of candidate genes and regulatory factors related to growth rate through hypothalamus transcriptome analyses in broiler chickens. BMC Genomics 21, 509 (2020). https://doi.org/10.1186/s12864-020-06884-5

Download citation

Keywords

  • Growth rate
  • Broilers
  • RNA-seq
  • Hypothalamus response