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

Milk-y Way: the impact of single-nucleotide polymorphisms on milk production traits in Kilis dairy goats

This study aimed to investigate the impact of single-nucleotide polymorphisms (SNPs) on milk production traits in Kilis dairy goats by analyzing the genotypes of POU1F1, PRLR, β -Lg, GH1, and GH2 genes and their association with lactation milk yield (LMY), lactation length (LL) and average daily milk yield (ADMY). Blood samples were collected from 227 goats, and genotyping was performed using polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP). The results revealed that the frequencies of the genotypes varied among the genes. The polymorphisms were found to be significantly linked with milk production traits. These findings suggest that SNPs of POU1F1/AluI, PRLR/RsaI and β -Lg/SacII are significantly associated with LMY and that the POU1F1-TC genotype, PRLR-TT genotype and β -Lg-AB genotype are associated with higher LMY and ADMY. Additionally, the POU1F1-TC genotype was found to have a longer LL. However, no significant association was found between the GH1 and GH2 genotypes and LMY, LL and ADMY. Overall, this study provides valuable insights into the genetic factors influencing milk production traits in Kilis dairy goats, which can be utilized for the selection of high-yielding animals in breeding programs.

Genome-wide association studies of lactation yields of milk, fat, protein and somatic cell score in New Zealand dairy goats

Background

Identifying associations between genetic markers and traits of economic importance will provide practical benefits for the dairy goat industry, enabling genomic prediction of the breeding value of individuals, and facilitating discovery of the underlying genes and mutations. Genome-wide association studies were implemented to detect genetic regions that are significantly associated with effects on lactation yields of milk (MY), fat (FY), protein (PY) and somatic cell score (SCS) in New Zealand dairy goats.

Methods

A total of 4,840 goats were genotyped with the Caprine 50 K SNP chip (Illumina Inc., San Diego, CA). After quality filtering, 3,732 animals and 41,989 SNPs were analysed assuming an additive linear model. Four GWAS models were performed, a single-SNP additive linear model and three multi-SNP BayesC models. For the single-SNP GWAS, SNPs were fitted individually as fixed covariates, while the BayesC models fit all SNPs simultaneously as random effects. A cluster of significant SNPs were used to define a haplotype block whose alleles were fitted as covariates in a Bayesian model. The corresponding diplotypes of the haplotype block were then fit as class variables in another Bayesian model.

Results

Across all four traits, a total of 43 genome-wide significant SNPs were detected from the SNP GWAS. At a genome-wide significance level, the single-SNP analysis identified a cluster of variants on chromosome 19 associated with MY, FY, PY, and another cluster on chromosome 29 associated with SCS. Significant SNPs mapped in introns of candidate genes (45%), in intergenic regions (36%), were 0-5 kb upstream or downstream of the closest gene (14%) or were synonymous substitutions (5%). The most significant genomic window was located on chromosome 19 explaining up to 9.6% of the phenotypic variation for MY, 8.1% for FY, 9.1% for PY and 1% for SCS.

Conclusions

The quantitative trait loci for yield traits on chromosome 19 confirms reported findings in other dairy goat populations. There is benefit to be gained from using these results for genomic selection to improve milk production in New Zealand dairy goats.

Characterization of milk composition, coagulation properties, and cheese-making ability of goats reared in extensive farms

The aims of this study were to explore the variability of milk composition, coagulation properties, and cheese-making traits of the Sarda goat breed, and to investigate the effects of animal and farm factors, and the geographic area (Central-East vs. South-West) of an insular region of Italy, Sardinia. A total of 570 Sarda goats reared in 21 farms were milk-sampled during morning milking. Individual milk samples were analyzed for composition, traditional milk coagulation properties (MCP), modeled curd-firming over time parameters (CF_t), and cheese-making traits (cheese yield, %CY; recovery of nutrients, %REC; daily cheese yield, dCY). Farms were classified into 2 categories based on milk energy level (MEL; high or low), defined according to the average net energy of milk daily produced by the lactating goats. Milk yield and composition were analyzed using a mixed model including the fixed effects of MEL, geographic area, days in milk, and parity, and the random effect of farm within MEL and geographic area. Data about MCP, CF_t, and the cheese-making process were analyzed using the same model, with the inclusion of the effects of animal and pendulum of the lactodynamograph instrument, allowing the measure of repeatability of these traits. Results showed that animal had greater influence on coagulation and cheese-making traits compared with farm effect. Days in milk influenced milk composition, whose changes partly reflected the modifications of %CY traits. Moreover, large differences were observed between primiparous and multiparous goats: primiparous goats produced less milk of better quality (higher fat, lower somatic cell and bacterial counts) and less cheese, but with higher recovery of fat and protein in the curd, compared with multiparous goats. The repeatability was very high, for both coagulation (84.0 to 98.8%) and cheese-making traits (89.7 to 99.9%). The effect of MEL was significant for daily productions of milk and cheese, coagulation time, and recovery of protein in the curd, which were better in high-MEL farms. As regards geographic area, milk composition and percentage cheese yield were superior in the Central-East area, whereas daily milk and cheese production and MCP were better in the South-West. This result was explainable by the phenomenon of crossbreeding Sarda goats with Maltese bucks, which occurred with greater intensity in the South-West than in the Central-East area of the island. The results provided by this study could be of great interest for the goat dairy sector. Indeed, the methods described in the present study could be applicable for other farming methods, goat breeds, and geographic areas. The collection of a wide range of phenotypes at individual animal level is fundamental for the characterization of local populations and can be used to guarantee breed conservation and the persistence of traditional farming systems, and to increase farmers' profit.

Whey protein polymorphisms in Sudanese goat breeds

The aim of the present study was to assess genetic variation that is characteristic for Sudanese goat breeds in the milk whey protein genes (LALBA and BLG). Four Sudanese goat breeds were screened for variability in LALBA and BLG genes at the DNA level by comparative sequencing of five animals per breed. Sixteen SNPs were identified in LALBA: seven in the upstream region, six synonymous, and three in the 3´-UTR. Three novel synonymous SNPs in exon 2 (ss5197800003, ss5197800012, and ss5197800004) were found in Nubian, Desert, and Nilotic, but not in Taggar goats. One SNP in the promoter of LALBA (rs642745519) modifies a predicted transcription factor binding site for Tcfe2a. The SNPs in the 3'-UTR (rs657915405, rs641559728, and rs664225585) affect predicted miRNA target sites. With respect to haplotypes in the exonic region, haplotype LALBA-A is most frequent in Nubian, Desert, and Nilotic goats, while haplotype LALBA-D is prevalent in Taggar goats. In BLG, 30 SNPs were detected: eight in the upstream gene region, two synonymous, 17 intronic, and three in the 3'-UTR. Among the 30 identified SNPs, 15 were novel. Four of these novel SNPs were located in the upstream gene region, one was synonymous, and ten were intronic. The novel synonymous SNP (ss5197800017), located in exon 2, was only found in Nubian and Nilotic goats. The SNPs ss5197800010 and rs635615192 in the promoter are located in predicted binding sites of transcription factors (M6097, Elk3, Elf5, and GABPA). Among seven haplotypes detected in the coding region, haplotype BLG-A is most frequent in Nubian and Nilotic goats while haplotype BLG-B is most frequent in Desert and Taggar goats. The high variability in regulatory gene regions among Sudanese goats could potentially affect the quality and yield of whey proteins in goat milk and provide a wide resource for genetic improvement of milk production and milk technology characteristics.

Genome-wide association study of conformation and milk yield in mixed-breed dairy goats

Identification of genetic markers that affect economically important traits is of high value from a biological point of view, enabling the targeting of candidate genes and providing practical benefits for the industry such as wide-scale genomic selection. This study is one of the first to investigate the genetic background of economically important traits in dairy goats using the caprine 50K single nucleotide polymorphism (SNP) chip. The aim of the project was to perform a genome-wide association study for milk yield and conformation of udder, teat, and feet and legs. A total of 137,235 milk yield records on 4,563 goats each scored for 10 conformation traits were available. Out of these, 2,381 goats were genotyped with the Illumina Caprine 50K BeadChip (Illumina Inc., San Diego, CA). A range of pseudo-phenotypes were used including deregressed breeding values and pseudo-estimated breeding values. Genome-wide association studies were performed using the multi-locus mixed model (MLMM) algorithm implemented in SNP & Variation Suite v7.7.8 (Golden Helix Inc., Bozeman, MT). A genome-wise significant [-log₁₀(P-value) > 5.95] SNP for milk yield was identified on chromosome 19, with additional chromosome-wise significant (-log₁₀(P-value) > 4.46] SNP on chromosomes 4, 8, 14, and 29. Three genome-wise significant SNP for conformation of udder attachment, udder depth, and front legs were identified on chromosome 19, and chromosome-wise SNP were found on chromosomes 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 21, 23, and 27. The proportion of variance explained by the significant SNP was between 0.4 and 7.0% for milk yield and between 0.1 and 13.8% for conformation traits. This study is the first attempt to identify SNP associated with milk yield and conformation in dairy goats. Two genome-wise significant SNP for milk yield and 3 SNP for conformation of udder attachment, udder depth, and front legs were found. Our results suggest that conformation traits have a polygenic background because, for most of them, we did not identify any quantitative trait loci with major effect.