Differential Expression of Salivary Proteins between Susceptible and Insecticide-Resistant Mosquitoes of Culex quinquefasciatus

Background: The Culex quinquefasciatus mosquito, a major pest and vector of filariasis and arboviruses in the tropics, has developed multiple resistance mechanisms to the main insecticide classes currently available in public health. Among them, the insensitive acetylcholinesterase (ace-1(R) allele) is widespread worldwide and confers cross-resistance to organophosphates and carbamates. Fortunately, in an insecticide-free environment, this mutation is associated with a severe genetic cost that can affect various life history traits. Salivary proteins are directly involved in human-vector contact during biting and therefore play a key role in pathogen transmission. Methods and Results: An original proteomic approach combining 2D-electrophoresis and mass spectrometry was adopted to compare the salivary expression profiles of two strains of C. quinquefasciatus with the same genetic background but carrying either the ace-1(R) resistance allele or not (wild type). Four salivary proteins were differentially expressed (> 2 fold, P < 0.05) in susceptible (SLAB) and resistant (SR) mosquito strains. Protein identification indicated that the D7 long form, a major salivary protein involved in blood feeding success, presented lower expression in the resistant strain than the susceptible strain. In contrast, three other proteins, including metabolic enzymes (endoplasmin, triosephosphate isomerase) were significantly over-expressed in the salivary gland of ace-1(R) resistant mosquitoes. A catalogue of 67 salivary proteins of C. quinquefasciatus sialotranscriptome was also identified and described. Conclusion: The "resistance"-dependent expression of salivary proteins in mosquitoes may have considerable impact on biting behaviour and hence on the capacity to transmit parasites/viruses to humans. The behaviour of susceptible and insecticide-resistant mosquitoes in the presence of vertebrate hosts and its impact on pathogen transmission urgently requires further investigation

Fast Mapping of Short Sequences with Mismatches, Insertions and Deletions Using Index Structures

With few exceptions, current methods for short read mapping make use of simple seed heuristics to speed up the search. Most of the underlying matching models neglect the necessity to allow not only mismatches, but also insertions and deletions. Current evaluations indicate, however, that very different error models apply to the novel high-throughput sequencing methods. While the most frequent error-type in Illumina reads are mismatches, reads produced by 454's GS FLX predominantly contain insertions and deletions (indels). Even though 454 sequencers are able to produce longer reads, the method is frequently applied to small RNA (miRNA and siRNA) sequencing. Fast and accurate matching in particular of short reads with diverse errors is therefore a pressing practical problem. We introduce a matching model for short reads that can, besides mismatches, also cope with indels. It addresses different error models. For example, it can handle the problem of leading and trailing contaminations caused by primers and poly-A tails in transcriptomics or the length-dependent increase of error rates. In these contexts, it thus simplifies the tedious and error-prone trimming step. For efficient searches, our method utilizes index structures in the form of enhanced suffix arrays. In a comparison with current methods for short read mapping, the presented approach shows significantly increased performance not only for 454 reads, but also for Illumina reads. Our approach is implemented in the software segemehl available at http://www.bioinf.uni-leipzig.de/Software/segemehl/

Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing

<p>Abstract</p> <p>Background</p> <p>The rapid evolution of 454 GS-FLX sequencing technology has not been accompanied by a reassessment of the quality and accuracy of the sequences obtained. Current strategies for decision-making and error-correction are based on an initial analysis by Huse <it>et al. </it>in 2007, for the older GS20 system based on experimental sequences. We analyze here the quality of 454 sequencing data and identify factors playing a role in sequencing error, through the use of an extensive dataset for Roche control DNA fragments.</p> <p>Results</p> <p>We obtained a mean error rate for 454 sequences of 1.07%. More importantly, the error rate is not randomly distributed; it occasionally rose to more than 50% in certain positions, and its distribution was linked to several experimental variables. The main factors related to error are the presence of homopolymers, position in the sequence, size of the sequence and spatial localization in PT plates for insertion and deletion errors. These factors can be described by considering seven variables. No single variable can account for the error rate distribution, but most of the variation is explained by the combination of all seven variables.</p> <p>Conclusions</p> <p>The pattern identified here calls for the use of internal controls and error-correcting base callers, to correct for errors, when available (e.g. when sequencing amplicons). For shotgun libraries, the use of both sequencing primers and deep coverage, combined with the use of random sequencing primer sites should partly compensate for even high error rates, although it may prove more difficult than previous thought to distinguish between low-frequency alleles and errors.</p

Impact of Chromatin Structures on DNA Processing for Genomic Analyses

Chromatin has an impact on recombination, repair, replication, and evolution of DNA. Here we report that chromatin structure also affects laboratory DNA manipulation in ways that distort the results of chromatin immunoprecipitation (ChIP) experiments. We initially discovered this effect at the Saccharomyces cerevisiae HMR locus, where we found that silenced chromatin was refractory to shearing, relative to euchromatin. Using input samples from ChIP-Seq studies, we detected a similar bias throughout the heterochromatic portions of the yeast genome. We also observed significant chromatin-related effects at telomeres, protein binding sites, and genes, reflected in the variation of input-Seq coverage. Experimental tests of candidate regions showed that chromatin influenced shearing at some loci, and that chromatin could also lead to enriched or depleted DNA levels in prepared samples, independently of shearing effects. Our results suggested that assays relying on immunoprecipitation of chromatin will be biased by intrinsic differences between regions packaged into different chromatin structures - biases which have been largely ignored to date. These results established the pervasiveness of this bias genome-wide, and suggested that this bias can be used to detect differences in chromatin structures across the genome

Whole Genome Sequencing Highlights Genetic Changes Associated with Laboratory Domestication of C. elegans

Defining the mutational landscape when individuals of a species grow separately and diverge over many generations can provide insights into trait evolution. A specific example of this involves studying changes associated with domestication where different lines of the same wild stock have been cultivated independently in different standard environments. Whole genome sequence comparison of such lines permits estimation of mutation rates, inference of genes' ancestral states and ancestry of existing strains, and correction of sequencing errors in genome databases. Here we study domestication of the C. elegans Bristol strain as a model, and report the genome sequence of LSJ1 (Bristol), a sibling of the standard C. elegans reference wild type N2 (Bristol). The LSJ1 and N2 lines were cultivated separately from shortly after the Bristol strain was isolated until methods to freeze C. elegans were developed. We find that during this time the two strains have accumulated 1208 genetic differences. We describe phenotypic variation between N2 and LSJ1 in the rate at which embryos develop, the rate of production of eggs, the maturity of eggs at laying, and feeding behavior, all the result of post-isolation changes. We infer the ancestral alleles in the original Bristol isolate and highlight 2038 likely sequencing errors in the original N2 reference genome sequence. Many of these changes modify genome annotation. Our study provides a starting point to further investigate genotype-phenotype association and offers insights into the process of selection as a result of laboratory domestication

Comparing De Novo Genome Assembly: The Long and Short of It

Recent advances in DNA sequencing technology and their focal role in Genome Wide Association Studies (GWAS) have rekindled a growing interest in the whole-genome sequence assembly (WGSA) problem, thereby, inundating the field with a plethora of new formalizations, algorithms, heuristics and implementations. And yet, scant attention has been paid to comparative assessments of these assemblers' quality and accuracy. No commonly accepted and standardized method for comparison exists yet. Even worse, widely used metrics to compare the assembled sequences emphasize only size, poorly capturing the contig quality and accuracy. This paper addresses these concerns: it highlights common anomalies in assembly accuracy through a rigorous study of several assemblers, compared under both standard metrics (N50, coverage, contig sizes, etc.) as well as a more comprehensive metric (Feature-Response Curves, FRC) that is introduced here; FRC transparently captures the trade-offs between contigs' quality against their sizes. For this purpose, most of the publicly available major sequence assemblers – both for low-coverage long (Sanger) and high-coverage short (Illumina) reads technologies – are compared. These assemblers are applied to microbial (Escherichia coli, Brucella, Wolbachia, Staphylococcus, Helicobacter) and partial human genome sequences (Chr. Y), using sequence reads of various read-lengths, coverages, accuracies, and with and without mate-pairs. It is hoped that, based on these evaluations, computational biologists will identify innovative sequence assembly paradigms, bioinformaticists will determine promising approaches for developing “next-generation” assemblers, and biotechnologists will formulate more meaningful design desiderata for sequencing technology platforms. A new software tool for computing the FRC metric has been developed and is available through the AMOS open-source consortium

RNA-Seq reveals large quantitative differences between the transcriptomes of outbreak and non-outbreak locusts

Outbreaks of locust populations repeatedly devastate economies and ecosystems in large parts of the world. The consequent behavioural shift from solitarious to gregarious and the concomitant changes in the locusts’ biology are of relevant scientific interest. Yet, research on the main locust species has not benefitted from recent advances in genomics. In this first RNA-Seq study on Schistocerca gregaria, we report two transcriptomes, including many novel genes, as well as differential gene expression results. In line with the large biological differences between solitarious and gregarious locusts, almost half of the transcripts are differentially expressed between their central nervous systems. Most of these transcripts are over-expressed in the gregarious locusts, suggesting positive correlations between the levels of activity at the population, individual, tissue and gene expression levels. We group these differentially expressed transcripts by gene function and highlight those that are most likely to be associated with locusts’ phase change either in a species-specific or general manner. Finally, we discuss our findings in the context of population-level and physiological events leading to gregariousness.M. Bakkali wishes to thank the Spanish Ministerio de Ciencia y Tecnología for the for the Ramón y Cajal fellowship and for the BFU2010-16438 grant that supported both this research and the FPI studentship to Rubén Martín Blázquez. We thank Mrs. Pernille Lavgesen for revision of the English language writing of this manuscript. We also thank the editor for the valuable comments on the manuscript

Integrating Communities of Practice in Technology Development Projects

Technology development projects usually benefit when knowledge and expertise are drawn from a variety of sources, including potential users. Orchestrating the involvement of people from disparate groups is a crucial task for project managers. It requires finding a balance between differentiation, when teams work in isolation, and integration, when groups come together to exchange knowledge. This article argues that a “community of practice” perspective can help project managers to achieve this balance, by drawing attention to the assumptions, interests, skills, and formal and tacit knowledge of the different groups involved. Successful integration can be achieved by ensuring that the developing technology is comprehensible to all the groups concerned, and making sure that it satisfies their various interests