Next Article in Journal
Epidemiology, Genetic Characterization, and Pathogenesis of Avian Influenza H5N8 Viruses Circulating in Northern and Southern Parts of Egypt, 2017–2019
Next Article in Special Issue
Comprehensive Proteomic Analysis Revealed a Large Number of Newly Identified Proteins in the Small Extracellular Vesicles of Milk from Late-Stage Lactating Cows
Previous Article in Journal
Development of Nanobodies Targeting Peste des Petits Ruminants Virus: The Prospect in Disease Diagnosis and Therapy
Previous Article in Special Issue
Physicochemical, Microbiological and Technological Properties of Red Deer (Cervus elaphus) Milk during Lactation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Genes Associated with Dairy Traits in Sarda Sheep

by
Michele Pazzola
1,
Giuseppe Massimo Vacca
1,*,
Pietro Paschino
1,
Giovanni Bittante
2 and
Maria Luisa Dettori
1
1
Department of Veterinary Medicine, University of Sassari, Via Vienna 2, 07100 Sassari, SS, Italy
2
Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale dell’Università 16, 35020 Legnaro, PD, Italy
*
Author to whom correspondence should be addressed.
Animals 2021, 11(8), 2207; https://doi.org/10.3390/ani11082207
Submission received: 21 June 2021 / Revised: 9 July 2021 / Accepted: 23 July 2021 / Published: 26 July 2021
(This article belongs to the Special Issue New Insights into the Milk)

Abstract

:

Simple Summary

In a population of 1112 individuals of Sarda sheep, we investigated the polymorphism of 45 single nucleotide polymorphisms (SNPs) in 15 different genes that have been previously investigated because of their relation to metabolism and innate immunity. For SNPs with a sufficient allele frequency, we also tested their influence on a wide panel of dairy traits. The genotyping evidenced some new genetic patterns, and seven SNPs had significant influences on the milk yield, composition and coagulation traits. The results may represent a starting point to acquire new knowledge in the field of dairy sheep sector and to improve the breeding schemes for the Sarda sheep breed.

Abstract

The aim of the present research was to analyze the variability of 45 SNPs from different genes involved in metabolism and innate immunity to perform an association analysis with the milk yield, composition and milk coagulation traits. A population of 1112 Sarda breed sheep was sampled. Genotyping was generated by a TaqMan Open ArrayTM. Thirty out of the 45 SNPs were polymorphic, and 12 displayed a minor allele frequency higher than 0.05. An association analysis showed that the variability at genes PRKAG3 and CD14 was significantly associated with the daily milk yield. The variability at PRKAG3 was also associated with the protein and casein content, somatic cell score and bacterial score. The variation at the PRKAA2 gene was associated with the milk lactose concentration. The SNPs at CD14 were also associated with the traditional milk coagulation properties, while the SNPs at GHR and GHRHR were associated with kSR, a derived coagulation parameter related to the rate of syneresis. The information provided here is new and increases our knowledge of genotype–phenotype interactions in sheep. Our findings might be useful in appropriate breeding schemes to be set up for the Sarda sheep breed, but these should be confirmed by further studies, possibly performed on independent populations.

1. Introduction

Sheep are farmed all over the world, because of their great adaptability and relatively high productions. Productive specialization is different among countries and geographical areas. In the Middle East and in the Mediterranean basin, sheep milk is of major importance [1]. In the European Union, sheep milk production has been markedly increased in the last decades, extending from 2 to about 3 million t in the period of 1960–2019, and dairy sheep farming is a relevant sector of agriculture in many areas, such as for Greece, the leading EU country with 944 thousand t of ewe milk produced in 2019, followed by Romania, Spain and Italy [2].
In dairy sheep farming, almost all milk is transformed into cheese; therefore, milk traits related to cheese-making, including milk coagulation properties (MCP), curd firming and syneresis, need to be evaluated, as they might impact the farm income [3]. The MCP are described as the rennet coagulation time (RCT), curd firming time (k20) and curd firmness after 30 min of analysis (a30), obtained by the lactodynamographic instrument [4]. The same analysis can be extended to 60 min with the additional single point measures of curd firmness at 45 (a45) and 60 (a60) minutes [5]. The data can be further analyzed by utilizing all the lactodynamographic records (one curd-firming data every 15 s) to set up the modeling equations of the coagulation process [6]. The derived parameters estimated by the modeling process are useful for understanding the entire milk coagulation process: the asymptotic potential curd firmness, CFP, the curd-firming instant rate constant, kCF, the syneresis instant rate constant, kSR, the maximum curd firmness value, CFmax, and the time to attain the maximum curd firmness, tmax [6]. Traditional and derived MCP have been studied in cattle and, later, in sheep to provide important information on the dynamics of coagulation, highlighting the similarities and differences between species [6,7].
The molecular information underlying the biological processes has often been useful for understanding important phenotypes in sheep, such as prolificacy [8] or muscular hypertrophy [9]. The path to understand the molecular basis of the exceptional meatiness of the Texel sheep went through the identification of a quantitative trait locus on sheep chromosome 2 after a whole-genome scan and, finally, to the identification of a SNP in the 3′UTR of the GDF8 (Growth Differentiation Factor 8) gene, which causes translational inhibition of the myostatin gene in sheep, thus contributing to muscular hypertrophy [9].
An association analysis, by correlating markers with phenotypes across a population, in addition to genome-wide association studies (based on whole-genome genotyping), and family-based linkage studies, might be useful to understand the role played by candidate genes or haplotypes in a population [10,11]. In the literature, many association studies investigating the ovine species (Ovis aries) are focused on growth, body weight and fat traits, as well as fertility and wool traits [12]. The research articles investigating sheep dairy traits are lower in number, and those mainly analyze casein and whey protein genes, as casein genes encode 80% of milk proteins (caseins), while whey proteins account for 20% of total milk proteins [13,14]. Noce et al. [13] investigated the polymorphism at the promoter regions of the four casein genes, in addition to the LALBA (alpha lactalbumin) and BLG (beta lactoglobulin) genes, and they highlighted significant associations with milk composition and coagulation properties in Sarda sheep [13]. Usai et al. [14] performed OvineSNP50 genotyping to map the genomic regions affecting milk traits and obtained a strong signal for the protein content on OAR chromosome 6, overlapping the region including the casein gene cluster in addition to OAR 4 for the fat content and OAR 11 for the milk yield traits.
Additionally, lactation is a complex phenomenon involving metabolic changes determined by many other genes, which ultimately, might affect the milk traits. The study by Bittante et al. [15] evidenced that, in sheep species, milk coagulation traits have moderate-to-high values of heritability, which can be exploitable in breeding programs. Hence, we hypothesize a role of some genes on the milk coagulation traits, based on their biological role of avoiding those genes that have already been studied and whose association with dairy traits has been demonstrated [12,13,14].
As regards the genes involved in cellular metabolism, PRKAA2 (Protein Kinase AMP-Activated Catalytic Subunit Alpha2), PRKAG3 (Protein Kinase AMP-Activated Non-Catalytic Subunit Gamma3) and the genes linked to the growth hormone system, such as GHRHR (Growth Hormone Releasing Hormone Receptor), GHR (Growth Hormone receptor) and IGF1 (Insulin Like Growth Factor 1), the PRKAA2 gene encodes a protein involved in the monitoring of the cellular energy status and, also, in response to stress [16]. The expression of PRKAA2 changes throughout the different stages of pregnancy and lactation in cattle, suggesting a crucial role in the complex processes of lipolysis and milk synthesis [17,18]. PRKAG3 is mainly studied for its role in regulating energy homeostasis and, also, in response to environmental or nutritional stress, mainly in skeletal muscle cells, where PRKAG3 increases the oxidation of fatty acids and glucose uptake to meet the muscle energy demands [19]. The PRKAG3 protein has been shown to affect the glycogen content in muscle cells and then the meat quality in livestock [20,21,22]. The growth hormone system plays an important role in protein and lipid metabolism and affects the galactopoiesis in dairy species [23]. We investigated the polymorphism at the B4GALT1 (β1,4-galactosyltransferase 1) and ABCG2 (ATP-binding cassette subfamily G member 2) genes. The β1,4-galactosyltransferase 1 protein is involved in glucose metabolism, as it forms a heterodimer with alpha lactalbumin in mammary epithelial cells. That generates the lactose synthase enzyme, which synthesizes lactose by combining one molecule of UDP-galactose with one of glucose [24]. The ABCG2 protein belongs to a transporter family, including the proteins involved in the transport of molecules across cell membranes [25]. The ABCG2 gene is expressed in different tissues, including the mammary gland epithelial cell during lactation, and plays an important role in active secretion of molecules across the blood–milk barrier, including drugs and toxins [26].
As regards the genes involved in the immune response and defense against pathogens, CD14 (Cluster of Differentiation 14 Molecule), TLR2 (Toll-Like Receptor 2), TLR4 (Toll-Like Receptor 4), and PI (Serpin Peptidase Inhibitor), whose variability in dairy sheep and correlation with milk traits has not been investigated yet, the CD14 gene encodes a molecule that is fundamental for innate immunity and acts against a range of pathogens. The CD14 protein is a surface antigen that is preferentially expressed on monocytes and macrophages, which contributes to mediate the response to bacterial lipopolysaccharide [16]. The TLR4 gene encodes a protein, a member of the Toll-like receptor family, involved in pathogen recognition and innate immunity. As the TLR4 receptor is involved in the cellular response to the lipopolysaccharides found in most Gram-negative bacteria, it also plays a role in mastitis [27]. We also evaluated the variability at the Neuropeptide FF receptor 2 (NPFFR2) gene. NPFFR2 is a member of the G protein-coupled receptor superfamily, membrane proteins activated by neuropeptides NPAF and NPFF [28]. The NPFF neuropeptides were shown to display anti-inflammatory activity in a mouse model and were linked to T-lymphocyte proliferation and the immune response [29].
In addition to the above-mentioned metabolism and immune response-related genes, we investigated the variability of the PLCE1 and TYRP1 genes. The PLCE1 (Phospholipase C Epsilon 1) gene encodes a phospholipase enzyme, which catalyzes the hydrolysis of phosphatidylinositol-4,5-bisphosphate, and produces the second messengers inositol 1,4,5-triphosphate (IP3) and diacylglycerol (DAG). The two second messengers can affect cell growth, differentiation, and gene expression [16]. The TYRP1 gene (Tyrosinase Related Protein 1) plays an important role in the synthesis of eumelanin and is associated with the coat color in domestic animals, including sheep [30].
The present study aimed to explore, in a population of 1112 dairy ewes of the Sarda breed, the variability of 45 SNPs located on fifteen genes responsible for metabolism and immune defense reactions in order to perform an association analysis between the genotyped SNPs and a wide range of dairy traits, including milk yield and composition, as well as traditional and model coagulation traits.

2. Materials and Methods

2.1. Farms, Animals, and Samples

For the present study, 1112 ewes in 23 different farms were sampled. Detailed features of the animals and farms were reported by Pazzola et al. [5]. The main aspects were the following: ewes were of Sarda breed, regularly recorded in the flock book, and daughters of 120 different sires (minimum 4, maximum 40 daughters per sire); parity order and stage of lactation of ewes were, respectively, between the first and ninth, and from the second to the seventh month after parturition; farms were of commercial type, located in the regional territory of Sardinia (Italy) and managed under semi-extensive farming systems.
At each farm, a variable group of ewes were sampled (minimum 22, maximum 89 ewes per farm) in a single sampling day. Individual milk samples were collected in 200-mL sterile plastic containers during the evening milking, which was performed by using manually operated milking machines. The daily milk yield (dMY) was calculated as the sum of the morning and evening milking of the same sampling day. Blood samples were individually collected in K3EDTA vacuum tubes (BD Vacutainer, Becton Dickinson, Milan, Italy) from the jugular vein. Milk and blood samples were transported at 4 °C to the laboratories to be stored and analyzed.

2.2. DNA Extracton and Genotyping

Blood samples were processed for DNA extraction within 24 h after sampling. Genomic DNA was extracted using the Gentra Puregene blood kit (Qiagen, Hilden, Germany). DNA concentration and purity were assessed using an Eppendorf BioPhotometer (Eppendorf, Hamburg, Germany).
We choose each SNP to be analyzed in the present experiment, respecting the technical requirements of the TaqMan® assay, and by using the Ensembl Genome Browser (https://www.ensembl.org/index.html; accessed date 30 March 2021) to localize and choose each SNP. The list of selected SNPs is reported in Supplementary Table S1. The selected SNPs and their flanking sequences of 60 bp (Supplementary Table S1) were submitted to the Custom TaqMan Assay Design Tool website (https://www5.appliedbiosystems.com/tools/cadt; accessed date 30 March 2021; Life Technologies Thermo Fisher Scientific, Waltham, MA, USA) to ascertain their suitability to be genotyped with the TaqMan Open Array multiplex assay.
The 12 K Flex QuantStudio instrument (Thermo Fisher Scientific, Waltham, MA, USA) was used for genotyping and Taqman Genotyper v.1.3 software (Applied Biosystems, Waltham, MA, USA) for visualizing the genotypes. Automatic allele calling was carried out using GeneCall software (Illumina, San Diego, CA, USA) with a CG threshold of 0.25 and a further analysis of the individual genotypes using SVS software version 7 (Golden Helix Inc., Bozeman, MT, USA). Samples or SNPs with call rates < 0.9 were removed from the dataset. The minor allele frequencies (MAF) observed and expected heterozygosity, consistency of the genotype distributions with the Hardy-Weinberg equilibrium, definition of the linkage disequilibrium (LD), and LD blocks were estimated by the Haploview software package [31].

2.3. Milk Analysis

Analysis of the milk samples was performed within 24 h after collection, as reported by Dettori et al. [32]. As regards the fat, protein, casein, lactose content, and pH, these were achieved by using the standard of the International Organization for Standardization and International Dairy Federation standard [33] and the MilkoScan FT6000 instrument (Foss Electric A/S, Hillerød, Denmark); the total bacterial count (TBC) was measured according to the ISO-IDF standard [34] by using a BactoScan FC150 instrument (Foss Electric); the somatic cell count (SCC) was measured according to the ISO-IDF method [35] by using a Fossomatic 5000 instrument (Foss Electric). The TBC and SCC were transformed into respective logarithmic scores to normalize the distribution: LBC (log-bacterial count = log10 (total bacterial count/1000), according to the ISO-IDF standard [36])) and SCS (somatic cell score = log2 (SCC × 10−5) + 3, according to Shook [37]).
As regards the coagulation traits, both the traditional milk coagulation properties (MCP) and curd firmness over time traits (CFt) were measured. The following MCPs were measured according to the method by McMahon and Brown [4] and the revisions for sheep species by Pazzola et al. [5] by using a Formagraph instrument (Foss Italia, Padova, Italy): RCT (rennet coagulation time in min); k20 (curd firming time in min); and a30, a45, and a60 (curd firmness 30, 45, and 60 min after rennet addition in mm). The following CFt traits were measured according to a method used for ovine species by Vacca et al. [6] by using the curd-firming individual point observations recorded by a Formagraph instrument: CFP (the maximum potential curd firmness at an infinite time in mm), kCF (curd-firming rate constant in % × min1), kSR (syneresis rate constant in % × min−1), CFmax (maximum curd firmness in mm), and tmax (time to attain CFmax in min).

2.4. Statistical Analysis

As regards the editing of the MCP and CFt traits from the 1112 milk samples, the values of RCT higher than 60 min, also called noncoagulating samples, were labeled as missing (n = 26); k20 was labeled as missing when higher than 5 min and for noncoagulating samples (n = 36); curd firmness was labeled as missing at 0 mm (a30: n = 29, a45: n = 24, a60: n = 26, CFP: n = 19, and CFmax: n = 19).
The association analysis between the genotypes of the polymorphic SNPs with the milk composition, MCP, and CFt traits was achieved using the MIXED procedure of SAS (version 9.4, SAS Inst. Inc., Cary, NC, USA), which was derived from our previous study by Dettori et al. [32]. The following model (1) was used to investigate one milk trait for each SNP at a time only for SNPs with a MAF higher than 0.05 and genotypes with a frequency higher than 0.01:
Yijklm = µ + Gi + Fj + Pk + DIMl + SIRE(G)m + eijklm
where Yijklm is the observed traits of the milk composition, MCP, and CFt; µ is the general mean; Gi is the fixed effect of the ith SNP genotype (i = 2 to 3 levels: the two homozygotes and the heterozygote); Fj is the fixed effect of the jth farm, which also includes animal management and feeding (j = 1–23 levels, i.e., the different farms where animals were reared); Pk is the fixed effect of the kth parity of the ewes (k = 4 levels: first, 455 sheep; second, 215 sheep; third, 212 sheep; and fourth or more parities, 230 sheep); DIMl is the fixed effect of the lth days in milking (l = 4 levels: level 1, ≤100 days, 210 sheep; level 2, 101–140 days, 381 sheep; level 3, 141–160 days, 133 sheep; and level 4, ≥161 days, 388 sheep); SIRE(G)m is the random effect of the mth sire (m = 124 different sires) nested within the genotype; and eijklm is the error random residual effect.
The association analysis between the milk traits and the LD block was performed by a further model (2) derived from model (1), with LDi (i = 2 to 3 levels) instead of Gi and SIRE(LD)m instead of SIRE(G)m, analyzing one milk trait for each LD block at a time.
The effects of the models were declared significant at p < 0.05, and multiple testing, for both models (1) and (2) and for the factors with more than two levels was performed using the Bonferroni method at α = 0.05.

3. Results

3.1. Allele Frequencies and Linkage Disequilibrium

The results of 45 SNPs’ genotyping are reported in Table 1. Thirty (33.3%) SNPs were polymorphic and, among these, 12 showed a MAF higher than 0.05 and were submitted to the association analysis. A linkage disequilibrium (LD) analysis was performed for all the SNPs mapping to the same chromosome. LD was highlighted in chromosome OAR:16 (Figure 1) within the GHR gene. The SNPs analyzed ranged from the 5′UTR region to exon 10, the last exon of the GHR gene. One LD block was located at exon 10 and included two htSNPs: rs408890407, synonymous mutation, and rs55631463, missense mutation. Three haplotypes were recorded with CT showing the highest frequency, 0.493, while CC and TT had 0.277 and 0.230, respectively.

3.2. Association Analysis

The descriptive statistics of the milk composition and coagulation traits achieved from the 1112 sampled ewes is reported in Table 2. The fixed effects included in models (1) and (2), i.e., farm, parity, and stage of lactation, showed high levels of significance, at p < 0.001 and p < 0.01, for almost all the milk traits (data not shown in the tables). The results about these effects have been widely discussed in previous papers by the same authors [5,13] using datasets connected to the present one. The proportion of variance explained by the random effect of the sire computed in both the models was always lower than 10% for almost all the SNPs and the LD block. The lowest values were recorded for the CFp, kCF, kSR, Tmax, and LBC, at about 0%, and the highest for the fat, protein, casein, and lactose, spanning from 5% to 9% (data not shown in the tables).
Supplementary Table S2 displays the results of the statistical analysis computed by model (1) regarding the genotype effects of the 12 polymorphic SNPs, with the MAF > 0.05, on the milk traits. The least square means of the milk yield and composition, MCP, and CFt, according to the genotypes with a frequency > 0.01, are reported in Table 3. SNP rs119102735 at PRKAA2 was the only investigated polymorphic locus affecting the lactose content, and the milk samples from ewes carrying genotype TT were characterized by a lower concentration than CC and CT. The SNP with the most conspicuous effects on the milk traits was rs159573167 at PRKAG3, with AA homozygous ewes producing the highest daily milk yield, lower concentrations of protein and casein, the lowest value of SCS and LBC, and the shortest rennet coagulation time. SNP rs160202315 at the TLR4 (Toll-Like Receptor 4) gene was associated with the milk composition and cell traits, milk samples from heterozygous TC ewes characterized by the best contents, i.e., higher, of protein and casein, and the highest scores of SCS and LBC. The effects of rs409504706 at GHRHR and rs55631463 at GHR were directed on the kSR trait, with TT and CC ewes, respectively, showing the fastest syneresis process. SNP rs428862267 at GHR affected the curd firmness, and AA ewes produced milk with the highest a30. SNP rs160087383 at CD14 had a highly significant effect on the milk yield (p < 0.001) and coagulation, and CC homozygous ewes showed both higher daily yields and better coagulation traits than TT, i.e., faster rennet coagulation times and larger curd firmness (p < 0.05): dMY 2012 vs. 1654 g/day, RCT 7.60 vs. 9.28 min, and a30 52.52 vs. 48.03 mm, respectively.
The results of the statistical analysis computed by model (2) regarding the genotype effect of the haplotype blocks at chromosome 16 on the milk traits are reported in Supplementary Table S3. No significant effect on the investigated milk traits was found for the LD block at GHR (SNPs rs408890407 and rs55631463).

4. Discussion

4.1. Genes Related to Sheep Metabolism

Important metabolic variations occur during pregnancy and lactation, which involve nutrient partitioning in order to guarantee productive functions while maintaining the animal’s well-being [38]. In the present study, we evaluated the variability of genes responsible for metabolic activities and their association with milk traits in a population of Sarda sheep.
AMPK, or 5′AMP-activated protein kinase, is an enzyme (EC 2.7.11.31) formed by three protein subunits (α, β, and γ) conserved from yeast to humans. The AMPK enzyme helps to maintain the energy homeostasis of the cell, activating glucose and fatty acid uptake and their oxidation when the AMP:ATP ratio varies [39]. The alpha subunit is encoded by the PRKAA2 gene and the gamma subunit by the PRKAG3 gene. In the present study, we highlighted that SNP rs119102735 of the PRKAA2 gene was associated with the milk lactose concentration. Based on the available literature, this is the first study to report a significant association between the polymorphisms at this gene and a trait regarding milk composition in dairy sheep. In accordance with our study, a recent finding by Huang et al. [40] elucidated that, in bovine mammary epithelial cells, the activation of AMP-activated protein kinase systems is linked with the reduction of lactose synthesis. The AMPK protein, composed of three subunits in humans, is expressed in different tissues, mainly in the skeletal muscle, heart, liver, and hypothalamus [41]. Hao et al. [42] did not reveal it by RNA-seq among the mRNAs expressed in the lactating mammary gland of sheep. Bionaz-Loor et al. [43] revealed, by using Real-Time PCR, that PRKAA2 is expressed in the mammary glands of lactating cattle, and the same authors did not reveal variations in the expression of this gene during the analyzed period. They observed that, although high-protein synthesis persisted during lactation, the level of mRNA synthesis for the AMPK subunits, or other genes corresponding to related biosynthetic pathways, did not vary. In sheep, Crisà et al. [44] analyzed a polymorphism in exon 4 of PRKAA2, which was not associated with changes in the synthesis of fatty acids in milk. AMPK is a molecule that functions as a metabolic switch, which, in the absence of ATP, blocks consumption and activates the catabolic pathways to produce ATP, and turns off the anabolic pathways, which consume ATP [45]. Since AMPK is involved in linking energy sensing to the metabolic pathways, it plays an important role at the whole-body level by translating endocrine communications into adapted metabolic responses [46]. For this reason, it is possible that the observed association between PRKAA2 and lactose might not to be directly linked to a variation of the AMPK expression or activity in the secretory cell of the mammary gland but, rather, to the regulation of the animal’s metabolism in the physiological phase of lactation.
We also genotyped three SNPs of the PRKAG3 gene and evidenced that SNP rs159573167 AA was associated with an increase of dMY in the sampled population (+17% in comparison with genotype GG; data not shown in the tables) and a decrease of the protein and casein contents. The variability at the PRKAG3 gene was investigated in New Zealand Suffolk sheep [47], and the variability and expression profile in skeletal muscle were investigated in five local Spanish breeds, where PRKAG3 was found to be highly expressed [48]. Here, for the first time, we highlighted an association between PRKGA3 and dMY and the protein and casein contents. Milk from PRKAG3 rs159573167 AA ewes also showed the lowest SCS and LBC scores, which was not directly explained by its function, and the rennet coagulation time (RCT). In the present study, the occurrence of lower concentrations of protein, somatic, and bacterial cells and higher daily yields found in PRKAG3 rs159573167AA ewes might be attributable, in addition to the genetic effect, to the so-called dilution effect, as the higher the milk yield, the lower the level of total solids and cells and vice versa [1].
We also genotyped twelve SNPs from genes IGF1, GHRHR, and GHR, which are part of the growth hormone pathway and play an important role in lipid, protein, and glucose metabolism [27]. We highlighted an association with curd firmness a30 (GHR rs428862267) and the syneresis instant rate constant, kSR (GHR rs55631463 and GHRHR rs409504706). Previous studies, carried out in a subset of the resource sheep population, have revealed an association of several SNPs at the GHR gene with the daily fat and protein yields [49] and with the kSR [50]. The kSR is a derived parameter depending on the speed of syneresis, which, in turn, is related to the lipid content of the milk [7].

4.2. Genes Related to Innate Immunity

The role of specific genes in the innate immune response to some pathogens in ruminants has been investigated [51], due to its potential positive impact in sustainable agriculture systems. Here, we evaluated the association of some genes involved in innate immunity with milk traits in dairy Sarda sheep. Genotype CC at SNP rs160087383 of the CD14 gene was associated with an increase of dMY. Li et al. [52] investigated the role of the CD14 gene in relation to the Gram-negative bacteria-induced innate immune response and revealed that the variability of this gene is associated with the morbidity of Chinese Holsteins. Thus, the CD14 gene, which is expressed in macrophages and plays an important role in a number of pathways related to the immune response, in our set of animals was found to be associated with the daily milk yield. The CD14 rs160087383 CC ewes were characterized by the highest daily milk yield and the concurrence of favorable coagulation traits, i.e., a faster RCT and a larger curd firmness a30. As observed also for the PRKAG3 gene, the milk samples from high-yielding ewes were characterized by the shortest, and consequently better, values of the rennet coagulation time. This positive feature is in accordance with a study using a dataset interconnected with the present one [5] that has evidenced this peculiar result regarding sheep species, because cow milk yields and RCT commonly show a positive (i.e., unfavorable) correlation [53].
In our investigation, SNP rs160202315 at the TLR4 gene was associated with both the protein and cellular traits. The significant effect of the TLR4 gene on the protein and casein contents is analogous to that recorded for cattle species [54], while effects on milk coagulation have been found for TLR2, encoding for a protein of the same family, which influences the curd-firming time in dairy cows [55]. Milk samples from TLR4 rs160202315 CC homozygous ewes were characterized by the lowest scores of SCS and LBC; this association can be explained in relation to the gene function, pathogen recognition, and innate immunity, as described above.
We have to keep in mind that the observed associations between the genotyped SNPs and dairy traits may be due to linkages with unknown causal variations.

5. Conclusions

The designed open array revealed some new genetic patterns and the related influences on the milk yield, composition, and coagulation traits. The results obtained in the present study from the association analysis evidenced several significant effects of the investigated SNPs on the milk traits, and these may represent an advancement for the dairy sheep sector. In particular, for the Sarda sheep breed, this information might be useful in appropriate breeding schemes, but these should be confirmed by further studies, possibly performed on independent populations.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11082207/s1: Table S1: Gene names, SNP ID, and context sequences of the 45 SNPs genotyped in the sampled population of Sarda sheep (n = 1112). Table S2: p-value and significance for the milk yield and composition, milk coagulation properties (MCP), and curd firmness over time traits (CFt), according to the effects of 12 polymorphic SNPs out of 45 investigated in Sarda sheep (n = 1112). Table S3: p-value and significance for the milk yield and composition, milk coagulation properties (MCP), and curd firmness over time traits (CFt), according to the effects of the block evidenced among the 45 SNPs investigated in Sarda sheep (n = 1112).

Author Contributions

Conceptualization, M.P., G.M.V., G.B. and M.L.D.; methodology G.M.V.; software, M.P., P.P. and M.L.D.; formal analysis and resources, M.P., P.P. and M.L.D.; writing—original draft preparation, M.P. and M.L.D.; supervision, M.P. and G.M.V.; and funding acquisition, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by fund “Fondazione di Sardegna 2016-L.R. 7/2007, year 2016” CUP J86C18000830005, University of Sassari, Sassari, Italy.

Institutional Review Board Statement

Ethical review and approval were waived for this study. No specific authorization from an animal ethics committee was required. DNA was isolated from blood samples that were collected by official veterinarians of the local health authorities (ASL) during scheduled control or eradication sanitary programs not linked with the present study. The milk samples were collected concurrently with official and performance controls of the flock book not directly linked with the present study. The sheep belonged to commercial private flocks and were included in this trial with the agreement of the farmers on a voluntary basis.

Data Availability Statement

Not applicable.

Acknowledgments

The authors sincerely thank the sheep breeders for giving free access to their farms, the Regional Farmers Associations of Sardinia (Associazione Regionale Allevatori, A.R.A., Cagliari, Italy) for the technical support in milk sampling and analysis, the veterinarians for sample collection, and Cinzia Ribeca for collaboration with the laboratory activity.

Conflicts of Interest

The authors declare no conflict of interest. The funder had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Ramos, M.; Juarez, M. Sheep Milk. In Encyclopedia of Dairy Sciences, 2nd ed.; Fuquay, J.W., Fox, P.F., McSweeney, P.L.H., Eds.; Academic Press: San Diego, CA, USA, 2011; pp. 494–502. [Google Scholar]
  2. Statistics Division of the Food and Agriculture Organization of the United FAOSTAT. 2021. Available online: http://www.fao.org/faostat/en/#data/QC (accessed on 6 June 2018).
  3. Boyazoglu, J.; Morand-Fehr, P. Mediterranean dairy sheep and goat products and their quality: A critical review. Small Rumin. Res. 2001, 40, 1–11. [Google Scholar] [CrossRef]
  4. McMahon, D.J.; Brown, R.J. Evaluation of Formagraph for Comparing Rennet Solutions. J. Dairy Sci. 1982, 65, 1639–1642. [Google Scholar] [CrossRef]
  5. Pazzola, M.; Dettori, M.L.; Cipolat-Gotet, C.; Cecchinato, A.; Bittante, G.; Vacca, G.M. Phenotypic factors affecting coagulation properties of milk from Sarda ewes. J. Dairy Sci. 2014, 97, 7247–7257. [Google Scholar] [CrossRef]
  6. Vacca, G.M.; Pazzola, M.; Dettori, M.L.; Pira, E.; Malchiodi, F.; Cipolat-Gotet, C.; Cecchinato, A.; Bittante, G. Modeling of coagulation, curd firming and syneresis of milk from Sarda. J. Dairy Sci. 2015, 98, 2245–2259. [Google Scholar] [CrossRef] [Green Version]
  7. Bittante, G.; Pellattiero, E.; Malchiodi, F.; Cipolat-Gotet, C.; Pazzola, M.; Vacca, G.M.; Schiavon, S.; Cecchinato, A. Quality traits and modelling of coagulation, curd firming and syneresis of sheep milk of Alpine breeds fed diets supplemented with rumen-protected conjugated fatty acid. J. Dairy Sci. 2014, 97, 4018–4028. [Google Scholar] [CrossRef]
  8. Davis, G.H. Major genes affecting ovulation rate in sheep. Genet. Sel. Evol. 2005, 37, S11–S23. [Google Scholar] [CrossRef] [Green Version]
  9. Clop, A.; Marcq, F.; Takeda, H.; Pirottin, D.; Tordoir, X.; Bibé, B.; Bouix, J.; Caiment, F.; Elsen, J.-M.; Eychenne, F.; et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat. Genet. 2006, 38, 813–818. [Google Scholar] [CrossRef]
  10. Kemper, K.E.; Daetwyler, H.D.; Visscher, P.; Goddard, M. Comparing linkage and association analyses in sheep points to a better way of doing GWAS. Genet. Res. 2012, 94, 191–203. [Google Scholar] [CrossRef] [Green Version]
  11. Barrett, J.; Hansoul, S.; Nicolae, D.; Cho, J.H.; Duerr, R.; Rioux, J.D.; Brant, S.R.; Silverberg, M.S.; Taylor, K.D.; Barmada, M.M.; et al. Genome-wide association defines more than 30 distinct susceptibility loci for Crohn’s disease. Nat. Genet. 2008, 40, 955–962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Gebreselassie, G.; Berihulay, H.; Jiang, L.; Ma, Y. Review on Genomic Regions and Candidate Genes Associated with Economically Important Production and Reproduction Traits in Sheep (Ovies aries). Animals 2019, 10, 33. [Google Scholar] [CrossRef] [Green Version]
  13. Noce, A.; Pazzola, M.L.M.; Dettori, M.; Amills, A.; Castelló, A.; Cecchinato, G.; Bittante, G.; Vacca, G.M. Variations at regu-latory regions of the milk protein genes are associated with milk traits and coagulation properties in the Sarda sheep. Anim. Genet. 2016, 47, 717–726. [Google Scholar] [CrossRef]
  14. Usai, M.G.; Casu, S.; Sechi, T.; Salaris, S.L.; Miari, S.; Sechi, S.; Carta, P.; Carta, A. Mapping genomic regions affecting milk traits in Sarda sheep by using the OvineSNP50 Beadchip and principal components to perform combined linkage and linkage dis-equilibrium analysis. Genet. Sel. Evol. 2019, 51, 65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Bittante, G.; Cipolat-Gotet, C.; Pazzola, M.; Dettori, M.L.; Vacca, G.M.; Cecchinato, A. Genetic analysis of coagulation proper-ties, curd firming modeling, milk yield, composition and acidity in Sarda dairy sheep. J. Dairy Sci. 2017, 100, 385–394. [Google Scholar] [CrossRef]
  16. The Human Gene Database Website. Genecards 2021. Available online: https://www.genecards.org (accessed on 18 March 2021).
  17. Khan, M.; Hosseini, A.; Burrell, S.; Rocco, S.; McNamara, J.; Loor, J. Change in subcutaneous adipose tissue metabolism and gene network expression during the transition period in dairy cows, including differences due to sire genetic merit. J. Dairy Sci. 2013, 96, 2171–2182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Xia, W.; Osorio, J.; Yang, Y.; Liu, D.; Jiang, M.F. Short communication: Characterization of gene expression profiles related to yak milk protein synthesis during the lactation cycle. J. Dairy Sci. 2018, 101, 11150–11158. [Google Scholar] [CrossRef] [PubMed]
  19. Winder, W.W.; Holmes, B.F.; Rubink, D.S.; Jensen, E.B.; Chen, M.; Holloszy, J.O. Activation of AMP-activated protein kinase increases mitochondrial enzymes in skeletal muscle. J. Appl. Physiol. 2000, 88, 2219–2226. [Google Scholar] [CrossRef] [Green Version]
  20. Yu, S.; Kim, J.; Chung, H.; Jung, K.; Lee, Y.; Yoon, D.; Lee, S.; Choi, I.; Bottema, C.; Sang, B.; et al. Molecular cloning and characterization of bovine PRKAG3 gene: Structure, expression and single nucleotide polymorphism detection. J. Anim. Breed. Genet. 2005, 122, 294–301. [Google Scholar] [CrossRef]
  21. Bongiorni, S.; Gruber, C.E.M.; Bueno, S.; Chillemi, G.; Ferre’, F.; Failla, S.; Moioli, B.; Valentini, A. Transcriptomic investigation of meat tenderness in two Italian cattle breeds. Anim. Genet. 2016, 47, 273–287. [Google Scholar] [CrossRef]
  22. Dall’Olio, S.; Scotti, E.; Costa, L.N.; Fontanesi, L. Effects of single nucleotide polymorphisms and haplotypes of the proteinkinase AMP-activated non-catalytic subunit gamma 3 (PRKAG3) gene onproduction, meat quality and carcass traits in Italian Large White pigs. Meat Sci. 2018, 136, 44–49. [Google Scholar] [CrossRef]
  23. Bergan-Roller, H.E.; Sheridan, M.A. The growth hormone signaling system: Insights into coordinating the anabolic and cata-bolic actions of growth hormone. Gen. Comp. Endocrinol. 2018, 258, 119–133. [Google Scholar] [CrossRef]
  24. Osorio, J.S.; Lohakare, J.; Bionaz, M. Biosynthesis of milk fat, protein, and lactose: Roles of transcriptional and posttranscrip-tional regulation. Physiol. Genom. 2016, 48, 231–256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Klein, I.; Sarkadi, B.; Váradi, A. An inventory of the human ABC proteins. Biochim. Biophys. Acta (BBA) Biomembr. 1999, 1461, 237–262. [Google Scholar] [CrossRef] [Green Version]
  26. Wassermann, L.; Halwachs, S.; Baumann, D.; Schaefer, I.; Seibel, P.; Honscha, W. Assessment of ABCG2-mediated transport of xenobiotics across the blood–milk barrier of dairy animals using a new MDCKII in vitro model. Arch. Toxicol. 2013, 87, 1671–1682. [Google Scholar] [CrossRef]
  27. Fleminger, G.; Heftsi, R.; Uzi, M.; Nissim, S.; Gabriel, L. Chemical and structural characterization of bacterially-derived casein peptides that impair milk clotting. Int. Dairy J. 2011, 21, 914–920. [Google Scholar] [CrossRef]
  28. Jhamandas, J.H.; Goncharuk, V. Role of neuropeptide FF in central cardiovascular and neuroendocrine regulation. Front. Endocrinol. 2013, 4, 8. [Google Scholar] [CrossRef] [Green Version]
  29. Sahana, G.; Guldbrandtsen, B.; Thomsen, B.; Holm, L.-E.; Panitz, F.; Brøndum, R.; Bendixen, C.; Lund, M. Genome-wide association study using high-density single nucleotide polymorphism arrays and whole-genome sequences for clinical mastitis traits in dairy cattle. J. Dairy Sci. 2014, 97, 7258–7275. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Paris, J.M.; Letko, A.A.; Häfliger, I.M.; Ammann, P.; Flury, C.; Drögemüller, C. Identification of two TYRP1 loss-of-function alleles in Valais Redsheep. Anim. Genet. 2019, 50, 778–782. [Google Scholar] [CrossRef]
  31. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [Green Version]
  32. Dettori, M.L.; Pazzola, M.; Petretto, E.; Vacca, G.M. Association Analysis between SPP1, POFUT1 and PRLR Gene Variation and Milk Yield, Composition and Coagulation Traits in Sarda Sheep. Animals 2020, 10, 1216. [Google Scholar] [CrossRef] [PubMed]
  33. ISO (International Organization for Standardization); IDF (International Dairy Federation). Milk and Liquid Milk Products: Determination of Fat, Protein, Casein, Lactose and pH Content; International Standard ISO 9622; IDF 141:2013; ISO: Geneva, Switzerland; IDF: Brussels, Belgium, 2013; Available online: https://www.iso.org (accessed on 5 April 2020).
  34. ISO (International Organization for Standardization); IDF (International Dairy Federation). Milk: Quantitative Determination of Bacteriological Quality—Guidance for Establishing and Verifying a Conversion Relationship between Routine Method Results and Anchor Method Results; International Standard ISO 21187; IDF 196:2004; ISO: Geneva, Switzerland; IDF: Brussels, Belgium, 2004; Available online: https://www.iso.org (accessed on 5 April 2020).
  35. ISO (International Organization for Standardization); IDF (International Dairy Federation). Milk: Enumeration of Somatic Cells—Part 2: Guidance on the Operation of Fluoro-Opto-Electronic Counters; International Standard ISO 13366-2; IDF IDF 148-2:2006; ISO: Geneva, Switzerland; IDF: Brussels, Belgium, 2006; Available online: https://www.iso.org (accessed on 5 April 2020).
  36. ISO (International Organization for Standardization); IDF (International Dairy Federation). Milk: Bacterial Count—Protocol for the Evaluation of Alternative; International Standard ISO 16297:2013; IDF 161:2013; ISO: Geneva, Switzerland; IDF: Brussels, Belgium, 2013; Available online: https://www.iso.org (accessed on 5 April 2020).
  37. Shook, G.E. Genetic Improvement of Mastitis Through Selection on Somatic Cell Count. Vet. Clin. N. Am. Food Anim. Pract. 1993, 9, 563–577. [Google Scholar] [CrossRef]
  38. Baumgard, L.; Collier, R.; Bauman, D. A 100-Year Review: Regulation of nutrient partitioning to support lactation. J. Dairy Sci. 2017, 100, 10353–10366. [Google Scholar] [CrossRef]
  39. Suter, M.; Riek, U.; Tuerk, R.; Schlattner, U.; Wallimann, T.; Neumann, D. Dissecting the role of 5’AMP for allosteric stimu-lation, activation, and deactivation of AMP-activated protein kinase. J. Biol. Chem. 2006, 281, 32207–32216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Huang, J.; Guesthier, M.-A.; Burgos, S. AMP-activated protein kinase controls lipid and lactose synthesis in bovine mammary epithelial cells. J. Dairy Sci. 2020, 103, 340–351. [Google Scholar] [CrossRef] [PubMed]
  41. Hardie, D.G.; Lin, S.-C. AMP-activated protein kinase—Not just an energy sensor. F1000Research 2017, 6, 1724. [Google Scholar] [CrossRef]
  42. Zhao, M.; Zhou, H.; Hickford, J.G.; Gong, H.; Wang, J.; Hu, J.; Liu, X.; Li, S.; Hao, Z.; Luo, Y. Transcriptome Profile Analysis of Mammary Gland Tissue from Two Breeds of Lactating Sheep. Genes 2019, 10, 781. [Google Scholar] [CrossRef] [Green Version]
  43. Bionaz, M.; Loor, J.J. Gene networks driving bovine mammary protein synthesis during the lactation cycle. Bioinform. Biol. Insights 2011, 5, 83–98. [Google Scholar] [CrossRef]
  44. Crisà, A.; Marchitelli, C.; Pariset, L.; Contarini, G.; Signorelli, F.; Napolitano, F.; Catillo, G.; Valentini, A.; Moioli, B. Exploring polymorphisms and effects of candidate genes on milk fat quality in dairy sheep. J. Dairy Sci. 2010, 93, 3834–3845. [Google Scholar] [CrossRef] [Green Version]
  45. Hardie, D.G. AMPK: A key regulator of energy balance in the single cell and the whole organism. Int. J. Obes. 2008, 32, S7–S12. [Google Scholar] [CrossRef] [Green Version]
  46. Nguyen, T.M.D. Role of AMPK in mammals reproduction: Specific controls and whole-body energy sensing. Comptes Rendus Biol. 2019, 342, 1–6. [Google Scholar] [CrossRef] [PubMed]
  47. Yang, G.; Zhou, H.; Wang, R.; Hickford, J. Variation in the ovine PRKAG3 gene. Gene 2015, 567, 251–254. [Google Scholar] [CrossRef]
  48. Noce, A.; Cardoso, T.F.; Manunza, A.; Martinez, A.; Cánovas, A.; Pons, A.; Bermejo, L.A.; Landi, V.; Sànchez, A.; Jordana, J.; et al. Expression patterns and genetic variation of the ovine skeletal muscle transcriptome of sheep from five Spanish meat breeds. Sci. Rep. 2018, 8, 10486. [Google Scholar] [CrossRef]
  49. Dettori, M.L.; Pazzola, M.; Paschino, P.; Amills, M.; Vacca, G.M. Association between the GHR, GHRHR and IGF1 gene pol-ymorphisms and milk yield and quality traits in Sarda sheep. J. Dairy Sci. 2018, 101, 9978–9986. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Dettori, M.L.; Pazzola, M.; Pira, E.; Stocco, G.; Vacca, G.M. Association between the GHR, GHRHR and IGF1 gene polymor-phisms and milk coagulation properties in Sarda sheep. J. Dairy Res. 2019, 86, 331–336. [Google Scholar] [CrossRef]
  51. Dhanasekaran, S.; Vignesh, A.R.; Raj, G.D.; Reddy, Y.K.M.; Raja, A.; Tirumurugaan, K.G. Comparative analysis of innate immune response following in vitro stimulation of sheep and goat peripheral blood mononuclear cells with bluetongue virus—serotype 23. Vet. Res. Commun. 2013, 37, 319–327. [Google Scholar] [CrossRef] [PubMed]
  52. Li, J.; Wang, Q.; Chen, F.; Wang, H.; Chen, J.; Wang, Z.; Huo, J.; Cai, Y. SNPs of CD14 change the mastitis morbidity of Chinese Holstein. Mol. Med. Rep. 2017, 16, 9102–9110. [Google Scholar] [CrossRef]
  53. Bittante, G.; Penasa, M.; Cecchinato, A. Invited review: Genetics and modeling of milk coagulation properties. J. Dairy Sci. 2012, 95, 6843–6870. [Google Scholar] [CrossRef]
  54. Viale, E.; Tiezzi, F.; Maretto, F.; De Marchi, M.; Penasa, M.; Cassandro, M. Association of candidate gene polymorphisms with milk technological traits, yield, composition, and somatic cell score in Italian Holstein-Friesian sires. J. Dairy Sci. 2017, 100, 7271–7281. [Google Scholar] [CrossRef]
  55. Cecchinato, A.; Chessa, S.; Ribeca, C.; Cipolat-Gotet, C.; Bobbo, T.; Casellas, J.; Bittante, G. Genetic variation and effects of candidate-gene polymorphisms on coagulation properties, curd firmness modeling and acidity in milk from Brown Swiss cows. Animal 2015, 9, 1104–1112. [Google Scholar] [CrossRef] [Green Version]
Figure 1. LD block defined by the SNPs genotyped in Sarda sheep (n = 1112), shown with haplotype frequencies on the right. LD plot of the SNPs on chromosome 16, including the GHR gene. Haploview plot of pairwise D′: red, D′ = 1.0 and logarithm of the odds (LOD) ≥ 2.0; white, D′ < 1.0 and LOD < 2.0. LD blocks are delimited by a black line. Haplotype-tagging SNPs are within black boxes. SNPs in bold are haplotype-tagging SNPs (htSNP).
Figure 1. LD block defined by the SNPs genotyped in Sarda sheep (n = 1112), shown with haplotype frequencies on the right. LD plot of the SNPs on chromosome 16, including the GHR gene. Haploview plot of pairwise D′: red, D′ = 1.0 and logarithm of the odds (LOD) ≥ 2.0; white, D′ < 1.0 and LOD < 2.0. LD blocks are delimited by a black line. Haplotype-tagging SNPs are within black boxes. SNPs in bold are haplotype-tagging SNPs (htSNP).
Animals 11 02207 g001
Table 1. The variability of 45 SNPs genotyped in Sarda sheep (n = 1112).
Table 1. The variability of 45 SNPs genotyped in Sarda sheep (n = 1112).
Gene SNP IDObsHPredHHWpv%GenMAFAlleles
PRKAA2: Protein Kinase AMP-Activated Catalytic Subunit Alpha2
rs1191027350.3190.3160.9196.60.197C:(T)
rs1597014430.0040.0041.0097.80.002G:(A)
rs1597014410.1610.1570.5786.20.086T:(C)
rs1597014400.0110.0190.0191.80.009G:(A)
PRKAG3: Protein Kinase AMP-Activated Non-Catalytic Subunit Gamma3
rs1595731400.0000.0001.001000.000G:G
rs1595731090.0010.0011.0099.00.000T:T
rs1595731670.0080.4350.0196.80.320G:(A)
B4GALT1: Beta-1,4-Galactosyltransferase 1
rs1601760290.0130.0180.9698.60.009G:(A)
rs1601760200.0000.0001.0099.90.000A:A
rs1601758210.0000.0001.001000.000A:A
rs1601758090.4190.4150.8579.50.294C:(T)
TLR4: Toll Like Receptor 4
rs1602023300.0000.0001.0099.30.000A:A
rs1602023210.0050.0051.0099.50.003A:(C)
rs1602023150.0530.4070.0197.50.284T:(C)
TYRP1: Tyrosinase Related Protein 1
rs4164172090.3230.3350.3196.60.212T:(C)
IGF1: Insulin Like Growth Factor 1
rs1598763900.3390.4590.0186.70.356G:(A)
rs1598763940.0020.0021.0099.50.001C:(G)
rs4010287810.0000.0001.0098.80.000G:G
GHRHR: Growth Hormone Releasing Hormone Receptor
rs4095047060.3140.3430.0195.50.220T:(G)
rs1617972460.0040.0041.0094.50.002C:(T)
CD14: Cluster of Differentiation 14 Molecule
rs1600873650.0010.0011.0098.30.000C:C
rs1600873710.0010.0080.0199.20.004G:(C)
rs1600873780.0000.0001.0099.00.000C:C
rs1600873830.0030.1210.0199.40.064T:(C)
ABCG2: ATP Binding Cassette Subfamily G Member 2
rs1599568450.0000.01.0099.30.000T:T
rs1599568850.0040.0050.0198.90.003A:(G)
rs1599569740.0020.0230.0187.40.011T:(G)
NPFFR2: Neuropeptide FF Receptor 2
rs1599805900.0340.0331.0091.50.017A:(C)
rs1599805930.0130.0150.1196.90.007T:(A)
CSNK1G1: Casein Kinase 1 Gamma 1
rs1603223860.0000.0001.00 1000.000G:G
GHR: Growth Hormone Receptor
rs4088904070.3260.3540.0196.30.230C:(T)
rs1611461640.0120.0121.0098.70.006T:(G)
rs556314630.3950.3990.7996.70.276T:(C)
rs4137760540.0100.0101.0099.50.005G:(A)
rs4042373210.0120.0121.0099.50.006C:(T)
rs4288622670.4480.4900.0194.40.430G:(A)
rs4258345830.0000.0001.001000.000A:A
TRL2: Toll Like Receptor 2
rs1620733180.0010.0011.0099.50.000C:C
PI: Serpin Peptidase Inhibitor
rs3975141700.0000.0001.001000.000A:A
rs3975141690.0010.0011.0099.50.000G:G
rs3975141680.0000.0001.0099.70.000T:T
rs3975141500.0020.0021.0098.90.001T:(C)
PLCE1: Phospholipase C Epsilon
rs1614731240.0610.0590.7399.00.030C:(T)
rs1614731260.0010.0030.0199.00.001T:(C)
rs1614731400.0060.0061.0090.90.003G:(A)
ObsH: observed heterozygosity, PredH: predicted heterozygosity, HWpv: Hardy-Weinberg test p-value, %Gen: percentage of genotyped samples, MAF: minor allele frequency, and minor alleles in brackets.
Table 2. Descriptive statistics of the milk yield and composition, milk coagulation properties (MCP), and curd firmness over time traits (CFt) from the sampled population of Sarda sheep (n = 1112).
Table 2. Descriptive statistics of the milk yield and composition, milk coagulation properties (MCP), and curd firmness over time traits (CFt) from the sampled population of Sarda sheep (n = 1112).
Milk TraitsMeanSDMinMax
Milk yield and composition
dMY (g/day)15978421815760
Fat (g/100 mL)6.441.191.9312.52
Protein (g/100 mL)5.410.653.528.24
Casein (g/100 mL)4.210.542.456.67
Lactose (g/100 mL)4.810.312.405.59
pH6.670.096.457.14
SCC (n × 1000/mL)12532993520,612
TBC (n × 1000/mL)25224872117,000
MCP
RCT (min)8.944.921.3052.30
k20 (min)1.980.791.3011.30
a30 (mm)49.4812.99070.00
a45 (mm)45.4615.17073.88
a60 (mm)41.8416.56077.98
CFt
CFP (mm)62.1316.20360.8
kCF (% × min−1)27.013.01.090.0
kSR (% × min−1)1.51.90.116.3
CFmax (mm)54.179.96075.96
tmax (min)30.5713.418.7560
dMY: daily milk yield; SCC: somatic cell count; TBC: total bacterial count; RCT: rennet coagulation time; k20: curd-firming time; a30, a45, and a60: curd firmness 30, 45, and 60 min after rennet addition; CFP: the maximum potential curd firmness after an infinite time; kCF: curd-firming rate constant; kSR: syneresis rate constant; CFmax: maximum curd firmness; and tmax: time to attain CFmax.
Table 3. Least square means of the milk yield and composition, the milk coagulation properties (MCP), and curd firmness over time traits (CFt), according to the genotypes of the SNPs investigated in the Sarda sheep (n = 1112).
Table 3. Least square means of the milk yield and composition, the milk coagulation properties (MCP), and curd firmness over time traits (CFt), according to the genotypes of the SNPs investigated in the Sarda sheep (n = 1112).
Milk Yield and CompositionMCPCFt
Gene SNP ID #Genotype and (n)dMYProteinCaseinLactoseSCSLBCRCTa30kSR
PRKAA2
rs119102735
CC (695)16815.394.204.79 b4.842.429.1248.471.46
CT (340)16575.364.174.79 b4.932.389.0847.941.65
TT (41)16795.424.234.68 a5.032.289.2645.471.46
PRKAG3
rs159573167
AA (346)1859 B5.24 a4.09 a4.844.32 A1.77 A8.03 a49.201.18
AG (11)1851 AB5.49 b4.29 b4.844.85 B2.22 AB10.79 b47.731.82
GG (724)1585 A5.45 b4.25 b4.765.17 B2.62 B9.77 b47.711.64
TLR4
rs160202315
CC (286)17475.21 a4.06 a4.854.26 a1.63 A8.3547.871.02
TC (58)17425.53 b4.32 b4.764.93 b2.50 B9.7748.381.95
TT (746)16375.44 b4.24 ab4.765.15 b2.61 B9.5048.521.71
GHRHR
rs409504706
GG (67)16455.334.154.765.322.319.7347.622.14 b
GT (334)16845.364.184.794.792.478.9247.861.58 ab
TT (662)16805.384.194.794.852.369.1448.431.42 a
CD14
rs160087383 §
CC (70)2012 B5.324.154.814.402.037.60 a52.52 b1.45
TT (1037)1654 A5.384.204.784.902.399.28 b48.03 a1.50
GHR
rs55631463
CC (85)17435.364.184.785.002.339.3648.611.11 A
CT (426)16455.374.194.774.822.419.3847.381.79 B
TT (567)16805.384.194.794.942.398.9848.701.39 A
GHR
rs428862267
AA (214)17375.394.194.774.932.278.8850.21 b1.28
GA (470)16485.374.194.804.802.328.8847.95 a1.61
GG (366)16755.364.174.794.902.479.5048.05 ab1.51
# SNPs out of the 45 investigated with MAP > 0.05, genotypes with frequencies > 0.01, and at least one significant p-value in one of the 18 investigated milk traits. § Genotype CT not computed because of the frequency < 0.01. dMY: daily milk yield in g/day. Protein, casein, and lactose in g/100 mL. SCS: somatic cell score = log2 (SCC × 10−5) + 3. LBC: logarithmic bacterial count = log10 total bacterial count (total bacterial count/1000). RCT: rennet coagulation time in min. a30: curd firmness 30 min after rennet addition in mm. kSR: syneresis rate constant in % × min−1. a,b Different lower-case letters for the same SNP differ significantly at p < 0.05. A,B Different upper-case letters for the same SNP differ significantly at p < 0.01. Least square means for the milk traits with statistical significance are in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pazzola, M.; Vacca, G.M.; Paschino, P.; Bittante, G.; Dettori, M.L. Novel Genes Associated with Dairy Traits in Sarda Sheep. Animals 2021, 11, 2207. https://doi.org/10.3390/ani11082207

AMA Style

Pazzola M, Vacca GM, Paschino P, Bittante G, Dettori ML. Novel Genes Associated with Dairy Traits in Sarda Sheep. Animals. 2021; 11(8):2207. https://doi.org/10.3390/ani11082207

Chicago/Turabian Style

Pazzola, Michele, Giuseppe Massimo Vacca, Pietro Paschino, Giovanni Bittante, and Maria Luisa Dettori. 2021. "Novel Genes Associated with Dairy Traits in Sarda Sheep" Animals 11, no. 8: 2207. https://doi.org/10.3390/ani11082207

APA Style

Pazzola, M., Vacca, G. M., Paschino, P., Bittante, G., & Dettori, M. L. (2021). Novel Genes Associated with Dairy Traits in Sarda Sheep. Animals, 11(8), 2207. https://doi.org/10.3390/ani11082207

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop