This work aimed to use 16S ribosomal RNA sequencing with the Illumina MiSeq platform to describe the milk microbiota from 50 healthy Assaf ewes. The global observed microbial community for clinically healthy milk samples analysed was complex and showed a vast diversity. The core microbiota of the sheep milk includes five genera: Staphylococcus, Lactobacillus, Corynebacterium, Streptococcus and Escherichia/Shigella. Although there are some differences, some of these genera are common with the microbiota core pattern of milk from other species, especially with dairy cows. The microbial composition of the studied samples, based on the definition of amplicon sequence variants, was analysed through a correlation network. A preliminary analysis by grouping the milk samples based on their somatic cell count (SCC), which is considered an indicator of subclinical mastitis (SM), showed certain differences for the core of the samples identified as SM. The differences in the microbiota diversity pattern among samples might also suggest that subclinical mastitis would be associated with the significant increase in some genera that are inhabitants of the mammary gland and a remarkable concomitant reduction in the microbial diversity. Additionally, we have also presented here a preliminary analysis to assess the impact of the sheep milk microbiome on SCC, as an indicator of subclinical mastitis. The results here reported provide a first characterization of the sheep milk microbiota and settle the basis for future studies in this field.

The global production of sheep milk is growing, and the main industrial use of sheep milk is cheese making. The Spanish Churra sheep breed is one of the most important native dairy breeds in Spain. The present study aimed to estimate genetic parameters for a wide range of traits influencing the cheese-making ability of Churra sheep milk. Using a total of 1,049 Churra ewes, we studied the following cheese-making traits: 4 traits related to milk coagulation properties (rennet coagulation time, curd-firming time, and curd firmness at 30 and 60 min after addition of rennet), 2 traits related to cheese yield (individual laboratory cheese yield and individual laboratory dried curd yield), and 3 traits measuring curd firmness over time (maximum curd firmness, time to attain maximum curd firmness, and syneresis). In addition, a list of milk traits, including the native pH of the milk and several milk production and composition traits (milk yield; the fat, protein, and dried extract percentages; and the somatic cell count), were also analyzed for the studied animals. After discarding the noncoagulating samples (only 3.7%), data of 1,010 ewes were analyzed with multiple-trait animal models by using the restricted maximum likelihood method to estimate (co)variance components, heritabilities, and genetic correlations. In general, the heritability estimates were low to moderate, ranging from 0.08 (for the individual laboratory dried curd yield trait) to 0.42 (for the fat percentage trait). High genetic correlations were found within pairs of related traits (i.e., 0.93 between fat and dried extract percentages, −0.93 between the log of the curd-firming time and curd firmness at 30 min, 0.70 between individual laboratory cheese yield and individual laboratory dried curd yield, and −0.94 between time to attain maximum curd firmness and syneresis). Considering all the information provided here, we suggest that in addition to the current consideration of the protein percentage trait for improving cheese yield traits, the inclusion of the pH of milk as a measured trait in the Churra dairy breeding program would represent an efficient strategy for improving the cheese-making ability of milk from this breed.

This study aimed to perform a GWAS to identify genomic regions associated with milk and cheese-making traits in Assaf and Churra dairy sheep breeds; second, it aimed to identify possible positional and functional candidate genes and their interactions through post-GWAS studies. For 2,020 dairy ewes from 2 breeds (1,039 Spanish Assaf and 981 Churra), milk samples were collected and analyzed to determine 6 milk production and composition traits and 6 traits related to milk coagulation properties and cheese yield. The genetic profiles of the ewes were obtained using a genotyping chip array that included 50,934 SNP markers. For both milk and cheese-making traits, separate single-breed GWAS were performed using GCTA software. The set of positional candidate genes identified via GWAS was subjected to guilt-by-association-based prioritization analysis with ToppGene software. Totals of 84 and 139 chromosome-wise significant associations for the 6 milk traits and the 6 cheese-making traits were identified in this study. No significant SNPs were found in common between the 2 studied breeds, possibly due to their genetic heterogeneity of the phenotypes under study. Additionally, 63 and 176 positional candidate genes were located in the genomic intervals defined as confidence regions in relation to the significant SNPs identified for the analyzed traits for Assaf and Churra breeds. After the functional prioritization analysis, 71 genes were identified as promising positional and functional candidate genes and proposed as targets of future research to identify putative causative variants in relation to the traits under examination. In addition, this multitrait study allowed us to identify variants that have a pleiotropic effect on both milk production and cheese-related traits. The incorporation of variants among the proposed functional and positional candidate genes into genomic selection strategies represent an interesting approach for achieving rapid genetic gains, specifically for those traits difficult to measure, such as cheese-making traits.

Transitioning from traditional to new genotyping technologies requires the development of bridging methodologies to avoid extra genotyping costs. This study aims to identify the optimum number of single nucleotide polymorphisms (SNPs) necessary to accurately impute microsatellite markers to develop a low-density SNP chip for parentage verification in the Assaf sheep breed. The accuracy of microsatellite marker imputation was assessed with three metrics: genotype concordance (C), genotype dosage (length r2), and allelic dosage (allelic r2), for all imputation scenarios tested (0.5–10 Mb microsatellite flanking SNP windows). The imputation accuracy for the three metrics analyzed for all haplotype lengths tested was higher than 0.90 (C), 0.80 (length r2), and 0.75 (allelic r2), indicating strong genotype concordance. The window with 2 Mb length provides the best accuracy for the imputation procedure and the design of an affordable low-density SNP chip for parentage testing. We additionally evaluated imputation performance under two null models, naive (imputing the most common allele) and random (imputing by randomly selecting the allele), which in comparison showed weak genotype concordances (0.41 and 0.15, respectively). Therefore, we describe a precise methodology in the present article to impute multiallelic microsatellite genotypes from a low-density SNP chip in sheep and solve the problem of parentage verification when different genotyping platforms have been used across generations.

Different studies have shown that polymorphisms in the sequence of genes coding for the milk proteins and milk fatty acids are associated with milk composition traits as well as with cheese-making traits. However, the lack of coincident results across sheep populations has prevented the use of this information in sheep breeding programs. The main objective of this study was to exploit the information derived from a total of 175 whole genome resequencing (WGR) datasets from 43 domestic sheep breeds and three wild sheep to evaluate the genetic diversity of 24 candidate genes for milk composition and identify genetic variants with a potential phenotypic effect. The functional annotation of the identified variants highlighted five single nucleotide polymorphisms (SNPs) predicted to have a high impact on the protein function and 42 missense SNPs with a putative deleterious effect. When comparing the allelic frequencies at these 47 polymorphisms with relevant functional effects between the genomes of Assaf and Churra sheep breeds, two missense deleterious variants were identified as potential markers associated to the milk composition differences found between the Churra and Assaf: XDH:92215727C>T and LALBA:137390760T>C. Future research is required to confirm the effect of the potential functionally relevant variants identified in the present study on milk composition and cheese-making traits.