JP6692226B2 - Mycotoxin analysis method and mycotoxin analyzer - Google Patents

Mycotoxin analysis method and mycotoxin analyzer Download PDF

Info

Publication number
JP6692226B2
JP6692226B2 JP2016126209A JP2016126209A JP6692226B2 JP 6692226 B2 JP6692226 B2 JP 6692226B2 JP 2016126209 A JP2016126209 A JP 2016126209A JP 2016126209 A JP2016126209 A JP 2016126209A JP 6692226 B2 JP6692226 B2 JP 6692226B2
Authority
JP
Japan
Prior art keywords
analysis
spectrum
mycotoxin
sample
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016126209A
Other languages
Japanese (ja)
Other versions
JP2018004252A (en
Inventor
望月 淳
望月  淳
忠生 國弘
忠生 國弘
宏夫 浜口
宏夫 浜口
安藤 正浩
正浩 安藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
TECHNOSURUGA LABORATORY CO.,LTD.
Original Assignee
TECHNOSURUGA LABORATORY CO.,LTD.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by TECHNOSURUGA LABORATORY CO.,LTD. filed Critical TECHNOSURUGA LABORATORY CO.,LTD.
Priority to JP2016126209A priority Critical patent/JP6692226B2/en
Publication of JP2018004252A publication Critical patent/JP2018004252A/en
Application granted granted Critical
Publication of JP6692226B2 publication Critical patent/JP6692226B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、振動分光分析法を用いるかび毒の分析方法及びかび毒の分析装置に関するものである。   TECHNICAL FIELD The present invention relates to a method for analyzing mycotoxin using vibrational spectroscopy and an apparatus for analyzing mycotoxin.

かびが産生する代謝産物のうちヒトや動物に対する毒性を有するものを、かび毒(マイコトキシン)と呼ぶ。かび毒は、ヒトや家禽、家畜などに対して急性又は慢性の生理的又は病理的障害を引き起こすことが知られている。そのため国際的には、コーデックス委員会(消費者の健康の保護、食品の公正な貿易の確保等を目的として、1963年に国連食糧農業機関(FAO)及び世界保健機関(WHO)により設置された国際的な政府間機関)がかび毒汚染の防止・低減に関する各種の規範等を定めている。また国内では、食品衛生法や農林水産省の各種管理ガイドライン等に基づく規範が定められている。   Among the metabolites produced by fungi, those that are toxic to humans and animals are called mycotoxins. Mycotoxin is known to cause acute or chronic physiological or pathological disorders in humans, poultry, livestock and the like. Therefore, internationally, the Codex Committee (established by the United Nations Food and Agriculture Organization (FAO) and the World Health Organization (WHO) in 1963 for the purpose of protecting consumer health and ensuring fair trade in food, etc. International intergovernmental organizations) have established various norms regarding the prevention and reduction of mycotoxin contamination. In Japan, norms are established based on the Food Sanitation Law and various management guidelines of the Ministry of Agriculture, Forestry and Fisheries.

食品や飼料の生産及び流通のさまざまな現場において、かび毒の許容値や規範の遵守状況が検査及び確認されなければならない。そのためには、できるだけ場所を選ばずにかび毒を簡便かつ迅速に分析する必要がある。分析は少なくとも定性的に行われる必要があり、可能であれば定量的に行われることが望ましい。また、同一試料に含まれる複数種のかび毒を並行してそれぞれ分析できるようにすることが、迅速かつ効率的な分析のために望ましい。   At various sites of food and feed production and distribution, compliance with acceptable values and norms of mycotoxin must be inspected and confirmed. For that purpose, it is necessary to analyze mycotoxins easily and quickly at any place. The analysis should be done at least qualitatively and preferably quantitatively. In addition, it is desirable for rapid and efficient analysis that multiple types of mycotoxin contained in the same sample can be analyzed in parallel.

従来のかび毒の分析は、一般的な化学分析法と同じく、目的物質(かび毒)を含む試料の採取、試料から溶媒への目的物質の抽出、抽出物の精製及び目的物質の検出という手順で行われる。かび毒の検出には、液体クロマトグラフィーやガスクロマトグラフィー等の方法が、蛍光検出器や質量分析計等の検出器と組み合わせて用いられている。しかし、これらの検出器は非常に高価であり、かつ、精製で除ききれなかった試料中の夾雑物により定量値が大きく増減する場合がある(例えば、非特許文献1参照。)。   As in the case of general chemical analysis methods, the conventional analysis of mycotoxins is the procedure of collecting a sample containing the target substance (mold poison), extracting the target substance from the sample to a solvent, purifying the extract and detecting the target substance. Done in. Methods such as liquid chromatography and gas chromatography are used in combination with detectors such as fluorescence detectors and mass spectrometers to detect mycotoxins. However, these detectors are very expensive, and the quantitative value may greatly increase or decrease due to impurities in the sample that cannot be removed by purification (see, for example, Non-Patent Document 1).

液体クロマトグラフィーやガスクロマトグラフィーを用いる従来技術では、測定結果が試料中の夾雑物の条件に左右されやすいという問題があり、また高価な検出器を食品等の生産や流通のさまざまな現場に持ち込むことは必ずしも現実的とはいえない。これに対して、比較的簡単な設備で分析を行える方法として酵素固定化免疫測定法(ELISA)が利用されているが、すべての特異性の高い抗原・抗体反応を1個のウェル内で行い、抗体に標識した酵素に基質を加えて反応、発色させ、その吸光度で結果を得ているために、夾雑物の影響を受けやすく擬陽性又は擬陰性の結果を得るなどの問題がある。   The conventional technique using liquid chromatography or gas chromatography has a problem that the measurement result is easily influenced by the condition of impurities in the sample, and also brings an expensive detector to various sites of production and distribution of food etc. This is not always realistic. On the other hand, enzyme-immobilized immunoassay (ELISA) is used as a method to perform analysis with relatively simple equipment, but all highly specific antigen / antibody reactions are performed in one well. Since a substrate is added to an enzyme labeled with an antibody to cause reaction and color development, and the result is obtained by its absorbance, there is a problem that it is easily affected by contaminants and a false positive or false negative result is obtained.

液体クロマトグラフィーやガスクロマトグラフィーに代わり、光学的手法を用いて煩雑な前処理を不要にするという技術が開示されている(例えば、特許文献1又は非特許文献2参照。)。特許文献1は、かび毒のような有害成分(危害要因)を含む測定対象物へ照射する光の励起波長と測定対象物から観測される蛍光の波長を段階的に変化させながら蛍光強度を測定して、蛍光指紋(又は、励起・蛍光マトリクス)と呼ばれるデータを取得し、当該蛍光指紋データに対して多変量解析を行うことにより危害要因を検知するという技術を開示する。   Instead of liquid chromatography or gas chromatography, a technique has been disclosed in which an optical method is used to eliminate the need for complicated pretreatment (for example, refer to Patent Document 1 or Non-Patent Document 2). Patent Document 1 measures the fluorescence intensity while gradually changing the excitation wavelength of the light irradiating the measurement target containing a harmful component (harm factor) such as mycotoxin and the wavelength of the fluorescence observed from the measurement target. Then, a technique is disclosed in which data called a fluorescent fingerprint (or excitation / fluorescence matrix) is acquired, and a hazard factor is detected by performing multivariate analysis on the fluorescent fingerprint data.

特許文献1に記載された発明によれば、液体クロマトグラフィーやガスクロマトグラフィーを用いる方法よりも迅速かつ容易に危害要因を検知することができるという効果を奏することが記載されている。   According to the invention described in Patent Document 1, it is described that the harmful factor can be detected more quickly and easily than the method using liquid chromatography or gas chromatography.

ただし、特許文献1に記載された方法は蛍光分光分析法の一種であるから、励起波長と蛍光波長の組合せ(波長条件)を変えながら、それらの波長条件ごとに光信号強度のデータを取得する必要がある(例えば明細書の段落「0069」に、5041波長条件が例示されている。)。また、有害物質の濃度の定量には、既知又は予め実測された検量線を用いる必要がある(明細書の段落「0075」)。これらを考慮すると、できるだけ場所を選ばず簡便かつ迅速にかび毒を分析できるようにするには、もう一歩の技術的前進が望まれるといってよい。   However, since the method described in Patent Document 1 is a kind of fluorescence spectroscopic analysis method, while changing the combination (wavelength condition) of the excitation wavelength and the fluorescence wavelength, the data of the optical signal intensity is acquired for each of these wavelength conditions. It is necessary (for example, 5041 wavelength condition is illustrated in paragraph “0069” of the description). In addition, it is necessary to use a known or previously measured calibration curve for quantifying the concentration of the harmful substance (paragraph “0075” in the specification). Taking these into consideration, another technological advance would be desirable in order to be able to analyze mycotoxins easily and quickly in any place.

特開2012−177606号公報JP, 2012-177606, A

田端節子 「講座カビ検査法10 カビ毒分析」 日本防菌防黴学会誌 Vol.43、No.10、pp.485−494、2015年Setsuko Tabata "Lecture Mold Inspection Method 10 Mycotoxin Analysis" Journal of Japanese Society of Antibacterial and Antifungal Vol. 43, No. 10, pp. 485-494, 2015 藤田、蔦、杉山「励起蛍光マトリクス計測を応用したデオキシニバレノールの新規判別法」 日本食品科学工学会誌 第55巻、第4号、pp.177−182、2008年4月Fujita, Tsuta, Sugiyama "A New Method for Deoxynivalenol Discrimination by Applying Excitation Fluorescence Matrix Measurement" Journal of Food Science and Technology, Japan Vol. 55, No. 4, pp. 177-182 April 2008

本願発明が解決しようとする課題は、できるだけ場所を選ばず簡便かつ迅速にかび毒を分析できるようにすることである。また望ましくは、同一試料に含まれる複数種のかび毒を並行してそれぞれ分析できるようにすることである。   The problem to be solved by the present invention is to enable simple and quick analysis of mycotoxin in any place as much as possible. It is also desirable that plural types of mycotoxin contained in the same sample can be analyzed in parallel.

そのためには、光学的手法の中でも、蛍光分光分析のように多くの波長条件に対して計測を繰り返すことを必要としない振動分光分析法(分子の振動に伴う双極子モーメントの変化によって生じる赤外線の吸収を観測する赤外吸収分光分析法、分極率の変化によって生じるラマン効果を観測するラマン分光分析法等)が有望である。振動分光分析法に用いる装置は小型化されたものを入手しやすいので、ある程度場所を選ばずに済む。なお、例えば試料中の夾雑物由来の蛍光のような背景光の影響をできるだけ抑えられることが望ましい。   For that purpose, among the optical methods, the vibrational spectroscopic analysis method that does not require repeated measurement for many wavelength conditions such as the fluorescence spectroscopic analysis (the infrared rays generated by the change of the dipole moment due to the vibration of the molecule) Infrared absorption spectroscopy for observing absorption, Raman spectroscopy for observing Raman effect caused by changes in polarizability, etc.) are promising. The equipment used for the vibrational spectroscopic analysis is easily available in a small size, so that it does not need to be selected in any place. It is desirable that the influence of background light such as fluorescence derived from impurities in the sample can be suppressed as much as possible.

試料からスペクトルのデータが得られたならば、目的物質を含む標準試料から予め取得しておいた参照スペクトルのデータも利用して多変量解析を行い、検量線を実測するのではなく仮想的に求める方法が有望である。このような方法は、本願発明者のうち2者の発明に係る特許第5780476号公報に開示されている。当該方法によれば複数の定量目的成分の組成比を同時に決定できるという効果の得られることが、特許第5780476号公報の明細書の段落「0011」に記載されている。すなわち、複数種のかび毒を並行してそれぞれ分析することが期待される。   Once the spectrum data is obtained from the sample, multivariate analysis is also performed using the reference spectrum data previously acquired from the standard sample containing the target substance, and the calibration curve is virtually measured instead of being actually measured. The way to ask is promising. Such a method is disclosed in Japanese Patent No. 5780476 related to the invention of two of the inventors of the present application. It is described in paragraph “0011” of the specification of Japanese Patent No. 5780476 that the effect of being able to simultaneously determine the composition ratio of a plurality of quantitative target components is obtained according to this method. In other words, it is expected to analyze multiple species of mycotoxin in parallel.

上述した課題を解決するため、本願発明に係るかび毒の分析方法は、ラマン分光分析法を用いるかび毒の分析方法において、分析の目的とする種類のかび毒が産生する可能性がある初期試料を採取し、前記採取した初期試料から、前記分析の目的とする種類のかび毒を含む可能性がある分析用試料を生成し、前記振動分光分析法を用いて、前記分析の目的とする種類のかび毒を含む標準試料から参照用スペクトルを取得し、前記分析用試料に対して、連続する複数の時間インターバルにわたってレーザ光を照射し、前記複数の時間インターバルごとに、前記分析用試料から散射されるラマン散乱光及び前記ラマン散乱光を除く背景光を集光すると共に、前記ラマン散乱光及び前記背景光を含むスペクトルを記録し、前記複数の時間インターバルごとに記録されたスペクトルから、多変量解析を用いて、前記ラマン散乱光の強度及び前記背景光の強度のそれぞれの時間依存性を抽出すると共に前記ラマン散乱光のスペクトルを前記背景光のスペクトルから分離する分離処理を行うことにより、前記分析用試料から分析用スペクトルを取得し、前記参照用スペクトルのデータ及び前記分析用スペクトルのデータを入力とする多変量解析を用いた演算処理を行うことにより、前記分析の目的とする種類のかび毒が前記初期試料に産生する場合において、前記分析用スペクトルから前記分析の目的とする種類のかび毒のスペクトルを抽出することを特徴とする。 In order to solve the above-mentioned problems, the method for analyzing mycotoxin according to the present invention is a method for analyzing mycotoxin using Raman spectroscopic analysis, in which the type of mycotoxin to be analyzed may be an initial sample. From the collected initial sample to produce an analytical sample that may contain a mold venom of the type intended for the analysis, and using the vibrational spectroscopic method, the type of the target of the analysis. Acquiring a reference spectrum from a standard sample containing fungal venom, irradiating the analysis sample with laser light over a plurality of consecutive time intervals, and emitting from the analysis sample at each of the plurality of time intervals. The Raman scattered light and the background light excluding the Raman scattered light are collected, a spectrum including the Raman scattered light and the background light is recorded, and the plurality of time intervals are recorded. From the spectrum recorded for each, using multivariate analysis, extract the time dependence of the intensity of the Raman scattered light and the intensity of the background light and the spectrum of the Raman scattered light from the spectrum of the background light. By performing a separation process to separate , obtain an analysis spectrum from the analysis sample, by performing an arithmetic process using multivariate analysis with the reference spectrum data and the analysis spectrum data as input In the case where the target type of mycotoxin in the analysis is produced in the initial sample, the spectrum of the target type of mycotoxin in the analysis is extracted from the analysis spectrum.

また、本願発明に係るかび毒の分析装置は、ラマン分光分析法を用いてかび毒の分析を行うことができるかび毒の分析装置において、分析の目的とする種類のかび毒が産生する可能性がある初期試料から生成された分析用試料から、前記振動分光分析法を用いて分析用スペクトルを取得することができる分光分析部と、前記分析用スペクトルのデータ及び前記分析の目的とする種類のかび毒を含む標準試料から前記振動分光分析法を用いて取得された参照用スペクトルのデータを入力として、多変量解析を用いた演算処理を行うことにより、前記分析の目的とする種類のかび毒が前記試料に産生する場合において前記分析の目的とする種類のかび毒のスペクトルを抽出することができる演算処理部とを備え、前記分光分析部は、レーザ光を照射された前記分析用試料が散射するラマン散乱光及び前記ラマン散乱光を除く背景光を集光して、前記ラマン散乱光及び前記背景光を含むスペクトルを記録することができる記録部と、連続する複数の時間インターバルにわたって前記分析用試料が連続してレーザ光を照射されたときに、前記複数の時間インターバルごとに記録されたスペクトルから、多変量解析を用いて、前記ラマン散乱光の強度及び前記背景光の強度のそれぞれの時間依存性を抽出すると共に前記ラマン散乱光のスペクトルを前記背景光のスペクトルから分離することにより、前記分析用スペクトルを取得することができる分離部とを備えたことを特徴とする。 Further, the mold venom analyzer according to the present invention is a mold venom analyzer that can perform analysis of the mold venom using Raman spectroscopy, and the type of mold venom to be analyzed may be produced. There is a spectroscopic analysis unit capable of acquiring an analysis spectrum by using the vibrational spectroscopic analysis method from an analysis sample generated from an initial sample, and data of the analysis spectrum and a target type of the analysis. By inputting the data of the reference spectrum obtained by using the vibrational spectroscopic analysis method from the standard sample containing the mycotoxin, by performing the arithmetic processing using the multivariate analysis, the mycotoxin of the kind of the analysis is performed. There and an arithmetic processing unit capable of extracting the spectrum of types of mycotoxin of interest of the analysis in the case of producing the sample, the spectroscopic analysis unit, the laser beam irradiation The Raman scattered light scattered by the sample for analysis and the background light excluding the Raman scattered light are condensed, and the recording section capable of recording the spectrum including the Raman scattered light and the background light is continuous. When the analysis sample is continuously irradiated with laser light over a plurality of time intervals, from the spectrum recorded for each of the plurality of time intervals, using multivariate analysis, the intensity of the Raman scattered light and the By separating each time dependence of the intensity of background light and separating the spectrum of the Raman scattered light from the spectrum of the background light, a separation unit capable of obtaining the analysis spectrum is provided. Characterize.

本願発明によれば、ある程度場所を選ばず簡便かつ迅速にかび毒を分析することができる。また、同一試料に含まれる複数種類のかび毒を並行してそれぞれ分析することができる。   According to the present invention, it is possible to analyze mycotoxin easily and quickly regardless of the location. Further, it is possible to analyze a plurality of types of mycotoxin contained in the same sample in parallel.

図1は実施例の手順を示すフローチャートである。FIG. 1 is a flowchart showing the procedure of the embodiment. 図2は実施例に係る分析装置のブロック図である。FIG. 2 is a block diagram of the analyzer according to the embodiment. 図3は分析用スペクトルの取得の手順をさらに詳しく示したフローチャートである。FIG. 3 is a flowchart showing in more detail the procedure for acquiring the spectrum for analysis. 図4はレーザ光照射とスペクトル記録の処理を説明するタイムチャートである。FIG. 4 is a time chart for explaining the processing of laser light irradiation and spectrum recording. 図5は時系列的に記録されたスペクトルを例示する図である。FIG. 5: is a figure which illustrates the spectrum recorded in time series. 図6は互いに分離されたかび毒(ゼアラレノン)由来のラマンスペクトル及び蛍光スペクトル並びにそれらの強度の時間依存性の分析結果を表す図である。FIG. 6 is a diagram showing a Raman spectrum and a fluorescence spectrum derived from fungal venom (zearalenone) separated from each other and analysis results of time dependence of their intensities. 図7は分析の条件を変えて図6と同じスペクトルが得られた結果を表す図である。FIG. 7 is a diagram showing the results of obtaining the same spectrum as in FIG. 6 by changing the analysis conditions. 図8は演算処理の手順をさらに詳しく示したフローチャートである。FIG. 8 is a flowchart showing the procedure of the arithmetic processing in more detail. 図9は演算処理の結果として抽出されたゼアラレノン由来のラマンスペクトルを表す図である。FIG. 9 is a diagram showing a Raman spectrum derived from zearalenone extracted as a result of the arithmetic processing. 図10はゼアラレノンの標品から得た参照用スペクトルの一例を示す図である。FIG. 10 is a diagram showing an example of a reference spectrum obtained from a preparation of zearalenone. 図11は演算処理の結果として抽出されたアフラトキシン4種由来のラマンスペクトルを表す図である。FIG. 11 is a diagram showing Raman spectra derived from four aflatoxin species extracted as a result of the arithmetic processing. 図12はアフラトキシン4種の標品から得た参照用スペクトルの一例を示す図である。FIG. 12 is a diagram showing an example of a reference spectrum obtained from four aflatoxin preparations. 図13は演算処理の結果として抽出されたデオキシニバレノール由来のラマンスペクトルを表す図である。FIG. 13 is a diagram showing a Raman spectrum derived from deoxynivalenol extracted as a result of the arithmetic processing. 図14はデオキシニバレノールの標品から得た参照用スペクトルの一例を示す図である。FIG. 14 is a diagram showing an example of a reference spectrum obtained from a standard sample of deoxynivalenol.

以下、図面を参照して、本願発明の実施例を説明する。以下に説明する本願発明の実施例は、食品等に目的物質(かび毒)を添加して初期試料を生成し、その初期試料から溶媒へ目的物質を抽出し、抽出物を精製して得た分析用試料から振動分光分析法の一つであるラマン分光分析法を用いて目的物質を検出するという手順で、食品等に添加したかび毒を分析するものである。図1は、実施例の手順を示すフローチャートである。   Embodiments of the present invention will be described below with reference to the drawings. Examples of the present invention described below were obtained by adding a target substance (mold poison) to foods and the like to produce an initial sample, extracting the target substance from the initial sample into a solvent, and purifying the extract. The method of detecting a target substance from a sample for analysis using Raman spectroscopic analysis, which is one of vibration spectroscopic methods, is used to analyze mycotoxin added to foods and the like. FIG. 1 is a flowchart showing the procedure of the embodiment.

上記の手順のうち、試料の生成及び採取から抽出物の精製まで(図1のステップS101及びS102)を、例えば以下のように実施する(本実施例はかび毒の分析方法を実験的に実施するものであるから、かび毒を含む試料の生成を予め行う。)。とうもろこし試料10.0gに4種(B1、B2、G1及びG2)のアフラトキシン混合物を10μg/kgの重量比で、デオキシニバレノールを1mg/kgの重量比で、又はゼアラレノンを1mg/kgの重量比で添加する(ステップS101)。   Of the above procedures, the steps from sample generation and collection to extract purification (steps S101 and S102 in FIG. 1) are carried out, for example, as follows (in this example, the method for analyzing mycotoxin is experimentally carried out). Therefore, the sample containing mycotoxin is generated in advance.). A mixture of four aflatoxins (B1, B2, G1 and G2) in a weight ratio of 10 μg / kg, deoxynivalenol in a weight ratio of 1 mg / kg, or zearalenone in a weight ratio of 1 mg / kg was added to 10.0 g of a corn sample. Add (step S101).

さらに超純水(MilliQ水(同一称呼のMILLI−Qは登録商標、以下同じ。))40mlを加えて振盪し(200rpm、30分)、さらにメタノール60mlを加えて振盪する(200rpm、30分)。これを静置した後、上清50mlを遠心分離して(7000rpm, 10分)、その上清5mlを窒素ガスで2mlに濃縮する。   Furthermore, 40 ml of ultrapure water (MilliQ water (the same name MILLI-Q is a registered trademark, the same applies below)) is added and shaken (200 rpm, 30 minutes), and further 60 ml of methanol is added and shaken (200 rpm, 30 minutes). .. After allowing this to stand, 50 ml of the supernatant is centrifuged (7,000 rpm, 10 minutes), and 5 ml of the supernatant is concentrated to 2 ml with nitrogen gas.

上記の2mlの濃縮液にPBS(リン酸緩衝生理食塩水)5mlを加え、かび毒濃縮カラム(Myco6in1+IACカラム)を用いてカラムトラップする。これを10mlの超純水(MilliQ水)10mlで洗浄し、メタノール1.5mlでかび毒をカラムより溶出させる。さらにカラムからメタノール1.5mlでかび毒を溶出させると共に、水溶性又は不溶性の別が異なるかび毒を溶出するため超純水(MilliQ)2mlを加える。以上のプロセスを経て得られた全溶出液を窒素ガスで乾固し、分光分析の対象とする分析用試料(符号10を付して表す。)を得る(ステップS102)。   5 ml of PBS (phosphate buffered saline) is added to the above 2 ml of concentrated solution, and column trap is performed using a fungal venom concentration column (Myco6in1 + IAC column). This is washed with 10 ml of ultrapure water (MilliQ water), and mycotoxin is eluted from the column with 1.5 ml of methanol. Further, the mold venom is eluted from the column with 1.5 ml of methanol, and 2 ml of ultrapure water (MilliQ) is added to elute the mold venom which is different in water solubility or insolubility. The whole eluate obtained through the above process is dried with nitrogen gas to obtain a sample for analysis (denoted by reference numeral 10) to be subjected to spectroscopic analysis (step S102).

ここでいったん図1から離れて、実施例に係るかび毒の分析装置100の構成を説明する。図2は、分析装置100のブロック図である。図2の左側に、ステップS102によって得られた分析用試料10を示す。分析装置100は、図1において一点鎖線で囲まれた構成を備える。それらの構成は、分光分析部110、演算処理部120及び制御部130である。分光分析部110は、受光部111、分光部112、記録部113及び分離部114を備える。演算処理部120は、生成部121及び解析部122を備える。   Here, the configuration of the mold venom analyzer 100 according to the embodiment will be described once with reference to FIG. FIG. 2 is a block diagram of the analyzer 100. The analysis sample 10 obtained in step S102 is shown on the left side of FIG. The analyzer 100 has a configuration surrounded by a chain line in FIG. Those components are the spectroscopic analysis unit 110, the arithmetic processing unit 120, and the control unit 130. The spectroscopic analysis section 110 includes a light receiving section 111, a spectroscopic section 112, a recording section 113, and a separating section 114. The arithmetic processing unit 120 includes a generation unit 121 and an analysis unit 122.

実施例では、レーザ光源を含む発光部20も使用される。分析装置100が発光部20を含む構成であってもよい。発光部20及び受光部111は、各種の商品化されたラマン分光分析装置の光学系に相当するように構成されることができる(ここでは詳しい説明を省略する。)。   In the embodiment, the light emitting unit 20 including the laser light source is also used. The analysis device 100 may include the light emitting unit 20. The light emitting unit 20 and the light receiving unit 111 can be configured to correspond to optical systems of various commercialized Raman spectroscopic analyzers (detailed description thereof will be omitted here).

上記の分析装置100の構成は、ハードウェア又はソフトウェアを含む多様な手段を用いて実現され得るものであって、必ずしも上記の各構成のとおりに明りょうに区画されたハードウェア又はソフトウェアのブロックが存在することを限定的に意味するものではない。図1において、矢印付きの実線は計測及び処理される信号の流れを示し、矢印付きの破線は制御信号の流れを示す。   The above-described configuration of the analysis device 100 can be realized by using various means including hardware or software, and a block of hardware or software clearly divided as in each of the above configurations is not necessarily required. It does not imply that it is present. In FIG. 1, solid lines with arrows show the flow of signals to be measured and processed, and broken lines with arrows show the flow of control signals.

発光部20は、分析用試料10等の対象物に対してレーザ光を照射することができる。受光部111は、レーザ照射の結果として上記の対象物から散射されるラマン散乱光及び蛍光その他の背景光を集光する。分光部112は、受光部111が集光した光を分析して、ラマンシフトを一方の軸(例えば横軸)及び光の強度を他方の軸(例えば縦軸)とする平面上にそのスペクトルをプロットすることができる。ラマンシフトを表す軸は、光の波長及び振動数のいずれかを単位とする軸に替えてもよい。記録部113は、上記のプロットされたスペクトルのデータを記録することができる。   The light emitting unit 20 can irradiate an object such as the analysis sample 10 with laser light. The light receiving unit 111 collects Raman scattered light, fluorescence and other background light scattered from the object as a result of laser irradiation. The spectroscopic unit 112 analyzes the light condensed by the light receiving unit 111 and displays the spectrum on a plane having Raman shift on one axis (for example, the horizontal axis) and light intensity on the other axis (for example, the vertical axis). Can be plotted. The Raman shift axis may be replaced with an axis whose unit is either the wavelength of light or the frequency of light. The recording unit 113 can record the data of the plotted spectrum.

分離部114は、例えばパーソナルコンピュータ(PC)に組み込まれたソフトウェアによって実現され、記録部113に記録されたスペクトルのデータに対して後述する多変量解析を含む分離処理を行うことにより、ラマン散乱光若しくは背景光の信号強度の時系列並びに信号対雑音比(SNR)が改善されたラマン散乱光を含む分析用スペクトルを得ることができる。   The separation unit 114 is realized by, for example, software installed in a personal computer (PC), and performs separation processing including multivariate analysis described below on the spectrum data recorded in the recording unit 113 to generate Raman scattered light. Alternatively, it is possible to obtain a spectrum for analysis including Raman scattered light in which the time series of the signal intensity of the background light and the signal-to-noise ratio (SNR) are improved.

演算処理部120の生成部121及び解析部122は、例えばPCに組み込まれたソフトウェアによって実現される。生成部121は上記の分析用スペクトル及び後述する参照用スペクトルを得て、行列として表すことができる解析対象データを生成することができる。解析部122は上記の解析対象データを得て、後述する多変量解析を含む演算処理を行うことにより、分析用試料10に含まれるかび毒のスペクトル及び信号強度プロファイル(一般に、検出されたかび毒の量に対して線形の関係にある。)を得ることができる。   The generation unit 121 and the analysis unit 122 of the arithmetic processing unit 120 are realized by software installed in a PC, for example. The generation unit 121 can obtain the analysis spectrum described above and the reference spectrum described later, and generate analysis target data that can be represented as a matrix. The analysis unit 122 obtains the above-mentioned data to be analyzed and performs a calculation process including a multivariate analysis described later to obtain a spectrum and a signal intensity profile (generally, detected mycotoxin of the mycotoxin contained in the analysis sample 10). , Which has a linear relationship with the quantity of.

制御部130は、例えばPCに組み込まれたソフトウェアによって実現される。制御部130は、上述した発光部20、分光分析部110及び演算処理部120に対して、データの取得と解析を行うのに必要な指示を与えると共に結果を受け取ることができる。   The control unit 130 is realized by software installed in a PC, for example. The control unit 130 can give an instruction necessary for performing data acquisition and analysis to the light emitting unit 20, the spectroscopic analysis unit 110, and the arithmetic processing unit 120, and can receive the result.

分析装置100の具体的な構成は、上述したものに限らない。それぞれの部分をハードウェア又はファームウェア主体に構成してもよく、ハードウェア又はファームウェアとソフトウェアを適宜組み合わせてもよく、本発明に係るプログラムを商品化された分光分析装置にインストールしたものであってもよい。   The specific configuration of the analysis device 100 is not limited to the above. Each part may be configured mainly by hardware or firmware, hardware or firmware and software may be appropriately combined, or the program according to the present invention may be installed in a commercialized spectroscopic analyzer. Good.

図1に戻り、実施例の手順をステップS103から説明する。このステップでは、図2において分析用試料10に代えて分析の目的とする種類のかび毒を含む標準試料(例えば、当該かび毒の標品自体又は標品を既知の溶媒に混合したもの)に対して発光部20からレーザ光を照射し、当該標準試料からのラマン散乱光を受光部111で受光する。受光されたラマン散乱光は分光部112で分光され、参照用スペクトルとして記録部113に記録される(ステップS103)。   Returning to FIG. 1, the procedure of the embodiment will be described from step S103. In this step, in place of the analytical sample 10 in FIG. 2, a standard sample containing a mold venom of the type to be analyzed (for example, the sample of the mold venom itself or a mixture of the sample and a known solvent) is used. On the other hand, laser light is emitted from the light emitting unit 20, and Raman scattered light from the standard sample is received by the light receiving unit 111. The received Raman scattered light is dispersed by the spectroscopic unit 112 and recorded in the recording unit 113 as a reference spectrum (step S103).

なお、標準試料が溶媒を用いて調製されている場合は、標準試料から得たスペクトルと溶媒のみから得たスペクトルの差分を参照用スペクトルとして、記録部113に記録してもよい。かび毒の標準試料から得た参照用スペクトルの例を、後述する図10等に示している。   When the standard sample is prepared using a solvent, the difference between the spectrum obtained from the standard sample and the spectrum obtained only from the solvent may be recorded in the recording unit 113 as a reference spectrum. An example of the reference spectrum obtained from the standard sample of mycotoxin is shown in FIG.

次にステップS104に移り、まず発光部20から分析用試料10に対してレーザ光を照射し、分析の目的とするかび毒に由来するラマン散乱光、及び当該かび毒以外の成分(夾雑物という。)に由来する蛍光をはじめとする背景光を受光部111で受光し、さらに分光部112で分光して得たスペクトルを、記録部113に記録する。このステップでは、かび毒由来のラマンスペクトルは分析用試料10に含まれる夾雑物由来の蛍光等の背景光のスペクトルに埋もれて、検出することが難しいという問題がある。   Next, the process proceeds to step S104. First, the light emitting unit 20 irradiates the analysis sample 10 with laser light, and Raman scattered light derived from the mycotoxin to be analyzed and components other than the mycotoxin (referred to as impurities). Background light such as fluorescence derived from (..) is received by the light receiving unit 111, and further, the spectrum obtained by being separated by the spectroscopic unit 112 is recorded in the recording unit 113. In this step, there is a problem that the Raman spectrum derived from mycotoxin is buried in the spectrum of background light such as fluorescence derived from contaminants contained in the analysis sample 10 and is difficult to detect.

図3は、ステップS104における処理の手順をさらに詳しく示したフローチャートである。試料中の夾雑物に由来する蛍光の強度は、当該試料に対するレーザ光の照射を続けると時間と共に減衰することが知られている(この現象を消光という。)。図3のフローチャートに示した処理においては、以下のように蛍光の消光現象を利用する。   FIG. 3 is a flowchart showing the processing procedure in step S104 in more detail. It is known that the intensity of fluorescence derived from impurities in a sample is attenuated with time when the sample is continuously irradiated with laser light (this phenomenon is called quenching). In the process shown in the flowchart of FIG. 3, the fluorescence quenching phenomenon is used as follows.

まず、発光部20から分析用試料10に対して、ある時間にわたり連続してレーザ光を照射する(ステップS301)。そしてレーザ光の照射時間を区切った時間インターバルごとに、分析用試料10から得られるスペクトル(ラマン散乱光及び背景光のスペクトルを含む。)を記録部113に記録する(ステップS302)。   First, the light emitting unit 20 continuously irradiates the analysis sample 10 with laser light for a certain period of time (step S301). Then, the spectrum (including the spectrum of the Raman scattered light and the background light) obtained from the analysis sample 10 is recorded in the recording unit 113 at each time interval that divides the irradiation time of the laser light (step S302).

制御部130は、上記のレーザ光の照射時間を連続する複数の時間インターバルに区切ることができる。記録部113は、時系列に従って(すなわち、制御部130から指示された複数の時間インターバルの順で)、分光部112がプロットしたスペクトルを記録することができる。   The control unit 130 can divide the laser light irradiation time into a plurality of continuous time intervals. The recording unit 113 can record the spectrum plotted by the spectroscopic unit 112 in chronological order (that is, in the order of a plurality of time intervals instructed by the control unit 130).

図4は、ステップS301及びS302の処理を説明するタイムチャートである。図の横軸は、時間軸である。時間軸上の”0”点において、発光部20から分析用試料10に対してレーザ光の照射を開始する。制御部130は、時間軸上のT、2T、…、NT(Tは時間、Nは自然数)の各タイミングにおいて、分析用試料10から得られるスペクトルを記録部113に記録させる。   FIG. 4 is a time chart explaining the processing of steps S301 and S302. The horizontal axis of the figure is the time axis. At the “0” point on the time axis, irradiation of laser light from the light emitting unit 20 to the analysis sample 10 is started. The control unit 130 causes the recording unit 113 to record the spectrum obtained from the analysis sample 10 at each timing of T, 2T, ..., NT (T is time, N is a natural number) on the time axis.

図5は、上述のように時系列的に記録されたスペクトルを例示する図である(分析用試料10は、トウモロコシにゼアラレノンを添加したサンプルのIMCカラム抽出定容液を固化したものである。励起光の波長は532nm、T=5秒、N=20とする。)。図5における各プロットは、時間の経過とともに上から下へ移行する。つまり、このスペクトル中の支配的な成分である夾雑物由来の蛍光の強度が時間の経過と共に低下していることがわかる。   FIG. 5: is a figure which illustrates the spectrum recorded in time series as mentioned above (The sample 10 for analysis is what solidified the IMC column extraction constant volume liquid of the sample which added zearalenone to corn. The wavelength of the excitation light is 532 nm, T = 5 seconds, and N = 20.) Each plot in FIG. 5 shifts from top to bottom with the passage of time. That is, it can be seen that the intensity of the fluorescence derived from the contaminant, which is the dominant component in this spectrum, decreases with the passage of time.

すべての時間インターバルについて記録されたスペクトルのデータ(以下、全スペクトルデータという。)は、例えば行の番号をラマンシフト軸上の位置(ピクセルという。)により、列の番号を上記の時系列上の順序により、それぞれ定めた行列の形で表すことができる。ピクセル数をMとすれば、全スペクトルデータはM行N列の行列(以下、符号Aで表す。)として表される。   The spectrum data recorded for all time intervals (hereinafter, referred to as all spectrum data) includes, for example, the row number by the position (pixel) on the Raman shift axis, and the column number by the above time series. It can be expressed in the form of a matrix defined by the order. Assuming that the number of pixels is M, all spectrum data are represented as a matrix of M rows and N columns (hereinafter referred to as symbol A).

次に、分離部114が行う全スペクトルデータの処理(分離処理)を説明する。例えば分析用試料10を構成する成分が、かび毒(ラマンスペクトルに寄与。)と夾雑物(蛍光スペクトルに寄与。)の2種類に分けられることが既知であるか又は推測される場合は、行列Aを、上記のそれぞれの成分のスペクトルデータを表す列ベクトルからなるM行2列の行列Wと、それぞれの成分の強度の時系列データを表す行ベクトルからなる2行N列の行列Hに因子分解することを考える。そうすると、行列Aは「数1」式、行列W及びHは「数2」式でそれぞれあらわされる。なお、「数1」式中の「実測スペクトル(kT)」は、図4の時間軸上のタイミング「kT」(1≦k≦N)において実測され記録されたスペクトルを表す。   Next, processing (separation processing) of all spectrum data performed by the separation unit 114 will be described. For example, when it is known or presumed that the components constituting the analytical sample 10 can be divided into two types, that is, mycotoxin (contribution to Raman spectrum) and contaminant (contribution to fluorescence spectrum), a matrix is formed. Factor A into a matrix W of M rows and 2 columns consisting of column vectors representing the spectrum data of each component and a matrix H of 2 rows N columns consisting of row vectors representing time series data of the intensity of each component. Think to disassemble. Then, the matrix A is expressed by the "equation 1" and the matrices W and H are expressed by the "equation 2", respectively. The "actually measured spectrum (kT)" in the equation 1 represents the spectrum actually measured and recorded at the timing "kT" (1≤k≤N) on the time axis of FIG.

Figure 0006692226
Figure 0006692226

Figure 0006692226
Figure 0006692226

上記の「数1」及び「数2」式から行列Wを求めることができれば、夾雑物由来の蛍光スペクトルから分離されると共に、時系列に沿って積分されることによりガウス雑音の影響が軽減されたかび毒由来のラマンスペクトルが得られる。「数1」式の解析的な解は得られないが、行列Wに初期値を与えて公知の最小二乗法により行列Hの近似解を求め、これからさらにWの近似解を求めるという演算を所定の収束条件が満たされるまで反復する交互最小二乗法により、近似解を得られることが知られている。   If the matrix W can be obtained from the above equations 1 and 2, the influence of Gaussian noise can be reduced by separating from the fluorescence spectrum derived from contaminants and integrating them in time series. Raman spectra derived from mycotoxin are obtained. Although an analytical solution of the equation (1) cannot be obtained, an initial value is given to the matrix W to obtain an approximate solution of the matrix H by a known least squares method, and then an approximate solution of W is further obtained. It is known that an approximate solution can be obtained by the alternating least squares method that iterates until the convergence condition of is satisfied.

上述した交互最小二乗法では、実測値を要素とする行列A及び近似値を要素とする行列W・Hの残差の二乗ノルム||A−W・H||を行列W又はHのいずれか一方について最小化することにより、行列W又はHのもう一方の次の段階の近似値を求める。この最小化演算をよりシンプルで高精度なものにするための数学的な工夫として、少なくとも行列Aの全要素が非負値で表される条件と、上記の残差の二乗ノルムに正則化項(例えばL1ノルム又はL2ノルム)を加算する拘束条件を課すことが知られている。また、ラマンスペクトルの強度の時系列データに時間依存性が顕著でないことから、行列Hの第1行の値を定数とする条件を課すことができる。 In the alternating least squares method described above, the squared norm of residuals || A−W · H || 2 of the matrix A whose actual measurement value is an element and the matrix W · H whose approximate value is an element are either the matrix W or H. One of the matrices W or H is approximated to the next stage by minimizing the other. As a mathematical device for making this minimization operation simpler and more accurate, at least all the elements of the matrix A are represented by non-negative values and the regularization term ( For example, it is known to impose a constraint condition of adding the L1 norm or the L2 norm). Moreover, since the time-series data of the intensity of the Raman spectrum does not have a remarkable time dependency, it is possible to impose a condition that the value in the first row of the matrix H is a constant.

図6は、図5に示した全スペクトルデータから上記の分離処理を経て求めた、「数2」式の行列Wの各列ベクトル及び行列Hの各行ベクトルの近似解を示す図である。図6の左側の図は、ラマンシフトを横軸に、夾雑物由来の蛍光スペクトル(番号1)及びかび毒由来のラマンスペクトル(番号2)を表している。図6の右側の図は、時系列上の測定番号(1≦N≦20)を横軸に、蛍光強度の時間依存性(番号1、実践)を表している(ラマンスペクトルの時間依存性(番号2、破線)はゼロと仮定した。)(ステップS303)。   FIG. 6 is a diagram showing an approximate solution of each column vector of the matrix W and each row vector of the matrix H of the “Equation 2” obtained by performing the above separation processing from all the spectrum data shown in FIG. The diagram on the left side of FIG. 6 shows the fluorescence spectrum (No. 1) derived from contaminants and the Raman spectrum (No. 2) derived from mycotoxin with the Raman shift as the horizontal axis. In the diagram on the right side of FIG. 6, the measurement number (1 ≦ N ≦ 20) on the time series is plotted on the horizontal axis, and the time dependence of fluorescence intensity (No. 1, practice) is shown (time dependence of Raman spectrum ( It is assumed that the number 2 and the broken line) are zero) (step S303).

分析用試料10に含まれる成分が例えば4種類に分けられることが既知であるか又は推測される場合は、行列Aを、それぞれの成分のスペクトルデータを表す列ベクトルからなるM行4列の行列Wと、それぞれの成分の強度の時系列データを表す行ベクトルからなる4行N列の行列Hに因子分解することを考える。そうすると、行列W及びHは「数3」式で表される。   When it is known or inferred that the components contained in the analytical sample 10 can be divided into, for example, four types, the matrix A is an M-row, 4-column matrix composed of column vectors representing the spectrum data of the respective components. Consider factorization into W and a matrix H of 4 rows and N columns consisting of row vectors representing time series data of intensity of each component. Then, the matrices W and H are represented by the formula "Equation 3".

Figure 0006692226
Figure 0006692226

「数1」及び「数3」式から、上記と同様の交互最小二乗法を用いて、行列W及びHの近似解を求めることができる。この場合にも、行列Aの全要素が非負値で表される条件と正則化項の拘束条件を用いる。また、行列W及びHの全要素が非負値で表される条件を課すことができる。   An approximate solution of the matrices W and H can be obtained from the equations 1 and 3 by using the alternating least squares method similar to the above. Also in this case, the condition that all the elements of the matrix A are represented by non-negative values and the constraint condition of the regularization term are used. Further, it is possible to impose a condition that all elements of the matrices W and H are represented by non-negative values.

図7は、図5に示した全スペクトルデータから上記の分離処理を経て求めた「数3」式の行列Wの各列ベクトル及び行列Hの各行ベクトルの近似解を示す図である。図7の左側の図は、ラマンシフトを横軸に、各成分(番号1ないし4)のスペクトルを表している。図7の右側の図は、時系列上の測定番号(1≦N≦20)を横軸に、各成分(番号1ないし4)の時間依存性を表している。番号1及び番号2の成分は、それぞれ図6の蛍光スペクトル(番号1)及びラマンスペクトル(番号2)とよく一致していることがわかる。このように分析用試料10に含まれる成分の種類が増えても、かび毒由来のラマンスペクトルを夾雑物由来の背景光から分離することができる(ステップS303)。   FIG. 7 is a diagram showing an approximate solution of each column vector of the matrix W and each row vector of the matrix H of the “Formula 3” obtained by performing the above separation processing from all the spectrum data shown in FIG. The diagram on the left side of FIG. 7 shows the spectrum of each component (numbers 1 to 4) with the Raman shift as the horizontal axis. The diagram on the right side of FIG. 7 shows the time dependence of each component (numbers 1 to 4) with the measurement number (1 ≦ N ≦ 20) on the time series as the horizontal axis. It can be seen that the components of No. 1 and No. 2 are in good agreement with the fluorescence spectrum (No. 1) and Raman spectrum (No. 2) of FIG. 6, respectively. Thus, even if the types of components contained in the analysis sample 10 increase, the Raman spectrum derived from mycotoxin can be separated from the background light derived from impurities (step S303).

以上の通り、分析用試料10に含まれる夾雑物由来の蛍光強度の時間依存性に着目して、かび毒由来のラマンスペクトルを分離して分析用スペクトルとして得ることができた(ステップS304、ステップS104)。続いて、ステップS103で得た参照用スペクトルとステップS104で得た分析用スペクトルを入力として、多変量解析を用いた演算処理を行うことにより、かび毒由来のラマンスペクトルを抽出する(ステップS105)。図8は、ステップS105における処理の手順をさらに詳しく示したフローチャートである。   As described above, paying attention to the time dependence of the fluorescence intensity derived from the contaminants contained in the analysis sample 10, the Raman spectrum derived from the mold must be separated and obtained as the analysis spectrum (step S304, step). S104). Subsequently, the Raman spectrum derived from mycotoxin is extracted by performing arithmetic processing using multivariate analysis with the reference spectrum obtained in step S103 and the analysis spectrum obtained in step S104 as inputs (step S105). .. FIG. 8 is a flowchart showing the procedure of the process in step S105 in more detail.

生成部121は、ステップS103で得られた参照用スペクトルを記録部113から読み出して複数通りの係数(実数)を乗算し、複数通りの積データを得る(ステップS801)。生成部121は、次にステップS104(又はS304)で得られた分析用スペクトルを分離部114から得て、上記の複数通りの積データとそれぞれ加算することにより、複数通りの和スペクトルのデータを得る(ステップS802)。上記の複数通りの和スペクトルのデータは、実空間において検量線を引くため分析対象試料に複数通りの濃度の標準試料を添加し、その都度測定されるスペクトルのデータに等価である。   The generation unit 121 reads the reference spectrum obtained in step S103 from the recording unit 113 and multiplies it by a plurality of coefficients (real numbers) to obtain a plurality of product data (step S801). Next, the generation unit 121 obtains the analysis spectrum obtained in step S104 (or S304) from the separation unit 114, and adds the analysis spectrum to each of the above-described plurality of product data to obtain a plurality of sum spectrum data. Get (step S802). The above-mentioned data of a plurality of sum spectra is equivalent to the data of a spectrum measured each time when a standard sample having a plurality of concentrations is added to a sample to be analyzed in order to draw a calibration curve in real space.

生成部121は、上記の複数通りの和スペクトルのデータに対して、ラマンシフト軸(又は対応する波長軸若しくは振動数軸)の向きにそれぞれ0回以上微分して得られたデータを列ベクトルとする行列を構成し、当該行列を解析対象データとして生成する(ステップS803)。解析対象データの行列は、ステップS302で得られた全スペクトルデータの行列と同様に、かび毒及びかび毒以外のそれぞれの成分のスペクトルデータを表す2列の列ベクトルからなる行列と、参照用スペクトルに乗算した係数(実空間における標準試料の濃度に相当する。)に対応するそれぞれの成分の強度プロファイルを表す2行の行ベクトルからなる行列の積に因子分解することができる。   The generation unit 121 differentiates the data of the above-described plurality of sum spectra in the direction of the Raman shift axis (or the corresponding wavelength axis or frequency axis) 0 times or more, and obtains the obtained data as a column vector. To form a matrix and generate the matrix as analysis target data (step S803). The matrix of the data to be analyzed is, like the matrix of all spectrum data obtained in step S302, a matrix composed of two column vectors representing spectrum data of each component other than mycotoxin and a spectrum for reference. Can be factored into a product of a matrix composed of two row vectors representing the intensity profile of each component corresponding to the coefficient (corresponding to the concentration of the standard sample in the real space) multiplied by.

解析部122は上記の解析対象データを得て、例えばかび毒以外の成分の強度プロファイルが一定である条件及び正則化項を加算する拘束条件を課した交互最小二乗法を用いることにより、かび毒由来のラマンスペクトル及び信号強度プロファイルを抽出することができる(特許第5780476号公報参照。)。   The analysis unit 122 obtains the above-described data to be analyzed and, for example, by using the alternate least squares method in which the intensity profile of components other than the mycotoxin is constant and the constraint condition for adding the regularization term is imposed, the mycotoxin is It is possible to extract the Raman spectrum and the signal intensity profile derived from it (see Japanese Patent No. 5780476).

図9は、ステップS104(S304)で得られたゼアラレノン由来の分析用スペクトル及びゼアラレノン標品から得た参照用スペクトルを用いて、ステップS801ないしS804の処理によって抽出されたゼアラレノン由来のラマンスペクトルを示す図である。図中の実線(最上段のプロット)は、蛍光から分離された後の分析用スペクトルである。図中の点線(最上段のすぐ下のプロット)は、ゼアラレノン以外の成分由来のスペクトルである。図中の一点鎖線(最下段のプロット)は、ゼアラレノン由来のラマンスペクトルである。   FIG. 9 shows a Raman spectrum derived from zearalenone extracted by the processes of steps S801 to S804 using the analysis spectrum derived from zearalenone obtained in step S104 (S304) and the reference spectrum obtained from the zearalenone preparation. It is a figure. The solid line (uppermost plot) in the figure is the analysis spectrum after being separated from the fluorescence. The dotted line in the figure (the plot immediately below the top row) is the spectrum derived from components other than zearalenone. The alternate long and short dash line (bottom plot) in the figure is the Raman spectrum derived from zearalenone.

図10は、ゼアラレノンの標品から得た参照用スペクトルの一例を示す。図9の最下段のプロットが図10のプロットと共通の特徴を示すことから、本実施例によりゼアラレノンのラマンスペクトルを抽出することができたことが確認される。   FIG. 10 shows an example of a reference spectrum obtained from a preparation of zearalenone. Since the plot at the bottom of FIG. 9 has the same characteristics as the plot in FIG. 10, it is confirmed that the Raman spectrum of zearalenone could be extracted by this example.

図11は、かび毒としてアフラトキシン4種(G1、G2、B1及びB2)を含む試料から抽出されたそれらのカビ毒由来のラマンスペクトルを、図9と同様に示す図である。図中の実線(最上段のプロット)は、蛍光から分離された後の分析用スペクトルである。図中の点線(最上段のすぐ下のプロット)は、アフラトキシン4種以外の成分由来のスペクトルである。図中の一点鎖線(最下段のプロット)4本は、アフラトキシン4種由来のラマンスペクトルである。図12のアフラトキシン4種の標品の参照用スペクトルと対比することにより、アフラトキシン4種のラマンスペクトルを抽出することができたことが確認される。これは、同一試料に含まれる複数種のかび毒を並行してそれぞれ分析し抽出した例である。   FIG. 11 is a diagram showing Raman spectra derived from those mold venoms extracted from a sample containing four types of aflatoxins (G1, G2, B1 and B2) as a mold toxicant, similar to FIG. 9. The solid line (uppermost plot) in the figure is the analysis spectrum after being separated from the fluorescence. The dotted line in the figure (the plot immediately below the top row) is the spectrum derived from components other than the four aflatoxins. Four dashed-dotted lines (lowermost plot) in the figure are Raman spectra derived from four aflatoxins. It is confirmed that the Raman spectra of the four types of aflatoxins could be extracted by comparing with the reference spectra of the four types of aflatoxin preparations in FIG. This is an example in which multiple types of mycotoxin contained in the same sample were analyzed and extracted in parallel.

