Error, reproducibility and sensitivity : a pipeline for data processing of Agilent oligonucleotide expression arrays

Background Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples. Results We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2% of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log2 units ( 6% of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators. Conclusions This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells

The SNAPSHOT study protocol : SNAcking, Physical activity, Self-regulation, and Heart rate Over Time

Peer reviewedPublisher PD

Serum levels of the angiogenic factor pleiotrophin in relation to disease stage in lung cancer patients

Pleiotrophin is a heparin-binding growth factor involved in the differentiation and proliferation of neuronal tissue during embryogenesis, and also secreted by melanoma and breast carcinoma cells. Pleiotrophin exhibits mitogenic and angiogenic properties and has been shown to influence the vascular supply, expansion and metastasis of tumour cells. Our aim was to study the serum and plasma concentrations of pleiotrophin and the classical angiogenic growth factor vascular endothelial growth factor. Using a specific ELISA-test we studied patients with small cell lung cancer (n=63), and patients with non-small cell lung cancer (n=22) in comparison to healthy control subjects (n=41). In most of the lung cancer patients (81%), we found serum levels of pleiotrophin above those of control subjects (P<0.001). Of the 63 small cell lung cancer patients in the study pleiotrophin serum levels were elevated in 55 cases (87%) and in 14 cases (63%) of the 22 non-small cell lung cancer patients. Pleiotrophin mean serum concentrations were 10.8-fold higher in the tumour patient group as compared to the control group (P<0.001). Furthermore, pleiotrophin serum levels correlated positively with the stage of disease and inversely with the response to therapy. Plasma vascular endothelial growth factor concentrations were elevated in only in 28.6% of small cell lung cancer and 45.5% of non-small cell lung cancer patients by an average of 2.3-fold. Quite strikingly, there was no apparent correlation between the plasma vascular endothelial growth factor concentration and the stage of disease. Our study suggests that pleiotrophin may be an early indicator of lung cancer and might be of use in monitoring the efficacy of therapy, which needs to be confirmed by larger studies

Validation of Plasmodium falciparum dUTPase as the target of 5'-tritylated deoxyuridine analogues with anti-malarial activity

BACKGROUND: Malaria remains as a major global problem, being one of the infectious diseases that engender highest mortality across the world. Due to the appearance of resistance and the lack of an effective vaccine, the search of novel anti-malarials is required. Deoxyuridine 5'-triphosphate nucleotido-hydrolase (dUTPase) is responsible for the hydrolysis of dUTP to dUMP within the parasite and has been proposed as an essential step in pyrimidine metabolism by providing dUMP for thymidylate biosynthesis. In this work, efforts to validate dUTPase as a drug target in Plasmodium falciparum are reported. METHODS: To investigate the role of PfdUTPase in cell survival different strategies to generate knockout mutants were used. For validation of PfdUTPase as the intracellular target of four inhibitors of the enzyme, mutants overexpressing PfdUTPase and HsdUTPase were created and the IC50 for each cell line with each compound was determined. The effect of these compounds on dUTP and dTTP levels from P. falciparum was measured using a DNA polymerase assay. Detailed localization studies by indirect immunofluorescence microscopy and live cell imaging were also performed using a cell line overexpressing a Pfdut-GFP fusion protein. RESULTS:Different attempts of disruption of the dut gene of P. falciparum were unsuccessful while a 3' replacement construct could recombine correctly in the locus suggesting that the enzyme is essential. The four 5'-tritylated deoxyuridine analogues described are potent inhibitors of the P. falciparum dUTPase and exhibit antiplasmodial activity. Overexpression of the Plasmodium and human enzymes conferred resistance against selective compounds, providing chemical validation of the target and confirming that indeed dUTPase inhibition is involved in anti-malarial activity. In addition, incubation with these inhibitors was associated with a depletion of the dTTP pool corroborating the central role of dUTPase in dTTP synthesis. PfdUTPase is mainly localized in the cytosol. CONCLUSION: These results strongly confirm the pivotal and essential role of dUTPase in pyrimidine biosynthesis of P. falciparum intraerythrocytic stages

Anonymisation of geographical distance matrices via Lipschitz embedding

BACKGROUND: Anonymisation of spatially referenced data has received increasing attention in recent years. Whereas the research focus has been on the anonymisation of point locations, the disclosure risk arising from the publishing of inter-point distances and corresponding anonymisation methods have not been studied systematically. METHODS: We propose a new anonymisation method for the release of geographical distances between records of a microdata file-for example patients in a medical database. We discuss a data release scheme in which microdata without coordinates and an additional distance matrix between the corresponding rows of the microdata set are released. In contrast to most other approaches this method preserves small distances better than larger distances. The distances are modified by a variant of Lipschitz embedding. RESULTS: The effects of the embedding parameters on the risk of data disclosure are evaluated by linkage experiments using simulated data. The results indicate small disclosure risks for appropriate embedding parameters. CONCLUSION: The proposed method is useful if published distance information might be misused for the re-identification of records. The method can be used for publishing scientific-use-files and as an additional tool for record-linkage studies

Gene expression profile of peripheral blood lymphocytes from renal cell carcinoma patients treated with IL-2, Interferon-α and dendritic cell vaccine

© The Author(s), 2012. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS One 7 (2012): e50221, doi:10.1371/journal.pone.0050221.Lymphocytes are a key component of the immune system and their differentiation and function are directly influenced by cancer. We examined peripheral blood lymphocyte (PBL) gene expression as a biomarker of illness and treatment effect using the Affymetrix Human Gene ST1 platform in patients with metastatic renal cell carcinoma (mRCC) who received combined treatment with IL-2, interferon-?-2a and dendritic cell vaccine. We examined gene expression, cytokine levels in patient serum and lymphocyte subsets as determined by flow cytometry (FCM). Pre-treatment PBLs from patients with mRCC exhibit a gene expression profile and serum cytokine profile consistent with inflammation and proliferation not found in healthy donors (HD). PBL gene expression from patients with mRCC showed increased mRNA of genes involved with T-cell and TREG-cell activation pathways, which was also reflected in lymphocyte subset distribution. Overall, PBL gene expression post-treatment (POST) was not significantly different than pre-treatment (PRE). Nevertheless, treatment related changes in gene expression (post-treatment minus pre-treatment) revealed an increased expression of T-cell and B-cell receptor signaling pathways in responding (R) patients compared to non-responding (NR) patients. In addition, we observed down-regulation of TREG-cell pathways post-treatment in R vs. NR patients. While exploratory in nature, this study supports the hypothesis that enhanced inflammatory cytotoxic pathways coupled with blunting of the regulatory pathways is necessary for effective anti-cancer activity associated with immune therapy. This type of analysis can potentially identify additional immune therapeutic targets in patients with mRCC.This work was supported by grants from the National Institutes of Health (RO1 CA5648, R21CA112761, P20RR016437, and P30CA023108)

Path and Ridge Regression Analysis of Seed Yield and Seed Yield Components of Russian Wildrye (Psathyrostachys juncea Nevski) under Field Conditions

The correlations among seed yield components, and their direct and indirect effects on the seed yield (Z) of Russina wildrye (Psathyrostachys juncea Nevski) were investigated. The seed yield components: fertile tillers m-2 (Y1), spikelets per fertile tillers (Y2), florets per spikelet- (Y3), seed numbers per spikelet (Y4) and seed weight (Y5) were counted and the Z were determined in field experiments from 2003 to 2006 via big sample size. Y1 was the most important seed yield component describing the Z and Y2 was the least. The total direct effects of the Y1, Y3 and Y5 to the Z were positive while Y4 and Y2 were weakly negative. The total effects (directs plus indirects) of the components were positively contributed to the Z by path analyses. The seed yield components Y1, Y2, Y4 and Y5 were significantly (P<0.001) correlated with the Z for 4 years totally, while in the individual years, Y2 were not significant correlated with Y3, Y4 and Y5 by Peason correlation analyses in the five components in the plant seed production. Therefore, selection for high seed yield through direct selection for large Y1, Y2 and Y3 would be effective for breeding programs in grasses. Furthermore, it is the most important that, via ridge regression, a steady algorithm model between Z and the five yield components was founded, which can be closely estimated the seed yield via the components

Geographical gradient of the <em>eIF4E</em> alleles conferring resistance to potyviruses in pea (<em>Pisum</em>) germplasm

<div><p>Background</p><p>The eukaryotic translation initiation factor 4E was shown to be involved in resistance against several potyviruses in plants, including pea. We combined our knowledge of pea germplasm diversity with that of the <i>eIF4E</i> gene to identify novel genetic diversity.</p><p>Methodology/Principal findings</p><p>Germplasm of 2803 pea accessions was screened for <i>eIF4E</i> intron 3 length polymorphism, resulting in the detection of four <i>eIF4E^A-B-C-S</i> variants, whose distribution was geographically structured. The <i>eIF4E^A</i> variant conferring resistance to the P1 PSbMV pathotype was found in 53 accessions (1.9%), of which 15 were landraces from India, Afghanistan, Nepal, and 7 were from Ethiopia. A newly discovered variant, <i>eIF4E^B</i>, was present in 328 accessions (11.7%) from Ethiopia (29%), Afghanistan (23%), India (20%), Israel (25%) and China (39%). The <i>eIF4E^C</i> variant was detected in 91 accessions (3.2% of total) from India (20%), Afghanistan (33%), the Iberian Peninsula (22%) and the Balkans (9.3%). The <i>eIF4E^S</i> variant for susceptibility predominated as the wild type. Sequencing of 73 samples, identified 34 alleles at the whole gene, 26 at cDNA and 19 protein variants, respectively. Fifteen alleles were virologically tested and 9 alleles (<i>eIF4E^{A-1-2-3-4-5-6-7}</i>, <i>eIF4E^B-1</i>, <i>eIF4E^C-2</i>) conferred resistance to the P1 PSbMV pathotype.</p><p>Conclusions/Significance</p><p>This work identified novel <i>eIF4E</i> alleles within geographically structured pea germplasm and indicated their independent evolution from the susceptible <i>eIF4E^S1</i> allele. Despite high variation present in wild <i>Pisum</i> accessions, none of them possessed resistance alleles, supporting a hypothesis of distinct mode of evolution of resistance in wild as opposed to crop species. The Highlands of Central Asia, the northern regions of the Indian subcontinent, Eastern Africa and China were identified as important centers of pea diversity that correspond with the diversity of the pathogen. The series of alleles identified in this study provides the basis to study the co-evolution of potyviruses and the pea host.</p></div