Susceptibility and Gene Interaction Study of the Angiotensin II Type 1 Receptor (AGTR1) Gene Polymorphisms with Non-Alcoholic Fatty Liver Disease in a Multi-Ethnic Population

Abstract

Angiotensin II type 1 receptor (AGTR1) has been reported to play a fibrogenic role in non-alcoholic fatty liver disease (NAFLD). In this study, five variants of the AGTR1 gene (rs3772622, rs3772627, rs3772630, rs3772633, and rs2276736) were examined for their association with susceptibility to NAFLD. Subjects made up of 144 biopsy-proven NAFLD patients and 198 controls were genotyped using TaqMan assays. The liver biopsy specimens were histologically graded and scored according to the method of Brunt. Single locus analysis in pooled subjects revealed no association between each of the five variants with susceptibility to NAFLD. In the Indian ethnic group, the rs2276736, rs3772630 and rs3772627 appear to be protective against NAFLD (p = 0.010, p = 0.016 and p = 0.026, respectively). Haplotype ACGCA is shown to be protective against NAFLD for the Indian ethnic subgroup (p = 0.03). Gene-gene interaction between the AGTR1 gene and the patatin-like phospholipase domain-containing 3 (PNPLA3) gene, which we previously reported as associated with NAFLD in this sample, showed a strong interaction between AGTR1 (rs3772627), AGTRI (rs3772630) and PNPLA3 (rs738409) polymorphisms on NAFLD susceptibility (p = 0.007). Further analysis of the NAFLD patients revealed that the G allele of the AGTR1 rs3772622 is associated with increased fibrosis score (p = 0.003). This is the first study that replicates an association between AGTR1 polymorphism and NAFLD, with further details in histological features of NAFLD. There is lack of evidence to suggest an association between any of the five variants of the AGTR1 gene and NAFLD in the Malays and Chinese. In the Indians, the rs2276736, rs3772630 and rs3772627 appear to protect against NAFLD. We report novel findings of an association between the G allele of the rs3772622 with occurrence of fibrosis and of the gene-gene interaction between AGTR1gene and the much-studied PNPLA3 gene.

About the Authors:

Shamsul Mohd Zain

* E-mail: [email protected]

Affiliation: The Pharmacogenomics Laboratory, Department of Pharmacology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Zahurin Mohamed

Affiliation: The Pharmacogenomics Laboratory, Department of Pharmacology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Sanjiv Mahadeva

Affiliation: Department of Medicine, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Sanjay Rampal

Affiliation: Julius Centre University of Malaya, Department of Social and Preventive Medicine, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Roma Choudhury Basu

Affiliation: Clinical Investigation Centre, University Malaya Medical Centre, Kuala Lumpur, Malaysia

Phaik-Leng Cheah

Affiliation: Department of Pathology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Agus Salim

Affiliation: Saw Swee Hock School of Public Health, National University of Singapore, Singapore, Singapore

Rosmawati Mohamed

Affiliation: Department of Medicine, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Introduction

The studies on non-alcoholic fatty liver disease (NAFLD) have rapidly increased since the term non-alcoholic steatohepatitis was first introduced by Ludwig et al [1]. NAFLD, defined as fatty accumulation of the liver without significant alcohol intake, ranges from simple steatosis to non-alcoholic steatohepatitis (NASH) and cirrhosis [2]. The histological features of NAFLD include steatosis, lobular inflammation, hepatocellular ballooning and fibrosis [3]. NAFLD is the hepatic manifestation of metabolic syndrome [4]–[8]. Although environmental factors, lifestyle and dietary intake underlie susceptibility to NAFLD, genetic factors play a significant upstream role [9]–[11].

The renin-angiotensin system acts on the hepatic stellate cells (HSC) to promote hepatic fibrosis [12]–[13]. The AGTR1 gene has been implicated with susceptibility to NAFLD [13]–[14]. To date, there is only one study of AGTR1 polymorphism and its association with the occurrence of NAFLD. The single nucleotide polymorphisms (SNPs) are located in the intron of AGTR1 on chromosome 3. In the study, five variants of AGTR1(rs3772622, rs3772627, rs3772630, rs3772633, and rs2276736) were revealed to be strongly associated with risk of NAFLD [15]. Furthermore, the deletion of angiotensin II type 1 receptor in animals was associated withdecreased hepatic steatosis. Thus, AGTR1 may be an important regulator of steatosis in the liver [16].

The data on AGTR1 gene polymorphisms and its association with NAFLD is currently limited to one study, therefore, we investigated the association between five previously reported SNPs of the AGTR1 gene with NAFLD susceptibility, in three major ethnic subgroups of the Malaysian population, namely the Malays, Chinese and Indians. These three ethnic groups, which are presumably of different genetic pools, provide a good setting to study possible ethnic differences in susceptibility to NAFLD. We further evaluated the association between the SNPs mentioned and the histological features of NAFLD, namely steatosis, lobular inflammation, hepatocellular ballooning and fibrosis. We previously reported in this same subjects, the association between PNPLA3 with NAFLD and disease severity [17]. In this study, we report new findings for the interaction between the AGTR1 gene and the PNPLA3 gene on susceptibility to NAFLD.

Methods

Ethics Statement

Subjects provided written informed consent after a full explanation of the research outline. The study protocol was reviewed and approved by the Medical Ethics Committee of University Malaya Medical Centre (UMMC).

Subjects

A case control study conducted at the UMMC recruited a total of 144 cases and 198 controls. All cases were biopsy-proven NAFLD patients based on increased echogenicity (compared to renal cortex) on ultrasound with or without abnormal alanine transferase (ALT). Ethnicity of the subjects were confirmed by affirmations of absence of mixed marriages for at least three generations.Exclusion criteria included alcohol consumption > 10g/day [18], hepatitis B or C infection, autoimmune hepatitis, those who had been on drugs that are known to cause steatosis, and Wilson's disease. All the biopsy specimens were on an average 1.5 cm long and contained six portal tracts. The biopsy specimens were graded according to the method of Brunt [19]–[20]. The controls were genetically unrelated healthy subjects presenting with body mass index (BMI) <25 kg/m², fasting plasma glucose <110 mg/dL, normal lipid profile and normal liver enzymes. NAFLD was excluded in the controls by ultrasonography according to the absence of the following criteria; diffuse increase in bright echoes in the liver parenchyma (exceeding that of renal cortex and spleen) with impaired visualization of the peripheral portal and hepatic vein borders, and/or loss of definition of the diaphragm and poor delineation of the portal and hepatic radicles.

Clinical and laboratory assessments

Venous blood samples were taken for the determination of haemoglobin A1c (HbA1c), total cholesterol, triglycerides, high-density lipoprotein cholesterol (HDL), low-density lipoprotein cholesterol (LDL), alanine transferase (ALT), aspartate aminotransferase (AST), and gamma-glutamyl transpeptidase (GGT) level in all subjects. The biochemical tests were analysed based on standard clinical laboratory methods. Other relevant parameters such as body mass index (BMI), waist circumference and blood pressure weredetermined according to standard protocol.

Genotyping

Blood samples were collected in EDTA tubesandgenomic DNA was prepared using the QiAamp DNA Mini Kit (Qiagen. Hilden, German) according to user protocol. The AGTR1 SNPs were genotyped using a predesigned TaqMan SNP genotyping assay (Applied Biosystems, Foster City, CA, USA) as previously described [17].

Statistical analyses

All values are presented as mean±standard deviation for continuous data and as percentages for categorical data. For normally distributed variables (age, BMI, waist circumference, ALT, HDL cholesterol, LDL cholesterol, total cholesterol, triglycerides, and systolic blood pressure), independentt-test was performed to determine the associations between the different grades of NAFLD. For variables that are not normally distributed (HbA1c, AST, GGT, diastolic blood pressure, steatosis grade, lobular inflammation, ballooning, and fibrosis), Mann-Whitney U test was performed. Hardy-Weinburg equilibirum (HWE) was checked for all the groups using a goodness-of fit χ2 test with one degree of freedom. A p value of less than 0.05 indicated a lack of agreement with HWE.

Association of allele was performed using logistic regression. Likelihood tests indicated a significant effect of ethnicity but no significant effect of age and gender. In order to avoid false discoveries due to population stratification, the association analysis was performed for each marker (SNP/haplotype) separately for each ethnic group. To come up with an overall estimate, the estimates from the different ethnic groups were combined using inverse variance weighting (IVW) technique, usually used for combining estimates in meta-analysis. The calibration and fit of the model were assessed using Hosmer Lemeshow goodness of fit and receiver operating characteristic (ROC) curves. Kruskal-Wallis test was performed for comparison of medians between three groups.

To investigate the influence of gene-gene interaction on NAFLD, Generalized Multifactor Dimensionality Reduction (GMDR) method was employed [21]. All possible interactions were tested using 10-fold cross validation with exhaustive search, which considers all possible variable combinations. GMDR provides a cross-validation consistency score which is a measure of the degree of consistency with which the selected interaction is identified as the best model among all possibilities considered. The testing balanced accuracy generated is a measure of the degree to which the interaction accurately predicts case-control status. Testing accuracy is a measure of the strength of gene-gene interaction with a power of 80% at an accuracy of 0.58–0.60, given a sample of <500 [22]. In addition, logistic regression model was performed to confirm the results of gene interaction analyses. Ethnicity was used as covariate in the gene-gene interaction analysis.

Univariate analysis of association of genotypes and histological ordinal variables was assessed using the Jonckheere-Terpstra test. Ordinal regression is perfomed for multivariate analysis of histological ordinal variables. To correct for testing for multiple histological parameters, the false discovery rate (FDR) was calculated using Benjamini-Hochberg procedure[23]. Analyses were performed using SPSS 16.0 (Chicago IL) with a two sided p< 0.05 considered to be statistically significant. Linkage disequilibrium and haplotype analyses for the five SNPs were performed using Haploview 4.2 program. The odds ratio for the haploypes was calculated using R software version 2.11.1.

We estimated that a sample size of 117 cases and 161 controls would provide 80% power at α of 0.05 with the following assumptions: the allele frequency range from 0.31-0.41, the baseline risk for the Malaysian population was 0.17, and the minimum detectable odds ratio was 2.0.

Results

Table 1 shows the demographic and clinical data of the NAFLD case (n = 144) and the control (n = 198) groups. There were 59 (41%) Malays, 54 (38%) Chinese and 31 (21%) Indians in the NAFLD group. Out of the 198 controls, 80 (41%) were Malays, 54 (27%) Chinese and 64 (32%) Indians. The differences in the parameters of BMI, HbA1c, LDL cholesterol, total cholesterol, triglycerides, AST, ALT and GGT between the two groups reflect the nature of the selection of the cases and controls.

[Figure omitted. See PDF.]

Table 1. Demographic and clinical data of the subjects.

https://doi.org/10.1371/journal.pone.0058538.t001

Table 2 further describes the demographic and clinical data of the NAFLD subjects grouped into simple steatosis (n = 33) and NASH (n = 111). The BMI, HbA1c, waist circumference, triglycerides, systolic, and diastolic blood pressure were significantly higher in the NASH group (p<0.05) as compared to the simple steatosis group.

[Figure omitted. See PDF.]

Table 2. Demographic and clinical data of the NAFLD patients.

https://doi.org/10.1371/journal.pone.0058538.t002

Genotypes and allele frequencies of AGTR1 polymorphisms

The genotypes of each SNP were in Hardy-Weinberg equilibrium for both NAFLD cases and controls, for pooled subjects as well as after stratification by ethnicity. The association tests between NAFLD and control subjects are shown in Table 3. None of the SNPs were associated with susceptibility to NAFLD. However, after ethnic stratification, in the Indian ethnic subgroup, the rs2276736, rs3772630 and rs3772627 were found to be protective against NAFLD (OR 0.40, 95% CI 0.20–0.81, p = 0.010; OR 0.43, 95% CI 0.22–0.86, p = 0.016; and OR 0.46, 95% CI 0.23–0.91, p = 0.026, respectively). We showed the sample power for the significant findings to be > 81%.

[Figure omitted. See PDF.]

Table 3. Association tests of AGTR1 SNPs and NAFLD.

https://doi.org/10.1371/journal.pone.0058538.t003

We then divided the NAFLD group into simple steatosis and NASH but did not find any association between the five SNPs with susceptibility to NASH in the pooled subjects (Table 3). Results by ethnicity group show that in the Indian ethnic subgroup, the three SNPs (rs2276736, rs3772630 and rs3772627) were also protective against NASH (OR 0.42, 95% CI 0.21–0.86, p = 0.017; OR 0.46, 95% CI 0.22–0.92, p = 0.029; and OR 0.49, 95% CI 0.24–0.98, p = 0.045, respectively).

