A gene cluster encoding human epidermis-type lipoxygenases at chromosome 17p13.1: cloning, physical mapping, and expression

Epidermis-type lipoxygenases, a distinct subclass within the multigene family of mammalian lipoxygenases (LOX), comprise recently discovered novel isoenzymes isolated from human and mouse skin including human 15-LOX-2, human and mouse 12R-LOX, mouse 8S-LOX, and mouse e-LOX-3. We have isolated the human homologue of mouse e-LOX-3. The cDNA of 3362 bp encodes a 711-amino-acid protein displaying 89% sequence identity with the mouse protein and exhibiting the same unusual structural feature, i.e., an extra segment of 41 amino acids, which can be located beyond the N-terminal beta-barrel domain at the surface of the C-terminal catalytic domain. The gene encoding e-LOX-3, ALOXE3, was found to be part of a gene cluster of approximately 100 kb on human chromosome 17p13.1 containing in addition the 12R-LOX gene, ALOX12B, the 15-LOX-2 gene, ALOX15B, and a novel 15-LOX pseudogene, ALOX15P. ALOXE3 and ALOX12B are arranged in a head-to-tail fashion separated by 8.5 kb. The genes are split into 15 exons and 14 introns spanning 22 and 15 kb, respectively. ALOX15P was found on the opposite DNA strand directly adjacent to the 3'-untranslated region of ALOX12B. ALOX15B is located in the same orientation 25 kb downstream of ALOX12B, and is composed of 14 exons and 13 introns spanning a total of 9.7 kb of genomic sequence. RT-PCR analysis demonstrated a predominant expression of ALOXE3, ALOX12B, and ALOX15B in skin.     

Human U1 small nuclear RNA genes: extensive conservation of flanking sequences suggests cycles of gene amplification and transposition.

The DNA immediately flanking the 164-base-pair U1 RNA coding region is highly conserved among the approximately 30 human U1 genes. The U1 multigene family also contains many U1 pseudogenes (designated class I) with striking although imperfect flanking homology to the true U1 genes. Using cosmid vectors, we now have cloned, characterized, and partially sequenced three 35-kilobase (kb) regions of the human genome spanning U1 homologies. Two clones contain one true U1 gene each, and the third bears two class I pseudogenes 9 kb apart in the opposite orientation. We show by genomic blotting and by direct DNA sequence determination that the conserved sequences surrounding U1 genes are much more extensive than previously estimated: nearly perfect sequence homology between many true U1 genes extends for at least 24 kb upstream and at least 20 kb downstream from the U1 coding region. In addition, the sequences of the two new pseudogenes provide evidence that class I U1 pseudogenes are more closely related to each other than to true genes. Finally, it is demonstrated elsewhere (Lindgren et al., Mol. Cell. Biol. 5:2190-2196, 1985) that both true U1 genes and class I U1 pseudogenes map to chromosome 1, but in separate clusters located far apart on opposite sides of the centromere. Taken together, these results suggest a model for the evolution of the U1 multigene family. We speculate that the contemporary family of true U1 genes was derived from a more ancient family of U1 genes (now class I U1 pseudogenes) by gene amplification and transposition. Gene amplification provides the simplest explanation for the clustering of both U1 genes and class I pseudogenes and for the conservation of at least 44 kb of DNA flanking the U1 coding region in a large fraction of the 30 true U1 genes.     

Characterization of the aldose reductase-encoding gene family in rat.

