Prerpints.org logo
Preprint
Article

Risk Factor Analysis and Genetic Parameter Estimation for Pre-weaning Mortality Traits in Boer, Spanish, and Crossbred Goat Kids

Altmetrics

Downloads

86

Views

50

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

08 March 2024

Posted:

11 March 2024

You are already at the latest version

Alerts
Abstract
The objectives of this study were to evaluate fixed risk factors associated with PWM and to estimate genetic parameters for PWM. A total of 927 birth records from a mixed population of purebred and crossbred Boer and Spanish goats born between 2016 and 2023 at the International Goat Research Center (IGRC) were used for this study. Four binary traits were studied: D0-3 (death within 3 days after birth), D4-60 (death between 4 and 60 days), D61-90 (death between 61 and 90 days), and D0-90 (death within 90 days). Logistic regression models were used to evaluate the risk factors associated with PWM traits. Bayesian threshold models and Gibbs sampling were used to estimate the genetic parameters. Birth weight, season, litter size, sex, dam age, breed, and heterosis were found to be significantly associated with at least one of the PWM traits. Heritability estimates were 0.263, 0.124, 0.080, and 0.207, for D0-3, D4-60, D61-90, and D0-90, respectively. The genetic correlations between the studied traits ranged from 0.892 (D0-3 and D0-90) to 0.999 (D0-3 and D61-90). These results suggest that PWM in goats is influenced by both non-genetic and genetic factors and can be reduced by management, genetic selection, and crossbreeding approaches.
Keywords: 
Subject: Biology and Life Sciences  -   Agricultural Science and Agronomy

1. Introduction

Goats are important livestock that produce meat, milk, fiber, and skin for human use. The global goat population has more than doubled in the past four decades, reaching about one billion goats [1]. Goats are hardy and adaptable animals that can survive in diverse climatic conditions and feed on a wide variety of vegetations. They are especially suited for small-scale and resource-poor farmers, both in the US and developing countries. In the US, the goat population rose by over 60% from 1.65 million in 1997 to 2.64 million in 2017, but then slightly decreased to 2.47 million in 2024 [2]. The US goat industry consists of three main segments: meat, dairy, and fiber, with the respective numbers of head in 2024 being 1.95 million, 415 thousand, and 100 thousand [2]. Despite the recent decline in US domestic goat production, the market demand for goat products remains high, driven by factors such as ethnic diversity and health awareness of consumers.
One of the major challenges that goat producers face is pre-weaning mortality (PWM), which is a significant source of economic loss and reduced productivity in goat farming. The PWM in goats varies depending on the breed, environment, and management of the animals and ranges between 7.8% and 46.0% [3,4]. The PWM is influenced by both non-genetic factors, such as birth weight, birth type, birth season, sex, parity, and age of dam [3,5], and genetic factors, such as breed, individual, maternal, and heterosis effects [6,7].
The estimation of genetic parameters for PWM traits is essential for understanding the genetic basis of PWM and designing effective breeding programs to reduce it. However, the literature on the genetic parameters of PWM in goats is scarce, and few studies have considered the effects of breed and heterosis on PWM [6,7]. Breed and heterosis effects are important to consider because they reflect the genetic diversity and adaptation of different goat populations to various environments and production systems [6]. Moreover, crossbreeding is a common practice in goat production to exploit the complementarity and hybrid vigor of different breeds.
Boer and Spanish goats are two of the most widely raised meat goat breeds in US, with different characteristics and performance. Boer goats originated from South Africa and are recognized for their high growth rate and carcass quality [8]. It has become the most popular meat goat in the US since its introduction in the 1990s. The Spanish goats are a landrace breed descended from the goats brought by Spanish explorers in the 1500s [9]. The Spanish goats are valued for their reproduction performance, maternal ability, and disease resistance [10,11]. In the US, there is a widespread trend to cross Boer with Spanish goats to combine the advantages of both breeds and increase the productivity and profitability of meat goat farming. However, there is limited information on the factors affecting PWM. The objectives of this study were to 1) evaluate fixed risk factors associated with PWM; and 2) estimate genetic parameters for PWM traits in a mixed population of purebred and crossbred Boer and Spanish goats. The evaluation of these factors and parameters can provide useful insights for goat producers and breeders to reduce PWM and improve profitability.

2. Materials and Methods

2.1. Animals and Data

The International Goat Research Center (IGRC) at Prairie View A&M University produced purebred and crossbred kids from Boer and Spanish goats. The crossbreeding program involved mating Boer sires with Spanish does (SB F1), and Spanish sires with Boer does (BS F1). The F1 female were then backcrossed with Boer sires (SBB) or Spanish sires (BSS). From 2016 to 2023, 927 kids were born from 233 does mated to 17 Boer and 16 Spanish sires. The breed compositions of the goat kids were: 265 Boer, 426 Spanish, 109 BS, 79 BS, and 30 BSS, and 18 SBB. Each kid born was assigned a unique identification number. Date of birth, sex, birth type, and body weights were recorded. The kids were separated from their respective dams between 28 and 138 days due to varied reasons, with an average of 80 days. High mortality was observed in two pre-weaning periods: day 0 to day 3, accounting for 55.9% of total pre-weaning loss, and day 61 to day 90, accounting for 15.3% of total pre-weaning loss. Four traits were defined for analysis: D0-3 (death within 3 days after birth), D4-60 (death between 4 and 60 days), D61-90 (death between 61 and 90 days), and D0-90 (death within 90 days). Goat kids that died before the trait start date (e.g., stillbirth) or left the herd before the trait end date were excluded from the analysis. Table 1 depicts the number and mortality of kids for each PWM trait. The risk factors considered were birth year, season (Winter: December – February, Spring: March – May), birth type (Single vs Multiple), sex (Female vs Male), birth weight, dam age, breed composition, and retained heterosis. We used the logarithm of birth year as a risk factor to test whether mortality decreased over time. Retained heterosis for each individual was calculated from the pedigree information as R H = 1 − ∑ i = 1 n P S i × P D i , where P S i and P D i indicate the probability that the sire and dam are from breed i , respectively.

2.2. Risk factor analysis

Logistic regression models with logit function implemented in R [12] were used to assess the effect of the risk factors on the PWM traits. The logistic regression model can be written as:
log ⁡ p 1 − p = b ' x ,
where p is the probability of death, b is a vector of coefficients, and x is a vector of risk factors. Odds ratio was calculated as O R = e x p ⁡ ( b ) and the 95 % confidence interval of the OR is calculated as C I = e x p ⁡ ( b ± 1.96 S E ) where SE is the standard error of b.

2.3. Genetic parameter estimation

We used liability-threshold models for the estimation of variance components and genetic parameters. The model assumes that the observed phenotypes are determined by an underlying continuous liability that follows a normal distribution. The equation that relates the liability to the observation can be written as:
y i = 0 i f   l i < t 1 i f   l i ≥ t ,
where y i is the observed phenotype for animal i , l i is the latent liability score for animal i , and t is the threshold value.
A univariate model was used for estimating the heritability for each of the PWM traits, and a bi-variate model was used for estimating genetic correlations between pairwise PWM traits. The univariate model can be written as:
l = X b + Z a + V m + W p e + e ,
where l is the vector of underlying liabilities; X, Z, V, and W are the design matrix for the fixed effect, additive genetic effect, maternal genetic effect, and maternal permanent environmental effect, respectively; b is the vector of fixed effect as described above; a, m, pe, and e are the vectors of additive genetic effect, maternal genetic effect, maternal permanent environmental effect, and residual effect, respectively. The prior distribution for the additive and maternal genetic effects is a m ~ M V N 0 , A ⨂ σ a 2 σ a m σ a m σ m 2 , where A is the numerator additive genetic relationship matrix calculated from the pedigree, which included 1116 individuals, with 55 sires and 262 dams; σ a 2 and σ m 2 are the additive and maternal genetic variances, respectively, and σ a m is the covariance between additive and maternal genetic effects.
The prior distributions for the maternal permanent environmental and residual effects are p e ~ M V N ( 0 , I σ p e 2 ) and e ~ M V N ( 0 ,   I σ e 2 ) , where I is the identity matrix, and σ p e 2 and σ e 2 are the variances for maternal permanent environmental and residual effects, respectively. The direct heritability is calculated as: h 2 = σ a 2 / σ a 2 + 2 σ a m + σ m 2 + σ p e 2 + σ e 2 .
The bivariate model can be written as:
l 1 l 2 = X 1 X 2 b 1 b 2 + Z 1 Z 2 a 1 a 2 + V 1 V 2 m 1 m 2 + W 1 W 2 p e 1 p e 2 + e 1 e 2 ,
where all the model items are as defined previously but with subscripts 1 and 2 denoting trait 1 and trait 2, respectively. The prior distribution for the additive and maternal genetic effects are: a 1 a 2 ~ M V N 0 , A ⊗ σ a 1 2 σ a 1 a 2 σ a 1 a 2 σ a 2 2 , and m 1 m 2 ~ M V N 0 , A ⊗ σ m 1 2 σ m 1 m 2 σ m 1 m 2 σ m 2 2 . The prior distribution for the maternal permanent environmental and residual effects are p e 1 p e 2 ~ M V N 0 , I σ p e 1 2 0 0 I σ p e 2 2 , and e 1 e 2 ~ M V N 0 , I σ e 1 2 0 0 I σ e 2 2 .
The parameters were estimated using Gibbs sampling algorithms implemented in the THRGIBBS1F90 program [13]. The Gibbs sampling process was run for 500,000 cycles, with a burn-in period of 75,000 cycles. The convergence and the burn-in period of the chain were determined by the Geweke convergence diagnostic statistic [14]. The remaining chain was thinned by keeping every 50th sample. Posterior means and standard deviations of the samples were computed using the POSTGIBBSF90 program [13].

3. Results

3.1. Logistic regression analysis of risk factors

