Prolificacy and Its Relationship with Age, Body Weight, Parity, Previous Litter Size and Body Linear Type Traits in Meat-type Goats.

Data on age and body weight at breeding, parity, previous litter size, days open and some descriptive body linear traits from 389 meat-type, prolific Black Bengal goats in Tripura State of India, were collected for 3 and 1/2 years (2007 to 2010) and analyzed using logistic regression model. The objectives of the study were i) to evaluate the effect of age and body weight at breeding, parity, previous litter size and days open on litter size of does; and ii) to investigate if body linear type traits influenced litter size in meat-type, prolific goats. The incidence of 68.39% multiple births with a prolificacy rate of 175.07% was recorded. Higher age (>2.69 year), higher parity order (>2.31), more body weight at breeding (>20.5 kg) and larger previous litter size (>1.65) showed an increase likelihood of multiple litter size when compared to single litter size. There was a strong, positive relationship between litter size and various body linear type traits like neck length (>22.78 cm), body length (>54.86 cm), withers height (>48.85 cm), croup height (>50.67 cm), distance between tuber coxae bones (>11.38 cm) and distance between tuber ischii bones (>4.56 cm) for discriminating the goats bearing multiple fetuses from those bearing a single fetus.

INTRODUCTION

Linear traits are currently used for classifying goat breeds (Yakubu et al., 2010) and thus suggested for the inclusion of linear dimensions in breeding programmes of goats (Wiggans and Hubbard, 2001). Where genetic evaluation has still limited use in developing countries, identification of some descriptive linear traits may be useful and farmers’ friendly tools for selecting goats with desirable characters.

Increasing kid production is of great interest to goat producers and thus, prolificacy is an economically important trait. Despite the wide use of evaluation for type in goats, little is known about the association between body conformation and some reproductive traits in this species. Linear type traits like physical strength, less sloped rump angle show positive relationships with litter size in dairy goats (Mellado et al., 2008). Although ultrasonographic scanning is used for the prediction of fetal number in goats (Abdelgafar et al., 2007), there is a dearth of knowledge on phenotypic descriptors to differentiate the pregnant goats bearing single or multiple fetuses under field conditions. Determination of fetal number during pregnancy in prolific goat breed is always helpful for differentiating the pregnant goats bearing single or multiple fetuses that aid in care and management to prevent pregnancy toxemia, dystocia and ensure the optimal birth weights and survivability of kids born to goats bearing multiple fetuses (Karen et al., 2006; Moallem et al., 2012). To date, there is an inadequate knowledge on the association between descriptive linear traits and the pregnant goats bearing single or multiple fetuses. A comprehensive description of phenotypic traits for differentiating the goats bearing single and multiple fetuses is thus urgently needed. Black Bengal goat is a small-sized, meat type, early maturing, non-seasonal and prolific breed (Acharya, 1982), in eastern and north eastern regions of India (Zeshmarani et al., 2007), Bangladesh (Rahman et al., 2006) and also in Pakistan (Khanum et al., 2007). This goat breed is famous for its quality of meat and leather with high fecundity. We thus aimed to evaluate the effect of age and body weight at breeding, parity, previous litter size and days open on litter size of does; and to test the hypothesis that body linear type traits might be associated to litter size in prolific goats.

MATERIALS AND METHODS

Study area

The study area included 23 villages in 3 districts (West Tripura, South Tripura and Dhalai) of Tripura state (longitude: 91°30′ E; latitude: 23°45′ N) in the north-eastern region of India. Two-stage stratified random sample survey (consisting of village as strata- 1 and animal within village as strata- 2) was conducted on 389 pregnant Black Bengal goats for the collection of monthly phenotypic data during 5 months period of pregnancy between June 2007 to December 2010 for 3 and 1/2 years. Care was taken to select villages located in isolation to minimize genetic linkage among the animals. The average annual rainfall for these areas is 2,100 mm. The climate is hot and humid with temperature ranging from 10°C in winter to 35°C in summer.

Experimental animals

In the villages, goats are reared in small groups of 2 to 5 goats per household. The animals are allowed for free grazing during the day on natural pasture, tree lopping and scavenging on kitchen waste when available. The female goats are mated naturally by the bucks available in the village. The study commenced with the primary visits to identify individual female goats and their owners. The goats having certain prominent physical characteristics like small size, deep body, straight, convex or concave head, short legs, thin and shiny hair coat of any single colour (black, white and brown) or mixed (combination of any two of the pure colours), erect or pendulous ears, curved or straight horn, short tail were identified as Black Bengal goats (Sahoo et al., 2005). The selected animals were free from diseases and had at least one previous kidding but less than 6 kiddings. Goats between their first and second month of pregnancy, as reported by the owners, were selected and identified by a neck tag with certain number. The animal’s feeds were supplemented with standard concentrate mixture at 200 gm per head per day, vitamins and mineral mixture and examined regularly by a veterinarian during the study period.