Although the enzyme aldose reductase (AR) is implicated in the development of tissue pathology in diabetes, the exact mechanism of this involvement remains unclear. To better understand the role that expression of the aldose reductase-encoding gene (ALR) may play in diabetic complications, we have begun to analyze the gene and its regulatory regions, and we present here the sequence of four ALR genes in the rat. The putative functional gene is 14.1 kb long, has ten exons which show perfect sequence identity to the rat lens AR RNA sequence, and nine introns with classical splice-site consensus sequences. Potential regulatory elements in the 5'-flanking region of this gene include a TATA box and two CCAAT boxes. Probing rat genomic Southern blots with a fragment from the first intron indicates that there is probably only one copy of this gene in the rat genome. The other three genes are processed pseudogenes which show approx. 90% identity to the rat lens AR RNA sequence, contain no introns, and have poly(A) regions at their 3' ends. Chromosomal localization studies show the presence of ALR genes on chromosomes 3, 4 and 6 in the rat with the putative functional gene mapped on chromosome 4.     

Presence of a differentially expressed U3A RNA variant in mouse. Structure and evolutive implications

A U3 RNA variant has been identified in mouse, the abundance of which relative to the previously characterized major form (U3B) appears to vary to a large extent depending upon the cell origin. Its partial sequence analysis shows that it is clearly related to the U3A form previously described in rat. Sequence comparisons suggest that the separation of the two forms of U3 genes now found in rat and mouse represent a relatively ancient event in rodent evolution. While mouse U3B RNA is encoded by four clustered genes, the U3A variant is encoded by a unique gene. Both mouse U3 RNAs differ substantially in primary structure (more than 10% divergence). Although rodent U3 RNAs exhibit a largely similar secondary structure, a specific difference between the A and B form can nevertheless be observed.     

Structural analysis of gene loci for rat U1 small nuclear RNA.

Four phage clones which hybridize with U1 small nuclear RNA were obtained from a rat gene library. Two clones contain a presumed pseudogene. A third clone includes two gene candidates that are co-linear with the rat U1-RNA, 3.6kb apart and in the opposite orientation. The two genes are surrounded by identical sequences of 491bp upstream and 178bp downstream. The upstream sequences do not contain a TATA box, but share many block homologies with those for the human U1-RNA gene(1-3). A 101bp "identifier (ID) sequence", which was reported to be specifically expressed in rat brain (4), is inserted immediately after the shared sequence downstream of one of the genes. In the fourth clone, there are two putative pseudogenes, which have one or three nucleotide changes, 3kb apart and in the same orientation. Southern blot analysis of total rat DNA reveals about 50 U1-RNA genes/pseudogenes in the genome.     

Organization of ribosomal RNA gene repeats of the mouse.

The organization of the ribosomal RNA (rRNA) genes of the mouse was determined by Southern blot hybridization using cloned rDNA fragments as probes, which could encompass the entire spacer region between two rRNA gene regions. The rRNA genes are organized into tandem repeats of nearly uniform length of about 44 kb. The heterogeneity detected in the nontranscribed spacer appears to be caused by its sequence rather than its length difference. At least three kinds of repetitive sequences are present in the non-transcribed spacer region; two of them are located 13 kb upstream from the 5'-end of 18S RNA gene and the other located 1 to 4 kb downstream from the 3'-end of 28S RNA gene.     

Concerted evolution within a trypsin gene cluster in Drosophila.

A cluster of four trypsin genes has previously been localized to cytological position 47D-F of the Drosophila melanogaster genome. One of these genes had been sequenced, and the presence of the other three genes was identified by cross-hybridization. Here, we present the DNA sequence of the entire genomic region encoding these four trypsin genes. In addition to the four previously inferred genes, we have identified a fifth trypsin-coding sequence located within this gene cluster. This new gene shows a high degree of sequence divergence (more than 30%) from the other four genes, although it retains all of the functional motifs that are characteristic of trypsin-coding sequences. In order to trace the molecular evolution of this gene cluster, we isolated and sequenced the homologous 7-kb region from the closely related species Drosophila erecta. A comparison of the DNA sequences between the two species provides strong evidence for the concerted evolution of some members of this gene family. Two genes within the cluster are evolving in concert, while a third gene appears to be evolving independently. The remaining two genes show an intermediate pattern of evolution. We propose a simple model, involving chromosome looping and gene conversion, to explain the relatively complex patterns of molecular evolution within this gene cluster.     

Gene Conversion Drives Within Genic Sequences: Concerted Evolution of Ribosomal RNA Genes in Bacteria and Archaea

Multiple copies of a given ribosomal RNA gene family undergo concerted evolution such that sequences of all gene copies are virtually identical within a species although they diverge normally between species. In eukaryotes, gene conversion and unequal crossing over are the proposed mechanisms for concerted evolution of tandemly repeated sequences, whereas dispersed genes are homogenized by gene conversion. However, the homogenization mechanisms for multiple-copy, normally dispersed, prokaryotic rRNA genes are not well understood. Here we compared the sequences of multiple paralogous rRNA genes within a genome in 12 prokaryotic organisms that have multiple copies of the rRNA genes. Within a genome, putative sequence conversion tracts were found throughout the entire length of each individual rRNA genes and their immediate flanks. Individual conversion events convert only a short sequence tract, and the conversion partners can be any paralogous genes within the genome. Interestingly, the genic sequences undergo much slower divergence than their flanking sequences. Moreover, genomic context and operon organization do not affect rRNA gene homogenization. Thus, gene conversion underlies concerted evolution of bacterial rRNA genes, which normally occurs within genic sequences, and homogenization of flanking regions may result from co-conversion with the genic sequence. Received: 31 March 2000 / Accepted: 15 June 2000     

Conservation of gene order, structure and sequence between three closely linked genes in Drosophila melanogaster and Drosophila virilis

In Drosophila melanogaster, the apparently unrelated genes anon-66Da, RpL14, and anon-66Db (from telomere to centromere) are located on a 5547 bp genomic fragment on chromosome arm 3L at cytological position 66D8. The three genes are tightly linked, and flanked by two relatively large genes with unknown function. We have taken a comparative genomic approach to investigate the evolutionary history of the three genes. To this end we isolated a Drosophila virilis 7.3 kb genomic fragment which is homologous to a 5.5 kb genomic region of D. melanogaster. Both fragments map to Muller's element D, namely to section 66D in D. melanogaster and to section 32E in D. virilis, and harbor the genes anon-66Da, RpL14, and anon-66Db. We demonstrate that the three genes exhibit a high conservation of gene topography in general and in detail. While most introns and intergenic regions reveal sequence divergences, there are, however, a number of interspersed conserved sequence motifs. In particular, two introns of the RpL14 gene contain a short, highly conserved 60 nt long sequence located at corresponding positions. This sequence represents a novel Drosophila small nucleolar RNA, which is homologous to human U49. Whereas DNA flanking the three genes shows no significant interspecies homologies, the 3'-flanking region in D. virilis contains sequences from the transposable element Penelope. The Penelope family of transposable elements has been shown to promote chromosomal rearrangements in the D. virilis species group. The presence of Penelope sequences in the D. virilis 7.3 kb genomic fragment may be indicative for a transposon-induced event of transposition which did not yet scramble the order of the three genes but led to the breakdown of sequence identity of the flanking DNA.     

Songbird genomics: analysis of 45 kb upstream of a polymorphic Mhc class II gene in red-winged blackbirds (Agelaius phoeniceus)

Here we present the sequence of a 45 kb cosmid containing a previously characterized poly-morphic Mhc class II B gene (Agph-DAB1) from the red-winged blackbird (Agelaius phoeniceus). We compared it with a previously sequenced cosmid from this species, revealing two regions of 7.5 kb and 13.0 kb that averaged greater than 97% similarity to each another, indicating a very recent shared duplication. We found 12 retroelements, including two chicken repeat 1 (CR1) elements, constituting 6.4% of the sequence and indicating a lower frequency of retroelements than that found in mammalian genomic DNA. Agph-DAB3, a new class II B gene discovered in the cosmid, showed a low rate of polymorphism and may be functional. In addition, we found a Mhc class II B gene fragment and three genes likely to be functional (encoding activin receptor type II, a zinc finger, and a putative gamma-filamin). Phylogenetic analysis of exon 2 alleles of all three known blackbird Mhc genes indicated strong clustering of alleles by locus, implying that large amounts of interlocus gene conversion have not occurred since these genes have been diverging. Despite this, interspecific comparisons indicate that all three blackbird Mhc genes diverged from one another less than 35 million years ago and are subject to concerted evolution in the long term. Comparison of blackbird and chicken Mhc promoter regions revealed songbird promoter elements for the first time. The high gene density of this cosmid confirms similar findings for the chicken Mhc, but the segment duplications and diversity of retroelements resembles mammalian sequences.     

Comparative human-mouse-rat sequence analysis of the ICAM gene cluster on HSA 19p13.2 and a 185-kb porcine region from SSC 2q

The human intercellular adhesion molecule gene (ICAM) cluster is located in a GC-rich and gene-rich region on HSA 19p13.2. We determined the complete DNA sequence of a 185-kb porcine bacterial artificial chromosome (BAC) clone containing parts of the ICAM gene cluster. We used the porcine sequence for a detailed comparative analysis between human, pig, mouse and rat. The 185 kb of porcine sequence covered 220 kb of homologous sequence in the human genome, which adds to the growing evidence that the porcine genome is somewhat smaller than the human genome. The genomic sequences of the four species showed a high level of conserved synteny and no rearrangements in gene order were observed. During evolution, the ICAM3 gene was inactivated by mutation in the mouse and rat genome, whereas it is still present in the human and pig genome. The loss of Icam3 in rodent genomes might be relevant for rodent-specific properties of the T-cell-mediated immune response. All the other investigated genes are conserved across all four investigated sequences.     

Structure and localization of the human SULT1B1 gene: neighborhood to SULT1E1 and a SULT1D pseudogene.

The soluble sulfotransferases are involved in the elimination of xenobiotics, the activation of procarcinogens, and the regulation of hormones. They comprise a gene superfamily (SULT). The structure and chromosomal location of nine human SULT genes are known. We have characterized a further gene, SULT1B1. Its structure is similar to that of other SULT1 genes. However, the total length of its eight exons and the introns (33.6 kb) is larger than that of other human SULT1 genes (4 to 21 kb). The SULT1B1 gene sequence is part of a sequence entry in the unfinished High-Throughput Genomic Sequences (HTGS) division of GenBank. However, the order and orientation of the SULT1B1 exons are not correct in this entry. SULT1B1 is located on chromosome 4q13.1, nearly 100 kb downstream of SULT1E1 on the same strand. The intervening sequence contains a SULT-like structure showing substantial homology to the mouse SULT1D1 cDNA recently described. However, in humans this structure represents a pseudogene (SULT1D1P) because of mutated splice donors/acceptors and in-frame stop codons in the sequence corresponding to exon II. This SULT gene cluster is located on the minus strand of chromosome 4 with SULT1B1 being closest to the centromer.     

Pig whey acidic protein gene is surrounded by two ubiquitously expressed genes

A 140-kb pig DNA fragment containing the whey acidic protein (WAP) gene cloned in a bacterial artificial chromosome (BAC344H5) has been shown to contain all of the cis-elements necessary for position-independent, copy-dependent and tissue-specific expression in transgenic mice. The insert from this BAC was sequenced. This revealed the presence of two other genes with quite different expression patterns in pig tissues and in transfected HC11 mouse mammary cells. The RAMP3 gene is located 15 kb upstream of the WAP gene in reverse orientation. The CPR2 gene is located 5 kb downstream of the WAP gene in the same orientation. The same locus organization was found in the human genome. The region between RAMP3 and CPR2 in the human genome contains a WAP gene-like sequence with several points of mutation which may account for the absence of WAP from human milk.