JP2019514360A - 超並列シークエンシングのためのdnaライブラリーを生成する方法及びキット - Google Patents

超並列シークエンシングのためのdnaライブラリーを生成する方法及びキット Download PDF

Info

Publication number
JP2019514360A
JP2019514360A JP2018553460A JP2018553460A JP2019514360A JP 2019514360 A JP2019514360 A JP 2019514360A JP 2018553460 A JP2018553460 A JP 2018553460A JP 2018553460 A JP2018553460 A JP 2018553460A JP 2019514360 A JP2019514360 A JP 2019514360A
Authority
JP
Japan
Prior art keywords
sequencing
primer
seq
section
library
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.)
Granted
Application number
JP2018553460A
Other languages
English (en)
Other versions
JP6899844B2 (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 JP2019514360A publication Critical patent/JP2019514360A/ja
Application granted granted Critical
Publication of JP6899844B2 publication Critical patent/JP6899844B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12NMICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09Recombinant DNA-technology
    • C12N15/10Processes for the isolation, preparation or purification of DNA or RNA
    • C12N15/1034Isolating an individual clone by screening libraries
    • C12N15/1082Preparation or screening gene libraries by chromosomal integration of polynucleotide sequences, HR-, site-specific-recombination, transposons, viral vectors
    • 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
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/6853Nucleic acid amplification reactions using modified primers or templates
    • C12Q1/6855Ligating adaptors
    • 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/155Modifications characterised by incorporating/generating a new priming site
    • 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
    • C12Q2535/00Reactions characterised by the assay type for determining the identity of a nucleotide base or a sequence of oligonucleotides
    • C12Q2535/122Massive parallel 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
    • C12Q2563/00Nucleic acid detection characterized by the use of physical, structural and functional properties
    • C12Q2563/143Magnetism, e.g. magnetic label
    • 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
    • C12Q2563/00Nucleic acid detection characterized by the use of physical, structural and functional properties
    • C12Q2563/149Particles, e.g. beads
    • 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
    • C12Q2563/00Nucleic acid detection characterized by the use of physical, structural and functional properties
    • C12Q2563/179Nucleic acid detection characterized by the use of physical, structural and functional properties the label being a nucleic acid

Abstract

a)WGAライブラリー汎用配列アダプターを含むフラグメントを含む一次WGA DNAライブラリーを準備する工程と、b)一次WGA DNAライブラリーを、少なくとも1つの第1プライマー及び少なくとも1つの第2プライマーを用いて再増幅する工程とを含み、少なくとも1つの第1プライマーが、少なくとも1つの第1シークエンシングアダプターと、少なくとも1つの第1シークエンシングバーコードと、WGAライブラリー汎用配列アダプター又はその逆相補鎖のいずれかにハイブリダイズする第1プライマー3’セクションとを含み、少なくとも1つの第2プライマーが、少なくとも1つの第1シークエンシングアダプターとは異なる少なくとも1つの第2シークエンシングアダプターと、WGAライブラリー汎用配列アダプター又はその逆相補体のいずれかにハイブリダイズする第2プライマー3’セクションとを含む、超並列シークエンシングライブラリー生成方法。【選択図】図3

Description

本発明は、全ゲノム増幅(WGA:Whole Genome Amplification)産物の全ゲノムシークエンシングのための超並列シークエンシングライブラリーを生成する方法及びキットに関する。特に、該方法は、確定的制限部位全ゲノム増幅(DRS−WGA:Deterministic Restriction-Site, Whole Genome Amplification)のDNA産物に適用することもできる。
このライブラリーは、ローパス全ゲノムシークエンシング及びゲノムワイドコピー数プロファイリングに有利に使用することができる。
単一細胞では、単純化のため及び/又はシークエンシング、SNP検出等を含む異なるタイプの遺伝分析を行うことを可能にするために、より多くのDNAが得られる全ゲノム増幅(WGA)を行うことが有益である。
確定的制限部位に基づくLM−PCRによるWGA(例えば特許文献1に記載される)が、当該技術分野で既知である(以下、単にDRS−WGAと称する)。DRS−WGAは、単一細胞の増幅により良好な解決法であり(非特許文献1を参照)、更には固定によるDNA分解により強い(非特許文献2を参照)ことが実証されている。
LM−PCRによるDRS−WGAの市販のキット(Ampli1(商標) WGAキット、Silicon Biosystems)が非特許文献3において使用されている。この研究では、単一細胞WGA材料に対してローパス全ゲノムシークエンシングによるコピー数分析が行われた。しかしながら、非特許文献3で用いられる標準ワークフローについては、i)WGAアダプターの消化、ii)DNA断片化、並びにiii)末端修復(EndRepair)、iv)A−テーリング、v)バーコード化アダプターライゲーション等の標準Illuminaワークフロー工程、加えてvi)バーコード化NGSライブラリーのサンプルプーリング及びvii)シークエンシングという通常の工程を含む幾つかの工程がIlluminaライブラリーの作製に必要とされた。非特許文献3(図5b)において示されるように、WBCは偽陽性と推定されるコピー数コールを殆ど提示しなかったが、CTCは概して、はるかに多くの異常を示した。
Ampli1(商標) WGAは、アレイ比較ゲノムハイブリダイゼーション(aCGH:array Comparative Genomic Hybridization)に適合する。実際に、幾つかのグループ(非特許文献4、非特許文献5)により、それが高分解能コピー数分析に好適であることが示されている。しかしながら、aCGH法は高コストかつ労働集約的であり、体細胞性コピー数変化(CNA)の検出のためのローパス全ゲノムシークエンシング(LPWGS)等の異なる方法が望ましい可能性がある。
Baslan et al(非特許文献6)は、DOP−PCR全ゲノム増幅から出発し、WGAアダプター消化、Illuminaアダプターのライゲーション、PCR増幅を含む幾つかの酵素工程を用いて全ゲノムコピー数プロファイリングを達成した。
非特許文献7は、染色体異常及び単一遺伝子疾患についての同時の着床前遺伝子診断へのMALBAC WGA(Yikon Genomics Inc)の使用を教示している。
特許文献2は、特別な縮重オリゴヌクレオチドプライミングPCR(DOP−PCR:Degenerate-Oligonucleotide-Priming-PCR)を採用するアプローチを教示している。特許文献2には種々のWGA法及び現行の技術水準の概要も掲載されている。特許文献2は、市販のキットPicoPlexのベースとなっている。
非特許文献8は、MALBAC及び多重置換増幅(MDA:Multiple Displacement Amplification)を含む幾つかのWGA法の性能比較を報告している。LPWGS及びWGSが非特許文献8に用いられている。ライブラリー調製がワークフローにより行われる。
DRS−WGAはアレイCGH(比較ゲノムハイブリダイゼーション)、中期CGH、及びヘテロ接合性の消失等の他の遺伝分析アッセイを用いた場合に、微量の顕微解剖FFPE材料からのコピー数プロファイルの分析にDOP−PCRよりも良好であることが示されている(非特許文献2、非特許文献9)。
特許文献3は、突然変異を導入、除去した、又は制限部位を変更した遺伝子座においてDRS−WGA産物から突然変異を検出する方法を教示している。
特許文献4は、非特許文献10に更に詳述及び報告されるような、細胞のゲノム完全性及び/又はDRS−WGA産物の品質を多重PCRにより評定する方法を教示している。
DRS−WGAは、一様かつバランスの取れた増幅に関して最良の結果をもたらすが、aCGH又は中期CGHに基づく現行のプロトコルは、困難及び/又は高コストである。ローパス全ゲノムシークエンシングが、幾つかのサンプルをaCGHよりも高い処理能力及び低いコストで分析するハイスループット方法として提案されている。しかしながら、WGA産物(DRS−WGA等)の超並列シークエンシングライブラリーを生成する既知の方法では、幾つかの酵素工程及び反応を含むプロトコルが依然として必要とされる。
上記に引用されるCTC分析への適用以外に、胚盤胞に対する出生前診断等の他の単一細胞分析用途、及び母体血から採取された循環胎児細胞についても、DRS−WGAの再現性及び品質と、ゲノムワイドコピー数多型(CNV:Copy-Number Variant)を分析する能力とを兼ね備えた、より簡素化された方法を行うことが望ましい場合がある。加えて、微量の細胞、FFPE又は組織生検材料から全ゲノムコピー数プロファイルを決定することも望ましい場合がある。
特許文献5は、ステムループアダプターオリゴを断片化ゲノムDNAに付加することを含む、シークエンシングのためのライブラリーを調製する方法を開示している。次いで、ループを切断することで、二本鎖アダプターによって隣接したゲノムフラグメントを生じる。次いで、バーコードを含むプライマーによりフラグメントを増幅し、Ion TorrentシークエンシングプラットホームでのDNAシークエンシングに使用する。
この方法は一連の欠点を有し、以下が最も重要な欠点である:
方法は、幾つかの反応及び幾つかの酵素が関わる多数の後続工程を伴う。
方法は、それ自体で単一細胞サンプルに由来するDNAに適用可能でない。
国際公開第2000/017390号 米国特許第8206913号(Kamberov et al, Rubicon Genomics) 国際公開第2014068519号(Fontana et al.) 国際公開第2015083121号(Klein et al.) 国際公開第2014/071361号
Lee YS, et al: Comparison of whole genome amplification methods for further quantitative analysis with microarray-based comparative genomic hybridization. Taiwan J Obstet Gynecol. 2008, 47(1):32-41. Stoecklein et al., SCOMP is superior to degenerated oligonucleotide primed-polymerase chain reaction for global amplification of minute amounts of DNA from microdissected archival tissue samples, Am J Pathol. 2002 Jul;161(1):43-51 Hodgkinson C.L. et al., Tumorigenicity and genetic profiling of circulating tumor cells in small-cell lung cancer, Nature Medicine 20, 897-903 (2014) Moehlendick B, et al. (2013) A Robust Method to Analyze Copy Number Alterations of Less than 100 kb in Single Cells Using Oligonucleotide Array CGH. PLoS ONE 8(6): e67031 Czyz ZT, et al (2014) Reliable Single Cell Array CGH for Clinical Samples. PLoS ONE 9(1): e85907 Optimizing sparse sequencing of single cells for highly multiplex copy number profiling, Genome Research, 25:1-11, April 9, 2015 Yan et al. Proc Natl Acad Sci U S A. 2015 Dec 29;112(52):15964-9 Hou et al., Comparison of variations detection between whole-genome amplification methods used in single-cell resequencing, GigaScience (2015) 4:37 Arneson et al., Comparison of whole genome amplification methods for analysis of DNA extracted from microdissected early breast lesions in formalin-fixed paraffin-embedded tissue, ISRN Oncol. 2012;2012:710692. doi: 10.5402/2012/710692. Epub 2012 Mar 14. Polzer et al. EMBO Mol Med. 2014 Oct 30;6(11):1371-86
本発明の一目的は、NGS(Next Generation Sequencing;次世代シークエンシング)ライブラリーをWGA産物から出発して簡素化された形で生成する方法を提供することである。特に、文献中に一般に報告されるよりも少ない酵素反応を含む方法を提供することが本発明の一目的である。
本発明の別の目的は、本発明によるライブラリー調製方法を用いてゲノムワイドコピー数プロファイルをWGA産物から出発して生成する方法を提供することである。
本発明の更なる目的は、上述の方法を行うためのキットを提供することである。作製されるライブラリーは選択シークエンシングプラットホーム、例えばIon Torrentプラットホーム又はIlluminaプラットホームに適合するのが好ましい。
本発明は、添付の特許請求の範囲に規定されるような全ゲノム増幅産物からの全ゲノムシークエンシングのための超並列シークエンシングライブラリーを生成する方法及びキットに関する。本発明は、本発明の方法により先に調製したライブラリーを用い、WGA産物から出発してゲノムワイドコピー数プロファイルを生成する方法に更に関する。
プライマー配列及び操作プロトコルも提供する。
ライブラリー生成反応は、幾つかのサンプルを同じNGSランで多重化するためのシークエンシングバーコードの導入を含むのが好ましい。好ましくは、WGAはDRS−WGAであり、ライブラリーを単一チューブ一段階PCR反応により生成する。
DRS−WGAにより生成されたDNAライブラリーに含まれる、本発明の第1の実施形態に使用される出発産物を示す図である。単一フラグメントを概略的にのみ示す。 MALBACにより生成されたDNAライブラリーに含まれる、本発明の第2の実施形態に使用される出発産物を示す図である。単一フラグメントを概略的にのみ示す。 図1に示されるようなDRS−WGAにより生成されたDNAライブラリーのフラグメントに適用され、Ion Torrent又はIlluminaシークエンシングプラットホームのようなシークエンシングプラットホームに適合するDNAライブラリーをもたらすように指向された、本発明による方法の再増幅工程の一実施形態を概略的に示す図である。 図1に示されるようなDRS−WGAのフラグメントに適用される、本発明に従って得られる再増幅反応工程、続いてフラグメントライブラリー選択を含むプロトコルワークフローを概略的に示す図である。本方法は、ILLUMINAシークエンシングプラットホームに適合するDNAライブラリーを直接提供する。 図4に示す工程に続いてDRS−WGAのフラグメントに適用される、本発明の方法の第3の実施形態に従って得られる最終一本鎖DNAライブラリーを概略的に示す図である。さらに、シークエンシングされた最終ssDNAライブラリー及び本発明に従って設計されたカスタムシークエンシングプライマーを図5に示す。DEPArrayシステム(Bolognesi et al.)によりFFPEからデジタル処理で選別された数百個の腫瘍細胞から出発して、DRS−WGAライブラリーが生成される。 本発明に従って調製され、PGMプラットホームによってシークエンシングされたDNAライブラリーに対してDEPArrayシステムによりFFPEからデジタル処理で選別された数百個の腫瘍細胞から出発して行われたローパス全ゲノムシークエンシングのシークエンシング結果を示す図である。 2つの異なる腫瘍細胞での本発明に従って調製されたDNAライブラリーに対するPGMプロトコルによって行ったローパス全ゲノムシークエンシングのシークエンシング結果を示す図である。 本発明に従って調製されたDNAライブラリーに対するILLUMINAプロトコル1によって行ったローパス全ゲノムシークエンシングのシークエンシング結果を示し、正常WBC細胞及び異常(腫瘍)細胞から得られた結果を比較する図である。 本発明に従って調製されたDNAライブラリーに対する本発明の一態様に従うILLUMINAプロトコル2によって行ったローパス全ゲノムシークエンシングのシークエンシング結果を示す図である。
定義
他に定義されていなければ、本明細書に使用される全ての技術用語及び科学用語は、本発明が関連する技術分野の当業者により一般に理解されるものと同じ意味を有するものである。本明細書に記載されるものと同様な又は均等な多くの方法及び材料を本発明の実施又は試験に使用することができるが、好ましい方法及び材料を下に記載している。他に言及されていなければ、本発明に使用される本明細書に記載の技法は当業者に既知の標準的な方法である。
「消化部位(DS:Digestion Site)」又は「制限部位(RS:Restriction Site)」という用語は、制限酵素がポリヌクレオチド鎖に沿って切り取る、制限酵素によって認識されるDNA分子に沿ったヌクレオチドの配列(通例、4塩基対〜8塩基対(bp)長)を意図するものである。
「切断部位」という用語は、制限酵素がヌクレオチド間のホスホジエステル結合を加水分解することによってヌクレオチドを切断する、ポリヌクレオチド鎖中の部位を意図するものである。
「アンプリコン」という用語は、PCR増幅によって生じるDNAの領域を意図するものである。
「DRS−WGAアンプリコン」又は短縮して「WGAアンプリコン」という用語は、ライゲートしたアダプターによって隣接した2つのRS間のDNA配列を含む、DRS−WGA中に増幅されたDNAフラグメントを意図するものである。
「元のDNA」という用語は、DRS−WGAによる増幅前のゲノムDNA(gDNA)を意図するものである。
「アダプター」又は「WGAアダプター」又は「WGA PCRプライマー」又は「WGAライブラリー汎用配列アダプター」という用語は、DRS−WGAの場合には、制限酵素の作用によって生成する各フラグメントにライゲートした付加的なオリゴヌクレオチド、又はMALBACの場合には、伸長及びPCRプロセスの結果としてWGA DNAライブラリーの各分子の5’セクションに存在する既知のポリヌクレオチド配列を意図するものである。
「コピー数変化(CNA:Copy Number Alteration)」という用語は、概して同じ個体のゲノムに関して規定される、ゲノム領域のコピー数の体細胞性変化を意図するものである。
「コピー数多型(CNV:Copy Number Variation)」という用語は、概して参照ゲノムに関して規定される、ゲノム領域のコピー数の生殖細胞系列変異を意図するものである。本明細書全体を通して、CNA及びCNVは、理論の殆どをどちらの状況にも当てはめることができるため、区別なく使用される場合がある。逆の指定がない限り、これらの用語の各々が両方の状況を指すことが意図されるものとする。
「超並列(Massive-parallel)次世代シークエンシング(NGS)」という用語は、空間的及び/又は時間的に分離され、クローン的に(clonally)シークエンシングされた(事前のクローン増幅(clonal amplification)を含む又は含まない)、DNA分子のライブラリーの作製を含む、DNAをシークエンシングする方法を意図するものである。例としては、Illuminaプラットホーム(Illumina Inc)、Ion Torrentプラットホーム(ThermoFisher Scientific Inc)、Pacific Biosciencesのプラットホーム、MinION(Oxford Nanopore Technologies Ltd)が挙げられる。
「標的配列」という用語は、元のDNA上の対象の領域を意図するものである。
「一次WGA DNAライブラリー(pWGAlib:Primary WGA DNA library)」という用語は、WGA反応によって得られるDNAライブラリーを意図するものである。
「多重アニーリング及びループ化による増幅サイクル法(MALBAC:Multiple Annealing and Looping Based Amplification Cycles)」という用語は、準線形全ゲノム増幅法を意図するものである(Zong et al., Genome-wide detection of single-nucleotide and copy-number variations of a single human cell, Science. 2012 Dec 21;338(6114):1622-6. doi: 10.1126/science.1229164.)。MALBACプライマーは、鋳型にハイブリダイズする8ヌクレオチドの3’ランダム配列、及び27ヌクレオチドの5’共通配列(GTG AGT GAT GGT TGA GGT AGT GTG GAG)を有する。最初の伸長後にセミアンプリコン(semiamplicons)を、相補的な5’末端及び3’末端を有する完全なアンプリコンをもたらす別の伸長の鋳型として使用する。数サイクルの準線形増幅の後、完全なアンプリコンを後続のPCRサイクルにより指数関数的に増幅することができる。
「DNAライブラリー精製」という用語は、DNAライブラリー材料を酵素、dNTP、塩、及び/又は所望のDNAライブラリーの一部ではない他の分子等の不要な反応成分から分離するプロセスを意図するものである。DNAライブラリー精製プロセスの例は、Beckman CoulterのAgencourt AMPure XPビーズ若しくは固相可逆的固定化(SPRI)ビーズ等の磁性ビーズを用いた技術、又はMerck MilliporeのAmiconスピンカラム等のスピンカラム精製による精製である。
「DNAライブラリーサイズ選択」という用語は、DNAライブラリーを構成する種々のフラグメントの塩基対分布を変更するプロセスを意図するものである。概して、或る特定の範囲内に含まれるDNAライブラリーの一部が実質的に保持されるのに対し、その範囲外のDNAライブラリー成分は実質的に廃棄される。DNAライブラリーサイズ選択プロセスの例は、電気泳動ゲル(例えば、ThermoFisher ScientificのE−gel)の切出し、又は磁性ビーズを用いた精製システム(例えば、Beckman CoulterのSPRIビーズ)による二重精製である。
「DNAライブラリー選択」という用語は、DNAライブラリー精製若しくはDNAライブラリーサイズ選択のいずれか、又はその両方を行うプロセスを意図するものである。
「NGS再増幅」という用語は、一次WGA DNAライブラリーの全て又は大部分が更に増幅されるPCR反応を意図するものである。NGSという用語は、本書全体を通して簡略化のために省略される場合があり、単に「再増幅」を参照されたい。
「シークエンシングアダプター(SA)」という用語は、DNA挿入断片のシークエンシングに役立つ1つ以上の分子を意図するものであり、各分子はポリヌクレオチド配列、官能基を含まないか、1つ以上含む場合がある。特に、シークエンサーがアウトプット配列を正確に生成するために超並列シークエンシングライブラリー中に存在する必要があるが、情報を保有しないポリヌクレオチド配列が意図される(非限定的な例:ssDNAをIlluminaシークエンシングの場合にはフローセル、若しくはIon Torrentシークエンシングの場合にはion−sphereにハイブリダイズさせるポリヌクレオチド配列、又は合成時解読(sequencing-by-synthesis)反応の開始に必要とされるポリヌクレオチド配列)。
「シークエンシングバーコード」という用語は、1つのシークエンサーリード内でシークエンシングする場合に、そのリードをそのバーコードと関連付けられた特定のサンプルに割り当てることを可能にするポリヌクレオチド配列を意図するものである。
「選択シークエンシングプラットホームに機能的な」という用語は、シークエンシングプロセス中にシークエンシングプラットホームに用いる必要があるポリヌクレオチド配列(例えば、バーコード又はシークエンシングアダプター)を意図するものである。
「ローパス全ゲノムシークエンシング」という用語は、1未満の平均シークエンシング深度での全ゲノムシークエンシングを意図するものである。
「平均シークエンシング深度」という用語は本明細書において、全参照ゲノムサイズで除算した、1サンプル当たりの、シークエンシングされ、参照ゲノムにマッピングされる塩基の総数を意図するものである。シークエンシング及びマッピングされる塩基の総数は、マッピングされるリードの数×平均リード長に近似することができる。
「二本鎖DNA(dsDNA)」という用語は、塩基対合則(AとT及びCとG)に従って、両者の窒素塩基の結合によって水素結合した2つの別個のポリヌクレオチド相補鎖を意図するものである。一本鎖DNA(ssDNA):DNAの2つの鎖は、2つの一本鎖DNA分子を形成し得る。すなわち、DNA分子はワトソン−クリック塩基対合によって連結した2つのssDNA分子で構成される。
「一本鎖DNA(ssDNA)」という用語は、例えば二本鎖DNAに由来する、又は相補的な一本鎖DNAと対合することができるポリヌクレオチド鎖、すなわち典型的ならせん型の2つのヌクレオチド鎖とは対照的に、一本鎖のみからなるポリヌクレオチドDNA分子を意図するものである。
「均等化する(equalizing)」は、1つ以上のサンプルの濃度を均等にするために調整する行為を意図するものである。
「正規化する(normalizing)」は、1つ以上のサンプルの濃度を所望の割合に対応するように調整する行為を意図するものである(均等化は、その割合が1である特別な例である)。本明細書では便宜上、「正規化する」及び「均等化する」という用語は、明らかに概念上同一であることから、区別せずに使用される。
「常磁性ビーズ」は、ストレプトアビジンコンジュゲート磁性ビーズ(例えば、Dynabeads(商標) MyOne(商標) Streptavidin C1、ThermoFisher Scientific)を意図するものである。常磁性ビーズのインキュベーションに言及する場合の「設計条件」という表現は、ストレプトアビジンコンジュゲート磁性ビーズを以下のバッファー:10mM Tris−HCl(pH7.5)、1mM EDTA、2M NaClで2回洗浄することに含まれる(consists in)活性化工程に必要とされる条件を指す。
ワークフロー
以下の表に本発明による幾つかの考え得るワークフローをまとめる。
プロセスインプット材料
本明細書全体で一次WGA DNAライブラリーに言及する。同じワークフローを、例えばdsDNA合成、又は標準WGAプライマーを用いた(例えば、可能な場合はAmpli1(商標) ReAmp/dsキット(Menarini Silicon Biosystems spa,Italy)を用いた)ライブラリー再増幅等の付加的なプロセスに更に供した一次WGA DNAライブラリーに適用することができる。便宜上、本明細書においては、これらの付加的なプロセスを考慮することなく一次WGA DNAライブラリーのみに言及する。この種の全てのインプットサンプルを、特許請求の範囲に報告されているものについても好適なサンプルインプットとして使用することができることが意図されるものとする。
初期精製
無視することができない量の一次WGAプライマーが一次WGAアウトプット産物中に存在する場合、DNAライブラリー精製を含む初期DNAライブラリー選択を行うことが有利である場合がある。実際、本発明によるプライマーは、一次WGAプライマー中に見られる共通配列に対応する配列を3’末端に有するため、無視することができない量の残留一次WGAプライマーの存在は、本発明による超並列シークエンシングライブラリーを得るために使用される再増幅プライマーと競合し、所望に応じて再増幅プライマー(複数の場合もある)(又はそれらの逆相補鎖(reverse complementary))を両末端に有するDNAライブラリー分子の収率を減少させる可能性がある。
サンプルにわたるリード数の均等化のための定量化
再増幅DNAライブラリーの量の変動が、共にプール及びシークエンシングされるサンプル間で比較的大きい場合、ライブラリーを等分し、各サンプルについてシークエンシングされるリード数を均等化するために、各サンプルからのDNAライブラリーの量を定量化することが有利である場合がある。
シークエンサーリード長及びWGAサイズピークの間のミスマッチは不正確な均等化を生じ得る
幾つかの超並列シークエンサー(Ion Torrentプラットホーム及びIlluminaプラットホームを含む)は、用いられる種々の化学作用に応じて50bp〜800bpに含まれるサイズ分布ピークを有するDNAフラグメント、例えば150bp、200bp、400bp、650bpに分布ピーキング(peaking)を有するDNAフラグメントのシークエンシングを用いる。pWGAlibサイズ分布は約1kbp等のより大きなフラグメントのピーク及び150bp、200bp、400bpのはるかに少量のDNAを有し得るため、所望の範囲の再増幅DNAライブラリー量の定量化は、所望のフラグメント範囲の事前サイズ選択なしに大量の再増幅DNAライブラリーに対して行った場合に不正確となる可能性がある。結果として、大量のDNA定量化及びプール中の様々なサンプルの均等化は、シークエンサーサイズ範囲内のフラグメントの数がライブラリー内のDNAフラグメントの分布の不正確さにより確率的に変動することから、1サンプル当たりの実際のリード数の比較的大きな変動を生じる可能性がある(このため、完全に均等化された総量のDNAライブラリーであっても、シークエンシングされるフラグメントの数の顕著な変動を生じる)。
均等化を改善するためにシークエンサーリード長内のDNAライブラリーの量を増大する
(再増幅前のサイズ選択による)
全DNAライブラリーに対するシークエンサーリード長範囲内のDNAライブラリーの量の割合を増大するために、一次WGA産物サイズ分布を変更する必要がある場合、DNAライブラリーサイズ選択を含む初期DNAライブラリー選択を行うことが有利である場合がある。
(選択的再増幅による)
代替的又は付加的には、所望の範囲のDNAライブラリーフラグメントの選択的増幅に好都合な条件下で再増幅反応を行うことが有利である場合がある。
ポリメラーゼの選択又は伸長サイクルの短縮による選択的再増幅
より短いフラグメントに好都合な反応条件は、より短いフラグメントを選択的に増幅するポリメラーゼを用いた再増幅PCR反応、又はより短い伸長段階によって長いフラグメントが完全長まで複製するのを防ぎ、不完全ライブラリーフラグメントを生成する初期PCRサイクルを含み得る。不完全ライブラリーフラグメントは再増幅プライマー(複数の場合もある)の3’セクションに逆相補的な3’末端部分を欠くため、このフラグメントは該再増幅プライマー(複数の場合もある)を用いた更なる複製工程から除外され、不完全フラグメントの指数関数的増幅が妨げられ、より長い一次WGA DNAライブラリーフラグメントに由来する線形の(サイクルによる)数の不完全増幅フラグメントしか生成することができない(consenting)。
表1によるワークフローの例
Wf1)は、Ion Torrent PGM(例えば、サンプルバーコードを必要としない単一サンプルを処理する314チップ)でのWGAライブラリーのLPWGSに適用することができる。2つのプライマーを用いた再増幅により、バーコードを用いない2つのシークエンシングアダプターの導入が可能となる。
Wf2)は、一次WGAの元のインプットサンプルが均一なタイプの非増幅材料、例えば同じ処理を受けた(例えば、新鮮な又は固定された)、非アポトーシス性の単一細胞に由来する場合に、Ion Torrent PGM又はIllumina MiSeqでの多数の一次WGAサンプルのLPWGSに適用することができる。このため、一次WGA収率が全体的におおむね同じであることから、定量化は必要とされない。バーコード化され、シークエンサーに適合したライブラリーをプールした後、サイズ選択して、シークエンサーリード長内の適切なサイズを有するフラグメントを単離し、精製し、シークエンシングする。サイズ選択をゲルによって行う場合、その後に精製が行われる。サイズ選択を例えば両側SPRIビーズ精製によって行う場合、得られるアウトプットは既に精製されており、更なる精製工程は必要でない。
Wf3)は、一次WGAの元のインプットサンプルが、例えば異なる元のDNAの品質を有する(一部が非アポトーシス性、一部がアポトーシス性であり、不均一なゲノム完全性指数を有する;Polzer et al. EMBO Mol Med. 2014 Oct 30;6(11):1371-86を参照されたい)、異なる処理を受けた(例えば、一部が新鮮で一部が固定された)、一部が単一細胞、一部が細胞プールの均一でないタイプの非増幅材料に由来する、Ion Torrent PGM又はIllumina MiSeqでの多数の一次WGAサンプルのLPWGSに適用することができる。このため、一次WGA収率がサンプル間で顕著に異なる可能性があることから、定量化が必要である。Wf2に対し、定量化が行われる。定量化の前に、例えば残留プライマー及びdNTP又はプライマー二量体が除去され、定量化に影響を及ぼさないことから、定量化工程の信頼性を高めるために精製することが有利である。
バーコード化され、シークエンサーに適合したライブラリーをプールした後、サイズ選択して、シークエンサーリード長内の適切なサイズを有するフラグメントを単離し、精製し、シークエンシングする。サイズ選択をゲルによって行う場合、その後に精製が行われる。サイズ選択を例えば両側SPRIビーズ精製によって行う場合、得られるアウトプットは既に精製されており、更なる精製工程は必要でない。
Wf4)は、Oxford Nanoporeシークエンシングのための超並列シークエンシングライブラリーの調製に適用することができる。Nanoporeはより長いリード長に対応することができるため、サイズ選択が不必要である可能性があり、ライブラリー内の実質的に全てのフラグメント長に対してシークエンシングを行うことができる。
Wf5)は、Oxford Nanoporeシークエンシングのための多重超並列シークエンシングライブラリーの調製に適用することができる。wf4に対し、再増幅プライマーは、より多くのサンプルを同じランで多重化するためのサンプルバーコードを更に含む。Nanoporeはより長いリード長に対応することができるため、サイズ選択が不必要である可能性がある。
Wf6)は、特殊アダプターの付加を必要としないOxford Nanoporeシークエンサーのための多重超並列シークエンシングライブラリーの調製に適用することができる。wf5に対し、再増幅プライマーはシークエンシングアダプターを含まず、より多くのサンプルを同じランで多重化するためのサンプルバーコードのみを含む。Nanoporeはより長いリード長に対応することができるため、サイズ選択が不必要である可能性がある。
Wf7)は、より短いリード長のシステム、例えば200bpケミストリー(chemistry:化学的性質)のシークエンシングを用いるIonProtonでの異なるDNA品質を有する不均一な処理後の均一でないサンプルから得られるDRS−WGA DNAライブラリーのシークエンシングのための多重超並列ライブラリーの調製に適用することができる。200bp前後の一次WGA DNAライブラリーの量は、一次WGA DNAライブラリー中の全DNAと比較して非常に少ないため、シークエンシングリード長外の全て又は実質的に全てのpWGAlibフラグメントを排除し、200bp前後のpWGAlibフラグメントを富化するサイズ選択を行うことが有利である場合がある。
次いで、IonProtonシステムに適合するバーコード及びシークエンシングアダプターを含む再増幅プライマーを用いて再増幅を行う。再増幅産物をこのように精製し、各サンプルについて定量化し、各サンプルバーコードについてリード数を均等化するために種々のサンプルの種々のアリコートを共にプールした後、シークエンシングしてLPWGSを行う。
上述のようなワークフローに含まれる工程の種々の組合せが、本発明の範囲から逸脱することなく可能であり、これが本明細書に開示されるような特別なプライマーによる一次WGA DNAライブラリーの再増幅に左右されることが当業者には明らかである。
WGA産物からの超並列シークエンシングライブラリー調製
本発明の第1の実施形態において、
a. 既知の5’WGA配列セクション(5SS)と、中央WGA配列セクション(MSS)と、該既知の5’WGA配列セクションに逆相補的な既知の3’WGA配列セクション(3SS)とを含むフラグメントを含む一次WGA DNAライブラリー(pWGAlib)を準備する工程であって、上記既知の5’WGA配列セクション(5SS)がWGAライブラリー汎用配列アダプターを含み、上記中央WGA配列セクション(MSS)が、WGA前の元の非増幅DNAのDNA配列に対応する挿入セクション(IS)を少なくとも含み、上記中央WGA配列が任意に、隣接5’中間セクション(F5)及び/又は隣接3’中間セクション(F3)を更に含む、工程と、
b. 上記一次WGA DNAライブラリーを、少なくとも1つの第1プライマー(1PR)及び少なくとも1つの第2プライマー(2PR)を用いて再増幅する工程と、
を含み、
上記少なくとも1つの第1プライマー(1PR)が、第1プライマー5’セクション(1PR5S)と第1プライマー3’セクション(1PR3S)とを含み、該第1プライマー5’セクション(1PR5S)が、少なくとも1つの第1シークエンシングアダプター(1PR5SA)と、該少なくとも1つの第1シークエンシングアダプター(1PR5SA)の3’位かつ上記第1プライマー3’セクション(1PR3S)の5’位に少なくとも1つの第1シークエンシングバーコード(1PR5BC)とを含み、上記第1プライマー3’セクション(1PR3S)が、上記既知の5’配列セクション(5SS)又は上記既知の3’配列セクション(3SS)のいずれかにハイブリダイズし、
上記少なくとも1つの第2プライマー(2PR)が、第2プライマー5’セクション(2PR5S)と第2プライマー3’セクション(2PR3S)とを含み、該第2プライマー5’セクション(2PR5S)が、上記少なくとも1つの第1シークエンシングアダプター(1PR5SA)とは異なる少なくとも1つの第2シークエンシングアダプター(2PR5SA)を含み、上記第2プライマー3’セクション(2PR3S)が、上記既知の5’配列セクション(5SS)又は上記既知の3’配列セクション(3SS)のいずれかにハイブリダイズする、方法が提供される。
既知の5’配列セクション(5SS)は、WGAライブラリー汎用配列アダプターからなるのが好ましい。一例として、DRS−WGA(Menarini Silicon BiosystemsのAmpli1(商標) WGAキット等)及びMALBAC(Yikon Genomics)が、本発明による方法のインプットについての要望により上記既知の5’配列セクションに逆相補的な既知の3’配列セクションを有するpWGAlibを生じる。
したがって、WGAライブラリー汎用配列アダプターは、好ましくはDRS−WGAライブラリー汎用配列アダプター(例えば、配列番号282)又はMALBACライブラリー汎用配列アダプター(例えば、配列番号283)、より好ましくはDRS−WGAライブラリー汎用配列アダプターである。
上記第2プライマー(2PR)が、上記少なくとも1つの第2シークエンシングアダプター(2PR5SA)の3’位かつ上記第2プライマー3’セクション(2PR3S)の5’位に少なくとも1つの第2シークエンシングバーコード(2PR5BC)を更に含むことが好ましい。
シークエンシングバーコードの存在から、ローパス全ゲノムシークエンシングの方法であって、
c. 複数のバーコード化された超並列シークエンシングライブラリーを準備すると共に、種々のシークエンシングバーコード(BC)を用いて得られるサンプルをプールする工程と、
d. プールしたライブラリーをシークエンシングする工程と、
を含む、方法が、本発明の一実施形態に従って実行される。
種々のシークエンシングバーコード(BC)を用いたサンプルをプールする工程が、
e.上記バーコード化された超並列シークエンシングライブラリーの各々におけるDNAを定量化する工程と、
f.バーコード化された超並列シークエンシングライブラリーの量を正規化する工程と、
を更に含む。
種々のシークエンシングバーコード(BC)を用いたサンプルをプールする工程が、少なくとも1つの選択塩基対範囲内に含まれるDNAフラグメントを選択する工程を更に含む。このような選択塩基対範囲は、下流のシークエンシングプラットホームの選択を考慮して異なる値を中心とする。例えばIlluminaシークエンシングプラットホームについては、塩基対範囲は650bp、好ましくは400bpを中心とする。他のシークエンシングプラットホーム、例えばIon Torrentについては、塩基対範囲は400bp、好ましくは200bp、より好ましくは150bp又は100bp又は50bpを中心とする。
本発明の更なる一実施形態によると、上記で言及されたローパス全ゲノムシークエンシングの方法は、第1シークエンシングアダプターと第2シークエンシングアダプターとの両方を含むDNAフラグメントを選択する工程を更に含む。
上記第1シークエンシングアダプターと上記第2シークエンシングアダプターとを含むDNAフラグメントを選択する工程は、超並列シークエンシングライブラリーを、例えば官能化常磁性ビーズに含まれる(consisting in)/それを含む少なくとも1つの固相に接触させることによって行うのが好ましい。本発明の方法の一実施形態において、上記常磁性ビーズがストレプトアビジンコーティングで官能化される。
本発明によるローパス全ゲノムシークエンシングの一方法において、少なくとも1つの第1プライマー(1PR)及び少なくとも1つの第2プライマー(2PR)の一方が5’末端でビオチン化され、選択フラグメントが上記ビーズから非ビオチン化ssDNAフラグメントを溶出させることで得られる。
図4から分かるように、上記の場合、再増幅されたWGA dsDNAライブラリーは、1)非ビオチン化dsDNAフラグメント、一方の鎖がビオチン化されたdsDNAフラグメント、及び両方の鎖がビオチン化されたdsDNAフラグメントを含む。本発明の方法は、
g.再増幅されたWGA dsDNAライブラリーを上記官能化常磁性ビーズと共に設計条件下でインキュベートすることにより、ビオチンと該官能化常磁性ビーズに割り当てられたストレプトアビジンとの間に共有結合を生じさせる工程と、
h.結合していない非ビオチン化dsDNAフラグメントを洗い流す工程と、
i.上記官能化常磁性ビーズから残留ssDNAフラグメントを溶出させる工程と、
の更なる工程を含む。
本発明はまた、超並列シークエンシングライブラリー調製キットであって、
第1プライマー5’セクション(1PR5S)と第1プライマー3’セクション(1PR3S)とを含み、該第1プライマー5’セクション(1PR5S)が、少なくとも1つの第1シークエンシングアダプター(1PR5SA)と、該少なくとも1つの第1シークエンシングアダプター(1PR5SA)の3’位かつ上記第1プライマー3’セクション(1PR3S)の5’位に少なくとも1つの第1シークエンシングバーコード(1PR5BC)とを含み、上記第1プライマー3’セクション(1PR3S)が、一次WGA DNAライブラリー(pWGAlib)のフラグメントのWGAライブラリー汎用配列アダプターを含む既知の5’配列セクション(5SS)又は該既知の5’配列セクションに逆相補的な既知の3’配列セクション(3SS)のいずれかにハイブリダイズし、上記フラグメントが上記既知の5’配列セクション(5SS)の3’かつ上記既知の3’配列セクション(3SS)の5’の中央配列セクション(MSS)を更に含む、1つの第1プライマー(1PR)、
第2プライマー5’セクション(2PR5S)と第2 3’セクション(2PR3S)とを含み、該第2プライマー5’セクション(2PR5S)が、上記少なくとも1つの第1シークエンシングアダプター(1PR5SA)とは異なる少なくとも1つの第2シークエンシングアダプター(2PR5SA)を含み、上記第2 3’セクションが、上記フラグメントの上記既知の5’配列セクション(5SS)又は上記既知の3’配列セクション(3SS)のいずれかにハイブリダイズする、1つの第2プライマー(2PR)、
を少なくとも含む、超並列シークエンシングライブラリー調製キットに関する。
特に、超並列シークエンシングライブラリー調製キットは、
a) 配列番号97(表2)のプライマー及び配列番号1〜配列番号96(表2)からなる群から選択される1つ以上のプライマー、又は、
b)配列番号194(表2)のプライマー及び配列番号98〜配列番号193(表2)からなる群から選択される1つ以上のプライマー、又は、
c)配列番号195〜配列番号202(表4)からなる群から選択される少なくとも1つのプライマー及び配列番号203〜配列番号214(表4)からなる群から選択される少なくとも1つのプライマー、又は、
d)配列番号215〜配列番号222(表6)からなる群から選択される少なくとも1つのプライマー及び配列番号223〜配列番号234(表6)からなる群から選択される少なくとも1つのプライマー、又は、
e)配列番号235〜配列番号242(表7)からなる群から選択される少なくとも1つのプライマー及び配列番号243〜配列番号254(表7)からなる群から選択される少なくとも1つのプライマー、又は、
f)配列番号259〜配列番号266(表8)からなる群から選択される少なくとも1つのプライマー及び配列番号267〜配列番号278(表8)からなる群から選択される少なくとも1つのプライマー、
を含む。
本発明の一実施形態によると、超並列シークエンシングライブラリー調製キットは、
最適なシングルリードシークエンシングプロセスが行われるように設計された、
配列番号235〜配列番号242(表7)からなる群から選択される少なくとも1つのプライマー、及び配列番号243〜配列番号254(表7)からなる群から選択される少なくとも1つのプライマー、配列番号255のカスタムシークエンシングプライマー、及び配列番号256若しくは配列番号258のプライマー、又は、
配列番号259〜配列番号266(表8)からなる群から選択される少なくとも1つのプライマー、配列番号267〜配列番号278(表8)からなる群から選択される少なくとも1つのプライマー、並びに配列番号279及び配列番号280のプライマー、
を含む。
上記キットは、選択シークエンシングプラットホームにおいて最適なペアエンドシークエンシングプロセスが行われるように設計された、配列番号257(表7)及び配列番号281(表8)から選択されるプライマーを更に含むことができる。
超並列シークエンシングライブラリー調製キットは、熱安定性DNAポリメラーゼを更に含むことが好ましい。
本発明はまた、最後に、ゲノムワイドコピー数プロファイリングの方法であって、
a. 上記のシークエンシングライブラリー調製キットを用いて開発されたDNAライブラリーをシークエンシングする工程と、
b. ゲノムの異なる領域にわたるシークエンシングリード密度を分析する工程と、
c. 上記ゲノムの領域についてのコピー数の値を、該領域におけるリード数を参照ゲノムの同じ領域において予測されたリード数に対して比較することによって決定する工程と、
を含む、方法に関する。
単一CTCからのローパス全ゲノムシークエンシング
LPWGSによるCNAプロファイリングは、aCGHが十分に明らかな結果を生じることができない場合がある、より低いゲノム完全性指数により寛容である。実際に、aCGHプローブはゲノムの定位置に設計される。これらの位置がDNAの架橋のために確率的に増幅することができない場合、対応するプローブはハイブリダイゼーション後に適切な量のシグナルを生成せず、試験DNAと参照DNAとのシグナル比にノイズピクセル(noisy pixel)を生じる。
これに対し、LPWGSを用いる場合、フラグメントはサイズ選択にのみ基づく。或る特定のフラグメントが例えばDNAの架橋又はアポトーシスによって誘導される切断のために確率的に増幅しない場合、依然として同じローパスビンに分けられる近接ゲノム領域における増幅に適した同じサイズの付加的なフラグメントが存在し得る。したがって、例えば図6〜図9に明らかに示されるように、シグナル対ノイズはライブラリーのゲノム完全性指数により柔軟である(resilient)。
DRS−WGAからの超並列シークエンシングライブラリー調製
サイズ選択は、シークエンシングライブラリーの塩基対サイズとしての実質的に同じ長さ(アダプター挿入の正味)のDRS−WGAフラグメントで構成される領域内のゲノムのサブサンプリングを意味する。
それにもかかわらず、これらのサブサンプリングは、コピー数多型コーリングの標準アルゴリズムを用いる場合であってもコピー数プロファイルの品質に影響を及ぼさないことが驚くべきことに見出された。
TTAA確定的制限部位を有するDRS−WGAを(Ampi1(商標) WGAキットとして)選択するのが有利である。このように、より短いフラグメントはゲノムの低GC含量領域でより高密度であり、フラグメント密度はより高いGC含量と負に相関する。
微量のデジタル処理で選別されたFFPE細胞からのローパス全ゲノムシークエンシング
DEPArrayシステム(Bolognesi et al.)を用いてFFPEからデジタル処理で選別された数百個の腫瘍細胞から出発して、DRS−WGAライブラリーを生成した。図6に示されるように、このライブラリーを用いて本発明によるIon/PGMのための超並列シークエンシングライブラリーを生成した。超並列ライブラリーを0.05未満の平均深度でシークエンシングした。
実施例1:
DRS−WGA後のIon Torrent PGMでのLPWGSのプロトコル
1)確定的制限部位全ゲノム増幅(DRS−WGA)
単一細胞DNAを、Ampli1(商標) WGAキット(Silicon Biosystems)を用い、製造業者の使用説明書に従って増幅した。
Ampli1(商標) WGAキットは、1つの単一細胞から得られたDNAからの全ゲノム増幅をもたらすように設計される。細胞溶解に続いて、DNAを制限酵素、好ましくはMseIで消化し、汎用アダプター配列をDNAフラグメントにライゲートする。増幅は、ゲノム全体に分布する200bp〜1000bp長の範囲サイズの生成される全てのフラグメントについて、単一特異的PCRプライマーによって媒介される。
2)WGA産物の再増幅
5μLのWGAにより増幅されたDNAを、5μLのNuclease−Free Waterを添加することによって希釈し、結合していないオリゴヌクレオチド、並びに過剰なヌクレオチド、塩及び酵素を除去するためにAgencourt AMPure XPシステム(Beckman Coulter)を用いて精製する。
ビーズを用いたDNA精製を以下のプロトコルに従って行った。18μLのビーズ(サンプル体積の1.8倍)を各サンプルに添加した。ビーズ及び反応産物を短時間のボルテックスにより混合した後、遠沈して液滴を収集した。次いで、混合反応物(mixed reactions)をRTで15分間、マグネット外で(off-magnet)インキュベートした後、DynaMag−96 Sideマグネット(Life Technologies)に移し、5分間静置した。上清を廃棄し、ビーズを150μLの新たに作製した80%EtOHで洗浄した。2回目のEtOH洗浄の後、ビーズをマグネット上で5分〜10分間乾燥させた。次いで、乾燥ビーズをマグネット外で15μLのLow TEバッファーに再懸濁し、10分間インキュベートし、続いてマグネット上で5分間のインキュベーションを行った。12マイクロリットルの溶出物を別のチューブに移し、続いて25ngのWGAにより精製されたDNAを含有する10μLアリコートを調製するために、Qubit(商標) 2.0 FluorometerでのdsDNA HSアッセイによって定量化した。
バーコード化再増幅を、Ampli1(商標) PCRキット(Menarini Silicon Biosystems)を用いて50μl容量で行った。各PCR反応物は、以下のもので構成される:5μlのAmpli1(商標) PCR Reaction Buffer(10倍)、1μLの配列番号1〜配列番号96の1つの25μMプライマー:
(1)(5’−CCATCTCATCCCTGCGTGTCTCCGACTCAG[BC−]AGTGGGATTCCTGCTGTCAGT−3’)
(ここで、[BC]=バーコード配列である)、1μLの配列番号97の25μMプライマー:
(2)(5’−CCTCTCTATGGGCAGTCGGTGATAGTGGGATTCCTGCTGTCAGT−3’)
、1.75μlのAmpli1(商標) PCR dNTP(10mM)、1.25μlのBSA、0.5 Ampli1(商標) PCR Taq Polymerase、37.5μlのAmpli1(商標)水及び25ngのWGAにより精製されたDNA。
Applied Biosystems(商標)の2720 Thermal Cyclerは以下のように設定した:95℃で4分間、95℃で30秒間、60℃で30秒間、72℃で2分間の1サイクル、95℃で30秒間、60℃で30秒間、72℃で2分間の10サイクル(20秒/サイクルで延長)、及び72℃で7分間の最終伸長。
図3に再増幅プロセスを概略的に示す。
バーコード化再増幅したWGA産物を1.8倍(90μl)のAMPure XPビーズを用いて精製し、上記の工程に従って35μlのLow TEバッファー中で溶出させた。
3)サイズ選択
Ion Torrentアダプターが設けられたフラグメントライブラリーに対応するバーコード化再増幅したWGA産物を、2100 Bioanalyzer(商標)(Agilent)でAgilentのDNA 7500キットによって適格にし、Qubit(商標) dsDNA HSアッセイキットを用いて定量化し、最終プールを得た。
異なるA−LIB−BC−Xアダプターを有する個々の7つのライブラリーを同じ量で合わせることによって等モルプールを作製し、42μLの最終容量で34ng/μLの濃度の最終プールを得た。プールの濃度をQubit(商標)法によって確認した。
Size Select 2%プレキャストアガロースゲル(Invitrogen)と組み合わせたE−Gel(商標) SizeSelect(商標)システムを、製造業者の使用説明書に従って対象のフラグメントのサイズ選択に使用した。
20μLの最終プールをE−gelの2つのレーンにロードし、サイズ標準(50bp DNA Ladder、Invitrogen)を用いて300bp〜400bpのセクション範囲をゲルから選択した。
サイズ選択後に、1.8倍(90μl)のAMPure XPビーズを用いてクリーンアップを行った。最終ライブラリーを上記の工程に従って30μlのLow TEバッファー中で溶出させ、2100 Bioanalyzer高感度チップ(Agilent Technologies)を用いて評価した。
4)Ion Torrent PGMシークエンシング
鋳型調製を、Ion PGM(商標) Hi−Q OT2キット−400bpのユーザーガイドに従って行った。
簡潔に述べると、ライブラリーフラグメントをIon Sphere Particle(ISP)上でエマルジョンPCRによりクローン的に増幅した後、鋳型陽性ISPについて富化した。PGMエマルジョンPCR反応を、Ion PGM(商標) Hi−Q OT2キット(Life Technologies)を用いて行い、エマルジョン及び増幅物(amplifications)を、Ion OneTouch Instrument(Life Technologies)を用いて生成した。回収後に、増幅されたライブラリーフラグメントを含有するISPをストレプトアビジンコーティング磁性ビーズに選択的に結合させ、空のISPを洗浄工程により除去し、鋳型陽性ISPが収集されるようにライブラリー鎖を変性させることによって富化を行った。
記載の富化工程は、Life Technologies ESシステム(Life Technologies)を用いて達成した。
Ion 318v2(商標)チップを、「Simplified Ion PGMTM Chip loading with the Ion PGMTM weighted chip bucket」プロトコル説明書(MAN0007517)に従ってロードした。
全てのサンプルを、Ion PGM(商標) Hi−Q(商標) Sequencingキット(Life Technologies)を用い、520フローランフォーマットを設定するIon Personal Genome Machine(PGM)(Life Technologies)によって処理した。
最後に、シークエンシングされたフラグメントを、それらの独自バーコードに基づいて特定のサンプルに割り当てた。
実施例2:
DRS−WGA後のIon Torrent ProtonでのLPWGSのプロトコル
1. 確定的制限部位全ゲノム増幅(DRS−WGA)
単一細胞DNAを、先の実施例に詳述されるようにAmpli1(商標) WGAキット(Menarini Silicon Biosystems)を用い、製造業者の使用説明書に従って増幅した。
2. 二本鎖DNA合成
5μLのWGAにより増幅されたDNAを、Ampli1(商標) ReAmp/dsキットを用い、製造プロトコルに従って二本鎖DNA(dsDNA)へと変換した。このプロセスにより、一本鎖DNA(ssDNA)分子のdsDNA分子への変換が確実となる。
3. dsDNA産物の精製
6μLのdsDNA合成産物を、44μLのNuclease−Free Waterを添加することで希釈し、結合していないオリゴヌクレオチド、並びに過剰なヌクレオチド、塩及び酵素を除去するためにAgencourt AMPure XPビーズ(Beckman Coulter)によって精製した。ビーズを用いたDNA精製を以下のプロトコルに従って行った。75μL(比率:サンプル体積の1.5倍)のAgencourt AMPure XPビーズを各50μlのサンプルに添加し、ボルテックスによって混合した。次いで、混合反応物を室温(RT)で15分間、マグネット外でインキュベートした後、溶液が透明になり、ペレットが形成されるまで(約5分間)マグネットプレート上に置いた。次いで、上清を除去し、ペレットをかき乱すことなく廃棄し(およそ5μlが残り得る)、チューブをマグネットプレート上に置いたままでビーズを150μLの新たに作製した70%EtOHで2回洗浄した。残留エタノール溶液を全てチューブの底から除去した後、ビーズペレットを短時間風乾した。22μLの10mM Tris Ultrapure(pH8.0)及び0.1mM EDTA(Low TE)バッファーを添加し、混合反応物を室温で2分間、マグネットプレート外でインキュベートし、続いてマグネットプレート上で5分間のインキュベーションを行った。20μLの溶出物を新たなチューブに移した。
そうでなければ、平均サイズ前後の一様な分布のフラグメントを生じさせるために代替工程3(下記)を用いた。
代替工程3)dsDNA産物の二重精製
SPRI選択は、次世代シークエンシングのためのフラグメントライブラリー調製の核酸サイズ選択を促進し、単純化するSPRIベースのケミストリーである。この工程を工程3に代替して行うことができた。6μLのdsDNA合成産物を、44μLのNuclease−Free Waterを添加することで希釈し、結合していないオリゴヌクレオチド、並びに過剰なヌクレオチド、塩及び酵素を除去し、平均サイズ前後のフラグメントの一様な分布を生じさせるためにSPRI選択ビーズ(Beckman Coulter)によって精製した。SPRIによるDNA精製を以下のプロトコルに従って行った。37.5μL(比率:サンプル体積の0.75倍)のSPRI選択ビーズを各50μlのサンプルに添加し、ボルテックスによって混合した。次いで、混合反応物をRTで15分間、マグネット外でインキュベートした後、溶液が透明になり、ペレットが形成されるまで(約5分間)マグネットプレート上に置いた。次いで、上清を回収し、新たなチューブに移した。2回目の精製を、37.5μLのSPRI選択ビーズを上清に添加して行い、ボルテックスによって混合した。次いで、混合反応物をRTで15分間、マグネット外でインキュベートした後、溶液が透明になり、ペレットが形成されるまで(約5分間)マグネットプレート上に置いた。次いで、上清を除去し、ペレットをかき乱すことなく廃棄し(およそ5μlが残り得る)、チューブをマグネットプレート上に置いたままでビーズを150μLの新たに作製した80%EtOHで2回洗浄した。残留エタノール溶液を全てチューブの底から除去した後、ビーズペレットを短時間風乾した。22μLのLow TEバッファーを添加し、混合反応物を室温で2分間、マグネット外でインキュベートし、続いてマグネットプレート上で5分間のインキュベーションを行った。20μLの溶出物を新たなチューブに移した。
4. バーコード化再増幅
バーコード化再増幅を、Ampli1(商標) PCRキット(Menarini Silicon Biosystems)を用いて50μl容量で行った。各PCR反応物は、以下のもので構成される:5μlのAmpli1(商標) PCR Reaction Buffer(10倍)、1μlの配列番号1〜配列番号96の1つの25μMプライマー:
(1)(5’−CCATCTCATCCCTGCGTGTCTCCGACTCAG[BC]AGTGGGATTCCTGCTGTCAGT−3’)
(ここで、[BC]=バーコード配列である)、1μlの配列番号97の25μMプライマー:
(2)(5’−CCTCTCTATGGGCAGTCGGTGATAGTGGGATTCCTGCTGTCAGT−3’)
、1.75μlのAmpli1(商標) PCR dNTP(10mM)、1.25μlのBSA、0.5 Ampli1(商標) PCR Taq Polymerase(FAST start)、37.5μlのAmpli1(商標)水及び2μlのdsにより精製されたDNA。これらは、Ion Torrent PGMに使用される同じプライマーであり、上記に示したIon TorrentライブラリーのためのNGS再増幅プライマー(PGM/ProtonのためのDRS WGA)の対応する表に報告する。
Applied Biosystems(商標)の2720 Thermal Cyclerは以下のように設定した:95℃で4分間、95℃で30秒間、60℃で30秒間、72℃で15秒間の11サイクル、その後72℃で30秒間の最終伸長。
5. バーコード化再増幅したdsDNA産物の精製
バーコード化再増幅したdsDNA産物を、1.5倍(75μl)の比率のAMPure XPビーズを用い、上記の工程3に従って精製し、35μlのLow TEバッファー中で溶出させた。溶出物を新たなチューブに移し、続いてQubit(商標) 2.0 FluorometerでdsDNA HSアッセイによって定量化し、最終等モルサンプルプールを得た。異なるA−LIB−BC−Xアダプターを有する各々のライブラリーを同じ量で合わせることによって等モルプールを作製し、42μLの最終容量で34ng/μLの濃度の最終プールを得た。
6. サイズ選択
Size Select 2%プレキャストアガロースゲル(Invitrogen)と組み合わせたE−Gel(商標) SizeSelect(商標)システムを、製造業者の使用説明書に従って対象のフラグメントのサイズ選択に使用した。
20μLの最終プールをE−gelの2つのレーンにロードし、サイズ標準(50bp DNA Ladder、Invitrogen)を用いて300bp〜400bpのセクション範囲をゲルから選択した。サイズ選択後に、上記の工程3に従って1.8倍(90μl)のAMPure XPビーズを用いてクリーンアップを行った。最終ライブラリーを30μlのLow TEバッファー中で溶出させた。
7. Ion Torrent Protonシークエンシング
精製工程後の等モルプールを、2100 Bioanalyzer(商標)(Agilent)でAgilentのDNA高感度キットによって適格にし、Qubit(商標) dsDNA HSアッセイキットを用いて定量化した。最後に、等モルプールを100pMの最終濃度まで希釈した。
鋳型調製をIon PI(商標) Hi−Q(商標) Chefのユーザーガイドに従って行った。Ion Chef(商標)システムは、Ion Proton(商標) Sequencerで用いられる自動化ハイスループット鋳型調製及びチップのローディングをもたらす。Ion Proton(商標) Sequencerにより、Ion Proton(商標) Hi−Q(商標)シークエンシングキット(Life Technologies)を用いて、Ion PI(商標)チップ上にロードされたライブラリーの自動化ハイスループットシークエンシングを行う。最後に、シークエンシングしたフラグメントを、それらの独自バーコードに基づいて特定のサンプルに割り当てた。
実施例3:
Illumina MiSeqでのローパス全ゲノムシークエンシングのプロトコル
プロトコル1
確定的制限部位全ゲノム増幅(DRS−WGA):
単一細胞DNAを、Ampli1(商標) WGAキット(Silicon Biosystems)を用い、製造業者の使用説明書に従って増幅した。5μLのWGAにより増幅されたDNAを、5μLのNuclease−Free Waterの添加によって希釈し、Agencourt AMPure XPシステムを用いて精製した(比率1.8倍)。DNAを12.5μL中で溶出させ、Qubit(商標) 2.0 FluorometerでのdsDNA HSアッセイによって定量化した。
バーコード化再増幅
バーコード化再増幅を、50μlの容量でAmpli1(商標) PCRキット(Menarini Silicon Biosystems)を用いて図4に概略的に示されるように行った。各PCR反応物は、以下のもので構成される:5μlのAmpli1(商標) PCR Reaction Buffer(10倍)、1μlの配列番号195〜配列番号202の1つのプライマー(25μM)、1μlの配列番号203〜配列番号214の1つのプライマー(25μM)、1.75μlのAmpli1(商標) PCR dNTP(10mM)、1.25μlのBSA、0.5 Ampli1(商標) PCR Taq Polymerase、25ngのWGAにより精製されたDNA、及び50μlの最終容量に達するまでのAmpli1(商標)水。
Applied Biosystems(商標)の2720 Thermal Cyclerは以下のように設定した:95℃で4分間、95℃で30秒間、60℃で30秒間、72℃で2分間の1サイクル、95℃で30秒間、60℃で30秒間、72℃で2分間の10サイクル(20秒/サイクルで延長)、及び72℃で7分間の最終伸長。
次いで、バーコード化再増幅したWGA産物(リストSEQ_IDs_ILL_PR1から選び出したIlluminaシークエンシングアダプター配列を含有する)を、2100 Bioanalyzer(商標)でAgilentのDNA 7500キットによって適格にし、Qubit(商標) 2.0 Fluorometerによって定量化した。
サイズ選択
次いで、ライブラリーを等モル濃度で合わせ、濃度28.6ng/μL及び最終容量100μLで得られるプールを、SPRIビーズを用いた二重精製によってサイズ選択した。簡潔に述べると、SPRIビーズをPCRグレード水で2倍希釈した。160μLの希釈SPRIビーズを100μlのプールに添加した。インキュベーション後に、25μLの上清を新たなバイアルに移し、30μLの希釈SPRIビーズを添加した。DNAを20μLのlow TE中で溶出させた。フラグメントサイズを2100 Bioanalyzer高感度チップ(Agilent Technologies)によって検証し、ライブラリー定量化を、Kapa Library定量化キットを用いるqPCRによって行った。
MiSeqシークエンシング
4nMのサイズ選択プールを、NaOH(0.1Nに等しいNaOH最終濃度)で5分間変性させた。次いで、変性させたサンプルをHT1で希釈して、1mM NaOH中の20pM変性ライブラリーを得た。570μLの20pM変性ライブラリー及び30μlの20pM変性PhiX対照をMiSeq(Illumina)でロードした。
150塩基のシングルエンドリードを、IlluminaのMiSeqのv3ケミストリーを用いて生成した。
配列番号リスト_第1プライマー_[ILLUMINA/DRS−WGA]プロトコル1
以下の表に、Illuminaプラットホームに適合するプライマーDRS−WGAの構造(5’→3’方向の配列、5’及び3’を除外)を示す。
以下の表に最終プライマー配列を報告する。
配列番号リスト_第1プライマー_[ILLUMINA/MALBAC]プロトコル1
以下の表に、Illuminaプラットホームに適合するプライマーMALBAC−WGAの構造(5’→3’方向の配列、5’及び3’を除外)を示す。
以下の表に最終プライマー配列を報告する。
プロトコル1の制限
Illuminaプロトコル1によって得られるライブラリーは、考え得る全てのP5/P7アダプター組合せ対を有する二本鎖pWGA libである。
フローセル内で同種のシークエンシングアダプター(P5/P5rc、P7rc/P7)を有するフラグメントによってもハイブリダイゼーションが生じたため、クラスター密度及び/又はクラスターの品質は、Illuminaプロトコル2の場合と比較して僅かに低くなる可能性がある。
プロトコル2
本発明による第2のプロトコルを例として提示する。本プロトコルは、pWGAlibから両方のシークエンシングアダプター(P5/P7)を含むフラグメントのみを選択し、同種のシークエンシングアダプター(P5/P5rc、P7rc/P7)を有するフラグメントを切り捨てることによってIlluminaフローセルにおけるクラスターの品質を高めるために有利であり得る。
図4に概略的に示されるプロトコル2(Illumina/DRS_WGA)のワークフローの説明
WGAにより増幅されたDNA産物は全て、長さが異なり、特定のタグ:個々のssDNA分子の各々の5’末端のLIB配列及び3’末端の相補的なLIB配列(図面中に青色で示される)を有する分子によって構成される。
本発明によると、両方の逆相補LIB配列がNGS Re−Amp(再増幅)プライマーの標的である。
それぞれ図面中で緑色−黄色−青色及び赤色−ピンク色−青色の2つのタイプのプライマー:LPb_DI_D50X(配列番号235〜配列番号242の範囲のプライマー)及びビオチン化プライマーLPb_DI_D70X(配列番号243〜配列番号254の範囲のプライマー)を設計した。
予測されるように、両方のタイプのプライマーがLIB配列及び相補的なLIB配列を結合することができ、実際のところ、図面に示されるように3つのタイプのアンプリコンがNGS Re−Ampプロセスによって生じる。
本発明による本プロトコルは、5’末端にビオチンタグを得るLPb_DI_D70X(図面中にP7rcアダプターとして示される)によってもたらされる。この特定のタグを用いて、図面に示されるようにビオチンタグを有しない唯一のフラグメント:
5’−P5−i5−LIB−挿入断片−LIB相補鎖−i7−P7−3’
をストレプトアビジンビーズにより選択する。
望ましい構成のssDNAを得るために(単純化のためにリードプライマーセクションを省略すると、望ましいssDNAは、5’−P5−i5−nnnnn−i7−P7−3’である)、プライマーは以下のようなものであるとする:
(1PR)5’−P5−i5−LIB−3’、及び、
(2PR)ビオチン−5’−P7rc−i7rc−LIB−3’(以下ではビオチンを説明の単純化のために省略するが、P7rcで始まるフラグメントの全てに、5’末端にのみ存在することが明らかである)。
再増幅により、以下のものが得られる:
開始:(WGA ssDNAフラグメントは全て、5’−LIB−nnn−LIBrc−3’として形成される)
伸長サイクルn=1):
1.5 5’−P5−i5−LIB−nnn−LIBrc−3’、
1.7 5’−P7rc−i7rc−LIB−nnn−LIBrc−3’
2^nのフラグメント(0%シークエンシング可能)
サイクルn=2):
2.5.5 5’−P5−i5−LIB−nnn−LIBrc−i5rc−P5rc−3’
2.5.7 5’−P7rc−i7rc−LIB−nnn−LIBrc−i5rc−P5rc−3’
2.7.5 5’−P5−i5−LIB−nnn−LIBrc−i7−P7−3’
2.7.7 5’−P7rc−i7rc−LIB−nnn−LIBrc−i7−P7−3’
2^n=4のフラグメント(25%のシークエンシング可能なフラグメント)
サイクルn=3):
2.5.5.5 5’−P5−i5−LIB−nnn−LIBrc−i5rc−P5rc−3’=2.5.5
2.5.5.7 5’−P7rc−i7rc−LIB−nnn−LIBrc−i5rc−P5rc−3’=2.5.7
2.5.7.5 5’−P5−i5−LIB−nnn−LIBrc−i7−P7−3’=2.7.5
2.5.7.7 5’−P7rc−i7rc−LIB−nnn−LIBrc−i7−P7−3’=2.7.7
2.7.5.5 5’−P5−i5−LIB−nnn−LIBrc−i5rc−P5rc−3’=2.5.5
2.7.5.7 5’−P7rc−i7rc−LIB−nnn−LIBrc−i5rc−P5rc−3’=2.5.7
2.7.7.5 5’−P5−i5−LIB−nnn−LIBrc−i7−P7−3’=2.7.5
2.7.7.7 5’−P7rc−i7rc−LIB−nnn−LIBrc−i7−P7−3’=2.7.7
2^n=8のフラグメント(25%のシークエンシング可能なフラグメント) シークエンシング可能なフラグメント=2^n/4=2^n/2^2=2^(n−2)
サイクルm) ... 2^(m−2) シークエンシング可能
最終的に、以下の4つのタイプのフラグメントが指数関数的増幅後に形成される。
2.5.5 5’−P5−i5−LIB−nnn−LIBrc−i5rc−P5rc−3’(→全てのビオチン化フラグメントが常磁性ビーズ上に保持される一方で第1の液体除去時に洗い流されるか、又は洗い流されない場合には、フローセル内の唯一の結合部位に連結するが、架橋増幅が起こらないことからシークエンシングクラスターを生成しない)
2.5.7 ビオチン−5’−P7rc−i7rc−LIB−nnn−LIBrc−i5rc−P5rc−3’(→ストレプトアビジンコーティングビーズによって除去される)
2.7.5 5’−P5−i5−LIB−nnn−LIBrc−i7−P7−3’(→シークエンシング可能)
2.7.7 ビオチン−5’−P7rc−i7rc−LIB−nnn−LIBrc−i7−P7−3’(→ストレプトアビジンコーティングビーズによって除去される)。
実施例4:
1.確定的制限部位全ゲノム増幅(DRS−WGA)
単一細胞DNAを、Ampli1(商標) WGAキット(Silicon Biosystems)を用い、製造業者の使用説明書に従って増幅した。
2.WGA産物の再増幅
5μLのWGAにより増幅されたDNAを、5μLのNuclease−Free Waterを添加することによって希釈し、結合していないオリゴ、並びに過剰なヌクレオチド、塩及び酵素を除去するためにAgencourt AMPure XPシステム(Beckman Coulter)を用いて精製する。
ビーズを用いたDNA精製を以下のプロトコルに従って行った。18μLのビーズ(サンプル体積の1.8倍)を各サンプルに添加した。ビーズ及び反応産物を短時間のボルテックスにより混合した後、遠沈して液滴を収集した。次いで、混合反応物を室温で15分間、マグネット外でインキュベートした後、DynaMag−96 Sideマグネット(Life Technologies)に移し、5分間静置した。上清を廃棄し、ビーズを150μLの新たに作製した80%EtOHで洗浄した。2回目のEtOH洗浄の後、ビーズをマグネット上で5分〜10分間乾燥させた。次いで、乾燥ビーズをマグネット外で15μLのLow TEバッファーに再懸濁し、10分間インキュベートし、続いてマグネット上で5分間のインキュベーションを行った。12マイクロリットルの溶出物を別のチューブに移し、続いて25ngのWGAにより精製されたDNAを含有する10μLアリコートを調製するために、Qubit(商標) 2.0 FluorometerでのdsDNA HSアッセイによって定量化した。
バーコード化再増幅を、50μlの容量でAmpli1(商標) PCRキット(Silicon Biosystems)を用いて行った。各PCR反応物は、以下のように構成される:5μlのAmpli1(商標) PCR Reaction Buffer(10倍)、1μLの配列番号235〜配列番号242の1つの25μMプライマー:
(3)5’AATGATACGGCGACCACCGAGATCTACAC[i5]GCTCTCCGTAGTGGGATTCCTGCTGTCAGTTAA3’)
1μLの配列番号243〜配列番号254の1つの25μMプライマー:
(4)(5’/Biosg/CAAGCAGAAGACGGCATACGAGAT[i7]GCTCACCGAAGTGGGATTCCTGCTGTCAGTTAA3’)
、1.75μlのAmpli1(商標) PCR dNTP(10mM)、1.25μlのBSA、0.5 Ampli1(商標) PCR Taq Polymerase、及び25ngのWGAにより精製されたDNA、及び37.5μlのAmpli1(商標)水。
Applied Biosystems(商標)の2720 Thermal Cyclerは以下のように設定した:95℃で4分間、95℃で30秒間、60℃で30秒間、72℃で2分間の1サイクル、95℃で30秒間、60℃で30秒間、72℃で2分間の10サイクル(20秒/サイクルで延長)、及び72℃で7分間の最終伸長。
3)サイズ選択
Illuminaアダプターが設けられたフラグメントライブラリーに対応するバーコード化再増幅したWGA産物を、2100 Bioanalyzer(商標)(Agilent)でAgilentのDNA 7500キットによって適格にし、Qubit(商標) dsDNA HSアッセイキットを用いて定量化し、最終プールを得た。
異なるLPb_DIデュアルインデックスアダプターを有する個々のライブラリーを同じ量で合わせることによって等モルプールを作製し、50μLの最終容量で35ng/μLの濃度の最終プールを得た。プールの濃度をQubit(商標)法によって確認した。
200bp〜1Kbのフラグメントセクション範囲を、SPRIビーズシステム(Beckman Coulter)をそれぞれR:0.47×及びL:0.85×の比率で用いる二重精製によって選択した。大きなDNAフラグメントを除去するために82μLの希釈SPRI(42μLのSPRIビーズ+42μLのPCRグレード水)を添加し、小さなDNAフラグメントを除去するために34.2μLの不希釈SPRIビーズを上清に添加した。
最終ライブラリーを50μlのLow TEバッファー中で溶出させ、2100 Bioanalyzer高感度チップ(Agilent Technologies)を用いて評価した。
4)異種P5/P7アダプター一本鎖ライブラリー選択
フラグメント選択を、P5/P7アダプターを含有する非ビオチン化DNAのみを解離させるためにDynabeads(商標) MyOne(商標) Streptavidin C1システムを用いて行ったが、これは熱又はNaOHをそれぞれ用いて得ることができた。2つの方法を下記に説明する。
1.5mlチューブ中の20μLのDynabeads(商標) MyOne(商標) Streptavidin C1を、1倍B&W溶液(10mM Tris−HCl(pH7.5);1mM EDTA;2M NaCl)で2回洗浄した。
50μLの分画プールライブラリーをDynabeads(商標) MyOne(商標) Streptavidin C1ビーズに添加し、完全に混合するために5分間隔でピペッティングしながら15分間インキュベートした。DNAでコーティングされたDynabeads(商標)を50μLの1倍SSC(0.15M NaCl、0.015Mクエン酸ナトリウム)で2回洗浄し、新たな50μLの1倍SSCを用いてビーズを再懸濁した。
95℃で5分間のインキュベーションの後、チューブをマグネットプレートに1分間置き、50μLの上清を新たなチューブに移し、氷上で5分間インキュベートした。
この点で、上清はP5/P7アダプターを有する非ビオチン化DNA鎖ライブラリーを含有する。
洗浄がより厳しくなるように、ストレプトアビジン選択手順を再度繰り返した。
熱を用いる代わりに、DNAでコーティングされたDynabeads(商標)の洗浄を、20μlの新たに調製した0.15M NaOHを用いて再懸濁することによって行うことができる。
室温で10分間のインキュベーションの後、チューブをマグネット内に1分〜2分間静置し、上清を新たなチューブに移した。
上清は非ビオチン化DNA鎖を含有する。一本鎖ライブラリーを、2.2μLの10倍TE(pH7.5)及び1.3μLの1.25M酢酸を添加することによって中和した。
Qubit(商標) ssDNAアッセイキットによって定量化された最終ライブラリー濃度は、25μMに相当する5ng/μLであった。
5)MiSeqシークエンシング
4nMの最終プールをNaOH(0.1Nに等しいNaOH最終濃度)で5分間変性させた。次いで、変性サンプルをHT1で希釈し、1mM NaOH中の20pM変性ライブラリーを得た。570μLの20pM変性ライブラリー及び30μlの20pM変性PhiX対照をMiSeq(Illumina)でロードした。
150塩基のシングルエンドリードを、IlluminaのMiSeqのv3ケミストリーを用い、標準リード1プライマー及び標準プライマーインデックス1を、それぞれ600μLの配列番号255のプライマー(カスタムリード1プライマー)及び600μLの配列番号256又は配列番号258のプライマー(カスタムプライマーインデックス1a(i7)及び1b(i7))と交換して生成した。
配列番号リスト_第1プライマー_[ILLUMINA_DRS_WGA]プロトコル2
以下の表に、Illuminaプロトコル2の最終プライマー配列を報告する。
配列番号リスト_第1プライマー_[ILLUMINA/MALBAC]プロトコル2
以下の表に、出発物質がWGA−MALBACライブラリーに由来する場合に適合したIllumina最終プライマー配列を報告する。
本発明によると、図4に示されるように両方のLIB逆相補鎖がNGS Re−Amp(再増幅)プライマーの標的である。さらに、リードが同じヌクレオチドで開始せず、シークエンシング性能が影響を受けるか、又はマトリックス、フェージング及びプレフェージングの算出のためのより多様なクラスターセットを確実にする高濃度スパイクイン(spike-in)の使用が回避される可能性があることから、ライブラリーの複雑さを増大するために、カスタムリード1シークエンシングプライマー(配列番号255)を設計した。図4に示されるように、カスタムリード1シークエンシングプライマー(配列番号255)はLIB配列を含有し、LIB逆相補配列に相補的である。
さらに、NGS Re−Amp(再増幅)産物は、インデックス1を読み取るためにIlluminaシステムに用いられるカノニカル配列を有さず、この理由から、インデックスi7の正確な読取りを可能にするためにカスタムシークエンシングプライマーインデックス1(i7)(配列番号256又は配列番号258)を使用する必要がある。カスタムシークエンシングプライマーインデックス1が逆相補的LIB配列を含有することに注目すべきである。
手順PGM/Proton及びIlluminaプロトコル1/2ワークフローを含む上記の全ての実施例は、上記表に挙げられる適合するプライマーMALBAC(配列番号98〜配列番号194及び配列番号215〜配列番号234及び配列番号259〜配列番号281)を用いて行うことができる。
データ分析
シークエンシングリードを、BWA MEMアルゴリズム(Li H. and Durbin R., 2010)を用いてhg19ヒト参照ゲノムにアラインメントした。PCR重複、二次/補助/QC不合格(not-passing-QC)アラインメント及びmultimapperリードを、Picard MarkDuplicates(http://broadinstitute.github.io/picard/)及びsamtools(Li H. et al, 2009)を用いてフィルタリングした。カバレッジ分析を、BEDTools(Quinlan A. et al, 2010)を用いて行った。
Control−FREEC(Boeva V. et al., 2011)アルゴリズムを用いて、対照サンプルを用いずにコピー数コールを得た。リード数をGC含量によって補正し、マッピング可能性(mappability)(uniqMatchオプション)及びウィンドウサイズを、変動係数=0.06を用いるソフトウェアによって決定した。倍数性を2に設定し、汚染調整は用いなかった。
CNVプロファイルのプロットを、図6〜図9に示されるようにカスタムpythonスクリプトを用いて得た。
本発明をAmpli1 WGAの方法のみを参照して説明したが、記載の技法は、当業者には明白であるように、変更すべき点を変更して、自己相補的5’及び3’領域を有するライブラリーを含む任意の他の種類のWGA(例えば、MALBAC)にも明らかに適用される。

Claims (27)

  1. 超並列シークエンシングライブラリーを生成する方法であって、
    a. 既知の5’配列セクション(5SS)と、中央配列セクション(MSS)と、該既知の5’配列セクションに逆相補的な既知の3’配列セクション(3SS)とを含むフラグメントを含む一次WGA DNAライブラリー(pWGAlib)を準備する工程であって、前記既知の5’配列セクション(5SS)がWGAライブラリー汎用配列アダプターを含み、前記中央配列セクション(MSS)が、WGA前の元の非増幅DNAのDNA配列に対応する挿入セクション(IS)を少なくとも含み、前記中央配列セクションが任意に、隣接5’中間セクション(F5)及び/又は隣接3’中間セクション(F3)を更に含む、工程と、
    b. 前記一次WGA DNAライブラリーを、少なくとも1つの第1プライマー(1PR)及び少なくとも1つの第2プライマー(2PR)を用いて再増幅する工程と、
    を含み、
    前記少なくとも1つの第1プライマー(1PR)が、第1プライマー5’セクション(1PR5S)と第1プライマー3’セクション(1PR3S)とを含み、該第1プライマー5’セクション(1PR5S)が、少なくとも1つの第1シークエンシングアダプター(1PR5SA)と、該少なくとも1つの第1シークエンシングアダプター(1PR5SA)の3’位かつ前記第1プライマー3’セクション(1PR3S)の5’位に少なくとも1つの第1シークエンシングバーコード(1PR5BC)とを含み、前記第1プライマー3’セクション(1PR3S)が、前記既知の5’配列セクション(5SS)又は前記既知の3’配列セクション(3SS)のいずれかにハイブリダイズし、
    前記少なくとも1つの第2プライマー(2PR)が、第2プライマー5’セクション(2PR5S)と第2プライマー3’セクション(2PR3S)とを含み、該第2プライマー5’セクション(2PR5S)が、前記少なくとも1つの第1シークエンシングアダプター(1PR5SA)とは異なる少なくとも1つの第2シークエンシングアダプター(2PR5SA)を含み、前記第2プライマー3’セクション(2PR3S)が、前記既知の5’配列セクション(5SS)又は前記既知の3’配列セクション(3SS)のいずれかにハイブリダイズする、方法。
  2. 前記第2プライマー(2PR)が、前記少なくとも1つの第2シークエンシングアダプター(2PR5SA)の3’位かつ前記第2プライマー3’セクション(2PR3S)の5’位に少なくとも1つの第2シークエンシングバーコード(2PR5BC)を更に含む、請求項1に記載の方法。
  3. 前記WGAライブラリー汎用配列アダプターが、DRS−WGAライブラリー汎用配列アダプター又はMALBACライブラリー汎用配列アダプターである、請求項1又は2に記載の方法。
  4. 前記WGAライブラリー汎用配列アダプターがDRS−WGAライブラリー汎用配列アダプターである、請求項3に記載の方法。
  5. 前記DRS−WGAライブラリー汎用配列アダプターが配列番号282であり、前記MALBACライブラリー汎用配列アダプターが配列番号283(MALBAC)である、請求項3又は4に記載の方法。
  6. ローパス全ゲノムシークエンシングの方法であって、
    c. 請求項1〜5のいずれか一項に記載の方法によって得られる複数のバーコード化された超並列シークエンシングライブラリーを準備すると共に、種々のシークエンシングバーコード(BC)を用いて得られるサンプルをプールする工程と、
    d. プールしたライブラリーをシークエンシングする工程と、
    を含む、方法。
  7. 種々のシークエンシングバーコード(BC)を用いたサンプルをプールする工程が、
    e)前記バーコード化された超並列シークエンシングライブラリーの各々におけるDNAを定量化する工程と、
    f)バーコード化された超並列シークエンシングライブラリーの量を正規化する工程と、
    を更に含む、請求項6に記載のローパス全ゲノムシークエンシングの方法。
  8. 種々のシークエンシングバーコード(BC)を用いたサンプルをプールする工程が、少なくとも1つの選択塩基対範囲を有するDNAフラグメントを選択する工程を更に含む、請求項7に記載のローパス全ゲノムシークエンシングの方法。
  9. 前記塩基対範囲が650bpを中心とする、請求項8に記載のローパス全ゲノムシークエンシングの方法。
  10. 前記塩基対範囲が400bpを中心とする、請求項8に記載のローパス全ゲノムシークエンシングの方法。
  11. 前記塩基対範囲が200bpを中心とする、請求項8に記載のローパス全ゲノムシークエンシングの方法。
  12. 前記塩基対範囲が150bpを中心とする、請求項8に記載のローパス全ゲノムシークエンシングの方法。
  13. 前記塩基対範囲が100bpを中心とする、請求項8に記載のローパス全ゲノムシークエンシングの方法。
  14. 前記塩基対範囲が50bpを中心とする、請求項8に記載のローパス全ゲノムシークエンシングの方法。
  15. 第1シークエンシングアダプターと第2シークエンシングアダプターとを含むDNAフラグメントを選択する工程を更に含む、請求項8〜14のいずれか一項に記載のローパス全ゲノムシークエンシングの方法。
  16. 前記第1シークエンシングアダプター及び前記第2シークエンシングアダプターを含むDNAフラグメントを選択する工程を、前記超並列シークエンシングライブラリーを少なくとも1つの固相に接触させることによって行う、請求項15に記載のローパス全ゲノムシークエンシングの方法。
  17. 前記少なくとも1つの固相が官能化常磁性ビーズを含む、請求項16に記載のローパス全ゲノムシークエンシングの方法。
  18. 前記常磁性ビーズがストレプトアビジンコーティングで官能化される、請求項17に記載のローパス全ゲノムシークエンシングの方法。
  19. 少なくとも1つの第1プライマー(1PR)及び少なくとも1つの第2プライマー(2PR)の一方が5’末端でビオチン化され、選択フラグメントが前記ビーズから非ビオチン化ssDNAフラグメントを溶出させることで得られる、請求項18に記載のローパス全ゲノムシークエンシングの方法。
  20. 前記少なくとも1つの第2プライマーが5’末端でビオチン化される、請求項19に記載のローパス全ゲノムシークエンシングの方法。
  21. g)再増幅されたWGA dsDNAライブラリーを前記官能化常磁性ビーズと共に設計条件下でインキュベートすることにより、ビオチンと該官能化常磁性ビーズに割り当てられたストレプトアビジンとの間に共有結合を生じさせる工程と、
    h)結合していない非ビオチン化dsDNAフラグメントを洗い流す工程と、
    i)前記官能化常磁性ビーズから残留ssDNAフラグメントを溶出させる工程と、
    の更なる工程を更に含む、請求項18〜20のいずれか一項に記載のローパス全ゲノムシークエンシングの方法。
  22. 超並列シークエンシングライブラリー調製キットであって、
    第1プライマー5’セクション(1PR5S)と第1プライマー3’セクション(1PR3S)とを含み、該第1プライマー5’セクション(1PR5S)が、少なくとも1つの第1シークエンシングアダプター(1PR5SA)と、該少なくとも1つの第1シークエンシングアダプター(1PR5SA)の3’位かつ前記第1プライマー3’セクション(1PR3S)の5’位に少なくとも1つの第1シークエンシングバーコード(1PR5BC)とを含み、前記第1プライマー3’セクション(1PR3S)が、一次WGA DNAライブラリー(pWGAlib)のフラグメントのWGAライブラリー汎用配列アダプターを含む既知の5’配列セクション(5SS)又は該既知の5’配列セクションに逆相補的な既知の3’配列セクション(3SS)のいずれかにハイブリダイズし、前記フラグメントが前記既知の5’配列セクション(5SS)の3’かつ前記既知の3’配列セクション(3SS)の5’の中央配列セクション(MSS)を更に含む、1つの第1プライマー(1PR)、
    第2プライマー5’セクション(2PR5S)と第2 3’セクション(2PR3S)とを含み、該第2プライマー5’セクション(2PR5S)が、前記少なくとも1つの第1シークエンシングアダプター(1PR5SA)とは異なる少なくとも1つの第2シークエンシングアダプター(2PR5SA)を含み、前記第2 3’セクションが、前記フラグメントの前記既知の5’配列セクション(5SS)又は前記既知の3’配列セクション(3SS)のいずれかにハイブリダイズする、1つの第2プライマー(2PR)、
    を少なくとも含む、超並列シークエンシングライブラリー調製キット。
  23. 超並列シークエンシングライブラリー調製キットであって、
    a. 配列番号97(表2)のプライマー及び配列番号1〜配列番号96(表2)からなる群から選択される1つ以上のプライマー、又は、
    b)配列番号194(表2)のプライマー及び配列番号98〜配列番号193(表2)からなる群から選択される1つ以上のプライマー、又は、
    c)配列番号195〜配列番号202(表4)からなる群から選択される少なくとも1つのプライマー及び配列番号203〜配列番号214(表4)からなる群から選択される少なくとも1つのプライマー、又は、
    d)配列番号215〜配列番号222(表6)からなる群から選択される少なくとも1つのプライマー及び配列番号223〜配列番号234(表6)からなる群から選択される少なくとも1つのプライマー、又は、
    e)配列番号235〜配列番号242(表7)からなる群から選択される少なくとも1つのプライマー及び配列番号243〜配列番号254(表7)からなる群から選択される少なくとも1つのプライマー、又は、
    f)配列番号259〜配列番号266(表8)からなる群から選択される少なくとも1つのプライマー及び配列番号267〜配列番号278(表8)からなる群から選択される少なくとも1つのプライマー、
    を含む、超並列シークエンシングライブラリー調製キット。
  24. 超並列シークエンシングライブラリー調製キットであって、
    最適なシングルリードシークエンシングプロセスが行われるように設計された、
    配列番号235〜配列番号242(表7)からなる群から選択される少なくとも1つのプライマー、配列番号243〜配列番号254(表7)からなる群から選択される少なくとも1つのプライマー、配列番号255のカスタムシークエンシングプライマー、及び配列番号256若しくは配列番号258のプライマー、又は、
    配列番号259〜配列番号266(表8)からなる群から選択される少なくとも1つのプライマー、配列番号267〜配列番号278(表8)からなる群から選択される少なくとも1つのプライマー、並びに配列番号279及び配列番号280のプライマー、
    を含む、超並列シークエンシングライブラリー調製キット。
  25. 選択シークエンシングプラットホームにおいて最適なペアエンドシークエンシングプロセスが行われるように設計された、配列番号257(表7)及び配列番号281(表8)から選択されるプライマーを更に含む、請求項24に記載の超並列シークエンシングライブラリー調製キット。
  26. 熱安定性DNAポリメラーゼを更に含む、請求項22〜25のいずれか一項に記載の超並列シークエンシングライブラリー調製キット。
  27. ゲノムワイドコピー数プロファイリングの方法であって、
    a. 請求項22又は23に記載のシークエンシングライブラリー調製キットを用いて開発されたDNAライブラリーをシークエンシングする工程と、
    b. ゲノムの異なる領域にわたるシークエンシングリード密度を分析する工程と、
    c. 前記ゲノムの領域についてのコピー数の値を、該領域におけるリード数を参照ゲノムの同じ領域において予測されたリード数に対して比較することによって決定する工程と、
    を含む、方法。
JP2018553460A 2016-04-15 2017-04-14 超並列シークエンシングのためのdnaライブラリーを生成する方法及びキット Active JP6899844B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
ITUA2016A002640A ITUA20162640A1 (it) 2016-04-15 2016-04-15 Metodo e kit per la generazione di librerie di dna per sequenziamento massivo parallelo
IT102016000039154 2016-04-15
PCT/EP2017/059075 WO2017178655A1 (en) 2016-04-15 2017-04-14 Method and kit for the generation of dna libraries for massively parallel sequencing

Publications (2)

Publication Number Publication Date
JP2019514360A true JP2019514360A (ja) 2019-06-06
JP6899844B2 JP6899844B2 (ja) 2021-07-07

Family

ID=56853703

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018553460A Active JP6899844B2 (ja) 2016-04-15 2017-04-14 超並列シークエンシングのためのdnaライブラリーを生成する方法及びキット

Country Status (9)

Country Link
US (2) US20200299680A1 (ja)
EP (1) EP3443115A1 (ja)
JP (1) JP6899844B2 (ja)
KR (1) KR102354422B1 (ja)
CN (1) CN109477245A (ja)
CA (1) CA3019714A1 (ja)
IT (1) ITUA20162640A1 (ja)
SG (2) SG10202010198VA (ja)
WO (1) WO2017178655A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3431611A1 (en) * 2017-07-21 2019-01-23 Menarini Silicon Biosystems S.p.A. Improved method and kit for the generation of dna libraries for massively parallel sequencing
CN108866155B (zh) * 2018-06-11 2022-07-26 中国农业科学院深圳农业基因组研究所 一种下一代测序文库的制备方法
IT201900013335A1 (it) * 2019-07-30 2021-01-30 Menarini Silicon Biosystems Spa Metodo per analizzare la perdita di eterozigosi (loh) a seguito di amplificazione totale del genoma basata su un sito di restrizione deterministico (drs-wga)
CN112662748B (zh) * 2019-10-15 2023-02-17 骏实生物科技(上海)有限公司 用于单细胞dna测序文库质量鉴定的引物对组合、方法和试剂盒
IT201900024159A1 (it) 2019-12-16 2021-06-16 Menarini Silicon Biosystems Spa Metodo e kit per amplificazione totale del genoma e analisi di molecole bersaglio in un campione biologico
CN111005074A (zh) * 2019-12-19 2020-04-14 江西海普洛斯医学检验实验室有限公司 一种基于illumina测序平台的DNA文库构建试剂盒、文库构建方法和应用
CN112592968B (zh) * 2020-12-27 2022-07-26 苏州科诺医学检验实验室有限公司 高通量测序用分子标签接头及其合成方法与应用
IT202100024101A1 (it) 2021-09-20 2023-03-20 Menarini Silicon Biosystems Spa Metodo per analizzare il grado di similarita' di almeno due campioni utilizzando amplificazione deterministica dell'intero genoma mediante siti di restrizione (drs-wga)
CN117721191A (zh) * 2024-02-07 2024-03-19 深圳赛陆医疗科技有限公司 基因测序方法、测序装置、可读存储介质和基因测序系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002526087A (ja) * 1998-09-18 2002-08-20 マイクロメット アーゲー 単一の細胞のdna増幅
JP2007523627A (ja) * 2003-01-29 2007-08-23 454 コーポレーション 核酸を増幅および配列決定する方法
WO2015148219A1 (en) * 2014-03-28 2015-10-01 Clarient Diagnostic Services, Inc. Accurate detection of rare genetic variants in next generation sequencing
WO2015171656A1 (en) * 2014-05-06 2015-11-12 Baylor College Of Medicine Methods of linearly amplifying whole genome of a single cell
JP2015533296A (ja) * 2012-11-05 2015-11-24 ルビコン ゲノミクス インコーポレイテッド バーコード化する核酸

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040209299A1 (en) * 2003-03-07 2004-10-21 Rubicon Genomics, Inc. In vitro DNA immortalization and whole genome amplification using libraries generated from randomly fragmented DNA
US8206913B1 (en) * 2003-03-07 2012-06-26 Rubicon Genomics, Inc. Amplification and analysis of whole genome and whole transcriptome libraries generated by a DNA polymerization process
WO2007073165A1 (en) 2005-12-22 2007-06-28 Keygene N.V. Method for high-throughput aflp-based polymorphism detection
AU2009226248B8 (en) * 2008-03-17 2014-10-02 Stichting Genetwister Ip Expression-linked gene discovery
CN102559856B (zh) * 2010-12-22 2014-03-12 深圳华大基因科技服务有限公司 去除测序文库中的载体片段的方法
KR102393608B1 (ko) * 2012-09-04 2022-05-03 가던트 헬쓰, 인크. 희귀 돌연변이 및 카피수 변이를 검출하기 위한 시스템 및 방법
ITTO20120962A1 (it) 2012-10-31 2014-05-01 Silicon Biosystems Spa Metodo e kit per rivelare una sequenza di dna bersaglio wild-type e/o mutata
CN105189780A (zh) * 2012-12-03 2015-12-23 以琳生物药物有限公司 核酸制备和分析的组合物和方法
US20140256571A1 (en) * 2013-03-06 2014-09-11 Life Technologies Corporation Systems and Methods for Determining Copy Number Variation
PT2881739T (pt) 2013-12-04 2017-08-11 Fraunhofer Ges Forschung Método e estojo para determinação da integridade do genoma e/ou da qualidade de uma biblioteca de sequências de adn obtida através de amplificação determinística do genoma inteiro por locais de restrição

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002526087A (ja) * 1998-09-18 2002-08-20 マイクロメット アーゲー 単一の細胞のdna増幅
JP2007523627A (ja) * 2003-01-29 2007-08-23 454 コーポレーション 核酸を増幅および配列決定する方法
JP2015533296A (ja) * 2012-11-05 2015-11-24 ルビコン ゲノミクス インコーポレイテッド バーコード化する核酸
WO2015148219A1 (en) * 2014-03-28 2015-10-01 Clarient Diagnostic Services, Inc. Accurate detection of rare genetic variants in next generation sequencing
WO2015171656A1 (en) * 2014-05-06 2015-11-12 Baylor College Of Medicine Methods of linearly amplifying whole genome of a single cell

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HOU YONG ET AL.: ""Comparison of variations detection between whole−", GIGASCIENCE, vol. 4, JPN6020042552, 2015, pages 1 - 16, ISSN: 0004382054 *
TSAI ISHENG J.ET AL.: ""Summarizing Specific Profiles in Illumina Sequenc", DNA RESEARCH, vol. 21, no. 3, JPN6020042549, 2014, pages 243 - 254, ISSN: 0004382055 *

Also Published As

Publication number Publication date
US20230250419A1 (en) 2023-08-10
SG10202010198VA (en) 2020-11-27
SG11201808675VA (en) 2018-11-29
CN109477245A (zh) 2019-03-15
KR102354422B1 (ko) 2022-01-21
ITUA20162640A1 (it) 2017-10-15
US20200299680A1 (en) 2020-09-24
CA3019714A1 (en) 2017-10-19
KR20190019049A (ko) 2019-02-26
WO2017178655A1 (en) 2017-10-19
EP3443115A1 (en) 2019-02-20
JP6899844B2 (ja) 2021-07-07

Similar Documents

Publication Publication Date Title
JP6899844B2 (ja) 超並列シークエンシングのためのdnaライブラリーを生成する方法及びキット
US20210071171A1 (en) Compositions and methods for targeted nucleic acid sequence enrichment and high efficiency library generation
US10072283B2 (en) Direct capture, amplification and sequencing of target DNA using immobilized primers
US11965157B2 (en) Compositions and methods for library construction and sequence analysis
CN107075513A (zh) 分离的寡核苷酸及其在核酸测序中的用途
JP2023126945A (ja) 超並列シークエンシングのためのdnaライブラリー生成のための改良された方法及びキット
KR20170133270A (ko) 분자 바코딩을 이용한 초병렬 시퀀싱을 위한 라이브러리 제조방법 및 그의 용도
EP3652345A2 (en) Sequencing method for genomic rearrangement detection
EP3565906B1 (en) Quantifying dna sequences
JP2022544779A (ja) ポリヌクレオチド分子の集団を生成するための方法
US20220411861A1 (en) A Multiplex Method of Preparing a Sequencing Library
Dilanyan Open-Source Custom Beads for Single-Cell Transcriptomics
JP2024035110A (ja) 変異核酸の正確な並行定量するための高感度方法
JP2024035109A (ja) 核酸の正確な並行検出及び定量のための方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200129

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200925

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201110

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210210

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: 20210608

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210615

R150 Certificate of patent or registration of utility model

Ref document number: 6899844

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150