説明

乳房炎耐性

本発明は、ウシ被験体における乳房炎耐性を判断する方法に関する。ここで、乳房炎耐性には無症候性乳房炎及び臨床型乳房炎の両方に対する耐性が含まれる。特に、本発明の方法は、ウシ被験体における乳房炎耐性を判断するための遺伝的マーカー及び/又は量的形質遺伝子座(QTL)の同定を含む。乳房炎耐性の決定は、特定のマイクロサテライト状態の解消を含む。さらに本発明は、乳房炎耐性に関連した遺伝的マーカーの検出用の診断キットに関する。本発明の方法及びキットは、育種目的でウシ被験体の選抜に利用され得る。したがって本発明は、より乳房炎になりにくい雌ウシを産生することとなる乳房炎耐性を有するウシ被験体を遺伝的に選抜する方法を提供する。

【発明の詳細な説明】
【技術分野】
【0001】
本発明は、ウシ被験体における乳房炎に対する耐性を判断する方法であって、ウシ染色体BTA9及びBTA11上に位置する少なくとも1つの遺伝的マーカーを検出することを含む、ウシ被験体における乳房炎に対する耐性を判断する方法に関する。さらに、本発明は、乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーの有無を検出する診断キットに関する。
【背景技術】
【0002】
乳房炎は、感染又は外傷から生じる雌ウシ(cow)の乳腺(mammary gland)又は乳房の炎症であり、乳房炎はウシにおける最も経済的に重要な疾患であると考えられている。本疾患は多様な物質によって生じるおそれがある。乳房炎の主な原因は、乳頭先端を介した微生物の乳腺侵入である。
【0003】
乳房炎は臨床型又は無症候性(sub-clinical)である可能性があり、無症候性の感染が臨床型の症状発現に先立つ。臨床型乳房炎(clinical mastitis)(CM)は、赤く且つ腫大した(swollen)乳腺、即ち、赤く腫大した乳房、凝結乳(clotted milk)の生産を通じて目視で検出され得る。一度検出されると、乳房炎の雌ウシ由来の乳は、全体の乳品質に影響を及ぼさないように、タンクから分離して保管される。
【0004】
無症候性乳房炎は、高い体細胞数(somatic cell counts)(SCS)、正常な又は上昇した体温、及び培養物の検査が陽性を示す乳試料を特徴とする乳房炎の一種である。したがって、無症候性乳房炎は目視で、乳房の腫大、又は腺若しくは産出される乳の観察によって検出することはできない。このため、農家には、無症候性乳房炎の雌ウシ由来の乳を他へそらすという選択肢がない。しかしながら、この乳は、非感染の雌ウシの乳よりも品質が劣っているため、タンク中の乳の残りを汚染するおそれがある。
【0005】
乳房炎は、体細胞数(SCS)を用いて検出され得る。ここでは、雌ウシ由来の乳の試料が体細胞(白血球)の存在について分析される。体細胞は、雌ウシの天然の防御機構の一部であり、乳房が感染すると細胞数は増加する。乳試料中の体細胞の数は、転球式粘度計(rolling-ball viscometer)及びクールター計数器によって間接的に推定され得る。
【0006】
乳房炎により、乳及び乳製品の量及び質が低下するので、農家及び乳業に経済的損失が生じる。したがって、ウシにおける乳房炎に対する耐性の遺伝的根拠を判断することが可能なことは、毎日の乳生産の観点だけではなく、乳房炎に対する耐性を有するウシ被験体を選抜する育種管理の観点においても、共に、乳業に対する非常に大きな経済的重要性がある。より乳房炎になりにくい雌ウシを産生する、耐性が改善したウシ被験体を遺伝的に選抜する方法が望ましい。
【0007】
遺伝的形質(genetic traits)に対する遺伝的決定因子(genetic determinants)を同定するための1つのアプローチは、過去の組換えの活用を目的とした連鎖不平衡(linkage disequilibrium)(LD)マッピングの使用であり、ヒト集団と比較して非常に長い染色体区域に渡って広がることが乳牛を含む幾つかの家畜(livestock)集団において示されている、(Farnir et al., 2000)。一度マッピングされると、量的形質遺伝子座(Quantitative Trait Locus)(QTL)は、マーカー利用選抜(marker assisted selection)において有益に利用され得る。
【0008】
連鎖不平衡
連鎖不均衡は、過去を遡る組換え事象(recombination events)を反映し、家系(families:ファミリー)内でのLDマッピングの使用はマッピングの分解能を増大させる。LDが存在するのは、集団中に観察されるハプロタイプが、各ハプロタイプにおける個々の遺伝的マーカーの頻度を相互に掛けることによって予測されるハプロタイプ頻度と一致しない場合である。この点で、ハプロタイプという用語は、相互に遺伝する傾向がある一染色体上に存在する一組の密接に連鎖した遺伝的マーカーを意味する。LDマッピングが効率的であるためには、遺伝的マーカーの密度が、所定の集団内でLDが広がる距離と調和する必要がある。多数の遺伝的マーカー(284個の常染色体マイクロサテライトマーカー)を使用する乳牛集団におけるLDの研究において、LDは、染色体内(intrachromosomal)マーカーについて数十センチモルガン(centimorgans)に渡って広がることが実証された(Farnir et al., 2000)。同様に、Georges, M(2000)の報告によれば、家畜における特定の表現型(phenotype)と連鎖する遺伝的マーカーの位置の信頼区間(confidence interval)は、典型的には20〜30cMである(おそらく500〜1000遺伝子に対応する)(Georges, M., 2000)。高い信頼度をもって興味のある特定の領域のマップを使用するために、連鎖不均衡の存在が考慮される。
【0009】
本発明において、臨床型乳房炎及び/又はSCSに関連した量的形質遺伝子座はウシ染色体BTA9で同定されており、ウシ被験体が乳房炎に対して耐性があるか否かを判断する方法を可能にする。
【発明の開示】
【発明が解決しようとする課題】
【0010】
乳房炎に対する耐性が増大したウシ被験体を選抜し、それにより乳房炎を患う動物に関連した経済的損失を避けることが可能になることは、ウシ産業内で経済的関心が大きい。乳房炎に対する耐性に関する遺伝的素因(genetic predisposition)は、本発明によって検出され得る。本発明は、乳房炎に対する耐性に関連し且つ/又はそれと連鎖する遺伝的マーカーに基づき、ウシ被験体において乳房炎に対する耐性を判断する方法を提供する。
【課題を解決するための手段】
【0011】
したがって、本発明の一態様は、ウシ被験体における乳房炎に対する耐性を判断するための方法であって、当該ウシ被験体からの試料中に、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を検出することを含み、当該少なくとも1つの遺伝的マーカーが、ウシ染色体BTA9上で、多型マイクロサテライトマーカーC6orf93及びinra084に隣接する領域、及び多型マイクロサテライトマーカーC6orf93及びinra084を含む領域、並びに/又はウシ染色体BTA11上で、多型マイクロサテライトマーカーHELMTT43及びBM3501に隣接する領域、及び多型マイクロサテライトマーカーHELMTT43及びBM3501を含む領域に位置し、当該少なくとも1つの遺伝的マーカーの有無が当該ウシ被験体又はその出生児の乳房炎耐性を示す方法に関する。
【0012】
本発明の第二の態様は、ウシ被験体において乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーの有無を検出する際に使用するための診断キットであって、少なくとも1つのオリゴヌクレオチド配列及びそれらの組合せを含み、当該ヌクレオチド配列が配列番号1〜配列番号192のいずれか及び/又はそれらの任意の組合せから選択される、診断キットに関する。
【発明を実施するための最良の形態】
【0013】
本発明は、乳牛における乳房炎耐性の遺伝的決定因子に関する。乳房炎、臨床型乳房炎及び無症候性乳房炎の発生は共に、乳業に対する実質的な経済的損失に関与する。したがって、乳房炎耐性に関する遺伝的素因を有するそれらのウシ被験体を同定することには、経済的な関心がある。かかる遺伝的素因を有するウシ被験体は、出生児に伝播し得る望ましい形質のキャリア(carriers)である。
【0014】
「ウシ被験体」という用語は任意の品種のウシを指し、成体動物であれ新生動物であれ、雌ウシ(cows)及び雄ウシ(bulls)の両方を含む意味である。この用語によって特定の年齢の動物が示されることはない。ウシ被験体の一例は、ホルスタイン種のひとつである。好ましい一実施形態において、ウシ被験体は、ホルスタインフリーシャンウシ集団のひとつである。別の実施形態において、ウシ被験体は、ホルスタインスワルトボント(Swartbont)ウシ集団のひとつである。別の実施形態において、ウシ被験体は、ドイツホルスタインシュワルツブントウシ集団のひとつである。別の実施形態において、ウシ被験体は、米国ホルスタインウシ集団のひとつである。一実施形態において、ウシ被験体は、レッドアンドホワイトホルスタイン種のひとつである。別の実施形態において、ウシ被験体は、ドイツホルスタインシュワルツブントウシ集団のひとつである。一実施形態において、ウシ被験体は、ホルスタイン種のひとつを含む、任意の家系のひとつである。好ましい一実施形態において、ウシ被験体は、デンマークレッド集団のひとつである。別の好ましい実施形態において、ウシ被験体は、フィンランドエアシャー集団のひとつである。さらに別の実施形態において、ウシ被験体は、スウェーデンレッドアンドホワイト(Swedish Red and White)集団のひとつである。さらなる実施形態において、ウシ被験体は、デンマークホルスタイン集団のひとつである。別の実施形態において、ウシ被験体は、スウェーデンレッドアンドホワイト集団のひとつである。さらに別の実施形態において、ウシ被験体は、ノルディックレッド(Nordic Red)集団のひとつである。
【0015】
本発明の一実施形態において、ウシ被験体は、スウェーデンレッドアンドホワイト、デンマークレッド、フィンランドエアシャー、ホルスタインフリーシャン、デンマークホルスタイン及びノルディックレッドから成る群から選択される。本発明の別の実施形態において、ウシ被験体は、フィンランドエアシャー及びスウェーデンレッドアンドホワイトウシから成る群から選択される。本発明の別の実施形態において、ウシ被験体は、フィンランドエアシャー及びスウェーデンレッドアンドホワイトウシから成る群から選択される。
【0016】
一実施形態において、ウシ被験体は表1aに示される品種の群から選択される。
【0017】
【表1a】

【0018】
一実施形態において、ウシ被験体は表1bに示される品種の群から選択される品種のひとつである。
【0019】
【表1b】

【0020】
一実施形態において、ウシ被験体は表1cに示される品種の群から選択される品種のひとつである。
【0021】
【表1c】

【0022】
「遺伝的マーカー」という用語は、ウシ染色体上のDNAの可変ヌクレオチド配列(variable nucleotide sequence)(多型(polymorphism))を指し、アレル同士を判別する。可変ヌクレオチド配列は、例えば増幅法(amplification methods)及び/又はサイズの差異の観察において特定のオリゴヌクレオチドを用いることにより、当業者に既知の方法によって同定することができる。しかしながら、可変ヌクレオチド配列は、シーケンシングにより、又は例えば制限断片長多型(restriction fragment length polymorphism)解析によっても検出され得る。可変ヌクレオチド配列は、欠失、挿入、反復及び/又は点変異(point mutation)によって表され得る。
【0023】
遺伝的マーカーの1つの型は、量的形質遺伝子座と連鎖するマイクロサテライトマーカーである。マイクロサテライトマーカーは、互いの後で反復される短い配列を指す。短い配列は、例えば1ヌクレオチド、例えば2ヌクレオチド、例えば3ヌクレオチド、例えば4ヌクレオチド、例えば5ヌクレオチド、例えば6ヌクレオチド、例えば7ヌクレオチド、例えば8ヌクレオチド、例えば9ヌクレオチド、例えば10ヌクレオチドである。しかしながら、変化が生じることがあり、反復の数は増加又は減少し得る。多型マイクロサテライトマーカーの特定の定義及び遺伝子座は、USDA遺伝子地図(genetic map)に見出され得る(Kappes et al., 1997;又は米国食肉動物研究所(http://www.marc.usda.gov/)へのリンクをたどることによる)。
【0024】
さらに、本発明の遺伝的マーカーのヌクレオチド配列は、ウシ被験体における乳房炎耐性に関する形質と遺伝的に連鎖する(genetically linked)と理解される。結果的に、本発明の方法により、遺伝的マーカーに隣接する領域、及び遺伝的マーカーを含むDNA領域のヌクレオチド配列から、多数の遺伝的マーカーが生成され得ることも理解される。
【0025】
「量的形質遺伝子座(QTL)」という用語は、特定形質(例えば、草丈)に関連するDNA領域である。必ずしもそれ自体が遺伝子である必要はないが、QTLは、問題の形質の基となる遺伝子に密接に連鎖する一続きの(stretches of)DNAである。
【0026】
「乳房炎」という用語は、本願において例えば高い体細胞スコア(SCS)を特徴とする無症候性乳房炎、及び臨床型乳房炎の両方を示すのに用いられる。
【0027】
「乳房炎耐性」及び「乳房炎に対する耐性」という用語は交換して用いることが可能であり、或るウシ被験体が他のウシ被験体より乳房炎になりにくいということに関係する。乳房炎に対する耐性に関連する遺伝的マーカーを判断するため本発明の如く多数のウシ被験体の解析を行う場合には、解析したウシ被験体における臨床型乳房炎及び/又は無症候性乳房炎の発生と連鎖した遺伝的マーカーの有無によって、乳房炎に対する耐性をほのめかす形質が観察され得る。乳房炎耐性には、ウシ被験体又はその出生児における乳房の健康に影響する形質に対する耐性、を包含することが理解される。したがって、雄ウシの乳房炎耐性は、その雌出生児によって身体的に明らかにされる。
【0028】
乳房炎耐性についてのスコア化
雄ウシの娘を、乳房炎耐性及びSCCについてスコアリングした。体細胞スコア(somatic cell score, SCS)は、乳記録スキーム(milk recording scheme)から得られるlog10変換した体細胞数値の平均として(10,000/mLで)定義した。この平均は出産後10日〜180日の期間でとった。息子の形質に対する推定育種価(estimated breeding values, EBV)を、家系構成を無視する単一形質の最良線形不偏予測(Best Linear Unbiased Prediction, BLUP)動物モデルを使用して算出した。これらのEBVをQTL解析に用いた。個々の形質で用いた娘の登記項目(registrations)は以下であった:
デンマークにおける臨床型乳房炎:第一出産前5日〜出産後50日の期間での臨床型乳房炎の治療症例、
スウェーデン及びフィンランドにおける臨床型乳房炎:第一出産前7日〜出産後150日の期間での臨床型乳房炎の治療症例、
デンマークにおけるSCS:第一出産後10日〜180日の期間での平均SCS、
スウェーデンにおけるSCS:第一出産後10日〜150日の期間での平均SCS、
フィンランドにおけるSCS:第一出産後10日〜305日の期間での平均SCS。
【0029】
本発明の一実施形態において、本明細書中に記載の方法及びキットは、乳房炎耐性に関する。本発明の別の実施形態において、本明細書中に記載の方法及びキットは、臨床型乳房炎に対する耐性に関する。別の実施形態において、本発明の方法及びキットは、体細胞数等によって検出される無症候性乳房炎に対する耐性に関する。さらに別の実施形態において、本発明の方法及びキットは主に、例えば、体細胞数によって検出される無症候性乳房炎に対する耐性と組み合わされた臨床型乳房炎に対する耐性、に関する。
【0030】
試料
本発明による方法はウシ被験体の試料を分析することを含む。ここで、当該試料は、当該方法において使用するためのウシ遺伝子材料を提供可能な任意の好適な試料であり得る。ウシ遺伝子材料は、例えば、血液試料、組織試料(例えば脾臓(spleen)、頬側塗抹(buccal smears))、体表(body surface)(毛又は爪)の切取り部、乳及び/又は精液(semen)から抽出、単離され、必要に応じて精製され得る。試料は新鮮でもよく、凍結してもよい。
【0031】
本発明の配列多型は、少なくとも1つのヌクレオチド差異、例えば少なくとも2つのヌクレオチド差異、例えば少なくとも3つのヌクレオチド差異、例えば少なくとも4つのヌクレオチド差異、例えば少なくとも5つのヌクレオチド差異、例えば少なくとも6つのヌクレオチド差異、例えば少なくとも7つのヌクレオチド差異、例えば少なくとも8つのヌクレオチド差異、例えば少なくとも9つのヌクレオチド差異、例えば10個のヌクレオチド差異を含む。ヌクレオチド差異は、ヌクレオチド差異、欠失及び/又は挿入、或いはそれらの任意の組合せを含む。
【0032】
グランドドーターデザイン
グランドドーターデザイン(grand daughter design)は、育種において広範に使用されている祖父についての、及び祖父の息子が出生児を出産した場合には当該息子についての、DNAベースのマーカーからデータを分析することを含む。DNAマーカーデータと共に使用されることとなる表現型データは、当該息子の娘由来である。かかる表現型データは、例えば乳生産の特徴、出産に関する特徴、肉質又は疾患であり得る。娘の一群は、その父親から一方のアレルを受け継ぐ一方、別の群の娘は、その父親から他方のアレルを受け継いでいる。2つの群からのデータを比較することにより、特定の染色体の断片が問題の形質に影響を及ぼす1つ又は複数の遺伝子を保有しているかどうかについての情報を得ることができる。QTLが染色体のこの断片内に存在するかどうかを結論付けることができる。
【0033】
グランドドーターデザインを実行するための前提条件は、詳細な表現型データが利用可能であることである。本発明において、かかるデータが入手可能となっている(http://www.lr.dk/kvaeg/diverse/principles.pdf)。
【0034】
QTLは、量的形質遺伝子座の短縮形である。個体に量的形質(quantitative traits)を与える遺伝子は、あたかも1つ又は複数の遺伝子が位置するように作用(act)する染色体の一片を観察することによって、間接的に見出され得る。
【0035】
対照的に、ある染色体上の多数のDNAマーカーが一方又は両方の親及びそれらの出生児について決定された場合、親からその出生児のうちの1頭又は複数頭へと伝えられた形質の情報を提供するために、DNAマーカーが直接使用され得る。このマーカーを使用して、DNAマーカーと連鎖した染色体の遺伝的履歴(genetic history)を算出し得る。
【0036】
染色体領域及びマーカー
BTAは、ウシ常染色体(Bos Taurus autosome)の短縮形である。本発明の一態様は、ウシ被験体における乳房炎に対する耐性を判断するための方法であって、当該ウシ被験体からの試料中で、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を検出することを含み、当該少なくとも1つの遺伝的マーカーが、ウシ染色体BTA9上で、多型マイクロサテライトマーカーC6orf93及びinra084に隣接する領域、及び多型マイクロサテライトマーカーC6orf93及びinra084を含む領域、並びに/又はウシ染色体BTA11上で、多型マイクロサテライトマーカーHELMTT43及びBM3501に隣接する領域、及び多型マイクロサテライトマーカーHELMTT43及びBM3501を含む領域に位置し、当該少なくとも1つの遺伝的マーカーの有無が当該ウシ被験体又はその出生児の乳房炎耐性を示す、方法に関する。本明細書に記載されるような連鎖不平衡の概念のために、本発明は、ウシ被験体において乳房炎に対する耐性を判断することにも関する。ここで、少なくとも1つの遺伝的マーカーが、乳房炎に対する耐性についてのウシの形質と連鎖する。ウシ被験体において乳房炎に対する耐性を判断するために、一つより多くの遺伝的マーカーが本発明において利用され得ると理解される。例えば少なくとも1つの遺伝的マーカーとは、精度が増大し得るように、少なくとも2つ以上の遺伝的マーカーの組合せであり得て、例えば少なくとも3つの遺伝的マーカー、例えば4つの遺伝的マーカー、例えば少なくとも5つの遺伝的マーカー、例えば6つの遺伝的マーカー、例えば少なくとも7つの遺伝的マーカー、例えば8つの遺伝的マーカー、例えば少なくとも9つの遺伝的マーカー、例えば10個の遺伝的マーカーである。
【0037】
少なくとも1つの遺伝的マーカーは、少なくとも1つのウシ染色体、例えば2つの染色体、例えば3つの染色体、例えば4つの染色体、例えば5つの染色体、及び/又は例えば6つの染色体上に位置し得る。少なくとも1つの遺伝的マーカーがウシ染色体9上に位置し得る。しかしながら、当該少なくとも1つの遺伝的マーカーは、異なる染色体上に位置するマーカーの組合せでもよい。
【0038】
少なくとも1つのマーカーは、本明細書中に示される表の個々のマーカーのいずれかから選択される。
【0039】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーc6orf93及びrgs17に隣接する領域、並びにマーカーc6orf93及びrgs17を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約69.35cM〜約79.8cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表1dで示されるマーカーの群から選択される。
【0040】
【表1d】

【0041】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上に位置する。一実施形態では、ウシ染色体BTA9上で、少なくとも1つの遺伝的マーカーは、マーカーc6orf93及びinra084に隣接する領域、並びにマーカーc6orf93及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約69.35cM〜約74.5cM(この分析に用いられる位置に従って)の領域に位置する。MARCマーカー地図に従って、遺伝的マーカーinra084の位置は90.98である。少なくとも1つの遺伝的マーカーは、表2で示されるマーカーの群から選択される。
【0042】
【表2】

【0043】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2251及びinra084に隣接する領域、並びにマーカーbms2251及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約71.3cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約86.58cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表3で示されるマーカーの群から選択される。
【0044】
【表3】

【0045】
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2251及びinra144に隣接する領域、並びにマーカーbms2251及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約71.3cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約86.58cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表4で示されるマーカーの群から選択される。
【0046】
【表4】

【0047】
さらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びinra084に隣接する領域、並びにマーカーbm7234及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表5で示されるマーカーの群から選択される。
【0048】
【表5】

【0049】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びinra144に隣接する領域、並びにマーカーbm7234及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表6で示されるマーカーの群から選択される。
【0050】
【表6】

【0051】
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びbms2819に隣接する領域、並びにマーカーbm7234及びbms2819を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約73.95cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表7で示されるマーカーの群から選択される。
【0052】
【表7】


【0053】
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びbm4208に隣接する領域、並びにマーカーbm7234及びbm4208を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約73.9cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.69cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表8で示されるマーカーの群から選択される。
【0054】
【表8】

【0055】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2819及びinra144に隣接する領域、並びにマーカーbms2819及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約73.95cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.98cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表9で示されるマーカーの群から選択される。
【0056】
【表9】

【0057】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2819及びinra084に隣接する領域、並びにマーカーbms2819及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約73.95cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.98cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表10で示されるマーカーの群から選択される。
【0058】
【表10】


【0059】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm4208及びinra144に隣接する領域、並びにマーカーbm4208及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約73.9cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.69cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表11で示されるマーカーの群から選択される。
【0060】
【表11】

【0061】
本発明のさらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーinra144及びinra084に隣接する領域、並びにマーカーinra144及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約74.2cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.98cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表12で示されるマーカーの群から選択される。
【0062】
【表12】


【0063】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2251及びbm7234に隣接する領域、並びにマーカーbms2251及びbm7234を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約71.3cM〜約72.3cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約86.58cM〜約88.136cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表13で示されるマーカーの群から選択される。
【0064】
【表13】


【0065】
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーEPM2A及びbm7234に隣接する領域、並びにマーカーEPM2A及びbm7234を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.1cM〜約72.3cM(この分析に用いられる位置に従って)の領域、及びMARCマーカー地図に従ってbm7234については約88.136cMに位置する。少なくとも1つの遺伝的マーカーは、表14で示されるマーカーの群から選択される。
【0066】
【表14】


【0067】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA9上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーinra144及びrgs17に隣接する領域、並びにマーカーinra144及びrgs17を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約74.2cM〜約79.8cM(この分析に用いられる位置に従って)の領域に位置し、MARCマーカー地図に従ったinra144の位置は90.98cMである。少なくとも1つの遺伝的マーカーは、表15で示されるマーカーの群から選択される。
【0068】
【表15】

【0069】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA9上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーinra084及びrgs17に隣接する領域、並びにマーカーinra084及びrgs17を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約74.5cM〜約79.8cM(この分析に用いられる位置に従って)の領域に位置し、MARCマーカー地図に従ったinra084の位置は90.68cMである。少なくとも1つの遺伝的マーカーは、表16で示されるマーカーの群から選択される。
【0070】
【表16】

【0071】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びBM3501に隣接する領域、並びにマーカーHELMTT43及びBM3501を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM〜約97.223cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表17で示されるマーカーの群から選択される。
【0072】
【表17−1】

【0073】
【表17−2】

【0074】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びINRA177に隣接する領域、並びにマーカーHELMTT43及びINRA177を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM〜約35.098cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表18で示されるマーカーの群から選択される。
【0075】
【表18】


【0076】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びMNB−70に隣接する領域、並びにマーカーHELMTT43及びMNB−70を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM〜約24.617cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表19で示されるマーカーの群から選択される。
【0077】
【表19】

【0078】
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB−40及びMNB−70に隣接する領域、並びにマーカーMNB−40及びMNB−70を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約24.617cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表20で示されるマーカーの群から選択される。
【0079】
【表20】


【0080】
さらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBP38及びINRA131に隣接する領域、並びにマーカーBP38及びINRA131を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約24.617cM〜約47.289cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表21で示されるマーカーの群から選択される。
【0081】
【表21】


【0082】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM2818及びINRA177に隣接する領域、並びにマーカーBM2818及びINRA177を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約30.009cM〜約35.098cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表22で示されるマーカーの群から選択される。
【0083】
【表22】


【0084】
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBMS1953及びBM2818に隣接する領域、並びにマーカーBMS1953及びBM2818を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約21.537cM〜約30.009cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表23で示されるマーカーの群から選択される。
【0085】
【表23】

【0086】
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びZAP70に隣接する領域、並びにマーカーHELMTT43及びZAP70を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM(MARCマーカー地図に従って)〜約5.4cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表24で示されるマーカーの群から選択される。
【0087】
【表24】

【0088】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーZAP70及びIL18RAに隣接する領域、並びにマーカーZAP70及びIL18RAを含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、(この分析に用いられる位置に従って)ウシ染色体BTA11上の約5.4cM〜約12.3cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表25で示されるマーカーの群から選択される。
【0089】
【表25】


【0090】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーINRA131及びBM6445に隣接する領域、並びにマーカーINRA131及びBM6445を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約47.289cM〜約61.570cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表26で示されるマーカーの群から選択される。
【0091】
【表26】

【0092】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM304及びBM7169に隣接する領域、並びにマーカーBM304及びBM7169を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約33.597cM〜約50.312cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表27で示されるマーカーの群から選択される。
【0093】
【表27】

【0094】
本発明のさらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM7169及びDIK5170に隣接する領域、並びにマーカーBM7169及びDIK5170を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約50.312cM〜約70.143cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表28で示されるマーカーの群から選択される。
【0095】
【表28】


【0096】
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM6445及びBMS1048に隣接する領域、並びにマーカーBM6445及びBMS1048を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約61.570cM〜約81.065cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表29で示されるマーカーの群から選択される。
【0097】
【表29】


【0098】
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMB110及びBMS2047に隣接する領域、並びにマーカーMB110及びBMS2047を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約68.679cM〜約78.457cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表30で示されるマーカーの群から選択される。
【0099】
【表30】


【0100】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びBM2818に隣接する領域、並びにマーカーIL18RA及びBM2818を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約30.009cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表31で示されるマーカーの群から選択される。
【0101】
【表31】

【0102】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM2818及びBM7169に隣接する領域、並びにマーカーBM2818及びBM7169を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約30.009cM〜約50.312cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表32で示されるマーカーの群から選択される。
【0103】
【表32】


【0104】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMAP4K4及びBM2818に隣接する領域、並びにマーカーMAP4K4及びBM2818を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約10.5cM(この分析に用いられる位置に従って)〜約30.009cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表33で示されるマーカーの群から選択される。
【0105】
【表33】


【0106】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びUMBTL103に隣接する領域、並びにマーカーIL18RA及びUMBTL103を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約23.829cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表34で示されるマーカーの群から選択される。
【0107】
【表34】

【0108】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB−40及びDIK2653に隣接する領域、並びにマーカーMNB−40及びDIK2653を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約20.135cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表35で示されるマーカーの群から選択される。
【0109】
【表35】

【0110】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM716及びDIK4637に隣接する領域、並びにマーカーBM716及びDIK4637を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約22.527cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表36で示されるマーカーの群から選択される。
【0111】
【表36】

【0112】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM716及びBMS2569に隣接する領域、並びにマーカーBM716及びBMS2569を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約21.082cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表37で示されるマーカーの群から選択される。
【0113】
【表37】


【0114】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBMS2325及びDIK4637に隣接する領域、並びにマーカーBMS2325及びDIK4637を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約21.082cM〜約22.527cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表38で示されるマーカーの群から選択される。
【0115】
【表38】


【0116】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びAUP1に隣接する領域、並びにマーカーIL18RA及びAUP1を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約17.6cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表39で示されるマーカーの群から選択される。
【0117】
【表39】

【0118】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びMNB−40に隣接する領域、並びにマーカーIL18RA及びMNB−40を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約19.440cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表40で示されるマーカーの群から選択される。
【0119】
【表40】

【0120】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB−40及びAUP1に隣接する領域、並びにマーカーMNB−40及びAUP1を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM(MARCマーカー地図に従って)〜約17.6cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表41で示されるマーカーの群から選択される。
【0121】
【表41】

【0122】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RAを含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表42で示される。
【0123】
【表42】

【0124】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB−40を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約16.0cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表43で示される。
【0125】
【表43】

【0126】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーAUP1を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約17.6cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表44で示される。
【0127】
【表44】

【0128】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーDIK4637及びUMBTL103に隣接する領域、並びにマーカーDIK4637及びUMBTL103を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約22.527cM〜約23.829cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表45で示されるマーカーの群から選択される。
【0129】
【表45】

【0130】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーDIK4637を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約22.527cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表46で示される。
【0131】
【表46】

【0132】
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーUMBTL103を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約23.829cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表47で示される。
【0133】
【表47】

【0134】
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、本明細書の他の箇所で記載されるように、BTA9の任意の領域及びマーカーがBTA11の任意の領域及びマーカーと組み合わされるマーカーの組合せである。
【0135】
プライマー
表50に、本発明により使用され得るプライマーを示す。表50で特定されるプライマー対は個々に使用してもよく、表50の1つ又は複数のプライマー対と併用してもよい。かかるプライマー又はプローブの設計は、当業者である分子生物学者には明らかであろう。かかるプライマーは、任意の好都合な長さ、例えば最大50塩基長、最大40塩基長、さらに好都合には最大30塩基長、例えば8〜25塩基長又は8〜15塩基長を有する。概して、かかるプライマーは、領域内の対応する野生型又はバリアント遺伝子座(variant locus)と完全に相補的な塩基配列(base sequences)を含むであろう。しかしながら、必要に応じて、過度にオリゴヌクレオチドプローブの判別力(discriminatory power)に影響を及ぼさないという条件で、1つ又は複数ののミスマッチ(mismatches)を導入してもよい。本発明のプライマー/プローブは、検出を容易にするための1つ又は複数の標識を保有し得る。
【0136】
一実施形態において、プライマー及び/又はプローブは、本明細書中に示されるようなマーカーによって描写される配列を含有する一塩基多型(single nucleotide polymorphism)とハイブリッド形成するサブシーケンス(subsequence)とのハイブリダイゼーション、及び/又はその増幅が可能である。
【0137】
本発明のプライマーヌクレオチド配列はさらに以下を含む:(a)遺伝的マーカー配列又はその相補配列若しくはRNA産物を含む核酸分子とハイブリッド形成する任意のヌクレオチド配列、ここでハイブリッド形成はストリンジェント条件下(例えば約45℃での6X塩化ナトリウム/クエン酸ナトリウム(SSC)中でのフィルター結合DNAとのハイブリダイゼーションと、その後の約50〜65℃で0.2XSSC/0.1%ドデシル(Dodecyl)硫酸ナトリウム(SDS)中での1回若しくは複数回の洗浄)、又は(b)高ストリンジェント条件下(例えば約45℃での6XSSC中でのフィルター結合核酸とのハイブリダイゼーションと、その後の約68℃で0.1XSSC/0.2% SDS中での1回若しくは複数回の洗浄)、又は当業者に明らかな他のハイブリダイゼーション条件下(例えば、Ausubel F.M. et al.編「分子生物学における現行のプロトコール(Current Protocols in Molecular Biology)」第I巻、Green Publishing Associates, Inc., and John Wiley & sons, Inc., New York, at pp. 6.3.1-6.3.6 and 2.10.3(1989)を参照されたい)、で行う。好ましくは、上記(a)及び(b)のヌクレオチド配列とハイブリッド形成する核酸分子は、遺伝的マーカー配列又はその相補配列若しくはRNA産物を含むゲノムDNA(genomic DNA)の核酸分子のコンプレメント(complement)を含むものである。
【0138】
本発明の核酸分子のうちの一つは、高ストリンジェント条件又はストリンジェント条件下で上記の核酸分子とハイブリッド形成する、デオキシオリゴヌクレオチド(「オリゴ(oligos)」)である。概して、14〜70ヌクレオチド長のプローブについては、融解温度(melting temperature)(TM)は、次式:
Tm(℃)=81.5+16.6(log[一価カチオン(monovalent cations)(モーラー(molar))])+0.41(%G+C)−(500/N)
(式中、Nはプローブ長である)を用いて算出される。ハイブリダイゼーションがホルムアミドを含有する溶液中で実行される場合、融解温度は、等式Tm(℃)=81.5+16.6(log[一価カチオン(モーラー(molar))])+0.41(%G+C)−(0.61%ホルムアミド)−(500/N)(式中、Nはプローブ長である)を使用して算出される。概して、ハイブリダイゼーションは、(DNA−DNAハイブリッドについては)Tmより約20〜25℃低い温度で、(RNA−DNAハイブリッドについては)Tmより10〜15℃低い温度で実行される。
【0139】
例示的な高ストリンジェント条件は、例えば、37℃(約14塩基オリゴについて)、48℃(約17塩基オリゴについて)、55℃(約20塩基オリゴについて)、及び60℃(約23塩基オリゴについて)での、6XSSC/0.05%ピロリン酸(pyrophosphate)ナトリウム中での洗浄を指す可能性がある。
【0140】
したがって、本発明はさらに、本発明の多型を検出するヌクレオチドプライマー又はプローブを提供する。評価は、少なくとも1つの核酸プライマー又はプローブの手段、例えばDNA、RNA又はペプチド核酸(peptide nucleic acid, PNA)若しくはロックド核酸(locked nucleic acid, LNA)等の核酸アナログ(analogue)のプライマー又はプローブによって実行され得る。
【0141】
本発明の一態様によれば、描写領域中の1つ又は複数の位置で多型を検出することが可能なアレル特異的オリゴヌクレオチドプローブが提供される。
【0142】
アレル特異的オリゴヌクレオチドプローブは、好ましくは5〜50ヌクレオチド、より好ましくは約5〜35ヌクレオチド、より好ましくは約5〜30ヌクレオチド、より好ましくは少なくとも9ヌクレオチドである。
【0143】
連鎖の決定
遺伝的マーカーが遺伝子材料中に存在するかどうかを検出するために、当業者に既知の標準的方法が適用され得る(例えば核酸増幅を用いて)。遺伝的マーカーが乳房炎耐性形質に遺伝的に連鎖するかどうかを判断するために、順列検定(permutation test)を適用するか(Doerge and Churchill, 1996)、又はピエフォ法(Piepho-method)を適用することができる(Piepho, 2001)。順列検定の原理はDoerge and Churchill(1996)によって十分に記載されている一方、ピエフォ法はPiepho(2001)によって十分に記載されている。回帰法を使用する家系内解析における有意な連鎖に、1変異当たり10000の順列が、順列検定を使用して作製された(Doerge and Churchill, 1996)。5%の染色体全体水準での閾値は、遺伝的マーカーと、乳房炎耐性及び体細胞数形質との間の連鎖についての有意な証拠であると見なされた。また、QTLは異なる雄親家系で確認された。分散成分(variance component)法を使用する家系間解析及び多形質解析(multi-trait analysis)については、ピエフォ法を使用して有意水準(significance level)を求めた(Piepho, 2001)。5%の染色体全体水準での閾値は、遺伝的マーカーと、乳房炎耐性及び体細胞数形質との間の連鎖についての有意な証拠であると見なされた。
【0144】
キット
本発明の別の態様は、ウシ被験体における乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーの有無を検出する際に使用するための診断キットであって、少なくとも1つのオリゴヌクレオチド配列及びそれらの組合せを含み、ヌクレオチド配列が配列番号1〜配列番号192のいずれか及び/又はそれらの任意の組合せから選択される、診断キットに関する。
【0145】
本発明によるウシ被験体についての乳房炎に対する耐性の遺伝的決定因子(genetic determinants)を確立するためのウシ被験体の遺伝子型同定(genotyping)は、DNA及び/又はRNAの解析に基づき得る。一例としては、本明細書中に記載の標準的なDNA抽出法を使用して提供され得るゲノムDNAである。ゲノムDNAは、例えば多型マーカー(polymorphic marker)領域に対応(相補的)するオリゴヌクレオチドプライマーを使用するポリメラーゼ連鎖反応等の、標準的な技法を使用して単離且つ増幅し得る。増幅反応前にDNAを精製するさらなる工程が含まれ得る。したがって、乳房炎耐性及び体細胞数特性を確立するための診断キットは、別個の包装中に、表xxに示される配列の群から選択される少なくとも1つのオリゴヌクレオチド配列、及びそれらの任意の組合せを含む。
[実施例]
【0146】
動物
動物材料は39個の父系性雄親家系(paternal sire families)を有するグランドドーターデザインから成り、出生児を調べた息子の総数は、4つの乳牛品種、即ちデンマークホルスタイン(DH)及びデンマークレッド(DR)、フィンランドエアシャー(FA)及びスウェーデンレッドアンドホワイト(SRB)で1513頭であった。これらの39個の家系は、5つのDH、9つのDR、11個のFA、及び14個のSRBの祖父家系から成る。祖父1頭当たりの息子の数は16〜161頭の範囲であり、平均家系サイズは38.8であった。
【0147】
ゲノムDNAの精製
以下のプロトコールに従って、精液からゲノムDNAを精製した。精液ストローを解凍後、ストローの両端を鋏で切り落とし、精液内容物を1.5ml容のエッペンドルフチューブに移した。0.9% NaCl 1mlを使用して、ストローを洗い、管中に流し入れた。次にチューブを2000rpmで5分間遠心分離し、その後上清(supernatant)を除去した。この洗浄工程を2回繰り返した。次に、300μlの緩衝液S(10mM Tris-HCl(pH8)、100mM NaCl、10mM EDTA(pH8)、0.5% SDS)、1M DTT 20μl及びプロナーゼ(20mg/ml)(Boehringer)20μlを、チューブに入れる。混合後、チューブをゆっくりと回転させながら一晩インキュベートし、飽和NaCl 180μlを添加後、15秒間激しく撹拌する。チューブを11000rpmで15分間遠心分離する。上清0.4mlを2ml容チューブに移し、96%エタノール1mlを添加し、チューブをゆっくりと回転させることによって混合を達成する。次にチューブを、11000rpmで10分間遠心分離する。液体を注ぎ出すことによって上清を除去し、70%エタノール(0.2ml)を用いてペレット(pellet)を洗浄し、再度11000rpmで10分間遠心分離する。エタノールを注ぎ出し、ペレットを乾燥させて、55℃で30分間、TE緩衝液0.5ml中に再懸濁する。
【0148】
増幅手順
TEMPアーゼ(GeneChoice)ポリメラーゼ及び供給元(GeneChoice)により提供される反応緩衝液Iを使用して、8μlの容量中でPCR反応を実行した。通常、各多重PCR(multiplex PCR)中には、5つの異なるマーカーが含まれる(DNA 1μl、TEMPアーゼ酵素0.1μl、0.2mM dNTP、1.2mM MgCl2、0.3μM各プライマー)。
【0149】
PCR混合物を、(TEMPアーゼについては)94℃で15分間、初期変性に付した。その後、試料をタッチダウンしながら10サイクル、サイクル処理した。即ち、温度を各サイクル(94℃、30秒間の変性、67℃、45秒間のアニーリング、72℃、30秒間の伸長(elongation))で1℃下げて、その後、試料を普通のPCR条件(94℃、30秒間の変性、58℃、45秒間のアニーリング、72℃、30秒間の伸長)で20サイクル、サイクル処理した。72℃で30分間の1サイクルでPCRサイクル処理を終結し、PCR機械を、「ずっと(for “ever”)」4℃に試料を冷却するようにプログラムした。
【0150】
マーカーの検出に使用されるプライマーのヌクレオチド配列を表50に示す。配列を5’末端から列記する。
【0151】
【表50−1】

【0152】
【表50−2】

【0153】
【表50−3】

【0154】
【表50−4】

【0155】
【表50−5】

【0156】
【表50−6】

【0157】
マーカー及びマップ
本研究のBTA9については、45個のマイクロサテライトマーカーが食肉動物研究所のウェブサイト(www.marc.usda.gov/genome/genome.html)から選択された。BTA9はHSA6qに対してオーソロガス(orthologous)であるので、HSA6qにそって28個の公開された遺伝子及びESTが(Ctgf、Vip、Vil2、Rgs17、Ros1、Slc16a10、Oprm、igf2r、Esr1、Deadc1、Pex3、C6orf93、Ifngr1、Shprh、Epm2a、AK094944、AK094379、Utrn、Tnf、plg、arid1b、Iama4、hivep2、C6orf055、CITED2、RP1−172K10、AIG1、GRM)、SNP及びマイクロサテライトスクリーニングのために選択された。連鎖解析によって、BTA9の緻密なマップを作製するために、本研究で同定された8つの新規のマイクロサテライトマーカー及び29個のSNPが、系統に渡って(across the pedigree)遺伝子型同定された。
【0158】
BTA11の連鎖マップにおける全37個のマーカーの中で、本研究では、30個のマイクロサテライトマーカーを食肉動物研究所のウェブサイト(http://www.marc.usda.gov/genome/cattle/cattle.html)から選択した。
【0159】
ラディエーションハイブリッド(radiation hybrid, RH)パネル情報
遺伝子及びMARCマイクロサテライトを含むマイクロサテライトをマッピングするために、特異的なプライマー対をウシの配列から設計した。染色体9に沿って、総数120個を超えるマーカーをウシRHパネルに用いた。65個の遺伝子及び34個のマイクロサテライトが増幅の成功、即ちウシDNAで適切なサイズのバンドを示し、ハムスターDNAでは増幅しなかった。これらは、3000−ラドパネルロスリン/ケンブリッジウシRHパネル(3000-rad panel Roslin/Cambridge bovine RH panel, Williams et al., 2002)で同定された。RH細胞株DNA 25ng、各プライマー0.5μM、200μM dNTP、3mM MgCl2、0.5UのTaqポリメラーゼ(BIOLine)を含有する全量20μlでPCR増幅を行った。反応条件は、94℃で3分間のタッチダウンから始めた後、93℃で30秒間、1サイクルごと0.5℃ずつ減少させる65〜45℃で30秒間のタッチダウン、及び72℃で1分間を40サイクル、そして72℃で5分間の最終伸長工程であった。PCR反応物を電気泳動し、PCR産物10μlを臭化エチジウム染色ミニゲル(2.5%アガロース)に充填し、2人の独立した観察者がアンプリコン(amplicons)の有無をスコア化した。複製物のパターンの間、又は異なる観察者のスコア間で幾らかの差異があった場合、PCR反応及びゲル電気泳動を繰り返した。幾つかのハイブリッドの結果が不明瞭なままであった場合、マーカーを廃棄した。
【0160】
ウシWGRHパネルでこれまでに同定されており、配置が既知(known assignments)のマーカーに対し、RHMAPPER(Soderlund et al., 1998)を用いて2点連鎖分析を行うことによって、ウシ染色体に対してマーカーを割り当てた(assigned)。それから、Williams et al(2002)によって記載されたように、Carthegeneソフトウェア(Schiex., 2002)を用いて、RHマップを構築した。ウシ染色体9上で、我々は、150個のマーカーによるラディエーションハイブリッドマップの情報を有する。
【0161】
CRIMAP2.4ソフトウェア(Green et al., 1990)を用いて、マーカーの序列及びマップ距離を推定した。我々は、MARCマップの序列に従ってマイクロサテライトマーカーを配置することによって我々の連鎖マップを構築することを開始した。次に、CRIMAPのBUILDオプションを作動させ、残りのマーカーを配置してから、新規のマイクロサテライトマーカー及びSNPマーカーを尤度が最も高い位置に挿入している。MARC(www.marc.usda.gov/genome/genome.html)、Ensembl(http://www.ensembl.org/Bos_taurus/index.html)、及びラディエーションハイブリッド(RH)情報を考慮し、マーカーの配置(emplacement)を再検討している。本発明によるBTA9のQTLマッピングに用いた最終の連鎖マップには、表51で列記されるように59個のマーカーが包含される。
【0162】
本発明によるBTA11のQTLマッピングに用いた最終の連鎖マップには、表52で列記されるような37個のマーカーが包含される。
【0163】
以下の表は、関連するQTLに用いたマーカーを示す。マーカーに対する任意のさらなる情報は、「http://www.marc.usda.gov/」、「http://www.ensembl.org/Bos_taurus/index. html」及び「http://www.ncbi.nih.gov/」に見出すことができる。
【0164】
【表51−1】

【0165】
【表51−2】

【0166】
【表52−1】


【0167】
【表52−2】

【0168】
表現型データ
雄ウシの娘を、乳房炎耐性及びSCSについてスコア化した。息子の形質に対する推定育種価(estimated breeding values, EBV)を、家系構造を無視する単一形質の最良線形不偏予測(Best Linear Unbiased Prediction, BLUP)動物モデルを使用して算出した。これらのEBVをQTL解析において使用した。個々の形質で使用される娘の登記項目(registrations)は以下の通り。
デンマークでの臨床型乳房炎:第一出産前5日〜出産後50日の期間での、臨床型乳房炎の治療症例。
スウェーデン及びフィンランドでの臨床型乳房炎:第一出産前7日〜出産後150日の期間での、臨床型乳房炎の治療症例。
デンマークでのSCS:第一出産後10日〜180日の期間での平均SCS。
スウェーデンでのSCS:第一出産後10日〜150日の期間での平均SCS。
フィンランドでのSCS:第一出産後10日〜305日の期間での平均SCS。
【実施例1】
【0169】
BTA9
統計分析
以下で示されるような多くの統計法が、乳房炎に関連する又はそれと連鎖する遺伝的マーカー、ひいては乳房炎耐性の決定に用いられた。
【0170】
QTL解析
QTLを同定するためには連鎖解析(LA)を用いるが、それは、偶然から予想されるより高頻度で系統内の疾患又は形質価に関連する染色体領域に、遺伝的マーカーを家系内でタイピングすることによって行う。かかる連鎖領域は、偶発的な遺伝的バリアントを包含する可能性が高い。データは一連のモデルを用いて解析した。最初に、すべての形質について多点回帰アプローチを用いた単一形質モデルを家系内で解析した。家系内で有意な効果を有する染色体は、家系間で見出されたQTLを検証し且つQTLを特徴付けるために、分散成分法(variance component method)を用いて解析した。
【0171】
回帰分析
マーカーの集団アレル頻度(allele frequencies)を、EMアルゴリズムを使用して推定した。その後、アレル頻度が誤差無く分かると仮定した。雄親におけるフェーズを出生児のマーカー型に基づいて決定した。その後、このフェーズが誤差無く分かると仮定した。各マップ位置での分離確率(segregation probabilities)は、Haldaneのマッピング関数を同時使用して、染色体上のすべてのマーカーからの情報を使用して算出した(Haldane, 1919)。表現型を分離確率に回帰した。有意性閾値(significance thresholds)を、順列検定を使用して算出した(Churchill and Doerge, 1994)。
【0172】
分散成分法(variance component method)
分散成分(VC)ベースの方法(Sorensen et al., 2003)を用いて、家系間連鎖解析を行った。VCを用いたLAでは、任意の2つの創始者ハプロタイプ(founder haplotypes)(Hs及びHm)のQTLアレル間の祖因同一性(Identity by descent, IBD)確率を0、即ち創始者ハプロタイプは関連しなかったと仮定する(Meuwissen et al., 2002)。雄親から父系性若しくは母系性のQTLアレルを受け継ぐ確率を計算するために(Freyer et al., 2004)、雄親ハプロタイプ及び息子の父系性遺伝ハプロタイプ(paternally inherited haplotypes)を用いた。そして、回帰アルゴリズム(Wang et al., 1995)を用いて、IBD行列を計算した。染色体に沿った各マーカーブラケット(
marker bracket)の中点(midpoint)でIBD行列を計算し、続く分散成分推定手順で用いた。QTLによって説明される相加的(additive)な遺伝的分散の総量の割合は、2σ2h/(2σ2h+σ2u)(式中、σ2h及びσ2uはそれぞれ、ハプロタイプ効果(haplotypes effect)及び相加的な多遺伝子性効果(polygenic effect)に関連する分散成分に対応する)と推定した。
【0173】
分散成分解析 単一形質単一QTL解析
各形質を、連鎖解析を使用して別個に解析した。完全なモデル(full model)は以下の通り表現され得る:
【0174】
【数1】

(式中、yはn個のEBVのベクトル(a vector of n EVBs)、Xは既知のデザイン行列、βは未知の固定効果のベクトル(ここでは単なる平均である)、Zは個体に関連する行列、uは相加的な多遺伝子性効果のベクトル、Wは個々の記録をそれぞれ未知の相加的なQTL効果と関連付ける既知の行列、qは個体の未知の相加的なQTL効果のベクトル、eは残余のベクトルである)。ランダム変数u、q及びeは、多変量であり(multivariate)、正規分布し且つ相互に独立すると仮定される(Lund et al., 2003)。
【0175】
多形質多QTL解析
多形質解析を行った。Lund et al., 2003に従って、モデル(1)をモデル(2)に記載されるように多形質多QTLモデルに拡張させることができる。
形質はnq個のQTLを用いた以下の線形混合モデルを使用してモデル化される:
【0176】
【数2】

(式中、yはt形質に対して記録されたn息子に対する観察値のベクトル(a vector of observations for n sons recorded on t traits)であり、μは全体の形質平均のベクトル(a vector of overall trait means)であり、Z及びWはそれぞれの息子の観察値を多遺伝子性効果及びQTL効果に関連付ける既知の行列であり、aは雄親及びその息子の多遺伝子性効果のベクトルであり、hiはi番目のQTLについての雄親及びその息子のQTLハプロタイプ効果のベクトルであり、eは残余のベクトルである)。ランダム変数a、hi及びeは多変量正規分布(multivariate normally distributed, MVN)であり且つ相互に無相関であると仮定される。特に、aは
【数3】

であり、hi
【数4】

であり、eは
【数5】

である。行列G、K及びEには、多遺伝子性効果、QTL効果及び残存効果のために、形質間の分散及び共分散(covariances)が含まれる。
記号
【数6】

は、クロネッカー積(Kronecker product)を表す。Aは多遺伝子性効果間の、共分散構造を示す相加的な関係行列(additive relationship matrix)であり、IBDiはi番目のQTLに対する効果間の、共分散構造を示す祖因同一性(IBD)行列であり、Iは単位行列(identity matrix)である。
【0177】
複合連鎖及び連鎖不平衡解析
複合連鎖及び連鎖不平衡解析(combined linkage and linkage disequilibrium analysis)において、任意の2つの創始者ハプロタイプのQTLアレル間のIBD確率を、Meuwissen and Goddard(2001)によって記載された方法を用いて計算した。この方法は、合体理論(coalescence theory)(Hudson, 1985)に基づき、隣接マーカーの状況同一性(identity-by-state, IBS)状態に依存しつつ、2つのハプロタイプが推定QTLにおいてIBDである確率を概算する。要するに、QTLでのIBD確率は、この位置を囲むアレルの周りのマーカーハプロタイプ(marker haplotypes)の類似性に基づいており、即ちこの位置近くに多くの(非)同一マーカーアレル(marker alleles)があれば、当該マップ位置における高い(低い)IBD確率が示唆される。IBD確率の実際の水準は有効集団サイズ(effective population size)、Neに影響される。両方のハプロタイプが共有するマーカーアレルを前提として、現世代と任意のベース世代(Tg世代前)との間の合体(coalescence)の確率を算出する(Hudson, 1985)。観察データからTg及びNeを推定するのは簡単ではない。シミュレーション研究は、QTL位置の推定がNe及びTgの選択に対して比較的鈍感であることを示す(Meuwissen and Goddard, 2000)。したがって、我々はTg=100及びNe=100の値を用いた。10個のマーカーのウィンドウを考慮し、IBD確率を計算した。我々は4マーカーウィンドウも用いて、4マーカーが10マーカーのハプロタイプによって既に同定されたピークを再現するのに十分であるかどうかを調べるために、LDLAピーク領域でのIBD確率を計算した。創始者ハプロタイプを機能的に明瞭に区別できるクラスターにグループ分けした。我々は、距離尺度(distance measure)として(1-IBDij)を用い、階層的クラスター形成アルゴリズム平均連鎖(hierarchical clustering algorithm average linkage)を適用し、すべての創始者ハプロタイプ間の遺伝的関係(genetic relationship)を表す有根樹状図(rooted dendrogram)を作製した。この樹状図を根元から下向きにスキャンし、すべての合体ハプロタイプ(coalescing haplotypes)が距離尺度(1-IBDij)<Tcを有するようなノード(nodes)に達するまで分枝を切断する。クラスターは、共通のノードに合体(coalesce)するハプロタイプの群として定義される。クラスター内のハプロタイプは同一のQTLアレルを有すると仮定し(IBD確率=1.0)、異なるクラスター由来のハプロタイプは異なるQTLアレルを有し、したがって独立していると考える(IBD確率=0)。そのため、連鎖不平衡情報に対応するIBD行列の上部は明瞭に区別できる創始者ハプロタイプに対応する単位行列である。息子の父系性ハプロタイプ(paternal haplotypes)における連鎖情報に対応するIBD行列の下部は、回帰アルゴリズムを用いて構築する(Wang et al., 1995)。IBD行列は、それぞれのマーカー区間(marker interval)の中点で計算し、続く分散成分推定手順で用いた。
【0178】
パラメーターの推定
平均情報被制限最大尤度アルゴリズム(average information restricted maximum likelihood algorithm)を用いて、分散成分を推定した(Jensen et al., 1997)。このモデルにおけるランダム効果に関する分散成分に関して、被制限尤度を最大にした。特定位置の格子(a grid of specific positions)に渡る一連の被制限尤度を最大にすることによって、QTL位置の被制限尤度のプロファイルが得られる(Sorensen et al., 2003)。染色体に沿って、それぞれのマーカーブラケットの中点でパラメーターを推定した。
【0179】
有意水準
ゲノム毎の第一種の過誤(genome-wise type I error)を制御するおおよその閾値水準を計算するための迅速な方法を使用して、分散成分解析についての有意性閾値を算出した(Piepho, 2001)。QTLの存在についての仮説検定は、尤度比検定(likelihood ratio test, LRT)統計、LRT=−2ln(Lreduced−Lfull)の漸近分布(asymptotic distribution)を基礎とした(式中、Lreduced及びLfullはそれぞれ、低減モデル(reduced model)及び完全モデル(full model)下での最大尤度であった)。低減モデルでは、解析する染色体についてのQTL効果を常に排除した。この方法は、順列手順(permutation procedures)の代替法であり、複雑な状況で適用可能である。この方法には、染色体に沿ったそれぞれの推定QTL位置からのLRT、染色体の数、LRTについての自由度(degree of freedom, df)(df=Hfullのパラメーター数−Hreducedのパラメーター数)、及び染色体毎の第一種の過誤発生率(chromosome-wise type I error rate)が必要である。染色体毎の5%の有意水準が有意であると見なした。
【0180】
BTA9の結果
表53に、BTA9についての回帰分析の結果を示す。図1〜図8は、回帰分析についてのQTLを示すグラフである。家系間QTL解析でQTLを検出するために、分散成分法を使用した。図9〜図16は、分散成分ベースの方法におけるQTLについてのLD、LDLA及びLDプロファイルを示す。図17〜図20はハプロタイプ効果を示す。
【0181】
デンマークレッド
家系内回帰分析によって、CM及びSCSについてのQTLがDR品種中の2つの家系で分離されることが分かった。2つの形質についてのQTLは同じ区間には位置していなかった。VC法を用いた家系間連鎖解析では、QTL効果は有意ではなかった。LDLA及びLD解析では、マーカーBMS2819及びINRA144間の74.08cMで高いQTLピークが観察された。LD解析におけるLRTピーク(13.6)は、複合LDLA解析で観察されたLRTピーク(8.51)よりも高かった。このQTLは、相加的な遺伝的分散の44%及びCMについての表現型分散の22%を説明した。初期設定としては(by default)、10個のマーカーのハプロタイプ(推定位置の各側で5つのマーカー)を用いて、或る位置のIBD確率を求めた。我々は、4つのマーカーのハプロタイプ(即ち推定位置の各側で2つのマーカー)も用いて、これらの4マーカー(BM4208−BMS2819−INRA144−INRA084)ハプロタイプ内に同様のLDLA/LDピークを観察した。この品種におけるこれらの4マーカーブラケット内では、SCSについてのLDLA複合ピークも観察された。IBD確率が0.90以上の雌親(dam)ハプロタイプを一緒にクラスター化した。DRにおいては、クラスター形成前では305個の創始者ハプロタイプが存在し、クラスター形成後クラスターは54個に減少した。5つのクラスターの頻度が5%を超え、最大クラスターの頻度は10%であり、また5つの雄親ハプロタイプが最大クラスターとクラスター形成した。これらの54個のハプロタイプ及びまた雄親から受け継いだハプロタイプについてのハプロタイプ効果が推定された。高及び低乳房炎耐性に関連したハプロタイプが同定された。
【0182】
フィンランドエアシャー
家系内回帰分析をフィンランドエアシャー家系で行った場合に、CM及びSCSに影響を及ぼすQTLが分離することが見出された。CMに対するQTLは58〜79cMの間に位置していた。SCSに影響を及ぼすQTLは、32〜44cM(LRT統計ピークは37cMにある)間に位置していた。LAプロファイル上、CMのQTLについての複合LDLAピークがマーカーBM4208−BMS2819−INRA144内に観察された。CMについて、同じ領域でLDピークも観察された。LAプロファイル上、SCsについての1つのLDLAピークがマーカーDIK2810及びDIK5364間の38cMで観察された。CMにおける全分散(total variance)の4%が74cMのQTLによって説明され、このQTLはFAにおいてSCSには効果を示さなかった。38cMでのQTLは、SCSにおいては全分散の18%を説明し、これはCMに対してほんのわずかしか効果がなかった。クラスター形成確率を0.90にすると、CMにおける最大LDLAピーク(即ちマーカーBM4208及びBMS2819間の区間中心(mid interval))では、442個の創始者ハプロタイプが38個のクラスターにグループ分けされた。頻度が5%を超えるクラスターが9つ存在していた。最大クラスターの頻度は14%であった。高及び低乳房炎耐性に関連したハプロタイプが同定された。
【0183】
スウェーデンレッドアンドホワイト
DR及びFAウシと同様に、スウェーデンレッドアンドホワイトウシでも、家系内回帰分析を行った場合に、CM及びSCSに影響を及ぼすQTLがBTA9上で分離することが見出された。CM及びSCSについてのQTLは共に同じ区間に位置していた。CMのQTLは家系間LA解析において有意であった(P<0.01)。SCSのQTLは家系間LA解析において有意ではなかった。SCSについてのLRTピークは73cMにあった。家系間解析におけるCMのQTLについての検定統計のピークは67.4cMにあり、QTL区間は59〜81cMであった。LRT統計は家系間LA解析では極めて有意であったが、SRBウシにおけるこのQTLについてのLAプロファイル上ではLDLAピークは全く観察されなかった。LA解析におけるLRT統計のピーク位置で、QTL分散はCM形質の全分散の25%であった。DR及びFA品種のCMのQTLについてのLDLAピークは、SRBで観察されたLAプロファイルに含まれる。CMについてのSRBデータにおいて、LDLAピークは全く観察されなかったが、DR及びFAにおいてピークが位置していたマーカー区間でSRBにおいて多くのクラスター形成が起こっていた。例えば、FAにおいて最大のLDLAピークが位置している、BM4208−BMS2819の区間中心(mid interval)で、及びDRにおいて最大のLDLAピークが位置している隣接区間(BMS2819−INRA144)で、SRBにおける全400個の創始者ハプロタイプの中で、それぞれ37個及び48個のクラスターが存在していた。
【0184】
デンマークホルスタイン
デンマークホルスタインウシにおいても、CM及びSCSに影響を及ぼすQTLの分離が家系内回帰分析で示された。CMのQTLは、VC法を用いた家系間LA解析で有意であった(P<0.01)。LAプロファイル上で小さなLDLAピークが観察されたが、DHにおけるQTLについては確証的なLDピークは見られなかった。LAにおけるCMのQTLについてのLRT統計は42.9cMで最大になり、そのLRT統計は10.6であった。QTL区間は29〜51cMと極めて大きく広がっていた。CMのQTLについての1つの小さなLDピークが、DR及びFA集団において74cMで観察されたLDピークと一致している。SCCのQTLの検定統計のピークは48.7cMにあり、区間は44〜58cMである。全分散のうち、各LAにおいて最大ピークをとるQTLによって説明される部分は、それぞれCMについて27%、及びSCSについて17%であった。CMについての最大LDピークは73.35cMにあり、この領域でDRについて高いLDピークが観察された。DHにおいては、SCSについてのLDピークは全く観察されなかった。
【0185】
品種間解析
品種内(within-breed)LA、LDLA及びLD解析によって、CMに影響を及ぼすQTLが、2つ以上の集団において74cM近辺で分離されることが分かった。したがって、この調査における様々な品種間のデータを組合せて、品種間QTL解析を行った。品種間のQTL解析の結果を表20、21及び22に示す。品種内解析を行った場合、DR及びFAウシにおけるCMのQTLについてのLDLAピークが隣接マーカー区間に位置していた。しかしながら、DR及びFAの複合データを解析した場合は、マーカーブラケット(BMS2819−INRA144)においてCMの高いLDLAピークが観察されたが、これはLDピークとも一致している。FA及びSRBデータの複合解析(combined analysis)によっては、LAプロファイル上でより高いLDLAピークは与えられなかったが、74cMの同じマーカー区間でLDピークが観察された。DR、FA及びSRBの複合データ解析によっても、LA上の同じ領域(即ちBM4208−BMS2819−INRA144)においてより高いLDLAピークが与えられた。LDピークも同じ位置に存在し、LA上のより高いLDLAピークを証明した。DR及びFAの連合解析(joint analysis)によって、SCSのQTLについて、マーカーDIK2810及びDIK5364間の38cMで高いLDLAピークが示された。また、FA及びSRBの複合解析においては、SCSのQTLについて同じ位置でLDLAピークが観察された。しかしながら、DR、FA及びSRBを一緒に解析すると、このLDLAピークは消失した。
【0186】
多形質解析
SCSは乳房炎耐性の指標形質(indicator trait)である。CMに関与する多くの遺伝子がSCSに対する効果も有していることが予測された。したがって、BTA9上で分離するQTLが両方の形質に対する多面発現効果(pleiotropic effect)を有しているかどうか、又はこのQTLが連鎖QTLであるかを調べるために、CM及びSCSの多形質解析を行った。DRデータの単一形質LDLA解析によっては、同じマーカー区間(即ちBMS2819及びINRA144間)でCM及びSCSのLDLAピークが示されたが、二形質の複合解析によっては、マーカーSLU2及びC6orf93間の69.1cMでLDLAピークが与えられた。品種内解析において、CM及びSCSの両方に影響を及ぼすQTLを有するモデルについて、FA、SRはLDLAピークを示さなかった。しかしながら、3つの品種DR、FA及びSRBを組合せ、二形質モデルで解析した場合、LRTの統計が19.7であるLDLAピークが、INRA144及びINRA084のマーカー区間で観察された。
【0187】
ハプロタイプ解析
上記のQTL細密(fine)マッピングの結果は、4マーカー領域(BM4208−BMS2819−INRA144−INRA084)内に、CMについて分離するQTLを指し示す。したがって、創始者ハプロタイプのクラスター形成及びハプロタイプ効果を、マーカーBMS2819及びINRA144間の中点で調べた。これは品種内で行われ、また3つの品種DR、FA及びSRB間でも行われた。これは、これらの3つの品種がその起源で関係しているためである。高及び低乳房炎耐性に関連したハプロタイプが同定された。
【0188】
【表53−1】

【0189】
【表53−2】

【0190】
【表54】

【0191】
【表55】

【0192】
【表56】

【実施例2】
【0193】
BTA11
統計分析
下記のような多くの統計法が、乳房炎、ひいては乳房炎耐性に関連した、又はそれと連鎖した遺伝的マーカーの決定に用いられた。
【0194】
QTL解析
偶然予測されるものよりも高頻度で、系統内の疾患又は形質価(trait values)に関連する染色体領域に対して、家系における遺伝的マーカーをタイピングすることによってQTLを同定するのに、連鎖解析(LA)を用いる。かかる連鎖領域は偶発的な遺伝的バリアントを含有する可能性が高い。データは一連のモデルで解析した。3つの相補的アプローチを用いた:(i)GDQTLソフトウェア(B. Guldbrandsten, 2005 personal communication)を用いる回帰ベースの方法(Haley and Knott, 1992)による半同胞家系内分離解析(within half-sib family segregation analysis)、(ii)分散成分法を用いる家系間連鎖解析、及び(iii)分散成分法を用いる複合連鎖不平衡連鎖解析(combined linkage disequilibrium linkage analysis, combined LDLA)。それぞれの形質についての雄親のQTL分離状態を決定するために、GDQTLを用いることによって、それぞれの家系を個々に解析した。それぞれの雄親についての染色体毎の有意水準を求めるのに順列検定を用いた(n=10000)(Churchill and Doerge, 1994)。次の工程は、形質及びQTL位置にかかわらず、QTLを分離する家系からのデータセットを組合せた、分散成分ベースの方法(Sorensen et al., 2003)を用いる家系間連鎖解析であった。Piepho(2001)によって提示された方法を用いて閾値を算出した。第3の工程は、すべての分離家系及び非分離家系を包含する複合LDLA解析(Lund et al., 2003)であった。連鎖QTLから多面発現性(pleiotropic)QTLを分けるために、多形質モデル及び多QTLモデルを解析した。QTLが、2つ以上の品種においてBTA11の同じ領域に分離することが観察された場合、品種間のデータを組合せてLDLA解析を行った。
【0195】
分散成分法
分散成分(VC)ベースの方法(Sorensen et al., 2003)を用いて、家系間連鎖解析を行った。VCを用いたLAにおいて、任意の2つの創始者ハプロタイプ(Hs及びHm)のQTLアレル間の祖因同一性(IBD)確率を0(即ち創始者ハプロタイプは関連がなかった)と仮定する(Meuwissen et al., 2002)。雄親から父系性又は母系性のQTLアレルを受け継ぐ確率を計算するために、雄親ハプロタイプ及び息子の父系性遺伝ハプロタイプを用いる(Freyer et al., 2004)、そして、再帰アルゴリズムを用いてIBD行列を計算した(Wang et al., 1995)。IBD行列は、染色体に沿った各マーカーブラケットの中点で計算し、続く分散成分推定手順で用いた。相加的な遺伝的分散の全体のうち、QTLによって説明される割合は、2σ2h/(2σ2h+σ2u)と推定した(式中、σ2h及びσ2uはそれぞれ、ハプロタイプ効果及び相加的な多遺伝子性効果に関連する分散成分に対応する)。
【0196】
分散成分解析 単一形質単一QTL解析
各形質は、連鎖解析を使用して別個に解析した。完全なモデルは以下の通り表現され得る:
【0197】
【数7】

(式中、yはn個のEBVのベクトル(a vector of n EVBs)、Xは既知のデザイン行列、βは未知の固定効果のベクトル(ここでは単なる平均である)、Zは個体に関連する行列、uは相加的な多遺伝子性効果のベクトル、Wは各個体の記録を未知の相加的なQTL効果と関連付ける既知の行列、qは個体の未知の相加的なQTL効果のベクトル、eは残余のベクトルである)。ランダム変数u、q及びeは、多変量であり、正規分布し且つ相互に独立すると仮定される(Lund et al., 2003)。
【0198】
多形質多QTL解析
多形質解析を行った。Lund et al., 2003に従って、モデル(1)をモデル(2)に記載されるように多形質多QTLモデルに拡張させることができる。
この形質はnq個のQTLで以下の線形混合モデルを用いてモデル化される:
【0199】
【数8】

(式中、yはt形質に対して記録されたn息子に対する観察値のベクトル(a vector of observations for n sons recorded on t traits)であり、μは全体の形質平均のベクトル(a vector of overall trait means)であり、Z及びWはそれぞれの息子の観察値を多遺伝子性効果及びQTL効果に対して関連付ける既知の行列であり、aは雄親及びそれらの息子の多遺伝子性効果のベクトルであり、hiはi番目のQTLについての雄親及びそれらの息子のQTLハプロタイプ効果のベクトルであり、eは残余のベクトルである)。ランダム変数a、hi及びeは多変量正規分布(multivariate normally distributed, MVN)であり且つ相互に無相関であると仮定される。特に、aは
【数9】

であり、hi
【数10】

であり、eは
【数11】

である。行列G、K及びEには、多遺伝子性効果、QTL効果及び残存効果のために、形質間で分散及び共分散が含まれる。
記号
【数12】

は、クロネッカー積(Kronecker product)を表す。Aは多遺伝子性効果間の、共分散構造を説明する相加的な関係行列(additive relationship matrix)であり、IBDiはi番目のQTLに対する効果間の、共分散構造を説明する祖因同一性(IBD)行列であり、Iは単位行列である。
【0200】
回帰分析
マーカーでの集団アレル頻度を、EMアルゴリズムを使用して推定した。その後、アレル頻度は誤差無く分かったと仮定した。雄親におけるフェーズは出生児のマーカー型に基づいて求めた。その後、このフェーズは誤差無く分かったと仮定した。各マップ位置での分離確率は、Haldaneのマッピング関数を同時使用して、染色体上のすべてのマーカーからの情報を使用して算出した(Haldane, 1919)。表現型を分離確率に回帰した。有意性閾値を、順列検定を使用して算出した(Churchil and Doerge, 1994)。
【0201】
パラメーターの推定
平均情報被制限最大尤度アルゴリズム(average information restricted maximum likelihood algorithm)を用いて、分散成分を推定した(Jensen et al., 1997)。このモデルにおけるランダム効果に関する分散成分に関して、被制限尤度を最大にした。特定位置の格子(a grid of specific positions)にわたる一連の被制限尤度を最大にすることによって、QTL位置の被制限尤度のプロファイルが得られる(Sorensen et al., 2003)。染色体に沿って、それぞれのマーカーブラケットの中点でパラメーターを推定した。相加的な遺伝的分散の全体のうち、QTLによって説明される割合は、2σ2h/(2σ2h+σ2u)と推定した(式中、σ2h及びσ2uはそれぞれ、ハプロタイプ効果及び相加的な多遺伝子性効果に関連する分散成分に対応する)。
【0202】
IBD確率の推定
連鎖解析:任意の2つの創始者ハプロタイプ(Hs及びHm)のQTLアレル間のIBD確率を0(即ち創始者ハプロタイプは関連しなかった)と仮定する(Meuwissen et al., 2002)。雄親から父系性又は母系性のQTLアレルを受け継ぐ確率を計算するために、雄親ハプロタイプ及び息子の父系性遺伝ハプロタイプを用いる、そして、回帰アルゴリズムを用いてIBD行列を計算した(Wang et al., 1995)。IBD行列は、染色体に沿った2cM区間毎で計算し、続く分散成分推定手順に用いた。
【0203】
複合連鎖及び連鎖不平衡解析
複合連鎖及び連鎖不平衡解析(combined linkage and linkage disequilibrium analysis)において、任意の2つの創始者ハプロタイプのQTLアレル間のIBD確率を、Meuwissen and Goddard(2001)によって記載された方法を用いて計算した。この方法は、合体理論(coalescence theory)(Hudson, 1985)に基づき、隣接マーカーの状況同一性(identity-by-state, IBS)状態に依存しつつ、2つのハプロタイプが推定QTLにおいてIBDである確率を概算する。要するに、QTLでのIBD確率は、この位置を囲むアレルの周りのマーカーハプロタイプの類似性に基づいており、即ちこの位置近くに多くの(非)同一マーカーアレルがあれば、当該マップ位置における高い(低い)IBD確率が示唆される。IBD確率の実際の水準は有効集団サイズ、Neに影響される。両方のハプロタイプが共有するマーカーアレルを前提として、現世代と任意のベース世代(Tg世代前)との間の合体(coalescence)の確率を算出する(Hudson, 1985)。観察データからTg及びNeを推定するのは簡単ではない。シミュレーション研究は、QTL位置の推定がNe及びTgの選択に対して比較的鈍感であることを示す(Meuwissen and Goddard, 2000)。したがって、我々はTg=100及びNe=100の値を用いた。10個のマーカーのウィンドウを考慮し、IBD確率を計算した。我々は異なるマーカー−ウィンドウ、例えば6マーカー、4マーカー等も用いて、少ないマーカーであっても、10マーカーのハプロタイプで検出したQTL分散を説明するのに十分であるかどうかを調べるために、LDLAピーク領域でのIBD確率を計算した。創始者ハプロタイプを明瞭に区別できるクラスターにグループ分けした。我々は、距離尺度として(1-IBDij)を用い、階層的クラスター形成アルゴリズム平均連鎖を適用し、すべての創始者ハプロタイプ間の遺伝的関係を表す有根樹状図を作製した。この樹状図を根元から下向きにスキャンし、すべての合体ハプロタイプ(coalescing haplotypes)がTcより小さい距離尺度(1-IBDij)を有するようなノードに達するまで分枝を切断する。クラスターは、共通のノードに合体(coalesce)するハプロタイプの群として定義される。クラスター内のハプロタイプは同一のQTLアレルを有すると仮定し(IBD確率=1.0)、異なるクラスター由来のハプロタイプは異なるQTLアレルを有し、したがって独立していると考える(IBD確率=0)。そのため、連鎖不平衡情報に対応するIBD行列の上部は明瞭に区別できる創始者ハプロタイプに対応する単位行列である。息子の父系性ハプロタイプにおける連鎖情報に対応するIBD行列の下部は、回帰アルゴリズムを用いて構築する(Wang et al., 1995)。IBD行列を、それぞれのマーカー区間の中点で計算し、続く分散成分推定手順で用いた。
【0204】
有意水準
ゲノム毎の第一種の過誤(genome-wise type I error)を制御するおおよその閾値水準を計算するための迅速な方法を使用して、分散成分解析についての有意性閾値を算出した(Piepho, 2001)。QTLの存在についての仮説検定は、尤度比検定(likelihood ratio test, LRT)統計、LRT=−2ln(Lreduced−Lfull)の漸近分布(asymptotic distribution)を基礎とした(式中、Lreduced及びLfullはそれぞれ、低減モデル(reduced model)及び完全モデル(full model)下での最大尤度であった)。低減モデルでは、解析する染色体についてのQTL効果を常に排除した。この方法は、順列手順(permutation procedures)の代替法であり、複雑な状況で適用可能である。この方法には、染色体に沿ったそれぞれの推定QTL位置からのLRT、染色体の数、LRTについての自由度(degree of freedom, df)(df=Hfullのパラメーター数−Hreducedのパラメーター数)、及び染色体毎の第一種の過誤発生率(chromosome-wise type I error rate)が必要である。染色体毎の5%の有意水準が有意であると見なした。
【0205】
BTA11の結果
表57はBTA11上の回帰分析の結果を示す。分散成分法を用いたLA解析の結果は表58に示す。LDLA結果は表59に示し、LD解析結果は表60に示す。図21(FA)、図23(SRB)及び図25(レッド混合)は臨床型乳房炎形質についての回帰分析のためのQTLを示すグラフであり、図22(FA)、図24(SRB)及び図26(レッド混合)は体細胞スコア形質についてのQTLを示すグラフである。家系間QTL解析でQTLを検出するのに、分散成分法を用いた。臨床型乳房炎について、LA、LDLA及びLDを用いた分散成分ベースの方法で得られたQTLプロファイルは、図27(FA)、図28(IBDのための4マーカーウィンドウを用いたFA)、及び図31(レッド混合)に示す。SCS形質についてのLA、LDLA及びLD解析におけるQTLプロファイルは、図29(FA)、図30(SRB)及び図32(レッド混合)に示す。4マーカーのハプロタイプを用いたフィンランドエアシャーにおける最大LDLAピークでの臨床型乳房炎形質に対する大きいクラスターの効果は図33に示す。
【0206】
フィンランドエアシャー
回帰分析を用いたBTA11についての8つのフィンランドエアシャー(FA)半同胞家系全データ間解析によって、5%の水準で有意な2つのQTLが得られた。臨床型乳房炎に影響を及ぼすQTLは11.3cMに位置し、体細胞スコアに影響を及ぼすQTLは64.1cMにあった。1つの家系が乳房炎のQTLで有意であった一方で、他の2つの家系は有意性閾値に達していた。これら3つの家系のQTL区間は重複していたが、広い領域に広がっていた。2つのフィンランドエアシャー家系はSCSのQTLについて有意であり、これら2つの家系におけるQTL位置は6cM離れていた。分散成分法を用いた臨床型乳房炎についての家系間連鎖解析の尤度比検定(LRT)の最大統計は、14.2cMで5.74であった。複合連鎖不平衡及び連鎖解析(LDLA)を行った場合、16.8cMではっきりしたQTLピークが存在し、マーカーMNB−40及びAUP1間のLRTは11.82であった。臨床型乳房炎についてのLD解析の最大LRTは、マーカーDIK4637及びUMBTL103間の20.6cMで観察されたが、MNB−40及びAUP1間のLD(LRT=3.18)についても幾らかの証拠が存在していた。臨床型乳房炎の分散のうち、最大LDLAピークにおけるQTLによって説明される部分は、全分散の15%であった。初期設定としては(by default)、IBD確率を推定するのに10マーカーウィンドウを用いた。4マーカーウィンドウでLDLA解析を繰り返した(図28)。マーカーAUP1及びBM716間の17.8cMで、はっきりしたQTLピーク(LRT=9.7)が観察された。この4マーカーの間の区間は2.1cMである。2つのフィンランドエアシャー雄親家系は、SCSのQTLについて55〜70cMの領域で分離していた。最も可能性のあるQTL位置は、すべてのFA家系を包含する多点回帰分析(multipoint regression analysis)において64.1cMであった。分散成分を用いたLA解析は、62.8cMで最大LRT(6.6)を有する。QTL区間は相当な幅を維持したままであった。この領域でLDが観察されなかったので、LDLA解析によってQTL区間を絞り込むことはできなかった。多面発現性QTLモデル(即ちBTA11上の臨床型乳房炎及びSCSの両方の形質に影響を及ぼすQTL)は収束しなかった。二連鎖QTLモデル(two linked QTL model)では、臨床型乳房炎のQTLは14.2cM、SCSのQTLは61.6cMに位置し、LRTは16.62であった。したがって、それぞれ1つが1つの形質に影響を及ぼす2つのQTLがBTA11上で分離されることを結論付けることができる。
【0207】
スウェーデンレッドアンドホワイト
スウェーデンレッドアンドホワイト(SRB)の家系間の多点回帰分析によって、SCSに影響を及ぼすQTLはBTA11上で分離し、QTLの最も可能性のある位置は61.2cMであることが示された。2つの家系は当該QTLについて有意であった。これらの2つの家系におけるQTLの考えられる位置は20cM離れていた(59.8及び40.4cM)。SRBにおいて家系間連鎖分析を分散成分法を用いて行った場合、SCSについてのQTL区間は非常に大きかった。QTL区間内のLDが不十分であったため、LDLA解析によってQTL区間を絞り込むことはできなかった。SCSに影響を及ぼす二連鎖QTLが30〜70cMの領域に位置して存在するかどうかを調べるのに、2-QTLモデルを使用した。SCSに影響を及ぼす61.2cMのQTLを固定し、SCSに影響を及ぼす別のQTLを求めてこの領域をスキャンした。しかしながら、この領域でSCSに影響を及ぼす第二のQTLについての痕跡は全くなかった。このQTLはSRBウシにおける臨床型乳房炎に対する多面発現効果を有していない。
【0208】
品種間解析
スウェーデンレッドアンドホワイト種はノルディックエアシャーウシ種に密接に関連している(Holmberg and Andersson-Eklund, 2004)。SCSに影響を及ぼすBTA11上のQTLは、FA(62.8cM)及びSRB(61.4cM)の両方において分離することが観察された。したがって、これら2つの品種のデータは、BTA11上のQTL細密マッピングのために組み合わされた。1つのデンマークレッド(DR)家系が、BTA11上のQTLについて臨床型乳房炎では56.0cMで、及びSCSでは66.9cMで分離して観察された。デンマークレッドウシはFA及びSRBにも歴史的に関連している。したがって、DR家系は、BTA11の連合解析(joint analysis)のために13のFA及びSRB家系に包含させた。分散成分を用いた臨床型乳房炎形質についての家系間連鎖解析は14.2cMで最大LRT(4.72)を有していた。レッドのデータの連合解析(joint Red data analysis)でFA内解析よりもLAピークがより小さい理由は、(BTA11の近位末端においてQTLについて分離しない)SRB及びDR家系を包含しているためであった。レッドの複合データについてのLDLAピークは16.8cMにあった(LRT=10.1)。レッドのデータにおいてLDの最大の証拠は18.2cMにあったが、マーカーMNB−40及びAUP1間の最大LDLAピーク(LRT=3.8)にLDの証拠が存在した。複合レッドデータ解析(combined Red data analysis)において、SCSに影響を及ぼすQTLは大きい区間を有していた(20cM)。分散成分を用いた連鎖解析のLRTは、62.4で14.6であった。最大のLDLAピークは、マーカーMS2177及びHELMTT44間の61.4cMにあった。LA解析よりもLDLA解析においてLRTが小さい理由は、SCSのQTL区間内でLDが不足しているためであった。家系内及び家系間の連鎖解析(回帰及び分散成分の両方)から、BTA11上でSCSのQTLが分離するという有力な証拠が存在したが、QTL区間内でLDが不足していたために、QTL位置を絞り込むことはできなかった。
【0209】
ハプロタイプ効果の推定
4マーカーウィンドウを用いたLDLA解析によって、マーカーAUP1及びBM716間の17.8cMに臨床型乳房炎のQTLが位置することが示された。FA、SRB及び1つのDR家系の連合解析(joint analysis)によって、この位置でのQTL情報は主にフィンランドエアシャー家系に由来していることが示された。したがって、マーカーAUP1及びBM716間の中点におけるクラスター及びその効果は、フィンランドエアシャーのみで調べた。この位置で、340個の創始者ハプロタイプが63個のクラスターに合体(coalesce)した。頻度が5%を超える8つのクラスターが存在し、最大のクラスターの頻度は9.7%であった。2つの祖父ハプロタイプを包含する32個のハプロタイプを有する1つのクラスターは、表現型の標準偏差が−0.13の推定効果を有していた。
【0210】
臨床型乳房炎及びSCSの形質についてのBTA11上のQTL細密マッピングによって、臨床型乳房炎に影響を及ぼす1つのQTLがフィンランドエアシャーウシにおいて分離し、SCSに影響を及ぼす1つのQTLがフィンランドエアシャー及びスウェーデンレッドアンドホワイトウシの両方で分離していることが確認された。LDLA解析によって、2.1cMの区間へ臨床型乳房炎が細密にマッピングされた。QTL区間内で連鎖不平衡が不足していたために、BTA11上のSCSに影響を与えるQTLは細密にマッピングすることができなかった。
ハプロタイプ解析
【0211】
上記のQTL細密マッピングの結果は、4マーカー領域(MNB−40−AUP1−BM716−DIK2653)内にCMについて分離するQTLを指し示す。したがって、創始者ハプロタイプのクラスター形成及びハプロタイプ効果を、マーカーAUP1及びBM716間の中点で検査した。また、上記のQTL細密マッピング結果は、4マーカー領域(BM304−INRA177−UMBTL20−RM96 INRA177)内のSCSについて分離するQTLを指し示す。したがって、創始者ハプロタイプのクラスター形成及びハプロタイプ効果を、それぞれ、マーカーINRA177及びUMBTL20間の中点で検査した。
【0212】
これは、FA内で行われた。高及び低乳房炎耐性に関連するハプロタイプが同定された。表61及び表62を参照されたい。
【0213】
【表57】


【0214】
【表58】

【0215】
【表59】

【0216】
【表60】

【0217】
【表61】


【表62】

【図面の簡単な説明】
【0218】
【図1】デンマークレッド家系の特徴である乳房炎耐性に関連したBTA9のゲノムスキャン(Genome scan)を示す図である。数字は、「血統書(herdbook)番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0219】
【図2】デンマークレッド家系の特徴である体細胞数に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0220】
【図3】フィンランドエアシャー家系の特徴である乳房炎耐性に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0221】
【図4】フィンランドエアシャー家系の特徴である体細胞数に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0222】
【図5a】スウェーデンレッドアンドホワイト家系の特徴である乳房炎耐性に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0223】
【図5b】スウェーデンレッドアンドホワイト家系の特徴である乳房炎耐性に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0224】
【図6a】スウェーデンレッドアンドホワイト家系の特徴である体細胞数に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0225】
【図6b】スウェーデンレッドアンドホワイト家系の特徴である体細胞数に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0226】
【図7】デンマークホルスタイン家系の特徴である乳房炎耐性に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0227】
【図8】デンマークホルスタイン家系の特徴である体細胞数に関連したBTA9のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0228】
【図9】デンマークレッドウシにおける74.08cMでのマーカーBMS2819及びINRA144間でLDLA/LDピークを示す乳房炎耐性形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0229】
【図10】デンマークレッドウシにおける74.075cMでのマーカーBMS2819及びINRA144間でLDLAピークを示す体細胞数形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0230】
【図11】フィンランドエアシャー種におけるマーカーBM4208及びINRA144間でLDLAピークを示す乳房炎耐性形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0231】
【図12】フィンランドエアシャー及びデンマークレッドウシの組合せにおけるマーカーBMS2819及びINRA144間でLDLA/LDピークを示す乳房炎耐性形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0232】
【図13】スウェーデンレッドアンドホワイトウシにおける60〜80cMでのマーカーBMS2819及びINRA144間でLAピークを示す乳房炎耐性形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0233】
【図14】マーカーBMS2819及びINRA084間でLDLAの証拠があるスウェーデンレッドアンドホワイトウシにおけるSCSについてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0234】
【図15】フィンランドエアシャー、デンマークレッドウシ及びスウェーデンレッドアンドホワイトウシの複合解析におけるマーカーBM4208及びINRA144間でLDLA/LDピークを示す乳房炎耐性形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0235】
【図16】デンマークホルスタインウシ乳房炎耐性形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0236】
【図17】デンマークレッドにおける74.08cMでの乳房炎耐性に対するハプロタイプ効果を示す図である。ハプロタイプ効果がy軸であり、ハプロタイプ数がx軸である。初めは、雌親ハプロタイプの大きいクラスターが存在し、続いて雄親ハプロタイプが存在する。
【0237】
【図18】フィンランドエアシャーにおける74.08cMでの乳房炎耐性に対するハプロタイプ効果を示す図である。ハプロタイプ効果がy軸であり、ハプロタイプ数がx軸である。初めは、雌親ハプロタイプの大きいクラスターが存在し、続いて雄親ハプロタイプが存在する。
【0238】
【図19】デンマークレッドアンドフィンランドエアシャーの組合せにおける74.08cMでの乳房炎耐性に対するハプロタイプ効果を示す図である。ハプロタイプ効果がy軸であり、ハプロタイプ数がx軸である。初めは、雌親ハプロタイプの大きいクラスターが存在し、続いて雄親ハプロタイプが存在する。
【0239】
【図20】デンマークレッド、フィンランドエアシャー及びスウェーデンレッドアンドホワイトの組合せにおける74.08cMでの乳房炎耐性に対するハプロタイプ効果を示す図である。ハプロタイプ効果がy軸であり、ハプロタイプ数がx軸である。初めは、雌親ハプロタイプの大きいクラスターが存在し、続いて雄親ハプロタイプが存在する。
【0240】
【図21】フィンランドエアシャー家系の特徴である臨床型乳房炎形質に関するBTA11のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0241】
【図22】フィンランドエアシャー家系の特徴である体細胞スコア形質に関するBTA11のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0242】
【図23】スウェーデンレッドアンドホワイト家系の特徴である臨床型乳房炎形質に関するBTA11のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0243】
【図24】スウェーデンレッドアンドホワイト家系の特徴である体細胞スコア形質に関するBTA11のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0244】
【図25】1つのデンマークレッド、8つのフィンランドエアシャー、及び4つのスウェーデンレッドアンドホワイトの家系の特徴である臨床型乳房炎形質に関するBTA11のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0245】
【図26】1つのデンマークレッド、8つのフィンランドエアシャー、及び4つのスウェーデンレッドアンドホワイトの家系の特徴である体細胞スコア形質に関するBTA11のゲノムスキャンを示す図である。数字は、「血統書番号」を指す。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、F値で表現されるQTL解析の検定統計を表す。
【0246】
【図27】フィンランドエアシャーウシにおけるLDLA/LDピークを示す臨床型乳房炎形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0247】
【図28】フィンランドエアシャーウシにおける、4マーカーのハプロタイプを有するLDLA/LDピークを示す臨床型乳房炎形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0248】
【図29】フィンランドエアシャーウシにおけるLA、LDLA及びLDプロファイルを示す体細胞スコア形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0249】
【図30】スウェーデンレッドアンドホワイトウシにおけるLA、LDLA及びLDプロファイルを示す体細胞スコア形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0250】
【図31】ノルディックレッド種の14家系の複合解析においてLA、LDLA及びLDプロファイルを示す臨床型乳房炎形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0251】
【図32】3つのノルディックレッドウシの14家系の複合解析においてQTLプロファイルを示す体細胞スコア形質についてのQTLプロファイルを示す図である。X軸は、本解析において利用される位置に従ってモルガンで表現される染色体の距離を表す。Y軸は、QTL解析の尤度比検定統計を表す。
【0252】
【図33】フィンランドエアシャーにおける17.8cMでの臨床型乳房炎に対する4マーカーのハプロタイプ効果を示す図である。ハプロタイプ効果がy軸であり、ハプロタイプ数がx軸である。初めは、雌親ハプロタイプの大きいクラスターが存在し、続いて雄親ハプロタイプが存在する。

【特許請求の範囲】
【請求項1】
ウシ被験体における乳房炎に対する耐性を判断する方法であって、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を該ウシ被験体からの試料中で検出することを含み、該少なくとも1つの遺伝的マーカーが、ウシ染色体BTA9上で、多型マイクロサテライトマーカーC6orf93及びinra084に隣接する領域、及び多型マイクロサテライトマーカーC6orf93及びinra084を含む領域、並びに/又はウシ染色体BTA11上で、多型マイクロサテライトマーカーHELMTT43及びBM3501に隣接する領域、及び多型マイクロサテライトマーカーHELMTT43及びBM3501を含む領域に位置し、該少なくとも1つの遺伝的マーカーの有無が、該ウシ被験体又はそれらの出生児の乳房炎耐性を示す、ウシ被験体における乳房炎に対する耐性を判断する方法。
【請求項2】
ウシ被験体における乳房炎に対する耐性を判断する方法であって、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を該ウシ被験体からの試料中で検出することを含み、該少なくとも1つの遺伝的マーカーが、ウシ染色体BTA9上で、多型マイクロサテライトマーカーC6orf93及びinra084に隣接する領域、及び多型マイクロサテライトマーカーC6orf93及びinra084を含む領域に位置し、該少なくとも1つの遺伝的マーカーの有無が該ウシ被験体又はそれらの出生児の乳房炎耐性を示す、ウシ被験体における乳房炎に対する耐性を判断する請求項1に記載の方法。
【請求項3】
前記少なくとも1つの遺伝的マーカーが乳房炎に対する耐性についてのウシの形質と連鎖する、請求項1に記載の方法。
【請求項4】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーC6orf93及びinra084に隣接する領域、及び遺伝的マーカーC6orf93及びinra084を含む領域に位置する、請求項1に記載の方法。
【請求項5】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2251及びinra084に隣接する領域、及び遺伝的マーカーbms2251及びinra084を含む領域に位置する、請求項1に記載の方法。
【請求項6】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbm7234及びinra144に隣接する領域、及び遺伝的マーカーbm7234及びinra144を含む領域に位置する、請求項1に記載の方法。
【請求項7】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2819及びinra144に隣接する領域、及び遺伝的マーカーbms2819及びinra144を含む領域に位置する、請求項1に記載の方法。
【請求項8】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2251及びinra144に隣接する領域、及び遺伝的マーカーbms2251及びinra144を含む領域に位置する、請求項1に記載の方法。
【請求項9】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2819及びinra084に隣接する領域、及び遺伝的マーカーbms2819及びinra084を含む領域に位置する、請求項1に記載の方法。
【請求項10】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbm7234及びbms2819に隣接する領域、及び遺伝的マーカーbm7234及びbms2819を含む領域に位置する、請求項1に記載の方法。
【請求項11】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbm7234及びbm4208に隣接する領域、及び遺伝的マーカーbm7234及びbm4208を含む領域に位置する、請求項1に記載の方法。
【請求項12】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーinra144及びrgs17に隣接する領域、及び遺伝的マーカーinra144及びrgs17を含む領域に位置する、請求項1に記載の方法。
【請求項13】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーinra144及びinra084に隣接する領域、及び遺伝的マーカーinra144及びinra084を含む領域に位置する、請求項1に記載の方法。
【請求項14】
前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーinra084及びrgs17に隣接する領域、及び遺伝的マーカーinra084及びrgs17を含む領域に位置する、請求項1に記載の方法。
【請求項15】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーbm4208及びinra144に隣接する領域、及び遺伝的マーカーbm4208及びinra144を含む領域に位置する、請求項1に記載の方法。
【請求項16】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーbm7234及びinra084に隣接する領域、及び遺伝的マーカーbm7234及びinra084を含む領域に位置する、請求項1に記載の方法。
【請求項17】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーEPM2A及びbm7234に隣接する領域、及び遺伝的マーカーEPM2A及びbm7234を含む領域に位置する、請求項1に記載の方法。
【請求項18】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーbms2251及びbm7234に隣接する領域、及び遺伝的マーカーbms2251及びbm7234を含む領域に位置する、請求項1に記載の方法。
【請求項19】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーC6orf93である、請求項1に記載の方法。
【請求項20】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDIK4986である、請求項1に記載の方法。
【請求項21】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーmm12e6である、請求項1に記載の方法。
【請求項22】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーPEX3である、請求項1に記載の方法。
【請求項23】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDEAD21である、請求項1に記載の方法。
【請求項24】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBMS2251である、請求項1に記載の方法。
【請求項25】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーEPM2Aである、請求項1に記載の方法。
【請求項26】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM7234である、請求項1に記載の方法。
【請求項27】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM4208である、請求項1に記載の方法。
【請求項28】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBMS2819である、請求項1に記載の方法。
【請求項29】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーINRA144である、請求項1に記載の方法。
【請求項30】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーINRA084である、請求項1に記載の方法。
【請求項31】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーrgs17である、請求項1に記載の方法。
【請求項32】
前記少なくとも1つの遺伝的マーカーが遺伝的マーカーの組合せである、請求項1に記載の方法。
【請求項33】
前記少なくとも1つの遺伝的マーカーが、BTA9の69.35cM〜74.5cMの領域に位置する、請求項1に記載の方法。
【請求項34】
前記少なくとも1つの遺伝的マーカーが、BTA9の71.3cM〜74.5cMの領域に位置する、請求項1に記載の方法。
【請求項35】
前記少なくとも1つの遺伝的マーカーが、BTA9の72.3cM〜74.5cMの領域に位置する、請求項1に記載の方法。
【請求項36】
前記少なくとも1つの遺伝的マーカーが、BTA9の73.95cM〜74.5cMの領域に位置する、請求項1に記載の方法。
【請求項37】
マイクロサテライトマーカーC6orf93を増幅するプライマー対が配列番号75及び配列番号76である、請求項19に記載の方法。
【請求項38】
マイクロサテライトマーカーPEX3を増幅するプライマー対が配列番号81及び配列番号82である、請求項22に記載の方法。
【請求項39】
マイクロサテライトマーカーEPM2Aを増幅するプライマー対が配列番号87及び配列番号88である、請求項25に記載の方法。
【請求項40】
マイクロサテライトマーカーrgs17を増幅するプライマー対が配列番号99及び配列番号100である、請求項31に記載の方法。
【請求項41】
ウシ被験体における乳房炎に対する耐性を判断する方法であって、少なくとも1つの遺伝的マーカーが、ウシ染色体11上で、多型マイクロサテライトマーカーHELMTT43及びBM3501に隣接する領域、及び多型マイクロサテライトマーカーHELMTT43及びBM3501を含む領域に位置し、該少なくとも1つの遺伝的マーカーが乳房炎耐性と連鎖し、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を前記被験体からの遺伝子材料中で検出することを包含し、前記少なくとも1つの遺伝的マーカーの存在が乳房炎耐性を示すこと、及び/又は乳房炎耐性を示す出生児を産むことを示す、ウシ被験体における乳房炎に対する耐性を判断する請求項1に記載の方法。
【請求項42】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーHELMTT43及びINRA177に隣接する領域、及び遺伝的マーカーHELMTT43及びINRA177を含む領域に位置する、請求項1に記載の方法。
【請求項43】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーHELMTT43及びMNB−70に隣接する領域、及び遺伝的マーカーHELMTT43及びMNB−70を含む領域に位置する、請求項1に記載の方法。
【請求項44】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーMNB−40及びMNB−70に隣接する領域、及び遺伝的マーカーMNB−40及びMNB−70を含む領域に位置する、請求項1に記載の方法。
【請求項45】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBP38及びINRA131に隣接する領域、及び遺伝的マーカーBP38及びINRA131を含む領域に位置する、請求項1に記載の方法。
【請求項46】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBM2818及びINRA177に隣接する領域、及び遺伝的マーカーBM2818及びINRA177を含む領域に位置する、請求項1に記載の方法。
【請求項47】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBMS1953及びBM2818に隣接する領域、及び遺伝的マーカーBMS1953及びBM2818を含む領域に位置する、請求項1に記載の方法。
【請求項48】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーHELMTT43及びZAP70に隣接する領域、及び遺伝的マーカーHELMTT43及びZAP70を含む領域に位置する、請求項1に記載の方法。
【請求項49】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーZAP70及びIL18RAに隣接する領域、及び遺伝的マーカーZAP70及びIL18RAを含む領域に位置する、請求項1に記載の方法。
【請求項50】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーINRA131及びBM6445に隣接する領域、及び遺伝的マーカーINRA131及びBM6445を含む領域に位置する、請求項1に記載の方法。
【請求項51】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBM304及びBM7169に隣接する領域、及び遺伝的マーカーBM304及びBM7169を含む領域に位置する、請求項1に記載の方法。
【請求項52】
前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBM7169及びDIK5170に隣接する領域、及び遺伝的マーカーBM7169及びDIK5170を含む領域に位置する、請求項1に記載の方法。
【請求項53】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーBM6445及びBMS1048に隣接する領域、及び遺伝的マーカーBM6445及びBMS1048を含む領域に位置する、請求項1に記載の方法。
【請求項54】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーMB110及びBMS2047に隣接する領域、及び遺伝的マーカーMB110及びBMS2047を含む領域に位置する、請求項1に記載の方法。
【請求項55】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーIL18RA及びBM2818に隣接する領域、及び遺伝的マーカーIL18RA及びBM2818を含む領域に位置する、請求項1に記載の方法。
【請求項56】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーBM2818及びBM7169に隣接する領域、及び遺伝的マーカーBM2818及びBM7169を含む領域に位置する、請求項1に記載の方法。
【請求項57】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーBMS2325及びDIK4637に隣接する領域、及び遺伝的マーカーBMS2325及びDIK4637を含む領域に位置する、請求項1に記載の方法。
【請求項58】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーIL18RA及びAUP1に隣接する領域、及び遺伝的マーカーIL18RA及びAUP1を含む領域に位置する、請求項1に記載の方法。
【請求項59】
前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーDIK4637及びUMBTL103に隣接する領域、及び遺伝的マーカーDIK4637及びUMBTL103を含む領域に位置する、請求項1に記載の方法。
【請求項60】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーIL18RAである、請求項1に記載の方法。
【請求項61】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーMNB−40である、請求項1に記載の方法。
【請求項62】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーUAP1である、請求項1に記載の方法。
【請求項63】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDIK4637である、請求項1に記載の方法。
【請求項64】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーUMBTL103である、請求項1に記載の方法。
【請求項65】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーINRA177である、請求項1に記載の方法。
【請求項66】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーUMBTL20である、請求項1に記載の方法。
【請求項67】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM716である、請求項1に記載の方法。
【請求項68】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDIK2653である、請求項1に記載の方法。
【請求項69】
前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM2818である、請求項1に記載の方法。
【請求項70】
前記少なくとも1つの遺伝的マーカーが遺伝的マーカーの一組合せである、請求項1に記載の方法。
【請求項71】
前記少なくとも1つの遺伝的マーカーが、BTA11の2.249cM〜97.223cMの領域に位置する、請求項1に記載の方法。
【請求項72】
前記少なくとも1つの遺伝的マーカーが、BTA11の2.249cM〜35.098cMの領域に位置する、請求項1に記載の方法。
【請求項73】
前記少なくとも1つの遺伝的マーカーが、BTA11の19.440cM〜23.829cMの領域に位置する、請求項1に記載の方法。
【請求項74】
マイクロサテライトマーカーIL18RAを増幅するプライマー対が配列番号125及び配列番号126である、請求項60に記載の方法。
【請求項75】
マイクロサテライトマーカーMNB−40を増幅するプライマー対が配列番号127及び配列番号128である、請求項61に記載の方法。
【請求項76】
マイクロサテライトマーカーDIK4637を増幅するプライマー対が配列番号141及び配列番号142である、請求項63に記載の方法。
【請求項77】
マイクロサテライトマーカーUMBTL103を増幅するプライマー対が配列番号143及び配列番号144である、請求項64に記載の方法。
【請求項78】
前記少なくとも1つの遺伝的マーカーが乳房炎を引き起こす遺伝子と連鎖する、請求項1に記載の方法。
【請求項79】
請求項1の方法によって乳房の健康特性を判断することを含む、育種目的でウシ被験体を選抜する方法。
【請求項80】
乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーのウシ被験体における有無を検出する際に用いる診断キットであって、少なくとも1つのオリゴヌクレオチド配列及びその組合せを包含し、該ヌクレオチド配列が配列番号1〜配列番号192のいずれか、及び/又はそれらの任意の組合せから選択される、乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーのウシ被験体における有無を検出する際に用いる診断キット。

【図1】
image rotate

【図2】
image rotate

【図3】
image rotate

【図4】
image rotate

【図5a】
image rotate

【図5b】
image rotate

【図6a】
image rotate

【図6b】
image rotate

【図7】
image rotate

【図8】
image rotate

【図9】
image rotate

【図10】
image rotate

【図11】
image rotate

【図12】
image rotate

【図13】
image rotate

【図14】
image rotate

【図15】
image rotate

【図16】
image rotate

【図17】
image rotate

【図18】
image rotate

【図19】
image rotate

【図20】
image rotate

【図21】
image rotate

【図22】
image rotate

【図23】
image rotate

【図24】
image rotate

【図25】
image rotate

【図26】
image rotate

【図27】
image rotate

【図28】
image rotate

【図29】
image rotate

【図30】
image rotate

【図31】
image rotate

【図32】
image rotate

【図33】
image rotate


【公表番号】特表2009−525733(P2009−525733A)
【公表日】平成21年7月16日(2009.7.16)
【国際特許分類】
【出願番号】特願2008−553616(P2008−553616)
【出願日】平成19年2月5日(2007.2.5)
【国際出願番号】PCT/DK2007/000058
【国際公開番号】WO2007/090399
【国際公開日】平成19年8月16日(2007.8.16)
【出願人】(508238358)オーフス ユニバーシティ (5)
【氏名又は名称原語表記】AARHUS UNIVERSITET
【出願人】(508238842)
【氏名又は名称原語表記】MTT AGRIFOOD RESEARCH FINLAND
【出願人】(508238277)エストニアン ユニバーシティ オブ ライフ サイエンシズ (1)
【氏名又は名称原語表記】ESTONIAN UNIVERSITY OF LIFE SCIENCES
【Fターム(参考)】