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

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.

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.

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.

Multiphasic mixed growth models for turkeys

Growth models are important for optimization of feed formulation and breeding programs in turkeys. The objectives of this study were 1) to compare sex and line differences in turkeys in parameter estimates of mono- and di-phasic Gompertz growth models, and 2) to evaluate mono and diphasic mixed Gompertz growth models to determine the variation in parameter estimates in a group of female line turkey toms. A total of 1,056 manually recorded weekly average body weight (BW) observations from male and female turkeys of a male and female line from weeks 1 to 20 were used for objective 1. Daily median values of automatically collected individual BW of female line turkey toms were used for objective 2 and random components associated with individual subject animals related to mature weight and/or timing of maximum gain during each phase were introduced in the Gompertz model. Growth curve shapes were different between male line toms, male line hens, female line toms, and female line hens (P < 0.001). However, inflection points were similar between male and female line toms and between male and female line hens (14.06 vs. 13.72 wk and 11.22 and 10.71 wk, respectively), while mature BW differed between lines by 6.49 and 3.81 kg for toms and hens, respectively. The normalized growth rate constant (growth rate constant corrected for mature weight) was around the same magnitude between male and female line toms (0.0031 vs. 0.0038, respectively), but slightly lower in male line hens compared to female line hens (0.0072 vs. 0.0091, respectively). Diphasic Gompertz models described growth better in all line × sex combinations compared to the monophasic models (P < 0.001) and mixed diphasic Gompertz models showed improved fit over mixed monophasic Gompertz models. The correlation structure of the random components identified that individuals with a higher mature weight had a later inflection point and lower growth rate coefficients. These results provide tools for improved breeding practices and a structure to evaluate the effects of dietary or environmental factors on growth trajectories.

Genetic parameters of feather corticosterone and fault bars and correlations with production traits in turkeys (Meleagris gallopavo)

Robustness can refer to an animal's ability to overcome perturbations. Intense selection for production traits in livestock has resulted in reduced robustness which has negative implications for livability as well as production. There is increasing emphasis on improving robustness through poultry breeding, which may involve identifying novel phenotypes that could be used in selection strategies. The hypothalamic-pituitary-adrenal (HPA) axis and associated hormones (e.g., corticosterone) participate in many metabolic processes that are related to robustness. Corticosterone can be measured non-invasively in feathers (FCORT) and reflects the average HPA axis activity over the feather growing period, however measurement is expensive and time consuming. Fault bars are visible feather deformities that may be related to HPA axis activity and may be a more feasible indicator trait. In this study, we estimated variance components for FCORT and fault bars in a population of purebred turkeys as well as their genetic and partial phenotypic correlations with other economically relevant traits including growth and efficiency, carcass yield, and meat quality. The estimated heritability for FCORT was 0.21 ± 0.07 and for the fault bar traits (presence, incidence, severity, and index) estimates ranged from 0.09 to 0.24. The genetic correlation of FCORT with breast weight, breast meat yield, fillet weight, and ultimate pH were estimated at -0.34 ± 0.21, -0.45 ± 0.23, -0.33 ± 0.24, and 0.32 ± 0.24, respectively. The phenotypic correlations of FCORT with breast weight, breast meat yield, fillet weight, drum weight, and walking ability were -0.16, -0.23, -0.18, 0.17, and 0.21, respectively. Some fault bar traits showed similar genetic correlations with breast weight, breast meat yield, and walking ability but the magnitude was lower than those with FCORT. While the dataset is limited and results should be interpreted with caution, this study indicates that selection for traits related to HPA axis activity is possible in domestic turkeys. Further research should focus on investigating the association of these traits with other robustness-related traits and how to potentially implement these traits in turkey breeding.

Genome-wide association study for reproduction-related traits in Chinese domestic goose

1. This study measured six reproduction traits in a Sichuan white goose population (209 individuals), including fertility, qualified egg rate, plasma concentrations of progesterone (P), follicle-stimulating hormone (FSH), prolactin (PRL) and oestrogen (E2).2. Whole-genome resequencing data from the same goose population (209 individuals) were used in a genome-wide association study (GWAS) utilising a mixed linear model to investigate the genes and genetic markers associated with reproduction traits. The frequency of the selected SNPs and haplotypes were determined using the Matrix-Assisted Laser Desorption Ionisation Time-Of-Flight Mass Spectrometry (MALDI-TOF MS) method.3. In total, 42 SNPs significantly associated with these traits were identified. A haplotype block was constructed based on five SNPs that were significantly associated with qualified egg rate, with individuals having the haplotype CCTTAAGGAA having the lowest qualified egg rate.4. In conclusion, these results provided potential markers for marker-assisted selection to improve goose reproductive performance and a basis for elucidating the genetics of goose reproduction.

Investigating inbreeding in the turkey (Meleagris gallopavo) genome

