RNA-seq was used to generate an extensive map of the Drosophila melanogaster transcriptome by broad sampling of 10 developmental stages. In total, 142.2 million uniquely mapped 64-100-bp paired-end reads were generated on the Illumina GA II yielding 3563 sequencing coverage. More than 95% of FlyBase genes and 90% of splicing junctions were observed. Modifications to 30% of FlyBase gene models were made by extension of untranslated regions, inclusion of novel exons, and identification of novel splicing events. A total of 319 novel transcripts were identified, representing a 2% increase over the current annotation. Alternate splicing was observed in 31% of D. melanogaster genes, a 38% increase over previous estimations, but significantly less than that observed in higher organisms. Much of this splicing is subtle such as tandem alternate splice sites.

SUMMARYObligatory cave species exhibit dramatic trait modifications such as eye reduction, loss of pigmentation and an increase in touch receptors. As molecular studies of cave adaptation have largely concentrated on vertebrate models, it is not yet possible to probe for genetic universalities underlying cave adaptation. We have therefore begun to study the strongly cave-adapted small carrion beetle Ptomaphagus hirtus. For over 100 years, this flightless signature inhabitant of Mammoth Cave, the world's largest known cave system, has been considered blind despite the presence of residual lens structures. By deep sequencing of the adult head transcriptome, we discovered the transcripts of all core members of the phototransduction protein machinery. Combined with the absence of transcripts of select structural photoreceptor and eye pigmentation genes, these data suggest a reduced but functional visual system in P. hirtus. This conclusion was corroborated by a negative phototactic response of P. hirtus in light/dark choice tests. We further detected the expression of the complete circadian clock gene network in P. hirtus, raising the possibility of a role of light sensation in the regulation of oscillating processes. We speculate that P. hirtus is representative of a large number of animal species with highly reduced but persisting visual capacities in the twilight zone of the subterranean realm. These can now be studied on a broad comparative scale given the efficiency of transcript discovery by next-generation sequencing.

One of the key advantages of using Drosophila melanogaster as a genetic model organism is the ability to conduct saturation mutagenesis screens to identify genes and pathways underlying a given phenotype. Despite the large number of genetic tools developed to facilitate downstream cloning of mutations obtained from such screens, the current procedure remains labor intensive, time consuming, and costly. To address this issue, we designed an efficient strategy for rapid identification of heterozygous mutations in the fly genome by combining rough genetic mapping, targeted DNA capture, and second generation sequencing technology. We first tested this method on heterozygous flies carrying either a previously characterized dac 5 or sens E2 mutation. Targeted amplification of genomic regions near these two loci was used to enrich DNA for sequencing, and both point mutations were successfully identified. When this method was applied to uncharacterized twr mutant flies, the underlying mutation was identified as a single-base mutation in the gene Spase18-21. This targeted-genomesequencing method reduces time and effort required for mutation cloning by up to 80% compared with the current approach and lowers the cost to <$1000 for each mutant. Introduction of this and other sequencing-based methods for mutation cloning will enable broader usage of forward genetics screens and have significant impacts in the field of model organisms such as Drosophila.

Copy number variation (CNV) contributes in phenotypically relevant ways to the genetic variability of many organisms. Cost-effective genomewide methods for identifying copy number variation are necessary to elucidate the contribution that these structural variants make to the genomes of model organisms. We have developed a novel approach for the identification of copy number variation by next generation sequencing. As a proof of concept our method has been applied to map the deletions of three Drosophila deficiency strains. We demonstrate that low sequence coverage is sufficient for identifying and mapping large deletions at kilobase resolution, suggesting that data generated from high-throughput sequencing experiments are sufficient for simultaneously analyzing many strains. Genomic DNA from two Drosophila deficiency stocks was barcoded and sequenced in multiplex, and the breakpoints associated with each deletion were successfully identified. The approach we describe is immediately applicable to the systematic exploration of copy number variation in model organisms and humans. S TRUCTURAL variation is known to contribute extensively to the genetic variability of humans, mammals, and many model organisms. One class of structural variant, termed copy number variation (CNV), includes deletions, duplications, insertions, and genomic rearrangements which affect the number of occurrences of a specific DNA sequence present in the genome . CNV is known to occur extensively in the Drosophila genome with functionally significant consequences (Bridges 1936;Dopman and Hartl 2007;Tibshirani and Wang 2008;Zhou et al. 2008). In one study of 15 Drosophila strains, as many as 10% of genes were observed to harbor CNVs (Emerson et al. 2008). Cryptic CNVs that affect the phenotype observed in a model organism have the potential to confound research on multiple levels. For example, a recent report indicates that terminal deletions on chromosome (chr) 2L are frequent among deficiency kit stocks with mutations on the second chromosome and that the associated deletion of lgl has distorted the results of several previous studies (Roegiers et al. 2009). Despite widespread existence of CNV, the biological consequences of this phenomenon remain largely unexplored due to the lack of efficient tools for detection and characterization.Until recently, comparative genomic hybridization with whole-genome tiling arrays (array-CGH) was the primary method for characterizing CNVs (Carter 2007); however, several limitations for this platform reduce its efficacy and efficiency. First, cross-hybridization and reliance on intensity scores lead to data that are difficult to interpret. Second, custom array design and optimization is labor intensive and costly. Third, array-CGH methods can only detect CNV, not other complex rearrangements such as balanced translocations and inversions. Finally, the overall cost of array-CGH methods is relatively high, particularly when high-resolution, whole-genome tiling arrays are employed.Direct sequencing usin...

Summary:The Evolutionary Trace Annotation (ETA) Server predicts enzymatic activity. ETA starts with a structure of unknown function, such as those from structural genomics, and with no prior knowledge of its mechanism uses the phylogenetic Evolutionary Trace (ET) method to extract key functional residues and propose a functionassociated 3D motif, called a 3D template. ETA then searches previously annotated structures for geometric template matches that suggest molecular and thus functional mimicry. In order to maximize the predictive value of these matches, ETA next applies distinctive specificity filters-evolutionary similarity, function plurality and match reciprocity. In large scale controls on enzymes, prediction coverage is 43% but the positive predictive value rises to 92%, thus minimizing false annotations. Users may modify any search parameter, including the template. ETA thus expands the ET suite for protein structure annotation, and can contribute to the annotation efforts of metaservers. Availability: The ETA Server is a web application available at