Survival analysis of productive life in Florida dairy goats using a Cox proportional hazards model.

Longevity is an economically important trait, since extending the functional life of a doe would allow us to keep the most productive females in the herd as long as possible, and this could result in the increased profitability of dairy farms. Thus, the objectives of this study were to determine the most important factors that influence the length of productive life (LPL) of female Florida goats and to estimate its genetic additive variance using a Cox proportional hazards model. The data consisted of 70,695 productive life records from 25,722 Florida females kidding between 2006 and 2020. A total of 19,495 does had completed their productive life while 6227 (24.2%) does had censored information. The pedigree contained information on 56,901 animals. The average censoring age and average failure age after first kidding for LPL were 36 and 47 months respectively. The model included, as time-independent effects, the age at first kidding and the interaction between herd, year and season of birth of the doe, and as time-dependent effects, the age at kidding, the interaction between herd, year and season of kidding, the within-herd class of milk production deviation, and the interaction between the lactation number and the stage of lactation. All fixed effects had a significant effect on LPL (p < 0.05). Does with older ages at the first kidding and an earlier age at kidding were at higher risk of being culled. A large difference among herds was observed in terms of culling risk, which highlighted the importance of adequate management practices. Also, high-producing does were less likely to be culled. The estimate of the additive genetic variance was 1.844 (in genetic standard deviation), with a heritability estimate of 0.58 ± 0.012. The results of this study are expected to contribute to the development of a genetic model for genetic evaluation of the length of the productive life of Spanish dairy goat breeds.

Molecular detection of Rickettsia in ectoparasites (Siphonaptera and Phthiraptera) of domestic and feral pigs from Argentina.

Rickettsioses are distributed among a variety of hematophagous arthropods, and represent an emergent threat. The presence of rickettsial bacteria in ectoparasites collected from pigs from Argentina is still unknown. This study investigated the presence and identity of Rickettsia spp. in fleas, Pulex irritans, and sucking lice, Haematopinus suis, of domestic and feral pigs, Sus scrofa, from Central-Northern Argentina, through the genes gltA and ompB. Rickettsial bacteria were detected in 50% of fleas and 24% of lice. The BLASTn analysis of the ompB gene fragments in P. irritans samples showed identities 99% and 100% with R. felis. Positive samples of H. suis were 99% similar with species from the spotted fever group, future amplifications of a more polymorphic fragment of the ompB gene will allow to corroborate the identity of the Rickettsia species present in these lice samples. The Rickettsia spp. reported in the present study are having eventually been associated with cases of human diseases, and the circulation of these agents in arthropods has already been reported in several countries. Therefore, the identification of circulating pathogenic agents, such as reported in this study, is crucial for development of preventive measures for the control of ectoparasite-borne rickettsiosis diseases. Further studies, using serology techniques, will be allow to explore the ability of pigs as a possible Rickettsia reservoir and its role as part of transmission cycle of Rickettsia spp. in the studied scenarios.

Genome-wide association study for feed efficiency in collective cage-raised rabbits under full and restricted feeding.

Feed efficiency (FE) is one of the most economically and environmentally relevant traits in the animal production sector. The objective of this study was to gain knowledge about the genetic control of FE in rabbits. To this end, GWASs were conducted for individual growth under two feeding regimes (full feeding and restricted) and FE traits collected from cage groups, using 114 604 autosome SNPs segregating in 438 rabbits. Two different models were implemented: (1) an animal model with a linear regression on each SNP allele for growth trait; and (2) a two-trait animal model, jointly fitting the performance trait and each SNP allele content, for FE traits. This last modeling strategy is a new tool applied to GWAS and allows information to be considered from non-genotyped individuals whose contribution is relevant in the group average traits. A total of 189 SNPs in 17 chromosomal regions were declared to be significantly associated with any of the five analyzed traits at a chromosome-wide level. In 12 of these regions, 20 candidate genes were proposed to explain the variation of the analyzed traits, including genes such as FTO, NDUFAF6 and CEBPA previously associated with growth and FE traits in monogastric species. Candidate genes associated with behavioral patterns were also identified. Overall, our results can be considered as the foundation for future functional research to unravel the actual causal mutations regulating growth and FE in rabbits.

