An independent confirmation of a quantitative trait locus for milk yield and composition traits on bovine chromosome 26

Several reports have demonstrated that bovine chromosome 26 (BTA26) harbours significant or suggestive quantitative trait loci (QTL) for milk production and composition traits in dairy cattle. Our previous study showed that a C/T substitution in the bovine TCF7L2 gene on BTA26 was significantly linked to QTL for protein yield (PY) in a Canadian dairy cattle population. Actually, this polymorphism was one of the markers derived from a genome‐wide screening of QTL for milk PY using an amplified fragment length polymorphism technique combined with a DNA pooling strategy. In the present study, 990 Holstein bulls with complete genotype and phenotype data from 14 sire families were analysed to confirm, if the QTL effects exist in other populations. Statistical analysis revealed that this marker was significantly associated with PY, protein percentage, milk yield and fat yield (FY) (p < 0.001) in the US Holstein population. These results indicate that this QTL region has a pleiotrophic effect on different milk traits and is portable in different populations.     

Multiple QTL on chromosome six in dairy cattle affecting yield and content traits

A granddaughter design containing five half-sib families from German Holstein–Friesian cattle was subjected to QTL analysis starting from the hypothesis of the existence of more than one QTL on chromosome BTA 6 affecting milk yield, fat yield, protein yield and content of fat and protein. The marker map consisted of 16 microsatellite markers with marker heterozygosity varying from 0.44 to 0.94. Two statistical methods were used: least squares (LS) and residual maximum likelihood (REML) allowing for two QTL simultaneously. The test statistics were calculated in steps of one cM along the chromosome. Significant QTL at the chromosome-wise 5% level according to the permutation test critical value were detected mainly in single families. The results were in conformance with the findings of several previous studies with approximate positions of putative QTL at 49 cM for milk yield, at 70 cM for fat and protein yield, and at 46 cM for protein content. Further QTL positions were suggested mostly for yield traits and protein content in the area of the casein gene cluster at 90…95 cM. The results of the two-QTL model analyses when using LS led to family specific inferences of a second QTL for fat yield and content of protein and fat, partly supported by the epistasis model.     

Quantitative trait loci for milk production and functional traits in two Danish Cattle breeds

Quantitative trait loci (QTL) in Danish Jersey and Danish Red cattle were independently mapped by least squares regression analysis. For Jersey breed, five grandsire families were genotyped for 186 markers on 16 chromosomes (BTAs). Eight traits analysed were milk yield (MY), fat percentage (FP), protein percentage (PP), clinical mastitis (CM), somatic cell score (SCS), maternal stillbirth, maternal calf size (MCS) and maternal calving difficulty. For Red breed, nine grandsire families were genotyped for 166 markers on 18 BTAs. Six traits analysed were MY, FP, PP, CM, SCS and female fertility. Nine and five QTL were detected in Jersey and Red breed, respectively, in across family tests. In Jersey breed, the results indicate QTL for CM and MCS on BTA 3. Additionally, there is an indication of QTL for MCS and FP on BTA 1 and a tentative evidence for a QTL for MY on BTA 26. There is a high risk of detected QTL being false positives. The detected QTL in Jersey breed indicate interesting results from a breeding perspective, but a practical application should await genome-wide association studies.     

Identification of candidate markers on bovine chromosome 14 (BTA14) under milk production trait quantitative trait loci in Holstein

The objective of this study was to identify single-nucleotide polymorphisms using a bovine chromosome 14 high-density SNP panel after accounting for the effect of DGAT1. Linkage disequilibrium information and sire heterozygosity were used to select markers for linkage analysis on bovine chromosome 14 for milk production traits in 321 Holstein animals. Results show putative milk peaks at 42 and 61 cM, both at p<0.10, a fat yield peak at 42 and 63 cM, both at p<0.05; a protein yield peak at 42 (p<0.01) and 84 cM (p<0.05); fat per cent peaks at 3 (p<0.01) and 29 cM (p<0.05), and a protein per cent peak at 4 cM (p<0.05). Once quantitative trait loci positions were established, allele substitution effects for all markers were evaluated using the same statistical model. Overlaying information between quantitative trait loci (QTL) and allele effect analysis enabled the identification (p<0.01) of 20 SNPs under the milk yield QTL, 2 under both of the fat yield peaks, 8 and 9 under the protein yield peaks, 2 and 6 for the fat per cent peaks and 5 for the protein per cent peak. One SNP in particular, ss61514555:A>C, showed association with 3 of the 5 traits: milk (p=1.59E-04), fat (p=6.88E-05) and protein yields (p=5.76E-05). Overall, combining information from linkage disequilibrium, sire heterozygosity and genetic knowledge of traits enabled the characterization of additional markers with significant associations with milk production traits.     

Dissection of genomic correlation matrices of US Holsteins using multivariate factor analysis

The aim of this study was to compare correlation matrices between direct genomic predictions for 31 traits at the genomic and chromosomal levels in US Holstein bulls. Multivariate factor analysis carried out at the genome level identified seven factors associated with conformation, longevity, yield, feet and legs, fat and protein content traits. Some differences were found at the chromosome level; variations in covariance structure on BTA 6, 14, 18 and 20 were interpreted as evidence of segregating QTL for different groups of traits. For example, milk yield and composition tended to join in a single factor on BTA 14, which is known to harbour the DGAT1 locus that affects these traits. Another example was on BTA 18, where a factor strongly correlated with sire calving ease and conformation traits was identified. It is known that in US Holstein, there is a segregating QTL on BTA18 influencing these traits. Moreover, a possible candidate gene for daughter pregnancy rate was suggested for BTA28. The methodology proposed in this study could be used to identify individual chromosomes, which have covariance structures that differ from the overall (whole genome) covariance structure. Such differences can be difficult to detect when a large number of traits are evaluated, and covariances may be affected by QTL that do not have large allele substitution effects.     

Mapping of quantitative trait loci for lactation persistency traits in German Holstein dairy cattle

A whole genome scan to map quantitative trait loci (QTL) for persistency of milk yield (PMY), persistency of fat yield (PFY), persistency of protein yield (PPY) and persistency of milk energy yield (PEY) was performed in a granddaughter design in the German Holstein dairy cattle population. The analysis included 16 paternal half‐sib families with a total of 872 bulls. The analysis was carried out for the first lactation and for the first three lactations combined using univariate weighted multimarker regression. Controlling the false discovery rate across traits and data sets at a level of 0.15 and treating the four persistency traits as different traits revealed 27 significant QTL. A total of 12 chromosomes showed significant QTL effects on a chromosomewise basis. The DGAT1 effect was highly significant for PPY and protein yield. A haplotype analysis using results of previous studies of the same design revealed a co‐segregation of various persistency QTL and QTL affecting health traits like dystocia and stillbirth and functional traits like non‐return rate 90 and somatic cell score.     

Joint analysis of the influence of CYP11B1 and DGAT1 genetic variation on milk production, somatic cell score, conformation, reproduction, and productive lifespan in German Holstein cattle

Recent publications indicate genetic variation in milk production traits on proximal BTA14, which cannot be explained solely with genetic variation in the DGAT1 gene. To elucidate these QTL effects, animals from a German Holstein granddaughter design (18 families, 1,291 sons) were genotyped for CYP11B1 (V30A) and DGAT1 (K232A) polymorphisms. Frequencies of alleles of maternal descent were estimated for CYP11B1(V) (0.776) and DGAT1(K) (0.549). Allele substitution effects (alpha/2) were first calculated for both alleles in separate models and then in a joint model. From the joint analysis, CYP11B1(V) effects on fat content (+0.04%) and protein content (+0.01%) were positive. Effects on milk yield (-82 kg), fat yield (-0.5 kg), and protein yield (-1.9 kg) were negative. Compared with the individual analysis, DGAT1(K) effects on fat content (+0.28%), protein content (+0.06%), and milk yield (-258 kg) were reduced; fat yield (+10.8 kg) was enhanced; and protein yield (-3.8 kg) was reduced. In the joint analysis, allele substitution effects of CYP11B1(V) and DGAT1(K) together explained more of the variation in milk production traits than DGAT1(K) alone. Further significant effects were found for CYP11B1(V) and DGAT1(K) among 6 reproduction traits and 14 conformational traits. These observations indicate a possible negative influence of DGAT1(K) on maternal nonreturn rate, and thus, on length of productive life.     

A primary screen of the bovine genome for quantitative trait loci affecting carcass and growth traits.

