Effectiveness of DArTseq markers application in genetic diversity and population structure of indigenous chickens in Eastern Province of Rwanda

BACKGROUND

The application of biotechnologies which make use of genetic markers in chicken breeding is developing rapidly. Diversity Array Technology (DArT) is one of the current Genotyping-By-Sequencing techniques allowing the discovery of whole genome sequencing. In livestock, DArT has been applied in cattle, sheep, and horses. Currently, there is no study on the application of DArT markers in chickens. The aim was to study the effectiveness of DArTSeq markers in the genetic diversity and population structure of indigenous chickens (IC) and SASSO in the Eastern Province of Rwanda.

METHODS

In total 87 blood samples were randomly collected from 37 males and 40 females of indigenous chickens and 10 females of SASSO chickens purposively selected from 5 sites located in two districts of the Eastern Province of Rwanda. Genotyping by Sequencing (GBS) using DArTseq technology was employed. This involved the complexity reduction method through digestion of genomic DNA and ligation of barcoded adapters followed by PCR amplification of adapter-ligated fragments.

RESULTS

From 45,677 DArTseq SNPs and 25,444 SilicoDArTs generated, only 8,715 and 6,817 respectively remained for further analysis after quality control. The average call rates observed, 0.99 and 0.98 for DArTseq SNPs and SilicoDArTs respectively were quite similar. The polymorphic information content (PIC) from SilicoDArTs (0.33) was higher than that from DArTseq SNPs (0.22). DArTseq SNPs and SilicoDArTs had 34.4% and 34% of the loci respectively mapped on chromosome 1. DArTseq SNPs revealed distance averages of 0.17 and 0.15 within IC and SASSO chickens respectively while the respective averages observed with SilicoDArTs were 0.42 and 0.36. The average genetic distance between IC and SASSO chickens was moderate for SilicoDArTs (0.120) compared to that of DArTseq SNPs (0.048). The PCoA and population structure clustered the chicken samples into two subpopulations (1 and 2); 1 is composed of IC and 2 by SASSO chickens. An admixture was observed in subpopulation 2 with 12 chickens from subpopulation 1.

CONCLUSIONS

The application of DArTseq markers have been proven to be effective and efficient for genetic relationship between IC and separated IC from exotic breed used which indicate their suitability in genomic studies. However, further studies using all chicken genetic resources available and large big sample sizes are required.

Modeling Environmental Conditions in Poultry Production: Computational Fluid Dynamics Approach

In recent years, computational fluid dynamics (CFD) has become increasingly important and has proven to be an effective method for assessing environmental conditions in poultry houses. CFD offers simplicity, efficiency, and rapidity in assessing and optimizing poultry house environments, thereby fueling greater interest in its application. This article aims to facilitate researchers in their search for relevant CFD studies in poultry housing environmental conditions by providing an in-depth review of the latest advancements in this field. It has been found that CFD has been widely employed to study and analyze various aspects of poultry house ventilation and air quality under the following five main headings: inlet and fan configuration, ventilation system design, air temperature-humidity distribution, airflow distribution, and particle matter and gas emission. The most commonly used turbulence models in poultry buildings are the standard k-ε, renormalization group (RNG) k-ε, and realizable k-ε models. Additionally, this article presents key solutions with a summary and visualization of fundamental approaches employed in addressing path planning problems within the CFD process. Furthermore, potential challenges, such as data acquisition, validation, computational resource requirements, meshing, and the selection of a proper turbulence model, are discussed, and avenues for future research (the integration of machine learning, building information modeling, and feedback control systems with CFD) are explored.

A Single-Step Genome-Wide Association Study for Semen Traits of Egyptian Buffalo Bulls

The present study aimed to contribute to the limited research on buffalo (Bubalus bubalis) semen traits by incorporating genomic data. A total of 8465 ejaculates were collected. The genotyping procedure was conducted using the Axiom® Buffalo Genotyping 90 K array designed by the Affymetrix Expert Design Program. After conducting a quality assessment, we utilized 67,282 SNPs genotyped in 192 animals. We identified several genomic loci explaining high genetic variance by employing single-step genomic evaluation. The aforementioned regions were located on buffalo chromosomes no. 3, 4, 6, 7, 14, 16, 20, 22, and the X-chromosome. The X-chromosome exhibited substantial influence, accounting for 4.18, 4.59, 5.16, 5.19, and 4.31% of the genomic variance for ejaculate volume, mass motility, livability, abnormality, and concentration, respectively. In the examined genomic regions, we identified five novel candidate genes linked to male fertility and spermatogenesis, four in the X-chromosome and one in chromosome no. 16. Additional extensive research with larger sample sizes and datasets is imperative to validate these findings and evaluate their applicability for genomic selection.

OMICs approaches and technologies for understanding low-high feed efficiency traits in chicken: implication to breeding

In poultry production, there has been a trend of continuous increase in cost of feed ingredients which represents the major proportion of the production costs. Feed costs can be reduced by improving feed efficiency traits which increase the possibility of using various indigestible feed sources and decrease the environmental impact of the enhanced poultry production. Therefore, feed efficiency has been used as one of the most important economic traits of selection in the breeding program of chickens. Recently, many OMICs experimental studies have been designed to characterize biological differences between the high and low feed efficiency chicken phenotypes. Biological complexity cannot be fully captured by main individual OMICs such as genomics, transcriptomics, proteomics and metabolomics. Therefore, researchers have combined multiple assays from the same set of samples to create multi-OMICs datasets. OMICs findings are crucial in improving existing approaches to poultry breeding. The current review aimed to highlight the components of feed efficiency and general OMICs approaches and technologies. Besides, individual and multi-OMICs based understanding of chicken feed efficiency traits and the application of the acquired knowledge in the chicken breeding program were addressed.

First large-scale genomic prediction in the honey bee

Genomic selection has increased genetic gain in several livestock species, but due to the complicated genetics and reproduction biology not yet in honey bees. Recently, 2970 queens were genotyped to gather a reference population. For the application of genomic selection in honey bees, this study analyzes the accuracy and bias of pedigree-based and genomic breeding values for honey yield, three workability traits, and two traits for resistance against the parasite Varroa destructor. For breeding value estimation, we use a honey bee-specific model with maternal and direct effects, to account for the contributions of the workers and the queen of a colony to the phenotypes. We conducted a validation for the last generation and a five-fold cross-validation. In the validation for the last generation, the accuracy of pedigree-based estimated breeding values was 0.12 for honey yield, and ranged from 0.42 to 0.61 for the workability traits. The inclusion of genomic marker data improved these accuracies to 0.23 for honey yield, and a range from 0.44 to 0.65 for the workability traits. The inclusion of genomic data did not improve the accuracy of the disease-related traits. Traits with high heritability for maternal effects compared to the heritability for direct effects showed the most promising results. For all traits except the Varroa resistance traits, the bias with genomic methods was on a similar level compared to the bias with pedigree-based BLUP. The results show that genomic selection can successfully be applied to honey bees.

Computational Fluid Dynamics Modeling of a Broiler House Microclimate in Summer and Winter

Simple Summary

Microclimate conditions in broiler housing are significant for maximizing poultry production and ensuring the welfare of the birds. In the present study, we modeled summer and winter microclimates in a mechanically ventilated broiler house. Validation of the simulated values was accomplished through comparison to field measurements. In visual simulations, the results were used to reconstruct microclimate conditions such as stagnant and stress zones of broiler houses. In conclusion, simulation techniques can be used as an alternative method for analyzing poultry house indoor environments.

Abstract

Appropriate microclimate conditions in broiler housing are critical for optimizing poultry production and ensuring the health and welfare of the birds. In this study, spatial variabilities of the microclimate in summer and winter seasons in a mechanically ventilated broiler house were modeled using the computational fluid dynamics (CFD) technique. Field measurements of temperature, relative humidity, and airspeeds were conducted in the house to compare the simulated results. The study identified two problems of high temperature in summer, which could result in bird heat stress and stagnant zones in winter, and simulated possible alternative solutions. In summer, if an evaporative cooling pad system was used, a decrease in temperature of approximately 3 °C could be achieved when the mean air temperature rose above 25 °C in the house. In winter, adding four 500-mm circulation fans of 20-m spacing inside the house could eliminate the accumulation of hot and humid air in the stagnant zones in the house. This study demonstrated that CFD is a valuable tool for adequate heating, ventilation, and air conditioning system design in poultry buildings.

Simulation studies to optimize genomic selection in honey bees

Background

With the completion of a single nucleotide polymorphism (SNP) chip for honey bees, the technical basis of genomic selection is laid. However, for its application in practice, methods to estimate genomic breeding values need to be adapted to the specificities of the genetics and breeding infrastructure of this species. Drone-producing queens (DPQ) are used for mating control, and usually, they head non-phenotyped colonies that will be placed on mating stations. Breeding queens (BQ) head colonies that are intended to be phenotyped and used to produce new queens. Our aim was to evaluate different breeding program designs for the initiation of genomic selection in honey bees.

