Correlated responses to long-term divergent selection for 8-week body weight in female White Plymouth Rock chickens: Sexual maturity

Molecular mechanisms of hypothalamic-pituitary-ovarian/thyroid axis regulating age at first egg in geese

Age at first egg (AFE) has consistently garnered interest as a crucial reproductive indicator within poultry production. Previous studies have elucidated the involvement of the hypothalamic-pituitary-ovarian (HPO) and hypothalamic-pituitary-thyroid (HPT) axes in regulating poultry sexual maturity. Concurrently, there was evidence suggesting a potential co-regulatory relationship between these 2 axes. However, as of now, no comprehensive exploration of the key pathways and genes responsible for the crosstalk between the HPO and HPT axes in the regulation of AFE has been reported. In this study, we conducted a comparative analysis of morphological differences and performed transcriptomic analysis on the hypothalamus, pituitary, thyroid, and ovarian stroma between normal laying group (NG) and abnormal laying group (AG). Morphological results showed that the thyroid index difference (D-) value (thyroid index D-value=right thyroid index-left thyroid index) was significantly (P < 0.05) lower in the NG than in the AG, while the ovarian index was significantly (P < 0.01) higher in the NG than in the AG. Furthermore, between NG and AG, we identified 99, 415, 167, and 1182 differentially expressed genes (DEGs) in the hypothalamus, pituitary, thyroid, and ovarian stroma, respectively. Gene ontology (GO) analysis highlighted that DEGs from 4 tissues were predominantly enriched in the "biological processes" category. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that 16, 14, 3, and 26 KEGG pathways were significantly enriched (P < 0.05) in the hypothalamus, pituitary, thyroid, and ovarian stroma. The MAPK signaling pathway emerged as the sole enriched pathway across all 4 tissues. Employing an integrated analysis of the protein-protein interaction (PPI) network and correlation analysis, we found GREB1 emerged as a pivotal component within the HPO axis to regulate estrogen-related signaling in the HPT axis, meanwhile, the HPT axis influenced ovarian development by regulating thyroid hormone-related signaling mainly through OPN5. Then, 10 potential candidate genes were identified, namely IGF1, JUN, ERBB4, KDR, PGF, FGFR1, GREB1, OPN5, DIO3, and THRB. These findings establish a foundation for elucidating the physiological and genetic mechanisms by which the HPO and HPT axes co-regulate goose AFE.

Physiological role of dietary energy in the sexual maturity: clues of body size, gonad development, and serum biochemical parameters of Chinese indigenous chicken

Sexual maturity is a crucial factor in the formation and development of poultry reproductive capacity. The nutritional status has been confirmed to play an important role in the regulation of sexual maturity. To investigate the effect of dietary energy levels on sexual maturity in chicken, diets with 3 energy levels (group L: 2,573 kcal/kg, group C: 2,836 kcal/kg, group H: 3,122 kcal/kg) were implemented to feed Guangyuan Gray chickens. During this trial, body weight, body size, organ development, sexual maturity, reproductive performance and blood biochemical parameters were monitored. The earlier sexual maturity was observed in group H, as well as a heavier first egg weight, larger interpubic distance and higher total cholesterol (T-CHO) content at sexual maturity. The dietary energy levels had no significant effect on body weight at first egg and egg production at 300 d of age. Although dietary energy levels had a significant effect on body weight, comb length, tibia length and girth, abdominal fat weight, oviduct weight and length, T-CHO, triglyceride (TG) content and estradiol (E2) level during the rearing period. No significant difference of gonadotropin releasing hormone (GnRH) and luteinizing hormone (LH) level among 3 groups was observed during the trial. The dietary energy levels had effects on mRNA expression of GnRH, estrogen receptor 1 (ESR1), estrogen receptor 2 (ESR2) in hypothalamus, gonadotropin inhibitory hormone receptor (GnIHR) in pituitary and luteinizing hormone receptor (LHR), ESR2 in ovary. The GnIHR/GnRHR ratio in pituitary was higher before sexual maturity and decreased at sexual maturity. The results of correlations analysis found that all the body size, carcass traits, serum biochemical parameters negatively correlated with age at first egg except for interpubic distance and serum blood glucose content. Collectively, dietary energy levels had effects on sexual maturity of chicken, which may be achieved by affecting body weight, gonad development, endocrine and the mRNA expression of genes related to hypothalamus-pituitary-gonad axis. These results further set our understanding of how dietary energy regulates sexual maturity.

Population dynamics of a long-term selection experiment in White Plymouth Rock chickens selected for low or high body weight

The population dynamics of 2 lines of chickens from a long-term (59 generations) selection experiment were assessed based on pedigree data. These lines were propagated from phenotypic selection for low and high 8-wk BW in White Plymouth Rock chickens. Our objective was to determine whether the 2 lines maintained similar population structures over the selection horizon to allow meaningful comparisons of their performance data. A complete pedigree of 31,909 individuals, consisting of 102 founders, 1,064 from the parental generation, and 16,245 low weight (LWS) and 14,498 high weight (HWS) select chickens, was available. Inbreeding (F) and average relatedness (AR) coefficients were computed. Average F per generation and AR coefficients were 1.3 (SD 0.8) % and 0.53 (SD 0.001) for LWS, and 1.5 (SD 1.1) % and 0.66 (SD 0.001) for HWS. Mean F for the entire pedigree was 0.26 (0.16) and 0.33 (0.19), and maximum F was 0.64 and 0.63, in LWS and HWS, respectively. Based on Wright's fixation index, at generation 59, substantial genetic differences were established between lines. The effective population size was 39 in LWS and 33 in HWS. The effective number of founders was 17 and 15, effective number of ancestors were 12 and 8, and genome equivalents were 2.5 and 1.9 in LWS and HWS, respectively. About 30 founders explained the marginal contribution to both lines. By generation 59, only 7 male and 6 female founders contributed to both lines. Moderately high levels of inbreeding and low effective population sizes were inevitable, as this was a closed population. However, effects on the fitness of the population were expected to be less substantial because founders were a combination of 7 lines. The effective numbers of founders and ancestors were relatively low compared to the actual number of founders, as few ancestors contributed to descendants. Based on these evaluations, it can be inferred that LWS and HWS had similar population structures. Comparisons of selection responses in the 2 lines therefore should be reliable.

Intestinal Microbiota and Host Cooperate for Adaptation as a Hologenome

Multiomic analyses reported here involved two lines of chickens, from a common founder population, that had undergone long-term selection for high (HWS) or low (LWS) 56-day body weight. In these lines that differ by around 15-fold in body weight, we observed different compositions of intestinal microbiota in the holobionts and variation in DNA methylation, mRNA expression, and microRNA profiles in the ceca. Insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) was the most upregulated gene in HWS ceca with its expression likely affected by the upregulation of expression of gga-miR-2128 and a methylated region near its transcription start site (388 bp). Correlation analysis showed that IGF2BP1 expression was associated with an abundance of microbes, such as Lactobacillus and Methanocorpusculum. These findings suggest that IGF2BP1 was regulated in the hologenome in adapting to long-term artificial selection for body weight. Our study provides evidence that adaptation of the holobiont can occur in the microbiome as well as in the epigenetic profile of the host.

IMPORTANCE The hologenome concept has broadened our perspectives for studying host-microbe coevolution. The multiomic analyses reported here involved two lines of chickens, from a common founder population, that had undergone long-term selection for high (HWS) or low (LWS) 56-day body weight. In these lines that differ by around 15-fold in body weight, we observed different compositions of intestinal microbiota in the holobionts, and variation in DNA methylation, mRNA expression, and microRNA profiles in ceca. The insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) was the most upregulated gene in HWS ceca with its expression likely affected by a methylated region near its transcription start site and the upregulation of expression of gga-miR-2128. Correlation analysis also showed that IGF2BP1 expression was associated with the abundance of microbes, such as Lactobacillus and Methanocorpusculum. These findings suggest that IGF2BP1 was regulated in the hologenome in response to long-term artificial selection for body weight. Our study shows that the holobiont may adapt in both the microbiome and the host's epigenetic profile.

Effect of age at first egg on reproduction performance and characterization of the hypothalamo-pituitary-gonadal axis in chickens

As an important indicator of sexual maturity of hens, age at first egg (AFE) is significantly associated with reproduction performance. In this study, 400 hens were divided into 6 groups based on AFE to analyze the difference of reproduction performance, reproduction hormone levels and the characterization of the hypothalamo-pituitary-gonadal axis in chickens. The results showed that the egg production of the hens in the late-maturing groups was significantly lower than that of the ones in other groups and the precocious hens had a lower egg production rate. The hens in late-maturing group had a lower fertility rate, LH levels and shorter duration of peak of egg production (PEP), the precocious hens had lower PRL levels. In addition, the characterization of the hypothalamo-pituitary-gonadal axis showed that the individuals with normal AFE had higher GNRH, GNRHR, ESR1, KITLG, and CYP11A1 expression levels than late-maturing and precocious individuals, which indicated that the chickens with normal AFE advantages on reproduction regulation system.

Haplotype Purging After Relaxation of Selection in Lines of Chickens that Had Undergone Long-Term Selection for High and Low Body Weight

Bi-directional selection for increased and decreased 56-day body weights (BW56) has been applied to two lines of White Plymouth Rock chickens—the Virginia high (HWS) and low (LWS) body weight lines. Correlated responses have been observed, including negative effects on traits related to fitness. Here, we use high and low body weight as proxies for fitness. On a genome-wide level, relaxed lines (HWR, LWR) bred from HWS and LWS purged some genetic variants in the selected lines. Whole-genome re-sequencing was here used to identify individual loci where alleles that accumulated during directional selection were purged when selection was relaxed. In total, 11 loci with significant purging signals were identified, five in the low (LW) and six in the high (HW) body weight lineages. Associations between purged haplotypes in these loci and BW56 were tested in an advanced intercross line (AIL). Two loci with purging signals and haplotype associations to BW56 are particularly interesting for further functional characterization, one locus on chromosome 6 in the LW covering the sour-taste receptor gene PKD2L1, a functional candidate gene for the decreased appetite observed in the LWS and a locus on chromosome 20 in the HW containing a skeletal muscle hypertrophy gene, DNTTIP1.

Genetic parameters for clutch and broodiness traits in turkeys (Meleagris Gallopavo) and their relationship with body weight and egg production

The objective of this study was to estimate phenotypic and genetic parameters for clutch and broodiness (BR) traits in turkeys and their relationship with body weight and egg production. Data on dam line hens was available and included: body weight at 18 wk of age (BW18), body weight at lighting (BWL, 29 to 33 wk), age at first egg (AFE), egg number (EN), rate of lay (RL), clutch length (CL), maximum clutch length (MCL), pause length (PL), maximum PL (MPL) and BR. BR was defined as the average number of consecutive pause days between clutches that was higher than the average PL per hen. Heritability estimates for BW18 and BWL were 0.50 and 0.53, respectively. The heritability for egg production, clutch, and pause traits varied from low (MPL = 0.15; BR = 0.15) to moderate (AFE = 0.22; EN = 0.28; RL = 0.29; CL = 0.21; MCL = 0.27; PL = 0.25). Genetic correlations were negative between body weight traits and EN (r_{g (BW18, EN)} = −0.27; r_{g(BWL, EN)} = −0.33) and CL (r_{g(BW18, CL)} = −0.40; r_{g(BWL, CL)} = −0.33). BR was negatively genetically correlated with EN (r_{g(BR, EN)} = −0.85) and CL (r_{g(BR, CL)} = −0.30), and positively genetically correlated with PL (r_{g(BR, PL)} = 0.93) and AFE (r_{g(BR, AFE)} = 0.21). EN had a positive (0.73) and a negative (−0.84) genetic correlation with CL and PL, respectively. Overall, the results of this study confirmed the negative (unfavorable) correlations between egg production and body weight. Despite unfavorable genetic and phenotypic correlations between egg production traits and those relating to BR, the inclusion of BR in a selection program through incorporation of clutch length traits and pause length traits is feasible. Integration of either clutch length traits or pause length traits in a selection index is likely to increase egg number while decreasing broodiness.

Yield Performance, Laying Behaviour Traits and Egg Quality of Purebred and Hybrid Hens Reared under Outdoor Conditions

Simple Summary

Nowadays, the consumer has a choice of table eggs depending on rearing system, egg size, eggshell and yolk colour, produced by a few worldwide hybrid strains. Many countries have a historic poultry tradition, with breeds used for domestic egg production. In Northern Italy, in the Po river plain, a relevant part of the maize production was used also for poultry feeding and many breeds existed in the last centuries. At present, the Veneto region shows an important poultry biodiversity with many breeds, whose egg production traits are not well known. A comparison between purebred and hybrid hens throughout the laying period was carried out in order to study yield performance, laying behavioural traits and egg quality. The knowledge of the variation of these traits is useful for managing the breeding of these genotypes and the offer to the consumers. The local breeds showed differences from hybrid hens both for yield performance and laying behaviour and for both egg size and external and internal quality.

Abstract

This study compared the yield performance, laying behavioural traits and egg quality of purebred and hybrid hens (from 28 until 44 weeks of age, considering four periods) reared under outdoor conditions. The four genotypes were reared on the same trial station, on four areas (one genotype/area), and under the same environmental conditions from hatching until the end of the trial. Italian dual-purpose purebred (Ermellinata di Rovigo—ER and Robusta maculata—RM) and hybrid (Hy-Line Brown—HB and Hy-Line White 36—HW) hens (flock size: 70 birds/genotype) were allowed outdoors (4 m²/bird, good pasture during the growing period and poor pasture throughout the laying period, according to the season) and indoors (0.20 m²/bird, five birds/individual nest) and fed commercial feed. Significant (p < 0.01) differences among genotypes were found. The hybrids showed a higher laying rate and hen-day edible egg mass, and a lower body weight than the purebreds. Broken and out-of-nest egg% were higher in RM and HW than ER and HB, respectively. Double-yolk egg% was higher in hybrids than in purebreds. The eggshell colour varied among brown eggshell ER, RM, and HB. The ER showed the lowest shape index. With aging, the yolk to albumen ratio linearly increased in all groups, eggshell% changed in ER, HW, RM (cubic) and in HB (linear). The purebreds showed meat spots% higher than blood spots; HW showed the lowest total inclusion%. In conclusion, according to an egg scoring evaluation (egg weight = medium-large size, yolk to albumen ratio = 0.5, total inclusions = none), HW showed a higher quality than HB and RM, and ER was intermediate. The RM hens showed the highest% of defective eggs, especially for overcrowding at nest, HB showed the lowest. Under outdoor conditions the laying behaviour of the purebred hens and the nest management are important factors for the saleable egg rate.

Gene Expression Profiling in Ovaries and Association Analyses Reveal HEP21 as a Candidate Gene for Sexual Maturity in Chickens

Simple Summary

