The impact of read length on quantification of differentially expressed genes and splice junction detection

BackgroundThe initial next-generation sequencing technologies produced reads of 25 or 36 bp, and only from a single-end of the library sequence. Currently, it is possible to reliably produce 300 bp paired-end sequences for RNA expression analysis. While read lengths have consistently increased, people have assumed that longer reads are more informative and that paired-end reads produce better results than single-end reads. We used paired-end 101 bp reads and trimmed them to simulate different read lengths, and also separated the pairs to produce single-end reads. For each read length and paired status, we evaluated differential expression levels between two standard samples and compared the results to those obtained by qPCR.ResultsWe found that, with the exception of 25 bp reads, there is little difference for the detection of differential expression regardless of the read length. Once single-end reads are at a length of 50 bp, the results do not change substantially for any level up to, and including, 100 bp paired-end. However, splice junction detection significantly improves as the read length increases with 100 bp paired-end showing the best performance. We performed the same analysis on two ENCODE samples and found consistent results confirming that our conclusions have broad application.ConclusionsA researcher could save substantial resources by using 50 bp single-end reads for differential expression analysis instead of using longer reads. However, splicing detection is unquestionably improved by paired-end and longer reads. Therefore, an appropriate read length should be used based on the final goal of the study.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0697-y) contains supplementary material, which is available to authorized users.     

Specific detection of DNA and RNA targets using a novel isothermal nucleic acid amplification assay based on the formation of a three-way junction structure

The formation of DNA three-way junction (3WJ) structures has been utilised to develop a novel isothermal nucleic acid amplification assay (SMART) for the detection of specific DNA or RNA targets. The assay consists of two oligonucleotide probes that hybridise to a specific target sequence and, only then, to each other forming a 3WJ structure. One probe (template for the RNA signal) contains a non-functional single-stranded T7 RNA polymerase promoter sequence. This promoter sequence is made double-stranded (hence functional) by DNA polymerase, allowing T7 RNA polymerase to generate a target-dependent RNA signal which is measured by an enzyme-linked oligosorbent assay (ELOSA). The sequence of the RNA signal is always the same, regardless of the original target sequence. The SMART assay was successfully tested in model systems with several single-stranded synthetic targets, both DNA and RNA. The assay could also detect specific target sequences in both genomic DNA and total RNA from Escherichia coli. It was also possible to generate signal from E.coli samples without prior extraction of nucleic acid, showing that for some targets, sample purification may not be required. The assay is simple to perform and easily adaptable to different targets.     

MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform

A multiple sequence alignment program, MAFFT, has been developed. The CPU time is drastically reduced as compared with existing methods. MAFFT includes two novel techniques. (i) Homo logous regions are rapidly identified by the fast Fourier transform (FFT), in which an amino acid sequence is converted to a sequence composed of volume and polarity values of each amino acid residue. (ii) We propose a simplified scoring system that performs well for reducing CPU time and increasing the accuracy of alignments even for sequences having large insertions or extensions as well as distantly related sequences of similar length. Two different heuristics, the progressive method (FFT-NS-2) and the iterative refinement method (FFT-NS-i), are implemented in MAFFT. The performances of FFT-NS-2 and FFT-NS-i were compared with other methods by computer simulations and benchmark tests; the CPU time of FFT-NS-2 is drastically reduced as compared with CLUSTALW with comparable accuracy. FFT-NS-i is over 100 times faster than T-COFFEE, when the number of input sequences exceeds 60, without sacrificing the accuracy.     

Sequence-specific detection using a universal probe system based on the formation of a four-way junction structure

We have developed a novel hybridization detection system using a universal probe based on the formation of a four-way junction (4WJ) structure. This methodology employs a combination of two sequence-specific probes and a universal quenching probe, and the same universal probe can be used for any target gene, allowing cost-effective assays. This 4WJ detection is ideal for extensive parallel identification of nucleic acids such as in multiplex polymerase chain reaction (PCR) systems. Compared with gel electrophoresis, this detection procedure is not only sensitive and rapid but also free of hazardous chemicals such as ethidium bromide. In addition, the 4WJ hybridization technology is more specific as an identifier than the size of a band on an agarose gel. We used a model multiplex PCR method that detected eight different virulence genes in Escherichia coli isolates, demonstrating that our 4WJ detection system is rapid, sensitive, and specific.     

A novel HPLC procedure for detection and quantification of aminoacetone, a precursor of methylglyoxal, in biological samples

