Phenotypic Characterization of Native Bovine Population in

Applied Science Reports
www.pscipub.com/ASR
E-ISSN: 2310-9440 / P-ISSN: 2311-0139
DOI: 10.15192/PSCP.ASR.2015.10.1.3844
App. Sci. Report.
10 (1), 2015: 38-44
© PSCI Publications
Phenotypic Characterization of Native Bovine Population in
Northern Tunisia
R. Baccouche1, B. Jemmali2, M. Haddad3, N. M’Hamdi1, M Ben Hamouda4
1. Institut National Agronomique de Tunisie
2. Laboratoire ADIPARA, Ecole Supérieure d'Agriculture de Mateur, Tunisie
3. Fondation Nationale d’Amélioration de la race Chevaline, Sidi Thabet, Tunisie
4. Institut National de la Recherche Agronomique de Tunisie
Corresponding author: R. Baccouche
Paper Information
ABSTRACT
This study aims to describe the native bovine population in Tunisia. This
Received: 17 December, 2014
description is based morphological traits. An inventory was conducted in
three governorates of North Tunisia (Bizerte, Béja and Manouba). A total
Accepted: 19 February, 2015
of 220 animals were described. Animals were individually described by 26
discontinuous descriptors. The results showed that the local cattle
Published: 20 April, 2015
population is formed of two types; la Blonde du Cap Bon and la Brune de
l’Atlas. It turns out that these animals have many similarities. The analysis
of variance (ANOVA) revealed eight significant variables with location
Citation
effect and seven with the gender effect. Multiple factor analysis (MFA)
generated two distinct groups of animals that differ only in the characters
Baccouche R, Jemmali B, Haddad M, Hamdi MN,
Hamouda MB. 2015. Phenotypic Characterization of
related to the horns. Gender and geographical location have not
Native Bovine Population in Northern Tunisia. Applied
contributed in this distribution. The discriminant factorial analysis (AFD),
Science Reports, 10(1), 38-44. Retrieved from
using the locality as discriminatory factor, did not allow to group
www.pscipub.com
individuals according to the region because they form a single group. The
(DOI:10.15192/PSCP.ASR.2015.10.1.3844)
hierarchical cluster analysis (HCA) has revealed three classes. The only
variables that were the cause of this dissimilarity are the profile of the back
and the orientation of the horns. It shows in the last analysis that the local
cattle population is diversified. This variability is neither gender nor the
geographic location implying that the phenotypic diversity within the local
cattle population is an individual variability.
© 2015 PSCI Publisher All rights reserved.
Key words: native bovine, population, phenotypic diversity, qualitative characters, ANOVA, MFA, HCA.
Introduction
Native cattle play social and cultural roles and are ranked as the most valued animal in Tunisian farming system.
It is adaptable and able to survive and perform under poor conditions. So, Tunisian cattle resulting from brown twig, type
brown of Atlas. In the early 18th century, they have faced the piroplasmosis bacilliform disease which led the colonial
authorities to import zebu breeds of India « Sindhi » and « Brahamas » and cross them to make the disease disappear. The
results of crosses were satisfying. Although this population is hardy and adapt easily to climate changes, its productivity
remain limited. This prompted the authorities after the independence to introduce specialized breeds to ensure selfsufficiency in milk and thus satisfy the consumer who has had changes in his eating habits. These crosses have affected
these cattle’s genetic structure, which is very confusing nowadays (Djemali and Berger, 1992). The downside these crosses
is that these specialized breeds have adjustment problems that will worsen as a result of climate changes scenarios (FAO,
2008). It is in this perspective that the need to conserve this population is introduced. Several studies indicate that the local
cattle population is formed of two races namely; the Blonde of Cap Bon, and the Brown of the Atlas which divides into
two branches: the gray branch and the tawny branch. The Blonde of Cap Bon is characterized by a white coat, the tawny
branch of the brown of the atlas by a tawny color, while the gray branch is characterized by a gray or black color. This
study aims to identify individuals in the population by phenotypic characterization and a short description of the breeding
system. A survey was conducted in three governorates (Bizerte, Béja and Manouba).
Materials and methods
The study was conducted in three districts namely; Bizerte, Béja, and Mannouba (Figure 1). A survey was
undertaken in 2013 and 2014. Collected information included 220 cows selected randomly for this study from the three
studied sites (Table 1). Structured questionnaires were used to collect information on the location of the farm, the breeds,
uses and preferences of breeds, adaptive, breeding and management.
App. Sci. Report. 10 (1), 2015: 38-44
Table 1. The collected sample’s number at different sites
District
Locality
Bizerte
Séjnène
Béja
Manouba
Total
Nevza
Amdoun
Location
Oued ezzen
Louka 2
Ouled Hmed
Araguib
Cap Sirat
M’chergua
Wled Ghanem
Gafaya
Number of animals
15
17
23
32
76
35
11
11
Total
163
46
11
220
Figure 1. Map of Tunisia indicating the study zones
Phenotypic characterization of animals
Photo shoots of each animal were established. Animals were individually described by the discontinuous
descriptors defined by the FAO (2012). There are 26 descriptors; Body hair, coat, color pattern, body hair coat color, body
skin color, hair type, muzzle color, eyelid color, hoof color, horn presence, horn attachment, horn shape, horn orientation,
horn color, ear shape, ear orientation, hump presence, hump size, hump position, hump shape, backline profile, dewlap
size, rump profile, navel flap, tail attachment, tail length, facial profile and penis sheath for males.
Statistical analysis
All data were processed using Microsoft Excel 2007 Microsoft Office. Many statistical analyzes were used to
describe and classes animals and identify the most significant characters in the differentiation. Indeed, a descriptive
statistical analysis was carried out at the beginning followed by analysis of variance (ANOVA), multiple factor analysis
(MFA), discriminant factorial analysis (AFD) and hierarchical cluster analysis (HCA). The XLSTAT software was used in
this study (Addinsoft, Paris, France, 2007).
The model used for ANOVA was: Y= μ + X+ ε
Where:
Y= the studied traits
µ = means of trait
X= the fixed effect
ε = residual error
Results and discussions
Descriptive analysis
Information on a total of 220 native cows were collected of which 198 females (90%) and 22 males (10%). The
average age of the animals is 60 ± 42 months. Table 2 shows the different qualitative traits studied with their modalities,
numbers and frequencies. We noticed a phenotypic similarity among animals of local population. Since for most studied
traits, there is always a modality that has a very high frequency compared to the others. Four coat colors are found; white,
gray, tawny and black. This confirms the existence of two ecotypes: la Blonde du Cap Bon and la brune de l’Atlas. Our
finding was in agreement with those of many studies in Tunisia (Bel hadj, 1972; Atti et al., 1987 and Ben Belgecem,
2005). The brune de l’Atlas is the most encountered (98.12%) population, which is mainly composed by tawny branch
subjects. Only four individuals are belongs to the Blonde du Cap Bon. 22.27% of the population studied had a medium
size hump with thoracic position and an erected orientation. 9% of individuals have a large navel. These two traits are
similar to those of Zebu of India.
39
App. Sci. Report. 10 (1), 2015: 38-44
Table 2. Morphological characterization of native cattle
Traits
Body hair coat color pattern
Color body hair coat
Body skin color
Hair type
Muzzle color
Eyelid color
Hoof color
Horn presence
Horn attachment
Horn shape
Horn orientation
Horn color
Ear shape
Ear orientation
Hump presence
Backline profile
Dewlap size
Rump profile
Navel flap
Tail attachment
Tail length
Facial profile
Modalities
pie
spotted
uniform
white
tawny
grey
black
Not pigmented
Pigmented
glossy/curly
glossy/straight
dull/curly
dull/straight
Not pigmented
pigmented
not pigmented
pigmented
not pigmented
pigmented
absent
present
One horn
Two horn
straight
curved
Lyre shape
lateral
backward
forward
upward
white
fawn
black
stamped
straight
rounded
lateral
drooping
erect
absent
present
straight
slopes up towards the rump
slopes down from withers
dipped
large
medium
small
sloping
roofy
flat
large
medium
small
low
high
short
long
medium
concave
convex
straight
Observations
15
4
201
4
163
26
27
2
2
2
69
10
139
164
56
48
172
14
206
5
215
5
210
44
155
16
17
5
33
160
35
9
146
25
50
170
205
7
8
171
49
106
4
63
47
40
117
63
189
17
14
20
33
167
114
106
39
92
89
132
17
71
Fréquency (%)
6,82
1,82
91,36
1,82
74,09
11,82
12,27
90,91
9,09
0,91
31,36
4,55
63,18
74,55
25,46
21,82
78,18
6,36
93,64
2,27
97,73
2,33
97,67
20,47
72,09
7,44
7,9
2,33
15,34
74,41
16,28
4,19
67.90
11,63
22,73
77,27
93,18
3,18
3,64
77,72
22,27
48,18
1,82
28,64
21,36
18,18
53,18
28,64
85,91
7,73
6,36
9,09
15
75,91
51,82
48,18
17,73
41,82
40,46
60
7,73
32,27
Analysis of variance
The probabilities associated with location effect are significant for some variables at α level of 5%, as shown in
Table 3. There are eight variables: Body hair coat color pattern, Body skin color, Hair type, Muzzle color, eyelid color,
horn orientation and color, hump presence, navel flap, and facial profile. The probabilities associated with the sex effect
are significant (P< 0.05) for seven descriptors (Table 2): presence, attachment, orientation, and color of horn, navel flap
and tail length.
40
App. Sci. Report. 10 (1), 2015: 38-44
Table 3. Results of one factor analysis of variance to a factor on morphological descriptors.
Location (ddl=7)
R²
F
0,072
2,347
0,061
1,977
0,094
3,144
0,186
6,937
0,097
3,251
0,107
3,613
0,052
1,663
0,026
0,816
0,028
0,885
0,058
1,848
0,083
2,735
0,082
2,699
0,048
1,540
0,034
1,060
0,111
3,772
0,037
1,164
0,032
1,007
0,012
0,364
0,088
2,923
0,059
1,899
0,038
1,187
0,092
3,083
Traits
Body hair coat color pattern
Color body hair coat
Body skin color
Hair type
Muzzle color
Eyelid color
Hoof color
Horn presence
Horn attachment
Horn shape
Horn orientation
Horn color
Ear shape
Ear orientation
Hump presence
Backline profile
Dewlap size
Rump Profile
Navel flap
Tail attachment
Tail length
Facial profile
Pr > F
0,025
0,060
0,004
< 0,0001
0,003
0,001
0,120
0,575
0,519
0,080
0,010
0,011
0,155
0,391
0,001
0,325
0,428
0,922
0,006
0,071
0,312
0,004
Gender (ddl=1)
R²
F
0,002
0,368
0,007
1,440
0,003
0,607
0,017
3,756
0,003
0,677
0,001
0,188
0,001
0,303
0,023
5,190
0,020
4,495
0,000
0,099
0,023
5,095
0,028
6,227
0,005
1,146
0,019
4,198
0,002
0,350
0,006
1,343
0,001
0,317
0,000
0,043
0,268
79,768
0,002
0,514
0,023
5,078
0,000
0,031
Pr > F
0,545
0,231
0,437
0,054
0,411
0,665
0,583
0,024
0,035
0,753
0,025
0,013
0,286
0,042
0,555
0,248
0,574
0,836
< 0,0001
0,474
0,025
0,860
R2: Total variance explained by Imada or sex effect; F: Fischer value; Pr: Probability
Multiple factorial analysis (MFA)
Multiple factorial analysis classify studied traits into two groups of variables; color and general conformation.
The results generated by this analysis provide several interpretable factors. Indeed, the eigen value of each variable was
examined which shows its contribution in explaining the distribution of animals. Based on the results presented in Table 4,
the first two axes explained 17.3% of the total inertia. The first axis is related to the overall conformation variables group
with high coordinates and contributions respectively of 0.9 and 62.12. As for the second dimension, it is linked to two
groups of variables namely color and general conformation with a contribution of 52.35 and 47.65.
An RV coefficient of the link between groups of variables was shown in Table 5. It shows the absence of
correlation between gender and locality with the two groups of variables. Figure 2 shows the spatial distribution of
individuals (cases) and morphology’s traits. Animals were clustered into three separate groups. The first opposes five
individuals (B103, B84, B88, B82, B54) to others, since they are distinguished by absence of horns. However, the second
axis contains the majority of animals (215) and it is constituted by other morphology’s traits. This analysis does not
provide enough information about the phenotypic diversity of the local cattle population. This is confirmed by the lower
eigen value of the two explicative axes (value???). A discriminant analysis seems to be very important to explain
phenotypic diversity.
Table 4. Eigen values of the MFA conducted on the morphological characteristics of native cattle, coordinates and group contribution
Eigenvalues
Variability (%)
cumulated
Coordinates
Color
General Conformation
Locality
Sex
Contribution
Color
General Conformation
Factor 1
1,481
9,748
9,748
Factor 2
1,146
7,545
17,293
Factor 3
0,886
5,833
23,126
0,561
0,920
0,042
0,023
0,600
0,546
0,115
0,009
0,611
0,275
0,007
0,000
37,878
62,122
52,350
47,650
69,012
30,988
Table 5. Matrix of the links between groups of variables: Coefficients RV
Group
Color
General Conformation
Locality
Sex
MFA
Color
1
0,185
0,069
0,023
0,822
General Conformation
locality
sex
1
0,044
0,073
0,713
1
0,015
0,075
1
0,059
41
App. Sci. Report. 10 (1), 2015: 38-44
Figure 2. Representation of the groups in the MFA on morphological data of native cattle
Discriminant analysis (AFD)
Figure 3 shows the distribution of individuals according to the locality factor. According to this analysis, the first
two axes have 84% of the total variability. There is an overlap between animals of different localities. The distribution is
very diffuse and animals form superposed groups where the dissimilarity between them is not obvious enough. Moreover,
we noticed that phenotypic diversity of the local population is not dependent on geographic location. This population is
very homogeneous.
Figure 3. Representation of variables and observations by AFD depending on the locality (F1-F2).
Hierarchical Cluster Analysis (HCA)
The goal of a classification is to constitute homogeneous groups. Figure 4 shows the dendrogram realized on the
different descriptors of the morphology. In fact, three separate groups of animals were identified. The first group contains
60 animals of different sex and age and belonging to different localities (Sejnene, Amdoun, Nefza, Gafaya). It is the same
for the second group, except that this one includes more animals of the order of 128. The third class consists of 32 animal
belonging only three locations namely; Sejnene, Amdoun, Nefza. The morphological characteristics that differentiate these
classes are shown in figure 5. In fact the first group is distinguished by the backline profile descriptor. Individuals in this
group have a sloping back up towards the rump or sloping down from withers. However, the second class is characterized
by the orientation of the horns. In fact, all the individuals in this group have horns forwardly directed. The third class is
similar to the first two except that individuals belonging to this class do not have the characteristics associated with the
profile of the backline or in the orientation of the horns.
42
App. Sci. Report. 10 (1), 2015: 38-44
Figure 4. Hierarchical classification of the animal groups in the native cattle population.
A1 : Body hair coat color pattern ;A2 : Body hair coat color ;A3 : Body skin color ; A4 : Hair type ; A5 : Muzzle
color ;A6 : Eyelid color ;A7 : Hoof color ;B1 : Horn presence ;B2 : Horn attachment ;B3 : Horn shape ;B4 : Horn
orientation ;B5 : Horn color ;C1 : Ear shape ;C2 : Ear orientation ;D1 : Hump presence ;D2 : Backline profile ;E :
Dewlap size ;F : Rump profile ;G : Navel flap ;H1 : Tail attachment;H2 : Tail length ;J : Facial profile
Figure 5. The morphological characters of each group of indigenous cattle CAH
Conclusion
The local cattle population is a very original source as it is in perfect harmony with the Tunisian environment.
The descriptive analysis of phenotypic variables showed the existence of predominant characters in the population, such as
the color of the horns, the facial profile and body skin color. This population is made of two races la Blonde du Cap Bon
and la Brune de l’Atlas that divides into two branches tawny and gray. It should be noted, that similar characters to those
of the zebu, such as the presence of the bump and the large navel size, still exists in Tunisia. The analysis of variance
43
App. Sci. Report. 10 (1), 2015: 38-44
revealed 8 significant variables with imada effect and 7 with the gender effect. Multiple Factor analysis (MFA) was
performed by grouping variables into two groups (color and general conformation). This analysis generated two distinct
groups of animals that differ only in the traits related to the horns. Gender and geographical location have not contributed
in this distribution. The discriminant analysis using the locality as discriminatory factor did not allow to group individuals
according to the locality because they form a single group. The hierarchical cluster analysis has revealed three classes. The
only variables that were the cause of this dissimilarity are the profile of the back and orientation of the horns. It is
concluded that the local cattle population is very heterogeneous phenotypically. This variability is neither gender nor the
geographic location implying that the phenotypic diversity within the local cattle population is an individual variability.
References
Atti N, Ben Dhia H.1987. Performances des bovins croisés en Tunisie. Annales de l’Institut National de Recherche Agronomique de Tunisie, 60. pp1-24.
Belhaj MT.1972. Essai d’amélioration génétique de la race bovine locale par le croisement d’absorption. Mémoire de fin d’études de troisième cycle
d’agronomie. Institut National Agronomique de Tunisie.85p.
Ben Belgacem H.2005. Evaluation zootechnique du programme de croisement de la vache locale dans la région de Séjnène. Mémoire de Mastère. Institut
National Agronomique de Tunisie.50p.
Casin P. 1999. Analyse des données et des panels de données. Boeck Université, Paris: 192p. [14-10-2013]. Available from internet<http[14-10-2013].
Available from internet<http://books.google.tn/books?id=dK2uCe2vZ4cC&dq=CAH+des+individu&hl=fr&source=gbs_navlinks_s>.
DGA.1931. Direction Générale de l’Agriculture, du Commerce et de la colonisation. Principales Productions Animales en Tunisie. 150p.
Djemali M, Berger PJ.1992. Yield and reproduction characteristics of Friesian cattle under North African conditions. Journal of Dairy Science. 75,pp
3568 - 3575.
FAO.2008. L’état des ressources zoogénétiques pour l’alimentation et l’agriculture dans le monde. Rome, .39p.
FAO.2012. Characterization of animal genetic resources. Rome, 2012, 87p.
Snayaa N.2003. Inventaire de la biodiversité bovine autochtone en Tunisie et développement d’un schéma de conservation. Mémoire de diplôme d’études
approfondies. Institut National Agronomique de Tunisie.72p.
44