Genetic parameters for linear conformation, stayability, performance and reproduction traits in German local Swabian-Hall landrace sows.

Improvements of longevity in sows are of overriding importance from an animal welfare as well as from an economic perspective, also in the local Swabian-Hall Landrace (SHL) breed. Consequently, the aim of the present study was a detailed assessment of linear conformation traits in the context of early predictions of sow stayability and to infer genetic (co)variance components between conformation traits with reproduction and performance traits. In this regard, we implemented a linear 9-point scoring system for conformation traits reflecting the body morphology and feet and leg characteristics in gilts. Linear conformation trait scoring included body length (BLEN), body height (BHEI), hind leg angularity (HLANG), hind leg pastern (HLPAS), hind cannon bone strength (HCBS), hind leg position (HLPOS), front leg position (FLPOS), claw length (CLEN) and muscling (MUSC) from 6042 SHL gilts before first insemination at the age of 180-200 days. For the same gilts, performance traits included average daily gain (ADG) and back fat thickness (BF) measured via ultrasound, reproduction traits included the number of piglets born alive (NBA) and the number of piglets weaned (NWEAN), and stayability was a binary trait reflecting sow survival after parity 2 (STAY_12), 3 (STAY_13) and 4 (STAY_14). For the estimation of variance components and heritabilities, single-trait linear animal models were defined for conformation, performance and reproduction traits, and single-trait generalized linear mixed models with a logit link function for binary stayability traits. Genetic covariances and correlations were inferred in series of bivariate runs for all combinations of conformation and remaining traits. The distribution of the assigned conformation scores reflected a Gaussian distribution, but a large fraction of records was assigned to the intermediate score classes 4, 5 and 6. Accordingly, the restricted utilization of the 9-point scale might explain the small genetic variances and small heritabilities for feet and leg traits in the range from 0.06 to 0.17. Heritabilities were larger for the linear traits reflecting body morphology, i.e., for BLEN (0.21) and BHEI (0.20). For conformation traits, the common litter environment explained up to 17% of the phenotypic trait variation. Genetic and especially phenotypic correlations between linear conformation traits with the reproduction traits NBA and NWEAN were close to zero. Heritabilities for STAY_12, STAY_13 and STAY_14 were 0.08, 0.07 and 0.05, respectively. Moderate genetic relationships were estimated between STAY_14 with some conformation traits. Especially high scores for BHEI and BLEN (i.e., long and large gilts) implied a decline in STAY_14 genetically (rg  = -0.24 and -0.53, respectively). Moderate genetic correlations were estimated between HLANG with STAY_14 (0.28), and between HCBS with STAY_12 (0.23). For most of the conformation traits with intermediate optimum, genetic correlations with STAY were close to zero, indicating improved longevity for gilts representing the population average with scores 4, 5 or 6, and suggesting the development of appropriate selection indices in this regard.

Multiple-trait and structural equation modelling approaches to infer genetic relationships between tail length and weight traits in Merinoland sheep.

Tail docking is routinely conducted in long-tailed sheep breeds to prevent flystrike infections, but it is not in agreement with legal guidelines and animal welfare issues. Selection on short tails is a sustainable alternative in this regard, but side effects on other breeding goal traits are unclear. In consequence, the present study aimed to estimate genetic parameters for tail length (TL) at birth, birth weight (BW), weaning weight (WW) and postweaning weight (PWW) at the slaughtering date considering single-trait (STM), multiple-trait (MTM) and structural equation models (SEM) with different random effects, and accordingly, different covariance structures. The SEM considered time-lagged recursive relationships among response variables in three different pathways. The first path pertained to the effect of TL on WW and of WW on PWW. The second path reflected the effect of BW on WW and of WW on PWW. The third path was the recursive effect of TL on PWW. The phenotypic data consisted of 2803 records for TL, 13,042 records for BW, 1556 records for WW and 3986 records for PWW from Merinoland lambs. Lambs were born in the period from 1995 to 2021 and kept at the university Gießen research station, Germany, with their naturally long tails. Genetic statistical model evaluation based on Bayesian and Akaike's information criteria suggested models simultaneously considering direct genetic, maternal genetic and maternal permanent environmental effects and respective covariances. For statistical models including the same random effects and covariance structures, SEM were superior over MTM. The direct heritability for TL from the best-fitting STM was 0.60 ± 0.08, indicating the potential for genetic reduction of tail length within a few generations. For growth traits, the direct heritabilities ranged from 0.16 ± 0.03 for BW to 0.31 ± 0.09 for PWW. The maternal heritabilities were 0.03 ± 0.03 for TL, 0.12 ± 0.02 for BW, 0.04 ± 0.03 for WW and 0.07 ± 0.03 for PWW, reflecting small, but the non-significant influence of uterine characteristics on the tail development. The direct genetic correlations between TL and all weight traits were positive and very similar to MTM and SEM but reflected antagonistic genetic relationships from a breeding perspective. Oppositely, the structural equation coefficients reflecting trait associations phenotypically were negative (favourable) for the time-lagged effects of TL on WW and on PWW. As an explanation, lambs with long and woolly tails have an increased risk for contamination with dirt and dust causing infections, which in turn impairs the body weight development. In conclusion, breeding on short tails should consider trait-associated environmental risk factors, for example, disease susceptibility, which can be mimicked via SEM approaches.

