CA2976580A1 - Methods for diagnosing infectious disease and determining hla status using immune repertoire sequencing - Google Patents
Methods for diagnosing infectious disease and determining hla status using immune repertoire sequencing Download PDFInfo
- Publication number
- CA2976580A1 CA2976580A1 CA2976580A CA2976580A CA2976580A1 CA 2976580 A1 CA2976580 A1 CA 2976580A1 CA 2976580 A CA2976580 A CA 2976580A CA 2976580 A CA2976580 A CA 2976580A CA 2976580 A1 CA2976580 A1 CA 2976580A1
- Authority
- CA
- Canada
- Prior art keywords
- cmv
- sequences
- tcr
- subject
- status
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
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
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6881—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for tissue or cell typing, e.g. human leukocyte antigen [HLA] probes
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
-
- 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/70—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving virus or bacteriophage
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/158—Expression markers
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/172—Haplotypes
Abstract
Methods are provided for predicting a subject's infection status using high-throughput T cell receptor sequencing to match the subject's TCR repertoire to a known set of disease-associated T cell receptor sequences. The methods of the present invention may be used to predict the status of several infectious agents in a single sample from a subject. Methods are also provided for predicting a subject's HLA status using high-throughput immune receptor sequencing.
Description
2 PCT/US2016/019343 METHODS FOR DIAGNOSING INFECTIOUS DISEASE AND DETERMINING
HLA STATUS USING IMMUNE REPERTOIRE SEQUENCING
CROSS REFERENCE TO RELATED APPLICATIONS
[0001] This application claims the benefit of U.S. Provisional Patent Application No.
62/120,249, filed on February 24, 2015, U.S. Provisional Patent Application No. 62/157,249, filed on May 5, 2015, and U.S. Provisional Patent Application No. 62/215,630, filed on September 8, 2015. The contents of each of these applications are herein incorporated by reference in their entirety.
DESCRIPTION OF THE TEXT FILE SUBMITTED ELECTRONICALLY
[0002] The contents of the text file submitted electronically herewith are incorporated herein by reference in their entirety: A computer readable format copy of the Sequence Listing (filename: ADBS 024 02W0 SeqList ST25, date recorded: February 24, 2016, file size 31,026 bytes).
BACKGROUND OF THE INVENTION
HLA STATUS USING IMMUNE REPERTOIRE SEQUENCING
CROSS REFERENCE TO RELATED APPLICATIONS
[0001] This application claims the benefit of U.S. Provisional Patent Application No.
62/120,249, filed on February 24, 2015, U.S. Provisional Patent Application No. 62/157,249, filed on May 5, 2015, and U.S. Provisional Patent Application No. 62/215,630, filed on September 8, 2015. The contents of each of these applications are herein incorporated by reference in their entirety.
DESCRIPTION OF THE TEXT FILE SUBMITTED ELECTRONICALLY
[0002] The contents of the text file submitted electronically herewith are incorporated herein by reference in their entirety: A computer readable format copy of the Sequence Listing (filename: ADBS 024 02W0 SeqList ST25, date recorded: February 24, 2016, file size 31,026 bytes).
BACKGROUND OF THE INVENTION
[0003] The cellular adaptive immune system conveys broad protection against infection by pathogens through the development of a vast and highly diverse repertoire of T
cell receptor (TCR) genes, which encode cell-surface T cell receptors with randomized antigen specificity.
The ability of a subject's adaptive immune system to adequately address an incipient infection relies on activation of an appropriate antigen-specific T-cell receptor (TCR). TCR-antigen interaction is mediated by the cell-surface presentation of foreign peptides by pathogen-infected cells in the context of major histocompatibility complex (MEW) proteins.
Specifically, CD8+ T cells recognize antigen in the context of MHC class I
proteins. Since MEW class I proteins are encoded by the human leukocyte antigen (HLA) loci A, B, and C, which are highly polymorphic, the antigen specificity of a TCR is modulated across individuals by HLA context.
cell receptor (TCR) genes, which encode cell-surface T cell receptors with randomized antigen specificity.
The ability of a subject's adaptive immune system to adequately address an incipient infection relies on activation of an appropriate antigen-specific T-cell receptor (TCR). TCR-antigen interaction is mediated by the cell-surface presentation of foreign peptides by pathogen-infected cells in the context of major histocompatibility complex (MEW) proteins.
Specifically, CD8+ T cells recognize antigen in the context of MHC class I
proteins. Since MEW class I proteins are encoded by the human leukocyte antigen (HLA) loci A, B, and C, which are highly polymorphic, the antigen specificity of a TCR is modulated across individuals by HLA context.
[0004] When an antigen has been encountered, activated T cells proliferate by clonal expansion and reside in the memory T cell compartment for many years as a clonal population of cells (clones) with identical-by-descent rearranged TCR genes (Arstila TP, et al. A direct estimate of the human alphabeta T cell receptor diversity.
Science 286: 958-961, 1999). Protection against future exposure to disease-causing pathogens is conferred by the ability of activated T cells to form long-lasting memory responses.
Science 286: 958-961, 1999). Protection against future exposure to disease-causing pathogens is conferred by the ability of activated T cells to form long-lasting memory responses.
[0005] The majority of TCR diversity resides in the B chain of the TCR
alpha/beta heterodimer. Each T cell clone is encoded by a single TCRB allele that has been randomly rearranged from the germ-line TCRB locus to form a mature TCRB gene. Immense diversity is generated by combining noncontiguous TCRB variable (V), diversity (D), and joining (J) region gene segments, which collectively encode the CDR3 region, the primary region of the TCRB locus for determining antigen specificity. Deletion and template-independent insertion of nucleotides during rearrangement at the V13-D13 and D13-J13 junctions further add to the potential diversity of receptors that can be encoded (Cabaniols JP, et al.
Most alpha/beta T
cell receptor diversity is due to terminal deoxynucleotidyl transferase. J Exp Med 194: 1385-1390, 2001). Typically, at a given point in time, a healthy adult expresses approximately 10 million unique TCRB chains on their 1012 circulating T cells (Robins HS, et al. (2009) Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells. Blood 114: 4099-4107). However, observing the same TCRB chain independently in two individuals is thousands of times more common than would be expected if all rearrangements were equally likely (Robins HS, et al. (2010) Overlap and effective size of the human CD8+
T cell receptor repertoire. Science Translational Medicine 2: 47ra64). It is expected that there are many TCRI3 sequences (especially those with few or no insertions) that are present in the naïve repertoires of most individuals and that these TCRI3 sequences will reliably proliferate upon exposure to their associated antigen (V. Venturi, et al., The molecular basis for public T-cell responses? Nature reviews. Immunology 8, 231-238 (2008); published online EpubMar (10.1038/nri2260)). This over-representation of specific TCRI3 sequence rearrangements in the naïve T cell repertoire forms the basis for public T
cell responses.
alpha/beta heterodimer. Each T cell clone is encoded by a single TCRB allele that has been randomly rearranged from the germ-line TCRB locus to form a mature TCRB gene. Immense diversity is generated by combining noncontiguous TCRB variable (V), diversity (D), and joining (J) region gene segments, which collectively encode the CDR3 region, the primary region of the TCRB locus for determining antigen specificity. Deletion and template-independent insertion of nucleotides during rearrangement at the V13-D13 and D13-J13 junctions further add to the potential diversity of receptors that can be encoded (Cabaniols JP, et al.
Most alpha/beta T
cell receptor diversity is due to terminal deoxynucleotidyl transferase. J Exp Med 194: 1385-1390, 2001). Typically, at a given point in time, a healthy adult expresses approximately 10 million unique TCRB chains on their 1012 circulating T cells (Robins HS, et al. (2009) Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells. Blood 114: 4099-4107). However, observing the same TCRB chain independently in two individuals is thousands of times more common than would be expected if all rearrangements were equally likely (Robins HS, et al. (2010) Overlap and effective size of the human CD8+
T cell receptor repertoire. Science Translational Medicine 2: 47ra64). It is expected that there are many TCRI3 sequences (especially those with few or no insertions) that are present in the naïve repertoires of most individuals and that these TCRI3 sequences will reliably proliferate upon exposure to their associated antigen (V. Venturi, et al., The molecular basis for public T-cell responses? Nature reviews. Immunology 8, 231-238 (2008); published online EpubMar (10.1038/nri2260)). This over-representation of specific TCRI3 sequence rearrangements in the naïve T cell repertoire forms the basis for public T
cell responses.
[0006] Public T-cell responses occur when T cells bearing identical T-cell receptors (TCRs) are observed to dominate the response to the same antigenic epitope in multiple individuals.
Many pathogenic antigens are known to induce such a public T cell response, in which a pathogenic antigen is targeted by the same T cell receptor sequence (and found to be immunodominant) in multiple individuals with specific HLA isotypes. H. Li et al., Determinants of public T cell responses. Cell research 22, 33-42 (2012);
published online EpubJan (10.1038/cr.2012.1); H. Li, et al., Recombinatorial biases and convergent recombination determine inter-individual TCRbeta sharing in murine thymocytes.
Journal of immunology (Baltimore, Md.: 1950) 189, 2404-2413 (2012); published online EpubSep 1 (10.4049/jimmuno1.1102087). In other words, public T cell responses are observed when the space of potential high-avidity TCRf3 chains that could bind to a particular antigen-MHC
complex includes one or more TCRf3 chains that also have a high likelihood of existing in the naive repertoire at any given time. Thus, sequences associated with a public T
cell response will only be intermittently present in the naive compartment of subjects that have not been exposed to a particular antigen; however, such clones should consistently appear in the T cell repertoire of subjects who have been exposed to the antigen, having undergone clonal expansion after antigen stimulation.
Many pathogenic antigens are known to induce such a public T cell response, in which a pathogenic antigen is targeted by the same T cell receptor sequence (and found to be immunodominant) in multiple individuals with specific HLA isotypes. H. Li et al., Determinants of public T cell responses. Cell research 22, 33-42 (2012);
published online EpubJan (10.1038/cr.2012.1); H. Li, et al., Recombinatorial biases and convergent recombination determine inter-individual TCRbeta sharing in murine thymocytes.
Journal of immunology (Baltimore, Md.: 1950) 189, 2404-2413 (2012); published online EpubSep 1 (10.4049/jimmuno1.1102087). In other words, public T cell responses are observed when the space of potential high-avidity TCRf3 chains that could bind to a particular antigen-MHC
complex includes one or more TCRf3 chains that also have a high likelihood of existing in the naive repertoire at any given time. Thus, sequences associated with a public T
cell response will only be intermittently present in the naive compartment of subjects that have not been exposed to a particular antigen; however, such clones should consistently appear in the T cell repertoire of subjects who have been exposed to the antigen, having undergone clonal expansion after antigen stimulation.
[0007] Previous work on public T cell responses has identified individual examples of public T cell responses to diseases (including CMV, EBV, influenza, multiple sclerosis and other disease malignancies and autoimmune conditions), (Venturi V, et al. (2008) J
Immunol 181:
7853-7862). These studies were limited by sequencing depth and the size of investigational cohorts. Typically, these public T cell responses were studied in the context of single antigens in a single HLA context, usually using antigen-MHC tetramers to purify antigen-specific T cells. However, because of technical limitations, only a relatively small number of public T cell responses have been identified to date. Not surprisingly, these results have been limited to the high frequency, easily observable public T cell responses that dominate the immune response to their target antigens. Therefore, these public responses represent only TCRI3 sequences that are common to the naive T cell compartment of nearly all adult subjects (i.e., available for antigen response), and misses both rare T cell rearrangements which are present in only some antigen-experienced patients due to rarity among naive T
cell sequences and T cell responses that are not immunodominant and therefore never comprise a sufficient fraction of the total T cell repertoire to be reliably observed with limited sequencing depth.
Immunol 181:
7853-7862). These studies were limited by sequencing depth and the size of investigational cohorts. Typically, these public T cell responses were studied in the context of single antigens in a single HLA context, usually using antigen-MHC tetramers to purify antigen-specific T cells. However, because of technical limitations, only a relatively small number of public T cell responses have been identified to date. Not surprisingly, these results have been limited to the high frequency, easily observable public T cell responses that dominate the immune response to their target antigens. Therefore, these public responses represent only TCRI3 sequences that are common to the naive T cell compartment of nearly all adult subjects (i.e., available for antigen response), and misses both rare T cell rearrangements which are present in only some antigen-experienced patients due to rarity among naive T
cell sequences and T cell responses that are not immunodominant and therefore never comprise a sufficient fraction of the total T cell repertoire to be reliably observed with limited sequencing depth.
[0008] Further, TCR-antigen interaction is mediated by the cell-surface presentation of foreign peptides by pathogen-infected cells in the context of major histocompatibility complex (MHC) class I proteins. Since MHC class I proteins are encoded by the human leukocyte antigen (HLA) loci A, B, and C, which are highly polymorphic, the antigen specificity of a TCR is modulated across individuals by HLA context.
[0009] The binding of T cell receptors to antigens is mediated by MHC
proteins, which present antigen on the surface of cells. MHC are encoded in humans by the HLA
loci, which are highly polymorphic. This polymorphism leads to heterogeneous T cell responses to the same antigen across individuals and to differential positive and negative selection of specific T cell receptor sequences during thymic training. Determination of an individual's HLA
alleles (HLA type) has several clinical applications. One example of a clinical application for HLA typing is testing an individual's suitability as a bone marrow transplant donor.
proteins, which present antigen on the surface of cells. MHC are encoded in humans by the HLA
loci, which are highly polymorphic. This polymorphism leads to heterogeneous T cell responses to the same antigen across individuals and to differential positive and negative selection of specific T cell receptor sequences during thymic training. Determination of an individual's HLA
alleles (HLA type) has several clinical applications. One example of a clinical application for HLA typing is testing an individual's suitability as a bone marrow transplant donor.
[0010] There is a need for improved methods to diagnose and/or predict an individual's status for an infectious disease, such as CMV, EBV, HPV, small pox, and others, with increased sensitivity and accuracy. Diagnostic methods are needed that harness the information about an individual's T cell receptor sequence profile, including the presence of public T cell clones, and assess the individual's infectious disease status based on the T cell receptor sequence profile. There is also a need for inferring an individual's HLA type based on the individual's T cell receptor sequence profile. The present invention fulfills these needs and provides further related advantages.
SUMMARY OF THE INVENTION
SUMMARY OF THE INVENTION
[0011] The present disclosure is based, in part, on methods of determining or predicting the presence or absence of one or more infectious disease agents in a subject of known or unknown infection status through immunological quantification techniques combined with mathematical modeling.
[0012] The present disclosure provides a method of predicting the presence or absence of an infection in a subject of unknown infection status. In one embodiment, the genomic DNA of a sample comprising T cells obtained from the subject is subjected to amplification and high throughput sequencing to determine a T cell receptor (TCR) profile comprising unique TCR
complementarity determining region 3 (CDR3) amino acid sequences. In a particular embodiment, the TCR profile is then compared against a database of previously identified diagnostic public T cell receptor sequences that are known to be statistically significantly associated with the infection. In one embodiment, a first score for the subject is then generated by determining the proportion of unique TCR sequences in the profile of the subject that match the public TCR sequences in the database. In this embodiment, the first score is input into an algorithm that compares the first score of the subject with the known infection statuses of a plurality of subjects of known infection status. In a further embodiment, an estimated probability of infection status is then determined for the subject as the algorithm output. In some embodiments, the method further comprises an initial step of obtaining a sample comprising T cells from the subject.
complementarity determining region 3 (CDR3) amino acid sequences. In a particular embodiment, the TCR profile is then compared against a database of previously identified diagnostic public T cell receptor sequences that are known to be statistically significantly associated with the infection. In one embodiment, a first score for the subject is then generated by determining the proportion of unique TCR sequences in the profile of the subject that match the public TCR sequences in the database. In this embodiment, the first score is input into an algorithm that compares the first score of the subject with the known infection statuses of a plurality of subjects of known infection status. In a further embodiment, an estimated probability of infection status is then determined for the subject as the algorithm output. In some embodiments, the method further comprises an initial step of obtaining a sample comprising T cells from the subject.
[0013] In particular embodiments, the present disclosure additionally provides a method of predicting the presence or absence of one or more viral infections in a subject of unknown infection status. In a further embodiment, the method comprises determining a profile of unique TCR sequences from a sample obtained from the subject (e.g., a sample comprising T
cells), and inputting these unique TCR sequences into one or more algorithms.
In one embodiment, the one or more algorithms are generated by determining at least 100,000 unique TCR sequences from each of a plurality of subjects of known infection status for each of the one or more infections. In some embodiments, the method further comprises an initial step of obtaining a sample from the subject. In some embodiments, unique TCR
sequences are statistically identified that correlate with the presence or absence of each of the one or more infections, thus generating a score which is predictive of the presence or absence of each of the one or more infections. In these embodiments, the scores are then input into a logistic regression model trained on each of the plurality of subjects of known infection status for each of the one or more infections, the output of which is the prediction as to whether the subject is either positive or negative for each of the one or more infections.
cells), and inputting these unique TCR sequences into one or more algorithms.
In one embodiment, the one or more algorithms are generated by determining at least 100,000 unique TCR sequences from each of a plurality of subjects of known infection status for each of the one or more infections. In some embodiments, the method further comprises an initial step of obtaining a sample from the subject. In some embodiments, unique TCR
sequences are statistically identified that correlate with the presence or absence of each of the one or more infections, thus generating a score which is predictive of the presence or absence of each of the one or more infections. In these embodiments, the scores are then input into a logistic regression model trained on each of the plurality of subjects of known infection status for each of the one or more infections, the output of which is the prediction as to whether the subject is either positive or negative for each of the one or more infections.
[0014] In some embodiments the one or more infections are from a cytomegalovirus (CMV), an Epstein-Barr virus (EBV), a Herpes simplex virus (HSV) or a small pox virus.
[0015] In particular embodiments, the present disclosure also provides a method for predicting the presence or absence of a CMV infection in a subject of unknown infection status. In one embodiment, the genomic DNA of a sample comprising T cells obtained from a subject is subjected to amplification and high throughput sequencing to determine a TCR
profile comprising unique CDR3 amino acid sequences. In some embodiments, the method further comprises an initial step of obtaining a sample comprising T cells from the subject. In further embodiments, the TCR profile is then compared against a database of previously identified diagnostic public T cell receptor sequences that known to be statistically significantly associated with CMV infection. In some embodiments, a CMV burden score for the subject is then generated by determining the proportion of unique TCR
sequences in the profile of the subject that match the public TCR sequences in the database. In further embodiments, the calculated CMV burden score is input into a logistic regression model that compares CMV burden and CMV infection status from a plurality of subjects of known CMV
infection status. In one embodiment, an estimated probability of CMV infection status is then determined for the subject as the logistic regression model output.
profile comprising unique CDR3 amino acid sequences. In some embodiments, the method further comprises an initial step of obtaining a sample comprising T cells from the subject. In further embodiments, the TCR profile is then compared against a database of previously identified diagnostic public T cell receptor sequences that known to be statistically significantly associated with CMV infection. In some embodiments, a CMV burden score for the subject is then generated by determining the proportion of unique TCR
sequences in the profile of the subject that match the public TCR sequences in the database. In further embodiments, the calculated CMV burden score is input into a logistic regression model that compares CMV burden and CMV infection status from a plurality of subjects of known CMV
infection status. In one embodiment, an estimated probability of CMV infection status is then determined for the subject as the logistic regression model output.
[0016] In some embodiments, the database of public T cell sequences are determined to be statistically associated by obtaining unique TCR sequences from a group of subjects with CMV infection and a group of subjects without CMV infection. In one embodiment, each unique TCR sequence is subjected to a one-tailed Fisher exact test based on the presence or absence of each TCR sequence in subjects of known CMV infection status. In one embodiment, the null hypothesis being that each TCR sequence is no more common in subjects with a CMV infection than in subjects without a CMV infection. In one embodiment, a nominal p value threshold is set, and the false discovery rate (FDR) is controlled by permutation of CMV infection status in each subject to generate an empirical null distribution of p-values. In a particular embodiment, a database is then generated that comprises T cell receptor sequences that are statistically significantly shared in subjects with CMV infection.
[0017] In some embodiments, the nominal p-value threshold is less than or equal to 1.0*10'4.
[0018] In some embodiments, the statistically identified unique TCR sequences that correlate with the presence or absence of CMV infection comprise one or more of SEQ ID
NOs: 1 to 142. In some embodiments, the statistically identified unique TCR sequences that correlate with the presence or absence of CMV infection comprise SEQ ID NOs: 1 to 142.
NOs: 1 to 142. In some embodiments, the statistically identified unique TCR sequences that correlate with the presence or absence of CMV infection comprise SEQ ID NOs: 1 to 142.
[0019] In some embodiments, the step of determining the at least 100,000 unique TCRf3 sequences from the sample includes amplifying the rearranged nucleic acids encoding the TCRO CDR3 region in a multiplex PCR reaction with a mixture of forward primers specific to TCR VP gene segments and reverse primers specific to TCR 43 gene segments.
In some embodiments, the reads of the amplified nucleic acids are sequenced, and the sequence reads are processed to remove errors in the primary sequence of each read and to compress the data. In particular embodiments, a nearest neighbor algorithm is applied to collapse the data into unique sequences by merging closely related sequences to remove both PCR
and sequencing errors.
In some embodiments, the reads of the amplified nucleic acids are sequenced, and the sequence reads are processed to remove errors in the primary sequence of each read and to compress the data. In particular embodiments, a nearest neighbor algorithm is applied to collapse the data into unique sequences by merging closely related sequences to remove both PCR
and sequencing errors.
[0020] The method also includes determining an HLA association for each unique TCR
sequence in the TCR profile. In some embodiments, the unique TCR sequence is associated with an HLA-A and/or an HLA-B allele.
sequence in the TCR profile. In some embodiments, the unique TCR sequence is associated with an HLA-A and/or an HLA-B allele.
[0021] The method of the invention also comprises steps for predicting a human leukocyte antigen (HLA) allele status of a subject, comprising (a) determining an immune receptor profile of unique T-cell receptor (TCR) rearranged DNA sequences for each of a plurality of subjects, each subject having a known HLA allele status; (b) categorizing the plurality of subjects based on (i) said known HLA allele status of the subject and (ii) a presence or absence in the subject's immune receptor profile of a feature comprising a unique TCR
rearranged DNA sequence; (c) determining a statistical score for the association between a set of features and a positive HLA allele status based on (b); (d) training a machine learning model using said set of features to define a set of classifiers for each HLA
allele status; (e) inputting one or more unique TCR rearranged DNA sequences of a subject with an unknown
rearranged DNA sequence; (c) determining a statistical score for the association between a set of features and a positive HLA allele status based on (b); (d) training a machine learning model using said set of features to define a set of classifiers for each HLA
allele status; (e) inputting one or more unique TCR rearranged DNA sequences of a subject with an unknown
22 PCT/US2016/019343 HLA allele status into said machine learning model to identify one or more features that match the set of classifiers; and (f) predicting an HLA allele status of said subject based on said one or more matched features.
[0022] The method comprises determining an immune receptor profile by determining the total number of unique TCR sequences and the frequency of each unique TCR
sequence. The method also comprises determining a statistical score comprising determining a p-value using a Fisher exact two-tailed test.
[0022] The method comprises determining an immune receptor profile by determining the total number of unique TCR sequences and the frequency of each unique TCR
sequence. The method also comprises determining a statistical score comprising determining a p-value using a Fisher exact two-tailed test.
[0023] In some embodiments, the method includes determining a cutoff p-value for identifying a set of features that are significantly associated with an HLA
allele status. In another embodiment, the method includes determining a false discovery rate (FDR) of the association of a feature with an HLA allele status. In other embodiments, the method includes determining a number of false-positive associations between said feature and said HLA allele status.
allele status. In another embodiment, the method includes determining a false discovery rate (FDR) of the association of a feature with an HLA allele status. In other embodiments, the method includes determining a number of false-positive associations between said feature and said HLA allele status.
[0024] In some embodiments, the method includes training a machine learning model by training a logistic regression model using said set of identified features and said known HLA
allele statuses of each subject. In one embodiment, the method includes training a machine learning model comprises performing a leave-one out cross validation method.
In another embodiment, the method comprises performing said leave-one out cross validation method for multiple rounds. In another embodiment, the prediction is at least 80%
accurate or is at least 90% accurate.
allele statuses of each subject. In one embodiment, the method includes training a machine learning model comprises performing a leave-one out cross validation method.
In another embodiment, the method comprises performing said leave-one out cross validation method for multiple rounds. In another embodiment, the prediction is at least 80%
accurate or is at least 90% accurate.
[0025] In certain embodiments, the TCR rearranged DNA sequence is a TCRA, TCRB, TCRG or TCRD rearranged DNA sequence. In other embodiments, the HLA allele is a HLA-A2 allele or a HLA-24 allele.
BRIEF DESCRIPTION OF THE DRAWINGS
BRIEF DESCRIPTION OF THE DRAWINGS
[0026] A better understanding of the novel features of the invention and advantages of the present invention will be obtained by reference to the following description that sets forth illustrative embodiments, in which the principles of the invention are utilized, and the accompanying drawings of which:
[0027] Figure (FIG.) 1 depicts eleven distinct CDR3P VDJ recombination events, while still yielding the same nucleotide sequence, and ultimately the same CDR3P amino acid sequence (SEQ ID NO:149). The nucleotide sequences on the left side of FIG. 1 are, in order from top to bottom, SEQ ID NO:143, SEQ ID NO: 144, SEQ ID NO: 144, SEQ ID NO: 144, SEQ
ID
NO: 144, SEQ ID NO: 144, SEQ ID NO:145, SEQ ID NO:146, SEQ ID NO:146, SEQ ID
NO:147, SEQ ID NO:148.
ID
NO: 144, SEQ ID NO: 144, SEQ ID NO:145, SEQ ID NO:146, SEQ ID NO:146, SEQ ID
NO:147, SEQ ID NO:148.
[0028] FIGs. 2A, 2B and 2C provide an overview of the method for predicting disease status of an individual, according to one embodiment of the invention. FIG 2A shows a dataset of peripheral blood samples from 640 healthy subjects (287 CMV- and 353 CMV+), which were analyzed by high-throughput TCR immune receptor profiling. In FIG. 2B, unique TCRf3 sequences were identified that were present in significantly more CMV+
subjects than CMV-subjects, controlling for false determination rate (FDR) by permutation of CMV
status.
Presence of these CMV-associated TCRf3 sequences was used to build a classification model.
The top panel of FIG. 2B depicts CASSLIGVSSYNEQFF (SEQ ID NO:12). In FIG. 2C, the classification model was tested using an exhaustive leave-one-out cross-validation, in which one sample was held out of the calculations, and the process was repeated from the beginning. The resulting classification model was used to predict the CMV
status of the holdout subject.
subjects than CMV-subjects, controlling for false determination rate (FDR) by permutation of CMV
status.
Presence of these CMV-associated TCRf3 sequences was used to build a classification model.
The top panel of FIG. 2B depicts CASSLIGVSSYNEQFF (SEQ ID NO:12). In FIG. 2C, the classification model was tested using an exhaustive leave-one-out cross-validation, in which one sample was held out of the calculations, and the process was repeated from the beginning. The resulting classification model was used to predict the CMV
status of the holdout subject.
[0029] FIG. 3A depicts a graphical representation of a subject data set, where 640 subjects were specifically phenotyped for HLA type and CMB status for eligibility as an HCT donor, and segregated between 287 CMV+ donors and 353 CMV- donors. There were roughly equal numbers of seropositive and seronegative samples for the investigation of public T cell responses.
[0030] FIG. 3B shows demographic characteristics of the subjects included in this study, classified by CMV status.
[0031] FIG. 4 depicts a test for CMV-association of TCRs, where the TCR amino acid sequence CASSLIGVSSYNEQFF (SEQ ID NO:12) was selected and identified in 27 of the CMV+ subjects, and only 2 of the CMV- subjects, with a p-value of 2.8E4)8.
[0032] FIG. 5 depicts exemplary TCR amino acid sequences identified in CMV+
subjects with a p-value threshold chosen at 10-4. In order from top to bottom, these are SEQ ID NOs:
1-5. One can see that at a lower p-value, larger numbers of CMB-associated TCRf3 sequences can be identified.
subjects with a p-value threshold chosen at 10-4. In order from top to bottom, these are SEQ ID NOs:
1-5. One can see that at a lower p-value, larger numbers of CMB-associated TCRf3 sequences can be identified.
[0033] FIG. 6 depicts the CMV burden, which is represented as the proportion of unique TCRs in each subject that are significantly CMV-associated versus those TCRs that are not significantly CMV-associated.
[0034] FIG. 7 depicts the cross-validation method, which leaves out one of the initial 640 samples, then the database of CMV-specific TCRf3 sequences and associated statistics is retrained in an unsupervised fashion to eliminate bias, and the CMV serostatus of the left-out sample is classified. FIG. 7 depicts CASSLIGVSSYNEQFF (SEQ ID NO:12).
[0035] FIG. 8 depicts the results of the cross-validation method shown for all subjects and cross-validated subjects. The data is presented as the area under receiver operating characteristic (AUROC) in the y-axis versus the p-value across the x-axis. The results of the cross-validation method are shown as a plot of the false positive rate versus the true positive rate. FIG. 8 (top graph) shows data for the classification performance of all and the cross-validation (CV) datasets for each p-value threshold, measured as the area under the ROC
curve (AUROC). The number above each set of data points corresponds to the number of CMV-associated TCRf3 identified at that p-value threshold, and the rectangle indicates the dataset selected for downstream analysis (p-value = 104). FIG. 8 (bottom graph) also shows a false discovery rate (FDR) estimated for each p-value threshold used in the identification of significantly CMV-associated TCRf3 sequences, using permutations of CMV
status. The best performance is seen at a p-value of 10-4, which corresponds to an estimated FDR of ¨20%, resulting in the identification of a set of 142 TCRf3 sequences that were significantly associated with positive CMV status (listed in Table 1).
curve (AUROC). The number above each set of data points corresponds to the number of CMV-associated TCRf3 identified at that p-value threshold, and the rectangle indicates the dataset selected for downstream analysis (p-value = 104). FIG. 8 (bottom graph) also shows a false discovery rate (FDR) estimated for each p-value threshold used in the identification of significantly CMV-associated TCRf3 sequences, using permutations of CMV
status. The best performance is seen at a p-value of 10-4, which corresponds to an estimated FDR of ¨20%, resulting in the identification of a set of 142 TCRf3 sequences that were significantly associated with positive CMV status (listed in Table 1).
[0036] FIG. 9 shows the ROC curves for both the all and the cross-validation datasets.
[0037] FIGs. 10A and 10B show HLA-restriction of CMV-associated TCRf3 sequences. FIG.
10A shows the distribution of HLA-A alleles in this cohort. FIG. 10B shows the distribution of HLA-B alleles in this cohort. FIGs. 10A and 10B depict the following sequences:
CASSLAPGATNEKLFF (SEQ ID NO:8), CASSLIGVSSYNEQFF (SEQ ID NO:12), CASSPSRNTEAFF (SEQ ID NO:73), CASSLQAGANEQFF (SEQ ID NO:119), and CASASANYGYTF (SEQ ID NO:118).
10A shows the distribution of HLA-A alleles in this cohort. FIG. 10B shows the distribution of HLA-B alleles in this cohort. FIGs. 10A and 10B depict the following sequences:
CASSLAPGATNEKLFF (SEQ ID NO:8), CASSLIGVSSYNEQFF (SEQ ID NO:12), CASSPSRNTEAFF (SEQ ID NO:73), CASSLQAGANEQFF (SEQ ID NO:119), and CASASANYGYTF (SEQ ID NO:118).
[0038] FIGs. 11A and 11B show the incidence of previously reported CMV-reactive TCRf3 sequences in this cohort. FIG. 11A shows the incidence of each such TCRf3 sequence in the cohort of 640 subjects plotted along the horizontal axis by decreasing total incidence, with the incidence in CMV+ subjects above the horizontal and the incidence in CMV-subjects below the horizontal. FIG. 11B shows a histogram of incidence of these TCRf3 sequences in the cohort of 640 subjects plotted for each group of sequences.
[0039] FIG. 12 shows the concordance of TCRB sequences in the cohort as compared to those in the literature. FIG. 12 depicts the following sequences:
CASSLAPGATNEKLFF
(SEQ ID NO:8), CASSLIGVSSYNEQFF (SEQ ID NO:12), CASSPSRNTEAFF (SEQ ID
NO:73), CASSLQAGANEQFF (SEQ ID NO:119), and CASASANYGYTF (SEQ ID
NO:118).
CASSLAPGATNEKLFF
(SEQ ID NO:8), CASSLIGVSSYNEQFF (SEQ ID NO:12), CASSPSRNTEAFF (SEQ ID
NO:73), CASSLQAGANEQFF (SEQ ID NO:119), and CASASANYGYTF (SEQ ID
NO:118).
[0040] FIG. 13 depicts an overview of one embodiment of the method. FIG.13 depicts CSARDRGIGNTIYF (SEQ ID NO:152).
[0041] FIG. 14 depicts feature selection using a two-tailed Fisher exact test to determine statistical significance of the association between a feature (unique TCRB
sequence) and HLA allele status (HLA-A2+ or HLA-A2-), according to an embodiment of the invention.
An exemplary list is shown of features (unique TCRB sequences) and the number of subjects who have a particular feature and whether the subject is positive or negative for the HLA-A2 allele. FIG. 14 depicts, in order from top to bottom, SEQ ID NOs:150-161.
sequence) and HLA allele status (HLA-A2+ or HLA-A2-), according to an embodiment of the invention.
An exemplary list is shown of features (unique TCRB sequences) and the number of subjects who have a particular feature and whether the subject is positive or negative for the HLA-A2 allele. FIG. 14 depicts, in order from top to bottom, SEQ ID NOs:150-161.
[0042] FIG. 15 depicts steps for feature selection, according to an embodiment of the invention. A p value is selected as a cutoff for identifying a set of "Feature TCRs" from the entire list of possible TCR sequences. Defining a p-value threshold and permuting the allele status across individuals provides an estimate of false discovery rate. This is performed for each HLA allele, resulting in a set of allele-associated TCRf3 sequences for each HLA allele.
A p value cutoff of p<10-4 and an FDR of 0.1 was used to identify 288 TCRB
sequences that are positively associated with HLA-A2. For each of the allele-associated TCRB
sequences, the frequency of the sequence is also determined in each subject. FIG. 15 depicts, in order from left to right, SEQ ID NOs:162 and 163.
A p value cutoff of p<10-4 and an FDR of 0.1 was used to identify 288 TCRB
sequences that are positively associated with HLA-A2. For each of the allele-associated TCRB
sequences, the frequency of the sequence is also determined in each subject. FIG. 15 depicts, in order from left to right, SEQ ID NOs:162 and 163.
[0043] FIG. 16 depicts a machine learning process for fitting a logistic regression model, according to an embodiment of the invention.
[0044] FIG. 17 depicts an exhaustive leave-one-out cross validation, according to an embodiment of the invention.
[0045] FIG. 18 shows the results of the cross-validation experiment and illustrates the accuracy of the method.
DETAILED DESCRIPTION OF THE INVENTION
DETAILED DESCRIPTION OF THE INVENTION
[0046] The invention provides methods for predicting and diagnosing infectious disease in a subject that are disease-specific (specific for each particular disease state), have a universal platform (do not require different processes/reagents for each disease state), and multiplexed (are able to assay multiple disease states simultaneously.
[0047] In one embodiment, the method includes steps for high-throughput immunosequencing of rearranged TCR genes in healthy subjects with known CMV
status.
Using the immunosequencing results, the method includes searching for TCRf3 sequences that are present in multiple subjects, and identifying a set of TCRf3 sequences that are significantly associated with positive CMV status. The method also includes calculating a p-value for the association of each TCRf3 sequence with CMV status using a Fisher exact test, controlling the false discovery rate (FDR) by permutation of the CMV status, and identifying a list of CMV-associated TCRf3 sequences (for a certain FDR and p-value). A
CMV score is calculated for each subject as the proportion of all that subject's TCRf3 sequences that are represented in the catalog of CMV-associated TCRf3 sequences. The CMV score is used to distinguish between CMV+ and CMV- subjects.
status.
Using the immunosequencing results, the method includes searching for TCRf3 sequences that are present in multiple subjects, and identifying a set of TCRf3 sequences that are significantly associated with positive CMV status. The method also includes calculating a p-value for the association of each TCRf3 sequence with CMV status using a Fisher exact test, controlling the false discovery rate (FDR) by permutation of the CMV status, and identifying a list of CMV-associated TCRf3 sequences (for a certain FDR and p-value). A
CMV score is calculated for each subject as the proportion of all that subject's TCRf3 sequences that are represented in the catalog of CMV-associated TCRf3 sequences. The CMV score is used to distinguish between CMV+ and CMV- subjects.
[0048] Infectious agents include pathogens, viruses, bacteria, parasites and/or microorganisms. In some embodiments, viruses include, but are not limited to, members of the herpes virus family (such as herpes simplex viruses 1 and 2, varicella-zoster virus, EBV
(Epstein-Barr virus), human cytomegalovirus (CMV), human herpesvirus 6, human herpesvirus 7, and Kaposi's sarcoma-associated herpesvirus), Hepatitis B virus (HBV), Hepatitis C virus (HCV), human immunodeficiency viruses (HIV) I and II, influenza A virus, influenza B virus, respiratory syncytial viruses (RSV) A and B, and human metapneumovirus (MPV). Other examples include human T-cell lymphocytotrophic virus, human papillomaviruses, orthomyxo viruses, paramyxo viruses, adenoviruses, corona viruses, rhabdo viruses, polio viruses, toga viruses, bunya viruses, arena viruses, rubella viruses, reo viruses, Norovirus, human metapneumovirus (MPV), West Nile virus, Yellow fever virus, Rabies virus, Rhinovirus, Rift Valley fever virus, Marburg virus, mumps virus, measles virus, human papilloma virus (HPV), Ebola virus, Colorado tick fever virus (CTFV), and/or rhinoviruses.
(Epstein-Barr virus), human cytomegalovirus (CMV), human herpesvirus 6, human herpesvirus 7, and Kaposi's sarcoma-associated herpesvirus), Hepatitis B virus (HBV), Hepatitis C virus (HCV), human immunodeficiency viruses (HIV) I and II, influenza A virus, influenza B virus, respiratory syncytial viruses (RSV) A and B, and human metapneumovirus (MPV). Other examples include human T-cell lymphocytotrophic virus, human papillomaviruses, orthomyxo viruses, paramyxo viruses, adenoviruses, corona viruses, rhabdo viruses, polio viruses, toga viruses, bunya viruses, arena viruses, rubella viruses, reo viruses, Norovirus, human metapneumovirus (MPV), West Nile virus, Yellow fever virus, Rabies virus, Rhinovirus, Rift Valley fever virus, Marburg virus, mumps virus, measles virus, human papilloma virus (HPV), Ebola virus, Colorado tick fever virus (CTFV), and/or rhinoviruses.
[0049] Other infectious organisms include Escherichia coli, Salmonella, Shigella, Campylobacter, Klebsiella, Pseudomonas, Listeria monocytogenes, Mycobacterium tuberculosis, Mycobacterium avium-intracellulare, Yersinia, Francisella, Pasteurella, Brucella, Clostridia, Bordetella pertussis, Bacteroides, Staphylococcus aureus, Streptococcus pneumonia, B-Hemolytic strep., Corynebacteria, Legionella, Mycoplasma, Ureaplasma, Chlamydia, Clostridium difficile, Gardnerella, Trichomonas vaginalis, Nei sseria gonorrhea, Neisseria meningitides, Hemophilus influenza, Enterococcus faecalis, Proteus vulgaris, Proteus mirabilis, Helicobacter pylori, Treponema palladium, Borrelia burgdorferi, Borrelia recurrentis, Rickettsial pathogens, Nocardia, Acitnomycetes and/or Acinetobacter.
[0050] In still other embodiments, fungal infectious agents include, but are not limited to, Cryptococcus neoformans, Blastomyces dermatitidis, Histoplasma capsulatum, Coccidioides immitis, Paracoccicioides brasiliensis, Candida albicans, Aspergillus fumigautus, Phycomycetes (Rhizopus), Sporothrix schenckii, Chromomycosis, and/or Maduromycosis.
[0051] In more embodiments, parasitic agents include, but are not limited to, Plasmodium falciparum, Plasmodium malaria, Plasmodium vivax, Plasmodium ovale, Onchoverva volvulus, Leishmania, Trypanosoma spp., Schistosoma spp., Entamoeba histolytica, Cryptosporidum, Giardia spp., Trichimonas spp., Balatidium coli, Wuchereria bancrofti, Toxoplasma spp., Enterobius vermicularis, Ascaris lumbricoides, Trichuris trichiura, Dracunculus medinesis, trematodes, Diphyllobothrium latum, Taenia spp., Pneumocystis carinii, and/or Necator americanis.
[0052] The major histocompatibility complex (MHC) is a set of cell surface molecules encoded by a large gene family which controls a major part of the immune system in all vertebrates. The major function of major histocompatibility complexes is to bind to peptide fragments derived from pathogens and display them on the cell surface for recognition by the appropriate T-cells.
[0053] The major histocompatibility complex (MHC) contains two types of polymorphic MHC genes, the class I and class II genes, which encode two groups of structurally distinct but homologous proteins, and other nonpolymorphic genes whose products are involved in antigen presentation.
[0054] The human MHC is called human leukocyte antigen (HLA). HLA proteins are encoded by genes of the MHC. HLA class I antigens include HLA-A, HLA-B, and HLA-C.
HLA class II antigens include HLA-DR, HLA-DQ, HLA-DP, HLA-DM, HLA-DOA, HLA-DOB. HLA's corresponding to MHC class I (A, B, and C) present peptides from inside the cell. HLA' s corresponding to MHC class II (DP, DM, DOA, DOB, DQ, and DR) present antigens from outside of the cell to T-lymphocytes. MHC molecules mediate the binding of a given T cell receptor to a given antigen, so MHC polymorphism across individuals modulates the T cell response to a given antigen. Clinical applications of HLA typing include vaccine trials, disease associations, adverse drug reactions, platelet transfusion, and transplantation of organs and stem cells.
HLA class II antigens include HLA-DR, HLA-DQ, HLA-DP, HLA-DM, HLA-DOA, HLA-DOB. HLA's corresponding to MHC class I (A, B, and C) present peptides from inside the cell. HLA' s corresponding to MHC class II (DP, DM, DOA, DOB, DQ, and DR) present antigens from outside of the cell to T-lymphocytes. MHC molecules mediate the binding of a given T cell receptor to a given antigen, so MHC polymorphism across individuals modulates the T cell response to a given antigen. Clinical applications of HLA typing include vaccine trials, disease associations, adverse drug reactions, platelet transfusion, and transplantation of organs and stem cells.
[0055] There are numerous alleles at each HLA gene locus (Al, A2, A3, etc.).
Each person inherits one complete set of HLA alleles (haplotype) from each parent, and this combination of encoded proteins constitutes a person's HLA type (e.g., different antigens of A23, A31, B7, B44, C7, C8, DR4, DR7, DQ2, DQ7, DP2, DP3). There are more than 50,000 different known HLA types.
Each person inherits one complete set of HLA alleles (haplotype) from each parent, and this combination of encoded proteins constitutes a person's HLA type (e.g., different antigens of A23, A31, B7, B44, C7, C8, DR4, DR7, DQ2, DQ7, DP2, DP3). There are more than 50,000 different known HLA types.
[0056] As used herein, adaptive immune receptor (AIR) refers to an immune cell receptor, e.g., a T cell receptor (TCR) or an Immunoglobulin (Ig) receptor found in mammalian cells.
In certain embodiments, the adaptive immune receptor is encoded by a TCRB, TCRG, TCRA, TCRD, IGH, IGK, and IGL gene or gene segment.
In certain embodiments, the adaptive immune receptor is encoded by a TCRB, TCRG, TCRA, TCRD, IGH, IGK, and IGL gene or gene segment.
[0057] The term "primer," as used herein, refers to an oligonucleotide sequence capable of acting as a point of initiation of DNA synthesis under suitable conditions.
Such conditions include those in which synthesis of a primer extension product complementary to a nucleic acid strand is induced in the presence of four different nucleoside triphosphates and an agent for extension (e.g., a DNA polymerase or reverse transcriptase) in an appropriate buffer and at a suitable temperature.
Such conditions include those in which synthesis of a primer extension product complementary to a nucleic acid strand is induced in the presence of four different nucleoside triphosphates and an agent for extension (e.g., a DNA polymerase or reverse transcriptase) in an appropriate buffer and at a suitable temperature.
[0058] As used herein, the term "gene" refers to the segment of DNA involved in producing a polypeptide chain, such as all or a portion of a TCR or Ig polypeptide (e.g., a CDR3-containing polypeptide); it includes regions preceding and following the coding region "leader and trailer" as well as intervening sequences (introns) between individual coding segments (exons), and can also include regulatory elements (e.g., promoters, enhancers, repressor binding sites and the like), and can also include recombination signal sequences (RSSs), as described herein.
[0059] The nucleic acids of the present embodiments, also referred to herein as polynucleotides, and including oligonucleotides, can be in the form of RNA or in the form of DNA, including cDNA, genomic DNA, and synthetic DNA. The DNA can be double-stranded or single-stranded, and if single stranded can be the coding strand or non-coding (anti-sense) strand. A coding sequence which encodes a TCR or an immunoglobulin or a region thereof (e.g., a V region, a D segment, a J region, a C region, etc.) for use according to the present embodiments can be identical to the coding sequence known in the art for any given TCR or immunoglobulin gene regions or polypeptide domains (e.g., V-region domains, CDR3 domains, etc.), or can be a different coding sequence, which, as a result of the redundancy or degeneracy of the genetic code, encodes the same TCR or immunoglobulin region or polypeptide.
[0060] Unless specific definitions are provided, the nomenclature utilized in connection with, and the laboratory procedures and techniques of, molecular biology, analytical chemistry, synthetic organic chemistry, and medicinal and pharmaceutical chemistry described herein are those well-known and commonly used in the art. Standard techniques can be used for recombinant technology, molecular biological, microbiological, chemical syntheses, chemical analyses, pharmaceutical preparation, formulation, and delivery, and treatment of patients.
[0061] Unless the context requires otherwise, throughout the present specification and claims, the word "comprise" and variations thereof, such as, "comprises" and "comprising"
are to be construed in an open, inclusive sense, that is, as "including, but not limited to." By "consisting of' is meant including, and typically limited to, whatever follows the phrase "consisting of" By "consisting essentially of' is meant including any elements listed after the phrase, and limited to other elements that do not interfere with or contribute to the activity or action specified in the disclosure for the listed elements. Thus, the phrase "consisting essentially of' indicates that the listed elements are required or mandatory, but that no other elements are required and can or cannot be present depending upon whether or not they affect the activity or action of the listed elements.
are to be construed in an open, inclusive sense, that is, as "including, but not limited to." By "consisting of' is meant including, and typically limited to, whatever follows the phrase "consisting of" By "consisting essentially of' is meant including any elements listed after the phrase, and limited to other elements that do not interfere with or contribute to the activity or action specified in the disclosure for the listed elements. Thus, the phrase "consisting essentially of' indicates that the listed elements are required or mandatory, but that no other elements are required and can or cannot be present depending upon whether or not they affect the activity or action of the listed elements.
[0062] In this specification and the appended claims, the singular forms "a,"
"an" and "the"
include plural references unless the content clearly dictates otherwise.
"an" and "the"
include plural references unless the content clearly dictates otherwise.
[0063] Reference throughout this specification to "one embodiment" or "an embodiment" or "an aspect" means that a particular feature, structure or characteristic described in connection with the embodiment is included in at least one embodiment of the present invention. Thus, the appearances of the phrases "in one embodiment" or "in an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
Furthermore, the particular features, structures, or characteristics can be combined in any suitable manner in one or more embodiments.
Methods of the Invention Cells
Furthermore, the particular features, structures, or characteristics can be combined in any suitable manner in one or more embodiments.
Methods of the Invention Cells
[0064] A sample containing lymphoid nucleic acid molecules (genomic DNA, cDNA
or alternatively, messenger RNA) from a subject can be obtained. The subject is a mammalian subject, such as a human.
or alternatively, messenger RNA) from a subject can be obtained. The subject is a mammalian subject, such as a human.
[0065] Lymphocytes (B cells and/or T cells) can be obtained from a biological sample, such as from a variety of tissue and biological fluid samples. These include but are not limited to bone marrow, thymus, lymph glands, lymph nodes, peripheral tissues and blood, or solid tissue samples. Any peripheral tissue can be sampled for the presence of B and T cells and is therefore contemplated for use in the methods described herein. Peripheral blood mononuclear cells (PBMC) are isolated by techniques known to those of skill in the art, e.g., by Ficoll-Hypaque density gradient separation. In certain embodiments, whole PBMCs are used for analysis.
Nucleic Acid Extraction
Nucleic Acid Extraction
[0066] Total genomic DNA can be extracted from cells by methods known to those of skill in the art. Examples include using the QIAamp DNA blood Mini Kit (QIAGEN ) or a Qiagen DNeasy Blood extraction Kit (Qiagen, Gaithersburg, MD, USA). The approximate mass of a single haploid genome is 3 pg. Preferably, at least 100,000 to 200,000 cells are used for analysis of diversity, i.e., about 0.6 to 1.2 [tg DNA from diploid T cells.
Using PBMCs as a source, the number of T cells can be estimated to be about 30% of total cells.
Alternatively, total nucleic acid can be isolated from cells, including both genomic DNA and mRNA. In other embodiments, cDNA is transcribed from mRNA and then used as templates for amplification.
Multiplex Quantitative PCR
Using PBMCs as a source, the number of T cells can be estimated to be about 30% of total cells.
Alternatively, total nucleic acid can be isolated from cells, including both genomic DNA and mRNA. In other embodiments, cDNA is transcribed from mRNA and then used as templates for amplification.
Multiplex Quantitative PCR
[0067] Multiplex quantitative PCR was performed as described herein and in Robins et at., 2009 Blood 114, 4099; Robins et at., 2010 Sci. Translat. Med. 2:47ra64; Robins et at., 2011 Immunol. Meth. doi:10.1016/j im.2011.09. 001; Sherwood et al. 2011 Sci.
Translat. Med.
3:90ra61; US 2012/0058902, US 2010/0330571, WO/2010/151416, WO/2011/106738 (PCT/US2011/026373), US 2015/0299785, W02012/027503 (PCT/US2011/049012), US
2013/0288237, US 9,181,590, US 9,181,591, and US 2013/0253842, which are incorporated by reference in their entireties. The present methods involve a multiplex PCR
method using a set of forward primers that specifically hybridize to V segments and a set of reverse primers that specifically hybridize to the J segments of a TCR or Ig locus, where a multiplex PCR
reaction using the primers allows amplification of all the possible VJ (and VDJ) combinations within a given population of T or B cells.
Translat. Med.
3:90ra61; US 2012/0058902, US 2010/0330571, WO/2010/151416, WO/2011/106738 (PCT/US2011/026373), US 2015/0299785, W02012/027503 (PCT/US2011/049012), US
2013/0288237, US 9,181,590, US 9,181,591, and US 2013/0253842, which are incorporated by reference in their entireties. The present methods involve a multiplex PCR
method using a set of forward primers that specifically hybridize to V segments and a set of reverse primers that specifically hybridize to the J segments of a TCR or Ig locus, where a multiplex PCR
reaction using the primers allows amplification of all the possible VJ (and VDJ) combinations within a given population of T or B cells.
[0068] Exemplary V segment primers and J segment primers are described in US
2012/0058902, US 2010/033057, WO/2010/151416, WO/2011/106738 (PCT/US2011/026373), US 2015/0299785, W02012/027503 (PCT/US2011/049012), US
2013/0288237, US 9,181,590, US 9,181,591, US 2013/0253842, WO 2013/188831 (PCT/U52013/045994), which are incorporated by reference in their entireties.
2012/0058902, US 2010/033057, WO/2010/151416, WO/2011/106738 (PCT/US2011/026373), US 2015/0299785, W02012/027503 (PCT/US2011/049012), US
2013/0288237, US 9,181,590, US 9,181,591, US 2013/0253842, WO 2013/188831 (PCT/U52013/045994), which are incorporated by reference in their entireties.
[0069] A multiplex PCR system can be used to amplify rearranged adaptive immune cell receptor loci. In certain embodiments, the CDR3 region is amplified from a TCRA, TCRB, TCRG or TCRD CDR3 region or similarly from an IgH or IgL (lambda or kappa) locus. A
plurality of V-segment and J-segment primers are used to amplify substantially all (e.g., greater than 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98% or 99%) rearranged adaptive immune receptor CDR3-encoding regions to produce a multiplicity of amplified rearranged DNA molecules. In certain embodiments, primers are designed so that each amplified rearranged DNA molecule is less than 600 nucleotides in length, thereby excluding amplification products from non-rearranged adaptive immune receptor loci.
plurality of V-segment and J-segment primers are used to amplify substantially all (e.g., greater than 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98% or 99%) rearranged adaptive immune receptor CDR3-encoding regions to produce a multiplicity of amplified rearranged DNA molecules. In certain embodiments, primers are designed so that each amplified rearranged DNA molecule is less than 600 nucleotides in length, thereby excluding amplification products from non-rearranged adaptive immune receptor loci.
[0070] In some embodiments, two pools of primers are used in a single, highly multiplexed PCR reaction. The "forward" pool of primers can include a plurality of V
segment oligonucleotide primers and the reverse pool can include a plurality of J
segment oligonucleotide primers. In some embodiments, there is a primer that is specific to (e.g., having a nucleotide sequence complementary to a unique sequence region of) each V region segment and to each J region segment in the respective TCR or Ig gene locus.
In other embodiments, a primer can hybridize to one or more V segments or J segments, thereby reducing the number of primers required in the multiplex PCR. In certain embodiments, the J-segment primers anneal to a conserved sequence in the joining ("J") segment.
segment oligonucleotide primers and the reverse pool can include a plurality of J
segment oligonucleotide primers. In some embodiments, there is a primer that is specific to (e.g., having a nucleotide sequence complementary to a unique sequence region of) each V region segment and to each J region segment in the respective TCR or Ig gene locus.
In other embodiments, a primer can hybridize to one or more V segments or J segments, thereby reducing the number of primers required in the multiplex PCR. In certain embodiments, the J-segment primers anneal to a conserved sequence in the joining ("J") segment.
[0071] Each primer can be designed such that a respective amplified DNA
segment is obtained that includes a sequence portion of sufficient length to identify each J segment unambiguously based on sequence differences amongst known J-region encoding gene segments in the human genome database, and also to include a sequence portion to which a J-segment-specific primer can anneal for resequencing. This design of V- and J-segment-specific primers enables direct observation of a large fraction of the somatic rearrangements present in the adaptive immune receptor gene repertoire within an individual.
segment is obtained that includes a sequence portion of sufficient length to identify each J segment unambiguously based on sequence differences amongst known J-region encoding gene segments in the human genome database, and also to include a sequence portion to which a J-segment-specific primer can anneal for resequencing. This design of V- and J-segment-specific primers enables direct observation of a large fraction of the somatic rearrangements present in the adaptive immune receptor gene repertoire within an individual.
[0072] In one embodiment, the present disclosure provides a plurality of V-segment primers and a plurality of J-segment primers. The plurality of V-segment primers and the plurality of J-segment primers amplify all or substantially all combinations of the V- and J-segments of a rearranged immune receptor locus. In some embodiments, the method provides amplification of substantially all of the rearranged AIR sequences in a lymphoid cell and is capable of quantifying the diversity of the TCR or IG repertoire of at least 106, 105, 104, or 103 unique rearranged AIR sequences in a sample. "Substantially all combinations" can refer to at least 80%, 85%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99% or more of all the combinations of the V- and J-segments of a rearranged immune receptor locus.
In certain embodiments, the plurality of V-segment primers and the plurality of J-segment primers amplify all of the combinations of the V- and J-segments of a rearranged immune receptor locus.
In certain embodiments, the plurality of V-segment primers and the plurality of J-segment primers amplify all of the combinations of the V- and J-segments of a rearranged immune receptor locus.
[0073] In general, a multiplex PCR system can use 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, or 25, and in certain embodiments, at least 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, or 39, and in other embodiments 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 65, 70, 75, 80, 85, or more forward primers, in which each forward primer specifically hybridizes to or is complementary to a sequence corresponding to one or more V region segments. The multiplex PCR
system also uses at least 2, 3, 4, 5, 6, or 7, and in certain embodiments, 8, 9, 10, 11, 12 or 13 reverse primers, or 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, or 25 or more primers, in which each reverse primer specifically hybridizes to or is complementary to a sequence corresponding to one or more J region segments. Various combinations of V and J segment primers can be used to amplify the full diversity of TCR and IG sequences in a repertoire.
For details on the multiplex PCR system, including primer oligonucleotide sequences for amplifying TCR and IG sequences, see, e.g., Robins et at., 2009 Blood 114, 4099; Robins et at., 2010 Sci.
Translat. Med. 2:47ra64; Robins et al., 20111 Immunol. Meth.
doi:10.1016/j.jim.2011.09.
001; Sherwood et al. 2011 Sci. Translat. Med. 3:90ra61; US 2012/0058902, US
2010/033057, WO/2010/151416, WO/2011/106738 (PCT/U52011/026373), US
2015/0299785, W02012/027503 (PCT/U52011/049012), US 2013/0288237, US
9,181,590, US 9,181,591, US 2013/0253842, WO 2013/188831 (PCT/U52013/045994), which is each incorporated by reference in its entirety.
system also uses at least 2, 3, 4, 5, 6, or 7, and in certain embodiments, 8, 9, 10, 11, 12 or 13 reverse primers, or 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, or 25 or more primers, in which each reverse primer specifically hybridizes to or is complementary to a sequence corresponding to one or more J region segments. Various combinations of V and J segment primers can be used to amplify the full diversity of TCR and IG sequences in a repertoire.
For details on the multiplex PCR system, including primer oligonucleotide sequences for amplifying TCR and IG sequences, see, e.g., Robins et at., 2009 Blood 114, 4099; Robins et at., 2010 Sci.
Translat. Med. 2:47ra64; Robins et al., 20111 Immunol. Meth.
doi:10.1016/j.jim.2011.09.
001; Sherwood et al. 2011 Sci. Translat. Med. 3:90ra61; US 2012/0058902, US
2010/033057, WO/2010/151416, WO/2011/106738 (PCT/U52011/026373), US
2015/0299785, W02012/027503 (PCT/U52011/049012), US 2013/0288237, US
9,181,590, US 9,181,591, US 2013/0253842, WO 2013/188831 (PCT/U52013/045994), which is each incorporated by reference in its entirety.
[0074] Oligonucleotides or polynucleotides that are capable of specifically hybridizing or annealing to a target nucleic acid sequence by nucleotide base complementarity can do so under moderate to high stringency conditions. In one embodiment, suitable moderate to high stringency conditions for specific PCR amplification of a target nucleic acid sequence can be between 25 and 80 PCR cycles, with each cycle consisting of a denaturation step (e.g., about 10-30 seconds (s) at greater than about 95 C), an annealing step (e.g., about 10-30s at about 60-68 C), and an extension step (e.g., about 10-60s at about 60-72 C), optionally according to certain embodiments with the annealing and extension steps being combined to provide a two-step PCR. As would be recognized by the skilled person, other PCR reagents can be added or changed in the PCR reaction to increase specificity of primer annealing and amplification, such as altering the magnesium concentration, optionally adding DMSO, and/or the use of blocked primers, modified nucleotides, peptide-nucleic acids, and the like.
[0075] A primer is preferably a single-stranded DNA. The appropriate length of a primer depends on the intended use of the primer but typically ranges from 6 to 50 nucleotides, or in certain embodiments, from 15-35 nucleotides in length. Short primer molecules generally require cooler temperatures to form sufficiently stable hybrid complexes with the template.
A primer need not reflect the exact sequence of the template nucleic acid, but must be sufficiently complementary to hybridize with the template. The design of suitable primers for the amplification of a given target sequence is well known in the art and described in the literature cited herein.
A primer need not reflect the exact sequence of the template nucleic acid, but must be sufficiently complementary to hybridize with the template. The design of suitable primers for the amplification of a given target sequence is well known in the art and described in the literature cited herein.
[0076] In some embodiments, the V- and J-segment primers are used to produce a plurality of amplicons from the multiplex PCR reaction. In certain embodiments, the amplicons range in size from 10, 20, 30, 40, 50, 75, 100, 200, 300, 400, 500, 600, 700, 800 or more nucleotides in length. In preferred embodiments, the amplicons have a size between 50-600 nucleotides in length.
[0077] According to non-limiting theory, these embodiments exploit current understanding in the art (also described above) that once an adaptive immune cell (e.g., a T or B lymphocyte) has rearranged its adaptive immune receptor-encoding (e.g., TCR or Ig) genes, its progeny cells possess the same adaptive immune receptor-encoding gene rearrangement, thus giving rise to a clonal population (clones) that can be uniquely identified by the presence therein of rearranged (e.g., CDR3-encoding) V- and J-gene segments that can be amplified by a specific pairwise combination of V- and J-specific oligonucleotide primers as herein disclosed.
[0078] In some embodiments, the V segment primers and J segment primers each include a second sequence at the 5'-end of the primer that is not complementary to the target V or J
segment. The second sequence can comprise an oligonucleotide having a sequence that is selected from (i) a universal adaptor oligonucleotide sequence, and (ii) a sequencing platform-specific oligonucleotide sequence that is linked to and positioned 5' to a first universal adaptor oligonucleotide sequence. Examples of universal adaptor oligonucleotide sequences can be pGEX forward and pGEX reverse adaptor sequences.
segment. The second sequence can comprise an oligonucleotide having a sequence that is selected from (i) a universal adaptor oligonucleotide sequence, and (ii) a sequencing platform-specific oligonucleotide sequence that is linked to and positioned 5' to a first universal adaptor oligonucleotide sequence. Examples of universal adaptor oligonucleotide sequences can be pGEX forward and pGEX reverse adaptor sequences.
[0079] In some embodiments, the resulting amplicons using the V-segment and J-segment primers described above include amplified V and J segments and the universal adaptor oligonucleotide sequences. The universal adaptor sequence can be complementary to an oligonucleotide sequence found in a tailing primer. Tailing primers can be used in a second PCR reaction to generate a second set of amplicons. In some embodiments, tailing primers can have the general formula:
5' -P- S¨B¨U¨ 3' (III),
5' -P- S¨B¨U¨ 3' (III),
[0080] wherein P comprises a sequencing platform-specific oligonucleotide,
[0081] wherein S comprises a sequencing platform tag-containing oligonucleotide sequence;
[0082] wherein B comprises an oligonucleotide barcode sequence and wherein said oligonucleotide barcode sequence can be used to identify a sample source, and
[0083] wherein U comprises a sequence that is complementary to the universal adaptor oligonucleotide sequence or is the same as the universal adaptor oligonucleotide sequence.
[0084] Additional description about universal adaptor oligonucleotide sequences, barcodes, and tailing primers are found in WO 2013/188831 (PCT/US13/45994), which is incorporated by reference in its entirety.
Amplification Bias Control
Amplification Bias Control
[0085] Multiplex PCR assays can result in a bias in the total numbers of amplicons produced from a sample, given that certain primer sets are more efficient in amplification than others.
To overcome the problem of such biased utilization of subpopulations of amplification primers, methods can be used that provide a template composition for standardizing the amplification efficiencies of the members of an oligonucleotide primer set, where the primer set is capable of amplifying rearranged DNA encoding a plurality of adaptive immune receptors (TCR or Ig) in a biological sample that comprises DNA from lymphoid cells.
To overcome the problem of such biased utilization of subpopulations of amplification primers, methods can be used that provide a template composition for standardizing the amplification efficiencies of the members of an oligonucleotide primer set, where the primer set is capable of amplifying rearranged DNA encoding a plurality of adaptive immune receptors (TCR or Ig) in a biological sample that comprises DNA from lymphoid cells.
[0086] In some embodiments, a template composition is used to standardize the various amplification efficiencies of the primer sets. The template composition can comprise a plurality of diverse template oligonucleotides of general formula (I):
5'-U1-B1-V-B2-R-J-B3-U2-3' (I)
5'-U1-B1-V-B2-R-J-B3-U2-3' (I)
[0087] The constituent template oligonucleotides are diverse with respect to the nucleotide sequences of the individual template oligonucleotides. The individual template oligonucleotides can vary in nucleotide sequence considerably from one another as a function of significant sequence variability among the large number of possible TCR or BCR variable (V) and joining (J) region polynucleotides. Sequences of individual template oligonucleotide species can also vary from one another as a function of sequence differences in Ul, U2, B
(B1, B2 and B3) and R oligonucleotides that are included in a particular template within the diverse plurality of templates.
(B1, B2 and B3) and R oligonucleotides that are included in a particular template within the diverse plurality of templates.
[0088] In certain embodiments, V is a polynucleotide comprising at least 20, 30, 60, 90, 120, 150, 180, or 210, and not more than 1000, 900, 800, 700, 600 or 500 contiguous nucleotides of an adaptive immune receptor variable (V) region encoding gene sequence, or the complement thereof, and in each of the plurality of template oligonucleotide sequences V
comprises a unique oligonucleotide sequence.
comprises a unique oligonucleotide sequence.
[0089] In some embodiments, J is a polynucleotide comprising at least 15-30, 31-60, 61-90, 91-120, or 120-150, and not more than 600, 500, 400, 300 or 200 contiguous nucleotides of an adaptive immune receptor joining (J) region encoding gene sequence, or the complement thereof, and in each of the plurality of template oligonucleotide sequences J
comprises a unique oligonucleotide sequence.
comprises a unique oligonucleotide sequence.
[0090] Ul and U2 can be each either nothing or each comprise an oligonucleotide having, independently, a sequence that is selected from (i) a universal adaptor oligonucleotide sequence, and (ii) a sequencing platform-specific oligonucleotide sequence that is linked to and positioned 5' to the universal adaptor oligonucleotide sequence.
[0091] Bl, B2 and B3 can be each either nothing or each comprise an oligonucleotide B that comprises a first and a second oligonucleotide barcode sequence of 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 25, 30, 35, 40, 45, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900 or 1000 contiguous nucleotides (including all integer values therebetween), wherein in each of the plurality of template oligonucleotide sequences B
comprises a unique oligonucleotide sequence in which (i) the first barcode sequence uniquely identifies the unique V oligonucleotide sequence of the template oligonucleotide and (ii) the second barcode sequence uniquely identifies the unique J oligonucleotide sequence of the template oligonucleotide.
comprises a unique oligonucleotide sequence in which (i) the first barcode sequence uniquely identifies the unique V oligonucleotide sequence of the template oligonucleotide and (ii) the second barcode sequence uniquely identifies the unique J oligonucleotide sequence of the template oligonucleotide.
[0092] R can be either nothing or comprises a restriction enzyme recognition site that comprises an oligonucleotide sequence that is absent from V, J, Ul, U2, Bl, B2 and B3.
[0093] Methods are used with the template composition for determining non-uniform nucleic acid amplification potential among members of a set of oligonucleotide amplification primers that are capable of amplifying productively rearranged DNA encoding one or a plurality of adaptive immune receptors in a biological sample that comprises DNA from lymphoid cells of a subject. The method can include the steps of: (a) amplifying DNA of a template composition for standardizing amplification efficiency of an oligonucleotide primer set in a multiplex polymerase chain reaction (PCR) that comprises: (i) the template composition (I) described above, wherein each template oligonucleotide in the plurality of template oligonucleotides is present in a substantially equimolar amount; (ii) an oligonucleotide amplification primer set that is capable of amplifying productively rearranged DNA encoding one or a plurality of adaptive immune receptors in a biological sample that comprises DNA
from lymphoid cells of a subject.
from lymphoid cells of a subject.
[0094] The primer set can include: (1) in substantially equimolar amounts, a plurality of V-segment oligonucleotide primers that are each independently capable of specifically hybridizing to at least one polynucleotide encoding an adaptive immune receptor V-region polypeptide or to the complement thereof, wherein each V-segment primer comprises a nucleotide sequence of at least 15 contiguous nucleotides that is complementary to at least one functional adaptive immune receptor V region-encoding gene segment and wherein the plurality of V-segment primers specifically hybridize to substantially all functional adaptive immune receptor V region-encoding gene segments that are present in the template composition, and (2) in substantially equimolar amounts, a plurality of J-segment oligonucleotide primers that are each independently capable of specifically hybridizing to at least one polynucleotide encoding an adaptive immune receptor J-region polypeptide or to the complement thereof, wherein each J-segment primer comprises a nucleotide sequence of at least 15 contiguous nucleotides that is complementary to at least one functional adaptive immune receptor J region-encoding gene segment and wherein the plurality of J-segment primers specifically hybridize to substantially all functional adaptive immune receptor J
region-encoding gene segments that are present in the template composition.
region-encoding gene segments that are present in the template composition.
[0095] The V-segment and J-segment oligonucleotide primers are capable of promoting amplification in said multiplex polymerase chain reaction (PCR) of substantially all template oligonucleotides in the template composition to produce a multiplicity of amplified template DNA molecules, said multiplicity of amplified template DNA molecules being sufficient to quantify diversity of the template oligonucleotides in the template composition, and wherein each amplified template DNA molecule in the multiplicity of amplified template DNA
molecules is less than 1000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80 or 70 nucleotides in length.
molecules is less than 1000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80 or 70 nucleotides in length.
[0096] The method also includes steps of: (b) sequencing all or a sufficient portion of each of said multiplicity of amplified template DNA molecules to determine, for each unique template DNA molecule in said multiplicity of amplified template DNA
molecules, (i) a template-specific oligonucleotide DNA sequence and (ii) a relative frequency of occurrence of the template oligonucleotide; and (c) comparing the relative frequency of occurrence for each unique template DNA sequence from said template composition, wherein a non-uniform frequency of occurrence for one or more template DNA sequences indicates non-uniform nucleic acid amplification potential among members of the set of oligonucleotide amplification primers.
molecules, (i) a template-specific oligonucleotide DNA sequence and (ii) a relative frequency of occurrence of the template oligonucleotide; and (c) comparing the relative frequency of occurrence for each unique template DNA sequence from said template composition, wherein a non-uniform frequency of occurrence for one or more template DNA sequences indicates non-uniform nucleic acid amplification potential among members of the set of oligonucleotide amplification primers.
[0097] Further description about bias control methods are provided in US
2013/0253842, US
9,150,905, US 2015/0203897, and WO 2013/169957 (PCT/U52013/040221), which are incorporated by reference in their entireties.
Sequencing
2013/0253842, US
9,150,905, US 2015/0203897, and WO 2013/169957 (PCT/U52013/040221), which are incorporated by reference in their entireties.
Sequencing
[0098] Sequencing can be performed using any of a variety of available high throughput single molecule sequencing machines and systems. Illustrative sequence systems include sequence-by-synthesis systems, such as the Illumina Genome Analyzer and associated instruments (Illumina HiSeq) (Illumina, Inc., San Diego, CA), Helicos Genetic Analysis System (Helicos BioSciences Corp., Cambridge, MA), Pacific Biosciences PacBio RS
(Pacific Biosciences, Menlo Park, CA), or other systems having similar capabilities.
Sequencing is achieved using a set of sequencing platform-specific oligonucleotides that hybridize to a defined region within the amplified DNA molecules. The sequencing platform-specific oligonucleotides are designed to sequence up amplicons, such that the V-and J- encoding gene segments can be uniquely identified by the sequences that are generated. See, e.g., US 2012/0058902; US 2010/033057; WO 2011/106738 (PCT/US2011/026373); US 2015/0299785; or WO 2012/027503 (PCT/US2011/049012), which is each incorporated by reference in its entirety.
(Pacific Biosciences, Menlo Park, CA), or other systems having similar capabilities.
Sequencing is achieved using a set of sequencing platform-specific oligonucleotides that hybridize to a defined region within the amplified DNA molecules. The sequencing platform-specific oligonucleotides are designed to sequence up amplicons, such that the V-and J- encoding gene segments can be uniquely identified by the sequences that are generated. See, e.g., US 2012/0058902; US 2010/033057; WO 2011/106738 (PCT/US2011/026373); US 2015/0299785; or WO 2012/027503 (PCT/US2011/049012), which is each incorporated by reference in its entirety.
[0099] In some embodiments, the raw sequence data is preprocessed to remove errors in the primary sequence of each read and to compress the data. A nearest neighbor algorithm can be used to collapse the data into unique sequences by merging closely related sequences, to remove both PCR and sequencing errors. See, e.g., US 2012/0058902; US
2010/033057;
WO 2011/106738 (PCT/U52011/026373); US 2015/0299785; or WO 2012/027503 (PCT/U52011/049012), which is each incorporated by reference in its entirety.
PCR template abundance estimation
2010/033057;
WO 2011/106738 (PCT/U52011/026373); US 2015/0299785; or WO 2012/027503 (PCT/U52011/049012), which is each incorporated by reference in its entirety.
PCR template abundance estimation
[00100] To estimate the average read coverage per input template in the multiplex PCR
and sequencing approach, a set of synthetic TCR (or BCR) templates (as described above) can be used, comprising each combination of VP and Jf3 gene segments. These synthetic molecules can be those described in general formula (I) above, and in US
2013/0253842, US
9,150,905, US 2015/0203897, and WO 2013/169957 (PCT/U52013/040221), which are incorporated by reference in their entireties.
and sequencing approach, a set of synthetic TCR (or BCR) templates (as described above) can be used, comprising each combination of VP and Jf3 gene segments. These synthetic molecules can be those described in general formula (I) above, and in US
2013/0253842, US
9,150,905, US 2015/0203897, and WO 2013/169957 (PCT/U52013/040221), which are incorporated by reference in their entireties.
[00101] These synthetic molecules can be included in each PCR reaction at very low concentration so that only some of the synthetic templates are observed. Using the known concentration of the synthetic template pool, the relationship between the number of observed unique synthetic molecules and the total number of synthetic molecules added to reaction can be simulated (this is very nearly one-to-one at the low concentrations that were used). The synthetic molecules allow calculation for each PCR reaction the mean number of sequencing reads obtained per molecule of PCR template, and an estimation of the number of T cells in the input material bearing each unique TCR rearrangement.
Discovery of diagnostic public T cell responses
Discovery of diagnostic public T cell responses
[00102] Given a large population of subjects with and without a given infectious disease, public T cell responses diagnostic for the disease state can be determined by applying the following statistical method. An immune receptor repertoire of unique T cell receptor sequences is determined, using the methods described above, for a group of subjects who have been identified as having or not having an infectious disease, such as CMV or small pox. It is possible, then, to determine which of those T cell receptor sequences are significantly more common (i.e., present in more individuals) among subjects with the disease state than in subjects without the disease state. When a common TCR
binds to an antigen, it is called a public T cell response. Public T cell responses are specific for a particular disease or antigen, are present in many individuals, and are encoded in a common format (specific rearranged receptor sequences) regardless of disease.
binds to an antigen, it is called a public T cell response. Public T cell responses are specific for a particular disease or antigen, are present in many individuals, and are encoded in a common format (specific rearranged receptor sequences) regardless of disease.
[00103] A one-tailed Fisher exact test is used, using presence or absence of T
cell receptor sequence in question vs. with or without disease state to construct a 2x2 contingency table (shown below) and with FDR (false discovery rate) controlled using an empirical distribution of null p-values determined by permutations of disease state. Since many clones are unique to a single subject (and consequently unique to either the positive or negative disease status classes), it is vital to control false discovery rate in feature selection to avoid over-fitting to the many spurious associations of unique TCRPs with positive infectious disease status. This process generates a list of T cell receptor sequences significantly more common in subjects with the disease state of interest.
Disease+ Disease-clone i present clone i not present N+ ¨ n,+ N_ Use of diagnostic public T cell responses to infer disease status
cell receptor sequence in question vs. with or without disease state to construct a 2x2 contingency table (shown below) and with FDR (false discovery rate) controlled using an empirical distribution of null p-values determined by permutations of disease state. Since many clones are unique to a single subject (and consequently unique to either the positive or negative disease status classes), it is vital to control false discovery rate in feature selection to avoid over-fitting to the many spurious associations of unique TCRPs with positive infectious disease status. This process generates a list of T cell receptor sequences significantly more common in subjects with the disease state of interest.
Disease+ Disease-clone i present clone i not present N+ ¨ n,+ N_ Use of diagnostic public T cell responses to infer disease status
[00104] Given a large population of subjects with and without a given infectious disease, and given a list of diagnostic public T cell responses generated as described above, it is possible to infer disease status in a subject of unknown disease status.
[00105] First, a 'disease burden' or quantitative measure of the presence and/or abundance of said diagnostic public T cell responses is calculated for each subject (those with known status as well as the subject of unknown status). The disease burden or quantitative measure is the proportion of unique T cell receptor sequences in each subject that are among the list of diagnostic public T cell responses.
[00106] Once this measure has been calculated, it is then determined whether the subject of unknown status has a disease burden consistent with known subjects who have the disease state of interest or with known subjects who do not have the disease state of interest. Herein, the method of comparison is to train a logistic regression model of (disease burden vs.
presence of disease state) in all subjects of known disease status, and use that model to assign to the subject of unknown status a probability of having the disease status of interest.
presence of disease state) in all subjects of known disease status, and use that model to assign to the subject of unknown status a probability of having the disease status of interest.
[00107] In some embodiments, the method comprises a classification model using an exhaustive leave-one-out cross-validation, in which the immune profile dataset from an individual is held out of the calculations, and the process is repeated from the beginning. In one aspect, cross-validation is used to assess the accuracy of a classification model. The resulting classification model is used to predict the infectious disease status of the holdout subject.
[00108] These methods can be performed for determining disease status for various types of disease, including but not limited to CMV, EBV, HPV, HIV, small pox, other infectious diseases, or even non-infectious diseases, like autoimmune diseases and malignancies.
HLA Typing
HLA Typing
[00109] HLA typing of bone marrow donor subjects was performed according to standard protocol by the Fred Hutchinson Cancer Research Center, and can be performed by methods known to those of skill in the art.
Method of Determining HLA Status of Unknown Subjects
Method of Determining HLA Status of Unknown Subjects
[00110] FIG. 13 shows an overview of the method for inferring HLA status. This method can be applied to any HLA allele type and any profile of TCR sequences, including TCRA, TCRB, TCRG, or TCRD gene sequences.
[00111] First, the HLA allele type is determined for each subject in a group, according to standard methods in the art for HLA typing. Then, a locus of the subject's immune repertoire is sequenced. For example, the TCRB locus is amplified and sequenced, using the methods described above.
[00112] For each subject, the total number of unique TCRB sequences and the frequency of each unique TCRB sequence are determined. For a given HLA allele, it is determined which TCRB sequences are significantly associated with the HLA allele. In one embodiment, for each unique TCRB sequence, it is determined how many subjects who are positive for an HLA allele have the TCRB sequence and how many subjects who are negative for the HLA allele have the TCRB sequence. In addition, it can be determined the number of subjects who are positive for an HLA allele and negative for the TCRB sequence and the number of subjects who are negative for the HLA allele and positive for the TCRB sequence.
The table below shows categorization of subjects by the presence or absence of a TCRB
sequence and the presence or absence of an HLA allele, HLA-A2.
HLA-A2+ HLA-A2-TCRf3 sequence i present n, n,_ TCRf3 sequence i not present N+ ¨ n, N_ ¨ n,_
The table below shows categorization of subjects by the presence or absence of a TCRB
sequence and the presence or absence of an HLA allele, HLA-A2.
HLA-A2+ HLA-A2-TCRf3 sequence i present n, n,_ TCRf3 sequence i not present N+ ¨ n, N_ ¨ n,_
[00113] FIG. 14 shows a list of exemplary unique TCRB sequences and the number of subjects who are positive or negative for an HLA-A2 allele that have a particular TCRB
sequence.
sequence.
[00114] A p-value is determined for the association of each TCRB sequence with an HLA
status using a Fisher exact test (two-tailed). The p value for association of each TCR with allele status using a Fisher exact test (two tailed) is calculated as follows:
(n+ +n_ (N++N_ ¨n+ ¨n_ Pr(table) = n+
(N++N
N
status using a Fisher exact test (two-tailed). The p value for association of each TCR with allele status using a Fisher exact test (two tailed) is calculated as follows:
(n+ +n_ (N++N_ ¨n+ ¨n_ Pr(table) = n+
(N++N
N
[00115] FIG. 14 also shows exemplary p values that are calculated for the association of a particular TCRB sequence with an HLA type (HLA-A2).
[00116] As shown in FIG. 15, a p value is selected as a cutoff for identifying a set of "Feature TCRS" from the entire list of possible TCR sequences. Defining a p-value threshold and permuting the allele status across individuals provides an estimate of false discovery rate. This is performed for each HLA allele, resulting in a set of allele-associated TCRf3 sequences for each HLA allele present in the training data. In FIG. 15, a p value cutoff of p<10-4 and an FDR of 0.1 is used to identify 288 TCRB sequences that are positively associated with HLA-A2. For each of the allele-associated TCRB sequences, the frequency of the sequence is also determined in each subject.
[00117] The feature selection step is followed by a machine learning process.
As shown in FIG. 16, for each HLA allele, a logistic regression model is trained using the set of feature vectors over all subjects, along with the known status for presence of that allele.
As shown in FIG. 16, for each HLA allele, a logistic regression model is trained using the set of feature vectors over all subjects, along with the known status for presence of that allele.
[00118] As shown in FIG. 17, an exhaustive leave-one-out cross validation is performed where one subject is removed from the analysis, and the HLA status of the subject is inferred based on feature selection and training from only the remaining subjects. The result is a set of classifiers (one for each HLA allele) that estimate the probability of positive status for each HLA allele, taking as input the feature vector for each allele.
[00119] Equipped with these classifiers, the HLA type of a new subject can be assessed by: 1. immune repertoire sequencing, 2. computing feature vectors for each allele, and 3.
defining a probability threshold for positive status that will be applied to the output from each classifier. In this manner, an HLA status can be inferred for a new subject with an unknown HLA status.
defining a probability threshold for positive status that will be applied to the output from each classifier. In this manner, an HLA status can be inferred for a new subject with an unknown HLA status.
[00120] In an embodiment, given the results of the feature selection scheme, a one-dimension allele score is determined as the fraction of an individual's unique TCRf3 sequences that appear in the set of TCRf3 sequences associated with a given allele. This single quantity, rather than a vector of features, may also be used for training a classifier and diagnosing novel individuals with similar accuracy.
[00121] Additionally, the specific classifier described (logistic regression) can be replaced with any of a number of other binary classifiers and gives substantially similar results. These include: k-nearest neighbors, random forests, artificial neural network, naive Bayes, and support vector machine.
EXAMPLES
EXAMPLES
[00122] Below are examples of specific embodiments for carrying out the present invention. The examples are offered for illustrative purposes only, and are not intended to limit the scope of the present invention in any way. Efforts have been made to ensure accuracy with respect to numbers used (e.g., amounts, temperatures, etc.), but some experimental error and deviation should, of course, be allowed for.
[00123] The practice of the present invention will employ, unless otherwise indicated, conventional methods of protein chemistry, biochemistry, recombinant DNA
techniques and pharmacology, within the skill of the art. Such techniques are explained fully in the literature. See, e.g., T.E. Creighton, Proteins: Structures and Molecular Properties (W.H.
Freeman and Company, 1993); A.L. Lehninger, Biochemistry (Worth Publishers, Inc., current addition); Sambrook, et al., Molecular Cloning: A Laboratory Manual (2nd Edition, 1989);
Methods In Enzymology (S. Colowick and N. Kaplan eds., Academic Press, Inc.);
Remington's Pharmaceutical Sciences, 18th Edition (Easton, Pennsylvania: Mack Publishing Company, 1990); Carey and Sundberg Advanced Organic Chemistry .3rd Ed. (Plenum Press) Vols A and B(1992).
Example 1: Identification of Cytomegalovirus (CMV)-Associated TCRs and Classification of CMV Status From Immunosequencing Data
techniques and pharmacology, within the skill of the art. Such techniques are explained fully in the literature. See, e.g., T.E. Creighton, Proteins: Structures and Molecular Properties (W.H.
Freeman and Company, 1993); A.L. Lehninger, Biochemistry (Worth Publishers, Inc., current addition); Sambrook, et al., Molecular Cloning: A Laboratory Manual (2nd Edition, 1989);
Methods In Enzymology (S. Colowick and N. Kaplan eds., Academic Press, Inc.);
Remington's Pharmaceutical Sciences, 18th Edition (Easton, Pennsylvania: Mack Publishing Company, 1990); Carey and Sundberg Advanced Organic Chemistry .3rd Ed. (Plenum Press) Vols A and B(1992).
Example 1: Identification of Cytomegalovirus (CMV)-Associated TCRs and Classification of CMV Status From Immunosequencing Data
[00124] The first aim of the study was to create a comprehensive catalog of CMV-specific T cell receptor sequences. The goal was to identify both common and rare, shared TCRI3 sequences enriched in CMV seropositive individuals relative to CMV
seronegative individuals. CMV was used as a model because it provides an ideal test bed for the development of T cell biomarkers. CMV serostatus is widely available and easy access to near-equal numbers of seropositive (cases) and seronegative (control) individuals among healthy adults provides good statistical power for a rigorous analysis.
seronegative individuals. CMV was used as a model because it provides an ideal test bed for the development of T cell biomarkers. CMV serostatus is widely available and easy access to near-equal numbers of seropositive (cases) and seronegative (control) individuals among healthy adults provides good statistical power for a rigorous analysis.
[00125] There are many possible TCR sequences generated during VDJ
recombination (-1012 easily accessible by VDJ recombination). As shown in FIG. 1, different VDJ
recombination events may result in the same nucleotide sequence. Moreover, it is possible that VDJ recombination can result in different nucleotide sequences that are translated into the same amino acid sequence. See Venturi, Price, Douek and Davenport 2008.
The molecular basis for public T-cell responses? Nature Reviews Immunology.
Certain TCRs are common rearrangements and shared by more than one individual. When a common TCR
binds to an antigen, it is called a public T cell response. Public T cell responses are specific for a particular disease or antigen, are present in many individuals, and are encoded in a common format (specific rearranged receptor sequences) regardless of disease.
recombination (-1012 easily accessible by VDJ recombination). As shown in FIG. 1, different VDJ
recombination events may result in the same nucleotide sequence. Moreover, it is possible that VDJ recombination can result in different nucleotide sequences that are translated into the same amino acid sequence. See Venturi, Price, Douek and Davenport 2008.
The molecular basis for public T-cell responses? Nature Reviews Immunology.
Certain TCRs are common rearrangements and shared by more than one individual. When a common TCR
binds to an antigen, it is called a public T cell response. Public T cell responses are specific for a particular disease or antigen, are present in many individuals, and are encoded in a common format (specific rearranged receptor sequences) regardless of disease.
[00126] In this study, the public T cell response to cytomegalovirus (CMV) was examined by sequencing of rearranged T cell receptor genes (TCRs) obtained from 640 subjects. CMV
is a life-long infection, usually asymptomatic, that afflicts most adults and elicits a robust memory T cell response.
is a life-long infection, usually asymptomatic, that afflicts most adults and elicits a robust memory T cell response.
[00127] FIGs. 2A, 2B and 2C provide an overview of the method, according to one embodiment of the invention. FIG 2A shows a dataset of peripheral blood samples from 640 healthy subjects (287 CMV- and 353 CMV+), which were analyzed by high-throughput TCR
profiling. In FIG. 2B, unique TCRf3 sequences were identified that were present in significantly more CMV+ subjects than CMV- subjects, controlling false determination rate (FDR) by permutation of CMV status. Presence of these CMV-associated TCR0 sequences was used to build a classification model. In FIG. 2C, the classification model was tested using exhaustive leave-one-out cross-validation, in which one sample was held out of the calculations, and the process was repeated from the beginning. The resulting classification model was used to predict the CMV status of the holdout subject.
profiling. In FIG. 2B, unique TCRf3 sequences were identified that were present in significantly more CMV+ subjects than CMV- subjects, controlling false determination rate (FDR) by permutation of CMV status. Presence of these CMV-associated TCR0 sequences was used to build a classification model. In FIG. 2C, the classification model was tested using exhaustive leave-one-out cross-validation, in which one sample was held out of the calculations, and the process was repeated from the beginning. The resulting classification model was used to predict the CMV status of the holdout subject.
[00128] Human peripheral blood samples were obtained from the Fred Hutchinson Cancer Research Center Research Cell Bank biorepository of healthy bone marrow donors under a protocol following written informed consent approved and supervised by the Fred Hutchinson Cancer Research Center Institutional Review Board. This biorepository houses an inventory of PBMC, B-LCL and DNA from hematopoietic cell transplant (HCT) patients, donors and family members along with cell lines and DNA derived from stem cell transplant patients, donors and selected family members from ¨7,800 stem cell transplant patients and 7,300 donors (from which samples were drawn). The DNA for this study was extracted from peripheral blood of HCT donors and extensively typed for HLA antigens and alleles along with other defining phenotypes needed for donor-patient matching or donor inclusion/exclusion. Among these are cytomegalovirus (CMV), Epstein-Barr virus (EBV), Herpes Simplex Virus (HSV), hepatitis, and diabetes mellitus.
[00129] As shown in FIG. 3A, 640 subjects were specifically phenotyped for HLA
type and CMV status for eligibility as a HCT donor, and segregated between 287 CMV+
donors and 353 CMV- donors. The CMV seropositive rate was approximately 45% (roughly equal numbers of seropositive and seronegative samples for the investigation of public T cell responses).
type and CMV status for eligibility as a HCT donor, and segregated between 287 CMV+
donors and 353 CMV- donors. The CMV seropositive rate was approximately 45% (roughly equal numbers of seropositive and seronegative samples for the investigation of public T cell responses).
[00130] FIG. 3B shows demographic characteristics of the subjects included in this study, classified by CMV status.
High-throughput TCRfl sequencing:
High-throughput TCRfl sequencing:
[00131] Genomic DNA was extracted from peripheral blood samples using the Qiagen DNeasy Blood extraction Kit (Qiagen, Gaithersburg, MD, USA). The CDR3 regions of rearranged TCRB genes were sequenced; the TCRB CDR3 region was defined according to the IMGT collaboration28' 20. TCRB CDR3 regions were amplified and sequenced using methods described above and in previously described protocols5' 29. The multiplexed PCR
method used a mixture of 60 forward primers specific to TCR VB gene segments and 13 reverse primers specific to TCR J0 gene segments. The resulting amplicons were sequenced using the methods described above. Reads of 87 bp were obtained using the Illumina HiSeq System. Raw HiSeq sequence data were preprocessed to remove errors in the primary sequence of each read, and to compress the data. A nearest neighbor algorithm was used to collapse the data into unique sequences by merging closely related sequences, to remove both PCR and sequencing errors.
method used a mixture of 60 forward primers specific to TCR VB gene segments and 13 reverse primers specific to TCR J0 gene segments. The resulting amplicons were sequenced using the methods described above. Reads of 87 bp were obtained using the Illumina HiSeq System. Raw HiSeq sequence data were preprocessed to remove errors in the primary sequence of each read, and to compress the data. A nearest neighbor algorithm was used to collapse the data into unique sequences by merging closely related sequences, to remove both PCR and sequencing errors.
[00132] In order to ensure adequate coverage of each T cell rearrangement, 8-10x sequence coverage, or ¨6-10 million sequencing reads per sample, was generated using approximately eight full sequencing runs. All sequencing reads were processed using a standardized bioinformatics pipeline to 1) demultiplex reads to specific samples, 2) eliminate low quality sequence and remove potential contaminants, 3) align and identify specific TCRB
V and J gene segments and CDR3 regions, 4) cluster highly similar sequences to account for PCR and sequencing errors, 5) normalize data to remove PCR amplification bias, 6) estimate total T-cell input, and 7) generate TCRB unique sequence counts and distributions.
V and J gene segments and CDR3 regions, 4) cluster highly similar sequences to account for PCR and sequencing errors, 5) normalize data to remove PCR amplification bias, 6) estimate total T-cell input, and 7) generate TCRB unique sequence counts and distributions.
[00133] Here, approximately 250,000 rearranged T cell receptor genes were sequenced from peripheral blood of each subject.
Feature selection:
Feature selection:
[00134] Across 640 subjects, there were 185,204 (+1¨ 84,171) unique TCRPs per subject, and 83,727,796 unique TCRPs in aggregate. Rather than attempting high dimensional CMV
classification using all unique TCRPs as potential features, a novel feature selection scheme was developed. Feature selection was selection of a particular amino acid sequence to identify common TCRs among individuals with CMV.
classification using all unique TCRPs as potential features, a novel feature selection scheme was developed. Feature selection was selection of a particular amino acid sequence to identify common TCRs among individuals with CMV.
[00135] Since many clones were unique to a single subject (and consequently unique to either the CMV+ or CMV¨ classes), it was vital to control false discovery rate in feature selection to avoid over-fitting to the many spurious associations of unique TCRPs with CMV
status.
status.
[00136] Each unique TCRf3 rearrangement, identified by V and J gene assignment and CDR3 amino acid sequence, was tested for CMV association. Each of these was subjected to a one-tailed Fisher exact test for its incidence in CMV¨ and CMV+ subjects.
Specifically, letting nu denote the number of subjects with CMV status j (with j ¨ or +) and clone i present, a p-valuep, was computed by performing Fisher's exact test on the contingency table.
CMV+ CMV-TCRP sequence i present TCRP sequence i not present N ¨ ni+ N_ ¨
where N+ and N_ denote the total number of subjects having each CMV status (CMV+ and CMV- respectively).
Specifically, letting nu denote the number of subjects with CMV status j (with j ¨ or +) and clone i present, a p-valuep, was computed by performing Fisher's exact test on the contingency table.
CMV+ CMV-TCRP sequence i present TCRP sequence i not present N ¨ ni+ N_ ¨
where N+ and N_ denote the total number of subjects having each CMV status (CMV+ and CMV- respectively).
[00137] To characterize a rejection region in the presence of many weakly dependent hypotheses (one for each unique TCR(3), CMV status assignments were randomly permuted 100 times, and statistics on the number of rejections at the nominal p-value threshold were calculated. Approximating the total fraction of true null hypotheses as unity, this allowed estimation of the false discovery rate (FDR) as the ratio of the mean number of rejections under permutation to the actual number of rejections.
[00138] As shown by example in FIG. 4, the TCR amino acid sequence "CASSLIGVSSYNEQFF" (SEQ ID NO:12) was selected and was identified in 27 CMV+
subjects and 2 CMV- subjects (p= 2.8E-08).
subjects and 2 CMV- subjects (p= 2.8E-08).
[00139] FIG. 5 shows exemplary TCR amino acid sequences identified in CMV+
subjects.
As shown in Table 1, at a selected p value of 104, 142 public T cell clones were identified as associated with CMV+ subjects. The p value was chosen at 10-4, but as one can see, at a lower p value, larger numbers of CMV associated TCRB sequences can be identified.
subjects.
As shown in Table 1, at a selected p value of 104, 142 public T cell clones were identified as associated with CMV+ subjects. The p value was chosen at 10-4, but as one can see, at a lower p value, larger numbers of CMV associated TCRB sequences can be identified.
[00140] Table 1: Each of the 142 TCRf3 sequences significantly associated with CMV
status (p<1x10-4, FDR-20%) in the cohort was included. Sequences were defined using the amino acid sequence of the CDR3 region along with V and J gene segments. The number of CMV+ and CMV- subjects in which each sequence was observed is given, as well as the p-value from a Fisher exact test for association with CMV status, and any HLA-A
or HLA-B
alleles that were significantly associated (p < 1x10-3) with the presence of each sequence.
Table 1. Exemplary CMV-Associated TCRI3 Sequences SEQ AA sequence V gene J gene CMV+/CMV- Fisher p value HLA
ID
association(s) NO:
1 CASSGQGAYEQYF TCRBV09-01 TCRBJ02-07*03 60/11 3.3142E-13 A29,B13 2 CASSPDRVGQETQYF TCRBV05-01*01 TCRBJ02-05*03 34/1 9.5395E-12 None 3 CASSIGPLEHNEQFF TCRBV19-01 TCRBJ02-01*03 30/0 1.4752E-11 Al,B8 4 CASSIEGNQPQHF TCRBV28-01*01 TCRBJ01-05*03 27/0 1.9463E-10 None CASSLVAGGRETQYF TCRBV05-06*01 TCRBJ02-05*03 51/11 2.3538E-6 CASSLEAEYEQYF TCRBV07-02*01 TCRBJ02-07*03 30/1 2.7085E-10 B8 7 CATSDGDEQFF TCRBV24 TCRBJ02-01*03 41/6 4.7339E-10 Al 8 CASSLAPGATNEKLFF TCRBV07-06*01 TCRBJ01-04*03 36/4 9.134E-10 A2 9 CASSRGRQETQYF TCRBV07-06*01 TCRBJ02-05*03 37/5 2.1095E-09 None CASSAGQGVTYEQYF TCRBV09-01 TCRBJ02-07*03 23/0 5.8832E-09 Al,B8 11 CASSLRREKLFF TCRBV05-06*01 TCRBJ01-04*03 28/2 1.2553E-08 A29,B50 12 CASSLIGVSSYNEQFF TCRBV07-09 TCRBJ02-01*03 27/2 2.7588E-08 A3,B7 13 CASSFPTSGQETQYF TCRBV07-09 TCRBJ02-05*03 27/2 2.7588E-08 A24 14 CASSPQRNTEAFF TCRBV04-03*01 TCRBJ01-01*03 27/2 2.7588E-08 A3,B7 CSVRDNFNQPQHF TCRBV29-01*01 TCRBJ01-05*03 21/0 3.1927E-08 None 16 CATSRDSQGSYGYTF TCRBV15-01*01 TCRBJ01-02*03 21/0 3.1927E-08 None 17 CASSPGDEQYF TCRBV25-01*01 TCRBJ02-07*03 24/1 3.7138E-08 B7 18 CASSQTGGRNQPQHF TCRBV12 TCRBJ01-05*03 26/2 6.0357E-08 19 CASSQNRGQETQYF TCRBV14-01*01 TCRBJ02-05*03 26/2 6.0357E-08 None CASSLVIGGDTEAFF TCRBV05-01*01 TCRBJ01-01*03 26/2 6.0357E-08 21 CASSFHGFNQPQHF TCRBV05-06*01 TCRBJ01-05*03 20/0 7.4139E-08 Al,B8 22 CASSRLAGGTDTQYF TCRBV07-03*01 TCRBJ02-03*03 50/16 7.5065E-08 Al,B8 23 CASSLPSGLTDTQYF TCRBV28-01*01 TCRBJ02-03*03 25/2 1.3144E-07 None 24 CATSRDTQGSYGYTF TCRBV15-01*01 TCRBJ01-02*03 19/0 1.7179E-07 None CATSDGDTQYF TCRBV24 TCRBJ02-03*03 51/18 2.2203E-07 Al,B8 26 CASSLVASGRETQYF TCRBV05-06*01 TCRBJ02-05*03 24/2 2.8483E-07 None 27 CASSIWGLDTEAFF TCRBV19-01 TCRBJ01-01*03 28/4 3.6207E-07 None 28 CASSPGDEQFF TCRBV25-01*01 TCRBJ02-01*03 28/4 3.6207E-07 B7 29 CASSPSTGTEAFF TCRBV05-06*01 TCRBJ01-01*03 18/0 3.9723E-07 None CSVEEDEGIYGYTF TCRBV29-01*01 TCRBJ01-02*03 18/0 3.9723E-07 None 31 CASSEIPNTEAFF TCRBV06-04 TCRBJ01-01*03 18/0 3.9723E-07 32 CASSQVPGQGDNEQFF TCRBV14-01*01 TCRBJ02-01*03 21/1 4.1428E-07 Al,B8 33 CASSPAGLNTEAFF TCRBV19-01 TCRBJ01-01*03 21/1 4.1428E-07 None 34 CASSLGLKGTQYF TCRBV12 TCRBJ02-05*03 21/1 4.1428E-07 None CASSGDRLYEQYF TCRBV02-01*01 TCRBJ02-07*03 23/2 6.1415E-07 None 36 CSVRDNYNQPQHF TCRBV29-01*01 TCRBJ01-05*03 23/2 6.1415E-07 None 37 CASSYGGLGSYEQYF TCRBV06-05*01 TCRBJ02-07*03 23/2 6.1415E-07 None 38 CASNRDRGRYEQYF TCRBV06-01*01 TCRBJ02-07*03 20/1 9.1836E-07 B45 39 CASMGGASYEQYF TCRBV27-01*01 TCRBJ02-07*03 20/1 9.1836E-07 A2 40 CASSLGVGPYNEQFF TCRBV07-02*01 TCRBJ02-01*03 20/1 9.1836E-07 B7 41 CASSLGGAGDTQYF TCRBV12 TCRBJ02-03*03 44/15 1.1802E-06 None 42 CATSRGTVSYEQYF TCRBV15-01*01 TCRBJ02-07*03 30/6 1.2303E-06 None 43 CATSDGETQYF TCRBV24 TCRBJ02-05*03 70/36 1.3031E-06 A1,B52 44 CASSEARGGVEKLFF TCRBV06-01*01 TCRBJ01-04*03 22/2 1.3173E-06 None 45 CASSLNRGQETQYF TCRBV14-01*01 TCRBJ02-05*03 19/1 2.0272E-06 None 46 CSVRDNHNQPQHF TCRBV29-01*01 TCRBJ01-05*03 19/1 2.0272E-06 None 47 CASSESGHRNQPQHF TCRBV10-02*01 TCRBJ01-05*03 16/0 2.1105E-06 None 48 CSASPGQGASYGYTF TCRBV20 TCRBJ01-02*03 16/0 2.1105E-06 None 49 CASSEARTRAFF TCRBV06-01*01 TCRBJ01-01*03 16/0 2.1105E-06 None 50 CASRPTGYEQYF TCRBV06-01*01 TCRBJ02-07*03 21/2 2.8099E-06 B39 51 CASSVTGGTDTQYF TCRBV09-01 TCRBJ02-03*03 74/41 2.8214E-06 A1,B8 52 CASSRLAASTDTQYF TCRBV07-03*01 TCRBJ02-03*03 23/3 3.1389E-06 A1,B8 53 CATSDSVTNTGELFF TCRBV24 TCRBJ02-02*03 18/1 4.4551E-06 54 CASSRNRESNQPQHF TCRBV06-05*01 TCRBJ01-05*03 18/1 4.4551E-06 None 55 CASSAQGAYEQYF TCRBV09-01 TCRBJ02-07*03 28/6 4.717E-06 None 56 CASSIQGYSNQPQHF TCRBV05-08*01 TCRBJ01-05*03 15/0 4.8494E-06 B8 57 CASSYNPYSNQPQHF TCRBV06-06 TCRBJ01-05*03 15/0 4.8494E-06 None 58 CASSLGHRDPNTGELFF TCRBV05-01*01 TCRBJ02-02*03 15/0 4.8494E-06 59 CASSTTGGDGYTF TCRBV19-01 TCRBJ01-02*03 26/5 5.5727E-06 None 60 CASSVLAGPTDTQYF TCRBV09-01 TCRBJ02-03*03 26/5 5.5727E-06 Al 61 CASSYRQETQYF TCRBV06-05*01 TCRBJ02-05*03 20/2 5.9592E-06 None 62 CASSSGQVYGYTF TCRBV05-06*01 TCRBJ01-02*03 22/3 6.4736E-06 Al,B8 63 CASGRDTYEQYF TCRBV02-01*01 TCRBJ02-07*03 17/1 9.7453E-06 None 64 CATSDSRTGGQETQYF TCRBV24 TCRBJ02-05*03 17/1 9.7453E-06 B8 65 CASSSPGRSGANVLTF TCRBV28-01*01 TCRBJ02-06*03 17/1 9.7453E-06 None 66 CASSYGGEGYTF TCRBV06-05*01 TCRBJ01-02*03 36/12 1.059E-05 None 67 CASSLAGVDYEQYF TCRBV07-09 TCRBJ02-07*03 25/5 1.0951E-05 B8 68 CASSLQGADTQYF TCRBV07-08*01 TCRBJ02-03*03 14/0 1.112E-05 None 69 CASSLEAENEQFF TCRBV07-02*01 TCRBJ02-01*03 14/0 1.112E-05 All,B8 70 CASSEAPSTSTDTQYF TCRBV02-01*01 TCRBJ02-03*03 14/0 1.112E-05 None 71 CASSLEGQQPQHF TCRBV28-01*01 TCRBJ01-05*03 14/0 1.112E-05 None 72 CASSLGHRDSSYEQYF TCRBV05-01*01 TCRBJ02-07*03 14/0 1.112E-05 A29,B57 73 CASSPSRNTEAFF TCRBV04-03*01 TCRBJ01-01*03 19/2 1.2561E-05 A3,B7 74 CSALGHSNQPQHF TCRBV20 TCRBJ01-05*03 19/2 1.2561E-05 None 75 CASSHRDRNYEQYF TCRBV07-09 TCRBJ02-07*03 19/2 1.2561E-05 A24 76 CASSPPGQGSDTQYF TCRBV18-01*01 TCRBJ02-03*03 28/7 1.3913E-05 None 77 CASSLQGYSNQPQHF TCRBV05-08*01 TCRBJ01-05*03 34/11 1.4711E-05 B8 78 CASSYVRTGGNYGYTF TCRBV06-05*01 TCRBJ01-02*03 31/9 1.5052E-05 None 79 CASSRDRNYGYTF TCRBV06-04 TCRBJ01-02*03 29/8 2.0024E-05 None 80 CASSTGTSGSYEQYF TCRBV06-01*01 TCRBJ02-07*03 16/1 2.1213E-05 B54 81 CASRSDSGANVLTF TCRBV06-04 TCRBJ02-06*03 16/1 2.1213E-05 None 82 CATSRVAGETQYF TCRBV15-01*01 TCRBJ02-05*03 24/5 2.1345E-05 B7 83 CASSEEGIQPQHF TCRBV02-01*01 TCRBJ01-05*03 22/4 2.4814E-05 None 84 CASSLGGPGDTQYF TCRBV12 TCRBJ02-03*03 22/4 2.4814E-05 None 85 CASSLVAAGRETQYF TCRBV05-06*01 TCRBJ02-05*03 13/0 2.5446E-05 B8 86 CASRGQGWDEKLFF TCRBV06-05*01 TCRBJ01-04*03 13/0 2.5446E-05 None 87 CASSLEGQGFGYTF TCRBV05-01*01 TCRBJ01-02*03 13/0 2.5446E-05 None 88 CASRDWDYTDTQYF TCRBV02-01*01 TCRBJ02-03*03 13/0 2.5446E-05 None 89 CASSRSGLAGNTGELFF TCRBV06 TCRBJ02-02*03 13/0 2.5446E-05 None 90 CASSPGQEAGANVLTF TCRBV05-01*01 TCRBJ02-06*03 13/0 2.5446E-05 None 91 CASSLGDRPDTQYF TCRBV11-02*02 TCRBJ02-03*03 13/0 2.5446E-05 None 92 CASSFPGGETQYF TCRBV11-01*01 TCRBJ02-05*03 18/2 2.6306E-05 None 93 CASSLETYGYTF TCRBV05-06*01 TCRBJ01-02*03 18/2 2.6306E-05 B78 94 CASSSGQVQETQYF TCRBV11-02*02 TCRBJ02-05*03 18/2 2.6306E-05 B51 95 CASSFDNYGYTF TCRBV05-04*01 TCRBJ01-02*03 18/2 2.6306E-05 None 96 CASSEGARQPQHF TCRBV10-02*01 TCRBJ01-05*03 18/2 2.6306E-05 None 97 CASSLTGGRNQPQHF TCRBV12 TCRBJ01-05*03 33/11 2.6324E-05 98 CASSLLWDQPQHF TCRBV05-05*01 TCRBJ01-05*03 20/3 2.6954E-05 None 99 CASSLFGTGGNTEAFF TCRBV05-06*01 TCRBJ01-01*03 20/3 2.6954E-05 B44 100 CASSISAGEAFF TCRBV19-01 TCRBJ01-01*03 20/3 2.6954E-05 None 101 CASSPPSGLTDTQYF TCRBV28-01*01 TCRBJ02-03*03 20/3 2.6954E-05 None 102 CASSPLSDTQYF TCRBV07-09 TCRBJ02-03*03 30/9 2.752E-05 None 103 CASSSRGTGELFF TCRBV28-01*01 TCRBJ02-02*03 25/6 3.3358E-05 None 104 CASSYAGDGYTF TCRBV06-05*01 TCRBJ01-02*03 31/10 3.6414E-05 B8 105 CASSDRGNTGELFF TCRBV04-01*01 TCRBJ02-02*03 31/10 3.6414E-05 A68,B41 106 CASSPGGTQYF TCRBV12 TCRBJ02-05*03 28/8 3.6921E-05 None 107 CASSLGDRAYNEQFF TCRBV05-06*01 TCRBJ02-01*03 28/8 3.6921E-05 None 108 CASSLRGSSYNEQFF TCRBV05-08*01 TCRBJ02-01*03 23/5 4.1254E-05 None 109 CASSLTASSYEQYF TCRBV05-01*01 TCRBJ02-07*03 23/5 4.1254E-05 None 110 CSASDHEQYF TCRBV20-01*01 TCRBJ02-07*03 23/5 4.1254E-05 None 111 CASSQGRHTDTQYF TCRBV14-01*01 TCRBJ02-03*03 15/1 4.5934E-05 A68 112 CASSRPGQGNTEAFF TCRBV12 TCRBJ01-01*03 26/7 4.8714E-05 None 113 CASSLVGDGYTF TCRBV07-08*01 TCRBJ01-02*03 36/14 4.9771E-05 None 114 CASSLGAGNQPQHF TCRBV28-01*01 TCRBJ01-05*03 29/9 4.9861E-05 None 115 CASSLTDTGELFF TCRBV11-02*02 TCRBJ02-02*03 29/9 4.9861E-05 None 116 CASSLTGGNSGNTIYF TCRBV07-02*01 TCRBJ01-03*03 19/3 5.4373E-05 B14 117 CAWRGTGNSPLHF TCRBV30-01*01 TCRBJ01-06*03 17/2 5.4712E-05 None 118 CASASANYGYTF TCRBV12 TCRBJ01-02*03 12/0 5.8109E-05 119 CASSLQAGANEQFF TCRBV07-02*01 TCRBJ02-01*03 12/0 5.8109E-05 None 120 CASSEEAGGSGYTF TCRBV06-01*01 TCRBJ01-02*03 12/0 5.8109E-05 None 121 CASRTGESGYTF TCRBV06-05*01 TCRBJ01-02*03 12/0 5.8109E-05 None 122 CASSGLNEQFF TCRBV06-01*01 TCRBJ02-01*03 12/0 5.8109E-05 None 123 CASSWDRDNSPLHF TCRBV25-01*01 TCRBJ01-06*03 12/0 5.8109E-05 B7 124 CASSIRTNYYGYTF TCRBV19-01 TCRBJ01-02*03 12/0 5.8109E-05 None 125 CSARSPEAFF TCRBV20-01*01 TCRBJ01-01*03 12/0 5.8109E-05 None 126 CASSRGTGATDTQYF TCRBV19-01 TCRBJ02-03*03 12/0 5.8109E-05 None 127 CASSPRVSNQPQHF TCRBV12 TCRBJ01-05*03 12/0 5.8109E-05 None 128 CAISESQDRGHEQYF TCRBV10-03*01 TCRBJ02-07*03 12/0 5.8109E-05 None 129 CASSLGRGYEKLFF TCRBV05-06*01 TCRBJ01-04*03 12/0 5.8109E-05 B62 130 CSVEVRGTDTQYF TCRBV29-01*01 TCRBJ02-03*03 12/0 5.8109E-05 None 131 CASRGQGAGELFF TCRBV02-01*01 TCRBJ02-02*03 30/10 6.4924E-05 None 132 CATSREGSGYEQYF TCRBV15-01*01 TCRBJ02-07*03 22/5 7.9023E-05 None 133 CASSLGWTEAFF TCRBV05-01*01 TCRBJ01-01*03 22/5 7.9023E-05 None 134 CASSLGSSSYNEQFF TCRBV07-02*01 TCRBJ02-01*03 35/14 8.4988E-05 None 135 CASSSAGADTQYF TCRBV07-09 TCRBJ02-03*03 14/1 9.8907E-05 B5 136 CASSERKNYGYTF TCRBV06-01*01 TCRBJ01-02*03 14/1 9.8907E-05 None 137 CASRDRDRVNTEAFF TCRBV06-01*01 TCRBJ01-01*03 14/1 9.8907E-05 A3 138 CASSRVGEQFF TCRBV03 TCRBJ02-01*03 14/1 9.8907E-05 None 139 CASSPRWQETQYF TCRBV27-01*01 TCRBJ02-05*03 14/1 9.8907E-05 B14 140 CASTPGDTIYF TCRBV25-01*01 TCRBJ01-03*03 14/1 9.8907E-05 B55
status (p<1x10-4, FDR-20%) in the cohort was included. Sequences were defined using the amino acid sequence of the CDR3 region along with V and J gene segments. The number of CMV+ and CMV- subjects in which each sequence was observed is given, as well as the p-value from a Fisher exact test for association with CMV status, and any HLA-A
or HLA-B
alleles that were significantly associated (p < 1x10-3) with the presence of each sequence.
Table 1. Exemplary CMV-Associated TCRI3 Sequences SEQ AA sequence V gene J gene CMV+/CMV- Fisher p value HLA
ID
association(s) NO:
1 CASSGQGAYEQYF TCRBV09-01 TCRBJ02-07*03 60/11 3.3142E-13 A29,B13 2 CASSPDRVGQETQYF TCRBV05-01*01 TCRBJ02-05*03 34/1 9.5395E-12 None 3 CASSIGPLEHNEQFF TCRBV19-01 TCRBJ02-01*03 30/0 1.4752E-11 Al,B8 4 CASSIEGNQPQHF TCRBV28-01*01 TCRBJ01-05*03 27/0 1.9463E-10 None CASSLVAGGRETQYF TCRBV05-06*01 TCRBJ02-05*03 51/11 2.3538E-6 CASSLEAEYEQYF TCRBV07-02*01 TCRBJ02-07*03 30/1 2.7085E-10 B8 7 CATSDGDEQFF TCRBV24 TCRBJ02-01*03 41/6 4.7339E-10 Al 8 CASSLAPGATNEKLFF TCRBV07-06*01 TCRBJ01-04*03 36/4 9.134E-10 A2 9 CASSRGRQETQYF TCRBV07-06*01 TCRBJ02-05*03 37/5 2.1095E-09 None CASSAGQGVTYEQYF TCRBV09-01 TCRBJ02-07*03 23/0 5.8832E-09 Al,B8 11 CASSLRREKLFF TCRBV05-06*01 TCRBJ01-04*03 28/2 1.2553E-08 A29,B50 12 CASSLIGVSSYNEQFF TCRBV07-09 TCRBJ02-01*03 27/2 2.7588E-08 A3,B7 13 CASSFPTSGQETQYF TCRBV07-09 TCRBJ02-05*03 27/2 2.7588E-08 A24 14 CASSPQRNTEAFF TCRBV04-03*01 TCRBJ01-01*03 27/2 2.7588E-08 A3,B7 CSVRDNFNQPQHF TCRBV29-01*01 TCRBJ01-05*03 21/0 3.1927E-08 None 16 CATSRDSQGSYGYTF TCRBV15-01*01 TCRBJ01-02*03 21/0 3.1927E-08 None 17 CASSPGDEQYF TCRBV25-01*01 TCRBJ02-07*03 24/1 3.7138E-08 B7 18 CASSQTGGRNQPQHF TCRBV12 TCRBJ01-05*03 26/2 6.0357E-08 19 CASSQNRGQETQYF TCRBV14-01*01 TCRBJ02-05*03 26/2 6.0357E-08 None CASSLVIGGDTEAFF TCRBV05-01*01 TCRBJ01-01*03 26/2 6.0357E-08 21 CASSFHGFNQPQHF TCRBV05-06*01 TCRBJ01-05*03 20/0 7.4139E-08 Al,B8 22 CASSRLAGGTDTQYF TCRBV07-03*01 TCRBJ02-03*03 50/16 7.5065E-08 Al,B8 23 CASSLPSGLTDTQYF TCRBV28-01*01 TCRBJ02-03*03 25/2 1.3144E-07 None 24 CATSRDTQGSYGYTF TCRBV15-01*01 TCRBJ01-02*03 19/0 1.7179E-07 None CATSDGDTQYF TCRBV24 TCRBJ02-03*03 51/18 2.2203E-07 Al,B8 26 CASSLVASGRETQYF TCRBV05-06*01 TCRBJ02-05*03 24/2 2.8483E-07 None 27 CASSIWGLDTEAFF TCRBV19-01 TCRBJ01-01*03 28/4 3.6207E-07 None 28 CASSPGDEQFF TCRBV25-01*01 TCRBJ02-01*03 28/4 3.6207E-07 B7 29 CASSPSTGTEAFF TCRBV05-06*01 TCRBJ01-01*03 18/0 3.9723E-07 None CSVEEDEGIYGYTF TCRBV29-01*01 TCRBJ01-02*03 18/0 3.9723E-07 None 31 CASSEIPNTEAFF TCRBV06-04 TCRBJ01-01*03 18/0 3.9723E-07 32 CASSQVPGQGDNEQFF TCRBV14-01*01 TCRBJ02-01*03 21/1 4.1428E-07 Al,B8 33 CASSPAGLNTEAFF TCRBV19-01 TCRBJ01-01*03 21/1 4.1428E-07 None 34 CASSLGLKGTQYF TCRBV12 TCRBJ02-05*03 21/1 4.1428E-07 None CASSGDRLYEQYF TCRBV02-01*01 TCRBJ02-07*03 23/2 6.1415E-07 None 36 CSVRDNYNQPQHF TCRBV29-01*01 TCRBJ01-05*03 23/2 6.1415E-07 None 37 CASSYGGLGSYEQYF TCRBV06-05*01 TCRBJ02-07*03 23/2 6.1415E-07 None 38 CASNRDRGRYEQYF TCRBV06-01*01 TCRBJ02-07*03 20/1 9.1836E-07 B45 39 CASMGGASYEQYF TCRBV27-01*01 TCRBJ02-07*03 20/1 9.1836E-07 A2 40 CASSLGVGPYNEQFF TCRBV07-02*01 TCRBJ02-01*03 20/1 9.1836E-07 B7 41 CASSLGGAGDTQYF TCRBV12 TCRBJ02-03*03 44/15 1.1802E-06 None 42 CATSRGTVSYEQYF TCRBV15-01*01 TCRBJ02-07*03 30/6 1.2303E-06 None 43 CATSDGETQYF TCRBV24 TCRBJ02-05*03 70/36 1.3031E-06 A1,B52 44 CASSEARGGVEKLFF TCRBV06-01*01 TCRBJ01-04*03 22/2 1.3173E-06 None 45 CASSLNRGQETQYF TCRBV14-01*01 TCRBJ02-05*03 19/1 2.0272E-06 None 46 CSVRDNHNQPQHF TCRBV29-01*01 TCRBJ01-05*03 19/1 2.0272E-06 None 47 CASSESGHRNQPQHF TCRBV10-02*01 TCRBJ01-05*03 16/0 2.1105E-06 None 48 CSASPGQGASYGYTF TCRBV20 TCRBJ01-02*03 16/0 2.1105E-06 None 49 CASSEARTRAFF TCRBV06-01*01 TCRBJ01-01*03 16/0 2.1105E-06 None 50 CASRPTGYEQYF TCRBV06-01*01 TCRBJ02-07*03 21/2 2.8099E-06 B39 51 CASSVTGGTDTQYF TCRBV09-01 TCRBJ02-03*03 74/41 2.8214E-06 A1,B8 52 CASSRLAASTDTQYF TCRBV07-03*01 TCRBJ02-03*03 23/3 3.1389E-06 A1,B8 53 CATSDSVTNTGELFF TCRBV24 TCRBJ02-02*03 18/1 4.4551E-06 54 CASSRNRESNQPQHF TCRBV06-05*01 TCRBJ01-05*03 18/1 4.4551E-06 None 55 CASSAQGAYEQYF TCRBV09-01 TCRBJ02-07*03 28/6 4.717E-06 None 56 CASSIQGYSNQPQHF TCRBV05-08*01 TCRBJ01-05*03 15/0 4.8494E-06 B8 57 CASSYNPYSNQPQHF TCRBV06-06 TCRBJ01-05*03 15/0 4.8494E-06 None 58 CASSLGHRDPNTGELFF TCRBV05-01*01 TCRBJ02-02*03 15/0 4.8494E-06 59 CASSTTGGDGYTF TCRBV19-01 TCRBJ01-02*03 26/5 5.5727E-06 None 60 CASSVLAGPTDTQYF TCRBV09-01 TCRBJ02-03*03 26/5 5.5727E-06 Al 61 CASSYRQETQYF TCRBV06-05*01 TCRBJ02-05*03 20/2 5.9592E-06 None 62 CASSSGQVYGYTF TCRBV05-06*01 TCRBJ01-02*03 22/3 6.4736E-06 Al,B8 63 CASGRDTYEQYF TCRBV02-01*01 TCRBJ02-07*03 17/1 9.7453E-06 None 64 CATSDSRTGGQETQYF TCRBV24 TCRBJ02-05*03 17/1 9.7453E-06 B8 65 CASSSPGRSGANVLTF TCRBV28-01*01 TCRBJ02-06*03 17/1 9.7453E-06 None 66 CASSYGGEGYTF TCRBV06-05*01 TCRBJ01-02*03 36/12 1.059E-05 None 67 CASSLAGVDYEQYF TCRBV07-09 TCRBJ02-07*03 25/5 1.0951E-05 B8 68 CASSLQGADTQYF TCRBV07-08*01 TCRBJ02-03*03 14/0 1.112E-05 None 69 CASSLEAENEQFF TCRBV07-02*01 TCRBJ02-01*03 14/0 1.112E-05 All,B8 70 CASSEAPSTSTDTQYF TCRBV02-01*01 TCRBJ02-03*03 14/0 1.112E-05 None 71 CASSLEGQQPQHF TCRBV28-01*01 TCRBJ01-05*03 14/0 1.112E-05 None 72 CASSLGHRDSSYEQYF TCRBV05-01*01 TCRBJ02-07*03 14/0 1.112E-05 A29,B57 73 CASSPSRNTEAFF TCRBV04-03*01 TCRBJ01-01*03 19/2 1.2561E-05 A3,B7 74 CSALGHSNQPQHF TCRBV20 TCRBJ01-05*03 19/2 1.2561E-05 None 75 CASSHRDRNYEQYF TCRBV07-09 TCRBJ02-07*03 19/2 1.2561E-05 A24 76 CASSPPGQGSDTQYF TCRBV18-01*01 TCRBJ02-03*03 28/7 1.3913E-05 None 77 CASSLQGYSNQPQHF TCRBV05-08*01 TCRBJ01-05*03 34/11 1.4711E-05 B8 78 CASSYVRTGGNYGYTF TCRBV06-05*01 TCRBJ01-02*03 31/9 1.5052E-05 None 79 CASSRDRNYGYTF TCRBV06-04 TCRBJ01-02*03 29/8 2.0024E-05 None 80 CASSTGTSGSYEQYF TCRBV06-01*01 TCRBJ02-07*03 16/1 2.1213E-05 B54 81 CASRSDSGANVLTF TCRBV06-04 TCRBJ02-06*03 16/1 2.1213E-05 None 82 CATSRVAGETQYF TCRBV15-01*01 TCRBJ02-05*03 24/5 2.1345E-05 B7 83 CASSEEGIQPQHF TCRBV02-01*01 TCRBJ01-05*03 22/4 2.4814E-05 None 84 CASSLGGPGDTQYF TCRBV12 TCRBJ02-03*03 22/4 2.4814E-05 None 85 CASSLVAAGRETQYF TCRBV05-06*01 TCRBJ02-05*03 13/0 2.5446E-05 B8 86 CASRGQGWDEKLFF TCRBV06-05*01 TCRBJ01-04*03 13/0 2.5446E-05 None 87 CASSLEGQGFGYTF TCRBV05-01*01 TCRBJ01-02*03 13/0 2.5446E-05 None 88 CASRDWDYTDTQYF TCRBV02-01*01 TCRBJ02-03*03 13/0 2.5446E-05 None 89 CASSRSGLAGNTGELFF TCRBV06 TCRBJ02-02*03 13/0 2.5446E-05 None 90 CASSPGQEAGANVLTF TCRBV05-01*01 TCRBJ02-06*03 13/0 2.5446E-05 None 91 CASSLGDRPDTQYF TCRBV11-02*02 TCRBJ02-03*03 13/0 2.5446E-05 None 92 CASSFPGGETQYF TCRBV11-01*01 TCRBJ02-05*03 18/2 2.6306E-05 None 93 CASSLETYGYTF TCRBV05-06*01 TCRBJ01-02*03 18/2 2.6306E-05 B78 94 CASSSGQVQETQYF TCRBV11-02*02 TCRBJ02-05*03 18/2 2.6306E-05 B51 95 CASSFDNYGYTF TCRBV05-04*01 TCRBJ01-02*03 18/2 2.6306E-05 None 96 CASSEGARQPQHF TCRBV10-02*01 TCRBJ01-05*03 18/2 2.6306E-05 None 97 CASSLTGGRNQPQHF TCRBV12 TCRBJ01-05*03 33/11 2.6324E-05 98 CASSLLWDQPQHF TCRBV05-05*01 TCRBJ01-05*03 20/3 2.6954E-05 None 99 CASSLFGTGGNTEAFF TCRBV05-06*01 TCRBJ01-01*03 20/3 2.6954E-05 B44 100 CASSISAGEAFF TCRBV19-01 TCRBJ01-01*03 20/3 2.6954E-05 None 101 CASSPPSGLTDTQYF TCRBV28-01*01 TCRBJ02-03*03 20/3 2.6954E-05 None 102 CASSPLSDTQYF TCRBV07-09 TCRBJ02-03*03 30/9 2.752E-05 None 103 CASSSRGTGELFF TCRBV28-01*01 TCRBJ02-02*03 25/6 3.3358E-05 None 104 CASSYAGDGYTF TCRBV06-05*01 TCRBJ01-02*03 31/10 3.6414E-05 B8 105 CASSDRGNTGELFF TCRBV04-01*01 TCRBJ02-02*03 31/10 3.6414E-05 A68,B41 106 CASSPGGTQYF TCRBV12 TCRBJ02-05*03 28/8 3.6921E-05 None 107 CASSLGDRAYNEQFF TCRBV05-06*01 TCRBJ02-01*03 28/8 3.6921E-05 None 108 CASSLRGSSYNEQFF TCRBV05-08*01 TCRBJ02-01*03 23/5 4.1254E-05 None 109 CASSLTASSYEQYF TCRBV05-01*01 TCRBJ02-07*03 23/5 4.1254E-05 None 110 CSASDHEQYF TCRBV20-01*01 TCRBJ02-07*03 23/5 4.1254E-05 None 111 CASSQGRHTDTQYF TCRBV14-01*01 TCRBJ02-03*03 15/1 4.5934E-05 A68 112 CASSRPGQGNTEAFF TCRBV12 TCRBJ01-01*03 26/7 4.8714E-05 None 113 CASSLVGDGYTF TCRBV07-08*01 TCRBJ01-02*03 36/14 4.9771E-05 None 114 CASSLGAGNQPQHF TCRBV28-01*01 TCRBJ01-05*03 29/9 4.9861E-05 None 115 CASSLTDTGELFF TCRBV11-02*02 TCRBJ02-02*03 29/9 4.9861E-05 None 116 CASSLTGGNSGNTIYF TCRBV07-02*01 TCRBJ01-03*03 19/3 5.4373E-05 B14 117 CAWRGTGNSPLHF TCRBV30-01*01 TCRBJ01-06*03 17/2 5.4712E-05 None 118 CASASANYGYTF TCRBV12 TCRBJ01-02*03 12/0 5.8109E-05 119 CASSLQAGANEQFF TCRBV07-02*01 TCRBJ02-01*03 12/0 5.8109E-05 None 120 CASSEEAGGSGYTF TCRBV06-01*01 TCRBJ01-02*03 12/0 5.8109E-05 None 121 CASRTGESGYTF TCRBV06-05*01 TCRBJ01-02*03 12/0 5.8109E-05 None 122 CASSGLNEQFF TCRBV06-01*01 TCRBJ02-01*03 12/0 5.8109E-05 None 123 CASSWDRDNSPLHF TCRBV25-01*01 TCRBJ01-06*03 12/0 5.8109E-05 B7 124 CASSIRTNYYGYTF TCRBV19-01 TCRBJ01-02*03 12/0 5.8109E-05 None 125 CSARSPEAFF TCRBV20-01*01 TCRBJ01-01*03 12/0 5.8109E-05 None 126 CASSRGTGATDTQYF TCRBV19-01 TCRBJ02-03*03 12/0 5.8109E-05 None 127 CASSPRVSNQPQHF TCRBV12 TCRBJ01-05*03 12/0 5.8109E-05 None 128 CAISESQDRGHEQYF TCRBV10-03*01 TCRBJ02-07*03 12/0 5.8109E-05 None 129 CASSLGRGYEKLFF TCRBV05-06*01 TCRBJ01-04*03 12/0 5.8109E-05 B62 130 CSVEVRGTDTQYF TCRBV29-01*01 TCRBJ02-03*03 12/0 5.8109E-05 None 131 CASRGQGAGELFF TCRBV02-01*01 TCRBJ02-02*03 30/10 6.4924E-05 None 132 CATSREGSGYEQYF TCRBV15-01*01 TCRBJ02-07*03 22/5 7.9023E-05 None 133 CASSLGWTEAFF TCRBV05-01*01 TCRBJ01-01*03 22/5 7.9023E-05 None 134 CASSLGSSSYNEQFF TCRBV07-02*01 TCRBJ02-01*03 35/14 8.4988E-05 None 135 CASSSAGADTQYF TCRBV07-09 TCRBJ02-03*03 14/1 9.8907E-05 B5 136 CASSERKNYGYTF TCRBV06-01*01 TCRBJ01-02*03 14/1 9.8907E-05 None 137 CASRDRDRVNTEAFF TCRBV06-01*01 TCRBJ01-01*03 14/1 9.8907E-05 A3 138 CASSRVGEQFF TCRBV03 TCRBJ02-01*03 14/1 9.8907E-05 None 139 CASSPRWQETQYF TCRBV27-01*01 TCRBJ02-05*03 14/1 9.8907E-05 B14 140 CASTPGDTIYF TCRBV25-01*01 TCRBJ01-03*03 14/1 9.8907E-05 B55
141 CASSENGGNQPQHF TCRBV06-01*01 TCRBJ01-05*03 14/1 9.8907E-05 None
142 CASSYPGETQYF TCRBV06-06 TCRBJ02-05*03 36/15 9.991E-05 None Dimensionality reduction and machine learning:
[00141] CMV burden was calculated for each subject. CMV burden is defined as the fraction of a subject's unique TCRf3s that are significantly CMV-associated, as shown in FIG.
6. The graph in FIG. 6 shows the distribution of CMV scores (i.e., the proportion of each subject's TCRf3 repertoire that matches our list of 142 CMV-associated TCRf3 sequences) among CMV+ and CMV- subjects.
[00142] This enabled fast training of a one-dimensional logistic regression classifier of CMV status. Exhaustive leave-one-out cross validation (including re-computation of CMV-associated clones) was performed, indicating high accuracy across a broad range of p-value thresholds, with performance degrading at a high FDR.
[00141] CMV burden was calculated for each subject. CMV burden is defined as the fraction of a subject's unique TCRf3s that are significantly CMV-associated, as shown in FIG.
6. The graph in FIG. 6 shows the distribution of CMV scores (i.e., the proportion of each subject's TCRf3 repertoire that matches our list of 142 CMV-associated TCRf3 sequences) among CMV+ and CMV- subjects.
[00142] This enabled fast training of a one-dimensional logistic regression classifier of CMV status. Exhaustive leave-one-out cross validation (including re-computation of CMV-associated clones) was performed, indicating high accuracy across a broad range of p-value thresholds, with performance degrading at a high FDR.
[00143] Each subject was removed from the dataset in turn, and the list of significantly CMV-associated TCRPs was re-calculated using the remaining subjects with known CMV
status. Each subject's CMV burden was then calculated as described above. A
one-dimensional logistic regression classifier was then trained on CMV burden vs.
CMV status, and the CMV burden of the subject which was held out was then calculated and input into the logistic regression. The output of this logistic regression model was the probability that the subject of unknown status was CMV+, which was then compared to the (known but not leveraged) status of the held-out subject to determine the accuracy of the hold-out classification. Classification of a subject with genuinely unknown CMV status proceeds likewise: sequence TCRP, calculate a CMV burden (as proportion of unique TCRPs from this subject present on the list of CMV-associated TCRf3s), and input this CMV
burden into a logistic regression model trained on subjects of known status, with an estimated probability of CMV-positivity as output.
status. Each subject's CMV burden was then calculated as described above. A
one-dimensional logistic regression classifier was then trained on CMV burden vs.
CMV status, and the CMV burden of the subject which was held out was then calculated and input into the logistic regression. The output of this logistic regression model was the probability that the subject of unknown status was CMV+, which was then compared to the (known but not leveraged) status of the held-out subject to determine the accuracy of the hold-out classification. Classification of a subject with genuinely unknown CMV status proceeds likewise: sequence TCRP, calculate a CMV burden (as proportion of unique TCRPs from this subject present on the list of CMV-associated TCRf3s), and input this CMV
burden into a logistic regression model trained on subjects of known status, with an estimated probability of CMV-positivity as output.
[00144] As shown in FIG. 7, the cross-validation method left out one of the initial 640 samples, then the database of CMV-specific TCRI3 sequences and associated statistics was retrained in an unsupervised fashion to eliminate bias, and the CMV serostatus of the left-out sample was classified.
[00145] The cross-validated results are shown in FIGs. 8 and 9 (shown for all subjects and cross-validated (CV) subjects). The accuracy of the model was shown to be at 89%, corresponding to a diagnostic odds ratio of ¨66. Assessing model performance across the range of specificity vs. sensitivity using an ROC curve, an area under the ROC
curve of 93%
was achieved.
curve of 93%
was achieved.
[00146] FIG. 8 (top graph) shows data for the classification performance of all and the cross-validation (CV) datasets for each p-value threshold, measured as the area under the ROC curve (AUROC). The number above each set of data points corresponds to the number of CMV-associated TCRf3 identified at that p-value threshold, and the rectangle indicates the dataset selected for downstream analysis (p-value = 10-4).
[00147] FIG. 8 (bottom graph) also shows a false discovery rate (FDR) estimated for each p-value threshold used in the identification of significantly CMV-associated TCRf3 sequences, using permutations of CMV status. The best performance is seen at a p-value of 104, which corresponds to an estimated FDR of ¨20%, resulting in the identification of a set of 142 TCRf3 sequences that were significantly associated with positive CMV
status (listed in Table 1). Using these conditions resulted in a good separation between the CMV+ and CMV-subjects in the cohort as measured by CMV score (See FIG. 6).
status (listed in Table 1). Using these conditions resulted in a good separation between the CMV+ and CMV-subjects in the cohort as measured by CMV score (See FIG. 6).
[00148] FIG. 9 shows the ROC curves for both the all and the cross-validation datasets.
The AUROC for the full dataset is 0.98, indicating that our approach resulted in an excellent classifier for CMV status. At the point of highest discriminating power, an accuracy of 0.89 and a diagnostic odds ratio of 66 in the cross-validation dataset was observed (achieved when classifying 86% of true positives correctly with a false positive rate of 8%).
Taken together, these data suggest that that presence of public T cell responses to CMV is highly correlated with CMV positive status.
The AUROC for the full dataset is 0.98, indicating that our approach resulted in an excellent classifier for CMV status. At the point of highest discriminating power, an accuracy of 0.89 and a diagnostic odds ratio of 66 in the cross-validation dataset was observed (achieved when classifying 86% of true positives correctly with a false positive rate of 8%).
Taken together, these data suggest that that presence of public T cell responses to CMV is highly correlated with CMV positive status.
[00149] Given that T cells recognize their cognate antigens in the context of MHC
molecules expressed by antigen presenting cells, it was tested whether the HLA-restriction of our CMV-associated TCRf3 sequences could be identified. A Fisher's exact test was performed on each CMV-associated TCRf3 sequence to determine if its presence was significantly associated with any of the HLA alleles observed in the cohort.
The association of 57 out of 142 CMV-associated TCRf3 sequences could be confidently assigned with at least one HLA allele, with a p-value cutoff of 1x103. Full results are presented in Table 1 (above) and FIG. 10.
molecules expressed by antigen presenting cells, it was tested whether the HLA-restriction of our CMV-associated TCRf3 sequences could be identified. A Fisher's exact test was performed on each CMV-associated TCRf3 sequence to determine if its presence was significantly associated with any of the HLA alleles observed in the cohort.
The association of 57 out of 142 CMV-associated TCRf3 sequences could be confidently assigned with at least one HLA allele, with a p-value cutoff of 1x103. Full results are presented in Table 1 (above) and FIG. 10.
[00150] FIGs. 10A and 10B show HLA-restriction of CMV-associated TCRf3 sequences.
FIG. 10A shows the distribution of HLA-A alleles in this cohort. FIG. 10B
shows the distribution of HLA-B alleles in this cohort. Each of the 142 CMV-associated TCRf3 sequences identified at p < 1x104 was tested for significant association with each HLA allele, with a p-value threshold of 1x103. Of these 142 CMV-associated TCRf3 sequences, 57 were significantly associated with an HLA-A and/or an HLA-B allele, and no sequences were significantly associated with more than a single allele from each locus.
Colored boxes and amino acid sequences indicate the 5 TCRf3 sequences identified in our study that had been previously identified. In 4 of the cases, the correct HLA association was recapitulated, and in the 5th case, there was no statistically-significant HLA association.
FIG. 10A shows the distribution of HLA-A alleles in this cohort. FIG. 10B
shows the distribution of HLA-B alleles in this cohort. Each of the 142 CMV-associated TCRf3 sequences identified at p < 1x104 was tested for significant association with each HLA allele, with a p-value threshold of 1x103. Of these 142 CMV-associated TCRf3 sequences, 57 were significantly associated with an HLA-A and/or an HLA-B allele, and no sequences were significantly associated with more than a single allele from each locus.
Colored boxes and amino acid sequences indicate the 5 TCRf3 sequences identified in our study that had been previously identified. In 4 of the cases, the correct HLA association was recapitulated, and in the 5th case, there was no statistically-significant HLA association.
[00151] A literature search was performed to identify previously reported CMV-reactive TCRf3 sequences, and 595 unique TCRf3 sequences were selected that had been identified by at least one previously published study 3, 10, 12-27. Of these, 30 unique TCRf3 sequences had previously been previously classified as public, or were reported in multiple studies. It was determined that 5 of these 30 public TCRf3 sequences unique were contained in the set of 142 CMV-associated TCRf3 sequences. Furthermore, an HLA association was identified for 4 of these 5 TCRf3 sequences, and in all 4 cases our HLA association agreed with the literature.
FIGs. 11A, 11B, and 12 provide more complete information on the prevalence in our dataset of previously identified CMV-reactive TCRf3 sequences, as described below.
FIGs. 11A, 11B, and 12 provide more complete information on the prevalence in our dataset of previously identified CMV-reactive TCRf3 sequences, as described below.
[00152] FIGs. 11A and 11B show the incidence of previously reported CMV-reactive TCRf3 sequences in this cohort. After a literature search, 565 'private' CMV-reactive TCRf3 sequences (reported in one individual from a single study) and 30 'public' CMV-reactive TCRf3 sequences (reported in multiple studies or in multiple individuals within a single study) were identified. FIG. 11A shows the incidence of each such TCRf3 sequence in the cohort of 640 subjects plotted along the horizontal axis by decreasing total incidence, with the incidence in CMV+ subjects above the horizontal and the incidence in CMV-subjects below the horizontal. Many previously-reported sequences were observed in the dataset, but most were seen in roughly equal number of CMV+ and CMV- subjects, which could be explained by receptor sequences with exceptionally high frequency in the naive repertoire, or could reflect cross-reactive receptors that bind to CMV antigens but also other common antigens.
FIG. 11B shows a histogram of incidence of these TCRf3 sequences in the cohort of 640 subjects plotted for each group of sequences. Most previously-reported CMV-reactive TCRf3 sequences were found in our dataset at appreciable levels, though only a handful were found disproportionately in CMV+ subjects at p < 1x104. 'Public' TRf3 sequences reported in the literature were considerably more common in our cohort.
FIG. 11B shows a histogram of incidence of these TCRf3 sequences in the cohort of 640 subjects plotted for each group of sequences. Most previously-reported CMV-reactive TCRf3 sequences were found in our dataset at appreciable levels, though only a handful were found disproportionately in CMV+ subjects at p < 1x104. 'Public' TRf3 sequences reported in the literature were considerably more common in our cohort.
[00153] FIG. 12 shows the concordance of TCRB sequences in the cohort as compared to those in the literature. Of the 142 TCRf3 sequences significantly associated with CMV status in the cohort with a p-value of less than 1x104, five TCRf3 sequences (defined by matching CDR3 amino acid sequence) have been previously reported in the literature as public clones.
FIG. 12 provides the CDR3 amino acid sequence, V and J genes, and HLA
association for each sequence and compares it to previous reports. As expected, the two sequences previously identified as public were seen in more subjects in our cohort than the three sequences previously reported only once. Concordance of V gene, J gene, and HLA allele association with those reported in the literature was very good. Five out of nine total comparisons had the same V gene, and the other 4 comparisons had the same V
gene subfamily. There were concordant J genes in 9 out of 9 comparisons, and an identical HLA
association in 8 out of 9 comparisons, with one sequence not significantly HLA-associated in this study.
FIG. 12 provides the CDR3 amino acid sequence, V and J genes, and HLA
association for each sequence and compares it to previous reports. As expected, the two sequences previously identified as public were seen in more subjects in our cohort than the three sequences previously reported only once. Concordance of V gene, J gene, and HLA allele association with those reported in the literature was very good. Five out of nine total comparisons had the same V gene, and the other 4 comparisons had the same V
gene subfamily. There were concordant J genes in 9 out of 9 comparisons, and an identical HLA
association in 8 out of 9 comparisons, with one sequence not significantly HLA-associated in this study.
[00154] Thus, immunosequencing can be used to determine CMV-specific T cell sequences in a large cohort and predict CMV status in new subjects with high accuracy. In summary, it has been demonstrated that information gleaned from rearranged T
cell receptors can be used to infer disease status based on the presence of public T cell responses; the only requirement is a large sample of pathogen-positive and -negative samples with which to identify these public T cell responses. Because high-throughput sequencing of T cell receptors captures all T cell responses equally, and these store immunological memory to all pathogens in a common format, reading T cell memory by looking for known public responses can be a viable strategy for simultaneously diagnosing a wide range of infectious agents using a single peripheral blood sample and a simple, unified assay.
These methods can be applied needed to acute infections, given that T cell memory persists for years, accounting for the fact that it is not known how public clones decay with time after an acute infection. Accordingly, this method can be used to assess multiple infections simultaneously, such as HPV, EBV, CMV and others. The method can also be used to predict or diagnose a non-infectious (e.g. autoimmune) disease.
cell receptors can be used to infer disease status based on the presence of public T cell responses; the only requirement is a large sample of pathogen-positive and -negative samples with which to identify these public T cell responses. Because high-throughput sequencing of T cell receptors captures all T cell responses equally, and these store immunological memory to all pathogens in a common format, reading T cell memory by looking for known public responses can be a viable strategy for simultaneously diagnosing a wide range of infectious agents using a single peripheral blood sample and a simple, unified assay.
These methods can be applied needed to acute infections, given that T cell memory persists for years, accounting for the fact that it is not known how public clones decay with time after an acute infection. Accordingly, this method can be used to assess multiple infections simultaneously, such as HPV, EBV, CMV and others. The method can also be used to predict or diagnose a non-infectious (e.g. autoimmune) disease.
[00155] The methods of the invention can also be used to quantify disease burden (e.g.
CMV reactivation post-transplant).
Example 2: Identification of HLA Type
CMV reactivation post-transplant).
Example 2: Identification of HLA Type
[00156] The first aim of the study was to create a comprehensive catalog of T
cell receptor sequences and associated HLA types.
HLA Typing:
cell receptor sequences and associated HLA types.
HLA Typing:
[00157] 640 subjects were phenotyped for HLA type according to standard protocol by the Fred Hutchinson Cancer Research Center.
Immune Repertoire Sequencing:
Immune Repertoire Sequencing:
[00158] Next, the TCRB sequences for each subject was determined by amplification and sequencing.
[00159] Genomic DNA was extracted from peripheral blood samples using the Qiagen DNeasy Blood extraction Kit (Qiagen, Gaithersburg, MD, USA). The CDR3 regions of rearranged TCRf3 genes were sequenced; the TCRf3 CDR3 region was defined according to the IMGT collaboration28' 20. TCRf3 CDR3 regions were amplified and sequenced using methods described above and in previously described protocols5' 29. The multiplexed PCR
method used a mixture of 60 forward primers specific to TCR VB gene segments and 13 reverse primers specific to TCR .10 gene segments. The resulting amplicons were sequenced using the methods described above. Reads of 87 bp were obtained using the Illumina HiSeq System. Raw HiSeq sequence data were preprocessed to remove errors in the primary sequence of each read, and to compress the data. A nearest neighbor algorithm was used to collapse the data into unique sequences by merging closely related sequences, to remove both PCR and sequencing errors.
method used a mixture of 60 forward primers specific to TCR VB gene segments and 13 reverse primers specific to TCR .10 gene segments. The resulting amplicons were sequenced using the methods described above. Reads of 87 bp were obtained using the Illumina HiSeq System. Raw HiSeq sequence data were preprocessed to remove errors in the primary sequence of each read, and to compress the data. A nearest neighbor algorithm was used to collapse the data into unique sequences by merging closely related sequences, to remove both PCR and sequencing errors.
[00160] In order to ensure adequate coverage of each T cell rearrangement, 8-10x sequence coverage, or ¨6-10 million sequencing reads per sample, was generated using approximately eight full sequencing runs. All sequencing reads were processed using a standardized bioinformatics pipeline to 1) demultiplex reads to specific samples, 2) eliminate low quality sequence and remove potential contaminants, 3) align and identify specific TCRB
V and J gene segments and CDR3 regions, 4) cluster highly similar sequences to account for PCR and sequencing errors, 5) normalize data to remove PCR amplification bias, 6) estimate total T-cell input, and 7) generate TCRB unique sequence counts and distributions.
V and J gene segments and CDR3 regions, 4) cluster highly similar sequences to account for PCR and sequencing errors, 5) normalize data to remove PCR amplification bias, 6) estimate total T-cell input, and 7) generate TCRB unique sequence counts and distributions.
[00161] Here, approximately 250,000 rearranged T cell receptor genes were sequenced from peripheral blood of each subject. For each subject, approximately 105 TCRB sequences were obtained, and the abundances of each unique sequence were quantified.
Feature Selection: To define association between an HLA allele and TCRB
sequence(s)
Feature Selection: To define association between an HLA allele and TCRB
sequence(s)
[00162] For each subject, the total number of unique TCRB sequences and the frequency of each unique TCRB sequence are determined. For a unique TCRB sequence, it is determined how many subjects who are positive for an HLA allele have the TCRB
sequence and how many subjects who are negative for the HLA allele have the TCRB
sequence. In addition, it can be determined the number of subjects who are positive for an HLA allele and negative for the TCRB sequence and the number of subjects who are negative for the HLA
allele and positive for the TCRB sequence. The table below shows categorization of subjects by the presence or absence of a TCRB sequence and the presence or absence of an HLA
allele, HLA-A2.
HLA-A2+ HLA-A2-TCRf3 sequence i present TCRf3 sequence i not present N+ ¨ n, N_
sequence and how many subjects who are negative for the HLA allele have the TCRB
sequence. In addition, it can be determined the number of subjects who are positive for an HLA allele and negative for the TCRB sequence and the number of subjects who are negative for the HLA
allele and positive for the TCRB sequence. The table below shows categorization of subjects by the presence or absence of a TCRB sequence and the presence or absence of an HLA
allele, HLA-A2.
HLA-A2+ HLA-A2-TCRf3 sequence i present TCRf3 sequence i not present N+ ¨ n, N_
[00163] FIG. 14 shows a list of exemplary unique TCRB sequences and the number of subjects who are positive or negative for an HLA-A2 allele that have a particular TCRB
sequence.
sequence.
[00164] A p-value is determined for the association of each TCRB sequence with an HLA
status using a Fisher exact test (two-tailed). The p value for association of each TCR with allele status using a Fisher exact test (two tailed) is calculated as follows:
(n++n_ (N++N_ -n+-n_) Pr(table) = n+
(N++N_ N
status using a Fisher exact test (two-tailed). The p value for association of each TCR with allele status using a Fisher exact test (two tailed) is calculated as follows:
(n++n_ (N++N_ -n+-n_) Pr(table) = n+
(N++N_ N
[00165] FIG. 14 shows exemplary p values that were calculated for the association of a particular TCRB sequence with an HLA type (HLA-A2).
[00166] As shown in FIG. 15, a p value is selected as a cutoff for identifying a set of "Feature TCRs" from the entire list of possible TCR sequences. Defining a p-value threshold and permuting the allele status across individuals provides an estimate of false discovery rate.
This is performed for each HLA allele, resulting in a set of allele-associated TCRf3 sequences for each HLA allele. In FIG. 15, a p value cutoff of p<10-4 and an FDR of 0.1 is used to identify 288 TCRB sequences that are positively associated with HLA-A2. For each of the allele-associated TCRB sequences, the frequency of the sequence is also determined in each subject.
This is performed for each HLA allele, resulting in a set of allele-associated TCRf3 sequences for each HLA allele. In FIG. 15, a p value cutoff of p<10-4 and an FDR of 0.1 is used to identify 288 TCRB sequences that are positively associated with HLA-A2. For each of the allele-associated TCRB sequences, the frequency of the sequence is also determined in each subject.
[00167] As shown in FIG. 15, a p value is selected as a cutoff for identifying a set of "Feature TCRS" from the entire list of possible TCR sequences. Defining a p-value threshold and permuting the allele status across individuals provides an estimate of false discovery rate. This is performed for each HLA allele, resulting in a set of allele-associated TCRf3 sequences for each HLA allele. In FIG. 15, a p value cutoff of p<10-4 and an FDR of 0.1 is used to identify 288 TCRB sequences that are positively associated with HLA-A2. For each of the allele-associated TCRB sequences, the frequency of the sequence is also determined in each subject.
[00168] The false discovery rate (FDR) was determined by permutation of allele status.
See Storey et al. Statistical significance for genomewide studies. PNAS, 100(6), pp. 9440-9445.
See Storey et al. Statistical significance for genomewide studies. PNAS, 100(6), pp. 9440-9445.
[00169] The feature selection step is followed by a machine learning process.
As shown in FIG. 16, for each HLA allele, a logistic regression model is trained using the set of feature vectors over all subjects, along with the known status for presence of that allele.
As shown in FIG. 16, for each HLA allele, a logistic regression model is trained using the set of feature vectors over all subjects, along with the known status for presence of that allele.
[00170] As shown in FIG. 17, an exhaustive leave-one-out cross validation is performed where one subject is removed from the analysis, and the HLA status of the subject is inferred based on feature selection and training from only the remaining subjects. The result is a set of classifiers (one for each HLA allele) that estimate the probability of positive status for each HLA allele, taking as input the feature vector for each allele.
[00171] Equipped with these classifiers, the HLA type of a new subject can be assessed by: 1. immune repertoire sequencing, 2. computing feature vectors for each allele, and 3.
defining a probability threshold for positive status that will be applied to the output from each classifier. In this manner, an HLA status can be inferred for a new subject with an unknown HLA status.
defining a probability threshold for positive status that will be applied to the output from each classifier. In this manner, an HLA status can be inferred for a new subject with an unknown HLA status.
[00172] FIG. 18 shows the results of a cross validation study. The leave-one out cross validation was performed for 78 rounds. For the HLA-A2 allele, the method described above resulted in identification of 288 feature TCRs and an accurate prediction of 41 subjects out of 43 as HLA-A2 positive. As described above, the HLA type of the subjects was known prior to performing the cross validation. Only 2/43 subjects were false-positives for HLA-A2.
Thus, the method was 96% accurate in predicting HLA-A2 presence in a subject based on the subject's TCRB sequence profile.
Thus, the method was 96% accurate in predicting HLA-A2 presence in a subject based on the subject's TCRB sequence profile.
[00173] FIG. 18 also shows the results for cross-validation of the HLA-A24 allele, which also had a 96% accuracy. 10 out of 13 subjects were accurately predicted to possess the HLA-A24 allele, and 65 out of 65 subjects were accurately predicted to as not have the HLA-A24 allele based on the identified TCRB features.
Example 3: HLA Study
Example 3: HLA Study
[00174] In another example, the public T-cell response to cytomegalovirus (CMV) was investigated by sequencing rearranged T cell receptors (TCRs) in 650 subjects (294 with and 356 without CMV). The concordance between ¨90 million unique TCRs and CMV
serostatus was assessed, focusing on identification of significant associations. In this study, 157 CMV-associated TCRs were identified at p < 10-3, FDR 0.15) Training a binary classifier on these features, it was predicted that CMV serostatus in a leave-one-out cross-validation procedure had a diagnostic odds ratio of 44. The classifier was also tested on a second independent cohort of 120 subjects with known CMV serostatus, yielding a diagnostic odds ratio of 49.
serostatus was assessed, focusing on identification of significant associations. In this study, 157 CMV-associated TCRs were identified at p < 10-3, FDR 0.15) Training a binary classifier on these features, it was predicted that CMV serostatus in a leave-one-out cross-validation procedure had a diagnostic odds ratio of 44. The classifier was also tested on a second independent cohort of 120 subjects with known CMV serostatus, yielding a diagnostic odds ratio of 49.
[00175] Next, the HLA-restriction of each CMV-associated TCR was investigated by assessing the over-representation of particular HLA types among the subjects that carry each CMV-associated TCR. Of 157 CMV-associated TCRs, 61 were HLA-associated at p <
10-3.
None of these were significantly associated with multiple HLA-A or HLA-B
alleles.
10-3.
None of these were significantly associated with multiple HLA-A or HLA-B
alleles.
[00176] There was substantial concordance between our data and previously published CMV- and HLA-associated TCRs. Most previously-reported public CMV-specific TCRs were seen in the data, although only 5/157 CMV-associated TCRs identified in this study have been previously reported. Of these, 4 were significantly HLA-restricted, and all four confirmed previous findings.
[00177] In addition, the association of TCRs with each HLA-A allele present in the cohort was investigated, with significant (p < 10-4) results for many higher frequency alleles. Binary classifier training resulted in high accuracy (-96%) prediction for these alleles, indicating that HLA type can be inferred from immunosequencing data.
[00178] In summary, this study demonstrated the validity of association studies using immunosequencing for detection and HLA-association of public T-cell responses to infection, and showed that assessing the presence of associated T-cell responses can serve as a powerful diagnostic classifier.
[00179] While the invention has been particularly shown and described with reference to a preferred embodiment and various alternate embodiments, it will be understood by persons skilled in the relevant art that various changes in form and details can be made therein without departing from the spirit and scope of the invention.
[00180] All references, issued patents and patent applications cited within the body of the instant specification are hereby incorporated by reference in their entirety, for all purposes.
REFERENCES
1. Neller, M.A., Burrows, J.M., Rist, M.J., Miles, J.J. St Burrows, S.R.
High frequency of herpesvirus-specific clonotypes in the human T cell repertoire can remain stable over decades with minimal turnover. J Virol 87, 697-700 (2013).
2. Li, H., Ye, C., Ji, G. St Han, J. Determinants of public T cell responses. Cell Res 22, 33-42 (2012).
3. Venturi, V., Price, D.A., Douek, D.C. St Davenport, M.P. The molecular basis for public T-cell responses? Nat Rev Immunol 8, 231-238 (2008).
4. Arstila, T.P. et al. A direct estimate of the human alphabeta T cell receptor diversity. Science 286, 958-961 (1999).
5. Robins, H.S. et al. Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells. Blood 114, 4099-4107 (2009).
6. Cabaniols, J.P., Fazilleau, N., Casrouge, A., Kourilsky, P. St Kanellopoulos, J.M. Most alpha/beta T cell receptor diversity is due to terminal deoxynucleotidyl transferase. J Exp Med 194, 1385-1390 (2001).
7. Davis, M.M. St Bjorkman, P.J. T-cell antigen receptor genes and T-cell recognition.
Nature 334 (1988).
8. Robins, H.S. et al. Overlap and effective size of the human CD8+ T cell receptor repertoire. Sci Transl Med 2, 47ra64 (2010).
9. Li, H. et al. Recombinatorial biases and convergent recombination determine interindividual TCRbeta sharing in murine thymocytes. J Immuno1189, 2404-2413 (2012).
10. Venturi, V. et al. TCR beta-chain sharing in human CD8+ T cell responses to cytomegalovirus and EBV. J Immuno1181, 7853-7862 (2008).
11. Gandhi, M.K. St Khanna, R. Human cytomegalovirus: clinical aspects, immune regulation, and emerging treatments. Lancet Infect Dis 4, 725-738 (2004).
12. Babel, N. et al. Clonotype analysis of cytomegalovirus-specific cytotoxic T
lymphocytes. J Am Soc Nephrol 20, 344-352 (2009).
13. lancu, E.M. et al. Clonotype selection and composition of human CD8 T
cells specific for persistent herpes viruses varies with differentiation but is stable over time./ Immuno1183, 319-331 (2009).
14. Khan, N., Cobbold, M., Keenan, R. St Moss, P.A. Comparative analysis of CD8+ T
cell responses against human cytomegalovirus proteins pp65 and immediate early 1 shows similarities in precursor frequency, oligoclonality, and phenotype.
J Infect Dis 185, 1025-1034 (2002).
15. Klarenbeek, P.L. et al. Deep sequencing of antiviral T-cell responses to HCMV and EBV in humans reveals a stable repertoire that is maintained for many years.
PLoS Pathog 8, e1002889 (2012).
16. Klinger, M. et al. Combining next-generation sequencing and immune assays: a novel method for identification of antigen-specific T cells. PLoS One 8, e74231 (2013).
17. Koning, D. et al. In vitro expansion of antigen-specific CD8(+) T cells distorts the T-cell repertoire. J Immunol Methods 405, 199-203 (2014).
18. Miconnet, I. et al. Large TCR diversity of virus-specific CD8 T cells provides the mechanistic basis for massive TCR renewal after antigen exposure./ Immunol 186, 7039-7049 (2011).
19. Price, D.A. et al. Avidity for antigen shapes clonal dominance in CD8+
T cell populations specific for persistent DNA viruses. J Exp Med 202, 1349-1361 (2005).
20. Retiere, C. et al. Generation of cytomegalovirus-specific human T-lymphocyte clones by using autologous B-Iymphoblastoid cells with stable expression of pp65 orlE1 proteins: a tool to study the fine specificity of the antiviral response.
J Virol 74, 3948-3952 (2000).
21. Scheinberg, P. et al. The transfer of adaptive immunity to CMV during hematopoietic stem cell transplantation is dependent on the specificity and phenotype of CMV-specific T cells in the donor. Blood 114, 5071-5080 (2009).
22. Schub, A., Schuster, I.G., Hammerschmidt, W. St Moosmann, A. CMV-specific TCR-transgenic T cells for immunotherapy. J Immunol 183, 6819-6830 (2009).
23. Trautmann, L. et al. Selection of T cell clones expressing high-affinity public TCRs within Human cytomegalovirus-specific CD8 T cell responses. J Immunol 175, 6123-6132 (2005).
24. Wynn, K.K. et al. Impact of clonal competition for peptide-MHC
complexes on the CD8+ T-cell repertoire selection in a persistent viral infection. Blood 111, 4292 (2008).
25. Peggs, K. et al. Characterization of human cytomegalovirus peptide-specific CD8(+) T-cell repertoire diversity following in vitro restimulation by antigen-pulsed dendritic cells. Blood 99, 213-223 (2002).
26. Weekes, M.P., Wills, M.R., Mynard, K., Carmichael, A.J. St Sissons, J.G. The memory cytotoxic T-lymphocyte (CTL) response to human cytomegalovirus infection contains individual peptide-specific CTL clones that have undergone extensive expansion in vivo. J Virol 73, 2099-2108 (1999).
27. Lim, A. et al. Frequent contribution of T cell clonotypes with public TCR features to the chronic response against a dominant EBV-derived epitope: application to direct detection of their molecular imprint on the human peripheral T cell repertoire./ Immuno1165, 2001-2011 (2000).
28. Yousfi Monod, M., Giudicelli, V., Chaume, D. St Lefranc, M.P.
IMGT/JunctionAnalysis: the first tool for the analysis of the immunoglobulin and T cell receptor complex V-J and V-D-J JUNCTIONs. Bioinformatics 20 Suppl 1, i379-385 (2004).
29. Carlson, C.S. et al. Using synthetic templates to design an unbiased multiplex PCR
assay. Nat Commun 4, 2680 (2013).
30. Storey, J.D. St Tibshirani, R. Statistical significance for genomewide studies. Proc Natl Acad Sci USA 100, 9440-9445 (2003).
REFERENCES
1. Neller, M.A., Burrows, J.M., Rist, M.J., Miles, J.J. St Burrows, S.R.
High frequency of herpesvirus-specific clonotypes in the human T cell repertoire can remain stable over decades with minimal turnover. J Virol 87, 697-700 (2013).
2. Li, H., Ye, C., Ji, G. St Han, J. Determinants of public T cell responses. Cell Res 22, 33-42 (2012).
3. Venturi, V., Price, D.A., Douek, D.C. St Davenport, M.P. The molecular basis for public T-cell responses? Nat Rev Immunol 8, 231-238 (2008).
4. Arstila, T.P. et al. A direct estimate of the human alphabeta T cell receptor diversity. Science 286, 958-961 (1999).
5. Robins, H.S. et al. Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells. Blood 114, 4099-4107 (2009).
6. Cabaniols, J.P., Fazilleau, N., Casrouge, A., Kourilsky, P. St Kanellopoulos, J.M. Most alpha/beta T cell receptor diversity is due to terminal deoxynucleotidyl transferase. J Exp Med 194, 1385-1390 (2001).
7. Davis, M.M. St Bjorkman, P.J. T-cell antigen receptor genes and T-cell recognition.
Nature 334 (1988).
8. Robins, H.S. et al. Overlap and effective size of the human CD8+ T cell receptor repertoire. Sci Transl Med 2, 47ra64 (2010).
9. Li, H. et al. Recombinatorial biases and convergent recombination determine interindividual TCRbeta sharing in murine thymocytes. J Immuno1189, 2404-2413 (2012).
10. Venturi, V. et al. TCR beta-chain sharing in human CD8+ T cell responses to cytomegalovirus and EBV. J Immuno1181, 7853-7862 (2008).
11. Gandhi, M.K. St Khanna, R. Human cytomegalovirus: clinical aspects, immune regulation, and emerging treatments. Lancet Infect Dis 4, 725-738 (2004).
12. Babel, N. et al. Clonotype analysis of cytomegalovirus-specific cytotoxic T
lymphocytes. J Am Soc Nephrol 20, 344-352 (2009).
13. lancu, E.M. et al. Clonotype selection and composition of human CD8 T
cells specific for persistent herpes viruses varies with differentiation but is stable over time./ Immuno1183, 319-331 (2009).
14. Khan, N., Cobbold, M., Keenan, R. St Moss, P.A. Comparative analysis of CD8+ T
cell responses against human cytomegalovirus proteins pp65 and immediate early 1 shows similarities in precursor frequency, oligoclonality, and phenotype.
J Infect Dis 185, 1025-1034 (2002).
15. Klarenbeek, P.L. et al. Deep sequencing of antiviral T-cell responses to HCMV and EBV in humans reveals a stable repertoire that is maintained for many years.
PLoS Pathog 8, e1002889 (2012).
16. Klinger, M. et al. Combining next-generation sequencing and immune assays: a novel method for identification of antigen-specific T cells. PLoS One 8, e74231 (2013).
17. Koning, D. et al. In vitro expansion of antigen-specific CD8(+) T cells distorts the T-cell repertoire. J Immunol Methods 405, 199-203 (2014).
18. Miconnet, I. et al. Large TCR diversity of virus-specific CD8 T cells provides the mechanistic basis for massive TCR renewal after antigen exposure./ Immunol 186, 7039-7049 (2011).
19. Price, D.A. et al. Avidity for antigen shapes clonal dominance in CD8+
T cell populations specific for persistent DNA viruses. J Exp Med 202, 1349-1361 (2005).
20. Retiere, C. et al. Generation of cytomegalovirus-specific human T-lymphocyte clones by using autologous B-Iymphoblastoid cells with stable expression of pp65 orlE1 proteins: a tool to study the fine specificity of the antiviral response.
J Virol 74, 3948-3952 (2000).
21. Scheinberg, P. et al. The transfer of adaptive immunity to CMV during hematopoietic stem cell transplantation is dependent on the specificity and phenotype of CMV-specific T cells in the donor. Blood 114, 5071-5080 (2009).
22. Schub, A., Schuster, I.G., Hammerschmidt, W. St Moosmann, A. CMV-specific TCR-transgenic T cells for immunotherapy. J Immunol 183, 6819-6830 (2009).
23. Trautmann, L. et al. Selection of T cell clones expressing high-affinity public TCRs within Human cytomegalovirus-specific CD8 T cell responses. J Immunol 175, 6123-6132 (2005).
24. Wynn, K.K. et al. Impact of clonal competition for peptide-MHC
complexes on the CD8+ T-cell repertoire selection in a persistent viral infection. Blood 111, 4292 (2008).
25. Peggs, K. et al. Characterization of human cytomegalovirus peptide-specific CD8(+) T-cell repertoire diversity following in vitro restimulation by antigen-pulsed dendritic cells. Blood 99, 213-223 (2002).
26. Weekes, M.P., Wills, M.R., Mynard, K., Carmichael, A.J. St Sissons, J.G. The memory cytotoxic T-lymphocyte (CTL) response to human cytomegalovirus infection contains individual peptide-specific CTL clones that have undergone extensive expansion in vivo. J Virol 73, 2099-2108 (1999).
27. Lim, A. et al. Frequent contribution of T cell clonotypes with public TCR features to the chronic response against a dominant EBV-derived epitope: application to direct detection of their molecular imprint on the human peripheral T cell repertoire./ Immuno1165, 2001-2011 (2000).
28. Yousfi Monod, M., Giudicelli, V., Chaume, D. St Lefranc, M.P.
IMGT/JunctionAnalysis: the first tool for the analysis of the immunoglobulin and T cell receptor complex V-J and V-D-J JUNCTIONs. Bioinformatics 20 Suppl 1, i379-385 (2004).
29. Carlson, C.S. et al. Using synthetic templates to design an unbiased multiplex PCR
assay. Nat Commun 4, 2680 (2013).
30. Storey, J.D. St Tibshirani, R. Statistical significance for genomewide studies. Proc Natl Acad Sci USA 100, 9440-9445 (2003).
Claims (24)
1. A method for predicting the presence or absence of a cytomegalovirus (CMV) infection in a subject of unknown infection status, the method comprising:
a) performing amplification and high throughput sequencing of genomic DNA
obtained from a sample comprising T cells obtained from the subject to determine a TCR
profile comprising unique TCR CDR3 amino acid sequences;
b) comparing the TCR profile with a database of previously identified diagnostic public T cell receptor sequences that are statistically significantly associated with CMV
infection;
c) generating a CMV burden score for the subject, wherein the CMV burden score is the proportion of unique TCR sequences in the profile of the subject that match the public TCR sequences in the database;
d) inputting the calculated CMV burden score from c) into a logistic regression model, wherein the logistic regression models compares CMV burden and CMV
infection status from a plurality of subjects of known CMV infection status; and e) determining an estimated probability of CMV infection status of the subject as the output of logistic regression model.
a) performing amplification and high throughput sequencing of genomic DNA
obtained from a sample comprising T cells obtained from the subject to determine a TCR
profile comprising unique TCR CDR3 amino acid sequences;
b) comparing the TCR profile with a database of previously identified diagnostic public T cell receptor sequences that are statistically significantly associated with CMV
infection;
c) generating a CMV burden score for the subject, wherein the CMV burden score is the proportion of unique TCR sequences in the profile of the subject that match the public TCR sequences in the database;
d) inputting the calculated CMV burden score from c) into a logistic regression model, wherein the logistic regression models compares CMV burden and CMV
infection status from a plurality of subjects of known CMV infection status; and e) determining an estimated probability of CMV infection status of the subject as the output of logistic regression model.
2. The method of claim 1, wherein the database of public T cell sequences are determined to be statistically significantly associated with CMV infection by:
a) obtaining unique TCR sequences from a group of subjects with CMV
infection and a group of subjects without CMV infection;
b) subjecting each unique TCR sequence to a one-tailed Fisher exact test based on presence or absence of each TCR sequence in subjects of known CMV infection status with the null hypothesis that each TCR sequence is no more common in subjects with a CMV
infection than in subjects without a CMV infection;
c) setting a nominal p value threshold;
d) controlling false discovery rate (FDR) by permutation of CMV infection status in each subject to generate an empirical null distribution of p-values; and e) generating a database of T cell receptor sequences that are statistically significantly shared in subjects with CMV infection.
a) obtaining unique TCR sequences from a group of subjects with CMV
infection and a group of subjects without CMV infection;
b) subjecting each unique TCR sequence to a one-tailed Fisher exact test based on presence or absence of each TCR sequence in subjects of known CMV infection status with the null hypothesis that each TCR sequence is no more common in subjects with a CMV
infection than in subjects without a CMV infection;
c) setting a nominal p value threshold;
d) controlling false discovery rate (FDR) by permutation of CMV infection status in each subject to generate an empirical null distribution of p-values; and e) generating a database of T cell receptor sequences that are statistically significantly shared in subjects with CMV infection.
3. The method of claim 2, wherein the nominal p-value threshold is less than or equal to 1.0*10 -4.
4. The method of claim 3, wherein the statistically identified unique TCR
sequences that correlate with the presence or absence of CMV infection comprise SEQ ID
NOs: 1-142.
sequences that correlate with the presence or absence of CMV infection comprise SEQ ID
NOs: 1-142.
5. The method of claim 1, wherein the step of determining TCR profile in step (b) comprises the steps of:
a) amplifying rearranged nucleic acids encoding the TCR CDR3 region in a multiplex PCR reaction with a mixture of forward primers specific to TCR V
gene segments and reverse primers specific to TCR J gene segments;
b) sequencing reads of the amplified nucleic acids;
c) processing the sequence reads to remove errors in the primary sequence of each read and to compress the data; and d) applying a nearest neighbor algorithm to collapse the data into unique sequences by merging closely related sequences to remove both PCR and sequencing errors.
a) amplifying rearranged nucleic acids encoding the TCR CDR3 region in a multiplex PCR reaction with a mixture of forward primers specific to TCR V
gene segments and reverse primers specific to TCR J gene segments;
b) sequencing reads of the amplified nucleic acids;
c) processing the sequence reads to remove errors in the primary sequence of each read and to compress the data; and d) applying a nearest neighbor algorithm to collapse the data into unique sequences by merging closely related sequences to remove both PCR and sequencing errors.
6. The method of claim 1, further comprising determining an HLA association for each unique TCR sequence in the TCR profile.
7. The method of claim 6, wherein said TCR sequence is associated with an HLA-A and/or an HLA-B allele.
8. A method for predicting the presence or absence of an infection in a subject of unknown infection status, the method comprising:
a) performing amplification and high throughput sequencing of genomic DNA
obtained from a sample comprising T cells obtained from the subject to determine a TCR
profile comprising unique CDR3 amino acid sequences;
b) comparing the TCR profile with a database of previously identified diagnostic public T cell receptor sequences that are statistically significantly associated with the infection;
c) generating a first score for the subject, wherein the first score is the proportion of unique TCR sequences in the profile of the subject that match the public TCR sequences in the database;
d) inputting the first score from c) into an algorithm, wherein the algorithm compares the first score of the subject and the infection status from a plurality of subjects of known infection status; and e) determining an estimated probability of infection status of the subject as the algorithm output.
a) performing amplification and high throughput sequencing of genomic DNA
obtained from a sample comprising T cells obtained from the subject to determine a TCR
profile comprising unique CDR3 amino acid sequences;
b) comparing the TCR profile with a database of previously identified diagnostic public T cell receptor sequences that are statistically significantly associated with the infection;
c) generating a first score for the subject, wherein the first score is the proportion of unique TCR sequences in the profile of the subject that match the public TCR sequences in the database;
d) inputting the first score from c) into an algorithm, wherein the algorithm compares the first score of the subject and the infection status from a plurality of subjects of known infection status; and e) determining an estimated probability of infection status of the subject as the algorithm output.
9. The method of claim 8, wherein the infection is from a cytomegalovirus (CMV), an Epstein-Barr virus (EBV), or a Herpes Simplex Virus (HSV).
10. A method for predicting the presence or absence of one or more viral infections in a subject of unknown infection status, the method comprising:
a) determining a profile of unique TCR sequences from a sample obtained from the subject;
b) inputting the unique TCR sequences from a) into one or more algorithms, wherein the one or more algorithms are generated by determining at least 10 5 unique TCR
sequences from each of a plurality of subjects of known infection status for each of the one or more infections and statistically identifying unique TCR sequences that correlate with the presence or absence of each of the one or more infections, to generate a score predictive of the presence or absence of each of the one or more infections; and c) inputting the scores from step b) into a logistic regression model trained on each of the plurality of subjects of known infection status for each of the one or more infections to predict whether the subject is either positive or negative for each of the one or more infections.
a) determining a profile of unique TCR sequences from a sample obtained from the subject;
b) inputting the unique TCR sequences from a) into one or more algorithms, wherein the one or more algorithms are generated by determining at least 10 5 unique TCR
sequences from each of a plurality of subjects of known infection status for each of the one or more infections and statistically identifying unique TCR sequences that correlate with the presence or absence of each of the one or more infections, to generate a score predictive of the presence or absence of each of the one or more infections; and c) inputting the scores from step b) into a logistic regression model trained on each of the plurality of subjects of known infection status for each of the one or more infections to predict whether the subject is either positive or negative for each of the one or more infections.
11. The method of claim 10, wherein the one or more viral infections are from a cytomegalovirus (CMV), an Epstein-Barr virus (EBV), a Herpes Simplex Virus (HSV) or a small pox virus.
12. A method for predicting a human leukocyte antigen (HLA) allele status of a subject, comprising:
(a) determining an immune receptor profile of unique T-cell receptor (TCR) rearranged DNA sequences for each of a plurality of subjects, each subject having a known HLA allele status;
(b) categorizing the plurality of subjects based on (i) said known HLA allele status of the subject and (ii) a presence or absence in the subject's immune receptor profile of a feature comprising a unique TCR rearranged DNA sequence;
(c) determining a statistical score for the association between a set of features and a positive HLA allele status based on (b);
(d) training a machine learning model using said set of features to define a set of classifiers for each HLA allele status;
(e) inputting one or more unique TCR rearranged DNA sequences of a subject with an unknown HLA allele status into said machine learning model to identify one or more features that match the set of classifiers; and (f) predicting an HLA allele status of said subject based on said one or more matched features.
(a) determining an immune receptor profile of unique T-cell receptor (TCR) rearranged DNA sequences for each of a plurality of subjects, each subject having a known HLA allele status;
(b) categorizing the plurality of subjects based on (i) said known HLA allele status of the subject and (ii) a presence or absence in the subject's immune receptor profile of a feature comprising a unique TCR rearranged DNA sequence;
(c) determining a statistical score for the association between a set of features and a positive HLA allele status based on (b);
(d) training a machine learning model using said set of features to define a set of classifiers for each HLA allele status;
(e) inputting one or more unique TCR rearranged DNA sequences of a subject with an unknown HLA allele status into said machine learning model to identify one or more features that match the set of classifiers; and (f) predicting an HLA allele status of said subject based on said one or more matched features.
13. The method of claim 12, wherein determining an immune receptor profile comprises determining the total number of unique TCR sequences and the frequency of each unique TCR sequence.
14. The method of claim 12, wherein determining a statistical score comprising determining a p-value using a Fisher exact two-tailed test.
15. The method of claim 14, further comprising determining a cutoff p-value for identifying a set of features that are significantly associated with an HLA
allele status.
allele status.
16. The method of claim 12, further comprising determining a false discovery rate (FDR) of the association of a feature with an HLA allele status.
17. The method of claim 16, further comprising determining a number of false-positive associations between said feature and said HLA allele status.
18. The method of claim 12, wherein training a machine learning model comprises training logistic regression model using said set of identified features and said known HLA
allele statuses of each subject.
allele statuses of each subject.
19. The method of claim 12, wherein training a machine learning model comprises performing a leave-one out cross validation method.
20. The method of claim 19, further comprising performing said leave-one out cross validation method for multiple rounds.
21. The method of claim 12, wherein said prediction is at least 80%
accurate.
accurate.
22. The method of claim 12, wherein said prediction is at least 90%
accurate.
accurate.
23. The method of claim 12, wherein said TCR rearranged DNA sequence is a TCRA, TCRB, TCRG or TCRD rearranged DNA sequence.
24. The method of claim 12, wherein said HLA allele is a HLA-A2 allele or a HLA-24 allele.
Applications Claiming Priority (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562120249P | 2015-02-24 | 2015-02-24 | |
US62/120,249 | 2015-02-24 | ||
US201562157249P | 2015-05-05 | 2015-05-05 | |
US62/157,249 | 2015-05-05 | ||
US201562215630P | 2015-09-08 | 2015-09-08 | |
US62/215,630 | 2015-09-08 | ||
PCT/US2016/019343 WO2016138122A1 (en) | 2015-02-24 | 2016-02-24 | Methods for diagnosing infectious disease and determining hla status using immune repertoire sequencing |
Publications (1)
Publication Number | Publication Date |
---|---|
CA2976580A1 true CA2976580A1 (en) | 2016-09-01 |
Family
ID=56789183
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA2976580A Abandoned CA2976580A1 (en) | 2015-02-24 | 2016-02-24 | Methods for diagnosing infectious disease and determining hla status using immune repertoire sequencing |
Country Status (6)
Country | Link |
---|---|
US (2) | US11047008B2 (en) |
EP (2) | EP3262196B1 (en) |
AU (2) | AU2016222788B2 (en) |
CA (1) | CA2976580A1 (en) |
ES (1) | ES2858306T3 (en) |
WO (1) | WO2016138122A1 (en) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2568509T3 (en) | 2009-01-15 | 2016-04-29 | Adaptive Biotechnologies Corporation | Adaptive immunity profiling and methods for the generation of monoclonal antibodies |
EP3388535B1 (en) | 2011-12-09 | 2021-03-24 | Adaptive Biotechnologies Corporation | Diagnosis of lymphoid malignancies and minimal residual disease detection |
JP6302847B2 (en) | 2012-03-05 | 2018-03-28 | アダプティヴ バイオテクノロジーズ コーポレーション | Determination of paired immunoreceptor chains from frequency matched subunits |
RU2631797C2 (en) | 2012-05-08 | 2017-09-26 | Эдэптив Байотекнолоджиз Корпорейшн | Compositions and methods of measurement and calibration of systematic mistake of amplification in multiplex pcr-reactions |
US9708657B2 (en) | 2013-07-01 | 2017-07-18 | Adaptive Biotechnologies Corp. | Method for generating clonotype profiles using sequence tags |
CA2941612A1 (en) | 2014-03-05 | 2015-09-11 | Adaptive Biotechnologies Corporation | Methods using randomer-containing synthetic molecules |
US11390921B2 (en) | 2014-04-01 | 2022-07-19 | Adaptive Biotechnologies Corporation | Determining WT-1 specific T cells and WT-1 specific T cell receptors (TCRs) |
US10066265B2 (en) | 2014-04-01 | 2018-09-04 | Adaptive Biotechnologies Corp. | Determining antigen-specific t-cells |
EP3132059B1 (en) | 2014-04-17 | 2020-01-08 | Adaptive Biotechnologies Corporation | Quantification of adaptive immune cell genomes in a complex mixture of cells |
EP3155426B1 (en) | 2014-06-13 | 2023-07-19 | Immudex ApS | General detection and isolation of specific cells by binding of labeled molecules |
WO2016069886A1 (en) | 2014-10-29 | 2016-05-06 | Adaptive Biotechnologies Corporation | Highly-multiplexed simultaneous detection of nucleic acids encoding paired adaptive immune receptor heterodimers from many samples |
US10246701B2 (en) | 2014-11-14 | 2019-04-02 | Adaptive Biotechnologies Corp. | Multiplexed digital quantitation of rearranged lymphoid receptors in a complex mixture |
EP3262196B1 (en) | 2015-02-24 | 2020-12-09 | Adaptive Biotechnologies Corp. | Method for determining hla status using immune repertoire sequencing |
SG11201707540QA (en) | 2015-03-16 | 2017-10-30 | Max-Delbrück-Centrum Für Molekulare Medizin In Der Helmholtz-Gemeinschaft | Method of detecting new immunogenic t cell epitopes and isolating new antigen-specific t cell receptors by means of an mhc cell library |
US11041202B2 (en) | 2015-04-01 | 2021-06-22 | Adaptive Biotechnologies Corporation | Method of identifying human compatible T cell receptors specific for an antigenic target |
US11514289B1 (en) * | 2016-03-09 | 2022-11-29 | Freenome Holdings, Inc. | Generating machine learning models using genetic data |
US10428325B1 (en) | 2016-09-21 | 2019-10-01 | Adaptive Biotechnologies Corporation | Identification of antigen-specific B cell receptors |
CN108070644B (en) * | 2016-11-08 | 2021-06-29 | 国家卫生计生委科学技术研究所 | Diagnosis system for gestational hypertension |
US11254980B1 (en) | 2017-11-29 | 2022-02-22 | Adaptive Biotechnologies Corporation | Methods of profiling targeted polynucleotides while mitigating sequencing depth requirements |
CA3095056A1 (en) * | 2018-04-13 | 2019-10-17 | Freenome Holdings, Inc. | Machine learning implementation for multi-analyte assay of biological samples |
CN111778355B (en) * | 2020-06-08 | 2023-05-05 | 璞晞(广州)生物免疫技术有限公司 | TCR sequence marker for hepatitis B virus infection and application thereof |
CN114203254B (en) * | 2021-12-02 | 2023-05-23 | 杭州艾沐蒽生物科技有限公司 | Method for analyzing immune characteristic related TCR based on artificial intelligence |
Family Cites Families (331)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5189147A (en) | 1984-06-13 | 1993-02-23 | Massachusetts Institute Of Technology | Meterodimeric T lymphocyte receptor antibody |
US4942124A (en) | 1987-08-11 | 1990-07-17 | President And Fellows Of Harvard College | Multiplex sequencing |
US5667967A (en) | 1990-05-01 | 1997-09-16 | The Board Of Trustees Of The Leland Stanford Junior University | T-cell receptor varible transcripts as disease related markers |
US5506126A (en) | 1988-02-25 | 1996-04-09 | The General Hospital Corporation | Rapid immunoselection cloning method |
WO1990004648A1 (en) | 1988-10-20 | 1990-05-03 | Alexander Alan Morley | Method for diagnosis of monoclonality in leukaemia and lymphoma |
US5143854A (en) | 1989-06-07 | 1992-09-01 | Affymax Technologies N.V. | Large scale photolithographic solid phase synthesis of polypeptides and receptor binding screening thereof |
US5336598A (en) | 1989-11-15 | 1994-08-09 | National Jewish Center For Immunology And Respiratory Medicine | Method for diagnosing a superantigen caused pathologial condition via assay of T-cells |
US5298396A (en) | 1989-11-15 | 1994-03-29 | National Jewish Center For Immunology And Respiratory Medicine | Method for identifying T cells disease involved in autoimmune disease |
US6916605B1 (en) | 1990-07-10 | 2005-07-12 | Medical Research Council | Methods for producing members of specific binding pairs |
EP0542810A1 (en) | 1990-08-02 | 1993-05-26 | B.R. Centre Limited | Methods for the production of proteins with a desired function |
IE76732B1 (en) | 1990-08-07 | 1997-11-05 | Becton Dickinson Co | One step test for absolute counts |
US5635354A (en) | 1991-01-09 | 1997-06-03 | Institut National De La Sante Et De La Recherche Medicale (Inserm) | Method for describing the repertoires of antibodies (Ab) and of T-cell receptors (TcR) of an individual's immune system |
JPH06502997A (en) | 1991-02-12 | 1994-04-07 | ルセル ユクラフ | Nucleotide sequence encoding the variable region of the β chain of human T lymphocyte receptor and corresponding peptide segment and use in diagnosis and therapy |
JP3080178B2 (en) | 1991-02-18 | 2000-08-21 | 東洋紡績株式会社 | Method for amplifying nucleic acid sequence and reagent kit therefor |
US6091000A (en) | 1991-03-15 | 2000-07-18 | Duke University | SCID mouse engrafted with human synovium tissue |
US20040018489A1 (en) | 1991-09-09 | 2004-01-29 | Ma Wupo | Detection of RNA |
US5674679A (en) | 1991-09-27 | 1997-10-07 | Amersham Life Science, Inc. | DNA cycle sequencing |
US5256542A (en) | 1992-03-09 | 1993-10-26 | Tanox Biosystems, Inc. | Selecting low frequency antigen-specific single B lymphocytes with correction for background noise |
US5213960A (en) | 1992-03-09 | 1993-05-25 | Tanox Biosystems, Inc. | Methods for selecting low frequency antigen-specific single B lymphocytes |
US5837447A (en) | 1992-04-15 | 1998-11-17 | Blood Center Research Foundation, Inc., The | Monitoring an immune response by analysis of amplified immunoglobulin or T-cell-receptor nucleic acid |
US5981176A (en) | 1992-06-17 | 1999-11-09 | City Of Hope | Method of detecting and discriminating between nucleic acid sequences |
US5925517A (en) | 1993-11-12 | 1999-07-20 | The Public Health Research Institute Of The City Of New York, Inc. | Detectably labeled dual conformation oligonucleotide probes, assays and kits |
WO1995028481A1 (en) | 1994-04-18 | 1995-10-26 | New York Society For The Ruptured And Crippled Maintaining The Hospital For Special Surgery | Conserved t-cell receptor sequences |
FR2724182B1 (en) | 1994-09-02 | 1996-12-13 | Pasteur Institut | OBTAINING A RECOMBINANT MONOCLONAL ANTIBODY FROM A HUMAN ANTI-RHESUS D MONOCLONAL ANTIBODY, ITS PRODUCTION IN INSECT CELLS, AND USES THEREOF |
US5846719A (en) | 1994-10-13 | 1998-12-08 | Lynx Therapeutics, Inc. | Oligonucleotide tags for sorting and identification |
US5604097A (en) | 1994-10-13 | 1997-02-18 | Spectragen, Inc. | Methods for sorting polynucleotides using oligonucleotide tags |
US5776737A (en) | 1994-12-22 | 1998-07-07 | Visible Genetics Inc. | Method and composition for internal identification of samples |
US6919434B1 (en) | 1995-02-20 | 2005-07-19 | Sankyo Co., Ltd. | Monoclonal antibodies that bind OCIF |
AU7437996A (en) | 1995-10-11 | 1997-04-30 | Leonard Adleman | Large scale dna sequencing by position sensitive hybridization |
WO1997013877A1 (en) | 1995-10-12 | 1997-04-17 | Lynx Therapeutics, Inc. | Measurement of gene expression profiles in toxicity determination |
US6087096A (en) | 1995-11-13 | 2000-07-11 | Dau; Peter C. | Method of intrafamily fragment analysis of the T cell receptor α and β chain CDR3 regions |
KR19990087749A (en) * | 1996-03-13 | 1999-12-27 | 시오노 요시히코 | Human T Cell Clones Specific to Rheumatoid Arthritis |
US6458530B1 (en) | 1996-04-04 | 2002-10-01 | Affymetrix Inc. | Selecting tag nucleic acids |
AU736742B2 (en) | 1996-06-03 | 2001-08-02 | Governors Of The University Of Alberta, The | Methods for detection of rearranged DNA |
ES2243393T3 (en) | 1996-06-04 | 2005-12-01 | University Of Utah Research Foundation | MONITORING OF HYBRIDIZATION DURING PCR. |
AU3878697A (en) | 1996-06-20 | 1998-02-02 | Cornell Research Foundation Inc. | Identification of abnormalities in the expression of t and cell antigen receptors as indicators of disease diagnosis, prognosis and therapeutic predictors |
CA2264925C (en) | 1996-09-06 | 2010-12-07 | Alain T. Luxembourg | Purification of antigen-specific t cells |
WO1998015644A2 (en) | 1996-09-27 | 1998-04-16 | The Chinese University Of Hong Kong | Parallel polynucleotide sequencing method |
US6228589B1 (en) | 1996-10-11 | 2001-05-08 | Lynx Therapeutics, Inc. | Measurement of gene expression profiles in toxicity determination |
GB9626815D0 (en) | 1996-12-23 | 1997-02-12 | Cemu Bioteknik Ab | Method of sequencing DNA |
JP2002503954A (en) | 1997-04-01 | 2002-02-05 | グラクソ、グループ、リミテッド | Nucleic acid amplification method |
DK1496120T3 (en) | 1997-07-07 | 2007-07-30 | Medical Res Council | In vitro sorting method |
US7572582B2 (en) | 1997-09-12 | 2009-08-11 | Exiqon A/S | Oligonucleotide analogues |
US6803019B1 (en) | 1997-10-15 | 2004-10-12 | Aclara Biosciences, Inc. | Laminate microstructure device and method for making same |
CA2307177C (en) | 1997-10-23 | 2004-06-29 | Exact Laboratories, Inc. | Methods for detecting contamination in molecular diagnostics using pcr |
US7351578B2 (en) | 1999-12-10 | 2008-04-01 | Invitrogen Corp. | Use of multiple recombination sites with unique specificity in recombinational cloning |
JP4262799B2 (en) | 1998-04-16 | 2009-05-13 | 平田機工株式会社 | Raw tire supply method |
DE19833738A1 (en) | 1998-07-27 | 2000-02-03 | Michael Giesing | Process for isolating cancer cells from cell-containing body fluids and kits for performing this process |
US6787308B2 (en) | 1998-07-30 | 2004-09-07 | Solexa Ltd. | Arrayed biomolecules and their use in sequencing |
DE19844931C1 (en) | 1998-09-30 | 2000-06-15 | Stefan Seeger | Procedures for DNA or RNA sequencing |
AR021833A1 (en) | 1998-09-30 | 2002-08-07 | Applied Research Systems | METHODS OF AMPLIFICATION AND SEQUENCING OF NUCLEIC ACID |
US6541608B1 (en) | 1999-02-23 | 2003-04-01 | Baylor College Of Medicine | T cell receptor Vβ-Dβ-Jβ sequence and methods for its detection |
US6300070B1 (en) | 1999-06-04 | 2001-10-09 | Mosaic Technologies, Inc. | Solid phase methods for amplifying multiple nucleic acids |
US6440706B1 (en) | 1999-08-02 | 2002-08-27 | Johns Hopkins University | Digital amplification |
CA2589418A1 (en) | 1999-08-24 | 2001-03-01 | Medarex, Inc. | Human ctla-4 antibodies and their uses |
US20040209314A1 (en) | 1999-09-06 | 2004-10-21 | Institut National De La Sante Et De La Recherche Medicale France | Means for detection and purification of CD8+ T lymphocyte populations specific to peptides presented in the context of HLA |
US6235483B1 (en) | 2000-01-31 | 2001-05-22 | Agilent Technologies, Inc. | Methods and kits for indirect labeling of nucleic acids |
US20040121364A1 (en) | 2000-02-07 | 2004-06-24 | Mark Chee | Multiplex nucleic acid reactions |
US7955794B2 (en) | 2000-09-21 | 2011-06-07 | Illumina, Inc. | Multiplex nucleic acid reactions |
CA2404568A1 (en) | 2000-03-31 | 2001-10-11 | University Of Southern California | Epigenetic sequences for esophageal adenocarcinoma |
AUPQ687600A0 (en) | 2000-04-13 | 2000-05-11 | Flinders Technologies Pty Ltd | A method of detection |
US20030207300A1 (en) | 2000-04-28 | 2003-11-06 | Matray Tracy J. | Multiplex analytical platform using molecular tags |
US6596492B2 (en) | 2000-07-11 | 2003-07-22 | Colorado State University Research Foundation | PCR materials and methods useful to detect canine and feline lymphoid malignancies |
US7567870B1 (en) | 2000-07-31 | 2009-07-28 | Institute For Systems Biology | Multiparameter analysis for predictive medicine |
AU2001284703B2 (en) | 2000-08-03 | 2007-03-22 | Therapeutic Human Polyclonals Inc. | Production of humanized antibodies in transgenic animals |
US6939451B2 (en) | 2000-09-19 | 2005-09-06 | Aclara Biosciences, Inc. | Microfluidic chip having integrated electrodes |
WO2002057491A2 (en) | 2000-10-24 | 2002-07-25 | The Board Of Trustees Of The Leland Stanford Junior University | Direct multiplex characterization of genomic dna |
US6778724B2 (en) | 2000-11-28 | 2004-08-17 | The Regents Of The University Of California | Optical switching and sorting of biological samples and microparticles transported in a micro-fluidic device, including integrated bio-chip devices |
US6947748B2 (en) | 2000-12-15 | 2005-09-20 | Adaptix, Inc. | OFDMA with adaptive subcarrier-cluster configuration and selective loading |
US7148040B2 (en) | 2001-02-20 | 2006-12-12 | University Of Georgia Research Foundation, Inc. | Method of rapid production of hybridomas expressing monoclonal antibodies on the cell surface |
US7265208B2 (en) | 2001-05-01 | 2007-09-04 | The Regents Of The University Of California | Fusion molecules and treatment of IgE-mediated allergic diseases |
WO2002094009A2 (en) | 2001-05-24 | 2002-11-28 | Guard Inc. | Methods for selecting and producing animals having a predicted level of immune response |
US20050260570A1 (en) | 2001-05-29 | 2005-11-24 | Mao Jen-I | Sequencing by proxy |
US7026121B1 (en) | 2001-06-08 | 2006-04-11 | Expression Diagnostics, Inc. | Methods and compositions for diagnosing and monitoring transplant rejection |
US6720144B1 (en) | 2001-06-21 | 2004-04-13 | Quest Diagnostics | Detection of clonal T-cell receptor-γ gene rearrangement by PCR/temporal temperature gradient gel electrophoresis (TTGE) |
WO2003008624A2 (en) | 2001-07-15 | 2003-01-30 | Keck Graduate Institute | Nucleic acid amplification using nicking agents |
WO2003020983A1 (en) | 2001-08-30 | 2003-03-13 | Virginia Commonwealth University | Allele specific pcr for genotyping |
ES2239246T3 (en) | 2001-08-31 | 2005-09-16 | Avidex Limited | SOLUBLE RECEIVER OF CELLS T. |
AU2002322518A1 (en) | 2001-08-31 | 2003-03-18 | Rosetta Inpharmactis Llc | Methods for preparing nucleic acid samples |
DK1436404T3 (en) | 2001-09-19 | 2010-03-08 | Alexion Pharma Inc | Manipulated templates and their use in single-primer amplification |
EP1451340B1 (en) | 2001-11-09 | 2014-01-08 | Life Technologies Corporation | Identification, monitoring and treatment of disease and characterization of biological condition using gene expression profiles |
GB2382137A (en) | 2001-11-20 | 2003-05-21 | Mats Gullberg | Nucleic acid enrichment |
GB0128153D0 (en) | 2001-11-23 | 2002-01-16 | Bayer Ag | Profiling of the immune gene repertoire |
WO2003052101A1 (en) | 2001-12-14 | 2003-06-26 | Rosetta Inpharmatics, Inc. | Sample tracking using molecular barcodes |
GB0130267D0 (en) | 2001-12-19 | 2002-02-06 | Neutec Pharma Plc | Focussed antibody technology |
WO2003059155A2 (en) | 2002-01-09 | 2003-07-24 | Maxx Genetech Co. Ltd | Method of detecting t-cell proliferation for diagnosis of diseases by gene array |
US7323306B2 (en) | 2002-04-01 | 2008-01-29 | Brookhaven Science Associates, Llc | Genome signature tags |
US7157274B2 (en) | 2002-06-24 | 2007-01-02 | Cytonome, Inc. | Method and apparatus for sorting particles |
AU2003299503A1 (en) | 2002-05-16 | 2004-06-15 | Vanderbilt University | Method for predicting autoimmune diseases |
EP1573636A2 (en) | 2002-07-01 | 2005-09-14 | Institut Pasteur | System, method, device, and computer program product for extraction, gathering, manipulation, and analysis of peak data from an automated sequencer |
EP1517995A4 (en) | 2002-07-03 | 2007-07-11 | Inst Scient Res Inc | Compositions and methods for the detection of human t cell receptor variable family gene expression |
CA2800113C (en) | 2002-07-12 | 2015-02-03 | The Johns Hopkins University | Reagents and methods for engaging unique clonotypic lymphocyte receptors |
AU2003261168A1 (en) | 2002-07-19 | 2004-02-09 | Althea Technologies, Inc. | Strategies for gene expression analysis |
US7157228B2 (en) | 2002-09-09 | 2007-01-02 | Bioarray Solutions Ltd. | Genetic analysis and authentication |
US7459273B2 (en) | 2002-10-04 | 2008-12-02 | Affymetrix, Inc. | Methods for genotyping selected polymorphism |
AU2003282612A1 (en) | 2002-10-11 | 2004-05-04 | The Regents Of The University Of California | A method for diagnosis and prognosis of multiple sclerosis |
ATE483822T1 (en) | 2002-10-11 | 2010-10-15 | Univ Erasmus | PRIMERS FOR NUCLEIC ACID AMPLIFICATION IN PCR-BASED CLONALITY STUDIES |
WO2004044209A1 (en) | 2002-11-13 | 2004-05-27 | Monoquant Pty Ltd | A method of detection |
US20060134704A1 (en) | 2002-11-14 | 2006-06-22 | Atsushi Muraguchi | Microwell array chip for detecting antigen-specific lymphocytes, method of detecting and method of manufacturing antigen-specific lymphocytes, and method of cloning antigen-specific lymphocyte antigen receptor genes |
AU2003300919A1 (en) | 2002-12-11 | 2004-06-30 | Coley Pharmaceutical Gmbh | 5' cpg nucleic acids and methods of use |
WO2004063706A2 (en) | 2003-01-08 | 2004-07-29 | Maxx Genetech Co., Ltd. | Method of detecting over-expression of t-cell receptor genes by real-time pcr |
US7244567B2 (en) | 2003-01-29 | 2007-07-17 | 454 Life Sciences Corporation | Double ended sequencing |
GB0304068D0 (en) | 2003-02-22 | 2003-03-26 | Avidex Ltd | Substances |
US20070042349A1 (en) | 2003-04-24 | 2007-02-22 | Ogle Brenda M | Methods for assessing biologic diversity |
AU2003902299A0 (en) | 2003-05-13 | 2003-05-29 | Flinders Medical Centre | A method of analysing a marker nucleic acid molecule |
WO2005020784A2 (en) | 2003-05-23 | 2005-03-10 | Mount Sinai School Of Medicine Of New York University | Surrogate cell gene expression signatures for evaluating the physical state of a subject |
US20050010030A1 (en) | 2003-07-02 | 2005-01-13 | Zang Jingwu Z. | T cell receptor CDR3 sequence and methods for detecting and treating rheumatoid arthritis |
EP1641809B2 (en) | 2003-07-05 | 2018-10-03 | The Johns Hopkins University | Method and compositions for detection and enumeration of genetic variations |
US7574305B2 (en) | 2003-07-15 | 2009-08-11 | Bioarray Solutions Ltd. | Concurrent optimization in selection of primer and capture probe sets for nucleic acid analysis |
US20060228350A1 (en) | 2003-08-18 | 2006-10-12 | Medimmune, Inc. | Framework-shuffling of antibodies |
US20050048498A1 (en) | 2003-08-29 | 2005-03-03 | Applera Corporation | Compositions, methods, and kits for assembling probes |
WO2005026686A2 (en) | 2003-09-09 | 2005-03-24 | Compass Genetics, Llc | Multiplexed analytical platform |
TWI333977B (en) | 2003-09-18 | 2010-12-01 | Symphogen As | Method for linking sequences of interest |
FR2863274B1 (en) | 2003-12-05 | 2012-11-16 | Commissariat Energie Atomique | METHOD FOR QUANTITATIVE EVALUATION OF REARRANGEMENT OR TARGETED GENETIC RECOMBINATION OF AN INDIVIDUAL AND APPLICATIONS THEREOF |
US20070238099A1 (en) | 2003-12-08 | 2007-10-11 | Cohen Irun R | Antigen Receptor Variable Region Typing |
DE60326052D1 (en) | 2003-12-15 | 2009-03-19 | Pasteur Institut | Determination of the repertoire of B lymphocyte populations |
US20080166718A1 (en) | 2003-12-15 | 2008-07-10 | Institut Pasteur | Repertoire determination of a lymphocyte B population |
WO2005059176A1 (en) | 2003-12-15 | 2005-06-30 | Institut Pasteur | Repertoire determination of a lymphocyte b population |
CA2556981C (en) | 2004-02-18 | 2015-10-13 | The Trustees Of Boston University | Method for detecting and quantifying rare mutations or polymorphisms |
US20060046258A1 (en) | 2004-02-27 | 2006-03-02 | Lapidus Stanley N | Applications of single molecule sequencing |
US20070161001A1 (en) | 2004-03-04 | 2007-07-12 | Dena Leshkowitz | Quantifying and profiling antibody and t cell receptor gene expression |
JP4480423B2 (en) | 2004-03-08 | 2010-06-16 | 独立行政法人科学技術振興機構 | Method for determining the presence or absence of expansion of immune cell clones |
DE102004016437A1 (en) | 2004-04-04 | 2005-10-20 | Oligene Gmbh | Method for detecting signatures in complex gene expression profiles |
US20050250147A1 (en) | 2004-05-10 | 2005-11-10 | Macevicz Stephen C | Digital profiling of polynucleotide populations |
JP2007536939A (en) | 2004-05-14 | 2007-12-20 | アモークス インコーポレーティッド | Immune cell biosensor and method of use thereof |
EP1598429A1 (en) | 2004-05-19 | 2005-11-23 | Amplion Ltd. | Detection of amplicon contamination during PCR exhibiting two different annealing temperatures |
GB0412973D0 (en) | 2004-06-10 | 2004-07-14 | Celltech R&D Ltd | Identification of antibody producing cells |
US20060020397A1 (en) | 2004-07-21 | 2006-01-26 | Kermani Bahram G | Methods for nucleic acid and polypeptide similarity search employing content addressable memories |
WO2006014869A1 (en) | 2004-07-26 | 2006-02-09 | Parallele Bioscience, Inc. | Simultaneous analysis of multiple genomes |
US20060094018A1 (en) | 2004-08-03 | 2006-05-04 | Bauer A R Jr | Discovery and a method for the early detection of pancreatic cancer and other disease conditions |
US7820382B2 (en) | 2004-08-03 | 2010-10-26 | Bauer A Robert | Method for the early detection of breast cancer, lung cancer, pancreatic cancer and colon polyps, growths and cancers as well as other gastrointestinal disease conditions and the preoperative and postoperative monitoring of transplanted organs from the donor and in the recipient and their associated conditions related and unrelated to the organ transplantation |
US7476724B2 (en) | 2004-08-05 | 2009-01-13 | Genentech, Inc. | Humanized anti-cmet antibodies |
CA2580412A1 (en) | 2004-09-13 | 2006-03-23 | Government Of The United States Of America, Represented By The Secretary , Department Of Health And Human Services | Compositions comprising t cell receptors and methods of use thereof |
US7170050B2 (en) | 2004-09-17 | 2007-01-30 | Pacific Biosciences Of California, Inc. | Apparatus and methods for optical analysis of molecules |
WO2006050138A2 (en) | 2004-10-29 | 2006-05-11 | Benaroya Research Institute At Virginia Mason | Methods of generating antigen-specific cd4+cd25+ regulatory t cells, compositions and methods of use |
US7966988B2 (en) | 2005-01-11 | 2011-06-28 | Exxonmobil Research And Engineering Company | Method for controlling soot induced lubricant viscosity increase |
US8029783B2 (en) | 2005-02-02 | 2011-10-04 | Genentech, Inc. | DR5 antibodies and articles of manufacture containing same |
FR2881436B1 (en) | 2005-02-03 | 2007-04-27 | Commissariat Energie Atomique | METHOD FOR DETERMINING THE DIVERSITY OF T LYMPHOCYTES IN A BIOLOGICAL SAMPLE |
US7393665B2 (en) | 2005-02-10 | 2008-07-01 | Population Genetics Technologies Ltd | Methods and compositions for tagging and identifying polynucleotides |
US20080280774A1 (en) | 2005-02-16 | 2008-11-13 | Wyeth | Methods and Systems for Diagnosis, Prognosis and Selection of Treatment of Leukemia |
US7537894B2 (en) | 2005-03-02 | 2009-05-26 | The University Of Chicago | Methods and kits for monitoring Barrett's metaplasia |
US20060216737A1 (en) | 2005-03-10 | 2006-09-28 | John Bodeau | Methods for multiplex amplification |
WO2006099604A2 (en) | 2005-03-16 | 2006-09-21 | Compass Genetics, Llc | Methods and compositions for assay readouts on multiple analytical platforms |
EP1877576B1 (en) | 2005-04-12 | 2013-01-23 | 454 Life Sciences Corporation | Methods for determining sequence variants using ultra-deep sequencing |
WO2006116155A2 (en) | 2005-04-21 | 2006-11-02 | The Regents Of The University Of California | A method for diagnosis and prognosis of multiple sclerosis subtypes |
US20060263789A1 (en) | 2005-05-19 | 2006-11-23 | Robert Kincaid | Unique identifiers for indicating properties associated with entities to which they are attached, and methods for using |
US7208795B2 (en) | 2005-05-24 | 2007-04-24 | Atmel Corporation | Low-cost, low-voltage single-layer polycrystalline EEPROM memory cell integration into BiCMOS technology |
US8445194B2 (en) | 2005-06-15 | 2013-05-21 | Callida Genomics, Inc. | Single molecule arrays for genetic and chemical analysis |
US20070020670A1 (en) | 2005-07-07 | 2007-01-25 | Hematologics, Inc. | Methods for detecting and confirming minimal disease |
US20070020640A1 (en) | 2005-07-21 | 2007-01-25 | Mccloskey Megan L | Molecular encoding of nucleic acid templates for PCR and other forms of sequence analysis |
GB0521521D0 (en) | 2005-10-21 | 2005-11-30 | Medical Res Council | Diagnostic methods and kits |
GB0522310D0 (en) | 2005-11-01 | 2005-12-07 | Solexa Ltd | Methods of preparing libraries of template polynucleotides |
US20070105165A1 (en) | 2005-11-04 | 2007-05-10 | Charles Goolsby | Composite profiles of cell antigens and target signal transduction proteins for analysis and clinical management of hematologic cancers |
US7375211B2 (en) | 2005-11-18 | 2008-05-20 | Kou Zhong C | Method for detection and quantification of T-cell receptor Vβ repertoire |
US8137936B2 (en) | 2005-11-29 | 2012-03-20 | Macevicz Stephen C | Selected amplification of polynucleotides |
EP1987162A4 (en) | 2006-01-23 | 2009-11-25 | Population Genetics Technologi | Nucleic acid analysis using sequence tokens |
WO2007087312A2 (en) | 2006-01-23 | 2007-08-02 | Population Genetics Technologies Ltd. | Molecular counting |
SG10201405158QA (en) | 2006-02-24 | 2014-10-30 | Callida Genomics Inc | High throughput genome sequencing on dna arrays |
US20100040606A1 (en) | 2006-03-06 | 2010-02-18 | Symphogen A/S | Recombinant polyclonal antibody for treatment of respiratory syncytial virus infections |
ES2546848T3 (en) | 2006-03-10 | 2015-09-29 | Epigenomics Ag | A method to identify a biological sample for methylation analysis |
EP2002017B1 (en) | 2006-04-04 | 2015-06-10 | Keygene N.V. | High throughput detection of molecular markers based on restriction fragments |
WO2007134220A2 (en) | 2006-05-11 | 2007-11-22 | University Of Maryland Biotechnology Institute | A general method for generating human antibody responses in vitro |
WO2007136518A2 (en) | 2006-05-17 | 2007-11-29 | Torrey Pines Institute For Molecular Studies | Treatment of autoimmune disorders |
AU2006345511B2 (en) | 2006-05-25 | 2013-03-21 | Institute For Advanced Study | Methods for identifying sequence motifs, and applications thereof |
US7833716B2 (en) | 2006-06-06 | 2010-11-16 | Gen-Probe Incorporated | Tagged oligonucleotides and their use in nucleic acid amplification methods |
CN101467082B (en) | 2006-06-12 | 2011-12-14 | 加利福尼亚太平洋生物科学公司 | Substrates for performing analytical reactions |
EP2589668A1 (en) | 2006-06-14 | 2013-05-08 | Verinata Health, Inc | Rare cell analysis using sample splitting and DNA tags |
US20100027896A1 (en) | 2006-06-28 | 2010-02-04 | Amir Geva | Automated application interaction using a virtual operator |
JP2009542220A (en) | 2006-06-30 | 2009-12-03 | アプライド バイオシステムズ, エルエルシー | Reversible terminator nucleotides and methods of use |
WO2008108803A2 (en) | 2006-07-13 | 2008-09-12 | Amaox, Ltd. | Immune cell biosensors and methods of using same |
US8394590B2 (en) | 2006-08-02 | 2013-03-12 | California Institute Of Technology | Capture agents and related methods and systems for detecting and/or sorting targets |
US20080274904A1 (en) | 2006-08-10 | 2008-11-06 | Niall Anthony Gormley | Method of target enrichment |
WO2008026927A2 (en) | 2006-08-30 | 2008-03-06 | Academisch Medisch Centrum | Process for displaying t- and b-cell receptor repertoires |
US8088379B2 (en) | 2006-09-26 | 2012-01-03 | The United States Of America As Represented By The Department Of Health And Human Services | Modified T cell receptors and related materials and methods |
WO2008039694A2 (en) | 2006-09-26 | 2008-04-03 | St. Jude Children's Research Hospital | Methods and compositions for monitoring t cell receptor diversity |
JP5026047B2 (en) | 2006-10-18 | 2012-09-12 | 株式会社膠原病研究所 | Method for identifying self-responsive T cells or T cell receptors involved in the development of autoimmune diseases and use thereof |
EP2102367A2 (en) | 2006-11-09 | 2009-09-23 | XDX, Inc. | Methods for diagnosing and monitoring the status of systemic lupus erythematosus |
US8262900B2 (en) | 2006-12-14 | 2012-09-11 | Life Technologies Corporation | Methods and apparatus for measuring analytes using large scale FET arrays |
US7862999B2 (en) | 2007-01-17 | 2011-01-04 | Affymetrix, Inc. | Multiplex targeted amplification using flap nuclease |
EP2121983A2 (en) | 2007-02-02 | 2009-11-25 | Illumina Cambridge Limited | Methods for indexing samples and sequencing multiple nucleotide templates |
AU2008221119B2 (en) | 2007-03-01 | 2013-11-14 | Les Laboratoires Servier | Method for cloning cognate antibodies |
WO2008147879A1 (en) | 2007-05-22 | 2008-12-04 | Ryan Golhar | Automated method and device for dna isolation, sequence determination, and identification |
WO2008150432A1 (en) | 2007-06-01 | 2008-12-11 | 454 Life Sciences Corporation | System and meth0d for identification of individual samples from a multiplex mixture |
WO2009015296A1 (en) | 2007-07-24 | 2009-01-29 | The Regents Of The University Of California | Microfabricated dropley generator |
JP2010534474A (en) | 2007-07-26 | 2010-11-11 | パシフィック バイオサイエンシーズ オブ カリフォルニア, インコーポレイテッド | Molecular redundant sequencing |
ITRM20070429A1 (en) | 2007-08-06 | 2009-02-07 | Uni Cattolica Del Sacro Cuor E | MEANS FOR DIAGNOSIS PREVENTION AND CARE OF RHEUMATOID ARTHRITIS. |
US20090042737A1 (en) | 2007-08-09 | 2009-02-12 | Katz Andrew S | Methods and Devices for Correlated, Multi-Parameter Single Cell Measurements and Recovery of Remnant Biological Material |
US8268564B2 (en) | 2007-09-26 | 2012-09-18 | President And Fellows Of Harvard College | Methods and applications for stitched DNA barcodes |
WO2009045898A2 (en) | 2007-09-28 | 2009-04-09 | Mayo Foundation For Medical Education And Research | Assessing t cell repertoires |
US7960116B2 (en) | 2007-09-28 | 2011-06-14 | Pacific Biosciences Of California, Inc. | Nucleic acid sequencing methods and systems |
KR101620642B1 (en) | 2007-11-07 | 2016-05-12 | 제넨테크, 인크. | Methods and compositions for assessing responsiveness of b-cell lymphoma to treatment with anti-cd40 antibodies |
EP2060637A1 (en) | 2007-11-14 | 2009-05-20 | Koninklijke Philips Electronics N.V. | Means and methods for detection of nucleic acids |
WO2009070767A2 (en) | 2007-11-28 | 2009-06-04 | Whitehead Institute For Biomedical Research | Systemic instigation systems to study tumor growth or metastasis |
EP2062982A1 (en) | 2007-11-26 | 2009-05-27 | ImmunID | Method for studying the V(D)J combinatorial diversity. |
CN101225441B (en) | 2007-12-05 | 2010-12-01 | 浙江大学 | Method for detecting genetic constitution of clone-specific T lymphocyte TCR BV CDR3 |
US20100021894A1 (en) | 2007-12-20 | 2010-01-28 | Northwestern University | Nanoparticle-Based Colorimetric Detection Of Cysteine |
US8621502B2 (en) | 2007-12-21 | 2013-12-31 | Microsoft Corporation | Obtaining user reactions to video |
US7767400B2 (en) | 2008-02-03 | 2010-08-03 | Helicos Biosciences Corporation | Paired-end reads in sequencing by synthesis |
EP2088432A1 (en) | 2008-02-11 | 2009-08-12 | MorphoSys AG | Methods for identification of an antibody or a target |
ES2600165T3 (en) | 2008-02-28 | 2017-02-07 | The Ohio State University Research Foundation | MiR-32 antagonists to increase the response of prostate cancer to apoptosis |
JP2011516033A (en) | 2008-02-28 | 2011-05-26 | ジ・オハイオ・ステイト・ユニバーシティ・リサーチ・ファウンデイション | MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid leukemia (AML) and uses thereof |
US20090226975A1 (en) | 2008-03-10 | 2009-09-10 | Illumina, Inc. | Constant cluster seeding |
TW200938840A (en) | 2008-03-12 | 2009-09-16 | Emo Biomedicine Corp | A method for in vitro study of immune modulation using pig blood cell |
US8143007B2 (en) | 2008-03-13 | 2012-03-27 | National Institute Of Immunology | Nested primer sets for amplifying mouse immunoglobulin variable gene segments |
CA2991818C (en) | 2008-03-28 | 2022-10-11 | Pacific Biosciences Of California, Inc. | Compositions and methods for nucleic acid sequencing |
EP2107040B1 (en) | 2008-03-31 | 2011-10-26 | Sony Deutschland Gmbh | A method of fabricating a membrane having a tapered pore |
AU2009244634B2 (en) | 2008-04-16 | 2014-12-18 | iRepertoire, Inc. | Method for evaluating and comparing immunorepertoires |
US8911948B2 (en) | 2008-04-30 | 2014-12-16 | Integrated Dna Technologies, Inc. | RNase H-based assays utilizing modified RNA monomers |
CN102171569A (en) | 2008-05-09 | 2011-08-31 | 杜克大学 | Autoantibodies in the detection and treatment of cancer |
US9068181B2 (en) | 2008-05-23 | 2015-06-30 | The General Hospital Corporation | Microfluidic droplet encapsulation |
DE102008025656B4 (en) | 2008-05-28 | 2016-07-28 | Genxpro Gmbh | Method for the quantitative analysis of nucleic acids, markers therefor and their use |
CN102112624A (en) | 2008-05-30 | 2011-06-29 | 麻省理工学院 | Compositions and methods for spatial separation and screening of cells |
WO2009151628A2 (en) | 2008-06-12 | 2009-12-17 | Gov't Of The Usa, As Represented By The Secretary, Department Of Health Human Services | Monitoring tcr-b to determine hiv therapy and disease progression |
EP2300823A4 (en) | 2008-06-25 | 2012-03-14 | Baylor Res Inst | Blood transcriptional signature of mycobacterium tuberculosis infection |
CA2769002C (en) | 2008-07-24 | 2018-05-08 | The Board Of Regents Of The University Of Texas System | Vh4 codon signature for multiple sclerosis |
WO2010011944A2 (en) | 2008-07-25 | 2010-01-28 | Wagner Richard W | Protein screeing methods |
TWI408888B (en) | 2008-08-07 | 2013-09-11 | Ind Tech Res Inst | Ultrasonic linear motor |
US8586310B2 (en) | 2008-09-05 | 2013-11-19 | Washington University | Method for multiplexed nucleic acid patch polymerase chain reaction |
US9156010B2 (en) | 2008-09-23 | 2015-10-13 | Bio-Rad Laboratories, Inc. | Droplet-based assay system |
CA3075139C (en) | 2008-09-23 | 2022-04-12 | Bio-Rad Laboratories, Inc. | Droplet-based assay system |
US8699361B2 (en) | 2008-09-30 | 2014-04-15 | Qualcomm Incorporated | Out-of-synchronization handling method and apparatus |
US20100137143A1 (en) | 2008-10-22 | 2010-06-03 | Ion Torrent Systems Incorporated | Methods and apparatus for measuring analytes |
US8546128B2 (en) | 2008-10-22 | 2013-10-01 | Life Technologies Corporation | Fluidics system for sequential delivery of reagents |
US8691510B2 (en) | 2008-11-07 | 2014-04-08 | Sequenta, Inc. | Sequence analysis of complex amplicons |
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 |
US8748103B2 (en) | 2008-11-07 | 2014-06-10 | Sequenta, Inc. | 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 |
AU2009311588B2 (en) | 2008-11-07 | 2016-07-28 | Adaptive Biotechnologies Corp. | Methods of monitoring conditions by sequence analysis |
US20140234835A1 (en) | 2008-11-07 | 2014-08-21 | Sequenta, Inc. | Rare clonotypes and uses thereof |
US9506119B2 (en) | 2008-11-07 | 2016-11-29 | Adaptive Biotechnologies Corp. | Method of sequence determination using sequence tags |
US8628927B2 (en) | 2008-11-07 | 2014-01-14 | Sequenta, Inc. | Monitoring health and disease status using clonotype profiles |
WO2010060051A2 (en) | 2008-11-21 | 2010-05-27 | Emory University | Systems biology approach predicts the immunogenicity of vaccines |
US8367330B2 (en) | 2008-12-22 | 2013-02-05 | Quest Diagnostics Investments Incorporated | Methods for detecting TCR-gamma gene rearrangement |
AU2009333489A1 (en) | 2008-12-30 | 2010-07-08 | Centocor Ortho Biotech Inc. | Serum markers predicting clinical response to anti-TNF antibodies in patients with ankylosing spondylitis |
ES2568509T3 (en) | 2009-01-15 | 2016-04-29 | Adaptive Biotechnologies Corporation | Adaptive immunity profiling and methods for the generation of monoclonal antibodies |
US9329170B2 (en) | 2009-01-20 | 2016-05-03 | The Board Of Trustees Of The Leland Stanford Junior University | Single cell gene expression for diagnosis, prognosis and identification of drug targets |
US20100323348A1 (en) | 2009-01-31 | 2010-12-23 | The Regents Of The University Of Colorado, A Body Corporate | Methods and Compositions for Using Error-Detecting and/or Error-Correcting Barcodes in Nucleic Acid Amplification Process |
US8673627B2 (en) | 2009-05-29 | 2014-03-18 | Life Technologies Corporation | Apparatus and methods for performing electrochemical reactions |
US8574835B2 (en) | 2009-05-29 | 2013-11-05 | Life Technologies Corporation | Scaffolded nucleic acid polymer particles and methods of making and using |
US20120058902A1 (en) | 2009-06-25 | 2012-03-08 | Livingston Robert J | Method of measuring adaptive immunity |
WO2010151416A1 (en) | 2009-06-25 | 2010-12-29 | Fred Hutchinson Cancer Research Center | Method of measuring adaptive immunity |
JP5829606B2 (en) | 2009-06-29 | 2015-12-09 | カリフォルニア・インスティテュート・オブ・テクノロジーCalifornia Institute Oftechnology | Isolation of unknown rearranged T cell receptors from single cells |
SG177688A1 (en) | 2009-07-17 | 2012-03-29 | Univ Brown | Detection of short rna sequences |
US8586359B2 (en) | 2009-07-28 | 2013-11-19 | Promising Furture, LLC | Compositions and methods of preparing alloreactive cytotoxic T cells |
WO2011035870A1 (en) | 2009-09-22 | 2011-03-31 | Roche Diagnostics Gmbh | Determination of kir haplotypes associated with disease |
GB0917094D0 (en) | 2009-09-29 | 2009-11-11 | Ucl Biomedica Plc | T-cell receptor |
US9043160B1 (en) | 2009-11-09 | 2015-05-26 | Sequenta, Inc. | Method of determining clonotypes and clonotype profiles |
US8835358B2 (en) | 2009-12-15 | 2014-09-16 | Cellular Research, Inc. | Digital counting of individual molecules by stochastic attachment of diverse labels |
US9315857B2 (en) | 2009-12-15 | 2016-04-19 | Cellular Research, Inc. | Digital counting of individual molecules by stochastic attachment of diverse label-tags |
US8545248B2 (en) | 2010-01-07 | 2013-10-01 | Life Technologies Corporation | System to control fluid flow based on a leak detected by a sensor |
KR101323827B1 (en) | 2010-01-08 | 2013-10-31 | 키스트 유럽 에프게엠베하 | Primers for diagnosing ankylosing spondylitis, and method for diagnosing ankylosing spondylitis using the same |
GB201000375D0 (en) | 2010-01-09 | 2010-02-24 | Univ Cardiff | T Cell clonotypes |
WO2011106738A2 (en) * | 2010-02-25 | 2011-09-01 | Fred Hutchinson Cancer Research Center | Use of tcr clonotypes as biomarkers for disease |
EP2367000A1 (en) * | 2010-03-04 | 2011-09-21 | Charité Universitätsmedizin Berlin | High throughput analysis of T-cell receptor repertoires |
JP6158080B2 (en) | 2010-05-06 | 2017-07-05 | アダプティヴ バイオテクノロジーズ コーポレーション | Health and disease status monitoring using chronotype profiles |
US9234240B2 (en) | 2010-05-07 | 2016-01-12 | The Board Of Trustees Of The Leland Stanford Junior University | Measurement and comparison of immune diversity by high-throughput sequencing |
US20130123120A1 (en) | 2010-05-18 | 2013-05-16 | Natera, Inc. | Highly Multiplex PCR Methods and Compositions |
EP2580353B1 (en) | 2010-06-11 | 2015-07-29 | Life Technologies Corporation | Alternative nucleotide flows in sequencing-by-synthesis methods |
SG186787A1 (en) | 2010-07-23 | 2013-02-28 | Esoterix Genetic Lab Llc | Identification of differentially represented fetal or maternal genomic regions and uses thereof |
WO2012017081A1 (en) | 2010-08-06 | 2012-02-09 | Ludwig-Maximilians-Universität München | Identification of t cell target antigens |
ME02810B (en) | 2010-09-20 | 2018-01-20 | Biontech Cell & Gene Therapies Gmbh | Antigen-specific t cell receptors and t cell epitopes |
WO2012038839A2 (en) | 2010-09-21 | 2012-03-29 | Population Genetics Technologies Ltd. | Increasing confidence of allele calls with molecular counting |
EP3447155A1 (en) | 2010-09-30 | 2019-02-27 | Raindance Technologies, Inc. | Sandwich assays in droplets |
EP2625295B1 (en) | 2010-10-08 | 2019-03-13 | President and Fellows of Harvard College | High-throughput immune sequencing |
EP2625320B1 (en) | 2010-10-08 | 2019-03-27 | President and Fellows of Harvard College | High-throughput single cell barcoding |
WO2012055929A1 (en) | 2010-10-26 | 2012-05-03 | Illumina, Inc. | Sequencing methods |
EP2635679B1 (en) | 2010-11-05 | 2017-04-19 | Illumina, Inc. | Linking sequence reads using paired code tags |
US9193997B2 (en) | 2010-12-15 | 2015-11-24 | The Board Of Trustees Of The Leland Stanford Junior University | Measuring and monitoring of cell clonality |
US20140057799A1 (en) | 2010-12-16 | 2014-02-27 | Gigagen | System and Methods for Massively Parallel Analysis of Nucleic Acids in Single Cells |
US9594870B2 (en) | 2010-12-29 | 2017-03-14 | Life Technologies Corporation | Time-warped background signal for sequencing-by-synthesis operations |
US10146906B2 (en) | 2010-12-30 | 2018-12-04 | Life Technologies Corporation | Models for analyzing data from sequencing-by-synthesis operations |
JP6211931B2 (en) | 2011-03-09 | 2017-10-11 | セル・シグナリング・テクノロジー・インコーポレイテツド | Methods and reagents for producing monoclonal antibodies |
US8759036B2 (en) | 2011-03-21 | 2014-06-24 | Affymetrix, Inc. | Methods for synthesizing pools of probes |
US9476095B2 (en) | 2011-04-15 | 2016-10-25 | The Johns Hopkins University | Safe sequencing system |
SI3415619T1 (en) | 2011-04-28 | 2021-04-30 | The Board of Trustees of the Leland Stanford Junior University Office of the General Counsel Building 170, Third Floor, Main Quad | Identification of polynucleotides associated with a sample |
SI3892295T1 (en) | 2011-05-24 | 2023-09-29 | BioNTech SE | Individualized vaccines for cancer |
WO2013033721A1 (en) | 2011-09-02 | 2013-03-07 | Atreca, Inc. | Dna barcodes for multiplexed sequencing |
WO2013036459A2 (en) | 2011-09-09 | 2013-03-14 | Sequenta, Inc. | Sequence-based measures of immune response |
US10385475B2 (en) | 2011-09-12 | 2019-08-20 | Adaptive Biotechnologies Corp. | Random array sequencing of low-complexity libraries |
WO2013055595A1 (en) | 2011-10-11 | 2013-04-18 | Sequenta, Inc. | Determining responsiveness of autoimmune patients to dmard treatment |
EP2768982A4 (en) | 2011-10-21 | 2015-06-03 | Adaptive Biotechnologies Corp | Quantification of adaptive immune cell genomes in a complex mixture of cells |
US20130137108A1 (en) | 2011-10-25 | 2013-05-30 | Anubhav Tripathi | Magnetic bead separation apparatus and method |
EP2773773B1 (en) | 2011-11-04 | 2017-01-11 | Adaptive Biotechnologies Corporation | T-cell receptor clonotypes shared among ankylosing spondylitis patients |
WO2013085855A1 (en) | 2011-12-05 | 2013-06-13 | Sequenta, Inc. | Clonotypes as biometric specimen tags |
EP3388535B1 (en) | 2011-12-09 | 2021-03-24 | Adaptive Biotechnologies Corporation | Diagnosis of lymphoid malignancies and minimal residual disease detection |
EP2788507A4 (en) | 2011-12-09 | 2015-07-08 | Sequenta Inc | Method of measuring immune activation |
US9499865B2 (en) | 2011-12-13 | 2016-11-22 | Adaptive Biotechnologies Corp. | Detection and measurement of tissue-infiltrating lymphocytes |
US20150031553A1 (en) | 2011-12-13 | 2015-01-29 | Sequenta, Inc. | Method of measuring immune activation |
JP2015504661A (en) | 2011-12-20 | 2015-02-16 | シーケンタ インコーポレイテッド | Monitoring of transformation from follicular lymphoma to diffuse large B-cell lymphoma by immunorepertoire analysis |
US20150031555A1 (en) | 2012-01-24 | 2015-01-29 | Gigagen, Inc. | Method for correction of bias in multiplexed amplification |
CA3133249C (en) | 2012-02-09 | 2023-07-25 | Memed Diagnostics Ltd. | Signatures and determinants for diagnosing infections and methods of use thereof |
EP3854873A1 (en) | 2012-02-17 | 2021-07-28 | Fred Hutchinson Cancer Research Center | Compositions and methods for accurately identifying mutations |
US11177020B2 (en) | 2012-02-27 | 2021-11-16 | The University Of North Carolina At Chapel Hill | Methods and uses for molecular tags |
WO2013131074A1 (en) | 2012-03-02 | 2013-09-06 | Diogenix, Inc. | Methods and reagents for evaluating autoimmune disease and determining antibody repertoire |
US20150038346A1 (en) | 2012-03-05 | 2015-02-05 | Sequenta, Inc. | Monitoring immune responsiveness to cancer vaccination |
JP6302847B2 (en) | 2012-03-05 | 2018-03-28 | アダプティヴ バイオテクノロジーズ コーポレーション | Determination of paired immunoreceptor chains from frequency matched subunits |
WO2013134302A1 (en) | 2012-03-05 | 2013-09-12 | Sequenta, Inc. | Monitoring immune responsiveness to cancer vaccination |
CA2867293C (en) | 2012-03-13 | 2020-09-01 | Abhijit Ajit PATEL | Measurement of nucleic acid variants using highly-multiplexed error-suppressed deep sequencing |
CN104395481A (en) | 2012-04-13 | 2015-03-04 | 赛昆塔公司 | Detection and quantitation of sample contamination in immune repertoire analysis |
WO2013158936A1 (en) | 2012-04-20 | 2013-10-24 | Sequenta, Inc | Monitoring immunoglobulin heavy chain evolution in b-cell acute lymphoblastic leukemia |
AU2013251701A1 (en) | 2012-04-24 | 2014-10-30 | Gen9, Inc. | Methods for sorting nucleic acids and multiplexed preparative in vitro cloning |
RU2631797C2 (en) | 2012-05-08 | 2017-09-26 | Эдэптив Байотекнолоджиз Корпорейшн | Compositions and methods of measurement and calibration of systematic mistake of amplification in multiplex pcr-reactions |
US20150247201A1 (en) | 2012-05-31 | 2015-09-03 | Malek Faham | Predicting relapse of chronic lymphocytic leukemia patients treated by allogeneic stem cell transplantation |
US20130324422A1 (en) | 2012-06-04 | 2013-12-05 | Sequenta, Inc. | Detecting disease-correlated clonotypes from fixed samples |
EP2859123A4 (en) | 2012-06-11 | 2015-12-16 | Sequenta Inc | Method of sequence determination using sequence tags |
WO2013188831A1 (en) | 2012-06-15 | 2013-12-19 | Adaptive Biotechnologies Corporation | Uniquely tagged rearranged adaptive immune receptor genes in a complex gene set |
WO2013192570A1 (en) | 2012-06-21 | 2013-12-27 | Gigagen, Inc. | System and methods for genetic analysis of mixed cell populations |
SG11201500313YA (en) | 2012-07-24 | 2015-02-27 | Sequenta Inc | Single cell analysis using sequence tags |
WO2014026031A1 (en) | 2012-08-10 | 2014-02-13 | Sequenta, Inc. | High sensitivity mutation detection using sequence tags |
US20140065629A1 (en) | 2012-08-29 | 2014-03-06 | Israel Barken | Methods of treating diseases |
CN114085896A (en) | 2012-09-24 | 2022-02-25 | Cb生物技术公司 | Multiplex pyrosequencing using non-interfering noise-canceling polynucleotides identification tags |
AU2013331135A1 (en) | 2012-10-19 | 2015-05-07 | Adaptive Biotechnologies Corp. | Monitoring diffuse large B-cell lymphoma from peripheral blood samples |
WO2014062945A1 (en) | 2012-10-19 | 2014-04-24 | Sequenta, Inc. | Monitoring clonotypes of plasma cell proliferative disorders in peripheral blood |
JP2015532123A (en) | 2012-10-22 | 2015-11-09 | シーケンタ・インコーポレイテッド | Monitoring of clonal type of mantle cell lymphoma in peripheral blood after immune transplantation |
AU2014218910A1 (en) | 2013-02-22 | 2015-08-20 | Adaptive Biotechnologies Corp. | Rare clonotypes and uses thereof |
US20140255944A1 (en) | 2013-03-08 | 2014-09-11 | Sequenta, Inc. | Monitoring treatment-resistant clones in lymphoid and myeloid neoplasms by relative levels of evolved clonotypes |
US20140255929A1 (en) | 2013-03-11 | 2014-09-11 | Sequenta, Inc. | Mosaic tags for labeling templates in large-scale amplifications |
SG10201707394PA (en) | 2013-03-15 | 2017-10-30 | Adaptive Biotechnologies Corp | Uniquely tagged rearranged adaptive immune receptor genes in a complex gene set |
US9708657B2 (en) | 2013-07-01 | 2017-07-18 | Adaptive Biotechnologies Corp. | Method for generating clonotype profiles using sequence tags |
US20160186260A1 (en) | 2013-07-26 | 2016-06-30 | Sequenta, Llc | Cancer vaccination with antigen evolution |
AU2014337063A1 (en) | 2013-10-18 | 2016-05-19 | Adaptive Biotechnologies Corp. | Predicting patient responsiveness to immune checkpoint inhibitors |
US20160340729A1 (en) | 2014-01-10 | 2016-11-24 | Adaptive Biotechnologies Corp. | Methods for defining and predicting immune response to allograft |
US20150218656A1 (en) | 2014-02-03 | 2015-08-06 | Adaptive Biotechnologies Corp. | Methods for detection and diagnosis of a lymphoid malignancy using high throughput sequencing |
CA2941612A1 (en) | 2014-03-05 | 2015-09-11 | Adaptive Biotechnologies Corporation | Methods using randomer-containing synthetic molecules |
US10066265B2 (en) | 2014-04-01 | 2018-09-04 | Adaptive Biotechnologies Corp. | Determining antigen-specific t-cells |
EP3132059B1 (en) | 2014-04-17 | 2020-01-08 | Adaptive Biotechnologies Corporation | Quantification of adaptive immune cell genomes in a complex mixture of cells |
WO2016069886A1 (en) | 2014-10-29 | 2016-05-06 | Adaptive Biotechnologies Corporation | Highly-multiplexed simultaneous detection of nucleic acids encoding paired adaptive immune receptor heterodimers from many samples |
US10246701B2 (en) | 2014-11-14 | 2019-04-02 | Adaptive Biotechnologies Corp. | Multiplexed digital quantitation of rearranged lymphoid receptors in a complex mixture |
EP3262196B1 (en) | 2015-02-24 | 2020-12-09 | Adaptive Biotechnologies Corp. | Method for 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 |
-
2016
- 2016-02-24 EP EP16756268.5A patent/EP3262196B1/en active Active
- 2016-02-24 EP EP19188465.9A patent/EP3591074A1/en active Pending
- 2016-02-24 AU AU2016222788A patent/AU2016222788B2/en active Active
- 2016-02-24 ES ES16756268T patent/ES2858306T3/en active Active
- 2016-02-24 US US15/553,040 patent/US11047008B2/en active Active
- 2016-02-24 CA CA2976580A patent/CA2976580A1/en not_active Abandoned
- 2016-02-24 WO PCT/US2016/019343 patent/WO2016138122A1/en active Application Filing
-
2021
- 2021-06-09 US US17/343,552 patent/US20210381050A1/en active Pending
-
2022
- 2022-06-29 AU AU2022204608A patent/AU2022204608A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
WO2016138122A1 (en) | 2016-09-01 |
EP3262196A4 (en) | 2018-11-21 |
EP3262196A1 (en) | 2018-01-03 |
EP3262196B1 (en) | 2020-12-09 |
ES2858306T3 (en) | 2021-09-30 |
US11047008B2 (en) | 2021-06-29 |
US20180037953A1 (en) | 2018-02-08 |
AU2022204608A1 (en) | 2022-07-21 |
AU2016222788B2 (en) | 2022-03-31 |
US20210381050A1 (en) | 2021-12-09 |
AU2016222788A1 (en) | 2017-08-31 |
EP3591074A1 (en) | 2020-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210381050A1 (en) | Methods for diagnosing infectious disease and determining hla status using immune repertoire sequencing | |
Norman et al. | Sequences of 95 human MHC haplotypes reveal extreme coding variation in genes other than highly polymorphic HLA class I and II | |
AU2016242967B2 (en) | Method of identifying human compatible T cell receptors specific for an antigenic target | |
Heather et al. | Dynamic perturbations of the T-cell receptor repertoire in chronic HIV infection and following antiretroviral therapy | |
AU2014232314B2 (en) | Uniquely tagged rearranged adaptive immune receptor genes in a complex gene set | |
EP2719774B1 (en) | Methods of monitoring conditions by sequence analysis | |
US20150299786A1 (en) | Uniquely tagged rearranged adaptive immune receptor genes in a complex gene set | |
WO2011106738A2 (en) | Use of tcr clonotypes as biomarkers for disease | |
WO2015106161A1 (en) | Methods for defining and predicting immune response to allograft | |
US11066705B2 (en) | Characterization of adaptive immune response to vaccination or infection using immune repertoire sequencing | |
JP2016516449A5 (en) | ||
US20220259670A1 (en) | Kit and method for analyzing single t cells | |
Emerson et al. | Immunosequencing reveals diagnostic signatures of chronic viral infection in T cell memory | |
Wagner et al. | Recruitment of highly functional SARS-CoV-2-specific CD8+ T cell receptors mediating cytotoxicity of virus-infected target cells in non-severe COVID-19 | |
Wei et al. | DRB genotyping in cynomolgus monkeys from China using polymerase chain reaction-sequence-specific primers | |
Tu | Recovery of T cell receptor variable sequences from 3'barcoded single-cell RNA sequencing libraries | |
EP4314337A1 (en) | Immune cell counting of sars-cov-2 patients based on immune repertoire sequencing | |
Tschochner et al. | Identifying Epstein-Barr virus EBNA-1 sequence variation using 454 FLX technology in Multiple Sclerosis patient samples | |
Nikula | Transcriptional Profiling of Organ‐Specific Autoimmunity in Human |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FZDE | Discontinued |
Effective date: 20220517 |
|
FZDE | Discontinued |
Effective date: 20220517 |