Molecular evolution of Drosophila and higher Diptera. I. Micro-complement fixation studies of a larval hemolymph protein

EVOLUTION OF MULTIPLE KINDS OF FEMALE SPERM-STORAGE ORGANS IN DROSOPHILA

Females of all species belonging to the family Drosophilidae have two kinds of sperm-storage organs: paired spherical spermathecae and a single elongate tubular seminal receptacle. We examined 113 species belonging to the genus Drosophila and closely allied genera and describe variation in female sperm-storage organ use and morphology. The macroevolutionary pattern of organ dysfunction and morphological divergence suggests that ancestrally both kinds of organs stored sperm. Loss of use of the spermathecae has evolved at least 13 times; evolutionary regain of spermathecal function has rarely if ever occurred. Loss of use of the seminal receptacle has likely occurred only once; in this case, all descendant species possess unusually elaborate spermathecae. Data further indicate that the seminal receptacle is the primary sperm-storage organ in Drosophila. This organ exhibits a pattern of strong correlated evolution with the length of sperm. The evolution of multiple kinds of female sperm-storage organs and the rapidly divergent and correlated evolution of sperm and female reproductive tract morphology are discussed.

Key enzymes and proteins of crop insects as candidate for RNAi based gene silencing

RNA interference (RNAi) is a mechanism of homology dependent gene silencing present in plants and animals. It operates through 21-24 nucleotides small RNAs which are processed through a set of core enzymatic machinery that involves Dicer and Argonaute proteins. In recent past, the technology has been well appreciated toward the control of plant pathogens and insects through suppression of key genes/proteins of infecting organisms. The genes encoding key enzymes/proteins with the great potential for developing an effective insect control by RNAi approach are actylcholinesterase, cytochrome P450 enzymes, amino peptidase N, allatostatin, allatotropin, tryptophan oxygenase, arginine kinase, vacuolar ATPase, chitin synthase, glutathione-S-transferase, catalase, trehalose phosphate synthase, vitellogenin, hydroxy-3-methylglutaryl coenzyme A reductase, and hormone receptor genes. Through various studies, it is demonstrated that RNAi is a reliable molecular tool which offers great promises in meeting the challenges imposed by crop insects with careful selection of key enzymes/proteins. Utilization of RNAi tool to target some of these key proteins of crop insects through various approaches is described here. The major challenges of RNAi based insect control such as identifying potential targets, delivery methods of silencing trigger, off target effects, and complexity of insect biology are very well illustrated. Further, required efforts to address these challenges are also discussed.

Perspective: female remating, operational sex ratio, and the arena of sexual selection in Drosophila species

As commonly observed among closely related species within a variety of taxa, Drosophila species differ considerably in whether they exhibit sexual dimorphism in coloration or morphology. Those Drosophila species in which male external sexual characters are minimal or absent tend, instead, to have exaggerated ejaculate traits such as sperm gigantism or seminal nutrient donations. Underlying explanations for the interspecific differences in the presence of external morphological sexual dimorphism versus exaggerated ejaculate traits are addressed here by examining the opportunity for sexual selection on males to occur before versus after mating in 21 species of Drosophila. Female remating frequency, an important component of the operational sex ratio, differs widely among Drosophila species and appears to dictate whether the arena of sexual selection is prior to, as opposed to after, copulation. Infrequent female mating results in fewer mating opportunities for males and thus stronger competition for receptive females that favors the evolution of male characters that maximize mating success. On the other hand, rapid female remating results in overlapping ejaculates in the female reproductive tract, such that ejaculate traits which enhance fertilization success are favored. The strong association between female remating frequency in a given species and the presence of sexually selected external versus internal male characters indicates that the relationship be examined in other taxa as well.

Paleo-demography of the Drosophila melanogaster subgroup: application of the maximum likelihood method

The species divergence times and demographic histories of Drosophila melanogaster and its three sibling species, D. mauritiana, D. simulans, and D. yakuba, were investigated using a maximum likelihood (ML) method. Thirty-nine orthologous loci for these four species were retrieved from DDBJ/EMBL/GenBank database. Both autosomal and X-linked loci were used in this study. A significant degree of rate heterogeneity across loci was observed for each pair of species. Most loci have the GC content greater than 50% at the third codon position. The codon usage bias in Drosophila loci is considered to result in the high GC content and the heterogenous rates across loci. The chi-square, G, and Fisher's exact tests indicated that data sets with 11, 23, and 9 pairs of DNA sequences for the comparison of D. melanogaster with D. mauritiana, D. simulans, and D. yakuba, respectively, retain homogeneous rates across loci. We applied the ML method to these data sets to estimate the DNA sequence divergences before and after speciation of each species pair along with their standard deviations. Using 1.6 x 10(-8) as the rate of nucleotide substitutions per silent site per year, our results indicate that the D. melanogaster lineage split from D. yakuba approximately 5.1 +/- 0.8 million years ago (mya), D. mauritiana 2.7 +/- 0.4 mya, and D. simulans 2.3 +/- 0.3 mya. It implies that D. melanogaster became distinct from D. mauritiana and D. simulans at approximately the same time and from D. yakuba no earlier than 10 mya. The effective ancestral population size of D. melanogaster appears to be stable over evolutionary time. Assuming 10 generations per year for Drosophila, the effective population size in the ancestral lineage immediately prior to the time of species divergence is approximately 3 x 10(6), which is close to that estimated for the extant D. melanogaster population. The D. melanogaster did not encounter any obvious bottleneck during the past 10 million years.

Sequence, structure and evolution of the ecdysone-inducible Lsp-2 gene of Drosophila melanogaster

The Lsp-2 gene encodes a major larval serum protein (hexamerin) of Drosophila melanogaster. Transcription of Lsp-2 is controlled by 20-hydroxyecdysone. Here we report the analysis of the structure of the Lsp-2 gene including the adjacent 5' and 3' sequences. In contrast to all other known hexamerin genes, Lsp-2 does not contain an intron. The Lsp-2 mRNA measures 2312 bases, as deduced from experimental determination of the transcription-start and stop sites and conceptual translation results in a 718 amino acid hexamerin subunit, including a 21-amino-acid signal peptide. While the calculated molecular mass of the native 697-amino-acid subunit is 83.5 kDa, mass spectrometry gave a value of 74.5 kDa. We detected in the Lsp-2 gene a 2052-bp antisense ORF that probably does not code for any protein. An unusual accumulation of rarely used codon triplets was found at the 5' and 3' ends of the Lsp-2 ORF. The calculated secondary structure matches well with that of arthropod hemocyanins. Electron micrographs show for LSP-2 hexamers a cubic shape, which can not be easily reconciled with its hexameric structure. Phylogenetic analysis revealed that LSP-2 diverged from the LSP-1 like hexamerins after separation of the Diptera from other insect orders.

Redundant cis-acting elements control expression of the Drosophila affinidisjuncta Adh gene in the larval fat body

The alcohol dehydrogenase (Adh) gene in the Hawaiian species of fruit fly, Drosophila affinidisjuncta, like the Adh genes from all Drosophila species analyzed, is expressed at high levels in the larval fat body via a larval-specific promoter. To identify the cis-acting elements involved in this highly conserved aspect of Adh gene expression, deleted D. affinidisjuncta genes were introduced into D. melanogaster by somatic transformation. Unlike previously described methods, this transformation system allows analysis of Adh gene expression specifically in the larval fat body. The arrangement of sequences influencing expression of the proximal promoter of this gene in the larval fat body differs markedly from that described for the Adh gene from the distant relative, D. melanogaster. Multiple redundant elements dispersed 5' and 3' to the gene, only some of which map to regions carrying evolutionarily conserved sequences, affect expression in the fat body. D. affinidisjuncta employs a novel mode of Adh gene regulation in which the proximal promoter is influenced by sequences having roles in expression of the distal promoter. This gene is also unique in that far upstream sequences can compensate for loss of sequences within 200 bp of the proximal RNA start site. Furthermore, expression is influenced in an unusual, context-dependent manner by a naturally-occurring 3' duplication of the proximal promoter--a feature found only in Hawaiian species.

The sequence and expression pattern of the Calliphora erythrocephala yolk protein A and B genes

The yolk protein genes (yps) are expressed in a temporal, tissue- and sex-specific fashion in Drosophila melanogaster. Here we report the sequence of two related genes in Calliphora erythrocephala. The predicted Calliphora yolk protein (YP) sequences are well conserved, especially at the C-terminal end when compared to those of D. melanogaster and Ceratitis capitata. Database searches with the Calliphora yolk protein B (CeYPB) sequence identify the vertebrate lipase similarity reported for the YPs of Drosophila and Ceratitis. Moreover, sequences with identity to divalent ion-binding sites were observed, which colocalized with putative tyrosine sulfation sites. Calliphora oogenesis differs from Drosophila in that it is cyclic in response to a meat feed. The Calliphora yp genes are expressed in the follicle cells of the egg chamber during vitellogenesis, as shown by in situ hybridization, and the yp message levels correlate with YP synthesis. The synthesis of the yp transcripts in ovaries of Calliphora occurs in the same pattern as that for ovarian transcripts in Drosophila. In the carcass, yp transcript levels are correlated with the production of a batch of eggs.

Molecular phylogeny of Drosophila based on ribosomal RNA sequences

Nucleotide sequences of 72 species of Drosophilidae were determined for divergent D1 and D2 domains (representing 200 and 341 nucleotides respectively in D. melanogaster) of large ribosomal RNA, using the rRNA direct sequencing method. Molecular phylogenetic trees were reconstructed using both distance and parsimony methods and the robustness of the nodes was evaluated by the bootstrap procedure. The trees obtained by these methods revealed four main lineages or clades which do not correspond to the taxonomical hierarchy. In our results, the genus Chymomyza is associated with the subgenus Scaptodrosophila of the genus Drosophila and their cluster constitutes the most ancient clade. The two other clades are constituted of groups belonging to the subgenus Sophophora of the genus Drosophila: the so-called Neotropical clade including the willistoni and saltans groups and the obscura-melanogaster clade itself split into three lineages: (1) obscura group + ananassae subgroup, (2) montium subgroup, and (3) melanogaster + Oriental subgroups. The fourth clade, the Drosophila one, contains three lineages. D. polychaeta, D. iri, and D. fraburu are branched together and constitute the most ancient lineage; the second lineage includes the annulimana, bromeliae, dreyfusi, melanica, mesophragmatica, repleta, robusta, and virilis groups. The third lineage is composed of the immigrans and the cardini, funebris, guaramunu, guarani, histrio, pallidipennis, quinaria, and tripunctata groups. The genera Samoaia, Scaptomyza, and Zaprionus are branched within the Drosophila clade. Although these four clades appear regularly in almost all tree calculations, additional sequencing will be necessary to determine their precise relationships.

The sequence-immunology correlation revisited: data for cetacean myoglobins and mammalian lysozymes

Quantitative microcomplement fixation tests employing rabbit antisera were done to compare immunologically 13 cetacean myoglobins and 15 mammalian lysozymes c of known amino acid sequence. In both cases there was a strong correlation between immunological distance (y) and percent sequence difference (x), as had been found for several other globular proteins. For myoglobin the relationship could be described by y = 10.5x and for lysozyme by y = 8.5x. The coefficients in both of these equations are appreciably higher than the values of 5.1-6.9 reported for three other vertebrate globular proteins (bird lysozyme c, mammalian ribonuclease, and mammalian serum albumin), and they imply that rabbit antisera to mammalian myoglobins and lysozymes are more sensitive to evolutionary substitutions. A strong inverse correlation (r = -0.95) was found when the slope of the line relating y to x for these five data sets was plotted against the percent sequence difference between the rabbit's own protein and the proteins immunized with. Specifically, the cetacean myoglobins on average differ in amino acid sequence from rabbit myoglobin by less than 13% and exhibit the steepest slope (10.5), while bird lysozyme sequences differ by nearly 40% from rabbit lysozyme and exhibit the shallowest slope (5.1).

