Investigating the footprint of post-domestication dispersal on the diversity of modern European, African and Asian goats.

BACKGROUND: Goats were domesticated in the Fertile Crescent about 10,000 years before present (YBP) and subsequently spread across Eurasia and Africa. This dispersal is expected to generate a gradient of declining genetic diversity with increasing distance from the areas of early livestock management. Previous studies have reported the existence of such genetic cline in European goat populations, but they were based on a limited number of microsatellite markers. Here, we have analyzed data generated by the AdaptMap project and other studies. More specifically, we have used the geographic coordinates and estimates of the observed (Ho) and expected (He) heterozygosities of 1077 European, 1187 African and 617 Asian goats belonging to 38, 43 and 22 different breeds, respectively, to find out whether genetic diversity and distance to Ganj Dareh, a Neolithic settlement in western Iran for which evidence of an early management of domestic goats has been obtained, are significantly correlated. RESULTS: Principal component and ADMIXTURE analyses revealed an incomplete regional differentiation of European breeds, but two genetic clusters representing Northern Europe and the British-Irish Isles were remarkably differentiated from the remaining European populations. In African breeds, we observed five main clusters: (1) North Africa, (2) West Africa, (3) East Africa, (4) South Africa, and (5) Madagascar. Regarding Asian breeds, three well differentiated West Asian, South Asian and East Asian groups were observed. For European and Asian goats, no strong evidence of significant correlations between Ho and He and distance to Ganj Dareh was found. In contrast, in African breeds we detected a significant gradient of diversity, which decreased with distance to Ganj Dareh. CONCLUSIONS: The detection of a genetic cline associated with distance to the Ganj Dareh in African but not in European or Asian goat breeds might reflect differences in the post-domestication dispersal process and subsequent migratory movements associated with the management of caprine populations from these three continents.

Detailed mineral profile of milk, whey, and cheese from cows, buffaloes, goats, ewes and dromedary camels, and efficiency of recovery of minerals in their cheese.

Milk and dairy products are important in the human diet not only for the macro nutrients, such as proteins and fats, that they provide, but also for the supply of essential micronutrients, such as minerals. Minerals are present in milk in soluble form in the aqueous phase and in colloidal form associated with the macronutrients of the milk. These 2 forms affect the nutritional functions of the minerals and their contribution to the technological properties of milk during cheese-making. The aim of the present work was to study and compare the detailed mineral profiles of dairy foods (milk, whey, and cheese) obtained from cows, buffaloes, goats, ewes and dromedary camels, and to analyze the recovery in the curd of the individual minerals according to a model cheese-making procedure applied to the milk of these 5 dairy species. The detailed mineral profile of the milk samples was obtained by inductively coupled plasma - optical emission spectroscopy (ICP - OES). We divided the 21 minerals identified in the 3 different matrices into essential macro- and micro-minerals, and environmental micro-minerals, and calculated the recovery of the individual minerals in the cheeses. The complete mineral profiles and the recoveries in the cheeses were then analyzed using a linear mixed model with Species and Food, and their interaction included as fixed effects, and Sample within Species as a random effect. The mineral profiles of each food matrix were then analyzed separately with a general linear model in which only the fixed effect of Species was included. The results showed that the species could be divided into 2 groups: those producing a more diluted milk characterized by a higher content of soluble minerals (in particular K), and those with a more concentrated milk with a higher colloidal mineral content in the skim of the milk (such as Ca and P). The recoveries of the minerals in the curd were in line with the initial content in the milk, and also highlighted the fact that the influence of the brine was not limited to the Na content but to its whole mineral makeup. These results provide valuable information for the evaluation of the nutritional and technological properties of milk, and for the uses made of the byproducts of cheese making from the milk of different species.

Cheese yield and nutrients recovery in the curd predicted by Fourier-transform spectra from individual sheep milk samples.

