Melanogaster Genome Research Articles

Discovering non-coding RNA elements in Drosophila 3' untranslated regions.

The Non-Coding RNA (ncRNA) elements in the 3' Untranslated Regions (3'-UTRs) are known to participate in the genes' post-transcriptional regulations. Inferring co-expression patterns of the genes through clustering these 3'-UTR ncRNA elements will provide invaluable insights for studying their biological functions. In this paper, we propose an improved RNA structural clustering pipeline. Benchmark of the new pipeline on Rfam data demonstrates over 10% performance improvements compared to the traditional hierarchical clustering pipeline. By applying the new clustering pipeline to 3'-UTRs of Drosophila melanogaster's genome, we have successfully identified 184 ncRNA clusters with 91.3% accuracy. One of these clusters corresponds to genes that are preferentially expressed in male Drosophila. Another cluster contains genes that are responsible for the functions of septate junction in epithelial cells. These discoveries encourage more studies on novel post-transcriptional regulation mechanisms.

Inducible Protein Traps with Dominant Phenotypes for Functional Analysis of theDrosophilaGenome

The Drosophila melanogaster genome has been extensively characterized, but there remains a pressing need to associate gene products with phenotypes, subcellular localizations, and interaction partners. A multifunctional, Minos transposon-based protein trapping system called Hostile takeover (Hto) was developed to facilitate in vivo analyses of endogenous genes, including live imaging, purification of protein complexes, and mutagenesis. The Hto transposon features a UAS enhancer with a basal promoter, followed by an artificial exon 1 and a standard 5' splice site. Upon GAL4 induction, exon 1 can splice to the next exon downstream in the flanking genomic DNA, belonging to a random target gene. Exon 1 encodes a dual tag (FLAG epitope and mCherry red fluorescent protein), which becomes fused to the target protein. Hto was mobilized throughout the genome and then activated by eye-specific GAL4; an F1 screen for abnormal eye phenotypes was used to identify inserts that express disruptive fusion proteins. Approximately 1.7% of new inserts cause eye phenotypes. Of the first 23 verified target genes, 21 can be described as regulators of cell biology and development. Most are transcription factor genes, including AP-2, CG17181, cut, klu, mamo, Sox102F, and sv. Other target genes [l(1)G0232, nuf, pum, and Syt4] make cytoplasmic proteins, and these lines produce diverse fluorescence localization patterns. Hto permits the expression of stable carboxy-terminal subfragments of proteins, which are rarely tested in conventional genetic screens. Some of these may disrupt specific cell pathways, as exemplified by truncated forms of Mastermind and Nuf.

Late Replication Domains Are Evolutionary Conserved in the Drosophila Genome

Drosophila chromosomes are organized into distinct domains differing in their predominant chromatin composition, replication timing and evolutionary conservation. We show on a genome-wide level that genes whose order has remained unaltered across 9 Drosophila species display late replication timing and frequently map to the regions of repressive chromatin. This observation is consistent with the existence of extensive domains of repressive chromatin that replicate extremely late and have conserved gene order in the Drosophila genome. We suggest that such repressive chromatin domains correspond to a handful of regions that complete replication at the very end of S phase. We further demonstrate that the order of genes in these regions is rarely altered in evolution. Substantial proportion of such regions significantly coincide with large synteny blocks. This indicates that there are evolutionary mechanisms maintaining the integrity of these late-replicating chromatin domains. The synteny blocks corresponding to the extremely late-replicating regions in the D. melanogaster genome consistently display two-fold lower gene density across different Drosophila species.

Comparison of Ultra-Conserved Elements in Drosophilids and Vertebrates