The logistic regression analysis revealed several significant risk factors associated with the four PWM traits in the current evaluated goat population. Table 2 presents the estimates of the effects, OR and their corresponding 95% CI, and P-values for the risk factors. Birth season was significantly (P < 0.05) associated with D0-3, with an estimated effect of 0.60 and an OR of 1.82. This indicates higher mortality during winter compared to spring. Birth weight showed a highly significant (P < 0.001) effect on D0-3, D4-60, and D0-90, and significant (P = 0.008) effect on D61-90, with estimates ranging from -0.33 to -0.44 and ORs ranging from 0.65 to 0.72, suggesting that increased birth weight is linked to reduced mortality risk. Birth type also had a significant (P < 0.05) effect on D0-3, with an estimate of -1.42 and an OR of 0.24, indicating higher mortality in multiple births than single births. Breed had significant (P < 0.05) effect on D4-60 and highly significant (P < 0.001) effect on other traits, with estimates ranging from -0.76 to -2.61 and ORs ranging from 0.07 to 0.47, suggesting reduced mortality in the Spanish breed compared to Boer breed. Retained heterosis was also significant (P = 0.009) for D0-3 and highly significant (P < 0.001) for D0-90, with estimates of -0.86 and -0.93 and ORs of 0.42 and 0.39, respectively, suggesting that heterosis contributed to reduced mortality. Dam age had a significant (P < 0.05) effect on D4-60 and D0-90, with estimates of 0.14 and 0.09 and ORs of 1.15 and 1.10, respectively, implying that older dams had higher mortality rates of their offspring. Sex had a highly significant (P < 0.001) effect on D0-90, with an estimate of -0.69 and an OR of 0.50, indicating lower mortality in females than males. Birth year has no significant (P > 0.05) effect on any of the traits, indicating no temporal trend in mortality in the studied population.

3.2. Estimates of genetic parameters

The posterior means and standard deviations of the variance components and the heritability estimates for the four PWM traits are presented in Table 3. The highest direct heritability was observed for D0-3 (0.263), followed by D0-90 (0.207), D4-60 (0.124), and D61-90 (0.080). The maternal heritability for the PWM traits ranged from 0.055 (D4-60) to 0.121 (D0-3), and the proportion of total variance explained by maternal permanent environmental effects varied from 0.027 (D61-90) to 0.043 (D0-3). The results indicate that there is substantial genetic variation for PWM in goats among individuals, and that genetic selection can improve the survival of the kids. The results also suggest that better maternal ability and care of the does can enhance the survivability of the kids.
The posterior means and standard deviations of the additive genetic correlations and phenotypic correlations among the four PWM traits are depicted in Table 4. The genetic correlations were very high for all pairs of traits, ranging from 0.892 (D0-3 and D0-90) to 0.999 (D0-3 and D61-90). These results indicate that there are polymorphic genes that affect multiple PWM traits and selection for one trait would have a strong effect on the others. However, the genetic correlations were not perfect, implying that some genes were specific to certain PWM traits. The phenotypic correlations were also high, ranging from 0.666 (D0-3 and D0-90) to 0.998 (D61-90 and D0-90). The lower phenotypic correlation between D0-3 and D0-90 and D4-60 and D0-90 indicate that environmental factors affect different PWM traits at different periods.

4. Discussion

4.1. Risk factors associated with pre-weaning mortality