Intra- and Interspecies RNA-Seq Based Variants in the Lactation Process of Ruminants.

The RNA-Seq data provides new opportunities for the detection of transcriptome variants' single nucleotide polymorphisms (SNPs) in various species and tissues. Herein, milk samples from two sheep breeds and two cow breeds were utilized to characterize the genetic variation in the coding regions in three stages (before-peak (BP), peak (P), and after-peak (AP)) of the lactation process. In sheep breeds Assaf and Churra, 100,462 and 97,768, 65,996 and 62,161, and 78,656 and 39,245 variants were observed for BP, P, and AP lactation stages, respectively. The number of specific variants was 59,798 and 76,419, 11,483 and 49,210, and 104,033 and 320,817 in cow breeds Jersy and Kashmiri, respectively, for BP, P, and AP stages. Via the transcriptome analysis of variation in regions containing QTL for fat, protein percentages, and milk yield, we detected a number of pathways and genes harboring mutations that could influence milk production attributes. Many SNPs detected here can be regarded as appropriate markers for custom SNP arrays or genotyping platforms to conduct association analyses among commercial populations. The results of this study offer new insights into milk production genetic mechanisms in cow and sheep breeds, which can contribute to designing suitable breeding systems for optimal milk production.

Estimation of direct and maternal genetic effects and annotation of potential candidate genes for weight and meat quality traits in a genotyped outdoor dual-purpose cattle breed.

With regard to potential applications of genomic selection in small numbered breeds, we evaluated genomic models and focused on potential candidate gene annotations for weight and meat quality traits in the local Rotes Höhenvieh (RHV) breed. Traits included 6,003 birth weights (BWT), 5,719 200 d-weights (200dw), 4,594 365 d-weights (365dw), and 547 records for intramuscular fat content (IMF). A total of 581,304 SNP from 370 genotyped cattle with phenotypic records were included in genomic analyses. Model evaluations focused on single- and multiple-trait models with direct and with direct and maternal genetic effects. Genetic relationship matrices were based on pedigree (A-matrix), SNP markers (G-matrix), or both (H-matrix). Genome-wide association studies (GWASs) were carried out using linear mixed models to identify potential candidate genes for the traits of interest. De-regressed proofs (DRP) for direct and maternal genetic components were used as pseudo-phenotypes in the GWAS. Accuracies of direct breeding values were higher from models based on G or on H compared to A. Highest accuracies (> 0.89) were obtained for IMF with multiple-trait models using the G-matrix. Direct heritabilities with maternal genetic effects ranged from 0.62 to 0.66 for BWT, from 0.45 to 0.55 for 200dW, from 0.40 to 0.44 for 365dW, and from 0.48 to 0.75 for IMF. Maternal heritabilities for BWT, 200dW, and 365dW were in a narrow range from 0.21 to 0.24, 0.24 to 0.27, and 0.21 to 0.25, respectively, and from 0.25 to 0.65 for IMF. Direct genetic correlations among body weight traits were positive and favorable, and very similar from different models but showed a stronger variation with 0.31 (A), -0.13 (G), and 0.45 (H) between BWT and IMF. In gene annotations, we identified 6, 3, 1, and 6 potential candidate genes for direct genetic effect on BWT, 200dW, 365dW, and IMF traits, respectively. Regarding maternal genetic effects, four (SHROOM3, ZNF609, PECAM1, and TEX2) and two (TMEM182 and SEC11A) genes were detected as potential candidate genes for BWT and 365dW, respectively. Potential candidate genes for maternal effect on IMF were GRHL2, FGA, FGB, and CTNNA3. As the most important finding from a practical breeding perspective, a small number of genotyped RHV cattle enabled accurate breeding values for high heritability IMF.

