The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura.

Nucleotide sequences of eight Est-5A and Est-5C genes corresponding to previously sequenced Est-5B genes in Drosophila pseudoobscura were determined to compare patterns of polymorphism and divergence among members of this small gene family. The three esterase genes were also sequenced from D. persimilis and D. miranda for interspecific comparisons. The data provide evidence that gene conversion between loci contributes to polymorphism and to the homogenization of the Est5 genes. For Est-5B, which encodes one of the most highly polymorphic proteins in Drosophila, 12% of the segregating amino acid variants appear to have been introduced via gene conversion from other members of the gene family. Interlocus gene conversion can also explain high sequence similarity, especially at synonymous sites, between Est-5B and Est-5A. Tests of neutrality using interspecific comparisons show that levels of polymorphism conform to neutral expectations at each Est-5 locus. However, McDonald-Kreitman tests based on intraspecific gene comparisons indicate that positive selection on amino acids has accompanied Est-5 gene duplication and divergence in D. pseudoobscura.     

Nucleotide sequence of the Xdh region in Drosophila pseudoobscura and an analysis of the evolution of synonymous codons

The nucleotide sequence of the Xdh region of Drosophila pseudoobscura is presented. The Xdh gene structure and organization are compared with the homologous region in D. melanogaster. This locus is shown to have similar organization in the two species, although an additional intron and three insertion/deletion events are described for the D. pseudoobscura coding region. The encoded proteins are predicted to have very similar charges and hydrophobic/hydrophilic domains even though 11% of the amino acids are different. A gene 5' to Xdh, putative l(3)s12, is suggested from sequence similarity between the species. Synonymous differences at the Xdh locus between the two species are analyzed using a new method described in the preceding paper by Lewontin. This analysis shows that synonymous positions within the Xdh locus are evolving at very different rates, being dependent on level of codon redundancy. A comparison of synonymous divergence between D. melanogaster and D. pseudoobscura in five additional genes reveals variation in the level of synonymous substitution.      

Conservation of alternative splicing and genomic organization of the myosin alkali light-chain (Mlc1) gene among Drosophila species

The Mlc1 gene of Drosophila melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. In larval muscle and tubular and abdominal muscles of adults, all of the six exons are included in the spliced mRNA, whereas, in the fibrillar indirect flight muscle of adult, exon 5 is excluded from the mRNA. We show that this tissue-specific pattern of alternative splicing of the Mlc1 pre-mRNA is conserved in D. simulans, D. pseudoobscura, and D. virilis. Isolation and sequencing of the Mlc1 genes from these three other Drosophila species have revealed that the overall organization of the genes is identical and that the genes have maintained a very high level of sequence identity within the coding region. Pairwise amino acid identities are 94%-99%, and there are no charge changes among the proteins. Total nucleotide divergence within the coding region of the four genes supports the accepted genealogy of these species, but the data indicate a significantly higher rate of amino acid replacement in the branch leading to D. pseudoobscura. A comparison of nucleotide substitutions in the coding portions of exon 5 and exon 6, which encode the alternative carboxyl termini of the two MLC1 isoforms, suggests that exon 5 is subject to greater evolutionary constraints than is exon 6. In addition to the coding sequences, there is significant sequence conservation within the 5' and 3' noncoding DNA and two of the introns, including one that flanks exon 5. These regions are candidates for cis- regulatory elements. Our results suggest that evolutionary constraints are acting on both the coding and noncoding sequences of the Mlc1 gene to maintain proper expression and function of the two MLC1 polypeptides.      

Molecular evolution of the Est-6 gene in Drosophila melanogaster: contrasting patterns of DNA variability in adjacent functional regions

