Genetic diversity of tyrosine hydroxylase (TH) and dopamine b-hydroxylase (DBH) genes in cattle breeds

Genetic diversity of tyrosine hydroxylase (TH) and dopamine b-hydroxylase (DBH) genes in cattle breeds

Diana Lelidett Lourenco-Jaramillo^I; Ana María Sifuentes-Rincón^I; Gaspar Manuel Parra-Bracamonte^I; Xochitl Fabiola de la Rosa-Reyna^I; Aldo Segura-Cabrera^II; Williams Arellano-Vera^I

^ILaboratorio de Biotecnología Animal, Centro de Biotecnología Genómica, Tamaulipas, México

^IILaboratorio de Biomedicina, Centro de Biotecnología Genómica, Tamaulipas, México

^{Send correspondence to}

ABSTRACT

DNA from four cattle breeds was used to re-sequence all of the exons and 56% of the introns of the bovine tyrosine hydroxylase (TH) gene and 97% and 13% of the bovine dopamine b-hydroxylase (DBH) coding and non-coding sequences, respectively. Two novel single nucleotide polymorphisms (SNPs) and a microsatellite motif were found in the TH sequences. The DBH sequences contained 62 nucleotide changes, including eight non-synonymous SNPs (nsSNPs) that are of particular interest because they may alter protein function and therefore affect the phenotype. These DBH nsSNPs resulted in amino acid substitutions that were predicted to destabilize the protein structure. Six SNPs (one from TH and five from DBH non-synonymous SNPs) were genotyped in 140 animals; all of them were polymorphic and had a minor allele frequency of > 9%. There were significant differences in the intra- and inter-population haplotype distributions. The haplotype differences between Brahman cattle and the three B. t. taurus breeds (Charolais, Holstein and Lidia) were interesting from a behavioural point of view because of the differences in temperament between these breeds.

Key words: behavior, B. t. indicus, B. t. taurus, catecholamine, temperament.

Basic cattle behavioural traits include social behaviour such as aggression and temperament (Buchenauer, 1999). These traits may have a direct economic impact and can be included in selection strategies (Mormède, 2005; Nkrumah et al., 2007; Core et al., 2009). Studies that have investigated an association between behaviour and production traits have found that nervous or aggressive animals have a decreased milk flow and yield (Rushen, cited by Hiendleder et al., 2003). Temperament, which is defined as the animal's response to handling by humans (Burrow and Dillon, 1997), has been studied in different breeds, and has been shown to affect growth, health, performance and carcass merit (Fox et al., 2004; Nkrumah et al., 2004, 2007; Core et al., 2009).

There are important temperament differences between Bos taurus taurus and Bos taurus indicus breeds and their crosses. Grandin (1980) found that the flight distances of Brahman cattle were longer than those of British breed cattle. Hearnshaw and Morris (1984) reported significant differences in the temperament of B. t. taurus- and B. t. indicus-derived breeds and among calves sired by B. t. indicus bulls; in contrast, there were no significant differences amongst calves sired by Hereford bulls.

Genes that encode different regulatory enzymes, transporters and receptors of the serotonin and dopamine signalling pathways may have a marked influence on genetic variation and phenotypic effects. Tyrosine hydroxylase (TH) is the rate-limiting enzyme for dopamine biosynthesis (Marsden, 2006). Genetic polymorphisms in the TH gene have been associated with different neurological and psychiatric conditions in humans (Giegling et al., 2009). Another target gene is dopamine b-hydroxylase (DBH), which is one of the most important enzymes for noradrenaline (NA) production. In humans, the genetic variation in this gene has been associated with maternal behaviour, feeding behaviour, impulsive behaviour, motivation, behaviour reinforcement and stimuli that produce stress (Hoebel 1985; Thomas and Palmiter 1997; Szczypka et al., 1999; Adriani et al., 2008; Geiger et al., 2009; Luksys et al., 2009).

The TH and DBH genes have also been studied in non-human mammals. Takeuchi et al. (2005) screened for polymorphisms in brain samples from ten unrelated Beagles and found two single nucleotide polymorphisms (SNPs) in DBH and four in TH. These SNPs were subsequently studied in 193 samples from five dog breeds. With the exception of one SNP from each gene, there was significant variation in the allelic distributions among the breeds, which suggested (based on known functions of these genes) that these SNPs could be used as candidate markers to evaluate behavioural traits in dogs.

In this work, we examined the genetic variation in the DBH and TH genes of cattle breeds that display different temperaments depending on their genetic background (i.e., B. t. taurus and B. t. indicus). Eleven DNA samples from four cattle breeds (three Holstein, two Charolais, three Brahman and three Lidia) were used to assess genetic variation in the DBH and TH genes. Whole blood DNA was extracted from these cattle breeds and from additional bovine samples using commercial Wizard® Genomic DNA purification kits (Promega Corporation, Madison, WI, USA). The TH and DBH gene sequences obtained from GenBank (accession numbers NC_007330 and NC_007309, respectively) were used to design primers with Primer Select version 7.0.0 software (Lasergene, DNASTAR, Madison, WI USA). Table S1 shows that most of the primers were located at intron-exon borders. For TH, these primers facilitated the resequencing of 100% and 54% of the coding and non-coding sequences, respectively, whereas for DBH they facilitated the resequencing of 97% of the coding sequences (11 exons) and 13% of the non-coding sequences, which included a complete intron. Polymerase chain reaction (PCR) amplifications were done in 25 mL volumes on a Thermocycler DNA Engine TETRAD 2 Peltier thermal cycler (MJ Research, Inc.,Waltham, MA, USA). The reaction mixtures contained 5-100 ng of genomic DNA, 1.5 or 2.5 mM MgCl2, 0.1 mM of each primer, 0.4 mM dNTPs and 2.5 U of GoTaq polymerase (Promega). A touchdown method was used and the amplification profile included an initial denaturation step of 95 °C for 10 min, five three-step cycles of 45 s at 95 °C, an annealing step for 45 s that started at 68 °C but decreased by 2 °C during each cycle and 45 s at 72 °C, and 25 three-step cycles of 45 s each at 95 °C, 65 °C and 72 °C. Bidirectional sequencing was done using the BigDye® Terminator procedure and an ABI PRISM 3100 Genetic Analyzer DNA sequencer (Applied Biosystems, Foster City, CA, USA). The sequences were assembled into contigs using Cap Win V6.0.6000 (Huang and Madan, 1999) and were aligned with ClustalX 2.0.8 (Larkin et al., 2007).

The resulting sequence information was used as an initial polymorphism screening step. We compared the data obtained from 11 animals belonging to different breeds with sequences deposited in the NCBI database. In this analysis, the presence of SNPs was determined by visually inspecting the sequence chromatograms. The SNPs were defined according to their presence in the screening population using the three expected genotypes.

Tyrosine hydroxylase and dopamine b-hydroxylase are critical enzymes in the dopamine and norepinephrine pathways. The resequencing of 4.3 kb and 3.2 kb of the TH and DBH bovine genes, respectively, of 11 animals from four breeds revealed a microsatellite motif and 64 SNPs (22 located in the exons and 42 located in the introns of both genes).

For the amino acid conservation analysis, the published B. t. taurus (AF118638), Equus caballus (AB029430), Homo sapiens (BC017174), Mus musculus (S50200), Rattus norvegicus (L12407) and Canis familiaris (Q68CI2) DBH sequences were aligned using ClustalX 2.0.8. The possible effects of the amino acid substitutions introduced by the eight nsSNPs were predicted using the IPTREE-STAB server (Huang et al., 2007). This software predicted the stability of the mutant proteins produced by the eight nsSNPs. The parameters used for the IPTREE-STAB server included a pH of 7 and a temperature of 25 °C.

Four restriction fragment length polymorphism (PCR-RFLP) and two amplification-created restriction site (PCR-ACRS) assays were designed to genotype one TH SNP and five of the eight DBH non-synonymous SNPs (nsSNPs) (Table 1). One hundred and forty animals from the four breeds (Lidia, n = 19; Charolais, n = 50; Brahman, n = 50; Holstein, n = 21) were genotyped using these assays.

The observed genotype frequencies were compared with the expected frequencies based on the Hardy-Weinberg equilibrium using Genepop version 4.0.10 (Rousset, 2008). The EM algorithm in the Arlequin version 3.11 software package was used to estimate the maximum likelihood haplotype frequencies (Excoffier et al., 2005) within and between the four tested populations. A correspondence analysis done with the program JMP_v3.2.2 (SAS Institute Inc., Cary, NC, USA) was used to assess the influence of the haplotypes in the four breeds.

The coding sequence of DBH contained one SNP per 85 bp. This value was higher than the average for human genes (one SNP per 185 bp) and also higher than the average for some highly polymorphic genes, such as the bovine myostatin gene which contains one SNP per 100 bp (Dunner et al., 2003). High polymorphism rates in genes with important biological roles may be related to the selection of de novo mutations because this mechanism increases the variability of the gene in question (Dunner et al., 2003). As expected, most of the SNPs in the bovine DBH gene were located in introns (66%); however, a significant number, were also found in exons. The latter included 13 synonymous SNPs and eight nsSNPs, six of which resulted in non-conservative amino acid substitutions (Table 2). Analysis of these nsSNPs showed that the resulting amino acid changes could affect the stability of the DBH protein structure. Protein destabilization is a common mechanism by which amino acid substitutions cause human diseases (Teng et al., 2010). The effect of these changes on gene function and phenotype needs to be assessed in specific experiments designed to address this question. To assess amino acid conservation in the DBH nsSNPs we compared the DBH amino acid sequences across five mammalian species (Mus musculus, Rattus norvegicus, Homo sapiens, Equus caballus and Canis familiaris). Two of the eight nsSNPs (rs109440428 and rs110234325) encoded amino acid residues that were conserved among these species. Notably, only one amino acid residue in the previously reported bovine DBH gene sequence was conserved among the five species examined.

The level of polymorphism in the bovine TH coding regions corresponded to one variant for every 733 bp. This low level of polymorphism could be related to the functional role of this enzyme as the limiting step in the entire dopamine production process. Currently, GenBank reports the presence of 25 polymorphisms in the bovine TH gene. However, none of the SNPs found here has been reported before, even in other non-human mammals. Takeuchi et al. (2005) reported that four variants were located in the coding region of canine TH; however, their locations were different from those of the SNPs present in the bovine gene.

In TH, the allelic frequencies of SNP rs109268356 revealed that the C allele predominated in the four genotyped bovine populations. Since this is a synonymous mutation it could be considered to be inconsequential for the primary protein structure and/or function; however, a growing number of studies have shown that synonymous mutations are also under evolutionary pressure and that they can be implicated in disease (Hunt et al., 2009). Until recently, bovine behavioural traits had not been associated with a particular gene. However, a whole genome study by Hiendleder et al. (2003) identified a QTL on chromosome 29 that was linked to temperament. The authors proposed that TH could be one of several candidate genes that influence temperament in cattle. Our results provide additional support for the hypothesis that the TH gene is a strong candidate associated with bovine behaviour.

The novel microsatellite identified in TH was located approximately 84 bp from the start of exon 5 and was characterised by a cytosine mononucleotide core motif. This motif exhibited an adenine interruption where the cytosine residues showed length differences of 6 to 13 bp and were repeated 2 to 9 times in the 5' to 3' direction. We analysed this motif at the structural level and sequenced this region not only in Lidia, Charolais, Brahman and Holstein animals but also in Beefmaster, Bradford and Nellore animals. The subsequent sequence analyses revealed interbreed variability but no individual intrabreed variation. The structure of the motif was AGC6AC9 AG in Charolais, Hereford and Holstein animals, AGC8AC4AG in Lidia and Beefmaster animals, AGC12AC2AG3 in Brahman and Bradford animals and AGC13AC3AG3 in Nellore animals.

The presence of the microsatellite rs109798407 AGC6/8/12/13A(C)2/3/4/9 AG motif in the bovine TH gene has not been reported before. The presence of this microsatellite is very interesting because there are several examples of polymorphic microsatellites within some genes, such as bovine growth hormone (bGH), calpastatin (CAST) and insulin-like growth factor 1 (IGF1). Genetic diversity analyses of these genes have shown an association between microsatellites and both productive traits and the allelic distribution and segregation between different genetic groups, including B. t. taurus, B. t. indicus and their crosses (Levéziel et al., 1994; Nonneman et al., 1999; Hale et al., 2000; Curi et al., 2005). As described here, the bovine TH microsatellite shows interbreed variation and contains at least four allelic forms with changes in their structure and length. Notably, the analysis of individuals from the eight breeds studied revealed no intrabreed differences. Although in mammals most of the repetitive sequences are located within introns the C-monomers are present at a frequency of 11.4%. In humans, the TH gene contains a four-nucleotide motif (TCATn) located in the first intron; this motif acts as a regulatory element in gene transcription (Meloni et al., 1998). Further studies are needed to establish whether the microsatellite has any functional effect in bovine TH transcription and to determine whether there is any correlation between variation in microsatellite structure and productive traits in cattle.

The only synonymous TH SNP and five of the eight DBH nsSNPs were genotyped in the animals of four breeds, which allowed their validation as SNPs based on their allelic frequencies (Table 2). Similarly, there were important differences in the genetic diversity parameters of the DBH nsSNPs among the populations tested. Although the Brahman breed had the highest number of polymorphic loci (five), the genetic variation was higher in Holstein and Charolais animals (32.6% and 20.9%, respectively). This increased genetic variation most likely reflected differences in the heterozygosity values, which were higher in the latter two breeds than in Brahman animals.

Table 3 shows the haplotype distributions calculated based on the five DBH nsSNP genotypes. There were significant differences (p £ 0.0001) in the haplotype distributions among the B. t. taurus breeds and between the B. t. taurus and B. t. indicus breeds. The correspondence analysis based on these frequencies showed the clear separation of only B. t. indicus from breeds with a B. t. taurus background (Figure 1). Some shared haplotypes were identified among the four breeds. The breeds with a B. t. taurus background shared haplotype I, which was also the haplotype with the highest frequency among the breeds. Haplotype I was characterized by the absence of nsSNPs in the sequence (Table 3). Haplotype II was the second most frequent haplotype and was present in Charolais, Brahman and Holstein animals. Some exclusive haplotypes were also observed in each breed, with Brahman animals having the highest number (six haplotypes); haplotype VII was the most frequent (55%) and was characterized by the presence of three nsSNPs. Notably, we only observed haplotype VI in this breed; haplotype VI was characterised by the presence of five nsSNPs (Table 3).

There were significant differences in the haplotype frequencies within and among the breeds studied. These differences were more evident among the three B. t. taurus breeds when compared with the Brahman breed, which was considered to be representative of the B. t. indicus genetic background. The implications of these differences could be discussed from a behavioural perspective since there are important differences in temperaments both within and among cattle breeds. Previous work has shown that when Brahman crosses are restrained in a squeeze chute they show higher cortisol levels than English breed crosses (Grandin, 1997). For the moment, we are unable to associate any genotype with differences in temperament. However, the genetic diversity observed among the breeds examined in this study provides a basis for investigations into such associations. Since the overall aim of such studies is to identify DNA markers for cattle temperament, this trait may have direct economic value, and temperament phenotypes could be included in future selection strategies (Mormède, 2005; Nkrumah et al., 2007; Core et al., 2009).

Acknowledgments

The authors thank the different herd owners for allowing us to sample their cattle. This work was supported by Research Grant SIP-IPN 20090817.

Supplementary Material

The following online material is available for this article:

Table S1 - Primer sequences used to analyse bovine TH and DBH genes.

This material is available as part of the online article from http://www.scielo.br/gmb.

Received: June 10, 2011

Accepted: February 15, 2012.

Associate Editor: Alexandre Rodrigues Caetano

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

References

• Adriani W, Boyer F, Gioiosa L, Macri S, Dreyer J and Laviola G (2008) Increased impulsive behavior and risk proneness following lentivirus-mediated dopamine transporter over-expression in rats nucleus accumbens. Neuroscience 159:47-58.
• Buchenauer D (1999) Genetics of behaviour in cattle. In: Fries R and Rubinsky A (eds) The Genetics of Cattle. CABI Publishing, Wallingford, pp 365-390.
• Burrow H and Dillon R (1997) Relationships between temperament and growth in a feedlot and commercial carcass traits of Bos indicus crossbreds. Aust J Exp Agric 37:407-411.
• Core S, Widowski T, Manson G and Miller S (2009) Eye white percentage as a predictor of temperament in beef cattle. J Anim Sci 87:2168-2174.
• Curi R, Oliveira H, Silveira A and Lopes C (2005) Effects of polymorphic microsatellites in the regulatory region of IGF-1 and GHR on growth and carcass traits in beef cattle. Anim Genet 36:58-62.
• Dunner S, Miranda M, Amigues Y, Cañón J, Georges M, Hanset R, Williams J and Ménissier F (2003) Haplotype diversity of the myostatin gene among beef cattle breeds. Genet Sel Evol 35:103-118.
• Fox J, Carstens G, Brown E, White M and Woods S (2004) Residual feed intake of growing bulls and relationships with temperament, fertility and performance traits. J Anim Sci 82:2-6.
• Excoffier L, Laval G and Schneider S (2005) Arlequin v. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online 1:47-50.
• Geiger B, Haburcak M, Avena N, Moyer M, Hoebel B and Pothos E (2009) Deficits of mesolimbic dopamine neurotransmission in rat dietary obesity. Neuroscience 159:1193-1199.
• Giegling I, Moreno-De-Luca D, Calati R, Hartmann A, Moller H, De Ronchi D, Rujescu D and Serretti A (2009) Tyrosine hydroxylase and DOPA decarboxylase gene variants in personality traits. Neuropsychobiology 59:23-27.
• Grandin T (1980) Livestock behaviour as related to handling facilities design. Int J Study Anim Probl 1:33-52.
• Grandin T (1997) Assessment of stress during handling and transport. J Anim Sci 75:249-257.
• Hale C, Herring W, Shibuya H, Lucy M, Lubahn D, Keisler D and Jonson G (2000) Decreased growth in Angus steers with a short TG-microsatellite allele in the P1 promoter of the growth hormone receptor gene. J Anim Sci 78:2099-2104.
• Hearnshaw H and Morris CA (1984) Genetic and environmental effect on temperament score in beef cattle. Aust J Agric Res 35:723-733.
• Hiendleder S, Thomsen H, Reinsch N, Bennewitz J, Leyhe-Horn B, Looft C, Xu N, Medjugorac I, Russ I, Kuhn C, et al. (2003) Mapping of QTL for body conformation and behavior in cattle. J Hered 94:496-506.
• Hoebel G (1985) Brain neurotransmitters in food and drug reward. Am J Clin Nutr 42:1133-1150.
• Huang X and Madan A (1999) CAP3: A DNA sequence assembly program. Genome Res 9:868-877.
• Huang L, Gromiha M and Ho S (2007) IPTREE-STAB: Interpretable decision tree based method for predicting protein stability changes upon mutations. Bioinformatics 23:1292-1293.
• Hunt R, Sauna Z, Ambudkar S, Gottesman M and Kimchi-Sarfaty C (2009) Silent (synonymous) SNPs: Should we care about them? In: Anton A. Komar (ed) Single Nucleotide Polymorphisms. 2nd edition. Humana Press, New Jersey, pp 23-39.
• Larkin M, Blackshields G, Brown N, Chenna R, McGettigan P, McWilliam H, Valentin F, Wallace I, Wilm A, Lopez R, et al. (2007) Clustal W and Clustal X ver. 2.0. Bioinformatics 23:2947-2948.
• Levéziel H, Rodellar C, Leroux C, Pepin L, Grohs C, Vaiman D, Mahé M, Martin P and Grosclaude F (1994) A microsatellite within the bovine ? -casein gene reveals a polymorphism correlating strongly with polymorphisms previously described at the protein as well as the DNA level. Anim Genet 25:223-228.
• Luksys G, Gerstner W and Sandi C (2009) Stress, genotype and norepinephrine in the prediction of mouse behavior using reinforcement learning. Nat Neurosci 12:1180-1186.
• Marsden C (2006) Dopamine: The rewarding years. Br J Pharmacol 147:S136-S144.
• Meloni R, Albanece V, Ravassard P, Treilhou F and Mallet J (1998) A tetranucleotide polymorphic microsatellite, located in the first intron of the tyrosine hydroxylase gene, acts as a transcription regulatory element in vitro. Hum Mol Genet 7:423-428.
• Mormède P (2005) Molecular genetics of behaviour: Research strategies and perspectives for animal production. Livestock Prod Sci 93:15-21.
• Nkrumah J, Basarab J, Price M, Okine E, Ammoura A, Guercio S, Hansen C, Li C, Benkel B and Moore S (2004) Different measures of energetic efficiency and their relationships with growth, feed intake, ultrasound, and carcass measurements in hybrid cattle. J Anim Sci 82:2451-2459.
• Nkrumah J, Crews H, Basarab J, Price M, Okine E, Wang Z and Li C (2007) Genetic and phenotypic relationships of feeding behavior and temperament with performance, feed efficiency, ultrasound, and carcass merit of beef cattle. J Anim Sci 85:2382-2390.
• Nonneman D, Kappes S and Koohmaraie M (1999) A polymorphic microsatellite in the promoter region of the bovine calpastatin gene. J Anim Sci 77:3114-3115.
• Rousset F (2008) Genepop'007: A complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Res 8:103-106.
• Szczypka M, Rainey M, Kim D, Alaynick W, Marck B, Matsumoto A and Palmiter R (1999) Feeding behavior in dopamine-deficient mice. Proc Natl Acad Sci USA 96:12138-12143.
• Takeuchi Y, Hashizume C, Chon E, Momozawa Y, Masuda K, Kikusui T and Mori Y (2005) Canine tyrosine hydroxylase (TH) gene and dopamine b-hydroxylase (DBH) gene: Their sequences, genetic polymorphisms, and diversities among five different dog breeds. J Vet Med Sci 67:861-867.
• Teng S, Srivastava S and Wang L (2010) Sequence feature-based prediction of protein stability changes upon amino acid substitutions. BMC Genomics 11(Suppl 2):S5.
• Thomas S and Palmiter R (1997) Impaired maternal behavior in mice lacking norepinephrine and epinephrine. Cell 91:583-592.

Publication Dates

• Publication in this collection
  19 June 2012
• Date of issue
  2012

History

• Received
  10 June 2011
• Accepted
  15 Feb 2012