Identification of candidate genes on the basis of SNP by time-lagged heat stress interactions for milk production traits in German Holstein cattle.

The aim of this study was to estimate genotype by time-lagged heat stress (HS) variance components as well as main and interaction SNP-marker effects for maternal HS during the last eight weeks of cow pregnancy, considering milk production traits recorded in the offspring generation. The HS indicator was the temperature humidity index (THI) for each week. A dummy variable with the code = 1 for the respective week for THI ≥ 60 indicated HS, otherwise, for no HS, the code = 0 was assigned. The dataset included test-day and lactation production traits from 14,188 genotyped first parity Holstein cows. After genotype quality control, 41,139 SNP markers remained for the genomic analyses. Genomic animal models without (model VC_nHS) and with in-utero HS effects (model VC_wHS) were applied to estimate variance components. Accordingly, for genome-wide associations, models GWA_nHS and GWA_wHS, respectively, were applied to estimate main and interaction SNP effects. Common genomic and residual variances for the same traits were very similar from models VC_nHS and VC_wHS. Genotype by HS interaction variances varied, depending on the week with in-utero HS. Among all traits, lactation milk yield with HS from week 5 displayed the largest proportion for interaction variances (0.07). For main effects from model GWA_wHS, 380 SNPs were suggestively associated with all production traits. For the SNP interaction effects from model GWA_wHS, we identified 31 suggestive SNPs, which were located in close distance to 62 potential candidate genes. The inferred candidate genes have various biological functions, including mechanisms of immune response, growth processes and disease resistance. Two biological processes excessively represented in the overrepresentation tests addressed lymphocyte and monocyte chemotaxis, ultimately affecting immune response. The modelling approach considering time-lagged genotype by HS interactions for production traits inferred physiological mechanisms being associated with health and immunity, enabling improvements in selection of robust animals.

Heat stress during late pregnancy and postpartum influences genetic parameter estimates for birth weight and weight gain in dual-purpose cattle offspring generations.

Impact of direct heat stress (HS) on genetic parameter estimates, i.e., HS close to the trait recording date, was verified in several previous studies conducted in dairy and beef cattle populations. The aim of the present study was to analyze the impact of time-lagged HS at different recording periods during late pregnancy (a.p.) and postpartum (p.p.) on genetic parameter estimates for birth weight (BWT) and weight gain traits (200 d- and 365 d-weight gain (200dg, 365dg)) in offspring of the dual-purpose cattle breed "Rotes Höhenvieh" (RHV). Furthermore, we estimated genetic correlations within traits across time-lagged climatic indicators, in order to proof possible genotype by environment interactions (G×E). Trait recording included 5,434 observations for BWT, 3,679 observations for 200dg and 2,998 observations for 365dg. Time-lagged climatic descriptors were classes for the mean temperature humidity index (mTHI) and number of HS days (nHS) from the following periods: 7 d-period a.p. (BWT), 56 d-period a.p., and 56 d-period p.p. (200dg and 365dg). Genetic parameters were estimated via 2-trait animal models, i.e., defining the same trait in different climatic environments as different traits. Genetic variances and heritabilities for all traits increased with increasing mTHI- and nHS-classes for all recording periods, indicating pronounced genetic differentiation with regard to time-lagged in utero HS and HS directly after birth. Similarly, in low mTHI- and nHS-classes indicating cold stress, genetic variances, and heritabilities were larger than for temperate climates. Genetic correlations substantially smaller than 0.80 indicating G × E were observed when considering same traits from mTHI- and nHS-classes in greater distance. Estimated breeding values (EBV) of the 10 most influential sires with the largest number of offspring records fluctuated across mTHI- and nHS-classes. Correlations between sire EBV for same traits from distant climatic classes confirmed the genetic correlation estimates. Sires displaying stable EBV with climatic alterations were also identified. Selection of those sires might contribute to improved robustness in the RHV outdoor population genetically.

Direct and maternal genetic components for body weight traits in Markhoz goat.