We have investigated nucleotide polymorphism at the esterase 6 gene (Est-6) gene, including the complete coding region (1686 bp), as well as the 5'-flanking (1183 bp) and 3'-flanking (193 bp) regions of the gene, in 30 strains of Drosophila melanogaster and in one strain of Drosophila simulans. The level of silent variation is similar in the coding and in the 3'-flanking region, but smaller in the 5'-flanking region. Strong linkage disequilibrium occurs within each region; and also, although less pronounced, between the 5'-flanking region and the rest of the gene, including the 3'-flanking region. We suggest that the pattern of nucleotide polymorphism of Est-6 may be shaped by: (1) directional and balancing selection acting on the promoter and the coding region; and (2) interactions between the two regions that involve variable degrees of hitchhiking. The patterns of linkage disequilibrium, as well as the statistics Z(nS) (Genetics 146 (1997)1197) and B and Q (Genet. Res. 74 (1999) 65), may be interpreted as there being multiple targets of selection within the gene. The previously reported Est-6 allozyme latitudinal clines may be accounted for by the interaction between selective processes in the promoter and coding regions.     

DNA polymorphism in the beta-Esterase gene cluster of Drosophila melanogaster

We have analyzed nucleotide polymorphism within a 5.3-kb region encompassing the functional Est-6 gene and the psiEst-6 putative pseudogene in 28 strains of Drosophila melanogaster and one of D. simulans. Two divergent sequence types were detected, which are not perfectly associated with Est-6 allozyme variation. The level of variation (pi) is very close in the 5'-flanking region (0.0059) and Est-6 gene (0.0057), but significantly higher in the intergenic region (0.0141) and putative pseudogene (0.0122). The variation in the 3'-flanking region is intermediate (0.0083). These observations may reflect different levels of purifying selection in the different regions. Strong linkage disequilibrium occurs within the region studied, with the largest values revealed in the putative pseudogene and 3'-flanking region. Moreover, recombination is restricted within psiEst-6. Gene conversion is detected both within and (to a lesser extent) between Est-6 and psiEst-6. The data indicate that psiEst-6 exhibits some characteristics that are typical of nonfunctional genes, while other characteristics are typically attributed to functional genes; the same situation has been observed in other pseudogenes (including Drosophila). The results of structural entropy analysis demonstrate higher structural ordering in Est-6 than in psiEst-6, in accordance with expectations if psiEst-6 is indeed a pseudogene. Taking into account that the function of psiEst-6 is not known (but could exist) and following the terminology of J. Brosius and S. J. Gould, we suggest that the term "potogene" may be appropriate for psiEst-6, indicating that it is a potential gene that may have acquired some distinctive but unknown function.     

Structure of the murine mb-1 gene encoding a putative sIgM-associated molecule

Genomic DNA clones containing the B cell-specific murine mb-1 gene were isolated and a 5.6-kb BamH I fragment was characterized. It is 5629 bp long and contains five exons: an exon containing the 5' untranslated and the coding sequence of the signal peptide, an exon of 294 bp, which contains most of the extracellular sequence of the MB-1 protein, a 119-bp long exon coding mainly for the transmembrane portion, and two exons of 69 bp and 427 bp encoding the cytoplasmic domain and the 3'-untranslated region, respectively. The mb-1 gene does not contain a "TATA box" found in many eukaryotic promoters. The 5'-flanking region has sequence stretches homologous to IgVH 5'-promoter regions and a bcl 2 intron sequence. It contains the decanucleotide sequence (ATGGCAAATA) almost identical to the octamer motif of IgVH promoters. A B cell-specific DNase I-hypersensitive site was found in the 3'-flanking region indicating that this region might be involved in B cell-specific expression of mb-1. Southern blot analysis of genomic liver DNA with the cloned mb-1 cDNA suggests the existence of another mb-1-related gene segment.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Evolutionary rearrangement of the amylase genomic regions between Drosophila melanogaster and Drosophila pseudoobscura

Two Drosophila pseudoobscura genomic clones have sequence similarity to the Drosophila melanogaster amylase region that maps to the 53CD region on the D. melanogaster cytogenetic map. The two clones with similarity to amylase map to sections 73A and 78C of the D. pseudoobscura third chromosome cytogenetic map. The complete sequences of both the 73A and 78C regions were compared to the D. melanogaster genome to determine if the coding region for amylase is present in both regions and to determine the evolutionary mechanism responsible for the observed distribution of the amylase gene or genes. The D. pseudoobscura 73A and 78C linkage groups are conserved with the D. melanogaster 41E and 53CD regions, respectively. The amylase gene, however, has not maintained its conserved linkage between the two species. These data indicate that amylase has moved via a transposition event in the D. melanogaster or D. pseudoobscura lineage. The predicted genes within the 73A and 78C regions show patterns of molecular evolution in synonymous and nonsynonymous sites that are consistent with previous studies of these two species.     

