Estimation of genetic parameters for feed efficiency traits using random regression models in dairy cattle.

Feed efficiency has become an increasingly important research topic in recent years. As feed costs rise and the environmental impacts of agriculture become more apparent, improving the efficiency with which dairy cows convert feed to milk is increasingly important. However, feed intake is expensive to measure accurately on large populations, making the inclusion of this trait in breeding programs difficult. Understanding how the genetic parameters of feed efficiency and traits related to feed efficiency vary throughout the lactation period is valuable to gain understanding into the genetic nature of feed efficiency. This study used 121,226 dry matter intake (DMI) records, 120,500 energy-corrected milk (ECM) records, and 98,975 metabolic body weight (MBW) records, collected on 7,440 first-lactation Holstein cows from 6 countries (Canada, Denmark, Germany, Spain, Switzerland, and the United States), from January 2003 to February 2022. Genetic parameters were estimated using a multiple-trait random regression model with a fourth-order Legendre polynomial for all traits. Weekly phenotypes for DMI were re-parameterized using linear regressions of DMI on ECM and MBW, creating a measure of feed efficiency that was genetically corrected for ECM and MBW, referred to as genomic residual feed intake (gRFI). Heritability (SE) estimates varied from 0.15 (0.03) to 0.29 (0.02) for DMI, 0.24 (0.01) to 0.29 (0.03) for ECM, 0.55 (0.03) to 0.83 (0.05) for MBW, and 0.12 (0.03) to 0.22 (0.06) for gRFI. In general, heritability estimates were lower in the first stage of lactation compared with the later stages of lactation. Additive genetic correlations between weeks of lactation varied, with stronger correlations between weeks of lactation that were close together. The results of this study contribute to a better understanding of the change in genetic parameters across the first lactation, providing insight into potential selection strategies to include feed efficiency in breeding programs.

Genomic analysis of perinatal sucking reflex in German Brown Swiss calves.

A healthy sucking reflex is essential for newborn calves to ensure sufficient colostrum intake in the first few hours postpartum. In recent decades, European Brown Swiss breeders have repeatedly reported that some calves lack the ability to consume colostrum directly after birth due to an absent sucking reflex. In this study, we collected the phenotypes of more than 5,500 German Brown Swiss calves and performed variance component estimation with sire threshold models using Markov chain Monte Carlo algorithms. The 50K (777K) genotypes of nearly 2,000 (200) calves were collected, and an imputation was performed for all 50K genotypes up to 777K. Genome-wide association studies (GWAS) for the trait sucking reflex were conducted for all 777K genotypes. Depending on the trait coding, a low heritability was estimated to range from 0.08 to 0.11. The GWAS results identified 34 trait-associated SNP on 6 different chromosomes. Post-GWAS analyses showed significant overrepresentation of Gene Ontologies for central nervous development and several regulative processes. Functional annotation clustering and pathway analysis revealed relations to lipid metabolism, immune and endocrine systems, and signal transduction. The results of this study suggest that breeding for an improved sucking reflex is possible but requires large data sets for the estimation of reliable breeding values (either large progeny testing groups or a large reference genome in a genomic selection program).

Linkage disequilibrium pattern and genome-wide association mapping for meat traits in multiple porcine F2 crosses.

In the present study, data from four F2 crosses were analysed and used to study the linkage disequilibrium (LD) structure within and across the crosses. Genome-wide association analyses (GWASes) for conductivity and dressing out meat traits were conducted using single-marker and Bayesian multi-marker models using the pooled data from all F2 crosses. Porcine F2 crosses generated from the distantly related founder breeds Wild Boar, Piétrain and Meishan, as well as from a porcine F2 cross from the closely related founder breed Piétrain and an F1 Large White × Landrace cross were pooled. A total of 2572 F2 animals were genotyped using a 62K SNP chip. The positions of the SNPs were based on genome assembly Sscrofa11.1. After post-alignment and genotype filtering, approximately 50K SNPs were usable for LD studies and GWASes. The main findings of the present study are that the breakdown of LD was faster in crosses from closely related founder breeds compared to crosses from distantly related founders. The fastest breakdown of LD was observed by pooling the data. Based on the single-marker results and LD structure, clusters and windows were built for 1-Mb intervals. For conductivity and dressing out, 183 and 191 nominal significant associations respectively and six and five clusters respectively were found. Dominance was important for conductivity, and considering dominance in GWASes improved the mapping signals. Most clear signals were found for conductivity on SSC6, 8 and 15 and for dressing out on SSC2 and 7. Considering dominance might contribute to the accuracy of genomic selection and serve as a guide for choosing mating pairs with good combining abilities. However, further research is needed to investigate if dominance is also important in crossbreed pig breeding schemes.

Investigations on the pattern of linkage disequilibrium and selection signatures in the genomes of German Piétrain pigs.

The aim of this study was to study the population structure, to characterize the LD structure and to define core regions based on low recombination rates among SNP pairs in the genome of Piétrain pigs using data from the PorcineSNP60 BeadChip. This breed is a European sire line and was strongly selected for lean meat content during the last decades. The data were used to map signatures of selection using the REHH test. In the first step, selection signatures were searched genome-wide using only core haplotypes having a frequency above 0.25. In the second step, the results from the selection signature analysis were matched with the results from the recently conducted genome-wide association study for economical relevant traits to investigate putative overlaps of chromosomal regions. A small subdivision of the population with regard to the geographical origin of the individuals was observed. The extent of LD was determined genome-wide using r(2) values for SNP pairs with a distance ≤5 Mb and was on average 0.34. This comparable low r(2) value indicates a high genetic diversity in the Piétrain population. Six REHH values having a p-value < 0.001 were genome-wide detected. These were located on SSC1, 2, 6 and 17. Three positional candidate genes with potential biological roles were suggested, called LOC100626459, LOC100626014 and MIR1. The results imply that for genome-wide analysis especially in this population, a higher marker density and higher sample sizes are required. For a number of nine SNPs, which were successfully annotated to core regions, the REHH test was applied. However, no selection signatures were found for those regions (p-value < 0.1).

Genetic parameters for feather pecking and aggressive behavior in a large F2-cross of laying hens using generalized linear mixed models.

Feather pecking and aggressive pecking is a well-known problem in egg production. In the present study, genetic parameters for 4 feather-pecking-related traits were estimated using generalized linear mixed models. The traits were bouts of feather pecking delivered (FPD), bouts of feather pecking received (FPR), bouts of aggressive pecking delivered (APD), and bouts of aggressive pecking received (APR). An F2-design was established from 2 divergent selected founder lines. The lines were selected for low or high feather pecking for 10 generations. The number of F2 hens was 910. They were housed in pens with around 40 birds. Each pen was observed in 21 sessions of 20 min, distributed over 3 consecutive days. An animal model was applied that treated the bouts observed within 20 min as repeated observations. An over-dispersed Poisson distribution was assumed for observed counts and the link function was a log link. The model included a random animal effect, a random permanent environment effect, and a random day-by-hen effect. Residual variance was approximated on the link scale by the delta method. The results showed a heritability around 0.10 on the link scale for FPD and APD and of 0.04 for APR. The heritability of FPR was zero. For all behavior traits, substantial permanent environmental effects were observed. The approximate genetic correlation between FPD and APD (FPD and APR) was 0.81 (0.54). Egg production and feather eating records were collected on the same hens as well and were analyzed with a generalized linear mixed model, assuming a binomial distribution and using a probit link function. The heritability on the link scale for egg production was 0.40 and for feather eating 0.57. The approximate genetic correlation between FPD and egg production was 0.50 and between FPD and feather eating 0.73. Selection might help to reduce feather pecking, but this might result in an unfavorable correlated selection response reducing egg production. Feather eating and feather pecking are genetically correlated and this needs further investigation.

Genome-wide association analysis for growth, muscularity and meat quality in Piétrain pigs.

Improvement in growth and meat quality is one of the main objectives in sire line pig breeding programmes. Mapping quantitative trait loci for these traits using experimental crosses and a linkage-based approach has been performed frequently in the past. The Piétrain breed often was involved as a founder breed to establish the experimental crosses. This breed was selected for muscularity and leanness but shows relatively poor meat quality. It is frequently used as a sire line breed. With the advent of genome-wide and dense SNP chips in pig genomic research, it is possible to also conduct genome-wide association studies within the Piétrain breed. In this study, around 500 progeny-tested sires were genotyped with 60k SNPs. Data filtering showed that around 48k SNPs were useable in this sample. These SNPs were used to conduct a genome-wide association study for growth, muscularity and meat quality traits. Because it is known that a mutation in the RYR1 gene located on chromosome 6 shows a major effect on meat quality, this mutation was included in the models. Single-marker and multimarker association analyses were performed. The results revealed between zero and eight significant associations per trait with P < 5 × 10(-5) . Of special interest are SNPs located on SSC6, SSC10 and SSC15.

A two-step approach to map quantitative trait loci for meat quality in connected porcine F(2) crosses considering main and epistatic effects.

The aim of this study was to map QTL for meat quality traits in three connected porcine F(2) crosses comprising around 1000 individuals. The three crosses were derived from the founder breeds Chinese Meishan, European Wild Boar and Pietrain. The animals were genotyped genomewide for approximately 250 genetic markers, mostly microsatellites. They were phenotyped for seven meat quality traits (pH at 45 min and 24 h after slaughter, conductivity at 45 min and 24 h after slaughter, meat colour, drip loss and rigour). QTL mapping was conducted using a two-step procedure. In the first step, the QTL were mapped using a multi-QTL multi-allele model that was tailored to analyse multiple connected F(2) crosses. It considered additive, dominance and imprinting effects. The major gene RYR1:g.1843C>T affecting the meat quality on SSC6 was included as a cofactor in the model. The mapped QTL were tested for pairwise epistatic effects in the second step. All possible epistatic effects between additive, dominant and imprinting effects were considered, leading to nine orthogonal forms of epistasis. Numerous QTL were found. The most interesting chromosome was SSC6. Not all genetic variance of meat quality was explained by RYR1:g.1843C>T. A small confidence interval was obtained, which facilitated the identification of candidate genes underlying the QTL. Epistasis was significant for the pairwise QTL on SSC12 and SSC14 for pH24 and for the QTL on SSC2 and SSC5 for rigour. Some evidence for additional pairwise epistatic effects was found, although not significant. Imprinting was involved in epistasis.

Mapping quantitative trait loci for metabolic and cytological fatness traits of connected F2 crosses in pigs.

In the present study 3 connected F(2) crosses were used to map QTL for classical fat traits as well as fat-related metabolic and cytological traits in pigs. The founder breeds were Chinese Meishan, European Wild Boar, and Pietrain with to some extent the same founder animals in the different crosses. The different selection history of the breeds for fatness traits as well as the connectedness of the crosses led to a high statistical power. The total number of F(2) animals varied between 694 and 966, depending on the trait. The animals were genotyped for around 250 genetic markers, mostly microsatellites. The statistical model was a multi-allele, multi-QTL model that accounted for imprinting. The model was previously introduced from plant breeding experiments. The traits investigated were backfat depth and fat area as well as relative number of fat cells with different sizes and 2 metabolic traits (i.e., soluble protein content as an indicator for the level of metabolic turnover and NADP-malate dehydrogenase as an indicator for enzyme activity). The results revealed in total 37 significant QTL on chromosomes 1, 2, 4, 5, 6, 7, 8, 9, 14, 17, and 18, with often an overlap of confidence intervals of several traits. These confidence intervals were in some cases remarkably small, which is due to the high statistical power of the design. In total, 18 QTL showed significant imprinting effects. The small and overlapping confidence intervals for the classical fatness traits as well as for the cytological and metabolic traits enabled positional and functional candidate gene identification for several mapped QTL.