Linkage studies of catechol-O-methyltransferase (COMT) and dopamine-beta-hydroxylase (DBH) cDNA expression levels
© Xing et al; licensee BioMed Central Ltd. 2007
Published: 18 December 2007
The COMT and DBH genes are physically located at chromosomes 22q11 and 9q34, respectively, and both COMT and DBH are involved in catecholamine metabolism and are strong candidates for certain psychiatric and neurological disorders. Although the genetic determinants for both enzymes' activities have been widely studied, their genetic involvement on gene mRNA expression levels remains unclear. In this study we performed quantitative linkage analysis of COMT and DBH cDNA expression levels, identifying transcriptional regulatory regions for both genes. Multiple Haseman-Elston regression was used to detect both additive and interactive effects between two loci. We found that the master transcriptional regulatory region 20q13 had an additive effect on the COMT expression level. We also found that chromosome 19p13 showed both additive and interactive effects with 9q34 on DBH expression level. Furthermore, a potential interaction between COMT and DBH was indicated.
Catechol-O-methyltransferase (COMT) (MIM 116790) catalyzes the transfer of a methyl group from S-adenosyl-methionine to catecholamines , and dopamine-beta-hydroxylase (DBH) (MIM 223360) catalyzes the conversion of dopamine to norepinephrine ; both enzymes are involved in the catecholamine metabolism system, in which abnormalities are hypothesized to be related to the development of psychiatric and neurological disorders such as schizophrenia and depression (for reviews, see [3, 4]). There are two distinct forms of COMT, a membrane-bound COMT (MB-COMT) and a soluble COMT (S-COMT) found in the cell cytoplasm ; the former is predominantly expressed in brain neurons, whereas the latter is predominantly expressed in other tissues such as liver and kidney. DBH is expressed in the vesicles of central noradrenergic and adrenergic neurons, peripheral noradrenergic neurons, and adrenomedullary neuronsecreotry cells . The COMT and DBH genes are physically located at chromosomes 22q11 and 9q34, respectively [7, 8].
There has been intensive research on the genetic determinants of plasma COMT and DBH activities, and findings so far strongly suggest that the structural genes COMT and DBH encoding these proteins are the major quantitative trait loci for their respective plasma activities; moreover, in both genes single-nucleotide polymorphisms (SNPs) have been identified that may account for most of the variation in the enzyme activity levels. On the common functional SNP rs165688, a G→A substitution in exon 4 of COMT produces a valine→methionine (Val/Met) substitution at codons 158 and 108 in the MB-COMT and S-COMT transcripts, respectively. This polymorphism has been shown to have a functional effect on enzyme activity; in particular, Val is a predominant factor that determines higher COMT activity . The SNP rs1611115 (-1021C→T) located in the 5' flanking region of the DBH gene accounts for 35% to 52% of the variance in plasma DBH activity in populations from diverse geographic origins .
Though the major effects of these common variants on enzyme activities have been replicated in independent samples, their effects on gene mRNA expression levels remain unclear. The effect of SNP rs165688 on COMT mRNA expression level has been controversial. Bray et al.  showed that it altered the COMT mRNA expression level; in contrast, several other studies failed to detect this association between COMT genotypes and expression level (for examples, see ). So far, studies on the determinants of COMT expression level have been limited to seeking cis-acting regulators, including SNP rs165688 and several other SNPs within COMT. To our knowledge, there is as yet no literature documenting effects of SNPs on DBH mRNA expression levels.
Phenotypic variation in the transcript level represents an intermediate stage between DNA sequence differences and the clinical expression of complex traits, and there have been studies suggesting a genetic contribution to the variation in gene expression level [13–15]. In the current study, we performed quantitative linkage analysis for both COMT and DBH expression phenotypes, identifying transcriptional regulatory regions for both genes. For COMT, we first identified cis-acting influences on the gene expression level, and then we performed a genome-wide linkage scan for trans-acting regulators across the genome conditional on the cis-acting effect at chromosome 22q11. For DBH, we similarly first identified the cis-acting influences on the gene expression level, and then, both conditional and unconditional on the cis-acting effect at chromosome 9q34, performed a linkage scan for trans-acting regulators specifically on chromosome 19, based on the evidence from previous linkage studies [16–18] that there might be a second locus near chromosome 19p13 acting jointly with DBH at chromosome 9q34 to influence DBH plasma activity.
where γ denotes variance (again, attenuated by recombination) due to the genetic interaction between locus 1 and locus 0, and a value significantly different from zero indicates an interactive effect of the two loci on gene expression level. Single-point IBD sharing was estimated by the program GENIBD. Both of the above programs are included in the Statistical Analysis for Genetic Epidemiology (S.A.G.E.) software .
Linkage of COMTexpression levels
Chromosome regions showing only additive trans-acting effectsa with 22q11 on COMT expression level
Physical distance (Mb)
p-Value for λ
4.4 × 10-3
3.4 × 10-4
4.1 × 10-3
2.5 × 10-3
3.7 × 10-3
1.8 × 10-3
2.2 × 10-3
2.0 × 10-7
2.1 × 10-4
3.1 × 10-3
Chromosome regions showing only interactive trans-acting effectsa with 22q11 on COMT expression level
Physical distance (Mb)
p-Value for γ
1.9 × 10-3
1.2 × 10-5
3.9 × 10-4
2.1 × 10-3
8.3 × 10-6
3.3 × 10-3
1.5 × 10-4
Chromosome regions showing both additive and interactive trans-acting effectsa with 22q11 on COMT expression level
Physical distance (Mb)
1.1 × 10-4
4.6 × 10-3
<1.0 × 10-8
3.5 × 10-2
2.3 × 10-2
2.9 × 10-5
1.5 × 10-3
5.0 × 10-7
1.7 × 10-2
1.8 × 10-3
1.1 × 10-3
4.7 × 10-2
5.7 × 10-3
3.5 × 10-3
4.7 × 10-2
2.9 × 10-5
1.3 × 10-4
6.6 × 10-3
1.4 × 10-4
1.1 × 10-2
1.0 × 10-3
5.9 × 10-3
4.7 × 10-3
2.6 × 10-3
4.8 × 10-2
4.4 × 10-3
4.8 × 10-4
8.6 × 10-3
Linkage of DBHexpression levels
Phenotypic variation in the transcript level of a gene involves genetic components; however, the transcriptional regulatory elements may or may not lie in the structural gene itself. Studies by Morley et al.  have suggested the existence of master regulatory regions, and we identified chromosome 20q13, nominated by them, as regulating COMT expression level. Note that Morley et al.  performed simple HE regression, while we performed multiple HE regression so as to condition on the IBD sharing at COMT, i.e., our analysis refined their original investigation and provided further evidence supporting chromosome 20q13 as a regulatory region. In the current data, overall there was a linkage signal at 22q11 for COMT expression level, but no signal at 9q34 for DBH expression level. Lack of an overall linkage signal at 9q34 for DBH could be due to there not being enough DBH polymorphism segregation in the sample, to there being less expression variation among individuals than between replicates within an individuals , from low levels of transcript abundance in view of its tissue-specific expression pattern (Dr. Joseph F. Cubells, personal communication), or to genetic heterogeneity, as also indicated in the earlier study of plasma DBH activity .
In this study we employed multiple HE regression to detect additive and interactive effects between two loci. In the three pedigrees selected for having a positive linkage signal at 9q34 for DBH expression level, we detected a boost of signal at chromosome 19p13, where previous studies [16–18] indicated the existence of a locus influencing DBH plasma activity. The multiple HE regression indicated both additive and interactive effects of the two loci, suggesting epistasis between them. For simplicity, and because of the limited amount of data available, we did not dissect the total genetic variance into additive and dominance components in this study; with more data the interaction could be characterized by more elaborate models , and even embedded in a general pedigree framework [25, 26]. We also noticed that chromosome 9q34, where DBH is located, showed an additive trans-acting effect with 22q11 on COMT expression level. This is not totally unexpected because both enzymes are involved in catecholamine metabolism; however, we could not find in the literature any genetic study on the interaction between them, which is surprising. Interaction between COMT and DBH could be a potential research topic.
In summary, we performed quantitative linkage analysis for both COMT and DBH expression phenotypes, identifying transcriptional regulatory regions for both genes. Multiple HE regression was used to detect additive and interactive effects between two loci. The master transcriptional regulatory region 20q13 had an additive effect on the COMT expression level. Chromosome regions 9q34 and 19p13 acted epistatically on DBH expression level. Potential interaction between COMT and DBH was indicated.
Some of the results in this paper were obtained by using the program package S.A.G.E., which is supported by a U.S. Public Health Service Resource Grant (RR03655) from the National Center for Research Resources. This work was supported in part by research grants GM-28356 from the National Institute of General Medical Sciences, Cancer Center Support Grant P30CAD43703 from the National Cancer Institute, and a Training grant (HL07567) from the National Heart, Lung and Blood Institute. CX was supported by a Merck Foundation fellowship. CX and QL were sponsored to attend GAW15 by the Endowment Sponsored Mentorship Program of Case Western Reserve University Graduate School.
This article has been published as part of BMC Proceedings Volume 1 Supplement 1, 2007: Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/1?issue=S1.
- Axelrod J, Tomchick R: Enzymatic O-methylation of epinephrine and other catechols. J Biol Chem. 1958, 233: 702-705.PubMedGoogle Scholar
- Kaufman S, Friedman S: Dopamine-beta-hydroxylase. Pharmacol Rev. 1965, 17: 71-100.PubMedGoogle Scholar
- Gogos JA, Gerber DJ: Schizophrenia susceptibility genes: emergence of positional candidates and future directions. Trends Pharmacol Sci. 2006, 27: 226-233. 10.1016/j.tips.2006.02.005.View ArticlePubMedGoogle Scholar
- Cubells JF, Zabetian CP: Human genetics of plasma dopamine beta hydroxylase activity: application to research in psychiatry and neurology. Psychopharmacology (Berl). 2004, 174: 463-476. 10.1007/s00213-004-1840-8.View ArticleGoogle Scholar
- Tenhunen J, Salminen M, Jalanko A, Ukkonen S, Ulmanen I: Structure of the rat catechol-O-methyltransferase gene: separate promoters are used to produce mRNAs for soluble and membrane-bound forms of the enzyme. DNA Cell Biol. 1993, 12: 253-263.View ArticlePubMedGoogle Scholar
- Kemper CM, O'Connor DT, Westlund KN: Immunocytochemical localization of dopamine-beta-hydroxylase in neurons of the human brain stem. Neuroscience. 1987, 23: 981-989. 10.1016/0306-4522(87)90173-4.View ArticlePubMedGoogle Scholar
- Grossman MH, Emanuel BS, Budarf ML: Chromosomal mapping of the human catechol-O-methyltransferase gene to 22q11.1-q11.2. Genomics. 1992, 12: 822-825. 10.1016/0888-7543(92)90316-K.View ArticlePubMedGoogle Scholar
- Craig SP, Buckle VJ, Lamouroux A, Mallet J, Craig IW: Localization of the human dopamine beta hydroxylase (DBH) gene to chromosome 9q34. Cytogenet Cell Genet. 1988, 48: 48-50.View ArticlePubMedGoogle Scholar
- Lotta T, Vidgren J, Tilgmann C, Ulmanen I, Melen K, Julkunen I, Taskinen J: Kinetics of human soluble and membrane-bound catechol-O-methyltransferase: a revised mechanism and description of the thermolabile variant of the enzyme. Biochemistry. 1995, 34: 4202-4210. 10.1021/bi00013a008.View ArticlePubMedGoogle Scholar
- Zabetian CP, Anderson GM, Buxbaum SG, Elston RC, Ichinose H, Nagatsu T, Kim KS, Kim CH, Malison RT, Gelernter J, Cubells JF: A quantitative-trait analysis of human plasma-dopamine beta-hydroxylase activity: evidence for a major functional polymorphism at the DBH locus. Am J Hum Genet. 2001, 68: 515-522. 10.1086/318198.View ArticlePubMed CentralPubMedGoogle Scholar
- Bray NJ, Buckland PR, Williams NM, Williams HJ, Norton N, Owen MJ, O'Donovan MC: A haplotype implicated in schizophrenia susceptibility is associated with reduced COMT expression in human brain. Am J Hum Genet. 2003, 73: 152-161. 10.1086/376578.View ArticlePubMed CentralPubMedGoogle Scholar
- Chen J, Lipska BK, Halim N, Ma QD, Matsumoto M, Melhem S, Kolachana BS, Hyde TM, Herman MM, Apud J, Egan MF, Kleinman JE, Weinberger DR: Functional analysis of genetic variation in catechol-O-methyltransferase (COMT): effects on mRNA, protein, and enzyme activity in postmortem human brain. Am J Hum Genet. 2004, 75: 807-821. 10.1086/425589.View ArticlePubMed CentralPubMedGoogle Scholar
- Yan H, Yuan W, Velculescu VE, Vogelstein B, Kinzler KW: Allelic variation in human gene expression. Science. 2002, 297: 1143-10.1126/science.1072545.View ArticlePubMedGoogle Scholar
- Cheung VG, Conlin LK, Weber TM, Arcaro M, Jen KY, Morley M, Spielman RS: Natural variation in human gene expression assessed in lymphoblastoid cells. Nat Genet. 2003, 33: 422-425. 10.1038/ng1094.View ArticlePubMedGoogle Scholar
- Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, Spielman RS, Cheung VG: Genetic analysis of genome-wide variation in human gene expression. Nature. 2004, 430: 743-747. 10.1038/nature02797.View ArticlePubMed CentralPubMedGoogle Scholar
- Wilson AF, Elston RC, Siervogel RM, Tran LD: Linkage of a gene regulating dopamine-beta-hydroxylase activity and the ABO blood group locus. Am J Hum Genet. 1988, 42: 160-166.PubMed CentralPubMedGoogle Scholar
- Wilson AF, Elston RC, Sellers TA, Bailey-Wilson JE, Gersting JM, Deen DK, Sorant AJM, Tran LD, Amos CI, Siervogel RM: Stepwise oligogenic segregation and linkage analysis illustrated with dopamine-betahydroxylase activity. Am J Med Genet. 1990, 35: 425-432. 10.1002/ajmg.1320350321.View ArticlePubMedGoogle Scholar
- Province M: Sequential methods of analysis for genome scans. Adv Genet. 2001, 42: 499-514.View ArticlePubMedGoogle Scholar
- Matise TC, Sachidanandam R, Clark AG, Kruglyak L, Wijsman E, Kakol J, Buyske S, Chui B, Cohen P, de Toma C, Ehm M, Glanowski S, He C, Heil J, Markianos K, McMullen I, Pericak-Vance MA, Silbergleit A, Stein L, Wagner M, Wilson AF, Winick JD, Winn-Deen ES, Yamashiro CT, Cann HM, Lai E, Holden AL: A 3.9-centimorgan-resolution human single-nucleotide polymorphism linkage map and screening set. Am J Hum Genet. 2003, 73: 271-284. 10.1086/377137.View ArticlePubMed CentralPubMedGoogle Scholar
- Haseman JK, Elston RC: The investigation of linkage between a quantitative trait and a marker locus. Behav Genet. 1972, 2: 3-19. 10.1007/BF01066731.View ArticlePubMedGoogle Scholar
- Shete S, Jacobs KB, Elston RC: Adding further power to the Haseman and Elston method for detecting linkage in larger sibships: weighting sums and differences. Hum Hered. 2003, 55: 79-85. 10.1159/000072312.View ArticlePubMedGoogle Scholar
- Liang K-Y, Zeger SL: Longitudinal data analysis using generalized linear models. Biometrika. 1986, 73: 13-22. 10.1093/biomet/73.1.13.View ArticleGoogle Scholar
- S.A.G.E. Statistical analysis for genetic epidemiology. Version 5.2. [http://darwin.case.edu/sage]
- Tiwari HK, Elston RC: Linkage of multilocus components of variance to polymorphic markers. Ann Hum Genet. 1997, 61: 253-261. 10.1017/S0003480097006167.View ArticlePubMedGoogle Scholar
- Wang T, Elston RC: Two-level Haseman-Elston regression for general pedigree data analysis. Genet Epidemiol. 2005, 29: 12-22. 10.1002/gepi.20075.View ArticlePubMedGoogle Scholar
- Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.View ArticlePubMed CentralPubMedGoogle Scholar
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.