Evolution of the Adh locus in the Drosophila willistoni group: the loss of an intron, and shift in codon usage

We report here the DNA sequence of the alcohol dehydrogenase gene (Adh) cloned from Drosophila willistoni. The three major findings are as follows: (1) Relative to all other Adh genes known from Drosophila, D. willistoni Adh has the last intron precisely deleted; PCR directly from total genomic DNA indicates that the deletion exists in all members of the willistoni group but not in any other group, including the closely related saltans group. Otherwise the structure and predicted protein are very similar to those of other species. (2) There is a significant shift in codon usage, especially compared with that in D. melanogaster Adh. The most striking shift is from C to U in the wobble position (both third and first position). Unlike the codon-usage-bias pattern typical of highly biased genes in D. melanogaster, including Adh, D. willistoni has nearly 50% G + C in the third position. (3) The phylogenetic information provided by this new sequence is in agreement with almost all other molecular and morphological data, in placing the obscura group closer to the melanogaster group, with the willistoni group farther distant but still clearly within the subgenus Sophophora.      

Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.      

Phylogeny of the subgenus sophophora (Diptera: drosophilidae) based on combined analysis of nuclear and mitochondrial sequences

Sequences from the nuclear (nu) alcohol dehydrogenase gene, the nu 28S ribosomal RNA locus, and the mitochondrial cytochrome oxidase II gene were used both individually and in combined analyses to infer the phylogeny of the subgenus Sophophora (Diptera: Drosophilidae). We used several optimality criteria, including maximum likelihood, maximum parsimony, and minimum evolution, to analyze these partitions to test the monophyly of the subgenus Sophophora and its four largest species groups, melanogaster, obscura, saltans, and willistoni. Our results suggest that the melanogaster and obscura species groups are each monophyletic and form a closely related clade. The Neotropical clade, containing the saltans and willistoni species groups, is also recovered, as previous studies have suggested. While the saltans species group is strongly supported as monophyletic, the results of several analyses indicate that the willistoni species group may be paraphyletic with respect to the saltans species group.     

Phylogeny of Drosophila and related genera: conflict between molecular and anatomical analyses

Drosophila species are extensively used in biological research; yet, important phylogenetic relationships within the genus and with related genera remain unresolved. The combined data for three genes (Adh, Sod, and Gpdh) statistically resolves outstanding issues. We define the genus Drosophila inclusively so as to include Scaptomyza and Zaprionus (considered distinct genera in the taxonomy of Wheeler, 1981) but excluding Scaptodrosophila. The genus Drosophila so defined is monophyletic. The subgenus Sophophora (including the melanogaster, obscura, and willistoni groups) is monophyletic and the sister clade to all other Drosophila subgenera. The Hawaiian Drosophila (including Scaptomyza) is a monophyletic group, but the subgenus Drosophila is not monophyletic, because the immigrans group is more closely related to the subgenus Hirtodrosophila than to other species of the subgenus Drosophila, such as the virilis and repleta groups.     

How Drosophila change their combs: the Hox gene Sex combs reduced and sex comb variation among Sophophora species

SUMMARY Identification of the events responsible for rapid morphological variation during evolution can help understand how developmental processes are changed by genetic modifications and thus produce diverse body features and shapes. Sex combs, a sexually dimorphic structure, show considerable variation in morphology and numbers among males from related species of Sophophora , a subgenus of Drosophila . To address which evolutionary changes in developmental processes underlie this diversity, we first analyzed the genetic network that controls morphogenesis of a single sex comb in the model D. melanogaster . We show that it depends on positive and negative regulatory inputs from proximo-distal identity specifying genes, including dachshund, bric à brac , and sex combs distal . All contribute to spatial regulation of the Hox gene Sex combs reduced (Scr ), which is crucial for comb formation. We next analyzed the expression of these genes in sexually dimorphic species with different comb numbers. Only Scr shows considerable expression plasticity, which is correlated with comb number variation in these species. We suggest that differences in comb numbers reflect changes of Scr expression in tarsus primordia, and discuss how initial comb formation could have occurred in an ancestral Sophophora fly following regulatory modifications of developmental programs both parallel to and downstream of Scr .     

A Functional and Structural Analysis of the Sex Combs Reduced Locus of Drosophila Melanogaster

We have undertaken a developmental genetic analysis of the homeotic gene Sex combs reduced (Scr) of Drosophila melanogaster by examining embryonic and adult phenotypes of mutations affecting Scr gene function. Molecular mapping of Scr breakpoint lesions has defined a segment of greater than 70 kb of DNA necessary for proper Scr gene function. This region is split by the fushi tarazu (ftz) gene, with lesions affecting embryonic Scr function molecularly mapping to the region proximal (5') to ftz and those exhibiting polyphasic semilethality predominantly mapping distal (3') to ftz. Gain-of-function mutations are associated with genomic rearrangements and map throughout the Scr locus. Our analysis has revealed that the Scr locus encompasses genetic elements that are responsible for functions in both the embryonic and larval to adult periods of development. From these studies, we conclude that Scr is a complex genetic locus with an extensive regulatory region that directs functions required for normal head and thoracic development in both the embryo and the adult and that the regulation of Scr during these two periods is distinct.     

