Aligning,analyzing, and visualizing sequences for antibody engineering: Automated recognition of immunoglobulin variable region features

The analysis and comparison of large numbers of immunoglobulin (Ig) sequences that arise during an antibody selection campaign can be time‐consuming and tedious. Typically, the identification and annotation of framework as well as complementarity‐determining regions (CDRs) is based on multiple sequence alignments using standardized numbering schemes, which allow identification of equivalent residues among different family members but often necessitate expert knowledge and manual intervention. Moreover, due to the enormous length variability of some CDRs the benefit of conventional Ig numbering schemes is limited and the calculation of correct sequence alignments can become challenging. Whereas, in principle, a well established set of rules permits the assignment of CDRs from the amino acid sequence alone, no currently available sequence alignment editor provides an algorithm to annotate new Ig sequences accordingly. Here we present a unique pattern matching method implemented into our recently developed ANTIC ALIgN editor that automatically identifies all hypervariable and framework regions in experimentally elucidated antibody sequences using so‐called “regular expressions.” By combination of this widely supported software syntax with the unique capabilities of real‐time aligning, editing and analyzing extended sets of amino acid and/or nucleotide sequences simultaneously on a local workstation, ANTIC ALIgN provides a powerful utility for antibody engineering. Proteins 2016; 85:65–71. © 2016 Wiley Periodicals, Inc.     

Father-to-Mother-to-Infant Transmission of HIV-1: Clonally Transmitted Isolate of Infant Mutates More Rapidly than That of the Mother and Rapidly Loses Reactivity with Neutralizing Antibody

The sequences of the V3 loop and surrounding regions of human immunodeficiency virus type-1 from a father-to-mother-to-infant trimmer were studied and the horizontal and vertical transmissions compared. The father's virus was variable for reactivity with neutralizing antibody and sequences of the V3 loop central core sequence. In contrast, the mother's viral sequences were much less diverse and reacted with a virus neutralizing antibody. The infant's viral sequences were also less diverse than those of the father, and N-glycosylation sites were conserved. By phylogenetic analysis, the major clone, of which V3-peptide reacted with the neutralizing antibody, was found to be transmitted from the mother to her infant; however, the mutated minor clones did not bind to the antibody. These findings suggest that both horizontal and vertical virus transmission were selective, and that the clonally transmitted virus in infants mutates more rapidly than viruses in the mother, to whom the virus was horizontally transmitted.     

PipTools: a computational toolkit to annotate and analyze pairwise comparisons of genomic sequences

