Chronic Kidney Disease: Combined Effects of Gene Polymorphisms of Tissue Inhibitors of Metalloproteinase 3, Total Urinary Arsenic, and Blood Lead Concentration

Abstract

The tissue inhibitor of metalloproteinase 3 (TIMP3) is known to be an anti-fibrotic factor. Arsenic, lead, and cadmium exposure and selenium intake may affect TIMP3 expression. The downregulation of TIMP3 expression is related to kidney fibrosis. Genotypes of TIMP3 are related to hypertension and cardiovascular diseases. Therefore, this study explored whether TIMP3 polymorphism is associated with hypertension-related chronic kidney disease (CKD). In addition, the combined effects of TIMP3 polymorphism and total urinary arsenic, blood lead and cadmium, and plasma selenium concentrations on CKD, were investigated. This was a case-control study, with 213 CKD patients and 423 age- and sex-matched controls recruited. Polymerase chain reaction-restriction fragment length polymorphism was used to determine TIMP3 gene polymorphisms. The concentrations of urinary arsenic species, plasma selenium, and blood lead and cadmium were measured. The odds ratio (OR) of CKD in the TIMP3rs9609643 GA/AA genotype was higher than that of the GG genotype at high levels of total urinary arsenic and blood lead; the OR and 95% confidence interval (CI) were 0.57 (0.31–1.05) and 0.52 (0.30–0.93), respectively, after multivariate adjustment. High blood lead levels tended to interact with the TIMP3rs9609643 GG genotype to increase the OR of CKD, and gave the highest OR (95% CI) for CKD of 5.97 (2.60–13.67). Our study supports a possible role for the TIMP3rs9609643 risk genotype combined with high total urinary arsenic or with high blood lead concentration to increase the OR of CKD.

1. Introduction

Chronic kidney disease (CKD) affects >10% of the world’s population and has emerged as one of leading causes of mortality worldwide [1]. Using an estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73 m² to define CKD, the prevalence of CKD in Taiwan was 11.9%, of which only 3.5% of patients were aware of their disease [2]. The incidence of end stage renal disease in Taiwan ranks first in the world [3]; therefore, exploring the etiology of CKD is an important issue in Taiwan.

Our recent study found that high plasma selenium concentrations significantly increased eGFR and decreased the odds ratio (OR) for CKD, and blood cadmium and lead concentrations, and total urinary arsenic concentration significantly decreased eGFR and increased OR for CKD [4]. Exposure to arsenic, lead and cadmium can cause tubular degeneration, fibrosis, hemorrhage and vacuolation in rat kidney tissue [5]. Studies have also found that arsenic, lead and cadmium can induce oxidative stress and cause nephrotoxicity [6,7,8]. Selenium can reduce the oxidative stress and fibrosis caused by these metals [5,9]. However, the mechanisms by which blood cadmium and lead, total urinary arsenic, and plasma selenium concentrations are associated with CKD have not been fully elucidated.

The tissue inhibitor of metalloproteinase 3 (TIMP3) is a physiological inhibitor of matrix metalloproteinases (MMPs). A disruption of the balance between MMPs and TIMPs can alter the stability and normal function in the extracellular matrix (ECM) and lead to abnormal tissue remodeling and homeostasis [10]. Of the four TIMPs, TIMP3 is the only one with an affinity for proteoglycans in the ECM [11], and it is also known to have anti-fibrotic effects [12]. The downregulation of TIMP3 may enhance the extent of tubule interstitial fibrosis (TIF) [12]. Recent studies reported that TIF was associated with CKD development and progression [13].

The increased expression of TIMP3 was observed in human kidney 2 epithelial cells under arsenic exposure [14]. Exposure to cadmium during pregnancy causes structural changes in fetal kidney tissue, which can be detected by increased levels of some kidney injury biomarkers in amniotic fluid such as albumin, osteopontin, vascular endothelial growth factor and TIMP1 [15]. One study found that with high blood lead concentration, both MMP2 and MMP9 were significantly increased, while TIMP2 was significantly decreased in the placenta of women [16]. A low-selenium diet may lead to decreased selenium content in adult rat kidneys, upregulation of MMP1 and MMP3 and downregulation of their inhibitors (TIMP1 and TIMP3), resulting in renal ultrastructural and ECM damage [17]. Exposure to arsenic, lead and cadmium may be positively or negatively associated with TIMP3, whereas selenium appears to be positively associated with TIMP3. However, results of current research are inconsistent.

The TIMP3 genes are located on chromosome 22q12.1 [18] and TIMP3 is a 24-kDa secreted protein that binds strongly to the ECM. A study of Chinese Han people found that TIMP3rs9619311 TC+CC or TIMP3rs2234921 AG+GG genotypes had a significantly higher risk of carotid plaque than those with the TT or AA genotypes, respectively [19]. A recent study reported that the TIMP3rs9619311 TT genotype had a significantly higher risk of essential hypertension than the TC+CC genotype [20]. One study demonstrated that TIMP3rs9609643 and TIMP3rs8136803 affect individual differences in breast cancer susceptibility and survival [21]. A study found that a significantly higher risk of colorectal cancer for TIMP3rs715521 AG+AA than GG genotype [22]. Current studies found that TIMP3 gene polymorphisms were associated with carotid plaques, hypertension, and cancer. Whether TIMP3 polymorphisms are associated with hypertension-related CKD remains to be explored. This study explored the association between TIMP3 genotypes and CKD. In addition, the combined effects of TIMP3 genotype and arsenic, lead or cadmium body burden and plasma selenium concentrations on CKD were evaluated.

2. Materials and Methods

2.1. Study Subjects