A primary genomic screen for quantitative trait loci (QTL) affecting carcass and growth traits was performed by genotyping 238 microsatellite markers on 185 out of 300 total progeny from a Bos indicus x Bos taurus sire mated to Bos taurus cows. The following traits were analyzed for QTL effects: birth weight (BWT), weaning weight (WW), yearling weight (YW), hot carcass weight (HCW), dressing percentage (DP), fat thickness (FT), marbling score (MAR), longissimus muscle area (LMA), rib bone (RibB), rib fat (RibF), and rib muscle (RibM), and the predicted whole carcass traits, retail product yield (RPYD), fat trim yield (FATYD), bone yield (BOYD), retail product weight (RPWT), fat weight (FATWT), and bone weight (BOWT). Data were analyzed by generating an F-statistic profile computed at 1-cM intervals for each chromosome by the regression of phenotype on the conditional probability of receiving the Brahman allele from the sire. There was compelling evidence for a QTL allele of Brahman origin affecting an increase in RibB and a decrease in DP on chromosome 5 (BTA5). Putative QTL at or just below the threshold for genome-wide significance were as follows: an increase in RPYD and component traits on BTA2 and BTA13, an increase in LMA on BTA14, and an increase in BWT on BTA1. Results provided represent a portion of our efforts to identify and characterize QTL affecting carcass and growth traits.     

A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls

The detection and mapping of genetic markers linked to quantitative trait loci (QTL) can be utilized to enhance genetic improvement of livestock populations. With the completion of the bovine genome sequence assembly, single nucleotide polymorphisms (SNP) assays spanning the whole bovine genome and research work on large scale identification, validation and analysis of genotypic variation in cattle has become possible. The objective of the present study was to perform a whole genome scan to identify and map QTL affecting milk production traits and somatic cell scores using linkage disequilibrium (LD) regression and 1536 SNP markers. Three and 18 SNP were found to be associated with only milk yield (MY) at a genome and chromosome wise significance (p < 0.05) level respectively. Among the 21 significant SNP, 16 were in a region reported to have QTL for MY in other dairy cattle populations and while the rest five were new QTL finding. Four SNP out of 21 are significant for the milk production traits (MY, fat yield, protein yield (PY), and milk contents) in the present study. Six and nine SNP were associated with PY at a genome and chromosome wise significant (p < 0.05) level respectively. Three and 17 SNP were found to be associated with FY at a genome and chromosome wise significant (p < 0.05) level. Five and seven SNP were mapped with somatic cell score at a genome and chromosome wise significant (p < 0.05) level respectively. The results of this study have revealed QTL for MY, PY, protein percentage, FY, fat percentage, somatic cell score and persistency of milk in the Canadian dairy cattle population. The chromosome regions identified in this study should be further investigated to potentially identify the causative mutations underlying the QTL.     

Effects of DGAT1 variants on milk production traits in German cattle breeds

Various QTL mapping experiments led to the detection of a QTL in the centromeric region of cattle chromosome 14 that had a major effect on the fat content of milk. Recently, the gene encoding diacylglycerol O-acyltransferase (DGAT1) was proposed to be a positional and functional candidate for this trait. This study investigated the effects of a nonconservative lysine to alanine (K232A) substitution in DGAT1, which very likely represents the causal mutation, on milk production traits. Existing granddaughter designs for Fleckvieh and German Holstein, the two major dairy/dual-purpose breeds in Germany, were used to estimate allele frequencies and gene substitution effects for milk, fat, and protein yield, as well as fat and protein content. A restriction fragment length polymorphism assay was applied to diagnose the K232A substitution in DGAT1. Estimates of the allele frequencies for the lysine-encoding variant were based on maternally inherited alleles in sons and amounted to 0.072 for Fleckvieh and 0.548 for German Holstein. Effects of DGAT1 variants on content traits were pronounced; estimates of the gene substitution effect for the lysine-encoding variant were 0.35 and 0.28% for fat content and 0.10 and 0.06% for protein content in Fleckvieh and German Holstein, respectively. Conversely, negative effects of the lysine variant of -242 to -180 kg for Fleckvieh and -260 to -320 kg for German Holstein were revealed for milk yield from first to third lactation, resulting in enhanced fat yield of 7.5 to 14.8 kg in Fleckvieh and 7.6 to 10.7 kg in German Holstein. For protein yield, however, mainly negative effects of -3.6 to 0.2 kg in Fleckvieh and -4.8 to -5.2 kg in German Holstein were observed. Pearson correlations between residuals of milk yield and content traits were decreased when omitting DGAT1 effects in the analysis, thereby indicating that DGAT1 contributes to negative correlations between these traits. Molecular tests allow for the direct selection among variants; however, the benefits of the alternative alleles depend on economic weights given to the different milk production traits in the breeding goal.     

Fine mapping of a calving QTL on Bos taurus autosome 18 in Holstein cattle

Decreased calving performance not only directly impacts the economic efficiency of dairy cattle farming but also influences public concern for animal welfare. Previous studies have revealed a QTL on Bos taurus autosome (BTA) 18 that has a large effect on calving traits in Holstein cattle. In this study, fine mapping of this QTL was performed using imputed high‐density SNP chip (HD) genotypes followed by imputed next‐generation sequencing (NGS) variants. BTA18 was scanned for seven direct calving traits in 6113 bulls with imputed HD genotypes. SNP rs136283363 (BTA18: 57 548 213) was consistently the most significantly associated SNP across all seven traits [e.g. p‐value = 2.04 × 10^?59 for birth index (BI)]. To finely map the QTL region and to explore pleiotropic effects, we studied NGS variants within the targeted region (BTA18: 57 321 450–57 625 355) for associations with direct calving traits and with three conformation traits. Significant variants were prioritized, and their biological relevance to the traits was interpreted. Considering their functional relationships with direct calving traits, SIGLEC12, CD33 and CEACAM18 were proposed as candidate genes. In addition, pleiotropic effects of this QTL region on direct calving traits and conformation traits were observed. However, the extent of linkage disequilibrium combined with the lack of complete annotation and potential errors in the Bos taurus genome assembly hampered our efforts to pinpoint the causal mutation.     

Detection of quantitative trait loci for growth and beef carcass fatness traits in a cross between Bos taurus (Angus) and Bos indicus (Brahman) cattle

This study was conducted to detect quantitative trait loci (QTL) affecting growth and beef carcass fatness traits in an experimental population of Angus and Brahman crossbreds. The three-generation mapping population was generated with 602 progeny from 29 reciprocal backcross and three F2 full-sib families, and 417 genetic markers were used to produce a sex-averaged map of the 29 autosomes spanning 2,642.5 Kosambi cM. Alternative interval-mapping approaches were applied under line-cross (LC) and random infinite alleles (RA) models to detect QTL segregating between and within breeds. A total of 35 QTL (five with genomewide significant and 30 with suggestive evidence for linkage) were found on 19 chromosomes. One QTL affecting yearling weight was found with genomewide significant evidence for linkage in the interstitial region of bovine autosome (BTA) 1, and an additional 19 QTL were detected with suggestive evidence for linkage under the LC model. Many of these QTL had a dominant (complete or overdominant) mode of gene action, and only a few of the QTL were primarily additive, which reflects the fact that heterosis for growth is known to be appreciable in crosses among Brahman and British breeds. Four QTL affecting growth were detected with genomewide significant evidence for linkage under the RA model on BTA 2 and BTA 6 for birth weight, BTA 5 for yearling weight, and BTA 23 for hot carcass weight. An additional 11 QTL were detected with suggestive evidence for linkage under the RA model. None of the QTL (except for yearling weight on BTA 5) detected under the RA model were found by the LC analyses, suggesting the segregation of alternate alleles within one or both of the parental breeds. Our results reveal the utility of implementing both the LC and RA models to detect dominant QTL and also QTL with similar allele frequency distributions within parental breeds.     

中国荷斯坦牛RPL23A和ACACB基因对产奶性状的遗传效应分析

作者所在团队前期通过奶牛乳腺上皮组织转录组测序及荷斯坦公牛全基因组重测序研究发现RPL23A和ACACB基因是奶牛乳蛋白和乳脂性状的候选功能基因,本研究旨在探究这两个基因是否对奶牛产奶性状具有显著遗传效应。以北京地区7个牧场的1059头中国荷斯坦母牛为试验群体,采集尾根静脉血并提取基因组DNA,通过飞行时间质谱方法检测SNP位点基因型,利用SAS9.4软件的MIXED过程进行关联分析。结果表明,RPL23A基因的SNP位点g.20146771C>T与第1泌乳期5个产奶性状达到显著或极显著关联(P=0.0001~0.0416),其优势等位基因为T;ACACB基因的g.63878254T>C位点与第1泌乳期产奶量、乳脂量和乳蛋白量呈极显著关联(P<0.01),其优势等位基因为C;g.63962768G>A位点与第1泌乳期产奶量、乳脂量、乳脂率和乳蛋白率关联显著或极显著(P=0.0001~0.0391),其优势等位基因为A。综上,RPL23A基因主要影响中国荷斯坦牛产奶量和乳蛋白,ACACB基因对产奶量和乳脂具有显著遗传效应,3个SNP位点可考虑作为遗传标记用于标记辅助选择培育奶牛高乳蛋白乳脂新品系和选育提高。     

