LongLangley1993.pdf

Natural Selection and the Origin of jingwei, a Chimeric Processed Functional Gene in Drosophila

Author(s): Manyuan Long and Charles H. Langley

Source: Science , Apr. 2, 1993, New Series, Vol. 260, No. 5104 (Apr. 2, 1993), pp. 91-95

Published by: American Association for the Advancement of Science

Stable URL: https://www.jstor.org/stable/2881129

JSTOR is a not-for-profit service that helps scholars, researchers, and students discover, use, and build upon a wide range of content in a trusted digital archive. We use information technology and tools to increase productivity and facilitate new forms of scholarship. For more information about JSTOR, please contact [email protected]. Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at https://about.jstor.org/terms

American Association for the Advancement of Science is collaborating with JSTOR to digitize, preserve and extend access to Science

This content downloaded from �������������71.57.134.161 on Wed, 29 Dec 2021 11:13:06 UTC�������������

All use subject to https://about.jstor.org/terms

REPORTS

Phosphoinositides could serve as specific membrane targets that bind proteins re- quired for the formation of transport vesi- cles, such as the polypeptides of the adaptor complex that link clathrin to the cytoplas- mic tails of certain transmembrane receptor proteins (for example, the mannose-6-phos- phate receptor) (28). Vesicles from rat adi- pocytes that contain the glucose transporter also contain PI 4-kinase, which may regulate the transport (fusion) of these vesicles with the plasma membrane in response to insulin (29). In addition to its potential role in signaling cell proliferation, PI 3-kinase asso- ciated with receptor protein tyrosine kinases at the plasma membrane may also take part in the endocytosis and down-regulation (ly- sosomal degradation) of these receptors. In this way, the duration and the magnitude of the growth signal might be modulated. The association of Vps34p with the membrane appears to be mediated by the product of another VPS gene, VPS 15. The VPS15 gene encodes a membrane-associated protein ki- nase (Vps15) (30, 31) that can be chem- ically cross-linked to Vps34p (21). This raises the possibility that the Vps15 and Vps34 proteins may function together as components of a signal transduction com- plex that regulates intracellular protein sorting decisions.

REFERENCES AND NOTES

1. L. C. Cantley etal., Cell 64, 281 (1991). 2. J. M. Backer et al., EMBO J. 11, 3469 (1992). 3. S. Soltoff, S. Rabin, L. Cantley, D. Kaplan, J. Biol.

Chem. 267,17472 (1992). 4. M. J. Berridge and R. F. Irvine, Nature 341, 197

(1989). 5. M. Whitman et al., ibid. 332, 644 (1988). 6. D. L. Lips et al., J. Biol. Chem. 264, 8759 (1989). 7. L. A. Serunian et al., ibid., p. 17809. 8. C. Carpenter et al., ibid. 265, 19704 (1990). 9. F. Shibasaki, Y. Homma, T. Takenawa, ibid. 266,

8108 (1991). 10. M. Otsu et al., Cell 65, 91 (1991). 11. I. D. Hiles etal., ibid. 70, 419 (1992). 12. C. A. Koch, D. Anderson, M. F. Moran, C. Ellis, T.

Pawson, Science 252, 668 (1991). 13. P. Hu et al., Mol. Cell. Biol. 12, 981 (1992). 14. C. J. McGlade et al., ibid., p. 991. 15. P. K. Herman and S. D. Emr, ibid. 10, 6742 (1990). 16. S. Kornfeld and 1. Mellman, Annu. Rev. Cell Biol.

5, 483 (1989). 17. D. J. Klionsky, P. K. Herman, S. D. Emr, Microbiol.

Rev. 54, 266 (1990). 18. J. E. Rothman and L Orci, Nature 355, 409

(1992). 19. T. Stevens et al., Cell 30, 439 (1982). 20. T. R. Graham and S. D. Emr, J. Cell Biol. 1 14, 207

(1991). 21. J. H. Stack, P. K. Herman, P. V. Schu, S. D. Emr,

EMBO J., in press. 22. K. R. Auger, C. L. Carpenter, L. C. Cantley, L.

Varticovski, J. Biol. Chem. 264, 20181 (1989). 23. G. Endemann, S. N. Dunn, L. C. Cantley, Bio-

chemistry26, 6845 (1987). 24. Genetics Computer Group sequence analysis

package; University of Wisconsin. 25. S. K. Hanks et al., Science 241, 42 (1988). 26. D. R. Knighton et al., ibid. 253, 407 (1991). 27. M. P. Sheetz and S. J. Singer, Proc. Natl. Acad.

Sci. U.S.A. 71, 4457 (1974). 28. B. M. Pearse and M. S. Robinson, Annu. Rev. Cell

Biol. 6, 151 (1990).

29. R. L. Del Vecchio and P. F. Pilch, J. Biol. Chem. 266,13278 (1991).

30. P. K. Herman et al., Cell 64, 425 (1991). 31. P. K. Herman, J. H. Stack, S. D. Emr, EMBO J. 10,

4049 (1991). 32. T. Kunkel, Proc. Natl. Acad. Sci. U.S.A. 82, 5463

(1985). 33. F. Sherman, G. R. Fink, L. W. Lawrence, Methods in

Yeast Genetics: A Laboratory Manual (Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, 1979).

34. M. Whitman et al., Nature 315, 239 (1985). 35. J. P. Walsh, K. K. Caldwell, P. W. Majerus, Proc.

Natl. Acad. Sci. U.S.A. 88, 9184 (1991). 36. L. J. Wickerham, J. Bacteriol. 52, 293 (1946). 37. K. R. Auger et al., Cell 57, 167 (1989). 38. We thank members of the Emr lab for discussions

and G. Huyer' for his help in optimizing the Pi 3-kinase assay. Supported by the Howard Hughes Medical Institute, a grant from the NSF (to S.D.E.), and by a fellowship from the Deutsche Forschungs- gemeinschaft (to P.V.S.). S.D.E. is an investigator of the Howard Hughes Medical Institute.

8 October 1992; accepted 13 January 1993

Natural Selection and the Origin of jingwei, a Chimeric Processed Functional Gene in Drosophila

Manyuan Long* and Charles H. Langley

The origin of new genes includes both the initial molecular events and subsequent pop- ulation dynamics. A processed Drosophila alcohol dehydrogenase (Adh) gene, previously thought to be a pseudogene, provided an opportunity to examine the two phases of the origin of a new gene. The sequence of the processed Adh messenger RNA became part of a new functional gene by capturing several upstream exons and introns of an unrelated gene. This novel chimeric gene, jingwei, differs from its parent Adh gene in both its pattern of expression and rate of molecular evolution. Natural selection participated in the origin and subsequent evolution of this gene.

How genes with novel functions evolve remains a fundamental and fascinating ques- tion. Gene duplications (1), exon shuffling, and processed genes (2) have been suggested as important sources of novel genes, but little is known about the evolutionary mech- anisms or the participation of natural selec- tion in their early history. Our analysis of the structure, expression, and evolution of a putative processed Adh pseudogene in Dro- sophila provided an opportunity to examine both the early molecular events and the evolutionary processes that created it. The jingwei (jgw) gene, as we named it (3), is a locus located on chromosome 3 in the Dro- sophila sibling species D. teissieri and D. yakuba. A part of jgw was initially observed to hybridize to the Adh probe (4). Further analysis suggested that in a single event, in the ancestor of the two species the Adh portion of this potential pseudogene was retrotransposed from an mRNA of the Adh locus on chromosome 2 (5). To understand the molecular population genetics of this potential pseudogene, we investigated its within-species DNA sequence variation. However, our results convince us that jgw is not a pseudogene and that the Adh-derived sequence is a part of a novel gene.

The Adh-derived portion was sequenced from ten jgw alleles of D. teissieri and of 20 jgw alleles of D. yakuba collected from natural populations (6). Figure 1 shows the

Section of Genetics, Section of Evolution and Ecology, and Center for Population Biology, University of Cali- fornia at Davis, Davis, CA 95616.

*To whom correspondence should be addressed.

distribution of nucleotide polymorphisms within species and the variation between species in a 765-bp segment that corre- sponds to the protein coding region of the Adh gene. Only one possible ancestral poly- morphism was apparent (site 782). A sum- mary of the DNA sequence variation in this region is presented (Table 1).

Unexpectedly, most polymorphisms were silent (eight out of ten in D. teissieri and 19 out of 21 in D. yakuba). If jgw were a pseudo- gene in which mutations had no phenotypic effect, most changes would be at replacement sites and there would be a reasonable frequen- cy of stop codons (7). No new stop codons were present. Another prediction of the pseu- dogene hypothesis is that a pseudogene has a higher overall level of nucleotide variation. Estimates of DNA sequence polymorphism in jgw were similar to that found in the Adh gene (Table 1) and are typical of many functional genes in Drosophila (-0.005) (8). A compar- ison of between-species divergence also re- vealed a significant bias toward silent versus replacement substitutions, and the degree of bias was smaller than for the within-species comparison. The bias was noted in the earlier study (5) but by itself was not large enough to convince the authors that this was not a pseudogene. In addition, insertions and dele- tions are abundant in the evolution of mam- malian pseudogenes (7). In jgw, no length polymorphism or divergence was observed in the coding region (9). These molecular pop- ulation genetic observations imply that the Adh-derived sequence is all or part of jgw, a new functional gene in D. teissier and D. yakuba. This conclusion motivated our search

SCIENCE * VOL. 260 * 2 APRIL 1993 91

This content downloaded from �������������71.57.134.161 on Wed, 29 Dec 2021 11:13:06 UTC�������������

All use subject to https://about.jstor.org/terms

for an RNA from jgw and a more detailed analysis of its structure.

Northem (RNA) analysis (10) of both total and polyadenylate [poly(A)J-selected RNAs from adults probed with jgw yielded only a single band that corresponded in size to that derived from the Adh gene (11). Because of the sequence similarity between the two genes and the abundance of the Adh mRNA, this result was not surprising. Two interpretations are possible. (i) The jgw

mRNA may be detectably abundant but obscured by the Adh signal if its size is similar to that of Adh. (ii) The jgw mRNA may be so rare (or nonexistent) that it was undetect- able by the Northem technique.

We used the polymerase chain reaction (PCR) in conjunction with reverse transcrip- tion (RT-PCR) (12) to identify RNAs with greater sensitivity and specificity than is pos- sible with Northem analysis. Several regions can be found in which Adh and the Adh-

derived portion of jgw differ by two or three contiguous substitutions. Under optimal PCR conditions, primers that contain these substi- tutions at their 3' ends specifically amplified jgw RNA. A jgw-specific amplification prod- uct was obtained from total RNA extracted from both species. Sequencing these products confirmed the characteristic substitutions of the jgw gene in the 571-bp amplified fragment from D. teissieri and the 321-bp amplified fragment from D. yakuba (13). These results demonstrated the presence of jgw-derived RNA in both species.

The single band in the 5' rapid amplifica- tion cDNA end (RACE) product and the identity of the two sequenced clones from this product indicated that jgw RNA has a distinct 5' end in D. teissieri (14). In both species, jgw appears to have captured more than 180 bp of additional exon (or exons) in the 5' direction (Fig. 2). The length of the 3' end of D. teissieri jgw RNA is similar to that of Adh RNA. Figure 2 also shows the sequence of these products. By extending the Adh reading frame in the 5' direction, the D. teissien jgw appears to have acquired a new start codon and potentially encodes a hybrid protein with an additional 77 amino acids added to the Adh- derived protein. A preliminary analysis of jgw in D. yakuba indicated that it contains a similar structure.

To determine the structure of the genomic region (or regions) containing the exon (or

Fig. 1. Polymorphic D. teis4eri D. yakuba and divergent differ- 29 a 1 2 3 4 56 7 9 90 1 23 45 6 7 8 9 3 0 1 1 1 2 1 3 1 4 1 5 16 ii 7 s 1 9 210 ences in the Adh-de- 2 2S A.T.A

231 T.A

rived portion of jgw. 232 C.G 2145 a.A.T . . – C

The nucleotide posi- 246 C.A 254 0.A

tions are indicated in 297 C.G 323..G.A the first column, num- 35 C.T T C

365.. . .

369 a.A bered according t ht 37

372 C.A reported from D.yaku- 37 A

387 T . . . . ba (5). The second and 388 a . . .C

389 C . . .T thirteenth columns show 390 a . . .A the consensus nucleo- 438 A .. .

450 A . . .T tides. The dots indicate 451 C . . . .

459 a . . .T that the nucleotides are 460 a c : : C identical to those of the 463 A . . . .

474 A . .. consensus. The first row 480 a . . .A

2 C . . .C. . . contains the numbers of 497 C . . .T

S00 T . . .C the sequenced alleles 525 C . . . .

526 a . . .C (1 0 from D. teissieri and 527 C . . .A 535 C …C .T

557 a Ak .' . A . 20 from D.yakuba).56 AAA

569 C.a So81 .A 599 a.a .. . .A A A AA 608 a.C 61 4 a.T 650 C.C . … .. T 680 C a .. C 719 a . . . . . A a .(Gc 72.8 C . . . . . C . …. .G. c 743 a. . . . . . . . . 749 T. . . .a a a T 776 C .. C 779 C. . . . . C . .. .. . T 782 A . . . . C . C . . . A . . . T

785 T. . . . . . . 812 T. . . . . . T. . . . . 863 a. . . . . . . T. . . . T T . T T T . . . T 878 a. . . . . . . A A'. A A A A A A . . . . A A 884 C. . . . . . T 887 a. . . . . . . A 'A 'A 'A 'A 'A ' . A 'A .A 893 a. . . . . . c 896 C. . . . . . C T .. . T' T . T T 917 A. . . . . . C 941 C. . . . . . A 951 A a i a a A 960 T. . . .. A 961 G. . . . . . C 974 C. . . . . . T 977 a T. . . . . .

Table 1. Summary of the variation within and between species. The numbers of segregating (or polymorphic) sites are classified as replacement (R) and silent (S). Estimates of 0 and vr are measures of within-population variation in terms of the parameter 4NpL per nucleotide, assuming the neutral theory of molecular evolution (20) where N is the effective population size and , is the mutation rate for selectively equivalent nucleotides. The values of divergence between species (average divergence per nucleotide) are the results of averaging pairwise comparisons between the alleles from the two species minus average within-species differences, subject to the multiple substitution correction according to the Jukes and Cantor model (28). The standard deviations of 0 and divergence estimates were calculated as described (29). Adh polymorphism data (19 of all the samples) of D. yakuba and D teissieri Adh sequences are from (5, 19).

Segregating sites (n) Gene 0 Ir

R S

Polymorphism within species

D. teissieri jgw 2 8 0.005 + 0.003 0.005 D. yakuba jgw 2 19 0.008 ? 0.003 0.006 D. yakuba Adh 0 18 0.006 ? 0.003 0.006

Divergence between species

jgw 0.041 ? 0.009 0.127 t 0.027 Adh 0.007 ? 0.004 0.059 ? 0.018

92 SCIENCE * VOL. 260 * 2 APRIL 1993

This content downloaded from �������������71.57.134.161 on Wed, 29 Dec 2021 11:13:06 UTC�������������

All use subject to https://about.jstor.org/terms

. REPORTS

exons) from which this new mRNA is de- rived, we carried out PCR-based genomic sequencing. Three introns (and three exons) were found 5' to the Adh-derived exon in both species (Fig. 2). The size and positions of the three introns are similar in the two spe- cies, and the standard intron-splicing signal (GT-AG) was found in each case except for the D. teissieri second intron (GT-AT).

Because the structure of jgw is similar in both species, it seems likely that the insert- ed exons of Adh recruited the observed three 5' exons (and introns) rather quickly. This interpretation is supported by the ob- servation that no silent substitution accu- mulated in the Adh-derived portion of jgw before speciation, although several replace- ment substitutions were common to jgw in both species. This suggests that there was relatively little time between the insertion event and the subsequent molecular evolu- tionary events: the recruitment of the addi- tional 5' components into a functional transcription unit, the substitution of sev- eral amino acids, and the speciation.

In spite of similarities of structure and early common ancestry, jgw appears to have evolved distinct pattems of expression in the two species. Figure 3 shows the pattem of jgw expression at different development stages de- termined by quantitative PCR technique (15). Quantitative PCR provides at least a relative comparison of abundance and con- firmed that the amount of jgw RNA was generally smaller than the amount of Adh. In D. teissieni, jgw showed sex-specific expression. The expected band was abundant in amplifi- cations of RNA from adult males, whereas no band was present in reactions from RNA of other stages or adult females. However, in D. yakuba reactions with RNA from all stages showed less expression than was seen with D. teissieri RNA. Larval stages Ll and L2 and

adult male RNA showed more expression than L3 and adult female RNA. Just as the protein sequences evolved rapidly, the se- quences controlling the regulation of jgw have also evolved rapidly. Further analysis will determine if the species-specific patterns of expression are the result of the divergence of sequences in and around the jgw transcription unit or whether these expression patterns reflect divergence at more distant, trans-act- ing regulatory loci.

Whereas a distinct mRNA was expressed from jgw in each species, the conclusion that jgw is a functional gene also rests on the evolutionary interpretation of the polymor- phism observed within each species. The extension of this evolutionary analysis in order to incorporate between-species diver- gence reveals that the early history of jgw was dominated by positive natural selection instead of the neutral mutation and genetic drift that are expected to characterize the dynamics of pseudogenes (16).

To explore the divergence of jgw in the period immediately after the Adh retro-

transposition but before the split of D. teissieri and D. yakuba, we compared the Adh-derived portions of all jgw alleles with the available Adh alleles in D. teissieri and D. yakuba and with Adh alleles from various members of the melanogaster subgroup (Fig. 4). The fixed nucleotide substitutions that distinguished the Adh-derived portion of all jgw alleles from Adh in all alleles of D. teissieri and D. yakuba were assumed to have been fixed in the ancestral species after the insertion. Eight such sites exist at positions 339, 576, 600, 791, 831, 861, 918, and 954. Surprisingly, all eight substitutions were amino acid replacements; no silent substitutions were found. By incorporating the available Adh sequences of the other six members of the melanogaster subgroup, we found that all eight substitutions were unique mutations that were fixed in the jgw lineage. Assuming the replacement rate of substitution of 0.5 x 10-9 per nucleotide per year [calculated from melanogaster sub- group Adh data (17)1, these eight substitu- tions (out of 572 possible replacement sites)

Fig. 2. Structure of the A lingwei B D. telssierl5' portion sequence of jgw jgw transcription unit 4-Captured region-* .4 Retrotransposed region . AG (13). (A) The structure p7 p P CTGTATATTAACCAATTT"TTCATAAzCGCrCGTCTTACAACCATrACrCTTCTAAAACGCACTCCTATTACTCGTGCCAAAGg of jgw deduced from P2 Exon the comparison of the tgagtctacaaaccctattagatcctgttcatatggtaattcattattaacaatgcgattcg cagAGCAGCACCGACAACACGCACT

sequenced D. teissieri AT : TMAMTMM CTACTACTAGoAhgtaagtaactcgtgggcaaaagtccaagtcagaattgcttatCtgCaaacccatoTCCTTTTCCAGooAxTCTGGTACC

RNA and genomic Intron 2 DNA. The four boxes T C GCATCATCCACCATGGGC.GGg taagtactatatgatgagttgtctagaxctcca

show the four exons in jgw, Intron and the solid and aaactcaagacttaccacctttgtcataaacctataaaaggattataaaagatatatatatatatattatttttcagTCCGGTTGTCAGTTG

hatched portions of the ATG GTG boxes are the putative Adh Exon4 (AAd-erive protein coding regions, Adh-derived and recruited, respectively. The open boxes are the putative D yakba 5 portion sequence of jgw

. …?TCTTTCGACCA GCTTAACCCcT TGTTCTGAAAAGgtgaggctaaaaatccctattagattctgttc untranslated portions. The lines among boxes are the three introns. The E I InTrAM I exon structure of the original Adh gene is also given, based on (5). The atatggtaattcactattaacaatgccattcgcaagTAAAAGCAGCACAGCCAACACGCACTCTACTACTCTTACTACAACAAgtaagtaa gray portions in four Adh exons are fused into one exon in the jgwgene. (B) Exon 2 Sequences of the 5' regions of jgwfrom D. teissieri and D. yakuba. Introns ctcgtgggcaaaagtccaagtcaggattgetaagtgcaaaccatgtccttttacagGGAATCTGGTTCCGTTTTTAGGTST" are shown in lowercase letters. The codon ATG underlined in the D. osctc&Txcc GGcGGG gtagtagtctagaactccaaaacttaccacatttgatataagcgattcaaatgaga

teissieri sequence is the hypothetical start codon. The 3' region of D. caaatacataactatccaggagatctgetccgggttctcttctcttcatatcgatcatattaatattgtctcagtgcagttgtcagttgcag

teissieri jgw is the same as that of the published sequence (5). An "n'3 ttcagcgagtaattgcatctctttaaatttacttatttgatcaaatcgaca gAGCCAGTGCGACCAAAAGGCGTTTT CC….

apparent polyadenylation signal, AATAAA, is at positions 1 123 to 1 128 of Exon 4 (AdhOned) D. teissieri.

Fig. 3. Developmental pattern of jgw expression. D. teissien The same method as described (12) was used Adult to amplify jgw RNA from total RNA extracted Li L2 L3 Female Male from the different development stages. The cDNA product from 300 ng of RNA was added

as a template for amplification because for the 571 range from 100 ng to 1000 ng of total RNA, the yield of jgw RNA PCR products was observed to be linearly related to the template RNA concen- tration (11). PCR cycles (30) were conducted at D. yakuba 940C (1 min) for denaturation, 600C (2 min) for annealing, and 720C (3 min) for polymerization.

The products of the PCR reactions were electro- 321 phoresed on a 1% agarose gel and stained with ethidium bromide. More Adh RNA was amplified from the same cDNA preparations, which sug-

gests that Adh is expressed in larger amounts than jgw (1 1). Molecular size markers are on the right in base pairs.

SCIENCE * VOL. 260 * 2 APRIL 1993 93

This content downloaded from �������������71.57.134.161 on Wed, 29 Dec 2021 11:13:06 UTC�������������

All use subject to https://about.jstor.org/terms

yielded an estimate of about 30 million years for the time of the last common ancestor of jgw and Adh. This conflicts with the age of the melanogaster subgroup, which is estimated to be 17 to 20 million years (18). The lack of silent substitutions im- plies that the insertion of jgw occurred close to the time of the speciation of D. teissieri and D. yakuba. The excess of replacement substitutions in the jgw line of descent between the time of its insertion and the speciation of D. teissieri and D. yakuba is consistent with the model that jgw respond- ed to positive natural selection and evolved a new function.

This analysis of the early history of jgw leads us to examine its dynamics after the divergence of the two species. One method to explore the role of natural selection at a single locus is to compare the replacement and silent substitutions found between species with the replacement and silent polymorphism found within species (19). The variation in the Adh-derived portion of jgw within and be- tween D. teissieri and D. yakuba is summarized (Fig. 4). A relative excess of (fixed) replace- ment substitutions over (fixed) silent substitu- tions between species (21:16) is apparent when compared to the proportion of replace- ment polymorphisms over silent polymor- phisms within a species (4:27). These results are inconsistent with the null hypothesis in the test proposed (19), which incorporates the assumption of the selective neutrality (X2 = 13.6, P < 0.001). Therefore, adaptive protein evolution remained an important force in the evolutionary history of jgw after the separation of the two species. Further evidence supporting this view of the diver- gence of jgw in D. teissieri and D. yakuba comes from a comparison of the captured 5' coding regions. Eleven out of the 15 differ- ences in 153 alignable coding sites cause amino acid replacements.

Neutral allele theory considers polymor-

phism to be a transient phase of molecular evolution (20). Under the assumptions of this theory, the level of the between-species divergence at different loci is positively cor- related with the level of within-species poly- morphism. This idea is embodied in a sim- ple, two-by-two statistical test (HKA test) (21) that calculates the expected amounts of polymorphism and divergence at two loci and can indicate whether the observed amounts are consistent with the neutral theory. Because we determined that there is an excess of replacement substitutions, we directed our attention to the silent variation in jgw and Adh. Figure 4 shows the observed numbers of silent substitutions and polymor- phisms in the Adh-derived portion of jgw. Whereas their within-species polymorphism is comparable, the between-species diver- gence of jgw is twofold greater than that of Adh (X2 = 7.12, P < 0.01).

The Adh gene in these two species has a very biased codon usage with a 84 to 87% G+C content at the third codon positions. Furthermore, the codon preference as shown by the within-species polymorphism data in D. yakuba (19) also shows a very high G+C content (85.7%). Nevertheless, the codon usage in the Adh-derived portion of jgw seems to be evolving to a smaller amount of G+C. The G+C contents at the third codon position summed over all silent poly- morphism sites in the two species are 70.0% for D. teissieri and 62.4% for D. yakuba, respectively. The G+C content at the third codon position in the fixed replacement sites between the two species is 60.0%. This evolution of a moderate G+C content of third sites is consistent with the higher rate of silent divergence for jgw and the general observation of a negative correlation be- tween the codon bias (and G+C content) and the rate of silent divergence (22).

The results presented reveal the molecu- lar, genetic, and evolutionary mechanisms

that participated in the origin of the chimer- ic processed functional gene jgw. The chi- meric nature of the expressed product of this novel gene is unique among processed genes, including those few that are functional, such as the second gene for rat insulin, the human PGK gene, and other examples that have been interpreted as outcomes of retroposi- tions (2, 23). It has been proposed that retrotransposition can be a major cause of intron loss during evolution (24), as it was with jgw. The evolution of jgw also demon- strates that retrotransposition can be a source of new, intron-containing genes in eukaryotic evolution. It is also unique among instances of exon-shuffling (25) in that it clearly involved the retrotransposi- tion of an mRNA. Experiments of cross- hybridization with a genomic Southem (DNA) blot and Northern analysis with the use of a probe generated from the captured portion of jgw indicated that the source of the captured portion of jgw was a duplication of an unrelated gene (26). The analysis of the early molecular evolution of jgw indicat- ed that jgw was functional from the begin- ning and experienced strong adaptive evolu- tion for what must have been a novel func- tion. Prevailing theories of the origin of new genes assume initial relaxation of selection (1). Our results provide a contrary interpre- tation in which natural selection is present throughout the origin of a new gene.

Fig. 4. Phylogenetic .,–

analysis of the evolution , 0/0 retrotransposition 8/0 of the Adh-derived por- tion of jingwei. The topol- / ogy reflects the as- / z/ , sumed evolutionary rela- ogy reflects the as- ~~~~~4/11 21/16 tion between jingwei and 11. 23.0 Adh (5). The small (stip-

pled) trees at the tips in- 0/18 18 8 J -2/8 2/19-{ }19 dicate that the detailed Adh Adh Adh jingwei jingwei phylogenies of the al- outgroup teissieri yakuba teissieri yakuba leles within species are

unknown. Numbers presented as fractions are the numbers of (fixed or polymorphic) amino acid replacement differences over the numbers of (fixed or polymorphic) silent differences. The Adh outgroup includes the sequences of D. melanogaster, D. simulans, D. sechellia, D. mauritiana, D. erecta, and D. Orena. These outgroup sequences were used to confirm that all functional Adh genes have the same DNA sequence for the eight codons that are distinct to jingwei (in both D. teissieri and D. yakuba) with one exception. At the eighth codon (nucleotide position 954) of a. orena Adh, there is a substitution that conferred an amino acid replacement. The italic numbers are the average numbers of silent substitutions and the numbers of silent polymorphisms under the HKA test, which also incorporates the standard assumptions of the neutral theory.

REFERENCES AND NOTES

1. S. Ohno, Evolution by Gene Duplication (Spring- er-Verlag, Berlin, 1970); T. Ohta, Genome 31, 304 (1989); J1 B. S. Haldane, The Causes of Evolution (Longmans, New York, 1932).

2. W. Gilbert, Nature 271, 501 (1978); G. F. Hollis, P. A. Hieter, 0. W. McBride, D. Swan, P. Leder, ibid. 296, 321 (1982); J. R. McCarrey and K. Thomas, ibid. 326, 501 (1987).

3. In an ancient legend from China (San Hai Jing), Jingwei, a daughter of the Emperor Yande (first Chinese emperor 3000 B.C.), tragically drowned while swimming in the East China Sea. Jingwei was then reincarnated as a beautiful bird that drops stones and wood into the sea in an attempt to fill it, thus preventing others from drowning. We used the name 'jingwei" because this gene avoided the usual fate of processed gene (death) and was "reincarnated" into a new structure with novel function.

4. C. H. Langley et al., Proc. Natl. Acad. Sci. U.S.A. 79, 5631 (1982).

5. P. Jeffs and M. Ashburner, Proc. R. Soc. London, Ser. B 244, 151 (1991).

6. The sequences were completely determined from both strands by the dideoxyribonucleotide chain- terminating method [F. Sanger et al., Proc. Natl. Acad. Sci. U.S.A. 74, 5463 (1977)] on single-strand- ed DNA templates prepared directly from the PCR products [R. G. Higuchi and H. Ochman, Nucleic Acids Res. 17, 5865 (1989); K. R. Saiki et al., Science 239, 487 (1988)]. The primers used in PCR are 5'-GACAGTGATATGAGATTGCCG-3' and 5'- GGAAGAATGTGAGTGTGCTTCG-3' for D. teissieri and 5'-CTCTAGAACACCAAAACTTACC-3' and 5'- GCATACATATTTGTAAAATGCAAG-3' for D. yaku- ba. Each sequence was from a single fly of a separate isofemale line that had been mated to its siblings for three generations. The two alleles or haplotypes in the five heterozygous flies were deter-

94 SCIENCE * VOL. 260 * 2 APRIL 1993

This content downloaded from �������������71.57.134.161 on Wed, 29 Dec 2021 11:13:06 UTC�������������

All use subject to https://about.jstor.org/terms

REPORTS

mined by the cloning and subsequent sequencing of the PCR products.

7. W.-H. Li, in Evolution of Duplicate Genes and Pseudogenes, M. Nei and R. N. Koehn, Eds. (Sinauer, Sunderland, MA, 1983), pp. 14-37.

8. C. H. Langley, in Population Biology of Genes and Molecules, N. T. Takahata and J. F. Crow, Eds. (Baifukan, Tokyo, 1990), pp. 75-91; M. Kreitman, in Evolution at the Molecular Level, R. K. Selander, A. G. Clark, T. S. Whittam, Eds. (Sinauer, Sunder- land, MA, 1991), pp. 204-221.

9. Only two length polymorphisms are observed in the 3' flanking region of D. yakuba.

10. J. C. Alwine etal., Proc. Natl. Acad. Sci. U.S.A. 74, 5350 (1977).

11. M. Long and C. H. Langley, unpublished data. 12. E. S. Kawasaki, in PCR Protocols, M. A. Innis et al.,

Eds. (Academic Press, New York, 1990), pp. 21-27. 13. Total RNA was extracted from male adult flies (D.

teissierl) and Li larvae (D. yakuba) as described [D. A. Goldberg et al., Cell 34, 59 (1982)]. We treated the RNA with ribonuclease-free deoxyribonuclease to avoid potential DNA contamination during the PCR experiment. For RT-PCR, the jgw-transcribed RNA sequences were amplified from cDNA as de-

scribed by Kawasaki (12) with the jgw-specific primers. The primers included the pair 5'-GGCAC- TCAATCCAAAGGTGTG-3' (P1) and 5'-GCCCA- AGTCCAGTTTCCAGAGT-3' (P2) for D. teissieri and the pair 5'-TCCTTGAGCAACAAGAACGTAA-3' (P3) and 5'-GTAGTTGACGGCGATGGTTGC-3' (P4) for D. yakuba. The PCR products were cloned into Bluescript SK+ and sequenced by the dideoxy- ribonucleotide chain-terminating method.

14. Drosophila teissieri 5' and 3' ends of RNA were amplified from cDNA with the RACE technique (27). For the 5' end of jgw RNA from D. teissieri, the primer used for primer extension was primer 5 (P5, 5'-GGTGGTCTCGGCAATGG-3'), which anneals to the RNA from both jgwand Adh. The cDNA gener- ated afterwards was tailed with deoxyadenosine triphosphate. Three primers were used for the sub- sequent PCR reactions. They included two artificial primers (27) [dT1 7-adapter and adapter primers that are able to anneal to the artificially added poly(A) tail of cDNA] and jgw-specific primer 6 (P6, 5'-TCACATCGTAGGGGTAGAAGGTGACGCA-3'). For the 3' end of jgw RNA in D. teissieri, the jgw-specific primer 7 (P7, 5'-CCAATTTGATTATAA- TGGCGCTTCG-3') or P1 in conjunction with the dT17-adapter primer was used to amplify the 3' portion of jgw RNA of D. teissieri. The 5' end of D. yakuba jgw RNA was amplified from cDNA with P7 and D. yakuba jgw-specific primer P4. Genomic DNA was amplified with the use of P7 and a down- stream primer; P6 was used for D. teissieri and primer 8 (P8, 5'-TCCAGACCAATGCCTCCCAGAC- CGGCAACGAAAATT-3') for D. yakuba. The posi- tions of the introns are determined by alignment of genomic and RNA sequences. Sequencing meth- ods of the PCR products included the direct se- quencing and the cloning-sequencing approaches as described (6).

15. J. Chelly, J.-C. Kaplan, P. Maire, S. Gautron, A. Kahn, Nature 333, 858 (1988); J. Singer-Sam et al., Nucleic Acids Res. 18, 1255 (1990).

16. W.-H. Li, T. Gojobori, M. Nei, Nature 292, 237 (1981).

17. D. T. Sullivan et al., Evol. Biol. 24, 107 (1990). 18. D. Lachaise et al., ibid. 22, 159 (1988). 19. J. H. McDonald and M. Kreitman, Nature 351, 652

(1991). 20. M. Kimura, The Neutral Theory of Molecular Evolu-

tion (Cambridge Univ. Press, Cambridge, 1983); G. A. Watterson, Theor. Popul. Biol. 7, 256 (1975); M. Nei and F. Tajima, Genetics 97,145 (1981).

21. R. R. Hudson etal., Genetics 116, 153 (1987). 22. E. N. Moriyama and T. Gojibori, ibid. 130, 855

(1992); P. Sharp and W.-H. Li, J. Mol. Evol. 28, 398 (1989).

23. M. B. Soares et al., Mol. Cell. Biol. 5, 2090 (1985); A. M. Weiner et al., Annu. Rev. Biochem. 55, 631

(1986); H.-H. M. Dahl et al., Genomics 8, 225 (1990); A. Ashworth et al., EMBO J. 9, 1529 (1990); J. Brosius, Science 251, 753 (1991 ).

24. G. R. Fink, Cell 49, 5 (1987).

25. R. L. Dorit, L. Schoenbach, W. Gilbert, Science 250, 1377 (1990); L. Banyai et al., FEBS Lett. 163, 37 (1983).

26. We named this non-Adh parent gene of jgw as Yande (ynd) in consistency with the Chinese

legend about jingwei (3, 11). 27. M. A. Frohman et al., Proc. Natl. Acad. Sci. U.S.A.

85, 8998 (1988). 28. M. Nei, Molecular Evolutionary Genetics (Colum-

bia Univ. Press, New York, 1987), pp. 276-279; T. H. Jukes and C. R. Cantor, in Mammalian Protein Metabolism, H. N. Munro, Ed. (Academic Press, New York, 1969), pp. 21-132.

29. W.-H. Li and D. Graur, Fundamentals of Molec- ular Evolution (Sinauer, Sunderland, MA, 1991); R. R. Hudson, in Oxford Surveys in Evolutionary Biology, D. Futuyma and J. Antonovics, Eds.

(Oxford Univ. Press, New York, 1990), vol. 7, pp. 1-40.

30. We thank J. Gillespie, A. Clark, M. Turelli, K. Burtis, and all members of our laboratory for helpful discussions and critical reading of the manuscript. D. Lachaise kindly provided the D. teissieri isofemale lines. We also thank P. Jeffs and M. Ashburner for sharing DNA sequences of

the Adh-derived portions of jgw in D. yakuba and D. teissieri before they published these se- quences (5). Supported by NSF (C.H.L.) and the Center for Population Biology of the University of California at Davis, a University of California at Davis Graduate Research Award, and a Jastro- Shields Research Scholarship (M.L.).

26 October 1992; accepted 4 February 1993

Modulation of Neuronal Migration by NMDA Receptors

Hitoshi Komuro and Pasko Rakic

The N-methyl-D-aspartate (NMDA) subtype of the glutamate receptor is essential for neuronal differentiation and establishment or elimination of synapses in a developing brain. The activity of the NMDA receptor has now been shown to also regulate the migration of granule cells in slice preparations of the developing mouse cerebellum. First, blockade of NMDA receptors by specific antagonists resulted in the curtailment of cell migration. Second, enhancement of NMDA receptor activity by the removal of magnesium or by the application of glycine increased the rate of cell movement. Third, increase of endogenous extracellular glutamate by inhibition of its uptake accelerated the rate of cell migration. These results suggest that NMDA receptors may play an early role in the regulation of calcium-dependent cell migration before neurons reach their targets and form synaptic contacts.

In the developing brain, most immature neurons migrate to their distant final desti- nations by extending their leading processes and translocating their soma through a ter- rain that is densely packed with previously generated neurons and their processes (1). This movement of immature neurons is es- sential for the establishment of normal cy- toarchitecture, synaptic connectivity, and function in the brain (2). In the cerebellum, granule cells migrate from the site of their origin in the germinal external granular layer toward the internal granular layer along the elongated processes of Bergmann glial cells (3). Recently Komuro and Rakic have dem- onstrated that the rate of granule cell move- ment across the molecular layer in the cere- bellum depends both on extracellular Ca2+ concentrations and on Ca2+ influx through N-type Ca2+ channels (4). However, the regulatory mechanism underlying this Ca2+- dependent process remains unknown. We have now examined the role of ionotropic receptors-NMDA, non-NMDA, GABAA, and GABAB (GABA is y-aminobutyric ac- id)-in granule cell migration; these recep- tors are expressed by immature granule cells (5) and can directly and indirectly affect

Section of Neurobiology, Yale University School of Medicine, New Haven, CT 06510.

Ca2+ influx and intracellular Ca2+ concen- trations (5, 6).

To examine whether NMDA, non- NMDA, GABAA, and GABAB receptors play a significant role in the migration of granule cells, we used slice preparations of the developing mouse cerebellum stained with a lipophilic carbocyanine dye [1,1'- dioctadecyl-3,3,3 ',3 '-tetramethylindocar- bocyanine perchlorate (DiI)] and a laser scanning confocal microscope (7). Postmi- totic granule cells in slice preparations mi- grate from the external granular layer to- ward the internal granular layer (Fig. 1) (8). Antagonists to these receptors were added to the culture medium in separate experiments. Blockade of the non-NMDA subtype of glutamate receptors [that is, kainate and AMPA (L-a-amino-3-hydroxy- 5-methyl-4-isoxazole propionic acid) recep- tors] by 6-cyano-7-nitroquinoxaline-2,3-di- one (CNQX) (9), GABAA receptors by bicuculline (6), and GABAB receptors by phaclofen (10) failed to alter the rate of cell migration (Fig. 2, A and B). However, blockade of the NMDA subtype of gluta- mate receptor by D-2-amino-5-phosphono- pentanoic acid (D-AP5) (11) or (+)-5- methyl- 10, 1 1-dihydro-5H-dibenzo[a,d]cy- clohepten-5, 10-imine hydrogen maleate (MK-801) (12) significantly decreased the

SCIENCE * VOL. 260 * 2 APRIL 1993 95

This content downloaded from �������������71.57.134.161 on Wed, 29 Dec 2021 11:13:06 UTC�������������

All use subject to https://about.jstor.org/terms