|
|
||||||||
Department of Dairy Science, University of Wisconsin, Madison 53706
1 Corresponding author: hkhatib{at}wisc.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: uterine milk protein gene differential allelic expression candidate gene productive life
| INTRODUCTION |
|---|
|
|
|---|
The UTMP gene was chosen for this study because previous studies of the UTMP region have shown an association with productive life (PL) and milk composition traits in dairy cattle. Heyen et al. (1999) reported an association between marker ILSTS054 and PL in the North American Holstein population. Khatib et al. (2005) reported a significant association between the protease inhibitor gene (PI) and PL and milk composition traits in Holstein dairy cattle. In addition, both PI and UTMP have protease-inhibition functions, and the genes are both located on chromosome 21, within a 321.6-kb genomic region. Thus, possible roles of UTMP in health traits and its proximity to PI were the motives for investigation of possible effects of UTMP on health traits in dairy cattle.
Recent studies in humans (Yan et al., 2002; Bray et al., 2003; Lo et al., 2003; Pastinen et al., 2003), in mice (Cowles et al., 2002), and in corn (Guo et al., 2004) have shown that alleles of nonimprinted genes are not expressed equally at the mRNA level, and such differential allelic expression (DAE) may be the basis for variation in disease susceptibility and in determining phenotypic diversity. Yan et al. (2002) examined single nucleotide polymorphisms (SNP) for 13 genes in 96 individuals and found significant differences in allelic variation in 6 out of the 13 genes studied, with 1.3- to 4.3-fold differences between alleles. In a survey of allelic variation among human genes, Pastinen et al. (2003) selected 129 genes based on their potential function in the immune system and metabolic disorders and identified 18% of the genes with deviations from the equimolar ratio between the 2 alleles. Lo et al. (2003) examined allele-specific gene expression of 1,063 transcribed SNP by using the Affymetrix HuSNP chip system. The authors found that among the 602 genes that were heterozygous, 326 showed preferential expression of 1 allele in at least 1 individual, and 170 of those genes showed a greater than 4-fold difference between the 2 alleles.
Accumulating evidence indicates that DAE is associated with phenotypic variability of quantitative and qualitative traits. For example, Hirota and colleagues (2004) studied the relationship between DAE of the human cytochrome P450 3A4 (CYP3A4) gene and metabolic activity in the liver. Individuals with a large difference in expression levels between the 2 alleles showed reduced total CYP3A4 mRNA and, hence, the total metabolic activity of the gene was reduced. That was the first study to demonstrate that CYP3A4 expression levels were affected by allele-specific expression levels.
In this study, we analyzed the expression pattern of UTMP in fetal and adult cattle tissues and investigated its association with production traits in 2 independent Holstein populations. In addition, we examined the relationship between DAE of UTMP and association with quantitative traits.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Reverse Transcription PCR
Reverse transcription (RT) of the bovine UTMP gene was performed using the primers UMP1 (ATGTCCCA CGGGAGAATGAA) and UMP2 (CCTCCTTCTGGGTG AGTTTG). The temperature cycles were as follows: 50°C for 30 min (reverse transcription); 95°C for 15 min, followed by 33 cycles of 94°C for 1 min, 58°C for 1 min, 72°C for 1 min; and a final extension at 72°C for 10 min. Primers ß-actinF (CAGCACAATGAAGAT CAAGATCATC) and ß-actinR (AAAGGGTGTAACGC AGCTAACAGT) were used to amplify 191 bp from the housekeeping gene ß-actin (GenBank accession number NM_173979) cDNA. The RT-PCR products were electrophoresed on a 1.5% agarose gel. The detectable amplification of any product was interpreted as a positive indication of UTMP expression; no formal quantitative test was performed. Primers ß-actinF and ß-actinR were used to amplify the housekeeping gene ß-actin as a control for the expression of UTMP in the various tissues.
Phenotypic Data
The association between UTMP variants and production traits was studied in 2 independent Holstein populations: the Cooperative Dairy DNA Repository (CDDR) population of AI sires and cows from the University of Wisconsin resource population (henceforth, UW resource population). The UW population was created to search (using a daughter design) for genetic markers in association with susceptibility to paratuberculosis. The 12 sire families of this population were chosen from a large number of candidate bulls with large numbers of daughters in production in 2000. Criteria for the final selection of the 12 bulls included having large numbers of daughters in production and relatively low pedigree relationships among the chosen bulls to more broadly sample the chromosomes of the US Holstein population. Phenotypic data for both analyses were obtained from the Animal Improvement Programs Laboratory (Beltsville, MD). The CDDR data were PTA for milk yield, milk protein and fat yields, milk protein and fat percentages, PL in months, and SCS from August 2005. The phenotypic data from the UW resource population comprised yield deviations for milk yield, milk protein and fat yields, and PL. Productive life is a longevity trait defined as a cows total lifetime months in milk with limits of 10 mo per lactation and 7 yr of age (VanRaden and Wiggans, 1995). Deviated records of PL were based on direct observations and also correlated traits measured earlier in life. The average reliability of PTA for PL in the CDDR population was 71%.
Polymorphism Detection and Genotyping
Semen samples from 28 Holstein sires and their 1,362 sons (23 to 104 sons per sire) were obtained from the CDDR, which is maintained by the USDA Bovine Functional Genomics Laboratory. In addition, blood samples were obtained from 913 cows from the UW resource population. For a detailed description of this population, see Cobanoglu et al. (2006). Genomic DNA was extracted from semen and blood samples using standard procedures. To detect SNP in the UTMP gene (GenBank accession number L22095), DNA pools were constructed from 220 bovine samples and amplified with the primers UTMP3 (GGCCCTACATCAAGCTG AGA) and UTMP4 (CTACTCAACTTGGGGGTTGA) as previously described (Leonard et al., 2005). The PCR products of the pooled DNA samples were sequenced and SNP were identified by visually inspecting sequence traces. For individual genotyping, primers UTMP3 and UTMP4 were used to amplify a 327-bp fragment of UTMP from genomic DNA. The PCR products were subjected to RFLP detection using the restriction enzyme BsrI that distinguishes alleles A and G of the SNP at position 1296. The digestion products were electrophoresed on a 3.0% agarose gel.
Allele Expression Quantification
Differential allelic expression of the UTMP gene was quantified by a sequencing-based approach and by single-nucleotide primer extension analysis. In the sequencing approach, RT-PCR products amplified from heterozygous individuals were sequenced according to standard procedures. Sequencing reactions of PCR amplicons were performed for 50 cycles at 96°C for 10 s, 58°C for 4 min, and a final extension at 72°C for 7 min. Data were analyzed using Applied Biosystems Sequencing Analysis (version 5.0; Applied Biosystems, Foster City, CA). The SNP were identified by visually inspecting each base in sequencing traces. Allelic variation was estimated by measuring the proportions of the peak heights of the 2 alternative alleles of the SNP. In the primer extension approach (Norton et al., 2002), RT-PCR products were purified from agarose gel using the GFX PCR DNA purification Kit (Amersham Biosciences, Piscataway, NJ). Primer extension reactions were prepared in a total volume of 10 µL containing 1 µL of purified RT-PCR product, 5 µL of SnaPshot Kit (Applied Biosystems), 0.02 µM of extension primer, and 1 µL of deionized water. The primer extension reactions were subjected to 25 cycles of 96°C for 10 s, 50°C for 5 s, and 60°C for 30 s. In a postextension treatment, reactions were treated with 1 unit of shrimp alkaline phosphatase at 37°C for 1 h followed by deactivation of the enzyme at 75°C for 15 min. The primer 1179.ext (GTGAAGGCAAAGGAGGTCCC) was used in the primer extension reactions. Samples were electrophoresed on a 3700 DNA sequencer (PE Applied Biosystems) and data were analyzed by using Genescan Analyzer version 2.5 software (PE Applied Biosystems). Allelic ratios were estimated by measuring the proportions of the peak heights of the 2 fluorescently labeled variants of the SNP found at the same position. No significance tests were performed.
Statistical Analysis
For the CDDR population, data relative to each trait were analyzed using an allele substitution model expressed as
![]() |
where yij is the PTA value relative to son j of sire i, µ is a general constant (intercept), si is the fixed effect of sire i, ß is the regression coefficient representing half of the allele substitution effect (
/2), xij is the number of G alleles (0, 1, or 2) on son j of sire i, and
ij is a residual term. Reliabilities of the sons PTA were incorporated as weights in the model to obtain weighted least squares estimates of the allele substitution effects.
Data relative to the daughters in the UW resource population were analyzed using 2 models. First, an allele substitution model similar to the CDDR described above was considered, which is expressed as
![]() |
where yijk represents the yield deviation of daughter j of sire i,
is an effect associated with Mycobacterium paratuberculosis infection status, and pij is an indicator variable assuming the values 0 and 1 for noninfected and infected cows, respectively. The index k indicates the genotype of the cows (k = AA, AG, or GG), and xk is the number of G alleles in the genotype k (xk = 0, 1, or 2, respectively). The remaining terms in the model are as previously defined. Mycobacterium paratuberculosis infection status was included in the model because the UW population was originally created to search for genetic markers associated with susceptibility to paratuberculosis.
In addition, a second model was used for each trait, which is described as
![]() |
where gk is the effect of the UTMP genotype k. Additive genetic effects were estimated as half the difference between the 2 homozygous groups, and dominance effects were estimated as the difference between the average of the 2 homozygous groups and the heterozygous group. The AA genotype was set as baseline for estimating the genotypic effects. All the analyses were implemented using the GLM procedure of SAS (SAS Institute, 1999).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
This is the first report on expression analysis of UTMP transcripts in a wide range of cattle tissues. Expression of UTMP was examined in a total of 198 tissues obtained from 10 fetuses and 17 pregnant cows. Transcripts of UTMP were found in all cotyledon tissues examined and to a lesser extent in ovary, pituitary, and spleen fetal tissues. The UTMP transcripts were not found in any other fetal tissues examined. To exclude the possibility of RT-PCR reaction failure in tissues where UTMP expression was not observed, the ß-actin gene was amplified for all tissues. Expression of the ßactin gene was detected in all examined tissues. However, we cannot exclude the possibility of variable expression of ß-actin among the different tissues or individuals examined. Although the expression of UTMP in fetal tissues was surprising, its presence in reproductive tissues (cotyledon and ovary) might shed more light on the role of this gene in maintenance of pregnancy and fetal survival. Roberts and Bazer (1988) suggested that nutrition of the conceptus, growth control, and immunosuppression of the local maternal system are possible functions of the UTMP.
The expression of UTMP has been mainly studied in sheep, with only a few studies on expression and localization of the bovine gene. In sheep, UTMP is the predominant protein present in the uterine fluid from d 30 of pregnancy until term (Moffatt et al., 1987); UTMP mRNA is not detected in endometrium until d 14 of pregnancy (Ing and Roberts, 1989). In addition, mRNA levels of ovine UTMP have been found to be affected by day of pregnancy (Stewart et al., 2000).
In this study, bovine UTMP was predominantly expressed in endometrium, ovary, and caruncle tissues in pregnant cows. Inconsistent amplification of RT-PCR products was observed in oocytes, lung, liver, and spleen tissues. On the other hand, heart, muscle, and pancreas tissues did not express UTMP in any examined individual. Thus, the predominant expression of UTMP in reproductive tissues is consistent with an important role for this protein in reproductive success as previously suggested for the ovine gene (Ing and Roberts, 1989). However, in contrast to the aforementioned studies, this is the first report on expression of UTMP in tissues other than endometrium.
In this study, we aimed to investigate the association between UTMP and milk production traits and PL in dairy cattle. To do this, we identified 2 synonymous SNP at positions 1179 (A/G) and 1296 (A/G) using the pooled DNA sequencing approach (Kwok et al., 1994; Khatib et al., 2006). The SNP at position 1296 was associated with a significant increase in PL. This finding parallels results previously obtained for the PI locus, which was also found to be associated with PL in dairy cattle (Heyen et al., 1999; Khatib et al., 2005). Taken together, these results provide strong evidence that the UTMP region is associated with the PL trait in dairy cattle. Unfortunately, these studies are not able to determine whether the effect observed is due to polymorphism in PI, UTMP, or another locus in that region. Additional research on this topic should be undertaken.
There is accumulating evidence that allelic variation in gene expression is correlated with phenotypic variability of quantitative and qualitative traits. In humans, only a few studies have addressed the possible roles of allele-specific expression in phenotypic variability among individuals and in disease susceptibility (Duno et al., 2004; Hirota et al., 2004; Murakami et al., 2004). For example, Duno et al. (2004) examined DAE of the CLCN1 gene, which causes myotonia congenita, an inherited disorder characterized by delayed skeletal muscle relaxation. Interestingly, they found that the expression level of the mutant allele was higher than that of the wild-type allele, and they suggested that differential expression of CLCN1 alleles might have roles in progression of myotonia congenital disease.
In this study, a polymorphism-based approach was used to analyze the allelic variation in expression of the bovine UTMP gene in reproductive tissues obtained from 9 heterozygous individuals. Reverse transcription-PCR analysis revealed that UTMP is expressed mainly in ovary, endometrium, and caruncle tissues. The DAE of UTMP in these tissues was measured as an allelic ratio of G/A alleles at SNP 1296. All heterozygous individuals examined displayed DAE, with allele G preferentially expressed in all endometrium tissues. Monoallelic expression (allele G) was observed in 2 of 3 ovary tissues and in 2 of 4 caruncle tissues. It is generally believed that genes displaying DAE would cause phenotypic variability across individuals. Given that UTMP has possible roles in health traits and based on the constant preferential expression of G allele at SNP 1296, it is plausible that this SNP is linked to improved PL.
Thus, the DAE of the bovine UTMP gene observed in this study is consistent with other studies on different genes that have also shown a correspondence between allelic variation in gene expression and phenotypic variability.
In summary, in this study we report that UTMP is primarily expressed in the bovine endometrium, ovary, and caruncle tissues. Also, we report that UTMP SNP 1296 shows significant association with PL in 2 independent Holstein cattle populations. This SNP could be exploited by marker-assisted selection for genetic improvement of PL within a breeding nucleus or commercial population; or by marker-assisted introgression for transferring desired alleles from a resource population to a commercial population.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication November 1, 2006. Accepted for publication January 9, 2007.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. Khatib, C. Maltecca, R. L. Monson, V. Schutzkus, X. Wang, and J. J. Rutledge The fibroblast growth factor 2 gene is associated with embryonic mortality in cattle J Anim Sci, September 1, 2008; 86(9): 2063 - 2067. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang, C. Maltecca, R. Tal-Stein, E. Lipkin, and H. Khatib Association of Bovine Fibroblast Growth Factor 2 (FGF2) Gene with Milk Fat and Productive Life: An Example of the Ability of the Candidate Pathway Strategy to Identify Quantitative Trait Genes J Dairy Sci, June 1, 2008; 91(6): 2475 - 2480. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Khatib, R. L. Monson, V. Schutzkus, D. M. Kohl, G. J. M. Rosa, and J. J. Rutledge Mutations in the STAT5A Gene Are Associated with Embryonic Survival and Milk Composition in Cattle J Dairy Sci, February 1, 2008; 91(2): 784 - 793. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |