Genomic selection for crossbred performance accounting for breed-specific effects.

BACKGROUND: Breed-specific effects are observed when the same allele of a given genetic marker has a different effect depending on its breed origin, which results in different allele substitution effects across breeds. In such a case, single-breed breeding values may not be the most accurate predictors of crossbred performance. Our aim was to estimate the contribution of alleles from each parental breed to the genetic variance of traits that are measured in crossbred offspring, and to compare the prediction accuracies of estimated direct genomic values (DGV) from a traditional genomic selection model (GS) that are trained on purebred or crossbred data, with accuracies of DGV from a model that accounts for breed-specific effects (BS), trained on purebred or crossbred data. The final dataset was composed of 924 Large White, 924 Landrace and 924 two-way cross (F1) genotyped and phenotyped animals. The traits evaluated were litter size (LS) and gestation length (GL) in pigs. RESULTS: The genetic correlation between purebred and crossbred performance was higher than 0.88 for both LS and GL. For both traits, the additive genetic variance was larger for alleles inherited from the Large White breed compared to alleles inherited from the Landrace breed (0.74 and 0.56 for LS, and 0.42 and 0.40 for GL, respectively). The highest prediction accuracies of crossbred performance were obtained when training was done on crossbred data. For LS, prediction accuracies were the same for GS and BS DGV (0.23), while for GL, prediction accuracy for BS DGV was similar to the accuracy of GS DGV (0.53 and 0.52, respectively). CONCLUSIONS: In this study, training on crossbred data resulted in higher prediction accuracy than training on purebred data and evidence of breed-specific effects for LS and GL was demonstrated. However, when training was done on crossbred data, both GS and BS models resulted in similar prediction accuracies. In future studies, traits with a lower genetic correlation between purebred and crossbred performance should be included to further assess the value of the BS model in genomic predictions.

Feasibility of pedigree recording and genetic selection in village sheep flocks of smallholder farmers.

Pedigree recording and genetic selection in village flocks of smallholder farmers have been deemed infeasible by researchers and development workers. This is mainly due to the difficulty of sire identification under uncontrolled village breeding practices. A cooperative village sheep-breeding scheme was designed to achieve controlled breeding and implemented for Menz sheep of Ethiopia in 2009. In this paper, we evaluated the reliability of pedigree recording in village flocks by comparing genetic parameters estimated from data sets collected in the cooperative village and in a nucleus flock maintained under controlled breeding. Effectiveness of selection in the cooperative village was evaluated based on trends in breeding values over generations. Heritability estimates for 6-month weight recorded in the village and the nucleus flock were very similar. There was an increasing trend over generations in average estimated breeding values for 6-month weight in the village flocks. These results have a number of implications: the pedigree recorded in the village flocks was reliable; genetic parameters, which have so far been estimated based on nucleus data sets, can be estimated based on village recording; and appreciable genetic improvement could be achieved in village sheep selection programs under low-input smallholder farming systems.

Novel insight into the genomic architecture of feed and nitrogen efficiency measured by residual energy intake and nitrogen excretion in growing pigs.

BACKGROUND: Improvement of feed efficiency in pigs is of great economical and environmental interest and contributes to use limited resources efficiently to feed the world population. Genome scans for feed efficiency traits are of importance to reveal the underlying biological causes and increase the rate of genetic gain. The aim of this study was to determine the genomic architecture of feed efficiency measured by residual energy intake (REI), in association with production, feed conversion ratio (FCR) and nitrogen excretion traits through the identification of quantitative trait loci (QTL) at different stages of growth using a three generation full-sib design population which originated from a cross between Pietrain and a commercial dam line. RESULTS: Six novel QTL for REI were detected explaining 2.7-6.1% of the phenotypic variance in REI. At growth from 60-90 kg body weight (BW), a QTL with a significant dominance effect was identified for REI on SSC14, at a similar location to the QTL for feed intake and nitrogen excretion traits. At growth from 90-120 kg BW, three QTL for REI were detected on SSC2, SSC4 and SSC7 with significant additive, imprinting and additive effects, respectively. These QTL (except for the imprinted QTL) were positionally overlapping with QTL for FCR and nitrogen excretion traits. During final growth (120-140 kg BW), a further QTL for REI was identified on SSC8 with significant additive effect, which overlapped with QTL for nitrogen excretion. During entire analysed growth (60-140 kg BW), a novel additive QTL for REI on SSC4 was observed, with no overlapping with QTL for any other traits considered. CONCLUSIONS: The occurrence of only one overlapping QTL of REI with feed intake suggests that only a small proportion of the variance in REI was explained by change in feed intake, whereas four overlapping QTL of REI with those of nitrogen excretion traits suggests that mostly underlying factors of feed utilisation such as metabolism and protein turnover were the reason for change in REI. Different QTL for REI were identified at different growth stages, indicating that different genes are responsible for efficiency in feed utilisation at different stages of growth.