Results from our logistic regression analysis showed that birth weight is consistently associated with PWM across all the traits studied. This finding is consistent with previous studies that reported a negative relationship between birth weight and mortality during the pre-weaning period [3,15,16]. According to [17], a 1-kg increase in kid birth weight reduces the probability of mortality by 32.5%. Luo et al. [4] also found that when the birth weight was 20% lower than the annual average, the PWM increased to 46%. Hailu et al. [18] also reported a maximum survival rate of 74% from kids greater than 3 kg at birth. These results suggest that birth weight is a crucial factor for survival of goat kids. Heavier kids are likely to have more sufficient energy reserves to maintain body heat and stronger immunity to resist diseases and environmental hazards [19,20]. Therefore, appropriate management, including feeding with colostrum, milk replacers, and brooding with artificial heating for kids with light weights is needed to improve their survival.
Birth type had a significant effect on trait D0-3 but not for the other traits from day 4 to day 90, indicating that multi-born kids are more vulnerable in early postnatal period than single-born kids. One likely reason for this is that multi-born kids tend to have lower birth weights than single-born kids, which may impair their ability to cope with environmental hazards or nurse colostrum from their mothers shortly after birth. In this study, the average birth weight was 3.65kg for single-born kids and 3.28 kg for multi-born kids. Moreover, multi-born kids may face competition from their siblings for accessing their mothers’ teats, which may further reduce their colostrum intake. Previous studies have also reported higher mortality rates for multi-born kids than single-born kids in various breeds [3,21,22,23,24]. This study did not separate the effects of twin-born and triplet-born effects due to the relatively low number of records. However, literature studies have reported that triplet-born kids had higher mortality than twin-born kids [22,23,24]. Therefore, it is recommended that multi-born kids, especially those with low birth weights, receive special care and management to improve their survival and growth performance.
Birth season had a significant influence on trait D0-3 and a marginal influence on trait D61-90. However, the direction of the effect was opposite for these two traits. Winter born kids had higher mortality for D0-3 and lower mortality for D61-90, while spring born kids had lower mortality for D0-3 and higher mortality for D61-90. The higher mortality rate in D0-3 for winter born kids may be attributed to the cold weather in winter. According to the historical weather data, the average temperature between December and February from 2016 and 2023 in Texas is 10.8°C (51.4°F). The higher mortality rate in D61-90 for spring born kids may be related to the hotter temperature in summer around weaning period. The average temperature between May and August in Texas is 28.3°C (82.9 °F). Previous literature studies have reported that both cold and hot weather can affect PWM. Luo et al. [4] found that the mortality for the first month after birth increases when temperature is below 11°C (52°F), while the mortality for the month between 2 to 4 after birth increases when temperature is greater than 26°C (79°F). Niverthikaa et al. [25] identified cold weather and heat stress as two of the major factors affecting PWM in non-native goat breeds reared in Ampara District in Sri Lanka. Snyman [22] from South Africa also found that season had a significant effect on kid mortality in Angora goats, with higher mortality rates in winter and summer than in autumn and spring. These results suggest that seasonal variations affect kid mortalities significantly and management practices should be adjusted accordingly to reduce the risks. Hailu et al. [18] also suggested that during spring season, does could have access to sufficient vegetation to produce milk for their young offspring compared to other seasons of the year.
Male kids had higher PWM than female kids, which is consistent with previous studies in different goat breeds [3,18,21,26]. The sex effect was more noticeable in later stages of pre-weaning than early stages. The reason for low survivability in male goat kids is unclear, but it may be associated with reduced thermoregulation and behaviors such as standing, udder seeking and sucking ability [3]. Male-linked mortality has also been reported in other species, such as piglets [27], sheep [28,29], and humans [30]. Further research is required to determine the causal factors affecting male kid PWM.
The PWM rate increased with the age of the dam in general. The higher mortality in older dams could be due to the decline in maternal ability, milk production, and immune function [5,23,31]. Our data showed that the mortality rate for kids born from dams aged 1, 2, 3, 4, 5, and 6+ years were: 28.6%, 10.7%, 23.5%, 29.9%, 29.7%, and 30.2%, respectively. The lowest mortality rate was observed in kids born from dams aged 2 years, which might be the optimal age for reproduction and lactation in goats. Kids born from dams aged 1 year had higher mortality rate than those born from dams aged 2 and 3 years. A possible reason is that these 1 year old does were bred before reached their optimal reproductive performance. It is worth noting that parity information was not recorded in our dataset. Other studies evaluated the effect of parity instead of dam age and consistently found that first-parity does tend to have higher mortality in their offspring [3,21,24]. First parity does may have low milk production and colostrum quality. They may also have less maternal experience and ability to care for their kids, especially if they have multiple births.
The Boer breed had higher mortality rates across all the pre-weaning traits studied, compared to the native Spanish breed. This is consistent with previous studies that showed the Boer breed had poor survival rates compared to other breeds in different countries and regions. For example, Tesema et al. [32] and Tessama et al. [33] reported that Boer breed had higher mortality rates than the local Ethiopian breeds. Harrison [34] showed that Boer and its crosses with Kiko had higher mortality rates compared to purebred Kiko breed in US. These results suggest that the Boer breed is less adapted to the environment and management system of the study areas, despite its higher growth rate. Crossbreeding can improve pre-weaning survival chances by introducing hybrid vigor. However, the effect of crossbreeding may vary depending on the breeds and environments involved. Pérez-Baena et al. [35] observed that Boer-crossed kids had significantly lower mortality than Murciano-Granadina purebred kids in Spain. Harrison [34] showed the crossbred Boer and Kiko had lower mortality than the purebred Boer, but still significantly higher mortality than the purebred Kiko. In this study, the heterosis effect reduced the probability of PWM by up to 61%, suggesting the potential of crossbreeding Boer with the local Spanish breed in US.

4.2. Genetic parameters for pre-weaning mortality traits

