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.      

A hyperconversion mechanism generates the chicken light chain preimmune repertoire

The chicken immunoglobulin light chain repertoire has been shown to be entirely derived from a single V lambda 1-J rearranged combination. The complete coding information of the lambda locus was determined: it comprises 25 V-hybridizing elements, all of which are pseudogenes, clustered in both orientations within 19 kb of DNA, starting 2.4 kb upstream of the V lambda 1 gene. Sequences of somatically rearranged V lambda 1 genes from embryonic and posthatching bursal cells show that diversification of light chain sequences occurs during ontogeny by a segmental gene conversion mechanism which takes place at a frequency of 0.05-0.1 per cell generation between the pseudogene pool and the unique rearranged functional V gene.     

Organization of variable region segments of the human immunoglobulin heavy chain: duplication of the D5 cluster within the locus and interchromosomal translocation of variable region segments.

We have studied the organization of variable region (V) genes of the human immunoglobulin heavy chain (H) by cosmid cloning. We isolated two independent immunoglobulin D5 clusters (D5-a and D5-b) from cosmid libraries of the human genome. Restriction maps of these two regions showed that downstream 15 kb portions of the 55 kb overlap were different although upstream 40 kb portions were almost identical. Four more D segments, (DM, DXP, DA and DK) were found around the D5 segment in the conserved region of each cluster. Nucleotide sequences of the corresponding D segments from each cluster were almost identical and they encoded potentially functional D regions. Analysis using human-rodent somatic cell hybrids demonstrated that both clusters were located in the immunoglobulin heavy chain (H) locus on chromosome 14, suggesting that the D5-a and D5-b regions evolved by internal duplication within this locus. We also isolated a 60 kb DNA region carrying four VH segments, designated as VH-F region, which was located on chromosome 16. Nucleotide sequences of the four VH segments were determined. Two of them encoded potentially functional VH segments, and the other two were pseudogenes. Some more VH segments were found to be located outside chromosome 14, by Southern blot hybridization of human-rodent hybrid cell DNAs. These results provide further evidence that the human VH locus has undergone recent reorganization.     

Simple DNA sequences in homologous flanking regions near immunoglobulin VH genes: a role in gene interaction?

Five closely related immunoglobulin VH genes (subgroup II) were compared by sequencing of several kb of DNA. In three of the genes homology greater than 75% was found along an area of 4 kb that includes the coding region. The homology in flanking regions is only slightly lower than that in the coding sequences. Two other genes, which are located on the same EcoRI fragment, show high homology to the first three genes in the coding and immediately flanking regions. In more distant flanking regions no homology is found with the first three genes. This indicates that their evolutionary history differs from that of the other three genes. A region of simple DNA sequence composed of repetitive TCC and TCA elements was found at a distance of approximately 380 bp upstream from the initiator ATG of these VH genes. This region is the site where the two sets of genes abruptly start to diverge. The structure of the simple DNA sequence in the various VH genes suggests that it may be involved in gene interaction. We propose that both simple DNA sequences and homology in flanking regions serve a function in the correction of VH genes, which seem to be rather free to diverge and drift into pseudogenes. A correction mechanism may help this gene family to maintain its two major features, multiplicity and diversity.     

Diverse organization of immunoglobulin VH gene loci in a primitive vertebrate.

The immunoglobulin (Ig) heavy chain variable (VH) gene family of Heterodontus francisci (horned shark), a phylogenetically distant vertebrate, is unique in that VH, diversity (DH), joining (JH) and constant region (CH) gene segments are linked closely, in multiple individual clusters. The V regions of 12 genomic (liver and gonad) DNA clones have been sequenced completely and three organization patterns are evident: (i) VH-D1-D2-JH-CH with unique 12/22 and 12/12 spacers in the respective D recombination signal sequences (RSSs); VH and JH segments have 23 nucleotide (nt) spacers, (ii) VHDH-JH-CH, an unusual germline configuration with joined VH and DH segments and (iii) VHDHJH-CH, with all segmental elements being joined. The latter two configurations do not appear to be pseudogenes. Another VH-D1-D2-JH-CH gene possesses a D1 segment that is flanked by RSSs with 12 nt spacers and a D2 segment with 22/12 spacers. Based on the comparison of spleen, VH+ cDNA sequences to a germline consensus, it is evident that both DH segments as well as junctional and N-type diversity account for Ig variability. In this early vertebrate, the Ig genes share unique properties with higher vertebrate T-cell receptor as well as with Ig and may reflect the structure of a common ancestral antigen binding receptor gene.     

Chicken immunoglobulin gamma-heavy chains: limited VH gene repertoire, combinatorial diversification by D gene segments and evolution of the heavy chain locus

