Scott Hunicke-Smith

High-throughput sequencing of the paired human immunoglobulin heavy and light chain repertoire

Each B-cell receptor consists of a pair of heavy and light chains. High-throughput sequencing can identify large numbers of heavy- and light-chain variable regions (VH and VL) in a given B-cell repertoire, but information about endogenous pairing of heavy and light chains is lost after bulk lysis of B-cell populations. Here we describe a way to retain this pairing information. In our approach, single B cells (>5 × 104 capacity per experiment) are deposited in a high-density microwell plate (125 pl/well) and lysed in situ. mRNA is then captured on magnetic beads, reverse transcribed and amplified by emulsion VH:VL linkage PCR. The linked transcripts are analyzed by Illumina high-throughput sequencing. We validated the fidelity of VH:VL pairs identified by this approach and used the method to sequence the repertoire of three human cell subsets—peripheral blood IgG+ B cells, peripheral plasmablasts isolated after tetanus toxoid immunization and memory B cells isolated after seasonal influenza vaccination.

Monoclonal antibodies isolated without screening by analyzing the variable-gene repertoire of plasma cells

Isolation of antigen-specific monoclonal antibodies (mAbs) and antibody fragments relies on high-throughput screening of immortalized B cells or recombinant antibody libraries. We bypassed the screening step by using high-throughput DNA sequencing and bioinformatic analysis to mine antibody variable region (V)-gene repertoires from bone marrow plasma cells (BMPC) of immunized mice. BMPCs, which cannot be immortalized, produce the vast majority of circulating antibodies. We found that the V-gene repertoire of BMPCs becomes highly polarized after immunization, with the most abundant sequences represented at frequencies between ∼1% and >10% of the total repertoire. We paired the most abundant variable heavy (VH) and variable light (VL) genes based on their relative frequencies, reconstructed them using automated gene synthesis, and expressed recombinant antibodies in bacteria or mammalian cells. Antibodies generated in this manner from six mice, each immunized with one of three antigens were overwhelmingly antigen specific (21/27 or 78%). Those generated from a mouse with high serum titers had nanomolar binding affinities.

Long-term neural and physiological phenotyping of a single human

Psychiatric disorders are characterized by major fluctuations in psychological function over the course of weeks and months, but the dynamic characteristics of brain function over this timescale in healthy individuals are unknown. Here, as a proof of concept to address this question, we present the MyConnectome project. An intensive phenome-wide assessment of a single human was performed over a period of 18 months, including functional and structural brain connectivity using magnetic resonance imaging, psychological function and physical health, gene expression and metabolomics. A reproducible analysis workflow is provided, along with open access to the data and an online browser for results. We demonstrate dynamic changes in brain connectivity over the timescales of days to months, and relations between brain connectivity, gene expression and metabolites. This resource can serve as a testbed to study the joint dynamics of human brain and metabolic function over time, an approach that is critical for the development of precision medicine strategies for brain disorders.

Rapid Construction and Characterization of Synthetic Antibody Libraries Without DNA Amplification

We report on a simple method to rapidly generate very large libraries of genes encoding mutant proteins without the use of DNA amplification, and the application of this methodology in the construction of synthetic immunoglobulin variable heavy (VH) and light (Vκ) libraries. Four high quality, chemically synthesized polynucleotides (90–140 bases) were annealed and extended using T4 DNA polymerase. Following electroporation, >109 transformants could be synthesized within 1 day. Fusion to β‐lactamase and selection on ampicillin resulted in 3.7 × 108 VH and 6.9 × 108 Vκ clones highly enriched for full‐length, in‐frame genes. High‐throughput 454 DNA sequencing of >250,000 VH and Vκ genes from the pre‐ and post‐selection libraries revealed that, in addition to the expected reduction in reading‐frame shifts and stop codons, selection for functional expression also resulted in a statistical decrease in the cysteine content. Apart from these differences, there was a good agreement between the expected and actual diversity, indicating that neither oligonucleotide synthesis nor biological constrains due to protein synthesis of VH/Vκ‐β‐lactamase fusions introduce biases in the amino acid composition of the randomized regions. This methodology can be employed for the rapid construction of highly diverse libraries with the near elimination of PCR errors in invariant regions. Biotechnol. Bioeng. 2010; 106: 347–357.

A capture-sequencing strategy identifies IRF8, EBF1, and APRIL as novel IGH fusion partners in B-cell lymphoma

The characterization of immunoglobulin heavy chain (IGH) translocations provides information on the diagnosis and guides therapeutic decisions in mature B-cell malignancies while enhancing our understanding of normal and malignant B-cell biology. However, existing methodologies for the detection of IGH translocations are labor intensive, often require viable cells, and are biased toward known IGH fusions. To overcome these limitations, we developed a capture sequencing strategy for the identification of IGH rearrangements at nucleotide level resolution and tested its capabilities as a diagnostic and discovery tool in 78 primary diffuse large B-cell lymphomas (DLBCLs). We readily identified IGH-BCL2, IGH-BCL6, IGH-MYC, and IGH-CCND1 fusions and discovered IRF8, EBF1, and TNFSF13 (APRIL) as novel IGH partners in these tumors. IRF8 and TNFSF13 expression was significantly higher in lymphomas with IGH rearrangements targeting these loci. Modeling the deregulation of IRF8 and EBF1 in vitro defined a lymphomagenic profile characterized by up-regulation of AID and/or BCL6, down-regulation of PRMD1, and resistance to apoptosis. Using a capture sequencing strategy, we discovered the B-cell relevant genes IRF8, EBF1, and TNFSF13 as novel targets for IGH deregulation. This methodology is poised to change how IGH translocations are identified in clinical settings while remaining a powerful tool to uncover the pathogenesis of B-cell malignancies.

High-throughput sequencing of human plasma RNA by using thermostable group II intron reverse transcriptases.

Next-generation RNA-sequencing (RNA-seq) has revolutionized transcriptome profiling, gene expression analysis, and RNA-based diagnostics. Here, we developed a new RNA-seq method that exploits thermostable group II intron reverse transcriptases (TGIRTs) and used it to profile human plasma RNAs. TGIRTs have higher thermostability, processivity, and fidelity than conventional reverse transcriptases, plus a novel template-switching activity that can efficiently attach RNA-seq adapters to target RNA sequences without RNA ligation. The new TGIRT-seq method enabled construction of RNA-seq libraries from <1 ng of plasma RNA in <5 h. TGIRT-seq of RNA in 1-mL plasma samples from a healthy individual revealed RNA fragments mapping to a diverse population of protein-coding gene and long ncRNAs, which are enriched in intron and antisense sequences, as well as nearly all known classes of small ncRNAs, some of which have never before been seen in plasma. Surprisingly, many of the small ncRNA species were present as full-length transcripts, suggesting that they are protected from plasma RNases in ribonucleoprotein (RNP) complexes and/or exosomes. This TGIRT-seq method is readily adaptable for profiling of whole-cell, exosomal, and miRNAs, and for related procedures, such as HITS-CLIP and ribosome profiling.

Next‐Generation Sequencing RNA‐Seq Library Construction

This unit presents protocols for construction of next‐generation sequencing (NGS) directional RNA sequencing libraries for the Illumina HiSeq and MiSeq from a wide variety of input RNA sources. The protocols are based on the New England Biolabs (NEB) small RNA library preparation set for Illumina, although similar kits exist from different vendors. The protocol preserves the orientation of the original RNA in the final sequencing library, enabling strand‐specific analysis of the resulting data. These libraries have been used for differential gene expression analysis and small RNA discovery and are currently being tested for de novo transcriptome assembly. The protocol is robust and applicable to a broad range of RNA input types and RNA quality, making it ideal for high‐throughput laboratories. Curr. Protoc. Mol. Biol. 106:4.21.1‐4.21.19.

Next‐Generation Sequencing Fragment Library Construction

All current next‐generation sequencing (NGS) platforms and applications require the sequencing library to have specific characteristics: in particular, size, size distribution, and 5′ and 3′ flanking sequences. This unit presents a robust protocol for converting a wide variety of input DNA samples into appropriate NGS libraries and discusses important considerations in experimental design, failure modes, and typical results. Curr. Protoc. Mol. Biol. 107:7.17.1‐7.17.16.