Pedigree-based estimation of covariance between dominance deviations and additive genetic effects in closed rabbit lines considering inbreeding and using a computationally simpler equivalent model.

Inbreeding generates covariances between additive and dominance effects (breeding values and dominance deviations). In this work, we developed and applied models for estimation of dominance and additive genetic variances and their covariance, a model that we call "full dominance," from pedigree and phenotypic data. Estimates with this model such as presented here are very scarce both in livestock and in wild genetics. First, we estimated pedigree-based condensed probabilities of identity using recursion. Second, we developed an equivalent linear model in which variance components can be estimated using closed-form algorithms such as REML or Gibbs sampling and existing software. Third, we present a new method to refer the estimated variance components to meaningful parameters in a particular population, i.e., final partially inbred generations as opposed to outbred base populations. We applied these developments to three closed rabbit lines (A, V and H) selected for number of weaned at the Polytechnic University of Valencia. Pedigree and phenotypes are complete and span 43, 39 and 14 generations, respectively. Estimates of broad-sense heritability are 0.07, 0.07 and 0.05 at the base versus 0.07, 0.07 and 0.09 in the final generations. Narrow-sense heritability estimates are 0.06, 0.06 and 0.02 at the base versus 0.04, 0.04 and 0.01 at the final generations. There is also a reduction in the genotypic variance due to the negative additive-dominance correlation. Thus, the contribution of dominance variation is fairly large and increases with inbreeding and (over)compensates for the loss in additive variation. In addition, estimates of the additive-dominance correlation are -0.37, -0.31 and 0.00, in agreement with the few published estimates and theoretical considerations.

Role of inbreeding depression, non-inbred dominance deviations and random year-season effect in genetic trends for prolificacy in closed rabbit lines.

In closed rabbit lines selected for prolificacy at the Polytechnic University of Valencia, genetic responses are predicted using BLUP. With a standard additive BLUP model and year-season (YS) effects fitted as fixed, genetic trends were overestimated compared to responses estimated using control populations obtained from frozen embryos. In these lines, there is a confounding between genetic trend, YS effects and inbreeding, and the role of dominance is uncertain. This is a common situation in data from reproductively closed selection lines. This paper fits different genetic evaluation models to data of these lines, aiming to identify the source of these biases: dominance, inbreeding depression and/or an ill-conditioned model due to the strong collinearity between YS, inbreeding and genetic trend. The study involved three maternal lines (A, V and H) and analysed two traits, total born (TB) and the number of kits at weaning (NW). Models fitting YS effect as fixed or random were implemented, in addition to additive genetic, permanent environment effects and non-inbred dominance deviations effects. When YS was fitted as a fixed effect, the genetic trends were overestimated compared to control populations, inbreeding had an apparent positive effect on litter size and the environmental trends were negative. When YS was fitted as random, the genetic trends were compatible with control populations results, inbreeding had a negative effect (lower prolificacy) and environmental trends were flat. The model fitting random YS, inbreeding and non-inbred dominance deviations yielded the following ratios of additive and dominance variances to total variance for NW: 0.06 and 0.01 for line A, 0.06 and 0.00 for line V and 0.01 and 0.08 for line H. Except for line H, dominance deviations seem to be of low relevance. When it is confounded with inbreeding as in these lines, fitting YS effect as random allows correct estimation of genetic trends.

Growth traits of four maternal lines of rabbits founded on different criteria: comparisons at foundation and at last periods after selection.

The objective was to compare growth traits in four maternal lines of rabbits (A, V, H and LP), with the aim of understanding the consequence of the different foundation and selection processes on the growth performance of the lines. The lines are currently in the 43th, 38th, 22th and 8th generations, respectively. Two comparisons were performed. One compared the values of the lines at their foundation, using the complete data set, the full pedigree and a two-trait analysis, including data on the selection criteria, litter size. The other comparisons were done during the last period when all the lines were housed together with the same feeding and management. The numbers of records were 323 208 for weaning weight, and 300 553 for slaughter weight and average daily Gain (from 46 708 l). The pedigree file included 346 638 animals. The second analysis used only the data corresponding to each period, and the analysis was conducted using a one-trait model. The model was the same as that defined for the comparisons at the foundation, but the additive effects were excluded. The H and LP lines showed highest values for all the traits compared. In the last periods, a good agreement was observed between the estimated differences, computed with the complete model and data set, or computed with an incomplete model and only data from the comparison period. At last periods, the differences were smaller than at foundation. The importance of the correlated response in growth after selection for litter size at weaning or the importance of a non-programmed intramating selection for the growth traits can explain the changes since foundation.

Effective population size and inbreeding depression on litter size in rabbits. A case study.

The purpose of this study is to use demographic and litter size data on four Spanish maternal lines of rabbits (A, V, H and LP), as a case study, in order to: (i) estimate the effective population size of the lines, as a measure of the rate of increase of inbreeding, and (ii) study whether the inbreeding effect on litter size traits depends on the pattern of its accumulation over time. The lines are being selected for litter size at weaning and are kept closed at the same selection nucleus under the same selection and management programme. The study considered 47,794 l and a pedigree of 14,622 animals. Some practices in mating and selection management allow an increase of the inbreeding coefficient lower than 0.01 per generation in these lines of around 25 males and 125 females. Their effective population size (Ne) was around 57.3, showing that the effect of selection, increasing the inbreeding, was counterbalanced by the management practices, intended to reduce the rate of inbreeding increase. The inbreeding of each individual was broken down into three components: old, intermediate and new inbreeding. The coefficients of regression of the old, intermediate and new inbreeding on total born (TB), number born alive (NBA) and number weaned (NW) per litter showed a decreasing trend from positive to negative values. Regression coefficients significantly different from zero were those for the old inbreeding on TB (6.79 ± 2.37) and NBA (5.92 ± 2.37). The contrast between the coefficients of regression between the old and new inbreeding were significant for the three litter size traits: 7.57 ± 1.72 for TB; 6.66 ± 1.73 for NBA and 5.13 ± 1.67 for NW. These results have been interpreted as the combined action of purging unfavourable genes and artificial selection favoured by the inbreeding throughout the generations of selection.

Novel phenotypes of feeding and social behaviour and their relationship with individual rabbit growth and feed efficiency.

Due to the lack of a recording system for individual consumption of group-housed rabbits, published studies about feeding behaviour are based on information recorded at the group- and not at the individual level and periods covering only a few days or, in some cases, only part of a day. Such information could be used to inform rabbit management systems but cannot be used for genetic selection. We aimed to generate and use information from a novel automated feeder for group-housed rabbits to identify new phenotypes for individual animals that could be incorporated into breeding programs to improve feed efficiency and social behaviour under different feeding regimens. At 39 d of age, rabbits from 15 batches were placed in cages and fed ad libitum to become used to the electronic feeder. From 42 to 58-59 d, one group of 1 086 rabbits was fed ad libitum (AL), while another group of 1 134 rabbits was fed on a restricted feeding schedule (R) by limiting the feeding time to the period between 1800 and 0600 h of the following day. We implemented a reliable multivariate method to remove anomalous feeding behaviour records. We then defined novel traits for feeding behaviour that apply to both types of feeding regimes, and for social behaviour that indicates an animal's rank within the cage hierarchy. We based these traits on feeder records and a biologically sound definition of a meal. Finally, we estimated the phenotypic correlations of those traits with growth and feed efficiency traits. Our findings demonstrate that variables about resource distribution among cage mates and an animal's priority for feed access were found to be good indicators of an animal's dominant or subordinate status within the cage. Based on results obtained in R animals (results were similar in AL animals), the most efficient animals were those that ate less frequently (phenotypic correlation with feed conversion ratio, rho = 0.6), and consumed smaller amounts per meal (rho = 0.7), spent less time at the feeder (rho = 0.4), and appeared to be subordinate, as they did not have priority access to the feeder (rho = -0.3), and had the smallest share of resources (range of rho = 0.2-0.6). We conclude that quantifying feeding and social behaviour traits can enhance the understanding of the mechanisms through which individuals exert their effects on the performance of their cage mates.

Replication and refinement of a quantitative trait locus influencing milk protein percentage on ovine chromosome 3.

