Evolution of yellow gene regulation and pigmentation in Drosophila

BACKGROUND: Changes in developmental gene expression are central to phenotypic evolution, but the genetic mechanisms underlying these changes are not well understood. Interspecific differences in gene expression can arise from evolutionary changes in cis-regulatory DNA and/or in the expression of trans-acting regulatory proteins, but few case studies have distinguished between these mechanisms. Here, we compare the regulation of the yellow gene, which is required for melanization, among distantly related Drosophila species with different pigment patterns and determine the phenotypic effects of divergent Yellow expression. RESULTS: Yellow expression has diverged among D. melanogaster, D. subobscura, and D. virilis and, in all cases, correlates with the distribution of black melanin. Species-specific Yellow expression patterns were retained in D. melanogaster transformants carrying the D. subobscura and D. virilis yellow genes, indicating that sequence evolution within the yellow gene underlies the divergence of Yellow expression. Evolutionary changes in the activity of orthologous cis-regulatory elements are responsible for differences in abdominal Yellow expression; however, cis-regulatory element evolution is not the sole cause of divergent Yellow expression patterns. Transformation of the D. melanogaster yellow gene into D. virilis altered its expression pattern, indicating that trans-acting factors that regulate the D. melanogaster yellow gene have also diverged between these two species. Finally, we found that the phenotypic effects of evolutionary changes in Yellow expression depend on epistatic interactions with other genes. CONCLUSIONS: Evolutionary changes in Yellow expression correlate with divergent melanin patterns and are a result of evolution in both cis- and trans-regulation. These changes were likely necessary for the divergence of pigmentation, but evolutionary changes in other genes were also required.     

Positive and negative DNA elements of the Drosophila grimshawi s18 chorion gene assayed in Drosophila melanogaster.

Germ line transformation has been used to map the cis regulatory DNA elements responsible for the precise and evolutionarily stable developmental expression of the s18 chorion gene. Constructs containing chimeric combinations of Drosophila melanogaster and D. grimshawi DNA regions, as well as D. grimshawi sequences alone, can direct expression in the follicular epithelium, in an s18-specific temporal and spatial pattern. The results indicate that both positive and negative regulatory elements can function when transferred from D. grimshawi to D. melanogaster. The first ca. 100 bp of the 5'-flanking DNA region constitute a minimal, developmentally regulated promoter, expression of which is inhibited by the next 100-bp DNA segment and activated by positive elements located further upstream. Expression of the minimal promoter can also be enhanced by more distant chorion regulatory elements, provided the inhibitory DNA segment is absent.     

Divergence of the Yellow Gene between Drosophila Melanogaster and D. Subobscura: Recombination Rate, Codon Bias and Synonymous Substitutions

The yellow (y) gene maps near the telomere of the X chromosome in Drosophila melanogaster but not in D. subobscura. Thus the strong reduction in the recombination rate associated with telomeric regions is not expected in D. subobscura. To study the divergence of a gene whose recombination rate differs between two species, the y gene of D. subobscura was sequenced. Sequence comparison between D. melanogaster and D. subobscura revealed several elements conserved in noncoding regions that may correspond to putative cis-acting regulatory sequences. Divergence in the y gene coding region between D. subobscura and D. melanogaster was compared with that found in other genes sequenced in both species. Both, yellow and scute exhibit an unusually high number of synonymous substitutions per site (p(s)). Also for these genes, the extent of codon bias differs between both species, being much higher in D. subobscura than in D. melanogaster. This pattern of divergence is consistent with the hitchhiking and background selection models that predict an increase in the fixation rate of slightly deleterious mutations and a decrease in the rate of fixation of slightly advantageous mutations in regions with low recombination rates such as in the y-sc gene region of D. melanogaster.     

The chorion genes of the medfly, Ceratitis capitata, I: Structural and regulatory conservation of the s36 gene relative to two Drosophila species.

We have used low stringency screening with the Drosophila melanogaster s36 chorion gene to recover its homologue from genomic and cDNA libraries of the medfly, Ceratitis capitata. The same gene has also been recovered from a genomic library of D. virilis. The medfly s36 gene shows similar developmental specificity as in Drosophila (early choriogenesis). It is also specifically amplified in ovarian follicles; this is the first report of chorion gene amplification outside the genus Drosophila. Alignments of s36 sequences from three species show that, in addition to its regulatory conservation, the s36 gene is extensively conserved in sequence, in a region corresponding to a central protein domain, and in short regions of 5' flanking DNA that might correspond to cis-regulatory elements.     

Molecular structure of a gypsy element of Drosophila subobscura (gypsyDs) constituting a degenerate form of insect retroviruses.

We have determined the nucleotide sequence of a 7.5 kb full-size gypsy element from Drosophila subobscura strain H-271. Comparative analyses were carried out on the sequence and molecular structure of gypsy elements of D.subobscura (gypsyDs), D.melanogaster (gypsyDm) and D.virilis (gypsyDv). The three elements show a structure that maintains a common mechanism of expression. ORF1 and ORF2 show typical motifs of gag and pol genes respectively in the three gypsy elements and could encode functional proteins necessary for intracellular expansion. In the three ORF1 proteins an arginine-rich region was found which could constitute a RNA binding motif. The main differences among the gypsy elements are found in ORF3 (env-like gene); gypsyDm encodes functional env proteins, whereas gypsyDs and gypsyDv ORF3s lack some motifs essential for functionality of this protein. On the basis of these results, while gypsyDm is the first insect retrovirus described, gypsyDs and gypsyDv could constitute degenerate forms of these retroviruses. In this context, we have found some evidence that gypsyDm could have recently infected some D.subobscura strains. Comparative analyses of divergence and phylogenetic relationships of gypsy elements indicate that the gypsy elements belonging to species of different subgenera (gypsyDs and gypsyDv) are closer than gypsy elements of species belonging to the same subgenus (gypsyDs and gypsyDm). These data are congruent with horizontal transfer of gypsy elements among different Drosophila spp.     

A 9.6 kb intervening sequence in D. virilis rDNA, and sequence homology in rDNA interruptions of diverse species of Drosophila and other diptera.

A large proportion of the 28S ribosomal RNA genes in Drosophila virilis are interrupted by a DNA sequence 9.6 kilobase pairs long. As regards both its presence and its position in the 28S gene (about two thirds of the way in), the D. virilis rDNA intervening sequence is similar to that found in D. melanogaster rDNA, but lengths differ markedly between the two species. Degrees of nucleotide sequence homology have been detected bewteen rDNA interruptions of the two species. This homology extends to putative rDNA intervening sequences in diverse higher diptera (other Drosophila species, the house fly and the flesh fly), but hybridization of cloned D. melanogaster and D. virilis rDNA interruption segments to DNA of several lower diptera has been negative. As is the case with melanogaster rDNA interruptions, segments of the virilis rDNA intervening sequence hybridize with non-rDNA components of the virilis genome, and interspecific homology may involve these non-rDNA sequences as well as rDNA interruptions. There is, however, evidence from buoyant density fractionation of DNA that the distributions of interruption-related sequences are distinct in D. melanogaster and D. virilis genomes. Moreover, thermal denaturation studies have indicated differing extents of homology between hybridizable sequences in D. virilis DNA and different segments of the D. melanogaster rDNA intervening sequence. We infer from our studies that rDNA intervening sequences are prevalent among higher diptera; that in the course of the evolution of these organisms, elements of the intervening sequences have been moderately to highly conserved; and that this conservation extends in at least two distantly related species of Drosophila to similar sequences found elsewhere in the genomes.     

Evidence for redundancy but not trans factor-cis element coevolution in the regulation of Drosophila Yp genes.

In Drosophila melanogaster and the endemic Hawaiian species D. grimshawi three Yolk protein (Yp) genes are expressed in a similar sex- and tissue-specific pattern. In contrast, DNA sequence comparisons of promoter/enhancer regions show low levels of similarity. We tested the functional significance of these observations by transforming D. melanogaster with the genomic region that includes the divergently transcribed D. grimshawi DgYp1 and DgYp2 genes; we found that the introduced genes were expressed in female fat body and in ovaries but not in males. Moreover, we found D. grimshawi proteins in the hemolymph and accumulating in ovaries. Using reporter constructs we showed that the intergenic region from D. grimshawi was sufficient to drive accurate expression, but some low level of ectopic expression was seen in males. Transforming D. melanogaster with constructs bearing deletions within the D. grimshawi intergenic region revealed only subtle effects in the overall level of expression, suggesting a high level of redundancy. Testing mutants in the sex-specific regulator doublesex revealed that it is capable of repressing the DgYp genes in males. Together, these data show that D. melanogaster trans-acting factors can regulate the in vivo pattern of DgYp expression and support the notion of a redundant and complex system of cis-acting elements.     

Organizational analysis of elav gene and functional analysis of ELAV protein of Drosophila melanogaster and Drosophila virilis.

Drosophila virilis genomic DNA corresponding to the D. melanogaster embryonic lethal abnormal visual system (elav) locus was cloned. DNA sequence analysis of a 3.8-kb genomic piece allowed identification of (i) an open reading frame (ORF) with striking homology to the previously identified D. melanogaster ORF and (ii) conserved sequence elements of possible regulatory relevance within and flanking the second intron. Conceptual translation of the D. virilis ORF predicts a 519-amino-acid-long ribonucleoprotein consensus sequence-type protein. Similar to D. melanogaster ELAV protein, it contains three tandem RNA-binding domains and an alanine/glutamine-rich amino-terminal region. The sequence throughout the RNA-binding domains, comprising the carboxy-terminal 346 amino acids, shows an extraordinary 100% identity at the amino acid level, indicating a strong structural constraint for this functional domain. The amino-terminal region is 36 amino acids longer in D. virilis, and the conservation is 66%. In in vivo functional tests, the D. virilis ORF was indistinguishable from the D. melanogaster ORF. Furthermore, a D. melanogaster ORF encoding an ELAV protein with a 40-amino-acid deletion within the alanine/glutamine-rich region was also able to supply elav function in vivo. Thus, the divergence of the amino-terminal region of the ELAV protein reflects lowered functional constraint rather than species-specific functional specification.     

Nucleotide variation at the no-on-transient A gene in Drosophila littoralis

The no-on-transient A (nonA) gene encodes a putative RNA-binding protein, and mutations in this gene are known to affect vision, male courtship song and viability in Drosophila melanogaster. Here we have sequenced the coding region of the nonA gene of Drosophila littoralis and compared it with those of Drosophila virilis and D. melanogaster. All portions of nonA appeared to be conserved between D. littoralis and D. virilis, while the 5' region of the gene of these two species showed high divergence from that of a more distantly-related species, D. melanogaster. The same was true for the glycine repeat regions. No significant deviation from neutrality was observed in the analysis of intraspecific nucleotide variation in 5' or 3' region of the nonA gene in D. littoralis population. Also, comparison of D. littoralis sequences with homologous sequence of D. virilis suggests that the gene is evolving neutrally in D. virilis group. Divergence of the 5' regions between D. virilis group species and D. melanogaster could be a result of positive selection, but this finding is obscured by the long divergence time of the species groups.     

Differentiation of Muller''s Chromosomal Elements D and E in the Obscura Group of Drosophila

Twenty-two markers located on Muller's elements D or E have been mapped by in situ hybridization in six species of the obscura group of Drosophila and in D. melanogaster. The obscura species can be grouped into a Palearctic cluster (D. subobscura, D. madeirensis and D. guanche) and a Nearctic one (D. pseudoobscura, D. persimilis and D. miranda). Eleven of the probes contain known genes: E74, Acp70A, Est5, hsp28/23, hsp83, emc, hsp70, Xdh, Acph-1, Cec and rp49. The remaining probes are recombinant phages isolated from a D. subobscura genomic library. All these markers hybridize to the putative homologous chromosome or chromosomal arm of elements D and E. Thus, these elements have conserved their genic content during species divergence. Chromosomal homologies proposed previously for each element among the species of the same cluster have been compared with the present results. The distribution of markers within each element has changed considerably as inferred from pairwise comparisons of obscura species included in the two different clusters. Only chromosomal segments defined by closely linked markers have been conserved: one such segment has been detected in element D and three in element E between D. subobscura and D. pseudoobscura.     

Identification and Characterization of the Cecropin Antibacterial Protein Gene Locus in Drosophila virilis

Cecropin is a type of antibacterial peptide that is synthesized in response to infection and has been characterized in many insect species and one mammal. The Cecropin locus of Drosophila melanogaster also contains the gene Andropin, which has been identified only in this species and encodes a male-specific antibacterial peptide. As a first step in studying the molecular evolution of the cecropin and andropin genes among Drosophila species, we have isolated genomic clones that cover the Cecropin locus in Drosophila virilis. The cloned region totals approximately 25 kb, within which a 9-kb fragment contains four cecropin genes and one pseudogene. All four genes have a high level of sequence homology to D. melanogaster Cecropin, about 80% identity in the coding regions, and the intron positions are conserved. As in D. melanogaster and other insects, κB-related cis-regulatory elements are found upstream of these cecropin genes. An Andropin-related sequence was not identified in D. virilis; however, genome Southern hybridizations suggest that Andropin-related sequences are present in at least the melanogaster species subgroup. Analysis of 19 insect cecropin genes identifies a common ancestral Cecropin before the divergence of Diptera and Lepidoptera. In addition, D. melanogaster and D. virilis can be identified by monophyletic clades for Cecropin. In contrast, the Lepidopteran species show polyphyletic relationships for duplicated cecropin genes. Received: 12 August 1996 / Accepted: 18 October 1996     

DNA Sequence Variation at the Period Locus Reveals the History of Species and Speciation Events in the Drosophila Virilis Group

The virilis phylad of the Drosophila virilis group consists of five closely related taxa: D. virilis, D. lummei, D. novamexicana, D. americana americana and D. americana texana. DNA sequences from a 2.1-kb pair portion of the period locus were generated in four to eight individuals from each of the five taxa. We found evidence of recombination and high levels of variation within species. We found no evidence of recent natural selection. Surprisingly there was no evidence of divergence between D. a. americana and D. a. texana, and they collectively appear to have had a large historical effective population size. The ranges of these two taxa overlap in a large hybrid zone that has been delineated in the eastern U.S. on the basis of the geographic pattern of a chromosomal fusion. Also surprisingly, D. novamexicana appears to consist of two distinct groups each with low population size and no gene flow between them.     

The nonA gene in Drosophila conveys species-specific behavioral characteristics

The molecular basis of species-specific differences in courtship behavior, a critical factor in preserving species boundaries, is poorly understood. Genetic analysis of all but the most closely related species is usually impossible, given the inviability of hybrids. We have therefore applied interspecific transformation of a single candidate behavioral locus, no-on-transient A (nonA), between Drosophila virilis and D. melanogaster, to investigate whether nonA, like the period gene, might encode species-specific behavioral information. Mutations in nonA can disrupt both visual behavior and the courtship song in D. melanogaster. The lovesong of nonA(diss) mutant males superficially resembles that of D. virilis, a species that diverged from D. melanogaster 40-60 mya. Transformation of the cloned D. virilis nonA gene into D. melanogaster hosts carrying a synthetic deletion of the nonA locus restored normal visual function (the phenotype most sensitive to nonA mutation). However, the courtship song of transformant males showed several features characteristic of the corresponding D. virilis signal, indicating that nonA can act as a reservoir for species-specific information. This candidate gene approach, together with interspecific transformation, can therefore provide a direct avenue to explore potential speciation genes in genetically and molecularly tractable organisms such as Drosophila.     

GEM, a cluster of repetitive sequences in the Drosophila subobscura genome

GEM is a new family of repetitive sequences detected in the D. subobscura genome. Two of the four described GEM elements encompass a heterogeneous central module, with no detectable ORF, flanked by two long inverted repeats. These elements are composed of a set of repetitive modules, which are inverted repeat (IR), direct repeat (DR), palindromic sequence (PS), long sequence (LS) and short sequence (SS). These five modules can be found either clustered or dispersed as single modules in the D. subobscura genome, in euchromatic and heterochromatic regions. In addition to the 3' region of Adh retrosequences, single IR and LS blocks were found associated with the promoter region of different genes, in particular, LS-like blocks have also been found associated with functional genes in D. melanogaster and D. virilis. Conversely, the DR block is highly similar to satellite DNAs from some other species of the obscura group. In addition, GEM elements share some structural features with IS elements described in different Drosophila species. It is likely that both GEM and IS sequences would be vestiges of an ancestral transposable element.     

Excess Polymorphism at the Adh Locus in DROSOPHILA MELANOGASTER

The evolutionary history of a region of DNA encompassing the Adh locus is studied by comparing patterns of variation in Drosophila melanogaster and its sibling species, D. simulans. An unexpectedly high level of silent polymorphism in the Adh coding region relative to the 5' and 3' flanking regions in D. melanogaster is revealed by a populational survey of restriction polymorphism using a four-cutter filter hybridization technique as well as by direct sequence comparisons. In both of these studies, a region of the Adh gene encompassing the three coding exons exhibits a frequency of polymorphism equal to that of a 4-kb 5' flanking region. In contrast, an interspecific sequence comparison shows a two-fold higher level of divergence in the 5' flanking sequence compared to the structural locus. Analysis of the patterns of variation suggest an excess of polymorphism within the D. melanogaster Adh locus, rather than lack of polymorphism in the 5' flanking region. An approach is outlined for testing neutral theory predictions about patterns of variation within and between species. This approach indicates that the observed patterns of variation are incompatible with an infinite site neutral model.