JP7597783B2 - バイサルファイトを使用しない全ゲノムメチル化解析 - Google Patents
バイサルファイトを使用しない全ゲノムメチル化解析 Download PDFInfo
- 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
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6813—Hybridisation assays
- C12Q1/6827—Hybridisation assays for detection of mutation or polymorphism
- C12Q1/683—Hybridisation assays for detection of mutation or polymorphism involving restriction enzymes, e.g. restriction fragment length polymorphism [RFLP]
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12P—FERMENTATION OR ENZYME-USING PROCESSES TO SYNTHESISE A DESIRED CHEMICAL COMPOUND OR COMPOSITION OR TO SEPARATE OPTICAL ISOMERS FROM A RACEMIC MIXTURE
- C12P19/00—Preparation of compounds containing saccharide radicals
- C12P19/26—Preparation of nitrogen-containing carbohydrates
- C12P19/28—N-glycosides
- C12P19/30—Nucleotides
- C12P19/34—Polynucleotides, e.g. nucleic acids, oligoribonucleotides
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6806—Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6813—Hybridisation assays
- C12Q1/6827—Hybridisation assays for detection of mutation or polymorphism
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Reaction characterised by the enzymatic activity
- C12Q2521/30—Phosphoric diester hydrolysing, i.e. nuclease
- C12Q2521/301—Endonuclease
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Reaction characterised by the enzymatic activity
- C12Q2521/50—Other enzymatic activities
- C12Q2521/531—Glycosylase
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Reactions characterised by treatment of reaction samples
- C12Q2523/10—Characterised by chemical treatment
- C12Q2523/115—Characterised by chemical treatment oxidising agents
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Reactions involving modified oligonucleotides, nucleic acids, or nucleotides
- C12Q2525/10—Modifications characterised by
- C12Q2525/117—Modifications characterised by incorporating modified base
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Reactions involving modified oligonucleotides, nucleic acids, or nucleotides
- C12Q2525/10—Modifications characterised by
- C12Q2525/191—Modifications characterised by incorporating an adaptor
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Oligonucleotides characterized by their use
- C12Q2600/154—Methylation 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の判定は、基礎研究に重要であるのみならず、診断及び治療を含む臨床用途にも有用である。
〔発明の概要〕
本開示は、全ゲノム塩基の分解能のメチローム解析方法を提供する。実施形態では、その方法は、DNA試料における5-メチルシトシン、5-ヒドロキシメチルシトシン、5-カルボキシルシトシン及び/または5-ホルミルシトシンのうちの1つ以上の位置を特定する。実施形態では、そのDNA試料は、全ゲノムを含む。本明細書に記載されている方法は、DNAのメチル化及びヒドロキシメチル化の解析であって、無修飾シトシンに影響を及ぼすことなく、塩基単位の分解能で、修飾シトシンを検出する穏和な反応を伴う解析を実現する。本発明で提供するのは、TETによる酸化と、ボラン誘導体(例えば、ピリジンボラン及び2-ピコリンボラン(pic-BH3))による還元を組み合わせることによって、5mC及び/または5hmCを特定するための改良型の方法であり、この方法は、本明細書では、TAPS(TETアシストピリジンボランシーケンシング)(表1)という。TAPSは、無修飾シトシンに影響を及ぼすことなく、高い感度及び特異性で、修飾を直接検出するとともに、本明細書に記載されているように、改変して、他のシトシン修飾を検出することができる。TAPSは、非破壊的であり、DNAを最長で10kbに保つものである。TAPSでは、バイサルファイトシーケンシングと比べて、マッピング率が高く、カバレッジが均一で、シーケンシングコストが低く、クオリティの上昇が可能となり、網羅性の向上した、より安価なメチローム解析が行われる。本明細書に記載されているように、5fC及び/または5caCを特定するには、酸化工程を用いない、この方法の変形形態を使用する。
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のいずれかの位置が得られる方法を提供する。
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の位置が得られる方法を提供する。
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のいずれかの位置が得られる方法を提供する。
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の位置が得られることと、
を含む方法を提供する。
a.5caC及び/または5fCを有する標的DNAを含むDNA試料を準備する工程と、
b.その5caC及び5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
c.その改変標的DNAを切断する工程と、
d.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
e.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCまたは5fCのいずれかの位置が得られる方法を提供する。
a.5caCを有する標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5fCに、ブロック基を付加する工程と、
c.その5caCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
d.その改変標的DNAを切断する工程と、
e.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
f.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCの位置が得られる方法を提供する。
a.5fCを有する標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5caCに、ブロック基を付加する工程と、
c.その5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
d.その改変標的DNAを切断する工程と、
e.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
f.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5fCの位置が得られる方法を提供する。
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の位置が得られる方法を提供する。
〔図面の簡単な説明〕
〔図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-ボラン反応である。1H NMRの結果は、2’-デオキシ-5,6-ジヒドロウリジンに関する過去の報告に従ったものである(I.Aparici-Espert et al.,J.Org.Chem.81,4031-4038(2016))。単一ヌクレオシドのpic-ボラン反応生成物の1H NMR(MeOH-d4、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-d4, 400 MHz)チャートである。δ ppm: 171.56 (CO), 153.54 (CO), 85.97 (CH), 83.86 (CH), 70.99 (CH), 61.92 (CH2), 36.04 (CH2), 35.46 (CH2), 30.49 (CH2).
〔図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-BH3によってジヒドロウラシル(DHU)に変換される。5fCは、O-エチルヒドロキシルアミン(EtONH2)のようなヒドロキシルアミン誘導体によってブロックしたが、その誘導体は、オキシムとなり、pic-ボランによる変換を阻止することになる。5caCは、EDC結合を介して、エチルアミンによってブロックして、pic-ボランによる変換をブロックするアミドに変換した。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。
〔図6-2〕5fC及び5caCを含むモデルDNAオリゴをpic-ボランによって、5fC及び5caCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5fC及び5caCは、pic-BH3によってジヒドロウラシル(DHU)に変換される。5fCは、O-エチルヒドロキシルアミン(EtONH2)のようなヒドロキシルアミン誘導体によってブロックしたが、その誘導体は、オキシムとなり、pic-ボランによる変換を阻止することになる。5caCは、EDC結合を介して、エチルアミンによってブロックして、pic-ボランによる変換をブロックするアミドに変換した。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。
〔図7-1〕5mC及び5hmCを含むモデルDNAオリゴをKRuO4及びpic-ボランによって、5hmCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5hmCは、βGTによってグルコースでブロックして、5gmCに変換できた。5mC、5hmC及び5gmCは、pic-ボランによって変換できなかった。5hmCは、KRuO4によって酸化させて、5fCにしてから、pic-ボランによってDHUに変換できた。各グラフの上に、MS(m/z)の計算値が示されており、ピークの左に、MSの観察値が示されている。
〔図7-2〕5mC及び5hmCを含むモデルDNAオリゴをKRuO4及びpic-ボランによって、5hmCをブロックした状態またはブロックしなかった状態で処理したものをMALDIによって特徴付けた結果である。5hmCは、βGTによってグルコースでブロックして、5gmCに変換できた。5mC、5hmC及び5gmCは、pic-ボランによって変換できなかった。5hmCは、KRuO4によって酸化させて、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のマッピングを行うものであり、バイサルファイトシーケンシングのような以前の方法の問題点を解消する。
別の態様では、本開示は、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試料における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の位置をもたらす工程と、
を含む方法を提供する。
本開示は、(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を含む試料から取った別々のアリコートであってよい。
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の位置が得られる方法を提供する。
一態様では、本開示は、DNA試料中の5caCまたは5fCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.5caC及び/または5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
f.その改変標的DNAを切断する工程と、
g.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
h.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCまたは5fCのいずれかの位置が得られる方法を提供する。
別の態様では、本開示は、DNA試料中の5caCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5fCに、ブロック基を付加する工程と、
c.その5caCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
g.その改変標的DNAを切断する工程と、
h.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
i.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5caCの位置が得られる方法を提供する。
別の態様では、本開示は、DNA試料中の5fCの特定方法であって、
a.標的DNAを含むDNA試料を準備する工程と、
b.そのDNA試料中の5caCに、ブロック基を付加する工程と、
c.その5fCをDHUに変換して、改変標的DNAを含む改変DNA試料をもたらす工程と、
d.その改変標的DNAを切断する工程と、
e.その切断された改変標的DNAに、アダプターDNA分子を付加する工程と、
f.その改変標的DNAの配列を検出する工程と、
を含み、
切断部位の存在により、その標的DNAにおける5fCの位置が得られる方法を提供する。
本開示の方法では、バイサルファイトシーケンシングのような方法に伴う大幅な分解を回避する穏和な酵素反応及び化学反応を利用する。すなわち、本発明の方法は、循環セルフリーDNAのような、インプットの少ない試料の解析、及びシングルセル解析に有用である。
本明細書に記載されている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)またはその誘導体である。
本明細書に記載されている方法は、DNA試料中の5caC及び/または5fCをDHUに変換する工程を含む。実施形態では、この工程は、DNA試料と、例えば、ボラン還元剤(ピリジンボラン、2-ピコリンボラン(pic-BH3)、ボラン、水素化ホウ素ナトリウム、シアノ水素化ホウ素ナトリウム及びトリアセトキシ水素化ホウ素ナトリウムなど)を含む還元剤とを接触させることを含む。好ましい実施形態では、その還元剤は、ピリジンボラン及び/またはpic-BH3である。
本明細書に記載されている方法は、変換工程(複数可)の前(すなわち、修飾シトシンをDHUに変換する1つまたは複数の工程の前)のDNA試料において修飾シトシン(5mC、5hmC、5fC及び/または5caC)があった位置に、DHUを含む改変標的DNAを切断する工程を含む。本明細書に記載されている切断工程は、DHUを含む改変標的DNAを特異的に切断する一方で、DHUを含まないDNAは、未切断または実質的に未切断のままにする。
実施形態では、本明細書に記載されている方法は、例えばUSER酵素によって切断した改変標的DNAに、アダプターDNA分子を付加する工程を含む。アダプターDNAまたはDNAリンカーは、化学合成した短い一本鎖または二本鎖のオリゴヌクレオチドであって、他のDNA分子の一方の末端または両方の末端にライゲーションできるオリゴヌクレオチドである。二本鎖のアダプターは、そのアダプターの各末端が、平滑末端、または5’もしくは3’オーバーハング(すなわち粘着末端)を有するように合成できる。DNAアダプターを、切断された改変標的DNAにライゲーションして、相補的なプライマーによるPCR増幅用、及び/またはクローニング及び/またはライブラリー作製(例えば次世代シーケンシングライブラリー)用の配列をもたらす。
本明細書に記載されている方法は、当該技術分野において知られている方法によって、改変標的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において、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-BH3)、ボラン、水素化ホウ素ナトリウム、シアノ水素化ホウ素ナトリウム及びトリアセトキシ水素化ホウ素ナトリウムからなる群のうちの1つ以上から選択されている。さらなる実施形態では、そのTET1酵素は、NgTet1、ヒトTET1またはマウスTet1であり、そのボラン還元剤は、ピリジンボラン及び/またはpic-BH3である。別の実施形態では、そのTET酵素は、mTET2またはその誘導体である。
〔実施例〕
実施例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)で精製した。
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCTCCCGAmCGCATGATCTGTACTTGATCGAChmCGTGCAACGATCGGAAGAGCACACGTCTGAACTCCAGTCACGCCAATATCTCGTATGCCGTCTTCTGCTTG
という配列(5’→3’)の最終生成物を得た。このモデルDNAの増幅用のPCRプライマーは、P5プライマーが5’-AATGATACGGCGACCACCGAG-3’、P7プライマーが5’-CAAGCAGAAGACGGCATACGAG-3’であった。
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCTAGGTGCGCTAAGTTCTAGATCGCCAACTGGTTGTGGCCTTAGCAGTCTmCGATCAGCTGmCTACTGTAmCGTAGCATCTATAGCCGGCTTGCTCTCTCTGCCTCTAGCAGCTGCTCCCTATAGTGAGTCGTATTAACGATCGGAAGAGCACACGTCTGAACTCCAGTCACGCCAATATCTCGTATGCCGTCTTCTGCTTG
であった。モデルDNAを増幅するためのPCRプライマーは、上に示したP5プライマー及びP7プライマーである。プライマー伸長用のビオチン標識プライマー配列は、このP7プライマーの5’末端に連結したビオチンである。T7 RNAポリメラーゼによって転写した後のRT-PCR用のPCRプライマーは、このP5プライマー、及び5’-TGCTAGAGGCAGAGAGAGCAAG-3’というRTプライマーであった。
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCTAGGTGCGCTAAGTTCTAGATCGCCAACTGGTTGTGGCCTTAGCAGTCTXGATCAGCTGCTACTGTACGTAGCATCTATAGCCGGCTTGCTCTCTCTGCCTCTAGCAGCTGCTCCCTATAGTGAGTCGTATTAACGATCGGAAGAGCACACGTCTGAACTCCAGTCACGCCAATATCTCGTATGCCGTCTTCTGCTTG
であり、この配列中、Xは、DHU、U、TまたはCである。このモデルDNAを増幅するためのPCRプライマーは、上記のP5プライマー及びP7プライマーである。
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で見ることができる。
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に列挙されている。
CACAGATGTCTGCCTGTTCATCCGCGTCCAGCTCGTTGAGTTTCTCCAGAAGCGTTAATGTCTGGCTTCTGATAAAGCGGGCCATGTTAAGGGCGGTTTTTTCCTGTTTGGTCACTGATGCCTCCGTGTAAGGGGGATTTCTGTTCATGGGGGTAATGATACCGATGAAACGAGAGAGGATGCTCACGATACGGGTTACTGATGATGAACATGCCCGGTTACTGGAACGTTGTGAGGGTAAACAACTGGCGGTATGGATGCGGCGGGACCAGAGAAAAATCACTCAGGGTCAATGCCAGCGCTTCGTTAATACAGATGTAGGTGTTCCACAGGGTAGCCAGCAGCATCCTGCGATGCAGATCCGGAACATAATGGTGCAGGGCGCTGACTTCCGCGTTTCCAGACTTTACGAAACACGGAAACCGAAGACCATTCATGTTGTTGCTCAGGTCGCAGACGTTTTGCAGCAGCAGTCGCTTCACGTTCGCTCGCGTATCGGTGATTCATTCTGCTAACCAGTAAGGCAACCCCGCCAGCCTAGCCGGGTCCTCAACGACAGGAGCACGATCATGCGCACCCGTGGGGCCGCCATGCCGGCGATAATGGCCTGCTTCTCGCCGAAACGTTTGGTGGCGGGACCAGTGACGAAGGCTTGAGCGAGGGCGTGCAAGATTCCGAATACCGCAAGCGACAGGCCGATCATCGTCGCGCTCCAGCGAAAGCGGTCCTCGCCGAAAATGACCCAGAGCGCTGCCGGCACCTGTCCTACGAGTTGCATGATAAAGAAGACAGTCATAAGTGCGGCGACGATAGTCATGCCCCGCGCCCACCGGAAGGAGCTGACTGGGTTGAAGGCTCTCAAGGGCATCGGTCGAGATCCCGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCGCTATCGGCTGAATTTGATTGCGAGTGAGATATTTATGCCAGCCAGCCAGACGCAGACGCGCCGAGACAGAACTTAATGGGCCCGCTAACAGCGCGATTTGCTGGTGACCCAATGCGACCAGATGCTCCACGCCCAGTCGCGTACCGTCTTCATGGGAGAAAATAATACTGTTGATGGGTGTCTGGTCAGAGACATCAAGAAATAACGCCGGAACATTAGTGCAGGCAGCTTCCACAGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCAGCCCACTGACGCGTTGCGCGAGAAGATTGTGCACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCT
120merのスパイクインコントロールの調製
120merのスパイクインコントロールをプライマー伸長法によって生成した。オリゴ配列及びプライマーは、表3に列挙されている。
N5mCNN及びN5hmCNNの配列を有する合成オリゴをアニーリング及び伸長方法によって生成した。オリゴ配列は、下記の表4に列挙されている。
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解析を行った。
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℃で保存した。
N末端Flagタグを有するmTET1CD触媒ドメイン(NM_001253857.2、4371~6392)をpcDNA3-Flagに、KpnI制限部位とBamH1制限部位の間にクローニングした。タンパク質の発現のために、1mgのプラスミドを1LのExpi293F(Gibco)細胞培養液に、1×106細胞/mLの密度でトランスフェクションし、細胞を48時間、37℃、170rpm及び5%CO2で成長させた。その後、細胞を遠心分離によって回収し、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によって検出されたことから、タンパク質がヌクレアーゼを含まないことが確認された。
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)によって、メーカーの指示に従って精製した。
Pic-BH3による還元:MeOH中の5Mのα-ピコリン-ボラン(pic-BH3、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)によって、メーカーの指示に従って精製した。
5hmCのブロック
50mMのHEPES緩衝液(pH8)、25mMのMgCl2、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)によって、メーカーの指示に従って精製した。
100mMのMES緩衝液(pH5.0)、10mMのO-エチルヒドロキシルアミン(Sigma-Aldrich)及び10μMの5fC DNAオリゴにおいて、5fCのブロックを2時間、37℃で行った。その生成物をMicro Bio-Spin 6 Columns(Bio-Rad)によって、メーカーの指示に従って精製した。
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 DNAオリゴ46μLを1MのNaOH2.5μLで30分、37℃で、振とう培養器において変性させてから、50mMのNaOH及び15mMの過ルテニウム酸カリウム(KRuO4、Sigma-Aldrich)を含む溶液1.5μLで1時間、氷上で酸化させた。その生成物をMicro Bio-Spin 6 Columnsによって、メーカーの指示に従って精製した。
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%超となるはずである。
TAPS前及びTAPS後のモデルDNA間で増幅曲線及び融解曲線を比較するために(図11)、1×LightCycler 480 High Resolution Melting Master Mix(Roche Diagnostics Corporation)、250nMのプライマー(フォワードプライマー:CCTGATGAAACAAGCATGTC及びリバースプライマー:CATTACTCACTTCCCCACTT)、ならびに3mMのMgSO4を含むPCRマスターミックス19μLに、DNA試料1ngを加えた。PCR増幅のために、初期変性工程を10分、95℃で行ってから、95℃での変性5分、カスタマイズしたアニーリング温度でのアニーリング5分、及び72℃での伸長5分というサイクルを40サイクル行った。最終工程には、95℃で1分、70℃で1分を含め、融解曲線(0.02℃刻み、毎回取得前に5分置いた)は65℃~95℃とした。
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産物を精製し、サンガーシーケンシング用に処理した。
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消化アッセイによって試験した。
ゼラチンコートプレートで、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%CO2で維持し、2日おきに継代した。
全ゲノムバイサルファイトシーケンシング(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用に、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のペアエンドシーケンシングを行った。
上記のようにして、全ゲノム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ビーズで精製した。
上記のようにして、全ゲノム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ビーズで精製した。
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ビーズで精製した。
ペアエンドリードを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)。参照により、本明細書に援用される)。
ペアエンドリードを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))。
カスタムの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)を用いて、改変塩基を検出した。
python3スクリプト(MF-phredder.py)によって、Illumina BaseSpaceからダウンロードしたような元のFASTQファイルから、ヌクレオチドの種類ごとのクオリティスコア統計を抽出した。
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を使用した。
塩基ごとの修飾リードの割合を、それぞれ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-BH3は、これまで知られていなかった還元的脱炭酸/脱アミノ反応によって、5fC及び5caCをDHUに容易に変換できることが見出された(図4)。その反応は、MALDIを用いて、単一のヌクレオシド及び複数のオリゴヌクレオチドの両方で、定量的であることが示された(図2~3及び6-1~7-2)。
方法
スパイクインコントロールの調製
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を調製した。
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の断片について選択した。
ゼラチンコートプレートで、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を抽出した。
mTet1触媒ドメイン(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ビーズで、メーカーの指示に従って精製した。
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%アガロースゲルに流した。
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のペアエンドシーケンシングを行った。
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のペアエンドシーケンシングを行った。
生のシーケンシングリードを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の割合によってサブサンプリングした。
生のシーケンシングリードをseqtk(https://github.com/lh3/seqtk)のtrimfq -b 2で処理して、各リードの左から2bpをトリミングした。Astair3.2.7を使用して、rrTAPSを処理した(8)。mm9ゲノムを参照として使用して、astair alignを用いて、クリーンしたリードをアラインメントした。メチル化CpGをastairコールによって抽出した。
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を選択した。
モデルDNAに対するeeTAPSの能力を示したので、eeTAPSを用いて、mESC gDNAにおけるCpGのメチル化をプロファイリングした(表11)。eeTAPSは、メチル化CpGを濃縮することになるので、コスト効率の良い方法であることが提議されている。実際、eeTAPSにおける断片の84.6%がC/Gとなったことを我々は見出した(図31A)。切断部位と、最も近いCpGとの距離に対するさらなる解析により、切断イベントの72.7%が、CpGで生じたことが特定された(図31B)。
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)。
インプットの少ない試料での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の適用が実行可能であることがさらに示された。
シーケンシング深度が、検出できる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の分解能での直接的なメチル化検出の可能性をもたらす。
Claims (16)
- 改変標的DNAの切断方法であって、
標的DNA中の5caC及び/または5fCをジヒドロウラシル(DHU)に変換し、1つ以上のDHUを含む改変標的DNAを提供する変換工程と、
1つ以上のDHUを含む前記改変標的DNAと、前記1つ以上のDHUの位置または前記1つ以上のDHUと隣接する位置で、前記改変標的DNAを切断する1つ以上のエンドヌクレアーゼとを接触させることとを含む、方法。 - 前記方法が、前記変換工程の前に、5mC及び/または5hmCを含む標的DNAを改変することをさらに含み、前記改変することが、前記標的DNAにおける5mC及び/または5hmCを5-カルボキシルシトシン(5caC)及び/または5-ホルミルシトシン(5fC)に変換する工程を含む、請求項1に記載の方法。
- 前記方法が、前記5mC及び/または5hmCを変換する工程の前に、5mC及び5hmCを含む前記標的DNAを改変することをさらに含み、前記改変することが、前記標的DNAにおける5hmCに、ブロック基を付加する工程を含む、請求項2に記載の方法。
- 前記5hmCに付加する前記ブロック基が、糖である、請求項3に記載の方法。
- 前記標的DNAにおける5mC及び/または5hmCを5caC及び/または5fCに変換する工程が、前記標的DNAと、10-11転座(TET)酵素とを接触させることを含む、請求項2から4のいずれか1項に記載の方法。
- 前記TET酵素が、ヒトTET1、ヒトTET2、ヒトTET3、マウスTET1、マウスTET2、マウスTET3、Naegleria TET(NgTET)、およびCoprinopsis cinerea TET(CcTET)からなる群から選択される、請求項5に記載の方法。
- 前記TET酵素が、NgTET、マウスTET2、またはヒトTET1である、請求項5または6に記載の方法。
- 前記5hmCを5caC及び/または5fCに変換する工程が、前記標的DNAを酸化剤と接触させることを含む、請求項2から4のいずれか1項に記載の方法。
- 前記酸化剤が、過ルテニウム酸カリウム、Cu(II)/TEMPO、ルテニウム酸カリウムまたは酸化マンガンである、請求項8に記載の方法。
- 前記5caC及び/または5fCをDHUに変換する工程が、前記標的DNAを還元剤と接触させることを含む、請求項1から9のいずれか1項に記載の方法。
- 前記還元剤が、ピリジンボラン、2-ピコリンボラン(pic-BH3)、ボラン、水素化ホウ素ナトリウム、シアノ水素化ホウ素ナトリウム及びトリアセトキシ水素化ホウ素ナトリウムからなる群から選択される、請求項10に記載の方法。
- 前記1つ以上のエンドヌクレアーゼが、アプリン酸/アピリミジン酸エンドヌクレアーゼ1(APE1)、ウラシルDNAグリコシラーゼ(UDG)、エンドヌクレアーゼIII、TmaエンドヌクレアーゼIII、TthエンドヌクレアーゼIV、エンドヌクレアーゼV、エンドヌクレアーゼVIII、ホルムアミドピリミジンDNAグリコシラーゼ(Fpg)及びhNEIL1からなる群のうちの1つ以上から選択される、請求項1から11のいずれか1項に記載の方法。
- 前記1つ以上のエンドヌクレアーゼが、ウラシルDNAグリコシラーゼ(UDG)及びDNAグリコシラーゼ-リアーゼエンドヌクレアーゼVIIIを組み合わせたものを含む、請求項1から12のいずれか1項に記載の方法。
- 前記切断された改変標的DNAに、アダプターDNA分子を付加する工程をさらに含む、請求項1から13のいずれか1項に記載の方法。
- 前記切断された改変標的DNAの配列を検出する工程をさらに含み、切断された部位の存在は、前記標的DNAにおける5caCまたは5fCの何れかの位置を提供する、請求項1から14のいずれか1項に記載の方法。
- 前記方法が、前記接触させる工程の前に、前記改変標的DNAのコピー数を増幅する工程をさらに含む、請求項1から15のいずれか1項に記載の方法。
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)
| 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)
| 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 |
-
2020
- 2020-07-08 EP EP20749939.3A patent/EP3997245B1/en active Active
- 2020-07-08 KR KR1020227004136A patent/KR20220038367A/ko active Pending
- 2020-07-08 CN CN202510225967.8A patent/CN120138085A/zh active Pending
- 2020-07-08 ES ES20749939T patent/ES2969617T3/es active Active
- 2020-07-08 MX MX2022000377A patent/MX2022000377A/es unknown
- 2020-07-08 EP EP23203846.3A patent/EP4306652B1/en active Active
- 2020-07-08 US US17/625,500 patent/US12410467B2/en active Active
- 2020-07-08 CN CN202080053099.5A patent/CN114174534B/zh active Active
- 2020-07-08 WO PCT/IB2020/056435 patent/WO2021005537A1/en not_active Ceased
- 2020-07-08 AU AU2020310613A patent/AU2020310613B2/en active Active
- 2020-07-08 BR BR112022000361A patent/BR112022000361A2/pt unknown
- 2020-07-08 JP JP2022500944A patent/JP7597783B2/ja active Active
-
2022
- 2022-01-06 IL IL289669A patent/IL289669A/en unknown
-
2024
- 2024-11-28 JP JP2024207194A patent/JP2025029007A/ja active Pending
Non-Patent Citations (2)
| 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 |