Use of DNA–Damaging Agents and RNA Pooling to Assess Expression Profiles Associated with BRCA1 and BRCA2 Mutation Status in Familial Breast Cancer Patients

A large number of rare sequence variants of unknown clinical significance have been identified in the breast cancer susceptibility genes, BRCA1 and BRCA2. Laboratory-based methods that can distinguish between carriers of pathogenic mutations and non-carriers are likely to have utility for the classification of these sequence variants. To identify predictors of pathogenic mutation status in familial breast cancer patients, we explored the use of gene expression arrays to assess the effect of two DNA–damaging agents (irradiation and mitomycin C) on cellular response in relation to BRCA1 and BRCA2 mutation status. A range of regimes was used to treat 27 lymphoblastoid cell-lines (LCLs) derived from affected women in high-risk breast cancer families (nine BRCA1, nine BRCA2, and nine non-BRCA1/2 or BRCAX individuals) and nine LCLs from healthy individuals. Using an RNA–pooling strategy, we found that treating LCLs with 1.2 µM mitomycin C and measuring the gene expression profiles 1 hour post-treatment had the greatest potential to discriminate BRCA1, BRCA2, and BRCAX mutation status. A classifier was built using the expression profile of nine QRT–PCR validated genes that were associated with BRCA1, BRCA2, and BRCAX status in RNA pools. These nine genes could distinguish BRCA1 from BRCA2 carriers with 83% accuracy in individual samples, but three-way analysis for BRCA1, BRCA2, and BRCAX had a maximum of 59% prediction accuracy. Our results suggest that, compared to BRCA1 and BRCA2 mutation carriers, non-BRCA1/2 (BRCAX) individuals are genetically heterogeneous. This study also demonstrates the effectiveness of RNA pools to compare the expression profiles of cell-lines from BRCA1, BRCA2, and BRCAX cases after treatment with irradiation and mitomycin C as a method to prioritize treatment regimes for detailed downstream expression analysis

Probe-level linear model fitting and mixture modeling results in high accuracy detection of differential gene expression

BACKGROUND: The identification of differentially expressed genes (DEGs) from Affymetrix GeneChips arrays is currently done by first computing expression levels from the low-level probe intensities, then deriving significance by comparing these expression levels between conditions. The proposed PL-LM (Probe-Level Linear Model) method implements a linear model applied on the probe-level data to directly estimate the treatment effect. A finite mixture of Gaussian components is then used to identify DEGs using the coefficients estimated by the linear model. This approach can readily be applied to experimental design with or without replication. RESULTS: On a wholly defined dataset, the PL-LM method was able to identify 75% of the differentially expressed genes within 10% of false positives. This accuracy was achieved both using the three replicates per conditions available in the dataset and using only one replicate per condition. CONCLUSION: The method achieves, on this dataset, a higher accuracy than the best set of tools identified by the authors of the dataset, and does so using only one replicate per condition

A novel approach to the clustering of microarray data via nonparametric density estimation

<p>Abstract</p> <p>Background</p> <p>Cluster analysis is a crucial tool in several biological and medical studies dealing with microarray data. Such studies pose challenging statistical problems due to dimensionality issues, since the number of variables can be much higher than the number of observations.</p> <p>Results</p> <p>Here, we present a general framework to deal with the clustering of microarray data, based on a three-step procedure: (i) gene filtering; (ii) dimensionality reduction; (iii) clustering of observations in the reduced space. Via a nonparametric model-based clustering approach we obtain promising results both in simulated and real data.</p> <p>Conclusions</p> <p>The proposed algorithm is a simple and effective tool for the clustering of microarray data, in an unsupervised setting.</p

EzArray: A web-based highly automated Affymetrix expression array data management and analysis system

<p>Abstract</p> <p>Background</p> <p>Though microarray experiments are very popular in life science research, managing and analyzing microarray data are still challenging tasks for many biologists. Most microarray programs require users to have sophisticated knowledge of mathematics, statistics and computer skills for usage. With accumulating microarray data deposited in public databases, easy-to-use programs to re-analyze previously published microarray data are in high demand.</p> <p>Results</p> <p>EzArray is a web-based Affymetrix expression array data management and analysis system for researchers who need to organize microarray data efficiently and get data analyzed instantly. EzArray organizes microarray data into projects that can be analyzed online with predefined or custom procedures. EzArray performs data preprocessing and detection of differentially expressed genes with statistical methods. All analysis procedures are optimized and highly automated so that even novice users with limited pre-knowledge of microarray data analysis can complete initial analysis quickly. Since all input files, analysis parameters, and executed scripts can be downloaded, EzArray provides maximum reproducibility for each analysis. In addition, EzArray integrates with Gene Expression Omnibus (GEO) and allows instantaneous re-analysis of published array data.</p> <p>Conclusion</p> <p>EzArray is a novel Affymetrix expression array data analysis and sharing system. EzArray provides easy-to-use tools for re-analyzing published microarray data and will help both novice and experienced users perform initial analysis of their microarray data from the location of data storage. We believe EzArray will be a useful system for facilities with microarray services and laboratories with multiple members involved in microarray data analysis. EzArray is freely available from <url>http://www.ezarray.com/</url>.</p

A new multitest correction (SGoF) that increases its statistical power when increasing the number of tests

<p>Abstract</p> <p>Background</p> <p>The detection of true significant cases under multiple testing is becoming a fundamental issue when analyzing high-dimensional biological data. Unfortunately, known multitest adjustments reduce their statistical power as the number of tests increase. We propose a new multitest adjustment, based on a sequential goodness of fit metatest (SGoF), which increases its statistical power with the number of tests. The method is compared with Bonferroni and FDR-based alternatives by simulating a multitest context via two different kinds of tests: 1) one-sample t-test, and 2) homogeneity G-test.</p> <p>Results</p> <p>It is shown that SGoF behaves especially well with small sample sizes when 1) the alternative hypothesis is weakly to moderately deviated from the null model, 2) there are widespread effects through the family of tests, and 3) the number of tests is large.</p> <p>Conclusion</p> <p>Therefore, SGoF should become an important tool for multitest adjustment when working with high-dimensional biological data.</p

A New Method for Isolation of Interstitial Fluid from Human Solid Tumors Applied to Proteomic Analysis of Ovarian Carcinoma Tissue

Major efforts have been invested in the identification of cancer biomarkers in plasma, but the extraordinary dynamic range in protein composition, and the dilution of disease specific proteins make discovery in plasma challenging. Focus is shifting towards using proximal fluids for biomarker discovery, but methods to verify the isolated sample's origin are missing. We therefore aimed to develop a technique to search for potential candidate proteins in the proximal proteome, i.e. in the tumor interstitial fluid, since the biomarkers are likely to be excreted or derive from the tumor microenvironment. Since tumor interstitial fluid is not readily accessible, we applied a centrifugation method developed in experimental animals and asked whether interstitial fluid from human tissue could be isolated, using ovarian carcinoma as a model. Exposure of extirpated tissue to 106 g enabled tumor fluid isolation. The fluid was verified as interstitial by an isolated fluid:plasma ratio not significantly different from 1.0 for both creatinine and Na+, two substances predominantly present in interstitial fluid. The isolated fluid had a colloid osmotic pressure 79% of that in plasma, suggesting that there was some sieving of proteins at the capillary wall. Using a proteomic approach we detected 769 proteins in the isolated interstitial fluid, sixfold higher than in patient plasma. We conclude that the isolated fluid represents undiluted interstitial fluid and thus a subproteome with high concentration of locally secreted proteins that may be detected in plasma for diagnostic, therapeutic and prognostic monitoring by targeted methods

A Bayesian Partition Method for Detecting Pleiotropic and Epistatic eQTL Modules

Studies of the relationship between DNA variation and gene expression variation, often referred to as “expression quantitative trait loci (eQTL) mapping”, have been conducted in many species and resulted in many significant findings. Because of the large number of genes and genetic markers in such analyses, it is extremely challenging to discover how a small number of eQTLs interact with each other to affect mRNA expression levels for a set of co-regulated genes. We present a Bayesian method to facilitate the task, in which co-expressed genes mapped to a common set of markers are treated as a module characterized by latent indicator variables. A Markov chain Monte Carlo algorithm is designed to search simultaneously for the module genes and their linked markers. We show by simulations that this method is more powerful for detecting true eQTLs and their target genes than traditional QTL mapping methods. We applied the procedure to a data set consisting of gene expression and genotypes for 112 segregants of S. cerevisiae. Our method identified modules containing genes mapped to previously reported eQTL hot spots, and dissected these large eQTL hot spots into several modules corresponding to possibly different biological functions or primary and secondary responses to regulatory perturbations. In addition, we identified nine modules associated with pairs of eQTLs, of which two have been previously reported. We demonstrated that one of the novel modules containing many daughter-cell expressed genes is regulated by AMN1 and BPH1. In conclusion, the Bayesian partition method which simultaneously considers all traits and all markers is more powerful for detecting both pleiotropic and epistatic effects based on both simulated and empirical data

RNA-seq Analysis Reveals That an ECF σ Factor, AcsS, Regulates Achromobactin Biosynthesis in Pseudomonas syringae pv. syringae B728a

Iron is an essential micronutrient for Pseudomonas syringae pv. syringae strain B728a and many other microorganisms; therefore, B728a has evolved methods of iron acquirement including the use of iron-chelating siderophores. In this study an extracytoplasmic function (ECF) sigma factor, AcsS, encoded within the achromobactin gene cluster is shown to be a major regulator of genes involved in the biosynthesis and secretion of this siderophore. However, production of achromobactin was not completely abrogated in the deletion mutant, implying that other regulators may be involved such as PvdS, the sigma factor that regulates pyoverdine biosynthesis. RNA-seq analysis identified 287 genes that are differentially expressed between the AcsS deletion mutant and the wild type strain. These genes are involved in iron response, secretion, extracellular polysaccharide production, and cell motility. Thus, the transcriptome analysis supports a role for AcsS in the regulation of achromobactin production and the potential activity of both AcsS and achromobactin in the plant-associated lifestyle of strain B728a