Search

Linkage disequilibrium (LD) and persistence of phase are fundamental approaches for exploring the genetic basis of economically important traits in cattle, including the identification of QTL for genomic selection and the estimation of effective population size (Ne) to determine the size of the training populations. In this study, we have used the Illumina BovineHD chip in 168 trios of 7 Spanish beef cattle breeds to obtain an overview of the magnitude of LD and the persistence of LD phase through the physical distance between markers. Also, we estimated the time of divergence based on the persistence of the LD phase and calculated past Ne from LD estimates using different alternatives to define the recombination rate. Estimates of average r² (as a measure of LD) for adjacent markers were close to 0.52 in the 7 breeds and decreased with the distance between markers, although in long distances, some LD still remained (0.07 and 0.05 for markers 200 kb and 1 Mb apart, respectively). A panel with a lower boundary of 38,000 SNP would be necessary to launch a successful within-breed genomic selection program. Persistence of phase, measured as the pairwise correlations between estimates of r in 2 breeds at short distances (10 kb), was in the 0.89 to 0.94 range and decreased from 0.33 to 0.52 to a range of 0.01 to 0.08 when marker distance increased from 200 kb to 1 Mb, respectively. The magnitude of the persistence of phase between the Spanish beef breeds was similar to those found in dairy breeds. For across-breed genomic selection, the size of the SNP panels must be in the range of 50,000 to 83,000 SNP. Estimates of past Ne showed values ranging from 26 to 31 for 1 generation ago in all breeds. The divergence among breeds occurred between 129 and 207 generations ago. The results of this study are relevant for the future implementation of within- and across-breed genomic selection programs in the Spanish beef cattle populations. Our results suggest that a reduced subset of the SNP panel would be enough to achieve an adequate precision of the genomic predictions. [ABSTRACT FROM AUTHOR]

The availability of SNP chips for massive genotyping has proven to be useful to genetically characterize populations of domestic cattle and to assess their degree of divergence. In this study, the Illumina BovineHD BeadChip genotyping array was used to describe the genetic variability and divergence among 7 important autochthonous Spanish beef cattle breeds. The within-breed genetic diversity, measured as the marker expected heterozygosity, was around 0.30, similar to other European cattle breeds. The analysis of molecular variance revealed that 94.22% of the total variance was explained by differences within individuals whereas only 4.46% was the result of differences among populations. The degree of genetic differentiation was small to moderate as the pairwise fixation index of genetic differentiation among breeds (FST) estimates ranged from 0.026 to 0.068 and the Nei's D genetic distances ranged from 0.009 to 0.016. A neighbor joining (N-J) phylogenetic tree showed 2 main groups of breeds: Pirenaica, Bruna dels Pirineus, and Rubia Gallega on the one hand and Avileña-Negra Ibérica, Morucha, and Retinta on the other. In turn, Asturiana de los Valles occupied an independent and intermediate position. A principal component analysis (PCA) applied to a distance matrix based on marker identity by state, in which the first 2 axes explained up to 17.3% of the variance, showed a grouping of animals that was similar to the one observed in the N-J tree. Finally, a cluster analysis for ancestries allowed assigning all the individuals to the breed they belong to, although it revealed some degree of admixture among breeds. Our results indicate large within-breed diversity and a low degree of divergence among the autochthonous Spanish beef cattle breeds studied. Both N-J and PCA groupings fit quite well to the ancestral trunks from which the Spanish beef cattle breeds were supposed to derive. [ABSTRACT FROM AUTHOR]

Demographic and pedigree analyses describe the structure and dynamics of livestock populations. We studied information recorded in the herdbooks of Asturiana de los Valles (AV; N = 458,806), Avileña-Negra Ibérica (ANI; N = 204,623), Bruna deis Pirineus (BP; N= 62,138), Morucha (Mo; N = 65,350), Pirenaica (Pi; N = 217,428), Retinta (Re; N = 135,300), and Rubia Gallega (RG; N= 235,511) beef breeds from their creation until 2009. All breeds have increased in the number of registered cows in recent years. In all breeds, herds do not behave as isolated entities and a high rate of exchange of breeding males between herds exists. A percentage of herds ( 12-52%) make some type of selection and sell bulls to other herds. There were large differences in average number of progeny per bull, ranging from 15.6 (AV) to 373.7 animals (RG, with a high incidence of AI). Generation interval estimates ranged from 4.7 (AV) to 7.6 (RG) yr in the sire pathway and from 5.95 (AV) to 7.8 (Mo) yr in the dam pathway. Density of pedigrees varied among breeds, with Pi, ANI, and Re having the more dense pedigrees, with average completeness indexes of more than 96% in the first generation and 80% when 6 generations were considered. A general increase in average inbreeding was observed in all breeds in the years analyzed. For animals bom in 2009, average inbreeding coefficients ranged from 0.6 (BP) to 7.2% (Re) when all animals were considered and from 3.6 (Pi) to 17.6% (BP) when only inbred animals were considered. Due to the lack of completeness of pedigrees in most populations, inbreeding coefficients may be considered as a lower bound of the true parameters. The proportion of inbred animals tended to increase in the periods analyzed in all breeds. Differences between inbreeding and coancestry rates (except in RG) suggest the presence of population structure. Effective population size (Ng) based on the inbreeding rate estimated by regression ranged from 43 to 378 for Re and BP, whereas N estimates based on coancestry were greater, with a range of 100 for RG to 9,985 for BP. These facts suggest that an adequate mating policy can help to monitor inbreeding so as not to lose genetic variability. Effective number of ancestors in 2009 for 6 of the breeds ranged from 42 (RG) to 220 (AV), with BP having much a greater value, and was lower than was the effective number of founders in all breeds, suggesting the existence of bottlenecks. [ABSTRACT FROM AUTHOR]

Although heteroskedasticity has been a main topic of interest in beef cattle during recent decades, the current availability of canalization models provided new insights for animal breeding programs. Within this context, birth BW (BWT) was analyzed in the Bruna dels Pirineus beef cattle breed by implementing canalization models that accounted for heterogeneous residual variances due to systematic, permanent environmental effects and additive genetic effects. Analyses were performed on BWT data from 8,130 calves born in 12 commercial breeding herds contributing to the yield recording scheme of the Bruna dels Pirineus breed. Analytical models accounted for direct additive genetic, permanent environmental, and 4 systematic effects (i.e., age of the dam, sex of the calf, birth type, and herd-year-season), and the same effects were evaluated as potential sources of variation in the residual term. Their relevance was checked by the deviance information criterion (DIC), and only residual additive genetic, permanent environmental, birth type, and herd-year-season remained in the operational model, all of them originating relevant reductions in the DIC parameter. Bruna dels Pirineus calves showed a moderate heritability of 0.30 (95% high posterior density, 0.19 to 0.40) for BWT; additional additive genetic variability was revealed in the residual term, this being positively correlated with the direct additive genetic component (0.44; 95% high posterior density, 0.37 to 0.54). Genetic trends were evaluated on both sources of additive genetic variance, and relevant patterns were identified in several herds. Although this breed did not evidence a homogeneous genetic trend for the whole population, herd-specific positive and negative trends were revealed, suggesting the plausibility of genetic selection for canalization on BWT in beef cattle breeds. These results must be viewed as a contribution to the canalization research field, providing relevant information for the breeding scheme of the Bruna dels Pirineus breed, as well as important insights about the genetic background of BWT for the beef industry worldwide. [ABSTRACT FROM AUTHOR]

The aim of this study was to assess the accuracy of ultrasound measurements for predicting carcass traits in 124 Spanish pascual-type lambs (13 to 16 kg carcass weight). Ultrasound images were taken transversal and longitudinal to the vertebral column and at thoracic (TV; between 12th and 13th ribs) and lumbar (LV; between first and second lumbar vertebrae) locations. Skin thickness, subcutaneous backfat thickness (BFT), and depth (DLD), width (WLD), and area (ALD) of longissimus dorsi were obtained with ImageJ 1.42q software. After slaughter, BFT (TV, 2.30 ± 0.06 mm; LV, 2.46 ± 0.06 mm), DLD (TV, 2.47 ± 0.03 cm; LV, 2.48 ± 0.03 cm), WLD (TV, 4.50 ± 0.04 cm; LV, 4.60 ± 0.04 cm), and ALD (TV, 9.96 ±0.12 cm²; LV, 10.19 ± 0.13 cm2) were directly measured on the lamb carcass. Correlations between ultrasound and direct carcass measurements were greater than 0.61 for DLD, WLD, and ALD (P < 0.05) whereas they fluctuated between 0.32 and 0.60 for BFT (P < 0.05); moreover, correlations were significantly (P < 0.05) greater for transversal than for longitudinal views. In a similar way, linear regression analyses suggested a moderate underestimation for BFT and lumbar DLD when using real-time ultrasound technologies whereas WLD, ALD, and thoracic DLD suffered from under-and overestimation for small and large values of carcass traits, respectively. After decomposing the mean square prediction error (MSPE) for the different ultrasound measurements, we found that the error due to disturbance contributed most to the MSPE followed by the error of central tendency and the error due to regression. The SE of prediction (SEP) was also calculated as an additional precision indicator, obtaining estimates less than that in previous studies with larger lambs. In conclusion, transversal ultrasound measurements at the thoracic and lumbar levels could be a useful tool for predicting DLD, WLD, and ALD in light lambs, perhaps suffering from worse prediction properties when focusing on BFT. This information could be of special relevance for light lamb producers worldwide, with a special emphasis in the Mediterranean basin where this kind of production system accounts for a large percentage of the sheep industry. [ABSTRACT FROM AUTHOR]

The aim of this study was to compare the goodness of fit of the threshold models with homoscedasticity or heteroscedasticity and the grouped data model for the analysis of calving ease in beef cattle by using a parametric bootstrap procedure. Field data included 8,205 records of the Bruna dels Pirineus beef cattle breed in the Pyrenean mountain areas of Catalonia (Spain). The actual distribution was 81.81% of calvings without assistance, 11.02% slightly assisted by the farmer, 5.12% strongly assisted by the farmer, 0.89% assisted by the veterinarian, and 1.16% cesarean, but these percentages were very different in the different herds. This can be explained partially by the different subjective way of scoring of each farmer. Primiparous cows had a greater (P < 0.001) difficulty calving than cows with 5 or more parities (11.74 vs. 4.49% of calvings strongly assisted by the farmer or the veterinarian and 2.8 vs. 0.65% cesarean). Male calves caused greater (P < 0.001) calving difficulty than females (7.71% of male calvings strongly assisted by the farmer or the veterinarian vs. 4.25% of females and 1.83% cesarean in males vs. 0.47% in females). The month and year of calving also had a strong influence on calving ease. These data were analyzed using 3 different models: the threshold models with homoscedasticity or heteroscedasticity and the grouped data model. The bootstrap comparison among models suggested that the threshold models, even allowing for heteroscedasticity, did not fit the herd effects well. In contrast, fitting deficiencies were not observed for the grouped data model in any factor. The variance of direct effect of the calf was estimated using the 3 models, and the heritability estimate ranged from 0.165 for the grouped data model to 0.185 for the hereroscedastic threshold model. This heritability was moderate, but it would justify the inclusion of direct effects of the calf on calving ease in the breeding objective. Overall, results highlighted the flexibility of the grouped data model for the analysis of discrete traits, like calving ease of beef calves. [ABSTRACT FROM AUTHOR]

