The Drosophila melanogaster transcriptome by paired-end RNA sequencing.

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 356× 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.

Rapid identification of heterozygous mutations in Drosophila melanogaster using genomic capture 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-genome-sequencing 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.

TRAFfic signals: High-throughput CAR discovery in NK cells reveals novel TRAF-binding endodomains that drive enhanced persistence and cytotoxicity.

Natural killer (NK) cells are a promising alternative therapeutic platform to CAR T cells given their favorable safety profile and potent killing ability. However, CAR NK cells suffer from limited persistence in vivo , which is, in part, thought to be the consequence of limited cytokine signaling. To address this challenge, we developed an innovative high-throughput screening strategy to identify CAR endodomains that could drive enhanced persistence while maintaining potent cytotoxicity. We uncovered a family of TRAF-binding endodomains that outperform benchmarks in primary NK cells along dimensions of persistence and cytotoxicity, even in low IL-2 conditions. This work highlights the importance of cell-type-specific cell therapy engineering and unlocks a wide range of high-throughput molecular engineering avenues in NK cells.

Global Gene Expression Analysis Reveals Complex Cuticle Organization of the Tribolium Compound Eye.

The red flour beetle Tribolium castaneum is a resource-rich model for genomic and developmental studies. To extend previous studies on Tribolium eye development, we produced transcriptomes for normal-eyed and eye-depleted heads of pupae and adults to identify differentially transcript-enriched (DE) genes in the visual system. Unexpectedly, cuticle-related genes were the largest functional class in the pupal compound eye DE gene population, indicating differential enrichment in three distinct cuticle components: clear lens facet cuticle, highly melanized cuticle of the ocular diaphragm, which surrounds the Tribolium compound eye for internal fortification, and newly identified facet margins of the tanned cuticle, possibly enhancing external fortification. Phylogenetic, linkage, and high-throughput gene knockdown data suggest that most cuticle proteins (CPs) expressed in the Tribolium compound eye stem from the deployment of ancient CP genes. Consistent with this, TcasCPR15, which we identified as the major lens CP gene in Tribolium, is a beetle-specific but pleiotropic paralog of the ancient CPR RR-2 CP gene family. The less abundant yet most likely even more lens-specific TcasCP63 is a member of a sprawling family of noncanonical CP genes, documenting a role of local gene family expansions in the emergence of the Tribolium compound eye CP repertoire. Comparisons with Drosophila and the mosquito Anopheles gambiae reveal a steady turnover of lens-enriched CP genes during insect evolution.

Phototransduction and clock gene expression in the troglobiont beetle Ptomaphagus hirtus of Mammoth cave.

Obligatory 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.

Evolutionary Trace Annotation Server: automated enzyme function prediction in protein structures using 3D templates.

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 function-associated 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 (http://mammoth.bcm.tmc.edu/eta/).

High-throughput multiplex sequencing to discover copy number variants in 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.