A first characterization of the microbiota-resilience link in swine.

BACKGROUND: The gut microbiome plays a crucial role in understanding complex biological mechanisms, including host resilience to stressors. Investigating the microbiota-resilience link in animals and plants holds relevance in addressing challenges like adaptation of agricultural species to a warming environment. This study aims to characterize the microbiota-resilience connection in swine. As resilience is not directly observable, we estimated it using four distinct indicators based on daily feed consumption variability, assuming animals with greater intake variation may face challenges in maintaining stable physiological status. These indicators were analyzed both as linear and categorical variables. In our first set of analyses, we explored the microbiota-resilience link using PERMANOVA, α-diversity analysis, and discriminant analysis. Additionally, we quantified the ratio of estimated microbiota variance to total phenotypic variance (microbiability). Finally, we conducted a Partial Least Squares-Discriminant Analysis (PLS-DA) to assess the classification performance of the microbiota with indicators expressed in classes. RESULTS: This study offers four key insights. Firstly, among all indicators, two effectively captured resilience. Secondly, our analyses revealed robust relationship between microbial composition and resilience in terms of both composition and richness. We found decreased α-diversity in less-resilient animals, while specific amplicon sequence variants (ASVs) and KEGG pathways associated with inflammatory responses were negatively linked to resilience. Thirdly, considering resilience indicators in classes, we observed significant differences in microbial composition primarily in animals with lower resilience. Lastly, our study indicates that gut microbial composition can serve as a reliable biomarker for distinguishing individuals with lower resilience. CONCLUSION: Our comprehensive analyses have highlighted the host-microbiota and resilience connection, contributing valuable insights to the existing scientific knowledge. The practical implications of PLS-DA and microbiability results are noteworthy. PLS-DA suggests that host-microbiota interactions could be utilized as biomarkers for monitoring resilience. Furthermore, the microbiability findings show that leveraging host-microbiota insights may improve the identification of resilient animals, supporting their adaptive capacity in response to changing environmental conditions. These practical implications offer promising avenues for enhancing animal well-being and adaptation strategies in the context of environmental challenges faced by livestock populations. Video Abstract.

Effect of extended heat stress in dairy cows on productive and behavioral traits.

This study evaluates the response of dairy cows to short and extended heat stressing conditions (from 1 to 28 days), as expressed in changes in their behavior. Due to climate change, heat stress and strong heat waves are frequently affecting the productivity and behavior of dairy cows. In the five years under study from 2018 to 2022, two were characterized by extremely strong heat waves occurring in the region analyzed in this study (Northern Italy). The dairy cattle farm involved in this study is located in Northern Italy and includes about 1 600 Holstein Friesian lactating dairy cows. Phenotypic data were provided by the Afimilk system and compromised behavioral and productive traits. Behavioral traits analyzed were activity, rest time, rest bouts, rest ratio, rest per bout and restlessness. Production traits were daily milk yield, average milking time, somatic cell count, fat percentage, protein percentage and lactose percentage. Climate data came from the National Aeronautics and Space Administration/Prediction of Worldwide Energy Resources database. Heat stress was analyzed considering Temperature-Humidity Index (THI) averaged over 28 different time windows of continuous heat stress. Results showed that rest time and milk yield were the two traits most affected by the increased THI. Rest time was immediately affected by high THI, showing a marked decrease already from 1d window and maintaining this all over the other windows. Furthermore, results show that rest time and rest ratio were only slightly negatively correlated with milk yield (-0.14 and -0.15). In addition, heat stress has a different effect depending on parity and lactation stages on the studied traits. In conclusion, the results indicate that heat stress increases activity and compromises milk production, rest time and milk quality traits. Results further suggest that rest time can be a better parameter than activity to describe the effects of heat stress on dairy cattle. The novel approach used in this study is based on the use of different time windows (up to 28 days) before the emergence of undesired THI and allows to identify the traits that are immediately influenced by the undesirable THI values and those that are influenced only after a prolonged heat stress period.

Persistence of autozygosity in crossbreds between autochthonous and cosmopolitan breeds of swine: a simulation study.