Methods

Stochastic simulations were conducted to evaluate the quality of the estimated breeding values. We developed a variation of the genomic relationship matrix to include genotypes of DPQ and tested different sizes of the reference population. The results were used to estimate genetic gain in the initial selection cycle of a genomic breeding program. This program was run over six years, and different numbers of genotyped queens per year were considered. Resources could be allocated to increase the reference population, or to perform genomic preselection of BQ and/or DPQ.

Results

Including the genotypes of 5000 phenotyped BQ increased the accuracy of predictions of breeding values by up to 173%, depending on the size of the reference population and the trait considered. To initiate a breeding program, genotyping a minimum number of 1000 queens per year is required. In this case, genetic gain was highest when genomic preselection of DPQ was coupled with the genotyping of 10–20% of the phenotyped BQ. For maximum genetic gain per used genotype, more than 2500 genotyped queens per year and preselection of all BQ and DPQ are required.

Conclusions

This study shows that the first priority in a breeding program is to genotype phenotyped BQ to obtain a sufficiently large reference population, which allows successful genomic preselection of queens. To maximize genetic gain, DPQ should be preselected, and their genotypes included in the genomic relationship matrix. We suggest, that the developed methods for genomic prediction are suitable for implementation in genomic honey bee breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00654-x.

Analyses of free amino acid and taste sensor traits in egg albumen and yolk revealed potential of value-added eggs in chickens

To create high-quality eggs by using different breed and feed materials, we investigated free amino acid contents and taste sensor traits using two chicken breeds (Rhode Island Red; RIR and Australorp; AUS) fed two feeds (mixed and fermented feeds). Two-way ANOVA revealed significant breed and feed main and interaction effects on albumen bitterness and a significant interaction effect on yolk bitterness. Albumen from RIR fed mixed feed and AUS fed fermented feed was higher bitterness, whereas yolk from those groups was lower bitterness. Significant breed effects were detected in four albumen amino acid traits (His, Met, Ile, and Lys) and a yolk His, whereas significant feed effects were found in 15 albumen amino acid traits (Asp, Glu, Ser, His, Gly, Thr, Ala, Tyr, Val, Met, Trp, Ile, Leu, Lys, and Pro) and a yolk cystine trait. Compared to albumen amino acids, yolk amino acids had limited effects by breed and feed. The present results suggest that genetic and nutritional factors can alter not only amino acid contents but also sensor values of bitterness, indicating that selecting the combination of breed and feed enable us to make amino acids enriched and taste added designer eggs in future.

Comparative Transcriptome Profiling of mRNA and lncRNA of Ovaries in High and Low Egg Production Performance in Domestic Pigeons (Columba livia)

Egg production performance is one of the most important economic traits in pigeon industry. However, little is known regarding how egg production performance is regulated by long non-coding RNAs (lncRNAs) in pigeons. To evaluate the lncRNAs and mRNAs in ovaries associated with egg production performance in domestic pigeons, high-throughput RNA sequencing of ovaries between high and low egg production performance groups were performed and analyzed in this study. A total of 34,346 mRNAs and 24,601 lncRNAs were identified, including 14,525 known lncRNAs and 10,076 novel lncRNAs, of which 811 mRNAs and 148 lncRNAs (P < 0.05) were significantly differentially expressed (DE) between the groups of high and low egg production performance. GO and KEGG annotation analysis indicated that the target genes of DE lncRNAs and DE mRNAs were related to cell differentiation, ATP binding and methylation. Moreover, we found that FOXK2, a target gene of lncRNA MSTRG.7894.4, was involved in regulating estrogen receptors. Our study provided a catalog of lncRNAs and mRNAs associated with egg production performance, and they deserve further study to deepen the understanding of biological processes in the ovaries of pigeons.

Genome-wide association mapping and accuracy of predictions for amoebic gill disease in Atlantic salmon (Salmo salar)

Amoebic gill disease (AGD) is a parasitic disease caused by the amoeba Paramoeba perurans, which colonizes the gill tissues and causes distress for the host. AGD can cause high morbidity and mortalities in salmonid and non-salmonid fish species. To understand the genetic basis of AGD and improve health status of farmed A. salmon, a population of ~ 6,100 individuals belonging to 150 full-sib families was monitored for development of AGD in the sea of Ireland. The population was followed for two rounds of AGD infections, and fish were gill scored to identify severity of disease in first (N = 3,663) and the second (N = 3,511) infection with freshwater treatment after the first gill-scoring. A subset of this gill-scored population (N = 1,141) from 119 full-sib families were genotyped with 57,184 SNPs using custom-made Affymetrix SNP-chip. GWAS analyses were performed which resulted in five significantly associated SNP variants distributed over chromosome 1, 2 and 5. Three candidate genes; c4, tnxb and slc44a4 were found within QTL region of chromosome 2. The tnxb and c4 genes are known to be a part of innate immune system, and may play a role in resistance to AGD. The gain in prediction accuracy obtained by involving genomic information was 9–17% higher than using traditional pedigree information.

Genetic Analysis of a Commercial Egg Laying Line Challenged With Newcastle Disease Virus

In low income countries, chickens play a vital role in daily life. They provide a critical source of protein through egg production and meat. Newcastle disease, caused by avian paramyxovirus type 1, has been ranked as the most devastating disease for scavenging chickens in Africa and Asia. High mortality among flocks infected with velogenic strains leads to a devastating loss of dietary protein and buying power for rural households. Improving the genetic resistance of chickens to Newcastle Disease virus (NDV), in addition to vaccination, is a practical target for improvement of poultry production in low income countries. Because response to NDV has a component of genetic control, it can be influenced through selective breeding. Adding genomic information to a breeding program can increase the amount of genetic progress per generation. In this study, we challenged a commercial egg-laying line with a lentogenic strain of NDV, measured phenotypic responses, collected genotypes, and associated genotypes with phenotypes. Collected phenotypes included viral load at 2 and 6 days post-infection (dpi), antibody levels pre-challenge and 10 dpi, and growth rates pre- and post-challenge. Six suggestive QTL associated with response to NDV and/or growth were identified, including novel and known QTL confirming previously reported associations with related traits. Additionally, previous RNA-seq analysis provided support for several of the genes located in or near the identified QTL. Considering the trend of negative genetic correlation between antibody and Newcastle Disease tolerance (growth under disease) and estimates of moderate to high heritability, we provide evidence that these NDV response traits can be influenced through selective breeding. Producing chickens that perform favorably in challenging environments will ultimately increase the supply of quality protein for human consumption.

Genome-wide association studies to identify quantitative trait loci affecting milk production traits in water buffalo

Water buffalo is the second largest resource of milk supply around the world, and it is well known for its distinctive milk quality in terms of fat, protein, lactose, vitamin, and mineral contents. Understanding the genetic architecture of milk production traits is important for future improvement by the buffalo breeding industry. The advance of genome-wide association studies (GWAS) provides an opportunity to identify potential genetic variants affecting important economical traits. In the present study, GWAS was performed for 489 buffaloes with 1,424 lactation records using the 90K Affymetrix Buffalo SNP Array (Affymetrix/Thermo Fisher Scientific, Santa Clara, CA). Collectively, 4 candidate single nucleotide polymorphisms (SNP) in 2 genomic regions were found to associate with buffalo milk production traits. One region affecting milk fat and protein percentage was located on the equivalent of Bos taurus autosome (BTA)3, spanning 43.3 to 43.8 Mb, which harbored the most likely candidate genes MFSD14A, SLC35A3, and PALMD. The other region on the equivalent of BTA14 at 66.5 to 67.0 Mb contained candidate genes RGS22 and VPS13B and influenced buffalo total milk yield, fat yield, and protein yield. Interestingly, both of the regions were reported to have quantitative trait loci affecting milk performance in dairy cattle. Furthermore, we suggest that buffaloes with the C allele at AX-85148558 and AX-85073877 loci and the G allele at AX-85106096 locus can be selected to improve milk fat yield in this buffalo-breeding program. Meanwhile, the G allele at AX-85063131 locus can be used as the favorable allele for improving milk protein percentage. Genomic prediction showed that the reliability of genomic estimated breeding values (GEBV) of 6 milk production traits ranged from 0.06 to 0.22, and the correlation between estimated breeding values and GEBV ranged from 0.23 to 0.35. These findings provide useful information to understand the genetic basis of buffalo milk properties and may play a role in accelerating buffalo breeding programs using genomic approaches.

Genetic aspects of feed efficiency and reduction of environmental footprint in broilers: a review

