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.     

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.     

Expression Pattern Diversity and Functional Conservation Between Retroposed PRAT Genes from Drosophila melanogaster and Drosophila virilis

Gene duplication by retrotransposition duplicates only the coding and untranslated regions of a gene and, thus, biases retroduplicated genes toward having different expression patterns from their parental genes. As such, genes duplicated by retrotransposition are more likely to develop novel expression domains. To explore this idea further, we used the Prat/Prat2 gene duplication in Drosophila as a case study to examine the aftermath of a retrotransposition event that resulted in both the parent and the child gene becoming essential for survival. We used the Gal4-UAS transgene system with EGFP as a reporter to determine the developmental expression patterns of Prat and Prat2 from D. melanogaster (DmPrat and DmPrat2) and Prat from D. virilis (DvPrat). We also tested the functional equivalence of the protein products of DmPrat and DmPrat2. We found that each of the proteins could rescue DmPrat mutations, showing that the requirement for both Prat and Prat2 in Drosophila is not simply due to differences in protein function. In contrast, we found that the DmPrat and DmPrat2 genes have developed nonoverlapping patterns of expression, which correlate with their respective loss-of-function phenotypes. We further found that DvPrat expression is similar to DmPrat during development but differs in adult gonads. Thus, the function of the Prat retrogene has not diverged in the D. melanogaster and D. virilis lineages, while some aspects of its expression pattern have evolved. Finally, we have identified promoter elements, conserved upstream of DmPrat and DvPrat, that this retrogene has acquired to drive its expression.     

Sequence and expression of the hunchback gene in Lucilia sericata: a comparison with other Dipterans

We have found that the hunchback (hb)gene from Lucilia sericata is conserved in its functional domains in comparison with related flies, although there is divergence in the protein outside these regions. The expression patterns of Lucilia hb in early embryos are broadly similar to other higher Dipterans. However, in the posterior region we report blastoderm and post-gastrulation expression patterns, which are diverged from Musca and Drosophila. These patterns are reminiscent of hb expression in more primitive insects and could be indicative of changes in the regulation of hb in Lucilia by the terminal system.     

Interspecific Comparison of the Transformer Gene of Drosophila Reveals an Unusually High Degree of Evolutionary Divergence

The transformer (tra) gene of Drosophila melanogaster occupies an intermediate position in the regulatory pathway controlling all aspects of somatic sexual differentiation. The female-specific expression of this gene's function is regulated by the Sex lethal (Sxl) gene, through a mechanism involving sex-specific alternative splicing of tra pre-mRNA. The tra gene encodes a protein that is thought to act in conjunction with the transformer-2 (tra-2) gene product to control the sex-specific processing of doublesex (dsx) pre-mRNA. The bifunctional dsx gene carries out opposite functions in the two sexes, repressing female differentiation in males and repressing male differentiation in females. Here we report the results from an evolutionary approach to investigate tra regulation and function, by isolating the tra-homologous genes from selected Drosophila species, and then using the interspecific DNA sequence comparisons to help identify regions of functional significance. The tra-homologous genes from two Sophophoran subgenus species, Drosophila simulans and Drosophila erecta, and two Drosophila subgenus species, Drosophila hydei and Drosophila virilis, were cloned, sequenced and compared to the D. melanogaster tra gene. This comparison reveals an unusually high degree of evolutionary divergence among the tra coding sequences. These studies also highlight a highly conserved sequence within intron one that probably defines a cis-acting regulator of the sex-specific alternative splicing event.     

Conservation of alternative splicing and genomic organization of the myosin alkali light-chain (Mlc1) gene among Drosophila species

The Mlc1 gene of Drosophila melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. In larval muscle and tubular and abdominal muscles of adults, all of the six exons are included in the spliced mRNA, whereas, in the fibrillar indirect flight muscle of adult, exon 5 is excluded from the mRNA. We show that this tissue-specific pattern of alternative splicing of the Mlc1 pre-mRNA is conserved in D. simulans, D. pseudoobscura, and D. virilis. Isolation and sequencing of the Mlc1 genes from these three other Drosophila species have revealed that the overall organization of the genes is identical and that the genes have maintained a very high level of sequence identity within the coding region. Pairwise amino acid identities are 94%-99%, and there are no charge changes among the proteins. Total nucleotide divergence within the coding region of the four genes supports the accepted genealogy of these species, but the data indicate a significantly higher rate of amino acid replacement in the branch leading to D. pseudoobscura. A comparison of nucleotide substitutions in the coding portions of exon 5 and exon 6, which encode the alternative carboxyl termini of the two MLC1 isoforms, suggests that exon 5 is subject to greater evolutionary constraints than is exon 6. In addition to the coding sequences, there is significant sequence conservation within the 5' and 3' noncoding DNA and two of the introns, including one that flanks exon 5. These regions are candidates for cis- regulatory elements. Our results suggest that evolutionary constraints are acting on both the coding and noncoding sequences of the Mlc1 gene to maintain proper expression and function of the two MLC1 polypeptides.      