Crossbreeding might be a valid strategy to valorize local pig breeds. Crossbreeding should reduce homozygosity and, as a consequence, yield hybrid vigor for fitness and production traits. This study aimed to quantify the persistence of autozygosity in terminal crossbred pigs compared with purebreds and, in turn, identify genomic regions where autozygosity's persistence would not be found. The study was based on genotyping data from 20 European local pig breeds and three cosmopolitan pig breeds used to simulate crossbred offspring. This study consisted of two steps. First, one hundred matings were simulated for each pairwise combination of the 23 considered breeds (for a total of 276 combinations), ignoring the sex of the parent individuals in order to generate purebred and crossbred matings leveraging all the germplasm available. Second, a few preselected terminal-maternal breed pairs were used to mimic a realistic terminal crossbreeding system: (i) Mora Romagnola (boars) or Cinta Senese (boars) crossed with Large White (sows) or Landrace (sows); (ii) Duroc (boars) crossed with Mora Romagnola (sows) or Cinta Senese (sows). Runs of homozygosity was used to estimate genome-wide autozygosity (FROH). Observed FROH was higher in purebreds than in crossbreds, although some crossbred combinations showed higher FROH than other purebred combinations. Among the purebreds, the highest FROH values were observed in Mora Romagnola and Turopolje (0.50 and 0.46, respectively). FROH ranged from 0.04 to 0.16 in the crossbreds Alentejana × Large White and Alentejana × Iberian, respectively. Persistence of autozygosity was found in several genomic segments harboring regions where quantitative trait loci (QTLs) were found in the literature. The regions were enriched in QTLs involved in fatty acid metabolism and associated with performance traits. This simulation shows that autozygosity persists in most breed combinations of terminal crosses. Results suggest that a strategy for crossbreeding is implemented when leveraging autochthonous and cosmopolitan breeds to obtain most of the hybrid vigor.

The potential of microbiota information to better predict efficiency traits in growing pigs fed a conventional and a high-fiber diet.

BACKGROUND: Improving pigs' ability to digest diets with an increased dietary fiber content is a lever to improve feed efficiency and limit feed costs in pig production. The aim of this study was to determine whether information on the gut microbiota and host genetics can contribute to predict digestive efficiency (DE, i.e. digestibility coefficients of energy, organic matter, and nitrogen), feed efficiency (FE, i.e. feed conversion ratio and residual feed intake), average daily gain, and daily feed intake phenotypes. Data were available for 1082 pigs fed a conventional or high-fiber diet. Fecal samples were collected at 16 weeks, and DE was estimated using near­infrared spectrometry. A cross-validation approach was used to predict traits within the same diet, for the opposite diet, and for a combination of both diets, by implementing three models, i.e. with only genomic (Gen), only microbiota (Micro), and both genomic and microbiota information (Micro+Gen). The predictive ability with and without sharing common sires and breeding environment was also evaluated. Prediction accuracy of the phenotypes was calculated as the correlation between model prediction and phenotype adjusted for fixed effects. RESULTS: Prediction accuracies of the three models were low to moderate (< 0.47) for growth and FE traits and not significantly different between models. In contrast, for DE traits, prediction accuracies of model Gen were low (< 0.30) and those of models Micro and Micro+Gen were moderate to high (> 0.52). Prediction accuracies were not affected by the stratification of diets in the reference and validation sets and were in the same order of magnitude within the same diet, for the opposite diet, and for the combination of both diets. Prediction accuracies of the three models were significantly higher when pigs in the reference and validation populations shared common sires and breeding environment than when they did not (P < 0.001). CONCLUSIONS: The microbiota is a relevant source of information to predict DE regardless of the diet, but not to predict growth and FE traits for which prediction accuracies were similar to those obtained with genomic information only. Further analyses on larger datasets and more diverse diets should be carried out to complement and consolidate these results.

Imputation accuracy from low- to medium-density SNP chips for US crossbred dairy cattle.

This study aimed at evaluating the quality of imputation accuracy (IA) by marker (IAm) and by individual (IAi) in US crossbred dairy cattle. Holstein × Jersey crossbreds were used to evaluate IA from a low- (7K) to a medium-density (50K) SNP chip. Crossbred animals, as well as their sires (53), dams (77), and maternal grandsires (63), were all genotyped with a 78K SNP chip. Seven different scenarios of reference populations were tested, in which some scenarios used different family relationships and others added random unrelated purebred and crossbred individuals to those different family relationship scenarios. The same scenarios were tested on Holstein and Jersey purebred animals to compare these outcomes against those attained in crossbred animals. The genotype imputation was performed with findhap (version 4) software (VanRaden, 2015). There were no significant differences in IA results depending on whether the sire of imputed individuals was Holstein and the dam was Jersey, or vice versa. The IA increased significantly with the addition of related individuals in the reference population, from 86.70 ± 0.06% when only sires or dams were included in the reference population to 90.09 ± 0.06% when sire (S), dam (D), and maternal grandsire genomic data were combined in the reference population. In all scenarios including related individuals in the reference population, IAm and IAi were significantly superior in purebred Jersey and Holstein animals than in crossbreds, ranging from 90.75 ± 0.06 to 94.02 ± 0.06%, and from 90.88 ± 0.11 to 94.04 ± 0.10%, respectively. Additionally, a scenario called SPB+DLD(where PB indicates purebread and LD indicates low density), similar to the genomic evaluations performed on US crossbred dairy, was tested. In this scenario, the information from the 5 evaluated breeds (Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey) genotyped with a 50K SNP chip and genomic information from the dams genotyped with a 7K SNP chip were combined in the reference population, and the IAm and IAi were 80.87 ± 0.06% and 80.85 ± 0.08%, respectively. Adding randomly nonrelated genotyped individuals in the reference population reduced IA for both purebred and crossbred cows, except for scenario SPB+DLD, where adding crossbreds to the reference population increased IA values. Our findings demonstrate that IA for US Holstein × Jersey crossbred ranged from 85 to 90%, and emphasize the significance of designing and defining the reference population for improved IA.

Inbreeding depression for producer-recorded udder, metabolic, and reproductive diseases in US dairy cattle.

This study leveraged a growing dataset of producer-recorded phenotypes for mastitis, reproductive diseases (metritis and retained placenta), and metabolic diseases (ketosis, milk fever, and displaced abomasum) to investigate the potential presence of inbreeding depression for these disease traits. Phenotypic, pedigree, and genomic information were obtained for 354,043 and 68,292 US Holstein and Jersey cows, respectively. Total inbreeding coefficients were calculated using both pedigree and genomic information; the latter included inbreeding estimates obtained using a genomic relationship matrix and runs of homozygosity. We also generated inbreeding coefficients based on the generational inbreeding for recent and old pedigree inbreeding, for different run-of-homozygosity length classes, and for recent and old homozygous-by-descent segment-based inbreeding. Estimates on the liability scale revealed significant evidence of inbreeding depression for reproductive-disease traits, with an increase in total pedigree and genomic inbreeding showing a notable effect for recent inbreeding. However, we found inconsistent evidence for inbreeding depression for mastitis or any metabolic diseases. Notably, in Holsteins, the probability of developing displaced abomasum decreased with inbreeding, particularly for older inbreeding. Estimates of disease probability for cows with low, average, and high inbreeding levels did not significantly differ across any inbreeding coefficient and trait combination, indicating that although inbreeding may affect disease incidence, it likely plays a smaller role compared with management and environmental factors.

Genomic predictions and GWAS for heat tolerance in pigs based on reaction norm models with performance records and data from public weather stations considering alternative temperature thresholds.

Genetic improvement of livestock productivity has resulted in greater production of metabolic heat and potentially greater susceptibility to heat stress. Various studies have demonstrated that there is genetic variability for heat tolerance and genetic selection for more heat tolerant individuals is possible. The rate of genetic progress tends to be greater when genomic information is incorporated into the analyses as more accurate breeding values can be obtained for young individuals. Therefore, this study aimed (1) to evaluate the predictive ability of genomic breeding values for heat tolerance based on routinely recorded traits, and (2) to investigate the genetic background of heat tolerance based on single-step genome-wide association studies for economically important traits related to body composition, growth and reproduction in Large White pigs. Pedigree information was available for 265,943 animals and genotypes for 8686 animals. The studied traits included ultrasound backfat thickness (BFT), ultrasound muscle depth (MDP), piglet weaning weight (WW), off-test weight (OTW), interval between farrowing (IBF), total number of piglets born (TNB), number of piglets born alive (NBA), number of piglets born dead (NBD), number of piglets weaned (WN) and weaning-to-estrus interval (IWE). The number of phenotypic records ranged from 6059 (WN) to 172,984 (TNB). Single-step genomic reaction norm predictions were used to calculate the genomic estimated breeding values for each individual. Predictions of breeding values for the validation population individuals were compared between datasets containing phenotypic records measured in the whole range of temperatures (WR) and datasets containing only phenotypic records measured when the weather station temperature was above 10°C (10C) or 15°C (15C), to evaluate the usefulness of these datasets that may better reflect the within-barn temperature. The use of homogeneous or heterogeneous residual variance was found to be trait-dependent, where homogeneous variance presented the best fit for MDP, BFT, OTW, TNB, NBA, WN and IBF, while the other traits (WW and IWE) had better fit with heterogeneous variance. The average prediction accuracy, dispersion and bias values considering all traits for WR were 0.36 ± 0.05, -0.07 ± 0.13 and 0.76 ± 0.10, respectively; for 10C were 0.39 ± 0.02, -0.05 ± 0.07 and 0.81 ± 0.05, respectively; and for 15C were 0.32 ± 0.05, -0.05 ± 0.11 and 0.84 ± 0.10, respectively. Based on the studied traits, using phenotypic records collected when the outside temperature (from public weather stations) was above 10°C provided better predictions for most of the traits. Forty-three and 62 candidate genomic regions were associated with the intercept (overall performance level) and slope term (specific biological mechanisms related to environmental sensitivity), respectively. Our results contribute to improve genomic predictions using existing datasets and better understand the genetic background of heat tolerance in pigs. Furthermore, the genomic regions and candidate genes identified will contribute to future genomic studies and breeding applications.

Longitudinal genomic analyses of automatically-recorded vaginal temperature in lactating sows under heat stress conditions based on random regression models.

BACKGROUND: Automatic and continuous recording of vaginal temperature (TV) using wearable sensors causes minimal disruptions to animal behavior and can generate data that enable the evaluation of temporal body temperature variation under heat stress (HS) conditions. However, the genetic basis of TV in lactating sows from a longitudinal perspective is still unknown. The objectives of this study were to define statistical models and estimate genetic parameters for TV in lactating sows using random regression models, and identify genomic regions and candidate genes associated with HS indicators derived from automatically-recorded TV. RESULTS: Heritability estimates for TV ranged from 0.14 to 0.20 over time (throughout the day and measurement period) and from 0.09 to 0.18 along environmental gradients (EG, - 3.5 to 2.2, which correspond to dew point values from 14.87 to 28.19 ËšC). Repeatability estimates of TV over time and along EG ranged from 0.57 to 0.66 and from 0.54 to 0.77, respectively. TV measured from 12h00 to 16h00 had moderately high estimates of heritability (0.20) and repeatability (0.64), indicating that this period might be the most suitable for recording TV for genetic selection purposes. Significant genotype-by-environment interactions (GxE) were observed and the moderately high estimates of genetic correlations between pairs of extreme EG indicate potential re-ranking of selection candidates across EG. Two important genomic regions on chromosomes 10 (59.370-59.998 Mb) and16 (21.548-21.966 Mb) were identified. These regions harbor the genes CDC123, CAMK1d, SEC61A2, and NUDT5 that are associated with immunity, protein transport, and energy metabolism. Across the four time-periods, respectively 12, 13, 16, and 10 associated genomic regions across 14 chromosomes were identified for TV. For the three EG classes, respectively 18, 15, and 14 associated genomic windows were identified for TV, respectively. Each time-period and EG class had uniquely enriched genes with identified specific biological functions, including regulation of the nervous system, metabolism and hormone production. CONCLUSIONS: TV is a heritable trait with substantial additive genetic variation and represents a promising indicator trait to select pigs for improved heat tolerance. Moderate GxE for TV exist, indicating potential re-ranking of selection candidates across EG. TV is a highly polygenic trait regulated by a complex interplay of physiological, cellular and behavioral mechanisms.

Domestication shapes the pig gut microbiome and immune traits from the scale of lineage to population.

Animal ecology and evolution have long been known to shape host physiology, but more recently, the gut microbiome has been identified as a mediator between animal ecology and evolution and health. The gut microbiome has been shown to differ between wild and domestic animals, but the role of these differences for domestic animal evolution remains unknown. Gut microbiome responses to new animal genotypes and local environmental change during domestication may promote specific host phenotypes that are adaptive (or not) to the domestic environment. Because the gut microbiome supports host immune function, understanding the effects of animal ecology and evolution on the gut microbiome and immune phenotypes is critical. We investigated how domestication affects the gut microbiome and host immune state in multiple pig populations across five domestication contexts representing domestication status and current living conditions: free-ranging wild, captive wild, free-ranging domestic, captive domestic in research or industrial settings. We observed that domestication context explained much of the variation in gut microbiome composition, pathogen abundances and immune markers, yet the main differences in the repertoire of metabolic genes found in the gut microbiome were between the wild and domestic genetic lineages. We also documented population-level effects within domestication contexts, demonstrating that fine scale environmental variation also shaped host and microbe features. Our findings highlight that understanding which gut microbiome and immune traits respond to host genetic lineage and/or scales of local ecology could inform targeted interventions that manipulate the gut microbiome to achieve beneficial health outcomes.

Genetic parameters for automatically-measured vaginal temperature, respiration efficiency, and other thermotolerance indicators measured on lactating sows under heat stress conditions.

BACKGROUND: Genetic selection based on direct indicators of heat stress could capture additional mechanisms that are involved in heat stress response and enable more accurate selection for more heat-tolerant individuals. Therefore, the main objectives of this study were to estimate genetic parameters for various heat stress indicators in a commercial population of Landrace × Large White lactating sows measured under heat stress conditions. The main indicators evaluated were: skin surface temperatures (SST), automatically-recorded vaginal temperature (TV), respiration rate (RR), panting score (PS), body condition score (BCS), hair density (HD), body size (BS), ear size, and respiration efficiency (Reff). RESULTS: Traits based on TV presented moderate heritability estimates, ranging from 0.15 ± 0.02 to 0.29 ± 0.05. Low heritability estimates were found for SST traits (from 0.04 ± 0.01 to 0.06 ± 0.01), RR (0.06 ± 0.01), PS (0.05 0.01), and Reff (0.03 ± 0.01). Moderate to high heritability values were estimated for BCS (0.29 ± 0.04 for caliper measurements and 0.25 ± 0.04 for visual assessments), HD (0.25 ± 0.05), BS (0.33 ± 0.05), ear area (EA; 0.40 ± 0.09), and ear length (EL; 0.32 ± 0.07). High genetic correlations were estimated among SST traits (> 0.78) and among TV traits (> 0.75). Similarly, high genetic correlations were also estimated for RR with PS (0.87 ± 0.02), with BCS measures (0.92 ± 0.04), and with ear measures (0.95 ± 0.03). Low to moderate positive genetic correlations were estimated between SST and TV (from 0.25 ± 0.04 to 0.76 ± 0.07). Low genetic correlations were estimated between TV and BCS (from - 0.01 ± 0.08 to 0.06 ± 0.07). Respiration efficiency was estimated to be positively and moderately correlated with RR (0.36 ± 0.04), PS (0.56 ± 0.03), and BCS (0.56 ± 0.05 for caliper measurements and 0.50 ± 0.05 for the visual assessments). All other trait combinations were lowly genetically correlated. CONCLUSIONS: A comprehensive landscape of heritabilities and genetic correlations for various thermotolerance indicators in lactating sows were estimated. All traits evaluated are under genetic control and heritable, with different magnitudes, indicating that genetic progress is possible for all of them. The genetic correlation estimates provide evidence for the complex relationships between these traits and confirm the importance of a sub-index of thermotolerance traits to improve heat tolerance in pigs.

Evaluating phenotypes associated with heat tolerance and identifying moderate and severe heat stress thresholds in lactating sows housed in mechanically or naturally ventilated barns during the summer under commercial conditions.

An accurate understanding of heat stress (HS) temperatures and phenotypes that indicate HS tolerance is necessary to improve swine HS resilience. Therefore, the study objectives were 1) to identify phenotypes indicative of HS tolerance, and 2) to determine moderate and severe HS threshold temperatures in lactating sows. Multiparous (4.10 ± 1.48) lactating sows and their litters (11.10 ± 2.33 piglets/litter) were housed in naturally ventilated (n = 1,015) or mechanically ventilated (n = 630) barns at a commercial sow farm in Maple Hill, NC, USA between June 9 and July 24, 2021. In-barn dry bulb temperatures (TDB) and relative humidity were continuously recorded for naturally ventilated (26.38 ± 1.21 °C and 83.38 ± 5.40%, respectively) and mechanically ventilated (26.91 ± 1.80 °C and 77.13 ± 7.06%, respectively) barns using data recorders. Sows were phenotyped between lactation days 11.28 ± 3.08 and 14.25 ± 3.26. Thermoregulatory measures were obtained daily at 0800, 1200, 1600, and 2000 h and included respiration rate, and ear, shoulder, rump, and tail skin temperatures. Vaginal temperatures (TV) were recorded in 10 min intervals using data recorders. Anatomical characteristics were recorded, including ear area and length, visual and caliper-assessed body condition scores, and a visually assessed and subjective hair density score. Data were analyzed using PROC MIXED to evaluate the temporal pattern of thermoregulatory responses, phenotype correlations were based on mixed model analyses, and moderate and severe HS inflection points were established by fitting TV as the dependent variable in a cubic function against TDB. Statistical analyses were conducted separately for sows housed in mechanically or naturally ventilated barns because the sow groups were not housed in each facility type simultaneously. The temporal pattern of thermoregulatory responses was similar for naturally and mechanically ventilated barns and several thermoregulatory and anatomical measures were significantly correlated with one another (P < 0.05), including all anatomical measures as well as skin temperatures, respiration rates, and TV. For sows housed in naturally and mechanically ventilated facilities, moderate HS threshold TDB were 27.36 and 26.69 °C, respectively, and severe HS threshold TDB were 29.45 and 30.60 °C, respectively. In summary, this study provides new information on the variability of HS tolerance phenotypes and environmental conditions that constitute HS in commercially housed lactating sows.

