"遺伝子の発現量を指標として用い、該遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、Mogat2、Alpl、Csf1、Klra2、Cdr2、Kcnq2、Prp15、Clec4e、Dusp13、Gpr77、Znf532、Oas1k、Per2、Upp1、Oas1a、Pth、Lrg1、Dnase2b、Vegfa、Hdc、Egr1、Clec4d、Alox5、Car2、Il1r2、Il8rb、Pla2g7、Dgat2、Mcpt10、Mcpt8、Mcpt8l2、Mcpt8l3、S100a9、Rab32、Emb、Fetub、Scgb3a1、Nlrp3、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、TC603221、RGD1565432、RGD1565374、RGD1565161、RGD1565983、RGD1310507、LOC690116、LOC691993、LOC689064、CF109441、A64P154430、A64P138390、A64P050137、ENSRNOT00000035459、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、Hbb、およびHbaa2からなる群から選択される少なくとも1種である、ストレスの評価方法。前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、およびHbbからなる群から選択される少なくとも1種である、請求項1に記載のストレスの評価方法。前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、Tspan8、Mmp9からなる群から選択される少なくとも1種である、請求項1又は2に記載のストレスの評価方法。前記遺伝子が、Tspan8、Mmp9、Add2、Eraf、Aqp1、Hbaa2、Hbb、およびAlas2の8種の遺伝子からなる群から選択される、請求項1に記載のストレスの評価方法。2種以上の前記遺伝子を指標として用いる、請求項1~4のいずれか1項に記載のストレスの評価方法。Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、Mogat2、Alpl、Csf1、Klra2、Cdr2、Kcnq2、Prp15、Clec4e、Dusp13、Gpr77、Znf532、Oas1k、Per2、Upp1、Oas1a、Pth、Lrg1、Dnase2b、Vegfa、Hdc、Egr1、Clec4d、Alox5、Car2、Il1r2、Il8rb、Pla2g7、Dgat2、Mcpt10、Mcpt8、Mcpt8l2、Mcpt8l3、S100a9、Rab32、Emb、Fetub、Scgb3a1、Nlrp3、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、TC603221、RGD1565432、RGD1565374、RGD1565161、RGD1565983、RGD1310507、LOC690116、LOC691993、LOC689064、CF109441、A64P154430、A64P138390、A64P050137、ENSRNOT00000035459、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、Hbb、およびHbaa2かからなる群から選択される少なくとも1種の遺伝子を含有する、ストレス評価マーカー。前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、およびHbbからなる群から選択される少なくとも1種である、請求項6に記載のストレス評価マーカー。前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、Tspan8、Mmp9からなる群から選択される少なくとも1種である、請求項6又は7に記載のストレス評価マーカー。前記遺伝子が、Tspan8、Mmp9、Add2、Eraf、Aqp1、Hbaa2、Hbb、およびAlas2の8種の遺伝子からなる群から選択される、請求項6に記載のストレス評価マーカー。2種以上の前記遺伝子を含有する、請求項6~9のいずれか1項に記載のストレス評価マーカー。少なくともTspan8、Mmp9、Add2、Eraf、Aqp1、Hbaa2、Hbb、およびAlas2の8種の遺伝子を含有する、請求項6~10のいずれか1項に記載のストレス評価マーカー。床敷飼育した後に、金網飼育することによる、ストレス負荷モデル動物の作成方法。床敷飼育を5日間以上行う、請求項12に記載の作成方法。金網飼育を5日間以上行う、請求項12又は13に記載の作成方法。モデル動物が、げっ歯類である、請求項12~14のいずれか1項に記載の作成方法。モデル動物が、ラット又はマウスである、請求項12~15のいずれか1項に記載の作成方法。請求項12~16のいずれか1項に記載の方法により作成される、ストレス負荷モデル動物。しかしながら、上記検査方法の何れも、対象者のストレス疲労状態、ストレス刺激の種類等は個人差が大きく、ストレス評価は簡便な方法かつ高精度で、個人的に持続して行うべきであるという課題が残された。一方、うつ病においても、脳内のモノアミン系の機能的変調があり、その変調は心身相関作用を通して、神経内分泌系、神経免疫系、自律神経系に少なからぬ影響を及ぼしていることが知られている。また、最近では、脳血流や脳内モノアミン受容体の変化なども指摘されているが、いずれも感度や再現性に問題が残る。また、従来の検査は、施行と評価には膨大な時間と労力が必要であり、簡便性という観点をも考慮すると、日常診療への応用はとうてい望むことが出来ないのが現状である。本発明者らは、上記課題を解決するために、白血球中の特定の遺伝子の発現量を定量することによって、新規なストレスの評価方法を提供することができることを見出し、本発明を完成した。本発明は、以下のとおりである。[1]遺伝子の発現量を指標として用い、該遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、Mogat2、Alpl、Csf1、Klra2、Cdr2、Kcnq2、Prp15、Clec4e、Dusp13、Gpr77、Znf532、Oas1k、Per2、Upp1、Oas1a、Pth、Lrg1、Dnase2b、Vegfa、Hdc、Egr1、Clec4d、Alox5、Car2、Il1r2、Il8rb、Pla2g7、Dgat2、Mcpt10、Mcpt8、Mcpt8l2、Mcpt8l3、S100a9、Rab32、Emb、Fetub、Scgb3a1、Nlrp3、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、TC603221、RGD1565432、RGD1565374、RGD1565161、RGD1565983、RGD1310507、LOC690116、LOC691993、LOC689064、CF109441、A64P154430、A64P138390、A64P050137、ENSRNOT00000035459、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、Hbb、およびHbaa2からなる群から選択される少なくとも1種である、ストレスの評価方法[2]前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、およびHbbからなる群から選択される少なくとも1種である、[1]に記載のストレスの評価方法。[3]前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、Tspan8、Mmp9からなる群から選択される少なくとも1種である、[1]又は[2]に記載のストレスの評価方法。[4]前記遺伝子が、Tspan8、Mmp9、Add2、Eraf、Aqp1、Hbaa2、Hbb、およびAlas2の8種の遺伝子からなる群から選択される、[1]に記載のストレスの評価方法。[5]2種以上の前記遺伝子を指標として用いる、[1]~[4]のいずれか1項に記載のストレスの評価方法。[6]Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、Mogat2、Alpl、Csf1、Klra2、Cdr2、Kcnq2、Prp15、Clec4e、Dusp13、Gpr77、Znf532、Oas1k、Per2、Upp1、Oas1a、Pth、Lrg1、Dnase2b、Vegfa、Hdc、Egr1、Clec4d、Alox5、Car2、Il1r2、Il8rb、Pla2g7、Dgat2、Mcpt10、Mcpt8、Mcpt8l2、Mcpt8l3、S100a9、Rab32、Emb、Fetub、Scgb3a1、Nlrp3、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、TC603221、RGD1565432、RGD1565374、RGD1565161、RGD1565983、RGD1310507、LOC690116、LOC691993、LOC689064、CF109441、A64P154430、A64P138390、A64P050137、ENSRNOT00000035459、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、Hbb、およびHbaa2かからなる群から選択される少なくとも1種の遺伝子を含有する、ストレス評価マーカー。[7]前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、およびHbbからなる群から選択される少なくとも1種である、[6]に記載のストレス評価マーカー。[8]前記遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、Tspan8、Mmp9からなる群から選択される少なくとも1種である、[6]又は[7]に記載のストレス評価マーカー。[9]前記遺伝子が、Tspan8、Mmp9、Add2、Eraf、Aqp1、Hbaa2、Hbb、およびAlas2の8種の遺伝子からなる群から選択される、[6]に記載のストレス評価マーカー。[10]2種以上の前記遺伝子を含有する、[6]~[9]のいずれか1項に記載のストレス評価マーカー。[11]少なくともTspan8、Mmp9、Add2、Eraf、Aqp1、Hbaa2、Hbb、およびAlas2の8種の遺伝子を含有する、[6]~[10]のいずれか1項に記載のストレス評価マーカー。[12]床敷飼育した後に、金網飼育することによる、ストレス負荷モデル動物の作成方法。[13]床敷飼育を5日間以上行う、[12]に記載の作成方法。[14]金網飼育を5日間以上行う、[12]又は[13]に記載の作成方法。[15]モデル動物が、げっ歯類である、[12]~[14]のいずれか1項に記載の作成方法。[16]モデル動物が、ラット又はマウスである、[12]~[15]のいずれか1項に記載の作成方法。[17][12]~[16]のいずれか1項に記載の方法により作成される、ストレス負荷モデル動物。本法は、通常の採血による数mLの血液をもとに解析可能であり、微侵襲的、簡便かつ日常的に行うことのできる検査方法である。以下、本発明を実施するための形態について以下詳細に説明する。なお、本発明は、以下の実施の形態に限定されるものではなく、その要旨の範囲内で種々変形して実施することができる。本発明は、遺伝子の発現量を指標として用い、該遺伝子が、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、Mogat2、Alpl、Csf1、Klra2、Cdr2、Kcnq2、Prp15、Clec4e、Dusp13、Gpr77、Znf532、Oas1k、Per2、Upp1、Oas1a、Pth、Lrg1、Dnase2b、Vegfa、Hdc、Egr1、Clec4d、Alox5、Car2、Il1r2、Il8rb、Pla2g7、Dgat2、Mcpt10、Mcpt8、Mcpt8l2、Mcpt8l3、S100a9、Rab32、Emb、Fetub、Scgb3a1、Nlrp3、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、TC603221、RGD1565432、RGD1565374、RGD1565161、RGD1565983、RGD1310507、LOC690116、LOC691993、LOC689064、CF109441、A64P154430、A64P138390、A64P050137、ENSRNOT00000035459、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、Hbb、およびHbaa2からなる群から選択される少なくとも1種である、ストレスの評価方法である。本発明において、ストレスの評価をするにあたり、上記111種の遺伝子のうち少なくとも1種の遺伝子発現量を測定することにより、ストレス負荷の有無を測定することができるが、以下のような方法によりストレスの負荷の有無を評価することができる。例えば、一遺伝子の発現量を測定して、ストレスの評価を行ってもよいし、複数遺伝子の発現量を同時に測定して、ストレスの評価を行ってもよい。また、111種遺伝子から選択される、1種又は複数種を含むマイクロチップを用いて、遺伝子発現量を測定してもよい。複数種の遺伝子を用いてストレスの評価を行う場合には、全ての遺伝子の発現量が増えている場合にはもちろん、発現量の亢進されている遺伝子の数が多いことをもってストレス負荷状態にあると評価することもできる。本発明のストレスの評価方法により、ストレスの有無を評価することができるが、かかる評価結果に基づいて、医師などによって、ヒトにおけるうつ病などの疾患に罹っているか否かを診断することができる。すなわち、本発明のストレスの評価方法においては、医師などが診断を行う際の基礎資料となる、血中のストレスマーカーである111種の遺伝子について、その少なくとも1種を用いてストレスの有無を評価するのであって、ヒトの病気を診断する方法には該当しない。ストレスの有無を評価するに当たり、かかる111種の遺伝子の少なくとも1種の遺伝子の発現量を測定する。かかる遺伝子の発現量について測定するサンプルは、特に限定されるものではないが、血液サンプルであることが好ましい。中でも、血液中の白血球における遺伝子の発現量を測定することが好ましい。遺伝子の発現量を測定し、ストレスの有無を評価するにあたっては、以下のような方法により行うことができる。サンプル中における遺伝子の発現量の絶対量を測定して、その多寡によりストレスの有無を測定することもできるが、データベース化した遺伝子発現量と対比することにより、ストレス評価を行うこともできる。なお、本発明においては、少なくとも1種の遺伝子の発現量を用いてストレス評価を行うが、うつ病患者やストレス負荷家畜からマーカー遺伝子群の発現量データを集め、データベース化し、そのデータベースと照合し、例えば、パターン認識技術を用いることによりストレス評価を行うこともできる。また、個人個体単位でマーカー発現量を定期的に測定記録しデータベース化し、そのデータベースと照合し、ストレス評価を行うこともできる。本発明における、遺伝子の発現量の測定方法としては、発現量が亢進していることを測定するために、mRNAの量を遺伝子の発現量として測定することもでき、また、遺伝子発現産物であるタンパクの発現変動を遺伝子の発現量として測定することもできる。遺伝子の発現量を転写産物としてのmRNAの量として測定する場合には、公知の遺伝子工学および分子生物学的技術に従い、当該分野で特定の遺伝子の発現を検知測定するために知られた手法、例えばinsituハイブリダイゼーション、ノーザンブロッティング、ドットブロット、RNaseプロテクションアッセイ、RT-PCR、Real-TimePCR等、当該技術分野で知られた任意の解析方法を用いることができる。また、遺伝子の発現変動を判定するための別の方法として、遺伝子の発現産物であるタンパク質量を定量する場合には、公知の遺伝子工学および分子生物学的技術に従い、当該分野で特定のタンパク質量を検知測定するために知られた手法、例えばウェスタンブロッティング、各種の免疫組織学的方法等、当該技術分野で知られた任意の解析方法を用いることができる。具体的には、遺伝子変動を判定するための手法としてタンパク質を定量する場合の例示として、Tspan8を例として説明すると、TSPAN8タンパク質を特異的に認識する抗体を用いて実施することができる。抗体はポリクローナル抗体またはモノクローナル抗体であり、それぞれTSPAN8タンパク質のエピトープに結合することができる全体分子、およびFab、Fab’2、Fv断片等が全て含まれる。このような抗体は、例えばポリクローナル抗体の場合には、タンパク質やその一部断片を免疫原として動物を免役した後、血清から得ることができる。あるいは、Tspan8遺伝子を組込んだ真核細胞用発現ベクターを注射や遺伝子銃によって、動物の筋肉や皮膚に導入した後、血清を採取することによって作製することができる。動物としては、マウス、ラット、ウサギ、ヤギ、ニワトリなどが用いられる。このような抗体を使用する方法は、例えば免疫染色、例えば組織あるいは細胞染色、免疫電子顕微鏡、イムノアッセイ、例えば競合型イムノアッセイまたは非競合型イムノアッセイで行うことができ、放射免疫測定法RIA、蛍光免疫測定法FIA、ルミネッセント免疫測定法LIA、酵素免疫測定法EIA、ELISAなどを用いることができ、B-F分離を行ってもよいし、あるいは行わないでその測定を行うことができる。好ましくはRIA、EIA、FIA、LIAであり、さらにサンドイッチ型アッセイが挙げられる。サンドイッチ型アッセイには、同時サンドイッチ型アッセイ、フォワードforwardサンドイッチ型アッセイあるいは逆サンドイッチ型アッセイなどであってもよい。Tspan8遺伝子の発現量をTSPAN8タンパク質を定量することによる方法を例示して説明したが、Tspan8遺伝子以外の遺伝子についても、同様の手法により遺伝子の発現量をタンパク質の定量により測定することもできる。本発明の評価方法においては、以上に例示した各方法の2以上を組み合わせて行うことによって、遺伝子の発現量を解析するにあたりより精度を高く判定することもできる。本発明においては、111種の遺伝子の発現量は、血液中の白血球画分における発現量を測定することが好ましい。本発明は、別の態様として、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、Mogat2、Alpl、Csf1、Klra2、Cdr2、Kcnq2、Prp15、Clec4e、Dusp13、Gpr77、Znf532、Oas1k、Per2、Upp1、Oas1a、Pth、Lrg1、Dnase2b、Vegfa、Hdc、Egr1、Clec4d、Alox5、Car2、Il1r2、Il8rb、Pla2g7、Dgat2、Mcpt10、Mcpt8、Mcpt8l2、Mcpt8l3、S100a9、Rab32、Emb、Fetub、Scgb3a1、Nlrp3、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、TC603221、RGD1565432、RGD1565374、RGD1565161、RGD1565983、RGD1310507、LOC690116、LOC691993、LOC689064、CF109441、A64P154430、A64P138390、A64P050137、ENSRNOT00000035459、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、Hbb、およびHbaa2からなる群から選択される少なくとも1種の遺伝子を含有する、ストレス評価マーカーである。ストレス評価マーカーとして、上記111種の遺伝子のうち少なくとも1種を含めば特に遺伝子が限定されるものではないが、111種の遺伝子それぞれを単独でマーカーとして用いてもよく、2種以上の遺伝子群をマーカーとして用いてもよい。複数の遺伝子をマーカーとして用いることにより、より信頼性の高いストレス評価マーカーとして用いることができる。本発明においては、ストレス評価マーカーとは、かかるマーカーの発現量を測定することにより、ストレス負荷の状態にあるか否かを判定することのできるマーカーであることを意味し、かかるマーカーとしては、111種の遺伝子から選択される少なくとも1種である。ストレス評価マーカーとして、2種以上を用いる場合には、キットのように、かかるマーカーを用いてもよく、同時に測定してもよく、一方を先行してその発現量を測定し、次いで、他のマーカーの発現量を測定することにより、ストレス評価用のマーカーとして用いてもよい。ストレスマーカーとしては、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、Schip1、Mycbpap、Aqp9、Ldhc、Degs2、Oas3、Slc22a4、Slc38a5、Slc43a3、Rsad2、Slc4a1、Wfdc15b、Plekhb1、Cd52、Anxa1、Rhd、Nrtn、Alox15、Scube1、Ube2c、Mrgprx3、Steap4、Selenbp1、Rrm2、S100a8、Cfd、Tuba3b、Slpi、Padi4、Fech、Bpgm、Nlrp12、Csf3r、Osgin1、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、RGD1311874、MGC72973、TC587327、LOC688574、Tspan8、Mmp9、Add2、Eraf、Aqp1、Alas2、およびHbbからなる群から選択される少なくとも1種を含有することが好ましく、Retnlg、Csta、Il1b、Camp、Ifitm6、Ifit1lb、Fcar、Pglyrp1、Mmp8、TC646750、ENSRNOT00000065312、ENSRNOT00000040211、ENSRNOT00000056434、Tspan8、Mmp9からなる群から選択される少なくとも1種を含有することがより好ましく、Retnlg、Csta、Il1b、Camp、Ifit1lb、Fcar、Pglyrp1、Mmp8、Tspan8、およびMmp9からなる群から選択される少なくとも1種を含有することがさらに好ましい。マイクロアレイのように、これら遺伝子を全てマーカーとして有していてもよく、また、その一部に遺伝子の少なくとも1種を含んでいてもよい。本発明の評価方法において用いられるストレス評価マーカーである111種の遺伝子は、本発明のもう1つの態様であるストレス負荷モデル動物の作成方法により見出されたものである。本発明において、ストレス負荷モデル動物の作成方法により、新規ストレス負荷モデル動物を作成し、ストレス負荷された動物における血液から白血球成分を抽出し、かかる血球成分中の遺伝子発現量を測定することにより見出された遺伝子群である。具体的には実施例に示すような実験により、ストレス検出用のマーカーとして選択された111種の遺伝子である。表1および表2に示す。表1及び表2は、各遺伝子についてのAccessionnumberとその機能と共に、白血球の遺伝子発現量をDNAマイクロアレイN数=5で分析した結果も示すものである。本発明におけるストレス負荷モデル動物の作成方法は、好ましくは、げっ歯類、中でも、マウスやラットに適用することで、新規ストレス負荷モデル動物を作成することができる。モデル動物を購入した際には、通常、輸送ストレス等がかかっていることが多い。そこで、本発明においては、まず、モデル動物を床敷条件で飼育して、かかるストレスを解消させる。床敷条件での飼育には、通常、5日以上飼育することにより、購入から輸送までの間にかかるストレスを解消させることができる。かかる飼育期間は、7日以上であることが好ましい。床敷条件での飼育とは、モデル動物を、例えば、夏目製作所製のポリカーボネートケージのようなケージに、ストレス負荷をより軽減する観点で、図1に示すように、ウッドチップやおがくずなどの床材や市販の床敷などを敷いた床敷ケージで飼育することを意味する。床敷条件で飼育することによりストレスの解消したモデル動物を、続いて、金網条件で飼育することによる、ストレス負荷モデル動物を作成することができる。床敷条件から、金網条件下でと、飼育環境が変わることにより、ストレス負荷のかかったモデル動物を作成することができる。金網条件での飼育には、通常、5日以上飼育することにより、ストレスから解消されたモデル動物にストレスを負荷することができる。かかる飼育期間は、7日以上であることが好ましい。金網条件での飼育は、長く行うと、モデル動物がストレス環境に順応してしまうため、好ましくは、10日以内で飼育することが好ましい。金網条件での飼育とは、モデル動物を、1.2cm角程度のメッシュの開いた金網ケージで飼育することを意味する。ストレス負荷モデル動物の作成の確認については、金網条件での飼育下のモデル動物と、床敷条件で飼育を続けたコントロールとしてのモデル動物との血液サンプル中の好中球数の違いや、ストレスマーカーとして知られた遺伝子の発現量を比較することによっても確認することができる。また、本発明により、111種の遺伝子マーカーが提供されているので、これら111種の遺伝子の発現量を比較することによっても、ストレス負荷モデル動物の作成を確認することができる。本発明においては、上記ストレス負荷モデル動物の作成方法により作成されるストレス負荷モデル動物をも提供する。以下、実施例を用いて、本発明をより詳細に説明するが、本発明はこれらの実施例に限定されるものではない。[実施例1]ストレス負荷ラットの作成1.試験方法購入した近交系の雄のラットFischer344の輸送ストレスを排除するために、7日間、床敷ケージ内でラットを飼育した。その後、同じ飼育部屋の別の床敷ケージ又は金網ケージへ、ラットを移し飼育した図1参照。それぞれのケージでの飼育7日目に採血し、得られたサンプル0.2mLを用いてヒト用血球分析装置XT-2000i、Sysmex製で血液成分解析を行った。床敷飼育については11匹の、金網飼育については12匹のラットの血液中の好中球の割合の測定結果およびその平均を表3に示す。図2は、血液中の好中球の割合の平均値をグラフ化したものであり、縦軸は、血液中の好中球の割合を表す。2.結果ラットに対して肉体的拘束などのストレスを加えると血液中の好中球の増加が観察されることが知られている新潟大学医学部医動物学教室のグループによって実証されてもいる:日本経済新聞19951124。また、日本補完代替医療学会誌Vol1No1pp3140や、臨床検査ガイド2007~2008文光堂出版なども参考になる。金網ケージ内で飼育したラットと床敷ケ-ジ内で飼育したラットとの血液成分を比較解析したところ、前者の白血球中好中球の割合は平均15.1であり後者は12.1であった。すなわち、金網飼育において、床敷飼育に比べ顕著に白血球中の好中球の割合が増加していたことから表3および図2参照、金網飼育により、ストレス負荷のかかったラットが得られた。[実施例2]CREB遺伝子およびIRF2遺伝子の白血球中での発現解析CREBcAMPResponseElementBindingprotein遺伝子およびIRF2InterferonRegulatoryFactor2遺伝子はうつ病やストレスと密接に関連していることが知られている。うつ病など精神疾患を対象に、末梢血白血球遺伝子発現を指標にして病因候補遺伝子を探索した結果、うつ病患者においてCREB遺伝子、つまりCREBmRNAの発現量が上昇することが報告されているJunichiIgaetalAlteredHDAC5andCREBmRNAexpressionsintheperipheralleukocytesofmajordepressionProgressinNeuroPsychopharmacologyandBiologicalPsychiatryVol31No3pp6286322007。また、健常人にストレスを負荷した際に、末梢血白血球由来のIRF2遺伝子、つまりIRF2mRNAの発現量が上昇することも報告されている特開2006-15号公報。そこで、実施例1で得られた血液サンプルを用いて白血球中CREB遺伝子およびIRF2遺伝子の発現変化の有無を確認した。<白血球の血液サンプルからの分離>床敷飼育したラット及び金網飼育したラットから採取した静脈血5mLを血液凝固阻止剤ヘパリンナトリウム注射液、100単位/mLを50μLで処理後、6mLの血球分離試薬モノポリ分離溶液、大日本住友製薬ファーマバイオメディカルの入った試験管へ移した。遠心分離1500rpm、20minを連続して3回、トータル60minを行い、白血球が含まれる白色の層を1.5mLチューブへ移した。生理食塩水を加え、遠心分離6000rpm、5minを3回を行って、白血球を沈殿させた。<白血球からのRNA抽出>次に、RNeasyminikitQiagenを用いて、添付説明書に従いtotalRNAの抽出を行った。具体的には以下のとおりである。BufferRLTにメルカプトエタノールM-Eを10μL/mLの濃度で添加した。遠心分離により沈殿させて得られた白血球サンプルにBufferRLTM-Eを700μL加えた。ボルテックスミキサーにより、物理的に細胞を破壊した。さらに、シリンジに取り付けた先の尖っていない注射針20-G、直径0.9mmに白血球サンプルを少なくとも5回通し、物理的に白血球を破壊した。得られたサンプルに600μLの70エタノール水溶液を添加し、ピペッティングによりよく混和した。600μLのサンプルを、2mLチューブ中にセットしたRNeasyスピンカラムに適用し、8000×g10,000rpm以上で15秒間遠心分離し、上清を除去した。350μLの洗浄液RW1bufferをRNeasyスピンカラムに添加し、8000xg10,000rpm以上で15秒間遠心分離し、上清を除去した。DNAのコンタミネーションを防ぐため、1サンプルあたりDNA分解酵素DNase20μLとバッファーRDD140μLの計160μLをカラム中央に適用し、室温で15分間静置した。350μLの洗浄液RW1bufferをRNeasyスピンカラムに添加し、8000xg10,000rpm以上で15秒間遠心分離し、上清を除去した。RNeasyスピンカラムに500μLの洗浄液RPEbufferを添加し、8000xg10,000rpm以上で15秒間遠心分離し、上清を除去した。再び、RNeasyスピンカラムに500μLの洗浄液RPEbufferを添加し、8000xg10,000rpm以上で2分間遠心操作した。RNeasyスピンカラムを新しい1.5mLチューブにセットした。蒸留水40μLを直接スピンカラムメンブレンに添加し、8000xg10,000rpm以上で1分間遠心操作を行ない、RNAを溶出した。新しい蒸留水40μLを用いて上記ステップを再度行ない、RNAを溶出した。抽出した各totalRNAサンプル50μLを濃度測定に用いた。残りのRNAサンプルは-80°Cにて凍結保存した。抽出したtotalRNAは吸光度測定にてRT-PCR解析に用いられる品質であるかどうかを確認した。<totalRNAからのcDNAの合成>逆転写反応用キットTOYOBO社製ReverTraAceを用いて、42°C、20分→99°C、5分→4°Cで逆転写反応を行った。逆転写反応用の試薬15sample分として、水150μL、buffer5xRT60μL、dNTP30μL、OligodT15μL、RNaseInhibitorTOYOBO15μL、ReverTraAseTOYOBO15μL、TemplateRNA15μLそれぞれ1μL117.5ngを用い、totalvol.20μLで実施した。<PCR法による遺伝子検出実験>CREB遺伝子およびIRF2遺伝子に対して特異的なプライマーを用いて、白血球由来cDNAを鋳型として、ExTaqを使用し、下記プロトコールに従いPCR反応を行った。PCRの反応液量totalvol.は1sampleあたり25μlで実施した。PCR用反応液30sample分として、DNAポリメラーゼ酵素Ex-Taq3.75μL、10xPCRbuffer75μL、dNTPs60μL、白血球由来templatecDNA60μL、滅菌水492μLを用い、混合液を23μLずつPCR用チューブへ分注し、それぞれの遺伝子検出用のプライマーDNAをそれぞれ2μLずつ添加した。PCR反応は、以下の条件で実施した。IRF2遺伝子のアニーリング温度は54°C、サイクル数は29、CREB遺伝子のアニーリング温度は56°C、サイクル数は30で実施した。RT-PCR解析で使用したCREB遺伝子およびIRF2遺伝子の特異的プライマーとして、CREBForwardPrimer5’GGCCTGCAGACATTAACCAT3’配列番号1ReversePrimer5’ATCCAGTCCATTTTCCACCA3’配列番号2IRF2ForwardPrimer5’ATCTGAGCGACCTTCCAAGA3’配列番号3ReversePrimer5’AGGACCGCATACTCAGGAGA3’配列番号4を用いた。PCR反応後、4アガロースゲル電気泳動を行い、エチジウムブロマイド染色で目的遺伝子の検出をした。検出にはUVトランスイルミネーターを用い、遺伝子バンドの強度を測定した。RT-PCR解析の結果、CREB遺伝子およびIRF2遺伝子は顕著に発現量が上昇していることが確認できた図3参照。CREB遺伝子およびIRF2遺伝子の発現量が上昇していることから、金網飼育により、ストレス負荷のかかったラットが得られたことが確認できた。なお、コントロールとして用いたGAPDH遺伝子のラットGAPDH遺伝子検出用プライマー配列は以下の通りである。ForwardPrimer5’TTCAATGGCACAGTCAAGGC3配列番号5ReversePrimer5’TCACCCCATTTGATGTTAGCG3配列番号6[実施例3]マイクロアレイ解析実施例1と同様にして、ラットを飼育して、金網飼育によりストレス負荷のかかったラット5個体および床敷飼育したラット5個体を得た。金網又は床敷上での飼育7日目に血液サンプルを採取した。実施例2と同様にして、RNeasyminikitQiagenを用いて、添付説明書に従いtotalRNAの抽出を行ったそれぞれの血液サンプルに由来するtotalRNAサンプルを、以下、ストレス負荷サンプル、対照サンプルという。。次に、DNAマイクロアレイ解析を行った。cDNAの合成、cRNAのラベルと増殖をLowInputQuickAmplabelingKitAgilentTechnologies社製を用いて行った。まず、totalRNA25ng、T7promoterprimer1.8μLを加え、totalで3.3μLになるようにNucleasefreewaterを加えた。RNASpikeInKit2カラー用のSpike-AをCy3ラベルする対照サンプルに、Spike-BをCy5ラベルするストレス負荷サンプルに2μL加えた。65°Cで10分間、インキュベート後、5分間、氷上に置き、熱変性させた。熱変性させたRNAにcDNAマスターミックス5xFirstStrandBuffer2.0μL、0.1MDTT1.0μL、10mMdNTPmix0.5μL、AffinityScript1.2μLを4.7μl加え、40°Cで2時間インキュベートした。チューブをウォーターバスから取り出し、70°Cで15分間インキュベートし、反応を停止した後、氷上で5分間冷却した。対照サンプルにCyanine3Transcriptionmixを、ストレス負荷サンプルにCyanine5Transcriptionmixを各反応チューブに加えた。Transcriptionマスターミックスとして、Nucleasefreewater0.75μL、5xTranscriptionbuffer3.2μL、0.1MDTT0.6μL、NTPmix1.0μL、T7RNApolymerasemix0.21μL、Cyanine3or5CTP0.24μLを1反応の必要量Ttotal6.0μLとして用いた。40°Cのウォーターバス中で、遮光しながら2時間インキュベートした。次に、RNeasyminispincolumusQiagen社を用いて、ラベル化cRNAを調整した。Nucleasefreewater84μL、RLTbuffer350μL、エタノール250μLを順次反応チューブに加えた。700μLのcRNAサンプルをカラムに移し、13,000rpmで1分間遠心した。フロースルーを別のチューブに移し、カラムをコレクションチューブに戻した。500μLのRPEBufferを加え、13000rpmで1分間遠心した。フロースルーを別のチューブに移し、カラムをコレクションチューブに戻した。500μLのRPEBufferを加え、13,000rpmで1分間遠心した。カラムを1.5mLチューブに移した。30μLのNucleasefreewaterを加え、1分間静置した。13,000rpmで1分間遠心し溶出した。次に、溶出した溶液状態のcRNAを分光光度計で濃度測定し、以下の溶液量となるようcRNAターゲット溶液の調整を行った。以下のように2xcRNAターゲット溶液を調整した。1チューブあたりの反応量を示す。Cyanine3ラベル化cRNA0.825μgCyanine5ラベル化cRNA0.825μg10xBlockingagent11.0μLNucleasefreewaterTotal52.8μL2xcRNAターゲット溶液52.8μLに25xフラグメンテーションバッファーを2.2μL加え、60°Cのウォーターバスで30分間インキュベートした。フラグメンテーションをストップさせるために、2xハイブリダイゼーションバッファーを55μL加えた。MicroarrayRatGE4x44kv32color、AgilentTechnologiest社製に100μLアプライし、65°Cで約17時間、ハイブリオーブンで攪拌しながらハイブリダイゼーションを行った。室温のAgilentgeneexpression洗浄バッファー1で1分間洗浄した。37°CのAgilentgeneexpression洗浄バッファー2で1分間洗浄した。スライドを引き上げて風乾した。乾燥後、Agilenttechnologiesmicroarrayscannerを用いて5μmの解像度でスキャンした。次に、変動した遺伝子を選出するために、microarrayscannerから出力された数値生データをExcelにより解析した。Controlフィーチャーは、ControlType→0にして、全て排除した。つまり、解析に使用しないポジティブコントロール、ネガティブコントロールを排除し、解析に使用するデータのみを表示した。なお、gIsFeatNonUnifOLGreenChannelCy3のフィーチャーに対するNonUniformityOutlier、gIsFeatPopnOLGreenChannelCy3のフィーチャーに対するPopulationOutlierのフラグは排除し、また、バックグラウンドより充分にシグナルの高いフィーチャーを集めるために、バックグラウンド±13xSDよりも低いシグナルフィーチャーを排除した。P値Pvalueは0.01以下、つまり99の信頼性をもって対照サンプルとの発現差があると判定されたフィーチャーのみを表示した。P値はいずれも10E-2未満であり、Retnlg遺伝子、Csta遺伝子、Il1b遺伝子、Camp遺伝子、Ifitm6遺伝子、Ifit1lb遺伝子、Fcar遺伝子、Pglyrp1遺伝子、Mmp8遺伝子、Schip1遺伝子、Mycbpap遺伝子、Aqp9遺伝子、Ldhc遺伝子、Degs2遺伝子、Oas3遺伝子、Slc22a4遺伝子、Slc38a5遺伝子、Slc43a3遺伝子、Rsad2遺伝子、Slc4a1遺伝子、Wfdc15b遺伝子、Plekhb1遺伝子、Cd52遺伝子、Anxa1遺伝子、Rhd遺伝子、Nrtn遺伝子、Alox15遺伝子、Scube1遺伝子、Ube2c遺伝子、Mrgprx3遺伝子、Steap4遺伝子、Selenbp1遺伝子、Rrm2遺伝子、S100a8遺伝子、Cfd遺伝子、Tuba3b遺伝子、Slpi遺伝子、Padi4遺伝子、Fech遺伝子、Bpgm遺伝子、Nlrp12遺伝子、Csf3r遺伝子、Osgin1遺伝子、Mogat2遺伝子、Alpl遺伝子、Csf1遺伝子、Klra2遺伝子、Cdr2遺伝子、Kcnq2遺伝子、Prp15遺伝子、Clec4e遺伝子、Dusp13遺伝子、Gpr77遺伝子、Znf532遺伝子、Oas1k遺伝子、Per2遺伝子、Upp1遺伝子、Oas1a遺伝子、Pth遺伝子、Lrg1遺伝子、Dnase2b遺伝子、Vegfa遺伝子、Hdc遺伝子、Egr1遺伝子、Clec4d遺伝子、Alox5遺伝子、Car2遺伝子、Il1r2遺伝子、Il8rb遺伝子、Pla2g7遺伝子、Dgat2遺伝子、Mcpt10遺伝子、Mcpt8遺伝子、Mcpt8l2遺伝子、Mcpt8l3遺伝子、S100a9遺伝子、Rab32遺伝子、Emb遺伝子、Fetub遺伝子、Scgb3a1遺伝子、Nlrp3遺伝子、TC646750遺伝子、ENSRNOT00000065312遺伝子、ENSRNOT00000040211遺伝子、ENSRNOT00000056434遺伝子、RGD1311874遺伝子、MGC72973遺伝子、TC587327遺伝子、LOC688574遺伝子、TC603221遺伝子、RGD1565432遺伝子、RGD1565374遺伝子、RGD1565161遺伝子、RGD1565983遺伝子、RGD1310507遺伝子、LOC690116遺伝子、LOC691993遺伝子、LOC689064遺伝子、CF109441遺伝子、A64P154430遺伝子、A64P138390遺伝子、A64P050137遺伝子、ENSRNOT00000035459遺伝子、Tspan8遺伝子、Mmp9遺伝子、Add2遺伝子、Eraf遺伝子、Aqp1遺伝子、Alas2遺伝子、Hbb遺伝子、およびHbaa2遺伝子において、ストレス負荷ラットにおける発現量が有意に増加していた。5シリーズ全てにおいて発現上昇していた遺伝子は、Retnlg遺伝子、Csta遺伝子、Il1b遺伝子、Camp遺伝子、Ifitm6遺伝子、Ifit1lb遺伝子、Fcar遺伝子、Pglyrp1遺伝子、Mmp8遺伝子、TC646750遺伝子、ENSRNOT00000065312遺伝子、ENSRNOT00000040211遺伝子、ENSRNOT00000056434遺伝子、Tspan8遺伝子、およびMmp9遺伝子であった。表5および表6にマイクロアレイ解析における測定結果の一部を例示として示す。マイクロアレイ解析においては、表1および表2において、changesintranscriptlevelとして、測定値の平均値としても記載しているが、かかる数値は、以下のようにして求めた値である。Tspan8を例にかかる数値の計算方法を示すと、表5および表6に示されるように、Tspan8遺伝子の発現量の変化は、FeatureNumberとして示す、25273、36630、16892、のデータとして各々サンプル数に応じて測定される。そして、まず、各サンプル毎の平均値を取る。すなわち、Tspan8のSample1のデータとしては、Tspan8のFeatureNumberが25273、36630、16892、であるSample1のデータの平均値を取っている。そうすることにより、Tspan8の測定結果として、Sample1~Sample5までの各平均値としての、遺伝子発現量の変化が確認される。ここで、Sample1~Sample5までの全ての各平均値において、発現上昇が確認できた遺伝子が、上記5シリーズ全てにおいて発現上昇していた遺伝子である。そして、Sample1~Sample5として得られた遺伝子発現量の変化を総平均値として求めた値が、上記表1および表2におけるchangesintranscriptlevelとして示されているものである。4シリーズにおいて発現上昇していた遺伝子は、Schip1遺伝子、Mycbpap遺伝子、Aqp9遺伝子、Ldhc遺伝子、Degs2遺伝子、Oas3遺伝子、Slc22a4遺伝子、Slc38a5遺伝子、Slc43a3遺伝子、Rsad2遺伝子、Slc4a1遺伝子、Wfdc15b遺伝子、Plekhb1遺伝子、Cd52遺伝子、Anxa1遺伝子、Rhd遺伝子、Nrtn遺伝子、Alox15遺伝子、Scube1遺伝子、Ube2c遺伝子、Mrgprx3遺伝子、Steap4遺伝子、Selenbp1遺伝子、Rrm2遺伝子、S100a8遺伝子、Cfd遺伝子、Tuba3b遺伝子、Slpi遺伝子、Padi4遺伝子、Fech遺伝子、Bpgm遺伝子、Nlrp12遺伝子、Csf3r遺伝子、Osgin1遺伝子、RGD1311874遺伝子、MGC72973遺伝子、TC587327遺伝子、LOC688574遺伝子、Add2遺伝子、Eraf遺伝子、Aqp1遺伝子、Alas2遺伝子、およびHbb遺伝子であった。3シリーズにおいて発現上昇していた遺伝子は、Mogat2遺伝子、Alpl遺伝子、Csf1遺伝子、Klra2遺伝子、Cdr2遺伝子、Kcnq2遺伝子、Prp15遺伝子、Clec4e遺伝子、Dusp13遺伝子、Fcar遺伝子、Gpr77遺伝子、Znf532遺伝子、Oas1k遺伝子、Per2遺伝子、Upp1遺伝子、Oas1a遺伝子、Pth遺伝子、Lrg1遺伝子、Dnase2b遺伝子、Vegfa遺伝子、Hdc遺伝子、Egr1遺伝子、Clec4d遺伝子、Alox5遺伝子、Car2遺伝子、Il1r2遺伝子、Il8rb遺伝子、Pla2g7遺伝子、Bpgm遺伝子、Dgat2遺伝子、Mcpt10遺伝子、Mcpt8遺伝子、Mcpt8l2遺伝子、Mcpt8l3遺伝子、S100a9遺伝子、Rab32遺伝子、Emb遺伝子、Fetub遺伝子、Scgb3a1遺伝子、Nlrp3遺伝子、TC603221遺伝子、RGD1565432遺伝子、RGD1565374遺伝子、RGD1565161遺伝子、RGD1565983遺伝子、RGD1310507遺伝子、LOC690116遺伝子、LOC691993遺伝子、LOC689064遺伝子、CF109441遺伝子、A64P154430遺伝子、A64P138390遺伝子、A64P050137遺伝子、ENSRNOT00000035459遺伝子、およびHbaa2遺伝子であった。[実施例4]111種の遺伝子の中から、無作為に、Tspan8、Mmp9、Add2、Eraf、Aqp1、Hbaa2、Hbb、およびAlas2の8種類の遺伝子を選択した。実施例1と同様にして作成したストレス負荷ラットの血液サンプル中の、遺伝子発現量をLightCycler480装置登録商標、ロシュ製を用いて、RT-PCRにより測定した。RT-PCR用プライマーとして表2に示すようにプライマーを、NCBIhttpwwwncbinlmnihgovから対象遺伝子のDNA情報を取得した後、プライマーデザインソフトPrimer3httpfrodowimiteduprimer3にDNA情報をアプライし、いくつか候補となるプライマーセットを選択し、そのうちから、対象遺伝子のイントロンを挟むようにプライマーセットを設計した。PCR反応は、以下の条件で実施した。対照サンプルおよびストレス負荷サンプルをそれぞれ5匹測定し、平均値をとりmRNAの相対定量法により遺伝子発現量を比較した図4参照。図4において、左側のカラムwは、床敷飼育でのラット対照サンプルの結果であり、該発現量を1としたときの、右側のカラムsには、金網飼育でのラットストレス負荷サンプルにおける相対発現量を示す。"