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

Specificity, Privacy, and Degeneracy in the CD4 T Cell Receptor Repertoire Following Immunization.

T cells recognize antigen using a large and diverse set of antigen-specific receptors created by a complex process of imprecise somatic cell gene rearrangements. In response to antigen-/receptor-binding-specific T cells then divide to form memory and effector populations. We apply high-throughput sequencing to investigate the global changes in T cell receptor sequences following immunization with ovalbumin (OVA) and adjuvant, to understand how adaptive immunity achieves specificity. Each immunized mouse contained a predominantly private but related set of expanded CDR3β sequences. We used machine learning to identify common patterns which distinguished repertoires from mice immunized with adjuvant with and without OVA. The CDR3β sequences were deconstructed into sets of overlapping contiguous amino acid triplets. The frequencies of these motifs were used to train the linear programming boosting (LPBoost) algorithm LPBoost to classify between TCR repertoires. LPBoost could distinguish between the two classes of repertoire with accuracies above 80%, using a small subset of triplet sequences present at defined positions along the CDR3. The results suggest a model in which such motifs confer degenerate antigen specificity in the context of a highly diverse and largely private set of T cell receptors

The new very small angle neutron scattering spectrometer at Laboratoire Leon Brillouin

The design and characteristics of the new very small angle neutron scattering spectrometer under construction at the Laboratoire Leon Brillouin is described. Its goal is to extend the range of scattering vectors magnitudes towards 2x10{-4} /A. The unique feature of this new spectrometer is a high resolution two dimensional image plate detector sensitive to neutrons. The wavelength selection is achieved by a double reflection supermirror monochromator and the collimator uses a novel multibeam design

Feature selection using a one dimensional naïve Bayes' classifier increases the accuracy of support vector machine classification of CDR3 repertoires.

MOTIVATION: Somatic DNA recombination, the hallmark of vertebrate adaptive immunity, has the potential to generate a vast diversity of antigen receptor sequences. How this diversity captures antigen specificity remains incompletely understood. In this study we use high throughput sequencing to compare the global changes in T cell receptor β chain complementarity determining region 3 (CDR3β) sequences following immunization with ovalbumin administered with complete Freund's adjuvant (CFA) or CFA alone. RESULTS: The CDR3β sequences were deconstructed into short stretches of overlapping contiguous amino acids. The motifs were ranked according to a one-dimensional Bayesian classifier score comparing their frequency in the repertoires of the two immunization classes. The top ranking motifs were selected and used to create feature vectors which were used to train a support vector machine. The support vector machine achieved high classification scores in a leave-one-out validation test reaching  : >90% in some cases. SUMMARY: The study describes a novel two-stage classification strategy combining a one-dimensional Bayesian classifier with a support vector machine. Using this approach we demonstrate that the frequency of a small number of linear motifs three amino acids in length can accurately identify a CD4 T cell response to ovalbumin against a background response to the complex mixture of antigens which characterize Complete Freund's Adjuvant. AVAILABILITY AND IMPLEMENTATION: The sequence data is available at www.ncbi.nlm.nih.gov/sra/?term¼SRP075893 The Decombinator package is available at github.com/innate2adaptive/Decombinator The R package e1071 is available at the CRAN repository https://cran.r-project.org/web/packages/e1071/index.html CONTACT: [email protected] information: Supplementary data are available at Bioinformatics online

The planform mobility of river channel confluences: Insights from analysis of remotely sensed imagery

This is the author accepted manuscript. The final version is available from the publisher via the DOI in this record.River channel confluences are widely acknowledged as important geomorphological nodes that control the downstream routing of water and sediment, and which are locations for the preservation of thick fluvial deposits overlying a basal scour. Despite their importance, there has been little study of the stratigraphic characteristics of river junctions, or the role of confluence morphodynamics in influencing stratigraphic character and preservation potential. As a result, although it is known that confluences can migrate through time, models of confluence geomorphology and sedimentology are usually presented from the perspective that the confluence remains at a fixed location. This is problematic for a number of reasons, not least of which is the continuing debate over whether it is possible to discriminate between scour that has been generated by autocyclic processes (such as confluence scour) and that driven by allocyclic controls (such as sea-level change). This paper investigates the spatial mobility of river confluences by using the 40-year record of Landsat Imagery to elucidate the styles, rates of change and areal extent over which large river confluence scours may migrate. On the basis of these observations, a new classification of the types of confluence scour is proposed and applied to the Amazon and Ganges-Brahmaputra-Meghna (GBM) basins. This analysis demonstrates that the drivers of confluence mobility are broadly the same as those that drive channel change more generally. Thus in the GBM basin, a high sediment supply, large variability in monsoonal driven discharge and easily erodible bank materials result in a catchment where over 80% of large confluences are mobile over this 40-year window; conversely this figure is < 40% for the Amazon basin. These results highlight that: i) the potential areal extent of confluence scours is much greater than previously assumed, with the location of some confluences on the Jamuna (Brahmaputra) River migrating over a distance of 20 times the tributary channel width; ii) extensive migration in the confluence location is more common than currently assumed, and iii) confluence mobility is often tied to the lithological and hydrological characteristics of the drainage basins that determine sediment yield.This work was funded by NERC grant NE/I023228/1 to Sambrook Smith, Bull, Nicholas and Best

Quantitative Characterization of the T Cell Receptor Repertoire of Naive 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

Cognitive impairment induced by delta9-tetrahydrocannabinol occurs through heteromers between cannabinoid CB1 and serotonin 5-HT2A receptors

Delta-9-tetrahydrocannabinol (THC), the main psychoactive compound of marijuana, induces numerous undesirable effects, including memory impairments, anxiety, and dependence. Conversely, THC also has potentially therapeutic effects, including analgesia, muscle relaxation, and neuroprotection. However, the mechanisms that dissociate these responses are still not known. Using mice lacking the serotonin receptor 5-HT2A, we revealed that the analgesic and amnesic effects of THC are independent of each other: while amnesia induced by THC disappears in the mutant mice, THC can still promote analgesia in these animals. In subsequent molecular studies, we showed that in specific brain regions involved in memory formation, the receptors for THC and the 5-HT2A receptors work together by physically interacting with each other. Experimentally interfering with this interaction prevented the memory deficits induced by THC, but not its analgesic properties. Our results highlight a novel mechanism by which the beneficial analgesic properties of THC can be dissociated from its cognitive side effects