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.     

Nucleotide variation of the Est-6 gene region in natural populations of Drosophila melanogaster

We have investigated nucleotide polymorphism in the Est-6 gene region in four samples of Drosophila melanogaster derived from natural populations of East Africa (Zimbabwe), Europe (Spain), North America (California), and South America (Venezuela). There are two divergent sequence types in the North and South American samples, which are not perfectly (North America) or not at all (South America) associated with the Est-6 allozyme variation. Less pronounced or no sequence dimorphism occurs in the European and African samples, respectively. The level of nucleotide diversity is highest in the African sample, lower (and similar to each other) in the samples from Europe and North America, and lowest in the sample from South America. The extent of linkage disequilibrium is low in Africa (1.23% significant associations), but much higher in non-African populations (22.59, 21.45, and 37.68% in Europe, North America, and South America, respectively). Tests of neutrality with recombination are significant in non-African samples but not significant in the African sample. We propose that demographic history (bottleneck and admixture of genetically different populations) is the major factor shaping the nucleotide patterns in the Est-6 gene region. However, positive selection modifies the pattern: balanced selection creates elevated levels of nucleotide variation around functionally important (target) polymorphic sites (RsaI-/RsaI+ in the promoter region and F/S in the coding region) in both African and non-African samples; and directional selection, acting during the geographic expansion phase of D. melanogaster, creates an excess of very similar sequences (RsaI- and S allelic lineages, in the promoter and coding regions, respectively) in the non-African samples.     

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster.

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.     

Nucleotide variation at the hypervariable esterase 6 isozyme locus of Drosophila simulans

Esterase 6 (Est-6/EST6) is polymorphic in both Drosophila melanogaster and D. simulans for two common allozyme forms, as well as for several other less common variants. Parallel latitudinal clines in the frequencies of the common EST6-F and EST6-S allozymes in these species have previously been interpreted in terms of a shared amino acid polymorphism that distinguishes the two variants and is subject to selection. Here we compare the sequences of four D. simulans Est-6 isolates and show that overall estimates of nucleotide heterozygosity in both coding and 5' flanking regions are more than threefold higher than those obtained previously for this gene in D. melanogaster. Nevertheless, the ratio of replacement to exon silent-site polymorphism in D. simulans is less than the ratio of replacement to silent divergence between D. simulans and D. melanogaster, which could be the result of increased efficiency of selection against replacement polymorphisms in D. simulans or to divergent selection between the two species. We also find that the amino acid polymorphisms separating EST6- F and EST6-S in D. simulans are not the same as those that separate these allozymes in D. melanogaster, implying that the shared clines do not reflect shared molecular targets for selection. All comparisons within and between the two species reveal a remarkable paucity of variation in a stretch of nearly 400 bp immediately 5' of the gene, indicative of strong selective constraint to retain essential aspects of Est-6 promoter function.      

香蕉果实成熟相关基因ACO1启动子区的克隆及其功能初探

根据已报道的香蕉课实表达ACC氧化酶基因（ACO1）的序列,用改进的接头连接PCR法从香蕉基因组中扩增并克隆了此基因5′旁侧区1526bp的片段,其中包含一个推测的TATA盒序列;与已公布的两个香蕉ACC氧化酶基因启动子序列（分别为934bp和1451bp）的相似性各为97．3％（Lopez-Gomez等）和88．8％（May和Kipp)。将4个含有不同大小启动子区的克隆片段与GUS基因编码区连接构建成嵌合成基因,通过基因枪轰击转入香蕉叶、根和果实的细胞后,瞬时表达结果表明不同大小的ACO1启动子区段都只在果实细胞中指导GUS基因表达,证明该启动子具有指导基因在果实中表达的功能,并推测负责果实特异性的顺式元件可能位于启动子近端0．7kb区段之内,在468至822的355bp区段内可能在与正控制有关的顺式元件。     

Nucleotide sequence of the B1-hordein gene of barley Hordeum vulgare L

The gene encoding B1 hordein of Hordeum vulgare (cv. Donetsky 4) was cloned and entirely sequenced. It contains no introns and codes for of 293 amino acids long polypeptide with molecular weight 33418. Our clone differs from the previously sequenced B1 hordein genes in some positions within the coding region (there are 4 nucleotide changes and a 12 bp deletion, as compared with the pB11 cDNA clone, and 5 nucleotide changes, as compared with the pBHP184 genomic clone). These changes result in polymorphism of amino acid sequences at 5 positions. 5'-flanking region contains putative regulatory and promoter sequences and differs from that of the pBHP184 clone in 3 positions.     

Molecular population genetics of the<Emphasis Type="Italic">β-esterase</Emphasis> gene cluster of<Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

We have investigated nucleotide polymorphism at theβ-esterase gene cluster including theEst-6 gene andψEst-6 putative pseudogene in four samples ofDrosophila melanogaster derived from natural populations of southern Africa (Zimbabwe), Europe (Spain), North America (USA: California), and South America (Venezuela). A complex haplotype structure is revealed in bothEst-6 andψEst-6. Total nucleotide diversity is twice inψEst-6 as inEst-6; diversity is higher in the African sample than in the non-African ones. Strong linkage disequilibrium occurs within theβ-esterase gene cluster in non-African samples, but not in the African one. Intragenic gene conversion events are detected withinEst-6 and, to a much greater extent, withinyEst-6; intergenic gene conversion events are rare. Tests of neutrality with recombination are significant for theβ-esterase gene cluster in the non-African samples but not significant in the African one. We suggest that the demographic history (bottleneck and admixture of genetically differentiated populations) is the major factor shaping the pattern of nucleotide polymorphism in theb-esterase gene cluster. However there are some ’footprints’ of directional and balancing selection shaping specific distribution of nucleotide polymorphism within the cluster. Intergenic epistatic selection betweenEst-6 andψEst-6 may play an important role in the evolution of theβ-esterase gene cluster preserving the putative pseudogene from degenerative destruction and reflecting possible functional interaction between the functional gene and the putative pseudogene.Est-6 andyEst-6 may represent an indivisible intergenic complex (‘intergene’) in which each single component (Est-6 orψEst-6) cannot separately carry out the full functional role.     

Nucleotide Polymorphism in the 5'' Promoter Region of Esterase 6 in Drosophila Melanogaster and Its Relationship to Enzyme Activity Variation

A 974-bp region immediately 5' of the esterase 6 gene was sequenced in 17 field derived third chromosome isoallelic lines. Twenty-three polymorphisms were identified, only two in the first 400 bp 5' but 16 in a 325-bp region from -494 to -819 bp. This distribution differs from previously published patterns in Drosophila simulans and D. mauritiana, where the first 800 bp are highly conserved. Fourteen common polymorphisms in the 325-bp region above are all in strong linkage disequilibrium with each other. Moreover, most of the haplotypes defined by the total of 23 polymorphisms fall into two groups that differ as a block at all 14 of these latter sites. Sequence differences between the two groups include some restriction sites that were scored in an earlier study of RFLPs and EST6 enzyme phenotypes among 42 isoallelic lines from the same population. By collating the two studies, we show that one haplotype group yields ~15% lower EST6 enzyme activity in adult males than the other. The promoter haplotypes show only weak disequilibrium with the esterase 6 fast/slow allozyme polymorphism, so it seems unlikely that previously reported latitudinal clines in the allozyme frequencies are due to their hitchhiking along with selection on the promoter difference.     

Human plasminogen activator inhibitor-1 gene. Promoter and structural gene nucleotide sequences

We have determined the nucleotide sequence of the human plasminogen activator inhibitor-1 (PAI-1) gene and significant stretches of DNA which extend into its 5'-and 3'-flanking DNA regions; a total sequence of 15,867 base pairs (bp) is presented. The sequenced 5'-flanking DNA (1,520 bp) contains the essential eukaryotic cis-type proximal regulatory elements CCAAT and TATAA; the more distal 5'-flanking DNA region, as well as some introns, contain sequence elements which share identities with known eukaryotic enhancer elements. A major finding is the identification of a large region of shared nucleotides (comprising of about 520 bp) between the 5'-flanking DNAs of PAI-1 and tissue-type plasminogen activator genes. The length of the PAI-1 5'-untranslated region was found to be 145 bp as determined by nuclease analysis. The remaining PAI-1 structural gene consists of amino acid coding regions (containing a total of 1,206 bp, coding for the 23 amino acids of the signal peptide and 379 amino acids of the mature PAI-1 protein), 8 intron regions (a total of 8,978 bp), and a long 3'-untranslated region of about 1,800 bp which contains several polyadenylation sites. Two types of repetitive DNA elements are located within the PAI-1 structural gene and flanking DNAs: we have found 12 Alu elements and 5 repeats of a long poly (Pur) element. These Alu-Pur elements may represent a subset of the more abundant Alu family of repetitive sequence elements.     

Nucleotide variation of the duplicated amylase genes in Drosophila kikkawai

We examined levels and patterns of the nucleotide polymorphism of the Amylase genes with a head-to-head duplication in Drosophila kikkawai. The levels of variation in D. kikkawai were comparable to those in Drosophila melanogaster. Tajima's test, Fu and Li's test, HKA test, and MK test did not show significant departure from neutrality. We found an excess of replacement changes in the within-locus class, representing polymorphism in one of the duplicated genes, compared with the between-locus class, representing polymorphism shared between the duplicated genes. Most replacement changes in the within-locus class were singletons. These results suggest that most replacement changes are deleterious. A contrasting evolutionary pattern, involving concerted evolution in the coding regions but differential evolution in the 5'-flanking regions, was observed. However, unlike the duplicated Amy genes of D. melanogaster, the coding regions of the duplicated genes in D. kikkawai tended to diverge. Using Ohta's model of the small multigene family, we found that recombination (interchromosomal equal crossing-over) rate was one order higher than gene conversion (unequal crossing-over) rate, resulting in a considerable but incomplete homogenization of the duplicated coding regions. Linkage disequilibria were found in the intron as well as within and around the regulatory cis-element sequences of one of the duplicated genes (Amy1). The possible causes of these linkage disequilibria were discussed.     

Molecular evolution of the 5'-flanking regions of the duplicated Amy genes in Drosophila melanogaster species subgroup

The nucleotide sequences of the 5'-flanking regions of the duplicated Amy genes in eight sibling species belonging to the melanogaster species subgroup are analyzed. In Drosophila melanogaster, a region of about 450 bp immediately upstream of the translation initiation site of the two paralogous genes (the proximal and distal genes) has sequence similarities. However, we could not detect any significant sequence similarity in the region more upstream than -450. This result indicates that the coding regions of the ancestral Amy gene were duplicated together with 450 bp of the 5'-flanking region as one unit. Multiple alignment of these 450-bp sequences in the proximal and distal genes of all eight species revealed a mosaic pattern of highly conserved and divergent regions. The conserved regions included almost all the putative regulatory elements identified in previous analyses of the sequences. A phylogenetic analysis of the aligned sequences shows that these 450-bp sequences are clustered into the proximal and the distal groups. As a whole, the divergence between groups in this region is very large in contrast to that in the coding regions. Based on the divergence between groups, the 450-bp region is divided into two subregions. We found that the ratios of the divergence between groups to that within groups differ in the two subregions. From these observations, we discuss a possibility of positive selection acting on the subregion immediately upstream of the Amy coding region to cause divergence of regulatory elements of the paralogous genes.      

The β-esterase Gene Cluster of Drosophila melanogaster: is ψEst-6 a Pseudogene, a Functional Gene, or both?

Pseudogenes have been defined as non-functional sequences of genomic DNA that are originally derived from functional genes, but exhibit degenerative features such as premature stop codons and frameshifts that prevent their expression. However, there is increasing evidence that pseudogenes are often evolutionarily conserved and may have retained some functional role or acquired new ones. Pseudogenes may exhibit non-functional features as well as functional ones. We investigate, as a model case, the beta-esterase gene cluster of Drosophila melanogaster that includes the Est-6 gene and the psiEst-6 putative pseudogene. We study four samples derived from natural populations of east Africa (Zimbabwe), Europe (Spain), North America (California), and South America (Venezuela). The level of nucleotide diversity is higher in Africa than in the non-African populations. There is twice more nucleotide diversity in psiEst-6 than in Est-6. Linkage disequilibrium within the beta-esterase gene cluster is strong in non-African samples, but much lower in Africa. The population recombination rate is the same for psiEst-6 and Est-6 in Africa, but significantly different in non-African samples. Intragenic gene conversion events are detected within Est-6 and, with much higher incidence, within psiEst-6; intergenic gene conversion events are rare. The extensive intragenic gene conversion within psiEst-6 can be explained by the invasion of retrotransposons that promote a form of homology-dependent gene conversion upon excision. Tests of neutrality with recombination are significant for the beta-esterase gene cluster in the non-African populations but not in Africa. The Est-6 gene sequences exhibit a well-known allozyme dimorphic structure. The sequences of psiEst-6 are also dimorphic in North and South America, but they do not correspond at all (South America) or only imperfectly (North America) to the Est-6 allozyme dimorphism. Sequence dimorphism is less pronounced in the European and African samples. We suggest that demographic history (bottleneck and admixture of genetically differentiated populations) is the major factor shaping the nucleotide pattern in the beta-esterase gene cluster. However, there are some clear indications of positive selection shaping the distribution of nucleotide polymorphism within the cluster. Intergenic epistatic selection may play an important role in the evolution of the beta-esterase gene cluster, preserving psiEst-6 from degenerative destruction and reflecting its functional interaction with Est-6. The Est-6 gene cluster of D. melanogaster represents an example of a functionally interacting complex ('intergene') in which two components (Est-6 and psiEst-6) or more are required to perform the final function.     

Genetic polymorphism at two linked loci,Sod and Est-6, in Drosophila melanogaster

We have examined the patterns of polymorphism at two linked loci, Sod and Est-6, separated by nearly 1000 kb on the left arm of chromosome 3 of Drosophila melanogaster. The evidence suggests that natural selection has been involved in shaping the polymorphisms. At the Sod locus, a fairly strong (s>0.01) selective sweep, started ≥2600 years ago, increased the frequency of a rare haplotype, F(A), to about 50% frequency in populations of Europe, Asia, and the Americas. More recently, an F(A) allele mutated to an S allele, which has increased to frequencies 5–15% in populations of Europe, Asia and North America. All S alleles are identical (or very nearly) in sequence and differ by one nucleotide substitution (which accounts for the F→S electrophoretic difference) from F(A) alleles. At the Est-6 locus, the evidence indicates both directional and balancing selection impacting separately the promoter and the coding regions of the gene, with linkage disequilibrium occurring within each region. Some linkage disequilibrium also exists between the two genes.