Gene-gene interaction analysis

The overall results, pooling all subjects, indicate a lack of evidence for an association betweenAGTR1 gene and NAFLD, which is in contrast with the positive findings by Yoneda et. al [15]. Our previous study showed a strong association of PNPLA3 gene with the occurence of NAFLD [17], and we were therefore curious to find out whether there is any gene-gene interaction between PNPLA3 and AGTR1. The results with the two-locus and three-locus models are shown in Table 4. The two-locus model has the highest testing accuracy while the three-locus model has the best cross-validation consistency. Investigation of the gene-gene interaction on the occurence of NAFLD, identified a significant effect of PNPLA3 and AGTR1 interaction (empirical p = 0.007). The best GMDR model suggested for the interaction was the three-locus model. The significant interaction was confirmed by logistic regression (permuted p = 0.017).

[Figure omitted. See PDF.]

Table 4. Best fitted gene-gene interaction model.

https://doi.org/10.1371/journal.pone.0058538.t004

Genotype association in additive model revealed that the risk genotypes (AGTR1 rs3772627TT, AGTR1 rs3772630AA and PNPLA3 rs738409GG) versus one risk factor (AGTR1 rs3772627TC, AGTR1 rs3772630AG and PNPLA3 rs738409CG) versus the protective genotypes (AGTR1 rs3772627CC, AGTR1 rs3772630GG and PNPLA3 rs738409CC), conferred enhanced risk for NAFLD (OR 2.23, 95% CI 1.16–4.31, p = 0.017). Since the protective effect was seen in the Indian ethnic subgroup, we additionally tested the association specifically in this group. The risk for NAFLD in this group was dramatically enhanced (OR 12.41, 95% CI 1.61–55.42), p = 0.016) following genotypes combinations. We showed that the gene-gene interaction with an effect size of 4.5 would provide a power of > 79%. We further evaluated the association of AGTR1 SNPs with NAFLD following stratification by PNPLA3 genotypes. The results indicated the absence of significant association (data not shown).

AGTR1 polymorphisms and liver histology

Of the five SNPs, rs3772622 was the only SNP to show association with liver histology in the NAFLD patients (Table 5). The GG genotype of this SNP is significantly associated with increased in fibrosis (p = 0.009) (Fig. 1). It has to be noted that the rs3772622 is not among the SNPs that have been found to be protective. To investigate whether the G allele of the rs3772622 is associated with risk of fibrosis, we performed univariate and multivariate tests. The G allele of this SNP was shown to be associated with presence of fibrosis (p = 0.003). Subjects with the G allele has a propensity of 2.18 times to develop fibrosis score ≥ 2. The association remains significant (p = 0.003) after adjustment of other histological attributes (Table 5). We also investigated the remaining SNPs and their association with histological features of NAFLD, however none of them showed a significant association (data not shown).

[Figure omitted. See PDF.]

Figure 1. Mean fibrosis score among the genotypes of rs3772622.

p value calculated from comparison between three groups using Kruskal-Wallis test.

https://doi.org/10.1371/journal.pone.0058538.g001

[Figure omitted. See PDF.]

Table 5. Association of G allele of rs3772622 with histological features in NAFLD patients.

https://doi.org/10.1371/journal.pone.0058538.t005

Linkage disequilibrium and haplotype analysis of AGTR1 SNPs

Linkage disequilibrium (LD) analysis revealed a strong LD between the five SNPs in overall subjects and after ethnic stratification except between rs3772633 and rs3772622. Variants rs2276736 and rs3772630 are in complete LD in the Chinese and Malays, however variants rs2276736 and rs3772627 are in complete LD in the Chinese only (Table 6).

[Figure omitted. See PDF.]

Table 6. Linkage disequilibrium coefficient among five SNPs of AGTR1 gene.

https://doi.org/10.1371/journal.pone.0058538.t006

Four haplotypes with frequencies of above 5% are presented in Table 7. In pooled subjects, haplotypes ATATG, GCGCA and ATATA are more frequent than the others (35.1%, 21.8% and 20.9%, respectively). Similar results are observed in the Malays (37.8%, 22.7% and 23.2%, respectively). Meanwhile in the Chinese, the most frequent haplotypes are ATATG (34.8%) and GCGCA (25.9%). The Indians on the other hand has high frequency of haplotypes ATATG (31.4%), ATATA (22.3%) and ACGCA (24.7%). Each of the following haplotypes, ATATG, GCGCA, ATATA and ACGCA is significantly different among the ethnic groups (p< 0.0001, for each haplotype respectively). Haplotype analysis revealed that haplotype ACGCA was protective against NAFLD in the Indians (OR 0.42, 95% CI 0.15-0.95, p = 0.03). Since rs2276736, rs3772630 and rs3772627 are presented in similar structure and some were in complete LD, we further analysed the haplotype among the three SNPs. We found that the Indians presenting with the CGC haplotype were protected (OR 0.41, 95% CI 0.21–0.84, p = 0.015) against NAFLD (data not shown).

[Figure omitted. See PDF.]

Table 7. Haplotype frequencies of AGTR1 SNPs.

https://doi.org/10.1371/journal.pone.0058538.t007

Discussion

Angiotensin II appears to have an important role in the progression to non-alcoholic steatohepatitis (NASH) [13]–[14]. Angiotensin II acts on angiotensin II type 1 receptor to activate hepatic stellate cells (HSCs). The activated HSCs will express the transforming growth factor-β 1 (TGF-β 1), a profibrogenic cytokine, as a counter response to liver injury [12]–[13]. The preventive role of AGTR1 blocker in the progression to NASH was succesfully observed in the animaland human studies [24]–[27]. The deletion of angiotensin II type 1 receptor in animal studies has been shown to reduce hepatic steatosis suggesting AGTR1 is an important regulator of hepatic steatosis [16].

In the present study, our data indicated that the five previously reported SNPs of the angiotensin II type I receptor (AGTR1) gene (rs3772622, rs3772627, rs3772630, rs3772633, and rs2276736) were not associated with susceptibility to non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH). This is in contrast with the Japanese study which showed that these AGTR1 polymorphisms were significantly associated with the occurence of NAFLD and NASH [15].

Among the SNPs studied, rs2276736, rs3772630 and rs3772627 were found to be protective against NAFLD and NASH in the Indian ethnic subgroup. Interestingly, these three SNPs share the same linkage disequilibrium structure. We revealed that haplotype ACGCA is protective against NAFLD in the Indians. The study by Yoneda et al reported that the protective haplotype in the Japanese was ATATG [15]. The discrepancy in the findings may be explained by ethnic differences between the populations studied. The three major ethnic subgroups in Malaysia, namely the Malay, Chinese and Indian are of the south Asian origin with the Chinese are of the Hans descendant of Southern China [28], Indians from Southern India [29], whereas the Malays are the indigenous ethnic group. There is a supportive evidence that the Japanese are from the north Asian (Northern China and Korea) origin but not the south Asian (Southern China and Southeast Asia) based on the classical marker polymorphisms, Y chormosome and mitochondrial DNA [30]–[31]. The difference in the association studies between the two population probably reflects the difference in the genetic pools of the Malaysian and the Japanese.

The results from our study indicated that there was no difference in levels of liver enzymes between patients with simple steatosis and those with NASH. This is in agreement with the results from Sorrentino et al [32] and Mofrad et al [33], but is in contrast with the findings of several other studies which showed differences in levels between the two groups [34]–[37]. The differences in these findings is probably due to the differences in the recruitment of subjects for the various studies, in which inclusion criteria ranges from morbid obese NAFLD to adults with NAFLD. Furthermore, we observed a significantly higher triglycerides levels in patients with NASH compared to those with simple steatosis. Triglyceride storage in hepatocytesmarks the intensity of exposure to potentially toxic fatty acids. As NAFLD progresses, excess free fatty acids are compensated with increase in triglycerides production [38].

Analysis of interaction between PNPLA3 and AGTR1 gene, revealed a strong interaction between AGTR1 (rs3772627), AGTRI (rs3772630) and PNPLA3 (rs738409) SNPs on NAFLD susceptibility. Compared to the Malays and Chinese, the Indians presented with the greatest susceptibilty to NAFLD in the PNPLA3 study [17]. Thus the PNPLA3can possibly mask the effect of AGTR1 on NAFLD among the Indians. Although both genes encode different products (PNPLA3 encodes adiponutrin while AGTR1 encodes angiotensin II type 1 receptor), the strong interaction observed suggested a possible interaction as both genes are involved in lipid accumulation. The action of angiotensin II on angiotensin II type 1 receptor (AGTR1) promotes lipid accumulation in the liver whilelack of AGTR1 reduces hepatic triglycerides [16]. Adiponutrin (PNPLA3) on the other hand is highly upregulated in response to feeding and downregulated in fasting state, displaying its role in lipid storage in the liver [39]. Of note, AGTR1 and PNPLA3 are correlatively upregulated in a state of excess lipid.

We further showed that the SNP rs3772622 was associated with occurence of fibrosis but not with the other histological parameters of NAFLD. The G allele of this SNP is associated with increased fibrosis score. The Japanese also found that SNP rs3772622 was associated with the presence of fibrosis [15], but in contrast to our findings, the G allele was found to be associated with decreased fibrosis score. The difference is probably explained by a higher minor allele frequency in our population than that of the Japanese. A greater systolic blood pressure mean value (136.00 mmHg) observed in the GG genotype compared to the AA genotype (128.65 mmHg) of the SNP rs3772633 could also contribute to the discrepancy. Data of the previous report showed that hypertensive patients have significantly higher fibrosis score than that of non-hypertensive [40]. The strong linkage disequilibrium observed between the five SNPs in the overall subjects and even after ethnic stratification indicates a potential tagSNP selection based on the pairwise LD. This approach will be useful in future in order to minimize the genotyping cost yet presenting the information that could reflect the role of the other relevant SNPs within the gene.

This present study has several limitations and strengths. Studies on NAFLD spectrum have been limited since this requires valid ascertainment and scoring of the different NAFLD stages, for which liver biopsy is the gold standard in providing such accurate diagnosis. Such biopsies are limited to situations when there is a definite clinical indication for the biopsy, and this reflects the relatively small sample size in our study. Unlike the diagnosis of NASH, diagnosis of fatty liver can be made based on clinical, radiological and laboratory findings, which therefore makes it slightly easier to carry out studies on NAFLD. Ethics consideration does not allow for biopsy in control subjects, thus no biopsies were performed in controls. In order to reduce likelihood of misclassifying the controls, a strict criteria are required. These criteria includes normal BMI, fasting plasma glucose, lipid profile and liver enzymes. The difference in the distribution of these parameters between the two extremes of the NAFLD disease spectrum, provides confidence in the selection process of NAFLD cases and controls. Mild steatosis among the controls cannot however be ruled out and this may result in a bias towards the null. The resulting association between AGTR1 gene and NAFLD was smaller than hypothesized (OR less than 2.0) and would have resulted in a loss of power for the relevant analysis. The small sample size after stratification into simple steatosis and NASH is another limitation for this study. A major strength of this study was the ability to compare the association between AGTR1 polymorphisms and NAFLD among three major ethnic groups in Asia. This study suggest that the rs2276736, rs3772630 and rs3772627 are protective against NAFLD among the Indians. This study is also the first to report the interaction of gene between the AGTR1 and the commonly reported PNPLA3 gene on NAFLD susceptibility.

In conclusion, the present study suggests that the AGTR1 variants rs2276736, rs3772630 and rs3772627 are protective against NAFLD among the Indians. The G allele of rs3772622 is significantly associated with presence of fibrosis. Among the Malays and Chinese, we failed to find any significant association between the five SNPs (rs3772622, rs3772627, rs3772630, rs3772633, and rs2276736) of AGTR1 gene and susceptibility to NAFLD. We report novel finding for the gene-gene interaction between AGTR1 and the much-studied PNPLA3 genes. A larger sample size and more ethnically diverse samples are needed to replicate these findings.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Association tests of AGTR1 SNPs and NAFLD.

https://doi.org/10.1371/journal.pone.0058538.s001

(DOCX)

Acknowledgments

We gratefully acknowledge the subjects for their participation in this study and to the staffs of UMMC for their assistance in subject recruitment.

Author Contributions

Conceived and designed the experiments: SMZ RM ZM. Performed the experiments: SMZ SM PLC RCB. Analyzed the data: SMZ SR AS. Contributed reagents/materials/analysis tools: ZM RM SM PLC. Wrote the paper: SMZ.

Citation: Zain SM, Mohamed Z, Mahadeva S, Rampal S, Basu RC, Cheah P-L, et al. (2013) Susceptibility and Gene Interaction Study of the Angiotensin II Type 1 Receptor (AGTR1) Gene Polymorphisms with Non-Alcoholic Fatty Liver Disease in a Multi-Ethnic Population. PLoS ONE 8(3): e58538. https://doi.org/10.1371/journal.pone.0058538

References

1. Ludwig J, Viggiano TR, McGill DB, Oh BJ (1980) Nonalcoholic steatohepatitis: Mayo Clinic experiences with a hitherto unnamed disease. Mayo Clin Proc 55: 434–438.

2. Charlton M (2004) Nonalcoholic fatty liver disease: a review of current understanding and future impact. Clin Gastroenterol Hepatol 2: 1048–1058.

3. Hubscher SG (2006) Histological assessment of non-alcoholic fatty liver disease. Histopathology 49: 450–465.

4. Adams LA, Waters OR, Knuiman MW, Elliott RR, Olynyk JK (2009) NAFLD as a risk factor for the development of diabetes and the metabolic syndrome: an eleven-year follow-up study. Am J Gastroenterol 104: 861–867.

5. Bertrand JB, Lefevre H, Prevot S, Perlemuter G (2011) Nonalcoholic fatty liver disease in a severely obese adolescent. An arguable liver biopsy. Arch Pediatr 18: 62–65.

6. Brookes MJ, Cooper BT (2007) Hypertension and fatty liver: guilty by association? J Hum Hypertens 21: 264–270.

7. Fatani S, Itua I, Clark P, Wong C, Naderali EK (2011) The effects of diet-induced obesity on hepatocyte insulin signaling pathways and induction of non-alcoholic liver damage. Int J Gen Med 4: 211–219.

8. Vanni E, Bugianesi E, Kotronen A, De Minicis S, Yki-Jarvinen H, et al. (2010) From the metabolic syndrome to NAFLD or vice versa? Dig Liver Dis 42: 320–330.

9. Day CP (2010) Genetic and environmental susceptibility to non-alcoholic fatty liver disease. Dig Dis 28: 255–260.

10. Merriman RB, Aouizerat BE, Bass NM (2006) Genetic influences in nonalcoholic fatty liver disease. J Clin Gastroenterol 40 Suppl 1S30–33.

11. Schwimmer JB, Celedon MA, Lavine JE, Salem R, Campbell N, et al. (2009) Heritability of nonalcoholic fatty liver disease. Gastroenterology 136: 1585–1592.

12. Li X, Meng Y, Wu P, Zhang Z, Yang X (2007) Angiotensin II and Aldosterone stimulating NF-kappaB and AP-1 activation in hepatic fibrosis of rat. Regul Pept 138: 15–25.

13. Yoshiji H, Kuriyama S, Yoshii J, Ikenaka Y, Noguchi R, et al. (2001) Angiotensin-II type 1 receptor interaction is a major regulator for liver fibrosis development in rats. Hepatology 34: 745–750.

14. Oakley F, Teoh V, Ching ASG, Bataller R, Colmenero J, et al. (2009) Angiotensin II activates I kappaB kinase phosphorylation of RelA at Ser 536 to promote myofibroblast survival and liver fibrosis. Gastroenterology 136: 2334–2344.

15. Yoneda M, Hotta K, Nozaki Y, Endo H, Uchiyama T, et al. (2009) Association between angiotensin II type 1 receptor polymorphisms and the occurrence of nonalcoholic fatty liver disease. Liver Int 29: 1078–1085.

16. Nabeshima Y, Tazuma S, Kanno K, Hyogo H, Chayama K (2009) Deletion of angiotensin II type I receptor reduces hepatic steatosis. J Hepatol 50: 1226–1235.

17. Zain SM, Mohamed R, Mahadeva S, Cheah PL, Rampal S, et al.. (2012) A multi-ethnic study of a PNPLA3 gene variant and its association with disease severity in non-alcoholic fatty liver disease. Hum Genet; doi:10.1007/s00439-012-1141-y

18. Ruhl CE, Everhart JE (2005) Joint effects of body weight and alcohol on elevated serum alanine aminotransferase in the United States population. Clin Gastroenterol Hepatol 3: 1260–1268.

19. Brunt EM, Kleiner DE, Wilson LA, Belt P, Neuschwander-Tetri BA (2011) Nonalcoholic fatty liver disease (NAFLD) activity score and the histopathologic diagnosis in NAFLD: distinct clinicopathologic meanings. Hepatology 53: 810–820.

20. Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, et al. (2005) Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology 41: 1313–1321.

21. Lou XY, Chen GB, Yan L, Ma JZ, Zhu J, et al. (2007) A generalized combinatorial approach for detecting gene-by-gene and gene-by-environment interactions with application to nicotine dependence. Am J Hum Genet 80: 1125–1137.

22. Chen GB, Xu Y, Xu HM, Li MD, Zhu J, et al. (2011) Practical and theoretical considerations in study design for detecting gene-gene interactions using MDR and GMDR approaches. PLoS One 6: e16981.

23. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol 57: 12.

24. Fujita K, Yoneda M, Wada K, Mawatari H, Takahashi H, et al. (2007) Telmisartan, an angiotensin II type 1 receptor blocker, controls progress of nonalcoholic steatohepatitis in rats. Dig Dis Sci 52: 3455–3464.

25. Hirose A, Ono M, Saibara T, Nozaki Y, Masuda K, et al. (2007) Angiotensin II type 1 receptor blocker inhibits fibrosis in rat nonalcoholic steatohepatitis. Hepatology 45: 1375–1381.

26. Moreno M, Gonzalo T, Kok RJ, Sancho-Bru P, van Beuge M, et al. (2010) Reduction of advanced liver fibrosis by short-term targeted delivery of an angiotensin receptor blocker to hepatic stellate cells in rats. Hepatology 51: 942–952.

27. Yokohama S, Yoneda M, Haneda M, Okamoto S, Okada M, et al. (2004) Therapeutic efficacy of an angiotensin II receptor antagonist in patients with nonalcoholic steatohepatitis. Hepatology 40: 1222–1225.

28. Hock SS (2007) Chinese Migration. The population of Malaysia. ISEAS Publishing.

29. Periasamy M (2007) Indian Migration Into Malaya and Singapore During the British Period. Biblioasia.

30. Hammer MF, Karafet TM, Park H, Omoto K, Harihara S, et al. (2006) Dual origins of the Japanese: common ground for hunter-gatherer and farmer Y chromosomes. J Hum Genet 51: 47–58.

31. Omoto K, Saitou N (1997) Genetic origins of the Japanese: a partial support for the dual structure hypothesis. Am J Phys Anthropol 102: 437–446.

32. Sorrentino P, Tarantino G, Conca P, Perrella A, Terracciano ML, et al. (2004) Silent non-alcoholic fatty liver disease-a clinical-histological study. J Hepatol 41: 751–757.

33. Mofrad P, Contos MJ, Haque M, Sargeant C, Fisher RA, et al. (2003) Clinical and histologic spectrum of nonalcoholic fatty liver disease associated with normal ALT values. Hepatology 37: 1286–1292.

34. Abrams GA, Kunde SS, Lazenby AJ, Clements RH (2004) Portal fibrosis and hepatic steatosis in morbidly obese subjects: A spectrum of nonalcoholic fatty liver disease. Hepatology 40: 475–483.

35. Dixon JB, Bhathal PS, O'Brien PE (2001) Nonalcoholic fatty liver disease: predictors of nonalcoholic steatohepatitis and liver fibrosis in the severely obese. Gastroenterology 121: 91–100.

36. Lima ML, Mourao SC, Diniz MT, Leite VH (2005) Hepatic histopathology of patients with morbid obesity submitted to gastric bypass. Obes Surg 15: 661–669.

37. Ong JP, Elariny H, Collantes R, Younoszai A, Chandhoke V, et al. (2005) Predictors of nonalcoholic steatohepatitis and advanced fibrosis in morbidly obese patients. Obes Surg 15: 310–315.

38. Jou J, Choi SS, Diehl AM (2008) Mechanisms of disease progression in nonalcoholic fatty liver disease. Semin Liver Dis 28: 370–379.

39. Hoekstra M, Li Z, Kruijt JK, Van Eck M, Van Berkel TJ, et al. (2010) The expression level of non-alcoholic fatty liver disease-related gene PNPLA3 in hepatocytes is highly influenced by hepatic lipid status. J Hepatol 52: 244–251.

40. Corey KE, Shah N, Misdraji J, Abu Dayyeh BK, Zheng H, et al. (2009) The effect of angiotensin-blocking agents on liver fibrosis in patients with hepatitis C. Liver Int. 29: 748–753.