A phylogeny of Drosophilidae using the Amyrel gene: questioning the Drosophila melanogaster species group boundaries

In this study, the phylogenetic relationships of 164 species of the family Drosophilidae are discussed, using the Amyrel gene, a member of the α -amylase multigene family. This study focuses on numerous species groups in the subgenera Sophophora and Drosophila of the genus Drosophila but also includes other closely related genera. Nucleotide data were analysed by several methods: maximum parsimony, neighbour joining, maximum likelihood and Bayesian inference. Heterogeneity of base composition (mainly low GC contents in the species groups willistoni and saltans ) has been addressed. In all analyses, the genus Drosophila appeared paraphyletic. The subgenus Sophophora clearly appeared to be a monophyletic group, showing well-resolved clades, with the Neotropical groups arising in a basal position. Here, it is proposed to raise the species subgroups ananassae and montium to the rank of species group, and to restrict the melanogaster species group to the melanogaster subgroup plus the 'Oriental' subgroups, among which the suzukii subgroup is polyphyletic. Some related genera such as Zaprionus , Liodrosophila , Scaptomyza and Hirtodrosophila are clustered with, or inside the subgenus Drosophila , which is therefore paraphyletic and should be reviewed.     

Cloning and transcriptional analysis of the segmentation gene fushi tarazu of Drosophila

In the course of studying the Antennapedia (Antp) locus, we found that one of the 3' Antp exons has weak cross-homology to another gene affecting segmentation, fushi tarazu (ftz; meaning "not enough segments"), which is 30 kb to the left of Antp. Homozygous ftz- embryos die before hatching and lack alternate body segments. The reduced number of segments results from the fusion of the anterior portion of one segment with the posterior portion of the next segment. The ftz gene encodes a single 1.9 kb poly(A)+ RNA expressed exclusively from the early blastoderm to gastrula stages of embryonic development. The structure of the ftz gene has been analyzed by S1 nuclease mapping and by restriction mapping of a cDNA clone. The ftz gene consists of two exons, and it is the 3' exon that cross-hybridizes with the 3' exon of Antp. The role of ftz in cell determination is discussed.     

Molecular evolution of Drosophila odorant receptor genes

A total of 752 odorant receptor (Or) genes, including pseudogenes, were identified in 11 Drosophila species and named after their orthologs in Drosophila melanogaster. The 813 Or genes, including 61 from D. melanogaster, were classified into 59 orthologous groups that are well supported by gene phylogeny. By reconciling with the gene family phylogeny, we estimated the number of gene duplication/loss events and intron gain/loss events in the species phylogeny. We found that these events are particularly frequent in Drosophila grimshawi, Drosophila willistoni, and obscura group. More than half of the duplicated genes stay as tandem arrays, whose size range from 2 to 8. These genes vary in sequence and some likely underwent positive selection, indicating that the gene duplication was important for flies to acquire new olfactory functions. We hypothesize that Or genes conferred the basic olfactory repertoire to ancestral flies before the speciation of the Drosophila and Sophophora subgenera about 40 Mya. This repertoire has been largely maintained in the current species, whereas lineage-specific gene duplication seems to have led to additional specialization in some species in response to specific ecological conditions.     

A comparative study of the short term cold resistance response in distantly related Drosophila species: the role of regucalcin and frost

The molecular basis of short term cold resistance (indexed as chill-coma recovery time) has been mostly addressed in D. melanogaster, where candidate genes (Dca (also known as smp-30) and Frost (Fst)) have been identified. Nevertheless, in Drosophila, the ability to tolerate short term exposure to low temperatures evolved several times independently. Therefore, it is unclear whether variation in the same candidate genes is also responsible for short term cold resistance in distantly related Drosophila species. It should be noted that Dca is a candidate gene for cold resistance in the Sophophora subgenus only, since there is no orthologous gene copy in the Drosophila subgenus. Here we show that, in D. americana (Drosophila subgenus), there is a north-south gradient for a variant at the 5' non-coding region of regucalcin (a Dca-like gene; in D. melanogaster the proteins encoded by the two genes share 71.9% amino acid identities) but in our D. americana F2 association experiment there is no association between this polymorphism and chill-coma recovery times. Moreover, we found no convincing evidence that this gene is up-regulated after cold shock in both D. americana and D. melanogaster. Size variation in the Fst PEST domain (putatively involved in rapid protein degradation) is observed when comparing distantly related Drosophila species, and is associated with short term cold resistance differences in D. americana. Nevertheless, this effect is likely through body size variation. Moreover, we show that, even at two hours after cold shock, when up-regulation of this gene is maximal in D. melanogaster (about 48 fold expression change), in D. americana this gene is only moderately up-regulated (about 3 fold expression change). Our work thus shows that there are important differences regarding the molecular basis of cold resistance in distantly related Drosophila species.     