Little is known about mutational variability in livestock, among which only a few mutations with relatively large effects have been reported. In this manuscript, mutational variability was analyzed in 1,765 litter size records from 404 Ripollesa ewes to characterize the magnitude of this genetic source of variation and check the suitability of including mutational effects in genetic evaluations of this breed. Threshold animal models accounting for additive genetic mutations were preferred to models without mutational contributions, with an average difference in the deviance information criterion of more than 5 units. Moreover, the statistical relevance of the additive genetic mutation term was checked through a Bayes factor approach, which showed that the models with mutational variability were 8.5 to 22.7 times more probable than the others. The mutational heritability (percentage of the phenotypic variance accounted for by mutational variance) was 0.6 or 0.9%, depending on whether genetic dominance effects were accounted for by the analytical model. The inclusion of mutational effects in the genetic model for evaluating litter size in Ripollesa ewes called for some minor modifications in the genetic merit order of the individuals evaluated, which suggested that the continuous uploading of new additive mutations could be taken into account to optimize the selection scheme. This study is the first attempt to estimate mutational variances in a livestock species and thereby contribute to better characterization of the genetic background of productive traits of interest. [ABSTRACT FROM AUTHOR]

Although inbreeding (F) is a topic of major concern in animal breeding, estimates of inbreeding depression are usually obtained by modeling the overall F coefficient of each individual, without considering that the recessive (deleterious) genetic load of a given population may be unevenly distributed among the founder genomes. The founder-specific partial F coefficient is calculated as the identity-by-descent probability at any given autosomal locus related to a particular founder and allows a more detailed analysis of inbreeding depression on productive traits. Within this context, birth BW data from 2,459 Ripollesa lambs were analyzed under a hierarchical animal model without F-related covariates (model 0), with inbreeding depression modeled by the overall F coefficient (model F1), or by the partial F coefficient of 9 founders that made a relevant contribution to the population inbreeding (model F2). A straightforward empirical Bayes factor (BF) was developed for testing statistical relevance of each F-related covariate, in which greater-than-1 values favored the model including the covariate. The deviance information criterion (DIC) clearly supported model F1 (5,767.8) rather than model 0 (5,771.2), suggesting that inbreeding depression had a relevant influence on birth BW data. The linear effect of inbreeding depression was statistically relevant in model F1 (BF = 2.52 × 1035), with lamb birth BW declining by -13.6 g with each 1% F increase. The quadratic effect of inbreeding depression was almost null in model F1 (BF = 0.02), as suggested by the reduction in DIC (5,766.9) when this effect was removed from model F1. On the other hand, model F2 provided a similar DIC (5,767.9) value, with this parameter decreasing to 5,764.7 when nonrelevant founder-specific inbreeding depression effects were removed. Substantial heterogeneity in founder-specific inbreeding depression was reported by model F2, in which estimates for 4 of the 9 founders did not differ from zero (BF between 0.05 and 0.42), whereas 5 founders originated moderate (-8.2 g for each 1% F increase; BF = 1.42) to large inbreeding depression (-96.2 g for each 1% F increase; BF = 8.80 × 1019). The substantial variability between founder estimates suggested that inbreeding depression effects may mainly be due to a few alleles with major deleterious effects. These results contribute valuable information that should help to achieve more accurate management of inbreeding in the Ripollesa breed. [ABSTRACT FROM AUTHOR]