Currently, optimization of feed efficiency is one of the main challenges in improvement programs of livestock and poultry genetics. The objective of this review is to present the genetic aspects of feed efficiency related traits in meat-type chicken and possible ways to reduce the environmental impact of poultry meat production with effective breeding. Basic measures of feed efficiency are defined and the genetic background of these traits, including a review of heritabilities is described. Moreover, a number of genomic regions and candidate genes determining feed efficiency traits of broilers that were detected over the past decades are described. Classical and genomic selection strategies for feed efficiency in the context of its relationships with other performance traits are discussed as well. Finally, future strategies to improve feed digestibility are described as it is expected that they will decrease wastes and greenhouse gas emission. Further genetic improvement of feed efficiency, should be examined jointly with appropriate feeding strategies in broilers.

Divergently expressed gene identification and interaction prediction of long noncoding RNA and mRNA involved in duck reproduction

Long noncoding RNAs (lncRNAs) and divergently expressed genes exist widely in different tissues of mammals and birds, in which they are involved in various biological processes. However, there is limited information on their role in the regulation of normal biological processes during differentiation, development, and reproduction in birds. In this study, whole transcriptome strand-specific RNA sequencing of the ovary from young ducks (60days), first-laying ducks (160days), and old ducks, i.e., ducks that stopped laying eggs (490days) was performed. The lncRNAs and mRNAs from these ducks were systematically analyzed and identified by duck genome sequencing in the three study groups. The transcriptome from the duck ovary comprised 15,011 protein-coding genes and 2905 lncRNAs; all the lncRNAs were identified as novel long noncoding transcripts. The comparison of transcriptome data from different study groups identified 2240 divergent transcription genes and 135 divergently expressed lncRNAs, which differed among the groups; most of them were significantly downregulated with age. Among the divergent genes, 38 genes were related to the reproductive process and 6 genes were upregulated. Further prediction analysis revealed that 52 lncRNAs were closely correlated with divergent reproductive mRNAs. More importantly, 6 remarkable lncRNAs were correlated significantly with the conversion of the ovary in different phases. Our results aid in the understanding of the divergent transcriptome of duck ovary in different phases and the underlying mechanisms that drive the specificity of protein-coding genes and lncRNAs in duck ovary.

Genome-wide association study for performance traits in chickens using genotype by sequencing approach

Performance traits are economically important and are targets for selection in breeding programs, especially in the poultry industry. To identify regions on the chicken genome associated with performance traits, different genomic approaches have been applied in the last years. The aim of this study was the application of CornellGBS approach (134,528 SNPs generated from a PstI restriction enzyme) on Genome-Wide Association Studies (GWAS) in an outbred F₂ chicken population. We have validated 91.7% of these 134,528 SNPs after imputation of missed genotypes. Out of those, 20 SNPs were associated with feed conversion, one was associated with body weight at 35 days of age (P < 7.86E-07) and 93 were suggestively associated with a variety of performance traits (P < 1.57E-05). The majority of these SNPs (86.2%) overlapped with previously mapped QTL for the same performance traits and some of the SNPs also showed novel potential QTL regions. The results obtained in this study suggests future searches for candidate genes and QTL refinements as well as potential use of the SNPs described here in breeding programs.

Genomic evaluation by including dominance effects and inbreeding depression for purebred and crossbred performance with an application in pigs

BACKGROUND

Improved performance of crossbred animals is partly due to heterosis. One of the major genetic bases of heterosis is dominance, but it is seldom used in pedigree-based genetic evaluation of livestock. Recently, a trivariate genomic best linear unbiased prediction (GBLUP) model including dominance was developed, which can distinguish purebreds from crossbred animals explicitly. The objectives of this study were: (1) methodological, to show that inclusion of marker-based inbreeding accounts for directional dominance and inbreeding depression in purebred and crossbred animals, to revisit variance components of additive and dominance genetic effects using this model, and to develop marker-based estimators of genetic correlations between purebred and crossbred animals and of correlations of allele substitution effects between breeds; (2) to evaluate the impact of accounting for dominance effects and inbreeding depression on predictive ability for total number of piglets born (TNB) in a pig dataset composed of two purebred populations and their crossbreds. We also developed an equivalent model that makes the estimation of variance components tractable.

RESULTS

For TNB in Danish Landrace and Yorkshire populations and their reciprocal crosses, the estimated proportions of dominance genetic variance to additive genetic variance ranged from 5 to 11%. Genetic correlations between breeding values for purebred and crossbred performances for TNB ranged from 0.79 to 0.95 for Landrace and from 0.43 to 0.54 for Yorkshire across models. The estimated correlation of allele substitution effects between Landrace and Yorkshire was low for purebred performances, but high for crossbred performances. Predictive ability for crossbred animals was similar with or without dominance. The inbreeding depression effect increased predictive ability and the estimated inbreeding depression parameter was more negative for Landrace than for Yorkshire animals and was in between for crossbred animals.

