Next Article in Journal
Oxytetracycline and Florfenicol Concentrations in Food-Additive Premixes Authorised for Broiler Chickens: Assessing Degree of Agreement with Manufacturers Labelling
Next Article in Special Issue
Does Small Ruminant Lentivirus Infection in Goats Predispose to Bacterial Infection of the Mammary Gland? A Preliminary Study
Previous Article in Journal
Estrus Detection Using Background Image Subtraction Technique in Tie-Stalled Cows
Previous Article in Special Issue
Genome-Wide SNP Analysis for Milk Performance Traits in Indigenous Sheep: A Case Study in the Egyptian Barki Sheep
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Candidate Gene Association Study for Economically Important Traits in Czech Dairy Goat Breeds

by
Michaela Brzáková
1,*,
Jana Rychtářová
2,
Jindřich Čítek
3,4 and
Zuzana Sztankóová
1
1
Department of Genetics and Breeding of Farm Animals, Institute of Animal Science, 104 00 Prague, Czech Republic
2
Department of Biology of Reproduction, Institute of Animal Science, 104 00 Prague, Czech Republic
3
Department of Genetics and Agricultural Biotechnologies, Faculty of Agriculture, University of South Bohemia, 370 05 Ceske Budejovice, Czech Republic
4
Department of Infectious Diseases and Preventive Medicine, Veterinary Research Institute, 621 00 Brno, Czech Republic
*
Author to whom correspondence should be addressed.
Animals 2021, 11(6), 1796; https://doi.org/10.3390/ani11061796
Submission received: 6 May 2021 / Revised: 9 June 2021 / Accepted: 10 June 2021 / Published: 16 June 2021

Abstract

:

Simple Summary

The milk production traits of goats are economically important. In the Czech Republic, goat milk is processed directly on farms and distributed as cheese. Although goat breeding is not a main focus of animal production in the Czech Republic, it is essential for the agricultural sector. A group of 14 SNPs (single-nucleotide polymorphisms) within four candidate genes (ACACA, BTN1A1, LPL, and SCD) were analysed in two Czech dairy goat breeds, White Shorthaired (WSH) goats and Brown Shorthaired (BSH) goats. The SNPs were significantly associated with milk-production traits (daily milk yield, protein, and fat percentage) and somatic cell count. This information may be useful for marker-assisted selection or related techniques to increase the accuracy of selection.

Abstract

Milk production is influenced by many factors, including genetic and environmental factors and their interactions. Animal health, especially udder health, is usually evaluated by the number of somatic cells. The present study described the effect of polymorphisms in the ACACA, BTN1A1, LPL, and SCD genes on the daily milk yield, fat, and protein percentages and somatic cell count. In this study, 590 White Shorthaired (WSH) and Brown Shorthaired (BSH) goats were included. SNP genotyping was performed by PCR-RFLP and multiplex PCR followed by SNaPshot minisequencing analysis. The linear mixed model with repeated measurement was used to identify the genetic associations between the studied genes/SNPs and chosen traits. All selected genes were polymorphic in the tested goat populations and showed significant associations with milk traits. Only BTN1A1 (SNP g.599 A > G) showed a significant association with the somatic cell score. After Bonferroni correction, a significant effect of LPL g.300G > A on daily milk yield and fat percentage, LPL g.185G > T on protein percentage, and LPL G50C, SCD EX3_15G > A, and SCD EX3_68A > G on fat percentage was found. The importance of environmental factors, such as the herd-year effect, month of milking, and lactation order on all milk performance indicators was confirmed.

1. Introduction

Milk and dairy products are essential for human nutrition. Goat milk is rich in minerals, vitamins, and bioactive components, is easily digestible, and contains fewer allergic proteins than cow milk. These characteristics also suggest the possible use of goat milk for therapeutic purposes [1,2]. In Europe, the small ruminant milk industry is not widespread because of the low number of animals and insufficient milk volume in goats compared to cows. However, the number of goats used for milk production is growing due to expanding demand. In many countries, somatic cell count (SCC), as an indicator of milk quality, plays a role in the milk industry [3,4,5]. A high SCC also negatively affects some flavour characteristics of cheese and ice cream due to a more intense caprine flavour [5,6,7].
Many studies have investigated the relationship between genetic background and milk performance [8,9,10,11,12]. Like cows, the heritabilities of milk yield in goats have been found to be low to moderate. E.g., a study in New Zealand (64,604 lactation records from 23,583 does distributed in 21 flocks) estimates heritabilities of 0.25 for milk yield, 0.24 for fat yield, 0.24 for protein yield, and 0.21 for SCS, suggesting the presence of useful heritable variation [8]. Multiple trait selections for these traits could improve the milk revenue produced from successive generations of New Zealand dairy goats. To a similar extent, a study in Germany estimated heritabilities of 0.15–0.31 for milk yield, 0.21–0.34 for fat content, 0.26–0.50 for protein content, and 0.10–0.17 for the persistence of milk yield [9]. The phenotypic and additive genetic correlations between the milk yield persistence and milk yield in kg were highly positive (0.52–0.72); similarly, the correlations between the protein and fat content were 0.45–0.55. The phenotypic correlations between the fat and protein content and the milk yield were negative, −0.13 to −0.26 and −0.21 to −0.36, respectively (n = 16,579 goats, 42,973 lactations). Others reported milk yield heritabilities of 0.10 to 0.29 (n = 529 goats, 15,509 milk yield test-day records) and a direct heritability for protein percentage of 0.441 (518 phenotypic records from the progeny of 48 sires and 131 dams) [10,11]. The genetic correlation between milk production and SCS can range broadly from −0.16 to 0.43, with a large standard error [12]. Non-genetic factors have been analysed as well. In Alpine goats in Croatia, the herd explained 24% of the variability in daily milk yield, 12% of that in fat content, and 9% of that in protein content; values for the herd test day were 17%, 29% and 30%, and those for the permanent environment were 16%, 3% and 5%, respectively [13].
The effect of major genes is the main research focus, similar to studies in other milking species. The consideration of the αS1-casein genotype may improve the model’s efficiency, translating into more accurate genetic parameters and breeding values [14]. In another study, 48 SNPs in αS1, αS2, β, and κ casein genes were included in the evaluation of genetic parameters [15]. Including genetic effects and relationships among these heritable biomarkers may improve the model efficiency, genetic parameters, and breeding values for milk yield and composition; this inclusion could also help optimise selection practices and profitability for components where technological application may be especially relevant for the cheese-making dairy sector [15]. In addition to the study of candidate loci, the genomic approach in the genome-wide association study has also been applied for the detection of genetic regions of interest [16]. They found a total of 43 genome-wide significant SNPs for lactation yields of milk (MY), fat (FY), protein (PY), and somatic cell score (SCS). A cluster of variants on chromosome 19 associated with MY, FY, and PY was identified, and another cluster on chromosome 29 associated with SCS. 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 [16].
Previous studies demonstrated that the biochemical pathways in mammary glands related to the biosynthesis and secretion of lipids, lactose, and proteins are regulated by complex gene networks [17,18]. Due to the significance of fatty acid biosynthesis, polymorphisms in the following genes were the focus of this study: butyrophilin (BTN1A1); acetyl-CoA carboxylase α (ACACA); lipoprotein lipase (LPL); and stearoyl-CoA desaturase (SCD). BTN1A1 is a milk fat globule membrane protein that plays a crucial role in lipid secretion and milk production [19]. The BTN1A1 gene is located at chromosome 23, has eight exons and seven coding exons, respectively, and a transcription length of 3256 bps. The gene codes for 526 amino acid (AA) residues (NCBI gene 100860762). In the transcript, 219 variant alleles were described. ACACA is the primary regulatory enzyme of fatty acid biosynthesis; it catalyses acetyl-CoA conversion to malonyl-CoA [18,19]. Fatty acids are essential for forming cell membranes and are used to synthesize fat for storage in adipose tissue or secretion into milk by the mammary gland [20]. The ACACA gene is located at chromosome 19, consists of 52 exons, the transcription length is 6990 bps, and codes for 2329 AA residues (NCBI gene 100861224). In all, 3555 variant alleles were found in the transcript.
The SCD gene, which is located on goat chromosome 26, has an important role in the cellular biosynthesis of monounsaturated fatty acids (MUFAs), because most of the conjugated linoleic acids are synthesised in the mammary gland by the action of SCD in circularizing vaccenic acid [21,22]. The SCD mRNA has been identified by Bernhard et al., 2001 [23] and Yahyaoui et al., 2003 [24]. There are many SNPs that have been described and identified in exon 3, intron 3, intron 4, exon 5 and 6, and a deletion of a nucleotide triplet in the 3′UTR [23,24,25,26]. The LPL gene is involved in the hydrolysation of triglycerides to glycerol and free fatty acids and in lipoprotein transportation [27,28]. It is synthesised in the mammary gland’s epithelial cells and influences the release of fatty acids in the mammary tissue [29]. LPL enzymes are encoded by the LPL gene, which consists of nine exons and eight introns, for a total of 3555 nucleotides. This gene encodes a protein containing 478 amino acids, JQ670882. A few SNPs have been described in goats: a missense mutation responsible for a Ser17Thr amino acid substitution at position 17 of the signal peptide (DQ370053:c.G50C), a C2094T (DQ370053) substitution in the 3′UTR of the gene, and a substitution ss522928251:C > T in intron 7 [21,30].
Our aim was to perform an association analysis of SNPs in the ACACA, BTN1A1, LPL, and SCD genes with milk production traits: daily milk yield (DMY), protein, and fat percentage (PP, FP), and SCC of goats on organic farms in the Czech Republic.

2. Materials and Methods

2.1. Ethical Approval

The experiment was carried out under Directive 2010/63/EU of the European Parliament and European Council of 22 September 2010, on the protection of animals used for scientific purposes.

2.2. Animals

In this study, a total of 590 animals were included. All individuals belonged to the two Czech national goat breeds: White Shorthaired (WSH) and Brown Shorthaired (BSH) goats or their crosses with other breeds. Thus, three breed groups were determined: purebred WSH, purebred BSH, and crossbreds of both breeds, where the proportion of WSH or BSH was 50% or higher. The experiment was performed on two farms that were located in the Czech Republic. Both farms specialize in organic goat milk production, which means that goats grazed and were fed only organic feed; no hormones, antibiotics, or similar substances were applied (except a form of veterinary treatment for individuals), no genetically modified organisms were included, and animals were kept under welfare conditions according to legislation in the European Commission Regulation 889/2008, and the European Commission Regulation 834/2007. The winter feed ration consisted of haylage at approximately 2 kg a day, hay ad libitum, and a grain mix, which was dosed during milking in the milking parlour with a total amount of 300 g a day. The summer feed consisted of grass at approximately 2 kg a day (loaded in the stable), hay ad libitum, and grain mix at 300 g a day. Goats were machine-milked twice a day.

2.3. Performance-Testing Database

Phenotype data of genotyped animals were obtained from the performance-testing database of the Czech-Moravian Breeders’ Association. In this study, the daily milk yield (DMY) in kg, milk protein (PP) and fat (FP) per cent, and somatic cell count (SCC) were considered. The analysis of DMY, PP, and FP was performed in the group of 590 goats belonging to the White Shorthair breed (n = 490), Brown Shorthair breed (n = 76), and crossbreeds between WSH, BSH, and Saanen (n = 24). Purebred individuals were approximately 75–100% pure. The goats were sired by 37 bucks; the average number of daughters per buck was 8 (minimum 1, maximum 35). They were on the first to eleventh lactation. In all, 8640 milking records from two farms were analysed. At farm A, 2241 repeated-milk records from 279 dairy goats in 2013, 2014, and 2016 were collected. Milk records from farm B included 311 dairy goats with 6399 milk records between 2010 and 2017. Milk samples were collected repeatedly during the milking seasons. There were 2 to 49 repeated records per goat, on average, there were 14 records (approximately four milk controls per year). Milk samples for DMY, PP and FP analysis were collected throughout the whole year as follows: January (n = 571), February (n = 533), March (n = 735), April (n = 1061), May (n = 1161), June (n = 1203), July (n = 1035), August (n = 1069), September (n = 950), October (n = 74), November (n = 57), and December (n = 191).
Samples from only one farm were analysed for somatic cell count, n = 146 goats, that is, White Shorthair goats n = 100, Brown Shorthair goats n = 38, and crossbreeds n = 8. The goats were sired by 27 bucks, and the average number of daughters was 5 (minimum 1, maximum 19). Milk samples were taken during the third to eleventh lactation. The number of repeated records for SCC was 857. Data were collected during 2016 (n = 728) and 2017 (n = 129). The analysis was conducted repeatedly throughout lactation. There were 2–11 repeated records per goat, on average 5 samples per goat, with approximately three controls per year. Milk samples for SCC analysis were collected throughout the milking season as follows: March (n = 26), April (n = 148), May (n = 148), June (n = 128), July (n = 163), August (n = 124), and September (n = 120).

2.4. DNA Extraction and SNP Genotyping

Blood samples from 590 animals (5 mL of each) were collected from the jugular vein and preserved in 0.5 mM EDTA (pH 8.0). Genomic DNA was extracted from blood using GeneAll, ExgeneTM, and a Clinic SV mini isolation kit (GeneAll Biotechnology cp., Ltd., Seoul, Korea; Bohemia Genetic Ltd., Prague, Czech Republic) according to the manufacturer’s recommendations. The SNPs analysed in our study are described in Table 1. SNPs in the BTN1A1 and LPL genes were genotyped according to previously described methods [30,31,32].
Individual SNPs of the ACACA and SCD loci were detected by primer extension analysis with the SNaPshot Commercial Kit (Applied Biosystem, Foster City, CA, USA). Primers used for the PCR, extension analysis and electropherogram of the SNaPshot product along with the GeneScanTM 120LIZTM size standard are given in Supplementary Table S1.
For the PCR reaction (512 bp fragment) of 3 SNPs in the 5′UTR [33], of ACACA locus, we used the set of primers, designed on the basis of the ovine sequence AJ292286 [34]. The PCR assay was performed in a 10 µL reaction mixture, consisting of 2 µL genomic DNA (10–100 ng), 1× PPP Master Mix (Top Bio Ltd., Prague, Czech Republic), and 0.5 µL (10 pmol/µL; 0.01 mM) of each primer (GENERI Biotech Ltd., Prague, Czech Republic), and sterile water up to volume. Thermal cycling conditions are presented in Supplementary Table S2.
A 536bp fragment of the SCD locus, at region exon3 and intron 3, was amplified by PCR with the following set of primers: F: 5′-TCCTAAgCTTATTCCAgCCCC-3′and R: 5′-gCCAgTCACTCAgAAgTACCC-3′, designed on the basis of the GenBank goat sequence (GenBank AH011188.2; AF422168.1) using Primer 3 software [35]. PCR assay was performed in a 20 µL reaction mixture consisting of 2 µL genomic DNA (10–100 ng), 1× PPP Master Mix (Top Bio Ltd., Prague, Czech Republic), of each primer (GENERI Biotech Ltd., Prague, Czech Republic) and sterile water up to volume. Thermal cycling conditions are presented in Supplementary Table S2.
The presence of fragments obtained in this phase was confirmed by gel-electrophoresis stained with ethidium bromide. The obtained PCR products were purified by using 1 unit of FastAP Thermosensitive Alkaline Phosphatase and Exonuclease I (Fermentas, Ltd., Prague, Czech Republic) to remove unwanted subproducts and incubated at 37 °C for 60 min, followed by 15 min at 85 °C.
The PEA assay utilises internal unlabelled primers which bind to a complementary PCR-generated template in the presence of AmpliTaq DNA Polymerase and fluorescently labelled ddNTPs. The polymerase extends the primer one nucleotide, adding a single ddNTP to its 3′ end. Primers were designed to allow size and colour discrimination between the different alleles (Table S1) and were optimised to be used simultaneously.
The single-base extension (SBE) reaction for ACACA locus was performed in a reaction mixture with final volume of 5.0 µL, containing 1.5 µL of purified multiplex PCR product, 1× extension primer mixture (0.01 mM concentrations): K-ACACA (1206) = 0.5 µL, K-ACACA (1255) = 0.5 µL, K-ACACA (1322) = 0.5 µL, 1.5 µL deionized water, and 2.0 µL of SNaPshot Multiplex Ready Reaction Mix (Applied Biosystems, Foster City, CA, USA). The single-base extension (SBE) reaction for SCD locus in positions EX3_15G > A, IVS3+46 C > T, IVS3+55A > G, EX3_68A > G and IVS3+105A > G was performed in a reaction mixture with a final volume of 6.0 µL, containing 1.5 µL of purified multiplex PCR product, 1× extension primer mixture (0.01 mM concentrations): K-EX3_15G > A = 0.5 µL, K-EX3_68A > G = 0.5 µL, K- IVS3+46C > T = 1.0 µL, K-IVS3 + 55A > G = 0.5 µL K-IVS3 + 105A > G = 0.5 µL, 0.7 µL deionized water, and 2.3 µL of SNaPshot Multiplex Ready Reaction Mix (Applied Biosystems, Foster City, CA, USA). Thermal cycling consisted of 25 cycles of denaturation at 96 °C for 10 s, primer annealing at 50 °C for 5 s, and primer extension at 60 °C for 30 s (Biometra Thermoblock: 050-801 TGradient 96, Biometra, Goettingen, Germany).
For electrophoretic detection, 0.5 µL of purified multiplex SBE reaction was mixed with 0.5 µL of GeneScan-120 LIZ size standard (Applied Biosystems, Foster City, CA, USA) and 9.0 µL of Hi–DiTMFormamide (Applied Biosystems, Foster City, CA, USA), following denaturation step at 95 °C for 5 min and analysed by capillary electrophoresis using the Applied Biosystems® 3130 Genetic Analyzer, an E5-Matrix Standard Set DS-02, a 36 cm capillary, and POP7 polymer (Applied Biosystems, Foster City, CA, USA). The results of genotyping were analysed and evaluated using GeneMapper v 3.5 software (Applied Biosystems, Foster City, CA, USA). The dye colour of the fragment was used to identify the nucleotide of interest (Figure S1).

2.5. Statistical Analysis

The dataset was edited, and unreliable data were removed. SCCs less than 13 and more than 9998 were removed from the analysis. To achieve an approximately normal distribution, the SCC was log-transformed into somatic cell score (SCS). The transformation was performed as follows:
SCS = log2 (SCC/100) + 3
where SCC is somatic cell count, which is expressed in thousands per 1 mL of milk.
Hardy–Weinberg equilibrium was tested by the χ2 test. The effect of gene polymorphisms on milk performance traits and SCS was analysed using the PROC MIXED procedure of SAS with repeated measurements [36]. Tested effects were considered statistically significant at p < 0.05, but biological importance was also considered. The following linear model was used for all traits (DMY, PP, FP, SCS):
Yijklmnop = µ + Gi + HYj + monthk + lacl + breedm + siren + goato + eijklmno ,
where Yijklmno = DMY, FP, PP, SCS; Gi = fixed effect of the genotype (class effect i = 1, 2, 3); HYj = combined fixed effect of herd-year (class effect j = 1, …, 11); monthk = fixed effect of the month of the year of milking (class effect l = 1, …, 12); lacl = fixed effect of the lactation order (class effect l = 1, …, 9 for DMY, FP and PP or l =1, …, 9 for SCS); breedm = fixed effect of the breed (class effect m = 1, …, 3); siren = random effect of the father of the goat; goato= permanent environment of the goat (repeated measurement); and eijklmno = random residual effect.
The post hoc comparison was performed by Scheffe´s method. A Bonferroni correction for multiple comparisons was applied to all significant associations. The correction factor was derived from the number of SNPs tested. The significance threshold (p < 0.05 and p < 0.01) was divided by the number of tests. Thus, Bonferroni-corrected significance levels of 0.05/13 = 0.004 and 0.01/13 = 0.0008 were applied.

3. Results and Discussion

3.1. Descriptive Statistics and Phenotypic Correlations

Genotype and allelic frequencies and the number of animals and records included in the analysis are shown in Table 2.
The descriptive statistics of the analysed traits are shown in Table 3.
Phenotypic correlations between milk traits (DMY, PP, FP) and SCS are shown in Table 4. The relationships of DMY with PP and FP were negative, which corresponds to the well-known dilution effect reflecting the reductions in fat and protein contents as milk yield increases [37]. Among milk components and SCC, positive correlations were found. Some authors found a negative and significant phenotypic correlation between the logarithm of somatic cell count and milk yield [38]. Milk protein content consistently showed a significant positive correlation to the logarithm of SCC. Their study showed a similar correlation between SCC and milk yield, or milk protein content of dairy goats’ milk as found in dairy cows’ milk. However, they stated the impossibility of employing commonly used physiological parameters for dairy cows in evaluating the mammary health status of dairy goats. According to other results, the goat milk composition did not change when milk SCC varied among three groups from 214,000 to 1450,000 cells/mL [39]. In ovine milk, the components can be expected to vary independently of milk SCC [40].

3.2. Environmental Factors

Environmental factors play an essential role in milk production and mastitis occurrence, so controlling environmental factors could permit or prevent animals from expressing their genetic potential [15]. In our analysis, we considered the following fixed effects: the combined effect of the herd-year (HY), month of milking, lactation order, and breed. All of the mentioned factors were significant for daily milk yield, protein, and fat percentages. The exception was an effect of the breed on DMY, where significance was observed only in models with a few SNPs LPL g.103G > A, g.185G > T, g.257C > T, and g.300G > A. The combined effect of HY explained 7% of the total of 26% variability explained by all tested factors. The HY effect included farm management, milking routines, milking frequency, and feed quality. The Czech Republic is situated in the middle of Europe. The climate is mild with four seasons (spring, summer, autumn, winter). Seasonal kidding usually starts in January in the winter season. After this, goat milk production increases with the growing needs of kids. The month of milking reflects changes during time as well as the nutritional condition of pasture which is rising with the increasing temperature, and also the variation of climate and phase of lactation throughout the year. Milk yield is also negatively correlated with milk fat and protein contents, so with increasing milk production, the milk fat and protein levels decrease [41]. The fixed effect of a month of milking explained approximately 10% of a total of 26% variability. The highest milk production was observed between the lactations. This is in agreement with many authors [41,42,43]. During this period, the highest milk yield is probably caused by the physical appearance, size, and quality of the udder. These morphological characteristics are affected by the breed and genotype of goats as well [44]. Contrary, younger goats tend to have a higher milk fat content than older goats [44]. Milk components were significantly affected by breed; for daily milk yield, the impact of breed varied. Differences between breeds were also confirmed by many authors [42,44].
The somatic cell count could be affected by numerous factors, such as milking routine, stage of lactation, lactation order, breed, feed quality, and health status [3,5,45]. In our analysis of the SCS, only the HY effect was significant, but not the effect of month of milking, lactation order, and breed. HY comprises farm management, feed quality, milking frequency, milking hygiene, pasture quality throughout the year, etc. This analysis was performed on only one farm over multiple years. Thus, the significance of the HY effect probably indicates the difference between the tested years. A trend of increasing somatic cell counts was observed throughout the year (from March to September); the differences among succeeding months were significant or near the significance threshold. Bergonier et al., 2003, claimed that higher rates of mastitis occurrence are observed at the beginning of machine milking and during the first third of lactation, but mastitis is rarely observed during drying-off or parturition [46]. The stage of lactation was not included in this study because of a lack of data. Nevertheless, the mastitis incidence in goats does not vary with the lactation stage, in contrast with cows [46]. Goat milk contains naturally higher SCCs than cow milk due to the apocrine secretory process in goats [47]. There is no consensus on whether the SCC is related to milk production (MY, FP, PP). Several authors claimed that a significantly lower fat content was found in goat milk infected with S. aureus than in noninfected milk. Even an SCC of approximately 3,300,000 cells/mL might not be connected with mastitis or pathological differences in the goat mammary gland [5].
Analysed farms were specialized in organic goat milk production, so the environmental conditions and herd management might differ in feeding routine, health management, and welfare issues compared to conventional goat farms. There is an assumption that there should be differences in SCC between organic and conventional farming. Goats from conventional farming could be exposed to various chemical agents such as disinfections, mycotoxins, pesticides, or residues of antibiotics. These agents affect the development and activity of the microbiological profile of milk. However, many studies on small ruminants and dairy cattle did not confirm differences between organic and conventional farms in milk quality parameters [48,49].

