IP Force 特許公報掲載プロジェクト 2022.1.31 β版

知財求人 - 知財ポータルサイト「IP Force」

▶ 株式会社ジェノミックスの特許一覧 ▶ 国立大学法人大阪大学の特許一覧

<>
  • 特許-三次元DNA構造相互作用分析方法 図1
  • 特許-三次元DNA構造相互作用分析方法 図2
  • 特許-三次元DNA構造相互作用分析方法 図3
  • 特許-三次元DNA構造相互作用分析方法 図4
  • 特許-三次元DNA構造相互作用分析方法 図5
  • 特許-三次元DNA構造相互作用分析方法 図6
  • 特許-三次元DNA構造相互作用分析方法 図7
  • 特許-三次元DNA構造相互作用分析方法 図8
  • 特許-三次元DNA構造相互作用分析方法 図9
  • 特許-三次元DNA構造相互作用分析方法 図10
  • 特許-三次元DNA構造相互作用分析方法 図11
  • 特許-三次元DNA構造相互作用分析方法 図12
< >
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-11-25
(45)【発行日】2024-12-03
(54)【発明の名称】三次元DNA構造相互作用分析方法
(51)【国際特許分類】
   C12Q 1/6869 20180101AFI20241126BHJP
   C12Q 1/686 20180101ALI20241126BHJP
   C12Q 1/34 20060101ALI20241126BHJP
【FI】
C12Q1/6869 Z
C12Q1/686 Z
C12Q1/34
【請求項の数】 12
(21)【出願番号】P 2020565236
(86)(22)【出願日】2020-01-10
(86)【国際出願番号】 JP2020000752
(87)【国際公開番号】W WO2020145405
(87)【国際公開日】2020-07-16
【審査請求日】2023-01-06
(31)【優先権主張番号】P 2019003811
(32)【優先日】2019-01-11
(33)【優先権主張国・地域又は機関】JP
(73)【特許権者】
【識別番号】506364400
【氏名又は名称】株式会社ステムリム
(73)【特許権者】
【識別番号】504176911
【氏名又は名称】国立大学法人大阪大学
(74)【代理人】
【識別番号】100145403
【弁理士】
【氏名又は名称】山尾 憲人
(74)【代理人】
【識別番号】100122301
【弁理士】
【氏名又は名称】冨田 憲史
(74)【代理人】
【識別番号】100170520
【弁理士】
【氏名又は名称】笹倉 真奈美
(72)【発明者】
【氏名】玉井 克人
(72)【発明者】
【氏名】新保 敬史
(72)【発明者】
【氏名】山崎 尊彦
【審査官】小倉 梢
(56)【参考文献】
【文献】eLife,2017年,Vol. 6:e21856,p. 1-35
【文献】bioRxiv,doi: https://doi.org/10.1101/193219,2017年,p. 1-28
【文献】Nature,2009年,Vol. 462,p. 58-64
【文献】Mol. Cytogenet.,2018年,Vol. 11:21,p. 1-10
(58)【調査した分野】(Int.Cl.,DB名)
C12Q 1/00 - 1/70
JSTPlus/JMEDPlus/JST7580(JDreamIII)
CAplus/MEDLINE/EMBASE/BIOSIS(STN)
(57)【特許請求の範囲】
【請求項1】
三次元DNA構造におけるDNA配列の相互作用を分析する方法であって、
(1)DNA切断ステップの前に、細胞核内におけるDNA配列の空間配置を固定するための架橋ステップと、
(2)前記DNA配列中のDNA鎖に接する任意のタンパク質を選択して、前記任意のタンパク質に一次抗体を結合させる一次抗体結合ステップと、前記一次抗体を介して前記任意のタンパク質にDNAを切断する酵素を結合させる酵素結合ステップとを含む、前記DNAを切断する酵素によって、前記任意のタンパク質の近傍のDNA鎖を特異的に切断することを特徴とするDNA切断ステップと、
(3)前記DNA切断ステップによる切断で生じたDNA末端をライゲーションするライゲーションステップと、
(4)前記ライゲーションされたDNA断片を選別する選別ステップと、
(5)前記選別ステップにて選別されたDNA断片に対し塩基配列を決定する配列決定ステップと、
(6)前記配列決定ステップの前にDNA断片に対しDNA構造の固定を解消する固定解消ステップとを含み、
前記DNA配列の相互作用を分析することを特徴とする三次元DNA構造相互作用分析方法。
【請求項2】
前記DNA配列の相互作用が、前記任意のタンパク質の近傍のDNA配列と、前記近傍のDNA配列に対してゲノム上の座位は離れているが空間的に近接するDNA配列との間の相互作用であることを特徴とする、請求項1記載の三次元DNA構造相互作用分析方法。
【請求項3】
前記架橋ステップは、ホルムアルデヒド、EGS(ethylene glycol bis(succinimidyl succinate))およびDSG(disuccinimidyl glutarate)からなる群より選択される薬剤により固定されることを特徴とする、請求項1または2に記載の三次元DNA構造相互作用分析方法。
【請求項4】
前記DNA切断ステップによる切断によって生じたDNA末端に標識を付加する標識付加ステップを含んでおり、
前記選別ステップは、ライゲーションされたDNA断片を前記標識に基づいて選別することを特徴とする、請求項1から3のいずれかに記載の三次元DNA構造相互作用分析方法。
【請求項5】
前記選別ステップは、前記ライゲーションステップにてライゲーションされたDNA断片を、より短いDNA断片に断片化した後に、選別を行うことを特徴とする、請求項1から4のいずれかに記載の三次元DNA構造相互作用分析方法。
【請求項6】
前記配列決定ステップは、選別されたDNA断片に対し、増幅を行うと共に増幅産物の塩基配列を決定する配列決定を行うことを特徴とする、請求項1から5のいずれかに記載の三次元DNA構造相互作用分析方法。
【請求項7】
前記選別ステップは、増幅されるDNA断片に対し、所定のアダプター配列を切断末端に付加することを特徴とする、請求項6に記載の三次元DNA構造相互作用分析方法。
【請求項8】
三次元DNA構造におけるDNA配列の相互作用を分析する方法であって、
(1)DNA切断ステップの前に、細胞核内におけるDNA配列の空間配置を固定するための架橋ステップと、
(2)前記DNA配列中のDNA鎖に接する任意のタンパク質を選択して、前記任意のタンパク質に一次抗体を結合させる一次抗体結合ステップと、前記一次抗体を介して前記任意のタンパク質にDNAを切断する酵素を結合させる酵素結合ステップとを含む、前記DNAを切断する酵素によって、前記任意のタンパク質の近傍のDNA鎖を特異的に切断することを特徴とするDNA切断ステップと、
(3)前記DNA切断ステップによる切断で生じたDNA末端をライゲーションするライゲーションステップと、
(4)前記ライゲーションされたDNA断片を捕捉する捕捉ステップと、
(5)前記捕捉ステップにより捕捉されたDNA断片に対しDNA構造の固定を解消する固定解消ステップと、
(6)固定が解消された後のDNA断片に対し、塩基配列を決定する配列決定ステップとを含み、
前記DNA配列の相互作用を分析することを特徴とする三次元DNA構造相互作用分析方法。
【請求項9】
三次元DNA構造におけるDNA配列の相互作用を分析する方法であって、
(1)DNA切断ステップの前に、細胞核内におけるDNA配列の空間配置を固定するための架橋ステップと、
(2)前記DNA配列中のDNA鎖に接する任意のタンパク質を選択して、前記任意のタンパク質に一次抗体を結合させる一次抗体結合ステップと、前記一次抗体を介して前記任意のタンパク質にDNAを切断する酵素を結合させる酵素結合ステップとを含む、前記DNAを切断する酵素によって、前記任意のタンパク質の近傍のDNA鎖を特異的に切断することを特徴とするDNA切断ステップと、
(3)前記DNA切断ステップによる切断で生じたDNA断片を捕捉する捕捉ステップと、
(4)前記捕捉ステップにより捕捉されたDNA断片の末端をライゲーションするライゲーションステップと、
(5)前記ライゲーションされたDNA断片に対しDNA構造の固定を解消する固定解消ステップと、
(6)固定が解消された後のDNA断片に対し塩基配列を決定する配列決定ステップとを含み、
前記DNA配列の相互作用を分析することを特徴とする三次元DNA構造相互作用分析方法。
【請求項10】
前記DNA配列の相互作用が、前記任意のタンパク質の近傍のDNA配列と、前記近傍のDNA配列に対してゲノム上の座位は離れているが空間的に近接するDNA配列との間の相互作用であることを特徴とする、請求項8または9のいずれかに記載の三次元DNA構造相互作用分析方法。
【請求項11】
前記架橋ステップは、ホルムアルデヒド、EGS(ethylene glycol bis(succinimidyl succinate))およびDSG(disuccinimidyl glutarate)からなる群より選択される薬剤により固定されることを特徴とする、請求項8から10のいずれかに記載の三次元DNA構造相互作用分析方法。
【請求項12】
前記配列決定ステップは、固定が解消された後のDNA断片に対し、増幅を行うと共に増幅産物の塩基配列を決定する配列決定を行うことを特徴とする、請求項8から11のいずれかに記載の三次元DNA構造相互作用分析方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、DNA配列の三次元構造における相互作用を分析する方法およびこれを実現する分析用キットに関し、特に、目的DNA配列に対して相互作用するDNA配列を迅速かつ的確に分析することのできる三次元DNA構造相互作用分析方法に関する。
【背景技術】
【0002】
近年、ヒトゲノムの配列が解読され、ヒトの遺伝情報が明らかとなったが、この情報は単に一次元の直線的な配列情報としてもたらされている。しかし、実際のゲノムDNAは、折り畳まれた複雑な三次元構造を形成して核内に収容されており、配列各部の空間的な配置(遠近)を一次元の直線配列から完全に予測することは困難である。つまり、DNA配列のある領域に着目すると、その領域に対し空間的に近傍に位置しているDNA配列は、一次元直線配列としては遠隔に位置する配列であるといったことが往々に生じている。
【0003】
また、近年、特定の遺伝子の発現は、三次元空間における配列間の相互作用にて制御されることが明らかとなりつつあり、配列の相互作用の解析法として染色体高次構造捕捉法(chromosome conformation capture;3C法)やこれが改良されたHiC法などが提案されている。
【0004】
3C法は、ゲノムの高次構造を維持したまま固定(架橋)してDNA切断をした後、空間的に近接する末端をライゲーション(近接ライゲーション)によって連結させる手法である。高次構造が固定された状態での近接ライゲーションは、空間的に近傍に位置する配列同士の結合となる。このため、連結によって生じるDNA配列は、ゲノム上の一次元の配列とは異なる配列が確認されることとなる。これまでに提案されているDNA配列の三次元の相互作用の分析方法はこの3C法がベースとなっている。
【0005】
しかしながら、3C法では、予め着目した配列間、即ち、既に明らかとなっている配列間の相互作用しか解析ができない。
【0006】
HiC法はその欠点を補うものであり、ゲノムワイドにDNA相互作用を解析する手法として提案され、未知の配列間の相互作用を分析できるよう構成されている。HiC法では、DNA切断により生じるDNA断片の末端をビオチン化した後に近接ライゲーションが行われる。ライゲーションされたDNA断片の一群はビオチン標識をマーカーとして回収され、更なる断片化を経て、PCR用のアダプターが付加され、増幅、シーケンスによる塩基配列の同定が実施される手法が基本となっている。
【0007】
一方で、HiC法では、理論上すべてのDNA三次元構造を検出するため、シーケンスコストが膨大になるうえ、非常に煩雑なコンピューター解析が必要となるため汎用性に問題を残している。そこで、低コスト化を目的に、特定配列とのハイブリダイゼーションを組み合わせた新たな手法が提案されている(特許文献1および非特許文献1参照)。
【0008】
また、DNAはクロマチンを形成しているので、クロマチンを断片化した後、クロマチン免疫沈降(ChIP)を利用して分析対象のクロマチン断片を回収し、クロマチン断片に基づいてDNA配列の相互作用を分析するChIA-PET法(特許文献2参照)が提案されている。
【先行技術文献】
【特許文献】
【0009】
【文献】特表2016-537029
【文献】特開2007-289152
【非特許文献】
【0010】
【文献】Stefan Schoenfelder, Biola-Maria Javierre, Mayra Furlan-Magaril, Steven W. Wingett, and Peter Fraser,Promoter Capture Hi-C: High-resolution, Genome-wide Profiling of Promoter Interactions,Journal of Visualized Experiments.(136), e57320,June 2018
【発明の概要】
【発明が解決しようとする課題】
【0011】
しかしながら、上述のHiC法は、実際上少なくとも10個以上の細胞が必要なサンプル数とされており、1細胞あたり10種類以上の断片のライゲーションが生じるとされ、作られるDNAペア(空間的に近接するDNA配列がライゲーションによって横並びとなったDNA断片)の自由度は、一説には、1018程度もの桁数に上るものとも推定されている。ここで、DNAペアの塩基配列の決定に先立って行われるPCRによる増幅では、GC%や制限酵素断片の長さの違いなど種々の要因から、増幅によって生成するDNAペアにバイヤスがかかることが知られている。また、得られたDNAペアの集合体(HiCライブラリ)のシーケンスは、通常のシーケンスに比べて得られるデータの複雑性が圧倒的に大きいため、解析結果に対し解釈の入り込む余地が大きく、分析結果の信頼性が、場合によっては低下しかねないという問題点を抱えている。
【0012】
特許文献1、非特許文献1に記載される手法は、従来の網羅的なHiC法を改良し、HiCライブラリに対し、DNAマイクロアレイでハイブリダイズを行うことで、注目領域特異的なシーケンスの実施を提案している。これによって、シーケンス負荷の軽減と、データ量の圧縮による解析負担軽減が図られている。しかし、DNAペアを作製するためのクロマチンの切断処理は、従来と同様にゲノムワイドに実施されている。従って、全体として膨大量のDNAペアが生じており、目的物以外に膨大な夾雑物が共存している状態となっている。このため、ミスハイブリダイズを回避できない。加えて、マイクロアレイに用いられるプローブには、配列設計段階で意図せずバイヤスがかかってしまうことが知られており、最終的に得られるデータが、十分に信頼性のあるものとならないという問題点があった。
【0013】
更に、クロマチン免疫沈降法は、制限酵素によるDNA切断にてクロマチンを断片化し、ヒストン等のタンパク質をターゲットに抗体を用いて免疫沈降を行うものである。免疫沈降に際しては、超音波処理や制限酵素処理によって断片化したクロマチン断片を液相へ可溶化させ、抗体免疫沈降を行うことで、DNA断片が採取され、分析用のライブラリソースが形成される。ところが、この可溶化工程(核画分の可溶化工程とも称される)において、不溶性画分がしばしば発生する。可溶化されないクロマチン断片は採取不能なため、ライブラリソースに取り込まれず、分析対象から除外されてしまう。結果として得られる解析データは実際とはズレたものになってしまいかねず、得られる結果の信頼性を低下させてしまっているという問題点があった。
【0014】
本発明は、上述した問題点に鑑みてなされたものであり、簡便かつ低コストで、解析データの信頼性に優れ、高精度にDNA配列相互作用を解析する手法およびこれを実現する分析用キットの提供を目的としている。
【課題を解決するための手段】
【0015】
[項目1]
三次元DNA構造におけるDNA配列の相互作用を分析する方法であって、
前記DNA配列中のDNA鎖に接する任意のタンパク質を選択して、前記任意のタンパク質の近傍のDNA鎖を特異的に切断するDNA切断ステップと、
前記DNA切断ステップによる切断で生じたDNA末端をライゲーションするライゲーションステップと、
前記ライゲーションされたDNA断片を選別する選別ステップと、
前記選別ステップにて選別されたDNA断片に対し塩基配列を決定する配列決定ステップと、
前期配列決定ステップの前にDNA断片に対しDNA構造の固定を解消する固定解消ステップとを含み、
前記DNA配列の相互作用を分析することを特徴とする三次元DNA構造相互作用分析方法。
【0016】
[項目2]
前記DNA切断ステップによる切断によって生じたDNA末端に標識を付加する標識付加ステップを含んでおり、
前記選別ステップは、ライゲーションされたDNA断片を前記標識に基づいて選別することを特徴とする、項目1記載の三次元DNA構造相互作用分析方法。
【0017】
[項目3]
前記選別ステップは、前記ライゲーションステップにてライゲーションされたDNA断片を、より短いDNA断片に断片化した後に、選別を行うことを特徴とする、項目1または2に記載の三次元DNA構造相互作用分析方法。
【0018】
[項目4]
前記配列決定ステップは、選別されたDNA断片に対し、増幅を行うと共に増幅産物の塩基配列を決定する配列決定を行うことを特徴とする、項目1から3のいずれかに記載の三次元DNA構造相互作用分析方法。
【0019】
[項目5]
前記選別ステップは、増幅されるDNA断片に対し、所定のアダプター配列を切断末端に付加することを特徴とする、項目4に記載の三次元DNA構造相互作用分析方法。
【0020】
[項目6]
三次元DNA構造におけるDNA配列の相互作用を分析する方法であって、
前記DNA配列中のDNA鎖に接する任意のタンパク質を選択して、前記任意のタンパク質の近傍のDNA鎖を特異的に切断するDNA切断ステップと、
前記DNA切断ステップによる切断で生じたDNA末端をライゲーションするライゲーションステップと、
前記ライゲーションされたDNA断片を捕捉する捕捉ステップと、
前記捕捉ステップにより捕捉されたDNA断片に対しDNA構造の固定を解消する固定解消ステップと、
固定が解消された後のDNA断片に対し、塩基配列を決定する配列決定ステップとを含み、
前記DNA配列の相互作用を分析することを特徴とする三次元DNA構造相互作用分析方法。
【0021】
[項目7]
三次元DNA構造におけるDNA配列の相互作用を分析する方法であって、
前記DNA配列中のDNA鎖に接する任意のタンパク質を選択して、前記任意のタンパク質の近傍のDNA鎖を特異的に切断するDNA切断ステップと、
前記DNA切断ステップによる切断で生じたDNA断片を捕捉する捕捉ステップと、
前記捕捉ステップにより捕捉されたDNA断片の末端をライゲーションするライゲーションステップと、
前記ライゲーションされたDNA断片に対しDNA構造の固定を解消する固定解消ステップと、
固定が解消された後のDNA断片に対し塩基配列を決定する配列決定ステップとを含み、
前記DNA配列の相互作用を分析することを特徴とする三次元DNA構造相互作用分析方法。
【0022】
[項目8]
前記配列決定ステップは、固定が解消された後のDNA断片に対し、増幅を行うと共に増幅産物の塩基配列を決定する配列決定を行うことを特徴とする、項目6または7に記載の三次元DNA構造相互作用分析方法。
【0023】
[項目9]
前記DNA切断ステップは、前記任意のタンパク質に一次抗体を結合させる一次抗体結合ステップと、前記一次抗体を介して前記任意のタンパク質にDNAを切断する酵素を結合させる酵素結合ステップとを含んでおり、
前記DNAを切断する酵素によって、前記任意のタンパク質の近傍のDNA鎖を特異的に切断することを特徴とする、項目1から8のいずれかに記載の三次元DNA構造相互作用分析方法。
【0024】
[項目10]
前記DNA切断ステップの前に、DNA構造を架橋によって固定化する架橋ステップを含む、項目1から9のいずれかに記載の三次元DNA構造相互作用分析方法。
【0025】
[項目11]
三次元DNA構造におけるDNA配列の相互作用を分析するために用いられる分析用キットであって、
任意のタンパク質に結合する一次抗体と、
前記一次抗体に結合するDNA切断酵素と、
前記DNA切断酵素による切断によって生じるDNA末端を識別するための標識と、
標識されるDNA末端をライゲーションする酵素とを含むことを特徴とする、三次元DNA構造相互作用分析キット。
【発明の効果】
【0026】
本発明の方法によれば、解析対象のDNA領域に対し選択的で局所的なDNA切断(DNA断片化)を行って解析サンプルが作製されるので、解析対象となるDNA領域にフォーカスされ夾雑物の少ない高品質な解析サンプルが得られる。ゆえに、準備するサンプル量を低減することができ、場合によっては圧倒的に少なくすることができる。
また、ゲノムワイドな解析に基づいて相互作用を特定するHiC法と比べた場合には、本発明の方法では、目的のDNA領域に絞って配列間の相互作用の解析が可能となるので、必要な解析結果を得るまでの時間を大幅に短縮し得、迅速かつ低コストでの解析が可能となる。
また、増幅を行う場合であっても、従来法のように多量の夾雑物が混在するサンプルでは、PCR増幅工程において、大量に存在する他の断片も増幅されるため、解析対象のいずれの断片についてもシーケンスデータを得るには、仕込みサンプル量を大量に準備せねばならないが、本発明では、夾雑物が少ないことから、PCR増幅を行ってもノイズの増大が抑制できるため、信頼性の高い解析結果を得ることができる。
【0027】
一実施形態において、本発明の方法を、DNA切断ステップが切断によって生じたDNA末端を選別するため標識を付加する標識付加ステップを含み、選別ステップにより、標識に基づいてライゲーションされたDNA断片を選別するようにした場合には、タンパク質に対する免疫沈降(クロマチン免疫沈降)によらず、DNA断片回収を成し得るので、クロマチン免疫沈降法で必要となる核画分の可溶化工程を不要とでき、不溶性画分の影響を排除できる。すなわち、全ゲノムが解析対象とされるHiChIP法等は従来のクロマチン免疫沈降を利用する解析手法であるが、核分画の可溶化工程において生じる不溶性画分の影響から、解析不能となる領域の発生を完全に回避することが困難である。このため、得られたシーケンス結果から相互作用を解析すると、場合によっては解析結果の信頼性が低下してしまう。しかし、本発明を、標識付加ステップを設けて構成した場合には、クロマチン免疫沈降法におけるような不溶性画分の影響を排除できるため、より、一層信頼性の高い解析結果を得ることができる。
【図面の簡単な説明】
【0028】
図1図1は、一例として、本発明の方法のフローチャートを示す。
図2図2は、H3K4me3が存在するクロマチン領域の相互作用のヒートマップを示す。縦・横軸はそれぞれ第1染色体から第23染色体を示し、遺伝子領域が並行している。図2中においては縦横の直線によって各染色体に対応する領域を区分しており、同染色体の同遺伝子領域は対角線で示す。検出された異なる遺伝子領域の相互作用は、点で表され、クロマチン領域の相互作用の度合いは、ヒートマップ上の点の密度に基づくグレースケールの濃淡で示されている。
図3図3は、一例として解析された、ゲノム上におけるH3K4me3の局在およびH3K4me3が存在するクロマチン領域間の相互作用を示す。最上部にはゲノム上の位置が数字で示されており、上列は、既知のCUT&RUN法(Skene et al., 2018)を実施して得た結果において解析されたH3K4me3の局在を示した分布図である。中列および下列は、本願実施例1により得られた結果であって、中列はH3K4me3の局在を示した分布図であり、下列にはH3K4me3が存在するクロマチン領域間の相互作用をループ状で示す。下列の下方には、第7染色体中の遺伝子(FBXL18、MIR589、ACTB)の位置が表示されており、遺伝子名称の右側にマーキングされた範囲が対応する遺伝子の存在位置を示している。
図4図4は、図3の部分的な拡大図である。
図5図5は、図3の部分的な拡大図である。
図6図6は、CTCFが存在するクロマチン領域の相互作用のヒートマップを示す。縦・横軸はそれぞれ第1染色体から第23染色体を示し、遺伝子領域が並行している。図6中においては縦横の直線によって各染色体に対応する領域を区分しており、同染色体の同遺伝子領域は対角線で示す。検出された異なる遺伝子領域の相互作用は、点で表され、クロマチン領域の相互作用の度合いは、ヒートマップ上のグレースケールの濃淡で示されている。
図7図7は、一例として解析された、ゲノム上におけるCTCFの局在およびCTCFが存在するクロマチン領域間の相互作用を示す。最上部にはゲノム上の位置が数字で示されており、上列は、本願実施例2により得られたCTCFの局在を示した分布図である。中列は、Jaspar Transcription Factor Databaseにより特定されたゲノム上に存在するCTCF結合サイトを示しており、ゲノム上の配列の方向性に応じ、一方の方向のものに下線を付して示している。下列は、本願実施例2により得られたCTCFが存在するクロマチン領域間の相互作用をループ状で示す。下列の下方には、第7染色体中の遺伝子(FBXL18、MIR589、ACTB)の位置が表示されており、遺伝子名称の右側にマーキングされた範囲が対応する遺伝子の存在位置を示している。
図8図8は、図7の部分的な拡大図である。
図9図9は、図7の部分的な拡大図である。
図10図10は、一例として解析された、ゲノム上におけるCTCFの局在およびCTCFが存在するクロマチン領域間の相互作用を示す。最上部にはゲノム上の位置が数字で示されており、上列は、既知のCUT&RUN法(Skene et al., 2018)を実施して得た結果において解析されたCTCFの局在を示した分布図である。中列は、Jaspar Transcription Factor Databaseにより特定されたゲノム上に存在するCTCF結合サイトを示しており、ゲノム上の配列の方向性に応じ、一方の方向のものに下線を付して示している。下列は、本願実施例3により得られたCTCFの局在を示した分布図と共に、CTCFの局在箇所の下に、本願実施例3により得られたCTCFが存在するクロマチン領域間の相互作用をループ状で示す。下列の更に下方には、第7染色体中の遺伝子(FBXL18、MIR589、ACTB)の位置が表示されており、遺伝子名称の右側にマーキングされた範囲が対応する遺伝子の存在位置を示している。
図11図11は、図10の部分的な拡大図である。
図12図12は、図10の部分的な拡大図である。
【発明を実施するための形態】
【0029】
本発明は、DNAの三次元構造に基づくDNA配列の相互作用を分析する方法およびこれを実現する分析用キットである。
【0030】
本発明において、DNAの三次元構造とは、DNA二重らせんが形成する高次構造、例えばループ及び折り畳みを有する構造を意味する。これらの構造体は、DNAのみから成ってもよく、又はタンパク質などの他の分子を含んでもよい。DNAとタンパク質複合体のクロマチンも含まれる概念である。クロマチンの主要な機能は、細胞内に収まるように、DNAを小さな体積にパッケージすることに加え、遺伝子発現、DNA複製及び修復を制御することといわれている。クロマチンの最も豊富なタンパク質成分は、DNAを凝縮させるヒストンである。
【0031】
近年、トポロジカルドメイン(以下、適宜「TD」と称す)と称される、大きなメガベースサイズの局所的クロマチン相互作用ドメインが同定されている。これらのドメインは、ヘテロクロマチンの広がりを制限するゲノムの領域と相関する。このドメインは、異なる細胞型にわたって安定であり、種を超えて高度に保存されている。トポロジカルドメインはまた、互いに相互作用する。本発明においては、このトポロジカルドメイン内のみならず、トポロジカルドメイン間のDNA配列相互作用についても分析対象としている。
【0032】
このような三次元構造に起因して、ゲノムDNAの配列各部の空間的な配置(遠近)を一次元の直線配列から完全に予測することは困難となっている。つまり、DNA配列のある領域に着目すると、その領域に対し空間的に近傍に位置しているDNA配列は、一次元直線配列としては遠隔に位置する配列であるといったことが往々に生じている。遺伝子発現などは、DNA配列の相互作用によって調整されると推定されており、DNA配列の空間的な遠近は重要な情報になっている。DNA配列の相互作用または三次元DNA構造相互作用は、例えば、三次元構造におけるDNA配列間の空間的な配置(例えば遠近)を表す。
【0033】
本発明の方法は、DNA鎖に結合、複合化または付着する各種タンパク質のうち、予め選択した任意のタンパク質の近傍のDNA鎖(目的のDNA領域)を特異的に切断する切断ステップと、その切断で生じたDNA末端をライゲーションすることで、(ゲノム上の座位は離れているが)空間的に近接するDNA配列が横並びとなった配列(DNAペア)を作製するライゲーションステップと、作製されたDNAペアを選別する選別ステップまたはライゲーションされたDNA断片もしくはDNA切断ステップによる切断で生じたDNA断片を(結合するタンパク質分子に基づいて)捕捉する捕捉ステップの少なくともいずれか一方と、固定解消ステップと、DNAペアを所望により増幅し、塩基配列を決定する配列決定ステップとを、任意の順序で、含むことができる。尚、目的のDNA領域とは、任意のタンパク質の近傍のDNA鎖の部分であって、解析対象となるDNA領域である。以下、各ステップについて説明する。
【0034】
なお、本発明は、切断ステップの前に、細胞内の高次構造を人為的操作によって固定する架橋ステップを設けて構成してもよい。
【0035】
(架橋ステップ)
架橋ステップは、細胞核内におけるDNA配列の空間配置を固定するための処置である。例えば、ホルムアルデヒドは、タンパク質と核酸、あるいは1のタンパク質と他の近隣のタンパク質とを架橋するために使用し得る。プロモータとエンハンサーとの関係等、相互作用する塩基配列は、多くの場合、特定のタンパク質と接しており、タンパク質を介して近接している。言い換えれば、DNAの高次構造は、多くの場合、タンパク質が関わって形成されている。このため、ホルムアルデヒドを用いてDNAとこれに接するタンパク質とを固定すること(架橋)でDNA配列の空間構造が固定される。
【0036】
本発明において使用し得るDNAを架橋する薬剤の例としては、ホルムアルデヒドに加え、例えば、EGS(ethylene glycol bis(succinimidyl succinate))、DSG(disuccinimidyl glutarate)が挙げられる。
【0037】
ホルムアルデヒドを用いれば、比較的短い距離、例えば約2Åを橋渡しする架橋を形成し、それによって、DNA配列間の密接な相互作用について解析するためのDNAペアを形成することができることとなる。
【0038】
ホルムアルデヒド架橋は、直接相互作用し合う二つの配列の検出に好適な手法であるが、ホルムアルデヒドはゼロレングス架橋剤と称されるものであり、極めて近い距離でコンタクトする配列しかDNAペアを作製することができない。このため、ホルムアルデヒドに代えて、EGS(16.1Å)又はDSG(7.7Å)など比較的長い架橋剤を用いれば、より複雑な高次構造を有する大型タンパク質複合体を捕捉でき、ホルムアルデヒドで架橋可能な距離よりも離れて座するDNA配列間の相互作用を検出するためのDNAペアを作製することができる。
【0039】
この架橋ステップは、例えば、氷上冷却から37℃の範囲の環境において、0.05vol%~4vol%のホルムアルデヒド溶液中で細胞をインキュベートすることによって行ってもよく、より好適には、0.5vol%~2.5vol%、更に好適には1vol%~2vol%のホルムアルデヒド溶液中でインキュベートすることによって行っても良い。インキュベート時間は30秒~16時間であり、好適には、5分~30分であり、より好適には10分~15分である。なお、このような架橋ステップを設けた場合は、高次構造を十分に固定してDNA配列の位置関係を的確に固定し、この後のステップにて得られるDNAペアの精度を向上させることができる。
【0040】
(切断ステップ)
切断ステップは、主に、予め選択したタンパク質(標的タンパク質)が接するDNA配列において当該タンパク質近傍のDNA鎖を、酵素を用いて切断するステップである。この切断ステップに際しては、酵素の透過性を高め、細胞核内に存在するDNAの切断を良好に行うべく、核膜の透過処理を実施しても良い。
【0041】
透過処理には、界面活性剤が用いられる。好適には、非イオン性または両イオン性の界面活性剤が用いられる。かかる界面活性剤の例としては、例えば、ジギトニン、オクチルフェノールエトキシレート類(商品名:トリトンX100、トリトンX114、ノニデットP40等)、n-ドデシル-β-D-マルトシド(DDM)、ポリオキシエチレンソルビタンモノラウレート(ツイン20)、ポリオキシエチレンソルビタンモノオレエート(ツイン80)、3-[(3-コールアミドプロピル)ジメチルアンモニオ]-1-プロパンスルホナート(商品名CHAPS)等が挙げられるがこれらに限られるものではない。界面活性剤を用いた透過処理の前または後において、液体窒素を用いて細胞を処理してもよい。
【0042】
この切断ステップにおいては、相互作用するDNA配列に接するタンパク質、例えば、相互作用する一方または双方のDNA配列に会合し、または複合体を形成し、あるいは付着して、これらのDNA配列に対し空間的に近傍に位置するタンパク質を標的とする。
【0043】
かかるタンパク質としては、クロマチン(DNA鎖)に接する各種のタンパク質を対象とすることができるものであり、分析対象に応じて適宜選定されるものである。その一例としては、例えば、DNAの三次元構造形成に密接に関与していることが知られているタンパク質として、CTCFが例示され、RNAの転写に関与するタンパク質としてRNAポリメラーゼII、ヌクレオソーム構造に必須であるヒストンH3,ヒストンH4、ヒストンH2A、ヒストンH2Bが挙げられるがこれらに限られるものでない。
【0044】
そして、選定されたタンパク質(標的タンパク質)に直接結合する一次抗体が準備される。かかる一次抗体としては、標的タンパク質に対応するモノクローナル抗体やポリクローナル抗体であるが、例えば、上記に例示したタンパク質に対応するものとして、CTCF、RNAポリメラーゼII、ヒストンH3,ヒストンH4、ヒストンH2A、ヒストンH2Bなどに対するモノクローナル抗体が用いられる。また、当該抗体は、必ずしも、完全な抗体である必要はなく、同様の機能を有する抗体の機能的断片、例えば標的タンパク質に結合する機能を有する断片であってもよい。
【0045】
切断ステップでは、上記標的タンパク質を捕捉した一次抗体を介して標的タンパク質に結合する酵素、換言すれば、一次抗体を介してDNA鎖に結合する酵素によって、DNA切断が行われる。なお、酵素が一次抗体を介してDNA鎖に結合するとは、DNA鎖に付着、結合する標的タンパク質に直接結合する一次抗体に、直接的または間接的に酵素が結合することを意味している。かかる酵素は、標的タンパク質の近傍位置のDNA鎖を切断する能力を有するものであれば、公知のものを適宜選択して使用する事ができ、例えば、ミクロコッカルヌクレアーゼ(MNase)、ディーエヌアーゼI(DNaseI)、トランスポザーゼ(transposase)などの各種エンドヌクレアーゼや、ヒンディースリー(HindIII)、MboI、DpnIIなどの各種制限酵素が例示できるがこれらに限られるものではない。好適には、配列非特異的に切断を行うエンドヌクレアーゼが用いられ、より好適には、MNaseが用いられる。なお、近傍とは標的タンパク質の位置から1000bp以内、例えば、900bp以内、800bp以内、700bp以内、600bp以内、500bp以内、400bp以内、300bp以内、200bp以内、好ましくは200bp以内を意味し、例えば、トランスポザーゼ、MNase、DNaseI、HindIII、MboI、DpnIIを用いれば、標的タンパク質の位置から1000bp以内を切断でき、また、MNase、DNaseI、HindIII、MboI、DpnIIを用いれば、標的タンパク質の位置から200bp以内を切断することができる。
【0046】
この切断ステップによって、切断酵素の特性に応じた塩基長のDNA断片を取得でき、例えばMNaseを用いた場合には、200bp程度(モノヌクレオソームサイズ程度)のDNA断片を取得できる。なお、得られるDNA断片長は、切断ステップの実施条件に応じて変化し、モノヌクレオソームサイズのものに、その倍数のものが含まれたものとなり得る。
【0047】
当該酵素が、DNA切断能を保有していても、抗体への結合能を本来的に有していない場合には、抗体結合能が酵素に付与された後に、切断ステップに用いられる。酵素に抗体結合能を付与する方法としては、各種方法を用いる事が可能であり、例えば、抗体のFc領域に高い親和性を示すタンパク質を、使用する酵素に融合した融合タンパク質とすることで抗体結合能を付与する事や、上述した一次抗体に特異的に結合する二次抗体に対し各種の化学反応で酵素を複合化する方法などが例示される。この抗体に親和性を示すタンパク質としては各種の公知のタンパク質を用いる事ができ、例えば、プロテインA、プロテインG、プロテインLが例示できるが、これらに限られるものではない。
【0048】
プロテインAなどのタンパク質と酵素との融合反応としては、公知方法を適宜用いる事ができ、例えば、過ヨウ素酸法、グルタルアルデヒド法、マレイミド法、NHS(N-ヒドロキシスクシンイミド)法を用いることができるがこの限りでない。
【0049】
また、生物工学的手法を用いて融合タンパク質を生成することもできる。具体的には、融合タンパク質を産生するベクター(プラスミド)に切断酵素の配列をコードする遺伝子を導入してベクターを構築し、又は、すでに構築されたベクターを、大腸菌に導入し、形質転換された大腸菌を培養すると共に、融合タンパク質の発現を誘導すれば、目的の融合タンパク質を生成することができる。
【0050】
例えば、プロテインAとMNaseとの融合タンパク質(pA-MN)であれば、公知のベクターpK19pA-MN等を用いて調整する事ができる。また、公知のプラスミドpET21a等にプロテインAおよびMNaseをコードするDNAを導入してpA-MNの発現ベクターを構築し、これを用いて調整できる。この場合においては、発現させたpA-MNを可溶性とすることができ、大腸菌を破砕することによってできた懸濁液を、遠心分離すれば、pA-MNは得られる上清中に含まれる。そして、この上清中からpA-MNをIgG-セファロースビーズ等にて吸着、洗浄、溶出の操作を経ることで回収することができる。
【0051】
本切断ステップで、DNA鎖に一次抗体を介して酵素を結合させる手順としては、1)DNA断片に一次抗体を加えた後、一次抗体に結合能を有する酵素を逐次的に添加し、一次抗体に酵素を結合させる。2)一次抗体と酵素とを予め結合させたものを加えることにより標的タンパク質の捕捉と共に一次抗体を介したDNA鎖への酵素の結合を行う。3)一次抗体添加後に二次抗体を加え、この二次抗体に対する結合能を有する酵素を加える。4)一次抗体添加後に、二次抗体と酵素とを予め結合させたものを加える等のいずれかの手順で実施されてよい。なお、一次抗体と酵素とを予め結合させたものをDNA断片に加える手法の場合には、項目9に記載の一次抗体結合ステップと酵素結合ステップとは一度に実施されることとなる。
【0052】
なお、上記の酵素には、DNA鎖の切断のタイミング(酵素活性)を容易に調整できるものが好適に用いられる。酵素活性の調整方法としては、反応温度、光照射、金属イオン濃度などが例示できる。より、好適には、反応のオンオフの切り替えが容易な金属イオンをエフェクターとするものが選択される。具体的には、例えば、カルシウムイオン依存性のものとしてMNaseやDNaseIが例示でき、マグネシウムイオン依存性のものとしてトランスポザーゼが例示できる。かかる酵素では、対応するイオンが添加されたタイミングでDNA切断が実行される。
【0053】
この切断ステップによって、予め選択した任意の標的タンパク質の近傍のDNA鎖が切断され、この後のライゲーションステップにて近接する切断端を結合させると、ゲノム上の配列とは異なる並びの配列となるDNAペアが作製される。
【0054】
本発明では、この切断ステップは、ゲノムワイドにランダムな切断を行うのではなく、標的タンパク質特異的に切断が行われる。このため、得られるDNAペアの一群(DNAライブラリ)は、従来のゲノムワイドにランダム切断して得られるDNAライブラリと比べて、高濃度に(原理上は、得られるDNAペア全て)所望の解析対象のDNA配列が存在するものとなる。つまり、ノイズとなる不要なDNA断片(ライゲーション副生成物、夾雑物)が極めて少ない高品質な分析サンプルを供し得ることとなる。
【0055】
従来のHiC法等では、シーケンス後のデータ解析で、意図しないDNAペアの情報を削除し解析データの精度向上を図っており、また、最終データの精度向上のためにシーケンスを膨大なリード数によって実行し、これを実現するために相当量の細胞数が必要となっている。一方、本発明は、シーケンスを行うDNAライブラリそのものの品質を向上させることで、データ信頼性を向上させるのである。このため、分析に必要な細胞数を抑制して、少量サンプルでの解析を可能とできる。また、増幅を行う場合であっても、分析サンプル中に含まれる目的物のDNAペアの割合が極めて高いため、夾雑物によるPCR時の増幅バイヤスの影響を抑制できる。更に、不要なDNA断片のリードを抑制でき、分析にかかる労力、コストを大幅に低減できる。
【0056】
なお、DNAペアを効率的に選別し濃縮するため、切断末端に、例えばビオチンを標識として付加する標識付加ステップを設けてもよい。本発明において、より好適には、免疫沈降法を用いずにDNAペアを濃縮するべく、この標識付加ステップを設ける事が望ましい。この標識付加ステップは、後述する選別ステップの前までに設けられればよい。このため、例えば、切断ステップに続けて設けてもよく、ライゲーションステップにおいてライゲーションと共に行われるように設けてもよい。このように、標識付加ステップを設ければ、ビオチン等の標識によって解析対象のDNAペアを濃縮できる。言い換えれば、免疫沈降法によらずに、必要なDNAペアを濃縮できるので、免疫沈降法において必要となる核画分の可溶化を行う必要がなく、形成されるDNAライブラリの品質をより一層向上させることができる。
【0057】
切断されたDNA末端のビオチン付加は、常法に従って行うことができる。例えば、切断されて生じた突出末端は、デオキシヌクレオチド三リン酸(dNTP)を基質に用いて伸張反応を行い平滑末端化することができるが、このdNTPをビオチン化しておけば、切断末端にビオチンを付加できる。ビオチン化には各種の市販のビオチン化試薬を用いる事できる。なお、酵素切断された切断末端が平滑末端であれば、一度突出末端とした後に、ビオチンを付加すればよい。
【0058】
なお、DNAペアの選別(濃縮)を可能とする物質であれば、切断末端に付加されるものはビオチンに限られるものではないが、DNA構造の固定解消にプロテアーゼ処理を行うことが好適であるところ、他の物質を選択する際においても、プロテアーゼで消化されない(プロテアーゼ耐性を有する)化合物が望ましい。かかる化合物としては、例えば、メチル化されたdATP等が例示される。このメチル化されたdATPを切断末端に付加すれば、抗メチル化dATP抗体でDNAペアを濃縮し得る。また、窒素の三重結合とアルキンとの環状化反応であるclick chemistryに反応する化合物を標識として用いてもよい。かかる化合物の一例として、例えばアルキンを持つ非天然塩基(5-エチニル-2’-デオキシウリジン(EdUなど))が例示できる。切断端に付加されたかかる化合物は、クリック反応によってアジド化合物と結合するため、アジド基を有するアジドビオチンやアジド基修飾フラッグタグなどと結合させることができ、その結果、これらビオチンやフラッグタグに基づいた濃縮を可能とし得る。また、ビオチン化dNTPに代えて、ビオチンが付加された短いDNA(リンカーDNA)を標識として用いてもよく、このリンカーDNAを介してライゲーション反応を行ってもよい。
【0059】
(ライゲーションステップ)
本ステップは、DNAリガーゼを用いてライゲーション反応を行うステップである。ライゲーションとは、2つのDNA断片をリン酸ジエステル結合により連結することを意味する。DNAリガーゼとしては、T4DNAリガーゼ、大腸菌DNAリガーゼ等の常温性DNAリガーゼを用いることができる。また、耐熱性DNAリガーゼを用いてもよい。このライゲーション反応によって、空間的に近位に位置する切断末端が連結することで、ゲノム上での並びでは遠方に位置する配列部分が横並びとなるDNAペアを構築できるのである。
【0060】
なお、このライゲーションステップにおいては、ライゲーションされなかったものを除去するため、エキソヌクレアーゼを用いて末端に存在するdNTP(ビオチン化dNTP)を消化しても良い。
【0061】
(固定解消ステップ)
本ステップは、ライゲーション反応を終了したサンプルに対し、プロテアーゼを用いることで、サンプル中に含有されるタンパク質を分解する工程である。これにより、DNAに結合するタンパク質、核膜などを含むサンプル中のタンパク質が分解される。相互作用するDNA配列は、標的タンパク質がDNA鎖に付着している結果として、近位に配置されるものとなっているが、ライゲーション反応によってDNAペアが形成された後は、タンパク質分子による固定が不要となるため、本ステップによりタンパク質を消化する。これによって、タンパク質によるDNA構造の固定(DNA配列の空間構造の固定)が解消される。このタンパク質消化に伴い、架橋ステップに基づく架橋による固定も解消される。これにより、相互作用するDNA配列が並びあって連結されたDNAペア(ライゲーション産物)の一群が得られる。なお、標識付加ステップによって付加した標識に基づいて後述する選別ステップを行う場合(クロマチン免疫沈降を行わない場合)には、本ステップによって核膜を消化することで、目的のサンプル(DNAペア)を核内から取り出すことができる。故に、核画分の可溶化工程を不要にできる。
なお、本ステップは、ライゲーションステップに引き続いて行うことが望ましいが、ライゲーションの後、PCRによる増幅操作または配列決定ステップの前までに設けられ、また、本発明において、クロマチン免疫沈降を行うステップが設けられる場合には、免疫沈降前に本ステップは設けられない。また、選別ステップにおいて断片化を行うことに伴いDNA構造の固定が解消される場合には、後述の選別ステップおける断片化処理が固定解消ステップに該当する。なお、かかる場合であっても、架橋ステップに基づく架橋の解消が必要な場合には、本ステップは固定解消ステップの一部を構成する。また、断片化に伴いDNA構造の固定が解消されるため本ステップによるDNA構造の固定解消が不要な場合であっても、核膜等の消化を行うべく本ステップ同様、プロテアーゼにて消化を行うステップを設けてもよい。
【0062】
(選別ステップ)
DNAペアを選択的に回収するステップである。本ステップにおいては、PCR効率、シーケンス効率を考慮して、ライゲーション産物として得たDNAペアの更なる断片化を行うことが好ましい。また、断片化と共にアダプター配列を断片末端に導入してもよい。
【0063】
一実施形態において、ライゲーションされたDNAペアを選別するステップは、ライゲーションされた当初配列の状態のものを選別することのみならず、ライゲーションの後、加工や処理がされたDNA断片や、その処理によって更に断片化されたDNA断片や、任意のDNA配列が付加されたDNA断片を選別することを含む。例えば、PCRステップを実施する場合には、ここで付加されるDNA配列は短いアダプター配列であってもよい。これにより、その後のPCRステップにおいて、短いアダプター配列の相補配列とインデックス配列とが含まれたプライマーを、DNAペア(DNA断片)に結合した短いアダプター配列にハイブリダイズさせ、増殖反応の実施と共にDNAペアにインデックス配列を付加することができる。また、PCRステップを行わない場合には、インデックス配列と短いアダプター配列との両方を備えたDNA配列をアダプター配列(フルアダプター)としてDNAペアにライゲーションさせ、選別ステップにて選別される配列が付加されたDNAを構築してもよい。例えば、アダプター配列(フルアダプター)やインデックス配列は、一般的な公知のシークエンサーシステムで用いられるフローセルに固定されたDNA断片と相補的なDNA配列を含んでよい。
この選別ステップでは、例えば、上記標識付加ステップにより、各ライゲーション産物の接合部に標識(ビオチン)が付加されていると、標識をターゲットとしてDNAペア(DNA断片)の回収を行う。この回収には、標識に対して親和性を示すタンパク質であれば、特に限定されず、市販の各種タンパク質を使用する事ができる。例えば、標識がビオチンである場合には、好適には、ビオチンに極めて高い親和性を有するアビジン、ストレプトアビジン、脱グリコシル化アビジンといったビオチン結合性タンパク質が挙げられる。また、当該ステップにおいては、これらのタンパク質がビーズなどの担持体に固定化されたものを用いて選別操作を行ってよい。本ステップで解析対象のDNAペアが選択的に回収され、DNAライブラリが構築される。また、例えば、標識に基づき回収したDNAペアで構築されたDNAライブラリに対し、オリゴヌクレオチドによるハイブリダイゼーションを行って、DNAペアの選別をさらに行ってもよい。
一実施形態において、DNAペア(DNA断片)を選別するステップは、ライゲーションの後に、電気泳動によってDNAペアを選別することを含む。例えば、アガロースゲルを用いて電気泳動し、DNA鎖長でライゲーション産物を選別する、標識を用いなくてよい手法が挙げられる。
【0064】
(配列決定ステップ)
上記にて構築されたDNAライブラリに対し、所望により、常法に従ってPCRを行い、各DNAペア(DNA断片)の増幅を行う。そして得られたPCR産物に対して、配列決定を行う。DNA断片に対し増幅を行う場合には、その増幅は、DNA断片の全長を増幅することのみを意味するのではなく、その一部を増幅することも含まれる。上記のPCRによって濃縮されたDNAペア(DNA断片)に対し、あるいは選別ステップで構築されたDNAライブラリ中のDNAペアに対し、前処理を行い、超並列DNA配列解析装置やサンガー法により、DNAペアの両端の塩基配列を解析し配列決定を行う。解析する塩基配列の長さは、1塩基からDNA断片全長でも良い。解析する塩基配列の長さは、任意に設定でき、例えば、数塩基以上、10塩基以上、15塩基以上、20塩基以上、30塩基以上であってよい。好適には、例えば、30-100塩基程度である。
【0065】
(解析方法)
得られたシーケンスデータは、従来のHiC法の解析手法を用いて解析できる。まず、得られた塩基配列に対し、両端の配列をそれぞれ独立に参照用ゲノム配列(マウスの細胞ならマウス、ヒト細胞ならヒト)にマッピングする。本法では、近接したDNAを人為的に結合させてDNA断片(DNAペア)を作製しているために、末端はそれぞれ独立にゲノムにマッピングできるが、DNA断片としてはマッピングできないことがある。例えば1番染色体のDNAと2番染色体のDNAが近接しており、本法によりDNA断片を作製した場合、得られたDNA断片は、一端に1番染色体の配列、他方に2番染色体の配列を有することになる。
【0066】
本法の解析においては、HiC法一般の実験と同様の解析手法が適用でき、こういった配列(参照ゲノム配列には存在し得ない人工配列)を選択的に抽出し、さらなる高次解析を行うのである。
【0067】
例えば、コンタクトマップ(接触確率行列)をヒートマップとして生成することによって、ゲノム上のある位置に存在するDNA配列と他の位置に存在するDNA配列との相互作用が可視化される。コンタクトマップは対称行列である。X座標とY座標に同じゲノム配列を並べ、シーケンスデータから位置Xiと位置Yjにマップされるペアエンドリードがあると、その(Xi,Yj)のカウントをひとつ増やす。したがって、各座標で値が大きい場合、接触確率が高い領域ペアであることを意味するものとなる。この座標の値は、これに応じた色の濃淡でヒートマップに表示することができる。
【0068】
(サンプル量)
本法は、従来法に比べて夾雑物の極めて低減された解析サンプルを得ることができる結果、解析に必要な細胞数を1~1010個程度とすることができ、好適には1~10個程度、より好適には、1~10個程度とすることができる。なお、細胞数の上限は必ずしも制限されず、1010個以上であっても良い。
【0069】
(本発明の方法の例示的な形態)
一実施形態において、本発明の方法は、従来のHiC法に代わる手法として用いることができるものである。かかる場合には、上述した、架橋ステップが適宜設けられた後、切断ステップ、標識付加ステップ、ライゲーションステップ、固定解消ステップ、選別ステップ、配列決定ステップを順次設けて構成される。本法によれば、上述したように固定解消ステップにて核膜の消化も行う事ができる。クロマチン免疫沈降法では、核画分の抽出が必要になるために、原理上不溶性画分の発生を避けられない。ゲノムワイドに制限酵素でクロマチンを切断し、DNA断片を生じさせる従来の免疫沈降法では、全DNA断片を可溶化させ核外に取り出すことが重要となる。不溶性画分の核内の残存は、データの誤差となるからである。一方、上記の標識付加ステップによって付加した標識に基づいて解析対象のDNA配列を選別するとすれば、ライゲーション処理の後にプロテアーゼを用いてDNAペアを選別することができ、(核から溶出しない不溶性画分がない)良好な分析サンプルを作製できることとなる。一実施形態においてクロマチン免疫沈降ステップを設け目的のDNA断片を補足する場合、目的のDNA断片のみが上清中に放出されればよいので、解析対象外の部分となる未消化の長鎖のDNA断片は、核内に不溶性画分として残存しても、得られるデータに影響がない。このため、従来の免疫沈降法において生じるような核の可溶化工程での不溶性画分に起因するデータの信頼性低下を抑止し得る。また、従来HiC法とは異なり、特定のタンパク質特異的にDNA切断を行うので、分析に供されるDNAライブラリにおいては、所望の解析対象のDNA配列が濃縮されており、ノイズとなる不要なDNA断片が極めて少ない高品質な分析サンプルを供し得ることとなる。
更には、得られたDNAペアのシーケンス結果からDNA切断部位を独立に解析すると、標的タンパク質の局在を検出することができる。すなわち、標的タンパク特異的にDNA切断を行う切断ステップを備えた本方法によれば、標的タンパク質によって介在される核内のDNA三次元構造と標的タンパク質の核内での局在の両方を一度に解析できるのである。
また、標的タンパク質特異的にDNA切断を行う切断ステップの前に、DNAに結合するタンパク質をビオチン化するタンパク質ビオチン化ステップを設けると共に、切断ステップの後に、核膜透過処理によって切断された断片を可溶化する可溶化ステップ(透析工程)を設けて、核内から、タンパク質DNA複合体を上清中に放出させ、これをストレプトアビジン等の担持ビーズを用いてプルダウンしてから、担持されたビーズ上にてライゲーションを行うようにライゲーションステップを構成し、更に、順次、固定解消ステップ、選別ステップ、配列決定ステップを実施してシーケンスデータを得てもよい。
【0070】
別の実施形態において、本発明の方法は、クロマチン免疫沈降を行うとして、核分画の可溶化ステップとクロマチン免疫沈降ステップとを備えて構成しても良い。なお、核分画の可溶化ステップとクロマチン免疫沈降ステップとは常法にしたがって実施できる。かかる場合には、一例として、上述した、架橋ステップが適宜設けられた後、切断ステップ、標識付加ステップ、ライゲーションステップ、核分画の可溶化ステップ、クロマチン免疫沈降ステップ、固定解消ステップ、選別ステップ、配列決定ステップを順次設けて構成することで、従来のHiChIPに代わるものとして、従来のHiChIPよりも、分析に供するDNAライブラリの作製段階において夾雑物を低減し得る手法を提供できる。この実施形態の方法の一例は、上記項目6に記載の三次元DNA構造相互作用分析方法に該当する。この実施形態の方法に含まれるクロマチン免疫沈降ステップは、上記項目6に記載の捕捉ステップに含まれる。
【0071】
別の実施形態において、本発明の方法は、上述した、架橋ステップが適宜設けられた後、切断ステップ、核分画の可溶化ステップ、クロマチン免疫沈降ステップ、ライゲーションステップおよび標識付加ステップ、固定解消ステップ、選別ステップ、配列決定ステップを順次設けて本発明を構成することで、ChIA-PETに代わるものとして、ChIA-PETよりも、分析に供するDNAライブラリの作製段階において夾雑物を低減し得る手法を提供できる。なお、本例においてはクロマチン免疫沈降ステップに基づいてDNAペアの濃縮が可能であるため選別ステップを不要としてもよい。また、本実施形態においても適宜追加のステップを設けてもよく、例えば、ライゲーションしたDNA断片に対し、適宜の再断片化や配列の付加を行うステップを設けることができ、かかるステップを経た後のDNA断片に対して増幅を行ってもよい。この実施形態の方法の一例は、上記項目7に記載の三次元DNA構造相互作用分析方法に該当する。この実施形態の方法に含まれるクロマチン免疫沈降ステップは、上記項目7に記載の捕捉ステップに含まれる。
【0072】
一実施形態において、本発明の分析用キットは、任意のタンパク質に結合する一次抗体と、一次抗体に結合するDNA切断酵素と、DNA切断酵素による切断によって生じるDNA末端を識別するための標識と、標識されるDNA末端をライゲーションする酵素とを含む。分析用キットに含まれる一次抗体、DNA切断酵素、標識およびDNA末端をライゲーションする酵素は、本発明のDNAの三次元構造に基づくDNA配列の相互作用を分析する方法において使用される材料である。
なお、一次抗体とDNA切断酵素とは、別々に含まれていてもよく、一次抗体とDNA切断酵素とが結合された状態で含まれていてもよい。
【実施例
【0073】
以下、本発明について、更に実施例に基づき詳細に説明するが、以下の実施例は本発明の一例であり、これらに限られるものでない。
【0074】
実施例において、以下の試薬が使用された。
(試薬)
PBS(Phosphate Buffered Saline; リン酸緩衝食塩水)
滅菌水(Nacalai Tesque社)
1M CaCl2 (Calcium Chloride)/滅菌水
1M DTT(Dithiothreitol)
1M Glycine/PBS
0.5% BSA/滅菌水(0.5wt.%ウシ血清アルブミン(BSA)水溶液)
75% Ethanol/滅菌水
Tween-20(Chem Cruz社)
DMSO (Dimethyl Sulfoxide)
5 M NaCl
1 M MgCl2
1 M Tris-HCl (pH7.4)(Sigma-Aldrich)
0.5 M EDTA (pH8.0) (Ambion社)
0.5 M EGTA (A.G. Scientific社)
EGTA(Sigma-Aldrich社
Triton(登録商標) X-100 (Nacalai Tesque社)
IGEPAL(登録商標) CA-630 (Sigma-Aldrich社)
Spermidine Trihydrochloride (Sigma-Aldrich社)
Digitonin (Nacalai Tesque社)
cOmpleteTM, EDTA-free Protease Inhibitor Cocktail Tablets (Roche社)
1mol/l-HEPES buffer solution(Nacalai Tesque社)
【0075】
RSB(Resuspension Buffer; 再懸濁バッファー)
組成:
10mM Tris-HCl(pH7.4)
10mM NaCl
3mM MgCl
0.1% Tween-20
滅菌水
Buenrostro JD, Wu B, Chang HY, Greenleaf WJ: ATAC‐seq: a method for assaying chromatin accessibility genome‐wide. Curr Protoc Mol Biol. 2015;109:21-9およびCorces MR, Trevino AE, Hamilton EG, Greenside PG, Sinnott-Armstrong NA, Vesuna S, Satpathy AT, Rubin AJ, Montine KS, Wu B, et al. An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues. Nat Methods. 2017;14:959-62を参考にして、RSBを調製した。
【0076】
5% Digitonin
組成:
50mg Digitonin(Nacalai Tesque社)
500μl DMSO
500μl 滅菌水
【0077】
Lysis Buffer
組成:
0.1% IGEPAL CA-630
0.01% Digitonin
RSB
Buenrostro JD, Wu B, Chang HY, Greenleaf WJ: ATAC‐seq: a method for assaying chromatin accessibility genome‐wide. Curr Protoc Mol Biol. 2015;109:21-9およびCorces MR, Trevino AE, Hamilton EG, Greenside PG, Sinnott-Armstrong NA, Vesuna S, Satpathy AT, Rubin AJ, Montine KS, Wu B, et al. An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues. Nat Methods. 2017;14:959-62を参考にして、Lysis Bufferを調製した。
【0078】
Wash Buffer
組成:
1M HEPESを1ml、5M NaClを1.5ml、2M Spermidineを12.5μl加え、50mlまで滅菌水を加え、さらにcOmpleteTM, EDTA-free Protease Inhibitor Cocktail Tablets (Roche社)を1 tablet加えることによって、Wash Bufferを調製した。
Skene et al., 2018. Targeted in situ genome-wide profiling with high efficiency for low cell numbers. Nature Protocols. 2018. Vol. 13. No. 5を参考にして、Wash Bufferを調製した。
【0079】
Digitonin Buffer
組成:
0.02% Digitonin(実施例3においては0.0375%)
Wash Buffer
Skene et al., 2018. Targeted in situ genome-wide profiling with high efficiency for low cell numbers. Nature Protocols. 2018. Vol. 13. No. 5を参考にして、Digitonin Bufferを調製した。
【0080】
Antibody Buffer
Digitonin Bufferを2ml、0.5M EDTA(pH8.0)を8μl加えることによって、Antibody Bufferを調製した。
Skene et al., 2018. Targeted in situ genome-wide profiling with high efficiency for low cell numbers. Nature Protocols. 2018. Vol. 13. No. 5を参考にして、Antibody Bufferを調製した。
【0081】
Mbuffer#2
組成:
50mM NaCl
10mM Tris-HCl(pH7.4)
10mM MgCl
滅菌水
Rando et al., 2016. Micro-C XL. Protocol Exchangeを参考にして、Mbuffer#2を調製した。
【0082】
10x NEBuffer#2
組成:
1M NaCl
200mM Tris-HCl(pH7.4)
200mM MgCl
滅菌水
Rando et al., 2016. Micro-C XL. Protocol Exchangeを参考にして、10x NEBuffer#2を調製した。
【0083】
10X T4 DNA ligase buffer(ATPなし)
組成:
500mM NaCl
583mM Tris-HCl (pH7.4)
200mM MgCl2
117mM DTT
Rando et al., 2016. Micro-C XL. Protocol Exchangeを参考にして、10X T4 DNA ligase bufferを調製した。
【0084】
1X Bind and Wash Buffer
組成:
5mM Tris-HCl, pH7.4
0.5mM EDTA
1M NaCl
滅菌水
Gabdank et al. 2016. A streamlined tethered chromosome conformation capture protocol. BMC Genomics. 2016. 17:274を参考にして、1X Bind and Wash Bufferを調製した。
【0085】
2X Bind and Wash Buffer
組成:
10mM Tris-HCl, pH7.4
1mM EDTA
2M NaCl
滅菌水
Gabdank et al. 2016. A streamlined tethered chromosome conformation capture protocol. BMC Genomics. 2016. 17:274を参考にして、2X Bind and Wash Bufferを調製した。
【0086】
(実施例1)
プロテインAとMNaseとの融合タンパク質(pA-MN)発現および精製
大腸菌株One Shot(登録商標) BL21(DE3)(invitrogen)をプラスミドpET11a-pA-MN(pET11aにpA-MNを発現するDNA配列を導入して構築したベクター)で形質転換し、100mg/Lのカルベニシリン(ナカライテスク社)を含有するLB培地(LB-Medium, Lennox)(フナコシ社)中で37℃で培養した。指数関数的に増殖している細胞に0.1mM イソプロピル-β-チオガラクトピラノシド(IPTG)を加え、培養温度条件を18℃に変更し培養を再開することによって、pA-MNの発現を誘導した。19時間の誘導期間後、500mlの培養物から細菌を回収し、5mM DTTを補った30ml TEN(10mM Tris-HCl[pH8.0]、1mM EDTA、および150mM NaCl)に回収した。1mgチキン卵白リゾチーム(Sigma社)を添加し、氷上で10分間インキュベートした。超音波ホモジナイザー(BIORUPTOR(登録商標)II Type6、ソニック・バイオ社)を用いて超音波処理を繰り返すことによって、細胞を破砕した。抽出物を4℃で遠心機(MX-301、TOMY社)中で15000×Gで20分間、遠心分離によって清澄化した。pA-MNの大部分は可溶性上清(S1)中に回収された。ビーズ(IgG-セファロース6 Fast Flow、GE healthcare社)を使用してS1からpA-MNを精製した。30mlの清澄化されたS1を600μlのIgG-セファロースビーズと共に16時間振とうしながらインキュベートした。エコノパッククロマトグラフィー用カラム(BIO-RAD社)に通して、ビーズを回収し、0.03% 界面活性剤(EMPIGEN(登録商標) BB 界面活性剤, SIGMA社)を補ったTENで2回、0.03% Empigen(登録商標)BBを含有するTEN500(500mM NaClを補ったTEN)で2回、そして5mM 酢酸アンモニウム(pH5.0)で1回洗浄した。pA-MNを0.5M 酢酸/酢酸アンモニウム(pH3.4)で溶出した。溶出液に1 mol/L NaOH水溶液を加え中和した。中和した溶出液と同量のグリセロールを添加し最終濃度50%に調製した。500mlの出発培養物あたり約1mgの高度に精製されたタンパク質を得た。長期間のアリコートは-80℃で、短期間のアリコートは-30℃で保存した。
【0087】
HEK293T細胞を、2vol%ウシ胎児血清(FBS)を含有するPBS(FBS-PBS)に懸濁した後、細胞数が10個となるようにチューブに注入し、室温で5分、遠心力x300gで遠心分離した。分離した細胞に、ホルムアルデヒド(商品名:Formaldehyde solution, 型番F8775, Sigma-Aldrich社)1wt.%を含有するFBS-PBS 1mlを加えて懸濁し、穏やかに攪拌した。その後、終濃度が250mMとなるように1MのグリシンPBS溶液の340μlを加え、室温5分間の混合操作を行い、続いて、氷上にて15分間保持しつつ適宜倒拌にて攪拌を行った。その後、4℃、遠心力×800g、5分間の遠心分離操作と5mlのPBSによる再懸濁とを複数回実施して洗浄を行った。なお、遠心分離操作では、適宜、分離した上清は除去される。最後の洗浄では、遠心分離後の上清(PBS)をできるだけ除去し、次いで、液体窒素にてサンプルを急冷した。
【0088】
上記手順で調整した細胞サンプルを溶解し、1mlのLysis Bufferにて再懸濁し、30分間、定期的に倒拌しつつ氷上に保持した。その後、1mlのRSBを加えてから、4℃、5分間、遠心力×800gで遠心分離した後、洗浄用バッファー(Wash Buffer)1.5mlにて洗浄を行った。洗浄用バッファーの添加と遠心分離操作とによる洗浄作業を複数回実施した後、4℃、5分間、遠心力×800gで遠心分離した沈殿物を回収し、抗H3K4me3抗体(型番ab8580、Abcam社)1μgを含むAntibody Buffer溶液500μlに再懸濁し、4℃で一晩、反応させた。抗H3K4me3抗体は、H3K4me3(4番目のリジンがトリメチル化されたヒストンH3)の抗体である。
【0089】
反応後のサンプルに対して、4℃、遠心力×500gで5分間の遠心分離操作と1mlのDigitonin Bufferでの再懸濁とを2回繰り返して行ってから、4℃、遠心力×500gで5分間の遠心分離操作の後、40μlのDigitonin Bufferにて再懸濁し、更に、Digitonin Bufferで1000倍希釈した(適宜に濃度調整した)pA-MNを加えて4℃にて1時間反応させた。その後、4℃、遠心力×500g、5分間の遠心分離操作と、1mlのDigitonin bufferでの再懸濁とを複数回実施した後、4℃、遠心力×500g、5分間で遠心分離した沈殿物を150μlのDigitonin bufferで再懸濁してから、100mM CaCl溶液を3μl添加して、37℃、20分反応させた。続いて、グリコールエーテルジアミン四酢酸(ethylene glycol tetraacetic acid、EGTA)を終濃度2mMになるように添加し、65℃、10分間処理をした。
【0090】
その後、4℃、遠心力×16000g、5分間の遠心分離操作と、これにより分離されたペレットを1mlのMbuffer#2で懸濁する操作と、を繰り返した後、4℃、遠心力×16000g、5分間の条件で遠心分離したペレットに、0.5wt.%ウシ血清アルブミン(BSA)水溶液を加えたMbuffer#2 (29μl Mbuffer#2+1μl 0.5%BSA/滅菌水)を30μl加えて再懸濁した。
【0091】
続いて、上記手順にて調整したサンプル30μlに対し、脱リン酸化酵素(商品名:Shrimp Alkaline Phosphatase (rSAP), New England BioLabs社)を2.5μl加え、37℃45分間の脱リン酸化処理を行った。反応終了時、脱リン酸化反応は、65℃5分間の熱処理にて停止させた。
【0092】
そして、10倍濃度に調整したNEBuffer#2(10x NEBuffer#2 )3μl、0.5wt.%BSA含有水溶液0.2μl、10mM ATP溶液(商品名 Adenosine 5’-Triphosphate (ATP)、New England BioLabs社)4.5μl、1MのDTT(商品名 UltraPureTM Dithiothreitol (Clelands Reagent)、Invitrogen/Life Technologies社)0.15μl、T4 DNAポリメラーゼ(商品名 T4 DNA Polymerase(5U/μl)、Thermo Scientific社)1.5μl、T4 ポリヌクレオチドキナーゼ(10U/μl)(商品名 T4 Polynucleotide Kinase、New England BioLabs社)3μlを加えて37℃で7分間反応させ、その後、10倍濃度のT4 DNA ligase buffer 6μl、0.5wt.%BSA含有水溶液0.3μl及び、ビオチン化dNTP(ビオチン化dATP(商品名:Biotin-14-dATP,Invitrogen社)、ビオチン化dCTP(商品名:Biotin-14-dCTP,Invitrogen社) がそれぞれ終濃度0.1mM、dGTP(商品名:dGTP,東洋紡)、dTTP(商品名:dTTP, 東洋紡))がそれぞれ終濃度0.05mMとなるよう加えて、溶液量約100μlとしたサンプルを、25℃25分、12℃15分の後、4℃保持とする温度条件で反応させた。そして、エチレンジアミン四酢酸(ethylenediaminetetraacetic acid、EDTA)を終濃度30mMとなるように加えて65℃で熱処理し、酵素活性を失活させた。これら一連の処理によってDNA切断末端にビオチンを付加した。
【0093】
そして、熱処理後のサンプル約110μlに、ATP添加10倍濃度T4 Ligase buffer(商品名:10X Buffer for T4 DNA ligase with 10mM ATP; New England BioLabs社) を250μl、0.5wt.%BSA含有水溶液を12.5μl、1MのMgClを3μl、T4 DNAリガーゼ(400U/μl)(商品名 T4 DNA Ligase、New England BioLabs社)を12.5μl加え、更に滅菌水を加えて全容量を2500μlとし、室温で1時間混合した。その後、4℃、遠心力×16000g、10分間で遠心分離操作を行い、ペレットを得た。
【0094】
10倍濃度で調整したNEBuffer#1(商品名 NEBuffer 1, New England BioLabs社)を10μlと、滅菌水89μl、エキソヌクレアーゼIII(Exonuclease III(100U/μl)、New England Biolabs社)1μlを、ペレットに加え、100μlのサンプル溶液とした。これを、37℃で5分間保持し、続いて、溶液の1/200容量のプロテアーゼ(商品名 Proteinase K (Recombinant) Solution、Nacalai Tesque社)を加え、55℃で3時間保持してサンプル溶液中のタンパク質を消化した(脱架橋)。
【0095】
続いて、サンプル溶液にフェノールクロロホルム溶液(商品名:フェノール:クロロホルム:イソアミルアルコール 25:24:1フェノール含量:約33w/w%、Nacalai Tesque社)100μlを加えて十分に混合してから、室温、遠心力×19800g、10分で遠心分離を行う操作を2回実施し、その後、上清を回収した。回収した上清の1/10容量の3M 酢酸ナトリウム溶液(商品名: 3mol/l Sodium Acetate Buffer Solution (pH 5.2)、Nacalai Tesque社)と、2.5倍容量の100%エタノールを入れて混合し、-80℃で1時間静置した。そのサンプルを4℃、遠心力×19800g、15分間の遠心分離操作を行ってDNAを沈殿させた後、上清を除いてから1mlの75%エタノールで洗浄し、再度4℃、遠心力×19800g、15分間の遠心分離操作を実施して上清を除去した。沈殿物を室温20分間で風乾し、50μlのBuffer EB(10 mM Tris・Cl、pH 8.5)(商品名Buffer EB Elution Buffer、Qiagen社)を加えて、37℃、15分保持してDNA溶液を調整した。そして、調整したDNA溶液を、更に精製した後、超微量分光光度計(商品名Qubit 3.0 Fluorometer、Invitrogen/ Life Technologies社)を用いてDNA量を確認した。
【0096】
そして、DNA量が10ngとなるようにDNA溶液をチューブに分注し、2倍濃度のTD buffer25μl(商品名 Tagment DNA Buffer (Nextera(登録商標) DNA Library Prep Kit、Illumina社)とトランスポソーム(Tagment DNA Enzyme 1 (TDE1)、Illumina社)1μlとを加えた後、溶液総量が50μlとなるまで滅菌水を追加した。これを55℃で10分間保持し、タグメンテーションを行った。反応後のサンプルを精製し、50μlのDNA溶液(Buffer EB溶液)とした。
【0097】
次に、ストレプトアビジンビーズ分散液(商品名DynabeadsTM(登録商標) M-280 Streptavidin、Invitrogen/ Thermo Fisher Scientific社)10μlを、500μlの1xBind and Wash bufferで2回洗浄した後、50μlの2xBind and Wash Bufferで懸濁した。この50μlのビーズ分散液を、上記で精製したDNA溶液とチューブ内で混合し、30分間、室温で保持した後、磁気スタンド(1.5ml用)に静置して、上清を除去し、ストレプトアビジンビーズを回収した。そして、回収したビーズを、0.1vol%の非イオン性界面活性剤溶液(商品名Triton X-100、Nacalai Tesque社)を含む500μlの1x Bind and wash bufferで1度洗浄してから、500μlのBuffer EBで更に洗浄した後、40μlのBuffer EBで懸濁して分散液とした。
【0098】
この分散液10μlに、滅菌水10μl、10μMのフォワード側プライマー(N707、Nextera(登録商標) XT Index Kit v.2, Illumina社)2.5μl、リバース側プライマー(S511 、商品名:Nextera(登録商標) XT Index Kit v.2, Illumina社)2.5μl、DNAポリメラーゼ(KAPA HiFi HotStart ReadyMix (2X)、Kapa Biosystems社)25μlを加えてから、PCRを行った。PCRの条件は、72℃5分間の予備加熱、98℃で3分間の初期変性後、98℃で10秒間の熱変性、63℃で30秒間のアニーリング反応、72℃で1分間の伸長反応を13サイクル行い、72℃で1分間の最終伸長反応を行った後、増幅産物を4℃に冷却した。
【0099】
なお、このPCRに先立って、DNA溶液の一部を用いて、リアルタイムPCRを実施し、そのプラトーから最適サイクル数は決定した。
【0100】
得られたPCR産物から、目的のDNAサイズ200-1000bpのDNA断片を回収するため、適宜のサイズセレクションを実施し、最終的に、150bp~1300bpのDNA断片が含まれたDNAライブラリ 10μl(DNA/Buffer EB溶液)を得た。
【0101】
この得られたDNAライブラリを用い、超並列DNA配列解析装置(装置名 NextSeq 500、Illumina社)を使用して、DNA断片(DNAペア)の両端の塩基配列を解析した。解析する塩基配列の長さは、36塩基とした。得られた塩基配列に対し、両端の配列をそれぞれ独立に参照用ヒトゲノム配列にマッピングすることで、DNA配列の相互作用が解析された。得られた結果を図2図3に示す。
【0102】
図2は、H3K4me3が存在するクロマチン領域の相互作用のヒートマップを示す。本発明の手法を確立するために、抗H3K4me3抗体を用いてDNA配列の相互作用解析を行い、H3K4me3が存在するクロマチン領域の相互作用をJuiceboxにより解析した。
図3は、一例として解析された、ゲノム上におけるH3K4me3の局在およびH3K4me3が存在するクロマチン領域の相互作用を示す。ACTB遺伝子領域付近(第7染色体)に局在するH3K4me3とそのヒストン翻訳後修飾が含まれるクロマチン領域間の相互作用をWashU EpiGenome Browserにより解析した。以上の結果より本発明の手法では別で実施したCUT&RUNと同様のH3K4me3の局在が示され、さらにH3K4me3が含まれるクロマチン領域の相互作用が確認された。また、これによって、H3K4me3が結合するDNA配列の相互作用が、同一の遺伝子およびその近傍の領域の間のみならず、離れた遺伝子領域との間においても生じていることが認められた。
図4は、図3の部分的な拡大図である。
図5は、図3の部分的な拡大図である。
【0103】
(実施例2)
ターゲットのタンパク質をCTCFとし、実施例1の抗H3K4me3抗体に代えて5μg抗CTCF抗体(ab70303 Abcam社)を用いたこと、また、リバース側プライマーをS511からS516(商品名:Nextera(登録商標) XT Index Kit v.2, Illumina社)に変更しPCRを13または14サイクルで行ったこと以外は、実施例1と同様にして、DNA配列相互作用の分析を行った。
結果を図6図7に示す。
【0104】
図6は、CTCFが存在するクロマチン領域の相互作用のヒートマップを示す。本発明の手法を確立するために、実施例1と同様に抗CTCF抗体を用いてDNA配列の相互作用解析を行い、CTCFが存在するクロマチン領域の相互作用をJuiceboxにより解析した。
図7は、一例として解析された、ゲノム上におけるCTCFの局在およびCTCFが存在するクロマチン領域の相互作用を示す。ACTB遺伝子領域付近(第7染色体)に局在するCTCFとその転写因子が含まれるクロマチン領域の相互作用をWashU EpiGenome Browserにより解析した。以上の結果、本発明の手法を用いてCTCFが結合するクロマチン領域の相互作用が確認できた。また、これによって、CTCFが結合するDNA配列の相互作用が、同一の遺伝子およびその近傍の領域の間のみならず、離れた遺伝子領域との間においても生じていることが認められた。
図8は、図7の部分的な拡大図である。
図9は、図7の部分的な拡大図である。
【0105】
(実施例3)
分析に使用されるサンプル数を106個の細胞数として、実施例1に基づいて、DNA配列相互作用の分析を行った。具体的には、以下のとおりである。
【0106】
HEK293T細胞を、2vol%ウシ胎児血清(FBS)を含有するPBS(FBS-PBS)に懸濁した後、細胞数が106個となるようにチューブに注入し、室温で5分、遠心力x300gで遠心分離した。分離した細胞に、ホルムアルデヒド(商品名:Formaldehyde solution, 型番F8775, Sigma-Aldrich社)1wt.%を含有するFBS-PBS 1mlを加えて懸濁し、穏やかに攪拌した。その後、終濃度が250mMとなるように1MのグリシンPBS溶液の340μlを加え、室温5分間の混合操作を行い、続いて、氷上にて15分間保持しつつ適宜倒拌にて攪拌を行った。その後、4℃、遠心力×800g、5分間の遠心分離操作と5mlのPBSによる再懸濁とを複数回実施して洗浄を行った。なお、遠心分離操作では、適宜、分離した上清は除去される。最後の洗浄では、遠心分離後の上清(PBS)をできるだけ除去し、次いで、液体窒素にてサンプルを急冷した。
【0107】
上記手順で調整した細胞サンプルを溶解し、1mlのLysis Bufferにて再懸濁し、30分間、定期的に倒拌しつつ氷上に保持した。その後、1mlのRSBを加えてから、4℃、5分間、遠心力×800gで遠心分離した後、洗浄用バッファー(Wash Buffer)1.5mlにて洗浄を行った。洗浄用バッファーの添加と遠心分離操作とによる洗浄作業を複数回実施した後、4℃、5分間、遠心力×800gで遠心分離した沈殿物を回収し、抗CTCF抗体(型番ab70303、Abcam社)1μgを含むAntibody Buffer溶液500μlに再懸濁し、4℃で一晩、反応させた。
【0108】
反応後のサンプルに対して、4℃、遠心力×500gで5分間の遠心分離操作と1mlのDigitonin Bufferでの再懸濁とを2回繰り返して行ってから、4℃、遠心力×500gで5分間の遠心分離操作の後、45μlのDigitonin Bufferにて再懸濁し、更に、Digitonin Bufferで1000倍希釈した(適宜に濃度調整した)pA-MNを加えて4℃にて1時間反応させた。その後、4℃、遠心力×500g、5分間の遠心分離操作と、1mlのDigitonin bufferでの再懸濁とを複数回実施した後、4℃、遠心力×500g、5分間で遠心分離した沈殿物を150μlのDigitonin bufferで再懸濁してから、100mM CaCl溶液を3μl添加して、37℃、20分反応させた。続いて、グリコールエーテルジアミン四酢酸(ethylene glycol tetraacetic acid、EGTA)を終濃度2mMになるように添加し、65℃、10分間処理をした。
【0109】
その後、4℃、遠心力×16000g、5分間の遠心分離操作と、これにより分離されたペレットを1mlのMbuffer#2で懸濁する操作と、を繰り返した後、4℃、遠心力×16000g、5分間の条件で遠心分離したペレットに、0.5wt.%ウシ血清アルブミン(BSA)水溶液を加えたMbuffer#2 (29μl Mbuffer#2+1μl 0.5%BSA/滅菌水)を30μl加えて再懸濁した。
【0110】
続いて、上記手順にて調整したサンプル30μlに対し、脱リン酸化酵素(商品名:Shrimp Alkaline Phosphatase (rSAP), New England BioLabs社)を2.5μl加え、37℃45分間の脱リン酸化処理を行った。反応終了時、脱リン酸化反応は、65℃5分間の熱処理にて停止させた。
【0111】
そして、10倍濃度に調整したNEBuffer#2(10x NEBuffer#2 )3μl、0.5wt.%BSA含有水溶液0.2μl、10mM ATP溶液(商品名 Adenosine 5’-Triphosphate(ATP)、New England BioLabs社)4.5μl、1MのDTT(商品名 UltraPureTM Dithiothreitol (Clelands Reagent)、Invitrogen/Life Technologies社)0.15μl、T4 DNAポリメラーゼ(商品名 T4 DNA Polymerase(5U/μl)、Thermo Scientific社)1.5μl、T4 ポリヌクレオチドキナーゼ(10U/μl)(商品名 T4 Polynucleotide Kinase、New England BioLabs社)3μlを加えて37℃で7分間反応させ、その後、10倍濃度のT4 DNA ligase buffer 6μl、0.5wt.%BSA含有水溶液0.3μl及び、ビオチン化dNTP(ビオチン化dATP(商品名:Biotin-14-dATP,Invitrogen社)、ビオチン化dCTP(商品名:Biotin-14-dCTP,Invitrogen社) がそれぞれ終濃度0.1mM、dGTP(商品名:dGTP,東洋紡)、dTTP(商品名:dTTP, 東洋紡))がそれぞれ終濃度0.05mMとなるよう加えて、溶液量約100μlとしたサンプルを、25℃25分、12℃15分の後、4℃保持とする温度条件で反応させた。そして、エチレンジアミン四酢酸(ethylenediaminetetraacetic acid、EDTA)を終濃度30mMとなるように加えて65℃で熱処理し、酵素活性を失活させた。これら一連の処理によってDNA切断末端にビオチンを付加した。
【0112】
そして、熱処理後のサンプル約110μlに、ATP添加10倍濃度T4 Ligase buffer(商品名:10X Buffer for T4 DNA ligase with 10mM ATP; New England BioLabs社) を250μl、0.5wt.%BSA含有水溶液を12.5μl、1MのMgClを3μl、T4 DNAリガーゼ(400U/μl)(商品名 T4 DNA Ligase、New England BioLabs社)を12.5μl加え、更に滅菌水を加えて全容量を2500μlとし、16℃で一晩混合した。その後、4℃、遠心力×16000g、10分間で遠心分離操作を行い、ペレットを得た。
【0113】
10倍濃度で調整したNEBuffer#1(商品名 NEBuffer 1, New England BioLabs社)を10μlと、滅菌水89μl、エキソヌクレアーゼIII(Exonuclease III(100U/μl)、New England Biolabs社)1μlを、ペレットに加え、100μlのサンプル溶液とした。これを、37℃で5分間保持し、続いて、溶液の1/200容量のプロテアーゼ(商品名 Proteinase K (Recombinant) Solution、Nacalai Tesque社)を加え、55℃で3時間保持してサンプル溶液中のタンパク質を消化した(脱架橋)。
【0114】
DNA溶液を、更に精製した後、超微量分光光度計(商品名Qubit 3.0 Fluorometer、Invitrogen/ Life Technologies社)を用いてDNA量を確認した。
【0115】
そして、DNA量が10ngとなるようにDNA溶液をチューブに分注し、2倍濃度のTD buffer25μl(商品名 Tagment DNA Buffer (Nextera(登録商標) DNA Library Prep Kit、Illumina社)とトランスポソーム(Tagment DNA Enzyme 1 (TDE1)、Illumina社)0.5μlとを加えた後、溶液総量が50μlとなるまで滅菌水を追加した。これを55℃で10分間保持し、タグメンテーションを行った。反応後のサンプルを精製し、50μlのDNA溶液(Buffer EB溶液)とした。
【0116】
次に、ストレプトアビジンビーズ分散液(商品名DynabeadsTM(登録商標) M-280 Streptavidin、Invitrogen/ Thermo Fisher Scientific社)10μlを、500μlの1xBind and Wash bufferで2回洗浄した後、50μlの2xBind and Wash Bufferで懸濁した。この50μlのビーズ分散液を、上記で精製したDNA溶液とチューブ内で混合し、30分間、室温で保持した後、磁気スタンド(1.5ml用)に静置して、上清を除去し、ストレプトアビジンビーズを回収した。そして、回収したビーズを、0.1vol%の非イオン性界面活性剤溶液(商品名Triton X-100、Nacalai Tesque社)を含む500μlの1x Bind and wash bufferで1度洗浄してから、500μlのBuffer EBで更に洗浄した後、10μlのBuffer EBで懸濁した。
【0117】
10μlに、滅菌水10μl、10μMのフォワード側プライマー(N704、Nextera(登録商標) XT Index Kit v.2, Illumina社)2.5μl、リバース側プライマー(S521、商品名:Nextera(登録商標) XT Index Kit v.2, Illumina社)2.5μl、DNAポリメラーゼ(KAPA HiFi HotStart ReadyMix (2X)、Kapa Biosystems社)25μlを加えてから、PCRを行った。PCRの条件は、72℃5分間の予備加熱、98℃で3分間の初期変性後、98℃で10秒間の熱変性、63℃で30秒間のアニーリング反応、72℃で1分間の伸長反応を11サイクル行い、72℃で1分間の最終伸長反応を行った後、増幅産物を4℃に冷却した。
【0118】
なお、このPCRに先立って、DNA溶液の一部を用いて、リアルタイムPCRを実施し、そのプラトーから最適サイクル数は決定した。
【0119】
得られたPCR産物から、目的のDNAサイズ200-1000bpのDNA断片を回収するため、適宜のサイズセレクションを実施し、最終的に、150bp~1300bpのDNA断片が含まれたDNAライブラリ 10μl(DNA/Buffer EB溶液)を得た。
【0120】
この得られたDNAライブラリを用い、超並列DNA配列解析装置(装置名 NextSeq 500、Illumina社)を使用して、DNA断片(DNAペア)の両端の塩基配列を解析した。解析する塩基配列の長さは、36塩基とした。得られた塩基配列に対し、両端の配列をそれぞれ独立に参照用ヒトゲノム配列にマッピングすることで、DNA配列の相互作用が解析された。得られた結果を図10に示す。
【0121】
図10は、一例として解析された、ゲノム上におけるCTCFの局在およびCTCFが存在するクロマチン領域間の相互作用を示す。最上部にはゲノム上の位置が数字で示されており、上列には、既知のCUT&RUN法(Skene et al., 2018)を別で実施して得た結果において解析されたCTCFの局在を示した分布図を、参照情報として示している。中列は、Jaspar Transcription Factor Databaseにより特定されたゲノム上に存在するCTCF結合サイトを示しており、ゲノム上の配列の方向性に応じ、一方の方向のものに下線を付して示している。下列は、本願実施例3により得られたCTCFの局在を示した分布図と共に、CTCFの局在箇所の下に、本願実施例3により得られたCTCFが存在するクロマチン領域間の相互作用をループ状で示す。下列の更に下方には、第7染色体中の遺伝子(FBXL18、MIR589、ACTB)の位置が表示されており、遺伝子名称の右側にマーキングされた範囲が対応する遺伝子の存在位置を示している。
図11は、図10の部分的な拡大図である。
図12は、図10の部分的な拡大図である。
図1
図2
図3
図4
図5
図6
図7
図8
図9
図10
図11
図12