A previous genome scan that was conducted in Spanish Churra sheep identified a significant quantitative trait locus (QTL) for milk protein percentage (PP) on chromosome 3 (OAR3), between markers KD103 and OARVH34. The aim of this study was to replicate these results and to refine the mapped position of this QTL. To accomplish this goal, we analysed 14 new half-sib families of Spanish Churra sheep including 1661 ewes from 29 different flocks. These animals were genotyped for 21 microsatellite markers mapping to OAR3. In addition to a classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) was performed with the aim of enhancing the resolution of the QTL mapping. The LA that was performed in this sheep population identified the presence of a highly significant QTL for PP near marker KD103 (P(c) < 0.001; P(exp) < 0.001). The phenotypic variance that was owing to the QTL was 2.74%. Two segregating families for the target QTL were identified in this population with QTL effect estimates of 0.47 and 0.95 SD. The LDLA identified the same QTL as the previous analyses with a high level of statistical significance (P = 9.184 E-11) and narrowed the confidence interval (CI) to a 13 cM region. These results confirm the segregation of the previously identified OAR3 QTL that influences PP in Spanish Churra sheep. Future research will aim to increase the marker density across the refined CI and to analyse the corresponding candidate genes to identify the allelic variant or variants that underlie this genetic effect.

The role of bovine causal genes underlying dairy traits in Spanish Churra sheep.

In dairy cattle, quantitative trait nucleotides (QTNs) underlying quantitative trait loci (QTL) for milk production traits have been identified in bovine DGAT1, GHR and ABCG2 genes. The SPP1 gene has also been proposed to be a regulator of lactation. In sheep, QTL underlying milk production traits have been reported only recently, and no proven QTN has been identified. Taking into account the close phylogenetic relationship between sheep and cattle, this study examined the possible effects of the aforementioned genes on sheep milk production traits. We first studied the genetic variability of the DGAT1, GHR, ABCG2 and SPP1 genes in 15 rams of the Spanish Churra dairy sheep breed. Second, we performed an association analysis between SNPs identified in these genes and three milk production traits recorded in a commercial population of Churra sheep. This analysis revealed only three significant associations at the nominal level (P-value <0.05) involving allelic variants of the ABCG2 gene, whereas no significant association was found for the DGAT1, GHR and SPP1 genes. When the Bonferroni correction was applied to take into account the multiple tests performed, none of the associations identified at the nominal level remained significant. Nevertheless, taking into account the high level of false-negative findings that can arise when applying the stringent Bonferroni correction, we think that our results provide a valuable primary assessment of strong candidate genes for milk traits in sheep.

A genome scan for quantitative trait loci affecting body conformation traits in Spanish Churra dairy sheep.

A genome scan for chromosomal regions influencing body conformation traits was conducted for a population of Spanish Churra dairy sheep following a daughter design. A total of 739 ewes from 11 half-sib sire families were included in the study. The ewes were scored for the 5 linear traits used in the breeding scheme of the Churra breed to assess body conformation: stature, rear legs-rear view, foot angle, rump width, and general appearance. All the animals, including the 11 sires, were genotyped for 181 microsatellite markers evenly distributed across the 26 sheep autosomes. Using the yield deviations of the raw scores adjusted for fixed factors as phenotypic measurements, a quantitative trait loci (QTL) analysis was performed on the basis of a multi-marker regression method. Seven suggestive QTL were identified on chromosomes Ovis aries (OAR)2, OAR5, OAR16, OAR23, and OAR26, but none reached a genome-wise significance level. Putative QTL were identified for all of the traits analyzed, except for general appearance score. The suggestive QTL showing the highest test statistic influenced rear legs-rear view and was localized on OAR16, close to the growth hormone receptor coding gene, GHR. Some of the putative linkage associations reported here are consistent with previously reported QTL in cattle for similar traits. To the best of our knowledge, this study provides the first report of QTL for body conformation traits in dairy sheep; further studies will be needed to confirm and redefine the linkage associations reported herein. It is expected that future genome-wide association analyses of larger families will help identify genes underlying these putative genetic effects and provide useful markers for marker-assisted selection of such functional traits.

The identification of QTL that affect the fatty acid composition of milk on sheep chromosome 11.

In this work, we analysed 11 genetic markers localized on OAR11 in a commercial population of Spanish Churra sheep to detect QTL that underlie milk fatty acid (FA) composition traits. Following a daughter design, we analysed 799 ewes distributed in 15 half-sib families. Eight microsatellite markers and three novel SNPs identified in two genes related to fatty acid metabolism, acetyl-CoA carboxylase alpha (ACACA) and fatty acid synthase (FASN), were genotyped in the whole population under study. The phenotypic traits considered in the study included 22 measurements related to the FA composition of the milk and three other milk production traits (milk protein percentage, milk fat percentage and milk yield). Across-family regression analysis revealed four significant QTL at the 5% chromosome-wise level influencing contents of capric acid (C10:0), lauric acid (C12:0), linoleic conjugated acid (CLA) and polyunsaturated fatty acids (PUFA) respectively. The peaks of the QTL affecting C10:0 and PUFA contents in milk map close to the FASN gene, which has been evaluated as a putative positional candidate for these QTL. The QTL influencing C12:0 content reaches its maximum significance at 58 cM, close to the gene coding for the glucose-dependent insulinotropic polypeptide. We were not able to find any candidate genes related to fat metabolism at the QTL influencing CLA content, which is located at the proximal end of the chromosome. Further research efforts will be needed to confirm and refine the QTL locations reported here.

Genetic determination of fatty acid composition in Spanish Churra sheep milk.

The objective of this study was to estimate the genetic variation of ovine milk fatty acid (FA) composition. We collected 4,100 milk samples in 14 herds from 976 Churra ewes sired mostly by 15 AI rams and analyzed them by gas-liquid chromatography for milk fatty acid composition. The studied traits were 12 individual FA contents (proportion in relation to the total amount of FA), 3 groups of fatty acids [saturated fatty acids (SFA), monounsaturated FA (MUFA), and polyunsaturated FA (PUFA)], and 2 FA ratios (n-6:n-3 and C18:2 cis-9,trans-11:C18:1 trans-11). In addition, percentages of fat and protein and daily milk yield were studied. For the analysis, repeatability animal models were implemented using Bayesian methods. In an initial step, univariate methods were conducted to test the hypothesis of the traits showing additive genetic determination. Deviance information criterion and Bayes factor were employed as model choice criteria. All the studied SFA showed additive genetic variance, but the estimated heritabilities were low. Among unsaturated FA (UFA), only C18:1 trans-11 and C18:2 cis-9,cis-12 showed additive genetic variation, their estimated heritabilities being [marginal posterior mean (marginal posterior SD)] 0.02(0.01) and 0.11(0.04), respectively. For the FA groups, only PUFA showed significant additive genetic variation. None of the studied ratios of FA showed additive genetic variation. In second multitrait analyses, genetic correlations between individual FA and production traits, and between groups of FA and ratios of FA and production traits, were investigated. Positive genetic correlations were estimated among medium-chain SFA, ranging from 0 to 0.85, but this parameter was close to zero between long-chain SFA (C16:0 and C18:0). Between long- and medium-chain SFA, estimated genetic correlations were negative, around -0.6. Among those UFA showing significant additive genetic variance, genetic correlations were close to zero. The estimated genetic correlations among all the investigated FA, milk yield, and fat and protein percentages were not different from zero. Our results suggest that low additive genetic variation is involved in the determination of the FA composition of milk fat in Churra sheep under current production conditions, which results in low values of heritabilities.

Detection of quantitative trait loci affecting the milk fatty acid profile on sheep chromosome 22: role of the stearoyl-CoA desaturase gene in Spanish Churra sheep.