Increase in methylglyoxal is thought to be involved in different pathological conditions. Deamination of aminoacetone by semicarbazide-sensitive amine oxidase (SSAO) leads to production of methylglyoxal. We have synthesized aminoacetone and developed a novel HPLC procedure for its quantitative determination. The urinary excretion of aminoacetone is approximately 20-30 microg/mouse/day, and the concentration is about 0.5 microg/g in mouse liver and small intestine. SSAO inhibitor increases aminoacetone levels in both tissues and urines. Results confirm that aminoacetone is an endogenous substrate for SSAO. However, data also indicate that deamination is not the only catabolic pathway for aminoacetone.     

TM-align: a protein structure alignment algorithm based on the TM-score

We have developed TM-align, a new algorithm to identify the best structural alignment between protein pairs that combines the TM-score rotation matrix and Dynamic Programming (DP). The algorithm is approximately 4 times faster than CE and 20 times faster than DALI and SAL. On average, the resulting structure alignments have higher accuracy and coverage than those provided by these most often-used methods. TM-align is applied to an all-against-all structure comparison of 10 515 representative protein chains from the Protein Data Bank (PDB) with a sequence identity cutoff <95%: 1996 distinct folds are found when a TM-score threshold of 0.5 is used. We also use TM-align to match the models predicted by TASSER for solved non-homologous proteins in PDB. For both folded and misfolded models, TM-align can almost always find close structural analogs, with an average root mean square deviation, RMSD, of 3 A and 87% alignment coverage. Nevertheless, there exists a significant correlation between the correctness of the predicted structure and the structural similarity of the model to the other proteins in the PDB. This correlation could be used to assist in model selection in blind protein structure predictions. The TM-align program is freely downloadable at http://bioinformatics.buffalo.edu/TM-align.     

Conformational impurity of disulfide proteins: detection, quantification, and properties

The conformations of native proteins are in principle, and in most cases, dictated by the law of thermodynamics. Accordingly, a native protein must always exist in equilibrium with a minor concentration of nonnative (denatured) conformational isomers even at nondenaturing conditions. The presence of an infinitesimal quantity of nonnative conformational isomers at physiological conditions is biologically relevant due to their propensity to aggregate, which is an underlying cause of many neurodegenerative diseases. However, their detection and quantification are inherently difficult. In this article, we describe a simple strategy using the technique of disulfide scrambling to identify and quantify such minute concentrations of nonnative isomers. It is demonstrated that even for small stable proteins such as epidermal growth factor and hirudin, approximately 1% of heterogeneous nonnative isomers coexist with the native proteins under physiological conditions.     

A novel fluorescence-based assay for the rapid detection and quantification of cellular deoxyribonucleoside triphosphates

Current methods for measuring deoxyribonucleoside triphosphates (dNTPs) employ reagent and labor-intensive assays utilizing radioisotopes in DNA polymerase-based assays and/or chromatography-based approaches. We have developed a rapid and sensitive 96-well fluorescence-based assay to quantify cellular dNTPs utilizing a standard real-time PCR thermocycler. This assay relies on the principle that incorporation of a limiting dNTP is required for primer-extension and Taq polymerase-mediated 5–3′ exonuclease hydrolysis of a dual-quenched fluorophore-labeled probe resulting in fluorescence. The concentration of limiting dNTP is directly proportional to the fluorescence generated. The assay demonstrated excellent linearity (R² > 0.99) and can be modified to detect between ∼0.5 and 100 pmol of dNTP. The limits of detection (LOD) and quantification (LOQ) for all dNTPs were defined as <0.77 and <1.3 pmol, respectively. The intra-assay and inter-assay variation coefficients were determined to be <4.6% and <10%, respectively with an accuracy of 100 ± 15% for all dNTPs. The assay quantified intracellular dNTPs with similar results obtained from a validated LC–MS/MS approach and successfully measured quantitative differences in dNTP pools in human cancer cells treated with inhibitors of thymidylate metabolism. This assay has important application in research that investigates the influence of pathological conditions or pharmacological agents on dNTP biosynthesis and regulation.     

The effects of alignment error and alignment filtering on the sitewise detection of positive selection

When detecting positive selection in proteins, the prevalence of errors resulting from misalignment and the ability of alignment filters to mitigate such errors are not well understood, but filters are commonly applied to try to avoid false positive results. Focusing on the sitewise detection of positive selection across a wide range of divergence levels and indel rates, we performed simulation experiments to quantify the false positives and false negatives introduced by alignment error and the ability of alignment filters to improve performance. We found that some aligners led to many false positives, whereas others resulted in very few. False negatives were a problem for all aligners, increasing with sequence divergence. Of the aligners tested, PRANK's codon-based alignments consistently performed the best and ClustalW performed the worst. Of the filters tested, GUIDANCE performed the best and Gblocks performed the worst. Although some filters showed good ability to reduce the error rates from ClustalW and MAFFT alignments, none were found to substantially improve the performance of PRANK alignments under most conditions. Our results revealed distinct trends in error rates and power levels for aligners and filters within a biologically plausible parameter space. With the best aligner, a low false positive rate was maintained even with extremely divergent indel-prone sequences. Controls using the true alignment and an optimal filtering method suggested that performance improvements could be gained by improving aligners or filters to reduce the prevalence of false negatives, especially at higher divergence levels and indel rates.     

KeR-EGI,a new index of gait quantification based on electromyography

PurposeTo define a new index of gait pathology in adults based on electromyographic data: the Ker-EGI for Kerpape-Rennes EMG-based Gait Index. The principle is similar to the one of Gait Deviation Index but using EMG profiles instead of joint angles. It first needs to build a database of healthy subjects gait to be able then to quantify the deviation of one peculiar patient’s gait from this typical behavior.MethodsNinety adults (59 healthy and 31 pathological) participated to this study. All pathological subjects had a diagnosis of central nervous system disorder. On each subject we collected the joint angles and the activation profile of seven muscles of each lower limb. Moreover, we recorded two videos (face and profile) of each patient to compute his/her Edinburgh Visual Gait Score (EVGS). Then for each patient, we computed the GGI (Gillette Gait Index), the GDI (Gait Deviation Index) and the Ker-EGI.ResultsCorrelation Ker-EGI and each of the three kinematical indices (GGI, GDI, EVGS) is fair to good (respectively R² = 0.62, 0.42, and 0.69).ConclusionKeR-EGI is a valid index to evaluate gait and is complementary to one of these three kinematical indices providing synthetic vision on patients’ motor control abilities.     

Label-free detection of DNA hybridization based on a novel functionalized conducting polymer

A new functionalized pyrrole monomer, 3-pyrrolylacrylic acid (PAA) was synthesized. It was used to prepare a copolymer with pyrrole, poly(Py-co-PAA), which was investigated by reflective FT-IR, UV-vis spectroscopy and cyclic voltammetry. A label-free DNA sensor was prepared based on a poly(Py-co-PAA) film. Hybridization with complementary and non-complementary DNA targets was studied by electrochemical impedance spectroscopy. Results show a significant increase in the charge-transfer resistance upon addition of complementary target. The impedance spectra were analyzed by using a modified Randles and Ershler equivalent circuit model. The change in charge-transfer resistance that was used as an index of sensor response was found to be linear with logarithmic target concentration in the range of 2 x 10(-9) to 2 x 10(-7)M. The detection limit was 0.98 nM.     

CASA: a server for the critical assessment of protein sequence alignment accuracy

SUMMARY: A public server for evaluating the accuracy of protein sequence alignment methods is presented. CASA is an implementation of the alignment accuracy benchmark presented by Sauder et al. (Proteins, 40, 6-22, 2000). The benchmark currently contains 39321 pairwise protein structure alignments produced with the CE program from SCOP domain definitions. The server produces graphical and tabular comparisons of the accuracy of a user's input sequence alignments with other commonly used programs, such as BLAST, PSI-BLAST, Clustal W, and SAM-T99. AVAILABILITY: The server is located at http://capb.dbi.udel.edu/casa.     

Chemical sensing by a novel electrical oscillator: detection and quantitation of polysaccharides in concanavalin A solutions

Electrical oscillations across two platinum electrodes connected to an external circuit and immersed in HEPES buffer solutions of concanavalin A (Con A) were measured at various concentrations of mannan, starch and dextran. The frequency of oscillations was found to change almost linearly with the logarithm of the concentration of mannan between 10(-4) and 1 w/v%; whereas the frequency remained nearly constant with dextran. With starch it slightly increased as the concentration increased. This method was suggested to be useful for quantitative and selective measurement of antigen-antibody reaction and also for detection of various cells with the specific binding site such as cancer cells.     

A point mutation in the conserved hexanucleotide at a yeast 5' splice junction uncouples recognition, cleavage, and ligation

We have constructed an actin-HIS4 gene fusion, such that expression of HIS4 requires proper splicing of the actin intron. Using this chimeric gene in an in vivo screen for splicing mutations, we have isolated a G to A transition in the fifth position of the yeast 5' consensus sequence/GTAPyGT. This mutation still allows the junction to be recognized by the splicing machinery, albeit inefficiently. Surprisingly, the fidelity of the 5' endonucleolytic cleavage is also reduced. This results in an incorrect cleavage 6 nucleotides 5' of the 5' junction, at the dinucleotide/AT. Cleavage at this abnormal site does not lead to the production of mature mRNA, although this species appears to be in a lariat structure. The behavior of this mutant argues that recognition of the 5' junction and subsequent cleavage are separable events and, furthermore, that requirements for 3' endonucleolytic cleavage may be more complex than previously imagined.