図13は、かび毒としてデオキシニバレノールを含む試料から抽出されたデオキシニバレノール由来のラマンスペクトルを、図9と同様に示す図である。図中の実線(最上段のプロット)は、蛍光から分離された後の分析用スペクトルである。図中の点線(最上段のすぐ下のプロット)は、デオキシニバレノール以外の成分由来のスペクトルである。図中の一点鎖線(最下段のプロット)は、デオキシニバレノール由来のラマンスペクトルである。図14のデオキシニバレノールの標品の参照用スペクトルと対比することにより、デオキシニバレノールのラマンスペクトルを抽出することができたことが確認される。   FIG. 13: is a figure which shows the Raman spectrum derived from deoxynivalenol extracted from the sample which contains deoxynivalenol as a mold poison similarly to FIG. The solid line (uppermost plot) in the figure is the analysis spectrum after being separated from the fluorescence. The dotted line in the figure (the plot immediately below the top row) is the spectrum derived from components other than deoxynivalenol. The alternate long and short dash line (bottom plot) in the figure is a Raman spectrum derived from deoxynivalenol. It is confirmed that the Raman spectrum of deoxynivalenol could be extracted by comparing with the reference spectrum of the standard sample of deoxynivalenol in FIG.

以上に説明した本発明の実施例によれば、試料中の夾雑物由来の蛍光等の背景光の影響を抑えて、かび毒由来のラマンスペクトルを抽出できることが確認された。実施例では、分析用スペクトルの取得(ステップS104)のためにステップS301ないしS304の処理を実行したが、試料の種類や特性によってはこれらの処理を省いて分析用スペクトルを取得してもよい。実施例ではラマン分光分析法を用いたが、図1のフローチャートに示した本発明の方法は、赤外吸収分光分析法を用いても実施することができる。   According to the examples of the present invention described above, it was confirmed that the Raman spectrum derived from mycotoxin can be extracted while suppressing the influence of background light such as fluorescence derived from contaminants in the sample. In the embodiment, the processes of steps S301 to S304 are executed to acquire the analysis spectrum (step S104), but depending on the type and characteristics of the sample, these processes may be omitted to acquire the analysis spectrum. Although Raman spectroscopy was used in the examples, the method of the present invention shown in the flowchart of FIG. 1 can also be carried out using infrared absorption spectroscopy.

実施例では図1のステップS102において抽出物を精製することを説明したが、かび毒を含む初期試料の特性によっては抽出物の精製を省き、初期試料をそのまま分析用試料にしてもよい。   Although the example explained that the extract is purified in step S102 in FIG. 1, the extract may be omitted depending on the characteristics of the initial sample containing the mycotoxin, and the initial sample may be directly used as the analytical sample.

10 分析用試料
20 発光部
100 分析装置
110 分光分析部
111 受光部
112 分光部
113 記録部
114 分離部
120 演算処理部
121 生成部
122 解析部
130 制御部
10 Analysis Sample 20 Light Emitting Part 100 Analyzing Device 110 Spectroscopic Analyzing Part 111 Light Receiving Part 112 Spectroscopic Part 113 Recording Part 114 Separing Part 120 Operation Processing Part 121 Generating Part 122 Analyzing Part 130 Control Part

Claims (8)

ラマン分光分析法を用いるかび毒の分析方法において、
分析の目的とする種類のかび毒が産生する可能性がある初期試料を採取し、
前記採取した初期試料から、前記分析の目的とする種類のかび毒を含む可能性がある分析用試料を生成し、
ラマン分光分析法を用いて、前記分析の目的とする種類のかび毒を含む標準試料から参照用スペクトルを取得し、
前記分析用試料に対して、連続する複数の時間インターバルにわたってレーザ光を照射し、
前記複数の時間インターバルごとに、前記分析用試料から散射されるラマン散乱光及び前記ラマン散乱光を除く背景光を集光すると共に、前記ラマン散乱光及び前記背景光を含むスペクトルを記録し、
前記複数の時間インターバルごとに記録されたスペクトルから、多変量解析を用いて、前記ラマン散乱光の強度及び前記背景光の強度のそれぞれの時間依存性を抽出すると共に前記ラマン散乱光のスペクトルを前記背景光のスペクトルから分離する分離処理を行うことにより、前記分析用試料から分析用スペクトルを取得し、
前記参照用スペクトルのデータ及び前記分析用スペクトルのデータを入力とする多変量解析を用いた演算処理を行うことにより、前記分析の目的とする種類のかび毒が前記初期試料に産生する場合において、前記分析用スペクトルから前記分析の目的とする種類のかび毒のスペクトルを抽出する
ことを特徴とするかび毒の分析方法。
In the method of analyzing mycotoxin using Raman spectroscopy,
Collecting an initial sample that may produce the type of mycotoxin for the purpose of analysis,
From the collected initial sample, generate an analytical sample that may contain a mold venom of the type intended for the analysis,
Using Raman spectroscopy, obtain a reference spectrum from a standard sample containing the type of mycotoxin of interest for the analysis,
Irradiating the analysis sample with laser light over a plurality of consecutive time intervals,
For each of the plurality of time intervals, while condensing the Raman scattered light scattered from the analysis sample and the background light excluding the Raman scattered light, record a spectrum containing the Raman scattered light and the background light,
From the spectra recorded for each of the plurality of time intervals, using multivariate analysis, the time dependence of the intensity of the Raman scattered light and the intensity of the background light and the spectrum of the Raman scattered light is extracted. By performing a separation process to separate from the spectrum of background light, to obtain an analysis spectrum from the analysis sample,
By performing arithmetic processing using the multivariate analysis with the data of the reference spectrum and the data of the analysis spectrum as an input, in the case where the type of mycotoxin of the type of the analysis is produced in the initial sample, A method for analyzing a mycotoxin, which comprises extracting a spectrum of a mycotoxin of a type intended for the analysis from the spectrum for analysis.
前記演算処理は、
前記参照用スペクトルのデータに複数通りの係数を乗算して得られた複数通りの積のデータのそれぞれを前記分析用スペクトルのデータに加算して得られた複数通りの和スペクトルのデータから、行列で表すことができる解析対象データを生成し、
多変量解析を用いて、前記解析対象データから前記分析の目的とする種類のかび毒のスペクトル及び信号強度プロファイルを抽出する
ことを特徴とする請求項1に記載のかび毒の分析方法。
The arithmetic processing is
From a plurality of sum spectrum data obtained by adding each of a plurality of product data obtained by multiplying the reference spectrum data by a plurality of coefficients, a matrix Generate analysis target data that can be expressed by
The method for analyzing mycotoxin according to claim 1, wherein a spectrum and signal intensity profile of the type of mycotoxin targeted for the analysis are extracted from the data to be analyzed using multivariate analysis.
前記分析の目的とする種類のかび毒に複数の種類がある場合において、
前記複数の種類のかび毒ごとに前記参照用スペクトルを取得し、
前記複数の分析の目的とする種類ごとに前記演算処理を行う
ことを特徴とする請求項1に記載のかび毒の分析方法。
In the case where there are multiple types of mycotoxin of the type of the analysis,
Obtaining the reference spectrum for each of the plurality of types of mycotoxin,
The method for analyzing a mycotoxin according to claim 1, wherein the arithmetic processing is performed for each target type of the plurality of analyses.
前記分離処理における多変量解析は、
前記記録されたスペクトルを、前記ラマン散乱光のスペクトル及び前記背景光のスペクトルをそれぞれ列ベクトルとして含む第1の行列と、前記複数の時間インターバルごとの順序に対応する前記ラマン散乱光の信号強度プロファイル及び前記複数の時間インターバルごとの順序に対応する前記背景光の信号強度プロファイルをそれぞれ行ベクトルとして含む第2の行列に分解して表したとき、
少なくとも前記記録されたスペクトルが非負値行列で表されるという条件と、誤差に正則化項を加算する拘束条件とを課した最小二乗法を用いて、前記第1の行列及び前記第2の行列を交互に求める演算を所定の収束条件が満たされるまで反復する
ことを特徴とする請求項1に記載のかび毒の分析方法。
Multivariate analysis in the separation process,
A first matrix containing the recorded spectrum as a column vector of the spectrum of the Raman scattered light and the spectrum of the background light as a column vector, and a signal intensity profile of the Raman scattered light corresponding to the order of the plurality of time intervals. And when decomposed into a second matrix including the signal intensity profile of the background light corresponding to the order of each of the plurality of time intervals as a row vector,
Using the least squares method that imposes at least the condition that the recorded spectrum is represented by a non-negative matrix and a constraint condition that adds a regularization term to an error, the first matrix and the second matrix The method for analyzing mycotoxin according to claim 1, wherein the calculation for alternately obtaining is repeated until a predetermined convergence condition is satisfied.
前記初期試料に対して前記分析の目的とする種類のかび毒に適合した溶媒を用いた抽出操作及び精製操作を行うことにより、前記分析用試料を生成することを特徴とする請求項1に記載のかび毒の分析方法。   The sample for analysis is produced by subjecting the initial sample to an extraction operation and a purification operation using a solvent suitable for the type of mycotoxin targeted for the analysis. Method for analysis of mycotoxin. 前記かび毒は、アフラトキシン4種(B1、B2、G1及びG2)、デオキシニバレノール並びにゼアラレノンのうち1種類以上を含むことを特徴とする請求項1ないし請求項5のいずれか1項に記載のかび毒の分析方法。   The mold according to any one of claims 1 to 5, wherein the mold venom contains one or more of four types of aflatoxins (B1, B2, G1 and G2), deoxynivalenol and zearalenone. Poison analysis method. 振動分光分析法を用いてかび毒の分析を行うことができるかび毒の分析装置において、
分析の目的とする種類のかび毒が産生する可能性がある初期試料から生成された分析用試料から、前記振動分光分析法を用いて分析用スペクトルを取得することができる分光分析部と、
前記分析用スペクトルのデータ及び前記分析の目的とする種類のかび毒を含む標準試料から前記振動分光分析法を用いて取得された参照用スペクトルのデータを入力として、多変量解析を用いた演算処理を行うことにより、前記分析の目的とする種類のかび毒が前記試料に産生する場合において前記分析の目的とする種類のかび毒のスペクトルを抽出する
ことができる演算処理部とを備え
前記分光分析部は、
レーザ光を照射された前記分析用試料が散射するラマン散乱光及び前記ラマン散乱光を除く背景光を集光して、前記ラマン散乱光及び前記背景光を含むスペクトルを記録することができる記録部と、
連続する複数の時間インターバルにわたって前記分析用試料が連続してレーザ光を照射されたときに、前記複数の時間インターバルごとに記録されたスペクトルから、多変量解析を用いて、前記ラマン散乱光の強度及び前記背景光の強度のそれぞれの時間依存性を抽出すると共に前記ラマン散乱光のスペクトルを前記背景光のスペクトルから分離することにより、前記分析用スペクトルを取得することができる分離部と
を備えたことを特徴とするかび毒の分析装置。
In a mold venom analyzer capable of analyzing mold venom using vibrational spectroscopy,
From an analytical sample produced from an initial sample in which the type of mycotoxin of interest may be produced, a spectroscopic analysis unit capable of obtaining an analytical spectrum using the vibrational spectroscopic analysis method,
The data of the spectrum for analysis and the data of the spectrum for reference obtained by using the vibrational spectroscopic analysis method from the standard sample containing the type of mycotoxin targeted for the analysis are input, and arithmetic processing using multivariate analysis is performed. By carrying out, a processing unit capable of extracting the spectrum of the target fungal venom of the analysis when the target fungal venom of the analysis is produced in the sample ,
The spectroscopic analysis unit,
A recording unit capable of recording a spectrum including the Raman scattered light and the background light by condensing the Raman scattered light scattered by the analysis sample irradiated with laser light and the background light excluding the Raman scattered light. When,
When the analysis sample is continuously irradiated with laser light over a plurality of continuous time intervals, from the spectrum recorded for each of the plurality of time intervals, using multivariate analysis, the intensity of the Raman scattered light. And separating the spectrum of the Raman scattered light from the spectrum of the background light while extracting the respective time dependence of the intensity of the background light, and a separation unit capable of acquiring the analysis spectrum . A mold venom analyzer characterized in that
前記演算処理部は、
前記参照用スペクトルのデータに複数通りの係数を乗算して得られた複数通りの積のデータのそれぞれを前記分析用スペクトルのデータに加算して得られた複数通りの和スペクトルのデータから、行列で表すことができる解析対象データを生成することができる生成部と、
多変量解析を用いて、前記解析対象データから前記分析の目的とする種類のかび毒のスペクトル及び信号強度プロファイルを抽出することができる解析部と
を備えたことを特徴とする請求項7に記載のかび毒の分析装置。
The arithmetic processing unit,
From a plurality of sum spectrum data obtained by adding each of a plurality of product data obtained by multiplying the reference spectrum data by a plurality of coefficients, a matrix A generation unit that can generate analysis target data that can be represented by
An analysis unit capable of extracting a spectrum and a signal intensity profile of a mold venom of a type intended for the analysis from the data to be analyzed by using a multivariate analysis. Mycotoxin analyzer.
JP2016126209A 2016-06-27 2016-06-27 Mycotoxin analysis method and mycotoxin analyzer Active JP6692226B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016126209A JP6692226B2 (en) 2016-06-27 2016-06-27 Mycotoxin analysis method and mycotoxin analyzer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016126209A JP6692226B2 (en) 2016-06-27 2016-06-27 Mycotoxin analysis method and mycotoxin analyzer