3.3. Associations between SNP, Milk Production Traits, and SCS

In total, 14 SNPs were included in the association analysis. SNP g.1255A > G of the ACACA gene was excluded because of monomorphism. The associations between 13 SNPs, milk production traits (MY, FP, PP), and SCS are shown in Table 5, and the differences between genotypes are shown in Table 6.
Polymorphisms other than those in the ACACA gene were associated with fat content in milk. Seven SNPs, namely, BTN1A1 g.599A > G; LPL g.185G > T, g.300G > A, and G50C, and SCD EX3_15G > A; EX3_68A > G, and IVS3+46 C > T, were detected as significant (Table 5). Significant differences between genotypes were found for five SNPs: BTN1A1 g.599A > G; LPL G50C; and SCD EX3_15G > A, EX3_68A > G, and IVS3+46 C > T (Table 6). After Bonferroni correction, only SNPs in the SCD and LPL genes showed a significant association with fat content in milk. The non-significance of ACACA polymorphisms is surprising, as other authors found an effect of the gene, especially on the milk fat content, but not for milk yield and protein content in goat and sheep milk [33,50], which keeps the gene promising for the future research [51]. The role of LPL in fat synthesis was also stated by other authors [52].
All analysed genes showed an association with protein content. A strong association (p < 0.01) was found between BTN1A1 g.599A > G and LPL g.185G > T. Other SNPs, such as ACACA g.1322T > C and g.257C > T and SCD IVS3+46 C > T, reached significant levels of p < 0.05 (Table 5). Significant differences between genotypes were found for ACACA g.1322T > C, BTN1A1 g.599A > G, LPL 257C > T, and SCD IVS3+46 C > T (Table 6). After Bonferroni correction, the number of statistically significant SNPs decreased from 5 SNPs to just LPL g.185G > T. However, for this SNP, the differences between genotypes were not significant. Only two SNPs were found to be associated with daily milk yield, ACACA g.1322T > C (p < 0.05) and LPL g.300G > A (p < 0.01) (Table 5). The differences between genotypes showed only slightly significant differences (Table 5). For the SCS, the only association was detected with the SNP BTN1A1 g.599A > G (p < 0.05). In this polymorphism, genotype AA was connected with the highest somatic cell score in milk, and significant differences were found between the AG and GG genotypes (p < 0.05). Unfortunately, this association was not significant after Bonferroni correction. No other significant associations were found for the SCS.
When comparing the effect of SNPs retained after Bonferroni correction and the differences among genotypes, only the LPL G50C and SCD EX_15G > A and EX3_68A > G polymorphisms showed a significant effect at the SNP and genotype levels. Additionally, other gene polymorphisms in goats are studied. For example, the AGPAT6 gene was found to significantly influence both the fat and protein contents and milk yield [53]. Additionally, the variability of the haplotypic sequences at the loci of the casein complex was studied [15]. The authors point out that a complete definition of the haplotypes in the casein complex in goats is difficult given the high genetic variability.
LPL, ACACA, and SCD ovine genes expression were found to be influenceable by diets [54]. The ability of nutrigenomic regulation of the transcription confirmed that these genes play a critical role in the regulation of lipid metabolism processes in sheep and could be associated with fatty acid profiles in milk and meat. Ovine LPL gene should also be studied due to its expression to microRNA-148a [55]. SCD and BTN1A1 genes were not found to be differentially expressed when comparing two Spanish sheep breeds [56].
Such studies like this are frequent in cattle. Often, the polymorphisms of genes for milk proteins are analysed [57]. The authors found a significant influence of the CSN1S2 gene on the milk protein content and also of the DGAT1 gene on the milk yield. Others found a significant association of the SCD1 bovine gene with fat and milk urea contents and of the ACACA gene with SCS [20]. Other analyses described the significant effect of DGAT1 and ACACA polymorphisms on the milk performance indicators [58,59]. DGAT1 and SCD1 are favourite subjects in dairy cattle [60]. More extensive studies were done, comprising eleven loci and 25 indels [61]. The BTN1A1 polymorphisms were studied even in water buffalo [62].

4. Conclusions

Concisely, some SNPs in the included genes showed association with milk traits but not with SCS. After more accurate Bonferroni correction, the significant associations of SNPs were only rare. Thus, our results support the generally accepted fact that environmental factors are more important than genetic for milk-production traits. However, along with population genetic analyses, the study of major genes in goats helps to better understand the genetic background of the milk-performance complex. This may contribute to future more effective selection in dairy goat populations.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11061796/s1, Table S1: Primers used for the amplification PCR product and extension primers of the goat Stearoyl-CoA desaturase, Acetyl-CoA carboxylase, cDNA, and annealing temperature, Table S2: Thermal cycling conditions of Stearoyl-CoA desaturase and Acetyl-CoA carboxylase. Figure S1: Electropherogram of SNPs in the Stearoyl-coenzyme A desaturase and Acetyl-CoA carboxylase α gene analysed with Gene Mapper software.

Author Contributions

Conceptualisation, M.B., J.R., J.Č., and Z.S.; methodology, J.R., Z.S., and M.B.; formal analysis, M.B. and J.Č.; writing—review and editing, M.B. and J.Č.; project administration, Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Ministry of Agriculture of the Czech Republic, institutional support MZE-RO0718 (V001), by the project NAZV QK1910320, and by the Grant Agency of South Bohemia University project GAJU028/2019/ Z.

Institutional Review Board Statement

All animal experiments were under the supervision of the Institutional Animal Care and Use Committee of the Faculty of Agriculture of the South Bohemia University, where the experiment was carried out approval number 22036/2019-MZE-18134.

Acknowledgments

