Theres Oakes

Computational analysis of stochastic heterogeneity in PCR amplification efficiency revealed by single molecule barcoding

The polymerase chain reaction (PCR) is one of the most widely used techniques in molecular biology. In combination with High Throughput Sequencing (HTS), PCR is widely used to quantify transcript abundance for RNA-seq, and in the context of analysis of T and B cell receptor repertoires. In this study, we combine DNA barcoding with HTS to quantify PCR output from individual target molecules. We develop computational tools that simulate both the PCR branching process itself, and the subsequent subsampling which typically occurs during HTS sequencing. We explore the influence of different types of heterogeneity on sequencing output, and compare them to experimental results where the efficiency of amplification is measured by barcodes uniquely identifying each molecule of starting template. Our results demonstrate that the PCR process introduces substantial amplification heterogeneity, independent of primer sequence and bulk experimental conditions. This heterogeneity can be attributed both to inherited differences between different template DNA molecules, and the inherent stochasticity of the PCR process. The results demonstrate that PCR heterogeneity arises even when reaction and substrate conditions are kept as constant as possible, and therefore single molecule barcoding is essential in order to derive reproducible quantitative results from any protocol combining PCR with HTS.

Dynamic Perturbations of the T-Cell Receptor Repertoire in Chronic HIV Infection and following Antiretroviral Therapy.

HIV infection profoundly affects many parameters of the immune system and ultimately leads to AIDS, yet which factors are most important for determining resistance, pathology, and response to antiretroviral treatment – and how best to monitor them – remain unclear. We develop a quantitative high-throughput sequencing pipeline to characterize the TCR repertoires of HIV-infected individuals before and after antiretroviral therapy, working from small, unfractionated samples of peripheral blood. This reveals the TCR repertoires of HIV+ individuals to be highly perturbed, with considerably reduced diversity as a small proportion of sequences are highly overrepresented. HIV also causes specific qualitative changes to the repertoire including an altered distribution of V gene usage, depletion of public TCR sequences, and disruption of TCR networks. Short-term antiretroviral therapy has little impact on most of the global damage to repertoire structure, but is accompanied by rapid changes in the abundance of many individual TCR sequences, decreases in abundance of the most common sequences, and decreases in the majority of HIV-associated CDR3 sequences. Thus, high-throughput repertoire sequencing of small blood samples that are easy to take, store, and process can shed light on various aspects of the T-cell immune compartment and stands to offer insights into patient stratification and immune reconstitution.

High-throughput sequencing of the T-cell receptor repertoire: pitfalls and opportunities

&NA; T‐cell specificity is determined by the T‐cell receptor, a heterodimeric protein coded for by an extremely diverse set of genes produced by imprecise somatic gene recombination. Massively parallel high‐throughput sequencing allows millions of different T‐cell receptor genes to be characterized from a single sample of blood or tissue. However, the extraordinary heterogeneity of the immune repertoire poses significant challenges for subsequent analysis of the data. We outline the major steps in processing of repertoire data, considering low‐level processing of raw sequence files and high‐level algorithms, which seek to extract biological or pathological information. The latest generation of bioinformatics tools allows millions of DNA sequences to be accurately and rapidly assigned to their respective variable V and J gene segments, and to reconstruct an almost error‐free representation of the non‐templated additions and deletions that occur. High‐level processing can measure the diversity of the repertoire in different samples, quantify V and J usage and identify private and public T‐cell receptors. Finally, we discuss the major challenge of linking T‐cell receptor sequence to function, and specifically to antigen recognition. Sophisticated machine learning algorithms are being developed that can combine the paradoxical degeneracy and cross‐reactivity of individual T‐cell receptors with the specificity of the overall T‐cell immune response. Computational analysis will provide the key to unlock the potential of the T‐cell receptor repertoire to give insight into the fundamental biology of the adaptive immune system and to provide powerful biomarkers of disease.

Quantitative Characterization of the T Cell Receptor Repertoire of Naïve and Memory Subsets Using an Integrated Experimental and Computational Pipeline Which Is Robust, Economical, and Versatile

