Identification of errors introduced during high throughput sequencing of the T cell receptor repertoire

Quantitative assessment of T cell repertoire recovery after hematopoietic stem cell transplantation

Delayed T-cell recovery and restricted T-cell receptor (TCR) diversity after allogeneic hematopoietic stem cell transplantation (allo-HSCT) are associated with increased risks of infection and cancer relapse. Technical challenges have limited faithful measurement of TCR diversity following allo-HSCT. Here we combined 5′-RACE PCR with deep sequencing, to quantify TCR diversity in 28 allo-HSCT recipients using a single oligonucleotide pair. Analysis of duplicate blood samples confirmed that the frequency of individual TCRs was accurately determined. After 6 months, cord blood graft recipients approximated the TCR diversity of healthy individuals, whereas recipients of T-cell-depleted peripheral blood stem cell grafts had a 28-fold and 14-fold lower CD4⁺ and CD8⁺ T-cell diversity, respectively. After 12 months, these deficiencies had improved for the CD4⁺, but not the CD8⁺ T-cell compartment. Overall, this method provides unprecedented views of T-cell repertoire recovery after allo-HSCT and may identify patients at high risk of infection or relapse.

Generalized DNA barcode design based on Hamming codes

The diversity and scope of multiplex parallel sequencing applications is steadily increasing. Critically, multiplex parallel sequencing applications methods rely on the use of barcoded primers for sample identification, and the quality of the barcodes directly impacts the quality of the resulting sequence data. Inspection of the recent publications reveals a surprisingly variable quality of the barcodes employed. Some barcodes are made in a semi empirical fashion, without quantitative consideration of error correction or minimal distance properties. After systematic comparison of published barcode sets, including commercially distributed barcoded primers from Illumina and Epicentre, methods for improved, Hamming code-based sequences are suggested and illustrated. Hamming barcodes can be employed for DNA tag designs in many different ways while preserving minimal distance and error-correcting properties. In addition, Hamming barcodes remain flexible with regard to essential biological parameters such as sequence redundancy and GC content. Wider adoption of improved Hamming barcodes is encouraged in multiplex parallel sequencing applications.

Rep-Seq: uncovering the immunological repertoire through next-generation sequencing

Recent scientific discoveries fuelled by the application of next-generation DNA and RNA sequencing technologies highlight the striking impact of these platforms in characterizing multiple aspects in genomics research. This technology has been used in the study of the B-cell and T-cell receptor repertoire. The novelty of immunosequencing comes from the recent rapid development of techniques and the exponential reduction in cost of sequencing. Here, we describe some of the technologies, which we collectively refer to as Rep-Seq (repertoire sequencing), to portray achievements in the field and to present the essential and inseparable role of next-generation sequencing to the understanding of entities in immune response. The large Rep-Seq data sets that should be available in the near future call for new computational algorithms to segue the transition from 'classic' molecular-based analysis to system-wide analysis. The combination of new algorithms with high-throughput data will form the basis for possible new clinical implications in personalized medicine and deeper understanding of immune behaviour and immune response.

Highly diverse TCRα chain repertoire of pre-immune CD8⁺ T cells reveals new insights in gene recombination

Although the T-cell receptor αδ (TCRαδ) locus harbours large libraries of variable (TRAV) and junctional (TRAJ) gene segments, according to previous studies the TCRα chain repertoire is of limited diversity due to restrictions imposed by sequential coordinate TRAV-TRAJ recombinations. By sequencing tens of millions of TCRα chain transcripts from naive mouse CD8(+) T cells, we observed a hugely diverse repertoire, comprising nearly all possible TRAV-TRAJ combinations. Our findings are not compatible with sequential coordinate gene recombination, but rather with a model in which contraction and DNA looping in the TCRαδ locus provide equal access to TRAV and TRAJ gene segments, similarly to that demonstrated for IgH gene recombination. Generation of the observed highly diverse TCRα chain repertoire necessitates deletion of failed attempts by thymic-positive selection and is essential for the formation of highly diverse TCRαβ repertoires, capable of providing good protective immunity.

High-throughput VDJ sequencing for quantification of minimal residual disease in chronic lymphocytic leukemia and immune reconstitution assessment

The primary cause of poor outcome following allogeneic hematopoietic cell transplantation (HCT) for chronic lymphocytic leukemia (CLL) is disease recurrence. Detection of increasing minimal residual disease (MRD) following HCT may permit early intervention to prevent clinical relapse; however, MRD quantification remains an uncommon diagnostic test because of logistical and financial barriers to widespread use. Here we describe a method for quantifying CLL MRD using widely available consensus primers for amplification of all Ig heavy chain (IGH) genes in a mixture of peripheral blood mononuclear cells, followed by high-throughput sequencing (HTS) for disease-specific IGH sequence quantification. To achieve accurate MRD quantification, we developed a systematic bioinformatic methodology to aggregate cancer clone sequence variants arising from systematic and random artifacts occurring during IGH-HTS. We then compared the sensitivity of IGH-HTS, flow cytometry, and allele-specific oligonucleotide PCR for MRD quantification in 28 samples collected from 6 CLL patients following allogeneic HCT. Using amplimer libraries generated with consensus primers from patient blood samples, we demonstrate the sensitivity of IGH-HTS with 454 pyrosequencing to be 10(-5), with a high correlation between quantification by allele-specific oligonucleotide PCR and IGH-HTS (r = 0.85). From the same dataset used to quantify MRD, IGH-HTS also allowed us to profile IGH repertoire reconstitution after HCT-information not provided by the other MRD methods. IGH-HTS using consensus primers will broaden the availability of MRD quantification in CLL and other B cell malignancies, and this approach has potential for quantitative evaluation of immune diversification following transplant and nontransplant therapies.

