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.     

Two independent duplications forming the Cyp307a genes in Drosophila

The conserved relationship between orthologs of many cytochrome P450 genes involved in ecdysone synthesis is not reflected in the evolution of the Drosophila Cyp307a genes. In Drosophila melanogaster Cyp307a1 (spook) and Cyp307a2 (spookier) both play essential roles in ecdysone synthesis and may possess biochemically redundant functions. Using phylogenetic analyses we show that the Drosophila Cyp307a genes were formed from two independent duplication events depicting a complicated evolutionary scenario. An initial duplication, from a Cyp307a2 ancestral gene produced the Cyp307a1 gene that has been maintained only in the Sophophoran subgenus. A second duplication in the Drosophila subgenus formed an additional paralog, Cyp307a3. Microsynteny is conserved for Cyp307a2 throughout the Drosophila species, but is not conserved between Cyp307a1 and Cyp307a3. These are located in different genomic positions in the Sophophora and Drosophila subgenera, respectively. Cyp307a3 appears to encode a functional gene product and is expressed in a different spatial and temporal manner to Cyp307a1. This suggests some level of functional divergence between the Cyp307a paralogs in different Drosophila species.     

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.     

Photoperiodic regulation of cold tolerance and expression levels of regucalcin gene in Drosophila montana

Temperature-induced plasticity of cold tolerance has been reported in many insect species, but cold tolerance can also be affected by changes in day (or night) length. In the present study we elucidate the direct and indirect effects of photoperiod on the cold tolerance of females of two Drosophila montana strains--one which possesses a robust photoperiodic diapause and another which does not. In the diapause-strain the time needed for recovery from chill coma showed a positive correlation with day length, but diapause itself played only a minor role in photoperiodic acclimation. The strain that was not able to enter to diapause as a response to day length also lacked photoperiodic cold acclimation ability indicating that this strain has deficiencies in its photoperiodic time measurement system. In the diapause-strain, the expression level of regucalcin gene was more than two times higher in diapausing than in non-diapausing females maintained in a single photoperiod, but day length per se did not cause significant changes in expression levels of this gene in either of the strains. In the non-diapausing strain this gene showed no expression changes in any comparison. Overall, the study shows that a decrease in day length can induce cold acclimation in D. montana, while changes in regucalcin expression are linked with photoperiodic diapause.     

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.      

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.     

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.     

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 .     

Chromosomal evolution of elements B and C in the Sophophora subgenus of Drosophila: evolutionary rate and polymorphism

The locations of 77 markers along the chromosomal elements B (41 markers) and C (36 markers) of Drosophila subobscura, D. pseudoobscura, and D. melanogaster were obtained by in situ hybridization on polytene chromosomes. In comparisons between D. subobscura and D. pseudoobscura, 10 conserved segments (accounting for 32% of the chromosomal length) were detected on element B and eight (17% of the chromosomal length) on element C. The fixation rate of paracentric inversions inferred by a maximum likelihood approach differs significantly between elements. Muller's element C (0.17 breakpoints/Mb/million years) is evolving two times faster than element B (0.08 breakpoints/Mb/million years). This difference in the evolutionary rate is paralleled by differences in the extent of chromosomal polymorphism in the corresponding lineages. Element C is highly polymorphic in D. subobscura, D. pseudoobscura, and in other obscura group species such as D. obscura and D. athabasca. In contrast, the level of polymorphism in element B is much lower in these species. The fixation rates of paracentric inversions estimated in the present study between species of the Sophophora subgenus are the highest estimates so far reported in the genus for the autosomes. At the subgenus level, there is also a parallelism between the high fixation rate and the classical observation that the species of the Sophophora subgenus tend to be more polymorphic than the species of the Drosophila subgenus. Therefore, the detected relationship between level of polymorphism and evolutionary rate might be a general characteristic of chromosomal evolution in the genus Drosophila.     

P elements in the saltans group of Drosophila: a new evaluation of their distribution and number of genomic insertion sites

Few are studies on P elements that have addressed the saltans group. These studies had shown that species from the cordata and elliptica subgroups were devoid of any discernible P homologous sequences, while species from the parasaltans, sturtevanti, and saltans subgroups all contain P element sequences. Our analyses showed the presence of one to 15 P element insertion sites in species of the saltans group, including Drosophila neocordata and Drosophila emarginata (cordata and elliptica subgroups, respectively). From these species, only those from the parasaltans, sturtevanti, and saltans subgroups harbor canonical P elements and, only those of the last two subgroups seem to harbor putative full-sized elements. Due to the low similarity of the sequences found in D. neocordata and D. emarginata to those earlier described, we suggest that these sequences might be rudimental P element derivatives that were present in the ancestral of the subgenus Sophophora.     

Molecular evolution of a duplication: the sex-peptide (Acp70A) gene region of Drosophila subobscura and Drosophila madeirensis

In Drosophila melanogaster, the Acp70A gene, which is involved in the postmating reactions of the female, is a single-copy gene. However, in Drosophila subobscura, the gene is duplicated and both copies are transcribed. To study the molecular evolution of the duplication, a 2.1- kb fragment encompassing both copies of the duplication was sequenced for 10 lines of D. subobscura and one line of Drosophila madeirensis. Estimates of the divergence between the two copies of the duplicated region and between the two species studied, D. subobscura and D. madeirensis, revealed that both copies of the Acp70a gene had evolved independently since their duplication. The ratio of nonsynonymous to silent divergence between copies was generally higher than one. The McDonald and Kreitman test revealed an excess of nonsynonymous changes fixed since the duplication and before the split of the D. subobscura and D. madeirensis lineages. These results point to natural selection driving protein evolution after the duplication. Specifically, adaptive evolution appears to have caused the initial differentiation between copies of the N-terminal parts of the proteins, while purifying selection could be responsible for the high conservation of the C- terminal parts.      

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.     

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.     

Molecular evolution of the ligands of the insulin-signaling pathway: dilp genes in the genus Drosophila

Drosophila melanogaster, unlike mammals, has seven insulin-like peptides (DILPS). In Drosophila, all seven genes (dilp1-7) are single copy in the 12 species studied, except for D. grimshawi with two tandem copies of dilp2. Our comparative analysis revealed that genes dilp1-dilp7 exhibit differential functional constraint, which is indicative of some functional divergence. Species of the subgenera Sophophora and Drosophila differ in some traits likely affected by the insulin-signaling pathway, such as adult body size. It is in the branch connecting the two subgenera that we found the footprint left by positive selection driving nonsynonymous changes at some dilp1 codons to fixation. Finally, the similar rate at which the two dilp2 copies of D. grimshawi have evolved since their duplication and the presence of a putative regulatory region highly conserved between the two paralogs would suggest that both copies were preserved either because of subfunctionalization or dose dependency rather than by the neofunctionalization of one of the two copies.     

A statistical analysis of nucleotide substitutions in the Drosophila Adh region reflects irregularities in molecular clocks

Substitutions rates are expected to be rather constant when a gene is compared between species. To analyze this feature, Ka/Ks ratios have been studied for Alcohol dehydrogenase (Adh) and Alcohol dehydrogenase duplication (Adh-dup) genes in Drosophila species. Adh Ka/Ks values are lower in intrasubgenus comparisons involving species of the Sophophora group than when these species are compared to the D. immigrans and S. lebanonensis, and this difference does not occur in the Adh-dup comparisons.