The human serglycin gene. Nucleotide sequence and methylation pattern in human promyelocytic leukemia HL-60 cells and T-lymphoblast Molt-4 cells.

The complete nucleotide sequence of the 16.7-kb human gene that encodes the peptide core (serglycin) of a secretory granule proteoglycan was determined, thus representing the first proteoglycan peptide core gene to be sequenced in its entirety. The exons, intron 1, and intron 2 comprised 7, 53, and 40% of the gene, respectively. Nineteen Alu-repetitive DNA sequences were interspersed in the gene, accounting for 28% of the total nucleotides in intron 1 and 40% of the nucleotides in intron 2. The nucleotide sequence was then used in an examination of the methylation pattern of the human serglycin gene in human promyelocytic leukemia HL-60 cells that contain serglycin mRNA and in T-lymphoblast Molt-4 cells that do not. With polymerase chain reaction methodology, 13 DNA probes of 250-880 base pairs in length were generated that corresponded to unique, non-Alu sequences spaced throughout the entire human serglycin gene. When blots containing genomic DNA digested with HpaII or MspI were examined with these genomic probes, it was discovered that the 5'-flanking region and intron 1 of the serglycin gene in HL-60 cells were both substantially less methylated than intron 2. In contrast, the entire serglycin gene in Molt-4 cells was highly methylated. Because hypomethylated genes generally are transcribed more efficiently than hypermethylated genes, the high level of serglycin mRNA in HL-60 cells probably is a consequence of the low level of methylation of intron 1 and the 5'-flanking region of the serglycin gene in these cells.     

Complete nucleotide sequence of a methylcholanthrene-inducible cytochrome P-450 (P-450d) gene in the rat

The rat cytochrome P-450d gene which is inducibly expressed by the administration of 3-methylcholanthrene (MC) has been cloned and analyzed for the complete nucleotide sequence. The gene is 6.9 kilobases long and is separated into 7 exons by 6 introns. The insertion sites of the introns in this gene are well-conserved as compared with those of another MC-inducible cytochrome P-450c gene, but are completely different from those of a phenobarbital-inducible cytochrome P-450e gene. The overall homologies in the coding nucleotide and deduced amino acid sequences were 75% and 68% between the two MC-inducible cytochrome P-450 genes, respectively. The similarity of the gene organization between cytochrome P-450d and P-450c as well as their homology in the deduced amino acid and the nucleotide sequences suggests that these two genes of MC-inducible cytochromes P-450 constitute a different subfamily than those of the phenobarbital-inducible one in the cytochrome P-450 gene family. In contrast with the notable sequence homology in the coding region of the two MC-inducible cytochromes P-450, all the introns and the 5'- and 3'-flanking regions of the two genes showed virtually no sequence homology between them except for several short DNA segments that are located in the promoter region and the first intron. The nucleotide sequences and the locations of these conserved short DNA segments in the two genes suggest that they may affect the expression of the genes. Middle repetitive sequence reported as ID or identifier sequence were found in and in the vicinity of the cytochrome P-450d gene.     

Structure of the human prealbumin gene

Using cloned human prealbumin cDNA as a probe, Southern blot hybridization of human genomic DNA revealed that the prealbumin gene consists of an unique, single-copy DNA. The nucleotide sequences of the entire human prealbumin gene, including both 581 base pairs of the 5'- and 95 base pairs of the 3'-flanking sequences, were determined. The gene spans about 7.0 kilobase pairs and consists of four exons and three introns. As in most eukaryotic genes, the consensus TATA and CAAT sequences are found 30 and 101 nucleotides, respectively, upstream from the putative cap site, and a polyadenylation signal sequence AA-TAAA is found in the 3'-untranslated region. Unexpectedly, two independent open reading frames provided with respective regulatory sequences were found within the gene: one in the first intron and the other in the third intron.