Skip to main content

Association between genome-wide copy number variation and arsenic-induced skin lesions: a prospective study

Abstract

Background

Exposure to arsenic in drinking water is a global health problem and arsenic-induced skin lesions are hallmark of chronic arsenic toxicity. We and others have reported germline genetic variations as risk factors for such skin lesions. The role of copy number variation (CNV) in the germline DNA in this regard is unknown.

Methods

From a large prospectively followed-up cohort, exposed to arsenic, we randomly selected 2171 subjects without arsenic-induced skin lesions at enrollment and genotyped their whole blood DNA samples on Illumina Cyto12v2.1 SNP chips to generate DNA copy number. Participants were followed up every 2 years for a total of 8 years, especially for the development of skin lesions. In Cox regression models, each CNV segment was used as a predictor, accounting for other potential covariates, for incidence of skin lesions.

Result

The presence of genomic deletion(s) in a number of genes (OR5J2, GOLGA6L7P, APBA2, GALNTL5, VN1R31P, PHKG1P2, SGCZ, ZNF658) and lincRNA genes (RP11-76I14.1, CTC-535 M15.2, RP11-73B2.2) were associated with higher risk [HR between 1.67 (CI 1.3-2.1) and 2.15 (CI 1.5-2.9) for different CNVs] for development of skin lesions independent of gender, age, and arsenic exposure. Some deletions had stronger effect in a specific gender (ZNF658 in males, SGCZ in females) and some had stronger effect in higher arsenic exposure (lincRNA CTD-3179P9.1) suggesting a possible gene-environment interaction.

Conclusion

This first genome-wide CNV study in a prospectively followed-up large cohort, exposed to arsenic, suggests that DNA deletion in several genes and lincRNA genes may predispose an individual to a higher risk of development of arsenic-induced skin lesions.

Peer Review reports

Background

Over 100 million individuals worldwide are exposed to arsenic through drinking water, including 28–57 million in Bangladesh [1] and 13 million in the United States [2]. Arsenic is a class-I human carcinogen, and chronic exposure to high levels (>300 μg/L) is associated with increased risk for a wide array of diseases including cancers of the lung [3], liver [4], bladder [5, 6], kidney [7, 8], as well as neurological [9], metabolic [10] and cardiovascular [11,12,13,14,15] diseases, skin lesions [16,17,18,19] and maternal health [20]. Chronic arsenic exposure through drinking water is associated with an increase in mortality [21]. Most arsenic-related cancers have a long latency period, but arsenic-induced skin lesions appear relatively early [6, 22]. Moreover, hyperkeratosis may be considered as precursors to arsenic-induced basal and squamous cell carcinoma of skin [23]. Smith and Steinmaus [24] have reviewed a large number of epidemiologic studies, mostly case–control, showing the association between arsenic exposure and skin lesions.

Previously, using case–control design in a Bangladeshi population, in the first genome wide association study (GWAS) in arsenic, our group found some single nucleotide polymorphisms (SNP) to be associated with arsenic metabolism [25]. Using clinical follow-up data of one of the largest cohorts exposed to arsenic through drinking water, our group also presented evidence that a higher risk of arsenic-induced skin lesions was found in the male gender, higher age and higher arsenic exposure [16]. We have clinically followed-up a large cohort in Bangladesh exposed to different levels of arsenic through drinking water [26, 27]. A large number of subjects from this cohort were randomly selected for GWAS and we have done SNP genotyping using oligonucleotide based arrays from whole blood DNA collected at baseline. We have previously demonstrated the utility of these oligonucleotide based arrays to detect and interpret copy number (CN) changes in clinical samples [28]. A copy number variant (CNV) is a term collectively used to describe gains or losses of DNA sequence >1 kb in length. These may have a direct effect on transcription and transcriptional regulation, which in turn may be a cause for disease susceptibility and phenotypic variation [29].

It may be noted that CNVs represent a large class of genomic variation that was not well studied in the past, but is now gaining the attention of many investigators [30,31,32]. CNVs have already been reported to be associated with autism [33], schizophrenia [34, 35] and Crohn’s disease [36, 37]. In a tumor tissue based study, DNA losses at chromosomes 1q21.1, 7p22.3, 9q12 and 19q13.31 have been reported in arsenic-related lung squamous cell carcinoma [38]. In this study, we evaluate whether there is any association of CNVs in germ line DNA in the development of arsenic-induced skin lesion. To our knowledge, we present the first paper addressing the role of germ line CNVs in the development of arsenic-induced skin lesions in a population exposed to arsenic through drinking water.

Methods

Study subjects

The Health Effects of Arsenic Longitudinal Study (HEALS) was designed to investigate the health effects of arsenic exposure through drinking water in a population-based sample of adults in Araihazar, Bangladesh [26]. The study methods have been described previously [26]. The study protocol was approved by the Institutional Review board of The University of Chicago, Columbia University, and the Bangladesh Medical Research Council. Informed consent was obtained from all participants. At the start of the study, we identified 12,050 eligible individuals for recruitment from the enumerated total of approximately 65,000 residents in the study area. Between October 2000 and May 2002, we recruited married individuals aged 18–75 years who had been residing in the study area for at least 5 years. A total of 11,746 men and women enrolled into the HEALS cohort. At the baseline interview, trained study physicians blinded to the arsenic concentrations in participant’s drinking water conducted in-person interviews and clinical evaluations including skin examination. They collected spot urine and blood samples from the participants according to a structured protocol. Participants were contacted for a follow-up examination at an interval of 2 years, which followed the same baseline protocol. For this study, we have utilized the biological samples collected at baseline, the clinical skin evaluation data at baseline and subsequent four biennial follow-ups (96 months). We randomly selected 2332 HEALS participants for genome-wide SNP genotyping. Among those 2332 HEALS participants, 2270 did not have any arsenic-induced skins lesion at the time of enrollment. In this study we only considered these 2270 HEALS participants with no prevalent skin lesions at baseline.

Arsenic-induced skin lesions

A structured protocol was used to ascertain arsenic-induced skin lesions by trained study physicians. The study physicians recorded the presence or absence of melanosis (hyperpigmentation), leukomelanosis (hypopigmentation) or keratosis (thickening of skin typically on the palms and soles) [16]. All the study physicians were specially trained to diagnose arsenic-induced skin lesions. We ascertained incident skin lesion cases in a prospective fashion using a structured protocol [16]. For the present study, “skin lesion” was classified as presence of any of these three or a combination of them.

Arsenic exposure

Well water arsenic (WAs) concentrations of all 5966 wells in the study area were measured by graphite furnace atomic absorption spectrometry, with a detection limit of 5 μg/L. Samples below the limit of detection were subsequently reanalyzed by inductively coupled plasma mass spectrometry, with a detection limit of 0.1 μg/L [39]. In our present study, the 25th percentile, 50th percentile and 75th percentile of well water arsenic were 12 μg/L, 56 μg/L and 142 μg/L respectively (see Additional file 1: Table S1). It may be noted that the 25th percentile was close to the WHO guideline for arsenic in drinking water (10 μg/L) and the 50th percentile was close to the Bangladesh national standard for arsenic in drinking water (50 μg/L). The urinary total arsenic concentration was measured by graphite furnace atomic absorption spectrometry [40]. Urinary creatinine was measured by a colorimetric method based on the Jaffe reaction described by Heinegard and Tiderstrom [41]. The urinary arsenic was measured from a spot urine collection. To take into account the hydration status, we used the urinary arsenic creatinine ratio (UACR) as measure of arsenic exposure. The log2-transformed UACR showed strong correlation to the log2-transformed well water arsenic concentration (r = 0.66, see Additional file 2: Figure S1).

Illumina SNP array

DNA was extracted from whole blood using the Flexigene kit (Qiagen, USA). Quantification was done using a NanoDrop 1000. According to Illumina protocol 250 ng of DNA was genotyped on Cyto12 v2.1 chips with 294,602 markers (289,773 SNP markers and 4829 copy number markers) and read on the BeadArray Reader. Image data was processed in GenomeStudio software V2010.3. After cluster generation, the genotype calls, B allele frequency and log2R ratio (LRR) were calculated. In GenomeStudio, the copy number (CN) is expressed as log2R ratio (LRR). For a particular locus, if a DNA sample has 2 copies (CN = 2), the ratio of signal intensity in a test sample to reference (which also should have CN = 2) would be 1 and thus log2 of the ratio (LRR) would be log2 1 = 0. In the same way, a sample with CN = 1 (intensity would be half compared to the reference) would have LRR = log2 0.5 = −1, whereas a sample with CN = 4 (expected intensity would be double the reference) would have LRR = log2 2 = 1.

Quality control (QC) and filtering

We excluded the markers in sex chromosomes (n = 17,442). From the remaining 277,160 autosomal markers (272,663 SNP and 4497 copy number markers), a total of 3064 SNP markers (only 1.12%) were excluded due to poor performance. The remaining 274,096 autosomal markers were considered for further analysis. For each sample, standard deviation (SD) of LRR of these 274,096 autosomal markers was calculated. Another 70 samples with SD of LRR >0.28 and 29 samples with call rate < 99.0% were excluded (see Additional file 3: Figure S2). Thus, we finally used high quality genomic data from a total of 2171 HEALS participants who did not have any arsenic-induced skin lesion at baseline and were prospectively followed-up for development of any arsenic-induced skin lesions. Characteristics of the study subjects are shown in Additional file 1: Table S1.

Genome-wide CN analysis

GenomeStudio generated LRR data was imported into Partek genomic suite and transformed to CN data in linear scale [42]. Standard Principal Component Analysis (PCA) and a sample histogram were generated as part of QC. After obtaining the CN value for each locus, to identify the genomic regions with amplification, normal CN or deletion, we used a genomic segmentation algorithm [42]. By “genomic region” in a particular sample we mean a stretch of DNA showing amplification or deletion. A genomic region with CN variation in one sample may or may not fully overlap with a genomic region in another sample. By amplification or deletion “segment”, we mean the stretch of amplified or deleted segment that is common in at least 5% of the samples in this paper.

For the data from the Cyto12 v2.1 chips, the genomic segmentation was done with a setting of a minimum of 6 markers, signal to noise 0.3, and p-value threshold of 0.001 for two neighboring regions having significantly differing means. A genomic region was considered as amplified if the geometric mean CN was >2.3 and a deletion if the mean was <1.7. We restricted the analysis for the autosomes only (2.8% of the segments had CN > 2.3, 22.12% had CN <1.7 and the 75.0% had a copy number between 1.7-2.3). The length of a genomic segment was calculated from the genomic location of the start and end SNP for that genomic segment. In this paper, we reported a genomic segment in a sample to have amplification (0:no amplification, 1:amplification) or deletion (0:no deletion, 1:deletion) only if it was at least 5 kb in size and the geometric mean of the CN within the genomic boundary of the segment for that particular sample was >2.3 or <1.7 respectively. We used the CN status for each segment as a binary predictor for development of arsenic-induced skin lesions (0: no skin lesion, 1: skin lesion) in survival analysis.

Statistical analysis

To compare the continuous variables (e.g. age, UACR, WAs, BMI,), we used one-way analysis of variance (ANOVA). For the categorical variables we used chi-square tests. We used both Kaplan-Meier curves and Cox Regression analysis. In survival analysis, the event was defined as any skin lesion detected during the follow-up visits. The time (months) was calculated from the enrollment to the first detection of any skin lesion (for “event”) and from the enrollment to the last follow-up (for the “censored”). We dichotomized the continuous variables by median value – age (0: age = <median 38, vs. 1: age > median), UACR (0: urinary arsenic creatinine ratio = <median 192 μg/g of creatinine vs. 1: >median). In Cox regression, we used the following model:

$$ \ln \left(\frac{H(t)}{H_0(t)}\right)={b}_1 Gene+{b}_2 Gender+{b}_3 Age+{b}_4 UACR $$

For the interaction models, we used the following:

$$ \ln \left(\frac{H(t)}{H_0(t)}\right)={b}_1 Gene+{b}_2 Gender+{b}_3 Age+{b}_4 UACR+{b}_5\times Gene\times Gender $$
$$ \ln \left(\frac{H(t)}{H_0(t)}\right)={b}_1 Gene+{b}_2 Gender+{b}_3 Age+{b}_4 UACR+{b}_5\times Gene\times UACR $$

Here, H(t)/H0(t) is the Hazard Ratio (HR). So in the model, if we use genomic deletion as a predictor (Gene), then the quantity exp.(b1) can be interpreted as the instantaneous relative risk of an event, at any time, for an individual with the genomic deletion present compared to an individual without the genomic deletion, given both the individuals are the same on all other covariates. For multiple testing, we used Bonferroni correction. The significance threshold was set at 4.4 E-5, which is 0.05/1135, the number of deletions identified and tested in this study.

Result

Among the 2171 HEALS subjects (m = 1032, f = 1139) without arsenical skin lesions at baseline, a total of 301 male (29%) and 115 female (10%) subjects developed skin lesions during the 8 years of follow-up. Kaplan-Meir plots by gender and arsenic exposure are shown in Additional file 4: Figure S3. The higher incidence of skin lesions among male subjects [HR 2.76 (CI 2.2 – 3.4)] compared to females of similar age and arsenic exposure, is consistent with our previous report on shorter follow-up on a larger number of subjects from the same cohort [27]. Higher age of the individual (reflecting duration of arsenic exposure) and higher Urinary Arsenic Creatinine Ratio (UACR) (reflecting the level of arsenic exposure) were also associated with higher risk [HR 2.97 (CI 2.3-3.7) and HR1.6 (CI 1.3-1.9) respectively] for the development of arsenic-induced skin lesions (see Additional file 4: Figure S3 and Additional file 5: Figure S4).

Structural variations were detected using Illumina Human Cyto12 v2.1 SNP chips. We identified a total of 1135 segments (at least 5 kb or longer in length), which showed CN loss <1.7 in at least 5% of the samples. We also found a total of 126 segments (at least 5 kb or longer in length), which showed CN gain >2.3 in at least 5% of the samples. Figure 1 shows an example of CNVs of a given region. The bottom-most panel (panel C) shows a deleted region in a number of samples. The genomic coordinates shown in x-axis correspond to the GRCh37/hg19 assembly. The length of deletion varies from sample to sample, however there are overlaps. In the middle panel (panel B), the y-axis shows the number of samples having deletion in that genomic region. So within that genomic region, there were five consecutive segments showing deletion. The top panel (panel A) shows annotation of the region as found in Refseq and Ensembl database.

Fig. 1
figure 1

Copy number variants detected in 8p22 region. Top panel (panel a) shows mapping of the region to Refseq and Ensembl database; middle panel (panel b): y-axis shows the number of samples showing deletion; lower panel (panel c) shows different samples having deletion of different length (shown on x-axis), but the overlap is clearly seen

To test the association between copy number loss/deletion and development of arsenic-induced skin lesions, we dichotomized the CN status of a segment as 0: no deletion and 1: deletion. Then, using Cox regression, we tested each of the 1135 segments for its potential role in the development of arsenic-induced skin lesions. We mentioned above that gender, age and UACR affect the risk of skin lesion. Therefore, along with the segments as a predictor, we also entered the gender (male vs. female), age (= < median 38 yrs. vs. >median) and UACR (= < median 192 μg/g of creatinine vs. >median) as covariates to see if the association(s) between segmental deletion and risk of skin lesion was independent of the covariates. Complete results from all these 1135 Cox regression models with HR and 95%CI of all the variables (genomic segment, gender, age and UACR) entered into the models are presented in Additional file 6: Table S2. We found a total of 24 segments covering 10 cytoband regions, deletion of which were significantly (Bonferroni p = <0.05) associated with higher risk of skin lesion development (HR ranging between 1.67 and 2.15 for different segments, see Table 1). The fact that multiple segments (mostly consecutive regions) within the same cytoband region (e.g. four successive segments in 5q34 cytoband region covering long intervening noncoding RNA (lincRNA) gene CTC-535 M15.2) were statistically significant for skin lesion risk, further strengthens our findings. Interestingly, out of these 10 cytoband regions, 3 of them (2q12.1, 5q34 and 7q11.21 shown in bold font in the Table 1) cover known lincRNA suggesting the possible significance of deletion of lincRNA as a risk factor for arsenic-induced skin lesions. For each CNV, the nearest gene is shown in Table 1. Our previous GWAS showed that two SNPs, rs9527 and rs11191659, were associated with arsenic metabolism [25]. However, conditioning on those two SNPs or smoking habit and chewing betel leaf did not change the effects of any of these genomic deletions on the development of skin lesions (see Additional file 7: Figure S5 and Additional file 8: Figure S6 respectively). Out of the 24 significant deletion segments shown in Table 1, 13 were reported in the database of genomic variants (DGV) (Table 1 shows the “distance to nearest reported CNV” = 0, for these regions), the 11 segments are not yet reported in the DGV (novel) (Table 1 shows the “distance to nearest reported CNV” > 0, for these regions), but these were within 82 kb distance from some other reported variants. Of the novel segments we found in the study that are associated with development of arsenic-induced skin lesions, 3 are from a known lincRNA RP11-76I14.1 (2q12.1), and 4 are from another known lincRNA CTC-535 M15.2 (5q34). The detail mappings are shown in Additional file 9: Figure S7.

Table 1 Cox regression analysis: Hazard Ratio (HR) for the development of arsenic-induced skin lesions by presence or absence of genomic deletions in different chromosomal locations. Combined analysis of total 2171 HEALS participants. Top genomic segments are sorted by chromosome

In the next step, in Cox regression models, in addition to the previous covariates, we also entered an interaction term “segment x gender” to find out if the deletion of any of the segment(s) affected the risk of skin lesion differently in male and female subjects. That also allowed us to identify segments significant in male and female subjects separately. Table 2 shows the segments that achieved statistical significance after Bonferroni correction for multiple testing for male subjects. For example, the association of 9p12 deletion covering intronic region of the gene ZNF658 and the development of arsenic-induced skin lesions is statistically significant in male subjects [HR 2.5, CI 1.7-3.7] and is stronger compared to that in female subjects [HR 1.3, CI 0.82-2.11] with interaction p = 0.03. This structural genomic variant has been reported in the DGV. Figure 2 shows the differential role of 15q13.3 deletion (APBA2 gene) in males and females for the development of arsenic-induced skin lesions. The detail mappings of these regions are shown in Additional file 10: Figure S8.

Table 2 Cox regression analysis: Hazard Ratio (HR) for the development of arsenic induced skin lesions by presence or absence of genomic deletions in different chromosomal locations in male subjects
Fig. 2
figure 2

Kaplan-Meier curves: x-axis represents time in months; deletion in chromosome 15q13.1 covering the gene APBA2 did not affect the development of arsenic-induced skin lesions in female subjects (p = 0.117, log rank test; shown on left), but significantly increased the risk in male subjects (p = 3.07E-05, log rank test; shown on right) exposed to arsenic through drinking water

Table 3 presents the segments with significant effect in female subjects. The association of 8p22 deletion covering intronic region of the SGCZ gene and the development of arsenic-induced skin lesion is stronger in female subjects [HR 2.4, CI 1.6-3.7] compared to that in males [HR 1.4, CI 1.02-1.96] with interaction p = 0.04. This structural genomic variant has been reported in the DGV.

Table 3 Cox regression analysis: Hazard Ratio (HR) for the development of arsenic induced skin lesions by presence or absence of genomic deletions in different chromosomal locations in female subjects

In the next step, in Cox regression analysis, in addition to gender, age, UACR, we also included an interaction term “segment x UACR” as predictors to find out if the deletion of any of the segment(s) affected the risk of skin lesion differently in subjects with high and low arsenic exposure (the Gene-Environment interaction). Results are presented in the Table 4. Most of these segments were in a lincRNA region. For example, among the group of individuals with high arsenic exposure (UACR > = median 192 μg/g of creatinine), those with deletion of 5q34 were at a 2.5 (CI 1.7-3.8) fold higher risk of skin lesion development compared to those without deletion. Among the group of individuals with low arsenic exposure, the corresponding HR was lower at 1.6 (CI 0.9-2.8). A similar effect was also seen for the deletion of chromosome 5q23.1 region (see Fig. 3). The detail mappings of these regions are shown in Additional file 11: Figure S9.

Table 4 Cox regression analysis: Hazard Ratio (HR) for the development of arsenic induced skin lesions by presence or absence of genomic deletions in different chromosomal locations in subjects with higher arsenic exposure (urinary arsenic creatinine ratio > = median value of 192 μg/g of creatinine)
Fig. 3
figure 3

Kaplan-Meier curves: x-axis represents time in months; deletion in chromosome 5q23.1 covering the lincRNA CTD-3179P9.1 did not affect the development of arsenic-induced skin lesions in subjects with lower UACR (< median 192 μg/g of creatinine, shown on left; p = 0.771, log rank test), but deletion of the same region significantly increased the risk in subjects with higher UACR (shown on right; p = 0.002, log rank test)

We tested if a deletion of any of these segments is associated with arsenic exposure. In that line, in logistic regression analysis for each of those segments, we used the segment (0: no deletion, 1: deletion) as the dependant variable; and for the independent variables we entered the measure of arsenic exposure (= < median vs. >median) along with gender (0: female, 1: male). Our data suggests that neither higher UACR, nor higher well water (as measure of intensity of arsenic exposure) were associated with higher prevalence of deletion for any of the segments. For many of the segments, however, deletion was more frequently found among the female subjects compared to the males (see Additional file 12: Table S3).

We also looked for if a copy number gain/amplification was associated with development of arsenic-induced skin lesions. Accordingly, we dichotomized the CN status of each segment as 0: no amplification and 1: amplification. In the Cox regression models, we calculated the HR for each amplification segment by entering the segment (0 vs. 1) as the predictor for arsenic-induced skin lesions, along with the covariates, gender (male vs. female), age (= < median 38 years vs. >median) and UACR (= < median vs. >median). None of the segments with amplification showed significantly higher risk for development of arsenic-induced skin lesions after Bonferroni correction for multiple testing.

Discussion

To our knowledge, this is the first large-scale genome-wide CN analysis to show from a prospectively followed-up cohort that the structural variation(s) in the germ line DNA may predispose an individual exposed to arsenic to develop arsenic-induced skin lesions. In tumor tissue, from arsenic-induced lung squamous cell carcinoma, deletion in chromosomal regions 1q21.1, 7p22.3, and 9q12 have been reported [38]. We were expecting to see some of the CNVs we found to overlap with those reported by Martinez et al. [38]. However, one of the explanations for the lack of overlap may be the fundamental difference between our study and the previous study. We have looked at CNV in blood DNA predisposing an individual to develop skin lesion. Our data does not suggest that the CNVs, we identified in the current study, were associated with arsenic exposure. On the other hand, the previous study focused at potentially “arsenic related” lung tissue specific and tumor-specific “somatic” CN change [38]. We have not yet examined the CN change in skin tissue from the arsenic-induced skin lesions biopsies. Hopefully in future we will be able to do that. But more importantly, as of now, we know very little about the functional significance of germline CNV; and have much to learn in chronic disease(s).

Previously, using case–control design, our group found some SNPs to be associated with arsenic metabolism [25]. Using clinical follow-up data in a larger cohort, our group also presented the evidence that higher risk of arsenic-induced skin lesion was found in male gender, increasing age and higher arsenic exposure [16, 27]. Now, we provide evidence that structural variation in the form of CN loss or deletion in certain genomic location(s) may have a role in the development of arsenic-induced skin lesions independent of gender, age, level of arsenic exposure and also independent of the SNPs related to arsenic metabolism. Our study is also the first to indirectly suggest the possible relationship between lincRNA and development of arsenic-induced skin lesions. The lincRNAs do not overlap exons of either protein-coding or other non-lincRNA types of genes. The role of lincRNAs are just being unveiled recently [43,44,45,46,47,48]. Thousands of lincRNAs are now known, however, many of their functions are still unknown [49].