Naive antibody gene-segment frequencies are heritable and unaltered by chronic lymphocyte ablation

A diverse antibody repertoire is essential for an effective adaptive immune response to novel molecular surfaces. Although past studies have observed common patterns of V-segment use, as well as variation in V-segment use between individuals, the relative contributions to variance from genetics, disease, age, and environment have remained unclear. Using high-throughput sequence analysis of monozygotic twins, we show that variation in naive V(H) and D(H) segment use is strongly determined by an individual's germ-line genetic background. The inherited segment-use profiles are resilient to differential environmental exposure, disease processes, and chronic lymphocyte depletion therapy. Signatures of the inherited profiles were observed in class switched germ-line use of each individual. However, despite heritable segment use, the rearranged complementarity-determining region-H3 repertoires remained highly specific to the individual. As it has been previously demonstrated that certain V-segments exhibit biased representation in autoimmunity, lymphoma, and viral infection, we anticipate our findings may provide a unique mechanism for stratifying individual risk profiles in specific diseases.

454 antibody sequencing - error characterization and correction

BACKGROUND

454 sequencing is currently the method of choice for sequencing of antibody repertoires and libraries containing large numbers (106 to 1012) of different molecules with similar frameworks and variable regions which poses significant challenges for identifying sequencing errors. Identification and correction of sequencing errors in such mixtures is especially important for the exploration of complex maturation pathways and identification of putative germline predecessors of highly somatically mutated antibodies. To quantify and correct errors incorporated in 454 antibody sequencing, we sequenced six antibodies at different known concentrations twice over and compared them with the corresponding known sequences as determined by standard Sanger sequencing.

RESULTS

We found that 454 antibody sequencing could lead to approximately 20% incorrect reads due to insertions that were mostly found at shorter homopolymer regions of 2-3 nucleotide length, and less so by insertions, deletions and other variants at random sites. Correction of errors might reduce this population of erroneous reads down to 5-10%. However, there are a certain number of errors accounting for 4-8% of the total reads that could not be corrected unless several repeated sequencing is performed, although this may not be possible for large diverse libraries and repertoires including complete sets of antibodies (antibodyomes).

CONCLUSIONS

The experimental test procedure carried out for assessing 454 antibody sequencing errors reveals high (up to 20%) incorrect reads; the errors can be reduced down to 5-10% but not less which suggests the use of caution to avoid false discovery of antibody variants and diversity.

High throughput cellular screens to interrogate the human T and B cell repertoires

There is growing consensus for an urgent need to develop efficient methods to analyze the human T and B cell response in order to understand its regulation, to unravel host-pathogen interactions and to design better vaccines. In this review we will consider different methods that have been developed to analyze human T and B cells and focus in particular on those based on high-throughput cell culture, such as memory B cell immortalization, plasma cell cultures and libraries of naïve or memory T cells. These culture-based methods can be used to interrogate with high efficiency the immune repertoires in order to characterize the frequency, function and specificity of naïve or memory T cells and to isolate potent and broadly neutralizing monoclonal antibodies.

Building a model: developing genomic resources for common milkweed (Asclepias syriaca) with low coverage genome sequencing

Background

Milkweeds (Asclepias L.) have been extensively investigated in diverse areas of evolutionary biology and ecology; however, there are few genetic resources available to facilitate and compliment these studies. This study explored how low coverage genome sequencing of the common milkweed (Asclepias syriaca L.) could be useful in characterizing the genome of a plant without prior genomic information and for development of genomic resources as a step toward further developing A. syriaca as a model in ecology and evolution.

Results

A 0.5× genome of A. syriaca was produced using Illumina sequencing. A virtually complete chloroplast genome of 158,598 bp was assembled, revealing few repeats and loss of three genes: accD, clpP, and ycf1. A nearly complete rDNA cistron (18S-5.8S-26S; 7,541 bp) and 5S rDNA (120 bp) sequence were obtained. Assessment of polymorphism revealed that the rDNA cistron and 5S rDNA had 0.3% and 26.7% polymorphic sites, respectively. A partial mitochondrial genome sequence (130,764 bp), with identical gene content to tobacco, was also assembled. An initial characterization of repeat content indicated that Ty1/copia-like retroelements are the most common repeat type in the milkweed genome. At least one A. syriaca microread hit 88% of Catharanthus roseus (Apocynaceae) unigenes (median coverage of 0.29×) and 66% of single copy orthologs (COSII) in asterids (median coverage of 0.14×). From this partial characterization of the A. syriaca genome, markers for population genetics (microsatellites) and phylogenetics (low-copy nuclear genes) studies were developed.

Conclusions

The results highlight the promise of next generation sequencing for development of genomic resources for any organism. Low coverage genome sequencing allows characterization of the high copy fraction of the genome and exploration of the low copy fraction of the genome, which facilitate the development of molecular tools for further study of a target species and its relatives. This study represents a first step in the development of a community resource for further study of plant-insect co-evolution, anti-herbivore defense, floral developmental genetics, reproductive biology, chemical evolution, population genetics, and comparative genomics using milkweeds, and A. syriaca in particular, as ecological and evolutionary models.