Systematic determination of patterns of gene expression during Drosophila embryogenesis

Enhancer-promoter communication at the Drosophila engrailed locus

Enhancers are often located many tens of kilobases away from the promoter they regulate, sometimes residing closer to the promoter of a neighboring gene. How do they know which gene to activate? We have used homing P[en] constructs to study the enhancer-promoter communication at the engrailed locus. Here we show that engrailed enhancers can act over large distances, even skipping over other transcription units, choosing the engrailed promoter over those of neighboring genes. This specificity is achieved in at least three ways. First, early acting engrailed stripe enhancers exhibit promoter specificity. Second, a proximal promoter-tethering element is required for the action of the imaginal disc enhancer(s). Our data suggest that there are two partially redundant promoter-tethering elements. Third, the long-distance action of engrailed enhancers requires a combination of the engrailed promoter and sequences within or closely linked to the promoter proximal Polycomb-group response elements. These data show that multiple mechanisms ensure proper enhancer-promoter communication at the Drosophila engrailed locus.

Essential requirement for RING finger E3 ubiquitin ligase Hakai in early embryonic development of Drosophila

Hakai is a RING finger type E3 ubiquitin ligase that is highly conserved in metazoans. Mammalian Hakai was shown to bind and ubiquitinate the intracellular domain of E-cadherin, and this activity is implicated in down-regulation of E-cadherin during v-Src-induced cellular transformation. To evaluate this model in vivo, we studied the function of the Drosophila homologue of Hakai. In cultured S2 cells, Drosophila Hakai and E-cadherin (Shotgun) formed a complex in a way distinct from the interaction described for mammalian counterparts. Hakai null mutants died during larval stages but this lethality could be offset by a HA-tagged Hakai construct. While zygotic Hakai function was dispensable for cell proliferation and differentiation in the wing disc epithelium, maternal Hakai mutants showed a variety of defects in epithelial integrity, including stochastic loss of E-cadherin expression and reduction of aPKC; defects in cell specification and cell migration were also observed. No increase of E-cadherin, however, was observed. Regulation of multiple target proteins under control of Hakai is, therefore, essential for early embryonic morphogenesis in Drosophila.

Unlocking the secrets of the genome

Despite the successes of genomics, little is known about how genetic information produces complex organisms. A look at the crucial functional elements of fly and worm genomes could change that.

The National Human Genome Research Institute's modENCODE project (the model organism ENCyclopedia Of DNA Elements) was set up in 2007 with the goal of identifying all the sequence-based functional elements in the genomes of two important experimental organisms, Caenorhabditis elegans and Drosophila melanogaster. Armed with modENCODE data, geneticists will be able to undertake the comprehensive molecular studies of regulatory networks that hold the key to how complex multicellular organisms arise from the list of instructions coded in the genome. In this issue, modENCODE team members outline their plan of campaign. Data from the project are to be made available on http://www.modencode.org and elsewhere as the work progresses.

Approaching Drosophila development through proteomic tools and databases: At the hub of the post-genomic era

The past decade has witnessed an explosion in the growth of proteomics. The completion of numerous genome sequences, the development of powerful protein analytical technologies, as well as the design of innovative bioinformatics tools have marked the beginning of a new post-genomic era. Proteomics, the large-scale analysis of proteins in an organism, organ or organelle encompasses different aspects: (1) the identification, analysis of post-translational modifications and quantification of proteins; (2) the study of protein-protein interactions; and (3) the functional analysis of interactome networks. Here, we briefly summarize the emerging analytical tools and databases that are paving the way for studying Drosophila development by proteomic approaches.

Stripy Ftz target genes are coordinately regulated by Ftz-F1

During development, cascades of regulatory genes act in a hierarchical fashion to subdivide the embryo into increasingly specified body regions. This has been best characterized in Drosophila, where genes encoding regulatory transcription factors form a network to direct the development of the basic segmented body plan. The pair-rule genes are pivotal in this process as they are responsible for the first subdivision of the embryo into repeated metameric units. The Drosophila pair-rule gene fushi tarazu (ftz) is a derived Hox gene expressed in and required for the development of alternate parasegments. Previous studies suggested that Ftz achieves its distinct regulatory specificity as a segmentation protein by interacting with a ubiquitously expressed cofactor, the nuclear receptor Ftz-F1. However, the downstream target genes regulated by Ftz and other pair-rule genes to direct segment formation are not known. In this study, we selected candidate Ftz targets by virtue of their early expression in Ftz-like stripes. This identified two new Ftz target genes, drumstick (drm) and no ocelli (noc), and confirmed that Ftz regulates a serotonin receptor (5-HT2). These are the earliest Ftz targets identified to date and all are coordinately regulated by Ftz-F1. Engrailed (En), the best-characterized Ftz/Ftz-F1 downstream target, is not an intermediate in regulation. The drm genomic region harbors two separate seven-stripe enhancers, identified by virtue of predicted Ftz-F1 binding sites, and these sites are necessary for stripe expression in vivo. We propose that pair-rule genes, exemplified by Ftz/Ftz-F1, promote segmentation by acting at different hierarchical levels, regulating first, other segmentation genes; second, other regulatory genes that in turn control specific cellular processes such as tissue differentiation; and, third, 'segmentation realizator genes' that are directly involved in morphogenesis.

Relationship between gene co-expression and sharing of transcription factor binding sites in Drosophila melanogaster

MOTIVATION

In functional genomics, it is frequently useful to correlate expression levels of genes to identify transcription factor binding sites (TFBS) via the presence of common sequence motifs. The underlying assumption is that co-expressed genes are more likely to contain shared TFBS and, thus, TFBS can be identified computationally. Indeed, gene pairs with a very high expression correlation show a significant excess of shared binding sites in yeast. We have tested this assumption in a more complex organism, Drosophila melanogaster, by using experimentally determined TFBS and microarray expression data. We have also examined the reverse relationship between the expression correlation and the extent of TFBS sharing.

RESULTS

Pairs of genes with shared TFBS show, on average, a higher degree of co-expression than those with no common TFBS in Drosophila. However, the reverse does not hold true: gene pairs with high expression correlations do not share significantly larger numbers of TFBS. Exception to this observation exists when comparing expression of genes from the earliest stages of embryonic development. Interestingly, semantic similarity between gene annotations (Biological Process) is much better associated with TFBS sharing, as compared to the expression correlation. We discuss these results in light of reverse engineering approaches to computationally predict regulatory sequences by using comparative genomics.

The genetic basis of transgressive ovary size in honeybee workers

Ovarioles are the functional unit of the female insect reproductive organs and the number of ovarioles per ovary strongly influences egg-laying rate and fecundity. Social evolution in the honeybee (Apis mellifera) has resulted in queens with 200-360 total ovarioles and workers with usually 20 or less. In addition, variation in ovariole number among workers relates to worker sensory tuning, foraging behavior, and the ability to lay unfertilized male-destined eggs. To study the genetic architecture of worker ovariole number, we performed a series of crosses between Africanized and European bees that differ in worker ovariole number. Unexpectedly, these crosses produced transgressive worker phenotypes with extreme ovariole numbers that were sensitive to the social environment. We used a new selective pooled DNA interval mapping approach with two Africanized backcrosses to identify quantitative trait loci (QTL) underlying the transgressive ovary phenotype. We identified one QTL on chromosome 11 and found some evidence for another QTL on chromosome 2. Both QTL regions contain plausible functional candidate genes. The ovariole number of foragers was correlated with the sugar concentration of collected nectar, supporting previous studies showing a link between worker physiology and foraging behavior. We discuss how the phenotype of extreme worker ovariole numbers and the underlying genetic factors we identified could be linked to the development of queen traits.

Motif composition, conservation and condition-specificity of single and alternative transcription start sites in the Drosophila genome

A map of transcription start sites across the Drosophila genome, providing insights into initiation patterns and spatiotemporal conditions.

Background

Transcription initiation is a key component in the regulation of gene expression. mRNA 5' full-length sequencing techniques have enhanced our understanding of mammalian transcription start sites (TSSs), revealing different initiation patterns on a genomic scale.

Results

To identify TSSs in Drosophila melanogaster, we applied a hierarchical clustering strategy on available 5' expressed sequence tags (ESTs) and identified a high quality set of 5,665 TSSs for approximately 4,000 genes. We distinguished two initiation patterns: 'peaked' TSSs, and 'broad' TSS cluster groups. Peaked promoters were found to contain location-specific sequence elements; conversely, broad promoters were associated with non-location-specific elements. In alignments across other Drosophila genomes, conservation levels of sequence elements exceeded 90% within the melanogaster subgroup, but dropped considerably for distal species. Elements in broad promoters had lower levels of conservation than those in peaked promoters. When characterizing the distributions of ESTs, 64% of TSSs showed distinct associations to one out of eight different spatiotemporal conditions. Available whole-genome tiling array time series data revealed different temporal patterns of embryonic activity across the majority of genes with distinct alternative promoters. Many genes with maternally inherited transcripts were found to have alternative promoters utilized later in development. Core promoters of maternally inherited transcripts showed differences in motif composition compared to zygotically active promoters.

Conclusions

Our study provides a comprehensive map of Drosophila TSSs and the conditions under which they are utilized. Distinct differences in motif associations with initiation pattern and spatiotemporal utilization illustrate the complex regulatory code of transcription initiation.

Expression of the apoptosis inducer gene head involution defective in primordial germ cells of the Drosophila embryo requires eiger, p53, and loki function

Nanos (Nos) is an evolutionarily conserved protein essential for the maintenance of primordial germ cells (PGCs). In Drosophila, the PGCs or pole cells express head involution defective (hid), which is required for caspase activation, but its translation is repressed by maternal Nos. In the absence of Nos activity, translation of hid mRNA into protein induces apoptosis in pole cells. However, it remains unclear how hid mRNA is regulated in pole cells. Here, we report that hid expression requires eiger (egr), a tumor necrosis factor ligand (TNF) homologue, which is induced in pole cells by decapentaplegic (dpp). In addition, we demonstrate that p53 and loki (lok), a damage-activated kinase known to be required for p53 phosphorylation, are both required for hid expression in pole cells. Since maternal lok mRNA is enriched in pole cells, it is possible that ubiquitously distributed p53 is activated in pole cells by maternal Lok. We propose that hid expression is activated in a pole cell-specific manner by loki/p53 and dpp/egr during embryogenesis.

Determination of gene expression patterns using high-throughput RNA in situ hybridization to whole-mount Drosophila embryos

We describe a high-throughput protocol for RNA in situ hybridization (ISH) to Drosophila embryos in a 96-well format. cDNA or genomic DNA templates are amplified by PCR and then digoxigenin-labeled ribonucleotides are incorporated into antisense RNA probes by in vitro transcription. The quality of each probe is evaluated before ISH using a RNA probe quantification (dot blot) assay. RNA probes are hybridized to fixed, mixed-staged Drosophila embryos in 96-well plates. The resulting stained embryos can be examined and photographed immediately or stored at 4 degrees C for later analysis. Starting with fixed, staged embryos, the protocol takes 6 d from probe template production through hybridization. Preparation of fixed embryos requires a minimum of 2 weeks to collect embryos representing all stages. The method has been used to determine the expression patterns of over 6,000 genes throughout embryogenesis.

Characterization of quantitative trait loci for the age of first foraging in honey bee workers

Identifying the basis of quantitative trait loci (QTL) remains challenging for the study of complex traits, such as behavior. The honey bee is a good model combining interesting social behavior with a high recombination rate that facilitates this identification. Several studies have focused on the pollen hoarding syndrome, identifying multiple QTL as the genetic basis of its behavioral components. One component, the age of first foraging, is central for colony organization and four QTL were previously described without identification of their genomic location. Enabled by the honey bee genome project, this study provides data from multiple experiments to scrutinize these QTL, including individual and pooled SNP mapping, sequencing of AFLP markers, and microsatellite genotyping. The combined evidence confirms and localizes two of the previous QTL on chromosome four and five, dismisses the other two, and suggests one novel genomic region on chromosome eleven to influence the age of first foraging. Among the positional candidates the Ank2, PKC, Erk7, and amontillado genes stand out due to corroborating functional evidence. This study thus demonstrates the power of combined, genome-based approaches to enable targeted studies of a manageable set of candidate genes for natural behavioral variation in the important, complex social trait "age of first foraging".

Regulation of glutamate receptor subunit availability by microRNAs

The efficacy of synaptic transmission depends, to a large extent, on postsynaptic receptor abundance. The molecular mechanisms controlling receptor abundance are poorly understood. We tested whether abundance of postsynaptic glutamate receptors (GluRs) in Drosophila neuromuscular junctions is controlled by microRNAs, and provide evidence that it is. We show here that postsynaptic knockdown of dicer-1, the endoribonuclease necessary for microRNA synthesis, leads to large increases in postsynaptic GluR subunit messenger RNA and protein. Specifically, we measured increases in GluRIIA and GluRIIB but not GluRIIC. Further, knockout of MiR-284, a microRNA predicted to bind to GluRIIA and GluRIIB but not GluRIIC, increases expression of GluRIIA and GluRIIB but not GluRIIC proportional to the number of predicted binding sites in each transcript. Most of the de-repressed GluR protein, however, does not appear to be incorporated into functional receptors, and only minor changes in synaptic strength are observed, which suggests that microRNAs primarily regulate Drosophila receptor subunit composition rather than overall receptor abundance or synaptic strength.

In situ analysis of gene expression in plants

In the post-genomic era, it is necessary to adapt methods for gene expression and functional analyses to more high-throughput levels of processing. mRNA in situ hybridization (ISH) remains a powerful tool for obtaining information regarding a gene's temporal and spatial expression pattern and can therefore be used as a starting point to define the function of a gene or a whole set of genes. We have deconstructed 'traditional' ISH techniques described for a range of organisms and developed protocols for ISH that adapt and integrate a degree of automation to standardized and shortened protocols. We have adapted this technique as a high-throughput means of gene expression analysis on wax-embedded plant tissues and also on whole-mount tissues. We have used wax-embedded wheat grains and Arabidopsis floral meristems and whole-mount Arabidopsis roots as test systems and show that it is capable of highly parallel processing.

Discovering structural cis-regulatory elements by modeling the behaviors of mRNAs

Gene expression is regulated at each step from chromatin remodeling through translation and degradation. Several known RNA-binding regulatory proteins interact with specific RNA secondary structures in addition to specific nucleotides. To provide a more comprehensive understanding of the regulation of gene expression, we developed an integrative computational approach that leverages functional genomics data and nucleotide sequences to discover RNA secondary structure-defined cis-regulatory elements (SCREs). We applied our structural cis-regulatory element detector (StructRED) to microarray and mRNA sequence data from Saccharomyces cerevisiae, Drosophila melanogaster, and Homo sapiens. We recovered the known specificities of Vts1p in yeast and Smaug in flies. In addition, we discovered six putative SCREs in flies and three in humans. We characterized the SCREs based on their condition-specific regulatory influences, the annotation of the transcripts that contain them, and their locations within transcripts. Overall, we show that modeling functional genomics data in terms of combined RNA structure and sequence motifs is an effective method for discovering the specificities and regulatory roles of RNA-binding proteins.

A bag-of-words approach for Drosophila gene expression pattern annotation

Background

Drosophila gene expression pattern images document the spatiotemporal dynamics of gene expression during embryogenesis. A comparative analysis of these images could provide a fundamentally important way for studying the regulatory networks governing development. To facilitate pattern comparison and searching, groups of images in the Berkeley Drosophila Genome Project (BDGP) high-throughput study were annotated with a variable number of anatomical terms manually using a controlled vocabulary. Considering that the number of available images is rapidly increasing, it is imperative to design computational methods to automate this task.

Results

We present a computational method to annotate gene expression pattern images automatically. The proposed method uses the bag-of-words scheme to utilize the existing information on pattern annotation and annotates images using a model that exploits correlations among terms. The proposed method can annotate images individually or in groups (e.g., according to the developmental stage). In addition, the proposed method can integrate information from different two-dimensional views of embryos. Results on embryonic patterns from BDGP data demonstrate that our method significantly outperforms other methods.

Conclusion

The proposed bag-of-words scheme is effective in representing a set of annotations assigned to a group of images, and the model employed to annotate images successfully captures the correlations among different controlled vocabulary terms. The integration of existing annotation information from multiple embryonic views improves annotation performance.

Dietary zinc absorption is mediated by ZnT1 in Drosophila melanogaster

Zinc is an essential nutritional factor involved in many key biological processes. However, the physiological function of zinc transporters at the organismal level is not well characterized. Early embryonic lethality of Znt1 knockout mice precludes functional analysis of the role of ZnT1 in dietary zinc absorption. Here, we report the identification and characterization of the Drosophila ZnT1 orthologue, dZnT1, for its role in Drosophila dietary zinc absorption. In cell culture, dZnT1 promoted zinc transport to reduce cytoplasmic zinc levels. Ubiquitous RNA interference of dZnT1 in Drosophila resulted in developmental arrest under restriction of dietary zinc, while dZnT1-overexpressing flies exhibited hypersensitivity to zinc. dZnT1 was prominently expressed in restricted regions of the midgut and exhibited a distribution on the basolateral membrane of the enterocytes. Gut-specific silencing of dZnT1 was sufficient to evoke lethality under zinc scarcity. Human ZnT1, but not ZnT7 or ZnT4, could rescue the zinc-acquiring defects caused by dZnT1 silencing. Taken together, our results proved that dZnT1 is a key zinc transporter in dietary zinc absorption, functioning by pumping zinc out of the enterocytes across the basolateral membrane. This study will be helpful in understanding the fundamental process of acquiring dietary zinc in higher eukaryotes.

Global implications of mRNA localization pathways in cellular organization

Genome expression profiling has led to the important realization that RNA molecules are more numerous and diverse than previously expected. One aspect of RNA biology that is just beginning to be fully appreciated is the extent to which mRNAs are transported to specific subcellular destinations before being translated, an exquisite mechanism for targeting proteins where they are required in the cell. While several excellent reviews have discussed the subject of mRNA localization, it is only in recent years that high-throughput technologies have been applied to address issues such as the extent and diversity of RNA localization events and mechanisms. This review focuses on these recent functional genomic approaches, their implications, and the new tools and methods that will be needed to further elucidate mRNA localization pathways.

Anterior-posterior positional information in the absence of a strong Bicoid gradient

The Bicoid (Bcd) transcription factor is distributed as a long-range concentration gradient along the anterior posterior (AP) axis of the Drosophila embryo. Bcd is required for the activation of a series of target genes, which are expressed at specific positions within the gradient. Here we directly tested whether different concentration thresholds within the Bcd gradient establish the relative positions of its target genes by flattening the gradient and systematically varying expression levels. Genome-wide expression profiles were used to estimate the total number of Bcd target genes, and a general correlation was found between the Bcd concentration required for activation and the positions where target genes are expressed in wild-type embryos. However, concentrations required for target gene activation in embryos with flattened Bcd were consistently lower than those present at each target gene's position in the wild-type gradient, suggesting that Bcd is in excess at every position along the AP axis. Also, several Bcd target genes were positioned in correctly ordered stripes in embryos with flattened Bcd, and we suggest that these stripes are normally regulated by interactions between Bcd and the terminal patterning system. Our findings argue strongly against the strict interpretation of the Bcd morphogen hypothesis, and support the idea that target gene positioning involves combinatorial interactions that are mediated by the binding site architecture of each gene's cis-regulatory elements.

Atlas of gene expression in the developing kidney at microanatomic resolution

Kidney development is based on differential cell-type-specific expression of a vast number of genes. While multiple critical genes and pathways have been elucidated, a genome-wide analysis of gene expression within individual cellular and anatomic structures is lacking. Accomplishing this could provide significant new insights into fundamental developmental mechanisms such as mesenchymal-epithelial transition, inductive signaling, branching morphogenesis, and segmentation. We describe here a comprehensive gene expression atlas of the developing mouse kidney based on the isolation of each major compartment by either laser capture microdissection or fluorescence-activated cell sorting, followed by microarray profiling. The resulting data agree with known expression patterns and additional in situ hybridizations. This kidney atlas allows a comprehensive analysis of the progression of gene expression states during nephrogenesis, as well as discovery of potential growth factor-receptor interactions. In addition, the results provide deeper insight into the genetic regulatory mechanisms of kidney development.

Identification of receptor-tyrosine-kinase-signaling target genes reveals receptor-specific activities and pathway branchpoints during Drosophila development

Receptor tyrosine kinases (RTKs) are an important family of signaling molecules with the unusual property that they are able to transduce their signals using the same downstream pathways. This has led to an unresolved debate as to whether individual receptors are interchangeable, or if each receptor can mediate specific downstream responses. To address this question, we have conducted a screen to identify target genes whose expression is differentially modulated by RTKs and their downstream pathway components. Using whole-mount in situ hybridization in Drosophila embryos exposed to constitutively active RTK pathway signaling, along with quantitative RT-PCR, we found that a significant fraction of target genes respond differentially in a spatial and/or quantitative manner. This includes differential responses to EGF receptor vs. fibroblast growth factor receptor signaling as well as to more downstream components such as Ras1 and pointed. We show that not only genes but also individual alternative transcripts can respond differently to signaling, and we present evidence that the differential responses can be mediated at the transcriptional level. Our results demonstrate that different RTKs can elicit distinct transcriptional responses, and the target genes obtained from our screen provide a valuable resource for further exploration of the mechanisms underlying this signaling specificity.

Image processing and analysis for quantifying gene expression from early Drosophila embryos

Correlation of quantities of transcriptional activators and repressors with the mRNA output of target genes is a central issue for modeling gene regulation. In multicellular organisms, both spatial and temporal differences in gene expression must be taken into account; this can be achieved by use of in situ hybridization followed by confocal laser scanning microscopy (CLSM). Here we present a method to correlate the protein levels of the short-range repressor Giant with lacZ mRNA produced by reporter genes using images of Drosophila blastoderm embryos taken by CLSM. The image stacks from CLSM are processed using a semiautomatic algorithm to produce correlations between the repressor levels and lacZ mRNA reporter genes. We show that signals derived from CLSM are proportional to actual mRNA levels. Our analysis reveals that a suggested parabolic form of the background fluorescence in confocal images of early Drosophila embryos is evident most prominently in flattened specimens, with intact embryos exhibiting a more linear background. The data extraction described in this paper is primarily conceived for analysis of synthetic reporter genes that are designed to decipher cis-regulatory grammar, but the techniques are generalizable for quantitative analysis of other engineered or endogenous genes in embryos.

Drosophila Gene Expression Pattern Annotation through Multi-Instance Multi-Label Learning

The Berkeley Drosophila Genome Project (BDGP) has produced a large number of gene expression patterns, many of which have been annotated textually with anatomical and developmental terms. These terms spatially correspond to local regions of the images; however, they are attached collectively to groups of images, such that it is unknown which term is assigned to which region of which image in the group. This poses a challenge to the development of the computational method to automate the textual description of expression patterns contained in each image. In this paper, we show that the underlying nature of this task matches well with a new machine learning framework, Multi-Instance Multi-Label learning (MIML). We propose a new MIML support vector machine to solve the problems that beset the annotation task. Empirical study shows that the proposed method outperforms the state-of-the-art Drosophila gene expression pattern annotation methods.

FlyEx, the quantitative atlas on segmentation gene expression at cellular resolution

The datasets on gene expression are the valuable source of information about the functional state of an organism. Recently, we have acquired the large dataset on expression of segmentation genes in the Drosophila blastoderm. To provide efficient access to the data, we have developed the FlyEx database (http://urchin.spbcas.ru/flyex). FlyEx contains 4716 images of 14 segmentation gene expression patterns obtained from 1579 embryos and 9,500,000 quantitative data records. Reference data are available for all segmentation genes in cycles 11-13 and all temporal classes of cycle 14A. FlyEx supports operations on images of gene expression patterns. The database can be used to examine the quality of data, analyze the dynamics of formation of segmentation gene expression domains, as well as to estimate the variability of gene expression patterns. Currently, a user is able to monitor and analyze the dynamics of formation of segmentation gene expression domains over the whole period of segment determination, that amounts to 1.5 h of development. FlyEx supports the data downloads and construction of personal reference datasets, that makes it possible to more effectively use and analyze data.

Integrating computational biology and forward genetics in Drosophila

Genetic screens are powerful methods for the discovery of gene-phenotype associations. However, a systems biology approach to genetics must leverage the massive amount of "omics" data to enhance the power and speed of functional gene discovery in vivo. Thus far, few computational methods for gene function prediction have been rigorously tested for their performance on a genome-wide scale in vivo. In this work, we demonstrate that integrating genome-wide computational gene prioritization with large-scale genetic screening is a powerful tool for functional gene discovery. To discover genes involved in neural development in Drosophila, we extend our strategy for the prioritization of human candidate disease genes to functional prioritization in Drosophila. We then integrate this prioritization strategy with a large-scale genetic screen for interactors of the proneural transcription factor Atonal using genomic deficiencies and mutant and RNAi collections. Using the prioritized genes validated in our genetic screen, we describe a novel genetic interaction network for Atonal. Lastly, we prioritize the whole Drosophila genome and identify candidate gene associations for ten receptor-signaling pathways. This novel database of prioritized pathway candidates, as well as a web application for functional prioritization in Drosophila, called Endeavour-HighFly, and the Atonal network, are publicly available resources. A systems genetics approach that combines the power of computational predictions with in vivo genetic screens strongly enhances the process of gene function and gene-gene association discovery.

Drosophila Gene Expression Pattern Annotation Using Sparse Features and Term-Term Interactions

The Drosophila gene expression pattern images document the spatial and temporal dynamics of gene expression and they are valuable tools for explicating the gene functions, interaction, and networks during Drosophila embryogenesis. To provide text-based pattern searching, the images in the Berkeley Drosophila Genome Project (BDGP) study are annotated with ontology terms manually by human curators. We present a systematic approach for automating this task, because the number of images needing text descriptions is now rapidly increasing. We consider both improved feature representation and novel learning formulation to boost the annotation performance. For feature representation, we adapt the bag-of-words scheme commonly used in visual recognition problems so that the image group information in the BDGP study is retained. Moreover, images from multiple views can be integrated naturally in this representation. To reduce the quantization error caused by the bag-of-words representation, we propose an improved feature representation scheme based on the sparse learning technique. In the design of learning formulation, we propose a local regularization framework that can incorporate the correlations among terms explicitly. We further show that the resulting optimization problem admits an analytical solution. Experimental results show that the representation based on sparse learning outperforms the bag-of-words representation significantly. Results also show that incorporation of the term-term correlations improves the annotation performance consistently.

Selective maintenance of Drosophila tandemly arranged duplicated genes during evolution

Genes occurring in conserved, tandemly-arrayed clusters in Drosophila melanogaster are co-expressed to a much higher extent than other duplicated genes.

Background

The physical organization and chromosomal localization of genes within genomes is known to play an important role in their function. Most genes arise by duplication and move along the genome by random shuffling of DNA segments. Higher order structuring of the genome occurs in eukaryotes, where groups of physically linked genes are co-expressed. However, the contribution of gene duplication to gene order has not been analyzed in detail, as it is believed that co-expression due to recent duplicates would obscure other domains of co-expression.

Results

We have catalogued ordered duplicated genes in Drosophila melanogaster, and found that one in five of all genes is organized as tandem arrays. Furthermore, among arrays that have been spatially conserved over longer periods than would be expected on the basis of random shuffling, a disproportionate number contain genes encoding developmental regulators. Using in situ gene expression data for more than half of the Drosophila genome, we find that genes in these conserved clusters are co-expressed to a much higher extent than other duplicated genes.

Conclusions

These results reveal the existence of functional constraints in insects that retain copies of genes encoding developmental and regulatory proteins as neighbors, allowing their co-expression. This co-expression may be the result of shared cis-regulatory elements or a shared need for a specific chromatin structure. Our results highlight the association between genome architecture and the gene regulatory networks involved in the construction of the body plan.

PROCESS FLOW FOR CLASSIFICATION AND CLUSTERING OF FRUIT FLY GENE EXPRESSION PATTERNS

The rapidly growing collection of fruit fly embryo images makes automated Image Segmentation and classification an indispensable requirement for a large-scale analysis of in situ hybridization (ISH) - gene expression patterns (GEP). We present here such an automated process flow for Segmenting, Classification, and Clustering large-scale sets of Drosophila melanogaster GEP that is capable of dealing with most of the complications implicated in the images.

A combinatorial code for pattern formation in Drosophila oogenesis

Two-dimensional patterning of the follicular epithelium in Drosophila oogenesis is required for the formation of three-dimensional eggshell structures. Our analysis of a large number of published gene expression patterns in the follicle cells suggests that they follow a simple combinatorial code based on six spatial building blocks and the operations of union, difference, intersection, and addition. The building blocks are related to the distribution of inductive signals, provided by the highly conserved epidermal growth factor receptor and bone morphogenetic protein signaling pathways. We demonstrate the validity of the code by testing it against a set of patterns obtained in a large-scale transcriptional profiling experiment. Using the proposed code, we distinguish 36 distinct patterns for 81 genes expressed in the follicular epithelium and characterize their joint dynamics over four stages of oogenesis. The proposed combinatorial framework allows systematic analysis of the diversity and dynamics of two-dimensional transcriptional patterns and guides future studies of gene regulation.

Cardiac expression of the Drosophila Transglutaminase (CG7356) gene is directly controlled by myocyte enhancer factor-2

The myocyte enhancer factor-2 (MEF2) family of transcription factors plays key roles in the activation of muscle structural genes. In Drosophila, MEF2 accumulates at high levels in the embryonic muscles, where it activates target genes throughout the mesoderm. Here, we identify the Transglutaminase gene (Tg; CG7356) as a direct transcriptional target of MEF2 in the cardiac musculature. Tg is expressed in cells forming the inflow tracts of the dorsal vessel, and we identify the enhancer responsible for this expression. The enhancer contains three binding sites for MEF2, and can be activated by MEF2 in tissue culture and in vivo. Moreover, loss of MEF2 function, or removal of the MEF2 binding sites from the enhancer, results in loss of Tg expression. These studies identify a new MEF2 target in the cardiac musculature. These studies provide a possible mechanism for the activation of transglutaminase genes.

DroID: the Drosophila Interactions Database, a comprehensive resource for annotated gene and protein interactions

BACKGROUND

Charting the interactions among genes and among their protein products is essential for understanding biological systems. A flood of interaction data is emerging from high throughput technologies, computational approaches, and literature mining methods. Quick and efficient access to this data has become a critical issue for biologists. Several excellent multi-organism databases for gene and protein interactions are available, yet most of these have understandable difficulty maintaining comprehensive information for any one organism. No single database, for example, includes all available interactions, integrated gene expression data, and comprehensive and searchable gene information for the important model organism, Drosophila melanogaster.

DESCRIPTION

DroID, the Drosophila Interactions Database, is a comprehensive interactions database designed specifically for Drosophila. DroID houses published physical protein interactions, genetic interactions, and computationally predicted interactions, including interologs based on data for other model organisms and humans. All interactions are annotated with original experimental data and source information. DroID can be searched and filtered based on interaction information or a comprehensive set of gene attributes from Flybase. DroID also contains gene expression and expression correlation data that can be searched and used to filter datasets, for example, to focus a study on sub-networks of co-expressed genes. To address the inherent noise in interaction data, DroID employs an updatable confidence scoring system that assigns a score to each physical interaction based on the likelihood that it represents a biologically significant link.

CONCLUSION

DroID is the most comprehensive interactions database available for Drosophila. To facilitate downstream analyses, interactions are annotated with original experimental information, gene expression data, and confidence scores. All data in DroID are freely available and can be searched, explored, and downloaded through three different interfaces, including a text based web site, a Java applet with dynamic graphing capabilities (IM Browser), and a Cytoscape plug-in. DroID is available at http://www.droidb.org.

The invertebrate B(0) system transporter, D. melanogaster NAT1, has unique d-amino acid affinity and mediates gut and brain functions

The CG3252 gene product, DmNAT1, represents the first Nutrient Amino acid Transporter cloned from Drosophila. It absorbs a broader set of neutral amino acids versus earlier characterized insect NATs and mammalian NATs-B(0) system transporters from the Sodium Neurotransmitter symporter Family (SNF, a.k.a. solute carrier family 6, SLC6). In addition to B(0)-specific l-substrates, DmNAT1 equally or more effectively transports d-amino acids with sub-millimolar affinities and 1:1 sodium:amino acid transport stoichiometry. DmNAT1 is strongly transcribed in the absorptive and secretory regions of the larval alimentary canal and larval brain, revealing its roles in the primary absorption and redistribution of large neutral l-amino acids as well as corresponding d-isomers. The absorption of d-amino acids via DmNAT1 may benefit the acquisition of fermented and symbiotic products, and may support the unique capacity of fruit fly larvae to utilize a diet with substitution of essential amino acids by d-isomers. It also suggests a remarkable adaptive plasticity of NAT-SLC6 mechanisms via alterations of a few identifiable sites in the substrate-binding pocket. The strong transcription in the brain suggests roles for DmNAT1 in neuronal nutrition and clearance of l-neutral amino acids from the fly brain. In addition, neuronal DmNAT1 may absorb synaptic d-serine and modulate NMDA receptor-coupled signal transduction. The characterization of the first invertebrate B(0)-like transporter extends the biological roles of the SLC6 family, revealing adaptations for the absorption of d-isomers of the essential amino acids. These findings suggest that some members of the NAT-SLC6 subfamily are evolving specific properties which contribute to nutrient symbiotic relationships and neuronal functions.

Cardiac gene regulatory networks in Drosophila

The Drosophila system has proven a powerful tool to help unlock the regulatory processes that occur during specification and differentiation of the embryonic heart. In this review, we focus upon a temporal analysis of the molecular events that result in heart formation in Drosophila, with a particular emphasis upon how genomic and other cutting-edge approaches are being brought to bear upon the subject. We anticipate that systems-level approaches will contribute greatly to our comprehension of heart development and disease in the animal kingdom.

Mapping the gene expression universe

Methods to globally survey gene expression provide valuable insights into gene function during development. In particular, comprehensive in situ hybridization studies have demonstrated that gene expression patterns are extraordinarily diverse and new imaging techniques have been introduced to capture these patterns with higher resolution at the tissue, cellular, and subcellular levels. The analysis of massive image databases can be greatly facilitated by computer vision techniques once annotated image sets reach the crucial mass sufficient to train the computer in pattern recognition. Ultimately, genome-wide atlases of gene expression during development will record gene activity in living animals with at least cellular resolution and in the context of morphogenetic events. These emerging datasets will lead to great advances in the field of comparative genomics and revolutionize our ability to decipher and model developmental processes for a variety of organisms.

The loci of evolution: how predictable is genetic evolution?

Is genetic evolution predictable? Evolutionary developmental biologists have argued that, at least for morphological traits, the answer is a resounding yes. Most mutations causing morphological variation are expected to reside in the cis-regulatory, rather than the coding, regions of developmental genes. This "cis-regulatory hypothesis" has recently come under attack. In this review, we first describe and critique the arguments that have been proposed in support of the cis-regulatory hypothesis. We then test the empirical support for the cis-regulatory hypothesis with a comprehensive survey of mutations responsible for phenotypic evolution in multicellular organisms. Cis-regulatory mutations currently represent approximately 22% of 331 identified genetic changes although the number of cis-regulatory changes published annually is rapidly increasing. Above the species level, cis-regulatory mutations altering morphology are more common than coding changes. Also, above the species level cis-regulatory mutations predominate for genes not involved in terminal differentiation. These patterns imply that the simple question "Do coding or cis-regulatory mutations cause more phenotypic evolution?" hides more interesting phenomena. Evolution in different kinds of populations and over different durations may result in selection of different kinds of mutations. Predicting the genetic basis of evolution requires a comprehensive synthesis of molecular developmental biology and population genetics.

The diversity of O-linked glycans expressed during Drosophila melanogaster development reflects stage- and tissue-specific requirements for cell signaling

Appropriate glycoprotein O-glycosylation is essential for normal development and tissue function in multicellular organisms. To comprehensively assess the developmental and functional impact of altered O-glycosylation, we have extensively analyzed the non-glycosaminoglycan, O-linked glycans expressed in Drosophila embryos. Through multidimensional mass spectrometric analysis of glycans released from glycoproteins by beta-elimination, we detected novel as well as previously reported O-glycans that exhibit developmentally modulated expression. The core 1 mucin-type disaccharide (Galbeta1-3GalNAc) is the predominant glycan in the total profile. HexNAcitol, hexitol, xylosylated hexitol, and branching extension of core 1 with HexNAc (to generate core 2 glycans) were also evident following release and reduction. After Galbeta1-3GalNAc, the next most prevalent glycans were a mixture of novel, isobaric, linear, and branched forms of a glucuronyl core 1 disaccharide. Other less prevalent structures were also extended with HexA, including an O-fucose glycan. Although the expected disaccharide product of the Fringe glycosyltransferase, (GlcNAcbeta1-3)fucitol, was not detectable in whole embryos, mass spectrometry fragmentation and exoglycosidase sensitivity defined a novel glucuronyl trisaccharide as GlcNAcbeta1-3(GlcAbeta1-4)fucitol. Consistent with the spatial distribution of the Fringe function, the GlcA-extended form of the Fringe product was enriched in the dorsal portion of the wing imaginal disc. Furthermore, loss of Fringe activity reduced the prevalence of the O-Fuc trisaccharide. Therefore, O-Fuc glycans necessary for the modulation of important signaling events in Drosophila are, as in vertebrates, substrates for extension beyond the addition of a single HexNAc.

Maternal Groucho and bHLH repressors amplify the dose-sensitive X chromosome signal in Drosophila sex determination

In Drosophila, XX embryos are fated to develop as females, and XY embryos as males, because the diplo-X dose of four X-linked signal element genes, XSEs, activates the Sex-lethal establishment promoter, SxlPe, whereas the haplo-X XSE dose leaves SxlPe off. The threshold response of SxlPe to XSE concentrations depends in part on the bHLH repressor, Deadpan, present in equal amounts in XX and XY embryos. We identified canonical and non-canonical DNA-binding sites for Dpn at SxlPe and found that cis-acting mutations in the Dpn-binding sites caused stronger and earlier Sxl expression than did deletion of dpn implicating other bHLH repressors in Sxl regulation. Maternal Hey encodes one such bHLH regulator but the E(spl) locus does not. Elimination of the maternal corepressor Groucho also caused strong ectopic Sxl expression in XY, and premature Sxl activation in XX embryos, but Sxl was still expressed differently in the sexes. Our findings suggest that Groucho and associated maternal and zygotic bHLH repressors define the threshold XSE concentrations needed to activate SxlPe and that they participate directly in sex signal amplification. We present a model in which the XSE signal is amplified by a feedback mechanism that interferes with Gro-mediated repression in XX, but not XY embryos.

A Drosophila Smyd4 homologue is a muscle-specific transcriptional modulator involved in development

Background

SET and MYND domain (Smyd) proteins are involved in the transcriptional regulation of cellular proliferation and development in vertebrates. However, the in vivo functions and mechanisms by which these proteins act are poorly understood.

Methodology/Principal Findings

We have used biochemical and genetic approaches to study the role of a Smyd protein in Drosophila. We identified eleven Drosophila genes that encode Smyd proteins. CG14122 encodes a Smyd4 homologue that we have named dSmyd4. dSmyd4 repressed transcription and recruited class I histone deacetylases (HDACs). A region of dSmyd4 including the MYND domain interacted directly with ∼150 amino acids at the N-termini of dHDAC1 and dHDAC3. dSmyd4 interacts selectively with Ebi, a component of the dHDAC3/SMRTER co-repressor complex. During embryogenesis dSmyd4 was expressed throughout the mesoderm, with highest levels in the somatic musculature. Muscle-specific RNAi against dSmyd4 resulted in depletion of the protein and lead to severe lethality. Eclosion is the final moulting stage of Drosophila development when adult flies escape from the pupal case. 80% of dSmyd4 knockdown flies were not able to eclose, resulting in late pupal lethality. However, many aspects of eclosion were still able to occur normally, indicating that dSmyd4 is likely to be involved in the development or function of adult muscle.

Conclusions/Significance

Repression of transcription by dSmyd4 and the involvement of this protein in development suggests that aspects of Smyd protein function are conserved between vertebrates and invertebrates.

Transcription factors bind thousands of active and inactive regions in the Drosophila blastoderm

Identifying the genomic regions bound by sequence-specific regulatory factors is central both to deciphering the complex DNA cis-regulatory code that controls transcription in metazoans and to determining the range of genes that shape animal morphogenesis. We used whole-genome tiling arrays to map sequences bound in Drosophila melanogaster embryos by the six maternal and gap transcription factors that initiate anterior-posterior patterning. We find that these sequence-specific DNA binding proteins bind with quantitatively different specificities to highly overlapping sets of several thousand genomic regions in blastoderm embryos. Specific high- and moderate-affinity in vitro recognition sequences for each factor are enriched in bound regions. This enrichment, however, is not sufficient to explain the pattern of binding in vivo and varies in a context-dependent manner, demonstrating that higher-order rules must govern targeting of transcription factors. The more highly bound regions include all of the over 40 well-characterized enhancers known to respond to these factors as well as several hundred putative new cis-regulatory modules clustered near developmental regulators and other genes with patterned expression at this stage of embryogenesis. The new targets include most of the microRNAs (miRNAs) transcribed in the blastoderm, as well as all major zygotically transcribed dorsal-ventral patterning genes, whose expression we show to be quantitatively modulated by anterior-posterior factors. In addition to these highly bound regions, there are several thousand regions that are reproducibly bound at lower levels. However, these poorly bound regions are, collectively, far more distant from genes transcribed in the blastoderm than highly bound regions; are preferentially found in protein-coding sequences; and are less conserved than highly bound regions. Together these observations suggest that many of these poorly bound regions are not involved in early-embryonic transcriptional regulation, and a significant proportion may be nonfunctional. Surprisingly, for five of the six factors, their recognition sites are not unambiguously more constrained evolutionarily than the immediate flanking DNA, even in more highly bound and presumably functional regions, indicating that comparative DNA sequence analysis is limited in its ability to identify functional transcription factor targets.

Cad74A is regulated by BR and is required for robust dorsal appendage formation in Drosophila oogenesis

Drosophila egg development is an established model for studying epithelial patterning and morphogenesis, but the connection between signaling pathways and egg morphology is still incompletely understood. We have identified a non-classical cadherin, Cad74A, as a putative adhesion gene that bridges epithelial patterning and morphogenesis in the follicle cells. Starting in mid-oogenesis, Cad74A is expressed in the follicle cells that contact the oocyte, including the border cells and most of the columnar follicle cells. However, Cad74A is repressed in two dorsolateral patches of follicle cells, which participate in the formation of tubular respiratory appendages. We show genetically that Cad74A is downstream of the EGFR and BMP signaling pathways and is repressed by the Zn-finger transcription factor Broad. The correlation of Cad74A repression in the cells that bend out of the plane of the follicular epithelium is preserved across Drosophila species and mutant backgrounds exhibiting a range of eggshell phenotypes. Complete removal of Cad74A from the follicle cells causes defects in dorsal appendage formation. Ectopic expression of Cad74A in the roof cells results in shortened, flattened appendages due to the hindered migration of the roof cells. Based on these results, we propose that Cad74A is part of the adhesive machinery that enables robust dorsal appendage formation, and as such provides a link between the patterning of the follicle cells and eggshell morphogenesis.

Prediction of gene expression in embryonic structures of Drosophila melanogaster

Understanding how sets of genes are coordinately regulated in space and time to generate the diversity of cell types that characterise complex metazoans is a major challenge in modern biology. The use of high-throughput approaches, such as large-scale in situ hybridisation and genome-wide expression profiling via DNA microarrays, is beginning to provide insights into the complexities of development. However, in many organisms the collection and annotation of comprehensive in situ localisation data is a difficult and time-consuming task. Here, we present a widely applicable computational approach, integrating developmental time-course microarray data with annotated in situ hybridisation studies, that facilitates the de novo prediction of tissue-specific expression for genes that have no in vivo gene expression localisation data available. Using a classification approach, trained with data from microarray and in situ hybridisation studies of gene expression during Drosophila embryonic development, we made a set of predictions on the tissue-specific expression of Drosophila genes that have not been systematically characterised by in situ hybridisation experiments. The reliability of our predictions is confirmed by literature-derived annotations in FlyBase, by overrepresentation of Gene Ontology biological process annotations, and, in a selected set, by detailed gene-specific studies from the literature. Our novel organism-independent method will be of considerable utility in enriching the annotation of gene function and expression in complex multicellular organisms.

The task of deciphering the complex transcriptional regulatory networks controlling development is one of the major current challenges for molecular biology. The problem is difficult, if not impossible, to solve without a detailed knowledge of the spatiotemporal dynamics of gene expression. Thus, to understand development, we need to identify and functionally characterize all players in regulatory networks. Data on gene expression dynamics obtained from whole transcriptome microarray experiments, combined with in situ hybridization mRNA localisation patterns for a subset of genes, may provide a route for predicting the localisation of gene expression for those genes for which in situ data has not been generated, as well as suggesting functional information for uncharacterised genes. Here, we report the development of one of the first methods for predicting the localisation of gene expression during Drosophila embryogenesis from microarray data. Pooling the subset of genes in the fly genome with in situ data to form functional units, localised in space and time for relevant developmental processes, facilitates the statement of a classification problem, which we address with machine-learning methods. Our approach promotes a richer annotation of biological function for genes in the absence of costly and time-consuming experimental analysis.

Inferring differentiation pathways from gene expression

Motivation: The regulation of proliferation and differentiation of embryonic and adult stem cells into mature cells is central to developmental biology. Gene expression measured in distinguishable developmental stages helps to elucidate underlying molecular processes. In previous work we showed that functional gene modules, which act distinctly in the course of development, can be represented by a mixture of trees. In general, the similarities in the gene expression programs of cell populations reflect the similarities in the differentiation path.

Results: We propose a novel model for gene expression profiles and an unsupervised learning method to estimate developmental similarity and infer differentiation pathways. We assess the performance of our model on simulated data and compare it with favorable results to related methods. We also infer differentiation pathways and predict functional modules in gene expression data of lymphoid development.

Conclusions: We demonstrate for the first time how, in principal, the incorporation of structural knowledge about the dependence structure helps to reveal differentiation pathways and potentially relevant functional gene modules from microarray datasets. Our method applies in any area of developmental biology where it is possible to obtain cells of distinguishable differentiation stages.

Availability: The implementation of our method (GPL license), data and additional results are available at http://algorithmics.molgen.mpg.de/Supplements/InfDif/

Contact:filho@molgen.mpg.de, schliep@molgen.mpg.de

Supplementary information:Supplementary data is available at Bioinformatics online.

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.

Automated annotation of Drosophila gene expression patterns using a controlled vocabulary

MOTIVATION

Regulation of gene expression in space and time directs its localization to a specific subset of cells during development. Systematic determination of the spatiotemporal dynamics of gene expression plays an important role in understanding the regulatory networks driving development. An atlas for the gene expression patterns of fruit fly Drosophila melanogaster has been created by whole-mount in situ hybridization, and it documents the dynamic changes of gene expression pattern during Drosophila embryogenesis. The spatial and temporal patterns of gene expression are integrated by anatomical terms from a controlled vocabulary linking together intermediate tissues developed from one another. Currently, the terms are assigned to patterns manually. However, the number of patterns generated by high-throughput in situ hybridization is rapidly increasing. It is, therefore, tempting to approach this problem by employing computational methods.

RESULTS

In this article, we present a novel computational framework for annotating gene expression patterns using a controlled vocabulary. In the currently available high-throughput data, annotation terms are assigned to groups of patterns rather than to individual images. We propose to extract invariant features from images, and construct pyramid match kernels to measure the similarity between sets of patterns. To exploit the complementary information conveyed by different features and incorporate the correlation among patterns sharing common structures, we propose efficient convex formulations to integrate the kernels derived from various features. The proposed framework is evaluated by comparing its annotation with that of human curators, and promising performance in terms of F1 score has been reported.

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.

Identification of genes influencing dendrite morphogenesis in developing peripheral sensory and central motor neurons

Background

Developing neurons form dendritic trees with cell type-specific patterns of growth, branching and targeting. Dendrites of Drosophila peripheral sensory neurons have emerged as a premier genetic model, though the molecular mechanisms that underlie and regulate their morphogenesis remain incompletely understood. Still less is known about this process in central neurons and the extent to which central and peripheral dendrites share common organisational principles and molecular features. To address these issues, we have carried out two comparable gain-of-function screens for genes that influence dendrite morphologies in peripheral dendritic arborisation (da) neurons and central RP2 motor neurons.

Results

We found 35 unique loci that influenced da neuron dendrites, including five previously shown as required for da dendrite patterning. Several phenotypes were class-specific and many resembled those of known mutants, suggesting that genes identified in this study may converge with and extend known molecular pathways for dendrite development in da neurons. The second screen used a novel technique for cell-autonomous gene misexpression in RP2 motor neurons. We found 51 unique loci affecting RP2 dendrite morphology, 84% expressed in the central nervous system. The phenotypic classes from both screens demonstrate that gene misexpression can affect specific aspects of dendritic development, such as growth, branching and targeting. We demonstrate that these processes are genetically separable. Targeting phenotypes were specific to the RP2 screen, and we propose that dendrites in the central nervous system are targeted to territories defined by Cartesian co-ordinates along the antero-posterior and the medio-lateral axes of the central neuropile. Comparisons between the screens suggest that the dendrites of peripheral da and central RP2 neurons are shaped by regulatory programs that only partially overlap. We focused on one common candidate pathway controlled by the ecdysone receptor, and found that it promotes branching and growth of developing da neuron dendrites, but a role in RP2 dendrite development during embryonic and early larval stages was not apparent.

Conclusion

We identified commonalities (for example, growth and branching) and distinctions (for example, targeting and ecdysone response) in the molecular and organizational framework that underlies dendrite development of peripheral and central neurons.

Bioimage informatics: a new area of engineering biology

In recent years, the deluge of complicated molecular and cellular microscopic images creates compelling challenges for the image computing community. There has been an increasing focus on developing novel image processing, data mining, database and visualization techniques to extract, compare, search and manage the biological knowledge in these data-intensive problems. This emerging new area of bioinformatics can be called ‘bioimage informatics’. This article reviews the advances of this field from several aspects, including applications, key techniques, available tools and resources. Application examples such as high-throughput/high-content phenotyping and atlas building for model organisms demonstrate the importance of bioimage informatics. The essential techniques to the success of these applications, such as bioimage feature identification, segmentation and tracking, registration, annotation, mining, image data management and visualization, are further summarized, along with a brief overview of the available bioimage databases, analysis tools and other resources.

Contact:pengh@janelia.hhmi.org

Supplementary information:Supplementary data are available at Bioinformatics online.

Regulation of the Drosophila lin-41 homologue dappled by let-7 reveals conservation of a regulatory mechanism within the LIN-41 subclade

Drosophila Dappled (DPLD) is a member of the RBCC/TRIM superfamily, a protein family involved in numerous diverse processes such as developmental timing and asymmetric cell divisions. DPLD belongs to the LIN-41 subclade, several members of which are micro RNA (miRNA) regulated. We re-examined the LIN-41 subclade members and their relation to other RBCC/TRIMs and dpld paralogs, and identified a new Drosophila muscle specific RBCC/TRIM: Another B-Box Affiliate, ABBA. In silico predictions of candidate miRNA regulators of dpld identified let-7 as the strongest candidate. Overexpression of dpld led to abnormal eye development, indicating that strict regulation of dpld mRNA levels is crucial for normal eye development. This phenotype was sensitive to let-7 dosage, suggesting let-7 regulation of dpld in the eye disc. A cell-based assay verified let-7 miRNA down-regulation of dpld expression by means of its 3'-untranslated region. Thus, dpld seems also to be miRNA regulated, suggesting that miRNAs represent an ancient mechanism of LIN-41 regulation.

Wiser (tsl): a recessive X-linked temperature-sensitive lethal mutation that affects the wings and the eyes in Drosophila melanogaster

In this study, we characterize a recessive X-linked temperature-sensitive mutation of the gene CG32711. The mutation, named wiser ( tsl ) (wings scalloped-eyes rough), was isolated from a dysgenic cross and is due to a natural P element insertion within the 5' regulatory region of the gene. Mutant (wiser ( tsl )) individuals exhibit wing notching, rough eyes, tarsal malformations and reduced life-span. At 29 degrees C they die at larval and late pharate stages or during eclosion. The CG32711 (wiser) gene is mainly expressed in the ventral midline cells, the peripheral neural system, the hemocytes and the tracheal system of embryos. It is also expressed in nurse cells of adult female ovaries. Our results show that the wiser gene is alternatively spliced generating two mRNAs, which share the same open reading frame, while western analysis identified two protein isoforms. Their expression pattern depends on the stage of development and the culture temperature. wiser ( tsl ) and wild-type individuals display different expression patterns of the two isoforms and this difference most probably accounts for the mutant phenotype. Our results indicate that wiser is a vital gene for the development of Drosophila melanogaster which has no orthologs outside the Drosophilidae.