The aim of the study was to estimate the (co)variance components and genetic trends of direct and maternal genetic effects in Markhoz goat using random regression (RR) models. A total of 16,283 body weight trait records from birth to 400 days of age were obtained from 3476 Markhoz kids. Kids were from 223 bucks and 1217 does raised in the Markhoz goat breeding station in Sanandaj, Iran. The fixed effects were gender, year of birth (1992-2014), age of dam (2 to 7 years old), and type of birth (single or twin). RR models based on Legendre polynomials of age at weight recording included direct and maternal genetic and permanent environmental effects. Heterogeneous residual variances were considered along age of recording. According to the Akaike information criterion (AIC), the Bayesian information criterion (BIC) and the likelihood ratio test (LRT), adding covariance between direct and maternal genetic effects, resulted in significant improvement in the level of fit. After examining models with different order of fits, the most suitable model consisted of the fourth-order polynomial for direct genetic effect and the third-order polynomial for maternal genetic, individual permanent environmental, and maternal permanent environmental effects. The direct heritabilities ranged from 0.21 at birth to 0.55 at 400 days of age. The largest (0.25) and the lowest (0.17) maternal heritabilities were estimated at birth and 60 days of age, respectively. Generally, negative direct-maternal genetic correlations were estimated between different ages, which might cause fluctuations of direct-maternal genetic trends along the years. It was concluded that selection at 100 days of age could reduce the impact of the antagonist relationship between maternal and direct genetic effects.

Maternal and direct genetic parameters for tail length, tail lesions, and growth traits in pigs.

Tail length and tail lesions are the major triggers for tail biting in pigs. Against this background, 2 datasets were analyzed to estimate genetic parameters for tail characteristics and growth traits. Dataset 1 considered measurements for trait tail length (T-LEN) and for the growth traits birth weight (BW), weaning weight (WW), postweaning weight (PWW), and average daily gain (ADG) from 9,348 piglets. Piglets were born in the period from 2015 to 2018 and kept on the university Gießen research station. Dataset 2 included 4,943 binary observations from 1,648 pigs from the birth years 2016 to 2019 for tail lesions (T-LES) as indicators for nail necrosis, tail abnormalities, or tail biting. T-LES were recorded at 30 ± 7 d after entry for rearing (T-Les-1), at 50 ± 7 d after entry for rearing (end of the rearing period, T-LES-2), and 130 ± 20 d after entry for rearing (end of fattening period, T-LES-3). Genetic statistical model evaluation for dataset 1 based on Akaike's information criterion and likelihood ration tests suggested multiple-trait animal models considering covariances between direct and maternal genetic effects. The direct heritability for T-LEN was 0.42 (±0.03), indicating the potential for genetic selection on short tails. The maternal genetic heritability for T-LEN was 0.05 (±0.04), indicating the influence of uterine characteristics on morphological traits. The negative correlation between direct and maternal effects for T-LEN of -0.35 (±0.13), as well as the antagonistic relationships (i.e., positive direct genetic correlations in the range from 0.03 to 0.40) between T-LEN with the growth traits BW, WW, PWW, and ADG, complicate selection strategies and breeding goal definitions. The correlations between direct effects for T-LEN and maternal effects for breeding goal traits, and vice versa, were positive but associated with a quite large SE. The heritability for T-LES when considering the 3 repeated measurements was 0.23 (±0.04) from the linear (repeatability of 0.30) and 0.21 (±0.06; repeatability of 0.29) from the threshold model. The breeding value correlations between T-LES-3 with breeding values from the repeatability models were quite large (0.74 to 0.90), suggesting trait lesion recording at the end of the rearing period. To understand all genetic mechanisms in detail, ongoing studies are focusing on association analyses between T-LEN and T-LES, and the identification of tail biting from an actor's perspective.

Time-lagged and acute impact of heat stress on production and fertility traits in the local dual-purpose cattle breed "Rotes Höhenvieh" under pasture-based conditions.