The objectives of this study were to explore the use of Fourier-transform infrared (FTIR) spectroscopy on individual sheep milk samples for predicting cheese-making traits, and to test the effect of the farm variability on their prediction accuracy. For each of 121 ewes from 4 farms, a laboratory model cheese was produced, and 3 actual cheese yield traits (fresh cheese, cheese solids, and cheese water) and 4 milk nutrient recovery traits (fat, protein, total solids, and energy) in the curd were measured. Calibration equations were developed using a Bayesian approach with 2 different scenarios: (1) a random cross-validation (80% calibration; 20% validation set), and (2) a leave-one-out validation (3 farms used as calibration, and the remaining one as validation set) to assess the accuracy of prediction of samples from external farms, not included in calibration set. The best performance was obtained for predicting the yield and recovery of total solids, justifying for the practical application of the method at sheep population and dairy industry levels. Performances for the remaining traits were lower, but still useful for the monitoring of the milk processing in the case of fresh curd and recovery of energy. Insufficient accuracies were found for the recovery of protein and fat, highlighting the complex nature of the relationships among the milk nutrients and their recovery in the curd. The leave-one-out validation procedure, as expected, showed lower prediction accuracies, as a result of the characteristics of the farming systems, which were different between calibration and validation sets. In this regard, the inclusion of information related to the farm could help to improve the prediction accuracy of these traits. Overall, a large contribution to the prediction of the cheese-making traits came from the areas known as "water" and "fingerprint" regions. These findings suggest that, according to the traits studied, the inclusion of water regions for the development of the prediction equation models is fundamental to maintain a high prediction accuracy. However, further studies are necessary to better understand the role of specific absorbance peaks and their contribution to the prediction of cheese-making traits, to offer reliable tools applicable along the dairy ovine chain.

Composition and aptitude for cheese-making of milk from cows, buffaloes, goats, sheep, dromedary camels, and donkeys.

Bovines produce about 83% of the milk and dairy products consumed by humans worldwide, the rest represented by bubaline, caprine, ovine, camelid, and equine species, which are particularly important in areas of extensive pastoralism. Although milk is increasingly used for cheese production, the cheese-making efficiency of milk from the different species is not well known. This study compares the cheese-making ability of milk sampled from lactating females of the 6 dairy species in terms of milk composition, coagulation properties (using lactodynamography), curd-firming modeling, nutrients recovered in the curd, and cheese yield (through laboratory model-cheese production). Equine (donkey) milk had the lowest fat and protein content and did not coagulate after rennet addition. Buffalo and ewe milk yielded more fresh cheese (25.5 and 22.9%, respectively) than cow, goat, and dromedary milk (15.4, 11.9, and 13.8%, respectively). This was due to the greater fat and protein contents of the former species with respect to the latter, but also to the greater recovery of fat in the curd of bubaline (88.2%) than in the curd of camelid milk (55.0%) and consequent differences in the recoveries of milk total solids and energy in the curd; protein recovery, however, was much more similar across species (from 74.7% in dromedaries to 83.7% in bovine milk). Compared with bovine milk, the milk from the other Artiodactyla species coagulated more rapidly, reached curd firmness more quickly (especially ovine milk), had a more pronounced syneresis (especially caprine milk), had a greater potential asymptotical curd firmness (except dromedary and goat milk), and reached earlier maximum curd firmness (especially caprine and ovine milk). The maximum measured curd firmness was greater for bubaline and ovine milk, intermediate for bovine and caprine milk, and lower for camelid milk. The milk of all ruminant species can be used to make cheese, but, to improve efficiency, cheese-making procedures need to be optimized to take into account the large differences in their coagulation, curd-firming, and syneresis properties.

Haplotype structure of MSTN, IGF1, and BMP2 genes in Tunisian goats (Capra hircus) and their association with morphometric traits.

The aim of the present study was to evaluate variability and haplotype structure of twenty-eight single nucleotide polymorphisms (SNPs) at myostatin (MSTN), insulin-like growth factor 1 (IGF1), and bone morphogenetic protein 2 (BMP2) genes. Association between the polymorphic SNPs and morphometric traits was performed on a population of 263 Tunisian goats. The SNPs analyzed were all polymorphic (except one), and the three genes had different haplotype structures. Significant association of SNPs at MSTN with head length was highlighted in the Tunisian goats. The variability at IGF1 gene was associated with body length, ear length, tail length, and chest depth. For BMP2, significant association was revealed with chest depth. Significant association was also detected between linkage disequilibrium (LD) block 2 at IGF1 with body length. These findings might play a potential role in gene-assisted programs.