Genomic regions associated with bovine milk fatty acids in both summer and winter milk samples.

BACKGROUND: In this study we perform a genome-wide association study (GWAS) for bovine milk fatty acids from summer milk samples. This study replicates a previous study where we performed a GWAS for bovine milk fatty acids based on winter milk samples from the same population. Fatty acids from summer and winter milk are genetically similar traits and we therefore compare the regions detected in summer milk to the regions previously detected in winter milk GWAS to discover regions that explain genetic variation in both summer and winter milk. RESULTS: The GWAS of summer milk samples resulted in 51 regions associated with one or more milk fatty acids. Results are in agreement with most associations that were previously detected in a GWAS of fatty acids from winter milk samples, including eight 'new' regions that were not considered in the individual studies. The high correlation between the -log10(P-values) and effects of SNPs that were found significant in both GWAS imply that the effects of the SNPs were similar on winter and summer milk fatty acids. CONCLUSIONS: The GWAS of fatty acids based on summer milk samples was in agreement with most of the associations detected in the GWAS of fatty acids based on winter milk samples. Associations that were in agreement between both GWAS are more likely to be involved in fatty acid synthesis compared to regions detected in only one GWAS and are therefore worthwhile to pursue in fine-mapping studies.

Response to genomic selection: the Bulmer effect and the potential of genomic selection when the number of phenotypic records is limiting.

BACKGROUND: Over the last ten years, genomic selection has developed enormously. Simulations and results on real data suggest that breeding values can be predicted with high accuracy using genetic markers alone. However, to reach high accuracies, large reference populations are needed. In many livestock populations or even species, such populations cannot be established when traits are difficult or expensive to record, or when the population size is small. The value of genomic selection is then questionable. METHODS: In this study, we compare traditional breeding schemes based on own performance or progeny information to genomic selection schemes, for which the number of phenotypic records is limiting. Deterministic simulations were performed using selection index theory. Our focus was on the equilibrium response obtained after a few generations of selection. Therefore, we first investigated the magnitude of the Bulmer effect with genomic selection. RESULTS: Results showed that the reduction in response due to the Bulmer effect is the same for genomic selection as for selection based on traditional BLUP estimated breeding values, and is independent of the accuracy of selection. The reduction in response with genomic selection is greater than with selection based directly on phenotypes without the use of pedigree information, such as mass selection. To maximize the accuracy of genomic estimated breeding values when the number of phenotypic records is limiting, the same individuals should be phenotyped and genotyped, rather than genotyping parents and phenotyping their progeny. When the generation interval cannot be reduced with genomic selection, large reference populations are required to obtain a similar response to that with selection based on BLUP estimated breeding values based on own performance or progeny information. However, when a genomic selection scheme has a moderate decrease in generation interval, relatively small reference population sizes are needed to obtain a similar response to that with selection on traditional BLUP estimated breeding values. CONCLUSIONS: When the trait of interest cannot be recorded on the selection candidate, genomic selection schemes are very attractive even when the number of phenotypic records is limited, because traditional breeding requires progeny testing schemes with long generation intervals in those cases.

Genome-wide association study of insect bite hypersensitivity in two horse populations in the Netherlands.