Data collection

Herd details were recorded on an individual data card consisted of number of females, kids and whether buck maintained in the herd, age of the female goats, parity, last date of kidding, previous litter size, current date/month of last breeding etc. Sixteen phenotypic variables measured included ear length, neck length, head rump length (HRL), head rump curve length (HRCL), body length (BL), sacral length (SL), body weight (BW), heart girth (HG), punch girth (PG), pelvic triangle area (PTA), withers height (WH), croup height (CH), clearance at sternum (CS), clearance at udder (CU), distance between tuber coxae bones (DTC) and distance between tuber ischii bones (DI). The length, height and circumference measurements (centimeters) were taken using a measuring tape. The measurements (centimeters) of DTC and DI were done using a calibrated wooden caliper. Body weight (kg) was taken using a spring balance. All measurements were collected from the goats at the third week of each month during five months pregnancy period in the morning before the animals were released for grazing. Measurements were recorded on an individual data sheet. Does that aborted during the study period were excluded from the experiment. Finally, 389 does completed the term and gave birth to live kid(s). Each individual doe was sampled once only. The numbers of single, twin and triplet pregnancies were 123, 240 and 26, respectively. Litter size (LZ) was defined as the number of kid(s) born by each goat and considered at 3 levels: single, twin and triplet births. Prolificacy was calculated as the percentage of number of kids born on total number of does delivered according to the following equation.

Statistical analysis

Data were analyzed in two steps. In the first step, age, parity, previous litter size, days open, body coat colour and data of sixteen phenotypic variables covering ear length, neck length, head rump length, head rump curve length, body length, sacral length, body weight, heart girth, punch girth, pelvic triangle area, withers height, croup height, clearance at sternum, clearance at udder, distance between tuber coxae bones and distance between tuber ischii bones collected on 389 does were analyzed using multiple stepwise logistic forward regressions (PROC LOGISTIC; SAS, 2012) to do the stepwise selection for identifying the prognostic factors for litter size. The model used was the logistic regression (LR) model:

where, litter size is the dependent variable; age, parity, previous litter size, days open, ear length, neck length, body weight, head lump length, head rump curve length, body length, sacral length, wither height, croup height, heart girth, punch girth, clearance at udder, clearance at sternum, pelvic triangle area, distance between tuber coxae bones, distance between tuber ischii bones are the predictor variables that ‘x’ represents; α = intercept of the model and β = the coefficient of predictor variables.

LR enters predictor variables in a stepwise manner; it will also fit specified models or perform forward selection of variables. At each step in the stepping process, an attempt was made to remove any insignificant variables from the model before adding a significant variable to the model. A significance level of 0.3 was required to allow a variable to enter into the model, and a significance level of 0.35 was required for a variable to stay in the model. The step selection was based on the maximum likelihood ratio. The evaluation of the logistic regression model was done by Hosmer and Lemeshow Goodness-of-Fit Test in the logistic procedure.

The response variable chooses litter size = 1 (when a pregnant goat carrying more than one kid) as against litter size = 0 (when a pregnant goat carrying one kid only) so that the probability of litter size is modeled. There were some missing data for some individuals during the five month period study. Hence, an average value of monthly data for each linear trait was considered for analysis. From stepwise logistic regressions, some traits were emerged as the prognostic factors that could influence litter size significantly.

In the second step, we wanted to find out certain point of measurement of various identified factors that could indicate the likelihood of multiple fetuses in does during pregnancy. An average value for each identified (stepwise logistic procedure) trait was determined for classifying the does into one of two categories of identified traits: those above or below average value for each trait. The null hypothesis that odds of prolificacy (does carrying two or more fetuses vs. does carrying single fetus) is the same in does classified as above or below average for different identified traits was tested by using logistic regression model (PROC LOGISTIC; SAS, 2012). The identified potential traits associated with litter size were included in the model as independent variables. The resulting values (odds ratios) were then interpreted as decisive factor. An odd ratio of ‘1’ indicates that the factor examined does not alter the chance of single or multiple fetuses. Odd ratios >1 imply a direct association with an increase likelihood of multiple fetuses. Odd ratios <1 imply an inverse association with the likelihood of multiple fetuses.

RESULTS AND DISCUSSION

Prolificacy

