WO2013155119A1 - Detection and quantitation of sample contamination in immune repertoire analysis - Google Patents
Detection and quantitation of sample contamination in immune repertoire analysis Download PDFInfo
- Publication number
- WO2013155119A1 WO2013155119A1 PCT/US2013/035857 US2013035857W WO2013155119A1 WO 2013155119 A1 WO2013155119 A1 WO 2013155119A1 US 2013035857 W US2013035857 W US 2013035857W WO 2013155119 A1 WO2013155119 A1 WO 2013155119A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- nucleic acids
- sequence
- recombined
- nucleic acid
- region
- Prior art date
Links
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/6848—Nucleic acid amplification reactions characterised by the means for preventing contamination or increasing the specificity or sensitivity of an amplification reaction
Definitions
- profiles of nucleic acids encoding immune molecules 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 (cited above); Boyd et al (cited above); He et al (cited above).
- sequence-based profiles are capable of much greater sensitivity than approaches based on size distributions of amplified CDR-encoding regions, sequence sampling by microarrays, hybridization kinetics curves from PCR amplicons, or other approaches, e.g. Morley et al, U.S.
- contaminating DNA may reduce the effective limit of detection in assays employing immune repertoire sequencing.
- Sources of contaminating DNA include assay reagents, equipment, operator handling, aersols, and the like, e.g. Urban et al, J. Forensic Sci., 45(6): 1307-1311 (2000); Kwok, pgs. 142-145, in Innis et al, Editors, PCR Protocols (Academic Press, 1990).
- Detection of minimal residual disease (MRD) of cancers is impacted by such contamination. Patients treated for many cancers often retain an MRD related to the cancer.
- clonotypes have been particularly useful in assessing MRD in leukemias and lymphomas, since such segment (referred to herein as "clonotypes") typically have unique sequences which may serve as molecular tags for their associated cancer cells, e.g. Van Dongen et al (cited above); Faham and Willis, U.S. patent publication 2011/0207134; and the like. Nevertheless, the sensitivity of such techniques is still limited by the presence of cross-over contamination from other individuals.
- the present invention is drawn to methods for detecting and quantifying contaminating nucleic acids in immune repertoire assays.
- 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 for determining a level of
- a tissue sample comprising T cells and/or B cells of an individual
- the method comprising the steps of: (a) obtaining a tissue sample from an individual, the tissue sample comprising nucleic acids from the individual and possibly from one or more other individuals, the nucleic acids including recombined nucleic acids from T cells and/or B cells and non-recombined nucleic acids; (b) generating a clonotype profile from nucleic acid of the tissue sample; (c) sequencing genetic markers at one or more genetic loci of the nucleic acids from the tissue sample to obtain proportions of the nucleic acid from the one or more other individuals and from the individual based on the genetic identification of the individual and the one or more other individuals from the genetic markers, wherein the genetic markers occur at loci present in both the recombined nucleic acids and the non-recombined nucleic acids; and (d) determining the level of contaminating nucleic acid in the clonotype profile as the
- the above method further includes the steps of (i) measuring total amount of said nucleic acids in said tissue sample; (ii) measuring total amount of said recombined nucleic acids in said tissue sample; (iii) sequencing genetic markers at one or more genetic loci in an excisable segment of nucleic acid present in said non-recombined nucleic acids and absent from said recombined nucleic acids to obtain proportions of said non-recombined nucleic acid from said individual and from said one or more other individuals; and (iv) determining said level of contaminating nucleic acid in said clonotype profile from a level of contaminating recombined nucleic acids in said nucleic acid of said tissue sample, the level of contaminating recombined nucleic acids being determined from the total amount of nucleic acids, the total amount of said recombined nucleic acids, said proportions of said nucleic acid from said individual and said one or more other individuals, and the proportions
- the invention is directed to a method for determining a limit of detection of a correlating clonotype in a clonotype profile of a tissue sample comprising T cells and/or B cells of an individual comprising the steps of: (a) obtaining a tissue sample from an individual, the tissue sample comprising nucleic acids from the individual and possibly from one or more other individuals, the nucleic acids including recombined nucleic acids from T cells and/or B cells and non-recombined nucleic acids; (b) generating a clonotype profile from nucleic acid of the tissue sample; (c) sequencing genetic markers at one or more genetic loci of the nucleic acids from the tissue sample to obtain proportions of the nucleic acid from the one or more other individuals and from the individual based on the genetic identification of the individual and the one or more other individuals from the genetic markers, wherein the genetic markers occur at loci present in both the recombined nucleic acids and the non-recombined
- FIG. 1 A diagrammatically illustrates a typical scheme for amplifying genetic markers for use in the invention.
- Fig. IB diagrammatically shows nucleic acid in a sample which comprises patient nucleic acid and contaminating nucleic acid, such as carry over DNA from one or more other patients.
- Fig. 1C illustrates the distribution of sequence reads among alleles of a genetic marker at a particular locus.
- Fig. ID diagrammatically shows DNA in a sample which comprises non-recombined nucleic acid (P n ) and recombined nucleic acid (P r ) of a patient and non-recombined nucleic acid (C n ) and recombined nucleic acid (C r ) contaminating nucleic acids, such as carry over DNA from one or more other patients.
- FIGS. 1E-1F illustrate an example of labeling by sampling to attach unique sequence tags to nucleic acid molecules.
- FIG. 1G illustrates an IgH transcript and sources of natural variability within it.
- Figs. 2A-2C show a two-staged PCR scheme for amplifying TCRP genes.
- FIG. 3A illustrates details of determining a nucleotide sequence of the PCR product of Fig. 2C.
- Fig. 3B illustrates details of another embodiment of determining a nucleotide sequence of the PCR product of Fig. 2C.
- 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 structure of V and J regions to improve base calls in the NDN region.
- FIG. 5 illustrates one implementation of the step of determining the sequence of a clonotype from sequence reads associated with the same sequence tag.
- 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
- the invention is directed to methods for detecting and quantifying nucleic acid contamination in a tissue sample of an individual that is used for generating a clonotype profile.
- the nucleic acid contamination of interest is that originating from other individuals, such as, from patient-to-patient or assay operator-to-patient, carry over contamination, where DNA containing clonotypes from one or more other individuals is inappropriately mixed with DNA of the individual for whom a measurement is intended.
- carry over contamination may lead to spurious estimates of the presence or quantity of clonotypes of great medical interest, such as a clonotype associated with a cancer, such as a leukemia.
- nucleic acid contamination from other individuals may be determined by measuring the amount or proportion of DNA having a genetic identity different from that of the individual of interest, i.e. the individual for whom the measurement or assay is intended.
- the genetic identity of nucleic acids in a tissue sample may be determined by conventional genetic identification assay that are based on short tandem repeats, single nucleotide polymorphisms, and the like, such as disclosed in the following references, which are incorporated by reference: Caskey et al, U.S. patent 5,364,759; Weber U.S. patent 5,075,217; Shumaker et al, Human Mutation, 7: 346-354 (1996); Sobrino et al, Forensic Sci.
- genetic markers in genomic DNA will have a polymorphic region (170) which may contain a SNP or STR along with upstream and downstream flanking regions (172) and (174), respectively, that may be used to locate the genetic marker (170) on the genome and provide primer binding sites for its amplification and analysis.
- upstream and downstream primers (176) and (180), respectively may be used to amplify genetic marker (170) in a two-stage PCR along with recombined nucleic acids from immune cells (as illustrated in Figs.
- primers (176) and (180) have tails (178) and (182), respectively, that have sequences that permit the same primers to be used in the amplification of genetic marker (170) as those used to amplify recombined nucleic acids in a second stage of a two-stage PCR.
- the tails may also contain barcodes or tags that identify the locus of the genetic marker (e.g. to identify it as being from locus 1 , 2, 3, or so on, in reference to Fig. 1A).
- primers (176) and (180) are selected so that their annealing and melting temperatures are approximately the same as those used to amplify the recombined nucleic acids (e.g. (212) of Fig. 2A) and so that the resulting amplicon has approximately the same length and GC content as the amplicons of the recombined nucleic acid molecules.
- nucleic acids (100) extracted from a tissue sample are represented as vertical lines (101) denoting genomes or fragments of genomes both from an individual of interest (102), such as a patient, and from other individuals, that is, contaminating nucleic acids (104). All of the individuals will have genetic markers at sites 1 , 2, 3, and 4, which have values Si, S 2 , S3, and S 4 (108) for the individual of interest, and values Si ', S 2 ', S3 ', and S 4 ' (1 10) for the contaminating DNA, which of course may include more than one value at each loci if the contamination is from more than one other individual.
- the values for the individual of interest may be distinguished from those of the other individuals that may be present in the tissue sample. If the genetic identity of the individual of interest based on these markers is not known beforehand, it may be determined from information generated in the assay. Assuming that DNA of the individual of interest will be the major fraction of nucleic acids, e.g. 80 percent, 90 percent, 95 percent, or the like, then the distribution of sequence reads corresponding to the different alleles of a genetic marker will be skewed to the allele corresponding to the individual of interest, i.e. the majority allele, as illustrated in Fig. 1C for genetic marker locus 1.
- Si genetic marker
- Sic genetic marker
- the allele having the highest number of sequence reads (152), Sic will correspond to the allele of the individual of interest.
- Some of the sequence reads of Sic (150) may be due to other individuals and some may be distributed among other alleles, e.g. (154). From independent counts of sequence reads from other loci, a system of linear equations may be set up and solved to give an estimate of the proportion of nucleic acids from the individual of interest and other individuals.
- Distributions of sequence reads among different alleles for each genetic marker can be used to identify the alleles which are characteristic of the individual of interest.
- a measure of the maximum sensitivity of a clonotype profile measurement e.g. for detecting a correlating clonotype, such as a leukemia clone
- line (155) in Fig. 1C is provided by the value of the largest non-patient allele, as illustrated by line (155) in Fig. 1C.
- any clonotype having a value greater than that of SID will be an accurate measurement of a patient clonotype.
- any clonotype present in a clonotype profile at a level less than that of the SID allele there exists a possibility that the clonotype does not originate from the patient.
- the invention provides a method for determining a limit of detection of a correlating clonotype in a clonotype profile of a tissue sample comprising T cells and/or B cells of an individual, wherein such method comprises steps: (a) obtaining a tissue sample from an individual, the tissue sample comprising nucleic acids from the individual and possibly from one or more other individuals, the nucleic acids including recombined nucleic acids from T cells and/or B cells and non-recombined nucleic acids; (b) generating a clonotype profile from nucleic acid of the tissue sample; (c) sequencing genetic markers at one or more genetic loci of the nucleic acids from the tissue sample to obtain proportions of the nucleic acid from the one or more other individuals and from the individual based on the genetic identification of the individual and the one or more other individuals from the genetic markers, wherein the genetic markers occur at loci present in both the recombined nucleic acids and the non-recombined nucle
- the ratio of cell types (and therefore the ratio of recombined DNA to non-recombined DNA) in a tissue sample may differ significantly between a sample from a patient, e.g. who is undergoing therapy, such as, cancer chemotherapy, and contamination from another individual, such as a healthy individual.
- further measurements can determine the different fractions of nucleic acid that are derived from different cells from a tissue sample that may include contaminating cells or nucleic acids.
- such fractions of nucleic acid may be referred to as patient (or individual of interest) non-recombined (P n ), patient recombined (P r ), contaminant non-recombined (C n ) and contaminant recombined (C r ).
- P n patient non-recombined
- P r patient recombined
- C n contaminant non-recombined
- C r contaminant recombined
- C r contaminant recombined
- An exemplary proportion of patient nucleic acid to contaminating nucleic acid may be approximately 90: 10.
- Each band is further divided into recombined nucleic acid (P r (124) and C r (128)) and non-recombined nucleic acid (P n (126) and C n (130)).
- P r /P n may be very different from the ratio C r /C n .
- the amount or concentration of total nucleic acid in a tissue sample may be determined by comparing the amount of a nucleic acid segment, such as a "house-keeping gene", e.g.
- GAPDH common to both recombined and non-recombined nucleic acids (referred to herein as a "common segment"), with one or more internal standards for such segment, or gene.
- the one or more internal standards are added in known amounts or concentrations to nucleic acid from a tissue sample prior to processing, for example, in a one- or two-stage PCR.
- total amount or concentration of recombined nucleic acid in a tissue sample may be determined by comparing the amount of sequence reads for clonotypes to the amount of sequence reads for one or more recombined sequence internal standards, e.g. V(D)J internal standards, added to nucleic acids from a tissue sample prior to processing.
- a set of recombined sequence internal standards are employed that are representative of the immune repertoire being analyzed, particularly in terms of length and composition.
- Proportions or amounts of non-recombined nucleic acid may be determined by measuring segments of nucleic acid that ordinarily would be excised during the somatic recombination process, i.e. V(D)J recombination, in which TCRs and/or BCRs are assembled. Such segments, referred to herein as “excisable segments,” are illustrated in Fig. ID by boxes (132) and (134) for the patient and contaminating nucleic acids, respectively. One or more of such segments are excised from the recombined nucleic acids as illustrated by lines (133) and (135), indicating the absence of boxes (132) and (134).
- V(D)J recombination for both TCRs and BCRs, at least two segments are always excised and lost: a first excised segment between the stretch of genome containing V region genes and the stretch of genome containing D region genes, and a second excised segment between the stretch of genome containing D region genes and the stretch of genome containing J region genes.
- the sequences of these segments are readily obtained from publicly available genome atlases, such as the NCBI Genome Browser
- V, D, and J-genes that encode human BCRs are found in or near the genomic region 14q32.33; a first excised segment between J genes and D genes has sequence located within coordinates -106,335K to -106,345K on chromosome 14 of the NCBI atlas of the human genome and a second excised segment between D genes and V genes has sequence located between the coordinates -106,406K to - 106,385K on chromosome 14 of the NCBI atlas of the human genome.
- Similar excisable segments exist for V, D, and J-genes that encode human TCRP chains, which are located on chromosome 7 at 7q32-35, e.g.
- Such an internal standard for the first and/or second excised segments is readily constructed from the above NCBI sequences using conventional techniques, such as discussed more fully below.
- the quantity, z (Nr* - Nr)/Nr*, referred to herein as "dark recombination,” provides a measure of recombined nucleic acids that are not determined in an assay.
- contaminating clonotypes may be detected and quantified by using sequence tags, as an alternative to genetic identity markers.
- each clonotype in a patient sample is determined by, and labeled with, a sequence tag, as described more fully below.
- sequence tags to determine whether a clonotype originated in the current sample or from another sample. This is accomplished by maintaining a record of sequence tags determined from each patient sample, then whenever a subsequent measurement is made the sequence tags of the current measurement are compared to those of prior measurements.
- Such records of sequence tags associated with clonotypes are conveniently maintained as electronic records on mass storage devices because of the large number of tag from each measurement and the ease of searching and comparing electronic records using conventional algorithms. If a match is found then the most likely cause is carry over contamination, provided that the populations of sequence tags employed in the measurements are sufficiently large.
- the same exemplary ratios of the size of sequence tag population to a clonotype population for labeling by sampling discussed above are applicable for detecting carry over contamination. In one embodiment, such ratio is 100: 1 or greater; in another embodiment, such ratio is 1000: 1 or greater.
- the invention comprises steps for obtaining and analyzing sequence data from a repertoire of immune molecules, such as T cell receptors (TCRs) or B cell receptors (BCRs) or defined fragments thereof, 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, from a DNA sequencer used to analyze the immune molecules.
- 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 extraction steps, sequencing chemistry, amplification chemistry, or the like.
- An aspect of the invention includes attaching a unique sequence tag to each clonotype in a sample to assist in determining whether sequence reads of such conjugates are derived from the same original clonotype.
- sequence tags are attached to the somatically recombined nucleic acid molecules to form tag- molecule conjugates wherein each recombined nucleic acid of such a conjugate has a unique sequence tag.
- attachment is made after nucleic acid molecules are extracted from a sample containing T cells and/or B cells.
- such unique sequence tags differ as greatly as possible from one another as determined by conventional distance measures for sequences, such as, a Hamming distance, or the like.
- a sequence tag of a conjugate remains far closer to its ancestral tag sequence than to that of any other tag sequence of a different conjugate. For example, if 16-mer sequence tags are employed and each such tag on a set of clonotypes has a Hamming distance of at least fifty percent, or eight nucleotides, from every other sequence tag on the clonotypes, then at least eight sequencing or amplification errors would be necessary to transform one such tag into another for a mis-read of a sequence tag (and the incorrect grouping of a sequence read of a clonotype with the wrong sequence tag).
- sequence tags are selected so that after attachment to recombined nucleic acids molecules to form tag-molecule conjugates, the Hamming distance between tags of the tag-molecule conjugates is a number at least twenty- five percent of the total length of such sequence tags (that is, each sequence tag differs in sequence from every other such tag in at least 25 percent of its nucleotides); in another embodiment, the Hamming distance between such sequence tags is a number at least 50 percent of the total length of such sequence tags.
- clonotype profiles are determined by the following steps: (a) obtaining a sample from an individual comprising T-cells and/or B-cells; (b) attaching sequence tags to molecules of recombined nucleic acids of T-cell receptor genes or immunoglobulin genes of the T-cells and/or B-cells to form tag-molecule conjugates, wherein substantially every molecule of the tag-molecule conjugates has a unique sequence tag; (c) amplifying the tag-molecule conjugates; (d) sequencing the tag-molecule conjugates; and (e) aligning sequence reads of like sequence tags to determine sequence reads corresponding to the same clonotypes of the repertoire.
- sequence tags are not only unique but also are sufficiently different from one another that the likelihood of even a large number of sequencing or amplification errors transforming one sequence tag into another would be close to zero.
- amplification of the tag-molecule conjugate is necessary for most sequencing technologies; however, whenever single-molecule sequencing technologies are employed an amplification step is optional.
- Single molecule sequencing technologies include, but are not limited to, single molecule real-time (SMRT) sequencing, nanopore sequencing, or the like, e.g. U.S. patents 7,313,308; 8,153,375; 7,907,800; 7,960,116; 8,137,569; Manrao et al, Nature Biotechnology, 4(8): 2685-2693 (2012); and the like.
- sequence tags are attached to recombined nucleic acid molecules of a sample by labeling by sampling, e.g. as disclosed by Brenner et al, U.S. patent 5,846,719; Brenner et al, U.S. patent 7,537,897; Macevicz, International patent publication WO
- polynucleotides of a population to be labeled are used to sample (by attachment, linking, or the like) sequence tags of a much larger population. That is, if the population of polynucleotides has K members (including replicates of the same polynucleotide) and the population of sequence tags has N members, then N»K.
- the size of a population of sequence tags used with the invention is at least 10 times the size of the population of clonotypes in a sample; in another embodiment, the size of a population of sequence tags used with the invention is at least 100 times the size of the population of clonotypes in a sample; and in another embodiment, the size of a population of sequence tags used with the invention is at least 1000 times the size of the population of clonotypes in a sample.
- a size of sequence tag population is selected so that substantially every clonotype in a sample will have a unique sequence tag whenever such clonotypes are combined with such sequence tag population, e.g.
- substantially every clonotype means at least 90 percent of such clonotypes will have a unique sequence tag; in other embodiments, substantially every clonotype means at least 99 percent of such clonotypes will have a unique sequence tag; in other embodiments, substantially every clonotype means at least 99.9 percent of such clonotypes will have a unique sequence tag.
- the number of T cells or B cells may be up to or about 1 million cells; thus, in some embodiments of the invention employing such samples, the number of unique sequence tags employed in labeling by sampling is at least 10 8 or in other embodiments at least 10 9 .
- sequence tags may be efficiently produced by combinatorial synthesis by reacting a mixture of all four nucleotide precursors at each addition step of a synthesis reaction, e.g. as disclosed in Church, U.S. patent 5,149,625, which is incorporated by reference.
- the number of sequence tags in a set of sequence tags made by such combinatorial synthesis is 4 k .
- a set of such sequence tags with k at least 14, or k in the range of about 14 to 18, is appropriate for attaching sequence tags to a 10 6 -member population of molecules by labeling by sampling.
- Sets of sequence tags with the above structure include many sequences that may introduce difficulties or errors while implementing the methods of the invention.
- the above combinatorially synthesized set of sequence tags includes many member tags with homopolymers segments that some sequencing
- the invention includes combinatorially synthesized sequence tags having structures that are efficient for particular method steps, such as sequencing.
- sequence tag structures efficient for sequencing-by-synthesis chemistries may be made by dividing the four natural nucleotides into disjoint subsets which are used alternatively in combinatorial synthesis, thereby preventing homopolymer segments above a given length. For example, let z be either A or C and x be either G or T, to give a sequence tag structure of
- i and j which may be the same or different, are selected to limit the size of any
- sequence tags may have lengths in the range of from 12 to 36 nucleotides; and in other embodiments, such sequence tags may have lengths in the range of from 12 to 24 nucleotides.
- other pairing of nucleotides may be used, for example, z is A or T and x is G or C; or z is A or G and x is T or C.
- z' be any combination of three of the four natural nucleotides and let x' be whatever nucleotide is not a z' (for example, z' is A, C or G, and x' is T). This gives a sequence tag structure as follows:
- attachment reactions may be used to attach unique tags to substantially every clonotype in a sample.
- such attachment is accomplished by combining a sample containing recombined nucleic acid molecules (which, in turn, comprise clonotype sequences) with a population or library of sequence tags so that members of the two populations of molecules can randomly combine and become associated or linked, e.g.
- clonotype sequences comprise linear single or double stranded polynucleotides and sequence tags are carried by reagent such as amplification primers, such as PCR primers, ligation adaptors, circularizable probes, plasmids, or the like.
- amplification primers such as PCR primers, ligation adaptors, circularizable probes, plasmids, or the like.
- Fig. IE and IF illustrate an attachment reaction comprising a PCR in which a population of sequence tags (T ls T 2 , T 3 ... T j , T j+ i ... T k , T k+ i ... T n _i, T n ) is incorporated into primers (1100).
- the population of sequence tags has a much greater size than that of recombined nucleic acid molecules (1102).
- the sequence tags are attached to the recombined nucleic acid molecules by annealing the primers to the nucleic acid molecules and extending the primers with a DNA polymerase in the first cycle of a PCR.
- the figure depicts how the recombined nucleic acid molecules select, or sample, a small fraction of the total population of sequence tags by randomly annealing to the primers by way of their common primer binding regions (1104), for example, in V region (1108). Since the primers (an therefore sequence tags) combine with the recombined nucleic acid sequence molecules randomly, there is a small possibility that the same sequence tag may be attached to different nucleic acid molecules; however, if the population of sequence tags is large as taught herein, then such possibility will be negligibly small so that substantially every recombined nucleic acid molecule will have a unique sequence tag attached.
- amplicon (1112) comprises the tag-molecule conjugates from the attachment reaction.
- Such immune molecules typically form an immune repertoire which comprises a very large set of very similar polynucleotides (e.g. >1000, but more usually from 100,000 to
- Fig. 1G illustrates diagrammatically a typical transcript of an IgH molecule (1120) 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 (1122) at a relative frequency roughly as indicated by curve (1128).
- 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 , e.g. as disclosed in Faham and Willis, U.S. patent publication 2011/0207134.
- 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, CA). 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.
- sequence reads (500) pursuant to the invention each comprise a copy of a sequence tag (502) and a copy of a clonotype (504). All sequence reads having the same sequence tag are assembled so that the nucleotides of each position in the clonotype portion can be compared.
- IgH-encoding sequences differ by no more than a single base, they will receive a distinct sequence tag, so that closely related IgH-encoding nucleic acids in a sample are not compared with one another in the clonotype determination process.
- errors in the sequence tags are not significant because the sequences of the sequence tags associated with clonotypes are so far apart in sequence space that a huge number of base changes could be sustained without one sequence tag becoming close in sequence space to any other sequence tag.
- 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.0x10 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).
- 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.
- Figs. 2A-2C illustrate exemplary steps of attaching unique sequence tags to recombined nucleic acid molecules in a two stage PCR.
- Population of recombined nucleic acid molecules (250) from a sample containing T-cells or B-cells are combined in a PCR mixture with forward and reverse primers (202) and (262).
- Primers (262) each comprise three regions: target annealing region (263) (which in this illustration is V region (206)); sequence tag (264); and primer binding region (265) for the second stage of the two-stage PCR.
- primers (262) comprise a mixture of target annealing regions to account for the diversity of V region sequence. Thus, every different primer is prepared with a sequence tag region.
- sequence tag element may be attached to C region primer (202) along with a primer binding region for the second PCR stage.
- recombined nucleic acid molecules (250) comprise constant, or C, region (203), J region (210), D region (208), and V region (206), which may represent a V(D)J segment encoding a CDR3 region of either a TCR or
- first stage amplicon (266) is produced with each member polynucleotide including sequence tag (270).
- polynucleotides of amplicon (266) are reamplified with new forward and reverse primers P5 (222) and P7 (220) which add further primer binding sites (224) and (223) for cluster formation using bridge PCR in a Solexa/Illumina sequencer.
- Primer P7 also include a secondary sequence tag (221) for optional multiplexing of samples in a single sequencing run.
- the secondary PCR amplicon (280) is produced with embedded P5 and P7 sequences by which bridge PCR may be carried out.
- the invention uses methods of labeling nucleic acids, such as fragments of genomic DNA, with unique sequence tags, which may include "mosaic tags," prior to amplification and sequencing.
- sequence tags are useful for identifying amplification and sequencing errors.
- Mosaic tags minimize sequencing and amplification artifacts due to inappropriate annealing, priming, hairpin formation, or the like, that may occur with completely random sequence tags of the prior art.
- mosaic tags are sequence tags that comprise alternating constant regions and variable regions, wherein each constant region has a position in the mosaic tag and comprises a predetermined sequence of nucleotides and each variable region has a position in the mosaic tag and comprises a predetermined number of randomly selected nucleotides.
- a 22-mer mosaic tag may have the following form:
- the sequence of b's, "***bbbbb*bb***bbb**b**" is selected to minimize the likelihood of having a perfect match in a genome of the organism making up the sample.
- all constant regions with the same position have the same length and all variable regions with the same position have the same length. This allows mosaic tags to be synthesized using partial combinatorial synthesis with conventional chemistries and instruments.
- mosaic tags comprise from 10 to 100 nucleotides, or from 12 to 80 nucleotides, or from 15 to 60 nucleotides. In some embodiments, mosaic tags comprise at least eight nucleotide positions with randomly selected nucleotides; in other embodiments, whenever mosaic tags have a length of at least 15 nucleotides, they comprise at least 12 nucleotide positions with randomly selected nucleotides. In another aspect, no variable region within a mosaic tag may have a length that is greater than seven nucleotides.
- mosaic tags may be used in the following steps: (i) preparing DNA templates from nucleic acids in a sample; (ii) labeling by sampling the DNA templates to form a multiplicity tag-template conjugates, wherein substantially every DNA template of a tag- template conjugate has a unique mosaic tag comprising alternating constant regions and variable regions, each constant region having a position in the mosaic tag and a length of from 1 to 10 nucleotides of a predetermined sequence and each variable region having a position in the mosaic tag and a length of from 1 to 10 randomly selected nucleotides, such that constant regions having the same positions have the same lengths and variable region having the same positions have the same lengths; (iii) amplifying the multiplicity of tag-template conjugates; (iv) generating a plurality of sequence reads for each of the amplified tag-template conjugates; and (v) determining a nucleotide sequence of each of the nucleic acids by determining a consensus nucleotide
- mosaic tags may be used in the following steps: (a) preparing single stranded DNA templates from nucleic acids in a sample; (b) labeling by sampling the single stranded DNA templates to form tag-template conjugates, wherein substantially every single stranded DNA template of a tag-template conjugate has a unique sequence tag (that is, a mosaic tag) having a length of at least 15 nucleotides and having the following form:
- the plurality of sequence reads is at least 10 4 ; in other embodiments, the plurality of sequence reads is at least 10 5 ; in still other embodiments, the plurality of sequence reads is at least 10 6 . In some embodiments, the total length of the above sequence tag is in the range of from 15 to 80 nucleotides.
- clonotypes of a sample are determined by first grouping sequence reads based on their sequence tags. Such grouping may be accomplished by conventional sequence alignment methods. Guidance for selecting alignment methods is available in Batzoglou, Briefings in Bio informatics, 6: 6-22 (2005), which is incorporated by reference. After sequence reads are assembled in groups corresponding to unique sequence tags, then the sequences of the associated clonotypes may be analyzed to determine the sequence of the clonotype from the sample.
- Fig. 5 illustrates an exemplary alignment and method from determining the sequence of a clonotype associated with a unique sequence tag.
- eleven sequence reads (500) are aligned by way of their respective sequence tags (502) after which nucleotides at each position of the clonotype portions of the sequence reads, indicated as 1, 2, 3, 4, ... n, are compared.
- nucleotides at position 6 (506) are t, t, g, t, t, t, t, t, t, c, t; that is, nine base calls are t's, one is "g” (508) and one is "c" (510).
- the correct base call of the clonotype sequence at a position is whatever the identity of the majority base is.
- the base call is "t", because it is the nucleotide in the majority of sequence reads at that position.
- other factors may be taken into account to determine a correct base call for a clonotype sequence, such as quality scores of the base calls of the sequence reads, identities of adjacent bases, or the like.
- clonotype profile comprising the abundances or frequencies of each different clonotype of a sample may be assembled.
- T-cells include, for example, cells that express T cell receptors (TCRs).
- B-cells include, for example, cells that express B cell receptors (BCRs).
- T-cells include helper T cells (effector T cells or Th cells), cytotoxic T cells (CTLs), memory T cells, and regulatory T cells, which may be distinguished by cell surface markers.
- a sample of T cells includes at least ⁇ , ⁇ 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, Bl cells, B2 cells, marginal-zone B cells, and follicular B cells.
- B-cells can express immunoglobulins (also referred to as antibodies or B cell receptors).
- a sample of B cells includes at least 1,000 B cells; but more typically, a sample includes at least 10,000 B cells, and more typically, at least 100,000 B cells.
- a sample includes a number of B cells in the range of from 1000 to 1,000,000 B cells.
- Samples used in the methods of the invention can come from a variety of tissues, including, for example, tumor tissue, blood and blood plasma, lymph fluid, cerebrospinal fluid surrounding the brain and the spinal cord, synovial fluid surrounding bone joints, and the like.
- the sample is a blood sample.
- the blood sample can be about 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, or 5.0 mL.
- the sample can be a tumor biopsy.
- the biopsy can be from, for example, from a tumor of the brain, liver, lung, heart, colon, kidney, or bone marrow.
- a biopsy can be an open biopsy, in which general anesthesia is used.
- the biopsy can be a closed biopsy, in which a smaller cut is made than in an open biopsy.
- the biopsy can be a core or incisional biopsy, in which part of the tissue is removed.
- the biopsy can be an excisional biopsy, in which attempts to remove an entire lesion are made.
- the biopsy can be a fine needle aspiration biopsy, in which a sample of tissue or fluid is removed with a needle.
- a sample or tissue sample includes nucleic acid, for example, DNA (e.g., genomic DNA) or RNA (e.g., messenger RNA).
- 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); Swarup et al, FEBS Lett., 581 : 795-799 (2007).
- the amount of RNA or DNA from a subject that can be analyzed includes varies widely. For example, DNA or RNA of a single cell may be all that is required for a calibration test (i.e.
- a correlating clonotype for a disease For generating a clonotype profile, sufficient nucleic acid must be in a sample to obtain a useful representation of an individual's immune receptor repertoire. More particularly, for generating a clonotype profile from genomic DNA at least 1 ng of total DNA from T cells or B cells (i.e. about 300 diploid genome equivalents) is extracted from a sample; in another embodiment, at least 2 ng of total DNA (i.e. about 600 diploid genome equivalents) is extracted from a sample; and in another embodiment, at least 3 ng of total DNA (i.e. about 900 diploid genome equivalents) is extracted from a sample.
- T cells or B cells i.e. about 300 diploid genome equivalents
- at least 2 ng of total DNA i.e. about 600 diploid genome equivalents
- at least 3 ng of total DNA i.e. about 900 diploid genome equivalents
- RNA is extracted so that at least 1000 transcripts are obtained which encode distinct TCRs, BCRs, or fragments thereof.
- the amount of RNA that corresponds to this limit varies widely from sample to sample depending on the fraction of lymphocytes in a sample, developmental stage of the lymphocytes, and the like.
- RNA used in methods of the invention may be either total RNA extracted from a tissue sample or polyA RNA extracted directly from a tissue sample or from total RNA extracted from a tissue sample.
- the above nucleic acid extractions may be carried out using commercially available kits, e.g. from Invitrogen (Carlsbad, CA), Qiagen (San Diego, CA), or like vendors. Guidance for extracting RNA is found in Liedtke et al, PCR Methods and
- a sample containing 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 .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 .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.
- 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 et al, Genome Research, 13: 954-964 (2003); and the like.
- WGA whole genome amplification
- MDA multiple displacement amplification
- Blood samples are of particular interest 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).
- PBMCs peripheral blood mononuclear cells
- Blood samples may range in volume from 100 to 10 mL; in one aspect, blood sample volumes are in the range of from 100 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, CA).
- DNeasy Blood & Tissue Kit Qiagen, Valencia, CA
- 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, CA), magnetically activated cell sorting (MACS)(Miltenyi Biotec, Auburn, CA), or the like.
- FACS fluorescently activated cell sorting
- MCS magnetically activated cell sorting
- recombined nucleic acids are present in the DNA of each individual's adaptive immunity cells as well as their associated RNA transcripts, so that 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.
- 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 a-chain and ⁇ -chain.
- the TCRa 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 are 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 (IgH) 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 TCR5 molecules (such as, molecules containing only D-J regions), and inactive IgK (for example, comprising Kde-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 TCRP 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.
- 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 R A 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.
- Known amounts of one or more internal standards to cDNA or genomic DNA can be added to an assay reaction to determine absolute quantities or concentrations of cDNA or genomic DNA samples of unknown quantity. By counting the number of molecules of the internal standard and comparing it 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 are well-known, e.g. Brenner et al, U.S. patent 7,537,897, which is incorporated herein by reference).
- 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 recombined immune molecules, such as T cell receptors or portions thereof or 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. patent 5,296,351; Gorski, U.S. patent 5,837,447; Dau, U.S. patent 6,087,096; Von Dongen et al, U.S.
- a step of generating a clonotype profile include steps of (a) amplifying a portion of T cell receptor genes and/or a portion of B cell receptor genes and (b) sequencing nucleic acids of the resulting amplicon.
- the number of amplicon nucleic acids sequenced may vary from application to application. For example, a clonotype profile to determine whether a leukemia patient is still in remission will be large so that the limit of detection of any tumor clones will be very low.
- the number of amplicon nucleic acids sequenced is at least 1000; in other embodiments, the number of amplicon nucleic acids sequenced is at least 10 4 ; in other embodiments, the number of amplicon nucleic acids sequenced is at least 10 5 .
- Such a generating step may also include further steps of coalescing sequence reads into clonotypes, enumerating or tabulating clonotypes, forming frequency distributions of clonotypes, identifying related subsets of clonotypes, displaying clonotype frequency information, and the like.
- 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, CA); 50 final reaction volume; 100 ng sample DNA; 10 pmol of each primer (subject to adjustments to balance amplification as described below); dNTPs at 200 ⁇ 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:
- 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.
- multiplex amplifications are carried out so that relative amounts of sequences in a starting population are 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 are not perfectly matched with it.
- the conditions are optimized so that each template can be amplified in a similar fashion irrespective of which primer 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 described in Faham and Willis, cited above) 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 are added to the ends of the sequences of the primary amplicon so that such sites are 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 are other considerations.
- One aspect of this is representation. This plays a role when the starting 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
- 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.
- a design goal for internal standards is providing a compound with the same amplification efficiency as the target molecule being quantified.
- Such a goal may be met by selecting an internal standard that includes the following properties: (a) homologous to target molecule (including having the same length and GC content); (b) same primer binding sites as the target molecule; (c) from a well-defined source (e.g., linearized plasmid DNA, purified PCR product, synthetic DNA, or the like); (d) readily detectable (e.g. a distinct sequence easily discernible even in presence of sequencing errors); (e) introduced into reaction in highly accurate concentrations.
- a well-defined source e.g., linearized plasmid DNA, purified PCR product, synthetic DNA, or the like
- readily detectable e.g. a distinct sequence easily discernible even in presence of sequencing errors
- Internal standards used with the invention include, but are not limited to, (a) internal standards for quantifying total nucleic acid in a sample (or more precisely, total nucleic acid, or total number of genome equivalents, in the starting material of an amplification reaction);
- internal standards for total recombined nucleic acids in a sample are designed with primers having the same melting and annealing temperatures as those used to amplify clonotypes.
- primers for second-stage amplification of the internal standards are identical to those used to amplify clonotypes.
- primers of internal standards for the embodiment illustrated may be the same as primer (202) and one (or several) of primers (212), and for a second amplification stage, primers of internal standards may be the same as primer (222) and primer (220).
- 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 are 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 5x10 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, TCRP molecules, TCRy molecules, complete TCR5 molecules, and incomplete TCR5 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 bp 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 bp per read.
- 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 as described in Faham and Willis, cited above).
- a sample is obtained that provides at least 0.5-1.0x10 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).
- 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.
- 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).
- telomeres 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 are used for generating the separate templates which, in turn, are combined and used to generate multiple sequence reads along the same strand.
- 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 (470, 472, and 474) which are 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 are 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 are 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 are selected. In other aspects, one or more sequence reads that in combination (but not separately) encompass a D and /or NDN region are sufficient.
- 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 Batzoglou, Briefings in Bio informatics, 6: 6-22 (2005), which is incorporated by reference.
- a tree search algorithm is employed, e.g.
- 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. 3A), so that the bases of the reads are in a reverse complementary relationship with one another.
- an overlap region e.g. 308 in Fig. 3A
- bases within such an overlap region are determined, at least in part, from such a reverse complementary relationship.
- 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.
- 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 a "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. 3A.
- 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. 3A).
- sequence reads are generated by extending sequencing primers, e.g. (302) and (310) in Fig. 3A, 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. 3A and 3B.
- 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 TCRP 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. 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 in Fig.
- 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.
- Figs. 4A and 4D 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, 414) all comprising the same NDN region and having different lengths encompassing successively larger portions (411, 413, 415) of V region (402).
- 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.
- One or more primers e.g. 435 and 437 in Fig. 4D
- One or more primers adjacent to the NDN region 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 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.
- amplicons 420, 422, and 4244, which may be distributed as single molecules on a solid surface, where they are further amplified by bridge PCR, or like technique.
- NDN regions Base calling in NDN regions (particularly of IgH chains) can be improved by using the codon structure 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 used to identify the recombined TCR or IgH sequences that have in frame codons from the V region through the NDN region and to the J region.
- 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
- 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 are 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. 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
- 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 are 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. Kits of the Invention
- kits comprising materials and reagents for carrying out the methods of the invention.
- kits of the invention comprise primers and recombined sequence internal standards for determining total amount or concentration of recombined nucleic acids in a tissue sample.
- kits of the invention include (a) primers and recombined sequence internal standards for determining total amount or concentration of recombined nucleic acids in a tissue sample and (b) primers and total nucleic acid internal standards for determining total amount or concentration of total nucleic acids in a tissue sample.
- such total nucleic acid internal standard is selected or derived from the group consisting of constitutively and universally expressed housekeeping genes; in a further such aspect, such total nucleic acid internal standard is selected or derived from the group consisting of the following genes: GAPDH, p 2 -microglobulin, 18S ribosomal R A, and ⁇ -actin.
- kits of the invention include (a) primers and recombined sequence internal standards for determining total amount or concentration of recombined nucleic acids in a tissue sample, (b) primers and total nucleic acid internal standards for determining total amount or concentration of total nucleic acids in a tissue sample, and (c) primers and first and/or second excised segment internal standards for determining total amount or concentration of nonrecombined nucleic acids in a tissue sample.
- kits of the invention include (a) primers and recombined sequence internal standards for determining total amount or concentration of recombined nucleic acids in a tissue sample, (b) primers and total nucleic acid internal standards for determining total amount or concentration of total nucleic acids in a tissue sample, (c) primers and first and/or second excised segment internal standards for determining total amount or concentration of nonrecombined nucleic acids in a tissue sample, and (d) primers for determining the value of alleles of genetic markers at one or more genetic loci located in both recombined nucleic acids and nonrecombined nucleic acids.
- kits of the invention include (a) primers and recombined sequence internal standards for determining total amount or concentration of recombined nucleic acids in a tissue sample, (b) primers and total nucleic acid internal standards for determining total amount or concentration of total nucleic acids in a tissue sample, (c) primers for determining the value of alleles of genetic markers at one or more genetic loci present in nonrecombined nucleic acids and absent in recombined nucleic acids, and (d) primers for determining the value of alleles of genetic markers at one or more genetic loci located in both recombined nucleic acids and nonrecombined nucleic acids.
- primers of a kit for amplifying a genetic marker may include for one or more of the genetic marker at least one forward primer specific to an upstream flanking region and at least one reverse primer specific to a downstream flanking region.
- primers of a kit may comprise a set of pairs of forward primer and reverse primer, wherein each forward primer has a tail and each reverse primer has a tail.
- tail means a portion of a primer that serves as a primer binding site in a second stage of a two-stage amplification reaction.
- a tail is not complementary to sequences in the nucleic acids being amplified, in particular, to those at or near a flanking region.
- the primer binding site of all the tails of forward primers of the pairs have the same sequence and the primer binding site of all the tails of reverse primers of the pairs have the same sequence.
- the sequences of the primer binding sites of the tails of the pairs of forward and reverse primers are different.
- 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 are "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 are 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.
- An amplification reaction may be a "realtime” 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.
- Crossing means a measure of the degree to which the distribution of clonotype abundances among clo no types of a repertoire is skewed to a single or a few
- 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.
- 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 lymphocyte which encodes an immune receptor or a portion thereof. More particularly, clonotype means a recombined nucleotide sequence of a T cell or B cell which encodes a T cell receptor (TCR) or B cell receptor (BCR), or a portion thereof.
- TCR T cell receptor
- BCR B cell receptor
- clonotypes may encode all or a portion of a VDJ rearrangement of IgH, a DJ rearrangement of IgH, a VJ rearrangement of IgK, a VJ rearrangement of IgL, a VDJ rearrangement of TCR ⁇ , a DJ rearrangement of TCR ⁇ , a VJ rearrangement of TCR a, a VJ rearrangement of TCR ⁇ , a VDJ rearrangement of TCR ⁇ , a VD rearrangement of TCR ⁇ , a Kde-V rearrangement, or the like.
- Clonotypes may also encode translocation breakpoint regions involving immune receptor genes, such as Bell -IgH or Bcll- IgH.
- clonotypes have sequences that are sufficiently long to represent or reflect the diversity of the immune molecules that they are derived from; consequently, clonotypes may vary widely in length. In some embodiments, clonotypes have lengths in the range of from 25 to 400 nucleotides; in other embodiments, clonotypes have lengths in the range of from 25 to 200 nucleotides.
- “Clonotype profile” means a listing of distinct clonotypes and their relative abundances that are derived from a population of lymphocytes. Typically, the population of lymphocytes are obtained from a tissue sample.
- the term "clonotype profile” is related to, but more general than, the immunology concept of immune "repertoire” as described in references, such as the following: 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 includes a wide variety of lists and abundances of rearranged immune receptor-encoding nucleic acids, which may be derived from selected subsets of lymphocytes (e.g. tissue-infiltrating lymphocytes, immunophenotypic subsets, or the like), or which may encode portions of immune receptors that have reduced diversity as compared to full immune receptors.
- lymphocytes e.g. tissue-infiltrating lymphocytes, immunophenotypic subsets, or the like
- clonotype profiles may comprise at least 10 3 distinct clonotypes; in other embodiments, clonotype profiles may comprise at least 10 4 distinct clonotypes; in other embodiments, clonotype profiles may comprise at least 10 5 distinct clonotypes; in other embodiments, clonotype profiles may comprise at least 10 6 distinct clonotypes. In such embodiments, such clonotype profiles may further comprise abundances or relative frequencies of each of the distinct clonotypes.
- a clonotype profile is a set of distinct recombined nucleotide sequences (with their abundances) that encode T cell receptors (TCRs) or B cell receptors (BCRs), or fragments thereof, respectively, in a population of lymphocytes of an individual, wherein the nucleotide sequences of the set have a one-to-one correspondence with distinct lymphocytes or their clonal subpopulations for substantially all of the lymphocytes of the population.
- nucleic acid segments defining clonotypes are 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, 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.
- a clonotype profile comprising human TCR ⁇ chains or fragments thereof comprises a number of distinct nucleotide sequences in the range of from 0.1 x 10 6 to 1.8 x 10 6 , or in the range of from 0.5 x 10 6 to 1.5 x 10 6 , or in the range of from 0.8 x 10 6 to 1.2 x 10 6 .
- a clonotype profile comprising human IgH chains or fragments thereof comprises a number of distinct nucleotide sequences in the range of from 0.1 x 10 6 to 1.8 x 10 6 , or in the range of from 0.5 x 10 6 to 1.5 x 10 6 , or in the range of from 0.8 x 10 6 to 1.2 x 10 6 .
- a clonotype profile of the invention comprises a set of nucleotide sequences encoding substantially all segments of the V(D)J region of an IgH chain.
- a clonotype profile 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 clonotype profile 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 clonotype profile 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 clonotype profile of the invention comprises a number of distinct nucleotide sequences that is substantially equivalent to the number of lymphocytes expressing a distinct IgH chain.
- a clonotype profile of the invention comprises a number of distinct nucleotide sequences that is substantially equivalent to the number of lymphocytes expressing a distinct TCR ⁇ chain.
- substantially equivalent means that with ninety- nine percent probability a clonotype profile 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 .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 .0001 percent or greater.
- clonotype profiles are derived from samples comprising from 10 5 to 10 7 lymphocytes. Such numbers of lymphocytes may be obtained from peripheral blood samples of from 1-10 mL.
- 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.
- Contamination means the presence in a tissue sample of one individual of nucleic acid from another individual.
- contamination means the presence of nucleic acid not originating from a patient which may affect the interpretation of a clonotype profile of the patient.
- Genetic identification means a unique correspondence between an individual and a set of values (or states) of genetic markers from one or more genetic loci of the individual.
- Genetic marker means a polymorphic segment of DNA at a genetic locus, which may be used to identify an individual.
- a genetic marker may be identified by its sequence or by adjacent or flanking sequences. Typically, a genetic marker can have a plurality of sequences, or values, in different individuals of a population. Exemplary genetic markers include, but are not limited to, short tandem repeats (STRs), single nucleotide polymorphisms (SNPs), and the like.
- the polymorphic segment of DNA may be genomic DNA or it may be reverse transcribed RNA. In one embodiment, the polymorphic segment is genomic DNA.
- a genetic marker for use with the invention is identified by amplification and sequencing using
- genetic markers are amplified and sequenced together with immune molecules during the process for generating a clonotype profile.
- Internal standard means a nucleic acid sequence that is processed in the same reaction as one or more target polynucleotides in order to permit absolute or relative quantification of the target polynucleotides in a sample.
- the reaction is an amplification reaction, such as PCR.
- An internal standard may be endogenous or exogenous. That is, an internal standard may occur naturally in the sample, or it may be added to the sample prior to a reaction.
- one or more exogenous internal standard sequences may be added to a reaction mixture in predetermined concentrations to provide a calibration to which an amplified sequence may be compared to determine the quantity of its corresponding target polynucleotide in a sample.
- Endogenous internal standards also referred to herein as "reference sequences,” are sequences natural to a sample that correspond to minimally regulated genes that exhibit a constant and cell cycle-independent level of transcription, e.g. Selvey et al, Mol. Cell Probes, 15: 307-311 (2001).
- Exemplary internal standards include, but are not limited to, sequences from the following genes: GAPDH, p 2 -microglobulin, 18S ribosomal RNA, and ⁇ -actin.
- Kit refers to any delivery system for delivering materials or reagents for carrying out a method of the invention.
- delivery systems include systems that allow for the storage, transport, or delivery of reaction reagents (e.g., primers, enzymes, internal standards, etc. in the appropriate containers) and/or supporting materials (e.g., buffers, written instructions for performing the assay etc.) from one location to another.
- reaction reagents e.g., primers, enzymes, internal standards, etc. in the appropriate containers
- supporting materials e.g., buffers, written instructions for performing the assay etc.
- kits include one or more enclosures (e.g., boxes) containing the relevant reaction reagents and/or supporting materials.
- Such contents may be delivered to the intended recipient together or separately.
- a first container may contain an enzyme for use in an assay, while a second container contains primers.
- 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 are 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 cell lymphoma (MCL), diffuse large B cell lymphoma (DLBCL), T cell lymphoma, or the like, e.g. Jaffe et al, Blood, 112: 4384-4399 (2008);
- Minimal residual disease means remaining cancer cells after treatment. The term is most frequently used in connection with treatment of lymphomas and leukemias.
- “Pecent homologous,” “percent identical,” or like terms used in reference to the comparison of a reference sequence and another 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 nano liters, e.g. 200 nL, to a few hundred ⁇ , e.g. 200 ⁇ .
- 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.
- 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 primers 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 et al, Anal. Biochem., 273: 221-228 (1999)(two-co lor 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, p 2 -microglobulin, ribosomal RNA, and the like.
- Polynucleotide refers to a linear polymer of nucleotide monomers and may be DNA or RNA. Monomers making up polynucleotides are capable of specifically binding to a natural polynucleotide by way of a regular pattern of monomer-to-monomer interactions, such as Watson-Crick type of base pairing, base stacking, Hoogsteen or reverse Hoogsteen types of base pairing, or the like. Such monomers and their internucleosidic linkages may be naturally occurring or may be analogs thereof, e.g., naturally occurring or non-naturally occurring analogs.
- Non-naturally occurring analogs may include PNAs, phosphorothioate internucleosidic linkages, bases containing linking groups permitting the attachment of labels, such as fluorophores, or haptens, and the like.
- Polynucleotides may comprise the four natural nucleosides (e.g.
- they may also comprise non-natural nucleotide analogs, e.g. including modified bases, sugars, or internucleosidic linkages. It is clear to those skilled in the art that where an enzyme has specific oligonucleotide or polynucleotide substrate requirements for activity (e.g.
- oligonucleotide refers to smaller polynucleotides, for example, having 3-60 monomeric units, or in some embodiments having from 12 to 60 monomeric units.
- a polynucleotide or oligonucleotides may be represented by a sequence of letters (upper or lower case), such as "ATGCCTG,” and it will be understood that the nucleotides are in 5' ⁇ 3' order from left to right and that "A” denotes deoxyadenosine, “C” denotes deoxycytidine, “G” denotes deoxyguanosine, and “T” denotes thymidine, “I” denotes deoxyinosine, "U” denotes uridine, unless otherwise indicated or obvious from context.
- A denotes deoxyadenosine
- C denotes deoxycytidine
- G denotes deoxyguanosine
- T denotes thymidine
- I denotes deoxyinosine
- U denotes uridine, unless otherwise indicated or obvious from context.
- 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 R A polymerase.
- the sequence of nucleotides added in the extension process is determined by the sequence of the template polynucleotide.
- primers are 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:
- 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 are 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.
- 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.
- Sequence tag (or “tag”) or “barcode” means an oligonucleotide that is attached to a polynucleotide or template molecule and is used to identify and/or track the polynucleotide or template in a reaction or a series of reactions.
- a sequence tag may be attached to the 3'- or 5'- end of a polynucleotide or template or it may be inserted into the interior of such polynucleotide or template to form a linear conjugate, sometime referred to herein as a "tagged polynucleotide,” or “tagged template,” or “tag-polynucleotide conjugate,” “tag-molecule conjugate,” or the like.
- Sequence tags may vary widely in size and compositions; the following references, which are incorporated herein by reference, provide guidance for selecting sets of sequence tags appropriate for particular embodiments: Brenner, U.S. patent 5,635,400; Brenner and Macevicz, U.S. patent 7,537,897; Brenner et al, Proc. Natl. Acad. Sci., 97: 1665-1670 (2000); Church et al, European patent publication 0 303 459; Shoemaker et al, Nature Genetics, 14: 450-456 (1996); Morris et al, European patent publication 0799897A1 ; Wallace, U.S. patent 5,981,179; and the like.
- Lengths and compositions of sequence tags can vary widely, and the selection of particular lengths and/or compositions depends on several factors including, without limitation, how tags are used to generate a readout, e.g. via a hybridization reaction or via an enzymatic reaction, such as sequencing; whether they are labeled, e.g. with a fluorescent dye or the like; the number of distinguishable oligonucleotide tags required to unambiguously identify a set of
- sequence tags can each have a length within a range of from 2 to 36 nucleotides, or from 4 to 30 nucleotides, or from 8 to 20 nucleotides, or from 6 to 10 nucleotides,
- sequence tag means a tree data structure for representing nucleotide sequences.
- 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 are 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.
- "Short tandem repeats" or "STRs” or “microsatellites” refer to genetic markers comprising a short (e.g., 1-5 nucleotide), tandemly repeated sequence motif. Microsatellites may contain repeat-motif interspersions, or "cryptically simple sequence” (Tautz, D. et al. (1986) Nature 322(6080):652-656). Such repeat-motif interspersions include simple repeat-motif interspersions wherein the micro satellite contains one or more interspersed repeats with the same length as the tandemly repeated sequence motif, but a different repeat sequence (Eichler, E. E. et al. (1994) Nat. Genet. 8:88-94; Eichler, E. E.
- tandemly repeated sequence motif is TGCC
- a simple repeat-motif interspersion may appear as follows: TGCC(TCTG) 2 (TGCC) 3 , wherein the interspersed repeat "TCTG” interrupts the repeat of the TGCC tandem repeated sequence motif.
- Repeat-motif interspersions also include more complex repeat-motif interspersions wherein the repeat motif interspersion is not the same length as the tandem repeated sequence motif.
- tandem repeated sequence motif is TGCC
- a complex repeat-motif interspersion may appear as follows: ⁇ $l HIL>TGC(TGCC) 2 , wherein the tandem repeated sequence motif is interrupted by TG and TGC.
- Other more complex repeat motif interspersions include the combination of the simple repeat- motif interspersion and the complex repeat-motif interspersion in the same micro satellite. For example, such a complex sequence repeat-motif interspersion may appear as follows:
- TGCC TGCC n (TCTG) o (TGCC)3TG(TGCC)3TGC(TGCC) 2 TGCCC(TGCC) p , wherein both forms of interspersed repeats interrupt the tandemly repeated sequence motif, TGCC.
- Microsatellites with and without interspersed repeats are encompassed by the terms "microsatellites” or "STR" as used herein.
Landscapes
- Chemical & Material Sciences (AREA)
- Organic Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Biochemistry (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Micro-Organisms Or Cultivation Processes Thereof (AREA)
Abstract
Description
Claims
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP13775514.6A EP2836611A4 (en) | 2012-04-13 | 2013-04-09 | Detection and quantitation of sample contamination in immune repertoire analysis |
SG11201406538VA SG11201406538VA (en) | 2012-04-13 | 2013-04-09 | Detection and quantitation of sample contamination in immune repertoire analysis |
CN201380031219.1A CN104395481A (en) | 2012-04-13 | 2013-04-09 | Detection and quantitation of sample contamination in immune repertoire analysis |
AU2013246050A AU2013246050B2 (en) | 2012-04-13 | 2013-04-09 | Detection and quantitation of sample contamination in immune repertoire analysis |
CA2869481A CA2869481A1 (en) | 2012-04-13 | 2013-04-09 | Detection and quantitation of sample contamination in immune repertoire analysis |
JP2015505861A JP2015512655A (en) | 2012-04-13 | 2013-04-09 | Detection and quantification of sample contamination in immune repertoire analysis |
Applications Claiming Priority (12)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261624002P | 2012-04-13 | 2012-04-13 | |
US61/624,002 | 2012-04-13 | ||
US201261658317P | 2012-06-11 | 2012-06-11 | |
US61/658,317 | 2012-06-11 | ||
US201261738277P | 2012-12-17 | 2012-12-17 | |
US61/738,277 | 2012-12-17 | ||
US201361768269P | 2013-02-22 | 2013-02-22 | |
US61/768,269 | 2013-02-22 | ||
US13/835,093 US9506119B2 (en) | 2008-11-07 | 2013-03-15 | Method of sequence determination using sequence tags |
US13/835,093 | 2013-03-15 | ||
US13/834,794 | 2013-03-15 | ||
US13/834,794 US9528160B2 (en) | 2008-11-07 | 2013-03-15 | Rare clonotypes and uses thereof |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013155119A1 true WO2013155119A1 (en) | 2013-10-17 |
Family
ID=49328112
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2013/035857 WO2013155119A1 (en) | 2012-04-13 | 2013-04-09 | Detection and quantitation of sample contamination in immune repertoire analysis |
Country Status (7)
Country | Link |
---|---|
EP (1) | EP2836611A4 (en) |
JP (1) | JP2015512655A (en) |
CN (1) | CN104395481A (en) |
AU (1) | AU2013246050B2 (en) |
CA (1) | CA2869481A1 (en) |
SG (1) | SG11201406538VA (en) |
WO (1) | WO2013155119A1 (en) |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015002908A1 (en) | 2013-07-01 | 2015-01-08 | Sequenta, Inc. | Large-scale biomolecular analysis with sequence tags |
US9043160B1 (en) | 2009-11-09 | 2015-05-26 | Sequenta, Inc. | Method of determining clonotypes and clonotype profiles |
US9150905B2 (en) | 2012-05-08 | 2015-10-06 | Adaptive Biotechnologies Corporation | Compositions and method for measuring and calibrating amplification bias in multiplexed PCR reactions |
US9181591B2 (en) | 2011-10-21 | 2015-11-10 | Adaptive Biotechnologies Corporation | Quantification of adaptive immune cell genomes in a complex mixture of cells |
CN105095687A (en) * | 2015-06-26 | 2015-11-25 | 南方科技大学 | Immune repertoire data analysis method and terminal |
US9217176B2 (en) | 2008-11-07 | 2015-12-22 | Sequenta, Llc | Methods of monitoring conditions by sequence analysis |
US9365901B2 (en) | 2008-11-07 | 2016-06-14 | Adaptive Biotechnologies Corp. | Monitoring immunoglobulin heavy chain evolution in B-cell acute lymphoblastic leukemia |
US9394567B2 (en) | 2008-11-07 | 2016-07-19 | Adaptive Biotechnologies Corporation | Detection and quantification of sample contamination in immune repertoire analysis |
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 |
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 |
US10066265B2 (en) | 2014-04-01 | 2018-09-04 | Adaptive Biotechnologies Corp. | Determining antigen-specific t-cells |
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 |
EP3344803A4 (en) * | 2015-09-01 | 2019-06-12 | The Regents of The University of California | Modular polypeptide libraries and methods of making and using same |
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 |
US11008609B2 (en) | 2017-09-01 | 2021-05-18 | Life Technologies Corporation | Compositions and methods for 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 |
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 |
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 |
US20220106633A1 (en) * | 2020-06-02 | 2022-04-07 | 10X Genomics, Inc. | Spatial transcriptomics for antigen-receptors |
US11608498B2 (en) | 2020-06-02 | 2023-03-21 | 10X Genomics, Inc. | Nucleic acid library methods |
US11613773B2 (en) | 2015-04-10 | 2023-03-28 | Spatial Transcriptomics Ab | Spatially distinguished, multiplex nucleic acid analysis of biological specimens |
US11618918B2 (en) | 2013-06-25 | 2023-04-04 | Prognosys Biosciences, Inc. | Methods and systems for determining spatial patterns of biological targets in a sample |
US11624063B2 (en) | 2020-06-08 | 2023-04-11 | 10X Genomics, Inc. | Methods of determining a surgical margin and methods of use thereof |
US11788122B2 (en) | 2011-04-13 | 2023-10-17 | 10X Genomics Sweden Ab | Methods of detecting analytes |
US11965213B2 (en) | 2019-05-30 | 2024-04-23 | 10X Genomics, Inc. | Methods of detecting spatial heterogeneity of a biological sample |
USRE50065E1 (en) | 2012-10-17 | 2024-07-30 | 10X Genomics Sweden Ab | Methods and product for optimising localised or spatial detection of gene expression in a tissue sample |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SG11201803646YA (en) * | 2015-11-04 | 2018-05-30 | Atreca Inc | Combinatorial sets of nucleic acid barcodes for analysis of nucleic acids associated with single cells |
CN105740650B (en) * | 2016-03-02 | 2019-04-05 | 广西作物遗传改良生物技术重点开放实验室 | A method of quick and precisely identifying high-throughput genomic data pollution sources |
US20210214776A1 (en) * | 2018-09-26 | 2021-07-15 | The University Of North Carolina At Chapel Hill | Compounds, compositions, and methods for improving assays |
AU2019402067A1 (en) * | 2018-12-17 | 2021-05-06 | Illumina, Inc. | Methods and means for preparing a library for sequencing |
CN110970091B (en) * | 2019-12-20 | 2023-05-23 | 北京优迅医学检验实验室有限公司 | Label quality control method and device |
CN111304309A (en) * | 2020-03-06 | 2020-06-19 | 上海韦翰斯生物医药科技有限公司 | Detection method for sequencing platform tag sequence pollution |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070172873A1 (en) * | 2006-01-23 | 2007-07-26 | Sydney Brenner | Molecular counting |
US20070264653A1 (en) * | 2006-03-10 | 2007-11-15 | Kurt Berlin | Method of identifying a biological sample for methylation analysis |
US20110207134A1 (en) * | 2008-11-07 | 2011-08-25 | Sequenta, Inc. | Monitoring health and disease status using clonotype profiles |
WO2011139371A1 (en) * | 2010-05-06 | 2011-11-10 | Sequenta, Inc. | Monitoring health and disease status using clonotype profiles |
WO2012083225A2 (en) * | 2010-12-16 | 2012-06-21 | Gigagen, Inc. | System and methods for massively parallel analysis of nycleic acids in single cells |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2307177C (en) * | 1997-10-23 | 2004-06-29 | Exact Laboratories, Inc. | Methods for detecting contamination in molecular diagnostics using pcr |
EP1598429A1 (en) * | 2004-05-19 | 2005-11-23 | Amplion Ltd. | Detection of amplicon contamination during PCR exhibiting two different annealing temperatures |
GB2483810B (en) * | 2008-11-07 | 2012-09-05 | Sequenta Inc | Methods for correlating clonotypes with diseases in a population |
-
2013
- 2013-04-09 EP EP13775514.6A patent/EP2836611A4/en not_active Withdrawn
- 2013-04-09 WO PCT/US2013/035857 patent/WO2013155119A1/en active Application Filing
- 2013-04-09 CN CN201380031219.1A patent/CN104395481A/en active Pending
- 2013-04-09 AU AU2013246050A patent/AU2013246050B2/en active Active
- 2013-04-09 CA CA2869481A patent/CA2869481A1/en not_active Abandoned
- 2013-04-09 JP JP2015505861A patent/JP2015512655A/en active Pending
- 2013-04-09 SG SG11201406538VA patent/SG11201406538VA/en unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070172873A1 (en) * | 2006-01-23 | 2007-07-26 | Sydney Brenner | Molecular counting |
US20070264653A1 (en) * | 2006-03-10 | 2007-11-15 | Kurt Berlin | Method of identifying a biological sample for methylation analysis |
US20110207134A1 (en) * | 2008-11-07 | 2011-08-25 | Sequenta, Inc. | Monitoring health and disease status using clonotype profiles |
WO2011139371A1 (en) * | 2010-05-06 | 2011-11-10 | Sequenta, Inc. | Monitoring health and disease status using clonotype profiles |
WO2012083225A2 (en) * | 2010-12-16 | 2012-06-21 | Gigagen, Inc. | System and methods for massively parallel analysis of nycleic acids in single cells |
Non-Patent Citations (3)
Title |
---|
ABBOTT ET AL.: "Design and use of signature primers to detect carry-over of amplified material", JOURNAL OF VIROLOGICAL METHODS, vol. 46, 1994, pages 51 - 59, XP023893356 * |
LOGAN ET AL.: "High-throughput VDJ sequencing for quantification of minimal residual disease in chronic lymphocytic leukemia and immune reconstitution assessment", PNAS, vol. 108, no. 52, 27 December 2011 (2011-12-27), pages 21194 - 21199, XP055170921 * |
See also references of EP2836611A4 * |
Cited By (78)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10519511B2 (en) | 2008-11-07 | 2019-12-31 | Adaptive Biotechnologies Corporation | 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 |
US9365901B2 (en) | 2008-11-07 | 2016-06-14 | Adaptive Biotechnologies Corp. | Monitoring immunoglobulin heavy chain evolution in B-cell acute lymphoblastic leukemia |
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 |
US11001895B2 (en) | 2008-11-07 | 2021-05-11 | Adaptive Biotechnologies Corporation | Methods of monitoring conditions by sequence analysis |
US9217176B2 (en) | 2008-11-07 | 2015-12-22 | Sequenta, Llc | Methods of monitoring conditions by sequence analysis |
US9228232B2 (en) | 2008-11-07 | 2016-01-05 | Sequenta, LLC. | Methods of monitoring conditions by sequence analysis |
US9512487B2 (en) | 2008-11-07 | 2016-12-06 | Adaptive Biotechnologies Corp. | Monitoring health and disease status using clonotype profiles |
US9347099B2 (en) | 2008-11-07 | 2016-05-24 | Adaptive Biotechnologies Corp. | Single cell analysis by polymerase cycling assembly |
US10155992B2 (en) | 2008-11-07 | 2018-12-18 | Adaptive Biotechnologies Corp. | Monitoring health and disease status using clonotype profiles |
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 |
US9394567B2 (en) | 2008-11-07 | 2016-07-19 | Adaptive Biotechnologies Corporation | Detection and quantification of sample contamination in immune repertoire analysis |
US10266901B2 (en) | 2008-11-07 | 2019-04-23 | Adaptive Biotechnologies Corp. | Methods of monitoring conditions by sequence analysis |
US9416420B2 (en) | 2008-11-07 | 2016-08-16 | 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 |
US11214793B2 (en) | 2009-06-25 | 2022-01-04 | Fred Hutchinson Cancer Research Center | Method of measuring adaptive immunity |
US9809813B2 (en) | 2009-06-25 | 2017-11-07 | Fred Hutchinson Cancer Research Center | Method of measuring adaptive immunity |
US9043160B1 (en) | 2009-11-09 | 2015-05-26 | Sequenta, Inc. | Method of determining clonotypes and clonotype profiles |
US11795498B2 (en) | 2011-04-13 | 2023-10-24 | 10X Genomics Sweden Ab | Methods of detecting analytes |
US11788122B2 (en) | 2011-04-13 | 2023-10-17 | 10X Genomics Sweden Ab | Methods of detecting analytes |
US10385475B2 (en) | 2011-09-12 | 2019-08-20 | Adaptive Biotechnologies Corp. | Random array sequencing of low-complexity libraries |
US9181590B2 (en) | 2011-10-21 | 2015-11-10 | Adaptive Biotechnologies Corporation | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US9279159B2 (en) | 2011-10-21 | 2016-03-08 | Adaptive Biotechnologies Corporation | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US9181591B2 (en) | 2011-10-21 | 2015-11-10 | 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 |
US9150905B2 (en) | 2012-05-08 | 2015-10-06 | Adaptive Biotechnologies Corporation | Compositions and method 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 |
US10894977B2 (en) | 2012-05-08 | 2021-01-19 | Adaptive Biotechnologies Corporation | Compositions and methods for measuring and calibrating amplification bias in multiplexed PCR reactions |
US10221461B2 (en) | 2012-10-01 | 2019-03-05 | Adaptive Biotechnologies Corp. | Immunocompetence assessment by adaptive immune receptor diversity and clonality characterization |
US11180813B2 (en) | 2012-10-01 | 2021-11-23 | Adaptive Biotechnologies Corporation | Immunocompetence assessment by adaptive immune receptor diversity and clonality characterization |
USRE50065E1 (en) | 2012-10-17 | 2024-07-30 | 10X Genomics Sweden Ab | Methods and product for optimising localised or spatial detection of gene expression in a tissue sample |
US10150996B2 (en) | 2012-10-19 | 2018-12-11 | Adaptive Biotechnologies Corp. | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US11618918B2 (en) | 2013-06-25 | 2023-04-04 | Prognosys Biosciences, Inc. | Methods and systems for determining spatial patterns of biological targets in a sample |
US11753674B2 (en) | 2013-06-25 | 2023-09-12 | Prognosys Biosciences, Inc. | Methods and systems for determining spatial patterns of biological targets in a sample |
US11821024B2 (en) | 2013-06-25 | 2023-11-21 | Prognosys Biosciences, Inc. | Methods and systems for determining spatial patterns of biological targets in a sample |
CN112375816A (en) * | 2013-07-01 | 2021-02-19 | 适应生物技术公司 | Large-scale biomolecular analysis with sequence tags |
US9708657B2 (en) | 2013-07-01 | 2017-07-18 | Adaptive Biotechnologies Corp. | Method for generating clonotype profiles using sequence tags |
CN105658812A (en) * | 2013-07-01 | 2016-06-08 | 适应生物技术公司 | Large-scale biomolecular analysis with 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 |
CN105658812B (en) * | 2013-07-01 | 2020-11-20 | 适应生物技术公司 | Large-scale biomolecular analysis with sequence tags |
WO2015002908A1 (en) | 2013-07-01 | 2015-01-08 | Sequenta, Inc. | Large-scale biomolecular analysis with sequence tags |
EP3486327A1 (en) * | 2013-07-01 | 2019-05-22 | Adaptive Biotechnologies Corporation | Large-scale biomolecular analysis with sequence tags |
JP2016523547A (en) * | 2013-07-01 | 2016-08-12 | シーケンタ,エルエルシー | Large-scale biomolecule analysis using sequence tags |
EP3017066A4 (en) * | 2013-07-01 | 2017-03-15 | Adaptive Biotechnologies Corporation | Large-scale biomolecular analysis with sequence tags |
US11248253B2 (en) | 2014-03-05 | 2022-02-15 | Adaptive Biotechnologies Corporation | Methods using randomer-containing synthetic molecules |
US11261490B2 (en) | 2014-04-01 | 2022-03-01 | Adaptive Biotechnologies Corporation | Determining antigen-specific T-cells |
US10435745B2 (en) | 2014-04-01 | 2019-10-08 | Adaptive Biotechnologies Corp. | Determining antigen-specific T-cells |
US10066265B2 (en) | 2014-04-01 | 2018-09-04 | Adaptive Biotechnologies Corp. | Determining antigen-specific t-cells |
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 |
US11613773B2 (en) | 2015-04-10 | 2023-03-28 | Spatial Transcriptomics Ab | Spatially distinguished, multiplex nucleic acid analysis of biological specimens |
US11739372B2 (en) | 2015-04-10 | 2023-08-29 | Spatial Transcriptomics Ab | Spatially distinguished, multiplex nucleic acid analysis of biological specimens |
CN105095687A (en) * | 2015-06-26 | 2015-11-25 | 南方科技大学 | Immune repertoire data analysis method and terminal |
EP3344803A4 (en) * | 2015-09-01 | 2019-06-12 | The Regents of The University of California | Modular polypeptide libraries and methods of making and using same |
US11560561B2 (en) | 2015-09-01 | 2023-01-24 | The Regents Of The University Of California | Modular polypeptide libraries and methods of making and using same |
EP4282979A3 (en) * | 2015-09-01 | 2024-02-28 | The Regents of the University of California | Modular polypeptide libraries and methods of making and using same |
EP3954772A1 (en) * | 2015-09-01 | 2022-02-16 | The Regents of The University of California | Modular polypeptide libraries and methods of making and using same |
US10428325B1 (en) | 2016-09-21 | 2019-10-01 | Adaptive Biotechnologies Corporation | Identification of antigen-specific B cell receptors |
US11008609B2 (en) | 2017-09-01 | 2021-05-18 | Life Technologies Corporation | Compositions and methods for immune repertoire sequencing |
US11254980B1 (en) | 2017-11-29 | 2022-02-22 | Adaptive Biotechnologies Corporation | Methods of profiling targeted polynucleotides while mitigating sequencing depth requirements |
US11965213B2 (en) | 2019-05-30 | 2024-04-23 | 10X Genomics, Inc. | Methods of detecting spatial heterogeneity of a biological sample |
US20220106633A1 (en) * | 2020-06-02 | 2022-04-07 | 10X Genomics, Inc. | Spatial transcriptomics for antigen-receptors |
US11608498B2 (en) | 2020-06-02 | 2023-03-21 | 10X Genomics, Inc. | Nucleic acid library methods |
US11840687B2 (en) | 2020-06-02 | 2023-12-12 | 10X Genomics, Inc. | Nucleic acid library methods |
US11845979B2 (en) | 2020-06-02 | 2023-12-19 | 10X Genomics, Inc. | Spatial transcriptomics for antigen-receptors |
US11859178B2 (en) | 2020-06-02 | 2024-01-02 | 10X Genomics, Inc. | Nucleic acid library methods |
US11692218B2 (en) * | 2020-06-02 | 2023-07-04 | 10X Genomics, Inc. | Spatial transcriptomics for antigen-receptors |
US11624063B2 (en) | 2020-06-08 | 2023-04-11 | 10X Genomics, Inc. | Methods of determining a surgical margin and methods of use thereof |
US11781130B2 (en) | 2020-06-08 | 2023-10-10 | 10X Genomics, Inc. | Methods of determining a surgical margin and methods of use thereof |
Also Published As
Publication number | Publication date |
---|---|
AU2013246050A1 (en) | 2014-10-23 |
CA2869481A1 (en) | 2013-10-17 |
AU2013246050B2 (en) | 2017-03-16 |
SG11201406538VA (en) | 2014-11-27 |
EP2836611A1 (en) | 2015-02-18 |
JP2015512655A (en) | 2015-04-30 |
EP2836611A4 (en) | 2015-12-30 |
CN104395481A (en) | 2015-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2013246050B2 (en) | Detection and quantitation of sample contamination in immune repertoire analysis | |
US9394567B2 (en) | Detection and quantification of sample contamination in immune repertoire analysis | |
AU2020210279B2 (en) | Large-scale biomolecular analysis with sequence tags | |
US9506119B2 (en) | Method of sequence determination using sequence tags | |
US9523129B2 (en) | Sequence analysis of complex amplicons | |
EP2959020B1 (en) | Method to select rare clonotypes | |
CA2859002C (en) | Detection and measurement of tissue-infiltrating lymphocytes | |
US8691510B2 (en) | Sequence analysis of complex amplicons | |
EP2859123A2 (en) | Method of sequence determination using sequence tags | |
US20140234835A1 (en) | Rare clonotypes and uses thereof | |
US20150167080A1 (en) | Method of determining clonotypes and clonotype profiles | |
WO2011139372A1 (en) | Sequence analysis of complex amplicons |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13775514 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2869481 Country of ref document: CA |
|
ENP | Entry into the national phase |
Ref document number: 2015505861 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2013246050 Country of ref document: AU Date of ref document: 20130409 Kind code of ref document: A |
|
REEP | Request for entry into the european phase |
Ref document number: 2013775514 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2013775514 Country of ref document: EP |