|
|
||||||||
1 Dipartimento di Produzioni Animali, Epidemiologia ed Ecologia, via Leonardo da Vinci 44, 10095 Grugliasco, Torino, Italy
2 Dipartimento di Scienze e Tecnologie Veterinarie per la Sicurezza Alimentare, Università degli Studi di Milano, 20134 Milano, Italy
3 Fondazione Parco Tecnologico Padano, Centro Ricerche e Studi Agroalimentari, 20090 Segrate (Milano), Italy
4 Dipartimento di Sanità e Benessere Animale, Strada Provinciale per Casamassima km 3, 70010 Valenzano (Bari), Italy
Corresponding author: P. Sacchi; e-mail: paola.sacchi{at}unito.it.
| ABSTRACT |
|---|
|
|
|---|
s1-casein (CSN1S1),
s2-casein, (CSN1S2), ß-casein (CSN2), and
-casein (CSN3) loci using several genomic techniques and milk protein analysis. Casein haplotype frequencies were estimated for each breed. Principal component analysis was carried out to highlight the relationship among breeds. Allele and haplotype distributions indicated considerable differences among breeds. The haplotype CSN1S1*F- CSN1S2*F-CSN3*D occurred in all breeds with frequencies >0.100 and was the most common haplotype in the Southern breeds. A high frequency of CSN1S1*0-CSN1S2*C-CSN3*A haplotype was found in Vallesana population (0.162). Principal component analysis clearly separated the Northern and Southern breeds by the first component. The variability of the caprine casein loci and variety of resulting haplotypes should be exploited in the future using specific breeding programs aiming to preserve biodiversity and to select goat genetic lines for specific protein production.
Key Words: goat casein polymorphism haplotype
Abbreviation key: AS-PCR = allele specific-PCR, CSN1S1 =
s1-casein locus, CSN1S2 =
s2-casein locus, CSN3 =
-casein locus, CSN2 = ß-casein locus, IEF = isoelectrofocusing, PCR-SSCP = PCR-single strand conformational polymorphism.
| INTRODUCTION |
|---|
|
|
|---|
s1-casein (CSN1S1), ß-casein (CSN2),
s2-casein (CSN1S2), and
-casein (CSN3) (Ferretti et al., 1990; Threadgill and Womack, 1990; Rijnkels et al., 1997). In goats, the entire casein gene cluster region spans about 250 kb on chromosome 6 (Hayes et al., 1993). Furthermore, CSN1S1 and CSN2 are only 12 kb apart and are convergently transcribed (Leroux and Martin, 1996).
In recent years, the genetic polymorphism of goat caseins has raised considerable research interest because goat casein polymorphisms are related to milk quality, composition, and technological properties (Martin et al., 2002). Researchers started with CSN1S1, which has 16 known co-dominant alleles. Alleles are associated with different rates of protein synthesis. In the 1980s, 7 genetic variants were known at the protein level: A, B, C, D, E, F, and 0 (Boulanger et al., 1984; Grosclaude et al., 1987; Brignon et al., 1989; Mahé and Grosclaude, 1989). Later on, more variants were identified both at the protein and DNA level and named as 02, G, B2, B3, B4 (Leroux et al., 1990; Martin and Leroux, 1994; Grosclaude and Martin, 1997); H, I, L (Chianese et al., 1997); M (Bevilacqua et al., 2002); and N (Ramunno et al., 2002). The B allele, now renamed B1, is considered ancestral for the species. The original 0 allele was renamed 01 to differentiate it from the second null allele (02) found by genomic studies. The D variant was traced back to the G variant (Martin and Leroux, 1994), characterized as the F allele, by an internal deletion. On the basis of the milk content of
S1-casein, the CSN1S1 variants can be classed into 4 groups: strong alleles (A, B1, B2, B3, B4, C, H, L, and M), producing almost 3.5 g/L of
S1-casein each; intermediate alleles (E and I; 1.1 g/L); weak alleles (F and G; 0.45 g/L); and null alleles (01, 02, and N), producing no
S1-casein (Grosclaude and Martin, 1997; Rando et al., 2000).
Three variants of CNS2 were found to be associated with a normal ß-casein content: A, B (Mahé and Grosclaude, 1993), and C (Neveu et al., 2002), the last allele differing for a single amino acid substitution (Ala177
Val177) from the A variant. Furthermore, 2 null CSN2 alleles were identified, both characterized by mutations responsible for premature stop codons in exon 7 (Ramunno et al., 1995: GenBank Accession number AJ011019; Persuy et al., 1999: GenBank Accession number AF172260). The mRNA analysis revealed that the transcript product amounts were almost 10 (Ramunno et al., 1995) and 100 (Persuy et al., 1999) times lower for the null alleles than for the A variant.
The genetic polymorphism of CSN1S2 was first described by Boulanger et al. (1984). At least 8 alleles have been identified. Alleles are associated with 3 synthesis levels; A, B (Boulanger et al., 1984), C (Bouniol et al., 1994), E (Lagonigro et al., 2001), F (Ramunno et al., 2001a), and G (Erhardt et al., 2002) are associated with a normal synthesis level, whereas D and 0 are associated with lower and null synthesis levels, respectively (Ramunno et al., 2001a, b).
Since the discovery of 2 caprine k-casein variants identifiable by isoelectrofocusing (IEF) (Di Luccia et al., 1990) and successively confirmed both at the protein and DNA level by Caroli et al. (2001), further CSN3 polymorphisms were detected by DNA analysis (Yahyaoui et al., 2001; Angiolillo et al., 2002; Yahyaoui et al., 2003; Jann et al., 2004). A total of 13 polymorphic sites were identified in the domestic goat (Jann et al., 2004), allowing the identification of 14 alleles corresponding to 11 protein variants, if taking into account the Capra pyrenaica variant (Yahyaoui et al., 2001), which was later observed in domestic goat breeds (Yahyaoui et al., 2003). Recently, 2 more polymorphisms were demonstrated by PCR-single strand conformational polymorphism (PCR-SSCP) analysis (Chessa et al., 2003), but they are not yet characterized. The number of caprine CSN3 polymorphisms and available sequences has led to some nomenclature problems, but new nomenclature has been proposed (Chessa et al., 2003; Jann et al., 2004), mainly based on the GenBank chronological order of the variants.
Because of the tight association among casein genes, the estimation of the relationship between casein variants and milk production traits can be improved by considering the entire casein haplotype instead of single gene typing. This was first suggested in cattle by Grosclaude et al. (1978) and was further developed in the same species (Aleandri et al., 1990; Ikonen et al., 2001; Boettcher et al., 2004; Caroli et al., 2004). The effect of an allele at a given locus could be confounded with the effect of linked alleles at another locus. Moreover, selection for an allele could automatically increase the frequency of linked alleles that may not result in a favorable effect on the trait of interest. Research focused at the haplotype level is necessary to detect important effects that could be used for the genetic improvement of goat breeds.
The aim of this work was to investigate the genetic structure at the casein gene cluster in 5 Italian goat breeds and to evaluate the casein haplotype variability within and among breeds.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Vallesana is reared in Northern Italy (Piedmont Region) and was originally imported from Switzerland as Walliser Schwarzhalsziege. The animals are medium sized and have black forequarters and white hindquarters. In 2002, the breed risk status was classified as critical because the total population size was around 100 head (http://dad.fao.org). Vallesana is considered a dual-purpose (milk and meat) breed, and the milk is mainly used for fresh cheese manufacturing. Average lactation milk production is about 260 kg in 210 d.
Roccaverano is a population found in Northern Italy. The animals are medium sized, brown or white, and polled. The breed is well adapted to the local environment and currently has a census of about 1000 individuals (http://dad.fao.org). Average lactation milk production is about 350 kg in 198 d.
Maltese was derived from Italian and North African breeds and is distributed in Southern Italy with an estimated population of 7000 head (http://dad.fao.org). The animals are lightweight and normally polled. Maltese hair is long and white or yellowish. Maltese have black marks on the ears and on both sides of the head and neck. Average lactation milk production is about 350 kg in 200 d.
Garganica goats are spread all over the Southern Apennine, and the breed is well adapted to the hilly environment. The animals are black and have screw-shaped horns. The current census indicates a population of about 19,000 breeding females (http://dad.fao.org). The typical lactation milk production is from 200 to 250 kg in 210 d. Garganica milk is characterized by high protein and fat content.
The Jonica breed resulted from crosses of the local population and Maltese. Jonica are almost exclusively raised in the southern area of Apulia. Jonica are adapted to live on arid soils. They are either polled or horned and have long, white hair. The number of breeding females is of about 1600 individuals (http://dad.fao.org). Average lactation milk production is about 300 kg in 215 d.
Genotyping
Genotyping was carried out at the protein level and at the DNA level to take into account both the phenotypic expression of the casein genes and the genetic structure of their nucleotide sequences. First, milk samples were typed by IEF according to Caroli et al. (2001). The DNA was extracted from blood or milk by the GFX Genomic Blood DNA Purification Kit (Amersham Biosciences, Uppsala, Sweden) and then typed using the following analyses:
CSN1S1: Allele-specific-PCR (AS-PCR) (Jànsa Pérez et al., 1994; Leroux et al., 1992); PCR-RFLP (Ramunno et al., 2000, 2002);
CSN1S2: PCR-RFLP (Ramunno et al., 1999, 2001a; Lagonigro et al., 2001);
CSN2: AS-PCR (Ramunno et al., 1995);
CSN3: PCR-SSCP (Chessa et al., 2003).
Statistical Analyses
Genepop (Raymond and Rousset, 1995) software was used to estimate allelic frequencies and to verify Hardy-Weinberg equilibrium. Significance levels were based on F-statistics and were from FSTAT (Goudet, 1995). Casein haplotype frequencies were estimated by using EH (Xie and Ott, 1993), considering only the alleles with frequencies >0.05. Principal component analysis was carried out on haplotype frequencies (SAS, 1990) to highlight the relationship among the breeds on the basis of the casein cluster variability.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
No polymorphism for CSN2 was observed in Roccaverano and Vallesana. Because milk IEF as well as the DNA typing methods do not allow the detection of the CSN2*C allele (Neveu et al., 2002), the A* allele in Table 1
could include this genetic variant, which is associated with a normal ß-casein content. In the Southern breeds, few heterozygous animals carried the null allele identified by Ramunno et al. (1995). No typing was done for the null allele characterized by Persuy et al. (1999), and it was found only in Creole and Pyrenean breeds. Interestingly, a Maltese goat heterozygous for CSN2*0 was homozygous for the F allele at CSN1S1 locus. This indicates the occurrence of a CSN1S1*F-CSN2*0 haplotype, which is associated with a very low casein level even if it is linked to a strong CSN1S2 allele. The casein content of this haplotype is around 2.5 g/L. This haplotype can be added to the 6 quantitative haplotypes observed at the goat calcium-sensitive casein loci (Table 2
) (Rando et al., 2000) on the basis of the approximate casein level (Ramunno et al., 1995; Grosclaude and Martin, 1997; Ramunno et al., 2001a, b). The new haplotype found in Maltese could be exploited in eventual breeding programs aiming at producing milk with specific nutritional properties (i.e., less caloric value) rather than milk used for cheese-making.
|
s2-casein content. The null allele was found in Maltese and Jonica at a very low frequency (0.007 and 0.005, respectively). The typing at both the protein and DNA level allowed us to assess that the CSN1S2*E allele, characterized at the DNA level by Lagonigro et al. (2001), matches with the G variant identified at the protein level by Erhardt et al. (2002). In the study at DNA level, only one nonsynonymous mutation was found in CSN1S2*E, leading the authors to hypothesize an amino acid substitution Pro197 to Arg197 in the mature protein (Lagonigro et al., 2001). The IEF migration pattern found in the study at protein level (Erhardt et al., 2002) clearly suggested the non-correspondence between the observed G variant and the hypothesized E variant. However, the coincidence assessed in the present study between CSN1S2*G, typed at the protein level, and CSN1S2*E, typed at the DNA level, indicates the necessity to analyze the CSN1S2*E protein to identify the modifications responsible for the IEF migration pattern. The usefulness of investigating milk protein loci not only at the DNA level but also taking into account their phenotypic expression is evident.
Extensive polymorphism was observed at CSN3. The most common allele was the CSN3*D according to the nomenclature proposed by Chessa et al. (2003). Frequencies ranged from 0.596 (Vallesana) to 0.776 (Garganica). The allele was referred to as B in Yahyaoui et al. (2003). It will be important to find an international agreement on CSN3 nomenclature in the future. The CSN3*A allele was also found in all breeds. This allele had variable frequencies ranging from 0.066 (Garganica) to 0.373 (Roccaverano). The CSN3*B, which is detectable at the protein level by IEF, results from the Gln44
Arg44 exchange (Caroli et al., 2001) and was observed in the Southern breeds and in Roccaverano. It had the highest frequencies in Maltese (0.179) and Garganica (0.132). Only the Northern populations carried the CNS3*C allele. The frequency was <0.05. In addition, 2 new SSCP patterns were detected, and the relative alleles were named X and Y, according to Chessa et al. (2003). The molecular characterization of the 2 variants is in progress.
Hardy-Weinberg equilibrium was verified for almost all loci within populations. Deviations were found for CSN1S1 and CSN1S2 (Vallesana), CSN1S2 (Maltese), and CSN3 (Jonica). The FST-statistics ranged from 0.013 (CSN2) to 0.074, with an overall value of 0.062, and were significantly >0 for all loci. The FIS-values were significantly >0 in 2 cases: for CSN1S2 in Maltese (0.169) and for CSN1S2 in Vallesana (0.136).
The casein haplotype frequencies are shown in Table 3
. Haplotypes with a frequency >0.05 in at least one breed are reported. The 18 most common haplotypes were identified. When considering each breed separately, frequencies >0.05 were found only for 6 (Vallesana), 7 (Roccaverano), 6 (Maltese), 4 (Jonica), and 9 (Garganica) haplotypes. The expected frequencies under the independence hypothesis are also shown in Table 3
. The
2 test was significantly different from the null hypothesis (no association) in Vallesana, Maltese, and Jonica, but no differences were found for Roccaverano and Garganica, indicating substantial linkage equilibrium in these breeds. These latter 2 breeds also had a higher haplotype number (7 and 9) compared with the other breeds. The haplotype FFD (in the order: CSN1S1*F-CSN1S2*F-CSN3*D) occurred in all breeds with a high frequency (>10%) and was the most common haplotype in the Southern breeds. The estimated frequencies of this haplotype were always much higher than expected frequencies. The distribution of the other haplotypes showed variability among the 5 breeds. In Vallesana and Roccaverano, ECD and FCD haplotypes were most prevalent, respectively. A high frequency of OCA haplotype was found in Vallesana (0.162), indicating that the null allele CSN1S1*0 (frequency of 0.175) is mainly present within the casein cluster. Few haplotypes had expected frequencies higher than 0.1: BFD and BAD in Garganica, FCD and ECD in Vallesana, and ACD in Jonica.
|
|
|
The evolutionary pathway of caprine CSN1S1, first proposed by Grosclaude et al. (1994), was slightly modified by Grosclaude and Martin (1997) and by Chianese et al. (1997). A putative B1 allele was considered ancestral, and 2 divergent lineages were suggested, one leading to A, G, 01, 02, I, and H variants (lineage A), and the other leading to B, L, F, C, and E variants (lineage B). Later on, Bevilacqua et al. (2002) included the M allele in the phylogeny, and they focused on the interallelic recombination event possibly responsible for this new variant.
As far as CSN1S2 is concerned, an evolutionary pathway can be proposed starting from the A variant and leading independently to the B, C, and F alleles, each one characterized by a different amino acid substitution with respect both to CSN1S2*A and to the bovine (Swissprot Accession number P02663) and ovine (P04654) sequences. Thus, CSN1S2*A may be considered as the ancestral variant. If goat CSN1S2*A, B, and C (P33049) and goat CSN1S2*E (CAC21704 sequences are compared, one can postulate that CSN1S2*E variant derives from CSN1S2*C, because the 2 variants share an Ile, instead of Lys, at position 167 of the mature protein.
A phylogeny for CSN3 was proposed by Yahyaoui et al. (2003) and by Jann et al. (2004). Apart from the conflicting nomenclature, both papers agree that CSN3*A appeared later in the evolutionary pathway and that 2 different lineages occurred.
Taking into account the amino acid differences among the protein variants resulting from the 18 casein clusters considered, the most likely ancestral haplotype could be the BAD (Table 5
), followed by BFD, which gives rise to the most common FFD. However, strong evidence of recombination among the casein genes as well as of further intragenic recombination events arise if one considers the possible evolutionary patterns involving the 18 haplotypes.
|
| CONCLUSIONS |
|---|
|
|
|---|
Arg44, revealed in CSN3, should be identified to clarify the functional importance of these substitutions. A screening of goat casein variability at the protein level is recommended for a more complete picture of the genetic polymorphisms in the casein genes.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication July 13, 2004. Accepted for publication November 29, 2004.
| REFERENCES |
|---|
|
|
|---|
-casein gene. J. Dairy Sci. 85:26792680.
S1-casein variant in the goat species. Eur. J. Biochem. 269:12931303.[Medline]
S1 et
S2 de la chèvre (Capra hircus). Genet. Sel. Evol. 16:157176.
S2-casein. Anim. Genet. 25:173177.[Medline]
S1-CN and characterization of those of its genetic variants which are synthesized at a high level,
S1-CN A, B and C. Protein Seq. Data Anal. 2:181188.[Medline]
-casein (CSN3) in different breeds and characterization at DNA level. Anim. Genet. 32:226230.[Medline]
-casein (CSN3) alleles in domestic goat by PCR-SSCP. J. Dairy Sci. 86:37263729.
S1-casein variants in goat milk. Pages 259267 in Milk Protein Polymorphism. International Dairy Federation, Brussels, Belgium.
-casein polymorphism in caprine milk. Sci. Tecn. Latt.-. Cas. 41:305314.
S2-casein (CSN1S2) and evidence for a further allele. Milchwissenschaft 57:137140.
S2 bovine: étroite liaison du locus
S2-Cn avec les loci
S1-Cn, ß-Cn et
-Cn, mise en évidence dun délétion dans le variant
S2-CnD. Ann. Genet. Sel. Anim. 10:313327.
S1-casein. Genet. Sel. Evol. 19:399412.
s1 caprine, ses effets, son évolution. INRA Prod. Anim. 7:319.
-casein (CSN3) gene from wild and domestic caprine species revealed by DNA sequencing. J. Dairy Res. 71:188195.[Medline]
S1-casein polymorphism in Spanish goat breeds. Small Rumin. Res. 20:215221.
S1-casein alleles. A potential tool in selection of individuals carrying alleles associated with a high level protein synthesis. Anim. Genet. 21:341351.[Medline]
S1-CnD, another allele associated with a decreased synthesis rate at the caprine
S1-casein locus. Genet. Sel. Evol. 21:127129.
This article has been cited by other articles:
![]() |
G. De la Torre, J. M. Serradilla, F. Gil Extremera, and M. R. Sanz Sampelayo Nutritional Utilization in Malaguena Dairy Goats Differing in Genotypes for the Content of {alpha}S1-Casein in Milk J Dairy Sci, June 1, 2008; 91(6): 2443 - 2448. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Caroli, F. Chiatti, S. Chessa, D. Rignanese, E. M. Ibeagha-Awemu, and G. Erhardt Characterization of the Casein Gene Complex in West African Goats and Description of a New {alpha}s1-Casein Polymorphism J Dairy Sci, June 1, 2007; 90(6): 2989 - 2996. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Caroli, F. Chiatti, S. Chessa, D. Rignanese, P. Bolla, and G. Pagnacco Focusing on the goat casein complex. J Dairy Sci, August 1, 2006; 89(8): 3178 - 3187. [Abstract] [Full Text] [PDF] |
||||
![]() |
E.-M. Prinzenberg, K. Gutscher, S. Chessa, A. Caroli, and G. Erhardt Caprine {kappa}-Casein (CSN3) Polymorphism: New Developments in Molecular Knowledge J Dairy Sci, April 1, 2005; 88(4): 1490 - 1498. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |