US20150167080A1 - Method of determining clonotypes and clonotype profiles - Google Patents
Method of determining clonotypes and clonotype profiles Download PDFInfo
- Publication number
- US20150167080A1 US20150167080A1 US13/196,885 US201113196885A US2015167080A1 US 20150167080 A1 US20150167080 A1 US 20150167080A1 US 201113196885 A US201113196885 A US 201113196885A US 2015167080 A1 US2015167080 A1 US 2015167080A1
- Authority
- US
- United States
- Prior art keywords
- sequence
- clonotype
- frequency candidate
- candidate clonotype
- clonotypes
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 124
- 238000012163 sequencing technique Methods 0.000 claims abstract description 93
- 150000007523 nucleic acids Chemical class 0.000 claims abstract description 80
- 102000039446 nucleic acids Human genes 0.000 claims abstract description 77
- 108020004707 nucleic acids Proteins 0.000 claims abstract description 77
- 238000004581 coalescence Methods 0.000 claims abstract description 12
- 239000002773 nucleotide Substances 0.000 claims description 56
- 125000003729 nucleotide group Chemical group 0.000 claims description 56
- 108091008874 T cell receptors Proteins 0.000 claims description 46
- 102000016266 T-Cell Antigen Receptors Human genes 0.000 claims description 45
- 210000003719 b-lymphocyte Anatomy 0.000 claims description 37
- 210000001744 T-lymphocyte Anatomy 0.000 claims description 34
- 210000004698 lymphocyte Anatomy 0.000 claims description 27
- 230000007423 decrease Effects 0.000 claims description 18
- 239000007787 solid Substances 0.000 claims description 18
- 102000018358 immunoglobulin Human genes 0.000 claims description 12
- 108060003951 Immunoglobulin Proteins 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 210000004369 blood Anatomy 0.000 claims description 11
- 239000008280 blood Substances 0.000 claims description 11
- 238000003786 synthesis reaction Methods 0.000 claims description 8
- 238000012217 deletion Methods 0.000 claims description 7
- 230000037430 deletion Effects 0.000 claims description 7
- 238000003780 insertion Methods 0.000 claims description 5
- 230000037431 insertion Effects 0.000 claims description 5
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 abstract description 18
- 238000005259 measurement Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 8
- 238000003556 assay Methods 0.000 abstract description 2
- 230000004069 differentiation Effects 0.000 abstract 1
- 230000003321 amplification Effects 0.000 description 85
- 238000003199 nucleic acid amplification method Methods 0.000 description 85
- 239000000523 sample Substances 0.000 description 75
- 238000006243 chemical reaction Methods 0.000 description 44
- 108091093088 Amplicon Proteins 0.000 description 43
- 102000006496 Immunoglobulin Heavy Chains Human genes 0.000 description 36
- 108010019476 Immunoglobulin Heavy Chains Proteins 0.000 description 36
- 210000004027 cell Anatomy 0.000 description 30
- 108020004414 DNA Proteins 0.000 description 24
- 238000003753 real-time PCR Methods 0.000 description 21
- 108091032973 (ribonucleotides)n+m Proteins 0.000 description 19
- 230000000875 corresponding effect Effects 0.000 description 19
- 108091028043 Nucleic acid sequence Proteins 0.000 description 18
- 230000000295 complement effect Effects 0.000 description 18
- 108091008875 B cell receptors Proteins 0.000 description 17
- 230000006870 function Effects 0.000 description 17
- 230000002441 reversible effect Effects 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 14
- 229920003266 Leaf® Polymers 0.000 description 13
- 230000035772 mutation Effects 0.000 description 12
- 108020004705 Codon Proteins 0.000 description 11
- 108091008915 immune receptors Proteins 0.000 description 10
- 102000027596 immune receptors Human genes 0.000 description 10
- 238000005215 recombination Methods 0.000 description 10
- 230000006798 recombination Effects 0.000 description 10
- 238000011160 research Methods 0.000 description 10
- 206010069754 Acquired gene mutation Diseases 0.000 description 9
- 230000037439 somatic mutation Effects 0.000 description 9
- 206010028980 Neoplasm Diseases 0.000 description 8
- 108091033319 polynucleotide Proteins 0.000 description 8
- 102000040430 polynucleotide Human genes 0.000 description 8
- 239000002157 polynucleotide Substances 0.000 description 8
- 239000000047 product Substances 0.000 description 8
- 108090000623 proteins and genes Proteins 0.000 description 8
- 108010047041 Complementarity Determining Regions Proteins 0.000 description 7
- 238000000137 annealing Methods 0.000 description 7
- 238000007796 conventional method Methods 0.000 description 7
- 238000003752 polymerase chain reaction Methods 0.000 description 7
- 239000002299 complementary DNA Substances 0.000 description 6
- 201000010099 disease Diseases 0.000 description 6
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 6
- 239000000463 material Substances 0.000 description 6
- 230000000392 somatic effect Effects 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- 108700028369 Alleles Proteins 0.000 description 5
- 206010025323 Lymphomas Diseases 0.000 description 5
- 238000003491 array Methods 0.000 description 5
- 230000037429 base substitution Effects 0.000 description 5
- 230000003247 decreasing effect Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 238000010790 dilution Methods 0.000 description 5
- 239000012895 dilution Substances 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 208000019420 lymphoid neoplasm Diseases 0.000 description 5
- 230000002829 reductive effect Effects 0.000 description 5
- -1 IgH molecules Chemical class 0.000 description 4
- 108091034117 Oligonucleotide Proteins 0.000 description 4
- 239000000427 antigen Substances 0.000 description 4
- 108091007433 antigens Proteins 0.000 description 4
- 102000036639 antigens Human genes 0.000 description 4
- 201000011510 cancer Diseases 0.000 description 4
- 239000007795 chemical reaction product Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000009396 hybridization Methods 0.000 description 4
- 208000032839 leukemia Diseases 0.000 description 4
- 238000007403 mPCR Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 108010014303 DNA-directed DNA polymerase Proteins 0.000 description 3
- 102000016928 DNA-directed DNA polymerase Human genes 0.000 description 3
- 108010006785 Taq Polymerase Proteins 0.000 description 3
- 238000007792 addition Methods 0.000 description 3
- 239000011324 bead Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 239000000839 emulsion Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000012634 fragment Substances 0.000 description 3
- 229920001519 homopolymer Polymers 0.000 description 3
- 210000002865 immune cell Anatomy 0.000 description 3
- 108020004999 messenger RNA Proteins 0.000 description 3
- 210000005259 peripheral blood Anatomy 0.000 description 3
- 239000011886 peripheral blood Substances 0.000 description 3
- 238000003757 reverse transcription PCR Methods 0.000 description 3
- 238000005096 rolling process Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 208000024893 Acute lymphoblastic leukemia Diseases 0.000 description 2
- 208000014697 Acute lymphocytic leukaemia Diseases 0.000 description 2
- 238000001712 DNA sequencing Methods 0.000 description 2
- 102000013463 Immunoglobulin Light Chains Human genes 0.000 description 2
- 108010065825 Immunoglobulin Light Chains Proteins 0.000 description 2
- 208000031671 Large B-Cell Diffuse Lymphoma Diseases 0.000 description 2
- 102000003960 Ligases Human genes 0.000 description 2
- 108090000364 Ligases Proteins 0.000 description 2
- TWRXJAOTZQYOKJ-UHFFFAOYSA-L Magnesium chloride Chemical compound [Mg+2].[Cl-].[Cl-] TWRXJAOTZQYOKJ-UHFFFAOYSA-L 0.000 description 2
- 208000031481 Pathologic Constriction Diseases 0.000 description 2
- 208000006664 Precursor Cell Lymphoblastic Leukemia-Lymphoma Diseases 0.000 description 2
- 238000012300 Sequence Analysis Methods 0.000 description 2
- 108700042075 T-Cell Receptor Genes Proteins 0.000 description 2
- JLCPHMBAVCMARE-UHFFFAOYSA-N [3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[5-(2-amino-6-oxo-1H-purin-9-yl)-3-[[3-[[3-[[3-[[3-[[3-[[5-(2-amino-6-oxo-1H-purin-9-yl)-3-[[5-(2-amino-6-oxo-1H-purin-9-yl)-3-hydroxyoxolan-2-yl]methoxy-hydroxyphosphoryl]oxyoxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxyoxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methyl [5-(6-aminopurin-9-yl)-2-(hydroxymethyl)oxolan-3-yl] hydrogen phosphate Polymers Cc1cn(C2CC(OP(O)(=O)OCC3OC(CC3OP(O)(=O)OCC3OC(CC3O)n3cnc4c3nc(N)[nH]c4=O)n3cnc4c3nc(N)[nH]c4=O)C(COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3CO)n3cnc4c(N)ncnc34)n3ccc(N)nc3=O)n3cnc4c(N)ncnc34)n3ccc(N)nc3=O)n3ccc(N)nc3=O)n3ccc(N)nc3=O)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)n3cc(C)c(=O)[nH]c3=O)n3cc(C)c(=O)[nH]c3=O)n3ccc(N)nc3=O)n3cc(C)c(=O)[nH]c3=O)n3cnc4c3nc(N)[nH]c4=O)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)O2)c(=O)[nH]c1=O JLCPHMBAVCMARE-UHFFFAOYSA-N 0.000 description 2
- 230000001594 aberrant effect Effects 0.000 description 2
- 239000002253 acid Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010804 cDNA synthesis Methods 0.000 description 2
- 230000009260 cross reactivity Effects 0.000 description 2
- 230000001351 cycling effect Effects 0.000 description 2
- 210000001151 cytotoxic T lymphocyte Anatomy 0.000 description 2
- 238000004925 denaturation Methods 0.000 description 2
- 230000036425 denaturation Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 206010012818 diffuse large B-cell lymphoma Diseases 0.000 description 2
- 238000001962 electrophoresis Methods 0.000 description 2
- 238000004128 high performance liquid chromatography Methods 0.000 description 2
- 238000012165 high-throughput sequencing Methods 0.000 description 2
- 229940072221 immunoglobulins Drugs 0.000 description 2
- 230000002401 inhibitory effect Effects 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 210000000265 leukocyte Anatomy 0.000 description 2
- 235000019689 luncheon sausage Nutrition 0.000 description 2
- 230000003211 malignant effect Effects 0.000 description 2
- 238000010369 molecular cloning Methods 0.000 description 2
- 238000007857 nested PCR Methods 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 239000000376 reactant Substances 0.000 description 2
- 239000011541 reaction mixture Substances 0.000 description 2
- 230000008707 rearrangement Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000000153 supplemental effect Effects 0.000 description 2
- 102000007469 Actins Human genes 0.000 description 1
- 108010085238 Actins Proteins 0.000 description 1
- 241001156002 Anthonomus pomorum Species 0.000 description 1
- 101100519158 Arabidopsis thaliana PCR2 gene Proteins 0.000 description 1
- 241000417436 Arcotheres Species 0.000 description 1
- 208000010839 B-cell chronic lymphocytic leukemia Diseases 0.000 description 1
- 241000894006 Bacteria Species 0.000 description 1
- 210000001266 CD8-positive T-lymphocyte Anatomy 0.000 description 1
- 206010065163 Clonal evolution Diseases 0.000 description 1
- 102000053602 DNA Human genes 0.000 description 1
- 102000012410 DNA Ligases Human genes 0.000 description 1
- 108010061982 DNA Ligases Proteins 0.000 description 1
- 230000004543 DNA replication Effects 0.000 description 1
- 108090000626 DNA-directed RNA polymerases Proteins 0.000 description 1
- 102000004163 DNA-directed RNA polymerases Human genes 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 102100031181 Glyceraldehyde-3-phosphate dehydrogenase Human genes 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 108700005091 Immunoglobulin Genes Proteins 0.000 description 1
- 102000012745 Immunoglobulin Subunits Human genes 0.000 description 1
- 108010079585 Immunoglobulin Subunits Proteins 0.000 description 1
- 208000030289 Lymphoproliferative disease Diseases 0.000 description 1
- 206010061240 Malignant lymphoid neoplasm Diseases 0.000 description 1
- 108700011259 MicroRNAs Proteins 0.000 description 1
- 208000034578 Multiple myelomas Diseases 0.000 description 1
- 102100030569 Nuclear receptor corepressor 2 Human genes 0.000 description 1
- 101710153660 Nuclear receptor corepressor 2 Proteins 0.000 description 1
- 108020005187 Oligonucleotide Probes Proteins 0.000 description 1
- 238000012408 PCR amplification Methods 0.000 description 1
- ZYFVNVRFVHJEIU-UHFFFAOYSA-N PicoGreen Chemical compound CN(C)CCCN(CCCN(C)C)C1=CC(=CC2=[N+](C3=CC=CC=C3S2)C)C2=CC=CC=C2N1C1=CC=CC=C1 ZYFVNVRFVHJEIU-UHFFFAOYSA-N 0.000 description 1
- 206010035226 Plasma cell myeloma Diseases 0.000 description 1
- 108091081021 Sense strand Proteins 0.000 description 1
- 108020004682 Single-Stranded DNA Proteins 0.000 description 1
- 108700042077 T-Cell Receptor beta Genes Proteins 0.000 description 1
- 206010042971 T-cell lymphoma Diseases 0.000 description 1
- 208000027585 T-cell non-Hodgkin lymphoma Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 150000007513 acids Chemical class 0.000 description 1
- 230000004721 adaptive immunity Effects 0.000 description 1
- 102000006707 alpha-beta T-Cell Antigen Receptors Human genes 0.000 description 1
- 108010087408 alpha-beta T-Cell Antigen Receptors Proteins 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 238000004873 anchoring Methods 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 210000001185 bone marrow Anatomy 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 239000006172 buffering agent Substances 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000010367 cloning Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010511 deprotection reaction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000000975 dye Substances 0.000 description 1
- 210000003162 effector t lymphocyte Anatomy 0.000 description 1
- 230000002255 enzymatic effect Effects 0.000 description 1
- 230000003325 follicular Effects 0.000 description 1
- 201000003444 follicular lymphoma Diseases 0.000 description 1
- 238000001502 gel electrophoresis Methods 0.000 description 1
- 238000013412 genome amplification Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 108020004445 glyceraldehyde-3-phosphate dehydrogenase Proteins 0.000 description 1
- XDDAORKBJWWYJS-UHFFFAOYSA-N glyphosate Chemical compound OC(=O)CNCP(O)(O)=O XDDAORKBJWWYJS-UHFFFAOYSA-N 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000003394 haemopoietic effect Effects 0.000 description 1
- 201000009277 hairy cell leukemia Diseases 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000002443 helper t lymphocyte Anatomy 0.000 description 1
- 239000000833 heterodimer Substances 0.000 description 1
- 230000003100 immobilizing effect Effects 0.000 description 1
- 210000000987 immune system Anatomy 0.000 description 1
- 230000002998 immunogenetic effect Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000009830 intercalation Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 210000001165 lymph node Anatomy 0.000 description 1
- 210000003563 lymphoid tissue Anatomy 0.000 description 1
- 229910001629 magnesium chloride Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 210000003826 marginal zone b cell Anatomy 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 210000001806 memory b lymphocyte Anatomy 0.000 description 1
- 210000003071 memory t lymphocyte Anatomy 0.000 description 1
- 239000000693 micelle Substances 0.000 description 1
- 239000002679 microRNA Substances 0.000 description 1
- 238000002493 microarray Methods 0.000 description 1
- 238000010899 nucleation Methods 0.000 description 1
- 238000001668 nucleic acid synthesis Methods 0.000 description 1
- 239000002777 nucleoside Substances 0.000 description 1
- 230000005257 nucleotidylation Effects 0.000 description 1
- 239000002751 oligonucleotide probe Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 229920001184 polypeptide Polymers 0.000 description 1
- 208000017805 post-transplant lymphoproliferative disease Diseases 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 230000037452 priming Effects 0.000 description 1
- 108090000765 processed proteins & peptides Proteins 0.000 description 1
- 102000004196 processed proteins & peptides Human genes 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000012175 pyrosequencing Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000011867 re-evaluation Methods 0.000 description 1
- 239000011535 reaction buffer Substances 0.000 description 1
- 102000005962 receptors Human genes 0.000 description 1
- 108020003175 receptors Proteins 0.000 description 1
- 238000010188 recombinant method Methods 0.000 description 1
- 210000003289 regulatory T cell Anatomy 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010839 reverse transcription Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 108020004418 ribosomal RNA Proteins 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000002864 sequence alignment Methods 0.000 description 1
- 238000011451 sequencing strategy Methods 0.000 description 1
- 239000007790 solid phase Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000010186 staining Methods 0.000 description 1
- 210000000130 stem cell Anatomy 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 239000001226 triphosphate Substances 0.000 description 1
- 235000011178 triphosphate Nutrition 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 238000010626 work up procedure Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/10—Sequence alignment; Homology search
-
- 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/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6881—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for tissue or cell typing, e.g. human leukocyte antigen [HLA] probes
-
- 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/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B10/00—ICT specially adapted for evolutionary bioinformatics, e.g. phylogenetic tree construction or analysis
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino 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
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/158—Expression markers
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/16—Primer sets for multiplex assays
Definitions
- nucleic acid sequences of large and complex populations of DNA and/or RNA e.g. Gloor et al, PLoS ONE 5(10): el 5406 (2010); Petrosino et al, Clinical Chemistry, 55(5): 856-866 (2009); Arstila et al, Science, 286: 958-961 (1999).
- profiles of nucleic acids encoding immune molecules, such as T cell or B cell receptors, or their components contain a wealth of information on the state of health or disease of an organism, so that the use of such profiles as diagnostic or prognostic indicators has been proposed for a wide variety of conditions, e.g. Faham and Willis, U.S.
- the present invention is drawn to methods for producing sequence-based profiles of complex nucleic acid populations.
- the invention is exemplified in a number of implementations and applications, some of which are summarized below and throughout the specification.
- the invention is directed to a method of forming a clonotype profile comprising the steps of (a) forming a data structure of somatically recombined immune molecules from sequence reads thereof, (b) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype, (c) removing the coalesced candidate clonotype from the data structure, and (d) repeating steps (b) and (c) until a clonotype profile is formed.
- such data structure is a sequence tree.
- the invention is directed to a method of generating a clonotype profile from an individual comprising the steps of: (a) spatially isolating individual molecules derived from a nucleic acid sample from T-cells and/or B-cells of the individual; (b) sequencing said spatially isolated individual molecules to produce a plurality of sequence reads each having portions old V region, an NDN region and a J region wherein for each such molecule at least one sequence lead encompasses a portion of the NDN region (c) forming from sequence reads encompassing at least a portion of an NDN region a data structure having elements representing candidate clonotypes, each element and its corresponding candidate clonotype having a frequency; (d) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype having a sequence of the highest frequency candidate clo
- the invention provides a method for determining clonotypes of an immune repertoire comprising the steps: (a) providing a set of sequence reads from a repertoire of recombined immune molecules each having a V region, an NDN region and a J region wherein for each such molecule at least one sequence read encompasses at least a portion of the NDN region of such molecule; (b) forming from sequence reads encompassing at least a portion of an NDN region a data structure having elements representing candidate clonotypes, each leaf and its corresponding candidate clonotype having a frequency; (c) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype having a sequence of the highest frequency candidate clonotype; (d) removing elements corresponding to the coalesced candidate clonotypes from the data structure; and
- the invention is directed to a method of determining clonotypes of an immune repertoire comprising the steps: (a) providing a set of sequence reads from a repertoire of recombined immune molecules each having portions of a V region, an NDN region and a J region wherein for each such molecule at least one sequence read encompasses a portion of the NDN region of such molecule; (b) forming from sequence reads encompassing portions of NDN regions a sequence tree having leaves representing candidate clonotypes, each leaf and its corresponding candidate clonotype having a frequency; (c) selecting a highest frequency candidate clonotype and identifying all said lesser frequency candidate clonotypes having a sequence difference therewith less than a predetermined difference value to form a coalescence subset; (d) coalescing with the highest frequency candidate clonotype any lesser frequency candidate clonotypes in the coalescence subset whenever such lesser frequency is below a predetermined frequency value to form a clonotype having
- the present invention provides methods for determining clonotypes and clonotype profiles from large sets of sequence data obtained by high throughput sequencing Of somatically recombined immune molecules.
- the invention implements the above methods by using one or more sequence trees for efficiently carrying out sequence comparisons necessary for distinguishing candidate clonotypes with genuine sequence differences from those with experimental or measurement errors, such as sequencing errors.
- FIG. 1A illustrates an IgH transcript and sources of natural variability within it.
- FIG. 1B illustrates concepts of clonotypes in sequence space and distances between closely related clonotypes.
- FIG. 1C is a flow chart illustrating one embodiment of a method for distinguishing genuinely different clonotypes from clonotypes that differ solely by sequencing errors (which should be coalesced).
- FIG. 1D illustrates the form of a numerical function used in one embodiment for determining whether or not to coalesce related clonotypes.
- FIGS. 1E-1F illustrate the use of sequence trees in the methods of the invention.
- FIGS. 2A-2B show a two-staged PCR scheme for amplifying TCR ⁇ genes.
- FIG. 3A illustrates a PCR product to be sequenced that was amplified using the scheme of FIGS. 2A-2B .
- FIG. 3B illustrates details of determining a nucleotide sequence of the PCR product of FIG. 3A .
- FIG. 3C illustrates details of another embodiment of determining a nucleotide sequence of the PCR product of FIG. 3A .
- FIG. 4A illustrates a PCR scheme for generating three sequencing templates from an IgH chain in a single reaction.
- FIGS. 4B-4C illustrates a PCR scheme for generating three sequencing templates from an IgH chain in three separate reactions after which the resulting amplicons are combined for a secondary PCR to add P5 and P7 primer binding sites.
- FIG. 4D illustrates the locations of sequence reads generated for an IgH chain.
- FIG. 4E illustrates the use of the codon stricture of V and J regions to improve base calls in the NDN region.
- the practice of the present invention may employ, unless otherwise indicated, conventional techniques and descriptions of molecular biology (including recombinant techniques), bioinformatics, cell biology, and biochemistry, which are within the skill of the art.
- conventional techniques include, but are not limited to, sampling and analysis of blood cells, nucleic acid sequencing and analysis, and the like. Specific illustrations of suitable techniques can be had by reference to the example herein below. However, other equivalent conventional procedures can, of course, also be used.
- Such conventional techniques and descriptions can be found in standard laboratory manuals such as Genome Analysis: A Laboratory Manual Series (Vols.
- the invention is directed to a method for obtaining and analyzing sequence data from a repertoire of immune molecules, such as T cell receptors (TCRs), to rapidly and efficiently determine a clonotype profile.
- Sequence data typically comprises a large collection of sequence reads, i.e. sequences of base calls and associated quality scores.
- a key challenge in constructing clonotype profiles is to rapidly and accurately distinguish sequence reads that contain genuine differences from those that contain errors from non-biological sources, such as the sequencing chemistry, amplification chemistry, or the like.
- An aspect of the invention is the use of sequence trees to represent and organize sequence data from immune molecules extracted from a sample.
- Such immune molecules typically form an immune repertoire which comprises a very laree set of very similar polynucleotides (e.g. >1000, but more usually from 100,000 to 1,000,000, or more) which arc relatively short in length (e.g. usually less than 300 bp).
- a very laree set of very similar polynucleotides e.g. >1000, but more usually from 100,000 to 1,000,000, or more
- arc relatively short in length e.g. usually less than 300 bp.
- sequence tree data structures e.g. lists, arrays, hash table, or the like.
- FIG. 1A illustrates diagrammatically a typical transcript of an IgH molecule ( 120 ) from which a clonotype profile is derived in accordance with some embodiments of the invention.
- Sources of natural sequence variability include modular recombination of the C, D, J and V segments from large sets carried by the genome, nucleotide additions and deletions to the ends of the D segment to produce the so-called “NDN” regions, and somatic hypermutation where substitutions are made randomly over the length of transcript ( 122 ) at a relative frequency roughly as indicated by curve ( 128 ).
- complex populations of such IgH and TCR transcripts are amplified and sequenced.
- one or both operations for IgH molecules are carried out by using redundant primers annealing to different sites in the V regions (described more fully below). This is particularly advantageous where a sequencing chemistry is employed that has a relatively high error rate or where such sequence variability is difficult or impossible to know beforehand. In the latter case, primer extension for amplification or generation of sequence reads takes place even if one or more primer binding sites are inoperable, or substantially inoperable, because of mismatches caused (for example) by one or more somatic mutations.
- a segment of recombined B cell nucleic acid is amplified by a PCR with a plurality of forward primers or a plurality of reverse primers to generate a nested set of templates (see FIGS. 4A and 4B and their descriptions below). Templates from such a set may be further amplified on a surface to form separate amplicons (e.g. by bridge PCR using a cBot instrument; Illumina, San Diego, Calif.).
- Templates from the same nested set may be associated with one another by sequence reads generated at their common ends. Nested sets of templates allow a sequencing chemistry with relative high error rates to be used to analyze longer sequences than otherwise would be possible, while at the same time maintaining high average quality scores over the entire length of the sequence. The nested sets also ensure that at least one sequence read is obtained from a V region even if it has been subjected to somatic hypermutation.
- sequencing chemistries may be used for analyzing highly variable nucleic acids, such as IgH molecules, that have error rates no better than the following: 0.2 percent of sequence reads contain at least one error in positions 1-50; 0.2-1.0 percent of sequence reads contain at least one error in positions 51-75; 0.5-1.5 percent of sequence reads contain at least one error in positions 76-100; and 1-5 percent of sequence reads contain at least one error in positions 101-125.
- the method of the invention includes steps for distinguishing clonotype sequences that are closely related and genuinely different from those that are closely related and the result of sequencing or other error.
- Constructing clonotypes from sequence read data depends in part on the sequencing method used to generate such data, as the different methods have different expected read lengths and data quality.
- a Solexa sequencer is employed to generate sequence read data for analysis.
- a sample is obtained that provides at least 0.5-1.0 ⁇ 10 6 lymphocytes to produce at least 1 million template molecules, which after optional amplification may produce a corresponding one million or more clonal populations of template molecules (or clusters).
- the Solexa approach such over sampling at the cluster level is, desirable so that each template sequence is determined with a large degree of redundancy to increase the accuracy of sequence determination.
- the sequence of each independent template is determined 10 times or more.
- different levels of redundancy may be used for comparable accuracy of sequence determination.
- coalescing of sequence reads depends on at least three factors: the number of sequences obtained for each of the clonotypes being compared; the number of bases at which they differ: and the sequencing quality score at the positions at which they are discordant.
- a likelihood ratio may be constructed and assessed that is based on the expected error rates and binomial distribution of errors. For example, two clonotypes, one with 150 reads and the other with 2 reads with one difference between them in an area of poor sequencing quality will likely be coalesced as they are likely to be generated by sequencing error.
- sequence reads are first converted into candidate clonotypes. Such a conversion depends on the sequencing platform employed. For platforms that generate high Q score long sequence reads, the sequence read or a portion thereof may. be taken directly as a candidate clonotype.
- candidate clonotypes are generated from collections of paired reads from multiple clusters, e.g. 10 or more, as mentioned above.
- the frequencies of candidate clonotypes may be plotted in sequence space, as illustrated in FIG. 1B , where such space is reduced to one dimension (the horizontal axis) for sake of illustration.
- the vertical axis gives the magnitude of each candidate clonotype's frequency, log(read count), or some like measure.
- candidate clonotypes are represented by the various symbols ( 130 ).
- whether two candidate clonotypes are coalesced depends on their respective frequencies or read counts (as noted above), the number of base differences between them (the more differences, the less likely is coalescence), and the quality scores of the bases at the locations where the respective sequences differ (higher quality scores makes coalescence less likely).
- Candidate clonotypes may be considered in the order of their respective frequencies.
- FIG. 1B shows candidate clonotype 1 ( 132 ), candidate clonotype 7 ( 134 ) and candidate clonotype 11 ( 136 ) as the three candidates with the highest three frequencies.
- each such candidate clonotype is other candidate clonotypes that are close in sequence, but with lesser frequencies, such as (i) for candidate clonotype 1 ( 132 ) there are candidate clonotype 2 ( 138 ) and the candidate clonotypes 3, 4, 5 and 6 enclosed by cone ( 140 ); for candidate clonotype 7 ( 134 ) there arc candidate clonotypes 8, 9 and 10 enclosed by cone ( 142 ); and (iii) for candidate clonotype 11, there is candidate clonotype 12 enclosed by cone ( 144 ).
- the cones represent likelihood boundaries within which a lesser frequency candidate clonotype would be coalesced with one of the higher frequency candidate clonotypes 1, 7 or 11.
- Such likelihood boundaries are functions of the frequency of the nearby candidate clonotypes (3, 4, 5 and 6 for 1; 8, 9 and 10 for 7; and 12 for 11) and their distances in sequence space from the respective higher frequency candidate clonotypes.
- Candidate clonotype 2 ( 138 ) is outside cone ( 140 ); thus, it would not be coalesced with candidate clonotype 1 ( 132 ).
- the likelihood (of coalesce) boundaries are shown as cones because candidate clones with higher frequencies are more likely to be genuinely different clonotypes than those of lower frequencies and multiple differences at lower frequencies are more likely to be errors than multiple differences at higher frequencies.
- the cloud of sequence reads surrounding each candidate cloncotype can be modeled using the binomial distribution and a simple, model for the probability of a single base error.
- This latter error model can be inferred from mapping V and J segments or from the clonotype finding algorithm itself, via self-consistency and convergence.
- a model is constructed for the probability of a given ‘cloud’ sequence Y with read count C2 and E errors (with respect to sequence X) being part of a true clonotype sequence X with perfect read count C1 under the null model that X is the only true clonotype in this region of sequence space.
- a decision is made whether or not to coalesce sequence Y into the clonotype X according the parameters C1, C2, and E.
- a max value C2 is pre-calculated for deciding to coalesce the sequence Y.
- the max values for C2 are chosen so that the probability of failing to coalesce Y under the null hypothesis that Y is part of clonotype X is less than some value P after integrating over all possible sequences Y with error E in the neighborhood of sequence X.
- the value P is controls the behavior of the algorithm and makes the coalescing more or less permissive.
- a sequence Y is not coalesced into clonotype X because its read count is above the threshold C2 for coalescing into clonotype X then it, becomes a candidate for seeding separate clonotypes (such as with candidate clonotype 2 ( 138 ) in FIG. 1B ).
- An algorithm implementing such principles would also make sure that any other sequences Y2, Y3, etc. which are ‘nearer’ to this sequence Y (that had been deemed independent of X) are not aggregated into X. This concept of ‘nearness’ includes both error counts with respect to Y and X and the absolute read count of X and Y, i.e.
- an algorithm proceeds in a top down fashion by starting with the sequence X with the highest read count. This sequence seeds the first clonotype. Neighboring sequences are either coalesced into this clonotype if their counts are below the precalculated thresholds (see above), or left alone if they are above the threshold or ‘closer’ to another sequence that was not coalesced. After searching all neighboring sequences within a maximum error count, the process of coalescing reads into clonotype X is finished. Its reads and all reads that have been coalesced into it are accounted for and removed from the list of reads available for making other clonotypes. The next sequence is then moved on to with the highest read count.
- Neighboring reads are coalesced into this clonotype as above and this process is continued until there are no more sequences with read counts above a given threshold, e.g. until all sequences with more than 1 count have been used as seeds for clonotypes.
- a further test may be added for determining whether to coalesce a candidate sequence Y into an existing clonotype X, which takes into account quality score of the relevant sequence reads.
- the average quality store(s) are determined for sequence(s) Y (averaged across all reads with sequence Y) were sequences Y and X differ. If the average score is above a predetermined value then it is more likely that the difference indicates a truly different clonotype that should not be coalesced and if the average score is below such predetermined value then it is more likely that sequence Y is caused by sequencing errors and therefore should be coalesced into X.
- All of the reads of a given sample are placed into the sequence tree. Each leaf nodes holds pointers to its associated reads.
- a unique sequence of a candidate clonotype is retrieved by traversing backwards in the tree from the leaf to the root node. The first sequence is placed into a simple tree with one root node and one leaf node that contains the full sequence of the read. Sequences arc next added one by one. For each added sequence either a new branch is formed at the last point of common sequence between the read and the existing tree or add the read to an existing leaf node if the tree already contains the sequence.
- Finding the highest read count sorting leaf nodes by read count allows one to find the leaf node (i.e. sequence) with the most reads, and successively lower numbers of reads; 2) Finding neighboring leafs: for any sequence all paths through the tree which have less than X errors with respect to this sequence are searchable.
- a path is started at the root and branch this path into separate paths proceeding along the tree. The current error count of each path as proceeding along the tree is noted. When the error count exceeds the max allowed errors the given path is terminated. In this way large parts of the tree arc pruned as early as possible. This is an efficient way of finding all paths.(i.e. all leafs) within X errors from any given sequence.
- a set of candidate clonotypes is obtained from sequence data obtained by sequencing recombined nucleic acids extracted from a sample of T cells or B cells.
- candidate clonotypes each include an NDN region and portions of V and J regions. These sequences are organized into a data structure ( 150 ), which may be a sequence tree.
- sequence trees may also be constructed for known V regions and known J regions.
- Sequence reacts making up a candidate clonotype may then be mapped, or aligned, to these known sequences via the sequence trees to, efficiently determine the most likely known V and J sequences of the candidate clonotypes.
- a data structure such as a sequence tree, is constructed for uses in a method for distinguishing genuine clonotypes from candidate clonotypes that contain experimental or measurement errors, such as sequencing errors.
- the candidate clonotype that has the highest frequency of occurrence among the current candidate clonotypes is selected ( 152 ) from the data structure, for example a sequence tree; in other words, HFCC k is the candidate clonotype with the highest number of copies, or read counts in cycle k.
- HFCC k is the candidate clonotype with the highest number of copies, or read counts in cycle k.
- LFCCs neighboring lesser frequency candidate clonotypes arc identified
- this identification is carried out using a sequence tree, which allows efficient sequence comparisons of relatively short ( ⁇ 300 bp) sequences.
- the comparisons, or sequence alignments are carried out using dynamic programming, e.g. as disclosed by Gusfield (cited above).
- dynamic programming is banded dynamic programming where sequences that differ from the selected HFCC by more than a predetermined distance are not considered, which speeds the computation.
- the candidates HFCC k and LFCC j may be compared on the basis of many different criteria or properties.
- candidate clonotypes arc compared on the basis of at least two properties: (i) frequency or read counts and (ii) sequence differences.
- candidate clonotypes are compared on the basis of at least three properties: (i) frequency or read counts, (ii) sequence differences, and (iii) quality scores or measures of the bases where differences occur.
- sequence differences include base substitutions; in another embodiment, sequence differences include base substitutions, deletions and insertions. The latter embodiment is especially applicable whenever sequence data is generated by sequencing-by-synthesis chemistries that do not employ terminators, such as 454 sequencers and Ion Torrent sequencers.
- comparisons of HFCCs and LFCCs may be implemented using a function, such as P(HFCC k , LFCC j , D, Q) shown in decision box ( 158 ), which depends on the quantities (i) through (iii) described above.
- Such a function may take many different forms, but generally the value of P changes with changes in (i), (ii) and (iii) as follows:
- the value of P preferably increases monotonically with the frequency of HFCC and the ratio of the frequency of HFCC to that of LFCC, such that the higher the ratio of the frequency of HFCC to that of LFCC, the higher the likelihood LFCC will be coalesced into HFCC.
- the value of P preferably decreases monotonically with degree to which the sequences of HFCC and LFCC differ, so that the greater the difference between HFCC and LFCC (e.g.
- the value of P preferably decreases monotonically with increasing quality scores of the locations where the sequences of HFCC and LFCC differ, so that for higher quality scores, the lower the likelihood LFCC will be coalesced with HFCC.
- the quality scores at the different locations may be combined in a variety of differ ways. In one embodiment, whenever there is a plurality of such differences, the plurality of quality scores is expressed as an average value, which may be either an unweighted average or a weighted average.
- FIG. 1D shows an exemplary function, P, computed for different quality values (curves a through e) for a given sequence difference.
- P computed for different quality values
- D D k
- D D(HFCC frequency)
- D D(HFCC frequency)
- D distance between HFCC k and LFCC j ; however, other distance measures may be used.
- D k is a non-increasing function of k: and in another embodiment, D k is a decreasing function of k.
- FIG. 1D illustrates computed values for one embodiment of P which relates read count thresholds for coalescing an LFCC given a read count of an HFCC for different quality scores, as described above. Curves “a” through “c” represent the relationships for different quality scores (with curve “a” corresponding to the highest quality score).
- path ( 192 ) (indicated by dashed line) in sequence tree ( 190 ) corresponds to HFCC ( 196 ), which is coalesced with LFCC ( 198 ).
- segment of path ( 192 ) in shaded area ( 199 ) is removed from sequence tree ( 190 )to give reduced sequence tree ( 197 ) shown in FIG. 1F , which is used in subsequent computations to find neighboring LFCC ( 154 ).
- clonotype determination is finished if a stopping criterion ( 170 ) is met.
- stopping criterion ( 170 ) is whether the last non-singleton candidate clonotype has been processed ( 152 ).
- stopping criterion ( 170 ) is whether the frequency or the read counts of the selected HFCC is below that corresponding to a single lymphocyte.
- an amplification step may result in each lymphocyte in a sample being represented by multiple copies of the same clonotype; thus, in one embodiment, whatever HFCC has a number of read counts below the number corresponding to a single lymphocyte, then the computation is stopped.
- such a number of read counts (or candidate clonotype copies) is 10; in another embodiment, such number is 20; in another embodiment, such a number is 30; in another embodiment, such a number is 40.
- the next HFCC is selected ( 172 ).
- the analytical steps summarized in the flow chart of FIG. 1C may be implemented in any suitable programming language, such as C, C++. Java, C#, Fortran, Pascal or the like.
- the above method for determining clonotypes and/or clonotype profiles comprises steps of (a) forming a data structure of recombined immune molecules from sequence reads obtained by high throughput nucleic acid sequencing, (b) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype, (c) removing the coalesced candidate clonotype from the data structure, and (d) repeatine, steps (b) and (c) until a clonotype profile is formed.
- the data structure is a sequence tree.
- the above method of determining clonotypes may be carried out by steps comprising: (a) providing a set of sequence reads from a repertoire of recombined immune molecules each having a V region, an NDN region and a J region wherein for each such molecule at least one sequence read encompasses at least a portion of the NDN region of such molecule; (b) forming from sequence reads encompassing at least a portion of an NDN region a sequence tree having leaves representing candidate clonotypes, each leaf and its corresponding candidate clonotype having a frequency; (c) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype having a sequence of the highest frequency candidate clonotype; (d) removing leaves corresponding to the coalesced candidate clonotypes from the sequence
- the step of forming further includes selecting a highest frequency candidate clonotype and identifying all said lesser frequency candidate clonotypes having a sequence difference therewith less than a predetermined difference value to form a coalescence subset.
- a predetermined difference value For the step of forming, one may limit the total number of LFCCs that must be compared for the coalescing operation (only ones within the predetermined difference value are considered).
- Such value is a process input depending on the application, e.g. the size of the repertoire, how much computing time is used, and so on.
- the function used for deciding whether to coalesce an HFCC with a LFCC can have a variety of forms. In one general aspect, for the step of coalescing, such a function may have the following properties.
- selection of a predetermined frequency value and a predetermined difference value is a design choice that depend on particular applications. Factors affecting such choices may include details of the biology, speed of implementation, and the like.
- Immune system repertoires may be obtained from samples of immune cells.
- immune cells can include T-cells and/or B-cells.
- T-cells include, for example, cells that express T cell receptors.
- T-cells include helper T cells (effector T cells or Th cells), cytotoxic T cells (CTLs), memory T cells, and regulatory T cells.
- a sample of T cells includes at least 1,000T cells; but more typically, a sample includes at least 10,000 T cells, and more typically, at least 100;000 T cells.
- a sample includes a number of T cells. In the range of from 1000 to 1,000,000 cells.
- a sample of immune cells may also comprise B cells.
- B-cells include, for example, plasma B cells, memory B cells, B1 cells, B2 cells, marginal-zone B cells, and follicular B cells.
- B-cells can express immunoglobulins (antibodies, B cell receptor).
- a sample of B cells includes at least 1,000 B cells; but more typically, a sample includes at least 10,009 B cells, and more typically, at least 100,000 B cells.
- a sample includes a number of Bcells in the range of from 1000 to 1,000,000 B cells.
- the sample can include nucleic acid, for example, DNA (e.g., genomic DNA or thitocliondrial DNA) or RNA (e.g., messenger RNA or microRNA).
- the nucleic acid can be cell-free DNA or RNA, e.g. extracted from the circulatory system, Vlassov et al, Curr. Mol. Med., 10: 142-165 (2010); Swamp et al, FEBS Lett., 581: 795-799 (2007).
- the amount of RNA or DNA from a subject that can be analyzed includes, for example, as low as a single cell in some applications (e.g., a calibration test) and as many as 10 million of cells or more translating to a range of DNA of 6 pg-60 ug, and RNA of approximately 1 pg-10 ug.
- a sample of lymphocytes is sufficiently large so that substantially every T cell or B cell with a distinct clonotype is represented therein, thereby forming a repertoire (as the term is used herein).
- a sample is taken that contains with a probability of ninety-nine percent every clonotype of a population present at a frequency of 0.001 percent or greater.
- a sample is taken that contains with a probability of ninety-nine percent every clonotype of a population present at a frequency of 0.0001 percent or greater.
- a sample of B cells or T cells includes at least a half million cells, and in another embodiment such sample includes at least one million cells.
- DNA from the material may be amplified by a non-biasing technique, such as whole genome amplification (WGA), multiple displacement amplification (MDA); or like technique, e.g. Hawkins et al, Curr. Opin. Biotech., 13: 65-67 (2002); Dean et al, Genome Research, 11: 1095-1099 (2001); Wang et al, Nucleic. Acids Research, 32: e76 (2004); Hosono ct al, Genome Research, 13: 954-964 (2003); and the like.
- WGA whole genome amplification
- MDA multiple displacement amplification
- Blood samples are of particular interest, especially in monitoring lymphoid neoplasms, such as lymphomas, leukemias, or the like; and may be obtained using conventional techniques, e.g. Innis et al, editors, PCR Protocols (Academic Press, 1990); or the like.
- white blood cells may be separated from blood samples using convention techniques, e.g. RosetteSep kit (Stem Cell Technologies, Vancouver, Canada).
- Blood samples may range in volume from 100 ⁇ L to 10 mL; in one aspect, blood sample volumes are in the range of from 200 100 ⁇ L to 2 mL.
- DNA and/or RNA may then be extracted from such blood sample using conventional techniques for use in methods of the invention, e.g.
- DNeasy Blood & Tissue Kit Qiagen, Valencia, Calif.
- subsets of white blood cells e.g. lymphocytes
- FACS fluorescently activated cell sorting
- MCS magnetically activated cell sorting
- RNA or DNA can be sequenced in the methods of the provided invention.
- a recombined sequence from a T-cell or B-cell encoding a T cell receptor or immunoglobulin molecule, or a portion thereof, is referred to as a clonotype.
- the DNA or RNA can correspond to sequences from T-cell receptor (TCR) genes or immunoglobulin (Ig) genes that encode antibodies.
- TCR T-cell receptor
- Ig immunoglobulin
- the DNA and RNA can correspond to sequences encoding ⁇ , ⁇ , ⁇ , or ⁇ chains of a TCR.
- the TCR is a heterodimer consisting of an ⁇ -chain and ⁇ -chain.
- the TCR ⁇ chain is generated by VJ recombination, and the ⁇ chain receptor is generated by V(D)J recombination.
- VJ recombination For the TCR ⁇ chain, in humans there arc 48 V segments, 2 D segments, and 13 J segments. Several bases may be deleted and others added (called N and P nucleotides) at each of the two junctions.
- N and P nucleotides Several bases may be deleted and others added (called N and P nucleotides) at each of the two junctions.
- the TCRs consist of ⁇ and ⁇ delta chains.
- the TCR ⁇ chain is generated by VJ recombination, and the TCR ⁇ chain is generated by V(D)J recombination (Kenneth Murphy, Paul Travers, and Mark Walport, Janeway's Immunology 7th edition, Garland Science, 2007, which is herein incorporated by reference in its entirety).
- the DNA and RNA analyzed in the methods of the invention can correspond to sequences encoding heavy chain immunoglobulins ago with constant regions ( ⁇ , ⁇ , ⁇ , ⁇ , or ⁇ ) or light chain immunoglobulins (IgK or IgL) with constant regions ⁇ or ⁇ .
- Each antibody has two identical light chains and two identical heavy chains.
- Each chain is composed of a constant (C) and a variable region.
- the variable region is composed of a variable (V), diversity (D), and joining (J) segments.
- V variable
- D diversity
- J joining
- Diversity in the light chain is generated in a similar fashion except that there is no D region so there is only VJ recombination. Somatic mutation often occurs close to the site of the recombination, causing the addition or deletion of several nucleotides, further increasing the diversity of heavy and light chains generated by B-cells. The possible diversity of the antibodies generated by a B-cell is then the product of the different heavy and light chains. The variable regions of the heavy and light chains contribute to form the antigen recognition (or binding) region or site. Added to this diversity is a process of somatic hypermutation which can occur after a specific response is mounted against some epitope.
- primers may be selected to generate amplicons of subsets of recombined nucleic acids extracted from lymphocytes. Such subsets may be referred to herein as “somatically rearranged regions.” Somatically rearranged regions may comprise nucleic acids from developing or from fully developed lymphocytes, where developing lymphocytes are cells in which rearrangement of immune genes has not been completed to form molecules having full V(D)J regions. Exemplary incomplete somatically rearranged regions include incomplete IgH molecules (such as, molecules containing only D-J regions), incomplete TCR ⁇ molecules (such as, molecules containing only D-J regions), and inactive IgK (for example, comprising Kdc-V regions).
- Adequate sampling of the cells is an important aspect of interpreting the repertoire data, as described further below in the definitions of “clonotype” and “repertoire”
- clonotype and “repertoire”
- one aspect of this invention is the development of methods to quantitate the number of input immune receptor molecules. This has been implemented this for TCR ⁇ and IgH sequences. In either case the same set of primers are used that are capable of amplifying all the different sequences.
- a real time PCR with the multiplex of primers is performed along with a standard with a known number of immune receptor copies.
- This real time PCR measurement can be made from the amplification reaction that will subsequently be sequenced or can be done on a separate aliquot of the same sample.
- DNA the absolute number of rearranged immune receptor molecules can be readily converted to number of cells (within 2 fold as some cells will have 2 rearranged copies of the specific immune receptor assessed and others will have one).
- cDNA the measured total number of rearranged molecules in the real time sample can be extrapolated to define the total number of these molecules used in another amplification reaction of the same sample.
- this method can be combined with a method to determine the total amount of RNA to define the number of rearranged immune receptor molecules in a unit amount (say 1 ⁇ g) of RNA assuming a specific efficiency of cDNA synthesis. If the total amount of cDNA is measured then the efficiency of cDNA synthesis need not be considered. If the number of cells is also known then the rearranged immune receptor copies per cell can be computed. If the number of cells is not known, one can estimate it from the total RNA as cells of specific type usually generate comparable amount of RNA. Therefore from the copies of rearranged immune receptor molecules per 1 ⁇ g one can estimate the number of these molecules per cell.
- Another approach that can be utilized is to add it known amount of unique immune receptor rearranged molecules with a known sequence, i.e. known amounts of one or more internal standards, to the cDNA or gnomic DNA from a sample of unknown quantity.
- a known sequence i.e. known amounts of one or more internal standards
- uch techniques for molecular counting arc well-known, e.g. Brenner et al, U.S. Pat. No. 7,537,897, which is incorporated herein by reference).
- Data from sequencing the added unique sequence can be used to distinguish the different possibilities if a real time PCR calibration is being used as well.
- amplicons of target populations of nucleic acids may be generated by a variety of amplification techniques.
- multiplex PCR is used to amplify members of a mixture of nucleic acids, particularly mixtures comprising recombined immune molecules such as T cell receptors, B cell receptors, or portions thereof.
- Guidance for carrying out multiplex PCRs of such immune molecules is found in the following references, which are incorporated by reference: Morley, U.S. Pat. No. 5,296,351; Gorski, U.S. Pat. No. 5,837,447; Dau, U.S. Pat. No. 6,087,096; Von Dongen et al, U.S.
- the step of spatially isolating individual nucleic acid molecules is achieved by carrying out a primary multiplex amplification of a preselected somatically rearranged region or portion thereof (i.e. target sequences) using forward and reverse primers that each have tails non-complementary to the target sequences to produce a first amplicon whose member sequences have common sequences at each end that allow further manipulation.
- such common ends may include primer binding sites for continued amplification using just a single forward primer and a single reverse primer instead of multiples of each, or for bridge amplification of individual molecules on a solid surface, or the like.
- Such common ends may be added in a single amplification as described above, or they may be added in a two-step procedure to avoid difficulties associated with manufacturing and exercising quality control over mixtures of long primers (e.g. 50-70 bases or more).
- the primary amplification is carried out as described above, except that the primer tails are limited in length to provide only forward and reverse primer binding sites at the ends of the sequences of the first amplicon.
- a secondary amplification is then carried out using secondary amplification primers specific to these primer binding sites to add further sequences to the ends of a second amplicon.
- the secondary amplification primers have tails non-complementary to the target sequences, which form the ends of the second amplicon and which may, be used in connection With sequencing the clonotypes of the second amplicon.
- such added sequences may include primer binding sites for generating sequence reads and primer binding sites for carrying out bridge PCR on a solid surface to generate clonal populations of spatially isolated individual molecules, for example, when Solexa-based sequencing is used.
- a sample of sequences from the second amplicon are disposed on a solid surface that has attached complementary oligonucleotides capable of annealing to sequences of the sample, idler which cycles of primer extension, denaturation, annealing are implemented until clonal populations of templates are formed.
- the size of the sample is selected so that (i) it includes an effective representation of clonotypes in the original sample, and (ii) the density of clonal populations on the solid surface is in a range that permits unambiguous sequence determination of clonotypes.
- TCR or BCR sequences or portions thereof can be amplified from nucleic acid in a multiplex reaction using at least one primer that anneals to the C region and one or more primers that can anneal to one or more V segments (as illustrated in FIGS. 2A-2B and FIGS. 4A-4B and discussed more fully below).
- the region to be amplified can include the full clonal sequence or a subset of the clonal sequence, including the V-D junction, D-J junction of an inummoglobulin or T-cell receptor gene, the full variable region of an immunoglobulin or T-cell receptor gene, the antigen recognition region, or a CDR, e.g., complementarity determining region 3 (CDR3).
- CDR3 complementarity determining region 3
- the TCR or immunoglobulin sequence can amplified using a primary and a secondary amplification step.
- Each of the different amplification steps can comprise different primers.
- the different primers can introduce sequence not originally present in the immune gene sequence.
- the amplification procedure can add new primer binding sites to the ends of the target sequences to convert a multiplex amplification to a singleplex amplification or the amplification procedure can add one or more tags to the 5′ and/or 3′ end of amplified TCR or immunoglobulin sequence (as illustrated in FIGS. 3A-3B ).
- the tag can be sequence that facilitates subsequent sequencing of the amplified DNA.
- the tag can be sequence that facilitates binding the amplified sequence to a solid support.
- the individual nucleic acid molecules can be isolated, optionally re-amplified, and then sequenced individually.
- Exemplary amplification protocols may be found in van Dongen et al, Leukemia, 17: 2257-2317 (2003) or van Dongen et al, U.S. patent publication 2006/0234234, which is incorporated by reference.
- an exemplary protocol is as follows: Reaction buffer: ABI Buffer II or ABI Gold Buffer (Life Technologies, San Diego, Calif.); 50 ⁇ L final reaction volume; 100 ng sample DNA; 10 pmol of each primer (subject to adjustments to balance amplification as described below); dNTPs at 200 ⁇ M final concentration; MgCl 2 at 1.5 mM final concentration (subject to optimization depending on target sequences and polymerase); Taq polymerase ( 1 - 2 U/tube); cycling conditions: preactivation 7 min at 95° C.; annealing at 60° C.; cycling times: 30 s denaturation; 30 s annealing; 30 s extension.
- Reaction buffer ABI Buffer II or ABI Gold Buffer (Life Technologies, San Diego, Calif.); 50 ⁇ L final reaction volume; 100 ng sample DNA; 10 pmol of each primer (subject to adjustments to balance amplification as described below); dNTPs at 200 ⁇ M final concentration; MgCl 2 at 1.5 mM final
- Polymerases that can be used for amplification in the methods of the invention are commercially available and include, for example, Taq polymerase, AccuPrime polymerase, or Pfu.
- the choice of polymerase to use can be based on whether fidelity or efficiency is preferred.
- Methods for isolation of nucleic acids from a pool include subcloning nucleic acid into DNA vectors and transforming bacteria (bacterial cloning), spatial separation of the molecules in two dimensions on a solid substrate (e.g., glass slide), spatial separation of the molecules in three dimensions in a solution within micelles (such as can be achieved using oil emulsions with or without immobilizing the molecules on a solid surface such as beads), or using microreaction chambers in, for example, microfluidic or nano-fluidic chips. Dilution can be used to ensure that on average a single molecule is present in a given volume, spatial region, bead, or reaction chamber.
- Real time PCR, picogreen staining, nanofluidic electrophoresis (e.g. LabChip) or UV absorption measurements can be used in an initial step to judge the functional amount of amplifiable material.
- multiplex amplifications are carried out so that relative amounts of sequences in a starting population arc substantially the same as those in the amplified population, or amplicon. That is, multiplex amplifications are carried out with minimal amplification bias among member sequences of a sample population. In one embodiment, such relative amounts are substantially the same if each relative amount in an amplicon is within five fold of its value in the starting sample. In another embodiment, such relative amounts are substantially the same if each relative amount in an amplicon is within two fold of its value in the starting sample. As discussed more fully below, amplification bias in PCR may be detected and corrected using conventional techniques so that a set of PCR primers may be selected for a predetermined repertoire that provide unbiased amplification of any sample.
- a multiplex amplification optionally uses all the V segments.
- the reaction is optimized to attempt to get amplification that maintains the relative abundance of the sequences amplified by different V segment primers.
- Some of the primers are related, and hence many of the primers may “cross talk,” amplifying templates that arc not perfectly matched with it.
- the conditions arc optimized so that each template can be amplified in a similar fashion irrespective of which primers amplified it. In other words if there are two templates, then after 1,000 fold amplification both templates can be amplified approximately 1,000 fold, and it does not matter that for one of the templates half of the amplified products carried a different primer because of the cross talk.
- the primer sequence is eliminated from the analysis, and hence it does not matter what primer is used in the amplification as long as the templates are amplified equally.
- amplification bias may be avoided by carrying out a two-stage amplification (as illustrated in FIGS. 2A-2B ) wherein a small number of amplification cycles are implemented in a first, or primary, stage using primers having tails non-complementary with the target sequences.
- the tails include primer binding sites that arc added to the ends of the sequences of the primary amplicon so that such sites arc used in a second stage amplification using only a single forward primer and a single reverse primer, thereby eliminating a primary cause of amplification bias.
- the primary PCR will have a small enough number of cycles (e.g. 5-10) to minimize the differential amplification by the different primers.
- the secondary amplification is done with one pair of primers and hence the issue of differential amplification is minimal.
- One percent of the primary PCR is taken directly to the secondary PCR. Thirty-five cycles (equivalent to 28 cycles without the 100 fold dilution step) used between the two amplifications were sufficient to show a robust amplification irrespective of whether the breakdown of cycles were: one cycle primary and 34 secondary or 25 primary and 10 secondary. Even though ideally doing only 1 cycle in the primary PCR may decrease the amplification bias, there arc other considerations.
- One aspect of this is representation. This plays a role when the stalling input amount is not in excess to the number of reads ultimately obtained.
- the 100 fold dilution between the 2 steps means that the representation is reduced unless the primary PCR amplification generated significantly more than 100 molecules. This indicates that a minimum 8 cycles (256 fold), but more comfortably 10 cycle (1,000 fold), may be used.
- the alternative to that is to take more than 1% of the primary PCR into the secondary but because of the high concentration of primer used in the primary PCR, a big dilution factor is can be used to ensure these primers do not interfere in the amplification and worsen the amplification bias between sequences.
- Another alternative is to add a purification or enzymatic step to eliminate the primers from the primary PCR to allow a smaller dilution of it. In this example, the primary PCR was 10 cycles and the second 25 cycles.
- any high-throughput technique for sequencing nucleic acids can be used in the method of the invention.
- such technique has a capability of generating in a cost-effective manner a volume of sequence data from which at least 1000 clonotypes can be determined, and preferably, from which at least 10,000 to 1,000,000 clonotypes can be determined.
- DNA sequencing techniques include classic dideoxy sequencing reactions (Sanger method) using labeled terminators or primers and gel separation in slab or capillary, sequencing by synthesis using reversibly terminated labeled nucleotides, pyrosequencing, 454 sequencing, allele specific hybridization to a library of labeled oligonucleotide probes, sequencing by synthesis using allele specific hybridization to a library of labeled clones that is followed by ligation, real time monitoring of the incorporation of labeled nucleotides during a polymerization step, polony sequencing, and SOLiD sequencing.
- Sequencing of the separated molecules has more recently been demonstrated by sequential or single extension reactions using polymerases or ligases as well as by single or sequential differential hybridizations with libraries of probes. These reactions have been performed on many clonal sequences in parallel including demonstrations in current commercial applications of over 100 million sequences in parallel. These sequencing approaches can thus be used to study the repertoire of T-cell receptor (TCR) and/or B-cell receptor (BCR).
- high-throughput methods of sequencing are employed that comprise a step of spatially isolating individual molecules on a solid surface where they are sequenced in parallel.
- Such solid surfaces may include nonporous surfaces (such as in Solexa sequencing, e.g.
- such methods comprise amplifying the isolated molecules either before or after they arc spatially isolated on a solid surface.
- Prior amplification may comprise emulsion-based amplification, such as emulsion PCR, or rolling circle amplification.
- emulsion-based amplification such as emulsion PCR, or rolling circle amplification.
- Solexa-based sequencing where individual template molecules are spatially isolated on a solid surface, after which they are amplified in parallel by bridge PCR to form separate clonal populations, or clusters, and then sequenced, as described in Bentley et al (cited above) and in manufacturer's instructions (e.g.
- individual molecules disposed and amplified on a solid surface form clusters in a density of at least 10 5 clusters per cm 2 ; or in a density of at least 5 ⁇ 10 5 per cm 2 ; or in a density of at least 10 6 clusters per cm 2 .
- sequencing chemistries are employed having relatively high error rates.
- the average quality scores produced by such chemistries are monotonically declining functions of sequence read lengths. In one embodiment, such decline corresponds to 0.5 percent of sequence reads have at least one error in positions 1-75; 1 percent of sequence reads have at least one error in positions 76-100; and 2 percent of sequence reads have at least one error in positions 101-125.
- a sequence-based clonotype profile of an individual is obtained using the following steps: (a) obtaining a nucleic acid sample from T-cells and/or B-cells of the individual; (b) spatially isolating individual molecules derived from such nucleic acid sample, the individual molecules comprising at least one template generated from a nucleic acid in the sample, which template comprises a somatically rearranged region or a portion thereof, each individual molecule being capable of producing at least one sequence read; (c) sequencing said spatially isolated individual molecules; and (d) determining abundances of different sequences of the nucleic acid molecules from the nucleic. acid sample to generate the clonotype profile.
- each of the somatically rearranged regions comprise a V region and a J region.
- the step of sequencing comprises bidirectionally sequencing each of the spatially isolated individual molecules to produce at least one forward sequence read and at least one reverse sequence read. Further to the latter embodiment, at least one of the forward sequence reads and at least one of the reverse sequence reads have an overlap region such that bases of such overlap region are determined by a reverse complementary relationship between such sequence reads.
- each of the somatically rearranged regions comprise a V region and a J region and the step of sequencing further includes determining a sequence of each of the individual nucleic acid molecules from one or more of its forward sequence reads and at least one reverse sequence read starting from a position in a J region and extending in the direction of its associated V region.
- individual molecules comprise nucleic acids selected from the group consisting of complete IgH molecules, incomplete IgH molecules, complete IgK complete, IgK inactive molecules, TCR ⁇ molecules, TCR ⁇ molecules, complete TCR ⁇ molecules, and incomplete TCR ⁇ molecules.
- the step of sequencing comprises generating the sequence reads having monotonically decreasing quality scores.
- monotonically decreasing quality scores are such that the sequence reads have error rates no better than the following: 0.2 percent of sequence reads contain at least one error in base positions 1 to 50, 0.2 to 1.0 percent of sequence reads contain at least one error in positions 51-75, 0.5 to 1.5 percent of sequence reads contain at least one error in positions 76-100.
- the above method comprises the following steps: (a) obtaining a nucleic acid sample from T-cells and/or B-cells of the individual; (b) spatially isolating individual molecules derived from such nucleic acid sample, the individual molecules comprising nested sets of templates each generated from a nucleic acid in the sample and each containing a somatically rearranged region or a portion thereof, each nested set being capable of producing a plurality of sequence reads each extending in the same direction and each starting from a different position on the nucleic acid from which the nested set was generated; (c) sequencing said spatially isolated individual molecules; and (d) determining abundances of different sequences of the nucleic acid molecules from the nucleic acid sample to generate the clonotype profile.
- the step of sequencing includes producing a plurality of sequence reads for each of the nested sets.
- each of the somatically rearranged regions comprise a V region and a J region, and each of the plurality of sequence reads starts from a different position in the V region and extends in the direction of its associated J region.
- the sequencing technique used in the methods of the invention generates sequences of least 1000 clonotypes per run; in another aspect, such technique generates sequences of at least 10,000 clonotypes per run; in another aspect, such technique generates sequences of at least 100,000 clonotypes per run; in another aspect, such technique generates sequences of at least 500,000 clonotypes per run; and in another aspect, such technique generates sequences of at least 1,000,000 clonotypes per run. In still another aspect, such technique generates sequences of between 100,000 to 1,000,000 clonotypes per run per individual sample.
- the sequencing technique used in the Methods of the provided invention can generate about 30 bp, about 40 bp, about 50 bp, about 60 bp, about 70 bp, about 80 bp, about 90 bp, about 100 bp, about 110, about 120 by per read, about 150 bp, about 200 bp, about 250 bp, about 300 bp, about 350 bp, about 400 bp, about 450 bp, about 500 bp, about 550 bp, or about 600 by per read.
- sequences of clonotypes may be determined by combining information from one or more sequence reads, for example, along the V(D)J regions of the selected chains.
- sequences of clonotypes are determined by combining information from a plurality of sequence reads.
- Such pluralities of sequence reads may include one or more sequence reads along a sense strand (i.e. “forward” sequence reads) and, one or more sequence reads along, its complementary strand (i.e. “reverse” sequence reads).
- FIGS. 4B-4C The generation of templates in separate reactions is illustrated in FIGS. 4B-4C .
- a sample containing IgH ( 400 ) is divided into three portions ( 472 , 474 , and 476 ) which arc added to separate PCRs using J region primers ( 401 ) and V region primers ( 404 , 406 , and 408 , respectively) to produce amplicons ( 420 , 422 and 424 , respectively).
- the latter amplicons are then combined ( 478 ) in secondary PCR ( 480 ) using P5 and P7 primers to prepare the templates ( 482 ) for bridge PCR and sequencing on an Illumina GA sequencer, or like instrument.
- Sequence reads of the invention may have a wide variety of lengths, depending in part on the sequencing technique being employed. For example, for some techniques, several trade-offs may arise in its implementation, for example, (i) the number and lengths of sequence reads per template and (ii) the cost and duration of a sequencing operation.
- sequence reads are in the range of from 20 to 400 nucleotides: in another embodiment, sequence reads are in a range of from 30 to 200 nucleotides; in still another embodiment, sequence reads arc in the range of from 30 to 120 nucleotides.
- 1 to 4 sequence reads are generated for determining the sequence of each clonotype; in another embodiment, 2 to 4 sequence reads arc generated for determining the sequence of each clonotype; and in another embodiment, 2 to 3 sequence reads are generated for determining the sequence of each clonotype.
- the numbers given are exclusive of sequence reads used to identify samples from different individuals.
- the lengths of the various sequence reads used in the embodiments described below may also vary based on the information that is sought to be captured by the read; for example, the starting location and length of a sequence read may be designed to provide the length of an NDN region as well as its nucleotide sequence: thus, sequence reads spanning the entire NDN region arc selected.
- sequences of clonotypes are determined in part by aligning sequence reads to one or more V region reference sequences and one or more J region reference sequences, and in part by base determination without alignment to reference sequences, such as in the highly variable NDN region.
- a variety of alignment algorithms may be applied to the sequence reads and reference sequences. For example, guidance for selecting alignment methods is available in Batzogloti, Briefings in Bioinformatics, 6: 6-22 (2005), which is incorporated by reference.
- a tree search algorithm is employed, e.g. as described generally in Gustield (cited above) and Cormen et al, Introduction to Algorithms, Third Edition (The MIT Press, 2009).
- an end of at least one forward read and an end of at least one reverse read overlap in an overlap region (e.g. 308 in FIG. 3B ), so that the bases of the reads arc in a reverse complementary relationship with one another.
- an overlap region e.g. 308 in FIG. 3B
- bases within such an overlap region arc determined, at least in part, from such a reverse complementary relationship.
- clonotypes of TCR ⁇ and IgH chains are determined by at least one sequence read starting in its J region and extending in the direction of its associated V region (referred to herein as “C read” ( 304 )) and at least one sequence read starting in its V region and extending in the direction of its associated J region (referred to herein as a “V read” ( 306 )).
- Overlap region ( 308 ) may or may not encompass the NDN region ( 315 ) as shown in FIG.
- Overlap region ( 308 ) may be entirely in the J region, entirely in the NDN region, entirely in the V region, or it may encompass a J region-NDN region boundary or a V region-NDN region boundary, or both such boundaries (as illustrated in FIG. 3B ).
- sequence reads arc generated by extending sequencing primers, e.g. ( 302 ) and ( 310 ) in FIG. 3B , with a polymerase in a sequencing-by-synthesis reaction, e.g. Metzger, Nature Reviews Genetics, 11: 31-46 (2010); Fuller et al, Nature Biotechnology, 27: 1013-1023 (2009).
- the binding sites for primers ( 302 ) and ( 310 ) are predetermined, so that they can provide a starting point or anchoring point for initial alignment and analysis, of the sequence reads.
- a C read is positioned so that it encompasses the D and/or NDN region of the TCR or IgH chain and includes a portion of the adjacent V region, e.g. as illustrated in FIGS. 3B and 3C .
- the overlap of the V read and the C read in the V region is used to align the reads with one another.
- such alignment of sequence reads is not necessary, e.g. with TCR ⁇ chains, so that a V read may only be long enough to identify the particular V region of a clonotype.
- Sequence read ( 330 ) is used to identify a V region, with or without overlapping another sequence read, and another sequence read ( 332 ) traverses the NDN region and is used to determine the sequence thereof. Portion ( 334 ) of sequence read ( 332 ) that extends into the V region is used to associate the sequence information of sequence read ( 332 ) with that of sequence read ( 330 ) to determine a clonotype.
- sequencing methods such as base-by-base approaches like the Solexa sequencing method, sequencing run time and reagent costs are reduced by minimizing the number of sequencing cycles in an analysis.
- FIG. 3C Sequence read ( 330 ) is used to identify a V region, with or without overlapping another sequence read, and another sequence read ( 332 ) traverses the NDN region and is used to determine the sequence thereof.
- Portion ( 334 ) of sequence read ( 332 ) that extends into the V region is used to associate the sequence information of sequence read ( 332 ) with that of sequence read ( 330
- amplicon ( 300 ) is produced with sample tag ( 312 ) to distinguish between clonotypes originating from different biological samples, e.g. different patients.
- Sample tag ( 312 ) may be identified by annealing a primer to primer binding region ( 316 ) and extending it ( 314 ) to produce a sequence read across tag ( 312 ), from which sample tag ( 312 ) is decoded.
- the IgH chain is more challenging to analyze than TCR ⁇ chain because of at least two factors: i) the presence of somatic mutations makes the mapping or alignment more difficult, and ii) the NDN region is larger so that it is often not possible to map a portion of the V segment to the C read.
- this problem is overcome by using a plurality of primer sets for generating V reads, which are located at different locations along the V region, preferably so that the primer binding sites are nonoverlapping and spaced apart, and with at least one primer binding site adjacent to the NDN region, e.g. in one embodiment from 5 to 50 bases from the V-NDN junction, or in another embodiment from 10 to 50 bases from the V-NDN junction.
- the redundancy of a plurality of primer sets minimizes the risk of failing to detect a clonotype due to a failure of one or two primers having binding sites affected by somatic mutations.
- the presence of at least one primer binding site adjacent to the NDN region makes it more likely that a V read will overlap with the C read and hence effectively extend the length of the C read.
- This allows for the generation of a continuous sequence that spans all sizes of NDN regions and that can also map substantially the entire V and J regions on both sides of the NDN region. Embodiments for carrying out such a scheme are illustrated in FIGS. 4A and 4D . In FIG.
- a sample comprising IgH chains ( 400 ) are sequenced by generating a plurality amplicons for each chain by amplifying the chains with a single set of J region primers ( 401 ) and a plurality (three shown) of sets of V region ( 402 ) primers ( 404 , 406 , 408 ) to produce a plurality of nested amplicons (e.g., 410 , 412 , 416 ) all comprising the same NDN region and having different lengths encompassing successively larger portions ( 411 , 413 , 415 ) of V region ( 402 ).
- nested amplicons e.g., 410 , 412 , 416
- the plurality of primer sets may be a number in the range of from 2 to 5. In another embodiment the plurality is 2-3; and still another embodiment the plurality is 3.
- the concentrations and positions of the primers in a plurality may vary widely. Concentrations of the V region primers may or may not be the same. In one embodiment, the primer closest to the NDN region has a higher concentration than the other primers of the plurality, e.g.
- a concentration ratio of 60:20:20 is used.
- One or more primers (e.g. 435 and 437 in FIG. 4B ) adjacent to the NDN region ( 444 ) may be used to generate one or more sequence reads (e.g. 434 and 436 ) that overlap the sequence read ( 442 ) generated by J region primer ( 432 ), thereby improving the quality of base calls in overlap region ( 440 ).
- Sequence reads from the plurality of primers may or may not overlap the adjacent downstream primer binding site and/or adjacent downstream sequence read.
- sequence reads proximal to the NDN region may be used to identify the particular V region associated with the clonotype.
- sequence reads proximal to the NDN region may be used to identify the particular V region associated with the clonotype.
- Such a plurality of primers reduces the likelihood of incomplete or failed amplification in case one of the primer binding sites is hypermutated during immunoglobulin development. It also increases the likelihood that diversity introduced by hypermutation of the V region will be capture in a clonotype sequence.
- a secondary PCR may be performed to prepare the nested amplicons for sequencing, e.g.
- amplicons 420 , 422 , and 424 ), which may be distributed as single molecules on a solid surface, where they are further amplified by bridge PCR, or like technique.
- NDN regions can be improved by using the codon stricture of the flanking J and V regions, as illustrated in FIG. 4E .
- “codon.structure” means the codons of the natural reading frame of segments of TCR or BCR transcripts or genes outside of the NDN regions, e.g. the V region, J region, or the like.
- amplicon ( 450 ) which is an enlarged view of the amplicon of FIG. 4B , is shown along with the relative positions of C read ( 442 ) and adjacent V read ( 434 ) above and the codon structures ( 452 and 454 ) of V region ( 430 ) and J region ( 446 ), respectively, below.
- bases in NDN region ( 456 ) are called (or identified) one base at a time moving from J region ( 446 ) toward V region ( 430 ) and in the opposite direction from V region ( 430 ) toward J region ( 446 ) using sequence reads ( 434 ) and ( 442 ).
- sequence reads ( 434 ) and ( 442 ) are sequence reads ( 434 ) and ( 442 ).
- the variants generated somatically only ones expressed are those whose J region and V region codon frames are in-frame with one another and remain in-frame through the NDN region. (Here the correct frames of the V and J regions arc determined from reference sequences). If an out-of-frame sequence is identified based one or more low quality base calls, the corresponding clonotype is flagged for re-evaluation or as a potential disease-related anomaly. If the sequence identified is in-frame and based on high quality base calls, then there is greater confidence that the corresponding clonotype has been correctly called.
- the invention includes a method of determining V(D)J-based clonotypes from bidirectional sequence reads comprising the steps of: (a) generating at least one J region sequence read that begins in a J region and extends into an NDN region and at least one V region sequence read that begins in the V regions and extends toward the NDN region such that the J region sequence read and the V region sequence read are overlapping in an overlap region, and the J region and the V region each have a codon structure; (b) determining whether the codon structure of the J region extended into the NDN region is in frame with the codon structure of the V region extended toward the NDN region.
- the step of generating includes generating at least one V region sequence read that begins in the V region and extends through the NDN region to the J region, such that the J region sequence read and the V region sequence read arc overlapping in an overlap region.
- IgH-based clonotypes that have undergone somatic hypermutation are determined as follows: A somatic mutation is defined as a sequenced base that is different from the corresponding base of a reference sequence (of the relevant segment, usually V, J or C) and that is present in a statistically significant number of reads.
- C reads may be used to find somatic mutations with respect to the mapped J segment and likewise V reads for the V segment. Only pieces of the C and V reads are used that are either directly mapped to J or V segments or that are inside the clonotype extension up to the NDN boundary.
- the NDN region is avoided and the same ‘sequence information’ is not used for mutation finding that was previously used for clonotype determination (to avoid erroneously classifying as mutations nucleotides that arc really just different recombined NDN regions).
- the mapped segment (major allele) is used as a scaffold and all reads are considered which have mapped to this allele during the read mapping phase.
- Each position of the reference sequences where at least one read has mapped is analyzed for somatic mutations.
- the criteria for accepting a non-reference base as a valid mutation include the following: 1) at least N reads with the given mutation base, 2) at least a given fraction N/M reads (where M is the total number of mapped reads at this base position) and 3) a statistical cut based on the binomial distribution, the average Q score of the N reads at the mutation base as well as the number (M-N) of reads with a non-mutation base.
- the above parameters arc selected so that the false discovery rate of mutations per clonotype is less than 1 in 1000, and more preferably, less than 1 in 10000.
- Phylogenic Clonotypes (Clans).
- cancers such as lymphoid neoplasms
- a single lymphocyte progenitor may give rise to many related lymphocyte progeny, each possessing and/or expressing a slightly different TCR or BCR, and therefore a different clonotype, due to cancer-related somatic mutation(s), such as base substitutions, aberrant rearrangements, or the like.
- phylogenic clones Cells producing such clonotypes are referred to herein as phylogenic clones, and a set of such related clones are referred to herein as a “clan.”
- clonotypes of phylogenic clones are referred to as phylogenic clonotypes and a set of phylogenic clonotypes may be referred to as a clan of clonotypes.
- methods of the invention comprise monitoring the frequency of a clan of clonotypes (i.e., the sum of frequencies of the constituent phylogenic clonotypes of the clan), rather than a frequency of an individual clonotype.
- Phylogenic clonotypes may be identified by one or more measures of relatedness to a parent clonotype. In one embodiment, phylogenic clonotypes may be grouped into the same clan by percent homology, as described more fully below. In another embodiment, phylogenic clonotypes are identified by common usage of V regions, J regions, and/or NDN regions.
- a clan may be defined by clonotypes having common J and ND regions but different V regions; or it may be defined by clonotypes having the same V and J regions (including identical base substitutions mutations) but with different NDN regions; or it may be defined by a clonotype-that has undergone one or more insertions and/or deletions of from 1-10 bases, or from 1-5 bases, or from 1-3 bases, to generate clan members.
- members of a clan are determined as follows.
- Clonotypes are assigned to the same clan if they satisfy the following criteria: i) they are mapped to the same V and J reference segments, with the mappings occurring at the same relative positions in the clonotype sequence, and ii) their NDN regions arc substantially identical. “Substantial” in reference to clan membership means that some small differences in the NDN region are allowed because somatic mutations may have occurred in this region. Preferably, in one embodiment, to avoid falsely calling a mutation in the NDN region, whether a base substitution is accepted as a cancer-related mutation depends directly on the size of the NDN region of the clan.
- a method may accept a clonotype as a clan member if it has a one-base difference from clan NDN sequence(s) as a cancer-related mutation if the length of the clan NDN sequence(s) is m nucleotides or greater, e.g. 9 nucleotides or greater, otherwise it is not accepted, or if it has a two-base difference from clan NDN sequence(s) as cancer-related mutations if the length of the clan NDN sequence(s) is n nucleotides or greater, e.g.
- members of a clan are determined using the following criteria: (a) V read maps to the same V region, (b) C read maps to the same J region, (c) NDN region substantially identical (as described above), and (d) position of NDN region between V-NDN boundary and J-NDN boundary is the same (or equivalently, the number of downstream base additions to D and the number of upstream base additions to D are the same).
- Clonotypes of a single sample may be grouped into clans and, clans from successive samples acquired at different times may be compared with one another.
- clans containing clonotypes correlated with a disease are identified from clonotypes of each sample and compared with that of the immediately previous sample to determine disease status, such as, continued remission, incipient relapse, evidence of further clonal evolution, or the like.
- PCR error is concentrated in some bases that were mutated in the early cycles of PCR. Sequencing error is expected to be distributed in many bases even though it is totally random as the error is likely to have some systematic biases. It is assumed that some bases will have sequencing error at a higher rate, say 5% (5 fold the average). Given these assumptions, sequencing error becomes the dominant type of error. Distinguishing PCR errors from the occurrence of highly related clonotypes will play a role in analysis. Given the biological significance to determining that there are two or more highly related clonotypes, a conservative approach to making such calls is taken.
- the detection of enough of the minor clonotypes so as to be sure with high confidence (say 99.9%) that there are more than one clonotype is considered. For example of clonotypes that are present at 100 copies/1,000,000, the minor variant is detected 14 or more times for it to be designated as an independent clonotype. Similarly, for clonotypes present at 1,000 copies/1,000,000 the minor variant can be detected 74 or more times to be designated as an independent clonotype. This algorithm can be enhanced by using the base quality score that is obtained with each sequenced base.
- the quality score can be used to decide the number of reads that need to be present to call an independent clonotype.
- the median quality score of the specific base in all the reads can be used, or more rigorously, the likelihood of being an error can be computed given the quality score of the specific base in each read, and then the probabilities can be combined (assuming independence) to estimate the likely number of sequencing error for that base.
- TCR ⁇ chains are analyzed.
- the analysis includes amplification, sequencing, and analyzing the TCR ⁇ sequences.
- One primer is complementary to a common sequence in C ⁇ 1 and C ⁇ 2, and there arc 34 V primers capable of amplifying all 48 V segments.
- C ⁇ 1 or C ⁇ 2 differ from each other at position 10 and 14 from the J/C junction.
- the primer for C ⁇ 1 and C ⁇ 2 ends at position 16 by and has no preference for C ⁇ 1 or C ⁇ 2.
- the 34 V primers are modified from an original set of primers disclosed in Van Dongen et al, U.S. patent publication 2006/0234234, which is incorporated herein by reference.
- the modified primers arc disclosed in Faham et al, U.S. patent publication 2010/0151471, which is also incorporated herein by reference.
- the Illumina Genome Analyzer is used to sequence the amplicon produced by the above primers.
- a two-stage amplification is performed on messenger RNA transcripts ( 200 ), as illustrated in FIGS. 2A-2B , the first stage employing the above primers and a second stage to add common primers for bridge amplification and sequencing.
- a primary PCR is performed using on one side a 20 by primer ( 202 ) whose 3′ end is 16 bases from the J/C junction ( 204 ) and which is perfectly complementary to C ⁇ 1( 203 ) and the two alleles of C ⁇ 2.
- primer set ( 212 ) contains primer sequences complementary to the different V region sequences ( 34 in one embodiment). Primers of set ( 212 ) also contain a non-complementary tail ( 214 ) that produces amplicon ( 216 ) having primer binding site ( 218 ) specific for P7 primers ( 220 ).
- amplicon ( 216 ) is formed that contains the highly diverse portion of the J(D)V region ( 206 , 208 , and 210 ) of the mRNA transcripts and common primer binding sites ( 203 and 218 ) for a secondary amplification to add a sample tag ( 221 ) and primers ( 220 and 222 ) for cluster formation by bridge PCR.
- a primer 222 in FIG.
- C10-17-P5 2B and referred to herein as “C10-17-P5”
- C10-17-P5 has at its 3′end the sequence of the 10 bases closest to the J/C junction, followed by 17 bp with the sequence of positions 15-31 from the J/C junction, followed by the PS sequence ( 224 ), which plays a role in cluster formation by bridge PCR in Solexa sequencing.
- the C10-17-P5 primer ( 222 ) anneals to the template generated from the first PCR, a 4 by loop (position 11-14) is created in the template, as the primer hybridizes to the sequence of the 10 bases closest to the J/C junction and bases at positions 15-31 from the J/C junction.
- the looping of positions 11-14 eliminates differential amplification of templates carrying C ⁇ 1 or C ⁇ 2.
- C10-17-P5 primer can be HPLC purified in order to ensure that all the amplified material has intact ends that can be efficiently utilized in the cluster formation.
- the length of the overhang on the V primers ( 212 ) is preferably 14 bp.
- the primary PCR is helped with a shorter overhang ( 214 ).
- the overhang in the V primer is used in the primary PCR as long as possible because the secondary PCR is priming from this sequence.
- a minimum size of overhang ( 214 ) that supports an efficient secondary PCR was investigated. Two series of V primers (for two different V segments) with overhang sires from 10 to 30 with 2 by steps were made. Using the appropriate synthetic sequences, the first PCR was performed with each of the primers in the series and gel electrophoresis was performed to show that all amplified.
- the primary PCR uses 34 different V primers ( 212 ) that anneal to V region ( 206 ) of RNA templates ( 200 ) and contain a common 14 by overhang on the 5′ tail.
- the 14 by is the partial sequence of one of the Illumina sequencing primers (termed the Read 2 primer).
- the secondary amplification primer ( 220 ) on the same side includes P7 sequence, a tag ( 221 ), and Read 2 primer sequence ( 223 ) (this primer is called Read2_tagX_P7).
- the P7 sequence is used for cluster formation.
- Read 2 primer and its complement are used for sequencing the V segment and the tag respectively.
- a set of 96 of these primers with tags numbered 1 through 96 arc created (sec below). These primers arc HPLC purified in order to ensure that all the amplified material has intact ends that can be efficiently utilized in the cluster formation.
- the second stage primer, C-10-17-P5 ( 222 , FIG. 2B ) has interrupted homology to the template generated in the first stage PCR.
- the efficiency of amplification using this primer has been validated.
- An alternative primer to C-10-17-P5, termed CsegP5 has perfect homology to the first stage C primer and a 5′ tail carrying P5.
- PCR templates was compared by performing real time PCR. In several replicates, it was found that PCR using the C-10-17-P5 primer had little or no difference in efficiency compared with PCR using the CsegP5 primer.
- Amplicon ( 300 ) resulting from the 2-stage amplification illustrated in FIGS. 2A-2B has the structure typically used with the Illumina sequencer as shown in FIG. 3A .
- Two primers that anneal to the outmost part of the molecule, Illumina primers P5 and P7 are used for solid phase amplification of the molecule (cluster formation).
- Three sequence reads arc done per molecule. The first read of 100 by is done with the C′ primer, which has a melting temperature that is appropriate for the Illumina sequencing process. The second read is 6 by long only and is solely for the purpose of identifying the sample tag. It is generated using a tag primer provided by the manufacturer (Illumina). The final read is the Read 2 primer, also provided by the manufacturer (Illumina). Using this primer, a 100 by read in the V segment is generated starting with the 1st PCR V primer sequence.
- “Aligning” means a method of comparing a test sequence, such as a sequence read, to one or more reference sequences to determine which reference sequence or which portion of a reference sequence is closest based on some sequence distance measure.
- An exemplary method of aligning nucleotide sequences is the Smith Waterman algorithm.
- Distance measures may include Hamming distance, Levenshtein distance, or the like. Distance measures may include a component related to the quality values of nucleotides of the sequences being compared.
- Amplicon means the product of a polynucleotide amplification reaction; that is, a clonal population of polynucleotides, which may be single stranded or double stranded, which are, replicated from one or more starting sequences.
- the one or more starting sequences may be one or more copies of the same sequence, or they may be a mixture of different sequences.
- amplicons are formed by the amplification of a single starting sequence. Amplicons may be produced by a variety of amplification reactions whose products comprise replicates of the one or more starting, or target, nucleic acids.
- amplification reactions producing amplicons arc “template-driven” in that base pairing of reactants, either nucleotides or oligonucleotides, have complements in a template polynucleotide that are required for the creation of reaction products.
- template-driven reactions are primer extensions with a nucleic acid polymerase or oligonucleotide ligations with a nucleic acid ligase.
- Such reactions include, but arc not limited to, polymerase chain reactions (PCRs), linear polymerase reactions, nucleic acid sequence-based amplification (NASBAs), rolling circle amplifications, and the like, disclosed in the following references that are incorporated herein by reference: Mullis et al, U.S. Pat.
- amplicons of the invention are produced by PCRs.
- An amplification reaction may be a “real-time” amplification if a detection chemistry is available that permits a reaction product to be measured as the amplification reaction progresses, e.g. “real-time PCR” described below, or “real-time NASBA” as described in Leone et al, Nucleic Acids Research, 26: 2150-2155 (1998), and like references.
- the term “amplifying” means performing an amplification reaction.
- a “reaction mixture” means a solution containing all the necessary reactants for performing a reaction, which may include, but not be limited to, buffering agents to maintain pH at a selected level during a reaction, salts, co-factors, scavengers, and the like.
- “Clonality” as used herein means a measure of the degree to which the distribution of clonotype abundances among clonotypes of a repertoire is skewed to a single or a few clonotypes. Roughly, clonality is an inverse measure of clonotype diversity. Many measures or statistics are available from ecology describing species-abundance relationships that may be used for clonality measures in accordance with the invention, e.g. Chapters 17 & 18, in Pielou, An Introduction to Mathematical Ecology, (Wiley-Interscience, 1969).
- a clonality measure used with the invention is a function of a clonotype profile (that is, the number of distinct clonotypes detected and their abundances), so that after a clonotype profile is measured, clonality may be computed from it to give a single number.
- a clonality measure is Simpson's measure, which is simply the probability that two randomly drawn clonotypes will be the same.
- Other clonality measures include information, based measures and McIntosh's diversity index, disclosed in Pielou (cited above).
- “Clonotype” means a recombined nucleotide sequence of a T cell or B cell encoding a T cell receptor (TCR) or B cell receptor (BCR); or a portion thereof.
- TCR T cell receptor
- BCR B cell receptor
- a collection of all the distinct clonotypes of a population of lymphocytes of an individual is a repertoire of such population, e.g. Arstila et al, Science, 286: 958-961 (1999); Yassai et al, Immunogenetics, 61: 493-502 (2009); Kedzierska et al, Mol. Immunol., 45(3): 607-618 (2008); and the like.
- clonotype profile is a tabulation of clonotypes of a sample of T cells and/or B cells (such as a peripheral blood sample containing such cells) that includes substantially all of the repertoire's clonotypes and their relative abundances.
- “Clonotype profile,” “repertoire profile,” and “repertoire” are used herein interchangeably. (That is, the term “repertoire,” as discussed more fully below, means a repertoire measured from a sample of lymphocytes).
- clonotypes comprise portions of an immunoglobulin heavy chain (IgH) or a TCR ⁇ chain.
- clonotypes may be based on other recombined molecules, such as immunoglobulin chains or TCR ⁇ chains, or portions thereof.
- Coding means treating two candidate clonotypes with sequence differences as the same by determining that such differences are due to experimental or measurement error and not due to genuine biological differences.
- a sequence of a higher frequency candidate clonotype is compared to that of a lower frequency candidate clonotype and if predetermined criteria are satisfied then the number of lower frequency candidate clonotypes is added to that of the higher frequency candidate clonotype and the lower frequency candidate clonotype is thereafter disregarded. That is, the read counts associated with the lower frequency candidate clonotype are added to those of the higher frequency candidate clonotype.
- CDRs complementarity determining regions
- T cell receptors and immunoglobulins each have three CDRs: CDR1 and CDR2 are found in the variable (V) domain, and CDR3 includes some of V, all of diverse (D) (heavy chains only) and joint (J), and some of the constant (C) domains.
- Data structure means an organization of information, usually in a computer or memory device, for better algorithm efficiency.
- Exemplary data structures include queues, stacks, linked lists, heaps, hash tables, arrays, trees, and the like.
- Data structures may have substructures that correspond to units of information or to subsets of related information. For example, arrays have rows and columns of entries; trees have nodes, branches, subtrees, and leaves; or the like.
- a data structure used herein is a sequence tree, an array or a hash table.
- Lymphoid neoplasm means an abnormal proliferation of lymphocytes that may be malignant or non-malignant.
- a lymphoid cancer is a malignant lymphoid neoplasm.
- Lymphoid neoplasms arc the result of, or are associated with, lymphoproliferative disorders, including but not limited to, follicular lymphoma, chronic lymphocytic leukemia (CLL), acute lymphocytic leukemia (ALL), hairy cell leukemia, lymphomas, multiple myeloma, post-transplant lymphoproliferative disorder, mantle lymphoma (MCL), diffuse large B cell lymphoma (DLBCL), T cell lymphoma, or the like; e.g.
- Percent homologous “percent identical,” or like terms used in reference to the comparison of a reference sequence and another sequence (“comparison sequence”) mean that in an optimal alignment between the two sequences, the comparison sequence is identical to the reference sequence in a number of subunit positions equivalent to the indicated percentage, the subunits being nucleotides for polynucleotide comparisons or amino acids for polypeptide comparisons.
- an “optimal alignment” of sequences being compared is one that maximizes matches between subunits and minimizes the number of gaps employed in constructing an alignment. Percent identities may be determined with commercially available implementations of algorithms, such as that described by Needleman and Wunsch, J. Mol.
- a polynucleotide having a nucleotide sequence at least 95 percent identical to a reference nucleotide sequence up to five percent of the nucleotides in the reference sequence may be deleted or substituted with another nucleotide, or a number of nucleotides up to five percent of the total number of nucleotides in the reference sequence may be inserted into the reference sequence.
- PCR Polymerase chain reaction
- PCR is a reaction for making multiple copies or replicates of a target nucleic acid flanked by primer binding sites, such reaction comprising one or more repetitions of the following steps: (i) denaturing the target nucleic acid, (ii) annealing primers to the primer binding sites, and (iii) extending the primers by a nucleic acid polymerase in the presence of nucleoside triphosphates.
- the reaction is cycled through different temperatures optimized for each step in a thermal cycler instrument.
- a double stranded target nucleic acid may be denatured at a temperature >90° C., primers annealed at a temperature in the range 50-75° C., and primers extended at a temperature in the range 72-78° C.
- PCR encompasses derivative forms of the reaction, including but not limited to, RT-PCR, real-time PCR, nested PCR, quantitative PCR, multiplexed PCR, and the like. Reaction volumes range from a few hundred nanoliters, e.g. 200 nL, to a few hundred ⁇ L, e.g. 200 ⁇ L.
- Reverse transcription PCR or “RT-PCR,” means a PCR that is preceded by a reverse transcription reaction that converts a target RNA to a complementary single stranded DNA, which is then amplified, e.g. Tecott et al, U.S. Pat. No. 5,168,038, which patent is incorporated herein by reference.
- Real-time PCR means a PCR for which the amount of reaction product, i.e. amplicon, is monitored as the reaction proceeds.
- Nested PCR means a two-stage PCR wherein the amplicon of a first PCR becomes the sample for a second PCR using a new set of primers, at least one of which binds to an interior location of the first amplicon.
- initial primers in reference to a nested amplification reaction mean the printers used to generate a first amplicon
- secondary primers mean the one or more primers used to generate a second, or nested, amplicon.
- Multiplexed PCR means a PCR wherein multiple target sequences (or a single target sequence and one or more reference sequences) are simultaneously carried out in the same reaction mixture, e.g. Bernard ct al, Anal. Biochem., 273: 221-228 (1999) (two-color real-time PCR). Usually, distinct sets of primers are employed for each sequence being amplified. Typically, the number of target sequences in a multiplex PCR is in the range of from 2 to 50, or from 2 to 40, or from 2 to 30. “Quantitative PCR” means a PCR designed to measure the abundance of one or more specific target sequences in a sample or specimen. Quantitative PCR includes both absolute quantitation and relative quantitation of such target sequences.
- Quantitative measurements are made using one or more reference sequences or internal standards that may be assayed separately or together with a target sequence.
- the reference sequence may be endogenous or exogenous to a sample or specimen, and in the latter case, may comprise one or more competitor templates.
- Typical endogenous reference sequences include segments of transcripts of the following genes: ⁇ -actin, GAPDH, ⁇ 2 -microglobulin, ribosomal RNA, and the like.
- Primer means an oligonucleotide, either natural or synthetic that is capable, upon forming a duplex with a polynucleotide template, of acting as a point of initiation of nucleic acid synthesis and being extended from its 3′ end along the template so that an extended duplex is formed.
- Extension of a primer is usually carried out with a nucleic acid polymerase, such as a DNA or RNA polymerase.
- the sequence of nucleotides added in the extension process is determined by the sequence of the template polynucleotide.
- primers arc extended by a DNA polymerase.
- Primers usually have a length in the range of from 14 to 40 nucleotides, or in the range of from 18 to 36 nucleotides.
- Primers are employed in a variety of nucleic amplification reactions, for example, linear amplification reactions using a single primer, or polymerase chain reactions, employing two or more primers.
- Guidance for selecting the lengths and sequences of primers for particular applications is well known to those of ordinary skill in the art, as evidenced by the following references that are incorporated by reference: Dieffenbach, editor, PCR Primer: A Laboratory Manual, 2 nd Edition (Cold Spring Harbor Press, New York, 2003).
- Quality score means a measure of the probability that a base assignment at a particular sequence location is correct.
- a variety methods are well known to those of ordinary skill for calculating quality scores for particular circumstances, such as, for bases called as a result of different sequencing chemistries, detection systems, base-calling algorithms, and so on.
- quality score values arc monotonically related to probabilities of correct base calling. For example, a quality score, or Q, of 10 may mean that there is a 90 percent chance that a base is called correctly, a Q of 20 may mean that there is a 99 percent chance that a base is called correctly, and so on.
- average quality scores decrease as a function of sequence read length, so that quality scores at the beginning of a sequence read are higher than those at the end of a sequence read, such declines being due to phenomena such as incomplete extensions, carry forward extensions, loss of template, loss of polymerase, capping failures, deprotection failures, and the like.
- “Repertoire”, or “immune repertoire” means a set of distinct recombined nucleotide sequences that encode T cell receptors (TCRs) or B cell receptors (BCRs), or fragments thereof, respectively; in a population Oft lymphocytes of an individual, wherein the nucleotide sequences of the set have a one-to-one correspondence with distinct lymphocytes or their subpopulations for substantially all of the lymphocytes of the population.
- a population of lymphocytes from which a repertoire is determined is taken from one or more tissue samples, such as one or more blood samples.
- clonotypes of a repertoire comprises any segment of nucleic acid common to a T cell or a B cell population which has undergone somatic recombination during the development of TCRs or BCRs, including normal or aberrant (e.g. associated with cancers) precursor molecules thereof, including, but not limited to, any of the following: an immunoglobulin heavy chain (IgH) or subsets thereof (e.g. an IgH variable region, CDR3 region, or the like), incomplete IgH molecules, an immunoglobulin light chain or subsets thereof (e.g.
- IgH immunoglobulin heavy chain
- subsets thereof e.g. an IgH variable region, CDR3 region, or the like
- incomplete IgH molecules e.g.
- nucleic acid segments defining clonotypes of a repertoire arc selected so that their diversity (i.e.
- the number of distinct nucleic acid sequences in the set is large enough so that substantially every T cell or B cell or clone thereof in an individual carries a unique nucleic acid sequence of such repertoire. That is, in accordance with the invention, a practitioner may select for defining clonotypes a particular segment or region of recombined nucleic acids that encode TCRs or BCRs that do not reflect the full diversity of a population of T cells or B cells; however, preferably, clonotypes are defined so that they do reflect the diversity of the population of T cells and/or B cells from which they are derived. That is, preferably each different clone of a sample has different clonotype.
- the population of lymphocytes corresponding to a repertoire may be circulating B cells, or may be circulating T cells, or may be subpopulations of either of the foregoing populations, including but not limited to, CD4+ T cells, or CD8+ T cells, or other subpopulations defined by cell surface markers, or the like.
- Such subpopulations may be acquired by taking samples from particular tissues, e.g.
- the population of lymphocytes corresponding to a repertoire may be derived from disease tissues, such as a tumor tissue, an infected tissue, or the like.
- a repertoire comprising human TCR ⁇ chains or fragments thereof comprises a number of distinct nucleotide sequences in the range of from 0.1 ⁇ 10 6 to 1.8 ⁇ 10 6 , or in the range of from 0.5 ⁇ 10 6 to 1.5 ⁇ 10 6 , or in the range of from 0.8 ⁇ 10 6 to 1.2 ⁇ 10 6 .
- a repertoire comprising human comprises a number of distinct nucleotide sequences in the range of from 0.1 ⁇ 10 6 to 1.8 ⁇ 10 6 , or in the range of from 0.5 ⁇ 10 6 to 1.5 ⁇ 10 6 , or in the range of from 0.8 ⁇ 10 6 to 1.2 ⁇ 10 6 .
- IgH chains or fragments thereof comprises a number of distinct nucleotide sequences in the range of from 0.1 ⁇ 10 6 to 1.8 ⁇ 10 6 , or in the range of from 0.5 ⁇ 10 6 to 1.5 ⁇ 10 6 , or in the range of from 0.8 ⁇ 10 6 to 1.2 ⁇ 10 6 .
- a repertoire of the invention comprises a set of nucleotide sequences encoding substantially all segments of the V(D)J region of an IgH chain.
- “substantially all” as used herein means every segment having a relative abundance of 0.001 percent or higher; or in another aspect, “substantially all” as used herein means every segment having a relative abundance of 0.0001 percent or higher.
- a repertoire of the invention comprises a set of nucleotide sequences that encodes substantially all segments of the V(D)J region of a TCR chain.
- a repertoire of the invention comprises a set of nucleotide sequences having lengths in the range of from 25-200 nucleotides and including segments of the V, D, and J regions of a TCR ⁇ chain.
- a repertoire of the invention comprises a set of nucleotide sequences having lengths in the range of from 25-200 nucleotides and including segments of the V, D, and J regions of an IgH chain.
- a repertoire of the invention comprises a number of distinct nucleotide sequences that is substantially equivalent to the number of lymphocytes expressing a distinct IgH chain. In another embodiment, a repertoire of the invention comprises a number of distinct nucleotide sequences that is substantially equivalent to the number of lymphocytes expressing a distinct TCR ⁇ chain. In still another embodiment, “substantially equivalent” means that with ninety-nine percent probability a repertoire of nucleotide sequences will include a nucleotide sequence encoding an IgH or TCR ⁇ or portion thereof carried or expressed by every lymphocyte of a population of an individual at a frequency of 0.001 percent or greater.
- “substantially equivalent” means that with ninety-nine percent probability a repertoire of nucleotide sequences will include a nucleotide sequence encoding an IgH or TCR ⁇ or portion thereof carried or expressed by every lymphocyte present at a frequency of 0.0001 percent or greater.
- the sets of clonotypes described in the foregoing two sentences are sometimes referred to herein as representing the “full repertoire” of IgH and/or TCR ⁇ sequences.
- a sufficiently large sample of lymphocytes is obtained so that such profile provides a reasonably accurate representation of a repertoire for a particular application.
- samples comprising from 10 5 to 10 7 lymphocytes are employed, especially when obtained from peripheral blood samples of from 1-10 mL.
- Sequence read means a sequence of nucleotides determined from a Sequence or stream of data generated by a sequencing technique, which determination is made, for example, by means; of base-calling software associated with the technique, e.g. base-calling software from a commercial provider of a DNA sequencing platform.
- a sequence read usually includes quality scores for each nucleotide in the sequence.
- sequence reads are made by extending a primer along a template nucleic acid, e.g. with a DNA polymerase or a DNA ligase. Data is generated by recording signals, such as optical, chemical (e.g. pH change), or electrical signals, associated with such extension. Such initial data is converted into a sequence read.
- a tree data structure of the invention is a rooted directed tree comprising nodes and edges that do not include cycles, or cyclical pathways. Edges from nodes of tree data structures of the invention are usually ordered. Nodes and/or edges arc structures that may contain, or be associated with, a value. Each node in a tree has zero or more child nodes, which by convention are shown below it in the tree. A node that has a child is called the child's parent node. A node has at most one parent. Nodes that do not have any children are called leaf nodes. The topmost node in a tree is called the root node. Being the topmost node, the root node will not have parents. It is the node at which operations on the tree commonly begin (although some algorithms begin with, the leaf nodes and work up ending at the root). All other nodes can be reached from it by following edges or links.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Organic Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biophysics (AREA)
- Biotechnology (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Genetics & Genomics (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medical Informatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Immunology (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Biochemistry (AREA)
- Microbiology (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Cell Biology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
Description
- This application is a continuation-in-part of copending U.S. application Ser. No. 13/100,389 filed 4 May 2011, which is a continuation-in-part of copending U.S. application Ser. No. 12/615,263 filed 9 Nov. 2009, and this application claims benefit of U.S. provisional applications 61/446,822 filed 25 Feb. 2011 and 61/455,743 filed 25 Oct. 2010, the foregoing applications being incorporated by reference in their entireties.
- Analysis of biological or medical samples often requires the determination of nucleic acid sequences of large and complex populations of DNA and/or RNA, e.g. Gloor et al, PLoS ONE 5(10): el 5406 (2010); Petrosino et al, Clinical Chemistry, 55(5): 856-866 (2009); Arstila et al, Science, 286: 958-961 (1999). In particular, profiles of nucleic acids encoding immune molecules, such as T cell or B cell receptors, or their components, contain a wealth of information on the state of health or disease of an organism, so that the use of such profiles as diagnostic or prognostic indicators has been proposed for a wide variety of conditions, e.g. Faham and Willis, U.S. patent publication 2010/0151471; Freeman et al, Genome Research, 19: 1817-1824 (2009); Boyd et al, Sci. Transl. Med., 1(12): 12ra23 (2009); He et al, Oncotarget (Mar. 8, 2011). Such sequence-based profiles arc capable of much greater sensitivity than approaches based on size distributions of amplified target nucleic acids, sequence sampling by microarrays, hybridization kinetics curves from PCR amplicons, or other approaches, e.g. Morley et al, U.S. Pat. No. 5,418,134; van Dongen et al, Leukemia, 17: 2257-2317 (2003); Ogle ct al, Nucleic Acids Research, 31: e139 (2003); Wang et al, BMC Genomics, 8: 329 (2007); Baum et al, Nature Methods, 3(11): 895,901 (2006). However, the efficient determination of clonotypes and clonotype profiles from sequence data poses challenges because of the size populations to be analyzed, the limited predictability of natural variability in the sequences extracted from samples, and noise introduced into the data by a host of sample preparation and measurement steps, e.g. Warren et al, Genome Research, 21(5): 790-797 (2011).
- In view of the potential importance of clonotype profiles for diagnostic and prognostic applications, it would be advantageous to many fields in medicine and biology, if methods were available for overcoming drawbacks of current methodologies for determining clonotypes and clonotype profiles from sequence data.
- The present invention is drawn to methods for producing sequence-based profiles of complex nucleic acid populations. The invention is exemplified in a number of implementations and applications, some of which are summarized below and throughout the specification.
- In one aspect, the invention is directed to a method of forming a clonotype profile comprising the steps of (a) forming a data structure of somatically recombined immune molecules from sequence reads thereof, (b) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype, (c) removing the coalesced candidate clonotype from the data structure, and (d) repeating steps (b) and (c) until a clonotype profile is formed. In one aspect, such data structure is a sequence tree.
- In another aspect, the invention is directed to a method of generating a clonotype profile from an individual comprising the steps of: (a) spatially isolating individual molecules derived from a nucleic acid sample from T-cells and/or B-cells of the individual; (b) sequencing said spatially isolated individual molecules to produce a plurality of sequence reads each having portions old V region, an NDN region and a J region wherein for each such molecule at least one sequence lead encompasses a portion of the NDN region (c) forming from sequence reads encompassing at least a portion of an NDN region a data structure having elements representing candidate clonotypes, each element and its corresponding candidate clonotype having a frequency; (d) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype having a sequence of the highest frequency candidate clonotype; (e) removing elements corresponding to the coalesced candidate clonotypes from the data structure; and (f) repeating steps (d) and (e) until clonotypes have been formed from all non-singleton lesser frequency candidate clonotypes, thereby generating the clonotype profile. In one aspect, the data structure is a sequence tree and the elements corresponding to candidate clonotypes are leaves of the sequence tree.
- In another aspect, the invention provides a method for determining clonotypes of an immune repertoire comprising the steps: (a) providing a set of sequence reads from a repertoire of recombined immune molecules each having a V region, an NDN region and a J region wherein for each such molecule at least one sequence read encompasses at least a portion of the NDN region of such molecule; (b) forming from sequence reads encompassing at least a portion of an NDN region a data structure having elements representing candidate clonotypes, each leaf and its corresponding candidate clonotype having a frequency; (c) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype having a sequence of the highest frequency candidate clonotype; (d) removing elements corresponding to the coalesced candidate clonotypes from the data structure; and (e) repeating steps (c) and (d) until a highest frequency of a lesser frequency candidate clonotype is below a predetermined stopping value. In one aspect, the data structure is a sequence tree and the elements corresponding to candidate clonotypes are leaves of the sequence tree.
- In still another aspect, the invention is directed to a method of determining clonotypes of an immune repertoire comprising the steps: (a) providing a set of sequence reads from a repertoire of recombined immune molecules each having portions of a V region, an NDN region and a J region wherein for each such molecule at least one sequence read encompasses a portion of the NDN region of such molecule; (b) forming from sequence reads encompassing portions of NDN regions a sequence tree having leaves representing candidate clonotypes, each leaf and its corresponding candidate clonotype having a frequency; (c) selecting a highest frequency candidate clonotype and identifying all said lesser frequency candidate clonotypes having a sequence difference therewith less than a predetermined difference value to form a coalescence subset; (d) coalescing with the highest frequency candidate clonotype any lesser frequency candidate clonotypes in the coalescence subset whenever such lesser frequency is below a predetermined frequency value to form a clonotype having a sequence of the highest frequency candidate clonotype; (e) removing leaves corresponding to the coalesced candidate clonotypes from the sequence tree; and (f) repeating steps (c) and (c) until clonotypes have been formed from all non-singleton lesser frequency candidate clonotypes.
- The present invention provides methods for determining clonotypes and clonotype profiles from large sets of sequence data obtained by high throughput sequencing Of somatically recombined immune molecules. In one aspect, the invention implements the above methods by using one or more sequence trees for efficiently carrying out sequence comparisons necessary for distinguishing candidate clonotypes with genuine sequence differences from those with experimental or measurement errors, such as sequencing errors.
- These above-characterized aspects, as well as other aspects, of the present invention are exemplified in a number of illustrated implementations and applications, some of which are shown in the figures and characterized in the claims section that follows. However, the above summary is not intended to describe each illustrated embodiment or every implementation of the invention.
- The novel features of the invention are set forth with particularity in the appended claims. A better understanding of the features and advantages of the present invention is obtained by reference to the following detailed description that sets forth illustrative embodiments, in which the principles of the invention are utilized, and the accompanying drawings of which:
-
FIG. 1A illustrates an IgH transcript and sources of natural variability within it. -
FIG. 1B illustrates concepts of clonotypes in sequence space and distances between closely related clonotypes. -
FIG. 1C is a flow chart illustrating one embodiment of a method for distinguishing genuinely different clonotypes from clonotypes that differ solely by sequencing errors (which should be coalesced). -
FIG. 1D illustrates the form of a numerical function used in one embodiment for determining whether or not to coalesce related clonotypes. -
FIGS. 1E-1F illustrate the use of sequence trees in the methods of the invention. -
FIGS. 2A-2B show a two-staged PCR scheme for amplifying TCRβ genes. -
FIG. 3A illustrates a PCR product to be sequenced that was amplified using the scheme ofFIGS. 2A-2B .FIG. 3B illustrates details of determining a nucleotide sequence of the PCR product ofFIG. 3A .FIG. 3C illustrates details of another embodiment of determining a nucleotide sequence of the PCR product ofFIG. 3A . -
FIG. 4A illustrates a PCR scheme for generating three sequencing templates from an IgH chain in a single reaction.FIGS. 4B-4C illustrates a PCR scheme for generating three sequencing templates from an IgH chain in three separate reactions after which the resulting amplicons are combined for a secondary PCR to add P5 and P7 primer binding sites.FIG. 4D illustrates the locations of sequence reads generated for an IgH chain.FIG. 4E illustrates the use of the codon stricture of V and J regions to improve base calls in the NDN region. - The practice of the present invention may employ, unless otherwise indicated, conventional techniques and descriptions of molecular biology (including recombinant techniques), bioinformatics, cell biology, and biochemistry, which are within the skill of the art. Such conventional techniques include, but are not limited to, sampling and analysis of blood cells, nucleic acid sequencing and analysis, and the like. Specific illustrations of suitable techniques can be had by reference to the example herein below. However, other equivalent conventional procedures can, of course, also be used. Such conventional techniques and descriptions can be found in standard laboratory manuals such as Genome Analysis: A Laboratory Manual Series (Vols. I-IV); PCR Primer: A Laboratory Manual; and Molecular Cloning: A Laboratory Manual (all from Cold Spring Harbor Laboratory Press); Gusfield, Algorithms on Strings, Trees, and Sequences (Cambridge University Press, 1997); and the like.
- In one aspect, the invention is directed to a method for obtaining and analyzing sequence data from a repertoire of immune molecules, such as T cell receptors (TCRs), to rapidly and efficiently determine a clonotype profile. Sequence data typically comprises a large collection of sequence reads, i.e. sequences of base calls and associated quality scores. A key challenge in constructing clonotype profiles is to rapidly and accurately distinguish sequence reads that contain genuine differences from those that contain errors from non-biological sources, such as the sequencing chemistry, amplification chemistry, or the like. An aspect of the invention is the use of sequence trees to represent and organize sequence data from immune molecules extracted from a sample. Such immune molecules typically form an immune repertoire which comprises a very laree set of very similar polynucleotides (e.g. >1000, but more usually from 100,000 to 1,000,000, or more) which arc relatively short in length (e.g. usually less than 300 bp). In one aspect of the invention, the inventors recognized and appreciated that these characteristics permitted the use of sequence tree data structures to efficiently compare sequence reads and candidate clonotypes as compared with analytical approaches using other data structures, e.g. lists, arrays, hash table, or the like.
- The enormous complexity of immune repertoires is well-known, e.g. Arstila et al, Science, 286: 958-961 (1999) and Warren et al (cited above).
FIG. 1A illustrates diagrammatically a typical transcript of an IgH molecule (120) from which a clonotype profile is derived in accordance with some embodiments of the invention. Sources of natural sequence variability include modular recombination of the C, D, J and V segments from large sets carried by the genome, nucleotide additions and deletions to the ends of the D segment to produce the so-called “NDN” regions, and somatic hypermutation where substitutions are made randomly over the length of transcript (122) at a relative frequency roughly as indicated by curve (128). In one aspect of the invention. complex populations of such IgH and TCR transcripts are amplified and sequenced. In one aspect one or both operations for IgH molecules are carried out by using redundant primers annealing to different sites in the V regions (described more fully below). This is particularly advantageous where a sequencing chemistry is employed that has a relatively high error rate or where such sequence variability is difficult or impossible to know beforehand. In the latter case, primer extension for amplification or generation of sequence reads takes place even if one or more primer binding sites are inoperable, or substantially inoperable, because of mismatches caused (for example) by one or more somatic mutations. Starting from promoter P (122) relative frequency shown by curve (128) climbs through leader region (124) to a maximum over the V(D)J region (126) of the transcript after which it drop to near zero. In one aspect of the invention, a segment of recombined B cell nucleic acid is amplified by a PCR with a plurality of forward primers or a plurality of reverse primers to generate a nested set of templates (seeFIGS. 4A and 4B and their descriptions below). Templates from such a set may be further amplified on a surface to form separate amplicons (e.g. by bridge PCR using a cBot instrument; Illumina, San Diego, Calif.). Templates from the same nested set may be associated with one another by sequence reads generated at their common ends. Nested sets of templates allow a sequencing chemistry with relative high error rates to be used to analyze longer sequences than otherwise would be possible, while at the same time maintaining high average quality scores over the entire length of the sequence. The nested sets also ensure that at least one sequence read is obtained from a V region even if it has been subjected to somatic hypermutation. In one embodiment, sequencing chemistries may be used for analyzing highly variable nucleic acids, such as IgH molecules, that have error rates no better than the following: 0.2 percent of sequence reads contain at least one error in positions 1-50; 0.2-1.0 percent of sequence reads contain at least one error in positions 51-75; 0.5-1.5 percent of sequence reads contain at least one error in positions 76-100; and 1-5 percent of sequence reads contain at least one error in positions 101-125. In view of the above, the method of the invention includes steps for distinguishing clonotype sequences that are closely related and genuinely different from those that are closely related and the result of sequencing or other error. - Constructing clonotypes from sequence read data depends in part on the sequencing method used to generate such data, as the different methods have different expected read lengths and data quality. In one approach, a Solexa sequencer is employed to generate sequence read data for analysis. In one embodiment, a sample is obtained that provides at least 0.5-1.0×106 lymphocytes to produce at least 1 million template molecules, which after optional amplification may produce a corresponding one million or more clonal populations of template molecules (or clusters). For most high throughput sequencing approaches, including the Solexa approach, such over sampling at the cluster level is, desirable so that each template sequence is determined with a large degree of redundancy to increase the accuracy of sequence determination. For Solexa-based implementations, preferably the sequence of each independent template is determined 10 times or more. For other sequencing approaches with different expected read lengths and data quality, different levels of redundancy may be used for comparable accuracy of sequence determination. Those of ordinary skill in the art recognize that the above parameters, e.g. sample size, redundancy, and the like, arc design choices related to particular applications.
- Reducing a set of reads for a given sample to a set of distinct clonotypes and recording the number of reads for each clonotype would be a trivial computational problem if sequencing technology was error free. However, in the presence of sequencing errors, each genuine clonotype is surrounded by a ‘cloud’ of reads with, varying numbers of errors with respect to the its sequence. The “cloud” of sequencing errors drops off in density as the distance increases from the clonotype in sequence space. A variety of algorithms are available for converting sequence reads into clonotypes. In one aspect, coalescing of sequence reads (that is, merging candidate clonotypes determined to have one or more sequencing errors) depends on at least three factors: the number of sequences obtained for each of the clonotypes being compared; the number of bases at which they differ: and the sequencing quality score at the positions at which they are discordant. A likelihood ratio may be constructed and assessed that is based on the expected error rates and binomial distribution of errors. For example, two clonotypes, one with 150 reads and the other with 2 reads with one difference between them in an area of poor sequencing quality will likely be coalesced as they are likely to be generated by sequencing error. On the other hand two clonotypes, one with 100 reads and the other with 50 reads with two differences between them are not coalesced as they are considered to be unlikely to be generated by sequencing error. In one embodiment of the invention, the algorithm described below may be used for determining clonotypes from sequence reads. Some of these concepts are illustrated in
FIG. 1B . In one aspect of the invention, sequence reads are first converted into candidate clonotypes. Such a conversion depends on the sequencing platform employed. For platforms that generate high Q score long sequence reads, the sequence read or a portion thereof may. be taken directly as a candidate clonotype. For platforms that generate lower Q score shorter sequence reads, some alignment and assembly steps may be required for converting a set of related sequence reads into a candidate clonotype. For example, for Solexa-based platforms, in some embodiments, candidate clonotypes are generated from collections of paired reads from multiple clusters, e.g. 10 or more, as mentioned above. The frequencies of candidate clonotypes may be plotted in sequence space, as illustrated inFIG. 1B , where such space is reduced to one dimension (the horizontal axis) for sake of illustration. The vertical axis gives the magnitude of each candidate clonotype's frequency, log(read count), or some like measure. In the figure, candidate clonotypes are represented by the various symbols (130). In accordance with one embodiment of the invention, whether two candidate clonotypes are coalesced depends on their respective frequencies or read counts (as noted above), the number of base differences between them (the more differences, the less likely is coalescence), and the quality scores of the bases at the locations where the respective sequences differ (higher quality scores makes coalescence less likely). Candidate clonotypes may be considered in the order of their respective frequencies.FIG. 1B shows candidate clonotype 1 (132), candidate clonotype 7 (134) and candidate clonotype 11 (136) as the three candidates with the highest three frequencies. Related to each such candidate clonotype are other candidate clonotypes that are close in sequence, but with lesser frequencies, such as (i) for candidate clonotype 1 (132) there are candidate clonotype 2 (138) and thecandidate clonotypes 3, 4, 5 and 6 enclosed by cone (140); for candidate clonotype 7 (134) therearc candidate clonotypes 8, 9 and 10 enclosed by cone (142); and (iii) for candidate clonotype 11, there is candidate clonotype 12 enclosed by cone (144). The cones represent likelihood boundaries within which a lesser frequency candidate clonotype would be coalesced with one of the higherfrequency candidate clonotypes 1, 7 or 11. Such likelihood boundaries are functions of the frequency of the nearby candidate clonotypes (3, 4, 5 and 6 for 1; 8, 9 and 10 for 7; and 12 for 11) and their distances in sequence space from the respective higher frequency candidate clonotypes. Candidate clonotype 2 (138) is outside cone (140); thus, it would not be coalesced with candidate clonotype 1 (132). Again, the likelihood (of coalesce) boundaries are shown as cones because candidate clones with higher frequencies are more likely to be genuinely different clonotypes than those of lower frequencies and multiple differences at lower frequencies are more likely to be errors than multiple differences at higher frequencies. - The cloud of sequence reads surrounding each candidate cloncotype can be modeled using the binomial distribution and a simple, model for the probability of a single base error. This latter error model can be inferred from mapping V and J segments or from the clonotype finding algorithm itself, via self-consistency and convergence. A model is constructed for the probability of a given ‘cloud’ sequence Y with read count C2 and E errors (with respect to sequence X) being part of a true clonotype sequence X with perfect read count C1 under the null model that X is the only true clonotype in this region of sequence space. A decision is made whether or not to coalesce sequence Y into the clonotype X according the parameters C1, C2, and E. For any given C1 and E a max value C2 is pre-calculated for deciding to coalesce the sequence Y. The max values for C2 are chosen so that the probability of failing to coalesce Y under the null hypothesis that Y is part of clonotype X is less than some value P after integrating over all possible sequences Y with error E in the neighborhood of sequence X. The value P is controls the behavior of the algorithm and makes the coalescing more or less permissive.
- If a sequence Y is not coalesced into clonotype X because its read count is above the threshold C2 for coalescing into clonotype X then it, becomes a candidate for seeding separate clonotypes (such as with candidate clonotype 2 (138) in
FIG. 1B ). An algorithm implementing such principles would also make sure that any other sequences Y2, Y3, etc. which are ‘nearer’ to this sequence Y (that had been deemed independent of X) are not aggregated into X. This concept of ‘nearness’ includes both error counts with respect to Y and X and the absolute read count of X and Y, i.e. it is modeled in the same fashion as the above model for the cloud of error sequences around clonotype X. In this way ‘cloud’ sequences can be properly attributed to their correct clonotype if they happen to be ‘near’ more than one clonotype. Thus, going toFIG. 1B , ifcandidate clonotype 2 is deemed to be genuinely distinct from candidate clonotype 1 (132), then a special routine, or subalgorithm, would provide a rule for determining which of candidate clonotypes 1 (132) and 2 (138), candidates 4 and 5, between 1 and 2, should be coalesced to (if either). - In one embodiment, an algorithm proceeds in a top down fashion by starting with the sequence X with the highest read count. This sequence seeds the first clonotype. Neighboring sequences are either coalesced into this clonotype if their counts are below the precalculated thresholds (see above), or left alone if they are above the threshold or ‘closer’ to another sequence that was not coalesced. After searching all neighboring sequences within a maximum error count, the process of coalescing reads into clonotype X is finished. Its reads and all reads that have been coalesced into it are accounted for and removed from the list of reads available for making other clonotypes. The next sequence is then moved on to with the highest read count. Neighboring reads are coalesced into this clonotype as above and this process is continued until there are no more sequences with read counts above a given threshold, e.g. until all sequences with more than 1 count have been used as seeds for clonotypes.
- As mentioned above, in another embodiment of the above algorithm, a further test may be added for determining whether to coalesce a candidate sequence Y into an existing clonotype X, which takes into account quality score of the relevant sequence reads. The average quality store(s) are determined for sequence(s) Y (averaged across all reads with sequence Y) were sequences Y and X differ. If the average score is above a predetermined value then it is more likely that the difference indicates a truly different clonotype that should not be coalesced and if the average score is below such predetermined value then it is more likely that sequence Y is caused by sequencing errors and therefore should be coalesced into X.
- Successful implementation of the above algorithm for coalescing candidate clonotypes is dependent upon having an efficient way of finding all sequences with less than E errors (i.e. less than some sequence distance measure) from some input sequence X. This problem is solved using a sequence tree. The implementation of such trees has sonic unusual features in that the nodes of the tree are not restricted to being single letters of the DNA sequences of the candidate clonotypes, as illustrated in
FIG. 1E . The nodes can have arbitrarily long sequences, which allows for a more efficient use of computer memory. - All of the reads of a given sample are placed into the sequence tree. Each leaf nodes holds pointers to its associated reads. A unique sequence of a candidate clonotype is retrieved by traversing backwards in the tree from the leaf to the root node. The first sequence is placed into a simple tree with one root node and one leaf node that contains the full sequence of the read. Sequences arc next added one by one. For each added sequence either a new branch is formed at the last point of common sequence between the read and the existing tree or add the read to an existing leaf node if the tree already contains the sequence. Having placed all the reads into the tree it is easy to use the tree for the following purposes: 1) Finding the highest read count: sorting leaf nodes by read count allows one to find the leaf node (i.e. sequence) with the most reads, and successively lower numbers of reads; 2) Finding neighboring leafs: for any sequence all paths through the tree which have less than X errors with respect to this sequence are searchable. A path is started at the root and branch this path into separate paths proceeding along the tree. The current error count of each path as proceeding along the tree is noted. When the error count exceeds the max allowed errors the given path is terminated. In this way large parts of the tree arc pruned as early as possible. This is an efficient way of finding all paths.(i.e. all leafs) within X errors from any given sequence.
- Features of the above concepts are illustrated in more detail in the flow chart of
FIG. 1C . A set of candidate clonotypes is obtained from sequence data obtained by sequencing recombined nucleic acids extracted from a sample of T cells or B cells. In one aspect, candidate clonotypes each include an NDN region and portions of V and J regions. These sequences are organized into a data structure (150), which may be a sequence tree. Not shown inFIG. 1C , as part of generating a set of candidate clonotypes, in one embodiment, sequence trees may also be constructed for known V regions and known J regions. Sequence reacts making up a candidate clonotype may then be mapped, or aligned, to these known sequences via the sequence trees to, efficiently determine the most likely known V and J sequences of the candidate clonotypes. Returning toFIG. 1C , once the candidate clonotypes are generated, a data structure, such as a sequence tree, is constructed for uses in a method for distinguishing genuine clonotypes from candidate clonotypes that contain experimental or measurement errors, such as sequencing errors. The candidate clonotype that has the highest frequency of occurrence among the current candidate clonotypes (HFCCk) is selected (152) from the data structure, for example a sequence tree; in other words, HFCCk is the candidate clonotype with the highest number of copies, or read counts in cycle k. Next, neighboring lesser frequency candidate clonotypes arc identified (LFCCs) (154); that is, candidate clonotypes within a distance of Dk arc identified. In one aspect of the invention, this identification is carried out using a sequence tree, which allows efficient sequence comparisons of relatively short (<300 bp) sequences. In one embodiment, the comparisons, or sequence alignments, are carried out using dynamic programming, e.g. as disclosed by Gusfield (cited above). In a further embodiment, such dynamic programming is banded dynamic programming where sequences that differ from the selected HFCC by more than a predetermined distance are not considered, which speeds the computation. The candidates HFCCk and LFCCj may be compared on the basis of many different criteria or properties. In one aspect, as mentioned above, candidate clonotypes arc compared on the basis of at least two properties: (i) frequency or read counts and (ii) sequence differences. In another aspect, as mentioned above, candidate clonotypes are compared on the basis of at least three properties: (i) frequency or read counts, (ii) sequence differences, and (iii) quality scores or measures of the bases where differences occur. In one embodiment, sequence differences include base substitutions; in another embodiment, sequence differences include base substitutions, deletions and insertions. The latter embodiment is especially applicable whenever sequence data is generated by sequencing-by-synthesis chemistries that do not employ terminators, such as 454 sequencers and Ion Torrent sequencers. Such sequencing approaches differentiate different sized homopolymer stretches by signal amplitude; thus, base-calling routines in such approaches are prone to insertion and deletion errors, because the difference in signal level from homopolymers differing by one nucleotide drops percipitiousty with increasing homopolymer size (that is, a 2-mer is readily distinguished from a 3-mer, but an 8-mer is almost indistinguishable from a 9-mer). In one aspect, comparisons of HFCCs and LFCCs may be implemented using a function, such as P(HFCCk, LFCCj, D, Q) shown in decision box (158), which depends on the quantities (i) through (iii) described above. Such a function may take many different forms, but generally the value of P changes with changes in (i), (ii) and (iii) as follows: The value of P preferably increases monotonically with the frequency of HFCC and the ratio of the frequency of HFCC to that of LFCC, such that the higher the ratio of the frequency of HFCC to that of LFCC, the higher the likelihood LFCC will be coalesced into HFCC. Likewise, the value of P preferably decreases monotonically with degree to which the sequences of HFCC and LFCC differ, so that the greater the difference between HFCC and LFCC (e.g. as measured by the minimal number of substitutions, insertions or deletions to change one to the other) the lower the likelihood LFCC will be coalesced with HFCC. Finally, the value of P preferably decreases monotonically with increasing quality scores of the locations where the sequences of HFCC and LFCC differ, so that for higher quality scores, the lower the likelihood LFCC will be coalesced with HFCC. When the sequences of HFCC and LFCC differ at more than one location, the quality scores at the different locations may be combined in a variety of differ ways. In one embodiment, whenever there is a plurality of such differences, the plurality of quality scores is expressed as an average value, which may be either an unweighted average or a weighted average.FIG. 1D shows an exemplary function, P, computed for different quality values (curves a through e) for a given sequence difference. As illustrated inFIG. 1D , whenever HFCC is at a level of about 200 read counts (170), then if the quality scores arc determined by curve (a), any LFCC with less than about 50 read counts (172) are coalesced into HFCC. The argument, D, of function P is a measure of the distance between the sequences HFCCk and LFCC; and its value may vary from cycle to cycle as an analysis progresses. (The indices “k” indicates that the values of constants with a “k” subscript may depend on the computational cycle, k.) In one embodiment, D=Dk, so that its value is a function of cycle, number. In another embodiment, D=D(HFCC frequency), so that its value is a function of the frequency of HFCC; independent of cycle number. For example, as the frequency of HFCC decreases, then distance, D, of candidates to be compared decreases. In one embodiment, D is a Hamming distance between HFCCk and LFCCj; however, other distance measures may be used. In one embodiment, Dk is a non-increasing function of k: and in another embodiment, Dk is a decreasing function of k. Decreasing the magnitude of D with increasing cycle number, or with decreasing frequency of HFCC, is advantageous in some embodiments because as a computation progresses to lower and lower frequency candidate clonotypes most such candidates arc singletons, so that sequence distance (rather than frequency difference) becomes the predominant comparison. By lowering D as the computation progresses, unproductive comparisons to distant low frequency candidate clonotypes are reduced, thereby speeding up the computation. Function P may be a complicated expression depending on the number of factors being considered.FIG. 1D illustrates computed values for one embodiment of P which relates read count thresholds for coalescing an LFCC given a read count of an HFCC for different quality scores, as described above. Curves “a” through “c” represent the relationships for different quality scores (with curve “a” corresponding to the highest quality score). - Returning to
FIG. 1C , if P<Pk, then LFCC; is not coalesced with HFCCk and another LFCC is selected (160). If P>Pk, then LFCC; is coalesced with HFCCk (162), in which case another LFCC is selected (166), unless there are no more LFCC left to evaluate (164). If there, are no more LFCC to evaluate (164), then the current HFCCk (including all of the LFCC's coalesced into it) is removed (168) from the data structure, such as the sequence tree. Such removal is illustrated in the simple sequence tree (190) ofFIGS. 1E-F . There, path (192) (indicated by dashed line) in sequence tree (190) corresponds to HFCC (196), which is coalesced with LFCC (198). After coalescence, the segment of path (192) in shaded area (199) is removed from sequence tree (190)to give reduced sequence tree (197) shown inFIG. 1F , which is used in subsequent computations to find neighboring LFCC (154). After such removal, clonotype determination is finished if a stopping criterion (170) is met. In one embodiment, stopping criterion (170) is whether the last non-singleton candidate clonotype has been processed (152). In another embodiment, stopping criterion (170) is whether the frequency or the read counts of the selected HFCC is below that corresponding to a single lymphocyte. In one aspect of the method of the invention, an amplification step may result in each lymphocyte in a sample being represented by multiple copies of the same clonotype; thus, in one embodiment, whatever HFCC has a number of read counts below the number corresponding to a single lymphocyte, then the computation is stopped. In some embodiments, such a number of read counts (or candidate clonotype copies) is 10; in another embodiment, such number is 20; in another embodiment, such a number is 30; in another embodiment, such a number is 40. If the stopping criterion is not met, then the next HFCC is selected (172). The analytical steps summarized in the flow chart ofFIG. 1C may be implemented in any suitable programming language, such as C, C++. Java, C#, Fortran, Pascal or the like. - In accordance with one aspect of the invention; the above method for determining clonotypes and/or clonotype profiles comprises steps of (a) forming a data structure of recombined immune molecules from sequence reads obtained by high throughput nucleic acid sequencing, (b) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype, (c) removing the coalesced candidate clonotype from the data structure, and (d) repeatine, steps (b) and (c) until a clonotype profile is formed. In one embodiment, the data structure is a sequence tree.
- In accordance with another aspect of the invention, the above method of determining clonotypes may be carried out by steps comprising: (a) providing a set of sequence reads from a repertoire of recombined immune molecules each having a V region, an NDN region and a J region wherein for each such molecule at least one sequence read encompasses at least a portion of the NDN region of such molecule; (b) forming from sequence reads encompassing at least a portion of an NDN region a sequence tree having leaves representing candidate clonotypes, each leaf and its corresponding candidate clonotype having a frequency; (c) coalescing with a highest frequency candidate clonotype any lesser frequency candidate clonotypes whenever such lesser frequency is below a predetermined frequency value and a sequence difference therebetween is below a predetermined difference value to form a clonotype having a sequence of the highest frequency candidate clonotype; (d) removing leaves corresponding to the coalesced candidate clonotypes from the sequence tree; and (e) repeating steps (c) and (d) until a highest frequency of a lesser frequency candidate clonotype is below a predetermined stopping value. In one embodiment, the step of forming further includes selecting a highest frequency candidate clonotype and identifying all said lesser frequency candidate clonotypes having a sequence difference therewith less than a predetermined difference value to form a coalescence subset. Thus, in such embodiment, one may limit the total number of LFCCs that must be compared for the coalescing operation (only ones within the predetermined difference value are considered). Such value is a process input depending on the application, e.g. the size of the repertoire, how much computing time is used, and so on. As mentioned above, the function used for deciding whether to coalesce an HFCC with a LFCC can have a variety of forms. In one general aspect, for the step of coalescing, such a function may have the following properties. It depends on frequencies of HFCC, LFCC, the sequence difference therebetween (which may be expressed as a conventional string difference measure, such as a Hamming distance) and quality scores of the one or more nucleotide locations where the HFCC and LFCC differ; such that the function (i) monotonically increases with increasing ratio of frequency of HFCC and frequency of LFCC, (ii) monotonically decreases with increasing sequence difference between HFCC and LFCC, and (iii) monotonically decreases with increasing quality scores of the one or more nucleotide locations. That is, in regard to property (iii), the surer one is that HFCC and LFCC are different (e.g., because there is a high level of confidence in the base calls), then the less likely they will be coalescenced.
- In the foregoing, selection of a predetermined frequency value and a predetermined difference value is a design choice that depend on particular applications. Factors affecting such choices may include details of the biology, speed of implementation, and the like.
- Complex populations of nucleic acids for analysis may arise from a variety of sources. Immune system repertoires may be obtained from samples of immune cells. For example, immune cells can include T-cells and/or B-cells. T-cells (T lymphocytes) include, for example, cells that express T cell receptors. T-cells include helper T cells (effector T cells or Th cells), cytotoxic T cells (CTLs), memory T cells, and regulatory T cells. In one aspect a sample of T cells includes at least 1,000T cells; but more typically, a sample includes at least 10,000 T cells, and more typically, at least 100;000 T cells. In another aspect, a sample includes a number of T cells. In the range of from 1000 to 1,000,000 cells. A sample of immune cells may also comprise B cells. B-cells include, for example, plasma B cells, memory B cells, B1 cells, B2 cells, marginal-zone B cells, and follicular B cells. B-cells can express immunoglobulins (antibodies, B cell receptor). As above, in one aspect a sample of B cells includes at least 1,000 B cells; but more typically, a sample includes at least 10,009 B cells, and more typically, at least 100,000 B cells. In another aspect, a sample includes a number of Bcells in the range of from 1000 to 1,000,000 B cells.
- The sample can include nucleic acid, for example, DNA (e.g., genomic DNA or thitocliondrial DNA) or RNA (e.g., messenger RNA or microRNA). The nucleic acid can be cell-free DNA or RNA, e.g. extracted from the circulatory system, Vlassov et al, Curr. Mol. Med., 10: 142-165 (2010); Swamp et al, FEBS Lett., 581: 795-799 (2007). In the methods of the provided invention, the amount of RNA or DNA from a subject that can be analyzed includes, for example, as low as a single cell in some applications (e.g., a calibration test) and as many as 10 million of cells or more translating to a range of DNA of 6 pg-60 ug, and RNA of approximately 1 pg-10 ug.
- As discussed more fully below (Definitions), a sample of lymphocytes is sufficiently large so that substantially every T cell or B cell with a distinct clonotype is represented therein, thereby forming a repertoire (as the term is used herein). In one embodiment, a sample is taken that contains with a probability of ninety-nine percent every clonotype of a population present at a frequency of 0.001 percent or greater. In another embodiment, a sample is taken that contains with a probability of ninety-nine percent every clonotype of a population present at a frequency of 0.0001 percent or greater. In one embodiment, a sample of B cells or T cells includes at least a half million cells, and in another embodiment such sample includes at least one million cells.
- Whenever a source of material from which a sample is taken is scarce, such as, clinical study samples, or the like, DNA from the material may be amplified by a non-biasing technique, such as whole genome amplification (WGA), multiple displacement amplification (MDA); or like technique, e.g. Hawkins et al, Curr. Opin. Biotech., 13: 65-67 (2002); Dean et al, Genome Research, 11: 1095-1099 (2001); Wang et al, Nucleic. Acids Research, 32: e76 (2004); Hosono ct al, Genome Research, 13: 954-964 (2003); and the like.
- 100401 Blood samples are of particular interest, especially in monitoring lymphoid neoplasms, such as lymphomas, leukemias, or the like; and may be obtained using conventional techniques, e.g. Innis et al, editors, PCR Protocols (Academic Press, 1990); or the like. For example, white blood cells may be separated from blood samples using convention techniques, e.g. RosetteSep kit (Stem Cell Technologies, Vancouver, Canada). Blood samples may range in volume from 100 μL to 10 mL; in one aspect, blood sample volumes are in the range of from 200 100 μL to 2 mL. DNA and/or RNA may then be extracted from such blood sample using conventional techniques for use in methods of the invention, e.g. DNeasy Blood & Tissue Kit (Qiagen, Valencia, Calif.). Optionally, subsets of white blood cells, e.g. lymphocytes, may be further isolated using conventional techniques, e.g. fluorescently activated cell sorting (FACS) (Becton Dickinson, San Jose, Calif.), magnetically activated cell sorting (MACS) (Miltenyi Biotec, Auburn, Calif.), or the like.
- Since the identifying recombinations are present in the DNA of each individual's adaptive immunity cell as well as their associated RNA transcripts, either RNA or DNA can be sequenced in the methods of the provided invention. A recombined sequence from a T-cell or B-cell encoding a T cell receptor or immunoglobulin molecule, or a portion thereof, is referred to as a clonotype. The DNA or RNA can correspond to sequences from T-cell receptor (TCR) genes or immunoglobulin (Ig) genes that encode antibodies. For example, the DNA and RNA can correspond to sequences encoding α, β, γ, or δ chains of a TCR. In a majority of T-cells, the TCR is a heterodimer consisting of an α-chain and β-chain. The TCRα chain is generated by VJ recombination, and the β chain receptor is generated by V(D)J recombination. For the TCRβ chain, in humans there arc 48 V segments, 2 D segments, and 13 J segments. Several bases may be deleted and others added (called N and P nucleotides) at each of the two junctions. In a minority of T-cells, the TCRs consist of γ and δ delta chains. The TCR γ chain is generated by VJ recombination, and the TCR δ chain is generated by V(D)J recombination (Kenneth Murphy, Paul Travers, and Mark Walport, Janeway's Immunology 7th edition, Garland Science, 2007, which is herein incorporated by reference in its entirety).
- The DNA and RNA analyzed in the methods of the invention can correspond to sequences encoding heavy chain immunoglobulins ago with constant regions (α, δ, ε, γ, or μ) or light chain immunoglobulins (IgK or IgL) with constant regions λ or κ. Each antibody has two identical light chains and two identical heavy chains. Each chain is composed of a constant (C) and a variable region. For the heavy chain, the variable region is composed of a variable (V), diversity (D), and joining (J) segments. Several distinct sequences coding for each type of these segments are present in the genome. A specific VDJ recombination event occurs during the development of a B-cell, marking that cell to generate a specific heavy chain. Diversity in the light chain is generated in a similar fashion except that there is no D region so there is only VJ recombination. Somatic mutation often occurs close to the site of the recombination, causing the addition or deletion of several nucleotides, further increasing the diversity of heavy and light chains generated by B-cells. The possible diversity of the antibodies generated by a B-cell is then the product of the different heavy and light chains. The variable regions of the heavy and light chains contribute to form the antigen recognition (or binding) region or site. Added to this diversity is a process of somatic hypermutation which can occur after a specific response is mounted against some epitope.
- As mentioned above, in accordance with the invention, primers may be selected to generate amplicons of subsets of recombined nucleic acids extracted from lymphocytes. Such subsets may be referred to herein as “somatically rearranged regions.” Somatically rearranged regions may comprise nucleic acids from developing or from fully developed lymphocytes, where developing lymphocytes are cells in which rearrangement of immune genes has not been completed to form molecules having full V(D)J regions. Exemplary incomplete somatically rearranged regions include incomplete IgH molecules (such as, molecules containing only D-J regions), incomplete TCRδ molecules (such as, molecules containing only D-J regions), and inactive IgK (for example, comprising Kdc-V regions).
- Adequate sampling of the cells is an important aspect of interpreting the repertoire data, as described further below in the definitions of “clonotype” and “repertoire” For example, starting with 1,000 cells creates a minimum frequency that the assay is sensitive to regardless of how many sequencing reads are obtained. Therefore one aspect of this invention is the development of methods to quantitate the number of input immune receptor molecules. This has been implemented this for TCRβ and IgH sequences. In either case the same set of primers are used that are capable of amplifying all the different sequences. In order to obtain an absolute number of copies, a real time PCR with the multiplex of primers is performed along with a standard with a known number of immune receptor copies. This real time PCR measurement can be made from the amplification reaction that will subsequently be sequenced or can be done on a separate aliquot of the same sample. In the case of DNA, the absolute number of rearranged immune receptor molecules can be readily converted to number of cells (within 2 fold as some cells will have 2 rearranged copies of the specific immune receptor assessed and others will have one). In the case of cDNA the measured total number of rearranged molecules in the real time sample can be extrapolated to define the total number of these molecules used in another amplification reaction of the same sample. In addition, this method can be combined with a method to determine the total amount of RNA to define the number of rearranged immune receptor molecules in a unit amount (say 1 μg) of RNA assuming a specific efficiency of cDNA synthesis. If the total amount of cDNA is measured then the efficiency of cDNA synthesis need not be considered. If the number of cells is also known then the rearranged immune receptor copies per cell can be computed. If the number of cells is not known, one can estimate it from the total RNA as cells of specific type usually generate comparable amount of RNA. Therefore from the copies of rearranged immune receptor molecules per 1μg one can estimate the number of these molecules per cell.
- One disadvantage of doing a separate real time PCR from the reaction that would be processed for sequencing is that there might be inhibitory effects that are different in the real time PCR from the other reaction as different enzymes, input DNA, and other conditions may be utilized. Processing the products of the real time PCR for sequencing would ameliorate this problem. However low copy number using real time PCR can be die to either low number of copies or to inhibitory effects, or other suboptimal conditions in the reaction.
- Another approach that can be utilized is to add it known amount of unique immune receptor rearranged molecules with a known sequence, i.e. known amounts of one or more internal standards, to the cDNA or gnomic DNA from a sample of unknown quantity. By counting the relative number of molecules that are obtained for the known added sequence compared to the rest of the sequences of the same sample, one can estimate the number of rearranged immune receptor molecules in the initial cDNA sample. (Such techniques for molecular counting arc well-known, e.g. Brenner et al, U.S. Pat. No. 7,537,897, which is incorporated herein by reference). Data from sequencing the added unique sequence can be used to distinguish the different possibilities if a real time PCR calibration is being used as well. Low copy number of rearranged immune receptor in the DNA (or cDNA) would create a high ratio between the number of molecules for the spiked sequence compared to the rest of the sample sequences. On the other hand, if the measured low copy number by real time PCR is due to inefficiency in the reaction, the ratio would not be high.
- As noted below, amplicons of target populations of nucleic acids may be generated by a variety of amplification techniques. In one aspect of the invention, multiplex PCR is used to amplify members of a mixture of nucleic acids, particularly mixtures comprising recombined immune molecules such as T cell receptors, B cell receptors, or portions thereof. Guidance for carrying out multiplex PCRs of such immune molecules is found in the following references, which are incorporated by reference: Morley, U.S. Pat. No. 5,296,351; Gorski, U.S. Pat. No. 5,837,447; Dau, U.S. Pat. No. 6,087,096; Von Dongen et al, U.S. patent publication 2006/0234234; European patent publication EP 1544308B1; and the like. The foregoing references describe the technique referred to as “spectratyping,” where a population of immune molecules are amplified by multiplex PCR after which the sequences of the resulting amplicon are physically separated, e.g. by electrophoresis, in order to determine whether there is a predominant size class. Such a class would indicate a predominant clonal population of lymphocytes-which, in turn, would be indicative of disease state. In spectratyping, it is important to select primers that display little or no cross-reactivity (i.e. that do not anneal to binding sites of other primers); otherwise there may be a false representation of size classes in the amplicon. In the present invention, so long as the nucleic acids of a population are uniformly amplified, cross-reactivity of primers is permissible because the sequences of the amplified nucleic acids are analyzed in the present invention, not merely their sizes. As described more fully below, in one aspect, the step of spatially isolating individual nucleic acid molecules is achieved by carrying out a primary multiplex amplification of a preselected somatically rearranged region or portion thereof (i.e. target sequences) using forward and reverse primers that each have tails non-complementary to the target sequences to produce a first amplicon whose member sequences have common sequences at each end that allow further manipulation. For example, such common ends may include primer binding sites for continued amplification using just a single forward primer and a single reverse primer instead of multiples of each, or for bridge amplification of individual molecules on a solid surface, or the like. Such common ends may be added in a single amplification as described above, or they may be added in a two-step procedure to avoid difficulties associated with manufacturing and exercising quality control over mixtures of long primers (e.g. 50-70 bases or more). In such a two-step process (described more fully below and illustrated in
FIGS. 4A-4B ), the primary amplification is carried out as described above, except that the primer tails are limited in length to provide only forward and reverse primer binding sites at the ends of the sequences of the first amplicon. A secondary amplification is then carried out using secondary amplification primers specific to these primer binding sites to add further sequences to the ends of a second amplicon. The secondary amplification primers have tails non-complementary to the target sequences, which form the ends of the second amplicon and which may, be used in connection With sequencing the clonotypes of the second amplicon. In one embodiment, such added sequences may include primer binding sites for generating sequence reads and primer binding sites for carrying out bridge PCR on a solid surface to generate clonal populations of spatially isolated individual molecules, for example, when Solexa-based sequencing is used. In this latter approach, a sample of sequences from the second amplicon are disposed on a solid surface that has attached complementary oligonucleotides capable of annealing to sequences of the sample, idler which cycles of primer extension, denaturation, annealing are implemented until clonal populations of templates are formed. Preferably, the size of the sample is selected so that (i) it includes an effective representation of clonotypes in the original sample, and (ii) the density of clonal populations on the solid surface is in a range that permits unambiguous sequence determination of clonotypes. - TCR or BCR sequences or portions thereof can be amplified from nucleic acid in a multiplex reaction using at least one primer that anneals to the C region and one or more primers that can anneal to one or more V segments (as illustrated in
FIGS. 2A-2B andFIGS. 4A-4B and discussed more fully below). The region to be amplified can include the full clonal sequence or a subset of the clonal sequence, including the V-D junction, D-J junction of an inummoglobulin or T-cell receptor gene, the full variable region of an immunoglobulin or T-cell receptor gene, the antigen recognition region, or a CDR, e.g., complementarity determining region 3 (CDR3). - The TCR or immunoglobulin sequence can amplified using a primary and a secondary amplification step. Each of the different amplification steps can comprise different primers. The different primers can introduce sequence not originally present in the immune gene sequence. For example, the amplification procedure can add new primer binding sites to the ends of the target sequences to convert a multiplex amplification to a singleplex amplification or the amplification procedure can add one or more tags to the 5′ and/or 3′ end of amplified TCR or immunoglobulin sequence (as illustrated in
FIGS. 3A-3B ). The tag can be sequence that facilitates subsequent sequencing of the amplified DNA. The tag can be sequence that facilitates binding the amplified sequence to a solid support. - After amplification of DNA from the genome (or amplification of nucleic acid in the form of cDNA by reverse transcribing RNA), the individual nucleic acid molecules can be isolated, optionally re-amplified, and then sequenced individually. Exemplary amplification protocols may be found in van Dongen et al, Leukemia, 17: 2257-2317 (2003) or van Dongen et al, U.S. patent publication 2006/0234234, which is incorporated by reference. Briefly, an exemplary protocol is as follows: Reaction buffer: ABI Buffer II or ABI Gold Buffer (Life Technologies, San Diego, Calif.); 50 μL final reaction volume; 100 ng sample DNA; 10 pmol of each primer (subject to adjustments to balance amplification as described below); dNTPs at 200 μM final concentration; MgCl2 at 1.5 mM final concentration (subject to optimization depending on target sequences and polymerase); Taq polymerase (1-2 U/tube); cycling conditions: preactivation 7 min at 95° C.; annealing at 60° C.; cycling times: 30 s denaturation; 30 s annealing; 30 s extension. Polymerases that can be used for amplification in the methods of the invention are commercially available and include, for example, Taq polymerase, AccuPrime polymerase, or Pfu. The choice of polymerase to use can be based on whether fidelity or efficiency is preferred.
- Methods for isolation of nucleic acids from a pool include subcloning nucleic acid into DNA vectors and transforming bacteria (bacterial cloning), spatial separation of the molecules in two dimensions on a solid substrate (e.g., glass slide), spatial separation of the molecules in three dimensions in a solution within micelles (such as can be achieved using oil emulsions with or without immobilizing the molecules on a solid surface such as beads), or using microreaction chambers in, for example, microfluidic or nano-fluidic chips. Dilution can be used to ensure that on average a single molecule is present in a given volume, spatial region, bead, or reaction chamber. Guidance for such methods of isolating individual nucleic acid molecules is found in the following references: Sambrook, Molecular Cloning: A Laboratory Manual (Cold Spring Harbor Laboratory Press, 2001s); Shendure et at, Science, 309: 1728-1732 (including supplemental matcrial) (2005); U.S. Pat. No. 6,300,070; Bentley ct al, Nature, 456:53-59 (including supplemental material) (2008); U.S. Pat. No. 7,323,305; Matsubara ct al, Biosensors & Bioelectronics, 20: 1482-1490 (2005): U.S. Pat. No. 6,753,147; and the like.
- Real time PCR, picogreen staining, nanofluidic electrophoresis (e.g. LabChip) or UV absorption measurements can be used in an initial step to judge the functional amount of amplifiable material.
- In one aspect, multiplex amplifications are carried out so that relative amounts of sequences in a starting population arc substantially the same as those in the amplified population, or amplicon. That is, multiplex amplifications are carried out with minimal amplification bias among member sequences of a sample population. In one embodiment, such relative amounts are substantially the same if each relative amount in an amplicon is within five fold of its value in the starting sample. In another embodiment, such relative amounts are substantially the same if each relative amount in an amplicon is within two fold of its value in the starting sample. As discussed more fully below, amplification bias in PCR may be detected and corrected using conventional techniques so that a set of PCR primers may be selected for a predetermined repertoire that provide unbiased amplification of any sample.
- In regard to many repertoires based on TCR or BCR sequences, a multiplex amplification optionally uses all the V segments. The reaction is optimized to attempt to get amplification that maintains the relative abundance of the sequences amplified by different V segment primers. Some of the primers are related, and hence many of the primers may “cross talk,” amplifying templates that arc not perfectly matched with it. The conditions arc optimized so that each template can be amplified in a similar fashion irrespective of which primers amplified it. In other words if there are two templates, then after 1,000 fold amplification both templates can be amplified approximately 1,000 fold, and it does not matter that for one of the templates half of the amplified products carried a different primer because of the cross talk. In subsequent analysis of the sequencing data the primer sequence is eliminated from the analysis, and hence it does not matter what primer is used in the amplification as long as the templates are amplified equally.
- In one embodiment, amplification bias may be avoided by carrying out a two-stage amplification (as illustrated in
FIGS. 2A-2B ) wherein a small number of amplification cycles are implemented in a first, or primary, stage using primers having tails non-complementary with the target sequences. The tails include primer binding sites that arc added to the ends of the sequences of the primary amplicon so that such sites arc used in a second stage amplification using only a single forward primer and a single reverse primer, thereby eliminating a primary cause of amplification bias. Preferably, the primary PCR will have a small enough number of cycles (e.g. 5-10) to minimize the differential amplification by the different primers. The secondary amplification is done with one pair of primers and hence the issue of differential amplification is minimal. One percent of the primary PCR is taken directly to the secondary PCR. Thirty-five cycles (equivalent to 28 cycles without the 100 fold dilution step) used between the two amplifications were sufficient to show a robust amplification irrespective of whether the breakdown of cycles were: one cycle primary and 34 secondary or 25 primary and 10 secondary. Even though ideally doing only 1 cycle in the primary PCR may decrease the amplification bias, there arc other considerations. One aspect of this is representation. This plays a role when the stalling input amount is not in excess to the number of reads ultimately obtained. For example, if 1,000,000 reads are obtained and starting with 1,000,000 input Molecules then taking only representation from 100,000 molecules to the secondary amplification would degrade the precision of estimating the relative abundance of the different species in the original sample. The 100 fold dilution between the 2 steps means that the representation is reduced unless the primary PCR amplification generated significantly more than 100 molecules. This indicates that a minimum 8 cycles (256 fold), but more comfortably 10 cycle (1,000 fold), may be used. The alternative to that is to take more than 1% of the primary PCR into the secondary but because of the high concentration of primer used in the primary PCR, a big dilution factor is can be used to ensure these primers do not interfere in the amplification and worsen the amplification bias between sequences. Another alternative is to add a purification or enzymatic step to eliminate the primers from the primary PCR to allow a smaller dilution of it. In this example, the primary PCR was 10 cycles and the second 25 cycles. - Any high-throughput technique for sequencing nucleic acids can be used in the method of the invention. Preferably, such technique has a capability of generating in a cost-effective manner a volume of sequence data from which at least 1000 clonotypes can be determined, and preferably, from which at least 10,000 to 1,000,000 clonotypes can be determined. DNA sequencing techniques include classic dideoxy sequencing reactions (Sanger method) using labeled terminators or primers and gel separation in slab or capillary, sequencing by synthesis using reversibly terminated labeled nucleotides, pyrosequencing, 454 sequencing, allele specific hybridization to a library of labeled oligonucleotide probes, sequencing by synthesis using allele specific hybridization to a library of labeled clones that is followed by ligation, real time monitoring of the incorporation of labeled nucleotides during a polymerization step, polony sequencing, and SOLiD sequencing. Sequencing of the separated molecules has more recently been demonstrated by sequential or single extension reactions using polymerases or ligases as well as by single or sequential differential hybridizations with libraries of probes. These reactions have been performed on many clonal sequences in parallel including demonstrations in current commercial applications of over 100 million sequences in parallel. These sequencing approaches can thus be used to study the repertoire of T-cell receptor (TCR) and/or B-cell receptor (BCR). In one aspect of the invention, high-throughput methods of sequencing are employed that comprise a step of spatially isolating individual molecules on a solid surface where they are sequenced in parallel. Such solid surfaces may include nonporous surfaces (such as in Solexa sequencing, e.g. Bentley et al, Nature, 456: 53-59 (2008) or Complete Genomics sequencing, e.g. Drmanac et al, Science, 327: 78-81 (2010)), arrays of wells, which may include bead- or particle-bound templates (such as with 454, e.g. Margulies et al, Nature, 437: 376-380 (2005) or Ion Torrent sequencing, U.S. patent publication 2010/0137143 or 2010/0304982), micromachined membranes (such as with SMRT sequencing, e.g. Eid et al, Science, 323: 133-138 (2009)), or bead arrays (as with SOLiD sequencing or polony sequencing, e.g. Kim et al, Science, 316: 1481-1414 (2007)). In another aspect, such methods comprise amplifying the isolated molecules either before or after they arc spatially isolated on a solid surface. Prior amplification may comprise emulsion-based amplification, such as emulsion PCR, or rolling circle amplification. Of particular interest is Solexa-based sequencing where individual template molecules are spatially isolated on a solid surface, after which they are amplified in parallel by bridge PCR to form separate clonal populations, or clusters, and then sequenced, as described in Bentley et al (cited above) and in manufacturer's instructions (e.g. TruSeq™ Sample Preparation Kit and Data Sheet, Illumina, Inc., San Diego, Calif., 2010); and further in the following references: U.S. Pat. Nos. 6,090,592; 6,300,070; 7,115,400; and EP0972081B1; which are incorporated by reference. In one embodiment, individual molecules disposed and amplified on a solid surface form clusters in a density of at least 105 clusters per cm2; or in a density of at least 5×105 per cm2; or in a density of at least 10 6 clusters per cm2. In one embodiment, sequencing chemistries are employed having relatively high error rates. In such embodiments, the average quality scores produced by such chemistries are monotonically declining functions of sequence read lengths. In one embodiment, such decline corresponds to 0.5 percent of sequence reads have at least one error in positions 1-75; 1 percent of sequence reads have at least one error in positions 76-100; and 2 percent of sequence reads have at least one error in positions 101-125.
- In one aspect, a sequence-based clonotype profile of an individual is obtained using the following steps: (a) obtaining a nucleic acid sample from T-cells and/or B-cells of the individual; (b) spatially isolating individual molecules derived from such nucleic acid sample, the individual molecules comprising at least one template generated from a nucleic acid in the sample, which template comprises a somatically rearranged region or a portion thereof, each individual molecule being capable of producing at least one sequence read; (c) sequencing said spatially isolated individual molecules; and (d) determining abundances of different sequences of the nucleic acid molecules from the nucleic. acid sample to generate the clonotype profile. In one embodiment, each of the somatically rearranged regions comprise a V region and a J region. In another embodiment, the step of sequencing comprises bidirectionally sequencing each of the spatially isolated individual molecules to produce at least one forward sequence read and at least one reverse sequence read. Further to the latter embodiment, at least one of the forward sequence reads and at least one of the reverse sequence reads have an overlap region such that bases of such overlap region are determined by a reverse complementary relationship between such sequence reads. In still another embodiment, each of the somatically rearranged regions comprise a V region and a J region and the step of sequencing further includes determining a sequence of each of the individual nucleic acid molecules from one or more of its forward sequence reads and at least one reverse sequence read starting from a position in a J region and extending in the direction of its associated V region. In another embodiment, individual molecules comprise nucleic acids selected from the group consisting of complete IgH molecules, incomplete IgH molecules, complete IgK complete, IgK inactive molecules, TCRβ molecules, TCRγ molecules, complete TCRδ molecules, and incomplete TCRδ molecules. In another embodiment, the step of sequencing comprises generating the sequence reads having monotonically decreasing quality scores. Further to the latter embodiment, monotonically decreasing quality scores are such that the sequence reads have error rates no better than the following: 0.2 percent of sequence reads contain at least one error in
base positions 1 to 50, 0.2 to 1.0 percent of sequence reads contain at least one error in positions 51-75, 0.5 to 1.5 percent of sequence reads contain at least one error in positions 76-100. In another embodiment, the above method comprises the following steps: (a) obtaining a nucleic acid sample from T-cells and/or B-cells of the individual; (b) spatially isolating individual molecules derived from such nucleic acid sample, the individual molecules comprising nested sets of templates each generated from a nucleic acid in the sample and each containing a somatically rearranged region or a portion thereof, each nested set being capable of producing a plurality of sequence reads each extending in the same direction and each starting from a different position on the nucleic acid from which the nested set was generated; (c) sequencing said spatially isolated individual molecules; and (d) determining abundances of different sequences of the nucleic acid molecules from the nucleic acid sample to generate the clonotype profile. In one embodiment, the step of sequencing includes producing a plurality of sequence reads for each of the nested sets. In another embodiment, each of the somatically rearranged regions comprise a V region and a J region, and each of the plurality of sequence reads starts from a different position in the V region and extends in the direction of its associated J region. - In one aspect, for each sample from an individual, the sequencing technique used in the methods of the invention generates sequences of least 1000 clonotypes per run; in another aspect, such technique generates sequences of at least 10,000 clonotypes per run; in another aspect, such technique generates sequences of at least 100,000 clonotypes per run; in another aspect, such technique generates sequences of at least 500,000 clonotypes per run; and in another aspect, such technique generates sequences of at least 1,000,000 clonotypes per run. In still another aspect, such technique generates sequences of between 100,000 to 1,000,000 clonotypes per run per individual sample.
- The sequencing technique used in the Methods of the provided invention can generate about 30 bp, about 40 bp, about 50 bp, about 60 bp, about 70 bp, about 80 bp, about 90 bp, about 100 bp, about 110, about 120 by per read, about 150 bp, about 200 bp, about 250 bp, about 300 bp, about 350 bp, about 400 bp, about 450 bp, about 500 bp, about 550 bp, or about 600 by per read.
- In one aspect of the invention, sequences of clonotypes (including but not limited to those derived from IgH, TCRα, TCRβ, TCRγ, TCRδ, and/or IgLK (IgK)) may be determined by combining information from one or more sequence reads, for example, along the V(D)J regions of the selected chains. In another aspect, sequences of clonotypes are determined by combining information from a plurality of sequence reads. Such pluralities of sequence reads may include one or more sequence reads along a sense strand (i.e. “forward” sequence reads) and, one or more sequence reads along, its complementary strand (i.e. “reverse” sequence reads). When multiple sequence reads are generated along the same strand, separate templates are first generated by amplifying sample molecules with primers selected for the different positions of the sequence reads. This concept is illustrated in
FIG. 4A where primers (404, 406 and 408) are employed to generate amplicons (410, 412, and 414, respectively) in a single reaction. Such amplifications may be carried out in the same reaction or in separate reactions. In one aspect, whenever PCR is employed, separate amplification reactions arc used for generating the separate templates which, in turn, are combined and used to generate multiple sequence reads along the same strand. This latter approach is preferable for avoiding the need to balance primer concentrations (and/or other reaction parameters) to ensure equal amplification of the multiple templates (sometimes referred to herein as “balanced amplification” or “unbias amplification”). The generation of templates in separate reactions is illustrated inFIGS. 4B-4C . There a sample containing IgH (400) is divided into three portions (472, 474, and 476) which arc added to separate PCRs using J region primers (401) and V region primers (404, 406, and 408, respectively) to produce amplicons (420, 422 and 424, respectively). The latter amplicons are then combined (478) in secondary PCR (480) using P5 and P7 primers to prepare the templates (482) for bridge PCR and sequencing on an Illumina GA sequencer, or like instrument. - Sequence reads of the invention may have a wide variety of lengths, depending in part on the sequencing technique being employed. For example, for some techniques, several trade-offs may arise in its implementation, for example, (i) the number and lengths of sequence reads per template and (ii) the cost and duration of a sequencing operation. In one embodiment, sequence reads are in the range of from 20 to 400 nucleotides: in another embodiment, sequence reads are in a range of from 30 to 200 nucleotides; in still another embodiment, sequence reads arc in the range of from 30 to 120 nucleotides. In one embodiment, 1 to 4 sequence reads are generated for determining the sequence of each clonotype; in another embodiment, 2 to 4 sequence reads arc generated for determining the sequence of each clonotype; and in another embodiment, 2 to 3 sequence reads are generated for determining the sequence of each clonotype. In the foregoing embodiments, the numbers given are exclusive of sequence reads used to identify samples from different individuals. The lengths of the various sequence reads used in the embodiments described below may also vary based on the information that is sought to be captured by the read; for example, the starting location and length of a sequence read may be designed to provide the length of an NDN region as well as its nucleotide sequence: thus, sequence reads spanning the entire NDN region arc selected. In other aspects, one or more sequence reads that in combination (but not separately) encompass a D and/or NDN region arc sufficient.
- In another aspect of the invention, sequences of clonotypes are determined in part by aligning sequence reads to one or more V region reference sequences and one or more J region reference sequences, and in part by base determination without alignment to reference sequences, such as in the highly variable NDN region. A variety of alignment algorithms may be applied to the sequence reads and reference sequences. For example, guidance for selecting alignment methods is available in Batzogloti, Briefings in Bioinformatics, 6: 6-22 (2005), which is incorporated by reference. In one aspect, whenever V reads or C reads (as mentioned above) are aligned to V and J region reference sequences, a tree search algorithm is employed, e.g. as described generally in Gustield (cited above) and Cormen et al, Introduction to Algorithms, Third Edition (The MIT Press, 2009).
- In another aspect, an end of at least one forward read and an end of at least one reverse read overlap in an overlap region (e.g. 308 in
FIG. 3B ), so that the bases of the reads arc in a reverse complementary relationship with one another. Thus, for example, if a forward read in the overlap region is “5′-acgttge”, then a reverse read in a reverse complementary relationship is “5′-gcaacgt” within the same overlap region. In one aspect, bases within such an overlap region arc determined, at least in part, from such a reverse complementary relationship. That is, a likelihood of a base call (or a related quality score) in a prospective overlap region is increased if it preserves,'or is consistent with, a reverse complementary relationship between the two sequence reads. In one aspect, clonotypes of TCR β and IgH chains (illustrated inFIG. 3B ) are determined by at least one sequence read starting in its J region and extending in the direction of its associated V region (referred to herein as “C read” (304)) and at least one sequence read starting in its V region and extending in the direction of its associated J region (referred to herein as a “V read” (306)). Overlap region (308) may or may not encompass the NDN region (315) as shown inFIG. 3B . Overlap region (308) may be entirely in the J region, entirely in the NDN region, entirely in the V region, or it may encompass a J region-NDN region boundary or a V region-NDN region boundary, or both such boundaries (as illustrated inFIG. 3B ). Typically, such sequence reads arc generated by extending sequencing primers, e.g. (302) and (310) inFIG. 3B , with a polymerase in a sequencing-by-synthesis reaction, e.g. Metzger, Nature Reviews Genetics, 11: 31-46 (2010); Fuller et al, Nature Biotechnology, 27: 1013-1023 (2009). The binding sites for primers (302) and (310) are predetermined, so that they can provide a starting point or anchoring point for initial alignment and analysis, of the sequence reads. In one embodiment, a C read is positioned so that it encompasses the D and/or NDN region of the TCR or IgH chain and includes a portion of the adjacent V region, e.g. as illustrated inFIGS. 3B and 3C . In one aspect, the overlap of the V read and the C read in the V region is used to align the reads with one another. In other embodiments, such alignment of sequence reads is not necessary, e.g. with TCRβ chains, so that a V read may only be long enough to identify the particular V region of a clonotype. This latter aspect is illustrated inFIG. 3C . Sequence read (330) is used to identify a V region, with or without overlapping another sequence read, and another sequence read (332) traverses the NDN region and is used to determine the sequence thereof. Portion (334) of sequence read (332) that extends into the V region is used to associate the sequence information of sequence read (332) with that of sequence read (330) to determine a clonotype. For some sequencing methods, such as base-by-base approaches like the Solexa sequencing method, sequencing run time and reagent costs are reduced by minimizing the number of sequencing cycles in an analysis. Optionally, as illustrated inFIG. 3B , amplicon (300) is produced with sample tag (312) to distinguish between clonotypes originating from different biological samples, e.g. different patients. Sample tag (312) may be identified by annealing a primer to primer binding region (316) and extending it (314) to produce a sequence read across tag (312), from which sample tag (312) is decoded. - The IgH chain is more challenging to analyze than TCRβ chain because of at least two factors: i) the presence of somatic mutations makes the mapping or alignment more difficult, and ii) the NDN region is larger so that it is often not possible to map a portion of the V segment to the C read. In one aspect of the invention, this problem is overcome by using a plurality of primer sets for generating V reads, which are located at different locations along the V region, preferably so that the primer binding sites are nonoverlapping and spaced apart, and with at least one primer binding site adjacent to the NDN region, e.g. in one embodiment from 5 to 50 bases from the V-NDN junction, or in another embodiment from 10 to 50 bases from the V-NDN junction. The redundancy of a plurality of primer sets minimizes the risk of failing to detect a clonotype due to a failure of one or two primers having binding sites affected by somatic mutations. In addition, the presence of at least one primer binding site adjacent to the NDN region makes it more likely that a V read will overlap with the C read and hence effectively extend the length of the C read. This allows for the generation of a continuous sequence that spans all sizes of NDN regions and that can also map substantially the entire V and J regions on both sides of the NDN region. Embodiments for carrying out such a scheme are illustrated in
FIGS. 4A and 4D . InFIG. 4A , a sample comprising IgH chains (400) are sequenced by generating a plurality amplicons for each chain by amplifying the chains with a single set of J region primers (401) and a plurality (three shown) of sets of V region (402) primers (404, 406, 408) to produce a plurality of nested amplicons (e.g., 410, 412, 416) all comprising the same NDN region and having different lengths encompassing successively larger portions (411, 413, 415) of V region (402). Members of a nested set may be grouped together after sequencing by noting the identify (or substantial identity) of their respective NDN, J and/or C regions, thereby allowing reconstruction of a longer V(D)J segment than would be the case otherwise for a sequencing platform with limited read length and/or sequence quality. In one embodiment, the plurality of primer sets may be a number in the range of from 2 to 5. In another embodiment the plurality is 2-3; and still another embodiment the plurality is 3. The concentrations and positions of the primers in a plurality may vary widely. Concentrations of the V region primers may or may not be the same. In one embodiment, the primer closest to the NDN region has a higher concentration than the other primers of the plurality, e.g. to insure that amplicons containing the NDN region are represented in the resulting amplicon. In a particular embodiment where a plurality a three primers is employed, a concentration ratio of 60:20:20 is used. One or more primers (e.g. 435 and 437 inFIG. 4B ) adjacent to the NDN region (444) may be used to generate one or more sequence reads (e.g. 434 and 436) that overlap the sequence read (442) generated by J region primer (432), thereby improving the quality of base calls in overlap region (440). Sequence reads from the plurality of primers may or may not overlap the adjacent downstream primer binding site and/or adjacent downstream sequence read. In one embodiment, sequence reads proximal to the NDN region (e.g. 436 and 438) may be used to identify the particular V region associated with the clonotype. Such a plurality of primers reduces the likelihood of incomplete or failed amplification in case one of the primer binding sites is hypermutated during immunoglobulin development. It also increases the likelihood that diversity introduced by hypermutation of the V region will be capture in a clonotype sequence. A secondary PCR may be performed to prepare the nested amplicons for sequencing, e.g. by amplifying with the P5 (401) and P7 (404, 406, 408) primers as illustrated to produce amplicons (420, 422, and 424), which may be distributed as single molecules on a solid surface, where they are further amplified by bridge PCR, or like technique. - Base calling in NDN regions (particularly of IgH chains) can be improved by using the codon stricture of the flanking J and V regions, as illustrated in
FIG. 4E . (As used herein, “codon.structure” means the codons of the natural reading frame of segments of TCR or BCR transcripts or genes outside of the NDN regions, e.g. the V region, J region, or the like.) There amplicon (450), which is an enlarged view of the amplicon ofFIG. 4B , is shown along with the relative positions of C read (442) and adjacent V read (434) above and the codon structures (452 and 454) of V region (430) and J region (446), respectively, below. In accordance with this aspect of the invention, after the codon structures (452 and 454) are identified by conventional alignment to the V and J reference sequences, bases in NDN region (456) are called (or identified) one base at a time moving from J region (446) toward V region (430) and in the opposite direction from V region (430) toward J region (446) using sequence reads (434) and (442). Under normal biological conditions, only the recombined TCR or IgH sequences that have in frame codons from the V region through the NDN region and to the J region are expressed as proteins. That is, of the variants generated somatically only ones expressed are those whose J region and V region codon frames are in-frame with one another and remain in-frame through the NDN region. (Here the correct frames of the V and J regions arc determined from reference sequences). If an out-of-frame sequence is identified based one or more low quality base calls, the corresponding clonotype is flagged for re-evaluation or as a potential disease-related anomaly. If the sequence identified is in-frame and based on high quality base calls, then there is greater confidence that the corresponding clonotype has been correctly called. Accordingly, in one aspect, the invention includes a method of determining V(D)J-based clonotypes from bidirectional sequence reads comprising the steps of: (a) generating at least one J region sequence read that begins in a J region and extends into an NDN region and at least one V region sequence read that begins in the V regions and extends toward the NDN region such that the J region sequence read and the V region sequence read are overlapping in an overlap region, and the J region and the V region each have a codon structure; (b) determining whether the codon structure of the J region extended into the NDN region is in frame with the codon structure of the V region extended toward the NDN region. In a further embodiment, the step of generating includes generating at least one V region sequence read that begins in the V region and extends through the NDN region to the J region, such that the J region sequence read and the V region sequence read arc overlapping in an overlap region. - Somatic Hypermutations. In one embodiment, IgH-based clonotypes that have undergone somatic hypermutation are determined as follows: A somatic mutation is defined as a sequenced base that is different from the corresponding base of a reference sequence (of the relevant segment, usually V, J or C) and that is present in a statistically significant number of reads. In one embodiment, C reads may be used to find somatic mutations with respect to the mapped J segment and likewise V reads for the V segment. Only pieces of the C and V reads are used that are either directly mapped to J or V segments or that are inside the clonotype extension up to the NDN boundary. In this way, the NDN region is avoided and the same ‘sequence information’ is not used for mutation finding that was previously used for clonotype determination (to avoid erroneously classifying as mutations nucleotides that arc really just different recombined NDN regions). For each segment type, the mapped segment (major allele) is used as a scaffold and all reads are considered which have mapped to this allele during the read mapping phase. Each position of the reference sequences where at least one read has mapped is analyzed for somatic mutations. In one embodiment, the criteria for accepting a non-reference base as a valid mutation include the following: 1) at least N reads with the given mutation base, 2) at least a given fraction N/M reads (where M is the total number of mapped reads at this base position) and 3) a statistical cut based on the binomial distribution, the average Q score of the N reads at the mutation base as well as the number (M-N) of reads with a non-mutation base. Preferably, the above parameters arc selected so that the false discovery rate of mutations per clonotype is less than 1 in 1000, and more preferably, less than 1 in 10000.
- Phylogenic Clonotypes (Clans). In cancers, such as lymphoid neoplasms, a single lymphocyte progenitor may give rise to many related lymphocyte progeny, each possessing and/or expressing a slightly different TCR or BCR, and therefore a different clonotype, due to cancer-related somatic mutation(s), such as base substitutions, aberrant rearrangements, or the like. Cells producing such clonotypes are referred to herein as phylogenic clones, and a set of such related clones are referred to herein as a “clan.” Likewise, clonotypes of phylogenic clones are referred to as phylogenic clonotypes and a set of phylogenic clonotypes may be referred to as a clan of clonotypes. In one aspect, methods of the invention comprise monitoring the frequency of a clan of clonotypes (i.e., the sum of frequencies of the constituent phylogenic clonotypes of the clan), rather than a frequency of an individual clonotype. Phylogenic clonotypes may be identified by one or more measures of relatedness to a parent clonotype. In one embodiment, phylogenic clonotypes may be grouped into the same clan by percent homology, as described more fully below. In another embodiment, phylogenic clonotypes are identified by common usage of V regions, J regions, and/or NDN regions. For example, a clan may be defined by clonotypes having common J and ND regions but different V regions; or it may be defined by clonotypes having the same V and J regions (including identical base substitutions mutations) but with different NDN regions; or it may be defined by a clonotype-that has undergone one or more insertions and/or deletions of from 1-10 bases, or from 1-5 bases, or from 1-3 bases, to generate clan members. In another embodiment, members of a clan are determined as follows. Clonotypes are assigned to the same clan if they satisfy the following criteria: i) they are mapped to the same V and J reference segments, with the mappings occurring at the same relative positions in the clonotype sequence, and ii) their NDN regions arc substantially identical. “Substantial” in reference to clan membership means that some small differences in the NDN region are allowed because somatic mutations may have occurred in this region. Preferably, in one embodiment, to avoid falsely calling a mutation in the NDN region, whether a base substitution is accepted as a cancer-related mutation depends directly on the size of the NDN region of the clan. For example, a method may accept a clonotype as a clan member if it has a one-base difference from clan NDN sequence(s) as a cancer-related mutation if the length of the clan NDN sequence(s) is m nucleotides or greater, e.g. 9 nucleotides or greater, otherwise it is not accepted, or if it has a two-base difference from clan NDN sequence(s) as cancer-related mutations if the length of the clan NDN sequence(s) is n nucleotides or greater, e.g. 20 nucleotides or greater, otherwise it is not accepted, In another embodiment, members of a clan are determined using the following criteria: (a) V read maps to the same V region, (b) C read maps to the same J region, (c) NDN region substantially identical (as described above), and (d) position of NDN region between V-NDN boundary and J-NDN boundary is the same (or equivalently, the number of downstream base additions to D and the number of upstream base additions to D are the same). Clonotypes of a single sample may be grouped into clans and, clans from successive samples acquired at different times may be compared with one another. In particular; in one aspect of the invention, clans containing clonotypes correlated with a disease, such as a lymphoid neoplasm, are identified from clonotypes of each sample and compared with that of the immediately previous sample to determine disease status, such as, continued remission, incipient relapse, evidence of further clonal evolution, or the like.
- It is expected that PCR error is concentrated in some bases that were mutated in the early cycles of PCR. Sequencing error is expected to be distributed in many bases even though it is totally random as the error is likely to have some systematic biases. It is assumed that some bases will have sequencing error at a higher rate, say 5% (5 fold the average). Given these assumptions, sequencing error becomes the dominant type of error. Distinguishing PCR errors from the occurrence of highly related clonotypes will play a role in analysis. Given the biological significance to determining that there are two or more highly related clonotypes, a conservative approach to making such calls is taken. The detection of enough of the minor clonotypes so as to be sure with high confidence (say 99.9%) that there are more than one clonotype is considered. For example of clonotypes that are present at 100 copies/1,000,000, the minor variant is detected 14 or more times for it to be designated as an independent clonotype. Similarly, for clonotypes present at 1,000 copies/1,000,000 the minor variant can be detected 74 or more times to be designated as an independent clonotype. This algorithm can be enhanced by using the base quality score that is obtained with each sequenced base. If the relationship between quality score and error rate is validated above, then instead of employing the conservative 5% error rate for all bases, the quality score can be used to decide the number of reads that need to be present to call an independent clonotype. The median quality score of the specific base in all the reads can be used, or more rigorously, the likelihood of being an error can be computed given the quality score of the specific base in each read, and then the probabilities can be combined (assuming independence) to estimate the likely number of sequencing error for that base. As a result, there are different thresholds of rejecting the sequencing error hypothesis for different bases with different quality scores. For example for a clonotype present at 1,000 copies/1,000,000 the minor variant is designated independent when it is detected 22 and 74 times if the probability of error were 0.01 and 0.05, respectively.
- In this example, TCRβ chains are analyzed. The analysis includes amplification, sequencing, and analyzing the TCRβ sequences. One primer is complementary to a common sequence in Cβ1 and Cβ2, and there arc 34 V primers capable of amplifying all 48 V segments. Cβ1 or Cβ2 differ from each other at
position 10 and 14 from the J/C junction. The primer for Cβ1 and Cβ2 ends at position 16 by and has no preference for Cβ1 or Cβ2. The 34 V primers are modified from an original set of primers disclosed in Van Dongen et al, U.S. patent publication 2006/0234234, which is incorporated herein by reference. The modified primers arc disclosed in Faham et al, U.S. patent publication 2010/0151471, which is also incorporated herein by reference. - The Illumina Genome Analyzer is used to sequence the amplicon produced by the above primers. A two-stage amplification is performed on messenger RNA transcripts (200), as illustrated in
FIGS. 2A-2B , the first stage employing the above primers and a second stage to add common primers for bridge amplification and sequencing. As shown inFIG. 2A , a primary PCR is performed using on one side a 20 by primer (202) whose 3′ end is 16 bases from the J/C junction (204) and which is perfectly complementary to Cβ1(203) and the two alleles of Cβ2. In the V region (206) of RNA transcripts (200), primer set (212) is provided which contains primer sequences complementary to the different V region sequences (34 in one embodiment). Primers of set (212) also contain a non-complementary tail (214) that produces amplicon (216) having primer binding site (218) specific for P7 primers (220). After a conventional multiplex PCR, amplicon (216) is formed that contains the highly diverse portion of the J(D)V region (206, 208, and 210) of the mRNA transcripts and common primer binding sites (203 and 218) for a secondary amplification to add a sample tag (221) and primers (220 and 222) for cluster formation by bridge PCR. In the secondary PCR, on the same side of the template, a primer (222 inFIG. 2B and referred to herein as “C10-17-P5”) is used that has at its 3′end the sequence of the 10 bases closest to the J/C junction, followed by 17 bp with the sequence of positions 15-31 from the J/C junction, followed by the PS sequence (224), which plays a role in cluster formation by bridge PCR in Solexa sequencing. (When the C10-17-P5 primer (222) anneals to the template generated from the first PCR, a 4 by loop (position 11-14) is created in the template, as the primer hybridizes to the sequence of the 10 bases closest to the J/C junction and bases at positions 15-31 from the J/C junction. The looping of positions 11-14 eliminates differential amplification of templates carrying Cβ1 or Cβ2. Sequencing is then done with a primer complementary to the sequence of the 10 bases closest to the J/C junction and bases at positions 15-31 from the J/C junction (this primer is called C′). C10-17-P5 primer can be HPLC purified in order to ensure that all the amplified material has intact ends that can be efficiently utilized in the cluster formation.) - In
FIG. 2A , the length of the overhang on the V primers (212) is preferably 14 bp. The primary PCR is helped with a shorter overhang (214). Alternatively, for the sake of the secondary PCR, the overhang in the V primer is used in the primary PCR as long as possible because the secondary PCR is priming from this sequence. A minimum size of overhang (214) that supports an efficient secondary PCR was investigated. Two series of V primers (for two different V segments) with overhang sires from 10 to 30 with 2 by steps were made. Using the appropriate synthetic sequences, the first PCR was performed with each of the primers in the series and gel electrophoresis was performed to show that all amplified. - As illustrated in
FIG. 2A , the primary PCR uses 34 different V primers (212) that anneal to V region (206) of RNA templates (200) and contain a common 14 by overhang on the 5′ tail. The 14 by is the partial sequence of one of the Illumina sequencing primers (termed theRead 2 primer). The secondary amplification primer (220) on the same side includes P7 sequence, a tag (221), andRead 2 primer sequence (223) (this primer is called Read2_tagX_P7). The P7 sequence is used for cluster formation. Read 2 primer and its complement are used for sequencing the V segment and the tag respectively. A set of 96 of these primers with tags numbered 1 through 96 arc created (sec below). These primers arc HPLC purified in order to ensure that all the amplified material has intact ends that can be efficiently utilized in the cluster formation. - As mentioned above, the second stage primer, C-10-17-P5 (222,
FIG. 2B ) has interrupted homology to the template generated in the first stage PCR. The efficiency of amplification using this primer has been validated. An alternative primer to C-10-17-P5, termed CsegP5, has perfect homology to the first stage C primer and a 5′ tail carrying P5. The efficiency of using C-10-17-P5 and CsegP5 in amplifying first stage. PCR templates was compared by performing real time PCR. In several replicates, it was found that PCR using the C-10-17-P5 primer had little or no difference in efficiency compared with PCR using the CsegP5 primer. - Amplicon (300) resulting from the 2-stage amplification illustrated in
FIGS. 2A-2B has the structure typically used with the Illumina sequencer as shown inFIG. 3A . Two primers that anneal to the outmost part of the molecule, Illumina primers P5 and P7 are used for solid phase amplification of the molecule (cluster formation). Three sequence reads arc done per molecule. The first read of 100 by is done with the C′ primer, which has a melting temperature that is appropriate for the Illumina sequencing process. The second read is 6 by long only and is solely for the purpose of identifying the sample tag. It is generated using a tag primer provided by the manufacturer (Illumina). The final read is theRead 2 primer, also provided by the manufacturer (Illumina). Using this primer, a 100 by read in the V segment is generated starting with the 1st PCR V primer sequence. - While the present invention has been described with reference to several particular example embodiments, those skilled in the art will recognize that many changes may be made thereto without departing from the spirit and scope of the present invention. The present invention is applicable to a variety of sensor implementations and other subject matter, in addition to those discussed above.
- Unless otherwise specifically defined herein, terms and symbols of nucleic acid chemistry, biochemistry, genetics, and molecular biology used herein follow those of standard treatises and texts in the field, e.g. Kornberg and Baker, DNA Replication, Second Edition (W. H. Freeman, New York, 1992); Lehninger, Biochemistry, Second Edition (Worth Publishers, New York, 1975); Strachan and Read, Human Molecular Genetics, Second Edition (Wiley-Liss, New York, 1999); Abbas et al, Cellular and Molecular Immunology, 6th edition (Saunders, 2007).
- “Aligning” means a method of comparing a test sequence, such as a sequence read, to one or more reference sequences to determine which reference sequence or which portion of a reference sequence is closest based on some sequence distance measure. An exemplary method of aligning nucleotide sequences is the Smith Waterman algorithm. Distance measures may include Hamming distance, Levenshtein distance, or the like. Distance measures may include a component related to the quality values of nucleotides of the sequences being compared.
- “Amplicon” means the product of a polynucleotide amplification reaction; that is, a clonal population of polynucleotides, which may be single stranded or double stranded, which are, replicated from one or more starting sequences. The one or more starting sequences may be one or more copies of the same sequence, or they may be a mixture of different sequences. Preferably, amplicons are formed by the amplification of a single starting sequence. Amplicons may be produced by a variety of amplification reactions whose products comprise replicates of the one or more starting, or target, nucleic acids. In one aspect, amplification reactions producing amplicons arc “template-driven” in that base pairing of reactants, either nucleotides or oligonucleotides, have complements in a template polynucleotide that are required for the creation of reaction products. In one aspect, template-driven reactions are primer extensions with a nucleic acid polymerase or oligonucleotide ligations with a nucleic acid ligase. Such reactions include, but arc not limited to, polymerase chain reactions (PCRs), linear polymerase reactions, nucleic acid sequence-based amplification (NASBAs), rolling circle amplifications, and the like, disclosed in the following references that are incorporated herein by reference: Mullis et al, U.S. Pat. Nos. 4,683,195; 4,965,188; 4,683,202; 4,800,159 (PCR); Gelfand ct al, U.S. Pat. No. 5,210,015 (real-time PCR with “taqman” probes); Wittwer et al, U.S. Pat. No. 6,174,670; Kacian et al, U.S. Pat. No. 5,399,491 (“NASBA”); Lizardi, U.S. Pat. No. 5,854,033; Aono et al, Japanese patent publ. JP 4-262799 (rolling circle amplification); and the like. In one aspect, amplicons of the invention are produced by PCRs. An amplification reaction may be a “real-time” amplification if a detection chemistry is available that permits a reaction product to be measured as the amplification reaction progresses, e.g. “real-time PCR” described below, or “real-time NASBA” as described in Leone et al, Nucleic Acids Research, 26: 2150-2155 (1998), and like references. As used herein, the term “amplifying” means performing an amplification reaction. A “reaction mixture” means a solution containing all the necessary reactants for performing a reaction, which may include, but not be limited to, buffering agents to maintain pH at a selected level during a reaction, salts, co-factors, scavengers, and the like.
- “Clonality” as used herein means a measure of the degree to which the distribution of clonotype abundances among clonotypes of a repertoire is skewed to a single or a few clonotypes. Roughly, clonality is an inverse measure of clonotype diversity. Many measures or statistics are available from ecology describing species-abundance relationships that may be used for clonality measures in accordance with the invention, e.g. Chapters 17 & 18, in Pielou, An Introduction to Mathematical Ecology, (Wiley-Interscience, 1969). In one aspect, a clonality measure used with the invention is a function of a clonotype profile (that is, the number of distinct clonotypes detected and their abundances), so that after a clonotype profile is measured, clonality may be computed from it to give a single number. One clonality measure is Simpson's measure, which is simply the probability that two randomly drawn clonotypes will be the same. Other clonality measures include information, based measures and McIntosh's diversity index, disclosed in Pielou (cited above).
- “Clonotype” means a recombined nucleotide sequence of a T cell or B cell encoding a T cell receptor (TCR) or B cell receptor (BCR); or a portion thereof. In one aspect, a collection of all the distinct clonotypes of a population of lymphocytes of an individual is a repertoire of such population, e.g. Arstila et al, Science, 286: 958-961 (1999); Yassai et al, Immunogenetics, 61: 493-502 (2009); Kedzierska et al, Mol. Immunol., 45(3): 607-618 (2008); and the like. As used herein, “clonotype profile,” or “repertoire profile,” is a tabulation of clonotypes of a sample of T cells and/or B cells (such as a peripheral blood sample containing such cells) that includes substantially all of the repertoire's clonotypes and their relative abundances. “Clonotype profile,” “repertoire profile,” and “repertoire” are used herein interchangeably. (That is, the term “repertoire,” as discussed more fully below, means a repertoire measured from a sample of lymphocytes). In one aspect of the invention, clonotypes comprise portions of an immunoglobulin heavy chain (IgH) or a TCR β chain. In other aspects of the invention, clonotypes may be based on other recombined molecules, such as immunoglobulin chains or TCRα chains, or portions thereof.
- “Coalescing” means treating two candidate clonotypes with sequence differences as the same by determining that such differences are due to experimental or measurement error and not due to genuine biological differences. In one aspect, a sequence of a higher frequency candidate clonotype is compared to that of a lower frequency candidate clonotype and if predetermined criteria are satisfied then the number of lower frequency candidate clonotypes is added to that of the higher frequency candidate clonotype and the lower frequency candidate clonotype is thereafter disregarded. That is, the read counts associated with the lower frequency candidate clonotype are added to those of the higher frequency candidate clonotype.
- “Complementarity determining regions” (CDRs) mean regions of an immunoglobulin (i.e., antibody) or T cell receptor where the molecule complements an antigen's conformation, thereby determining the molecule's specificity and contact with a specific antigen. T cell receptors and immunoglobulins each have three CDRs: CDR1 and CDR2 are found in the variable (V) domain, and CDR3 includes some of V, all of diverse (D) (heavy chains only) and joint (J), and some of the constant (C) domains.
- “Data structure” means an organization of information, usually in a computer or memory device, for better algorithm efficiency. Exemplary data structures include queues, stacks, linked lists, heaps, hash tables, arrays, trees, and the like. Data structures may have substructures that correspond to units of information or to subsets of related information. For example, arrays have rows and columns of entries; trees have nodes, branches, subtrees, and leaves; or the like. In one aspect, a data structure used herein is a sequence tree, an array or a hash table.
- “Lymphoid neoplasm” means an abnormal proliferation of lymphocytes that may be malignant or non-malignant. A lymphoid cancer is a malignant lymphoid neoplasm. Lymphoid neoplasms arc the result of, or are associated with, lymphoproliferative disorders, including but not limited to, follicular lymphoma, chronic lymphocytic leukemia (CLL), acute lymphocytic leukemia (ALL), hairy cell leukemia, lymphomas, multiple myeloma, post-transplant lymphoproliferative disorder, mantle lymphoma (MCL), diffuse large B cell lymphoma (DLBCL), T cell lymphoma, or the like; e.g. Jaffe et al, Blood, 112:4384-4399 (2008); Swerdlow et al, WHO Classification of Tumours of Haematopoietic and Lymphoid Tissues (e. 4th) (IARC Press, 2008).
- “Percent homologous,” “percent identical,” or like terms used in reference to the comparison of a reference sequence and another sequence (“comparison sequence”) mean that in an optimal alignment between the two sequences, the comparison sequence is identical to the reference sequence in a number of subunit positions equivalent to the indicated percentage, the subunits being nucleotides for polynucleotide comparisons or amino acids for polypeptide comparisons. As used herein, an “optimal alignment” of sequences being compared is one that maximizes matches between subunits and minimizes the number of gaps employed in constructing an alignment. Percent identities may be determined with commercially available implementations of algorithms, such as that described by Needleman and Wunsch, J. Mol. Biol., 48: 443-453 (1970) (“GAP” program of Wisconsin Sequence Analysis Package, Genetics Computer Group, Madison, Wis.), or the like. Other software packages in the art for constructing alignments and calculating percentage identity or other measures of similarity include the “BestFit” program, based'on the algorithm of Smith and Waterman, Advances in Applied Mathematics, 2: 482-489 (1081) (Wisconsin Sequence Analysis Package, Genetics Computer Group, Madison, Wis.). In other words, for example, to obtain a polynucleotide having a nucleotide sequence at least 95 percent identical to a reference nucleotide sequence, up to five percent of the nucleotides in the reference sequence may be deleted or substituted with another nucleotide, or a number of nucleotides up to five percent of the total number of nucleotides in the reference sequence may be inserted into the reference sequence.
- “Polymerase chain reaction,” or “PCR,” means a reaction for the in vitro amplification of specific DNA sequences by the simultaneous primer extension of complementary strands of DNA. In other words, PCR is a reaction for making multiple copies or replicates of a target nucleic acid flanked by primer binding sites, such reaction comprising one or more repetitions of the following steps: (i) denaturing the target nucleic acid, (ii) annealing primers to the primer binding sites, and (iii) extending the primers by a nucleic acid polymerase in the presence of nucleoside triphosphates. Usually, the reaction is cycled through different temperatures optimized for each step in a thermal cycler instrument. Particular temperatures, durations at each step, and rates of change between steps depend on many factors well-known to those of ordinary skill in the art, e.g. exemplified by the references: McPherson et al, editors, PCR: A Practical Approach and PCR2: A Practical Approach (IRL Press, Oxford, 1991 and 1995, respectively). For example, in a conventional PCR using Taq DNA polymerase, a double stranded target nucleic acid may be denatured at a temperature >90° C., primers annealed at a temperature in the range 50-75° C., and primers extended at a temperature in the range 72-78° C. The term “PCR” encompasses derivative forms of the reaction, including but not limited to, RT-PCR, real-time PCR, nested PCR, quantitative PCR, multiplexed PCR, and the like. Reaction volumes range from a few hundred nanoliters, e.g. 200 nL, to a few hundred μL, e.g. 200 μL. “Reverse transcription PCR,” or “RT-PCR,” means a PCR that is preceded by a reverse transcription reaction that converts a target RNA to a complementary single stranded DNA, which is then amplified, e.g. Tecott et al, U.S. Pat. No. 5,168,038, which patent is incorporated herein by reference. “Real-time PCR” means a PCR for which the amount of reaction product, i.e. amplicon, is monitored as the reaction proceeds. There are many forms of real-time PCR that differ mainly in the detection chemistries used for monitoring the reaction product, e.g. Gelfand et al, U.S. Pat. No. 5,210,015 (“taqman”); Wittwer et al, U.S. Pat. Nos. 6,174,670 and 6,569,627 (intercalating dyes); Tyagi et al, U.S. Pat. No. 5,925,517 (molecular beacons); which patents are incorporated herein by reference. Detection chemistries for real-time PCR are reviewed in Mackay et al, Nucleic Acids Research, 30: 1292-1305 (2002), which is also incorporated herein by reference. “Nested PCR” means a two-stage PCR wherein the amplicon of a first PCR becomes the sample for a second PCR using a new set of primers, at least one of which binds to an interior location of the first amplicon. As used herein, “initial primers” in reference to a nested amplification reaction mean the printers used to generate a first amplicon, and “secondary primers” mean the one or more primers used to generate a second, or nested, amplicon. “Multiplexed PCR” means a PCR wherein multiple target sequences (or a single target sequence and one or more reference sequences) are simultaneously carried out in the same reaction mixture, e.g. Bernard ct al, Anal. Biochem., 273: 221-228 (1999) (two-color real-time PCR). Usually, distinct sets of primers are employed for each sequence being amplified. Typically, the number of target sequences in a multiplex PCR is in the range of from 2 to 50, or from 2 to 40, or from 2 to 30. “Quantitative PCR” means a PCR designed to measure the abundance of one or more specific target sequences in a sample or specimen. Quantitative PCR includes both absolute quantitation and relative quantitation of such target sequences. Quantitative measurements are made using one or more reference sequences or internal standards that may be assayed separately or together with a target sequence. The reference sequence may be endogenous or exogenous to a sample or specimen, and in the latter case, may comprise one or more competitor templates. Typical endogenous reference sequences include segments of transcripts of the following genes: β-actin, GAPDH, β2-microglobulin, ribosomal RNA, and the like. Techniques for quantitative PCR are well-known to those of ordinary skill in the art, as exemplified in the following references that are incorporated by reference: Freeman et al, Biotechniques, 26: 112-126 (1999); Becker-Andre et al, Nucleic Acids Research, 17: 9437-9447 (1989); Zimmerman et al, Biotechniques, 21: 268-279 (1996); Diviacco et al, Gene, 122: 3013-3020 (1992); Becker-Andre et al, Nucleic Acids Research, 17: 9437-9446 (1989); and the like.
- “Primer” means an oligonucleotide, either natural or synthetic that is capable, upon forming a duplex with a polynucleotide template, of acting as a point of initiation of nucleic acid synthesis and being extended from its 3′ end along the template so that an extended duplex is formed. Extension of a primer is usually carried out with a nucleic acid polymerase, such as a DNA or RNA polymerase. The sequence of nucleotides added in the extension process is determined by the sequence of the template polynucleotide. Usually primers arc extended by a DNA polymerase. Primers usually have a length in the range of from 14 to 40 nucleotides, or in the range of from 18 to 36 nucleotides. Primers are employed in a variety of nucleic amplification reactions, for example, linear amplification reactions using a single primer, or polymerase chain reactions, employing two or more primers. Guidance for selecting the lengths and sequences of primers for particular applications is well known to those of ordinary skill in the art, as evidenced by the following references that are incorporated by reference: Dieffenbach, editor, PCR Primer: A Laboratory Manual, 2nd Edition (Cold Spring Harbor Press, New York, 2003).
- “Quality score” means a measure of the probability that a base assignment at a particular sequence location is correct. A variety methods are well known to those of ordinary skill for calculating quality scores for particular circumstances, such as, for bases called as a result of different sequencing chemistries, detection systems, base-calling algorithms, and so on. Generally, quality score values arc monotonically related to probabilities of correct base calling. For example, a quality score, or Q, of 10 may mean that there is a 90 percent chance that a base is called correctly, a Q of 20 may mean that there is a 99 percent chance that a base is called correctly, and so on. For some sequencing platforms, particularly those using sequencing-by-synthesis chemistries, average quality scores decrease as a function of sequence read length, so that quality scores at the beginning of a sequence read are higher than those at the end of a sequence read, such declines being due to phenomena such as incomplete extensions, carry forward extensions, loss of template, loss of polymerase, capping failures, deprotection failures, and the like.
- “Repertoire”, or “immune repertoire”, means a set of distinct recombined nucleotide sequences that encode T cell receptors (TCRs) or B cell receptors (BCRs), or fragments thereof, respectively; in a population Oft lymphocytes of an individual, wherein the nucleotide sequences of the set have a one-to-one correspondence with distinct lymphocytes or their subpopulations for substantially all of the lymphocytes of the population. In one aspect, a population of lymphocytes from which a repertoire is determined is taken from one or more tissue samples, such as one or more blood samples. A member nucleotide sequence of a repertoire is referred to herein as a “clonotype.” In one aspect, clonotypes of a repertoire comprises any segment of nucleic acid common to a T cell or a B cell population which has undergone somatic recombination during the development of TCRs or BCRs, including normal or aberrant (e.g. associated with cancers) precursor molecules thereof, including, but not limited to, any of the following: an immunoglobulin heavy chain (IgH) or subsets thereof (e.g. an IgH variable region, CDR3 region, or the like), incomplete IgH molecules, an immunoglobulin light chain or subsets thereof (e.g. a variable region, CDR region, or the like), T cell receptor a chain or subsets thereof, T cell receptor β chain or subsets thereof (e.g. variable region, CDR3, V(D)J region, or the like), a CDR (including CDR1, CDR2 or CDR3, of either TCRs or BCRs, or combinations of such CDRs), V(D)J regions of either TCRs or BCRs, hypermutated regions of IgH variable regions. or the like. In one aspect, nucleic acid segments defining clonotypes of a repertoire arc selected so that their diversity (i.e. the number of distinct nucleic acid sequences in the set) is large enough so that substantially every T cell or B cell or clone thereof in an individual carries a unique nucleic acid sequence of such repertoire. That is, in accordance with the invention, a practitioner may select for defining clonotypes a particular segment or region of recombined nucleic acids that encode TCRs or BCRs that do not reflect the full diversity of a population of T cells or B cells; however, preferably, clonotypes are defined so that they do reflect the diversity of the population of T cells and/or B cells from which they are derived. That is, preferably each different clone of a sample has different clonotype. (Of course, in some applications, there will be multiple copies of one or more particular clonotypes within a profile, such as in the case of samples from leukemia or lymphoma patients). In other aspects of the invention, the population of lymphocytes corresponding to a repertoire may be circulating B cells, or may be circulating T cells, or may be subpopulations of either of the foregoing populations, including but not limited to, CD4+ T cells, or CD8+ T cells, or other subpopulations defined by cell surface markers, or the like. Such subpopulations may be acquired by taking samples from particular tissues, e.g. bone marrow, or lymph nodes, or the like, or by sorting or enriching cells from a sample (such as peripheral blood) based on one or more cell surface markers, size, morphology, or the like. In still other aspects, the population of lymphocytes corresponding to a repertoire may be derived from disease tissues, such as a tumor tissue, an infected tissue, or the like. In one embodiment, a repertoire comprising human TCR β chains or fragments thereof comprises a number of distinct nucleotide sequences in the range of from 0.1×106to 1.8×106, or in the range of from 0.5×106 to 1.5×106, or in the range of from 0.8×106 to 1.2×106. In another embodiment, a repertoire comprising human. IgH chains or fragments thereof comprises a number of distinct nucleotide sequences in the range of from 0.1×106 to 1.8×106, or in the range of from 0.5×106to 1.5×106, or in the range of from 0.8×106to 1.2×106. In a particular embodiment, a repertoire of the invention comprises a set of nucleotide sequences encoding substantially all segments of the V(D)J region of an IgH chain. In one aspect, “substantially all” as used herein means every segment having a relative abundance of 0.001 percent or higher; or in another aspect, “substantially all” as used herein means every segment having a relative abundance of 0.0001 percent or higher. In another particular embodiment, a repertoire of the invention comprises a set of nucleotide sequences that encodes substantially all segments of the V(D)J region of a TCR chain. In another embodiment, a repertoire of the invention comprises a set of nucleotide sequences having lengths in the range of from 25-200 nucleotides and including segments of the V, D, and J regions of a TCR β chain. In another embodiment, a repertoire of the invention comprises a set of nucleotide sequences having lengths in the range of from 25-200 nucleotides and including segments of the V, D, and J regions of an IgH chain. In another embodiment, a repertoire of the invention comprises a number of distinct nucleotide sequences that is substantially equivalent to the number of lymphocytes expressing a distinct IgH chain. In another embodiment, a repertoire of the invention comprises a number of distinct nucleotide sequences that is substantially equivalent to the number of lymphocytes expressing a distinct TCR β chain. In still another embodiment, “substantially equivalent” means that with ninety-nine percent probability a repertoire of nucleotide sequences will include a nucleotide sequence encoding an IgH or TCR β or portion thereof carried or expressed by every lymphocyte of a population of an individual at a frequency of 0.001 percent or greater. In still another embodiment, “substantially equivalent” means that with ninety-nine percent probability a repertoire of nucleotide sequences will include a nucleotide sequence encoding an IgH or TCR β or portion thereof carried or expressed by every lymphocyte present at a frequency of 0.0001 percent or greater. The sets of clonotypes described in the foregoing two sentences are sometimes referred to herein as representing the “full repertoire” of IgH and/or TCRβ sequences. As mentioned above, when measuring or generating a clonotype profile (or repertoire profile), a sufficiently large sample of lymphocytes is obtained so that such profile provides a reasonably accurate representation of a repertoire for a particular application. In one aspect, samples comprising from 105 to 107 lymphocytes are employed, especially when obtained from peripheral blood samples of from 1-10 mL.
- “Sequence read” means a sequence of nucleotides determined from a Sequence or stream of data generated by a sequencing technique, which determination is made, for example, by means; of base-calling software associated with the technique, e.g. base-calling software from a commercial provider of a DNA sequencing platform. A sequence read usually includes quality scores for each nucleotide in the sequence. Typically, sequence reads are made by extending a primer along a template nucleic acid, e.g. with a DNA polymerase or a DNA ligase. Data is generated by recording signals, such as optical, chemical (e.g. pH change), or electrical signals, associated with such extension. Such initial data is converted into a sequence read.
- “Sequence tree” means a tree data structure for representing nucleotide sequences. In one aspect, a tree data structure of the invention is a rooted directed tree comprising nodes and edges that do not include cycles, or cyclical pathways. Edges from nodes of tree data structures of the invention are usually ordered. Nodes and/or edges arc structures that may contain, or be associated with, a value. Each node in a tree has zero or more child nodes, which by convention are shown below it in the tree. A node that has a child is called the child's parent node. A node has at most one parent. Nodes that do not have any children are called leaf nodes. The topmost node in a tree is called the root node. Being the topmost node, the root node will not have parents. It is the node at which operations on the tree commonly begin (although some algorithms begin with, the leaf nodes and work up ending at the root). All other nodes can be reached from it by following edges or links.
Claims (46)
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/196,885 US9043160B1 (en) | 2009-11-09 | 2011-08-02 | Method of determining clonotypes and clonotype profiles |
US14/611,878 US20150252419A1 (en) | 2009-11-09 | 2015-02-02 | Method of determining clonotypes and clonotype profiles |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/615,263 US8236503B2 (en) | 2008-11-07 | 2009-11-09 | Methods of monitoring conditions by sequence analysis |
US45574310P | 2010-10-25 | 2010-10-25 | |
US201161446822P | 2011-02-25 | 2011-02-25 | |
US13/100,389 US8691510B2 (en) | 2008-11-07 | 2011-05-04 | Sequence analysis of complex amplicons |
US13/196,885 US9043160B1 (en) | 2009-11-09 | 2011-08-02 | Method of determining clonotypes and clonotype profiles |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/100,389 Continuation-In-Part US8691510B2 (en) | 2008-11-07 | 2011-05-04 | Sequence analysis of complex amplicons |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/611,878 Division US20150252419A1 (en) | 2009-11-09 | 2015-02-02 | Method of determining clonotypes and clonotype profiles |
Publications (2)
Publication Number | Publication Date |
---|---|
US9043160B1 US9043160B1 (en) | 2015-05-26 |
US20150167080A1 true US20150167080A1 (en) | 2015-06-18 |
Family
ID=53176440
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/196,885 Active 2031-01-07 US9043160B1 (en) | 2009-11-09 | 2011-08-02 | Method of determining clonotypes and clonotype profiles |
US14/611,878 Abandoned US20150252419A1 (en) | 2009-11-09 | 2015-02-02 | Method of determining clonotypes and clonotype profiles |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/611,878 Abandoned US20150252419A1 (en) | 2009-11-09 | 2015-02-02 | Method of determining clonotypes and clonotype profiles |
Country Status (1)
Country | Link |
---|---|
US (2) | US9043160B1 (en) |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150275296A1 (en) * | 2014-04-01 | 2015-10-01 | Sequenta, Inc. | Determining antigen-specific t-cells |
US9279159B2 (en) | 2011-10-21 | 2016-03-08 | Adaptive Biotechnologies Corporation | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US9347099B2 (en) | 2008-11-07 | 2016-05-24 | Adaptive Biotechnologies Corp. | Single cell analysis by polymerase cycling assembly |
US9365901B2 (en) | 2008-11-07 | 2016-06-14 | Adaptive Biotechnologies Corp. | Monitoring immunoglobulin heavy chain evolution in B-cell acute lymphoblastic leukemia |
US9371558B2 (en) | 2012-05-08 | 2016-06-21 | Adaptive Biotechnologies Corp. | Compositions and method for measuring and calibrating amplification bias in multiplexed PCR reactions |
US9416420B2 (en) | 2008-11-07 | 2016-08-16 | Adaptive Biotechnologies Corp. | Monitoring health and disease status using clonotype profiles |
US9499865B2 (en) | 2011-12-13 | 2016-11-22 | Adaptive Biotechnologies Corp. | Detection and measurement of tissue-infiltrating lymphocytes |
US9506119B2 (en) | 2008-11-07 | 2016-11-29 | Adaptive Biotechnologies Corp. | Method of sequence determination using sequence tags |
US9512487B2 (en) | 2008-11-07 | 2016-12-06 | Adaptive Biotechnologies Corp. | Monitoring health and disease status using clonotype profiles |
US9528160B2 (en) | 2008-11-07 | 2016-12-27 | Adaptive Biotechnolgies Corp. | Rare clonotypes and uses thereof |
US9708657B2 (en) | 2013-07-01 | 2017-07-18 | Adaptive Biotechnologies Corp. | Method for generating clonotype profiles using sequence tags |
US9809813B2 (en) | 2009-06-25 | 2017-11-07 | Fred Hutchinson Cancer Research Center | Method of measuring adaptive immunity |
US9824179B2 (en) | 2011-12-09 | 2017-11-21 | Adaptive Biotechnologies Corp. | Diagnosis of lymphoid malignancies and minimal residual disease detection |
US10077478B2 (en) | 2012-03-05 | 2018-09-18 | Adaptive Biotechnologies Corp. | Determining paired immune receptor chains from frequency matched subunits |
US10150996B2 (en) | 2012-10-19 | 2018-12-11 | Adaptive Biotechnologies Corp. | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US10221461B2 (en) | 2012-10-01 | 2019-03-05 | Adaptive Biotechnologies Corp. | Immunocompetence assessment by adaptive immune receptor diversity and clonality characterization |
US10246701B2 (en) | 2014-11-14 | 2019-04-02 | Adaptive Biotechnologies Corp. | Multiplexed digital quantitation of rearranged lymphoid receptors in a complex mixture |
US10323276B2 (en) | 2009-01-15 | 2019-06-18 | Adaptive Biotechnologies Corporation | Adaptive immunity profiling and methods for generation of monoclonal antibodies |
US10385475B2 (en) | 2011-09-12 | 2019-08-20 | Adaptive Biotechnologies Corp. | Random array sequencing of low-complexity libraries |
US10392663B2 (en) | 2014-10-29 | 2019-08-27 | Adaptive Biotechnologies Corp. | Highly-multiplexed simultaneous detection of nucleic acids encoding paired adaptive immune receptor heterodimers from a large number of samples |
US10428325B1 (en) | 2016-09-21 | 2019-10-01 | Adaptive Biotechnologies Corporation | Identification of antigen-specific B cell receptors |
US11041202B2 (en) | 2015-04-01 | 2021-06-22 | Adaptive Biotechnologies Corporation | Method of identifying human compatible T cell receptors specific for an antigenic target |
US11047008B2 (en) | 2015-02-24 | 2021-06-29 | Adaptive Biotechnologies Corporation | Methods for diagnosing infectious disease and determining HLA status using immune repertoire sequencing |
US11066705B2 (en) | 2014-11-25 | 2021-07-20 | Adaptive Biotechnologies Corporation | Characterization of adaptive immune response to vaccination or infection using immune repertoire sequencing |
WO2021173502A1 (en) * | 2020-02-28 | 2021-09-02 | 10X Genomics, Inc. | Systems and methods for identifying adaptive immune cell clonotypes |
US11248253B2 (en) | 2014-03-05 | 2022-02-15 | Adaptive Biotechnologies Corporation | Methods using randomer-containing synthetic molecules |
US11254980B1 (en) | 2017-11-29 | 2022-02-22 | Adaptive Biotechnologies Corporation | Methods of profiling targeted polynucleotides while mitigating sequencing depth requirements |
US11390921B2 (en) | 2014-04-01 | 2022-07-19 | Adaptive Biotechnologies Corporation | Determining WT-1 specific T cells and WT-1 specific T cell receptors (TCRs) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111755075B (en) * | 2019-03-28 | 2023-09-29 | 深圳华大生命科学研究院 | Method for filtering sequence pollution among high-throughput sequencing samples of immune repertoire |
Family Cites Families (127)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1990004648A1 (en) | 1988-10-20 | 1990-05-03 | Alexander Alan Morley | Method for diagnosis of monoclonality in leukaemia and lymphoma |
US5336598A (en) | 1989-11-15 | 1994-08-09 | National Jewish Center For Immunology And Respiratory Medicine | Method for diagnosing a superantigen caused pathologial condition via assay of T-cells |
US5298396A (en) | 1989-11-15 | 1994-03-29 | National Jewish Center For Immunology And Respiratory Medicine | Method for identifying T cells disease involved in autoimmune disease |
US5635354A (en) | 1991-01-09 | 1997-06-03 | Institut National De La Sante Et De La Recherche Medicale (Inserm) | Method for describing the repertoires of antibodies (Ab) and of T-cell receptors (TcR) of an individual's immune system |
US5837447A (en) | 1992-04-15 | 1998-11-17 | Blood Center Research Foundation, Inc., The | Monitoring an immune response by analysis of amplified immunoglobulin or T-cell-receptor nucleic acid |
EP0753060A4 (en) | 1994-04-18 | 1999-11-24 | New York Society | Conserved t-cell receptor sequences |
US6090592A (en) | 1994-08-03 | 2000-07-18 | Mosaic Technologies, Inc. | Method for performing amplification of nucleic acid on supports |
US5698396A (en) | 1995-06-07 | 1997-12-16 | Ludwig Institute For Cancer Research | Method for identifying auto-immunoreactive substances from a subject |
US6087096A (en) | 1995-11-13 | 2000-07-11 | Dau; Peter C. | Method of intrafamily fragment analysis of the T cell receptor α and β chain CDR3 regions |
KR19990087749A (en) | 1996-03-13 | 1999-12-27 | 시오노 요시히코 | Human T Cell Clones Specific to Rheumatoid Arthritis |
WO1997046706A1 (en) | 1996-06-03 | 1997-12-11 | University Of Alberta | Methods for detection of rearranged dna |
WO1998001738A2 (en) | 1996-06-20 | 1998-01-15 | Cornell Research Foundation, Inc. | Identification of abnormalities in the expression of t and b cell antigen receptors as disease indicators |
PT937251E (en) | 1996-09-06 | 2007-01-31 | Ortho Mcneil Pharm Inc | Purification of antigen-specific t cells |
GB9626815D0 (en) | 1996-12-23 | 1997-02-12 | Cemu Bioteknik Ab | Method of sequencing DNA |
US6210910B1 (en) | 1998-03-02 | 2001-04-03 | Trustees Of Tufts College | Optical fiber biosensor array comprising cell populations confined to microcavities |
DE19844931C1 (en) | 1998-09-30 | 2000-06-15 | Stefan Seeger | Procedures for DNA or RNA sequencing |
AUPQ687600A0 (en) | 2000-04-13 | 2000-05-11 | Flinders Technologies Pty Ltd | A method of detection |
US6596492B2 (en) | 2000-07-11 | 2003-07-22 | Colorado State University Research Foundation | PCR materials and methods useful to detect canine and feline lymphoid malignancies |
US7567870B1 (en) | 2000-07-31 | 2009-07-28 | Institute For Systems Biology | Multiparameter analysis for predictive medicine |
US6720144B1 (en) | 2001-06-21 | 2004-04-13 | Quest Diagnostics | Detection of clonal T-cell receptor-γ gene rearrangement by PCR/temporal temperature gradient gel electrophoresis (TTGE) |
US20030096277A1 (en) | 2001-08-30 | 2003-05-22 | Xiangning Chen | Allele specific PCR for genotyping |
EP1436404B1 (en) | 2001-09-19 | 2009-11-11 | Alexion Pharmaceuticals, Inc. | Engineered templates and their use in single primer amplification |
JP5236856B2 (en) | 2001-11-09 | 2013-07-17 | ライフ テクノロジーズ コーポレーション | Identification, observation and treatment of diseases and identification of biological states using gene expression profiles |
GB0128153D0 (en) | 2001-11-23 | 2002-01-16 | Bayer Ag | Profiling of the immune gene repertoire |
WO2003059155A2 (en) | 2002-01-09 | 2003-07-24 | Maxx Genetech Co. Ltd | Method of detecting t-cell proliferation for diagnosis of diseases by gene array |
WO2004046098A2 (en) | 2002-05-16 | 2004-06-03 | Vanderbilt University | Method for predicting autoimmune diseases |
AU2003247129A1 (en) | 2002-07-01 | 2004-01-19 | Centre National De La Recherche Scientifique (Cnrs) | System, method, device, and computer program product for extraction, gathering, manipulation, and analysis of peak data from an automated sequencer |
EP1517995A4 (en) | 2002-07-03 | 2007-07-11 | Inst Scient Res Inc | Compositions and methods for the detection of human t cell receptor variable family gene expression |
AU2003261168A1 (en) | 2002-07-19 | 2004-02-09 | Althea Technologies, Inc. | Strategies for gene expression analysis |
US20050009096A1 (en) | 2002-10-11 | 2005-01-13 | The Regents Of The University Of California | Method for diagnosis and prognosis of multiple sclerosis |
DK1549764T3 (en) | 2002-10-11 | 2011-01-24 | Univ Erasmus | Primers for nucleic acid amplification in PCR-based clonality studies |
US20060147925A1 (en) | 2002-11-13 | 2006-07-06 | Morley Alexander A | Method of detection |
ATE536548T1 (en) | 2002-11-14 | 2011-12-15 | Vivalis | MICROWELL ARRAY CHIP FOR DETECTING AN ANTIGEN-SPECIFIC LYMPHOCYTE, METHOD FOR DETECTING AN ANTIGEN-SPECIFIC LYMPHOCYTE AND METHOD FOR CLONING AN ANTIGEN RECEPTOR GENE OF AN ANTIGEN-SPECIFIC LYMPHOCYTE |
WO2004063706A2 (en) | 2003-01-08 | 2004-07-29 | Maxx Genetech Co., Ltd. | Method of detecting over-expression of t-cell receptor genes by real-time pcr |
WO2004096985A2 (en) | 2003-04-24 | 2004-11-11 | Mayo Foundation For Medical Education And Research | Methods for assessing biologic diversity |
US20070105105A1 (en) | 2003-05-23 | 2007-05-10 | Mount Sinai School Of Medicine Of New York University | Surrogate cell gene expression signatures for evaluating the physical state of a subject |
US20050010030A1 (en) | 2003-07-02 | 2005-01-13 | Zang Jingwu Z. | T cell receptor CDR3 sequence and methods for detecting and treating rheumatoid arthritis |
TWI333977B (en) | 2003-09-18 | 2010-12-01 | Symphogen As | Method for linking sequences of interest |
FR2863274B1 (en) | 2003-12-05 | 2012-11-16 | Commissariat Energie Atomique | METHOD FOR QUANTITATIVE EVALUATION OF REARRANGEMENT OR TARGETED GENETIC RECOMBINATION OF AN INDIVIDUAL AND APPLICATIONS THEREOF |
WO2005053603A2 (en) | 2003-12-08 | 2005-06-16 | Yeda Research And Development Co. Ltd. | Antigen receptor variable region typing |
US20080166718A1 (en) | 2003-12-15 | 2008-07-10 | Institut Pasteur | Repertoire determination of a lymphocyte B population |
ATE421998T1 (en) | 2003-12-15 | 2009-02-15 | Pasteur Institut | DETERMINATION OF THE REPERTOIRE OF B-LYMPHOCYTE POPULATIONS |
EP1704251A1 (en) | 2003-12-15 | 2006-09-27 | Institut Pasteur | Repertoire determination of a lymphocyte b population |
US20060046258A1 (en) | 2004-02-27 | 2006-03-02 | Lapidus Stanley N | Applications of single molecule sequencing |
WO2005084134A2 (en) | 2004-03-04 | 2005-09-15 | Dena Leshkowitz | Quantifying and profiling antibody and t cell receptor gene expression |
DE102004016437A1 (en) | 2004-04-04 | 2005-10-20 | Oligene Gmbh | Method for detecting signatures in complex gene expression profiles |
JP2007536939A (en) | 2004-05-14 | 2007-12-20 | アモークス インコーポレーティッド | Immune cell biosensor and method of use thereof |
US20060094018A1 (en) | 2004-08-03 | 2006-05-04 | Bauer A R Jr | Discovery and a method for the early detection of pancreatic cancer and other disease conditions |
US7820382B2 (en) | 2004-08-03 | 2010-10-26 | Bauer A Robert | Method for the early detection of breast cancer, lung cancer, pancreatic cancer and colon polyps, growths and cancers as well as other gastrointestinal disease conditions and the preoperative and postoperative monitoring of transplanted organs from the donor and in the recipient and their associated conditions related and unrelated to the organ transplantation |
WO2006031221A1 (en) | 2004-09-13 | 2006-03-23 | Government Of The United States Of America, Represented By The Secretary, Department Of Health And Human Services | Compositions comprising t cell receptors and methods of use thereof |
FR2881436B1 (en) | 2005-02-03 | 2007-04-27 | Commissariat Energie Atomique | METHOD FOR DETERMINING THE DIVERSITY OF T LYMPHOCYTES IN A BIOLOGICAL SAMPLE |
US20080280774A1 (en) | 2005-02-16 | 2008-11-13 | Wyeth | Methods and Systems for Diagnosis, Prognosis and Selection of Treatment of Leukemia |
EP1877576B1 (en) | 2005-04-12 | 2013-01-23 | 454 Life Sciences Corporation | Methods for determining sequence variants using ultra-deep sequencing |
WO2006116155A2 (en) | 2005-04-21 | 2006-11-02 | The Regents Of The University Of California | A method for diagnosis and prognosis of multiple sclerosis subtypes |
US7375211B2 (en) | 2005-11-18 | 2008-05-20 | Kou Zhong C | Method for detection and quantification of T-cell receptor Vβ repertoire |
MX2008011280A (en) | 2006-03-06 | 2008-09-12 | Symphogen As | Recombinant polyclonal antibody for treatment of respiratory syncytial virus infections. |
PT2963127T (en) | 2006-04-04 | 2017-10-06 | Keygene Nv | High throughput detection of molecular markers based on restriction fragments |
WO2007136518A2 (en) | 2006-05-17 | 2007-11-29 | Torrey Pines Institute For Molecular Studies | Treatment of autoimmune disorders |
US7553949B2 (en) | 2006-06-30 | 2009-06-30 | Applied Biosystems, Llc | Reversible terminator nucleotides and methods of use |
WO2008108803A2 (en) | 2006-07-13 | 2008-09-12 | Amaox, Ltd. | Immune cell biosensors and methods of using same |
US20080274904A1 (en) | 2006-08-10 | 2008-11-06 | Niall Anthony Gormley | Method of target enrichment |
WO2008026927A2 (en) | 2006-08-30 | 2008-03-06 | Academisch Medisch Centrum | Process for displaying t- and b-cell receptor repertoires |
WO2008039694A2 (en) | 2006-09-26 | 2008-04-03 | St. Jude Children's Research Hospital | Methods and compositions for monitoring t cell receptor diversity |
JP5026047B2 (en) | 2006-10-18 | 2012-09-12 | 株式会社膠原病研究所 | Method for identifying self-responsive T cells or T cell receptors involved in the development of autoimmune diseases and use thereof |
US8148067B2 (en) | 2006-11-09 | 2012-04-03 | Xdx, Inc. | Methods for diagnosing and monitoring the status of systemic lupus erythematosus |
JP5669397B2 (en) | 2007-03-01 | 2015-02-12 | シムフォゲン・アクティーゼルスカブSymphogen A/S | Methods for cloning cognate antibodies |
WO2008147879A1 (en) | 2007-05-22 | 2008-12-04 | Ryan Golhar | Automated method and device for dna isolation, sequence determination, and identification |
US8454906B2 (en) | 2007-07-24 | 2013-06-04 | The Regents Of The University Of California | Microfabricated droplet generator for single molecule/cell genetic analysis in engineered monodispersed emulsions |
ITRM20070429A1 (en) | 2007-08-06 | 2009-02-07 | Uni Cattolica Del Sacro Cuor E | MEANS FOR DIAGNOSIS PREVENTION AND CARE OF RHEUMATOID ARTHRITIS. |
JP5547071B2 (en) | 2007-08-09 | 2014-07-09 | セルラ・インコーポレイテッド | Method and apparatus for associating multi-parameter single cell measurements and recovery of residual biological material |
WO2009045898A2 (en) | 2007-09-28 | 2009-04-09 | Mayo Foundation For Medical Education And Research | Assessing t cell repertoires |
WO2009070767A2 (en) | 2007-11-28 | 2009-06-04 | Whitehead Institute For Biomedical Research | Systemic instigation systems to study tumor growth or metastasis |
EP2062982A1 (en) | 2007-11-26 | 2009-05-27 | ImmunID | Method for studying the V(D)J combinatorial diversity. |
US7767400B2 (en) | 2008-02-03 | 2010-08-03 | Helicos Biosciences Corporation | Paired-end reads in sequencing by synthesis |
JP2011516033A (en) | 2008-02-28 | 2011-05-26 | ジ・オハイオ・ステイト・ユニバーシティ・リサーチ・ファウンデイション | MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid leukemia (AML) and uses thereof |
JP2011517283A (en) | 2008-02-28 | 2011-06-02 | ジ・オハイオ・ステイト・ユニバーシティ・リサーチ・ファウンデイション | MicroRNA-based methods and compositions for prognosis and treatment of prostate related disorders |
US20090226975A1 (en) | 2008-03-10 | 2009-09-10 | Illumina, Inc. | Constant cluster seeding |
EP2100970B1 (en) | 2008-03-13 | 2017-05-10 | National Institute of Immunology | Ig genes specific oligonucleotides and uses thereof |
US9012148B2 (en) | 2008-04-16 | 2015-04-21 | Jian Han | Method for evaluating and comparing immunorepertoires |
JP2012500964A (en) | 2008-05-09 | 2012-01-12 | デューク ユニバーシティ | Autoantibodies in cancer detection and treatment |
US9068181B2 (en) | 2008-05-23 | 2015-06-30 | The General Hospital Corporation | Microfluidic droplet encapsulation |
CN102112624A (en) | 2008-05-30 | 2011-06-29 | 麻省理工学院 | Compositions and methods for spatial separation and screening of cells |
WO2009151628A2 (en) | 2008-06-12 | 2009-12-17 | Gov't Of The Usa, As Represented By The Secretary, Department Of Health Human Services | Monitoring tcr-b to determine hiv therapy and disease progression |
TW201022492A (en) | 2008-06-25 | 2010-06-16 | Baylor Res Inst | Blood transcriptional signature of mycobacterium tuberculosis infection |
US8394583B2 (en) | 2008-07-24 | 2013-03-12 | The Board Of Regents Of The University Of Texas System | VH4 codon signature for multiple sclerosis |
AU2009297108B2 (en) | 2008-09-23 | 2015-02-12 | Bio-Rad Laboratories, Inc. | Droplet-based assay system |
US9156010B2 (en) | 2008-09-23 | 2015-10-13 | Bio-Rad Laboratories, Inc. | Droplet-based assay system |
TW201027105A (en) | 2008-10-31 | 2010-07-16 | Qualcomm Inc | Using magnetometer with a positioning system |
US9528160B2 (en) | 2008-11-07 | 2016-12-27 | Adaptive Biotechnolgies Corp. | Rare clonotypes and uses thereof |
US8748103B2 (en) | 2008-11-07 | 2014-06-10 | Sequenta, Inc. | Monitoring health and disease status using clonotype profiles |
EP2719774B8 (en) | 2008-11-07 | 2020-04-22 | Adaptive Biotechnologies Corporation | Methods of monitoring conditions by sequence analysis |
US9394567B2 (en) | 2008-11-07 | 2016-07-19 | Adaptive Biotechnologies Corporation | Detection and quantification of sample contamination in immune repertoire analysis |
US8691510B2 (en) | 2008-11-07 | 2014-04-08 | Sequenta, Inc. | Sequence analysis of complex amplicons |
US8628927B2 (en) | 2008-11-07 | 2014-01-14 | Sequenta, Inc. | Monitoring health and disease status using clonotype profiles |
US20140234835A1 (en) | 2008-11-07 | 2014-08-21 | Sequenta, Inc. | Rare clonotypes and uses thereof |
US9506119B2 (en) | 2008-11-07 | 2016-11-29 | Adaptive Biotechnologies Corp. | Method of sequence determination using sequence tags |
US9365901B2 (en) | 2008-11-07 | 2016-06-14 | Adaptive Biotechnologies Corp. | Monitoring immunoglobulin heavy chain evolution in B-cell acute lymphoblastic leukemia |
AU2010206843B2 (en) | 2009-01-20 | 2015-01-29 | The Board Of Trustees Of The Leland Stanford Junior University | Single cell gene expression for diagnosis, prognosis and identification of drug targets |
SG10201403451QA (en) | 2009-06-25 | 2014-09-26 | Hutchinson Fred Cancer Res | Method of measuring adaptive immunity |
GB201000375D0 (en) | 2010-01-09 | 2010-02-24 | Univ Cardiff | T Cell clonotypes |
WO2011106738A2 (en) | 2010-02-25 | 2011-09-01 | Fred Hutchinson Cancer Research Center | Use of tcr clonotypes as biomarkers for disease |
EP3456847A1 (en) | 2010-05-06 | 2019-03-20 | Adaptive Biotechnologies Corporation | Monitoring solid transplant rejection using clonotype profiles |
US9234240B2 (en) | 2010-05-07 | 2016-01-12 | The Board Of Trustees Of The Leland Stanford Junior University | Measurement and comparison of immune diversity by high-throughput sequencing |
DK2625295T3 (en) | 2010-10-08 | 2019-06-11 | Harvard College | HIGH-THROUGHPUT-IMMUNE SEQUENCING |
WO2012061832A1 (en) | 2010-11-05 | 2012-05-10 | Illumina, Inc. | Linking sequence reads using paired code tags |
US9193997B2 (en) | 2010-12-15 | 2015-11-24 | The Board Of Trustees Of The Leland Stanford Junior University | Measuring and monitoring of cell clonality |
WO2013036459A2 (en) | 2011-09-09 | 2013-03-14 | Sequenta, Inc. | Sequence-based measures of immune response |
US10385475B2 (en) | 2011-09-12 | 2019-08-20 | Adaptive Biotechnologies Corp. | Random array sequencing of low-complexity libraries |
US20140256592A1 (en) | 2011-10-11 | 2014-09-11 | Sequenta, Inc. | Determining responsiveness of autoimmune patients to dmard treatment |
ES2619713T3 (en) | 2011-11-04 | 2017-06-26 | Adaptive Biotechnologies Corporation | Clonotypes of T-cell receptors shared among patients with ankylosing spondylitis |
US20140336059A1 (en) | 2011-12-05 | 2014-11-13 | Sequenta, Inc. | Clonotypes as biometric specimen tags |
WO2013086462A1 (en) | 2011-12-09 | 2013-06-13 | Sequenta, Inc. | Method of measuring immune activation |
US20150031553A1 (en) | 2011-12-13 | 2015-01-29 | Sequenta, Inc. | Method of measuring immune activation |
US9499865B2 (en) | 2011-12-13 | 2016-11-22 | Adaptive Biotechnologies Corp. | Detection and measurement of tissue-infiltrating lymphocytes |
JP2015504661A (en) | 2011-12-20 | 2015-02-16 | シーケンタ インコーポレイテッド | Monitoring of transformation from follicular lymphoma to diffuse large B-cell lymphoma by immunorepertoire analysis |
EP3372694A1 (en) | 2012-03-05 | 2018-09-12 | Adaptive Biotechnologies Corporation | Determining paired immune receptor chains from frequency matched subunits |
WO2013134302A1 (en) | 2012-03-05 | 2013-09-12 | Sequenta, Inc. | Monitoring immune responsiveness to cancer vaccination |
WO2013155119A1 (en) | 2012-04-13 | 2013-10-17 | Sequenta, Inc. | Detection and quantitation of sample contamination in immune repertoire analysis |
WO2013158936A1 (en) | 2012-04-20 | 2013-10-24 | Sequenta, Inc | Monitoring immunoglobulin heavy chain evolution in b-cell acute lymphoblastic leukemia |
WO2013181428A2 (en) | 2012-05-31 | 2013-12-05 | Sequenta, Inc. | Predicting relapse of chronic lymphocytic leukemia patients treated by allogeneic stem cell transplantation |
US20130324422A1 (en) | 2012-06-04 | 2013-12-05 | Sequenta, Inc. | Detecting disease-correlated clonotypes from fixed samples |
EP2859123A4 (en) | 2012-06-11 | 2015-12-16 | Sequenta Inc | Method of sequence determination using sequence tags |
CA2878694A1 (en) | 2012-07-24 | 2014-01-30 | Sequenta, Inc. | Single cell analysis using sequence tags |
WO2014026031A1 (en) | 2012-08-10 | 2014-02-13 | Sequenta, Inc. | High sensitivity mutation detection using sequence tags |
AU2014218910A1 (en) | 2013-02-22 | 2015-08-20 | Adaptive Biotechnologies Corp. | Rare clonotypes and uses thereof |
US20140255944A1 (en) | 2013-03-08 | 2014-09-11 | Sequenta, Inc. | Monitoring treatment-resistant clones in lymphoid and myeloid neoplasms by relative levels of evolved clonotypes |
US20140255929A1 (en) | 2013-03-11 | 2014-09-11 | Sequenta, Inc. | Mosaic tags for labeling templates in large-scale amplifications |
-
2011
- 2011-08-02 US US13/196,885 patent/US9043160B1/en active Active
-
2015
- 2015-02-02 US US14/611,878 patent/US20150252419A1/en not_active Abandoned
Cited By (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9528160B2 (en) | 2008-11-07 | 2016-12-27 | Adaptive Biotechnolgies Corp. | Rare clonotypes and uses thereof |
US9523129B2 (en) | 2008-11-07 | 2016-12-20 | Adaptive Biotechnologies Corp. | Sequence analysis of complex amplicons |
US10155992B2 (en) | 2008-11-07 | 2018-12-18 | Adaptive Biotechnologies Corp. | Monitoring health and disease status using clonotype profiles |
US9365901B2 (en) | 2008-11-07 | 2016-06-14 | Adaptive Biotechnologies Corp. | Monitoring immunoglobulin heavy chain evolution in B-cell acute lymphoblastic leukemia |
US10519511B2 (en) | 2008-11-07 | 2019-12-31 | Adaptive Biotechnologies Corporation | Monitoring health and disease status using clonotype profiles |
US9416420B2 (en) | 2008-11-07 | 2016-08-16 | Adaptive Biotechnologies Corp. | Monitoring health and disease status using clonotype profiles |
US10246752B2 (en) | 2008-11-07 | 2019-04-02 | Adaptive Biotechnologies Corp. | Methods of monitoring conditions by sequence analysis |
US9506119B2 (en) | 2008-11-07 | 2016-11-29 | Adaptive Biotechnologies Corp. | Method of sequence determination using sequence tags |
US10760133B2 (en) | 2008-11-07 | 2020-09-01 | Adaptive Biotechnologies Corporation | Monitoring health and disease status using clonotype profiles |
US10865453B2 (en) | 2008-11-07 | 2020-12-15 | Adaptive Biotechnologies Corporation | Monitoring health and disease status using clonotype profiles |
US10266901B2 (en) | 2008-11-07 | 2019-04-23 | Adaptive Biotechnologies Corp. | Methods of monitoring conditions by sequence analysis |
US9347099B2 (en) | 2008-11-07 | 2016-05-24 | Adaptive Biotechnologies Corp. | Single cell analysis by polymerase cycling assembly |
US9512487B2 (en) | 2008-11-07 | 2016-12-06 | Adaptive Biotechnologies Corp. | Monitoring health and disease status using clonotype profiles |
US10323276B2 (en) | 2009-01-15 | 2019-06-18 | Adaptive Biotechnologies Corporation | Adaptive immunity profiling and methods for generation of monoclonal antibodies |
US9809813B2 (en) | 2009-06-25 | 2017-11-07 | Fred Hutchinson Cancer Research Center | Method of measuring adaptive immunity |
US11905511B2 (en) | 2009-06-25 | 2024-02-20 | Fred Hutchinson Cancer Center | Method of measuring adaptive immunity |
US11214793B2 (en) | 2009-06-25 | 2022-01-04 | Fred Hutchinson Cancer Research Center | Method of measuring adaptive immunity |
US10385475B2 (en) | 2011-09-12 | 2019-08-20 | Adaptive Biotechnologies Corp. | Random array sequencing of low-complexity libraries |
US9279159B2 (en) | 2011-10-21 | 2016-03-08 | Adaptive Biotechnologies Corporation | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US9824179B2 (en) | 2011-12-09 | 2017-11-21 | Adaptive Biotechnologies Corp. | Diagnosis of lymphoid malignancies and minimal residual disease detection |
US9499865B2 (en) | 2011-12-13 | 2016-11-22 | Adaptive Biotechnologies Corp. | Detection and measurement of tissue-infiltrating lymphocytes |
US10077478B2 (en) | 2012-03-05 | 2018-09-18 | Adaptive Biotechnologies Corp. | Determining paired immune receptor chains from frequency matched subunits |
US10894977B2 (en) | 2012-05-08 | 2021-01-19 | Adaptive Biotechnologies Corporation | Compositions and methods for measuring and calibrating amplification bias in multiplexed PCR reactions |
US10214770B2 (en) | 2012-05-08 | 2019-02-26 | Adaptive Biotechnologies Corp. | Compositions and method for measuring and calibrating amplification bias in multiplexed PCR reactions |
US9371558B2 (en) | 2012-05-08 | 2016-06-21 | Adaptive Biotechnologies Corp. | Compositions and method for measuring and calibrating amplification bias in multiplexed PCR reactions |
US11180813B2 (en) | 2012-10-01 | 2021-11-23 | Adaptive Biotechnologies Corporation | Immunocompetence assessment by adaptive immune receptor diversity and clonality characterization |
US10221461B2 (en) | 2012-10-01 | 2019-03-05 | Adaptive Biotechnologies Corp. | Immunocompetence assessment by adaptive immune receptor diversity and clonality characterization |
US10150996B2 (en) | 2012-10-19 | 2018-12-11 | Adaptive Biotechnologies Corp. | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US9708657B2 (en) | 2013-07-01 | 2017-07-18 | Adaptive Biotechnologies Corp. | Method for generating clonotype profiles using sequence tags |
US10526650B2 (en) | 2013-07-01 | 2020-01-07 | Adaptive Biotechnologies Corporation | Method for genotyping clonotype profiles using sequence tags |
US10077473B2 (en) | 2013-07-01 | 2018-09-18 | Adaptive Biotechnologies Corp. | Method for genotyping clonotype profiles using sequence tags |
US11248253B2 (en) | 2014-03-05 | 2022-02-15 | Adaptive Biotechnologies Corporation | Methods using randomer-containing synthetic molecules |
US10435745B2 (en) | 2014-04-01 | 2019-10-08 | Adaptive Biotechnologies Corp. | Determining antigen-specific T-cells |
US20150275296A1 (en) * | 2014-04-01 | 2015-10-01 | Sequenta, Inc. | Determining antigen-specific t-cells |
US10066265B2 (en) * | 2014-04-01 | 2018-09-04 | Adaptive Biotechnologies Corp. | Determining antigen-specific t-cells |
US11261490B2 (en) | 2014-04-01 | 2022-03-01 | Adaptive Biotechnologies Corporation | Determining antigen-specific T-cells |
US11390921B2 (en) | 2014-04-01 | 2022-07-19 | Adaptive Biotechnologies Corporation | Determining WT-1 specific T cells and WT-1 specific T cell receptors (TCRs) |
US10392663B2 (en) | 2014-10-29 | 2019-08-27 | Adaptive Biotechnologies Corp. | Highly-multiplexed simultaneous detection of nucleic acids encoding paired adaptive immune receptor heterodimers from a large number of samples |
US10246701B2 (en) | 2014-11-14 | 2019-04-02 | Adaptive Biotechnologies Corp. | Multiplexed digital quantitation of rearranged lymphoid receptors in a complex mixture |
US11066705B2 (en) | 2014-11-25 | 2021-07-20 | Adaptive Biotechnologies Corporation | Characterization of adaptive immune response to vaccination or infection using immune repertoire sequencing |
US11047008B2 (en) | 2015-02-24 | 2021-06-29 | Adaptive Biotechnologies Corporation | Methods for diagnosing infectious disease and determining HLA status using immune repertoire sequencing |
US11041202B2 (en) | 2015-04-01 | 2021-06-22 | Adaptive Biotechnologies Corporation | Method of identifying human compatible T cell receptors specific for an antigenic target |
US10428325B1 (en) | 2016-09-21 | 2019-10-01 | Adaptive Biotechnologies Corporation | Identification of antigen-specific B cell receptors |
US11254980B1 (en) | 2017-11-29 | 2022-02-22 | Adaptive Biotechnologies Corporation | Methods of profiling targeted polynucleotides while mitigating sequencing depth requirements |
WO2021173502A1 (en) * | 2020-02-28 | 2021-09-02 | 10X Genomics, Inc. | Systems and methods for identifying adaptive immune cell clonotypes |
Also Published As
Publication number | Publication date |
---|---|
US9043160B1 (en) | 2015-05-26 |
US20150252419A1 (en) | 2015-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9043160B1 (en) | Method of determining clonotypes and clonotype profiles | |
US10526650B2 (en) | Method for genotyping clonotype profiles using sequence tags | |
US9523129B2 (en) | Sequence analysis of complex amplicons | |
US8691510B2 (en) | Sequence analysis of complex amplicons | |
US10385475B2 (en) | Random array sequencing of low-complexity libraries | |
EP2959020B1 (en) | Method to select rare clonotypes | |
US20150031553A1 (en) | Method of measuring immune activation | |
US9365901B2 (en) | Monitoring immunoglobulin heavy chain evolution in B-cell acute lymphoblastic leukemia | |
US9528160B2 (en) | Rare clonotypes and uses thereof | |
US9394567B2 (en) | Detection and quantification of sample contamination in immune repertoire analysis | |
JP6085249B2 (en) | Sequence analysis of complex amplicons | |
US9506119B2 (en) | Method of sequence determination using sequence tags | |
US20140349883A1 (en) | Monitoring transformation of follicular lymphoma to diffuse large b-cell lymphoma by immune repertoire analysis | |
US9499865B2 (en) | Detection and measurement of tissue-infiltrating lymphocytes | |
US20140342360A1 (en) | Method of measuring immune activation | |
US20140356339A1 (en) | Sequence-based measures of immune response | |
US20140234835A1 (en) | Rare clonotypes and uses thereof | |
US20150299800A1 (en) | Monitoring mantle cell lymphoma clonotypes in peripheral blood after immunotransplant | |
US20140255944A1 (en) | Monitoring treatment-resistant clones in lymphoid and myeloid neoplasms by relative levels of evolved clonotypes | |
US20140256592A1 (en) | Determining responsiveness of autoimmune patients to dmard treatment | |
WO2013158936A1 (en) | Monitoring immunoglobulin heavy chain evolution in b-cell acute lymphoblastic leukemia |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SEQUENTA, INC., CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:MOORHEAD, MARTIN;FAHAM, MALEK;WILLIS, THOMAS;REEL/FRAME:026696/0222 Effective date: 20110802 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: SEQUENTA, LLC, WASHINGTON Free format text: MERGER;ASSIGNOR:SEQUENTA, INC.;REEL/FRAME:036503/0732 Effective date: 20150107 |
|
FEPP | Fee payment procedure |
Free format text: PAT HOLDER CLAIMS SMALL ENTITY STATUS, ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: LTOS); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
AS | Assignment |
Owner name: ADAPTIVE BIOTECHNOLOGIES CORP., WASHINGTON Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SEQUENTA, LLC;REEL/FRAME:037560/0788 Effective date: 20160119 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2551); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 4 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2552); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 8 |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |