Background There is an increasing need to account for genotype-by-environment (G × E) interactions in livestock breeding programs to improve productivity and animal welfare across environmental and management conditions. This is even more relevant for pigs because selection occurs in high-health nucleus farms, while commercial pigs are raised in more challenging environments. In this study, we used single-step homoscedastic and heteroscedastic genomic reaction norm models (RNM) to evaluate G × E interactions in Large White pigs, including 8686 genotyped animals, for reproduction (total number of piglets born, TNB; total number of piglets born alive, NBA; total number of piglets weaned, NW), growth (weaning weight, WW; off-test weight, OW), and body composition (ultrasound muscle depth, MD; ultrasound backfat thickness, BF) traits. Genetic parameter estimation and single-step genome-wide association studies (ssGWAS) were performed for each trait. Results The average performance of contemporary groups (CG) was estimated and used as environmental gradient in the reaction norm analyses. We found that the need to consider heterogeneous residual variance in RNM models was trait dependent. Based on estimates of variance components of the RNM slope and of genetic correlations across environmental gradients, G × E interactions clearly existed for TNB and NBA, existed for WW but were of smaller magnitude, and were not detected for NW, OW, MD, and BF. Based on estimates of the genetic variance explained by the markers in sliding genomic windows in ssGWAS, several genomic regions were associated with the RNM slope for TNB, NBA, and WW, indicating specific biological mechanisms underlying environmental sensitivity, and dozens of novel candidate genes were identified. Our results also provided strong evidence that the X chromosome contributed to the intercept and slope of RNM for litter size traits in pigs. Conclusions We provide a comprehensive description of G × E interactions in Large White pigs for economically-relevant traits and identified important genomic regions and candidate genes associated with GxE interactions on several autosomes and the X chromosome. Implementation of these findings will contribute to more accurate genomic estimates of breeding values by considering G × E interactions, in order to genetically improve the environmental robustness of maternal-line pigs.

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single-and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from −0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Properly quantifying environmental heat stress (HS) is still a major challenge in livestock breeding programs, especially as adverse climatic events become more common. The definition of critical periods and climatic variables to be used as the environmental gradient is a key step for genetically evaluating heat tolerance (HTol). Therefore, the main objectives of this study were to define the best critical periods and environmental variables (ENV) to evaluate HT and estimate variance components for HT in Large White pigs. The traits included in this study were 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). Seven climatic variables based on public weather station data were compared based on three criteria, including the following: (1) strongest G×E estimate as measured by the slope term, (2) ENV yielding the highest theoretical accuracy of the genomic estimated breeding values (GEBV), and (3) variable yielding the highest distribution of GEBV per ENV. Relative humidity (for BFT, MDP, NBD, WN, and WW) and maximum temperature (for OTW, TNB, NBA, IBF, and IWE) are the recommended ENV based on the analyzed criteria. The acute HS (average of 30 days before the measurement date) is the critical period recommended for OTW, BFT, and MDP in the studied population. For WN, WW, IBF, and IWE, a period ranging from 34 days prior to farrowing up to weaning is recommended. For TNB, NBA, and NBD, the critical period from 20 days prior to breeding up to 30 days into gestation is recommended. The genetic correlation values indicate that the traits were largely (WN, WW, IBF, and IWE), moderately (OTW, TNB, and NBA), or weakly (MDP, BFT, and NBD) affected by G×E interactions. This study provides relevant recommendations of critical periods and climatic gradients for several traits in order to evaluate HS in Large White pigs. These observations demonstrate that HT in Large White pigs is heritable, and genetic progress can be achieved through genetic and genomic selection.

analysis of pig growth curves combining pedigree and genomic information, Livestock Science,

Genome association analyses have been successful in identifying quantitative trait loci (QTLs) for pig body weights measured at a single age. However, when considering the whole weight trajectories over time in the context of genome association analyses, it is important to look at the markers that affect growth curve parameters. The easiest way to consider them is via the two-step method, in which the growth curve parameters and marker effects are estimated separately, thereby resulting in a reduction of the statistical power and the precision of estimates. One efficient solution is to adopt nonlinear mixed models (NMM), which enables a joint modeling of the individual growth curves and marker effects. Our aim was to propose a genome association analysis for growth curves in pigs based on NMM as well as to compare it with the traditional two-step method. In addition, we also aimed to identify the nearest candidate genes related to significant SNP (single nucleotide polymorphism) markers. The NMM presented a higher number of significant SNPs for adult weight (A) and maturity rate (K), and provided a direct way to test SNP significance simultaneously for both the A and K parameters. Furthermore, all significant SNPs from the two-step method were also reported in the NMM analysis. The ontology of the three candidate genes (SH3BGRL2, MAPK14, and MYL9) derived from significant SNPs (simultaneously affecting A and K) allows us to make inferences with regards to their contribution to the pig growth process in the population studied.

This research aimed to estimate genetic parameters for milk yield and type traits [withers height (WH), croup height (CH), body length (BL), croup length (CL), iliac width (ILW), ischial width (ISW), and thoracic circumference] in Murrah buffaloes and to identify genomic regions related to type traits by applying a single-step genome-wide association study. Data used to estimate the genetic parameters consisted of 601 records of milk yield in the first lactation and the aforementioned type traits. For the single-step genome-wide association study, 322 samples genotyped with a 90K Axiom Buffalo Genotyping array (Thermo Fisher Scientific, Santa Clara, CA) were used. Bivariate analysis revealed that heritability for milk yield (kg) at 305 d was 0.31 ± 0.11, whereas it ranged from 0.22 ± 0.07 to 0.34 ± 0.09 for the studied conformation traits. Based on the percentages of genetic variance explained by windows of 10 markers, there were 16 genomic regions explaining more than 0.5% of the variance for WH, CH, BL, CL, ILW, ISW, and thoracic circumference. Between those regions, 4 were associated with more than 1 trait, suggesting pleiotropic roles for some genes of Bos taurus autosome (BTA) 12 on CL and WH, BTA13 on ISW and ILW, BTA23 on CH and BL, and BTA28 on ISW and BL. Most of these regions coincide with known quantitative trait loci for milk traits. Thus, further studies based on sequence data will help to validate the association of this region with type traits and likely identify the causal mutations.

We compared different Bayesian models to handle censored data for genetic parameters estimation of age at first calving (AFC) in Brazilian Brahman cattle. Data from females with AFC above 1825 days of age were assumed to have failed to calve and were considered as censored records. Data including information of 53,703 cows were analyzed through the following methods: conventional linear model method (LM), which consider only uncensored records; simulation method (SM), in which the data were augmented by drawing random samples from positive truncated normal distributions; penalty method (PM), in which a constant of 21 days was added to censored records; and the bivariate threshold-linear method (TLcens). The LM was the most suited for genetic evaluation of AFC in Brazilian Brahman cattle based on the predictive ability evaluation through cross-validation analysis. The similar results for LM and PM regarding Spearman correlations, and the higher percentages of selected animals in common, indicated that there was not relevant reranking of animals when censored records were used. In summary, the heritability estimates for AFC ranged from 0.09 (TLcens) to 0.20 (LM). Given its poor predictive performance, the SM is not recommended for handling censored records for genetic evaluation of AFC.

RESUMOObjetivou-se com este estudo estimar parâmetros genéticos para o número total de leitões nascidos (NTLN), número de leitões nascidos vivos (NLNV) e número de leitões vivos aos cinco dias de idade (NLV5) com modelos de regressão aleatória e averiguar melhor modelagem da variância residual na avaliação das trajetórias genéticas do tamanho da leitegada de fêmeas Landrace e Large White. Os dados utilizados foram provenientes de uma granja de melhoramento genético de suínos e continham 2.388 observações de fêmeas Landrace e 2.325 de Large White. Os modelos de melhor ajuste para o NTLN e NLV5 foram os que consideraram a variância residual homogênea e, para NLNV, o modelo com quatro classes de variâncias residuais foi o mais adequado (BIC). Para Landrace, o efeito materno não foi significativo. O modelo que incluiu o efeito materno e quatro classes de variância residual foi o que apresentou melhor ajuste para NTLN na raça Large White, sendo os modelos sem efeito materno e com variância residual homogênea os mais adequados para NLNV e NLV5. As herdabilidades estimadas variaram de baixas a altas (0,08-0,34; 0,04-0,29 e 0,05-0,21 na raça Landrace e 0,16-0,30; 0,10-0,37 e 0,09-0,32 na Large White, para NTLN, NLNV e NLV5, respectivamente). A alta correlação de posto entre os valores genéticos do NLNV e NLV5 sugere que não há necessidades do controle do NLV5 nesse programa de melhoramento genético. Maiores ganhos podem ser obtidos pela seleção no NLNV de fêmeas primíparas, em função da diminuição do intervalo de gerações.Palavras-chave: leitões, correlação, herdabilidade, modelo ABSTRACT This study aimed to estimate genetic parameters for total number of piglets born (NTLN), number of piglets born alive (NLNV) and number of piglets alive at five days of age (NLV5) using random regression models and to evaluate the best way for modelling the residual variance in the description of the genetic trajectories of litter size in Landrace and Large White breeds. The data came from a swine breed improvement program, and a total of 2388 and 2325 litter size records of Landrace and Large White, respectively were used in the analyses. The models considering homogenous residual variance showed the best goodness of fit for NTLN and NLV5 and the model with four classes of residual variances was most appropriate for NLNV (BIC). In the Landrace breed the maternal effect was not significant. The model including maternal effect and four classes of residual variance adequately described NTLN of Large White breed and the models without maternal effectand with homogeneous residual variance were most appropriate to describe NLNV and NLV5. The estimated heritability for NTLN,