Adh and Adh-dup sequences of Drosophila lebanonensis and D. immigrans: interspecies comparisons

We have cloned and sequenced the Adh genomic region of Drosophila lebanonensis (subgenus Scaptodrosophila) and D. immigrans (subgenus Drosophila). This region, which contains Adh, encoding the alcohol dehydrogenase enzyme, and Adh-dup (duplicate of Adh), has been compared with the same fragment from D. subobscura (subgenus Sophophora). Even though the flanking regions and introns of both genes have been affected by high substitution rates, the consensus sequences have been clearly identified. Although the overall homology of the coding regions was 76-78% among the species compared, there were differences in the exon distribution of the nucleotide substitutions when Adh or Adh-dup were compared, thus showing that these two genes differ in their evolutionary pattern.

Homeotic gene Antennapedia mRNA contains 5'-noncoding sequences that confer translational initiation by internal ribosome binding

The Antennapedia (Antp) homeotic gene of Drosophila melanogaster has two promoters, P1 and P2. The resulting Antp mRNAs contain 1512-nucleotide (P1) and 1727-nucleotide (P2) 5'-noncoding regions, composed of exons A, B, D, and E (P1) or exons C, D, and E (P2), respectively. Multiple AUG codons are present in exons A, B, and C. We have found that 252-nucleotide exon D, common to mRNAs from both transcription units and devoid of AUG codons, can mediate initiation of translation by internal ribosome binding in cultured cells. Many mRNAs in Drosophila contain long 5'-noncoding regions with apparently unused AUG codons, suggesting that internal ribosome binding may be a common mechanism of translational initiation, and possibly its regulation, in Drosophila.

Complementary DNA-DNA hybridization in Drosophila

We have performed DNA-DNA hybridization experiments among several species of Drosophila using the evolutionarily conserved portion of the genome representing sequences coding for amino acids of proteins. This was done by using as tracer, radioactively labeled complementary DNA that was reverse transcribed from adult mRNA. We show that this procedure extends phylogenetically the distance over which the technique can be applied to fast-evolving groups such as Drosophila. The major phylogenetic conclusions are (1) the subgenus Sophophora is a monophyletic lineage; (2) within Sophophora the melanogaster subgroup is closer to the obscura group than either group is to the willistoni group; (3) the subgenus Drosophila is complex with most major lineages originating deep in the phylogeny; the subgenus may not be monophyletic; (4) as with most groups classically placed in Drosophila, the Hawaiian Drosophila originate early, supporting the notion that this lineage is older than the extant islands; and (5) the virilis/repleta lineage is monophyletic within Drosophila.

The phylogenetic relationships of flies in the family Drosophilidae deduced from mtDNA sequences

The phylogenetic relationships of several taxa from representative genera, subgenera, groups, and subgroups in the Drosophilidae were examined using sequences from a 905-bp mtDNA fragment. Conventional cloning and sequencing techniques were used to obtain nucleotide sequences. In addition, polymerase chain reaction primers were designed for the rapid amplification and sequencing of this region for the species examined in the Drosophilidae. Phylogenetic analysis was done by cladistic techniques. Because of the coding nature of the 905-bp mtDNA fragment, several separate analyses of these sequences were performed. The genera Scaptomyza and Hirtodrosophila occupy ancestral branching positions in the molecular phylogeny. The genera Chymomyza and Zaprionus have intermediate branching positions, while the subgenera Drosophila and Sophophora are in the most derived position in the molecular phylogeny. Within the subgenus Sophophora, there is little resolution using these sequences, while within the subgenus Drosophila, D. melanica, D. funebris, and D. pinicola form a clade in a derived part of the phylogeny, with D. robusta and D. immigrans branching in an intermediate position in the phylogeny. D. mercatorum, a member of the repleta species group, occupies an ancestral position in the molecular phylogeny.