Genetic parameter estimates for PWM traits in goats are limited in literature. Previous studies have reported low heritability estimates for this trait, ranging from 0.03 to 0.10 in different goat breeds and populations. For example, Synman [22] estimated a heritability of 0.04 for PWM in South African Angora goats. Rout et al. [36] estimated a heritability of 0.03 in Jamunapari goats in India. Josiane et al. [37] estimated a heritability of 0.04 in indigenous goats in Burundi. Tesema et al. [38] estimated a heritability of 0.10 for kid survival at 3 months of age in a crossbred population of Boer and Central Highland goats in Ethiopia. Our estimates of heritability were higher than those reported in literature. This could be due to the genetic diversity of population we studied, which was consisted of Boer, Spanish, and their crossbreds. Another possible explanation is that most of the previous studies used field data, which may have more measurement errors and unaccounted environmental factors. However, the number of records used in our study was still relatively small, which resulted in large standard deviations of the estimates. Therefore, more data are needed to obtain more accurate estimates of the genetic parameters for PWM in our population.
Heritability estimates for different periods of PWM in goats are lacking in literature. Several studies conducted in sheep have reported low heritability estimates in different pre-weaning periods. Sallam [39] reported that heritability estimates ranged from 0.011 to 0.020 for D0-3, D4-60, D61-90, and D0-90 in Barki lamb in Egypt. Brien et al. [40] found similar results for D0-3 and D4-85 in different sheep breeds in Australia, with heritability estimates of 0.014 and 0.010, respectively. Vostry and Milerski [41] reported heritability estimates for D0-1 (0.024) and D2-14 (0.033) in different sheep breeds in the Czech Republic. Our results suggest that the direct heritability was highest for D0-3 (0.263) and lowest for D61-90 (0.080), which is reasonable as more random environmental effects will likely affect mortality in later stages. The results suggest that selecting for D0-3 could be a viable strategy for improving PWM. One advantage of this strategy is that D0-3 has the highest heritability and therefore more genetic gain can be achieved than selection for other traits. Furthermore, more phenotypic records can be used for D0-3 as some records will be missing for traits that are recorded later as animals leave the herds for varied reasons, such as slaughter or sale.
Our results showed that maternal genetic effects explained less variation than additive genetics, but more than maternal permanent environmental effects. Previous studies have reported different maternal effects for PWM in goats. For example, Rout et al. [36] estimated a maternal permanent environmental effect of 0.01 in Jamunapari goats. Josiane et al. [37] estimated a maternal genetic heritability of 0.04, similar to the direct heritability, in indigenous goats in Burundi. Tesema et al. [38] estimated a maternal permanent environmental effect of 0.06 in Boer-Central Highland crossbred goats. The combined maternal effects in our study, both genetic and environmental, accounted for around 8-16% of the phenotypic variance. This suggests that the does’ selection and management are important for their kids’ survival.
Genetic correlations among different time periods of PWM for goats have not been reported in literature. However, several studies have estimated these correlations for other livestock species. For example, Hansen et al. [42] estimated a genetic correlation of 0.73 between D1-14 and D15-60 in Danish Holstein cattle. Sallam [39] reported moderate to high genetic correlations among PWM traits, ranged from 0.72 (D0-3 and D4-60) to 0.95 (D4-60 and D61-90) in Barki lamb in Egypt. Ibi et al. [43] estimated a moderately high genetic correlation of 0.69 between D15-60 and D61-180, 0.50 between D1-14 and D15-60, but low genetic correlation of 0.06 between D1-14 and D61-180 in Japanese Black beef cattle. The positively high genetic correlations between pairs of the traits in our study suggest that mortality at different pre-weaning periods may share some common genes or genetic factors. However, these correlations may vary depending on the population and environment. Therefore, more studies are needed to estimate the genetic correlations among PWM in different time periods for different goat populations.

5. Conclusions

Pre-weaning mortality in goats is affected by both genetic and non-genetic factors. Our study identified birth season, birth type, sex, birth weight, dam age, breed, and heterosis as significant risk factors for this trait. There are substantial genetic variations for PWM among individuals and maternal effects. To reduce PWM in goats, it is essential to account for all the factors that influence this economically important trait. Better management practices and genetic approaches such as selective breeding and crossbreeding that exploit breed and heterosis effects can be used to improve pre-weaning survival of goat kids. The reason for lower survivability in male goat kids compared to female goat kids is unclear, necessitating further research to pinpoint the primary cause of observed higher mortality rate in male goat kids.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work and approved it for publication. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the USDA National Institute of Food and Agriculture Evans-Allen 1890 Research Formula Program under Section 1445.

Institutional Review Board Statement

The study was conducted in the International Goat Research Center (IGRC) at Prairie View A&M University. All procedures involving animals were reviewed and approved by the Institutional Animal Care and Use Committee of the Prairie View A&M University (IACUC protocol number: AUP 2022-041).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available from the corresponding authors upon reasonable request. The data are not publicly available due to privacy or ethical restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Utaaker, K.S.; Chaudhary, S.; Kifleyohannes, T.; Robertson, L.J. Global goat! Is the expanding goat population an important reservoir of Cryptosporidium? Front. Vet. Sci. 2021, 8, 648500. [Google Scholar] [CrossRef] [PubMed]
  2. USDA-NASS. Sheep and Goat Inventory Report. Available online: https://usda.library.cornell.edu/concern/publications/000000018 (accessed on 28 February 2024).
  3. Chauhan, I.S.; Misra, S.S.; Kumar, A.; Gowane, G.R. Survival analysis of mortality in pre-weaning kids of Sirohi goat. Animal. 2019, 13, 2896–2902. [Google Scholar] [CrossRef]
  4. Luo, N.; Wang, J.; Hu, Y.; Zhao, Z.; Zhao, Y.; Chen, X. Cold and heat climatic variations reduce indigenous goat birth weight and enhance pre-weaning mortality in subtropical monsoon region of China. Trop. Anim. Health Prod. 2020, 52, 1385–1394. [Google Scholar] [CrossRef]
  5. Tesema, Z.; Tilahun, M.; Deribe, B.; Lakew, M.; Belayneh, N.; Zegeye, A.; Aychew, D. Effect of non-genetic factors on pre-weaning growth, survivability and prolificacy of Central Highland x Boer crossbred goats in North Eastern Ethiopia. Livest. Res. Rural Dev. 2017, 29, 11–17. [Google Scholar]
  6. Mavrogenis, A.P.; Antoniades, A.; Koumas, A. Genetic and environmental effects on kid mortality in Cyprus goats. Small Rumin. Res. 2018, 160, 1–5. [Google Scholar]
  7. Snyman, M.A. Genetic parameters for pre-weaning mortality in South African Angora goats. Small Rumin. Res. 2019, 176, 1–6. [Google Scholar]
  8. Casey, N.H.; Van Niekerk, W.A. The Boer goat. I. Origin, adaptability, performance testing, reproduction and milk production. Small Rumin. Res. 1988, 1, 291–302. [Google Scholar] [CrossRef]
  9. Glimp, H.A. Meat goat production and marketing. J. Anim. Sci. 1995, 73, 291–295. [Google Scholar] [CrossRef]
  10. Browning Jr, R.; Leite-Browning, M.L.; Byars Jr, M. Reproductive and health traits among Boer, Kiko, and Spanish meat goat does under humid, subtropical pasture conditions of the southeastern United States. J. Anim. Sci. 2011, 89, 648–660. [Google Scholar] [CrossRef]
  11. Pellerin, A.N.; Browning, R. Comparison of Boer, Kiko, and Spanish meat goat does for stayability and cumulative reproductive output in the humid subtropical southeastern United States. BMC Vet. Res. 2012, 8, 1–9. [Google Scholar] [CrossRef]
  12. R Core Team. R: A language and environment for statistical computing; R Foundation for Statistical Computing: Vienna, Austria, 2021.
  13. Tsuruta, S.; Misztal, I. THRGIBBS1F90 for estimation of variance components with threshold-linear model. In: Proceedings of 8th World Congress on Genetics Applied to Livestock Production, Belo Horizonte, Brazil, 13-18 August 2006.
  14. Geweke, J. Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments. Research Department Staff Report 148. Federal Reserve Bank of Minneapolis, 1991. [Google Scholar]
  15. Husain, S.S.; Horst, P.; Islam, A.B.M.M. Effect of different factors on pre-weaning survivability of Black Bengal kids. Small Rumin. Res. 1995, 18, 1–5. [Google Scholar] [CrossRef]
  16. Deribe, G.; Abebe, G.; Tegegne, A. Non-genetic factors influencing reproductive traits and pre-weaning mortality of lambs and kids under smallholder management, Southern Ethiopia. J. Anim. Plant Sci. 2014, 24, 413–417. [Google Scholar]
  17. Yitagesu, E.; Alemnew, E. Mortality rate of Boer, Central Highland goat and their crosses in Ethiopia: Nonparametric survival analysis and piecewise exponential model. Vet. Sci. 2022, 8, 2183–2193. [Google Scholar] [CrossRef]
  18. Hailu, D.; Mieso, G.; Nigatu, A.; Fufa, D.; Gamada, D. The effect of environmental factors on preweaning survival rate of Borana and Arsi-Bale kids. Small Rumin. Res. 2006, 66, 291–294. [Google Scholar] [CrossRef]
  19. Casellas, J.; Caja, G.; Such, X.; Piedrafita, J. Survival analysis from birth to slaughter of Ripollesa lambs under semi-intensive management. J. Anim. Sci. 2007, 85, 512–517. [Google Scholar] [CrossRef]
  20. Abdelqader, A.; Irshaid, R.; Tabbaa, M.J.; Abuajamieh, M.; Titi, H.; Al-Fataftah, A.R. Factors influencing Awassi lambs survivorship under fields conditions. Livest. Sci. 2017, 199, 1–6. [Google Scholar] [CrossRef]
  21. Singh, M.K.; Rai, B.; Sharma, N. Factors affecting survivability in Jamunapari kids under semi-intensive management system. Indian J. Anim. Sci. 2008, 78, 178–181. [Google Scholar]
  22. Snyman, M.A. Factors affecting preweaning kid mortality in South African Angora goats. S. Afr. J. Anim. Sci. 2010, 40, 54–64. [Google Scholar] [CrossRef]
  23. Subramaniyan, M.; Thanga, T.V.; Subramanian, M.; Senthilnayagam, H. Factors affecting pre-weaning survivability of kids in an organized goat farm. Int. J. Livest. Res. 2016, 6, 83–92. [Google Scholar] [CrossRef]
  24. Singh, M.K.; Pourouchottamane, R.; Singh, S.P.; Kumar, R.A.; Sharma, N.I.; Kumar, A.K.; Dass, G.O.; Pundir, R.K. Non-genetic factors affecting pre-weaning growth and survival rate in Barbari kids under semi-intensive management system. Indian J. Anim. Sci. 2022, 92, 1081–1087. [Google Scholar] [CrossRef]
  25. Niverthika, N.; Fouzi, M.N.M.; Nikzaad, R.M.; Thaiuba, A. Factors affecting pre-weaning mortality in non-native goat breeds reared in Ampara. Int. J. Res. 2023, 133, 155–161. [Google Scholar] [CrossRef]
  26. Barbind, R.P.; B. I. Dandewar. Pre-weaning mortality pattern in osmanabadi crossbred goats. Indian J. Anim. Res. 2004, 38, 75–76. [Google Scholar]
  27. Baxter, E.M.; Jarvis, S.; Palarea-Albaladejo, J.; Edwards, S.A. The weaker sex? The propensity for male-biased piglet mortality. PLoS One, 2012, 7, e30318. [Google Scholar] [CrossRef]
  28. Southey, B.R.; Rodriguez-Zas, S.L.; Leymaster, K.A. Discrete time survival analysis of lamb mortality in a terminal sire composite population. J. Anim. Sci., 2003, 81, 1399–1405. [Google Scholar] [CrossRef]
  29. Sawalha, R.M.; Conington, J.; Brotherstone, S.; Villanueva, B. Analyses of lamb survival of Scottish Blackface sheep. Animal, 2007, 1, 151–157. [Google Scholar] [CrossRef]
  30. Kraemer, S. The fragile male. Bmj, 2000, 321, 1609–1612. [Google Scholar] [CrossRef]
  31. Ruvuga, P.R.; Maleko, D.D. Dairy goats’ management and performance under smallholder farming systems in Eastern Africa: the systematic review and meta-analysis. Trop. Anim. Health Prod. 2023, 55, 255. [Google Scholar] [CrossRef]
  32. Tesema, E.Y.; Alemnew, E. Morbidity and Mortality Rates of Boer and Central Highland Goats at Ataye Crossbreeding Program Research Site: Non-parametric Survival Analysis and Piecewise Exponential Model. 2021. Available online at https://www.researchsquare.com/article/rs-177769/v1 (accessed on February 28, 2024).
  33. Tessema, T.; Teshome, D.; Kumsa, S. Evaluation of growth and survival performances of Pure Borana goat and their crosses with Boer goats. Int. J. Genet. Mol. Biol. 2022, 14, 1–8. [Google Scholar]
  34. Harrison, A. Reproductive and Progeny Performance of Pure and Composite Breeds of Meat Goat Over Two Breeding Cycles. Doctoral dissertation, Tuskegee University, Tuskegee, Alabama, 2023.
  35. Pérez-Baena, I.; Jarque-Durán, M.; Gómez, E.A.; Díaz, J.R.; Peris, C. Terminal crossbreeding of murciano-granadina goats to boer bucks: effects on reproductive performance of goats and growth of kids in artificial rearing. Animals. 2021, 11, 986. [Google Scholar] [CrossRef]
  36. Rout, P.K.; Matika, O.; Kaushik, R.; Dige, M.S.; Dass, G.; Singh, M.K.; Bhusan, S. Genetic analysis of growth parameters and survival potential of Jamunapari goats in semiarid tropics. Small Rumin. Res. 2018, 165, 124–130. [Google Scholar] [CrossRef]
  37. Josiane, M.; Gilbert, H.; Johann, D. Genetic parameters for growth and kid survival of indigenous goat under smallholding system of burundi. Animals. 2020, 10, 135. [Google Scholar] [CrossRef]
  38. Tesema, Z.; Alemayehu, K.; Kebede, D.; Getachew, T.; Deribe, B.; Taye, M.; Tilahun, M.; Kefale, A.; Belayneh, N.; Yizengaw, L. Genetic analysis of survival potential of Boer x Central Highland goats under semi-intensive management. Small Rumin. Res. 2020, 193, 106253. [Google Scholar] [CrossRef]
  39. Sallam, A.M. Risk factors and genetic analysis of pre-weaning mortality in Barki lambs. Livest. Sci., 2019, 230, 103818. [Google Scholar] [CrossRef]
  40. Brien, F.D.; Hebart, M.L.; Smith, D.H.; Edwards, J.H.; Greeff, J.C.; Hart, K.W.; Refshauge, G.; Bird-Gardiner, T.L.; Gaunt, G.; Behrendt, R.; Robertson, M.W. Opportunities for genetic improvement of lamb survival. Anim. Prod. Sci., 2010, 50, 1017–1025. [Google Scholar] [CrossRef]
  41. Vostrý, L.; Milerski, M. Genetic and non-genetic effects influencing lamb survivability in the Czech Republic. Small Rumin. Res., 2013, 113, 47–54. [Google Scholar] [CrossRef]
  42. Hansen, M.; Madsen, P.; Jensen, J.; Pedersen, J.; Christensen, L.G. Genetic parameters of postnatal mortality in Danish Holstein calves. J. Dairy Sci. 2003, 86, 1807–1817. [Google Scholar] [CrossRef]
  43. Ibi, T.; Kahi, A.K.; Hirooka, H. Genetic parameters of postnatal mortality and birth weight in J apanese B lack calves. Anim. Sci. J. 2015, 86, 25–30. [Google Scholar] [CrossRef]
Table 1. Number and mortality of goat kids for the pre-weaning mortality traits1.
Table 1. Number and mortality of goat kids for the pre-weaning mortality traits1.
Trait Kids alive (N) Kids died (N) Total kids (N) Mortality (%)
D0-3 828 99 927 10.68
D4-60 691 51 742 6.87
D61-90 556 27 583 4.63
D0-90 556 177 733 24.15
1 D0-3: death within 3 days after birth; D4-60: death between 4 and 60 days after birth; D61-90: death between 61 and 90 days after birth; D0-90: death within 90 days after birth.
Table 2. Estimates of the effect (b ± SE), odds ratio (OR), 95% confidence interval (CI), and P-values for the risk factors associated with the pre-weaning mortality traits1.
Table 2. Estimates of the effect (b ± SE), odds ratio (OR), 95% confidence interval (CI), and P-values for the risk factors associated with the pre-weaning mortality traits1.
Trait Risk factor b ± SE OR (95% CI) P-value
D0-3 Birth year 0.14 ± 0.22 1.15 (0.75, 1.80) 0.535
Season (Winter) 0.60 ± 0.24 1.82 (1.13, 2.93) 0.013
Sex (Female) -0.23 ± 0.23 0.79 (0.51, 1.24) 0.306
Birth weight -0.33 ± 0.08 0.72 (0.61, 0.85) <0.001
Birth type (Single) -1.42 ± 0.61 0.24 (0.06, 0.68) 0.019
Dam age 0.06 ± 0.05 1.06 (0.96, 1.17) 0.253
Breed (Spanish) -1.53 ± 0.27 0.22 (0.13,0.37) <0.001
Heterosis -0.86 ± 0.33 0.42 (0.22,0.79) 0.009
D4-60 Birth year -0.05 ± 0.27 0.96 (0.57, 1.67) 0.869
Season (Winter) 0.05 ± 0.34 1.05 (0.54, 2.02) 0.874
Sex (Female) -0.52 ± 0.31 0.59 (0.32, 1.09) 0.094
Birth weight -0.39 ± 0.12 0.68 (0.54, 0.86) 0.001
Birth type (Single) -0.39 ± 0.55 0.68 (0.20, 1.77) 0.476
Dam age 0.14 ± 0.06 1.15 (1.02, 1.30) 0.019
Breed (Spanish) -0.76 ± 0.36 0.47 (0.23, 0.96) 0.035
Heterosis -0.45 ± 0.43 0.64 (0.26, 1.44) 0.297
D61-90 Birth year 0.15 ± 0.43 1.17 (0.52, 2.90) 0.724
Season (Winter) -1.12 ± 0.58 0.33 (0.09, 0.94) 0.055
Sex (Female) -0.82 ± 0.46 0.44 (0.18, 1.07) 0.072
Birth weight -0.39 ± 0.15 0.67 (0.50, 0.90) 0.008
Birth type (Single) -0.01 ± 0.67 0.99 (0.22, 3.26) 0.990
Dam age 0.00 ± 0.11 1.00 (0.79, 1.24) 0.984
Breed (Spanish) -2.61 ± 0.57 0.07 (0.02, 0.21) <0.001
Heterosis -0.77 ± 0.69 0.46 (0.10, 1.64) 0.261
D0-90 Birth year 0.00 ± 0.19 1.00 (0.69, 1.46) 0.993
Season (Winter) 0.08 ± 0.21 1.09 (0.72, 1.64) 0.686
Sex (Female) -0.69 ± 0.20 0.50 (0.34, 0.74) <0.001
Birth weight -0.44 ± 0.08 0.65 (0.56, 0.75) <0.001
Birth type (Single) -0.70 ± 0.37 0.50 (0.23, 0.99) 0.060
Dam age 0.09 ± 0.04 1.10 (1.01, 1.19) 0.030
Breed (Spanish) -1.67 ± 0.24 0.19 (0.12, 0.30) <0.001
Heterosis -0.93 ± 0.27 0.39 (0.23, 0.66) <0.001
1 D0-3: death within 3 days after birth; D4-60: death between 4 and 60 days after birth; D61-90: death between 61 and 90 days after birth; D0-90: death within 90 days after birth.
Table 3. Posterior estimates of means (and standard deviations) for variance components and heritability values for the pre-weaning mortality traits1,2.
Table 3. Posterior estimates of means (and standard deviations) for variance components and heritability values for the pre-weaning mortality traits1,2.
Trait σ a 2 σ m 2 σ m p e 2 σ e 2 h d 2 h m 2 c 2
D0-3 0.028 (0.023) 0.012 (0.006) 0.004 (0.004) 0.058 (0.013) 0.263 (0.190) 0.121 (0.059) 0.043 (0.037)
D4-60 0.009 (0.007) 0.004 (0.002) 0.002 (0.001) 0.053 (0.005) 0.124 (0.098) 0.055 (0.030) 0.027 (0.021)
D61-90 0.004 (0.003) 0.003 (0.002) 0.001 (0.001) 0.036 (0.003) 0.080 (0.071) 0.073 (0.036) 0.027 (0.024)
D0-90 0.038 (0.024) 0.012 (0.007) 0.006 (0.005) 0.121 (0.016) 0.207 (0.119) 0.068 (0.038) 0.034 (0.029)
1  σ a 2 : additive genetic variance; σ m 2 : maternal genetic variance; σ m p e 2 : maternal permanent environmental variance; σ e 2 : residual variance; h d 2 : direct heritability; h m 2 : maternal genetic heritability; c 2 : maternal environmental variance ratio. 2 D0-3: death within 3 days after birth; D4-60: death between 4 and 60 days after birth; D61-90: death between 61 and 90 days after birth; D0-90: death within 90 days after birth.
Table 4. Posterior estimates of means (and standard deviations) for genetic correlations (above diagonal) and phenotypic correlations (below diagonal) among the pre-weaning mortality traits1.
Table 4. Posterior estimates of means (and standard deviations) for genetic correlations (above diagonal) and phenotypic correlations (below diagonal) among the pre-weaning mortality traits1.
Trait D0-3 D4-60 D61-90 D0-90
D0-3 0.998 (0.003) 0.999 (0.002) 0.892 (0.123)
D4-60 0.996 (0.004) 0.998 (0.004) 0.941 (0.081)
D61-90 0.996 (0.002) 0.996 (0.003) 0.997 (0.010)
D0-90 0.666 (0.029) 0.843 (0.020) 0.998 (0.005)
1 D0-3: death within 3 days after birth; D4-60: death between 4 and 60 days after birth; D61-90: death between 61 and 90 days after birth; D0-90: death within 90 days after birth.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated