Evolution of the human alpha-amylase multigene family through unequal, homologous, and inter- and intrachromosomal crossovers

Human amylase haplotypes differ from each other by different numbers of a long direct repeat unit of approximately 100 kb, encompassing two complete salivary amylase genes and one amylase pseudogene lacking the first three exons. The two salivary genes are part of a 75-kb-long inverted repeat. Two short sequences, hybridizing with a probe containing exons 1-3, were found in the central part of the inverted repeat. Sequencing showed that these fragments, designated r, contain exon 3 sequences. We present evidence that these r-fragments and the pseudogene most likely are remnants of the same ancestral pancreatic gene. We determined the orientation of the exon 3 sequences present in the r-fragment and show that an inversion can explain their origination. Hybridization studies, using random fragments from the intergenic region of the AMY gene cluster as probes, enabled us to detect more extended homologous regions in this cluster than were found previously on the basis of restriction maps only. Together, these results allow us to present a model for the evolution of the human amylase multigene family by a number of consecutive events involving inter- and intrachromosomal crossovers.     

Retrotransposons and the evolution of mammalian gene expression

Transposable elements, and retroviral-like elements in particular, are a rich potential source of genetic variation within a host's genome. Many mutations of endogenous genes in phylogenetically diverse organisms are due to insertion of elements that affect gene expression by altering the normal pattern of regulation. While few such associations are known to have been maintained over time, two recently elucidated examples suggest transposable elements may have a significant impact in evolution of gene expression. The first example, concerning the mouse sex-limited protein (Slp), clearly establishes that ancient retroviral enhancer sequences now confer hormonal dependence on the adjacent gene. The second example shows that within the human amylase gene family, salivary specific expression has arisen due to inserted sequences, deriving perhaps from a conjunction of two retrotransposable elements.     

Evolutionary implications of the human aldolase-A, -B, -C, and -pseudogene chromosome locations.

The aldolase genes represent an ancient gene family with tissue-specific isozymic forms expressed only in vertebrates. The chromosomal locations of the aldolase genes provide insight into their tissue-specific and developmentally regulated expression and evolution. DNA probes for the human aldolase-A and -C genes and for an aldolase pseudogene were used to quantify and map the aldolase loci in the haploid human genome. Genomic hybridization of restriction fragments determined that all the aldolase genes exist in single copy in the haploid human genome. Spot-blot analysis of sorted chromosomes mapped human aldolase A to chromosome 16, aldolase C to chromosome 17, the pseudogene to chromosome 10; it previously had mapped the aldolase-B gene to chromosome 9. All loci are unlinked and located on to two pairs of morphologically similar chromosomes, a situation consistent with tetraploidization during isozymic and vertebrate evolution. Sequence comparisons of expressed and flanking regions support this conclusion. These locations on similar chromosome pairs correctly predicted that the aldolase pseudogene arose when sequences from the aldolase-A gene were inserted into the homologous aldolase location on chromosome 10.     

Three transposed elements in the intron of a human VK immunoglobulin gene.

Two gene segments coding for the variable region of human immunoglobulin light chains of the kappa type (VK genes, ref. 2) were found to have unusual structures. The two genes which are called A6 and A22 are located in duplicated gene clusters. Their restriction maps are very similar. About 4 kb of the A22 gene region were sequenced. It turned out that the intron contains an insert with the characteristics of a transposed element. The inserted DNA of 1.2 kb length contains imperfect direct and inverted repeats at its ends; at the insertion site a duplication of five nucleotides was found. Within the inserted DNA one copy each of an Alu element and of the simple sequence motif (T-G)17 were identified. Also these two repetitive sequences are themselves flanked by short direct repeats. The major inserted DNA has no significant homology to published human nucleic acid sequences. The whole structure is interpreted best by assuming a sequential insertion of the three elements. The coding region of the VK gene itself has several mutations which by themselves would render it a pseudogene; we assume that the insertion event(s) occurred prior to the mutations. According to mapping and hybridization data A6 is very similar to A22.     

Concerted evolution of human amylase genes.

Cosmid clones containing 250 kilobases of genomic DNA from the human amylase gene cluster have been isolated. These clones contain seven distinct amylase genes which appear to comprise the complete multigene family. By sequence comparison with the cDNAs, we have identified two pancreatic amylase genes and three salivary amylase genes. Two truncated pseudogenes were also recovered. Intergenic distances of 17 to 22 kilobases separate the amylase gene copies. Within the past 10 million years, duplications, gene conversions, and unequal crossover events have resulted in a very high level of sequence similarity among human amylase gene copies. To identify sequence elements involved in tissue-specific expression and hormonal regulation, the promoter regions of the human amylase genes were sequenced and compared with those of the corresponding mouse genes. The promoters of the human and mouse pancreatic amylase genes are highly homologous between nucleotide -160 and the cap site. Two sequence elements thought to influence pancreas-specific expression of the rodent genes are present in the human genes. In contrast, similarity in the 5' flanking sequences of the salivary amylase genes is limited to several short sequence elements whose positions and orientations differ in the two species. Some of these sequence elements are also associated with other parotid-specific genes and may be involved in their tissue-specific expression. A glucocorticoid response element and a general enhancer element are closely associated in several of the amylase promoters.     

Conserved sequences in both coding and 5'' flanking regions of mammalian opal suppressor tRNA genes.

The rabbit genome encodes an opal suppressor tRNA gene. The coding region is strictly conserved between the rabbit gene and the corresponding gene in the human genome. The rabbit opal suppressor gene contains the consensus sequence in the 3' internal control region but like the human and chicken genes, the rabbit 5' internal control region contains two additional nucleotides. The 5' flanking sequences of the rabbit and the human opal suppressor genes contain extensive regions of homology. A subset of these homologies is also present 5' to the chicken opal suppressor gene. Both the rabbit and the human genomes also encode a pseudogene. That of the rabbit lacks the 3' half of the coding region. Neither pseudogene has homologous regions to the 5' flanking regions of the genes. The presence of 5' homologies flanking only the transcribed genes and not the pseudogenes suggests that these regions may be regulatory control elements specifically involved in the expression of the eukaryotic opal suppressor gene. Moreover the strict conservation of coding sequences indicates functional importance for the opal suppressor tRNA genes.     

The opcA and ΨopcB regions in Neisseria: genes, pseudogenes, deletions, insertion elements and DNA islands

Previous data have indicated that the opc gene encoding an immunogenic invasin is specific to Neisseria meningitidis (Nm) and is lacking in Neisseria gonorrhoeae (Ng). The data presented here show that Nm and Ng both contain two paralogous opc-like genes, opcA, corresponding to the former opc gene, and (psi)opcB, a pseudogene. The predicted OpcA and OpcB proteins possess transmembrane regions with conserved non-polar faces but differ extensively in four of the five surface-exposed loops. Gonococcal OpcA was expressed weakly under in vitro conditions, and it is unknown whether these bacteria can express this protein at high levels. Analysis of the sequences flanking opcA and (psi)opcB revealed a framework of conserved housekeeping genes interspersed with DNA islands. These regions also contained several pseudogenes, deletions and IS elements, attesting to considerable genome plasticity. Both opcA and (psi)opcB are located on DNA islands that have probably been imported from unrelated bacteria. A third island encodes the dcmD/dcrD R/M genes in Ng versus a small open reading frame in most strains of Nm. Rare strains of Nm were identified in which the R/M island has been imported. DNA islands in Nm and Ng seem to have been acquired by recombination via conserved flanking housekeeping genes rather than by insertion of mobile genetic elements.     

Evolution of human cytoplasmic actin gene sequences: chromosome mapping and structural characterizations of three cytoplasmic actin-like pseudogenes including one on the Y chromosome

Eight recombinant phage clones containing cytoplasmic actin-like gene sequences have been isolated from a human genomic library for structural characterization. Kpn I family repeat sequences flank six of these actin genes isolated, and Alu family repeats are scattered throughout the DNA inserts of all eight phage clones. Three of these genes are γ actin-like, and the other five are β actin-like. The complete nucleotide sequence analysis of one β and one γ actin-like genes and their flanking regions demonstrates that they both are processed pseudogenes. Using unique DNA sequences flanking these two pseudogenes as hybridization probes for human-mouse somatic cell hybrid DNAs, we have mapped the two actin pseudogenes on human chromosomes 8 and 3, respectively. We have also determined the DNA sequence of a human Y chromosome-linked, processed actin pseudogene. The different values of sequence divergence of these processed pseudogenes and their functional counterparts allow us to estimate the time of generation of the pseudogenes. The results suggest that the cDNA insertion events generating the human cytoplasmic actin-like pseudogenes have occurred at significantly different times during the evolution of primates, after their separation from other mammalian species.     

A dimorphic Alu Sb-like insertion in COL3A1 is ethnic-specific

Alu elements are a class of repetitive DNA sequences found throughout the human genome that are thought to be duplicated via an RNA intermediate in a process termed retroposition. Recently inserted Alu elements are closely related, suggesting that they are derived from a single source gene or closely related source genes. Analysis of the type III collagen gene (COL3A1) revealed a polymorphic Alu insertion in intron 8 of the gene. The Alu insertion in the COL3A1 gene had a high degree of nucleotide identity to the Sb family of Alu elements, a family of older Alu elements. The Alu sequence was less similar to the consensus sequence for the PV or Sb2 subfamilies, subfamilies of recently inserted Alu elements. These data support the observations that at least three source genes are active in the human genome, one of which is distinct from the PV and Sb2 subfamilies and predates either of these two subfamilies. Appearance of the Alu insertion in different ethnic populations suggests that the insertion may have occurred in the last 100,000 years. This Alu insert should be a useful marker for population studies and for marking COL3A1 alleles.     