Chicken meat and egg productions are essential for human beings. Sexual maturity is important for both egg production and meat flavor. It is necessary to elucidate the genetic mechanism of chicken sexual maturity. In current study, we used digital gene expression (DGE) RNA-sequencing analysis to investigate differential expression of genes in pre-pubertal and post-pubertal ovaries in two different sub-breeds of chicken with different onsets of sexual maturity. After the analysis of RNA-sequencing data, numerous differentially expressed genes were found in both comparisons (32 day old, early-sexual-maturity pre-laying hens (P-F-O1) vs. 103 day old early-sexual-maturity laying hens (P-F-O2), and 32 day old late-sexual-maturity pre-laying hens (L-F-O1) vs. 153 day old late-sexual-maturity pre-laying hens (L-F-O2)). With the bioinformatic analysis, hen egg protein 21 kDa (HEP21) was chosen as the candidate gene to conduct following experiment. The variations in HEP21 were screened and association analyses between rs315156783 and reproductive traits were investigated in fifth-generation Ningdu Yellow chickens from a closely bred population. These results demonstrated that HEP21 is a candidate gene for sexual maturity and ovary development in chickens. However, the underlying mechanism of how HEP21 regulates chicken sexual maturity needs further focused studies.

Abstract

The age of onset of sexual maturity is an important reproductive trait in chickens. In this study, we explored candidate genes associated with sexual maturity and ovary development in chickens. We performed DGE RNA-sequencing analyses of ovaries of pre-laying (P-F-O1, L-F-O1) and laying (P-F-O2, L-F-O2) hens of two sub-breeds of Ningdu Yellow chicken. A total of 3197 genes were identified in the two comparisons, and 966 and 1860 genes were detected exclusively in comparisons of P-F-O1 vs. P-F-O2 and L-F-O1 vs. L-F-O2, respectively. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that genes involved in transmembrane signaling receptor activity, cell adhesion, developmental processes, the neuroactive ligand–receptor interaction pathway, and the calcium signaling pathway were enriched in both comparisons. Genes on these pathways, including growth hormone (GH), integrin subunit beta 3 (ITGB3), thyroid stimulating hormone subunit beta (TSHB), prolactin (PRL), and transforming growth factor beta 3 (TGFB3), play indispensable roles in sexual maturity. As a gene unique to poultry, hen egg protein 21 kDa (HEP21) was chosen as the candidate gene. Differential expression and association analyses were performed. RNA-seq data and qPCR showed that HEP21 was significantly differentially expressed in pre-pubertal and pubertal ovaries. A total of 23 variations were detected in HEP21. Association analyses of single nucleotide polymorphisms (SNPs) in HEP21 and reproductive traits showed that rs315156783 was significantly related to comb height at 84 and 91 days. These results indicate that HEP21 is a candidate gene for sexual maturity in chickens. Our results contribute to a more comprehensive understanding of sexual maturity and reproduction in chickens.

Glucose Tolerance and Plasma Non-Esterified Fatty Acid Levels in Chickens Selected for Low Body Weight, Red Junglefowl, and their Reciprocal Cross

Responses of an individual to food deprivation, such as a 16-h fast, are complex, and are influenced by environmental and genetic factors. Domestication is an ongoing process during which adaptations to changing environments occur over generations. Food deprivation by their caretakers is less for domestic chickens than for their junglefowl ancestors. Unlike domestic chicken, the junglefowl adapted over generations to periods of food deprivation, which may be reflected in differences in metabolic responses to brief periods without food. Here, we compared the blood glucose and plasma levels of non-esterified fatty acids (NEFA) among four populations when deprived of feed for 16 h. The four populations included a domestic White Rock experimental line (LWS) maintained for generations under ad libitum feeding, adult red junglefowl (RJF), and a reciprocal cross of the lines. Although there were significant differences in adult (31-week) body weight between the RJF (683 g) and LWS (1282 g), with the weight of F₁ crosses being intermediate, the amount of abdominal fat relative to body weight was similar for all populations. Patterns for blood glucose responses to a glucose bolus after a 16-h fast were similar for the initial and final points in the parental and cross populations. However, RJF reached their peak faster than LWS, with the reciprocal cross intermediate to the parental populations. Plasma NEFA concentrations were higher after the 16-h fast than in fed states, with no population differences for the fasting state. However, in the fed state, NEFA levels were lesser for LWS than for others, which was reflected further in percentage change from fed to fasted. This larger change in LWS suggests differences in mobilization of energy substrates and implies that during domestication or development of the LWS line, thresholds for responses to acute stressors may have increased.

