Sequences In Drosophila Research Articles

Drosophila Polymorphism Database (DPDB)A Portal for Nucleotide Polymorphism in Drosophila

As a growing number of haplotypic sequences from re-sequencing studies are now accumulating for Drosophila in the main primary sequence databases, collectively they can now be used to describe the general pattern of nucleotide variation across species and genes of this genus. The Drosophila Polymorphism Database (DPDB) is a secondary database that provides a collection of all well-annotated polymorphic sequences in Drosophila together with their associated diversity measures and options for re-analysis of the data that greatly facilitate both multi-locus and multi-species diversity studies in one of the most important group of model organisms. Here we describe the state-of-the-art of the DPDB database and provide a step-by-step guide to all its searching and analytic capabilities. Finally, we illustrate its usefulness through selected examples. DPDB is freely available at http://dpdb.uab.cat.

Drosophila Omi, a mitochondrial-localized IAP antagonist and proapoptotic serine protease

Although essential in mammals, in flies the importance of mitochondrial outer membrane permeabilization for apoptosis remains highly controversial. Herein, we demonstrate that Drosophila Omi (dOmi), a fly homologue of the serine protease Omi/HtrA2, is a developmentally regulated mitochondrial intermembrane space protein that undergoes processive cleavage, in situ, to generate two distinct inhibitor of apoptosis (IAP) binding motifs. Depending upon the proapoptotic stimulus, mature dOmi is then differentially released into the cytosol, where it binds selectively to the baculovirus IAP repeat 2 (BIR2) domain in Drosophila IAP1 (DIAP1) and displaces the initiator caspase DRONC. This interaction alone, however, is insufficient to promote apoptosis, as dOmi fails to displace the effector caspase DrICE from the BIR1 domain in DIAP1. Rather, dOmi alleviates DIAP1 inhibition of all caspases by proteolytically degrading DIAP1 and induces apoptosis both in cultured cells and in the developing fly eye. In summary, we demonstrate for the first time in flies that mitochondrial permeabilization not only occurs during apoptosis but also results in the release of a bona fide proapoptotic protein.

Human DNA Replication-related Element Binding Factor (hDREF) Self-association via hATC Domain Is Necessary for Its Nuclear Accumulation and DNA Binding

We previously demonstrated that hDREF, a human homologue of Drosophila DNA replication-related element binding factor (dDREF), is a DNA-binding protein predominantly distributed with granular structures in the nucleus. Here, glutathione S-transferase pulldown and chemical cross-linking assays showed that the carboxyl-terminal hATC domain of hDREF, highly conserved among hAT transposase family members, possesses self-association activity. Immunoprecipitation analyses demonstrated that hDREF self-associates in vivo, dependent on hATC domain. Moreover, analyses using a series of hDREF mutants carrying amino acid substitutions in the hATC domain revealed that conserved hydrophobic amino acids are essential for self-association. Immunofluorescence studies further showed that all hDREF mutants lacking self-association activity failed to accumulate in the nucleus. Self-association-defective hDREF mutants also lost association with endogenous importin beta1. Moreover, electrophoretic gel-mobility shift assays revealed that the mutations completely abolished the DNA binding activity of hDREF. These results suggest that self-association of hDREF via the hATC domain is necessary for its nuclear accumulation and DNA binding. We also found that ZBED4/KIAA0637, another member of the human hAT family, also self-associates, again dependent on the hATC domain, with deletion resulting in loss of efficient nuclear accumulation. Thus, hATC domains of human hAT family members appear to have conserved functions in self-association that are required for nuclear accumulation.

Drosophila retinophilin contains MORN repeats and is conserved in humans

The function of conserved novel human genes can be efficiently addressed in genetic model organisms. From a collection of genes expressed in the Drosophila visual system, cDNAs expressed in vertebrates were identified and one similar to a novel human gene was chosen for further investigation. The results reported here characterize the Drosophila retinophilin gene and demonstrate that a similar gene is expressed in the human retina. The Drosophila and human retinophilin sequences are 50% identical, and they share an additional 16% conserved substitutions. Examination of the cDNA and genomic sequence indicates that it corresponds to the gene CG10233 of the annotated genome and predicts a 22.7 kDa protein. Polyclonal antibodies generated to a predicted retinophilin peptide recognize an antigen in Drosophila photoreceptor cells. The retinophilins encode 4 copies of a repeat associated with a Membrane Occupation and Recognition Nexus (MORN) function first discovered in junctophilins, which may interact with the plasma membrane. These results therefore show that Drosophila retinophilin is expressed in fly photoreceptor cells, demonstrate that a conserved human gene is expressed in human retina, and suggest that a mutational analysis of the Drosophila gene would be valuable.

Expression, Purification, and Characterization of a Predicted Drosophila Glutaminyl Cyclase

Glutaminyl cyclase (QC) is responsible for the N-terminal conversion of glutamine to pyroglutamate in many peptides. Though the N-terminal pyroglutamate is found in many different peptides in virtually all species, including Drosophila, QC activity has not been linked to this peptide modification in Drosophila. Sequence analyses revealed putative protein sequences in Drosophila with greater than 50% sequence homology to human QC. To evaluate the enzymatic activity of the predicted Drosophila QC sequence a bacterial fusion protein was designed which contained the predicted amino acid sequence for Drosophila QC along with a purification tag. Expression of the fusion protein and subsequent purification revealed a protein of the expected size of 44kD. The 44kD protein was evaluated using two previously published QC activity assays. Both assays revealed QC activity in the fusion protein similar to human QC as well similar inhibition with published QC inhibitors suggesting the Drosophila protein is a QC. Further analyses will be completed to compare the structural properties of the Drosophila QC to the human QC.

Indelign: a probabilistic framework for annotation of insertions and deletions in a multiple alignment

A quantitative study of molecular evolutionary events such as substitutions, insertions and deletions from closely related genomes requires (1) an accurate multiple sequence alignment program and (2) a method to annotate the insertions and deletions that explain the 'gaps' in the alignment. Although the former requirement has been extensively addressed, the latter problem has received little attention, especially in a comprehensive probabilistic framework. Here, we present Indelign, a program that uses a probabilistic evolutionary model to compute the most likely scenario of insertions and deletions consistent with an input multiple alignment. It is also capable of modifying the given alignment so as to obtain a better agreement with the evolutionary model. We find close to optimal performance and substantial improvement over alternative methods, in tests of Indelign on synthetic data. We use Indelign to analyze regulatory sequences in Drosophila, and find an excess of insertions over deletions, which is different from what has been reported for neutral sequences. The Indelign program may be downloaded from the website http://veda.cs.uiuc.edu/indelign/ Supplementary material is available at Bioinformatics online.

Detection of a gypsy-like sequence in the genome of the cat flea Ctenocephalides felis (Bouché 1835)

The gypsy virus of Drosophila melanogaster is the first identified retrovirus of insects. Its infectious properties have been reported, and it is able to cause diverse mutations of the phenotype in the fruit fly. As a typical endogenous retrovirus, it is transmitted vertically as part of the host genome, but horizontal transmission has also been shown for the D. melanogaster gypsy virus. Using polymerase chain reaction with degenerate primers designed from the gypsy element of Drosophila subobscura, a gypsy-like sequence was amplified from the genome of the cat flea Ctenocephalides felis in a pool of laboratory and field strains. High identities with other gypsy-like sequences of Drosophila spp. can be detected. In contrast to some of these gypsy sequences, the 410-bp fragment of C. felis possesses two deletions, which are localised in the env gene region. Because of these deletions and the lack of virus particles in C. felis, the criteria for host infectivity are not fulfilled. Nevertheless, the general detection of a gypsy-like sequence in C. felis and the capability of the ectoparasite to act as transmitter of viral diseases can be of importance for the cat flea's position in the epidemiology of diseases, which are accompanied by a higher rate of mutation.

