Lazar G - Effects of The Environmental Stress On Two Fish Populations Revealed by Statistical and Spectral Analysis
Lazar G - Effects of The Environmental Stress On Two Fish Populations Revealed by Statistical and Spectral Analysis
Lazar G - Effects of The Environmental Stress On Two Fish Populations Revealed by Statistical and Spectral Analysis
net/publication/259641327
CITATIONS READS
9 318
9 authors, including:
Some of the authors of this publication are also working on these related projects:
Health and environmental risks related to the groundwater and soil contamination with organochlorine pesticides and heavy metals View project
All content following this page was uploaded by Gabriel Lazar on 05 April 2014.
Abstract
The aquatic species are permanently subjected to various stresses through the ecosystem and the food habit. Elucidating the
impact at the cellular level allows us to understand the ways in which the polluting elements and their derivatives can disrupt the
metabolism of aquatic organisms triggering apoptosis in individual cells or necrosis affecting groups of neighbouring cells. The
aim of this preliminary study is to highlight the influence of the environmental factors on two fish population by statistical and
spectroscopic methods, allowing the observation of the changes in structural components of cell’s membranes. In many studies, an
accurate method for biomolecules quantification is necessary for the purpose of selecting optimum species and environmental
living conditions. Biological samples were collected from two different ecosystems: the Oituz River placed in a region less
affected by urban and industrial pollution, and the Mures River, located in an area with a higher degree of pollution. Different
samples from brain and muscle and theirs etheric extract were analyzed using ATR-FTIR and UV-VIS spectrometry. Detailed
analysis of the different spectral regions after spectral decomposition reveals differences between the two fish populations. The
infrared spectroscopy was assisted by statistical analyses in order to develop a calibration model and to reveal the mainly spectral
region that correlate with the etheric extract.
Received: September, 2011; Revised final: January 2012; Accepted: January, 2012
Author to whom all correspondence should be addressed: e-mail: [email protected]; Phone: +40 234542411; Fax: +40 234545753
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
protein when the pollutants penetrate the cell through each organ system at the molecular level that is
different transport mechanisms, like endocytosis affected.
(Moisescu et al., 2008). When the production rate of FTIR spectroscopy was applied to monitor the
reactive oxygen and nitrogen species exceeds the rate structural changes in lipids and in amino acids of
of decomposition induced by antioxidants, the proteins (Iwaki et al., 2003), to study the damage and
oxidative stress appears, modifying the metabolism of coordination effects of heavy metals on photosystem
cellular macromolecules. II proteins (Nahar and Tajmir-Riahi, 1996), to analyze
The common methods for monitoring the role of the metallothioneins (MT) (Coyle et al.,
parameters like oxidation of lipids in fish population 2002; Formigari et al., 2007; Vergani et al., 2005), to
are based on chemical analysis. The oxidative study the enzymatic activities in the detoxification
deterioration of lipids in fish oil were monitored using mechanism in response to the presence of different
reflectance NIR and absorbance MIR spectroscopy pollutants in the aquatic environment (Henczová et
accompanied by PLS regression. Hydrolytic al., 2008). Other advantages of FTIR spectroscopy
deterioration of lipids from fish oil, which seriously were to detect abnormal changes in cells and tissues
affects the use of oil and industry storage with minimal preparation under in vivo conditions
requirements, may be successfully monitored by during minimal invasive investigation (Krafft et al.,
spectroscopy and chemometric methods in fishmeal 2004).
industry (Cozzolino et al., 2005). Additionally, it was The infrared spectroscopy has been used by
also possible to achieve reliable estimates of fat Laurens and Wolfrum (2011) to develop the
content using the middle infrared signal intensity multivariate chemometric model for estimating levels
from 1744 cm-1 (Sun, 2009). of spiked neutral and polar fats from micro algae.
Other methods currently available for Benefit of infrared spectroscopic techniques is the
screening lipid content in different sample are based nature straight, fast and non-destructive method of
on fluorescent lipophilic dyes such as Nile red, thin screening for specific characteristics of a large
layer chromatography (TLC), and mass spectroscopy number of samples bypassing the need for long and
(MS). UV-VIS Diffuse Reflectance Spectroscopy laborious wet chemical analyses.
(Arsene et al., 2010), Fourier Transform Infrared The calibration and validation model provides
(FTIR) (Harja et al., 2010) are non-invasive a quick way and high throughput to determinate the
alternative techniques to study the molecular structure content of lipid, offering an alternative to other
and composition modifications concerning a large methods. A multivariate partial least squares (PLS)
variety of samples. A mathematical model to predict like chemometric approach was used to correlate the
the oxidative values of Menhaden fish oil, the main spectral changes during increasing
artificial neural network analysis (ANN) (Gosav et phospholipids’ and triglyceride content in algal
al., 2007), was used by Klaypradit to evaluate biomass collected and to predict the levels of
oxidative quality with ATR technique (Klaypradit et exogenously added lipids of trilaurin as a triglyceride
al., 2010). ATR spectroscopy along with ANN proves and phosphatidylcholine as a phospholipids (Laurens
a rapid technique instead of the conventional method and Wolfrum, 2011).
for prediction of fish lipids oxidation. Levels of antioxidant protection and lipid
Heavy metals (Hlihor and Gavrilescu, 2009; peroxidation were evaluated in different species of
Stefanescu et al., 2011) present in some industrial marine organisms (Almeida et al., 2004). Mussels
wastewaters, despite the essential role of some of the exposed to Cd, Cu, Fe, Pb for different time periods
metals in a number of enzymatic processes are increased lipid peroxidation that could share in low
potentially toxic to freshwater environments levels of reduced glutathione (GSH) and glutathione
(Teodosiu et al., 2009) because of toxic effects on the peroxidase (GPX) activity and depletion of GSH and
receiver, influencing the performance of biological caused GPX activity increased. The correlation
treatment processes (Aksu, 2005). Another FTIR coefficients between the enzyme phospholipid
spectroscopy application was to detect the damaging hydroperoxide glutathione peroxidase PHGPx and
effects of the heavy metals on carp (Cyprinus carpio malonaldehyde levels were performed with statistical
L.), silver carp (Hypopthalmichtys molitrix V.) and software. Negative correlation suggested that enzyme
wels (Silurus glanis L.). To detect the damaging PHGPx was used as a biomarker of the toxicity of
effects of the heavy metals, the investigated fish new potential pollutants related with exposure to the
isoenzymes showed that the Cu2+ was a stronger mussels. Peroxidase lipid is correlated with
inhibitor than Pb2+ (Henczová et al., 2008). antioxidant activity classical systems (Almeida et al.,
The high Cu2+ concentrations appear to be 2004).
toxic to fish, especially if it accumulates beyond the Therefore, correlating macroscopic with
chelating capacity of the cell, generating oxidative microscopic and structural properties is a key element
stress and cellular damage (Gaetke and Chow, 2003; in models development. Chemometrical processing of
Henczová et al., 2008). The continuous exposure of measured data allows the detection of the most
marine organisms to Pb2+ low concentrations may important spectral features, which are crucial for
result in their bioaccumulation and subsequent ranking the fish samples into several groups
transfer to humans through the food chain, owed for according to their specific geographical region.
110
Effects of the environmental stress on two fish populations revealed by statistical and spectral analysis
In the presented study such chemometric using a ZnSe reflectance cell. Samples were carefully
analysis methods are performed using Unscrambler X homogenized and were pressed against the ZnSe cell
10 and SPSS 17.0 software to obtain information with prior to scanning.
minimum error from a reduced number of fish A total of 32 scans were taken for each
samples whose macroscopic characteristics are spectrum, for three replicate samples. The spectra
known. The uses of chemometric method in were collected in the range of 4000 to 600 cm−1 at 1
combination with FTIR spectroscopy are cm−1 resolution and data were exported using Opus
recommended for the discrimination of two fish 6.0 software. Infrared measurements were performed
populations. The main objective of this article is to on samples taken from muscles and brain, carefully
interpret the IR and UV-VIS absorption spectra with homogenized. Background spectra, which were
deconvoluted absorption peaks highlighting the collected under identical conditions, were subtracted
differences between the fish samples. The second from the brain and muscle sample spectra
objective is to develop a rapid calibration method to automatically. All measurements were made at the
predict the percent of lipid extracts in different fish laboratory temperature 22°C. Spectra were recorded
species. on samples before extraction, after extraction of fatty
acids and on extracted samples.
2. Material and methods The extraction was performed using a
microwave equipment Milestone ETHOS lab station,
2.1. Biological material a high-performance microwave-accelerated system
with advanced process control. This method provides
Biological samples were collected from two the solvent extraction of fatty acid from the samples
different ecological environments: in a closed vessel device using temperature control
- Oituz River, Oituz collection point: GPS microwave heating. Diethyl ether was used as a
Coordinates: N: 46.15063 E: 26.45176 Altitude: 200 solvent. Extraction takes place at 35°C and 800 W for
m; covered area: 527 m2; type of course: linear; 20 minutes.
structure of the substrate: gravel, sand, silt, riparian; UV-Vis spectra were recorded only for
vegetation: shrubs, herbaceous plants; aquatic extracted samples using a Varian Cary 100
vegetation: periphyton; physicochemical Spectrophotometer. The molecular absorption spectra
characteristics: air/water temperature: of etheric extract from muscle and brain of the
29.10°C/25.30°C, pH: 8.55; fishing date: 03.09.2008; investigated fish population were measured in the
fishing time: 17 h 10 '; region of 200-400 nm. The absorbance values at 200
- Mureş River, Mureş Basin, Vidrasău collection wavelengths were used as spectral variables for
point: GPS coordinates: N: 46.28877 E: 24.25147 statistical analysis.
Altitude: 285 m; covered area: 228 m2; type of
course: sinuous; structure of the substrate: gravel, 2.3. Spectral processing of middle infrared and UV-
sand; riparian vegetation: trees, shrubs, herbaceous VIS experimental data
plants, reed; aquatic vegetation: periphyton;
physicochemical characteristics: air/water The FTIR instrument provides software for
temperature: 28.10°C/20.10°C, pH: 8.278; fishing spectral subtraction, smoothing, derivative and
date: 10.09.2009; fishing time: 16 h 50 '. deconvolution analysis. The availability of software
The biological material was collected by plays an important role in analyzing of subtle changes
electroshock narcosis, using an electric fishing in infrared biomolecules (Barth and Haris, 2009). For
machine with double insulation EGF 5000 (EFKO - FTIR models the mathematical pre-treatment is an
Elektrofischfanggerate GmbH, Germany). The improvement over the models using the raw spectra.
collected material was fixed with formalin 4%. In our study the FTIR spectral frequency has
Conservation has been achieved in technical alcohol, been achieved with the Opus Bruker software. The
water and glycerin in a ratio of 1: 2: 0.25 (Nicoară et obtained data were exported to Origin 8.0 software.
al., 2006). The mathematical treatment of spectral data was
In Romania the Squalius cephalus L., 1758 made before making the calibration models for
(chub) species lives almost exclusively in small rivers chemometric analyzed. This pre-treatment can help to
and rarely in ponds. In Europe chub live north of the reduce spectral variation due to instrument or sample
Pyrenees, Alps, Dinaric and Balkan. The Barbus variability (Laurens and Wolfrum, 2011).
petenyi Heckel, 1847 species lives only in rivers and Attention is required because the mathematical
streams from mountains and the upper region of hills. pre-treatment of the raw spectra affects the quality of
Both species live in rocky rivers, fast and cold, and in the prediction models. All spectra were smoothed
some streams where summer heat is strong, but only with Savitzky-Golay method at 9 points of window
in the mountains. and 2 polynomial orders. The areas of the absorption
intensity of peaks which are directly related to the
2.2. FTIR and UV-VIS spectroscopic method concentration of the molecules were calculated using
the base line correction (Venkataramana et al., 2010).
FTIR attenuated total reflectance spectra were The fitting curve for FTIR and UV-VIS
collected with a Tensor 27 Bruker spectrophotometer, spectra data were made with Gaussian peaks type.
111
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
UV-VIS spectra were smoothed with Savitzky-Golay measuring parameters and electrical noise. Among
method at 15 points of window and polynomial of these factors, the greatest influence is given by
order 3. R2 value of fitting curve for all samples were changes of concentrations of mixture constituents.
≈ 0.9999 value and Reduced Chi Sq ≈ 8·10-5. The Thus, using the complex and sophisticated
value of R2 and Reduced Chi Sq were represented on mathematical methods, which explanation exceeds
each graph. this paper, based on calibration sample spectra helps
achieving a more limited set of principal components
2.4. Multivariate calibration or PCs factors. Representation of the samples in the
PCs space gives us a set of scores and loadings for
Generally, constituents of a sample have each PCs. So, based on calibration sample spectra,
numerous bands that are partially or completely were achieved a more restricted set of main
overlapping. However, based on Lambert - Beer law components or factors that include 90-99% of
it can be taken into account that mixtures’ absorbance variables values.
is additive. Therefore, for any two components m and
n, the equations bases of the calibration process are 2.4.2. Partial Least Squares (PLS)
Eqs. (1) and (2) (Neamtu and Anoaica, 2006). The Partial Least Squares (PLS) method
performs correlations between spectrum
X i ki mCm ki nCn (1) decomposition and concentration of mixture’s
components. The method has two different
X j k j m C m k j n Cn (2) algorithms: PLS1 and PLS2. PLS1 calculates one set
of spectral factors and scores for each mixture’s
constituent of interest. PLS2 provides a spectral
where: kλ is a constant which includes the extinction
decomposition for all constituents in the mixture. This
coefficient and the specific optical path way for each
means that in this case, it will be obtained a set of
component m and n, and λi, λj are the wavelengths
spectral factors and one of scores for whole set of
corresponding to absorption bands and Cm, n is the
calibration spectra. The main advantages of the PLS
concentration of m and n constituents that have
methods are: the decomposition and regression takes
numerous partially or totally overlap bands.
place in one step, derived spectral factors are directly
Using data measured at multiple wavelengths
linked to constituents of the mixture.
for a great number of samples with known properties,
The method can be applied to highly complex
there is a possibility to minimize the errors, so the
mixtures, in which it is enough to know only the
measurement results are almost near the real value.
concentration of calibration constituents. The
Specific calculation methods consist of writing the
calibration model can be used to predict contaminated
linear algebra matrix of the system equations
samples with constituents that were not present in the
characterizing the complex sample analyzed in the
calibration process. The advantage of these methods
infrared or in the ultraviolet-visible region.
is the possibility of using the wavelengths in the
Multivariate spectrometry, name also
entire range in the calibration method, thus achieving
chemometric method, is based on the specific spectral
an effect of maintenance on the experimental results,
mathematical dependence of absorption intensity of
enhancing the accuracy of the calibration model.
elements, radicals, functional groups or molecules
from simple compounds or complex mixtures of
2.5. Correlations and regression analysis
compounds and their concentrations. This method
applies particularly when the theory available at that
The chemometric method is a less expensive
time is not sufficient for solving problems that have
approach because it is not necessary to work with the
many unknown variables in relationships.
corresponding standards to confirm the outputs of
The chemometric strategy consists in
qualitative analysis performed by a chosen spectral
following steps:
method. It is not necessary to assign all spectral bands
- data acquisition;
or peaks. In this study, the correlations between the
- generate a mathematical model that is based on
percent of lipid extract and the relative intensity area
statistical methods or neural networks;
of middle infrared spectral regions were investigated.
- interpretation of model parameters achieved;
For this purpose the Bivariate Correlations procedure
- applying the model to a new situation.
was computed and Pearson's correlation coefficient
was calculated with SPSS 17.0 software.
2.4.1. Principal Component Analysis PCA
These variables are negatively or positively
The method of Principal Components Analysis
correlated and the correlation is significant at the 0.05
(PCA) underlies chemometric analysis of several
level. General aim of multiple regressions is to
quantitative and qualitative spectrophotometric
observe the relationship between several dependent
methods.
and independent variable.
Note that a sample spectrum is obtained as a
In this case the dependent variables are:
result of various factors changes such as constituents
phosphatidylcholine, deformation of carbohydrates,
of mixture, interactions between constituents, change
sphingolipids&oligosaccharides, PO2- stretching
of external conditions, variations of instruments
symmetric mainly phospholipids, PO2- stretching
112
Effects of the environmental stress on two fish populations revealed by statistical and spectral analysis
asymmetric mainly phospholipids, C=O stretching mathematical pre treatment. The position of MIR
(symmetric) fatty acids and amino acids, CH2 bending peaks presented in literature, the calculated position
mainly lipids, C=O stretching esther phospholipids, of main absorption bands, and their assignments are
CH2 symmetric stretching in saturated fatty acid presented in Table 1.
chain, CH3 symmetric stretching mainly in proteins, 3.2. FTIR spectra samples after mathematical
C-H stretching of C-H in methine groups, CH2 processing
asymmetric stretching in saturated fatty acid chain,
CH3 asymmetric stretching in saturated fatty acid Lipids play a major role as energy provider for
chain regions infrared bands and the independent organism exposed to stress conditions (Jagadeesan et
variable is etheric extract percent. Regression al., 2005). The infrared spectrum of proteins is
statistics and related plots for different curves of characterized by a set of absorption region named as
estimation regression models were made. A separate amide I (1618-1390 cm-1), amide II (1540-1550 cm-1)
model for each dependent variable was realized. and amide III (1240-1310 cm-1) (Naumann et al.,
2009; Sun, 2010). For muscle and brain before and
3. Results and discussions after extraction, and for extract samples, the protein
bands were founded at 1633-1687 cm-1 (Amide I),
3.1. Spectral band assignments of brain and muscle 1541-1545 (Amide II) and 1278-1312 cm-1 (Amide
sample III). Generally, the muscle tissue (Fig. 1) stores
energy rich molecules like glycogen which is a
FTIR spectroscopy was applied to monitor the glucose polymer. In our study, the glycogen
structural changes in lipids and in amino acids of absorption due to C-O, C-C, C-O-H, C-O-C
proteins (Iwaki et al., 2003), to study the damage and deformation of carbohydrates was assigned to 1028-
coordination effects of heavy metals on photo system 1045 cm-1 region.
II proteins (Nahar and Tajmir-Riahi, 1996), to analyze
the role of the metallothioneins (MT) (Coyle et al.,
2002; Formigari et al., 2007; Vergani et al., 2005), to
study the enzymatic activities in the detoxification
mechanism in response to the presence of different
pollutants in the aquatic environment (Henczová et
al., 2008). Other uses of FTIR spectroscopy were to
detect abnormal changes in cells and tissues with
minimal preparation under in vivo conditions during
minimal invasive investigation (Krafft et al., 2004).
Traditional lipid and fatty acid analyses
require relatively large amounts of biomass (>1 g of
dry biomass), it takes a large amount of time and it is
not recommended for the analysis of a large number
of samples. A high number of variables can have a Fig. 1. FTIR spectra of Oituz chub muscle before extraction
potential impact on the lipid content and composition (solid line) and Mures chub muscle
of the fish, influencing the shape of IR absorption before extraction (dot line)
spectra. For fish IR spectra Laurence and Wolfrum
(2011) describes three distinct absorption bands The zones with an important lipid contribution
which are mainly characteristic for lipids. Two of for our study were assigned at 1052-1240 cm-1, 1718-
them: the CH3/CH2 group is identified from 3025 to 1740 cm-1 and 2853-2963 cm-1. The cellular
2954 cm−1 and the C=O esters group is assigned from components are usually not independent, but interact
1746 to 1654 cm−1 in middle infrared region (MIR). with each other. Such examples are the lipoproteins
The hydroxyl and phosphate groups from the and the biological membranes. Therefore,
phospholipids, the last band, can be found in 1200– mathematics treatment of spectral data is required to
500 cm−1 regions. Arrondo and Goni, 1998 shows that separate the contributions of each cellular component.
phospholipids’ infrared band in the spectrum splits In order to highlight the possible differences at
into three different absorption bands corresponding to molecular level between the two fish populations
the polar (1000-1300 cm-1), interfacial (1700-1750 provided by variations of ecosystems, the infrared
cm-1) and hydrophobic (2800-3100 cm-1) moieties of spectra of these samples were analyzed in detail,
the lipid. using five regions of interest. The regions were
In our previous studies (Lazar et al., 2002; deconvoluted in individual peaks, using information
Lazar et al., 2008) a detailed deconvolution of C-H from Table 1. Figs. 2 and 3 show obtained results
simple and double bonds in hydrophobic region is from multi-peak analysis for four of these regions.
presented. In the present study, the complex spectrum Comparing the two samples (Mures and Oituz
which contains several bands arising from chub muscle before extraction), significant
contributions of various functional groups belonging differences in regions 1000-1150 cm-1 (Figs. 2a and
to carbohydrates, proteins and lipids was analyzed. 3a), 1500-1800 cm-1 (Figs. 2c and 3c) and 2800-3000
For accurate data, all spectra were undergoing to cm-1 (Figs. 2d and 3d) can be observed.
113
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
Table 1. FTIR spectral band assignments of brain and muscle samples before, after extraction and for extract samples
114
Effects of the environmental stress on two fish populations revealed by statistical and spectral analysis
a) Multi-peak analysis for 950-1150 cm-1 region b) Multi-peak analysis for 1150-1350 cm-1 region
c) Multi-peak analysis for 1450-1800 cm-1 region d) Multi-peak analysis for 2800-3000 cm-1 region
Fig. 2. Multi-peak analysis for four zones (a- 950-1150 cm-1 ; b-1150-1350 cm-1 ; c- 1450-1800 cm-1; d- 2800-3000 cm-1)
of FTIR spectra of Mures chub muscle before extraction
a) Multi-peak analysis for 950-1150 cm-1 region b) Multi-peak analysis for 1150-1350 cm-1 region
c) Multi-peak analysis for 1450-1800 cm-1 region d) Multi-peak analysis for 2800-3000 cm-1 region
Fig. 3. Multi-peak analysis for four zones (a- 950-1150 cm-1 ; b-1150-1350 cm-1 ; c- 1450-1800 cm-1; d- 2800-3000 cm-1)
of FTIR spectra of Oituz chub muscle before extraction
115
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
So, these differences between two fish A comparison of the absorption spectrum of
populations are mainly attributed to lipid bands, in Oituz and Mures chub brain etheric extract revealed
concordance with studies of Mecozzi et al. (2009). the bathochromic shift for Mures samples.
116
Effects of the environmental stress on two fish populations revealed by statistical and spectral analysis
Table 3. List of FTIR band assignment 3.5. Matrix plot analysis of various infrared spectral
regions
In order to confirm the differences between the
Absorption band (cm-1) Functional group infrared spectra of analyzed samples, Unscrambler X
software was used to obtain the matrix plot for
970 phosphatidylcholine
different samples pairs. The results are presented in
deformation of carbohydrates,
1016-1045 Fig. 6. It can be noticed clear differences between
mainly Glycogen absorption
1052-1061 sphingolipids & oligosaccharides samples before extraction (Figs. 6a, 6b, 6c, 6e) and
PO2- stretching symmetric, mainly between etheric extract samples (Fig. 6f). No
1079-1088 significant differences were observed between Mures
phospholipids
PO2- stretching asymmetric, and Oituz chub brain after extraction (Fig. 6d). This
1225-1240
mainly phospholipids confirms the hypothesis that lipids are the main
C=O stretching symmetric of components of cell membrane correlated with the
1395 COO-, fatty acids and amino environmental characteristics.
acids
CH2 bending, mainly lipids with 3.6. Statistical analysis
1453-1456
little contributions from protein
C=O stretching, ester
1708-1740
phospholipids
3.6.1. Correlation loadings
CH2 symmetric stretching in The regions of interest in our infrared spectral
2853 data were found analyzing the variables using the
saturated fatty acid chains
CH3 symmetric stretching, mainly Principal Component Analysis (PCA) method. The
2865-2871 data matrix of infrared spectra was decomposed into
in proteins
CH2 asymmetric stretching in scores, loadings and residuals. These results indicate
2918-2937
saturated fatty acid chains the degrees of correlation between different FTIR
CH3 asymmetric stretching in spectral regions.
2950-2963
saturated fatty acid chains
Table 4. Measured fish’s mean dimensions and comparison between the relative intensity of major IR characteristic bands
117
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
a. Oituz Barbus petenyi brain before extraction and Oituz Barbus b. Mures Barbus petenyi muscle before extraction and Oituz
petheny muscle before extraction Barbus petenyi muscle before extraction
c. Mures chub muscle before extraction and Oituz chub before d. Mures chub brain after extraction and Oituz chub brain after
extraction extraction
e. Mures chub muscle before extraction and Oituz chub before f. Mures chub muscle etheric extract and Oituz chub etheric
extraction CH2 CH3 bonds region extract CH2 CH3 bonds region
The lowest residual variance is found for 6 the positive correlation between C=O stretching
PCs, but we chose to work with a model consisting of esther phospholipids, CH2 symmetric stretching in
4 PCs. The loading plot, also called map of variables saturated fatty acid chain, CH2 asymmetric stretching
(Fig. 7), displays information about the variables in in saturated fatty acid chain and fingerprint region.
the PCA model. Generally, the variables from the right sided
The plot shows that variables from circle are the variables with the higher degree of
fingerprint region (1000-1300 cm-1) have an extreme correlation. This means that the mentioned variables
position to the right of the plot along PC1, indicating are principal contribution in analyzed infrared
a strong correlation each other. Other observation is spectral data.
118
Effects of the environmental stress on two fish populations revealed by statistical and spectral analysis
3.6.2. Regression analysis dimensions are closely positive correlated. Also, the
The regression methods establish relations correlations founded from Fig. 7 are confirmed by
between two blocks of variables. First, it was made a this correlated map of variables. The regression
regression model with PLS2 of following groups of coefficients (BW) are used to calculate the response
variables: dimensions of samples and relative value from the infrared spectral relative intensity at
intensity of different infrared peaks (Fig. 8). The different region. The size of the coefficients gives a
predicted versus measured plot displays the sign of which variables have a significant impact
predictions done during calibration. In this case, the assessment on the variables response. The BW
R-square of calibration model was 0.902126. The coefficients are used for interpretation of the
loading weights plot can be studied from the regression model. In this case, the CH2 asymmetric
regression model, to find correlating variables. It is stretching in saturated fatty acid chain and CH3
observed that the dimension of samples (length, asymmetric stretching in saturated fatty acid chain are
width, head length, height, thickness and weight) are relevant for macroscopic properties, like dimensions
strongly negatively correlated with PO2- stretching data (Fig. 10).
symmetric mainly phospholipids, PO2- stretching
asymmetric mainly phospholipids, C=O stretching 3.6.3. Correlation analyzed with SPSS software
(symmetric) fatty acids and amino acids, CH2 bending The Bivariate Correlations from SPSS
mainly lipids corresponds to 1000-1500 cm-1 region software calculates the Pearson's correlation
(Fig. 9). coefficient. The Pearson's correlation coefficient is
Dimensions are correlated with fingerprint significant if it is greater than 0.05.
region (1000-1500 cm-1). As is expected, fish’s
Fig. 8. Predicted versus measured plot for regression model calibration between dimensions of samples and relative intensity of
different infrared peaks
119
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
Fig. 9. Correlated map of the variables: dimensions of samples and infrared spectral relative intensity at different region
4. Conclusions
120
Effects of the environmental stress on two fish populations revealed by statistical and spectral analysis
Fig 11. Regression plot between C=O stretching esther and etheric extract percent
Fig 12. Regression plot between CH2 symmetric stretching in saturated fatty acid chain and etheric extract percent
Fig 13. Regression plot between CH2 asymmetric stretching in saturated fatty acid chain and etheric extract percent
121
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
Table 5. The Pearson's correlations of etheric extract percent and each functional group from infrared region
Table 6. ANOVA summary table for C=O stretching esther, CH2 symmetric stretching in saturated fatty acid chain and CH2
asymmetric stretching in saturated fatty acid chain variables
In our study, the chemometrical discover the relations between spectral, biological and
characterization of fish population’s sample, the chemical sample properties. Mathematical
Correlation Analysis and the Regression Analysis as pretreatments can contribute to reduce spectral
the main multivariate techniques made it possible to variation (Laurens and Wolfrum, 2011).
122
Effects of the environmental stress on two fish populations revealed by statistical and spectral analysis
The important spectral regions were assigned Cozzolinoa D., Murraya I., Chreeb A., Scaifec J.R., (2005),
using the regression coefficients, indicating which Multivariate determination of free fatty acids and
infrared regions of the spectra are the main moisture in fish oils by partial least-squares regression
and near-infrared spectroscopy, LWT Food Science and
responsible in calibration model of macroscopic and
Technology, 38, 821-828.
spectral variables. Dreissig I, Machill S., Salzer R., Krafft C., (2009),
The most important differences between the Quantification of brain lipids by FTIR spectroscopy
muscle sample’s spectra from River and Oituz and partial least squares regression, Spectrochimica
populations were attributed to the polar (1000-1300 Acta Part A, 71, 2069-2075.
cm-1) and hydrophobic (2850-2960 cm-1) moieties of Gaetke L.M., Chow C.K., (2003), Copper toxicity,
the lipid components. oxidative stress and antioxidant nutrients, Toxicology,
In the case of brain extract from River and 189, 147-163.
Oituz populations, significant differences appear only Gosav S., Praisler M., Dorohoi D.O., (2007), ANN expert
for interfacial (1700-1750 cm-1) moieties of the lipid system screening for illicit amphetamines using
molecular descriptors, Journal of Molecular Structure,
components. The infrared relative intensity positively 834–836, 188-194.
correlated with etheric extract percent parameter at Formigari A., Irato P., Santon A., (2007), Zinc, antioxidant
1708-1740 cm-1, 2853 cm-1 and 2918-2937 cm-1 was systems and metallothionein in metal mediated-
revealed by statistical analysis of the infrared spectral apoptosis: biochemical and cytochemical aspects,
bands. Note that in the correlation matrix, all infrared Comparative Biochemistry and Physiology Part C:
variables are dependent. Toxicology & Pharmacology, 146, 443-459.
The proposed goals of this study, to use Harja M., Ciocîntă R.C.1, Barbuta M., Lăcrămioara Rusu
attenuated total reflectance (ATR)-FTIR and UV-Vis L., Simion A.I., Bistricianu I.L., (2010), CaCO3
spectroscopy as an alternative method compared with Controllable synthesis by double exchange method
using CaCl2 residual solutions, Environmental
chemical method for evaluating the difference Engineering and Management Journal, 9, 1571-1577.
between two fish population and to apply Henczová M., Deér A.K., Filla A., Komlósi V., Mink J.,
chemometric analysis to predict the etheric extract (2008), Effects of Cu2+ and Pb2+ on different fish
values were achieved. species: Liver cytochrome P450-dependent
In conclusion, the use of spectral information, monooxygenase activities and FTIR spectra,
correlated with statistical methods, allows Comparative Biochemistry and Physiology, Part C,
highlighting the environmental influence on lipid 148, 53-60.
metabolism, for the same fish species from different Hlihor R.M., Gavrilescu M., (2009), Removal of some
hydrographic areas. environmentally relevant heavy metals using low-cost
natural sorbents, Environmental Engineering and
Management Journal, 8, 353-372.
References Jagadeesan G., Kavitha A.V., Subashini J., (2005), FT-IR
Study of the influence of Tribulus terrestris on mercury
Aksu Z., (2005), Application of biosorption for the removal intoxicated mice, Mus musculus liver, Tropical
of organic pollutants: A review, Process Biochemistry, Biomedicine, 22, 15-22.
40, 997-1026. Klaypradit W., Kerdpiboon S., Singh R.K., (2010),
Almeida E.A., Miyamoto S., Bainy A.C.D., Medeiros Application of Artificial Neural Networks to Predict
M.H.G., Mascio P., (2004), Protective effect of the Oxidation of Menhaden Fish Oil Obtained from
phospholipid hydroperoxide glutathione peroxidase Fourier Transform Infrared Spectroscopy Method,
(PHGPx) against lipid peroxidation in mussels Perna Food Bioprocess Technology, 4, 475-480.
perna exposed to different metals, Marine Pollution Krafft C., Stephan B. Sobottka S.B., Schackert G., Salzer
Bulletin, 49, 386-392. R., (2004), Analysis of human brain tissue, brain
Arrondo J.L.R., Goni F.M., (1998), Infrared studies of tumors and tumor cells by infrared spectroscopic
protein-induced perturbation of lipids in lipoproteins mapping, Analyst, 129, 921-925.
and membranes, Chemistry and Physics of Lipids, 96, Lankmayr E., Mocak J., Serdt K., Balla B., Wenz T.,
53-68. Bandoniene D., Gfrerer M., Wagner S., (2004),
Arsene D., Catrinescu C., Drăgoi B., Teodosiu C., (2010), Chemometrical classification of pumpkin seed oils
Catalytic wet hydrogen peroxide oxidation of 4- using UV–Vis, NIR and FTIR spectra, Journal of
chlorophenol over iron-exchanged clays, Biochemical and Biophysical Methods, 61, 95-106.
Environmental Engineering and Management Journal, Laurens L.M.L., Wolfrum E.J., (2011), Feasibility of
9, 7-16. Spectroscopic Characterization of Algal Lipids:
Barth A., Haris P.I. (2009), Biological and Biomedical Chemometric Correlation of NIR and FTIR Spectra
Infrared Spectroscopy, In: Advances in Biomedical with Exogenous Lipids in Algal Biomass, Bioenergy
Spectroscopy, Volume 2, IOS Press BV Nieuwe Research, 4, 22-35.
Hemweg 6B 1013 BG Amsterdam Netherlands. Lazar G., Vascan I., Lazar I., Stamate M., (2002), Some
Brereton R.G., (2007), Biological and Medical Applications properties of sputtered a-C:H, Journal of Non-
of Chemometrics, In: Applied Chemometrics for Crystalline Solids, 299–302, 835-839.
Scientists, John Wiley & Sons Ltd., The Atrium, Lazar G., Bouchet-Fabre B., Zellama K., (2008), Structural
Southern Gate, Chichester, West Sussex PO19 8SQ, properties of nitrogenated amorphous carbon films:
England, 287-317. Influence of deposition temperature and
Coyle P., Philcox J.C., Carey L.C., Rofe A.M., (2002), radiofrequency discharge power, Journal of Applied
Metallothionein: The multipurpose protein, Cellular Physics, 104, 073534, doi.org/10.1063/1.2951622.
and Molecular Life Sciences, 59, 627-647.
123
Lazar et al./Environmental Engineering and Management Journal 11 (2012), 1, 109-124
Mecozzi M., Pietroletti M., Gallo V., Conti M.E., (2009), Palaniappan PL.RM., Nishanth T., Renju V.B., (2010),
Formation of incubated marine mucilages investigated Bioconcentration of zinc and its effect on the
by FTIR and UV–VIS spectroscopy and supported by biochemical constituents of the gill tissues of Labeo
two-dimensional correlation analysis, Marine rohita: An FT-IR study, Infrared Physics and
Chemistry, 116, 18-35. Technology, 53, 103-111.
Moisescu M.G., Leveque P., Bertrand J.R., Kovacs E., Mir Sun D.W., (2009), Infrared Spectroscopy for Food Quality
L.M., (2008), Microscopic observation of living cells Analyses and Control, 1th Edition, Elsevier Inc.,
during their exposure to modulated electromagnetic Oxford, UK.
fields, Bioelectrochemistry, 74, 9-15. Stefanescu I.A., Gavrila L., Mocanu R.D., Olariu R.I.,
Movasaghia Z., Rehmanb S., ur Rehman I., (2008), Fourier Arsene C., (2011), Bioremediation perspective of
Transform Infrared (FTIR) Spectroscopy of Biological Bacillus Megaterium towards heavy metals in
Tissues, Applied Spectroscopy Reviews, 43, 134-179. environments enriched with phosphogypsum,
Nahar S., Tajmir-Riahi H.A., (1996), Complexation of Chemistry Magazine, 62, 245-249.
heavy-metal cations Hg, Cd, and Pb with proteins of Teodosiu C., Cojocariu C., Musteret C.P., Dascalescu I.G.,
Psii—Evidence for metal–sulfur binding and protein Caraene I., (2009), Assessment of human and natural
conformational transition by FTIR spectroscopy, impacts over water quality in the Prut river basin,
Journal of Colloid and Interface Science, 178, 648- Romania, Environmental Engineering and
656. Management Journal, 8, 1439-1450.
Naumann D., Fabian H., Lasch P., (2009), FTIR Venkataramana G.V., Kumar J. K., Prasad A.G.D., Karimi
Spectroscopy of Cells, Tissues and Body Fluids, In: P., (2010), Fourier Transform Infrared spectroscopic
Advances in Biomedical Spectroscopy, vol. 2, IOS study on liver of freshwater fish Oreochromis
Press BV, Amsterdam, Netherlands. Mossambicus, Romanian Journal of Biophysics, 20,
Neamtu J., Anoaica P.G., (2006), Applications of 315-322.
Electromagnetic Radiations in The Medical Field, (in Vergani L., Grattarola M., Borghi C., Dondero F., Viarengo
Romanian), University of Medicine Publishing House, A., (2005), Fish and molluscan metallothioneins, FEBS
Craiova, Romania. Journal, 272, 6014-6023.
Nicoară M., Ureche D., Plăvan G., Erhan M., (2006), Wolkers W.F., (2009), Biomedical FTIR Spectroscopy of
Comparative analyses of the food spectra of the fish Lipids In: Biological and Biomedical Infrared
populations living in the Rivers Oituz and Casin, Spectroscopy, Barth A. and Haris P.I. (Eds.), IOS
Annals Scientific of the University "Al.I. Cuza" Iaşi, Press, 272-287.
Section I: Animal biology, LII, 87-93. Yoshida S., Miyazaki M., Sakai K., Takeshita M., Yuasa S.,
Palaniappan PL.RM., Vijayasundaram V., (2008), FTIR Sato A., Kobayashi T., Watanabe S., Okuyama H.
study of arsenic induced biochemical changes on the (1997), Fourier Transform Infrared spectroscopic
liver tissues of fresh water fingerlings LABEO analysis of rat brain microsomal membranes modified
ROHITA, Romanian Journal of Biophysics, 18, 135- by dietary fatty acids: possible correlation with altered
144. learning behavior, Biospectroscopy, 3, 281–290.
124