Evolution of immunoglobulin VH pseudogenes in chickens

In chickens, there is a single functional gene (VH1) coding for the heavy chain variable region of immunoglobulins, and immunoglobulin diversity is generated by gene conversion of the VH1 gene by many variable region pseudogenes (psi VH's) that exist on the 5' side of the VH1 gene. To understand the evolution of this unique genetic system, we conducted statistical analyses of VH1 and psi VH genes together with functional VH genes from other higher vertebrate species. The results indicate, first, that chicken VH genes are all closely related to one another and were derived relatively recently from an ancestral gene belonging to one of the three major groups of VH genes in higher vertebrates. Second, the rate of nonsynonymous substitution is slightly higher than that of synonymous substitution in the complementarity- determining regions (CDRs), which suggests that diversity-enhancing selection has operated in the CDRs even for pseudogenes. However, both the rates of synonymous and nonsynonymous substitution are higher in the CDRs than in the framework regions (FRs), apparently because of an interaction between positive selection and meiotic gene conversion in the CDRs. Third, a dot matrix analysis of the psi VH genes and genomic diversity (D) genes has indicated that the 3' end of psi VH genes is attached by D-gene-like sequences, and this region of psi VH genes has high similarity with D gene sequences. This suggests that V and D genes were fused at some point of evolutionary time and this fused element multiplied by gene duplication. Finally, two alternative hypotheses of explaining the evolution of the chicken VH gene system are presented.      

Analysis of VH gene diversity in rainbow trout (Oncorhynchus mykiss): both nonsynonymous and synonymous nucleotide changes are more frequent in CDRs than in FRs

 Forty-six immunoglobulin VH gene sequences of rainbow trout were compiled to analyze the extent of variations and the frequency of nucleotide changes in CDRs and FRs. The results show that the frequency of nonsynonymous (amino acid replacing) changes (Ka) are on average 4.9 times higher in complementarity determining regions (CDRs) than in FRs, thus contributing more diversity in CDRs. Unexpectedly, however, the frequency of synonymous (silent) changes (Ks) show the same tendency: it was 5.3 times higher in CDRs than in framework regions (FRs). The distribution of Ks/Ka ratios of each comparison shows no segregation between CDRs and FRs. The same analysis applied to five germline VH genes of Heterodontus francisci shows the same result as was found with the rainbow trout. In contrast, the results from mouse data show that, while the CDR/FR ratio for Ka is much higher (7.4), the CDR/FR ratio for Ks is only slightly higher (1.8). The distribution of Ks/Ka ratios in mouse indicates clear segregation between CDRs and FRs. This suggests that CDR germline diversity is largely generated by gene conversion in VHs of lower vertebrates such as rainbow trout or shark. This mechanism might be advantageous to lower vertebrates in generating V gene diversity faster than other mechanisms such as point mutation and selection. Received: 10 July 1996 / Revised: 9 September 1996     

Relative contributions of germline gene variation and somatic mutation to immunoglobulin diversity in the mouse

The relative contributions of germline gene variation and somatic mutation to immunoglobulin diversity were studied by comparing germline gene sequences with their rearranged counterparts for the mouse VH, V kappa, and V lambda genes. The mutation rate at the amino acid level was estimated to be 7.0% in the first and second complementarity- determining regions (CDRs) and 2.0% in the framework regions (FRs). The difference in the mutation rate at the nucleotide level between the CDRs and FRs was of the same order of magnitude as that for the amino acid level. Analysis of amino acid diversity or nucleotide diversity indicated that the contribution of somatic mutation to immunoglobulin diversity is approximately 5%. However, the contribution of somatic mutation to the number of different amino acid sequences of immunoglobulins is much larger than that estimated by the analysis of amino acid diversity, and more than 90% of the different immunoglobulins seem to be generated by somatic mutation. Examination of the pattern of nucleotide substitution has suggested that clonal selection after somatic mutation may not be as strong as generally believed.      

A statistical examination of hypervariability in complementarity-determining regions of immunoglobulins

To determine the relative importance of gene conversion followed by natural selection and of natural selection for point mutation in generating variability in immunoglobulins, the numbers of synonymous and nonsynonymous substitutions in immunoglobulin sequences of various subgroups were estimated for complementarity-determining regions (CDRs) and for framework regions (FRs). Both the number of synonymous substitutions and the number of nonsynonymous substitutions in the CDR were found to exceed the corresponding numbers in the FR. Therefore, gene conversion is likely to be an important mechanism for providing variability in the CDR of immunoglobulins. The correlation coefficients between the number of synonymous substitutions and the number of nonsynonymous substitutions and between the substitution number in the CDR and that in the FR were found to be very low. Again, gene conversion is thought to be responsible for this finding.     

Specific features of immunoglobulin VH genes of the Antarctic teleost Trematomus bernacchii

The somatic recombination of different germline-encoded gene segments constitutes a principal source of antibody diversity. In order to investigate the diversity in recombined gene segments encoding the immunoglobulin heavy chain of the Antarctic teleost Trematomus bernacchii, a VH library was constructed by 5'-RACE (rapid amplification of cDNA ends) using RNA isolated from the spleen of an individual specimen. Analysis of cDNA sequences of 45 rearranged VH/D/JH segments revealed specific features, such as: high number of repeats, up to 8 bp long, and palindromic sequences, especially in CDRs (complementary determining regions); occurrence of the RGYW consensus, known as mutational hot spot, higher than in other species. Sixty-four percent of single base substitutions was found within this motif. In addition, the usage of serine codons showed a clear bias for AGY in CDRs, particularly in CDR2, and for TCN in FRs (framework regions). In CDRs, the frequency of non-synonymous changes was higher than that of synonymous changes. Diversity generated by insertions/deletions occurred more often than in other species; inserted bases were often repeats of adjacent bases. In particular the CDR2 showed the highest length variability as compared to other species. Alignment of VH sequences indicated that also the gene conversion mechanism may contribute to generating diversity. These data indicate a CDR mutability higher than in other species and provide some insights into the hypermutational events that may also occur in teleosts.     

Construction and expression of a reshaped VH domain against human CD28 molecules

Compared with the amino acid sequence of a mouse anti-human CD28 VH domain antibody, the two most homologous sequences of human antibodies were pulled out from Genbank. One of them was used as the main template for the framework regions of the reshaped VH domain. While the original mouse antibody CDRs were inserted into the human acceptor FRs, some residues in human acceptor FRs, which were different from those of the original mouse FRs in corresponding positions, were then determined or, alternatively, mutagenized to their conservative properties in kappa classification. Based on the amino acid sequences of the designed VH domain, the nucleotide sequence was deduced by using E. coli bias codons. The sequence was split into ten 30 to 60 nucleotide fragments for synthesizing, then annealed and amplified by overlap PCR. Taq DNA polymerase was used in a buffer with high Mg2+ concentration to induce more random mutations, both in FRs and CDRs. A phage display library was constructed by cloning these PCR products. After three rounds of panning, several reshaped VH with high antigen binding activity were obtained. One of them had the same CDR amino acid sequences as that of the original mouse VH domain. Further study showed that it retained a high antigen binding affinity after being expressed in E. coli BL21 (DE3).     

Adaptive evolution of variable region genes encoding an unusual type of immunoglobulin in camelids

A typical immunoglobulin (Ig) molecule is composed of four polypeptide chains: two identical heavy (H) chains and two identical light (L) chains. This tetrameric structure is conserved in almost all jawed vertebrate species. However, it has been discovered that camels and llamas (family: Camelidae) possess a type of dimeric Ig that consists of two H chains only. These H chains do not associate with L chains, and they do not have the first constant region (CH1), which is present in the conventional Ig. In spite of these changes, the dimeric Ig maintains the normal immune function. To understand the evolution of the dimeric Ig, we studied the phylogenetic relationships of the variable region (V(H)H) genes of the dimeric Ig from Camelidae and those (V(H)) of the conventional Ig from mammals. The results showed that the V(H)H genes form a monophyletic cluster within one of the mammalian V(H) groups, group C. We examined the type of selective force in complementarity-determining regions (CDRs) and framework regions (FRs) by comparing the rate of synonymous (dS) and nonsynonymous (dN) substitutions. We found that the results obtained from V(H)H genes were similar to those from V(H) genes in that CDRs showed an excess of dN over dS (indicating positive selection), whereas the reverse was true for FRs (purifying selection). However, when the extent of positive selection or purifying selection was investigated at each codon site, three major differences between V(H)H and V(H) genes were found. That is, very different types of selective force were observed between V(H)H and V(H) genes (1) at the sites that contact the L chain in the conventional Ig, (2) at the sites that interact with the CH1 region in the conventional Ig, and (3) in the H1 loop. Our findings suggest that adaptive evolution has occurred in the functionally important sites of the V(H)H genes to maintain the normal immune function in the dimeric Ig.     

Comprehensive analysis and characterization of the TCR α chain sequences in the common marmoset

The common marmoset (Callithrix jacchus) is useful as a nonhuman primate model of human diseases. Although the marmoset model has great potential for studying autoimmune diseases and immune responses against pathogens, little information is available regarding the genes involved in adaptive immunity. Here, we identified one TCR α constant (TRAC), 46 TRAJ (joining), and 35 TRAV (variable) segments from marmoset cDNA. Marmoset TRAC, TRAJ, and TRAV shared 80%, 68–100%, and 79–98% identity with their human counterparts at the amino acid level, respectively. The amino acid sequences were less conserved in TRAC than in TCRβ chain constant (TRBC). Comparative analysis of TRAV between marmosets and humans showed that the rates of synonymous substitutions per site (d _S ) were not significantly different between the framework regions (FRs) and complementarity determining regions (CDRs), whereas the rates of nonsynonymous substitutions per site (d _N ) were significantly lower in the FRs than in CDRs. Interestingly, the d _N values of the CDRs were greater for TRBV than TRAV. These results suggested that after the divergence of Catarrhini from Platyrrhini, amino acid substitutions were decreased in the FRs by purifying selection and occurred more frequently in CDRβ than in CDRα by positive selection, probably depending on structural and functional constraints. This study provides not only useful information facilitating the investigation of adaptive immunity using the marmoset model but also new insight into the molecular evolution of the TCR heterodimer in primate species.     

Concerted evolution of the immunoglobulin VH gene family

With the aim of understanding the concerted evolution of the immunoglobulin VH multigene family, a phylogenetic tree for the DNA sequences of 16 mouse and five human germ line genes was constructed. This tree indicates that all genes in this family have undergone substantial evolutionary divergence. The most closely related genes so far identified in the mouse genome seem to have diverged about 6 million years (MY) ago, whereas the most distantly related genes diverged about 300 MY ago. This suggests that gene duplication caused by unequal crossing-over or gene conversion occurs very slowly in this gene family. The rate of occurrence of gene duplication in the VH gene family has been estimated to be 5 x 10(-7) per gene per year, which seems to be at least about 100 times lower than that for the rRNA gene family. This low rate of concerted evolution in the VH gene family helps retain intergenic genetic variability that in turn contributes to antibody diversity. Because of accumulation of destructive mutations, however, about one-third of the mouse and human VH genes seem to have become nonfunctional. Many of these pseudogenes have apparently originated recently, but some of them seem to have existed in the genome for more than 10 MY. The rate of nucleotide substitution for the complementarity-determining regions (CDRs) is as high as that of pseudogenes. This suggests that there is virtually no purifying selection operating in the CDRs and that germ line mutations are effectively used for generating antibody diversity.      

SAIVGeM: spreadsheet analysis of immunoglobulin VH gene mutations

The analysis of mutations in immunoglobulin heavy chain variable (IGHV) region genes is a tedious process when performed by hand on multiple sequences. This report describes a set of linked Microsoft Excel files that perform several common analyses on large numbers of IGHV sequences. The spreadsheet analysis of immunoglobulin VH gene mutations (SAIVGeM) package determines the distribution of mutations among each nucleotide, the nature of the mutation at both the nucleotide and amino acid level, the frequency of mutation in the A/G G C/T A/T (RGYW) hotspot motifs of both strand polarity, and the distribution of replacement and silent mutations among the complementarity determining regions (CDRs) and the framework regions (FRs) of the immunoglobulin gene as defined by either the Kabat or IMGT conventions. These parameters are summarized and graphically presented where appropriate. In addition, the SAIVGeM package analyzes those mutations that occur in third positions of redundant codons. Because any nucleotide change in these positions is inherently silent, these positions can be used to study the mutational spectra without biases from the selection of protein structure.     

抗人黑色素瘤单链抗体基因的克隆和分泌型表达

应用ＲＴ－ＰＣＲ技术从分泌抗人黑色素瘤单克隆抗体的杂交瘤细胞ＨＢ８７６０中克隆了抗体轻、重链可变区基因,然后用（Ｇｌｙ４Ｓｅｒ）３连接肽基因将ＶＨ、ＶＬ连接成ＳｃＦｖ基因,并进行了序列测定．计算机分析表明ＶＨ,ＶＬ均符合小鼠抗体可变区的特征,为功能性重排的抗体可变区基因．ＶＨ、ＶＬ、ｌｉｎｋｅｒ拼接正确．ＳｃＦｖ基因全长７２９ｂｐ,其中ＶＨ基因长３６０ｂｐ,编码１２０个氨基酸,ＶＬ基因长３２４ｂｐ,编码１０８个氨基酸．在噬菌粒表达载体ｐＣＡＮＴＡＢ５Ｅ中表达了可溶性的ＳｃＦｖ蛋白,表达产物经流式细胞仪检测可特异地与黑色素瘤细胞结合,不与肝癌、胃癌及良性黑痣细胞结合     

Mutation pattern of paired immunoglobulin heavy and light variable domains in chronic lymphocytic leukemia B cells

B-cell chronic lymphocytic leukemia (CLL) patients display leukemic clones bearing either germline or somatically mutated immunoglobulin heavy variable (IGHV ) genes. Most information on CLL immunoglobulins (Igs), such as the definition of stereotyped B-cell receptors (BCRs), was derived from germline unmutated Igs. In particular, detailed studies on the distribution and nature of mutations in paired heavy- and light-chain domains of CLL clones bearing mutated Igs are lacking. To address the somatic hyper-mutation dynamics of CLL Igs, we analyzed the mutation pattern of paired IGHV-diversity-joining (IGHV-D-J ) and immunoglobulin kappa/lambda variable-joining (IGK/LV-J ) rearrangements of 193 leukemic clones that displayed ≥ 2% mutations in at least one of the two immunoglobulin variable (IGV ) genes (IGHV and/or IGK/LV ). The relationship between the mutation frequency in IGHV and IGK/LV complementarity determining regions (CDRs) and framework regions (FRs) was evaluated by correlation analysis. Replacement (R) mutation frequency within IGK/LV chain CDRs correlated significantly with mutation frequency of paired IGHV CDRs in λ but not κ isotype CLL clones. CDRs of IGKV-J rearrangements displayed a lower percentage of R mutations than IGHVs. The frequency/pattern of mutations in kappa CLL Igs differed also from that in κ-expressing normal B cells described in the literature. Instead, the mutation frequency within the FRs of IGHV and either IGKV or IGLV was correlated. Notably, the amount of diversity introduced by replaced amino acids was comparable between IGHVs and IGKVs. The data indicate a different mutation pattern between κ and λ isotype CLL clones and suggest an antigenic selection that, in κ samples, operates against CDR variation.     

β淀粉样肽单克隆抗体可变区基因的5′RACE扩增及序列分析

目的:克隆并分析抗β淀粉样肽单克隆抗体轻链与重链可变区基因。方法:从分泌抗β淀粉样肽单克隆抗体的杂交瘤细胞株A8中提取总RNA,根据恒定区序列设计基因特异性引物,通过5′RACE法扩增抗体的轻链和重链可变区基因,测定并分析可变区基因序列,并克隆入pMD18-T载体。结果:重链可变区基因序列全长450bp,编码150个氨基酸残基;轻链可变区基因序列全长429bp,编码143个氨基酸残基。在GeneBank中对氨基酸序列进行比对分析,二者均符合小鼠IgG可变区基因的特征。根据Kabat法则对A8抗体轻链和重链可变区氨基酸序列基因进行分析并确定了3个抗原互补决定区(CDR)、4个框架区(FR)和信号肽。结论:通过5′RACE法得到了抗β淀粉样肽单克隆抗体轻链与重链可变区基因,为进一步研究抗体三维结构,以及对该抗体进行人源化改造奠定了基础。     

一种构建改形单域抗体的方法

为验证一种构建改形单域抗体的实用新方法,与以往方法不同的是,该方法不需要对抗体进行空间结构模拟,以确定人源抗体的FRs接受序列及在人源FRs接受序列中哪些氨基酸残基需要突变,并且该方法将抗体的改形与亲和力成熟于同一过程完成,利用该方法构建了改形抗CD28重链单域抗体,根据一种鼠源抗CD28重链单域抗体的氨基酸序列,于GenBank中查得两条与之最同源的人源抗体序列,利用其中一条的FRs作为改形抗体的主框架进行改形构建,将鼠源抗体的CDR区插入到人源FR区后,对人源FR区的一些氨基酸残基进行替换突变,替换的氨基酸残基数及替换原则主要是根据对查到的人源抗体序列,鼠源抗体序列,以及这些序列与Kabat分类中的种属序列进行的比较,为了增加改形抗体基因的多样性,对要被替换的氨基酸残基在基因合成中采用简并的方式,使要被替换的氨基酸残基和替换的氨基酸残基都有机会出现,二者出现的几率各为50%,同时,在将大小不同的合成核苷酸片段采用重叠PCR扩增以获得完整改形抗体基因时,采用高Mg^2 浓度下和使用TaqDNA聚合酶,以进一步随机引入突变,利用重叠PCR产物构建了一个噬菌体抗体库,经过3轮淘选后,获得了几个具有较高免疫学活性的改形抗体,对其中的两个抗体进行了进一步研究,将两个抗体的基因在大肠杆菌BL21（DE3）中表达,复性后的表达蛋白仍具有较高的免疫学活性,结果表明该方法是有效可行的。     

抗炭疽保护性抗原单克隆抗体可变区基因的克隆及序列分析

目的:克隆并分析抗重组炭疽芽孢杆菌保护性抗原(rPA)单克隆抗体4B2、5E1和2A8轻、重链可变区基因。方法:从分泌抗rPA特异性单抗的杂交瘤细胞株4B2、5E1和2A8中提取总RNA;利用RT-PCR技术克隆抗rPA单抗4B2、5E1和2A8的VL、VH基因,并将其连入pMD18-T载体中进行测序分析。结果:4B2VL基因全长378bp,VH基因全长414bp;5E1VL基因全长420bp,VH基因全长426bp;2A8VL基因全长384bp,VH基因全长342bp。通过分析,这些基因均符合小鼠IgGV区基因的特征,含有4个框架区、3个抗原互补决定区,而且抗体特征性的半胱氨酸残基的数量和位置也正确。结论:克隆了抗rPA单抗4B2、5E1和2A8的VL、VH基因,为下一步构建多种形式的基因工程抗体奠定了良好的基础。     

A fully synthetic human Fab antibody library based on fixed VH/VL framework pairings with favorable biophysical properties

This report describes the design, generation and testing of Ylanthia, a fully synthetic human Fab antibody library with 1.3E+11 clones. Ylanthia comprises 36 fixed immunoglobulin (Ig) variable heavy (VH)/variable light (VL) chain pairs, which cover a broad range of canonical complementarity-determining region (CDR) structures. The variable Ig heavy and Ig light (VH/VL) chain pairs were selected for biophysical characteristics favorable to manufacturing and development. The selection process included multiple parameters, e.g., assessment of protein expression yield, thermal stability and aggregation propensity in fragment antigen binding (Fab) and IgG1 formats, and relative Fab display rate on phage. The framework regions are fixed and the diversified CDRs were designed based on a systematic analysis of a large set of rearranged human antibody sequences. Care was taken to minimize the occurrence of potential posttranslational modification sites within the CDRs. Phage selection was performed against various antigens and unique antibodies with excellent biophysical properties were isolated. Our results confirm that quality can be built into an antibody library by prudent selection of unmodified, fully human VH/VL pairs as scaffolds.     

Enhanced Evolvability in Immunoglobulin V Genes Under Somatic Hypermutation

Darwinian theory requires that mutations be produced in a nonanticipatory manner; it is nonetheless consistent to suggest that mutations that have repeatedly led to nonviable phenotypes would be introduced less frequently than others—if under appropriate genetic control. Immunoglobulins produced during infection acquire point mutations that are subsequently selected for improved binding to the eliciting antigen. We and others have speculated that an enhancement of mutability in the complementarity-determining regions (CDR; where mutations have a greater chance of being advantageous) and/or decrement of mutability in the framework regions (FR; where mutations are more likely to be lethal) may be accomplished by differential codon usage in concert with the known sequence specificity of the hypermutation mechanism. We have examined 115 nonproductively rearranged human Ig sequences. The mutation patterns in these unexpressed genes are unselected and therefore directly reflect inherent mutation biases. Using a χ² test, we have shown that the number of mutations in the CDRs is significantly higher than the number of mutations found in the FRs, providing direct evidence for the hypothesis that mutations are preferentially targeted into the CDRs.     

Limited diversity of the immunoglobulin heavy chain variable domain of the emerald rockcod Trematomus bernacchii

To investigate the diversity of the immunoglobulin heavy chain variable domain of the cold adapted teleost Trematomus bernacchii, 45 cDNA clones, containing complete or partial sequences of rearranged VH/D/JH segments, were analysed. Clones were isolated from a spleen library constructed by 5' RACE or from an expression library previously constructed and immunoscreened with rabbit anti- T. bernacchii Ig heavy chain antibodies. VH sequences shared, on average, 79.9% nucleotide identity and defined only two gene families referred to as Trbe VH I and Trbe VH II, the latter comprising 89% of the VH sequences analysed in this study. A Southern blot analysis, performed with family specific probes, revealed that there are at least 25 genomic VH genes. A phylogenetic tree showed that Trbe VH I clustered with VH genes belonging to group D and Trbe VH II with those of group C. Four putative distinct D segments were found to contribute to the diversity of CDR3, which showed a high glycine content. The Shannon analysis revealed that FRs are very highly conserved. Of CDRs, CDR2 exhibits a mean entropy value higher than CDR1, contributing to variability in a significant manner. Moreover, eight distinct JH segments were identified. These findings provide several clues suggesting a limited diversity of the VH genes in the Antarctic teleost T. bernacchii.     

Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila.

Evolutionary analysis of codon bias in Drosophila indicates that synonymous mutations are not neutral, but rather are subject to weak selection at the translation level. Here we show that the effectiveness of natural selection on synonymous sites is strongly correlated with the rate of recombination, in accord with the nearly neutral hypothesis. This correlation, however, is apparent only in genes encoding short proteins. Long coding regions have both a lower codon bias and higher synonymous substitution rates, suggesting that they are affected less efficiently by selection. Therefore, both the length of the coding region and the recombination rate modulate codon bias. In addition, the data indicate that selection coefficients for synonymous mutations must vary by a minimum of one or two orders of magnitude. Two hypotheses are proposed to explain the relationship among the coding region length, the codon bias, and the synonymous divergence and polymorphism levels across the range of recombination rates in Drosophila. The first hypothesis is that selection coefficients on synonymous mutations are inversely related to the total length of the coding region. The second hypothesis proposes that interference among synonymous mutations reduces the efficacy of selection on these mutations. We investigated this second hypothesis by carrying out forward simulations of weakly selected mutations in model populations. These simulations show that even with realistic recombination rates, this interference, which we call the "small-scale" Hill-Robertson effect, can have a moderately strong influence on codon bias.     

抗人CD3改形单链抗体的构建、表达及活性测定

ＣＤ３单抗通过多种途径有效地同体的免疫状态,在临床应用中具有极大的潜力。为克服鼠源单抗用于临床的局限性,拟采用抗体工程技术研制抗人ＣＤ３改形单链抗体。首先,将鼠源ＣＤ３单抗ＯＫＴ３轻重链ＣＤＲｓ分别移植到人源抗体ＬＳ１轻链和Ｎｄ重链的框架中,经计算机模拟其空间构象,进行残基替换,确定ＣＤ３改形ＶＬ、ＶＨ氨基酸序列,化学合成改形ＶＬ、ＶＨ基因,将其分别插入至载体ｐＲＯＨ８０中,构建成抗人ＣＤ３改形单