Bidirectional Selection for Body Weight on Standing Genetic Variation in a Chicken Model

Experimental populations of model organisms provide valuable opportunities to unravel the genomic impact of selection in a controlled system. The Virginia body weight chicken lines represent a unique resource to investigate signatures of selection in a system where long-term, single-trait, bidirectional selection has been carried out for more than 60 generations. At 55 generations of divergent selection, earlier analyses of pooled genome resequencing data from these lines revealed that 14.2% of the genome showed extreme differentiation between the selected lines, contained within 395 genomic regions. Here, we report more detailed analyses of these data exploring the regions displaying within- and between-line genomic signatures of the bidirectional selection applied in these lines. Despite the strict selection regime for opposite extremes in body weight, this did not result in opposite genomic signatures between the lines. The lines often displayed a duality of the sweep signatures, where an extended region of homozygosity in one line, in contrast to mosaic pattern of heterozygosity in the other line. These haplotype mosaics consisted of short, distinct haploblocks of variable between-line divergence, likely the results of a complex demographic history involving bottlenecks, introgressions and moderate inbreeding. We demonstrate this using the example of complex haplotype mosaicism in the growth1 QTL. These mosaics represent the standing genetic variation available at the onset of selection in the founder population. Selection on standing genetic variation can thus result in different signatures depending on the intensity and direction of selection.

Genomic signatures of 60 years of bidirectional selection for 8-week body weight in chickens

Sixty years, constituting 60 generations, have passed since the founding of the Virginia body weight lines, an experimental population of White Plymouth Rock chickens. Using a stringent breeding scheme for divergent 8-week body weight, the lines, which originated from a common founder population, have responded to bidirectional selection with an approximate 15-fold difference in the selected trait. They provide a model system to study the genetics of complex traits in general and the influences of artificial selection on quantitative genetic architectures in particular. As we reflect on the 60th anniversary of the initiation of the Virginia body weight lines, there is opportunity to discuss the findings obtained using different analytical and experimental genetic and genomic strategies and integrate them with a recent pooled genome resequencing dataset. Hundreds of regions across the genome show differentiation between the 2 lines, reinforcing previous findings that response to selection relied on standing variation across many genes and giving insights into the haplotype complexity underlying regions associated with body weight.

Growth patterns for three generations of an intercross between red junglefowl and chickens selected for low body weight

Growth is a complex and dynamic process that may be measured at a specific point or over a period of time. Compared was the growth of male and female chickens over a three-generation period. Involved were red junglefowl (RJF; Gallus gallus), a line of White Plymouth Rock chickens (LWS; Gallus gallus domesticus) selected for low body weight, and their reciprocal F1 and F2 crosses. In both sexes, Gompertz's description of growth showed that RJF had significantly lower asymptotes, earlier inflection points, and faster growth rates than LWS. Heterosis for these measures was positive for asymptote and negative for growth rate and inflection point. The RJF commenced egg production at a significantly younger age and lower body weight than LWS. Although F1 and F2 reciprocal crosses were similar for body weight and for age at first egg, the F1 reciprocal crosses began lay at significantly younger ages than the F2 crosses and parental lines. When viewed on a physiological basis where age and body weight were simultaneously standardized, both parental lines and reciprocal F1 and F2 crosses had differing rapid and lag growth phases. Overall, sexual dimorphism increased in all populations from hatch to sexual maturity. The LWS males had a longer growth period consistent with their female counterparts who became sexually mature at older ages. Comprehensively, these results indicate additive and nonadditive genetic variation for distinct growth patterns and changes in resource allocation strategies over time.