elPrep: high-performance preparation of sequence alignment/map files for variant calling

elPrep is a high-performance tool for preparing sequence alignment/map files for variant calling in sequencing pipelines. It can be used as a replacement for SAMtools and Picard for preparation steps such as filtering, sorting, marking duplicates, reordering contigs, and so on, while producing identical results. What sets elPrep apart is its software architecture that allows executing preparation pipelines by making only a single pass through the data, no matter how many preparation steps are used in the pipeline. elPrep is designed as a multithreaded application that runs entirely in memory, avoids repeated file I/O, and merges the computation of several preparation steps to significantly speed up the execution time. For example, for a preparation pipeline of five steps on a whole-exome BAM file (NA12878), we reduce the execution time from about 1: 40 hours, when using a combination of SAMtools and Picard, to about 15 minutes when using elPrep, while utilising the same server resources, here 48 threads and 23GB of RAM. For the same pipeline on whole-genome data (NA12878), elPrep reduces the runtime from 24 hours to less than 5 hours. As a typical clinical study may contain sequencing data for hundreds of patients, elPrep can remove several hundreds of hours of computing time, and thus substantially reduce analysis time and cost

Targeted knock-down of miR21 primary transcripts using snoMEN vectors induces apoptosis in human cancer cell lines

We have previously reported an antisense technology, 'snoMEN vectors', for targeted knock-down of protein coding mRNAs using human snoRNAs manipulated to contain short regions of sequence complementarity with the mRNA target. Here we characterise the use of snoMEN vectors to target the knock-down of micro RNA primary transcripts. We document the specific knock-down of miR21 in HeLa cells using plasmid vectors expressing miR21-targeted snoMEN RNAs and show this induces apoptosis. Knock-down is dependent on the presence of complementary sequences in the snoMEN vector and the induction of apoptosis can be suppressed by over-expression of miR21. Furthermore, we have also developed lentiviral vectors for delivery of snoMEN RNAs and show this increases the efficiency of vector transduction in many human cell lines that are difficult to transfect with plasmid vectors. Transduction of lentiviral vectors expressing snoMEN targeted to pri-miR21 induces apoptosis in human lung adenocarcinoma cells, which express high levels of miR21, but not in human primary cells. We show that snoMEN-mediated suppression of miRNA expression is prevented by siRNA knock-down of Ago2, but not by knock-down of Ago1 or Upf1. snoMEN RNAs colocalise with Ago2 in cell nuclei and nucleoli and can be co-immunoprecipitated from nuclear extracts by antibodies specific for Ago2

Transcriptional responses of ecologically diverse drosophila species to larval diets differing in relative sugar and protein ratios

We utilized three ecologically diverse Drosophila species to explore the influence of ecological adaptation on transcriptomic responses to isocaloric diets differing in their relative proportions of protein to sugar. Drosophila melanogaster, a cosmopolitan species that breeds in decaying fruit, exemplifies individuals long exposed to a Western diet higher in sugar, while the natural diet of the cactophilic D. mojavensis, is much lower in carbohydrates. Drosophila arizonae, the sister species of D. mojavensis, is largely cactophilic, but also utilizes rotting fruits that are higher in sugars than cacti. We exposed third instar larvae for 24 hours to diets either (1) high in protein relative to sugar, (2) diets with equal amounts of protein and sugar, and (3) diets low in protein but high in sugar. As we predicted, based upon earlier interspecific studies of development and metabolism, the most extreme differences in gene expression under different dietary conditions were found in D. mojavensis followed by D. arizonae. No differential expression among diets was observed for D. melanogaster, a species that survives well under all three conditions, with little impact on its metabolism. We suggest that these three species together provide a model to examine individual and population differences in vulnerability to lifestyle-associated health problems such as metabolic syndrome and diabetes

The Participation of Calponin in the Cross Talk between 20-Hydroxyecdysone and Juvenile Hormone Signaling Pathways by Phosphorylation Variation

20-hydroxyecdysone (20E) and juvenile hormone (JH) signaling pathways interact to mediate insect development, but the mechanism of this interaction is poorly understood. Here, a calponin homologue domain (Chd) containing protein (HaCal) is reported to play a key role in the cross talk between 20E and JH signaling by varying its phosphorylation. Chd is known as an actin binding domain present in many proteins including some signaling proteins. Using an epidermal cell line (HaEpi), HaCal was found to be up-regulated by either 20E or the JH analog methoprene (JHA). 20E induced rapid phosphorylation of HaCal whereas no phosphorylation occurred with JHA. HaCal could be quickly translocated into the nuclei through 20E or JH signaling but interacted with USP1 only under the mediation of JHA. Knockdown of HaCal by RNAi blocked the 20E inducibility of USP1, PKC and HR3, and also blocked the JHA inducibility of USP1, PKC and JHi. After gene silencing of HaCal by ingestion of dsHaCal expressed by Escherichia coli, the larval development was arrested and the gene expression of USP1, PKC, HR3 and JHi were blocked. These composite data suggest that HaCal plays roles in hormonal signaling by quickly transferring into nucleus to function as a phosphorylated form in the 20E pathway and as a non-phosphorylated form interacting with USP1 in the JH pathway to facilitate 20E or JH signaling cascade, in short, by switching its phosphorylation status to regulate insect development

A Novel Ecdysone Receptor Mediates Steroid-Regulated Developmental Events during the Mid-Third Instar of Drosophila

The larval salivary gland of Drosophila melanogaster synthesizes and secretes glue glycoproteins that cement developing animals to a solid surface during metamorphosis. The steroid hormone 20-hydroxyecdysone (20E) is an essential signaling molecule that modulates most of the physiological functions of the larval gland. At the end of larval development, it is known that 20E—signaling through a nuclear receptor heterodimer consisting of EcR and USP—induces the early and late puffing cascade of the polytene chromosomes and causes the exocytosis of stored glue granules into the lumen of the gland. It has also been reported that an earlier pulse of hormone induces the temporally and spatially specific transcriptional activation of the glue genes; however, the receptor responsible for triggering this response has not been characterized. Here we show that the coordinated expression of the glue genes midway through the third instar is mediated by 20E acting to induce genes of the Broad Complex (BRC) through a receptor that is not an EcR/USP heterodimer. This result is novel because it demonstrates for the first time that at least some 20E-mediated, mid-larval, developmental responses are controlled by an uncharacterized receptor that does not contain an RXR-like component

Identification and Characterization of Two Novel RNA Viruses from Anopheles gambiae Species Complex Mosquitoes

Mosquitoes of the Anopheles gambiae complex display strong preference for human blood-meals and are major malaria vectors in Africa. However, their interaction with viruses or role in arbovirus transmission during epidemics has been little examined, with the exception of O'nyong-nyong virus, closely related to Chikungunya virus. Deep-sequencing has revealed different RNA viruses in natural insect viromes, but none have been previously described in the Anopheles gambiae species complex. Here, we describe two novel insect RNA viruses, a Dicistrovirus and a Cypovirus, found in laboratory colonies of An. gambiae taxa using small-RNA deep sequencing. Sequence analysis was done with Metavisitor, an open-source bioinformatic pipeline for virus discovery and de novo genome assembly. Wild-collected Anopheles from Senegal and Cambodia were positive for the Dicistrovirus and Cypovirus, displaying high sequence identity to the laboratory-derived virus. Thus, the Dicistrovirus (Anopheles C virus, AnCV) and Cypovirus (Anopheles Cypovirus, AnCPV) are components of the natural virome of at least some anopheline species. Their possible influence on mosquito immunity or transmission of other pathogens is unknown. These natural viruses could be developed as models for the study of Anopheles-RNA virus interactions in low security laboratory settings, in an analogous manner to the use of rodent malaria parasites for studies of mosquito anti-parasite immunity