A performance data set from 376 Ripollesa purebred ewes of the experimental flock of the Universitat Autônoma of Barcelona was analyzed using a bivariate Bayesian threshold-linear model. The data set contained 1,598 litter size records and 1,699 days-to-lambing records. The model included the additive genetic effect of each animal and 3 nongenetic sources of variation: ewe age, year of lambing, and the permanent environmental effect characterized by the ewe. The flock was phenotypically selected for litter size since 1986, and replacement ewes and rams were selected from the progeny of the more prolific ewes, which had at least 3 deliveries recorded. The phenotypic trend for litter size was positive, whereas days to lambing followed an unclear pattern. Both traits had low heritabilities; 0.13 for litter size and 0.11 for days to lambing. Response to selection was evaluated through (a) the average breeding value of the ewe lambs chosen annually, and (b) the average breeding value of the overall flock. The first measurement suggested a positive trend for litter size, although it showed important oscillations. On the other hand, the average breeding value for the overall flock showed a stable positive tendency after yr 4 of selection, with estimates clearly different from zero after yr 11 of selection. A significant increase in the incidence of multiple births was observed, with a mode of approximately 10%. The correlated response in days to lambing did not show a significant trend. The effect of year of lambing also positively influenced both litter size and days to lambing, although important oscillations were observed between years. Results indicated that litter size in sheep can be effectively improved through phenotypic selection, even in small flocks; moreover, days to lambing could also be genetically improved, given the estimate obtained for its heritability. [ABSTRACT FROM AUTHOR]

The analysis focused on model fitting of 2 ewe reproductive traits, litter size, and days to lambing (interval between the introduction of the ram into the flock and the subsequent parturition of the ewes). The experimental data set of the Universitat Autònoma of Barcelona flock was used, including 1,598 records of litter size and 1,699 records of days to lambing from 376 Ripollesa ewes between 1986 and 2005. Univariate and bivariate models were considered as beginning points with linear or threshold approximation for litter size. Model fitting was evaluated in terms of goodness-of-fit and predictive ability, using the mean square error and the correlation between phenotypic and predicted records (ρyŷ) as reference parameters. The bivariate model was preferable for both variables, minimizing mean square error and maximizing ρyŷ. A threshold approximation for litter size was preferable over a linear approximation. Models were also compared with a simulation study, comparing the correlation coefficient between simulated and predicted breeding values (ρa,â). The bivariate threshold model was favored, with a ρyŷ of 0.677 and 0.834 for litter size and days to lambing, respectively. Correlation coefficients between simulated and predicted breeding values in the bivariate linear model were reduced slightly to 0.651 and 0.831, respectively, and they were lowest with linear univariate models (0.642 and 0.802). Although the bivariate models for ewe litter size and days to lambing were more accurate than the univariate models, the threshold approaches showed a greater advantage under the bivariate model. For the purpose of genetic evaluation of litter size in sheep, use of the threshold-linear model seems justified. In the Ripollesa breed, the evaluation of litter size can benefit from recording birth weight. [ABSTRACT FROM AUTHOR]

The objective of this study was to analyze the association between the haplotypes of the prion protein (PrP) locus and several reproductive and lamb weight traits in Ripollesa sheep. Prion protein genotypes were available for a total of 310 sheep (7 rams, 114 ewes, and 189 lambs), all of them belonging to the purebred Ripollesa flock of the Universitat Autònoma of Barcelona, for which all sheep had a known pedigree. In addition, the genotype of 24 historical descendants of the previously genotyped adult individuals was reconstructed, provided that both parents were homozygous for PrP haplotypes. Only 3 haplotypes (ARR, ARQ, and ARH) were observed in the PrP locus of the sheep sampled. Reproductive traits included conception rate and litter size, whereas birth BW and 90-d BW were the lamb weight traits studied. The additive effect of PrP haplotypes was analyzed through Bayesian animal threshold and linear models, for reproduction and weight traits, respectively. Ewe reproductive data belonged to 89 ewes that gave 492 conception rate records and 440 litter size records. Analyses of BW at birth and at 90 d of age were made on 323 and 164 lamb records, respectively. No associations between PrP haplotypes and conception rate and BW traits were observed. For litter size, the effect of the ARH haplotype was greater than that of the ARQ haplotype. Differences between ARH and ARR haplotypes also suggested an advantage for the ARH. As a whole, our results indicated that the selection favorable to increase litter size in Ripollesa ewes may also increase the ARH haplotype frequency, which contradicts the recommendations of the current European Union legislation aiming to increase the genetic resistance to scrapie. As a consequence, scrapie genotyping needs to be included as a new selection criterion in the breed. [ABSTRACT FROM AUTHOR]