Duplicative and conservative transpositions of larval serum protein 1 genes in the genus Drosophila

Interspecific comparative molecular analyses of transposed genes and their flanking regions can help to elucidate the time, direction, and mechanism of gene transposition. In the Drosophila melanogaster genome, three Larval serum protein 1 (Lsp1) genes (alpha, beta and gamma) are present and each of them is located on a different chromosome, suggesting multiple transposition events. We have characterized the molecular organization of Lsp1 genes in D. buzzatii, a species of the Drosophila subgenus and in D. pseudoobscura, a species of the Sophophora subgenus. Our results show that only two Lsp1 genes (beta and gamma) exist in these two species. The same chromosomal localization and genomic organization, different from that of D. melanogaster, is found in both species for the Lsp1beta and Lsp1gamma genes. Overall, at least two duplicative and two conservative transpositions are necessary to explain the present chromosomal distribution of Lsp1 genes in the three Drosophila species. Clear evidence for implication of snRNA genes in the transposition of Lsp1beta in Drosophila has been found. We suggest that an ectopic exchange between highly similar snRNA sequences was responsible for the transposition of this gene. We have also identified the putative cis-acting regulatory regions of these genes, which seemingly transposed along with the coding sequences.     

Convergent intron loss of MRP1 in Drosophila and mosquito species

Previous study revealed that the MRP1 gene ortholog DMRP1/CG6214 of Drosophila melanogaster contains 12 exons in the coding region. In the current study, the genes of DMRP1/CG6214 from D. melanogaster and Drosophila virilis were compared, and the result indicated that D. virilis had an extra intron located in exon 2, implying that intron loss or gain might have occurred at this locus. To track the evolution of the extra intron (Intron Z), orthologous nucleotide sequences of 37 arthropod species were cloned or annotated. Based on phylogenetic analysis, we found that Intron Z should present in the common ancestor of arthropod species, more than 420 Ma. In addition, we found that Sophophora subgenus species and mosquito (Culex pipiens) lost Intron Z independently, showing evolutionary convergence.     

Evidence for Horizontal Transmission of the P Transposable Element between Drosophila Species

Several studies have suggested that P elements have rapidly spread through natural populations of Drosophila melanogaster within the last four decades. This observation, together with the observation that P elements are absent in the other species of the melanogaster subgroup, has lead to the suggestion that P elements may have entered the D. melanogaster genome by horizontal transmission from some more distantly related species. In an effort to identify the potential donor in the horizontal transfer event, we have undertaken an extensive survey of the genus Drosophila using Southern blot analysis. The results showed that P-homologous sequences are essentially confined to the subgenus Sophophora. The strongest P hybridization occurs in species from the closely related willistoni group. A wild-derived strain of D. willistoni was subsequently selected for a more comprehensive molecular examination. As part of the analysis, a complete P element was cloned and sequenced from this line. Its nucleotide sequence was found to be identical to the D. melanogaster canonical P, with the exception of a single base substitution at position 32. When the cloned element was injected into D. melanogaster embryos, it was able to both promote transposition of a coinjected marked transposon and induce singed-weak mutability, thus demonstrating its ability to function as an autonomous element. The results of this study suggest that D. willistoni may have served as the donor species in the horizontal transfer of P elements to D. melanogaster.     

Distinct activation patterns of EGF receptor signaling in the homoplastic evolution of eggshell morphology in genus Drosophila

Homoplasy is a phenomenon in which organisms in different phylogenetic groups independently acquire similar traits. However, it is largely unknown how developmental mechanisms are altered to give rise to homoplasy. In the genus Drosophila, all species of the subgenus Sophophora, including Drosophila (D.) melanogaster, have eggshells with two dorsal appendages (DAs); most species in the subgenus Drosophila, including D. virilis, and in the subgenus Dorsilopha, have four-DAs. D. melanica belongs to the Drosophila subgenus, but has two-DAs, and phylogenetic analyses suggest that it acquired this characteristic independently. The patterning of the DAs is tightly regulated by epidermal growth factor receptor (EGFR) signaling in D. melanogaster. Previous studies suggested that a change in the EGFR signal activation pattern could have led to the divergence in DA number between D. melanogaster and D. virilis. Here, we compared the patterns of EGFR signal activation across the Drosophila subgenera by immunostaining for anti-activated MAP kinase (MAPK). Our analysis revealed distinct patterns of EGFR signal activation in each subgenus that was consistent with their phylogenetic relationship. In addition, the number of DAs always corresponded to the number of EGFR signaling activation domains in two, three, and four-DA species. Despite their common two-DA characteristic, the EGFR signaling activation pattern in D. melanica diverged significantly from that of species in the subgenus Sophophora. Our results suggest that acquisition of the homoplastic two-DA characteristic could be explained by modifications of the EGFR signaling system in the genus Drosophila that occurred independently and at least twice during evolution.