(57)【要約】“妊娠受容能をもつ”卵母細胞を同定する遺伝学的手段を提供する。これらの手段は、妊娠受容能をもつ卵母細胞に由来する卵丘細胞に特徴的レベルで（アップレギュレートまたはダウンレギュレートされて）発現する１以上の遺伝子の発現レベルを検出することを含む。この特徴的な遺伝子発現レベル、すなわち本明細書中で“妊娠シグナチャー”と呼ぶパターンを用いて、生存可能な妊娠の発生を妨害または阻止する原因状態、たとえば閉経前状態、他のホルモン機能障害、卵巣機能障害、卵巣嚢胞、癌または他の細胞増殖障害、自己免疫障害などの状態を伴う対象を同定することもできる。好ましい態様において、妊娠シグナチャーはＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰのうち１以上を含むであろう。【選択図】図１
2013510575.tif 000002
【特許請求の範囲】
【請求項１】
下記の段階を含む、受精した際に生存可能な妊娠を生じることができる卵母細胞を同定する非侵襲的方法：
（ｉ）雌性ドナーからの妊娠受容能について検査すべき卵母細胞またはそのドナーの他の卵母細胞に随伴する少なくとも１個の卵丘細胞を入手する；
（ｉｉ）その少なくとも１個の卵丘細胞による少なくとも１つの遺伝子の発現をアッセイする；その遺伝子の発現は、その卵丘細胞に随伴する卵母細胞が、受精させて適切な子宮環境内へ移植した際に生存可能な妊娠をもたらす能力と相関し、その際、遺伝子は表４中のもの、および／またはＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰ、またはそれらのオルソログ、スプライスバリアントもしくは対立遺伝子バリアントから選択される；そして
（ｉｉｉ）その少なくとも１つの遺伝子の発現レベルを、妊娠受容能をもつ卵母細胞に随伴する卵丘細胞による特徴的な発現レベルと比較したものに基づいて、その卵母細胞またはその雌性ドナーに由来する他の卵母細胞が、受精させて適切な子宮環境内へ移植した際に生存可能な妊娠をもたらすことが潜在的に可能であるかどうかを同定する。
【請求項２】
その発現を検出すべき少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される、請求項１に記載の方法。
【請求項３】
その方法により、それらの遺伝子のうち２以上の発現を検出する、請求項２に記載の方法。
【請求項４】
その方法により、それらの遺伝子のうち３以上の発現を検出する、請求項２に記載の方法。
【請求項５】
その方法により、それらの遺伝子のうち４以上の発現を検出する、請求項２に記載の方法。
【請求項６】
その方法により、それらの遺伝子のうち５以上の発現を検出する、請求項２に記載の方法。
【請求項７】
その方法により、それらの遺伝子のうち６以上の発現を検出する、請求項２に記載の方法。
【請求項８】
その方法により、それらの遺伝子のうち７以上の発現を検出する、請求項２に記載の方法。
【請求項９】
その方法により、それらの遺伝子のうち８以上の発現を検出する、請求項２に記載の方法。
【請求項１０】
卵母細胞が哺乳動物の卵母細胞である、請求項１に記載の方法。
【請求項１１】
卵母細胞がヒトの卵母細胞である、請求項１０に記載の方法。
【請求項１２】
卵母細胞がヒト以外の霊長類の卵母細胞である、請求項１０に記載の方法。
【請求項１３】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰから選択される少なくとも５つの遺伝子の発現を測定する、請求項１に記載の方法。
【請求項１４】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰから選択される少なくとも６つの遺伝子の発現を測定する、請求項１３に記載の方法。
【請求項１５】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰから選択される少なくとも７つの遺伝子の発現を測定する、請求項１４に記載の方法。
【請求項１６】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰから選択される少なくとも８つの遺伝子の発現を測定する、請求項１５に記載の方法。
【請求項１７】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰから選択される少なくとも９つの遺伝子の発現を測定する、請求項１４に記載の方法。
【請求項１８】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する少なくとも１０の遺伝子の発現を測定する、請求項１４に記載の方法。
【請求項１９】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する少なくとも１５の遺伝子の発現を同定する、請求項１に記載の方法。
【請求項２０】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する少なくとも２０の遺伝子の発現を同定する、請求項１に記載の方法。
【請求項２１】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する、表４中のものから選択される少なくとも２０〜５０の遺伝子の発現をアッセイする、請求項１に記載の方法。
【請求項２２】
その発現が生存可能な妊娠を潜在的にもたらす卵母細胞の能力と相関する、表４中のものから選択される少なくとも５０〜１００の遺伝子の発現をアッセイする、請求項１に記載の方法。
【請求項２３】
遺伝子の発現をアッセイする方法が差示的遺伝子発現をモニターする方法を使用する、請求項１に記載の方法。
【請求項２４】
その方法が、インデクシングディファレンシャルディスプレイ逆転写酵素ポリメラーゼ連鎖反応（ＤＤＲＴ−ＰＣＲ）を指標とすることを含む、請求項１４に記載の方法。
【請求項２５】
卵母細胞が最高で３５歳のヒト女性から入手される、請求項１に記載の方法。
【請求項２６】
卵母細胞が少なくとも２５歳のヒト女性から入手される、請求項１に記載の方法。
【請求項２７】
卵母細胞が少なくとも３０歳のヒト女性から入手される、請求項１６に記載の方法。
【請求項２８】
卵母細胞が少なくとも３５歳のヒト女性から入手される、請求項１８に記載の方法。
【請求項２９】
卵母細胞が少なくとも４０歳のヒト女性から入手される、請求項１９に記載の方法。
【請求項３０】
その少なくとも１つの遺伝子の異常な発現が、閉経、癌、卵巣機能障害、卵巣嚢胞、自己免疫障害およびホルモン機能障害から選択される状態と相関する、請求項１に記載の方法。
【請求項３１】
下記の段階を含む、受精処置の有効性を評価する方法：
（ｉ）ヒト女性を想定受精増強処置により処置する；
（ｉｉ）処置後にそのヒト女性から卵母細胞およびそれに随伴する卵丘細胞を入手し、そしてその卵母細胞に随伴する少なくとも１つの卵丘細胞による、表４中に含まれるものならびにさらにＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰ、またはそれらのオルソログおよび対立遺伝子バリアントもしくはスプライスバリアントから選択される少なくとも１つの遺伝子の発現を測定する；そして
（ｉｉｉ）その卵母細胞に随伴する細胞によるその少なくとも１つの遺伝子の発現レベルを、正常卵母細胞または妊娠卵母細胞に随伴する卵丘細胞によるその遺伝子の特徴的な発現レベルと比較したものに基づいて、その処置が有効であるかどうかを評価する。
【請求項３２】
受精処置がホルモン療法を含む、請求項３１に記載の方法。
【請求項３３】
対象が閉経期にあり、処置がホルモン置換療法を含む、請求項３１に記載の方法。
【請求項３４】
遺伝子発現をリアルタイムポリメラーゼ連鎖反応（ＲＴ−ＰＣＲ）により検出する、請求項１に記載の方法。
【請求項３５】
遺伝子発現をＲＴ−ＰＣＲにより検出する、請求項３１に記載の方法。
【請求項３６】
遺伝子発現をインデクシングディファレンシャルディスプレイ逆転写酵素ポリメラーゼ連鎖反応（ＤＤＲＴ−ＰＣＲ）を指標とすることにより差示検出する、請求項１または３１に記載の方法。
【請求項３７】
遺伝子発現結果が卵丘細胞からのＲＮＡを用いて得られる、請求項３６に記載の方法。
【請求項３８】
対象において受精潜在性を評価する方法であって、妊娠潜在性を評価しようとする卵母細胞またはその対象から採取した他の卵母細胞に随伴する卵丘細胞による、表４中のもの、ならびにＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰ、またはそれらのオルソログ、スプライスバリアントもしくは対立遺伝子バリアントから選択される特異的妊娠シグナチャー遺伝子の発現レベルを検出し、その発現レベルを、生存可能な妊娠をもたらすことができる卵母細胞に随伴する卵丘細胞によるそれらの遺伝子の特徴的な発現レベルと比較し；そして、その卵丘細胞が妊娠受容能をもつ卵母細胞に特徴的なレベルで１以上の妊娠シグナチャー遺伝子を発現しているかどうかに基づいて、その対象が潜在的に“妊娠受容能をもつ”かどうかを判定することを含む方法。
【請求項３９】
少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される少なくとも２つの遺伝子を含む、請求項３８に記載の方法。
【請求項４０】
少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される少なくとも３つの遺伝子を含む、請求項３８に記載の方法。
【請求項４１】
少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される少なくとも４つの遺伝子を含む、請求項３８に記載の方法。
【請求項４２】
少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される少なくとも５つの遺伝子を含む、請求項３８に記載の方法。
【請求項４３】
少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される少なくとも６つの遺伝子を含む、請求項３８に記載の方法。
【請求項４４】
少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される少なくとも７つの遺伝子を含む、請求項３８に記載の方法。
【請求項４５】
少なくとも１つの遺伝子が、ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される８、９、１０、１１、１２、１３、または１４すべての遺伝子を含む、請求項３８に記載の方法。
【請求項４６】
その方法が、インデクシングディファレンシャルディスプレイ逆転写酵素ポリメラーゼ連鎖反応（ＤＤＲＴ−ＰＣＲ）を指標とすることにより差示的遺伝子発現をモニターする、請求項４５に記載の方法。
【請求項４７】
“妊娠シグナチャー”ポリペプチドに特異的に結合する抗体を用いて遺伝子発現を検出する、請求項１に記載の方法。
【請求項４８】
ドナー卵丘細胞における、表４からの遺伝子のセット、および／またはＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰの発現を用いて、異数性から選択される妊娠を排除する卵母細胞異常性を同定する、請求項１に記載の方法。
【請求項４９】
受容能をもつ卵母細胞または受容能をもつ胚を選択するための方法であって、卵母細胞および胚を囲む卵丘細胞において、表４中の遺伝子、および／またはＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰから選択される１以上の遺伝子の発現レベルを測定する段階を含む方法。
【請求項５０】
受容能をもつ卵母細胞を選択するための方法であって、表４中の遺伝子のうち少なくとも５つまたはＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる遺伝子のうち少なくとも５つの発現レベルを測定する段階を含む方法。
【請求項５１】
受容能をもつ胚を選択するための方法であって、胚を囲む卵丘細胞において少なくとも１０の遺伝子の発現レベルを測定する段階を含み、その際、それらの遺伝子が表４中のもの、またはＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群の遺伝子から選択される方法。
【請求項５２】
さらに、試料中の遺伝子の発現レベルを対照と比較することからなる段階を含み、その際、試料と対照の間の遺伝子の発現レベルの差の検出はその卵母細胞または胚が受容能をもつかどうかの指標となる、請求項４９〜５１のいずれか１項に記載の方法。
【請求項５３】
さらに、卵丘細胞に随伴する卵母細胞を受精させることを含む、前記請求項のいずれか１項に記載の方法。
【請求項５４】
さらに、卵母細胞を雌性ホストの子宮に移植することを含む、請求項５３に記載の方法。
【請求項５５】
雌性ホストは卵母細胞と異なる種のものである、請求項５４に記載の方法。
【請求項５６】
雌性ホストは卵母細胞と同じ種のものである、請求項５５に記載の方法。
【請求項５７】
雌性ホストはそれから卵母細胞を得た個体である、請求項５６に記載の方法。
【請求項５８】
雌性ホストはそれから卵母細胞を得た個体と異なる個体である、請求項５６に記載の方法。
【請求項５９】
妊娠受容能について検査すべき卵母細胞に随伴する少なくとも１つの卵丘細胞の遺伝子発現と、妊娠受容能をもつ卵母細胞に随伴する卵丘細胞による特徴的な発現レベルとの比較を、重み付き投票法、ベイズ複合共変量法、対角線形識別法、最近接重心法、ｋ−最近接法、シュランケン重心法、支援ベクターマシーン、複合共変量法、およびそのいずれかの組合わせからなる群から選択される方法を用いて実施する、前記請求項のいずれか１項に記載の方法。
【請求項６０】
妊娠受容能について検査すべき卵母細胞に随伴する少なくとも１つの卵丘細胞の遺伝子発現と、妊娠受容能をもつ卵母細胞に随伴する卵丘細胞による特徴的な発現レベルとの比較を、重み付き投票法を用いて実施する、請求項５９に記載の方法。
【請求項６１】
妊娠受容能について検査すべき卵母細胞に随伴する少なくとも１つの卵丘細胞の、遺伝子ＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰ、またはそれらのオルソログ、スプライスバリアントもしくは対立遺伝子バリアントの遺伝子発現を、妊娠受容能をもつ卵母細胞に随伴する卵丘細胞による特徴的な発現レベルと比較することを含む、請求項５９または６０に記載の方法。
【請求項６２】
さらに、雌性ドナーに由来する卵母細胞を受精させて適切な子宮環境内へ移植した際に生存可能な妊娠を潜在的にもたらすことができるかどうかを示すインディケーターを作成することを含む、前記請求項のいずれか１項に記載の方法。
【請求項６３】
インディケーターをレポートとして提供する、請求項６２に記載の方法。
【請求項６４】
インディケーターを電子ディスプレイ上に表示する、請求項６２に記載の方法。
【請求項６５】
インディケーターを電子コミュニケーションとして提供する、請求項６２に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
関連出願の引照
[0001]  本出願は、US provisional application Serial No. 61/388,296、2010年9月30日出願; US provisional application Serial No. 61/387,313および61/387,286、両方とも2010年月9日28出願; US provisional application Serial No. 61/360,556、2010年7月1日出願およびUS provisional application Serial No. 61/259,783、2009年11月10日出願に基づく優先権を主張する。本出願はUS Serial No 11/584,580、2006年10月23日出願にも関係し、これはUS Serial No. 11/437,797、2006年月5日22出願の一部継続出願であり、これはUS Serial No. 11/091,883、2005年3月29日出願の一部継続出願であり、これはprovisional application No. 60/556,875、2004年3月29日出願に基づく優先権を主張する。これらの出願すべてを全体として本明細書に援用する。
【０００２】
発明の分野
[0002]  本発明は、卵丘細胞（cumulus cell）におけるその発現が、その卵丘細胞に随伴する卵母細胞または同じドナーから得られる卵母細胞が妊娠受容能をもつ（pregnancy competent）かどうか、すなわち体外受精に際して生存可能な妊娠をもたらすことができるかどうかと相関する、２２７のヒト遺伝子群、および好ましい１４遺伝子のセットを同定する。さらに本発明は、卵丘細胞における発現が卵母細胞の受容能と相関するこれら２２７の遺伝子の同定および好ましい１４遺伝子のセットの同定をもたらした遺伝子発現検出法および統計分析法を提供する。
【０００３】
[0003]  この知見に基づいて、本発明は、卵母細胞に随伴するかまたは同じドナーに由来する卵丘細胞によるこれら２２７の遺伝子のうち１以上またはこれら１４の遺伝子のうち１以上の発現レベルを検出することによってＩＶＦ法に用いるのに潜在的に適したヒト卵母細胞を同定するための方法および検査キットを提供する。さらにこの知見に基づいて、本発明はさらに、受精させて適切な子宮環境へ移植した際に生存可能な妊娠をもたらす可能性がより高いヒト卵母細胞を同定するための検査キットを提供する。卵丘細胞におけるその発現が妊娠潜在性と相関する２２７の遺伝子のセットは後記の表４に含まれる。さらに、好ましい１４遺伝子のセットは表１２中にあり、ＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる。
【０００４】
[0004]  以上に基づいて、本発明はさらに、雌性対象から単離した卵母細胞から得た卵丘細胞における、表４に含まれるこれら２２７の特異的遺伝子または好ましい１４遺伝子のセットのうち１以上の発現を分析することにより、年齢（閉経）、原疾患状態または薬物療法による卵巣機能障害の結果として妊娠機能障害をもつ雌性対象、好ましくはヒト女性を同定するための遺伝学的方法を提供する。
【０００５】
[0005]  本発明はまた、受精処置またはホルモン処置を受けた雌性対象の卵丘細胞によるこれら２２７または１４の特異的遺伝子のうち１以上の発現にその処置が及ぼす影響を評価することにより、想定受精処置またはホルモン処置の有効性を評価する方法を提供する。
【背景技術】
【０００６】
[0006]  発明の背景
[0007]  現在、ある雌性対象が、“妊娠受容能をもつ”卵母細胞、すなわち自然または人為的手段で受精した際に胚を生成でき、それが次いで適切な子宮環境へ移行した際に生存可能な子孫をもたらすことができる卵母細胞を産生するかどうかを同定するための信頼できる商業的に入手される遺伝学的方法または遺伝学以外の方法はない。むしろ、一般的な受精能評価方法は、ホルモンレベル、卵母細胞の数および質の視診、雌性生殖系臓器の外科的または非侵襲的（ＭＲＩ）検査などに基づいて受精能を評価する。しばしば、ある女性が生存可能な妊娠の形成に問題をもつ場合、長期間、たとえば１年以上後には、診断は“説明できない”受精問題となる可能性があり、その女性は単純に試みを続けるか、あるいは他の選択肢、たとえば養子または代理出産を求めるようにアドバイスされる。
【０００７】
[0008]  一部は、妊娠受容能をもつ卵母細胞を同定するための手段がないためであろうが、人工生殖技術（assisted reproductive technology）（ＡＲＴ）の成功率、すなわち体外受精（in vitro fertilization）（ＩＶＦ）の試みによる妊娠率および出産率は依然として低い。主観的な形態学的パラメーターが今でもＩＶＦおよびＩＣＳＩプログラムに使用する健全な胚を選択するための主な基準である。しかし、そのような基準では胚の受容能を実際に予測することはできない。多数の研究が、幾つかの異なる形態学的基準の組合わせによってより厳密に胚が選択されることを示している。胚の選択のための形態学的基準は移植する当日に評価され、原則として初期胚の卵割（受精後２５〜２７時間目）、２日目、３日目または５日目の卵割球の数およびサイズ、４または８細胞期のフラグメンテーション率および多核化の存在に基づく(Fenwick et al., Hum Reprod, 17, 407-12. (2002)。
【０００８】
[0009]  最近の研究により、受精用の卵母細胞の選択は得られる胚すべてをそれらの質に関係なく移植するのと比較してＡＲＴの結果を改善しないことが示された(La Sala et al., Fertil Steril. (2008))。
【０００９】
[0010]  ＩＶＦ成功率を高め、新鮮置換のための胚の数を減らして多胎妊娠率を低下させるために、最高の着床潜在性をもつ生存可能な胚を同定することが求められている。これらすべての理由で、胚選択のための幾つかのバイオマーカーが現在調べられている(Haouzi et al., Gynecol Obstet Fertil, 36, 730-742. (2008); He et al., Nature, 444, 12-3. (2006))。
【００１０】
[0011]  妊娠をもたらす胚は、そうではない胚と比較してそれらの代謝プロフィールが異なるので、幾つかの研究は非侵襲的な胚培養培地の評価により検出できる分子シグナチャーの同定を試みている(Brison et al.,. Hum Reprod, 19, 2319-24. (2004); Gardner et al., Fertil Steril, 76, 1175-80. (2001); Sakkas and Gardner, Curr Opin Obstet Gynecol, 17, 283-8 (2005); Seli et al., Fertil Steril, 88, 1350-7. (2007); Zhu et al. Fertil Steril.(2007)。
【００１１】
[0012]  ゲノミクスも、胚の発生に際しての遺伝学的機能および細胞機能の重要な知識を提供しつつある。McKenzie et al., Hum Reprod, 19, 2869-74. (2004); Feuerstein et al, Hum Reprod, 22, 3069-77は、卵丘細胞中の幾つかの遺伝子、たとえばシクロオキシゲナーゼ２（ＣＯＸ２）の発現が卵母細胞および胚の質の指標になると報告した。さらに、グレムリン１（Gremlin 1）（ＧＲＥＭ１）、ヒアルロン酸シンターゼ２（ＨＡＳ２）、ステロイド形成急性調節タンパク質（steroidogenic acute regulatory protein）（ＳＴＡＲ）、ステアロイル−補酵素Ａデサチュラーゼ１および５（ＳＣＤ１および５）、アンフィレギュリン（amphiregulin）（ＡＲＥＧ）、ならびにペントラキシン３（pentraxin 3）（ＰＴＸ３）も胚の質と正の相関性をもつと報告されている(Zhang et al., Fertil Steril, 83 Suppl 1, 1169-79. (2005))。より最近になって、ヒト卵丘細胞中のグルタチオンペルオキシダーゼ３（ＧＰＸ３）、ケモカイン受容体４（ＣＸＣＲ４）、サイクリンＤ２（ＣＣＮＤ２）およびカテニン  デルタ１（ＣＴＮＮＤ１）は胚の質と逆の相関性をもつことが、胚の発生に際しての初期卵割速度に基づいて示された(van Montfoort et al., (2008) Mol Hum Reprod, 14, 157-68. (2008))。
【００１２】
[0013]  Cillo et al., Reprod. 134: 645-50 (2007)も、特定の卵丘細胞遺伝子、すなわちＨＡＳ２、ＧＲＥＭ１およびＰＴＸ３の発現と、卵母細胞の質および胚の発生との相関性を示唆している。さらに、Assidi et al. Biol. Reprod. 79(2) 209-222 (2008)は、特定の卵丘細胞遺伝子、すなわちＥＧＦＲ、ＣＤ４４、ＨＡＳ２、ＰＴＳＧ２およびＢＴＣの発現と、卵母細胞の質およびそれからの胚の発生との相関性を示唆している。さらに、Bettegowda et al., Biol. Reprod. 79(2): 301-309 (2008)は、特定のプロテイナーゼカテプシン遺伝子の発現と、ウシ卵母細胞の質およびそれからの子孫の発生との相関性を示唆している。
【００１３】
[0014]  さらに、Zhangらに最近付与された特許(2009年8月11日)は、卵丘細胞におけるペントラキシン３およびＢＣＬ−２メンバーを検出して卵母細胞の質を評価することを権利請求している。US20040058975、2004年3月25日公開も、ＥＰ２受容体および／またはシクロオキシゲナーゼＣＯＸ−２の拮抗が卵丘細胞の増殖および卵母細胞の発育を促進すると教示している。
【００１４】
[0015]  また、初期卵割は妊娠を予測するための信頼性のあるバイオマーカーであることが示されている(Lundin et al, Hum Reprod, 16, 2652-7. (2001); Van Montfoort et al., Hum Reprod, 19, 2103-8 (2004); Yang et al, Fertil Steril, 88, 1573-8 (2007))が、妊娠予後と相関する卵丘細胞の遺伝子発現プロフィールについてはほとんど報告されていない(ただし、Assou et al., Mol Hum Reprod. 2008 Dec; 14(12): 711-9. Epub 2008 Nov 21を参照)。
【先行技術文献】
【特許文献】
【００１５】
【特許文献１】Zhang et al. (2009年8月11日)
【特許文献２】US20040058975
【非特許文献】
【００１６】
【非特許文献１】Fenwick et al., Hum Reprod, 17, 407-12. (2002)
【非特許文献２】La Sala et al., Fertil Steril. (2008)
【非特許文献３】Haouzi et al., Gynecol Obstet Fertil, 36, 730-742. (2008)
【非特許文献４】He et al., Nature, 444, 12-3. (2006)
【非特許文献５】Brison et al.,. Hum Reprod, 19, 2319-24. (2004)
【非特許文献６】Gardner et al., Fertil Steril, 76, 1175-80. (2001)
【非特許文献７】Sakkas and Gardner, Curr Opin Obstet Gynecol, 17, 283-8 (2005)
【非特許文献８】Seli et al., Fertil Steril, 88, 1350-7. (2007)
【非特許文献９】Zhu et al. Fertil Steril. (2007)
【非特許文献１０】McKenzie et al., Hum Reprod, 19, 2869-74. (2004)
【非特許文献１１】Feuerstein et al, Hum Reprod, 22, 3069-77
【非特許文献１２】Zhang et al., Fertil Steril, 83 Suppl 1, 1169-79. (2005)
【非特許文献１３】van Montfoort et al., (2008) Mol Hum Reprod, 14, 157-68.(2008)
【非特許文献１４】Cillo et al., Reprod. 134: 645-50 (2007)
【非特許文献１５】Assidi et al. Biol. Reprod. 79(2) 209-222 (2008)
【非特許文献１６】Bettegowda et al., Biol. Reprod. 79(2): 301-309 (2008)
【非特許文献１７】Lundin et al, Hum Reprod, 16, 2652-7. (2001)
【非特許文献１８】Van Montfoort et al., Hum Reprod, 19, 2103-8 (2004)
【非特許文献１９】Yang et al, Fertil Steril, 88, 1573-8 (2007)
【非特許文献２０】Assou et al., Mol Hum Reprod. 2008 Dec; 14(12): 711-9
【非特許文献２１】Epub 2008 Nov 21
【発明の概要】
【発明が解決しようとする課題】
【００１７】
[0016]  したがって、以上にもかかわらず、ＩＶＦ法に用いるのに適した卵母細胞を同定し、女性における遺伝子に基づく受精問題を同定するための、代替となる、より予測性のある方法を提供することがきわめて望ましいであろう。特に、卵丘細胞によるその発現が卵母細胞の妊娠受容能と相関する他の遺伝子およびバイオマーカーの同定、ならびにそれらを用いた検査キットおよびアッセイ法はＩＶＦ法の結果を向上させることができるので、きわめて望ましいであろう。
【００１８】
[0017]  これらの方法および検査キットはさらに、卵母細胞関連の受精問題を伴う女性の同定を提供するであろう；これは、そのような受精問題は妊娠を排除する他の健康問題、たとえば癌、閉経状態、ホルモン機能障害、卵巣嚢胞、または他の原疾患もしくは健康関連問題と相関している可能性があるので望ましい。
【課題を解決するための手段】
【００１９】
[0018]  本発明は、受容能をもつ卵母細胞を選択するための方法であって、表４中の２２７の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のいずれかの発現レベルを測定する段階を含む方法に関する。
【００２０】
[0019]  本発明はまた、受容能をもつ胚を選択するための方法であって、卵母細胞を囲む卵丘細胞における特異的遺伝子の発現レベルを測定する段階を含み、その際、遺伝子は表４中のものまたはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子である方法に関する。
【００２１】
[0020]  本発明はまた、受容能をもつ卵母細胞または受容能をもつ胚を選択するための方法であって、卵母細胞または胚を囲む卵丘細胞において表４中の２２７の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子から選択される１以上の遺伝子の発現レベルを測定する段階を含む方法に関する。
【００２２】
[0021]  これらの遺伝子のうち１以上の異常な発現レベルは、早期の胚分裂停止に起因する受容能のない卵母細胞または胚の予測に役立つ。
[0022]  前記のように、女性ドナーの卵丘細胞によるこれらの遺伝子の発現レベルは、その卵丘細胞に随伴するかまたは同じ対象に由来する卵母細胞が自然または人為的手段により受精した際に“妊娠受容能をもつ”可能性と相関することが見出された。これらの遺伝子および発現レベルは、本出願人らが“妊娠シグナチャー”と呼ぶものを構成する。さらに、妊娠シグナチャーは前記に示した本出願人の先の出願に開示した１以上の遺伝子をさらに含むことができる。
【００２３】
[0023]  “妊娠シグナチャー”を構成する遺伝子または対応するポリペプチドのうち１以上の発現検出レベルに基づいて、ある個体がその個体の卵母細胞をＩＶＦ法に使用するのに潜在的に不適切にする遺伝子関連の受精問題をもつかどうかを判定する新規方法を提供することは、本発明の関連する目的である。妊娠シグナチャーを構成する遺伝子および遺伝子産物は、この場合も好ましくは表４に含まれるものから選択され、および／またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される。
【００２４】
[0024]  本発明の他の目的は、下記の段階を含む、雌性受精処置の有効性を評価する方法を提供することである：
（ｉ）“生存可能な妊娠”をもつことを阻止または阻害する問題を伴うと推定される雌性対象を処置する；
（ｉｉ）その受精処置後にその雌性対象から少なくとも１つの卵母細胞およびそれに随伴する細胞を単離する；
（ｉｉｉ）単離した卵母細胞に随伴する少なくとも１つの卵丘細胞を単離し、“妊娠受容能をもつ”卵母細胞に特徴的な発現レベルで発現する表４中のものから選択される少なくとも１つの遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰから選択される少なくとも１つの遺伝子の発現レベルを検出する；そして
（ｉｖ）処置の結果としてその遺伝子が“妊娠受容能をもつ”卵母細胞に特徴的な発現レベルで発現しているかどうかに基づいて、その受精処置の推定有効性を判定する。
【００２５】
[0025]  本発明の他の特定の目的は、妊娠シグナチャーを構成する１以上の遺伝子の発現を調節することにより不妊症を処置するための新規方法を提供することである。これらの方法は、表４に含まれる遺伝子、またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰのうち１以上、ならびにそれらのスプライスバリアントまたは対立遺伝子バリアントの発現に作動または拮抗する化合物を投与することを含む。
【００２６】
[0026]  本発明の他の目的は、想定受精処置の有効性を評価するための動物モデルを提供することであり、ヒト以外の動物、たとえばヒト以外の霊長類の妊娠受容能をもつ卵母細胞に随伴する卵丘細胞に特徴的なレベルで発現する遺伝子を同定し；そのヒト以外の動物において想定受精処置がそれらの遺伝子の発現レベルに及ぼした作用、すなわちその処置の結果として、妊娠受容能をもつ卵母細胞に随伴する卵丘細胞にみられる遺伝子発現レベル（“妊娠シグナチャー”）をより良好に模倣した遺伝子発現レベルが、表４中の２２７の遺伝子のうち１以上、またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる１４の遺伝子群のうち１以上において生じるかどうかに基づいて、その処置の有効性を評価することを含む。
【図面の簡単な説明】
【００２７】
【図１ａ】[0027]  図１ａ〜ｃは、ＦおよびＮの全試料（６５試料）、トレイニングセット（training set）の試料（３３試料）および検証セット（validation set）の試料（３２試料）のクラスタリングを、全遺伝子を用いて別個に示す（ａ：全試料，ｂ：トレイニングセット，ｃ：検証セット）。類似度の距離（metric of similarity）としてピアソンの相関（Pearson's correlation）を用いる平均連鎖法（average linkage method）を利用し、行正規化（row normalization）に従って、試料を階層的クラスタリング法（hierarchical clustering）でクラスタリングする；Sneath, P. (1973) Numerical taxonomy; the principles and practice of numerical classification. W. H. Freeman, San Francisco, CA USA。図１：全遺伝子を用いたＮおよびＦ試料のクラスタリング（ａ：全試料，ｂ：トレイニングセット，ｃ：検証セット）。
【図１ｂ】[0027]  図１ａ〜ｃは、ＦおよびＮの全試料（６５試料）、トレイニングセット（training set）の試料（３３試料）および検証セット（validation set）の試料（３２試料）のクラスタリングを、全遺伝子を用いて別個に示す（ａ：全試料，ｂ：トレイニングセット，ｃ：検証セット）。類似度の距離（metric of similarity）としてピアソンの相関（Pearson's correlation）を用いる平均連鎖法（average linkage method）を利用し、行正規化（row normalization）に従って、試料を階層的クラスタリング法（hierarchical clustering）でクラスタリングする；Sneath, P. (1973) Numerical taxonomy; the principles and practice of numerical classification. W. H. Freeman, San Francisco, CA USA。図１：全遺伝子を用いたＮおよびＦ試料のクラスタリング（ａ：全試料，ｂ：トレイニングセット，ｃ：検証セット）。
【図１ｃ】[0027]  図１ａ〜ｃは、ＦおよびＮの全試料（６５試料）、トレイニングセット（training set）の試料（３３試料）および検証セット（validation set）の試料（３２試料）のクラスタリングを、全遺伝子を用いて別個に示す（ａ：全試料，ｂ：トレイニングセット，ｃ：検証セット）。類似度の距離（metric of similarity）としてピアソンの相関（Pearson's correlation）を用いる平均連鎖法（average linkage method）を利用し、行正規化（row normalization）に従って、試料を階層的クラスタリング法（hierarchical clustering）でクラスタリングする；Sneath, P. (1973) Numerical taxonomy; the principles and practice of numerical classification. W. H. Freeman, San Francisco, CA USA。図１：全遺伝子を用いたＮおよびＦ試料のクラスタリング（ａ：全試料，ｂ：トレイニングセット，ｃ：検証セット）。
【図２】[0028]  図２ａ〜ｂは、１１８０の記述遺伝子（descriptive gene）を用いたＮおよびＦ試料のクラスタリングを示す（ａ：トレイニングセット，ｂ：検証セット）。その結果は、トレイニングセット（Ｆ対Ｎ）で同定されたｔ−検定に基づく差次発現を示す（多重仮説検定のためにボンフェロニ補正（Bonferroni correction）してｐ＜０．０５）。得られた“記述遺伝子”と呼ばれる１１８０の遺伝子を用いてトレイニングセットおよび検証セットを個別にクラスタリングした。
【図３】[0029]  図３は、２２７の予測遺伝子（predictor gene）を用いたＮおよびＦ試料のクラスタリングを示す（ａ：トレイニングセット，ｂ：検証セット）。このデータはトレイニングセットにおいて不適正予測された試料のみがクラスタリングにおいても誤置される（misplaced）ことを明らかにするが、検証セットクラスタリングにおけるＦ試料とＮ試料の混合挙動は重み付き投票法（weighted voting approach）によりなされた貢献を強調する。
【図４】[0030]  図４は、Ａ）表４中の予測遺伝子セット（ＰＧ’ｓ）に有意性を指定する方法、Ｂ）ＰＧ’ｓを改良する方法、およびＣ）最終的な予測遺伝子セットをさらに分析する方法を模式的に示す。
【図５】[0031]  図５は、各遺伝子について、その遺伝子が４０の数値をもつ試料の数を示す。結果は“古い”３５試料について計算されている。図５には、各遺伝子について４０の数値をもつ試料の数を、本発明の遺伝子（“Ｈａｓａｎ遺伝子”と表示した１９６の遺伝子）およびＴＬＤＡ上の全３７９遺伝子（“全遺伝子”と表示）について別個にプロットしたものを示す。
【図６】[0032]  図６は、４０の数値をもつ遺伝子の数を各試料について示す。結果は“古い”３５試料について計算されている：各遺伝子について、その遺伝子が４０の数値をもつ試料の数を示す。結果は“新しい”１４試料について計算されている。
【図７】[0033]  図７は、各遺伝子について、その遺伝子が４０の数値をもつ試料の数を示す。結果は“新しい”１４試料について計算されている。
【図８】[0034]  図８は、４０の数値をもつ遺伝子の数を各試料について示す。結果は“新しい”１４試料について計算されている。
【図９】[0035]  図９は、下記の因子に基づく遺伝子の分布を示す：その遺伝子が属するグループ（ＰまたはＡ）；ＴＬＤＡとマイクロアレイにおけるその遺伝子のアップ／ダウンレギュレーションの一致（方向が同じであれば１０、そうでなければ−１０）；その遺伝子が４０の数値をもつ試料の数。この分析は、種々の個数の異常値を除外して、スケーリングした（scaled）数値およびスケーリングしていない（unscaled）数値について個別に実施される。
【発明を実施するための形態】
【００２８】
[0036]  発明の詳細な記述
[0037]  本発明をより詳細に述べる前に、以下の定義を示す。本出願においてその他のすべての語および句は、本発明の分野の当業者が解釈するようにそれらの通常の意味により解釈すべきである。
【００２９】
[0038]  “妊娠受容能をもつ卵母細胞”は、自然または人為的手段により受精した場合にそれが適切な子宮環境内で構成されると生存可能な妊娠をもたらすことができる雌性生殖細胞または卵子を表わす。
【００３０】
[0039]  用語“受容能をもつ胚”は、同様に妊娠に導く高い着床率をもつ胚を表わす。用語“高い着床率”は、子宮内へ移植した際に胚が子宮環境内に着床して生存可能な胎児を生じ、その妊娠を終結させる措置または事象がなければそれが発育して生存可能な子孫になる潜在性を意味する。
【００３１】
[0040]  “生存可能な妊娠”は、受精した卵母細胞が適切な子宮環境内に収容された際に発育すること、およびそれが発育して生存可能な胎児になり、その妊娠を終結させる措置または事象がなければそれが発育して生存可能な子孫になることを表わす。
【００３２】
[0041]  “卵丘細胞”は、卵母細胞を囲む細胞塊に含まれる細胞を表わす。これは“卵母細胞に随伴する細胞”の例である。これらの細胞は、受精した際に“妊娠受容能をもつ”卵母細胞をもたらすのに必要なそれの栄養要求その他の要求の一部を卵母細胞に提供することに関与すると考えられている(Buccione, R., Schroeder, A.C., and Eppig, J.J. (1990). Interactions between somatic cells and germ cells throughout mammalian oogenesis. Biol Reprod 43, 543-547.)。
【００３３】
[0042]  “差示的遺伝子発現”は、目的組織、本明細書中では好ましくは卵母細胞に随伴する細胞、たとえば卵丘細胞内で、その発現が変動する遺伝子を表わす。
[0043]  “リアルタイムＲＴ−ＰＣＲ”は、試料中の特定のＲＮＡ転写体の増幅と定量が同時にできる方法またはそれに用いる装置を表わす。
【００３４】
[0044]  “マイクロアレイ分析”は、特定の試料、たとえば組織試料または細胞試料中の、特異的遺伝子の発現レベルの定量を表わす。
[0045]  本明細書中で“妊娠シグナチャー”は、表４中の特異的遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子、ならびにそれらのオルソログ、スプライスバリアントまたは対立遺伝子バリアントから選択される１以上の遺伝子またはそれらによりコードされるポリペプチドの正常な発現レベルを表わす；その際これらの遺伝子またはポリペプチドは、正常な卵丘細胞において、これらの１以上の遺伝子またはポリペプチドを特徴的レベルで発現する卵丘細胞に随伴する卵母細胞は生存可能な妊娠を生じる可能性がより高いという可能性と相関するレベルで発現する。
【００３５】
[0046]  本明細書中で、特定の検出した発現核酸配列またはポリペプチドに関して“卵丘遺伝子の特徴的発現レベル”は、特定の遺伝子またはポリペプチドが、正常な卵丘細胞に随伴する卵丘細胞または正常もしくは発生受容能をもつ卵母細胞に随伴するものにみられるレベルと実質的に類似するレベルで発現することを意味する。
【００３６】
[0047]  “実質的に類似する”は、個々の遺伝子の発現レベルが、正常な卵丘細胞による遺伝子またはポリペプチドの発現の検出レベルの好ましくは±１〜５倍の範囲内、より好ましくは±１〜３倍の範囲内、よりさらに好ましくは±１〜１．５倍の範囲内、最も好ましくは±１．０〜１．３、１．０〜１．２、１．０〜１．２倍の範囲内にあることを意味する。
【００３７】
[0048]  本発明によれば、卵母細胞は自然サイクル、変異させた自然サイクル、またはｃＩＶＦもしくはＩＣＳＩのために刺激したサイクルから生じたものであってもよい。用語“自然サイクル”は、雌または女性が卵母細胞を産生する自然サイクルを表わす。用語“変異させた自然サイクル”は、組換えＦＳＨまたはｈＭＧに関連するＧｎＲＨアンタゴニストによる緩和な卵巣刺激下で、雌または女性が１または２個の卵母細胞を産生するプロセスを表わす。用語“刺激したサイクル”は、組換えＦＳＨまたはｈＭＧに関連するＧｎＲＨアゴニストまたはアンタゴニストによる緩和な卵巣刺激下で、雌または女性が１以上の卵母細胞を産生するプロセスを表わす。
【００３８】
[0049]  “適切な妊娠シグナチャーをもつかまたは妊娠受容能をもつと判定された卵母細胞または卵丘細胞”は、卵母細胞またはその卵母細胞に随伴する卵丘細胞、あるいは同じ対象から被験卵丘細胞とほぼ同じ時点で（０〜６か月以内）得られた卵母細胞であって、表４中のものまたはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子、またはそれらのオルソログ、スプライスバリアントもしくは対立遺伝子バリアントのうち少なくとも１つの遺伝子またはそれらによりコードされるポリペプチドを正常な卵丘細胞が発現するレベルに特徴的な様式で発現すると判定されたものを表わす。表４に含まれるもの、またはそれらの対立遺伝子バリアントもしくはスプライスバリアントのうち、好ましくは少なくとも２もしくは３つの遺伝子、より好ましくは少なくとも３〜１０の遺伝子、１０〜５０の遺伝子、さらには最高１００またはそれ以上の遺伝子が特徴的な様式で発現する。多数の、たとえば１０以上のオーダーの遺伝子の発現を対象の遺伝学に基づくアッセイで評価する場合、適切な妊娠シグナチャーは、全部または実質的に全部、すなわち少なくとも７０〜８０％の検出遺伝子が正常な卵丘細胞に特徴的な様式で発現することを意味すると理解すべきである。たとえば、１０の遺伝子の発現を検出する場合、それらの遺伝子のうち少なくとも７、８または９の遺伝子が、好ましくは正常な卵丘細胞、すなわち正常な胚および生存可能な妊娠を生じることができる卵母細胞に随伴する卵丘細胞と一致するレベルで発現するであろう。
【００３９】
[0050]  一般に、妊娠シグナチャーに関して、特徴的な発現レベルは少なくとも３〜５、５〜１０、１０〜２０、可能ならば少なくとも５０〜１００の遺伝子についてみられ、それらは“妊娠受容能をもつ”卵母細胞を囲む卵丘細胞に特徴的なレベルで発現する。これは、分析される卵丘細胞内で遺伝子が発現するレベルおよび遺伝子発現の分布を含むものとする。
【００４０】
[0051]  “妊娠シグナチャー遺伝子”は、正常または“妊娠受容能をもつ”卵母細胞に随伴する卵丘細胞が特徴的なレベルで発現する遺伝子を表わす。これらの遺伝子は表４中に含まれ、さらにＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子、ならびにそれらのオルソログ、スプライスバリアントおよび対立遺伝子バリアントを含む。この表中に、遺伝子はそれらの名称および寄託番号で表示される。本発明はさらにこれらの遺伝子の対立遺伝子バリアントおよびスプライスバリアントならびにオルソログの検出を含むことを理解すべきである。
【００４１】
[0052]  “妊娠シグナチャー遺伝子またはポリペプチドの発現を検出するのに適したプローブ”は、転写された遺伝子または対応するポリペプチドの発現を特異的に検出する核酸配列（単数または複数）またはリガンド、たとえば抗体を表わす。好ましい態様において、発現はリアルタイムＰＣＲ検出法の使用により選択される。
【００４２】
[0053]  “ＩＶＦ”は、体外受精を表わす。
[0054]  用語“古典的な体外受精”または“ｃＩＶＦ”は、卵母細胞を体外、すなわちインビトロで精子により受精させるプロセスを表わす。ＩＶＦは、体内受胎が失敗した場合の主要な受精処置である。用語“細胞質内精子注入（intracytoplasmic sperm injection）”または“ＩＣＳＩ”は、１個の精子を卵母細胞に直接注入する体外受精法を表わす。この方法は、雄性不妊因子を克服するために採用されるのが最も一般的であるが、それは卵母細胞に精子が容易に透過できない場合にも、また時には特に精子の供与を伴う体外受精法として採用できる。
【００４３】
[0055]  “透明帯”は、卵母細胞の最も外側の領域を表わす。
[0056]  “差示発現遺伝子を検出するための方法”は、発現した遺伝子転写体またはコードされたポリペプチドを特異的に検出するプローブを用いて差示的遺伝子発現を定量評価するためのいずれか既知の方法を包含する。そのような方法の例には、下記のものが含まれる：インデクシング（indexing）ディファレンシャルディスプレイ逆転写ポリメラーゼ連鎖反応（differential display reverse transcription polymerase chain reaction）（ＤＤＲＴ−ＰＣＲ；Mahadeva et al, 1998, J. Mol. Biol. 284: 1391-1318; WO 94/01582；サブトラクティブｍＲＮＡハイブリダイゼーション（参照：Advanced Mol. Biol.; R.M. Twyman (1999) Bios Scientific Publishers, Oxford, p. 334、核酸アレイまたはマイクロアレイの使用(参照：Nature Genetics, 1999, vol. 21, Suppl. 1061)、および遺伝子発現の連続分析（serial analysis of gene expression）（ＳＡＧＥ）、たとえばValculesev et al, Science (1995) 270: 484-487を参照)、およびリアルタイムＰＣＲ（ＲＴ−ＰＣＲ）。たとえば、卵母細胞における転写遺伝子の差示的レベルは、ノーザンブロット法および／またはＲＴ−ＰＣＲの使用により検出できる。
【００４４】
[0057]  好ましい方法はＣＲＬ増幅プロトコルであり、これは本出願人の先の出願に開示した新規な全ＲＮＡ増幅プロトコルであって鋳型スイッチングＰＣＲとＴ７を基礎とする増幅法を組み合わせたものを表わす。このプロトコルは、数個の細胞または限定された全ＲＮＡしかない試料に好適である。
【００４５】
[0058]  好ましくは“妊娠シグナチャー”遺伝子は、ＤＮＡチップ、たとえばｃＤＮＡ配列を含むフィルターアレイ、またはｃＤＮＡもしくはインサイチュー合成オリゴヌクレオチド配列を含むガラスチップへの、ＲＮＡまたはＤＮＡのハイブリダイゼーションにより検出される。フィルター付きアレイは、一般に高存在量または中等度存在量の遺伝子用として、より良好である。ＤＮＡチップは低存在量の遺伝子を検出できる。例示態様において、ヒトゲノム由来の遺伝子またはそのサブセットを含む試料をAffymetrix GeneChipで検査することができる。
【００４６】
[0059]  あるいは、妊娠シグナチャー遺伝子によりコードされるポリペプチドまたはそれに結合する抗体を含むポリペプチドアレイを作成して検出および診断に使用できる。
[0060]  “ＥＡＳＥ”は遺伝子個体学（gene ontology）プロトコルであり、遺伝子のリストから、あるカテゴリー内にその番号のサブグループ遺伝子がみられる確率に基づいて各遺伝子に指定した機能カテゴリーに基づき、マイクロアレイ上に現われるそのカテゴリーからの遺伝子の頻度を考慮してサブグループを形成する。
【００４７】
[0061]  以上に基づいて、本発明は、ある雌性、好ましくはヒトまたはヒト以外の哺乳動物が“妊娠受容能をもつ”卵母細胞を産生するかどうか、あるいは特定の卵母細胞が妊娠受容能をもつかどうかを検出する新規方法を提供する。この方法は、“妊娠受容能をもつ”卵母細胞に随伴する（それらを囲む）卵丘細胞が特徴的なレベルで発現する表４中の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち１以上の発現レベルを検出することを伴う；すなわち、これらの卵母細胞は自然または人為的手段（ＩＶＦ）により受精して適切な子宮環境内へ移行した際に、その妊娠が何らかの事象または措置、たとえば外科的介入またはホルモン介入により終結されない限り、生存可能な妊娠、すなわち生存可能な胎児に発育して最終的に子孫になる胚をもたらすことができる。
【００４８】
[0062]  本明細書中に記載するように、本発明者らは、卵丘細胞に発現する、胚の潜在性および妊娠予後のバイオマーカーである一組の遺伝子を決定した。本発明者らは、卵母細胞を囲む卵丘細胞の遺伝子発現プロフィールが種々の妊娠予後と相関することを立証し、これにより妊娠に向けて発育する卵母細胞の特異的発現シグナチャーを同定することができた。本発明者らの結果は、卵母細胞を囲む卵丘細胞の分析は卵母細胞選択のための非侵襲的方法であることを示す。
【００４９】
[0063]  表４中の予測遺伝子のセットおよび表１２において同定した１４遺伝子のセットは、既知のヒト遺伝子である。しかし、（卵丘細胞における）これらの遺伝子の発現はこれまで卵母細胞の受容能または胚の発生と関連づけられていなかった。したがって本発明は、受容能をもつ卵母細胞を選択するための方法であって、卵母細胞を囲む卵丘細胞における特異的遺伝子の発現レベルを測定する段階を含み、その際、遺伝子は表４中の２２７の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち少なくとも１つを含む方法に関する。
【００５０】
[0064]  本発明の方法は、さらに試料中の遺伝子の発現レベルを対照と比較することからなる段階を含み、その際、試料と対照の間の遺伝子の発現レベルの差の検出は、その卵母細胞が受容能をもつかどうかの指標となる。対照は、受容能をもつ卵母細胞に随伴する卵丘細胞を含む試料、または未受精（unfertilized）卵母細胞に随伴する卵丘細胞を含む試料で構成されてもよい。
【００５１】
[0065]  本発明の方法は、好ましくはヒト女性に適用できるが、他の哺乳動物（たとえば、霊長類、イヌ、ネコ、ブタ、ウシなど）にも適用できる。
[0066]  本発明の方法は、体外受精処置の有効性を評価するために特に適切である。したがって本発明はまた、雌性対象において管理された卵巣過剰刺激（controlled ovarian hyperstimulation）（ＣＯＳ）プロトコルの有効性を評価するための下記を含む方法に関する：ｉ）雌性対象から少なくとも１個の卵母細胞をそれの卵丘細胞と共に入手し；ｉｉ）本発明の方法によりその卵母細胞が受容能をもつ卵母細胞であるかどうかを判定する。
【００５２】
[0067]  次いで、その方法の後に、受容能をもつ卵母細胞を発生専門医が選択し、そしてたとえば古典的な体外受精（ｃＩＶＦ）プロトコルを用いて、または細胞質内精子注入（ＩＣＳＩ）プロトコル下において、体外でそれらを受精させる。
【００５３】
[0068]  本発明の他の目的は、管理された卵巣過剰刺激（ＣＯＳ）プロトコルの有効性をモニターするための下記を含む方法に関する：ｉ）自然、変異または刺激サイクル下で女性から少なくとも１個の卵母細胞をそれの卵丘細胞と共に単離し；ｉｉ）本発明の方法によりその卵母細胞が受容能をもつ卵母細胞であるかどうかを判定し；ｉｉｉ）そしてＣＯＳ処置の有効性を、それにより受容能をもつ卵母細胞が生成するかどうかに基づいてモニターする。
【００５４】
[0069]  ＣＯＳ処置は、組換えＦＳＨまたはｈＭＧに関連するＧｎＲＨアゴニストまたはアンタゴニストからなる群から選択される少なくとも１種類の有効成分に基づくものであってもよい。
【００５５】
[0070]  本発明はまた、受容能をもつ胚を選択するための方法であって、表４中の２２７の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち少なくとも１つの発現レベルを測定する段階を含む方法に関する。
【００５６】
[0071]  本発明の方法はさらに、試料中の遺伝子の発現レベルを対照と比較することからなる段階を含み、その際、試料と対照の間の遺伝子の発現レベルの差の検出は、その胚が受容能をもつかどうかの指標となる。対照は、生存可能な胎児を生じる胚に随伴する卵丘細胞を含む試料、または生存可能な胎児を生じない胚に随伴する卵丘細胞を含む試料で構成されてもよい。
【００５７】
[0072]  本発明の方法は胚の形態学的考察からの独立をもたらすことが注目される。２つの胚が同じ形態学的観点をもつけれども本発明の方法によれば妊娠に至る着床率が異なるという可能性がある。
【００５８】
[0073]  本発明の方法は、好ましくはヒト女性に適用できるが、他の哺乳動物（たとえば、霊長類、イヌ、ネコ、ブタ、ウシなど）にも適用できる。
[0074]  本発明はまた、ある胚が受容能をもつ卵母細胞であるかどうかを判定するための方法であって、胚を囲む卵丘細胞中の４５の遺伝子の発現レベルを測定することを含み、その際、それらの遺伝子は表４中の２２７の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち少なくとも１つを含む方法に関する。
【００５９】
[0075]  本発明はまた、ある胚が受容能をもつ胚であるかどうかを判定するための下記を含む方法に関する：ｉ）卵母細胞をそれの卵丘細胞と共に入手し；ｉｉ）その卵母細胞を体外受精させ；そしてｉｉｉ）本発明の方法により段階ｉ）の卵母細胞が受容能をもつ卵母細胞であるかどうかを判定することにより、段階ｉｉ）から得られる胚が受容能をもつかどうかを判定する。
【００６０】
[0076]  本発明はまた、受容能をもつ卵母細胞または受容能をもつ胚を選択するための方法であって、卵母細胞または胚を囲む卵丘細胞において、表４中の２２７の遺伝子のうち少なくとも１つまたはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子から選択される１以上の遺伝子の発現レベルを測定する段階を含む方法に関する。選択したこれらの遺伝子のうち１以上の異常な発現は、受容能をもたない卵母細胞または胚、着床できない不能な胚、あるいは早期胚停止が原因で受容能をもたない卵母細胞または胚を予測するのに役立つ。
【００６１】
[0077]  本発明の方法は、雌性の妊娠予後を向上させるために特に適切である。したがって本発明はまた、雌性の妊娠予後を向上させるための下記を含む方法に関する：ｉ）本発明の方法を実施することにより受容能をもつ胚を選択し；ｉｉｉ）段階ｉ）で選択した胚をその雌性の子宮内に移植し、その際、雌性は卵母細胞のドナーであってもよく、そうでなくてもよい。
【００６２】
[0078]  したがって前記の方法は、発生専門医が妊娠予後の潜在性が乏しい胚を子宮に移植するのを避けるのに役立つであろう。前記の方法は、着床および妊娠を導くことができる受容能をもつ胚を選択することにより、多胎妊娠を避けるのにも特に適切である。
【００６３】
[0079]  前記のすべての場合、それらの方法は、卵丘細胞および胚の遺伝子発現プロフィールと妊娠予後の関係を表わす。
[0080]  本発明の遺伝子の発現レベルを測定するための方法：
[0081]  “妊娠シグナチャー”中の遺伝子、すなわち表４中の２２７の遺伝子のうち少なくとも１つ、またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち少なくとも１つの発現レベルの測定は、多様な手法で実施できる。一般に、測定する発現レベルは相対的発現レベルである。
【００６４】
[0082]  より好ましくは、測定は、試料を選択試薬、たとえばプローブ、プライマーまたはリガンドと接触させ、これにより試料中に最初にあった目的とするポリペプチドまたは核酸の存在を検出し、あるいは量を測定することを含む。接触は、いずれか適切な器具、たとえばプレート、マイクロタイターディッシュ、試験管、ウェル、ガラス、カラムなどにおいて実施できる。特定の態様において、接触は試薬でコートした支持体、たとえば核酸アレイまたは特異的リガンドのアレイ上で実施される。支持体は固体または半固体支持体、たとえばガラス、プラスチック、ナイロン、紙、金属、ポリマーなどで構成される適切な支持体のいずれであってもよい。支持体は多様な形態およびサイズのもの、たとえばスライド、膜、ビーズ、カラム、ゲルなどであってもよい。接触は、試薬と試料の核酸またはポリペプチドとの間で形成される検出可能な複合体、たとえば核酸ハイブリッドまたは抗体−抗原複合体に適切な、いずれかの条件下で行なうことができる。
【００６５】
[0083]  好ましい態様において、発現レベルはｍＲＮＡの量を測定することにより測定できる。
[0084]  ｍＲＮＡの量を測定するための方法は当技術分野で周知である。たとえば、試料（たとえば、患者から調製した細胞または組織）中に含有される核酸をまず標準法に従って、たとえば細胞溶解酵素または化学物質溶液を用いて抽出し、あるいは製造業者の指示に従って核酸結合樹脂により抽出する。抽出したｍＲＮＡを次いでハイブリダイゼーション（たとえば、ノーザンブロット分析）および／または増幅（たとえば、ＲＴ−ＰＣＲ）により検出する。好ましくは、定量または半定量ＲＴ−ＰＣＲが好ましい。リアルタイム定量または半定量ＲＴ−ＰＣＲが特に有利である。他の増幅方法には、リガーゼ連鎖反応（ligase chain reaction）（ＬＣＲ）、転写仲介増幅（transcription-mediated amplification）（ＴＭＡ）、鎖置換増幅（strand displacement amplification）（ＳＤＡ）および核酸配列に基づく増幅（nucleic acid sequence based amplification）（ＮＡＳＢＡ）が含まれる。
【００６６】
[0085]  少なくとも１０のヌクレオチドをもち、かつ本発明の目的ｍＲＮＡに対して配列の相補性または相同性を示す核酸は、ハイブリダイゼーションプローブまたは増幅プライマーとして有用である。そのような核酸は同一である必要はないが、一般に匹敵するサイズの相同領域に対して少なくとも約８０％同一、好ましくは少なくとも８５％同一、よりさらに好ましくは９５〜９５％同一であると理解される。特定の態様において、核酸をハイブリダイゼーションの検出のための適宜な手段、たとえば検出可能な標識と組み合わせて用いるのが有利である。蛍光性、放射性、酵素性、または他のリガンド（たとえば、アビジン／ビオチン）を含めた広範な適宜な指示薬が当技術分野で知られている。
【００６７】
[0086]  プローブは一般に、１０〜１０００、たとえば１０〜８００、より好ましくは１５〜７００、一般に２０〜５００ヌクレオチドの長さの一本鎖核酸を含む。プライマーは一般に、より短い１０〜２５ヌクレオチドの長さの一本鎖核酸であり、増幅すべき目的核酸と完全またはほぼ完全に一致するように設計される。プローブおよびプライマーは、それらがハイブリダイズする核酸に“特異的”である；すなわち、それらは好ましくは高緊縮ハイブリダイゼーション条件（最高融解温度Ｔｍに相当する；たとえば５０％ホルムアミド、５×または６×ＳＣＣ。ＳＣＣは０．１５Ｍ  ＮａＣｌ、０．０１５Ｍクエン酸Ｎａである）下でハイブリダイズする。前記のハイブリダイゼーションおよび検出法に用いる核酸プライマーまたはプローブをキットとして組み込むことができる。そのようなキットは、コンセンサスプライマーおよび分子プローブを含む。好ましいキットは、増幅が起きたかどうかを判定するために必要な構成要素も含む。キットは、たとえばＰＣＲ緩衝液および酵素；陽性対照配列、反応制御プライマー；ならびに特異的配列を増幅および検出するための指示も含むことができる。
【００６８】
[0087]  特定の態様において、本発明の方法は卵丘細胞から抽出した全ＲＮＡを用意し、それらのＲＮＡを増幅させ、そして特異的プローブにハイブリダイズさせる段階を含む；より具体的には定量または半定量ＲＴ−ＰＣＲによる。
【００６９】
[0088]  他の好ましい態様において、発現レベルをＤＮＡチップ分析により測定する。そのようなＤＮＡチップまたは核酸マイクロアレイは支持体に化学結合した種々の核酸プローブから構成され、支持体はマイクロチップ、ガラススライドまたはマイクロスフェアサイズのビーズであってもよい。マイクロチップは、ポリマー、プラスチック、樹脂、多糖類、シリカもしくはシリカ系材料、カーボン、金属、無機ガラス、またはニトロセルロースで作成することができる。プローブは、核酸、たとえばｃＤＮＡまたはオリゴヌクレオチドを含むことができ、それらは約１０〜約６０塩基対であってもよい。発現レベルを測定するために、被験対象からの試料を場合によりまず逆転写したものを標識し、ハイブリダイゼーション条件下でマイクロアレイと接触させて、マイクロアレイ表面に付着したプローブ配列に相補的な標的核酸との間に複合体を形成させる。標識付きのハイブリダイズした複合体を次いで検出し、定量または半定量することができる。標識化は種々の方法で、たとえば放射性標識法または蛍光標識法を用いて達成できる。マイクロアレイハイブリダイゼーション法の多数の変法を当業者は利用できる(たとえば、Hoheisel, Nature Reviews, Genetics, 2006, 7: 200-210による総説を参照)。
【００７０】
[0089]  これに関して、本発明はさらに、表４中の２２７の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち少なくとも１つに対して特異的な核酸を保有する固体支持体を含むＤＮＡチップを提供する。
【００７１】
[0090]  これらの遺伝子の発現レベルを測定するための他の方法には、それらの遺伝子がコードするタンパク質の量の測定が含まれる。
[0091]  そのような方法は、試料中に存在するマーカータンパク質と選択的に相互作用しうる結合パートナーと試料を接触させることを含む。結合パートナーは一般に抗体であり、それはポリクローナルまたはモノクローナルであってもよく、好ましくはモノクローナルである。
【００７２】
[0092]  タンパク質の存在は、標準的な電気泳動法および免疫診断法（イムノアッセイ法、たとえば競合、直接反応、またはサンドイッチタイプのアッセイ法を含む）を用いて検出できる。そのようなアッセイ法には下記のものが含まれるが、それらに限定されない：ウェスタンブロット法；凝集試験；酵素標識および仲介型のイムノアッセイ、たとえばＥＬＩＳＡ；ビオチン／アビジンタイプのアッセイ法；ラジオイムノアッセイ；免疫電気泳動；免疫沈降法など。これらの反応は一般に、標識、たとえば蛍光標識、化学発光標識、放射性標識、酵素標識または色素分子の解明、あるいは抗原およびそれと反応した抗体（単数または複数）の間の複合体の形成を検出するための他の方法を含む。
【００７３】
[0093]  上記のアッセイ法は一般に、液相中の結合していないタンパク質を、抗原−抗体複合体が結合した固相支持体から分離することを伴う。本発明の実施に使用できる固体支持体には、ニトロセルロース（たとえば、膜またはマイクロタイターウェルの形のもの）；ポリ塩化ビニル（たとえば、シートまたはマイクロタイターウェル）；ポリスチレンラテックス（たとえば、ビーズまたはマイクロタイタープレート）；ポリフッ化ビニリデン；ジアゾ化紙；ナイロン膜；活性化ビーズ、磁気反応性ビーズなどの支持体が含まれる。より具体的にはＥＬＩＳＡ法を使用でき、その際、マイクロタイタープレートのウェルを被験タンパク質に対する抗体でコートする。次いで、マーカータンパク質を含有するかまたは含有する疑いのある生体試料を、コートされたウェルに添加する。抗体−抗原複合体が形成しうるのに十分な期間のインキュベーション後、プレート（単数または複数）を洗浄して、結合していない部分を除去し、検出可能な状態に標識した二次結合分子を添加する。二次結合分子を捕獲されているいずれかの試料マーカータンパク質と反応させ、プレートを洗浄し、そして二次結合分子の存在を当技術分野で周知の方法を用いて検出する。
【００７４】
[0094]  あるいは、免疫組織化学的（ＩＨＣ）方法が好ましい場合がある。ＩＨＣは特に試料または組織検体中の標的をインサイチュー検出する方法を提供する。試料の全体的な細胞統合性はＩＨＣにおいて維持されるので、目的とする標的の存在と位置の両方を検出できる。一般に試料をホルマリンで固定し、パラフィンに包埋し、切断して染色およびその後の光学顕微鏡による検査のための切片にする。現在のＩＨＣ法は、直接標識化、または二次抗体を基礎とするかもしくはハプテンを基礎とする標識化のいずれかを用いる。既知のＩＨＣシスチムの例には、たとえばEnVision(商標) (DakoCytomation)、Powervision(登録商標) (Immunovision,アリゾナ州スプリングデール)、NBA(商標)キット(Zymed Laboratories Inc.,カリフォルニア州サウスサンフランシスコ)、ヒストファイン(HistoFine) (登録商標) (ニチレイ社,日本東京)が含まれる。
【００７５】
[0095]  特定の態様において、組織切片（たとえば、卵丘細胞を含む試料）を、目的遺伝子によりコードされるタンパク質に対する抗体と共にインキュベートした後、スライドまたは他の支持体上にマウントすることができる。次いで、適切な固体支持体にマウントしたこの試料の顕微鏡検査を行なうことができる。顕微鏡写真を作成するためには、試料を含む切片をガラススライドまたは他の平坦な支持体上にマウントし、選択的染色により目的タンパク質の存在を強調することができる。
【００７６】
[0096]  したがって、ＩＨＣ試料はたとえば下記を含むことができる：（ａ）卵丘細胞を含む標本、（ｂ）それらの細胞の固定および包埋、ならびに（ｃ）それらの細胞試料中の目的タンパク質の検出。ある態様において、ＩＨＣ染色法はたとえば下記の段階を含むことができる：組織の切取りおよびトリミング、固定、脱水、パラフィン浸潤、薄い切片への切断、ガラススライド上へのマウンティング、ベーキング、脱パラフィン、再水和、抗原回復（retrieval）、遮断段階、一次抗体の適用、洗浄、二次抗体（場合によっては適切な検出可能な標識に結合したもの）の適用、洗浄、対比染色、および顕微鏡検査。
【００７７】
[0097]  本発明はまた、前記の方法を実施するためのキットに関するものであり、その際、キットは卵母細胞または胚が受容能をもつかどうかの指標となる表４中の２２７の遺伝子またはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち少なくとも１つの発現レベルを測定するための手段を含む。
【００７８】
[0098]  卵丘細胞における表４中の２２７の遺伝子および１４遺伝子のセットの発現が移植に際しての卵母細胞の受容能および胚の発生と相関することを本発明者らがいかにして決定したかについての以下の記載および実施例により、本発明をさらに説明する。ただし、決してこれらの例および記載が本発明の範囲を限定するとは解釈すべきでない。
【実施例】
【００７９】
[0099]  これを行なった手段の例を以下に詳細に記載する。
[00100]  本発明を実施するための第１観点は、女性および潜在的に他の哺乳動物の妊娠シグナチャーを構成する遺伝子を同定することを伴い、妊娠受容能をもつかまたはもたない女性ドナーから採集した卵丘細胞における遺伝子の発現を同定および比較することにより達成された。これは、ドナー女性の種々のヒト卵母細胞から卵丘細胞を採集し、そして患者に１または２個の想定受精卵を移植することにより行なわれた。これらの患者を次いで移植の結果に基づいて、完全妊娠、部分妊娠および妊娠なしに基づき３グループに分けた。この方法に用いた卵母細胞それぞれについて、５４，０００以上の転写体を含むAffymetrix HG 133 Plus 2アレイを用いてその卵母細胞を囲む少なくとも１つの卵丘細胞の転写プロフィールを決定した。患者が表１に定める除外基準のいずれにも該当しない場合にのみ試験に含めた。
【００８０】
表１：患者の除外基準
【００８１】
【表１】
2013510575.tif 000003
【００８２】
[00101]  より具体的には、卵母細胞が健全な幼児を産生する能力の予測に役立つ遺伝子シグナチャーを見出すために、本発明者らは５４，０００以上の転写体を含むAffymetrix HG 133 Plus 2アレイを用いて、その卵母細胞を囲む少なくとも１つの卵丘細胞のトランスクリプトームをプロファイリングした。個々の卵丘試料からの全ＲＮＡをPicoPure RNA単離キット(Molecular Devices,カリフォルニア州サニーベール)により単離した。少量のＲＮＡをマイクロアレイ分析に必要なレベルにまで確実に一貫して増幅するのを保証する自社開発プロトコルにより試料ＲＮＡを増幅した(Kocabas, et al., Proc Natl Acad Sci U S A, 103, 14027-14032 (2006))。
【００８３】
[00102]  得られた増幅ＲＮＡ（ａＲＮＡ）をAffymetrixアレイにハイブリダイズさせた。どの胚移植も妊娠成功をもたらさなった３６の試料（不成功（No success）に関するＮと表示）およびすべての移植が妊娠成功をもたらした３０の試料（完全妊娠（Full success）に関するＦと表示）を用いた。妊娠成功に影響を及ぼす既知の混同因子はなく、年齢またはＩＶＦサイクル数など関連の臨床パラメーターはＦとＮのグループ間で有意差がなかった。
【００８４】
[00103]  品質管理（ＱＣ）パラメーターを６５試料すべてについて、製造業者(Affymetrix)から無料で入手できるExpression Console（商標）(EC)ソフトウェアにより計算した。スケーリングファクター（全体的チップ強度の２％調整平均(trimmed mean)を同等化するために必要な係数）、存在(present)と呼出されるプローブセットのパーセント、スパイクおよびラベリングコントロールについての３’−５’比、ならびにハウスキーピング遺伝子を含めたすべてのＱＣパラメーターが、すべての試料について許容範囲内であった（製造業者のガイドラインの記載に従った）。妊娠成功に影響を及ぼす既知の混同因子はなく、年齢またはＩＶＦサイクル数など関連の臨床パラメーターはＦとＮのグループ間で有意差がなかった（ｔ−検定  ｐ＞０．０５）（表１を参照）。許容のための追加の基準には、多嚢胞性卵巣症候群（Polycystic Ovarian Syndrome）（ＰＣＯＳ）が存在しないこと、化学療法歴または腹部もしくは骨盤への放射線照射歴がないこと、＞４ｃｍの壁内または粘膜下類線維腫が存在しないこと、また雄性側には、精巣生検歴がないこと、および精子数＞５００万であることが含まれていた。
【００８５】
[00104]  予測モデルの確実性を証明するために、ＦおよびＮ試料をトレイニングセットと検証セットにランダムに分けた。目標は、トレイニングセットに発現した予測遺伝子セットを見出し、次いで予測モデルの開発に用いていない検証セットにおいてこれらの予測遺伝子の性能を検定することであった。この戦略は（すべての試料をシグナチャーの開発に用いるのとは異なり）過剰当てはめ（over-fitting）を防ぎ、予測シグナチャーの頑健性（robustness）の評価を得ることができる(Nevins, J.R. and Potti, A. (2007) Mining gene expression profiles: expression signatures as cancer phenotypes, Nat Rev Genet, 8, 601-609.)。
【００８６】
[00105]  表２に示すように、３３試料（１５Ｆ；１８Ｎ）をトレイニングセットに用い、３２試料（１５Ｆ；１８Ｎ）を検証セットに用いた。
[00106]  トレイニングセットおよび検証セットに用いた試料を表２に示す。
【００８７】
表２．予測の目的でトレイニングセットおよび検証セットに用いた試料
【００８８】
【表２】
2013510575.tif 000004
【００８９】
[00107]  図１に、本発明者らは、全遺伝子を用いた、ＦおよびＮの全試料（６５試料）、トレイニングセットの試料（３３試料）および検証セットの試料（３２試料）のクラスタリングを示す。類似度の距離（metric of similarity）としてピアソンの相関（Pearson's correlation）を用いる平均連鎖法（average linkage method）を利用し、行正規化（row normalization）に従って、試料を階層的クラスタリング法（hierarchical clustering）でクラスタリングする；Sneath, P. (1973) Numerical taxonomy; the principles and practice of numerical classification. W. H. Freeman, San Francisco, CA USA。
【００９０】
[00108]  クラスタリング過程で、チップ上の完全転写プロファイリング、すなわち５４Ｋ＋転写体すべてを用いる。３データセットすべてのクラスタリングが、妊娠成功に基づく分離がないことを指摘している。これは、発現結果を成功と相関させる表現型特異的遺伝子の同定を見出して最終的に予測遺伝子セットを同定するためには、管理下での学習分析が必要であることを示唆する。
【００９１】
[00109]  成功と相関する遺伝子を見出すために、ｔ−検定に基づく差次発現を示す遺伝子（多重仮説検定のためにボンフェロニ補正（Bonferroni correction）してｐ＜０．０５）をトレイニングセット（Ｆ対Ｎ）において同定した。得られた“記述遺伝子（descriptive gene）”と呼ばれる１１８０の遺伝子を用いてトレイニングセットおよび検証セットを個別にクラスタリングした（図２）。その結果は、特にトレイニングセットにおいてＮとＦの試料の分離に成功したことを示す；このセットの試料は、表４中の妊娠シグナチャー遺伝子ならびにＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＫＲＡＳ、ＮＣＡＭ１、ＮＤＮＬ２、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子を同定するために用いられたからである。
【００９２】
[00110]  次いで、これらの１１８０の遺伝子を用いて１点排除交差検証法（leave-one-out-cross-validation）（Ｌ１ＯＸＶ）をトレイニングセットにおいて実施した。この方法では、予測遺伝子セット中の最初の遺伝子数（複数）、たとえばＰを固定する。次いでトレイニングセット中の１試料を排除し、残りの試料を用いてＮとＦを識別するトップＰ個の遺伝子（複数）を計算する。これらＰ個の遺伝子（複数）を用いて、排除される試料をＮまたはＦと予測する。この操作をトレイニングセットの３３試料すべてについて繰り返す（一度に１つを排除する）。適正予測の総数をトレイニングセットにおける予測子の精度として挙げる。
【００９３】
[00111]  Ｌ１ＯＸＶ過程で種々の数値のＰ、すなわち予測遺伝子数を試みて、良好なＬ１ＯＸＶ予測精度を示すものについて、これらの遺伝子を検証セットに適用する。検証セットにおいて適正に予測された試料数を検証セットにおける予測精度としてレポートする。高いトレイニング精度および検証精度が得られる最小のＰ、すなわち精度グラフが平坦になるＰを予測遺伝子セットとしてレポートする。用いた予測アルゴリズムは、重み付き投票法（Weighted Voting）(Golub et aL, Molecular classification of cancer: class discovery and class prediction by gene expression monitoring, Science, 286, 531-537(1999)、ベイズ複合共変量法（Bayesian Compound Covariate）;（Wright et 1., Proc Natl Acad Sci U S A, 100, 9991-9996 (1999))；対角線形識別法（Diagonal Linear Discriminant）（Dudoit et aL, citation, 97, 77-87)；最近接重心法（Nearest Centroid）、ｋ−最近接法（k-Nearest Neighbors）(Golub et al., Molecular classification of cancer: class discovery and class prediction by gene expression monitoring, Science, 286, 531-537 (1999))；シュランケン重心法(Shrunken Centroids）(Tibshirani et al., Proc Natl Acad Sci U S A, 99, 6567-6572 (2002))；支援ベクターマシーン（Support Vector Machines）(Radmacher, et al., J Comput Biol, 9, 505-511 (2002))、および複合共変量法（Compound Covariate）(Hedenfalk et aL, N Engl J Med, 344, 539-548 (2001))であった。
【００９４】
[00112]  本発明者らの分析では重み付き投票法が２２７遺伝子の予測子セットについて最良の予測を行ない、９７％のＬ１ＯＸＶ精度（３２／３３の適正予測）および８８％（２８／３２適正予測）の検証セット精度を与えた。“予測遺伝子”と呼ぶこれら２２７の遺伝子は、フィッシャー検定（Fisher's test）を用いてトレイニングセットと検証セットの両方に関して有意（ｐ＜０．０５）の予測結果を与えた。予測結果を表３に示す。
【００９５】
[00113]  表３．トレイニングおよび検証データセットについての予測結果
【００９６】
【表３】
2013510575.tif 000005
【００９７】
[00114]  図３に、２２７の予測遺伝子リストを用いたトレイニングセットおよび検証セットのクラスタリングを示す。トレイニングセットにおいて不適正予測された試料のみはクラスタリングにおいても誤置される（misplaced）が、検証セットクラスタリングにおけるＦおよびＮ試料の混合挙動は重み付き投票法によりなされた貢献を強調することに注目されたい。
【００９８】
[00115]  以上に基づいて、本発明者らは、卵丘細胞におけるその発現を評価できる一組のヒト遺伝子（表４に含まれる）を同定し、その卵丘細胞に随伴する卵母細胞（または被験卵丘細胞を単離するために用いた同じドナーに由来する卵母細胞）の妊娠潜在性を評価するために、正常な卵母細胞に相関するそれらの遺伝子の発現レベル、すなわち生存可能な妊娠を生じることができる卵母細胞に随伴するものと比較した。追加試料の分析によって、その発現を卵丘細胞において評価することができる好ましい１４遺伝子のセット（表１２に含まれる）を同定し、その卵丘細胞に随伴する卵母細胞（または被験卵丘細胞を単離するために用いた同じドナーに由来する卵母細胞）の妊娠潜在性を評価するために、正常な卵母細胞に相関するそれらの遺伝子の発現レベル、すなわち生存可能な妊娠を生じることができる卵母細胞に随伴するものと比較した。
【００９９】
[00116]  これらの遺伝子のリストは臨床関連のものであるが、より最適な妊娠シグナチャー遺伝子セットを同定するために、これら２２７および１４の予測遺伝子のリストを機能分析し、さらに検証して、さらに改良できると本発明者らは期待する。
【０１００】
[00117]  それに関して、この表４に含まれる予測遺伝子のリストおよび表１２に含まれる１４遺伝子のセットを改良するために使用できる下記の手法を以下に記載する。
[00118]  本発明者らはこれら２２７または１４の予測遺伝子（ＰＧ’ｓ）に対する有意性を、２つの戦略を用いて指定する：ｉ）クラスラベルを並べ替え(permute)、先に適用したものと同じ方法を用いて最適な２２７遺伝子の予測子を同定し、この摂動(perturbed)データセットを用いてそれらの性能を計算する；ｉｉ）元のデータセットを用いて、ランダムに選択した２２７遺伝子の予測子の性能を検定する。上記の戦略を用いて得た結果に対するＰＧ’Ｓの性能比較を用いて、ＰＧ’Ｓの有意性を評価する。ＰＧ’Ｓを改良するために、完全データセットを別のトレイニングセットと検証セットに分け、それぞれの区画(split)について最適な予測遺伝子セットを計算する。これらの予測子を相互およびＰＧ’Ｓと比較して、一貫して良好な予測子として現われる遺伝子からなる最終的な予測遺伝子セットを得る。次いでこの最終セットをそれの有意性および精度について評価する。最後に、さらにこの最終的な予測遺伝子セットを下記について分析する：ｉ）有意に過剰呈示された（over-represented）ＧＯカテゴリーを見出すための遺伝子個体学（Gene Ontology）（ＧＯ）機能分類、ｉｉ）既知の生物学的経路、遺伝子ネットワーク、疾患経路、小分子と薬物標的の相互作用における関与、ｉｉｉ）共通／個別の転写調節エレメントを同定するためのプロモーター領域分析、およびｍｉＲＮＡ標的分析、ならびにｉｖ）それらの局在性、分泌潜在性、およびライゲーション特性。以下において、本発明者らはこの全作業流れを詳細に説明し、それを図４にまとめる。
【０１０１】
[00119]  本発明者らは、ハイスループット転写プロファイリングおよびプロテオミクスデータセットの両方を用いて糖尿病、脊髄異形成症候群、慢性肝疾患、および種々の形態の癌の予測バイオマーカーを見出すのに、本明細書に定める方法を以前に適用した(参照: Aivado et al., Proc Natl Acad Sci U S A, 104, 1307-1312 2007); Jones et al., Clin Cancer Res, 11, 5730-5739 (2007); Jones, et al., J Urol, 179, 730-736 (2007); Out et al. Diabetes Care, 30, 638-643 (2007); Out et al., J Biol Chem, 282, 11197-11204 (2007); Prall et al., Int J Hematol, 89, 173-187(2009); Spentzos et al., J Clin Oncol, 23, 7911-7918 (2005)、およびZinkin et al, Clin Cancer Res, 14, 470-477 (2008)。これらの研究すべてにおいてこれらのバイオマーカーは有意の予測性をもつことが認められ、この予測子セットを定めるために用いられていなかった新しいデータセットについて実験と解析の両方により独立して検証された。
【０１０２】
[00120]  予測アルゴリズム
[00121]  “信号対ノイズ比”（ＳＮＲ）は、遺伝子ｇの予測値を評価するために用いられる(Golub et al., Science, 286, 531-537 (1999).)。μＦ（ｇ）およびμＮ（ｇ）を、それぞれＦ（妊娠成功）およびＮ（妊娠失敗）試料グループの遺伝子ｇの平均値とする。同様に、σＦ（ｇ）およびσＮ（ｇ）を、それぞれＦおよびＮ試料グループの遺伝子ｇの標準偏差とする。ＳＮＲ（ｇ）＝［μＦ（ｇ）−μＮ（ｇ）］／［σＦ（ｇ）＋σＮ（ｇ）］と定義する。この算術式は両グループについて理想遺伝子発現ベクター付近の近隣ＲＭを規定し、ここでＭ＝｜Ｆ｜＋｜Ｎ｜、すなわちデータセット中の試料の総数である。ＳＮＲはいずれかのグループで著しく偏った発現を伴う遺伝子を処罰し（punish）、遺伝子のメンバーシップについての符号付きランキング法（signed ranking method）を提供する。この場合、大きな正の数値はＦグループについて良好な予測子の指標となり、大きな負の数値（絶対値で）はＮグループについて良好な予測子の指標となる。本発明者らはまた、理想化発現パターンと特定遺伝子ｇの間の相関性の境界をＢ（ｇ）＝［μＦ（ｇ）＋μＮ（ｇ）］／２と定義する。
【０１０３】
[00122]  この方法で、Ｐ個の遺伝子の予測遺伝子セットＧ＝｛ｇ１，ｇ２，．．．，ｇＰ｝、ＦおよびＮ試料のグループ、ならびに予測すべき新しい試料Ｓを評価する。ｇｉ（１＜ｉ＜Ｐ）の投票（vote）をＶｉ＝ＳＮＲ（ｇｉ）［Ｓ（ｇｉ）−Ｂ（ｇｉ）］と定義する；ここで、Ｓ（ｇｉ）はＳ中の遺伝子ｇｉの信号値を表わす。Ｖｉは、Ｓ（ｇｉ）がいかに良好にＦおよびＮ試料のｇｉの“挙動”に関係しているかを表わす。Ｖｉが正であればｇｉに基づいてＳはＦであると予測されると結論し、Ｖｉが負であればｇｉはＳをＮであると予測する。予測子セット中のすべての遺伝子について繰り返して、Ｐ個の投票を求め、ＶＦをすべての正の投票の和とし、ＶＮをすべての負の投票の和とする。絶対値でＶＦがＶＮより大きければ、試料ＳをＦと予測する；そうでなければ、ＳをＮと予測する。本発明者らの先の試験で、トレイニングセットを用いて頑健な予測遺伝子が得られ、これを独立した検証セットについて試験した。
【０１０４】
[00123]  予測遺伝子セットの有意性
[00124]  頑健な予測遺伝子セットを得るために、試料をトレイニングセットと検証セットにランダムに分ける。本発明者らは予め重み付き投票法を用いて交差検証を排除して（Ｌ１ＯＸＶ）、２２７の予測遺伝子のセットを得た。この分析法で、ＴおよびＶを、それぞれトレイニングセットおよび検証セットに用いたＦおよびＮ試料の数とした。これに関して、トレイニンググループの試料のうち１つを排除し、トレイニングセット中のＴ−ｌ  Ｆ試料およびＮ試料中で最大ＳＮＲをもつＰ個の遺伝子を見出した。これらＰ個の遺伝子を用いて、前記の方法により、排除した試料をＦまたはＮのいずれであるか予測する。適正予測の数がトレイニングセットに関するＬ１ＯＸＶ予測精度になり、最高予測精度をもつ最小のＰを続いて検証セットについて試験する。本発明者らの予備試験で、Ｐは２２７であることが見出され、その予測精度はそれぞれトレイニングセットおよび検証セットについて９７％および８８％であった。
【０１０５】
[00125]  これら２２７の予測遺伝子の統計的有意性も評価する。これを評価するために、２段階戦略を採用する。第１相では、トレイニングセットに用いた試料のクラスラベルをランダムにシャフルする；すなわち、ランダムサブセットのＦ試料をＮと呼び、同様にランダムサブセットのＮ試料をＦと呼ぶ。この摂動データセットを用いて、トレイニングセットに関して２２７の遺伝子の重み付き投票法によるＬ１ＯＸＶ精度を計算する。この設定で、並べ替えたクラスラベルを用いて最高ＳＮＲ値をもつ２２７の遺伝子を選択する。クラスラベル並べ替えをＢ回実施し、ｂ回目の並べ替えＬ１ＯＸＶ精度をＡｂとする。本発明者らの元の２２７遺伝子の予測子セットのＬ１ＯＸＶ精度の有意性をｐ＝１／ＢΣＩ（Ａｂ＞９７％）と評価し、ここで、Ｉは指示関数であり、Ａｂ＞９７％であれば１、その他の場合は０の数値と仮定する。同様に、Ａ’ｂを、検証セットに関してｂ回目の並べ替えを用いて得た２２７遺伝子の予測精度とする。次いで、検証セットに関して本発明者らの２２７遺伝子の予測子セットの予測精度有意性をｐ＝１／ＢΣΙ（Α’ｂ＞８８％）と評価する。
【０１０６】
[00126]  有意性を評価するための第２相では、元のクラスラベルを維持するが、ランダムに２２７遺伝子のセットを取り出し、トレイニングセットに関するそれらのＬ１ＯＸＶ精度および検証セットにおける予測精度を評価する。このランダム選択をＢ回実施し、ＡｂおよびＡ’ｂを、それぞれｂ回目の選択についてのトレイニングセットに関するＬ１ＯＸＶ精度および検証セットにおける予測精度とする。同様に、Ｌ１ＯＸＶ精度についての有意性をｐ＝１／ＢΣＩ（Ａｂ＞９７％）として計算し、検証セットの予測精度についての有意性をｐ＝１／ＢΣΙ（Α’ｂ＞８８％）として計算する。両相においてＢ＝１０００である。
【０１０７】
[00127]  予測遺伝子セットのさらなる改良
[00128]  ２２７遺伝子の予測子セットをさらに改良するために、同様にトレイニングセットおよび検証セットの種々の区画を使用し、得られる遺伝子セットの重なりを調べる。良好に規定された予測子を見出すために、改良プロセスでは本発明者らの全データセットをトレイニングセットと検証セットに半分に分けるのではなく；むしろ７５％の試料をトレイニングセットに用い、２５％の試料を検証セットに用いる。さらに、Ｌ１ＯＸＶの代わりに１０倍交差検証戦略を適応する。この場合、トレイニングセットに関して予測子を見出す際、一度に１つの試料を排除するのではなく；代わりにトレイニングセット中の試料の総数の１０％を排除する。各反復に際して残り９０％の試料を用いて予測子セットを計算し、この予測子セットで、排除した１０％の試料を予測する。次いで適正予測の総数を用いて予測精度を計算する。
【０１０８】
[00129]  各データ区画について、検証セットに関する予測精度および予測遺伝子セットの有意性を前記に従って評価する。次いで各区画で見出した予測遺伝子セットを相互に、および元の２２７遺伝子の予測子セットと交差させて、一貫して良好な予測子として現われる遺伝子を見出す。このために、交差検証法を用いて、全データセットに関してより改良された予測遺伝子セットを形成し、それの予測精度および有意性を計算する。
【０１０９】
[00130]  クラスタリングおよび機能分析
[00131]  算術平均付き重みなし対グループ法（Unweighted Pair Group Method with Arithmetic-mean）（ＵＰＧＭＡ）（Sneath, P. (1973) Numerical taxonomy; the principles and practice of numerical classification. W. H. Freeman, San Francisco, CA USA）、類似度木（similarity tree）構築のために用いた階層的クラスタリング（hierarchical clustering）法、および主成分分析（principal components analysis）（ＰＣＡ）(Otu, et al., J Biol Chem, 282, 11197-11204 (2007))を用いて、試料および遺伝子のクラスタリングを実施する。階層的クラスタリング法では、平均連鎖クラスタリング（average linkage clustering）を適用する前に発現データマトリックスを各遺伝子について行正規化（row-normalized）し、ピアソンの相関（Pearson's correlation）を距離の尺度として用いる。元の分散を可能な限り多く維持しながら多変量データ対象物をより低次元空間内へ射影するＰＣＡでは、各試料を平均ゼロおよび標準偏差１まで正規化する。
【０１１０】
[00132]  機能分析は、それ以上の調査の根拠となる目的遺伝子リスト中に遺伝子個体学（Gene Ontology）（ＧＯ）カテゴリーを見出すことからなる(Ashburner, et al., Nat Genet, 25, 25-29 (2000))。発現分析の系統的探索法（Expression Analysis Systematic Explorer）（ＥＡＳＥ）は、そのセットにおいて過剰呈示され（over-represented）、したがってさらに関心がもたれる可能性のある、生物学的関連カテゴリーを同定する(Hosack et al., Genome Biol, 4, R70) (2003)。
【０１１１】
[00133]  これを達成するために、ＥＡＳＥは各プローブをＧＯカテゴリーに付随するEntrez Gene identifier (Maglott, et al., Nucleic Acids Res, 35, D26-31 (2001))にマッピングする。GO Consortium (Geneontology. org)は、各遺伝子を（適用可能な場合は）３つのＧＯカテゴリー中の１以上のクラスに指定する：生物学的機能、細胞プロセス、および分子機能。ＥＡＳＥは入力遺伝子リスト中の過剰呈示されているＧＯカテゴリーを、フィッシャー正解確率（Fisher exact probabilities）のジャックナイフ反復リサンプリング（jacknife iterative resampling）を用いて同定し、ボンフェロニ多重検定補正（Bonferroni multiple testing correction）する。“ＥＡＳＥスコア”はジャックナイフフィッシャー正解確率の分布の上限である。少ない数の遺伝子を含むカテゴリーは標準重みに達しない（under-weighted）ので、ＥＡＳＥスコアはフィッシャー正解検定より頑健である。ＥＡＳＥ分析は予測遺伝子リストをチップ上の全遺伝子に対して検定し、ＥＡＳＥスコアは入力リスト中のＧＯカテゴリーの過剰呈示の尤度（likelihood）について計算される。過剰呈示は、与えられた入力リスト中に、分布がランダムであれば出現するであろうと予想されたより頻繁に現われる特定のＧＯカテゴリー、たとえば細胞周期、に属する遺伝子のグループを記述する。ＥＡＳＥスコアは有意性レベルであり、小さいＥＡＳＥスコアほど過剰呈示における信頼性が増すことの指標となる。０．０５以下のＥＡＳＥスコアをもつＧＯカテゴリーを有意に過剰呈示しているとして選択する。
【０１１２】
[00134]  経路分析
[00135]  経路分析、機能強化分析（functional enrichment analysis）、または遺伝子セット分析は、マイクロアレイ試験で有意に調節される前決定した遺伝子セットまたはクラスに注目する。Ingenuity Software Knowledge Base （ＩＫＢ）(カリフォルニア州レッドウッドシティー)を用いて、基礎となる転写調節を最も良く説明するネットワークおよび経路を同定する。ＩＫＢは、頑健な相互作用データベースを提供する科学文献の手動管理（manual curation）に基づく、遺伝子および／または遺伝子産物間の相互作用を利用する。目的とする遺伝子リストが分析されると、結果を２方法で考察する：ｉ）入力遺伝子リストおよび限定的な拡張を用いて形成されたネットワーク、ｉｉ）入力遺伝子のサブセットを有意に宿す既知の生物学的経路。両方の結果を、最終遺伝子ネットワークと部分的に相互作用することが知られている薬物、低分子代謝産物、機能、および疾患に関してさらに分析することができる。この方法で、本発明者らの予測遺伝子セットによって最も良く説明される経路を同定し、これらのネットワークに対するさらなる摂動を推論する。さらに、前決定した正準経路（canonical pathway）に必ずしも属さない新しい相互作用ネットワークを作成することができる。
【０１１３】
[00136]  プロモーターおよびｍｉＲＮＡ標的の分析
[00137]  本発明者らの予測シグナチャー中の遺伝子の５ｋｂ上流を調べ、同定したこれらの領域を、プロモーター分析および相互作用ツール（promoter analysis and interaction tool）ＰＡＩＮＴ(Vadigepalli et al, OMICS, 7, 235-252 (2003))およびｏＰＯＳＳＵＭ(Ho Sui et al., Nucleic Acids Res, 33, 3154-3164 (2005))を用いて分析する。これらの方法は、入力遺伝子リストのプロモーター領域のモデルを作成し、これらの領域の転写調節エレメント（trascription regulatory element）（ＴＲＥ）を同定する。次いでそれらの結果をＴＲＥの過剰呈示について分析し、基礎となるマイクロアレイデータに関して遺伝子調節に重要な潜在的転写因子を見出す。
【０１１４】
[00138]  予測遺伝子のもうひとつの潜在的調節機序はｍｉＲＮＡによるものであり、これは主に転写後抑制またはｍＲＮＡ分解によって配列特異的な様式で遺伝子発現を調節する。ＴａｒｇｅｔＳｃａｎＳ(http://www.targetscan.org/)を用いて、予測遺伝子についてｍｉＲＮＡ標的部位を同定することができる。ＴａｒｇｅｔＳｃａｎＳはＴａｒｇｅｔｓｃａｎの改良バージョンであり、ＵＴＲを２〜８塩基のｍｉＲＮＡに対する完全ワトソン−クリック相補セグメントについて検索し、各ｍｉＲＮＡ標的部位相互作用についてＲＮＡｅｖａｌによりフォールディング自由エネルギーＧを計算し、次いで各ＵＴＲにＺスコアを指定する(Lewis et al., Cell, 115, 787-798 (2003))。こうして、入力遺伝子リストにおけるｍｉＲＮＡ標的が、ｍｉＲＮＡ結合部位のコンピューター処理により同定される。
【０１１５】
[00139]  分泌タンパク質および膜リガンドの予測
[00140]  細胞内でのタンパク質の局在性および輸送は、それらのアミノ酸配列中の固有のシグナルにより支配されている。本発明者らの予測遺伝子リストを、これに関してそれらがコードするタンパク質により分析することができる。タンパク質のアノテーションがそれの局在性、分泌ポテンシャル、ライゲーションその他の特性について情報を伝えるのに十分な場合、標的をそれ以後の試験のために単離するには手動でのフィルタリングで十分であろう。実質的なアノテーションがない場合、高精度コンピューターアルゴリズムＴａｒｇｅｔＰを用いてタンパク質の上記特性を予測することができる；これは、神経回路網を利用し、タンパク質のＮ末端の特徴を考慮に入れる(Emanuelsson et al. Nat Protoc, 2, 953-971 (2007))。分泌タンパク質を予測する場合、タンパク質の分泌ポテンシャルを低下させる膜貫通ドメインをもつものをフィルタリングで除去することが役立つ。これらのタンパク質における膜貫通ドメインの予測は、予測のために隠れマルコフ・モデル（hidden Markov model）を用いる独立バージョンのＴＭＨＭＭを用いて行なわれる(Sonnhammer et al., Proc Int Conf Intell Syst Mol Biol, 6, 175-182 (1998))。ＰＳＯＲＴ９に規定された方法によりこの分析を補足することもできる（Horton et al, Nucleic Acids Res, 35, W585-587 (2007))；これは、それの予測戦略において、ｋ−最近接分類法、およびギャップなし多重アラインメントから得られる重み付きマトリックスを使用する。
【０１１６】
[00141]  別の戦略
[00142]  さらに、場合により２方法で本発明者らの分析を向上させることができる。
[00143]  信号値の計算および正規化：モデルを基礎とする方法、たとえばｄＣｈｉｐ(Li et al., Genome Biol, 2, RESEARCH0032 (2001).)およびＲＭＡ(Irizarry et al., Nucleic Acids Res, 31, el5. (2001))を用いて遺伝子チップデータを正規化および要約することができ、これらはAffymetrix MAS 5.0正規化(Hubbell et al., Bioinformatics, 18, 1585-1592.(2002))より優れていることが示されたが、本発明者らはこれらの方法の使用を避けた；それらはベースライン試料に依存し、新しい試料の追加に適応しないからである(Barash et al., Bioinformatics, 20, 839-846 (2004))。言い換えると、新しい試料をデータセットに追加するので、それ以前のモデルに基づく信号分析は無効になり、したがって完全な分析ワークフローが繰り返される。本発明者らのデータセットは絶えず増大し、以前の所見を新たに追加した試料において検証する必要があるので、本発明者らはＭＡＳ  ５．０、すなわち新しい試料の追加に対して頑健な信号値計算および正規化方法を用いる。AffymetrixからのExpression Console中のProbe Logarithmic Intensity Error（ＰＬＩＥＲ）ワークフローと呼ばれる新しい方法(www.affymetrix.com)は、四分位点正規化（quintile normalization）を含み、プローブアフィニティー、経験的プローブ性能の利用によって、また低濃度および高濃度付近の誤りを適宜処理することによって、改良された信号値を提供する(Katz et al., BMC Bioinformatics, 7, 464 (2006)。ＰＬＩＥＲは新しい試料の追加に対してＭＡＳ  ５．０ほど頑健ではないが、それはモデルを基礎とする方法と比較してより良好な性能を備えており、したがって信号値計算および正規化のための本発明者らの別の戦略である。
【０１１７】
[00144]  予測アルゴリズム：
[00145]  本発明者らの先の試験において最良の性能をもつ予測アルゴリズムであった重み付き投票法のほかに、本発明者らの結果が本出願に述べる成功基準に一致しない場合は下記の予測アルゴリズムを試みる計画を立てる：ベイズ複合共変量法（Bayesian Compound Covariate）(Dudoit et al., Journal of the American Statistical Association, 97, 77-87 (2002); Wright et al., Proc Natl Acad Sci U S A, 100, 9991-9996 (2003))、対角線形識別法（Diagonal Linear Discriminant）、最近接重心法（Nearest Centroid）、ｋ−最近接法（k-Nearest Neighbors）(Golub et al., Science, 286, 531-537 (1999)、シュランケン重心法（Shrunken Centroids）(Tibshirani et al., Proc Natl Acad Sci U S A, 99, 6567-6572 (2002)、支援ベクターマシーン（Support Vector Machines）(Radmacher et al., J Comput Biol, 9, 505-511 (2002))、および複合共変量法（Compound Covariate）(Hedenfalk et al., N Engl J Med, 344,539-548 (2001))。
【０１１８】
[00146]  ｑＲＴ−ＰＣＲを基礎とするTaqman（商標）Low Density Arrays（ＬＤＡ）を用いる予測遺伝子セットの検証
[00147]  Taqman ＬＤＡ(Applied Biosystems)は、３８４の遺伝子をごく少量のＲＮＡから高い適合度で同時ｑＲＴ−ＰＣＲ分析できるマイクロフルイドプレート(microfluidic plate）である。単一の特注ＬＤＡを用いて２セットの遺伝子を検証する：１．）予備マイクロアレイ試験で同定したＰＧ’Ｓ、および２．）本発明者らの共同研究者Dagan Wells（オックスフォード大学）が提供したセット。独立した試験で、Wellsらは卵母細胞に異数性で随伴する卵丘細胞に発現する遺伝子のセットを同定した。細胞遺伝学研究により、平均年齢３２の女性において１／４の卵母細胞が異数体であることが明らかになった(Fragouli Cytogenet Genome Res; 114:30-38 (2006)。これらの遺伝子セットの両方を合わせることによって、より最適でより強力な妊娠シグナチャーを同定することができる。
【０１１９】
[00148]  それに基づいて、本発明者らの特注ＬＤＡで処理するために先にマイクロアレイ分析を行なった６５の組から２５のＦ（完全妊娠成功）および２４のＮ（妊娠成功なし）卵丘細胞試料をランダムに選択する。これらのｃＤＮＡ（マイクロアレイ分析から残ったもの）を調製し、ABI 7900HTマシーンでＬＤＡ当たり１試料をApplied Biosystems' Taqman LDA指示に従って処理する。各転写体の絶対定量を実施する。得られた増幅データを7000 System SDSソフトウェア(Applied Biosystems)により分析し、各遺伝子の強度値を対照遺伝子に対して正規化する。
【０１２０】
[00149]  マイクロアレイ結果との比較
[00150]  スピアマンの順位相関（Spearman's rank correlation）およびマン−ホイットニーＵ検定（Mann-Whitney U test）を適用して、試料全体にわたってマイクロアレイおよびＬＤＡにより得た各遺伝子につき発現レベルを比較する。パラメトリックモデルにおける分布仮定は２プラットホームでは有効でない可能性があり、かつ基礎となる手法および正規化法の差のため２プラットホーム間の信号値を比較できない可能性があるので、これらのノンパラメトリック法および順位に基づく方法を採用する。試験した遺伝子それぞれについて、ＬＤＡに用いたＦとＮの試料間の変化倍率も計算する。この変化の方向がマイクロアレイ結果にみられるものと一致し、２プラットホーム間で遺伝子発現プロフィールが有意には変化しないままであり（ｐ＞０．０５）、相関性が高い（ｒ＞０．９）ならば、これらの遺伝子を検証されたとラベル付けする。ＬＤＡ実験に用いた試料に対する検証された遺伝子の予測力を、ＬＤＡとマイクロアレイ両方の信号値を別個に用いて評価する。同様に、検証された遺伝子の予測精度を、マイクロアレイ信号値を用いた完全データセットを用いて評価する。これらの予測プロセスで、試料をトレイニングセットと検証セットに分け（７５％−２５％の区画）、そしてトレイニングセットについてモデルを構築してそれを検証セットについて検査することにより、予測精度を計算する。この予測戦略において、前記の重み付き投票法を適用する。
【０１２１】
[00151]  別の戦略
[00152]  予測遺伝子がＬＤＡプラットホームにおいて前記のように検証されない場合、または検証された遺伝子を用いて良好な予測結果が得られない場合、予測戦略を改変し、予測遺伝子シグナチャーの頑健性を犠牲にすることによって、より高い精度および妊娠予後と良好に相関する遺伝子を得ることができる。この目的のために、元のデータセットをトレイニングセットと検証セットに分けずに予測遺伝子シグナチャーを計算する。代わりに、完全なマイクロアレイデータセットを用いて予測子を構築し、前記に従ってそれの精度を計算し、それの統計的有意性を評価する。この予測戦略に際して、重み付き投票法を用いて完全データセットに適用する１０倍排除交差検証法（leave-ten-fold-out cross validation）を用いる。
【０１２２】
[00153]  他の試料への妊娠予測子の適用
[00154]  １５のＦ（完全妊娠成功）および１５のＮ（妊娠成功なし）の新しい卵丘細胞試料を採集する。前記に従ってＲＮＡを単離する。ABI's High Capacity cDNA Reverse Transcriptionキットを用いて逆転写を完了させ、特注の予備増幅プール（ＡＢＩ）を用いてｃＤＮＡの増幅を完了させる。予備増幅およびｑＲＴ−ＰＣＲサイクリングは、ABI 7900HTマシーンにより、ＬＤＡ当たり１試料で、Applied Biosystems' Taqman LDA指示に従って逐次行なわれる。各転写体の絶対定量を行なう。得られた増幅データを7000 System SDSソフトウェア(Applied Biosystems)により分析し、各遺伝子の強度値を対照遺伝子に対して正規化する。
【０１２３】
[00155]  ＬＤＡプラットホームにおいて検証した最終的な予測遺伝子を用いて、前向き研究における試料からなる完全データセットに関して予測モデルを構築する。次いで、後向き研究で得た信号値を用いて、このモデルを新しい試料に適用する。この予測戦略で、前記に述べた重み付き投票法を適用する。
【０１２４】
[00156]  別の戦略
[00157]  場合によっては、より厳密な信号値を得るために、後向き研究で実施した実験を繰り返すか、あるいは本発明者らの検証努力の統計的パワーを潜在的に増強するために、用いる試料数を増加してもよい。
【０１２５】
[00158]  したがって、以上に基づいて、好ましい態様において本発明方法は、妊娠受容能をもつ卵母細胞を産生する女性対象または産生しない女性対象を、表４に含まれる差示発現する遺伝子のセットもしくは表１２に含まれる１４遺伝子のセットまたは対応するコードされたポリペプチドの発現レベルに基づいて同定するために用いられる。ただし本発明方法は、ヒト以外の動物、たとえばイヌ、コウモリ、ウシ、ウマを含めた他の哺乳動物、鳥類、両生類、爬虫類などにも適用できる。たとえば、本発明を用いて、すなわち表４中の妊娠シグナチャー遺伝子もしくは表１２に含まれる１４遺伝子のセットまたはその相補体のうち１以上の発現を調節する化合物をスクリーニングすることにより、想定雌性受精処置の研究のための動物モデルを誘導することができる。さらに本発明を用いて、妊娠受容能をもつ卵母細胞を産生する能力を排除または阻害する異常、たとえば卵巣機能障害、卵巣嚢胞、閉経前もしくは閉経状態、癌、自己免疫障害、ホルモン機能障害、細胞増殖障害、または妊娠受容能をもつ卵母細胞の発生を阻害もしくは排除する他の状態を伴う雌性対象を同定することができる。
【０１２６】
[00159]  たとえば、特異的な妊娠シグナチャー遺伝子を特徴的な発現レベルで発現しない対象をスクリーニングして、それらが妊娠受容能をもつ卵母細胞を産生するのを排除する基礎健康状態を伴うかどうかを評価する。特に、そのような対象をスクリーニングして、それらが閉経の徴候を示しているかどうか、癌、自己免疫疾患または卵巣異常、たとえば卵巣嚢胞を伴うかどうか、あるいは“妊娠受容能をもつ”卵母細胞を産生するのを排除する可能性のある他の健康状態、たとえばホルモン障害、アレルギー性障害などを伴うかどうかを評価する。
【０１２７】
[00160]  さらに本発明を用いて、ヒトまたはヒト以外の雌性対象における想定雌性受精処置の有効性を評価することができる。本質的に、そのような方法は下記を含むであろう：雌性対象、好ましくは女性を、想定妊娠増強処置により処置し、処置後に少なくとも１個の卵母細胞および随伴する周囲卵丘細胞をその女性から単離し、場合によりさらに処置前に少なくとも１個の卵母細胞および随伴する周囲細胞をその女性から単離し、単離した卵母細胞それぞれから少なくとも１個の卵丘細胞を単離し；妊娠受容能をもつ卵母細胞に随伴する（を囲む）卵丘細胞が特徴的レベルで発現するかまたは特徴的レベルで発現しない少なくとも１つの遺伝子の発現レベルを検出し；その処置が、妊娠受容能をもつ卵母細胞を囲む卵丘細胞に（処置しない場合より）より特徴的なレベルで少なくとも１つの妊娠シグナチャー遺伝子を発現する卵丘細胞を生じるかどうかに基づいて、その想定受精処置の有効性を評価する。前記のように、雌性ヒト対象が好ましいが、本発明方法はヒト以外の雌性動物、たとえばヒト以外の雌性霊長類、または想定ヒト受精処置の評価のための他の適切な動物モデルにおいて、想定受精処置の有効性を評価するために使用できる。
【０１２８】
[00161]  さらにまた本発明は、体外または体内受精処置の有効性を向上させるために使用できる。特に、“妊娠受容能をもたない”または未熟であることが認められた卵母細胞を、体内または体外受精の前、途中または後に、“妊娠シグナチャー”遺伝子であると同定された遺伝子がコードする１以上の遺伝子産物、たとえばホルモン、増殖因子、分化因子などを含有する培地で培養することができる。本質的にこれらの遺伝子産物の存在は、通常は“妊娠受容能をもつ”卵母細胞を囲む卵丘細胞が発現して、それらが普通は卵母細胞に栄養供給し、これによりこれらの卵母細胞が受精した際に生存可能な妊娠をもたらす能力が促進される栄養性遺伝子産物の不足を補うはずである。
【０１２９】
[00162]  あるいは、それらの妊娠シグナチャー遺伝子がコードする１以上の遺伝子産物またはそれらの遺伝子の発現を調節する化合物を、妊娠受容能をもつ卵母細胞を産生しないことが本発明の方法により見出された対象に投与してもよい。そのような投与は、非経口、たとえば静脈内、筋肉内、皮下への注射によるか、あるいは経口または経皮投与によるものであってもよい。あるいは、これらの遺伝子産物を標的部位、たとえば雌性の卵巣または子宮環境に局所投与してもよい。たとえば、雌性対象の子宮または卵巣に、“妊娠シグナチャー”遺伝子がコードする１以上の遺伝子産物またはそのような遺伝子の調節剤を持続送達する薬物送達デバイスを埋め込むことができる。
【０１３０】
[00163]  したがって、一般に本発明は、排卵時の卵子、好ましくはヒトの卵子に由来する卵丘細胞（好ましくは卵丘細胞）が発現する遺伝子であってそれの発現レベルはその卵子が自然受精または人工受精して適切な子宮環境へ移行した際に生存可能な妊娠を生じる能力と相関する遺伝子の、同一性および相対量に関する同定および特性解明を伴う。
【０１３１】
[00164]  １例示態様において、本発明は、卵母細胞の受容能またはある個体が受容能をもつ卵母細胞を産生するかどうかを評価する手段として前記に例示した核酸またはタンパク質の既知の検査法を用いて、卵丘細胞による表４中のいずれかの組合わせ、すなわち少なくとも１、２、３、４、５、６、１０、５０、．．．１００、２００．．．．２２７の遺伝子、または表１２中の１４のいずれかの組合わせの遺伝子の発現を、正常な卵丘細胞が発現するレベルと比較してアッセイする。
【０１３２】
[00165]  ただし、例示態様における本発明は表４中の遺伝子および表１２の遺伝子のいずれかの組合わせを選択するが、好ましい態様において本発明は、表１２中の少なくとも２つの遺伝子、より好ましくは少なくとも３つ、およびこれらの遺伝子１４すべての発現を、卵母細胞の受容能を評価する手段としてアッセイする。さらに本発明方法は、あるいは、差示発現する卵丘細胞が発現する表４および表１２中のものから選択される１以上の遺伝子であってそれらの発現が卵母細胞の受容能と相関する遺伝子の発現レベルを、その発現が同様に予測に役立つことが認められた他の遺伝子と共にモニターすることにより実施できる。
【０１３３】
[00166]  表４：ヒト卵丘細胞が差示発現する、妊娠受容能に関連するヒト遺伝子
【０１３４】
【表４−１】
2013510575.tif 000006
【０１３５】
【表４−２】
2013510575.tif 000007
【０１３６】
【表４−３】
2013510575.tif 000008
【０１３７】
【表４−４】
2013510575.tif 000009
【０１３８】
【表４−５】
2013510575.tif 000010
【０１３９】
【表４−６】
2013510575.tif 000011
【０１４０】
【表４−７】
2013510575.tif 000012
【０１４１】
【表４−８】
2013510575.tif 000013
【０１４２】
【表４−９】
2013510575.tif 000014
【０１４３】
【表４−１０】
2013510575.tif 000015
【０１４４】
【表４−１１】
2013510575.tif 000016
【０１４５】
表５：差示発現データを検証するために用いたヒト遺伝子
【０１４６】
【表５−１】
2013510575.tif 000017
【０１４７】
【表５−２】
2013510575.tif 000018
【０１４８】
【表５−３】
2013510575.tif 000019
【０１４９】
【表５−４】
2013510575.tif 000020
【０１５０】
【表５−５】
2013510575.tif 000021
【０１５１】
【表５−６】
2013510575.tif 000022
【０１５２】
【表５−７】
2013510575.tif 000023
【０１５３】
[00167]  表４に示すように、本発明者らはＩＶＦ成功の指標となる２２７の遺伝子および４１２８の差示発現する遺伝子を同定した。ＴＬＤＡチップ上にこれらの遺伝子のうち１９９を呈示し、１９６のユニークなものを呈示した（Ｂ３ＧＮＴ７、ＡＴＰ８Ａ１およびＤＮＡＪＣ１５は２プローブセットにより呈示される）。これら１９６のうち１４１は予測セットに属し（Ｅｘｃｅｌファイルの“ｇｒｏｕｐ（グループ）”の欄にＰにより呈示する）、５５は差示発現セットに属する（差示発現遺伝子を見出すために全試料を用いたので、“Ａｌｌ”にちなんでＡにより呈示する）。
【０１５４】
[00168]  ＴＬＤＡ実験の出力は閾値サイクル（Ｃｔ）値であり、これは“蛍光が閾値を通過する分画サイクル数（fractional cycle number）”である。これらの値は対数尺度（底２）であり、発現と反比例する。たとえば、遺伝子ＸおよびＹに関するＣｔ値をそれぞれ１８および１９とする。これはＸを１８サイクル目で“良好に”“検出する”ことを意味し、これはＹ（１９サイクル目で良好に検出される）と比較して２倍大きい発現レベルとなる。したがって、ＸとＹの間の古典的変化倍率（Ｘ／Ｙ）を計算するためにはＸをＹから差し引いてそれを２の冪数とする必要がある。
【０１５５】
[00169]  アレイ全体にわたる正規化を提供するスパイクイン対照（spike-in control）はない。おおよその正規化を行なうために、全遺伝子のＣｔ値を試料全体について考慮し、最も安定な発現を選択する。“安定な”遺伝子（そのプレート内で）のＣｔ値を各遺伝子のＣｔ値から差し引くことにより、プレート（試料）内の各遺伝子についてデルタＣｔ（ｄＣｔ）値を計算する。Ｓｔａｔｍｉｎｅｒ分析により、“安定な”遺伝子はＧＡＰＤＨであることが見出された。
【０１５６】
[00170]  本発明者らの“予測分析に用いた古い”試料からＴＬＤＡに用いた試料を下記の表６に挙げる。表中の“＊”は、先の分析で不適正予測された試料に相当する。
[00171]  表６．先の分析からＴＬＤＡ実験に用いた試料
【０１５７】
【表６】
2013510575.tif 000024
【０１５８】
したがって、１９のＦ試料および１６のＮ試料、合計３５がある。１４の新しい試料についても実験した：表７に示す１１のＮ、３のＦ。
表７：ＴＬＤＡ実験に用いた新しい試料
【０１５９】
【表７】
2013510575.tif 000025
【０１６０】
[00172]  データの質
[00173]  “検出”呼出しの遺伝子について分布を検査した。遺伝子の発現レベルが決定されなければ、それに４０の数値を指定する。図５に、各遺伝子について４０の数値を伴う試料の数を、本発明者らの遺伝子（“Ｈａｓａｎ遺伝子”と表示した１９６の遺伝子）およびＴＬＤＡ上の全３７９遺伝子（“全遺伝子”と表示）に対して個別にプロットしたものを示す。この図には、各遺伝子についてその遺伝子が４０の数値をもつ試料の数を示す。この図では、全試料について結果を計算している。
【０１６１】
[00174]  本発明の遺伝子について全体として、６８６０（１９６＊３５）の数値のうち２２３７の数値が４０であり、これは３２．６％の測定で検出されない遺伝子があることを示す。若干の遺伝子は若干の試料中にほとんど存在しないと予想されるので、本発明の定量分析では４０の数値をもつ遺伝子があっても許容される。これに対し、本発明方法では３５試料すべてに検出されない遺伝子があることは許容されない。たとえば、本発明の遺伝子について、３０以上の試料に検出されない遺伝子が６２ある。全遺伝子を考慮すると、４０の数値を伴う測定の割合は約２７．４％である：１３２３０（３７８＊３５）のうち３６３１。
【０１６２】
[00175]  図６に、４０と呼び出された遺伝子の数を各試料について示す。これは、ある意味において、各試料について得た検出レベルを示す。この図に、４０の数値を伴う１１７のＨａｓａｎ遺伝子および２１２の全遺伝子をもつ９＿０７２４０７のように、最適下の全体的検出を伴う若干の試料がみられる。試料全体において４０と呼び出された遺伝子の数の平均±標準偏差値は、Ｈａｓａｎ遺伝子については６３．９±１４．５、全遺伝子については１０３．７±２８．６である。
【０１６３】
[00176]  さらに、１４の“新しい”試料について同様な分析を行なった。その際、Ｈａｓａｎ遺伝子においては２７４４（１４＊１９６）の数値のうち５４１が４０であり（測定の１９．７％に相当）、全遺伝子においては５２９２（１４＊３７８）の数値のうち８９６が４０である（測定の１６．９％に相当）ことが認められる。これらは“古い”試料においてみられるものより１０％以上良好である。図７および８に、それぞれ遺伝子および試料に関してデータを分析した場合の、４０の数値を伴う試料の数および遺伝子の数を示す。
【０１６４】
[00177]  図７は、各試料について、遺伝子が４０の数値をもつ試料の数を示す。結果は１４の試料について計算されている。
[00178]  本発明の遺伝子リストについて、４０の数値をもつ試料が１２以上ある２５の遺伝子がある。図８に、各試料について４０の数値をもつ遺伝子の数を示す；結果は１４の“新しい”試料について計算されている。
【０１６５】
[00179]  その際、全体的検出が不良である１試料１７＿１００７０９＿２２がみられる；これは、古いデータセットにおいて４０の数値をもつ、卵丘細胞に差示発現する“Ｈａｓａｎ遺伝子”と呼ぶ８５のヒト遺伝子のセットおよび１４０の遺伝子をもつ。新しいデータセット中では、試料全体において４０と呼び出された遺伝子の数の平均±標準偏差値は、Ｈａｓａｎ遺伝子については３８．６±１３．７、全遺伝子については６４．０±２２．８であり、これは前のデータセットと比較して改良されている。
【０１６６】
[00180]  これらの結果は、新しい試料についてのＴＬＤＡ検出は古い試料と比較して改良されており、遺伝子に対するある種のフィルタリングが必要であることを示唆する。このフィルタリングは、遺伝子が４０の発現値をもつ、すなわち検出されない試料の数を基礎とすることができる。前記に述べたように、たとえば前の試料については試料のグループに関係なく８５％を超える試料においてその遺伝子が検出されない多数の遺伝子がある。
【０１６７】
[00181]  データ分析
[00182]  ４０の発現値をもつデータ点の取扱い、変化倍率の計算、および対数化した値を用いるか用いないかなど、考慮すべきさまざまなことがある。以下においてそのようなことに対処して、有望な解決策を提示する。
【０１６８】
[00183]  スケーリング：２セットの出力が得られる：Ｃｔ値（対数化した発現レベル）およびｄＣｔ値；その際、与えられた試料について、各遺伝子のｄＣ値はＧＡＰＤＨのＣｔ値をその遺伝子のＣｔから差し引くことにより計算される。Ｃｔ値は対数であるので、これは各遺伝子の発現値をＧＡＰＤＨの発現値で割ることに相当する。言い換えると、それはある遺伝子とＧＡＰＤＨの間の変化倍率である。これらの数値を用いて進めることは、ＧＡＰＤＨに対する遺伝子の変化倍率に基づいてグループ間の変化倍率を計算することを意味する。ＧＡＰＤＨはアレイに用いる内在対照のひとつではないので、ＴＬＤＡに用いるスパイクイン対照はなく、対数目盛りでの小さな変動が実際の数値において大きな差を示す場合があり、ある程度注意してこれを調べる。それにもかかわらず、スケーリングした数値とスケーリングしていない数値の両方を用いた分析を提示する。このレポートの残りの部分について、スケーリングしていない数値は増幅ファイルで得たＣｔ値を表わし、スケーリングした数値はＧＡＰＤＨを差し引くことにより得たｄＣｔ値を表わす。
【０１６９】
[00184]  変化倍率：
[00185]  ２つの試料ＡおよびＢがあり、これらの試料における遺伝子Ｘの発現値がそれぞれａＸおよびｂＸとする。ＴＬＤＡ出力においてみられるもの（Ｃｔ値）は、ｌｏｇ（ａＸ）およびｌｏｇ（ｂＸ）である。これら２つの試料間の変化倍率を計算したい場合、Ｃｔ値を引算し、それを２の冪数とする。すなわち、ＦＣ＝２  ｌｏｇ（ａＸ）−ｌｏｇ（ｂＸ）。その理由は、下記の規則である：ｌｏｇｐ−ｌｏｇｑ＝ｌｏｇ（ｐ／ｑ）、および２  １ｏｇ２ｐ＝ｐ。ただし、Ｃｔ値は逆転しており、すなわち小さい数値ほど大きい表示を意味するので、このＦＣは変化倍率Ｂ／Ａを与える。たとえば、Ａにおいて１０．８、Ｂにおいて１２．３のＣｔ値がみられるとすれば、これはこの遺伝子がＡにおいてアップレギュレートされており、Ｂ／Ａの変化倍率は２  １０．８−  １２．３＝２  −１．５＝０．３５であることを意味する。言い換えると、この遺伝子はＡにおいて１／０．３５＝２．８倍、アップレギュレートされている。この点に到達するもうひとつの方法は、まずＣｔ値を対数解除し（unlog）、次いで既知の方法でＦＣを計算する；ただし、方向が逆転しており、すなわちＣｔワールドでは小さいほど大きいことを意味する。したがって、Ａについての発現レベル＝２  １０．８＝１７８２、Ｂについての発現レベル＝２  １２．３＝５０４２、およびＦＣはＢ／Ａ＝１７８２／５０４２＝０．３５となる。
【０１７０】
[00186]  １未満のＦＣ値は解釈が難しいので、行なうことはそれらを逆転させ、負の記号を付けることである。上記の例については、Ｂ／ＡについてのＦＣが０．３５であると言う代わりに、Ｂ／ＡについてのＦＣが−１／０．３５＝−２．８であると言う。すべての計算において、常にＦ値をＮ値から差し引く（ｌｏｇ目盛りを用いたとすれば）か、あるいはＮ値をＦ値で割って（ｌｏｇ解除した値を用いたとすれば）ＦＣをＦ／Ｎについて計算した。上記に例示したように、１未満のＦＣを表現するためには負の数値を用いた。
【０１７１】
[00187]  単純なＦＣを計算することはさほど複雑ではないとしても、さらに考慮すべきことがある。上記の例は２つの試料のみを含み、すなわちそれは各グループに１つの試料があるとみることができる。本発明の場合のように各グループに１つより多い試料があればどうであろうか（１６Ｎ，１９Ｆ）？Ｃｔ値を平均すれば、実際に発現レベルの幾何平均が得られる。次いで２グループのＣｔ値の平均を引算し、次いで２の冪数とすると、これは２グループの幾何平均を割ることによりＦＣを計算することを意味する。その理由は、下記の規則である：a logX = logX aおよびlogp + logq = log(pq)。
【０１７２】
[00188]  一例をあげると、発現レベルａ、ｂおよびｃがグループＮにあり、ｄ、ｅ、ｆおよびｇがグループＦにあると仮定する。ＴＬＤＡ出力にみられるのはｌｏｇａ、ｌｏｇｂなどである。ＦＣ（Ｆ／Ｎ）を計算するために、Ｆにおける平均値をＮにおける平均値から差し引き、次いで２の冪数とすると、下記が得られる。
【０１７３】
[00189]  Ｎにおける平均値 = 1/3 [loga + logb + logc] = l/3 1og[abc] = log(abc) l/3
[00190]  Ｆにおける平均値 = 1/4 [logd + loge + logf + logg] = l/4 1og[defg] = log(defg) l/4
[00191]  FC (F/N) = 2 ＾ [log(abc) l/3 - log(defg) l/4] = 2 ＾ {log[(abc) l/3 / (defg) l/4]}= (abc) l/3 / (defg) l/4
[00192]  数字の幾何平均はそれらの積のｎ乗根であることを思い起こされたい。したがって、我々は常に対数解除した数値を用いて作業することを選択している。すなわち、本発明者らはまずＣｔ値を２の冪数とし、次いで分析を行なった。
【０１７４】
[00193]  ４０：  ４０は、検出されなかった遺伝子を呈示するのに十分なほど高いとみなされた任意Ｃｔ値である。しかし、それを４０ではなく４２に設定すれば、結果はすべて変化するであろう。したがって、本発明者らはまず４０ではないすべての数値を調べてそれらをランク付けすることによりこれを解決した。Ｈａｓａｎ遺伝子について、これは４６２３の数値をランク付けすることに相当する。次いでこれらの遺伝子のうち底２％、すなわち最底の９２遺伝子を調べた；それらの平均および標準偏差を計算すると、３７．９および０．８となった。次いで各４０を区間［３７．９−０．８、３７．９＋０．８］間でランダムに選択した数字に置き換えた。
【０１７５】
[00194]  異常値：発現レベルを手動で調べる場合、与えられた遺伝子について異常値として行動する試料がしばしばみられる。これを克服するために、ＦＣを計算する際にグループ（ＮまたはＦ）において最高および最低の発現レベルを除外した。同様に、各グループにおいて最高２つおよび最低２つを除外することによりこの操作を繰り返した。
【０１７６】
[00195]  結論として、以上の統計的方法を用いることにより、本発明者らは、卵丘細胞による表４中の遺伝子（この遺伝子セットはＡＢＣＡ６、ＤＤＩＴ４、ＤＵＳＰ１、ＧＰＲ１７３Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３およびＴＥＲＦ２ＩＰ、ならびにさらに６つの遺伝子ＫＲＡＳ、ＮＣＡＭｌ、ＯＬＦＭＬ３、ＰＴＰＲＡおよびＳＤＦ４を含む）の発現レベルは、それに随伴する卵母細胞または同じ女性ドナーからの卵母細胞が生存可能な妊娠を生じる能力と相関することを見出した。したがって、卵丘細胞によるこれらの遺伝子のうち１以上の発現を検出する方法は、それに随伴する卵母細胞または同じ女性ドナーからの卵母細胞がＩＶＦ法に使用するのに適切であるかどうかを判定するために使用できる。
【０１７７】
[00196]  これらの結果を確認するために、以下に記載するように１４Ｎ対１７Ｆの比較および１２Ｎと１５Ｆの比較を行なった。
[00197]  遺伝子のフィルタリング：大部分の試料中に検出されない若干の遺伝子がある。極端な場合、３５試料すべてにおいて４０のＣｔ値を与えられる“Ｈａｓａｎ遺伝子”と表示される２０の遺伝子がある。これらの遺伝子を採用すると予測値にはならない。一般的方法として、本発明者らは２５以上の試料中に検出されない遺伝子を除外した。
【０１７８】
[00198]  結果
[00199]  以上の方法を用いて、本発明者らは異常値の除外に基づいて３つのデータセットを作成した；ここで異常値は、ある遺伝子について、与えられたあるグループにおいて最高または最低の発現を伴う試料（単数または複数）である：ｉ）試料を除外しなかった（Ｅｘｃｅｌファイル名に“除外なし”により表示）、ｉｉ）各グループ（ＮおよびＦ）において最高および最低の発現値を伴う試料を除外し、すなわち１４Ｎおよび１７Ｆの試料を用いてＦＣを計算する（Ｅｘｃｅｌファイル名に“１つを除外”により表示）、ｉｉｉ）各グループ（ＮおよびＦ）において最高のトップ２つおよび最低のトップ２つの発現値を除外し、すなわち１２Ｎおよび１５Ｆの試料を用いてＦＣを計算する（Ｅｘｃｅｌファイル名に“２つを除外”により表示）。
【０１７９】
[00200]  各データセットについて、４０を前記に説明したように扱った；すなわち、それらを、検出したすべての値の最低２％の平均値±それらの数値の標準偏差のうちのランダムな数値で置き換えた。次いで発現レベルを対数解除し、各データセットについて２つのファイルを作成した：スケーリングしたもの（与えられた試料について、各遺伝子の発現レベルをその試料中のＧＡＰＤＨの発現レベルで割る）およびスケーリングしていないもの（ＧＡＰＤＨまたは他の何らかのスケーリングを適用しない）。
【０１８０】
[00201]  各Ｅｘｃｅｌファイル（合計６；３データセットそれぞれについて２つ）につき、下記の欄も含めた：“グループ”：Ｐ（遺伝子は予測シグナチャーからのもの）  Ａ（遺伝子は差示発現分析からのもの）  トレイニング試料および検証試料、すなわちマイクロアレイ分析のすべての試料を使用；“カウント４０”：遺伝子が数値４０と推定される試料の数を示す；“ＦＣ  ＴＬＤＡ”：ＦＣを平均の比として計算；“ＦＣ  Ａｆｆｙ”：ＦＣはマイクロアレイから得られる；“一致”：増加の方向がＴＬＤＡとＡｆｆｙで同じであるかどうかを示す；一致の指標として１０、そうでなければ−１０を用いた。
【０１８１】
[00202]  図９に、ＴＬＤＡとマイクロアレイの間で一致する遺伝子の分布を示す。これらの結果を表にしたものが表８にあり、その際、一致レベルをマイクロアレイ結果からの遺伝子のアップ／ダウンレギュレーションに基づいて分ける。ここにはすべての遺伝子を用いる；すなわち遺伝子が４０であることが示された試料の数に基づくフィルタリングを適用していない。
【０１８２】
[00203]  下記の表８に、アップレギュレーションの方向がＴＬＤＡとマイクロアレイ分析の間で一致を示す遺伝子の数を示す；ＴＬＤＡとマイクロアレイ分析の間のアップレギュレーションの方向。Ｓ：スケーリングしたもの，Ｕ：スケーリングしていないもの．Ｐ：予測遺伝子リスト．Ａ：全試料を一緒に分析した際に得られたマイクロアレイデータセット中の遺伝子．Ｔ：総遺伝子．遺伝子をフィルタリングしなかった。
【０１８３】
[00204]  表８
【０１８４】
【表８】
2013510575.tif 000026
【０１８５】
[00205]  これらの結果は、各グループにおいてトップの１または２つの異常値を除外すると、スケーリングしたデータを用いて最良の全体的一致が６１％で達成されることを示唆する。さらに、２５以上の試料において“検出されない”遺伝子をフィルタリングで除去することによりこの分析を繰り返した。１１以上の試料において検出される遺伝子を用いてＡｆｆｙとＴＬＤＡの間の一致の分布を示す結果を表９にまとめる。
【０１８６】
【表９】
2013510575.tif 000027
【０１８７】
[00206]  両グループにおいていずれの側にもトップの１または２つの異常値を除外し、かつスケーリングした数値を用いた場合は、全体として６６％が一致するので、これらの結果はさらに最適化される。同様に興味深いのは、スケーリングしていない数値では、Ｆグループにおいてほぼすべての遺伝子がアップレギュレーションを示すことである。
【０１８８】
[00207]  予測
[00208]  予測分析には、２５以上の試料において“検出されない”遺伝子をフィルタリングで除去した場合、ＡｆｆｙとＴＬＤＡの間で一致する遺伝子（表９）を用いた。本発明者らは、重み付き投票法を用いる１点排除交差検証法（leave one out cross validation）を適用した。スケーリングしていないデータを用いた場合、６遺伝子で最良の結果が得られた（２７／３５＝７７％の予測精度）。スケーリングしたデータの場合、最良の予測精度は（２２／３５＝６３％）であった。これらの遺伝子を、３５試料を用いたＴＬＤＡにおけるそれらの変化倍率値（Ｆ／Ｎ）およびＡｆｆｙ、ならびに“新しい”１４試料におけるこれらの遺伝子の変化倍率値（Ｆ／Ｎ）と共に表１０に示す。
【０１８９】
【表１０】
2013510575.tif 000028
【０１９０】
[00209]  ３５試料についての予測結果を表１１に示す。不適正に予測された試料に影を付けてある。合計で３つのＮおよび５つのＦ試料が不適正に予測された。
【０１９１】
【表１１】
2013510575.tif 000029
【０１９２】
【表１２】
2013510575.tif 000030
【０１９３】
[00210]  結論として、以上の統計的方法を用いて、本発明者らは、卵丘細胞による表４中の１以上の遺伝子、より好ましくはＡＢＣＡ６、ＮＣＡＭ１、ＯＬＦＭＬ３、ＰＴＰＲＡ、ＳＤＦ４、ＧＰＲ１３７Ｂ、ＤＤＩＴ４、ＤＵＳＰｌ、ＧＰＲ１３７Ｂ、ＩＤＵＡ、ＫＣＴＤ５、ＮＤＮＬ２、ＳＬＣ２６Ａ３、およびＴＥＲＦ２ＩＰからなる群から選択される１４の遺伝子のうち１以上の発現レベルは、それに随伴する卵母細胞または同じ女性ドナーからの卵母細胞が生存可能な妊娠を生じる能力と相関することを見出した。したがって、卵丘細胞によるこれらの遺伝子のうち１以上の発現を検出する方法は、それに随伴する卵母細胞または同じ女性ドナーからの卵母細胞がＩＶＦ法に使用するのに適切であるかどうかを判定するために、またＩＶＦ法に使用するのに不適切な卵母細胞を生じる状態を伴う個体を同定するために、また受精処置の成功をモニターするために使用できる。
【０１９４】
[00211]  参考文献
[00212]  本明細書全体を通して、種々の参考文献が本発明に関係する技術水準について記載している。これらの参考文献の開示内容を本発明の開示に引用することにより本明細書に援用する。
【図１ａ】
2013510575.tif 000031
【図１ｂ】
2013510575.tif 000032
【図１ｃ】
2013510575.tif 000033
【図２】
2013510575.tif 000034
【図３】
2013510575.tif 000035
【図４】
2013510575.tif 000036
【図５】
2013510575.tif 000037
【図６】
2013510575.tif 000038
【図７】
2013510575.tif 000039
【図８】
2013510575.tif 000040
【図９】
2013510575.tif 000041
【国際調査報告】
2013510575.tif 000042
2013510575.tif 000043
2013510575.tif 000044
2013510575.tif 000045