Bouabid Badaoui

Integrating Y-Chromosome, Mitochondrial, and Autosomal Data to Analyze the Origin of Pig Breeds

We have investigated the origin of swine breeds through the joint analysis of mitochondrial, microsatellite, and Y-chromosome polymorphisms in a sample of pigs and wild boars with a worldwide distribution. Genetic differentiation between pigs and wild boars was remarkably weak, likely as a consequence of a sustained gene flow between both populations. The analysis of nuclear markers evidenced the existence of a close genetic relationship between Near Eastern and European wild boars making it difficult to infer their relative contributions to the gene pool of modern European breeds. Moreover, we have shown that European and Far Eastern pig populations have contributed maternal and paternal lineages to the foundation of African and South American breeds. Although West African pigs from Nigeria and Benin exclusively harbored European alleles, Far Eastern and European genetic signatures of similar intensity were detected in swine breeds from Eastern Africa. This region seems to have been a major point of entry of livestock species in the African continent as a result of the Indian Ocean trade. Finally, South American creole breeds had essentially a European ancestry although Asian Y-chromosome and mitochondrial haplotypes were found in a few Nicaraguan pigs. The existence of Spanish and Portuguese commercial routes linking Asia with America might have favored the introduction of Far Eastern breeds into this continent.

Short communication: Effect of αS1-casein (CSN1S1) and κ-casein (CSN3) genotypes on milk composition in Murciano-Granadina goats

The effects of the caprine alpha(S1)-casein (CSN1S1) polymorphisms on milk quality have been widely demonstrated. However, much less is known about the consequences of the kappa-casein (CSN3) genotype on milk composition in goats. Moreover, the occurrence of interactions between CSN3 and CSN1S1 genotypes has not been investigated. In this study, an association analysis between CSN1S1 and CSN3 genotypes and milk quality traits was performed in 89 Murciano-Granadina goats. Total milk yield as well as total protein, fat, solids-not-fat, lactose, alpha(S1)-casein (CSN1S1), and alpha(S2)-casein (CSN1S2) contents were recorded every other month during a whole lactation (316 observations). Data analysis using a linear mixed model for repeated observations revealed no interaction between the CSN1S1 and CSN3 genotypes. With regard to the effect of the CSN3 locus, AB and BB genotypes were significantly associated with higher levels of total casein and protein content compared with the AA CSN3 genotype. In strong contrast with French breeds, the CSN1S1 genotype did not affect protein, casein, and fat concentrations in Murciano-Granadina goats. These results highlight the importance of taking into consideration the CSN3 genotype when performing selection for milk composition in dairy goats.

Mitochondrial DNA diversity and origins of South and Central American goats.

We have analysed the genetic diversity of South and Central American (SCA) goats by partially sequencing the mitochondrial control region of 93 individuals with a wide geographical distribution. Nucleotide and haplotype diversities reached values of 0.020 +/- 0.00081 and 0.963 +/- 0.0012 respectively. We have also observed a rather weak phylogeographic structure, with almost 69% of genetic variation included in the within-breed variance component. The topology of a median-joining network analysis including 286 European, Iberian, Atlantic and SCA mitochondrial sequences was very complex, with most of the haplotypes forming part of independent small clusters. SCA sequences showed a scattered distribution throughout the network, and clustering with Spanish and Portuguese sequences occurred only occasionally, not allowing the distinguishing of a clear Iberian signature. Conversely, we found a prominent cluster including Canarian, Chilean, Argentinian and Bolivian mitochondrial haplotypes. This result was independently confirmed by constructing a Bayesian phylogenetic tree (posterior probability of 0.97). Sharing of mitochondrial haplotypes by SCA and Canarian goats suggests that goat populations from the Atlantic archipelagos, where Spanish and Portuguese ships en route to the New World used to stow food and supplies, participated in the foundation of SCA caprine breeds.

Effects of α s1-casein (CSN1S1) and κ-casein (CSN3) genotypes on milk coagulation properties in Murciano-Granadina goats.

The effects of the caprine α s1-casein (CSN1S1) polymorphisms on milk quality and cheese yield have been widely studied in French and Italian goat breeds. Much less is known about the consequences of κ-casein (CSN3) genotype on the technological and coagulation properties of goat milk. In the current study, we have performed an association analysis between polymorphisms at the goat CSN1S1 and CSN3 genes and milk coagulation (rennet coagulation time, curdling rate and curd firmness) and technological (time to cutting of curd and cheese yield) properties. In this analysis, we have included 193 records from 74 Murciano-Granadina goats (with genotypes constituted by different combinations of alleles B, E and F of the gene CSN1S1 and alleles A and B of the gene CSN3) distributed in three herds, which were collected bimonthly during a whole lactation. Data analysis, using a linear mixed model for repeated observations, revealed significant associations between CSN1S1 genotypes and the rate of the curdling process. In this way, milk from EE goats had a significantly higher curdling rate than milk from BB individuals (P<0·05). Contrary to previous experiments performed in French breeds, cheese yield was not significantly different in BB, EE and EF goats. Moreover, we have shown that CSN3 genotype has a significant effect on the rennet coagulation time (BB>AB, P<0·05) but not on cheese yield. No interaction between the CSN1S1 and CSN3 genotypes was observed.

Identification of a single nucleotide polymorphism at intron 16 of the caprine acyl-coenzyme A: diacylglycerol acyltransferase 1 (DGAT1) gene.

The DGAT1 gene encodes a microsomal enzyme that catalyses the only committed step in triacylglycerol synthesis by joining diacylglycerol and fatty acyl coenzyme A. In cattle, a K232A substitution in the DGAT1 molecule has a significant effect on enzyme activity and milk fat content. The prominent role of this gene in lipid metabolism led us to undertake the structural characterization of DGAT1 in goats. In this way, we have sequenced a 1552 bp fragment of the goat DGAT1 cDNA, which encompasses most of the coding sequence (from exon 1 to 17), and a genomic fragment covering exons 12 to 17. Multiple alignment of the goat DGAT1 sequences revealed the existence of a single nucleotide polymorphism (SNP) involving a T to C substitution at intron 16. We optimized a primer extension based genotyping method that allowed us to determine that the C variant is a minority allele with frequencies ranging from 0.062 (Murciano-Granadina) to 0.109 (Malagueña). This SNP, although not expected to have any functional effect, might be useful as a genetic marker in association studies to detect additional DGAT1 polymorphisms which might influence fat milk content and other traits of economic interest.

Detecting the existence of gene flow between Spanish and North African goats through a coalescent approach

Human-driven migrations are one of the main processes shaping the genetic diversity and population structure of domestic species. However, their magnitude and direction have been rarely analysed in a statistical framework. We aimed to estimate the impact of migration on the population structure of Spanish and African goats. To achieve this goal, we analysed a dataset of 1,472 individuals typed with 23 microsatellites. Population structure of African and Spanish goats was moderate (mean FST = 0.07), with the exception of the Canarian and South African breeds that displayed a significant differentiation when compared to goats from North Africa and Nigeria. Measurement of gene flow with Migrate-n and IMa coalescent genealogy samplers supported the existence of a bidirectional gene flow between African and Spanish goats. Moreover, IMa estimates of the effective number of migrants were remarkably lower than those calculated with Migrate-n and classical approaches. Such discrepancies suggest that recent divergence, rather than extensive gene flow, is the main cause of the weak population structure observed in caprine breeds.

Technical note: Advantages and limitations of authenticating Palmera goat dairy products by pyrosequencing the melanocortin 1 receptor (MC1R) gene

Inferring the breed of origin of dairy products can be achieved through molecular analysis of genetic markers with a population-specific pattern of segregation. The goal of the current work was to generate such markers in goats by resequencing several pigmentation genes [melanocortin 1 receptor (MC1R), v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (KIT), tyrosinase (TYR), and tyrosinase-related protein 2 (TYRP2)]. This experiment revealed 10 single nucleotide polymorphisms (SNP), including 5 missense mutations and 1 nonsense mutation. These markers were genotyped in 560 goats from 18 breeds originally from Italy, the Iberian Peninsula, the Canary Islands, and North Africa. Although the majority of SNP segregated at moderate frequencies in all populations (including 2 additional markers that were used as a source of information), we identified a c.764G>A SNP in MC1R that displayed highly divergent allelic frequencies in the Palmera breed compared with the Majorera and Tinerfeña breeds from the Canary Islands. Thus, we optimized a pyrosequencing-based technique that allowed us to estimate, very accurately, the allele frequencies of this marker in complex DNA mixtures from different individuals. Once validated, we applied this method to generating breed-specific DNA profiles that made it possible to detect fraudulent cheeses in which Palmero cheese was manufactured with milk from Majorera goats. One limitation of this approach, however, is that it cannot be used to detect illegal manufacturing where Palmero dairy products are produced by mixing milk from Palmera and Majorera goats, because the c.764G>A SNP segregates in both breeds.

Identification of c.483C>T polymorphism in the caprine tyrosinase-related protein 1 (TYRP1) gene

Tyrosinase-related protein 1 (TYRP1) has been shown to play a fundamental role in pigmentation both in human and mouse. In this work, we aimed to characterize the variability of the caprine TYRP1 gene and investigate its segregation in a wide array of goat breeds. By partially sequencing the coding region of the TYRP1 gene in 18 individuals from eight different breeds, we were able to identify a synonymous nucleotide substitution at exon 3 (c.483C>T). An extensive survey of Iberian and Balearic (N=175), Italian (N=99), Swiss (N=54), Asian (N=14), Canarian (N=92) and North African (N=117) goats with different coat colours was carried out. We found that the C-allele has a different distribution in European vs African breeds, being almost fixed in the latter. Moreover, the C-allele showed an increased frequency in white coated breeds (Girgentana, Grigia Molisana, Blanca de Rasquera and Saanen) when compared with those displaying a dark pigmentation (Cilentana Nera, Azpi Gorri and Murciano-Granadina). This could be due to genetic drift, migration and other factors associated with the demographic history of breeds under analysis or to a genetic hitchhiking event (c.483C>T frequencies would be shaped by a neighbouring causal mutation differentially selected in white and black goats). More refined studies will be needed to distinguish between these two alternative explanations.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.