A novel family of potentially mobile DNA elements encoding site-specific gene-integration functions: integrons

A family of novel mobile DNA elements is described, examples of which are found at several independent locations and encode a variety of antibiotic resistance genes. The complete elements consist of two conserved segments separated by a segment of variable length and sequence which includes inserted antibiotic resistance genes. The conserved segment located 3' to the inserted resistance genes was sequenced from Tn21 and R46, and the sequences are identical over a region of 2026 bases, which includes the sulphonamide resistance gene sull, and two further open reading frames of unknown function. The complete sequences of both the 3' and 5' conserved regions of the DNA element have been determined. A 59-base sequence element, found at the junctions of inserted DNA sequences and the conserved 3' segment, is also present at this location in the R46 sequence. A copy of one half of this 59-base element is found at the end of the sull gene, suggesting that sull, though part of the conserved region, was also originally inserted into an ancestral element by site-specific integration. Inverted or direct terminal repeats or short target site duplications, both of which are characteristics of class I and class II transposons, are not found at the outer boundaries of the elements described here. Furthermore, the conserved regions do not encode any proteins related to known transposition proteins, except the DNA integrase encoded by the 5' conserved region which is implicated in the gene insertion process. Mobilization of this element has not been observed experimentally; mobility is implied from the identification of the element in at least four independent locations, in Tn21, R46 (IncN), R388 (IncW) and Tn1696. The definitive features of these novel elements are (i) that they include site-specific integration functions (the integrase and the insertion site); (ii) that they are able to acquire various gene units and act as an expression cassette by supplying the promoter for the inserted genes. As a consequence of acquiring different inserted genes, the element exists in a variety of forms which differ in the number and nature of the inserted genes. This family of elements appears formally distinct from other known mobile DNA elements and we propose the name DNA integration elements, or integrons.     

Evolution of a multigene family that encodes the Kunitz chymotrypsin inhibitor in winged bean: a possible intermediate in the generation of a new gene with a distinct pattern of expression

Winged bean Kunitz chymotrypsin inhibitor (WCI) accumulates in an organ-specific and temporally regulated manner. The protein is encoded by a multigene family that includes at least four putative inhibitor-coding genes and three pseudogenes. The structure of the WCI genes indicates that an insertion at a 5′ proximal site occurred after duplication of the ancestral WCI gene and that several gene conversion events subsequently contributed to the evolution of this gene family. Analysis of the promoter activity of the 5′ regions of the WCI genes in transgenic tobacco showed that only the 5′ regions of the WCI-3a and WCI-3b genes, which encode the major WCI protein in winged bean, promoted the organ-specific and temporally regulated expression of a reporter gene. The 5′ region of a pseudogene, the WCI-P1 gene which contains frameshift mutations, exhibited constitutive promoter activity in tobacco, an indication that the 5′ region of the WCI-P1 gene might spontaneously have acquired new regulatory sequences during evolution. Since gene conversion is a relatively frequent event and since the homology between the WCI-P1 and WCI-3a/b genes is disrupted at a 5′ proximal site by remnants of an inserted sequence, the WCI-P1 gene appears to be a possible intermediate that could be converted into a new functional gene with a distinct pattern of expression by a single gene-conversion event. Received: 17 April 1996 / Accepted: 23 October 1996     

A human dihydrofolate reductase intronless pseudogene with an Alu repetitive sequence: multiple DNA insertions at a single chromosomal site

A dihydrofolate reductase (DHFR) pseudogene, hDHFR-psi 3 has been isolated from a human genomic DNA fragment library. Sequence analysis of this gene revealed a lack of introns and the presence of a tract of nine adenines, 90 bp downstream from the end of the coding sequence. These features suggest that hDHFR-psi 3 was derived from a processed RNA molecule that has been converted into DNA and inserted into a chromosome, analogous to the origin of three intronless human DHFR genes previously described. An interesting feature of hDHFR-psi 3 is the presence of a member of the Alu moderately repetitive DNA sequence family within the DHFR coding region. This Alu element is flanked by a 16 bp directly repeated DNA segment derived from DHFR coding sequences. The Alu element apparently has been inserted into the intronless DHFR pseudogene and thus, there have been two insertions at a single chromosomal locus. The hDHFR-psi 3 contains only the 3' half of the DHFR coding sequence. Immediately upstream from the directly repeated sequence before the Alu element is an adenine-rich tract. The DNA farther upstream is moderately repetitive and is related to neither DHFR nor Alu DNA sequence. Therefore, it seems possible that a third insertion has occurred at the same site further disrupting the hDHFR coding sequences.     

