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.     

Analysis of hABC1 gene 5' end: additional peptide sequence, promoter region, and four polymorphisms

Evidence linking mutations in ATP-binding-cassette transporter gene 1 (ABC1) to Tangier disease suggests it functions in the active transport of free cholesterol out of cells. Since its mRNA level is regulated in response to cellular cholesterol stores it is of interest to explore its promoter response elements, and to investigate polymorphisms for their contributions to the prevalence of low levels of HDL in the population that promotes premature coronary heart disease. Investigation of the 5' end of the gene by 5' RACE analysis revealed 455 nucleotides additional to published sequences, and predicts another 60 amino acid N-terminal residues, resulting in a 2261-residue protein. Protein sequence analysis predicts a membrane-spanning region and possible signal peptide. The 5' flanking region was located by a Human Research Project BLAST search. This region contains regulatory elements that potentially control ABC1 gene expression. In addition to numerous SP1 binding sites there are four putative sterol regulatory elements (SREs). Our studies uncovered three single nucleotide substitution polymorphisms, one in the promoter region and two in the 5' untranslated region (5'UTR), plus an insertion/deletion polymorphism.     

Intraspecific nuclear DNA variation in Drosophila

We have summarized and analyzed all available nuclear DNA sequence polymorphism studies for three species of Drosophila, D. melanogaster (24 loci), D. simulans (12 loci), and D. pseudoobscura (5 loci). Our major findings are: (1) The average nucleotide heterozygosity ranges from about 0.4% to 2% depending upon species and function of the region, i.e., coding or noncoding. (2) Compared to D. simulans and D. pseudoobscura (which are about equally variable), D. melanogaster displays a low degree of DNA polymorphism. (3) Noncoding introns and 3' and 5' flanking DNA shows less polymorphism than silent sites within coding DNA. (4) X-linked genes are less variable than autosomal genes. (5) Transition (Ts) and transversion (Tv) polymorphisms are about equally frequent in non-coding DNA and at fourfold degenerate sites in coding DNA while Ts polymorphisms outnumber Tv polymorphisms by about 2:1 in total coding DNA. The increased Ts polymorphism in coding regions is likely due to the structure of the genetic code: silent changes are more often Ts's than are replacement substitutions. (6) The proportion of replacement polymorphisms is significantly higher in D. melanogaster than in D. simulans. (7) The level of variation in coding DNA and the adjacent noncoding DNA is significantly correlated indicating regional effects, most notably recombination. (8) Surprisingly, the level of polymorphism at silent coding sites in D. melanogaster is positively correlated with degree of codon usage bias. (9) Three proposed tests of the neutral theory of DNA polymorphisms have been performed on the data: Tajima's test, the HKA test, and the McDonald-Kreitman test. About half of the loci fail to conform to the expectations of neutral theory by one of the tests. We conclude that many variables are affecting levels of DNA polymorphism in Drosophila, from properties of nucleotides to population history and, perhaps, mating structure. No simple, all encompassing explanation satisfactorily accounts for the data.      

Isolation and characterization of cDNA clones for human skeletal muscle alpha actin.

Two cDNA libraries corresponding to polyA+ RNA from human adult skeletal muscle have been constructed by cloning in the PstI site of pBR322. Skeletal alpha actin cDNA clones have been isolated and characterized. Three of these plasmids have overlapping inserts which together contain the complete 5' non-coding and protein-coding region and part of the 3' untranslated region. Determination of the sequence of the cloned cDNA confirms the complete conservation in human of the amino-acid sequence of skeletal alpha actin compared to the rabbit or rat proteins. The 5' untranslated region, but not the 3' untranslated region, shows good homology with the corresponding one in the rat gene. Analysis of changes at silent sites within the protein-coding region suggests that the divergence of skeletal and cardiac alpha actin took place much earlier than the mammalian radiation. The plasmids described here have been used as probes to detect the homologous gene among the about thirty actin sequences present in the human genome.     

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.      

Nucleotide sequence comparison of the rp49 gene region between Drosophila subobscura and D. melanogaster

A 1.6-kb fragment encompassing the rp49 gene, which codes for a ribosomal protein, has been cloned and sequenced in Drosophila subobscura. The rp49 coding region has accumulated 46 nucleotide differences out of 402 bp since D. subobscura diverged from D. melanogaster. Forty-three percent of the effectively silent sites have changed since both species diverged. Both silent and replacement differences are distributed at random between the two exons of the gene. The frequency of silent differences in exons does not differ from that observed in the 5' leader sequence and in the intron. The frequency of silent differences in exon and intron sites is much greater than the number of amino acid replacement differences. This observation indicates strong purifying selection against amino acid replacements.      

Molecular cloning of rat glucose-dependent insulinotropic peptide (GIP).

A cDNA clone encoding glucose-dependent insulinotropic peptide (GIP) was identified that consisted of 34 bp of 5' untranslated sequence, an open reading frame of 432 bp and 115 bp in the 3' untranslated region. The deduced amino acid sequence revealed a 144 amino acid preprohormone consisting of a 43 amino acid N-terminal extension including a signal peptide, a 42 amino acid hormone, and a 59 amino acid C-terminal extension. Rat GIP differs from the human hormone by two amino acid substitutions: arginine for histidine at position 18 and leucine for isoleucine at position 40. A single mRNA from small intestine of approximately 800 bases was identified on Northern blot analysis in equivalent amounts in proximal and distal small intestine.     

Polymorphisms in the human DNA repair gene XPF.

DNA sequence polymorphisms were sought in the coding region and at the exon-intron boundaries of the human XPF gene, which plays a role in nucleotide excision repair. Based on a survey of 38 individuals, we found six single nucleotide polymorphisms, one in the 5' non-coding region of the XPF gene, and five in the 2751 bp coding region. At each site, the frequency of the rarer allele varies from about 0.01 to over 0.38. Except for the 5' non-coding and one coding sequence polymorphism, the rarer alleles for the remaining four polymorphisms were found only in heterozygotes. Of the five polymorphisms in the coding region, one is silent, one results in a conserved amino acid difference, and the remaining three result in non-conserved amino acid differences. Because of its biological function in nucleotide excision repair, functionally significant XPF gene polymorphisms are candidates for influencing cancer susceptibility and overall genetic stability. Nucleotide sequence diversity estimates for XPF are similar to the lipoprotein lipase and beta-globin genes.     

糜子抗旱节水相关基因PmMYB的克隆及表达分析

根据在糜子抗旱节水分子基础研究中获得的一个糜子MYB基因的EST序列, 以其序列及水稻MYB18基因的序列为基础设计引物, 扩增得到1 739 bp的全长基因组序列。序列分析表明, 其包含121 bp(347~467 bp)和93 bp(599~691 bp)的两个内含子, 3个外显子; 全长cDNA序列为1 525 bp, 其中3′非翻译区为212 bp, 5′非翻译区为41 bp, 编码区为1 272 bp, 共编码424个氨基酸, C-端存在一个丝氨酸(Ser, S)丰富区。该基因具有两个典型的MYB类转录因子基因的DNA结合区(DNA-binding domain), 分别为13~63、66~114位氨基酸, 属于典型的R2R3-MYB转录因子。对其与水稻、玉米、火炬松、拟南芥、辣椒、陆地棉、大麦及茄子等9种植物的MYB基因的R2、R3重复区的氨基酸序列多重比较, 表明R2R3重复序列在植物中具有较高的保守性; 基于氨基酸序列的编码区系统进化树分析表明, 不同植物的MYB基因遗传分化很大, 序列相似性为32%~84%, 其中糜子MYB基因与水稻的MYB18相似程度最高(84%), 与大麦和玉米的相似性分别为46%和41%。通过半定量RT-PCR对其表达模式分析表明, 该基因在水分胁迫和干旱后复水条件下上调表达, 与糜子抗旱节水紧密相关。该基因的克隆为进一步探讨利用该基因改良其他植物的抗旱节水性奠定了良好的基础。     

Gene for selenium-dependent glutathione peroxidase maps to human chromosomes 3, 21 and X

A human glutathione peroxidase cDNA has been used as a probe to hybridize to DNAs isolated from human - rodent somatic cell hybrids that have segregated human chromosomes. A 609 bp probe which contains the entire coding region hybridizes to human chromosomes 3, 21 and Xp. Fragments of the cDNA coding sequence and of the 3' untranslated region were also used as probes. These fragments hybridized to each of the three chromosomes with the same efficiency, suggesting similarity between the loci, whereas an intronic probe detected only the gene on chromosome 3. The general organization of each gene was determined from the hybridization data. The data suggest that the locus on chromosome 3 is a functional gene containing a single intron and a pattern of restriction sites identical to those found in the cDNA coding sequence. The data also suggest that the sequences on chromosomes X and 21 have equal conservation of the 3' untranslated and coding sequences but do not contain introns, providing evidence that the latter two sequences are processed pseudogenes. A simple two allele polymorphism in PvuII digests was detected at the locus on chromosome 21.     

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.     

Isolation and sequence analysis of a full-length cDNA for human M creatine kinase

A full length cDNA for human M creatine kinase has been isolated and sequenced. The cDNA contains 77 bp of 5' untranslated, 338 bp of 3' untranslated sequence and the entire coding region (1146 bp) for human M creatine kinase. The M creatine kinases from different species share considerable sequence homology within the coding region (77-91%) and in amino acid sequence (82-97%). Little or no sequence homology is observed in the 3' untranslated sequence of the mammalian M creatine kinases, although canine and human creatine kinase share overall 80% sequence homology in 5' untranslated sequence. A unique 8 bp sequence was identified in the 5' untranslated regions of mammalian M creatine kinase but is not present in B creatine kinase cDNA. The degree of sequence conservation observed implies an evolutionary constraint on M creatine kinase structure beyond that which would be expected for the maintenance of enzymatic function.     

Sequences of mouse immunoglobulin light chain genes before and after somatic changes.

We have determined the nucleotide sequences of the germ line gene as well as a corresponding somatically mutated and rearranged gene coding for a mouse immunoglobulin lambdaI type light chain. These sequencing studies were carried out on three Eco RI-DNA fragments which had been cloned from BALB/c mouse embryos or a lambdaI chainsecreting myeloma, H2020. The embryonic DNA clone Ig 99lambda contains two protein-encoding segments, one for the majority of the hydrophobic leader (L) and the other for the rest of the leader and the variable (V) region of the lambda0 chain (Cohn et al., 1974); these segments are separated by a 93 base pair (bp) intervening sequence (I-small). The coding of the V region ends with His at residue 97. The second embryonic DNA clone Ig 25lambda includes a 39 bp DNA segment (J) coding for the rest of the conventionally defined V region (that is, up to residue 110), and also contains the sequence coding for the constant (C) region approximately 1250 untranslated bp (I-large) away from the J sequence. The J sequence is directly linked with the V-coding sequence in the myeloma DNA clone, Ig 303lambda, which has the various DNA segments arranged in the following order: 5' untranslated region, L, l-small, V linked with J, l-large, C, 3' untranslated sequence. The lg 303lambda V DNA sequence codes for the V region synthesized by the H2020 myeloma and is different from the lg 99lambda V DNA sequence by only two bases. No silent base change was observed between the two DNA clones for the entire sequence spanning the 5' untranslated regions and the V-coding segments. These results confirm the previously drawn conclusion that an active complete lambdaI gene arises by somatic recombination that takes place at the ends of the V-coding DNA segment and the J sequence. No sequence homology was observed at or near the sites of the recombination.     

Inference of positive and negative selection on the 5' regulatory regions of Drosophila genes

Both positive selection and negative selection have been shown to drive the evolution of coding regions. It is of interest to know if the corresponding 5' regions of genes may be subjected to selection of comparable intensities. For such a comparison, we chose the Accessory gland protein (Acp) genes as our test case. About 700 bp and 600 bp for the 5' and coding regions, respectively, of eight previously unstudied genes were sequenced from 21 isogenic lines of D. melanogaster and one line from D. simulans. The ratio of divergence at the amino-acid replacement sites (A) over that at the synonymous sites (S) was twice the ratio for common polymorphism. Interestingly, the 5' region shows the same trend, with the 5'/S divergence ratio being 1.8 times higher than the 5'/S ratio for common polymorphism. There are several possible explanations for the 5'/S ratios, including demography, negative selection, and positive selection. Under normal conditions, positive selection is the most likely explanation. If that is true, about 45 to 50 percent of all fixed differences at both the replacement and 5' sites were adaptive, even though the substitution rate in the former is only half that of the latter (K(A)/K(S) approximately 0.3 vs. K(5')/K(S) approximately 0.6). As previous analyses have indicated, the inclusion of slightly deleterious polymorphism confounds the inference of positive selection. The analysis of published polymorphism data covering 97 verified 5' regions of Drosophila suggests more pronounced selective constraint on the 5' untranslated region and the core promoter (together corresponding to approximately 200 bp in this data set) when compared to the more distal portion of the 5' region of genes.     

Transpecific polymorphisms in an inversion linked esterase locus in Drosophila buzzatii

Nucleotide variation was studied in a 1.1 kb section of the coding region of an Esterase gene (Est-A) that maps in the center of the segments rearranged by polymorphic inversions in the cactophilic Drosophila buzzatii. We examine 30 homozygous second-chromosome lines differing in gene arrangement and three D. koepferae isofemale lines as outgroups. Our data show that Est-A is a highly polymorphic gene at both synonymous and replacement sites. Significant departures from homogeneity in the distribution of the ratio of silent polymorphism to divergence predicted by the neutral theory reveals a local excess of silent polymorphism. This is consistent with the presence of two apparent narrow peaks of elevated silent polymorphism surrounding nonconservative amino acid substitutions. These polymorphisms as well as others at synonymous and nonsynonymous sites are shared with D. koepferae. We suggest that the presence of shared nucleotide polymorphisms is probably due to interspecific gene flow and/or balancing selection acting on replacement variants and/or to a decreased probability of loss of ancestral polymorphisms caused by linkage to an adaptive inversion polymorphism. Recurrent mutation and persistence of neutral ancestral polymorphisms cannot, however, be ruled out. The analysis of the distribution of nucleotide variation among the three chromosomal arrangements sampled reveals that derived arrangements (J and JZ(3)) are less polymorphic than the ancestral ST, and that the widely distributed ST and J arrangements are genetically differentiated. However, a significant number of polymorphisms are shared between arrangements, suggesting frequent exchange either from gene conversion or from double crossovers in heterokaryotypes. Finally, our present results in combination with data of sequence variation at the breakpoints of inversion J suggest that this old gene arrangement has risen in frequency in relatively recent times.