Higher risk of development of arsenic-induced skin lesions among male individuals is reported by others as well [27, 50, 51]. As expected, higher UACR or higher well water arsenic content was related to higher risk of skin lesions. Reduction in arsenic exposure increases the odds that an individual with skin lesions would recover or show less severe lesions within 10 years [52]. Historically, in Bangladesh, the arsenic contamination in drinking water started after 1971 when the digging of deep wells started with an assumption that deep tube wells in rural areas would provide safe drinking water in terms of bacterial contamination [53]. Among the younger subjects (age = <median 38 yrs), the effect of higher UACR on development of skin lesions (HR2.6, CI 1.7-4.0) was stronger than its effect among the older subjects (HR1.39, CI 1.1-1.7). In most of the households, by culture, the male subjects (husbands) were older than the female subjects (wives), and male subjects usually consume higher volume of water than the female (3.8 L/day SD1.4 vs. 3.1 L/day SD 1.1, p = 3.3 × 10−35, in the present study). So even though they were consuming the water from the same source, perhaps the cumulative arsenic dosage was higher in males due to higher consumption and higher age (longer duration of exposure). But even after controlling for age, UACR, genetic markers (SNP and CNV), smoking and betel use, the male gender was strongly related to skin lesion development.

The structural variants, we are reporting in this paper to have significant association with the development of arsenic-induced skin lesions, have a frequency between 5% and 21% among the study population. However, for some of these variants that are previously reported in the DGV [31], the frequency was much different. For example, the deletion in 7q11.21 region: the deletion frequency ranges from 1 in 29,084 in a case–control study for developmental delay using arrays [54] to 35 in 2504 in a sequence-based study [32]; amplification/gain was reported as high as 121/270 [55]. These differences may be due to different population, different test platform, different disease entity of interest etc. The significant germ line CNVs we found in this study, to be related to arsenic-induced skin lesions, do not match with the CNVs in arsenic related lung cancer tissue [38].

To achieve genome-wide coverage, by design, the Illumina Cyto12 chip or most of the other commercially available SNP chips actually interrogate mostly SNPs in the inter-genic regions. Therefore, many of the GWAS hit SNPs are in fact far away from a gene. It is important to note the fact that, most of the statistically significant regions (20 out of 24) with deletion we found in this study, were overlapping or were within a gene. In other words, the list of significant deletion regions was enriched in genic regions. Only a few (4 out of 24) were located within inter-genic regions.

The utility of oligonucleotide based SNP chips for the detection of CN change and its interpretation in clinical samples was demonstrated by our group in the past [28]. In the past, we also validated the CN changes detected by similar SNP chips by comparing them with the results from a Luminex based multiplex assay for our other study (see Additional file 13: Figure S10). Our current study clearly suggests that there is some role of structural change in the genome (in the form of CN loss/deletion) in the development of arsenic-induced skin lesions, independent of the known clinical factors/parameters like age, sex, UACR level as well as the genotypes known to affect the arsenic metabolism.

One of the major strengths of this study is the long-term prospective follow-up and reasonably large sample size. The HRs of deletion(s) of these genomic segments for skin lesions may be slightly lower than the HRs of age, sex, but are not negligible. Over the last decade, we have been following up one of the largest cohorts exposed to arsenic [26] and depending on future availability of funding support, we have the opportunity to check the reproducibility of this novel finding in a larger and independent set of subjects.

We acknowledge the fact that the microarray platform used in this study (Cyto12) is not ultra-high density, the majority of the markers are intergenic and/or intronic in location and it does not have many markers in exonic regions. Functional characterization of the CNVs will be done in a future study. We did not have the source of RNA samples for these individuals to run gene expression to confirm the effects of deletion of lincRNA regions in the genomic DNA samples. For future studies, we will focus on this issue as well as tissue specificity.

Conclusion

Our genome wide CN analysis study of a prospectively followed-up cohort of arsenic exposure from drinking water suggests that individuals with CNVs in several genomic locations are predisposed to higher risk of development of arsenic-induced skin lesions. A few of these CNVs include lincRNA gene regions. The findings need to be replicated in another independent study and functional characterization would be needed to better understand the underlying genetic mechanism.

Abbreviations

ANOVA:

Analysis of variance

CI:

Confidence interval

CN:

Copy number

CNV:

Copy number variation

GWAS:

Genome wide association study

HEALS:

Health Effects of Arsenic Longitudinal Study

HR:

Hazard Ratio

lincRNA:

Long intervening noncoding RNA

LRR:

Log R ratio

QC:

Quality control

SD:

Standard deviation

SNP:

Single nucleotide polymorphism

UARC:

Urinary arsenic creatinine ratio

WAs:

Water arsenic

References

  1. Smith AH, Lingas EO, Rahman M. Contamination of drinking-water by arsenic in Bangladesh: a public health emergency. Bull World Health Organ. 2000;78(9):1093–103.

    CAS  Google Scholar 

  2. United States Environmental Protection Agency: Office of Water (4606) drinking water standard for arsenic. EPA 815-F-00-015. In.; 2001.

  3. Celik I, Gallicchio L, Boyd K, Lam TK, Matanoski G, Tao X, Shiels M, Hammond E, Chen L, Robinson KA, et al. Arsenic in drinking water and lung cancer: a systematic review. Environ Res. 2008;108(1):48–55.

    Article  CAS  Google Scholar 

  4. Liu J, Waalkes MP. Liver is a target of arsenic carcinogenesis. Toxicol Sci. 2008;105(1):24–32.

    Article  CAS  Google Scholar 

  5. Mink PJ, Alexander DD, Barraj LM, Kelsh MA, Tsuji JS. Low-level arsenic exposure in drinking water and bladder cancer: a review and meta-analysis. Regul Toxicol Pharmacol. 2008;52(3):299–310.

    Article  CAS  Google Scholar 

  6. Steinmaus CM, Ferreccio C, Romo JA, Yuan Y, Cortes S, Marshall G, Moore LE, Balmes JR, Liaw J, Golden T, et al. Drinking water arsenic in northern chile: high cancer risks 40 years after exposure cessation. Cancer Epidemiol Biomark Prev. 2013;22(4):623–30.

    Article  CAS  Google Scholar 

  7. Chen CJ, Chen CW, Wu MM, Kuo TL. Cancer potential in liver, lung, bladder and kidney due to ingested inorganic arsenic in drinking water. Br J Cancer. 1992;66(5):888–92.

    Article  CAS  Google Scholar 

  8. Yuan Y, Marshall G, Ferreccio C, Steinmaus C, Liaw J, Bates M, Smith AH. Kidney cancer mortality: fifty-year latency patterns related to arsenic exposure. Epidemiology. 2010;21(1):103–8.

    Article  Google Scholar 

  9. Vahidnia A, van der Voet GB, de Wolff FA. Arsenic neurotoxicity--a review. Hum Exp Toxicol. 2007;26(10):823–32.

    Article  CAS  Google Scholar 

  10. Martin E, Gonzalez-Horta C, Rager J, Bailey KA, Sanchez-Ramirez B, Ballinas-Casarrubias L, Ishida MC, Gutierrez-Torres DS, Hernandez Ceron R, Viniegra Morales D, et al. Metabolomic characteristics of arsenic-associated diabetes in a prospective cohort in Chihuahua, Mexico. Toxicol Sci. 2015;144(2):338–46.

    Article  CAS  Google Scholar 

  11. Hall EM, Acevedo J, Lopez FG, Cortes S, Ferreccio C, Smith AH, Steinmaus CM. Hypertension among adults exposed to drinking water arsenic in northern Chile. Environ Res. 2017;153:99–105.

    Article  CAS  Google Scholar 

  12. Jiang J, Liu M, Parvez F, Wang B, Wu F, Eunus M, Bangalore S, Newman JD, Ahmed A, Islam T, et al. Association between arsenic exposure from drinking water and Longitudinal change in blood pressure among HEALS cohort participants. Environ Health Perspect. 2015;123(8):806–12.

    Google Scholar 

  13. Moon KA, Guallar E, Umans JG, Devereux RB, Best LG, Francesconi KA, Goessler W, Pollak J, Silbergeld EK, Howard BV, et al. Association between exposure to low to moderate arsenic levels and incident cardiovascular disease. A prospective cohort study. Ann Intern Med. 2013;159(10):649–59.

    Google Scholar 

  14. Tsuji JS, Perez V, Garry MR, Alexander DD. Association of low-level arsenic exposure in drinking water with cardiovascular disease: a systematic review and risk assessment. Toxicology. 2014;323:78–94.

    Article  CAS  Google Scholar 

  15. Wu F, Jasmine F, Kibriya MG, Liu M, Cheng X, Parvez F, Islam T, Ahmed A, Rakibuz-Zaman M, Jiang J, et al. Interaction between arsenic exposure from drinking water and genetic polymorphisms on cardiovascular disease in Bangladesh: a prospective case-cohort study. Environ Health Perspect. 2015;123(5):451–7.

    CAS  Google Scholar 

  16. Ahsan H, Chen Y, Parvez F, Zablotska L, Argos M, Hussain I, Momotaj H, Levy D, Cheng Z, Slavkovich V, et al. Arsenic exposure from drinking water and risk of premalignant skin lesions in Bangladesh: baseline results from the health effects of arsenic Longitudinal study. Am J Epidemiol. 2006;163(12):1138–48.

    Article  Google Scholar 

  17. Flanagan SV, Johnston RB, Zheng Y. Arsenic in tube well water in Bangladesh: health and economic impacts and implications for arsenic mitigation. Bull World Health Organ. 2012;90(11):839–46.

    Article  Google Scholar 

  18. Karagas MR, Gossai A, Pierce B, Ahsan H. Drinking water arsenic contamination, skin lesions, and malignancies: a systematic review of the global evidence. Curr Environ Health Rep. 2015;2(1):52–68.

    Article  CAS  Google Scholar 

  19. Mayer JE, Goldman RH. Arsenic and skin cancer in the USA: the current evidence regarding arsenic-contaminated drinking water. Int J Dermatol. 2016;55(11):e585–91.

    Article  CAS  Google Scholar 

  20. Kile ML, Rodrigues EG, Mazumdar M, Dobson CB, Diao N, Golam M, Quamruzzaman Q, Rahman M, Christiani DC. A prospective cohort study of the association between drinking water arsenic exposure and self-reported maternal health symptoms during pregnancy in Bangladesh. Environ Health. 2014;13(1):29.

    Article  Google Scholar 

  21. Argos M, Kalra T, Rathouz PJ, Chen Y, Pierce B, Parvez F, Islam T, Ahmed A, Rakibuz-Zaman M, Hasan R, et al. Arsenic exposure from drinking water, and all-cause and chronic-disease mortalities in Bangladesh (HEALS): a prospective cohort study. Lancet. 2010;376(9737):252–8.

    Article  CAS  Google Scholar 

  22. Haque R, Mazumder DN, Samanta S, Ghosh N, Kalman D, Smith MM, Mitra S, Santra A, Lahiri S, Das S, et al. Arsenic in drinking water and skin lesions: dose–response data from West Bengal, India. Epidemiology. 2003;14(2):174–82.

    Google Scholar 

  23. NationalResearchCouncilSubcommittee. In: Arsenic in drinking water: 2001 update. Washington DC: 2001 by the National Academy of Sciences; 2001.

  24. Smith AH, Steinmaus CM. Health effects of arsenic and chromium in drinking water: recent human findings. Annu Rev Public Health. 2009;30:107–22.

    Article  Google Scholar 

  25. Pierce BL, Kibriya MG, Tong L, Jasmine F, Argos M, Roy S, Paul-Brutus R, Rahaman R, Rakibuz-Zaman M, Parvez F, et al. Genome-wide association study identifies chromosome 10q24.32 variants associated with arsenic metabolism and toxicity phenotypes in Bangladesh. PLoS Genet. 2012;8(2):e1002522.

    Article  CAS  Google Scholar 

  26. Ahsan H, Chen Y, Parvez F, Argos M, Hussain AI, Momotaj H, Levy D, van Geen A, Howe G, Graziano J. Health effects of arsenic Longitudinal study (HEALS): description of a multidisciplinary epidemiologic investigation. J Expo Sci Environ Epidemiol. 2006;16(2):191–205.

    Article  CAS  Google Scholar 

  27. Argos M, Kalra T, Pierce BL, Chen Y, Parvez F, Islam T, Ahmed A, Hasan R, Hasan K, Sarwar G, et al. A prospective study of arsenic exposure from drinking water and incidence of skin lesions in Bangladesh. Am J Epidemiol. 2011;174(2):185–94.

    Article  Google Scholar 

  28. Jasmine F, Rahaman R, Dodsworth C, Roy S, Paul R, Raza M, Paul-Brutus R, Kamal M, Ahsan H, Kibriya MG. A genome-wide study of cytogenetic changes in colorectal cancer using SNP microarrays: opportunities for future personalized treatment. PLoS One. 2012;7(2):e31968.

    Article  CAS  Google Scholar 

  29. Zhang J, Feuk L, Duggan GE, Khaja R, Scherer SW. Development of bioinformatics resources for display and analysis of copy number and other structural variants in the human genome. Cytogenet Genome Res. 2006;115(3–4):205–14.

    Article  CAS  Google Scholar 

  30. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C. Detection of large-scale variation in the human genome. Nat Genet. 2004;36(9):949–51.

    Article  CAS  Google Scholar 

  31. MacDonald JR, Ziman R, Yuen RK, Feuk L, Scherer SW. The database of genomic variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 2014;42(Database issue):D986–92.

    Article  CAS  Google Scholar 

  32. Mills RE, Walter K, Stewart C, Handsaker RE, Chen K, Alkan C, Abyzov A, Yoon SC, Ye K, Cheetham RK, et al. Mapping copy number variation by population-scale genome sequencing. Nature. 2011;470(7332):59–65.

    Article  CAS  Google Scholar 

  33. Pinto D, Pagnamenta AT, Klei L, Anney R, Merico D, Regan R, Conroy J, Magalhaes TR, Correia C, Abrahams BS, et al. Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010;466(7304):368–72.

    Article  CAS  Google Scholar 

  34. McCarthy SE, Makarov V, Kirov G, Addington AM, McClellan J, Yoon S, Perkins DO, Dickel DE, Kusenda M, Krastoshevsky O, et al. Microduplications of 16p11.2 are associated with schizophrenia. Nat Genet. 2009;41(11):1223–7.

    Article  CAS  Google Scholar 

  35. Stefansson H, Rujescu D, Cichon S, Pietilainen OP, Ingason A, Steinberg S, Fossdal R, Sigurdsson E, Sigmundsson T, Buizer-Voskamp JE, et al. Large recurrent microdeletions associated with schizophrenia. Nature. 2008;455(7210):232–6.

    Article  CAS  Google Scholar 

  36. Craddock N, Hurles ME, Cardin N, Pearson RD, Plagnol V, Robson S, Vukcevic D, Barnes C, Conrad DF, Giannoulatou E, et al. Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls. Nature. 2010;464(7289):713–20.

    Article  CAS  Google Scholar 

  37. McCarroll SA, Huett A, Kuballa P, Chilewski SD, Landry A, Goyette P, Zody MC, Hall JL, Brant SR, Cho JH, et al. Deletion polymorphism upstream of IRGM associated with altered IRGM expression and Crohn’s disease. Nat Genet. 2008;40(9):1107–12.

    Article  CAS  Google Scholar 

  38. Martinez VD, Buys TP, Adonis M, Benitez H, Gallegos I, Lam S, Lam WL, Gil L. Arsenic-related DNA copy-number alterations in lung squamous cell carcinomas. Br J Cancer. 2010;103(8):1277–83.

    Article  CAS  Google Scholar 

  39. Cheng Z, Zheng Y, Mortlock R, Van Geen A. Rapid multi-element analysis of groundwater by high-resolution inductively coupled plasma mass spectrometry. Anal Bioanal Chem. 2004;379(3):512–8.

    Article  CAS  Google Scholar 

  40. Nixon DE, Mussmann GV, Eckdahl SJ, Moyer TP. Total arsenic in urine: palladium-persulfate vs nickel as a matrix modifier for graphite furnace atomic absorption spectrophotometry. Clin Chem. 1991;37(9):1575–9.

    CAS  Google Scholar 

  41. Heinegard D, Tiderstrom G. Determination of serum creatinine by a direct colorimetric method. Clin Chim Acta. 1973;43(3):305–10.

    Article  CAS  Google Scholar 

  42. Downey T. Analysis of a multifactor microarray study using Partek genomics solution. Methods Enzymol. 2006;411:256–70.

    Article  CAS  Google Scholar 

  43. Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F, et al. Landscape of transcription in human cells. Nature. 2012;489(7414):101–8.

    Article  CAS  Google Scholar 

  44. Esteller M. Non-coding RNAs in human disease. Nat Rev Genet. 2011;12(12):861–74.

    Article  CAS  Google Scholar 

  45. Huarte M. The emerging role of lncRNAs in cancer. Nat Med. 2015;21(11):1253–61.

    Article  CAS  Google Scholar 

  46. Ling H, Vincent K, Pichler M, Fodde R, Berindan-Neagoe I, Slack FJ, Calin GA. Junk DNA and the long non-coding RNA twist in cancer genetics. Oncogene. 2015;34(39):5003–11.

    Article  CAS  Google Scholar 

  47. Ulitsky I, Bartel DP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154(1):26–46.

    Article  CAS  Google Scholar 

  48. Xiong XD, Ren X, Cai MY, Yang JW, Liu X, Yang JM. Long non-coding RNAs: an emerging powerhouse in the battle between life and death of tumor cells. Drug Resist Updat. 2016;26:28–42.

    Article  Google Scholar 

  49. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.

    Article  Google Scholar 

  50. Guha Mazumder DN, Haque R, Ghosh N, De BK, Santra A, Chakraborty D, Smith AH. Arsenic levels in drinking water and the prevalence of skin lesions in West Bengal, India. Int J Epidemiol. 1998;27(5):871–7.

    Article  CAS  Google Scholar 

  51. Tondel M, Rahman M, Magnuson A, Chowdhury IA, Faruquee MH, Ahmad SA. The relationship of arsenic levels in drinking water and the prevalence rate of skin lesions in Bangladesh. Environ Health Perspect. 1999;107(9):727–9.

    Article  CAS  Google Scholar 

  52. Seow WJ, Pan WC, Kile ML, Baccarelli AA, Quamruzzaman Q, Rahman M, Mahiuddin G, Mostofa G, Lin X, Christiani DC. Arsenic reduction in drinking water and improvement in skin lesions: a follow-up study in Bangladesh. Environ Health Perspect. 2012;120(12):1733–8.

    CAS  Google Scholar 

  53. Yunus FM, Khan S, Chowdhury P, Milton AH, Hussain S, Rahman M. A review of groundwater arsenic contamination in Bangladesh: the millennium development goal era and beyond. Int J Environ Res Public Health. 2016;13(2):215.

    Article  Google Scholar 

  54. Coe BP, Witherspoon K, Rosenfeld JA, van Bon BW, Vulto-van Silfhout AT, Bosco P, Friend KL, Baker C, Buono S, Vissers LE, et al. Refining analyses of copy number variation identifies specific genes associated with developmental delay. Nat Genet. 2014;46(10):1063–71.

    Article  CAS  Google Scholar 

  55. Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, et al. Global variation in copy number in the human genome. Nature. 2006;444(7118):444–54.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors express their thanks to all the participants of the Health Effects of Arsenic Longitudinal Study (HEALS). Authors also thank all the study physicians and health workers for their dedicated works that made the study a success.

Funding

This work was supported by the National Institutes of Health grants U01 CA122171, P30 CA 014599, P42ES010349, R01CA102484, and R01CA107431. However the freedoms to design, conduct, interpret, and publish research is not compromised by any controlling sponsor.

Availability of data and materials

Data explaining the results is provided as additional files. The genetic dataset is available to bona fide researchers through reasonable request to the PI of the grant – Prof. Habibul Ahsan (habib@uchicago.edu) for the purpose of confirming the study results presented in this paper.

Author information

Authors and Affiliations

Authors

Contributions

MGK, FJ, HA made contributions to conception, design, data interpretation, statistical analysis, and drafting. FP, MA, AA, TI, MRZ, JHG made contribution for epidemiological data collection and analysis. FJ, SR, RPB, JS, VS contributed for data acquisitions. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Muhammad G. Kibriya.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Institutional Review board of The University of Chicago, Columbia University, and the Bangladesh Medical Research Council. Informed consent was obtained from all participants.

Consent for publication

Not applicable.

Competing interests

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

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Characteristics of HEALS participants selected for this study. None of them had arsenic induced skin lesion at baseline. (XLS 21 kb)

Additional file 2: Figure S1.

Correlation of log2 transformed urinary arsenic creatinine ratio (UACR) and log2 transformed well water arsenic (WAs) concentration (r = 0.66, p = 3.7 E-280). (PPT 108 kb)