Invasion of the human albumin-alpha-fetoprotein gene family by Alu, Kpn, and two novel repetitive DNA elements

The human albumin-alpha-fetoprotein genomic domain contains 13 repetitive DNA elements randomly distributed throughout the symmetrical structures of these genes. These repeated sequences are located at different sites within the two genes. The human albumin gene contains five Alu elements within four of its 14 intervening sequences. Two of these repeats are located in intron 2, and the remaining three are located in introns 7, 8, and 11. The human alpha-fetoprotein gene contains three of these Alu elements, one in intron 4 and the remaining two in the 3'-untranslated region. In addition, the human alpha-fetoprotein gene contains a Kpn repeat and two classes of novel repeats that are absent from the human albumin gene. Six of the Alu elements within the two genes are bound by short direct repeats that harbor five base substitutions in 120 possible positions (60 bp times 2 termini). The absence of Alu repeats from analogous positions in rodents indicates that these repeats invaded the albumin-alpha-fetoprotein domain less than 85 Myr ago (the time of mammalian radiation). Furthermore, considering the conservation of terminal repeats flanking the Alu sequences of the albumin-alpha-fetoprotein domain (0.042 changes per site), we submit that the average time of Alu insertion into this gene family could have been as recently as 15-30 Myr ago.     

Cloning and molecular characterization of the chalcone synthase multigene family of Petunia hybrida

Chalcone synthase-encoding genes (chs) in Petunia hybrida comprise a multigene family. Some of the chs genes have been grouped into a subfamily, based upon their strong cross-hybridization and tight genomic linkage. From genomic libraries eight 'complete' chs genes, two chs gene 5'-fragments and two chs gene 3'-fragments have been isolated. The nucleotide sequence of six complete chs genes is presented and discussed in relation to their evolutionary origin and expression in different tissues. Each member of the family consists of two exons separated by an intron of variable size and sequence, which is located at a conserved position. The chs gene fragments represent single exons. Homology between non-linked chs genes is approx. 80% at the DNA level and restricted to protein-coding sequences. Homology between subfamily members (which are tightly linked) is higher (90-99%) and extends into untranslated regions of the gene, strengthening the view that they arose by recent gene duplications. The chsD gene contains a mutated translation stop codon, suggesting that this is an inactive (pseudo)gene. None of the other members of the gene family exhibits characteristics of a pseudogene, indicating that if gene inactivation has occurred during their evolution, it must characteristics of a pseudogene, indicating that if gene inactivation has occurred during their evolution, it must have been a recent event. Homology at the protein level between some (expressed) chs genes is surprisingly low. The possibility that these genes encode proteins with slightly different enzymatic activities is discussed.     

Gene regulation and evolution in the chorion locus of Bombyx mori. Structural and developmental characterization of four eggshell genes and their flanking DNA regions

We report on the detailed structural and developmental characterization of four chorion genes and a truncated pseudogene located within a 9.5 X 10(3) base chromosomal segment. These genes belong to the A and B multigene families and, like previously characterized moth chorion genes, are arranged in tightly linked pairs, which are divergently transcribed (A/B.L11 and A/B.L12). On the basis of their high degree of sequence divergence, the A genes define two distinct subfamilies, while the more homologous B genes represent different copies of the same gene type. The A.L11 and B.L11 introns are much longer, in each case because of a single inserted DNA segment that is missing from A.L12 or B.L12. The 2.1 X 10(3) base insertion in A.L11 is the first retrovirus-like transposable element characterized in Bombyx mori. The very short 5' flanking sequences of A/B.L11 and A/B.L12 (277 and 276 base-pairs) are distinct as shown by hybridization but both recur in additional chorion gene pairs, forming two respective classes that are expressed during distinctly different developmental periods. The divergently transcribed genes of each pair, which border the same 5' flanking sequence, are expressed co-ordinately, during the same developmental period. Detailed comparisons of the 5' flanking regions, and of the corresponding region of the Drosophila s15-1 chorion gene, revealed numerous, very short sequence elements that are shared. One such element, T-C-A-C-G-T, is also associated with all five sequenced Drosophila chorion genes. Some elements are repeated in a dyad symmetrical pattern, i.e. are associated with each of the two genes in a pair, while others, including T-C-A-C-G-T, occur only once per 5' flanking region, and, if functionally important, would presumably act bi-directionally on both genes of the pair.