JP7638309B2 - 低減した増幅バイアスによるハイスループット単一細胞シークエンシング - Google Patents
低減した増幅バイアスによるハイスループット単一細胞シークエンシング Download PDFInfo
- Publication number
- JP7638309B2 JP7638309B2 JP2023025703A JP2023025703A JP7638309B2 JP 7638309 B2 JP7638309 B2 JP 7638309B2 JP 2023025703 A JP2023025703 A JP 2023025703A JP 2023025703 A JP2023025703 A JP 2023025703A JP 7638309 B2 JP7638309 B2 JP 7638309B2
- Authority
- JP
- Japan
- Prior art keywords
- cells
- nucleic acid
- nuclei
- indexed
- amplification
- 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
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
- C12N15/1003—Extracting or separating nucleic acids from biological samples, e.g. pure separation or isolation methods; Conditions, buffers or apparatuses therefor
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
- C12N15/1034—Isolating an individual clone by screening libraries
- C12N15/1065—Preparation or screening of tagged libraries, e.g. tagged microorganisms by STM-mutagenesis, tagged polynucleotides, gene tags
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
- C12N15/1034—Isolating an individual clone by screening libraries
- C12N15/1093—General methods of preparing gene libraries, not provided for in other subgroups
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N9/00—Enzymes; Proenzymes; Compositions thereof; Processes for preparing, activating, inhibiting, separating or purifying enzymes
- C12N9/10—Transferases (2.)
- C12N9/12—Transferases (2.) transferring phosphorus containing groups, e.g. kinases (2.7)
- C12N9/1241—Nucleotidyltransferases (2.7.7)
- C12N9/1247—DNA-directed RNA polymerase (2.7.7.6)
-
- 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
-
- 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/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C40—COMBINATORIAL TECHNOLOGY
- C40B—COMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
- C40B40/00—Libraries per se, e.g. arrays, mixtures
- C40B40/04—Libraries containing only organic compounds
- C40B40/06—Libraries containing nucleotides or polynucleotides, or derivatives thereof
-
- C—CHEMISTRY; METALLURGY
- C40—COMBINATORIAL TECHNOLOGY
- C40B—COMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
- C40B50/00—Methods of creating libraries, e.g. combinatorial synthesis
- C40B50/06—Biochemical methods, e.g. using enzymes or whole viable microorganisms
-
- C—CHEMISTRY; METALLURGY
- C40—COMBINATORIAL TECHNOLOGY
- C40B—COMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
- C40B50/00—Methods of creating libraries, e.g. combinatorial synthesis
- C40B50/08—Liquid phase synthesis, i.e. wherein all library building blocks are in liquid phase or in solution during library creation; Particular methods of cleavage from the liquid support
- C40B50/10—Liquid phase synthesis, i.e. wherein all library building blocks are in liquid phase or in solution during library creation; Particular methods of cleavage from the liquid support involving encoding steps
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B82—NANOTECHNOLOGY
- B82Y—SPECIFIC USES OR APPLICATIONS OF NANOSTRUCTURES; MEASUREMENT OR ANALYSIS OF NANOSTRUCTURES; MANUFACTURE OR TREATMENT OF NANOSTRUCTURES
- B82Y5/00—Nanobiotechnology or nanomedicine, e.g. protein engineering or drug delivery
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biochemistry (AREA)
- Biotechnology (AREA)
- General Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Plant Pathology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Structural Engineering (AREA)
- Nanotechnology (AREA)
- Medical Informatics (AREA)
- Pharmacology & Pharmacy (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Preparation Of Compounds By Using Micro-Organisms (AREA)
- Bipolar Transistors (AREA)
Description
この出願は2018年5月17日に出願の米国特許仮出願第62/673,023号および2019年3月21日に出願の米国特許仮出願第62/821,864号の権益を主張し、その各々は参照により完全に本明細書に組み込まれる。
政府資金提供
現代の単一細胞ゲノムシークエンシング技術は、2つの主要な限界を有する。第1に、ほとんどの方法は個々の細胞を区分することを必要とし、それはスループットを制限する可能性がある。第2に、ほとんどの増幅方法はPCRベースであり、したがって、指数関数的増幅バイアスを被る。第1の問題を解決するために、我々および同僚は、単一細胞の核酸内容を特異的にタグ付けし、それによって各連続した回のインデクシングによりスループットにおける指数関数的増加を可能にするために、数回のスプリット-プール分子バーコード化を実行する、単一細胞コンビナトリアルインデクシング(「sci-」)を開発した。sci-方法は、多数の単一細胞において、クロマチンアクセス性(sci-ATAC-seq)、トランスクリプトーム(sci-RNA-seq)、ゲノム(sci-DNA-seq)、メチローム(sci-MET)、染色体高次構造(sci-Hi-C)をプロファイリングするように首尾よく開発されている(Cao et al., 2017, Science 357:661-667;Cusanovich et al., 2015, Science, 348:910-914;Mulqueen et al., 2018, Nat. Biotechnol. 36:428-431;Ramani et al., 2017, Nat. Methods 14:263-266;Vitak et al., 2017, Nat. Methods 14:302-308)。第2の問題を解決するために、T7ベースの転写による線形増幅は、単一細胞アッセイとの関係で以前に利用された潜在的解法を提供する(Eberwine et al., 1992; Proceedings of the National Academy of Sciences 89:3010-3014;Hashimshony et al., 2012, Cell Rep. 2:666-673;Sos et al., 2016, Genome Biolol., 17:20)。例えば、最近、Chen et al.は、ゲノムを断片化し、インビトロ転写(IVT)のためにT7 RNAプロモーターを同時に挿入するためにTn5トランスポゾンを使用するトランスポゾン挿入による線形増幅を開発した(「LIANTI」)。DNA鋳型から生成されるRNAコピーは、さらなる増幅のために鋳型の役割をすることができない。したがって、全てのコピーは、元のDNA鋳型に直接的に由来する。指数関数的増幅を回避することによって、LIANTIは均一性を維持し、配列のエラーを最小にする。しかし、この方法は各単一細胞からの連続的ライブラリー調製を必要とするので、ロースループットである(Chen et al., 2017, Science 356:189-194)。
本明細書で使用される用語は、特に明記しない限り、関連する分野におけるそれらの通常の意味をとることが理解される。本明細書で使用されるいくつかの用語およびそれらの意味は、本明細書に示される。
本開示は例えば以下を提供する。
[項1]
複数の単一の核または細胞からの核酸を含むシークエンシングライブラリーを調製するための方法であって、
第1の複数のコンパートメントに複数の単離された核または細胞を提供する工程であって、
各コンパートメントは単離された核または細胞のサブセットを含み、
核または細胞は核酸断片を含む、工程と;
線形増幅媒介物を前記細胞または核に導入する工程と;
線形増幅によって前記核酸断片を増幅する工程と;
核または細胞の各サブセットを処理してインデックス付きの核または細胞を生成する工程であって、
前記処理は単離された前記核または細胞に存在する核酸断片に第1のコンパートメント特異的インデックス配列を加えて単離された核または細胞に存在するインデックス付きの核酸をもたらすことを含み、
前記処理は、ライゲーション、プライマー伸長、ハイブリダイゼーション、増幅または転位を含む工程と;
インデックス付きの前記核または細胞を合わせてプールされたインデックス付きの核または細胞を生成し、それによって複数の前記核または細胞からシークエンシングライブラリーを生成する工程と
を含む方法。
[項2]
前記増幅が前記処理の前に起こる、項1に記載の方法。
[項3]
前記処理が前記増幅の前に起こる、項1に記載の方法。
[項4]
複数の単一の核または細胞からの核酸を含むシークエンシングライブラリーを調製するための方法であって、
複数の単離された核または細胞を提供する工程であって、
核または細胞は核酸断片を含む、工程と;
線形増幅媒介物を単離された前記核または細胞に導入する工程と;
単離された前記核または細胞を第1の複数のコンパートメントに分配する工程であって、
各コンパートメントは単離された核または細胞のサブセットを含む工程と;
線形増幅によって前記核酸断片を増幅する工程と;
単離された核または細胞の各サブセットを処理してインデックス付きの核または細胞を生成する工程であって、
前記処理は単離された前記核または細胞に存在する核酸断片に第1のコンパートメント特異的インデックス配列を加えて単離された核または細胞に存在するインデックス付きの核酸をもたらすことを含み、
前記処理は、ライゲーション、プライマー伸長、増幅または転位を含む工程と;
インデックス付きの前記核を合わせてプールされたインデックス付きの核または細胞を生成し、それによって複数の前記核または細胞からシークエンシングライブラリーを生成する工程と
を含む方法。
[項5]
複数の単一の核または細胞からの核酸を含むシークエンシングライブラリーを調製するための方法であって、
複数の単離された核または細胞を第1の複数のコンパートメントに提供する工程であって、
各コンパートメントは単離された核または細胞のサブセットを含み、
核または細胞は核酸断片を含む、工程と;
核または細胞の各サブセットを処理してインデックス付きの核または細胞を生成する工程であって、
前記処理は、単離された前記核または細胞に存在する核酸断片に、(i)第1のコンパートメント特異的インデックス配列を加えて単離された核または細胞に存在するインデックス付きの核酸をもたらすこと、および(ii)線形増幅媒介物によって認識されるヌクレオチド配列を加えることを含み、
前記処理は、ライゲーション、プライマー伸長、ハイブリダイゼーション、増幅または転位を含む工程と;
線形増幅媒介物を前記細胞または核に導入する工程と;
線形増幅によって前記核酸断片を増幅する工程と;
インデックス付きの前記核または細胞を合わせてプールされたインデックス付きの核または細胞を生成し、それによって複数の前記核または細胞からシークエンシングライブラリーを生成する工程と
を含む方法。
[項6]
前記線形増幅媒介物がファージRNAポリメラーゼまたは線形増幅プライマーを含む、項1、4または5のいずれか一項に記載の方法。
[項7]
前記核酸断片がT7プロモーターを含み、前記ファージRNAポリメラーゼがT7 RNAポリメラーゼを含む、項6に記載の方法。
[項8]
前記線形増幅媒介物を導入することが、前記線形増幅媒介物を単離された前記核または細胞に存在する核酸断片に加えることを含む、項6に記載の方法。
[項9]
各コンパートメントの複数の単離された前記核または細胞を所定の条件に曝露させることをさらに含む、項1または5に記載の方法。
[項10]
前記曝露の後に複数の前記細胞から核を単離することをさらに含む、項9に記載の方法。
[項11]
複数の単離された前記核または細胞を所定の条件に曝露させることをさらに含む、項4に記載の方法。
[項12]
単離された前記核の完全性を維持しつつヌクレオソームが枯渇された核を生成する条件に単離された前記核を付すことをさらに含む、項1、4または5のいずれか一項に記載の方法。
[項13]
前記処理が、
各サブセットをトランスポソーム複合体と接触させる工程であって、
各コンパートメントの中の前記トランスポソーム複合体は、他のコンパートメントの中の第1のインデックス配列と異なる第1のインデックス配列を含む工程と;
前記サブセットの中の核酸を複数の核酸に断片化し、前記第1のインデックス配列を前記核酸の少なくとも1つの鎖に組み込んでインデックス付きの前記核酸を含むインデックス付きの前記核または細胞を生成する工程と
を含む、項1、4または5のいずれか一項に記載の方法。
[項14]
前記処理が、
各サブセットを逆転写酵素および単離された前記核の中のRNA分子にアニールするプライマーと接触させる工程であって、各コンパートメントの中の前記プライマーは、他のコンパートメントの中の第1のインデックス配列と異なる第1のインデックス配列を含み、インデックス付きの前記核酸を含むインデックス付きの前記核または細胞を生成する工程
を含む、項1、4または5のいずれか一項に記載の方法。
[項15]
前記接触させることが特異的ヌクレオチド配列にアニールする標的特異的プライマーをさらに含む、項14に記載の方法。
[項16]
前記第1のコンパートメント特異的インデックス配列を加える前記処理が、遍在配列を含むヌクレオチド配列を前記核酸断片に加え、次に前記第1のコンパートメント特異的インデックス配列を前記核酸断片に加える2工程法を含む、項1、4または5のいずれか一項に記載の方法。
[項17]
前記加えることが前記遍在配列を含むトランスポソーム複合体を含む、項16に記載の方法。
[項18]
前記処理が、第1のインデックスを単離された前記核または細胞に存在するDNA核酸に加えること、第1のインデックスを単離された前記核または細胞に存在するRNA核酸に加えること、またはその組合せを含む、項1、4または5のいずれか一項に記載の方法。
[項19]
前記の第1のインデックス配列をRNA核酸に加えることが、
各サブセットを逆転写酵素および単離された前記核または細胞の中のRNA分子にアニールするプライマーと接触させる工程であって、
各コンパートメントの中の前記プライマーは、前記第1のコンパートメント特異的インデックス配列を含み、インデックス付きの前記核酸を含むインデックス付きの前記核または細胞を生成する工程
を含む、項18に記載の方法。
[項20]
前記の第1のインデックス配列をDNA核酸に加えることが、
各サブセットをトランスポソーム複合体と接触させる工程であって、
各コンパートメントの中の前記トランスポソーム複合体は前記第1のコンパートメント特異的インデックス配列を含む工程と;
前記サブセットの中の核酸を複数の核酸に断片化し、前記第1のコンパートメント特異的インデックス配列を前記核酸の少なくとも1つの鎖に組み込んでインデックス付きの前記核酸を含むインデックス付きの前記核または細胞を生成する工程と
を含む、項18に記載の方法。
[項21]
各コンパートメントにおいてDNA核酸に加えられる前記第1のインデックス配列およびRNA核酸に加えられる前記第1のインデックス配列が同一である、項19に記載の方法。
[項22]
各コンパートメントにおいてDNA核酸に加えられる前記第1のインデックス配列およびRNA核酸に加えられる前記第1のインデックス配列が同一でない、項19に記載の方法。
[項23]
前記核酸断片の指数関数的増幅をさらに含み、前記指数関数的増幅は特異的ヌクレオチド配列にアニールする標的特異的プライマーを含む、項1、4または5のいずれか一項に記載の方法。
[項24]
前記合わせることの後に、
プールされたインデックス付きの前記核または細胞のサブセットを第2の複数のコンパートメントに分配する工程と;
第2のコンパートメント特異的インデックス配列をインデックス付きの核酸に導入して二重インデックス付きの核酸を含む二重インデックス付きの核または細胞を生成する工程であって、
前記導入することは、ライゲーション、プライマー伸長、増幅または転位を含む工程と
をさらに含む、項1、4または5のいずれか一項に記載の方法。
[項25]
二重インデックス付きの前記核を合わせてプールされた二重インデックス付きの核または細胞を生成する工程と、
プールされた二重インデックス付きの前記核または細胞のサブセットを第3の複数のコンパートメントに分配する工程と;
第3のコンパートメント特異的インデックス配列をインデックス付きの核酸に導入して三重インデックス付きの核酸を含む三重インデックス付きの核または細胞を生成する工程であって、
前記導入することは、ライゲーション、プライマー伸長、増幅または転位を含む工程と
をさらに含む、項24に記載の方法。
[項26]
メチル化分析のためにインデックス付きの前記核または細胞を処理してメチル化分析のために適する核酸断片を生成することをさらに含む、項1、4または5のいずれか一項に記載の方法。
[項27]
インデックス付きの前記核または細胞を近接性ライゲーションに付してクロマチン高次構造の分析のために適する核酸断片を生成することをさらに含む、項1、4または5のいずれか一項に記載の方法。
[項28]
前記シークエンシングライブラリーの前記核酸断片を増幅してDNAナノボールを生成することをさらに含む、項1、4または5のいずれか一項に記載の方法。
[項29]
前記コンパートメントがウェルまたは液滴を含む、項1~28のいずれか一項に記載の方法。
[項30]
前記第1の複数のコンパートメントの各コンパートメントが50から100,000,000個の核または細胞を含む、項1~29のいずれか一項に記載の方法。
[項31]
前記第2の複数のコンパートメントの各コンパートメントが50から100,000,000個の核または細胞を含む、項1~29のいずれか一項に記載の方法。
[項32]
複数の増幅部位を含む表面を提供する工程であって、
前記増幅部位が遊離の3’末端を有する付着した一本鎖の捕捉オリゴヌクレオチドの少なくとも2つの集団を含む工程と、
複数の増幅部位を生成するのに適する条件の下で増幅部位を含む前記表面をインデックス付きの前記断片と接触させ、各々は複数のインデックスを含む個々の断片からのアンプリコンのクローン集団を含む工程と
をさらに含む、項1~31のいずれか一項に記載の方法。
[項33]
複数の単一細胞からの核酸を含むシークエンシングライブラリーを調製する方法であって、
(a)複数の細胞からの単離された核を提供する工程と;
(b)単離された前記核の完全性を維持しつつヌクレオソームが枯渇された核を生成する化学処理に単離された前記核を付す工程と;
(c)ヌクレオソームが枯渇された前記核のサブセットを第1の複数のコンパートメントに分配し、各サブセットをトランスポソーム複合体と接触させる工程であって、各コンパートメントの中の前記トランスポソーム複合体がトランスポザーゼおよび他のコンパートメントの中の第1のインデックス配列と異なる第1のインデックス配列を含む工程と;
(d)ヌクレオソームが枯渇された核の前記サブセットの中の核酸を複数の核酸断片に断片化し、前記第1のインデックス配列を前記核酸断片の少なくとも1つの鎖に組み込んでインデックス付きの核酸断片を含むインデックス付きの核を生成する工程であって、インデックス付きの前記核酸断片は前記トランスポザーゼに付着したままである工程と;
(d)インデックス付きの前記核を合わせてプールされたインデックス付きの核を生成する工程と;
(e)プールされたインデックス付きの前記核のサブセットを第2の複数のコンパートメントに分配し、各サブセットをヘアピンライゲーション二重鎖と、インデックス付きの核酸断片の片方または両方の末端への前記ヘアピンライゲーション二重鎖のライゲーションに適する条件の下で接触させて二重インデックス付きの核酸断片をもたらす工程であって、前記ヘアピンライゲーション二重鎖が他のコンパートメントの中の第2のインデックス配列と異なる第2のインデックス配列を含む工程と;
(f)二重インデックス付きの前記核を合わせてプールされたインデックス付きの核を生成する工程と;
(g)プールされた二重インデックス付きの前記核のサブセットを第3の複数のコンパートメントに分配する工程と;
(h)二重インデックス付きの前記核を溶解させる工程と;
(i)他のコンパートメントの中の第3のインデックス配列と異なる第3のインデックス配列を含むように二重インデックス付きの前記核断片を処理する工程と;
(j)前記三重インデックス断片を合わせ、それによって複数の前記単一細胞からの全ゲノム核酸を含むシークエンシングライブラリーを生成する工程と
を含む方法。
本明細書で提供される方法は、細胞または複数の細胞から単離した核を提供することを含む。細胞および核は、任意の試料、例えば、任意の生物体(複数可)に、および生物体(複数可)の任意の細胞型または任意の組織に由来することができる。一実施形態では、細胞は、生殖細胞、例えば精子細胞または卵細胞であってよい。一実施形態では、組織は、生殖組織、例えば精巣上体であってよい。一実施形態では、細胞または核は、がんまたは病的な組織に由来することができる。本方法は、細胞を解離させることおよび/または核を単離することをさらに含むことができる。核を細胞から単離するための方法は当業者に公知であり、ルーチンである。核または細胞の数は、少なくとも2つであってよい。上限は、本明細書に記載される方法の他の工程で使用される装置(例えば、マルチウェルプレート)の実際の限界に依存する。使用することができる核または細胞の数は限定するものではなく、数10億を数えることができる。例えば、一実施形態では、核または細胞の数は、100,000,000以下、10,000,000以下、1,000,000,000以下、100,000,000以下、10,000,000以下、1,000,000以下、100,000以下、10,000以下、1,000以下、500以下、または50以下であってよい。1つまたは複数の試料を提供することができる。例えば、試料は、1生物体からの1つの細胞型または組織であってよい。試料を同定し、次に組み合わせるために、本明細書に記載されるインデクシング方法を使用して、複数の試料、例えば、1生物体からの異なる細胞型、2つ以上の生物体からの1つの細胞型もしくは組織、または2つ以上の生物体からの異なる細胞型もしくは組織を第1のインデックスで別々にインデックス付けすることができる。当業者は、一部の実施形態では、各核の中の核酸分子は、生物体の全体の遺伝子相補体(生物体の全ゲノムとも呼ばれる)を表し、イントロンおよびエクソン配列の両方、ならびにプロモーターおよびエンハンサー配列などの非コード調節配列を含むゲノムDNA分子であることを認める。
本明細書で提供される方法は、核、例えばヌクレオソームが枯渇された核または細胞のサブセットを、複数のコンパートメントに分配することを含む。本方法は、単離した核または細胞の集団(本明細書ではプールとも呼ばれる)がサブセットに分割される、複数の分配工程を含むことができる。一般的に、単離した核または細胞のサブセットのプールから複数のコンパートメントへの分配は、単離した核または細胞のサブセットに存在する核酸断片へのインデックスの付加の前に起こる。したがって、本方法は、プールされた単離した核または細胞をとってそれらを分配する少なくとも1つの「スプリットアンドプール」工程を含み、ここで、「スプリットアンドプール」工程の数は、核酸断片に加えられる異なるインデックスの数に依存することができる。インデックスを付けた後、十分な数のインデックスが核酸断片に加えられるまで、必要に応じてサブセットをプールし、サブセットに分割し、インデックスを付け、再びプールすることができる。
一実施形態では、単離した核または細胞の中のDNA核酸、例えば染色体および/またはプラスミドを核酸断片に断片化するために、単離した核または細胞の処理を使用することができる。配列決定をする標的核酸が核または細胞に存在するDNAに由来するとき、処理が一般的に必要である。しかし、一部の実施形態では、RNA分子は断片化する必要がしばしばないので、配列決定をする標的核酸が核または細胞に存在するRNA(例えば、mRNAおよび/または非コードRNA)に由来するとき、処理は適宜である。核または細胞の中の核酸を処理することは、ヌクレオチド配列を処理によって生成される核酸断片の片方または両方の末端に一般的に加え、ヌクレオチド配列は1つまたは複数の遍在配列を含むことができ、および一般的に含む。遍在配列は、例えば、ライゲーション、プライマー伸長または増幅の以降の工程による核酸断片へのインデックスなどの別のヌクレオチド配列の付加のためのプライマーとして使用することができるヌクレオチド配列をアニールする以降の工程における「ランディングパッド」として使用することができる。そのようなプライマーのヌクレオチド配列は、インデックス配列を含んでもよい。核または細胞の中の核酸を処理することは、処理によって生成される核酸断片の片方または両方の末端に1つまたは複数の特異な分子識別子を加えることができる。
タグまたはバーコードとも呼ばれるインデックス配列は、特定の核酸が存在したコンパートメントに特徴的なマーカーとして有益である。したがって、インデックスは、特定のコンパートメントに存在する標的核酸の各々に付着している核酸配列タグであり、その存在は、単離した核または細胞の集団が本方法の特定の段階に存在したコンパートメントを示すか、またはそれを同定するのに使用される。インデックスの核酸断片への付加は、異なるコンパートメントに分配される単離した核または細胞のサブセットで達成される。
本明細書で提供される方法は、核酸断片の線形増幅を含む。ほとんどの増幅方法はPCRベースであり、したがって、指数関数的増幅バイアスを被る。本明細書で使用される線形増幅は指数関数的増幅バイアスを低減または排除することができ、それによってより良い均一性および低減された配列エラーにつながる。全ゲノム増幅を利用する全ての単一細胞ゲノム研究方法では、増幅産物はコンパートメント(例えばウェルまたは液滴)に入れられ、バーコードが直接的または間接的に増幅産物に付着している。このように、1コンパートメントにつき単一の細胞だけが存在し、スループットが制限され、費用が増加する。本発明の特異な態様は、単一のコンパートメントの中で指数関数的増幅バイアスなしで複数の単一細胞ライブラリーを増幅することができることである。単一細胞からのライブラリーは、特異な単一細胞ごとの特異な1つまたは複数のバーコードに基づいて割り当てることができる。
一実施形態では、処理および/またはインデクシング工程の間のヌクレオチドの付加は、断片の固定化およびシークエンシングにおいて有益である遍在配列を加える。別の実施形態では、インデックス付きの核酸断片は、核酸断片の固定化およびシークエンシングにおいて有益な遍在配列を加えるために、さらに処理することができる。当業者は、コンパートメントが液滴である実施形態において、核酸断片を固定化するための配列は適宜であることを認める。一実施形態では、断片の固定化およびシークエンシングにおいて有益な遍在配列の組込みは、同一の遍在アダプター(「ミスマッチアダプター」とも呼ばれ、その一般的フィーチャーは、Gormleyら、米国特許第7,741,463号、およびBignellら、米国特許第8,053,192号に記載される)をインデックス付きの核酸断片の5’および3’末端にライゲーションすることを含む。一実施形態では、遍在アダプターは、アレイの上にインデックス付きの核酸断片を固定化するための配列を含む、シークエンシングに必要な全ての配列を含む。
シークエンシングのために、複数のインデックス付きの断片を調製することができる。例えば、三重インデックス付きの断片のライブラリーが生成される実施形態では、三重インデックス付きの断片は、シークエンシングの前に一般的に固定化および/または増幅によって富化される(図2、ブロック21)。1つまたは複数の供給源からのインデックス付きの断片を基質に付着させるための方法は、当技術分野で公知である。一実施形態では、インデックス付きの断片は、そのインデックス付きの断片に特異性を有する複数の捕捉オリゴヌクレオチドを使用して富化され、捕捉オリゴヌクレオチドは、固体基質の表面に固定化することができる。例えば、捕捉オリゴヌクレオチドは遍在の結合性の対の第1のメンバーを含むことができ、ここで、結合性の対の第2のメンバーは固体基質の表面に固定化される。同様に、固定化された二重インデックス付きの断片を増幅する方法には、限定されずに、架橋増幅および動態学的排除が含まれる。シークエンシングの前に固定化および増幅する方法は、例えば、Bignellら(米国特許第8,053,192号)、Gundersonら(WO2016/130704)、Shenら(米国特許第8,895,249号)およびPipenburgら(米国特許第9,309,502号)に記載される。
表面へのインデックス付きの断片の付着の後、固定化され、増幅されたインデックス付きの断片の配列が決定される。シークエンシングは、任意の適するシークエンシング技術を使用して実行することができ、鎖再合成を含む、固定化および増幅されたインデックス付きの断片の配列を決定する方法は当技術分野で公知であり、例えば、Bignellら(米国特許第8,053,192号)、Gundersonら(WO2016/130704)、Shenら(米国特許第8,895,249号)およびPipenburgら(米国特許第9,309,502号)に記載される。
実施形態1。複数の単一の核または細胞からの核酸を含むシークエンシングライブラリーを調製するための方法であって:第1の複数のコンパートメントに複数の単離された核または細胞を提供する工程であって、各コンパートメントは単離された核または細胞のサブセットを含み、核または細胞は核酸断片を含む工程;
線形増幅媒介物を細胞または核に導入する工程;
線形増幅によって核酸断片を増幅する工程;
核または細胞の各サブセットを処理してインデックス付きの核または細胞を生成する工程であって、処理は、単離された核または細胞に存在する核酸断片に第1のコンパートメント特異的インデックス配列を加えて単離された核または細胞に存在するインデックス付きの核酸をもたらすことを含み、処理は、ライゲーション、プライマー伸長、ハイブリダイゼーション、増幅または転位を含む工程;および
インデックス付きの核または細胞を合わせてプールされたインデックス付きの核または細胞を生成し、それによって複数の核または細胞からシークエンシングライブラリーを生成する工程を含む方法。
核酸断片を含む複数の単離された核または細胞を提供する工程;
線形増幅媒介物を単離された核または細胞に導入する工程;
第1の複数のコンパートメントに単離された核または細胞を分配する工程であって、各コンパートメントは単離された核または細胞のサブセットを含む工程;
線形増幅によって核酸断片を増幅する工程;
単離された核または細胞の各サブセットを処理してインデックス付きの核または細胞を生成する工程であって、処理は、単離された核または細胞に存在する核酸断片に第1のコンパートメント特異的インデックス配列を加えて単離された核または細胞に存在するインデックス付きの核酸をもたらすことを含み、処理は、ライゲーション、プライマー伸長、増幅または転位を含む工程;
インデックス付きの核を合わせてプールされたインデックス付きの核または細胞を生成し、それによって複数の核または細胞からシークエンシングライブラリーを生成する工程を含む方法。
第1の複数のコンパートメントに複数の単離された核または細胞を提供する工程であって、各コンパートメントは単離された核または細胞のサブセットを含み、核または細胞は核酸断片を含む工程;
核または細胞の各サブセットを処理してインデックス付きの核または細胞を生成する工程であって、処理は、単離された核または細胞に存在する核酸断片に、(i)第1のコンパートメント特異的インデックス配列を加えて単離された核または細胞に存在するインデックス付きの核酸をもたらすこと、および(ii)線形増幅媒介物によって認識されるヌクレオチド配列を加えることを含み、処理は、ライゲーション、プライマー伸長、ハイブリダイゼーション、増幅または転位を含む工程;
線形増幅媒介物を細胞または核に導入する工程;
線形増幅によって核酸断片を増幅する工程;および
インデックス付きの核または細胞を合わせてプールされたインデックス付きの核または細胞を生成し、それによって複数の核または細胞からシークエンシングライブラリーを生成する工程を含む方法。
各サブセットをトランスポソーム複合体と接触させる工程であって、各コンパートメントの中のトランスポソーム複合体は、他のコンパートメントの中の第1のインデックス配列と異なる第1のインデックス配列を含む工程;および
サブセットの中の核酸を複数の核酸に断片化し、第1のインデックス配列をその核酸の少なくとも1つの鎖に組み込んでインデックス付きの核酸を含むインデックス付きの核または細胞を生成する工程を含む、実施形態1~12のいずれか1つの方法。
各サブセットをトランスポソーム複合体と接触させる工程であって、各コンパートメントの中のトランスポソーム複合体は、第1のコンパートメント特異的インデックス配列を含む工程;および
サブセットの中の核酸を複数の核酸に断片化し、第1のコンパートメント特異的インデックス配列をその核酸の少なくとも1つの鎖に組み込んでインデックス付きの核酸を含むインデックス付きの核または細胞を生成する工程を含む、実施形態1~19のいずれか1つの方法。
プールされたインデックス付きの核または細胞のサブセットを、第2の複数のコンパートメントに分配する工程;および
第2のコンパートメント特異的インデックス配列をインデックス付きの核酸に導入して二重インデックス付きの核酸を含む二重インデックス付きの核または細胞を生成する工程であって、導入することは、ライゲーション、プライマー伸長、増幅または転位を含む工程をさらに含む、実施形態1~23のいずれか1つの方法。
プールされた二重インデックス付きの核または細胞のサブセットを第3の複数のコンパートメントに分配する工程;および
第3のコンパートメント特異的インデックス配列をインデックス付きの核酸に導入して三重インデックス付きの核酸を含む三重インデックス付きの核または細胞を生成する工程であって、導入することは、ライゲーション、プライマー伸長、増幅または転位を含む工程をさらに含む、実施形態1~24のいずれか1つの方法。
複数の増幅部位を生成するのに適する条件の下で増幅部位を含む表面をインデックス付きの断片と接触させ、各々は複数のインデックスを含む個々の断片からのアンプリコンのクローン集団を含む工程をさらに含む、実施形態1~31のいずれか1つの方法。
(a)複数の細胞からの単離された核を提供する工程;
(b)単離された核の完全性を維持しつつヌクレオソームが枯渇された核を生成する化学処理に単離された核を付す工程;
(c)ヌクレオソームが枯渇された核のサブセットを第1の複数のコンパートメントに分配し、各サブセットをトランスポソーム複合体と接触させる工程であって、各コンパートメントの中のトランスポソーム複合体がトランスポザーゼおよび他のコンパートメントの中の第1のインデックス配列と異なる第1のインデックス配列を含む工程;
(d)ヌクレオソームが枯渇された核のサブセットの中の核酸を複数の核酸断片に断片化し、第1のインデックス配列を核酸断片の少なくとも1つの鎖に組み込んでインデックス付きの核酸断片を含むインデックス付きの核を生成する工程であって、インデックス付きの核酸断片はトランスポザーゼに付着したままである工程;
(d)インデックス付きの核を合わせてプールされたインデックス付きの核を生成する工程;
(e)プールされたインデックス付きの核のサブセットを第2の複数のコンパートメントに分配し、各サブセットをヘアピンライゲーション二重鎖と、インデックス付きの核酸断片の片方または両方の末端へのヘアピンライゲーション二重鎖のライゲーションに適する条件の下で接触させて二重インデックス付きの核酸断片をもたらす工程であって、ヘアピンライゲーション二重鎖が他のコンパートメントの中の第2のインデックス配列と異なる第2のインデックス配列を含む工程;
(f)二重インデックス付きの核を合わせてプールされたインデックス付きの核を生成する工程;
(g)プールされた二重インデックス付きの核のサブセットを第3の複数のコンパートメントに分配する工程;
(h)二重インデックス付きの核を溶解させる工程;
(i)他のコンパートメントの中の第3のインデックス配列と異なる第3のインデックス配列を含むように二重インデックス付きの核断片を処理する工程;および
(j)三重インデックス断片を合わせ、それによって複数の単一の細胞からの全ゲノム核酸を含むシークエンシングライブラリーを生成する工程
を含む方法。
単一細胞ゲノムシークエンシングのための従来の方法は、均一性およびスループットに関して制限される。ここでは、「sci-L3」、単一細胞コンビナトリアルインデクシング(「sci」)および線形(「L」)増幅を組み合わせたハイスループットの高カバレージ単一細胞シークエンシング方法を記載する。sci-L3方法は、スループットの指数関数的増加を可能にしつつ増幅バイアスを最小にする、一方向性3レベル(「3」)インデクシングスキームを採択する。我々は、単一細胞全ゲノムシークエンシング(「sci-L3-WGS」)、標的化ゲノムシークエンシング(「sci-L3-標的-seq」)ならびにゲノムおよびトランスクリプトームの共アッセイ(「sci-L3-RNA/DNA」)の概念実証を通したsci-L3フレームワークの一般化の可能性を実証する。F1ハイブリッド雄マウスからの>10,000個の精子および精子前駆体のゲノムをプロファイリングし、86,786個の交差をマッピングし、雄減数分裂における全ゲノム均等染色体分離の事例を含む稀な染色体誤分離事象を特徴付けるために、我々はsci-L3-WGSを適用する。CRISPR摂動およびゲノム安定性の測定を連結するために、組換え状況を完全に特徴付けるために、およびハイスループットの高カバレージ単一細胞ゲノムシークエンシングを必要とする他の目標に、sci-L3アッセイを適用することができることが期待される。
現代の単一細胞ゲノムシークエンシング技術は、2つの主要な限界を有する。第1に、ほとんどの方法は個々の細胞を区分することを必要とし、それはスループットを制限する。第2に、ほとんどの増幅方法はPCRベースであり、したがって、指数関数的増幅バイアスを被る。第1の問題を解決するために、我々および同僚は、単一細胞コンビナトリアルインデクシング(「sci-」)を開発し、ここでは、単一細胞の核酸内容を特異的にタグ付けするために数回のスプリット-プール分子バーコード化を実行し、それによって各連続した回のインデクシングによりスループットの指数関数的増加を可能にする。sci-方法は、多数の単一細胞において、クロマチンアクセス性(sci-ATAC-seq)、トランスクリプトーム(sci-RNA-seq)、ゲノム(sci-DNA-seq)、メチローム(sci-MET)、染色体高次構造(sci-Hi-C)をプロファイリングするように首尾よく開発されている(Cao et al., 2017;Cusanovich et al., 2015;Mulqueen et al., 2018;Ramani et al., 2017;Vitak et al., 2017)。第2の問題を解決するために、T7ベースの転写による線形増幅は、単一細胞アッセイとの関係で以前に利用された潜在的解法を提供する(Eberwine et al., 1992;Hashimshony et al., 2012;Sos et al., 2016)。例えば、最近、Chenらは、ゲノムを断片化し、インビトロ転写(IVT)のためにT7 RNAプロモーターを同時に挿入するためにTn5トランスポゾンを使用する、トランスポゾン挿入による線形増幅を開発した(「LIANTI」)。DNA鋳型から生成されるRNAコピーは、さらなる増幅のために鋳型の役割をすることができない。したがって、全てのコピーは、元のDNA鋳型に直接的に由来する。指数関数的増幅を回避することによって、LIANTIは均一性を維持し、配列のエラーを最小にする。しかし、この方法は、各単一細胞からの連続的ライブラリー調製を必要とするので、ロースループットである(Chen et al., 2017)。
スループットを増加させつつ増幅バイアスを最小にする潜在的技術的経路は、「sci」および「LIANTI」方法を単に組み合わせることであろう。しかし、T7プロモーターがTn5トランスポゾンを通して挿入されるLIANTIの分子構造は、わずか2回だけの細胞バーコード化の機会を与え、それは、スループットを1実験につき数千の単一細胞に制限するだろう。それは、ゲノムDNAのプロファイリングにさらに制限される(Chen et al., 2017;Sos et al., 2016)。sci-L3の開発では、単一細胞コンビナトリアルインデクシング、線形増幅およびライゲーションによりT7プロモーターを導入することによる3回の細胞バーコード化(「3レベル」)を統合した(図3A)。「sci」および「LIANTI」を単に組み合わせることと比べて、sci-L3アプローチはいくつかの主要な利点を有する。第1に、潜在的スループットは、3レベルインデクシングによってかなり低減された費用で1実験につき1,000,000個を超える細胞に指数関数的に増加する(Cao et al., 2019)。第2には、単一細胞バーコード化の一方向性の性質は、CRISPR摂動および生じるゲノム不安定性を連結すること、ならびに多数の単一細胞にわたって特異的ゲノム遺伝子座を配列決定することが望ましい他の適用を可能にする全ゲノムシークエンシング(「WGS」)に加えて、標的化シークエンシング(「標的-seq」)にsci-L3を容易に変換することを可能にする。第3に、一般化可能である線形増幅およびハイスループット細胞バーコード化スキームとして、sci-L3は、sci-L3ベースの単一細胞RNA/DNA共アッセイの我々の概念実証によってここで実証される通り、プロトコールの小さい改変によって他の単一細胞アッセイおよび共アッセイに適応させる柔軟性を提供する。
sci-L3-WGSおよびsci-L3-標的-seqの概念実証
sci-L3-WGSおよびsci-L3-標的-seqの3レベルコンビナトリアルインデクシングおよび増幅スキームは、図3Aに示す:(i)細胞をホルムアルデヒドで固定し、ヌクレオソームはSDSによって枯渇させる(Vitak et al., 2017)。結果として生じる核は、次に24個のウェルに均一に分配する。(ii)インデックス付きのTn5の挿入(「タグメンテーション)」)により1回目のバーコードを24ウェルの各々の中に加える。Tn5トランスポゾンがバーコードなしのT7プロモーターを含有するLIANTIと異なって、スペーサー配列がバーコードの5’側に含まれ、それは以降のライゲーション工程のための「ランディングパッド」の役目をする(Tn5トランスポゾン設計の詳細については、図4および実施例2、「sci-L3-WGSおよびsci-L3-標的-seqの方法および分子設計」のセクションを参照する)。(iii)核の全てをプールし、64個の新しいウェルに均一に再分配する;2回目のバーコードをライゲーションによって加え、それは両方のバーコードの外側に配置されるT7プロモーター配列を含む。(iv)核の全てを再度一緒にプールし、蛍光標示式細胞分取(FACS)サイトメトリーによって選別し、1ウェルにつき最高300細胞で最終回のウェルに分配する。異なる倍数性の核をDAPI(4’,6-ジアミジノ-2-フェニルインドール)染色によってゲーティングし、富化することができることに留意する。さらに、単純な希釈は、損失率を低減することができる、FACSの代替手段である。(v)選別された核を溶解し、in situギャップ伸長に付して二重鎖T7プロモーターを形成する。これの後に、IVT、逆転写(RT)および第2鎖合成(SSS)が続き、ゲノムを線形増幅する。3回目のバーコードを特異な分子識別子(UMI)と一緒にSSS工程の間に加えて、個々のIVT転写物をタグ付けする。(vi)それらの起源細胞を規定する3つのバーコードを各々含有する二重鎖DNA分子(図3B、上)は、目標が単一細胞WGS(例えば、ライゲーション(図3B、中央)またはタグメンテーションによって配列アダプターを付け加える)である場合は従来のライブラリー構築方法に、または、目標が単一細胞標的化DNA-seq(例えば、プライマーの1つが標的特異的であるPCR工程を加える(図3B、下))である場合はわずかに改変された方法に適合する。
プロトコールへの小さい改変によってsci-L3-WGSスキームを分子生物学の他の態様に潜在的に適合させることができると認識した。これを実証するために、sci-L3-RNA/DNA共アッセイの概念実証実験を実行した。簡潔には、sci-L3-WGSにおけるように1回目のDNAバーコード化はTn5挿入によって実行されるが、1回目のRNAバーコード化を同時に実行し、逆転写を通してバーコードおよびUMI含有ポリTプライマーによってmRNAをタグ付けする(図6A)。Tn5挿入およびRTプライマーの両方は2回目のバーコードならびにT7プロモーターのライゲーションを媒介することができるオーバーハングを有し、3レベルインデクシングおよび以降のIVTベースの線形増幅を、sci-L3-WGSとほとんど同一の様式で効果的に可能にする(図6A~6B、実施例2、「sci-L3-RNA/DNA共アッセイの方法および分子設計」のセクション)。概念実証として、マウス細胞を2つのヒト細胞系からの細胞と一緒に混合し、sci-L3-RNA/DNA共アッセイを実行した。大多数の細胞のために、RNA(5.2%の衝突率)およびDNA(6.6%の衝突率)の両方について、読み取りデータはマウスまたはヒトゲノムにマッピングされた(図6C~6D)。さらに、良好な共アッセイと一貫して、細胞の100%はそれらのRNAおよびDNAプロファイルによって同じ種標識を割り当てられた。さらなるチェックとして、それらのRNAプロファイルに基づいてt-SNE空間の中のヒト細胞に可視化した。予想通り、それらは2つのクラスターに分離した。Y染色体の存在の有無に基づく個々の細胞の標識化は、96.5%の精度でクラスターをBJ細胞(雄)またはHEK293T細胞(雌)(図6E)に対応すると首尾一貫して同定した。
正常な有糸分裂細胞分裂では、二倍体染色体は複製を経て4コピーのDNAを生成し、姉妹染色分体は相反的な娘細胞に分離する。娘細胞は各々母および父から引き継いだDNA配列の1コピーを受け取り、動原体近位の配列でほとんどいつでもヘテロ接合性を維持する(図7A)。稀に、染色体は染色体相同体の間で有糸分裂交差を経るが、2つの組み換えられた染色分体が異なる娘細胞に分離する場合、それは交差に対し動原体遠位の配列でヘテロ接合性喪失(LOH)の二倍体細胞を時にはもたらすことができる(図7B~C)。
不妊性(B6×Spret)交雑からのM2細胞における染色体分離
上記の通りに得られた、不妊性(B6×Spretus)F1雄の副睾丸からの細胞における減数分裂を分析しようと先ず努めた。2つのsci-L3-WGS実験にわたって、2,689個(>10kの生の読み取りデータを有する2,919個の選別された細胞の92%)および4,239個(>30kの生の読み取りデータを有する4,497個の選別された細胞の94%)の単一細胞のゲノムをプロファイリングした。特異的マッピングの読み取りデータの数を、図5Fに示す。2つのライブラリーのための1.6×および1.4×のシークエンシング深度で(詳細は図5にある)、それぞれ0.7%および1.4%の中間ゲノムカバレージに対応する、1細胞あたり約70kおよび約144kの特異なTn5部位の中間値を得た。
我々は、均等分離が種内(B6×Cast)F1雄の妊性後代においてもMIの間に起こるかどうか疑問に思った。上に示す通り、この交雑からの副睾丸はほとんど完全に1C成熟精子からなる。したがって、全精巣から2C二次精母細胞のために富化した。次に、副睾丸および精巣の両方からの細胞で、sci-L3-WGSを実行した。
これをより正式に考慮するために、各実験からのデータを3つの二項分布のベイズの有限混合物にあてはめた。詳細は、実施例2、「非1C細胞の3集団をあてはめるための有限混合物モデル」のセクションおよび図12に提供され、重要な結論はここで要約される。1番目に、種内(B6×Cast)F1雄の精巣(すなわち、バーコード第2群)からの非1C細胞は、還元(28%)対均等(2%)分離する細胞サブセット、ならびにおそらく1Cダブレット(69%)を含むと推定される(図12B)。この割合は、種間(B6×Spret)F1雄からのM2細胞では異なり、これらは、還元(66%)対均等(14%)分離する細胞サブセット、ならびにおそらく1Cダブレット(20%)を含むと推定される(図12C)。これらの分析は、不妊性(B6×Spret)の交雑が、還元分離より均等分離にバイアスされる細胞のかなりより高い割合を有するという結論を裏付ける。
次に、交差事象のゲノム相関を調査することを目指した。全体として、我々は、(B6×Spret)交雑からの19,601個の交差ブレークポイントを有する1,663個の1C細胞および4,184個の交差ブレークポイントを有する240個のM2細胞、ならびに、(B6×Cast)交雑からの60,755個の交差ブレークポイントを有する5,547個の1C細胞および2,246個の交差ブレークポイントを有する115個のM2細胞を分析した。我々の知る限りでは、これは、哺乳動物の減数分裂に関連して同定された交差事象の数に関して先例のないデータセットである。
交差ホットニスを調節するゲノムフィーチャー
より微細なスケールで交差の分配を評価するために、各マウス染色体に沿って「ホットニスマップ」を生成するために、全ての交差事象を崩壊させた。これらのマップを、減数分裂のDSBホットスポットを最も高い分解能で同定する、一本鎖DNAシークエンシング(SSDS)マップ(Brick et al., 2018;Smagulova et al., 2011, 2016)およびSpo11オリゴヌクレオチド複合体マップ(Lange et al., 2016)と先ず比較した(図15A)。これらの2つのマッピング方法からのB6株におけるDSBマップは、100kbウィンドウ(ロー=0.87、p<2e-308)に沿ってお互いと強く相関する。我々の1CおよびM2細胞の交差パイルアップはお互いと相関するが((B6×Spret)交雑ではロー=0.67、(B6×Cast)交雑ではロー=0.55、両方でp<2e-308、図15B~C)、両方ともDSBマップから逸脱する。関連して、ホットスポット特異化のための主要プレーヤーであるPRDM9遺伝子は、分かれたマウス系統の間で、さらにマウスの亜種の間でさえ、異なるモチーフに結合するように発展した(Davies et al., 2016;Gregorova et al., 2018)。我々は、実施例2、「交差ホットニスに及ぼすPRDM9の影響」のセクションにおける2つの交雑の間の差に及ぼすその潜在的影響を検討する。
1CおよびM2細胞は交差パイルアップにおいて広く類似しているように見えるが(図15)、単一細胞のサブセットにおける交差分配に影響するフィーチャーに何らかの構造があるかどうかについて疑問に思った。これを探究するために、78個のフィーチャーの各々のために各単一細胞について交差関連の情報を集めた(実施例2、「バイオインフォマティクスおよび統計分析の方法」のセクション)。次に、各行を1つの単一細胞とし、各列を1つの要約されたフィーチャー値としたマトリックスの上で、主成分分析(PCA)を使用した。(B6×Spret)交雑については、最初の2つの主成分(PC)は分散の26%を捕捉し、(B6×Cast)交雑については、PC1およびPC3は分散の17%を捕捉する。両交雑では、1CおよびM2細胞は、これらのPCによって2つのクラスターに分離される。図21および図22に、これらのPCに投影される各フィーチャーをプロットする。交差の染色体分配、単為生殖染色体および染色体四分位における交差の位置は、1CおよびM2細胞の分離を促進するようであるフィーチャーである。
最後に、我々は、交差位置の予測モデルを構築するために、ここで観察される多数の事象を活用することを目指した。具体的には、1は同じ地帯長分配からのゲノムから採取した交差地帯であり0はランダム地帯である、二元応答の線形モデルを構築した(詳細は、実施例2、「バイオインフォマティクスおよび統計分析の方法」のセクションにある)。BMA分析におけるのと同じ76個のフィーチャーを使用して、(B6×Spret)交雑のために、0.73の平均レシーバーオペレーター曲線(ROC)の曲線下面積(AUC)のヘルドアウトデータで、交差地帯を予測することができる。BMAによって同定される高い包含確率(MIP>0.5)の25変数のサブセットで、0.72の類似の平均AUCを達成する(図16E)。同様に、(B6×Cast)交雑については、全てのフィーチャーまたはMIP>0.5の25個のフィーチャーのサブセットを使用するとき、0.85の平均AUCを達成する(図16F)。
ここでは、3レベルの単一細胞コンビナトリアルインデクシングおよび線形増幅を組み合わせたフレームワーク、sci-L3を記載する。我々は、sci-L3が単一細胞全ゲノムシークエンシング(sci-L3-WGS)、単一細胞標的化DNAシークエンシング(sci-L3-標的-seq)ならびにゲノムおよびトランスクリプトームの単一細胞共アッセイ(sci-L3-RNA/DNA)に適用可能であることを実証する。sci-L3-WGSにより、少なくとも数万および潜在的に数百万の単一細胞のゲノムを2日間の実験で、10kの細胞では1細胞につき0.14ドル、および1Mの細胞では1細胞につき0.008ドルのライブラリー構築費用で処理することができる。sci-L3-WGSのスループットは、「イン-チューブ」LIANTIなどの線形増幅に基づく代替の単一細胞WGS方法より数桁高い(Chen et al., 2017)。各単一細胞から回収される特異分子の数を、低い数千(Pellegrino et al., 2018)または低い数万(Vitak et al., 2017)から数十万にそれはさらに向上させる。
バーコード第2群からの(B6×Cast)雑種から回収された非1C細胞は、1Cダブレット、均等分離の方へバイアスされているようである細胞、および還元分離の方へバイアスされているようである細胞を含む。それらの相対的な割合を定量化するために、0.01、0.48および0.95の均等分離染色体の確率ならびに0.28、0.69および0.02の混合割合で、3つの二項分布の混合物にデータをあてはめた(図12A)。対照的に、バーコード第1群からの非1C細胞を3つの二項分布の混合物に同様にあてはめようとするとき、我々は0.46、0.5および0.53(全て0.5に近い)の均等分離染色体の確率ならびに0.24、0.44および0.31の混合割合を得る(図12B)。
(B6×Spret)交雑からの19,601個の交差ブレークポイントを有する1,663個の1C細胞および4,184個の交差ブレークポイントを有する240個のM2細胞、ならびに(B6×Cast)交雑からの60,755個の交差ブレークポイントを有する5,547個の1C細胞および2,246個の交差ブレークポイントを有する115個のM2細胞に基づき、我々は、染色体全域にわたって減数分裂交差の分配を先ず考慮した。交差密度は、1Mbあたりの1回の分裂あたりの1細胞あたりの交差の平均数に2(1C細胞)または1(M2細胞)を掛け算したものとここで規定される。(B6×Spret)交雑では、1C細胞における染色体サイズと交差密度の間で強力な負の相関が観察された(図13A、r=-0.66、p=0.002)。以前の知見(Lange et al., 2016)と一貫して、この相関はSpo11オリゴヌクレオチド複合体密度(r=-0.46、p<0.05)によって部分的に説明されるだけであり、これは、より小さい染色体がより多くのDSBを維持し、それらのDSBは交差を生成する可能性がより高いことを示唆する。この負の相関は、M2細胞においてさらにより強力である(図13B、r=-0.83、p=1e-5)。図10A~Bでは、我々は1細胞あたりの1染色体あたりの複数の交差の場合を単一の事象と考え、それは負の相関をさらに強化する(1C細胞ではr=-0.87、p=2e-6;M2細胞ではr=-0.91、p=8e-8)。これらの観察は、より小さい染色体が交差に対して、特に1回の細胞分裂につき少なくとも1つの交差を有することに対してよりホットであることを示唆する。(B6×Cast)交雑において、同じ傾向が観察される(図14A~D)。1C細胞は種間および種内交雑について1細胞あたりの1染色体につき平均0.62個および0.58個の交差をそれぞれ有したが、M2細胞は1細胞あたりの1染色体につき平均0.92個および1.03個を有した(図13C~D、10C~D)。種間M2細胞における交差率は、2%の配列分散にもかかわらず、B6近交系マウスの4C精母細胞においてMlh1焦点によって測定した交差数(Froenicke et al., 2002)よりわずか9%低い。1C細胞における交差率は、単一のヒト精子シークエンシングにおいて観察されるより45%低い(Lu et al., 2012;Wang et al., 2012)。後者の差は、マウス染色体の端部動原体の性質に大部分よるかもしれない。種間(B6×Spret)交雑は(B6×Cast)交雑と比較して1Cにおいて検出される交差のより高い平均数を有するが(p=7e-26、マン-ホイットニー検定法)、M2細胞における交差の平均数はより低い(p=2e-10)。あらゆる染色体に交差を有する19個の常染色体の全てを還元分離させたM2細胞の割合は、(B6×Spret)交雑(41/80または51%)より(B6×Cast)交雑(60/91または66%)で高いことに注意し(p=0.06、フィッシャーの直接確率検定)、それは後者の不妊に寄与する可能性がある。
ゲノム全体にわたって染色体に沿って交差ブレークポイントを積み上げることによる交差ホットニスマップに基づいて(図15)、おそらくCast PRDM9対立遺伝子がF1雑種において半優性である結果として、種内(B6×Cast)交雑では、交差ホットニスがB6雄よりもCast雄においてマッピングされたDSBホットドメインとよりよく相関する(それぞれ、ロー=0.28および0.12、p<2e-308およびp=1e-83)ことを見出した。相関は、(B6×Cast)F1動物においてマッピングされたDSBホットドメインでより強い(ロー=0.3、p<2e-308)。(B6×Spret)交雑では、PRDM9コンセンサス結合性部位の浸食は、Spo11オリゴヌクレオチド複合体マップによって規定される4種類のDSBホットスポットをもたらす:B6とSpretの間で保存される、「対称」ホットスポットと呼ばれるもの、B6またはSpretにだけ存在する、「非対称」ホットスポットと呼ばれるもの、およびいずれの種においてもPRDM9結合性部位を含有しないもの。4種類のDSBホットドメインの全ては、(B6×Spret)交雑からの交差と十分に相関しない(B6にマッピングされた全てのSpo11ホットスポットを使用した場合ロー=0.13、p=4e-87;「対称ホットスポット」だけを使用する場合ロー=0.11、p=3e-63)。1つの可能性は、(B6×Spret)交雑におけるDSB部位がSpretPRDM9対立遺伝子によって強く支配され、そのため、B6系統バックグラウンドにマッピングされたDSBホットスポットから交差部位が予測されないということである。
我々は逆分離の高い発生率を、特に種間(B6×Spret)交雑において観察した。我々は、以下について推測する:1)何が未熟な動原体コヒーシン分離を引き起こす可能性があるか、2)1つの交差が適切な還元分離に十分であるかどうかについて、および3)MIにおける均等分離はどんな結果をもたらす可能性があるか。
交差ホットニスは連続帯であり、多くの因子によって形づくられる。(B6×Cast)交雑における交差は、2つの親系統のための個々のマップにおけるより、F1交雑にマッピングされる減数分裂のDSBホットスポットと強く相関し、これは、新規の減数分裂ホットスポットがF1雑種において形成することができるという以前の知見に基づいて予想される(Smagulova et al., 2016)。(B6×Spret)交雑では、交差は、弱いがSpo11切断と正に相関する。Spo11マップは、B6対立遺伝子のPRDM9タンパク質に結合したPRDM9部位だけを説明することに留意すべきであり、PRDM9のSpretコピーが異なる部位に結合して、我々の分析で説明されていない新しい減数分裂DSBホットスポットを形成する可能性がある。減数分裂交差と正に相関していることが観察されるゲノムフィーチャーには、GCに富む領域(酵母減数分裂における場合も同様である(Petes, 2001;Petes and Merker, 2002))、系統(Lilue et al., 2018)、遺伝子体、偽遺伝子転写物、CTCF結合性部位、複製ドメイン(Marchal et al., 2018)、DNAトランスポゾン、サテライトDNAならびにH3K4me1、H3K27me3およびH3K36me3を含むヒストン改変のサブセット(Mu et al., 2017)の間のCNVの増加が含まれる。興味深いことに、雄生殖細胞における有糸分裂から減数分裂への転換の調節に関与するDmrt6の結合性部位(Zhang et al., 2014)は、減数分裂の交差ホットニスと強く相関している。注目すべきことに減数分裂交差と負に相関しているゲノムフィーチャーには、3’UTR、LINEおよび低い複雑性のDNAが含まれる。rDNAが減数分裂交差に対して極めて冷めている(Petes and Botstein, 1977)酵母の場合と異なり、マウスrDNAは交差を抑制しないようである。これらのゲノムフィーチャーで、(B6×Spret)および(B6×Cast)でそれぞれ0.73および0.85の精度によって、真の減数分裂交差開始部位をマウスゲノムの中のランダムに採取された地帯から区別することができ、(B6×Cast)交雑における0.85の予測精度は25個のゲノムフィーチャーのサブセットに適用される。様々なフィーチャーがモデル化アプローチの間でほとんど一貫してふるまうが、さらなる実験なしでいかなる因果律も割り当てることができないことを我々は強調する。
sci-L3-WGSおよびsci-L3-標的-seqの方法および分子設計
単一細胞の調製およびヌクレオソーム枯渇
細胞懸濁液は、シャーレからのトリプシン処理または組織からの均質化によって調製される。ワシントン大学IACUCが承認したプロトコールに従って、雄F1マウスはCO2および続く頚椎脱臼によって安楽死させた。雄生殖細胞の単離のために、10%FBSを追加した1×PBSの1mlの中で管をスライスし、室温で15分間、組織をインキュベートすることによって精巣上体を解剖した。インキュベーションの後、ピペット操作によって細胞懸濁液を収集した。精巣上体から単離された細胞は、(B6×Spret)交雑の実験のために、および、(B6×Cast)交雑の成熟精子(「バーコード第3群」)の供給源としても使用した。(B6×Cast)交雑のための2C細胞のための富化方法としての全精巣からの核の単離のために、精巣細胞を先ず1%ホルムアルデヒドによって架橋させ、低張性緩衝液を使用して核を抽出した。次に、主にDAPIシグナルに基づくDNA含有量によって、1Cおよび2C核をFACS選別した。培養したヒトおよびマウス細胞を4℃で5分間、550gでペレット化し、雄生殖細胞は4℃で10分間、2400gでペレット化する。
トランスポゾンDNAオリゴはリン酸化された2本の鎖の両方の5’で合成され、1つはLIANTIおよびNexteraにおけるのと同様にTn5挿入のために必要とされ(5’/リン酸/CTGTCTCTTATACACATCT、IDT、PAGE精製(配列番号1))、他は、ライゲーションのために必要とされる(5’/リン酸/GTCTTG XXXXXXXX[1回目のバーコード]AGATGTGTATAAGAGACAG、IDT、標準の脱塩(配列番号2))。アニール緩衝液(10mMトリス-HCl pH8.0、50mM NaCl、1mM EDTA、pH8.0)の中での段階的な冷却(95℃5分間、-0.1℃/サイクル、9秒間/サイクル、25℃まで700サイクル)による1:1のアニーリングの後に、5’オーバーハングを有するTn5二重鎖を1.5μMに希釈する。次に7.2μLの保存緩衝液(50%グリセロールを含む1xTE)を12μL~1μMのTn5トランスポザーゼ(Lucigen、TNP92110)に加え、0.79μLの希釈されたトランスポザーゼを室温で30分間、0.4μLの1.5μM Tn5二重鎖と一緒にインキュベートする。トランスポソームは、0.2μMの最終濃度に二量体化する。トランスポソーム複合体は、最高1年の間、-20℃で安定して保存することができる。1回目に24ウェルのバーコード化のために24個の反応を立ち上げたが、適用によってはより多くのウェルが望ましいかもしれない。新しい生物適用ごとに、検査実験のためにトランスポソームを先ず0.1μMまでさらに希釈する。特異な読み取りデータの数およびライブラリー複雑性はより最適でないが(図5)、低分解能でのマッピングのために使用できる。
次に、lo-bind96ウェルプレートの各ウェルの中に1.5μLの核を20,000/μL濃度で分配し、6.5μLのH2Oおよび0.7μLの50mM MgCl2を加える(溶解緩衝液の中のEDTAは3.24mMの最終濃度である)。上で調製される1.2μLのトランスポソームを各ウェルに加え、次にプレートを55℃で20分間インキュベートする(サーモミキサーが推奨されるが、要求されていない)。次に5μLの停止液(40mM EDTAおよび1mMスペルミジン)を加え、核をトラフにプールする。溶解緩衝液のさらなる1mLを核懸濁液に加えてからペレット化する。上清を注意深く除去した後に、312μLの再懸濁緩衝液(24μLの10mM dNTP、48μLの10×タグメンテーション緩衝液[50mM MgCl2、100mMトリス-HCl pH 8.0]、96μLのH2O、144μLの溶解緩衝液)に核を再懸濁し、新しいlo-bind96ウェルプレートの各ウェルに4.7μLの核ミックスを分配する。ヘアピンライゲーション二重鎖(1.CAAGAC2.Y’Y’Y’Y’Y’Y’Y’[2回目のバーコードの逆相補体]3.CAGGAGCGAGCTGCATCCC4.AATTTAATACGACTCACTATA5.GGGATGCAGCTCGCTCCTG6.YYYYYYY[2回目のバーコード](配列番号3))は、Tn5トランスポゾン二重鎖と同様にプレアニールされ、1.5μMまで希釈される。ライゲーション二重鎖は5つのエレメントを含有することに留意する:1)Tn5の上のライゲーションアダプターの逆相補体;2)2回目のバーコードの逆相補体;3)第2の鎖の合成(SSS)プライマーの逆相補体;4)T7プロモーター、これはヘアピンのループ領域であることに留意する;5)T7転写を強化するためのGGGから開始する第2の鎖の合成(SSS)プライマー領域(図4Bの中の「sp2」);6)2回目のバーコード(図4Bの中の「bc2」)。これらの二重鎖の0.8μLを、核懸濁液を有する64ウェルの各々に加え、1.18μLのライゲーションミックス(0.6μLの10×NEB T4リガーゼ緩衝液、0.48μLのPEG-4000、0.1μLのT4 DNAリガーゼ[Thermo EL0011])を各ウェルに加え、20℃で30分間インキュベートする。ライゲーションの後、ループ状の構造はLIANTIのそれに類似し、RT工程(下で議論する)での効率を促進する点に、および、両ラウンドのバーコードがT7プロモーターの3’に存在し、したがって増幅された分子に含まれる点に留意する。ライゲーション反応は、4μLの停止液を加えることによって中止される。細胞は、次に新しいトラフ(約630μL)にプールされ、5μg/mLの最終濃度のDAPIで染色され、細胞選別の前に加えられる3μLの溶解緩衝液によって各新しいウェルの中に100~300で選別される。FACSによる各選別事象がノズルのサイズによって約3~5nLのFACS緩衝液に関連する点に留意する、我々は、塩濃度を低くしておくために各ウェルに加えられる液体の総量を<1μLに保つことを推奨する。
次に、細胞溶解のために、75℃で45分間インキュベートし、4℃に冷却し、新たに希釈したQiagenプロテアーゼ(最終濃度2mg/mL)によって55℃で8時間処理することによって、各ウェルの中の合計3.5~4μLの選別された核にとりかかる。75℃で30分間インキュベートすることによって、プロテアーゼを次に熱不活性化する。細胞溶解物は、-80℃で保存することができる。以降の増幅工程がRNAを含み、時間に敏感であるので、各実験のために32個以下のウェルの試料(約9600個の単一細胞)を処理することを推奨する。ギャップ伸長(図4C)のために、2μLのH2O、0.7μLの10×タグメンテーション緩衝液、0.35μLの10mM dNTPおよび鎖置換活性を有する0.35μLのBst WarmStart2.0ポリメラーゼの混合物を加えることによって、鎖置換活性を有するポリメラーゼを使用し、68℃で5分間インキュベートする。ライゲーションが両末端で成功した場合、両側にT7プロモーターがあって二重鎖は対称的であるが、ライゲーションが一方の末端だけで成功した場合は、破線ボックスの中の領域は一方の側で欠落している点に留意する。分子間のライゲーションは、一般的に非効率的である。分子間のライゲーションの必要性を最小にするためにプレアニールされたヘアピンループを含めたが、2つの分子(ヘアピンループのない3つの代わりに)はなお互いを見出す必要がある。ライゲーション効率が50%である場合、両末端にライゲーションを有する率は25%であるが、いずれかの末端にライゲーションを有する率は75%である。RT工程の後半では、ライゲーションの成功は一方の末端だけで必要とされることを我々は示す。ギャップ伸長の後、2μLのH2O、2μLのT7 Polミックスおよび10μLのrNMPミックス(NEB、HiScribe(商標)T7 Quick High Yield RNA Synthesisキット)を加えることによって、20μLのT7インビトロ転写システムをアセンブルする。混合物を37℃で10~16時間インキュベートする。
2.2μLの0.5M EDTAを加えることによって、転写を終了する。増幅されたRNA分子は、次にRCC-5(Zymo Research、R1016)で精製し、18μLの0.1×TEで溶出する。最初に0.6μLのRNA RTプライマー(rArGrArUrGrUrGrUrArUrArArGrArGrArCrArG、IDT(配列番号4))、2μLの10mM dNTPおよび0.5μLのSUPERase*In(商標)RNアーゼ阻害剤(20U/μL、Thermo Fisher AM2696)を加えることによって、30μLのRTシステムをアセンブルする。次に、変性および二次構造の除去のために70℃で1分間および90℃で20秒間インキュベートし、氷上で急冷する。6μLの5×RT緩衝液、1.5μLの0.1M DTT、1μLのSUPERase*In(商標)および1μLのSSIVによるRTのために、SuperScript(商標)IV逆転写酵素(SSIV、Thermo Fisher 18090050)を使用する。RT反応は、55℃で15分間、60℃で10分間、65℃で12分間、70℃で8分間、75℃で5分間および80℃で10分間インキュベートする。反応を室温まで冷却させてから、0.5μLのRNアーゼH(NEB)および0.3μLのRNアーゼA(Life Technologies、AM2270)を加え、37℃で30分間インキュベートする。図4Eは、RT工程の間の2つのシナリオを示すことに留意する:1)両末端のライゲーションが成功した場合、LIANTIにおける通りRTは折り返しループによって初回刺激される可能性がある;2)一方の末端だけでライゲーションが成功した場合、変性工程の前に加えられるRNA RTプライマーによってRTが初回刺激される可能性がある。過剰なRNAプライマーおよびRNA転写物は、cDNA合成の後に分解される。最後に、27μLのH2O、20μLの5×Q5緩衝液、20μLのQ5GCエンハンサー、1μLのQ5ポリメラーゼおよび1μLのSSSプライマー(NNNN[UMI]ZZZZZZ[3回目のバーコード]GGGATGCAGCTCGCTCCTG、IDT、標準の脱塩(配列番号5)を加えることによって、Q5 DNAポリメラーゼによって第2の鎖を合成する。結果として生じる二本鎖DNAはDCC-5(ZymoResearch、D4014)で精製することができ、シークエンシングアダプターを加えるための最小限の3サイクルのPCRによってNEBNextウルトラIIなどのライブラリー調製キットで続行することができる。
単一細胞調製およびヌクレオソーム枯渇
下に示す差以外はsci-L3-WGSにおけるのと同じプロトコールで、細胞懸濁液を調製する。HEK293T、BJ-5taおよび3T3細胞をシャーレからトリプシン処理し、1M/mLの細胞濃度で、室温で10分間、1×PBS中の2%PFAで固定した。以降のクエンチング(グリシンによる)、洗浄、核単離(0.1%IGEPALによる)、ヌクレオソーム枯渇(xSDS方法)工程は、1%Superase-Inを全ての溶解緩衝液および1xNEBuffer2.1に加えること以外はsci-L3-WGSと同一である。核は、1μLにつき20,000個の核で、1%Superase-Inを含む溶解緩衝液に再懸濁させる。
単一細胞ゲノム増幅構成成分については、トランスポソームの設計およびアセンブリーは、sci-L3-WGSと同一である。
次に、RT工程の準備のために、lo-bind 96ウェルプレートの各ウェルの中に1.5μLの核を20,000/μL濃度で分配し、0.2μLのH2O、0.3μLの50mM MgCl2(溶解緩衝液中のEDTAを中和するために)、0.25μLの10mM dNTPおよび上記の1μLの25μM RTプライマーを加える。次に、二次構造を除去し、氷上で速やかに失活させるために、核混合物を55℃で5分間インキュベートする。次に、1μLの5×RT緩衝液、0.03μLの100mMDTT(溶解緩衝液からのDTTがある点に留意する、最終濃度5mM)、0.25μLのSSIV、0.25μLのRNaseOUT(Thermo Fisherカタログ番号10777019)を加え、RT反応のために25℃で1分間、37℃で1分間、42℃で1分間、50℃で1分間、55℃で15分間インキュベートする。次に、各ウェルに0.4μLのMgCl2および3.52μLのH2O、および上で調製した1.2μLのトランスポソームを加える。細胞溶解後までの全ての以降の工程は、sci-L3-WGSと同一である。
部分的NEBNext読み取りデータ1プライマーを5’オーバーハングとした(CACGACGCTCTTCCGATCT NNNNNNN(配列番号9))ギャップ伸長のために、ランダム七量体を使用する。1μLの20μMオリゴを加え、DNAを変性させるために95℃で3分間インキュベートし、オリゴがアニールするために室温まで徐々に冷却する(約5分間)。次に、2μLのH2O、0.8μLの10×NEBuffer2、0.4μLの10mM dNTP、0.4μLクレノウ断片(3’→5’エキソ-、NEB M0212S)を加え、30℃で8分間および75℃で10分間インキュベートする。ギャップ伸長の後、20μLのT7インビトロ転写システムを同じsci-L3-WGSプロトコールによってアセンブルする。
異なるオリゴ配列以外、全ての工程はsci-L3-WGSと同一である。IVTの後のRT工程で、0.6μLのRNA RTプライマーを使用する代わりに、0.6μLのNEBNext読み取りデータ1プライマー(AATGATACGGCGACCACCG AGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT、IlluminaシークエンシングのP5末端、IDT(配列番号10))を使用する。SSSプライマーのために、シークエンシングアダプターを加えるために、AAGCAGAAGACGGCATACGAGAT[P7末端]NNNN[UMI2]Z’Z’Z’Z’Z’Z’[3回目のバーコード]CGTCTCTAC GGGATGCAGCTCGCTCCTG(配列番号11)を使用する。結果として生じる二本鎖DNAは今やIlluminaシークエンシングのためのP5およびP7の両末端を含有し、1.1×AmpureXPビーズで精製してシークエンシングに進行することができることに注意すべきである。シークエンシングアダプターを加えるためのライブラリー調製工程およびsci-L3-WGSにおける最小限の3サイクルのPCRは、共アッセイのために不要である。
(B6×Spret)交雑
2つの別個の実験において、それぞれ70日齢および88日齢の(B6×Spret)F1雄からの6つおよび3つの副睾丸から単離した細胞をプールし、1%ホルムアルデヒドで固定した。各実験で、ヌクレオソーム枯渇の後、1ウェルにつき30,000個の細胞を分配し、1回目のバーコードを加えるために24ウェルにわたってin situのインデックス付きTn5挿入を実行した。次に全ての細胞をプールし、ライゲーションによって2回目のバーコードおよびT7プロモーターを加えるために、これらを64個のウェルに再分配した。全ての細胞を再びプールした後に、細胞混合物を1:6に分割し、大部分の細胞(6/7)をFACS選別し、残り(1/7)を希釈した。結果として生じたウェルは、1ウェルにつき100~360個の細胞を含有し、概算4~11%の衝突率であった。
6つの精巣から、約12M個の1C円形精子および約0.5M個の2C細胞を回収した。しかし、1C細胞の>20倍多い数のために、2C細胞のために選別した集団において多くの1C細胞がなお見出された(図8F)。2C細胞を富化しようと試みたsci-L3-WGS実験の1つでは、精巣上体から約160k個の精子、約160k個の1C円形精子および約70k個の2C細胞をタグメントしたと推定し、sci-L3-WGSのFACS工程の間に2C細胞のためにさらに富化した(図8G)。しかし、2回の富化にもかかわらず、1C細胞はなお優勢であった。
読み取りデータ処理、アラインメントおよびSNVコーリング
ベースコールは、インデックスにおけるエラーのために1つのミスマッチを可能にしたbcl2fastqによってfastqファイルに変換した。次に、非多重化するために、カスタマイズしたシェルスクリプト「sci_lianti_v2.sh」(パイソンスクリプトおよびR Markdownファイルは、「sci_lianti_inst.tar.gz」として別々にアップロードされる;全ての主要なおよび補助的図面を生成するための中間データファイルを含有するRパッケージは、以下のリンク:https://drive.google.com/file/d/19NFubouHrahZ8WoblL-tcDrrTlIZEpJh/view?usp=sharingを通してダウンロードし、インストールすることができる)を使用し、それは以下の工程のためにパイソンスクリプトまたはNGSツールを呼び出す:1)全ての単一細胞コンビナトリアルバーコードが読み取りデータ1(R1)にあるように、読み取りデータ対を命じる;2)3回目の(SSS、6nt、エラーは許されない)バーコードを非多重化し、転写物のためのバーコードおよびUMIの両方を読み取りデータ名に付け、3回目のバーコードでライブラリーを分割する。100~300個の単一細胞を含有する個々のライブラリーの3回目のバーコードによる分割のために、全ての以降の工程は並行して実行されることに留意すべきである;3)R1における1回目(Tn5、8nt、1つのエラーが許される)および2回目(ライゲーション、7nt、1つのエラーが許される)のバーコードを分割するためにcutadaptを使用し、エラーはLevenshtein距離によって計算し、両ラウンドのバーコードを読み取りデータ名に付ける。この工程は対形成末端モードで実行される、すなわち、R1が正しいバーコードおよびスペーサー構造を有しない場合は、対形成読み取りデータ2(R2)は棄却される;4)cutadaptを使用してR2をクリーンアップする;5)対形成末端モードでbwa memを有するhg19またはmm10ゲノムに整列させる(Li and Durbin, 2009)。バーコード衝突を調査する実験のために、hg19およびmm10の連結された参照を使用し、特異的に整列させた読み取りデータを使用してヒトまたはマウスゲノムへの相対的なマッピング率を決定する;6)読み取りデータ名に付けられた1回目および2回目のバーコードを使用してbamファイルを単一細胞bamファイルに分割する;7)bedtool(Quinlan and Hall, 2010)によってbamファイルをbedファイルに変換し、R1またはR2が同じエンドポイントを共有する場合、特異な挿入部位を決定する。特異なTn5挿入部位は、読み取りデータ対の両末端が異なる必要性がある断片と規定される;8)「lianti」パッケージ(https://github.com/lh3/lianti/blob/master/pileup.c)(Chen et al., 2017)の「パイルアップ」機能を使用して、対立遺伝子-awareモードで変異体を呼び出す。この工程で、SNPコールが最終vcfファイルに含まれるために各SNP位置の深度の閾値がバルクファイルにおいて超えられることだけが必要であり、したがって、変異体がバルクファイルの中でヘテロ接合のSNPとして存在する限り、REFおよびALT対立遺伝子の生の数値が単一細胞カラムに含まれるように、各単一細胞bamファイルと合わせたバルクbamファイル(全ての約6900個の単一細胞のsamtool併合によって生成される(Chen et al., 2017;Li and Durbin, 2009)、30×を超える)を含めることに留意する。これは、新規のSNPコーリングクエスチョンを遺伝子タイピングクエスチョンに変換することによって、単一細胞における低深度シークエンシングによる高い偽陰性率の問題を回避する;9)Spretのための参照SNP vcfファイル(Mouse Genome ProjectからダウンロードされるSPRET_EiJ.mgp.v5.snps.dbSNP142.vcf.gz)によって各単一細胞においてSNP品質に関して呼び出されるSNVを注釈する。注釈したSNPファイルは、以降の交差ブレークポイント分析のための入力として次に使用する。
所与のSNP部位の遺伝子型は、参照および代替の対立遺伝子を裏づける読み取りデータの数を比較することによって決定される。1C細胞では、交差位置は、3つの状態、参照、代替物およびヘテロ接合体の隠れマルコフモデルをあてはめることによって決定される。
この工程はHMM出力からrdsファイルをとり、単為生殖染色体コールを生成する。
この工程は、図10、13および14に示す減数分裂交差の染色体レベルの特徴を生成する。
均等分離する染色体のそれらの確率を表しているp1、p2、p3によってそれぞれパラメータ化される、3つの二項分布の混合物にデータをあてはめた。これらの3つの二項分布の相対寄与は、長さ3のベクターシータによって表される。Rパッケージrstan(http://mc-stan.org/users/interfaces/rstan)をθのために先行する均一なDirichlet:θ~Dir(K=3、α=1)およびpのために先行するベータ:p~ベータ(a=5、b=5)と使用して、それらの後方分布から試料を描画することによって、p1、p2、p3ならびにθを推定する。モデル仕様に関するさらなる詳細については、Stanファイルmt_mixture_model.stanを参照する。
我々は、様々なゲノムエレメントに関して、以前のゲノム研究からの、ならびにダウンロードされたgff3フォーマットのマウスアノテーションファイルおよびUCSCゲノムブラウザーからのRepeatMasker(https://genome.ucsc.edu/cgi-bin/hgTables)からのデータセットを処理した。mm9に基づくデータセットを、mm10に先ず上げる。これらのデータセットは、2つのカテゴリーにだいたい分類される:bedフォーマットの数値データまたはbedGraphフォーマットの様々な遺伝子もしくは後成的マークのシグナル。細胞クラスタリングおよび予測モデル化では、交差地帯は異なる長さを有する。細胞クラスタリング分析のために各単一細胞における全ての交差から合計される配列の全量を分けることによって数値データを正規化し、および極めて短い地帯が過剰に重み付けされないように、各交差地帯またはランダム採取地帯のために地帯長プラス1kbを分けることによって正規化する。1kbを加えることが大いに余分の配列を含まないように、中間の地帯長は150kbであることに留意すべきである。様々なマークの連続的シグナルを有するデータセットについては、交差またはランダム地帯を横切るマークの平均シグナルをとる。交差パイルアップデータセットについては、均一サイズの100kbのウィンドウを使用したので、数値データを使用するときに地帯長を正規化しなかった。
それらのブレークポイントフィーチャーに基づいて1CおよびM2細胞の分離を2Dで可視化するために、主成分分析を使用する。各単一細胞について、3タイプに対応する合計78個のフィーチャーの交差関連情報を集めた。最初のタイプとして、各細胞の各染色体のために交差または全染色体LOH事象の数を単に計算した。第2のタイプとして、GC含有量、配列分散、クロマチンマークの強度等のフィーチャーのために、各細胞における交差ブレークポイントの中間値を計算した。第3のタイプとして、各細胞における、遺伝子体、末端反復配列(LTR)、交差ブレークポイントと重なったLINEエレメントなどのゲノムエレメントの正規化された数値を計算した。
Claims (9)
- 複数の細胞からの核酸を含むシークエンシングライブラリーを調製するための方法であって、以下:
(a)複数の細胞からの単離された核を提供する工程;
(b)第1の複数のコンパートメントに単離された前記核を分配する工程であって、ここで、各コンパートメントが単離された核のサブセットを含む、工程、および各サブセットをトランスポソーム複合体と接触させる工程であって、ここで、各コンパートメントの中の前記トランスポソーム複合体が、トランスポザーゼ、他のコンパートメントの中の第1のインデックス配列と異なる第1のインデックス配列、および少なくとも1つの遍在配列を含む、工程;
(c)ヌクレオソームが枯渇された核の前記サブセットの中の核酸を複数の核酸断片に断片化し、前記第1のインデックス配列を前記核酸断片の少なくとも1つの鎖に組み込んでインデックス付きの核酸断片を含むインデックス付きの核を生成する工程;
(d)インデックス付きの前記核を合わせてプールされたインデックス付きの核を生成する工程;
(e)プールされたインデックス付きの前記核のサブセットを第2の複数のコンパートメントに分配する工程;
(f)各コンパートメントの中のインデックス付きの前記核酸断片に第2のインデックス配列を組み込んで二重インデックス付きの核酸断片を生成する工程であって、ここで、前記第2のインデックス配列が各コンパートメントに特異である、工程;
(g)二重インデックス付きの前記核を合わせて、プールされた二重インデックス付きの核を生成する工程;
(h)プールされた二重インデックス付きの前記核のサブセットを第3の複数のコンパートメントに分配し、二重インデックス付きの前記核を溶解して、線形増幅により二重インデックス付きの核酸断片を増幅する工程;
(i)各コンパートメントの中の二重インデックス付きの前記核酸断片に第3のインデックス配列を組み込んで三重インデックス付きの核酸断片を生成する工程であって、ここで各コンパートメントの中の前記第3のインデックス配列が、他のコンパートメントの中の第1および第2のインデックス配列と異なり、かつ各コンパートメントの中の前記第3のインデックス配列が、他のコンパートメントの中の第3のインデックス配列と異なる、工程;
(j)三重インデックス付きの前記核酸断片を合わせてプールされた三重インデックス付きの核酸断片を生成し、それにより複数の前記核からのシークエシングライブラリーを産生する工程、
とを含む、方法。 - 単離された前記核が、ヌクレオソーム無しである、請求項1に記載の方法。
- 単離された前記核の完全性を維持しつつヌクレオソームが枯渇された核を生成する条件に、工程(b)の前に、単離された前記核を付すことをさらに含む、請求項1に記載の方法。
- 前記遍在配列が、遍在プライマーを含む、請求項1に記載の方法。
- 工程(c)のインデックス付きの前記核酸断片が、前記トランスポザーゼに付着したままであり、そのため同じゲノムDNA分子に由来する核酸断片が物理的に連結されたままである、請求項1に記載の方法。
- 前記トランスポザーゼをインデックス付きの前記核酸断片から解離させることをさらに含む、請求項5に記載の方法。
- 工程(f)が、インデックス付きの核酸断片の片方または両方の末端へのヘアピンライゲーション二重鎖のライゲーションに適する条件の下で各サブセットをヘアピンライゲーション二重鎖と接触させて、二重インデックス付きの核酸断片をもたらすことを含み、ここで、前記ヘアピンライゲーション二重鎖が、前記第2のインデックス配列を含む、請求項1に記載の方法。
- 工程(h)の後および工程(i)の前に、二重インデックス付きの前記核酸断片の人為操作をさらに含み、ここで前記人為操作が、ギャップ伸長、インビトロ転写、逆転写、またはそれらの組み合わせを含む、請求項1に記載の方法。
- 工程(i)の組み込むことが、第2の鎖の合成を含む、請求項1に記載の方法。
Applications Claiming Priority (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201862673023P | 2018-05-17 | 2018-05-17 | |
| US62/673,023 | 2018-05-17 | ||
| US201962821864P | 2019-03-21 | 2019-03-21 | |
| US62/821,864 | 2019-03-21 | ||
| PCT/US2019/032966 WO2019222688A1 (en) | 2018-05-17 | 2019-05-17 | High-throughput single-cell sequencing with reduced amplification bias |
| JP2019569897A JP7234146B2 (ja) | 2018-05-17 | 2019-05-17 | 低減した増幅バイアスによるハイスループット単一細胞シークエンシング |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019569897A Division JP7234146B2 (ja) | 2018-05-17 | 2019-05-17 | 低減した増幅バイアスによるハイスループット単一細胞シークエンシング |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2023075143A JP2023075143A (ja) | 2023-05-30 |
| JP7638309B2 true JP7638309B2 (ja) | 2025-03-03 |
Family
ID=68541159
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019569897A Active JP7234146B2 (ja) | 2018-05-17 | 2019-05-17 | 低減した増幅バイアスによるハイスループット単一細胞シークエンシング |
| JP2023025703A Active JP7638309B2 (ja) | 2018-05-17 | 2023-02-22 | 低減した増幅バイアスによるハイスループット単一細胞シークエンシング |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019569897A Active JP7234146B2 (ja) | 2018-05-17 | 2019-05-17 | 低減した増幅バイアスによるハイスループット単一細胞シークエンシング |
Country Status (13)
| Country | Link |
|---|---|
| US (2) | US11981891B2 (ja) |
| EP (1) | EP3794141A4 (ja) |
| JP (2) | JP7234146B2 (ja) |
| KR (2) | KR102447811B1 (ja) |
| CN (2) | CN111201329B (ja) |
| AU (3) | AU2019270185A1 (ja) |
| BR (1) | BR112020005982A2 (ja) |
| CA (2) | CA3206377A1 (ja) |
| IL (1) | IL271454B2 (ja) |
| MX (2) | MX2019015013A (ja) |
| RU (1) | RU2744175C1 (ja) |
| SG (1) | SG11202001827UA (ja) |
| WO (1) | WO2019222688A1 (ja) |
Families Citing this family (36)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106987629B (zh) * | 2017-03-31 | 2021-01-15 | 上海市第一妇婴保健院 | 一种在单细胞水平上检测基因组上核小体排布的方法 |
| IL271454B2 (en) | 2018-05-17 | 2025-04-01 | Illumina Inc | Rapid single-cell genetic sequencing with low Aggregation bias |
| 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 |
| JP2023508792A (ja) * | 2019-12-19 | 2023-03-06 | イルミナ インコーポレイテッド | ハイスループット単一細胞ライブラリー、並びに製造方法及び使用方法 |
| WO2021146187A1 (en) * | 2020-01-13 | 2021-07-22 | Fluent Biosciences Inc. | Emulsion based drug screening |
| JP7730449B2 (ja) | 2020-01-13 | 2025-08-28 | イルミナ インコーポレイテッド | 単一細胞遺伝子プロファイリングのための方法およびシステム |
| CN115698282A (zh) | 2020-01-13 | 2023-02-03 | 福路伦特生物科学公司 | 单细胞测序 |
| WO2021147873A1 (en) * | 2020-01-21 | 2021-07-29 | Thunderbio Innovation Ltd | Method for amplifying a target nucleic acid |
| EP4107286A4 (en) * | 2020-02-19 | 2024-04-24 | Sequencing Health, Inc. | METHODS AND SYSTEMS FOR PROCESSING NUCLEIC ACIDS |
| US20230212656A1 (en) * | 2020-02-20 | 2023-07-06 | The Regents Of The University Of California | Methods of spatially resolved single cell sequencing |
| US11866782B2 (en) | 2020-03-16 | 2024-01-09 | Fluent Biosciences Inc. | Multi-omic analysis in monodisperse droplets |
| AU2021240092A1 (en) * | 2020-03-20 | 2022-10-27 | Mission Bio, Inc. | Single cell workflow for whole genome amplification |
| WO2021195169A1 (en) | 2020-03-24 | 2021-09-30 | Fluent Biosciences Inc. | Viral detection using template emulsification |
| IL298821A (en) * | 2020-06-09 | 2023-02-01 | Illumina Inc | A method for increasing the throughput of sequencing libraries |
| CA3184751A1 (en) * | 2020-06-12 | 2021-12-16 | President And Fellows Of Harvard College | Compositions and methods for dna methylation analysis |
| EP4182477A4 (en) | 2020-07-15 | 2024-09-11 | Fluent Biosciences Inc. | MULTI-LEVEL LIGATION OLIGOS |
| EP4196600A4 (en) * | 2020-08-14 | 2024-09-11 | Factorial Diagnostics, Inc. | PREPARATION OF IN SITU LIBRARIES FOR SEQUENCING |
| CN112251422B (zh) * | 2020-10-21 | 2024-04-19 | 华中农业大学 | 含独特分子标签序列的转座酶复合体及其应用 |
| CN112309500B (zh) * | 2020-10-30 | 2024-08-30 | 广州序科码生物技术有限责任公司 | 一种基于单细胞测序数据唯一片段序列捕获方法 |
| CN114686564B (zh) * | 2020-12-31 | 2023-10-24 | 深圳华大生命科学研究院 | 一种适用于微量样本的dna样本制备方法 |
| WO2022141366A1 (zh) * | 2020-12-31 | 2022-07-07 | 北京寻因生物科技有限公司 | 含有移码碱基的细胞标签微珠及其制备方法和应用 |
| EP4027567A1 (en) | 2021-01-08 | 2022-07-13 | Comcast Cable Communications, LLC | Beam control for repetitions |
| MX2023009046A (es) * | 2021-02-04 | 2023-08-10 | Illumina Inc | Generacion de lectura unida indexada prolongada en globulos unidos con transposomas. |
| CA3211172A1 (en) * | 2021-03-31 | 2022-10-06 | Illumina, Inc. | Methods of preparing directional tagmentation sequencing libraries using transposon-based technology with unique molecular identifiers for error correction |
| US20220325277A1 (en) * | 2021-04-13 | 2022-10-13 | Inscripta, Inc. | Genomic edit detection at the single cell level |
| WO2023019024A2 (en) * | 2021-08-13 | 2023-02-16 | Arizona Board Of Regents On Behalf Of Arizona State University | A method for single-cell dna sequencing via in situ genomic amplification and combinatorial barcoding |
| WO2023079047A1 (en) * | 2021-11-04 | 2023-05-11 | Universal Diagnostics S.A | Systems and methods for preparing biological samples for genetic sequencing |
| US20240218355A1 (en) * | 2021-11-09 | 2024-07-04 | Chengdu Boe Optoelectronics Technology Co., Ltd. | Dna end repair reagent and kit thereof, dna library construction kit, and method for constructing dna library |
| CN114807305A (zh) * | 2022-04-13 | 2022-07-29 | 首都医科大学附属北京口腔医院 | 一种构建原核生物单细胞rna测序文库的方法 |
| CN114990146A (zh) * | 2022-06-30 | 2022-09-02 | 昆明理工大学 | 一种用于染色质分析的多任务工具包及其制备方法 |
| CN115064210B (zh) * | 2022-07-27 | 2022-11-18 | 北京大学第三医院(北京大学第三临床医学院) | 一种鉴定二倍体胚胎细胞中染色体交叉互换位置的方法及应用 |
| US12534721B2 (en) | 2022-10-06 | 2026-01-27 | Illumina, Inc. | Quantitative detection and analysis of moleculesip |
| CN115537408B (zh) * | 2022-10-08 | 2025-04-18 | 厦门大学 | 一种单细胞多组学文库及其构建方法 |
| CN118727162B (zh) * | 2024-06-05 | 2025-04-04 | 首都医科大学附属北京口腔医院 | 单细胞4c文库构建方法和检测方法 |
| EP4686761A1 (en) * | 2024-07-31 | 2026-02-04 | 4basebio UK Ltd | Linear dna products produced from multiple amplified dna products |
| CN120673842A (zh) * | 2025-05-08 | 2025-09-19 | 内蒙古农业大学 | 一种绒山羊绒厚性状的基因组选择育种的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017506877A (ja) | 2013-12-30 | 2017-03-16 | アトレカ インコーポレイテッド | 核酸バーコードを用いた単一細胞と関連づけられた核酸の分析方法 |
| US20180023119A1 (en) | 2016-07-22 | 2018-01-25 | Illumina, Inc. | Single cell whole genome libraries and combinatorial indexing methods of making thereof |
| JP2020523011A (ja) | 2017-06-07 | 2020-08-06 | オレゴン ヘルス アンド サイエンス ユニバーシティ | メチル化配列決定のための単一細胞全ゲノムライブラリー |
Family Cites Families (102)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4683195A (en) | 1986-01-30 | 1987-07-28 | Cetus Corporation | Process for amplifying, detecting, and/or-cloning nucleic acid sequences |
| US4683202A (en) | 1985-03-28 | 1987-07-28 | Cetus Corporation | Process for amplifying nucleic acid sequences |
| US4882268A (en) | 1985-12-24 | 1989-11-21 | Massachusetts Institute Of Technology | Method for determining tissue of origin and degree of malignancy of tumor cells |
| CA1323293C (en) | 1987-12-11 | 1993-10-19 | Keith C. Backman | Assay using template-dependent nucleic acid probe reorganization |
| CA1341584C (en) | 1988-04-06 | 2008-11-18 | Bruce Wallace | Method of amplifying and detecting nucleic acid sequences |
| AU3539089A (en) | 1988-04-08 | 1989-11-03 | Salk Institute For Biological Studies, The | Ligase-based amplification method |
| US5130238A (en) | 1988-06-24 | 1992-07-14 | Cangene Corporation | Enhanced nucleic acid amplification process |
| ATE144556T1 (de) | 1988-06-24 | 1996-11-15 | Amgen Inc | Verfahren und mittel zum nachweis von nukleinsäuresequenzen |
| ATE138106T1 (de) | 1988-07-20 | 1996-06-15 | David Segev | Verfahren zur amplifizierung und zum nachweis von nukleinsäuresequenzen |
| US5185243A (en) | 1988-08-25 | 1993-02-09 | Syntex (U.S.A.) Inc. | Method for detection of specific nucleic acid sequences |
| WO1991006678A1 (en) | 1989-10-26 | 1991-05-16 | Sri International | Dna sequencing |
| CA2035010C (en) | 1990-01-26 | 1996-12-10 | Keith C. Backman | Method of amplifying target nucleic acids applicable to both polymerase and ligase chain reactions |
| US5573907A (en) | 1990-01-26 | 1996-11-12 | Abbott Laboratories | Detecting and amplifying target nucleic acids using exonucleolytic activity |
| US5223414A (en) | 1990-05-07 | 1993-06-29 | Sri International | Process for nucleic acid hybridization and amplification |
| US5455166A (en) | 1991-01-31 | 1995-10-03 | Becton, Dickinson And Company | Strand displacement amplification |
| WO1995021271A1 (en) | 1994-02-07 | 1995-08-10 | Molecular Tool, Inc. | Ligase/polymerase-mediated genetic bit analysistm of single nucleotide polymorphisms and its use in genetic analysis |
| US5677170A (en) | 1994-03-02 | 1997-10-14 | The Johns Hopkins University | In vitro transposition of artificial transposons |
| WO1995025180A1 (en) | 1994-03-16 | 1995-09-21 | Gen-Probe Incorporated | Isothermal strand displacement nucleic acid amplification |
| US5846719A (en) | 1994-10-13 | 1998-12-08 | Lynx Therapeutics, Inc. | Oligonucleotide tags for sorting and identification |
| US5750341A (en) | 1995-04-17 | 1998-05-12 | Lynx Therapeutics, Inc. | DNA sequencing by parallel oligonucleotide extensions |
| GB9620209D0 (en) | 1996-09-27 | 1996-11-13 | Cemu Bioteknik Ab | Method of sequencing DNA |
| GB9626815D0 (en) | 1996-12-23 | 1997-02-12 | Cemu Bioteknik Ab | Method of sequencing DNA |
| JP2002503954A (ja) | 1997-04-01 | 2002-02-05 | グラクソ、グループ、リミテッド | 核酸増幅法 |
| US6969488B2 (en) | 1998-05-22 | 2005-11-29 | Solexa, Inc. | System and apparatus for sequential processing of analytes |
| US6245974B1 (en) | 1997-08-06 | 2001-06-12 | North Carolina State University | Matrix attachment regions |
| AR021833A1 (es) | 1998-09-30 | 2002-08-07 | Applied Research Systems | Metodos de amplificacion y secuenciacion de acido nucleico |
| DE19849348A1 (de) * | 1998-10-26 | 2000-04-27 | Univ Ludwigs Albert | Linear Amplification mediated PCR (=LAM PCR) |
| US6274320B1 (en) | 1999-09-16 | 2001-08-14 | Curagen Corporation | Method of sequencing a nucleic acid |
| US7955794B2 (en) | 2000-09-21 | 2011-06-07 | Illumina, Inc. | Multiplex nucleic acid reactions |
| US7582420B2 (en) | 2001-07-12 | 2009-09-01 | Illumina, Inc. | Multiplex nucleic acid reactions |
| US7611869B2 (en) | 2000-02-07 | 2009-11-03 | Illumina, Inc. | Multiplexed methylation detection methods |
| US7001792B2 (en) | 2000-04-24 | 2006-02-21 | Eagle Research & Development, Llc | Ultra-fast nucleic acid sequencing device and a method for making and using the same |
| ATE377093T1 (de) | 2000-07-07 | 2007-11-15 | Visigen Biotechnologies Inc | Sequenzbestimmung in echtzeit |
| AU2002227156A1 (en) | 2000-12-01 | 2002-06-11 | Visigen Biotechnologies, Inc. | Enzymatic nucleic acid synthesis: compositions and methods for altering monomer incorporation fidelity |
| AR031640A1 (es) | 2000-12-08 | 2003-09-24 | Applied Research Systems | Amplificacion isotermica de acidos nucleicos en un soporte solido |
| US7057026B2 (en) | 2001-12-04 | 2006-06-06 | Solexa Limited | Labelled nucleotides |
| US7399590B2 (en) | 2002-02-21 | 2008-07-15 | Asm Scientific, Inc. | Recombinase polymerase amplification |
| US8030000B2 (en) | 2002-02-21 | 2011-10-04 | Alere San Diego, Inc. | Recombinase polymerase amplification |
| WO2004018497A2 (en) | 2002-08-23 | 2004-03-04 | Solexa Limited | Modified nucleotides for polynucleotide sequencing |
| CA2498764C (en) | 2002-09-20 | 2015-11-10 | New England Biolabs, Inc. | Helicase dependent amplification of nucleic acids |
| CA2527947A1 (en) | 2003-06-06 | 2004-12-16 | The Board Of Trustees Of The University Of Illinois | Methods for assessing the invasive potential of a cell employing chromatin analysis |
| WO2005003304A2 (en) | 2003-06-20 | 2005-01-13 | Illumina, Inc. | Methods and compositions for whole genome amplification and genotyping |
| GB0321306D0 (en) | 2003-09-11 | 2003-10-15 | Solexa Ltd | Modified polymerases for improved incorporation of nucleotide analogues |
| EP1701785A1 (en) | 2004-01-07 | 2006-09-20 | Solexa Ltd. | Modified molecular arrays |
| CA2579150C (en) | 2004-09-17 | 2014-11-25 | Pacific Biosciences Of California, Inc. | Apparatus and method for analysis of molecules |
| WO2006064199A1 (en) | 2004-12-13 | 2006-06-22 | Solexa Limited | Improved method of nucleotide detection |
| EP1888743B1 (en) | 2005-05-10 | 2011-08-03 | Illumina Cambridge Limited | Improved polymerases |
| EP3257949A1 (en) | 2005-06-15 | 2017-12-20 | Complete Genomics Inc. | Nucleic acid analysis by random mixtures of non-overlapping fragments |
| GB0514936D0 (en) | 2005-07-20 | 2005-08-24 | Solexa Ltd | Preparation of templates for nucleic acid sequencing |
| US7405281B2 (en) | 2005-09-29 | 2008-07-29 | Pacific Biosciences Of California, Inc. | Fluorescent nucleotide analogs and uses therefor |
| GB0522310D0 (en) | 2005-11-01 | 2005-12-07 | Solexa Ltd | Methods of preparing libraries of template polynucleotides |
| CN101432439B (zh) | 2006-02-24 | 2013-07-24 | 考利达基因组股份有限公司 | Dna阵列上的高通量基因组测序 |
| SG170028A1 (en) | 2006-02-24 | 2011-04-29 | Callida Genomics Inc | High throughput genome sequencing on dna arrays |
| WO2007107710A1 (en) | 2006-03-17 | 2007-09-27 | Solexa Limited | Isothermal methods for creating clonal single molecule arrays |
| CA2648149A1 (en) | 2006-03-31 | 2007-11-01 | Solexa, Inc. | Systems and devices for sequence by synthesis analysis |
| US8343746B2 (en) | 2006-10-23 | 2013-01-01 | Pacific Biosciences Of California, Inc. | Polymerase enzymes and reagents for enhanced nucleic acid sequencing |
| US7910354B2 (en) | 2006-10-27 | 2011-03-22 | Complete Genomics, Inc. | Efficient arrays of amplified polynucleotides |
| US8349167B2 (en) | 2006-12-14 | 2013-01-08 | Life Technologies Corporation | Methods and apparatus for detecting molecular interactions using FET arrays |
| US8262900B2 (en) | 2006-12-14 | 2012-09-11 | Life Technologies Corporation | Methods and apparatus for measuring analytes using large scale FET arrays |
| GB2457851B (en) | 2006-12-14 | 2011-01-05 | Ion Torrent Systems Inc | Methods and apparatus for measuring analytes using large scale fet arrays |
| WO2008093098A2 (en) | 2007-02-02 | 2008-08-07 | Illumina Cambridge Limited | Methods for indexing samples and sequencing multiple nucleotide templates |
| CN101802223A (zh) | 2007-08-15 | 2010-08-11 | 香港大学 | 用于高通量亚硫酸氢盐dna-测序的方法和组合物及其用途 |
| WO2010003132A1 (en) | 2008-07-02 | 2010-01-07 | Illumina Cambridge Ltd. | Using populations of beads for the fabrication of arrays on surfaces |
| US20100137143A1 (en) | 2008-10-22 | 2010-06-03 | Ion Torrent Systems Incorporated | Methods and apparatus for measuring analytes |
| US9080211B2 (en) | 2008-10-24 | 2015-07-14 | Epicentre Technologies Corporation | Transposon end compositions and methods for modifying nucleic acids |
| EP2370599B1 (en) | 2010-01-19 | 2015-01-21 | Verinata Health, Inc | Method for determining copy number variations |
| WO2012061832A1 (en) | 2010-11-05 | 2012-05-10 | Illumina, Inc. | Linking sequence reads using paired code tags |
| US8829171B2 (en) | 2011-02-10 | 2014-09-09 | Illumina, Inc. | Linking sequence reads using paired code tags |
| US9074251B2 (en) | 2011-02-10 | 2015-07-07 | Illumina, Inc. | Linking sequence reads using paired code tags |
| US8951781B2 (en) | 2011-01-10 | 2015-02-10 | Illumina, Inc. | Systems, methods, and apparatuses to image a sample for biological or chemical analysis |
| EP2714938B1 (en) | 2011-05-27 | 2017-11-15 | President and Fellows of Harvard College | Methods of amplifying whole genome of a single cell |
| US9926552B2 (en) | 2011-06-06 | 2018-03-27 | Cornell University | Microfluidic device for extracting, isolating, and analyzing DNA from cells |
| EP2718465B1 (en) | 2011-06-09 | 2022-04-13 | Illumina, Inc. | Method of making an analyte array |
| CA2859660C (en) | 2011-09-23 | 2021-02-09 | Illumina, Inc. | Methods and compositions for nucleic acid sequencing |
| CN103060924B (zh) * | 2011-10-18 | 2016-04-20 | 深圳华大基因科技有限公司 | 微量核酸样本的文库制备方法及其应用 |
| EP3305400A3 (en) | 2011-10-28 | 2018-06-06 | Illumina, Inc. | Microarray fabrication system and method |
| EP3366348B1 (en) | 2012-01-16 | 2023-08-23 | Greatbatch Ltd. | Emi filtered co-connected hermetic feedthrough, feedthrough capacitor and leadwire assembly for an active implantable medical device |
| CN204832037U (zh) | 2012-04-03 | 2015-12-02 | 伊鲁米那股份有限公司 | 检测设备 |
| US8895249B2 (en) | 2012-06-15 | 2014-11-25 | Illumina, Inc. | Kinetic exclusion amplification of nucleic acid libraries |
| RU2529784C2 (ru) * | 2012-12-26 | 2014-09-27 | Закрытое акционерное общество "Геноаналитика" | Технология определения анеуплоидии методом секвенирования |
| US9683230B2 (en) | 2013-01-09 | 2017-06-20 | Illumina Cambridge Limited | Sample preparation on a solid support |
| US9512422B2 (en) | 2013-02-26 | 2016-12-06 | Illumina, Inc. | Gel patterned surfaces |
| DK2970951T3 (da) | 2013-03-13 | 2019-05-13 | Illumina Inc | Fremgangsmåder til nukleinsyresekventering |
| US9328382B2 (en) | 2013-03-15 | 2016-05-03 | Complete Genomics, Inc. | Multiple tagging of individual long DNA fragments |
| WO2014145078A1 (en) | 2013-03-15 | 2014-09-18 | Verinata Health, Inc. | Generating cell-free dna libraries directly from blood |
| SG11201508985VA (en) | 2013-05-23 | 2015-12-30 | Univ Leland Stanford Junior | Transposition into native chromatin for personal epigenomics |
| DK3017065T3 (en) | 2013-07-01 | 2018-11-26 | Illumina Inc | Catalyst-free Surface functionalization and polymer grafting |
| ES2764096T3 (es) | 2013-08-19 | 2020-06-02 | Abbott Molecular Inc | Bibliotecas de secuenciación de próxima generación |
| CA2936564C (en) | 2014-01-07 | 2022-10-18 | Fundacio Privada Institut De Medicina Predictiva I Personalitzada Del Cancer | Methods for generating double stranded dna libraries and sequencing methods for the identification of methylated cytosines |
| US9677132B2 (en) | 2014-01-16 | 2017-06-13 | Illumina, Inc. | Polynucleotide modification on solid support |
| RU2571855C2 (ru) | 2014-05-12 | 2015-12-20 | Государственное научное учреждение Всероссийский научно-исследовательский институт ветеринарной вирусологии и микробиологии | Набор олигонуклеотидных праймеров для получения библиотеки полноразмерных генов, кодирующих иммунодоминантные белки р17 (d11l), р30 (cp204l), р54 (e183l), р60 (cp530r), р72 (b646l) вируса африканской чумы свиней |
| CN112430641A (zh) | 2014-06-30 | 2021-03-02 | 亿明达股份有限公司 | 使用单侧转座的方法和组合物 |
| US20160061517A1 (en) | 2014-08-29 | 2016-03-03 | Black & Veatch Holding Company | Dual mixed refrigerant system |
| RU2736728C2 (ru) | 2014-10-17 | 2020-11-19 | Иллумина Кембридж Лимитед | Транспозиция с сохранением сцепления генов |
| WO2016066586A1 (en) | 2014-10-31 | 2016-05-06 | Illumina Cambridge Limited | Novel polymers and dna copolymer coatings |
| US10900065B2 (en) | 2014-11-14 | 2021-01-26 | University Of Washington | Methods and kits for labeling cellular molecules |
| GB2532790B (en) | 2014-11-28 | 2017-10-04 | Porvair Filtration Group Ltd | Modified pipette tips for chromatin immunoprecipitation assay |
| CN107406890B (zh) * | 2015-02-10 | 2023-07-18 | 亿明达股份有限公司 | 用于分析细胞组分的方法和组合物 |
| WO2017156336A1 (en) | 2016-03-10 | 2017-09-14 | The Board Of Trustees Of The Leland Stanford Junior University | Transposase-mediated imaging of the accessible genome |
| US10240196B2 (en) * | 2016-05-27 | 2019-03-26 | Agilent Technologies, Inc. | Transposase-random priming DNA sample preparation |
| US10400235B2 (en) * | 2017-05-26 | 2019-09-03 | 10X Genomics, Inc. | Single cell analysis of transposase accessible chromatin |
| IL271454B2 (en) | 2018-05-17 | 2025-04-01 | Illumina Inc | Rapid single-cell genetic sequencing with low Aggregation bias |
-
2019
- 2019-05-17 IL IL271454A patent/IL271454B2/en unknown
- 2019-05-17 JP JP2019569897A patent/JP7234146B2/ja active Active
- 2019-05-17 US US16/415,969 patent/US11981891B2/en active Active
- 2019-05-17 AU AU2019270185A patent/AU2019270185A1/en not_active Abandoned
- 2019-05-17 EP EP19803332.6A patent/EP3794141A4/en active Pending
- 2019-05-17 KR KR1020207000938A patent/KR102447811B1/ko active Active
- 2019-05-17 SG SG11202001827UA patent/SG11202001827UA/en unknown
- 2019-05-17 BR BR112020005982-0A patent/BR112020005982A2/pt unknown
- 2019-05-17 CA CA3206377A patent/CA3206377A1/en active Pending
- 2019-05-17 CA CA3067435A patent/CA3067435C/en active Active
- 2019-05-17 KR KR1020227032955A patent/KR102640255B1/ko active Active
- 2019-05-17 WO PCT/US2019/032966 patent/WO2019222688A1/en not_active Ceased
- 2019-05-17 MX MX2019015013A patent/MX2019015013A/es unknown
- 2019-05-17 RU RU2019144286A patent/RU2744175C1/ru active
- 2019-05-17 CN CN201980003422.5A patent/CN111201329B/zh active Active
- 2019-05-17 CN CN202510085190.XA patent/CN119955902A/zh active Pending
- 2019-12-11 MX MX2023010271A patent/MX2023010271A/es unknown
-
2022
- 2022-04-26 AU AU2022202739A patent/AU2022202739B2/en active Active
-
2023
- 2023-02-22 JP JP2023025703A patent/JP7638309B2/ja active Active
-
2024
- 2024-05-03 US US18/654,545 patent/US20240376460A1/en active Pending
- 2024-09-20 AU AU2024219947A patent/AU2024219947A1/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017506877A (ja) | 2013-12-30 | 2017-03-16 | アトレカ インコーポレイテッド | 核酸バーコードを用いた単一細胞と関連づけられた核酸の分析方法 |
| US20180023119A1 (en) | 2016-07-22 | 2018-01-25 | Illumina, Inc. | Single cell whole genome libraries and combinatorial indexing methods of making thereof |
| JP2020523011A (ja) | 2017-06-07 | 2020-08-06 | オレゴン ヘルス アンド サイエンス ユニバーシティ | メチル化配列決定のための単一細胞全ゲノムライブラリー |
Non-Patent Citations (4)
| Title |
|---|
| CHEN C. et al.,Science,356 (2017),p.189-194 |
| CUSANOVICH D. A. et al.,Science,Vol.348 Issue 6237 (2015),p.910-914 |
| VITAK S. A. et al.,Nature Methods,Vol.14 No.3 (2017),p.302-308, ONLINE METHODS |
| YIN Y. et al.,Molecular Cell,76 (2019 Nov),p.676-690 |
Also Published As
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7638309B2 (ja) | 低減した増幅バイアスによるハイスループット単一細胞シークエンシング | |
| JP7564158B2 (ja) | 単細胞全ゲノムライブラリおよびそれを作成する組み合わせインデックス付加方法 | |
| CN110997932B (zh) | 用于甲基化测序的单细胞全基因组文库 | |
| JP2025108660A (ja) | ハイスループット単一細胞ライブラリー、並びに製造方法及び使用方法 | |
| RU2833615C2 (ru) | Высокопроизводительное секвенирование одиночной клетки со сниженной ошибкой амплификации | |
| HK40031087B (zh) | 具有减少的扩增偏倚的高通量单细胞测序 | |
| HK40111421A (zh) | 用於甲基化测序的单细胞全基因组文库 | |
| HK40031087A (en) | High-throughput single-cell sequencing with reduced amplification bias | |
| HK40027672B (zh) | 用於甲基化测序的单细胞全基因组文库 | |
| NZ749719B2 (en) | Single cell whole genome libraries and combinatorial indexing methods of making thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230323 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230323 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20240502 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240507 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20240806 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20241002 |
|
| 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: 20250121 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20250218 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7638309 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |