JP2009525733A5 - - Google Patents
Download PDFInfo
- Publication number
- JP2009525733A5 JP2009525733A5 JP2008553616A JP2008553616A JP2009525733A5 JP 2009525733 A5 JP2009525733 A5 JP 2009525733A5 JP 2008553616 A JP2008553616 A JP 2008553616A JP 2008553616 A JP2008553616 A JP 2008553616A JP 2009525733 A5 JP2009525733 A5 JP 2009525733A5
- Authority
- JP
- Japan
- Prior art keywords
- genetic
- marker
- markers
- genetic marker
- microsatellite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000003550 marker Substances 0.000 claims description 384
- 230000002068 genetic Effects 0.000 claims description 374
- 241000283690 Bos taurus Species 0.000 claims description 243
- 210000000349 Chromosomes Anatomy 0.000 claims description 194
- 208000004396 Mastitis Diseases 0.000 claims description 136
- 229920002393 Microsatellite Polymers 0.000 claims description 60
- 102100000256 AUP1 Human genes 0.000 claims description 13
- 108060000731 AUP1 Proteins 0.000 claims description 13
- 108060004491 LTV1 Proteins 0.000 claims description 13
- 102100015408 LTV1 Human genes 0.000 claims description 13
- 229920001850 Nucleic acid sequence Polymers 0.000 claims description 13
- 229920000272 Oligonucleotide Polymers 0.000 claims description 10
- 238000009007 Diagnostic Kit Methods 0.000 claims description 8
- 102100005283 ZAP70 Human genes 0.000 claims description 8
- 101700083482 ZAP70 Proteins 0.000 claims description 8
- 101700081256 EPM2A Proteins 0.000 claims description 7
- 102000033191 EPM2A Human genes 0.000 claims description 6
- 230000001488 breeding Effects 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 5
- 210000000481 Breast Anatomy 0.000 claims description 4
- 101700009067 PEX3 Proteins 0.000 claims description 3
- 102100007525 PEX3 Human genes 0.000 claims 2
- 241000894006 Bacteria Species 0.000 claims 1
- 102100012384 UAP1 Human genes 0.000 claims 1
- 101700002661 UAP1 Proteins 0.000 claims 1
- 238000004458 analytical method Methods 0.000 description 178
- 238000003657 Likelihood-ratio test Methods 0.000 description 44
- 230000000694 effects Effects 0.000 description 42
- 239000002773 nucleotide Substances 0.000 description 29
- 125000003729 nucleotide group Chemical group 0.000 description 29
- 210000001082 somatic cell Anatomy 0.000 description 26
- 229920003013 deoxyribonucleic acid Polymers 0.000 description 24
- 239000011159 matrix material Substances 0.000 description 23
- 210000004080 Milk Anatomy 0.000 description 16
- 239000000654 additive Substances 0.000 description 16
- 230000000996 additive Effects 0.000 description 16
- 235000013336 milk Nutrition 0.000 description 16
- 239000008267 milk Substances 0.000 description 16
- 239000000203 mixture Substances 0.000 description 16
- 239000000523 sample Substances 0.000 description 15
- 230000003234 polygenic Effects 0.000 description 13
- 238000000611 regression analysis Methods 0.000 description 13
- 238000000034 method Methods 0.000 description 10
- 108020004707 nucleic acids Proteins 0.000 description 10
- 150000007523 nucleic acids Chemical class 0.000 description 10
- 230000000875 corresponding Effects 0.000 description 8
- 235000013365 dairy product Nutrition 0.000 description 8
- 230000003321 amplification Effects 0.000 description 7
- 238000009396 hybridization Methods 0.000 description 7
- 238000003199 nucleic acid amplification method Methods 0.000 description 7
- FAPWRFPIFSIZLT-UHFFFAOYSA-M sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 7
- 230000002759 chromosomal Effects 0.000 description 6
- 230000000295 complement Effects 0.000 description 6
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 6
- 241000042032 Petrocephalus catostoma Species 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 239000002131 composite material Substances 0.000 description 5
- 201000010099 disease Diseases 0.000 description 5
- 239000000047 product Substances 0.000 description 5
- 238000005204 segregation Methods 0.000 description 5
- 229920000160 (ribonucleotides)n+m Polymers 0.000 description 4
- 238000007476 Maximum Likelihood Methods 0.000 description 4
- 210000000582 Semen Anatomy 0.000 description 4
- DBMJMQXJHONAFJ-UHFFFAOYSA-M Sodium laurylsulphate Chemical compound [Na+].CCCCCCCCCCCCOS([O-])(=O)=O DBMJMQXJHONAFJ-UHFFFAOYSA-M 0.000 description 4
- 238000004581 coalescence Methods 0.000 description 4
- 238000007405 data analysis Methods 0.000 description 4
- 230000000337 effect on mastitis Effects 0.000 description 4
- 230000004983 pleiotropic Effects 0.000 description 4
- 241000894007 species Species 0.000 description 4
- 238000007619 statistical method Methods 0.000 description 4
- 108020005187 Oligonucleotide Probes Proteins 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 210000004027 cells Anatomy 0.000 description 3
- 230000036425 denaturation Effects 0.000 description 3
- 238000004925 denaturation Methods 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- ZHNUHDYFZUAESO-UHFFFAOYSA-N formamide Chemical compound NC=O ZHNUHDYFZUAESO-UHFFFAOYSA-N 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000002751 oligonucleotide probe Substances 0.000 description 3
- 239000011780 sodium chloride Substances 0.000 description 3
- 239000010902 straw Substances 0.000 description 3
- 239000006228 supernatant Substances 0.000 description 3
- 102100009763 MAP4K4 Human genes 0.000 description 2
- 101710038965 MAP4K4 Proteins 0.000 description 2
- 210000004293 Mammary Glands, Human Anatomy 0.000 description 2
- TWRXJAOTZQYOKJ-UHFFFAOYSA-L MgCl2 Chemical compound [Mg+2].[Cl-].[Cl-] TWRXJAOTZQYOKJ-UHFFFAOYSA-L 0.000 description 2
- 229920002224 Peptide nucleic acid Polymers 0.000 description 2
- 238000010171 animal model Methods 0.000 description 2
- 238000000137 annealing Methods 0.000 description 2
- 150000001768 cations Chemical class 0.000 description 2
- 230000003247 decreasing Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 201000009910 diseases by infectious agent Diseases 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 244000144972 livestock Species 0.000 description 2
- 238000007403 mPCR Methods 0.000 description 2
- 235000013372 meat Nutrition 0.000 description 2
- 238000002844 melting Methods 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 239000008188 pellet Substances 0.000 description 2
- 238000001558 permutation test Methods 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- LENZDBCJOHFCAS-UHFFFAOYSA-N 2-amino-2-(hydroxymethyl)propane-1,3-diol Chemical compound OCC(N)(CO)CO LENZDBCJOHFCAS-UHFFFAOYSA-N 0.000 description 1
- 108060002080 ADAT2 Proteins 0.000 description 1
- 102100001573 AIG1 Human genes 0.000 description 1
- 108060000247 AIG1 Proteins 0.000 description 1
- 229920000936 Agarose Polymers 0.000 description 1
- 229920002287 Amplicon Polymers 0.000 description 1
- 210000004369 Blood Anatomy 0.000 description 1
- 206010006242 Breast enlargement Diseases 0.000 description 1
- 206010006262 Breast inflammation Diseases 0.000 description 1
- 101700026049 CCN2 Proteins 0.000 description 1
- 101710015494 CITED2 Proteins 0.000 description 1
- 102100010290 CITED2 Human genes 0.000 description 1
- 241001466804 Carnivora Species 0.000 description 1
- 241000699800 Cricetinae Species 0.000 description 1
- 238000007400 DNA extraction Methods 0.000 description 1
- 239000003155 DNA primer Substances 0.000 description 1
- 101700054181 ESR1 Proteins 0.000 description 1
- 101710040068 EZR Proteins 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- ZMMJGEGLRURXTF-UHFFFAOYSA-N Ethidium bromide Chemical compound [Br-].C12=CC(N)=CC=C2C2=CC=C(N)C=C2[N+](CC)=C1C1=CC=CC=C1 ZMMJGEGLRURXTF-UHFFFAOYSA-N 0.000 description 1
- 229920002760 Expressed sequence tag Polymers 0.000 description 1
- 210000004907 Glands Anatomy 0.000 description 1
- 101710030947 IFNGR1 Proteins 0.000 description 1
- 206010022114 Injury Diseases 0.000 description 1
- 210000000265 Leukocytes Anatomy 0.000 description 1
- 210000002264 Mammary Glands, Animal Anatomy 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 210000000282 Nails Anatomy 0.000 description 1
- 108010059712 Pronase Proteins 0.000 description 1
- 101700022364 RGS17 Proteins 0.000 description 1
- 101710027587 ROS1 Proteins 0.000 description 1
- -1 RP1-172K10 Proteins 0.000 description 1
- 101700056561 SHPRH Proteins 0.000 description 1
- 210000000952 Spleen Anatomy 0.000 description 1
- 108010006785 Taq Polymerase Proteins 0.000 description 1
- FQENQNTWSFEDLI-UHFFFAOYSA-J Tetrasodium pyrophosphate Chemical compound [Na+].[Na+].[Na+].[Na+].[O-]P([O-])(=O)OP([O-])([O-])=O FQENQNTWSFEDLI-UHFFFAOYSA-J 0.000 description 1
- 239000007984 Tris EDTA buffer Substances 0.000 description 1
- HRXKRNGNAMMEHJ-UHFFFAOYSA-K Trisodium citrate Chemical compound [Na+].[Na+].[Na+].[O-]C(=O)CC(O)(CC([O-])=O)C([O-])=O HRXKRNGNAMMEHJ-UHFFFAOYSA-K 0.000 description 1
- 101700003320 VIP Proteins 0.000 description 1
- QQEHDZXXCDSAFE-JBSAMAPISA-N [(3S,8R,9S,10R,13S,14S,17R)-17-acetyl-6-chloro-3-hydroxy-10,13-dimethyl-1,2,3,8,9,11,12,14,15,16-decahydrocyclopenta[a]phenanthren-17-yl] acetate Chemical compound C1=C(Cl)C2=C[C@@H](O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@@](C(C)=O)(OC(=O)C)[C@@]1(C)CC2 QQEHDZXXCDSAFE-JBSAMAPISA-N 0.000 description 1
- 230000001058 adult Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- KCXVZYZYPLLWCC-UHFFFAOYSA-N edta Chemical compound OC(=O)CN(CC(O)=O)CCN(CC(O)=O)CC(O)=O KCXVZYZYPLLWCC-UHFFFAOYSA-N 0.000 description 1
- 230000029578 entry into host Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001502 gel electrophoresis Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000000813 microbial Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000037125 natural defense Effects 0.000 description 1
- 238000003752 polymerase chain reaction Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 239000011535 reaction buffer Substances 0.000 description 1
- 238000007894 restriction fragment length polymorphism technique Methods 0.000 description 1
- 239000001509 sodium citrate Substances 0.000 description 1
- 229940048086 sodium pyrophosphate Drugs 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000000392 somatic Effects 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 235000019818 tetrasodium diphosphate Nutrition 0.000 description 1
- 239000001577 tetrasodium phosphonato phosphate Substances 0.000 description 1
- 238000010257 thawing Methods 0.000 description 1
- 210000001519 tissues Anatomy 0.000 description 1
- 239000011778 trisodium citrate Substances 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Description
本発明は、ウシ被験体における乳房炎に対する耐性を判断する方法であって、ウシ染色体BTA9及びBTA11上に位置する少なくとも1つの遺伝的マーカーを検出することを含む、ウシ被験体における乳房炎に対する耐性を判断する方法に関する。さらに、本発明は、乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーの有無を検出する診断キットに関する。
乳房炎は、感染又は外傷から生じる雌ウシ(cow)の乳腺(mammary gland)又は乳房の炎症であり、乳房炎はウシにおける最も経済的に重要な疾患であると考えられている。本疾患は多様な物質によって生じるおそれがある。乳房炎の主な原因は、乳頭先端を介した微生物の乳腺侵入である。
乳房炎は臨床型又は無症候性(sub-clinical)である可能性があり、無症候性の感染が臨床型の症状発現に先立つ。臨床型乳房炎(clinical mastitis)(CM)は、赤く且つ腫大した(swollen)乳腺、即ち、赤く腫大した乳房、凝結乳(clotted milk)の生産を通じて目視で検出され得る。一度検出されると、乳房炎の雌ウシ由来の乳は、全体の乳品質に影響を及ぼさないように、タンクから分離して保管される。
無症候性乳房炎は、高い体細胞数(somatic cell counts)(SCS)、正常な又は上昇した体温、及び培養物の検査が陽性を示す乳試料を特徴とする乳房炎の一種である。したがって、無症候性乳房炎は目視で、乳房の腫大、又は腺若しくは産出される乳の観察によって検出することはできない。このため、農家には、無症候性乳房炎の雌ウシ由来の乳を他へそらすという選択肢がない。しかしながら、この乳は、非感染の雌ウシの乳よりも品質が劣っているため、タンク中の乳の残りを汚染するおそれがある。
乳房炎は、体細胞数(SCS)を用いて検出され得る。ここでは、雌ウシ由来の乳の試料が体細胞(白血球)の存在について分析される。体細胞は、雌ウシの天然の防御機構の一部であり、乳房が感染すると細胞数は増加する。乳試料中の体細胞の数は、転球式粘度計(rolling-ball viscometer)及びクールター計数器によって間接的に推定され得る。
乳房炎により、乳及び乳製品の量及び質が低下するので、農家及び乳業に経済的損失が生じる。したがって、ウシにおける乳房炎に対する耐性の遺伝的根拠を判断することが可能なことは、毎日の乳生産の観点だけではなく、乳房炎に対する耐性を有するウシ被験体を選抜する育種管理の観点においても、共に、乳業に対する非常に大きな経済的重要性がある。より乳房炎になりにくい雌ウシを産生する、耐性が改善したウシ被験体を遺伝的に選抜する方法が望ましい。
遺伝的形質(genetic traits)に対する遺伝的決定因子(genetic determinants)を同定するための1つのアプローチは、過去の組換えの活用を目的とした連鎖不平衡(linkage
disequilibrium)(LD)マッピングの使用であり、ヒト集団と比較して非常に長い染色体区域に渡って広がることが乳牛を含む幾つかの家畜(livestock)集団において示されている、(Farnir et al., 2000)。一度マッピングされると、量的形質遺伝子座(Quantitative Trait Locus)(QTL)は、マーカー利用選抜(marker assisted selection)において有益に利用され得る。
disequilibrium)(LD)マッピングの使用であり、ヒト集団と比較して非常に長い染色体区域に渡って広がることが乳牛を含む幾つかの家畜(livestock)集団において示されている、(Farnir et al., 2000)。一度マッピングされると、量的形質遺伝子座(Quantitative Trait Locus)(QTL)は、マーカー利用選抜(marker assisted selection)において有益に利用され得る。
連鎖不平衡
連鎖不均衡は、過去を遡る組換え事象(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)。高い信頼度をもって興味のある特定の領域のマップを使用するために、連鎖不均衡の存在が考慮される。
連鎖不均衡は、過去を遡る組換え事象(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)。高い信頼度をもって興味のある特定の領域のマップを使用するために、連鎖不均衡の存在が考慮される。
本発明において、臨床型乳房炎及び/又はSCSに関連した量的形質遺伝子座はウシ染色体BTA9で同定されており、ウシ被験体が乳房炎に対して耐性があるか否かを判断する方法を可能にする。
乳房炎に対する耐性が増大したウシ被験体を選抜し、それにより乳房炎を患う動物に関連した経済的損失を避けることが可能になることは、ウシ産業内で経済的関心が大きい。乳房炎に対する耐性に関する遺伝的素因(genetic predisposition)は、本発明によって検出され得る。本発明は、乳房炎に対する耐性に関連し且つ/又はそれと連鎖する遺伝的マーカーに基づき、ウシ被験体において乳房炎に対する耐性を判断する方法を提供する。
したがって、本発明の一態様は、ウシ被験体における乳房炎に対する耐性を判断するための方法であって、当該ウシ被験体からの試料中に、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を検出することを含み、当該少なくとも1つの遺伝的マーカーが、ウシ染色体BTA9上で、多型マイクロサテライトマーカーC6orf93及びinra084に挟まれ、且つ、多型マイクロサテライトマーカーC6orf93及びinra084を含む領域、並びに/又はウシ染色体BTA11上で、多型マイクロサテライトマーカーHELMTT43及びBM3501に挟まれ、且つ、多型マイクロサテライトマーカーHELMTT43及びBM3501を含む領域に位置し、当該少なくとも1つの遺伝的マーカーの有無が当該ウシ被験体又はその出生児の乳房炎耐性を示す方法に関する。
本発明の第二の態様は、ウシ被験体において乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーの有無を検出する際に使用するための診断キットであって、少なくとも1つのオリゴヌクレオチド配列及びそれらの組合せを含み、当該ヌクレオチド配列が配列番号1〜配列番号192のいずれか及び/又はそれらの任意の組合せから選択される、診断キットに関する。
本発明は、乳牛における乳房炎耐性の遺伝的決定因子に関する。乳房炎、臨床型乳房炎及び無症候性乳房炎の発生は共に、乳業に対する実質的な経済的損失に関与する。したがって、乳房炎耐性に関する遺伝的素因を有するそれらのウシ被験体を同定することには、経済的な関心がある。かかる遺伝的素因を有するウシ被験体は、出生児に伝播し得る望ましい形質のキャリア(carriers)である。
「ウシ被験体」という用語は任意の品種のウシを指し、成体動物であれ新生動物であれ、雌ウシ(cows)及び雄ウシ(bulls)の両方を含む意味である。この用語によって特定の年齢の動物が示されることはない。ウシ被験体の一例は、ホルスタイン種のひとつである。好ましい一実施形態において、ウシ被験体は、ホルスタインフリーシャンウシ集団のひとつである。別の実施形態において、ウシ被験体は、ホルスタインスワルトボント(Swartbont)ウシ集団のひとつである。別の実施形態において、ウシ被験体は、ドイツホルスタインシュワルツブントウシ集団のひとつである。別の実施形態において、ウシ被験体は、米国ホルスタインウシ集団のひとつである。一実施形態において、ウシ被験体は、レッドアンドホワイトホルスタイン種のひとつである。別の実施形態において、ウシ被験体は、ドイツホルスタインシュワルツブントウシ集団のひとつである。一実施形態において、ウシ被験体は、ホルスタイン種のひとつを含む、任意の家系のひとつである。好ましい一実施形態において、ウシ被験体は、デンマークレッド集団のひとつである。別の好ましい実施形態において、ウシ被験体は、フィンランドエアシャー集団のひとつである。さらに別の実施形態において、ウシ被験体は、スウェーデンレッドアンドホワイト(Swedish Red and White)集団のひとつである。さらなる実施形態において、ウシ被験体は、デンマークホルスタイン集団のひとつである。別の実施形態において、ウシ被験体は、スウェーデンレッドアンドホワイト集団のひとつである。さらに別の実施形態において、ウシ被験体は、ノルディックレッド(Nordic Red)集団のひとつである。
本発明の一実施形態において、ウシ被験体は、スウェーデンレッドアンドホワイト、デンマークレッド、フィンランドエアシャー、ホルスタインフリーシャン、デンマークホルスタイン及びノルディックレッドから成る群から選択される。本発明の別の実施形態において、ウシ被験体は、フィンランドエアシャー及びスウェーデンレッドアンドホワイトウシから成る群から選択される。本発明の別の実施形態において、ウシ被験体は、フィンランドエアシャー及びスウェーデンレッドアンドホワイトウシから成る群から選択される。
一実施形態において、ウシ被験体は表1aに示される品種の群から選択される。
一実施形態において、ウシ被験体は表1bに示される品種の群から選択される品種のひとつである。
一実施形態において、ウシ被験体は表1cに示される品種の群から選択される品種のひとつである。
「遺伝的マーカー」という用語は、ウシ染色体上のDNAの可変ヌクレオチド配列(variable nucleotide sequence)(多型(polymorphism))を指し、アレル同士を判別する。可変ヌクレオチド配列は、例えば増幅法(amplification methods)及び/又はサイズの差異の観察において特定のオリゴヌクレオチドを用いることにより、当業者に既知の方法によって同定することができる。しかしながら、可変ヌクレオチド配列は、シーケンシングにより、又は例えば制限断片長多型(restriction fragment length polymorphism)解析によっても検出され得る。可変ヌクレオチド配列は、欠失、挿入、反復及び/又は点変異(point mutation)によって表され得る。
遺伝的マーカーの1つの型は、量的形質遺伝子座と連鎖するマイクロサテライトマーカーである。マイクロサテライトマーカーは、互いの後で反復される短い配列を指す。短い配列は、例えば1ヌクレオチド、例えば2ヌクレオチド、例えば3ヌクレオチド、例えば4ヌクレオチド、例えば5ヌクレオチド、例えば6ヌクレオチド、例えば7ヌクレオチド、例えば8ヌクレオチド、例えば9ヌクレオチド、例えば10ヌクレオチドである。しかしながら、変化が生じることがあり、反復の数は増加又は減少し得る。多型マイクロサテライトマーカーの特定の定義及び遺伝子座は、USDA遺伝子地図(genetic map)に見出され得る(Kappes et al., 1997;又は米国食肉動物研究所(http://www.marc.usda.gov/)へのリンクをたどることによる)。
さらに、本発明の遺伝的マーカーのヌクレオチド配列は、ウシ被験体における乳房炎耐性に関する形質と遺伝的に連鎖する(genetically linked)と理解される。結果的に、本発明の方法により、複数の遺伝的マーカーに挟まれ、且つ、当該遺伝的マーカーを含むDNA領域のヌクレオチド配列から、多数の遺伝的マーカーが生成され得ることも理解される。
「量的形質遺伝子座(QTL)」という用語は、特定形質(例えば、草丈)に関連するDNA領域である。必ずしもそれ自体が遺伝子である必要はないが、QTLは、問題の形質の基となる遺伝子に密接に連鎖する一続きの(stretches of)DNAである。
「乳房炎」という用語は、本願において例えば高い体細胞スコア(SCS)を特徴とする無症候性乳房炎、及び臨床型乳房炎の両方を示すのに用いられる。
「乳房炎耐性」及び「乳房炎に対する耐性」という用語は交換して用いることが可能であり、或るウシ被験体が他のウシ被験体より乳房炎になりにくいということに関係する。乳房炎に対する耐性に関連する遺伝的マーカーを判断するため本発明の如く多数のウシ被験体の解析を行う場合には、解析したウシ被験体における臨床型乳房炎及び/又は無症候性乳房炎の発生と連鎖した遺伝的マーカーの有無によって、乳房炎に対する耐性をほのめかす形質が観察され得る。乳房炎耐性には、ウシ被験体又はその出生児における乳房の健康に影響する形質に対する耐性、を包含することが理解される。したがって、雄ウシの乳房炎耐性は、その雌出生児によって身体的に明らかにされる。
乳房炎耐性についてのスコア化
雄ウシの娘を、乳房炎耐性及び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。
雄ウシの娘を、乳房炎耐性及び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。
本発明の一実施形態において、本明細書中に記載の方法及びキットは、乳房炎耐性に関する。本発明の別の実施形態において、本明細書中に記載の方法及びキットは、臨床型乳房炎に対する耐性に関する。別の実施形態において、本発明の方法及びキットは、体細胞数等によって検出される無症候性乳房炎に対する耐性に関する。さらに別の実施形態において、本発明の方法及びキットは主に、例えば、体細胞数によって検出される無症候性乳房炎に対する耐性と組み合わされた臨床型乳房炎に対する耐性、に関する。
試料
本発明による方法はウシ被験体の試料を分析することを含む。ここで、当該試料は、当該方法において使用するためのウシ遺伝子材料を提供可能な任意の好適な試料であり得る。ウシ遺伝子材料は、例えば、血液試料、組織試料(例えば脾臓(spleen)、頬側塗抹(buccal smears))、体表(body surface)(毛又は爪)の切取り部、乳及び/又は精液(semen)から抽出、単離され、必要に応じて精製され得る。試料は新鮮でもよく、凍結してもよい。
本発明による方法はウシ被験体の試料を分析することを含む。ここで、当該試料は、当該方法において使用するためのウシ遺伝子材料を提供可能な任意の好適な試料であり得る。ウシ遺伝子材料は、例えば、血液試料、組織試料(例えば脾臓(spleen)、頬側塗抹(buccal smears))、体表(body surface)(毛又は爪)の切取り部、乳及び/又は精液(semen)から抽出、単離され、必要に応じて精製され得る。試料は新鮮でもよく、凍結してもよい。
本発明の配列多型は、少なくとも1つのヌクレオチド差異、例えば少なくとも2つのヌクレオチド差異、例えば少なくとも3つのヌクレオチド差異、例えば少なくとも4つのヌクレオチド差異、例えば少なくとも5つのヌクレオチド差異、例えば少なくとも6つのヌクレオチド差異、例えば少なくとも7つのヌクレオチド差異、例えば少なくとも8つのヌクレオチド差異、例えば少なくとも9つのヌクレオチド差異、例えば10個のヌクレオチド差異を含む。ヌクレオチド差異は、ヌクレオチド差異、欠失及び/又は挿入、或いはそれらの任意の組合せを含む。
グランドドーターデザイン
グランドドーターデザイン(grand daughter design)は、育種において広範に使用されている祖父についての、及び祖父の息子が出生児を出産した場合には当該息子についての、DNAベースのマーカーからデータを分析することを含む。DNAマーカーデータと共に使用されることとなる表現型データは、当該息子の娘由来である。かかる表現型データは、例えば乳生産の特徴、出産に関する特徴、肉質又は疾患であり得る。娘の一群は、その父親から一方のアレルを受け継ぐ一方、別の群の娘は、その父親から他方のアレルを受け継いでいる。2つの群からのデータを比較することにより、特定の染色体の断片が問題の形質に影響を及ぼす1つ又は複数の遺伝子を保有しているかどうかについての情報を得ることができる。QTLが染色体のこの断片内に存在するかどうかを結論付けることができる。
グランドドーターデザイン(grand daughter design)は、育種において広範に使用されている祖父についての、及び祖父の息子が出生児を出産した場合には当該息子についての、DNAベースのマーカーからデータを分析することを含む。DNAマーカーデータと共に使用されることとなる表現型データは、当該息子の娘由来である。かかる表現型データは、例えば乳生産の特徴、出産に関する特徴、肉質又は疾患であり得る。娘の一群は、その父親から一方のアレルを受け継ぐ一方、別の群の娘は、その父親から他方のアレルを受け継いでいる。2つの群からのデータを比較することにより、特定の染色体の断片が問題の形質に影響を及ぼす1つ又は複数の遺伝子を保有しているかどうかについての情報を得ることができる。QTLが染色体のこの断片内に存在するかどうかを結論付けることができる。
グランドドーターデザインを実行するための前提条件は、詳細な表現型データが利用可能であることである。本発明において、かかるデータが入手可能となっている(http://www.lr.dk/kvaeg/diverse/principles.pdf)。
QTLは、量的形質遺伝子座の短縮形である。個体に量的形質(quantitative traits)を与える遺伝子は、あたかも1つ又は複数の遺伝子が位置するように作用(act)する染色体の一片を観察することによって、間接的に見出され得る。
対照的に、ある染色体上の多数のDNAマーカーが一方又は両方の親及びそれらの出生児について決定された場合、親からその出生児のうちの1頭又は複数頭へと伝えられた形質の情報を提供するために、DNAマーカーが直接使用され得る。このマーカーを使用して、DNAマーカーと連鎖した染色体の遺伝的履歴(genetic history)を算出し得る。
染色体領域及びマーカー
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個の遺伝的マーカーである。
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個の遺伝的マーカーである。
少なくとも1つの遺伝的マーカーは、少なくとも1つのウシ染色体、例えば2つの染色体、例えば3つの染色体、例えば4つの染色体、例えば5つの染色体、及び/又は例えば6つの染色体上に位置し得る。少なくとも1つの遺伝的マーカーがウシ染色体9上に位置し得る。しかしながら、当該少なくとも1つの遺伝的マーカーは、異なる染色体上に位置するマーカーの組合せでもよい。
少なくとも1つのマーカーは、本明細書中に示される表の個々のマーカーのいずれかから選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーc6orf93及びrgs17に挟まれ、且つ、マーカーc6orf93及びrgs17を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約69.35cM〜約79.8cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表1dで示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上に位置する。一実施形態では、ウシ染色体BTA9上で、少なくとも1つの遺伝的マーカーは、マーカーc6orf93及びinra084に挟まれ、且つ、マーカーc6orf93及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約69.35cM〜約74.5cM(この分析に用いられる位置に従って)の領域に位置する。MARCマーカー地図に従って、遺伝的マーカーinra084の位置は90.98である。少なくとも1つの遺伝的マーカーは、表2で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2251及びinra084に挟まれ、且つ、マーカーbms2251及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約71.3cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約86.58cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表3で示されるマーカーの群から選択される。
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2251及びinra144に挟まれ、且つ、マーカーbms2251及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約71.3cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約86.58cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表4で示されるマーカーの群から選択される。
さらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びinra084に挟まれ、且つ、マーカーbm7234及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表5で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びinra144に挟まれ、且つ、マーカーbm7234及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表6で示されるマーカーの群から選択される。
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びbms2819に挟まれ、且つ、マーカーbm7234及びbms2819を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約73.95cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表7で示されるマーカーの群から選択される。
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm7234及びbm4208に挟まれ、且つ、マーカーbm7234及びbm4208を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.3cM〜約73.9cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約88.136cM〜約90.69cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表8で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2819及びinra144に挟まれ、且つ、マーカーbms2819及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約73.95cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.98cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表9で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2819及びinra084に挟まれ、且つ、マーカーbms2819及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約73.95cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.98cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表10で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbm4208及びinra144に挟まれ、且つ、マーカーbm4208及びinra144を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約73.9cM〜約74.2cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.69cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表11で示されるマーカーの群から選択される。
本発明のさらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーinra144及びinra084に挟まれ、且つ、マーカーinra144及びinra084を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約74.2cM〜約74.5cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約90.98cM〜約90.98cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表12で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーbms2251及びbm7234に挟まれ、且つ、マーカーbms2251及びbm7234を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約71.3cM〜約72.3cM(この分析に用いられる位置に従って)、及びMARCマーカー地図に従って約86.58cM〜約88.136cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表13で示されるマーカーの群から選択される。
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーEPM2A及びbm7234に挟まれ、且つ、マーカーEPM2A及びbm7234を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約72.1cM〜約72.3cM(この分析に用いられる位置に従って)の領域、及びMARCマーカー地図に従ってbm7234については約88.136cMに位置する。少なくとも1つの遺伝的マーカーは、表14で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA9上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーinra144及びrgs17に挟まれ、且つ、マーカーinra144及びrgs17を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約74.2cM〜約79.8cM(この分析に用いられる位置に従って)の領域に位置し、MARCマーカー地図に従ったinra144の位置は90.98cMである。少なくとも1つの遺伝的マーカーは、表15で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA9上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上で、マーカーinra084及びrgs17に挟まれ、且つ、マーカーinra084及びrgs17を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA9上の約74.5cM〜約79.8cM(この分析に用いられる位置に従って)の領域に位置し、MARCマーカー地図に従ったinra084の位置は90.68cMである。少なくとも1つの遺伝的マーカーは、表16で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びBM3501に挟まれ、且つ、マーカーHELMTT43及びBM3501を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM〜約97.223cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表17で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーはウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びINRA177に挟まれ、且つ、マーカーHELMTT43及びINRA177を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM〜約35.098cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表18で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びMNB-70に挟まれ、且つ、マーカーHELMTT43及びMNB-70を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM〜約24.617cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表19で示されるマーカーの群から選択される。
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB-40及びMNB-70に挟まれ、且つ、マーカーMNB-40及びMNB-70を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約24.617cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表20で示されるマーカーの群から選択される。
さらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBP38及びINRA131に挟まれ、且つ、マーカーBP38及びINRA131を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約24.617cM〜約47.289cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表21で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM2818及びINRA177に挟まれ、且つ、マーカーBM2818及びINRA177を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約30.009cM〜約35.098cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表22で示されるマーカーの群から選択される。
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBMS1953及びBM2818に挟まれ、且つ、マーカーBMS1953及びBM2818を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約21.537cM〜約30.009cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表23で示されるマーカーの群から選択される。
別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーHELMTT43及びZAP70に挟まれ、且つ、マーカーHELMTT43及びZAP70を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約2.249cM(MARCマーカー地図に従って)〜約5.4cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表24で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーZAP70及びIL18RAに挟まれ、且つ、マーカーZAP70及びIL18RAを含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、(この分析に用いられる位置に従って)ウシ染色体BTA11上の約5.4cM〜約12.3cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表25で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーINRA131及びBM6445に挟まれ、且つ、マーカーINRA131及びBM6445を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約47.289cM〜約61.570cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表26で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM304及びBM7169に挟まれ、且つ、マーカーBM304及びBM7169を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約33.597cM〜約50.312cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表27で示されるマーカーの群から選択される。
本発明のさらに別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM7169及びDIK5170に挟まれ、且つ、マーカーBM7169及びDIK5170を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約50.312cM〜約70.143cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表28で示されるマーカーの群から選択される。
さらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM6445及びBMS1048に挟まれ、且つ、マーカーBM6445及びBMS1048を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約61.570cM〜約81.065cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表29で示されるマーカーの群から選択される。
またさらなる実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMB110及びBMS2047に挟まれ、且つ、マーカーMB110及びBMS2047を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約68.679cM〜約78.457cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表30で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びBM2818に挟まれ、且つ、マーカーIL18RA及びBM2818を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約30.009cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表31で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM2818及びBM7169に挟まれ、且つ、マーカーBM2818及びBM7169を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約30.009cM〜約50.312cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表32で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMAP4K4及びBM2818に挟まれ、且つ、マーカーMAP4K4及びBM2818を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約10.5cM(この分析に用いられる位置に従って)〜約30.009cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表33で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びUMBTL103に挟まれ、且つ、マーカーIL18RA及びUMBTL103を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約23.829cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表34で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB-40及びDIK2653に挟まれ、且つ、マーカーMNB-40及びDIK2653を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約20.135cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表35で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM716及びDIK4637に挟まれ、且つ、マーカーBM716及びDIK4637を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約22.527cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表36で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBM716及びBMS2569に挟まれ、且つ、マーカーBM716及びBMS2569を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM〜約21.082cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表37で示されるマーカーの群から選択される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーBMS2325及びDIK4637に挟まれ、且つ、マーカーBMS2325及びDIK4637を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約21.082cM〜約22.527cMの領域に位置する。少なくとも1つの遺伝的マーカーは、表38で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びAUP1に挟まれ、且つ、マーカーIL18RA及びAUP1を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約17.6cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表39で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RA及びMNB-40に挟まれ、且つ、マーカーIL18RA及びMNB-40を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)〜約19.440cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表40で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB-40及びAUP1に挟まれ、且つ、マーカーMNB-40及びAUP1を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約19.440cM(MARCマーカー地図に従って)〜約17.6cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは、表41で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーIL18RAを含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約12.3cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表42で示される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーMNB-40を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約16.0cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表43で示される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーAUP1を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約17.6cM(この分析に用いられる位置に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表44で示される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーDIK4637及びUMBTL103に挟まれ、且つ、マーカーDIK4637及びUMBTL103を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約22.527cM〜約23.829cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表45で示されるマーカーの群から選択される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーDIK4637を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約22.527cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表46で示される。
本発明の一実施形態において、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上に位置する。一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上で、マーカーUMBTL103を含む領域に位置する。本発明の一実施形態では、少なくとも1つの遺伝的マーカーは、ウシ染色体BTA11上の約23.829cM(MARCマーカー地図に従って)の領域に位置する。少なくとも1つの遺伝的マーカーは表47で示される。
本発明の別の実施形態において、少なくとも1つの遺伝的マーカーは、本明細書の他の箇所で記載されるように、BTA9の任意の領域及びマーカーがBTA11の任意の領域及びマーカーと組み合わされるマーカーの組合せである。
プライマー
表50に、本発明により使用され得るプライマーを示す。表50で特定されるプライマー対は個々に使用してもよく、表50の1つ又は複数のプライマー対と併用してもよい。かかるプライマー又はプローブの設計は、当業者である分子生物学者には明らかであろう。かかるプライマーは、任意の好都合な長さ、例えば最大50塩基長、最大40塩基長、さらに好都合には最大30塩基長、例えば8〜25塩基長又は8〜15塩基長を有する。概して、かかるプライマーは、領域内の対応する野生型又はバリアント遺伝子座(variant locus)と完全に相補的な塩基配列(base sequences)を含むであろう。しかしながら、必要に応じて、過度にオリゴヌクレオチドプローブの判別力(discriminatory power)に影響を及ぼさないという条件で、1つ又は複数ののミスマッチ(mismatches)を導入してもよい。本発明のプライマー/プローブは、検出を容易にするための1つ又は複数の標識を保有し得る。
表50に、本発明により使用され得るプライマーを示す。表50で特定されるプライマー対は個々に使用してもよく、表50の1つ又は複数のプライマー対と併用してもよい。かかるプライマー又はプローブの設計は、当業者である分子生物学者には明らかであろう。かかるプライマーは、任意の好都合な長さ、例えば最大50塩基長、最大40塩基長、さらに好都合には最大30塩基長、例えば8〜25塩基長又は8〜15塩基長を有する。概して、かかるプライマーは、領域内の対応する野生型又はバリアント遺伝子座(variant locus)と完全に相補的な塩基配列(base sequences)を含むであろう。しかしながら、必要に応じて、過度にオリゴヌクレオチドプローブの判別力(discriminatory power)に影響を及ぼさないという条件で、1つ又は複数ののミスマッチ(mismatches)を導入してもよい。本発明のプライマー/プローブは、検出を容易にするための1つ又は複数の標識を保有し得る。
一実施形態において、プライマー及び/又はプローブは、本明細書中に示されるようなマーカーによって描写される配列を含有する一塩基多型(single nucleotide polymorphism)とハイブリッド形成するサブシーケンス(subsequence)とのハイブリダイゼーション、及び/又はその増幅が可能である。
本発明のプライマーヌクレオチド配列はさらに以下を含む:(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)を含むものである。
本発明の核酸分子のうちの一つは、高ストリンジェント条件又はストリンジェント条件下で上記の核酸分子とハイブリッド形成する、デオキシオリゴヌクレオチド(「オリゴ(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℃低い温度で実行される。
例示的な高ストリンジェント条件は、例えば、37℃(約14塩基オリゴについて)、48℃(約17塩基オリゴについて)、55℃(約20塩基オリゴについて)、及び60℃(約23塩基オリゴについて)での、6XSSC/0.05%ピロリン酸(pyrophosphate)ナトリウム中での洗浄を指す可能性がある。
したがって、本発明はさらに、本発明の多型を検出するヌクレオチドプライマー又はプローブを提供する。評価は、少なくとも1つの核酸プライマー又はプローブの手段、例えばDNA、RNA又はペプチド核酸(peptide nucleic acid, PNA)若しくはロックド核酸(locked nucleic acid, LNA)等の核酸アナログ(analogue)のプライマー又はプローブによって実行され得る。
本発明の一態様によれば、描写領域中の1つ又は複数の位置で多型を検出することが可能なアレル特異的オリゴヌクレオチドプローブが提供される。
アレル特異的オリゴヌクレオチドプローブは、好ましくは5〜50ヌクレオチド、より好ましくは約5〜35ヌクレオチド、より好ましくは約5〜30ヌクレオチド、より好ましくは少なくとも9ヌクレオチドである。
連鎖の決定
遺伝的マーカーが遺伝子材料中に存在するかどうかを検出するために、当業者に既知の標準的方法が適用され得る(例えば核酸増幅を用いて)。遺伝的マーカーが乳房炎耐性形質に遺伝的に連鎖するかどうかを判断するために、順列検定(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%の染色体全体水準での閾値は、遺伝的マーカーと、乳房炎耐性及び体細胞数形質との間の連鎖についての有意な証拠であると見なされた。
遺伝的マーカーが遺伝子材料中に存在するかどうかを検出するために、当業者に既知の標準的方法が適用され得る(例えば核酸増幅を用いて)。遺伝的マーカーが乳房炎耐性形質に遺伝的に連鎖するかどうかを判断するために、順列検定(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%の染色体全体水準での閾値は、遺伝的マーカーと、乳房炎耐性及び体細胞数形質との間の連鎖についての有意な証拠であると見なされた。
キット
本発明の別の態様は、ウシ被験体における乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーの有無を検出する際に使用するための診断キットであって、少なくとも1つのオリゴヌクレオチド配列及びそれらの組合せを含み、ヌクレオチド配列が配列番号1〜配列番号192のいずれか及び/又はそれらの任意の組合せから選択される、診断キットに関する。
本発明の別の態様は、ウシ被験体における乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーの有無を検出する際に使用するための診断キットであって、少なくとも1つのオリゴヌクレオチド配列及びそれらの組合せを含み、ヌクレオチド配列が配列番号1〜配列番号192のいずれか及び/又はそれらの任意の組合せから選択される、診断キットに関する。
本発明によるウシ被験体についての乳房炎に対する耐性の遺伝的決定因子(genetic determinants)を確立するためのウシ被験体の遺伝子型同定(genotyping)は、DNA及び/又はRNAの解析に基づき得る。一例としては、本明細書中に記載の標準的なDNA抽出法を使用して提供され得るゲノムDNAである。ゲノムDNAは、例えば多型マーカー(polymorphic marker)領域に対応(相補的)するオリゴヌクレオチドプライマーを使用するポリメラーゼ連鎖反応等の、標準的な技法を使用して単離且つ増幅し得る。増幅反応前にDNAを精製するさらなる工程が含まれ得る。したがって、乳房炎耐性及び体細胞数特性を確立するための診断キットは、別個の包装中に、表xxに示される配列の群から選択される少なくとも1つのオリゴヌクレオチド配列、及びそれらの任意の組合せを含む。[実施例]
動物
動物材料は39個の父系性雄親家系(paternal sire families)を有するグランドドーターデザインから成り、出生児を調べた息子の総数は、4つの乳牛品種、即ちデンマークホルスタイン(DH)及びデンマークレッド(DR)、フィンランドエアシャー(FA)及びスウェーデンレッドアンドホワイト(SRB)で1513頭であった。これらの39個の家系は、5つのDH、9つのDR、11個のFA、及び14個のSRBの祖父家系から成る。祖父1頭当たりの息子の数は16〜161頭の範囲であり、平均家系サイズは38.8であった。
動物材料は39個の父系性雄親家系(paternal sire families)を有するグランドドーターデザインから成り、出生児を調べた息子の総数は、4つの乳牛品種、即ちデンマークホルスタイン(DH)及びデンマークレッド(DR)、フィンランドエアシャー(FA)及びスウェーデンレッドアンドホワイト(SRB)で1513頭であった。これらの39個の家系は、5つのDH、9つのDR、11個のFA、及び14個のSRBの祖父家系から成る。祖父1頭当たりの息子の数は16〜161頭の範囲であり、平均家系サイズは38.8であった。
ゲノム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中に再懸濁する。
以下のプロトコールに従って、精液からゲノム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中に再懸濁する。
増幅手順
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各プライマー)。
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各プライマー)。
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℃に試料を冷却するようにプログラムした。
マーカーの検出に使用されるプライマーのヌクレオチド配列を表50に示す。配列を5'末端から列記する。
マーカー及びマップ
本研究の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)遺伝子型同定された。
本研究の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)遺伝子型同定された。
BTA11の連鎖マップにおける全37個のマーカーの中で、本研究では、30個のマイクロサテライトマーカーを食肉動物研究所のウェブサイト(http://www.marc.usda.gov/genome/cattle/cattle.html)から選択した。
ラディエーションハイブリッド(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反応及びゲル電気泳動を繰り返した。幾つかのハイブリッドの結果が不明瞭なままであった場合、マーカーを廃棄した。
遺伝子及び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反応及びゲル電気泳動を繰り返した。幾つかのハイブリッドの結果が不明瞭なままであった場合、マーカーを廃棄した。
ウシWGRHパネルでこれまでに同定されており、配置が既知(known assignments)のマーカーに対し、RHMAPPER(Soderlund et al., 1998)を用いて2点連鎖分析を行うことによって、ウシ染色体に対してマーカーを割り当てた(assigned)。それから、Williams et al(2002)によって記載されたように、Carthegeneソフトウェア(Schiex., 2002)を用いて、RHマップを構築した。ウシ染色体9上で、我々は、150個のマーカーによるラディエーションハイブリッドマップの情報を有する。
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個のマーカーが包含される。
本発明によるBTA11のQTLマッピングに用いた最終の連鎖マップには、表52で列記されるような37個のマーカーが包含される。
以下の表は、関連するQTLに用いたマーカーを示す。マーカーに対する任意のさらなる情報は、「http://www.marc.usda.gov/」、「http://www.ensembl.org/Bos_taurus/index. html」及び「http://www.ncbi.nih.gov/」に見出すことができる。
表現型データ
雄ウシの娘を、乳房炎耐性及び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。
雄ウシの娘を、乳房炎耐性及び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。
BTA9
統計分析
以下で示されるような多くの統計法が、乳房炎に関連する又はそれと連鎖する遺伝的マーカー、ひいては乳房炎耐性の決定に用いられた。
統計分析
以下で示されるような多くの統計法が、乳房炎に関連する又はそれと連鎖する遺伝的マーカー、ひいては乳房炎耐性の決定に用いられた。
QTL解析
QTLを同定するためには連鎖解析(LA)を用いるが、それは、偶然から予想されるより高頻度で系統内の疾患又は形質価に関連する染色体領域に、遺伝的マーカーを家系内でタイピングすることによって行う。かかる連鎖領域は、偶発的な遺伝的バリアントを包含する可能性が高い。データは一連のモデルを用いて解析した。最初に、すべての形質について多点回帰アプローチを用いた単一形質モデルを家系内で解析した。家系内で有意な効果を有する染色体は、家系間で見出されたQTLを検証し且つQTLを特徴付けるために、分散成分法(variance component method)を用いて解析した。
QTLを同定するためには連鎖解析(LA)を用いるが、それは、偶然から予想されるより高頻度で系統内の疾患又は形質価に関連する染色体領域に、遺伝的マーカーを家系内でタイピングすることによって行う。かかる連鎖領域は、偶発的な遺伝的バリアントを包含する可能性が高い。データは一連のモデルを用いて解析した。最初に、すべての形質について多点回帰アプローチを用いた単一形質モデルを家系内で解析した。家系内で有意な効果を有する染色体は、家系間で見出されたQTLを検証し且つQTLを特徴付けるために、分散成分法(variance component method)を用いて解析した。
回帰分析
マーカーの集団アレル頻度(allele frequencies)を、EMアルゴリズムを使用して推定した。その後、アレル頻度が誤差無く分かると仮定した。雄親におけるフェーズを出生児のマーカー型に基づいて決定した。その後、このフェーズが誤差無く分かると仮定した。各マップ位置での分離確率(segregation probabilities)は、Haldaneのマッピング関数を同時使用して、染色体上のすべてのマーカーからの情報を使用して算出した(Haldane,
1919)。表現型を分離確率に回帰した。有意性閾値(significance thresholds)を、順列検定を使用して算出した(Churchill and Doerge, 1994)。
マーカーの集団アレル頻度(allele frequencies)を、EMアルゴリズムを使用して推定した。その後、アレル頻度が誤差無く分かると仮定した。雄親におけるフェーズを出生児のマーカー型に基づいて決定した。その後、このフェーズが誤差無く分かると仮定した。各マップ位置での分離確率(segregation probabilities)は、Haldaneのマッピング関数を同時使用して、染色体上のすべてのマーカーからの情報を使用して算出した(Haldane,
1919)。表現型を分離確率に回帰した。有意性閾値(significance thresholds)を、順列検定を使用して算出した(Churchill and Doerge, 1994)。
分散成分法(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σ2 h/(2σ2 h+σ2 u)(式中、σ2 h及びσ2 uはそれぞれ、ハプロタイプ効果(haplotypes effect)及び相加的な多遺伝子性効果(polygenic effect)に関連する分散成分に対応する)と推定した。
分散成分(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σ2 h/(2σ2 h+σ2 u)(式中、σ2 h及びσ2 uはそれぞれ、ハプロタイプ効果(haplotypes effect)及び相加的な多遺伝子性効果(polygenic effect)に関連する分散成分に対応する)と推定した。
分散成分解析 単一形質単一QTL解析
各形質を、連鎖解析を使用して別個に解析した。完全なモデル(full model)は以下の通り表現され得る:
各形質を、連鎖解析を使用して別個に解析した。完全なモデル(full model)は以下の通り表現され得る:
多形質多QTL解析
多形質解析を行った。Lund et al., 2003に従って、モデル(1)をモデル(2)に記載されるように多形質多QTLモデルに拡張させることができる。
形質はnq個のQTLを用いた以下の線形混合モデルを使用してモデル化される:
多形質解析を行った。Lund et al., 2003に従って、モデル(1)をモデル(2)に記載されるように多形質多QTLモデルに拡張させることができる。
形質はnq個のQTLを用いた以下の線形混合モデルを使用してモデル化される:
i
は
複合連鎖及び連鎖不平衡解析
複合連鎖及び連鎖不平衡解析(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)
al., 1995)。IBD行列は、それぞれのマーカー区間(marker interval)の中点で計算し、続く分散成分推定手順で用いた。
複合連鎖及び連鎖不平衡解析(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)
al., 1995)。IBD行列は、それぞれのマーカー区間(marker interval)の中点で計算し、続く分散成分推定手順で用いた。
パラメーターの推定
平均情報被制限最大尤度アルゴリズム(average information restricted maximum likelihood algorithm)を用いて、分散成分を推定した(Jensen et al., 1997)。このモデルにおけるランダム効果に関する分散成分に関して、被制限尤度を最大にした。特定位置の格子(a grid of specific positions)に渡る一連の被制限尤度を最大にすることによって、QTL位置の被制限尤度のプロファイルが得られる(Sorensen et al., 2003)。染色体に沿って、それぞれのマーカーブラケットの中点でパラメーターを推定した。
平均情報被制限最大尤度アルゴリズム(average information restricted maximum likelihood algorithm)を用いて、分散成分を推定した(Jensen et al., 1997)。このモデルにおけるランダム効果に関する分散成分に関して、被制限尤度を最大にした。特定位置の格子(a grid of specific positions)に渡る一連の被制限尤度を最大にすることによって、QTL位置の被制限尤度のプロファイルが得られる(Sorensen et al., 2003)。染色体に沿って、それぞれのマーカーブラケットの中点でパラメーターを推定した。
有意水準
ゲノム毎の第一種の過誤(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%の有意水準が有意であると見なした。
ゲノム毎の第一種の過誤(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%の有意水準が有意であると見なした。
BTA9の結果
表53に、BTA9についての回帰分析の結果を示す。図1〜図8は、回帰分析についてのQTLを示すグラフである。家系間QTL解析でQTLを検出するために、分散成分法を使用した。図9〜図16は、分散成分ベースの方法におけるQTLについてのLD、LDLA及びLDプロファイルを示す。図17〜図20はハプロタイプ効果を示す。
表53に、BTA9についての回帰分析の結果を示す。図1〜図8は、回帰分析についてのQTLを示すグラフである。家系間QTL解析でQTLを検出するために、分散成分法を使用した。図9〜図16は、分散成分ベースの方法におけるQTLについてのLD、LDLA及びLDプロファイルを示す。図17〜図20はハプロタイプ効果を示す。
デンマークレッド
家系内回帰分析によって、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個のハプロタイプ及びまた雄親から受け継いだハプロタイプについてのハプロタイプ効果が推定された。高及び低乳房炎耐性に関連したハプロタイプが同定された。
家系内回帰分析によって、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個のハプロタイプ及びまた雄親から受け継いだハプロタイプについてのハプロタイプ効果が推定された。高及び低乳房炎耐性に関連したハプロタイプが同定された。
フィンランドエアシャー
家系内回帰分析をフィンランドエアシャー家系で行った場合に、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%であった。高及び低乳房炎耐性に関連したハプロタイプが同定された。
家系内回帰分析をフィンランドエアシャー家系で行った場合に、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%であった。高及び低乳房炎耐性に関連したハプロタイプが同定された。
スウェーデンレッドアンドホワイト
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個のクラスターが存在していた。
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個のクラスターが存在していた。
デンマークホルスタイン
デンマークホルスタインウシにおいても、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ピークは全く観察されなかった。
デンマークホルスタインウシにおいても、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ピークは全く観察されなかった。
品種間解析
品種内(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ピークは消失した。
品種内(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ピークは消失した。
多形質解析
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のマーカー区間で観察された。
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のマーカー区間で観察された。
ハプロタイプ解析
上記のQTL細密(fine)マッピングの結果は、4マーカー領域(BM4208-BMS2819-INRA144-INRA084)内に、CMについて分離するQTLを指し示す。したがって、創始者ハプロタイプのクラスター形成及びハプロタイプ効果を、マーカーBMS2819及びINRA144間の中点で調べた。これは品種内で行われ、また3つの品種DR、FA及びSRB間でも行われた。これは、これらの3つの品種がその起源で関係しているためである。高及び低乳房炎耐性に関連したハプロタイプが同定された。
上記のQTL細密(fine)マッピングの結果は、4マーカー領域(BM4208-BMS2819-INRA144-INRA084)内に、CMについて分離するQTLを指し示す。したがって、創始者ハプロタイプのクラスター形成及びハプロタイプ効果を、マーカーBMS2819及びINRA144間の中点で調べた。これは品種内で行われ、また3つの品種DR、FA及びSRB間でも行われた。これは、これらの3つの品種がその起源で関係しているためである。高及び低乳房炎耐性に関連したハプロタイプが同定された。
BTA11
統計分析
下記のような多くの統計法が、乳房炎、ひいては乳房炎耐性に関連した、又はそれと連鎖した遺伝的マーカーの決定に用いられた。
統計分析
下記のような多くの統計法が、乳房炎、ひいては乳房炎耐性に関連した、又はそれと連鎖した遺伝的マーカーの決定に用いられた。
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解析を行った。
偶然予測されるものよりも高頻度で、系統内の疾患又は形質価(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解析を行った。
分散成分法
分散成分(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σ2 h/(2σ2 h+σ2 u)と推定した(式中、σ2 h及びσ2 uはそれぞれ、ハプロタイプ効果及び相加的な多遺伝子性効果に関連する分散成分に対応する)。
分散成分(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σ2 h/(2σ2 h+σ2 u)と推定した(式中、σ2 h及びσ2 uはそれぞれ、ハプロタイプ効果及び相加的な多遺伝子性効果に関連する分散成分に対応する)。
分散成分解析 単一形質単一QTL解析
各形質は、連鎖解析を使用して別個に解析した。完全なモデルは以下の通り表現され得る:
各形質は、連鎖解析を使用して別個に解析した。完全なモデルは以下の通り表現され得る:
多形質多QTL解析
多形質解析を行った。Lund et al., 2003に従って、モデル(1)をモデル(2)に記載されるように多形質多QTLモデルに拡張させることができる。
この形質はnq個のQTLで以下の線形混合モデルを用いてモデル化される:
多形質解析を行った。Lund et al., 2003に従って、モデル(1)をモデル(2)に記載されるように多形質多QTLモデルに拡張させることができる。
この形質はnq個のQTLで以下の線形混合モデルを用いてモデル化される:
回帰分析
マーカーでの集団アレル頻度を、EMアルゴリズムを使用して推定した。その後、アレル頻度は誤差無く分かったと仮定した。雄親におけるフェーズは出生児のマーカー型に基づいて求めた。その後、このフェーズは誤差無く分かったと仮定した。各マップ位置での分離確率は、Haldaneのマッピング関数を同時使用して、染色体上のすべてのマーカーからの情報を使用して算出した(Haldane, 1919)。表現型を分離確率に回帰した。有意性閾値を、順列検定を使用して算出した(Churchil and Doerge, 1994)。
マーカーでの集団アレル頻度を、EMアルゴリズムを使用して推定した。その後、アレル頻度は誤差無く分かったと仮定した。雄親におけるフェーズは出生児のマーカー型に基づいて求めた。その後、このフェーズは誤差無く分かったと仮定した。各マップ位置での分離確率は、Haldaneのマッピング関数を同時使用して、染色体上のすべてのマーカーからの情報を使用して算出した(Haldane, 1919)。表現型を分離確率に回帰した。有意性閾値を、順列検定を使用して算出した(Churchil and Doerge, 1994)。
パラメーターの推定
平均情報被制限最大尤度アルゴリズム(average information restricted maximum likelihood algorithm)を用いて、分散成分を推定した(Jensen et al., 1997)。このモデルにおけるランダム効果に関する分散成分に関して、被制限尤度を最大にした。特定位置の格子(a grid of specific positions)にわたる一連の被制限尤度を最大にすることによって、QTL位置の被制限尤度のプロファイルが得られる(Sorensen et al., 2003)。染色体に沿って、それぞれのマーカーブラケットの中点でパラメーターを推定した。相加的な遺伝的分散の全体のうち、QTLによって説明される割合は、2σ2 h/(2σ2 h+σ2 u)と推定した(式中、σ2 h及びσ2 uはそれぞれ、ハプロタイプ効果及び相加的な多遺伝子性効果に関連する分散成分に対応する)。
平均情報被制限最大尤度アルゴリズム(average information restricted maximum likelihood algorithm)を用いて、分散成分を推定した(Jensen et al., 1997)。このモデルにおけるランダム効果に関する分散成分に関して、被制限尤度を最大にした。特定位置の格子(a grid of specific positions)にわたる一連の被制限尤度を最大にすることによって、QTL位置の被制限尤度のプロファイルが得られる(Sorensen et al., 2003)。染色体に沿って、それぞれのマーカーブラケットの中点でパラメーターを推定した。相加的な遺伝的分散の全体のうち、QTLによって説明される割合は、2σ2 h/(2σ2 h+σ2 u)と推定した(式中、σ2 h及びσ2 uはそれぞれ、ハプロタイプ効果及び相加的な多遺伝子性効果に関連する分散成分に対応する)。
IBD確率の推定
連鎖解析:任意の2つの創始者ハプロタイプ(Hs及びHm)のQTLアレル間のIBD確率を0(即ち創始者ハプロタイプは関連しなかった)と仮定する(Meuwissen et al., 2002)。雄親から父系性又は母系性のQTLアレルを受け継ぐ確率を計算するために、雄親ハプロタイプ及び息子の父系性遺伝ハプロタイプを用いる、そして、回帰アルゴリズムを用いてIBD行列を計算した(Wang et al., 1995)。IBD行列は、染色体に沿った2cM区間毎で計算し、続く分散成分推定手順に用いた。
連鎖解析:任意の2つの創始者ハプロタイプ(Hs及びHm)のQTLアレル間のIBD確率を0(即ち創始者ハプロタイプは関連しなかった)と仮定する(Meuwissen et al., 2002)。雄親から父系性又は母系性のQTLアレルを受け継ぐ確率を計算するために、雄親ハプロタイプ及び息子の父系性遺伝ハプロタイプを用いる、そして、回帰アルゴリズムを用いてIBD行列を計算した(Wang et al., 1995)。IBD行列は、染色体に沿った2cM区間毎で計算し、続く分散成分推定手順に用いた。
複合連鎖及び連鎖不平衡解析
複合連鎖及び連鎖不平衡解析(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行列を、それぞれのマーカー区間の中点で計算し、続く分散成分推定手順で用いた。
複合連鎖及び連鎖不平衡解析(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行列を、それぞれのマーカー区間の中点で計算し、続く分散成分推定手順で用いた。
有意水準
ゲノム毎の第一種の過誤(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%の有意水準が有意であると見なした。
ゲノム毎の第一種の過誤(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%の有意水準が有意であると見なした。
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に示す。
表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に示す。
フィンランドエアシャー
回帰分析を用いた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上で分離されることを結論付けることができる。
回帰分析を用いた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上で分離されることを結論付けることができる。
スウェーデンレッドアンドホワイト
スウェーデンレッドアンドホワイト(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ウシにおける臨床型乳房炎に対する多面発現効果を有していない。
スウェーデンレッドアンドホワイト(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ウシにおける臨床型乳房炎に対する多面発現効果を有していない。
品種間解析
スウェーデンレッドアンドホワイト種はノルディックエアシャーウシ種に密接に関連している(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位置を絞り込むことはできなかった。
スウェーデンレッドアンドホワイト種はノルディックエアシャーウシ種に密接に関連している(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位置を絞り込むことはできなかった。
ハプロタイプ効果の推定
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の推定効果を有していた。
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の推定効果を有していた。
臨床型乳房炎及びSCSの形質についてのBTA11上のQTL細密マッピングによって、臨床型乳房炎に影響を及ぼす1つのQTLがフィンランドエアシャーウシにおいて分離し、SCSに影響を及ぼす1つのQTLがフィンランドエアシャー及びスウェーデンレッドアンドホワイトウシの両方で分離していることが確認された。LDLA解析によって、2.1cMの区間へ臨床型乳房炎が細密にマッピングされた。QTL区間内で連鎖不平衡が不足していたために、BTA11上のSCSに影響を与えるQTLは細密にマッピングすることができなかった。
ハプロタイプ解析
ハプロタイプ解析
上記のQTL細密マッピングの結果は、4マーカー領域(MNB-40-AUP1-BM716-DIK2653)内にCMについて分離するQTLを指し示す。したがって、創始者ハプロタイプのクラスター形成及びハプロタイプ効果を、マーカーAUP1及びBM716間の中点で検査した。また、上記のQTL細密マッピング結果は、4マーカー領域(BM304-INRA177-UMBTL20-RM96 INRA177)内のSCSについて分離するQTLを指し示す。したがって、創始者ハプロタイプのクラスター形成及びハプロタイプ効果を、それぞれ、マーカーINRA177及びUMBTL20間の中点で検査した。
これは、FA内で行われた。高及び低乳房炎耐性に関連するハプロタイプが同定された。表61及び表62を参照されたい。
Claims (80)
- ウシ被験体における乳房炎に対する耐性を判断する方法であって、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を該ウシ被験体からの試料中で検出することを含み、該少なくとも1つの遺伝的マーカーが、ウシ染色体BTA9上で、多型マイクロサテライトマーカーC6orf93及びinra084に挟まれ、且つ、多型マイクロサテライトマーカーC6orf93及びinra084を含む領域、並びに/又はウシ染色体BTA11上で、多型マイクロサテライトマーカーHELMTT43及びBM3501に挟まれ、且つ、多型マイクロサテライトマーカーHELMTT43及びBM3501を含む領域に位置し、該少なくとも1つの遺伝的マーカーの有無が、該ウシ被験体又はそれらの出生児の乳房炎耐性を示す、ウシ被験体における乳房炎に対する耐性を判断する方法。
- ウシ被験体における乳房炎に対する耐性を判断する方法であって、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を該ウシ被験体からの試料中で検出することを含み、該少なくとも1つの遺伝的マーカーが、ウシ染色体BTA9上で、多型マイクロサテライトマーカーC6orf93及びinra084に挟まれ、且つ、多型マイクロサテライトマーカーC6orf93及びinra084を含む領域に位置し、該少なくとも1つの遺伝的マーカーの有無が該ウシ被験体又はそれらの出生児の乳房炎耐性を示す、ウシ被験体における乳房炎に対する耐性を判断する請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが乳房炎に対する耐性についてのウシの形質と連鎖する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーC6orf93及びinra084に挟まれ、且つ、遺伝的マーカーC6orf93及びinra084を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2251及びinra084に挟まれ、且つ、遺伝的マーカーbms2251及びinra084を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbm7234及びinra144に挟まれ、且つ、遺伝的マーカーbm7234及びinra144を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2819及びinra144に挟まれ、且つ、遺伝的マーカーbms2819及びinra144を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2251及びinra144に挟まれ、且つ、遺伝的マーカーbms2251及びinra144を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbms2819及びinra084に挟まれ、且つ、遺伝的マーカーbms2819及びinra084を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbm7234及びbms2819に挟まれ、且つ、遺伝的マーカーbm7234及びbms2819を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーbm7234及びbm4208に挟まれ、且つ、遺伝的マーカーbm7234及びbm4208を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーinra144及びrgs17に挟まれ、且つ、遺伝的マーカーinra144及びrgs17を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーinra144及びinra084に挟まれ、且つ、遺伝的マーカーinra144及びinra084を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の遺伝的マーカーinra084及びrgs17に挟まれ、且つ、遺伝的マーカーinra084及びrgs17を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーbm4208及びinra144に挟まれ、且つ、遺伝的マーカーbm4208及びinra144を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーbm7234及びinra084に挟まれ、且つ、遺伝的マーカーbm7234及びinra084を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーEPM2A及びbm7234に挟まれ、且つ、遺伝的マーカーEPM2A及びbm7234を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーbms2251及びbm7234に挟まれ、且つ、遺伝的マーカーbms2251及びbm7234を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーC6orf93である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDIK4986である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーmm12e6である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーPEX3である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDEAD21である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBMS2251である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーEPM2Aである、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM7234である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM4208である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBMS2819である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーINRA144である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーINRA084である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーrgs17である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが遺伝的マーカーの組合せである、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の69.35cM〜74.5cMの領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の71.3cM〜74.5cMの領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の72.3cM〜74.5cMの領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA9の73.95cM〜74.5cMの領域に位置する、請求項1に記載の方法。
- マイクロサテライトマーカーC6orf93を増幅するプライマー対が配列番号75及び配列番号76である、請求項19に記載の方法。
- マイクロサテライトマーカーPEX3を増幅するプライマー対が配列番号81及び配列番号82である、請求項22に記載の方法。
- マイクロサテライトマーカーEPM2Aを増幅するプライマー対が配列番号87及び配列番号88である、請求項25に記載の方法。
- マイクロサテライトマーカーrgs17を増幅するプライマー対が配列番号99及び配列番号100である、請求項31に記載の方法。
- ウシ被験体における乳房炎に対する耐性を判断する方法であって、少なくとも1つの遺伝的マーカーが、ウシ染色体11上で、多型マイクロサテライトマーカーHELMTT43及びBM3501に挟まれ、且つ、多型マイクロサテライトマーカーHELMTT43及びBM3501を含む領域に位置し、該少なくとも1つの遺伝的マーカーが乳房炎耐性と連鎖し、乳房炎耐性を示す少なくとも1つの形質と連鎖する少なくとも1つの遺伝的マーカーの有無を前記被験体からの遺伝子材料中で検出することを包含し、前記少なくとも1つの遺伝的マーカーの存在が乳房炎耐性を示すこと、及び/又は乳房炎耐性を示す出生児を産むことを示す、ウシ被験体における乳房炎に対する耐性を判断する請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーHELMTT43及びINRA177に挟まれ、且つ、遺伝的マーカーHELMTT43及びINRA177を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーHELMTT43及びMNB-70に挟まれ、且つ、遺伝的マーカーHELMTT43及びMNB-70を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーMNB-40及びMNB-70に隣接する領域、及び遺伝的マーカーMNB-40及びMNB-70を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBP38及びINRA131に挟まれ、且つ、遺伝的マーカーBP38及びINRA131を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBM2818及びINRA177に挟まれ、且つ、遺伝的マーカーBM2818及びINRA177を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBMS1953及びBM2818に挟まれ、且つ、遺伝的マーカーBMS1953及びBM2818を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーHELMTT43及びZAP70に挟まれ、且つ、遺伝的マーカーHELMTT43及びZAP70を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーZAP70及びIL18RAに挟まれ、且つ、遺伝的マーカーZAP70及びIL18RAを含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーINRA131及びBM6445に挟まれ、且つ、遺伝的マーカーINRA131及びBM6445を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBM304及びBM7169に挟まれ、且つ、遺伝的マーカーBM304及びBM7169を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の遺伝的マーカーBM7169及びDIK5170に挟まれ、且つ、遺伝的マーカーBM7169及びDIK5170を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーBM6445及びBMS1048に挟まれ、且つ、遺伝的マーカーBM6445及びBMS1048を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーMB110及びBMS2047に挟まれ、且つ、遺伝的マーカーMB110及びBMS2047を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーIL18RA及びBM2818に挟まれ、且つ、遺伝的マーカーIL18RA及びBM2818を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーBM2818及びBM7169に挟まれ、且つ、遺伝的マーカーBM2818及びBM7169を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーBMS2325及びDIK4637に挟まれ、且つ、遺伝的マーカーBMS2325及びDIK4637を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーIL18RA及びAUP1に挟まれ、且つ、遺伝的マーカーIL18RA及びAUP1を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、遺伝的マーカーDIK4637及びUMBTL103に挟まれ、且つ、遺伝的マーカーDIK4637及びUMBTL103を含む領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーIL18RAである、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーMNB-40である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーUAP1である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDIK4637である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーUMBTL103である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーINRA177である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーUMBTL20である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM716である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーDIK2653である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーがマイクロサテライトマーカーBM2818である、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが遺伝的マーカーの一組合せである、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の2.249cM〜97.223cMの領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の2.249cM〜35.098cMの領域に位置する、請求項1に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが、BTA11の19.440cM〜23.829cMの領域に位置する、請求項1に記載の方法。
- マイクロサテライトマーカーIL18RAを増幅するプライマー対が配列番号125及び配列番号126である、請求項60に記載の方法。
- マイクロサテライトマーカーMNB-40を増幅するプライマー対が配列番号127及び配列番号128である、請求項61に記載の方法。
- マイクロサテライトマーカーDIK4637を増幅するプライマー対が配列番号141及び配列番号142である、請求項63に記載の方法。
- マイクロサテライトマーカーUMBTL103を増幅するプライマー対が配列番号143及び配列番号144である、請求項64に記載の方法。
- 前記少なくとも1つの遺伝的マーカーが乳房炎を引き起こす遺伝子と連鎖する、請求項1に記載の方法。
- 請求項1の方法によって乳房の健康特性を判断することを含む、育種目的でウシ被験体を選抜する方法。
- 乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーのウシ被験体における有無を検出する際に用いる診断キットであって、少なくとも1つのオリゴヌクレオチド配列及びその組合せを包含し、該ヌクレオチド配列が配列番号1〜配列番号192のいずれか、及び/又はそれらの任意の組合せから選択される、乳房炎に対する耐性に関連した少なくとも1つの遺伝的マーカーのウシ被験体における有無を検出する際に用いる診断キット。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DKPA200600164 | 2006-02-06 | ||
DKPA200700147 | 2007-01-30 | ||
PCT/DK2007/000058 WO2007090399A2 (en) | 2006-02-06 | 2007-02-05 | Qtls for mastitis resistance in cattle |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2009525733A JP2009525733A (ja) | 2009-07-16 |
JP2009525733A5 true JP2009525733A5 (ja) | 2012-11-08 |
Family
ID=38110214
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008553616A Pending JP2009525733A (ja) | 2006-02-06 | 2007-02-05 | 乳房炎耐性 |
Country Status (13)
Country | Link |
---|---|
US (1) | US20100003675A1 (ja) |
EP (1) | EP2069531B1 (ja) |
JP (1) | JP2009525733A (ja) |
AR (1) | AR059326A1 (ja) |
AU (1) | AU2007214118B2 (ja) |
CA (1) | CA2677523A1 (ja) |
DK (1) | DK2069531T3 (ja) |
ES (1) | ES2425105T3 (ja) |
NZ (1) | NZ571106A (ja) |
PL (1) | PL2069531T3 (ja) |
PT (1) | PT2069531E (ja) |
SI (1) | SI2069531T1 (ja) |
WO (1) | WO2007090399A2 (ja) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8010403B2 (en) * | 2006-12-29 | 2011-08-30 | American Express Travel Related Services Company, Inc. | System and method for targeting transaction account product holders to receive upgraded transaction account products |
WO2010020252A1 (en) * | 2008-08-19 | 2010-02-25 | Viking Genetics Fmba | Methods for determining a breeding value based on a plurality of genetic markers |
WO2010108498A1 (en) * | 2009-03-26 | 2010-09-30 | Aarhus Universitet | Genetic markers for mastitis resistance |
RU2015111205A (ru) * | 2012-08-28 | 2016-10-27 | Орхус Университет | Генетические маркеры устойчивости к маститу |
CN104611425B (zh) * | 2015-01-15 | 2017-08-08 | 中国农业大学 | 一种与中国荷斯坦奶牛金葡菌乳房炎抗性相关的snp分子标记及其应用 |
-
2007
- 2007-02-05 WO PCT/DK2007/000058 patent/WO2007090399A2/en active Application Filing
- 2007-02-05 CA CA002677523A patent/CA2677523A1/en not_active Abandoned
- 2007-02-05 PL PL07702476T patent/PL2069531T3/pl unknown
- 2007-02-05 SI SI200731258T patent/SI2069531T1/sl unknown
- 2007-02-05 AU AU2007214118A patent/AU2007214118B2/en not_active Ceased
- 2007-02-05 PT PT77024768T patent/PT2069531E/pt unknown
- 2007-02-05 ES ES07702476T patent/ES2425105T3/es active Active
- 2007-02-05 JP JP2008553616A patent/JP2009525733A/ja active Pending
- 2007-02-05 EP EP07702476.8A patent/EP2069531B1/en not_active Not-in-force
- 2007-02-05 DK DK07702476.8T patent/DK2069531T3/da active
- 2007-02-05 NZ NZ571106A patent/NZ571106A/en not_active IP Right Cessation
- 2007-02-05 US US12/223,689 patent/US20100003675A1/en not_active Abandoned
- 2007-02-06 AR ARP070100476A patent/AR059326A1/es unknown
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Otsu et al. | Cosegregation of porcine malignant hyperthermia and a probable causal mutation in the skeletal muscle ryanodine receptor gene in backcross families | |
Levin et al. | Genetic map of the chicken Z chromosome using random amplified polymorphic DNA (RAPD) markers | |
Cushwa et al. | Applications of the random amplified polymorphic DNA (RAPD) assay for genetic analysis of livestock species | |
US20070186298A1 (en) | Follistatin gene as a genetic marker for reproductive and performance traits in pigs | |
JP2015526099A (ja) | 乳房炎抵抗性の遺伝子マーカ | |
Ferreira et al. | Genetic diversity and population structure of different varieties of Morada Nova hair sheep from Brazil. | |
AU2009275988B2 (en) | A genetic marker test for Brachyspina and fertility in cattle | |
US7785778B2 (en) | Porcine polymorphisms and methods for detecting them | |
JP2009525733A5 (ja) | ||
JP2009525733A (ja) | 乳房炎耐性 | |
US20060166244A1 (en) | DNA markers for increased milk production in cattle | |
US20060172329A1 (en) | DNA markers for cattle growth | |
CN100406572C (zh) | 用于鉴定牛脊椎畸形综合征携带者的遗传试验 | |
US20080118914A1 (en) | Follistatin gene as a genetic marker for first parity litter size in pigs | |
WO2010108498A1 (en) | Genetic markers for mastitis resistance | |
KR101787260B1 (ko) | 종우 판별을 위한 키트 및 종우 판별 방법 | |
AU2004276248A1 (en) | Adrenergic receptor SNP for improved milking characteristics | |
Kumar et al. | Molecular genetic characterization of four White Leghorn strains using microsatellite markers | |
US20070209084A1 (en) | Adrenergic Receptor SNP for Improved Milking Characteristics | |
Mishra et al. | Application of advanced molecular marker technique for improvement of animal: A critical | |
AU2379399A (en) | Identification of animals resistant to nematode parasite infection | |
MARIANA et al. | Incidence of BLAD and DUMPS carriers in Romanian cattle breeds | |
EP2149611A1 (en) | A genetic marker test for brachyspina and fertility in cattle | |
Soattin | The use of molecular markers for analyzing genes and genomes of livestock | |
Kirkpatrick | Refined mapping of twinning rate QTL on bovine chromosome 5 and analysis of 2 IGF1 as a positional candidate gene 3 4 Eui-Soo Kim, Xianwei Shi, Ozden Cobanoglu, Kent Weigel, P. Jeffrey Berger, and 5 |