CONCLUSIONS

Methodological developments led to closed-form estimators of inbreeding depression, variance components and correlations that can be easily interpreted in a quantitative genetics context. Our results confirm that genetic correlations of breeding values between purebred and crossbred performances within breed are positive and moderate. Inclusion of dominance in the GBLUP model does not improve predictive ability for crossbred animals, whereas inclusion of inbreeding depression does.

Imputation of genotypes in Danish purebred and two-way crossbred pigs using low-density panels

Background

Genotype imputation is commonly used as an initial step in genomic selection since the accuracy of genomic selection does not decline if accurately imputed genotypes are used instead of actual genotypes but for a lower cost. Performance of imputation has rarely been investigated in crossbred animals and, in particular, in pigs. The extent and pattern of linkage disequilibrium differ in crossbred versus purebred animals, which may impact the performance of imputation. In this study, first we compared different scenarios of imputation from 5 K to 8 K single nucleotide polymorphisms (SNPs) in genotyped Danish Landrace and Yorkshire and crossbred Landrace-Yorkshire datasets and, second, we compared imputation from 8 K to 60 K SNPs in genotyped purebred and simulated crossbred datasets. All imputations were done using software Beagle version 3.3.2. Then, we investigated the reasons that could explain the differences observed.

Results

Genotype imputation performs as well in crossbred animals as in purebred animals when both parental breeds are included in the reference population. When the size of the reference population is very large, it is not necessary to use a reference population that combines the two breeds to impute the genotypes of purebred animals because a within-breed reference population can provide a very high level of imputation accuracy (correct rate ≥ 0.99, correlation ≥ 0.95). However, to ensure that similar imputation accuracies are obtained for crossbred animals, a reference population that combines both parental purebred animals is required. Imputation accuracies are higher when a larger proportion of haplotypes are shared between the reference population and the validation (imputed) populations.

Conclusions

The results from both real data and pedigree-based simulated data demonstrate that genotype imputation from low-density panels to medium-density panels is highly accurate in both purebred and crossbred pigs. In crossbred pigs, combining the parental purebred animals in the reference population is necessary to obtain high imputation accuracy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0134-4) contains supplementary material, which is available to authorized users.

Using multiple regression, Bayesian networks and artificial neural networks for prediction of total egg production in European quails based on earlier expressed phenotypes

The prediction of total egg production (TEP) potential in poultry is an important task to aid optimized management decisions in commercial enterprises. The objective of the present study was to compare different modeling approaches for prediction of TEP in meat type quails (Coturnix coturnix coturnix) using phenotypes such as weight, weight gain, egg production and egg quality measurements. Phenotypic data on 30 traits from two lines (L1, n=180; and L2, n=205) of quail were modeled to predict TEP. Prediction models included multiple linear regression and artificial neural network (ANN). Moreover, Bayesian network (BN) and a stepwise approach were used as variable selection methods. BN results showed that TEP is independent from other earlier expressed traits when conditioned on egg production from 35 to 80 days of age (EP1). In addition, the prediction accuracy was much lower when EP1 was not included in the model. The best predictive model was ANN, after feature selection, showing prediction correlations of r=0.792 and r=0.714 for L1 and L2, respectively. In conclusion, machine learning methods may be useful, but reasonable prediction accuracies are obtained only when partial egg production measurements are included in the model.

Evaluation of measures of correctness of genotype imputation in the context of genomic prediction: a review of livestock applications

In livestock, many studies have reported the results of imputation to 50k single nucleotide polymorphism (SNP) genotypes for animals that are genotyped with low-density SNP panels. The objective of this paper is to review different measures of correctness of imputation, and to evaluate their utility depending on the purpose of the imputed genotypes. Across studies, imputation accuracy, computed as the correlation between true and imputed genotypes, and imputation error rates, that counts the number of incorrectly imputed alleles, are commonly used measures of imputation correctness. Based on the nature of both measures and results reported in the literature, imputation accuracy appears to be a more useful measure of the correctness of imputation than imputation error rates, because imputation accuracy does not depend on minor allele frequency (MAF), whereas imputation error rate depends on MAF. Therefore imputation accuracy can be better compared across loci with different MAF. Imputation accuracy depends on the ability of identifying the correct haplotype of a SNP, but many other factors have been identified as well, including the number of genotyped immediate ancestors, the number of animals with genotypes at the high-density panel, the SNP density on the low- and high-density panel, the MAF of the imputed SNP and whether imputed SNP are located at the end of a chromosome or not. Some of these factors directly contribute to the linkage disequilibrium between imputed SNP and SNP on the low-density panel. When imputation accuracy is assessed as a predictor for the accuracy of subsequent genomic prediction, we recommend that: (1) individual-specific imputation accuracies should be used that are computed after centring and scaling both true and imputed genotypes; and (2) imputation of gene dosage is preferred over imputation of the most likely genotype, as this increases accuracy and reduces bias of the imputed genotypes and the subsequent genomic predictions.

Genomic evaluation of both purebred and crossbred performances

Background

For a two-breed crossbreeding system, Wei and van der Werf presented a model for genetic evaluation using information from both purebred and crossbred animals. The model provides breeding values for both purebred and crossbred performances. Genomic evaluation incorporates marker genotypes into a genetic evaluation system. Among popular methods are the so-called single-step methods, in which marker genotypes are incorporated into a traditional animal model by using a combined relationship matrix that extends the marker-based relationship matrix to non-genotyped animals. However, a single-step method for genomic evaluation of both purebred and crossbred performances has not been developed yet.

Results

An extension of the Wei and van der Werf model that incorporates genomic information is presented. The extension consists of four steps: (1) the Wei van der Werf model is reformulated using two partial relationship matrices for the two breeds; (2) marker-based partial relationship matrices are constructed; (3) marker-based partial relationship matrices are adjusted to be compatible to pedigree-based partial relationship matrices and (4) combined partial relationship matrices are constructed using information from both pedigree and marker genotypes. The extension of the Wei van der Werf model can be implemented using software that allows inverse covariance matrices in sparse format as input.

Conclusions

A method for genomic evaluation of both purebred and crossbred performances was developed for a two-breed crossbreeding system. The method allows information from crossbred animals to be incorporated in a coherent manner for such crossbreeding systems.

Genomic selection: Status in different species and challenges for breeding

Technical advances and development in the market for genomic tools have facilitated access to whole-genome data across species. Building-up on the acquired knowledge of the genome sequences, large-scale genotyping has been optimized for broad use, so genotype information can be routinely used to predict genetic merit. Genomic selection (GS) refers to the use of aggregates of estimated marker effects as predictors which allow improved individual differentiation at young age. Realizable benefits of GS are influenced by several factors and vary in quantity and quality between species. General characteristics and challenges of GS in implementation and routine application are described, followed by an overview over the current status of its use, prospects and challenges in important animal species. Genetic gain for a particular trait can be enhanced by shortening of the generation interval, increased selection accuracy and increased selection intensity, with species- and breed-specific relevance of the determinants. Reliable predictions based on genetic marker effects require assembly of a reference for linking of phenotype and genotype data to allow estimation and regular re-estimation. Experiences from dairy breeding have shown that international collaboration can set the course for fast and successful implementation of innovative selection tools, so genomics may significantly impact the structures of future breeding and breeding programmes. Traits of great and increasing importance, which were difficult to improve in the conventional systems, could be emphasized, if continuous availability of high-quality phenotype data can be assured. Equally elaborate strategies for genotyping and phenotyping will allow tailored approaches to balance efficient animal production, sustainability, animal health and welfare in future.

The value of resequence data for poultry breeding: a primary layer breeder perspective

Poultry breeding companies are facing a new paradigm. Since 2004, extensive resources have been developed to increase understanding of the fundamental biology of the chicken. The chicken genome has been sequenced and revised twice, millions of novel DNA variants have been identified, and new tools have been created that allow rapid and inexpensive detection of these DNA variations. These developments have led to the establishment of molecular-based breeding programs within major poultry breeding companies that are revolutionizing the primary poultry breeding industries. Costs of sequencing continue to drop and are predicted to eventually reach the point where it is feasible to sequence the entire genome of elite birds before selection. There are multiple challenges to be resolved before this information can be fully incorporated into a breeding program. These include handling and analyzing the extremely large data sets generated, understanding which genes, variants, or both are relevant for commercial production traits, development of new bio-informatic tools, and integration of molecular information with traditional breeding programs. The novel variation identified within elite commercial lines will lead to enhancements in commercial breeding programs. Applications of this information include whole genomic selection, parentage identification, trait association studies, and quality control.

Genome-enabled prediction of quantitative traits in chickens using genomic annotation

BACKGROUND

Genome-wide association studies have been deemed successful for identifying statistically associated genetic variants of large effects on complex traits. Past studies have found enrichment of trait-associated SNPs in functionally annotated regions, while depletion was reported for intergenic regions (IGR). However, no systematic examination of connections between genomic regions and predictive ability of complex phenotypes has been carried out.

RESULTS

In this study, we partitioned SNPs based on their annotation to characterize genomic regions that deliver low and high predictive power for three broiler traits in chickens using a whole-genome approach. Additive genomic relationship kernels were constructed for each of the genic regions considered, and a kernel-based Bayesian ridge regression was employed as prediction machine. We found that the predictive performance for ultrasound area of breast meat from using genic regions marked by SNPs was consistently better than that from SNPs in IGR, while IGR tagged by SNPs were better than the genic regions for body weight and hen house egg production. We also noted that predictive ability delivered by the whole battery of markers was close to the best prediction achieved by one of the genomic regions.

CONCLUSIONS

Whole-genome regression methods use all available quality filtered SNPs into a model, contrary to accommodating only validated SNPs from exonic or coding regions. Our results suggest that, while differences among genomic regions in terms of predictive ability were observed, the whole-genome approach remains as a promising tool if interest is on prediction of complex traits.

Applied animal genomics: results from the field

Genomic selection (GS) is the use of statistical methods to estimate the genetic merit of a genotyped animal based on prediction equations derived from large ancestral populations with both phenotypes and genotypes. It has revolutionized the dairy cattle breeding industry and has been implemented with varying degrees of success in other animal breeding programs, including swine, poultry, and beef cattle. The findings of empirical field studies applying GS to the breeding sectors of these main animal protein industries are reviewed. Several translational considerations must be addressed before implementing GS in genetic improvement programs. These include determining and obtaining economically relevant phenotypes and determining the optimal size of the training population, cost-effective genotyping strategies, the practicality of field implementation, and the relative costs versus the benefits of the realized rate of genetic gain. GS may additionally change the optimal breeding scheme design, and studies that address this consideration are also reviewed briefly.

Epigenetics and phenotypic variability: some interesting insights from birds

Little is known about epigenetic mechanisms in birds with the exception of the phenomenon of dosage compensation of sex chromosomes, although such mechanisms could be involved in the phenotypic variability of birds, as in several livestock species. This paper reviews the literature on epigenetic mechanisms that could contribute significantly to trait variability in birds, and compares the results to the existing knowledge of epigenetic mechanisms in mammals. The main issues addressed in this paper are: (1) Does genomic imprinting exist in birds? (2) How does the embryonic environment influence the adult phenotype in avian species? (3) Does the embryonic environment have an impact on phenotypic variability across several successive generations? The potential for epigenetic studies to improve the performance of individual animals through the implementation of limited changes in breeding conditions or the addition of new parameters in selection models is still an open question.

Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo)

Background

The turkey (Meleagris gallopavo) is an important agricultural species and the second largest contributor to the world’s poultry meat production. Genetic improvement is attributed largely to selective breeding programs that rely on highly heritable phenotypic traits, such as body size and breast muscle development. Commercial breeding with small effective population sizes and epistasis can result in loss of genetic diversity, which in turn can lead to reduced individual fitness and reduced response to selection. The presence of genomic diversity in domestic livestock species therefore, is of great importance and a prerequisite for rapid and accurate genetic improvement of selected breeds in various environments, as well as to facilitate rapid adaptation to potential changes in breeding goals. Genomic selection requires a large number of genetic markers such as e.g. single nucleotide polymorphisms (SNPs) the most abundant source of genetic variation within the genome.

Results

Alignment of next generation sequencing data of 32 individual turkeys from different populations was used for the discovery of 5.49 million SNPs, which subsequently were used for the analysis of genetic diversity among the different populations. All of the commercial lines branched from a single node relative to the heritage varieties and the South Mexican turkey population. Heterozygosity of all individuals from the different turkey populations ranged from 0.17-2.73 SNPs/Kb, while heterozygosity of populations ranged from 0.73-1.64 SNPs/Kb. The average frequency of heterozygous SNPs in individual turkeys was 1.07 SNPs/Kb. Five genomic regions with very low nucleotide variation were identified in domestic turkeys that showed state of fixation towards alleles different than wild alleles.

Conclusion

The turkey genome is much less diverse with a relatively low frequency of heterozygous SNPs as compared to other livestock species like chicken and pig. The whole genome SNP discovery study in turkey resulted in the detection of 5.49 million putative SNPs compared to the reference genome. All commercial lines appear to share a common origin. Presence of different alleles/haplotypes in the SM population highlights that specific haplotypes have been selected in the modern domesticated turkey.