This was a hospital-based case-control study. From September 2005 to September 2011, 214 CKD patients and 423 age- and sex-matched healthy controls were recruited at Taipei Medical University Hospital and Taipei Wanfang Medical Center [23]. This study was approved by the Institutional Review Board of Taipei Medical University (N202101029). All study subjects were interviewed by questionnaires and biological samples were collected after they provided their informed consent.

Based on blood urea nitrogen, serum creatinine, and proteinuria, the Modification of Diet in Renal Disease formula was used to calculate eGFR by nephrologists from Taipei Medical University Hospital and Taipei Wanfang Hospital to determine the different stages CKD patients: eGFR (mL/min/1.73 m²) = 186.3 × serum creatinine^−1.154 × age^−0.203 × 1.212 (if patient is black) × 0.742 (if female) [24].

2.2. Interview and Bio-Specimen Collection

Study subjects were interviewed using a structured questionnaire by well-trained interviewers. The contents of the questionnaire included sociodemographic data; lifestyle such as cigarette smoking habit and consumption of alcohol, coffee, and tea; analgesic usage; and disease history.

An EDTA-vacuum syringe was used to collect 5–8 mL of blood, and the buffy coats were separated for DNA extraction and analysis of TIMP3rs9619311, TIMP3rs11547635, TIMP3rs715572, TIMP3rs9609643, TIMP3rs8136803, and TIMP3rs2234921 genotypes. Red blood cells were separated for measurement of lead and cadmium concentrations, and plasma was separated for measurement of selenium concentrations.

2.3. Arsenic, Cadmium, Lead, and Selenium Measurement

To ensure absence of arsenobetaine or arsenocholine (less toxic than inorganic arsenic and its methylated metabolites), high performance liquid chromatography was used to separate urinary arsenic species: arsenite (As^III), arsenate (As^V), and its metabolites, monomethylarsonic acid (MMA^V) and dimethylarsinic acid (DMA^V). Concentration of arsenic species was determined by hydride generator linked with atomic absorption spectrometry [25]. Plasma selenium and blood lead and cadmium concentrations were analyzed by inductively coupled plasma mass spectrometry [4]. If the experimental value was lower than the detection limit, the data analysis was carried out at the half-of-detection-limit concentration. The determination method, detection limit, reliability, and validity are shown in Supplementary Table S1. The sum of As^III, As^V, MMA^V, and DMA^V concentrations was termed the total urinary arsenic concentration.

2.4. Determination of TIMP3 Gene Polymorphisms

Genomic DNA was extracted by digestion with proteinase K followed by phenol and chloroform. The Agena Bioscience MassARRAY System was used according to the manufacturer’s instructions to determine the TIMP3rs9619311, TIMP3rs11547635, TIMP3rs715572, TIMP3rs9609643, TIMP3rs8136803 and TIMP3rs2234921 genotypes.

2.5. Statistical Analysis

Continuous variables are presented as mean ± standard deviation, while categorical variables are presented as frequencies (percentages). The chi-square test was used to analyze the distribution of categorical variables among the groups of the subjects, and to test whether the TIMP3 genotypes of the control group fitted a Hardy-Weinberg equilibrium. The Wilcoxon rank-sum test was conducted to compare the continuous variables between CKD cases and controls. Multiple logistic regression was used to evaluate the associations between TIMP3 genotypes and CKD by estimating OR and 95% confidence interval (CI). All models were adjusted for confounders including age, sex, and educational level; consumption of alcohol, coffee, and tea; analgesic usage; and disease histories of diabetes and hypertension. All data were analyzed using SAS 9.4 software (SAS Institute, Cary, NC, USA). A two-sided p-value < 0.05 was considered significant.

3. Results

Table 1 shows sociodemographic characteristics, lifestyle, and disease history for CKD cases and controls. The percentage of educational level above high school was significantly higher in the controls than in CKD cases. There was no difference in the proportion of cigarette smoking between the two groups. The proportion of frequent or occasional consumption of alcohol, tea and coffee was significantly higher in controls than in CKD cases. The proportion of CKD patients who routinely used analgesics was higher than in controls. Significantly more CKD patients had hypertension and diabetes than controls.

Table 2 presents the association between TIMP3 polymorphisms and CKD. The gene polymorphisms of TIMP3rs9619311, TIMP3rs11547635, TIMP3rs715572, TIMP3rs9609643, TIMP3rs8136803 and TIMP3rs2234921 were not associated with CKD.

Table 3 compares the total urinary arsenic, blood cadmium and lead, and plasma selenium concentrations between CKD and control groups. The total urinary arsenic and blood cadmium and lead concentrations were significantly higher, while plasma selenium levels were significantly lower in CKD cases than in controls.

We analyzed the association between TIMP3 genotypes and CKD after stratifying for total urinary arsenic, blood cadmium and lead, and plasma selenium concentrations (Table 4). The median values of total urinary arsenic, blood cadmium and lead, and plasma selenium concentrations in the control group were used as cut-off points for stratification analysis. The OR of CKD in the TIMP3rs9609643 GA/AA genotype was significantly lower than that in GG genotype at high total urinary arsenic and blood lead concentrations, but not at low total urinary arsenic and blood lead concentrations. The OR of CKD in TIMP3rs8136803 GT/TT genotype was lower than that in GG genotype under low total urinary arsenic, and the OR of CKD in the TIMP3rs8136803GT/TT genotype was higher than in the GG genotype under low blood lead concentration. Thus, TIMP3rs9609643 and TIMP3rs8136803 may interact with total urinary arsenic or blood lead concentrations. Therefore, the combined effects of TIMP3rs9609643 and TIMP3rs8136803, and total urinary arsenic or and blood lead concentrations, on CKD were subsequently analyzed. However, under blood cadmium and plasma selenium stratification, no association between TIMP3 polymorphisms and CKD was observed.