Publications (2)

Publication Number Publication Date
JP2018004252A JP2018004252A (en) 2018-01-11
JP6692226B2 true JP6692226B2 (en) 2020-05-13

Family

ID=60946108

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016126209A Active JP6692226B2 (en) 2016-06-27 2016-06-27 Mycotoxin analysis method and mycotoxin analyzer

Country Status (1)

Country Link
JP (1) JP6692226B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2572662B (en) * 2018-10-05 2020-06-03 Res & Innovation Uk Raman spectrometer
CN109342602A (en) * 2018-11-21 2019-02-15 中国计量大学 Aflatoxin detection device in a kind of liquid
JP7262345B2 (en) * 2019-08-30 2023-04-21 一般財団法人電力中央研究所 Material inspection device, material inspection method and material inspection program
CN112649414A (en) * 2020-11-30 2021-04-13 安徽大学 In-vivo detection method of fusarenol toxoid in wheat head with scab lesion
JP7440870B2 (en) 2021-01-22 2024-02-29 株式会社島津製作所 How to analyze data determined by two variables

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1263301A4 (en) * 1999-10-18 2005-06-08 Kerry Scott Lane Method and system for assay and removal of harmful toxins during processing of tobacco products
US9662047B2 (en) * 2010-08-05 2017-05-30 Massachusetts Institute Of Technology Portable raman diagnostic system
JP6080577B2 (en) * 2013-02-06 2017-02-15 日本分光株式会社 Noise removal method and noise removal apparatus
JP5780476B1 (en) * 2014-09-05 2015-09-16 株式会社分光科学研究所 Spectral quantification method, spectral quantification apparatus and program

Also Published As

Publication number Publication date
JP2018004252A (en) 2018-01-11

Similar Documents

Publication Publication Date Title
JP6692226B2 (en) Mycotoxin analysis method and mycotoxin analyzer
Yaseen et al. Raman imaging for food quality and safety evaluation: Fundamentals and applications
Cogdill et al. Single-kernel maize analysis by near-infrared hyperspectral imaging
Mitsutake et al. Raman imaging spectroscopy: history, fundamentals and current scenario of the technique
Pizarro et al. Determination of the peroxide value in extra virgin olive oils through the application of the stepwise orthogonalisation of predictors to mid-infrared spectra
EP2828643B1 (en) Method and apparatus for characterising samples by measuring light scattering and fluorescence
US20140085630A1 (en) Spectroscopic apparatus and methods for determining components present in a sample
Montemurro et al. Highly sensitive quantitation of pesticides in fruit juice samples by modeling four-way data gathered with high-performance liquid chromatography with fluorescence excitation-emission detection
Yande et al. Detection of pesticides on navel orange skin by surface-enhanced Raman spectroscopy coupled with Ag nanostructures
CN113056672A (en) Information processing apparatus, control method for information processing apparatus, program, calculation apparatus, and calculation method
Gabrielsson et al. Recent developments in multivariate calibration
Montemurro et al. Automatic generation of photochemically induced excitation-emission-kinetic four-way data for the highly selective determination of azinphos-methyl in fruit juices
US20080154549A1 (en) Noise-Component Removing Method
CN112098340A (en) Turquoise identification method based on hyperspectral imaging technology and assembly line process
Ong et al. Determination of aflatoxin B1 level in rice (Oryza sativa L.) through near-infrared spectroscopy and an improved simulated annealing variable selection method
Hara et al. Effect of Raman exposure time on the quantitative and discriminant analyses of carotenoid concentrations in intact tomatoes
Deore et al. FTIR combined with chemometric tools—a potential approach for early screening of grazers in microalgal cultures
JP6630944B2 (en) Raman scattered light observation method, Raman scattered light observation device, and program
Joshi et al. Au nanogap SERS substrate for the carbaryl pesticide determination in juice and milk using chemomterics
Caponigro et al. Single-drop technique for lactose prediction in dry milk on metallic surfaces: Comparison of Raman, FT–NIR, and FT–MIR spectral imaging
US11841373B2 (en) Information processing apparatus, method for controlling information processing apparatus, and program
Monago-Maraña et al. Second-order calibration in combination with fluorescence fibre-optic data modelling as a novel approach for monitoring the maturation stage of plums
Tseng et al. Internet-enabled near-infrared analysis of oilseeds
WO2023123329A1 (en) Method and system for extracting net signal in near-infrared spectrum
Zhang et al. Non‐destructive detection and recognition of pesticide residue levels on cauliflowers using visible/near‐infrared spectroscopy combined with chemometrics

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190522

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200129

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200129

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200323

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200406

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200414

R150 Certificate of patent or registration of utility model

Ref document number: 6692226

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250