Diversity and dynamics of the Drosophila transcriptome.

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Diversity of miRNAs, siRNAs, and piRNAs across 25 Drosophila cell lines.

We expanded the knowledge base for Drosophila cell line transcriptomes by deeply sequencing their small RNAs. In total, we analyzed more than 1 billion raw reads from 53 libraries across 25 cell lines. We verify reproducibility of biological replicate data sets, determine common and distinct aspects of miRNA expression across cell lines, and infer the global impact of miRNAs on cell line transcriptomes. We next characterize their commonalities and differences in endo-siRNA populations. Interestingly, most cell lines exhibit enhanced TE-siRNA production relative to tissues, suggesting this as a common aspect of cell immortalization. We also broadly extend annotations of cis-NAT-siRNA loci, identifying ones with common expression across diverse cells and tissues, as well as cell-restricted loci. Finally, we characterize small RNAs in a set of ovary-derived cell lines, including somatic cells (OSS and OSC) and a mixed germline/somatic cell population (fGS/OSS) that exhibits ping-pong piRNA signatures. Collectively, the ovary data reveal new genic piRNA loci, including unusual configurations of piRNA-generating regions. Together with the companion analysis of mRNAs described in a previous study, these small RNA data provide comprehensive information on the transcriptional landscape of diverse Drosophila cell lines. These data should encourage broader usage of fly cell lines, beyond the few that are presently in common usage.

The developmental transcriptome of Drosophila melanogaster.

Drosophila melanogaster is one of the most well studied genetic model organisms; nonetheless, its genome still contains unannotated coding and non-coding genes, transcripts, exons and RNA editing sites. Full discovery and annotation are pre-requisites for understanding how the regulation of transcription, splicing and RNA editing directs the development of this complex organism. Here we used RNA-Seq, tiling microarrays and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events, and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. These data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.

Diversification of doublesex function underlies morph-, sex-, and species-specific development of beetle horns.

Sex-specific trait expression is frequently associated with highly variable, condition-dependent expression within sexes and rapid divergence among closely related species. Horned beetles are an excellent example for studying the molecular basis of these phenomena because horn morphology varies markedly among species, between sexes, and among alternative, nutritionally-cued morphs within sexes. In addition, horns lack obvious homology to other insect traits and provide a good opportunity to explore the molecular basis of the rapid diversification of a novel trait within and between species. Here we show that the sex-determination gene doublesex (dsx) underlies important aspects of horn development, including differences between sexes, morphs, and species. In male Onthophagus taurus, dsx transcripts were preferentially expressed in the horns of the large, horned morph, and RNAi-mediated knockdown of dsx dramatically altered male horn allometry by massively reducing horn development in large males, but not in smaller males. Conversely, dsx RNAi induced ectopic, nutrition-sensitive horn development in otherwise hornless females. Finally, in a closely related species (Onthophagus sagittarius) that has recently evolved a rare reversed sexual dimorphism, dsx RNAi revealed reversed as well as novel dsx functions despite an overall conservation of dsx expression. This suggests that rapid evolution of dsx functions has facilitated the transition from a regular sexual dimorphism to a reversed sexual dimorphism in this species. Our findings add beetle horns to existing examples of a close relationship between dsx and sexual trait development, and suggest that dsx function has been coopted to facilitate both the evolution of environmentally-cued intrasexual dimorphisms and rapid species divergences in a novel trait.

The spread of a transposon insertion in Rec8 is associated with obligate asexuality in Daphnia.

Although transitions from sexual to asexual reproduction are thought to have important evolutionary consequences, little is known about the mechanistic underpinnings of these changes. The cyclical parthenogen Daphnia pulex is a powerful model in which to address these issues because female-limited meiosis suppression can be transmitted to sexual individuals via males, providing the opportunity for genetic dissection of the trait. A previous study identified genomic regions differentiating obligately asexual females from their sexual counterparts, and a candidate gene within one such region, encoding the meiotic cohesin Rec8, is the subject of this investigation. The D. pulex genome contains three Rec8 loci, all of which are quite polymorphic. However, at one of the loci, all obligately asexual clones carry an allele containing an identical upstream insertion of a transposable element as well as a frameshift mutation, both of which are completely absent from sexual lineages. The low level of variation within the insertion allele across all asexual lineages suggests that this element may be in the process of spreading through the species, and abrogation or modification of Rec8 function is possibly responsible for converting meiotically reproducing lineages into obligate asexuals.

