Evolution of a Vκ gene family

To examine the evolution of multigene families we have selected as an example an immunoglobulin light chain variable region subgroup (V24) which has been extensively characterized in inbred mice (Mus musculus domesticus). Homologous genes have been isolated and sequenced from Mus pahari, a genetically and geographically isolated species believed to be the oldest living representative of the genus. Southern blot analysis using probes corresponding to individual genes in this subgroup reveals changes in the overall size of the family occurring at the level of individual genes but not at the level of the entire family. Nucleotide sequence analysis indicates an absence of regulatory sequences such as the CAT and TATA boxes 5 to the coding region, but a decanucleotide sequence involved in light chain expression is highly conserved. Within coding regions highly complex patterns of variation are seen which appear to reflect quite different selective pressures on various subregions of the coding sequence. Complementarity determining regions (CDR) are conserved to different extents, with the first CDR region in all family members being among the most conserved segments of the molecule. Conservation is similarly variable among framework segments, indicating complex and variable evolutionary pressures not only at the level of individual genes or their products but also at subregions within homologous molecules.     

Genetic polymorphisms of gene conversion within the duplicated human α-globin loci

Summary Of the 645 Japanese subjects studied, we have identified 10 individuals heterozygous for a chromosome with the triplicated -globin loci. The frequency of the triple -loci was 0.008 in this population, while that of the single -locus, i.e., -thalassemia 2 gene, might be lower than 0.0008. Analysis of haplotypes using particular RsaI site polymorphism in the -globin gene complex strongly suggests that the triple loci may have had multiple origins in this population.     

Human Vκ immunoglobulin gene number: Implications for the origin of antibody diversity

To assess the relative contributions of germline versus somatically mutated genes in the human immune system, we have examined the size of the kappa light-chain variable region (V_κ) gene pool. Two cloned kappa subgroup 1 (V_κ1) gene probes detected the same family of 15 to 20 crosshybridizing restriction fragments in human DNA, whereas flanking region probes detected fewer hybridizing fragments. Most of the hybridizing bands represent single-copy genes, as judged by a “gene titration” experiment. Furthermore, the number of hybridization bands is a good estimate of the haploid gene number, since we observed little polymorphism of restriction sites in the V_κ locus of eight unrelated people. A cloned V_κ3 probe hybridized to essentially the same 15–20 genes in human DNA as the V_κ1 probes. These results strongly suggest that a discrete family of 15–20 genes constitutes a large proportion of the V genes from three of the four V_κ subgroups. The small number of V_κ genes in the human genome supports the idea that somatic mutation plays a major role in the origin of antibody diversity in man.     

Altitudinal variation at duplicated β-globin genes in deer mice: effects of selection, recombination, and gene conversion

Spatially varying selection on a given polymorphism is expected to produce a localized peak in the between-population component of nucleotide diversity, and theory suggests that the chromosomal extent of elevated differentiation may be enhanced in cases where tandemly linked genes contribute to fitness variation. An intriguing example is provided by the tandemly duplicated β-globin genes of deer mice (Peromyscus maniculatus), which contribute to adaptive differentiation in blood-oxygen affinity between high- and low-altitude populations. Remarkably, the two β-globin genes segregate the same pair of functionally distinct alleles due to a history of interparalog gene conversion and alleles of the same functional type are in perfect coupling-phase linkage disequilibrium (LD). Here we report a multilocus analysis of nucleotide polymorphism and LD in highland and lowland mice with different genetic backgrounds at the β-globin genes. The analysis of haplotype structure revealed a paradoxical pattern whereby perfect LD between the two β-globin paralogs (which are separated by 16.2 kb) is maintained in spite of the fact that LD within both paralogs decays to background levels over physical distances of less than 1 kb. The survey of nucleotide polymorphism revealed that elevated levels of altitudinal differentiation at each of the β-globin genes drop away quite rapidly in the external flanking regions (upstream of the 5' paralog and downstream of the 3' paralog), but the level of differentiation remains unexpectedly high across the intergenic region. Observed patterns of diversity and haplotype structure are difficult to reconcile with expectations of a two-locus selection model with multiplicative fitness.     

Repetitive extragenic palindromic sequences, mRNA stability and gene expression: evolution by gene conversion? A review

Repetitive extragenic palindromic (REP) sequences are highly conserved inverted repeats present in up to 1000 copies on the Escherichia coli chromosome. We have shown both in vivo and in vitro that REP sequences can stabilize upstream mRNA by blocking the processive action of 3'----5' exonucleases. In a number of operons, mRNA stabilization by REP sequences plays an important role in the control of gene expression. Furthermore, differential mRNA stability mediated by the REP sequences can be responsible for differential gene expression within polycistronic operons. Despite the key role of REP sequences in mRNA stability and gene expression in a number of operons, several lines of evidence suggest that this is unlikely to be the primary reason for the exceptionally high degree of sequence conservation between REP sequences. Other possible functions for REP sequences are discussed. We propose that REP sequences may be a prokaryotic equivalent of 'selfish DNA' and that gene conversion may play a role in the evolution and maintenance of REP sequences.     

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.     

Evolution of the human sarcomeric-actin genes: Evidence for units of selection within the 3′ untranslated regions of the mRNAs

Summary The complete 3 untranslated region (3UTR) sequence of the human skeletal-actin gene has been compared with the corresponding regions of the rat and chicken skeletal-actin genes. This comparison reveals that the skeletal-actin 3UTR is composed of conserved and nonconserved segments. By using genomic Southern transfer blots and thermal stability (T_m) measurements, we found that the cardiac-actin gene 3UTR also consists of conserved and nonconserved segments. Comparison of human andXenopus laevis cardiac-actin mRNA sequences confirms the presence of a region of high similarity in the 3UTR. We conclude that subsegments of the 3UTRs of both skeletal- and cardiac-actin genes of birds and mammals are under considerable selective pressure. This suggests that these conserved sequences may have functional roles in actin-gene expression or regulation, and that these roles might be different for each actin isoform.     