The detrimental effects of increased homozygosity due to inbreeding have prompted the development of methods to reduce inbreeding. The detection of runs of homozygosity (ROH), or contiguous stretches of homozygous marker genotypes, can be used to describe and quantify the level of inbreeding in an individual. The estimation of inbreeding coefficients can be calculated based on pedigree information, ROH, or the genomic relationship matrix. The aim of this study was to detect and describe ROH in the turkey genome and compare estimates of pedigree-based inbreeding coefficients (F_PED) with genomic-based inbreeding coefficients estimated from ROH (F_ROH) and the genomic relationship matrix (F_GRM). A total of 2,616,890 pedigree records were available. Of these records, 6,371 genotyped animals from three purebred turkey (Meleagris gallopavo) lines between 2013 and 2019 were available, and these were obtained using a dense single nucleotide polymorphism array (56,452 SNPs). The overall mean length of detected ROH was 2.87 ± 0.29 Mb with a mean number of 84.87 ± 8.79 ROH per animal. Short ROH with lengths of 1 to 2 Mb long were the most abundant throughout the genome. Mean ROH coverage differed greatly between chromosomes and lines. Considering inbreeding coefficient means across all lines, genomic derived inbreeding coefficients (F_ROH = 0.27; F_GRM = 0.32) were higher than coefficients estimated from pedigree records (F_PED = 0.14). Correlations between F_ROH and F_PED, F_ROH and F_GRM, and F_PED and F_GRM ranged between 0.19 to 0.31, 0.68 to 0.73, and 0.17 to 0.30, respectively. Additionally, correlations between F_ROH from different lengths and F_PED substantially increased with ROH length from -0.06 to 0.33. Results of the current research, including the distribution of ROH throughout the genome and ROH-derived inbreeding estimates, can provide a more comprehensive description of inbreeding in the turkey genome. This knowledge can be used to evaluate genetic diversity, a requirement for genetic improvement, and develop methods to minimize inbreeding in turkey breeding programs.

Accuracy of breeding values for production traits in turkeys (Meleagris gallopavo) using recursive models with or without genomics

Background

Knowledge about potential functional relationships among traits of interest offers a unique opportunity to understand causal mechanisms and to optimize breeding goals, management practices, and prediction accuracy. In this study, we inferred the phenotypic causal networks among five traits in a turkey population and assessed the effect of the use of such causal structures on the accuracy of predictions of breeding values.

Methods

Phenotypic data on feed conversion ratio, residual feed intake, body weight, breast meat yield, and walking score in addition to genotype data from a commercial breeding population were used. Causal links between the traits were detected using the inductive causation algorithm based on the joint distribution of genetic effects obtained from a standard Bayesian multiple trait model. Then, a structural equation model was implemented to infer the magnitude of causal structure coefficients among the phenotypes. Accuracies of predictions of breeding values derived using pedigree- and blending-based multiple trait models were compared to those obtained with the pedigree- and blending-based structural equation models.

Results

In contrast to the two unconditioned traits (i.e., feed conversion ratio and breast meat yield) in the causal structures, the three conditioned traits (i.e., residual feed intake, body weight, and walking score) showed noticeable changes in estimates of genetic and residual variances between the structural equation model and the multiple trait model. The analysis revealed interesting functional associations and indirect genetic effects. For example, the structural coefficient for the path from body weight to walking score indicated that a 1-unit genetic improvement in body weight is expected to result in a 0.27-unit decline in walking score. Both structural equation models outperformed their counterpart multiple trait models for the conditioned traits. Applying the causal structures led to an increase in accuracy of estimated breeding values of approximately 7, 6, and 20% for residual feed intake, body weight, and walking score, respectively, and different rankings of selection candidates for the conditioned traits.

Conclusions

Our results suggest that structural equation models can improve genetic selection decisions and increase the prediction accuracy of breeding values of selection candidates. The identified causal relationships between the studied traits should be carefully considered in future turkey breeding programs.

Discovering lethal alleles across the turkey genome using a transmission ratio distortion approach

Deviation from Mendelian inheritance expectations (transmission ratio distortion, TRD) has been observed in several species, including the mouse and humans. In this study, TRD was characterized in the turkey genome using both allelic (specific- and unspecific-parent TRD) and genotypic (additive- and dominance-TRD) parameterizations within a Bayesian framework. In this study, we evaluated TRD for 23 243 genotyped Turkeys across 56 393 autosomal SNPs. The analyses included 500 sires, 2013 dams and 11 047 offspring (trios). Three different haplotype sliding windows of 4, 10 and 20 SNPs were used across the autosomal chromosomes. Based on the genotypic parameterizations, 14 haplotypes showed additive and dominance TRD effects highlighting regions with a recessive TRD pattern. In contrast, the allelic model uncovered 12 haplotype alleles with the allelic TRD pattern which showed an underrepresentation of heterozygous offspring in addition to the absence of homozygous animals. For regions with the allelic pattern, only one particular region showed a parent-specific TRD where the penetrance was high via the dam, but low via the sire. The gene set analysis uncovered several gene ontology functional terms, Reactome pathways and several Medical Subject Headings that showed significant enrichment of genes associated with TRD. Many of these gene ontology functional terms (e.g. mitotic spindle assembly checkpoint, DRM complex and Aneuploidy), Reactome pathways (e.g. Mismatch repair) and Medical Subject Headings (e.g. Adenosine monophosphate) are known to be related to fertility, embryo development and lethality. The results of this study revealed potential novel candidate lethal haplotypes, functional terms and pathways that may enhance breeding programs in Turkeys through reducing mortality and improving reproduction rate.