Sequence conservation between species is useful both for locating coding regions of genes and for identifying functional noncoding segments. Hence interspecies alignment of genomic sequences is an important computational technique. However, its utility is limited without extensive annotation. We describe a suite of software tools, PipTools, and related programs that facilitate the annotation of genes and putative regulatory elements in pairwise alignments. The alignment server PipMaker uses the output of these tools to display detailed information needed to interpret alignments. These programs are provided in a portable format for use on common desktop computers and both the toolkit and the PipMaker server can be found at our Web site (http://bio.cse.psu.edu/). We illustrate the utility of the toolkit using annotation of a pairwise comparison of the mouse MHC class II and class III regions with orthologous human sequences and subsequently identify conserved, noncoding sequences that are DNase I hypersensitive sites in chromatin of mouse cells.     

Modeling the feasibility of whole genome shotgun sequencing using a pairwise end strategy

In pairwise end sequencing, sequences are determined from both ends of random subclones derived from a DNA target. Sufficiently similar overlapping end sequences are identified and grouped into contigs. When a clone's paired end sequences fall in different contigs, the contigs are connected together to form scaffolds. Increasingly, the goals of pairwise strategies are large and highly repetitive genomic targets. Here, we consider large-scale pairwise strategies that employ mixtures of subclone sizes. We explore the properties of scaffold formation within a hybrid theory/simulation mathematical model of a genomic target that contains many repeat families. Using this model, we evaluate problems that may arise, such as falsely linked end sequences (due either to random matches or to homologous repeats) and scaffolds that terminate without extending the full length of the target. We illustrate our model with an exploration of a strategy for sequencing the human genome. Our results show that, for a strategy that generates 10-fold sequence coverage derived from the ends of clones ranging in length from 2 to 150 kb, using an appropriate rule for detecting overlaps, we expect few false links while obtaining a single scaffold extending the length of each chromosome.     

The heterochromatin of grasshoppers from the Caledia captiva species complex. I. Sequence evolution and conservation in a highly repeated DNA family

The restriction enzyme TaqI digests 0.2% of the genomic DNA from the grasshopper Caledia captiva to a family of sequences 168 bp in length (length of consensus sequence). The sequence variation of this "Taq family" of repeat units was examined among four races from C. captiva to assay the pattern of evolution within this highly repeated DNA. The Taq-family repeats are located in C-banded heterochromatin on at least one member of each homologous pair of chromosomes; the locations range from centromeric to telomeric. Thirty-nine cloned repeats isolated from two population 1A individuals along with 11 clones from seven populations taken from three of the races demonstrated sequence variation at 72 positions. Pairwise comparisons of the cloned repeats, both within an individual and between different races, indicate that levels of intraspecific divergence, as measured by reproductive incompatibility, do not correlate with sequence divergence among the 168-bp repeats. A number of subsequences within the repeat remain unchanged among all 50 clones; the longest of these is 18 bp. That the same 18-bp subsequence is present in all clones examined is a finding that departs significantly (P less than 0.01) from what would be expected to occur at random. Two other cloned repeats, from a reproductively isolated race of C. captiva, have sequences that show 56% identity with this 18-bp conserved region. An analysis showed that the frequency of occurrence of an RsaI recognition site within the 168- bp repeat in the entire Taq family agreed with that found in the cloned sequences. These data, along with a partial sequence for the entire Taq family obtained by sequencing uncloned repeats, suggest that the consensus sequence from the cloned copies is representative of this highly repeated family and is not a biased sample resulting from the cloning procedure. The 18-bp conserved sequence is part of a 42-bp sequence that possesses dyad symmetry typical of protein-binding sites. We speculate that this may be significant in the evolution of the Taq family of sequences.      

Nucleotide sequence determination and secondary structure of Xenopus U3 snRNA

Using a combination of RNA sequencing and construction of cDNA clones followed by DNA sequencing, we have determined the primary nucleotide sequence of U3 snRNA in Xenopus laevis and Xenopus borealis. This molecule has a length of 219 nucleotides. Alignment of the Xenopus sequences with U3 snRNA sequences from other organisms reveals three evolutionarily conserved blocks. We have probed the secondary structure of U3 snRNA in intact Xenopus laevis nuclei using single-strand specific chemical reagents; primer extension was used to map the positions of chemical modification. The three blocks of conserved sequences fall within single-stranded regions, and are therefore accessible for interaction with other molecules. Models of U3 snRNA function are discussed in light of these data.     

Construction and characterization of single-chain variable fragment antibody library derived from germline rearranged immunoglobulin variable genes

Antibody repertoires for library construction are conventionally harvested from mRNAs of immune cells. To examine whether germline rearranged immunoglobulin (Ig) variable region genes could be used as source of antibody repertoire, an immunized phage-displayed scFv library was prepared using splenocytic genomic DNA as template. In addition, a novel frame-shifting PCR (fsPCR) step was introduced to rescue stop codon and to enhance diversity of the complementarity-determining region 3 (CDR3). The germline scFv library was initially characterized against the hapten antigen phenyloxazolone (phOx). Sequence analysis of the phOx-selective scFvs indicated that the CDRs consisted of novel as well as conserved motifs. In order to illustrate that the diversity of CDR3 was increased by the fsPCR step, a second scFv library was constructed using a single scFv clone L3G7C as a template. Despite showing similar binding characteristics towards phOx, the scFv clones that were obtained from the L3G7C-derived antibody library gave a lower non-specific binding than that of the parental L3G7C clone. To determine whether germline library represented the endogenous immune status, specific scFv clones for nucleocapsid (N) protein of SARS-associated coronavirus (SCoV) were obtained both from naïve and immunized germline scFv libraries. Both libraries yielded specific anti-N scFvs that exhibited similar binding characteristics towards recombinant N protein, except the immunized library gave a larger number of specific anti-N scFv, and clones with identical nucleotide sequences were found. In conclusion, highly diversified antibody library can be efficiently constructed using germline rearranged immunoglobulin variable genes as source of antibody repertoires and fsPCR to diversify the CDR3.     

S-type lectins occur also in invertebrates: High conservation of the carbohydrate recognition domain in the lectin genes from the marine sponge Geodia cydonium

The marine sponge Geodia cydonium contains several lectins.The main component, called lectin-1, is composed of three tofour identical subunits. The subunits of the lectins were clonedfrom a cDNA library; two clones were obtained. From the deducedaa sequence of one clone, LECT-1, a mol. wt of 15 313 Da iscalculated; this value is in good agreement with mass spectrometricanalysis of 15 453 25 Da. The sequence of another clone, LECT-2,was analysed and the aa sequence was deduced (15 433 Da). Thetwo subunits have a framework sequence of 38 conserved aa whichare characteristic for the carbohydrate-binding site of vertebrateS-type lectins. Clustering of lectin sequences of various speciesfollowing their pairwise comparison establishes a dendrogram,which reveals that the sponge lectin could be considered asthe ancestor for vertebrate S-type lectins. Geodia cydonium lectin sponges S-type lectin     

Molecular cloning of a multifunctional chicken protein acting as the prolyl 4-hydroxylase beta-subunit, protein disulphide-isomerase and a cellular thyroid-hormone-binding protein. Comparison of cDNA-deduced amino acid sequences with those in other species.

Several recent studies indicate that a single polypeptide may act as the beta-subunit of prolyl 4-hydroxylase, the enzyme protein disulphide-isomerase and a cellular thyroid-hormone-binding protein. We report here the isolation and characterization of cDNA clones encoding this multifunctional protein in the chicken. All the coding sequences were determined on the basis of nucleotide sequencing of five cDNA clones and amino acid sequencing of the N-terminal end of the chicken beta-subunit. The processed polypeptide contains 493 amino acid residues, the size of the respective mRNA being about 2.7 kb. The chicken beta-subunit cDNA sequences were 78% homologous to the previously reported human beta-subunit cDNA sequences at the nucleotide level and 85% homologous at the amino acid level. The homology of the chicken beta-subunit sequences to those reported for bovine thyroid-hormone-binding protein and rat protein disulphide-isomerase was also 85% at the amino acid level. Primary-structure comparisons between the four species indicated that the two proposed active sites of protein disulphide-isomerase, the two Trp-Cys-Gly-His-Cys-Lys sequences, are located within highly conserved regions, which are also homologous to the active sites of a number of thioredoxins. The middle of the polypeptide has an additional conserved region 100 amino acid residues in length in which the degree of homology between the four species is 94% at the amino acid level. This long conserved region may also be important for some of the multiple functions of the protein. The four extreme C-terminal amino acids of the polypeptide in all four species are Lys-Asp-Glu-Leu, a sequence that has been suggested to function as a signal for the retention of a protein in the endoplasmic reticulum.     

Initial analysis of tandemly repetitive sequences in the genome of Zhikong scallop (Chlamys farreri Jones et Preston).

Tandemly repetitive sequences are widespread in all eukaryotic genomes, but data on tandem repeats are limited in Zhikong scallop (Chlamys farreri). In the present study, paired-end sequencing of 2016 individual fosmid clones resulted in 3646 sequences. A total of 2,286,986 bp of genomic sequences were generated, representing approximately 1.84 per thousand of the Zhikong scallop genome. Using tandem repeats finder (TRF) software, a total of 2500 tandem repeats were found, including 313 satellites, 1816 minisatellites and 371 microsatellites. The cumulative length of tandem repeats was 552,558 bp, accounting for 24.16% of total length. Specifically, the length of microsatellites, minisatellites and satellites was 9425, 336,001 and 207,132 bp, accounting for 1.71, 60.81 and 37.49% of the length of tandem repeats, and 0.41, 14.69 and 9.06% of total length, respectively. The detailed information on the characteristic of all repeat units was also represented, which will provide a useful resource for physical mapping and better utilization of the existing genomic information in Zhikong scallop.     

The 10 kb Drosophila virilis 28S rDNA intervening sequence is flanked by a direct repeat of 14 base pairs of coding sequence

Most repeat units of rDNA in Drosophila virilis are interrupted in the 28S rRNA coding region by an intervening sequence about 10 kb in length; uninterrupted repeats have a length of about 11 kb. We have sequenced the coding/intervening sequence junctions and flanking regions in two independent clones of interrupted rDNA, and the corresponding 28S rRNA coding region in a clone of uninterrupted rDNA. The intervening sequence is terminated at both ends by a direct repeat of a fourteen nucleotide sequence that is present once in the corresponding region of an intact gene. This is a phenomenon associated with transposable elements in other eukaryotes and in prokaryotes, and the Drosophila rDNA intervening sequence is discussed in this context. We have compared more than 200 nucleotides of the D. virilis 28S rRNA gene with sequences of homologous regions of rDNA in Tetrahymena pigmentosa (Wild and Sommer, 1980) and Xenopus laevis (Gourse and Gerbi, 1980): There is 93% sequence homology among the diverse species, so that the rDNA region in question (about two-thirds of the way into the 28S rRNA coding sequence) has been very highly conserved in eukaryote evolution. The intervening sequence in T. pigmentosa is at a site 79 nucleotides upstream from the insertion site of the Drosophila intervening sequence.     

Functional inhibition of transitory proteins by intrabody-mediated retention in the endoplasmatic reticulum

Intrabodies are recombinantly expressed intracellular antibody fragments that can be used to specifically bind and inhibit the function of cellular proteins of interest. Intrabodies can be targeted to various cell compartments by attaching an appropriate localization peptide sequence to them. An efficient strategy with a high success rate is to anchor intrabodies in the endoplasmatic reticulum where they can inhibit transitory target proteins by binding and preventing them to reach their site of action. Intrabodies can be assembled from antibody gene fragments from various sources into dedicated expression vectors. Conventionally, antibody cDNA sequences are derived from selected hybridoma cell clones that express antibodies with the desired specificity. Alternatively, appropriate clones can be isolated by affinity selection from an antibody in vitro display library. Here an evaluation of endoplasmatic reticulum targeted intrabodies with respect to other knockdown approaches is given and the characteristics of various intrabody expression vectors are discussed. A step by step protocol is provided that was repeatedly used to construct intrabodies derived from diverse antibody isotypes producing hybridoma cell clones. The inactivation of the cell surface receptor neural cell adhesion molecule (NCAM) by a highly efficacious novel endoplasmatic reticulum-anchored intrabody is demonstrated.     

A Novel Human scFv Library with Non-Combinatorial Synthetic CDR Diversity

The present work describes the construction and validation of a human scFv library with a novel design approach to synthetic complementarity determining region (CDR) diversification. The advantage of synthetic antibody libraries includes the possibility of exerting fine control over factors like framework sequences, amino acid and codon usage, and CDR diversity. However, random combinatorial synthesis of oligonucleotides for CDR sequence diversity also produces many clones with unnatural sequences and/or undesirable modification motifs. To alleviate these issues, we designed and constructed a novel semi-synthetic human scFv library with non-combinatorial, pre-designed CDR diversity and a single native human framework each for heavy, kappa, and lambda chain variable domains. Next-generation sequencing analysis indicated that the library consists of antibody clones with highly nature-like CDR sequences and the occurrence of the post-translational modification motifs is minimized. Multiple unique clones with nanomolar affinity could be isolated from the library against a number of target antigens, validating the library design strategy. The results demonstrate that it is possible to construct a functional antibody library using low, non-combinatorial synthetic CDR diversity, and provides a new strategy for the design of antibody libraries suitable for demanding applications.     

Amplicon –Based Metagenomic Analysis of Mixed Fungal Samples Using Proton Release Amplicon Sequencing

Next generation sequencing technology has revolutionised microbiology by allowing concurrent analysis of whole microbial communities. Here we developed and verified similar methods for the analysis of fungal communities using a proton release sequencing platform with the ability to sequence reads of up to 400 bp in length at significant depth. This read length permits the sequencing of amplicons from commonly used fungal identification regions and thereby taxonomic classification. Using the 400 bp sequencing capability, we have sequenced amplicons from the ITS1, ITS2 and LSU fungal regions to a depth of approximately 700,000 raw reads per sample. Representative operational taxonomic units (OTUs) were chosen by the USEARCH algorithm, and identified taxonomically through nucleotide blast (BLASTn). Combination of this sequencing technology with the bioinformatics pipeline allowed species recognition in two controlled fungal spore populations containing members of known identity and concentration. Each species included within the two controlled populations was found to correspond to a representative OTU, and these OTUs were found to be highly accurate representations of true biological sequences. However, the absolute number of reads attributed to each OTU differed among species. The majority of species were represented by an OTU derived from all three genomic regions although in some cases, species were only represented in two of the regions due to the absence of conserved primer binding sites or due to sequence composition. It is apparent from our data that proton release sequencing technologies can deliver a qualitative assessment of the fungal members comprising a sample. The fact that some fungi cannot be amplified by specific “conserved” primer pairs confirms our recommendation that a multi-region approach be taken for other amplicon-based metagenomic studies.     

抗人AFP抗体噬菌体呈现文库的构建筛选及基因序列测定

从人ＡＦＰ免疫小鼠脾细胞ｍＲＮＡ中扩增出全套抗体Ｖ区基因并随机拼接为ＳｃＦｖ基因,构建全套ＳｃＦｖ基因噬菌体呈现文库,经两轮ｐａｎｎｉｎｇ筛选富集,一次从９４个单个重组噬菌体克隆中筛选到１４个具有ＡＦＰ结合活性的克隆,测定两个阳性克隆中ＳｃＦｖ基因的核苷酸序列,获得一个Ｖ＿Ｈ基因和两个Ｖ＿Ｋ基因,其基因序列分别与鼠ＩｇＶ＿ＨＪ５５８、Ｖ＿ｋ－ＯＸ１和Ｖ＿ｋ４／５家族同源性最高;推导出的氨基酸序列中均含有抗体Ｖ区特征性的两个恒定的半胱氨酸残基、具有明确的三个ＣＤＲ和四个ＦＲ序列,表明这三个基因均系新发现的功能性鼠抗体Ｖ区基因序列。     

从大容量噬菌体抗体库中筛选抗Acr蛋白人源单链抗体

目的:从天然的大容量噬菌体抗体库中筛选特异的抗结核分枝杆菌晶体蛋白( alpha-crystallin Acr)的人源抗体.方法:以结核分枝杆菌Acr蛋白包被免疫管,通过对噬菌体抗体库进行4轮“吸附-洗脱-扩增”的过程从大容量抗体库中筛选特异性抗结核分枝杆菌Acr蛋白的抗体,并对可变区序列进行了测序分析.将特异性的噬菌体抗体感染HB2151菌,经IPTG诱导表达,制备了抗结核分枝杆菌Acr蛋白的可溶性单链抗体;对其序列和抗原结合活性进行分析鉴定.结果:经过4轮筛选,获得了43个与结核分枝杆菌Acr蛋白结合的阳性克隆,其中29个特异结合的克隆;测序分析有26不同的可变区片段;通过可溶性单链抗体(scFv)表达筛选到14株特异性结合Acr蛋白的可溶性单链抗体克隆;经过基因测序,分析了可变区基因的亚群.成功制备了可溶性单链抗体.Westren blotting分析证实筛选的人源单链抗体能与天然蛋白结合.结论:利用单链大容量抗体库获得抗结核分枝杆菌Acr蛋白的噬菌体抗体并且成功制备抗结核分枝杆菌Acr天然蛋白的可溶性单链抗体,为今后的研究和应用奠定基础.     

Modelling of the combining sites of three anti-lysozyme monoclonal antibodies and of the complex between one of the antibodies and its epitope.

Models of the antigen combining sites of three monoclonal antibodies, which recognise different but overlapping epitopes within the 'loop' region of hen egg lysozyme (HEL), have been generated from the cDNA sequences of their Fv regions (the VL and VH domains) and the known crystal structures of immunoglobulin fragments. The alpha-carbon backbone of the structurally conserved framework region has been derived from the IgG myeloma protein NEW, and models for the hypervariable loop regions have been selected on the basis of length and maximum sequence homology. The model structures have been refined by energy minimisation. Both the size and chemical nature of the predicted combining site models correlate broadly with the epitope boundaries previously determined by affinity studies. A model of the complex formed between one antibody and the corresponding lysozyme epitope is described, and contact residues are identified for subsequent testing by oligonucleotide-directed site-specific mutagenesis.     

Heterogeneity of the ribosomal DNA episome in strains and species of Entamoeba

Ribosomal DNA sequences in several species of the genus Entamoeba are highly repeated and display restriction fragment-length polymorphism (RFLP), which has been used to identify species and differentiate strains. However, the continuous variability of the non-transcribed repeat sequences in the ribosomal episome hinders an accurate typification. Looking for more reliable markers, we used DNA probes containing conserved sequences in the ribosomal episome — coding regions for the 16S and 5.8S rRNAs and transcribed spacers flanking the rDNA sequences, and the coding region for the 3 end of the 26S rRNA — to analyse hybridization patterns from five cloned pathogenic strains of Entamoeba histoiytica, two strains of the also pathogenic Entamoeba invadens and the non-pathogenic Laredo strain of Entamoeba moshkovskii. Our results provide reliable bases for the differentation of clones, strains and species of Entamoeba and the reconstruction of E. histolytica episomes. Differences in the number and length of rDNA-containing DNA fragments, previously observed by other investigators and confirmed by us, can be better defined by the present analysis.