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.      

Tissue-Specific Expression Phenotypes of Hawaiian Drosophila Adh Genes in Drosophila Melanogaster Transformants

Interspecific differences in the tissue-specific patterns of expression displayed by the alcohol dehydrogenase (Adh) genes within the Hawaiian picture-winged Drosophila represent a rich source of evolutionary variation in gene regulation. Study of the cis-acting elements responsible for regulatory differences between Adh genes from various species is greatly facilitated by analyzing the behavior of the different Adh genes in a homogeneous background. Accordingly, the Adh gene from Drosophila grimshawi was introduced into the germ line of Drosophila melanogaster by means of P element-mediated transformation, and transformants carrying this gene were compared to transformants carrying the Adh genes from Drosophila affinidisjuncta and Drosophila hawaiiensis. The results indicate that the D. affinidisjuncta and D. grimshawi genes have relatively higher levels of expression and broader tissue distribution of expression than the D. hawaiiensis gene in larvae. All three genes are expressed at similar overall levels in adults, with differences in tissue distribution of enzyme activity corresponding to the pattern in the donor species. However, certain systematic differences between Adh gene expression in transformants and in the Hawaiian Drosophila are noted along with tissue-specific position effects in some cases. The implications of these findings for the understanding of evolved regulatory variation are discussed.     

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.     

Phylogeny of the Drosophila immigrans species group (Diptera: Drosophilidae) based on Adh and Gpdh sequences

The immigrans species group in the Drosophilinae is one of the representative species groups of Drosophila in East Asia. Although this group constitutes a significant part of the drosophilid fauna in the Old World, only a few species have been analyzed in previous molecular phylogenetic studies. To study the phylogeny of the immigrans group, we analyzed the nucleotide sequences of two nuclear genes, alcohol dehydrogenase (Adh) and glycerol-3-phosphate dehydrogenase (Gpdh), for 36 drosophilid species, including 12 species of the immigrans group. In the resultant phylogenetic trees, 10 species of the immigrans group (D. immigrans, D. formosana, D. ruberrima, D. albomicans, D. nasuta, D. neonasuta, D. pallidifrons, D. hypocausta, D. neohypocausta, D. siamana) consistently formed a clade (the immigrans group proper), although the phylogeny within this clade did not exactly correspond to the classification of species subgroups. However, D. annulipes and D. quadrilineata, both of which belong to the quadrilineata subgroup of the immigrans group, were not included in the immigrans group proper. Furthermore, we obtained the unexpected result that D. annulipes was included in a clade comprising Scaptomyza and Hawaiian Drosophila, together with D. maculinotata of the funebris group, although the phylogenetic relationships within this clade remain uncertain and need to be substantiated with further studies. Thus, according to the present study, the immigrans group is polyphyletic.     

DNA divergence in and around the alcohol dehydrogenase locus in five closely related species of Hawaiian Drosophila

The alcohol dehydrogenase (Adh) region from five planitibia subgroup species of Hawaiian picture-wing Drosophila has been cloned. A total of 15 kb of DNA in and around the Adh gene has been compared among the five species. Genetic distances were calculated to determine evolutionary relationships. These distances agree with previous distances determined by protein polymorphism and DNA hybridization techniques and can be interpreted in terms of specific island colonization and speciation (founder) events over the past 5 Myr. Examination of the restriction maps of the cloned Adh region from the five species shows many instances of small deletions, insertion of a transposable element in D. heteroneura, and the existence of a highly variable region on the 3' side of the Adh gene. Clustering relationships and rates of DNA change are calculated and compared with the relationship found for other species of Drosophila.      

The relationship between third-codon position nucleotide content, codon bias, mRNA secondary structure and gene expression in the drosophilid alcohol dehydrogenase genes Adh and Adhr

To gain insights into the relationship between codon bias, mRNA secondary structure, third-codon position nucleotide distribution, and gene expression, we predicted secondary structures in two related drosophilid genes, Adh and Adhr, which differ in degree of codon bias and level of gene expression. Individual structural elements (helices) were inferred using the comparative method. For each gene, four types of randomization simulations were performed to maintain/remove codon bias and/or to maintain or alter third-codon position nucleotide composition (N3). In the weakly expressed, weakly biased gene Adhr, the potential for secondary structure formation was found to be much stronger than in the highly expressed, highly biased gene Adh. This is consistent with the observation of approximately equal G and C percentages in Adhr ( approximately 31% across species), whereas in Adh the N3 distribution is shifted toward C (42% across species). Perturbing the N3 distribution to approximately equal amounts of A, G, C, and T increases the potential for secondary structure formation in Adh, but decreases it in Adhr. On the other hand, simulations that reduce codon bias without changing N3 content indicate that codon bias per se has only a weak effect on the formation of secondary structures. These results suggest that, for these two drosophilid genes, secondary structure is a relatively independent, negative regulator of gene expression. Whereas the degree of codon bias is positively correlated with level of gene expression, strong individual secondary structural elements may be selected for to retard mRNA translation and to decrease gene expression.     

Rates of DNA change and phylogeny from the DNA sequences of the alcohol dehydrogenase gene for five closely related species of Hawaiian Drosophila

The sequence of 1.6 kb of DNA surrounding the alcohol dehydrogenase (Adh) gene from five species of the Planitibia subgroup of the Hawaiian picture-winged Drosophila, with estimated divergence times of 0.4-5.1 Myr, has been determined. The gene trees which were found by using the sequence divergence from different regions of the sequences are generally in accord with the phylogeny proposed for these species when chromosomal inversions and island of origin are used. One of the species (D. picticornis) appears to be more distant from the other species in this group than they are from a member of the Grimshawi group (D. affinidisjuncta) which is chromosomally more distant. Two of the species (D. differens and D. plantibia) show heterogeneity in the nucleotide changes in the Adh coding region, heterogeneity which is interpreted to be due to a gene conversion or recombination after hybridization between the two species. The minimal rate of nucleotide substitution of synonymous nucleotides and of nontranscribed nucleotides downstream from the coding region is estimated as 1.5 x 10(-8) and 1.1 x 10(-8) substitutions/nucleotide/year, respectively. This rate is two to three times the maximal rate estimated for mammalian synonymous substitutions.     

Isolation and characterization of the genomic region from Drosophila kuntzei containing the Adh and Adhr genes

The nucleotide sequences of the Adh and Adhr genes of Drosophila kuntzei were derived from combined overlapping sequences of clones isolated from a genomic library and from cloned PCR and inverse-PCR fragments. Only a proximal promoter was detected upstream of the Adh gene, indicating that D. kuntzei Adh is regulated by a one-promoter system. Further upstream of the Adh structural gene, an adult enhancer region (AAE) was found that contains most of the regulatory sequences described for AAEs of other Drosophila species. Analysis of the ADH protein showed an amino acid change from valine to threonine in the active site at position 189 which is also found in D. funebris but is otherwise unique among DROSOPHILA: This difference alone may be responsible for the very low ADH activity found in this species and may cause a difference in substrate usage pattern. Codon bias in Adh and Adhr was comparable and found to be very low compared with other species. Phylogenetic analysis showed that D. kuntzei is closest related to D. funebris and D. immigrans. The time of divergence between D. kuntzei and D. funebris was estimated to be 14.2-20.2 Myr and that between D. kuntzei-D. funebris and D. immigrans to be 30.8-44.0 Myr. An analysis of the genetic variation in the Adh gene and upstream sequences of four European strains showed that this gene was highly variable. Overall nucleotide diversity (pi) was 0.0139, which is two times higher than that in D. melanogaster.     

Intergroup Phylogenies in Drosophila as Determined by Comparisons of Salivary Banding Patterns

A salivary gland chromosome phylogeny is presented which summarizes the evolutionary relationships of twenty-two species belonging to the sub-genus Drosophila, and members of the twelve species groups: D. melanica, D. repleta, D. carbonaria, D. polychaeta, D. annulimana, D. robusta, D. carsoni, D. virilis, D. funebris and the "picture-wing," D. mimica and D. crassifemur groups (of Hawaii).-Photographic salivary chromosome maps were prepared for all twenty-two species studied. While the chromosomes of different species belonging to the same group can usually be homologized almost completely, so that construction of intragroup phylogenies is easy, chromosomes of members of different groups are so modified structurally that in most cases only short sections can be fully homologized, and these in only one or two chromosome elements.-Broadly homologous chromosome elements were compared for three species at a time, and on the basis of overlapping homologous sections, or overlapping inversions included within homologous sections, the trio of chromosomes, and the species to which they belonged can often be arranged in a two-step phylogenetic series. Detection of many such ordered trios permits construction of a single phylogenetic scheme encompassing all species.-D. nigromelanica, of the D. melanica group is found to be chromosomally intermediate between the rest of its group and the species belonging to other groups, suggesting that it is the most nearly ancestral member of its group. When trios of species including D. nigromelanica and members of two other species groups are compared, it is found that in twelve of fourteen such comparisons the chromosomes of D. nigromelanica are structurally intermediate between those of the members of the other two species groups, indicating the central position of D. nigromelanica in the phylogeny as a whole.-Available cytological evidence indicates that among the nine continental groups studied, it is the D. robusta group which is chromosomally closest to the Hawaiian "picture-wing" groups. Among the members of the Hawaiian groups it is D. primaeva and D. attigua which are found to be closest to the continental species. This finding tends to confirm the earlier conclusion of Carson and Stalker, based on different evidence, that the above two species were in an ancestral position in the Hawaiian phylogeny.-The relationship of the D. robusta and D. melanica groups demonstrated in this paper, the phylogenies within each of these two groups earlier worked out by Narayanan and by Stalker, and the present geographical distributions of the species within them, require that at least three Asiatic-New World migrations must have occurred during the evolution of the two groups.     

Fungal diversity associated with hawaiian Drosophila host plants

Hawaiian Drosophila depend primarily, sometimes exclusively, on specific host plants for oviposition and larval development, and most specialize further on a particular decomposing part of that plant. Differences in fungal community between host plants and substrate types may establish the basis for host specificity in Hawaiian Drosophila. Fungi mediate decomposition, releasing plant micronutrients and volatiles that can indicate high quality substrates and serve as cues to stimulate oviposition. This study addresses major gaps in our knowledge by providing the first culture-free, DNA-based survey of fungal diversity associated with four ecologically important tree genera in the Hawaiian Islands. Three genera, Cheirodendron, Clermontia, and Pisonia, are important host plants for Drosophila. The fourth, Acacia, is not an important drosophilid host but is a dominant forest tree. We sampled fresh and rotting leaves from all four taxa, plus rotting stems from Clermontia and Pisonia. Based on sequences from the D1/D2 domain of the 26S rDNA gene, we identified by BLAST search representatives from 113 genera in 13 fungal classes. A total of 160 operational taxonomic units, defined on the basis of ≥97% genetic similarity, were identified in these samples, but sampling curves show this is an underestimate of the total fungal diversity present on these substrates. Shannon diversity indices ranged from 2.0 to 3.5 among the Hawaiian samples, a slight reduction compared to continental surveys. We detected very little sharing of fungal taxa among the substrates, and tests of community composition confirmed that the structure of the fungal community differed significantly among the substrates and host plants. Based on these results, we hypothesize that fungal community structure plays a central role in the establishment of host preference in the Hawaiian Drosophila radiation.     

Patterns of natural selection at the Alcohol dehydrogenase gene of Drosophila americana

Similar outcomes are often observed in species exposed to similar selective regimes, but it is unclear how often the same mechanism of adaptive evolution is followed. Here we present an analysis of selection affecting sequence variation in the Alcohol dehydrogenase (Adh) gene of Drosophila americana, a species endemic to a large climate range that has been colonized by D. melanogaster. Unlike D. melanogaster, there is no evidence of selection on allozymes of ADH across the sampled range. This indicates that if there has been a similar adaptive response to climate in D. americana, it is not within the coding region of Adh. Instead, analyses of a combined dataset containing 86 alleles of Adh reveal purifying selection on the Adh gene, especially within its intron sequences. Frequency spectra of derived unpreferred variants at synonymous sites indicate that these sites are affected by weak purifying selection, but the deviation from neutrality is less drastic than observed for derived variants in noncoding introns. This contrast further supports the notion that noncoding sites in Drosophila are often subject to stronger selection pressures than synonymous sites.     

The molecular evolution of the alcohol dehydrogenase and alcohol dehydrogenase-related genes in the Drosophila melanogaster species subgroup

The DNA sequences of the Adh genes of three members of the Drosophila melanogaster species subgroup have been determined. This completes the Adh sequences of the eight species of this subgroup. Two species, D. yakuba and D. teissieri, possess processed Adh pseudogenes. In all of the species of the subgroup, a gene of unknown function, Adhr, is located about 300 bp 3' to Adh. Although this gene is experiencing a higher rate of synonymous substitution than Adh, it is more constrained at the amino acid level. Phylogenetic relationships between all eight members of the melanogaster subgroup have been analyzed using a variety of methods. All analyses suggested that the D. yakuba and D. teissieri pseudogenes have a single common ancestor, rather than evolving independently in each species, and that D. melanogaster is the sister species to D. simulans, D. sechellia, and D. mauritiana. The evolutionary relationships of the latter three species remain equivocal.      

Evolution of alcohol dehydrogenase genes in peonies (Paeonia): phylogenetic relationships of putative nonhybrid species

Alcohol dehydrogenase genes were amplified by PCR, cloned, and sequenced from 11 putative nonhybrid species of the angiosperm genus Paeonia. Sequences of five exons and six intron regions of the Adh gene were used to reconstruct the phylogeny of these species. Two paralogous genes, Adh1A, and Adh2, were found; an additional gene, Adh1B, is also present in section Moutan. Phylogenetic analyses of exon sequences of the Adh genes of Paeonia and a variety of other angiosperms imply that duplication of Adh1 and Adh2 occurred prior to the divergence of Paeonia species and was followed by a duplication resulting in Adh1A and Adh1B. Concerted evolution appears to be absent between these paralogous loci. Phylogenetic analysis of only the Paeonia Adh exon sequences, positioning the root of the tree between the paralogous genes Adh1 and Adh2, suggests that the first evolutionary split within the genus occurred between the shrubby section Moutan and the other two herbaceous sections Oneapia and Paeonia. Restriction of Adh1B genes to section Moutan may have resulted from deletion of Adh1B from the common ancestor of sections Oneapia and Paeonia. A relative-rate test was designed to compare rates of molecular change among lineages based on the divergence of paralogous genes, and the results indicate a slower rate of evolution within the shrubby section Moutan than in section Oneapia. This may be responsible for the relatively long branch length of section Oneapia and the short branch length between section Moutan and the other two sections found on the Adh, ITS (nrDNA), and matK (cpDNA) phylogenies of the genus. Adh1 and Adh2 intron sequences cannot be aligned, and we therefore carried out separate analyses of Adh1A and Adh2 genes using exon and intron sequences together. The Templeton test suggested that there is not significant incongruence among Adh1A, ITS, and matK data sets, but that these three data sets conflict significantly with Adh2 sequence data. A combined analysis of Adh1A, ITS, and matK sequences produced a tree that is better resolved than that of any individual gene, and congruent with morphology and the results of artificial hybridization. It is therefore considered to be the current best estimate of the species phylogeny. Paraphyly of section Paeonia in the Adh2 gene tree may be caused by longer coalescence times and random sorting of ancestral alleles.      

The structure of the Adh locus of Drosophila mettleri: an intermediate in the evolution of the Adh locus in the repleta group of Drosophila.

Members of species of the mulleri and hydei subgroups of the repleta group of Drosophila have duplicate Adh genes. The Adh regions of D. mojavensis, D. mulleri, and D. hydei contain three genes--a pseudogene, Adh-2, and Adh-1--arranged 5' to 3'. To understand the evolution of the triplicate Adh structure, we have cloned and sequenced the Adh locus of D. mettleri. This region consists of a 5' pseudogene and a 3' functional Adh gene. On the basis of the structure and nucleotide sequence comparisons of Adh genes of D. mettleri and other species, we propose that an initial duplication of the ancestral Adh gene generated two Adh genes arranged in tandem. The more 5' Adh gene became a pseudogene, while the more 3' gene remained functional through all the developmental stages. A second duplication of this 3' gene resulted in Adh regions with three genes--a pseudogene, Adh-2, and Adh-1.     

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.     

Enzyme polymorphism and species discrimination in fruit flies of the genus Dacus (Tephritidae).

Sympatric populations of D. tryoni and D. neohumeralis are difficult to completely distinguish taxonomically. Using five pigmentation characters, each of some taxonomic value, a small proportion of individuals cannot be assigned to either species nor definitely classified as hybrids. To aid in species discrimination and hybrid identification gene frequencies in natural populations were estimated at three polymorphic protein loci, an alcohol dehydrogenase (Adh), an octanol dehydrogenase (Odh) and an esterase (E-2). Samples of flies were taken from four sites spread over 1200 miles along the Australian eastern coast. Within each species allelic frequencies at each locus were largely the same at all localities. Consistent differences in gene frequencies between species occurred at all three loci, strongly supporting the hypothesis of two distinct gene pools. The Adh locus best discriminated between species with a unique allele occurring in D. neohumeralis at a frequency of 0-85. None of the loci showed complete differentiation and hence it was not possible to find a quick and easy method to distinguish the species nor to detect field hybrids. Directional selection of laboratory populations for a change in callus colour (the best pigmentation character for separating the species) indicated that at the Adh and E-2 loci frequencies of major alleles were not genetically associated with major genes for callus colour. Thus genotype determination at these loci when considered together with pigmentation characters may be valuable taxonomically for further distinguishing between the species.     

Evolution of the Transposable Element Mariner in Drosophila Species

The distribution of the transposable element mariner was examined in the genus Drosophila. Among the eight species comprising the melanogaster species subgroup, the element is present in D. mauritiana, D. simulans, D. sechellia, D. yakuba and D. teissieri, but it is absent in D. melanogaster, D. erecta and D. orena. Multiple copies of mariner were sequenced from each species in which the element occurs. The inferred phylogeny of the elements and the pattern of divergence were examined in order to evaluate whether horizontal transfer among species or stochastic loss could better account for the discontinuous distribution of the element among the species. The data suggest that the element was present in the ancestral species before the melanogaster subgroup diverged and was lost in the lineage leading to D. melanogaster and the lineage leading to D. erecta and D. orena. This inference is consistent with the finding that mariner also occurs in members of several other species subgroups within the overall melanogaster species group. Within the melanogaster species subgroup, the average divergence of mariner copies between species was lower than the coding region of the alcohol dehydrogenase (Adh) gene. However, the divergence of mariner elements within species was as great as that observed for Adh. We conclude that the relative sequence homogeneity of mariner elements within species is more likely a result of rapid amplification of a few ancestral elements than of concerted evolution. The mariner element may also have had unequal mutation rates in different lineages.