Molecular evolution in Drosophila and the higher Diptera II. A time scale for fly evolution

Independent origins of new sex-linked chromosomes in the melanica and robusta species groups of Drosophila

Background

Recent translocations of autosomal regions to the sex chromosomes represent important systems for identifying the evolutionary forces affecting convergent patterns of sex-chromosome heteromorphism. Additions to the sex chromosomes have been reported in the melanica and robusta species groups, two sister clades of Drosophila. The close relationship between these two species groups and the similarity of their rearranged karyotypes motivates this test of alternative hypotheses; the rearranged sex chromosomes in both groups are derived through a common origin, or the rearrangements are derived through at least two independent origins. Here we examine chromosomal arrangement in representatives of the melanica and the robusta species groups and test these alternative hypotheses using a phylogenetic approach.

Results

Two mitochondrial and two nuclear gene sequences were used to reconstruct phylogenetic relationships of a set of nine ingroup species having fused and unfused sex chromosomes and representing a broad sample of both species groups. Different methods of phylogenetic inference, coupled with concurrent cytogenetic analysis, indicate that the hypothesis of independent origins of rearranged sex chromosomes within each species group is significantly more likely than the alternative hypothesis of a single common origin. An estimate tightly constrained around 8 My was obtained for the age of the rearranged sex chromosomes in the melanica group; however, a more loosely constrained estimate of 10–15 My was obtained for the age of the rearrangement in the robusta group.

Conclusion

Independent acquisition of new chromosomal arms by the sex chromosomes in the melanica and robusta species groups represents a case of striking convergence at the karyotypic level. Our findings indicate that the parallel divergence experienced by newly sex-linked genomic regions in these groups represents an excellent system for studying the tempo of sex chromosome evolution.

Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs

MicroRNA (miRNA) genes give rise to small regulatory RNAs in a wide variety of organisms. We used computational methods to predict miRNAs conserved among Drosophila species and large-scale sequencing of small RNAs from Drosophila melanogaster to experimentally confirm and complement these predictions. In addition to validating 20 of our top 45 predictions for novel miRNA loci, the large-scale sequencing identified many miRNAs that had not been predicted. In total, 59 novel genes were identified, increasing our tally of confirmed fly miRNAs to 148. The large-scale sequencing also refined the identities of previously known miRNAs and provided insights into their biogenesis and expression. Many miRNAs were expressed in particular developmental contexts, with a large cohort of miRNAs expressed primarily in imaginal discs. Conserved miRNAs typically were expressed more broadly and robustly than were nonconserved miRNAs, and those conserved miRNAs with more restricted expression tended to have fewer predicted targets than those expressed more broadly. Predicted targets for the expanded set of microRNAs substantially increased and revised the miRNA-target relationships that appear conserved among the fly species. Insights were also provided into miRNA gene evolution, including evidence for emergent regulatory function deriving from the opposite arm of the miRNA hairpin, exemplified by mir-10, and even the opposite strand of the DNA, exemplified by mir-iab-4.

Following the Chromosome Path to the Garden of the Genome

I have been fascinated by chromosomes for longer than I care to mention; their beautiful structure, cell-type-specific changes in morphology, and elegant movements delight me. Shortly before I began graduate study, the development of nucleic acid hybridization made it possible to compare two nucleic acids whether or not their sequences were known. From this stemmed a progression of development in tools and techniques that continues to enhance our understanding of how chromosomes function. As my PhD project I contributed to this progression by developing in situ hybridization, a technique for hybridization to nucleic acids within their cellular context. Early studies with this technique initiated several lines of research, two of which I describe here, that I have pursued to this day. First, analysis of RNA populations by hybridization to polytene chromosomes (a proto-microarray-type experiment) led us to characterize levels of regulation during heat shock beyond those recognizable by puffing studies. We found also that one still-undeciphered major heat shock puff encodes a novel set of RNAs for which we propose a regulatory role. Second, localization of various multicopy DNA sequences has suggested roles for them in chromosome structure: Most recently we have found that Drosophila telomeres consist of and are maintained by special non-LTR (long terminal repeat) retrotransposons.

The evolutionary significance of cis-regulatory mutations

For decades, evolutionary biologists have argued that changes in cis-regulatory sequences constitute an important part of the genetic basis for adaptation. Although originally based on first principles, this claim is now empirically well supported: numerous studies have identified cis-regulatory mutations with functionally significant consequences for morphology, physiology and behaviour. The focus has now shifted to considering whether cis-regulatory and coding mutations make qualitatively different contributions to phenotypic evolution. Cases in which parallel mutations have produced parallel trait modifications in particular suggest that some phenotypic changes are more likely to result from cis-regulatory mutations than from coding mutations.

Comparative genomic analysis of allatostatin-encoding (Ast) genes in Drosophila species and prediction of regulatory elements by phylogenetic footprinting

The role of the YXFGLa family of allatostatin (AST) peptides in dipterans is not well-established. The recent completion of sequencing of genomes for multiple Drosophila species provides an opportunity to study the evolutionary variation of the allatostatins and to examine regulatory elements that control gene expression. We performed comparative analyses of Ast genes from seven Drosophila species (Drosophila melanogaster, Drosophila simulans, Drosophila ananassae, Drosophila yakuba, Drosophila pseudoobscura, Drosophila mojavensis, and Drosophila grimshawi) and used phylogenetic footprinting methods to identify conserved noncoding motifs, which are candidates for regulatory regions. The peptides encoded by the Ast precursor are nearly identical across species with the exception of AST-1, in which the leading residue may be either methionine or valine. Phylogenetic footprinting predicts as few as 3, to as many as 17 potential regulatory sites depending on the parameters used during analysis. These include a Hunchback motif approximately 1.2 kb upstream of the open reading frame (ORF), overlapping motifs for two Broad-complex isoforms in the first intron, and a CF2-II motif located in the 3'-UTR. Understanding the regulatory elements involved in Ast expression may provide insight into the function of this neuropeptide family.

Exploratory activity in Drosophila requires the kurtz nonvisual arrestin

When Drosophila adults are placed into an open field arena, they initially exhibit an elevated level of activity followed by a reduced stable level of spontaneous activity. We have found that the initial elevated component arises from the fly's interaction with the novel arena since: (1) the increased activity is independent of handling prior to placement within the arena, (2) the fly's elevated activity is proportional to the size of the arena, and (3) the decay in activity to spontaneous levels requires both visual and olfactory input. These data indicate that active exploration is the major component of elevated initial activity. There is a specific requirement for the kurtz nonvisual arrestin in the nervous system for both the exploration stimulated by the novel arena and the mechanically stimulated activity. kurtz is not required for spontaneous activity; kurtz mutants display normal levels of spontaneous activity and average the same velocities as wild-type controls. Inhibition of dopamine signaling has no effect on the elevated initial activity phase in either wild-type or krz(1) mutants. Therefore, the exploratory phase of open field activity requires kurtz in the nervous system, but is independent of dopamine's stimulation of activity.

Drosophila CK2 regulates lateral-inhibition during eye and bristle development

Lateral inhibition is critical for cell fate determination and involves the functions of Notch (N) and its effectors, the Enhancer of Split Complex, E(spl)C repressors. Although E(spl) proteins mediate the repressive effects of N in diverse contexts, the role of phosphorylation was unclear. The studies we describe implicate a common role for the highly conserved Ser/Thr protein kinase CK2 during eye and bristle development. Compromising the functions of the catalytic (alpha) subunit of CK2 elicits a rough eye and defects in the interommatidial bristles (IOBs). These phenotypes are exacerbated by mutations in CK2 and suppressed by an increase in the dosage of this protein kinase. The appearance of the rough eye correlates, in time and space, to the specification and refinement of the 'founding' R8 photoreceptor. Consistent with this observation, compromising CK2 elicits supernumerary R8's at the posterior margin of the morphogenetic furrow (MF), a phenotype characteristic of loss of E(spl)C and impaired lateral inhibition. We also show that compromising CK2 elicits ectopic and split bristles. The former reflects the specification of excess bristle SOPs, while the latter suggests roles during asymmetric divisions that drive morphogenesis of this sensory organ. In addition, these phenotypes are exacerbated by mutations in CK2 or E(spl), indicating genetic interactions between these two loci. Given the centrality of E(spl) to the repressive effects of N, our studies suggest conserved roles for this protein kinase during lateral inhibition. Candidates for this regulation are the E(spl) repressors, the terminal effectors of this pathway.

Basal relationships in the Drosophila melanogaster species group

The Drosophila melanogaster species group is a popular model for evolutionary studies due to its morphological and ecological diversity and its inclusion of the model species D. melanogaster. However, phylogenetic relationships among major lineages within this species group remain controversial. In this report, the phylogeny of 10 species representing each of the well-supported monophyletic clades in the melanogaster group was studied using the sequences of 14 loci that together comprise 9493 nucleotide positions. Combined Bayesian analysis using gene-specific substitution models produced a 100% credible set of two trees. In the strict consensus of these trees, the ananassae subgroup branches first in the melanogaster species group, followed by the montium subgroup. The remaining lineages form a monophyletic clade in which D. ficusphila and D. elegans branch first, followed by D. biarmipes, D. eugracilis, and the melanogaster subgroup. This strongly supported phylogeny resolves most basal relationships in the melanogaster species group, and provides a framework that can be extended in the future to encompass more species.

Wingless signaling in a large insect, the blowflyLucilia sericata: A beautiful example of evolutionary developmental biology

Blowflies are the primary facultative agent in causing myiasis of domestic sheep in the whole world and, at the same time, it is an important tool for forensic medicine. Surprisingly, and in contrast to its importance, almost no data regarding the embryology and molecular markers are known for this insect. In this report, we present a detailed description of the blowfly Lucilia sericata embryogenesis and of imaginal disc development. The embryogenesis of Lucilia strongly resembles that of Drosophila, despite their apparent size difference. Moreover, imaginal disc development appears to be equally well conserved. Through cloning, expression, and functional studies, we show that the Lucilia Wingless (Wg) protein is highly conserved between the two species. We further show that parasegments are established in Lucilia, however, engrailed expression shows a more dynamic expression pattern than expected in comparison to Drosophila. Over-expression of Lucilia Wingless in Drosophila shows wingless-like wing phenotypes, suggesting that Lucilia Wingless blocks the signalling activity of Drosophila Wingless. Upon injection of wg dsRNA, we observe a "lawn of denticle" phenotype, closely resembling that of Drosophila. Due to the large size of the insect, the distance over which Wingless exerts signalling activity is up to three times larger than in Drosophila, yet the consequences are very similar. Our data demonstrate long-range wingless signaling mechanisms adapted for patterning large domains of naked cuticle and suggest signaling properties of Lucilia Wingless that are distinct from those of Drosophila Wingless.

Two retrotransposons maintain telomeres in Drosophila

Telomeres across the genus Drosophila are maintained, not by telomerase, but by two non-LTR retrotransposons, HeT-A and TART, that transpose specifically to chromosome ends. Successive transpositions result in long head-to-tail arrays of these elements. Thus Drosophila telomeres, like those produced by telomerase, consist of repeated sequences reverse transcribed from RNA templates. The Drosophila repeats, complete and 5'-truncated copies of HeT-A and TART, are more complex than telomerase repeats; nevertheless, these evolutionary variants have functional similarities to the more common telomeres. Like other telomeres, the Drosophila arrays are dynamic, fluctuating around an average length that can be changed by changes in the genetic background. Several proteins that interact with telomeres in other species have been found to have homologues that interact with Drosophila telomeres. Although they have hallmarks of non-LTR retrotransposons, HeT-A and TART appear to have a special relationship to Drosophila. Their Gag proteins are efficiently transported into diploid nuclei where HeT-A Gag recruits TART Gag to chromosome ends. Gags of other non-LTR elements remain predominantly in the cytoplasm. These studies provide intriguing evolutionary links between telomeres and retrotransposable elements.

Mitochondrial voltage-dependent anion channel gene family in Drosophila melanogaster: complex patterns of evolution, genomic organization, and developmental expression

Voltage-dependent anion channels (VDACs), also known as mitochondrial porins, are a family of small pore-forming proteins of the mitochondrial outer membrane found in all eukaryotes. VDACs play important roles in the regulated flux of metabolites between the cytosolic and mitochondrial compartments, energy metabolism, and apoptosis. Annotation of the genome sequence of Drosophila melanogaster revealed three genes (CG17137, CG31722-A, and CG31722-B) with homology to porin, the previously described Drosophila VDAC. Molecular analysis reveals a complex pattern of organization and expression. The genomic organization of these four genes and sequence comparisons with other insect VDAC homologs indicate that this gene family evolved through a mechanism of duplication and divergence from an ancestral VDAC gene during the radiation of the genus Drosophila. CG17137, CG31722-A, and CG31722-B are expressed in a male-specific pattern on both transcriptional and translational levels, while porin is equally expressed in both male and female flies. Additionally, CG31722-A and CG31722-B are expressed as a dicistronic transcript. Western blot analysis and immunofluorescence microscopy confirm that these proteins localize to the mitochondrion. Further expression analysis showed that CG17137 and CG31722-B are abundant in testes, while porin is ubiquitously expressed. While porin, CG17137, and CG31722-B are expressed to different degrees during embryogenesis, all of these proteins are dramatically reduced relative to cytochrome c content during larvogenesis. These studies illustrate a complex genomic organization and spatiotemporal pattern of expression for Drosophila VDACs as well as an evolutionary history consistent with either a partitioning of VDAC functions or an acquisition of novel functions among isoforms.