Evolution of the variable gene segments and recombination signal sequences of the human T-cell receptor α/δ locus

The T-cell receptor (TCR) and loci are particularly interesting because of their unique genomic structure, in that the gene segments for each locus are interspersed. The origin of this remarkable gene segment arrangement is obscure. In this report, we investigated the evolution of the TCR and variable loci and their respective recombination signal sequences (RSSs). Our phylogenetic analyses divided the and variable gene segments into two major groups each with distinguishing motifs in both the framework and complementarity determining regions (CDRs). Sequence analyses revealed that TCR variable segments share similar CDR2 sequences with immunoglobulin light chain variable segments, possibly revealing similar evolutionary histories. Maximum likelihood analysis of the region on Chromosome 14q11.2 containing the loci revealed two possible ancestral TCR / variable segments, TRDV2 and TRAV1-1/1-2, respectively. Maximum parsimony revealed different evolutionary patterns between the variable segment and RSS of the same variable gene arguing for dissimilar evolutionary origins. Two models could account for this difference: a V(D)J recombination activity involving embedded heptamer-like motifs in the germline genome, or, more plausibly, an unequal sister chromatid crossing-over. Either mechanism would have resulted in increased diversity for the adaptive immune system.     

Phylogeny of flowering plants by the chloroplast genome sequences: in search of a “lucky gene”

One of the most complicated remaining problems of molecular-phylogenetic analysis is choosing an appropriate genome region. In an ideal case, such a region should have two specific properties: (i) results of analysis using this region should be similar to the results of multigene analysis using the maximal number of regions; (ii) this region should be arranged compactly and be significantly shorter than the multigene set. The second condition is necessary to facilitate sequencing and extension of taxons under analysis, the number of which is also crucial for molecular phylogenetic analysis. Such regions have been revealed for some groups of animals and have been designated as "lucky genes". We have carried out a computational experiment on analysis of 41 complete chloroplast genomes of flowering plants aimed at searching for a "lucky gene" for reconstruction of their phylogeny. It is shown that the phylogenetic tree inferred from a combination of translated nucleotide sequences of genes encoding subunits of plastid RNA polymerase is closest to the tree constructed using all protein coding sites of the chloroplast genome. The only node for which a contradiction is observed is unstable according to the different type analyses. For all the other genes or their combinations, the coincidence is significantly worse. The RNA polymerase genes are compactly arranged in the genome and are fourfold shorter than the total length of protein coding genes used for phylogenetic analysis. The combination of all necessary features makes this group of genes main candidates for the role of "lucky gene" in studying phylogeny of flowering plants.     

The misidentification of bacterial genes as human cDNAs: was the human D-1 tumor antigen gene acquired from bacteria?

The initial analysis of the draft copy of the human genome sequence revealed the presence of several genes that were proposed to have been directly transferred from bacteria. We investigated the human D-1 antigen as a potential lateral transfer event. We report that although the human D-1 antigen seems to be an excellent candidate for lateral transfer, it is a contaminating bacterial sequence present in a human cDNA library that was included in the human genome analysis. Furthermore, several other genes present in the publicly available databases that were included in the analysis of the human genome are also likely contaminating bacterial sequences present in cDNA libraries.     

Sheep α-globin gene sequences: Implications for their concerted evolution and for the down-regulation of the 3′ genes

In sheep as in man and most other mammals, there are two -globin genes (^I and ^II), which are expressed at different levels, the upstream gene being the most efficient. In -globin gene triplication and quadruplication, this trend is confirmed, i.e., the -chain output of the downstream genes progressively decreases. In this study, we have determined the complete sequence of the cDNAs and of both the introns in a triple- haplotype in which each gene could be recognized for the presence of distinct alleles. The sequence analysis reveals that the bodies of the three -globin genes are essentially identical (99.9% homology) and moreover indicates that the down-regulation of additional -globin genes in sheep is not the effect of sequence variation from the Cap to the Poly(A) addition sites. This striking similarity among -genes is higher than that seen in other mammals and is probably sustained by particularly efficient mechanisms of gene conversion and cross-over fixation. Correspondence to: Dr. M.S. Ristaldi     

Evolution of Mhc class II B genes in Darwin's finches and their closest relatives: birth of a new gene

The 15 extant species of Darwin's finches on the Galápagos and Cocos Islands are the products of an unfinished adaptive radiation from a founder flock of birds related to the South American species Tiaris obscura. Molecular characterization of their major histocompatibility complex ( Mhc) class II B genes has revealed the existence of several related groups of sequences (presumably encoded in distinct loci) from which one (group 5) stands out because of its low divergence over extended time periods. Analysis of group 5 exon 2 and intron 2 sequences has revealed that the encoding locus apparently arose 2-3 million years ago in the Tiaris group of South and Central American Thraupini. The locus shows no evidence of inactivation, but displays a very low degree of polymorphism, both in terms of number of alleles and genetic distances between alleles. Some of the polymorphism, however, appears to be trans-specific. All the observed intergenic differences can be explained by point mutations and most of the exon 2 changes represent non-synonymous substitutions, although the rate of non-synonymous and synonymous substitutions appears to be the same. The origin of the new locus is explained by the birth-and-death model of Mhc evolution with two important extensions. First, the ancestor of the group 5 genes may have arisen without new gene duplication and second, the birth of the new group may have been brought about by a switch from balancing to directional selection. The ancestor of the group 5 genes may have been a classical class II B allele (one of many) which directional selection fixed in the ancestral population and drove into the category of nonclassical genes.