BACKGROUND: Insect bite hypersensitivity is a common allergic disease in horse populations worldwide. Insect bite hypersensitivity is affected by both environmental and genetic factors. However, little is known about genes contributing to the genetic variance associated with insect bite hypersensitivity. Therefore, the aim of our study was to identify and quantify genomic associations with insect bite hypersensitivity in Shetland pony mares and Icelandic horses in the Netherlands. METHODS: Data on 200 Shetland pony mares and 146 Icelandic horses were collected according to a matched case-control design. Cases and controls were matched on various factors (e.g. region, sire) to minimize effects of population stratification. Breed-specific genome-wide association studies were performed using 70 k single nucleotide polymorphisms genotypes. Bayesian variable selection method Bayes-C with a threshold model implemented in GenSel software was applied. A 1 Mb non-overlapping window approach that accumulated contributions of adjacent single nucleotide polymorphisms was used to identify associated genomic regions. RESULTS: The percentage of variance explained by all single nucleotide polymorphisms was 13% in Shetland pony mares and 28% in Icelandic horses. The 20 non-overlapping windows explaining the largest percentages of genetic variance were found on nine chromosomes in Shetland pony mares and on 14 chromosomes in Icelandic horses. Overlap in identified associated genomic regions between breeds would suggest interesting candidate regions to follow-up on. Such regions common to both breeds (within 15 Mb) were found on chromosomes 3, 7, 11, 20 and 23. Positional candidate genes within 2 Mb from the associated windows were identified on chromosome 20 in both breeds. Candidate genes are within the equine lymphocyte antigen class II region, which evokes an immune response by recognizing many foreign molecules. CONCLUSIONS: The genome-wide association study identified several genomic regions associated with insect bite hypersensitivity in Shetland pony mares and Icelandic horses. On chromosome 20, associated genomic regions in both breeds were within 2 Mb from the equine lymphocyte antigen class II region. Increased knowledge on insect bite hypersensitivity associated genes will contribute to our understanding of its biology, enabling more efficient selection, therapy and prevention to decrease insect bite hypersensitivity prevalence.

Long-term response to genomic selection: effects of estimation method and reference population structure for different genetic architectures.

BACKGROUND: Genomic selection has become an important tool in the genetic improvement of animals and plants. The objective of this study was to investigate the impacts of breeding value estimation method, reference population structure, and trait genetic architecture, on long-term response to genomic selection without updating marker effects. METHODS: Three methods were used to estimate genomic breeding values: a BLUP method with relationships estimated from genome-wide markers (GBLUP), a Bayesian method, and a partial least squares regression method (PLSR). A shallow (individuals from one generation) or deep reference population (individuals from five generations) was used with each method. The effects of the different selection approaches were compared under four different genetic architectures for the trait under selection. Selection was based on one of the three genomic breeding values, on pedigree BLUP breeding values, or performed at random. Selection continued for ten generations. RESULTS: Differences in long-term selection response were small. For a genetic architecture with a very small number of three to four quantitative trait loci (QTL), the Bayesian method achieved a response that was 0.05 to 0.1 genetic standard deviation higher than other methods in generation 10. For genetic architectures with approximately 30 to 300 QTL, PLSR (shallow reference) or GBLUP (deep reference) had an average advantage of 0.2 genetic standard deviation over the Bayesian method in generation 10. GBLUP resulted in 0.6% and 0.9% less inbreeding than PLSR and BM and on average a one third smaller reduction of genetic variance. Responses in early generations were greater with the shallow reference population while long-term response was not affected by reference population structure. CONCLUSIONS: The ranking of estimation methods was different with than without selection. Under selection, applying GBLUP led to lower inbreeding and a smaller reduction of genetic variance while a similar response to selection was achieved. The reference population structure had a limited effect on long-term accuracy and response. Use of a shallow reference population, most closely related to the selection candidates, gave early benefits while in later generations, when marker effects were not updated, the estimation of marker effects based on a deeper reference population did not pay off.

Genome-wide association of milk fatty acids in Dutch dairy cattle.