Pairwise analysis of combined effects of high total urinary arsenic or blood lead levels and TIMP3 risk genotype is shown in Table 5. The OR of CKD was significantly increased in dose-response with no risk factor, one risk factor, or both risk factors. We observed 3.75-fold increased odds (95% CI 1.89–7.45) of CKD cases carrying the TIMP3rs9609643 GG genotype and high levels of blood lead (>37.44 μg/L) compared to controls. The p-value for the interaction term of TIMP3rs9609643 and blood lead concentration was 0.027, and it appeared that TIMP3rs9609643 had a multiplicative interaction with blood lead on CKD; however, the significance disappeared with the multivariate adjustment. Other interactions were not significant.

4. Discussion

We found the distributions of TIMP3rs9619311, TIMP3rs11547635, TIMP3rs715572, TIMP3rs9609643, TIMP3rs8136803 and TIMP3rs2234921 between CKD cases and controls did not differ. However, subjects with the TIMP3rs9609643 GA/AA genotype had a marginally or significantly lower OR of CKD than those with GG genotype at high total urinary arsenic or high blood lead concentration. The combined effect of TIMP3rs9609643 or TIMP3rs8136803 risk genotypes and high urinary total arsenic or high blood lead concentrations gradually increased the OR for CKD with increasing risk factors.

The pathogenesis of CKD is complex and may be caused by a combination of environmental, genetic, and other factors. According to our study, the TIMP3rs9609643 GA/AA genotype was significantly associated with lower risk of CKD compared with GG genotype at high total urinary arsenic or at high blood lead concentration. There are few studies on TIMP3rs9609643; of three studies in China, one explored the correlation between TIMP3rs9609643 and thoracic aortic dissection [26], one explored its relationship with osteoarthritis [27], and one explored the correlation with high myopia [28], but all had negative results. However, one study found that women with the TIMP3rs9609643 AA genotype were 60% less likely to develop breast cancer than women with the GG genotype (OR 0.4, 95% CI 0.2–1.0) [21].

However, we observed that the TIMP3rs8136803 GT/TT genotype decreased the risk of CKD at low total arsenic levels while increasing the risk of CKD at low blood lead concentrations, compared with the GG genotype. A current study found that the genotype and allele of TIMP3rs8136803 significantly differed between with and without primary open-angle glaucoma; the frequency of TIMP3rs8136803 genotype GG in primary open-angle glaucoma was higher than that in controls [29]. However, the TIMP3rs8136803 genotype was not associated with osteoarthritis [27] and one study indicated that women with the TIMP3rs8136803 TT genotype were five times more likely to develop breast cancer than those with the GG genotype (OR 5.1, 95% CI 1.1–24.3) [21]. In addition, breast cancer cases with TIMP3rs8136803 TT were almost four times more likely to have reduced disease-free survival and there was a trend toward reduced overall survival compared to the GG genotype [21].

A Chinese study found that TIMP3rs2234921 and TIMP3rs9619311 were associated with mixture plaque [19]. Another study pointed out that TIMP3rs9619311 is related to essential hypertension [20]. Studies found that TIMP3rs9619311 was associated with hepatocellular carcinoma [30] and colorectal cancer [31]. A recent study in Taiwan found that TIMP3rs9619311 was associated with survival in cervical cancer [32]. Additionally, TIMP3rs715572 was associated with colorectal cancer [22] and survival of adenocarcinoma of the gastroesophageal junction [33]. However, TIMP3rs9619311, TIMP3rs11547635, TIMP3rs715572 and TIMP3rs2234921 were not associated with CKD in our study, and there are also some studies with similar results to ours [34,35,36,37]. At present, there are few studies on the relationship between TIMP3 genotype and disease and they have inconsistent results, so further investigation is needed. The functional relevance of these polymorphisms is also unclear. They may directly affect the expression or activity of TIMP3, or be markers for other functionally relevant variants, which requires further investigation.

Significantly gradually increased ORs for CKD with increasing risk factors (high total urinary arsenic concentration, high blood lead concentration, and TIMP3rs9609643 and TIMP3rs8136803 risk genotype) were observed in this study. This may be because exposure to arsenic, lead and cadmium induces oxidative stress and fibrosis, resulting in nephrotoxicity [6,7,8]. High concentrations of lead [16] or total urinary arsenic [38] reduced the expression of TIMP, leading to an imbalance in the MMPs/TIMPs ratio, favoring proteolytic enzyme activity and leading to generation of tissue abnormalities. In contrast, some studies found that long-term exposure to arsenic, possibly due to downregulation of TIMP3, and TIMP3 deficiency may lead to oxidative stress [39], resulting in increased renal fibrosis [12,14]. However, the role of TIMP3 in CKD remains unclear. Although the functional significance of the polymorphisms of TIMP3rs9609643 and TIMP3rs8136803 is unknown, some of the associations identified in our study support a possible role for these polymorphisms as, when combined with high total urinary arsenic or high blood lead concentration, they increased the OR of CKD.

This study had some limitations. The small number of homozygous individuals with rare alleles may have produced unstable OR estimates. Further studies with larger sample size are needed to improve the precision of point estimates when assessing TIMP3 polymorphisms and environmental metals exposure in relation to CKD. The analysis of six TIMP3 polymorphisms may not represent all the gene functions. Our study did not analyze gene polymorphisms regulating TIMP3 expression. Further studies should be conducted to assess the function of TIMP3 and its associated gene polymorphisms to determine their role in CKD development.

5. Conclusions

The risk of CKD related to high levels of blood lead or high levels of total urinary arsenic was modified by TIMP3rs9609643 GA/AA genotypes. High blood lead levels tended to interact with the TIMP3rs9609643 risk genotype to increase the risk of CKD. We recommend future studies of the levels of serum TIMP3, to determine the relevant mechanism regarding the relationships between CKD and TIMP3 polymorphisms, and environmental metals exposure.

Author Contributions

Conceptualization, formal analysis and methodology, Y.-C.L., Y.-L.H., H.-S.S., S.-L.H. and Y.-M.H., writing—original draft, Y.-C.L. and Y.-M.H., funding acquisition, Y.-M.H., writing—review and editing, Y.-M.H., visualization, Y.-M.H. and Y.-L.H. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki, and approved by approved by the Institutional Review Board of Taipei Medical University (N202101029). All study subjects were interviewed by questionnaires and biological samples were collected after they provided their informed consent.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Conflicts of Interest

The authors declare they have no actual or potential competing financial interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Tables

Table 1

Sociodemographic characteristics, lifestyle, and disease histories for the CKD cases and controls.

Variables	CKD Cases(N = 214)	Controls(N = 423)	p
Age	64.99 ± 13.54	64.33 ± 12.58	0.486 a
Gender
Male	132 (61.68)	260 (61.47)	0.958 a
Female	82 (38.32)	163 (38.53)
Educational level			<0.0001 b
Illiterate/elementary school	88 (41.12)	96 (22.70)
Junior/senior high school	71 (33.18)	146 (34.52)
College and above	55 (25.70)	181 (42.79)
Cigarette smoking			0.678 b
Non-smoker	157 (73.36)	307 (72.58)
Former smoker	32 (14.95)	73 (17.26)
Current smoker	35 (11.68)	43 (10.17)
Alcohol consumption			<0.0001 b
Never	177 (82.71)	270 (63.83)
Frequently	33 (15.42)	65 (15.37)
Occasional	4 (1.87)	88 (20.80)
Coffee consumption			<0.0001 b
Never	167 (78.04)	218 (51.05)
Frequently	28 (13.08)	102 (24.11)
Occasional	19 (8.88)	103 (24.35)
Tea consumption			<0.0001 b
Never	119 (55.61)	149 (34.89)
Frequently	69 (32.24)	168 (39.72)
Occasional	26 (12.15)	106 (25.06)
Analgesic usage			<0.0001 b
No/yes as-needed basis	188 (87.85)	404 (95.52)
Yes, routinely	26 (12.15)	19 (4.48)
Diabetes			<0.0001 b
No	129 (60.28)	379 (89.60)
Yes	85 (39.72)	44 (10.40)
Hypertension			<0.0001 b
No	92 (42.99)	296 (69.98)
Yes	122 (57.01)	127 (30.02)

Abbreviations: CKD, chronic kidney disease. Values are expressed as mean ± standard deviation or number (%) of cases and controls. ^a Wilcoxon rank-sum test. ^b χ² test.

Table 2

Associations between TIMP3 gene polymorphisms and CKD.

TIMP3 Genotypes	CKD Cases	Controls	Age-Gender AdjustedORs (95% CI)	Multivariate AdjustedORs (95% CI) a
rs11547635 C > T
CC	96 (45.50)	195 (46.54)	1.00	1.00
CT	89 (42.18)	181 (43.20)	0.99 (0.70–1.42)	1.09 (0.72–1.67)
TT	25 (12.32)	43 (10.26)	1.23 (0.71–2.12)	1.12 (0.57–2.21)
rs2234921 A > G
AA	181 (84.58)	348 (83.05)	1.00	1.00
AG	30 (14.02)	66 (15.75)	0.87 (0.55–1.39)	0.81 (0.47–1.42)
GG	3 (1.40)	5 (1.19)	1.12 (0.26–4.78)	0.92 (0.19–4.40)
rs715572 G > A
GG	89 (41.59)	179 (42.82)	1.00	1.00
GA	104 (48.60)	189 (45.22)	1.11 (0.78–1.58)	1.16 (0.76–1.76)
AA	21 (9.81)	50 (11.96)	0.85 (0.48–1.50)	0.95 (0.48–1.90)
rs9609643 G > A
GG	166 (77.57)	301 (71.67)	1.00	1.00
GA	46 (21.50)	109 (25.95)	0.77 (0.52–1.14)	0.74 (0.46–1.18)
AA	2 (0.93)	10 (2.38)	0.36 (0.08–1.66)	0.22 (0.03–1.65)
rs9619311 T > C
TT	181 (84.98)	348 (82.86)	1.00	1.00
TC	29 (13.62)	67 (15.95)	0.83 (0.52–1.32)	0.80 (0.46–1.41)
CC	3 (1.41)	5 (1.19)	1.11 (0.26–4.74)	0.92 (0.19–4.44)
rs8136803 G > T
GG	189 (88.32)	381 (90.07)	1.00	1.00
GT	23 (10.75)	37 (8.75)	1.25 (0.72–2.17)	0.88 (0.46–1.68)
TT	2 (0.93)	5 (1.18)	0.80 (0.15–4.15)	0.69 (0.12–4.12)

Abbreviations: CKD, chronic kidney disease; TIMP3, tissue inhibitor of metalloproteinase 3; OR, odds ratio; CI, confidence interval. Seven participants were missing for TIMP3rs11547635; four were missing for TIMP3rs2234921 and TIMP3rs9619311; three were missing for TIMP3rs9609643; and five were missing for TIMP3rs715572. ^a Adjusted for age, sex, educational level, analgesic usage, disease histories of diabetes and hypertension, and alcohol, coffee, and tea consumption. Multiple logistic regression models were used to calculate the association between TIMP3 genotypes and CKD.

Table 3

Total urinary arsenic, blood cadmium and lead, and plasma selenium concentrations for CKD cases and controls.

	CKD Cases (N = 214)			Controls (N = 423)
Variables	Median	First Quartile	Third Quartile	Median	First Quartile	Third Quartile	p
Total urinary arsenic(μg/g creatinine)	22.54	15.76	34.32	16.04	10.64	26.00	<0.0001
Blood cadmium (μg/L)	1.66	1.18	2.65	1.04	0.68	1.50	<0.0001
Blood lead (μg/dL)	63.65	41.44	88.15	37.44	25.56	52.64	<0.0001
Plasma selenium (μg/L)	185.68	147.80	223.45	217.85	182.60	253.15	<0.0001

p-value was tested by Wilcoxon rank-sum test.

Table 4

Associations between TIMP3 gene polymorphisms and CKD stratified by total urinary arsenic and blood lead concentrations.

	Total Urinary Arsenic > 16.04 μg/g Creatinine			Total Urinary Arsenic ≤ 16.04 μg/g Creatinine
TIMP3 Genotypes	CKD Cases/Controls	Age-Sex Adjusted ORs (95% CI)	Multivariate Adjusted ORs (95% CI) a	CKD Cases/Controls	Age-Sex Adjusted ORs (95% CI)	Multivariate Adjusted ORs (95% CI) a
rs9609643 G > A
GG	124/143	1.00	1.00	42/158	1.00	1.00
GA/AA	34/67	0.55 (0.34–0.90) *	0.57 (0.31–1.05) +	14/52	0.98 (0.49–1.95)	0.97 (0.45–2.12)
rs8136803 G > T
GG	137/195	1.00	1.00	52/186	1.00	1.00
GT/TT	21/16	1.86 (0.94–3.71) +	1.54 (0.67–3.54)	4/26	0.56 (0.19–1.70)	0.31 (0.09–1.09) +
	Blood lead > 37.44 μg/dL			Blood lead ≤ 37.44 μg/dL
rs9609643 G > A
GGGA/AA	134/13836/74	1.000.50 (0.31–0.79) **	1.000.52 (0.30–0.93) *	32/16212/45	1.001.36 (0.65–2.85)	1.000.89 (0.37–2.17)
rs8136803 G > T
GG	153/189	1.00	1.00	36/192	1.00	1.00
GT/TT	17/23	0.91 (0.47–1.76)	0.55 (0.24–1.26)	8/19	2.28 (0.92–5.63) +	2.11 (0.74–1.51)

Abbreviations: TIMP3, tissue inhibitor of metalloproteinase 3; OR, odds ratio; CI, confidence interval. Seven participants were missing for TIMP3rs11547635; four were missing for TIMP3rs2234921 and TIMP3rs9619311; three were missing for TIMP3rs9609643; and five were missing for TIMP3rs715572. ^a Adjusted for age, sex, educational level, analgesic usage, disease histories of diabetes and hypertension, and alcohol, coffee, and tea consumption. ⁺ 0.05 < p < 0.1, ^* p < 0.05, ^** p < 0.01.

Table 5

Effects of interactions of TIMP3 gene polymorphisms with total urinary arsenic and blood lead concentrations on CKD.

Metals	TIMP3 Genotypes	CKD Cases/Controls	Age-Gender AdjustedORs (95% CI)	Multivariate AdjustedORs (95% CI) a
Total urinary arsenic (μg/g creatinine)	rs9609643 G > A
≤16.04	GA/AA	14/52	1.00 &,***	1.00 &,***
≤16.04	GG	42/158	0.98 (0.50–1.95)	1.14 (0.53–2.43)
>16.04	GA/AA	34/67	1.91 (0.93–3.93) +	1.87 (0.82–4.28)
>16.04	GG	124/143	3.29 (1.73–6.24) ***	3.10 (1.51–6.40) **
		Synergy index	2.57 (0.59–11.31)	2.08 (0.47–9.22)
		p interaction	0.197	0.761
Total urinary arsenic (μg/g creatinine)	rs8136803 G > T
≤16.04	GT/TT	4/26	1.00 &,***	1.00 &,***
≤16.04	GG	52/186	1.80 (0.60–5.38)	3.21 (0.96–10.75) +
>16.04	GG	137/195	4.58 (1.56–13.44) **	6.85 (2.08–10.75) **
>16.04	GT/TT	21/16	8.55 (2.48–29.48) ***	10.42 (2.60–41.66) ***
		Synergy index	1.72 (0.73–4.05)	1.17 (0.48–2.85)
		p interaction	0.979	0.797
Blood lead (μg/dL)	rs9609643 G > A
≤37.44	GA/AA	12/45	1.00 &,***	1.00 &,***
≤37.44	GG	32/163	0.74 (0.35–1.55)	1.06 (0.44–2.52)
>37.44	GA/AA	36/74	1.86 (0.88–3.96)	3.13 (1.27–7.72) *
>37.44	GG	134/138	3.75 (1.89–7.45) ***	5.97 (2.60–13.67) ***
		Synergy index	4.58 (0.45–47.15)	2.27 (0.87–5.91)
		p interaction	0.027	0.137
Blood lead (μg/dL)	rs8136803 G > T
≤37.44	GG	36/192	1.00 &,***	1.00 &,***
≤37.44	GT/TT	8/19	2.28 (0.93–5.61) +	1.93 (0.69–5.42)
>37.44	GT/TT	17/23	4.02 (1.95–8.30) ***	3.32 (1.41–7.82) **
>37.44	GG	153/189	4.43 (2.91–6.67) ***	5.76 (3.44–9.66) ***
		Synergy index	0.80 (0.36–1.75)	1.46 (0.53–4.04)
		p interaction	0.647	0.563

Abbreviations: TIMP3, tissue inhibitor of metalloproteinase 3; OR, odds ratio; CI, confidence interval. ^& tested for linear trend. ⁺ 0.05 < p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001. ^a Adjusted for age, sex, educational level, analgesic usage, disease histories of diabetes and hypertension, and alcohol, coffee, and tea consumption.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph20031886/s1, Table S1. The validity and reliability of urinary arsenic species, plasma selenium, and red blood cell lead and cadmium.

References

1. Kovesdy, C.P. Epidemiology of chronic kidney disease: An update 2022. Kidney Int. Suppl.; 2022; 12, pp. 7-11. [DOI: https://dx.doi.org/10.1016/j.kisu.2021.11.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35529086]

2. Wen, C.P.; Cheng, T.Y.D.; Tsai, M.K.; Chang, Y.C.; Chan, H.T.; Tsai, S.P.; Chiang, P.H.; Hsu, C.C.; Sung, P.K.; Hsu, Y.H. et al. All-cause mortality attributable to chronic kidney disease: A prospective cohort study based on 462 293 adults in Taiwan. Lancet; 2008; 371, pp. 2173-2182. [DOI: https://dx.doi.org/10.1016/S0140-6736(08)60952-6]

3. Saran, R.; Robinson, B.; Abbott, K.C.; Agodoa, L.Y.; Albertus, P.; Ayanian, J.; Balkrishnan, R.; Bragg-Gresham, J.; Cao, J.; Chen, J.L.T. et al. US Renal Data System 2016 Annual Data Report: Epidemiology of kidney disease in the United States. Am. J. Kidney Dis.; 2017; 69, (Suppl. S1), pp. S1-S688. [DOI: https://dx.doi.org/10.1053/j.ajkd.2016.12.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28236831]

4. Wu, C.-Y.; Wong, C.-S.; Chung, C.-J.; Wu, M.-Y.; Huang, Y.-L.; Ao, P.-L.; Lin, Y.-F.; Lin, Y.-C.; Shiue, H.-S.; Su, C.-T. et al. The association between plasma selenium and chronic kidney disease related to lead, cadmium and arsenic exposure in a Taiwanese population. J. Hazard. Mater; 2019; 375, pp. 224-232. [DOI: https://dx.doi.org/10.1016/j.jhazmat.2019.04.082] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31075550]

5. Riaz, M.A.; Nisa, Z.U.; Anjum, M.S.; Butt, H.; Mehmood, A.; Riaz, A.; Akhtar, A.B.T. Assessment of metals induced histopathological and gene expression changes in different organs of non-diabetic and diabetic rats. Sci. Rep.; 2020; 10, 5897. [DOI: https://dx.doi.org/10.1038/s41598-020-62807-0]

6. Gungor, H.; Kara, H. Effects of selenium, zinc, insulin and metallothionein on cadmium-induced oxidative stress and metallothionein gene expression levels in diabetic rats. J. Basic Clin. Physiol. Pharmacol.; 2020; 31, 20190198. [DOI: https://dx.doi.org/10.1515/jbcpp-2019-0198] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32114524]

7. Shi, Y.; Tian, C.; Yu, X.; Fang, Y.; Zhao, X.; Zhang, X.; Xia, D. Protective Effects of Smilax glabra Roxb. Against Lead-Induced Renal Oxidative Stress, Inflammation and Apoptosis in Weaning Rats and HEK-293 Cells. Front. Pharmacol.; 2020; 11, 556248. [DOI: https://dx.doi.org/10.3389/fphar.2020.556248]

8. Dutta, S.; Saha, S.; Mahalanobish, S.; Sadhukhan, P.; Sil, P.C. Melatonin attenuates arsenic induced nephropathy via the regulation of oxidative stress and inflammatory signaling cascades in mice. Food Chem. Toxicol.; 2018; 118, pp. 303-316. [DOI: https://dx.doi.org/10.1016/j.fct.2018.05.032]

9. Le, G.; Yang, L.; Du, H.; Hou, L.; Ge, L.; Sylia, A.; Muhmood, A.; Chen, X.; Han, B.; Huang, K. Combination of zinc and selenium alleviates ochratoxin A-induced fibrosis via blocking ROS-dependent autophagy in HK-2 cells. J. Trace Elem Med. Biol.; 2022; 69, 126881. [DOI: https://dx.doi.org/10.1016/j.jtemb.2021.126881]

10. Takawale, A.; Sakamuri, S.S.; Kassiri, Z. Extracellular matrix communication and turnover in cardiac physiology and pathology. Compr. Physiol.; 2015; 5, pp. 687-719.

11. Brew, K.; Nagase, H. The tissue inhibitors of metalloproteinases (TIMPs): An ancient family with structural and functional diversity. Biochim. Biophys Acta; 2010; 1803, pp. 55-71. [DOI: https://dx.doi.org/10.1016/j.bbamcr.2010.01.003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20080133]

12. Kassiri, Z.; Oudit, G.Y.; Kandalam, V.; Awad, A.; Wang, X.; Ziou, X.; Maeda, N.; Herzenberg, A.M.; Scholey, J.W. Loss of TIMP3 enhances interstitial nephritis and fibrosis. J. Am. Soc. Nephrol.; 2009; 20, pp. 1223-1235. [DOI: https://dx.doi.org/10.1681/ASN.2008050492] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19406980]

13. Zhou, D.; Fu, H.; Zhang, L.; Zhang, K.; Min, Y.; Xiao, L.; Lin, L.; Bastacky, S.I.; Liu, Y. Tubule-Derived Wnts Are Required for Fibroblast Activation and Kidney Fibrosis. J. Am. Soc. Nephrol.; 2017; 28, pp. 2322-2336. [DOI: https://dx.doi.org/10.1681/ASN.2016080902] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28336721]

14. Chang, Y.W.; Singh, K.P. Arsenic induces fibrogenic changes in human kidney epithelial cells potentially through epigenetic alterations in DNA methylation. J. Cell Physiol.; 2019; 234, pp. 4713-4725. [DOI: https://dx.doi.org/10.1002/jcp.27244] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30191986]

15. Jacobo-Estrada, T.; Cardenas-Gonzalez, M.; Santoyo-Sánchez, M.; Parada-Cruz, B.; Uria-Galicia, E.; Arreola-Mendoza, L.; Barbier, O. Evaluation of kidney injury biomarkers in rat amniotic fluid after gestational exposure to cadmium. J. Appl. Toxicol.; 2013; 6, pp. 1183-1193. [DOI: https://dx.doi.org/10.1002/jat.3286] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26815315]

16. Gonzalez-Puebla, E.; Gonzalez-Horta, C.; Infante-Ramirez, R.; Sanin, L.H.; Levario-Carrillo, M.; Sanchez-Ramirez, B. Altered expressions of MMP-2, MMP-9, and TIMP-2 in placentas from women exposed to lead. Hum. Exp. Toxicol.; 2012; 31, pp. 662-670. [DOI: https://dx.doi.org/10.1177/0960327111431706] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22241630]

17. Han, J.; Liang, H.; Yi, J.; Tan, W.; He, S.; Wu, X.; Shi, X.; Ma, J.; Guo, X. Selenium deficiency induced damages and altered expressions of metalloproteinases and their inhibitors (MMP1/3, TIMP1/3) in the kidneys of growing rats. J. Trace Elem. Med. Biol.; 2016; 34, pp. 1-9. [DOI: https://dx.doi.org/10.1016/j.jtemb.2015.11.003]

18. Apte, S.S.; Mattei, M.G.; Olsen, B.R. Cloning of the cDNA encoding human tissue inhibitor of metalloproteinases-3 (TIMP-3) and mapping of the TIMP3 gene to chromosome 22. Genomics; 1994; 19, pp. 86-90. [DOI: https://dx.doi.org/10.1006/geno.1994.1016]

19. Zheng, Z.; He, X.; Zhu, M.; Jin, X.; Li, C.; Zhu, F.; Lv, C.; Li, W.; Hu, X.; Wang, W. et al. Tissue inhibitor of the metalloproteinases-3 gene polymorphisms and carotid plaque susceptibility in the Han Chinese population. Int. J. Neurosci.; 2018; 128, pp. 920-927. [DOI: https://dx.doi.org/10.1080/00207454.2018.1436544]

20. Chen, C.; Yang, M.; Dou, L.P.; Ling, D.M.; Huang, S. Association of the interaction between the rs9619311 and rs402007 polymorphisms and smoking with essential hypertension in Chinese Han population. Medicine; 2020; 99, e20552. [DOI: https://dx.doi.org/10.1097/MD.0000000000020552]

21. Peterson, N.B.; Beeghly-Fadiel, A.; Gao, Y.-T.; Long, J.; Cai, Q.; Shu, X.-O.; Zheng, W. Polymorphisms in tissue inhibitors of metalloproteinases-2 and -3 and breast cancer susceptibility and survival. Int. J. Cancer; 2009; 125, pp. 844-850. [DOI: https://dx.doi.org/10.1002/ijc.24405]

22. Wang, N.; Zhou, S.; Fang, X.-C.; Gao, P.; Qiao, Q.; Wu, T.; He, X.-L. MMP-2, -3 and TIMP-2, -3 polymorphisms in colorectal cancer in a Chinese Han population: A case-control study. Gene; 2020; 730, 144320. [DOI: https://dx.doi.org/10.1016/j.gene.2019.144320] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31887330]

23. Chen, W.-J.; Huang, Y.-L.; Shiue, H.-S.; Chen, T.-W.; Lin, Y.-F.; Huang, C.-Y.; Lin, Y.-C.; Han, B.-C.; Hsueh, Y.-M. Renin-angiotensin-aldosterone system related gene polymorphisms and urinary total arsenic is related to chronic kidney disease. Toxicol. Appl. Pharmacol.; 2014; 279, pp. 95-102. [DOI: https://dx.doi.org/10.1016/j.taap.2014.05.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24907556]

24. Levey, A.S.; Bosch, J.P.; Lewis, J.B.; Greene, T.; Rogers, N.; Roth, D. A more accurate method to estimate glomerular filtration rate from serum creatinine: A new prediction equation. Modification of Diet in Renal Disease Study Group. AnnInternMed; 1999; 130, pp. 461-470.

25. Hsueh, Y.M.; Huang, Y.L.; Huang, C.C.; Wu, W.L.; Chen, H.M.; Yang, M.H.; Lue, L.C.; Chen, C.J. Urinary levels of inorganic and organic arsenic metabolites among residents in an arseniasis-hyperendemic area in Taiwan. J. Toxicol. Environ. Health A; 1998; 54, pp. 431-444. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9661909]

26. Wang, X.L.; Liu, O.; Qin, Y.W.; Zhang, H.J.; Lv, Y. Association of the polymorphisms of MMP-9 and TIMP-3 genes with thoracic aortic dissection in Chinese Han population. Acta Pharmacol. Sin.; 2014; 35, pp. 351-355. [DOI: https://dx.doi.org/10.1038/aps.2013.179]

27. Tong, Z.; Liu, Y.; Chen, B.; Yan, L.; Hao, D. Association between MMP3 and TIMP3 polymorphisms and risk of osteoarthritis. Oncotarget; 2017; 8, pp. 83563-83569. [DOI: https://dx.doi.org/10.18632/oncotarget.18745]

28. Leung, K.H.; Yiu, W.C.; Yap, M.K.; Ng, P.W.; Fung, W.Y.; Sham, P.C.; Yip, S.P. Systematic investigation of the relationship between high myopia and polymorphisms of the MMP2, TIMP2, and TIMP3 genes by a DNA pooling approach. Invest. Ophthalmol. Vis. Sci.; 2011; 52, pp. 3893-3900. [DOI: https://dx.doi.org/10.1167/iovs.11-7286]

29. Ji, M.L.; Jia, J. Correlations of TIMP2 and TIMP3 gene polymorphisms with primary open-angle glaucoma. Eur. Rev. Med. Pharmacol. Sci.; 2019; 23, pp. 5542-5547.

30. Tsai, H.-T.; Hsieh, M.-J.; Chiou, H.-L.; Lee, H.-L.; Hsin, M.-C.; Liou, Y.-S.; Yang, C.-C.; Yang, S.-F.; Kuo, W.-H. TIMP-3 -1296 T>C and TIMP-4 -55 T>C gene polymorphisms play a role in the susceptibility of hepatocellular carcinoma among women. Tumour. Biol.; 2014; 35, pp. 8999-9007. [DOI: https://dx.doi.org/10.1007/s13277-014-2170-z]

31. Banday, M.Z.; Sameer, A.S. Strong association of tissue inhibitor of metalloproteinase (TIMP)-2 and -3 promoter single nucleotide polymorphisms with risk of colorectal cancer in ethnic Kashmiri population-a case control study. Biosci. Rep.; 2019; 39, BSR20190478. [DOI: https://dx.doi.org/10.1042/BSR20190478] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30988064]

32. Chung, J.-F.; Chen, C.-L.; Nassef, Y.; Shiu, B.-H.; Wang, C.-H.; Kuo, F.-H.; Hsiao, Y.-H.; Yang, S.-F.; Wang, P.-H. Effect of tissue inhibitor of metalloproteinases-3 genetics polymorphism on clinicopathological characteristics of uterine cervical cancer patients in Taiwan. Int. J. Med. Sci.; 2022; 19, pp. 1013-1022. [DOI: https://dx.doi.org/10.7150/ijms.72378]

33. Bashash, M.; Shah, A.; Hislop, G.; Treml, M.; Bretherick, K.; Janoo-Gilani, R.; Leach, S.; Le, N.; Bajdik, C.; Brooks-Wilson, A. Genetic polymorphisms at TIMP3 are associated with survival of adenocarcinoma of the gastroesophageal junction. PLoS ONE; 2013; 8, e59157. [DOI: https://dx.doi.org/10.1371/journal.pone.0059157]

34. Chang, J.-H.; Lai, T.-C.; Yang, P.-J.; Shih, P.-C.; Yang, Y.-C.; Lee, K.-L.; Liu, T.-C.; Tsao, T.; Yang, S.-F.; Chien, M.-H. Associations of TIMP-3 Genetic Polymorphisms with EGFR Statuses and Cancer Clinicopathologic Development in Lung Adenocarcinoma Patients. Int. J. Mol. Sci.; 2020; 21, 8023. [DOI: https://dx.doi.org/10.3390/ijms21218023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33126605]

35. Song, G.; Yan, J.; Zhang, Q.; Li, G.; Chen, Z.J. Association of tissue inhibitor of metalloproteinase gene polymorphisms and unexplained recurrent spontaneous abortions in Han Chinese couples. Eur. J. Obstet. Gynecol. Reprod. Biol.; 2014; 181, pp. 84-88. [DOI: https://dx.doi.org/10.1016/j.ejogrb.2014.07.013] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25128867]

36. Su, C.W.; Huang, Y.W.; Chen, M.K.; Su, S.C.; Yang, S.F.; Lin, C.W. Polymorphisms and Plasma Levels of Tissue Inhibitor of Metalloproteinase-3: Impact on Genetic Susceptibility and Clinical Outcome of Oral Cancer. Medicine; 2015; 94, e2092. [DOI: https://dx.doi.org/10.1097/MD.0000000000002092]

37. Wu, Z.; Chen, H.; Pan, L.; Yu, W.; Lou, C.; Chen, J.; He, D. Effect of TIMP2/TIMP3 genes on the risk of osteosarcoma in Zhejiang population. Medicine; 2021; 100, e24818. [DOI: https://dx.doi.org/10.1097/MD.0000000000024818] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33725949]

38. Josyula, A.B.; Poplin, G.S.; Kurzius-Spencer, M.; McClellen, H.E.; Kopplin, M.J.; Stürup, S.; Lantz, R.C.; Burgess, J.L. Environmental arsenic exposure and sputum metalloproteinase concentrations. Environ. Res.; 2006; 102, pp. 283-290. [DOI: https://dx.doi.org/10.1016/j.envres.2006.01.003]

39. Fujii, T.; Duarte, S.; Lee, E.; Ke, B.; Busuttil, R.W.; Coito, A.J. Tissue Inhibitor of Metalloproteinase 3 Deficiency Disrupts the Hepatocyte E-Cadherin/beta-Catenin Complex and Induces Cell Death in Liver Ischemia/Reperfusion Injury. Liver Transpl.; 2020; 26, pp. 113-126. [DOI: https://dx.doi.org/10.1002/lt.25667]