Ebbie: automated analysis and storage of small RNA cloning data using a dynamic web server

BACKGROUND: DNA sequencing is used ubiquitously: from deciphering genomes[1] to determining the primary sequence of small RNAs (smRNAs) [2-5]. The cloning of smRNAs is currently the most conventional method to determine the actual sequence of these important regulators of gene expression. Typical smRNA cloning projects involve the sequencing of hundreds to thousands of smRNA clones that are delimited at their 5' and 3' ends by fixed sequence regions. These primers result from the biochemical protocol used to isolate and convert the smRNA into clonable PCR products. Recently we completed a smRNA cloning project involving tobacco plants, where analysis was required for ~700 smRNA sequences[6]. Finding no easily accessible research tool to enter and analyze smRNA sequences we developed Ebbie to assist us with our study. RESULTS: Ebbie is a semi-automated smRNA cloning data processing algorithm, which initially searches for any substring within a DNA sequencing text file, which is flanked by two constant strings. The substring, also termed smRNA or insert, is stored in a MySQL and BlastN database. These inserts are then compared using BlastN to locally installed databases allowing the rapid comparison of the insert to both the growing smRNA database and to other static sequence databases. Our laboratory used Ebbie to analyze scores of DNA sequencing data originating from an smRNA cloning project[6]. Through its built-in instant analysis of all inserts using BlastN, we were able to quickly identify 33 groups of smRNAs from ~700 database entries. This clustering allowed the easy identification of novel and highly expressed clusters of smRNAs. Ebbie is available under GNU GPL and currently implemented on CONCLUSION: Ebbie was designed for medium sized smRNA cloning projects with about 1,000 database entries [6-8].Ebbie can be used for any type of sequence analysis where two constant primer regions flank a sequence of interest. The reliable storage of inserts, and their annotation in a MySQL database, BlastN[9] comparison of new inserts to dynamic and static databases make it a powerful new tool in any laboratory using DNA sequencing. Ebbie also prevents manual mistakes during the excision process and speeds up annotation and data-entry. Once the server is installed locally, its access can be restricted to protect sensitive new DNA sequencing data. Ebbie was primarily designed for smRNA cloning projects, but can be applied to a variety of RNA and DNA cloning projects[2,3,10,11]

Differences in the bovine milk whey proteome between early pregnancy and the estrous cycle

peer-reviewedCurrent bovine pregnancy detection methods are not reliable until at least day 28 post artificial insemination (AI). The bovine estrous cycle is approximately 21 days; consequently, producers miss an opportunity to rebreed at the next estrous event. Therefore, commercial interest exists for the discovery of novel biomarkers of pregnancy which could reliably detect pregnancy status at or before day 21 of pregnancy. The objective of the present study was to use liquid chromatography tandem mass spectrometry (LC-MS/MS) to perform a global, label-free, proteomics study on (i) milk whey and (ii) extracellular vesicle (EV) enriched milk whey samples, from day 21 of pregnancy, compared with day 21 of the estrous cycle, in order to identify potential protein biomarkers of early pregnancy. The estrous cycles of 10 dairy cows were synchronized, they went through one (control) estrous cycle and these cows were artificially inseminated during the following estrus. These cows were confirmed pregnant by ultrasound scanning. Milk whey samples were collected on day 21 of the estrous cycle and on day 21 post AI. Milk whey samples and EV enriched milk whey samples were analyzed by LC-MS/MS and subsequent analyzes of the label-free quantitative data was performed in MaxQuant and Perseus. Four proteins (APOB, SPADH1, PLIN2 and LPO) were differentially expressed between the proteomes of milk whey from day 21 of pregnancy and day 21 of the estrous cycle (P < 0.05). Ten proteins (PIGR, PGD, QSOX1, MUC1, SRPRA, MD2, GAPDH, FOLR1, GPRC5B and HHIPL2) were differentially expressed between the proteomes of EV enriched milk whey from day 21 of pregnancy and day 21 of the estrous cycle (P < 0.05). These proteins are potential milk whey biomarkers of early pregnancy

Meta-analysis of small RNA-sequencing errors reveals ubiquitous post-transcriptional RNA modifications

Recent advances in DNA-sequencing technology have made it possible to obtain large datasets of small RNA sequences. Here we demonstrate that not all non-perfectly matched small RNA sequences are simple technological sequencing errors, but many hold valuable biological information. Analysis of three small RNA datasets originating from Oryza sativa and Arabidopsis thaliana small RNA-sequencing projects demonstrates that many single nucleotide substitution errors overlap when aligning homologous non-identical small RNA sequences. Investigating the sites and identities of substitution errors reveal that many potentially originate as a result of post-transcriptional modifications or RNA editing. Modifications include N1-methyl modified purine nucleotides in tRNA, potential deamination or base substitutions in micro RNAs, 3′ micro RNA uridine extensions and 5′ micro RNA deletions. Additionally, further analysis of large sequencing datasets reveal that the combined effects of 5′ deletions and 3′ uridine extensions can alter the specificity by which micro RNAs associate with different Argonaute proteins. Hence, we demonstrate that not all sequencing errors in small RNA datasets are technical artifacts, but that these actually often reveal valuable biological insights to the sites of post-transcriptional RNA modifications

Platelet-Derived Microparticles From Obese Individuals: Characterization of Number, Size, Proteomics, and Crosstalk With Cancer and Endothelial Cells

Rationale: Obesity is a risk factor for atherothrombosis and various cancers. However, the mechanisms are not yet completely clarified.Objectives: We aimed to verify whether the microparticles (MPs) released from thrombin-activated platelets differed in obese and non-obese women for number, size, and proteomics cargo and the capacity to modulate in vitro the expression of (i) genes related to the epithelial to mesenchymal transition (EMT) and the endothelial to mesenchymal transition (EndMT), and (ii) cyclooxygenase (COX)-2 involved in the production of angiogenic and inflammatory mediators.Methods and Results: MPs were obtained from thrombin activated platelets of four obese and their matched non-obese women. MPs were analyzed by cytofluorimeter and protein content by liquid chromatography-mass spectrometry. MPs from obese women were not different in number but showed increased heterogeneity in size. In obese individuals, MPs containing mitochondria (mitoMPs) expressed lower CD41 levels and increased phosphatidylserine associated with enhanced Factor V representing a signature of a prothrombotic state. Proteomics analysis identified 44 proteins downregulated and three upregulated in MPs obtained from obese vs. non-obese women. A reduction in the proteins of the α-granular membrane and those involved in mitophagy and antioxidant defenses-granular membrane was detected in the MPs of obese individuals. MPs released from platelets of obese individuals were more prone to induce the expression of marker genes of EMT and EndMT when incubated with human colorectal cancer cells (HT29) and human cardiac microvascular endothelial cells (HCMEC), respectively. A protein, highly enhanced in obese MPs, was the pro-platelet basic protein with pro-inflammatory and tumorigenic actions. Exclusively MPs from obese women induced COX-2 in HCMEC.Conclusion: Platelet-derived MPs of obese women showed higher heterogeneity in size and contained different levels of proteins relevant to thrombosis and tumorigenesis. MPs from obese individuals presented enhanced capacity to cause changes in the expression of EMT and EndMT marker genes and to induce COX-2. These effects might contribute to the increased risk for the development of thrombosis and multiple malignancies in obesity.Clinical Trial Registration:www.ClinicalTrials.gov, identifier NCT01581801

Searching for a Stochastic Background of Gravitational Waves with LIGO

The Laser Interferometer Gravitational-wave Observatory (LIGO) has performed the fourth science run, S4, with significantly improved interferometer sensitivities with respect to previous runs. Using data acquired during this science run, we place a limit on the amplitude of a stochastic background of gravitational waves. For a frequency independent spectrum, the new limit is $\Omega_{\rm GW} < 6.5 \times 10^{-5}$. This is currently the most sensitive result in the frequency range 51-150 Hz, with a factor of 13 improvement over the previous LIGO result. We discuss complementarity of the new result with other constraints on a stochastic background of gravitational waves, and we investigate implications of the new result for different models of this background.Comment: 37 pages, 16 figure

Upper limit map of a background of gravitational waves

We searched for an anisotropic background of gravitational waves using data from the LIGO S4 science run and a method that is optimized for point sources. This is appropriate if, for example, the gravitational wave background is dominated by a small number of distinct astrophysical sources. No signal was seen. Upper limit maps were produced assuming two different power laws for the source strain power spectrum. For an f^-3 power law and using the 50 Hz to 1.8 kHz band the upper limits on the source strain power spectrum vary between 1.2e-48 Hz^-1 (100 Hz/f)^3 and 1.2e-47 Hz^-1 (100 Hz /f)^3, depending on the position in the sky. Similarly, in the case of constant strain power spectrum, the upper limits vary between 8.5e-49 Hz^-1 and 6.1e-48 Hz^-1. As a side product a limit on an isotropic background of gravitational waves was also obtained. All limits are at the 90% confidence level. Finally, as an application, we focused on the direction of Sco-X1, the closest low-mass X-ray binary. We compare the upper limit on strain amplitude obtained by this method to expectations based on the X-ray luminosity of Sco-X1.Comment: 11 pages, 9 figures, 2 table

Upper limit map of a background of gravitational waves

We searched for an anisotropic background of gravitational waves using data from the LIGO S4 science run and a method that is optimized for point sources. This is appropriate if, for example, the gravitational wave background is dominated by a small number of distinct astrophysical sources. No signal was seen. Upper limit maps were produced assuming two different power laws for the source strain power spectrum. For an f^-3 power law and using the 50 Hz to 1.8 kHz band the upper limits on the source strain power spectrum vary between 1.2e-48 Hz^-1 (100 Hz/f)^3 and 1.2e-47 Hz^-1 (100 Hz /f)^3, depending on the position in the sky. Similarly, in the case of constant strain power spectrum, the upper limits vary between 8.5e-49 Hz^-1 and 6.1e-48 Hz^-1. As a side product a limit on an isotropic background of gravitational waves was also obtained. All limits are at the 90% confidence level. Finally, as an application, we focused on the direction of Sco-X1, the closest low-mass X-ray binary. We compare the upper limit on strain amplitude obtained by this method to expectations based on the X-ray luminosity of Sco-X1.Comment: 11 pages, 9 figures, 2 table

Structure-based assessment and network analysis of targeting 14-3-3 proteins in prostate cancer

Abstract Developing combination therapy for castrate-resistant prostate cancer (CRPC) may require exploiting new drug targets outside androgen receptor and PI3K / AKT / mTOR signal transduction pathways implicated in prostate cancer (PCa) progression. One such possible new target is YWHAZ of the 14-3-3 protein family as this gene has prognostic significance for metastatic CRPC patients. However, there are no small molecules targeting YWHAZ commercially available. Hence, we explored whether the small molecule BV02 targeting another 14-3-3 protein family member SFN also binds to YWHAZ. Using advanced docking algorithms we find that BV02 docks many other 14-3-3 family members. In addition, the amphipathic groove where drug binding occurs also has a high binding affinity for other drugs used to treat PCa such as docetaxel. The proteome of metastatic PCa models (LNCaP clone FGC and PC-3) was perturbed as a result of BV02 treatment. Through data integration of three proteomics data sets we found that BV02 modulates numerous protein-protein interactions involving 14-3-3 proteins in our PCa models

Characterizing multiple exogenous and endogenous small RNA populations in parallel with subfemtomolar sensitivity using a streptavidin gel-shift assay

Here we present a simple and inexpensive gel-shift assay for the detection and quantification of small RNAs. The assay is at least 5–10 times more sensitive than a conventional Northern, and is highly scalable. Total RNA is first size purified to enrich the desired size range, phosphatase treated, and then radiolabeled to high specific activity using polynucleotide kinase. The resulting RNA stock is then hybridized to an excess of biotinylated DNA probe oligonucleotide, prior to mixing with streptavidin and loading on a native gel. The amount of supershifted material was proportional to the amount of labeled target RNA in the sample. We applied this method to verify sequencing data originally obtained from a four-point comparison study on the effect of endogenous expression of HC-Pro on Y-satellite/cucumber mosaic virus infection in tobacco plants. The results of the streptavidin gel-shift assay were consistent with the concentrations of small RNA infected plants inferred by our original cloning data, and rapidly provided information about the relative concentration of a number of viral and endogenous small RNAs. Further straightforward improvements to this simple methodology might be expected to improve the methods sensitivity by as much as another 10-fold