Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7597783B2 - バイサルファイトを使用しない全ゲノムメチル化解析 - Google Patents
[go: Go Back, main page]

JP7597783B2 - バイサルファイトを使用しない全ゲノムメチル化解析 - Google Patents

バイサルファイトを使用しない全ゲノムメチル化解析 Download PDF

Info

Publication number
JP7597783B2
JP7597783B2 JP2022500944A JP2022500944A JP7597783B2 JP 7597783 B2 JP7597783 B2 JP 7597783B2 JP 2022500944 A JP2022500944 A JP 2022500944A JP 2022500944 A JP2022500944 A JP 2022500944A JP 7597783 B2 JP7597783 B2 JP 7597783B2
Authority
JP
Japan
Prior art keywords
dna
taps
target dna
5hmc
5cac
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2022500944A
Other languages
English (en)
Other versions
JP2022540453A5 (ja
JP2022540453A (ja
Inventor
ソン,チュンシャオ
チェン,ジンフェイ
シエイカ-ジエリンスカ,パウリナ
リウ,イビン
Original Assignee
ルドウイグ インスティテュート フォー キャンサー リサーチ エルティーディー
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ルドウイグ インスティテュート フォー キャンサー リサーチ エルティーディー filed Critical ルドウイグ インスティテュート フォー キャンサー リサーチ エルティーディー
Publication of JP2022540453A publication Critical patent/JP2022540453A/ja
Publication of JP2022540453A5 publication Critical patent/JP2022540453A5/ja
Priority to JP2024207194A priority Critical patent/JP2025029007A/ja
Application granted granted Critical
Publication of JP7597783B2 publication Critical patent/JP7597783B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813Hybridisation assays
    • C12Q1/6827Hybridisation assays for detection of mutation or polymorphism
    • C12Q1/683Hybridisation assays for detection of mutation or polymorphism involving restriction enzymes, e.g. restriction fragment length polymorphism [RFLP]
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12PFERMENTATION OR ENZYME-USING PROCESSES TO SYNTHESISE A DESIRED CHEMICAL COMPOUND OR COMPOSITION OR TO SEPARATE OPTICAL ISOMERS FROM A RACEMIC MIXTURE
    • C12P19/00Preparation of compounds containing saccharide radicals
    • C12P19/26Preparation of nitrogen-containing carbohydrates
    • C12P19/28N-glycosides
    • C12P19/30Nucleotides
    • C12P19/34Polynucleotides, e.g. nucleic acids, oligoribonucleotides
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6806Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813Hybridisation assays
    • C12Q1/6827Hybridisation assays for detection of mutation or polymorphism
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2521/00Reaction characterised by the enzymatic activity
    • C12Q2521/30Phosphoric diester hydrolysing, i.e. nuclease
    • C12Q2521/301Endonuclease
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2521/00Reaction characterised by the enzymatic activity
    • C12Q2521/50Other enzymatic activities
    • C12Q2521/531Glycosylase
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2523/00Reactions characterised by treatment of reaction samples
    • C12Q2523/10Characterised by chemical treatment
    • C12Q2523/115Characterised by chemical treatment oxidising agents
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2525/00Reactions involving modified oligonucleotides, nucleic acids, or nucleotides
    • C12Q2525/10Modifications characterised by
    • C12Q2525/117Modifications characterised by incorporating modified base
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2525/00Reactions involving modified oligonucleotides, nucleic acids, or nucleotides
    • C12Q2525/10Modifications characterised by
    • C12Q2525/191Modifications characterised by incorporating an adaptor
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/154Methylation markers

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Molecular Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Genetics & Genomics (AREA)
  • Microbiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biotechnology (AREA)
  • General Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Immunology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • General Chemical & Material Sciences (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Saccharide Compounds (AREA)
  • Enzymes And Modification Thereof (AREA)

Description

発明の詳細な説明
〔技術分野〕
本開示は、全ゲノムにおいて、5-メチルシトシン、5-ヒドロキシメチルシトシン、5-カルボキシルシトシン及び/または5-ホルミルシトシンの位置を特定する方法を提供する。
〔背景技術〕
5-メチルシトシン(5mC)及び5-ヒドロキシメチルシトシン(5hmC)は、哺乳動物ゲノムに見られる2つの主要エピジェネティックマークである。5hmCは、10-11転座(TET)ファミリーのジオキシゲナーゼによって、5mCから生成される。Tetはさらに、5hmCを酸化して、5-ホルミルシトシン(5fC)及び5-カルボキシルシトシン(5caC)にすることができ、哺乳動物ゲノムにおいて、5caCは、存在量が5mC及び5hmCよりもかなり少ない(5hmCの量の10分の1~100分の1)。5mC及び5hmCは連動して、遺伝子調節から正常な発生に至る広範な生物学的プロセスで重要な役割を果たす。DNAの異常なメチル化及びヒドロキシメチル化は、様々な疾患と関連付けられており、がんの特徴として広く認められている。したがって、DNA配列における5mC及び5hmCの判定は、基礎研究に重要であるのみならず、診断及び治療を含む臨床用途にも有用である。
5fC及び5caCは、5mCの最終酸化誘導体の2つであり、塩基除去修復経路におけるチミンDNAグリコシラーゼ(TDG)によって、無修飾シトシンに変換できる。すなわち、5fC及び5caCは、胚発生で重要な役割を果たす能動的脱メチル化プロセスにおいて、2つの重要な鍵となる中間体である。5fC及び5caCは、これらの状況において確認されており、5mCの脱メチル化がほぼ完了したことの指標として機能し得る。5fC及び5caCは、特異的タンパク質との結合、ならびにRNAポリメラーゼIIの速度及び特異性に対する作用といった追加の機能も果たし得る。
DNAのメチル化及びヒドロキシメチル化の解析方法のうち、現行のゴールドスタンダードであるとともに、最も広く使われている方法は、バイサルファイトシーケンシング(BS)、ならびに、Tetアシストバイサルファイトシーケンシング(TAB-Seq)及び酸化バイサルファイトファイトシーケンシング(oxBS)といったその派生方法である。同様に、バイサルファイトシーケンシングは、全ゲノムDNAのメチル化のアッセイ方法のうち、最も確立された方法である。これらの方法のいずれも、バイサルファイト処理を用いて、非メチル化シトシンをウラシルに変換する一方で、5mC及び/または5hmCは、インタクトなままにする。バイサルファイト処理したDNAのPCR増幅(ウラシルをチミンとして読み取る)を通じて、各シトシンの修飾情報を1塩基の分解能で得ることができる(CをTに変換することで、非メチル化シトシンの位置が得られる)。しかしながら、バイサルファイトシーケンシングには、主な欠点が少なくとも2つある。第1に、バイサルファイト処理は、所要の酸性条件及び温熱条件下で、脱プリン化が原因で、DNAの90%超が分解される過酷な化学反応である点である。この分解により、その用途は、循環セルフリーDNA及びシングルセルシーケンシングを含む臨床試料のように、インプットの少ない試料に厳しく制限される。第2に、バイサルファイトシーケンシングは、無修飾シトシンのチミンへの完全な変換に依存する点である。無修飾シトシンは、ヒトゲノムにおいて、全シトシンのおよそ95%を占める。これらのすべての位置がチミンに変換されると、配列複雑度が大幅に低下し、シーケンシングクオリティが下がり、マッピング率が低くなり、ゲノムカバレッジが不均一となり、シーケンシングコストが上昇する。バイサルファイトシーケンシング法は、無修飾シトシンのチミンへの変換が不完全なことにより、5mC及び5hmCの偽検出も発生しやすい。
全ゲノム塩基の分解能のメチロームシーケンシングにより、DNAメチル化は最も網羅的に解析可能となるが、高額なシーケンシングコストにより、用途が制限される場合が多い。シーケンシングコストを削減するために、Reduced Representation Bisulfite Sequencing(RRBS)を用いてもよい。しかしながら、RRBSは、所定の配列コンテクストにおけるごく一部のCpG部位しか対象としないので、メチル化の網羅的な状況は得られない。したがって、低コストで、mCpGのカバレッジを向上させるアプローチに対するニーズが存在する。
〔発明の概要〕
本開示は、全ゲノム塩基の分解能のメチローム解析方法を提供する。実施形態では、その方法は、DNA試料における5-メチルシトシン、5-ヒドロキシメチルシトシン、5-カルボキシルシトシン及び/または5-ホルミルシトシンのうちの1つ以上の位置を特定する。実施形態では、そのDNA試料は、全ゲノムを含む。本明細書に記載されている方法は、DNAのメチル化及びヒドロキシメチル化の解析であって、無修飾シトシンに影響を及ぼすことなく、塩基単位の分解能で、修飾シトシンを検出する穏和な反応を伴う解析を実現する。本発明で提供するのは、TETによる酸化と、ボラン誘導体(例えば、ピリジンボラン及び2-ピコリンボラン(pic-BH))による還元を組み合わせることによって、5mC及び/または5hmCを特定するための改良型の方法であり、この方法は、本明細書では、TAPS(TETアシストピリジンボランシーケンシング)(表1)という。TAPSは、無修飾シトシンに影響を及ぼすことなく、高い感度及び特異性で、修飾を直接検出するとともに、本明細書に記載されているように、改変して、他のシトシン修飾を検出することができる。TAPSは、非破壊的であり、DNAを最長で10kbに保つものである。TAPSでは、バイサルファイトシーケンシングと比べて、マッピング率が高く、カバレッジが均一で、シーケンシングコストが低く、クオリティの上昇が可能となり、網羅性の向上した、より安価なメチローム解析が行われる。本明細書に記載されているように、5fC及び/または5caCを特定するには、酸化工程を用いない、この方法の変形形態を使用する。
全ゲノムTAPS(wgTAPS)では、全ゲノムバイサルファイトシーケンシング(WGBS)と比べて、シーケンシングコストが半減する。しかしながら、特に、5mC及び5hmCが、哺乳動物ゲノムにおける全シトシン残基の約4%を占めるに過ぎないため、ショートリード全ゲノムシーケンシングによって生成されるリードの65%~80%前後が、メチル化CpG部位(mCpG)を全く含まないことを考慮すると、全ゲノムシーケンシングのコストは、多くのプロジェクトにとって、依然として法外である。シーケンシングコストを削減するために、バイサルファイト処理の前に、CpGに富む領域を制限酵素によって濃縮するReduced Representation Bisulfite Sequencing(RRBS)が、広く用いられている方法である。しかしながら、RRBSは、所定の配列コンテクストにおけるごく一部のCpG部位しかカバーしないので、メチル化の網羅的な状況は得られない。したがって、低コストで、mCpGのカバレッジを向上させるアプローチに対するニーズが存在する。一実施形態では、本開示は、費用効果が良いゲノムワイドなメチル化解析を1CpGの分解能で行うための新たな方策をもたらす、TAPS法の改変型(本明細書では、エンドヌクレアーゼ濃縮TAPS(eeTAPS)という)を提供する。別の実施形態では、本開示は、5mC、5hmC、5fC及び/または5caCを検出するのに使用できる、TAPS法の変形形態の改変型(例えば、TAPSβ及びCAPS)を提供する。
すなわち、本開示は、全ゲノムにおいて、5-メチルシトシン、5-ヒドロキシメチルシトシン、5-カルボキシルシトシン及び/または5-ホルミルシトシンのうちの1つ以上を特定するための、低コストな方法を提供する。本明細書に記載されている方法は、例えばTETアシストピリジンボラン処理によって、修飾シトシン(5mC、5hmC、5fC、5caC)をジヒドロウラシル(DHU)に変換してから、そのDHU部位をエンドヌクレアーゼによって切断して、DNA断片を作製した後、その断片をシーケンシングライブラリーにすることをベースとする。断片化されなかったゲノムDNAは、直接シーケンシングすることはできず、DHUに変換される修飾シトシンが存在する場合のみ、ゲノムDNAは切断されて、DNA断片となった後、その断片の各末端(修飾シトシンの位置を示す)により、シーケンシングできる。その切断された断片をシーケンシングすることによって、元のDNA試料におけるメチル化シトシンを塩基単位の分解能で特定できる。さらに、高メチル化シトシンは、低メチル化部位よりも、切断される傾向が高いので、この方法を用いて、DNAメチル化の半定量を行うことができる。
一態様では、本開示は、DNA試料における5mCまたは5hmCの特定方法であって、
a.5mC及び/または5hmCを有する標的DNAを含むDNA試料を準備する工程と、
b.そのDNAを改変する工程(以下の工程を含む)と、
i.そのDNA試料中の5mC及び5hmCを5-カルボキシルシトシン(5caC)及び/または5fCに変換する工程
ii.その5caC及び/または5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5mCまたは5hmCのいずれかの位置が得られる方法を提供する。
別の態様では、本開示は、DNA試料における5-メチルシトシン(5mC)の特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料を改変する工程(以下の工程を含む)と、
i.そのDNA試料中の5-ヒドロキシメチルシトシン(5hmC)に、ブロック基を付加する工程
ii.そのDNA試料中の5mCを5-カルボキシルシトシン(5caC)及び/または5-ホルミルシトシン(5fC)に変換する工程
iii.その5caC及び/または5fCをジヒドロウラシル(DHU)に変換して、改変標的DNAを含む改変DNA試料をもたらす工程
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5mCの位置が得られる方法を提供する。
別の態様では、本開示は、DNA試料における5mCまたは5hmCの特定方法であって、
a.5mC及び/または5hmCを有する標的DNAを含むDNA試料を準備する工程と、
b.そのDNAを改変する工程(以下の工程を含む)と、
i.そのDNA試料中の5mC及び5hmCを5-カルボキシルシトシン(5caC)及び/または5fCに変換する工程
ii.その5caC及び/または5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5mCまたは5hmCのいずれかの位置が得られる方法を提供する。
別の態様では、本開示は、DNA試料において、5mCを同定するとともに、5hmCを同定する方法であって、
a.そのDNA試料中の5mCを特定すること(以下の工程を含む)と、
i.5mC及び/または5hmCを有する標的DNAを含む第1のDNA試料を準備する工程
ii.その第1の試料中のDNAを改変する工程(以下の工程を含む)
1.その第1のDNA試料中の5-ヒドロキシメチルシトシン(5hmC)に、ブロック基を付加する工程
2.その第1のDNA試料中の5mCを5caC及び/または5fCに変換する工程
3.その5caC及び/または5fCをDHUに変換して、改変標的DNAを含む第1の改変DNA試料をもたらす工程
iii.その改変標的DNAを切断する工程
iv.その切断された改変標的DNAに、アダプターDNA分子を付加する工程
v.その改変標的DNAの配列を検出する工程であって、切断部位の存在により、その標的DNAにおける5mCの位置が得られる工程
b.そのDNA試料中の5mCまたは5hmCを特定すること(以下の工程を含む)と、
i.5mC及び/または5hmCを有する標的DNAを含む第2のDNA試料を準備する工程
ii.その第2の試料中のDNAを改変する工程(以下の工程を含む)
1.その第2のDNA試料中の5mC及び5hmCを5caC及び/または5fCに変換する工程
2.その5caC及び/または5fCをDHUに変換して、改変標的DNAを含む第2の改変DNA試料をもたらす工程
iii.その改変標的DNAを切断する工程
iv.その切断された改変標的DNAに、アダプターDNA分子を付加する工程
v.その第2の試料に由来する改変標的DNAの配列を検出する工程であって、切断部位の存在により、その標的DNAにおける5mCまたは5hmCのいずれかの位置が得られる工程
c.工程(a)及び(b)の結果を比較することであって、工程(b)では存在するが、工程(a)では存在しない切断部位により、その標的DNAにおける5hmCの位置が得られることと、
を含む方法を提供する。
実施形態では、そのDNA試料中の5hmCに付加するブロック基は、糖である。実施形態では、その糖は、天然の糖または修飾糖、例えば、グルコースまたは修飾グルコースである。実施形態では、ブロック基は、グルコシルトランスフェラーゼ酵素、例えば、T4バクテリオファージβ-グルコシルトランスフェラーゼ(βGT)及びT4バクテリオファージα-グルコシルトランスフェラーゼ(αGT)、ならびにそれらの誘導体及び類似体の存在下で、DNA試料を、糖に連結されたUDP、例えば、UDP-グルコース、または修飾グルコースに連結されたUDPと接触させることによって、5hmCに付加する。
実施形態では、DNA試料中の5mCを5caC及び/または5fCに変換する工程、ならびにDNA試料中の5mC及び5hmCを5caC及び/または5fCに変換する工程はそれぞれ、そのDNA試料を10-11転座(TET)酵素と接触させることを含む。さらなる実施形態では、そのTET酵素は、ヒトTET1、ヒトTET2、ヒトTET3、マウスTet1、マウスTet2、マウスTet3、Naegleria TET(NgTET)、Coprinopsis cinerea(CcTET)、及びこれらの誘導体または類似体のうちの1つ以上である。実施形態では、そのTET酵素は、NgTETまたはその誘導体である。実施形態では、そのTET酵素は、マウスmTet1(mTet1CD)またはその誘導体である。別の実施形態では、そのTET酵素は、ヒトTET2(hTET2)またはその誘導体である。
別の態様では、本開示は、DNA試料中の5caCまたは5fCの特定方法であって、
a.5caC及び/または5fCを有する標的DNAを含むDNA試料を準備する工程と、
b.その5caC及び5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCまたは5fCのいずれかの位置が得られる方法を提供する。
別の態様では、本開示は、DNA試料中の5caCの特定方法であって、
a.5caCを有する標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5fCに、ブロック基を付加する工程と、
c.その5caCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
d.その改変標的DNAを切断する工程と、
e.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
f.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCの位置が得られる方法を提供する。
実施形態では、そのDNA試料中の5fCに、ブロック基を付加する工程は、そのDNAをアルデヒド反応性化合物(例えば、ヒドロキシルアミン誘導体(O-エチルヒドロキシルアミンなど)、ヒドラジン誘導体及びヒドラジド誘導体を含む)と接触させることを含む。
別の態様では、本開示は、DNA試料中の5fCの特定方法であって、
a.5fCを有する標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5caCに、ブロック基を付加する工程と、
c.その5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
d.その改変標的DNAを切断する工程と、
e.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
f.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5fCの位置が得られる方法を提供する。
実施形態では、そのDNA試料中の5caCに、ブロック基を付加する工程は、そのDNA試料をカルボン酸誘導体化試薬(例えば、1-エチル-3-(3-ジメチルアミノプロピル)カルボジイミド(EDC)を含む)、及び(ii)アミン(エチルアミンなど)、ヒドラジンまたはヒドロキシルアミン化合物と接触させることを含む。
別の態様では、本開示は、DNA試料中の5hmCの特定方法であって、
a.5hmCを有する標的DNAを含むDNA試料を準備することと、
b.その試料中のDNAを改変すること(以下の工程を含む)と、
i.そのDNA試料中の5hmCを5caC及び/または5fCに変換する工程
ii.その5caC及び/または5fCをジヒドロウラシル(DHU)に変換して、改変標的DNAを含む改変DNA試料をもたらす工程
c.その改変標的DNAを切断することと、
d.その切断された改変標的DNAに、アダプターDNA分子を付加することと、
e.その改変標的DNAの配列を検出することと、
を含み、
切断部位の存在により、その標的DNAにおける5hmCの位置が得られる方法を提供する。
実施形態では、5hmCを5caC及び/または5fCに変換する工程は、そのDNA試料を酸化剤と接触させることを含む。さらなる実施形態では、その酸化剤は、過ルテニウム酸カリウム、Cu(II)/TEMPO、ルテニウム酸カリウムまたは酸化マンガンである。
実施形態では、そのDNA試料は、ゲノムDNAを含む。実施形態では、そのDNA試料は、ピコグラム量のDNAを含む。実施形態では、そのDNA試料は、DNAを約1pg~約900pg、DNAを約1pg~約500pg、DNAを約1pg~約100pg、DNAを約1pg~約50pg、DNAを約1~約10pg、DNAを約200pg未満、約100pg未満、DNAを約50pg未満、DNAを約20pg未満、及びDNAを約5pg未満含む。別の実施形態では、そのDNA試料は、ナノグラム量のDNAを含む。実施形態では、そのDNA試料は、DNAを約1~約500ng、DNAを約1~約200ng、DNAを約1~約100ng、DNAを約1~約50ng、DNAを約1ng~約10ng、DNAを約1ng~約5ng、DNAを約100ng未満、DNAを約50ng未満、DNAを約5ng、またはDNAを約2ng未満含む。実施形態では、そのDNA試料は、循環セルフリーDNA(cfDNA)を含む。実施形態では、そのDNA試料は、マイクログラム量のDNAを含む。
実施形態では、5caC及び/または5fCをDHUに変換する工程は、そのDNA試料を還元剤(例えば、ピリジンボラン、2-ピコリンボラン(pic-BH)、tert-ブチルアミンボラン、ボラン、水素化ホウ素ナトリウム、シアノ水素化ホウ素ナトリウム及びトリアセトキシ水素化ホウ素ナトリウムを含む)と接触させることを含む。好ましい実施形態では、その還元剤は、pic-BH及び/またはピリジンボランである。
実施形態では、その改変標的DNAを切断する工程は、改変標的DNAを含む改変DNA試料と、DHUの存在により、改変標的DNAを特異的に切断する1つ以上のエンドヌクレアーゼとを接触させることによって、DHUを含む改変標的DNAを特異的に切断することを含む。実施形態では、そのエンドヌクレアーゼは、DNA N-グリコシラーゼ活性及びAPリアーゼ活性を有する2官能性DNAエンドヌクレアーゼであり、例えば、TmaエンドヌクレアーゼIII、エンドヌクレアーゼVIII、ホルムアミドピリミジンDNAグリコシラーゼ(Fpg)及び/またはhNEIL1が挙げられる。実施形態では、その改変標的DNAは、ウラシル特異的除去試薬(USER)を用いて切断する。そのUSER酵素は、ウラシルDNAグリコシラーゼ(UDG)及びエンドヌクレアーゼVIIIを組み合わせたものを含む。改変標的DNAを切断するために、単独でまたは組み合わせて使用できる他の酵素は、アプリン酸/アピリミジン酸エンドヌクレアーゼ1(APE1)、UDG、エンドヌクレアーゼIII、TmaエンドヌクレアーゼIII、TthエンドヌクレアーゼIV、エンドヌクレアーゼV、エンドヌクレアーゼVIII、Fpg及び/またはhNEIL1である。
実施形態では、上記の方法は、その切断された改変標的DNAのサイズ選別を行う工程をさらに含む。実施形態では、上記の方法は、その改変標的DNAのコピー数を増幅する工程をさらに含む。実施形態では、この増幅工程は、その改変標的DNAの配列を検出する工程の前に行う。その改変標的DNAのコピー数を増幅する工程は、ポリメラーゼ連鎖反応(PCR)、プライマー伸長法及び/またはクローニングを行うことによって行ってよい。
実施形態では、その改変標的DNAの配列を決定する工程は、チェーンターミネーションシーケンシング、マイクロアレイ、ハイスループットシーケンシング及び制限酵素解析を含む。実施形態では、その改変標的DNAの配列を検出する工程は、次世代シーケンシング法を含む。
一態様では、本開示は、改変標的DNAの切断方法であって、1つ以上のDHU(すなわちDHU残基)を含む改変標的DNAと、その1つ以上のDHUの位置またはその1つ以上のDHUと隣接する位置で、その改変標的DNAを切断する1つ以上のエンドヌクレアーゼとを接触させること含む方法を提供する。その1つ以上のエンドヌクレアーゼは、例えば、本明細書に記載されているいずれかのタイプのエンドヌクレアーゼまたはそれらを組み合わせたものであってよい。実施形態では、その改変標的DNAにおける1つ以上のDHUは、例えば、本明細書に記載されている方法によって、5mC及び/または5hmCから派生させたものである。実施形態では、その改変標的DNAにおける1つ以上のDHUは、例えば、本明細書に記載されている方法によって、5caC及び/または5fCから派生させたものである。
実施形態では、その方法は、接触工程の前に、その改変標的DNAをもたらすために、5mC及び/または5hmCを含む標的DNAを改変することをさらに含み、その改変は、(i)例えば、本明細書に記載されている方法によって、その標的DNAにおける5mC及び5hmCを5-カルボキシルシトシン(5caC)及び/または5-ホルミルシトシン(5fC)に変換する工程と、(ii)例えば、本明細書に記載されている方法によって、その5caC及び/または5fCをジヒドロウラシル(DHU)に変換する工程を含む。
実施形態では、その方法は、接触工程の前に、改変標的DNAをもたらすために、5mC及び5hmCを含む標的DNAを改変することをさらに含み、その改変は、(i)例えば、本明細書に記載されている方法によって、その標的DNAにおける5hmCに、ブロック基を付加する工程と、(ii)例えば、本明細書に記載されている方法によって、その標的DNAにおける5mCを5-カルボキシルシトシン(5caC)及び/または5-ホルミルシトシン(5fC)に変換する工程と、(iii)例えば、本明細書に記載されている方法によって、その5caC及び/または5fCをジヒドロウラシル(DHU)に変換する工程を含む。
実施形態では、その方法は、接触工程の前に、改変標的DNAをもたらすために、5caC及び/または5fCを含む標的DNAを改変することをさらに含み、その改変は、例えば、本明細書に記載されている方法によって、その5caC及び/または5fCをジヒドロウラシル(DHU)に変換する工程を含む。
実施形態では、その方法は、接触工程の前に、改変標的DNAをもたらすために、5caC及び/または5fCを含む標的DNAを改変することをさらに含み、その改変は、(i)例えば、本明細書に記載されている方法によって、その標的DNAにおける5fCにブロック基を付加する工程と、(ii)例えば、本明細書に記載されている方法によって、その5caCをジヒドロウラシル(DHU)に変換する工程を含む。
実施形態では、その方法は、接触工程の前に、改変標的DNAをもたらすために、5caC及び/または5fCを含む標的DNAを改変することをさらに含み、その改変は、(i)例えば、本明細書に記載されている方法によって、その標的DNAにおける5caCに、ブロック基を付加する工程と、(ii)例えば、本明細書に記載されている方法によって、その5fCをジヒドロウラシル(DHU)に変換する工程を含む。
実施形態では、その方法は、接触工程の前に、改変標的DNAをもたらすために、5mC及び/または5hmCを含む標的DNAを改変することをさらに含み、その改変は、(i)例えば、本明細書に記載されている方法によって、その標的DNAにおける5mCを5-カルボキシルシトシン(5caC)及び/または5-ホルミルシトシン(5fC)に変換する工程と、(ii)例えば、本明細書に記載されている方法によって、その5caC及び/または5fCをジヒドロウラシル(DHU)に変換する工程を含む。
〔図面の簡単な説明〕
〔図1〕ボラン含有化合物のスクリーニングである。11merのオリゴヌクレオチド(「オリゴ」)における5caCのDHUへの変換について、ボラン含有化合物をスクリーニングし、変換率は、MALDIによって推定した。2-ピコリンボラン(pic-ボラン)、ボランピリジン及びtert-ブチルアミンボランは、5caCをDHUに完全に変換できたが、エチレンジアミンボラン及びジメチルアミンボランでは、変換率は30%前後であった。ジクロロヘキシルアミンボラン、モルホリンボラン、4-メチルモルホリンボラン及びトリメチルアミンボランでは、検出可能な生成物は、測定されなかった(不検出)。水素化ホウ素ナトリウム及びトリ(アセトキシ)水素化ホウ素ナトリウムのような他の還元剤は、酸性培地中で速やかに分解されて、変換が不完全となった。シアノ水素化ホウ素ナトリウムは、酸性条件下では、シアン化水素が形成される可能性があることから、使用しなかった。完全な変換、低い毒性及び高い安定性から、Pic-ボラン及びピリジンボランを選択した。
〔図2〕DNAオリゴにおけるPic-ボラン反応である。(A)5caCを含む11merのモデルDNAをpic-ボランで処理したものをMALDIによって特徴付けたものである。各グラフの上に、質量(m/z)の計算値が示されており、ピークの左に、質量の観察値が示されている。(B)dC及び各種シトシン誘導体の変換率をHPLC-MS/MSによって定量した。データは、3つのレプリケートの平均±標準偏差としてされている。
〔図3A〕単一ヌクレオシドのpic-ボラン反応である。H NMRの結果は、2’-デオキシ-5,6-ジヒドロウリジンに関する過去の報告に従ったものである(I.Aparici-Espert et al.,J.Org.Chem.81,4031-4038(2016))。単一ヌクレオシドのpic-ボラン反応生成物のH NMR(MeOH-d、400MHz)チャートである。δ ppm: 6.28 (t, 1H, J = 7 Hz), 4.30 (m, 1H), 3.81 (m, 1H), 3.63 (m, 2H), 3.46 (m, 2H), 2.65 (t, 2H, J = 6 Hz), 2.20 (m, 1H), 2.03 (m, 1H).
〔図3B〕単一ヌクレオシドのpic-ボラン反応である。13C NMRの結果は、2’-デオキシ-5,6-ジヒドロウリジンに関する過去の報告に従ったものである(I.Aparici-Espert et al.,J.Org.Chem.81,4031-4038(2016))。単一ヌクレオシドのpic-ボラン反応生成物の13C NMR (MeOH-d, 400 MHz)チャートである。δ ppm: 171.56 (CO), 153.54 (CO), 85.97 (CH), 83.86 (CH), 70.99 (CH), 61.92 (CH), 36.04 (CH), 35.46 (CH), 30.49 (CH).
〔図4A〕5caCのDHUへの、ボランによる変換と、5caCのDHUへのボラン反応の機序案を示す図である。
〔図4B〕5fCのDHUへの、ボランによる変換と、5fCのDHUへのボラン反応の機序案を示す図である。
〔図5〕(A)TAPS法が、5mC及び5hmCの両方をDHUに変換し、そのDHUが、複製時にはチミンとして機能することを示す図である。(B)TAPS、TAPSβ及びCAPSの方法の概説である。
〔図6-1〕5fC及び5caCを含むモデルDNAオリゴをpic-ボランによって、5fC及び5caCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5fC及び5caCは、pic-BHによってジヒドロウラシル(DHU)に変換される。5fCは、O-エチルヒドロキシルアミン(EtONH)のようなヒドロキシルアミン誘導体によってブロックしたが、その誘導体は、オキシムとなり、pic-ボランによる変換を阻止することになる。5caCは、EDC結合を介して、エチルアミンによってブロックして、pic-ボランによる変換をブロックするアミドに変換した。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。
〔図6-2〕5fC及び5caCを含むモデルDNAオリゴをpic-ボランによって、5fC及び5caCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5fC及び5caCは、pic-BHによってジヒドロウラシル(DHU)に変換される。5fCは、O-エチルヒドロキシルアミン(EtONH)のようなヒドロキシルアミン誘導体によってブロックしたが、その誘導体は、オキシムとなり、pic-ボランによる変換を阻止することになる。5caCは、EDC結合を介して、エチルアミンによってブロックして、pic-ボランによる変換をブロックするアミドに変換した。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。
〔図7-1〕5mC及び5hmCを含むモデルDNAオリゴをKRuO及びpic-ボランによって、5hmCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5hmCは、βGTによってグルコースでブロックして、5gmCに変換できた。5mC、5hmC及び5gmCは、pic-ボランによって変換できなかった。5hmCは、KRuOによって酸化させて、5fCにしてから、pic-ボランによってDHUに変換できた。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。
〔図7-2〕5mC及び5hmCを含むモデルDNAオリゴをKRuO及びpic-ボランによって、5hmCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5hmCは、βGTによってグルコースでブロックして、5gmCに変換できた。5mC、5hmC及び5gmCは、pic-ボランによって変換できなかった。5hmCは、KRuOによって酸化させて、5fCにしてから、pic-ボランによってDHUに変換できた。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。
〔図8〕制限酵素消化によって、TAPSが、5mCをTに効果的に変換できることが示された。(A)TAPSによる配列変化を確認するための制限酵素消化アッセイの図である。(B)TAPSによる、CからTへの変換を確認するためのTaqαI消化試験である。TAPSは、TaqαI制限部位を有するとともに、完全にメチル化された5つのCpG部位(5mC)を含む222bpのモデルDNAと、その非メチル化コントロール(C)で行った。5mC、C及びC TAPSに示されているように、PCR増幅した222bpのモデルDNAは、TaqαIによって、約160bp及び約60bpの断片に切断できる。5mC-TAPSのレーンに示されているように、メチル化DNAでTAPSを行った後、T(mC)GA配列は、TTGAに変換され、TaqαI消化によっては切断されなくなる。
〔図9〕222bpのモデルDNA及びmESC gDNAでのTAPSである。(A)TAPSの前(5mC、C)及びTAPSの後(5mC TAPS、C TAPS)における、完全にメチル化された5つのCpG部位を含む222bpのモデルDNAと、その非メチル化コントロールのサンガーシーケンシングの結果である。5mCのみが、TAPS法によってTに変換される。(B)NgTET1による酸化後及びpic-ボランによる還元後のmESC gDNAコントロールにおける相対的な改変レベルの、HPLC-MS/MSによる定量結果である。データは、3つのレプリケートの平均±標準偏差としてされている。
〔図10A〕TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。氷浴で冷却する前(の222bpの非メチル化DNA、222bpのメチル化DNA及びmESC gDNAのアガロースゲル画像である。TAPS後、検出可能なDNA分解は、観察されず、DNAは、二本鎖を保ち、冷却せずに見ることができた。バイサルファイト変換では、分解が起こり、DNAは、一本鎖となり、氷上で冷やした後のみに見ることができた。
〔図10B〕TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。氷浴で冷却した後(B)の222bpの非メチル化DNA、222bpのメチル化DNA及びmESC gDNAのアガロースゲル画像である。TAPS後、検出可能なDNA分解は、観察されず、DNAは、二本鎖を保ち、冷却せずに見ることができた。バイサルファイト変換では、分解が起こり、DNAは、一本鎖となり、氷上で冷やした後のみに見ることができた。
〔図10C〕TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。TAPS及びバイサルファイトで処理した、様々な断片長のmESC gDNAであって、氷上で冷やす前(左パネル)及び氷上で冷やした後(右パネル)のmESC gDNAのアガロースゲル画像である。TAPS後のDNAは、二本鎖を保ち、ゲル上で、直接見ることができた。バイサルファイト処理の方が、試料の損傷及び断片化が大きく、DNAは、一本鎖となり、氷上で冷やした後のみに見ることができた。図15に示されているように、TAPSによる変換は、すべてのgDNAにおいて、断片の長さにかかわらず、完全なものであった。
〔図10D〕TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。TAPSの前及び後の222bpのモデルDNAのアガロースゲル画像であり(独立して3回反復)、その反応後、検出可能な分解は見られなかった。
〔図11〕TAPSの前及び後のモデルDNA間での増幅曲線及び融解曲線の比較である。qPCRアッセイによって、TAPSの前及び後のモデルDNAで、増幅曲線において小さな違いが見られた。メチル化DNA(5mC)の融解曲線は、TAPSの後、低い温度の方にシフトしたことから、Tmが低下する、CからTへの変換の可能性が示されたが、非メチル化DNA(C)では、シフトは見られなかった。
〔図12〕サンガーシーケンシングによって示されたように、TAPS、TAPSβ及びCAPSの後、CからTへの完全な変換が誘導された。本明細書に記載されているようにして、1つのメチル化CpG部位及び1つのヒドロキシメチル化CpG部位を含むモデルDNAを調製した。本明細書に記載されているように、NgTET1による酸化及びピリジンボランによる還元のプロトコールに従って、TAPSによる変換を行った。5hmCのブロック、NgTET1による酸化及びピリジンボランによる還元のプロトコールに従って、TAPSβによる変換を行った。5hmCの酸化及びピリジンボランによる還元のプロトコールに従って、CAPSによる変換を行った。変換後、変換されたDNA試料1ngをTaq DNAポリメラーゼによってPCR増幅し、サンガーシーケンシング用に処理した。TAPSは、5mC及び5hmCの両方をTに変換した。TAPSβは、5mCを選択的に変換し、それに対し、CAPSは、5hmCを選択的に変換した。これらの3つの方法のいずれでも、無修飾シトシン及びその他の塩基は変換されなかった。
〔図13A〕TAPSは、様々なDNAポリメラーゼ及びRNAポリメラーゼと適合しており、サンガーシーケンシングによって示されているように、CからTへの完全な変換を誘導する。ポリメラーゼ試験用の、メチル化CpG部位を含むモデルDNAとプライマー配列は、本明細書に記載されている。TAPS処理後、5mCは、DHUに変換された。KAPA HiFi Uracil Plusポリメラーゼ、Taqポリメラーゼ及びVentエキソポリメラーゼは、DHUをTとして読み取るので、PCR後、CからTへの完全な変換を誘導する。あるいは、Klenow断片、Bst DNAポリメラーゼ及びphi29 DNAポリメラーゼを含め、ビオチン標識プライマー及び等温ポリメラーゼで、プライマー伸長を行った。新たに合成したDNA鎖をDynabeads MyOne Streptavidin C1によって分離してから、Taqポリメラーゼを用いたPCRによって増幅し、サンガーシーケンシング用に処理した。T7 RNAポリメラーゼは、DHUを効率的にバイパスして、DHU部位の反対側にアデニンを挿入できたが、これは、RT-PCR及びサンガーシーケンシングによって証明されている。
〔図13B〕特定の他の市販されているポリメラーゼは、DHUを含むDNAを効率的には増幅しなかった。TAPS処理後、5mCは、DHUに変換された。KAPA HiFi Uracil Plusポリメラーゼ及びTaqポリメラーゼは、DHUをTとして読み取るので、CからTへの完全な変換を誘導することになる。KAPA HiFiポリメラーゼ、Pfuポリメラーゼ、Phusionポリメラーゼ及びNEB Q5ポリメラーゼを含め、特定の他の市販されているポリメラーゼでは、CからTへの変換が低度であったか、または変換されなかったことが観察された(不図示)。
〔図14〕DHUでは、T及びCと比べて、PCRバイアスが見られない。本明細書に記載されているように、対応するDNAオリゴを用いて、1つのDHU/U/T/C修飾を含むモデルDNAを合成した。qPCR反応に基づき、DHU/U/T/C修飾を有する各モデルDNAの標準曲線をプロットし、モデルDNAインプットは、1:10の段階希釈を行った(0.1pg~1ng。すべてのqPCR実験をトリプリケートで行った)。ExcelのSLOPE関数によって、対数濃度(ng)の値とCt平均値との回帰の傾きを計算した。効率(%)=(10^(-1/傾き)-1)×100%という式を用いて、PCR効率を計算した。増幅係数=10^(-1/傾き)という式を用いて、増幅係数を計算した。DHU、TまたはCの修飾を有するモデルDNAでのPCR効率は、ほぼ同じであったことから、DHUは、通常の塩基として読み取ることができ、PCRバイアスの原因とならないと見られることが示された。
〔図15〕TAPSは、DNA断片の長さにかかわらず、5mCをTに完全に変換した。(A)TaqαI-消化アッセイのアガロースゲル画像により、すべての試料において、DNA断片の長さにかかわらず、5mCからTへの完全な変換が確認された。TAPS後、λゲノムに由来する194bpのモデル配列をPCR増幅し、TaqαI酵素で消化した。変換されなかった試料から増幅されたPCR産物は、切断できたが、TAPS処理試料で増幅された産物は、インタクトなままであったことから、制限部位の喪失、すなわち、5mCからTへの完全な変換が示唆されている。(B)ゲルバンドの定量によって、CからTへの変換パーセンテージを推定したところ、試験したすべてのDNA断片長で、100%となった。
〔図16〕異なるTAPS条件での変換及び偽陽性である。mTet1及びピリジンボランを組み合わせたところ、NgTET1またはpic-ボランを用いた他の条件と比べて、メチル化Cの変換率が最も高くなり(96.5%。完全にCpGがメチル化されたλDNAで計算した)と、無修飾のCの変換率が最も低くなった(0.23%、2kbの無修飾スパイクインで計算した)。バーの上に示されているのは、試験したすべてのシトシン部位の変換率±標準誤差である。
〔図17〕短いスパイクインにおける変換率である。5mC及び5hmCを含む(A)120mer-1及び(B)120mer-2である。両方の鎖の5mC及び5hmCの部位で、ほぼ完全な変換が行われた。実際の配列は、修飾の状態とともに、上下に示されている。
〔図18A〕全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。TAPSで処理したDNAにおける5mC及び5hmCの変換率である。左:既知の位置においてメチル化またはヒドロキシメチル化された合成スパイクイン(CpN)である。
〔図18B〕全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。無修飾の2kbのスパイクインから得られた、TAPSの偽陽性率である。
〔図18C〕全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。100万個のシミュレーションリードを1つのIntel Xeon CPUの1つのコアで処理した時の、TAPS及びWGBSの総実行時間である。
〔図18D〕全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。シーケーシングしたすべてのリードペア(トリミング後)のうち、ゲノムにマッピングされた割合である。
〔図18E〕全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。シーケンシングしたすべてのリードペアにおける第1及び第2のリードにおける1塩基あたりのシーケンシングクオリティスコアをIllumina BaseSpaceによって記録したものである。上は、TAPSである。下は、WGBSである。
〔図19〕TAPSでは、WGBSよりも均一なカバレッジが得られ、カバーされなかった位置がWGBSよりも少なかった。WGBSとTAPSの間で、(A)すべての塩基及び(B)CpG部位におけるカバレッジ深度を比較したものであり、両方の鎖で計算した。「TAPS(ダウンサンプリング済み)」では、マッピングされたすべてのTAPSリードから取ったランダムなリードを選択して、カバレッジ中央値が、WGBSのカバレッジ中央値と一致するようにした。最後のビンには、カバレッジが50×超の位置が示されている。
〔図20〕すべての染色体にわたる修飾レベルの分布である。マウス染色体全体にわたる、100kbのウィンドウにおける改変レベルの平均を、CpGのカバレッジによって重み付けして、Gaussian加重移動平均フィルターを用いて、ウィンドウサイズ10で平滑化したものである。
〔図21A〕TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。すべてのマウスCpGアイランド(20個のウィンドウにビニングした)及び4kbpのフランキング領域(サイズの等しい50個のウィンドウにビニングした)におけるシーケンシングカバレッジ深度の平均である。シーケンシング深度の差を考慮するために、マッピングされたすべてのTAPSリードをダウンサンプリングして、ゲノム全体にわたる、マッピングされたWGBSリードの数の中央値と適合させた。
〔図21B〕TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。TAPSのみ、TAPS及びWGBSの両方、またはWGBSのみによって、少なくとも3本のリードによってカバーされたCpG部位である。
〔図21C〕TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。少なくとも3本のリードによってカバーされたとともに、TAPSのみ、TAPS及びWGBS、またはWGBSのみによって検出された修飾レベルが0.1超であるCpG部位の数である。
〔図21D〕TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。TAPS及びWGBSにおける修飾レベル(%)の染色体分布の例である。マウス染色体4における、100kbのウィンドウ当たりの修飾CpGの割合の平均を、Gaussian加重移動平均フィルターを用いて、ウィンドウサイズ10で平滑化したものである。
〔図21E〕TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。TAPS及びWGBSの両方において、少なくとも3本のリードによってカバーされたCpG部位の数を表すヒートマップであり、各方法によって測定した場合の修飾レベルによって分類したものである。コントラストを向上させるために、両方の方法で無修飾のCpGを含む第1のビンは、カラースケールから除外し、星印で示されている。
〔図22〕CpGアイランド周辺の修飾レベルである。CpGアイランド(20個のウィンドウにビニングした)及び4kbpのフランキング領域(サイズの等しい50個のウィンドウにビニングした)における修飾レベルの平均である。カバレッジが3リード未満のビンは無視した。
〔図23〕TAPSでは、WGBSよりもカバレッジ-修飾バイアスが小さい。すべてのCpG部位を、そのカバレッジに従ってビニングし、各ビンにおいて、修飾値の平均(●)及び中央値(黒三角)が、WGBS(A)及びTAPS(B)に関して示されている。最後のビンには、100個超のリードによってカバーされたCpG部位が示されている。線は、データポイントにわたる直線近似を表している。
〔図24〕dsDNA及びssDNAライブラリー調製キットを用いて、インプットの少ないgDNA及びセルフリーDNA TAPSを調製した。(A)dsDNAライブラリーキットNEBNext Ultra IIまたはKAPA HyperPrepキットを用いて、わずか1ngのマウス胚性幹細胞(mESC)ゲノムDNA(gDNA)による、シーケンシングライブラリーの構築が成功した。ssDNAライブラリーキットAccel-NGS Methyl-Seqキットを用いて、インプットDNA量をさらに、(B)0.01ngのmESC gDNAまたは(C)1ngのセルフリーDNAまで減少させた。
〔図25〕dsDNA KAPA HyperPrepライブラリー調製キットを用いて、インプットの少ないgDNA及びセルフリーDNAのTAPSライブラリーを調製した。KAPA HyperPrepキットを用いて、1ngほどの少ない(A)mESC gDNA及び(B)セルフリーDNAで、シーケンシングライブラリーの構築が成功した。セルフリーDNAは、血漿ヌクレアーゼ消化により、160bp前後(ヌクレオソームサイズ)のシャープな長さ分布を有する。ライブラリーの構築後、そのDNAは、約300bpとなり、それは、Bにおけるシャープなバンドである。
〔図26〕高品質なセルフリーDNA TAPSである。(A)TAPSで処理したcfDNAにおける5mCの変換率である。(B)TAPSで処理したcfDNAにおける偽陽性率である。(C)シーケンシングされたすべてのリードペアのうち、ゲノムにユニークにマッピングされた割合である。(D)シーケンシングされたすべてのリードペアのうち、ゲノムにユニークにマッピングされ、PCR重複リードを除去した後の割合である。CHG及びCHHは、非CpGコンテクストである。
〔図27〕TAPSは、遺伝的バリアントを検出できる。メチル化(MOD、上の行)及びC→T SNP(下の行)では、元の上鎖(OT)/元の下鎖(OB)(左のカラム)、ならびにOTと相補的な鎖(CTOT)及びOBと相補的な鎖(CTOB)(右のカラム)において、別々の塩基分布が見られた。
〔図28〕TAPS変換生成物のエンドヌクレアーゼ切断である。(A)TAPS変換の前及び後のmESC gDNAの、異なるエンドヌクレアーゼによる消化の結果である。TAPS変換では、メチル化シトシンに代わって、DHUが導入される。Endo VIIIはエンドヌクレアーゼVIII、Endo IVはエンドヌクレアーゼIV、Tma Endo IIIはエンドヌクレアーゼIII、Tth Endo IVはエンドヌクレアーゼIV、Nth Endo IIIはNthエンドヌクレアーゼIII、Endo VはエンドヌクレアーゼVである。(B)TAPSで処理したmESC gDNAであって、USER消化前、USER消化後、及びサイズ選別後のmESC gDNAの代表的な画像である。(C)BS-seq及びeeTAPSの両方によって測定したすべてのCpGにおけるメチル化レベルを示す散布図である。
〔図29A〕エンドヌクレアーゼ濃縮TAPS(eeTAPS)である。eeTAPSの概略図(上図)と、CGメチル化レベルの、コンピューターによる測定値(下図)である。まず、TAPSで5-メチルシトシン(mC)をジヒドロウラシル(DHU)に変換してから、USER消化を通じて濃縮した。その後、サイズ選別を行ったDNA断片を調製して、シーケンシングライブラリーにし、PCRによって増幅した。リードのアラインメント後、各CpG部位で切断されるリードの数を、各CpG部位で切断されたリードの総数、または各CpG部位をカバーするリードの総数で除したものとして、CGメチル化レベルを計算した。
〔図29B〕4kbのモデルDNAでのeeTAPSの検証である。トラックは、上から下の順で、バイサルファイトシーケンシング(BS-seq)、eeTAPS、及びeeTAPSコントロールで測定したメチル化レベルを示している。eeTAPSコントロールでは、TAPSで変換せずに、USER酵素を用いて、DNAを消化した。
〔図30A〕mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS、eeTAPS及びrrTAPSのゲノムの断片化方法を示す図である。wgTAPSでは、ゲノムDNAは、ランダムに断片化されるのに対して、断片化は、eeTAPSでは、メチル化CpG(mCG)部位で、rrTAPSでは、CCGG部位で特異的に行われる。
〔図30B〕mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS、eeTAPS及びrrTAPSによってカバーされたCpG部位のパーセンテージを、クロマチンの異なる特徴と重ね合わせたものを示すバープロットである。これらのクロマチンの特徴は、以前の研究(20)で定義されている。
〔図30C〕mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS及びeeTAPSによって、1CpGの分解能で求めたメチル化レベルの相関性を示すヒートマップである。wgTAPS及びeeTAPSの両方において、メチル化レベルを16群に分けた。カラーは、所定の間隔でのCpGの数を示している。wgTAPSでのカバレッジが5超の部位のみを考慮に入れた。プロットの上部には、ピアソンの相関係数が示されている。
〔図30D〕mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS、rrTAPS及びeeTAPSにおける、検出されたmCpG部位の重複を示すベンプロットである。wgTAPS及びeeTAPSの両方において、メチル化レベルが、メチル化レベルの第1四分位超のCpG部位を、メチル化CpGとして定義した。rrTAPSでは、wgTAPSと同じメチル化カットオフ値を用いた。(wgTAPS及びrrTAPSでは、メチル化レベルが0.5超のCpG部位を、mCpGとして定義し、eeTAPSでは、メチル化レベルが0.28超のCpG部位をmCpGとして定義した。)
〔図31〕(A)シーケンシングした断片の末端におけるヌクレオチド頻度である。(B)切断部位と、その部位に最も近いCpGとの距離の分布を示すバープロットである。(C)wgTAPS及びrrTAPSによって、1CpGの分解能で求めたメチル化レベルの相関性を示すヒートマップである。wgTAPS及びrrTAPSの両方において、メチル化レベルを16群に分けた。カラーは、所定の間隔でのCpGの数を示している。wgTAPSでのカバレッジが5超であり、rrTAPSでのカバレッジが5超である部位のみを考慮に入れた。ピアソンの相関係数は、0.92である。(D)eeTAPSのレプリケートで検出されたmCpG部位の重複であり、レプリケートは、同じ深度になるようにサブサンプリングして、mCpGを検出した。
〔図32〕異なるゲノム機構におけるメチル化をeeTAPSでプロファイリングしたものである。(A)mESCの1番染色体において、eeTAPS(上の線)及びwgTAPS(下の線)の両方によって測定したメチル化レベルである。100kbのウィンドウを使用し、移動平均値をRにおいて、bw=10で、movAvg2関数によって計算した。(B)eeTAPS及びwgTAPS間での、染色体ビンにおけるメチル化レベルの相関性を示す密度プロットである。100kbのウィンドウを用いて、各ビンにおいて、メチル化レベルの平均を計算した。プロットの上部には、ピアソンの相関係数が示されている。(C)eeTAPS及びwgTAPSにおける、CpGアイランド(CGI)及び4kbのフランキング領域にわたるメチル化レベルの平均である。(D)eeTAPS及びwgTAPS間での、CpGアイランドにおけるメチル化レベルの相関性を示す密度プロットである。プロットの上部には、ピアソンの相関係数が示されている。(E)eeTAPS及びwgTAPSでの、クロマチンの異なる特徴におけるメチル化レベルの相関性を示す密度プロットである。これらのクロマチンの特徴は、ヒストンマーカーによって以前に定義されたものであり(Bogu,G.K.,Vizan,P.,Stanton,L.W.,Beato,M.,Di Croce,L.and Marti-Renom,M.A.(2015)Chromatin and RNA Maps Reveal Regulatory Long Noncoding RNAs in Mouse. Mol Cell Biol,36,809-819)、(F)にも示されている。プロットの上部には、ピアソンの相関係数を示した。(F)クロマチンのすべての特徴にわたるメチル化レベルの分布をeeTAPS(下のバー)及びwgTAPS(上のバー)によって測定したものを示すボックスプロットである。
〔図33〕wgTAPS及びeeTAPSにおける、転写開始部位(TSS)周辺のメチル化分布の平均である。遺伝子は、GSE72855のデータセットに従い、それらの発現レベルによって分類した。
〔図34〕インプットの少ない試料での、eeTAPSによる解析である。(A)eeTAPSを用いて、1ng、10ng、50ng、200ngのmESC gDNAインプットで検出されたmCpGの数(wgTAPSによって特定されたmCpGのうちの数)である。200ngのmESCでは、リードを2分の1にダウンサンプリングして、インプットの少ない試料に合わせて、シーケンス深度を適合させた。0.28超のメチル化レベルという基準を用いて、mCpGを割り当て、eeTAPSでは、1超の切断数をmCpGとして割り当て、wgTAPSでは、0.5超のメチル化レベルをmCpGとして割り当てた。バーの上に示されているパーセンテージは、検出されたmCpGのパーセンテージである(wgTAPSで検出されたmCpGは、真として定義する)。(B)異なるインプットレベルを用いて、マウスゲノムにわたり、eeTAPSで測定したメチル化分布を示すヒートマップである。各染色体を100kbのウィンドウに分け、それらのウィンドウは、ヒートマップの行によって表されている。メチル化レベルは、100kbの各ウィンドウにおいて、メチル化CpG部位の数を、カバーされたCpG部位の総数で除したものとして定義した。(C)インプットの少ない試料間での、200ngのインプット試料に対するメチル化の相関性を示す密度プロットである。メチル化レベルは、(B)に示されているように、100kbのウィンドウにおいて、全マウスゲノムにわたり計算した。ピアソンの相関係数は、各プロットについて示されている。
〔図35〕eeTAPSシーケンシング深度の解析結果である。(A)1~10×のシーケンシング深度でリードをサンプリングしたときに検出されるメチル化CpGの数である。上に示されているパーセンテージは、eeTAPSによって検出されたmCpGのパーセンテージである(wgTAPSで検出されたmCpGを真と定義する)。(B)1~10×のシーケンシング深度でリードをサンプリングしたときの、100kbのウィンドウにおける、全マウスゲノム(上の線)及びCpGアイランド(CGI)(下の線)におけるメチル化の相関性である。
〔発明を実施するための形態〕
本開示は、バイサルファイトを使用せずに、塩基単位の分解能で、全ゲノムDNAを含むDNA試料中の標的DNAにおけるシトシン修飾を特定する方法を提供する。本明細書に記載されている方法は、PCT/US2019/012627(参照により、その全体が本明細書に援用される)に記載されている方法に対する改良点を含み、そのPCT特許には、バイサルファイトを使用せずに、塩基単位の分解能で、配列における5mC及び5hmCを検出する方法(TAPSという)を含む方法が記載されている。TAPSは、無修飾シトシンに影響を及ぼすことなく、5mC及び5hmCを直接かつ定量的に、塩基単位の分解能で検出するための穏和な酵素反応及び化学反応からなる。本開示は、無修飾シトシンに影響を及ぼすことなく、5fC及び5caCを塩基単位の分解能で検出するための改良型の方法も提供する。すなわち、本発明で提供する方法は、5mC、5hmC、5fC及び5caCのマッピングを行うものであり、バイサルファイトシーケンシングのような以前の方法の問題点を解消する。
本明細書に記載されている方法は、例えばTetアシストピリジンボラン処理によって、修飾シトシン(5mC、5hmC、5fC、5caC)をジヒドロウラシル(DHU)に変換してから、DHU部位を、(例えばUSER(ウラシル特異的除去試薬)によって)切断して、DNA断片を作製した後、その断片をシーケンシングライブラリーにすることをベースとする。断片化されなかったゲノムDNAは、直接シーケンシングすることはできない。DHUに変換される修飾シトシンが存在する場合のみ、DNAは切断されて、DNA断片となった後、その断片の各末端(修飾シトシンの位置を示す)により、シーケンシングできる。すなわち、例えば、eeTAPSは、5mC/5hmCのDHUへの変換により、切断された断片のみをシーケンシングすることによって、全ゲノムTAPS(WGTAPS)のコストを削減する。切断された断片をシーケンシングすることによって、元のDNA試料におけるメチル化シトシンを塩基単位の分解能で特定できる。
5mCまたは5hmCを(一緒に)特定する方法
別の態様では、本開示は、DNA試料における5mCまたは5hmCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNAを改変する工程(以下の工程を含む)と、
i.そのDNA試料中の5mC及び5hmCを5-カルボキシルシトシン(5caC)及び/または5fCに変換する工程
ii.その5caC及び/または5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5mCまたは5hmCのいずれかの位置が得られる方法を提供する。
標的DNAにおける5mCまたは5hmCの特定方法の実施形態では、その方法は、標的DNAにおいて、5mCまたは5hmCの改変が認められた各位置における、5mCまたは5hmCの改変頻度に関する半定量的な尺度を提供する。5mCまたは5hmCを特定するこの方法は、5mC及び5hmCの位置をもたらすが、これらの2つのシトシン修飾の区別は行わない。それどころか、5mC及び5hmCの両方をDHUに変換する。
5mCの特定方法
一態様では、本開示は、DNA試料における5-メチルシトシン(5mC)の特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNAを改変する工程(以下の工程を含む)と、
i.そのDNA試料中の5-ヒドロキシメチルシトシン(5hmC)に、ブロック基を付加する工程
ii.そのDNA試料中の5mCを5-カルボキシルシトシン(5caC)及び/または5-ホルミルシトシン(5fC)に変換する工程
iii.その5caC及び/または5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程であって、切断部位の存在により、その標的DNAにおける5mCの位置をもたらす工程と、
を含む方法を提供する。
標的DNAにおける5mCの特定方法の実施形態では、その方法は、標的DNAにおいて、5mCの改変が認められた各位置における、5mCの改変頻度に関する半定量的な尺度をもたらす。
5hmCを含めずに、標的DNAにおける5mCを特定するために、試料中の5hmCをブロックして、5hmCが、5caC及び/または5fCに変換されないようにする。本開示の方法では、試料DNAにおける5hmCは、ブロック基を5hmCに付加することによって、次の工程で反応しないようになる。一実施形態では、そのブロック基は、修飾糖を含む糖、例えば、グルコースまたは6-アジド-グルコース(6-アジド-6-デオキシ-D-グルコース)である。その糖ブロック基は、1つ以上のグルコシルトランスフェラーゼ酵素の存在下で、DNA試料をウリジン二リン酸(UDP)-糖と接触させることによって、5hmCのヒドロキシメチル基に付加する。
実施形態では、そのグルコシルトランスフェラーゼは、T4バクテリオファージβ-グルコシルトランスフェラーゼ(βGT)、T4バクテリオファージα-グルコシルトランスフェラーゼ(αGT)、ならびにこれらの誘導体及び類似体である。βGTは、DNAにおいて、β-D-グルコシル(グルコース)残基をUDP-グルコースから5-ヒドロキシメチルシトシン残基に移転する化学反応を触媒する酵素である。
ブロック基が、例えばグルコースであるとは、グルコース部分(例えば、β-D-グルコシル残基)が5hmCに付加されて、グルコシル5-ヒドロキシメチルシトシンが得られることを指す。その糖ブロック基は、グルコシルトランスフェラーゼ酵素の基質であるとともに、後で、5hmCが5caC及び/または5fCに変換されないようにブロックするいずれの糖または修飾糖であることもできる。続いて、TET酵素を用いた酸化によるような、本発明で提供する方法によって、DNA試料中の5mCを5caC及び/または5fCに変換する工程を行う。そして、ボランによる酸化のような、本発明で提供する方法によって、その5caC及び/または5fCのDHUへの変換を行う。
5mCを特定するとともに、5hmCを特定する方法
本開示は、(i)本明細書に記載されている第1のDNA試料で、5mCを特定する方法を行い、(ii)本明細書に記載されている第2のDNA試料で、5mCまたは5hmCを特定する方法を行うことによって、標的DNAにおいて、5mCを特定するとともに、5hmCを特定する方法を提供する。(i)によって、5mCの位置が得られる。(i)及び(ii)の結果を比較することによって、(i)では存在するが、(ii)では存在しない切断部位によって、標的DNAにおける5hmCの位置が得られる。実施形態では、第1及び第2のDNA試料は、同じDNA試料に由来する。例えば、第1及び第2の試料は、解析対象のDNAを含む試料から取った別々のアリコートであってよい。
5mC及び5hmC(ブロックされていない)を5fC及び5caCに変換してから、DHUに変換するので、DNA試料に存在するすべての5fC及び5caCが、5mC及び/または5hmCとして検出されることになる。しかしながら、通常の条件下では、ゲノムDNAにおける5fC及び5caCのレベルが極めて低いことを考慮すると、DNA試料中のメチル化及びヒドロキシメチル化を解析するときには、上記の事象は、許容されることが多いであろう。その5fC及び5caCのシグナルは、例えば、5fCにおいては、ヒドロキシルアミンの結合によって、5caCにおいてEDCのカップリングによって、5fC及び5caCを保護して、DHUに変換されないようにすることによって除去できる。
上記の方法は、5mCの位置と、5mCまたは5hmC(組み合わさったもの)の位置とを比較することを通じて、標的DNAにおける5hmCの位置を特定する。あるいは、標的DNAにおける5hmCの改変位置を直接測定できる。すなわち、一態様では、本開示は、DNA試料中の5hmCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNAを改変する工程(以下の工程を含む)と、
i.そのDNA試料中の5hmCを5caC及び/または5fCに変換する工程
ii.その5caC及び/または5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5hmCの位置が得られる方法を提供する。
実施形態では、その5hmCを5fCに変換する工程は、そのDNAと、例えば、過ルテニウム酸カリウム(KRuO)(Science.2012,33,934-937及びWO2013017853(参照により、本明細書に援用される)に記載されているようなもの)またはCu(II)/TEMPO(過塩素酸銅(II)及び2,2,6,6-テトラメチルピペリジン-1-オキシル(TEMPO))(Chem.Commun.,2017,53,5756-5759及びWO2017039002(参照により、本明細書に援用される)に記載されているようなもの)と、を接触させることによって、5hmCを酸化して、5fCにすることを含む。使用できる他の酸化剤は、ルテニウム酸カリウム及び/または酸化マンガンである。続いて、本明細書に開示されている方法によって、例えばボラン反応によって、DNA試料中の5fCをDHUに変換する。
5caCまたは5fCの特定方法
一態様では、本開示は、DNA試料中の5caCまたは5fCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.5caC及び/または5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
f.その改変標的DNAを切断する工程と、
g.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
h.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCまたは5fCのいずれかの位置が得られる方法を提供する。
5fCまたは5caCを特定するこの方法は、5fCまたは5caCの位置をもたらすが、これらの2つのシトシン修飾の区別は行わない。むしろ、5fC及び5caCの両方をDHUに変換し、そのDHUを、本明細書に記載されている方法によって検出する。
5caCの特定方法
別の態様では、本開示は、DNA試料中の5caCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5fCに、ブロック基を付加する工程と、
c.その5caCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
g.その改変標的DNAを切断する工程と、
h.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
i.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCの位置が得られる方法を提供する。
標的DNAにおける5caCの特定方法の実施形態では、その方法は、標的DNAにおいて、5caCの改変が認められた各位置における、5caCの改変頻度に関する半定量的な尺度をもたらす。
この方法では、5fCをブロックし(5mC及び5hmCは、DHUに変換されない)、標的DNAにおける5caCを特定可能にする。実施形態では、DNA試料中の5fCに、ブロック基を付加する工程は、そのDNAと、アルデヒド反応性化合物(例えば、ヒドロキシルアミン誘導体、ヒドラジン誘導体及びヒドラジド誘導体を含む)とを接触させることを含む。ヒドロキシルアミン誘導体としては、ヒドロキシルアミン、ヒドロキシルアミン塩酸塩、ヒドロキシルアンモニウム酸の硫酸塩、リン酸ヒドロキシルアミン、O-メチルヒドロキシルアミン、O-ヘキシルヒドロキシルアミン、O-ペンチルヒドロキシルアミン、O-ベンジルヒドロキシルアミン、及び特に、O-エチルヒドロキシルアミン(EtONH2)、O-アルキル化ヒドロキシルアミンまたはO-アリール化ヒドロキシルアミン、これらの酸または塩が挙げられる。ヒドラジン誘導体としては、N-アルキルヒドラジン、N-アリールヒドラジン、N-ベンジルヒドラジン、Ν,Ν-ジアルキルヒドラジン、N,N-ジアリールヒドラジン、Ν,Ν-ジベンジルヒドラジン、N,N-アルキルベンジルヒドラジン、Ν,Ν-アリールベンジルヒドラジン及びN,N-アルキルアリールヒドラジンが挙げられる。ヒドラジド誘導体としては、-トルエンスルホニルヒドラジド、N-アシルヒドラジド、Ν,Ν-アルキルアシルヒドラジド、N,N-ベンジルアシルヒドラジド、Ν,Ν-アリールアシルヒドラジド、N-スルホニルヒドラジド、N,N-アルキルスルホニルヒドラジド、Ν,Ν-ベンジルスルホニルヒドラジド及びN,N-アリールスルホニルヒドラジドが挙げられる。
5fCの特定方法
別の態様では、本開示は、DNA試料中の5fCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5caCに、ブロック基を付加する工程と、
c.その5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
d.その改変標的DNAを切断する工程と、
e.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
f.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5fCの位置が得られる方法を提供する。
標的DNAにおける5fCの特定方法の実施形態では、その方法は、標的DNAにおいて、5fCの改変が認められた各位置における、5fCの改変頻度に関する半定量的な尺度をもたらす。
そのDNA試料中の5caCに、ブロック基を付加する工程は、(i)そのDNA試料と、カップリング剤、例えば、1-エチル-3-(3-ジメチルアミノプロピル)カルボジイミド(EDC)またはΝ,Ν’-ジシクロヘキシルカルボジイミド(DCC)などのカルボジイミド誘導体のようなカルボン酸誘導体化試薬とを接触させ、(ii)そのDNA試料と、アミン、ヒドラジンまたはヒドロキシルアミン化合物とを接触させることによって行うことができる。すなわち、例えば、DNA試料をEDCで処理してから、ベンジルアミン、エチルアミンまたはその他のアミンで処理して、例えば5caCをブロックして、pic-BHによってDHUに変換されないようにするアミドを形成させることによって、5caCをブロックできる。EDC触媒によって、5caCをカップリングする方法は、WO2014165770に記載されており、参照により、本明細書に援用される。
DNA試料/標的DNA
本開示の方法では、バイサルファイトシーケンシングのような方法に伴う大幅な分解を回避する穏和な酵素反応及び化学反応を利用する。すなわち、本発明の方法は、循環セルフリーDNAのような、インプットの少ない試料の解析、及びシングルセル解析に有用である。
実施形態では、DNA試料は、ピコグラム量のDNAを含む。実施形態では、DNA試料は、DNAを約1pg~約900pg、DNAを約1pg~約500pg、DNAを約1pg~約100pg、DNAを約1pg~約50pg、DNAを約1~約10pg、約200pg未満、DNAを約100pg未満、DNAを約50pg未満、DNAを約20pg未満、及びDNAを約5pg未満含む。別の実施形態では、DNA試料は、ナノグラム量のDNAを含む。本明細書に開示されている方法で用いる試料のDNAは、単一の細胞に由来するDNAまたはバルクDNA試料を含め、いずれの量であることもできる。実施形態では、本発明の方法は、DNAを約1~約500ng、DNAを約1~約200ng、DNAを約1~約100ng、DNAを約1~約50ng、DNAを約1~約10ng、DNAを約2~約5ng、DNAを約100ng未満、DNAを約50ng未満、DNAを5ng未満、及びDNAを2ng未満含むDNA試料で行うことができる。実施形態では、そのDNA試料は、マイクログラム量のDNAを含む。
DNA試料を準備するとは、本明細書で使用する場合、DNA試料をいずれかの供給源から、直接または間接的のいずれかで得ることを指す。本明細書に記載されている方法で用いるDNA試料は、例えば、体液、組織試料、器官、細胞小器官または単一の細胞を含むいずれかの供給源に由来するものであってよい。実施形態では、DNA試料は、循環セルフリーDNA(セルフリーDNAまたはcfDNA)であり、このDNAは、血液中に見られるDNAであり、細胞内には存在しない。cfDNAは、当該技術分野において知られている方法を用いて、血液または血漿から単離できる。cfDNAを単離するには、例えばCirculating DNA Kit(Qiagen)を含む市販のキットが利用可能である。DNA試料は、濃縮工程に起因してもよく、その工程としては、抗体免疫沈降、クロマチン免疫沈降、制限酵素消化ベースの濃縮、ハイブリダイゼーションベースの濃縮または化学標識ベースの濃縮が挙げられるが、これらに限らない。
標的DNAは、シトシン修飾(すなわち、5mC、5hmC、5fC及び/または5caC)を有するいずれのDNAであってもよく、そのDNAとしては、組織、器官、細胞及び細胞小器官から精製したDNA断片またはゲノムDNAが挙げられるが、これらに限らない。標的DNAは、試料中の単一のDNA分子であることができ、またはシトシン修飾を有する試料(またはそのサブセット)中のDNA分子集団全体であってもよい。標的DNAは、供給源に由来する天然型のDNAであることができ、または例えば、断片化、修復及びシーケンシング用のアダプターとのライゲーションによって、ハイスループットシーケンシングを行える状態に事前に変換されていることもできる。すなわち、本明細書に記載されている方法を用いて、(例えば、個別の標的の配列を求めることによって)個別に解析したり、または(例えば、ハイスループットシーケンシング法または次世代シーケンシング法によって)集団で解析したりできる標的DNA配列のライブラリーを作製できるように、標的DNAは、複数のDNA配列を含むことができる。
標的DNAを含むDNA試料は、モネラ界(細菌)、原生生物界、菌類界、植物界及び動物界の生物から得ることができる。DNA試料は、患者もしくは対象、環境試料、または該当する生物から得てよい。実施形態では、DNA試料は、細胞、細胞群、体液、組織試料、器官及び/または細胞小器官から抽出するか、精製するか、または取り出す。好ましい実施形態では、試料DNAは、全ゲノムDNAである。
5mC及び5hmCの5caC及び/または5fCへの変換
本明細書に記載されているeeTAPS法など、本発明で提供する方法の実施形態は、5mC及び5hmC(または5hmCをブロックする場合には、5mCのみ)を5caC及び/または5fCに変換する工程を含む。実施形態では、この工程は、DNA試料を10-11転座(TET)酵素と接触させることを含む。TET酵素は、5mCのN5メチル基への酸素分子の移転を触媒する酵素のファミリーであり、この移転により、5-ヒドロキシメチルシトシン(5hmC)が形成される。TETはさらに、5hmCを酸化して5fCにし、5fCを酸化して5caCを形成させるのを触媒する(図5Aを参照されたい)。本明細書に記載されている方法で有用なTET酵素としては、ヒトTET1、ヒトTET2、ヒトTET3、マウスTet1、マウスTet2、マウスTet3、Naegleria TET(NgTET)、Coprinopsis cinerea(CcTET)、及びこれらの誘導体または類似体のうちの1つ以上が挙げられる。実施形態では、そのTET酵素は、NgTETまたはその誘導体である。別の実施形態では、そのTET酵素は、ヒトTET1(hTET1)またはその誘導体である。実施形態では、そのTET酵素は、マウスTet1またはその誘導体(mTet1CD)である。別の実施形態では、そのTET酵素は、ヒトTET2(hTET2)またはその誘導体である。
5caC及び/または5fCのDHUへの変換
本明細書に記載されている方法は、DNA試料中の5caC及び/または5fCをDHUに変換する工程を含む。実施形態では、この工程は、DNA試料と、例えば、ボラン還元剤(ピリジンボラン、2-ピコリンボラン(pic-BH)、ボラン、水素化ホウ素ナトリウム、シアノ水素化ホウ素ナトリウム及びトリアセトキシ水素化ホウ素ナトリウムなど)を含む還元剤とを接触させることを含む。好ましい実施形態では、その還元剤は、ピリジンボラン及び/またはpic-BHである。
改変標的DNAを切断する工程
本明細書に記載されている方法は、変換工程(複数可)の前(すなわち、修飾シトシンをDHUに変換する1つまたは複数の工程の前)のDNA試料において修飾シトシン(5mC、5hmC、5fC及び/または5caC)があった位置に、DHUを含む改変標的DNAを切断する工程を含む。本明細書に記載されている切断工程は、DHUを含む改変標的DNAを特異的に切断する一方で、DHUを含まないDNAは、未切断または実質的に未切断のままにする。
DHUを含む改変標的DNAを切断する工程は、DHUを含む改変標的DNAと、その改変標的DNAを特異的に切断する1つ以上のDNAエンドヌクレアーゼとを接触させることによって行うことができる。実施形態では、そのDNAエンドヌクレアーゼのうちの1つ以上は、DNA N-グリコシラーゼ活性及びAPリアーゼ活性を有する2官能性DNAエンドヌクレアーゼであり、そのDNAエンドヌクレアーゼとしては、例えば、TmaエンドヌクレアーゼIII、エンドヌクレアーゼVIII、ホルムアミドピリミジンDNAグリコシラーゼ(Fpg)及び/またはhNEIL1が挙げられる。実施形態では、DHUを含む改変標的DNAは、ウラシル特異的除去試薬(USER)によって切断する。USER酵素は、ウラシルDNAグリコシラーゼ(UDG)及びDNAグリコシラーゼ-リアーゼエンドヌクレアーゼVIIIを組み合わせたものを含む。改変標的DNAを切断するのに使用できる他の酵素は、アプリン酸/アピリミジン酸エンドヌクレアーゼ1(APE1)、エンドヌクレアーゼIII(Endo III)、TmaエンドヌクレアーゼIII、TthエンドヌクレアーゼIV、エンドヌクレアーゼV、エンドヌクレアーゼVIII、Fpg及びhNEIL1のうちの1つ以上である。
実施形態では、DHUを含む改変標的DNAを切断する工程は、House CH,Miller SL Hydrolysis of dihydrouridine and related compounds.Biochemistry.1996;35(1):315-320に記載されているように、その改変標的DNAを酸性pH及び/または加熱条件に暴露することを含む。実施形態では、その切断工程は、その改変標的DNAを少なくとも70℃、75℃、80℃、85℃、90℃、95℃、100℃、105℃もしくは110℃の温度、及び/または7以上、7.5以上、8以上、8.5以上、9以上、9.5以上、10以上、10.5以上、11以上、11.5以上、12以上、12.5以上もしくは13以上のpHに暴露することを含む。
切断された改変標的DNAへのアダプターDNA分子の付加
実施形態では、本明細書に記載されている方法は、例えばUSER酵素によって切断した改変標的DNAに、アダプターDNA分子を付加する工程を含む。アダプターDNAまたはDNAリンカーは、化学合成した短い一本鎖または二本鎖のオリゴヌクレオチドであって、他のDNA分子の一方の末端または両方の末端にライゲーションできるオリゴヌクレオチドである。二本鎖のアダプターは、そのアダプターの各末端が、平滑末端、または5’もしくは3’オーバーハング(すなわち粘着末端)を有するように合成できる。DNAアダプターを、切断された改変標的DNAにライゲーションして、相補的なプライマーによるPCR増幅用、及び/またはクローニング及び/またはライブラリー作製(例えば次世代シーケンシングライブラリー)用の配列をもたらす。
アダプターを、切断された標的DNAにライゲーションする前に、例えば、5’リン酸基を有する平滑末端を作製する末端修復によって、その切断されたDNAの末端をライゲーション用に調製してよい。平滑末端をアダプターへのライゲーションに使用でき、またはライゲーションの前に、例えばテーリング反応によって、オーバーハングを作製できる。テーリングは、ノンテンプレートヌクレオチドを平滑な二本鎖DNA分子の3’末端に付加する酵素的な方法である。3’末端のAテーリング(すなわち、3’末端へのdAの付加)を用いて、アダプターと、相補的なdTオーバーハングとのライゲーションを促進できる。
実施形態では、切断された改変標的DNAに、DNAアダプター分子を付加する工程の前または後のいずれかに、その切断された標的DNAのサイズ選択を行う。実施形態では、そのサイズ選別は、切断された標的DNAにDNAアダプターを付加した後に行う。サイズ選別は、当該技術分野において知られている方法によって行うことができ、その方法としては、固相逆固定化(SPRI)常磁性ビーズ(例えば、AMPure XPビーズを使用するもの)が挙げられるが、これらに限らない。
改変標的DNAのコピー数の増幅
本明細書に記載されている方法は、当該技術分野において知られている方法によって、改変標的DNAのコピー数を増幅する(増加させる)工程を任意に含んでよい。改変標的DNAがDNAであるときには、そのコピー数は、例えば、PCR、クローニング及びプライマー伸長法によって増加できる。特定の標的DNA配列に特異的なプライマーを用いるPCRによって、個々の標的DNAのコピー数を増幅できる。あるいは、標準的な技法によって、DNAベクターにクローニングすることによって、複数の異なる改変標的DNA配列を増幅できる。実施形態では、複数の異なる改変標的DNA配列のコピー数をPCRによって増加させて、次世代シーケンシング用のライブラリーを作製し、その場合、例えば、二本鎖のアダプターDNAを事前に、試料DNA(または改変試料DNA)にライゲーションし、そのアダプターDNAと相補的なプライマーを用いて、PCRを行う。
次世代シーケンシングライブラリーの作製
切断された改変標的DNAにアダプターDNA分子を付加したら、その改変標的DNAのコピー数を(例えばPCRによって)増幅して、次世代シーケンシング用のライブラリーDNA配列を作製できる。PCR用のプライマーは、切断された標的DNAに事前にライゲーションしたアダプターDNAに対応する(アダプターDNAと相補的な)配列を有する。本発明で提供する方法(試薬、工程及びその順序を含む)により、ハイスループットな次世代シーケンシング法を用いてシーケンシングできるDNA配列のライブラリーを作製可能になる。
改変標的DNAの切断部位の検出
本明細書に開示されている方法の実施形態では、その方法は、切断された改変標的DNAの配列を検出する工程を含む。その改変標的DNAは、未改変の標的DNAにおいて、5mC、5hmC、5fC及び5caCの1つ以上が存在していた位置に、DHUを含む。DHUを含む改変標的DNAは、DHU感受性のエンドヌクレアーゼによる消化を含む、本明細書に記載されている方法によって切断する。続いて、切断された断片をシーケンシングライブラリーに変換でき、その各断片の開始点及び終点は、修飾シトシン(5mC、5hmC、5fCまたは5caC)の部位に対応する。これにより、メチル化CpG部位をゲノムワイドに濃縮可能となり、その一方で、メチル化されていないゲノムの大部分が除去される。したがって、切断部位を特定する方法であって、当該技術分野において知られているいずれかの方法によって、シトシン修飾を検出できる。このような方法としては、サンガーシーケンシング法、マイクロアレイ法及び次世代シーケンシング法のようなシーケンシング法が挙げられる。
キット
加えて、本開示は、標的DNAにおける5mC及び5hmCを特定するためのキットを提供する。このようなキットは、本明細書に記載されている方法によって5mC及び5hmCを特定するための試薬を含む。そのキットは、本明細書に記載されている方法によって、5caCを特定するための試薬、及び5fCを特定するための試薬も含んでよい。実施形態では、そのキットは、TET酵素、ボラン還元剤、及び本発明の方法を行うための説明を含む。さらなる実施形態では、そのTET酵素は、TET1もしくはTET2(またはこれらの誘導体)であり、そのボラン還元剤は、ピリジンボラン、2-ピコリンボラン(pic-BH)、ボラン、水素化ホウ素ナトリウム、シアノ水素化ホウ素ナトリウム及びトリアセトキシ水素化ホウ素ナトリウムからなる群のうちの1つ以上から選択されている。さらなる実施形態では、そのTET1酵素は、NgTet1、ヒトTET1またはマウスTet1であり、そのボラン還元剤は、ピリジンボラン及び/またはpic-BHである。別の実施形態では、そのTET酵素は、mTET2またはその誘導体である。
実施形態では、そのキットはさらに、5hmCブロック基及びグルコシルトランスフェラーゼ酵素を含む。さらなる実施形態では、その5hmCブロック基は、ウリジン二リン酸(UDP)-糖であり、その糖は、グルコースまたはグルコース誘導体であり、そのグルコシルトランスフェラーゼ酵素は、T4バクテリオファージβ-グルコシルトランスフェラーゼ(βGT)、T4バクテリオファージα-グルコシルトランスフェラーゼ(αGT)、ならびにこれらの誘導体及び類似体である。
実施形態では、そのキットはさらに、過ルテニウム酸カリウム(KRuO4)、Cu(II)/TEMPO(過塩素酸銅(II)及び2,2,6,6-テトラメチルピペリジン-1-オキシル(TEMPO))、ルテニウム酸カリウム、ならびに酸化マンガンのうちの1つ以上から選択した酸化剤を含む。
実施形態では、そのキットは、DNA試料における5fCをブロックするための試薬を含む。実施形態では、そのキットは、例えば、本明細書に記載されているようなヒドロキシルアミン誘導体、ヒドラジン誘導体及びヒドラジド誘導体を含むアルデヒド反応性化合物を含む。実施形態では、そのキットは、本明細書に記載されているような、5caCをブロックするための試薬を含む。
実施形態では、そのキットは、DNAを単離するための試薬を含む。実施形態では、そのキットは、インプットの少ないDNAを試料から単離するため、例えば、cfDNAを血液、血漿または血清から単離するための試薬を含む。実施形態では、そのキットは、ゲノムDNAを単離するための試薬を含む。
実施形態では、そのキットは、本明細書に記載されているように、DHUを含む改変標的DNAを切断するための酵素を1つ以上含む。実施形態では、そのキットは、本明細書に記載されているようなアダプターDNA分子を含む。加えて、そのキットは、切断された改変標的DNAに、そのアダプターDNA分子をライゲーションするための酵素を含んでもよい。
〔実施例〕
実施例1:TAPS及びWGTAPS
方法
モデルDNAの調製
MALDI及びHPLC-MS/MS試験用のDNAオリゴ:C、5mC及び5hmCを有するDNAオリゴヌクレオチド(「オリゴ」)をIntegrated DNA Technologies(IDT)から購入した。すべての配列及び修飾を図6-1及び図6-2並びに図7-1及び図7-2に見ることができた。5fCを有するDNAオリゴをCテーリング法によって合成した。すなわち、DNAオリゴ5’-GTCGACCGGATC-3’及び5’-TTGGATCCGGTCGACTT-3’をアニーリングしてから、5-ホルミル-2’-dCTP(Trilink Biotech)及びKlenow断片(3’→5’エキソヌクレアーゼ活性欠損変異体)(New England Biolabs)と、NEBuffer2において2時間、37℃でインキュベートした。その生成物をBio-Spin P-6 Gel Columns(Bio-Rad)で精製した。
Expedite 8900 DNA Synthesis Systemを用いて、標準的なホスホロアミダイト(Sigma)及び5-カルボキシ-dC-CEホスホロアミダイト(Glen Research)によって、5caCを有するDNAオリゴを合成した。その後の脱保護及び精製は、Glen-Pak Cartridges(Glen Research)をメーカーの指示に従って用いて行った。Voyager-DE MALDI-TOF(マトリックス支援レーザー脱離イオン化飛行時間型)Biospectrometry Workstationによって、精製オリゴヌクレオチドの特徴付けを行った。
変換試験用の222bpのモデルDNA:CpG部位を5個含む222bpのモデルDNAを作製するために、Taq DNAポリメラーゼ(New England Biolabs)を用いて、バクテリオファージλDNA(Thermo Fisher)をPCR増幅し、AMPure XPビーズ(Beckman Coulter)によって精製した。プライマー配列は、フォワードプライマーが5’-CCTGATGAAACAAGCATGTC-3’、リバースプライマーが5’-CAUTACTCACUTCCCCACUT-3’である。PCR産物の逆鎖におけるウラシル塩基をUSER酵素(New England Biolabs)によって除去した。続いて、1×NEBuffer2、0.64mMのS-アデノシルメチオニン及び20UのM.SssI CpGメチルトランスフェラーゼ(New England Biolabs)を含む溶液20μlにおいて、精製PCR産物100ngを2時間、37℃でメチル化してから、20分、熱不活化を65℃で行った。メチル化された222bpのモデルDNAをAMPure XPビーズによって精製した。
サンガーシーケンシングによる、TAPS、TAPSβ及びCAPSの検証用のモデルDNA:5mMのTris-Cl(pH7.5)、5mMのMgCl及び50mMのNaClを含むアニーリング緩衝液中で、5mC部位を1つ及び5hmC部位を1つ含む34bpのDNAオリゴを他のDNAオリゴとアニーリングしてから、400UのT4リガーゼ(NEB)を含む反応液中で、25℃において1時間ライゲーションし、1.8×AMPure XPビーズによって精製した。
ライゲーション反応後、USER酵素によってウラシルリンカーを除去して、
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCTCCCGACGCATGATCTGTACTTGATCGAChmCGTGCAACGATCGGAAGAGCACACGTCTGAACTCCAGTCACGCCAATATCTCGTATGCCGTCTTCTGCTTG
という配列(5’→3’)の最終生成物を得た。このモデルDNAの増幅用のPCRプライマーは、P5プライマーが5’-AATGATACGGCGACCACCGAG-3’、P7プライマーが5’-CAAGCAGAAGACGGCATACGAG-3’であった。
ポリメラーゼ試験及びサンガーシーケンシング用のモデルDNA:ポリメラーゼ試験及びサンガーシーケンシング用のモデルDNAは、異なるDNAオリゴを使用した以外は、上記と同じライゲーション方法で調製した。
最終生成物の配列(5’→3’)は、
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCTAGGTGCGCTAAGTTCTAGATCGCCAACTGGTTGTGGCCTTAGCAGTCTCGATCAGCTGCTACTGTACGTAGCATCTATAGCCGGCTTGCTCTCTCTGCCTCTAGCAGCTGCTCCCTATAGTGAGTCGTATTAACGATCGGAAGAGCACACGTCTGAACTCCAGTCACGCCAATATCTCGTATGCCGTCTTCTGCTTG
であった。モデルDNAを増幅するためのPCRプライマーは、上に示したP5プライマー及びP7プライマーである。プライマー伸長用のビオチン標識プライマー配列は、このP7プライマーの5’末端に連結したビオチンである。T7 RNAポリメラーゼによって転写した後のRT-PCR用のPCRプライマーは、このP5プライマー、及び5’-TGCTAGAGGCAGAGAGAGCAAG-3’というRTプライマーであった。
PCRバイアス試験用のモデルDNA:PCRバイアス試験用のモデルDNAは、異なるDNAオリゴを使用した以外は、上記と同じライゲーション方法で調製した。
最終生成物の配列(5’→3’)は、
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCTAGGTGCGCTAAGTTCTAGATCGCCAACTGGTTGTGGCCTTAGCAGTCTXGATCAGCTGCTACTGTACGTAGCATCTATAGCCGGCTTGCTCTCTCTGCCTCTAGCAGCTGCTCCCTATAGTGAGTCGTATTAACGATCGGAAGAGCACACGTCTGAACTCCAGTCACGCCAATATCTCGTATGCCGTCTTCTGCTTG
であり、この配列中、Xは、DHU、U、TまたはCである。このモデルDNAを増幅するためのPCRプライマーは、上記のP5プライマー及びP7プライマーである。
メチル化されたバクテリオファージλゲノムDNAの調製
Mg2+不含緩衝液(10mMのTris-Cl(pH8.0)、50mMのNaCl及び10mMのEDTA)中の0.64mMのSAM及び0.8U/μlのM.SssI酵素を含む反応液50μL中で、メチル化されていないバクテリオファージλDNA(Promega)1μgを2時間、37℃でメチル化した。続いて、M.SssI酵素0.5μL及びSAM1μLを加え、その反応物をさらに2時間、37℃でインキュベートした。その後、メチル化されたDNAを1×AMPure XPビーズで精製した。完全なメチル化を確認するために、すべての手順をNEB Buffer2において繰り返した。続いて、DNAのメチル化をHpaII消化アッセイで検証した。CutSmart緩衝液(NEB)中の2UのHpaII酵素(NEB)を含む反応液10μL中で、50ngのメチル化DNA及び非メチル化DNAを1時間、37℃で消化した。消化生成物を1%アガロースゲルに、未消化のλDNAコントロールとともに流した。そのアッセイ後、非メチル化λDNAは消化されたが、メチル化λDNAは、インタクトなままであったことから、CpGのメチル化が完全かつ成功したことが確認された。λDNAの配列は、GenBank-EMBLアクセッション番号J02459で見ることができる。
2kbの無修飾スパイクインコントロールの調製
pNIC28-Bsa4プラスミド(Addgene、カタログ番号26103)から、1ngのDNA鋳型、0.5μMのプライマー、1UのPhusion High-Fidelity DNA Polymerase(Thermo Fisher)を含む反応液中で、2kbのスパイクインコントロール(2kb-1、2、3)をPCR増幅した。PCRプライマー配列は、表2に列挙されている。
PCR産物をZymo-Spinカラムで精製した。2kbの無修飾コントロール配列(5’→3’)は、以下のとおりであった。
CACAGATGTCTGCCTGTTCATCCGCGTCCAGCTCGTTGAGTTTCTCCAGAAGCGTTAATGTCTGGCTTCTGATAAAGCGGGCCATGTTAAGGGCGGTTTTTTCCTGTTTGGTCACTGATGCCTCCGTGTAAGGGGGATTTCTGTTCATGGGGGTAATGATACCGATGAAACGAGAGAGGATGCTCACGATACGGGTTACTGATGATGAACATGCCCGGTTACTGGAACGTTGTGAGGGTAAACAACTGGCGGTATGGATGCGGCGGGACCAGAGAAAAATCACTCAGGGTCAATGCCAGCGCTTCGTTAATACAGATGTAGGTGTTCCACAGGGTAGCCAGCAGCATCCTGCGATGCAGATCCGGAACATAATGGTGCAGGGCGCTGACTTCCGCGTTTCCAGACTTTACGAAACACGGAAACCGAAGACCATTCATGTTGTTGCTCAGGTCGCAGACGTTTTGCAGCAGCAGTCGCTTCACGTTCGCTCGCGTATCGGTGATTCATTCTGCTAACCAGTAAGGCAACCCCGCCAGCCTAGCCGGGTCCTCAACGACAGGAGCACGATCATGCGCACCCGTGGGGCCGCCATGCCGGCGATAATGGCCTGCTTCTCGCCGAAACGTTTGGTGGCGGGACCAGTGACGAAGGCTTGAGCGAGGGCGTGCAAGATTCCGAATACCGCAAGCGACAGGCCGATCATCGTCGCGCTCCAGCGAAAGCGGTCCTCGCCGAAAATGACCCAGAGCGCTGCCGGCACCTGTCCTACGAGTTGCATGATAAAGAAGACAGTCATAAGTGCGGCGACGATAGTCATGCCCCGCGCCCACCGGAAGGAGCTGACTGGGTTGAAGGCTCTCAAGGGCATCGGTCGAGATCCCGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCGCTATCGGCTGAATTTGATTGCGAGTGAGATATTTATGCCAGCCAGCCAGACGCAGACGCGCCGAGACAGAACTTAATGGGCCCGCTAACAGCGCGATTTGCTGGTGACCCAATGCGACCAGATGCTCCACGCCCAGTCGCGTACCGTCTTCATGGGAGAAAATAATACTGTTGATGGGTGTCTGGTCAGAGACATCAAGAAATAACGCCGGAACATTAGTGCAGGCAGCTTCCACAGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCAGCCCACTGACGCGTTGCGCGAGAAGATTGTGCACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCT
120merのスパイクインコントロールの調製
120merのスパイクインコントロールをプライマー伸長法によって生成した。オリゴ配列及びプライマーは、表3に列挙されている。
簡潔に述べると、120mer-1のスパイクインでは、5mMのTris-Cl(pH7.5)、5mMのMgCl及び50mMのNaClを含むアニーリング緩衝液中で、3μMのオリゴを10μMのプライマーとアニーリングした。120mer-2のスパイクインでは、5μMのオリゴを7.5μMのプライマーとアニーリングした。0.4μMのdNTP(120mer-1:dATP/dGTP/dTTP/dhmCTP、120mer-2:dATP/dGTP/dTTP/dCTP)及び5UのKlenowポリメラーゼ(New England Biolabs)を含むNEB緩衝液2において、プライマー伸長を1時間、37℃で行った。反応後、スパイクインコントロールをZymo-Spinカラム(Zymo Research)で精製した。続いて、NEB緩衝液2中に0.64mMのSAM及び0.8U/μlのM.SssI酵素を含む反応液50μL中で、120merのスパイクインコントロールを2時間、37℃でメチル化し、Zymo-Spinカラムによって精製した。使用したすべてのスパイクイン配列は、https://figshare.com/s/80c3ab713c261262494bからダウンロードできる。
N5mCNN及びN5hmCNNを有する合成スパイクインの生成
N5mCNN及びN5hmCNNの配列を有する合成オリゴをアニーリング及び伸長方法によって生成した。オリゴ配列は、下記の表4に列挙されている。
簡潔に述べると、5mMのTris-Cl(pH7.5)、5mMのMgCl及び50mMのNaClを含むアニーリング緩衝液中で、10μMのN5mCNN及びN5hmCNNのオリゴ(IDT)を合わせてアニーリングした。0.4mMのdNTP(dATP/dGTP/dTTP/dCTP)及び5UのKlenowポリメラーゼ(NEB)を含むNEB緩衝液2中で、伸長を1時間、37℃で行った。反応後、スパイクインコントロールをZymo-Spinカラム(Zymo Research)で精製した。
N5mCNN及びN5hmCNNを有する合成スパイクイン(5’→3’):
GAAGATGCAGAAGACAGGAAGGATGAAACACTCAGGCGCACGCTGGCATNmCNNGACAAACCACAAGAACAGGCTAGTGAGAATGAAGGGATATGTTTGTAAGATGGTCNNGNATCTTGGGTTGTGTGGTGGATGTTGGCGTTGGTGGGTTTCAGAGTTGG
相補鎖(5’→3’):CCAACTCTGAAACCCACCAACGCCAACATCCACCACACAACCCAAGATNhmCNNGACCATCTTACAAACATATCCCTTCATTCTCACTAGCCTGTTCTTGTGGTTTGTCNNGNATGCCAGCGTGCGCCTGAGTGTTTCATCCTTCCTGTCTTCTGCATCTTC
DNAの消化及びHPLC-MS/MS解析
DNA試料を2UのNuclease P1(Sigma-Aldrich)及び10nMのデアミナーゼ阻害剤エリスロ-9-アミノ-β-ヘキシル-α-メチル-9H-プリン-9-エタノール塩酸塩(Sigma-Aldrich)で消化した。一晩37℃でインキュベート後、その試料をさらに、6Uのアルカリホスファターゼ(Sigma-Aldrich)及び0.5UのホスホジエステラーゼI(Sigma-Aldrich)で3時間、37℃で処理した。消化されたDNA溶液をAmicon Ultra-0.5mL 10K遠心フィルター(Merck Millipore)でろ過して、タンパク質を除去し、HPLC-MS/MS解析を行った。
HPLC-MS/MS解析は、6495B Triple Quadrupole Mass Spectrometer(Agilent)と接続した1290 Infinity LC Systems(Agilent)で行った。ZORBAX Eclipse Plus C18カラム(2.1×150mm、1.8マイクロメートル、Agilent)を使用した。そのカラム温度を40℃に保ち、溶媒系は、10mMの酢酸アンモニウムを含む水(pH6.0、溶媒A)及び水-アセトニトリル(60/40(v/v)、溶媒B)であり、流速は0.4mL/分であった。グラジエントは、溶媒B0%で0~5分、溶媒B0→5.63%で5~8分、溶媒B5.63%で8~9分、溶媒B5.63→13.66%で9~16分、溶媒B13.66→100%で16~17分、溶媒B100%で17~21分、溶媒B100→0%で21~24.3分、溶媒B0%で24.3~25分であった。定量のために、MSの動的多重反応モニタリングモード(dMRM)を使用した。供給源依存性パラメーターは、ガス温度230℃、ガス流束14L/分、ネブライザー40psi、シースガス温度400℃、シースガス流束11L/分、正イオンモードにおけるキャピラリー電圧1500V、ノズル電圧0V、高圧RF110V及び低圧RF80V(いずれも正イオンモード)であった。フラグメンター電圧は、すべての化合物において380Vであったとともに、他の化合物依存性パラメーターは、表5にまとめられているとおりであった。
NgTET1の発現及び精製
Hisタグ付加NgTET1タンパク質(GG739552.1)をコードするpRSET-Aプラスミドを設計し、Invitrogenから購入した。以前に説明されたものに、いくつかの改変を加えて(J.E.Pais et al.,Biochemical characterization of a Naegleria TET-like oxygenase and its application in single molecule sequencing of 5-methylcytosine.Proc.Natl.Acad.Sci.U.S.A.112,4316-4321(2015)。参照により、本明細書に援用される)、タンパク質をE.coli BL21(DE3)株で発現させて、精製した。簡潔に述べると、一晩、小規模培養して得た細菌を、タンパク質の発現のために、OD600が0.7~0.8になるまで、LB培地において37℃及び200rpmで成長させた。続いて、培養液を室温まで冷却し、標的タンパク質の発現を0.2mMのイソプロピル-β-d-1-チオガラクトピラノシド(IPTG)で誘導した。細胞を、さらに18時間、18℃及び180rpmで維持した。その後、細胞を回収し、20mMのHEPES(pH7.5)、500mMのNaCl、1mMのDTT、20mMのイミダゾール、1μg/mLのロイペプチン、1μg/mLのペプスタチンA及び1mMのPMSFを含む緩衝液に再懸濁させた。細胞をEmulsiFlex-C5高圧ホモジナイザーで破砕し、溶解液を遠心分離によって1時間、30,000×g及び4℃で清澄化した。回収した上清をNi-NTA樹脂に充填し、20mMのHEPES(pH7.5)、500mMのイミダゾール、2MのNaCl、1mMのDTTを含む緩衝液で、NgTET1タンパク質を溶出させた。続いて、回収した画分をHiLoad 16/60 Sdx 75(20mMのHEPES(pH7.5)、2MのNaCl、1mMのDTT)で精製した。次に、NgTET1を含む画分を回収し、20mMのHEPES(pH7.0)、10mMのNaCl、1mMのDTTを含む緩衝液に、緩衝液を交換し、HiTrap HP SPカラムに入れた。純粋なタンパク質を塩濃度勾配で溶出させ、回収し、20mMのTris-Cl(pH8.0)、150mMのNaCl及び1mMのDTTを含む最終緩衝液に、緩衝液を交換した。続いて、タンパク質を130μMまで濃縮し、グリセロール(30%(v/v))と混合し、アリコートを-80℃で保存した。
mTET1CDの発現及び精製
N末端Flagタグを有するmTET1CD触媒ドメイン(NM_001253857.2、4371~6392)をpcDNA3-Flagに、KpnI制限部位とBamH1制限部位の間にクローニングした。タンパク質の発現のために、1mgのプラスミドを1LのExpi293F(Gibco)細胞培養液に、1×10細胞/mLの密度でトランスフェクションし、細胞を48時間、37℃、170rpm及び5%COで成長させた。その後、細胞を遠心分離によって回収し、50mMのTris-Cl(pH=7.5)、500mMのNaCl、1×cOmplete Protease Inhibitor Cocktail(Sigma)、1mMのPMSF、1%Triton X-100を含む溶解緩衝液に再懸濁させ、氷で20分インキュベートした。続いて、細胞溶解液を遠心分離によって30分、30000×g及び4℃で清澄化した。回収した上清をANTI-FLAG M2 Affinity Gel(Sigma)で精製し、20mMのHEPES(pH=8.0)、150mMのNaCl、0.1mg/mLの3×Flagペプチド(Sigma)、1×cOmplete Protease Inhibitor Cocktail(Sigma)、1mMのPMSFを含む緩衝液で、純粋なタンパク質を溶出させた。回収した画分を濃縮し、20mMのHEPES(pH=8.0)、150mMのNaCl及び1mMのDTTを含む最終緩衝液に、緩衝液を交換した。濃縮したタンパク質をグリセロール(30%(v/v))と混合し、液体窒素で凍結し、アリコートを-80℃で保存した。組み換えmTET1CDの活性及び品質をMALDI質量分析解析によって調べた。このアッセイに基づくと、組み換えmTET1CDは、十分な活性を有し、5mCの5caCへの酸化を触媒できる。試験したモデルオリゴの、いくらかの有意な消化が、MALDIによって検出されたことから、タンパク質がヌクレアーゼを含まないことが確認された。
TETによる酸化
NgTET1による酸化:222bpのモデルDNAオリゴをTetによって酸化するために、50mMのMOPS緩衝液(pH6.9)、100mMの硫酸アンモニウム鉄(II)、1mMのa-ケトグルタル酸塩、2mMのアスコルビン酸、1mMのジチオスレイトール(DTT)、50mMのNaCl及び5μMのNgTETを含む溶液20μl中で、222bpのDNA100ngを1時間、37℃でインキュベートした。その後、0.4UのプロテイナーゼK(New England Biolabs)をその反応混合物に加え、30分、37℃でインキュベートした。その生成物をZymo-Spinカラム(Zymo Research)によって、メーカーの指示に従って精製した。
ゲノムDNAをNgTET1によって酸化するために、50mMのMOPS緩衝液(pH6.9)、100mMの硫酸アンモニウム鉄(II)、1mMのa-ケトグルタル酸塩、2mMのアスコルビン酸、1mMのジチオスレイトール、50mMのNaCl及び5μMのNgTET1を含む溶液50μl中で、ゲノムDNA500ngを1時間、37℃でインキュベートした。その後、4UのプロテイナーゼK(New England Biolabs)をその反応混合物に加え、30分、37℃でインキュベートした。その生成物を1.8×Ampureビーズで、メーカーの指示に従って精製した。
mTET1による酸化:50mMのHEPES緩衝液(pH8.0)、100μMの硫酸アンモニウム鉄(II)、1mMのa-ケトグルタル酸塩、2mMのアスコルビン酸、1mMのジチオスレイトール、100mMのNaCl、1.2mMのATP及び4μMのmTET1CDを含む反応液50μl中で、ゲノムDNA100ngを80分、37℃でインキュベートした。その後、0.8UのプロテイナーゼK(New England Biolabs)をその反応混合物に加え、1時間、50℃でインキュベートした。その生成物をBio-Spin P-30 Gel Column(Bio-Rad)及び1.8×AMPure XPビーズで、メーカーの指示に従って精製した。
ボランによる還元
Pic-BHによる還元:MeOH中の5Mのα-ピコリン-ボラン(pic-BH、Sigma-Aldrich)25μLと、3Mの酢酸ナトリウム溶液(pH5.2、Thermo Fisher)5μLをDNA試料20μLに加え、60℃で1時間インキュベートした。その生成物を、222bpの場合にはZymo-Spinカラム(Zymo Research)によって、メーカーの指示に従って精製し、またはオリゴの場合にはMicro Bio-Spin 6 Columns(Bio-Rad)によって、メーカーの指示に従って精製した。
あるいは、2-ピコリン-ボラン(pic-ボラン、Sigma-Aldrich)100mgをDMSO187μLに溶解して、3.26M前後の溶液を得た。各反応では、pic-ボラン溶液25μLと、3Mの酢酸ナトリウム溶液(pH5.2、Thermo Fisher)5μLをDNA試料20μLに加え、3時間、70℃でインキュベートした。その生成物を、ゲノムDNAの場合にはZymo-Spinカラムによって、DNAオリゴの場合にはMicro Bio-Spin 6 Columns(Bio-Rad)によって、メーカーの指示に従って精製した。
ピリジンボランによる還元:600mMの酢酸ナトリウム溶液(pH=4.3)及び1Mのピリジンボランを含む反応液50μL中で、水35μL中の酸化DNA50~100ngをEppendorf ThermoMixerにおいて16時間、37℃及び850rpmで還元した。その生成物をZymo-Spinカラムによって精製した。
単一ヌクレオシドのpic-ボラン反応:MeOH中の3.26Mの2-ピコリン-ボラン(pic-ボラン、Sigma-Aldrich)500μLと、3Mの酢酸ナトリウム溶液(pH5.2、Thermo Fisher)500μLを2’-デオキシシチジン-5-カルボン酸ナトリウム塩(Berry&Associates)10mgに加えた。その混合物を1時間、60℃で攪拌した。その生成物をHPLCによって精製して、純粋な化合物を白色泡状物として得た。高分解能MS (Q-TOF) m/z [M + Na]+ C14Naに対する計算値: 253.0800;実測値: 253.0789.
5hmCのブロック
50mMのHEPES緩衝液(pH8)、25mMのMgCl、200μMのウリジン二リン酸グルコース(UDP-Glc、New England Biolabs)、10UのβGT(Thermo Fisher)及び10μMの5hmC DNAオリゴを含む溶液20μl中で、5hmCのブロックを1時間、37℃で行った。その生成物をMicro Bio-Spin 6 Columns(Bio-Rad)によって、メーカーの指示に従って精製した。
5fCのブロック
100mMのMES緩衝液(pH5.0)、10mMのO-エチルヒドロキシルアミン(Sigma-Aldrich)及び10μMの5fC DNAオリゴにおいて、5fCのブロックを2時間、37℃で行った。その生成物をMicro Bio-Spin 6 Columns(Bio-Rad)によって、メーカーの指示に従って精製した。
5caCのブロック
75mMのMES緩衝液(pH5.0)、20mMのN-ヒドロキシスクシンイミド(NHS、Sigma-Aldrich)、20mMの1-(3-ジメチルアミノプロピル)-3-エチルカルボジイミド塩酸塩(EDC、Fluorochem)及び10μMの5caC DNAオリゴにおいて、5caCのブロックを37℃で0.5時間行った。続いて、Micro Bio-Spin 6 Columns(Bio-Rad)をメーカーの指示に従って用いて、緩衝液を100mMのリン酸ナトリウム(pH7.5)、150mMのNaClに交換した。10mMのエチルアミン(Sigma-Aldrich)をそのオリゴに加え、1時間、37℃でインキュベートした。その生成物をMicro Bio-Spin 6 Columns(Bio-Rad)によって、メーカーの指示に従って精製した。
5hmCの酸化
5hmC DNAオリゴ46μLを1MのNaOH2.5μLで30分、37℃で、振とう培養器において変性させてから、50mMのNaOH及び15mMの過ルテニウム酸カリウム(KRuO、Sigma-Aldrich)を含む溶液1.5μLで1時間、氷上で酸化させた。その生成物をMicro Bio-Spin 6 Columnsによって、メーカーの指示に従って精製した。
TaqαIアッセイによるTAPS変換の検証
TaqαI制限部位(TCGA)を含む標的領域をPCR増幅し、その後、TaqαI消化を行うことによって、TAPS後の5mCの変換を試験した。例えば、我々のTAPSライブラリーにおける5mCの変換は、CpGがメチル化されたλDNAスパイクインコントロールから増幅される単一のTaqαI制限部位を含む194bpのアンプリコンに基づき試験できる。194bpのアンプリコンから増幅されたPCR産物をTaqαI制限酵素で消化し、消化産物を2%アガロースゲルで調べる。未変換のコントロールDNAで増幅されたPCR産物は、TaqαIによって消化され、そのゲル上に、2本のバンドが示される。TAPSで変換された試料では、CからTへの変換により、制限部位は喪失するので、194bpのアンプリコンは、インタクトのままとなる。全体的な変換レベルは、消化されたゲルバンド及び消化されなかったゲルバンドの定量に基づき評価でき、成功したTAPS試料では、そのレベルは、95%超となるはずである。
簡潔に述べると、対応するプライマーを用いて、Taq DNAポリメラーゼ(New England Biolabs)によって、変換させたDNA試料をPCR増幅した。そのPCR産物を4単位のTaqαI制限酵素(New England Biolabs)とともに、1×CutSmart緩衝液(New England Biolabs)において、30分、65℃でインキュベートし、2%アガロースゲル電気泳動によって調べた。
定量ポリメラーゼ連鎖反応(qPCR)
TAPS前及びTAPS後のモデルDNA間で増幅曲線及び融解曲線を比較するために(図11)、1×LightCycler 480 High Resolution Melting Master Mix(Roche Diagnostics Corporation)、250nMのプライマー(フォワードプライマー:CCTGATGAAACAAGCATGTC及びリバースプライマー:CATTACTCACTTCCCCACTT)、ならびに3mMのMgSOを含むPCRマスターミックス19μLに、DNA試料1ngを加えた。PCR増幅のために、初期変性工程を10分、95℃で行ってから、95℃での変性5分、カスタマイズしたアニーリング温度でのアニーリング5分、及び72℃での伸長5分というサイクルを40サイクル行った。最終工程には、95℃で1分、70℃で1分を含め、融解曲線(0.02℃刻み、毎回取得前に5分置いた)は65℃~95℃とした。
他のアッセイでは、qPCRは、1×Fast SYBR Green Master Mix(Thermo Fisher)、200nMのフォワードプライマー及びリバースプライマーを含むPCRマスターミックス19μLに、所要量のDNA試料を加えることによって行った。PCR増幅のために、初期変性工程を20秒、95℃で行ってから、95℃での変性3秒、アニーリング20秒、及び60℃での伸長というサイクルを40サイクル行った。
HpaII-qPCRアッセイによる、mESC gDNAにおけるCCGGのメチル化レベルの検証
mESC gDNA1μgを、反応液50μL中の50単位のHpaII(NEB、50単位/μL)及び1×CutSmart緩衝液と16時間、37℃でインキュベートした。コントロールの反応では、HpaIIを加えなかった。プロテイナーゼK1μLを反応物に加え、40℃で30分インキュベートしてから、プロテイナーゼKを10分、95℃で不活化した。特異的なCCGGの位置について、対応するプライマーセット(表9に列挙されている)を用いて、上記のようなqPCRアッセイによって、HpaIIで消化した試料またはコントロール試料のCt値を測定した。
サンガーシーケンシング
エキソヌクレアーゼI及びエビ由来アルカリホスファターゼ(New England Biolabs)またはZymo-SpinカラムによってPCR産物を精製し、サンガーシーケンシング用に処理した。
異なる長さの断片に対するDNA損傷試験
mESCゲノムDNAに0.5%のCpGメチル化λDNAをスパイクインし、断片化しないままにするか、またはCovaris M220という装置で超音波処理し、AMPure XPビーズで、500~1kbまたは1kb~3kbでサイズ選別を行った。DNA200ngをmTET1CDで1回酸化させ、上記のようなピリジンボラン複合体で還元するか、またはEpiTect Bisulfite Kit(Qiagen)をメーカーのプロトコールに従って用いて変換した。TAPS及びバイサルファイトによる変換の前及び後のDNA10ngを1%アガロースゲルに流した。可視化するために、バイサルファイト変換を行ったゲルを10分、氷浴で冷やした。TAPS試料における5mCの変換を、上記のように、TaqαI消化アッセイによって試験した。
mESCの培養及びゲノムDNAの単離
ゼラチンコートプレートで、15%FBS(Gibco)、2mMのL-グルタミン(Gibco)、1%非必須アミノ酸(Gibco)、1%ペニシリン/ストレプトアビジン(Gibco)、0.1mMのβ-メルカプトエタノール(Sigma)、1000単位/mLのLIF(Millipore)、1μMのPD0325901(Stemgent)及び3μMのCHIR99021(Stemgent)を添加したダルベッコ改変イーグル培地(DMEM)(Invitrogen)中で、マウスESC(mESC)E14を培養した。培養液を37℃及び5%COで維持し、2日おきに継代した。
ゲノムDNAの単離のために、細胞を遠心分離によって5分、1000×g及び室温で回収した。Quick-DNA Plus Kit(Zymo Research)をメーカーのプロトコールに従って用いて、DNAを抽出した。
TAPS及びWGBS用のmESC gDNAの調製
全ゲノムバイサルファイトシーケンシング(WGBS)用に、mESC gDNAに、0.5%の非メチル化λDNAをスパイクインした。全ゲノムTAPS用には、mESC gDNAに、0.5%のメチル化λDNA及び0.025%の2kbの無修飾スパイクインコントロールをスパイクインした。DNA試料をCovaris M220という装置によって断片化し、AMPure XPビーズで、200~400bpでサイズ選別を行った。加えて、AMPure XPビーズによるサイズ選別後、TAPS用のDNAに、0.25%のN5mCNN及びN5hmCNNコントロールオリゴをスパイクインした。
全ゲノムバイサルファイトシーケンシング
全ゲノムバイサルファイトシーケンシング(WGBS)には、0.5%の非メチル化バクテリオファージλDNAをスパイクインした断片化mESC gDNA200ngを使用した。KAPA HyperPrepキット(Kapa Biosystems)をメーカーのプロトコールに従って用いて、末端修復反応及びA-テーリング反応、ならびにメチル化アダプター(NextFlex)のライゲーションを行った。その後、EpiTect Bisulfite Kit(Qiagen)をIlluminaのプロトコールに従って用いて、DNAに対してバイサルファイト変換を行った。KAPA Hifi Uracil Plus Polymerase(Kapa Biosystems)を用いて、最終的なライブラリーを6サイクル増幅し、1×Ampureビーズで精製した。NextSeq500シーケンサー(Illumina)で、NextSeq High Output Kitを15%のPhiXというコントロールライブラリースパイクインとともに用いて、WGBSシーケンシングライブラリーに対して、80bpのペアエンドシーケンシングを行った。
全ゲノムTAPS
全ゲノムTAPS用に、0.5%のメチル化λDNA及び0.025%の2kbの無修飾スパイクインコントロールをスパイクインした断片化mESC gDNA100ngを使用した。KAPA HyperPrepキットをメーカーのプロトコールに従って用いて、末端修復反応及びA-テーリング反応、ならびにIllumina Multiplexingアダプターのライゲーションを行った。ライゲーションしたDNAを、上記のプロトコールに従って、mTET1CDで2回酸化させてから、ピリジンボランで還元した。KAPA Hifi Uracil Plus Polymeraseを用いて、最終的なシーケンシングライブラリーを5サイクル増幅し、1×Ampureビーズで精製した。NextSeq500シーケンサー(Illumina)で、1つのNextSeq High Outputキットを用いて、1%のPhiXコントロールライブラリースパイクインとともに、全ゲノムTAPSシーケンシングライブラリーに対して、80bpのペアエンドシーケンシングを行った。
dsDNAライブラリー調製キットを用いたローインプット全ゲノムTAPS
上記のようにして、全ゲノムTAPS用に調製したmESC gDNAをローインプット全ゲノムTAPSに使用した。簡潔に述べると、NgTET1を1回、上記のプロトコールに従って用いて、mESC gDNAを100ng、10ng及び1ng含む試料を酸化させた。NEBNext Ultra II(New England Biolabs)またはKAPA HyperPrepキットをメーカーのプロトコールに従って用いて、末端修復反応及びA-テーリング反応、ならびにライゲーションを行った。その後、上記のようにして、DNAに対してpic-ボラン反応を行った。KAPA Hifi Uracil Plus Polymeraseを用いて、変換されたライブラリーを増幅し、1×Ampureビーズで精製した。
ssDNAライブラリー調製キットを用いたローインプット全ゲノムTAPS
上記のようにして、全ゲノムTAPS用に調製したmESC gDNAをローインプット全ゲノムTAPSに使用した。簡潔に述べると、mESC gDNAを100ng、10ng、1ng、100pg及び10pg含む試料をNgTET1で1回酸化させ、上記のようにして、pic-ボランで還元した。Accel-NGS Methyl-Seq DNA Library Kit(Swift Biosciences)をメーカーのプロトコールに従って用いて、シーケンシングライブラリーを調製した。KAPA Hifi Uracil Plus Polymeraseを用いて、最終的なライブラリーを6サイクル(100ng)、9サイクル(10ng)、13サイクル(1ng)、16サイクル(100pg)及び21サイクル(10pg)増幅し、0.85×Ampureビーズで精製した。
別の実験において、上記のようにして、全ゲノムTAPS用に調製したmESC gDNAをローインプット全ゲノムTAPSに使用した。簡潔に述べると、KAPA HyperPrepキットをメーカーのプロトコールに従って用いて、mESC gDNAを100ng、10ng及び1ng含む試料を末端修復反応及びA-テーリング反応に用いて、Illumina Multiplexingアダプターにライゲーションした。続いて、ライゲーションした試料を、上記のプロトコールに従って、mTET1CDで1回酸化させてから、ピリジンボランで還元した。KAPA Hifi Uracil Plus Polymeraseを用いて、変換されたライブラリーを5サイクル(100ng)、8サイクル(10ng)及び13サイクル(1ng)増幅し、1×AMPure XPビーズで精製した。
セルフリーDNA TAPS
10ng及び1ngのセルフリーDNA試料から、セルフリーDNA TAPS試料を調製した。簡潔に述べると、上記のようにして、試料をNgTET1で1回酸化させ、pic-ボランで還元した。Accel-NGS Methyl-Seq DNA Library Kit(Swift Biosciences)をメーカーのプロトコールに従って用いて、シーケンシングライブラリーを調製した。KAPA Hifi Uracil Plus Polymeraseを用いて、最終的なライブラリーを9サイクル(10ng)及び13サイクル(1ng)増幅し、0.85×Ampureビーズで精製した。
別の実験において、全ゲノムTAPSに関して上記されているようにして、セルフリーDNA TAPS試料を10ng及び1ngのセルフリーDNA試料から調製した。簡潔に述べると、KAPA HyperPrepキットをメーカーのプロトコールに従って用いて、セルフリーDNA試料を末端修復反応及びA-テーリング反応に使用して、Illumina Multiplexingアダプターにライゲーションした。続いて、ライゲーションした試料を、上記のプロトコールに従って、mTET1CDで1回酸化させてから、ピリジンボランで還元した。KAPA Hifi Uracil Plus Polymeraseを用いて、変換されたライブラリーを7サイクル(10ng)及び13サイクル(1ng)増幅し、1×AMPure XPビーズで精製した。
WGBSデータ処理
ペアエンドリードをFASTQとしてIllumina BaseSpaceからダウンロードし、その後、Trim Galore! v0.4.4(https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/)でクオリティトリミングを行った。トリミング後に、少なくとも1つのリードが35bpよりも短いリードペアを除去した。--no_overlapオプションを用いたBismark v0.19(F.Krueger,S.R.Andrews,Bismark:Bismark:a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 27,1571-1572(2011)、参照により本明細書に援用される)を使用して、mm9バージョンのマウスゲノム、λファージ及びPhiX(Illumina iGENOMESに由来する配列)を組み合わせたゲノムに、トリミングしたリードをマッピングした。「three-C」フィルターを用いて、非変換率が極端なリードを除去した。Picard v1.119(http://broadinstitute.github.io/picard/)のMarkDuplicatesを用いて、PCRデュプリケートをコールした。マッピングアーチファクトが発生しやすいことが知られている領域をダウンロードし(https://sites.google.com/site/anshulkundaje/projects/blacklists)、さらなる解析から除外した(E.P.Consortium,An integrated encyclopedia of DNA elements in the human genome.Nature 489,57-74(2012)。参照により、本明細書に援用される)。
TAPSデータの事前処理
ペアエンドリードをIllumina BaseSpaceからダウンロードし、その後、Trim Galore! v0.4.4でクオリティトリミングを行った。トリミング後に、少なくとも1つのリードが35bpよりも短いリードペアを除去した。BWA mem v.0.7.15(H.Li,R.Durbin,Fast and accurate short read alignment with Burrows-Wheeler transform.transform. Bioinformatics25,1754-1760(2009)。参照により、本明細書に援用される)をデフォルトパラメーターで用いて、スパイクイン配列、λファージ及びmm9バージョンのマウスゲノムを組み合わせたゲノムに、トリミングしたリードをマッピングした。マッピングアーチファクトが発生しやすいことが知られている領域をダウンロードし(https://sites.google.com/site/anshulkundaje/projects/blacklists)、さらなる解析から除外した(E.P.Consortium,Nature489,57-74(2012))。
TAPSにおいて変換された塩基の検出
カスタムのpython3スクリプト(MF-filter.py)を用いて、アラインメントしたリードを元の上鎖(OT)及び元の下鎖(OB)に分けた。続いて、OT及びOBで別々に、Picard MarkDuplicatesを用いて、PCRデュプリケートを除去した。デュプリケートが除去されたマッピング済みのOTリード及びOBリードで別々に、BamUtil clipOverlap(https://github.com/statgen/bamUtil)を用いて、リードペアにおける重複セグメントを除去した。続いて、samtools mpileup及びカスタムのpython3スクリプト(MF-caller_MOD.py)を用いて、改変塩基を検出した。
TAPS及びWGBSのシーケンシングクオリティ解析
python3スクリプト(MF-phredder.py)によって、Illumina BaseSpaceからダウンロードしたような元のFASTQファイルから、ヌクレオチドの種類ごとのクオリティスコア統計を抽出した。
TAPS及びWGBSのカバレッジ解析
Bedtools v2.25 genomecov(A.R.Quinlan,I.M.Hall,BEDTools:a flexible suite of utilities for comparing genomic features.Bioinformatics 26,841-842(2010)。参照により、本明細書に援用される)で、塩基ごとのゲノムカバレッジファイルを作成した。TAPS及びWGBSの間で、相対的なカバレッジ分布を比較するために、samtools viewの-sオプションを用いて、WGBSにおける対応するカバレッジ中央値に対して、TAPSリードをサブサンプリングした。WGBS及びサブサンプリング済みのTAPSにおけるカバレッジを比較する解析では、TAPS及びWGBSの両方のbamファイルで、clipOverlapを使用した。
TAPS及びWGBSによって測定されたシトシン修飾の解析
塩基ごとの修飾リードの割合を、それぞれBismarkのアウトプット及びMF-caller_MOD.pyのアウトプットから計算した。Bedtools intersectを用いて、積集合を求め、統計解析結果及び図をR及びMatlabで作成した。IGV v2.4.6(J.T.Robinson et al.,Integrative genomics viewer.Nat.Biotechnol.29,24-26(2011)。参照により、本明細書に援用される)を用いて、ゲノム領域を可視化した。CGI周囲のカバレッジ及び修飾レベルをプロットするために、mm9におけるすべてのCGI座標をUCSCゲノムブラウザーからダウンロードし、20個のウィンドウにビニングし、両側において、サイズが80bpのウィンドウ最大50個分だけ広げた(ただし、次のCGIまでの距離の半分に達しないことを条件とした)。Bedtools mapを用いて、各ビンにおいて、(CpGにおける)平均修飾レベル、及び(すべての塩基、両方の鎖における)カバレッジを計算した。各ビンの値を再び平均化し、その後、Matlabでプロットした。
データ処理時間のシミュレーション
(-p -ss NS50 --errfree --minQ 15 -k 0 -nf 0 -l 75 -c 1000000 -m 240 -s 0 -ir 0 -ir2 0 -dr 0 -dr2 0 -sam -rs 10というパラメーターで、)λファージゲノムに基づき、ART42を用いて、合成ペアエンドシーケンシングリードをシミュレートした。その後、すべてのCpGの位置の50%を、修飾されたものとしてマークし、カスタムのpython3スクリプトを用いて、TAPS(修飾塩基を変換する)またはWGBS(無修飾塩基を変換する)のいずれかとして、2つのライブラリーを作製した。その後、本書における各方法に使用するパイプラインに従って、それらのリードを処理した。Linuxのコマンド時間によって、処理時間を測定した。すべての解析工程は、1つのIntel Xeon CPU及び250GBメモリーを用いて、シングルスレッドモードで行った。
結果及び考察
pic-BHは、これまで知られていなかった還元的脱炭酸/脱アミノ反応によって、5fC及び5caCをDHUに容易に変換できることが見出された(図4)。その反応は、MALDIを用いて、単一のヌクレオシド及び複数のオリゴヌクレオチドの両方で、定量的であることが示された(図2~3及び6-1~7-2)。
11merの5caC含有DNAオリゴをモデルとして使用して、5caCと反応できる化学物質をスクリーニングした(マトリックス支援レーザー脱離/イオン化質量分析(MALDI)によってモニタリングした)。特定のボラン含有化合物が、5caCオリゴと効率的に反応して、分子量が41Da減少したことがわかった(図1及び2)。さらなる試験用に、ピリジンボラン及びその誘導体2-ピコリンボラン(pic-ボラン)を選択した。それらは、市販されているとともに、環境に優しい還元剤であるからである。
単一の5caCヌクレオシドでの反応を繰り返し、ピリジンボラン及びpic-ボランが、5caCをジヒドロウラシル(DHU)に変換することを確認した(図3、4B)。興味深いことに、ピリジンボラン及びpic-ボランは、どうやら還元的脱炭酸/脱アミノのような機序を通じて、5fCをDHUに変換することもわかった(図4C、6-1及び6-2)。両方の反応の詳細な機序は、まだ明らかになっていない。そのDNAオリゴでのボラン反応に対する、HPLC-MS/MSによる定量解析から、pic-ボランが、98%前後の効率性で、5caC及び5fCをDHUに変換するとともに、非メチル化シトシン、5mCまたは5hmCに対する活性を有さないことが確認される(図2B)。
ウラシル誘導体として、DHUは、DNAポリメラーゼ及びRNAポリメラーゼの両方によって、チミンとして認識されることができる。したがって、ボランによる還元を用いて、5caCのTへの変換及び5fCのTへの変換の両方を誘導できるとともに、ボランによる還元を、5fC及び5caCの、塩基単位の分解能でのシーケンシングに用いることができ、我々は、このシーケンシングをピリジンボランシーケンシング(「PS」)と称した(表6)。5fCTへのTへの、ボランによる還元は、ヒドロキシルアミンの結合(C.X.Song et al.,Genome-wide profiling of 5-formylcytosine reveals its roles in epigenetic priming.Cell 153,678-691(2013)。参照により、本明細書に援用される)を通じて、5caCのTへの、ボランによる還元は、EDCのカップリング(X.Lu et al.,Chemical modification-assisted bisulfite sequencing(CAB-Seq)for 5-carboxylcytosine detection in DNA.J.Am.Chem.Soc.135,9315-9317(2013)。参照により、本明細書に援用される)を通じてブロックできる(図6-1及び図6-2)。このブロッキングにより、PSを用いて、5fCまたは5caCを特異的にシーケンシング可能になる(表6)。
さらに、Tetアシストピリジンボランシーケンシング(「TAPS」)という本発明のプロセスでは、TET酵素を用いて、5mC及び5hmCを酸化させて5caCにし、その後、5caCをボランによって還元することができる(図5A~B、表1)。TAPSは、5mC及び5hmCのC→T変換を誘導できるので、5mC及び5hmCを塩基単位の分解能で検出するのに使用できる。
加えて、本明細書でTAPSβと称するプロセスでは、β-グルコシルトランスフェラーゼ(βGT)は、5hmCをグルコースで標識することによって、5hmCがTETによって酸化されないように(M.Yu et al.,Base-resolution analysis of 5-hydroxymethylcytosine in the mammalian genome.Cell 149,1368-1380(2012))、かつボランによって還元されないように保護でき(図7-1及び図7-2)、5mCのみを選択的にシーケンシング可能になる(図5B、表1)。その後、TAPSの測定結果からTAPSβの結果を差し引くことによって、5hmC部位を推定できる。あるいは、TETの代わりに、これまで酸化バイサルファイトファイトシーケンシング(oxBS)で用いてきた試薬である過ルテニウム酸カリウム(KRuO)(M.J.Booth et al.,Quantitative Sequencing of 5-Methylcytosine and 5-Hydroxymethylcytosine at Single-Base Resolution.Science 336,934-937(2012))を化学的酸化剤として用いて、5hmCを特異的に酸化させて、5fCにすることができる(図7-1及び図7-2)。このアプローチは、本明細書では、化学アシストピリジンボランシーケンシング(「CAPS」)と称し、このアプローチを用いて、5hmCを特異的にシーケンシングできる(図5B、表1)。したがって、TAPS及び関連方法は原則的に、シトシンの4つのエピジェネティックな修飾をすべてシーケンシングするための網羅的なシステムを提供できる(図5B、表1、表6)。
TAPS単独では、ゲノムにおける既存の5fC及び5caCも検出されることになる。しかしながら、通常の条件下では、ゲノムDNAにおける5fC及び5caCのレベルが極めて低いことを考慮すると、上記の事象は許容されるであろう。特定の条件下で、5fC及び5caCのシグナルを完全に除去したい場合には、ヒドロキシルアミンの結合によって5fCを保護し、EDCのカップリングによって5caCを保護することによって、容易に除去して、DHUへの変換を防ぐこともできる。
5mC及び5hmCを塩基単位の分解能でマッピングするための現行の標準的かつ最も広く用いられている方法であるバイサルファイトシーケンシングと比較して、TAPSの性能を評価した。Naegleria TET様オキシゲナーゼ(NgTET1)及びマウスTet1(mTet1)は両方とも、in vitroにおいて、5mCを効率的に酸化させて、5caCにできるので、これらを使用した。5mCのTへの変換を確認するために、完全にメチル化されたCpG部位を含むモデルDNAにTAPSを適用したところ、制限酵素消化(図8A~B)及びサンガーシーケンシング(図9A)によって実証されたように、TAPSが、5mCをTに効果的に変換できることが示された。TAPSβ及びCAPSも、サンガーシーケンシングによって検証した(図12)。
TAPSは、マウス胚性幹細胞(mESC)に由来するゲノムDNA(gDNA)にも適用した。HPLC-MS/MSによる定量により、予想どおり、5mCが、mESC gDNAにおけるシトシン修飾の98.5%を占め、残りが、5hmC(1.5%)、ならびに微量の5fC及び5caCで構成され、DHUを有さないことが示された(図9B)。TETによる酸化後、シトシン修飾の約96%が酸化されて、5caCになり、3%が酸化されて、5fCになった(図9B)。ボランによる還元後、シトシン修飾の99%超が、DHUに変換された(図9B)。これらの結果から、TETによる酸化及びボランによる還元の両方が、ゲノムDNAに対して効率的に機能することが示されている。
TETによる酸化及びボランによる還元の両方とも、穏和な反応であり、バイサルファイトと比べて、顕著なDNA分解は見られず(図10A~D)、これにより、DNAの回収率が高くなる。バイサルファイトシーケンシングを上回る別の顕著な利点は、TAPSが、非破壊的であり、DNAを最長で10kbの長さに維持できる点である(図10C)。さらに、DNAは、TAPS後、二本鎖を保ち(図10A~C)、その変換は、DNAの長さの影響を受けない(図15A~B)。
加えて、DHUは、天然の塩基に近いので、様々なDNAポリメラーゼ、及び等温DNAポリメラーゼまたはRNAポリメラーゼと適合し(図13A~B)、PCR中に、T/Cと比べて、バイアスが見られない(図14)。
mESC gDNAの2つの試料に対して、全ゲノムシーケンシングを行ったが、比較のために、一方は、TAPSを用いて変換し、もう一方は、標準的な全ゲノムバイサルファイトシーケンシング(WGBS)を用いて変換した。
TAPSの正確性を評価するために、完全に無修飾であるか、(上記の方法を用いて)CpGメチルトランスフェラーゼ(M.SssI)またはGpCメチルトランスフェラーゼ(M.CviPI)を用いてin vitroでメチル化したかのいずれかである異なる長さのスパイクインを付加した。5mC及び5hmC含む短いスパイクイン(120mer-1及び120mer-2)では、CpGコンテクスト及び非CpGコンテクストの両方の状況において、両方の鎖に対する両方の修飾に関して、ほぼ完全な変換が観察された(図17A~B)。
WGBSにおけるgDNA200ngに対して、TAPSには、gDNA100ngを使用した。TAPSの正確性を評価するために、3つの異なる種類のスパイクインコントロールを付加した。すべてのCpGが完全にメチル化されたλDNAを使用して、偽陰性率(5mCの非変換率)を推定し、2kbの無修飾アンプリコンを使用して、偽陽性率(無修飾Cの変換率)を推定し、いずれかの他の塩基(メチル化CはN5mCNN、ヒドロキシメチル化CはN5hmCNN)に囲まれたメチル化C及びヒドロキシメチル化Cの両方を含む合成オリゴスパイクインを使用して、異なる配列コンテクストにおける5mC及び5hmCでの変換率を比較した。mTet1及びピリジンボランを組み合わせた場合に、5mC変換率が最も高くなり(λスパイクインでは96.5%、合成スパイクインでは97.3%)、無修飾Cの変換率が最も低くなった(0.23%)(図18A~B及び図16)。2.7%~3.5%の偽陰性率、及びわずか0.23%の偽陽性率は、バイサルファイトシーケンシングに匹敵し、最近の研究では、9個の市販バイサルファイトキットの偽陰性率の平均は1.7%であり、偽陽性率は0.6%であった(Holmes,E.E.et al.Performance evaluation of kits for bisulfite-conversion of DNA from tissues,cell lines,FFPE tissues,aspirates,lavages,effusions,plasma,serum,and urine.PLoS One 9,e93933(2014))。合成スパイクインにより、TAPSが、5mC及び5hmCの両方に対して十分に機能すること、ならびにTAPSが、非CpGコンテクストにおいては、多少劣るのみであることが示唆されている。5hmCにおける変換は、5mCよりも8.2%低く、非CpGコンテクストにおける変換は、CpGコンテクストよりも11.4%低い(図18A)。
WGBSデータでは、アラインメント工程及び修飾コール工程の両方において、特別なソフトウェアが必要となる。これに対して、我々の処理パイプラインでは、標準的なゲノムアライナー(bwa)の後に、我々が「asTair」と呼ぶカスタムの修飾コールツールを使用する。シミュレーションしたWGBSリード及びTAPSリード(同じ部分メチル化ソース配列に由来する)を処理したところ、TAPS/asTairは、WGBS/Bismarkよりも3倍速かった(図18C)。
ほぼすべてのシトシンがチミンに変換されることにより、WGBSライブラリーは、極めて歪んだヌクレオチド組成を特徴とし、この組成は、Illuminaシーケンシングに悪影響を及ぼし得る。その結果、WGBSリードでは、TAPSと比べて、シトシン/グアニン塩基対において、実質的に、より低いシーケンシングクオリティスコアが見られた(図18E)。ヌクレオチド組成バイアスを相殺するために、少なくとも10~20%のPhiX DNA(塩基バランスの取れたコントロールライブラリー)を一般的に、WGBSライブラリーに加える(例えば、IlluminaのWhole-Genome Bisulfite Sequencing on the HiSeq 3000/HiSeq 4000 Systemsを参照されたい)。したがって、我々は、WGBSライブラリーに15%PhiXを添加した。この結果、BSで変換されたリードの情報内容の低下、及びバイサルファイト処理による、DNAの分解と組み合わさって、WGBSにおけるマッピング率は、TAPSと比べて有意に低下した(図18D及び表7)。
すなわち、同じシーケンシングコスト(NextSeq High Outputで1回)において、TAPSの平均深度は、WGBSの平均深度を上回った(TAPSでは21倍、WGBSでは13.1倍。表8)。さらに、TAPSでは、カバーされなかった領域はわずかで、全体では、WGBSと同じシーケンシング深度にダウンサンプリングした後でも、より均一なカバレッジ分布が見られた(四分位範囲:TAPSでは9、WGBSでは11。図19A及び表8)。
例えば、特に、CpGアイランド(CGI)は、WGBSとTAPSの間でシーケンシング深度の差について調整した時でも、概ね、TAPSの方が多くカバーされ(図21A)、その一方で、いずれにおいても、CGI内の脱メチル化は同程度であった(図22)。さらに、WGBSでは、高度にカバーされたCpG部位において、修飾レベルが低下したわずかなバイアスが見られたが(図23A)、我々の結果により、TAPSでは、修飾-カバレッジのバイアスがほんのわずかであることが示唆されている(図23B)。これらの結果から、TAPSでは、WGBSと比べて、シーケンシングクオリティが劇的に向上すると同時に、シーケンシングコストが効果的に半減されることが示されている。
TAPSのゲノムカバレッジの向上及び均一化により、少なくとも3本のリードによってカバーされたCpG部位の数が増加した。TAPSでは、マウスゲノムにおける合計43,205,316個のCpG部位の88.3%が、このレベルでカバーされたのに対し、WGBSでは、わずか77.5%であった(図21B及び19B)。TAPS及びWGBSにおいて、染色体領域全体にわたって、メチル化測定値の相関性が高くなった(図21D及び図20)。ヌクレオチドごとに基づくと、いずれの方法においても、少なくとも3本のリードによって、32,755,271個のCpG位置がカバーされた(図21B)。我々は、これらの部位内で、「修飾CpG」を、修飾レベルが少なくとも10%であるすべてのCpG位置として定義した(L.Wen et al.,Whole-genome analysis of 5-hydroxymethylcytosine and 5-methylcytosine at base resolution in the human brain.Genome Biology 15,49(2014))。この閾値を用いたところ、TAPSとWGBSの間で、CpGの95.8%において、修飾状態が一致した。少なくとも3本のリードによってカバーされたとともに、WGBSにおいて、修飾されていることが見出されたすべてのCpGの98.5%が、TAPSによって、修飾されたものとしてリコールされたことから、WGBSとTAPSの間での一致が良好であることが示された(図21C)。WGBS及びTAPSの両方で、少なくとも3本のリードによってカバーされたCpGごとに修飾レベルを比較したところ、TAPSとWGBSの間で良好な相関性が観察された(ピアソンのr=0.63、p<2e-16、図21E)。注目すべきことに、TAPSでは、高度に修飾されたCpG位置であって、WGBSによっては見落とされたCpG位置のサブセットが特定された(図21E、右下隅)。我々はさらに、直交性制限消化処理及びリアルタイムPCRアッセイを用いて、これらのCpGのうちの7個を検証し、それらのすべてが、完全にメチル化及び/またはヒドロキシメチル化されていることを確認した(表9)。
勘案すると、これらの結果によって、TAPSは、WGBSに直接置き換わることができ、実際に、WGBSよりも網羅的にメチロームが調べられることが示されている。
最後に、インプットの少ないDNAで、TAPSを試験したところ、TAPSが、1ngほどの少ないgDNAでも機能し、場合によっては、1細胞レベルに近い10pgのgDNAまで減少しても機能することが示された。TAPSは、1ngの循環セルフリーDNAでも有効に機能する。これらの結果から、インプットの少ないDNA及び臨床用途でのTAPSの潜在能力が示されている(図24A~C、図25A~B)。
1~2mlの血漿から採取した1つの健全な試料、1つのバレット食道(バレット)試料及び1つの膵臓癌試料に由来する3つの循環セルフリーDNA試料(cfDNA)で、TAPSを試験した。標準的なTAPSプロトコールは、下記のとおりであり、各試料を約10倍のカバレッジでシーケンシングした。cfDNA TAPSの結果を解析したところ、TAPSにより、高い5mC変換率(図26A)、低い偽陽性率(無修飾シトシンの変換、図26B)、高いマッピング率(図26C)及び低いPCR重複率(図26D)を含め、インプットの少ないcfDNAから、バルクなゲノムDNAの場合と同じ高品質なメチロームシーケンシングが行われることが示された。これらの結果から、cfDNAからの疾患診断に関して、TAPSに検出力があることが示されている。
TAPSは、メチル化を、C→Tの遺伝的バリアントまたは一塩基多型(SNP)から識別することもできるので、遺伝的バリアントを検出することができる。TAPSでは、メチル化及びC→TのSNPでは、異なるパターンが得られる。すなわち、メチル化では、元の上鎖(OT)/元の下鎖(OB)においてはT/Gリード、OTと相補的な鎖(CTOT)及びOBと相補的な鎖(CTOB)においてはA/Cリードが得られるのに対して、C→TのSNPでは、OT/OB及び(CTOB/CTOT)においてT/Aリードが得られる(図27)。これにより、1回の実験及びシーケンシングランで、メチル化情報及び遺伝的バリアント、すなわち変異の両方を得る際のTAPSの有用性がさらに向上する。本明細書に開示されているTAPS法のこの能力により、ゲノム解析とエピジェニック解析が統合され、標準的な全ゲノムシーケンシング(WGS)を行う必要がなくなることによって、シーケンシングコストが実質的に削減される。
要約すると、我々は、バイサルファイトを使用せずに、シトシンのエピジェニックな修飾を塩基単位の分解能でシーケンシングするPS由来の一連の方法を開発し、全メチロームシーケンシングにおけるTAPSの有用性を実証した。TAPSは、穏和な酵素反応及び化学反応を用いて、塩基単位の分解能で、高い感度及び特異性で、無修飾シトシンに影響を及ぼすことなく、5mC及び5hmCを直接検出することによって、半分のシーケンシングコストで、高品質かつより完全なメチロームをもたらすという点で、バイサルファイトシーケンシングよりも優れる。したがって、TAPSは、DNAのメチルシトシン及びヒドロキシメチルシトシンの解析における新たなスタンダードとして、バイサルファイトシーケンシングに置き換わる可能性がある。最近報告された、バイサルファイトを使用しない5fCシーケンシング方法(B.Xia et al.,Bisulfite-free,base-resolution analysis of 5-formylcytosine at the genome scale.Nat.Methods 12,1047-1050(2015)、C.Zhu et al.,Single-Cell 5-Formylcytosine Landscapes of Mammalian Early Embryos and ESCs at Single-Base Resolution.Cell Stem Cell 20,720-731(2017))において、シトシンに大規模な改変を導入するのとは異なり、TAPSでは、修飾シトシンを、天然塩基に近いDHUに変換し、そのDHUを一般的なポリメラーゼによって、Tとして「読み取る」ことができ、TAPSは、PCRを使用しないDNAシーケンシングと適合する可能性がある。TAPSは、様々な下流解析(パイロシーケンシング、メチル化感受性PCR、制限消化処理、MALDI質量分析、マイクロアレイ及び全ゲノムシーケンシングが挙げられるが、これらに限らない)と適合する。TAPSでは、長いDNAを保持できるので、マッピングしにくい特定の領域を調べるために、SMRTシーケンシング及びナノポアシーケンシングのようなロングリードシーケンシング技術と組み合わせると、極めて有用であり得る。プルダウン法をTAPSと組み合わせて、シーケンシングコストをさらに削減し、塩基単位の分解能の情報を低分解能のアフィニティーベースマップに加えることも可能である。本実施例では、TAPSが、コスト、複雑さ、及び解析に必要な時間を低減しながら、常法の用途で、WGBSに直接置き換わる可能性があることが示された。これにより、学術的な研究及び臨床診断において、エピジェネティックな解析がさらに広く採用可能になるであろう。
実施例2:エンドヌクレアーゼ濃縮TAPS(eeTAPS)
方法
スパイクインコントロールの調製
1ngのDNA鋳型、0.5μMのプライマー及び1×Phusion High-Fidelity PCR Master Mix with HF Buffer(Thermo Scientific)を含む反応液中で、pNIC28-Bsa4プラスミド(Addgene、カタログ番号26103)をPCR増幅することによって、4kbのスパイクインコントロールを調製した。プライマー配列は、表10に列挙されている。そのPCR産物をZymo-ICカラム(Zymo Research)によって精製し、反応液50μL中で、HpaII Methyltransferase(New England Biolabs)によって2時間、37℃でメチル化した。1×AMPure XPビーズ(Beckman Coulter)をメーカーのプロトコールに従って用いて、メチル化された生成物を精製した。以前に説明されたようにして(Wu,H.,Wu,X.J.and Zhang,Y.(2016)Base-resolution profiling of active DNA demethylation using MAB-seq and caMAB-seq.Nat Protoc,11:1081-1100)、非メチル化λ-DNA(Promega)をM.SssI酵素(New England Biolabs)でメチル化することによって、CpGが完全にメチル化されたλ-DNAを調製した。
担体DNAの調製
1ngのDNA鋳型、0.5μMのプライマー及び1×Phusion High-Fidelity PCR Master Mix with HF Buffer(Thermo Scientific)を含む反応液中で、pNIC28-Bsa4プラスミド(Addgene、カタログ番号26103)をPCR増幅することによって、担体DNAを調製した。プライマー配列は、表10に列挙されている。そのPCR産物をZymo-ICカラム(Zymo Research)によって精製し、Covaris M220によって断片化し、0.9×AMPure XPビーズで精製して、200~500bpの断片について選択した。
mESCの培養及びゲノムDNAの単離
ゼラチンコートプレートで、15%FBS(Gibco)、2mMのl-グルタミン(Gibco)、1%非必須アミノ酸(Gibco)、1%ペニシリン/ストレプトアビジン(Gibco)、0.1mMのβ-メルカプトエタノール(Sigma)、1,000単位ml-1の白血病阻止因子(Millipore)、1μMのPD0325901(Stemgent)及び3μMのCHIR99021(Stemgent)を添加したDMEM(Invitrogen)中で、E14 mESCを培養した。培養液を37℃及び5%CO2で維持し、2日おきに継代した。ゲノムDNAの単離のために、細胞を遠心分離によって5分、1000g及び室温で回収した。Quick-DNA Plus Kit(Zymo Research)をメーカーのプロトコールに従って用いて、DNAを抽出した。
mTet1CDの発現及び精製
mTet1触媒ドメイン(mTet1CD)の発現及び精製は、上記のようにして行った。
mTet1CDによる酸化
50mMのHEPES緩衝液(pH8.0)、100μMの硫酸アンモニウム鉄(II)、1mMのα-ケトグルタル酸塩、2mMのアスコルビン酸、1mMのジチオスレイトール、100mMのNaCl、1.2mMのATP及び4μMのmTet1CDを含む反応液50μl中で、0.5%のメチル化λ-DNA及び0.025%の2kbの無修飾DNAコントロールをスパイクインしたmESC gDNA200ngを80分、37℃で酸化させた。その後、0.8UのプロテイナーゼK(New England Biolabs)をその反応混合物に加え、1時間、50℃でインキュベートした。その生成物をBio-Spin P-30 Gel Column(Bio-Rad)及び1.8×AMPure XPビーズで、メーカーの指示に従って精製した。
DHUを消化するエンドヌクレアーゼのスクリーニング
mESC gDNA1μgを、上記のようにして、mTet1CDによって酵素的に酸化させた。その後、600mMの酢酸ナトリウム溶液(pH4.3)及び1Mのピリジンボランを含む反応液50μl中で、水35μl中の酸化DNAをEppendorf ThermoMixerにおいて、16時間、37℃及び850r.p.m.で還元した。Zymo-Spinカラムを用いて、その生成物を精製した。続いて、TAPSで変換したDNA、または未変換のDNA40ngを、USER(カタログ番号M5505S)、エンドヌクレアーゼIV(カタログ番号M0304S)、TmaエンドヌクレアーゼIII(カタログ番号M0291S)、エンドヌクレアーゼV(カタログ番号M0305S)、UDG(カタログ番号M0280S)、TthエンドヌクレアーゼIV(カタログ番号M0294S)、Fpg(カタログ番号M0240S)、エンドヌクレアーゼIII(Nth)(カタログ番号M0268S)、エンドヌクレアーゼVIII(カタログ番号M0299S)、APE1(カタログ番号M0282S)という酵素(すべてNew England Biolabs製)によって、メーカーのプロトコールに従って消化した。消化産物を1.8×AMPure XPビーズで、メーカーの指示に従って精製し、各生成物10ngを2%アガロースゲルに流した。
eeTAPS
mESCゲノムDNA(200ng、50ng、10ngまたは1ng)に、CCGGの配列コンテクストにおいてメチル化した4kbの0.05%コントロールをスパイクインし、上記のようにして、mTet1CDによって酸化させた。その後、600mMの酢酸ナトリウム溶液(pH4.3)及び1Mのピリジンボランを含む反応液50μl中で、水35μl中の酸化DNA試料をEppendorf ThermoMixerにおいて、16時間、37℃及び850r.p.m.で還元した。Zymo-Spinカラムを用いて、その生成物を精製した。CutSmart緩衝液中の2UのUSER酵素(New England Biolabs)を含む反応液20μL中で、変換した試料を1時間、37℃で消化し、0.35×~1×AMPure XPビーズでサイズ選別を行った。KAPA HyperPrepキットをメーカーのプロトコールに従って用いて、末端修復反応及びA-テーリング反応、ならびにIllumina Multiplexingアダプターのライゲーションを行った。コントロールライブラリーを調製するために、スパイクインコントロールを含む未変換のmESC gDNA200ngを、上記のようにして、USER酵素によって消化し、サイズ選別を行い、ライブラリーの構築に使用した。KAPA HiFi HotStart ReadyMixを用いて、最終的なシーケンシングライブラリーを6サイクル(200ngのインプットの場合)、8サイクル(50ngのインプットの場合)、10サイクル(10ngのインプットの場合)または14サイクル(1ngのインプットの場合)増幅し、0.35×~1×AMPure XPビーズでサイズ選別を行った。NextSeq500シーケンサー(Illumina)で、他のシーケンシングライブラリーと併せて、最終的なライブラリーに対して、80bpのペアエンドシーケンシングを行った。
rrTAPS
mESC gDNA1μgに、CpGがメチル化された1%のλをスパイクし、Fast digest Msp1酵素(Thermo Scientific)によって、反応液50μL中で30分、37℃で消化した。消化されたDNAをフェノール/クロロホルム沈殿法によって精製した。NEBNext(登録商標)Ultra(商標)II DNA Library Prep Kitをメーカーのプロトコールに従って用いて、末端修復反応及びA-テーリング反応、ならびにIllumina Multiplexingアダプターのライゲーションを行った。続いて、ライゲーションしたライブラリーを1.6×AMPure XPビーズで精製し、1%アガロースゲルに流した。100~400bpのDNA断片を切り出し、Monarch(登録商標)DNA Gel Extraction Kitによって、メーカーのプロトコールに従って精製した。アダプターにライゲーションした試料を担体DNA100ngにスパイクし、上記のようにして、mTet1CDによって二重酸化させた。600mMの酢酸ナトリウム溶液(pH4.3)及び1Mのピリジンボランを含む反応液50μl中で、水35μl中の酸化DNAをEppendorf ThermoMixerにおいて、16時間、37℃及び850r.p.m.で還元した。Zymo-Spinカラムを用いて、その生成物を精製した。最終的なシーケンシングライブラリーをKAPA HiFi ウラシル(+)Master Mixで6サイクル増幅し、1×AMPure XPビーズで精製した。NextSeq 500シーケンサー(Illumina)で、他のシーケンシングライブラリーと併せて、最終的なライブラリーに対して、80bpのペアエンドシーケンシングを行った。
eeTAPSのデータ解析
生のシーケンシングリードをTrimGalore(https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/)で処理して、--paired--length35というパラメーターで、アダプター及び品質のトリミングを行った。bwa mem 0.7.17-r1188(Li,H.and Durbin,R.(2009)Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics,25,1754-1760)をデフォルトパラメーターで用いて、クリーンしたリードをアラインメントした。4kbのモデルDNAでは、2,627~6,911のpNIC28-Bsa4配列を参照として使用した。mESC gDNAでは、mm9ゲノムを参照として用いた。適正にマッピングされたリードペアのみ(83または99として割り当てられたフラッグを有するリード1)を抽出して、bedtools v2.27.1(Quinlan,A.R.and Hall,I.M.(2010) BEDTools: a flexible suite of utilities for comparing genomic features.Bioinformatics,26,841-842)を用いて、その断片全体のエンドポイント及びリードスルーの両方でカバレッジを計算し、切断の割合を計算するときには、非切断部位も考慮に入れた。eeTAPSを解析するための詳しい計算パイプラインは、https://gitlab.com/jfeicheng/userenrichに見ることができる。eeTAPSにおいて、2つのテクニカルレプリケートをシーケンシングした。シーケンス深度がeeTAPSに及ぼす作用を解析するときに、2つのレプリケートに由来するアラインメントファイルをマージしてから、samtools view(Li,H.,Handsaker,B.,Wysoker,A.,Fennell,T.,Ruan,J.,Homer,N.,Marth,G.,Abecasis,G.,Durbin,R.and Genome Project Data Processing,S.(2009)The Sequence Alignment/Map format and SAMtools.Bioinformatics,25,2078-2079)を用いて、0.1~1の割合によってサブサンプリングした。
rrTAPSのデータ解析
生のシーケンシングリードをseqtk(https://github.com/lh3/seqtk)のtrimfq -b 2で処理して、各リードの左から2bpをトリミングした。Astair3.2.7を使用して、rrTAPSを処理した(8)。mm9ゲノムを参照として使用して、astair alignを用いて、クリーンしたリードをアラインメントした。メチル化CpGをastairコールによって抽出した。
mESCにおけるwgTAPS、eeTAPS及びrrTAPSの比較
wgTAPSデータをGSE112520(Liu,Y.B.,Siejka-Zielinska,P.,Velikova,G.,Bi,Y.,Yuan,F.,Tomkova,M.,Bai,C.S.,Chen,L.,Schuster-Bockler,B.and Song,C.X.(2019)Bisulfite-free direct detection of 5-methylcytosine and 5-hydroxymethylcytosine at base resolution.Nat Biotechnol,37,424-429)からダウンロードした。少なくとも4本のリードでカバーされたCpG部位のみを、カバーされたCpG部位とみなした。CpGメチル化レベルが、すべてのCpGメチル化レベルの第1四分位(wgTAPSでは0.5、eeTAPSでは0.28)よりも高いという基準に従って、メチル化CpG部位の数を定義した。bedtoolsを用いて、ゲノムを、重複しない100kbのウィンドウに分割した。CpG island trackをhttp://hgdownload.soe.ucsc.edu/goldenPath/mm9/database/cpgIslandExt.txt.gzからダウンロードした。Gene annotation fileをttp://hgdownload.soe.ucsc.edu/goldenPath/mm9/database/refGene.txt.gzからダウンロードした。wgTAPSでは、メチル化の平均値を用いて、各ウィンドウにおいて、メチル化を割り当てた。eeTAPSでは、切断割合が0.28超のCpG部位を、メチル化された部位とし、このカットオフ未満の部位は、非メチル化部位とし、そして、各ビンのメチル化レベルをメチル化CpG数/(メチル化CpG数+非メチル化CpG数)として測定した。e14 mESC細胞株の発現データをGEOエントリーGSE72855(Neri,F.,Rapelli,S.,Krepelova,A.,Incarnato,D.,Parlato,C.,Basile,G.,Maldotti,M.,Anselmi,F.and Oliviero,S.(2017)Intragenic DNA methylation prevents spurious transcription initiation.Nature,543,72-77)から得て、そのデータを用いて、遺伝子を、その発現レベルに従って4つの群に分けた。
結果
eeTAPSの開発
TAPS反応後のシーケンシングのために、メチル化CpG部位を濃縮するために、DHUを含むTAPS産物を特異的に切断するエンドヌクレアーゼを特定した。DHUまたは構造の似ているヌクレオチド(ウラシル、5-ヒドロキシメチルウラシル、ジヒドロチミン)を消化する能力が知られている10個の市販のエンドヌクレアーゼを試験した。USER、エンドヌクレアーゼVIII、エンドヌクレアーゼIII、及びFpgを含むヌクレアーゼは、TAPSで変換されたDNAを切断したが、APE1及びUDGのような他のヌクレアーゼは、TAPSで変換されたDNAを実質的に切断しなかった(図28A)。USERは、TAPSで変換されたDNAの切断効率が最も高かったとともに、変換されなかったDNAに対する影響が最も小さかったので(図28A)、USERを選択した。
続いて、TAPSによる変換をUSER消化と組み合わせて、メチル化された配列を濃縮した。まず、マウス胚性幹細胞(mESC)に由来する、断片化されていないゲノムDNA(gDNA)をTAPSで変換し、USERで消化した。切断によって、100bp~10kbの範囲のDNA断片が得られた(図28B)。どうやら、短い断片ほど、密集してメチル化された領域に対応し、長い断片は、ゲノムのうち、まばらにメチル化された部分に対応するようである。断片化したDNAのサイズ選別を行って、メチル化状態が中程度となるように、200bp~1kbの断片を保持し、Illuminaシーケンシングライブラリーを調製した(図29A、図28B)。メチル化CpG部位を特定及び定量するために、計算パイプラインを作成した。各CpG部位で切断されるリードの数を、各CpG部位で切断されるかまたは各CpG部位をカバーするリードの総数で除した数として、メチル化レベルを計算した(図29A)。
eeTAPSの性能を評価するために、HpaIIメチルトランスフェラーゼによってメチル化したCCGG部位におけるすべてのCpGで、4kbのスパイクインモデルDNAを調製し、このメチルトランスフェラーゼによっては、オフターゲットの非CCGG部位においても、いくらか低いレベルでのCpGのメチル化が行われた。モデルDNAにおいては、eeTAPSによるメチル化とバイサルファイトによるメチル化の間での一致性が優れていたことから(ピアソンの相関係数(r)=0.98)(図28C)、DNAメチル化レベルの定量における、eeTAPSの検出力が裏付けられた(図29B)。一方で、USER酵素を用いて、TAPSで変換されなかった4kbのモデルDNAを消化したコントロール試料では、有意なメチル化を有するCpGは検出されなかったことから(図29B)、メチル化の検出における、eeTAPSの高い特異性が示された。勘案すると、これらの結果から、eeTAPSが、4kbのモデルDNAにおいて、DNAのメチル化状態を正確に導き出せることが示された。
mESCでのeeTAPS
モデルDNAに対するeeTAPSの能力を示したので、eeTAPSを用いて、mESC gDNAにおけるCpGのメチル化をプロファイリングした(表11)。eeTAPSは、メチル化CpGを濃縮することになるので、コスト効率の良い方法であることが提議されている。実際、eeTAPSにおける断片の84.6%がC/Gとなったことを我々は見出した(図31A)。切断部位と、最も近いCpGとの距離に対するさらなる解析により、切断イベントの72.7%が、CpGで生じたことが特定された(図31B)。
この点をさらに示すために、eeTAPSをwgTAPS及びrrTAPSと比較した(図30A、表11)。まず、3つのすべての方法でカバーされるCpG部位の数を比較した(カバレッジが3超のCpGを、カバーされるCpGとみなした)。wgTAPS及びeeTAPSでは、CpG部位の大半がカバーされた(wgTAPSでは、19.3M個の部位、全CpGの92.1%、eeTAPSでは15.2M個の部位、全CpGの74.2%)のに対して、rrTAPSでカバーされたのは、約1.5M個の部位(全CpGの7.2%)に過ぎなかった(図30B)。これらのアッセイによってカバーされるゲノム領域をさらに比較するために、カバーされた部位を、様々なゲノム領域(Bogu,G.K.,Vizan,P.,Stanton,L.W.,Beato,M.,Di Croce,L.and Marti-Renom,M.A.(2015)Chromatin and RNA Maps Reveal Regulatory Long Noncoding RNAs in Mouse.Mol Cell Biol,36,809-819)にマッピングした。遠位の調節エレメントにおけるメチル化のような、遺伝子間のメチル化も、遺伝子調節に重要であると知られている。wgTAPS及びeeTAPSは、類似の広範な特徴分布を共有し、カバーされたCpG部位の大半が、ヘテロクロマチンに位置する(wgTAPSでは65.6%、eeTAPSでは71%)のに対して、rrTAPSは、プロモーター領域に偏る(カバーされたCpGの46%)ことを我々は見出した(図30B)。1CpGの分解能では、eeTAPS及びwgTAPSにおいて、良好な相関性が見られ(r=0.56、図30C)、それに対して、rrTAPS及びwgTAPSにおいては、非常に高い相関性が見られた(r=0.92、図31C)。eeTAPSでは、メチル化が過剰推定され、これは、非メチル化CpGを有するDNA断片の方が長くなり、あまりうまく増幅されないことによる可能性が高い。
次に、異なるアッセイでカバーされるメチル化CpG部位を比較した。eeTAPS及びwgTAPSでは、メチル化CpG(mCpG)部位として定義される部位に関して、高い一致性が見られる(wgTAPSによって検出された82.4%のmCpG部位が、eeTAPSによっても検出される。図30D)が、rrTAPSは、約20.0%のmCpGを検出するのみである(図30D)。さらに、eeTAPSでは、高い再現性が見られ、レプリケートにおいて、81%のmCpGが観察された(図31D)。勘案すると、これらの解析により、eeTAPSが、全ゲノム規模で、mCpG部位を正確かつ確実に検出でき、1CpGの分解能で、メチル化を測定するための強力な半定量的ツールであり得ることが裏付けられている。
ゲノム機構でのeeTAPS及びwgTAPSの比較
eeTAPSとwgTAPSの間で、異なるゲノム機構にわたるメチル化パターンを比較した。領域内のメチル化レベルを定量するために、wgTAPSでは、メチル化の平均値を使用し、eeTAPSでは、検出されたCpG部位の総数に対する、メチル化CpGの割合を使用した。eeTAPS及びwgTAPSにおいて、相関性の高い、染色体ワイドなメチル化パターンが見られたが、eeTAPSでは、メチル化レベルが過剰に推定された(図32A、B)。CpGアイランド(CGI)は、DNAのメチル化が非常に少ないことが知られており、これらは、eeTAPS及びwgTAPSの両方に反映されている(図32C)。eeTAPS及びwgTAPSを用いて、CGIでのメチル化レベルを測定したものの相関性は0.81だったことから、eeTAPSは、様々な機構において、CpGのメチル化状態を正確に捉えることができることがさらに示されている(図32D)。
以前の研究により、プロモーター領域におけるDNAのメチル化は概ね、遺伝子の発現と逆相関することが明らかにされている。我々は、遺伝子を、その発現レベルに従って4つの群に分け、転写開始部位(TSS)の4kb上流から4kb下流までのメチル化の平均をプロットした。我々は、eeTAPS及びwgTAPSの両方を用いて、高度に発現する遺伝子が、メチル化レベルの低い傾向があり、発現レベルの低い遺伝子が、メチル化レベルの高いことを見出した(図33)。我々は、以前に定義されたようにして(Bogu,G.K.,Vizan,P.,Stanton,L.W.,Beato,M.,Di Croce,L.and Marti-Renom,M.A.(2015)Chromatin and RNA Maps Reveal Regulatory Long Noncoding RNAs in Mouse.Mol Cell Biol,36,809-819)、クロマチンの異なる特徴においても、メチル化分布を比較した(図32E、F)。以前の研究と一致して、ヘテロクロマチン領域は、高度にメチル化される一方で、ユークロマチンにおけるプロモーター領域は大抵の場合、CGのメチル化が少ない(図32F)。
インプットの少ない試料に対するeeTAPSの適用
インプットの少ない試料でのeeTAPSの性能を評価するために、我々は、eeTAPSに、1ng、10ng及び50ngのmESC gDNAをそれぞれ適用した。200ngのmESC DNA試料では、シーケンシングリードを2分の1にダウンサンプリングして、インプットの少ない試料のシーケンシング深度と適合させた。eeTAPSにおいて、1ngのDNAを用いても、wgTAPSによって特定されたmCpG部位の27%が再現されることを我々は見出した。50ngのmESC DNAを使用したところ、そのパーセンテージは、47%まで上昇した(図34A)。これらのインプットの少ない試料によって、全ゲノムのメチル化プロファイルをさらに比較するために、我々は、ゲノムを100kbのウィンドウにビニングし、各ビンにおいて、メチル化レベルの平均を計算した(図34B)。これらのインプットの少ない試料の間で、一貫性の高いメチル化プロファイルが観察されたことから(2分の1の200ngでのeeTAPSと比較した1ngでのr=0.88、20ngでのr=0.92、50ngでのr=0.95。図34C)、インプットの少ないDNA試料に対するeeTAPSの適用が実行可能であることがさらに示された。
シーケンシング深度がeeTAPSに及ぼす作用
シーケンシング深度が、検出できるmCpG部位の総数に及ぼす作用を評価するために、我々は、eeTAPSをダウンサンプリングし、その性能を評価した。検出されるmCpG部位の総数は、シーケンシング深度が深いほど増加した(図35A)。とはいえ、4×(70M個のリード)のシーケンシング深度では、74%のmCpG部位を検出するのに成功できた(4×eeTAPSでも、wgTAPSで検出された14.9M個のmCpG部位のうち、10.9M個の部位がmCpGとして定義された)。染色体及びCGIの全体にわたるメチル化の相関性の点でも、同様の傾向が観察され(図35B)、CGIにおけるピアソンの相関係数は、4×カバレッジでは、0.83に達した(図35B)。したがって、eeTAPSが、WGBSと比べて、低いシーケンシングコストで、グローバルなメチル化プロファイルを正確にもたらすことができることを我々は実証した。
考察
wgTAPSは、網羅性が最も高い定量的かつ塩基単位の分解能の全ゲノムメチル化解析を行うことができた。しかしながら、全ゲノムシーケンシングの高いコストと、大量の処理データにより、多くのプロジェクトにおける、メチル化解析の広範な用途は依然として制限を受ける。メチル化CpG部位は、哺乳動物のゲノムにおけるほんの一部であるので、メチル化の状態を調べるには、全ゲノムシーケンシングは、データ効率が最も高いアプローチというわけではない。コスト効率の良いアプローチは、シーケンシングによってさらに解析する際に、メチル化CpGを含む領域のみを特異的に選択するものであろう。CpGリッチな領域の制限酵素消化濃縮に基づくReduced-Representation Sequencingを行い、その後にバイサルファイトシーケンシングを行うことは、費用効果の高いメチローム解析アプローチであるが、この方法では、ゲノムにおけるCpG部位のごく一部がカバーされたに過ぎなかった(Meissner,A.,Gnirke,A.,Bell,G.W.,Ramsahoye,B.,Lander,E.S.and Jaenisch,R.(2005)Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis.Nucleic Acids Res,33,5868-5877)。TAPSは、Reduced-Representationのアプローチと適合しており、rrTAPSが、ゲノムのサブセットにおいて、特にCGIにおいて、メチル化を正確に定量できることを我々は実証した。遺伝子プロモーターにおけるCpGのメチル化の、十分に確立された生物学的意義に加えて、遺伝子間のDNAメチル化は、細胞運命決定及び腫瘍形成に関与する可能性から、広範な研究で注目されている。濃縮アプローチをゲノムワイドなCpG部位まで広げるために、我々はTAPSの利点をさらに用いて、5mCをDHUに直接変換し、それにより、DHU感受性のエンドヌクレアーゼによる誘導によって、これらの特異的な改変塩基での切断を可能にした。これらの断片の選択的な濃縮をシーケンシングと組みわせることを通じて、eeTAPSにより、CpGのメチル化をゲノムワイドな規模で検出可能になることを我々は実証した。TAPSによってメチル化シトシンが直接検出されるため、このような方策が可能である。従来の抗体ベースの濃縮方法とは異なり、eeTAPSは、1CpGの分解能での直接的なメチル化検出の可能性をもたらす。
eeTAPSを用いて、費用効果の良い形式で、1CpGの分解能で、ゲノムワイドなメチル化の特徴を捉えることができ、それにより、rrTAPS及びwgTAPS間のギャップが埋まることを我々は実証した。複数の異なるゲノム機構にわたる、eeTAPSによるメチル化プロファイルは、wgTAPSを用いて得られるプロファイルとの相関性が高かった。さらに、eeTAPSは、わずか70M個のリードで、wgTAPSによって検出されたメチル化CpGの74%を検出できる。eeTAPSの考え得る制約(メチル化CpG部位の間隔のばらつきに起因する)は、単一のCpG部位でのメチル化レベルの半定量的な測定である可能性がある。この研究では、我々は、概念実証として、200bp~1kbの断片を選択した。しかしながら、それでも、wgTAPSとeeTAPSとの相関係数は0.56ほどの良好さである。さらに、TAPS反応の穏和な性質を足場として、インプットの少ないDNA試料でのメチル化の検出において、eeTAPSが、費用効果が高い有望なプロトコールであることも我々は示した。
ボラン含有化合物のスクリーニングである。11merのオリゴヌクレオチド(「オリゴ」)における5caCのDHUへの変換について、ボラン含有化合物をスクリーニングし、変換率は、MALDIによって推定した。2-ピコリンボラン(pic-ボラン)、ボランピリジン及びtert-ブチルアミンボランは、5caCをDHUに完全に変換できたが、エチレンジアミンボラン及びジメチルアミンボランでは、変換率は30%前後であった。ジクロロヘキシルアミンボラン、モルホリンボラン、4-メチルモルホリンボラン及びトリメチルアミンボランでは、検出可能な生成物は、測定されなかった(不検出)。水素化ホウ素ナトリウム及びトリ(アセトキシ)水素化ホウ素ナトリウムのような他の還元剤は、酸性培地中で速やかに分解されて、変換が不完全となった。シアノ水素化ホウ素ナトリウムは、酸性条件下では、シアン化水素が形成される可能性があることから、使用しなかった。完全な変換、低い毒性及び高い安定性から、Pic-ボラン及びピリジンボランを選択した。 DNAオリゴにおけるPic-ボラン反応である。(A)5caCを含む11merのモデルDNAをpic-ボランで処理したものをMALDIによって特徴付けたものである。各グラフの上に、質量(m/z)の計算値が示されており、ピークの左に、質量の観察値が示されている。(B)dC及び各種シトシン誘導体の変換率をHPLC-MS/MSによって定量した。データは、3つのレプリケートの平均±標準偏差としてされている。 単一ヌクレオシドのpic-ボラン反応である。H NMRの結果は、2’-デオキシ-5,6-ジヒドロウリジンに関する過去の報告に従ったものである(I.Aparici-Espert et al.,J.Org.Chem.81,4031-4038(2016))。単一ヌクレオシドのpic-ボラン反応生成物のH NMR(MeOH-d、400MHz)チャートである。δ ppm: 6.28 (t, 1H, J = 7 Hz), 4.30 (m, 1H), 3.81 (m, 1H), 3.63 (m, 2H), 3.46 (m, 2H), 2.65 (t, 2H, J = 6 Hz), 2.20 (m, 1H), 2.03 (m, 1H). 単一ヌクレオシドのpic-ボラン反応である。13C NMRの結果は、2’-デオキシ-5,6-ジヒドロウリジンに関する過去の報告に従ったものである(I.Aparici-Espert et al.,J.Org.Chem.81,4031-4038(2016))。単一ヌクレオシドのpic-ボラン反応生成物の13C NMR (MeOH-d, 400 MHz)チャートである。δ ppm: 171.56 (CO), 153.54 (CO), 85.97 (CH), 83.86 (CH), 70.99 (CH), 61.92 (CH), 36.04 (CH), 35.46 (CH), 30.49 (CH). 5caCのDHUへの、ボランによる変換と、5caCのDHUへのボラン反応の機序案を示す図である。 5fCのDHUへの、ボランによる変換と、5fCのDHUへのボラン反応の機序案を示す図である。 (A)TAPS法が、5mC及び5hmCの両方をDHUに変換し、そのDHUが、複製時にはチミンとして機能することを示す図である。(B)TAPS、TAPSβ及びCAPSの方法の概説である。 5fC及び5caCを含むモデルDNAオリゴをpic-ボランによって、5fC及び5caCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5fC及び5caCは、pic-BHによってジヒドロウラシル(DHU)に変換される。5fCは、O-エチルヒドロキシルアミン(EtONH)のようなヒドロキシルアミン誘導体によってブロックしたが、その誘導体は、オキシムとなり、pic-ボランによる変換を阻止することになる。5caCは、EDC結合を介して、エチルアミンによってブロックして、pic-ボランによる変換をブロックするアミドに変換した。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。 5fC及び5caCを含むモデルDNAオリゴをpic-ボランによって、5fC及び5caCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5fC及び5caCは、pic-BHによってジヒドロウラシル(DHU)に変換される。5fCは、O-エチルヒドロキシルアミン(EtONH)のようなヒドロキシルアミン誘導体によってブロックしたが、その誘導体は、オキシムとなり、pic-ボランによる変換を阻止することになる。5caCは、EDC結合を介して、エチルアミンによってブロックして、pic-ボランによる変換をブロックするアミドに変換した。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。 5mC及び5hmCを含むモデルDNAオリゴをKRuO及びpic-ボランによって、5hmCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5hmCは、βGTによってグルコースでブロックして、5gmCに変換できた。5mC、5hmC及び5gmCは、pic-ボランによって変換できなかった。5hmCは、KRuOによって酸化させて、5fCにしてから、pic-ボランによってDHUに変換できた。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。 5mC及び5hmCを含むモデルDNAオリゴをKRuO及びpic-ボランによって、5hmCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5hmCは、βGTによってグルコースでブロックして、5gmCに変換できた。5mC、5hmC及び5gmCは、pic-ボランによって変換できなかった。5hmCは、KRuOによって酸化させて、5fCにしてから、pic-ボランによってDHUに変換できた。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。 制限酵素消化によって、TAPSが、5mCをTに効果的に変換できることが示された。(A)TAPSによる配列変化を確認するための制限酵素消化アッセイの図である。(B)TAPSによる、CからTへの変換を確認するためのTaqαI消化試験である。TAPSは、TaqαI制限部位を有するとともに、完全にメチル化された5つのCpG部位(5mC)を含む222bpのモデルDNAと、その非メチル化コントロール(C)で行った。5mC、C及びC TAPSに示されているように、PCR増幅した222bpのモデルDNAは、TaqαIによって、約160bp及び約60bpの断片に切断できる。5mC-TAPSのレーンに示されているように、メチル化DNAでTAPSを行った後、T(mC)GA配列は、TTGAに変換され、TaqαI消化によっては切断されなくなる。 222bpのモデルDNA及びmESC gDNAでのTAPSである。(A)TAPSの前(5mC、C)及びTAPSの後(5mC TAPS、C TAPS)における、完全にメチル化された5つのCpG部位を含む222bpのモデルDNAと、その非メチル化コントロールのサンガーシーケンシングの結果である。5mCのみが、TAPS法によってTに変換される。(B)NgTET1による酸化後及びpic-ボランによる還元後のmESC gDNAコントロールにおける相対的な改変レベルの、HPLC-MS/MSによる定量結果である。データは、3つのレプリケートの平均±標準偏差としてされている。 TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。氷浴で冷却する前(の222bpの非メチル化DNA、222bpのメチル化DNA及びmESC gDNAのアガロースゲル画像である。TAPS後、検出可能なDNA分解は、観察されず、DNAは、二本鎖を保ち、冷却せずに見ることができた。バイサルファイト変換では、分解が起こり、DNAは、一本鎖となり、氷上で冷やした後のみに見ることができた。 TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。氷浴で冷却した後(B)の222bpの非メチル化DNA、222bpのメチル化DNA及びmESC gDNAのアガロースゲル画像である。TAPS後、検出可能なDNA分解は、観察されず、DNAは、二本鎖を保ち、冷却せずに見ることができた。バイサルファイト変換では、分解が起こり、DNAは、一本鎖となり、氷上で冷やした後のみに見ることができた。 TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。TAPS及びバイサルファイトで処理した、様々な断片長のmESC gDNAであって、氷上で冷やす前(左パネル)及び氷上で冷やした後(右パネル)のmESC gDNAのアガロースゲル画像である。TAPS後のDNAは、二本鎖を保ち、ゲル上で、直接見ることができた。バイサルファイト処理の方が、試料の損傷及び断片化が大きく、DNAは、一本鎖となり、氷上で冷やした後のみに見ることができた。図15に示されているように、TAPSによる変換は、すべてのgDNAにおいて、断片の長さにかかわらず、完全なものであった。 TAPSでは、バイサルファイトと比べて、有意なDNA分解が起こらない。TAPSの前及び後の222bpのモデルDNAのアガロースゲル画像であり(独立して3回反復)、その反応後、検出可能な分解は見られなかった。 TAPSの前及び後のモデルDNA間での増幅曲線及び融解曲線の比較である。qPCRアッセイによって、TAPSの前及び後のモデルDNAで、増幅曲線において小さな違いが見られた。メチル化DNA(5mC)の融解曲線は、TAPSの後、低い温度の方にシフトしたことから、Tmが低下する、CからTへの変換の可能性が示されたが、非メチル化DNA(C)では、シフトは見られなかった。 サンガーシーケンシングによって示されたように、TAPS、TAPSβ及びCAPSの後、CからTへの完全な変換が誘導された。本明細書に記載されているようにして、1つのメチル化CpG部位及び1つのヒドロキシメチル化CpG部位を含むモデルDNAを調製した。本明細書に記載されているように、NgTET1による酸化及びピリジンボランによる還元のプロトコールに従って、TAPSによる変換を行った。5hmCのブロック、NgTET1による酸化及びピリジンボランによる還元のプロトコールに従って、TAPSβによる変換を行った。5hmCの酸化及びピリジンボランによる還元のプロトコールに従って、CAPSによる変換を行った。変換後、変換されたDNA試料1ngをTaq DNAポリメラーゼによってPCR増幅し、サンガーシーケンシング用に処理した。TAPSは、5mC及び5hmCの両方をTに変換した。TAPSβは、5mCを選択的に変換し、それに対し、CAPSは、5hmCを選択的に変換した。これらの3つの方法のいずれでも、無修飾シトシン及びその他の塩基は変換されなかった。 TAPSは、様々なDNAポリメラーゼ及びRNAポリメラーゼと適合しており、サンガーシーケンシングによって示されているように、CからTへの完全な変換を誘導する。ポリメラーゼ試験用の、メチル化CpG部位を含むモデルDNAとプライマー配列は、本明細書に記載されている。TAPS処理後、5mCは、DHUに変換された。KAPA HiFi Uracil Plusポリメラーゼ、Taqポリメラーゼ及びVentエキソポリメラーゼは、DHUをTとして読み取るので、PCR後、CからTへの完全な変換を誘導する。あるいは、Klenow断片、Bst DNAポリメラーゼ及びphi29 DNAポリメラーゼを含め、ビオチン標識プライマー及び等温ポリメラーゼで、プライマー伸長を行った。新たに合成したDNA鎖をDynabeads MyOne Streptavidin C1によって分離してから、Taqポリメラーゼを用いたPCRによって増幅し、サンガーシーケンシング用に処理した。T7 RNAポリメラーゼは、DHUを効率的にバイパスして、DHU部位の反対側にアデニンを挿入できたが、これは、RT-PCR及びサンガーシーケンシングによって証明されている。 特定の他の市販されているポリメラーゼは、DHUを含むDNAを効率的には増幅しなかった。TAPS処理後、5mCは、DHUに変換された。KAPA HiFi Uracil Plusポリメラーゼ及びTaqポリメラーゼは、DHUをTとして読み取るので、CからTへの完全な変換を誘導することになる。KAPA HiFiポリメラーゼ、Pfuポリメラーゼ、Phusionポリメラーゼ及びNEB Q5ポリメラーゼを含め、特定の他の市販されているポリメラーゼでは、CからTへの変換が低度であったか、または変換されなかったことが観察された(不図示)。 DHUでは、T及びCと比べて、PCRバイアスが見られない。本明細書に記載されているように、対応するDNAオリゴを用いて、1つのDHU/U/T/C修飾を含むモデルDNAを合成した。qPCR反応に基づき、DHU/U/T/C修飾を有する各モデルDNAの標準曲線をプロットし、モデルDNAインプットは、1:10の段階希釈を行った(0.1pg~1ng。すべてのqPCR実験をトリプリケートで行った)。ExcelのSLOPE関数によって、対数濃度(ng)の値とCt平均値との回帰の傾きを計算した。効率(%)=(10^(-1/傾き)-1)×100%という式を用いて、PCR効率を計算した。増幅係数=10^(-1/傾き)という式を用いて、増幅係数を計算した。DHU、TまたはCの修飾を有するモデルDNAでのPCR効率は、ほぼ同じであったことから、DHUは、通常の塩基として読み取ることができ、PCRバイアスの原因とならないと見られることが示された。 TAPSは、DNA断片の長さにかかわらず、5mCをTに完全に変換した。(A)TaqαI-消化アッセイのアガロースゲル画像により、すべての試料において、DNA断片の長さにかかわらず、5mCからTへの完全な変換が確認された。TAPS後、λゲノムに由来する194bpのモデル配列をPCR増幅し、TaqαI酵素で消化した。変換されなかった試料から増幅されたPCR産物は、切断できたが、TAPS処理試料で増幅された産物は、インタクトなままであったことから、制限部位の喪失、すなわち、5mCからTへの完全な変換が示唆されている。(B)ゲルバンドの定量によって、CからTへの変換パーセンテージを推定したところ、試験したすべてのDNA断片長で、100%となった。 異なるTAPS条件での変換及び偽陽性である。mTet1及びピリジンボランを組み合わせたところ、NgTET1またはpic-ボランを用いた他の条件と比べて、メチル化Cの変換率が最も高くなり(96.5%。完全にCpGがメチル化されたλDNAで計算した)と、無修飾のCの変換率が最も低くなった(0.23%、2kbの無修飾スパイクインで計算した)。バーの上に示されているのは、試験したすべてのシトシン部位の変換率±標準誤差である。 短いスパイクインにおける変換率である。5mC及び5hmCを含む(A)120mer-1及び(B)120mer-2である。両方の鎖の5mC及び5hmCの部位で、ほぼ完全な変換が行われた。実際の配列は、修飾の状態とともに、上下に示されている。 全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。TAPSで処理したDNAにおける5mC及び5hmCの変換率である。左:既知の位置においてメチル化またはヒドロキシメチル化された合成スパイクイン(CpN)である。 全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。無修飾の2kbのスパイクインから得られた、TAPSの偽陽性率である。 全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。100万個のシミュレーションリードを1つのIntel Xeon CPUの1つのコアで処理した時の、TAPS及びWGBSの総実行時間である。 全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。シーケーシングしたすべてのリードペア(トリミング後)のうち、ゲノムにマッピングされた割合である。 全ゲノムバイサルファイトシーケンシング(WGBS)よりも改善されたTAPSのシーケンシングクオリティ。シーケンシングしたすべてのリードペアにおける第1及び第2のリードにおける1塩基あたりのシーケンシングクオリティスコアをIllumina BaseSpaceによって記録したものである。上は、TAPSである。下は、WGBSである。 TAPSでは、WGBSよりも均一なカバレッジが得られ、カバーされなかった位置がWGBSよりも少なかった。WGBSとTAPSの間で、(A)すべての塩基及び(B)CpG部位におけるカバレッジ深度を比較したものであり、両方の鎖で計算した。「TAPS(ダウンサンプリング済み)」では、マッピングされたすべてのTAPSリードから取ったランダムなリードを選択して、カバレッジ中央値が、WGBSのカバレッジ中央値と一致するようにした。最後のビンには、カバレッジが50×超の位置が示されている。 すべての染色体にわたる修飾レベルの分布である。マウス染色体全体にわたる、100kbのウィンドウにおける改変レベルの平均を、CpGのカバレッジによって重み付けして、Gaussian加重移動平均フィルターを用いて、ウィンドウサイズ10で平滑化したものである。 TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。すべてのマウスCpGアイランド(20個のウィンドウにビニングした)及び4kbpのフランキング領域(サイズの等しい50個のウィンドウにビニングした)におけるシーケンシングカバレッジ深度の平均である。シーケンシング深度の差を考慮するために、マッピングされたすべてのTAPSリードをダウンサンプリングして、ゲノム全体にわたる、マッピングされたWGBSリードの数の中央値と適合させた。 TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。TAPSのみ、TAPS及びWGBSの両方、またはWGBSのみによって、少なくとも3本のリードによってカバーされたCpG部位である。 TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。少なくとも3本のリードによってカバーされたとともに、TAPSのみ、TAPS及びWGBS、またはWGBSのみによって検出された修飾レベルが0.1超であるCpG部位の数である。 TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。TAPS及びWGBSにおける修飾レベル(%)の染色体分布の例である。マウス染色体4における、100kbのウィンドウ当たりの修飾CpGの割合の平均を、Gaussian加重移動平均フィルターを用いて、ウィンドウサイズ10で平滑化したものである。 TAPS及びWGBSによるゲノムワイドなメチローム測定値の比較である。TAPS及びWGBSの両方において、少なくとも3本のリードによってカバーされたCpG部位の数を表すヒートマップであり、各方法によって測定した場合の修飾レベルによって分類したものである。コントラストを向上させるために、両方の方法で無修飾のCpGを含む第1のビンは、カラースケールから除外し、星印で示されている。 CpGアイランド周辺の修飾レベルである。CpGアイランド(20個のウィンドウにビニングした)及び4kbpのフランキング領域(サイズの等しい50個のウィンドウにビニングした)における修飾レベルの平均である。カバレッジが3リード未満のビンは無視した。 TAPSでは、WGBSよりもカバレッジ-修飾バイアスが小さい。すべてのCpG部位を、そのカバレッジに従ってビニングし、各ビンにおいて、修飾値の平均(●)及び中央値(黒三角)が、WGBS(A)及びTAPS(B)に関して示されている。最後のビンには、100個超のリードによってカバーされたCpG部位が示されている。線は、データポイントにわたる直線近似を表している。 dsDNA及びssDNAライブラリー調製キットを用いて、インプットの少ないgDNA及びセルフリーDNA TAPSを調製した。(A)dsDNAライブラリーキットNEBNext Ultra IIまたはKAPA HyperPrepキットを用いて、わずか1ngのマウス胚性幹細胞(mESC)ゲノムDNA(gDNA)による、シーケンシングライブラリーの構築が成功した。ssDNAライブラリーキットAccel-NGS Methyl-Seqキットを用いて、インプットDNA量をさらに、(B)0.01ngのmESC gDNAまたは(C)1ngのセルフリーDNAまで減少させた。 dsDNA KAPA HyperPrepライブラリー調製キットを用いて、インプットの少ないgDNA及びセルフリーDNAのTAPSライブラリーを調製した。KAPA HyperPrepキットを用いて、1ngほどの少ない(A)mESC gDNA及び(B)セルフリーDNAで、シーケンシングライブラリーの構築が成功した。セルフリーDNAは、血漿ヌクレアーゼ消化により、160bp前後(ヌクレオソームサイズ)のシャープな長さ分布を有する。ライブラリーの構築後、そのDNAは、約300bpとなり、それは、Bにおけるシャープなバンドである。 高品質なセルフリーDNA TAPSである。(A)TAPSで処理したcfDNAにおける5mCの変換率である。(B)TAPSで処理したcfDNAにおける偽陽性率である。(C)シーケンシングされたすべてのリードペアのうち、ゲノムにユニークにマッピングされた割合である。(D)シーケンシングされたすべてのリードペアのうち、ゲノムにユニークにマッピングされ、PCR重複リードを除去した後の割合である。CHG及びCHHは、非CpGコンテクストである。 TAPSは、遺伝的バリアントを検出できる。メチル化(MOD、上の行)及びC→T SNP(下の行)では、元の上鎖(OT)/元の下鎖(OB)(左のカラム)、ならびにOTと相補的な鎖(CTOT)及びOBと相補的な鎖(CTOB)(右のカラム)において、別々の塩基分布が見られた。 TAPS変換生成物のエンドヌクレアーゼ切断である。(A)TAPS変換の前及び後のmESC gDNAの、異なるエンドヌクレアーゼによる消化の結果である。TAPS変換では、メチル化シトシンに代わって、DHUが導入される。Endo VIIIはエンドヌクレアーゼVIII、Endo IVはエンドヌクレアーゼIV、Tma Endo IIIはエンドヌクレアーゼIII、Tth Endo IVはエンドヌクレアーゼIV、Nth Endo IIIはNthエンドヌクレアーゼIII、Endo VはエンドヌクレアーゼVである。(B)TAPSで処理したmESC gDNAであって、USER消化前、USER消化後、及びサイズ選別後のmESC gDNAの代表的な画像である。(C)BS-seq及びeeTAPSの両方によって測定したすべてのCpGにおけるメチル化レベルを示す散布図である。 エンドヌクレアーゼ濃縮TAPS(eeTAPS)である。eeTAPSの概略図(上図)と、CGメチル化レベルの、コンピューターによる測定値(下図)である。まず、TAPSで5-メチルシトシン(mC)をジヒドロウラシル(DHU)に変換してから、USER消化を通じて濃縮した。その後、サイズ選別を行ったDNA断片を調製して、シーケンシングライブラリーにし、PCRによって増幅した。リードのアラインメント後、各CpG部位で切断されるリードの数を、各CpG部位で切断されたリードの総数、または各CpG部位をカバーするリードの総数で除したものとして、CGメチル化レベルを計算した。 4kbのモデルDNAでのeeTAPSの検証である。トラックは、上から下の順で、バイサルファイトシーケンシング(BS-seq)、eeTAPS、及びeeTAPSコントロールで測定したメチル化レベルを示している。eeTAPSコントロールでは、TAPSで変換せずに、USER酵素を用いて、DNAを消化した。 mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS、eeTAPS及びrrTAPSのゲノムの断片化方法を示す図である。wgTAPSでは、ゲノムDNAは、ランダムに断片化されるのに対して、断片化は、eeTAPSでは、メチル化CpG(mCG)部位で、rrTAPSでは、CCGG部位で特異的に行われる。 mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS、eeTAPS及びrrTAPSによってカバーされたCpG部位のパーセンテージを、クロマチンの異なる特徴と重ね合わせたものを示すバープロットである。これらのクロマチンの特徴は、以前の研究(20)で定義されている。 mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS及びeeTAPSによって、1CpGの分解能で求めたメチル化レベルの相関性を示すヒートマップである。wgTAPS及びeeTAPSの両方において、メチル化レベルを16群に分けた。カラーは、所定の間隔でのCpGの数を示している。wgTAPSでのカバレッジが5超の部位のみを考慮に入れた。プロットの上部には、ピアソンの相関係数が示されている。 mESC DNAでのwgTAPS、eeTAPS及びrrTAPSの比較である。wgTAPS、rrTAPS及びeeTAPSにおける、検出されたmCpG部位の重複を示すベンプロットである。wgTAPS及びeeTAPSの両方において、メチル化レベルが、メチル化レベルの第1四分位超のCpG部位を、メチル化CpGとして定義した。rrTAPSでは、wgTAPSと同じメチル化カットオフ値を用いた。(wgTAPS及びrrTAPSでは、メチル化レベルが0.5超のCpG部位を、mCpGとして定義し、eeTAPSでは、メチル化レベルが0.28超のCpG部位をmCpGとして定義した。) (A)シーケンシングした断片の末端におけるヌクレオチド頻度である。(B)切断部位と、その部位に最も近いCpGとの距離の分布を示すバープロットである。(C)wgTAPS及びrrTAPSによって、1CpGの分解能で求めたメチル化レベルの相関性を示すヒートマップである。wgTAPS及びrrTAPSの両方において、メチル化レベルを16群に分けた。カラーは、所定の間隔でのCpGの数を示している。wgTAPSでのカバレッジが5超であり、rrTAPSでのカバレッジが5超である部位のみを考慮に入れた。ピアソンの相関係数は、0.92である。(D)eeTAPSのレプリケートで検出されたmCpG部位の重複であり、レプリケートは、同じ深度になるようにサブサンプリングして、mCpGを検出した。 異なるゲノム機構におけるメチル化をeeTAPSでプロファイリングしたものである。(A)mESCの1番染色体において、eeTAPS(上の線)及びwgTAPS(下の線)の両方によって測定したメチル化レベルである。100kbのウィンドウを使用し、移動平均値をRにおいて、bw=10で、movAvg2関数によって計算した。(B)eeTAPS及びwgTAPS間での、染色体ビンにおけるメチル化レベルの相関性を示す密度プロットである。100kbのウィンドウを用いて、各ビンにおいて、メチル化レベルの平均を計算した。プロットの上部には、ピアソンの相関係数が示されている。(C)eeTAPS及びwgTAPSにおける、CpGアイランド(CGI)及び4kbのフランキング領域にわたるメチル化レベルの平均である。(D)eeTAPS及びwgTAPS間での、CpGアイランドにおけるメチル化レベルの相関性を示す密度プロットである。プロットの上部には、ピアソンの相関係数が示されている。(E)eeTAPS及びwgTAPSでの、クロマチンの異なる特徴におけるメチル化レベルの相関性を示す密度プロットである。これらのクロマチンの特徴は、ヒストンマーカーによって以前に定義されたものであり(Bogu,G.K.,Vizan,P.,Stanton,L.W.,Beato,M.,Di Croce,L.and Marti-Renom,M.A.(2015)Chromatin and RNA Maps Reveal Regulatory Long Noncoding RNAs in Mouse. Mol Cell Biol,36,809-819)、(F)にも示されている。プロットの上部には、ピアソンの相関係数を示した。(F)クロマチンのすべての特徴にわたるメチル化レベルの分布をeeTAPS(下のバー)及びwgTAPS(上のバー)によって測定したものを示すボックスプロットである。 wgTAPS及びeeTAPSにおける、転写開始部位(TSS)周辺のメチル化分布の平均である。遺伝子は、GSE72855のデータセットに従い、それらの発現レベルによって分類した。 インプットの少ない試料での、eeTAPSによる解析である。(A)eeTAPSを用いて、1ng、10ng、50ng、200ngのmESC gDNAインプットで検出されたmCpGの数(wgTAPSによって特定されたmCpGのうちの数)である。200ngのmESCでは、リードを2分の1にダウンサンプリングして、インプットの少ない試料に合わせて、シーケンス深度を適合させた。0.28超のメチル化レベルという基準を用いて、mCpGを割り当て、eeTAPSでは、1超の切断数をmCpGとして割り当て、wgTAPSでは、0.5超のメチル化レベルをmCpGとして割り当てた。バーの上に示されているパーセンテージは、検出されたmCpGのパーセンテージである(wgTAPSで検出されたmCpGは、真として定義する)。(B)異なるインプットレベルを用いて、マウスゲノムにわたり、eeTAPSで測定したメチル化分布を示すヒートマップである。各染色体を100kbのウィンドウに分け、それらのウィンドウは、ヒートマップの行によって表されている。メチル化レベルは、100kbの各ウィンドウにおいて、メチル化CpG部位の数を、カバーされたCpG部位の総数で除したものとして定義した。(C)インプットの少ない試料間での、200ngのインプット試料に対するメチル化の相関性を示す密度プロットである。メチル化レベルは、(B)に示されているように、100kbのウィンドウにおいて、全マウスゲノムにわたり計算した。ピアソンの相関係数は、各プロットについて示されている。 eeTAPSシーケンシング深度の解析結果である。(A)1~10×のシーケンシング深度でリードをサンプリングしたときに検出されるメチル化CpGの数である。上に示されているパーセンテージは、eeTAPSによって検出されたmCpGのパーセンテージである(wgTAPSで検出されたmCpGを真と定義する)。(B)1~10×のシーケンシング深度でリードをサンプリングしたときの、100kbのウィンドウにおける、全マウスゲノム(上の線)及びCpGアイランド(CGI)(下の線)におけるメチル化の相関性である。

Claims (16)

  1. 改変標的DNAの切断方法であって、
    標的DNA中の5caC及び/または5fCをジヒドロウラシル(DHU)に変換し、1つ以上のDHUを含む改変標的DNAを提供する変換工程と、
    1つ以上のDHUを含む前記改変標的DNAと、前記1つ以上のDHUの位置または前記1つ以上のDHUと隣接する位置で、前記改変標的DNAを切断する1つ以上のエンドヌクレアーゼとを接触させることとを含む、方法。
  2. 前記方法が、前記変換工程の前に、5mC及び/または5hmCを含む標的DNAを改変することをさらに含み、前記改変することが、前記標的DNAにおける5mC及び/または5hmCを5-カルボキシルシトシン(5caC)及び/または5-ホルミルシトシン(5fC)に変換する工程を含む、請求項1に記載の方法。
  3. 前記方法が、前記5mC及び/または5hmCを変換する工程の前に、5mC及び5hmCを含む前記標的DNAを改変することをさらに含み、前記改変することが、前記標的DNAにおける5hmCに、ブロック基を付加する工程を含む、請求項2に記載の方法。
  4. 前記5hmCに付加する前記ブロック基が、糖である、請求項3に記載の方法。
  5. 前記標的DNAにおける5mC及び/または5hmCを5caC及び/または5fCに変換する工程が、前記標的DNAと、10-11転座(TET)酵素とを接触させることを含む、請求項2から4のいずれか1項に記載の方法。
  6. 前記TET酵素が、ヒトTET1、ヒトTET2、ヒトTET3、マウスTET1、マウスTET2、マウスTET3、Naegleria TET(NgTET)、およびCoprinopsis cinerea TET(CcTET)からなる群から選択される、請求項5に記載の方法。
  7. 前記TET酵素が、NgTET、マウスTET2、またはヒトTET1である、請求項5または6に記載の方法。
  8. 前記5hmCを5caC及び/または5fCに変換する工程が、前記標的DNAを酸化剤と接触させることを含む、請求項2から4のいずれか1項に記載の方法。
  9. 前記酸化剤が、過ルテニウム酸カリウム、Cu(II)/TEMPO、ルテニウム酸カリウムまたは酸化マンガンである、請求項8に記載の方法。
  10. 前記5caC及び/または5fCをDHUに変換する工程が、前記標的DNAを還元剤と接触させることを含む、請求項1から9のいずれか1項に記載の方法。
  11. 前記還元剤が、ピリジンボラン、2-ピコリンボラン(pic-BH)、ボラン、水素化ホウ素ナトリウム、シアノ水素化ホウ素ナトリウム及びトリアセトキシ水素化ホウ素ナトリウムからなる群から選択される、請求項10に記載の方法。
  12. 前記1つ以上のエンドヌクレアーゼが、アプリン酸/アピリミジン酸エンドヌクレアーゼ1(APE1)、ウラシルDNAグリコシラーゼ(UDG)、エンドヌクレアーゼIII、TmaエンドヌクレアーゼIII、TthエンドヌクレアーゼIV、エンドヌクレアーゼV、エンドヌクレアーゼVIII、ホルムアミドピリミジンDNAグリコシラーゼ(Fpg)及びhNEIL1からなる群のうちの1つ以上から選択される、請求項1から11のいずれか1項に記載の方法。
  13. 前記1つ以上のエンドヌクレアーゼが、ウラシルDNAグリコシラーゼ(UDG)及びDNAグリコシラーゼ-リアーゼエンドヌクレアーゼVIIIを組み合わせたものを含む、請求項1から12のいずれか1項に記載の方法。
  14. 前記切断された改変標的DNAに、アダプターDNA分子を付加する工程をさらに含む、請求項1から13のいずれか1項に記載の方法。
  15. 前記切断された改変標的DNAの配列を検出する工程をさらに含み、切断された部位の存在は、前記標的DNAにおける5caCまたは5fCの何れかの位置を提供する、請求項1から14のいずれか1項に記載の方法。
  16. 前記方法が、前記接触させる工程の前に、前記改変標的DNAのコピー数を増幅する工程をさらに含む、請求項1から15のいずれか1項に記載の方法。
JP2022500944A 2019-07-08 2020-07-08 バイサルファイトを使用しない全ゲノムメチル化解析 Active JP7597783B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2024207194A JP2025029007A (ja) 2019-07-08 2024-11-28 バイサルファイトを使用しない全ゲノムメチル化解析

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201962871444P 2019-07-08 2019-07-08
US62/871,444 2019-07-08
PCT/IB2020/056435 WO2021005537A1 (en) 2019-07-08 2020-07-08 Bisulfite-free, whole genome methylation analysis

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2024207194A Division JP2025029007A (ja) 2019-07-08 2024-11-28 バイサルファイトを使用しない全ゲノムメチル化解析

Publications (3)

Publication Number Publication Date
JP2022540453A JP2022540453A (ja) 2022-09-15
JP2022540453A5 JP2022540453A5 (ja) 2023-07-14
JP7597783B2 true JP7597783B2 (ja) 2024-12-10

Family

ID=71895032

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2022500944A Active JP7597783B2 (ja) 2019-07-08 2020-07-08 バイサルファイトを使用しない全ゲノムメチル化解析
JP2024207194A Pending JP2025029007A (ja) 2019-07-08 2024-11-28 バイサルファイトを使用しない全ゲノムメチル化解析

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2024207194A Pending JP2025029007A (ja) 2019-07-08 2024-11-28 バイサルファイトを使用しない全ゲノムメチル化解析

Country Status (11)

Country Link
US (1) US12410467B2 (ja)
EP (2) EP3997245B1 (ja)
JP (2) JP7597783B2 (ja)
KR (1) KR20220038367A (ja)
CN (2) CN120138085A (ja)
AU (1) AU2020310613B2 (ja)
BR (1) BR112022000361A2 (ja)
ES (1) ES2969617T3 (ja)
IL (1) IL289669A (ja)
MX (1) MX2022000377A (ja)
WO (1) WO2021005537A1 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021161192A1 (en) 2020-02-11 2021-08-19 The Chancellor, Masters And Scholars Of The University Of Oxford Targeted, long-read nucleic acid sequencing for the determination of cytosine modifications
US20230357833A1 (en) * 2020-09-14 2023-11-09 Ludwig Institute For Cancer Research Ltd Cytosine modification analysis
CN117881795A (zh) * 2021-08-17 2024-04-12 因美纳有限公司 用于鉴定甲基化胞嘧啶的方法和组合物
JP2024540168A (ja) * 2021-11-02 2024-10-31 ガーダント ヘルス, インコーポレイテッド 品質管理方法
CN114085894A (zh) * 2021-11-03 2022-02-25 翌圣生物科技(上海)股份有限公司 核酸甲基化胞嘧啶转化方法
WO2023107453A1 (en) * 2021-12-06 2023-06-15 Chan Zuckerberg Biohub, Inc. Method for combined genome methylation and variation analyses
EP4499859A1 (en) * 2022-03-31 2025-02-05 Illumina, Inc. Compositions including aqueous amine borane complexes and polynucleotides, and methods of using the same to detect methylcytosine or hydroxymethylcytosine
CN115161400B (zh) * 2022-06-16 2023-06-20 温州医科大学 一种dna甲基化检测方法及所用试剂
EP4296372A1 (en) * 2022-06-22 2023-12-27 Universitat Pompeu Fabra (UPF) Method to detect and discriminate cytosine modifications
CN115820580A (zh) * 2022-08-31 2023-03-21 武汉爱博泰克生物科技有限公司 活性重组hTET2蛋白的表达及纯化方法
EP4594482A1 (en) * 2022-09-30 2025-08-06 Illumina, Inc. Cytidine deaminases and methods of use in mapping modified cytosine nucleotides
US20260103701A1 (en) * 2022-10-04 2026-04-16 Exact Sciences Innovation Ltd. Tet-assisted pyridine borane sequencing
CN115404270B (zh) * 2022-10-31 2023-01-31 臻和(北京)生物科技有限公司 一种dna甲基化测序文库甲基化转化率的评估方法、应用、终端设备及存储介质
KR20260005347A (ko) 2023-04-28 2026-01-09 이그잭트 사이언시스 이노베이션 리미티드 조작된 tet 효소 및 후성유전학 및 차세대 서열분석(ngs), 예컨대 tet-보조 피리딘 보란 서열분석(taps)에서의 용도
KR20250147902A (ko) * 2024-04-01 2025-10-14 (주) 제노텍 Dna 메틸화 분석 방법 및 그 활용
WO2025235365A1 (en) * 2024-05-08 2025-11-13 Freenome Holdings, Inc. Methods and systems for improved methylation sequencing

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7582420B2 (en) 2001-07-12 2009-09-01 Illumina, Inc. Multiplex nucleic acid reactions
WO2008091541A1 (en) * 2007-01-22 2008-07-31 University Of Vermont And State Agricultural College Molecular accessibility assay
US9115386B2 (en) 2008-09-26 2015-08-25 Children's Medical Center Corporation Selective oxidation of 5-methylcytosine by TET-family proteins
WO2011127136A1 (en) 2010-04-06 2011-10-13 University Of Chicago Composition and methods related to modification of 5-hydroxymethylcytosine (5-hmc)
EP2694686B2 (en) 2011-04-06 2023-07-19 The University of Chicago COMPOSITION AND METHODS RELATED TO MODIFICATION OF 5-METHYLCYTOSINE (5mC)
US20130157261A1 (en) * 2011-06-01 2013-06-20 The Methodist Hospital Research Institute Compositions and Methods for Quantitative Histology, Calibration of Images in Fluorescence Microscopy, and ddTUNEL Analyses
RU2612902C2 (ru) 2011-07-29 2017-03-13 Кембридж Эпидженетикс Лимитед Способы обнаружения модификации нуклеотидов
US9267117B2 (en) 2012-03-15 2016-02-23 New England Biolabs, Inc. Mapping cytosine modifications
CN104955960A (zh) 2012-11-30 2015-09-30 剑桥表现遗传学有限公司 用于修饰的核苷酸的氧化剂
WO2014165770A1 (en) 2013-04-05 2014-10-09 The University Of Chicago Single-base resolution sequencing of 5-formylcytosine (5fc) and 5-carboxylcytosine (5cac)
US20160194696A1 (en) 2013-08-09 2016-07-07 New England Biolabs, Inc. Detecting, Sequencing and/or Mapping 5-Hydroxymethylcytosine and 5-Formylcytosine at Single-Base Resolution
EP3280793A4 (en) 2015-04-06 2018-10-03 The Regents of the University of California Methods for determining base locations in a polynucleotide
WO2017039002A1 (ja) 2015-09-04 2017-03-09 国立大学法人東京大学 5-ヒドロキシメチルシトシン酸化剤及び5-ヒドロキシメチルシトシン解析方法
US11162139B2 (en) 2016-03-02 2021-11-02 Shanghai Epican Genetech Co. Ltd. Method for genomic profiling of DNA 5-methylcytosine and 5-hydroxymethylcytosine
WO2017176630A1 (en) 2016-04-07 2017-10-12 The Board Of Trustees Of The Leland Stanford Junior University Noninvasive diagnostics by sequencing 5-hydroxymethylated cell-free dna
EP3452591B1 (en) * 2016-05-02 2023-08-16 Encodia, Inc. Macromolecule analysis employing nucleic acid encoding
NZ765943A (en) 2018-01-08 2022-10-28 Ludwig Inst For Cancer Res Ltd Bisulfite-free, base-resolution identification of cytosine modifications
WO2019160994A1 (en) 2018-02-14 2019-08-22 Bluestar Genomics, Inc. Methods for the epigenetic analysis of dna, particularly cell-free dna
US20260103701A1 (en) 2022-10-04 2026-04-16 Exact Sciences Innovation Ltd. Tet-assisted pyridine borane sequencing

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DIZDAROGLU, M. et al,Substrate Specificity of the Escherichia coli Endonuclease III: Excision of Thymine- and Cytosine-Derived Lesions in DNA Produced by Radiation-Generated Free Radicals,Biochemistry,1993年,Vol.32,pp.12105-12111
LIU, Y. et al,Bisulfite-free direct detection of 5-methylcytosine and 5-hydroxymethylcytosine at base resolution,nature biotechnology,2019年,Vol.37,pp.424-429

Also Published As

Publication number Publication date
CN114174534B (zh) 2025-03-25
KR20220038367A (ko) 2022-03-28
WO2021005537A1 (en) 2021-01-14
CN114174534A (zh) 2022-03-11
AU2020310613A1 (en) 2022-02-17
EP3997245A1 (en) 2022-05-18
BR112022000361A2 (pt) 2022-05-10
IL289669A (en) 2022-03-01
US20220275424A1 (en) 2022-09-01
EP4306652A3 (en) 2024-03-27
EP3997245B1 (en) 2023-10-18
MX2022000377A (es) 2022-03-25
JP2022540453A (ja) 2022-09-15
ES2969617T3 (es) 2024-05-21
EP4306652B1 (en) 2026-03-25
EP4306652A2 (en) 2024-01-17
CN120138085A (zh) 2025-06-13
US12410467B2 (en) 2025-09-09
CA3146212A1 (en) 2021-01-14
AU2020310613B2 (en) 2026-02-05
JP2025029007A (ja) 2025-03-05

Similar Documents

Publication Publication Date Title
JP7597783B2 (ja) バイサルファイトを使用しない全ゲノムメチル化解析
JP7539948B2 (ja) シトシン修飾の、亜硫酸水素塩非含有、塩基分解能特定
CA3146212C (en) Bisulfite-free, whole genome methylation analysis
HK40103747A (en) Bisulfite-free, whole genome methylation analysis
HK40069820B (en) Bisulfite-free, whole genome methylation analysis
HK40069820A (en) Bisulfite-free, whole genome methylation analysis
EA047457B1 (ru) Модификация цитозина без применения бисульфита с разрешением до одного основания
NZ793130A (en) Bisulfite-free, base-resolution identification of cytosine modifications

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230706

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230706

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240521

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240821

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20241029

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241128

R150 Certificate of patent or registration of utility model

Ref document number: 7597783

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150