BACKGROUND: Identifying genomic regions, and preferably individual genes, responsible for genetic variation in milk fat composition of bovine milk will enhance the understanding of biological pathways involved in fatty acid synthesis and may point to opportunities for changing milk fat composition via selective breeding. An association study of 50,000 single nucleotide polymorphisms (SNPs) was performed for even-chain saturated fatty acids (C4:0-C18:0), even-chain monounsaturated fatty acids (C10:1-C18:1), and the polyunsaturated C18:2cis9,trans11 (CLA) to identify genomic regions associated with individual fatty acids in bovine milk. RESULTS: The two-step single SNP association analysis found a total of 54 regions on 29 chromosomes that were significantly associated with one or more fatty acids. Bos taurus autosomes (BTA) 14, 19, and 26 showed highly significant associations with seven to ten traits, explaining a relatively large percentage of the total additive genetic variation. Many additional regions were significantly associated with the fatty acids. Some of the regions harbor genes that are known to be involved in fat synthesis or were previously identified as underlying quantitative trait loci for fat yield or content, such as ABCG2 and PPARGC1A on BTA 6; ACSS2 on BTA 13; DGAT1 on BTA 14; ACLY, SREBF1, STAT5A, GH, and FASN on BTA 19; SCD1 on BTA26; and AGPAT6 on BTA 27. CONCLUSIONS: Medium chain and unsaturated fatty acids are strongly influenced by polymorphisms in DGAT1 and SCD1. Other regions also showed significant associations with the fatty acids studied. These additional regions explain a relatively small percentage of the total additive genetic variance, but they are relevant to the total genetic merit of an individual and in unraveling the genetic background of milk fat composition. Regions identified in this study can be fine mapped to find causal mutations. The results also create opportunities for changing milk fat composition through breeding by selecting individuals based on their genetic merit for milk fat composition.

Genetic and phenotypic parameter estimates for body weights and egg production in Horro chicken of Ethiopia.

A breeding program has been established in 2008 to improve productivity of Horro chicken, an indigenous population in the western highlands of Ethiopia. The pedigree descended from 26 sires and 260 dams. Body weights were measured every 2 weeks from hatch to 8 weeks then every 4 weeks for the next 8 weeks. Egg production was recorded to 44 weeks of age for one generation. Genetic parameters were estimated using animal model fitted with common environmental effects for growth traits and ignoring common environment for egg production traits. Direct heritabilities ranged from low (0.15 ± 0.08), for body weight at 6 weeks, to moderate (0.40 ± 0.23), for hatch weight. Heritabilities of common environmental effects on growth were high at hatch (0.39 ± 0.10) and remained low afterwards. Age at first egg showed a very low heritability (0.06 ± 0.15). Heritabilities of egg numbers in the first, second, third, and fourth months of laying were 0.32 (±0.13), 0.20 (±0.16), 0.56 (±0.15), and 0.25 (±0.14), respectively. Heritabilities of cumulative of monthly records of egg numbers were from 0.24 ± 0.16 (for the first 2 months, EP12) to 0.35 ± 0.16 (over the 6 months, EP16). Body weight at 16 weeks of age (BW16) has a strong genetic correlation with the cumulative of monthly records: 0.92 (with EP12), 0.69 (with EP36), and 0.73 (with EP16). Besides their strong association, BW16 and EP16 showed higher heritability, relative to their respective trait categories. These two traits seemed to have common genes and utilizing them as selection traits would be expected to improve both egg production and growth performance of local chicken. However, the standard errors of estimates in this study were mostly high indicating that the estimates have low precision. Parameter estimations based on more data are needed before applying the current results in breeding programs.

Sensitivity of methods for estimating breeding values using genetic markers to the number of QTL and distribution of QTL variance.

The objective of this simulation study was to compare the effect of the number of QTL and distribution of QTL variance on the accuracy of breeding values estimated with genomewide markers (MEBV). Three distinct methods were used to calculate MEBV: a Bayesian Method (BM), Least Angle Regression (LARS) and Partial Least Square Regression (PLSR). The accuracy of MEBV calculated with BM and LARS decreased when the number of simulated QTL increased. The accuracy decreased more when QTL had different variance values than when all QTL had an equal variance. The accuracy of MEBV calculated with PLSR was affected neither by the number of QTL nor by the distribution of QTL variance. Additional simulations and analyses showed that these conclusions were not affected by the number of individuals in the training population, by the number of markers and by the heritability of the trait. Results of this study show that the effect of the number of QTL and distribution of QTL variance on the accuracy of MEBV depends on the method that is used to calculate MEBV.