Climate change and the associated increase in global temperatures have a well-described negative impact on swine production. Therefore, improving swine heat stress resilience is of utmost importance to reduce the deleterious effects of heat stress on swine health, performance, and welfare. Genomic selection for heat stress resilience may be a viable strategy to improve swine productivity in a changing climate. However, identifying environmental conditions that constitute heat stress and deriving novel traits that can be easily collected on farm and provide accurate and precise predictions of heat stress tolerance is a necessary step. The present study demonstrated that housing conditions had a limited influence on heat stress tolerance phenotypes, several anatomical and thermoregulatory measures were correlated, and housing conditions impacted heat stress threshold temperatures. Results from this study may be applied to large-scale phenotyping initiatives to develop or refine genomic selection indexes for heat stress resilience in pigs.

Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle.

Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (FPED), marker-based using a genomic relationship matrix (FGRM), and segment-based using ROH (FROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using FPED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using FGRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using FROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.

Gut microbiota and host genetics contribute to the phenotypic variation of digestive and feed efficiency traits in growing pigs fed a conventional and a high fiber diet.

BACKGROUND: Breeding pigs that can efficiently digest alternative diets with increased fiber content is a viable strategy to mitigate the feed cost in pig production. This study aimed at determining the contribution of the gut microbiota and host genetics to the phenotypic variability of digestive efficiency (DE) traits, such as digestibility coefficients of energy, organic matter and nitrogen, feed efficiency (FE) traits (feed conversion ratio and residual feed intake) and growth traits (average daily gain and daily feed intake). Data were available for 791 pigs fed a conventional diet and 735 of their full-sibs fed a high-fiber diet. Fecal samples were collected at 16 weeks of age to sequence the V3-V4 regions of the 16S ribosomal RNA gene and predict DE with near-infrared spectrometry. The proportions of phenotypic variance explained by the microbiota (microbiability) were estimated under three OTU filtering scenarios. Then, microbiability and heritability were estimated independently (models Micro and Gen) and jointly (model Micro+Gen) using a Bayesian approach for all traits. Breeding values were estimated in models Gen and Micro+Gen. RESULTS: Differences in microbiability estimates were significant between the two extreme filtering scenarios (14,366 and 803 OTU) within diets, but only for all DE. With the intermediate filtering scenario (2399 OTU) and for DE, microbiability was higher (> 0.44) than heritability (< 0.32) under both diets. For two of the DE traits, microbiability was significantly higher under the high-fiber diet (0.67 ± 0.06 and 0.68 ± 0.06) than under the conventional diet (0.44 ± 0.06). For growth and FE, heritability was higher (from 0.26 ± 0.06 to 0.44 ± 0.07) than microbiability (from 0.17 ± 0.05 to 0.35 ± 0.06). Microbiability and heritability estimates obtained with the Micro+Gen model did not significantly differ from those with the Micro and Gen models for all traits. Finally, based on their estimated breeding values, pigs ranked differently between the Gen and Micro+Gen models, only for the DE traits under both diets. CONCLUSIONS: The microbiota explained a significant proportion of the phenotypic variance of the DE traits, which was even larger than that explained by the host genetics. Thus, the use of microbiota information could improve the selection of DE traits, and to a lesser extent, of growth and FE traits. In addition, our results show that, at least for DE traits, filtering OTU is an important step and influences the microbiability.

Exploring methods to summarize gut microbiota composition for microbiability estimation and phenotypic prediction in swine.

The microbial composition resemblance among individuals in a group can be summarized in a square covariance matrix and fitted in linear models. We investigated eight approaches to create the matrix that quantified the resemblance between animals based on the gut microbiota composition. We aimed to compare the performance of different methods in estimating trait microbiability and predicting growth and body composition traits in three pig breeds. This study included 651 purebred boars from either breed: Duroc (n = 205), Landrace (n = 226), and Large White (n = 220). Growth and body composition traits, including body weight (BW), ultrasound backfat thickness (BF), ultrasound loin depth (LD), and ultrasound intramuscular fat (IMF) content, were measured on live animals at the market weight (156 ± 2.5 d of age). Rectal swabs were taken from each animal at 158 ± 4 d of age and subjected to 16S rRNA gene sequencing. Eight methods were used to create the microbial similarity matrices, including 4 kernel functions (Linear Kernel, LK; Polynomial Kernel, PK; Gaussian Kernel, GK; Arc-cosine Kernel with one hidden layer, AK1), 2 dissimilarity methods (Bray-Curtis, BC; Jaccard, JA), and 2 ordination methods (Metric Multidimensional Scaling, MDS; Detrended Correspondence analysis, DCA). Based on the matrix used, microbiability estimates ranged from 0.07 to 0.21 and 0.12 to 0.53 for Duroc, 0.03 to 0.21 and 0.05 to 0.44 for Landrace, and 0.02 to 0.24 and 0.05 to 0.52 for Large White pigs averaged over traits in the model with sire, pen, and microbiome, and model with the only microbiome, respectively. The GK, JA, BC, and AK1 obtained greater microbiability estimates than the remaining methods across traits and breeds. Predictions were made within each breed group using four-fold cross-validation based on the relatedness of sires in each breed group. The prediction accuracy ranged from 0.03 to 0.18 for BW, 0.08 to 0.31 for BF, 0.21 to 0.48 for LD, and 0.04 to 0.16 for IMF when averaged across breeds. The BC, MDS, LK, and JA achieved better accuracy than other methods in most predictions. Overall, the PK and DCA exhibited the worst performance compared to other microbiability estimation and prediction methods. The current study shows how alternative approaches summarized the resemblance of gut microbiota composition among animals and contributed this information to variance component estimation and phenotypic prediction in swine.

Gut microbiota has received significant research attention in farm animals because of its close relationship with host performance. We chose eight approaches to create a square covariance matrix that characterizes the relationship among animals based on their gut microbiota composition. Then, we fitted this information with linear models to evaluate the proportion of phenotypic variance explained by gut microbiota composition and predict host growth and body composition traits in three pig breeds. We found that different matrices had varying performance in predicting host phenotypes, but the results highly depended on the trait and breed considered in the prediction. Our findings highlight possible alternative approaches to incorporate gut microbiome data in regression models and emphasize the value of gut microbiome data in better understanding complex traits in pigs with diverse genetic backgrounds.

Genotyping and phenotyping strategies for genetic improvement of meat quality and carcass composition in swine.

BACKGROUND: Meat quality and composition traits have become valuable in modern pork production; however, genetic improvement has been slow due to high phenotyping costs. Combining genomic information with multi-trait indirect selection based on cheaper indicator traits is an alternative for continued cost-effective genetic improvement. METHODS: Data from an ongoing breeding program were used in this study. Phenotypic and genomic information was collected on three-way crossbred and purebred Duroc animals belonging to 28 half-sib families. We applied different methods to assess the value of using purebred and crossbred information (both genomic and phenotypic) to predict expensive-to-record traits measured on crossbred individuals. Estimation of multi-trait variance components set the basis for comparing the different scenarios, together with a fourfold cross-validation approach to validate the phenotyping schemes under four genotyping strategies. RESULTS: The benefit of including genomic information for multi-trait prediction depended on the breeding goal trait, the indicator traits included, and the source of genomic information. While some traits benefitted significantly from genotyping crossbreds (e.g., loin intramuscular fat content, backfat depth, and belly weight), multi-trait prediction was advantageous for some traits even in the absence of genomic information (e.g., loin muscle weight, subjective color, and subjective firmness). CONCLUSIONS: Our results show the value of using different sources of phenotypic and genomic information. For most of the traits studied, including crossbred genomic information was more beneficial than performing multi-trait prediction. Thus, we recommend including crossbred individuals in the reference population when these are phenotyped for the breeding objective.

Use of Milk Infrared Spectral Data as Environmental Covariates in Genomic Prediction Models for Production Traits in Canadian Holstein.

The purpose of this study was to provide a procedure for the inclusion of milk spectral information into genomic prediction models. Spectral data were considered a set of covariates, in addition to genomic covariates. Milk yield and somatic cell score were used as traits to investigate. A cross-validation was employed, making a distinction for predicting new individuals' performance under known environments, known individuals' performance under new environments, and new individuals' performance under new environments. We found an advantage of including spectral data as environmental covariates when the genomic predictions had to be extrapolated to new environments. This was valid for both observed and, even more, unobserved families (genotypes). Overall, prediction accuracy was larger for milk yield than somatic cell score. Fourier-transformed infrared spectral data can be used as a source of information for the calculation of the 'environmental coordinates' of a given farm in a given time, extrapolating predictions to new environments. This procedure could serve as an example of integration of genomic and phenomic data. This could help using spectral data for traits that present poor predictability at the phenotypic level, such as disease incidence and behavior traits. The strength of the model is the ability to couple genomic with high-throughput phenomic information.

Use of Host Feeding Behavior and Gut Microbiome Data in Estimating Variance Components and Predicting Growth and Body Composition Traits in Swine.

The purpose of this study was to investigate the use of feeding behavior in conjunction with gut microbiome sampled at two growth stages in predicting growth and body composition traits of finishing pigs. Six hundred and fifty-one purebred boars of three breeds: Duroc (DR), Landrace (LR), and Large White (LW), were studied. Feeding activities were recorded individually from 99 to 163 days of age. The 16S rRNA gene sequences were obtained from each pig at 123 ± 4 and 158 ± 4 days of age. When pigs reached market weight, body weight (BW), ultrasound backfat thickness (BF), ultrasound loin depth (LD), and ultrasound intramuscular fat (IMF) content were measured on live animals. Three models including feeding behavior (Model_FB), gut microbiota (Model_M), or both (Model_FB_M) as predictors, were investigated. Prediction accuracies were evaluated through cross-validation across genetic backgrounds using the leave-one-breed-out strategy and across rearing environments using the leave-one-room-out approach. The proportions of phenotypic variance of growth and body composition traits explained by feeding behavior ranged from 0.02 to 0.30, and from 0.20 to 0.52 when using gut microbiota composition. Overall prediction accuracy (averaged over traits and time points) of phenotypes was 0.24 and 0.33 for Model_FB, 0.27 and 0.19 for Model_M, and 0.40 and 0.35 for Model_FB_M for the across-breed and across-room scenarios, respectively. This study shows how feeding behavior and gut microbiota composition provide non-redundant information in predicting growth in swine.

Genotype by Environment Interaction and Selection Response for Milk Yield Traits and Conformation in a Local Cattle Breed Using a Reaction Norm Approach.

Local breeds are often reared in various environmental conditions (EC), suggesting that genotype by environment interaction (GxE) could influence genetic progress. This study aimed at investigating GxE and response to selection (R) in Rendena cattle under diverse EC. Traits included milk, fat, and protein yields, fat and protein percentage, and somatic cell score, three-factor scores and 24 linear type traits. The traits belonged to 11,085 cows (615 sires). Variance components were estimated in a two-step reaction norm model (RNM). A single trait animal model was run to obtain the solutions of herd-EC effect, then included in a random regression sire model. A multivariate response to selection (R) in different EC was computed for traits under selection including beef traits from a performance test. GxE accounted on average for 10% of phenotypic variance, and an average rank correlation of over 0.97 was found between bull estimated breeding values (EBVs) by either including or not including GxE, with changing top ranks. For various traits, significantly greater genetic components and R were observed in plain farms, loose housing rearing system, feeding total mixed ration, and without summer pasture. Conversely, for beef traits, a greater R was found for mountain farms, loose housing, hay-based feeding and summer pasture.

Exploring the role of gut microbiota in host feeding behavior among breeds in swine.

BACKGROUND: The interplay between the gut microbiota and feeding behavior has consequences for host metabolism and health. The present study aimed to explore gut microbiota overall influence on feeding behavior traits and to identify specific microbes associated with the traits in three commercial swine breeds at three growth stages. Feeding behavior measures were obtained from 651 pigs of three breeds (Duroc, Landrace, and Large White) from an average 73 to 163 days of age. Seven feeding behavior traits covered the information of feed intake, feeder occupation time, feeding rate, and the number of visits to the feeder. Rectal swabs were collected from each pig at 73 ± 3, 123 ± 4, and 158 ± 4 days of age. DNA was extracted and subjected to 16 S rRNA gene sequencing. RESULTS: Differences in feeding behavior traits among breeds during each period were found. The proportion of phenotypic variances of feeding behavior explained by the gut microbial composition was small to moderate (ranged from 0.09 to 0.31). A total of 21, 10, and 35 amplicon sequence variants were found to be significantly (q-value < 0.05) associated with feeding behavior traits for Duroc, Landrace, and Large White across the three sampling time points. The identified amplicon sequence variants were annotated to five phyla, with Firmicutes being the most abundant. Those amplicon sequence variants were assigned to 28 genera, mainly including Christensenellaceae_R-7_group, Ruminococcaceae_UCG-004, Dorea, Ruminococcaceae_UCG-014, and Marvinbryantia. CONCLUSIONS: This study demonstrated the importance of the gut microbial composition in interacting with the host feeding behavior and identified multiple archaea and bacteria associated with feeding behavior measures in pigs from either Duroc, Landrace, or Large White breeds at three growth stages. Our study provides insight into the interaction between gut microbiota and feeding behavior and highlights the genetic background and age effects in swine microbial studies.

Genetic diversity and population history of eight Italian beef cattle breeds using measures of autozygosity.

In the present study, GeneSeek GGP-LDv4 33k single nucleotide polymorphism chip was used to detect runs of homozygosity (ROH) in eight Italian beef cattle breeds, six breeds with distribution limited to Tuscany (Calvana, Mucca Pisana, Pontremolese) or Sardinia (Sarda, Sardo Bruna and Sardo Modicana) and two cosmopolitan breeds (Charolais and Limousine). ROH detection analyses were used to estimate autozygosity and inbreeding and to identify genomic regions with high frequency of ROH, which might reflect selection signatures. Comparative analysis among breeds revealed differences in length and distribution of ROH and inbreeding levels. The Charolais, Limousine, Sarda, and Sardo Bruna breeds were found to have a high frequency of short ROH (~ 15.000); Calvana and Mucca Pisana presented also runs longer than 16 Mbp. The highest level of average genomic inbreeding was observed in Tuscan breeds, around 0.3, while Sardinian and cosmopolitan breeds showed values around 0.2. The population structure and genetic distances were analyzed through principal component and multidimensional scaling analyses, and resulted in a clear separation among the breeds, with clusters related to productive purposes. The frequency of ROH occurrence revealed eight breed-specific genomic regions where genes of potential selective and conservative interest are located (e.g. MYOG, CHI3L1, CHIT1 (BTA16), TIMELESS, APOF, OR10P1, OR6C4, OR2AP1, OR6C2, OR6C68, CACNG2 (BTA5), COL5A2 and COL3A1 (BTA2)). In all breeds, we found the largest proportion of homozygous by descent segments to be those that represent inbreeding events that occurred around 32 generations ago, with Tuscan breeds also having a significant proportion of segments relating to more recent inbreeding.