Detecting the footprint of selection on the genomes of Murciano-Granadina goats.

Artificial selection is one of the major forces modifying the genetic composition of livestock populations. Identifying genes under selection could be useful to elucidate their impact on phenotypic variation. We aimed to identify genomic regions targeted by selection for dairy and pigmentation traits in Murciano-Granadina goats. Performance of a selection scan based on the integrated haplotype score test in a population of 1183 Murciano-Granadina goats resulted in the identification of 77 candidate genomic regions/SNPs. The most significant selective sweeps mapped to chromosomes 1 (69.86 Mb), 4 (41.80-49.95 Mb), 11 (65.74 Mb), 12 (31.24 and 52.51 Mb), 17 (34.76-37.67 Mb), 22 (31.75 Mb), and 26 (26.69-31.05 Mb). By using previously generated RNA-Seq data, we built a catalogue of 6414 genes that are differentially expressed across goat lactation (i.e. 78 days post-partum, early lactation; 216 days post-partum, late lactation; 285 days post-partum, dry period). Interestingly, 183 of these genes mapped to selective sweeps and several of them display functions related with lipid, protein, and carbohydrate metabolism, insulin signaling, cell proliferation, as well as mammary development and involution. Of particular interest are the CSN3 and CSN1S2 genes, which encode two major milk proteins. Additionally, we found three pigmentation genes (GLI3, MC1R, and MITF) co-localizing with selective sweeps. Performance of a genome-wide association study and Sanger sequencing and TaqMan genotyping experiments revealed that the c.801C>G (p.Cys267Trp) polymorphism in the melanocortin 1 receptor (MC1R) gene is the main determinant of the black (GG or GC genotypes) and brown (CC genotypes) colorations of Murciano-Granadina goats.

Nonparametric analysis of casein complex genes' epistasis and their effects on phenotypic expression of milk yield and composition in Murciano-Granadina goats.

Improving knowledge on the causative polymorphisms or genes regulating the expression of milk quantitative and qualitative traits and their interconnections plays a major role in dairy goat breeding programs and genomic research. This information enables optimization of predictive and selective tools, to obtain better-performing animals to help satisfy market demands more efficiently. Goat milk casein proteins (αS1, αS2, ß, and κ) are encoded by 4 loci (CSN1S1, CSN1S2, CSN2, and CSN3) clustered within 250 kb on chromosome 6. Among the statistical methods used to identify epistatic interactions in genome-wide qualitative association studies (GWAS), gene-based methods have recently grown in popularity due to their better statistical power and biological interpretability. However, most of these methods make strong assumptions about the magnitude of the relationships between SNP and phenotype, limiting statistical power. Thus, the aims of this study were to quantify the epistatic relationships among 48 SNP in the casein complex on the expression of milk yield and components (fat, protein, dry matter, lactose, and somatic cells) in Murciano-Granadina goats, to explain the qualitative nature of the SNP used to quantify the genotypes produced as a result. Categorical principal component analysis (CATPCA) was used to delimit and group the number of SNP studied depending on their implications in the explanation of milk yield and components variability. Afterward, nonlinear canonical correlation analysis was used to identify relationships among and within the SNP groups detected by CATPCA. Our results suggest that 79.65% of variability in the traits evaluated may be ascribed to the epistatic relationships across and within 7 SNP groups. Two partially overlapping groups of epistatically interrelated SNP were detected: one group of 21 SNP, explaining 57.56% of variability, and another group of 20 SNP, explaining 42.43% (multiple fit ≥ 0.1). Additionally, SNP18, 32, and 36 (CSN1S2, CSN1S1, and CSN2 loci, respectively) were the most significant SNP to explain intragroup epistatic variability (component loading > |0.5|). Conclusively, milk yield and quality may not only depend on the specific casein gene pool of individuals, but may also be relevantly conditioned by the relationships set across and within such genes. Hence, studying epistasis in isolation may be crucial to optimize selective practices for economically important dairy traits.

Modelling for the inheritance of multiple births and fertility in endangered equids: Determining risk factors and genetic parameters in donkeys (Equus asinus).

Multiple births or twinning in equids are dangerous, undesirable situations that compromise the life of the dam and resulting offspring. However, embryo vitrification and freezing techniques take advantage of individuals whose multiple ovulations allow flushing more fertilised embryos from the oviduct to be collected, increasing the productivity and profitability of reproductive techniques. Embryo preservation is especially important in highly endangered populations such as certain donkey (Equus asinus) breeds; for which conventional reproductive techniques have previously been deemed inefficient. For instance, becoming an effective alternative to artificial insemination with frozen semen to preserve the individuals' genetic material. The objective of this study was to examine the historical foaling records of Andalusian donkeys to estimate prevalence, risk factors, phenotypic and genetic parameters for multiple births, assessing the cumulative foal number born per animal, maximum foal number per birth and multiple birth number per animal. We designed a Bayesian General Animal Mixed Model with single records considering the 'fixed' effects of birth year, birth season, birth month, sex, farm, location, and husbandry system. Age was considered and included as a linear and quadratic covariate. Gibbs sampling reported heritability estimates ranging from 0.18 ±â€¯0.101 to 0.24 ±â€¯0.078. Genetic and phenotypic correlations ranged from 0.496 ±â€¯0.298 to 0.846 ±â€¯0.152 and 0.206 ±â€¯0.063 to 0.607 ±â€¯0.054, respectively. Predicted breeding values obtained enable the potential selection against/for these traits, offering a new perspective for donkey breeding and conservation.

Low genome-wide homozygosity in 11 Spanish ovine breeds.

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (FROH ) are generally small (FROH  < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.

Effect of different olive oil-derived antioxidants (hydroxytyrosol and 3,4-dihydroxyphenylglycol) on the quality of frozen-thawed ram sperm.

The aim of the present study was to evaluate the effect of the addition of different concentrations of two olive oil-derived antioxidants, hydroxytyrosol (3,4-dihydroxyphenylethanol, HT) and 3,4-dihydroxyphenylglycol (DHPG), on ovine semen during the freezing-thawing process. Sperm was collected, pooled and diluted with commercial extenders and then divided into aliquots supplemented with different concentrations (10 µg/ml, 30 µg/ml, 50 µg/ml and 70 µg/ml) of HT, DHPG and a mixture (MIX) of both antioxidants. A control group, without antioxidant, was also prepared. Sperm motility, viability, acrosome integrity, mitochondrial membrane potential and lipid peroxidation (LPO) were assessed. The results showed that frozen-thawed ram spermatozoa exhibited lower values for motility, membrane integrity, acrosome and mitochondrial membrane potential than fresh samples (P ≤ 0.01). However, when antioxidants were added, thawed spermatozoa exhibited relatively low LPO, recording values similar to fresh spermatozoa; by contrast, the control group of frozen-thawed spermatozoa without antioxidants exhibited significantly higher LPO (P ≤ 0.01). The addition of a HT+DHPG mixture (MIX) had a negative impact on sperm membrane and acrosome integrity, suggesting that a pure antioxidant supplementation has the potential to offer superior results. In conclusion, HT and DHPG exhibited a positive effect on the frozen-thawed spermatozoa inasmuch as they reduced the LPO. These olive oil-derived antioxidants have the potential to improve frozen-thawed sperm quality, although further studies should be carried out to analyse the antioxidant effect at different times after thawing.

Expression patterns and genetic variation of the ovine skeletal muscle transcriptome of sheep from five Spanish meat breeds.

The goal of the current study is to analyse the gene expression profile of the ovine skeletal muscle as well as to characterize the genetic variation of transcripts expressed in such tissue. This aim has been achieved by sequencing the longissimus dorsi transcriptomes of 50 sheep distributed in five pools representing the Canaria de Pelo, Roja Mallorquina, Gallega, Xisqueta and Ripollesa Spanish autochthonous breeds. Approximately, 363 million reads per pool have been produced and 71.9-82.9% have been successfully mapped to the ovine genome in a paired-end mode (2 × 75 bp). The 200 most expressed muscle transcripts (≈1% of the total transcript count) account for 51% (Canaria de Pelo) to 67% (Gallega) of the total ovine skeletal muscle mRNA expression. These highly expressed genes play key roles in pathways related with striated muscle contraction, gluconeogenesis, glycolysis, citric acid cycle and respiratory electron transport. RNA-Sequencing of muscle transcripts has also revealed that ~72% of the SNPs detected with this approach are shared by at least two pools, and 10% of them segregate in the five pools under analysis. Most of the substitutions detected by RNA-Seq are synonymous or missense and only a minority are predicted to have consequences on protein function.

Dissection of ancestral genetic contributions to Creole goat populations.

Goats have played a key role as source of nourishment for humans in their expansion all over the world in long land and sea trips. This has guaranteed a place for this species in the important and rapid episode of livestock expansion triggered by Columbus' arrival in the Americas in the late 1400s. The aims of this study are to provide a comprehensive perspective on genetic diversity in American goat populations and to assess their origins and evolutionary trajectories. This was achieved by combining data from autosomal neutral genetic markers obtained in more than two thousand samples that encompass a wide range of Iberian, African and Creole goat breeds. In general, even though Creole populations differ clearly from each other, they lack a strong geographical pattern of differentiation, such that populations of different admixed ancestry share relatively close locations throughout the large geographical range included in this study. Important Iberian signatures were detected in most Creole populations studied, and many of them, particularly the Cuban Creole, also revealed an important contribution of African breeds. On the other hand, the Brazilian breeds showed a particular genetic structure and were clearly separated from the other Creole populations, with some influence from Cape Verde goats. These results provide a comprehensive characterisation of the present structure of goat genetic diversity, and a dissection of the Iberian and African influences that gave origin to different Creole caprine breeds, disentangling an important part of their evolutionary history. Creole breeds constitute an important reservoir of genetic diversity that justifies the development of appropriate management systems aimed at improving performance without loss of genomic diversity.

Storage temperature and sucrose concentrations affect ram sperm quality after vitrification.

The aim of this study was to analyse the characteristics of ram spermatozoa subjected to varying concentrations of sucrose, and the influence of storage temperature (22°C or 5°C) prior to vitrification. Ejaculated semen was diluted in TCFEY (tris-citric acid-fructose 20% egg yolk), and two aliquots were prepared at a final concentration of 100×106spz/ml, one maintained at room temperature (22°C) and the other at 5°C. In the first experiment, the toxicity of sucrose diluents on the sperm was analysed; sperm samples at different temperatures were diluted (1:2) in TCF-BSA 2% (control) or in the same extender supplemented with various sucrose concentrations (0.4M, 0.6M and 0.8M). The effects of vitrification were studied in the second experiment, where sperm samples were mixed with different concentrations of cryoprotectants (sucrose) and vitrified by being plunged directly into liquid nitrogen. In both experiments, the sperm quality was assessed by measuring motility, morphology, membrane functionality (HOST), viability, acrosome integrity and DNA fragmentation. The toxicity test revealed significant differences (p≤0.05) when different sucrose concentrations were used; lower total and progressive motility, normal morphology and membrane functionality were noted when sucrose concentration was higher, compared to the control treatment. Samples maintained at room temperature showed significantly (p≤0.05) higher viability than samples stored at 5°C. In contrast, although the quality of vitrified sperm was drastically decreased in comparison with fresh sperm, sucrose was associated with greater total motility, viability and membrane functionality. This improvement was closely linked to the temperature at which the sperm had been previously maintained, showing higher values when sperm was stored at 5°C. The main conclusions to be drawn from the study are therefore that sucrose shows promising potential as a cryoprotectant, and storing samples at 5°C is linked to improved sperm quality following vitrification.

A model to infer the demographic structure evolution of endangered donkey populations.

Stemming from The Worldwide Donkey Breeds Project, an initiative aiming at connecting international researchers and entities working with the donkey species, molecularly tested pedigree analyses were carried out to study the genetic diversity, structure and historical evolution of the Andalusian donkey breed since the 1980s to infer a model to study the situation of international endangered donkey breeds under the remarkably frequent unknown genetical background status behind them. Demographic and genetic variability parameters were evaluated using ENDOG (v4.8). Pedigree completeness and generation length were quantified for the four gametic pathways. Despite mean inbreeding was low, highly inbred animals were present in the pedigree. Average coancestry, relatedness, and non-random mating degree trends were computed. The effective population size based on individual inbreeding rate was about half when based on individual coancestry rate. Nei's distances and equivalent subpopulations number indicated differentiated farms in a highly structured population. Although genetic diversity loss since the founder generations could be considered small, intraherd breeding policies and the excessive contribution of few ancestors to the gene pool could lead to narrower pedigree bottlenecks. Long average generation intervals could be considered when reducing inbreeding. Wright's fixation statistics indicated slight inbreeding between farms. Pedigree shallowness suggested applying new breeding strategies to reliably estimate descriptive parameters and control the negative effects of inbreeding, which could indeed, mean the key to preserve such valuable animal resources avoiding the extinction they potentially head towards, making the present model become an international referent when assessing endangered donkey populations.

The legacy of Columbus in American horse populations assessed by microsatellite markers.

Criollo horse populations descend from horses brought from the Iberian Peninsula over the period of colonization (15th to 17th century). They are spread throughout the Americas and have potentially undergone genetic hybridization with other breeds in the recent past. In this study, 25 autosomal microsatellites were genotyped in 50 horse breeds representing Criollo populations from 12 American countries (27 breeds), breeds from the Iberian Peninsula (19), one breed each from France and Morocco and two cosmopolitan horse breeds (Thoroughbred and Arabian). The genetic relationships among breeds identified five clusters: Celtic; Iberian; North American with Thoroughbred influence; most Colombian breeds; and nearly all other Criollo breeds. The group of "all other Criollo breeds" had the closest genetic relationship with breeds originating from the Iberian Peninsula, specifically with the Celtic group. For the whole set of Criollo breeds analysed, the estimated genetic contribution from other breeds was approximately 50%, 30% and 20% for the Celtic, Iberian and Arab-Thoroughbred groups, respectively. The spatial distribution of genetic diversity indicates that hotspots of genetic diversity are observed in populations from Colombia, Ecuador, Brazil, Paraguay and western United States, possibly indicating points of arrival and dispersion of Criollo horses in the American continent. These results indicate that Criollo breeds share a common ancestry, but that each breed has its own identity.

Population structure of eleven Spanish ovine breeds and detection of selective sweeps with BayeScan and hapFLK.

The goals of the current work were to analyse the population structure of 11 Spanish ovine breeds and to detect genomic regions that may have been targeted by selection. A total of 141 individuals were genotyped with the Infinium 50 K Ovine SNP BeadChip (Illumina). We combined this dataset with Spanish ovine data previously reported by the International Sheep Genomics Consortium (N = 229). Multidimensional scaling and Admixture analyses revealed that Canaria de Pelo and, to a lesser extent, Roja Mallorquina, Latxa and Churra are clearly differentiated populations, while the remaining seven breeds (Ojalada, Castellana, Gallega, Xisqueta, Ripollesa, Rasa Aragonesa and Segureña) share a similar genetic background. Performance of a genome scan with BayeScan and hapFLK allowed us identifying three genomic regions that are consistently detected with both methods i.e. Oar3 (150-154 Mb), Oar6 (4-49 Mb) and Oar13 (68-74 Mb). Neighbor-joining trees based on polymorphisms mapping to these three selective sweeps did not show a clustering of breeds according to their predominant productive specialization (except the local tree based on Oar13 SNPs). Such cryptic signatures of selection have been also found in the bovine genome, posing a considerable challenge to understand the biological consequences of artificial selection.

Conservation priorities of Iberoamerican pig breeds and their ancestors based on microsatellite information.

Criollo pig breeds are descendants from pigs brought to the American continent starting with Columbus second trip in 1493. Pigs currently play a key role in social economy and community cultural identity in Latin America. The aim of this study was to establish conservation priorities among a comprehensive group of Criollo pig breeds based on a set of 24 microsatellite markers and using different criteria. Spain and Portugal pig breeds, wild boar populations of different European geographic origins and commercial pig breeds were included in the analysis as potential genetic influences in the development of Criollo pig breeds. Different methods, differing in the weight given to within- and between-breed genetic variability, were used in order to estimate the contribution of each breed to global genetic diversity. As expected, the partial contribution to total heterozygosity gave high priority to Criollo pig breeds, whereas Weitzman procedures prioritized Iberian Peninsula breeds. With the combined within- and between-breed approaches, different conservation priorities were achieved. The Core Set methodologies highly prioritized Criollo pig breeds (Cr. Boliviano, Cr. Pacifico, Cr. Cubano and Cr. Guadalupe). However, weighing the between- and within-breed components with FST and 1-FST, respectively, resulted in higher contributions of Iberian breeds. In spite of the different conservation priorities according to the methodology used, other factors in addition to genetic information also need to be considered in conservation programmes, such as the economic, cultural or historical value of the breeds involved.

Genetic relationships among American donkey populations: insights into the process of colonization.

This study presents the first insights into the genetic diversity and structure of the American donkey metapopulation. The primary objectives were to detect the main structural features underlying variability among American donkey populations, identify boundaries between differentiated gene pools, and draw the main colonization pathways since the introduction of donkeys into America in the 15th century. A panel of 14 microsatellite markers was applied for genotyping 350 American donkeys from 13 countries. The genetic structure of this metapopulation was analysed using descriptive statistics and Bayesian model-based methods. These populations were then compared to a database containing information on 476 individuals from 11 European breeds to identify the most likely ancestral donor populations. Results showed the presence of two distinct genetic pools, with confluence of the two in Colombia. The southern pool showed a unique genetic signature subsequent to an older founder event, but lacked any significant influence of modern gene flow from Europe. The northern pool, conversely, may have retained more ancestral polymorphisms and/or have experienced modern gene flow from Spanish breeds. The Andalusian and, to a lesser extent, the Catalan breeds have left a more pronounced footprint in some of the American donkey populations analysed.

Genetic parameters of traits associated with the growth curve in Segureña sheep.

This paper studies the genetic importance of growth curve parameters and their relevance as selection criteria in breeding programmes of Segureño sheep. Logistic and Verhulst growth functions were chosen for their best fit to BW/age in this breed; the first showed the best general fit and the second the best individual fit. Live weights of 41 330 individuals from the historical archives of the National Association of Segureña Sheep Breeders were used in the analysis. The progeny of 1464 rams and 27 048 ewes were used to study the genetic and phenotypic parameters of growth curve parameters and derived traits. Reproductive management in the population consists in controlled natural mating inside every herd, with a minimum of 15% of the females fertilized by artificial insemination with fresh semen; with the purpose being the herd genetic connections, all herd genealogies are screened with DNA markers. Estimates of growth curve parameters from birth to 80 days were obtained for each individual and each function by the non-linear regression procedure using IBM SPSS statistics (version 21) with the Levenberg-Marquart estimation method. (Co)variance components and genetic parameters were estimated by using the REML/Animal model methodology. The heritability of mature weight was estimated as 0.41±0.042 and 0.38±0.021 with the logistic and Verhulst models, respectively, and the heritability of other parameters ranged from 0.41 to 0.62 and 0.37 to 0.61, with the models, respectively. A negative genetic correlation between mature weight and rate of maturing was found.

A mitochondrial analysis reveals distinct founder effect signatures in Canarian and Balearic goats.

In the course of human migrations, domestic animals often have been translocated to islands with the aim of assuring food availability. These founder events are expected to leave a genetic footprint that may be recognised nowadays. Herewith, we have examined the mitochondrial diversity of goat populations living in the Canarian and Balearic archipelagos. Median-joining network analysis produced very distinct network topologies for these two populations. Indeed, a majority of Canarian goats shared a single ancestral haplotype that segregated in all sampled islands, suggesting a single founder effect followed by a stepping-stone pattern of diffusion. This haplotype also was present in samples collected from archaeological assemblies at Gran Canaria and Lanzarote, making evident its widespread distribution in ancient times. In stark contrast, goats from Majorca and Ibiza did not share any mitochondrial haplotypes, indicating the occurrence of two independent founder events. Furthermore, in Majorcan goats, we detected the segregation of the mitochondrial G haplogroup that has only been identified in goats from Egypt, Iran and Turkey. This finding suggests the translocation of Asian and/or African goats to Majorca, possibly as a consequence of the Phoenician and Carthaginian colonisations of this island.

Genetic characterization of Uruguayan Pampa Rocha pigs with microsatellite markers.

In this study, we genetically characterized the Uruguayan pig breed Pampa Rocha. Genetic variability was assessed by analyzing a panel of 25 microsatellite markers from a sample of 39 individuals. Pampa Rocha pigs showed high genetic variability with observed and expected heterozygosities of 0.583 and 0.603, respectively. The mean number of alleles was 5.72. Twenty-four markers were polymorphic, with 95.8% of them in Hardy Weinberg equilibrium. The level of endogamy was low (FIS = 0.0475). A factorial analysis of correspondence was used to assess the genetic differences between Pampa Rocha and other pig breeds; genetic distances were calculated, and a tree was designed to reflect the distance matrix. Individuals were also allocated into clusters. This analysis showed that the Pampa Rocha breed was separated from the other breeds along the first and second axes. The neighbour-joining tree generated by the genetic distances DA showed clustering of Pampa Rocha with the Meishan breed. The allocation of individuals to clusters showed a clear separation of Pampa Rocha pigs. These results provide insights into the genetic variability of Pampa Rocha pigs and indicate that this breed is a well-defined genetic entity.

Characterization of commercial and biological growth curves in the Segureña sheep breed.

Non-linear models were analysed to describe both the biological and commercial growth curves of the Segureña sheep, one of the most important Spanish breeds. We evaluated Brody, von Bertalanffy, Verhulst, logistic and Gompertz models, using historical data from the National Association of Segureña Sheep Breeders (ANCOS). These records were collected between 2000 and 2013, from a total of 129 610 weight observations ranging from birth to adulthood. The aim of this research was to establish the mathematical behaviour of body development throughout this breed's commercial life (birth to slaughter) and biological life (birth to adulthood); comparison between both slopes gives important information regarding the best time for slaughter, informs dietary advice according to animals' needs, permits economical predictions of productions and, by using the curve parameters as selection criteria, enables improvements in growth characteristics of the breed. Models were fitted according to the non-linear regression procedure of statistical package SPSS version19. Model parameters were estimated using the Levenberg-Marquardt algorithm. Candidate models were compared using the determinative coefficient, mean square error, number of iterations, Akaike information coefficient and biological coherence of the estimated parameters. The von Bertalanffy and logistic models were found to be best suited to the biological and commercial growth curves, respectively, for both sexes. The Brody equation was found to be unsuitable for studying the commercial growth curve. Differences between the parameters in both sexes indicate a strong impact of sexual dimorphism on growth. This can emphasize the value of the highest growth rate for females, indicating that they reach maturity earlier.

The biodiversity and genetic structure of Balearic sheep breeds.

The Balearic sheep breeds, Mallorquina, Menorquina, Roja Mallorquina, Ibicenca and one possible new genetic group, Formentera, constitute a unique genetic resource in the Mediterranean farming landscape, displaying high genetic diversity levels and being well differentiated among themselves and with respect to the continental sheep breeds. We used a microsatellite panel of markers to study genetic diversity and relationships with other Spanish breeds. The results reported in this study have important implications for the use, conservation and breeding of Balearic sheep stocks. A mean number of 7.59 alleles was found among the Balearic sheep breeds for the microsatellites scored. The whole mean value of observed heterozygosity amounted to 0.62, whereas the expected heterozygosity value was 0.69, suggesting the presence of a great degree of genetic variability, although a significant deficit of heterozygotes was detected for some markers. Genetic distance estimates showed that Balearic sheep are differentiated from the other Spanish breeds and in particular, from the Merino type. The Ibicenca breed showed the highest distance value from other breeds. The neighbour-net method of analysis clustered the Roja Mallorquina, Menorquina and Mallorquina breeds. The Structure results clearly demonstrated the genetic differentiation among the four Balearic sheep breeds, with the Ibicenca and Formentera races joined, with slight migration among them. Few external genetic influences from the Spanish mainland breeds were detected.

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.

Genetic diversity analysis of the Uruguayan Creole cattle breed using microsatellites and mtDNA markers.

The Uruguayan Creole cattle population (N = 600) is located in a native habitat in south-east Uruguay. We analyzed its genetic diversity and compared it to other populations of American Creole cattle. A random sample of 64 animals was genotyped for a set of 17 microsatellite loci, and the D-loop hyper-variable region of mtDNA was sequenced for 28 calves of the same generation. We identified an average of 5.59 alleles per locus, with expected heterozygosities between 0.466 and 0.850 and an expected mean heterozygosity of 0.664. The polymorphic information content ranged from 0.360 to 0.820, and the global FIS index was 0.037. The D-loop analysis revealed three haplotypes (UY1, UY2 and UY3), belonging to the European matriline group, with a haplotype diversity of 0.532. The history of the population, changes in the effective population size, bottlenecks, and genetic drift are possible causes of the genetic variability patterns that we detected.