Genotype by Environment Interactions (G*E) of Chickens Tested in Ethiopia Using Body Weight as a Performance Trait.

Ethiopia is a developing nation that could highly benefit from securing food via improved smallholder poultry farming. To support farmer and breeding decisions regarding which chicken strain to use in which Ethiopian environment, G*E analyses for body weight (BW) of growing male and female chickens were conducted. Research questions were (1) if a G*E is present for BW and (2) which strain performs best in which environment in terms of predicted BW. Analyses were performed using predicted BW at four different ages (90, 120, 150, and 180 days) of five strains (Horro, Koekoek, Kuroiler, Sasso-Rhode Island Red (S-RIR), and Sasso) tested in five Ethiopian regions (Addis Ababa, Amhara, Oromia, South Region, and Tigray) that are part of three Agro-Ecological Zones (AEZ) (cool humid, cool sub-humid, and warm semi-arid). The indigenous Horro strain was used as a control group to compare four other introduced tropically adapted strains. The dataset consisted of 999 female and 989 male farm-average BW measurements. G*E was strongly present (p < 0.001) for all combinations of strain and region analyzed. In line with previous research, Sasso was shown to have the highest predicted BW, especially at an early age, followed by Kuroiler. Horro had the lowest predicted BW at most ages and in most regions, potentially due to its young breeding program. The highest predicted BW were observed in Tigray, Oromia, and Amhara regions, which are in the main part of the cool sub-humid AEZ.

Comparison of linkage disequilibrium estimated from genotypes versus haplotypes for crossbred populations.

BACKGROUND: Linkage disequilibrium (LD) is commonly measured based on the squared coefficient of correlation [Formula: see text] between the alleles at two loci that are carried by haplotypes. LD can also be estimated as the [Formula: see text] between unphased genotype dosage at two loci when the allele frequencies and inbreeding coefficients at both loci are identical for the parental lines. Here, we investigated whether [Formula: see text] for a crossbred population (F1) can be estimated using genotype data. The parental lines of the crossbred (F1) can be purebred or crossbred. METHODS: We approached this by first showing that inbreeding coefficients for an F1 crossbred population are negative, and typically differ in size between loci. Then, we proved that the expected [Formula: see text] computed from unphased genotype data is expected to be identical to the [Formula: see text] computed from haplotype data for an F1 crossbred population, regardless of the inbreeding coefficients at the two loci. Finally, we investigated the bias and precision of the [Formula: see text] estimated using unphased genotype versus haplotype data in stochastic simulation. RESULTS: Our findings show that estimates of [Formula: see text] based on haplotype and unphased genotype data are both unbiased for different combinations of allele frequencies, sample sizes (900, 1800, and 2700), and levels of LD. In general, for any allele frequency combination and [Formula: see text] value scenarios considered, and for both methods to estimate [Formula: see text], the precision of the estimates increased, and the bias of the estimates decreased as sample size increased, indicating that both estimators are consistent. For a given scenario, the [Formula: see text] estimates using haplotype data were more precise and less biased using haplotype data than using unphased genotype data. As sample size increased, the difference in precision and biasedness between the [Formula: see text] estimates using haplotype data and unphased genotype data decreased. CONCLUSIONS: Our theoretical derivations showed that estimates of LD between loci based on unphased genotypes and haplotypes in F1 crossbreds have identical expectations. Based on our simulation results, we conclude that the LD for an F1 crossbred population can be accurately estimated from unphased genotype data. The results also apply for other crosses (F2, F3, Fn, BC1, BC2, and BCn), as long as (selected) individuals from the two parental lines mate randomly.

Species and Phenotypic Distribution Models Reveal Population Differentiation in Ethiopian Indigenous Chickens.

Smallholder poultry production dominated by indigenous chickens is an important source of livelihoods for most rural households in Ethiopia. The long history of domestication and the presence of diverse agroecologies in Ethiopia create unique opportunities to study the effect of environmental selective pressures. Species distribution models (SDMs) and Phenotypic distribution models (PDMs) can be applied to investigate the relationship between environmental variation and phenotypic differentiation in wild animals and domestic populations. In the present study we used SDMs and PDMs to detect environmental variables related with habitat suitability and phenotypic differentiation among nondescript Ethiopian indigenous chicken populations. 34 environmental variables (climatic, soil, and vegetation) and 19 quantitative traits were analyzed for 513 adult chickens from 26 populations. To have high variation in the dataset for phenotypic and ecological parameters, animals were sampled from four spatial gradients (each represented by six to seven populations), located in different climatic zones and geographies. Three different ecotypes are proposed based on correlation test between habitat suitability maps and phenotypic clustering of sample populations. These specific ecotypes show phenotypic differentiation, likely in response to environmental selective pressures. Nine environmental variables with the highest contribution to habitat suitability are identified. The relationship between quantitative traits and a few of the environmental variables associated with habitat suitability is non-linear. Our results highlight the benefits of integrating species and phenotypic distribution modeling approaches in characterization of livestock populations, delineation of suitable habitats for specific breeds, and understanding of the relationship between ecological variables and quantitative traits, and underlying evolutionary processes.

An evaluation of inbreeding measures using a whole-genome sequenced cattle pedigree.

The estimation of the inbreeding coefficient (F) is essential for the study of inbreeding depression (ID) or for the management of populations under conservation. Several methods have been proposed to estimate the realized F using genetic markers, but it remains unclear which one should be used. Here we used whole-genome sequence data for 245 individuals from a Holstein cattle pedigree to empirically evaluate which estimators best capture homozygosity at variants causing ID, such as rare deleterious alleles or loci presenting heterozygote advantage and segregating at intermediate frequency. Estimators relying on the correlation between uniting gametes (FUNI) or on the genomic relationships (FGRM) presented the highest correlations with these variants. However, homozygosity at rare alleles remained poorly captured. A second group of estimators relying on excess homozygosity (FHOM), homozygous-by-descent segments (FHBD), runs-of-homozygosity (FROH) or on the known genealogy (FPED) was better at capturing whole-genome homozygosity, reflecting the consequences of inbreeding on all variants, and for young alleles with low to moderate frequencies (0.10 < . < 0.25). The results indicate that FUNI and FGRM might present a stronger association with ID. However, the situation might be different when recessive deleterious alleles reach higher frequencies, such as in populations with a small effective population size. For locus-specific inbreeding measures or at low marker density, the ranking of the methods can also change as FHBD makes better use of the information from neighboring markers. Finally, we confirmed that genomic measures are in general superior to pedigree-based estimates. In particular, FPED was uncorrelated with locus-specific homozygosity.

Effect of genomic selection and genotyping strategy on estimation of variance components in animal models using different relationship matrices.

BACKGROUND: The traditional way to estimate variance components (VC) is based on the animal model using a pedigree-based relationship matrix (A) (A-AM). After genomic selection was introduced into breeding programs, it was anticipated that VC estimates from A-AM would be biased because the effect of selection based on genomic information is not captured. The single-step method (H-AM), which uses an H matrix as (co)variance matrix, can be used as an alternative to estimate VC. Here, we compared VC estimates from A-AM and H-AM and investigated the effect of genomic selection, genotyping strategy and genotyping proportion on the estimation of VC from the two methods, by analyzing a dataset from a commercial broiler line and a simulated dataset that mimicked the broiler population. RESULTS: VC estimates from H-AM were severely overestimated with a high proportion of selective genotyping, and overestimation increased as proportion of genotyping increased in the analysis of both commercial and simulated data. This bias in H-AM estimates arises when selective genotyping is used to construct the H-matrix, regardless of whether selective genotyping is applied or not in the selection process. For simulated populations under genomic selection, estimates of genetic variance from A-AM were also significantly overestimated when the effect of genomic selection was strong. Our results suggest that VC estimates from H-AM under random genotyping have the expected values. Predicted breeding values from H-AM were inflated when VC estimates were biased, and inflation differed between genotyped and ungenotyped animals, which can lead to suboptimal selection decisions. CONCLUSIONS: We conclude that VC estimates from H-AM are biased with selective genotyping, but are close to expected values with random genotyping.VC estimates from A-AM in populations under genomic selection are also biased but to a much lesser degree. Therefore, we recommend the use of H-AM with random genotyping to estimate VC for populations under genomic selection. Our results indicate that it is still possible to use selective genotyping in selection, but then VC estimation should avoid the use of genotypes from one side only of the distribution of phenotypes. Hence, a dual genotyping strategy may be needed to address both selection and VC estimation.

Benefits of testing in both bio-secure and production environments in genomic selection breeding programs for commercial broiler chicken.

BACKGROUND: A breeding program for commercial broiler chicken that is carried out under strict biosecure conditions can show reduced genetic gain due to genotype by environment interactions (G × E) between bio-secure (B) and commercial production (C) environments. Accuracy of phenotype-based best linear unbiased prediction of breeding values of selection candidates using sib-testing in C is low. Genomic prediction based on dense genetic markers may improve accuracy of selection. Stochastic simulation was used to explore the benefits of genomic selection in breeding schemes for broiler chicken that include birds in both B and C for assessment of phenotype. RESULTS: When genetic correlations ([Formula: see text]) between traits measured in B and C were equal to 0.5 and 0.7, breeding schemes with 15, 30 and 45% of birds assessed in C resulted in higher genetic gain for performance in C compared to those without birds in C. The optimal proportion of birds phenotyped in C for genetic gain was 30%. When the proportion of birds in C was optimal and genotyping effort was limited, allocating 30% of the genotyping effort to birds in C was also the optimal genotyping strategy for genetic gain. When [Formula: see text] was equal to 0.9, genetic gain for performance in C was not improved with birds in C compared to schemes without birds in C. Increasing the heritability of traits assessed in C increased genetic gain significantly. Rates of inbreeding decreased when the proportion of birds in C increased because of a lower selection intensity among birds retained in B and a reduction in the probability of co-selecting close relatives. CONCLUSIONS: If G × E interactions ([Formula: see text] of 0.5 and 0.7) are strong, a genomic selection scheme in which 30% of the birds hatched are phenotyped in C has larger genetic gain for performance in C compared to phenotyping all birds in B. Rates of inbreeding decreased as the proportion of birds moved to C increased from 15 to 45%.

Genomic prediction of survival time in a population of brown laying hens showing cannibalistic behavior.

BACKGROUND: Mortality due to cannibalism causes both economic and welfare problems in laying hens. To limit mortality due to cannibalism, laying hens are often beak-trimmed, which is undesirable for animal welfare reasons. Genetic selection is an alternative strategy to increase survival and is more efficient by taking heritable variation that originates from social interactions into account, which are modelled as the so-called indirect genetic effects (IGE). Despite the considerable heritable variation in survival time due to IGE, genetic improvement of survival time in laying hens is still challenging because the detected heritable variation of the trait with IGE is still limited, ranging from 0.06 to 0.26, and individuals that are still alive at the end of the recording period are censored. Furthermore, survival time records are available late in life and only on females. To cope with these challenges, we tested the hypothesis that genomic prediction increases the accuracy of estimated breeding values (EBV) compared to parental average EBV, and increases response to selection for survival time compared to a traditional breeding scheme. We tested this hypothesis in two lines of brown layers with intact beaks, which show cannibalism, and also the hypothesis that the rate of inbreeding per year is lower for genomic selection than for the traditional breeding scheme. RESULTS AND DISCUSSION: The standard deviation of genomic prediction EBV for survival time was around 22 days for both lines, indicating good prospects for selection against mortality in laying hens with intact beaks. Genomic prediction increased the accuracy of the EBV by 35 and 32 % compared to the parent average EBV for the two lines. At the current reference population size, predicted response to selection was 91 % higher when using genomic selection than with the traditional breeding scheme, as a result of a shorter generation interval in males and greater accuracy of selection in females. The predicted rate of inbreeding per generation with truncation selection was substantially lower for genomic selection than for the traditional breeding scheme for both lines. CONCLUSIONS: Genomic selection for socially affected traits is a promising tool for the improvement of survival time in laying hens with intact beaks.

Indirect genetic effects contribute substantially to heritable variation in aggression-related traits in group-housed mink (Neovison vison).

BACKGROUND: Since the recommendations on group housing of mink (Neovison vison) were adopted by the Council of Europe in 1999, it has become common in mink production in Europe. Group housing is advantageous from a production perspective, but can lead to aggression between animals and thus raises a welfare issue. Bite marks on the animals are an indicator of this aggressive behaviour and thus selection against frequency of bite marks should reduce aggression and improve animal welfare. Bite marks on one individual reflect the aggression of its group members, which means that the number of bite marks carried by one individual depends on the behaviour of other individuals and that it may have a genetic basis. Thus, for a successful breeding strategy it could be crucial to consider both direct (DGE) and indirect (IGE) genetic effects on this trait. However, to date no study has investigated the genetic basis of bite marks in mink. RESULT AND DISCUSSION: A model that included DGE and IGE fitted the data significantly better than a model with DGE only, and IGE contributed a substantial proportion of the heritable variation available for response to selection. In the model with IGE, the total heritable variation expressed as the proportion of phenotypic variance (T2) was six times greater than classical heritability (h2). For instance, for total bite marks, T2 was equal to 0.61, while h2 was equal to 0.10. The genetic correlation between direct and indirect effects ranged from 0.55 for neck bite marks to 0.99 for tail bite marks. This positive correlation suggests that mink have a tendency to fight in a reciprocal way (giving and receiving bites) and thus, a genotype that confers a tendency to bite other individuals can also cause its bearer to receive more bites. CONCLUSION: Both direct and indirect genetic effects contribute to variation in number of bite marks in group-housed mink. Thus, a genetic selection design that includes both direct genetic and indirect genetic effects could reduce the frequency of bite marks and probably aggression behaviour in group-housed mink.