Associations between variants of the HALgene and milk production traits in Chinese Holstein cows
BMC Genetics volume 15, Article number: 125 (2014)
The histidine ammonia-lyse gene (HAL) encodes the histidine ammonia-lyase, which catalyzes the first reaction of histidine catabolism. In our previous genome-wide association study in Chinese Holstein cows to identify genetic variants affecting milk production traits, a SNP (rs41647754) located 357 bp upstream of HAL, was found to be significantly associated with milk yield and milk protein yield. In addition, the HAL gene resides within the reported QTLs for milk production traits. The aims of this study were to identify genetic variants in HAL and to test the association between these variants and milk production traits.
Fifteen SNPs were identified within the regions under study of the HAL gene, including three coding mutations, seven intronic mutations, one promoter region mutation, and four 3′UTR mutations. Nine of these identified SNPs were chosen for subsequent genotyping and association analyses. Our results showed that five SNP markers (ss974768522, ss974768525, ss974768531, ss974768533 and ss974768534) were significantly associated with one or more milk production traits. Haplotype analysis showed that two haplotype blocks were significantly associated with milk yield and milk protein yield, providing additional support for the association between HAL variants and milk production traits in dairy cows (P < 0.05).
Our study shows evidence of significant associations between SNPs within the HAL gene and milk production traits in Chinese Holstein cows, indicating the potential role of HAL variants in these traits. These identified SNPs may serve as genetic markers used in genomic selection schemes to accelerate the genetic gains of milk production traits in dairy cattle.
The dairy industry plays a key role in contributing to agricultural economy and it is the main financial source of dairy farmers. Hence an improvement in milk production traits continues to be the most profitable breeding goal. Although classical selection approaches have obtained great genetic gains of farm animal performance, the implementation of selection programs containing molecular information could facilitate the improvement in selection response of animals .
With the progress in DNA-based marker technology, the identification of genomic regions (i.e. quantitative trait loci, QTLs) that associated with complex traits has been available . Incorporating the detected QTLs into genetic evaluation makes it possible to evidently improve the accuracy of selection, thus accelerating the genetic improvement of animal productivity. Over the last decades, numerous investigations regarding identification of QTLs for milk production traits in dairy cattle have been reported -. Advances in detecting causal genes for complex traits are slow as the linkage mapping results in large confidence intervals ,. Particularly, the region that a QTL is mapped might contain too many plausible candidate genes. Until recently, only few candidate genes responsible for variation in complex traits were discovered and functionally confirmed in dairy cattle, i.e., the diacylglycerol acyltransferase 1 (DGAT1) gene  and the growth hormone receptor (GHR) gene .
High throughput SNP genotyping technology has made it feasible to screen for mutations associated with complex traits by genome-wide association study (GWAS). The advantages of GWAS are an increase in the power to detect underlying variants in complex traits and in simplifying the discovery of causal variants . The GWAS has been widely recognized as a valid approach for gene discovery and successfully identified many genes associated with economically important traits in several domestic animal species -.
In our previous GWAS study in Chinese Holstein dairy cows to identify genetic variants affecting milk production traits, we not only identified some previously reported functional genes but also discovered 105 SNPs that significantly associated with milk production traits . Among these SNPs, a SNP, BFGL-NGS-110018 (P = 1.50E-07, n = 1815), located 357 bp upstream of the histidine ammonia-lyse (HAL) gene was highly significant in our GWAS results with effects on milk yield and milk protein yield. The HAL gene located on bovine chromosome (BTA) 5 encodes a cytosolic enzyme known as histidase or histidinase, which catalyzes the first reaction of histidine catabolism, the nonoxidative deamination of histidine to urocanic acid . Histidine is the first-limiting amino acid for milk protein synthesis in dairy cows, the infusion of which was related to significant increases in milk and milk protein yields . In addition, the HAL gene is located within the reported QTLs for milk production traits -. These findings strongly suggest that HAL is a promising positional and functional candidate gene for milk production traits in dairy cattle.
In this study, we explored genetic variants of the HAL gene and tested the association between HAL variants and milk production traits in Chinese Holstein cows, and performed functional prediction analysis of the identified variants. Significant associations were observed between a subset of HAL variants and the milk production traits of interest. Our study firstly identified the genetic variants within the HAL gene in dairy cattle and the information of these SNPs markers will be available for genomic selection in dairy cattle breeding programs.
The whole study protocols for collection of the tissue samples of experimental individuals were reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of China Agricultural University.
Animals and phenotypic data
Thirteen unrelated sires and their corresponding 638 daughters with known milk production phenotypes were collected to construct the study population. Family size ranges from 5 to 136 daughters with an average of about 49 daughters per sire. Daughters were from 15 Chinese Holstein cattle farms in Beijing, China, where regular and standard performance testing (Dairy Herd Improvement, DHI) has been conducted since 1999. The official up to date estimated breeding values (EBVs) of five milk production traits, including first lactation 305-d milk yield (MY), first lactation 305-d fat yield (FY), first lactation 305-d protein yield (PY), first lactation average 305-d fat percentage (FP), and first lactation average 305-d protein percentage (PP) were considered as phenotypic observations used for the subsequent analysis. These EBVs were calculated by the official Dairy Data Center of China based on the genetic parameters that estimated via the complete DHI data of Chinese Holstein cattle population.
Genomic DNA extraction
Genomic DNA was isolated from blood sample of the daughters by QIAamp DNA blood kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions and from semen sample of the sires using the routine salt-out procedure. The quality and quantity of extracted DNA were measured using the spectrophotometer and gel electrophoresis.
SNP identification and genotyping
A DNA pool was constructed by mixing equal amounts of DNA from the 13 unrelated sires to identify variants of the HAL gene. All exons and their adjacent intronic sequences were targeted for selective amplification by PCR. Sixteen pairs of nucleotide primers (Additional file 1: Table S1) targeting the regions of interest were designed using Primer Premier 5.0  based on the reference sequence NC_007303.5 and synthesized by Sangon Biotech Co. Ltd. (Shanghai, China).
PCR amplifications were performed in a total volume of 25 μL. The following reagents were used for amplification: 50-100 ng of genomic DNA, 10 pmol of each primer, 5 mM of dNTP mix, 2.5 μL of 10× PCR buffer, 0.625 U Taq DNA Polymerase (Takara Biotechnology (Dalian) Co., Ltd., Dalian, China). The PCR reaction conditions were as follows: 95°C for 5 min, followed by 35 cycles of 94°C for 30’sec, annealing from 46°C to 58°C for 30’sec, 72°C for 35’sec, and a final extension at 72°C for 10 min. The PCR products were directly sequenced using an ABI3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA).
We chose nine out of these identified SNPs focusing primarily on SNPs located in the coding regions for subsequent genotyping in 638 Chinese Holstein cows using the iPLEX MassARRAY system (Sequenom Inc., San Diego, California, USA).
To study the effects of HAL variants on milk production traits, we performed a single locus-based regression analysis. To this end, we utilized the same analytical method as employed in our previous GWAS study . A linear mixed regression model was adapted as follows:
where y is the vector of EBVs of all daughters, μ is the overall mean, b is the regression coefficient of EBVs on SNP genotypes. x is the vector of the SNP genotype predictors set as 0, 1 or 2, representing the three genotypes 11, 12 and 22 respectively (setting 1 as the minor frequency allele). a is the vector of residual polygenetic effects with (where A is the additive genetic relationship matrix according to the pedigree data regarding all experimental individuals and is the additive variance), Z is the incidence matrix of a, and e is the vector of residual errors with (where W is a diagonal matrix with the diagonal elements equal to 1/RELij, of which RELij is the reliability of EBV of daughter i in family j, and is the residual error variance). For each SNP, the estimate of b and the corresponding sampling variance were obtained by solving the mixed model equations. Afterwards, a Wald chi-squared statistic with df = 1 was applied to determine whether the SNPs were associated with the milk production traits considered in this study.
Pairwise linkage disequilibrium (LD) was calculated and plotted by using Haploview . Haplotype block structure was defined for marker pairs showing D’ > 0.75, with the fraction of strong linkage disequilibrium informative comparisons set at 0.9. For haplotype-based association analysis, considering multiple loci in high LD, we extended the original haplotype trend regression (HTR) approach  by allowing random effects in the regression model. The haplotype trend regression model with polygenic random effects, as described in our previous work , was as following:
where y is the vector of EBV of all experimental daughters, μ is the overall mean. 1 is the n-dimensional vector with all elements equal to 1, h is the haplotype fixed effect vector with elements hi (i = 1, 2, ..., k) being the effect of haplotype of all k distinct haplotypes within the haplotype block. X is the n × k indicator matrix with the same pattern as that in , a is the vector of the residual polygenetic effects with (where A is the additive genetic relationship matrix based on the pedigree data regarding all experimental individuals and is the additive variance), Z is the incidence matrix of a, and e is the vector of residual errors with (where is the residual error variance and the weight matrix W is a diagonal matrix with each diagonal element equal to the reciprocal of the reliability of the EBV of individual i). With regard to each haplotype block, the estimate of the haplotype effects vector h and the corresponding sampling variance Var(h'h) were acquired via solving the mixed model equations, and a Wald chi-squared statistic with df = k-1 (k is the number of distinct haplotypes in the experimental populations within the region of haplotype block investigated) was constructed to examine whether the haplotype block was associated with the trait.
In the haplotype analysis, the haplotype with frequency > 5% was regarded as a distinguishable haplotype, and those haplotypes each with relative frequency < 5% were pooled into a single group. For both analyses, false discovery rate (FDR) correction was performed to correct for multiple testing. An adjusted P-value of ≤ 0.05 was considered statistically significant.
Hardy-Weinberg equilibrium (HWE) test was performed on each SNP. We calculated the SNP allele frequencies by determining their proportions in the samples and then estimated the expected genotype numbers using the expected genotype frequencies under HWE. Then, a goodness-of-fit test (Chi-square) was used to compare the expected and observed genotype numbers, choosing a significance level of 0.05.
Quantification of gene expression
Total RNA was extracted from the mammary gland tissues using TRIzol Reagent (Life technology, Carlsbad, CA, USA) following the manufacturer’s protocols. The tissue samples were collected from eight genetically unrelated Chinese Holstein cows in late lactation (approximately 300 days) after slaughter within 30 min and frozen in liquid nitrogen and then stored at -80°C. The quality of extracted RNA was assessed by electrophoresis in 1% agarose gel and quantified with the NanoDropTM 2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). RNA was then purified and reversely transcribed into cDNA using PrimerScript® RT reagent Kit with gDNA Eraser (Takara Biotechnology (Dalian) Co., Ltd.) following the manufacturer’s instructions. Quantities of mRNA were then measured with real-time quantitative PCR (qPCR) using a LightCycler® 480 Real-Time PCR System (Roche, Hercules, CA, USA). The qPCR assays were performed with a volume of 20 μL containing 10 μL SYBR Green Mixture, 7 μL deionized water, 1 μL template of cDNA, 1 μL of each primer. The reaction conditions were as follows: 95°C for 5 min, 45 cycles of 95°C for 10’sec, 60°C for 10’sec, 72°C for 10’sec. The glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH) was used as a reference to normalize gene expression. Primers for the amplification of both HAL and GAPDH genes are presented in Additional file 2: Table S2. Each qPCR assay was carried out in triplicate. The relative gene expression was calculated by using the 2-ΔΔCt method .
Bioinformatics analysis of identified SNPs
The potential impact of the amino acid transition on the structure or function of the protein was predicted using the two web server tools SIFT (http://sift.jcvi.org/www/SIFT.html) and PolyPhen (http://genetics.bwh.harvard.edu/pph2/). The search for potential transcription factor binding site (TFBS) for polymorphisms around promoter region was performed by using the online tool TFSEARCH (http://www.cbrc.jp/research/db/TFSEARCH). Effects of 3′UTR variants on microRNA binding sites were searched with STarMir .
SNPs within the HALgene in Chinese Holstein cows
A total of fifteen SNPs in the HAL gene were discovered in this study (Table 1). Among these identified SNPs, seven were found within introns and one within the promoter region. In addition, seven SNPs were exonic, of which two were synonymous substitutions and one was a non-synonymous switch resulting in an amino acid substitution (Gly to Glu) at residue 228, and the remaining four were located in the 3′UTR.
SNPs in coding parts of exons are considered to be more likely to be functional than those in other genic regions. We therefore chose nine SNPs mainly located in the coding regions for genotyping (Figure 1). Locations and allele frequencies of the nine SNPs are shown in Table 2.
Compared with the SNP loci contained in the 700 K BovineHD SNP array, none of these identified SNPs are currently present. The details of all identified SNPs have been submitted to dbSNP (http://www.ncbi.nlm.nih.gov/SNP/) and will be publicly available (submission numbers from ss974768522 to ss974768536) in dbSNP Build (B142).
Our findings showed that significant associations exist between five markers (ss974768522, ss974768525, ss974768531, ss974768533 and ss974768534) and one or more milk production traits (FDR P < 0.05). Out of these associations, the synonymous SNP ss974768525 was associated with FY, with the genotype CC showing higher levels in FY. Two 3′UTR variants, ss974768533 and ss974768534, showed significant associations with PP with the genotype CC linked to higher increases in PP. The promoter variant ss974768522 was significantly associated with PY with the heterozygous genotype CT related to higher increases in this trait. The intronic mutation ss974768531 was also found to be associated with MY, FY, PY and FP. A marginal significant association between SNP ss974768527 and PY was observed (P = 0.0447, not significant after correction for multiple testing). Results of SNP marker associations with milk production traits are presented in Additional file 3: Table S3.
Two haplotype blocks were constructed (Figure 2). The block 1 included 6 SNPs, which formed 5 major haplotypes (> 5%) in our resource population. The block 2 consisted of 3 major haplotypes. The constitution of the haplotypes and their frequencies are listed in Table 3. The haplotype blocks containing these SNPs were also tested for associations with milk production traits. Both the two blocks were associated with MY and PY (P < 0.05) (Table 3). These findings produced additional support for the association between these SNP markers and milk production traits.
Expression analysis of the HALgene
As the possible effects of non-synonymous mutations on gene expression by altering mRNA stability , we measured HAL expression in mammary gland tissues with different genotypes at the locus of SNP ss974768527. Results indicated that individuals with the GG genotype have greater expression level (n = 4, mean relative expression = 0.331) than those with the GA genotype (n = 4, mean relative expression = 0.014) (Figure 3a). We also analyzed the possible effects of SNP ss974768522 located in the regulatory region on gene expression and observed that individuals with the CC genotype have higher expression level (n = 4, mean relative expression = 0.368) than those with the CT genotype (n = 4, mean relative expression = 0.046) (Figure 3b). As the SNP ss974768527 was found to be a non-significantly associated mutation, the possible effects of SNP ss974768522 may be responsible for the observed differential gene expression results. For both SNPs, as the relatively lower frequency of the homozygous mutant genotype, we failed to obtain all genotypes in the tested samples.
Prediction of potentially functional SNPs
Effects of the detected non-synonymous mutation ss974768527 were predicted to be “tolerated” and “benign” by the SIFT and PolyPhen programs respectively, implying that the related protein structure and function may not be affected.
We also examined the impact of promoter variant ss974768522 on the binding site for transcriptional factors. Results showed that an 8 bp sequence containing the C allele was within the transcription factor myeloid zinc finger 1 (MZF1) binding motif, which was disrupted by the presence of the T allele (Figure 4).
In addition, the STarMir analysis showed that none of the four SNPs within 3′UTR is situated within the binding site for mature microRNAs in cattle.
In this study, based on the reported QTLs for milk production traits - and our previous GWAS findings , we explored the genetic variants of the HAL gene and tested their association with milk production traits in Chinese Holstein cows. Our results demonstrated that significant associations exist between HAL variants and milk production traits, which provide important SNP marker information that can be considered for genetic improvement in dairy breeding schemes.
The identified non-synonymous SNP ss974768527 caused an amino acid substitution (Gly to Glu) at residue 228 of the protein. These two amino acids are both polar amino acids having aliphatic side chain groups. This substitution does not seem to influence the structure and function of the related protein, which is further supported by the prediction of functional SNPs. Therefore this non-functional change may explain in part why no significant association was observed between this SNP and milk traits under study.
As synonymous mutations do not change the amino acid sequence of the encoded proteins, this type of genetic variation has generally been considered as silent mutation with no effect on protein structure and function. However, accumulating evidence indicates that synonymous SNPs could have functional consequences . Considerable studies have identified the significant contribution of synonymous SNPs to human disease -. Synonymous SNPs could influence protein expression and enzymatic activity through modifying mRNA stability ,. In addition, synonymous codons may have different codon usage, and therefore lead to varied gene expression . According to the codon usage table of cattle (http://www.kazusa.or.jp/codon/), the two synonymous SNPs detected here resulted in moderate alterations in codon usage frequency, for SNP (ss974768523) with AAT changed to AAC from 14.7 to 21.4 per thousand, for SNP (ss974768525) with ATC changed to ATT from 23.3 to 14.6 per thousand. Therefore, the two synonymous mutations identified in HAL might have the potential to influence gene expression and thus function of gene product.
One of the most important mechanisms by which SNPs in the regulatory region of a gene can affect transcription rate is by abolishing or introducing a TFBS . In this study, the C/T substitution in the promoter region (ss974768522) of HAL may lead to the disappearance of binding sites of transcription factor MZF1, which pertains to the kruppel family of zinc-finger proteins. It has been identified as a bi-functional transcription regulator, repressing transcription in non-hematopoietic cells and activating transcription in hematopoietic cells . Therefore, the possible effects on MZF1 may be responsible for the observed differential gene expression results.
In the present study, the EBVs of daughters were treated as phenotypic observations for association analyses. In addition to EBVs, de-regressed EBVs and yield deviation of individual were also often considered as phenotypes in the association study relevant to milk production traits ,. Considering the effects of such phenotypes on QTL mapping and marker-assisted selection, none of them has a remarkably greater advantage over the others . Consequently, we only used the EBVs for association analysis. It is also notable that we fulfilled the association analyses with the method of SNP by SNP individually by regressing the phenotypes of a single trait on either the genotypes of a SNP. As the multi-marker analysis containing multiple marker information was generally thought more powerful than single marker analysis , we further performed the haplotype-based association analysis to investigate the association of HAL variants with milk production traits.
SNP-based association analysis revealed that five out of the nine tested SNPs were associated with the milk production traits under analysis. Haplotype-based association analysis provided further evidence for these associations. As these identified SNPs are currently absent in the 700 K BovineHD SNP array, our findings will be helpful for the enrichment of genetic marker information used in genomic selection program of dairy cattle. The information of these identified SNPs could be integrated into the panel of high density SNP arrays used in genomic selection of dairy cattle breeding schemes to improve the frequency of the genetic marker which is positively correlated with the milk production traits of interest.
In summary, we firstly explored genetic variants within the HAL gene and uncovered associations between the tested SNPs and milk production traits in Chinese Holstein cows. Results of this study could be used as a first step to explore HAL variants potentially responsible for variation in milk production traits and further studies will be required to replicate these association findings in different cattle populations. The information of these variants will be available for genomic selection to accelerate the genetic gains in Chinese Holstein cattle breeding programs.
Carta A, Casu S, Salaris S: Invited review: current state of genetic improvement in dairy sheep. J Dairy Sci. 2009, 92 (12): 5814-5833. 10.3168/jds.2009-2479.
Liu Z, Cordes J: DNA marker technologies and their applications in aquaculture genetics. Aquaculture. 2004, 238 (1): 1-37. 10.1016/j.aquaculture.2004.05.027.
Ashwell MS, Heyen DW, Sonstegard TS, Van Tassell CP, Da Y, VanRaden PM, Ron M, Weller JI, Lewin HA: Detection of quantitative trait loci affecting milk production, health, and reproductive traits in Holstein cattle. J Dairy Sci. 2004, 87 (2): 468-475. 10.3168/jds.S0022-0302(04)73186-0.
Chen HY, Zhang Q, Yin CC, Wang CK, Gong WJ, Mei G: Detection of quantitative trait loci affecting milk production traits on bovine chromosome 6 in a Chinese Holstein population by the daughter design. J Dairy Sci. 2006, 89 (2): 782-790. 10.3168/jds.S0022-0302(06)72140-3.
Kolbehdari D, Wang Z, Grant JR, Murdoch B, Prasad A, Xiu Z, Marques E, Stothard P, Moore SS: A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls. J Anim Breed Genet. 2009, 126 (3): 216-227. 10.1111/j.1439-0388.2008.00793.x.
Andersson L, Georges M: Domestic-animal genomics: deciphering the genetics of complex traits. Nat Rev Genet. 2004, 5 (3): 202-212. 10.1038/nrg1294.
Goddard ME, Hayes BJ: Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet. 2009, 10 (6): 381-391. 10.1038/nrg2575.
Grisart B, Coppieters W, Farnir F, Karim L, Ford C, Berzi P, Cambisano N, Mni M, Reid S, Simon P, Spelman R, Georges M, Snell R: Positional candidate cloning of a QTL in dairy cattle: identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition. Genome Res. 2002, 12 (2): 222-231. 10.1101/gr.224202.
Blott S, Kim JJ, Moisio S, Schmidt-Kuntzel A, Cornet A, Berzi P, Cambisano N, Ford C, Grisart B, Johnson D, Karim L, Simon P, Snell R, Spelman R, Wong J, Vilkki J, Georges M, Farnir F, Coppieters W: Molecular dissection of a quantitative trait locus: a phenylalanine-to-tyrosine substitution in the transmembrane domain of the bovine growth hormone receptor is associated with a major effect on milk yield and composition. Genetics. 2003, 163 (1): 253-266.
Hirschhorn JN, Daly MJ: Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005, 6 (2): 95-108. 10.1038/nrg1521.
Jiang L, Liu J, Sun D, Ma P, Ding X, Yu Y, Zhang Q: Genome wide association studies for milk production traits in Chinese Holstein population. PLoS One. 2010, 5 (10): e13661-10.1371/journal.pone.0013661.
Fontanesi L, Schiavo G, Galimberti G, Calo DG, Scotti E, Martelli PL, Buttazzoni L, Casadio R, Russo V: A genome wide association study for backfat thickness in Italian Large White pigs highlights new regions affecting fat deposition including neuronal genes. BMC Genomics. 2012, 13: 583-10.1186/1471-2164-13-583.
Xie L, Luo C, Zhang C, Zhang R, Tang J, Nie Q, Ma L, Hu X, Li N, Da Y, Zhang X: Genome-wide association study identified a narrow chromosome 1 region associated with chicken growth traits. PLoS One. 2012, 7 (2): e30910-10.1371/journal.pone.0030910.
Suchi M, Sano H, Mizuno H, Wada Y: Molecular cloning and structural characterization of the human histidase gene (HAL). Genomics. 1995, 29 (1): 98-104. 10.1006/geno.1995.1219.
Huhtanen P, Vanhatalo A, Varvikko T: Effects of abomasal infusions of histidine, glucose, and leucine on milk production and plasma metabolites of dairy cows fed grass silage diets. J Dairy Sci. 2002, 85 (1): 204-216. 10.3168/jds.S0022-0302(02)74069-1.
Olsen HG, Gomez-Raya L, Vage DI, Olsaker I, Klungland H, Svendsen M, Adnoy T, Sabry A, Klemetsdal G, Schulman N, Kramer W, Thaller G, Ronningen K, Lies S: A genome scan for quantitative trait loci affecting milk production in Norwegian dairy cattle. J Dairy Sci. 2002, 85 (11): 3124-3130. 10.3168/jds.S0022-0302(02)74400-7.
Viitala SM, Schulman NF, de Koning DJ, Elo K, Kinos R, Virta A, Virta J, Maki-Tanila A, Vilkki JH: Quantitative trait loci affecting milk production traits in Finnish Ayrshire dairy cattle. J Dairy Sci. 2003, 86 (5): 1828-1836. 10.3168/jds.S0022-0302(03)73769-2.
Awad A, Russ I, Emmerling R, Forster M, Medugorac I: Confirmation and refinement of a QTL on BTA5 affecting milk production traits in the Fleckvieh dual purpose cattle breed. Anim Genet. 2010, 41 (1): 1-11. 10.1111/j.1365-2052.2009.01957.x.
Lalitha S: Primer premier 5. Biotech Software & Internet Report: The Computer Software Journal for Scient. 2000, 1 (6): 270-272. 10.1089/152791600459894.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, Ehm MG: Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered. 2002, 53 (2): 79-91. 10.1159/000057986.
Wang H, Jiang L, Liu X, Yang J, Wei J, Xu J, Zhang Q, Liu JF: A Post-GWAS Replication Study Confirming the PTK2 Gene Associated with Milk Production Traits in Chinese Holstein. PLoS One. 2013, 8 (12): e83625-10.1371/journal.pone.0083625.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Rennie W, Liu C, Carmack CS, Wolenc A, Kanoria S, Lu J, Long D, Ding Y: STarMir: a web server for prediction of microRNA binding sites. Nucleic Acids Res. 2014, 42 (Web Server issue): W114-W118. 10.1093/nar/gku376.
Vasilopoulos Y, Cork MJ, Teare D, Marinou I, Ward SJ, Duff GW, Tazi-Ahnini R: A nonsynonymous substitution of cystatin A, a cysteine protease inhibitor of house dust mite protease, leads to decreased mRNA stability and shows a significant association with atopic dermatitis. Allergy. 2007, 62 (5): 514-519. 10.1111/j.1398-9995.2007.01350.x.
Sauna ZE, Kimchi-Sarfaty C: Understanding the contribution of synonymous mutations to human disease. Nat Rev Genet. 2011, 12 (10): 683-691. 10.1038/nrg3051.
Macaya D, Katsanis SH, Hefferon TW, Audlin S, Mendelsohn NJ, Roggenbuck J, Cutting GR: A synonymous mutation in TCOF1 causes Treacher Collins syndrome due to mis-splicing of a constitutive exon. Am J Med Genet A. 2009, 149A (8): 1624-1627. 10.1002/ajmg.a.32834.
Ramser J, Ahearn ME, Lenski C, Yariz KO, Hellebrand H, von Rhein M, Clark RD, Schmutzler RK, Lichtner P, Hoffman EP, Meindl A, Baumbach-Reardon L: Rare missense and synonymous variants in UBE1 are associated with X-linked infantile spinal muscular atrophy. Am J Hum Genet. 2008, 82 (1): 188-193. 10.1016/j.ajhg.2007.09.009.
Ho PA, Kuhn J, Gerbing RB, Pollard JA, Zeng R, Miller KL, Heerema NA, Raimondi SC, Hirsch BA, Franklin JL, Lange B, Gamis AS, Alonzo TA, Meshinchi S: WT1 synonymous single nucleotide polymorphism rs16754 correlates with higher mRNA expression and predicts significantly improved outcome in favorable-risk pediatric acute myeloid leukemia: a report from the children’s oncology group. J Clin Oncol. 2011, 29 (6): 704-711. 10.1200/JCO.2010.31.9327.
Capon F, Allen MH, Ameen M, Burden AD, Tillman D, Barker JN, Trembath RC: A synonymous SNP of the corneodesmosin gene leads to increased mRNA stability and demonstrates association with psoriasis across diverse ethnic groups. Hum Mol Genet. 2004, 13 (20): 2361-2368. 10.1093/hmg/ddh273.
Nackley AG, Shabalina SA, Tchivileva IE, Satterfield K, Korchynskyi O, Makarov SS, Maixner W, Diatchenko L: Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science. 2006, 314 (5807): 1930-1933. 10.1126/science.1131262.
Gao R, Reece K, Sissung T, Reed E, Price DK, Figg WD: The ERCC1 N118N polymorphism does not change cellular ERCC1 protein expression or platinum sensitivity. Mutat Res. 2011, 708 (1-2): 21-27. 10.1016/j.mrfmmm.2011.01.002.
Garcia-Roves PM, Osler ME, Holmstrom MH, Zierath JR: Gain-of-function R225Q mutation in AMP-activated protein kinase gamma3 subunit increases mitochondrial biogenesis in glycolytic skeletal muscle. J Biol Chem. 2008, 283 (51): 35724-35734. 10.1074/jbc.M805078200.
Hromas R, Davis B, Rauscher FJ, Klemsz M, Tenen D, Hoffman S, Xu D, Morris JF: Hematopoietic transcriptional regulation by the myeloid zinc finger gene, MZF-1. Curr Top Microbiol Immunol. 1996, 211: 159-164.
VanRaden PM, Wiggans GR: Derivation, calculation, and use of national animal model information. J Dairy Sci. 1991, 74 (8): 2737-2746. 10.3168/jds.S0022-0302(91)78453-1.
Jairath L, Dekkers JC, Schaeffer LR, Liu Z, Burnside EB, Kolstad B: Genetic evaluation for herd life in Canada. J Dairy Sci. 1998, 81 (2): 550-562. 10.3168/jds.S0022-0302(98)75607-3.
Thomsen H, Reinsch N, Xu N, Looft C, Grupe S, K hn C, Brockmann GA, Schwerin M, Leyhe-Horn B, Hiendleder S, Erhardt G, Medjugorac I, Russ I, Forster M, Brenig B, Reinhardt F, Reents R, Blumel J, Averdunk G, Kalm E: Comparison of estimated breeding values, daughter yield deviations and de-regressed proofs within a whole genome scan for QTL. J Anim Breed Genet. 2001, 118 (6): 357-370. 10.1046/j.1439-0388.2001.00302.x.
Martin ER, Lai EH, Gilbert JR, Rogala AR, Afshari AJ, Riley J, Finch KL, Stevens JF, Livak KJ, Slotterbeck BD, Slifer SH, Warren LL, Michael Conneally P, Schmechel DE, Purvis I, Pericak-Vance MA, Roses AD, Vance JM: SNPing away at complex diseases: analysis of single-nucleotide polymorphisms around APOE in Alzheimer disease. Am J Hum Genet. 2000, 67 (2): 383-394. 10.1086/303003.
The authors appreciate the financial support provided by the National Major Development Program of Transgenic Breeding (2014ZX0800953B), the National High Technology Research and Development Program of China (863 Program 2011AA100302 and 2013AA102503), the National Natural Science Foundations of China (31272419), the Chinese Universities Scientific Fund (2014XJ003), and Program for Changjiang Scholar and Innovation Research Team in University (IRT1191).
The authors declare that they have no competing interests.
JL designed the study and revised the manuscript. HW and LJ prepared DNA samples for SNP identification and genotyping and participated in the manuscript drafting. WW performed RNA extraction and qPCR assays. SZ, ZY and QZ participated in the data analysis and provided suggestions for the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Table S3.: Associations of SNP marker with estimated breeding values (EBVs) of milk production traits in dairy cows. (DOCX 25 KB)
About this article
Cite this article
Wang, H., Jiang, L., Wang, W. et al. Associations between variants of the HALgene and milk production traits in Chinese Holstein cows. BMC Genet 15, 125 (2014). https://doi.org/10.1186/s12863-014-0125-4
- HAL gene
- Single nucleotide polymorphism
- Milk production traits
- Chinese Holstein cows