|Home | About | Journals | Submit | Contact Us | Français|
Exploring the capabilities of instrumental techniques for discriminating n-3 rich oils derived from animals is a very important though much neglected area that was emphasized more than 100 years ago. In this study the potential of gas chromatography (GC) for discriminating full fatty acid methyl ester (FAME) profiles from fish (cod liver and salmon) and marine mammal (seal and whale) oils is evaluated by means of principal component analysis (PCA). The FAME profiles from plant oils such as rapeseed, linseed and soy oils and seven different brands of n-3 supplements are also used in the discrimination process. The results from the PCA plots can reliably distinguish between plant, n-3 supplements, fish and marine mammal oils. By removing the contribution of the n-3 supplements and plant oils it is possible to discriminate between types of fish and marine animal oils. GC offers a rapid, simple and convenient means of discriminating oils from different species, brands and grades.
The importance of developing techniques aiming at detecting adulteration of fish oils has been emphasized since the late nineteenth early twentieth century when a great scarcity of cod liver oil accompanied by famine prices on the market brought about adulteration of genuine cod liver oil with low-grade shark oil [1, 2]. An overview of the different instrumental techniques used for the analysis of plant, fish and marine mammal oils in studies where the terms discrimination, adulteration, classification, profiling, differentiation, authentication or characterization have been a vital and important component of the various studies revealed that the focus has been on plant oils [3–41] by using mainly different types of chromatography, infrared absorption, nuclear magnetic resonance and mass spectrometry techniques. Studies on fish and marine mammal oils are scarce [41–47]. In addition, the literature overview revealed that multivariate techniques, especially principal component analysis (PCA), were the preferred data analysis approaches used in the majority of these oil discrimination studies [3, 5, 8, 10–14, 16, 17, 19–22, 25–28, 30, 34–38, 40, 42–44, 46, 47] and that the potential of simple techniques such as gas chromatography (GC) has not been explored yet in the discrimination of fish and marine mammal oils by using the conventional fatty acid methyl ester (FAME) profiles.
The previous observations seem to indicate that the issue about the importance of techniques for detecting adulteration in n-3 rich oils derived from animals has been overlooked by practitioners during the course of a century and it is a topic of contemporary relevance that needs attention.
The main goal of this study is to evaluate the capability of GC to discriminate between n-3 rich oils derived from fish and marine mammals of different species, brands and grades by using the FAME profiles and PCA. The GC FAME profiles from plant oils such as rapeseed, linseed and soy oils and seven different brands of n-3 supplements are also used in the discrimination process. The discrimination between and within animal oils is studied by using three different data analysis strategies: firstly, the analysis of the full FAME profiles from plant, supplement and animal oils; secondly, the analysis of selected FAME profiles from plant, supplement and animal oils with levels higher than 0.5% of the total composition; thirdly, the analysis of the full FAME profiles from animal oils. It must be mentioned that discrimination studies of n-3 rich oils derived from fish and marine mammals have not been previously reported.
Sodium hydroxide, hexane, methanol, boron trifluoride in methanol (20% w/v) and chloroform were purchased from Merck (Darmstadt, Germany). Butylated hydroxytoluene (BHT) and boron trichloride in methanol (14%) were purchased from Sigma–Aldrich Co., USA. FAME pure standards and also model mixture standards 2A and 2B (18:0, 18:1n-9, 18:2n-6, 18:3n-3, 20:4n-6), 3A (18:2n-6, 18:3n-3, 20:4n-6, 22:6n-3), 4A (6:0, 8:0, 10:0, 12:0, 14:0) 6A (16:0, 18:0, 20:0, 22:0, 24:0), 7A (16:1n-7, 16:1n-9, 20:1n-9, 22:1n-11, 24:1n-9) and 14A (13:0, 15:0, 17:0, 19:0, 21:0) were purchased from Nu-Chek Prep (Elysian, MN, USA). Nonadecanoic acid methyl ester (19:0) internal standard was from Fluka (Buchs, Switzerland). De-ionized water was purified in a Milli-Q system (Milli-Q system Millipore, Milford, MA, USA). The fish oils were cod liver oil from Peter Möller, Lysaker, Norway and salmon oil from Havnegater, Sortland, Norway. The two brands of harp seal oils (Pagophilus groenlandicus) were from Rieber Skinn A/S, Bergen, Norway (two refined samples from different batches, designated as RSA1 and RSA2, and one crude sample designated as CSA were provided) and from JFM Sunile A/S, Os, Norway (one refined sample designated as RSB was provided). Whale oil (Balaenoptera acutorostrata) conventionally (WC) and molecularly (WM) distilled were from Myklebust Trading AS, Myklebost, Norway. The plant oils analyzed were soy oil (Mills DA, Sofienberg, Norway), linseed and rapeseed oils (Kinsarvik Naturkost, Bergen, Norway). The seven commercial n-3 supplements obtained from a local pharmacy were Fri Flyt (Vesterålens Naturprodukter AS, Sortland, Norway), Natur-Omega (Naturhuset AS, Vøyenenga, Norway), Møllers dobbel (Peter Möller, Lysaker, Norway), Pikasol (Axellus A/S, Oslo, Norway), Omega-3 Forte (Vitamed, Sarpsborg, Norway), Omega-3 høykonsentrert (Sunkost, Oslo, Norway), El Dorado (Probio Nutraceuticals, Tromsø, Norway). The supplements were designated as K1, K2, K3, K4, K5, K6 and K7, respectively.
The FAME preparation protocol has been published elsewhere . Briefly, 50 mg of the sample are mixed with 2 ml BF3/CH3OH and 5 mg of the 19:0 internal standard. The mixture is heated at 100 °C for 1 h and cooled down to room temperature. Aliquots of 1 ml of hexane and 2 ml of H2O are added, vortex-mixed for 15 s, placed in a centrifuge at 3,000 rpm for 2 min and the FAME are then extracted from the upper hexane phase. Depending on the fat content the sample is either concentrated under nitrogen or diluted with hexane and subsequently subjected to GC analysis.
The FAME for every kind of oil (fish, marine mammal and plant) were prepared in triplicate and for the n-3 supplements in duplicate.
Analysis of the FAME was performed on a Perkin-Elmer AutoSystem XL gas chromatograph (Perkin-Elmer, Norwalk, Connecticut) equipped with a liquid autosampler and a flame ionization detector. The FAME samples were analyzed on a CP-Sil 88 capillary column (50 m × 0.32 mm ID 0.2 μm film thickness, Varian, Courtaboeuf, France). Data collection was performed by the Perkin-Elmer TotalChrom Data System software version 6.3. The temperature program was as follows: the oven temperature was held at 60 °C for 1 min, ramped to 160 °C at 25 °C/min, held at 160 °C for 28 min, ramped to 190 °C at 25 °C/min, held at 190 °C for 17 min, ramped to 220 °C at 25 °C/min and finally held at 220 °C for 10 min. Direct on-column injection was used. The injector port temperature was ramped instantaneously from 50 to 250 °C and the detector temperature was 250 °C. The carrier gas was ultra-pure helium at a pressure of 82 KPA. The analysis time was 60 min. This time interval was sufficient to detect FAME with chains from 10 to 24 carbons in length. The FAME peaks were identified by comparison of their retention times with the retention times of highly purified FAME standards. Since the concentration of the internal standard (19:0) is known and its recovery in the different oils was constant, the comparison of the various FAME peak areas with 19:0 peak area is used to calculate the concentration of the FAME in the various oils and supplements.
Principal component analysis (PCA) is a well documented multivariate method for reducing the dimensionality of a data set by rotating and constructing orthogonal linear combinations of the original variables and projecting the maximum variability onto a new axis also known as principal components (PCs). The results are classified according to the level of information produced by the various combinations of the original variables in a way that the first component (PC1) is the major axis of the points in the p-dimensional space that accounts for the largest variation in the data and consequently it contains the most possible information. The second component (PC2) is perpendicular to PC1 and it defines the next largest amount of variation.
The results as presented in the various tables (Tables 1, ,2,2, ,3,3, ,44 and and5)5) were combined (according to the various approaches to be discussed) and arranged in m × n data matrices where m represents every prepared oil sample and n represents every analyzed fatty acid. The matrices submitted to PCA are standardized by subtracting their means and dividing by their standard deviations in order to construct linear combinations of the predictor variables n that contains the greatest variance. The PCA score and loading plots of the FAME profiles from the various oils were computed with the software package Statgraphics Plus 5.1 (Statistical Graphics Corp.).
The fish, marine mammal and plant oils were prepared in triplicate and the n-3 supplements in duplicate. The triplicate and duplicate lipid profiles from the various injected oil samples, expressed as mg FAME/g sample, are presented in Tables 1, ,2,2, ,3,3, ,44 and and5.5. The individual profiles were arranged in a data matrix consisting of 47 rows representing the various analyzed oils with their respective replicates and 34 columns representing the individual FAME detected by GC. The 34 individual FAME profiles were: 14:0, 14:1n-9, 15:0, 16:0, 16:1n-9, 16:1n-7, 17:0, 16:2n-4, 18:0, 16:3n-3, 18:1n-11, 18:1n-9, 18:1n-7, 16:4n-3, 18:2n-6, 20:0, 18:3n-3, 20:1n-11, 20:1n-9, 20:1n-7, 18:4n-3, 20:2n-6, 20:3n-6, 22:0, 20:3n-3, 20:4n-6, 22:1n-11, 22:1n-9, 20:4n-3, 20:5n-3, 24:0, 24:1n-9, 22:5n-3, 22:6n-3.
The 47 × 34 matrix was submitted to PCA as a data exploration technique and a total of six components (PCs) were extracted and grouped in decreasing order of variance. The first component (PC1) which explains 41.91% of the total variation can be used to discriminate the oils according to their nature as is shown in Fig. 1. A plot of the scores of the two first components (Fig. 2), which explain 66.35% of the data variation, differentiates basically the same number of groups and sub-groups found in Fig. 1. Besides, the loadings of the two first components, were plotted to investigate the relationship between the various FAME (Fig. 3).
The six PCs computed by using the 47 × 34 matrix were plotted against each other to produce two and three dimensional PC scores graphs and consequently explore the capability of these PCs to discriminate confidently within the marine animal oils. Unfortunately, clear and well-defined patterns that allow one to differentiate the various oils and grades were not observed in any of the graphical representations, hence a data reduction, based on selected FAME profiles, was implemented.
FAME data reduction has been used in the discrimination of oils derived from one fish species (cod liver oil) by selecting the 15 FAME with levels higher than 1% of the total composition . Considering that in the present study the 34 FAME or variables are given in mg/g and arranged in columns for PCA purposes (47 × 34), it was decided to discard all the FAME columns with average values <5 mg/g (<0.5% of the total averaged FAME profile), in that way 19 FAME profiles were retained (14:0, 16:0, 16:1n-7, 18:0, 18:1n-11, 18:1n-9, 18:1n-7, 18:2n-6, 18:3n-3, 20:1n-11, 20:1n-9, 18:4n-3, 20:4n-6, 22:1n-11, 20:4n-3, 20:5n-3, 24:1n-9, 22:5n-3, 22:6n-3). A new data matrix of size 47 × 19 was submitted to PCA and three PCs, explaining 85.29% of the total data variability, were computed. The three PCs were used to generate various two and pseudo-three dimensional score plots which essentially showed a clear discrimination between plant, supplements and marine animal oils as already observed in Figs. 1, ,22.
The contribution of the supplement and plant oils was removed as an alternative approach to discriminate the various marine oils. A 24 × 34 matrix was constructed by using the data presented in Tables 1, ,22 and submitted to PCA. A considerable percentage of the variability of this matrix (97.58%) was explained by four extracted PCs (49.80, 22.65, 20.08 and 5.05%). The score plot of PC1 and PC4 (Fig. 4) indicates that it is possible to discriminate within the four different types of animal oils by using these two components.
Further assessment of the capability of GC for discriminating between pure and adulterated oils, was achieved using conventionally distilled whale oil (WC) debased by cod liver oil in the proportions 90:10, 80:20 and 50:50. The FAME profiles of the mixtures prepared in triplicate (Table 5) along with the profiles of the marine animal oils described in Tables 1 and and22 were submitted to PCA and four PCs extracted, explaining 47.05, 25.08, 20.19 and 5.36% of the data variability, respectively.
A graph of PC1 for the various oil samples (Fig. 1) shows that different kinds of oils can be basically differentiated along the PC1 axis. Specifically, the scores of the analyzed plant and animal oils grouped themselves at opposite ends of the PC1 axis while the scores of six supplements (K1 and K3–7) in the middle of this axis. The scores for K2 overlaps those from rapeseed oil, hence it is likely that this supplement contains this particular oil. In addition, the supplements as well as the animals exhibit some sub-groups which discriminate supplement K3 in the former and salmon oil in the latter. These latent sub-structures are attributed to the consistently high levels of 16:0 and 16:1n-7 in supplement K3 (on average three times higher compared to the others supplements of this group) and the high levels of 18:1n-9, 18:2n-6 and 18:3n-3 in the salmon oil (on average 1.7, 4.7 and 5.0 times higher than the rest of marine animal oils, respectively).
The score plot of the first two PC (Fig. 2) which explained 41.91 and 24.44%, respectively, of the total variability of the different types of oil FAME profiles allowed concluding that PC1 discriminates in effect between animal and plant oils while PC2 differentiates between supplements and plant oils. Besides, this score plot revealed that in addition to rapeseed oil, supplement K2 also contains linseed oil due to the proximity of their scores. This proximity was constantly observed when the scores of any of the six PCs were plotted against each other. The presence of rapeseed and linseed oils in the composition of K2 was confirmed by searching the webpage of the manufacturer of this particular supplement. The PC1 versus PC2 plot (Fig. 2) revealed some variability in the individual scores for K7 which could be attributed to experimental errors, indicating the importance of replication in discrimination studies.
The plot of the loadings of the two first components, expressing the relationship between the various FAME (Fig. 3) showed the lack of correlation between 20:5n-3 and 18:1n-9 on the PC2 axis, indicating that none of the plant oils studied in the present investigation are present in the composition of supplements K1, K3, K4, K5, K6 and K7. These observations are in agreement with the various manufacturers who have reported some special developed oils (name not disclosed), refined fish oil (from non-specified origin), gelatin from pork, etc. among the various constituents of their supplements rather than plant oils.
The superimposition of the three main clusters from Fig. 2 on the loading plot (Fig. 3) demonstrates unequivocally that the observed inverse correlation between 20:5n-3 (positive PC2 loading value) and 18:1n-9 (negative PC2 loading value) is responsible for the discrimination between supplements (with scores highly associated to 20:5n-3) and plant oils (with scores highly associated to 18:1n-9). Similarly, the lack of correlation between 16:1n-9, 14:0, 20:1n-9, 16:2n-4 (positive PC1 loading values) and 22:0, 18:2n-6 (negative PC1 loading values) is responsible for discriminating between animal and plant oils.
The discrimination patterns of the various two and pseudo-three dimensional score plots obtained after performing PCA on selected FAME profiles were basically similar to those obtained by using the full FAME profiles (Figs. 1, ,2).2). However, the graphs consistently misclassified supplement K2 as containing a mixture of the three plant oils (soy, linseed and rapeseed oil) while in fact only linseed and rapeseed oil are present in this particular supplement. This result indicates that PCA on full FAME profiles outperforms the proposed data reduction approach for supplement classification. The various plots generated with the aforementioned three PCs were also unable to establish a clear distinction between the various species, brands and grades of animal oils; hence a further discrimination study was carried out by using only the full FAME profiles derived from marine oils.
The PC1 and PC4 score plot (Fig. 4) shows that PC1 can discriminate between the four different types of animal oils, namely seal oil, cod liver oil, whale oil and salmon oil and that PC4 can discriminate effectively within every animal oil species and their alleged qualities. For instance, it is observed that the two different batches of refined seal oils from manufacturer-A (designated as RSA1 and RSA2 in Fig. 4) display positive and negative PC4 score values, respectively, indicating some differences between them. In addition, the crude and the refined seal oils from the same manufacturer and designated as CSA and RSA2 in Fig. 4, exhibit negative scores values, indicating a correlation between these two oils regardless their alleged qualities. The discrimination between the seal oils from different manufacturers namely Rieber Skinn (SA) and JFM Sunile (SB) is observed in Fig. 4. The variables responsible for the discrimination between the manufactures are visualized by means of the PC1 and PC4 loading plot (Fig. 5). The higher contents of 16:1n-9, 18:1n-7, 20:1n-7 and the lower contents of 16:2n-4 in SB compared to SA are the main fatty acids that contribute to distinguish SA from SB. Similarly, Fig. 4 shows the unmistakably differentiation between molecularly distilled whale oil (WM) from conventionally distilled whale oil (WC) which is mainly due to the slightly higher levels of 20:1n-9, the lower levels of 22:6n-3 and 18:4n-3 as well as the lack of 16:4n-3 in WM compared to WC. The clear-cut distinction between cod liver and salmon oil along the PC1 axis (Fig. 4) is mainly attributed to the contribution of several fatty acids such as 18:1n-11 and 16:1n-7 at negative PC1 values and 20:2n-6 and 18:3n-3 at positive PC1 values as indicated in Fig. 5.
The PC1and PC4 score plot (Fig. 6) demonstrates that it is possible to discriminate adulterated from genuine whale oil samples. The information retained by PC1 is mainly connected with the nature of the various oils. For instance, five specific regions can be visualized along the PC1 axis (salmon, whale, whale + cod liver, cod liver and seal) while the information retained by PC4 is mainly connected with the discrimination within animal oil species, their alleged qualities and the various proportions of cod liver oil used to adulterate genuine whale oil.
In conclusion, the different approaches used in the discrimination process indicated that PCA on the full FAME profiles is the best strategy to discriminate between the various oils considered in this study. Considering that n-3 rich oils derived from animals are highly regarded as alternative medicines worldwide, the potentiality of unexplored single or coupled techniques for authentication and discrimination of these kinds of oils should be investigated to prevent fraudulent practices.
The European Commission in the context of the Erasmus Mundus Program and The Norwegian Research Council (SIP project NRF 173534/I30) are gratefully acknowledged for financial support of Y.Z and Z.D, respectively.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.