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

Megan Scholtens,Andrew Jiang,Ashley Smith,Mathew Littlejohn,Klaus Lehnert,Russell Snell,Nicolas Lopez-Villalobos,Dorian Garrick,Hugh Blair

DOI: https://doi.org/10.1186/s40104-020-00453-2

2020-05-25

Journal of Animal Science and Biotechnology

Abstract:Abstract 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.

agriculture, dairy & animal science

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to identify the genetic regions significantly associated with milk production traits in New Zealand dairy goats through genome - wide association studies (GWAS). Specifically, the authors hope to find gene loci significantly associated with the following four milk production traits by analyzing single - nucleotide polymorphism (SNP) markers:

1. **Milk Yield (MY)**: The total amount of milk.
2. **Fat Yield (FY)**: The yield of milk fat.
3. **Protein Yield (PY)**: The yield of milk protein.
4. **Somatic Cell Score (SCS)**: Used to assess the number of somatic cells in milk, indirectly reflecting udder health.

### Research Background

In the dairy goat industry, identifying genetic markers associated with economic traits can provide important tools for breeding, thereby improving the accuracy of individual breeding value prediction and helping to discover potential functional genes and mutations. This information can help accelerate the breeding process, shorten the generation interval, and increase the genetic gain rate.

### Research Methods

The researchers genotyped 4,840 dairy goats using the Illumina Caprine 50K SNP chip. After quality filtering, they finally analyzed 3,732 animals and 41,989 SNPs. To detect the genetic regions significantly associated with the above four milk production traits, the researchers adopted four different GWAS models:

- **Single - SNP additive linear model**: Each SNP was fitted separately as a fixed covariate.
- **Three multi - SNP BayesC models**: All SNPs were simultaneously fitted as random effects.

In addition, the researchers also defined haplotype blocks by clustering significant SNPs and incorporated the alleles of these haplotype blocks as covariates into the Bayesian model for further analysis.

### Main Results

1. **Discovery of significant SNPs**:
- A total of 43 genome - wide significant SNPs were found among all four traits.
- Single - SNP analysis found a SNP cluster on chromosome 19 significantly associated with MY, FY, and PY, and a SNP cluster on chromosome 29 significantly associated with SCS.

2. **The most significant gene window**:
- A gene window on chromosome 19 explained up to 9.6% of the MY phenotypic variation, 8.1% of the FY phenotypic variation, 9.1% of the PY phenotypic variation, and 1% of the SCS phenotypic variation.

3. **Annotation of candidate genes**:
- Significant SNPs were mainly located in introns (45%), intergenic regions (36%), and within 0 - 5 kb upstream or downstream of the nearest gene (14%).

### Conclusion

This study confirmed the correlation between the quantitative trait loci (QTL) on chromosome 19 and milk production traits, which is consistent with previous research results in other dairy goat populations. The research results provide valuable genetic information for improving the milk production of New Zealand dairy goats through genomic selection.

### Formula Representation

The statistical tests and model formulas involved in the paper are as follows:

- **Single - SNP GWAS model**:
\[
Y_i=\beta_0 + \beta_j X_{ij}+e_i
\]
where \(Y_i\) is the adjusted phenotype of individual \(i\), \(\beta_0\) is the intercept, \(\beta_j\) is the effect of the \(j\) - th SNP, \(X_{ij}\) is the genotype coding of individual \(i\) at the \(j\) - th SNP, and \(e_i\) is the residual.

- **BayesC model**:
\[
Y_i = \mu+\sum_{j = 1}^{p}u_j x_{ij}+e_i
\]
where \(\mu\) is the overall mean, \(u_j\) is the effect of the \(j\) - th SNP, \(x_{ij}\) is the genotype coding of individual \(i\) at the \(j\) - th SNP, and \(e_i\) is the residual.

Through these methods, the researchers successfully identified the genetic regions significantly associated with milk production traits in New Zealand dairy goats, providing an important basis for further genomic selection and breeding programs.