Chemical variability of artificial stone powders in relation to their health effects

The occurrence of highly severe silica-related diseases among the resin- and silica-based artificial stone workers was claimed, associated to an extremely short latency. High levels of exposure and intrinsic properties of AS are thought to modulate the development of silicosis and auto-immune diseases. This study compares parent materials and processed dusts, to shed light on changes of AS occurring in the manufacturing process, through an XRF, EPR and XAS investigation. We point out the extremely wide variability of the materials, the occurrence of chemical signatures impressed by the processing techniques, and the unprecedented generation of stable radicals associated to the lysis of the Si-O chemical bond inside the resin coated respirable crystalline silica. These results suggest that the AS processing in industrial stone workshops can create respirable dusts with peculiar physical and chemical properties, to be correlated to the observed clinical evidences.

. Number of cases, maximum, minimum and median at% values from XRF analyses.
Na, Mg, Al, K, Ca, Ti are frequently found, while Cl and Fe show seldom occurrences. Na, Mg, K and Ca exhibit ranges fully covering both the trace and minor element fields, whereas P, S, Cl, Fe, Co, Cu and Zr are trace contaminants. The elements not shown in Fig. 1 never exceed the field of the minor elements, being mostly trace contaminants. Zn occurs as minor component only in three samples (7A, 7B and 7C). Moreover, data concerning some other elements (Na, Mg, Ca and Ti, and, partly, P, S, Co and Cu) suggest that the chemistry of these three samples is markedly different from that of all other analysed samples (Fig. 1). Since their anomalous behaviour can affect the join multivariate structure of the data 26,27 the subsequent multivariate analysis was performed without including them in the database (SI, Section E). The dendrogram (Fig. 2) points to two main groups of samples. A and B samples cluster in different groups (with the exception of the 1A sample), whereas C samples do not exhibit preferential grouping (with A or B) nor group alone. Considering that C samples are parent of both A and B, these latter being linked to different stone processing, the dendrogram points to the occurrence of a certain "chemical signature" impressed by the processing. The cluster analysis was performed on log-centred data by using the squared Euclidean distance as similarity measure and the Ward method to link cases 28 , (SI section C).
A further step in the analysis of the dataset variability has been carried out using the biplot methodology for compositional data 24,29 (SI, sec. C). Having the longest rays (Fig. 3), Ca, Na, Co and Ti explain most of the variability of the XRF dataset, with respect to the compositional barycentre. These rays define quadrants where samples of different groups are well discriminated. In particular, the rays of Na and Ca (plus Mg, Cl, Fe) head towards samples of groups A, whereas rays of Co (plus Cu, Zr, Cr, Al, K, Fe) towards group B and Ti (plus Si, P, S, Cl) towards group C. Indeed, the A, B and C groups appear characterised by a different compositional signature. Closer detail on the relationships among elements in the biplot has been provided by opportune sets of 3 variables (sub-compositions), representable in a ternary diagram (SI, section F). In the following, two specific  www.nature.com/scientificreports www.nature.com/scientificreports/ hypotheses about the origin of contamination of samples during processing have been tested exploiting specific ternary diagrams: • Contamination occurring at the interface between water and CS (wet processing) • Contamination occurring at the interface between CS and machine-tool materials (dry processing) The first hypothesis was tested through the study of the subcomposition Ca-Fe-Si. The ternary diagram (Fig. 4a) points to an increase in Ca associated with an increase in the ratio Fe/Si, as shown by the PC1 direction, which explains ~90% of the data variability. Accordingly, the internal relationship among these three variables appears a feature of the whole system. Thus, the diagram of Fig. 4a supports the fact that a part of Fe is added to the system together with Ca by wet processing.
Concerning the second hypothesis, the dry process may result in a sample contamination accounting for the relative hardness of the minerals/materials involved. Under this scheme, Ti is contributed by rutile, Fe by tools and Si by CS. The diagram of the Si-Ti-Fe sub-composition (Fig. 4b), the PC1 of which explains ~83% of the data variability, points to a light increase in Si associated with a high variability in the Fe/Ti ratio. These results indicate an alternative presence of either Ti or Fe, in agreement with the hardnesses of rutile and steel, respectively. The higher the rutile content, at the expenses of the CS content, the softer is the material, and thus the Fe content decreases. In order to assess to which extent Fe can be considered a proxy of the contamination by working-tools, we used an indirect method, considering Co as an element of special steel. The ternary diagram of the sub-composition Si-Ti-Co (Fig. 4c), with PC1 capturing the 82% of the data variability, reveals a close similarity to that of Fig. 4b. Thus, we can safely assume that at least a part of the Fe on the processed samples is provided as an impurity during the processing step, and that the amounts provided by the wet and dry processes may not be characterised by an equivalent speciation.
Na is one of the most variable elements in the compositional dataset, suggesting possible contributions by different sources (Fig. 3). Figure 4d depicts the ternary diagram of the Na-Al-K field, chosen to assess if the association between these three elements (commonly found associated in alkaline feldspars) was occurring. The observed trend confirms that Al and K, whose ratio is almost constant, can be attributed, as hypothesized, to feldspars, contained in the parent materials (in agreement with the XRD findings). Conversely, Na is provided by at least two sources, the most relevant of which is not linked to the parent materials, and it is likely associated to the wet processing.
EPR spectroscopy: transition metal ions. The intrinsic heterogeneity of the samples, and especially the variability of the chemical composition, has a counterpart in the results of the EPR investigations. Considering only the analytically detected elements, at least Al, Si, Ti, Cr, Mn, Fe, Co, Cu can occur as EPR active species (transition metal ions, TMI, and/or inorganic radicals). The spectral variability, observed in the panoramic EPR spectra of all 21 samples (in SI, Section G), can be ascribed to this circumstance. However, most of the main experimental features can be attributed to the Fe speciation. Signals from three different Fe species have been observed: (a) Fe(III) as isolated ion, occurring in a rhombic coordination, revealed by the typical narrow EPR  Fig. 5b: in this case, the extremely broad width, and the occurrence of a zero-field absorption, unambiguously indicate the excitation of magnetic resonant modes 30 ; (c) Fe occurring in superparamagnetic particles, whose crystal size is small enough to behave as single-domain magnetic particles, as in the case of Fig. 5c; these species are characterised by moderately broad signals (ΔH < 200 mT), whose centre changes towards lower magnetic field values with decreasing temperature 30 . Only the EPR signal of Mn(II), among other TMI possibly occurring in the samples, has been frequently observed (Fig. 5d). This species, with its isotropic hyperfine structure, is likely to occur in a carbonate environment 31 . The overall findings of the EPR investigation are summarised in the Table 3.
Certain Fe species apparently show a general trend over the set of investigated samples. The most striking trend concerns the permanent magnetic species. These, in fact, are detected in A and B, but not in C samples. This suggests a link with the treatments undergone by the A and B samples: worked materials are contaminated by magnetic Fe species during the cutting processes. Conversely, no evidence of a significant discrimination among the A and B samples was observed. Signals attributed to superparamagnetic Fe species are observed in all C samples, and in most of the treated A and B samples, without evidence of specific trends. However, the data in Table 3 suggest that the superparamagnetic species, when abundantly present in the C samples, are always found also in the treated samples. This fact could be interpreted considering these species as isolated from the silica framework, and, being softer 32 , transferred to the treated dusts. Fe(III) ions are only seldom observed. Mn(II) is never detected in the raw C samples, while being preferentially found in the wet A samples. The fact that A samples are characterized by a robust enrichment in Ca(II) (see §3.2) and the attribution of the Mn(II) spectrum to calcium carbonate, point to a definite relationship between these two elements. We thus infer that Mn(II) contamination occurs, during the wet cutting, through the water used in abating procedures. EPR spectroscopy: radical species. All samples reveal a noticeable spectrum due to the h Al radical 33 .
This spectrum consists of two multiplets, centred at B = 337 mT and at B = 346 mT, due to the superhyperfine interaction of the unpaired electron, located onto an O − anion, with neighbouring H and Al nuclei (Fig. 6a). This radiogenic radical, very common in quartz, is located inside the crystals, with the proton hosted in the channels parallel to [001] axis in the α-quartz structure. The presence of the h Al spectrum is observed almost independently on the type of process followed by the investigated A and B samples. Thus, it cannot be considered as an efficient proxy of changes in the radical speciation during AS processing. It is worth mentioning that no h Al signal is detected in the samples of the series 4, in full agreement with the almost complete absence of quartz ( §3.1 and SI, Section H).
All A and B samples, but those of the series 4, present the ubiquitary evidence of a further radical species, hereafter called as R, superimposed to the h Al spectrum. Its presence, however, is never observed in the #C samples. Thus, the R radical can be considered as a product of the AS processing. From the detailed EPR spectra (e.g. Fig. 6a), only an approximate definition of its parameters can be gained, i.e. its position, corresponding to a value of g = 2.0032(2) and its very narrow width, of the order of 0.15 mT. No evidence of hyperfine structure of www.nature.com/scientificreports www.nature.com/scientificreports/ the signal is observed. A closer detail of the R species was obtained through the EEPR investigation: from this technique, the experimental spectra of both the 5B and 5C samples reveal, as previously stated, the presence of the asymmetric spectrum of the h Al radical (the hyperfine structure of this species is unresolved due to the used echo sequence 34 , Fig. 6b). The spectrum of the 5B sample presents a further signal, due to the R species. The simultaneous presence of both spectra points to the attribution of the R radical to a species with the same S = 1/2 spin state as the h Al centre. A better visualization of the signal due to the R species was obtained by subtracting the spectrum of the 5C sample from that of the 5B one (Fig. 6b). The residual signal intensity is thus fully attributable to the R species. A spectral fit carried out assuming a Lorentzian line shape indicates that the spectrum is characterized by a slight anisotropy in the g values (2.0043(2) and 2.0032(2) for the parallel and perpendicular components, respectively) associated to a larger difference in the line width (1.04(3) and 0.20(1) mT, respectively). The absence of hyperfine interaction was confirmed. On this basis, we can tentatively attribute this signal to an unpaired electron not chemically bonded to magnetic nuclei, with g values compatible with either an inorganic radical as e.g. the E centre in quartz (i.e. an unpaired electron localized on a tri-coordinated Si atom) 35 , or to an organic radical, likely generated during the thermal and mechanical of sample processing.
The radical speciation of the samples belonging to the series 4 appears to be different. In the 4C sample, a very broad and undefined signal at g ≈ 2 (labelled as "D") is registered. However, the very bad signal-to-noise ratio prevents any kind of more accurate characterization of its position and shape. Despite the bad spectral quality, the EEPR spectrum of the D species (Fig. 6c) reveals a slightly broad signal, centred at g = 2.0040(2), and a line

Magnetic Fe species
Superparam. Fe species Fe(III) Other TMI www.nature.com/scientificreports www.nature.com/scientificreports/ width of 1.9(1) mT. The very poor signal-to-noise ratio allows its attribution to a radical species, originally present in the polyester resin; the absence of h Al radicals is in fact linked to the absence of quartz. This attribution is also in agreement with the fact that cristobalite, a high temperature polymorph of silica, crystallizes without intrinsic radicals 35 . When sample 4C undergoes processing, EPR (Fig. 6d) and EEPR (Fig. 6e) indicate that a new radical species is formed (labelled as "X"). This species is observed on both 4A and 4B samples. In fact, the EEPR of Fig. 6e shows that the sharp line of the species X is overlapped on the broader line of the D species. The X line position is consistent with a g = 2.0030(2) value, very similar to that of the R radical. However, the absence of any appreciable Zeeman anisotropy is supporting the attribution to a different species. The line width of the spectrum of the radical X is 0.45 (2) mT. This result suggests further considerations. As far as the radical speciation is concerned, the quartz-bearing and quartz-free samples subjected to processing exhibit different results, so the final radical speciation is closely linked to the mineralogical composition of the original material. In contrast, no, or very weak, apparent relationships are found with the sample chemistry (i.e. with respect to the chemical signatures described §3.2) or in the type (wet, dry) of processing. Moreover, the resin in its original state may contain some radicals, but these are very few, especially if compared with h Al , whose concentration has been estimated in 5*10 14 defect/mol 33 , and they seem substantially unaffected by the sample processing. Thus, we attribute the newly formed R and X species to the interaction between the surfaces of mineral/resins and of processing tools. Apparently, all A and B samples but those of the series #4 present a common FT pattern (Fig. 7a): they mark the presence of different groups of peaks. Those at ~7-9, ~11-12 and ~23-27 MHz are due to the interaction of the h Al radical with the neighbouring H nuclei, within the channel of the quartz structure 33 . Conversely, the additional peak at ~14-15 MHz can be attributed to the Larmor resonance of free H nuclei, i.e. to several nuclei that weakly interact with the paramagnetic centre through dipolar interaction. This signal is attributed to the R species. If hosted inside the quartz crystal, or inside the resin, the R species would exhibit strong dipolar and probably isotropic hyperfine interactions: thus, we can reasonably state that the R species is located at the interface between the resin and the crystal, i.e. it is a surface species, of the type of the E centres 35 .
Concerning the FT spectra of the 4A and 4B samples, no appreciable nuclear (proton) modulation patterns are observed (Fig. 7b). Accordingly, also the X radical does not belong to the resin. However, no information about the local surrounding of the X species was gained. Thus, considering the main species present in the 4B and 4C samples, i.e. cristobalite, rutile, and resin, we can exclude only the last one. The two tentative attributions can be surface radicals of either cristobalite, or rutile. One has to consider, that for this second attribution, the weak peaks in the FT pattern of Fig. 7b, occurring at 0.35 MHz and 1.30 MHz could be attributed to hyperfine interaction with 47 Ti and 49 Ti isotopes, with quadrupolar contribution. X-ray absorption spectroscopy. The analysis of the XANES and EXAFS regions highlights that Fe is present in multiple species in all the studied samples. The co-presence of different oxidized phases, together with metallic Fe, was therefore always taken into account. Linear combination fits (LCF) of the XANES region were performed employing several Fe-bearing minerals as standards (Section J in SI). The ubiquitous presence of Fe in multiple phases does not allow one to identify unambiguously all the host components and LCF results should be regarded as merely indicative of the "kind" of Fe-bearing phases.
Data from LCF of the XANES were then used as a starting point to sort out the most probable components contributing to the overall experimental EXAFS spectra. The quantitative analysis of the EXAFS data was carried www.nature.com/scientificreports www.nature.com/scientificreports/ out considering the possible co-presence of metallic Fe together with other oxidized phases. Whenever the presence of metallic Fe was highlighted by the LCF results, fits were performed by initially assigning two different amplitude factors to the "Fe-oxide" phase and metallic Fe. After the achievement of a suitable model for the oxidized phase, a single amplitude factor was used. The metallic and oxidized phases were then related to the overall amplitude by means of a proportional factor and their sum constrained to be 1, thus obtaining an estimation of the ratio between metallic and total Fe content in a manner similar to that described by Di Benedetto et al. 36 ( Table 4).
The EXAFS signal of the samples #6A, shown as example in Fig. 8a,b, is dominated by an oxidized component; however, the comparison with the metallic Fe spectrum suggests the presence of moderate amounts of this species. Indeed, the quantitative EXAFS analysis yielded a metallic Fe fraction of about 30(10) % of the total Fe content with a double distance 2.48(3)-2.88(3) Å for metallic Fe, typical of the bulk bcc metal. Figure 8c shows the metallic Fe fraction thus obtained by the EXAFS analysis on all the studied samples; such values are in good agreement, especially if considering the complexity of the studied system, with those obtained by the LCF analysis. The most evident feature is that none of the raw samples shows the presence of metallic Fe (Fig. 8c), while the mechanically processed samples, both in wet and dry conditions, show in almost every series the presence of variable amounts of metallic Fe, with contents reaching up to the 60% of the total Fe content.
As for the oxidized portion of the spectrum, Fig. 8d shows the results obtained on the I-shell bond distance; in a few cases, the EXAFS analysis of the raw samples required the use of multiple paths in order to properly model the data. Second shell contributions are inconclusive since the presence of multiple Fe-bearing phases makes unambiguous information hard to infer. For the same reason we are prevented from getting information on the host phases by the investigation of the first shell only: indeed, Fe could be present in at least 4 different coordination environments, other than that pertaining to its metallic form, since the presence of both tetrahedral and octahedral Fe at both 2+ and 3+ oxidation states, cannot be excluded. However, the obtained data allow us to infer some information about the effects of the mechanical treatments. The Fe-O bond length of the treated samples rarely shows distances similar to those belonging to the raw fragments (Fig. 8d). In general, average Fe-O bond distances in C samples seem to be shorter than those of the A and B samples. The above-mentioned modifications can be ascribed to both a change in the Fe 2+ /Fe 3+ ratio and to a variation of Fe coordination; unfortunately, in such complex systems, it is not possible to infer more detailed information about the effects of the process.

Discussion
The results highlight an extreme variability of the trace elements composition, already attributed to an intrinsic consequence of the stone assemblages from CS, pigments and resins 19 . Nevertheless, PCA was able to sort out internal trends in the dataset structure, pointing out that: (1) A specific chemical contamination of the parent AS occurs depending on the wet or dry processing operated.
(2) The dry polishing of AS induces a contamination determined by the working-tools composition (special steels) 37 ; (3) The wet processed materials include several elements compatible with the ion content of the water; namely, the unexpectedly large Ca content observed in almost all the investigated wet samples is certainly linked to the water chemical composition. In fact, the wet processing involves recycled water from a closed circuit in all considered firms. The chemical signature of water is driven by the continuous contact with the abated sediment (which is only seldom removed).
The processing of AS results in detectable changes of the Fe speciation. Namely, the EXAFS analysis proves that, during both wet and dry processing, considerable amounts of metallic Fe (likely provided by working-tools) are added to the studied materials, identifiable by EPR as permanent magnetic phases. Moreover, by delving into the Fe-O bond distance, it is possible to notice that the modifications induced by the mechanical treatments extend to all the detected Fe forms by modifying Fe coordination and/or Fe oxidation state.
The process of dust contamination verified by both XRF and XAS analyses has to be considered as extremely relevant for toxicity studies: Pavan et al. 19 , in fact, have already pointed out that the presence of numerous redox active transition species coupled to RCS during respiration can have a relevant role in modulating (and especially improving) reactivity in HO˙ generation. The present data confirm that hetero-ions are provided to the final respirable dusts by the processing itself, and that difference can also occur between the wet and dry processing. Accordingly, this evidence can be linked to the high and variable toxicity of these dusts.
The AS processing is also effective in changing the radical speciation. At the initial stage are observed only h Al , radiogenic radicals of quartz, and resin radicals, these latter being very diluted. The manufacturing process does not annihilate h Al , but it generates new radicals (R and X), not observed in previous studies on CS bearing materials. These radicals are specifically linked to the quartz presence/absence and not to the type of processing (wet/dry). ESEEM shows that only weak dipolar interactions with H nuclei are operative, suggesting that R and X radicals are located neither in the inorganic portion, nor in the resin, but rather at their interface.
Our data suggest that the operating temperature during the process is confined in the range between 260 °C and 380 °C (values at which the Ti and h Al centres are annihilated, respectively) 35 . This range compares well with simulations and experimental determinations reported by several authors [38][39][40][41] , and is lower than the degradation temperature of most polyester resins 42,43 . Unfortunately, we have not enough spectroscopic information to extend this consideration to the wet process too. Nevertheless, similar, or lower, T values can be expected. This information allows establishing that the resins effectively persist 19 , but they do not react chemically to form additional radical species. ( Table 4. Main parameters for EXAFS analysis. Notes: S 0 2 = Amplitude reduction factor, R = refined path distance, N = path degeneracy, σ 2 = Debye-Waller factor. * = Fixed parameters. Metallic Fe (Fe 0 ) has been fitted according to the BCC structure using the correlated Debye model in order to compute Debye-Waller factors for each path 45  www.nature.com/scientificreports www.nature.com/scientificreports/ Thus, we point out that the main role played by resins is the protection from annihilation of the inorganic radicals here described. The fact that radical species appear protected, in full agreement with the findings of Pavan et al. 19 , confirms the anomalous behaviour of the interaction between silica dusts arising from AS, if compared with other conventional sources of exposure to RCS.

Conclusion
The information provided by the present multianalytical investigation of artificial stone samples allows us to complete the results provided by Pavan et al. 19 . We confirm the extremely wide variability of the raw AS materials. The processed dusts clearly differ in the chemical signature from the parent AS, and they also differ in response to the operated dry/wet process. Moreover, significant differences in the Fe and radical speciation are observed between raw and processed materials. From these considerations, workers operating at different tasks may be exposed to dusts having different features, and the toxicity of the AS may be linked not only to internal properties but also to modifications externally induced by the processing. A specific task linked to the variable exposure for ASW is represented by the large amount of cristobalite, found in samples #4: this phase can play a confusing role in the evaluation of the degree of toxicity of the ASW cohorts 44 . The stability of the radicals observed in this study, and the persistence of the resin after the processing could suggest that resin, coating the RCS particles, is able to protect surface radicals for a limited time during the interaction with the lung lining fluids. This mechanism can thus be correlated to the higher toxicity of the AS powder with respect to the conventional silica bearing materials, whose unprotected radicals could be annihilated before reaching the lung tissues. Nevertheless, the fact that at least a part of the contaminants is available for the interaction with tissues immediately after the exposure (i.e. the metal contamination provided by the type of processing), while suggesting that complex paths of toxicity may arise, seems in line with the observation provided by Hoy et al. 1 , concerning the role of high level exposures during short periods in a working day.