BackgroundAdapter trimming is a prerequisite step for analyzing next-generation sequencing (NGS) data when the reads are longer than the target DNA/RNA fragments. Although typically used in small RNA sequencing, adapter trimming is also used widely in other applications, such as genome DNA sequencing and transcriptome RNA/cDNA sequencing, where fragments shorter than a read are sometimes obtained because of the limitations of NGS protocols. For the newly emerged Nextera long mate-pair (LMP) protocol, junction adapters are located in the middle of all properly constructed fragments; hence, adapter trimming is essential to gain the correct paired reads. However, our investigations have shown that few adapter trimming tools meet both efficiency and accuracy requirements simultaneously. The performances of these tools can be even worse for paired-end and/or mate-pair sequencing.ResultsTo improve the efficiency of adapter trimming, we devised a novel algorithm, the bit-masked k-difference matching algorithm, which has O(kn) expected time with O(m) space, where k is the maximum number of differences allowed, n is the read length, and m is the adapter length. This algorithm makes it possible to fully enumerate all candidates that meet a specified threshold, e.g. error ratio, within a short period of time. To improve the accuracy of this algorithm, we designed a simple and easy-to-explain statistical scoring scheme to evaluate candidates in the pattern matching step. We also devised scoring schemes to fully exploit the paired-end/mate-pair information when it is applicable. All these features have been implemented in an industry-standard tool named Skewer (https://sourceforge.net/projects/skewer). Experiments on simulated data, real data of small RNA sequencing, paired-end RNA sequencing, and Nextera LMP sequencing showed that Skewer outperforms all other similar tools that have the same utility. Further, Skewer is considerably faster than other tools that have comparative accuracies; namely, one times faster for single-end sequencing, more than 12 times faster for paired-end sequencing, and 49% faster for LMP sequencing.ConclusionsSkewer achieved as yet unmatched accuracies for adapter trimming with low time bound.

Dicer enzymes process virus-specific double-stranded RNA (dsRNA) into small interfering RNAs (siRNAs) to initiate specific antiviral defense by related RNA interference (RNAi) pathways in plants, insects, nematodes, and mammals. Antiviral RNAi in Caenorhabditis elegans requires Dicer-related helicase 1 (DRH-1), not found in plants and insects but highly homologous to mammalian retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs), intracellular viral RNA sensors that trigger innate immunity against RNA virus infection. However, it remains unclear if DRH-1 acts analogously to initiate antiviral RNAi in C. elegans. Here, we performed a forward genetic screen to characterize antiviral RNAi in C. elegans. Using a mapping-by-sequencing strategy, we uncovered four loss-of-function alleles of drh-1, three of which caused mutations in the helicase and C-terminal domains conserved in RLRs. Deep sequencing of small RNAs revealed an abundant population of Dicer-dependent virus-derived small interfering RNAs (vsiRNAs) in drh-1 single and double mutant animals after infection with Orsay virus, a positive-strand RNA virus. These findings provide further genetic evidence for the antiviral function of DRH-1 and illustrate that DRH-1 is not essential for the sensing and Dicer-mediated processing of the viral dsRNA replicative intermediates. Interestingly, vsiRNAs produced by drh-1 mutants were mapped overwhelmingly to the terminal regions of the viral genomic RNAs, in contrast to random distribution of vsiRNA hot spots when DRH-1 is functional. As RIG-I translocates on long dsRNA and DRH-1 exists in a complex with Dicer, we propose that DRH-1 facilitates the biogenesis of vsiRNAs in nematodes by catalyzing translocation of the Dicer complex on the viral long dsRNA precursors.

The genus Dacus is one of the most economically important tephritid fruit flies. The first complete mitochondrial genome (mitogenome) of Dacus species – D. longicornis was sequenced by next-generation sequencing in order to develop the mitogenome data for this genus. The circular 16,253 bp mitogenome is the typical set and arrangement of 37 genes present in the ancestral insect. The mitogenome data of D. longicornis was compared to all the published homologous sequences of other tephritid species. We discovered the subgenera Bactrocera, Daculus and Tetradacus differed from the subgenus Zeugodacus, the genera Dacus, Ceratitis and Procecidochares in the possession of TA instead of TAA stop codon for COI gene. There is a possibility that the TA stop codon in COI is the synapomorphy in Bactrocera group in the genus Bactrocera comparing with other Tephritidae species. Phylogenetic analyses based on the mitogenome data from Tephritidae were inferred by Bayesian and Maximum-likelihood methods, strongly supported the sister relationship between Zeugodacus and Dacus.

Leaf chlorosis induced by plant virus infection has a short fluorescence lifetime, which reflects damaged photosynthetic complexes and degraded chloroplasts. Plant viruses often induce chlorosis and necrosis, which are intimately related to photosynthetic functions. Chlorophyll fluorescence lifetime measurement is a valuable noninvasive tool for analyzing photosynthetic processes and is a sensitive indicator of the environment surrounding the fluorescent molecules. In this study, our central goal was to explore the effect of viral infection on photosynthesis by employing chlorophyll fluorescence lifetime imaging (FLIM), steady-state fluorescence, non-photochemical quenching (NPQ), transmission electron microscopy (TEM), and pigment analysis. The data indicated that the chlorophyll fluorescence lifetime of chlorotic leaves was significantly shorter than that of healthy control leaves, and the fitted short lifetime component of chlorophyll fluorescence of chlorotic leaves was dominant. This dominant short lifetime component may result from damage to the structure of thylakoid, which was confirmed by TEM. The NPQ value of chlorotic leaves was slightly higher than that of healthy green leaves, which can be explained by increased neoxanthin, lutein and violaxanthin content relative to chlorophyll a. The difference in NPQ is slight, but FLIM can provide simple and direct characterization of PSII structure and photosynthetic function. Therefore, this technique shows great potential as a simple and rapid method for studying mechanisms of plant virus infection.

Graphical abstract