DNA-Barcoded pMHC Tetramers for Detection of Single Antigen-Specific T Cells by Digital PCR
Antigen-specific T cells are found at low frequencies in circulation but carry important diagnostic information as liquid biomarkers in numerous biomedical settings, such as monitoring the efficacy of vaccines and cancer immunotherapies. To enable detection of antigen-specific T cells with high sensitivity, we develop peptide-MHC (pMHC) tetramers labeled with DNA barcodes to detect single T cells by droplet digital PCR (ddPCR). We show that site-specific conjugation of DNA via photocleavable linkers allows barcoded tetramers to stain T cells with similar avidity compared to conventional fluorescent tetramers, and efficient recovery of barcodes by light with no loss in cell viability. We design an orthogonal panel of DNA-barcoded tetramers to simultaneously detect multiple antigen-specific T cell populations, including from a mouse model of viral infection, and discriminate single cancer-specific T cells with high diagnostic sensitivity and specificity. This approach of DNA-barcoding can be broadened to encompass additional rare cells for monitoring immunological health at the single cell level.
<!>DNA Barcode Conjugation to StvC<!>Barcode Release and Cell Viability<!>DNA Barcode Analysis by ddPCR<!>Results and Discussion
<p>Rare cells circulating at low frequencies in blood (e.g., 1–10 cells in 1 ml blood)1, such as stem cells2,3, circulating tumors cells4,5, and antigen-specific immune cells6,7, have important applications in basic research and translation8,9, including use as liquid biomarkers for monitoring disease. Low frequency T cell subsets are particularly important and carry a wealth of information; for example, the number of naïve T cells predicts the strength of an immune response10,11, the dynamics of T cell expansion and contraction correlate with pathogen infection and clearance12–16, and the quantity of circulating antigen-specific T cells is a key biomarker of efficacy for cancer immunotherapies including engineered T cells and checkpoint blockade17–21. Despite their importance, it remains technically challenging to detect antigen-specific T cells, requiring expression of engineered peptide MHC (pMHC) multimers22, and complex analytical platforms capable of multiplexed analysis such as multi-color combinatorial flow cytometry23,24, next generation sequencing25, or cytometry by time of flight (CyTOF)26. Although powerful, these approaches require deep technical expertise and target cells cannot be recovered after analysis, precluding routine use for biomedical diagnostics.</p><p>Here we develop pMHC tetramers labeled with photo-labile DNA barcodes to allow detection of rare viral-specific and cancer-specific T cells by droplet digital PCR (ddPCR). Analytical platforms based on DNA leverage the unique properties of nucleic acids, such as base-pairing, to achieve ultrasensitive or highly multiplexed detection. These approaches include the use of protein-DNA conjugates for self-assembled cellular arrays27,28 and engineered tissues29, and dynamic DNA gates for multiplexed cell sorting30. DNA barcodes have also been incorporated into lipid nanoparticle libraries for drug delivery screens in vivo31, and used for ultrasensitive detection of clinical biomarkers with nanoparticles32,33, immunoPCR34–36, and proximity ligation37,38. By labeling pMHC tetramers with photo-labile DNA barcodes, our approach allows staining of antigen-specific T cells in a complex sample, harvest of barcodes by light while preserving cell viability, and ultrasensitive detection of single T cells by droplet digital PCR.</p><!><p>DNA barcodes containing a 5' NH2, an internal Cy5, and an internal photocleavable linker (see Table S1) were purchased from Integrated DNA Technologies with HPLC purification and reconstituted to 100μM in 1X PBS. StvC was produced as described in Supporting Information. For bioconjugation, StvC was reduced with 10mM TCEP and then reacted with 50-fold excess MHPH (Solulink), and amine-terminated DNA was reacted with 25-fold excess S-4FB (Solulink) for 4 hours. Excess linker was removed by buffer exchanging into citrate buffer (50mM sodium citrate, 150mM NaCl, pH 6) using Amicon spin filters (Millipore). Functionalized DNA and StvC were combined at a 1:1 ratio and reacted overnight before purifying on a Superdex 200 Increase 10/300 GL column using an AKTA Pure FPLC (GE Healthcare). Conjugation was verified by SDS-PAGE followed by Coomassie staining.</p><p>To compare binding avidity of DNA-barcoded to conventional fluorophore-labeled pMHC tetramers, tetramers were synthesized by reacting biotinylated Db-GP10025–33 pMHC with StvC conjugated with either Cy5 or DNA barcode-Cy5 at a 4:1 molar ratio. Splenocytes isolated from pmel TCR transgenic mice (Jackson Labs) were co-stained with anti-mouse CD8 (Clone: 53–6.7, BD Biosciences) and tetramer at various concentrations, and Cy5 intensity was analyzed by flow cytometry using a BD Accuri C6. Mean Cy5 intensity in the tetramer positive population was used as a measure of tetramer binding, and binding data was fit to a one site model using GraphPad Prism 6.0.</p><!><p>pMHC tetramers were generated as described above using PE-conjugated streptavidin (Invitrogen) or StvC conjugated with DNA barcode containing an internal photocleavable linker and a 3' FAM (see Table S1). Splenocytes were co-stained with anti-mouse CD8 (Clone: 53–6.7, BD Biosciences), tetramer-PE, and tetramer-DNA-FAM, and then tetramer positive cells were enriched using anti-PE microbeads (Miltenyi Biotec) and LS columns (Miltenyi Biotec) according to the manufacturer's instructions. Eluted cells were resuspended in 100μl PBS + 1% BSA in 1.5ml Eppendorf tubes, placed on ice, and irradiated for various durations in a UVP CL-1000L Crosslinker (365nm). FAM and PE signal was measured by flow cytometry using a BD Accuri C6. To measure cell viability after UV irradiation, unstained splenocytes were irradiated under UV light for various durations, stained using a Annexin V Apoptosis Detection Kit (BD Biosciences) according to the manufacturer's instructions, and then analyzed using a BD Accuri C6.</p><!><p>To check for orthogonality between multiple barcodes and primer sets, free barcodes were serially diluted 1:10 in nuclease-free water to 10−13 M followed by a single 1:5 dilution to 2×10−14 M. A QX200 EvaGreen Supermix master mix (Bio-Rad) was created containing 125nM primers, and 5μl of each diluted barcode was combined with 20μl of each primer master mix. 20μl of this reaction was used to create droplets using a QX200 Droplet Generator (Bio-Rad) according to the manufacturer's instructions. Droplets were thermal cycled in a Eppendorf Mastercycler using the following conditions: 95°C for 5 min, 40 cycles of 95°C for 30 sec and 60°C for 1 min, 4°C for 5 min, 90°C for 5 min, and hold at 4°C. Droplets were then analyzed on a QX200 Droplet Reader (Bio-Rad).</p><p>For cell studies, splenocytes from OT1, P14, and/or pmel TCR transgenic mice (Jackson Labs) or from mice infected with LCMV Armstrong were stained with DNA barcoded tetramers. Barcodes were isolated by irradiating stained cells in a UVP CL-1000L Crosslinker (365nm) on ice for 30 min. Cells were spun down at 1000xg for 5 min and the supernatant containing barcode was analyzed by ddPCR as described above. See Supporting Information for more detailed experimental conditions.</p><!><p>Conventional pMHC tetramers are produced by combining biotinylated pMHC monomers with streptavidin to produce a four-fold valent construct. To label tetramers with DNA barcodes, we expressed recombinant streptavidin engineered with C-terminal cysteines (StvC)39 to introduce free thiol groups for site-specific conjugation. We reacted DNA barcodes containing a 5' amine and a 3' Cy5 fluorescent reporter (Table S1) with StvC using heterobifunctional linkers (MHPH (maleimide-to-hydrazine) and S-4FB (NHS ester-to-aldehyde)) to form a stable bis-arylhydrazone bond as described previously28 (Figure 1A). The resulting StvC-DNA conjugates were purified by size-exclusion chromatography (Figure 1B) and the coupling efficiency was quantified by gel mobility shift assays (~1.4 DNA barcodes per StvC molecule; Figure 1C). To analyze tetramer binding avidity, we compared the ability of DNA-barcoded tetramers to stain cancer-specific T cells isolated from transgenic pmel mice, which recognize the melanoma antigen gp100 (Db-GP10025–33), to control fluorescent tetramers made by reacting StvC with Cy5 maleimide (Figure 1D). By quantifying the mean fluorescent intensities at tetramer staining concentrations ranging from 0–600 nM, we observed nearly identical binding curves and calculated apparent dissociation constants (Kd,app)40 of 32 nM and 12 nM for DNA-barcoded and fluorescent tetramers respectively (Figure 1E). Together these results demonstrated that DNA-barcoded tetramers, made by site-specific coupling to StvC, stain antigen-specific T cells with binding similar to standard fluorescent tetramers.</p><p>We next sought to incorporate photo-cleavable linkers to allow harvest of DNA barcodes from the surface of target T cells without affecting cell viability. We synthesized DNA barcodes containing a photocleavable 1-(2-nitrophenyl) ethyl linker, which is sensitive to cleavage by long-wavelength UV light (~365nm). To quantify cleavage efficiency on the surface of T cells, we co-stained pmel splenocytes with photocleavable tetramers labeled with FAM, and conventional tetramers labeled with PE as a control (Figure 2A). We found that co-stained T cells initially exhibited fluorescence in both FAM and PE channels as expected (0 minute sample), but after 30 minutes of UV-light irradiation, we detected barcode cleavage by loss of FAM fluorescence (Figure 2B). By contrast, the intensity of PE fluorescence remained at identical levels after photo-cleavage, indicating that the ability of tetramers to bind to T cells was not affected by exposure to UV-light. We further stained splenocytes with Annexin V and propidium iodide to quantify the percentage of apoptotic and necrotic cells, respectively, following photo-cleavage. We observed no increase in the frequency of Annexin V+/PI− or Annexin V+/PI+ cells post-irradiation, indicating that our barcode isolation workflow did not affect cell viability (Figure 2C). Taken together, these results showed that DNA barcodes can be efficiently recovered from the surface of target cells by UV light while preserving cell viability.</p><p>We next designed a multiplex set of DNA barcodes (A–C) containing orthogonal primer binding regions to minimize potential amplification biases during barcode quantification by ddPCR (Table S1); barcodes contained 50% G/C content with primer melting temperatures between 63–64°C and were screened for minimal secondary structure in silico using NUPACK41. ddPCR uses a water-in-oil emulsion to dilute target DNA until a single molecule is isolated within a nanoliter droplet. Only droplets containing a strand of target DNA are amplified during thermal cycling, and the final frequency of amplified droplets allows the target DNA concentration to be calculated by Poisson statistics. In samples containing a matched barcode-primer set, target counts were amplified to greater than 103 copies per μl compared to unmatched controls (<5 copies per μl), and the coefficient of variance of the barcode-primer sets was less than 3% (Figure 3A). These results were within the range of values reported by previous studies designing orthogonal primers for unbiased multiplexed PCR42. We synthesized a panel of three pMHC tetramers to simultaneously detect T cells specific for the model immunogen ovalbumin (Kb-OVA257–264; Barcode A), the viral antigen LCMV glycoprotein (Db-LCMV GP33–41; Barcode B), and the human melanoma antigen gp100 (Db-GP10025–33; Barcode C). In samples containing a mixture of different frequencies of transgenic T cells specific for these antigens (OT1, P14, and pmel T cells respectively), we stained T cells en masse, and quantified barcodes released by UV-light cleavage by ddPCR (Figure 3B,C). For each antigen-specific T cell population, barcode counts were highly correlated (R2 = 0.95–0.99) with T cell frequency within the population as measured by flow cytometry (Figure 3D, Figure S1). Importantly, we found that detection of each strain was barcode sequence independent as switching the barcode associated with each tetramer resulted in similar counts by ddPCR (Figure S2). Because TCR transgenic mice express high affinity, monoclonal TCRs, we then tested if this approach could be used to monitor T cell dynamics during an endogenous polyclonal immune response, where T cells may have a broad range of avidities to different antigens derived from a single pathogen43. We infected mice with the Armstrong strain of lymphocytic choriomeningitis virus (LCMV), isolated splenocytes 0, 3, 8, and 148 days post-infection, and performed barcode analysis to detect T cells specific for the LCMV antigens GP33, GP276, and NP396 (Figure 4A). We found that barcode counts correlated well with frequencies as measured by flow cytometry (Figure 4B), which validated this approach for multiplexed analysis of T cells.</p><p>To determine whether barcoded tetramers could resolve single T cells, we first quantified the sensitivity of ddPCR by assaying samples containing free DNA barcodes at a copy number range of 7.5×101–1.5×105 (or 5×10−14–2.5×10−17 M in 5 μl), which we selected to correspond to the average number of surface TCRs per cell44. Under these conditions, barcode counts were linearly quantified across the entire range of concentrations (R2 = 0.99) down to 0.025 fM (Figure S3), demonstrating the potential of ddPCR to detect free barcodes at levels expected from single cells. We verified this observation using a ladder of P14 T cells sorted by fluorescence-activated cell sorting (FACS) and stained with barcoded tetramer and similarly found a linear relationship between cell number and barcode concentration (R2=0.99), even for low cell counts (Figure 5A). Sensitivity was not reduced in the presence of high background, as we could detect P14 T cells spiked into control B6 splenocytes at frequencies as low as 5 P14 T cells in 1 million B6 splenocytes (~0.01% when gated on CD8 positive cells) (Figure S4). To determine the diagnostic sensitivity and specificity of single T cell detection, we quantified barcode counts from pmel splenocytes co-stained with barcoded GP100 and control GP33 tetramers, and found that single cell barcode counts for the pmel tetramer were significantly elevated compared to control P14 tetramers (Figure 5B). By receiver operating characteristic (ROC) analysis – which is used to determine the sensitivity (true positive rate) and specificity (true negative rate) of a diagnostic classifier – DNA-barcoded tetramers discriminated single antigen-specific T cells with an area under the curve (AUC) accuracy of 0.91, which is significantly higher than a random classifier (AUC = 0.5, perfect classifier AUC = 1.0) (Figure 5C). These results show that staining by DNA-barcoded tetramers followed by ddPCR analysis can be used to measure the frequency of multiple antigen-specific T cell populations from biological samples with single-cell sensitivity.</p><p>Here we develop DNA barcoded pMHC tetramers for analysis of antigen-specific T cells by ddPCR. We show that barcoded tetramers stain T cells equivalent to their fluorescent counterparts and demonstrate that barcodes are non-destructively isolated from the surface of cells by UV irradiation. Using ddPCR, we show multiplexed detection of T cells against three different antigens from an LCMV-infected mouse and achieve single cell sensitivity. Importantly, this barcoding approach is compatible with advances in tetramer staining protocols45 and the development of higher order multimers46 for detection of T cells expressing low affinity TCRs. In the future, the library of DNA barcodes can be expanded to screen peptides for TCR reactivity47 and to monitor the dynamics of T cells through naïve, effector, and memory states during the course of disease.</p>