The sides A and B of water sphere hit by 1 MeV electrons from the opposite directions got the maximum dose to compare with the sides C and D (Fig. 3a). In view of such a dose distribution, a 3–4 mm thick layer of water sphere from sides A and B were irradiated with 1 MeV electrons. According to that, а 1-mm thick skin was removed from the tubers, and a 1-mm thick layer of the tuber was peeled for the subsequent extraction.
In water parallelepiped, the layers located close to beryllium window of the X-ray tube received the maximum dose (Fig. 3b). Further, the dose plummeted as photons were penetrating deeper than 1 mm into the phantom.
Development of a chemical test. The task was to differentiate between potato tubers irradiated with different doses and the non-irradiated samples. Our attempt to use the classical fingerprinting approach for this purpose was not successful: the samples did exhibit intrinsic fluorescence, but its intensity did not change as a result of irradiation. The technique of adding fluorophores to potato samples was not successful either: the signal of the added dye did not change with the dose.
We had to turn to a more sophisticated fluorescent fingerprinting technique which we suggested recently. This strategy relies on the effect of analytes on two types of indicator reactions: (1) catalytic oxidation of a dye with H2O239 and (2) the formation of fluorescent aggregates of the dye with the analyte and a surfactant33,40. By the analytes we imply reactive chemical components of the sample whose concentration could have changed as a result of radiation treatment. In system (1), the analytes are supposed to bind with the transition metal ion (Cu2+ or Pd2+) that catalyzes the oxidation of a carbocyanine dye with hydrogen peroxide. As a consequence, the oxidation reaction rate changes, which is tracked by the visible light absorption and by the fluorescence intensity change of the carbocyanine dye in solution39.
In system (2), the reactive analytes form a “hydrophobic ion pair”40 with an oppositely charged surfactant. For the biological species supposedly containing in the sample, which are for the most part anionic, a cationic surfactant cetyltrimethylammonium bromide (CTAB) was used. The hydrophobic ion pairs were found33 to incorporate hydrophobic dyes, which is accompanied by fluorescence enhancement. The types (1) and (2) indicator reactions were studied by us earlier, and the concentration conditions for the present research were selected based on the previous studies35,39. In reactions of type (2) we employed hydrophobic dyes prone to forming fluorescent aggregates with a surfactant and an oppositely charged analyte (dyes 1 and 4)33,35,39. In type (1) reactions we used less hydrophobic and more water-soluble carbocyanines that could be easily oxidized by H2O2 (dyes 2–4)35,39. The concentrations of the oxidant and metal ion-catalyst were chosen to provide a reaction rate convenient for observation (1 M hydrogen peroxide, 1 mM Cu2+). The general strategy is schemed in Fig. 4.
The development of signal over time was recorded photographically in the visible and near-IR region (Fig. 5, each three adjacent columns represent a sample). It can be seen that the differences between samples are not very pronounced. There was a question, whether these differences were sufficient to distinguish between the absorbed doses. To solve this problem, the images were digitized and subjected to chemometric treatment. The indicator reactions and concrete images used in chemometric treatment are listed in Supplementary Table S2. The results were presented as scores plots allowing to more explicitly visualize of the differences between the samples.
Discrimination of tubers irradiated by an electron beam. For the set of experiments No 1, the different doses were not discriminated by using principal component analysis (PCA). We turned to discriminant analysis (LDA), which is a supervised technique, i.e. the software is informed about which signals belong to which class (irradiation dose). Before performing the LDA procedure, all data were divided into training set and validation set. The training set is used by the LDA software to establish the relationship between the sample and its dose (create the model), and the validation set is used to verify the quality of this model. As a compromise, we chose a validation set size as 28–30% of all data (14–15 samples out of the total amount of 50). A larger number of validation samples would shorten the training set and impede the discrimination quality.
LDA yields a 4-dimensional score plot in the coordinates of factors F1–F4. The potato samples irradiated with different doses are represented as groups of points in this plot. The accuracy of discrimination by LDA technique is estimated as the percentage of correctly assigned validation samples. For example, the points in Fig. 6 correspond to 36 training and 14 validation samples. Since the discrimination is not perfect, it is difficult to assign some of the validation points to specific groups of training samples. The software considered such assignment as correct if the Mahalanobis distance between the validation point and the training sample group of the correct dose was minimal among the distances to the other groups. In other words, the correctly assigned validation point should be located closer to the group of points of the dose it originated from. The distances for such disputed points are illustrated in Fig. 6. Accuracy was calculated as the ratio of the number of correctly assigned validation points to the total number of validation points; the software presented this result as confusion tables (Table 1).
Overall, for the data demonstrated in Fig. 6, the tubers irradiated with an unknown dose can be assigned to one of the doses (10; 100; 1,000 and 10,000 Gy) with 100% accuracy using the basic data set.
Table 1
An example of confusion matrix for the validation points obtained by the XLSTAT software. Incorrectly predicted doses are shown in bold font
True
dose, Gy
|
Predicted dose, Gy
|
Total No
of points
|
Accuracy*
|
0
|
10
|
100
|
1,000
|
10,000
|
0
|
1
|
0
|
1
|
0
|
0
|
2
|
50,0%
|
10
|
0
|
3
|
0
|
0
|
0
|
3
|
100,0%
|
100
|
0
|
0
|
3
|
0
|
0
|
3
|
100,0%
|
1,000
|
0
|
0
|
0
|
2
|
1
|
3
|
66,7%
|
10,000
|
0
|
0
|
0
|
0
|
3
|
3
|
100,0%
|
Total No of points
|
1
|
3
|
4
|
2
|
4
|
14
|
85,7%
|
* Number of accurately predicted values / total number of validation points. |
The red arcs in graph (a) connect the validation points and the barycentres of the nearby groups corresponding to certain absorbed doses shown as black squares; the Mahalanobis distances in the space of factors F1–F4 are shown in red numbers beside the arcs. Validation procedure can be considered, for example, for a validation point located between the doses of 10 kGy and 100 Gy; this point is seemingly closer to the 10 kGy group; however, in the other dimension (F1–F3) that point is at a greater distance from the 10 kGy group. Quantitatively, the 4D distances from this point to 10 kGy and 100 Gy groups are 8.9 and 6.1 units, respectively (shown by the red arcs), which implies that the point belongs to the 100 Gy group. Similar situations can be viewed for the other validation points.
Table 2
Accuracy of discrimination of potato tubers irradiated with electron beam for various data sets of experiment No 1
Data set used for treatment
|
Number of data columns*
|
Accuracy, %**
|
Basic data set
|
22
|
100
|
Aggregation reaction of dye 1 removed
|
20
|
78
|
Oxidation reaction of dye 3 removed
|
16
|
78
|
Oxidation reaction of dye 2 removed (samples with sulfite)
|
16
|
86
|
Oxidation reaction of dye 2 removed (samples with ascorbate)
|
14
|
86
|
Only dye 2 reactions are used
|
12
|
86
|
Only NIR images are used
|
12
|
80
|
Columns with the highest standard deviations are used***
|
7
|
64
|
Data are selected by the largest visual difference between photographs
|
5
|
57
|
* Each column contains data for one indicator reaction at a certain reaction time (concrete characteristics of the data cloumns are given in Supplementary Table S2). |
** Percentage of validation points correctly assigned to their groups (automatically calculated based on the comparison of Mahalanobis distances to the groups. The results were presented as tables similar to Table 1). |
*** Standard deviations were calculated for the columns as measures of data diversity between samples. |
Reduction of data. Discrimination of potatoes irradiated with different doses was based on 3–4 indicator reactions. In order to reduce the number of experiments, an attempt was made to find out whether it was possible to reduce the amount of data. If an indicator reaction (based on dyes 1, 2 or 3) was excluded from the basic data set of 22 columns to give a reduced set of 14–20 columns, the discrimination accuracy dropped to 78–85% (Table 2). The same accuracy was observed if only dye 2 reaction was used or only near IR images were treated (12 columns). Further reduction of data resulted in substantial deterioration of accuracy. The corresponding score plots are given in Supplementary Fig. S1, a–f.
Effect of potato variety. Experiment No 1 was conducted with the potatoes of one variety. However, a situation may arise where the tubers irradiated with known doses belong to one variety, while the tubers irradiated with an unknown dose are of a different variety. For this reason, in experiment No 2 we studied two different potato varieties (X and Y), three tubers of each. The indicator reactions used in experiment No 2 were the same as those in experiment No 1. The dose of 10 Gy was not used in experiment No 2 in order to reduce the number of rows in the data table.
When only the samples of variety X were considered, 100% accurate discrimination of the samples according to the doses was achieved (Fig. 7,a). For variety Y, the 100 Gy and 1,000 Gy samples were poorly separated (Fig. 7,b). Due to that, all doses were correctly assigned only in 93% cases (the average of 5 validation runs, each run was performed with 48 training and 12 random validation samples). For both varieties it was possible to confidently distinguish between the irradiated and non-irradiated samples.
When all data were treated without the distinction of X and Y varieties, the groups of points were also partly overlapping (Fig. 7,с). The accuracy of determination of the dose for the combined varieties was 89% (the result of validation procedure repeated 5 times with 25 validation and 95 training samples). Consequently, the discrimination of doses can be achieved even for the tubers of different varieties, though it can be less efficient that within one variety.
Discrimination of tubers irradiated by X-rays. The same dose of electron beam radiation can be more efficient than X-rays since it can destroy microorganisms more efficiently, as it was found in the study of irradiation of turkey meat13. Nevertheless, we explored the feasibility of discriminating between the X-ray doses absorbed by the potato samples. Two tubers of one variety were studied. Technically, the whole tubers could not be treated, for which reason the potato pieces sized 6 mm × 6 mm × 15 mm were placed in the Eppendorf tubes and irradiated during 50 s to 21 min, which corresponded to doses of 0, 100, 1000, and 5000 Gy. After irradiation, the samples were extracted with water similarly to the samples irradiated with an electron beam, and after 24 h of extraction, the same indicator reactions were carried out in the 96-well plate. Each sample was measured in 5 parallels, which gave totally 40 samples (2 tubers × 4 doses × 5 parallels). The data were collected 3 times during 20 min for the redox reaction, which totally yielded 24 data columns. The whole data set was used for the discriminant analysis treatment.
The results shown in Fig. 7,d demonstrate that the X-ray doses can be completely discriminated. Validation was performed 5 times with a random validation set (each time with 32 training and 8 validation samples). Out of the total number of 40 validation samples, 38 were assigned correctly, which corresponds to the discrimination accuracy of 95%.
Advantages, limitations and prospects. A limitation of the suggested chemical fingerprinting method is the necessity to analyze the reference samples of known composition along with the unknown ones. In other words, the samples irradiated with known doses and control samples should be treated simultaneously with the unknowns. Theoretically, the measurements of known samples could be performed in advance, but the real day-to-day reproducibility of the optical signal in the indicator reactions is not perfect, and the results obtained on one day cannot be reliably used on another day. This is a consequence of the nature of the reactions used in this method and a price to pay for their sensitivity.
In this study, we extracted the samples on the day of irradiation. Some additional research effort is needed to find out how long after exposure it is still possible to determine the dose (for example, if the sample was irradiated a week or two ago).
It is clear that the properties of potato samples may strongly depend on the variety, conditions of growth and other factors, which can eventually alter the discrimination efficiency. This work only shows the fundamental possibility of determining the order of the absorbed dose. Other foods should be also studied in this regard.
The fingerprinting technique does not allow us to reveal the nature of compounds formed in the irradiated samples that influence the indicator reaction rate and allow for the discrimination of samples. However, this knowledge is unnecessary for solving practical tasks.
The advantage of this method over chromatography techniques is its relative simplicity with no spectral or other sophisticated instrumentation involved except for the photo cameras and red light LED source. In this study we used synthesized carbocyanines but the protocol can be later adapted to commercially available dyes. The accuracy of discrimination can be further improved by implementing new indicator reactions that would be more sensitive to the composition of irradiated samples.