Error analysis of deep sequencing of phage libraries: peptides censored in sequencing

Structure, Biology, and Applications of Filamentous Bacteriophages

The closely related Escherichia coli Ff filamentous phages (f1, fd, and M13) have taken a fantastic journey over the past 60 years, from the urban sewerage from which they were first isolated, to their use in high-end technologies in multiple fields. Their relatively small genome size, high titers, and the virions that tolerate fusion proteins make the Ffs an ideal system for phage display. Folding of the fusions in the oxidizing environment of the E. coli periplasm makes the Ff phages a platform that allows display of eukaryotic surface and secreted proteins, including antibodies. Resistance of the Ffs to a broad range of pH and detergents facilitates affinity screening in phage display, whereas the stability of the virions at ambient temperature makes them suitable for applications in material science and nanotechnology. Among filamentous phages, only the Ffs have been used in phage display technology, because of the most advanced state of knowledge about their biology and the various tools developed for E. coli as a cloning host for them. Filamentous phages have been thought to be a rather small group, infecting mostly Gram-negative bacteria. A recent discovery of more than 10 thousand diverse filamentous phages in bacteria and archaea, however, opens a fascinating prospect for novel applications. The main aim of this review is to give detailed biological and structural information to researchers embarking on phage display projects. The secondary aim is to discuss the yet-unresolved puzzles, as well as recent developments in filamentous phage biology, from a viewpoint of their impact on current and future applications.

Restriction of the Global IgM Repertoire in Antiphospholipid Syndrome

The typical anti-phospholipid antibodies (APLA) in the anti-phospholipid syndrome (APS) are reactive with the phospholipid-binding protein β2GPI as well as a growing list of other protein targets. The relation of APLA to natural antibodies and the fuzzy set of autoantigens involved provoked us to study the changes in the IgM repertoire in APS. To this end, peptides selected by serum IgM from a 7-residue linear peptide phage display library (PDL) were deep sequenced. The analysis was aided by a novel formal representation of the Igome (the mimotope set reflecting the IgM specificities) in the form of a sequence graph. The study involved women with APLA and habitual abortions (n=24) compared to age-matched clinically healthy pregnant women (n=20). Their pooled Igomes (297 028 mimotope sequences) were compared also to the global public repertoire Igome of pooled donor plasma IgM (n=2 796 484) and a set of 7-mer sequences found in the J regions of human immunoglobulins (n=4 433 252). The pooled Igome was represented as a graph connecting the sequences as similar as the mimotopes of the same monoclonal antibody. The criterion was based on previously published data. In the resulting graph, identifiable clusters of vertices were considered related to the footprints of overlapping antibody cross-reactivities. A subgraph based on the clusters with a significant differential expression of APS patients' mimotopes contained predominantly specificities underrepresented in APS. The differentially expressed IgM footprints showed also an increased cross-reactivity with immunoglobulin J regions. The specificities underexpressed in APS had a higher correlation with public specificities than those overexpressed. The APS associated specificities were strongly related also to the human peptidome with 1 072 mimotope sequences found in 7 519 human proteins. These regions were characterized by low complexity. Thus, the IgM repertoire of the APS patients was found to be characterized by a significant reduction of certain public specificities found in the healthy controls with targets representing low complexity linear self-epitopes homologous to human antibody J regions.

Expanding the Chemical Diversity of Genetically Encoded Libraries

The power of ribosomes has increasingly been harnessed for the synthesis and selection of molecular libraries. Technologies, such as phage display, yeast display, and mRNA display, effectively couple genotype to phenotype for the molecular evolution of high affinity epitopes for many therapeutic targets. Genetic code expansion is central to the success of these technologies, allowing researchers to surpass the intrinsic capabilities of the ribosome and access new, genetically encoded materials for these selections. Here, we review techniques for the chemical expansion of genetically encoded libraries, their abilities and limits, and opportunities for further development. Importantly, we also discuss methods and metrics used to assess the efficiency of modification and library diversity with these new techniques.

Compositional Bias in Naïve and Chemically-modified Phage-Displayed Libraries uncovered by Paired-end Deep Sequencing

Understanding the composition of a genetically-encoded (GE) library is instrumental to the success of ligand discovery. In this manuscript, we investigate the bias in GE-libraries of linear, macrocyclic and chemically post-translationally modified (cPTM) tetrapeptides displayed on the M13KE platform, which are produced via trinucleotide cassette synthesis (19 codons) and NNK-randomized codon. Differential enrichment of synthetic DNA {S}, ligated vector {L} (extension and ligation of synthetic DNA into the vector), naïve libraries {N} (transformation of the ligated vector into the bacteria followed by expression of the library for 4.5 hours to yield a "naïve" library), and libraries chemically modified by aldehyde ligation and cysteine macrocyclization {M} characterized by paired-end deep sequencing, detected a significant drop in diversity in {L} → {N}, but only a minor compositional difference in {S} → {L} and {N} → {M}. Libraries expressed at the N-terminus of phage protein pIII censored positively charged amino acids Arg and Lys; libraries expressed between pIII domains N1 and N2 overcame Arg/Lys-censorship but introduced new bias towards Gly and Ser. Interrogation of biases arising from cPTM by aldehyde ligation and cysteine macrocyclization unveiled censorship of sequences with Ser/Phe. Analogous analysis can be used to explore library diversity in new display platforms and optimize cPTM of these libraries.