Prices for Ingenuity

Prices for Ingenuity

Cyclic nucleotide phosphodiesterases in Drosophila melanogaster.

Cyclic nucleotide PDEs (phosphodiesterases) are important enzymes that regulate intracellular levels of cAMP and cGMP. In the present study, we identify and characterize novel PDEs in the genetic model, Drosophila melanogaster. The Drosophila genome encodes five novel PDE genes in addition to dunce. Predicted PDE sequences of Drosophila show highly conserved critical domains when compared with human PDEs. Thus PDE-encoding genes of D. melanogaster are CG14940-PDE1C, CG8279-PDE6beta, CG5411-PDE8A, CG32648-PDE9 and CG10231-PDE11. Reverse transcriptase-PCRs of adult tissues reveal widespread expression of PDE genes. Drosophila Malpighian (renal) tubules express all the six PDEs: Drosophila PDE1, dunce (PDE4), PDE6, PDE8, PDE9 and PDE11. Antipeptide antibodies were raised against PDE1, PDE6, PDE9 and PDE11. Verification of antibody specificity by Western blotting of cloned and expressed PDE constructs allowed the immunoprecipitation studies of adult Drosophila lysates. Biochemical characterization of immunoprecipitated endogenous PDEs showed that PDE1 is a dual-specificity PDE (Michaelis constant Km for cGMP: 15.3+/-1 microM; Km cAMP: 20.5+/-1.5 microM), PDE6 is a cGMP-specific PDE (Km cGMP: 37+/-13 microM) and PDE11 is a dual-specificity PDE (Km cGMP: 6+/-2 microM; Km cAMP: 18.5+/-5.5 microM). Drosophila PDE1, PDE6 and PDE11 display sensitivity to vertebrate PDE inhibitors, zaprinast (IC50 was 71+/-39 microM for PDE1, 0.65+/-0.015 microM for PDE6 and 1.6+/-0.5 microM for PDE11) and sildenafil (IC50 was 1.3+/-0.9 microM for PDE1, 0.025+/-0.005 microM for PDE6 and 0.12+/-0.06 microM for PDE11). We provide the first characterization of a cGMP-specific PDE and two dual-specificity PDEs in Drosophila, and show a high degree of similarity in structure and function between human and Drosophila PDEs.

Genomic Heterogeneity of Background Substitutional Patterns in Drosophila melanogaster

Mutation is the underlying force that provides the variation upon which evolutionary forces can act. It is important to understand how mutation rates vary within genomes and how the probabilities of fixation of new mutations vary as well. If substitutional processes across the genome are heterogeneous, then examining patterns of coding sequence evolution without taking these underlying variations into account may be misleading. Here we present the first rigorous test of substitution rate heterogeneity in the Drosophila melanogaster genome using almost 1500 nonfunctional fragments of the transposable element DNAREP1_DM. Not only do our analyses suggest that substitutional patterns in heterochromatic and euchromatic sequences are different, but also they provide support in favor of a recombination-associated substitutional bias toward G and C in this species. The magnitude of this bias is entirely sufficient to explain recombination-associated patterns of codon usage on the autosomes of the D. melanogaster genome. We also document a bias toward lower GC content in the pattern of small insertions and deletions (indels). In addition, the GC content of noncoding DNA in Drosophila is higher than would be predicted on the basis of the pattern of nucleotide substitutions and small indels. However, we argue that the fast turnover of noncoding sequences in Drosophila makes it difficult to assess the importance of the GC biases in nucleotide substitutions and small indels in shaping the base composition of noncoding sequences.

Selection on Codon Usage for Error Minimization at the Protein Level

Given the structure of the genetic code, synonymous codons differ in their capacity to minimize the effects of errors due to mutation or mistranslation. I suggest that this may lead, in protein-coding genes, to a preference for codons that minimize the impact of errors at the protein level. I develop a theoretical measure of error minimization for each codon, based on amino acid similarity. This measure is used to calculate the degree of error minimization for 82 genes of Drosophila melanogaster and 432 rodent genes and to study its relationship with CG content, the degree of codon usage bias, and the rate of nucleotide substitution. I show that (i) Drosophila and rodent genes tend to prefer codons that minimize errors; (ii) this cannot be merely the effect of mutation bias; (iii) the degree of error minimization is correlated with the degree of codon usage bias; (iv) the amino acids that contribute more to codon usage bias are the ones for which synonymous codons differ more in the capacity to minimize errors; and (v) the degree of error minimization is correlated with the rate of nonsynonymous substitution. These results suggest that natural selection for error minimization at the protein level plays a role in the evolution of coding sequences in Drosophila and rodents.

Looking at the bigger picture

Hobert et al. ([Hobert et al., 2004][1]) have made a number of criticisms on our paper ([Vogel et al., 2003][2]). In the following paragraphs we give our replies to these criticisms. In a number of cases, the comments do provide useful corrections to the paper. Nevertheless, the major conclusions of

Evolutionary history and mode of the amylase multigene family in Drosophila.

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 +/- 1.93 and 12.38 +/- 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.

EST Analysis of the Cnidarian Acropora millepora Reveals Extensive Gene Loss and Rapid Sequence Divergence in the Model Invertebrates

A significant proportion of mammalian genes are not represented in the genomes of Drosophila, Caenorhabditis or Saccharomyces, and many of these are assumed to have been vertebrate innovations. To test this assumption, we conducted a preliminary EST project on the anthozoan cnidarian, Acropora millepora, a basal metazoan. More than 10% of the Acropora ESTs with strong metazoan matches to the databases had clear human homologs but were not represented in the Drosophila or Caenorhabditis genomes; this category includes a surprising diversity of transcription factors and metabolic proteins that were previously assumed to be restricted to vertebrates. Consistent with higher rates of divergence in the model invertebrates, three-way comparisons show that most Acropora ESTs match human sequences much more strongly than they do any Drosophila or Caenorhabditis sequence. Gene loss has thus been much more extensive in the model invertebrate lineages than previously assumed and, as a consequence, some genes formerly thought to be vertebrate inventions must have been present in the common metazoan ancestor. The complexity of the Acropora genome is paradoxical, given that this organism contains apparently few tissue types and the simplest extant nervous system consisting of a morphologically homogeneous nerve net.

The advantages of the ITS2 region of the nuclear rDNA cistron for analysis of phylogenetic relationships of insects: a Drosophila example

We examined the utility for phylogenetic reconstruction of the second internal transcribed spacer (ITS2), lying between the nuclear 5.8S gene and the gene for large subunit ribosomal RNA, using sequences of Ceratitis, Bactrocera, Musca, and Drosophila. We aligned and analyzed 13 sequences from GenBank and 11 new sequences from diverse species of Drosophila. Derivation of the secondary structure of the ITS2, the RNA transcript folding pattern required for transcript processing into functional RNA units, revealed the facets of sequence conservation common to all the sequences, that then allowed alignment of all the genera. The resultant tree, though including only a sparse representation of the enormous Drosophila diversity, conforms generally with the consensus of all prior phylogenetic reconstructions, using eight other nuclear and mitochondrial gene regions; where species representation is greater, as in the melanogaster subgroup of the Sophophora subgenus representatives, it conforms exactly. The paradigm ITS2 secondary structure presented can now be used to assess the genus more thoroughly, since its base pairing pattern makes alignment of sequences obvious. In addition, it shows that these insects share the ITS2 secondary structure characteristics of the other major animal groups as well as the green line of eukaryote evolution. The relatively short (<400 bp) ITS2 region seems ideal for reconstructing evolutionary relationships at the levels of species, genera, and perhaps even higher.

The role of DNA sequence in centromere formation

Recent analyses of centromeric sequences in human, Drosophila and Arabidopsis, have provided new insights into the importance of DNA sequence for centromere function.

The regulation of yp3 expression in the Drosophila melanogaster fat body.

The regulation of the Drosophila melanogaster yolk protein genes 1 and 2 have been well characterised. Cis-acting DNA elements and trans-acting factors regulating ovarian fat body and sex-specific expression have been identified. In this paper we have analysed the regulation of yolk protein 3, which is separated from the other two genes on the X-chromosome. We have separated sex-specific control from fat body control in some constructs in transgenic flies. We propose that the organisation of the regulatory elements in yp3 differs from yp1 and yp2 for control of fat body expression and that it closely resembles the regulation of a reporter gene using Musca and Calliphora yp promoter enhancer sequences in transgenic Drosophila.

Human-fly homologies

HOMOPHILA is a searchable database of human disease genes, obtained from the Online Mendelian Inheritance in Man database, that have significant homology at the protein level to genes contained in FlyBase, the current Drosophila sequence database.

Identification of the β‐dystroglycan binding epitope within the C‐terminal region of α‐dystroglycan

Dystroglycan is a receptor for extracellular matrix proteins that plays a crucial role during embryogenesis in addition to adult tissue stabilization. A precursor product of a single gene is post-translationally cleaved to form two different subunits, alpha and beta. The extracellular alpha-dystroglycan is a membrane-associated, highly glycosylated protein that binds to various extracellular matrix molecules, whereas the transmembrane beta-dystroglycan binds, via its cytosolic domain, to dystrophin and many other proteins. alpha- and beta-Dystroglycan interact tightly but noncovalently. We have previously shown that the N-terminal region of beta-dystroglycan, beta-DG(654-750), binds to the C-terminal region of murine alpha-dystroglycan independently from glycosylation. Preparing a series of deleted recombinant fragments and using solid-phase binding assays, the C-terminal sequence of alpha-dystroglycan containing the binding epitope for beta-dystroglycan has been defined more precisely. We found that a region of 36 amino acids, from position 550-585, is required for binding the extracellular region, amino acids 654-750 of beta-dystroglycan. Recently, a dystroglycan-like gene was identified in Drosophila that showed a moderate degree of conservation with vertebrate dystroglycan (31% identity, 48% similarity). Surprisingly, the Drosophila sequence contains a region showing a higher degree of identity and conservation (45% and 66%) that coincides with the 550-585 sequence of vertebrate alpha-dystroglycan. We have expressed this Drosophila dystroglycan fragment and measured its binding to the extracellular region of vertebrate (murine) beta-dystroglycan (Kd = 6 +/- 1 microM). These data confirm the proper identification of the beta-dystroglycan binding epitope and stress the importance of this region during evolution. This finding might help the rational design of dystroglycan-specific binding drugs, that could have important biomedical applications.

Transcriptional repression by Pax5 (BSAP) through interaction with corepressors of the Groucho family.

Pax5 (BSAP) functions as both a transcriptional activator and repressor during midbrain patterning, B-cell development and lymphomagenesis. Here we demonstrate that Pax5 exerts its repression function by recruiting members of the Groucho corepressor family. In a yeast two-hybrid screen, the groucho-related gene product Grg4 was identified as a Pax5 partner protein. Both proteins interact cooperatively via two separate domains: the N-terminal Q and central SP regions of Grg4, and the octapeptide motif and C-terminal transactivation domain of Pax5. The phosphorylation state of Grg4 is altered in vivo upon Pax5 binding. Moreover, Grg4 efficiently represses the transcriptional activity of Pax5 in an octapeptide-dependent manner. Similar protein interactions resulting in transcriptional repression were also observed between distantly related members of both the Pax2/5/8 and Groucho protein families. In agreement with this evolutionary conservation, the octapeptide motif of Pax proteins functions as a Groucho-dependent repression domain in Drosophila embryos. These data indicate that Pax proteins can be converted from transcriptional activators to repressors through interaction with corepressors of the Groucho protein family.