Evolutionary conservation of Ceratitis capitata transformer gene function

Transformer functions as a binary switch gene in the sex determination and sexual differentiation of Drosophila melanogaster and Ceratitis capitata, two insect species that separated nearly 100 million years ago. The TRA protein is required for female differentiation of XX individuals, while XY individuals express smaller, presumably nonfunctional TRA peptides and consequently develop into adult males. In both species, tra confers female sexual identity through a well-conserved double-sex gene. However, unlike Drosophila tra, which is regulated by the upstream Sex-lethal gene, Ceratitis tra itself is likely to control a feedback loop that ensures the maintenance of the female sexual state. The putative CcTRA protein shares a very low degree of sequence identity with the TRA proteins from Drosophila species. However, in this study we show that a female-specific Ceratitis Cctra cDNA encoding the putative full-length CcTRA protein is able to support the female somatic and germline sexual differentiation of D. melanogaster XX; tra mutant adults. Although highly divergent, CcTRA can functionally substitute for DmTRA and induce the female-specific expression of both Dmdsx and Dmfru genes. These data demonstrate the unusual plasticity of the TRA protein that retains a conserved function despite the high evolutionary rate. We suggest that transformer plays an important role in providing a molecular basis for the variety of sex-determining systems seen among insects.

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.

Conservation of an inhibitor of the epidermal growth factor receptor, Kekkon1, in dipterans

Regulation of epidermal growth factor receptor (EGFR) signaling requires the concerted action of both positive and negative factors. While the existence of numerous molecules that stimulate EGFR activity has been well documented, direct biological inhibitors appear to be more limited in number and phylogenetic distribution. Kekkon1 (Kek1) represents one such inhibitor. Kek1 was initially identified in Drosophila melanogaster and appears to be absent from vertebrates and the invertebrate Caenorhabditis. To further investigate Kek1's function and evolution, we identified kek1 orthologs within dipterans. In D. melanogaster, kek1 is a transcriptional target of EGFR signaling during oogenesis, where it acts to attenuate receptor activity through an inhibitory feedback loop. The extracellular and transmembrane portion of Kek1 is sufficient for its inhibitory activity in D. melanogaster. Consistent with conservation of its role in EGFR signaling, interspecies comparisons indicate a high degree of identity throughout these regions. During formation of the dorsal-ventral axis Kek1 is expressed in dorsal follicle cells in a pattern that reflects the profile of receptor activation. D. virilis Kek1 (DvKek1) is also expressed dynamically in the dorsal follicle cells, supporting a conserved role in EGFR signaling. Confirming this, biochemical and transgenic assays indicate that DvKek1 is functionally interchangeable with DmKek1. Strikingly, we find that the cytoplasmic domain contains a region with the highest degree of conservation, which we have implicated in EGFR inhibition and dubbed the Kek tail (KT) box.

A highly conserved candidate chemoreceptor expressed in both olfactory and gustatory tissues in the malaria vector Anopheles gambiae

Anopheles gambiae is a highly anthropophilic mosquito responsible for the majority of malaria transmission in Africa. The biting and host preference behavior of this disease vector is largely influenced by its sense of smell, which is presumably facilitated by G protein-coupled receptor signaling [Takken, W. & Knols, B. (1999) Annu. Rev. Entomol. 44, 131-157]. Because of the importance of host preference to the mosquitoes' ability to transmit disease, we have initiated studies intended to elucidate the molecular mechanisms underlying olfaction in An. gambiae. In the course of these studies, we have identified a number of genes potentially involved in signal transduction, including a family of candidate odorant receptors. One of these receptors, encoded by GPRor7 (hereafter referred to as AgOr7), is remarkably similar to an odorant receptor that is expressed broadly in olfactory tissues and has been identified in Drosophila melanogaster and other insects [Krieger, J., Klink, O., Mohl, C., Raming, K. & Breer, H. (2003) J. Comp. Physiol. A 189, 519-526; Vosshall, L. B., Amrein, H., Morozov, P. S., Rzhetsky, A. & Axel, R. (1999) Cell 96, 725-736]. We have observed AgOr7 expression in olfactory and gustatory tissues in adult An. gambiae and during several stages of the mosquitoes' development. Within the female adult peripheral chemosensory system, antiserum against the AgOR7 polypeptide labels most sensilla of the antenna and maxillary palp as well as a subset of proboscis sensilla. Furthermore, AgOR7 antiserum labeling is observed within the larval antenna and maxillary palpus. These results are consistent with a role for AgOr7 in both olfaction and gustation in An. gambiae and raise the possibility that AgOr7 orthologs may also be of general importance to both modalities of chemosensation in other insects.

Exploring the evolutionary history of the alcohol dehydrogenase gene (Adh) duplication in species of the family tephritidae

In the olive fruit fly Bactrocera oleae and the med fly Ceratitis capitata previous studies have shown the existence of two Adh genes in each species. This observation, in combination with the former finding that various Drosophila species of virilis and repleta group encode two isozymes of ADH which are the result of a gene duplication, challenged us to address a scenario dealing with the evolutionary history of the Adh gene duplication in Tephritidae. In our lab we proceeded to the cloning and sequence analysis of Adh genes from more tephritid species, a prerequisite for further study of this issue. Here we show that phylogenetic trees produced from either the nucleotide or the amino acid sequences of 14 tephritid Adh genes consisted of two main clusters, with Adh sequences of the same "type" grouping together (i.e., Adh1 sequences form a cluster and Adh2 sequences form a second one), as expected if there was one duplication event before speciation within the family Tephritidae. We used the amount of divergence between the two isozymic forms of Adh of the species carrying both Adh1 and Adh2 genes to obtain an estimate of the age of the duplication event. Interestingly, our data again support the hypothesis that the duplication of an ancestral Adh single gene in the family Tephritidae occurred before the emergence of the genera Bactrocera and Ceratitis, thus suggesting that Adh duplication was based on a prespeciation rather than a postspeciation event that might have involved two independent duplication events, one in each of the two genera.

Bifunctional killing activity encoded by conserved reaper proteins