Evolutionary conservation of homeodomain-binding sites and other sequences upstream and within the major transcription unit of the Drosophila segmentation gene engrailed.

The engrailed (en) gene functions throughout Drosophila development and is expressed in a succession of intricate spatial patterns as development proceeds. Normal en function relies on an extremely large cis-acting regulatory region (70 kilobases). We are using evolutionary conservation to help identify en sequences important in regulating patterned expression. Sequence comparison of 2.6 kilobases upstream of the en coding region of D. melanogaster and D. virilis (estimated divergence time, 60 million years) showed that 30% of this DNA occurs in islands of near perfect sequence conservation. One of these conserved islands contains binding sites for homeodomain-containing proteins. It has been shown genetically that homeodomain-containing proteins regulate en expression. Our data suggested that this regulation may be direct. The remaining conserved islands may contain binding sites for other regulatory proteins.     

Conservation of gene order, structure and sequence between three closely linked genes in Drosophila melanogaster and Drosophila virilis

In Drosophila melanogaster, the apparently unrelated genes anon-66Da, RpL14, and anon-66Db (from telomere to centromere) are located on a 5547 bp genomic fragment on chromosome arm 3L at cytological position 66D8. The three genes are tightly linked, and flanked by two relatively large genes with unknown function. We have taken a comparative genomic approach to investigate the evolutionary history of the three genes. To this end we isolated a Drosophila virilis 7.3 kb genomic fragment which is homologous to a 5.5 kb genomic region of D. melanogaster. Both fragments map to Muller's element D, namely to section 66D in D. melanogaster and to section 32E in D. virilis, and harbor the genes anon-66Da, RpL14, and anon-66Db. We demonstrate that the three genes exhibit a high conservation of gene topography in general and in detail. While most introns and intergenic regions reveal sequence divergences, there are, however, a number of interspersed conserved sequence motifs. In particular, two introns of the RpL14 gene contain a short, highly conserved 60 nt long sequence located at corresponding positions. This sequence represents a novel Drosophila small nucleolar RNA, which is homologous to human U49. Whereas DNA flanking the three genes shows no significant interspecies homologies, the 3'-flanking region in D. virilis contains sequences from the transposable element Penelope. The Penelope family of transposable elements has been shown to promote chromosomal rearrangements in the D. virilis species group. The presence of Penelope sequences in the D. virilis 7.3 kb genomic fragment may be indicative for a transposon-induced event of transposition which did not yet scramble the order of the three genes but led to the breakdown of sequence identity of the flanking DNA.     

Evolution of the glucose dehydrogenase gene in Drosophila

The glucose dehydrogenase genes (Gld) of Drosophila melanogaster, of D. pseudoobscura, and of D. virilis have been isolated and compared with each other in order to identify conserved and divergent aspects of their structure and expression. The exon/intron structure of Gld is conserved. The Gld mRNAs are similar, with a range of 2.6-2.8 kb among the three species. All three species exhibit peaks of Gld expression during every major developmental stage, although considerable variation in the precise timing of these peaks exists between species. Interspecific gene transfer experiments demonstrate that the regulation and function of the D. pseudoobscura Gld is similar enough to the homologous gene in D. melanogaster to substitute for its essential role in the eclosion process. Comparison of the putative promoter sequences has identified both shared and divergent sequence elements which are likely responsible, respectively, for the conserved and divergent patterns of expression observed. The entire coding sequences of the pseudoobscura and melanogaster Gld genes are presented and shown to encode a 612-amino-acid pre-protein. The inferred amino acid sequences are 92% conserved between the two species. In general the intronic regions of Gld are unusually well conserved.     

Genomic Organization and Evolution of Alternative Exons in a Drosophila Calcium Channel Gene

The genomic organization of a gene coding for an α1 subunit of a voltage-gated calcium channel of Drosophila melanogaster (Dmca1A) was determined. Thirty-four exons, distributed over 45 kb of genomic sequence, have been identified and mapped, including exons in three regions involved in alternative splicing and new sites potentially involved in RNA editing. The comparison of the intron/exon boundaries of this channel with a mammalian counterpart shows that the genomic structure of these two genes has remained fairly similar during evolution, with more than half of the Drosophila intron positions being perfectly conserved compared to the human channel. Phylogenetic analysis of the mutually exclusive alternative exons revealed that they have diverged considerably. It is suggested that this divergence, rather than reflecting evolutionary age, is the likely result of accelerated rates of evolution following duplication.     

Microevolutionary divergence pattern of the segmentation gene hunchback in Drosophila

To study the microevolutionary processes shaping the evolution of the segmentation gene hunchback (hb) from Drosophila melanogaster, we cloned and sequenced the gene from 12 isofemale lines representing wild-type populations of D. melanogaster, as well as from the closely related species Drosophila sechellia, Drosophila orena, and Drosophila yakuba. We find a relatively low degree of sequence variation in D. melanogaster (theta = 0.0017), which is, however, consistent with its chromosomal location in a region of low recombination. Tests of neutrality do not reject a neutral-evolution model for the whole region. However, pairwise tests with different subregions indicate that there is a relative excess of polymorphic sites in the leader and the intron. Codon usage pattern analysis shows a particularly biased codon usage in the highly conserved regions, which is in line with the hypothesis that selection on translational accuracy is the driving force behind such a bias. A comparison of the expression pattern of hb in different sibling species of D. melanogaster reveals some regulatory changes in D. yakuba, which could be interpreted as changes in the timing of secondary expression domains.     

Promoter analysis of the Drosophila melanogaster gene encoding the pterin 4alpha-carbinolamine dehydratase.

Pterin-4alpha-carbinolamine dehydratase (PCD) is a key enzyme in the regeneration pathway of tetrahydrobiopterin. Previously, we isolated and reported the Drosophila melanogaster gene encoding PCD. In the present study, we isolated and characterized the Drosophila virilis gene encoding PCD. The Drosophila virilis PCD gene has two introns and an open reading frame to encode a protein of 101 amino acids. The amino acid sequence of Drosophila virilis PCD shows a 83% homology to that of the Drosophila melanogaster PCD protein. From the alignment of the nucleotide sequence in the 5'-flanking region of the Drosophila melanogaster and Drosophila virilis PCD genes, we found four conserved sequences. Using a transient transfection assay, we showed that one of the conserved sequences (-127 to approximately -115) is critical for expression, also the minimal promoter region between -127 and +51 is necessary for the efficient expression of Drosophila melanogaster PCD.     

Spatial regulation of the gap gene giant during Drosophila development

We describe the regulated expression of the segmentation gene giant (gt) during early embryogenesis. The gt protein is expressed in two broad gradients in precellular embryos, one in anterior regions and the other in posterior regions. Double immunolocalization studies show that the gt patterns overlap with protein gradients specified by the gap genes hunchback (hb) and knirps (kni). Analysis of all known gap mutants, as well as mutations that disrupt each of the maternal organizing centers, indicate that maternal factors are responsible for initiating gt expression, while gap genes participate in the subsequent refinement of the pattern. The maternal morphogen bicoid (bcd) initiates the anterior gt pattern, while nanos (nos) plays a role in the posterior pattern. Gene dosage studies indicate that different thresholds of the bcd gradient might trigger hb and gt expression, resulting in overlapping but noncoincident patterns of expression. We also present evidence that different concentrations of hb protein are instructive in defining the limits of kni and gt expression within the presumptive abdomen. These results suggest that gt is a bona fide gap gene, which acts with hb, Krüppel and kni to initiate striped patterns of gene expression in the early embryo.     

Codon Usage Bias and Base Composition of Nuclear Genes in Drosophila

The nuclear genes of Drosophila evolve at various rates. This variation seems to correlate with codon-usage bias. In order to elucidate the determining factors of the various evolutionary rates and codon-usage bias in the Drosophila nuclear genome, we compared patterns of codon-usage bias with base compositions of exons and introns. Our results clearly show the existence of selective constraints at the translational level for synonymous (silent) sites and, on the other hand, the neutrality or near neutrality of long stretches of nucleotide sequence within noncoding regions. These features were found for comparisons among nuclear genes in a particular species (Drosophila melanogaster, Drosophila pseudoobscura and Drosophila virilis) as well as in a particular gene (alcohol dehydrogenase) among different species in the genus Drosophila. The patterns of evolution of synonymous sites in Drosophila are more similar to those in the prokaryotes than they are to those in mammals. If a difference in the level of expression of each gene is a main reason for the difference in the degree of selective constraint, the evolution of synonymous sites of Drosophila genes would be sensitive to the level of expression among genes and would change as the level of expression becomes altered in different species. Our analysis verifies these predictions and also identifies additional selective constraints at the translational level in Drosophila.