Additional file 3: Figure S2.

QC of samples by standard deviation of Log R Ratio (LRR) and by SNP call rate of the array. (PPT 163 kb)

Additional file 4: Figure S3.

Kaplan-Meir plots show (a) male subjects were at higher risk of developing arsenic-induced skin lesion than the female subjects (p = 4.5 E-35, log rank test; shown on left side) exposed to arsenic through drinking water; and (b) subjects with higher age (>median 38 years) were also at higher risk for development of skin lesion (p = 3.8 E-40; shown on right side) than those who were younger. X-axis represents time to event (months of follow-up after enrollment). (PPT 76 kb)

Additional file 5: Figure S4.

Kaplan-Meir plots show (a) subjects with higher UACR (>median 192 μg/g of creatinine) were at higher risk of developing arsenic-induced skin lesion than those with lower UACR (p = 0.001, log rank test; shown on left side); (b) categorization by well water arsenic (WAs) also showed similar effect – higher risk in subjects drinking water with higher arsenic concentration (>median 56 μg/L) compared to those drinking water with lower arsenic concentration (p = 1.97 E-09, log rank test; shown on right side). X-axis represents time to event (months of follow-up after enrollment). (PPT 77 kb)

Additional file 6: Table S2.

Results from all the 1135 Cox regression models with HR and 95%CI of all the variables (genomic segment, gender, age and UACR). (XLS 1116 kb)

Additional file 7: Figure S5.

Conditioning on Arsenic metabolism SNPs did not show any effect on the HR of the genomic segments (adjusted for gender, age & UACR). (PPT 138 kb)

Additional file 8: Figure S6.

Adjusting for smoking habit or use of betel leaf did not show any effect on the HR of the genomic segments (adjusted for gender, age & UACR). (PPT 141 kb)

Additional file 9: Figure S7.

Detail mapping of the CNV regions (presented in Table 1), predisposing to significantly higher risk for development of arsenic-induced skin lesions. (PDF 201 kb)

Additional file 10: Figure S8.

Detail mapping of the CNV regions showing interaction with gender (presented in Tables 2 and 3) for higher risk of developing arsenic-induced skin lesions. (PDF 234 kb)

Additional file 11: Figure S9.

Detail mapping of the CNV regions showing interaction with higher arsenic exposure (presented in Table 4) for higher risk of developing arsenic-induced skin lesions. (PDF 239 kb)

Additional file 12: Table S3.

Frequency of the significant deleted segments in female and male subjects. (XLS 24 kb)

Additional file 13: Figure S10.

Correlation between copy number of a genomic region derived from oligonucleotide SNP chip and fluorescent intensity derived from Luminex based assay. (PPT 159 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kibriya, M.G., Jasmine, F., Parvez, F. et al. Association between genome-wide copy number variation and arsenic-induced skin lesions: a prospective study. Environ Health 16, 75 (2017). https://doi.org/10.1186/s12940-017-0283-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12940-017-0283-8

Keywords