Metazoan genomes contain many ultra-conserved elements (UCEs), long sequences identical between distant species. In this study we identified UCEs in drosophilid and vertebrate species with a similar level of phylogenetic divergence measured at protein-coding regions, and demonstrated that both the length and number of UCEs are larger in vertebrates. The proportion of non-exonic UCEs declines in distant drosophilids whilst an opposite trend was observed in vertebrates. We generated a set of 2,126 Sophophora UCEs by merging elements identified in several drosophila species and compared these to the eutherian UCEs identified in placental mammals. In contrast to vertebrates, the Sophophora UCEs are depleted around transcription start sites. Analysis of 52,954 P-element, piggyBac and Minos insertions in the D. melanogaster genome revealed depletion of the P-element and piggyBac insertions in and around the Sophophora UCEs. We examined eleven fly strains with transposon insertions into the intergenic UCEs and identified associated phenotypes in five strains. Four insertions behave as recessive lethals, and in one case we observed a suppression of the marker gene within the transgene, presumably by silenced chromatin around the integration site. To confirm the lethality is caused by integration of transposons we performed a phenotype rescue experiment for two stocks and demonstrated that the excision of the transposons from the intergenic UCEs restores viability. Sequencing of DNA after the transposon excision in one fly strain with the restored viability revealed a 47 bp insertion at the original transposon integration site suggesting that the nature of the mutation is important for the appearance of the phenotype. Our results suggest that the UCEs in flies and vertebrates have both common and distinct features, and demonstrate that a significant proportion of intergenic drosophila UCEs are sensitive to disruption.

OnTheFly: a database ofDrosophila melanogastertranscription factors and their binding sites

We present OnTheFly (http://bhapp.c2b2.columbia.edu/OnTheFly/index.php), a database comprising a systematic collection of transcription factors (TFs) of Drosophila melanogaster and their DNA-binding sites. TFs predicted in the Drosophila melanogaster genome are annotated and classified and their structures, obtained via experiment or homology models, are provided. All known preferred TF DNA-binding sites obtained from the B1H, DNase I and SELEX methodologies are presented. DNA shape parameters predicted for these sites are obtained from a high throughput server or from crystal structures of protein–DNA complexes where available. An important feature of the database is that all DNA-binding domains and their binding sites are fully annotated in a eukaryote using structural criteria and evolutionary homology. OnTheFly thus provides a comprehensive view of TFs and their binding sites that will be a valuable resource for deciphering non-coding regulatory DNA.

Transposable element 'roo' attaches to nuclear matrix of the Drosophila melanogaster.

The genome of eukaryotes is organized into structural units of chromatin loops. This higher order organization is supported by a nuclear skeleton called the nuclear matrix. The genomic DNA associated with the nuclear matrix is called the matrix associated region (MAR). Only a few genome-wide screens have been attempted, although many studies have characterized locusspecific MAR DNA sequences. In this study, a MAR DNA library was prepared from the Drosophila melanogaster Meigen (Diptera: Drosophilidae) genome. One of the sequences identified as a MAR was from a long terminal repeat region of ‘roo’ retrotransposon (roo MAR). Sequence analysis of roo MAR showed its distribution across the D. melanogaster genome. roo MAR also showed high sequence similarity with a previously identified MAR in Drosophila, namely the ‘gypsy’ retrotransposon. Analysis of the genes flanking roo MAR insertions in the Drosophila genome showed that genes were co-ordinately expressed. The results from the present study in D. melanogaster suggest this sequence plays an important role in genome organization and function. The findings point to an evolutionary role of retrotransposons in shaping the genomic architecture of eukaryotes.

Change point analysis of histone modifications reveals epigenetic blocks linking to physical domains

Histone modification is a vital epigenetic mechanism for transcriptional control in eukaryotes. High-throughput techniques have enabled whole-genome analysis of histone modifications in recent years. However, most studies assume one combination of histone modification invariantly translates to one transcriptional output regardless of local chromatin environment. In this study we hypothesize that, the genome is organized into local domains that manifest similar enrichment pattern of histone modification, which leads to orchestrated regulation of expression of genes with relevant bio- logical functions. We propose a multivariate Bayesian Change Point (BCP) model to segment the Drosophila melanogaster genome into consecutive blocks on the basis of combinatorial patterns of histone marks. By modeling the sparse distribution of histone marks across the chromosome with a zero-inflated Gaussian mixture, our partitions capture local BLOCKs that manifest relatively homogeneous enrichment pattern of histone modifications. We further characterized BLOCKs by their transcription levels, distribution of genes, degree of co-regulation and GO enrichment. Our results demonstrate that these BLOCKs, although inferred merely from histone modifications, reveal strong relevance with physical domains, which suggest their important roles in chromatin organization and coordinated gene regulation.

‘Escaping’ the X chromosome leads to increased gene expression in the male germline of Drosophila melanogaster

Genomic analyses of Drosophila species suggest that the X chromosome presents an unfavourable environment for the expression of genes in the male germline. A previous study in D. melanogaster used a reporter gene driven by a testis-specific promoter to show that expression was greatly reduced when the gene was inserted onto the X chromosome as compared with the autosomes. However, a limitation of this study was that only the expression regulated by a single, autosomal-derived promoter was investigated. To test for an increase in expression associated with 'escaping' the X chromosome, we analysed reporter gene expression driven by the promoters of three X-linked, testis-expressed genes (CG10920, CG12681 and CG1314) that were inserted randomly throughout the D. melanogaster genome. In all cases, insertions on the autosomes showed significantly higher expression than those on the X chromosome. Thus, even genes whose regulation has adapted to the X-chromosomal environment show increased male germline expression when relocated to an autosome. Our results provide direct experimental evidence for the suppression of X-linked gene expression in the Drosophila male germline that is independent of gene dose.

Characterisation of class I and II α-mannosidases from Drosophila melanogaster

Homology searches indicated that up to five class I α-mannosidases (glycohydrolase family 47) and eight class II α-mannosidases (glycohydrolase family 38) are encoded by the fruitfly (Drosophila melanogaster) genome. Selected example mannosidases were expressed in secreted form using the yeast Pichia pastoris. A number of characteristics of these enzymes were determined with p-nitrophenyl-α-mannoside as substrate; particularly striking were the low optima (pH 5) of three class II mannosidases most closely related to known lysosomal mannosidases and the distinct Co(II)-requirement of a mannosidase previously named ManIIb. Some of the recombinant mannosidases were demonstrably active towards oligomannosidic glycans, specifically, the Co(II)-requiring ManIIb, two 'acidic' mannosidases and the class I mas-1 mannosidase. Other than previous characterisations of the well-known Golgi mannosidase II, this is the first study summarising various properties of recombinant mannosidases from the fruitfly.

Optimized Position Weight Matrices in Prediction of Novel Putative Binding Sites for Transcription Factors in the Drosophila melanogaster Genome

Position weight matrices (PWMs) have become a tool of choice for the identification of transcription factor binding sites in DNA sequences. DNA-binding proteins often show degeneracy in their binding requirement and thus the overall binding specificity of many proteins is unknown and remains an active area of research. Although existing PWMs are more reliable predictors than consensus string matching, they generally result in a high number of false positive hits. Our previous study introduced a promising approach to PWM refinement in which known motifs are used to computationally mine putative binding sites directly from aligned promoter regions using composition of similar sites. In the present study, we extended this technique originally tested on single examples of transcription factors (TFs) and showed its capability to optimize PWM performance to predict new binding sites in the fruit fly genome. We propose refined PWMs in mono- and dinucleotide versions similarly computed for a large variety of transcription factors of Drosophila melanogaster. Along with the addition of many auxiliary sites the optimization includes variation of the PWM motif length, the binding sites location on the promoters and the PWM score threshold. To assess the predictive performance of the refined PWMs we compared them to conventional TRANSFAC and JASPAR sources. The results have been verified using performed tests and literature review. Overall, the refined PWMs containing putative sites derived from real promoter content processed using optimized parameters had better general accuracy than conventional PWMs.

Mitochondrial glutamate carriers from Drosophila melanogaster: Biochemical, evolutionary and modeling studies

The mitochondrial carriers are members of a family of transport proteins that mediate solute transport across the inner mitochondrial membrane. Two isoforms of the glutamate carriers, GC1 and GC2 (encoded by the SLC25A22 and SLC25A18 genes, respectively), have been identified in humans. Two independent mutations in SLC25A22 are associated with severe epileptic encephalopathy. In the present study we show that two genes (CG18347 and CG12201) phylogenetically related to the human GC encoding genes are present in the D. melanogaster genome. We have functionally characterized the proteins encoded by CG18347 and CG12201, designated as DmGC1p and DmGC2p respectively, by overexpression in Escherichia coli and reconstitution into liposomes. Their transport properties demonstrate that DmGC1p and DmGC2p both catalyze the transport of glutamate across the inner mitochondrial membrane. Computational approaches have been used in order to highlight residues of DmGC1p and DmGC2p involved in substrate binding. Furthermore, gene expression analysis during development and in various adult tissues reveals that CG18347 is ubiquitously expressed in all examined D. melanogaster tissues, while the expression of CG12201 is strongly testis-biased. Finally, we identified mitochondrial glutamate carrier orthologs in 49 eukaryotic species in order to attempt the reconstruction of the evolutionary history of the glutamate carrier function. Comparison of the exon/intron structure and other key features of the analyzed orthologs suggests that eukaryotic glutamate carrier genes descend from an intron-rich ancestral gene already present in the common ancestor of lineages that diverged as early as bilateria and radiata.

The use ofDrosophilaS2 cells in R&D and bioprocessing

Drosophila Schneider 2 (S2) cells have been available for approximately 40 years. Their use has intensified over the past 15 years: resolution of the whole Drosophila melanogaster genome and the amenability of S2 cells for siRNA-based studies are some of the reasons for their growing use. This review covers recent publications on use of S2 cells for research and manufacturing and points to some possible future developments in their use in the vaccine field. Relatively few groups have systematically developed the system to enable expression of challenging proteins. They demonstrated that these cells can constitute a robust, efficient protein expression system, with specific advantages such as homogeneous glycosylation profile; reproducibility between production runs; options for cultivation modes, including perfusion; and no cell lysis, leading to relatively low levels of contaminating host cell proteins. The platform has shown to be particularly well adapted for the production of challenging viral, malaria and immunotherapy antigens.

Studying Recombination with High-Throughput Sequencing: An Educational Primer for Use with “Fine-Scale Heterogeneity in Crossover Rate in the garnet-scalloped Region of the Drosophila melanogaster X chromosome”

An article by Singh and colleagues in this issue of GENETICS quantifies variation in recombination rate across a small region of the Drosophila melanogaster genome, providing an opportunity for instructors of genetics to introduce or reinforce important concepts such as recombination and recombination rate variation, genome sequencing, and sequence features of the genome. Additional background information, a detailed explanation of the methods used in this study, and discussion questions are provided.

Prediction of Drosophila melanogaster gene function using Support Vector Machines

BackgroundWhile the genomes of hundreds of organisms have been sequenced and good approaches exist for finding protein encoding genes, an important remaining challenge is predicting the functions of the large fraction of genes for which there is no annotation. Large gene expression datasets from microarray experiments already exist and many of these can be used to help assign potential functions to these genes. We have applied Support Vector Machines (SVM), a sigmoid fitting function and a stratified cross‐validation approach to analyze a large microarray experiment dataset from Drosophila melanogaster in order to predict possible functions for previously un‐annotated genes. A total of approximately 5043 different genes, or about one‐third of the predicted genes in the D. melanogaster genome, are represented in the dataset and 1854 (or 37%) of these genes are un‐annotated.Results39 Gene Ontology Biological Process (GO‐BP) categories were found with precision value equal or larger than 0.75, when recall was fixed at the 0.4 level. For two of those categories, we have provided additional support for assigning given genes to the category by showing that the majority of transcripts for the genes belonging in a given category have a similar localization pattern during embryogenesis. Additionally, by assessing the predictions using a confidence score, we have been able to provide a putative GO‐BP term for 1422 previously un‐annotated genes or about 77% of the un‐annotated genes represented on the microarray and about 19% of all of the un‐annotated genes in the D. melanogaster genome.ConclusionsOur study successfully employs a number of SVM classifiers, accompanied by detailed calibration and validation techniques, to generate a number of predictions for new annotations for D. melanogaster genes. The applied probabilistic analysis to SVM output improves the interpretability of the prediction results and the objectivity of the validation procedure.

An analysis and prediction of nucleosome positioning based on information content

Nucleosome positioning plays a key role in the regulation of many biological processes. In this study, the statistical difference of information content was investigated in nucleosome and linker DNA regions across eukaryotic organisms. By analyzing the information redundancy, D k , in Saccharomyces cerevisiae, Drosophila melanogaster, and Caenorhabditis elegans genomes, the short-range dominance of nucleotide correlation in nucleosome and linker DNA regions was confirmed. Significant difference of the D k value between the nucleosome and linker DNA regions was also found. The underlying reason for many successful oligonucleotide-based predictions of nucleosome positioning in eukaryotic model organisms may be attributed to the short-range dominance of nucleotide correlation in the nucleosome and linker DNA regions. When applying power spectrum analysis to the nucleosome and linker DNA regions, some obvious differences in sequence periodic signals were observed. The parameter F k was introduced to describe particular base correlation. Furthermore, the support vector machine combining F k was used to classify nucleosome and linker DNA regions in Homo sapiens, Oryzias latipes, C. elegans, Candida albicans, and S. cerevisiae. Independent test demonstrated that a good performance can be achieved by using this algorithm. This result further revealed that base correlation information has an important role in nucleosome positioning.

Strong Mutational Bias Toward Deletions in the Drosophila melanogaster Genome Is Compensated by Selection

Insertions and deletions (collectively indels) obviously have a major impact on genome evolution. However, before large-scale data on indel polymorphism became available, it was difficult to estimate the strength of selection acting on indel mutations. Here, we analyze indel polymorphism and divergence in different compartments of the Drosophila melanogaster genome: exons, introns of different lengths, and intergenic regions. Data on low-frequency polymorphisms indicate that 0.036–0.039 short (1–30 nt) insertion mutations and 0.085–0.092 short deletion mutations, with mean lengths 3.23 and 4.78, respectively, occur per single-nucleotide substitution. The excess of short deletion over short insertion mutations implies that indel mutations of these lengths should lead to a loss of approximately 0.30 nt per single-nucleotide replacement. However, polymorphism and divergence data show that this deletion bias is almost completely compensated by selection: Negative selection is stronger against deletions, whereas insertions are more likely to be favored by positive selection. Among the inframe low-frequency polymorphic mutations in exons, long introns, and intergenic regions, selection prevents a larger fraction of deletions (80–87%, depending on the type of the compartment) than of insertions (70–82%) or single-nucleotide substitutions (49–73%), from reaching high frequencies. The corresponding fractions were the lowest in short introns: 66%, 47%, and 15%, respectively, consistent with the weakest selective constraint in them. The McDonald–Kreitman test shows that 32–46% of the deletions and 60–73% of the insertions that were fixed in the recent evolution of D. melanogaster are adaptive, whereas this fraction is only 0–29% for single-nucleotide substitutions.

Anent the Genomics of Spermatogenesis in Drosophila melanogaster

An appreciable fraction of the Drosophila melanogaster genome is dedicated to male fertility. One approach to characterizing this subset of the genome is through the study of male-sterile mutations. We studied the relation between vital and male-fertility genes in three large autosomal regions that were saturated for lethal and male-sterile mutations. The majority of male-sterile mutations affect genes that are exclusively expressed in males. These genes are required only for male fertility, and several mutant alleles of each such gene were encountered. A few male-sterile mutations were alleles of vital genes that are expressed in both males and females. About one-fifth of the genes in Drosophila melanogaster show male-specific expression in adults. Although some earlier studies found a paucity of genes on the X chromosome showing male-biased expression, we did not find any significant differences between the X chromosome and the autosomes either in the relative frequencies of mutations to male sterility or in the frequencies of genes with male-specific expression in adults. Our results suggest that as much as 25% of the Drosophila genome may be dedicated to male fertility.

Turning gold into ‘junk’: transposable elements utilize central proteins of cellular networks

The numerous discovered cases of domesticated transposable element (TE) proteins led to the recognition that TEs are a significant source of evolutionary innovation. However, much less is known about the reverse process, whether and to what degree the evolution of TEs is influenced by the genome of their hosts. We addressed this issue by searching for cases of incorporation of host genes into the sequence of TEs and examined the systems-level properties of these genes using the Saccharomyces cerevisiae and Drosophila melanogaster genomes. We identified 51 cases where the evolutionary scenario was the incorporation of a host gene fragment into a TE consensus sequence, and we show that both the yeast and fly homologues of the incorporated protein sequences have central positions in the cellular networks. An analysis of selective pressure (Ka/Ks ratio) detected significant selection in 37% of the cases. Recent research on retrovirus-host interactions shows that virus proteins preferentially target hubs of the host interaction networks enabling them to take over the host cell using only a few proteins. We propose that TEs face a similar evolutionary pressure to evolve proteins with high interacting capacities and take some of the necessary protein domains directly from their hosts.

The Drosophila early ovarian transcriptome provides insight to the molecular causes of recombination rate variation across genomes

BackgroundEvidence in yeast indicates that gene expression is correlated with recombination activity and double-strand break (DSB) formation in some hotspots. Studies of nucleosome occupancy in yeast and mice also suggest that open chromatin influences the formation of DSBs. In Drosophila melanogaster, high-resolution recombination maps show an excess of DSBs within annotated transcripts relative to intergenic sequences. The impact of active transcription on recombination landscapes, however, remains unexplored in a multicellular organism. We then investigated the transcription profile during early meiosis in D. melanogaster females to obtain a glimpse at the relevant transcriptional dynamics during DSB formation, and test the specific hypothesis that DSBs preferentially target transcriptionally active genomic regions.ResultsOur study of transcript profiles of early- and late-meiosis using mRNA-seq revealed, 1) significant differences in gene expression, 2) new genes and exons, 3) parent-of-origin effects on transcription in early-meiosis stages, and 4) a nonrandom genomic distribution of transcribed genes. Importantly, genomic regions that are more actively transcribed during early meiosis show higher rates of recombination, and we ruled out DSB preference for genic regions that are not transcribed.ConclusionsOur results provide evidence in a multicellular organism that transcription during the initial phases of meiosis increases the likelihood of DSB and give insight into the molecular determinants of recombination rate variation across the D. melanogaster genome. We propose that a model where variation in gene expression plays a role altering the recombination landscape across the genome could provide a molecular, heritable and plastic mechanism to observed patterns of recombination variation, from the high level of intra-specific variation to the known influence of environmental factors and stress conditions.

75 OnTheFly database – structural basis to study TF’s DNA-binding specificity

We describe a systematic determination of Drosophila melanogaster transcription factor DNA-binding specificities. We annotated and classified all transcription factors (TFs) predicted in the Drosophila melanogaster genome (Pfreundt, U. et al. 2009) and collected the known preferred DNA binding sites of the TFs based on the B1H (Zhu, L.J. et al. 2011), DNaseI (Bergman C.M. et al. 2005), and SELEX (Slattery, M. et al. 2011) methods. Then, we identified the sequence and shape preferences for all DNA-binding proteins (Kuziemko A et al. 2011) and also characterized the shapes of their preferred DNA binding sites using structural models. The identification of the preferred DNA binding sites, and their shapes, for all DNA binding proteins will provide an unprecedented and extremely valuable database for anyone attempting to decipher noncoding regulatory DNA. Specifically, we showed that using structural criteria such as the width of the minor groove (Rohs R. et al. 2010), we distinguished DNA sequences bound by proteins which possess a Homeodomain from other proteins that possess a ZNF-C2H2 domain or ETS domains (pictured bellow from right to left). Furthermore, based on key TF–DNA interactions from the template structure stored in PDB, we superpose query DNA structure onto the template. Therefore, we obtain a homology model, where TF is from the template and DNA is from the query. This model will provide potential structural basis to study TF’s DNA-binding specificity. The insights from such a study could help in selecting the best DNA candidates to be bound by a certain TF for experimental testing. Finding such pairs will help to characterize the unique properties of protein–DNA interfaces and identify new drug target sites.