A DNA adductome analysis revealed a reduction in the global level of C5-hydroxymethyl-2′-deoxycytidine in the non-tumoral upper urinary tract mucosa of urothelial carcinoma patients
Genes and Environment volume 43, Article number: 52 (2021)
DNA adducts, covalent modifications to DNA due to exposure to specific carcinogens, cause the mispairing of DNA bases, which ultimately results in DNA mutations. DNA methylation in the promoter region, another type of DNA base modification, alters the DNA transcription process, and has been implicated in carcinogenesis in humans due to the down-regulation of tumor suppressor genes. Difficulties are associated with demonstrating the existence of DNA adducts or chemically modified bases in the human urological system. Apart from aristolochic acid-DNA adducts, which cause urothelial carcinoma and endemic nephropathy in a particular geographical area (Balkan), limited information is currently available on DNA adduct profiles in renal cell carcinoma and upper urinary tract urothelial carcinoma, including renal pelvic cancer and ureteral cancer.
To elucidate the significance of DNA adducts in carcinogenesis in the urothelial system, we investigated 53 DNA adducts in the non-tumoral renal parenchyma and non-tumoral renal pelvis of patients with renal cell carcinoma, upper urinary tract urothelial carcinoma, and other diseases using liquid chromatography coupled with tandem mass spectrometry. A comparative analysis of tissue types, the status of malignancy, and clinical characteristics, including lifestyle factors, was performed.
C5-Methyl-2′-deoxycytidine, C5-hydroxymethyl-2′-deoxycytidine (5hmdC), C5-formyl-2′-deoxycytidine, 2′-deoxyinosine, C8-oxo-2′-deoxyadenosine, and C8-oxo-2′-deoxyguanosine (8-OHdG) were detected in the renal parenchyma and renal pelvis. 8-OHdG was more frequently detected in the renal pelvis than in the renal cortex and medulla (p = 0.048 and p = 0.038, respectively). 5hmdC levels were significantly lower in the renal pelvis of urothelial carcinoma patients (n = 10) than in the urothelium of patients without urothelial carcinoma (n = 15) (p = 0.010). Regarding 5hmdC levels in the renal cortex and medulla, Spearman’s rank correlation test revealed a negative correlation between age and 5hmdC levels (r = − 0.46, p = 0.018 and r = − 0.45, p = 0.042, respectively).
The present results revealed a reduction of 5hmdC levels in the non-tumoral urinary tract mucosa of patients with upper urinary tract urothelial carcinoma. Therefore, the urothelial cell epithelia of patients with upper urinary tract cancer, even in non-cancerous areas, may be predisposed to urothelial cancer.
The initiation of carcinogenesis involves the induction of irreversible genetic alterations, including mutations, transversions, transitions, and/or deletions in DNA. One of the causes of DNA mutations is DNA adducts, which are covalent modifications to DNA due to exposure to specific carcinogens. Another cause is modified DNA bases, which are enzymatically generated in the human body . Some DNA base modifications in promoter regions (adducts in a broad sense), such as C5-methyl-2′-deoxycytidine (5mdC), alter the DNA transcription process, which regulates gene expression. Promoter methylation in tumor suppressor genes has been implicated in carcinogenesis and the poor prognosis of some cancers [2,3,4,5]. DNA adducts and modified DNA bases, which can be detected in analytical chemistry using tissue, urine, and blood, have been examined as biomarkers in individual assessments of the risk of developing cancer. Therefore, the analysis of DNA adducts and modified DNA bases may eventually lead to the early detection of cancer .
DNA adducts in human tissue have been investigated as the cause of DNA mutations in the lungs and gastrointestinal tract, which are directly exposed to potential procarcinogens and mutagens derived from environmental substances, such as tobacco smoke, polluted air, alcohol, and dietary carcinogens. Thus, the DNA lesions caused by these environmental constituents have been detected in the lungs and gastrointestinal tract [7,8,9,10]. The urinary system, which is constantly exposed to urine containing hydrophilic phase II metabolites, is another organ that is exposed to many environmental mutagens. Urinary C8-oxo-2′-deoxyguanosine (8-OHdG) has been investigated as a biomarker of oxidative DNA damage . However, limited information is currently available on other DNA adducts in human urological organs, including the kidneys, ureter, and bladder [12,13,14,15,16]. Aristolochic acid (AA)-DNA adducts and their role in the development of urothelial carcinoma (UC), renal cell carcinoma (RCC), and nephropathy in endemic areas and Chinese herb users  have recently been attracting increasing attention [18,19,20,21,22]. The mechanisms underlying carcinogenesis induced by AA have been investigated in vitro and in vivo . Smoking has been identified as the most important risk factor for UC, particularly bladder cancer . DNA adducts induced by aromatic amines, such as 4-aminobiphenyl from tobacco smoke, have been detected in the non-tumor bladder tissues of patients with bladder cancer and interactions with genetic polymorphisms were reported [12, 13].
The roles of the hydroxymethylation of cytosine, detected as C5-hydroxymethyl-2′-deoxycytidine (5hmdC) or 5-hydroxymethylcytosine (5hmC), in cancer are currently being investigated. A previous study demonstrated that 5hmdC levels were lower in malignant tumor tissue, including clear cell RCC and UC, than in normal tissue . 5hmC functions as a cis-regulatory element in transcription factor-binding regions. Furthermore, it may associate with histone modifications to alter the configuration of chromatin and modulate alternative splicing . 5hmC has also been shown to maintain the integrity of stalled replication forks and has been linked to resistance to poly (ADP-ribose) polymerase inhibitors and cisplatin .
Although various types of known and unknown noncanonical deoxyribonucleosides may be involved in urinary carcinogenesis, their role remains unclear [12, 16, 25]. A more detailed understanding of DNA adducts and modified DNA bases, particularly in the human renal parenchyma and renal pelvis, will provide important insights into the etiologies and carcinogenesis of RCC and upper urinary tract urothelial carcinoma (UTUC). Therefore, our research has focused on DNA adducts and modified DNA bases in the non-tumoral tissue of urological organs.
We previously detected lipid peroxidation-induced DNA adducts in the non-tumoral tissue of gastric cancer patients . This finding prompted us to investigate whether these adducts are also present in urological organs. RCC originates from the renal cortex and UTUC from the urothelium of the renal pelvis and ureter. Therefore, we herein compared the DNA adduct profiles of the renal parenchyma obtained from patients with and without RCC and those with and without UTUC.
MATERIALS and METHODS
Sample tissues were obtained from consecutive cases of RCC, UTUC, or other benign diseases, including oncocytoma and urinary lithiasis, that underwent nephrectomy and nephroureterectomy at Hamamatsu University Hospital. The renal cortex, renal medulla, and urothelial mucosa of the renal pelvis were collected as the non-tumoral tissue of resected kidneys. Representative cases are shown in Fig. 1. Tissues were also collected from autopsy cases of duodenal ulcers, germ cell tumors, interstitial pneumonia, pleomorphic sarcoma, and lung cancer. Tissues were obtained from 49 patients (26 renal cortex tissue samples from 14 RCC patients [the cortex-RCC group] and 12 non-RCC patients [the cortex-non-RCC group], including 8 UTUC cases, 2 non-malignant cases, one ovarian cancer case, and one sarcoma case; 21 renal medulla tissue samples from 12 RCC patients [the medulla-RCC group] and 9 non-RCC patients [the medulla-non-RCC group], including 8 UTUC cases and one non-malignant case; 25 samples of the urothelium of the renal pelvis from 10 UTUC patients [the pelvis-UTUC group] and 15 non-UTUC patients [the pelvis-non-UTUC group], including 10 RCC cases and 5 cases without urological malignancies). Clinicopathological features and tissue types are shown in Table 1. In several cases, two or three portions were collected from individual patients. For example, sample No.1 from the renal cortex and sample No.29 from the renal medulla were taken from the same patient. Details on the kidney samples collected, including clinicopathological features, are shown in Supplementary Table S1. Kidney samples obtained immediately after resection or autopsy were frozen in liquid nitrogen and stored at − 80 °C in a deep freezer until DNA extraction. The locations of sampled portions in the kidneys were photographed to ensure that tumorous tissues were not contaminated in the sampled tissues. All procedures to collect tissue samples were performed according to the principles guided by the Japanese Society of Pathology .
DNA extraction from kidney tissue samples and enzymatic digestion
DNA extraction was conducted on all tissue samples using a Gentra Puregene™ Tissue Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions. Using a median tissue mass of 110 mg (range, 27–296 mg), we purified a median of 93.4 μg DNA (range, 41.3–526.3 μg). Purified DNA was quantified by measuring absorbance at 260 and 280 nm with a Nano Drop ND-1000 spectrophotometer (Thermo Fisher Scientific, Tokyo, Japan). To prevent the oxidation of DNA, deferoxamine was added to the DNA solvent at a final concentration of 0.1 mM. DNA was prepared for liquid chromatography coupled with tandem mass spectrometry (LC/MS) using an 8-OHdG Assay Preparation Reagent Set (Wako, Osaka, Japan). Ten micrograms of DNA was transferred to a 1.5-mL Eppendorf tube and diluted to 52.5 μL using distilled water. After samples had been incubated at 98 °C for 2 min and placed on ice for 5 min, 6.7 μL of Acetic Acid Buffer and 3.5 μL of Nuclease P1 were added to the mixture, which was then incubated at 37 °C for 30 min. Following the addition of 7.0 μL of Tris Buffer and 0.35 μL of Alkaline Phosphatase Solution, samples were incubated at 37 °C for 30 min. The solution was centrifuged at 15000×g at 20 °C for 15 min using Nanosep 10 K (Pall Life Science, Michigan, USA) to filter enzymatically hydrolyzed deoxyribonucleoside monomers.
DNA adduct identification and semi-quantification by LC/MS
A 4000 QTRAP mass spectrometer (SCIEX, Framingham, MA, USA) with an Acquity UPLC system (Waters, Milford, MA, USA) was used to analyze DNA adducts. An aliquot of a digested DNA sample (1.1 μg/7.5 μL) was injected and separated using an Acquity UPLC HSS T3 column (2.1 × 100 mm, Waters) with mobile phase A: 0.02% (v/v) acetic acid in water and mobile phase B: methanol at 0.2 mL/min with an isocratic and linear gradient elution program (3% B from 0 to 5 min; 3–15% B from 5 to 10 min; 15–80% B from 10 to 25 min; 80% B from 25 to 30 min). Other experimental conditions were set as follows: ion source temperature, 400 °C; ion spray voltage, 5500 V; declustering potential, 51 V; entrance potential, 15 V; collision energy, 18 eV; collision cell exit potential, 6 V. Selected reaction monitoring in the positive ion mode was designed to mainly detect the neural loss of 2′-deoxyribose from positively ionized 2′-deoxynucleoside adducts ([M + H]+ → [M + H-116]+). To enhance detection, we limited and scheduled the dwell time for counting ions at the column retention time of standard compounds according to each mass transition. We visually evaluated chromatograms and identified known DNA adducts and modifications, which are described in Supplementary Table S2. Distilled water and a “blank sample” in which distilled water was used instead of DNA in the enzymatic digestion procedure were used in every measurement to monitor the background. A “DNA sample”, consisting of 24.5% 2′-deoxyadenosine (dA), 24.5% 2′-deoxythymidine (dT), 20.5% 2′-deoxyguanosine (dG), and 20.5% 2′-deoxycytidine (dC) at 1.0 μg/mL, and a “standard sample” containing all the standard compounds we prepared (Supplementary Table S2), were used to identify each DNA adduct and monitor day-to-day variance.
Acquired wiff files were analyzed using Skyline software v184.108.40.206 (https://skyline.ms). The peak boundary of each adduct was defined as the minimum and maximum column retention times based on five measurements using the “standard sample”. Chromatograms were visually evaluated to exclude apparent peaks with a deviation of more than 0.2 min from the authentic retention time of the “standard sample” as well as narrower peaks, which were suspected to be artifacts. We set the detection limit as the peak area, which was less than the mean values of five “blank sample” measurements. We surveyed 53 types of DNA adducts and modifications and characterized the following 6 molecules: 5mdC, 5hmdC, C5-formyl-2′-deoxycytidine (5fdC), 2′-deoxyinosine (dI), C8-oxo-2′-deoxyadenosine (8-OHdA), and 8-OHdG, from kidney samples. Raw data on the areas of these 6 adducts and 4 isotopologues of deoxyribonucleosides are shown in Supplementary Table S3.
Raw area data were normalized by the quantity of the naturally-occurring isotopologues of deoxyribonucleosides detected in each sample, such as 13C2-dA, 13C2-dT, 13C2-dG, and 13C2-dC. To adjust intra-day and day-to-day differences, the areas of these four isotopologues in kidney samples were divided by the median value of the corresponding isotopologue in all measurements. The median of the ratios of these four isotopologues in a sample was defined as the “correction coefficient” of each sample. In our experiments, the “correction coefficient” ranged between 0.43 and 1.60. To normalize the area data of kidney samples, we divided each area of kidney samples by the “correction coefficient” and calculated the “normalized area” (Supplementary Table S4).
To convert “normalized area” data into a molar ratio, we used “mol/area ratio” calculated by the molecular weight of each DNA adduct and the average area of the standard sample with a defined amount. The “mol/area ratios” of 5mdC, 5hmdC, 5fdC, dI, 8-OHdA, and 8-OHdG were 3.58 × 10− 19, 5.10 × 10− 19, 6.55 × 10− 19, 3.95 × 10− 19, 1.34 × 10− 18, and 8.58 × 10− 19, respectively. The molar quantity of DNA adducts in kidney samples was acquired as a product of the “normalized area” and “mol/area ratio” and was divided by the number of molecules of intact deoxyribonucleoside calculated from the genome DNA amount and composition with GC contents 41% to acquire the molar ratio (Supplementary Table S5). Data were shown with three-digit accuracy.
Collection of clinical data
Electronic medical records were reviewed to obtain variables including age, sex, weight, laboratory data, and pathological diagnosis. Information on smoking habits and a history of taking herbal medicine was also retrospectively collected. The estimated glomerular filtration rate (eGFR) was calculated using the Cockcroft-Gault equation: in males, eGFR = (140 - age) × weight (kg)/72 × plasma creatinine (mg/dL); in females, eGFR = 0.85 × (140 - age) × weight (kg)/72 × plasma creatinine (mg/dL).
The number of samples in which each DNA adduct was detected was compared by Fisher’s exact test between 3 tissue types: the renal cortex, renal medulla, and renal pelvis. The Mann–Whitney U test was performed to compare the molar ratio of DNA adducts in the following pairs: the cortex-RCC group and cortex-non-RCC group, the medulla-RCC group and medulla-non-RCC group, and the pelvis-UTUC group and pelvis-non-UTUC group. The relationships between each DNA adduct and clinical characteristics, including age, sex, smoking history, herb usage, eGFR, and the status of malignant tumors were assessed. Continuous variables, including age and eGFR, were analyzed by Spearman’s rank correlation test and other dichotomous variables were examined using the Mann–Whitney U test. Patients were classified as non-smokers [smoking history (−)] and ex-smokers and current smokers as smokers [smoking history (+)]. In pelvis-UTUC group, extended analyses were performed to compare the quantity of 5hmdC between solitary UTUC and multiple UTUC, and between UTUC with history of bladder cancer and UTUC without it using the Mann–Whitney U test. All calculations were performed using SPSS (IBM Corp. Released 2019. IBM SPSS Statistics for Windows, Version 26.0. Armonk, NY: IBM Corp). The significance of differences in the present study was set at p < 0.05 and all reported p values were two-sided.
Figure 2 shows representative chromatograms of 5hmdC in the “standard sample”, “blank sample”, and kidney samples of every group. After adopting the defined threshold, the prevalence of 5mdC, 5hmdC, 5fdC, dI, 8-OHdA, and 8-OHdG is shown in Table 2. 5mdC, 5hmdC, and dI were detected in all kidney samples. As shown in Supplementary Table S2, although C5-carboxy-2′-deoxycytidine (5cadC) was investigated, its levels were below the threshold. Fisher’s exact test revealed differences in the prevalence of the 6 molecules between tissue types, with 8-OHdG being more frequently detected in the renal pelvis than in the renal cortex and medulla (p = 0.048 and p = 0.038, respectively). The p values of the other pairs were > 0.05 (Supplementary Table S6). The prevalence of 5mdC, 5hmdC, 5fdC, dI, 8-OHdA, and 8-OHdG was also compared among the following cohort pairs: the cortex-RCC group and cortex-non-RCC group, medulla-RCC group and medulla-non-RCC group, and pelvis-UTUC group and pelvis-non-UTUC group. However, the p values of these pairs were > 0.05 (Supplementary Table S6).
Table 3 shows the p values of the Mann-Whitney U test on the relationships between clinical characteristics and DNA adducts in each tissue. Herb usage was associated with 5mdC and 5hmdC levels in the renal cortex. Although the number of herb users was very small in our cohort, as shown in the dot plots in Fig. 3, 5mdC and 5hmdC levels appeared to be reduced in the renal cortex of these patients. No significant differences were observed between smoking history and oxidative adducts, including 8-OHdA and 8-OHdG.
5hmdC levels differed between the pelvis-UTUC group and pelvis-non-UTUC group (p = 0.010). The median and range of the molar ratio of 5hmdC in the pelvis-UTUC group were 1.2 × 10− 3 and 6.5 × 10− 4 to 2.2 × 10− 3, respectively, while those in the pelvis-non-UTUC group were 2.0 × 10− 3 and 6.7 × 10− 4 to 2.6 × 10− 3, respectively. Furthermore, we compared the quantity of 5hmdC according to the two clinical features, tumor multicentricity of UTUC and history of bladder cancer in the patients with UTUC and the p values were 0.61 and 0.27, respectively (Fig. 4).
Figure 5 shows dot plots of the molar ratio of DNA adducts in the groups examined. The levels of the DNA adducts tested did not significantly differ between the cortex-RCC group and cortex-non-RCC group or between the medulla-RCC group and medulla-non-RCC group. RCC was not associated with the accumulation of specific DNA adducts in the renal cortex or medulla. Descriptive statistics, including the median, range, and standard deviation, of each DNA adduct according to the three tissue types are shown in Supplementary Table S7.
A relationship was observed between age and 5hmdC levels. Table 4 and Fig. 6 shows the relationships between DNA adducts, age, and eGFR. Specifically, age and eGFR correlated with 5hmdC levels in the renal cortex (r = − 0.46, p = 0.018 and r = 0.64, p < 0.001, respectively). Age also correlated with 5hmdC levels in the renal medulla (r = − 0.45, p = 0.042). However, neither 5mdC nor 5fdC correlated with age or eGFR in the renal cortex, renal medulla, or renal pelvis.
We also examined the relationships between all 6 DNA adducts (Supplementary Table S8) and found a negative correlation between 5mdC and dI in the renal cortex and medulla, but not in the renal pelvis (Fig. 7).
The present results revealed the presence of 5mdC, 5hmdC, 5fdC, dI, 8-OHdA, and 8-OHdG in the renal cortex, renal medulla, and urothelium of the renal pelvis. In our system, the adductomic approach theoretically facilitates the detection of multiple adducts other than those listed above. We previously identified several lipid peroxidation-induced DNA adducts, including butanone-etheno-2′-deoxyadenosine, in the gastric mucosa of gastric cancer patients . We examined the chromatograms of all urological tissues for the retention times and m/z of these lipid peroxidation-derived DNA adducts; however, they were not found. Therefore, DNA adduct profiles appear to markedly differ between the stomach and urological organs .
dI is produced by cadmium and triptolide, which cause renal injury [29, 30]. Regarding carcinogenesis, Wang et al. reported the potential of ultraperformance LC/MS to detect lipid and purine metabolites in urine as differential markers of urothelial cancer from RCC, and among the potential biomarkers examined, they identified dI in urine . However, in the present study, dI quantified in non-tumoral urinary tissue did not correlate with a malignant tumor status or eGFR.
8-OHdG and 8-OHdA are representative markers of DNA oxidative damage during oxidative stress . Among urological cancers, a high level of 8-OHdG was detected in the RCC tissue and urine of patients with schistosomiasis-associated squamous cell carcinoma in the bladder [14, 16]. Recent studies proposed the use of 8-OHdG in urine to diagnose gastric carcinoma and lung cancer [33, 34]; however, only a few studies examined 8-OHdG levels in tissue as a potential biomarker of urological cancer in clinical settings [14, 15]. In the present study, 8-OHdG levels did not significantly differ between UTUC and non-UTUC patients. Oxidative agents, such as tobacco, directly damage the lungs and gastrointestinal system and are excreted in urine; therefore, 8-OHdG levels in urine may be useful as a biomarker of these oxidative agents. In the present study, 8-OHdG levels did not correlate with the urological cancer status, which implies that the urinary system is only indirectly exposed to tobacco carcinogens via urine. We also investigated the relationship between oxidative DNA adducts in urinary tissues and smoking habits because tobacco smoke has been shown to induce oxidative DNA damage in the lungs ; however, the results obtained revealed no significant differences in 8-OHdG or 8-OHdA levels between smokers and non-smokers. Using leukocyte DNA, Asami et al. and Besaratinia et al. investigated the relationship between 8-OHdG and smoking habits, but reported inconsistent findings [36, 37]. This discrepancy may have been due to leukocytes being indirectly exposed to tobacco smoke, in contrast to direct exposure in the lungs. The relationship between urinary 8-OHdG levels and smoking habits assessed by LC/MS is also controversial: Andreoli et al. reported a positive relationship , whereas Harman et al. found a negative relationship . In the present study, oxidation-induced DNA adducts were not associated with smoking habits. The reasons for this result in our cohort may have been because we did not examine smoking habits in detail, such as quantity and duration, age is a confounding factor , and exposure to smoke in the urinary system is indirect. Further studies to compare DNA adducts in blood, urine, and urinary tissue and quantitatively evaluate smoking are needed in order to confirm the carcinogenicity of smoke-induced DNA adducts in the urinary system.
In the present study, 5hmdC levels in the non-tumoral renal pelvis were significantly lower in the UTUC group than in the non-UTUC group. However, no significant differences were observed in its precursor, 5mdC or its product, 5fdC, and another oxidation product, 5cadC was not detected. 5hmdC levels were previously reported to be reduced in many cancer tissues, including UC and premalignant lesions [25, 40,41,42,43]. 5mdC is oxidized and sequentially converted to 5hmdC, 5fdC, and 5cadC by ten-eleven translocation family proteins (TETs 1, 2, and 3). We previously reported that TET family protein levels were markedly lower in the tumoral portion than in the non-tumoral portion of the stomach , with TET1 playing an important role in the oxidation of 5mdC to 5hmdC. Foksinski et al. found a positive correlation between 5hmdC and TET1 transcript levels in malignant cell lines . Furthermore, Dziaman et al. showed that 5hmdC and TET1 expression levels were significantly lower in colon cancer and colon adenoma, a premalignant lesion, than in normal colon tissue . In addition, the silencing of TET1 reduced 5hmC levels and promoted cell proliferation by down-regulating the expression of AJAP1 in bladder cancer cells . Another reason for the reduced levels of 5hmdC in malignant tissue may be a high cell proliferation rate. Jin et al. showed that 5hmC levels were reduced in Ki-67-positive cells using an immunohistochemical analysis ; the loss of 5hmC in tumors was attributed to an elevated rate of cell proliferation, which appeared to result in passive decreases in 5hmC levels in human tumor tissue. Previous studies also reported a negative linear correlation between the cell proliferation rate and global level of 5hmC in several mammalian tissues [46, 47]. As described above, 5hmdC levels are reduced by a high cell proliferation rate and low TET1 expression. We hereafter discuss the effects of reductions in 5hmdC levels in the non-tumoral tissue of UTUC patients, not RCC patients. Reductions in 5hmdC levels in the non-tumoral portion of patients with malignant tumors have been reported in other diseases. We previously demonstrated that the global levels of 5hmdC were lower in the non-tumor gastric mucosa of gastric cancer patients than in control gastric mucosa . Other pathological changes in non-tumoral tissue have also been detected. Li et al. showed the clonal expansion of and somatic mutations in the normal urothelium of the non-tumoral tissue of UC patients . Therefore, reductions in 5hmdC levels in UTUC patients may be a factor contributing to differences in the characteristics of cancer, particularly local recurrence. Since local recurrence after kidney-sparing surgery has been reported in 37–52% of UTUC patients , major guidelines recommend nephroureterectomy in many cases . In contrast, local recurrence is rarer in patients with RCC, at a rate of 1–9% after partial nephrectomy and thermal ablation, while surgery robot-assisted partial nephrectomy recently achieved 0.5% [51, 52]. Therapeutic strategies for these two carcinomas markedly differ because of the risk of local recurrence. In the present study, we revealed reductions in 5hmdC levels in the non-tumoral tissue of UTUC patients, but not RCC patients. This result may reflect clinical differences in the two malignancies and, thus, reductions in the tissue levels of 5hmdC may be used for the risk stratification of the recurrence of UTUC. Additionally, we analyzed the relationships between the quantity of 5hmdC and the two clinical features, tumor multicentricity and the history of bladder cancer which are known as the prognostic factor of intravesical recurrence  in Fig. 4. However, the significant differences were not observed. A limitation of the present study is the small sample size. 5fdC, 8-OHdA, and 8-OHdG are rare DNA adducts in kidney tissue and the sample size were comparatively small in extended analysis about multicentricity and history of bladder cancer in patients with UTUC; therefore, investigations on a larger sample size are needed to reach concrete conclusions. Another limitation is that the control group included patients with other malignant tumors. For example, the cortex-non-RCC group included 8 UTUC patients, one ovarian cancer patient, and one sarcoma patient; therefore, the comparisons conducted in the present study were not necessarily between malignant tumor patients and non-tumor patients.
In the present study, herb usage appeared to contribute to reductions in 5hmdC levels in the renal cortex (Table 3). We did not extensively investigate all potentially detectable DNA adducts and their possible origins, particularly AA, in terms of urothelial carcinogenesis and kidney toxicity. The FDA issued a drug alert of the risk of AA in some of the herbs prescribed  and the herbs containing AA were not officially approved for prescription by Japanese government when outbreaks of AA related nephropathy occurred due to misidentification or private importation a decade ago in Japan . In the present cohort, a large proportion of patients used Chinese herbal medicine; however, to the best of our knowledge and based on a detailed examination of clinical records, herbs such as kakkonto, daikenchuto, rikkunshito, daiokanzoto, and mashiningan, which are commonly prescribed in Japanese clinical practice, do not contain AA. We cannot conclude whether AA contributes to urothelial diseases in Japanese patients due to small number of patients investigated herein. Therefore, further studies on a larger number of patients are needed to clarify the effects of herb use on DNA methylation and hydroxymethylation.
In the present study, the status of RCC did not correlate with 5hmdC levels in the renal parenchyma. However, negative correlations were observed between age and 5hmdC levels in the renal cortex and medulla, as shown in Table 4 and Fig. 6. Therefore, the relationship between 5hmdC levels and aging appears to differ among organs. Previous studies reported age-dependent decreases in the global levels of 5hmdC in human T cells and the cerebral cortex [56, 57], whereas increases were detected in the cerebellum and sperm [58, 59]. Based on these findings, the mechanisms underlying 5hmdC reductions with aging warrant further study.
As shown in Fig. 6, eGFR positively correlated with 5hmdC levels, namely, 5hmdC levels decreased with the deterioration of renal function. Using a mouse system, Huang et al. reported down-regulation of the mRNA expression of TET1 and TET2 after ischemia-reperfusion injury, and subsequent reduction of global levels of 5hmdC in the kidney . Furthermore, chronic hypoxia was shown to impair the activity of TET and alter the demethylation status in the mouse kidney . Therefore, the low expression of TET may also lead to 5hmdC reductions with decreases in eGFR in humans.
Investigations to obtain a more comprehensive understanding of DNA adducts and modifications in urological systems are still in the burgeoning stage. The more extensive identification of other DNA adducts in more cases is needed. And the present study used LC/MS to detect the DNA adducts which is not widely available in a clinical practice. To apply the findings of DNA adducts acquired by LC/MS to a clinical setting, the investigation which evaluates the DNA adducts by immunohistochemistry staining is required.
5hmdC levels in the normal renal pelvis were significantly lower in UTUC patients than in non-UTUC patients. Pathological changes also occurred in the normal urothelial mucosa, which may explain the higher risk of local recurrence in UTUC patients than in RCC patients.
Availability of data and materials
All processed files are presented in additional files. Raw data files are available upon request.
Renal cell carcinoma
Upper urinary tract urothelial carcinoma
estimated glomerular filtration rate
2′-deoxycytidine containing two carbon-13 atoms
2′-deoxythymidine containing two carbon-13 atoms
2′-deoxyadenosine containing two carbon-13 atoms
2′-deoxyguanosine containing two carbon-13 atoms
N7-pyridyloxobutyl-2′-deoxyguanosine containing four deuterium atoms
O2-pyridyloxobutyl-2′-deoxythymidine containing four deuterium atoms
Weinberg RA. The biology of cancer. 2nd ed: Garland science; 2013.
Yu J, Ni M, Xu J, Zhang H, Gao B, Gu J, et al. Methylation profiling of twenty promoter-CpG islands of genes which may contribute to hepatocellular carcinogenesis. BMC Cancer. 2002;2:29.
Bearzatto A, Conte D, Frattini M, Zaffaroni N, Andriani F, Balestra D, et al. p16(INK4A) Hypermethylation detected by fluorescent methylation-specific PCR in plasmas from non-small cell lung cancer. Clin Cancer Res. 2002;8(12):3782–7.
Lee MG, Huh JS, Chung SK, Lee JH, Byun DS, Ryu BK, et al. Promoter CpG hypermethylation and downregulation of XAF1 expression in human urogenital malignancies: implication for attenuated p53 response to apoptotic stresses. Oncogene. 2006;25(42):5807–22.
Martinez R, Setien F, Voelter C, Casado S, Quesada MP, Schackert G, et al. CpG island promoter hypermethylation of the pro-apoptotic gene caspase-8 is a common hallmark of relapsed glioblastoma multiforme. Carcinogenesis. 2007;28(6):1264–8.
Totsuka Y, Watanabe M, Lin Y. New horizons of DNA adductome for exploring environmental causes of cancer. Cancer Sci. 2021;112(1):7–15.
Kanaly RA, Hanaoka T, Sugimura H, Toda H, Matsui S, Matsuda T. Development of the adductome approach to detect DNA damage in humans. Antioxid Redox Signal. 2006;8(5–6):993–1001.
Matsuda T, Tao H, Goto M, Yamada H, Suzuki M, Wu Y, et al. Lipid peroxidation-induced DNA adducts in human gastric mucosa. Carcinogenesis. 2013;34(1):121–7.
Chou PH, Kageyama S, Matsuda S, Kanemoto K, Sasada Y, Oka M, et al. Detection of lipid peroxidation-induced DNA adducts caused by 4-oxo-2(E)-nonenal and 4-oxo-2(E)-hexenal in human autopsy tissues. Chem Res Toxicol. 2010;23(9):1442–8.
Ohnishi I, Iwashita Y, Matsushita Y, Ohtsuka S, Yamashita T, Inaba K, et al. Mass spectrometric profiling of DNA adducts in the human stomach associated with damage from environmental factors. Genes Environ. 2021;43(1):12.
Kasai H. What causes human cancer? Approaches from the chemistry of DNA damage. Genes Environ. 2016;38:19.
Kumondai M, Hosono H, Orikasa K, Arai Y, Arai T, Sugimura H, et al. CYP2A13 genetic polymorphisms in relation to the risk of bladder Cancer in Japanese smokers. Biol Pharm Bull. 2016;39(10):1683–6.
Kumondai M, Hosono H, Orikasa K, Arai Y, Arai T, Sugimura H, et al. Genetic polymorphisms of CYP2A6 in a case-control study on bladder Cancer in Japanese smokers. Biol Pharm Bull. 2016;39(1):84–9.
Miyake H, Hara I, Kamidono S, Eto H. Prognostic significance of oxidative DNA damage evaluated by 8-hydroxy-2′-deoxyguanosine in patients undergoing radical nephrectomy for renal cell carcinoma. Urology. 2004;64(5):1057–61.
Miyake H, Hara I, Kamidono S, Eto H. Oxidative DNA damage in patients with prostate cancer and its response to treatment. J Urol. 2004;171(4):1533–6.
Gouveia MJ, Santos J, Brindley PJ, Rinaldi G, Lopes C, Santos LL, et al. Estrogen-like metabolites and DNA-adducts in urogenital schistosomiasis-associated bladder cancer. Cancer Lett. 2015;359(2):226–32.
Arlt VM, Pfohl-Leszkowicz A, Cosyns J-P, Schmeiser HH. Analyses of DNA adducts formed by ochratoxin a and aristolochic acid in patients with Chinese herbs nephropathy. Mutat Res. 2001;494(1–2):143–50.
Hoang ML, Chen CH, Chen PC, Roberts NJ, Dickman KG, Yun BH, et al. Aristolochic Acid in the etiology of renal cell carcinoma. Cancer Epidemiol Biomark Prev. 2016;25(12):1600–8.
Nortier JL, Martinez MC, Schmeiser HH, Arlt VM, Bieler CA, Petein M, et al. Urothelial carcinoma associated with the use of a Chinese herb (Aristolochia fangchi). N Engl J Med. 2000;342(23):1686–92.
Schmeiser HH, Bieler CA, Wiessler M, Van Ypersele De Strihou C, Cosyns JP. Detection of DNA adducts formed by aristolochic acid in renal tissue from patients with Chinese herbs nephropathy. Cancer Res. 1996;56(9):2025–8.
Schmeiser HH, Kucab JE, Arlt VM, Phillips DH, Hollstein M, Gluhovschi G, et al. Evidence of exposure to aristolochic acid in patients with urothelial cancer from a Balkan endemic nephropathy region of Romania. Environ Mol Mutagen. 2012;53(8):636–41.
Jelakovic B, Karanovic S, Vukovic-Lela I, Miller F, Edwards KL, Nikolic J, et al. Aristolactam-DNA adducts are a biomarker of environmental exposure to aristolochic acid. Kidney Int. 2012;81(6):559–67.
Luan Y, Xing G, Ren J, Gu J. Role of hepatic cytochrome P450 enzymes in the detoxication of aristolochic acid I; effects on DNA adduct, mutation, and tumor formation. Genes Environ. 2015;37:11.
Burger M, Catto JW, Dalbagni G, Grossman HB, Herr H, Karakiewicz P, et al. Epidemiology and risk factors of urothelial bladder cancer. Eur Urol. 2013;63(2):234–41.
Munari E, Chaux A, Vaghasia AM, Taheri D, Karram S, Bezerra SM, et al. Global 5-Hydroxymethylcytosine levels are profoundly reduced in multiple genitourinary malignancies. PLoS One. 2016;11(1):e0146302.
Shi DQ, Ali I, Tang J, Yang WC. New insights into 5hmC DNA modification: generation, Distribution and Function. Front Genet. 2017;8:100.
Kharat SS, Ding X, Swaminathan D, Suresh A, Singh M, Sengodan SK, et al. Degradation of 5hmC-marked stalled replication forks by APE1 causes genomic instability. Sci Signal. 2020;13(645):eaba8091.
Kanai Y, Nishihara H, Miyagi Y, Tsuruyama T, Taguchi K, Katoh H, et al. The Japanese Society of Pathology Guidelines on the handling of pathological tissue samples for genomic research: standard operating procedures based on empirical analyses. Pathol Int. 2018;68(2):63–90.
Huang W, Liu C, Xie L, Wang Y, Xu Y, Li Y. Integrated network pharmacology and targeted metabolomics to reveal the mechanism of nephrotoxicity of triptolide. Toxicol Res (Camb). 2019;8(6):850–61.
Guan T, Xin Y, Zheng K, Wang R, Zhang X, Jia S, et al. Metabolomics analysis of the effects of quercetin on renal toxicity induced by cadmium exposure in rats. Biometals. 2021;34(1):33–48.
Wang Z, Liu X, Liu X, Sun H, Guo Z, Zheng G, et al. UPLC-MS based urine untargeted metabolomic analyses to differentiate bladder cancer from renal cell carcinoma. BMC Cancer. 2019;19(1):1195.
Kalam MA, Haraguchi K, Chandani S, Loechler EL, Moriya M, Greenberg MM, et al. Genetic effects of oxidative DNA damages: comparative mutagenesis of the imidazole ring-opened formamidopyrimidines (Fapy lesions) and 8-oxo-purines in simian kidney cells. Nucleic Acids Res. 2006;34(8):2305–15.
Borrego S, Vazquez A, Dasi F, Cerda C, Iradi A, Tormos C, et al. Oxidative stress and DNA damage in human gastric carcinoma: 8-Oxo-7′8-dihydro-2′-deoxyguanosine (8-oxo-dG) as a possible tumor marker. Int J Mol Sci. 2013;14(2):3467–86.
Loft S, Svoboda P, Kawai K, Kasai H, Sorensen M, Tjonneland A, et al. Association between 8-oxo-7,8-dihydroguanine excretion and risk of lung cancer in a prospective study. Free Radic Biol Med. 2012;52(1):167–72.
Asami S, Manabe H, Miyake J, Tsurudome Y, Hirano T, Yamaguchi R, et al. Cigarette smoking induces an increase in oxidative DNA damage, 8-hydroxydeoxyguanosine, in a central site of the human lung. Carcinogenesis. 1997;18(9):1763–6.
Asami S, Hirano T, Yamaguchi R, Tomioka Y, Itoh H, Kasai H. Increase of a type of oxidative DNA damage, 8-hydroxyguanine, and its repair activity in human leukocytes by cigarette smoking. Cancer Res. 1996;56(11):2546–9.
Besaratinia A, Van Schooten FJ, Schilderman PA, De Kok TM, Haenen GR, Van Herwijnen MH, et al. A multi-biomarker approach to study the effects of smoking on oxidative DNA damage and repair and antioxidative defense mechanisms. Carcinogenesis. 2001;22(3):395–401.
Andreoli R, Mutti A, Goldoni M, Manini P, Apostoli P, De Palma G. Reference ranges of urinary biomarkers of oxidized guanine in (2′-deoxy)ribonucleotides and nucleic acids. Free Radic Biol Med. 2011;50(2):254–61.
Harman SM, Liang L, Tsitouras PD, Gucciardo F, Heward CB, Reaven PD, et al. Urinary excretion of three nucleic acid oxidation adducts and isoprostane F (2) alpha measured by liquid chromatography-mass spectrometry in smokers, ex-smokers, and nonsmokers. Free Radic Biol Med. 2003;35(10):1301–9.
Pfeifer GP, Kadam S, Jin SG. 5-hydroxymethylcytosine and its potential roles in development and cancer. Epigenetics Chromatin. 2013;6(1):10.
Dziaman T, Gackowski D, Guz J, Linowiecka K, Bodnar M, Starczak M, et al. Characteristic profiles of DNA epigenetic modifications in colon cancer and its predisposing conditions-benign adenomas and inflammatory bowel disease. Clin Epigenetics. 2018;10:72.
Jin SG, Jiang Y, Qiu R, Rauch TA, Wang Y, Schackert G, et al. 5-Hydroxymethylcytosine is strongly depleted in human cancers but its levels do not correlate with IDH1 mutations. Cancer Res. 2011;71(24):7360–5.
Du C, Kurabe N, Matsushima Y, Suzuki M, Kahyo T, Ohnishi I, et al. Robust quantitative assessments of cytosine modifications and changes in the expressions of related enzymes in gastric cancer. Gastric Cancer. 2015;18(3):516–25.
Foksinski M, Zarakowska E, Gackowski D, Skonieczna M, Gajda K, Hudy D, et al. Profiles of a broad spectrum of epigenetic DNA modifications in normal and malignant human cell lines: proliferation rate is not the major factor responsible for the 5-hydroxymethyl-2′-deoxycytidine level in cultured cancerous cell lines. PLoS One. 2017;12(11):e0188856.
Yan YL, Huang ZN, Zhu Z, Cui YY, Li MQ, Huang RM, et al. Downregulation of TET1 promotes bladder Cancer cell proliferation and invasion by reducing DNA Hydroxymethylation of AJAP1. Front Oncol. 2020;10:667.
Gackowski D, Zarakowska E, Starczak M, Modrzejewska M, Olinski R. Tissue-specific differences in DNA modifications (5-Hydroxymethylcytosine, 5-Formylcytosine, 5-Carboxylcytosine and 5-Hydroxymethyluracil) and their interrelationships. PLoS One. 2015;10(12):e0144859.
Bachman M, Uribe-Lewis S, Yang X, Williams M, Murrell A, Balasubramanian S. 5-Hydroxymethylcytosine is a predominantly stable DNA modification. Nat Chem. 2014;6(12):1049–55.
Li R, Du Y, Chen Z, Xu D, Lin T, Jin S, et al. Macroscopic somatic clonal expansion in morphologically normal human urothelium. Science. 2020;370(6512):82–9.
Cutress ML, Stewart GD, Zakikhani P, Phipps S, Thomas BG, Tolley DA. Ureteroscopic and percutaneous management of upper tract urothelial carcinoma (UTUC): systematic review. BJU Int. 2012;110(5):614–28.
Roupret M, Babjuk M, Burger M, Capoun O, Cohen D, Comperat EM, et al. European Association of Urology guidelines on upper urinary tract urothelial carcinoma: 2020 update. Eur Urol. 2021;79(1):62–79.
Kriegmair MC, Bertolo R, Karakiewicz PI, Leibovich BC, Ljungberg B, Mir MC, et al. Systematic review of the Management of Local Kidney Cancer Relapse. Eur Urol Oncol. 2018;1(6):512–23.
Yoon YE, Lee HH, Kim KH, Park SY, Moon HS, Lee SR, et al. Focal therapy versus robot-assisted partial nephrectomy in the management of clinical T1 renal masses: a systematic review and meta-analysis. Medicine (Baltimore). 2018;97(45):e13102.
Azemar MD, Comperat E, Richard F, Cussenot O, Roupret M. Bladder recurrence after surgery for upper urinary tract urothelial cell carcinoma: frequency, risk factors, and surveillance. Urol Oncol. 2011;29(2):130–6.
Aristolochic Acid: Safety alert. In: Safety Alerts & Advisories. U.S. Food and Drug Administration. 2015. http://wayback.archive-it.org/7993/20171114232638/https://www.fda.gov/Food/RecallsOutbreaksEmergencies/SafetyAlertsAdvisories/ucm095272.htm. Accessed 11 May 2021.
Tanaka A, Nishida R, Yoshida T, Koshikawa M, Goto M, Kuwahara T. Outbreak of Chinese herb nephropathy in Japan: are there any differences from Belgium? Intern Med. 2001;40(4):296–300.
Truong TP, Sakata-Yanagimoto M, Yamada M, Nagae G, Enami T, Nakamoto-Matsubara R, et al. Age-dependent decrease of DNA Hydroxymethylation in human T cells. J Clin Exp Hematop. 2015;55(1):1–6.
Dzitoyeva S, Chen H, Manev H. Effect of aging on 5-hydroxymethylcytosine in brain mitochondria. Neurobiol Aging. 2012;33(12):2881–91.
Szulwach KE, Li X, Li Y, Song CX, Wu H, Dai Q, et al. 5-hmC-mediated epigenetic dynamics during postnatal neurodevelopment and aging. Nat Neurosci. 2011;14(12):1607–16.
Jenkins TG, Aston KI, Cairns BR, Carrell DT. Paternal aging and associated intraindividual alterations of global sperm 5-methylcytosine and 5-hydroxymethylcytosine levels. Fertil Steril. 2013;100(4):945–51.
Huang N, Tan L, Xue Z, Cang J, Wang H. Reduction of DNA hydroxymethylation in the mouse kidney insulted by ischemia reperfusion. Biochem Biophys Res Commun. 2012;422(4):697–702.
Gu Y, Chen J, Zhang H, Shen Z, Liu H, Lv S, et al. Hydrogen sulfide attenuates renal fibrosis by inducing TET-dependent DNA demethylation on klotho promoter. FASEB J. 2020;34(9):11474–87.
We acknowledge Drs. Masako Suzuki, Takuya Kitamoto, Minako Kondo, and Chiharu Uchida in the Advanced Research Facilities and Services, the Preeminent Medical Photonics Education & Research Center, Hamamatsu University School of Medicine for their excellent support on the maintenance and maneuvering of instruments. We also acknowledge Dr. Satoshi Baba and Ms. Juri Takahashi as well as other members of the Department of Diagnostic Pathology, Hamamatsu University Hospital.
This work was supported by JSPS KAKENHI Grant Numbers JP22659072, JP24659161, JP26670187, and JP16K15256, AMED Grant Numbers JP19ck0106264 and JP20ck0106545, the Smoking Research Foundation, and a HUSM Grant-in-Aid by Hamamatsu University School of Medicine.
Ethics approval and consent to participate
All procedures were performed in accordance with the Declaration of Helsinki and the use of residual pathological tissue in this research was approved by the Institutional Review Board of Hamamatsu University School of Medicine (No. 20–011).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Clinicopathological characteristics of cases. Supplementary Table S2: Standard compounds of DNA adducts and mass transition. Supplementary Table S3: Raw area data on DNA adducts in kidney samples. Supplementary Table S4: “Correction coefficient” and “normalized area” data on DNA adducts. Supplementary Table S5: The molar ratio of DNA adducts. Supplementary Table S6: Differences in the prevalence of DNA adducts by tissue types and the tumor status. Supplementary Table S7: Descriptive statistics of DNA adducts according to tissue types. Supplementary Table S8: A correlation analysis of DNA adducts.
About this article
Cite this article
Matsushita, Y., Iwashita, Y., Ohtsuka, S. et al. A DNA adductome analysis revealed a reduction in the global level of C5-hydroxymethyl-2′-deoxycytidine in the non-tumoral upper urinary tract mucosa of urothelial carcinoma patients. Genes and Environ 43, 52 (2021). https://doi.org/10.1186/s41021-021-00228-9