Latitudinal clines for nucleotide polymorphisms in the Esterase 6 gene of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Previous studies have found non-neutral patterns of nucleotide polymorphism in the promoter and coding regions of Est6 in D. melanogaster. Coding region polymorphism peaks around two closely linked replacement differences associated with the EST6-F/EST6-S allozyme polymorphism. The promoter contains two common, highly diverged haplotype groups, P1 and P7, that differentially affect Est6 expression. Allozyme studies have also revealed latitudinal clines in EST6-F and EST6-S frequencies that recur across continents. Here we analyse nucleotide polymorphisms across the promoter and the region of peak coding sequence polymorphism in 10 Australian populations along a 25° latitudinal gradient in order to examine the basis for the allozyme clines. As with the earlier studies, we find an excess of intermediate to high frequency variants in both the P1/P7 region and around the two EST6-F/EST6-S replacements in some populations. The two EST6-F/EST6-S replacement polymorphisms show latitudinal clines whereas the P1 and P7 groups of promoter haplotypes do not. However the strongest clines are for three co-segregating silent site polymorphisms in a 4 bp stretch at the 3′ end of the sequenced region. Monte Carlo simulations show that the clines for those three sites can explain all others in the data but none of the others can explain those three. Thus the allozyme clines may not reflect selection on either the P1/P7 polymorphism or the two replacements previously associated with the EST6-F/EST-S difference.     

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.     

Nucleotide Variation and Conservation at the Dpp Locus, a Gene Controlling Early Development in Drosophila

A study of polymorphism and species divergence of the dpp gene of Drosophila has been made. Eighteen lines from a population of D. melanogaster were sequenced for 5200 bp of the Hin region of the gene, coding for the dpp polypeptide. A comparison was made with sequence from D. simulans. Ninety-six silent polymorphisms and three amino acid replacement polymorphisms were found. The overall silent polymorphism (0.0247) is low, but haplotype diversity (0.0066 for effectively silent sites and 0.0054 for all sites) is in the range found for enzyme loci. Amino acid variation is absent in the N-terminal signal peptide, the C-terminal TGF-β peptide and in the N-terminal half of the pro-protein region. At the nucleotide level there is strong conservation in the middle half of the large intron and in the 3' untranslated sequence of the last exon. The 3' untranslated conservation, which is perfect for 110 bp among all the divergent species, is unexplained. There is strong positive linkage disequilibrium among polymorphic sites, with stretches of apparent gene conversion among originally divergent sequences. The population apparently is a migration mixture of divergent clades.     

Associations of esterase 6 allozyme and activity variation with reproductive fitness inDrosophila melanogaster

Previous studies have shown that the esterase 6 (EST6) enzyme ofD. melanogaster is mainly produced in the sperm ejaculatory duct of the adult male and comparisons of wild-type males with laboratory null mutants have suggested that the enzyme plays a role in reproductive fitness. In this study we have compared 18 field-derived lines each isoallelic forEst6 for differences in five components of male reproductive fitness. No consistent fitness differences were found among lines differing in respect of the two major allozyme classes EST6-F and EST6-S, despite other evidence that these two classes are not selectively equivalent in the field. However, differences in reproductive fitness were found among lines differing in the minor mobility variants that segregate within EST6-F and EST6-S. A failure to distinguish among these minor forms may explain the discrepancies in previous studies on the effects of the major EST6 allozymes on reproductive fitness. The most significant associations we have found between EST6 and reproductive fitness were due to variation in EST6 activity levels. Male EST6 activity levels were found to be positively correlated with their time to first mating, negatively correlated with the numbers of eggs laid and progeny produced by their mates, and negatively correlated with the frequency with which their mates remate. We conclude that some EST6 variants differ in components of male reproductive fitness operative in laboratory cultures. However, the evidence for fitness differences is stronger for variants affecting the amount, rather than the structure of the enzyme, and the direction of the differences varies between some of the fitness components tested.     

Historical Selection, Amino Acid Polymorphism and Lineage-Specific Divergence at the G6pd Locus in Drosophila Melanogaster and D. Simulans

The nucleotide diversity across 1705 bp of the G6pd gene is studied in 50 Drosophila melanogaster and 12 D. simulans lines. Our earlier report contrasted intraspecific polymorphism and interspecific differences at silent and replacement sites in these species. This report expands the number of European and African lines and examines the pattern of polymorphism with respect to the common A/B allozymes. In D. melanogaster the silent nucleotide diversity varies 2.8-fold across localities. The B allele sequences are two- to fourfold more variable than the derived A allele, and differences between allozymes are twice as among B alleles. There is strong linkage disequilibrium across the G6pd region. In both species the level of silent polymorphism increases from the 5' to 3' ends, while there is no comparable pattern in level of silent site divergence or fixation. The neutral model is not rejected in either species. Using D. yakuba as an outgroup, the D. melanogaster lineage shows a twofold greater rate of silent fixation, but less than half the rate of amino acid replacement. Lineage-specific differences in mutation fixation are inconsistent with neutral expectations and suggest the interaction of species-specific population size differences with both weakly advantageous and deleterious selection.     

Inferring the Evolutionary Histories of the Adh and Adh-Dup Loci in Drosophila Melanogaster from Patterns of Polymorphism and Divergence

The DNA sequences of 11 Drosophila melanogaster lines are compared across three contiguous regions, the Adh and Adh-dup loci and a noncoding 5' flanking region of Adh. Ninety-eight of approximately 4750 sites are segregating in the sample, 36 in the 5' flanking region, 38 in Adh and 24 in Adh-dup. Several methods are presented to test whether the patterns and levels of polymorphism are consistent with neutral molecular evolution. The analysis of within- and between-species polymorphism indicates that the region is evolving in a nonneutral and complex fashion. A graphical analysis of the data provides support for a hypothesized balanced polymorphism at or near position 1490, site of the amino acid replacement difference between Adhf and Adhs. The Adh-dup locus is less polymorphic than Adh and all 24 of its polymorphisms occur at low frequency--suggestive of a recent selective substitution in the Adh-dup region. Adhs alleles form two distinct evolutionary lineages that differ one from another at a total of nineteen sites in the Adh and Adh-dup loci. The polymorphisms are in complete linkage disequilibrium. A recombination experiment failed to find evidence for recombination suppression between the two allelic classes. Two hypotheses are presented to account for the widespread distribution of the two divergent lineages in natural populations. Natural selection appears to have played an important role in governing the overall patterns of nucleotide variation across the two-gene region.     

Nucleotide polymorphism in the Est6 promoter,which is widespread in derived populations of Drosophila melanogaster,changes the level of Esterase 6 expressed in the male ejaculatory duct

Previous analysis of an Australian population of D. melanogaster revealed two predominant Est6 promoter haplotypes, P1 and P7. These haplotypes, which differ at 14 sites over a 325-bp region, are associated with a 15-20% difference in male EST6 activity. Here we show that the P1/P7 sequence difference causes the male activity variation by recreating the activity difference among >60 independently transformed lines containing representative P1 or P7 promoter alleles fused to an identical Est6 coding region. Furthermore we find that the whole fly difference reflects about a twofold difference in EST6 activity in the anterior sperm ejaculatory duct. EST6 activity variation in this tissue is known to affect reproductive fitness. Using a combination of RFLP analysis and DNA sequencing, we show that P1 and P7 are predominant in six populations from America, Asia, and Australia, albeit less frequent in a population from the presumptively ancestral east African range of the species. The sequence data show significant departures from neutral expectations for the derived American and Australian populations but not the presumptively ancestral Zimbabwean population. Thus the P1/P7 difference could be a major source of adaptively significant EST6 activity variation through much of the now cosmopolitan range of D. melanogaster.     

Population variation in linkage disequilibrium across the COMT gene considering promoter region and coding region variation

Catechol-O-methyl transferase (COMT) catalyzes the first step in one of the major pathways in the degradation of catecholamines. The COMT gene on chromosome 22 has been considered a candidate gene for many neuropsychiatric disorders, in part because an exon 4 single nucleotide polymorphism (SNP) in COMT causes an amino acid substitution associated with significantly altered enzyme activity. This functional variant, detected as an NlaIII restriction site polymorphism (RSP), is polymorphic in populations from around the world. A four-site haplotype spanning 28 kb effectively encompasses COMT. This haplotype is comprised of two novel polymorphisms [a tetranucleotide short tandem repeat polymorphism (STRP) in intron 1 and a HindIII RSP at the 5' end of COMT], the NlaIII site, and another previously published site - a BglI RSP at the 3' end of the gene. Overall linkage disequilibrium (LD) for this haplotype is strong and significant in 32 population samples from around the world. Conditional probabilities indicate that, in spite of moderate to strong disequilibrium in most non-African populations, the NlaIII site, although often used for prediction, would not always be a reliable predictor of allelic variation at the other sites. Because other functional variation might exist, especially regulatory variation, these findings indicate that haplotypes would be more effective indicators of possible involvement of COMT in disease etiology.     

Correlation between patterns of DNase I-hypersensitive sites and upstream promoter activity of the human epsilon-globin gene at different stages of erythroid development.

DNA 5' to the human epsilon-globin gene exhibits unique patterns of DNase I-hypersensitive sites (DHS) in three human erythroleukemic cell lines which represent the embryonic (K562), fetal (HEL), and adult (KMOE) stages of erythroid development. We have mapped 10 epsilon-globin DHS in K562 cells, in which the epsilon-globin gene is maximally active. Major sites are located -11.7, -10.5, -6.5, -2.2 kilobase pairs (kbp) and -200 base pairs (bp) upstream of the gene and directly over the major cap site. Minor sites are located -5.5, -4.5, and -1.48 kbp and -900 bp upstream of the cap site. In HEL cells, in which the epsilon-globin gene is expressed at extremely low levels, the -11.7-, -10.5-, -5.5-, -4.5-, and -2.2-kbp DHS are no longer detectable; the -200-bp site is approximately 300-fold less sensitive to DNase I; and the -1.48-kbp, -900-bp, and major cap site DHS are 3- to 4-fold less sensitive. Only the DHS located -6.5 kbp relative to the major cap site is detectable at all three stages of erythroid development, including KMOE cells in which epsilon-globin synthesis is undetectable. We suggest that this site may be implicated in maintaining the entire beta-globin cluster in an active chromatin conformation. The five DHS downstream of the -6.5-kbp element possess associated promoters. Thus two distinct types of DHS exist--promoter positive and promoter negative. In HEL cells, all the upstream promoters are inactivated, although the -1.48-kbp and -900- and -200-bp DHS are still present. This suggests that the maintenance of DHS and regulation of their associated promoters occur by independent mechanisms. The inactivation of the upstream promoters in HEL cells while the major cap site remains active represents a unique pattern of expression and suggests that HEL cells possess regulatory factors which specifically down regulate the epsilon-globin upstream promoters.     

Insertion elements and deletion formation in a halophilic archaebacterium.

Deletion events that occur spontaneously in 36-kilobase-pair (kbp) plasmid pHH4 from the archaebacterium Halobacterium halobium were investigated. Four different deletion derivatives with sizes ranging from 5.7 to 17 kbp were isolated. Three of these deletion variants derived from pHH4 (pHH6 [17 kbp], pHH7 [16 kbp], and pHH8 [6.3 kbp]), whereas the 5.7-kbp plasmid pHH9 derived from pHH6. Strains containing pHH6, pHH7, or pHH9 each lacked the parental plasmid pHH4, while pHH8 occurred at a 1:1 ratio together with pHH4. Common to all of these plasmids was the 5.7-kbp region of pHH9 DNA. The regions containing the fusion site in the deletion derivatives were investigated and compared with the corresponding area of the parental plasmid. Each deletion occurred exactly at the terminus of an insertion element. In pHH6 and pHH7, a halobacterial insertion element (ISH2) was located at the deletion site. The DNA fused to ISH2 displayed a 7-base-pair (bp) (pHH7) or 10-bp (pHH6) sequence homology to the inverted repeat of ISH2. In the two smaller plasmids, pHH8 and pHH9, an ISH27 element was located at the deletion site. Most likely, all of these smaller plasmids resulted from an intramolecular transposition event. The ISH27 insertion sequence contains a 16-bp terminal inverted repeat and duplicates 5 bp of target DNA during the transposition with the specificity 5'ANNNT3'. Four ISH27 copies were analyzed, and two ISH27 element types were identified that have approximately 85% sequence similarity. The ISH27 insertion elements constitute a family which is related to the ISH51 family characterized for H. volcanii, another halophilic archaebacterium.     

Functional organization of the Harvey murine sarcoma virus genome.

The comparative infectivity of Harvey murine sarcoma virus (Ha-MuSV) DNA for NIH 3T3 cells was determined for supercoiled Ha-MuSV DNA molecularly cloned in lambda phage and pBR322 at its unique EcoRI site (which is located near the middle of the 6-kilobase pair [kbp] unintegrated linear viral DNA) and for two cloned subgenomic fragments: one was 3.8 kbp and lacked about 1 kbp from each side of the EcoRI site, and the second did not contain the 3 kbp of the unintegrated linear viral DNA located on the 3' side of the EcoRI site. Each subgenomic DNA induced foci of transformed cells, but with a lower relative efficiency then genomic DNA. Transfection with intact vector Ha-MuSV DNA yielded results similar to those obtained after separation of Ha-MuSV DNA from vector DNA. Cells lines were then derived from individual foci transformed with each type of viral DNA. Focus-forming virus was recovered from transformed cells after superinfection with a helper-independent virus, but the efficiency varied by several orders of magnitude. For several transformed lines, the efficiency of recovery of focus-forming virus was correlated with the structure of the Ha-MuSV DNA in the cells before superinfection. When 32P-labeled Ha-MuSV DNA probes specific for sequences on either the 3' or 5' side of the EcoRI site were used to analyze the viral RNA in the transformed cell lines, all lines were found to hybridize with the 5' probe, but some lines did not hybridize with the 3' probe. The transformed lines contained high levels of the Ha-MuSV-coded p21 or its associated GDP-binding activity. We conclude that the transforming region and the sequences that code for the viral p21 protein are both located within the 2 kilobases closest to the 5' end of the Ha-MuSV genome.     

Nuclease hypersensitivity in the beta-globin gene region of K562 cells

We have investigated chromatin structure in the beta-globin gene region of the K562 human erythroleukemic cell line by using S1 and DNase I nuclease sensitivity assays. Despite the lack of beta-globin gene expression in these cells, we find nuclease-hypersensitive sites to these enzymes in its 5' and 3' flanking regions in K562 chromatin. This result is in contrast to previous reports in which no hypersensitive sites were found in the immediate vicinity of this gene. In the 3' region, one major hypersensitive site at 0.9 kpb 3' and three minor hypersensitive sites at 0.7 kbp, 0.5 kbp 3' and 0.2 kbp 5' of the polyadenylation site were observed; these sites are very similar to those found in fetal liver and adult bone marrow cells in which the beta-globin gene is expressed. We find hypersensitive sites to both enzymes in the 5' region of the beta-globin gene: a major site 0.8 kbp 5' to the cap site, and two minor sites 1.2 and 1.5 kbp 5' to the cap site. The -0.8 kbp site is also present in plasmids containing the beta-globin gene. Our results suggest that the lack of beta-globin gene expression may be related to the lack of hypersensitivity sites in the immediate (150 bp) 5' flanking region of the beta-globin gene, as occurs in other active globin genes.