cDNA clones encoding the variable and constant regions of chicken immunoglobulin (Ig) gamma-chains were obtained from spleen cDNA libraries. Southern blots of kidney DNA show that the variable region sequences of eight cDNA clones reveal the same set of bands corresponding to approximately 30 cross-hybridizing VH genes of one subgroup. Since the VH clones were randomly selected, it is likely that the bulk of chicken H-chains are encoded by a single VH subgroup. Nucleotide sequence determinations of two cDNA clones reveal VH, D, JH and the constant region. The VH segments are closely related to each other (83% homology) as expected for VH or the same subgroup. The JHs are 15 residues long and differ by one amino acid. The Ds differ markedly in sequence (20% homology) and size (10 and 20 residues). These findings strongly indicate multiple (at least two) D genes which by a combinatorial joining mechanism diversify the H-chains, a mechanism which is not operative in the chicken L-chain locus. The most notable among the chicken Igs is the so-called 7S IgG because its H-chain differs in many important aspects from any mammalian IgG. The sequence of the C gamma cDNA reported here resolves this issue. The chicken C gamma is 426 residues long with four CH domains (unlike mammalian C gamma which has three CH domains) and it shows 25% homology to the chicken C mu. The chicken C gamma is most related to the mammalian C epsilon in length, the presence of four CH domains and the distribution of cysteines in the CH1 and CH2 domains. We propose that the unique chicken C gamma is the ancestor of the mammalian C epsilon and C gamma subclasses, and discuss the evolution of the H-chain locus from that of chicken with presumably three genes (mu, gamma, alpha) to the mammalian loci with 8-10 H-chain genes.     

Molecular structures of human argininosuccinate synthetase pseudogenes. Evolutionary and mechanistic implications

In the human genome there is one expressed gene for argininosuccinate synthetase and 14 pseudogenes. A cDNA coding for human argininosuccinate synthetase was used to screen a human genomic library. Twenty-five unique genomic clones were isolated and extensively characterized. At least seven clones represented processed argininosuccinate synthetase pseudogenes that lost the introns in the expressed gene. Restriction mapping demonstrated that these processed pseudogenes were located in distinct regions of the human genome. Complete nucleotide sequences of two processed pseudogenes, psi AS-1 and psi AS-3, and a partial sequence of psi AS-7 were determined. Both psi AS-1 and psi AS-3 had an adenine-rich region at their 3' end and were flanked by distinct imperfect direct repeats. A comparison of these pseudogene sequences to that of the cDNA demonstrated that psi AS-1 and psi AS-3 were 93% homologous to the cDNA, whereas psi AS-7 was 89% homologous to the cDNA. Therefore, it is estimated that psi AS-1 and psi AS-3 were created 10-11 million years ago, whereas psi AS-7 arose approximately 21 million years ago. We have estimated the evolutionary rate for the expressed argininosuccinate synthetase gene based on the sequences of psi AS-1 and psi AS-3. These data indicate that the expressed argininosuccinate synthetase gene is evolving at a rate similar to that of the beta-globin gene and much faster than the alpha-tubulin gene. Furthermore, a comparison of the sequences of psi AS-1 and psi AS-3 suggests the possibility that these pseudogenes arose from a common intermediate.     

Organization and rearrangement of immunoglobulin M genes in the amphibian Xenopus.

Sequences of immunoglobulin (Ig) cDNA clones of Xenopus laevis show that at least three different VH families are expressed in association with different JH elements and different isotypes of Ig constant regions. In genomic Southern blot analysis, the VH probes for each family hybridize to a distinct set of multiple DNA fragments. In contrast, the genomic JH elements and the IgM constant region gene are localized in a single DNA fragment of approximately 15 kb. Genomic VH elements contain regulatory sequences similar to those in VH genes of shark, fish and mammals and have a leader peptide sequence that contains an intron; they encode the VH region until residue 95 and have heptamer--23-bp--nonamer motifs similar to the rearrangement signal sequences (RSS) in all other vertebrate VH elements. The six genomic JH elements so far sequenced have a nonamer--23-bp--heptamer motif at their 5' end. These RSS motifs imply the existence of DH elements. The comparison of cDNA clones that contain similar constant regions but different VH regions or JH elements suggest rearrangement events. This is shown by Southern blot analysis of erythrocyte and B cell DNA with a JH probe. Thus, the overall organization of the Xenopus Ig gene locus is similar to that of mammals but strikingly different from shark.     

Identification of rabbit genomic Ig-VH pseudogenes that could serve as donor sequences for latent allotype expression

Synthetic DNA oligomers specific for the VHa allotypes of rabbit Ig genes have been used to identify latent allotypic sequences in homozygous a1 and a2 rabbits. Two Ig VH pseudogenes containing latent a3 regions have been cloned from the genome of a homozygous a2 rabbit. Analysis of the regions associated with allotype expression indicates that these two pseudogenes contain VHa- sequences in framework region 1 (FR1) and VHa3 sequences in FR3. One gene has undergone an unusual rearrangement with a third VH gene, deleting their intervening sequences and recombining in FR3 with sequences 5' to the leader exon. Our results demonstrate the presence of latent VH sequences in the genomic DNA of normal rabbits and suggest that a mechanism such as gene conversion is responsible for expression of genetically-unexpected Ig VH genes.     

Organization, sequence and nuclease hypersensitivity of repetitive elements flanking the chicken apoVLDLII gene: extended sequence similarity to elements flanking the chicken vitellogenin gene.

