Next Article in Journal
Consequences Matter: Compassion in Conservation Means Caring for Individuals, Populations and Species
Previous Article in Journal
Partnering with Taxidermists for Improved Chronic Wasting Disease Surveillance
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 9
Issue 12
10.3390/ani9121114
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats

by
Yanghai Zhang
1,
Ke Wang
1,
Jinwang Liu
2,3,
Haijing Zhu
2,3,
Lei Qu
2,3,
Hong Chen
1,
Xianyong Lan
1,
Chuanying Pan
1,* and
Xiaoyue Song
2,3,*
1
Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, China
2
Shaanxi Provincial Engineering and Technology Research Center of Cashmere Goats, Yulin University, Yulin, Shaanxi 719000, China
3
Life Science Research Center, Yulin University, Yulin, Shaanxi 719000, China
*
Authors to whom correspondence should be addressed.
Animals 2019, 9(12), 1114; https://doi.org/10.3390/ani9121114
Submission received: 15 November 2019 / Revised: 8 December 2019 / Accepted: 9 December 2019 / Published: 11 December 2019
(This article belongs to the Section Small Ruminants)
Browse Figure
Versions Notes

Abstract

:

Simple Summary

In the field of animal breeding, selection using molecular genetic markers, such as insertions/deletions (indels), is a novel method to reduce the generation interval and more accurately improve performance. The casein alpha s1 (CSN1S1) gene is one of the major genes regulating the milk performance of mammals, and it is also relevant to body development. Previous studies have shown that an 11-bp indel in this gene was strongly associated with goat litter size, however, its effects on milk performance and body measurement traits have not been studied. In this study, this indel was genotyped in three Chinese goat breeds, and further analysis results revealed that this indel polymorphism was significantly associated with milk performance in dairy goats, as well as affecting the body measurement traits in all three breeds. These findings hint that the 11-bp indel within CSN1S1 gene could be utilized for effectively selecting goats with excellent performances in terms of production and breeding.

Abstract

The casein alpha s1 (CSN1S1) gene encodes α-s1 casein, one of the proteins constituting milk, which affects milk performance, as well as improving the absorption of calcium and bone development in mammals. A previous study found that an 11-bp insertion/deletion (indel) of this gene strongly affected litter size in goats. However, to our knowledge, the relationships between this polymorphism and the milk performance and body measurement traits of goats have not been reported. In this paper, the previously identified indel has been recognized in three Chinese goat breeds, namely the Guanzhong dairy goat (GZDG; n = 235), Shaanbei white cashmere goat (SBWC; n = 1092), and Hainan black goat (HNBG; n = 278), and the following three genotypes have been studied for all of the breeds: insertion/insertion (II), deletion/deletion (DD), and insertion/deletion (ID). The allele frequencies analyzed signified that the frequencies of the “D” allele were higher (47.8%–65.5%), similar to the previous report, which indicates that this polymorphism is genetically stable in different goat breeds. Further analysis showed that this indel was markedly associated with milk fat content, total solids content, solids-not-fat content, freezing point depression, and acidity in GZDG (p < 0.05), and also affected different body measurement traits in all three breeds (p < 0.05). The goats with II genotypes had superior milk performance, compared with the others; however, goats with DD genotypes had better body measurement sizes. Hence, it may be necessary to select goats with an II or DD genotype, based on the desired traits, while breeding. Our study provides information on the potential impact of the 11-bp indel polymorphism of the CSN1S1 gene for improving the milk performance and body measurement traits in goats.

1. Introduction

The goat industry is one of the oldest animal husbandry industries, with a high potential; however, the actual efficiency of goat breeding is not high, and, at present, the demand for goat products exceeds the supply [1,2,3,4,5]. Thus, a pressing question that has been proposed is how to improve the production performance of goats [6,7]. Traditional breeding, for example, using estimated breeding values (EBVs), has been extremely successful in observing the genetic progress. However, the aim of breeding programs is to predict the genetic merit of the next generation and to improve the genetic performance [8]. Based on the identification of genetic DNA markers, such as single nucleotide polymorphisms (SNPs) and insertions/deletions (Indels), novel breeding methods were built, including marker-assisted selection (MAS) and genomic selection (GS), which would significantly reduce the generation interval and increase accuracy [8,9]. At present, some candidate genes and DNA markers have been identified for improving reproductive and body measurement traits in livestock [10,11,12,13], but more crucial and powerful candidate genes require identification in order to comprehensively improve production traits in goat breeding.
Casein gene family (CSNs) is linked as a CSN1S1–CSN2–CSN1S2–CSN3 gene cluster encoding the following four caseins: αs1, β, αs2, and κ. Milk proteins are constituted of whey proteins and above-mentioned four kinds of caseins, where caseins account for the majority of the milk [14]. Among these casein genes, the casein alpha s1 (CSN1S1) gene encodes α-s1 casein, which is mainly expressed in breast tissues, especially in the lactating breasts of females (https://www.proteinatlas.org/). This type of casein is also the most abundant casein in bovine and caprine milk, and can provide sufficient amino acids for animal growth and development [15]. Therefore, this gene exerts a crucial role in milk yield and composition. In previous studies, the goat CSN1S1 gene has shown at least 18 alleles, including A, B1, B2, B3, B4, B′, C, D, E, F, G, H, I, L, M, N, O1, and O2 [16,17,18], and many studies have shown that these polymorphisms affect caprine milk performance in different goat breeds. The N and F alleles of CSN1S1 have been shown to increase the milk yield in the Xinong Saanen dairy goat [19]; the BB and AB genotypes have been associated with the highest milk protein and fat in Sarda goats [20]; and some polymorphisms have been shown to affect the fat content, protein content, and coagulation properties of goat milk [21,22]. In other ruminants, such as cows [23], sheep [24], and buffalo [25], CSN1S1 polymorphisms have also been shown to significantly affect the milk production traits.
On the other hand, casein has been considered to have a high nutritional value, with the potential to improve bone development and growth (as it contains almost all of the necessary amino acids), as well as increasing the absorption of calcium and other micro-elements in the infant gut by complexing and solubilizing calcium ions [14,15]. Animal experiments in vivo have convincingly testified that casein or casein phosphopeptides (CPPs; originating from the digestion of casein) can promote calcium uptake and bone formation, as well as increase bone strength and quality, and body weight [26,27,28]. Therefore, this gene may play a crucial role in the body size, structure, and development of goats.
Furthermore, CSN1S1 may even play an important role, in terms of reproductive traits. The BMPR1B gene (also called FecB) is a famous reproduction-related gene, as it can significantly regulate the ovulation rate in sheep [29,30]. As far back as 1994, Montgomery et al. determined that the FecB gene is located on chromosome six in sheep, and identified genetic linkages with platelet-derived growth factor receptor-alpha (PDGFRA), CSN1S1, and other casein genes [31]. According to the current genome sequence, the number of bases between the BMPR1B, and PDGFRA and CSN1S1 that differ by 42.7 and 60.5 Mb in sheep (NC_040257.1), respectively; these are 40.0 and 55.6 Mb in goat (NC_030813.1), respectively. This information builds the basis for the conclusions from 1994 [31]. Furthermore, Lan et al. found that the polymorphisms in other casein family genes (CSN3 and CSN1S2) had significant positive effects on the litter size of the Xinong Saanen dairy goat [32]. Therefore, we speculated that the CSN1S1 gene may be linked to the BMPR1B gene in goats, and thus, the variations within CSN1S1 might affect the regulation of fecundity [1].
As a powerful candidate gene in livestock breeding, our previous study found that an 11-bp indel within the CSN1S1 gene strongly affected the litter size in a larger sample of the Shaanbei white cashmere goat (SBWC) [1]. However, the relationships between this polymorphism and the milk performance and body measurement traits of goats have not yet been studied, to our knowledge. Therefore, based on our preliminary work, this study is aimed at unveiling the association of the 11-bp indel of CSN1S1, and the milk and body measurement traits in Chinese goats, which would contribute to selecting a herd with a high productivity using the MAS method in goat breeding.

2. Materials and Methods

2.1. Sample and Data Collection

All of the animal experiments in this study conformed to the relevant laws and guidelines and were approved by the Faculty of Animal Policy and Welfare Committee of Northwest A&F University (NWAFAC1008).
A total of 1605 ear samples (small ear notches) were randomly collected from the following three Chinese indigenous goat breeds: Guanzhong dairy goats (GZDG; n = 235), Shaanbei white cashmere goats (SBWC; n = 1092), and Hainan black goats (HNBG; n = 278). The dairy-use breed GZDG were adult healthy female goats (2–4 years old) from a local farm in Sanyuan county of Shaanxi Province [33]. These goats were unrelated, and the milk performance indices were recorded by Lan et al. in 2013 [33]. Based on a 305-day lactation period, different lactation data were collected from the GZDG [33]. The milk samples were collected twice a day (05:00 and 17:00) when the dairy goats were milked. Then, the milk samples were measured for the following eight milk performance indices: milk fat content (%), protein content (%), total solids (TS) content (%), solids-not-fat (SNF) content (%), lactose content (%), density, freezing point depression (FPD; °C), and acidity (pH) [33]. Each milk performance index was taken as the average of two values (05:00 and 17:00). Additionally, the body measurement traits of the GZDG were collected, including the body height, body length, height at hip cross, chest circumference, chest width, hip width, and cannon circumference.
The meat and cashmere dual-purpose SBWC goats were adult females (2–3 years old), and reared in the Shaanbei white cashmere goat breeding farm in Yulin city of Shaanxi Province [1,34]. Ear samples were randomly collected in these populations between 2016 and 2018 [1,34]. These goats were unrelated and healthy. For this breed, we measured seven values of body size, including the body height, body length, height at hip cross, chest circumference, chest depth, chest width, and cannon circumference [34]. In addition, the adult female HNBGs (2–3 years old) were raised at a local farm in Zanzhou country of Hainan Province [35]. This breed is a famous meat-type of goat breed; therefore, the data of the body measurement traits were recorded, including the body height, body length, body weight, chest circumference, chest depth, chest width, hip width, and cannon circumference [35].

2.2. DNA Extraction and Indel Genotyping

The genomic DNA samples were extracted from the ear tissue samples using the high-salt extraction method [36]. A Nanodrop 1000 Spectrometer (Thermo Scientific, Waltham, MA, USA) was used to detect the DNA purity and quality. Then, the extracted DNA were diluted to 10 ng/μL and stored at −20 °C.
According to the previous report [1], a pair of primers (F: 5′-GCTGGAAGCAGTTCGTCA-3′; R: 5′-GGGTTGATAGCCTTGTATGTT-3′) was synthesized to detect the 11-bp indel within the CSN1S1 gene in all of the goats. The PCR amplification processes used were the same as in the previous study, including the primers, PCR reaction mixture, and amplification condition. If the PCR products migrated at the same position on the agarose gel when they were genotyped by 3.0% agarose gel electrophoresis, as previously reported, the genotype of each goat could be visualized. For example, one band appearing at 170 bp meant that the individual was of a homozygote insertion type (insertion/insertion, II), one band appearing at 159 bp meant that the individual was of a homozygote deletion type (deletion/deletion, DD), and two bands appearing in a lane meant that the individual was of a heterozygote type (insertion/deletion, ID). Meanwhile, the PCR products were confirmed by sequencing and compared to the reference genome (NC_030813.1; Tsingke company, Xi’an, Shaanxi, China).

2.3. Statistical Analysis

The genotypes, allele frequencies, and a test for the Hardy–Weinberg equilibrium (HWE) were calculated using the SHEsis platform (http://analysis.bio-x.cn/myAnalysis.php) [37]. The population indices (polymorphism information content, PIC; homozygosity, Ho; heterozygosity, He; effective allele numbers, Ne) were calculated following Nei’s methods [38].
The general linear model (GLM) of the analysis of variance (ANOVA) was used to analyze the relationships between the different genotypes and economically important traits using the SPSS 23.0 statistical software package (IBM, New York, NY, USA). The statistical model for the milk performance traits used was as follows: Yijk = µ + Li + Gj + (LG)ij + εijk, where Yijk is the observation of the milk performance indices, µ is the overall mean of each index, Li is the effect due to ith lactation, Gj is the fixed effect of jth genotype, (LG)ij is interaction between the ith lactation and the jth genotype, and εijk is the random error [33]. Meanwhile, for the body measurement traits, the GLM was calculated as follows: Yil = μ + Ai + Gl + εil, where Yil represents the observations of the body measurement traits, μ is the population mean, Ai is the fixed effect of age, Gl is the fixed effect of genotype, and εil is the random error [39,40]. The model excluded the effects of farm, sex, and season of birth, which had no significant effects on the variation of traits in this population [40]. The different breeds were raised in their respective farms and analyzed separately. All of the data were presented as the mean ± standard error (S.E.), and p < 0.05 was considered to be significant.

3. Results

3.1. Genetic Parameter of 11-bp Indel of CSN1S1 in Three Chinese Goat Breeds

Based on the gel electrophoresis pattern and sequence chromatograms (Figure 1), an 11-bp indel of the CSN1S1 gene was detected in all of the tested goats. The polymorphism was present in the following three variants: the homozygote insertion type (II, 170 bp), the homozygote deletion type (DD, 159 bp), and the heterozygote type (ID, 170 bp and 159 bp; Figure 1A). After sequence alignment with the updated reference genome (NC_030813.1), it was found that the polymorphism was an insertion mutation, and the position changed from 6654–6664 to 8045–8055 (Figure 1B). Hence, we redefined this polymorphism as NC_030813.1: g.8045–8055insTTTCCGTAATG, compared to the previous name (NC_030813.1: g.6654–6664delTTTCCGTAATG) given by Wang et al. [1]. Moreover, in this study, the heterozygote type “ID” actually showed three bands on the agarose gel (Figure 1A), as the top band was a heteroduplex [41,42].
Subsequently, the genotype and allele frequencies of this polymorphism in the three goat breeds were calculated (Table 1). The minor allele frequencies of this indel in GZDG, SBWC, and HNBG were 0.345, 0.478, and 0.448, respectively. The frequency of the ID genotype was higher than the other genotypes in SBWC and HNBG; however, in the GZDG population, the highest genotype frequency (0.481) was the homozygote deletion type “DD”. Based on the frequency numbers of the different genotypes, the population indices (Ho, He, Ne, and PIC), indicating the degrees of genetic variation, were calculated (Table 1). The Ho, He, and Ne values of this polymorphism were high in all three populations; meanwhile, the PIC value showed that all of the breeds had moderate genetic diversity (0.25 < PIC < 0.5) at this indel locus of the CSN1S1 gene. In addition, the χ2 test results demonstrated that no breed was in accordance with HWE (p < 0.05; Table 1).

3.2. Association between the 11-bp Indel and Milk Performance Traits in the Guanzhong Dairy Goat

To assess the effect of the 11-bp indel on the milk performance traits, the milk samples were collected and eight milk physicochemical indices were measured in the Guanzhong dairy goat (GZDG) [33]. The analysis result demonstrated that this indel was significantly associated with five of the indices in dairy goat, including milk fat content (%), TS content (%), SNF content (%), FPD (°C), and acidity (pH; Table 2; p < 0.05). Goats with II or ID genotypes had higher values of milk physicochemical indices compared with individuals with the DD genotype (p < 0.05).

3.3. Association Analysis of the 11-bp Indel with Body Measurement Traits in Three Goat Breeds

The body measurement traits can directly reflect the skeleton structure, body size, and growth development, which are related to the physiological function and the production performance of goats [35]. Here, the various body sizes were collected to analyze the relationships between this polymorphism and the body measurement traits in different breeds, including the meat-use breeds (HNBG and SBWC) and the dairy-use breed (GZDG). The association analysis results revealed that the 11-bp indel was significantly associated with different body measurement traits in different goat breeds (Table 3). The indel was significantly related to the BH of all three of the breeds (p < 0.05). Two body measurement traits—body height and body length—were related in the GZDG population (p < 0.05); six kinds of body measurement traits—body height, body length, height at hip cross, chest circumference, chest width, and cannon circumference—were associated with this polymorphism in SBWC goats (p < 0.05); and, in the HNBG population, three different body measurement traits—body height, chest depth, and body weight—were associated with the polymorphism (p < 0.05). An interesting tendency that was observed was that female goats with DD or ID genotypes had greater body measurement traits than those of other genotypes in the three goat breeds (p < 0.05).

4. Discussion

The CSN1S1 gene encodes α-s1 casein, one of the proteins constituting milk that can improve calcium uptake in mammals, which affects milk performance and growth development [14,15,23,26]. Many studies have shown that CSN1S1 polymorphisms can affect the milk performance in different goat breeds [19,20,21,22]. Our team has previously identified an 11-bp indel within the CSN1S1 gene, which strongly affected litter size in a large sample of SBWC goats (n = 3047) [1]. However, whether this polymorphism affects the milk performance and body measurement traits of goats has not yet been investigated. In this paper, we detected this indel in three local goat breeds, and showed the following three genotypes: II (homozygote insertion type), ID (heterozygote type), and DD (homozygote deletion type). As the reference caprine genome sequence of CSN1S1 (NC_030813.1) was updated, the polymorphism was redefined as NC_030813.1: g.8045–8055insTTTCCGTAATG. In three goat breeds, the average frequency of the “D” allele was 56.2%, which was similar to that reported in a previous report on SBWC goats (54.0%) by Wang et al. [1]. However, in SBWC goats, the frequency of the “D” allele was lower than that of the previous study, which could be due to the sample sizes in this experiment being smaller. Furthermore, no breed was in accordance with the HWE, suggesting that the populations of the three breeds were not in dynamic equilibrium, due to migration, genetic drift, and artificial selection [43]. Overall, this polymorphic mutation had a high frequency and was present in all three goat breeds from different regions, indicating that this polymorphic mutation is genetically stable in goats and shows potential as a candidate DNA marker in molecular breeding progress.
Subsequently, the association analysis uncovered that the strong influences of the 11-bp indel of CSN1S1 on the milk performance traits in dairy goats (Table 2). According to the results, dairy goats of the II or ID genotypes demonstrated better milk performance traits, particularly in terms of milk protein content (%), TS content (%), and SNF content (%; p < 0.05). Higher values of these traits indicate that a goat’s milk contains more dry matter, especially proteins, because there is no difference in milk fat content (%) and lactose content (%) among individuals of different genotypes (p > 0.05). The reason for this difference may be due to the fact that the CSN1S1 gene encodes α-s1 casein, which is the richest protein in milk [14]. Therefore, the polymorphism may affect the expression of the protein, resulting in changes in milk protein content and other components of dairy goats’ milk. The increase in the milk protein improves the nutritional value of milk and is also beneficial in the cheese-making process [44,45]. However, this may also have some negative effects. Some literature indicates that the variants led to increased milk components, but reduced the milk yield [46,47,48]. Regrettably, in this study, because of the lack of data on the milk yield of dairy goats, the effect of this indel on this aspect cannot be studied. In the future, more comprehensive data about milk performance will need to be collected in order to explore the impact of this 11-bp indel.
In addition, as caseins could promote bone formation and body development by increasing intestinal calcium absorption and bioavailability [14,26,27,28,49], the effect of this 11-bp indel of the CSN1S1 gene on the body measurement traits of goats is worth studying. The association analysis results demonstrated that this polymorphism significantly affected the body measurement traits of goats and affected the different body measurement sizes in different goat breeds (Table 3; p < 0.05). These differences might be due to the diversity of goat breeds, and individuals with DD or ID genotypes had better body measurement traits, compared with the II type, in all three breeds (p < 0.05). Therefore, for selecting the body measurement trait of goats, the DD genotype of the 11-bp indel within the CSN1S1 gene is the dominant type.
Intriguingly, the dominant genotype of milk performance (Table 2) and litter size [1] in goats was the II genotype for the 11-bp indel polymorphism within the CSN1S1 gene, but the dominant genotype of the body measurement trait is the DD genotype and individuals with the II genotype had lower body measurement sizes (Table 3). In general, female animals, after parturition and during lactation, need to provide more nutrients to their offspring, and so their body conditions decline [50,51,52]. For goat breeds dedicated to milk production, their body sizes are not as big as that of meat-producing goats, because it is less important to increase body size instead of milk production. This may explain why the dominant genotypes are the opposite in the two traits of goats. Therefore, according to the existing analysis results, the genotype selection for breeding can be carried out based on the use of goats, as follows: if the milk performance needs to be improved, for example, in dairy goat breeds, goats with the II genotype of the 11-bp indel in the CSN1S1 gene should be selected; on the other hand, for improving body measurement traits in goats, such as in meat-use or meat–cashmere dual purpose goat breeds, the DD genotype should be chosen. However, as the number of samples used to study milk performance traits in this experiment was not large and only one breed type was used (only GZDG), no decisive conclusion can be drawn. It is necessary to increase the sample size and collect more production performance indices with different goat breeds for further conclusions to be drawn.
In addition, considering that this indel is an intron variation and often has no direct effect on gene expression, it might be linked with other polymorphisms in the exon or the transcriptional regulation region to exert indirect effects [53,54]. In fact, in our previous study on litter size, we explored whether this 11-bp indel was linked to the polymorphisms of the BMPR1B gene in goats; however, the Q279R mutation of BMPR1B was not found in SBWC population [1]. Therefore, it is possible that this polymorphism affects litter size by it linking to other potential mutations of the BMPR1B gene or other genes in goats [1]. Unfortunately, at present, most studies have only focused on the relationships between this gene and the milk performance traits, and only one study has reported on the associations between genetic variants of the CNS1S1 gene and different economically important traits in goats [55]. The author found that some SNPs within exons 9–11 significantly affected the milk yield, litter size, and body measurement traits in the Xinong Sannen dairy goats [55]. Therefore, we suspect that this indel might link with these SNPs to control the production traits of goats. In addition, the intronic polymorphisms might also change the transcription and translation processes of the gene by binding trans-acting factors [56,57], altering the gene splicing pattern [58,59], or mediating the biogenesis of the microRNA [60]. The possible mechanisms of this 11-bp indel polymorphism (i.e., those affecting the economically important traits) still need to be explored.

5. Conclusions

An 11-bp indel in the CSN1S1 gene significantly affected the milk performance and body measurement traits in goats, and the dominant genotypes were opposite in the two traits. Therefore, the genotype selection of breeding in the MAS program can be carried out based on the use of goats, as follows: if the milk performance needs to be improved, for example, in dairy goat breeds, goats with the II genotype of this polymorphism should be selected; on the other hand, for improving the body measurement traits into goats, such as in meat-use breeds, the DD genotype should be chosen.

Author Contributions

conceptualization, Y.Z., C.P., and X.S.; formal analysis, Y.Z. and K.W.; funding acquisition, C.P. and X.S.; investigation, K.W., J.L., H.Z., and L.Q.; methodology, H.C. and X.L.; project administration, X.L., C.P., and X.S; resources, X.L., J.L., H.Z., L.Q., and X.S.; software, Y.Z. and K.W.; supervision, C.P. and X.S.; validation, C.P. and X.S.; visualization, C.P. and X.S.; writing (original draft), Y.Z. and K.W.; writing (review and editing), Y.Z., C.P., and X.S.

Funding

This research was funded by the National Natural Science Foundation of China (31660650; 31172184), and the Provincial Key Projects of Shaanxi (2014KTDZ02-01).

Acknowledgments

We greatly thanked the staff of the Guanzhong dairy goat breeding farm (Shaanxi, China), the Shaanbei Cashmere goat breeding farm (Shaanxi, China), and the Hainan black goat breeding farm (Hainan, China) for samples and data collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, K.; Yan, H.; Xu, H.; Yang, Q.; Zhang, S.; Pan, C.; Chen, H.; Zhu, H.; Liu, J.; Qu, L.; et al. A novel indel within goat casein alpha S1 gene is significantly associated with litter size. Gene 2018, 671, 161–169. [Google Scholar] [CrossRef] [PubMed]
  2. Zhang, B.; Chang, L.; Lan, X.; Asif, N.; Guan, F.; Fu, D.; Li, B.; Yan, C.; Zhang, H.; Zhang, X.; et al. Genome-wide definition of selective sweeps reveals molecular evidence of trait-driven domestication among elite goat (Capra species) breeds for the production of dairy, cashmere, and meat. Gigascience 2018, 7, giy105. [Google Scholar] [CrossRef] [PubMed]
  3. Peacock, C. Goats—A pathway out of poverty. Small. Rumin. Res. 2005, 60, 179–186. [Google Scholar] [CrossRef]
  4. Dubeuf, J.P.; Morand-Fehr, P.; Rubino, R. Situation, changes and future of goat industry around the world. Small. Rumin. Res. 2004, 51, 165–173. [Google Scholar] [CrossRef]
  5. FAO Commission on Genetic Resources for Food and Agriculture Assessments. The Second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture; FAO: Rome, Italy, 2015; Available online: http://www.fao.org/publications/sowangr/en/ (accessed on 18 September 2019).
  6. Kang, Z.; Jiang, E.; Wang, K.; Pan, C.; Chen, H.; Yan, H.; Zhu, H.; Liu, J.; Qu, L.; Lan, X. Goat membrane associated ring-CH-type finger 1 (MARCH1) mRNA expression and association with litter size. Theriogenology 2019, 128, 8–16. [Google Scholar] [CrossRef]
  7. Yang, Q.; Zhang, S.; Li, J.; Wang, X.; Peng, K.; Lan, X.; Pan, C. Development of a touch-down multiplex PCR method for simultaneously rapidly detecting three novel insertion/deletions (indels) within one gene: An example for goat GHR gene. Anim. Biotechnol. 2019, 30, 366–371. [Google Scholar] [CrossRef]
  8. Jonas, E.; de Koning, D.J. Genomic selection needs to be carefully assessed to meet specific requirements in livestock breeding programs. Front. Genet. 2015, 6, 49. [Google Scholar] [CrossRef] [Green Version]
  9. Zhu, H.; Zhang, Y.; Bai, Y.; Yang, H.; Yan, H.; Liu, J.; Shi, L.; Song, X.; Li, L.; Dong, S.; et al. Relationship between SNPs of POU1F1 gene and litter size and growth traits in Shaanbei white cashmere goats. Animals 2019, 9, 114. [Google Scholar] [CrossRef] [Green Version]
  10. Chen, M.; Wang, J.; Liu, N.; Cui, W.; Dong, W.; Xing, B.; Pan, C. Pig SOX9: Expression profiles of Sertoli cell (SCs) and a functional 18 bp indel affecting testis weight. Theriogenology 2019, 138, 94–101. [Google Scholar] [CrossRef]
  11. Kang, Z.; Zhang, S.; He, L.; Zhu, H.; Wang, Z.; Yan, H.; Huang, Y.; Dang, R.; Lei, C.; Chen, H.; et al. A 14-bp functional deletion within the CMTM2 gene is significantly associated with litter size in goat. Theriogenology 2019, 139, 49–57. [Google Scholar] [CrossRef]
  12. Wang, X.; Yang, Q.; Zhang, S.; Zhang, X.; Pan, C.; Chen, H.; Zhu, H.; Lan, X. Genetic Effects of single nucleotide polymorphisms in the goat GDF9 gene on prolificacy: True or false positive? Animals 2019, 9, 886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ren, T.; Li, W.; Liu, D.; Liang, K.; Wang, X.; Li, H.; Jiang, R.; Tian, Y.; Kang, X.; Li, Z. Two insertion/deletion variants in the promoter region of the QPCTL gene are significantly associated with body weight and carcass traits in chickens. Anim. Genet. 2019, 50, 279–282. [Google Scholar] [CrossRef] [PubMed]
  14. Rezaei, R.; Wu, Z.; Hou, Y.; Bazer, F.W.; Wu, G. Amino acids and mammary gland development: Nutritional implications for milk production and neonatal growth. J. Anim. Sci. Biotechnol. 2016, 7, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Juvarajah, T.; Wan-Ibrahim, W.I.; Ashrafzadeh, A.; Othman, S.; Hashim, O.H.; Fung, S.Y.; Abdul-Rahman, P.S. Human milk fat globule membrane contains hundreds of abundantly expressed and nutritionally beneficial proteins that are generally lacking in caprine milk. Breastfeed. Med. 2018, 13, 631–637. [Google Scholar] [CrossRef] [PubMed]
  16. Grosclaude, F.; Martin, P. Casein polymorphism in the goat. In Milk Protein Polymorphism; International Dairy Federation, Ed.; Special Issue, 9702; International Dairy Federation: Brussels, Belgium, 1997; pp. 241–253. [Google Scholar]
  17. Martin, P.; Szymanowska, M.; Zwierzchowski, L.; Leroux, C. The impact of genetic polymorphisms on the protein composition of ruminant milks. Reprod. Nutr. Dev. 2002, 42, 433–459. [Google Scholar] [CrossRef]
  18. Guan, D.; Mármol-Sánchez, E.; Cardoso, T.F.; Such, X.; Landi, V.; Tawari, NR.; Amills, M. Genomic analysis of the origins of extant casein variation in goats. J. Dairy Sci. 2019, 102, 5230–5241. [Google Scholar] [CrossRef]
  19. Yue, X.P.; Zhang, X.M.; Wang, W.; Ma, R.N.; Deng, C.J.; Lan, X.Y.; Chen, H.; Li, F.; Xu, X.R.; Ma, Y.; et al. The CSN1S1 N and F alleles identified by PCR-SSCP and their associations with milk yield and composition in Chinese dairy goats. Mol. Biol. Rep. 2011, 38, 2821–2825. [Google Scholar] [CrossRef]
  20. Vacca, G.M.; Dettori, M.L.; Piras, G.; Manca, F.; Paschino, P.; Pazzola, M. Goat casein genotypes are associated with milk production traits in the Sarda breed. Anim. Genet. 2014, 45, 723–731. [Google Scholar] [CrossRef]
  21. Mestawet, T.A.; Girma, A.; Ådnøy, T.; Devold, T.G.; Vegarud, G.E. Newly identified mutations at the CSN1S1 gene in Ethiopian goats affect casein content and coagulation properties of their milk. J. Dairy Sci. 2013, 96, 4857–4869. [Google Scholar] [CrossRef]
  22. Carillier-Jacquin, C.; Larroque, H.; Robert-Granié, C. Including αs1 casein gene information in genomic evaluations of French dairy goats. Genet. Sel. Evol. 2016, 48, 54. [Google Scholar] [CrossRef] [Green Version]
  23. 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] [PubMed] [Green Version]
  24. Noce, A.; Pazzola, M.; Dettori, M.L.; Amills, M.; Castelló, A.; Cecchinato, A.; Bittante, G.; Vacca, G.M. Variations at regulatory 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] [PubMed]
  25. Bonfatti, V.; Giantin, M.; Gervaso, M.; Coletta, A.; Dacasto, M.; Carnier, P. Effect of CSN1S1-CSN3 (αS1-κ-casein) composite genotype on milk production traits and milk coagulation properties in Mediterranean water buffalo. J. Dairy Sci. 2012, 95, 3435–3443. [Google Scholar] [CrossRef] [PubMed]
  26. Masarwi, M.; Gabet, Y.; Dolkart, O.; Brosh, T.; Shamir, R.; Phillip, M.; Gat-Yablonski, G. Skeletal effect of casein and whey protein intake during catch-up growth in young male Sprague-Dawley rats. Br. J. Nutr. 2016, 116, 59–69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Liu, G.; Miao, J.; Sun, S.; Luo, Z.; Xia, Z.; Guo, B.; Liu, F.; Chen, Y.; Ying, D.; Miao, B.; et al. The effect of dietary magnesium and caseinphosphopeptides on bone metabolism in rats. Food. Funct. 2017, 8, 4487–4495. [Google Scholar] [CrossRef] [PubMed]
  28. Liu, G.; Sun, S.; Guo, B.; Miao, B.; Luo, Z.; Xia, Z.; Ying, D.; Liu, F.; Guo, B.; Tang, J.; et al. Bioactive peptide isolated from casein phosphopeptides promotes calcium uptake in vitro and in vivo. Food. Funct. 2018, 9, 2251–2260. [Google Scholar] [CrossRef] [PubMed]
  29. Mulsant, P.; Lecerf, F.; Fabre, S.; Schibler, L.; Monget, P.; Lanneluc, I.; Pisselet, C.; Riquet, J.; Monniaux, D.; Callebaut, I.; et al. Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Mérino ewes. Proc. Natl. Acad. Sci. USA. 2001, 98, 5104–5109. [Google Scholar] [CrossRef] [Green Version]
  30. Tang, J.; Hu, W.; Di, R.; Liu, Q.; Wang, X.; Zhang, X.; Zhang, J.; Chu, M. Expression analysis of the prolific candidate genes, BMPR1B, BMP15, and GDF9 in Small Tail Han ewes with three fecundity (FecB gene) genotypes. Animals 2018, 8, 166. [Google Scholar] [CrossRef] [Green Version]
  31. Montgomery, G.W.; Lord, E.A.; Penty, J.M.; Dodds, K.G.; Broad, T.E.; Cambridge, L.; Sunden, S.L.; Stone, R.T.; Crawford, A.M. The Booroola fecundity (FecB) gene maps to sheep chromosome 6. Genomics 1994, 22, 148–153. [Google Scholar] [CrossRef]
  32. Lan, X.Y.; Chen, H.; Pan, C.Y.; Li, R.B.; Li, X.C.; Fang, X.T. Relationship between polymorphisms of CSN3, CSN1S2 and β-lg genes and litter sizes of Xinong Saanen dairy goat. Sci. Agric. Sin. 2005, 11, 2333–2338. [Google Scholar]
  33. Lan, X.Y.; Zhao, H.Y.; Li, Z.J.; Zhou, R.; Pan, C.Y.; Lei, C.Z.; Chen, H. Exploring the novel genetic variant of PITX1 gene and its effect on milk performance in dairy goats. J. Integr. Agric. 2013, 12, 118–126. [Google Scholar] [CrossRef]
  34. Wang, X.Y.; Yang, Q.; Wang, K.; Zhang, S.H.; Pan, C.Y.; Chen, H.; Qu, L.; Yan, H.L.; Lan, X.Y. A novel 12-bp indel polymorphism within the GDF9 gene is significantly associated with litter size and growth traits in goats. Anim. Genet. 2017, 48, 735–736. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, S.; Jiang, E.; Wang, K.; Zhang, Y.; Yan, H.; Qu, L.; Chen, H.; Lan, X.; Pan, C. Two insertion/deletion variants within SPAG17 gene are associated with goat body measurement traits. Animals 2019, 9, 379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Aljanabi, S.M.; Martinez, I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic. Acids. Res. 1997, 25, 4692–4693. [Google Scholar] [CrossRef]
  37. Li, Z.Q.; Zhang, Z.; He, Z.D.; Tang, W.; Li, T.; Zeng, Z.; He, L.; Shi, YY. A partition-ligation-combination-subdivision EM algorithm for haplotype inference with multiallelic markers: Update of the SHEsis (http://analysis.bio-x.cn). Cell. Res. 2009, 19, 519–523. [Google Scholar] [CrossRef]
  38. Nei, M. Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. USA. 1973, 70, 3321–3323. [Google Scholar] [CrossRef] [Green Version]
  39. Wang, X.; Yang, Q.; Wang, K.; Yan, H.; Pan, C.; Chen, H.; Liu, J.; Zhu, H.; Qu, L.; Lan, X. Two strongly linked single nucleotide polymorphisms (Q320P and V397I) in GDF9 gene are associated with litter size in cashmere goats. Theriogenology 2019, 125, 115–121. [Google Scholar] [CrossRef]
  40. Lan, X.Y.; Pan, C.Y.; Chen, H.; Zhang, C.L.; Li, J.Y.; Zhao, M.; Lei, C.Z.; Zhang, A.L.; Zhang, L. An AluI PCR-RFLP detecting a silent allele at the goat POU1F1 locus and its association with production traits. Small. Rumin. Res. 2007, 73, 8–12. [Google Scholar] [CrossRef]
  41. Ren, F.; Yu, S.; Chen, R.; Lv, X.; Pan, C. Identification of a novel 12-bp insertion/deletion (indel) of iPS-related Oct4 gene and its association with reproductive traits in male piglets. Anim. Reprod. Sci. 2017, 178, 55–60. [Google Scholar] [CrossRef]
  42. Nagamine, C.M.; Chan, K.; Lau, Y.F. A PCR antifact generation of heteroduplexes. Am. J. Hum. Genet. 1989, 45, 337–339. [Google Scholar]
  43. Zhang, X.Y.; Wu, X.F.; Jia, W.C.; Pan, C.Y.; Li, X.C.; Lei, C.Z.; Chen, H.; Lan, X.Y. Novel nucleotide variations, haplotypes structure and associations with growth related traits of goat AT motif-binding factor (ATBF1) gene. Asian Australas. J. Anim. 2015, 28, 1394–1406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Młynek, K.; Oler, A.; Zielińska, K.; Tkaczuk, J.; Zawadzka, W. The effect of selected components of milk and ripening time on the development of the hardness and melting properties of cheese. Acta. Sci. Pol. Technol. Aliment. 2018, 17, 133–140. [Google Scholar] [PubMed]
  45. Pazzola, M.; Stocco, G.; Dettori, M.L.; Bittante, G.; Vacca, G.M. Effect of goat milk composition on cheesemaking traits and daily cheese production. J. Dairy. Sci. 2019, 102, 3947–3955. [Google Scholar] [CrossRef] [PubMed]
  46. Li, Y.; Zhou, H.; Cheng, L.; Edwards, G.R.; Hickford, J.G. Effect of DGAT1 variant (K232A) on milk traits and milk fat composition in outdoor pasture-grazed dairy cattle. N. Z. J. Agric. Res. 2019, 1–13. [Google Scholar] [CrossRef]
  47. Xiang, R.; MacLeod, I.M.; Bolormaa, S.; Goddard, M.E. Genome-wide comparative analyses of correlated and uncorrelated phenotypes identify major pleiotropic variants in dairy cattle. Sci. Rep. 2017, 7, 9248. [Google Scholar] [CrossRef]
  48. Vacca, G.M.; Dettori, M.L.; Balia, F.; Luridiana, S.; Mura, M.C.; Carcangiu, V.; Pazzola, M. Sequence polymorphisms at the growth hormone GH1/GH2-N and GH2-Z gene copies and their relationship with dairy traits in domestic sheep (Ovis aries). Mol. Biol. Rep. 2013, 40, 5285–5294. [Google Scholar] [CrossRef]
  49. Bennett, T.; Desmond, A.; Harrington, M.; McDonagh, D.; Fitzgerald, R.J.; Flynn, A.; Cashman, K.D. The effect of high intakes of casein and casein phosphopeptide on calcium absorption in the rat. Br. J. Nutr. 2000, 83, 673–680. [Google Scholar] [CrossRef] [Green Version]
  50. Kim, I.H.; Suh, G.H. Effect of the amount of body condition loss from the dry to near calving periods on the subsequent body condition change, occurrence of postpartum diseases, metabolic parameters and reproductive performance in Holstein dairy cows. Theriogenology 2003, 60, 1445–1456. [Google Scholar] [CrossRef]
  51. Roche, J.R.; Meier, S.; Heiser, A.; Mitchell, M.D.; Walker, C.G.; Crookenden, M.A.; Riboni, M.V.; Loor, J.J.; Kay, J.K. Effects of precalving body condition score and prepartum feeding level on production, reproduction, and health parameters in pasture-based transition dairy cows. J. Dairy Sci. 2015, 98, 7164–7182. [Google Scholar] [CrossRef] [Green Version]
  52. Smith, G.L.; Friggens, N.C.; Ashworth, C.J.; Chagunda, M.G.G. Association between body energy content in the dry period and post-calving production disease status in dairy cattle. Animal 2017, 11, 1590–1598. [Google Scholar] [CrossRef] [Green Version]
  53. Abecassis, I.; Sedgewick, A.J.; Romkes, M.; Buch, S.; Nukui, T.; Kapetanaki, M.G.; Vogt, A.; Kirkwood, J.M.; Benos, P.V.; Tawbi, H. PARP1 rs1805407 increases sensitivity to PARP1 inhibitors in cancer cells suggesting an improved therapeutic strategy. Sci. Rep. 2019, 9, 3309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Ganapathipillai, S.; Laval, G.; Hoffmann, I.S.; Castejon, A.M.; Nicod, J.; Dick, B.; Frey, F.J.; Frey, B.M.; Cubeddu, L.X.; Ferrari, P. CYP11B2-CYP11B1 haplotypes associated with decreased 11 beta-hydroxylase activity. J. Clin. Endocrinol. Metab. 2005, 90, 1220–1225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Yu, J. Studies on Genetic Diversity of CSN1S1, BMPR-1B, BMP15 and IGFBP3 genes in two dairy goat breeds. Master’s Thesis, Northwest A&F University, Xianyang, China, 2006. [Google Scholar]
  56. Farias, F.H.G.; Dahlqvist, J.; Kozyrev, S.V.; Leonard, D.; Wilbe, M.; Abramov, S.N.; Alexsson, A.; Pielberg, G.R.; Hansson-Hamlin, H.; Andersson, G.; et al. A rare regulatory variant in the MEF2D gene affects gene regulation and splicing and is associated with a SLE sub-phenotype in Swedish cohorts. Eur. J. Hum. Genet. 2019, 27, 432–441. [Google Scholar] [CrossRef] [PubMed]
  57. Markljung, E.; Jiang, L.; Jaffe, J.D.; Mikkelsen, T.S.; Wallerman, O.; Larhammar, M.; Zhang, X.; Wang, L.; Saenz-Vash, V.; Gnirke, A.; et al. ZBED6, a novel transcription factor derived from a domesticated DNA transposon regulates IGF2 expression and muscle growth. PLoS. Biol. 2009, 7, e1000256. [Google Scholar] [CrossRef] [PubMed]
  58. Mishra, S.K.; Thakran, P. Intron specificity in pre-mRNA splicing. Curr. Genet. 2018, 64, 777–784. [Google Scholar] [CrossRef]
  59. Hashiramoto, A.; Konishi, Y.; Murayama, K.; Kawasaki, H.; Yoshida, K.; Tsumiyama, K.; Tanaka, K.; Mizuhara, M.; Shiotsuki, T.; Kitamura, H.; et al. A variant of death-receptor 3 associated with rheumatoid arthritis interferes with apoptosis-induction of T cell. J. Biol. Chem. 2018, 293, 1933–1943. [Google Scholar] [CrossRef] [Green Version]
  60. Hecker, M.; Boxberger, N.; Illner, N.; Fitzner, B.; Schröder, I.; Winkelmann, A.; Dudesek, A.; Meister, S.; Koczan, D.; Lorenz, P.; et al. A genetic variant associated with multiple sclerosis inversely affects the expression of CD58 and microRNA-548ac from the same gene. PLoS. Genet. 2019, 15, e1007961. [Google Scholar] [CrossRef]
Animals 09 01114 g001 550
Figure 1. Visualization and sequencing of the 11-bp indel under study in the goat casein alpha s1 (CSN1S1) gene. (A) The electrophoresis pattern of the 11-bp indel within the goat CSN1S1 gene. Lanes 1, 4, and 5 show the homozygote deletion type (DD); lane 3 shows the homozygote insertion type (II); lanes 2 and 6 show the heterozygote type (ID); lane 7 is a negative control; and lane 8 is a DNA ladder marker. (B) The sequence chromatograms of the 11-bp indel in goat CSN1S1. The upper and lower panels show the sequence chromatograms of the II and DD genotypes, respectively.
Figure 1. Visualization and sequencing of the 11-bp indel under study in the goat casein alpha s1 (CSN1S1) gene. (A) The electrophoresis pattern of the 11-bp indel within the goat CSN1S1 gene. Lanes 1, 4, and 5 show the homozygote deletion type (DD); lane 3 shows the homozygote insertion type (II); lanes 2 and 6 show the heterozygote type (ID); lane 7 is a negative control; and lane 8 is a DNA ladder marker. (B) The sequence chromatograms of the 11-bp indel in goat CSN1S1. The upper and lower panels show the sequence chromatograms of the II and DD genotypes, respectively.
Animals 09 01114 g001
Table
Table 1. The genotype and allele frequencies and population indices of the 11-bp indel of CSN1SI gene in three different goat breeds.
Table 1. The genotype and allele frequencies and population indices of the 11-bp indel of CSN1SI gene in three different goat breeds.
BreedsTotal SizesGenotype Frequencies
(N)		Allele FrequenciesPopulation Parametersχ2
(p Value)
IIIDDDIDHoHeNePIC
GZDG2350.170
(40)		0.429
(82)		0.481
(113)		0.3450.6550.5480.4521.8240.35012.173
(p = 4.85 × 10−4)
SBWC10920.295
(322)		0.454
(496)		0.251
(274)		0.5220.4780.5010.4991.9960.3758.809
(p = 0.003)
HNBG2780.255
(71)		0.385
(107)		0.360
(100)		0.4480.5520.5050.4951.9780.37213.670
(p = 2.18 × 10−4)
Note: GZDG—Guanzhong dairy goat; SBWC—Shaanbei white cashmere goat; HNBG—Hainan black goat; Ho—homozygosity; He—heterozygosity; Ne—effective allele numbers; PIC—polymorphism information content.
Table
Table 2. Association between the indel in the CSN1S1 gene and milk performance indices in the Guanzhong dairy goats.
Table 2. Association between the indel in the CSN1S1 gene and milk performance indices in the Guanzhong dairy goats.
Milk Performance IndicesGenotypes (N)p Value
II (11)ID (35)DD (46)
Fat content (%)4.61 ± 0.404.36 ± 0.203.86 ± 0.150.103
Protein content (%)3.93 a ± 0.213.73 ab ± 0.123.38 b ± 0.070.006
TS content (%)13.73 a ± 0.5713.38 ab ± 0.3412.44 b ± 0.210.018
SNF content (%)8.99 a ± 0.198.91 a ± 0.138.46 b ± 0.080.003
Lactose content (%)3.88 ± 0.094.02 ± 0.083.95 ± 0.040.481
Density1028.80 ± 0.451029.25 ± 0.461028.27 ± 0.270.144
FPD (°C)0.50a ± 0.010.50 a ± 0.010.48 b ± 0.010.043
Acidity (pH)7.49a ± 0.516.79 ab ± 0.336.22 b ± 0.190.046
Note: TS—total solids; SNF—solids-not-fat; FPD—freezing point depression. Values with different letters (a and b) within the same row differ significantly at p < 0.05.
Table
Table 3. Relationships between the 11-bp indel of the CSN1S1 gene and body measurement traits in the three goat breeds.
Table 3. Relationships between the 11-bp indel of the CSN1S1 gene and body measurement traits in the three goat breeds.
BreedsBody Measurement TraitsGenotypes (N)p Value
II (11)ID (35)DD (46)
GZDGBody height (cm)67.98 c ± 3.0969.05 b ± 2.8971.26 a ± 0.900.009
Body length (cm)75.54 b ± 0.9877.24 a ± 1.1877.60 a ± 0.680.043
Height at hip cross (cm)68.98 ± 2.9169.33 ± 1.1170.81 ± 0.950.571
Chest circumference (cm)89.18 ± 0.9488.28 ± 0.9888.92 ± 0.820.841
Chest width (cm)28.70 ± 0.8527.81 ± 0.5529.17 ± 0.490.180
Hip width (cm) 21.94 ± 1.4123.54 ± 0.5922.26 ± 0.670.293
Cannon circumference (cm)11.34 ± 0.4010.89 ± 0.2110.30 ± 0.270.099
II (247)ID (380)DD (195)
SBWCBody height (cm)56.62 c ± 0.2958.11 b ± 0.1960.12 a ± 0.331.61 × 10−16
Body length (cm)66.42 b ± 0.3167.95 a ±0.2368.62 a ± 0.331.23 × 10−6
Height at hip cross (cm)59.41 c ± 0.2960.66 b ± 0.2161.97 a ± 0.304.22 × 10−9
Chest circumference (cm)84.97 b ± 0.4986.78 a ± 0.3586.93 a ± 0.500.003
Chest width (cm)19.59 b ± 0.1820.36 a ±0.1420.48 a ±0.18 3.90 × 10−4
Chest depth (cm)30.61 ± 01.1630.15 ± 0.1930.46 ± 0.180.851
Cannon circumference (cm)8.13 b ± 0.058.32 a ± 0.038.26 a ± 0.050.003
II (26)ID (39)DD (36)
HNBGBody height (cm)51.50 b ± 0.8953.63 a ± 0.6453.98 a ± 0.650.049
Body length (cm)54.71 ± 1.1456.50 ± 0.7557.55 ± 0.640.069
Body weight (kg)26.79 b ± 1.4829.35 ab ± 1.1731.34 a ± 1.000.043
Chest circumference (cm)71.04 ± 1.3973.00 ± 1.1174.40 ± 0.900.137
Chest width (cm)15.10 ± 0.3315.61 ± 0.3315.51 ± 0.260.521
Chest depth (cm)25.83 b ± 0.4726.67 ab ± 0.3427.28 a ± 0.320.033
Hip width (cm)13.26 ± 0.3513.51 ± 0.2614.10 ± 0.170.069
Cannon circumference (cm)7.92 ± 0.138.08 ± 0.108.01 ± 0.090.621
Note: GZDG—Guanzhong dairy goat; SBWC—Shaanbei white cashmere goat; HNBG—Hainan black goat. Values with different letters (a and b) within the same row differ significantly at p < 0.05.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Wang, K.; Liu, J.; Zhu, H.; Qu, L.; Chen, H.; Lan, X.; Pan, C.; Song, X. An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats. Animals 2019, 9, 1114. https://doi.org/10.3390/ani9121114

AMA Style

Zhang Y, Wang K, Liu J, Zhu H, Qu L, Chen H, Lan X, Pan C, Song X. An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats. Animals. 2019; 9(12):1114. https://doi.org/10.3390/ani9121114

Chicago/Turabian Style

Zhang, Yanghai, Ke Wang, Jinwang Liu, Haijing Zhu, Lei Qu, Hong Chen, Xianyong Lan, Chuanying Pan, and Xiaoyue Song. 2019. "An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats" Animals 9, no. 12: 1114. https://doi.org/10.3390/ani9121114

APA Style

Zhang, Y., Wang, K., Liu, J., Zhu, H., Qu, L., Chen, H., Lan, X., Pan, C., & Song, X. (2019). An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats. Animals, 9(12), 1114. https://doi.org/10.3390/ani9121114

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