Production objectives and trait preferences of village poultry producers of Ethiopia: implications for designing breeding schemes utilizing indigenous chicken genetic resources.

To generate information essential for the implementation of breeding schemes suitable for village poultry producers in Ethiopia, a survey was conducted aimed at defining the socioeconomic characteristics of the production environments in different geographic regions, understanding the important functions of chickens, identifying farmers' choice of chicken breeds and the underlying factors that determine the choice of genetic stock used. The survey included both questionnaire survey and a participatory group discussion. A total of 225 households (45 households from each of five Woredas) were interviewed. The questionnaire was designed to collect data covering general information on village poultry production such as socio-management characteristics, production objectives, population structure, breed choice and trait preferences, market preferences of specific traits, and farmers' selection practices. The participatory farmers' discussions were designed to involve stakeholders in defining the breeding objective "traits" and deriving their relative importance in the production environment based on the different functions of chickens and "traits" identified in the interviews. The results showed that production of eggs for consumption is the principal function of chickens in most regions followed by the use as source of income and meat for home consumption. The production system in all geographic regions studied revealed similar features generally characterized by extensive scavenging management, absence of immunization programs, increased risk of exposure of birds to disease and predators, and reproduction entirely based on uncontrolled natural mating and hatching of eggs using broody hens. Farmers' ratings of indigenous chickens with respect to modern breeds showed the highest significance of the adaptive traits in general, and the superior merits of indigenous chickens to high yielding exotic breeds in particular. Adaptation to the production environment was the most important attribute of chickens in all the study areas. The high significance attributed to reproduction traits indicates the need for maintaining broody behavior and high level of hatchability while breeding for improved productivity of indigenous chickens for village conditions. The market price of chickens is primarily dictated by weight, but farmers rated growth (males) and number of eggs followed by growth (females) as the production traits they would like the most to be improved. Therefore, the ultimate breeding goal should be to develop a dual-purpose breed based on indigenous chicken genetic resources with any of the comb types other than single for all the regions studied having the most preferred white body plumage for farmers in the Amhara region and red body plumage for those in Oromia, Benshangul-Gumuz, and Southern regions.

Multilevel selection 1: Quantitative genetics of inheritance and response to selection.

Interaction among individuals is universal, both in animals and in plants, and substantially affects evolution of natural populations and responses to artificial selection in agriculture. Although quantitative genetics has successfully been applied to many traits, it does not provide a general theory accounting for interaction among individuals and selection acting on multiple levels. Consequently, current quantitative genetic theory fails to explain why some traits do not respond to selection among individuals, but respond greatly to selection among groups. Understanding the full impacts of heritable interactions on the outcomes of selection requires a quantitative genetic framework including all levels of selection and relatedness. Here we present such a framework and provide expressions for the response to selection. Results show that interaction among individuals may create substantial heritable variation, which is hidden to classical analyses. Selection acting on higher levels of organization captures this hidden variation and therefore always yields positive response, whereas individual selection may yield response in the opposite direction. Our work provides testable predictions of response to multilevel selection and reduces to classical theory in the absence of interaction. Statistical methodology provided elsewhere enables empirical application of our work to both natural and domestic populations.

Multilevel selection 2: Estimating the genetic parameters determining inheritance and response to selection.

