"ユーカリ属植物のゲノムDNA中のシンナモイルCoAリダクターゼをコードする遺伝子領域内の複数の一塩基多型に基づいてユーカリ属植物の樹種を識別することを含む、ユーカリ属植物の種間雑種の樹種を識別する方法。前記ユーカリ属植物のゲノムDNA中のNS領域及び葉緑体ETS領域内の複数のSNPに基づいて上記ユーカリ属植物の種間雑種の樹種を識別することを行い、この識別結果と、上記CCR遺伝子領域のSNPを利用する識別結果とを組み合わせて該ユーカリ属植物の種間雑種の樹種を決定することをさらに含む、請求項1に記載の方法。前記方法が、前記ユーカリ属植物の種間雑種のゲノムDNA中の1又は2以上のSNPを含む領域を、各SNPを挟む該ゲノムDNAの配列に相補的な約20~約30塩基の配列からなるプライマーを用いるポリメラーゼ連鎖反応により増幅する工程、その増幅産物を配列決定し1又は2以上のSNPの各塩基を決定する工程、並びにSNPの塩基と樹種との関係を示すデータに基づいて該種間雑種の樹種を決定する工程を含む、請求項1又は2に記載の方法。識別可能なユーカリ属植物の樹種が、ユーカリグロブラス、ユーカリカマルドレンシス、ユーカリブラシアーナ、ユーカリテレティコルニスEtereticornis、ユーカリグランディス、ユーカリサリグナEsaligna、ユーカリデグルプタ、ユーカリペリータ、及びユーカリユーロフィラからなる群から選択される、請求項1~3のいずれか1項に記載の方法。前記NS領域が、ユーカリグランディスのゲノムDNA配列配列番号31の約180000番目の塩基から約181000番目の塩基までの領域に相当する領域であり、上記葉緑体ETS領域が、該ゲノムDNA配列配列番号31の約193500番目の塩基から約194000番目の塩基までの領域に相当する領域であり、及び上記CCR領域が、ユーカリグロブラスのゲノムDNA配列配列番号32に相当する領域である、請求項2~4のいずれか1項に記載の方法。前記SNPの塩基と樹種との関係を示すデータが、下記の表1NS及びETS領域を利用する場合及び表2CCR領域を利用する場合に示されるものである、ただし表中のスラッシュを挟む塩基は、各アレルの塩基を表す、請求項3に記載の方法。前記CCR領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーが、P220配列番号22、CCR02配列番号23、P223配列番号26、CCR001配列番号24、P221配列番号27、P225配列番号25、P219配列番号28、及びP217配列番号29のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される、請求項3~6のいずれか1項に記載の方法。前記NS領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーが、NS2配列番号2、NS3配列番号3、NS4配列番号4、NS6、NS11、及びNS13配列番号14のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される、請求項3~6のいずれか1項に記載の方法。前記ETS領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーが、ETS2配列番号19及びETS4配列番号21のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される、請求項3~6のいずれか1項に記載の方法。ユーカリ属植物の種間雑種の樹種を複数の一塩基多型を利用して識別するためのPCRプライマーセットを含むキットであって、該PCRプライマーセットが、ユーカリ属植物ゲノムのNS領域及びETS領域内のSNP1~SNP8並びにCCR領域内のSNP9~SNP17の各々を含む領域を増幅するためのプライマーを含み、該SNP1~SNP8がそれぞれ、ユーカリグランディスのNS及びETS領域のゲノムDNA配列配列番号30の180133番目、180769番目、193840番目、193651番目、193650番目、193624番目、193591番目、193545番目に相当する位置に存在し、該SNP9~SNP17がそれぞれ、ユーカリグロブラスのCCR領域のゲノムDNA配列配列番号31の989番目、1163番目、1210番目、1223番目、1846番目、1878番目、2312番目、2428番目、2548番目に相当する位置に存在し、該プライマーが、該SNP1~SNP17の各SNPを挟む該ゲノムDNAの塩基配列、又は該塩基配列と90以上の同一性を有する塩基配列、に相補的な約20~約30塩基長の配列からなるセンス及びアンチセンスプライマーのセットを含む、上記キット。前記PCRプライマーセットが、下記の~:SN領域内のSNP1について、増幅プライマーがNS11配列番号12NS13配列番号14の組合せ、及び又はNS11配列番号12NS2配列番号2の組合せ、及び又はNS11NS6の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;SN領域内のSNP2について、増幅プライマーがNS11配列番号12NS13配列番号14の組合せ、及び又はNS3配列番号3NS4配列番号4の組合せ、及び又はNS11NS6の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい並びに、ETS領域内のSNP3~SNP8の各々について、増幅プライマーがETS4配列番号21ETS2配列番号19の組合せ、及び又はETS18S配列番号16ETSMyrtF配列番号17の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP9について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP10について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP11について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP12について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP13について、増幅プライマーがCCR001配列番号24P219配列番号28の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP14について、増幅プライマーがCCR001配列番号24P219配列番号28の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP15について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP16について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;並びに、CCR領域内のSNP17について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;を含む、請求項10に記載のキット。ユーカリ属植物の種間雑種の樹種を複数の一塩基多型を利用して識別するための、SNP1~SNP8、SNP9~17と樹種との関係を示す、請求項6に示される表1及び表2の組合せの使用方法。上記の背景技術に記載されたように、植物の種間雑種または種内雑種を正確にかつ簡単に識別または判別するための手法が非常に少ないことが分かった。ユーカリ属植物についてもこれらの交雑種のなかに成長性、病害性などの特性の点で優良な品種が存在することが認められるため、そのような優良品種を選抜するためにひとつひとつ栽培して各交雑種の特性を調べるやり方では多大の時間と労力を要することは明らかである。実際に、例えば主要樹種であるユーカリグロブラスEglobulusおよびユーカリカマルドレンシスEcamaldulensisの種間雑種には、両者の利点を兼ね備えた新品種が得られる可能性があるため、形態的には雑種と思われる候補木が遺伝的にも雑種であるかを識別する手法を開発する必要性がある。また、ラオスをはじめとする東南アジア植林地では、外部から導入したユーカリ雑種クローンあるいは将来目標にあげている種間雑種の両親樹種が何であるかを調べる必要性が高まっているし、今後も東南アジアを中心とした国々で様々なクローンあるいは雑種が導入されることが予想され、樹種識別ツールは将来にわたって有用であると予想される。このため、本発明は、ユーカリ属植物の種間雑種を複数のSNPを利用して正確かつ容易に識別する方法を提供することを目的とする。本発明者らは、ユーカリ属由来の多数の樹種のサンプルを種々の産地から集め、全ゲノム中で近縁種間のSNPが存在すると考えられる3領域NS領域、ETS領域及びCCR領域を中心に解析して樹種間で差のある17種のSNPを見出し、そして、これらのSNPでユーカリ属9樹種を4~8グループに分類できることを見出し、本発明を完成させた。なお、NS領域とETS領域の8種のSNPを利用する樹種の識別については、この出願と同じ出願人による出願である特願2009295174号にも記載されている。したがって、本発明は、以下の特徴を有する。ユーカリ属植物のゲノムDNA中のシンナモイルCoAリダクターゼをコードする遺伝子領域内の複数の一塩基多型に基づいてユーカリ属植物の樹種を識別することを含む、ユーカリ属植物の種間雑種の樹種を識別する方法。上記ユーカリ属植物のゲノムDNA中のNS領域及び葉緑体ETS領域内の複数のSNPに基づいて上記ユーカリ属植物の種間雑種の樹種を識別することを行い、この識別結果と、上記CCR遺伝子領域のSNPを利用する識別結果とを組み合わせて該ユーカリ属植物の種間雑種の樹種を決定することをさらに含む、上記に記載の方法。上記方法が、前記ユーカリ属植物の種間雑種のゲノムDNA中の1又は2以上のSNPを含む領域を、各SNPを挟む該ゲノムDNAの配列に相補的な約20~約30塩基の配列からなるプライマーを用いるポリメラーゼ連鎖反応により増幅する工程、その増幅産物を配列決定し1又は2以上のSNPの各塩基を決定する工程、並びにSNPの塩基と樹種との関係を示すデータに基づいて該種間雑種の樹種を決定する工程を含む、上記又はに記載の方法。識別可能なユーカリ属植物の樹種が、ユーカリグロブラス、ユーカリカマルドレンシス、ユーカリブラシアーナ、ユーカリテレティコルニスEtereticornis、ユーカリグランディス、ユーカリサリグナEsaligna、ユーカリデグルプタ、ユーカリペリータ、及びユーカリユーロフィラからなる群から選択される、上記~のいずれかに記載の方法。上記NS領域が、ユーカリグランディスのゲノムDNA配列配列番号31の約180000番目の塩基から約181000番目の塩基までの領域に相当する領域であり、上記葉緑体ETS領域が、該ゲノムDNA配列配列番号31の約193500番目の塩基から約194000番目の塩基までの領域に相当する領域であり、及び上記CCR領域が、ユーカリグロブラスのゲノムDNA配列配列番号32に相当する領域である、上記~のいずれかに記載の方法。上記SNPの塩基と樹種との関係を示すデータが、下記の表1NS及びETS領域を利用する場合及び表2CCR領域を利用する場合に示されるものである、ただし表中のスラッシュを挟む塩基は、各アレルの塩基を表す、上記に記載の方法。上記CCR領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーが、P220配列番号22、CCR02配列番号23、P223配列番号26、CCR001配列番号24、P221配列番号27、P225配列番号25、P219配列番号28、及びP217配列番号29のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される、上記~のいずれかに記載の方法。上記NS領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーが、NS2配列番号2、NS3配列番号3、NS4配列番号4、NS6、NS11、及びNS13配列番号14のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される、上記~のいずれかに記載の方法。上記ETS領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーが、ETS2配列番号19及びETS4配列番号21のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される、上記~のいずれかに記載の方法。ユーカリ属植物の種間雑種の樹種を複数の一塩基多型を利用して識別するためのPCRプライマーセットを含むキットであって、該PCRプライマーセットが、ユーカリ属植物ゲノムのNS領域及びETS領域内のSNP1~SNP8並びにCCR領域内のSNP9~SNP17の各々を含む領域を増幅するためのプライマーを含み、該SNP1~SNP8がそれぞれ、ユーカリグランディスのNS及びETS領域のゲノムDNA配列配列番号30の180133番目、180769番目、193840番目、193651番目、193650番目、193624番目、193591番目、193545番目に相当する位置に存在し、該SNP9~SNP17がそれぞれ、ユーカリグロブラスのCCR領域のゲノムDNA配列配列番号31の989番目、1163番目、1210番目、1223番目、1846番目、1878番目、2312番目、2428番目、2548番目に相当する位置に存在し、該プライマーが、該SNP1~SNP17の各SNPを挟む該ゲノムDNAの塩基配列、又は該塩基配列と90以上の同一性を有する塩基配列、に相補的な約20~約30塩基長の配列からなるセンス及びアンチセンスプライマーのセットを含む、上記キット。上記PCRプライマーセットが、下記の~:SN領域内のSNP1について、増幅プライマーがNS11配列番号12NS13配列番号14の組合せ、及び又はNS11配列番号12NS2配列番号2の組合せ、及び又はNS11NS6の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;SN領域内のSNP2について、増幅プライマーがNS11配列番号12NS13配列番号14の組合せ、及び又はNS3配列番号3NS4配列番号4の組合せ、及び又はNS11NS6の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい並びに、ETS領域内のSNP3~SNP8の各々について、増幅プライマーがETS4配列番号21ETS2配列番号19の組合せ、及び又はETS18S配列番号16ETSMyrtF配列番号17の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP9について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP10について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP11について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP12について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP13について、増幅プライマーがCCR001配列番号24P219配列番号28の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP14について、増幅プライマーがCCR001配列番号24P219配列番号28の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP15について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP16について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;並びに、CCR領域内のSNP17について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;からなる、上記に記載のキット。ユーカリ属植物の種間雑種の樹種を複数の一塩基多型を利用して識別するための、SNP1~SNP8、SNP9~17と樹種との関係を示す上記に示される表1及び表2の組合せの使用方法。本発明により、ユーカリ属植物の種間雑種の樹種を実験室レベルで容易にかつ正確に識別することが可能になり、これにより、各雑種をひとつひとつ栽培してその特性を調べる必要性がなくなった。すなわち、約1000種以上も存在するといわれるユーカリ属植物の天然交雑種のなかから優良品種を従来法と比べてより効率よく正確に選抜することを本発明の方法が可能にする。樹種識別法ユーカリ属植物は、主にオーストラリア原産であるが、今日では、東南アジア、オーストラリア、南米などの地域に植林されるようになり、世界的に有用な産業造林樹種のひとつであり、特に、紙パルプ、製材、木炭、オイルなどとして利用されている。この植物には、多数の種間交雑種が存在する。種間交雑には、天然交雑、人工交配、その他の交雑が含まれる。天然交雑は、分布域と開花期が重なった場合に起こり、例えばユーカリペリータ×ユーカリブラシアーナ、ユーカリグロブラス×ユーカリカマルドレンシスなどでこの種の交雑が認められている。人工交配は、優良木同士の人工交配であり、例えばユーカリユーロフィラ×ユーカリグランディスなどで雑種育種が試みられている。多数存在する交雑種について、現在の種間雑種の判別法は形態的な指標によるものであり、より早期により確実に判別する方法としてDNAマーカーを利用する方法が求められてきた。このような状況のなかで、本発明は、ユーカリ属植物の種間雑種の樹種を複数のSNPを利用して識別する方法であって、該雑種のゲノムDNAを鋳型にしてSNP1~SNP17の各々を含む領域を増幅する工程、増幅産物を配列決定しSNP1~SNP17の各塩基を決定する工程、並びに表1及び表2に基づいて該雑種の樹種を判定する又は、該雑種の親種を特定する工程を含んでなる方法を提供する。本発明の方法によって識別可能な又はグループ分け可能なユーカリ属植物の樹種は、ユーカリグロブラス、ユーカリカマルドレンシス、ユーカリブラシアーナ、ユーカリテレティコルニスEtereticornis、ユーカリグランディス、ユーカリサリグナEsaligna、ユーカリデグルプタ、ユーカリペリータ、及びユーカリユーロフィラからなる9種である。本発明によれば、ユーカリ属植物の上記樹種のゲノムDNA中で特定されたSNP1~17からなる17種のSNPを利用して雑種の樹種又は親種を識別又はグループ分けすることができる。これらのSNPは3領域、すなわちNS領域、ETS領域及びCCR領域に存在しており、具体的には、NS領域に2種のSNP、ETS領域に6種のSNP、CCR領域に9種のSNPがそれぞれ存在している。これらのSNPでユーカリ属9樹種を、NS領域及びETS領域の8種のSNPを利用するとき5つのグループに、CCR領域の9種のSNPを利用するとき4つのグループに、そしてこれらの17種のSNPを利用するとき8つのグループに分類できる。すなわち、SNP1~SNP8を使用するとき、ユーカリ属植物の上記9樹種は以下の5つの群にグループ分け可能である。ユーカリグロブラス、ユーカリカマルドレンシス又はユーカリブラシアーナ、又はユーカリテレティコルニスEtereticornisユーカリグランディス、又はユーカリサリグナEsalignaユーカリデグルプタ、及びユーカリペリータ又はユーカリユーロフィラ。また、SNP9~SNP17を使用するとき、ユーカリ属植物の上記9樹種は以下の4つの群a群~d群にグループ分け可能である。ユーカリグロブラス、ユーカリテレティコルニスEtereticornis、ユーカリサリグナEsaligna、ユーカリカマルドレンシス、ユーカリブラシアーナ、ユーカリグランディス、ユーカリデグルプタ、ユーカリペリータ又はユーカリユーロフィラ。さらにまた、SNP1~SNP8とSNP9~SNP17の両方を組み合せて使用するとき、ユーカリ属植物の上記9樹種は以下の8つの群にグループ分け可能である。ユーカリグロブラス、ユーカリカマルドレンシスユーカリブラシアーナ、ユーカリテレティコルニスEtereticornis、ユーカリグランディス、ユーカリサリグナEsaligna、ユーカリデグルプタ、及びユーカリペリータ又はユーカリユーロフィラ。したがって、本発明の方法により上記の17種のSNPを利用することによって、ユーカリ属植物の9樹種のうち7種の樹種を識別することが可能であるし、また、残りの2樹種についても、SNP1SNP3SNP5SNP6SNP7SNP8SNP11SNP15SNP17の少なくとも1つによる判定結果に基づいて、場合によっては各樹種の識別が可能になるか、或いはそうでない場合には、その2樹種のいずれかであると判定可能である。上記のSNP1~SNP8のうち、SNP1とSNP2はユーカリゲノムDNAのNS領域に存在し、SNP3~SNP8はユーカリゲノムDNAの葉緑体ETS領域に存在する。ユーカリゲノム配列に基づくSNP1~SNP8の位置は次のとおりである。SNP1scaffold522180133番目T又はCSNP2scaffold522180769番目A又はGSNP3scaffold522193840番目A又はGSNP4scaffold522193651番目G又はASNP5scaffold522193650番目T又はCSNP6scaffold522193624番目G又はCSNP7scaffold522193591番目T又はCSNP8scaffold522193545番目A又はGここで、scaffold522は、ユーカリグランディスのゲノムDNA断片の番号を表す。httpeucalyptusdbbiupaczaは、Eucspressoが提供するユーカリ遺伝子発現データベースである。また、上記のSNP9~SNP17は、ユーカリゲノムDNAのCCR領域に存在する。このCCR領域の配列情報は、例えばユーカリグロブラスのCCR遺伝子の塩基配列についてGenBankAY656821配列番号31に記載されている。この塩基配列に基づくと、SNP9~SNP17は、それぞれ989番目、1163番目、1210番目、1223番目、1846番目、1878番目、2312番目、2428番目、2548番目に位置する。本発明の方法では、まず、樹種が不明のユーカリ属植物の組織例えば、葉、根、茎などからゲノムDNAを抽出する。具体的には、組織を粉砕し、フェノールクロロホルム法、グアニジンイソチオシアネート酸性フェノール法などのDNA抽出法を用いてDNAを抽出し、エタノールを加え沈殿するDNAを回収する。DNAの抽出は、市販のDNA抽出キットを使用して実施することができる。次に、抽出されたゲノムDNA検査すべき目的の領域を含むDNA断片でもよい。を鋳型にしてSNP1~SNP2、SNP3~SNP8、SNP9~SNP17の1又は2以上のSNPの各々を含む領域をポリメラーゼ連鎖反応によって増幅し、その後、各増幅産物に対し配列決定用プライマーを使用してPCRを行って相補鎖DNAを合成し、その塩基配列を決定し、これによってSNP1~SNP17の1又は2以上、好ましくはすべての各塩基を決定する。上記の9樹種のユーカリ属植物における上記のSNPを含むゲノム領域をPCRによって増幅するためのプライマーセットは、該ユーカリ属植物ゲノムの上記のNS領域及びETS領域内のSNP1~SNP8並びに上記のCCR領域内のSNP9~SNP17の各々を含む領域を増幅するためのプライマーを含み、該SNP1~SNP8がそれぞれ、ユーカリグランディスのNS及びETS領域のゲノムDNA配列配列番号30の180133番目、180769番目、193840番目、193651番目、193650番目、193624番目、193591番目、193545番目に相当する位置に存在し、該SNP9~SNP17がそれぞれ、ユーカリグロブラスのCCR領域のゲノムDNA配列配列番号31の989番目、1163番目、1210番目、1223番目、1846番目、1878番目、2312番目、2428番目、2548番目に相当する位置に存在し、該プライマーが、該SNP1~SNP17の各SNPを挟む該ゲノムDNAの塩基配列、又は、該塩基配列と90以上、好ましくは95以上、さらに好ましくは98以上の同一性を有する塩基配列、に相補的な約20~約30塩基長の任意の配列からなるセンス及びアンチセンスプライマーセットを含む。このようなプライマーセットを使用して上記SNPを含む領域を増幅したときの増幅産物のサイズは、少なくとも40~70塩基長であり、好ましくは少なくとも100塩基長であり、例えば200~1000塩基長又はそれ以上であってもよい。上記のプライマーの例は以下のとおりである。上記CCR領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーの例は、P220配列番号22、CCR02配列番号23、P223配列番号26、CCR001配列番号24、P221配列番号27、P225配列番号25、P219配列番号28、及びP217配列番号29のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される。上記NS領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーの例は、NS2配列番号2、NS3配列番号3、NS4配列番号4、NS6、NS11、及びNS13配列番号14のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される。上記ETS領域のSNPの各々を含む領域をPCRによって増幅するためのプライマーは、ETS2配列番号19及びETS4配列番号21のプライマー、並びに、これらのプライマーの各塩基配列と90以上の同一性を有する塩基配列からなるプライマーからなる群から選択される。増幅工程のプライマーの組合せは、SNP1~SNP8及びSNP9~SNP17の各塩基を含む領域を増幅するのに適したユーカリゲノム配列由来のセンスプライマー及びアンチセンスプライマーの組合せであり、例えば、下記の~に示されるプライマーの組み合わせである。SN領域内のSNP1について、増幅プライマーがNS11配列番号12NS13配列番号14の組合せ、及び又はNS11配列番号12NS2配列番号2の組合せ、及び又はNS11NS6の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;SN領域内のSNP2について、増幅プライマーがNS11配列番号12NS13配列番号14の組合せ、及び又はNS3配列番号3NS4配列番号4の組合せ、及び又はNS11NS6の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい並びに、ETS領域内のSNP3~SNP8の各々について、増幅プライマーがETS4配列番号21ETS2配列番号19の組合せ、及び又はETS18S配列番号16ETSMyrtF配列番号17の組合せであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP9について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP10について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP11について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP12について、増幅プライマーがPP20配列番号22P223配列番号26の組み合わせ又はCCR002配列番号23P221配列番号27の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP13について、増幅プライマーがCCR001配列番号24P219配列番号28の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP14について、増幅プライマーがCCR001配列番号24P219配列番号28の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP15について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;CCR領域内のSNP16について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい;並びに、CCR領域内のSNP17について、増幅プライマーがP225配列番号25P217配列番号29の組み合わせであり、但し各組合せのプライマーは、独立に、前記配列番号の各塩基配列と90以上の同一性を有する塩基配列からなるプライマーであってもよい。上記プライマーは、通常、塩基数17~30、好ましくは20~30、より好ましくは20~25のサイズのDNAからなり、また、自然界で植物の突然変異が起こりうることを考慮して、それぞれの上記配列番号の塩基配列と90以上、好ましくは95以上、より好ましくは98以上の同一性を有する塩基配列からなるプライマーを使用することができる。変異は、1、2又は3個の塩基の置換、欠失又は挿入若しくは、付加を含む。配列に関する同一性は、2つの塩基配列を、最大の一致率となるようにギャップを導入して又はギャップを導入しないで整列アラインメントしたときに、ギャップを含めての塩基の総数に対する同一塩基数の割合を表す。同一性は、BLAST、FASTAなどの公知のアルゴリズムを使用して決定することができるAltschulSFetalNucleicAcidsRes2533893402McGinnisSMaddenTLNucleicAcidsRes32W20W25等。PCRは、変性、アニーリング及び伸長反応を1サイクルとして約20~40サイクルを繰り返すことによって行うことができる。変性は、二本鎖DNAを加熱によって解離して一本鎖にする、すなわちセンス鎖とアンチセンス鎖に分離する、ステップであり、そのための条件は、例えば約94~96°Cの温度で約10~60秒間二本鎖DNAを含む溶液を加熱処理することを含む。アニーリングは、センス鎖及びアンチセンス鎖の各鋳型DNAに、それらに対応するプライマー、すなわちフォワードプライマー又はリバースプライマー、を結合するステップであり、そのための条件は、鋳型DNAとプライマーを含む溶液を約50~65°Cで約10~60秒間加熱処理することを含む。伸長反応は、鋳型DNAにアニーリングしたプライマーを起点にして相補的DNA鎖を合成するステップであり、そのための条件は、耐熱性DNAポリメラーゼ及びdNTPsの存在下で、プライマーがアニーリングした鋳型DNAを含む溶液を約60~80°Cで約10~60秒間加熱処理することを含む。PCRの緩衝液は、dNTPs及びMgCl2を含むPCR用緩衝液であり、伸長反応に際して、このなかにさらに、耐熱性DNAポリメラーゼ、例えばTaqポリメラーゼThermusaquatics由来やPfuポリメラーゼPyrococcusfuriosus由来など、が添加される。PCRは市販のサーマルサイクラーを使用して実施することによって自動的なDNA増幅が可能になる。サーマルサイクラーは、例えば宝酒造、AppliedBiosystems、BioRadLabsなどから市販されている。本発明の方法では、増幅されたDNAを回収、精製し、さらに配列決定のために、精製DNAを鋳型にしてPCRを行い、増幅産物のDNA配列を決定する。DNAの回収、精製のために、ポリアクリルアミドゲル電気泳動、アガロースゲル電気泳動などの電気泳動の使用のほかに、市販のPCR産物精製キット、例えばQIAquickPCRpurificationkitQIAGEN、などを使用すると便利である。また、塩基配列の決定は、ジデオキシ法を利用した自動DNAシークエンサーを使用する方法などによって行うことができる。Sanger法では、上記のPCR条件ただし、4種類の基質のほかに、蛍光標識された特異的競合阻害剤を加える。で鋳型DNAに相補的なDNAを合成するが、この合成の際に基質と阻害剤の間で競合反応が起こり、阻害剤が取り込まれると反応が停止し、生成する相補鎖DNAはその3末端がいずれもddNTPNATG又はCで停止した種々の長さのDNA断片として得られ、これをゲル電気泳動でサイズ分離し、それを読み取って塩基配列を決定する。読み取りに際しては、レーザー光を励起し、その蛍光を検出するなどの手法を使用できる。キット及び使用本発明はさらに、ユーカリ属植物の種間雑種の樹種を複数の一塩基多型を利用して識別するためのPCRプライマーセットを提供する。該PCRプライマーセットは、NS領域及びETS領域内のSNP1~SNP8及びCCR領域内のSNP9~SNP17の各々を含む領域を増幅するための上記のプライマー並びに、必要に応じて、増幅産物を配列決定するためのプライマーを含み、好ましい例として上記の~のプライマーを含むセットである。本発明はさらに、上記のプライマー、好ましくは上記の~のプライマーセットを含む、ユーカリ属植物の種間雑種の樹種を識別するためのキットを提供する。キットには、上記プライマーセットのほかに、PCRのために必要な試薬類、例えば、PCRバッファー、耐熱性ポリメラーゼ、dNTPsなどの試薬、鋳型DNA調製用試薬、例えばゲノムDNA抽出用溶媒、制限酵素などの試薬、並びに使用説明書手順や操作法を含む。を含めることができる。上記のようにして決定されたテストサンプルのSNP1~SNP8及びSNP9~SNP17の各塩基について、表1及び表2に示されたそれぞれSNP1~SNP8及びSNP9~SNP17の各SNP塩基と樹種との関係に基づいて、後述の実施例の樹種の識別例のようにして、検査すべき交雑種の親種の樹種を決定することができる。例えば、テストサンプルのSNP1~SNP8の塩基がそれぞれTCGGAGGATCGCCCAAであれば、そのサンプルの樹種は、Ecamaldulensis×Edeglputaであると推定される。また、例えば、ユーカリカマルドレンシスとユーカリブラシアーナは、SNP1~SNP8を利用しても識別がむずかしいが、SNP1~SNP8とSNP9~SNP17例えば、SNP15及びSNP17の組合せを利用することによって、それらの識別が可能になる。したがって、本発明には、ユーカリ属植物の種間雑種の樹種を複数の一塩基多型を利用して識別するための、SNPと樹種との関係を示す表1及び表2の組合せの使用方法も包含される。この表に記載されるSNP樹種の関係を利用することによって、上記のとおり、ユーカリ属植物の種間雑種の樹種を識別することが可能になる。本発明を以下の実施例によってさらに説明するが、本発明の範囲はこれらの実施例によって制限されないものとする。[実施例1]ユーカリ属の樹種間で差のあるSNPの探索ユーカリ属植物の全ゲノム中で近縁種間のSNPが存在すると考えられる2領域を中心に解析し、それらの領域に含まれる複数のSNPを組み合わせることにより、未知の樹種の識別の信頼性を高めることを目的とする。これらの2領域すなわち、領域1及び領域2を、EcamaldulensisとEglobulusを例に以下に示す。ここで、領域1は、SNP1とSNP2が存在する領域であり、scaffold522として知られるEgrandisのゲノムDNA配列配列番号30の約180000番目の塩基から約181000番目の塩基までの領域に相当する領域であり、領域2は、SNP3~SNP8の各SNPが存在する領域であり、該ゲノムDNA配列の約193500番目の塩基から約194000番目の塩基までの領域に相当する領域である。SNP1~SNP8のゲノム上の各位置は、雑種によって幾分異なる場合もあると予想されるが、これらの領域は比較的保存性が高いため、当業者であれば、SNP周辺の配列との対比から各SNPの位置を特定することは容易である。材料発見したSNPが個体間ではなく樹種間で差のあるものであることを確かめるため、それぞれの樹種で大量のサンプルについて調査した。解析に使用したサンプルは、表3に示した。方法ユーカリ属全DNAの抽出乾燥したユーカリ属の葉約50ngをジルコニアボールとともにマイクロチューブミキサーで粉砕し、東洋紡Plantgenomeを用いてゲノムDNAを抽出した。得られたDNAは滅菌水に溶解した。標的領域の増幅得られたユーカリ属のゲノムDNAを鋳型としてPCR増幅を行った。PCRはゲノムDNA50ng、プライマーNS0NS9ETS4ETS2各終濃度03μM、KODPlusver2東洋紡1Unit、dATP、dTTP、dGTP、dCTP各200μM、MgCl215mMを加えた緩衝液50μl中で行った。反応条件はサーマルサイクラーDice、宝酒造社製中で、鋳型DNAを94°Cで2分間変性させ、その後、変性:94°Cで10秒間、アニーリング:58°Cで30秒間、伸長反応:68°Cで1分間を1サイクルとして30サイクル繰り返した。増幅及び配列決定に使用したプライマーは、次の表4及び表5に示したとおりである。上記の表4及び表5に示したプライマーの塩基配列は、配列番号22~24のDNAの塩基配列に基づいて設計された。ここで、配列番号22は、Eucalyptuslehmanniiの18SrRNA遺伝子の部分配列登録番号AM235528を示し、配列番号23は、EucalyptuscamaldulensisのexternaltranscribedspacerETSと18SrRNA遺伝子の部分配列登録番号DQ352528を示し、配列番号24は、Eucalyptusperrinianaの葉緑体externaltranscribedspacerETSの部分配列登録番号AM489907を示す。PCR産物の精製得られたPCR産物をQIAquickPCRpurificationkitQIAGENを用いて精製した。配列決定反応精製したサンプルを塩基長に応じて10~20ng、プライマーNS0NS2NS3NS4NS5NS6NS9、ETS2ETS4各終濃度016μM、BigDyeTerminatorv31CycleSequencingKitABIを用いてシークエンス反応を行った。反応条件はサーマルサイクラーDice、宝酒造社製中で、鋳型DNAを96°Cで1分間変性させたのち、変性:96°Cで10秒間、アニーリング:50°Cで5秒間、伸長反応:60°Cで1分間を1サイクルとして25サイクル繰り返した。精製得られたサンプルをBigDyeXterminatorABIを用いて精製した。配列決定GeneticAnalyzer3130xlABIを用いて配列決定を行った。解析条件はBDxFASTseqPOP7BDv3を用いた。得られたデータは変異検出ソフトSeqScapeABIを用いて解析した。結果それぞれ2領域を解析した結果、樹種間で差のあるSNPが7つ見つかった上記表1。これらのSNPでユーカリ属9樹種を5つのグループに分類できる可能性が高い表6。上記表1及び表6中、1~8の数字はそれぞれSNP1~SNP8を表す。また、例えばAAGという表記は、AA又はAGを意味する。また、同様に、GAGAという表記は、GG、GA、AG又はAAのいずれかを意味する。なお、領域1及び領域2の括弧内のNSは18SrRNA領域、ETSは葉緑体externaltranscribedspacerをそれぞれ表す。SNP1~SNP8のユーカリゲノム上の位置について、SNP1とSNP2はユーカリゲノムDNAの領域1に、SNP3~SNP8はユーカリゲノムDNAの領域2にそれぞれ由来し、ユーカリグランディスのゲノム配列配列番号30に基づくSNP1~SNP8の位置は、上で説明したとおりである。また、SNP1~SNP8の各塩基をPCRによって特定するために、SNP1~SNP8を含む領域を増幅するプライマー、及び増幅産物を配列決定するためのプライマーを次のように決定した。SNP1NS11NS13又はNS11NS2又はNS11NS6で増幅し、NS11及びNS2を用いて配列決定する。SNP2NS11NS13又はNS3NS4又はNS11NS6で増幅し、NS3及びNS4を用いて配列決定する。SNP3~SNP8:ETS4ETS2又はETS18SETSMyrtFで増幅し、ETS2ETS4、あるいは、ETS18SETSMyrtFで配列決定する。[実施例2]ユーカリ属樹種の識別例供試材料ユーカリカマルドレンシスと考えられるユーカリDH2の葉ユーカリカマルドレンシス×ユーカリデグルプタと考えられるユーカリK7の葉方法1DNA抽出:乾燥した葉約20mgをMagExtractorPlantGenomeTOYOBO用いて抽出した。2PCR:KODplusver2TOYOBOを用いて対象領域を増幅した。対象領域はNSETSの2領域でそれぞれに対応したプライマーNS11NS13もしくはNS11NS2NS3NS4ETS4ETS2を用いた。PCR条件は変性94°C10secアニール60°C30sec伸長68°C3060sec30cycleとした。PCR後、マルチスクリーンμ96Milliporeを用いて精製した。3配列決定:GeneticAnalyzer3130xlABIを用いて配列解析を行った。配列決定プライマーはそれぞれSNP1NS11NS2SNP2NS3NS4SNP3SNP8ETS4ETS2を用いた。結果結果を表7樹種unkownDH2及びunkownK7に示す。DH2は、ユーカリカマルドレンシスである可能性が高いことが確認できた。また、K7は、ユーカリカマルドレンシス×ユーカリデグルプタの雑種である可能性が高いことが確認できた。[実施例3]ユーカリ属の樹種間で差のあるSNPの探索この実施例では、実施例1に記載の手法と同様の手法が使用された。ユーカリ属植物9樹種ユーカリグロブラス、ユーカリカマルドレンシス、ユーカリブラシアーナ、ユーカリテレティコルニスEtereticornis、ユーカリグランディス、ユーカリサリグナEsaligna、ユーカリデグルプタ、ユーカリペリータ、及びユーカリユーロフィラについて、それらの核ゲノムDNA中のCCRCinnamoylCoAReductaseをコードする遺伝子領域の塩基配列を比較した。その結果、CCR遺伝子領域に9ヶ所のSNPが存在することを見出した。これら9つのSNPのそれぞれを、該遺伝子の5側からSNP9SNP10SNP11SNP12SNP13SNP14SNP15SNP16SNP17と称した。なお、例えばユーカリグロブラスのCCR遺伝子領域の配列は、例えばGenBankのデータベースから入手可能であり、AY656821の登録番号が付与されている。この配列配列番号31を利用して、配列が不明のユーカリ属植物のCCR遺伝子領域の増幅をPCR法で行い、該遺伝子領域を配列決定した。各樹種、少なくとも3個体以上についてCCR遺伝子領域の配列決定を行い、樹種内で共通している配列を決定した。さらに樹種間で変異があるSNPを決定した。なお、各個体、各樹種の塩基配列の比較はGENETYXVer905株式会社ゼネティックスを用いて行った。上記のSNP9~SNP17の1又は2つ以上のSNPを含む領域をPCRによって増幅するためのプライマーは次のとおりである。プライマーの位置関係は図1に示されている。また、増幅用プライマー及び塩基配列配列番号は、P220配列番号22、CCR02配列番号23、CCR001配列番号24、P225配列番号25、P223配列番号26、P221配列番号27、P219配列番号28、P217配列番号29である。これらの配列、及びプライマーの組合せフォワードとリバースは、表8に示すとおりである。このような組合せによってSNP9~SNP17のすべてを増幅することができる。増幅産物は、配列決定プライマーを用いて配列決定され、各SNPの塩基が特定された。その結果は、表2上記に示されたとおりである。SNP9~SNP17によって、ユーカリ属9樹種を4つのグループに分けることができた表9。さらに実施例1で決定されたSNP1~SNP8と、本実施例で決定されたSNP9~SNP17とを組合せて雑種の各SNPのアレル塩基を対比することによって、上記9樹種を8つのグループに分けることが可能であった表9。[実施例4]雑種の識別例EcamxEgloの識別例実施例2と同様に領域1NSのSNPと領域2ETSのSNPを利用することでグループAグロブラスとグループBの雑種であることが判るが、樹種の特定は出来なかった。しかし、実施例3と同様にCCRこの場合、CCR4、5、6、7のSNPを利用することで、グロブラスとカマルドレンシスの雑種であることが特定できた表10,11。EgraxEcamの識別例実施例2と同様に領域1のSNPと領域2のSNPを利用することでグループBとグループCの雑種であることが判る。しかし実施例3と同様にCCRこの場合、CCR1、4、5、6、7、8のSNPを利用することで、グランディスとカマルドレンシスの雑種であることが特定できた表10,11。"