Abstract
The VH replacement process is a RAG-mediated secondary recombination in which the variable region of a rearranged VHDJH is replaced by a different germline VH gene. In almost all human and mouse VH genes, two sequence features appear to be crucial for VH replacement. First, an embedded heptamer, which is located near the 3′ end of the rearranged VH gene, serves as a cryptic recombination signal sequence (3′cRSS) for the VH replacement process. Second, a short stretch of nucleotides located downstream of the 3′cRSS serve as a footprint of the original VH region, frequently encoding charged amino acids. In this review, we show that both of these two features are conserved in the VH genes of all jawed vertebrates, which suggests that the VH replacement process may be a conserved mechanism.
Introduction
In vertebrates, the largely diversified repertoire in the variable domain of immunoglobulins is initially generated through V(D)J recombination in developing B lymphocytes. V(D)J recombination is a somatic, cell type-specific process in which the separated germline variable (V), diversity (D) (only for immunoglobulin heavy chain, IgH), and joining (J) gene segments are assembled to form an exon that encodes the functional variable domain. The V(D)J recombination process is also a site-specific recombination that depends on the recognition, binding, and cleavage of recombination signal sequences (RSSs) by the enzymes encoded by the recombination-activating genes 1 and 2 (RAG-1 and RAG-2). RSSs flank all potentially functional V, D, and J gene segments. The consensus RSS consists of a conversed heptamer (CACAGTG) and an A-rich nonamer (ACAAAAACC), which are separated by approximately 12- or 23-bp spacers. V(D)J recombination occurs efficiently between two RSSs with different spacer lengths (Jung et al., ). Because the immunoglobulin gene loci of the IgH genes in tetrapods and teleosts are organized in a “translocon” fashion, a quasi-random selection and combination of the individual V, D, and J segments from the germline repertoire generates a large diversity of the antibody specificities (Flajnik, ; Nemazee, ). A further increase in this diversity is provided by the imprecise processing of the coding region junctions, including the deletion of nucleotides in the coding end and the addition of non-templated (N) nucleotides (Lieber et al., ).
Due to the stochastic nature of the V(D)J recombination process, B cell receptors (BCRs) that recognize autoantigens can trigger B cell central tolerance via anergy, clonal deletion, and receptor editing. Whereas anergy and clonal deletion inactivate or clear the self-reactive clones, receptor editing in immature B cells allows the cells to continue their immunoglobulin gene rearrangements to alter the specificity of their BCR (Nemazee, ). The secondary rearrangement occurs at the immunoglobulin kappa and lambda (Igκ and Igλ) loci, as previously demonstrated by many studies. However, the ongoing rearrangement of the IgH gene locus was once considered impossible because the primary rearrangement deletes the entire D locus and thus no D segments with the appropriate 5′ and 3′ 12-RSS sequences remain to join with new VH and JH segments (Gay et al., ; Tiegs et al., ; Prak and Weigert, ). Originally discovered in studies of murine pre-B cell lines, a type of secondary rearrangement called VH replacement, has now been documented as a receptor editing mechanism for the IgH gene in an increasing number of studies using knock-in mice and human normal and transformed B cells (Kleinfield et al., ; Reth et al., ; Chen et al., ; Cascalho et al., ; Zhang et al., ; Koralov et al., ). The VH replacement process closely resembles the mechanism of V(D)J recombination; this process uses a normal RSS of an upstream germline VH segment and a cryptic RSS (cRSS) that is embedded close to the 3′ end of the rearranged V exon to mediate a VH-to-VHDJH recombination. In addition, VH replacement usually produces junctional diversity or leads to frame-shifts in vivo (Covey et al., ; Chen et al., ; Koralov et al., ).
Conservation of the 3′cRSS and a downstream charged amino acid-encoding nucleotide sequence in the VH genes of human and mouse
In almost all known human germline VH genes (47/51), the cRSS is composed of a heptamer (TACTGTG) in the opposite orientation to the RSS of the germline VH segments. In addition, no conserved nonamer similar to the consensus nonamer is located upstream of the heptamer (Covey et al., ; Radic and Zouali, ). Similar conserved heptamers have been identified in more than 60% of the mouse VH nucleotide sequences that are available in GenBank (Chen et al., ). Some studies suggested that the VH replacement process is a RAG-mediated recombination process because of the detection of the double-stranded DNA breaks at the cRSS and the extrachromosomal DNA circles. Zhang et al. provided further evidence that the recombinant RAG-1/RAG-2 proteins can cleave the cRSS in vitro (Covey et al., ; Usuda et al., ; Zhang et al., ). Furthermore, many additional 3′ cryptic recombination signal sequence (3′cRSS)-like motifs that only contain the most conserved trinucleotide of the heptamer, 5′CAC (or 3′GTG), in both orientations of the coding region of the VH gene have been considered to play a role in VH gene revision, which is a second receptor replacement mechanism that occurs in germinal center B cells that may have undergone clonal expansion in response to antigen stimulation (Itoh et al., ; Wilson et al., ). Some predicted cRSSs that are initiated by the CAC motifs have been found to support detectable levels of recombination in extrachromosomal recombination assays (Davila et al., ). Therefore, any heptamer that contains a CAC motif at its 5′ end may have the potential to act as a cRSS for secondary rearrangement.
During each round of VH replacement, the recipient VH may leave a short stretch of nucleotides downstream of the 3′cRSS as a footprint. The analysis of the VH replacement footprints (the residual 3′ sequences of the replaced VH at the V-D junctions) in natural human IgH sequences by Zhang et al. indicated that the footprints frequently contribute charged amino acids to the IgH CDR3 region, regardless of the reading frame. In addition, 80% of the amino acids encoded by the 3′ end of human VH genes in all three reading frames are highly charged (Zhang et al., ). In the mouse, the arginine (Arg)-encoding AGA codon was also found at the 3′ end of most VH genes (Koralov et al., ). Previous studies have indicated that somatic mutations to Arg are common in the majority of high-affinity anti-dsDNA antibodies generated in autoimmune mice (Radic et al., ). Because the germline D genes and the normal VH-D and D-JH junctions of the IgH gene in the human and mouse rarely encode charged amino acids, the antibodies that contain VH replacement footprints may have a tendency to become autoreactive (Zhang et al., ). In addition, antibodies containing an Arg-rich CDR3 are negatively selected in a mouse strain in which the IgH repertoire is generated by VH replacement, although the level of anti-DNA antibodies in the sera of these mutant mice is still elevated (Koralov et al., ). A similar observation was recently made in humans. In systematic lupus erythematosus (SLE) patients, the frequency of VH replacement is significantly higher than in healthy individuals, and more than half of the autoreactive antibodies are encoded by VH replacement products with CDR3 regions that are rich in charged amino acids (Fan, ).
The cRSS near the 3′ end of VH genes and the charged amino acid-encoding nucleotide sequence following the 3′cRSS are conserved in both human and mouse. However, the conservation of these two features is not comprehensive to all six groups of jawed vertebrates (cartilaginous fishes, teleosts, amphibians, reptiles, birds, and mammals). Because the genomic organization of the VH genes in cartilaginous fishes and birds does not provide an advantageous condition for VH replacement (McCormack et al., ; Dooley and Flajnik, ), we will present a detailed analysis of the VH genes in the other four classes of jawed vertebrates, including six mammals (mouse, Norway rat, guinea pig, rabbit, African elephant, and gray short-tailed opossum), two reptiles (painted turtle and anole lizard), one amphibian (western clawed frog), and three teleosts (zebrafish, Atlantic salmon, and channel catfish), to determine whether these two features have been conserved throughout the evolution of jawed vertebrates.
Conservation of the 3′cRSS in the functional VH genes of different vertebrates
In our analysis, the functional germline VH sequences are available from the IMGT database (www.imgt.org) (for mouse and Norway rat), Ensembl genome database (www.ensembl.org) (for western clawed frog, painted turtle, and anole lizard) and other references (Ros et al., ; Danilova et al., ; Bengten et al., ; Wang et al., ; Yasuike et al., ; Guo et al., , ).
Regarding only the canonical heptamer (TACTGTG) of the 3′cRSS, the percentage of VH genes with an embedded 3′cRSS varies widely among the listed species, from zero in the rabbit to 90.5% in the opossum (Table 1). If only those heptamers that contain the critical 3′ GTG (NNNNGTG) are considered to be functional 3′cRSSs, the 3′cRSS is present in more than 65% of the VH genes in all analyzed species except the channel catfish; in addition, the percentage of VH genes with this sequence is higher than 85% in most mammals (except the Norway rat), western clawed frog, zebrafish, and Atlantic salmon (Table 1). The first two nucleotides (GT) of the GTG motif of the 3′cRSS arise from the TGT codon for cysteine 104 (Cys104, IMGT numbering) and the third nucleotide (G) belonging to the following codon. More than 50% (67/122) of the heptamers that do not contain the GTG motif retain the GT nucleotides. In the majority of the germline VH genes from the most species analyzed below, the amino acid that follows Cys104 is alanine (Ala), which is encoded by a GCN codon (mouse 88/107, Norway rat 79/117, guinea pig 74/89, rabbit 12/12, African elephant 42/48, opossum 19/21, painted turtle 60/68, anole lizard 59/71, western clawed frog 31/38, zebrafish 32/33, and Atlantic salmon 43/50). Fanning et al. speculated that the 3′cRSS reflects the conservation of Cys104, which is critical for the structure of the H chain (Fanning et al., ). Our analysis supports this hypothesis, but the preference of the TGT codon for Cys104 and that of the following GCN codon for Ala is also important.
Table 1
| Species | VH family | Type of the heptamer at the 3′ end of VH genesa | References | ||
|---|---|---|---|---|---|
| TACTGTG (CACTGTG) (%) | NNNNGTG (%) | Others (%) | |||
| MAMMALS | |||||
| Mouseb(Mus musculus) | IGHV1 (J558) | 56.9 (29/51) | 94.1 (48/51) | 5.9 (3/51) | IMGT database |
| IGHV2 (Q52) | 100 (8/8) | 100 (8/8) | − | ||
| IGHV3 (36–60) | 100 (6/6) | 100 (6/6) | − | ||
| IGHV4 (X-24) | 100 (1/1) | 100 (1/1) | − | ||
| IGHV5 (7183) | 90 (9/10) | 90 (9/10) | 10 (1/10) | ||
| IGHV6 (J606) | − | − | 100 (5/5) | ||
| IGHV7 (S107) | 100 (3/3) | 100 (3/3) | − | ||
| IGHV8 (3609) | 100 (6/6) | 100 (6/6) | − | ||
| IGHV9 (VGAM3-8) | − | 100 (4/4) | − | ||
| IGHV10 (VH10) | 100 (2/2) | 100 (2/2) | − | ||
| IGHV11 (CP3) | − | − | 100 (2/2) | ||
| IGHV12 (CH27) | 100 (1/1) | 100 (1/1) | − | ||
| IGHV13 (3609N) | 50 (1/2) | 50 (1/2) | 50 (1/2) | ||
| IGHV14 (SM7) | 50 (2/4) | 50 (2/4) | 50 (2/4) | ||
| IGHV15 (VH15A) | 100 (1/1) | 100 (1/1) | − | ||
| IGHV16 | − | 100 (1/1) | − | ||
| Total | 64.5 (69/107) | 86.9 (93/107) | 13.1 (14/107) | ||
| Norway rat (Rattus norvegicus) | IGHV1 | 41.7 (10/24) | 87.5 (21/24) | 12.5 (3/24) | IMGT database |
| IGHV2 | 28.1 (9/32) | 59.4 (19/32) | 40.6 (13/32) | ||
| IGHV3 | 100 (4/4) | 100 (4/4) | − | ||
| IGHV4 | 50 (1/2) | 50 (1/2) | 50 (1/2) | ||
| IGHV5 | 71.4 (15/21) | 80.9 (17/21) | 19.1 (4/21) | ||
| IGHV6 | 12.5 (1/8) | 12.5 (1/8) | 87.5 (7/8) | ||
| IGHV7 | 66.7 (4/6) | 66.7 (4/6) | 33.3 (2/6) | ||
| IGHV8 | 87.5 (7/8) | 87.5 (7/8) | 12.5 (1/8) | ||
| IGHV9 | − | 75 (3/4) | 25 (1/4) | ||
| IGHV10 | − | − | 100 (2/2) | ||
| IGHV11 | 40 (2/5) | 40 (2/5) | 60 (3/5) | ||
| IGHV12 | − | − | 100 (1/1) | ||
| Total | 45.3 (53/117) | 67.5 (79/117) | 32.5 (38/117) | ||
| Guinea pig (Cavia porcellus) | VH1 | 90.4 (19/21) | 95.2 (20/21) | 4.8 (1/21) | Guo et al., |
| VH2 | 93.8 (15/16) | 100 (16/16) | − | ||
| VH3 | 78.8 (41/52) | 84.6 (44/52) | 15.4 (8/52) | ||
| Total | 84.3 (75/89) | 89.9 (80/89) | 10.1 (9/89) | ||
| Rabbit (Oryctolagus cuniculus) | IGHV1 | − | 100 (12/12) | − | Ros et al., |
| Total | − | 100 (12/12) | − | ||
| African elephant (Loxodonta Africana) | VH1 | 66.7 (2/3) | 66.7 (2/3) | 33.3 (1/3) | Guo et al., |
| VH2 | 100 (2/2) | 100 (2/2) | − | ||
| VH3 | 83.3 (5/6) | 83.3 (5/6) | 16.7 (1/6) | ||
| VH4 | 85.3 (29/34) | 88.2 (30/34) | 11.8 (4/34) | ||
| VH5 | 100 (2/2) | 100 (2/2) | − | ||
| VH7 | 100 (1/1) | 100 (1/1) | − | ||
| Total | 85.4 (41/48) | 87.5 (42/48) | 12.5 (6/48) | ||
| Gray short-tailed opossum (Monodelphis domestica) | VH1 | 94.7 (18/19) | 94.7 (18/19) | 5.3 (1/19) | Wang et al., |
| VH2 | − | 100 (1/1) | − | ||
| VH3 | 100 (1/1) | 100 (1/1) | − | ||
| Total | 90.5 (19/21) | 95.2 (20/21) | 4.8 (1/21) | ||
| REPTILES | |||||
| Painted turtlec (Chrysemys picta) | VH1 | 89.7 (26/29) | 96.6 (28/29) | 3.4 (1/29) | Ensembl database |
| VH2 | 10 (1/10) | 10 (1/10) | 90 (9/10) | ||
| VH3 | 100 (2/2) | 100 (2/2) | − | ||
| VH4 | − | − | 100 (1/1) | ||
| VH5 | 100 (3/3) | 100 (3/3) | − | ||
| VH6 | − | − | 100 (2/2) | ||
| VH7 | 50 (3/6) | 50 (3/6) | 50 (3/6) | ||
| VH8 | 100 (1/1) | 100 (1/1) | − | ||
| VH9 | 75 (3/4) | 75 (3/4) | 25 (1/4) | ||
| VH10 | 100 (9/9) | 100 (9/9) | − | ||
| VH11 | − | 100 (1/1) | − | ||
| Total | 70.6 (48/68) | 75 (51/68) | 25 (17/68) | ||
| Anole lizardd (Anolis carolinensis) | Group I | 46.7 (7/15) | 86.7 (13/15) | 13.3 (2/15) | Ensembl database |
| Group II | 35.3 (6/17) | 76.5 (13/17) | 23.5 (4/17) | ||
| Group III | 53.3 (8/15) | 80 (12/15) | 20 (3/15) | ||
| Group IV | 65 (13/20) | 100 (20/20) | − | ||
| Ac VH38265 | 100 (1/1) | 100 (1/1) | − | ||
| Ac VH338184 | 100 (1/1) | 100 (1/1) | − | ||
| Ac VH377057 | − | − | 100 (1/1) | ||
| Ac VH405628 | − | − | 100 (1/1) | ||
| Total | 50.7 (36/71) | 84.5 (60/71) | 15.5 (11/71) | ||
| AMPHIBIANS | |||||
| Western clawed froge (Xenopus tropicalis) | VH1 | 90.9 (10/11) | 100 (11/11) | − | Ensembl database |
| VH2 | 83.3 (5/6) | 83.3 (5/6) | 16.7 (1/6) | ||
| VH3 | − | 20 (1/5) | 80 (4/5) | ||
| VH4 | 50 (1/2) | 100 (2/2) | − | ||
| VH5 | 100 (3/3) | 100 (3/3) | − | ||
| VH6 | 100 (1/1) | 100 (1/1) | − | ||
| VH8 | 100 (6/6) | 100 (6/6) | − | ||
| VH9 | − | 100 (1/1) | − | ||
| VH10 | 100 (1/1) | 100 (1/1) | − | ||
| VH11 | 100 (2/2) | 100 (2/2) | − | ||
| Total | 76.3 (29/38) | 86.8 (33/38) | 13.2 (5/38) | ||
| TELEOSTS | |||||
| Zebrafish (Danio rerio) | IGHV1 | 100 (4/4) | 100 (4/4) | − | Danilova et al., |
| IGHV2 | 100 (3/3) | 100 (3/3) | − | ||
| IGHV3 | − | 100 (1/1) | − | ||
| IGHV4 | 85.7 (6/7) | 85.7 (6/7) | 14.3 (1/7) | ||
| IGHV5 | − | 66.7 (2/3) | 33.3 (1/3) | ||
| IGHV6 | 100 (1/1) | 100 (1/1) | − | ||
| IGHV7 | 100 (1/1) | 100 (1/1) | − | ||
| IGHV8 | 75 (3/4) | 100 (4/4) | − | ||
| IGHV9 | 50 (2/4) | 50 (2/4) | 50 (2/4) | ||
| IGHV10 | − | 100 (1/1) | − | ||
| IGHV11 | 100 (2/2) | 100 (2/2) | − | ||
| IGHV13 | − | 100 (1/1) | − | ||
| IGHV14 | 100 (1/1) | 100 (1/1) | − | ||
| Total | 69.7 (23/33) | 87.9 (29/33) | 12.1 (4/33) | ||
| Atlantic salmon (Salmo salar) | IGHV1 | 63.6 (7/11) | 100 (11/11) | − | Yasuike et al., |
| IGHV2 | 100 (3/3) | 100 (3/3) | − | ||
| IGHV3 | − | 100 (1/1) | − | ||
| IGHV4 | − | 80 (4/5) | 20 (1/5) | ||
| IGHV6 | 12.5 (1/8) | 100 (8/8) | − | ||
| IGHV7 | 100 (2/2) | 100 (2/2) | − | ||
| IGHV8 | 100 (10/10) | 100 (10/10) | − | ||
| IGHV9 | − | 100 (1/1) | − | ||
| IGHV10 | 100 (1/1) | 100 (1/1) | − | ||
| IGHV12 | 100 (1/1) | 100 (1/1) | − | ||
| IGHV15 | 100 (2/2) | 100 (2/2) | − | ||
| IGHV16 | 100 (4/4) | 100 (4/4) | − | ||
| IGHV17 | − | 100 (1/1) | − | ||
| Total | 62 (31/50) | 98 (49/50) | 2 (1/50) | ||
| Channel catfish (Ictalurus punctatus) | VH1 | 100 (2/2) | 100 (2/2) | − | Bengten et al., |
| VH2 | − | − | 100 (1/1) | ||
| VH3 | 20 (1/5) | 20 (1/5) | 80 (4/5) | ||
| VH5 | − | − | 100 (1/1) | ||
| VH6 | − | − | 100 (2/2) | ||
| VH7 | − | − | 100 (5/5) | ||
| VH9 | 100 (1/1) | 100 (1/1) | − | ||
| VH11 | − | 50 (1/2) | 50 (1/2) | ||
| VH12 | − | − | 100 (2/2) | ||
| VH14 | 100 (1/1) | 100 (1/1) | − | ||
| Total | 22.7 (5/22) | 27.3 (6/22) | 72.7 (16/22) | ||
Frequency of the 3′cRSS heptamer in the functional germline VH genes of 12 species.
The denominator inside the parentheses represents the total number of functional VH genes in the corresponding VH family, and the numerator represents the number of functional VH genes that contain the specified type of heptamer in the corresponding VH family. If a VH gene had more than one allele, only the 01 allele was used in the analysis of the 3′cRSS.
All of the functional VH genes were obtained from the VH repertoire of the C57BL/6 strain (imgt.org/IMGTrepertoire).
The functional germline VH genes of the painted turtle were identified from the genome scaffolds JH585110, JH585065, and JH585278 (pre.ensembl.org/Chrysemys_picta_bellii).
The functional germline VH genes of the anole lizard were identified from the genome scaffold GL343491 (ensembl.org/Anolis_carolinensis). The groups of the anole lizard VH genes are classified according to the phylogenetic study by Gambon Deza et al. (). Ac VH338184, Ac VH377057, and Ac VH405628 are three functional germline VH genes that are not included in the following phylogenetic analysis.
The functional germline VH genes of the western clawed frog were identified from the genome scaffolds GL173803, GL173909, and GL173564 (ensembl.org/Xenopus_tropicalis). Two VH genes from VH4 family are not included in the following phylogenetic analysis.
To further determine whether the maintenance of the 3′cRSS is driven by an evolutionary force, we used a phylogenetic analysis to classify the VH sequences from all 12 species into eight groups. The A, B, and C groups include the mammalian clans I, II, and III, respectively, as well as a few of the VH genes from painted turtle, anole lizard, and western clawed frog. Group D only contains the reptile VH genes. Group E consists of the VH genes from painted turtle, anole lizard, western clawed frog, and teleosts; group F contains VH genes from western clawed frog and teleosts. All VH genes from group G and H belong to the teleosts (Figure 1). We then calculated the frequency of the 3′cRSS in each VH group. As shown in Table 2, group A to F possess a high proportion of VH genes that contain the canonical heptamer motif (TACTGTG > 55%). By contrast, most VH genes from two teleost-specific VH groups, G and H, do not contain the canonical heptamer motifs (TACTGTG < 40%). However, all eight groups, regardless of the divergence time, contain high the proportion of VH genes that contain the critical 3′ GTG in the heptamer motif (NNNNGTG > 69%). Because the VH genes are subjected to divergent evolution by the birth-and-death process (Ota and Nei, ) and six VH groups (A–F) contain genes that have persisted for a long time in living species from different classes, the high frequency of the 3′cRSS in VH groups that evolved relatively later suggests that the maintenance of the 3′cRSS was positively selected during the evolution of the VH genes.
Figure 1
Table 2
| Group | VH family | Type of the 3′cRSS heptamer | |
|---|---|---|---|
| TACTGTG (CACTGTG) (%) | NNNNGTG (%) | ||
| Group A | Mus musculus IGHV1 (J558), IGHV9 (VGAM3-8), IGHV14 (SM7), IGHV15 (VH15A) | 56.9 (66/116) | 90.5 (105/116) |
| Rattus norvegicus IGHV1, IGHV9 | |||
| Loxodonta africana VH1, VH5, VH7 | |||
| Chrysemys picta VH8, VH9, VH10,VH11 | |||
| Anolis carolinensis Ac VH38265 | |||
| Xenopus tropicalis VH5, VH9, VH11 | |||
| Group B | Mus musculus IGHV2 (Q52), IGHV3 (36–60), IGHV8 (3609), IGHV12 (CH27) | 67.2 (127/189) | 78.3 (148/189) |
| Rattus norvegicus IGHV2, IGHV3, IGHV8, IGHV12 | |||
| Cavia porcellus VH1, VH2 | |||
| Loxodonta Africana VH2, VH4 | |||
| Chrysemys picta VH2, VH3, VH4, VH6 | |||
| Anolis carolinensis Group II | |||
| Xenopus tropicalis VH2, VH3, VH6, VH8 | |||
| Group C | Mus musculus IGHV4 (X-24), IGHV5 (7183), IGHV6 (J606), IGHV7 (S107), IGHV10 (VH10), IGHV11 (CP3), IGHV13 (3609N), IGHV16 | 67.3 (138/205) | 79.5 (163/205) |
| Rattus norvegicus IGHV4, IGHV5, IGHV6, IGHV7, IGHV10, IGHV11 | |||
| Cavia porcellus VH3 | |||
| Oryctolagus cuniculus IGHV1 | |||
| Loxodonta Africana VH3 | |||
| Monodelphis domestica VH1, VH2, VH3 | |||
| Chrysemys picta VH1 | |||
| Anolis carolinensis Group III | |||
| Group D | Chrysemys picta VH5 | 55.6 (10/18) | 88.9 (16/18) |
| Anolis carolinensis Group I | |||
| Group E | Chrysemys picta VH7 | 68.5 (50/73) | 87.7 (64/73) |
| Anolis carolinensis Group IV | |||
| Xenopus tropicalis VH1 | |||
| Danio rerio IGHV3, IGHV6, IGHV9, IGHV10, IGHV11, IGHV14 | |||
| Salmo salar IGHV1, IGHV8 | |||
| Ictalurus punctatus VH3 | |||
| Group F | Xenopus tropicalis VH10 | 94.1 (16/17) | 94.1 (16/17) |
| Danio rerio IGHV4 | |||
| Salmo salar IGHV2, IGHV15, IGHV16 | |||
| Group G | Danio rerio IGHV5, IGHV7, IGHV8 | 36.7 (11/30) | 76.7 (23/30) |
| Salmo salar IGHV3, IGHV6, IGHV7, IGHV10, IGHV12, IGHV17 | |||
| Ictalurus punctatus VH2, VH5, VH6, VH9, VH12, VH14 | |||
| Group H | Danio rerio IGHV1, IGHV2, IGHV13 | 39.1 (9/23) | 69.6 (16/23) |
| Salmo salar IGHV4, IGHV9 | |||
| Ictalurus punctatus VH1, VH7, VH11 | |||
Frequency of the 3′cRSS heptamer in the functional germline VH genes of the eight groups.
Conservation of the downstream charged amino acid-encoding nucleotide sequence in the functional VH genes of different vertebrates
The germline VH genes from 11 species were analyzed to determine the frequency of charged amino acids that are encoded by the nucleotide sequence following the 3′cRSS in three reading frames (Figure 2). The length of the nucleotide sequence following the 3′cRSS is usually 7 nt in seven tetrapod species and 9 nt in teleosts. Therefore, the number of amino acids that are encoded by this sequence in three reading frames is usually one, two, and two in tetrapods and two, three, and two in teleosts (Figure 2A). Due to the high percentage of A and G nucleotides in this sequence in all 11 species, the frequency of charged amino acids that are encoded by all three reading frames is higher than the random frequency (14/64, ~22%) and the frequency of the charged amino acids encoded by functional germline DH genes (Figure 2B). In addition, it is noteworthy that the frequency of charged amino acids encoded by reading frame I is greater than 60% for all 11 species (Figure 2B). Moreover, most of the charged amino acids are positively charged (Norway rat 108/115, guinea pig 65/69, rabbit 12/12, African elephant 28/30, opossum 17/17, painted turtle 65/70, anole lizard 71/81, western clawed frog 28/34, zebrafish 35/45, Atlantic salmon 45/63, and channel catfish 27/38). The reading frame I is the only correct reading frame that can ensure the encoding of a natural H chain protein; reading frame II and III might also be found in VH replacement footprints if the primary rearrangement is non-functional; thus, the VH replacement process in all 11 species should be prone to generate CDR3s that are rich in positively charged amino acids if this mechanism is conserved in jawed vertebrates.
Figure 2
Conclusion
The main conclusion from the present analysis is that both the 3′cRSS and the charged amino acid-encoding nucleotide sequence following the 3′cRSS are conserved among different classes of vertebrates, which suggests that the VH replacement may be a conserved mechanism in all jawed vertebrates. However, additional experimental evidence from species other than human and mouse are needed to support this hypothesis. The biological function of the VH replacement process in non-mammalian vertebrates is worth careful and thorough study.
Conflict of interest statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Statements
Acknowledgments
This work was supported by the National Science Fund for Distinguished Young Scholars (30725029), the National Basic Research Program of China (973 Program-2010CB945300), and the National Natural Science foundation of China (31272433).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
BengtenE.QuiniouS.HikimaJ.WaldbieserG.WarrG. W.MillerN. W.et al. (2006). Structure of the catfish IGH locus: analysis of the region including the single functional IGHM gene. Immunogenetics58, 831–844. 10.1007/s00251-006-0139-9
2
CascalhoM.WongJ.WablM. (1997). VH gene replacement in hyperselected B cells of the quasimonoclonal mouse. J. Immunol. 159, 5795–5801.
3
ChenC.NagyZ.PrakE. L.WeigertM. (1995). Immunoglobulin heavy chain gene replacement: a mechanism of receptor editing. Immunity3, 747–755.
4
CoveyL. R.FerrierP.AltF. W. (1990). VH to VHDJH rearrangement is mediated by the internal VH heptamer. Int. Immunol. 2, 579–583.
5
DanilovaN.BussmannJ.JekoschK.SteinerL. A. (2005). The immunoglobulin heavy-chain locus in zebrafish: identification and expression of a previously unknown isotype, immunoglobulin Z. Nat. Immunol. 6, 295–302. 10.1038/ni1166
6
DavilaM.LiuF.CowellL. G.LiebermanA. E.HeikampE.PatelA.et al. (2007). Multiple, conserved cryptic recombination signals in VH gene segments: detection of cleavage products only in pro B cells. J. Exp. Med. 204, 3195–3208. 10.1084/jem.20071224
7
DooleyH.FlajnikM. F. (2006). Antibody repertoire development in cartilaginous fish. Dev. Comp. Immunol. 30, 43–56. 10.1016/j.dci.2005.06.022
8
FanR. (2009). The Potential Role of VH Replacement in Editing and Generating Autoreactive Antibodies. Ph.D, The University of Alabama.
9
FanningL.BertrandF. E.SteinbergC.WuG. E. (1998). Molecular mechanisms involved in receptor editing at the Ig heavy chain locus. Int. Immunol. 10, 241–246. 10.1093/intimm/10.2.241
10
FlajnikM. F. (2002). Comparative analyses of immunoglobulin genes: surprises and portents. Nat. Rev. Immunol. 2, 688–698. 10.1038/nri889
11
Gambon DezaF.Sanchez EspinelC.Magadan MompoS. (2009). The immunoglobulin heavy chain locus in the reptile Anolis carolinensis. Mol. Immunol. 46, 1679–1687. 10.1016/j.molimm.2009.02.019
12
GayD.SaundersT.CamperS.WeigertM. (1993). Receptor editing: an approach by autoreactive B cells to escape tolerance. J. Exp. Med. 177, 999–1008. 10.1084/jem.177.4.999
13
GuoY.BaoY.MengQ.HuX.RenL.LiN.et al. (2012). Immunoglobulin genomics in the guinea pig (Cavia porcellus). PLoS ONE7:e39298. 10.1371/journal.pone.0039298
14
GuoY.BaoY.WangH.HuX.ZhaoZ.LiN.et al. (2011). A preliminary analysis of the immunoglobulin genes in the African elephant (Loxodonta africana). PLoS ONE6:e16889. 10.1371/journal.pone.0016889
15
ItohK.MeffreE.AlbesianoE.FarberA.DinesD.SteinP.et al. (2000). Immunoglobulin heavy chain variable region gene replacement as a mechanism for receptor revision in rheumatoid arthritis synovial tissue B lymphocytes. J. Exp. Med. 192, 1151–1164. 10.1084/jem.192.8.1151
16
JungD.GiallourakisC.MostoslavskyR.AltF. W. (2006). Mechanism and control of V(D)J recombination at the immunoglobulin heavy chain locus. Annu. Rev. Immunol. 24, 541–570. 10.1146/annurev.immunol.23.021704.115830
17
KleinfieldR.HardyR. R.TarlintonD.DanglJ.HerzenbergL. A.WeigertM. (1986). Recombination between an expressed immunoglobulin heavy-chain gene and a germline variable gene segment in a Ly 1+ B-cell lymphoma. Nature322, 843–846. 10.1038/322843a0
18
KoralovS. B.NovobrantsevaT. I.KonigsmannJ.EhlichA.RajewskyK. (2006). Antibody repertoires generated by VH replacement and direct VH to JH joining. Immunity25, 43–53. 10.1016/j.immuni.2006.04.016
19
LieberM. R.MaY.PannickeU.SchwarzK. (2003). Mechanism and regulation of human non-homologous DNA end-joining. Nat. Rev. Mol. Cell Biol. 4, 712–720. 10.1038/nrm1202
20
McCormackW. T.TjoelkerL. W.ThompsonC. B. (1991). Avian B-cell development: generation of an immunoglobulin repertoire by gene conversion. Annu. Rev. Immunol. 9, 219–241. 10.1146/annurev.iy.09.040191.001251
21
NemazeeD. (2006). Receptor editing in lymphocyte development and central tolerance. Nat. Rev. Immunol. 6, 728–740. 10.1038/nri1939
22
OtaT.NeiM. (1994). Divergent evolution and evolution by the birth-and-death process in the immunoglobulin VH gene family. Mol. Biol. Evol. 11, 469–482.
23
PrakE. L.WeigertM. (1995). Light chain replacement: a new model for antibody gene rearrangement. J. Exp. Med. 182, 541–548. 10.1084/jem.182.2.541
24
RadicM. Z.MackleJ.EriksonJ.MolC.AndersonW. F.WeigertM. (1993). Residues that mediate DNA binding of autoimmune antibodies. J. Immunol. 150, 4966–4977.
25
RadicM. Z.ZoualiM. (1996). Receptor editing, immune diversification, and self-tolerance. Immunity5, 505–511. 10.1016/S1074-7613(00)80266-6
26
RethM.GehrmannP.PetracE.WieseP. (1986). A novel VH to VHDJH joining mechanism in heavy-chain-negative (null) pre-B cells results in heavy-chain production. Nature322, 840–842. 10.1038/322840a0
27
RosF.PuelsJ.ReichenbergerN.Van SchootenW.BuelowR.PlatzerJ. (2004). Sequence analysis of 0.5 Mb of the rabbit germline immunoglobulin heavy chain locus. Gene330, 49–59. 10.1016/j.gene.2003.12.037
28
TiegsS. L.RussellD. M.NemazeeD. (1993). Receptor editing in self-reactive bone marrow B cells. J. Exp. Med. 177, 1009–1020. 10.1084/jem.177.4.1009
29
UsudaS.TakemoriT.MatsuokaM.ShirasawaT.YoshidaK.MoriA.et al. (1992). Immunoglobulin V gene replacement is caused by the intramolecular DNA deletion mechanism. EMBO J. 11, 611–618.
30
WangX.OlpJ. J.MillerR. D. (2009). On the genomics of immunoglobulins in the gray, short-tailed opossum Monodelphis domestica. Immunogenetics61, 581–596. 10.1007/s00251-009-0385-8
31
WeiZ.WuQ.RenL.HuX.GuoY.WarrG. W.et al. (2009). Expression of IgM, IgD, and IgY in a reptile, Anolis carolinensis. J. Immunol. 183, 3858–3864. 10.4049/jimmunol.0803251
32
WilsonP. C.WilsonK.LiuY. J.BanchereauJ.PascualV.CapraJ. D. (2000). Receptor revision of immunoglobulin heavy chain variable region genes in normal human B lymphocytes. J. Exp. Med. 191, 1881–1894. 10.1084/jem.191.11.1881
33
YasuikeM.De BoerJ.Von SchalburgK. R.CooperG. A.McKinnelL.MessmerA.et al. (2010). Evolution of duplicated IgH loci in Atlantic salmon, Salmo salar. BMC Genomics11:486. 10.1186/1471-2164-11-486
34
ZhangZ.BurrowsP. D.CooperM. D. (2004). The molecular basis and biological significance of VH replacement. Immunol. Rev. 197, 231–242. 10.1111/j.0105-2896.2004.0107.x
35
ZhangZ.ZemlinM.WangY. H.MunfusD.HuyeL. E.FindleyH. W.et al. (2003). Contribution of Vh gene replacement to the primary B cell repertoire. Immunity19, 21–31. 10.1016/S1074-7613(03)00170-5
36
ZhaoY.Pan-HammarstromQ.YuS.WertzN.ZhangX.LiN.et al. (2006). Identification of IgF, a hinge-region-containing Ig class, and IgD in Xenopus tropicalis. Proc. Natl. Acad. Sci. U.S.A. 103, 12087–12092. 10.1073/pnas.0600291103
Summary
Keywords
immunoglobulin, VH replacement, 3′cRSS, VH replacement footprint, charged amino acid, vertebrate
Citation
Sun Y, Liu Z, Li Z, Lian Z and Zhao Y (2012) Phylogenetic conservation of the 3′ cryptic recombination signal sequence (3′cRSS) in the VH genes of jawed vertebrates. Front. Immun. 3:392. doi: 10.3389/fimmu.2012.00392
Received
04 September 2012
Accepted
05 December 2012
Published
20 December 2012
Volume
3 - 2012
Edited by
Zhixin (Jason) Zhang, University of Nebraska Medical Center, USA
Reviewed by
Biao Zheng, Baylor College of Medicine and University of Texas Medical School, USA; Zhixin (Jason) Zhang, University of Nebraska Medical Center, USA
Copyright
© 2012 Sun, Liu, Li, Lian and Zhao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Yi Sun and Yaofeng Zhao, State Key Laboratory of Agrobiotechnology, College of Biological Sciences, China Agricultural University, No.2 Yuanmingyan west Road, Haidian District, Beijing 100193, China. e-mail: sunyi@cau.edu.cn; yaofengzhao@cau.edu.cn
This article was submitted to Frontiers in B Cell Biology, a specialty of Frontiers in Immunology.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.