Filamentous Phage: Structure and Biology

Ff filamentous phage (fd, M13 and f1) of Escherichia coli have been the workhorse of phage display technology for the past 30 years. Dominance of Ff over other bacteriophage in display technology stems from the titres that are about 100-fold higher than any other known phage, efficacious transformation ensuring large library size and superior stability of the virion at high temperatures, detergents and pH extremes, allowing broad range of biopanning conditions in screening phage display libraries. Due to the excellent understanding of infection and assembly requirements, Ff phage have also been at the core of phage-assisted continual protein evolution strategies (PACE). This chapter will give an overview of the Ff filamentous phage structure and biology, emphasizing those properties of the Ff phage life cycle and virion that are pertinent to phage display applications.

PHASTpep: Analysis Software for Discovery of Cell-Selective Peptides via Phage Display and Next-Generation Sequencing

Next-generation sequencing has enhanced the phage display process, allowing for the quantification of millions of sequences resulting from the biopanning process. In response, many valuable analysis programs focused on specificity and finding targeted motifs or consensus sequences were developed. For targeted drug delivery and molecular imaging, it is also necessary to find peptides that are selective—targeting only the cell type or tissue of interest. We present a new analysis strategy and accompanying software, PHage Analysis for Selective Targeted PEPtides (PHASTpep), which identifies highly specific and selective peptides. Using this process, we discovered and validated, both in vitro and in vivo in mice, two sequences (HTTIPKV and APPIMSV) targeted to pancreatic cancer-associated fibroblasts that escaped identification using previously existing software. Our selectivity analysis makes it possible to discover peptides that target a specific cell type and avoid other cell types, enhancing clinical translatability by circumventing complications with systemic use.

Facile, high quality sequencing of bacterial genomes from small amounts of DNA

Sequencing bacterial genomes has traditionally required large amounts of genomic DNA (~1 μg). There have been few studies to determine the effects of the input DNA amount or library preparation method on the quality of sequencing data. Several new commercially available library preparation methods enable shotgun sequencing from as little as 1 ng of input DNA. In this study, we evaluated the NEBNext Ultra library preparation reagents for sequencing bacterial genomes. We have evaluated the utility of NEBNext Ultra for resequencing and de novo assembly of four bacterial genomes and compared its performance with the TruSeq library preparation kit. The NEBNext Ultra reagents enable high quality resequencing and de novo assembly of a variety of bacterial genomes when using 100 ng of input genomic DNA. For the two most challenging genomes (Burkholderia spp.), which have the highest GC content and are the longest, we also show that the quality of both resequencing and de novo assembly is not decreased when only 10 ng of input genomic DNA is used.

Identification of target-binding peptide motifs by high-throughput sequencing of phage-selected peptides

High-throughput sequencing was previously applied to phage-selected peptides in order to gain insight into the abundance and diversity of isolated peptides. Herein we developed a procedure to efficiently compare the sequences of large numbers of phage-selected peptides for the purpose of identifying target-binding peptide motifs. We applied the procedure to analyze bicyclic peptides isolated against five different protein targets: sortase A, urokinase-type plasminogen activator, coagulation factor XII, plasma kallikrein and streptavidin. We optimized sequence data filters to reduce biases originating from the sequencing method and developed sequence correction algorithms to prevent identification of false consensus motifs. With our strategy, we were able to identify rare target-binding peptide motifs, as well as to define more precisely consensus sequences and sub-groups of consensus sequences. This information is valuable to choose peptide leads for drug development and it facilitates identification of epitopes. We furthermore show that binding motifs can be identified after a single round of phage selection. Such a selection regimen reduces propagation-related bias and may facilitate application of phage display in non-specialized laboratories, as procedures such as bacterial infection, phage propagation and purification are not required.

Impact of next-generation sequencing error on analysis of barcoded plasmid libraries of known complexity and sequence

Barcoded vectors are promising tools for investigating clonal diversity and dynamics in hematopoietic gene therapy. Analysis of clones marked with barcoded vectors requires accurate identification of potentially large numbers of individually rare barcodes, when the exact number, sequence identity and abundance are unknown. This is an inherently challenging application, and the feasibility of using contemporary next-generation sequencing technologies is unresolved. To explore this potential application empirically, without prior assumptions, we sequenced barcode libraries of known complexity. Libraries containing 1, 10 and 100 Sanger-sequenced barcodes were sequenced using an Illumina platform, with a 100-barcode library also sequenced using a SOLiD platform. Libraries containing 1 and 10 barcodes were distinguished from false barcodes generated by sequencing error by a several log-fold difference in abundance. In 100-barcode libraries, however, expected and false barcodes overlapped and could not be resolved by bioinformatic filtering and clustering strategies. In independent sequencing runs multiple false-positive barcodes appeared to be represented at higher abundance than known barcodes, despite their confirmed absence from the original library. Such errors, which potentially impact barcoding studies in an application-dependent manner, are consistent with the existence of both stochastic and systematic error, the mechanism of which is yet to be fully resolved.