O. Vidal

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.

Genetics of serum and muscle lipids in pigs

Pork meat is one of the most important sources of animal protein in the human diet. Its nutritional properties are partly determined by intramuscular fat content and composition, with existing general consensus about the detrimental effects of cholesterol and saturated fat on cardiovascular health in humans. Because of their physiological resemblance, pigs can be also used as a valuable animal model to study the genetics of human diseases such as atherosclerosis, obesity and dyslipidaemias. Heritability estimates and QTL maps of porcine muscle and serum lipid traits evidence that a considerable amount of genetic variance determining these phenotypes exists, but its molecular basis remains mostly unknown. The recent advent of high-throughput genotyping and sequencing technologies has revolutionised the field of animal genomics. With these powerful tools, finding needles in the genomic haystack has become increasingly feasible. However, these methodological advances should not be deemed as magic bullets. The goal of identifying the many polymorphisms that shape the variability of lipid phenotypes is so challenging that success can be achieved only under the scope of large international consortia.

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.

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.

Nucleotide Sequence and Polymorphism of the Pig Acyl Coenzyme A Synthetase Long-Chain 1 (ACSL1) Gene

We have sequenced 3,013 bp of the pig acyl coenzyme A long-chain synthetase 1 (ACSL1) gene. Sequence analysis allowed us to identify three conserved elements in the predicted amino acid sequence, two of them related to the ATP/AMP signature motif and the third involved in enzyme catalysis and fatty acid substrate specificity. In addition, we have identified five C → T and one G → A transition polymorphisms located in exon 16 (SNPe16), exon 17 (SNPe17) and the 3′ UTR (SNPa to d), which have been genotyped in 143 pigs from the Landrace, Large White, Piétrain, Iberian, and Duroc breeds.

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.

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.

Conservation of Goat Populations from Southwestern Europe Based on Molecular Diversity Criteria

Goat farming plays a key role in agricultural activity and in maintaining forest lands in Southwestern Europe. Remarkably, the Iberian Peninsula represents nearly 25% of the European goat census. Goat husbandry is often associated with low input production systems and uses selective breeding programs, which are less advanced than those employed in other livestock. Native goat breeds are very well adapted to produce in marginal areas under extensive conditions. Loss of their genetic diversity could have important economic, ecological and scientific implications as well as social consequences. Several methodologies have been developed to preserve the genetic diversity of single populations, but additional problems arise when a group of breeds, i.e., subpopulations, is considered in conservation programs. The conservation priority of a breed depends on its contribution to the overall genetic diversity of the species, in terms of the intrinsic genetic variation that it harbors and also of its relationship with other breeds. However, the estimation of the contributions of each of these two components to overall genetic diversity cannot be easily assessed. Besides, conservation goals in the short-term (avoidance of inbreeding) and long-term (adaptation to future environmental changes) should be considered when taking conservative decisions. A comprehensive analysis of Iberian goat breeds has been carried out to evaluate conservation priorities based on methodologies that account for within- or between-breed genetic diversity, or combinations of both. Based on genetic distinctiveness, breeds such as Palmera, Formentera, and Blanca Celtiberica were prioritized, whereas the maximum priority was assigned to Florida, Pirenaica, Retinta, and Moncaina breeds when focusing on within-breed diversity. Overall, combined approaches showed very little variation among breeds reflecting a history of extensive gene flow, partly due to transhumance and recent divergence. The main conclusion of our study is that these statistical analyses are useful, but conservation decisions must take into account other factors in addition to strict genetic diversity classification.