Analysis of nuclease hypersensitivity of regions flanking the estrogen-dependent, chicken apoVLDLII gene has revealed an hepatic, DNaseI hypersensitive site whose sensitivity is influenced by both the developmental stage and sex of the bird. The site is located 3.0kb upstream from the gene, in a block of middle repetitive elements. Contact hybridization studies indicate that the block consists of contiguous copies of two elements with reiteration frequencies of 500-1000 and 10,000-30,000 copies per haploid genome. Sequencing of 1.8kb spanning the repeats has revealed that the higher frequency element is a member of the CR1 family. The adjacent lower frequency repeat can also be found next to another member of the CR1 family located in the 3' flanking region of the vitellogenin gene. The hypersensitive site has been mapped to one of the two most highly conserved regions of the CR1 element. This region displays homology with a silencer sequence recently identified in a CR1 element flanking the chicken lysozyme gene.     

Molecular cloning and characterization of two sets of alpha-theta genes in the rat alpha-like globin gene cluster

The highly heterogeneous rat hemoglobin system was investigated at the gene level. Two regions of the alpha-like globin gene cluster from a Wistar rat were isolated. Four lambda Dash recombinant clones carrying rat alpha-like globin genes were localized on two distinct gene regions. A region of approximately 16kb was found to contain the 5'-IIalpha1-psi theta 1-3' loci, and another of approximately 24kb the 5'-IIalpha2-psi theta2-psiI alpha3-3' loci. Both IIalpha1 and IIalpha2 are considered to be active, coding the IIalpha-globin chain. The nt sequences of IIalpha1 and IIalpha2 are identical except for six nt in the non-coding region. The psiI alpha3 locus is a truncated pseudogene. The putative promoter region of an alpha-like globin gene is joined directly to the third exon, homologous to that of Ialpha-globin cDNA. psi theta1 and psi theta2 are also pseudogenes, as evidenced by several deletions located in the protein-coding regions of these loci. The psi theta1 and psi theta2 loci exhibit extensive homology, but the restriction maps of these genes and their flanking regions differ considerably. Genomic Southern blot analyses of the total liver DNA from six rats showed the existence of three theta-globin-related genes, including psi theta1 and psi theta2. These results indicate that the two gene regions investigated are not allelic variants, but may be generated by block duplication. This is the first report of the existence of rodent theta-globin genes.     

Germ line maintenance of the pseudogene donor pool for somatic immunoglobulin gene conversion in chickens.

Somatic immunoglobulin diversity is generated in avian species by sequential gene conversion of variable (V) gene segments of the immunoglobulin heavy- and light-chain loci during B-cell development. The germ line pools of donor sequence information for somatic V-region gene conversion are found in families of V pseudogenes, located 5' of the single functional V gene of each locus. The sequence relationships among the pseudogenes (psi VL) and functional VL1 gene of the chicken light-chain alleles in three inbred strains were compared to determine the extent of diversity within the germ line pseudogene cluster. Numerous differences were observed. For example, compared with the previously reported CB allele and the G4 allele, the S3 allele contains two intact pseudogenes between psi VL16 and psi VL18. These two adjacent psi VL gene segments (psi VL17a and psi VL17b) could have given rise to the psi VL17 segment of the G4 and CB alleles by homologous recombination. The majority of other sequence polymorphisms among the psi VL alleles appear to be the result of meiotic gene conversion. The incidence of untemplated mutations within psi VL segments is significantly lower than the incidence of mutation within the pseudogene flanking regions. Together with the observations that most psi VL segments have open reading frames and lack stop codons, these data support the hypothesis that the psi VL cluster resembles a functional multigene family maintained by evolutionary selection for its functional role in generating somatic antibody diversity. Meiotic gene conversion events within the psi VL cluster serve both to introduce diversity by the exchange of short segments between family members and to prevent the accumulation of random mutations.     

Conserved regulatory elements in the promoter region of the N-CAM gene.

Genomic clones containing 5'-flanking sequences, the first exon, and the entire first intron from the chicken N-CAM gene were characterized by restriction mapping and DNA sequencing. A > 600-bp segment that includes the first exon is very G + C-rich and contains a large proportion of CpG dinucleotides, suggesting that it represents a CpG island. SP-1 and AP-1 consensus elements are present, but no TATA- or CCAAT-like elements were found within 300 bp upstream of the first exon. Comparison of the chicken promoter region sequence with similar regions of the human, rat, and mouse N-CAM genes revealed that some potential regulatory elements including a "purine box" seen in mouse and rat N-CAM genes, one of two homeodomain binding regions seen in mammalian N-CAM genes, and several potential SP-1 sites are not conserved within this region. In contrast, high CpG content, a homeodomain binding sequence, an SP-1 element, an octomer element, and an AP-1 element are conserved in all four genes. The first intron of the chicken gene is 38 kb, substantially smaller than the corresponding intron from mammalian N-CAM genes. Together with previous studies, this work completes the cloning of the chicken N-CAM gene, which contains at least 26 exons distributed over 85 kb.