The percent of Black Bengal does showing different birth type based on kidding size is presented in Table 1. The litter size proportions for single, twin and triplet were 31.61%, 61.70% and 6.69%, respectively. The present survey recorded the birth of 681 kids from 389 kidding does, averaging 1.75 kids per doe and registering prolificacy rate of 175.07%. The average litter size being 1.75 in this study is quite comparable with some world prolific goat breeds including Nubian, Pygmy, American Alpine, French Alpine, Saanen and Toggenburg with the average litter size of 2.0, 1.9, 1.9, 1.7, 1.7 and 1.6, respectively (Amoah et al., 1996), suggesting that Black Bengal goats are prolific goat breed. Our findings on 68.39% multiple births, of which 61.70% were twinning agree with the previous record of Hassan et al. (2007) who reported the incidence of 56.3% twinning births in Black Bengal goats in neighboring Bangladesh.

Table 1

Indian Black Bengal goats with different litter size groups and prolificacy

	Litter size			Total kids born	Average kid born/doe	Prolificacy (%)
	Single	Twin	Triplet
Number of animal	123	240	26	681	1.75	175.07
% of animal	31.61	61.70	6.69

Factors influencing litter size

Table 2 provides a summary of various descriptive linear traits including age, parity, previous litter size and days open included in the analysis. A limited work on the effect of age, body weight, parity and descriptive linear traits on prolificacy has been reported in small ruminants, especially goats. The result of stepwise logistic regression analysis is presented in Table 3 that shows the factors influencing litter size in pregnant does. Withers height was the best discriminating factor between does carrying multiple fetuses and does bearing single fetus during pregnancy. Distance between tuber coxae bones, neck length, body length, croup height, body weight and age were the next important factors for multiple litter size during pregnancy. Larger litter size was strongly influenced by higher parity and higher previous litter size. The results of the Hosmer and Lemeshow test register the goodness-of-fit of the selected logistic regression model where 0.1216 is a probability value (Table 3). Higher age, body weight, parity and previous litter size have positive influence on multiple births in Black Bengal goats. These findings are in general agreement with the previous reports on does (Constantinou, 1989; Mellado et al., 1991; Amoah et al., 1996) and ewes (Laster et al., 1972; Smith, 1985; Gaskins et al., 2005). Though some genes have been identified for influencing litter size in goats (An et al., 2009; Chu et al., 2011; Feng et al., 2011), certain body weight along with desirable body condition at certain age and parity may be necessary for optimum metabolic requirements which in turn influence hypophyseal-pituitary-gonadal axis for more ovulation that ultimately determine the number of successful fertilization and subsequently litter size (Perry et al., 1991; Tummaruk et al., 2007). In earlier studies, body condition score, higher live weight, age, physical strength and less slopped rump angles were found to be associated with an increase possibility of multiple births in does (Constantinou, 1989; Mellado et al., 2008) and in ewes (Hall et al., 1994; Gonzalez et al., 1997; Gaskins et al., 2005; Aliyari et al., 2012). It may be assumed that the does with higher previous litter size may have potential to give multiple births in subsequent kidding. Goats with longer neck length have showed to give multiple births. Since neck length may also be adjusted with age and body weight, it can’t be ruled out whether neck length is a heritable character for the incidence of multiple births in Black Bengal goat. On the other hand, different linear traits have been suggested to be valuable tools in breeding programmes for some likelihood of productive and reproductive traits in goats (Wiggans and Hubbard, 2001), in sheep (Snowder, 2002) and in dairy cattle (Royal et al., 2002), because of the moderate to high heritability of some linear traits (Luo et al., 1997).

Table 2

Descriptive statistics for different body linear traits including age, parity in meat-type, prolific Black Bengal goats in India

Trait	Mean±SD
Age (yr)	2.69±1.27
Parity (No.)	2.31±1.45
Previous litter size (No.)	1.65±0.67
Days open (d)	104.11±41.54
Ear length (cm)	12.40±1.23
Neck length (cm)	22.78±3.32
Body weight (kg)	20.50±4.68
Head rump length (cm)	57.15±5.16
Head rump curve length (cm)	62.14±4.86
Body length (cm)	54.86±6.60
Sacral length (cm)	12.80±1.26
Withers height (cm)	48.85±3.35
Croup height (cm)	50.67±3.26
Heart girth (cm)	62.70±5.68
Punch girth (cm)	73.14±6.85
Clearance at udder (cm)	22.77±2.68
Clearance at sternum (cm)	24.43±2.38
Pelvic triangle area (square cm)	462.65±138.71
Distance between tuber coxae bones (cm)	11.38±1.18
Distance between tuber ischii bones (cm)	4.56±0.54

Table 3

Summary of stepwise logistic selection of various traits influencing litter size in meat-type, prolific Black Bengal goats in India

