Organ Metallome Processed with Chemometric Methods Enable the Determination of Elements that May Serve as Markers of Exposure to Iron Oxide Nanoparticles in Male Rats

The systemic influence of iron oxide nanoparticles on the elemental homeostasis of key organs was examined in male rats. In tissues taken at different intervals from nanoparticles injection, the dynamics of elemental changes was analyzed. The organ metallome was studied using total reflection X-ray fluorescence. The obtained data were processed with advanced cluster and discriminant analyses—to classify the tissues according to their organs of origin and to distinguish accurately the nanoparticle-treated and normal rats. Additionally, in the case of liver and heart, it was possible to determine the elements of highest significance for different treatments, which may serve as markers of exposure to iron oxide nanoparticles. Electronic supplementary material The online version of this article (10.1007/s12011-020-02104-z) contains supplementary material, which is available to authorized users.


Introduction
Nanoparticles (NPs) are defined as particles with two or three dimensions in the range from 1 to 100 nm. Their properties depend on the shape, chemical composition, and size which are comparable to biological molecules and subcellular structures [1,2]. Among different nanomaterials, iron oxide nanoparticles (IONPs) have the greatest potential to serve for biomedical purposes. Because of their magnetic properties, IONPs can be used as contrast agents in magnetic resonance imaging, and their ability to interact with the external magnetic field gradients makes them potentially ideal drug carriers [2][3][4][5][6]. Furthermore, in response to an external magnetic field, they can also induce local hyperthermia in pathologically changed tissues [2,4]. Thus, the use of IONPs may improve effectiveness of medical diagnostics and therapy and/or reduce their negative side effects [2,4,6].
For successful development of IONPs medical application, it is necessary to monitor their potential toxicity considering relationships between the shape, size, or chemical properties of the surface of the nanomaterials including their in vivo behavior [1,7,8]. It is important to verify the final fate of NP-based drug carriers after the drug release at the site of action and to determine whether and how the nanomaterials themselves are metabolized, degraded, and/or effectively removed from the body [1,3].
The essential stage of the in vivo investigations is animal tests, which allow to follow the systemic impact of NPs on a living organism [9,10]. The evidence gained through such a study points out that IONPs may accumulate within organs such as liver, brain, heart, or kidneys, leading to transient histopathological and functional changes [11][12][13][14][15][16][17]. One of the possible ways to obtain information about the biodistribution and biokinetics of IONPs as well as the range of their effects is exploration of elemental changes occurring in organs of NPs-treated animals [9,[18][19][20][21][22]. In combination with multivariate statistical analyses, this approach may allow to determine elemental markers of biological effects following exposure to nanoparticles.
Multivariate methods improve the analysis of multidimensional data, allowing useful comparisons to be made from the elemental composition of the samples [23,24]. The most important ones for consideration here include discriminant analysis and cluster analysis, the main goals of which are simplification of the data structure, grouping observations, and investigation of the relationship between variables [25][26][27].
In this work, the utility of multivariate methods for the determination of the elemental markers of exposure to low doses of polyethylene glycol-coated IONPs (PEG-IONPs) was investigated. For this purpose the main organs taken from three groups of rats treated with PEG-IONPs and normal animals were examined. Three different time intervals counting from NPs administration were studied, namely 2 h, 24 h, and 7 days from the injection. Cluster analysis was used for unsupervised classification of organs and examined groups of animals. Similar two approaches were used in the case of discriminant analysis. However, the latter classification was conducted into the populations defined a priori, before analysis. Additional purpose of the discriminant analysis was determination of the original variables (elements) playing the greatest role in the process of groups discrimination.

Experimental Animals and Sample Preparation
The tissues subjected to elemental analysis were obtained from male Wistar rats, which came from the colony of the Department of Neuroanatomy, Institute of Zoology and Biomedical Research, Jagiellonian University in Krakow. All procedures related to the use of animals were carried out according to the agreement no. 121/2015 of the First Local Ethical Committee on Animal Testing at the Jagiellonian University in Krakow and were performed in accordance with the international standards. The animals were grown at the constant temperature (20 ± 2°C), with unrestricted access to food (Labofeed) and maintaining 12 h/12 h light/dark cycles. On the 60th day of their postnatal development, 24 rats were divided into four equinumerous groups and three of them were injected through the tail vein with 1 ml of PEG-IONPs solution (described below), while the controls (group N, n = 6) were treated with the same volume of saline. PEG-coated magnetic iron (II, III) oxide nanoparticles (Sigma-Aldrich 747,408) were dispersed in 15% mannitol (Baxter), and the final concentration of Fe in the injected solution was 8.14 ppm. The hydrodynamic diameter and zeta potential of NPs (both values determined in mannitol solution after removing of the aggregates/agglomerates) were equal to 35 nm and − 98 mV, respectively.
Three groups of NPs-treated animals were anesthetized and perfused with saline 2 h (2H group, n = 5 1 ), 24 h (24H group, n = 6) and 7 days (7D group, n = 6) from the injection. The selected organs (brain, heart, kidneys, liver, spleen, and muscle) were taken from the animal body immediately after the perfusion. The elemental analysis was done on the liquid organ samples which were obtained in the process of digestion in high purity 65% nitric acid. This was carried out in the Speedwave 4 microwave digestion system (BERGHOF). As an internal standard, gallium solution in concentration of 10 ppm was used at 0.3 ml per 1 ml of digested tissue sample. A detailed protocol of the whole procedure can be found elsewhere [20,22].
In the in vivo studies of IONPs toxicity that have been carried out so far, animals were exposed to relatively high doses of Fe. They equaled from 0.8 to even 12 mg Fe per kg of body mass and were equal or much higher than doses used clinically in humans (~mgs of Fe per kg of body weight) [28][29][30][31]. In our experiment, each animal obtained 0.00814 mg of Fe which may seem little. However, if we do simple calculations taking into account (1) the typical Fe contents in wet mass of rat organs [32], (2) the masses of examined organs in adult rats, and (3) the mass of administered iron, we can easily show that the injected Fe constitutes from around 1% (for the liver) to even 35% (for the kidney) of total Fe mass in the organ and, the same, might affect the elemental homeostasis of the rat organism. In the context, we need also to take into account the route of administration, as in our study, NPs were injected directly to blood and not administered orally or intraperitoneally.

Elemental Analysis
The elemental analysis of tissue samples was performed using the total reflection X-ray fluorescence (TXRF) method in the Laboratory of X-ray Methods of the Institute of Physics at the Jan Kochanowski University in Kielce. S2 PICOFOX (Bruker Nano) automatic analytical instrument was used for this purpose. The details concerning the used apparatus and measurement conditions were shown in our previous papers [20,22].
The analysis of TXRF spectra allowed a quantitative analysis of the following elements in the measured tissues: P, S, K, Ca, Fe, Cu, Zn, and Se. Their concentrations in the liver, heart, brain, spleen, kidneys, and muscles in each of the examined animal were the subject of further statistical study including those using multivariate methods.
The calculations necessary to obtain information about the concentrations of elements present in the samples were done in two steps. First, the concentration C i of the element i in the measured liquid sample was calculated according to Formula (1). In turn, the final elemental concentration C x i within organ was calculated, taking into account dilution factor d and the liquid to organ mass conversion factor f x what was described by Eq. (2). where: C IS − concentration of the internal standard (Ga) in the liquid sample (ppm), N IS − net pulse number for the internal standard in the sample spectrum (a.u.), N i − net pulse number for the element i in the sample spectrum (a.u.), s i − relative sensitivity for the element i. where: liquid to organ mass conversion factor for organ x (a.u.).

Multivariate Statistical Analysis
Two multivariate methods, namely cluster and discriminant analysis, were used for advanced statistical evaluation of the elemental data. Cluster analysis is typically used to search for patterns in a data set by grouping the observations into clusters in such a way that objects within each cluster are similar while the clusters are dissimilar to each other. Two approaches are used for this purpose-hierarchical and divisive methods. In the case of hierarchical clustering, in the beginning, each observation is treated as a single cluster. Further combining of observations into clusters (agglomeration) occurs with the increase of the degree of non-similarity of objects until one cluster containing all observations is received. In divisive methods, this process is reversed, but in both cases, the result is a structure called a dendrogram [25].
There are several hierarchical methods of clustering (nearest neighbor, farthest neighbor, weighted pair-group average method, etc.); however, the most popular is Ward's method using an analysis of variance approach to create the links. The sums of squares of deviations within the clusters are minimized, and at each level of combining of clusters, a pair creating the cluster with the least possible diversity is chosen [25]. Ward's method was applied also in this work.
Grouping of objects into clusters is possible, thanks to the dissimilarity function being the measure of the distance between the observations. The most commonly used distance measures are the Euclidean or squared Euclidean distances. The Euclidean distance between the two vectors x = (x 1 , x 2 , …, x p ) and y = (y 1 , y 2 , …, y p ) is defined as [25]: The main task of the discriminant analysis is to confirm the existence of differences between two or more groups to which the objects were assigned. To elucidate these differences, linear functions of the variables (discriminant functions) are used. Discriminant functions are determined in such a way that the external variability (between populations) is significantly greater than the internal variability (between the objects within the same populations). In order to check which discriminant functions are necessary to differentiate groups, Wilks' Λ statistics, which can be approximated with a chi-square distribution, is used. In turn, F Fisher' statistics corresponds to partial Wilks' Λ, which is utilized to determine the variables playing the greatest role in the process of discrimination [25].
The choice of methods of multivariate analysis in this paper was motivated by the nature of the samples tested. In our study, each object was described by a set of features which constituted the concentrations of particular elements.
It is necessary to mention that the data set used within this study has been partially utilized before in the papers of Matusiak et al. [22] and Skoczen et al. [20]. In the study of Matusiak et al., the changes in the contents of 4 elements, namely Ca, Fe, Cu, and Zn, were examined only in the liver of animals exposed to NPs action and controls, while in the work of Skoczen et al., the anomalies of the same 4 elements were measured for the remaining organs (kidneys, spleen, heart, and brain). In this paper, for the first time, we examined the whole data set obtained during the animal experiment, which now include also information about changes occurring in the accumulation of P, S, K, and Se. Besides the analysis of additional elements and tissue (muscle), the obtained results were based on a different research approach in which all the data were treated as one set of variables. The used solution was justified by the European Union directive on the use of laboratory animals, which obliges researchers to comply with the 3Rs principle. One of its key elements is the reduction of the number of animals while maintaining the statistical significance of subsequent analyses. In our work, this principle was respected, trying to use all information from the collected till now material. Repetition of the experiment on new animals in order to obtain a relatively small, complementary material that can be obtained from the same samples would be inconsistent with the recommended research standards.

Results
The result of the quantitative elemental analysis were the concentrations of elements (P, S, K, Ca, Fe, Cu, Zn, and Se) in the main internal organs (liver, heart, brain, spleen, kidneys, and muscles marked as Li, He, Br, Sp, Ki, Mu, respectively). As the first step, the median concentrations of elements were calculated and the differences between PEG-IONPs-treated animals and controls were tested using non-a parametric U Mann-Whitney test at the 5% level of statistical significance. The obtained median values together with their uncertainties calculated as interquartile spans are presented in Table 1. The statistically significant differences between rats injected with NPs and normal animals were also marked there.
The only organs sharing statistically significant differences in Fe levels were the kidneys ( Table 1). As a result of PEG-IONPs injection, the level of Fe increased in kidneys after 2 h and 7 days from the injection. Nonetheless, the treatment of animals with NPs, led to many other significant elemental abnormalities of the examined organs, which were described and discussed in details in our previous papers [20,22].
In this paper, the data on accumulation of elements in organs have been treated as a set of features describing each animal (object) and, for the purpose of their simultaneous comparisons, two multivariate methods (cluster and discriminant analysis) were applied. They were used in two approaches, for classification based on the examined organs and classification due to the analyzed groups of animals. The data were explored using the STATISTICA software (version 7.1).

Evaluation of the Statistical Significance of Differences in the Elemental Composition of Main Organs
The first approach allowed to check whether samples of the examined organs can be distinguished based on their elemental composition. The analyses were carried out for each animal group (N, 2H, 24H, 7D) separately.

Cluster Analysis
Cluster analysis was performed on the concentrations of elements after the data standardization using the STATISTICA software. The hierarchical Ward method of agglomeration and the Euclidean distance as a measure of the distance between the observations were applied. The results are presented in the form of dendrograms in Fig. 1. For all the studied animal groups, the samples from particular organs form distinct clusters. However, in the case of spleen, the allocation of all samples into one cluster was observed only in the 24H group, while for the N, 2H and 7D groups, some of the spleen samples showed a greater similarity to the samples originating from the brain. In addition, the shapes of dendrograms for animals from the N and 24H groups as well as the 2H and 7D groups presented pairwise similarities. In the case of the first pair, two distinct clusters could be observed. The first one included muscle, brain, and spleen, while the second one-heart, kidney, and liver samples. For the 2H and 7D groups, three distinct clusters could be noticed instead of two. The first of them included samples from muscle and heart, indicating similarity of their elemental composition. The second cluster consisted of brain and spleen samples, while the third contained kidney and liver tissues. It is worth to mention that the classification of heart samples together with muscle samples is the main differentiator between the second (2H and 7D) and the first pair of groups (N and 24H).

Discriminant Analysis
First, the analysis was done for the N group using a standard method in which all primary variables were taken into account. Five discriminant functions were distinguished and their characteristics are presented in Table 2.
The data presented in Table 2 confirmed the statistical significance of all the obtained discriminant functions. Therefore, their raw canonical coefficients are shown in Table 3.
Based on the data from Table 3, five significant discriminant functions could be calculated and each object could be characterized by them instead of eight original variables. The dependences describing how the first two discriminant functions depend on raw primary variables were presented in Eqs.
(4) and (5). As one can see it, the discriminant functions are sums of the original variables weighted by respective raw coefficients of discriminant function.
By replacing the primary variables with the D 1 and D 2 coordinates of the discriminant functions, a scatterplot of observations was obtained (Fig. 2a).
In order to determine the elements of the highest significance for the model, the partial Wilks lambda, F test of partial Wilks' lambda, and its p value were used. All these parameters are shown in Table 4.
Based on the data presented in Table 4, the elements of the greatest importance in the process of organ distinguishing in the case of the N group were determined.  These were P, S, K, Fe, Cu, and Se; the elements for which the lowest partial Wilks lambdas, largest values of F Fisher statistics, and p values less than 0.05 were observed. Further analysis was carried out for these six elements, and the newly calculated raw canonical coefficients are presented in Table 5 and used to create a scatterplot of observations in the space of discriminant variables shown in Fig. 2b. Despite the reduced number of variables in the model, different organs taken from normal animals are still very well separated in the space of discriminant variables (Fig. 2b).
A typically used measure of differences between the studied groups in the space of discriminant variables is the squared Mahalanobis distance between their centroids-meaning the centers of gravity of the observations in this space. The squared Mahalanobis distances between different organs representing the N group together with the F Fisher statistics and its p values are included in Table 6. All the p values presented there are less than 5%, which confirms the statistical significance of the calculated Mahalanobis distances and significance of the differences between the six analyzed organs for the N group.
An analogous approach was used for the 2H, 24H, and 7D groups. In all cases, it was shown that the examined organs were well separated in the space of discriminant variables. The elements contributing the most to the organ discrimination in the 2H group were S, K, Fe, Cu, Zn, and Se; in the 24H group P, S, Ca, Fe, Cu, and Se; and in the 7D group, P, S, K, Fe, Cu, and Se. Excluding of the redundant elements from the model allowed still good differentiation of examined organs in the space of discriminant variables what one can see from Figs. S1-S3 presented in the Supplementary materials. What is more, these qualitative results have been confirmed by statistics related to the squared Mahalanobis distances.
Discriminant analyses for N, 2H, 24H, and 7D groups are characterized in Table 7 by numbers of both all and statistically significant discriminant functions and the elements of the highest significance for the models determined based on partial Wilks' lambdas.   Accordingly to Table 7, all determined discriminant functions were statistically significant in the case of N, 2H, and 24H groups. Regardless of the studied animal population, the elements such as S, Fe, Cu, and Se played an important role in the process of organ differentiation. However, in order to obtain good separation of cases in the space of discriminant variables, depending on the studied group, one had to take into account P and K for N and 7D groups, K and Zn for 2H group, and P and Ca for 24H group, additionally.

Evaluation of Differences in the Elemental Composition Between Organs from IONPs-Treated and Normal Rats
The second methodical approach allowed to check whether a tissue sample of given main organ could be assigned to a specific group of IONPs-treated animals or to untreated controls based on the analysis of its elemental composition.

Cluster Analysis
Cluster analysis was performed on the standardized concentrations of elements in the main internal organs. For this purpose, the Ward's hierarchical method of agglomeration and the Euclidean distance as a measure of dissimilarity between objects were applied. The results obtained for each organs are presented in the form of dendrograms and compared in Fig. 3.
The cluster analysis carried out for the examined organs indicated that the differences in elemental composition of samples from the studied groups were too small to enable the separation of the population in the Ward agglomeration process. However, some regularity can be observed for the liver, for which the majority of cases representing four studied groups of animals were classified into separate clusters.

Discriminant Analysis
At the beginning, the analysis was done for the liver using a standard method in which all primary variables were taken into account. Three statistically significant discriminant functions were distinguished and their characteristics are presented in Table 8.
As one can see it from Table 8, all separated discriminant functions were statistically significant. Therefore, for all of them, raw canonical coefficients are calculated and presented in the Table 9.
Based on the data included in Table 9, three discriminant functions could characterize the object instead of eight original variables. Equations (6) and (7) define how the first two discriminant functions depend on raw primary variables. These equations were used to obtain the scatterplot of observations in the space of discriminant functions which is presented in Fig. 4a.  Figure 4a shows that livers taken from normal and 2H rats were well separated in the space of discriminant variables  from those taken from two remaining experimental groups. Thus, standard discriminant analysis based on concentrations of elements in livers allowed the differentiation of IONPstreated rats from controls. In Table 10, the parameters used to determine the elements of the highest significance for the model in the case of liver were presented. They indicate Ca, Cu, and Zn as the elements most important for distinguishing the livers taken from IONPs-treated animals. Therefore, further analysis was conducted only for these variables and the new raw canonical coefficients of discriminant functions were presented in Table 11.
The scatterplot of observations in the space of the first and the second discriminant functions is presented in Fig. 4b. Exclusion of the redundant variables (elements) from the model improved the separation of the cases representing the groups 24H and 7D but worsen the differentiation of the cases from the groups N and 2H.
The squared Mahalanobis distances between livers representing four examined groups together with the F Fisher statistics and their p values are included in Table 12. All the p values presented there are less than 5%, which confirms the statistical significance of the calculated Mahalanobis distances.
Analogically, the discriminant analysis was done based on the data concerning the elemental composition of other organs under analysis. The results obtained for kidneys are presented in Tables S1-S3 and Fig. S4 of Supplementary materials, for heart in Tables S4-S7 and Fig. S5, while for brain, spleen, and muscles in Tables S8, S9, and S10, respectively.
Due to the correlation occurring between Ca and Zn concentrations found for kidneys, it was not possible to include all the primary variables in the discriminant analysis. As the contribution of Ca in the process of group differentiation was greater than that for Zn (see the data presented in Table 13), the concentration of this element was the subject of further statistical calculations.
For both the kidneys and heart, three discriminant functions were distinguished but only one of them was statistically significant (Tables S1 and S4). On the basis of raw canonical coefficients of discriminant functions (Tables S2 and S5), two new discriminant variables were created for the two organs and used to present the observations in the space of discriminant functions (Figs. S4 for kidneys and S5 for heart).
For the kidneys (Fig. S4), first discriminant function gave good separation of control animals (N group) from those injected with IONPs. In turn, consideration of a second discriminant function did not contribute to a better differentiation of the examined animal groups. Table S3 shows that none of the primary variables contribute significantly to the discrimination process. Therefore, statistical analysis for the kidneys was terminated at this stage.
In the heart case, discriminant analysis based on all the primary variables allowed satisfying separation of all animal groups in the space of discriminant variables. Two elements, namely K and Se, were determined as the most relevant for the process of the distinguishing the examined populations. However, removal of the redundant variables from the model significantly worsened classification of cases compared with the analysis carried out for all primary variables what can be seen in Fig. S5. The squares of the Mahalanobis distance were, therefore, not calculated.
The analyses performed for the brain, spleen, and muscles did not allow to determine any statistically significant discriminant functions (Tables S8-S10) neither primary variables relevant for the process of group differentiation.

Discussion and Conclusions
As can be seen from the existing literature, multivariate methods, also known as chemometrics, are frequently used in the analysis of elemental data especially in the field of life sciences [24,[33][34][35][36][37][38][39][40][41]. For example, hierarchical cluster analysis was successfully applied for the prediction of cancerrelated Cu-binding proteins, which can serve as a source for mechanistic-molecular studies of Cu-dependent processes in cancer [42]. The method was also used for quantitative imaging of the influence of chronic Mn exposure on metal accumulation in hippocampal formation [43]. In turn, discriminant analysis was adapted inter alia, for characterization of metal profiles within serum in Alzheimer's disease [34] and gouty arthritis [44] as well as identification of metal biomarkers of different cancers in human tissues [39,45,46]. In the literature, there are also many examples of the chemometric method application in the research utilizing X-ray fluorescencebased techniques [24,36,38,47,48]. Both multivariate methods, which were used in the study, proved that the elemental composition of examined organs is so unique that enables to correctly classify them into the groups that were separated during analysis or defined a priori. The shape of the dendrograms obtained for examined groups showed big similarity between animals from which organs were collected 2 h and 7 days following IONPs injection. This is in agreement with the data presented in Table 1 and also with our previous results showing that these time intervals are connected with the most intensive elemental Fig. 3 Dendrograms presenting the results of unsupervised classification of animals based on the elemental composition of the liver (a), kidneys (b), brain (c), heart (d), muscles (e), and spleen (f) anomalies occurring in the internal organs [20,22]. Moreover, the discriminant analysis showed that, regardless of the animal population studied, S, Fe, Cu, and Se are the elements that play important roles in the process of organ categorization.
Unsupervised classification (cluster analysis) of cases according to the animal groups showed some regularities only for samples taken from the liver. For this organ, the rats representing different populations were usually categorized to separate clusters. Such patterns were not observed for other organs including those for which significant elemental abnormalities occurred as a result of exposure to NPs. Better results were obtained using discriminant analysis. The scatterplots of observations were obtained for the liver, kidneys, and heart. For other examined tissues, nonstatistically significant discriminant functions were determined. For the three mentioned above organs, our previous investigations showed significant elemental anomalies in NPs-treated rats [20,22] compared with controls, and such a result is probably connected with the functions performed by these organs. The discriminant analysis carried out for the liver, kidneys, and heart using all primary variables allowed good distinguishing of the control group from IONPs-treated animals. Therefore, information on elemental composition of one of the three mentioned organs supported by this multivariate method could be used for verification whether the animals were or were not treated with NPs.
Based on the partial Wilks lambdas, the elements of highest significance in the process of group discrimination were found. In the case of liver, these were Cu, Ca, and Zn. As one can notice from Table 1, the level of Cu decreased while those of Ca and Zn increased in livers of PEG-IONPs-treated rats. Regulation of the iron metabolism is one of the main liver tasks [49]. The liver is responsible for the iron storage in a form of ferritin, regulation of this element movement through the production of peptide hepcidin, and the synthesis of major iron metabolism proteins, such as transferrin and ceruloplasmin [50]. The observed temporary (in 2H group) and longterm anomalies of liver Cu level could be an effect of ceruloplasmin release to blood in response to the increased iron level and are in agreement with our earlier results demonstrating increased concentration of Cu in rat serum after PEG-IONPs administration [22].
Also, Zn and Ca concentrations were changing in the liver as a consequence of the treatment with PEG-IONPs. The levels of both elements increased between 1st and 7th day from the injection what could indicate the occurrence of inflammatory processes being a consequence of Fe-induced oxidative stress [22]. Zn together with Cu is the part of Cu,Zn-superoxide dismutase (Cu,Zn-SOD)-an oxidoreductase enzyme responsible for rapid two-step dismutation of the toxic superoxide radical to molecular oxygen and hydrogen peroxide [51]. Because of its function, Cu,Zn-SOD is a crucial component of the cellular response to oxidative stress [51]. In our study, simultaneous increase of Cu and Zn concentrations within the liver was not observed. Therefore, an elevation in the level of Cu,Zn-SOD enzyme in the organ occurring in response to iron-induced oxidative stress cannot be confirmed. On the other hand, this cannot be ruled out either, because the increase in copper level associated with Cu,Zn-SOD could be masked by the release of copper, probably bound to ceruloplasmin, from the liver which was observed within our previous study and was described in details elsewhere [20,22].
Significant increase of Ca content, observed for the liver between 1st and 7th day after PEG-IONPs administration could point at the role of Ca in signaling of increased production of reactive oxygen species [52]. In turn, the disturbances in the calcium economy resulting from persistent oxidative stress could lead to the loss of efficiency of adaptive mechanisms and, consequently, to the appearance of apoptosis [53].  Kidneys were another organ for which discriminant analysis allowed to distinguish the control group from IONPstreated rats. However, in contrary to liver, the determination of elemental markers of NPs exposure was not possible for them. Kidneys are responsible for removal of small unnecessary or toxic particles from the organism [54]. However, this ability is strictly connected with the size of these elements. It was shown that nanoparticles with the diameter higher than 10 nm cannot be removed from the organism by the kidneys [21,55]. Our previous results are in agreement with these data showing the redistribution process of PEG-IONPs [20].
In the case of heart, the discriminant analysis allowed both to differentiate normal rats from those injected with NPs and to determine the elements significant for the discrimination process, namely K and Se. Despite the lack of statistically significant changes in the iron level occurring in hearts of animals subjected to PEG-IONPs, the intense flow of blood with NPs through the organ may trigger long-term anomalies manifesting, among others, in elemental changes of the organ [20,56]. In the recent literature, there is no unambiguous information about the nature of the observed changes. Some authors indicate a negative effect on heart related to exposure to IONPs [57], while others show their protective nature [58,59]. Therefore, further research in this area should be carried out.
Potassium is an element of great importance for proper functioning of the nervous and muscular systems with particular emphasis on myocardial work. The right level of potassium lowers blood pressure, has positive effect on the endothelium of the vessels, and also reduces the production of free radicals what decrease the risk of brain stroke [59,60]. What is more, the extracellular potassium concentration influences the electrical stability of the heart [60]. Increases of extracellular potassium ions concentration in serum are responsible for many early electrophysiological changes, and as a result can result in cardiac arrhythmias in both physiological and ischemic myocardium [60]. Therefore, determination of K as the element important for differentiation of normal rats from those treated with NPs could indicate some pathological processes occurring in heart. Such a conclusion is in agreement with our previous observation, showing increased Ca levels in the organ of IONPs-treated animals [20]. Ca plays an important role in the regulation of secretion and release of hormones into the bloodstream. Moreover, it influences the transmission of nerve impulses and muscle contractility [61]. Also, proper functioning of the vascular system and myocardium as well   as coagulation process are associated with the occurrence of calcium ions [62,63]. The second element proposed as the biomarker of PEG-IONPs action on heart was Se. The role of selenium in living organism is mainly connected with regulation of the inflammatory response, antioxidant properties, the proliferation/ differentiation of immune cells, and optimal functioning of the cardiovascular system [64]. Selenium in a form of selenoproteins regulates various signaling processes by influencing the redox homeostasis and cellular Ca 2+ influx [65,66]. The determination of Se as an element important for hearts discrimination together with the fact that the level of Ca was disturbed in the organ seem to confirm the existence of oxidative stress in rats treated with PEG-IONPs.
Reassuming the obtained results, it could be concluded that the applied chemometric techniques enabled correct classification of tissues according to the organ of origin and, what is more significant, discrimination of rats subjected to NPs action from normal animals. Additionally, in the case of liver and heart, it was possible to determine the elements of the highest significance for the differentiation process which may be the candidates for the markers of exposure to nanomaterials which could potentially be applied in diagnostic procedures or medical treatments.
This study should be continued. Especially, additional technique providing the direct evidence for the deposition of PEG-IONPs (e.g., TEM) in the organs should be used. Neither TEM analysis nor histological staining was used in this study, as both techniques require appropriate and different protocols of biological material preparation. While, the reliable elemental studies are possible only using unfixed tissues. The tissue which, after a short perfusion with saline, was not fixed but only frozen, cannot be the subject of immunolabeling. Also for TEM, a special fixative (glutaraldehyde with paraformaldehyde) must be used. However, as we have shown it in our previous paper, the fixation may strongly affect the elemental composition of the tissue [67]. Therefore, we decided to carry out the elemental analysis on unfixed samples what, unfortunately, strongly limited the possibilities of their further examinations.
On the other hand, both histological and TEM study are usually limited to very low amount of the biological material. This is typically a few-micrometer-thick slice of the organ per animal. Taking into account such small amount of the sample at the beginning stage of the study, we could easily miss the areas of NPs accumulation or action. Therefore, at the first stage of the research, the general elemental anomalies occurring within organs were examined.
Our future research should also be focused on the evaluation of the blood levels of the examined elements, as such data could make easier the interpretation of the obtained statistical analysis results. Such data were not presented in this paper because of the poor optimization of the procedure of serum preparation for the analysis of Fe and other elements. The methods we used till now did not allow us to avoid the hemolysis of erythrocytes during blood centrifugation and as consequence, to obtain the reliable results of elemental analysis of serum samples.
Acknowledgments Special thanks for Dr. Aldona Kubala-Kukus and MSc Ilona Stabrawa from the Jan Kochanowski University in Kielce for their great help in performing the TXRF measurements.
Funding Information This work was partially financed by the Faculty of Physics and Applied Computer Science AGH UST statutory tasks No.11.11.220.01/3 within subsidy of Ministry of Science and Higher Education.

Compliance with Ethical Standards
All procedures related to the use of animals were carried out according to the agreement no. 121/2015 of the First Local Ethical Committee on Animal Testing at the Jagiellonian University in Krakow and were performed in accordance with the international standards.

Conflict of Interest
The authors declare that they have no conflict of interest. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.