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

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

The sedimentology of river confluences

This is the author accepted manuscript. The final version is available from Wiley via the DOI in this record.Channel confluences are key nodes within large river networks, and yet surprisingly little is known about their spatial and temporal evolution. Moreover, because confluences are associated with vertical scour that typically extends to several times the mean channel depth, the deposits associated with such scours should have a high preservation potential within the rock record. Paradoxically, such scours are rarely observed, and their preservation and sedimentological interpretation are poorly understood. The present study details results from a physically‐based morphodynamic model that is applied to simulate the evolution and alluvial architecture of large river junctions. Boundary conditions within the model were defined to approximate the junction of the Ganges and Jamuna rivers, Bangladesh, with the model output being supplemented by geophysical datasets collected at this junction. The numerical simulations reveal several distinct styles of sedimentary fill that are related to the morphodynamic behaviour of bars, confluence scour downstream of braid bars, bend scour and major junction scour. Comparison with existing, largely qualitative, conceptual models reveals that none of these can be applied simply, although elements of each are evident in the deposits generated by the numerical simulation and observed in the geophysical data. The characteristics of the simulated scour deposits are found to vary according to the degree of reworking caused by channel migration, a factor not considered adequately in current conceptual models of confluence sedimentology. The alluvial architecture of major junction scours is thus characterized by the prevalence of erosion surfaces in conjunction with the thickest depositional sets. Confluence scour downstream of braid bar and bend scour sites may preserve some large individual sets, but these locations are typically characterized by lower average set thickness compared to major junction scour and by a lack of large‐scale erosional surfaces. Areas of deposition not related to any of the specific scour types highlighted above record the thinnest depositional sets. This variety in the alluvial architecture of scours may go some way towards explaining the paradox of ancient junction scours, that while abundant large scours are likely in the rock record, they have been reported rarely. The present results outline the likely range of confluence sedimentology and will serve as a new tool for recognizing and interpreting these deposits in the ancient fluvial record.This work was funded by a UK Natural Environment Research Council award to Sambrook Smith (NE/I023228/1), Bull (NE/I023864/1) and Nicholas (NE/I023120/1)

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

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

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

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 inevitable youthfulness of known high-redshift radio galaxies

Radio galaxies can be seen out to very high redshifts, where in principle they can serve as probes of the early evolution of the Universe. Here we show that for any model of radio-galaxy evolution in which the luminosity decreases with time after an initial rapid increase (that is, essentially all reasonable models), all observable high-redshift radio-galaxies must be seen when the lobes are less than 10^7 years old. This means that high-redshift radio galaxies can be used as a high-time-resolution probe of evolution in the early Universe. Moreover, this result helps to explain many observed trends of radio-galaxy properties with redshift [(i) the `alignment effect' of optical emission along radio-jet axes, (ii) the increased distortion in radio structure, (iii) the decrease in physical sizes, (iv) the increase in radio depolarisation, and (v) the increase in dust emission] without needing to invoke explanations based on cosmology or strong evolution of the surrounding intergalactic medium with cosmic time, thereby avoiding conflict with current theories of structure formation.Comment: To appear in Nature. 4 pages, 2 colour figures available on request. Also available at http://www-astro.physics.ox.ac.uk/~km