Sheep milk fat contains several components that may provide human health benefits, such as monounsaturated fatty acids and conjugated linoleic acid (CLA). Most of the CLA in ruminant milk is synthesized in the mammary gland by the action of the enzyme stearoyl-CoA desaturase (SCD) on circulating vaccenic acid (trans-11 C18:2; VA). Previous studies have found significant associations between polymorphisms in the SCD gene and the fatty acid composition of ruminant products, including sheep milk. Based on this, we performed a quantitative trait loci (QTL) analysis of an ovine chromosome (22) that harbors the SCD gene for effects on milk fatty acid composition traits and classical milk production traits. We identified a suggestive QTL influencing the CLA/VA ratio with the maximum statistic at position 26 cM of the studied chromosome, whereas the SCD gene has been mapped to position 41.6 cM. The individual introduction of 4 SCD single nucleotide polymorphisms in the QTL model did not cause a reduction of the variance explained by the QTL, which suggests that the SCD gene is not directly responsible for the detected effect in the Churra population studied herein. This conclusion was supported by the lack of any significant association identified between the 4 SCD single nucleotide polymorphisms and the CLA/VA ratio. This association analysis suggested a possible effect of the SCD gene on milk fat percentage in Churra sheep. An independent confirmation of these primary results will be required before attempting its practical implementation in selection programs.

Anticipating the Impact of the USMLE Step 1 Pass/Fail Scoring Decision on Underrepresented-in-Medicine Students.

Three-digit United States Medical Licensing Examination (USMLE) Step 1 scores have assumed an outsized role in residency selection decisions, creating intense pressure for medical students to obtain a high score on this exam. In February 2020, the Federation of State Medical Boards and the National Board of Medical Examiners announced that Step 1 would transition to pass/fail scoring beginning in 2022.The authors discuss the potential advantages and disadvantages of the pass/fail scoring change for underrepresented-in-medicine (UiM) trainees. UiM students may benefit from this change because it reduces the effect of an inequitable exam; helps correct for students who attend medical schools with a curriculum heavier on nontested formative elements; and decreases stress, improves quality of life, and undermines imposter syndrome. However, this change may also precipitate unforeseen challenges, such as increased discrimination toward UiM trainees, an increase in high-stakes test failures due to a reduced focus on preparing for standardized exams, or the development of new (e.g., subject exams) or overreliance on existing (e.g., school ranking) metrics that would make UiM residency candidates less competitive.To enhance UiM representation in the future health care workforce, it is imperative that national organizations (e.g., accrediting, licensing, regulatory, professional, honor, student, and faculty), hospitals, residency programs, and patient advocacy groups undertake a shared, rigorous approach in assessing the impact of the pass/fail scoring change on UiM applicants' selection to specialty and subspecialty residencies.

Indirect genetic effects on the relationships between production and feeding behaviour traits in growing Duroc pigs.

Performance and feeding behaviour traits in growing pigs could be affected by social interaction effects when animals are raised in group. So, properly knowing the genetic correlations between direct and social interaction effects among performance and feeding behaviour traits could improve the accuracy of the genetic evaluations. Our aim was to explore the role of feeding behaviour traits (FBT) and indirect genetic effects (IGEs) in the genetic evaluations of growing pigs. Thus, genetic parameters were estimated for production traits (PT): average daily gain, average daily feed consumption, feed conversion ratio and backfat thickness; as well as for FBT: average daily feeding rate, average daily feeding frequency, average daily occupation time and average daily time between consecutive visits. Traits were recorded in 1144 Duroc pigs during the fattening period. Two bivariate models were fitted: classic animal model and an animal model fitting IGE. Estimations were done following Bayesian procedures. Heritability estimates obtained with classic animal model for all studied traits were medium-high. The additional heritable variation captured by IGE supposed that the ratios of total genetic variance to phenotypic variance (T2) were higher than the heritability estimates obtained with the classic model, except for occupation time trait, when a lower value (0.20 ± 0.19) was estimated. This is due to a high and negative correlation between IGE and direct genetic effects (DGEs) of this particular trait (-0.78 ± 0.27). Results from classic animal model do not evidence a clear role of FBT to improve the accuracy of breeding value predictions for PT; only average daily feeding rate seems to show a positive correlation (around 0.50 to 0.60) with average daily gain, average daily feed consumption and backfat thickness. However, when IGE model was fitted, the number of estimates of genetic correlations between FBT and PT showing a relevant magnitude increased, generally for the correlations between IGE of FBT and DGE of PT; or particularly for the correlations between IGE of average daily feeding frequency, and the IGE of all the PT, except average daily gain. Thus, in evaluations using the animal model with IGE fitted, the inclusion of FBT could aid the improvement of the accuracy of breeding value predictions for PT. This is a consequence of the improved genetic relationships between traits that can be fitted when considering such models.

Genetic parameters of sow feed efficiency during lactation and its underlying traits in a Duroc population.

As a result of the genetic selection for prolificacy and the improvements in the environment and farms management, litter size has increased in the last few years so that energy requirements of the lactating sow are greater. In addition, selection for feed efficiency of growing pigs is also conducted in maternal lines, and this has led to a decrease in appetite and feed intake that is extended to the lactation period, so the females are not able to obtain the necessary energy and nutrients for milk production and they mobilize their energetic reserves. When this mobilization is excessive, reproductive and health problems occur which ends up in an early sow culling. In this context, it has been suggested to improve feed efficiency at lactation through genetic selection. The aim of this study is to know, in a Duroc population, the genetic determinism of sow feed efficiency during lactation and traits involved in its definition, as well as genetic and environmental associations between them. The studied traits are daily lactation feed intake (dLFI), daily sow weight balance (dSWB), backfat thickness balance (BFTB), daily litter weight gain (dLWG), sow residual feed intake (RFI) and sow restricted residual feed intake (RRFI) during lactation. Data corresponded to 851 parities from 581 Duroc sows. A Bayesian analysis was performed using Gibbs sampling. A four-trait repeatability animal model was implemented including the systematic factors of batch and parity order, the standardized covariates of sow weight (SWf) and litter weight (LWs) at farrowing for all traits and lactation length for BFTB. The posterior mean (posterior SD) of heritabilities were: 0.09 (0.03) for dLFI, 0.37 (0.07) for dSWB, 0.09 (0.03) for BFTB, 0.22 (0.05) for dLWG, 0.04 (0.02) for RFI and null for RRFI. The genetic correlation between dLFI and dSWB was high and positive (0.74 (0.11)) and null between dLFI and BFTB. Genetic correlation was favourable between RFI and dLFI and BFTB (0.71 (0.16) and -0.69 (0.18)), respectively. The other genetic correlations were not statistically different from zero. The phenotypic correlations were low and positive between dLFI and dSWB (0.27 (0.03), dSWB and BFTB (0.25 (0.04)), and between dLFI and dLWG (0.16 (0.03)). Therefore, in the population under study, the improvement of the lactation feed efficiency would be possible either using RFI, which would not have unfavourable correlated effects, or through an index including its component traits.

Genetic variability of functional longevity in five rabbit lines.

The objectives of this study were to analyse the differences in the genetic determination of functional longevity in five Spanish lines of rabbits and to check how different systematic factors might affect this genetic determination. Four of the lines were maternal (lines A, V, H and LP), these lines were established selecting base generation animals according to different criteria, but in the subsequent generations all of them were selected for litter size at weaning. The other is the paternal line R, this line was constituted by selecting animals with an outstanding daily growth rate. The trait analysed, length of productive life, was the time in days between the date of the first positive pregnancy test and the date of culling or death of a doe. Four models extended from the Cox proportional hazard model were used to analyse data of each line separately and jointly. The complete model (Model 1) included the fixed effect of year-season (YS) combination, positive palpation order (OPP), that is, reproductive cycle, physiological status of the doe (PS) at service and number of kits born alive (NBA) in each kindling as time-dependent factors. The inbreeding coefficient was fitted as a continuous covariate and the animal's additive genetic effect was also fitted to the model (Model 1). The other models were identical to Model 1 but excluding OPP (Model 2) or PS (Model 3) or NBA (Model 4), which were explored to assess the consequence on additive variance estimates of not correcting for these animal-dependent factors. Estimated effective heritabilities of longevity were 0.07 ± 0.03, 0.03 ± 0.02, 0.14 ± 0.09, 0.05 ± 0.04, 0.02 ± 0.01 and 0.04 ± 0.01 for lines A, V, H, LP, R and for the merged data set, respectively. Removing the PS from the model led to an increase in the estimated additive genetic variance in all lines (0.17 ± 0.05, 0.05 ± 0.03, 0.29 ± 0.19, 0.29 ± 0.20, 0.07 ± 0.04 and 0.05 ± 0.02 for lines A, V, H, LP, R and the merged data set, respectively). The highest hazard of death and/or culling was observed during the first two parities and decreased as the order of parity progressed. Does non-pregnant-non-lactating had the highest risk of death or culling. The does that had zero kits born alive incurred the highest risk, and this risk decreased as the NBA increased. In conclusion, the consideration of longevity as selection criterion for the studied rabbit lines is not recommended.

Norway rat (Rattus norvegicus) ectoparasites in livestock production systems from central Argentina: Influencing factors on parasitism.

Haematophagous ectoparasites are worldwide vectors of many zoonotic bacterial diseases, both emerging and re-emerging, whose incidences are rising. Livestock development alters different environmental characteristics such as the microclimate of a site, changing the availability, density and susceptibility of the hosts to pathogens and vectors, indirectly influencing the spread and persistence of a disease within an ecosystem. The Norway rat (Rattus norvegicus), the most abundant vertebrate pest species found on livestock farms from Argentina, is a reservoir for several important zoonotic bacteria and may harbor ectoparasite species, which act as their vectors. Even though the Norway rat is widely known for its role as an ectoparasite host, the ecological characteristics of their ectoparasite communities and the related factors with parasitism on livestock farms have never been described. In the present study, we describe the ectoparasite community in Norway rats from central Argentina livestock farms, while also depicting the influencing factors on both ectoparasite occurrence and abundance. Ectoparasites were collected from rats captured in 20 sites from Buenos Aires province, between the winter of 2016 and the summer of 2018. A total of 1441 ectoparasite individuals were collected from 159 Norway rat individuals [Total ectoparasite prevalence = 69.2%; Mean ectoparasite specimen abundance (± CI) = 9.06 ±â€¯2.32 ectoparasite individuals per rat; Mean ectoparasite specimen intensity (±CI) = 13.10 ±â€¯3.08 ectoparasite individuals per infested rat found]. Ectoparasite assemblage consisted of four cosmopolitan species, recognized for their sanitary relevance: mites (Laelapidae: Laelaps nuttalli and Laelaps echidninus), lice (Polyplacidae: Polyplax spinulosa) and fleas (Pulicidae: Xenopsylla cheopis). We observed higher Norway rat abundance in sites related to higher ectoparasite occurrence and abundance frequencies on the rats. Additionally, ectoparasites were more abundant on rats in warm seasons and on male individuals, over female rats. Moreover, the geographical location of the studied sites influenced the ectoparasite assemblage structure observed on the rats. This study broadens the knowledge on the role of Norway rats as zoonotic ectoparasites hosts and analyzes the drivers influencing ectoparasite occurrence and abundance on the most populated region of Argentina, which is also the region with the most intensive livestock farming. Therefore, this survey may assist in evaluating potential risks for humans and generate effective sanitary control strategies for ectoparasite-borne infectious diseases.

Genetic determination of the onset of heat stress on daily milk production in the US Holstein cattle.

Existence of individual variation in the onset of heat stress for daily milk yield of dairy cows was assessed. Data included 353,376 test-day records of 38,383 first-parity Holsteins from a random sample of US herds. Three hierarchical models were investigated. Model 1 inferred the value of a temperature-humidity index (THI) at which mean yield began to decline as well as the extent of that decline. Model 2 assumed individual variation in yield decline beyond a common THI threshold. Model 3 additionally assumed individual variation for the onset of heat stress. Deviance information criteria indicated the superiority of model 3 over model 2. For model 2, genetic correlation between milk yield in the absence of heat stress and the THI threshold for heat stress was -0.4 (0.11) [marginal posterior mean (marginal posterior standard deviation)]. For model 3, genetic correlations were -0.53 (0.05) between milk yield and THI threshold and -0.62 (0.08) between milk yield and yield decay beyond the THI threshold. Total standard deviation (sum of additive genetic and permanent environmental standard deviations) for the THI threshold was 3.95 (0.06), and more than half of that variation had an additive genetic origin [56% (5%)]. Because of the high genetic correlation [0.95 (0.03)] between yield decay and THI threshold with model 3, using only one of them as a selection criterion for heat tolerance would modify the other in the desired direction.