Genome-wide analysis of promoter architecture in Drosophila melanogaster.

Core promoters are critical regions for gene regulation in higher eukaryotes. However, the boundaries of promoter regions, the relative rates of initiation at the transcription start sites (TSSs) distributed within them, and the functional significance of promoter architecture remain poorly understood. We produced a high-resolution map of promoters active in the Drosophila melanogaster embryo by integrating data from three independent and complementary methods: 21 million cap analysis of gene expression (CAGE) tags, 1.2 million RNA ligase mediated rapid amplification of cDNA ends (RLM-RACE) reads, and 50,000 cap-trapped expressed sequence tags (ESTs). We defined 12,454 promoters of 8037 genes. Our analysis indicates that, due to non-promoter-associated RNA background signal, previous studies have likely overestimated the number of promoter-associated CAGE clusters by fivefold. We show that TSS distributions form a complex continuum of shapes, and that promoters active in the embryo and adult have highly similar shapes in 95% of cases. This suggests that these distributions are generally determined by static elements such as local DNA sequence and are not modulated by dynamic signals such as histone modifications. Transcription factor binding motifs are differentially enriched as a function of promoter shape, and peaked promoter shape is correlated with both temporal and spatial regulation of gene expression. Our results contribute to the emerging view that core promoters are functionally diverse and control patterning of gene expression in Drosophila and mammals.

The transcriptional diversity of 25 Drosophila cell lines.

Drosophila melanogaster cell lines are important resources for cell biologists. Here, we catalog the expression of exons, genes, and unannotated transcriptional signals for 25 lines. Unannotated transcription is substantial (typically 19% of euchromatic signal). Conservatively, we identify 1405 novel transcribed regions; 684 of these appear to be new exons of neighboring, often distant, genes. Sixty-four percent of genes are expressed detectably in at least one line, but only 21% are detected in all lines. Each cell line expresses, on average, 5885 genes, including a common set of 3109. Expression levels vary over several orders of magnitude. Major signaling pathways are well represented: most differentiation pathways are "off" and survival/growth pathways "on." Roughly 50% of the genes expressed by each line are not part of the common set, and these show considerable individuality. Thirty-one percent are expressed at a higher level in at least one cell line than in any single developmental stage, suggesting that each line is enriched for genes characteristic of small sets of cells. Most remarkable is that imaginal disc-derived lines can generally be assigned, on the basis of expression, to small territories within developing discs. These mappings reveal unexpected stability of even fine-grained spatial determination. No two cell lines show identical transcription factor expression. We conclude that each line has retained features of an individual founder cell superimposed on a common "cell line" gene expression pattern.

Fine scale analysis of gene expression in Drosophila melanogaster gonads reveals Programmed cell death 4 promotes the differentiation of female germline stem cells.

BACKGROUND: Germline stem cells (GSCs) are present in the gonads of Drosophila females and males, and their proper maintenance, as well as their correct differentiation, is essential for fertility and fecundity. The molecular characterization of factors involved in maintenance and differentiation is a major goal both in Drosophila and stem cell research. While genetic studies have identified many of these key factors, the use of genome-wide expression studies holds the potential to greatly increase our knowledge of these pathways. RESULTS: Here we report a genome-wide expression study that uses laser cutting microdissection to isolate germline stem cells, somatic niche cells, and early differentiating germ cells from female and male gonads. Analysis of this data, in association with two previously published genome-wide GSC data sets, revealed sets of candidate genes as putatively expressed in specific cell populations. Investigation of one of these genes, CG10990 the Drosophila ortholog of mammalian Programmed cell death 4 (Pdcd4), reveals expression in female and male germline stem cells and early differentiating daughter cells. Functional analysis demonstrates that while it is not essential for oogenesis or spermatogenesis, it does function to promote the differentiation of GSCs in females. Furthermore, in females, Pdcd4 genetically interacts with the key differentiation gene bag of marbles (bam) and the stem cell renewal factor eIF4A, suggesting a possible pathway for its function in differentiation. CONCLUSIONS: We propose that Pdcd4 promotes the differentiation of GSC daughter cells by relieving the eIF4A-mediated inhibition of Bam.