The authors are grateful for the financial support of the Ministry of Agriculture of the Czech Republic. Additionally, we would like to thank the Sheep and Goat Breeders´ Association and goat farmers for providing the performance-recording data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lad, S.S.; Aparnathi, K.D.; Mehta, B.; Velpula, S. Goat Milk in Human Nutrition and Health—A Review. Int. J. Curr. Microbiol. App. Sci. 2017, 6, 1781–1792. [Google Scholar] [CrossRef] [Green Version]
  2. Sonu, K.S. Compositional and therapeutic signatures of goat milk: A review. Int. J. Chem. Stud. 2020, 8, 1013–1019. [Google Scholar] [CrossRef] [Green Version]
  3. Kuchtík, J.; Králíčková, Š.; Zapletal, D.; Węglarzy, K.; Šustová, K.; Skrzyżala, I. Changes in physico-chemical characteristics, somatic cell count and fatty acid profile of Brown Short-haired goat milk during lactation. Anim. Sci. Pap. Rep. 2015, 33, 71–83. [Google Scholar]
  4. Králíčková, Š.; Kuchtík, J.; Filipčík, R.; Lužová, T.; Šustová, K. Effect of chosen factors on milk yield, basic composition and somatic cell count of organic milk of brown short-haired goats. Acta Univ. Agric. Silvic. Mendel. Brun. 2013, 61, 99–105. [Google Scholar] [CrossRef] [Green Version]
  5. Türkmen, N.; Kanca, H.; Gursoy, A. Effects of Somatic Cell Count in Goat Milk on Some Physical, Chemical and Sensory Properties of Vanilla Ice Cream. Atatürk Üniversitesi Vet. Bil. Derg. 2018, 13, 135–140. [Google Scholar] [CrossRef] [Green Version]
  6. Sánchez-Macías, D.; Morales-delaNuez, A.; Torres, A.; Hernández-Castellano, L.E.; Jiménez-Flores, R.; Castro, N.; Argüello, A. Effects of addition of somatic cells to caprine milk on cheese quality. Int. Dairy J. 2013, 29, 61–67. [Google Scholar] [CrossRef] [Green Version]
  7. Silanikove, N.; Merin, U.; Shapiro, F.; Leitner, G. Subclinical mastitis in goats is associated with upregulation of nitric oxide-derived oxidative stress that causes reduction of milk antioxidative properties and impairment of its quality. J. Dairy Sci. 2014, 97, 3449–3455. [Google Scholar] [CrossRef] [PubMed]
  8. Scholtens, M.R.; Lopez-Villalobos, N.; Garrick, D.; Blair, H.; Lehnert, K.; Snell, R. Genetic parameters for total lactation yields of milk, fat, protein, and somatic cell score in New Zealand dairy goats. Anim. Sci. J. 2020, 91, e13310. [Google Scholar] [CrossRef]
  9. Wolber, M.R.; Hamann, H.; Herold, P. Continuous and extended milking in dairy goats 2. Communication: Genetic analysis of dairy traits. Züchtungskunde 2019, 91, 129–140. [Google Scholar]
  10. Nogueira, B.C.F.; Oliveira, H.R.; Souza, N.O.; Junqueira, V.S.; Rodrigues, M.T.; Silva, F.F.; Brito, L.F. Random regression models for genetic evaluation of milk yield in the second lactation of tropical dairy goats. Livest. Res. Rural Dev. 2020, 32, 26. [Google Scholar]
  11. Verma, M.; Dige, M.S.; Kaushik, R.; Gautam, D.; De, S.; Rout, P.K. Milk composition traits in Jamunapari goats: Genetic parameters and effect of CSN1S1 gene variants. Int. J. Dairy Technol. 2020, 73, 12–21. [Google Scholar] [CrossRef]
  12. Biffani, S.; Tiezzi, F.; Fresi, P.; Stella, A.; Minozzi, G. Genetic parameters of weeping teats in Italian Saanen and Alpine dairy goats and their relationship with milk production and somatic cell score. J. Dairy Sci. 2020, 103, 9167–9176. [Google Scholar] [CrossRef] [PubMed]
  13. Špehar, M.; Mulc, D.; Barać, Z.; Mioč, B.; Kasap, A. Estimation of genetic parameters for dairy traits using repeatability animal models based on test-day data for Alpine goat in Croatia. Mljekarstvo 2019, 69, 215–221. [Google Scholar] [CrossRef]
  14. Inostroza, M.G.P.; Landi, V.; González, F.J.N.; Jurado, J.M.L.; Bermejo, J.V.D.; Álvarez, J.F.; Martínez, M.D.A. Integrating Casein Complex SNPs Additive, Dominance and Epistatic Effects on Genetic Parameters and Breeding Values Estimation for Murciano-Granadina Goat Milk Yield and Components. Genes 2020, 11, 309. [Google Scholar] [CrossRef] [Green Version]
  15. Inostroza, M.G.P.; González, F.J.N.; Landi, V.; Jurado, J.M.L.; Bermejo, J.V.D.; Álvarez, J.F.; Martínez, M.D.A. Bayesian Analysis of the Association between Casein Complex Haplotype Variants and Milk Yield, Composition, and Curve Shape Parameters in Murciano-Granadina Goats. Animals 2020, 10, 1845. [Google Scholar] [CrossRef]
  16. Scholtens, M.; Jiang, A.; Smith, A.; Littlejohn, M.; Lehnert, K.; Snell, R., II; Lopez-Villalobos, N.; Garrick, D.; Blair, H. Genome-wide association studies of lactation yields of milk, fat, protein and somatic cell score in New Zealand dairy goats. J. Anim. Sci. Biotechnol. 2020, 11, 55. [Google Scholar] [CrossRef]
  17. Andres, A.C.; Djonov, V. The mammary gland vasculature revisited. J. Mammary Gland Biol. Neoplasia 2010, 15, 319–328. [Google Scholar] [CrossRef] [Green Version]
  18. Rudolph, M.C.; Neville, M.C.; Anderson, S.M. Lipid synthesis in lactation: Diet and the fatty acid switch. J. Mammary Gland Biol. Neoplasia 2007, 12, 269–281. [Google Scholar] [CrossRef]
  19. Qu, Y.; Liu, Y.; Ma, L.; Sweeney, S.; Lan, X.; Chen, Z.; Li, Z.; Lei, C.; Chen, H. Novel SNPs of butyrophilin (BTN1A1) and milk fat globule epidermal growth factor (EGF) 8 (MFG-E8) are associated with milk traits in dairy goat. Mol. Biol. Rep. 2011, 38, 371–377. [Google Scholar] [CrossRef]
  20. Cecchinato, A.; Ribeca, C.; Chessa, S.; Cipolat-Gotet, C.; Maretto, F.; Casellas, J.; Bittante, G. Candidate gene association analysis for milk yield, composition, urea nitrogen and somatic cell scores in Brown Swiss cows. Animal 2014, 8, 1062–1070. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Crepaldi, P.; Nicoloso, L.; Coizet, B.; Milanesi, E.; Pagnacco, G.; Fresi, P.; Dimauro, C. Association of acetyl-coenzyme A carboxylase α, stearoyl-coenzyme A desaturase, and lipoprotein lipase genes with dairy traits in Alpine goats. J. Dairy Sci. 2013, 96, 1856–1864. [Google Scholar] [CrossRef]
  22. Bauman, D.E.; Mather, I.H.; Wall, R.J.; Lock, A.L. Major advances associated with the biosynthesis of milk. J. Dairy Sci. 2006, 89, 1235–1243. [Google Scholar] [CrossRef]
  23. Bernard, L.; Leroux, C.; Hayeds, H.; Gautier, M.; Chilliard, Y.; Martin, P. Characterization of the caprine stearoyl-CoA desaturase gene and its mRNA showing an unusually long 3′-UTR sequence arising from a single exon. Gene 2001, 281, 53–61. [Google Scholar] [CrossRef]
  24. Yahyaoui, M.H.; Vaiman, D.; Sánchez, A.; Folch, J.M. Mapping of the goat stearoyl coenzyme A desaturase gene to chromosome 26. Anim. Genet. 2003, 34, 474. [Google Scholar] [CrossRef] [PubMed]
  25. Chen, Z.; Sun, J.; Li, Z.; Lan, X.; Zhang, C.; Qu, Y.; Liu, Y.; Fang, X.; Lei, C.; Chen, H. Novel SNPs in the caprine stearoyl-CoA desaturase (SCD) and decorin (DCN) genes that are associated with growth traits in Chine goat breeds. Mol. Biol. Rep. 2011, 38, 3121–3127. [Google Scholar] [CrossRef]
  26. Zhang, C.-L.; Gao, X.-Y.; Shao, R.-Y.; Wang, Y.-H.; Fang, X.-T.; Chen, H. Stearoyl-Co A Desaturase (SCD) gene polymorphism in goat breeds. Biochem. Genet. 2010, 48, 822–828. [Google Scholar] [CrossRef]
  27. Crisa, A.; Marchitelli, C.; Pariset, L.; Gontarini, G.; Signarelli, F.; Napolitano, R.; Catillo, G.; Valentini, A.; Moili, 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] [PubMed] [Green Version]
  28. Hidayati, C.; Sumantri, R.; Noor, R.; Priyanto, R.; Rahayu, S. Single nucleotide polymorphisms of lipoprotein lipase gene and its association with marbling quality in local sheeps. J. Indones. Trop. Anim. Agric. 2015, 40, 1–10. [Google Scholar] [CrossRef]
  29. Moioli, B.; Contarini, G.; Avalli, A.; Gatillo, G.; Orru, L.; De Matteis, G.; Masoero, G.; Napolitano, F. Short communication: Effect of stearoyl-coenzyme A desaturase polymorphism on fatty acid composition of milk. J. Dairy Sci. 2007, 90, 3553–3558. [Google Scholar] [CrossRef] [PubMed]
  30. Badaoui, B.; Serradilla, J.M.; Tomàs, A.; Urrutia, B.; Ares, J.L.; Carrizosa, J.; Sànchez, A.; Jordana, J.; Amills, M. Identification of two polymorphisms in the goat lipoprotein lipase gene and their association with milk production traits. J. Dairy Sci. 2007, 90, 3012–3017. [Google Scholar] [CrossRef]
  31. Kyselova, J.; Sztankoova, Z.; Svitakova, A.; Melcova, S.; Krejcova, M. Effect of single nucleotide polymorphisms of BTN1A1 gene on the milk production traits of Czech dairy goats. In Proceedings of the 34th International Society for Animal Genetics Conference—ISAG 34th, Xi’an, China, 27 July–1 August 2014. [Google Scholar]
  32. Rychtářová, J.; Sztankóová, Z.; Hofmannová, M.; Vostrý, L.; Milerski, M. Characterization of the new genetic variant in the caprine lipoprotein lipase gene. Small Rumin. Res. 2020, 182, 5–10. [Google Scholar] [CrossRef]
  33. Signorrelli, F.; Napolitano, F.; De Mateis, G.; Scatà, M.C.; Gennaro, C.; Tripaldi, C.; Moioli, B. Identification of novel single nucleotide polymorphisms in promotor III of the Acetyl-CoA carboxylase-α gene in goats affecting milk production traits. J. Hered. 2009, 100, 386–389. [Google Scholar]
  34. Sztankóová, Z.; Rychtářová, J.; Kyselová, J.; Czerneková, V. Simultaneous genotyping of 4 SNPs in promoter III of the ovine ACACA. Small. Rum. Res. 2016, 138, 25–30. [Google Scholar] [CrossRef]
  35. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef] [Green Version]
  36. SAS Institute Inc. Base SAS 9.4 Procedures Guide: Statistical Procedures, 2nd ed.; Statistical Analysis System Institute Inc.: Cary, NC, USA, 2013. [Google Scholar]
  37. Pulina, G.; Macciotta, N.; Nudda, A. Milk composition and feeding in the Italian dairy sheep. Ital. J. Anim. Sci. 2010, 4 (Suppl. 1), 5–14. [Google Scholar] [CrossRef] [Green Version]
  38. Ying, C.; Wang, H.T.; Hsu, J.T. Relationship of somatic cell count, physical, chemical and enzymatic properties to the bacterial standard plate count in dairy goat milk. Livest. Prod. Sci. 2002, 74, 63–77. [Google Scholar] [CrossRef]
  39. Chen, S.; Wang, J.; Kessel, J.; Ren, F.; Zeng, S. Effect of somatic cell count in goat milk on yield, sensory quality and fatty acid profile of semisoft cheese. J. Dairy Sci. 2010, 93, 1345–1354. [Google Scholar] [CrossRef] [Green Version]
  40. Jaeggi, J.J.; Govindasamy-Lucey, S.; Berger, Y.M.; Johnson, M.E.; McKusick, B.C.; Thomas, D.L.; Wendorff, W.L. Hard ewe’s milk cheese manufactured from milk of three different groups of somatic cell counts. J. Dairy Sci. 2003, 86, 3082–3089. [Google Scholar] [CrossRef]
  41. Goetsch, A.L.; Zeng, S.S.; Gipson, T.A. Factors affecting goat milk production and quality. Small Rumin. Res. 2011, 101, 55–63. [Google Scholar] [CrossRef]
  42. Vacca, G.M.; Stocco, G.; Dettori, M.L.; Pira, E.; Bittante, G.; Pazzola, M. Milk yield, quality, and coagulation properties of 6 breeds of goats: Environmental and individual variability. J. Dairy. Sci. 2018, 101, 7236–7247. [Google Scholar] [CrossRef]
  43. Ciappesoni, G.; Přibyl, J. Factors affecting goat milk yield and its composition. Czech J. Anim. Sci. 2004, 11, 465–473. [Google Scholar] [CrossRef] [Green Version]
  44. Idowu, S.T.; Adewumi, O.O. Genetic and non-genetic factors affecting yield and milk composition in goats. J. Adv. Dairy. Res. 2017, 5, 2. [Google Scholar] [CrossRef] [Green Version]
  45. Csanádi, J.; Fenyvessy, J.; Bohata, S. Somatic cell count of milk from different goat breeds. Acta Univ. Sapientiae Aliment. 2015, 845–854. [Google Scholar] [CrossRef] [Green Version]
  46. Bergonier, D.; de Cremoux, R.; Rupp, R.; Lagriffoul, G.; Berthelot, X. Mastitis of dairy small ruminants. Vet. Res. 2003, 34, 689–716. [Google Scholar] [CrossRef] [Green Version]
  47. Park, Y.W.; Guo, M. Goat milk Products: Types of products, manufacturing technology, chemical composition, and marketing. In Handbook of Milk of Non-Bovine Mammals; Park, Y.W., Haenlein, G.F.W., Eds.; Blackwell Publishing: Ames, IA, USA; Oxford, UK, 2006; p. 449. [Google Scholar]
  48. Malissiova, E.; Tzora, A.; Katsioulis, A.; Hatzinikou, M.; Tsakalof, A.; Arvanitoyannis, I.S.; Govaris, A.; Hadjichristodoulou, C. Relationship between production conditions and milk gross composition in ewe´s and goat´s organic and conventional farms in central Greece. Dairy Sci. Technol. 2015, 95, 437–450. [Google Scholar] [CrossRef] [Green Version]
  49. Kouřimská, L.; Legarová, V.; Panovská, Z.; Pánek, J. Quality of cow’s milk from organic and conventional farming. Czech J. Food Sci. 2014, 32, 398–405. [Google Scholar] [CrossRef] [Green Version]
  50. Moioli, B.; Scata, M.C.; De Matteis, G.; Annicchiarico, G.; Catillo, G.; Napolitano, F. The ACACA gene is a potential candidate gene for fat content in sheep milk. Anim. Genet. 2013, 44, 601–603. [Google Scholar] [CrossRef]
  51. Di Gerlando, R.; Mastrangelo, S.; Tortorici, L.; Tolone, M.; Sutera, A.M.; Sardina, M.T.; Portolano, B. Full-length sequencing and identification of novel polymorphisms in the ACACA gene of Valle del Belice sheep breed. J. Genet. 2017, 96, 591–597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Zhao, W.-S.; Hu, S.-L.; Yu, K.; Wang, H.; Wang, W.; Loor, J.; Luo, J. Lipoprotein Lipase, Tissue Expression and Effects on Genes Related to Fatty Acid Synthesis in Goat Mammary Epithelial Cells. Int. J. Mol. Sci. 2014, 15, 22757–22771. [Google Scholar] [CrossRef] [Green Version]
  53. He, C.; Wang, C.; Chang, Z.H.; Guo, B.L.; Li, R.; Yue, X.P.; Lan, X.Y.; Chen, H.; Lei, C.Z. AGPAT6 polymorphism and its association with milk traits of dairy goats. Genet. Mol. Res. 2011, 10, 2747–2756. [Google Scholar] [CrossRef]
  54. Ropka-Molik, K.; Knapik, J.; Pieszka, M.; Szmatola, T.; Piórkowska, K. Nutritional modification of SCD, ACACA and LPL gene expressions in different ovine tissues. Arch. Anim. Breed. 2017, 60, 243–250. [Google Scholar] [CrossRef] [Green Version]
  55. Jin, X.; Hao, Z.; Zhao, M.; Shen, J.; Ke, N.; Song, Y.; Qiao, L.; Lu, Y.; Hu, L.; Wu, X.; et al. MicroRNA-148a Regulates the Proliferation and Differentiation of Ovine Preadipocytes by Targeting PTEN. Animals 2021, 11, 820. [Google Scholar] [CrossRef] [PubMed]
  56. Suárez-Vega, A.; Gutiérrez-Gil, B.; Arranz, J.J. Transcriptome expression analysis of candidate milk genes affecting cheese-related traits in 2 sheep breeds. J. Dairy Sci. 2016, 99, 6381–6390. [Google Scholar] [CrossRef]
  57. Ardicli, S.; Soyudal, B.; Samli, H.; Dincel, D.; Balci, F. Effect of STAT1, OLR1, CSN1S1, CSN1S2, and DGAT1 genes on milk yield and composition traits of Holstein breed. R. Bras. Zootec. 2018, 47, e20170247. [Google Scholar] [CrossRef] [Green Version]
  58. Fontanesi, L.; Calo, D.G.; Galimberti, G.; Negrini, R.; Marino, R.; Nardone, A.; Ajmone-Marsan, P.; Russo, V. A candidate gene association study for nine economically important traits in Italian Holstein cattle. Anim. Genet. 2014, 45, 576–580. [Google Scholar] [CrossRef]
  59. Fontanesi, L.; Scotti, E.; Samorè, A.B.; Bagnato, A.; Russo, V. Association of 20 candidate gene markers with milk production and composition traits in sires of Reggiana breed, a local dairy cattle population. Liv. Sci. 2015, 176, 14–21. [Google Scholar] [CrossRef]
  60. Mao, Y.J.; Chen, R.J.; Chang, L.L.; Chen, Y.; Ji, D.J.; Wu, X.X.; Shi, X.K.; Wu, H.T.; Zhang, M.R.; Yang, Z.P.; et al. Effects of SCD1- and DGAT1-genes on production traits of Chinese Holstein cows located in the Delta Region of Yangtze River. Liv. Sci. 2012, 145, 280–286. [Google Scholar] [CrossRef]
  61. Jiang, J.; Liu, L.; Gao, Y.; Shi, L.; Li, Y.; Liang, W.; Sun, D. Determination of genetic associations between indels in 11 candidate genes and milk composition traits in Chinese Holstein population. BMC Genet. 2019, 20, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Wu, C.; Liu, L.; Huo, J.; Li, L.; Miao, Y. Isolation, bioinformatic and tissue expression analysis of a novel water buffalo gene BTN1A1. Buffalo Bull. 2014, 33, 449–461. [Google Scholar]
