Genome-wide association analysis of egg production performance in chickens across the whole laying period

Integrative analysis of the ovarian metabolome and transcriptome of the Yaoshan chicken and its improved hybrids

Introduction: Laying performance is a key factor affecting production efficiency in poultry, but its molecular mechanism is still indistinct. In this study, Yaoshan chickens, a local breed in Guizhou, China, and merchant chickens (GYR) with higher egg yield after the three-line cross improvement hybridization of Yaoshan chickens were used as animal samples. Methods: To explore the regulatory mechanism of the diversities in laying performance, RNA-seq and ultra-performance liquid chromatographytandem mass spectrometry (UPLC-MS/MS) were used to describe the transcriptional and metabolic profiles of the ovaries of Yaoshan and GYR chickens. Results: At the transcriptional level, 288 differentially expressed genes were upregulated in Yaoshan chickens and 353 differentially expressed genes were upregulated in GYR chickens. In addition, GSEA showed that ECM-receptor interactions and the TGF-β signaling pathway were restrained, resulting in increased egg production in GYR chickens. Furthermore, the upregulation of thiamine and carnitine was identified by metabolomic analysis to facilitate the laying performance of hens. Finally, comprehensive analyses of the transcriptome and metabolome found that thiamine and carnitine were negatively correlated with ECM-receptor interactions and the TGF-β signaling pathway, which jointly regulate the laying performance of Yaoshan chickens and GYR chickens. Discussion: Taken together, our research delineates differences in the transcriptional and metabolic profiles of the ovaries of Yaoshan and GYR chickens during the peak egg production period and provides new hypotheses and clues for further research on poultry egg production performance and the improvement of economic benefits.

Investigating the genetic determination of duration-of-fertility trait in breeding hens

The duration-of-fertility (DF), which was defined as the number of days when breeding hens lay fertile eggs following copulation or artificial insemination (AI), is an important economic trait in chick production when it has strong effects on fertile egg output and production costs. Little is known about the underlying genes and molecular markers related to DF trait to date. Here, we measured the DF of 701 Chinese Jinghong hens and 408 Jingfen hens. The DF showed high individual variability and potential for genetic improvement. Then, 192 Jinghong breeding hens were provided for a genome-wide association study, 27 SNPs respectively located in three genomic linkage regions (GGA1:41Kb; GGA3:39Kb and GGA8:39Kb) were suggested to be significantly associated with DF. Particularly, 6 of these 27 SNPs were further verified to be associated with DF in the 701 Jinghong and 408 Jingfen hens using PCR-RFLP genotyping method. These 27 SNPs were also mapped to 7 genes according to their genomic position. Furtherly, 5 of these 7 genes were tested using qPCR. Results show that the CYP2D6, WBP2NL, ESR1 and TGFBR3 mRNA expression levels of hens with long DF were significantly higher than the hens with short DF (P < 0.05). Overall, findings in our research provide new insight into the genetic basis of duration-of-fertility in breeding hens while providing new clues for further functional validation on the DF-related genetic regulation mechanism and improvement of DF through chicken breeding.

Genome-wide association study identified candidate genes for egg production traits in the Longyan Shan-ma duck

Egg production is an important economic trait in layer ducks and understanding the genetics basis is important for their breeding. In this study, a genome-wide association study (GWAS) for egg production traits in 303 female Longyan Shan-ma ducks was performed based on a genotyping-by-sequencing strategy. Sixty-two single nucleotide polymorphisms (SNPs) associated with egg weight traits were identified (P < 9.48 × 10-5), including 8 SNPs at 5% linkage disequilibrium (LD)-based Bonferroni-corrected genome-wide significance level (P < 4.74 × 10-6). One hundred and nineteen SNPs were associated with egg number traits (P < 9.48 × 10-5), including 13 SNPs with 5% LD-based Bonferroni-corrected genome-wide significance (P < 4.74 × 10-6). These SNPs annotated 146 target genes which contained known candidate genes for egg production traits, such as prolactin and prolactin releasing hormone receptor. This study identified that these associated genes were significantly enriched in egg production-related pathways (P < 0.05), such as the oxytocin signaling, MAPK signaling, and calcium signaling pathways. It was notable that 18 genes were differentially expressed in ovarian tissues between higher and lower egg production in Shan-ma ducks. The identified potential candidate genes and pathways provide insight into the genetic basis underlying the egg production trait of layer ducks.

Genome-Wide Association Analysis and Genetic Parameters for Egg Production Traits in Peking Ducks

Egg production traits are crucial in the poultry industry, including age at first egg (AFE), egg number (EN) at different stages, and laying rate (LR). Ducks exhibit higher egg production capacity than other poultry species, but the genetic mechanisms are still poorly understood. In this study, we collected egg-laying data of 618 Peking ducks from 22 to 66 weeks of age and genotyped them by whole-genome resequencing. Genetic parameters were calculated based on SNPs, and a genome-wide association study (GWAS) was performed for these traits. The SNP-based heritability of egg production traits ranged from 0.09 to 0.54. The GWAS identified nine significant SNP loci associated with AFE and egg number from 22 to 66 weeks. These loci showed that the corresponding alleles were positively correlated with a decrease in the traits. Moreover, three potential candidate genes (ENSAPLG00020011445, ENSAPLG00020012564, TMEM260) were identified. Functional enrichment analyses suggest that specific immune responses may have a critical impact on egg production capacity by influencing ovarian function and oocyte maturation processes. In conclusion, this study deepens the understanding of egg-laying genetics in Peking duck and provides a sound theoretical basis for future genetic improvement and genomic selection strategies in poultry.

Whole genome resequencing reveals the adaptability of native chickens to drought, tropical and frigid environments in Xinjiang

Chickens exhibit extensive genetic diversity and are distributed worldwide. Different chicken breeds have evolved to thrive in diverse environmental conditions. However, research on the genetic mechanisms underlying chicken adaptation to extreme environments, such as tropical, frigid and drought-prone regions, remains limited. In this study, we conducted whole-genome sequencing of 240 individuals from six native chicken breeds in Xinjiang, China, as well as 4 publicly available chicken breeds inhabiting regions with varying annual precipitations, temperatures, and altitudes. Our analysis revealed several genetic variants among the examined breeds. Furthermore, we investigated the genetic diversity and population structure of breeds residing in extreme drought and temperature environments by comparing them. Notably, native chicken breeds exhibited different genetic diversity and population structures. Moreover, we identified candidate genes associated with chicken adaptability to the environment, such as CORO2A, CTNNA3, AGMO, GRID2, BBOX1, COL3A1, INSR, SOX5, MAP2 and PLPPR1. Additionally, pathways such as lysosome, cysteine and methionine metabolism, glycosaminoglycan degradation, and Wnt signaling may be play crucial roles in regulating chicken adaptation to drought environments. Overall, these findings contribute to our understanding of the genetic mechanisms governing chicken adaptation to extreme environments, and also offer insights for enhancing the resilience of chicken breeds to different climatic conditions.

Genetic Analysis of Egg Production Traits in Luhua Chickens: Insights from a Multi-Trait Animal Model and a Genome-Wide Association Study

Egg production plays a pivotal role in the economic viability of hens. To analyze the genetic rules of egg production, a total of 3151 Luhua chickens were selected, the egg production traits including egg weight at first laying (Start-EW), egg weight at 43 weeks (EW-43), egg number at 43 weeks (EN-43), and total egg number (EN-All) were recorded. Then, the effects of related factors on egg production traits were explored, using a multi-trait animal model for genetic parameter estimation and a genome-wide association study (GWAS). The results showed that body weight at first egg (BWFE), body weight at 43 weeks (BW-43), age at first egg (AFE), and seasons had significant effects on the egg production traits. Start-EW and EW-43 had moderate heritability of 0.30 and 0.21, while EN-43 and EN-All had low heritability of 0.13 and 0.16, respectively. Start-EW exhibited a robust positive correlation with EW-43, while Start-EW was negatively correlated with EN-43 and EN-All. Furthermore, gene ontology (GO) results indicated that Annexin A2 (ANXA2) and Frizzled family receptor 7 (FZD7) related to EW-43, Cyclin D1 (CCND1) and A2B adenosine receptor (ADORA2B) related to EN-All, and have been found to be mainly involved in metabolism and growth processes, and deserve more attention and further study. This study contributes to accelerating genetic progress in improving low heritability egg production traits in layers, especially in Luhua chickens.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Genome-wide association studies of egg production traits by whole genome sequencing of Laiwu Black chicken

Compared to high-yield commercial laying hens, Chinese indigenous chicken breeds have poor egg laying capacity due to the lack of intensive selection. However, as these breeds have not undergone systematic selection, it is possible that there is a greater abundance of genetic variations related to egg laying traits. In this study, we assessed 5 egg number (EN) traits at different stages of the egg-laying period: EN1 (from the first egg to 23 wk), EN2 (from 23 to 35 wk), EN3 (from 35 to 48 wk), EN4 (from the first egg to 35 wk), and EN5 (from the first egg to 48 wk). To investigate the molecular mechanisms underlying egg number traits in a Chinese local chicken breed, we conducted a genome-wide association study (GWAS) using data from whole-genome sequencing (WGS) of 399 Laiwu Black chickens. We obtained a total of 3.01 Tb of raw data with an average depth of 7.07 × per individual. A total of 86 genome-wide suggestive or significant single-nucleotide polymorphisms (SNP) contained within a set of 45 corresponding candidate genes were identified and found to be associated with stages EN1-EN5. The genes vitellogenin 2 (VTG2), lipase maturation factor 1 (LMF1), calcium voltage-gated channel auxiliary subunit alpha2delta 3 (CACNA2D3), poly(A) binding protein cytoplasmic 1 (PABPC1), programmed cell death 11 (PDCD11) and family with sequence similarity 213 member A (FAM213A) can be considered as the candidate genes associated with egg number traits, due to their reported association with animal reproduction traits. Noteworthy, results suggests that VTG2 and PDCD11 are not only involved in the regulation of EN3, but also in the regulation of EN5, implies that VTG2 and PDCD11 have a significant influence on egg production traits. Our study offers valuable genomic insights into the molecular genetic mechanisms that govern egg number traits in a Chinese indigenous egg-laying chicken breed. These findings have the potential to enhance the egg-laying performance of chickens.

Genetic parameter estimation and molecular foundation of chicken egg-laying trait

The age of first egg (AFE) in chicken can affect early and even life-time egg production performance to some extent, and therefore is an important economic trait that affects production efficiency. To better understand the genetic patterns of AFE and other production traits including body weight at first egg (BWA), first egg weight (FEW), and total egg number from AFE to 58 wk of age (total-EN), we recorded the production performance of 2 widely used layer breeds, white leghorn (WL) and Rhode Island Red (RIR) and estimated genetic parameters based on pedigree and production data. The results showed that the heritability of AFE in both breeds ranged from 0.4 to 0.6, and AFE showed strong positive genetic and phenotypic correlations to BWA as well as FEW, while showing strong negative genetic and phenotypic correlations with total-EN. Furtherly, by genome-wide association analysis study (GWAS), we identified 12 and 26 significant SNPs to be related to AFE in the 2-layer breeds, respectively. A total of 18 genes were identified that could affect AFE based on the significant SNP annotations obtained, but there were no gene overlapped in the 2 breeds indicating the genetic foundation of AFE could differ from breed to breed. Our results provided a deeper understanding of genetic patterns and molecular basement of AFE in different breeds and could help in the selection of egg production traits.

New insights into the genetic loci related to egg weight and age at first egg traits in broiler breeder

Egg weight (EW) and age at first egg (AFE) are economically important traits in breeder chicken production. The genetic basis of these traits, however, is far from understood, especially for broiler breeders. In this study, genetic parameter estimation, genome-wide association analysis, meta-analysis, and selective sweep analysis were carried out to identify genetic loci associated with EW and AFE in 6,842 broiler breeders. The study found that the heritability of EW ranged from 0.42 to 0.44, while the heritability of AFE was estimated at 0.33 in the maternal line. Meta-analysis and selective sweep analysis identified two colocalized regions on GGA4 that significantly influenced EW at 32 wk (EW32W) and at 43 wk (EW43W) with both paternal and maternal lines. The genes AR, YIPF6, and STARD8 were located within the significant region (GGA4: 366.86-575.50 kb), potentially affecting EW through the regulation of follicle development, cell proliferation, and lipid transfer etc. The promising genes LCORL and NCAPG were positioned within the significant region (GGA4:75.35-75.42 Mb), potentially influencing EW through pleiotropic effects on growth and development. Additionally, 3 significant regions were associated with AFE on chromosomes GGA7, GGA19, and GGA27. All of these factors affected the AFE by influencing ovarian development. In our study, the genomic information from both paternal and maternal lines was used to identify genetic regions associated with EW and AFE. Two genomic regions and eight genes were identified as the most likely candidates affecting EW and AFE. These findings contribute to a better understanding of the genetic basis of egg production traits in broiler breeders and provide new insights into future technology development.

TAC3 regulates GnRH/gonadotropin synthesis in female chickens

Neurokinin B (NKB), a peptide encoded by the tachykinin 3 (TAC3), is critical for reproduction in all studied species. However, its potential roles in birds are less clear. Using the female chicken (c-) as a model, we showed that cTAC3 is composed of five exons with a full-length cDNA of 787 bp, which was predicted to generate the mature NKB peptide containing 10 amino acids. Using cell-based luciferase reporter assays, we demonstrated that cNKB could effectively and specifically activate tachykinin receptor 3 (TACR3) in HEK293 cells, suggesting its physiological function is likely achieved via activating cTACR3 signaling. Notably, cTAC3 and cTACR3 were predominantly and abundantly expressed in the hypothalamus of hens and meanwhile the mRNA expression of cTAC3 was continuously increased during development, suggesting that NKB-TACR3 may emerge as important components of the neuroendocrine reproductive axis. In support, intraperitoneal injection of cNKB could significantly promote hypothalamic cGnRH-Ι, and pituitary cFSHβ and cLHβ expression in female chickens. Surprisingly, cTAC3 and cTACR3 were also expressed in the pituitary gland, and cNKB treatment significantly increased cLHβ and cFSHβ expression in cultured primary pituitary cells, suggesting cNKB can also act directly at the pituitary level to stimulate gonadotropin synthesis. Collectively, our results reveal that cNKB functionally regulate GnRH/gonadotropin synthesis in female chickens.

Genetic parameters, reciprocal cross differences, and age-related heterosis of egg-laying performance in chickens

BACKGROUND

Egg-laying performance is economically important in poultry breeding programs. Crossbreeding between indigenous and elite commercial lines to exploit heterosis has been an upward trend in traditional layer breeding for niche markets. The objective of this study was to analyse the genetic background and to estimate the heterosis of longitudinal egg-laying traits in reciprocal crosses between an indigenous Beijing-You and an elite commercial White Leghorn layer line. Egg weights were measured for the first three eggs, monthly from 28 to 76 weeks of age, and at 86 and 100 weeks of age. Egg quality traits were measured at 32, 54, 72, 86, and 100 weeks of age. Egg production traits were measured from the start of lay until 43, 72, and 100 weeks of age. Heritabilities and phenotypic and genetic correlations were estimated. Heterosis was estimated as the percentage difference of performance of a crossbred from that of the parental average. Reciprocal cross differences were estimated as the difference between the reciprocal crossbreds as a percentage of the parental average.

RESULTS

Estimates of heritability of egg weights ranged from 0.29 to 0.75. Estimates of genetic correlations between egg weights at different ages ranged from 0.72 to 1.00. Estimates of heritability for cumulative egg numbers until 43, 72, and 100 weeks of age were around 0.15. Estimates of heterosis for egg weight and cumulative egg number increased with age, ranging from 1.0 to 9.0% and from 1.4 to 11.6%, respectively. From 72 to 100 weeks of age, crossbreds produced more eggs per week than the superior parent White Leghorn (3.5 eggs for White Leghorn, 3.8 and 3.9 eggs for crossbreds). Heterosis for eggshell thickness ranged from 2.7 to 6.6% when using Beijing-You as the sire breed. No significant difference between reciprocal crosses was observed for the investigated traits, except for eggshell strength at 54 weeks of age.

CONCLUSIONS

The heterosis was substantial for egg weight and cumulative egg number, and increased with age, suggesting that non-additive genetic effects are important in crossbreds between the indigenous and elite breeds. Generally, the crossbreds performed similar to or even outperformed the commercial White Leghorns for egg production persistency.

Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.

Genome-Wide Association Study of Egg Production Traits in Shuanglian Chickens Using Whole Genome Sequencing

Egg production is the most important economic trait in laying hens. To identify molecular markers and candidate genes associated with egg production traits, such as age at first egg (AFE), weight at first egg (WFE), egg weight (EW), egg number (EN), and maximum consecutive egg laying days (MCD), a genome-wide analysis by whole genome sequencing was performed in Shuanglian chickens. Through whole genome sequencing and quality control, a total of 11,006,178 SNPs were obtained for further analysis. Heritability estimates ranged from moderate to high for EW (0.897) and MCD (0.632), and from low to moderate (0.193~0.379) for AFE, WFE, and EN. The GWAS results showed 11 genome-wide significant SNPs and 23 suggestive significant SNPs were identified to be associated with EN, MCD, WFE, and EW. Linkage disequilibrium analysis revealed twenty-seven SNPs associated with EN were located in a 0.57 Mb region on GGA10, and clustered into five blocks. Through functional annotation, three candidate genes NEO1, ADPGK, and CYP11A1, were identified to be associated with EN, while the S1PR4, LDB2, and GRM8 genes was linked to MCD, WFE, and EW, respectively. These findings may help us to better understand the molecular mechanisms underlying egg production traits in chickens and contribute to genetic improvement of these traits.

Research Note: Genetic parameters for egg production and clutch-related traits in indigenous Beijing-You chickens

Egg products from indigenous chickens have growing market shares as consumers are pursuing differentiation in egg consumption. The genetic improvement in egg production performance of those breeds is crucial for increasing the economic profit. This study aimed to estimate genetic parameters for egg production and clutch-related traits in indigenous Beijing-You chickens for understanding the genetic architecture and exploring proper biological traits for selection. Data on traits including age at first egg (AFE), egg number (EN), average clutch length (ACL), maximum clutch length (MCL), number of clutches (NC) and pauses (NP), and average pause length (APL) were collected from 4 generations of purebred Beijing-You chickens based on the 43-wk and 66-wk of individual egg production record. The heritabilities, genetic and phenotypic correlations were analyzed by the DMU software with the restricted maximum likelihood method in a multivariate animal model. The results showed that the AFE of Beijing-You chickens was 174.45 d of age, and its heritability was as high as 0.62. The heritability was 0.26 for EN43 and 0.18 for EN66. The clutch traits including ACL, MCL, NC, and NP were moderate to high heritable (h2 = 0.15-0.39), but APL was very low heritable (h2 = 0.05). Genetic correlations were high between AFE and EN (rG(AFE, EN43) = -0.79, rG(AFE, EN66) = -0.39), whereas low between AFE and ACL (rG(AFE, ACL43) = -0.08, rG(AFE, ACL66) = 0.01) and MCL (rG(AFE, MCL) = -0.07). EN had higher correlations with ACL (rG(EN43, ACL43) = 0.59, rG(EN66, ACL66) = 0.40) than that with MCL (rG(EN43, MCL43) = 0.56, rG(EN66, MCL66) = 0.32). The heritability for ACL43 (h2 = 0.38) was higher than that for MCL43 (h2 = 0.33). ACL43 had a positive correlation with EN66 (rG(ACL43, EN66) = 0.62). These results indicated that the egg production of whole laying period could be improved by early selection for AFE and ACL at the same time in Beijing-You chickens.

Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds

BACKGROUND

The genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes.

RESULTS

Using the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) F_ST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF).

CONCLUSION

Based on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

Whole-genome sequencing identifies potential candidate genes for egg production traits in laying ducks (Anas platyrhynchos)

Egg production traits are economically important in laying ducks. Genetic molecular mechanisms and candidate genes underlying these traits remain unclear. In this study, whole genome variants were identified through whole-genome resequencing using three high-egg producing (HEN) and three low-egg producing (LEN) laying ducks. The gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genome (KEGG) pathways for the genes of common differential variants between HEN and LEN ducks were determined. Frizzled class receptor 6 (FZD6) was further genotyped using the Sequenom MassARRAY iPLEX platform. The association of FZD6 gene polymorphisms with 73 egg production and weight traits in 329 female ducks were estimated. A total of 65,535 single nucleotide polymorphisms (SNPs) and 4,702 indels were identified across the genome. Fourteen GO terms and 14 KEGG pathways were determined for the genes of common differential variants, including MAPK signaling, Wnt signaling, melanogenesis and calcium signaling pathways, which are key functional pathways for poultry egg production reported in previous reports. Further analysis showed that 27 SNPs of FZD6 were associated with three early egg production of duck and egg weight traits, including egg production at 17 weeks (EP17), 18 weeks (EP18) and 19 weeks (EP19) and egg weight at 59 weeks (EW59). The FZD6 should be considered a novel candidate gene for egg production traits in laying ducks.

Integration of transcriptome sequencing and whole genome resequencing reveal candidate genes in egg production of upright and pendulous-comb chickens

Egg production performance plays an important role in the poultry industry across the world. Previous studies have shown a great difference in egg production performance between pendulous-comb (PC) and upright-comb (UC) chickens. However, there are no reports to identify potential candidate genes for egg production in PC and UC chickens. In the present study, 1,606 laying chickens were raised, and the egg laid by individual chicken was collected for 100 d. Moreover, the expression level of estrogen and progesterone hormones was measured at the start-laying and peak-laying periods of hens. Besides, 4 PC and 4 UC chickens were selected at 217 d of age to perform transcriptome sequencing (RNA-seq) and whole genome resequencing (WGS) to screen the potential candidate genes of egg production. The results showed that PC chicken demonstrated better egg production performance (P < 0.05) and higher estrogen and progesterone hormone expression levels than UC chicken (P < 0.05). RNA-seq analysis showed that 341 upregulated and 1,036 downregulated differentially expressed genes (DEGs) were identified in the ovary tissues of PC and UC chickens. These DEGs were mainly enriched in protein-related, lipid-related, and nucleic acids-related biological processes including ribosome, peptide biosynthetic process, lipid transport terms, and catalytic activity acting on RNA which can significantly affect egg production in chicken. The enrichment results of WGS analysis were consistent with RNA-seq. Further, joint analysis of WGS and RNA-seq data was utilized to screen 30 genes and CAMK1D, CLSTN2, MAST2, PIK3C2G, TBC1D1, STK3, ADGRB3, and PPARGC1A were identified as potential candidate genes for egg production in PC and UC chickens. In summary, our study provides a wealth of information for a better understanding of the genetic and molecular mechanism for the future breeding of PC and UC chickens for egg production.

Investigation of egg production curve in ostrich using nonlinear functions

In most countries, ostrich farming is considered a developing branch of the efficient poultry industry. The profitability of ostrich farm requires specific consideration of productions features such as the female fertility, egg production, hatchability, and growth performance. Hence, this study aimed to fit nonlinear functions to describe the ostrich egg production pattern to achieve the most appropriate and recommendable mathematical function for future studies. For this purpose, 14,507 daily records of 184 female ostriches in 5 production seasons (periods) during 2016 to 2021 were used. Five nonlinear functions including Incomplete gamma (Wood function), Corrected gamma (McNally), nonlinear Logistic (Yang), Logistic (Nelder), and Lokhorst were fitted for modeling the egg production curve in ostrich. The goodness of fit criteria's including Mean Square Error (MSE), Likelihood Ratio Test (LRT), Akaike Information Criterion (AIC), and Bayesian Information Criterion (BIC) were used to evaluate and selection of the best function. The results indicated that the Wood and the McNally functions with a slight difference in all fitting criteria were the best-fitted functions and the Yang function with the highest values of MSE, LRT, AIC, BIC, were the most inappropriate function to describe the ostrich egg production curve. The McNally and the Wood can be recommended as appropriate functions to describe egg production during 5 production seasons in the studied ostrich flock.

Advanced Practical Strategies to Enhance Table Egg Production

The global demand for table eggs has increased exponentially due to the growing human population. To meet this demand, major advances in hen genetics, nutrition, and husbandry procedures are required. Developing cost-effective and practically applicable strategies to improve egg production and quality is necessary for the development of egg industry worldwide. Consumers have shown a strong desire regarding the improvement of hens' welfare and egg quality. They also become interested in functional and designer foods. Modifications in the nutritional composition of laying hen diets significantly impact egg nutritional composition and quality preservation. According to previous scientific research, enriched egg products can benefit human health. However, producers are facing a serious challenge in optimizing breeding, housing, and dietary strategies to ensure hen health and high product quality. This review discussed several practical strategies to increase egg production, quality, and hens' welfare. These practical strategies can potentially be used in layer farms for sustainable egg production. One of these strategies is the transition from conventional to enriched or cage-free production systems, thereby improving bird behavior and welfare. In addition, widely use of plant/herbal substances as dietary supplements in layers' diets positively impacts hens' physiological, productive, reproductive, and immunological performances.

A significant quantitative trait locus on chromosome Z and its impact on egg production traits in seven maternal lines of meat-type chicken

BACKGROUND

Egg production is economically important in the meat-type chicken industry. To better understand the molecular genetic mechanism of egg production in meat-type chicken, genetic parameter estimation, genome-wide association analyses combined with meta-analyses, Bayesian analyses, and selective sweep analyses were performed to screen single nucleotide polymorphisms (SNPs) and other genetic loci that were significantly associated with egg number traits in 11,279 chickens from seven material lines.

RESULTS

Yellow-feathered meat-type chickens laid 115 eggs at 43 weeks of age and white-feathered chickens laid 143 eggs at 60 weeks of age, with heritability ranging from 0.034-0.258. Based on meta-analyses and selective sweep analyses, one region (10.81-13.05 Mb) on chromosome Z was associated with egg number in all lines. Further analyses using the W2 line was also associated with the same region, and 29 SNPs were identified that significantly affected estimation of breeding value of egg numbers. The 29 SNPs were identified as having a significant effect on the egg number EBV in 3194 birds in line W2. There are 36 genes in the region, with glial cell derived neurotrophic factor, DAB adaptor protein 2, protein kinase AMP-activated catalytic subunit alpha 1, NAD kinase 2, mitochondrial, WD repeat domain 70, leukemia inhibitory factor receptor alpha, complement C6, and complement C7 identified as being potentially affecting to egg number. In addition, three SNPs (rs318154184, rs13769886, and rs313325646) associated with egg number were located on or near the prolactin receptor gene.

CONCLUSION

Our study used genomic information from different chicken lines and populations to identify a genomic region (spanning 2.24 Mb) associated with egg number. Nine genes and 29 SNPs were identified as the most likely candidate genes and variations for egg production. These results contribute to the identification of candidate genes and variants for egg traits in poultry.

Genome-Wide Association Study and Selective Sweep Analysis Reveal the Genetic Architecture of Body Weights in a Chicken F2 Resource Population

Rapid growth is one of the most important economic traits in broiler breeding programs. Identifying markers and genes for growth traits may not only benefit marker-assisted selection (MAS)/genomic selection (GS) but also provide important information for understanding the genetic architecture of growth traits in broilers. In the present study, an F₂ resource population derived from a cross between the broiler and Baier yellow chicken (a Chinese local breed) was used and body weights from 1 to 12 weeks of age [body weight (BW) 1–BW12)] were measured. A total of 519 F₂ birds were genome re-sequenced, and a combination of genome-wide association study (GWAS) and selective sweep analysis was carried out to characterize the genetic architecture affecting chicken body weight comprehensively. As a result, 1,539 SNPs with significant effects on body weights at different weeks of age were identified using a genome-wide efficient mixed-model association (GEMMA) package. These SNPs were distributed on chromosomes 1 and 4. Besides, windows under selection identified for BW1–BW12 varied from 1,581 to 2,265. A total of 42 genes were also identified with significant effects on BW1–BW12 based on both GWAS and selective sweep analysis. Among these genes, diacylglycerol kinase eta (DGKH), deleted in lymphocytic leukemia (DLEU7), forkhead box O17 (FOXO1), karyopherin subunit alpha 3 (KPNA3), calcium binding protein 39 like (CAB39L), potassium voltage-gated channel interacting protein 4 (KCNIP4), and slit guidance ligand 2 (SLIT2) were considered as important genes for broiler growth based on their basic functions. The results of this study may supply important information for understanding the genetic architecture of growth traits in broilers.

Genome-Wide Association Study of Egg-Laying Traits and Egg Quality in LingKun Chickens

Egg production is the most important trait of laying hens. To identify molecular markers and candidate genes associated with egg production and quality, such as body weight at first oviposition (BWF), the number of eggs produced in 500 days (EN500), egg weight (EW), egg shell thickness (EST), egg shell strength (ESS), and Haugh unit (HU), a genome-wide analysis was performed in 266 LingKun Chickens. The results showed that thirty-seven single nucleotide polymorphisms (SNPs) were associated with all traits (p < 9.47 × 10⁻⁸, Bonferroni correction). These SNPs were located in close proximity to or within the sequence of the thirteen candidate genes, such as Galanin And GMAP Prepropeptide (GAL), Centromere Protein (CENPF), Glypican 2 (GPC2), Phosphatidylethanolamine N-Methyltransferase (PEMT), Transcription Factor AP-2 Delta (TFAP2D), and Carboxypeptidase Q (CPQ) gene related to egg-laying and Solute Carrier Family 5 Member 7 (SLC5A7), Neurocalcin Delta (NCALD), Proteasome 20S Subunit Beta 2 (PSMB2), Slit Guidance Ligand 3 (SLIT3), and Tubulin Tyrosine Ligase Like 7 (TTLL7) genes related to egg quality. Interestingly, one of the genes involved in bone formation (SLIT3) was identified as a candidate gene for ESS. Our candidate genes and SNPs associated with egg-laying traits were significant for molecular breeding of egg-laying traits and egg quality in LingKun chickens.

Host-genotype-dependent cecal microbes are linked to breast muscle metabolites in Chinese chickens

In chickens, the effect of host genetics on the gut microbiota is not fully understood, and the extent to which the heritable gut microbes affect chicken metabolism and physiology is still an open question. Here, we explored the interactions among chicken genetics, the cecal microbiota and metabolites in breast muscle from ten chicken breeds in China. We found that different chicken breeds displayed distinct cecal microbial community structures and functions, and 15 amplicon sequence variants (ASVs) were significantly associated with host genetics through different genetic loci, such as those related to the intestinal barrier function. We identified five heritable ASVs significantly associated with 53 chicken muscle metabolites, among which the Megamonas probably affected lipid metabolism through the production of propionate. Our study revealed that the chicken genetically associated cecal microbes may have the potential to affect the bird’s physiology and metabolism.

•

The cecal microbiota are different among ten chicken breeds

•

The chicken genetics influences the cecal microbiota structures and functions

•

The chicken heritable cecal microbes are associated with muscle metabolites

•

Megamonas may affect lipid metabolism by the production of propionate

Using egg production longitudinal recording to study the genetic background of resilience in purebred and crossbred laying hens

Background

There is growing interest in using genetic selection to obtain more resilient farm animals (i.e. that are minimally affected by disturbances or rapidly recover from them). The aims of this study were to: (i) estimate the genetic parameters of resilience indicator traits based on egg production data, (ii) assess whether these traits are genetically correlated in purebreds and crossbreds, and (iii) assess the genetic correlations of these traits with egg production (EP) as total number of eggs between 25 and 83 weeks. Purebred hens (33,825 from a White Leghorn (WA) line and 34,397 from a Rhode Island (BD) line were housed in individual cages, while crossbred hens were housed in collective cages of 6 to 8 paternal half-sibs (12,852 WA and 3898 BD crossbred groups, where the name of the group refers to the line used as the sire). Deviations of a hen’s weekly egg production from the average of the corresponding batch were calculated. Resilience indicator traits investigated were the natural logarithm of the variance (LNVAR), the skewness (SKEW), and the lag-one autocorrelation (AUTO-R) of these deviations.

Results

In both purebred lines, EP was estimated to be lowly heritable (WA: 0.11 and BD: 0.12). Resilience indicators were also estimated to be lowly heritable in both lines (LNVAR: 0.10 and 0.12, SKEW: 0.04 and 0.02, AUTO-R: 0.06 and 0.08 in WA and BD, respectively). In both crossbred groups, EP, AUTO-R, and SKEW were estimated to be less heritable than in purebreds (EP: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 ≤ 0.07; and resilience indicator traits: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 ≤ 0.03), while LNVAR had an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 estimate that was similar to or higher in crossbreds (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 ranged from 0.13 to 0.21) than in purebreds. In both purebreds and crossbreds, resilience indicator traits were estimated to have favorable genetic correlations with EP and between each other. For all traits and in both lines, estimates of genetic correlations between purebreds and crossbreds (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc) differed from 1 and ranged from 0.16 to 0.63.

Conclusions

These results show that selection for resilience based on EP data can be considered in breeding programs for layers. Genetic improvement of resilience in crossbreds can be achieved by using information on purebreds, but would be greatly enhanced by the integration of information on crossbreds in breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00716-8.

Genome-Wide Association Studies and Haplotype-Sharing Analysis Targeting the Egg Production Traits in Shaoxing Duck

Age at first egg (AFE) and egg number (EN) are economically important traits related to egg production, as they directly influence the benefits of the poultry industry, but the molecular genetic research that affects those traits in laying ducks is still sparse. Our objective was to identify the genomic regions and candidate genes associated with AFE, egg production at 43 weeks (EP43w), and egg production at 66 weeks (EP66w) in a Shaoxing duck population using genome-wide association studies (GWASs) and haplotype-sharing analysis. Single-nucleotide polymorphism (SNP)-based genetic parameter estimates showed that the heritability was 0.15, 0.20, and 0.22 for AFE, EP43w, and EP66w, respectively. Subsequently, three univariate GWASs for AFE, EP43w, and EP66w were carried out independently. Twenty-four SNPs located on chromosome 25 within a 0.01-Mb region that spans from 4.511 to 4.521 Mb were associated with AFE. There are two CIs that affect EP43w, i.e., twenty-five SNPs were in strong linkage disequilibrium region spanning from 3.186 to 3.247 Mb on chromosome 25, a region spanning from 4.442 to 4.446 Mb on chromosome 25, and two interesting genes, ACAD8 and THYN1, that may affect EP43w in laying ducks. There are also two CIs that affect EP66w, i.e., a 2.412-Mb region that spans from 127.497 to 129.910 Mb on chromosome 2 and a 0.355-Mb region that spans from 4.481 to 4.837 Mb on chromosome 29, and CA2 and GAMT may be the putative candidate genes. Our study also found some haplotypes significantly associated with these three traits based on haplotype-sharing analysis. Overall, this study was the first publication of GWAS on egg production in laying ducks, and our findings will be helpful to provide some candidate genes and haplotypes to improve egg production performance based on breeding in laying duck. Additionally, we learned from a method called bootstrap test to verify the reliability of a GWAS with small experimental samples that users can access at https://github.com/xuwenwu24/Bootstrap-test.

Genome-wide association studies for growth traits in broilers

BACKGROUND

The identification of markers and genes for growth traits may not only benefit for marker assist selection /genomic selection but also provide important information for understanding the genetic foundation of growth traits in broilers.

RESULTS

In the current study, we estimated the genetic parameters of eight growth traits in broilers and carried out the genome-wide association studies for these growth traits. A total of 113 QTNs discovered by multiple methods together, and some genes, including ACTA1, IGF2BP1, TAPT1, LDB2, PRKCA, TGFBR2, GLI3, SLC16A7, INHBA, BAMBI, APCDD1, GPR39, and GATA4, were identified as important candidate genes for rapid growth in broilers.

CONCLUSIONS

The results of this study will provide important information for understanding the genetic foundation of growth traits in broilers.

Genetics and breeding of a black-bone and blue eggshell chicken line. 2. Laying patterns and egg production in two consecutive generations

The BG line was originated from the cross between 2 Chinese indigenous breeds, Dongxiang blue eggshell, and Jiangshan black-bone, and has been bred to combine dark heavy black-bone body and high production of blue-shell eggs, into single dual-purpose line. Full-pedigree hens from 2 generations, G4 (n = 441) and G5 (n = 464), were reared in the same single-cage laying facility in 2019–2020 and 2020–2021, respectively. Starting from the first egg of each hen, its daily egg production was recorded until 300 days-of-age. Up to 7 "no-egg" days were considered normal laying breaks between clutches, whereas laying cessation of 8 or more days was considered Pause, and the laying pattern of each hen was assigned either with Pause or No-Pause. The other traits included PsDays: number of Pause days; AFE: age at first egg; EN300: eggs laid until 300 d; %L300: total laying rate (EN300/[300-AFE]); %Lnet: net laying rate (EN300/[300-AFE-PsDays]); ClLng: average clutch length; EW200 and EW300: average egg weight around 200 d and 300 d. Estimates of heritability (h²) of each trait, and phenotypic and genetic correlation between traits, were calculated in each generation using the animal model. Heritability estimates were calculated also by regressing the means of full-sib G5 hens on their G4 parents' means. Mean overall laying rate of all G4 hens was low (%L300 = 57%) because 53% of them had Pause in their laying pattern. In G5, incidence of Pause was higher (75%) due to a 3-wk cold stress, with mean %L300 = 54%. However, significant estimates of heritability and genetic correlations suggest that selection for low PsDays will reduce the incidence of Pause in BG hens and elevate the line's mean laying rate towards %L300 = 70%, as the No-Pause hens in G5. PsDays-free laying rate (%Lnet) was found to be highly correlated with the significantly heritable (h²≈0.4) clutch length (ClLng). Selection index combining the genetically independent low PsDays and high ClLng is expected to maximize egg production improvement in the BG line, and in similar populations derived from indigenous breeds.

Genome-Wide Association Study of Growth Performance and Immune Response to Newcastle Disease Virus of Indigenous Chicken in Rwanda

A chicken genome has several regions with quantitative trait loci (QTLs). However, replication and confirmation of QTL effects are required particularly in African chicken populations. This study identified single nucleotide polymorphisms (SNPs) and putative genes responsible for body weight (BW) and antibody response (AbR) to Newcastle disease (ND) in Rwanda indigenous chicken (IC) using genome-wide association studies (GWAS). Multiple testing was corrected using chromosomal false detection rates of 5 and 10% for significant and suggestive thresholds, respectively. BioMart data mining and variant effect predictor tools were used to annotate SNPs and candidate genes, respectively. A total of four significant SNPs (rs74098018, rs13792572, rs314702374, and rs14123335) significantly (p ≤ 7.6E-5) associated with BW were identified on chromosomes (CHRs) 8, 11, and 19. In the vicinity of these SNPs, four genes such as pre-B-cell leukaemia homeobox 1 (PBX1), GPATCH1, MPHOSPH6, and MRM1 were identified. Four other significant SNPs (rs314787954, rs13623466, rs13910430, and rs737507850) all located on chromosome 1 were strongly (p ≤ 7.6E-5) associated with chicken antibody response to ND. The closest genes to these four SNPs were cell division cycle 16 (CDC16), zinc finger, BED-type containing 1 (ZBED1), myxovirus (influenza virus) resistance 1 (MX1), and growth factor receptor bound protein 2 (GRB2) related adaptor protein 2 (GRAP2). Besides, other SNPs and genes suggestively (p ≤ 1.5E-5) associated with BW and antibody response to ND were reported. This work offers a useful entry point for the discovery of causative genes accountable for essential QTLs regulating BW and antibody response to ND traits. Results provide auspicious genes and SNP-based markers that can be used in the improvement of growth performance and ND resistance in IC populations based on gene-based and/or marker-assisted breeding selection.

Genome-Wide Association Study Using Whole-Genome Sequencing Identifies a Genomic Region on Chromosome 6 Associated With Comb Traits in Nandan-Yao Chicken

Comb traits have potential economic value in the breeding of indigenous chickens in China. Identifying and understanding relevant molecular markers for comb traits can be beneficial for genetic improvement. The purpose of this study was to utilize genome-wide association studies (GWAS) to detect promising loci and candidate genes related to comb traits, namely, comb thickness (CT), comb weight (CW), comb height, comb length (CL), and comb area. Genome-wide single-nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) in 300 Nandan-Yao chickens were detected using whole-genome sequencing. In total, we identified 134 SNPs and 25 INDELs that were strongly associated with the five comb traits. A remarkable region spanning from 29.6 to 31.4 Mb on chromosome 6 was found to be significantly associated with comb traits in both SNP- and INDEL-based GWAS. In this region, two lead SNPs (6:30,354,876 for CW and CT and 6:30,264,318 for CL) and one lead INDEL (a deletion from 30,376,404 to 30,376,405 bp for CL and CT) were identified. Additionally, two genes were identified as potential candidates for comb development. The nearby gene fibroblast growth factor receptor 2 (FGFR2)-associated with epithelial cell migration and proliferation-and the gene cytochrome b5 reductase 2 (CYB5R2)-identified on chromosome 5 from INDEL-based GWAS-are significantly correlated with collagen maturation. The findings of this study could provide promising genes and biomarkers to accelerate genetic improvement of comb development based on molecular marker-assisted breeding in Nandan-Yao chickens.

Novel genome reveals susceptibility of popular gamebird, the red-legged partridge (Alectoris rufa, Phasianidae), to climate change

We produced a high-quality de novo genome assembly of the red-legged partridge A. rufa, the first reference genome of its genus, by utilising novel 10× Chromium technology. The estimated genome size was 1.19 Gb with an overall genome heterozygosity of 0.0022; no runs of homozygosity were observed. In total, 21,589 protein coding genes were identified and assigned to 16,772 orthologs. Of these, 201 emerged as unique to Alectoris and were enriched for positive regulation of epithelial cell migration, viral genome integration and maturation. Using PSMC analysis, we inferred a major demographic decline commencing ~140,000 years ago, consistent with forest expansion and reduction of open habitats during the Eemian interglacial. Present-day populations exhibit the historically lowest genetic diversity. Besides implications for management and conservation, this genome also promises key insights into the physiology of these birds with a view to improving poultry husbandry practices.

Identification of candidate genomic regions for chicken egg number traits based on genome-wide association study

Background

Since the domestication of chicken, various breeds have been developed for food production, entertainment, and so on. Compared to indigenous chicken breeds which generally do not show elite production performance, commercial breeds or lines are selected intensely for meat or egg production. In the present study, in order to understand the molecular mechanisms underlying the dramatic differences of egg number between commercial egg-type chickens and indigenous chickens, we performed a genome-wide association study (GWAS) in a mixed linear model.

Results

We obtained 148 single nucleotide polymorphisms (SNPs) associated with egg number traits (57 significantly, 91 suggestively). Among them, 4 SNPs overlapped with previously reported quantitative trait loci (QTL), including 2 for egg production and 2 for reproductive traits. Furthermore, we identified 32 candidate genes based on the function of the screened genes. These genes were found to be mainly involved in regulating hormones, playing a role in the formation, growth, and development of follicles, and in the development of the reproductive system. Some genes such as NELL2 (neural EGFL like 2), KITLG (KIT ligand), GHRHR (Growth hormone releasing hormone receptor), NCOA1 (Nuclear receptor coactivator 1), ITPR1 (inositol 1, 4, 5-trisphosphate receptor type 1), GAMT (guanidinoacetate N-methyltransferase), and CAMK4 (calcium/calmodulin-dependent protein kinase IV) deserve our attention and further study since they have been reported to be closely related to egg production, egg number and reproductive traits. In addition, the most significant genomic region obtained in this study was located at 48.61–48.84 Mb on GGA5. In this region, we have repeatedly identified four genes, in which YY1 (YY1 transcription factor) and WDR25 (WD repeat domain 25) have been shown to be related to oocytes and reproductive tissues, respectively, which implies that this region may be a candidate region underlying egg number traits.

Conclusion

Our study utilized the genomic information from various chicken breeds or populations differed in the average annual egg number to understand the molecular genetic mechanisms involved in egg number traits. We identified a series of SNPs, candidate genes, or genomic regions that associated with egg number, which could help us in developing the egg production trait in chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07755-3.

Comparison of Selection Signatures between Korean Native and Commercial Chickens Using 600K SNP Array Data

Korean native chickens (KNCs) comprise an indigenous chicken breed of South Korea that was restored through a government project in the 1990s. The KNC population has not been developed well and has mostly been used to maintain purebred populations in the government research institution. We investigated the genetic features of the KNC population in a selection signal study for the efficient improvement of this breed. We used 600K single nucleotide polymorphism data sampled from 191 KNCs (NG, 38; NL, 29; NR, 52; NW, 39; and NY, 33) and 54 commercial chickens (Hy-line Brown, 10; Lohmann Brown, 10; Arbor Acres, 10; Cobb, 12; and Ross, 12). Haplotype phasing was performed using EAGLE software as the initial step for the primary data analysis. Pre-processed data were analyzed to detect selection signals using the ‘rehh’ package in R software. A few common signatures of selection were identified in KNCs. Most quantitative trait locus regions identified as candidate regions were associated with traits related to reproductive organs, eggshell characteristics, immunity, and organ development. Block patterns with high linkage disequilibrium values were observed for LPP, IGF11, LMNB2, ERBB4, GABRB2, NTM, APOO, PLOA1, CNTN1, NTSR1, DEF3, CELF1, and MEF2D genes, among regions with confirmed selection signals. NL and NW lines contained a considerable number of selective sweep regions related to broilers and layers, respectively. We recommend focusing on improving the egg and meat traits of KNC NL and NW lines, respectively, while improving multiple traits for the other lines.

Genomic prediction for growth using a low-density SNP panel in dromedary camels

For thousands of years, camels have produced meat, milk, and fiber in harsh desert conditions. For a sustainable development to provide protein resources from desert areas, it is necessary to pay attention to genetic improvement in camel breeding. By using genotyping-by-sequencing (GBS) method we produced over 14,500 genome wide markers to conduct a genome- wide association study (GWAS) for investigating the birth weight, daily gain, and body weight of 96 dromedaries in the Iranian central desert. A total of 99 SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.002). Genomic breeding values (GEBVs) were estimated with the BGLR package using (i) all 14,522 SNPs and (ii) the 99 SNPs by GWAS. Twenty-eight SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.001). Annotation of the genomic region (s) within ± 100 kb of the associated SNPs facilitated prediction of 36 candidate genes. The accuracy of GEBVs was more than 0.65 based on all 14,522 SNPs, but the regression coefficients for birth weight, daily gain, and body weight were 0.39, 0.20, and 0.23, respectively. Because of low sample size, the GEBVs were predicted using the associated SNPs from GWAS. The accuracy of GEBVs based on the 99 associated SNPs was 0.62, 0.82, and 0.57 for birth weight, daily gain, and body weight. This report is the first GWAS using GBS on dromedary camels and identifies markers associated with growth traits that could help to plan breeding program to genetic improvement. Further researches using larger sample size and collaboration of the camel farmers and more profound understanding will permit verification of the associated SNPs identified in this project. The preliminary results of study show that genomic selection could be the appropriate way to genetic improvement of body weight in dromedary camels, which is challenging due to a long generation interval, seasonal reproduction, and lack of records and pedigrees.

Changes in the Control of the Hypothalamic-Pituitary Gonadal Axis Across Three Differentially Selected Strains of Laying Hens (Gallus gallus domesticus)

Genetic selection for earlier sexual maturation and extended production cycles in laying hens has significantly improved reproductive efficiency. While limited emphasis has been placed on the underlying physiological changes, we hypothesize that modifications in the control of the hypothalamic-pituitary gonadal (HPG) axis have occurred. Thus, three strains of White leghorn derivatives were followed from hatch to 100 weeks of age (woa), including Lohmann LSL-lite (n = 120) as current commercial hens, heritage Shaver White leghorns (n = 100) as 2000s commercial equivalents, and Smoky Joe hens (n = 68) as 1960s commercial equivalents. Body weight (BW) and egg production were monitored, and blood samples were collected throughout to monitor estradiol (E₂) concentrations. Tissue samples were collected at 12, 17, 20, 25, 45, 60, 75, and 100 woa to capture changes in mRNA levels of key genes involved in the HPG axis and monitor ovarian follicular pools. All hens, regardless of strain, age or photoperiod laid their first egg within a 64-gram BW window and, as E₂ levels increased prior to photostimulation (PS) in Lohmann and Shaver hens, a metabolic trigger likely induced sexual maturation. However, increased levels of Opsin 5 (OPN5) were observed during the maturation period. Although an elevation in gonadotrophin-releasing hormone I (GnRH-I) mRNA levels was associated with early maturation, no changes in gonadotrophin-inhibitory hormone (GnIH) mRNA levels were observed. Nonetheless, a significant shift in pituitary sensitivity to GnRH was associated with maturation. Throughout the trial, Lohmann, Shaver, and Smoky Joe hens laid 515, 417, and 257 eggs, respectively (p < 0.0001). Results show that the extended laying persistency in Lohmann hens was supported by sustained pituitary sensitivity to GnRH-I, recurrent elevations in follicle-stimulating hormone (FSH) mRNA levels, and five cyclical elevations in E₂ levels. This was also associated with a consistently higher pool of small white ovarian follicles. In summary, our results demonstrate first that, regardless of photoperiodic cues, meeting a specific narrow body weight threshold is sufficient to initiate sexual maturation in Leghorn chicken derivatives. Furthermore, recurrent increases in E₂ and FSH may be the key to sustain extended laying period, allowing modern layers to double their reproductive capacity compared to their 1960s-counterparts.

Genome-Wide Association Study-Based Identification of SNPs and Haplotypes Associated With Goose Reproductive Performance and Egg Quality

Geese are one of the most economically important waterfowl. However, the low reproductive performance and egg quality of geese hinder the development of the goose industry. The identification and application of genetic markers may improve the accuracy of beneficial trait selection. To identify the genetic markers associated with goose reproductive performance and egg quality traits, we performed a genome-wide association study (GWAS) for body weight at birth (BBW), the number of eggs at 48 weeks of age (EN48), the number of eggs at 60 weeks of age (EN60) and egg yolk color (EYC). The GWAS acquired 2.896 Tb of raw sequencing data with an average depth of 12.44× and identified 9,279,339 SNPs. The results of GWAS showed that 26 SNPs were significantly associated with BBW, EN48, EN60, and EYC. Moreover, five of these SNPs significantly associated with EN48 and EN60 were in a haplotype block on chromosome 35 from 4,512,855 to 4,541,709 bp, oriented to TMEM161A and another five SNPs significantly correlated to EYC were constructed in haplotype block on chromosome 5 from 21,069,009 to 21,363,580, which annotated by TMEM161A, CALCR, TFPI2, and GLP1R. Those genes were enriched in epidermal growth factor-activated receptor activity, regulation of epidermal growth factor receptor signaling pathway. The SNPs, haplotype markers, and candidate genes identified in this study can be used to improve the accuracy of marker-assisted selection for the reproductive performance and egg quality traits of geese. In addition, the candidate genes significantly associated with these traits may provide a foundation for better understanding the mechanisms underlying reproduction and egg quality in geese.

Detection of CNV in the SH3RF2 gene and its effects on growth and carcass traits in chickens

BACKGROUND

The SH3RF2 gene is a protein-coding gene located in a quantitative trait locus associated with body weight, and its deletion has been shown to be positively associated with body weight in chickens.

RESULTS

In the present study, CNV in the SH3RF2 gene was detected in 4079 individuals from 17 populations, including the "Gushi ×Anka" F2 resource population and populations of Chinese native chickens, commercial layers, and commercial broilers. The F2 resource population was then used to investigate the genetic effects of the chicken SH3RF2 gene. The results showed that the local chickens and commercial layers were all homozygous for the wild-type allele. Deletion mutation individuals were detected in all of the commercial broiler breeds except Hubbard broiler. A total of, 798 individuals in the F2 resource group were used to analyze the effects of genotype (DD/ID/II) on chicken production traits. The results showed that CNV was associated with 2-, 6-, 10-, and 12-week body weight (P = 0.026, 0.042, 0.021 and 0.039 respectively) and significantly associated with 8-week breast bone length (P = 0.045). The mutation was significantly associated with 8-week body weight (P = 0.007) and 4-week breast bone length (P = 0.010). CNV was significantly associated with evisceration weight, leg muscle weight, carcass weight, breast muscle weight and gizzard weight (P = 0.032, 0.033, 0.045, 0.004 and 0.000, respectively).

CONCLUSIONS

CNV of the SH3RF2 gene contributed to variation in the growth and weight gain of chickens.

Genome-Wide Association Study of Body Weight Traits in Chinese Fine-Wool Sheep

Body weight is an important economic trait for sheep and it is vital for their successful production and breeding. Therefore, identifying the genomic regions and biological pathways that contribute to understanding variability in body weight traits is significant for selection purposes. In this study, the genome-wide associations of birth, weaning, yearling, and adult weights of 460 fine-wool sheep were determined using resequencing technology. The results showed that 113 single nucleotide polymorphisms (SNPs) reached the genome-wide significance levels for the four body weight traits and 30 genes were annotated effectively, including AADACL3, VGF, NPC1, and SERPINA12. The genes annotated by these SNPs significantly enriched 78 gene ontology terms and 25 signaling pathways, and were found to mainly participate in skeletal muscle development and lipid metabolism. These genes can be used as candidate genes for body weight in sheep, and provide useful information for the production and genomic selection of Chinese fine-wool sheep.