Systems-biology dissection of eukaryotic cell growth.

A recent article in BMC Biology illustrates the use of a systems-biology approach to integrate data across the transcriptome, proteome and metabolome of budding yeast in order to dissect the relationship between nutrient conditions and cell growth.

Gene discovery in the horned beetle Onthophagus taurus.

BACKGROUND: Horned beetles, in particular in the genus Onthophagus, are important models for studies on sexual selection, biological radiations, the origin of novel traits, developmental plasticity, biocontrol, conservation, and forensic biology. Despite their growing prominence as models for studying both basic and applied questions in biology, little genomic or transcriptomic data are available for this genus. We used massively parallel pyrosequencing (Roche 454-FLX platform) to produce a comprehensive EST dataset for the horned beetle Onthophagus taurus. To maximize sequence diversity, we pooled RNA extracted from a normalized library encompassing diverse developmental stages and both sexes. RESULTS: We used 454 pyrosequencing to sequence ESTs from all post-embryonic stages of O. taurus. Approximately 1.36 million reads assembled into 50,080 non-redundant sequences encompassing a total of 26.5 Mbp. The non-redundant sequences match over half of the genes in Tribolium castaneum, the most closely related species with a sequenced genome. Analyses of Gene Ontology annotations and biochemical pathways indicate that the O. taurus sequences reflect a wide and representative sampling of biological functions and biochemical processes. An analysis of sequence polymorphisms revealed that SNP frequency was negatively related to overall expression level and the number of tissue types in which a given gene is expressed. The most variable genes were enriched for a limited number of GO annotations whereas the least variable genes were enriched for a wide range of GO terms directly related to fitness. CONCLUSIONS: This study provides the first large-scale EST database for horned beetles, a much-needed resource for advancing the study of these organisms. Furthermore, we identified instances of gene duplications and alternative splicing, useful for future study of gene regulation, and a large number of SNP markers that could be used in population-genetic studies of O. taurus and possibly other horned beetles.

Programed cell death shapes the expression of horns within and between species of horned beetles.

Holometabolous insects provide an excellent opportunity to study both the properties of development as well as their evolution and diversification across taxa. Here we investigate the developmental basis and evolutionary diversification of secondary trait loss during development in the expression of beetle horns, a novel and highly diverse class of secondary sexual traits. In many species, horn growth during late larval development is followed by a period of dramatic remodeling during the pupal stage, including the complete resorption of horns in many cases. Here we show that programed cell death plays an important and dynamic role in the secondary resorption of pupal horn primordia during pupal development. Surprisingly, the degree of cell death mediated horn resorption depended on species, sex, and body region, suggesting the existence of regulatory mechanisms that can diversify quickly over short phylogenetic distances. More generally, our results illustrate that secondary, differential loss of structures during development can be a powerful mechanism for generating considerable morphological diversity both within and between species.

EST and microarray analysis of horn development in Onthophagus beetles.

BACKGROUND: The origin of novel traits and their subsequent diversification represent central themes in evo-devo and evolutionary ecology. Here we explore the genetic and genomic basis of a class of traits that is both novel and highly diverse, in a group of organisms that is ecologically complex and experimentally tractable: horned beetles. RESULTS: We developed two high quality, normalized cDNA libraries for larval and pupal Onthophagus taurus and sequenced 3,488 ESTs that assembled into 451 contigs and 2,330 singletons. We present the annotation and a comparative analysis of the conservation of the sequences. Microarrays developed from the combined libraries were then used to contrast the transcriptome of developing primordia of head horns, prothoracic horns, and legs. Our experiments identify a first comprehensive list of candidate genes for the evolution and diversification of beetle horns. We find that developing horns and legs show many similarities as well as important differences in their transcription profiles, suggesting that the origin of horns was mediated partly, but not entirely, by the recruitment of genes involved in the formation of more traditional appendages such as legs. Furthermore, we find that horns developing from the head and prothorax differ in their transcription profiles to a degree that suggests that head and prothoracic horns are not serial homologs, but instead may have evolved independently from each other. CONCLUSION: We have laid the foundation for a systematic analysis of the genetic basis of horned beetle development and diversification with the potential to contribute significantly to several major frontiers in evolutionary developmental biology.

A machine-learning approach to combined evidence validation of genome assemblies.

MOTIVATION: While it is common to refer to 'the genome sequence' as if it were a single, complete and contiguous DNA string, it is in fact an assembly of millions of small, partially overlapping DNA fragments. Sophisticated computer algorithms (assemblers and scaffolders) merge these DNA fragments into contigs, and place these contigs into sequence scaffolds using the paired-end sequences derived from large-insert DNA libraries. Each step in this automated process is susceptible to producing errors; hence, the resulting draft assembly represents (in practice) only a likely assembly that requires further validation. Knowing which parts of the draft assembly are likely free of errors is critical if researchers are to draw reliable conclusions from the assembled sequence data. RESULTS: We develop a machine-learning method to detect assembly errors in sequence assemblies. Several in silico measures for assembly validation have been proposed by various researchers. Using three benchmarking Drosophila draft genomes, we evaluate these techniques along with some new measures that we propose, including the good-minus-bad coverage (GMB), the good-to-bad-ratio (RGB), the average Z-score (AZ) and the average absolute Z-score (ASZ). Our results show that the GMB measure performs better than the others in both its sensitivity and its specificity for assembly error detection. Nevertheless, no single method performs sufficiently well to reliably detect genomic regions requiring attention for further experimental verification. To utilize the advantages of all these measures, we develop a novel machine learning approach that combines these individual measures to achieve a higher prediction accuracy (i.e. greater than 90%). Our combined evidence approach avoids the difficult and often ad hoc selection of many parameters the individual measures require, and significantly improves the overall precisions on the benchmarking data sets.

Correction to: DNA copy number evolution in Drosophila cell lines.

Following publication of the original article [1], the authors reported the following errors.

Profiling sex-biased gene expression during parthenogenetic reproduction in Daphnia pulex.

BACKGROUND: Sexual reproduction is a core biological function that is conserved throughout eukaryotic evolution, yet breeding systems are extremely variable. Genome-wide comparative studies can be effectively used to identify genes and regulatory patterns that are constrained to preserve core functions from those that may help to account for the diversity of animal reproductive strategies. We use a custom microarray to investigate gene expression in males and two reproductive stages of females in the crustacean Daphnia pulex. Most Daphnia species reproduce by cyclical parthenogenesis, alternating between sexual and clonal reproduction. Both sex determination and the switch in their mode of reproduction is environmentally induced, making Daphnia an interesting comparative system for the study of sex-biased and reproductive genes. RESULTS: Patterns of gene expression in females and males reveal that 50% of assayed transcripts show some degree of sex-bias. Female-biased transcription is enriched for translation, metabolic and regulatory genes associated with development. Male-biased expression is enriched for cuticle and protease function. Comparison with well studied arthropods such as Drosophila melanogaster and Anopheles gambiae suggests that female-biased patterns tend to be conserved, whereas male-biased genes are evolving faster in D. pulex. These findings are based on the proportion of female-biased, male-biased, and unbiased genes that share sequence similarity with proteins in other animal genomes. CONCLUSION: Some transcriptional differences between males and females appear to be conserved across Arthropoda, including the rapid evolution of male-biased genes which is observed in insects and now in a crustacean. Yet, novel patterns of male-biased gene expression are also uncovered. This study is an important first step towards a detailed understanding of the genetic basis and evolution of parthenogenesis, environmental sex determination, and adaptation to aquatic environments.

Sampling Daphnia's expressed genes: preservation, expansion and invention of crustacean genes with reference to insect genomes.

BACKGROUND: Functional and comparative studies of insect genomes have shed light on the complement of genes, which in part, account for shared morphologies, developmental programs and life-histories. Contrasting the gene inventories of insects to those of the nematodes provides insight into the genomic changes responsible for their diversification. However, nematodes have weak relationships to insects, as each belongs to separate animal phyla. A better outgroup to distinguish lineage specific novelties would include other members of Arthropoda. For example, crustaceans are close allies to the insects (together forming Pancrustacea) and their fascinating aquatic lifestyle provides an important comparison for understanding the genetic basis of adaptations to life on land versus life in water. RESULTS: This study reports on the first characterization of cDNA libraries and sequences for the model crustacean Daphnia pulex. We analyzed 1,546 ESTs of which 1,414 represent approximately 787 nuclear genes, by measuring their sequence similarities with insect and nematode proteomes. The provisional annotation of genes is supported by expression data from microarray studies described in companion papers. Loci expected to be shared between crustaceans and insects because of their mutual biological features are identified, including genes for reproduction, regulation and cellular processes. We identify genes that are likely derived within Pancrustacea or lost within the nematodes. Moreover, lineage specific gene family expansions are identified, which suggest certain biological demands associated with their ecological setting. In particular, up to seven distinct ferritin loci are found in Daphnia compared to three in most insects. Finally, a substantial fraction of the sampled gene transcripts shares no sequence similarity with those from other arthropods. Genes functioning during development and reproduction are comparatively well conserved between crustaceans and insects. By contrast, genes that were responsive to environmental conditions (metal stress) and not sex-biased included the greatest proportion of genes with no matches to insect proteomes. CONCLUSION: This study along with associated microarray experiments are the initial steps in a coordinated effort by the Daphnia Genomics Consortium to build the necessary genomic platform needed to discover genes that account for the phenotypic diversity within the genus and to gain new insights into crustacean biology. This effort will soon include the first crustacean genome sequence.

Automated liquid handling and high-throughput preparation of polymerase chain reaction-amplified DNA for microarray fabrication.

Genome-wide studies of gene expression and transcription factor-binding sites using DNA microarrays are leading to new systems level insights. The massively parallel nature of microarrays presents technical challenges: fabricating high-quality microarrays at the front end and data analysis and interpretation downstream. A principal challenge in fabricating microarrays is preparation of the DNA samples. This is particularly the case for polymerase chain reaction-amplified DNA samples. The challenge is to scale up efficiently to high-throughput preparation of tens of thousands of DNA samples while ensuring a uniform high quality. This chapter outlines strategic considerations, including automated liquid handling and workflow development to maximize efficiency, and quality control (QC) measures to ensure uniform quality. The protocols are presented with commentary to illustrate their logic and specific techniques. These principles and techniques are extensible to other high-throughput molecular biological applications.

Troubleshooting microarray hybridizations.

Microarray experiments are being performed more widely than ever before, but even seasoned investigators can experience technical problems with hybridizations. This chapter provides guidelines for recognizing, rectifying, and avoiding common trouble areas. Specifically, it addresses frequent complications related to artifacts of printing, RNA sample preparation and quality, fluorophore labeling, hybridization conditions, and posthybridization washes. Emphasis is placed on investigating problems though a combination of appropriate controls and image analysis, where diagnostic plots of data quality are used to illustrate characteristics of acceptable and unsatisfactory hybridizations. This chapter also discusses resources available to microarray users hoping to improve the sensitivity and specificity of their experiments.

Sex determination signals control ovo-B transcription in Drosophila melanogaster germ cells.

Nonautonomous inductive signals from the soma and autonomous signals due to a 2X karyotype determine the sex of Drosophila melanogaster germ cells. These two signals have partially overlapping influences on downstream sex determination genes. The upstream OVO-B transcription factor is required for the viability of 2X germ cells, regardless of sexual identity, and for female germline sexual identity. The influence of inductive and autonomous signals on ovo expression has been controversial. We show that ovo-B is strongly expressed in the 2X germ cells in either a male or a female soma. This indicates that a 2X karyotype controls ovo-B expression in the absence of inductive signals from the female soma. However, we also show that female inductive signals positively regulate ovo-B transcription in the 1X germ cells that do not require ovo-B function. Genetic analysis clearly indicates that inductive signals from the soma are not required for ovo-B function in 2X germ cells. Thus, while somatic inductive signals and chromosome karyotype have overlapping regulatory influences, a 2X karyotype is a critical germline autonomous determinant of ovo-B function in the germline.

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.