Breed of goat affects the prediction accuracy of milk coagulation properties using Fourier-transform infrared spectroscopy.

The prediction of traditional goat milk coagulation properties (MCP) and curd firmness over time (CFt) parameters via Fourier-transform infrared (FTIR) spectroscopy can be of significant economic interest to the dairy industry and can contribute to the breeding objectives for the genetic improvement of dairy goat breeds. Therefore, the aims of this study were to (1) explore the variability of milk FTIR spectra from 4 goat breeds (Camosciata delle Alpi, Murciano-Granadina, Maltese, and Sarda), and to assess the possible discriminant power of milk FTIR spectra among breeds, (2) assess the viability to predict coagulation traits by using milk FTIR spectra, and (3) quantify the effect of the breed on the prediction accuracy of MCP and CFt parameters. In total, 611 individual goat milk samples were used. Analysis of variance of measured MCP and CFt parameters was carried out using a mixed model including the farm and pendulum as random factors, and breed, parity, and days in milk as fixed factors. Milk spectra for each goat were collected over the spectral range from wavenumber 5,011 to 925 × cm-1. Discriminant analysis of principal components was used to assess the ability of FTIR spectra to identify breed of origin. A Bayesian model was used to calibrate equations for each coagulation trait. The accuracy of the model and the prediction equation was assessed by cross-validation (CRV; 80% training and 20% testing set) and stratified CRV (SCV; 3 breeds in the training set, one breed in the testing set) procedures. Prediction accuracy was assessed by using coefficient of determination of validation (R2VAL), the root mean square error of validation (RMSEVAL), and the ratio performance deviation. Moreover, measured and FTIR predicted traits were compared in the SCV procedure by assessing their least squares means for the breed effect, Pearson correlations, and variance heteroscedasticity. Results showed the feasibility of using FTIR spectra and multivariate analyses to correctly assign milk samples to their breeds of origin. The R2VAL values obtained with the CRV procedure were moderate to high for the majority of coagulation traits, with RMSEVAL and ratio performance deviation values increasing as the coagulation process progresses from rennet addition. Prediction accuracy obtained with the SCV were strongly influenced by the breed, presenting general low values restricting a practical application. In addition, the low Pearson correlation coefficients of Sarda breed for all the traits analyzed, and the heteroscedastic variances of Camosciata delle Alpi, Murciano-Granadina, and Maltese breeds, further indicated that it is fundamental to consider the differences existing among breeds for the prediction of milk coagulation traits.

Variability of the caprine whey protein genes and their association with milk yield, composition and renneting properties in the Sarda breed: 2. The BLG gene.

The variability of the promoter region and the 3'UTR (exon-7) of the BLG gene, encoding the ß-lactoglobulin, was investigated by sequencing in 263 lactating Sarda goats in order to assess its association with milk traits. Milk traits included: milk yield, fat, total protein and lactose content, pH, daily fat and protein yield (DFPY), freezing point, milk energy, somatic cell count, total microbial mesophilic count, rennet coagulation time (RCT), curd firming rate (k20) and curd firmness (a30). A total of 7 polymorphic sites were detected and the sequence analysed was given accession number KM817769. Only three SNPs (c.-381C>T, c.-323C>T and c.*420C>A) had minor allele frequency higher than 0.05. The effects of farm, stage of lactation and the interaction farm × stage of lactation significantly influenced all the milk traits (P T and c.*420C>A (P T (P < 0.001). The c.-381TT homozygous goats showed lower pH, RCT and k20 than c.-381CT (P < 0.05). In conclusion the polymorphism of the goat BLG gene did not affect the total protein content of the Sarda goat milk, and only weakly influenced RCT and k20. On the other hand, an interesting effect on milk yields and DFPY emerged in two SNPs. This information might be useful in dairy goat breeding programs.

Variability of the caprine whey protein genes and their association with milk yield, composition and renneting properties in the Sarda breed. 1. The LALBA gene.