The survivability from birth to slaughter of 1,487 Ripollesa lambs with a preslaughter overall mortality of 9.6% was studied under the proportional hazards framework, assuming a Weibull distribution for the baseline hazards function. A sire frailty model was fitted, with the common environment received by the lamb as an additional random source of variation. Common environment was considered time-dependent and was characterized by the dam and the contemporary lamb group during the preweaning and fattening periods, respectively. Only 3 fixed effects were statistically significant: the linear and quadratic effects of birth weight (P < 0.001), the relative position of the delivery within the lambing season (P < 0.001), and the presence of stillbirths or mummified fetuses within the litter (P < 0.05). Birth type and parity of the ewe were significant only when birth weight was removed from the model (P < 0.001 and P < 0.05, respectively). Nevertheless, the model including birth weight became preferable according to the Akaike's information criterion. Survivability dramatically decreased with extreme birth weights, although it reached a survival probability greater than 93.5% within the 3.3 to 5.4 kg range, indicating an optimum birth weight range of Ripollesa lambs for survival purposes. The hazard ratio (HR) increased for births occurring within the last third of the lambing period (HR = 1.70; P c 0.05), as well as for primiparous ewes that lambed in December and January (HR = 5.36; P < 0.001). Survival probability decreased for lambs born from litters with 1 or more stillbirths or mummified fetuses (HR = 1.61; P < 0.05). The variance component estimated for sire variance (0.07) was clearly lower than that of the common environment (1.87), with a heritability estimate of 0.027. [ABSTRACT FROM AUTHOR]

Given that correct assumptions on the baseline survival function are determinant for the validity of further inferences, specific tools to test the fit of a model to real data become essential in proportional hazards models. In this sense, we have proposed a parametric bootstrap to test the fit of survival models. Monte Carlo simulations are used to generate new data sets from the estimates obtained through the assumed models, and then bootstrap intervals can be established for the survival function along the time space studied. Significant fitting deficiencies are revealed when the real survival function is not included within the bootstrap interval. We tested this procedure in a survival data set of Bruna dels Pirineus beef calves, assuming 4 parametric models (exponential, Weibull, exponential time-dependent, Weibull time-dependent) and the Cox's semiparametric model. Fitting deficiencies were not observed for the Cox's model and the exponential time-dependent model, whereas the Weibull time-dependent model suffered from moderate overestimation at different ages. Thus, the exponential time-dependent model appears to be preferable because of its correct fit for survival data of beef calves and its smaller computational and time requirements. Exponential and Weibull models were completely rejected due to the continuous over- and underestimation of the survival probability reported. Results here highlighted the flexibility of parametric models with time-dependent effects, achieving a fit comparable to nonparametric models. [ABSTRACT FROM AUTHOR]

The Bayes factor (BF) procedure was applied to examine the additive genetic component of several physiological and vitality variables for newborn pigs. Nine variables were studied: heart rate, arterial oxygen saturation, rectal temperature (all at birth and 60 mm later), birth weight, interval between birth and first teats contact, and interval between birth and first colostrum intake. The available numbers of data ranged from 288 (heart rate at 60 mm) to 839 records (birth weight) from F2 Iberian × Meishan newborn pigs. We compared a model with zero heritability (nonheritable) with the one where the additive genetic background was included. The BF was used to discriminate between both candidate models. Very strong evidence of genetic background was detected for heart rate 60 mm after birth (BF = 48.90), and strong evidence was detected for rectal temperature at birth (BF = 13.82). Posterior modes (means) of heritabilities were 0.29 (0.32) and 0.40 (0.39), respectively. In addition, substantial evidence of absence of genetic background was detected for arterial oxygen saturation at birth. [ABSTRACT FROM AUTHOR]

Iberian × Meishan F2 piglet's preweaning survivability was analyzed using categorical data regression procedures within the proportional hazards assumption. A frailty sire model was assumed with the litter effect treated as an additional random source of variation. Moreover, the relative birth weight within litter and the litter effect were considered time-dependent covariates that changed their values in the second day of life due to cross fostering carried out to standardize litters. Six variables had a significant effect on survivability: birth weight (P < 0.01), relative birth weight within litter (P < 0.001), rectal temperature 60 min after birth (P < 0.01), type of presentation at birth (P < 0.05), presence of stillbirths (P < 0.001), and presence of mummified fetuses (P < 0.001). Small piglets (<0.98 kg) suffered a high hazard ratio (6.57; P < 0.001), with this variable being clearly lower for the rest of birth weight categories. Piglets that were small in relation to their siblings (relative birth weight within litter) also suffered an increased death risk, with a hazard ratio of 1.81 (P < 0.05), which was similar to animals with posterior presentations at birth (hazard ratio = 1.80; P < 0.05). Piglets with a rectal temperature lower than 35.4°C 60 min after birth showed the highest hazard ratio (7.18; P < 0.01). Furthermore, the presence of mummified fetuses decreased the survivability of the remaining siblings, with a hazard ratio of 2.03 (P < 0.01), as did the presence of stillbirths (hazard ratio = 3.55; P < 0.001). The inclusion of the two random effects allowed us to estimate the mode of the joint posterior density of the sire variance (0.08) and the litter variance (1.98). The estimated heritability of preweaning survival reached a value of 0.03. We conclude that piglet survival involves several systematic influences related to birth weight, thermoregulatory ability, and injuries suffered during gestation and farrowing. The genetic variance was small compared with those generated by the common environment, for which the genetic improvement of piglet survival seems difficult. [ABSTRACT FROM AUTHOR]

Eight physiological and vitality variables related to the first minutes of the life of Iberian × Meishan F2 piglets were analyzed. Variables included rectal temperature, arterial oxygen saturation, heart rate (all were monitored at birth and 60 min later), time to reach the udder, and time to the first colostrum intake. Litter effect within boar, a random source of variation, influenced all variables, whereas the random boar effect was not significant for heart rate and arterial oxygen saturation at birth and the time to reach the udder. Birth weight influenced rectal temperature at birth (P < 0.01) and 60 min later (P < 0.001), arterial oxygen saturation i h after birth (P < 0.05), heart rate 1 h after birth (P < 0.001), time to reach the udder (P < 0.05), and time to the first suckle (P < 0.01), whereas viability score influenced rectal temperatures (at birth and 60 min later; P < 0.001), and time to reach the udder (P < 0.001) and to suckle (P < 0.001). Finally, the order of birth showed significant effects for rectal temperature I h after birth (P < 0.001) and time to first colostrum intake (P < 0.001). Correlation coefficients between physiological variables were generally low, with the exception of the ones for rectal temperature 1 h after birth and arterial oxygen saturation i h after birth (0.38; P < 0.001), and for rectal temperature at birth and rectal temperature 1 h after birth (0.34; P < 0.001). Times to reach the udder and to suckle were highly correlated (0.67; P < 0.001) and also were moderately and negatively correlated with rectal temperature i h after birth (-0.36 and -0.38 respectively; P < 0.001). Heart rates at birth and 1 h later, as well as arterial oxygen saturation at birth and I h later, were not correlated, showing that the values at birth do not necessarily provide information about the physiological status of the pig 1 h after birth. The recording of physiological variables, birth weight, and viability score may be useful to identify weak piglets quickly and to establish palliative measures. [ABSTRACT FROM AUTHOR]