Molecular evolution and phylogeny of the Drosophila virilis species group as inferred by two-dimensional electrophoresis

Systematic relationship among the 12 species of the Drosophila virilis species group, and Drosophila robusta, were investigated by the use of two-dimensional electrophoresis (2-DE). A total of 389 protein characters (about 200 loci) were scored and analyzed both phylogenetically and phenetically. The resulting phylogeny was found to be largely concordant with the current views of evolution among these species based on other independent morphological, chromosomal, electrophoretic, and immunological data sets, although some notable differences were observed. The 2-DE data also appeared to be useful for constructing a molecular clock to date the absolute times of divergence among the species. It appears from this analysis that the evolution of the major clades within the species group occurred about 20 million years ago. Previous suggestions that the rate of molecular evolution was different between the virilis and montana phylads was not confirmed. The technique of 2-DE seems to be an excellent tool for reconstructing phylogenies and should be particularly valuable for examining relatively closely related species.

ADH and phylogenetic relationships of Drosophila lebanonesis (Scaptodrosophila)

Increasing data on Drosophila alcohol dehydrogenase (ADH) sequences have made it possible to calculate the rate of amino acid replacement per year, which is 1.7 x 10(-9). This value makes this protein suitable for reconstructing phylogenetic relationships within the genus for those species for which no molecular data are available such as Scaptodrosophila. The amino acid sequence of Drosophila lebanonensis is compared to all of the already known Drosophila ADHs, stressing the unique characteristic features of this protein such as the conservation of an initiating methionine at the N-terminus, the unique replacement of a glycine by an alanine at a very conserved position in the NAD domain of all dehydrogenases, the lack of a slow-migrating peptide, and the total conservation of the maximally hydrophilic peptide. The functional significance of these features is discussed. Although the percent amino acid identity of the ADH molecule in Drosophila decreases as the number of sequences compared increases, the conservation of residue type in terms of size and hydrophobocity for the ADH molecule is shown to be very high throughout the genus Drosophila. The distance matrix and parsimony methods used to establish the phylogenetic relationships of D. lebanonensis show that the three subgenera, Scaptodrosophila, Drosophila, and Sophophora separated at approximately the same time.

Molecular phylogeny of the subgenus Sophophora of Drosophila derived from large subunit of ribosomal RNA sequences

RNA sequencing has been used to assess the relationships among species of the subgenus Sophophora of the genus Drosophila. Two divergent domains, D1 and D2, of the large ribosomal RNA (28S), totalling 550 nucleotides have been sequenced using the rRNA direct sequencing method. A tree has been reconstructed from the neighbor-joining algorithm and the confidence intervals were evaluated by the bootstrap procedure. Results have shown that the branching of the willistoni and saltans groups of the subgenus Sophophora is very ancient and probably predates that of the subgenus Drosophila. The other groups and subgroups of Sophophora are clustered in three main lineages: 1) the melanogaster and oriental subgroups; 2) the montium subgroup; 3) the ananassae subgroup of the melanogaster group clustered with the fima and obscura groups. Thus, in comparison with our results, several taxa of various ranks appear paraphyletic (the genus Drosophila, the subgenus Sophophora and the melanogaster group). Our biochemical phylogeny is only in partial agreement with the pattern of Throckmorton's radiations as well as with classical taxonomy, both based on morphological data.

Nucleotide sequence of the genomic region encoding alcohol dehydrogenase in Drosophila affinidisjuncta

The DNA sequence of a 3886-bp genomic region containing the alcohol dehydrogenase (Adh) gene from Drosophila affinidisjuncta, and the RNA sequences of the D. affinidisjuncta Adh transcripts, are presented. These data support the conclusion that two Adh promoters generate distinct, developmentally regulated Adh transcripts. Correlations between these sequences and the transcription map are discussed. Comparisons between these and equivalent data from D. melanogaster are also presented. We note the following observations: (1) Except at the extreme 3' end, the two genes are identically organized. (2) Drosophila Adh protein accumulates amino acid replacements at the rate of approximately 0.5 per million years. (3) Among the non-protein-coding DNA sequences, putative homologies occur in the two promoter regions.

Molecular evolution among some Drosophila species groups as indicated by two-dimensional electrophoresis

The evolutionary and phylogenetic relationships of seven Drosophila species groups (represented by D. melanogaster, D. mulleri, D. mercatorum, D. robusta, D. virilis, D. immigrans, D. funebris, and D. melanica) were investigated by the use of two-dimensional electrophoresis. The resulting phylogeny is congruent with the current views of evolution among these groups based on morphological characters and immunological distances. Previous studies indicated that the ability of one-dimensional electrophoresis to resolve relationships between distantly related taxa extended to about the Miocene [25 million years (Myr) ago], but the present study demonstrates that two-dimensional electrophoresis is a useful indicator of phylogeny even back to the Paleocene (65 Myr ago). In addition, two-dimensional electrophoresis is shown to be a useful technique for detecting slowly evolving structural proteins such as actins and tropomyosins.

Biochemical phylogeny of the eight species in the Drosophila melanogaster subgroup, including D. sechellia and D. orena

The phylogenetic relationships of the eight species of the Drosophila melanogaster subgroup are examined on the basis of genetic variation at 33 putative enzyme loci. Values of Nei's genetic distance (ds) range from 0·28 to 1·74. D. sechellia appears closer to D. simulans than to D. mauritiana, the two former being the most closely related. D. orena is quite distantly related to D. erecta (ds = 1). Genetic differentiation supports the existence of three main lineages within the melanogaster subgroup and the yakuba-teissieri pair appears to be closer to the melanogaster lineage than to the erecta-orena one. Inferences of the times of species divergence from allozyme data are made and their agreement to other estimates is discussed.

Tempo and mode of sequence evolution in mitochondrial DNA of Hawaiian Drosophila

Sequence comparisons were made for up to 667 bp of DNA cloned from 14 kinds of Hawaiian Drosophila and five other dipteran species. These sequences include parts of the genes for NADH dehydrogenase (subunits 1, 2, and 5) and rRNA (from the large ribosomal subunit). Because the times of divergence among these species are known approximately, the sequence comparisons give insight into the evolutionary dynamics of this molecule. Transitions account for nearly all of the differences between sequences that have diverged by less than 2%: for these sequences the mean rate of divergence appears to be about 2%/Myr. In comparisons involving greater divergence times and greater sequence divergence, relatively more of the sequence differences are due to transversions. Specifically, the fraction of these differences that are counted as transversions rises from an initial value of less than 0.1 to a plateau value of nearly 0.6. The time required to reach half of the plateau value, about 10 Myr, is similar to that for mammalian mtDNA. The mtDNAs of flies and mammals are also alike in the shape of the curve relating the percentage of positions at which there are differences in protein-coding regions to the time of divergence. For both groups of animals, the curve has a steep initial slope ascribable to fast accumulation of synonymous substitutions and a shallow final slope resulting from the slow accumulation of substitutions causing amino acid replacements. However, the percentage of all sites that can experience a high rate of substitution appears to be only about 8% for fly mtDNA compared to about 20% for mammalian mtDNA.(ABSTRACT TRUNCATED AT 250 WORDS)

DNA sequence comparison among closely related Drosophila species in the mulleri complex