Evaluation of association between polymorphism within the thyroglobulin gene and milk production traits in dairy cattle

In dairy cattle, many studies have reported quantitative trait loci (QTL) on the centromeric end of chromosome 14 that affect milk production traits. One of the candidate genes in this QTL region – thyroglobulin (TG) – was previously found to be significantly associated with marbling in beef cattle. Thus, based on QTL studies in dairy cattle and because of possible effects of this gene on fat metabolism, we investigated the association of TG with milk yield and composition in Holstein dairy cattle. A total of 1279 bulls from the Cooperative Dairy DNA Repository Holstein population were genotyped for a single nucleotide polymorphism in TG used previously in beef cattle studies. Analysis of 29 sire families showed no significant association between TG variants and milk production traits. Within‐sire family analysis suggests that TG is neither the responsible gene nor a genetic marker in association with milk production traits.     

Bovine quantitative trait loci analysis for growth, carcass, and meat quality traits in an F2 population from a cross between Japanese Black and Limousin

A genome-wide scan for QTL affecting economically important traits in beef production was performed using an F(2) resource family from a Japanese Black x Limousin cross, where 186 F(2) animals were measured for growth, carcass, and meat-quality traits. All family members were genotyped for 313 informative microsatellite markers that spanned 2,382 cM of bovine autosomes. The centromeric region of BTA2 contained significant QTL (i.e., exceeding the genome-wide 5% threshold) for 5 carcass grading traits [LM area, beef marbling standards (BMS) number, luster, quality grade, and firmness), 8 computer image analysis (CIA) traits [LM lean area, ratio of fat area (RFA) to LM area, LM area, RFA to musculus (M.) trapezius area, M. trapezius lean area, M. semispinalis lean area, RFA to M. semispinalis area, and RFA to M. semispinalis capitis area], and 5 meat quality traits (contents of CP, crude fat, moisture, C16:1, and C18:2 of LM). A significant QTL for withers height was detected at 80.3 cM on BTA5. We detected significant QTL for the C14:0 content in backfat and C14:0 and C14:1 content in intermuscular fat around the 62.3 to 71.0 cM region on BTA19 and for C14:0, C14:1, C18:1, and C16:0 content and ratio of total unsaturated fatty acid content to total SFA content in intramuscular fat at 2 different regions on BTA19 (41.1 cM for C14:1 and 62.3 cM for the other 4 traits). Overall, we identified 9 significant QTL regions controlling 27 traits with genome-wide significance of 5%; of these, 22 traits exceeded the 1% genome-wide threshold. Some of the QTL affecting meat quality traits detected in this study might be the same QTL as previously reported. The QTL we identified need to be validated in commercial Japanese Black cattle populations.     

Non-additive effects on milk production in Czech dairy cows

Crossbreeding effects on milk production traits of Czech dual‐purpose and dairy cattle breeds were estimated. Nearly 370 000 cows with known gene proportions from Czech Pied, Ayrshire or Holstein cattle were selected from the national milk recording data base. Single‐trait animal models were calculated for milk, fat and protein yield, fat and protein content. The model of Dickerson including additive, additive maternal, heterotic and recombination effects was used for the part of the animal model describing the crossbreeding effects in all calculations. For milk yield, the additive genetic effect (defined as deviation from Czech Pied cattle) was 850–900 kg for Holstein and 240–480 kg for Ayrshire. The maternal effects were low and negative. Low significant positive heterotic effects were observed being up to approximately 100 kg for Czech Pied × Holstein. The recombination effects were negative and statistically significant for Czech Pied × Holstein. The results for fat and protein yield were similar to the results for milk yield. For fat and protein content, nearly no statistically significant crossbreeding effects were found.     

Searching for DNA markers for milk production and composition on chromosome 6 in sheep

Several milk protein polymorphisms are potential tools for selection in dairy ruminants. However, research results for dairy sheep are not as conclusive as those for goats or cattle and are often controversial. The main objective of this study was to find and later use molecular genetic markers in selection to improve milk production and milk composition in Awassi ewes. Chromosome 6 was chosen because several studies have reported the presence of significant quantitative trait loci (QTL) affecting milk production traits on ovine and bovine chromosome 6. Altogether, genotypes for 13 microsatellite loci were determined for 258 ewes, which were purebred Awassi or Awassi-Merino crosses. Phenotypic data were lactation yield of milk, milk fat, protein and lactose (kg), average milk protein and fat percentage and average somatic cell count. Five out of the 13 microsatellites showed significant association with at least one of the examined traits.     

Theoretical efficiency of multiple-trait quantitative trait loci-assisted selection

The effectiveness of five selection methods for genetic improvement of net merit comprising trait 1 of low heritability (h² = 0.1) and trait 2 of high heritability (h² = 0.4) was examined: (i) two‐trait quantitative trait loci (QTL)‐assisted selection; (ii) partial QTL‐assisted selection based on trait 1; (iii) partial QTL‐assisted selection based on trait 2; (iv) QTL‐only selection; and (v) conventional selection index without QTL information. These selection methods were compared under 72 scenarios with different combinations of the relative economic weights, the genetic correlations between traits, the ratio of QTL variance to total genetic variance of the trait, and the ratio of genetic variances between traits. The results suggest that the detection of QTL for multiple‐trait QTL‐assisted selection is more important when the index traits are negatively correlated than when they are positively correlated. In contrast to literature reports that single‐trait marker‐assisted selection (MAS) is the most efficient for low heritability traits, this study found that the identified QTL of the low heritability trait contributed negligibly to total response in net merit. This is because multiple‐trait QTL‐assisted selection is designed to maximize total net merit rather than the genetic response of the individual index trait as in the case of single‐trait MAS. Therefore, it is not economical to identify the QTL of the low heritability traits for the improvement of total net merit. The efficient, cost‐effective selection strategy is to identify the QTL of the moderate or high heritability traits of the QTL‐assisted selection index to facilitate total economic returns. Detection of the QTL of the low h² traits for the QTL‐assisted index selection is justified when the low h² traits have high negative genetic correlation with the other index traits and/or when both economic weights and genetic variances of the low h² traits are larger as compared to the other index traits of higher h². This study deals with theoretical efficiency of QTL‐assisted selection, but the same principle applies to SNP‐based genomic selection when the proportion of the genetic variance ‘explained by the identified QTLs’ in this study is replaced by ‘explained by SNPs’.     

中国荷斯坦牛PIN1基因对产奶性状的遗传效应分析

前期研究通过荷斯坦公牛全基因组重测序鉴定到17个奶牛产奶性状候选功能基因,其中,肽基脯氨酸顺反异构酶基因(PIN1)参与甘油三酯代谢、甘油磷脂代谢以及mTOR信号通路,且位于产奶量和乳蛋白量性状QTL区间。为进一步系统分析PIN1基因是否对奶牛产奶性状具有遗传效应,本实验基于40头公牛的基因组DNA混池,采用PCR产物直接测序法对PIN1基因的全部编码区以及上下游调控区2000 bp进行扫描,在内含子2检测到1个SNP位点7:g.14432394G>A,A、G等位基因频率分别为0.4797和0.5203。采用靶向测序基因型技术对北京地区987头中国荷斯坦母牛进行个体基因型检测,对SNP位点7:g.14432394G>A与5个产奶性状进行关联分析。结果表明:在第1泌乳期,SNP 7:g.14432394G>A与产奶量、乳脂量、乳蛋白量和乳蛋白率呈显著或极显著关联(P=0.0001~0.0493);在第2泌乳期,SNP与产奶量、乳脂量、乳脂率和乳蛋白量呈显著或极显著关联(P=0.0001~0.0104);SNP位点7:g.14432394G>A对产奶量、乳脂量、乳蛋白量和乳蛋白率的加性效应或等位基因替代效应均达到显著或极显著。综上,PIN1基因对中国荷斯坦牛的产奶量和乳蛋白、乳脂性状具有显著遗传效应,可作为遗传标记用于基因组选择,以加快遗传进展。