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

BeadArray Expression Analysis Using Bioconductor

Illumina whole-genome expression BeadArrays are a popular choice in gene profiling studies. Aside from the vendor-provided software tools for analyzing BeadArray expression data (GenomeStudio/BeadStudio), there exists a comprehensive set of open-source analysis tools in the Bioconductor project, many of which have been tailored to exploit the unique properties of this platform. In this article, we explore a number of these software packages and demonstrate how to perform a complete analysis of BeadArray data in various formats. The key steps of importing data, performing quality assessments, preprocessing, and annotation in the common setting of assessing differential expression in designed experiments will be covered

Consistent phenological shifts in the making of a biodiversity hotspot: the Cape flora

Background The best documented survival responses of organisms to past climate change on short (glacial-interglacial) timescales are distributional shifts. Despite ample evidence on such timescales for local adaptations of populations at specific sites, the long-term impacts of such changes on evolutionary significant units in response to past climatic change have been little documented. Here we use phylogenies to reconstruct changes in distribution and flowering ecology of the Cape flora - South Africa's biodiversity hotspot - through a period of past (Neogene and Quaternary) changes in the seasonality of rainfall over a timescale of several million years. Results Forty-three distributional and phenological shifts consistent with past climatic change occur across the flora, and a comparable number of clades underwent adaptive changes in their flowering phenology (9 clades; half of the clades investigated) as underwent distributional shifts (12 clades; two thirds of the clades investigated). Of extant Cape angiosperm species, 14-41% have been contributed by lineages that show distributional shifts consistent with past climate change, yet a similar proportion (14-55%) arose from lineages that shifted flowering phenology. Conclusions Adaptive changes in ecology at the scale we uncover in the Cape and consistent with past climatic change have not been documented for other floras. Shifts in climate tolerance appear to have been more important in this flora than is currently appreciated, and lineages that underwent such shifts went on to contribute a high proportion of the flora's extant species diversity. That shifts in phenology, on an evolutionary timescale and on such a scale, have not yet been detected for other floras is likely a result of the method used; shifts in flowering phenology cannot be detected in the fossil record

Butyrate Transcriptionally Enhances Peptide Transporter PepT1 Expression and Activity

Background: PepT1, an intestinal epithelial apical di/tripeptide transporter, is normally expressed in the small intestine and induced in colon during chronic inflammation. This study aimed at investigating PepT1 regulation by butyrate, a short-chain fatty acid produced by commensal bacteria and accumulated inside inflamed colonocyte. Results: We found that butyrate treatment of human intestinal epithelial Caco2-BBE cells increased human PepT1 (hPepT1) promoter activity in a dose- and time-dependent manner, with maximal activity observed in cells treated with 5 mM butyrate for 24 h. Under this condition, hPepT1 promoter activity, mRNA and protein expression levels were increased as assessed by luciferase assay, real-time RT-PCR and Western blot, respectively. hPepT1 transport activity was accordingly increased by,2.5-fold. Butyrate did not alter hPepT1 mRNA half-life indicating that butyrate acts at the transcriptional level. Molecular analyses revealed that Cdx2 is the most important transcription factor for butyrate-induced increase of hPepT1 expression and activity in Caco2-BBE cells. Butyrate-activated Cdx2 binding to hPepT1 promoter was confirmed by gel shift and chromatin immunoprecipitation. Moreover, Caco2-BBE cells overexpressing Cdx2 exhibited greater hPepT1 expression level than wild-type cells. Finally, treatment of mice with 5 mM butyrate added to drinking water for 24 h increased colonic PepT1 mRNA and protein expression levels, as well as enhanced PepT1 transport activity in colonic apical membranes vesicles

QTL meta-analysis of root traits in Brassica napus under contrasting phosphorus supply in two growth systems

A high-density SNP-based genetic linkage map was constructed and integrated with a previous map in the Tapidor x Ningyou7 (TNDH) Brassica napus population, giving a new map with a total of 2041 molecular markers and an average marker density which increased from 0.39 to 0.97 (0.82 SNP bin) per cM. Root and shoot traits were screened under low and ‘normal’ phosphate (Pi) supply using a ‘pouch and wick’ system, and had been screened previously in an agar based system. The P-efficient parent Ningyou7 had a shorter primary root length (PRL), greater lateral root density (LRD) and a greater shoot biomass than the P-inefficient parent Tapidor under both treatments and growth systems. Quantitative trait loci (QTL) analysis identified a total of 131 QTL, and QTL meta-analysis found four integrated QTL across the growth systems. Integration reduced the confidence interval by ~41%. QTL for root and shoot biomass were co-located on chromosome A3 and for lateral root emergence were co-located on chromosomes A4/C4 and C8/C9. There was a major QTL for LRD on chromosome C9 explaining ~18% of the phenotypic variation. QTL underlying an increased LRD may be a useful breeding target for P uptake efficiency in Brassica

Cancer pharmacogenetics

The large number of active combination chemotherapy regimens for most cancers has led to the need for better information to guide the \u27standard\u27 treatment for each patient. In an attempt to individualise therapy, pharmacogenetics and pharmacogenomics (a polygenic approach to pharmacogenetic studies) encompass the search for answers to the hereditary basis for interindividual differences in drug response. This review will focus on the results of studies assessing the effects of polymorphisms in drug-metabolising enzymes and drug targets on the toxicity and response to commonly used chemotherapy drugs. In addition, the need for polygenic pharmacogenomic strategies to identify patients at risk for adverse drug reactions will be highlighted

Domain-and species-specific monoclonal antibodies recognize the Von Willebrand Factor-C domain of CCN5

The CCN family of proteins typically consists of four distinct peptide domains: an insulin-like growth factor binding protein-type (IGFBP) domain, a Von Willebrand Factor C (VWC) domain, a thrombospondin type 1 repeat (TSP1) domain, and a carboxy-terminal (CT) domain. The six family members participate in many processes, including proliferation, motility, cell-matrix signaling, angiogenesis, and wound healing. Accumulating evidence suggests that truncated and alternatively spliced isoforms are responsible for the diverse functions of CCN proteins in both normal and pathophysiologic states. Analysis of the properties and functions of individual CCN domains further corroborates this idea. CCN5 is unique among the CCN family members because it lacks the CT-domain. To dissect the domain functions of CCN5, we are developing domain-specific mouse monoclonal antibodies. Monoclonal antibodies have the advantages of great specificity, reproducibility, and ease of long-term storage and production. In this communication, we injected mixtures of GST-fused rat CCN5 domains into mice to generate monoclonal antibodies. To identify the domains recognized by the antibodies, we constructed serial expression plasmids that express dual-tagged rat CCN5 domains. All of the monoclonal antibodies generated to date recognize the VWC domain, indicating it is the most highly immunogenic of the CCN5 domains. We characterized one particular clone, 22H10, and found that it recognizes mouse and rat CCN5, but not human recombinant CCN5. Purified 22H10 was successfully applied in Western Blot analysis, immunofluorescence of cultured cells and tissues, and immunoprecipitation, indicating that it will be a useful tool for domain analysis and studies of mouse-human tumor models