DNA hybridization was used to establish DNA sequence relationships among seven Drosophila species. Single-copy DNA was isolated from four species within the Drosophila mulleri complex, D. mojavensis, D. arizonensis, D. ritae and D. starmeri. These single-copy DNAs were used as tracers to be hybridized with each other and one additional member of the mulleri complex, D. aldrichi, a member of a closely related complex, D. hydei, and a distantly related species, D. melanogaster. Two methods have been used to determine the relatedness between these species: (1) the extent of duplex formed as measured by binding to hydroxyapatite and (2) the thermal stability of the duplexed DNA. Moderately repetitive DNA was purified from these species and used similarly to determine the divergence of this family of sequences. The rate of nucleotide substitution was estimated to be 0.2 +/-, 0.1% base pair change per million years for both single-copy and middle-repetitive DNAs. The size of the D. arizonensis genome, a representative of the mulleri complex, was calculated to be 2.2 X 10(8) base pairs from its kinetic complexity similar to that of D. hydei. The relative amounts (18%) and average reiteration frequency (100 copies) of the middle-repetitive DNA are similar for all Drosophila species studied. Finally, the data are presented in a phylogenetic tree.

Mitochondrial DNA evolution in the melanogaster species subgroup of Drosophila

Detailed restriction maps (40 cleavage sites on average) of mitochondrial DNAs (mtDNAs) from the eight species of the melanogaster species subgroup of Drosophila were established. Comparison of the cleavage sites allowed us to build a phylogenetic tree based on the matrix of nucleotide distances and to select the most parsimonious network. The two methods led to similar results, which were compared with those in the literature obtained from nuclear characters. The three chromosomally homosequential species D. simulans, D. mauritiana, and D. sechellia are mitochondrially very related, but exhibit complex phylogenetic relationships. D. melanogaster is their closest relative, and the four species form a monophyletic group (the D. melanogaster complex), which is confirmed by the shared unusual length of their mt genomes (18-19 kb). The other four species of the subgroup (D. yakuba, D. teissieri, D. erecta, and D. orena) are characterized by a much shorter mt genome (16-16.5 kb). The monophyletic character of the D. yakuba complex, however, is questionable. Two species of this complex, D. yakuba and D. teissieri, are mitochondrially indistinguishable (at the level of our investigation) in spite of their noticeable allozymic and chromosomal divergence. Finally, mtDNA distances were compared with the nuclear-DNA distances thus far established. These sequences seem to evolve at rather similar rates, the mtDNA rate being barely double that of nuclear DNA.

Sarcophagid larval proteins: partial sequence homologies among three cuticle proteins and related structures of drosophilids

Three structural proteins from the larval cuticle of Sarcophaga bullata have been sequenced at the amino terminus for 30-40 residues. We observed a high degree of homology with related proteins of Drosophila melanogaster, based on the previous findings of M. Snyder, J. Hirsh, and N. Davidson [(1981) Cell 25:165-177]. S. bullata protein SC1 had 65% homology with Drosophila isolate CP1, and SC6 showed 49% homology with CPX and 54% with CP2a. The three sarcophagid polypeptides also resembled each other with respect to mapped products of tryptic cleavage. The sites of posttranslational arylation required for puparium formation, namely histidyl and lysyl residues, were asymmetrically distributed in the sarcophagid samples. In SC1 the bulk of the loci of putative crosslinks lay beyond the 43-residue fragment. In SC6 half the histidines fell within the first 25% of the primary chain.

Ancient origin for Hawaiian Drosophilinae inferred from protein comparisons

Immunological comparisons of a larval hemolymph protein enabled us to build a tree relating major groups of drosophiline flies in Hawaii to one another and to continental flies. The tree agrees in topology with that based on internal anatomy. Relative rate tests suggest that evolution of hemolymph proteins has been about as fast in Hawaii as on continents. Since the absolute rate of evolution of hemolymph proteins in continental flies is known, one can erect an approximate time scale for Hawaiian fly evolution. According to this scale, the Hawaiian fly fauna stems from a colonist that landed on the archipelago about 42 million years ago-i.e., before any of the present islands harboring drosophilines formed. This date fits with the geological history of the archipelago, which has witnessed the sequential rise and erosion of many islands during the past 70 million years. We discuss the bearing of the molecular time scale on views about rates of organismal evolution in the Hawaiian flies.

A model of immunological distances in systematics

While immunological distances among taxa have had wide use in systematics, there has been some doubt about their utility because of the observed non-metricity of such distance matrices. A model is presented here relating observed immunological distance to the actual number of antigenic site differences between taxa. This model accounts for the observed departures of these distances from the metric condition of reciprocity and triangle inequality. Based upon the model, two procedures are suggested for the transformation of immunological distances to metric distances appropriate for phylogenetic analysis. The model implies that the usual scaling adjustments applied to the immunological distance matrix are inappropriate; however, the same transformation applied instead to an initial similarity matrix will solve a scaling problem. Non-reciprocity of the distances is shown to remain a problem independent of this initial scaling problem. It is suggested that further transformation of these re-scaled distances may be obtained through an extension of the ADCLUS procedure developed in psychology. This approach suggests a general strategy for a transformation to metric distances, given a particular model of non-metricity for the data.

Molecular evolution in Drosophila and the higher Diptera II. A time scale for fly evolution

In this paper, we examine first the steadiness of the rate of evolutionary change in a larval hemolymph protein, LHP, in numerous Drosophila species. We estimated amino acid sequence divergence from immunological distances measured with the quantitative microcomplement fixation technique. Using tests not depending on knowledge of absolute times of divergence, we estimated the variance of the rate of evolutionary change to be at least 4 times as large as that for a process resembling radioactive decay. Thus, the rate of evolution of this protein is as uniform as that of vertebrate proteins. Our analysis indicates no acceleration of protein evolution in the lineages leading to Hawaiian drosophilines. Second, we give an explicit description of a procedure for calculating the absolute value of the mean rate of evolutionary change in this protein. This procedure is suggested for general use in calculating absolute rates of molecular evolution. The mean rate of evolution of LHP is about 1.2 immunological distance units per million years, which probably corresponds to a unit evolutionary period of 4 million years; LHP thus evolves at a rate comparable to that of mammalian hemoglobins. Finally, we utilize the calibrated rate of LHP evolution to derive a time scale of evolution in the Drosophilidae and higher Diptera.

Evolution of single-copy DNA and the ADH gene in seven drosophilids

Single-copy DNA was isolated from Drosophila melanogaster and hybridized with total genomic DNA of D. melanogaster, D. mauritiana, D. simulans, D. pseudoobscura, D. willistoni, D. hydei and D. virilis. The duplexes were thermally eluted from hydroxyapatite and the data used to assess the relatedness of each species to D. melanogaster. The general pattern of relatedness was similar to that predicted by morphological methods but with some notable exceptions. The rate of nucleotide substitution was estimated to be greater than 0.66% of bases per million years. An unexpected, rapidly evolving component of D. melanogaster single-copy DNA was identified. The relatedness of these species was also studied with respect to the gene coding for alcohol dehydrogenase (ADH). The ADH gene, previously cloned from D. melanogaster (Goldberg 1980), was hybridized with Southern blots of genomic digests of the seven species. The intensity and position of the hybridizing bands suggest the amount of divergence of the gene. Divergence was quantitated by reassociation of a fragment of the cloned ADH gene with total DNA of the seven drosophilids and thermal elution of the resultant duplexes from hydroxyapatite. The ADH gene was isolated from genomic clone libraries of D. melanogaster, D. simulans and D. mauritiana and further studied by comparison of position of restriction sites. Species relationships deduced from comparison of total single-copy DNA and the ADH gene were consistent, demonstrating that a single gene can reflect divergence of the entire genome.