Climate change causes rising temperatures and extreme weather events worldwide, with possible detrimental time-lagged and acute impact on production and functional traits of cattle kept in outdoor production systems. The aim of the present study was to infer the influence of mean daily temperature humidity index (mTHI) and number of heat stress days (nHS) from different recording periods on birth weight (BWT), 200 d- and 365 d-weight gain (200 dg, 365 dg) of calves, and on the probability of stillbirth (SB), and calving interval (CINT) of their dams. Data recording included 4,362 observations for BWT, 3,136 observations for 200 dg, 2,502 observations for 365 dg, 9,293 observations for the birth status, and 2,811 observations for CINT of the local dual-purpose cattle breed "Rotes Höhenvieh" (RHV). Trait responses on mTHI and nHS were studied via generalized linear mixed model applications with identity link functions for Gaussian traits (BWT, 200 dg, 365 dg, CINT) and logit link functions for binary SB. High mTHI and high nHS before autumn births had strongest detrimental impact on BWT across all antepartum- (a.p.) periods (34.4 ± 0.79 kg maximum). Prolonged CINT was observed when cows suffered heat stress (HS) before or after spring calvings, with maximum length of 391.6 ± 3.82 d (56 d a.p.-period). High mTHI and high nHS during the 42 d- and 56 d a.p.-period implied increased probabilities for SB. We found a significant (P < 0.05) seasonal effect on SB in model 3 across all a.p.-periods, with the highest probability in autumn (maximum of 5.4 ± 0.82% in the 7 d a.p.-period). Weight gains of calves (200 dg and 365 dg) showed strongest HS response for mTHI and nHS measurements from the long-term postnatal periods (42 d- and 56 d-periods), with minimum 200 dg of 194.2 ± 4.15 kg (nHS of 31 to 42 d in the 42 d-period) or minimum 365 dg of 323.8 ± 3.82 kg (mTHI ≥ 60 in the 42 d-period). Calves born in summer, combined with high mTHI or high nHS pre- or postnatal, had lower weight gains, compared with calves born in other calving seasons or under cooler conditions. Highest BWT, weight gains, and shortest CINT mostly were detected under cool to moderate climate conditions for mTHI, and small to moderate nHS. Results indicate acute and time-lagged HS effects and address possible HS-induced epigenetic modifications of the bovine genome across generations and limited acclimatization processes to heat, especially when heat occurs during the cooler spring and autumn months.

The role of agro-ecological factors and transboundary transhumance in shaping the genetic diversity in four indigenous cattle populations of Benin.

The indigenous cattle population of Benin is a diverse mix of taurine and hybrid breeds shaped by diverse ecological and climatic conditions with eight agro-ecological zones (AEZ). Presumably, the taurine breeds face current endangerment due to ongoing indicine introgression following climate change and transboundary transhumance. The aim of the study was to investigate the genetic diversity and population structure of the indigenous breeds Lagune, Somba, Pabli and Borgou considering spatial agro-ecological and socio-economic factors (transhumance) based on 50k SNP and microsatellite data. Among the four sampled breeds, six genetic clusters were identified using model-free (discriminant analysis of principal components) and model-based (TESS and ADMIXTURE) methods separating taurine from hybrid breeds. Results based on an extension with publicly available historic SNP data sets from taurine and indicine West African cattle and additional outgroups provided additional insight into changes of genetic structure in the sampled breeds over time. Both taurine breeds, Somba and Lagune, showed a stable foundation but also spatially limited partial indicine introgression associated with transhumance leading to high genetic diversity. In addition, we found evidence for spatial diversity and changes in genetic structure over time in the Borgou breed in comparison of our samples with the historic samples which could be explained by potential continuous indicine introgression into the Borgou breed in two sample regions. Results for the Pabli breed do not conclusively point to full absorbance by the Borgou in comparison with all available Borgou samples. Further research is needed in this regard.

Variation in autosomal and sex-linked genetic effects for growth traits in Markhoz goat using multivariate animal models.

The present study was conducted to estimate autosomal and sex-linked genetic parameters by restricted maximum likelihood method using four different multivariate models for growth traits in Markhoz goat. Data were collected over a period of 23 years (1992-2014) from the breeding station of Markhoz goat. The studied growth traits were birth weight (BW), weaning weight (WW), 6-month weight (6WT), 9-month weight (9WT), and 12-month weight (12WT). The best model was determined by Akaike's information criterion. Based on the best model, direct autosomal heritabilities of 0.18, 0.10, 0.12, 0.22, and 0.24 were obtained for BW, WW, 6WT, 9WT, and 12WT, respectively. On the other hand, estimates of direct sex-linked heritability for BW, WW, 6WT, 9WT, and 12WT were 0.02, 0.01, 0.002, 0.03, and 0.04, respectively. The maternal permanent environmental component was remarkable for growth traits. Direct autosomal additive genetic correlations among the traits were positive and varied from 0.11 between BW and 9WT to 0.95 between BW and WW. The positive correlations were also estimated for direct sex-linked additive genetic effect except between BW-6WT (- 0.07) and WW-6WT (-0.19). The result of this study indicated that sex-linked chromosome effect had more considerable influence at post-weaning traits.