Step	Trait entered	DF	Score Chi-square	Pr.> Chi square1
1	Withers height	1	46.0132	<0.0001
2	Distance between tuber coxae bones	1	13.1773	0.0003
3	Body weight	1	2.5079	0.1133
4	Neck length	1	3.6537	0.0559
5	Body length	1	2.5711	0.1088
6	Previous litter size	1	1.3312	0.2586
7	Distance between tuber ischii bones	1	1.368	0.2422
8	Croup height	1	2.3705	0.1236
9	Parity	1	1.608	0.2048
10	Clearance at udder	1	1.9078	0.1672
11	Age	1	2.0224	0.155

The Pr.>Chi square is 0.1216 at 8 d.f. and value of Chi square is 12.7277 as per Hosmer and Lemeshow Goodness-of-Fit test.

Table 4 lists the factors that influenced likelihood of litter size >1. Body weight (>20.5 kg), withers height (>48.85 cm) and parity (>2.31) were associated with four times increased likelihood of multiple litter size. Previous litter size (>1.65), distance between tuber coxae bones (>11.38 cm), neck length (>22.78 cm) registered three times more chances of multiple litter size. Does with >2.69 years age, >54.86 cm body length, >50.67 cm croup height and >4.56 cm distance between tuber ischii bones were two times more likely to have multiple litter size. The odd ratio (<1) for clearance at udder measurement indicated an inverse association of such factor with the likelihood of multiple fetuses. The increased likelihood of bearing multiple fetuses for certain increased measurements in body length, withers height, croup height, distance between tuber coxae bones and distance between tuber ischii bones in pregnant goats could be due to enlargement of the body size during pregnancy to make room for twin or triplet fetuses. Physical strength has been implicated in an increase likelihood of multiple births in does (Mellado et al., 2008) and ewes (Hall et al., 1994). An inverse relationship between clearance at udder measurement and the likelihood of multiple fetuses suggests downward movement of udder in pregnant goats carrying multiple fetuses. Thus, phenotypic variation could be quite attractive for screening the goats bearing single or multiple fetuses.

Table 4

Summary of logistic regression model of the influence of various traits associated with litter size during pregnancy in meat-type, prolific Black Bengal goats in India

Item	N	Litter size (Mean±SD)	Odds ratio estimate	95% wald confidence limits	p-value
Withers height (cm)
<48.85	186	1.55±0.54	1.000	2.532–6.352	<0.0001
>48.85	200	1.92±0.52	4.003
Distance between tuber coxae bones (cm)
<11.38	187	1.57±0.53	1.000	2.126–5.265	<0.0001
>11.38	199	1.90±0.54	3.346
Body weight (cm)
<20.50	192	1.55±0.53	1.000	2.753–7.088	<0.0001
>20.50	194	1.94±0.52	4.417
Neck length (cm)
<22.78	175	1.56±0.53	1.000	2.102–5.149	<0.0001
>22.78	211	1.89±0.54	3.290
Body length (cm)
<54.86	163	1.56±0.51	1.000	1.806–4.361	<0.0001
>54.86	223	1.87±0.56	2.806
Previous litter size (number)
<1.65	158	1.56±0.56	1.000	2.169–5.293	<0.0001
>1.65	228	1.87±0.52	3.389
Distance between tuber ischii bones (cm)
<4.56	179	1.62±0.53	1.000	1.332–3.179	0.0011
>4.56	207	1.85±0.56	2.058
Croup height (cm)
<50.67	192	1.61±0.56	1.000	1.778–4.353	<0.0001
>50.67	194	1.87±0.53	2.782
Parity (number)
<2.31	235	1.60±0.55	1.000	2.466–6.890	<0.0001
>2.31	151	1.95±0.51	4.122
Clearance at udder (cm)
<22.77	184	1.76±0.54	1.000	0.521–1.231	0.3114
>22.77	202	1.72±0.58	0.801
Age (yr)
<2.69	239	1.65±0.57	1.000	1.724–4.560	<0.0001
>2.69	147	1.89±0.51	2.804

IMPLICATIONS

A clear distinction between the pregnant goats bearing single and multiple fetuses is very necessary to the goat keepers for taking proper care of the goats carrying multiple fetuses and achieving economical benefit out of prolificacy trait. The present study generates the basic measurements of descriptive linear traits in meat type, prolific breed like Black Bengal goats. The older does (>2.69 years) with the history of larger previous litter size (>1.65), higher parity (>2.31) and more body weight (>20.5 kg) are associated with an increase likelihood of multiple litter size. This study also indicates a positive relationship between litter size and some body linear type traits. Neck length (>22.78 cm), body length (>54.86 cm), withers height (>48.85 cm), croup height (>50.67 cm), distance between tuber coxae bones (>11.38 cm) and distance between tuber ischii bones (>4.56 cm) show an increased chance of carrying multiple fetuses. Hence, these body linear traits could be used for discriminating the goats bearing multiple fetuses from the goats bearing single fetus.