Table 1. Analysed SNPs and methods used.
Table 1. Analysed SNPs and methods used.
GeneGeneBank Access. No.SNPsRegionAAMethodReference
Acetyl-CoA carboxylase α (ACACA)AJ292286g.1206C > T3′UTR PEA a[33]
AJ292286g.1255A > G5′UTR PEA
AJ292286g.1322T > C3′UTR PEA
Butyrophilin (BTN1A1)NM001285618.1g.599A > GExon4Glu184/LysPCR-RFLP[31]
Lipoprotein lipase (LPL)KP261023g.103G > Asignal peptideGly/ArgPEA[32]
KP261023g.185G > Tintron I PEA[32]
KP261023g.257C > Tintron I PEA[32]
KP261023g.300G > Aintron I PEA[32]
DQ370053G50Csignal peptideSer36/ThrPEA[30]
Stearoyl-coenzyme A desaturase (SCD)AF422168.1EX3_15G > AExon3Val109/MetPEA[25]
AF422168.1EX3_68A > GExon3Arg/ArgPEA[25]
AF422168.1IVS3+46C > TIntron3 PEAPresent work
AF422168.1IVS3+55A > G Intron3 PEA[25]
AF422168.1IVS3+105A > GIntron3dataPEAPresent work
a SNPs of loci were genotyped by using multiplex Primer Extension Analysis (Supplementary Data Table S1, Figure S1), AA—Amino Acid.
Table 2. Genotype and allelic frequencies and the number of animals and records included in the analysis.
Table 2. Genotype and allelic frequencies and the number of animals and records included in the analysis.
GeneGene Bank Access. No.SNPsGenotypesNFrequencyAlleleFrequencyχ2NmilkNscc
Acetyl-CoA carboxylase α (ACACA)AJ292286g.1206C > TCC1810.522C0.722.7543675490
CT820.401 1885288
TT220.077T0.28 47964
Acetyl-CoA carboxylase α (ACACA)AJ292286g.1322T > C CT190.067C0.030.23741049
TT2660.933T0.97 5629793
Butyrophilin (BTN1A1)NM001285618.1g.599A > GAA70.023A0.150.01914014
AG760.249 1553199
GG2220.728G0.85 4559636
Lipoprotein lipase (LPL)KP261023g.103G > AGG2650.892G0.940.3603851494
GA300.101 36148
AA20.007A0.06 4212
Lipoprotein lipase (LPL)KP261023g.185G > TGG2100.707G0.850.3652793352
GT820.276 1414202
TT50.017T0.15 470
Lipoprotein lipase (LPL)KP261023g.257C > TTT110.037T0.220.44414015
CT1090.367 1669241
CC1770.596C0.78 2445298
Lipoprotein lipase (LPL)KP261023g.300G > AGG2190.737G0.860.2573285453
GA700.236 88895
AA80.027A0.14 816
CC30.010C0.120.221889
Lipoprotein lipase (LPL)DQ370053G50CCG680.224 1410210
GG2320.766G0.88 4752631
Stearoyl-coenzyme A desaturase (SCD)AF422168EX3_15G > AAA1820.591A0.760.7913808512
AG1030.334 2049262
GG230.075G0.24 52483
Stearoyl-coenzyme A desaturase (SCD)AF422168EX3_68A > GAA2560.831A0.910.0515393691
AG500.162 948158
GG20.006G0.09 408
Stearoyl-coenzyme A desaturase (SCD)AF422168 CC2340.760C0.862.8234851665
IVS3+46 C > TCT620.201 1333172
TT120.039T0.14 19720
Stearoyl-coenzyme A desaturase (SCD)AF422168 AA230.075A0.1515.899 **52483
IVS3+55A > G AG480.156 899142
GG2370.769G0.85 4958632
Stearoyl-coenzyme A desaturase (SCD)AF422168IVS3+105A > GAA250.081A0.250.77054686
AG1070.347 2127283
GG1760.571G0.75 3708488
N—number of animals in milk traits analysis (the number of animals for SCC analysis was smaller); Nmilk—number of samples with MDY, PP, FP performance; Nscs—number of samples with SCS; ** significant at p < 0.01.
Table 3. Descriptive statistics of goat milking traits of analysed data (mean ± standard deviation).
Table 3. Descriptive statistics of goat milking traits of analysed data (mean ± standard deviation).
TraitNMean ± SDMinMax
Daily milk yield (L)86402.94 ± 1.0040.86.4
Milk fat percentage (%)86403.06 ± 0.6442.05.0
Milk protein percentage (%)86403.02 ± 0.3242.24.2
SCC 18571353.52 ± 1608.24019.09625.0
SCS 28575.87 ± 1.6940.609.59
1 SCC—somatic cell count in thousands per 1 mL of milk, 2 SCS = log2 (SCC/100) + 3, SCC—somatic cell count.
Table 4. Pearson correlation coefficients between all analysed traits.
Table 4. Pearson correlation coefficients between all analysed traits.
TraitDMYPPFP
PP−0.205 ** ± 0.98
FP−0.154 ** ± 0.660.401 ** ± 0.61
SS−0.418 ** ± 1.540.154 ** ± 1.670.113 ** ± 1.68
DMY—daily milk yield; PP—milk protein %; FP—milk fat %; SCS = somatic cell score; ** p < 0.01.
Table 5. Significance of fixed effects for milk traits and somatic cell score—individual SNPs.
Table 5. Significance of fixed effects for milk traits and somatic cell score—individual SNPs.
Daily Milk Yield (L) SCS
SNPsSNPSNP BCHerd-YearMonth of MilkingLactation OrderBreedSNPSNP BCHerd-YearMonth of MilkingLactation OrderBreed
g.1206C > Tnsns******nsnsns**nsnsns
g.1322T > C *ns******nsnsns**nsnsns
g.599A > Gnsns******ns*ns**nsnsns
g.103G > Ansns*******nsns**nsnsns
g.185G > Tnsns*******nsns**nsnsns
g.257C > Tnsns*******nsns**nsnsns
g.300G > A**********nsns**nsnsns
G50Cnsns******nsnsns**nsnsns
EX3_15G > Ansns******nsnsns**nsnsns
EX3_68A > Gnsns******nsnsns**nsnsns
IVS3+46 C > Tnsns******nsnsns**nsnsns
IVS3+55A > G nsns******nsnsns**nsnsns
IVS3+105A > Gnsns******nsnsns**nsnsns
Protein (%) Fat (%)
g.1206C > Tnsns*******nsns********
g.1322T > C *ns*******nsns********
g.599A > G**ns*********ns********
g.103G > Ansns********nsns********
g.185G > T************ns********
g.257C > T*ns*******nsns*******
g.300G > Ansns******************
G50Cnsns******************
EX3_15G > Ansns*****************
EX3_68A > Gnsns******************
IVS3+46C > T*ns*********ns*******
IVS3+55A > G nsns*******nsns********
IVS3+105A > Gnsns*********ns*******
SCS—somatic cell score; SNP BCp-value after Bonferroni correction of SNP in each row are significant as computed in model with effects of individual SNP, herd-year, month, lactation order, and breed * significant at p < 0.05; ** significant at p < 0.01; ns not significant.
Table 6. Daily milk yield, protein, and fat percentage, and somatic cell score (LSM ± SE) and differences among genotypes.
Table 6. Daily milk yield, protein, and fat percentage, and somatic cell score (LSM ± SE) and differences among genotypes.
SNPsGenotypeDMYPPFPSCS
ACACA g.1206C > TCC2.50 ± 0.1123.16 ± 0.0373.23 ± 0.0606.25 ± 0.461
CT2.55 ± 0.1103.14 ± 0.0363.20 ± 0.0596.24 ± 0.448
TT2.50 ± 0.1393.15 ± 0.0433.16 ± 0.0727.16 ± 0.660
ACACA g.1322T > C CT2.75 ± 0.135 a3.09 ± 0.043 a3.19 ± 0.0716.04 ± 0.725
TT2.52 ± 0.105 a3.15 ± 0.035 a3.21 ± 0.0586.36 ± 0.413
BTN1A1 g.599A > GAA2.74 ± 0.1773.13 ± 0.0513.09 ± 0.0887.07 ± 1.087
AG2.51 ± 0.1083.17 ± 0.034 A3.24 ± 0.057 a7.04 ± 0.477 a
GG2.53 ± 0.1063.12 ± 0.034 A3.17 ± 0.056 a6.08 ± 0.418 a
LPL g.103G > AGG2.70 ± 0.1193.13 ± 0.0453.18 ± 0.071 A6.21 ± 0.588
GA3.07 ± 0.1623.10 ± 0.0563.33 ± 0.093 A,B5.92 ± 0.837
AA3.07 ± 0.2673.18 ± 0.0822.90 ± 0.142 B5.99 ± 1.378
LPL g.185G > TGG2.67 ± 0.1153.11 ± 0.0473.14 ± 0.0776.23 ± 0.622
GT2.72 ± 0.1133.14 ± 0.0453.18 ± 0.0746.10 ± 0.632
TTNot est.Not est.Not est.Not est.
LPL g.257C > TTT2.81 ± 1773.27 ± 0.062 A,B3.18 ± 0.1076.67 ± 1.179
CT2.69 ± 0.1123.15 ± 0.046 A3.20 ± 0.077 A6.07 ± 0.601
CC2.69 ± 0.1173.11 ± 0.047 B3.10 ± 0.079 A6.26 ± 0.658
LPL g.300G > AGG2.69 ± 0.1113.13 ± 0.0443.16 ± 0.0746.22 ± 0.596
GA2.70 ± 0.1243.15 ± 0.0483.21 ± 0.0815.94 ± 0.699
AA2.80 ± 0.3003.32 ± 0.0913.08 ± 0.1637.02 ± 1.854
LPL G50CCC2.61 ± 0.2333.05 ± 0.0662.70 ± 0.114 A,B6.24 ± 1.324
CG2.59 ± 0.1163.13 ± 0.0373.22 ± 0.061 B6.53 ± 0.508
GG2.51 ± 0.1053.15 ± 0.0343.19 ± 0.06 A6.33 ± 0.418
SCD EX3_15G > AAA2.52 ± 0.1053.15 ± 0.0343.21 ± 0.057 A6.33 ± 0.415
AG2.53 ± 0.1103.13 ± 0.0353.14 ± 0.059 A6.30 ± 0.494
GG2.38±0.1303.16 ± 0.0403.24 ± 0.0687.00 ± 0.679
SCD EX3_68A > GAA2.51 ± 0.1043.15 ± 0.0343.21 ± 0.056 A6.36 ± 0.414
AG2.53 ± 0.1203.13 ± 0.0383.14 ± 0.062 B6.36 ± 0.592
GG2.18 ± 0.2773.20 ± 0.0783.55 ± 0.136 A,B6.62 ± 1.738
SCD IVS3+46C > TCC2.52 ± 1033.15 ± 0.034 A3.21 ± 0.056 A6.33 ± 0.412
CT2.52 ± 0.1123.15 ± 0.038 B3.15 ± 0.060 A6.59 ± 0.519
TT2.43 ± 0.1573.05 ± 0.047 A,B3.15 ± 0.0817.00 ± 0.916
SCD IVS3+55A > G AA2.35 ± 0.1343.17 ± 0.0403.26 ± 0.0676.92 ± 0.693
AG2.50 ± 0.1253.14 ± 0.0383.18 ± 0.0646.10±0.609
GG2.52 ± 0.1073.15 ± 0.0343.20 ± 0.0566.34 ± 0.410
SCD IVS3+105A > GAA2.33 ± 0.1283.16 ± 0.0403.27 ± 0.067 A6.64 ± 0.665
AG2.51 ± 0.1103.14 ± 0.0353.16 ± 0.059 A6.16 ± 0.493
GG2.54 ± 0.1053.15 ± 0.0343.21 ± 0.0576.38 ± 0.416
DMY—daily milk yield, PP—milk protein percentage, FP—fat protein percentage, SCS—somatic cell score. a,b Differences between genotypes with the same letters in the same column are significant at p < 0.05. A,B Differences between genotypes with the same letters in the same column are significant at p < 0.01.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brzáková, M.; Rychtářová, J.; Čítek, J.; Sztankóová, Z. A Candidate Gene Association Study for Economically Important Traits in Czech Dairy Goat Breeds. Animals 2021, 11, 1796. https://doi.org/10.3390/ani11061796

AMA Style

Brzáková M, Rychtářová J, Čítek J, Sztankóová Z. A Candidate Gene Association Study for Economically Important Traits in Czech Dairy Goat Breeds. Animals. 2021; 11(6):1796. https://doi.org/10.3390/ani11061796

Chicago/Turabian Style

Brzáková, Michaela, Jana Rychtářová, Jindřich Čítek, and Zuzana Sztankóová. 2021. "A Candidate Gene Association Study for Economically Important Traits in Czech Dairy Goat Breeds" Animals 11, no. 6: 1796. https://doi.org/10.3390/ani11061796

APA Style

Brzáková, M., Rychtářová, J., Čítek, J., & Sztankóová, Z. (2021). A Candidate Gene Association Study for Economically Important Traits in Czech Dairy Goat Breeds. Animals, 11(6), 1796. https://doi.org/10.3390/ani11061796

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