The present disclosure relates generally to methods for targeted hybrid capture of rearranged T cell receptors. More particularly, some embodiments relate to a method for direct and quantitative, error-corrected counting of genomic sequences. Some embodiments also relate to specific counts of T cell populations that are present in a sample.
T cells are integral mediators of the adaptive immune response in vertebrate organisms. They control the production of antibodies by co-stimulating B cells, and they mediate direct clearance of pathogen-infected and physiologically-defective cells by direct physical engagement between the T cell and the distressed target cell. The cell-to-cell interaction between T cells and targets is undeniably complex, yet central to the process is engagement of T cell receptors (TCRs) found on the surface of the T cell surface and major histocompatibility complex (MHC) molecules displayed on the surface of target cells. The genes encoding TCRs are assembled from a pre-existing array of possible gene segments that are present as germline sequences in all cells. During T cell development, this array is assembled by site-specific recombinases into potential T cell receptor sequences (TCRs). Those cells that produce a functional TCR that does not recognize self eventually mature and become part an individual's T cell repertoire.
The introduction of therapies that rely on the stimulation of innate T cells to treat cancers has garnered well-deserved attention. Some treated patients have experienced complete and durable responses for disease indications that previously had dismal survival prognoses. The current goal of clinical research is to understand how these T cells become activated. Similarly, in the context of clinical therapy there remains a need to determine if and when efficacious T cell populations become mobilized in the eradication of cancerous tissues.
It is therefore an aspect of this disclosure to provide methods for profiling adaptive immune response genes in a sample.
Some embodiments provided herein relate to methods of identifying a rearranged adaptive immune response gene. In some embodiments, the method comprises: obtaining a sample comprising genomic DNA; isolating genomic DNA from the sample; capturing a rearranged adaptive immune response gene from the isolated genomic DNA by sequential hybridization; amplifying the second extended sequence; and/or sequencing the second extended sequence. In some embodiments, the sequential hybridization comprises: hybridizing the genomic DNA with a first set of probes specific to a first portion of the rearranged adaptive immune response gene to generate a hybridized sequence; extending the first set of probes to generate a first extended sequence; purifying or isolating the first extended sequence; hybridizing the purified first extended sequence with a second set of probes specific to a second portion of the rearranged adaptive immune response gene; and/or extending the second set of probes to generate a second extended sequence.
In some embodiments, the sample is obtained from a tissue or a biofluid. In some embodiments, the sample is obtained from a tumor tissue, a region proximal to a tumor tissue, an organ tissue, peripheral tissue, lymph, urine, cerebral spinal fluid, a buffy coat isolate, whole blood, peripheral blood, bone marrow, amniotic fluid, breast milk, plasma, serum, aqueous humor, vitreous humor, cochlear fluid, saliva, stool, sweat, vaginal secretions, semen, bile, tears, mucus, sputum, and/or vomit. In some embodiments, the sample comprises adaptive immune cells. In some embodiments, the sample comprises one or more immune cells, such as T cells.
In some embodiments, the rearranged adaptive immune response gene is encoded by the T cell receptor (TCR) alpha gene (TRA), the TCR beta gene (TRB), the TCR delta gene (TRD), the TCR gamma gene (TRG), the antibody heavy chain gene (IGH), the kappa light chain antibody gene (IGK), and/or the lambda light chain antibody gene (IGL).
In some embodiments, the first portion of the rearranged adaptive immune response gene is a CDR3-encoding region, comprising a V, D, or J region of the rearranged adaptive immune response gene. In some embodiments, the first extended sequence is copied with T4 DNA polymerase and T4 gene 32 protein.
In some embodiments, extending is performed in a solution containing polyethylene glycol (PEG). In some embodiments, the PEG has an average molecular weight of 8000 Daltons (PEG8000). In some embodiments, PEG is present in an amount of 2-40% w/v, such as 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10, 15, 20, 25, 30, 35, or 40% w/v, or an amount within a range defined by any two of the aforementioned values.
In some embodiments, the method further comprises fragmenting and end-repairing the genomic DNA prior to sequential hybridization. In some embodiments, the method further comprises ligating an amplification adaptor to the first extended sequence. In some embodiments, the amplifying is performed by polymerase chain reaction (PCR).
In some embodiments, the first set of probes comprises J region sequences of human TCR alpha (TRA), human TCR beta (TRB), human TCR gamma (TRG), human TCR delta (TRG), a human antibody heavy chain (IGH), a human kappa light chain antibody (IGK), and/or a human lambda light chain antibody (IGL). In some embodiments, the first set of probes comprises V region sequences of human TRA, human TRB, human TRG, human TRD, human IGH, human IGK, and/or human IGL. In some embodiments, the second set of probes comprises J region sequences of human TRA, human TRB, human TRG, human TRD, human IGH, human IGK, and/or human IGL. In some embodiments, the second set of probes comprises V region sequences of human TRA, human TRB, human TRG, human TRD, human IGH, human IGK, and/or human IGL.
In some embodiments, the first set of probes comprises a DNA sequence tag for identification of specific clones. In some embodiments, the DNA sequence tag is a nucleic acid sequence from including 2-10 nucleic acids, such as 2, 3, 4, 5, 6, 7, 8, 9, or 10 nucleic acids selected at random. In some embodiments, the DNA sequence tag includes a sequence of NN, NNN, NNNN, NNNNN, NNNNNN, NNNNNNN, NNNNNNNN, NNNNNNNNN, or NNNNNNNNNN, wherein N is A, T, G, or C. In some embodiments, the DNA sequence tags, the first and second set of probes, and the captured sequences are all used in informatic identification of clones. In some embodiments, the sample comprises a plurality of rearranged genomic sequences.
In some embodiments, the method further comprises determining the frequency of specific T cell clones, B cell clones, or both in the sample to determine a T cell immune repertoire, a B cell repertoire, or both in the sample. In some embodiments, the method further comprises profiling circulating nucleic acids, TCR repertoire, and/or Ab repertoire in a whole blood sample. In some embodiments, the profiling comprises a determination of the characteristics of a population of nucleic acids, TCR repertoire, and/or Ab repertoire in a sample.
In some embodiments, the method further comprises assessing both circulating nucleic acid and immune repertoire from a single whole blood sample. In some embodiments, an amount of single cell genomic DNA is increased by whole genome amplification prior to analysis. In some embodiments, single cell analysis is used to identify pairing between alpha and beta chain TCR within a single cell. In some embodiments, the first set of probes comprises a nucleic acid having at least 90% sequence identity to any sequence defined by any one or more of SEQ ID NOs: 62-128. In some embodiments, the second set of probes comprises a nucleic acid having at least 90% sequence identity to any sequence defined by any one or more of SEQ ID NOs: 129-227.
Embodiments provided herein relate to methods for profiling adaptive immune response genes in a sample, including determination of adaptive immune response gene repertoires in a sample.
TCRs are a unique signature for each T cell, and therefore the determination of TCR repertoires provides direct insight into the activities of the adaptive immune response. There are several other clinical applications of TCR profiling that include minimal residual disease monitoring in T cell lymphomas, individual response to vaccines meant to stimulate the adaptive immune system, and adaptive immune responses to infectious diseases.
As shown in
Methods have been described in which TCR-specific PCR primers are used amplify and sequence rearranged TCR segments from genomic DNA (Robins H, Desmarais C, Matthis J, Livingston R, Andriesen J, Reijonen H, et al. Ultra-sensitive detection of rare T cell clones. J Immunol Methods. 2012 Jan. 31; 375(1-2):14-9, expressly incorporated herein by reference in its entirety). Several commercially-available methods take advantage of the fact that rearranged TCR are expressed as messenger RNAs, and they use RNA-seq methods to monitor TCR repertoires (e.g. Immunoverse from Archer Dx, Immune repertoire-seq from CD-Genomics, Full-Length V(D)J Sequences from 10× genomics). The use of molecular identifiers has been used to provide error-correction and a quantitative framework for analysis (Shugay M, Britanova O V, Merzlyak E M, Turchaninova M A, Mamedov I Z, Tuganbaev T R, et al. Towards error-free profiling of immune repertoires. Nat Methods. 2014 June; 11(6):653-5, expressly incorporated herein by reference in its entirety). Both genomic PCR and mRNA profiling, even with molecular tags, are indirect measurements of T cell repertoires. The genomic methods rely on multiplex PCR and are subject to amplification biases. Moreover, they lack error-correcting strategies and are therefore prone to over-estimates of TCR diversity. Expression-based methods measure TCR expression levels rather than T cell populations, and the well-established observation that TCR expression is governed by T cell activation (Paillard F, Sterkers G, and Vaquero C. Transcriptional and post-transcriptional regulation of TCR, CD4 and CD8 gene expression during activation of normal human T lymphocytes. EMBO J. 1990 June; 9(6): 1867-1872, expressly incorporated herein by reference in its entirety) is likely to provide a distorted view of T cell populations. This is a particularly critical consideration in the context of oncology where the efficacy of immune checkpoint inhibitors relies on a pre-existing population of inactive but potentially responsive tumor-specific killer T cells.
Some embodiments provided herein relate to a method to tag, retrieve, and/or quantify TCR repertoires. The next generation sequencing (NGS) readout is an accurate census of T cells that are present in the analysis sample. The method utilizes targeted hybrid capture technology. In the current context, tagged capture probes are used to retrieve and copy one of the partner gene segments that is rearranged to a functional TCR gene in T cells. Notably, this first capture step captures all possible gene segments, including the vast majority that is not rearranged in cells other than T cells. In a second capture step, probes specific for the other partner gene segment, which are brought in close proximity to the first partner during TCR gene development, are used to retrieve rearranged TCR genes from the initial library. In some embodiments, the method of using two capture steps is referred to herein as “sequential capture.” In some embodiments, this method provides readouts of the highly-diverse, antigen-binding CDR3 regions as a signature of individual T cells. Importantly, the TCR repertoires collected from within one individual over short periods of time may be highly similar while the repertoires collected from different individuals may differ substantially. In some embodiments, the method is both reproducible and specific.
In some embodiments, sequential capture (e.g., comprising the aforementioned two capture steps) may be used for determination of adaptive immune response gene repertoires of adaptive immune systems that undergo gene rearrangements. In some embodiments, for example, sequential capture may be used with TCR alpha and TCR beta gene targets for determination of TCR repertoires. However, the methods described herein may be used on other targets, such as other TCRs (e.g. gamma and delta chains) present on T cells that generally inhabit the digestive system. Antibody-producing B cells also possess repertoires of genes produced by genomic rearrangement. In some embodiments, methods described herein are applicable to profiling of these cell populations as well.
In some embodiments, the method of immune repertoire profiling is conducted on circulating alpha and beta chain bearing T cells. In some embodiments, the method of immune repertoire profiling is conducted on antibody producing B cells and gastric T cell delta gamma repertoires. In some embodiments, the method of immune repertoire profiling is nucleic acid hybridization and capture based. Significantly, the methods described herein differ from other profiling methods, which are PCR based. The methods described herein may use PCR to amplify DNA, but “sequential hybridization” with a first probe to one end of the TCR gene (for example, the J region or the V region), enrichment of these clones, and a second probe for the other end of the TCR (J→V, or V→J) of the enriched clones differentiates the present disclosure from standard techniques.
In addition, in some embodiments, the method for immune repertoire profiling is a genomic method that interrogates genomic DNA. In contrast, other commercially available technologies rely on mRNA transcript analysis, where mRNA is converted to cDNA and then enriched by specific PCR primers. One problem with these standards techniques is that clinicians care about T cell populations rather than expression levels of TCRs. Another issue that these standard techniques present is inaccurate test results. By way of example, consider a system having two populations of T cells, where one population is fighting off an infection. This population would be transcribing TCR message at a furious rate. The other population can fight off cancer, but the tumor is down-regulating its response. This population is making TCR message in minute quantities. If the TCR repertoire is profiled based on messenger RNA, a false conclusion would be that there are far more infection fighting cells than cancer fighting cells, even though in reality they are equal populations.
Some embodiments provided herein relate to methods that quantitatively analyze or count individual T cell clones by introducing a tag at the first hybridization step. This tag persists throughout the hybridization, capture, and sequencing steps and is used in post-sequence analysis to count T cell clones. The methods provided herein are not amenable to standard PCR-based profiling methods.
In some embodiments, these tags serve a purpose of eliminating false TCR clones. Using PCR only, it is not possible to tell the difference between a true positive clone that is rare versus a false positive clone that is the result of an error, such as a sequencing error. These false positive clones are particularly troublesome in the face of next-generation sequencing that generates millions of sequences. With the significant amount of data that is generated, errors can create functional TCR sequences that were not actually present in the biological sample being analyzed. However, the methods described herein using tags allow for identification of related sequences that arise by post-sample, error-driven processes.
Quantitative analysis of T cell clones is important for profiling T cell repertoire, and changes thereof. For example, profiling the T cell repertoire before and after an immunotherapy administration is useful for monitoring efficacy during treatment. Without wishing to be bound by theory, but by way of example, many of the newest class of immunotherapies rely on stimulating a preexisting set of TCR clones that have been inactivated by immune checkpoint molecules, such as PD-L1. By blocking the influence of PD-L1 (for example, with monoclonal antibodies), it is possible to activate the anti-tumor T cell repertoire. The course of therapy can be followed by profiling the T cell repertoire before and after administration of the PD-L1 checkpoint inhibitor. The methods described herein are useful for monitoring efficacy during methods of therapy, such as methods of treatment or inhibition of diseases such as cancer, which is valuable because some tumors respond to activation and others do not.
While still not wishing to be bound by theory, each DNA:DNA hybridization reaction is independent of a different reaction that involves a different set of sequences. By extension, it is possible to conduct thousands of probe:genomic-target capture steps simultaneously within a single reaction vessel, as long as each reaction is a simple bimolecular complex. Still further, methods described herein, including the capture methods are capable of capturing and removing TCRs, Ab-producing genes, MHC genes, tumor-related cancer genes and other adaptive immune response genes in a single reaction. In contrast, PCR-based methods rely only on the specificity of a trimolecular hybridization in which the genomic fragment, the first primer, and the second primer all specifically interact on the same genomic sequence. PCR is a far more complex reaction because subtle interactions between highly concentrated PCR primers can dominate the hybridization outcome. Thus, multiplex PCR systems are very limited and complex. The hybridization-based methods described herein operate on fundamentally different principles than existing multiplex PCR methods.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as is commonly understood by one of ordinary skill in the art. All patents, applications, published applications and other publications referenced herein are expressly incorporated by reference in their entireties unless stated otherwise. In the event that there is a plurality of definitions for a term herein, those in this section prevail unless stated otherwise.
As used herein, the term “adaptive immune system” has its ordinary meaning as understood in light of the specification, and refers to highly specialized, systemic cells and processes that eliminate pathogenic challenges. The cells of the adaptive immune system are a type of leukocyte, called a lymphocyte. B cells and T cells are the major types of lymphocytes.
As used herein, the term “immune cell” has its ordinary meaning as understood in light of the specification, and refers to cells that play a role in the immune response. Immune cells are of hematopoietic origin, and include lymphocytes, such as B cells and T cells; natural killer (NK) cells; myeloid cells, such as monocytes, macrophages, eosinophils, mast cells, basophils, and/or granulocytes.
As used herein, the term “T cell” has its ordinary meaning as understood in light of the specification, and includes CD4+ T cells and CD8+ T cells. The term T cell also includes T helper 1 type T cells, T helper 2 type T cells, T helper 17 type T cells and/or inhibitory T cells. The term “antigen presenting cell” includes antigen presenting cells (e.g., B lymphocytes, monocytes, dendritic cells, and/or Langerhans cells), as well as, other antigen presenting cells (e.g., keratinocytes, endothelial cells, astrocytes, fibroblasts, and/or oligodendrocytes). Some embodiments provided herein relate to providing or administering T cells to subjects in need of an immune response. Some embodiments provided herein relate to profiling of T cell compartments. The sorting of T cells using surface-specific markers coupled to fluorescence-activated cell sorting is a fundamental technology in immunological research. As used herein, the term “T cell compartments” has its ordinary meaning as understood in light of the specification, and refers to specific sets of T cells that all have the same surface markers.
As used herein, the term “immune response” has its ordinary meaning as understood in light of the specification, and includes T cell mediated and/or B cell mediated immune responses that are influenced by modulation of T cell co-stimulation. Exemplary immune responses include T cell responses, e.g., cytokine production, and/or cellular cytotoxicity. In addition, the term immune response includes immune responses that are indirectly affected by T cell activation, e.g., antibody production (humoral responses) and/or activation of cytokine responsive cells, e.g., macrophages. In the adaptive immune response, antigens are recognized by hypervariable molecules, such as antibodies or TCRs, which are expressed with sufficiently diverse structures to be able to recognize any antigen. Antibodies can bind to any part of the surface of an antigen. TCRs, however, are restricted to binding to short peptides bound to class I or class II molecules of the major histocompatibility complex (MHC) on the surface of APCs. TCR recognition of a peptide/MHC complex triggers activation (clonal expansion) of the T cell.
As used herein, “T cell receptor (TCR)” has its ordinary meaning as understood in light of the specification, and refers to a T cell receptor or a T cell antigen receptor, or a receptor expressed on a cell membrane of a T cell that regulates an immune system, and recognizes an antigen. There are α chain, β chain, γ chain and δ chain, constituting an αβ or γδ dimer. A TCR consisting of the former combination is called an αβ TCR and a TCR consisting of the latter combination is called a γδ TCR. T cells having such TCRs are called αβ T cell or γδ T cell. The structure is very similar to a Fab fragment of an antibody produced by a B cell, and recognizes an antigen molecule bound to an MHC molecule. Since a TCR gene of a mature T cell has undergone gene rearrangement, an individual has a diverse TCR and is able to recognize various antigens. A TCR further binds to an invariable CD3 molecule present in a cell membrane to form a complex. CD3 has an amino acid sequence called the ITAM (immunoreceptor tyrosine-based activation motif) in an intracellular region. This motif is considered to be involved in intracellular signaling. Each TCR chain is composed of a variable section (V) and a constant section (C). The constant section penetrates through the cell membrane and has a short cytoplasm portion. The variable section is present extracellularly and binds to an antigen-MHC complex. The variable section has three regions called a hypervariable section or a complementarity determining region (CDR), which binds to an antigen-MHC complex. The three CDRs are each called CDR1, CDR2, and CDR3. For a TCR, CDR1 and CDR2 are considered to bind to an MHC, while CDR3 is considered to bind to an antigen. Gene rearrangement of a TCR is similar to the process for a B cell receptor known as an immunoglobulin. In gene rearrangement of an αβ TCR, VDJ rearrangement of a β chain is first performed and then VJ rearrangement of an α chain is performed. Since a gene of a δ chain is deleted from a chromosome in rearrangement of an α chain, a T cell having an αβ TCR would not simultaneously have a γδ TCR. In contrast, in a T cell having a γδ TCR, a signal mediated by this TCR suppresses expression of a β chain. Thus, a T cell having a γδ TCR would not simultaneously have an αβ TCR.
As used herein, “B cell receptor (BCR)” has its ordinary meaning as understood in light of the specification, and is also called a B cell receptor or B cell antigen receptor and refers to those composed of an Igα/Igβ (CD79a/CD79b) heterodimer (α/β) conjugated with a membrane-bound immunoglobulin (mIg). An mIg subunit binds to an antigen to induce aggregation of the receptors, while an α/β subunit transmits a signal to the inside of a cell. BCRs, when aggregated, are understood to quickly activate Lyn, Blk, and Fyn of Src family kinases as in Syk and Btk of tyrosine kinases. Results greatly differ depending on the complexity of BCR signaling, the results including survival, resistance (allergy; lack of hypersensitivity reaction to antigen) or apoptosis, cell division, differentiation into antibody-producing cell or memory B cell and the like. Several hundred million types of T cells with a different TCR variable region sequence are produced and several hundred million types of B cells with a different BCR (or antibody) variable region sequence are produced. Individual sequences of TCRs and BCRs vary due to an introduced mutation or rearrangement of the genomic sequence. Thus, it is possible to obtain a clue for antigen specificity of a T cell or a B cell by determining a genomic sequence of TCR/BCR or a sequence of an mRNA (cDNA).
As used herein, “V region” has its ordinary meaning as understood in light of the specification, and refers to a variable section (V) of a variable region of a TCR chain or a BCR chain. As used herein, “D region” has its ordinary meaning as understood in light of the specification, and refers to a D region of a variable region of a TCR chain or a BCR chain. As used herein, “J region” has its ordinary meaning as understood in light of the specification, and refers to a J region of a variable region of a TCR chain or a BCR chain. As used herein, “C region” has its ordinary meaning as understood in light of the specification, and refers to a constant section (C) region of a TCR chain or a BCR chain.
The combinatorial joining of V and J segments in α chains and V, D and J segments in β chains produces a large number of possible molecules, thereby creating a diversity of TCRs. Diversity is also achieved in TCRs by alternative joining of gene segments. In contrast to Ig, β and δ gene segments can be joined in alternative ways. RSS flanking gene segments in β and δ gene segments can generate VJ and VDJ in the β chain, and VJ, VDJ, and VDDJ on the δ chain. As in the case of Ig, diversity is also produced by variability in the joining of gene segments. Some embodiments provided herein relate to gene segments, including T cell receptor alpha chain V region (TRAV), T cell receptor beta chain V region (TRBV) T cell receptor alpha chain J region (TRAJ), or T cell receptor beta chain J region (TRBJ).
In some embodiments, adaptive immune response genes may include TCR alpha gene (TRA), the TCR beta gene (TRB), the TCR delta gene (TRD), the TCR gamma gene (TRG), the antibody heavy chain gene (IGH), the kappa light chain antibody gene (IGK), and/or the lambda light chain antibody gene (IGL).
As used herein, the term “rearranged” has its ordinary meaning as understood in light of the specification, and refers to a configuration of a heavy chain or light chain immunoglobulin locus wherein a V segment is positioned immediately adjacent to a D-J or J segment in a conformation encoding essentially a complete VH and VL domain, respectively. A rearranged immunoglobulin gene locus can be identified by comparison to germline DNA; a rearranged locus will have at least one recombined heptamer/nonamer homology element.
As used herein, the term “unrearranged” or “germline configuration” in reference to a V segment has its ordinary meaning as understood in light of the specification, and refers to the configuration wherein the V segment is not recombined so as to be immediately adjacent to a D or J segment.
The term “gene” has its ordinary meaning as understood in light of the specification, and includes the segment of DNA involved in producing a polypeptide chain. Specifically, a gene includes, without limitation, regions preceding and following the coding region, such as the promoter and 3′-untranslated region, respectively, as well as intervening sequences (introns) between individual coding segments (exons). As used herein, “genomic DNA” refers to chromosomal DNA, as opposed to complementary DNA copied from an RNA transcript. “Genomic DNA”, as used herein, may be all of the DNA present in a single cell, or may be a portion of the DNA in a single cell.
The term “nucleic acid” or “polynucleotide” has its ordinary meaning as understood in light of the specification, and includes deoxyribonucleotides or ribonucleotides and polymers thereof in either single- or double-stranded form. Unless specifically limited, the term encompasses nucleic acids containing known analogues of natural nucleotides that have similar binding properties as the reference nucleic acid and are metabolized in a manner similar to naturally occurring nucleotides. Unless otherwise indicated, a particular nucleic acid sequence also implicitly encompasses conservatively modified variants thereof (e.g., degenerate codon substitutions), alleles, orthologs, SNPs, and complementary sequences as well as the sequence explicitly indicated. Specifically, degenerate codon substitutions may be achieved by generating sequences in which the third position of one or more selected (or all) codons is substituted with mixed-base and/or deoxyinosine residues (Batzer et al., Nucleic Acid Res., 19:5081 (1991); Ohtsuka et al., J. Biol. Chem., 260:2605-2608 (1985); Rossolini et al., Mol. Cell. Probes, 8:91-98 (1994), each of which is expressly incorporated herein by reference in its entirety). The term nucleic acid is used interchangeably with gene, cDNA, and mRNA encoded by a gene.
As used herein, the terms “nucleic acid” and “polynucleotide” are interchangeable and has its ordinary meaning as understood in light of the specification, and refer to any nucleic acid, whether composed of phosphodiester linkages or modified linkages such as phosphotriester, phosphoramidate, siloxane, carbonate, carboxymethylester, acetamidate, carbamate, thioether, bridged phosphoramidate, bridged methylene phosphonate, bridged phosphoramidate, bridged phosphoramidate, bridged methylene phosphonate, phosphorothioate, methylphosphonate, phosphorodithioate, bridged phosphorothioate and/or sulfone linkages, or combinations of such linkages. The terms “nucleic acid” and “polynucleotide” has its ordinary meaning as understood in light of the specification, and also specifically include nucleic acids composed of bases other than the five biologically occurring bases (adenine, guanine, thymine, cytosine and uracil).
As used herein, the term “antibody” has its ordinary meaning as understood in light of the specification, and includes whole antibodies and any antigen binding fragment (i.e., “antigen-binding portion”) or single chain thereof. An “antibody” refers to a glycoprotein comprising at least two heavy (H) chains and two light (L) chains inter-connected by disulfide bonds, or an antigen binding portion thereof. Each heavy chain is comprised of a heavy chain variable region (abbreviated herein as VH) and a heavy chain constant region. The heavy chain constant region is comprised of three domains, CH1, CH2 and CH3. Each light chain is comprised of a light chain variable region (abbreviated herein as VL) and a light chain constant region. The light chain constant region is comprised of one domain, CL. The VH and VL regions can be further subdivided into regions of hypervariability, termed complementarity determining regions (CDR), interspersed with regions that are more conserved, termed framework regions (FR). Each VH and VL is composed of three CDRs and four FRs, arranged from amino-terminus to carboxy-terminus in the following order: FR1, CDR1, FR2, CDR2, FR3, CDR3, and FR4. The variable regions of the heavy and light chains contain a binding domain that interacts with an antigen.
As used herein, “CDR3” has its ordinary meaning as understood in light of the specification, and refers to the third complementarity-determining region (CDR). In this regard, CDR is a region that directly contacts an antigen and undergoes a particularly large change among variable regions, and is referred to as a hypervariable region. Each variable region of a light chain and a heavy chain has three CDRs (CDR1-CDR3) and 4 FRs (FR1-FR4) surrounding the three CDRs. Because a CDR3 region is considered to be present across V region, D region and J region, it is considered as an important key for a variable region, and is thus used as a subject of analysis. As used herein, “front of CDR3 on a reference V region” refers to a sequence corresponding to the front of CDR3 in a V region targeted by the present disclosure. As used herein, “end of CDR3 on a reference J” refers to a sequence corresponding to the end of CDR3 in a J region targeted by the present disclosure.
As used herein, the term “antigen-binding portion” of an antibody (or simply “antibody portion”), has its ordinary meaning as understood in light of the specification, and refers to one or more fragments of an antibody that retain the ability to specifically bind to an antigen (e.g., PD-1, PD-L1, and/or PD-L2). It has been shown that the antigen-binding function of an antibody can be performed by fragments of a full-length antibody. Examples of binding fragments encompassed within the term “antigen-binding portion” of an antibody include (i) a Fab fragment, a monovalent fragment consisting of the VH, VL, CL and CH1 domains; (ii) a F(ab′)2fragment, a bivalent fragment comprising two Fab fragments linked by a disulfide bridge at the hinge region; (iii) a Fd fragment consisting of the VH and CH1 domains; (iv) a Fv fragment consisting of the VH and VL domains of a single arm of an antibody, (v) a dAb fragment, which consists of a VH domain; and (vi) an isolated complementarity determining region (CDR) or (vii) a combination of two or more isolated CDRs which may optionally be joined by a synthetic linker.
As used herein, the term “variant” has its ordinary meaning as understood in light of the specification, and refers to a polynucleotide (or polypeptide) having a sequence substantially similar to a reference polynucleotide (or polypeptide). In the case of a polynucleotide, a variant can have deletions, substitutions, additions of one or more nucleotides at the 5′ end, 3′ end, and/or one or more internal sites in comparison to the reference polynucleotide. Similarities and/or differences in sequences between a variant and the reference polynucleotide can be detected using conventional techniques known in the art, for example polymerase chain reaction (PCR) and hybridization techniques. Variant polynucleotides also include synthetically derived polynucleotides, such as those generated, for example, by using site-directed mutagenesis. Generally, a variant of a polynucleotide, including, but not limited to, a DNA, can have at least 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, or 99% or more sequence identity to the reference polynucleotide as determined by sequence alignment programs known by skilled artisans. In the case of a polypeptide, a variant can have deletions, substitutions, additions of one or more amino acids in comparison to the reference polypeptide. Similarities and/or differences in sequences between a variant and the reference polypeptide can be detected using conventional techniques known in the art, for example Western blot. Generally, a variant of a polypeptide, can have at least 60%, 65%, 70%, 75%, 80%, 85%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99% or more sequence identity to the reference polypeptide as determined by sequence alignment programs known by skilled artisans.
As used herein, the term “profile” has its ordinary meaning as understood in light of the specification, and includes any set of data that represents the distinctive features or characteristics associated with a tumor, tumor cell, and/or cancer. The term encompasses a “nucleic acid profile” that analyzes one or more genetic markers, a “protein profile” that analyzes one or more biochemical or serological markers, and combinations thereof. Examples of nucleic acid profiles include, but are not limited to, a genotypic profile, gene copy number profile, gene expression profile, DNA methylation profile, and combinations thereof. Non-limiting examples of protein profiles include a protein expression profile, protein activation profile, and combinations thereof. For example, a “genotypic profile” includes a set of genotypic data that represents the genotype of one or more genes associated with a tumor, tumor cell, and/or cancer. Similarly, a “gene copy number profile” includes a set of gene copy number data that represents the amplification of one or more genes associated with a tumor, tumor cell, and/or cancer. Likewise, a “gene expression profile” includes a set of gene expression data that represents the mRNA levels of one or more genes associated with a tumor, tumor cell, and/or cancer. In addition, a “DNA methylation profile” includes a set of methylation data that represents the DNA methylation levels (e.g., methylation status) of one or more genes associated with a tumor, tumor cell, and/or cancer. Furthermore, a “protein expression profile” includes a set of protein expression data that represents the levels of one or more proteins associated with a tumor, tumor cell, and/or cancer. Moreover, a “protein activation profile” includes a set of data that represents the activation (e.g., phosphorylation status) of one or more proteins associated with a tumor, tumor cell, and/or cancer.
As used herein, “repertoire of a variable region” refers to a collection of V(D)J regions created in any manner by gene rearrangement in a TCR or BCR. The terms such as TCR repertoire and BCR repertoire are used, which are also called, for example, T cell repertoire, B cell repertoire or the like in some cases. For instance, “T cell repertoire” refers to a collection of lymphocytes characterized by expression of a T cell receptor (TCR) serving an important role in antigen recognition. A change in a T cell repertoire provides a significant indicator of an immune status in a physiological condition and disease condition. In some embodiments provided herein, a repertoire determination may include determination of a T cell immune repertoire, a B cell repertoire, circulating nucleic acids repertoire, TCR repertoire, and/or Ab repertoire.
The term “identifying” has its ordinary meaning as understood in light of the specification, and refers to assessing, determining, or ascertaining the presence, absence, identity, quality, and/or quantity of an endpoint of interest. For example, identifying a rearranged adaptive immune response gene may refer to a determination of the presence and/or quantity of an adaptive immune response gene in a sample, including a determination of the identity of the adaptive immune response gene.
The term “sample” has its ordinary meaning as understood in light of the specification, and includes any biological specimen obtained from a subject. Samples include, without limitation, a biofluid, whole blood, peripheral blood, plasma, serum, red blood cells, white blood cells (e.g., peripheral blood mononuclear cells), saliva, urine, stool, sweat, tears, vaginal secretions, nipple aspirate, amniotic fluid, breast milk, semen, bile, mucus, sputum, vomit, lymph, fine needle aspirate, cerebrospinal fluid, a buffy coat isolate, aqueous humor, vitreous humor, cochlear fluid, any other bodily fluid, bone marrow, a tissue sample, a tumor tissue, a region proximal to a tumor tissue, an organ tissue, peripheral tissue, and/or cellular extracts thereof. In some embodiments, the sample is whole blood or a fractional component thereof such as plasma, serum, or a cell pellet.
Each T cell has a unique T cell receptor (TCR). The TCRs are protein dimers on the cell surface—either α and β chains in the case of circulating T cells or γ and δ chains in T cells localized to the gut (there are yet more chains expressed during development).
There are at least two goals to immune repertoire profiling. First, a determination the unique sequences of TCRs. The CDR3 regions are the protein segments that give each T cell its unique recognition specificity. The CDR3 coding sequence is created when V regions join with J regions. Occasionally, a small D region may exist between the V and J regions. The join between V and J is error prone by design, such that when these segments are fused, there is an intentional process where random DNA bases are inserted. This process further elaborates the TCR diversity. In some embodiments, the methods provided herein provide a determination of the DNA sequences of the V-J region across many different T cells.
Second, a count of T cell clones is determined. During an infection, certain T cell clones (as defined by their TCRs) are expanded because they are effective against an invader. Counting the numbers of each clone, even if they have the same TCR, provides a profile of the TCRs.
When genomic DNA is isolated from a sample, such as from a whole blood sample that contains T cells, for example, a molecular DNA tag is added to each genomic fragment before amplification of the genomic DNA. In this way, each TCR gene has a unique tag. Even if the TCR sequence is the same, the tag allows distinguishing of clones from different T cells versus those that are replicates from the same cell.
Normally all of the V segments and J segments are separated from one another by large, intervening genomic sequences. Only in adaptive immune response genes, such as TCR genes or antibody encoding genes, are the V and J sequences brought together in close proximity. By selecting for short genomic fragments that have both a V region and a J region on the same fragment, it is possible to enrich for functional TCR genes. A short genomic fraction can include a fraction of less than about 400 base pairs, such as less than 400, less than 350, less than 300, less than 250, less than 200, less than 150, less than 100, less than 90, less than 80, less than 70, less than 60, less than 50, or less than 40 base pairs or within a range defined by any two of the aforementioned values. Enrichment of a functional TCR gene is achieved by a sequential hybridization strategy in which all J regions are retrieved with J region specific probes. A majority of the sequences may be unrearranged, germline J segments. Following amplification of this J region enriched clone pool, fragments that also contain V regions are retrieved from the initial J pool using V region specific probes.
Some embodiments provided herein relate to methods and systems for target hybrid capture-based TCR enrichment.
1. Extraction of genomic DNA from a sample. The sample is obtained from a tumor tissue, a region proximal to a tumor tissue, an organ tissue, peripheral tissue, lymph, urine, cerebral spinal fluid, a buffy coat isolate, whole blood, peripheral blood, bone marrow, amniotic fluid, breast milk, plasma, serum, aqueous humor, vitreous humor, cochlear fluid, saliva, stool, sweat, vaginal secretions, semen, bile, tears, mucus, sputum, or vomit, or any other specimen thought to contain T cells. Genomic DNA is extracted by methods known in the art, including, for example, salting-out methods, organic extraction methods, cesium chloride density gradient methods, anion-exchange methods, and silica-based methods (Green, M. R. and Sambrook J., 2012, Molecular Cloning (4th ed.), Cold Spring Harbor, N.Y.: Cold Spring Harbor Laboratory Press).
2. Fragmentation of genomic DNA to an average size of about 300 bp or 300 bp followed by end repair. Because V and J regions are normally separated by large distances (>1000 bp) in unrearranged genomes and they only move into close proximity (<100 bp) in rearranged TCR genes, this fragmentation and the subsequent demand that a fragment have both a J region and a V region heavily enriches for TCR-encoding genes. Fragmentation can be performed by standard fragmentation techniques, including, for example, shearing, sonication, or enzymatic digestion; including restriction digests, as well as other methods or combinations of these approaches. In particular embodiments, any method known in the art for fragmenting DNA can be employed with the present disclosure.
3. As shown in
J region probes may be 89 nt in length. They may include a 45 nt tail that is complementary to biotinylated oligo 588 (e.g., SEQ ID NO: 232). This may be followed currently by a 4 nt random sequence (NNNN). More specific and longer sequences may be used. The 40 nt J region probes may be a combination of the J coding region that comes after the conserved triplet codon for F (inclusive of the F triplet). However, the J coding region is short, so these probes also include the genomic sequences found just 3′ of the J coding regions.
The J probes may have a tail sequence that is annealed to a complementary, biotinylated sequence (e.g., 588 J-probe complement, GGTAGTGTAGACTTAAGCGGCTATAGGGACTGGTCATCGTCATCG/3BioTEG/, SEQ ID NO: 232, Table 3). The biotin moiety is used for purification by attachment of the probe:genomic DNA complex to streptavidin-coated magnetic beads.
TCR J probes (
4. As shown in
5. An amplification segment is ligated to J region clones and subsequently PCR amplified (
6. As shown in
V region probes may be 101 nt long (
TCR V probes may include a 45 nucleotide tail sequence, followed by a tag of random nucleotides (e.g., NNNN), wherein N is A, T, C, or G, and wherein the tag can be 2-10 nucleotides in length, such as 2, 3, 4, 5, 6, 7, 8, 9, or 10 nucleotides in length, following by a J region probe sequences, as shown in the table below.
7. The annealed V region probes are extended. This copy of a copy is what actually is sequenced, following amplification with V probe and J probe specific primers. The temporary adaptor is lost.
8. As shown in
The clones are first amplified with primers that both add the sequences required for Illumina sequencing and that “index” each sample so that samples may be analyzed together. Indexing is achieved by amplifying each sample with a unique primer pair. Once the clones are amplified, they are sequenced in three separate steps using the specific sequencing primers. One PCR primer (CAC3 FLFP, oligo 568 AATGATACGGCGACCACCGAGATCTACACGTGACTGGCACGGGAGTTGATCCTG GTTTTCAC, SEQ ID NO: 229) is common to all samples. The other primer (chosen from oligos 607-638, SEQ ID NOs: 236-267) is unique to a sample and it marks each independent sample with its own “index.” In
9. Informatics analysis is then performed on the sequenced clones. Embedded in the sequencing data is the T cell repertoire. “Repertoire” in this case means a quantitative listing of all observed V-CDR3-J sequences. The ID tags were added in order to enable a count of different T cells with the same TCRs as two different events. This is important when assessing an immune response, for example, a T cell response directed against a tumor that is stimulated by immunotherapy.
The overall T cell repertoire data from a single sample is large. For example, in one microgram of whole blood DNA, about 5000 different TCR alpha chain and 5000 different TCR beta chain sequences may be present. One microgram of human genomic DNA has about 167,000 diploid genomes and about 5% of the genomes present are from T cells, it is reasonable to expect to count about 8000 unique T cells (unique α+β TCRs) per analyzed sample. Many times, the exact sequence is observed multiple times, and one function of post-sequence analysis is to condense these into a unique, consensus TCR.
In some embodiments, a data analysis, including a heatmap of TCRs, may be recognizable within a person's samples that are collected at intervals of weeks. Thus, in some embodiments, the T cell repertoires are reasonably stable over time. They can shift dramatically in response to an infection, a sickness, or in response to immune checkpoint blocker therapy in a cancer patient. In addition, in some embodiments, the heatmaps between different individuals are different from one another.
The primary objective of TCR analysis is counting. Each legitimate sequence is derived from a unique T cell, and the end result is census of all the T cells present in one microgram of whole blood genomic DNA.
Because each α chain is derived from the pairwise combination of 45 possible V regions and 54 possible J regions—representing a total of 2430 possible combinations—classifying the population based on the number of independent α chain clones of a particular V region that is joined to a specific J region in a table format provides a practical overview of the T cell population. Similarly, there are 45 possible β chain V regions and 12 possible β chain J regions—a total of 540 possibilities—that are also amenable to graphical display if provided in table format.
At least four elements may be taken into consideration for counting purposes. These include: 1) the J probe UMI—the first four bases of READ2; 2) the J probe sequence—the last 20 bases of READ2 (in some instances this 20 base sequence is not unique and therefore two or three α chain sequences are condensed together); 3) the V probe sequence—bases 5 through 14 of READ1 (this is the identifier that uniquely tags each V region probe; and 4) the CDR3 sequence (for example, bases 60-69 of READ1)
In addition, there are at least two kinds of artifacts in the data. The artifacts may include: 1) clones generated by probe-probe interactions, reads derived from these clones may be short and have terminal vector sequence (e.g. GCCGTCTTCTGCTTG; SEQ ID NO: 268) or they may possess J probe ACC4 primer sequences (e.g. GGTAGTGTAGACTTA; SEQ ID NO: 269). These artifacts add clones that should not be counted; and 2) clones lost because of single base read errors. The classification system described herein may include 30 error-free bases (20 for J and 10 for V) for a clone to be counted. Analyses that tolerate mismatches may increase the number of clones that are currently removed from counting consideration.
An additional artifact may occur with abundant unoccupied probes. The 3′ to 5′ exonuclease activity of T4 DNA polymerase is capable of generating a blunt end on these molecules, which then becomes a substrate for ligation to the P1 adaptor sequence (
Additional alternatives are disclosed in further detail in the following examples, which are not in any way intended to limit the scope of the claims.
Genomic DNA samples collected from various sources were purified using the Oragene saliva collection kit. The oligonucleotides that enable post-processing suppressive PCR, full-length amplification and sequencing are shown in
The probes are shown in
Four gDNAs (F, S, C and L) were diluted to 20 ng/μL in 150 μL final volume. The samples were sonicated to 500 bp and 125 μL was purified with 125 μL of beads. The starting material and purified, fragmented gDNA for each sample was run on a gel shown in
For capture, 10 μL of gDNA sample was heated to 98° C. for 2 minutes (to achieve strand dissociation) and cooled on ice. 5 μL of 4× bind and 5 μL of the 49 probe tagged V2 probe pool (probes listed in
For processing and adaptor ligation, 100 μL of T4 mix was made that contained: 60 μL water, 10 μL NEB “CutSmart” buffer, 15 μL 50% PEG8000, 10 μL 10 mM ATP, 1 μL 1 mM dNTP blend, 1 μL T4 gene 32 protein (NEB), and 0.5 μL T4 DNA polymerase (NEB). 25 μL of this mix was added to each of the four samples and incubated at 20° C. for 15 minutes followed by a 70° C. incubation for 10 minutes to heat inactivate the T4 polymerase. Following this 1.25 μL of adaptor (10 μM in ligation strand, pre-annealed) and 1.25 μL of HC T4 DNA ligase were added. This mixture was further incubated at 22° C. for 30 minutes and 65° C. for 10 minutes.
Here, one attractive feature of library free is that processed complexes are, at least in theory, still attached to beads. Beads were pulled from the ligation buffer and washed once with 200 μL of TEzero. The complexes were then resuspended in 2 μL. For amplification, the idea is to use single primer amplification in a 20 μL volume to both amplify target fragments and to enrich for long genomic fragments over probe “stubs”. Following this, a larger volume PCR reaction with full length primers will be used to create a “sequence-ready” library.
A Q5-based, single primer PCR amplification buffer was made by combining 57 μL water, 20 μL 5×Q5 reaction buffer, 10 μL of single primer 117 (see list above), 2 μL of 10 mM dNTPs, and 1 μL of Q5 hot start polymerase. Eighteen μL was added to each tube followed by amplification for 20 cycles (98° C.-30 seconds; 98° C.-10 seconds, 69° C.-10 seconds, 72° C.-10 seconds for 20 cycles; 10° C. hold). Following this, the beads were pulled out and the 20 μL of pre-amp supernatant was transferred to 280 μL of PCR mix that contained 163.5 μL water, 60 μL 5×Q5 buffer, 15 μL of forward primer 118 (10 μM), 15 μM of reverse primer 119 (10 μM), 6 μL of 10 mM dNTPs, 13.5 μL of EvaGreen+ROX dye blend (1.25 parts EG to 1 part ROX), and 3 μL of Q5 hot start polymerase (adding the dye to all reactions was unintended). Two of 100 μL aliquots were amplified by conventional PCR (98° C.-10 seconds, 69° C.-10 seconds, 72° C.-10 seconds) and quadruplicate ten μL aliquots were amplified under qPCR conditions. The amplification plot shown in
qPCR capture assays were used to determine whether gene specific targets were captured and selectively amplified. The target regions for various assays are shown in Table 2.
For qPCR analysis, genomic DNA from sample F at 10 ng/μL (2 μL is added to 8 μL of PCR mix to give a final volume and concentration of 10 μL and 2 ng/μL, respectively) was used as control. Purified processed material from the F and S samples was diluted to 0.01 ng/μL=10 pg/μL and 2 μL was added to each 8 μL PCR reaction to give a final concentration of 2 pg/μL. These are more or less standard qPCR assay conditions to evaluate any capture reaction. The results are shown in
To this point, library-free was a collection of promising-looking smears. The qPCR data indicates that the technology is in fact very effective at retrieving the targeted genomic regions and at leaving off-target regions behind (Assays 6, 8). The fold purifications, often >500,000-fold, are directly comparable to our SOP technology.
The results from the preliminary investigation described in Example 1 were sufficiently compelling for investigation of the enzymatic requirements for complex processing. The design of experiment is shown in Table 3.
To make capture complexes for analysis, twelve identical reactions were created. Ten μL of 135 ng/μL sonicated gDNA was melted, annealed with tagged V2 probe, stuck to strep coated beads, washed and resuspended in TEzero as described above. Five hundred μL of processing master mix was prepared by combining 270 μL water, 50 μL 10× CutSmart buffer, 50 μL of 10 mM ATP, 75 μL of 50% PEG8000, and 5 μL of 10 mM dNTPs. This buffer was divided into 10 of 90 μL aliquots (duplicate tests were performed) and enzyme was added in the amounts described above (per 90 μL of master mix was added 1 μL of T4 gene 32 protein, 0.5 μL of T4 polymerase, 5 μL of adaptor and/or 5 μL of HC T4 ligase). Following T4 fill-in and ligation as described above, the complexes were washed free of processing mix in TEzero and resuspended in 2 μL TEzero. Complexes were resuspended in 20 μL final volume each of single primer amplification mix and amplified for 20 cycles as described above. The beads were then pulled aside using a magnet and the 20 μL clarified amplification was diluted into 180 μL of full-length F+R (118+119) PCR amplification mix. Fifty μL was pulled aside for qPCR analysis and the remaining 150 μL was split in two and amplified by conventional PCR. The 50 μL qPCR samples were mixed with 2.5 μL of dye blend and 10 μL aliquots were monitored by fluorescence change. The traces of this experiment are shown in
The methods described in Examples 1 and 2 were used to produce a DNA sequencing library with the four Coriell samples. Each one of the four samples was coded with an individual index code in the final PCR step. The creation of such libraries highlights that library-free methods demand that all samples in a collection be processed separately, which is undesirable. The final library constituents (shown separately prior to pooling) are shown in the gel image in
The library-free methods worked well for CNV analyses. Unique read counts for the X-linked gene PLP1 were normalized to the autosomal loci KRAS and MYC and the plot of these data is shown in
The following example demonstrates the profiling of one microgram of genomic DNA. This genomic DNA can be isolated from whole blood cells, from the buffy coat, from peripheral blood mononuclear cells, or from other samples and tissues as described herein. In reality, all of these are similar sources of nucleated leukocytes that include T cells that have alpha and beta chain TCRs. The steps described in this protocol are illustrated in
The adaptor for this Example was made from oligos 596 (J-probe-part, CCGCTTAAGTCTACACTAC/3ddC/, SEQ ID NO: 233) and 597 (J-probe-lig, /5Phos/GGTAGTGTAGACTTAAGCGGCTATAGG, SEQ ID NO: 234). 20 μL of each oligo was combined in 160 μL of TEzero+25 mM NaCl to generate a duplex with a final concentration of 10 μM.
The PCR primer for this experiment was oligo 489 (ACC4_27, CCTATAGCCGCTTAAGTCTACACTACC, SEQ ID NO: 228). 50 μL of oligo 489 was combined with 450 μL of TEzero to obtain 10 μM PCR primer.
The following oligonucleotides were also used, as described below: 568 PCR Primer post V-hyb (SEQ ID NO 229); 571 Forward Sequencing Primer (SEQ ID NO: 230); 573 Reverse Sequencing Primer (SEQ ID NO: 231); and 606 Index Sequencing Primer (SEQ ID NO: 235).
In separate reactions, 130 μL of gDNA was sonicated from patient samples VSC7-2, 7-3, 7-4 and 7-5 to 300 bp. 125 μL of sonicated gDNA was added to 150 μL of beads. The mixture was washed twice with 70% EtOH. The pellets were resuspended in 50 μL TEZ. 1000 ng of sonicated gDNA was added to a new tube. Standard end repair was performed (ST1, ST2). Each end repaired sample was captured with: 12.5 μL of 1.0 nM TRAJ Probe+12.5 μL of 1.0 nM TRBJ Probe. The mixture was heated to 98° C. for 2 minutes, and 112.5 μL of hybridization buffer was added. Run on O/N at 65° C. hybridization.
Following hybridization, the mixture was washed as followed. 150 μL of the hybridization reactions was mixed with 40 μL of washed MyOne streptavidin beads in 1 mL TT. The mixture was incubated for 30 minutes with occasional mixing. Beads were pulled out and resuspended in 400 μL TT. Two 200 μL aliquots were separated in PCR strip tubes. The beads were pulled down and resuspended in 200 μL per tube wash buffer, incubated at 45° C. for 5 minutes, pulled out and resuspended in 200 μL TEzero, followed by pulled out and resuspended in 20 μL per tube TEzero.
For T4 extension, 80 μL of T4 mix containing 52.5 μL water, 10 μL 10× CutSmart buffer, 15 μL 50% PEG8000, 1 μL of 10 mM dNTPs, 1 μL T4 Gene 32 protein, and 0.5 μL T4 DNA polymerase was prepared. The mixture was incubated at 20° C. for 15 minutes followed by 70° C. for 10 minutes. The beads were pulled out and resuspended in 200 μL TEzero, pulled out and resuspended in 50 μL TEzero. 20 μL of adaptor was added and 30 μL of standard ligation cocktail (10 μL 10× ligation buffer, 15 μL 50% PEG8000, 5 μL T4 DNA ligation buffer) was added. The standard ligation protocol was run (60 minutes at 20° C., followed by 10 minutes at 65° C.).
The beads were pulled out and resuspended in 20 μL TEzero. 80 μL of “C+P” PCR mix: 50 μL 2× master blend, 10 μL TCR PCR primer 489 (SEQ ID NO: 228), and 20 μL water was added. The sequence was amplified for 5 cycles.
The beads were pulled out, and 60 μL of supernatant was added to 240 μL post C+P PCR mix: 120 μL 2× master blend, 24 μL TCR primer 489 (SEQ ID NO: 228), and 96 μL water. The amplification was monitored by qPCR.
All samples were amplified for 10 cycles (regardless of qPCR results). The beads were purified, and resuspended in 20 μL H2O for a total of 40 μL H20. Each 40 μL sample was captured by adding: 10 μL of 1.0 nM TRAV Probe+10 μL of 1.0 nM TRBV Probe. The mixture was heated to 98° C. for 2 minutes. 90 μL of hybridization buffer was added, and run on O/N 65° C. hybridization.
The mixture was washed post hybridization by combining 150 μL hybridization reactions with 40 μL of washed MyOne streptavidin beads in 1 mL TT. The mixture was incubated for 30 minutes with occasional mixing. The beads were pulled out and resuspended in 400 μL TT. Two 200 μL aliquots were split in PCR strip tubes. The beads were pulled out, resuspended in 200 μL per tube wash buffer, and incubated at 45° C. for 5 minutes. The beads were pulled out and resuspended in 200 μL TEzero, and then pulled out and resuspended in 20 μL per tube TEzero.
80 μL of “C+P” PCR mix was added: 50 μL 2× master blend, 10 μL TCR PCR primer 568 (SEQ ID NO: 229), 10 μL TCR PCR index primer, and 20 μL water. The mixture was amplified for 5 cycles, the beads pulled out, and 60 μL of supernatant was added to 240 μL post C+P PCR mix: 120 μL 2× master blend, 12 μL TCR PCR primer 568 (SEQ ID NO: 229), 12 μL TCR PCR index primer (including index primers 607 (SEQ ID NO: 236), 608 (SEQ ID NO: 237), 623 (SEQ ID NO: 252), and 624 (SEQ ID NO: 253) for patient samples 7-2, 7-3, 7-4 and 7-5, respectively), and 96 μL water. Amplification was monitored by qPCR. Beads were purified by resuspending in 20 μL TEZ for a total of 40 μL TEZ.
Follow standard MiSeq protocol. Use the following primers in the corresponding MiSeq wells. Primer 571 FTCSP (SEQ ID NO: 23) to 18 Primer; 606 ITCSP (SEQ ID NO: 235) to 19 Primer; and 573 RTCSP (SEQ ID NO: 231) to 20 Primer.
The raw output from the Illumina MiSeq run produced approximately 8 million sequencing reads, about 2 million reads per patient sample after parsing the data using the sample index information. The data for each patient was filtered in several steps that included: discarding reads that did not have a legitimate V region or J region probe sequence; discarding reads that did not have a protein coding open reading frame in the CDR3 region between the V and the J probes (Importantly, the observed distribution of CDR3 sequence lengths (average=36 bases for alpha chains and 39 bases for beta chains) was concordant with previous literature reports); identifying redundant reads into a single, consensus TCR “unique sequence”; classifying unique read sets into alpha or beta chains; classifying alpha unique reads or beta unique reads according to their V and J regions; counting the number of TCRs in each V/J intersection (pixel); and presenting the population distribution of TCRs in patient series 7-2 through 7-5 in heat maps.
Approximately 5000 unique alpha and 5000 unique beta TCR sequences were observed in each sample (the range was 3217 to 7684 unique sequences). An example of a heat map for one alpha chain sample is shown in
One microgram of human genomic DNA is the equivalent of about 150,000 diploid genomes, or, in other words, representative of 150,000 cells. In whole blood, roughly 4-7% of nucleated cells are T cells. Therefore, the expectation is that 6000 to 10,500 unique TCRs in each sample should be observed. The observed density of about 5000 unique TCRs is consistent with this expectation, especially when the fact that cancer patients are often immunosuppressed by therapy is taken into account. The TCR repertoire produced by the methods provided herein is likely to reflect a snapshot of the peripheral, circulating T cells present in a sample. Modifying J probe tags will expand the detection of redundant clones and on profiling of the tumor infiltrating T cells in resected tumor tissue.
Development of the method requires several iterations that were not initially obvious from a priori consideration of the assay. The method has significant clinical utility in applications such as infectious disease monitoring and assessment of the efficacy of immune-oncology therapies.
It is to be understood that the description, specific examples and data, while indicating exemplary embodiments, are given by way of illustration and are not intended to limit the various embodiments of the present disclosure. Various changes and modifications within the present disclosure will become apparent to the skilled artisan from the description and data contained herein, and thus are considered part of the various embodiments of this disclosure.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2020/038474 | 6/18/2020 | WO |
Number | Date | Country | |
---|---|---|---|
62887938 | Aug 2019 | US |