Association of dopaminergic pathway gene polymorphisms with chronic renal insufficiency among Asian Indians with type-2 diabetes

Background Genetic markers conferring susceptibility to diabetes specific renal disease remains to be identified for early prediction and development of effective drugs and therapies. Inconsistent results obtained from analysis of genes from classical pathways generate need for examination of unconventional genetic markers having role in regulation of renal function. Experimental and clinical evidences suggest that dopamine is an important natriuretic hormone. Therefore, various genes involved in regulation of dopamine bioavailability could play a role in diabetic chronic renal insufficiency (CRI). We investigated the contribution of 12 polymorphisms from five Dopaminergic pathway genes to CRI among type-2 diabetic Asian Indian subjects. Methods Genetic association of 12 polymorphisms (SNPs) from five genes namely-dopamine receptor-1 (DRD1), DRD2, DRD3, DRD4, andcatechol-O-methyltransferase (COMT) with diabetic CRI was investigated using a case-control approach. Logistic regression analysis was carried out to correlate various clinical parameters with genotypes, and to study pair wise interactions between SNPs of different genes. Results SNPs -141 ins/del C and G>A (1 kb upstream from exon 2) in DRD2 gene showed significant allelic and genotypic association. Allele -141 insC and genotype -141 insC/insC of -141 ins/del C polymorphism, and allele A of G>A SNP were found to be predisposing to CRI. Our result of allelic and genotypic association of -141 insC/delC SNP was also reflected in the haplotypic association. Heterozygous genotype of polymorphism 900 ins/del C in COMT gene was predisposing towards CRI. Conclusion Some polymorphisms in DRD2 and COMT genes are significantly associated with susceptibility to CRI in the Asian Indian population which, if confirmed would be consistent with a suggested role of dopamine metabolism in disease occurrence.


Background
A large number of genetic association studies for unraveling the genetic basis of renal disease in diabetes have been carried out using polymorphisms in candidate genes from various postulated biochemical and metabolic pathways. These include renin-angiotensin-aldosterone system (RAAS), protein kinase C (PKC)-polyol, advanced glycation end (AGE) product, hexosamine and oxidative stress. However the results are confusing and identification of genetic susceptibility to chronic renal insufficiency (CRI) continues to be a challenge. This suggests the need for examining other unconventional genetic markers having role in regulation of normal physiological functions of kidney, especially under hyperglycemic condition.
Dopamine is an important natriuretic hormone involved in the regulation of blood pressure, and salt-water reabsorption in kidney. While angiotensin II, another hormone involved in natriuresis, stimulates sodium reabsorption through activation of Na-K-ATPase and Na/ H exchanger, dopamine acts antagonistically and inhibits reabsorption by inhibiting Na-K-ATPase and Na/H exchanger, regulating sodium transport and excretion in the proximal tubule [1][2][3]. Moreover, there are ample empirical evidences to support that dopamine through specific dopamine receptors influences various components of renin angiotensin aldosterone system (RAAS). While D1 like receptors stimulate renin release, D3 receptor inhibits it [4]. D3 receptor of dopamine decreases angiotensin II type 1 receptor (AT1) expression and angiotensin II binding to AT1 in renal proximal tubules [3,4]. Disruption of D3 receptor gene in mice induces a renin dependent form of hypertension [5]. Such interactions between dopamine and components of RAAS could modulate the potential effects of RAAS, which besides regulating vascular tone and Na ion balance in kidney also influences cytokine mediated inflammatory pathway in diabetic kidney disease [6][7][8][9]. In addition, administration of nitecapone, an inhibitor of dopamine metabolizing enzyme catechol-O-methyltransferase (COMT; EC 2.1.1.6), has been found to effectively abolish glomerular hypertension and reduce progression to glomerulosclerosis by inhibiting Na-K-ATPase in a dopamine dependent fashion [10,11]. These findings reiterate a potential role of dopaminergic pathway genes in conferring susceptibility to diabetic CRI through modulation of dopamine levels.
Several polymorphisms in dopaminergic pathway genes have been identified and extensively investigated for their role in hypertension and schizophrenia [12][13][14][15][16][17]. There is only one report on association of D3 receptor of dopamine with diabetic nephropathy but in individuals with type-1 diabetes [18]. Due to significant role of genes from dopaminergic pathway in kidney function, in this pilot study, we investigated the susceptibility conferred by polymorphisms in dopaminergic pathway genes to CRI among individuals with type-2 diabetes, hitherto unexplored.

Subjects
Complete clinical and demographic characteristics of the study population have been previously reported [19]. Briefly, in this retrospective case-control analysis consecutive subjects suffering from type-2 diabetes with CRI (cases, CRI, n = 196), and diabetics without any evidence of diabetic kidney disease (controls, DM, n = 225) were recruited from the outpatient departments of the four participating medical institutions situated across India. The study was approved by respective institutional ethics committee and an informed consent obtained from all the participants. For CRI, the inclusion criteria were subjects with type 2 diabetes for ≥2 years, plus two or more of the following: serum creatinine ≥150 umol/l, urinary albumin excretion rate (UAER) > 200 mg/l, and presence of retinopathy. Retinopathy was diagnosed by either fundoscopic examination or fluoroangiographic study. 10 ml venous blood was collected from each individual included in the study for biochemical and genetic analysis. Biochemical analyses to determine fasting glucose, glycated hemoglobin, serum creatinine, triglycerides, total cholesterol, and albumin were carried out at the respective centres using automated analyzers. Using serum creatinine as a surrogate marker we calculated glomerular filtration rate (GFR) by using the online modified diet in renal disease (MDRD) calculator. Patients with drug induced nephrotoxic damage or secondary causes of albuminuria such as obstructive renal disease, renal stone disease and acute urinary tract infection were excluded from this group. Normoalbuminuric (AER < 20 mg/l) type-2 diabetes subjects of ≥ 10 years duration of diabetes (mean duration 17.07 ± 6.69 years) were recruited as control diabetes subjects. Normoalbuminuria among control (DM) subjects was determined by timed urine collections. An aliquot of blood from the four centres was transported to the genetics laboratory for DNA isolation. DNA was isolated from the lymphocytes using the conventional phenolchloroform organic extraction method and used for genetic analysis.
Genotyping was done using either polymerase chain reaction (PCR)-restriction fragment length polymorphisms (RFLP) or PCR-length polymorphism (LP) approach following the method described in the respective quoted references. Genotypic profiles obtained for each of the polymorphisms are presented in Additional file 1. The digested PCR products were resolved on 2-3% agarose gel stained with ethidium bromide.

Statistical Analysis
Comparison of all clinical variables between DM and CRI subjects were carried out by χ 2 test for nominal variables or t-test for continuous variables. Continuous variables where skewed distribution was observed were compared by Mann Whitney U test and values are reported as mean and range. Hardy-Weinberg equilibrium (HWE) was tested for each of the genetic marker. Allelic and genotypic associations of SNPs/VNTRs were evaluated by Pearson's χ 2 test followed by odds ratio and 95% CI computation.
Power of the sample size for each of the SNPs was calculated using PAWE software version 1.2 [25,26]. Haplotype analysis was performed using PHASE-standard analysis version 2.0.2 [27,28]. Chi-square values were derived from a series of 2 × 2 contingency tables based on the frequency of each haplotype versus all others between the DM and the CRI groups. Linear and logistic regression analyses were carried out to correlate various clinical parameters with genotypes and to study pair wise interactions between SNPs of different genes. P values < 0.05 were subject to Bonferroni's correction.

Clinical analysis
Clinical characteristics of the two patient groups are presented in Table 1.

Genetic analysis: Efficiency of PCR based genotyping was approximately 98% with at least two additional attempts for genotyping of failed samples COMT
No allelic association of any of the four SNPs in this gene with diabetic CRI was observed ( Table 2). Only the heterozygote genotype of 900 ins/del C SNP was significantly associated with CRI (P = 0.018, OR 1.66, 95%CI 1.09-2.52). Haplotypes generated using genotypic profiles of the four SNPs in COMT gene are presented in Table 3. We observed significant association of three haplotypes in COMT gene, which were retained even after Bonferroni's correction.

DRD1
No allelic or genotypic association of the -48 A>G polymorphism with CRI in our sample set was seen ( Table 2).

DRD2
Of the three polymorphisms analysed in this gene, we observed a highly significant allelic and genotypic association of -141 ins/del C and G>A (1 kb upstream from exon 2) SNPs. T>C (10 kb downstream from exon 8) SNP did not show any allelic or genotypic association with diabetic kidney disease ( Table 2). Out of the eight possible haplotypes in our sample set, three (-141 insC-A-T, -141 insC-A-C, and -141 delC-A-C) were significantly associated with CRI (Table 4).

DRD3
No allelic or genotypic association of the most commonly investigated Ser9Gly SNP with CRI was observed in this study (Table 2).

DRD4
None of the three polymorphisms tested in the gene showed significant genotypic or allelic association with CRI. Out of 21 possible haplotypes in our case-control population, only haplotype 2-1-3 [(-120) no duplication -(-521) C -VNTR 467] was found to be marginally associated with CRI. However, after Bonferroni's correction the haplotype did not retain significance (Additional file 2).    27-8.84) and AA (P = 0.013; OR 1.88, 95%CI 1.14-3.10) of SNPs -141 ins/del C, and G>A (1 kb upstream from exon 2) respectively of DRD2 gene to be predisposing to CRI. Pair-wise interactions between different polymorphisms included in this study were tested using multiple logistic regression analysis. No significant interaction between any of these polymorphisms/ genes was observed.

Discussion
In addition to its classical role as a neurotransmitter, dopamine also acts as a natriuretic hormone that has implications in kidney function. Dopamine down regulation causes Na retention and renal tissue injury [11]. In addition dopamine also modulates the activity and availability of Ang II, the key molecule of renin-angiotensinaldosterone system in kidney. An altered interplay of dopamine and Ang II could disturb the vascular tone, sodium ion balance, and enhance cytokine-mediated injury to renal tissue. Thus, as mentioned earlier, genetic variations in genes controlling the bio-availability of dopamine could be pathological. While various studies have implicated genetic variations in COMT and dopamine receptor genes in essential hypertension [29][30][31][32][33][34][35], and one study has implicated COMT (based on a functional assay using a COMT inhibitor) in diabetic kidney disease in rats [10], an association of dopaminergic pathway gene polymorphisms with kidney disease among  a Order of SNPs in the DRD2 haplotypes: -141 ins/del C -G>A (Intron 1) -T>C (Intron 7) b Significant after Bonferroni's correction α = 0.0063 *n = total number **f = frequency individuals with type-2 diabetes has not been tested in any population so far. This is the first report on the role of multiple SNPs in dopaminergic pathway genes conferring susceptibility to diabetic CRI.
Out of four SNPs tested in COMT gene only 900 ins/del C showed a strong trend towards genotypic association (Table 2), with insC/delC genotype being predisposing (OR 1.66, CI 1.09-2.52) towards CRI. However, in the absence of any evidence for the functional status of this polymorphism at present, it would be difficult to substantiate its role in CRI. Val158Met polymorphism has been the most widely implicated for the activity status of COMT. Presence of methionine instead of valine in the protein makes it heat-labile and leads to fourfold reduction in activity [36]. In addition, different studies indicate that this functional polymorphism accounts for most of the variation in peripheral COMT activity [37]. Though, there was no allelic or genotypic association of this SNP ( Table 2) in our study, in haplotypic combinations this polymorphism seems to be important. Haplotype 1-1-2-2 [A-C-A (Met)-del C] is predisposing and haplotype 1-1-1-2 [A-C-A (Val)-del C], with only one allele different at Val158Met locus, is protective to CRI. 158Val, the high activity allele of the enzyme, is the predominant allele in our population and is required for regulating normal physiological levels of dopamine. Replacement of Val allele by Met in the haplotype may increase dopamine levels, which is an antagonist to Ang II, and may alter vascular tone and sodium ion balance in kidney. The only other significantly associated haplotype in COMT gene, 2-1-2-1 (G-C-A-InsC), is protective towards CRI (Table 3). In the absence of any information on functional status of the other three polymorphisms in this gene, it seems difficult to assess the protective nature of this haplotype.
Out of eight polymorphisms in four dopamine receptor genes (D1-D4), only SNPs -141 ins/del C and G>A (1 kb upstream from exon 2) in DRD2 gene, showed significant allelic and genotypic association (   (Table 4). Based on functional analysis it has been reported that -141 delC allele has decreased promoter strength as compared to the -141 ins C allele [38]. Considering the above mentioned facts, it could be argued that -141 ins C allele may lead to excess bio-availability of dopamine, which may alter vascular tone and sodium ion balance in kidney. Function of G>A SNP (I Kb upstream from exon 2) in the gene is not known and thus it's difficult to comment on our observation on this SNP.
There are a few study limitations. Present study is a candidate gene based case-control study and most candidate gene case-control studies of complex traits to date have been disappointing. In diabetic nephropathy alone, many initial reports of significant associations have not withstood the test of replication in other cohorts. Such inconsistencies stem primarily from the very low pre-test probability that any given gene contributes to the susceptibility of a complex trait despite apriori hypotheses based on cell, tissue, or animal model experiments. Other reasons include the use of underpowered sample sizes, multiple testing, phenotypic heterogeneity, poor phenotype characterization, selection bias, population stratification, and incomplete knowledge of the complete set of allelic variants in the region of a candidate gene. Our study involves a small sample size. Harrington has commented on standards of evidence in the genomic age and suggests that in a case control study the sample size should be >1000 subjects to arrive at worthwhile conclusions [39]. However, only a very few investigators have access to such a large data and inspite of being a multi-site study the present study could not achieve such large numbers. On the other hand the present study fulfils most of the criteria of a good genetic association study suggested by Hattersley et al [40]. The diagnosis of diabetic renal disease in the present study is based on multiple criteria. Diabetes can produce a large variety of renal lesions and it is likely that molecular and genetic pathways involved in each of these are different. This is the first investigation of association of polymorphisms in the dopaminergic pathway genes with CRI. Significant allelic, genotypic, and haplotypic associations observed with diabetic CRI in this study seem interesting but extrapolation of the results demands replication in a larger sample set.

Conclusion
It is concluded that polymorphisms -141 ins/del C and G>A (1 kb upstream from exon 2) in DRD2, and 900 ins/ del C in COMT genes are associated with susceptibility to diabetic CRI, which if confirmed by replication studies and/or expression analysis would be consistent with a role of dopamine metabolism in disease occurrence. Given that some of the genetic variants were associated with susceptibility and others with protection could be an important indication of the mechanisms of CRI itself. The use of the apposite polymorphisms for prediction or prognosis could well be a subordinate step to that of the confirmation of the observation and of understanding the disease itself.