The T cell receptor (TCR) repertoire can provide a personalized biomarker for infectious and non-infectious diseases. We describe a protocol for amplifying, sequencing, and analyzing TCRs which is robust, sensitive, and versatile. The key experimental step is ligation of a single-stranded oligonucleotide to the 3′ end of the TCR cDNA. This allows amplification of all possible rearrangements using a single set of primers per locus. It also introduces a unique molecular identifier to label each starting cDNA molecule. This molecular identifier is used to correct for sequence errors and for effects of differential PCR amplification efficiency, thus producing more accurate measures of the true TCR frequency within the sample. This integrated experimental and computational pipeline is applied to the analysis of human memory and naive subpopulations, and results in consistent measures of diversity and inequality. After error correction, the distribution of TCR sequence abundance in all subpopulations followed a power law over a wide range of values. The power law exponent differed between naïve and memory populations, but was consistent between individuals. The integrated experimental and analysis pipeline we describe is appropriate to studies of T cell responses in a broad range of physiological and pathological contexts.

The T Cell Response to the Contact Sensitizer Paraphenylenediamine Is Characterized by a Polyclonal Diverse Repertoire of Antigen-Specific Receptors

Paraphenylenediamine (PPD) is a common component of hair dyes and black henna tattoos and can cause skin sensitization and allergic contact dermatitis (ACD). The cutaneous inflammatory reaction associated with ACD is driven by both CD4+ and CD8+ T cells. However, the characteristics of such responses with respect to clonal breadth and magnitude are poorly defined. In this study, we have characterized the in vitro recall response of peripheral blood T cells prepared from PPD-allergic individuals to a PPD–human serum albumin (HSA) conjugate (PPD–HSA). Quantitative high throughput sequencing was used to characterize the changes in the repertoire of T cell receptor (TCR) α and β genes after exposure to antigen in vitro. The PPD conjugate induced expansion of T cells carrying selected TCRs, with around 800 sequences (around 1%) being 8 or more times as abundant after culture than before. The expanded sequences showed strong skewing of V and J usage, consistent with an antigen-driven clonal expansion. The complementarity-determining region 3 sequences of the expanded TCRs could be grouped into several families of related amino acid sequence, but the overall diversity of the expanded sample was not much less than that of a random sample of the same size. The results suggest a model in which PPD–HSA conjugate stimulates a broad diversity of TCRs, with a wide range of stimulation strengths, which manifest as different degrees of in vitro expansion.

Sequence and primer independent stochastic heterogeneity in PCR amplification efficiency revealed by single molecule barcoding

The polymerase chain reaction (PCR) is one of the most widely used techniques in molecular biology. In combination with High Throughput Sequencing (HTS), PCR is widely used to quantify transcript abundance for RNA-seq and especially in the context of analysis of T cell and B cell receptor repertoires. In this study, we combine molecular DNA barcoding with HTS to quantify PCR output from individual target molecules. Our results demonstrate that the PCR process exhibits very significant unexpected heterogeneity, which is independent of the sequence of the primers or target, and independent of bulk experimental conditions. The mechanistic origin of this heterogeneity is not clear, but simulations suggest that it must derive from inherited differences between different DNA molecules within the reaction. The results illustrate that single molecule barcoding is important in order to derive reproducible quantitative results from any protocol which combines PCR with HTS.

Urine-derived lymphocytes as a non-invasive measure of the bladder tumor immune microenvironment

Despite the advances in cancer immunotherapy, only a fraction of patients with bladder cancer exhibit responses to checkpoint blockade, highlighting a need to better understand drug resistance and identify rational immunotherapy combinations. However, accessibility to the tumor prior and during therapy is a major limitation in understanding the immune tumor microenvironment (TME). Herein, we identified urine-derived lymphocytes (UDLs) as a readily accessible source of T cells in 32 patients with muscle invasive bladder cancer (MIBC). We observed that effector CD8+ and CD4+ cells and regulatory T cells within the urine accurately map the immune checkpoint landscape and T cell receptor repertoire of the TME. Finally, an increased UDL count, specifically high expression of PD-1 (PD-1hi) on CD8+ at the time of cystectomy, was associated with a shorter recurrence-free survival. UDL analysis represents a dynamic liquid biopsy that is representative of the bladder immune TME that may be used to identify actionable immuno-oncology (IO) targets with potential prognostic value in MIBC.