Interactions among individuals are universal, both in animals and in plants and in natural as well as domestic populations. Understanding the consequences of these interactions for the evolution of populations by either natural or artificial selection requires knowledge of the heritable components underlying them. Here we present statistical methodology to estimate the genetic parameters determining response to multilevel selection of traits affected by interactions among individuals in general populations. We apply these methods to obtain estimates of genetic parameters for survival days in a population of layer chickens with high mortality due to pecking behavior. We find that heritable variation is threefold greater than that obtained from classical analyses, meaning that two-thirds of the full heritable variation is hidden to classical analysis due to social interactions. As a consequence, predicted responses to multilevel selection applied to this population are threefold greater than classical predictions. This work, combined with the quantitative genetic theory for response to multilevel selection presented in an accompanying article in this issue, enables the design of selection programs to effectively reduce competitive interactions in livestock and plants and the prediction of the effects of social interactions on evolution in natural populations undergoing multilevel selection.

Estimating relatedness between individuals in general populations with a focus on their use in conservation programs.

Relatedness estimators are widely used in genetic studies, but effects of population structure on performance of estimators, criteria to evaluate estimators, and benefits of using such estimators in conservation programs have to date received little attention. In this article we present new estimators, based on the relationship between coancestry and molecular similarity between individuals, and compare them with existing estimators using Monte Carlo simulation of populations, either panmictic or structured. Estimators were evaluated using statistical criteria and a diversity criterion that minimized relatedness. Results show that ranking of estimators depends on the population structure. An existing estimator based on two-gene and four-gene coefficients of identity performs best in panmictic populations, whereas a new estimator based on coancestry performs best in structured populations. The number of marker alleles and loci did not affect ranking of estimators. Statistical criteria were insufficient to evaluate estimators for their use in conservation programs. The regression coefficient of pedigree relatedness on estimated relatedness (beta2) was substantially lower than unity for all estimators, causing overestimation of the diversity conserved. A simple correction to achieve beta2 = 1 improves both existing and new estimators. Using relatedness estimates with correction considerably increased diversity in structured populations, but did not do so or even decreased diversity in panmictic populations.

On the detection of imprinted quantitative trait loci in experimental crosses of outbred species.

In this article, the quantitative genetic aspects of imprinted genes and statistical properties of methods to detect imprinted QTL are studied. Different models to detect imprinted QTL and to distinguish between imprinted and Mendelian QTL were compared in a simulation study. Mendelian and imprinted QTL were simulated in an F2 design and analyzed under Mendelian and imprinting models. Mode of expression was evaluated against the H(0) of a Mendelian QTL as well as the H(0) of an imprinted QTL. It was shown that imprinted QTL might remain undetected when analyzing the genome with Mendelian models only. Compared to testing against a Mendelian QTL, using the H(0) of an imprinted QTL gave a higher proportion of correctly identified imprinted QTL, but also gave a higher proportion of false inference of imprinting for Mendelian QTL. When QTL were segregating in the founder lines, spurious detection of imprinting became more prominent under both tests, especially for designs with a small number of F1 sires.

Dairy cattle production in Europe.

The European Union (EU) is a major player on world markets for most dairy products and produces the largest single share of the global market. Dairying is one of the most profitable sectors of EU agriculture. Milk yields per cow have increased steadily in every member state between 1985 and 1997. Overall EU dairy production continues to follow a trend towards increased intensification on a smaller number of larger, more specialised production units. In this paper we highlight the main characteristics of dairy cattle production in the EU, as well as their impact on cow and herd reproductive performances. AI is widely used in areas where there is high input:high output farming and farmers' primary objective is to have each cow produce a calf every 12 months without any particular calving pattern. The use of embryo technologies is increasing and in some areas (e.g. IVP embryo transfers) Europe leads the world. In most areas of Europe, because of both year round calving and ethical or consumer concerns, hormones are mainly used to treat reproductive diseases and not for pharmaceutical control of breeding.

Factors affecting commercial application of embryo technologies in dairy cattle in Europe--a modelling approach.

Reproductive techniques have a major impact on the structure of breeding programmes, the rate of genetic gain and dissemination of genetic gain in populations. This manuscript reviews the impact of reproductive technologies on the underlying components of genetic gain and inbreeding with special reference to the role of female reproductive technology. Evaluation of alternative breeding schemes should be based on genetic gain while constraining inbreeding. Optimum breeding schemes can be characterised by: decreased importance of sib information; increased accuracy at the expense of intensity; and a factorial mating strategy. If large-scale embryo cloning becomes feasible, this will have a small impact on the rate of genetic gain but will have a large impact on the structure of breeding programmes.

The prospects of selection for social genetic effects to improve welfare and productivity in livestock.

Social interactions between individuals living in a group can have both positive and negative effects on welfare, productivity, and health of these individuals. Negative effects of social interactions in livestock are easier to observe than positive effects. For example, laying hens may develop feather pecking, which can cause mortality due to cannibalism, and pigs may develop tail biting or excessive aggression. Several studies have shown that social interactions affect the genetic variation in a trait. Genetic improvement of socially-affected traits, however, has proven to be difficult until relatively recently. The use of classical selection methods, like individual selection, may result in selection responses opposite to expected, because these methods neglect the effect of an individual on its group mates (social genetic effects). It has become clear that improvement of socially-affected traits requires selection methods that take into account not only the direct effect of an individual on its own phenotype but also the social genetic effects, also known as indirect genetic effects, of an individual on the phenotypes of its group mates. Here, we review the theoretical and empirical work on social genetic effects, with a focus on livestock. First, we present the theory of social genetic effects. Subsequently, we evaluate the evidence for social genetic effects in livestock and other species, by reviewing estimates of genetic parameters for direct and social genetic effects. Then we describe the results of different selection experiments. Finally, we discuss issues concerning the implementation of social genetic effects in livestock breeding programs. This review demonstrates that selection for socially-affected traits, using methods that target both the direct and social genetic effects, is a promising, but sometimes difficult to use in practice, tool to simultaneously improve production and welfare in livestock.

Risk factors for insect bite hypersensitivity in Friesian horses and Shetland ponies in The Netherlands.

Insect bite hypersensitivity (IBH) is an equine skin allergy caused by bites of Culicoides spp. and impacts on the welfare of affected horses. The aim of this study was to identify and quantify risk factors for IBH. Data from 3453 Friesian horse mares and 7074 Shetland pony mares scored for IBH by inspectors during obligatory foal inspections were analysed using breed-specific multivariable logistic regression models. The combined effect of month and year of scoring, Province and inspector were significantly associated with IBH in both breeds. In Shetland pony mares, withers height and coat colour were also significantly associated with IBH, while body condition had a nearly significant effect. The outcomes from this study on risk factors might contribute to the development of more efficient measures to reduce the prevalence of IBH.

Genetic variation in vitamin B-12 content of bovine milk and its association with SNP along the bovine genome.

Vitamin B-12 (also called cobalamin) is essential for human health and current intake levels of vitamin B-12 are considered to be too low. Natural enrichment of the vitamin B-12 content in milk, an important dietary source of vitamin B-12, may help to increase vitamin B-12 intake. Natural enrichment of the milk vitamin B-12 content could be achieved through genetic selection, provided there is genetic variation between cows with respect to the vitamin B-12 content in their milk. A substantial amount of genetic variation in vitamin B-12 content was detected among raw milk samples of 544 first-lactation Dutch Holstein Friesian cows. The presence of genetic variation between animals in vitamin B-12 content in milk indicates that the genotype of the cow affects the amount of vitamin B-12 that ends up in her milk and, consequently, that the average milk vitamin B-12 content of the cow population can be increased by genetic selection. A genome-wide association study revealed significant association between 68 SNP and vitamin B-12 content in raw milk of 487 first-lactation Dutch Holstein Friesian cows. This knowledge facilitates genetic selection for milk vitamin B-12 content. It also contributes to the understanding of the biological mechanism responsible for the observed genetic variation in vitamin B-12 content in milk. None of the 68 significantly associated SNP were in or near known candidate genes involved in transport of vitamin B-12 through the gastrointestinal tract, uptake by ileum epithelial cells, export from ileal cells, transport through the blood, uptake from the blood, intracellular processing, or reabsorption by the kidneys. Probably, associations relate to genes involved in alternative pathways of well-studied processes or to genes involved in less well-studied processes such as ruminal production of vitamin B-12 or secretion of vitamin B-12 by the mammary gland.