Drosophila activators of apoptosis mapping to the Reaper region function, in part, by antagonizing IAP proteins through a shared RHG motif. We isolated Reaper from the Blowfly L. cuprina, which triggered extensive apoptosis in Drosophila cells. Conserved regions of Reaper were tested in the context of GFP fusions and a second killing activity, distinct from the RHG, was identified. A 20 amino-acid peptide, designated R3, conferred targeting to a focal compartment and promoted membrane blebbing. Killing by the R3 fragment did not correlate with translational suppression or with reduced DIAP1 levels. Likewise, R3-induced cell deaths were only modestly suppressed by silencing of Dronc and involved no detectable association with DIAP1. Instead, a second IAP-binding domain, distinct from the R3, was identified at the C terminus of Reaper that bound to DIAP1 but failed to trigger apoptosis. Collectively, these findings are inconsistent with single effector models for cell killing by Reaper and suggest, instead, that Reaper encodes conserved bifunctional death activities that propagate through distinct effector pathways.

Retrotransposons provide an evolutionarily robust non-telomerase mechanism to maintain telomeres

Telomere molecular biology is far more complex than originally thought. Understanding biological systems is aided by study of evolutionary variants, and Drosophila telomeres are remarkable variants. Drosophila lack telomerase and the arrays of simple repeats generated by telomerase in almost all other organisms; instead, Drosophila telomeres are long tandem arrays of two non-LTR retrotransposons, HeT-A and TART. These are the first transposable elements found to have a bona fide role in cell structure, revealing an unexpected link between telomeres and what is generally considered to be parasitic DNA. In addition to providing insight into the cellular functions performed by telomeres, analysis of HeT-A and TART is providing insight into the evolution of chromosomes, retrotransposons, and retroviruses. Recent studies show that retrotransposon telomeres constitute a robust system for maintaining chromosome ends. These telomeres are now known to predate the separation of extant Drosophila species, allowing ample time for elements and hosts to coevolve interesting mechanisms.

Cytogenetic evidence for asexual evolution of bdelloid rotifers

DNA sequencing has shown individual bdelloid rotifer genomes to contain two or more diverged copies of every gene examined and has revealed no closely similar copies. These and other findings are consistent with long-term asexual evolution of bdelloids. It is not entirely ruled out, however, that bdelloid genomes consist of previously undetected pairs of sequences so similar as to be identical over the regions sequenced, as might result if bdelloids were highly inbred sexual diploids or polyploids. Here, we employ fluorescent in situ hybridization with cosmid probes to determine the copy number and chromosomal distribution of the heat shock gene hsp82 and adjacent sequences in the bdelloid Philodina roseola. We conclude that the four copies identified by sequencing are the only ones present and that each is on a separate chromosome. Bdelloids therefore are not highly homozygous sexually reproducing diploids or polyploids.

Recruitment of the Proneural Gene scute to the Drosophila Sex-Determination Pathway

In flies, scute (sc) works with its paralogs in the achaete-scute-complex (ASC) to direct neuronal development. However, in the family Drosophilidae, sc also acquired a role in the primary event of sex determination, X chromosome counting, by becoming an X chromosome signal element (XSE)—an evolutionary step shown here to have occurred after sc diverged from its closest paralog, achaete (ac). Two temperature-sensitive alleles, scsisB2 and scsisB3, which disrupt only sex determination, were recovered in a powerful F1 genetic selection and used to investigate how sc was recruited to the sex-determination pathway. scsisB2 revealed 3′ nontranscribed regulatory sequences likely to be involved. The scsisB2 lesion abolished XSE activity when combined with mutations engineered in a sequence upstream of all XSEs. In contrast, changes in Sc protein sequence seem not to have been important for recruitment. The observation that the other new allele, scsisB3, eliminates the C-terminal half of Sc without affecting neurogenesis and that scsisB1, the most XSE-specific allele previously available, is a nonsense mutant, would seem to suggest the opposite, but we show that housefly Sc can substitute for fruit fly Sc in sex determination, despite lacking Drosophilidae-specific conserved residues in its C-terminal half. Lack of synergistic lethality among mutations in sc, twist, and dorsal argue against a proposed role for sc in mesoderm formation that had seemed potentially relevant to sex-pathway recruitment. The screen that yielded new sc alleles also generated autosomal duplications that argue against the textbook view that fruit fly sex signal evolution recruited a set of autosomal signal elements comparable to the XSEs.

HeT-A elements in Drosophila virilis: retrotransposon telomeres are conserved across the Drosophila genus

Drosophila melanogaster telomeres are composed of two retrotransposons, HeT-A and TART. Drosophila virilis has recently been shown to have telomere-specific TART elements with many of the characteristics of their D. melanogaster homologues. We now report identification of the second telomere-specific retrotransposon, HeT-A, from D. virilis. These results show that HeT-A and TART have been maintaining telomeres in Drosophila for more than the 60 million years that separate D. melanogaster and D. virilis. All Drosophila species and stocks studied have both of these telomeric elements, suggesting that the elements collaborate, an assumption supported by evidence from D. melanogaster that their Gag proteins interact. Although the HeT-A sequence evolves at a high rate, the element retains the unusual structural features that characterize all HeT-A homologues. These features may be involved in the role of HeT-A at the telomere. The Gag protein from HeT-Avir is as much like TART Gag from other species as it is like HeT-A Gag, suggesting that these Gags are evolving under similar constraints, probably to maintain appropriate interactions with host telomeres and possibly to allow collaborative interactions like those seen in D. melanogaster. In addition, we have identified a chimeric element, Uvir, carrying a pol coding sequence only distantly related to sequences thus far found in any telomere arrays.

The evolution of enzyme specificity in Fasciola spp

Fasciola spp., commonly known as liver fluke, are significant trematode parasites of livestock and humans. They secrete several cathepsin L-like cysteine proteases, some of which differ in enzymatic properties and timing of expression in the parasite's life cycle. A detailed sequence and evolutionary analysis is presented, based on 18 cathepsin L-like enzymes isolated from Fasciola spp. (including a novel clone identified in this study). The enzymes form a monophyletic group which has experienced several gene duplication events over the last approximately 135 million years, giving rise to the present-day enzymatic repertoire of the parasite. This timing of these duplications appears to correlate with important points in the evolution of the mammalian hosts. Furthermore, the dates suggest that Fasciola hepatica and Fasciola gigantica diverged around 19 million years ago. A novel analysis, based on the pattern of amino acid diversity, was used to identify sites in the enzyme that are predicted to be subject to positive adaptive evolution. Many of these sites occur within the active site cleft of the enzymes, and hence would be expected to lead to differences in substrate specificity. Using homology modeling, with reference to previously obtained biochemical data, we are able to predict S2 subsite specificity for these enzymes: specifically those that can accommodate bulky hydrophobic residues in the P2 position and those that cannot. A number of other positions subject to evolutionary pressure and potentially significant for enzyme function are also identified, including sites anticipated to diminish cystatin binding affinity.

RNA editing and regulation of Drosophila 4f-rnp expression by sas-10 antisense readthrough mRNA transcripts

We have previously described an example of extensively A-to-G edited cDNA derived from adult heads of the fruitfly Drosophila melanogaster. In that study, the source of the predicted antisense RNA pairing strand for template recognition by dADAR editase was not identified, and the biological significance of the observed hyperediting was not known. Here, we address each of these questions. 4f-rnp and sas-10 are closely adjacent X-linked genes located on opposite DNA strands that produce convergent transcripts. We show that developmentally regulated antisense sas-10 readthrough mRNA arises by activation of an upstream promoter P2 during the late embryo stage of fly development. The sas-10 readthrough transcripts pair with 4f-rnp mRNA to form double-stranded molecules, as indicated by A-to-G editing observed in both RNA strands. It would be predicted that perfect RNA duplexes would be targeted for modification/degradation by enzyme pathways that recognize double-stranded RNAs, leading to decline in 4f-rnp mRNA levels, and this is what we observe. The observation using quantitative RT-PCR that sas-10 readthrough and 4f-rnp transcript levels are inversely related suggests a role for the antisense RNA in posttranscriptional regulation of 4f-rnp gene expression during development. Potential molecular mechanisms that could lead to this result are discussed, one of which is targeted transcript degradation via the RNAi pathway. Insofar as the dADAR editase and RNAi pathways are known to be constitutive in this system, it is likely that control of antisense RNA transcription is the rate-limiting factor. The results provide insight into roles of naturally occurring antisense RNAs in regulation of eukaryotic gene expression.

The Drosophila hairy RNA localization signal modulates the kinetics of cytoplasmic mRNA transport

In several Drosophila cell types, mRNA transport depends on microtubules, the molecular motor dynein and trans-acting factors including Egalitarian and Bicaudal-D. However, the molecular basis of transcript recognition by the localization machinery is poorly understood. Here, we characterize the features of hairy pair-rule RNA transcripts that mediate their apical localization, using in vivo injection of fluorescently labelled mRNAs into syncytial blastoderm embryos. We show that a 121-nucleotide element within the 3'-untranslated region is necessary and sufficient to mediate apical transport. The signal comprises two essential stem-loop structures, in which double-stranded stems are crucial for localization. Base-pair identities within the stems are not essential, but can contribute to the efficiency of localization, suggesting that specificity is mediated by higher-order structure. Using time-lapse microscopy, we measure the kinetics of localization and show that impaired localization of mutant signals is due to delayed formation of active motor complexes and, unexpectedly, to slower movement. These findings, and those from co-injecting wild-type and mutant RNAs, suggest that the efficiency of molecular motors is modulated by the character of their cargoes.

Transcription control of a gene for Drosophila transcription factor, DREF by DRE and cis-elements conserved between Drosophila melanogaster and virilis

A DNA replication-related element (DRE)-binding factor (DREF) has been revealed to be an important transcription factor for activating promoters of cell proliferation and differentiation related genes. The amino acid sequences of DREF are conserved in evolutionary separate Drosophila species, Drosophila melanogaster (Dm) and Drosophila virilis (Dv) in three regions. In the present study, evidence was obtained that there are several highly conserved regions in the 5' flanking region between the DmDREF and DvDREF genes. Band mobility shift assays using oligonucleotides corresponding to these conserved regions revealed that specific trans-acting factors can bind to at least three regions -554 to -543 (5'-TTTGTTCTTGCG), -81 to -70 (5'-GCCCACGTGGCT) and +225 to +234 (5'-GCAATCAGTG). Using a transient luciferase expression assay, we demonstrated that the region -554 to -543 functions as a negative regulatory element for DmDREF promoter activity, while the regions -77 to -70 (5'-ACGTGGCT) and +225 to +236 (5'-GCAATCAGTGTT) function as positive regulatory elements. In previous studies, we observed that expression of the homeodomain protein Zerknüllt (Zen) represses PCNA gene transcription, by reducing the DNA binding activity of DREF. Here we show Zen downregulates DREF gene promoter activity through action on the region between +241 and +254 (5'-AGAATACTCAACA). In addition, the DmDREF promoter contains five DREs. Using a double stranded RNA-mediated interference method, we generated evidence that expression of DmDREF could be auto-regulated by DREF through the third DRE located at +211 to +218. In living flies we obtained results consistent with those obtained in vitro and in cultured cells. The study thus indicates that DmDREF is effectively regulated via highly conserved regions between the DmDREF and DvDREF promoters, suggesting the existence of common regulatory factors, and that DmDREF can be positively regulated by itself via the third DRE located in its most highly conserved region.

The Musca domestica larval hexamerin is composed of multiple, similar polypeptides

The Musca domestica larval hexamerin (MdHex-L) is a hexameric glycoprotein with an apparent native molecular weight of 500 kDa. Seven different cDNAs that encode MdHex-L subunits were cloned and sequenced. Furthermore, amino acid sequences of isolated subunits were determined by the Edman degradation method and compared to the conceptual translation products derived from the cloned cDNAs. The obtained data indicate the existence of multiple forms of MdHex-L subunits and that these multiple forms may be grouped into three categories according to their percentages of nucleotide sequence identity.

Transposon telomeres are widely distributed in the Drosophila genus: TART elements in the virilis group

Telomeres of most animals, plants, and unicellular eukaryotes are made up of tandem arrays of repeated DNA sequences produced by the enzyme telomerase. Drosophila melanogaster has an unusual variation on this theme; telomeres consist of tandem arrays of sequences produced by successive transpositions of two non-LTR retrotransposons, HeT-A and TART. To explore the phylogenetic distribution of these variant telomeres, we have looked for TART homologues in a distantly related Drosophila species, virilis. We have found elements that, despite many differences in nucleotide sequence, retain significant amino acid similarity to TART from D. melanogaster. These D. virilis TART elements have features that characterize TART elements in D. melanogaster: (i) they are found in tandem arrays on chromosome ends, (ii) they are not found in euchromatin, and (iii) they produce both sense and antisense transcripts, with the antisense RNA being in excess. The D. virilis TART elements have one surprising feature: both of the ORFs contain long stretches of the trinucleotide repeat CAX, encoding polyglutamine (with a few interspersed histidines). These long polyglutamine stretches are conserved in the three D. virilis elements sequenced. They do not interrupt any domains of known function in the TART proteins and are not seen in TART proteins from other species. Comparison of the D. virilis and D. melanogaster telomeres suggests that the retrotransposon mechanism of telomere maintenance may have arisen before the separation of the genus Drosophila.

Transcriptional heterochrony of scute and changes in bristle pattern between two closely related species of blowfly

Temporal shifts in the expression of regulatory genes, relative to other events taking place during development, can result in changes in morphology. Such transcriptional heterochrony can introduce dramatic morphological changes that involve rather few genetic events and so has the potential to cause rapid changes during evolution. We have shown previously that stereotyped species-specific bristle patterns on the notum of higher Diptera correlate with changes in the spatial regulation of scute expression. scute encodes a proneural gene required for the development of sensory bristle precursors and is expressed before pupation in discrete domains on the presumptive notum at sites where the macrochaete precursors arise. Thus, for Ceratitis capitata and Calliphora vicina, species separated from Drosophila melanogaster by about 80 and 100 million years respectively, the domains of sc expression differ. In all three species, a second phase of ubiquitous sc expression, after pupation, precedes formation of the microchaete precursors. Here, we describe sc expression in Phormia terranovae, a species belonging to the family Calliphoridae that is closely related to C. vicina. We find that spatial regulation is almost identical between P. terranovae and C. vicina, in spite of their different bristle patterns. The timing of sc expression differs, however, between the two. The first spatially restricted phase of expression is slightly delayed and the second ubiquitous phase remarkably accelerated, such that there is a period of overlap. As a result, the last precursors from the first phase of expression arise at the same time as the first precursors from the second phase of expression and are morphologically indistinguishable from the late-arising microchaetes. These observations illustrate the power of developmental heterochrony in bringing about rapid morphological change.

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.

Default repression and Notch signaling: Hairless acts as an adaptor to recruit the corepressors Groucho and dCtBP to Suppressor of Hairless

The DNA-binding transcription factor Suppressor of Hairless [Su(H)] functions as an activator during Notch (N) pathway signaling, but can act as a repressor in the absence of signaling. Hairless (H), a novel Drosophila protein, binds to Su(H) and has been proposed to antagonize N signaling by inhibiting DNA binding by Su(H). Here we show that, in vitro, H directly binds two corepressor proteins, Groucho (Gro) and dCtBP. Reduction of gro or dCtBP function enhances H mutant phenotypes and suppresses N phenotypes in the adult mechanosensory bristle. This activity of gro is surprising, because it is directed oppositely to its traditionally defined role as a neurogenic gene. We find that Su(H)-H complexes can bind to DNA with high efficiency in vitro. Furthermore, a H-VP16 fusion protein causes dominant-negative phenotypes in vivo, a result consistent with the proposal that H functions in transcriptional repression. Taken together, our findings indicate that "default repression" of N pathway target genes by an unusual adaptor/corepressor complex is essential for proper cell fate specification during Drosophila peripheral nervous system development.

A complex array of DNA-binding proteins required for pairing-sensitive silencing by a polycomb group response element from the Drosophila engrailed gene

Regulatory DNA from the Drosophila gene engrailed causes silencing of a linked reporter gene (mini-white) in transgenic Drosophila. This silencing is strengthened in flies homozygous for the transgene and has been called "pairing-sensitive silencing." The pairing-sensitive silencing activities of a large fragment (2.6 kb) and a small subfragment (181 bp) were explored. Since pairing-sensitive silencing is often associated with Polycomb group response elements (PREs), we tested the activities of each of these engrailed fragments in a construct designed to detect PRE activity in embryos. Both fragments were found to behave as PREs in a bxd-Ubx-lacZ reporter construct, while the larger fragment showed additional silencing capabilities. Using the mini-white reporter gene, a 139-bp minimal pairing-sensitive element (PSE) was defined. DNA mobility-shift assays using Drosophila nuclear extracts suggested that there are eight protein-binding sites within this 139-bp element. Mutational analysis showed that at least five of these sites are important for pairing-sensitive silencing. One of the required sites is for the Polycomb group protein Pleiohomeotic and another is GAGAG, a sequence bound by the proteins GAGA factor and Pipsqueak. The identity of the other proteins is unknown. These data suggest a surprising degree of complexity in the DNA-binding proteins required for PSE function.

Evolution of 5S rRNA gene families in Drosophila

In Drosophila virilis, the three clusters of 5S rRNA genes on chromosome 5 comprise two different gene families (B and C), which differ profoundly in the organization of their spacer sequences. While C-type genes, which are found in two of the clusters, exhibit a true repetitive character, the B-type genes of the third cluster are each embedded in completely different genomic environments. Southern blots of genomic DNA of different D. virilis subspecies, D. hydei and D. melanogaster probed with 5S rRNA gene spacer and coding sequences demonstrate the specificity of C-type sequences for the D. virilis species group. The comparative analysis of flanking sequences of 5S rRNA genes of D. virilis, members of the D. melanogaster species subgroup and of the blowfly Calliphora erythrocephala reveals the existence of conserved sequence motifs both in the 5' upstream and 3' downstream flanking regions. Their possible roles in the control of expression and processing of the 5S rRNA precursor molecule are discussed.

A comparative genomic analysis of two distant diptera, the fruit fly, Drosophila melanogaster, and the malaria mosquito, Anopheles gambiae

Genome evolution entails changes in the DNA sequence of genes and intergenic regions, changes in gene numbers, and also changes in gene order along the chromosomes. Genes are reshuffled by chromosomal rearrangements such as deletions/insertions, inversions, translocations, and transpositions. Here we report a comparative study of genome organization in the main African malaria vector, Anopheles gambiae, relative to the recently determined sequence of the Drosophila melanogaster genome. The ancestral lines of these two dipteran insects are thought to have separated approximately 250 Myr, a long period that makes this genome comparison especially interesting. Sequence comparisons have identified 113 pairs of putative orthologs of the two species. Chromosomal mapping of orthologous genes reveals that each polytene chromosome arm has a homolog in the other species. Between 41% and 73% of the known orthologous genes remain linked in the respective homologous chromosomal arms, with the remainder translocated to various nonhomologous arms. Within homologous arms, gene order is extensively reshuffled, but a limited degree of conserved local synteny (microsynteny) can be recognized.

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.

Interaction of gap genes in the Drosophila head: tailless regulates expression of empty spiracles in early embryonic patterning and brain development

Unlike gap genes in the trunk region of Drosophila embryos, gap genes in the head were presumed not to regulate each other's transcription. Here, we show that in tailless (tll) loss-of-function mutants the empty spiracles (ems) expression domain in the head expands, whereas it retracts in tll gain-of-function embryos. We have identified a 304bp element in the ems-enhancer which is sufficient to drive expression in the head and brain and which contains two TLL and two BCD binding sites. Transgenic reporter gene lines containing mutations of the TLL binding sites demonstrate that tll directly inhibits the expression of ems in the early embryonic head and the protocerebral brain anlage. These results are the first demonstration of direct transcriptional regulation between gap genes in the head.

Role of CtBP in transcriptional repression by the Drosophila giant protein

The giant protein is a short-range transcriptional repressor that refines the expression pattern of gap and pair-rule genes in the Drosophila blastoderm embryo. Short-range repressors including knirps, Krüppel, and snail utilize the CtBP cofactor for repression, but it is not known whether a functional interaction with CtBP is a general property of all short-range repressors. We studied giant repression activity in a CtBP mutant and find that this cofactor is required for giant repression of some, but not all, genes. While targets of giant such as the even-skipped stripe 2 enhancer and a synthetic lacZ reporter show clear derepression in the CtBP mutant, another giant target, the hunchback gene, is expressed normally. A more complex situation is seen with regulation of the Krüppel gene, in which one enhancer is repressed by giant in a CtBP-dependent manner, while another is repressed in a CtBP-independent manner. These results demonstrate that giant can repress both via CtBP-dependent and CtBP-independent pathways, and that promoter context is critical for determining giant-CtBP functional interaction. To initiate mechanistic studies of the giant repression activity, we have identified a minimal repression domain within giant that encompasses residues 89-205, including an evolutionarily conserved region bearing a putative CtBP binding motif.

The functions of the multiproduct and rapidly evolving dec-1 eggshell gene are conserved between evolutionarily distant species of Drosophila

The Drosophila dec-1 gene encodes multiple proteins that are required for female fertility and proper eggshell morphogenesis. Genetic and immunolocalization data suggest that the different DEC-1 proteins are functionally distinct. To identify regions within the proteins with potential biological significance, we cloned and sequenced the D. yakuba and D. virilis dec-1 homologs. Interspecies comparisons of the predicted translation products revealed rapidly evolving sequences punctuated by blocks of conserved amino acids. Despite extensive amino acid variability, the proteins produced by the different dec-1 homologs were functionally interchangeable. The introduction of transgenes containing either the D. yakuba or the D. virilis dec-1 open reading frames into a D. melanogaster DEC-1 protein null mutant was sufficient to restore female fertility and wild-type eggshell morphology. Normal expression and extracellular processing of the DEC-1 proteins was correlated with the phenotypic rescue. The nature of the conserved features highlighted by the evolutionary comparison and the molecular resemblance of some of these features to those found in other extracellular proteins suggests functional correlates for some of the multiple DEC-1 derivatives.

Rapid evolution of a cyclin A inhibitor gene, roughex, in Drosophila

The recent sequencing of the complete genome of the fruit fly Drosophila melanogaster has yielded about 30% of the predicted genes with no obvious counterparts in other organisms. These rapidly evolving genes remain largely unexplored. Here, we present evidence for a striking variability in an important Drosophila cell cycle regulator encoded by the gene roughex (rux) in closely related fly species. The unusual level of Rux protein variability indicates that there are very low overall constraints on amino acid substitutions. Despite the lack of sequence similarity, certain common features, including the presence of a C-terminal nuclear localization signal and a functionally important N-terminal RXL cyclin-binding motif, exist between Rux and cyclin-dependent kinase inhibitors of the Cip/Kip family. These results indicate that even some genes involved in key regulatory processes in eukaryotes evolve at extremely high rates.

The Twisted gastrulation family of proteins, together with the IGFBP and CCN families, comprise the TIC superfamily of cysteine rich secreted factors

AIMS

To analyse the similarities between the Twisted gastrulation (TSG) proteins known to date; in addition, to determine phylogenetic relations among the TSG proteins, and between the TSGs and other protein families--the CCN (for example, CCN2 (CTGF), CCN1 (CYR61), and CCN3 (NOV)) and IGFBP (insulin-like growth factor binding protein) families.

METHODS

TBLASTN and FASTA3 were used to identify new tsg genes and relatives of the TSG family. The sequences were aligned with ClustalW. The predictions of sites for signal peptide cleavage, post-translational modifications, and putative protein domains were carried out with software available at various databases. Unrooted phylogenetic trees were calculated using the UPGMA method.

RESULTS

Several tsg genes from vertebrates and invertebrates were compared. Alignment of protein sequences revealed a highly conserved family of TSG proteins present in both vertebrates and invertebrates, whereas the slightly less well conserved IGFBP and CCN proteins are apparently present only in vertebrates. The TSG proteins display strong homology among themselves and they are composed of a putative signal peptide at the N-terminus followed by a cysteine rich (CR) region, a conserved domain devoid of cysteines, a variable midregion, and a C-terminal CR region. The most striking similarity between the TSGs and the IGFBP and CCN proteins occurs in the N-terminal conserved cysteine rich domain and the characteristic 5' cysteine rich domain(s), spacer region, and 3' cysteine rich domain structure.

CONCLUSION

The family of highly conserved TSG proteins, together with the IGFBP and CCN families, constitute an emerging multigene superfamily of secreted cysteine rich factors. The TSG branch of the superfamily appears to pre-date the others because it is present in all species examined, whereas the CCN and IGFBP genes are found only in vertebrates.

Salivary gland determination in Drosophila: a salivary-specific, fork head enhancer integrates spatial pattern and allows fork head autoregulation

In the early Drosophila embryo, a system of coordinates is laid down by segmentation genes and dorsoventral patterning genes. Subsequently, these coordinates must be interpreted to define particular tissues and organs. To begin understanding this process for a single organ, we have studied how one of the first salivary gland genes, fork head (fkh), is turned on in the primordium of this organ, the salivary placode. A placode-specific fkh enhancer was identified 10 kb from the coding sequence. Dissection of this enhancer showed that the apparently homogeneous placode is actually composed of at least four overlapping domains. These domains appear to be developmentally important because they predict the order of salivary invagination, are evolutionarily conserved, and are regulated by patterning genes that are important for salivary development. Three dorsoventral domains are defined by EGF receptor (EGFR) signaling, while stripes located at the anterior and posterior edges of the placode depend on wingless signaling. Further analysis identified sites in the enhancer that respond either positively to the primary activator of salivary gland genes, SEX COMBS REDUCED (SCR), or negatively to EGFR signaling. These results show that fkh integrates spatial pattern directly, without reference to other early salivary gland genes. In addition, we identified a binding site for FKH protein that appears to act in fkh autoregulation, keeping the gene active after SCR has disappeared from the placode. This autoregulation may explain how the salivary gland maintains its identity after the organ is established. Although the fkh enhancer integrates information needed to define the salivary placode, and although fkh mutants have the most extreme effects on salivary gland development thus far described, we argue that fkh is not a selector gene for salivary gland development and that there is no master, salivary gland selector gene. Instead, several genes independently sense spatial information and cooperate to define the salivary placode.

Isolation and characterization of the Xanthine dehydrogenase gene of the Mediterranean fruit fly, Ceratitis capitata

Xanthine dehydrogenase (XDH) is a member of the molybdenum hydroxylase family of enzymes catalyzing the oxidation of hypoxanthine and xanthine to uric acid. The enzyme is also required for the production of one of the major Drosophila eye pigments, drosopterin. The XDH gene has been isolated in many species representing a broad cross section of the major groups of living organisms, including the cDNA encoding XDH from the Mediterranean fruit fly Ceratitis capitata (CcXDH) described here. CcXDH is closely related to other insect XDHs and is able to rescue the phenotype of the Drosophila melanogaster XDH mutant, rosy, in germline transformation experiments. A previously identified medfly mutant, termed rosy, whose phenotype is suggestive of a disruption in XDH function, has been examined for possible mutations in the XDH gene. However, we find no direct evidence that a mutation in the CcXDH gene or that a reduction in the CcXDH enzyme activity is present in rosy medflies. Conclusive studies of the nature of the medfly rosy mutant will require rescue by germline transformation of mutant medflies.

Characterization of muscle actin genes in Drosophila virilis reveals significant molecular complexity in skeletal muscle types

Actin is a ubiquitous and highly conserved eukaryotic protein required for cell motility and locomotion. In this manuscript, we characterize the four muscle actin genes of the insect Drosophila virilis and demonstrate strong similarities between the D. virilis genes and their homologues in Drosophila melanogaster; intron locations are conserved, and there are few amino acid differences between homologues. We also found strong conservation in temporal expression patterns of the muscle actin genes--the homologues of the D. melanogaster genes Act57B and Act87E are expressed throughout the life cycle, whereas the other two D. virilis genes, homologous to Act79B and Act88F are specific to pupal and adult stages. In situ hybridization revealed that each D. virilis gene is expressed in a unique pattern in the muscles of the thorax and abdomen. These muscle-specific patterns of actin isoforms suggest a greater physiological diversity for the adult muscles of insects than has been appreciated to date from their categorization into fibrillar, tubular (non-fibrillar) and supercontractile muscle types.

Divergent structure and function of the bicoid gene in Muscoidea fly species

We have investigated the evolution of the bicoid (bcd) gene in fly species of the Muscoidea Superfamily. We obtained the complete bcd sequence from the housefly Musca domestica and found polymorphism in the coding region among Musca strains. In addition to Musca, we cloned most of the bcd coding sequences from two blowfly species Calliphora vicina and Lucilia sericata. The 5' and 3' regulatory regions flanking the Musca bcd gene are widely diverged in sequence from Drosophila; however, some important sequence motifs identified in Drosophila bcd are present. The predicted RNA secondary structures of the 3' UTRs are similar, despite sequence divergence. Comparison of Bicoid (Bcd) proteins shows a serine-rich domain of unknown function is present in the Muscoidea species, but is absent in other species. The in vivo function of bcd in Musca was tested by RNAi to mimic loss of function phenotype. We obtained a head defect phenotype similar to weak bcd alleles of Drosophila. Although our comparisons initially suggest functional conservation between species, closer inspection reveals significant differences. Divergence of structural motifs, such as regulatory elements in flanking regions and conservation of protein domains in some species but not in others, points to functional divergence between species. We suggest that the larger embryonic size in Muscoidea species restricts the morphogenetic activity of a weak Bcd activator, which has evolved a more specialized role in head determination and lost some functions in thoracic development.

Spatial and temporal control of RNA stability

Maternally encoded RNAs and proteins program the early development of all animals. A subset of the maternal transcripts is eliminated from the embryo before the midblastula transition. In certain cases, transcripts are protected from degradation in a subregion of the embryonic cytoplasm, thus resulting in transcript localization. Maternal factors are sufficient for both the degradation and protection components of transcript localization. Cis-acting elements in the RNAs convert transcripts progressively (i) from inherently stable to unstable and (ii) from uniformly degraded to locally protected. Similar mechanisms are likely to act later in development to restrict certain classes of transcripts to particular cell types within somatic cell lineages. Functions of transcript degradation and protection are discussed.

Phylogenetic utility of mitochondrial COI and nuclear Gpdh genes in Drosophila

Phylogenetic utility of the mitochondrial COI (cytochrome oxidase subunit I) and nuclear Gpdh (glycerol-3-phosphate dehydrogenase) genes was studied in the Drosophila melanogaster species group. The rate of substitution was higher in the COI gene than in the Gpdh gene. In addition, multiple substitutions, not only for transitional but also for transversional substitutions, occurred faster in the COI gene. None of the trees obtained using the COI gene supported the well-established monophyly of the ananassae subgroup. In addition, the incongruence length difference test, Templeton test, and partitioned Bremer support revealed that the trees based on the COI data are considerably different from those based on the Gpdh and the combined data set. Thus, the COI gene did not show good phylogenetic performance in the melanogaster group. The present analyses based on the Gpdh gene and the combined data set revealed that the ananassae subgroup branched off first in the melanogaster group followed by the montium subgroup and further by the melanogaster subgroup in contrast to the most recent phylogenetic hypothesis based on Amy multigenes.

Molecular organization of the Drosophila melanogaster Adh chromosomal region in D. repleta and D. buzzatii, two distantly related species of the Drosophila subgenus

The molecular organization of a 1.944-Mb chromosomal region of Drosophila melanogaster around the Adh locus has been analyzed in two repleta group species: D. repleta and D. buzzatii. The extensive genetic and molecular information about this region in D. melanogaster makes it a prime choice for comparative studies of genomic organization among distantly related species. A set of 26 P1 phages from D. melanogaster were successfully hybridized using fluorescence in-situ hybridization (FISH) to the salivary gland chromosomes of both repleta group species. The results show that the Adh region is distributed in D. repleta and D. buzatii over six distant sites of chromosome 3, homologous to chromosomal arm 2L of D. melanogaster (Muller's element B). This observation implies a density of 2.57 fixed breakpoints per Mb in the Adh region and suggests a considerable reorganization of this chromosomal element via the fixation of paracentric inversions. Nevertheless, breakpoint density in the Adh region is three times lower than that estimated for D. repleta chromosome 2, homologous to D. melanogaster 3R (Muller's element E). Differences in the rate of evolution among chromosomal elements are seemingly persistent in the Drosophila genus over long phylogenetic distances.