The 5' flanking region and 3' UTR of the caprine LALBA gene were analysed by SSCP and sequencing. A total of nine SNPs were detected: three in the promoter region, two were synonymous coding SNPs at exon-1, and four SNPs were in exon-4, within the 3'UTR. The nucleotide changes located in the promoter region (c.-358T>C, c.-163G>A, c.-121T>G) were genotyped by SSCP in 263 Sarda goats to evaluate their possible effect on milk yield, composition and renneting properties. We observed an effect of the three SNPs on milk yield and lactose content. Genotypes TT and CT at c.-358T>C (P A (P C and c.-121T>G were part of transcription factors binding sites, potentially involved in modulating the LALBA gene expression. The LALBA genotype affected renneting properties (P < 0.001), as heterozygotes c.-358CT and c.-163GA were characterised by delayed rennet coagulation time and curd firming time and the lowest value of curd firmness. The present investigation increases the panel of SNPs and adds new information about the effects of the caprine LALBA gene polymorphism.

The sheep growth hormone gene polymorphism and its effects on milk traits.

Growth hormone (GH) is encoded by the GH gene, which may be single copy or duplicate in sheep. The two copies of the sheep GH gene (GH1/GH2-N and GH2-Z) were entirely sequenced in one 106 ewes of Sarda breed, in order to highlight sequence polymorphisms and investigate possible association between genetic variants and milk traits. Milk traits included milk yield, fat, protein, casein and lactose percentage. We evidenced 75 nucleotide changes. Transcription factor binding site prediction revealed two sequences potentially recognised by the pituitary-specific transcription factor POU1FI at the GH1/GH2-N gene, which were lost at the promoter of GH2-Z, which might explain the different tissues of expression of GH1/GH2-N (pituitary) and GH2-Z (placenta). Significant differences in milk traits were observed among genotypes at polymorphic loci only for the GH2-Z gene. Sheep with homozygote genotype ss748770547 CC had higher fat percentage (P < 0.01) than TT. SNP ss748770547 was part of a potential transcription factor binding site for C/EBP alpha (CCAAT/Enhancer Binding Protein), which is involved in the regulation of adipogenesis and adipoblast differentiation. SNP ss748770547, located in the GH2-Z gene 5' flanking region, may be a causal mutation affecting milk fat content. These findings might contribute to the knowledge of the sheep GH locus and might be useful in selection processes in sheep.

Evolution of satellite DNA sequences in two tribes of Bovidae: A cautionary tale.

Two clones, Bt1 from Bos taurus and Om1 from Ovis orientalis musimon, were used as probes for hybridization on genomic DNA and on metaphase chromosomes in members of Bovini and Caprini tribes. Bt1 and Om1 are sequences respectively belonging to the 1.715 and 1.714 DNA satellite I families. Southern blots and fluorescence in situ hybridization experiments showed completely coherent results: the Bovini probe Bt1 hybridized only to members of the Bovini tribe and not to members of Caprini. Likewise, the Caprini probe Om1 hybridized only to members of the Caprini tribe and not to members of Bovini. Hybridization signals were detected in the heterochromatic regions of every acrocentric autosome, except for two pairs of autosomes from Capra hircus that did not show hybridization to probe Om1. No signal was detected on X and Y chromosomes or on bi-armed autosomes. Remarkably, probe Om1 showed almost 100% homology with a bacterial sequence reported in Helicobacter pylori.

Identification of novel SNPs in the Sarda breed goats POU1F1 gene and their association with milk productive performance.

The aim of the study was to detect polymorphism in the POU1F1 gene in Sarda breed goat, as well as to establish if SNPs could be associated with milk productive traits. The research was conducted on 129 Sarda breed goats from 4 to 5 years old, multiparous, lactating and in their third to fifth lactation. We report nine exonic and seven non-coding regions SNPs within the Sarda goat POU1F1 gene, namely, Ex 1 61 G>C; Ex 1 108 G>A; Ex 3 C>T; Ex 3 92 C>T; Ex 4 110 A>G; Ex 5 34 G>A resulting in Arg213Lys change; IVS4 641 G>A, IVS4 643 A>C, IVS4 659 G>A, IVS4 677 A>C, IVS4 G699Del, IVS4 709 C>G, Ex 6 17 G>A resulting in Arg228Ser change, Ex 6 58 G>T, Ex 6 172 T>C, 3'UTR 110 T>C. A statistically significant association was found between genotype TT, in position 17 of the exon 6 (3.1 % of frequency), and increased milk yield (P < 0.01) while genotype GT (25.6 % of frequency) was associated with a higher fat content. Genotype TT in position 58 of the exon 6 (3.9 % of frequency) was found to be associated with a higher fat (P < 0.01) and protein content (P < 0.05). Twenty-eight haplotypes were detected, but no significant association between the haplotypes and the milk production traits have been found. Our data, as well as providing new SNPs extending the POU1F1 gene characterization, evidence a relationship between polymorphism and milk production traits in Sarda goat breed.

Development of a RNA extraction method from milk for gene expression study in the mammary gland of sheep.

The aim of the study was to develop a reliable method for the RNA extraction from milk of Sarda sheep breed and to highlight if the extracted RNA can be used for expression study on mammary genes involved in milk fat synthesis using RT-qPCR. The main result is that a sample of 150 ml of milk provides an optimal amount of RNA (73.5 µg/ml). The highest RNA concentration has been found in the samples analysed within 4 h after collection. The RNA extracted was positively correlated to the number of somatic cells (P < 0.001). The efficiency of the extraction method was confirmed by the results obtained from qPCR which showed a Ct value, for SREBPF1 gene of 26.8 ± 0.15. This research demonstrated that the high-quality of the RNA obtained is suited to use for studies of mammary genes expression in sheep, avoiding any damage caused by mammary gland biopsy.

Effect of CSN1S1 gene polymorphism and stage of lactation on milk yield and composition of extensively reared goats.

The effect of CSN1S1 genotype and lactation stage on milk yield and composition were investigated in 80 extensively reared goats. Milk yield was recorded in early, mid and late lactation and individual milk samples were collected to determine: fat, protein, lactose and casein content, pH, freezing point, somatic cell count (SCC) and total microbic mesophilic count (TMC). Relative casein composition and amino acid profile were quantified by HPLC. Fatty acid profile was measured by gas-chromatography. Genotype did not affect milk yield, while this trait was significantly affected by lactation stage (P < 0.01). CSN1S1 BB goats produced significantly higher protein and casein percentages (P < 0.05). αs1-casein (CN) was significantly higher in BB and AB goats than AF and BF, showing intermediate values in AA goats (P < 0.01). The protein percentage and the αs1 and αs2-CN fractions were not affected by lactation stage, while the casein content and the ß and κ-CN significantly increased throughout lactation (P < 0.01). C4 : 0 and C6 : 0 were not affected by genotype, while C8 : 0 and C10 : 0 were higher in the AA goats than BB; most of the long chain FA were higher in BB than AA goats. MUFA and PUFA increased in late lactation. In addition, BB goats showed higher essential amino acids, resulting in an optimal composition from the nutritional point of view, when compared with AA goats. The increase of MUFA, PUFA, essential and cis-FA in late lactation indicate that the lipid composition of goat's milk, with the progress of lactation, tends to improve its nutritional value.

Growth hormone gene variability and its effects on milk traits in primiparous Sarda goats.

The aim of this research was to investigate variability in each of the five exons of the caprine growth hormone (gGH) gene, in order to establish the possible relationships with milk traits in Sarda breed goat. The general linear model procedure was used to analyse the effects of the single strand conformation (SSCP) profiles on milk traits of 100 lactating goats. Analysis of conformational polymorphism at exons 1-5 revealed a total of 25 differing banding patterns. Sequencing revealed 21 nucleotide changes (compared with GenBank D00476): 14 were polymorphic and 7 monomorphic; 19 in exonic regions, 5 of which were nonsynonymous. A SNP upstream of the transcription initiation codon (c.-3A>G) and an indel (c.*29_30insC) in the 3'UTR, were detected. Alignment of 4 cloned sequences including the entire gGH gene led to the identification of 22 nucleotide variations within the intron regions, including two indels. Association analysis revealed that each exon, except exon-1, affected milk yield, exons 1 and 3 influenced milk fat percentage, and all exons, except exon-2, had an effect on protein percentage, supporting previous results in livestock. The variability detected at the caprine GH gene might provide useful information for the phylogeny of ruminants and, more importantly, have implications on the biological function of the growth hormone and on those traits resulting from its physiological action, including milk production and composition. The caprine GH gene may become a useful molecular marker for a more effective genetic selection for milk production traits in goats.

Variations in Casein Genes Are Associated with Milk Protein and Fat Contents in Sarda Goats (Capra hircus), with an Important Role of CSN1S2 for Milk Yield.

This work aimed to assess the variability of casein genes in a population of 153 bucks and 825 lactating does of the Sarda breed, and to perform association analysis between polymorphic sites and milk yield and composition traits. To genotype the casein genes, we chose an SNP panel including 44 SNPs mapping to the four casein genes CSN1S1, CSN2, CSN1S2, and CSN3. Genotyping (made by KASP™ genotyping assay, based on competitive allele-specific PCR) revealed the high variability of the Sarda goat, and haplotype analysis revealed linkage disequilibrium (LD) between CSN1S1 and CSN2 genes, in addition to two LD blocks within the CSN1S2 and two LD blocks within the CSN3 gene, in bucks and does. Association analysis revealed that variability at all four casein genes was associated with milk protein content, total solids, and milk energy. The three Ca-sensitive casein genes were associated with lipid content, and CSN1S2 showed a unique pattern, with intron variants associated with milk yield, in addition to milk pH, NaCl, and SCS (Somatic Cell Score). This information might prove useful in selection schemes and in future investigations aiming to better understand the biology of lactation, and the direct link between genotype and phenotype.

Phenotypic Analysis of Fourier-Transform Infrared Milk Spectra in Dairy Goats.

The infrared spectrum of bovine milk is used to predict many interesting traits, whereas there have been few studies on goat milk in this regard. The objective of this study was to characterize the major sources of variation in the absorbance of the infrared spectrum in caprine milk samples. A total of 657 goats belonging to 6 breeds and reared on 20 farms under traditional and modern dairy systems were milk-sampled once. Fourier-transform infrared (FTIR) spectra were taken (2 replicates per sample, 1314 spectra), and each spectrum contained absorbance values at 1060 different wavenumbers (5000 to 930 × cm-1), which were treated as a response variable and analyzed one at a time (i.e., 1060 runs). A mixed model, including the random effects of sample/goat, breed, flock, parity, stage of lactation, and the residual, was used. The pattern and variability of the FTIR spectrum of caprine milk was similar to those of bovine milk. The major sources of variation in the entire spectrum were as follows: sample/goat (33% of the total variance); flock (21%); breed (15%); lactation stage (11%); parity (9%); and the residual unexplained variation (10%). The entire spectrum was segmented into five relatively homogeneous regions. Two of them exhibited very large variations, especially the residual variation. These regions are known to be affected by the absorbance of water, although they also exhibited wide variations in the other sources of variation. The average repeatability of these two regions were 45% and 75%, whereas for the other three regions it was about 99%. The FTIR spectrum of caprine milk could probably be used to predict several traits and to authenticate the origin of goat milk.

Analysis of polymorphism within POU1F1 gene in relation to milk production traits in dairy Sarda sheep breed.

The ovine POU1F1 gene is localized on chromosome 1 and it contains five introns and six exons. In different mammalian species some mutations in different exons are associated with different production traits. The aim of our research was to study the POU1F1 gene nucleotide sequence to detect possible polymorphisms and their relationships with milk productive traits in Sarda breed sheep. The study had been conducted on 140 ewes, 4 or 5 years old coming from a farm located in Sardinia. All the animals were multiparous, lactating and in their third to fifth lactation. Individual milk yield had been recorded monthly and for each sample fat, protein, casein, lactose, and somatic cell count values were analysed. A jugular blood sample was collected from each ewe to perform genomic DNA extraction. PCR, SSCP and sequencing analysis were carried out to examine the six exons to highlight possible SNPs. One-way ANOVA was used to analyse association of variants with milk yield and/or its composition. Two novel SNP were found: 121 C>T in the 5'UTR of the fourth intron fragment and 249 G>A in the 3'UTR of the sixth exon fragment. The statistical analysis did not shown association between milk productive traits and the found polymorphisms. However, further investigations about the promoter region or the prophet genes, like the PROP-1, could clarify its exact role in regulating the productive traits in sheep.

Mitochondrial DNA diversity of the Sardinian local cattle stock.

The aim of this research was to characterize the genetic diversity of the Sarda (Sa, n = 131), Sardo Bruna (SB, n = 44) and Sardo Modicana (SM, n = 26) cattle breeds, reared in the island of Sardinia (Italy). A portion of the mitochondrial DNA hypervariable region was sequenced, in order to identify a potential signature of African introgression. The FST coefficients among populations ranged between 0.056 for Sa vs SB and 0.167 for SB vs SM. AMOVA analysis indicated there was a significant differentiation of the three breeds, although most of diversity was gathered at the within-breed level. The Median Joining Network of the Sardinian sequences showed a potential founder effect signature. A MJ network including Sardinian cattle plus African, Italian, Iberian and Asian sequences, revealed the presence of haplogroup T3, already detected in Sa cattle, and the presence of Hg T1 and Hg T1'2'3, in Sa and SB. The presence of a private haplotype belonging to haplogroup T1, which is characteristic of African taurine breeds, may be due to the introgression of Sardinian breeds with African cattle, either directly (most probable source: North African cattle) or indirectly (through a Mediterranean intermediary already introgressed with African blood).

A polymorphism at the melatonin receptor 1A (MTNR1A) gene in Sarda ewes affects fertility after AI in the spring.

The effect of MTNR1A gene polymorphisms on the fertility rate after AI in Sarda sheep was evaluated in 600 lactating adult ewes. Genomic DNA was subjected to amplification of the MTNR1A gene exon II. Amplicons were digested with restriction endonuclease MnlI. Ten samples from each genotype were sequenced. A polymorphism was detected (A612G) and ewes were determined to be +/+, +/- or -/- for the allele. Allelic frequency was 0.77 for the+allele and 0.23 for the - allele. The frequency of the +/+, +/- and -/- genotypes was 68, 19 and 13%, respectively. On 16 May 2009, 60 ewes from each genotype group were synchronised using intravaginal sponges containing 40 mg fluorogestone acetate for 14 days. At sponge removal, the ewes were administered 350 IU pregnant mare's serum gonadotropin and were then inseminated, 54-56 h later, with 400 × 10(6) spermatozoa. Pregnancies were confirmed 50 days after AI using transabdominal ultrasonography. Lambing dates and the number of newborn lambs were recorded within 155 days after AI. Conception and lambing rate were higher for ewes with the +/+ and +/- genotypes compared with those with the -/- genotype (P<0.01). In conclusion, there was a positive correlation between MTNR1A allele polymorphisms the reproductive response following synchronisation and AI in the spring.

Effects of different storage conditions, the farm and the stage of lactation on renneting parameters of goat milk investigated using the Formagraph method.

This study evaluated the effect of storage on renneting properties of goat milk investigated using the Formagraph method. Milk samples from 169 goats in three farms (F1, F2 and F3) were analysed during an entire lactation (45, 75, 105, 135 and 165 days in milking DIM), to obtain renneting parameters, both from fresh milk and after storage with Bronopol and freezing at -20°C and -80°C. As regards fresh milk, mean values of clotting time were between 12·51 (45 DIM) and 13·29 min (105 DIM and F2), the curd firming time between 1·77 (45 DIM) and 2·15 min (F1) and curd firmness between 42·09 (165 DIM) and 49·55 mm (45 DIM). No statistical difference was recorded after storage. After regression analysis, all prediction models showed significance value at P<0·001 with the highest R2 value for clotting time, 0·710 (fresh vs. frozen milk at -20°C), and the lowest for clot firmness, 0·281 (fresh vs. frozen milk at -80°C). Results demonstrated that assessment of goat milk coagulation properties using the Formagraph method is also achievable after freezing or Bronopol addition.