Detection of genetic association and functional polymorphisms of UGDH affecting milk production trait in Chinese Holstein cattle
- Qing Xu†2,
- Gui Mei†1, 3,
- Dongxiao Sun1Email author,
- Qin Zhang1,
- Yuan Zhang1,
- Cengceng Yin1,
- Huiyong Chen1,
- Xiangdong Ding1 and
- Jianfeng Liu1
© Xu et al.; licensee BioMed Central Ltd. 2012
Received: 15 August 2012
Accepted: 30 October 2012
Published: 2 November 2012
We previously localized a quantitative trait locus (QTL) on bovine chromosome 6 affecting milk production traits to a 1.5-Mb region between BMS483 and MNB-209 via genome scanning followed by fine mapping.
Totally 15 genes were mapped within such linkage region through bioinformatic analysis of the cattle-human comparative map and bovine genome assembly. Of them, the UDP-glucose dehydrogenase (UGDH) was suggested as a potential positional candidate gene for milk production traits based on its corresponding physiological and biochemical functions and genetic effects. By sequencing all the coding exons and the untranslated regions in UGDH with pooled DNA of 8 sires represented the separated families detected in our previous studies, a total of ten SNPs were identified and genotyped in 1417 Holstein cows of 8 separation families. Individual SNP-based association analysis revealed 4 significant associations of SNP Ex1-1, SNP Int3-1, SNP Int5-1, and SNP Ex12-3 with milk yield (P < 0.05), and 2 significant associations of SNP Ex1-1 and SNP Ex12-3 with protein yield (P < 0.05). Furthermore, our haplotype-based association analyses indicated that haplotypes G-C-C, formed by SNP Ex12-2-SNP Int11-1-SNP Ex11-1, T-G, formed by SNP Int9-3-SNP Int9-2, and C-C, formed by SNP Int5-1-SNP Int3-1, are significantly associated with protein percentage (F=4.15; P=0.0418) and fat percentage (F=5.18~7.25; P=0.0072~0.0231). Finally, by using an in vitro expression assay, we demonstrated that the A allele of SNP Ex1-1 and T allele of SNP Ex11-1of UGDH significantly decreases the expression of UGDH by 68.0% at the RNA, and 50.1% at the protein level, suggesting that SNP Ex1-1 and Ex11-1 represent two functional polymorphisms affecting expression of UGDH and may partly contributed to the observed association of the gene with milk production traits in our samples.
Taken together, our findings strongly indicate that UGDH gene could be involved in genetic variation underlying the QTL for milk production traits.
KeywordsDairy cattle BTA6 Positional candidate Milk traits UGDH Function validation
Both theoretical and simulation studies agree that application of gene-assisted selection has the potential to increase the rate of genetic gain by pre-selecting young candidate bulls prior to progeny testing in dairy cattle. This kind of selection is based on the identification of genes that may affect the traits of interest. In dairy cattle, since the first report on quantitative trait locus (QTL) mapping by Georges et al., extensive QTL mapping has been implemented to detect QTLs with major effect on the milk production traits. So for, a total number of 1,651 QTLs for milk yield and composition traits have been reported via genome scan based on marker-QTL linkage analysis (http://www.animalgenome.org/QTLdb/cattle.html, April 29, 2012). With the completion of genome sequencing of cattle and development of comparative map of human/bovine, it is feasible to identify the genes or other functional elements, which are responsible for a quantitative trait with comparative candidate positional cloning strategy. In recent years, based on fine mapping results and positional cloning tools, several positional candidate SNPs have been confirmed to be true QTL, including the AA to GC dinucleotide substitution (K232A) in the exon8 of DGAT1 on BTA 14 and an F279 missense mutation of GHR on BTA20 for fat percentage, a mutation in the regulatory element of OPN on BTA6 and a Y581S missense mutation of ABCG2 on BTA6 for milk fat and protein concentration.
Even though QTLs for milk production traits have been found in almost 29 chromosomes, BTA6 is one of the most concerned chromosomes in QTL mapping for milk yield and content. Of the 1651 reported QTLs for milk production traits, 193 are on BTA6. In our previous study in a daughter-design Chinese Holstein population consisting of 26 sire families, a QTL for milk yield, fat yield, protein yield and fat percentage was detected around the microsatellite marker BM470 with a confidence interval of 4 cM on BTA6 by using both linear regression and variance component approaches[8, 9]. This region was also found to harbor one or more QTLs for milk production traits by several previous independent studies[10–13]. Thus, further fine mapping has been employed with 15 microsatellites in the region of 14.3 cM between markers BMS690 and BM4528 in 918 daughters of 8 segregating families identified by Chen et al.[8, 9]. With haplotype sharing based LD and single marker regression mapping, a QTL with significant effects on milk yield and milk composition was determined between BMS483 and MNB209, which spans a narrow genetic distance of 0.6 cM and a physical distance of 1.5 Mb. In the present study, by using a high-resolution whole-genome cattle-human comparative map[15–17] and the bovine genome assembly Btau 4.0 (http://www.ncbi.nlm.nih.gov), a total of 15 known genes have been mapped within such narrow linkage region. Out of them, we herein focused on UGDH as the most plausible comparative functional candidate gene affecting milk production traits. The enzyme UGDH can convert UDP-glucose to UDP-glucuronic acid, a critical component of the glycosaminoglycans, hyaluronan, chondroitin sulfate, and heparan sulfate, thus promotes normal cellular growth, embryogenesis and adult organism physiology[18, 19]. UGDH is also implicated in the progression of epithelial cancers such as colon, breast, and prostate[20–22]. The purpose of the study was to determine the genetic effect of the UGDH gene on milk yield and component traits in dairy cows and regulation by functional polymorphisms.
The blood samples and frozen semen were collected along with the regular quarantine inspection of the farms and breeding station, so no ethical approval was required for this study.
Animals and phenotypic data
A total of 1417 Chinese Holstein cows were selected from 8 sire families with 67–540 daughters in each family. Such 8 sires represent all the segregating families detected in our previous QTL mapping studies[8, 9, 14]. Estimated breeding values (EBV) for five milk production traits (i.e., milk yield, fat yield, protein yield, fat percentage, and protein percentage over 305 days) were provided by the Dairy Data Processing Center, Dairy Association of China (DAC) which was calculated with a random regression multiple traits test-day model based on 6,980,000 test-day records of 585,121 Holstein cows collected from 1993 to 2008 in China.
Genomic DNA was extracted from whole blood samples of cows and frozen semen of the 8 bulls by a standard phenol-chloroform method and stored at −20°C.
Positional candidate cloning, SNP discovery and genotyping
Within the 1.5-Mb region between microsatellites BMS483 and MNB-209 on BTA 6 where harbored a QTL for milk production traits[8, 9, 14], positional candidate cloning was implemented to mine the known functional genes by bioinformatic analysis of a high-resolution whole-genome cattle-human comparative map[15–17] and the bovine genome assembly Btau 4.0 (http://www.ncbi.nlm.nih.gov). Then, based on Btau 4.0, full-length cDNA sequences of candidate genes were obtained and totally 88 pairs of primers were designed with Primer3.0 program to amplify the entire coding region of each candidate gene to identify potential polymorphisms (primers were not shown with the exception of UGDH).
PCR primers for amplification and identification of single nucleotide polymorphisms in UGDH
Forward primer sequence
Reverse primer sequence
Amplicon size (bp)
GAA TGA ATC CTC CGC CT
TTCCTG TGG TGA TGT TCG
GGG CTA CTG ATG ATG TT
CTC TGT CTG GGT TCT TC
AGC AGA AAG TAT TCG TCG
CAG CTT AGC AGC AGG TaG
ACC AAC CCA CCC ATT GCC T
CGG ACA CGA CTG AGC TAC T
AGT TGA TTA GTT TCT TGC
TAA TAA TTG AGT TCC TGG
TTT GGA GAA GTT GGC TAA GG
TGC TAT GCT TCT GAA TGG C
CGT TCA TAT CCT CAC AAG A
TGG GGA AAT CGT TAC AGT A
TTC AAC CAG ACC ATT CCA CT
TGC GTT CAT AAT CCA GTT CC
TAA GGT GAA GTG ACG GAA GC
GGA GCA ACT AAG CCT GTG TG
GGG AATGTA GAC TGA CTC
GCT ACT CAT CCA ATC ACT
Single SNP analysis
Associations between individual SNP and the five milk production traits were estimated using the MIXED procedure of the SAS 8.02 software based on the model y = μ + G + a + e, where y is the EBV of cows, μ is the general mean, G is the SNP genotype effect, a is the random polygenic effect with distribution of N(0, A [σ a 2]) (A is the additive genetic relationship matrix among all individuals, σ a 2 is the additive genetic variance), and e is a random residual with distribution of N(0, W σ a 2) (W is a diagonal matrix with the diagonal elements equal to 1/REL ij , REL ij is the reliability of the EBV of daughter j in family i, and σ a 2 is the residual error variance).
Pair-wise linkage disequilibrium (LD) between all SNP pairs and haplotype blocks were estimated using the program Haploview4.0. Then, haplotype reconstruction of each cow was carried out within the blocks via expectation maximization (EM) algorithm. Association between each haplotype and the five milk production traits was performed based on the model y ij = μ + s i + bh ij + e ij , where s i is the fixed effect of sire i, b is the regression coefficient, hij is an indicator variable with a value 0, 1 or 2 to indicate the copy number of the haplotype carried by the individual.
Both analyses (e.g. Single SNP analysis and Haplotype analysis) were corrected for multiple testing using false discroverage rate (FDR) methods. We declared a significant SNP or haplotype if the corresponding FDR value < 0.05.
cDNA synthesis, plasmid construction and mutagenesis
On the basis of association analysis results and the nucleotide location of each SNP within the gene, we investigated three SNPs with two located in exon 1 (SNPs exon1-1 and exon1-2), one in exon 11 (SNP exon11-1). The complete coding region of UGDH was synthesized based on the sequence of cattle full-length cDNA at NCBI (GenBank Accession No. NM: 174211). To make cloning step easier, we added a restriction endonuclease enzyme digestion site to the 5′ and 3′end of the synthesized gene, respectively. The synthesized DNA fragment was double-digested with restriction enzymes Bam H I and Xho I and then cloned into the digested pcDNA3.1/myc-HisA expression vector (Clontech, Mountain View, CA) with the same enzymes; this yielded the first plasmid construct, named as pcDNA3-GAC (g-a-c). The other seven plasmid constructs (i.e., a-a-c, g-c-c, g-a-t, a-c-c, g-c-t, a-a-t, and a-c-t) were generated from the g-a-c plasmid using a site-directed mutagenesis kit purchased from Stratagene Inc. (La Jolla, CA). All eight plasmid constructs were verified by automated DNA sequencing in both directions.
Cell culture and transfection
COS-7 cells (an African Green Monkey SV40-transf'd kidney fibroblast cell line) were purchased from the American Tissue Cell Culture Inc. (Manassas, VA), and cultured in RPMI-1640 supplemented with 10% fetal bovine serum (Invitrogen, Carlsbad, CA) at 37°C in 5% CO2. The culture medium was changed every other day. The cells were transfected by Lipofectamine 2000 (Invitrogen, Carlsbad, CA) according to the manufacturer’s protocol. Briefly, a day prior to transfection, the cells were trypsinized and plated at 4 × 104 cells per well on a 48-well plate. Subsequently, 0.4 μg of plasmid DNA containing the appropriate structure was diluted in 20 μl of Opti-MEM I (Invitrogen, Carlsbad, CA) Reduced Serum Medium and mixed with an equal volume of Opti-MEM I Reduced Serum Medium containing 1 μl of Lipofectamine 2000. The mixture was incubated for 20 min at room temperature and added directly to the cells in a 48-well plate. The cells were cultured for an additional 48 hours prior to harvesting for extraction of RNA and protein for real-time RT-PCR and Western blotting, respectively.
RNA isolation and real-time quantitative RT-PCR
To evaluate the regulatory effect of SNPs exon1-1 (G/A), exon1-2 (A/C) and exon11-1 (C/T) on the expression of UGDH, eight plasmids were constructed corresponding to different allele combinations of exon1-1 (G/A), exon1-2 (A/C) or exon11-1 (C/T). Total RNA from cell transfected by eight plasmids respectively and control cells was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA). One microgram of total RNA was reverse transcribed in a final volume of 20 ul containing 4 ul of 5 first-strand buffer (250 mM Tris–HCl, pH 8.3; 375 mM KCl; 15 mM MgCl2), 10 mM DTT, 0.5 mM each dNTP, 40 U RNaseOUT™, 1 ul of 50 nM random hexamers, and 200 U of Superscript II RNase H- reverse transcriptase (Invitrogen). The TaqMan® specific probes and primers used for real-time PCR were designed according to the above synthesized UGDH fragment and myc sequences in pcDNA3.1/myc-HisA expression vector were: forward primer 5′-CTCCTTCTGGTGAAATTCCAAAGT-3′, reverse primer 5′-TCTAGACTCGAGTACTCTGGGTTTCTT-3′, and the probe 5′-AGTCTTCAGGATATGTGCCCAA-3′. Amplification of 2 ul of cDNA was carried out in a total volume of 20 ul according to the manual of TaqMan® Gene Expression Assays (Applied Biosystems). The mRNA level of each gene was determined using a calibration method and normalized to that of 18S rRNA in each sample. In addition, the RNA of transfected cells was treated with RNA-free DNase I (Ambion Inc, Foster City, CA) to remove any potential remaining plasmid DNA from the RNA sample prior to reverse transcription. All the primers and probe were synthesized by Applied Biosystems Inc. (Foster, CA).
Total protein was extracted from individual transfected cell from an independent experiment by homogenization in RIPA buffer (50 mM Tris–HCl, pH 8.0; 150 mM NaCl, 1% NP-40, 0.5% sodium deoxycholate, and 0.1% SDS), and the protein concentration was determined by the Bio-Rad Protein Assay (Bio-Rad Laboratories, Hercules, CA). Ten micrograms of total protein was separated by 10% SDS-PAGE followed by transfer to nitrocellulose membranes (0.45 um) at 25V overnight at 4°C. The membrane was first incubated in a blocking buffer (5% nonfat milk and 0.2% Tween 20) for 1.5h at room temperature and then 1.5h at room temperature in the blocking buffer containing cattle anti-His antibody (dilution 1:5000; Abcam Inc, Cambridge, MA). After three washes in TBST (10 mM Tris–HCl, pH 8.0; 0.15 M NaCl; 0.2% Tween 20) for 10 min each, the membranes were exposed to horseradish peroxidase-conjugated secondary antibody at room temperature for 1.5h, and then exposed to X-ray film. After hybridization with the antibody of interest, membranes were stripped and re-probed with antibody to β-Actin (dilution 1:5000; Abcam Inc, Cambridge, MA), which was used for normalization of the protein content of each sample. Then, the films were scanned for quantitative analysis with ImageQuant 5.2 (Molecular Dynamics, Sunnyvale, CA).
Positional candidate cloning
With bioinformatic analysis of whole-genome cattle-human comparative map and the bovine genome assembly Btau 4.0, totally 13 annotated genes and 2 putative genes were mapped within the previously identified 1.5Mb linkage region on BTA6 (Mei et al., 2009), which were compose of TBC1D1, KLF3, TLR10, TLR1, TLR6, LOC511583, TMEM156, LOC528668, WDR19, LOC514842, LOC514925, RPL9, LIAS and UGDH. Of them, five genes including KLF3, KLHL5, RFC1, WDR19 and UGDH, were suggested as candidates for milk production traits based on their corresponding physiological and biochemical functions in human, so that the entire coding region of each gene was screened to identify any potential polymorphisms by direct sequencing of pooled DNA from 8 segregating sire families. As a result, a total of 64 SNPs were revealed and their associations with five milk yield and milk composition traits were determined (data of KLF3, KLHL5, RFC1 and WDR19 were not shown). We herein considered UGDH as the most plausible comparative functional candidate gene affecting milk production traits.
SNPs identification and selection
Information for 10 SNPs in UGDH and allele frequencies in Chinese Holstein population
Position (BTA 6)
Primer and probe sequence
CAAACAAATG [A/C] AACAAAACCA
AAATCACTCC [C/T] CTCTTACTCC
TCAGACTCTC [G/A] ATTCCTCCTG
GAGAAATAAG [T/C] CCTTGGCTTA
AGCCATAGCA [T/G] TACTTAGTAT
TGGTTGCAGC [A/C] CTCGGTGCCT
GAGCCGAGGA [G/A] ATAGAAGTGG
Association between individual SNP and milk production traits
Significant UGDH SNPs associated with milk production traits in Chinese Holstein population (P<0.05)
Milk yield, kg LSMEAN ± SE
Protein yield, kg LSMEAN ± SE
Associations between haplotypes and milk production traits
Associations of UGDH haplotypes with milk production traits in Chinese Holstein population
Analysis of UGDH and its mutant Ex1-1 (G/A), Ex1-2 (A/C) and Ex11-1 (C/T).
After normalization to the corresponding 18S rRNA of each sample, we found all mutant constructs except g-c-c showed significantly lower mRNA expression than the construct g-a-c (~52.0%; P < 0.05). Moreover, three mutant constructs (i.e., a-a-c, a-c-c, and a-c-t) were observed to be highly significant (P < 0.01), where UGDH mRNA was down-regulated by 63%, 62%, and 68%, respectively (Figure3B). At the protein level, as shown in Figure3D, the seven mutant constructs showed generally lower expression than the construct g-a-c (−15.8~−51.1%) and among them, constructs a-a-t, a-c-c, a-a-c, g-c-c, and a-c-t reached the significance level of 0.05, with their UGDH protein concentration being reduced by 37.8%,51.1%, 34.7%, 26.5%, and 50.1%, respectively.
These results indicate that at mRNA level, the mutation of Ex1-2 has little impact on the expression of UGDH and that all differences observed among the eight constructs were attributable primarily to SNPs Ex1-1 and Ex11-1, for which the A allele of SNP Ex1-1 and T allele of SNP Ex11-1 caused lower expression than the G and C allele respectively. But at protein level, the two alleles of SNP Ex1-1 play more important role in the expression of UGDH than those of SNPs Ex1-2 and Ex11-1, where the G allele of SNP Ex1-1 corresponds to higher protein concentration.
In this study, we provide a strong evidence for the significant associations of UGDH, with 2 milk production traits in a Chinese Holstein cattle population. Moreover, we demonstrated that the SNPs Ex1-1 and Ex11-1 in exonic regions of UGDH are two functional polymorphisms in which A allele of SNP Ex1-1 and T allele of SNP Ex11-1 cause lower expression level of UGDH compared with the G and C allele respectively. Together, our findings not only confirmed our early finding that a QTL located in the 1.5-Mb region between BMS483 and MNB-209 for milk yield and protein yield, but also strongly suggest that the SNPs Ex1-1 and Ex11-1 in UGDH might be QTN responsible for this QTL.
Previously, we mapped a QTL near BMS470 with effects on milk production traits with daughter design using 14 microsatellites on bovine chromosome 6 in a Chinese Holstein population involving 26 paternal half-sib families with 2356 daughters[8, 9]. By increasing the marker density with 17 microsatellite markers spanning from BMS690 to BMS4528 near BMS470 and genotyping them for 918 daughters from 8 sires families which had been proved to contain segregating QTLs affecting milk production traits in the previous studies[8, 9], we further fine mapped this QTL in region of BMS483 and MNB-209 by employing linkage disequilibrium (LD) mapping approaches. On the base of these studies, the UGDH gene within this interval was selected as a potential positional candidate gene for further study by genetic association and functional analysis as reported in this study. Association analysis of individual SNPs indicated that the SNP Ex1-1 located in the 5` UTR (c.exon1-1 G>A) was significantly associated, even after multiple testing correction of FDR, with milk yield (P = 0.003). Moreover, haplotype analysis revealed two significant associations of two percentage traits with haplotypes formed by SNP Ex12-2- SNP Int11-1- SNP Ex11-1 and SNP Int5-1- SNP Int3-1, respectively. Given the aforementioned association results of SNP Ex1-1 and SNP Ex11-1 with milk production traits and their locations within the gene, we hypothesized that these two SNPs and the other SNP Ex1-2 (located in the 5` UTR of UGDH and near to the SNP Ex1-1) contribute to the observed association of UGDH with milk production traits. By expressing plasmid constructs containing different alleles for above three SNPs in COS 7 cells, we found that the plasmid construct containing the A allele of SNP Ex1-1 and T allele of SNP Ex11-1 produced significantly lower expression of mRNA than the G and C allele respectively. At protein level, the two alleles of SNP Ex1-1 play more important role on the expression of UGDH than those of SNPs Ex1-2 and SNP Ex11-1, where the G allele of SNP Ex1-1 corresponds to higher protein concentration. Thus, we concluded that SNP Ex1-1 is involved in the regulation of the expression both at the transcription and protein level and SNP Ex11-1 is involved in the regulation of the expression of UGDH at the transcription level Although further study is required to determine how these two DNA polymorphisms modulates the expression of UGDH, they are most likely responsible, at least part, for the observed effect of above QTL affecting milk yield and composition in dairy cattle used in current study. As for SNP Ex1-2, in agreement with the results of the individual association analysis on this SNP, we did find no significant difference between the two alleles at either the RNA or the protein level. Although a non-synonymous SNP in a gene can change its function(s), causing the gene to be associated with a disease, some recent studies found that synonymous SNPs without altered coding sequences also can affect the gene’s function. These SNPs may alter the structure of the substrate and inhibitor interaction site by a rare codon in their polymorphism or affect gene expression by altering mRNA secondary structure.
Based on the previously identified 1.5-Mb region between markers BMS483 and MNB-209 on BTA6 including a QTL for milk production traits, the UGDH gene was identified and shown to be associated milk yield and milk composition. Further in vitro expression assay demonstrated that the SNPs in exon1 and exon11 of UGDH represent two functional polymorphisms affecting expression of UGDH and might partly contributed to the observed association with milk production traits.
Bos Taurus Automosome
Dairy association of China
AcylCoA: diacylglycerol acyltransferase
Growth hormone receptor
Single nucleotide polymorphism
Quantitative trait locus
This work was supported by the National Natural Science Foundation of China (31072016), Beijing InnovationTeam of Technology System in the National Dairy Industry, Key Development of New Transgenic Breeds Program of China (2009ZX08009-156B), National “948” Project (2011-G2A), National Key Technologies R & D Program (2011BAD28B02, 10ZCKFNC02700) and Agricultural Science Technology Achievement Transformation Project (20110409002). The authors would also like to thank the Dairy Association of China for providing semen samples and pedigree information of bull.
- Khatkar MS, Thomson PC, Tammen I, Raadsma HW: Quantitative trait loci mapping in dairy cattle: review and meta-analysis. Genet Sel Evol. 2004, 36: 163-190. 10.1186/1297-9686-36-2-163.PubMed CentralView ArticlePubMed
- Georges M, Nielsen D, Mackinnon M, Mishra A, Okimoto R, Pasquino AT, Sargeant LS, Sorensen A, Steele MR, Zhao X: Mapping quantitative trait loci controlling milk production in dairy cattle by exploiting progeny testing. Genetics. 1995, 139: 907-920.PubMed CentralPubMed
- Watanabe M, Nikaido M, Tsuda TT, Kobayashi T, Mindell D, Cao Y, Okada N, Hasegawa M: New candidate species most closely related to penguins. Gene. 2006, 378: 65-73.View ArticlePubMed
- Grisart B, Farnir F, Karim L, Cambisano N, Kim JJ, Kvasz A, Mni M, Simon P, Frère JM, Coppieters W, Georges M: Genetic and functional confirmation of the causality of the DGAT1 K232A quantitative trait nucleotide in affecting milk yield and composition. Proc Natl Acad Sci USA. 2004, 101: 2398-2403. 10.1073/pnas.0308518100.PubMed CentralView ArticlePubMed
- 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: 253-266.PubMed CentralPubMed
- Schnabel RD, Kim JJ, Ashwell MS, Sonstegard TS, Van Tassell CP, Connor EE, Taylor JF: Fine-mapping milk production quantitative trait loci on BTA6: analysis of the bovine osteopontin gene. Proc Natl Acad Sci USA. 2005, 102: 6896-6901. 10.1073/pnas.0502398102.PubMed CentralView ArticlePubMed
- Cohen-Zinder M, Seroussi E, Larkin DM, Loor JJ, Everts-van der Wind A, Lee JH, Drackley JK, Band MR, Hernandez AG, Shani M, Lewin HA, Weller JI, Ron M: Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle. Genome Res. 2005, 15: 936-944. 10.1101/gr.3806705.PubMed CentralView ArticlePubMed
- Chen HY, Zhang Q, Wang CK, Shu J, Mei G, Yin CC, Hu F, Xu JJ, Gong WJ, Li HJ, Qiu XT: Mapping QTLs on BTA6 affecting milk production traits in a Chinese Holstein population. Chin Sci Bull. 2005, 50: 1737-1742. 10.1360/982005-725.View Article
- 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: 782-790. 10.3168/jds.S0022-0302(06)72140-3.View ArticlePubMed
- Zhang Q, Boichard D, Hoeschele I, Ernst C, Eggen A, murkve B, Pfister-Genskow M, Witte LA, Grignola FE, Uimari P, Thaller G, Bishop MD: Mapping quantitative trait loci for milk production and health of dairy cattle in a large outbred pedigree. Genetics. 1998, 149: 1959-1973.PubMed CentralPubMed
- Velmala RJ, Vilkki HJ, Elo KT, de Koning DJ, Maki-Tanila AV: A search for quantitative trait loci for milk production traits on chromosome 6 in Finnish Ayrshire cattle. Anim Genet. 1999, 30: 136-143. 10.1046/j.1365-2052.1999.00435.x.View ArticlePubMed
- Szyda J, Liu Z, Reinhardt F, Reents R: Estimation of quantitative trait loci parameters for milk production traits in German Holstein dairy cattle population. J Dairy Sci. 2005, 88: 356-367. 10.3168/jds.S0022-0302(05)72695-3.View ArticlePubMed
- Harder B, Bennewitz J, Reinsch N, Thaller G, Thomsen H, Kuhn C, Schwerin M, Erhardt G, Förster M, Reinhardt F, Kalm E: Mapping of quantitative trait loci for lactation persistency traits in German Holstein dairy cattle. J Anim Breed Genet. 2006, 123: 89-96. 10.1111/j.1439-0388.2006.00577.x.View ArticlePubMed
- Mei G, Yin C, Ding X, Zhang Q: Fine mapping quantitative trait loci affecting milk production traits on bovine chromosome 6 in a Chinese Holstein population. J Genet Genomics. 2009, 36: 653-660. 10.1016/S1673-8527(08)60157-5.View ArticlePubMed
- Weikard R, Goldammer T, Laurent P, Womack JE, Kuehn C: A gene-based high-resolution comparative radiation hybrid map as a framework for genome sequence assembly of a bovine chromosome 6 region associated with QTL for growth, body composition, and milk performance traits. BMC Genomics. 2006, 7: 53-10.1186/1471-2164-7-53.PubMed CentralView ArticlePubMed
- Everts-van der Wind A, Kata SR, Band MR, Rebeiz M, Larkin DM, Everts RE, Green CA, Liu L, Natarajan S, Goldammer T, Lee JH, McKay S, Womack JE, Lewin HA: A 1463 gene cattle-human comparative map with anchor points defined by human genome sequence coordinates. Genome Res. 2004, 14: 1424-1437. 10.1101/gr.2554404.PubMed CentralView ArticlePubMed
- Everts-van der Wind A, Larkin DM, Green CA, Elliott JS, Olmstead CA, Chiu R, Schein JE, Marra MA, Womack JE, Lewin HA: A high-resolution whole-genome cattle-human comparative map reveals details of mammalian chromosome evolution. Proc Natl Acad Sci USA. 2005, 102: 18526-18531. 10.1073/pnas.0509285102.View ArticlePubMed
- Trochon V, Mabilat C, Bertrand P, Legrand Y, Smadja-Joffe F, Soria C, Delpech B, Lu H: Evidence of involvement of CD44 in endothelial cell proliferation, migration and angiogenesis in vitro. Int J Cancer. 1996, 66: C664-668. 10.1002/(SICI)1097-0215(19960529)66:5<664::AID-IJC14>3.0.CO;2-4.View Article
- Iozzo RV, Murdoch AD: Proteoglycans of the extracellular environment: clues from the gene and protein side offer novel perspectives in molecular diversity and function. FASEB J. 1996, 10: 598-614.PubMed
- Ropponen K, Tammi M, Parkkinen J, Eskelinen M, Tammi R, Lipponen P, Agren U, Alhava E, Kosma VM: Tumor cell-associated hyaluronan as an unfavorable prognostic factor in colorectal cancer. Cancer Res. 1998, 58: 342-347.PubMed
- Auvinen P, Tammi R, Parkkinen J, Tammi M, Agren U, Johansson R, Hirvikoski P, Eskelinen M, Kosma VM: Hyaluronan in peritumoral stroma and malignant cells associates with breast cancer spreading and predicts survival. Am J Pathol. 2000, 156: 529-536. 10.1016/S0002-9440(10)64757-8.PubMed CentralView ArticlePubMed
- Aaltomaa S, Lipponen P, Tammi R, Tammi M, Viitanen J, Kankkunen JP, Kosma VM: Strong Stromal Hyaluronan Expression Is Associated with PSA Recurrence in Local Prostate Cancer. Urol Int. 2002, 69: 266-272. 10.1159/000066123.View ArticlePubMed
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMed
- Ding X, Zhang Q, Flury C, Simianer H: Haplotype reconstruction and estimation of haplotype frequencies from nuclear families with only one parent available. Hum Hered. 2006, 62: 12-19. 10.1159/000095598.View ArticlePubMed
- Winer J, Jung CK, Shackel I, Williams PM: Development and validation of real-time quantitative reverse transcriptase-polymerase chain reaction for monitoring gene expression in cardiac myocytes in vitro. Anal Biochem. 1999, 270: 41-9. 10.1006/abio.1999.4085.View ArticlePubMed
- Kanda A, Chen W, Othman M, Branham KE, Brooks M, Khanna R, He S, Lyons R, Abecasis GR, Swaroop A: A variant of mitochondrial protein LOC387715/ARMS2, not HTRA1, is strongly associated with age-related macular degeneration. Proc Natl Acad Sci USA. 2007, 104: 16227-16232. 10.1073/pnas.0703933104.PubMed CentralView ArticlePubMed
- Kimchi-Sarfaty C, Oh JM, Kim IW, Sauna ZE, Calcagno AM, Ambudkar SV, Gottesman MM: A "silent" polymorphism in the MDR1 gene changes substrate specificity. Science. 2007, 315: 525-528. 10.1126/science.1135308.View ArticlePubMed
- 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: 1930-1933. 10.1126/science.1131262.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.