Probability fold change: a robust computational approach for identifying differentially expressed gene lists

Eight potential biomarkers for distinguishing between lung adenocarcinoma and squamous cell carcinoma

Lung adenocarcinoma (LADC) and squamous cell carcinoma (LSCC) are the most common non-small cell lung cancer histological phenotypes. Accurate diagnosis distinguishing between these two lung cancer types has clinical significance. For this study, we analyzed four Gene Expression Omnibus (GEO) datasets (GSE28571, GSE37745, GSE43580, and GSE50081). We then imported the datasets into the Gene-Cloud of Biotechnology Information online platform to identify genes differentially expressed in LADC and LSCC. We identified DSG3 (desmoglein 3), KRT5 (keratin 5), KRT6A (keratin 6A), KRT6B (keratin 6B), NKX2-1 (NK2 homeobox 1), SFTA2 (surfactant associated 2), SFTA3 (surfactant associated 3), and TMC5 (transmembrane channel-like 5) as potential biomarkers for distinguishing between LADC and LSCC. Receiver operating characteristic curve analysis suggested that KRT5 had the highest diagnostic value for discriminating between these two cancer types. Using the PrognoScan online survival analysis tool and the Kaplan-Meier Plotter, we found that high KRT6A or KRT6B levels, or low NKX2-1, SFTA3, or TMC5 levels correlated with unfavorable prognoses in LADC patients. Further studies will be needed to verify our findings in additional patient samples, and to elucidate the mechanisms of action of these potential biomarkers in non-small cell lung cancer.

Probabilistic strain optimization under constraint uncertainty

Background

An important step in strain optimization is to identify reactions whose activities should be modified to achieve the desired cellular objective. Preferably, these reactions are identified systematically, as the number of possible combinations of reaction modifications could be very large. Over the last several years, a number of computational methods have been described for identifying combinations of reaction modifications. However, none of these methods explicitly address uncertainties in implementing the reaction activity modifications. In this work, we model the uncertainties as probability distributions in the flux carrying capacities of reactions. Based on this model, we develop an optimization method that identifies reactions for flux capacity modifications to predict outcomes with high statistical likelihood.

Results

We compare three optimization methods that select an intervention set comprising up- or down-regulation of reaction flux capacity: CCOpt (Chance constrained optimization), DetOpt (Deterministic optimization), and MCOpt (Monte Carlo-based optimization). We evaluate the methods using a Monte Carlo simulation-based method, MCEval (Monte Carlo Evaluations). We present two case studies analyzing a CHO cell and an adipocyte model. The flux capacity distributions required for our methods were estimated from maximal reaction velocities or elementary mode analysis. The intervention set selected by CCOpt consistently outperforms the intervention set selected by DetOpt in terms of tolerance to flux capacity variations. MCEval shows that the optimal flux predicted based on the CCOpt intervention set is more likely to be obtained, in a probabilistic sense, than the flux predicted by DetOpt. The intervention sets identified by CCOpt and MCOpt were similar; however, the exhaustive sampling required by MCOpt incurred significantly greater computational cost.

Conclusions

Maximizing tolerance to variable engineering outcomes (in modifying enzyme activities) can identify intervention sets that statistically improve the desired cellular objective.

Reliable and sample saving gene expression analysis approach for diagnostic tool development

This work answers the question of whether it is necessary to hybridize individual instead of pooled RNA samples on microarrays for screening gene targets suitable as diagnostic tools for radiation exposure scenarios, while at the same time meeting comparable microarray quality criteria. For developing new clinical diagnostic tools, a two-stage study design was employed in five projects. At first, pooled and not individual RNA samples were hybridized on microarrays for screening purposes. Potential gene candidates were selected based on their fold-change only. This was followed by a validation/quantification step using individual RNA samples and quantitative RT-PCR. Quality criteria from the screening approach with pooled RNA samples were compared with published data from the MicroArray Quality Control (MAQC) consortium that hybridized each reference RNA sample separately and established quality criteria for microarrays. When comparing both approaches, only insignificant differences for quality criteria such as false positives, sensitivity, specificity, and overall agreement were found. However, material, costs, and time were drastically reduced when hybridizing pooled RNA and gene targets applicable for clinical diagnostic purposes could be successfully selected. In search of new diagnostic tools for radiation exposure scenarios, the two stage study design using either pooled or individual RNA samples on microarrays shows comparable quality criteria, but the RNA pooling approach saves unique material, costs, and efforts and successfully selects gene targets that can be used for the desired diagnostic purposes.

CDS: a fold-change based statistical test for concomitant identification of distinctness and similarity in gene expression analysis

The problem of identifying differential activity such as in gene expression is a major defeat in biostatistics and bioinformatics. Equally important, however much less frequently studied, is the question of similar activity from one biological condition to another. The fold-change, or ratio, is usually considered a relevant criterion for stating difference and similarity between measurements. Importantly, no statistical method for concomitant evaluation of similarity and distinctness currently exists for biological applications. Modern microarray, digital PCR (dPCR), and Next-Generation Sequencing (NGS) technologies frequently provide a means of coefficient of variation estimation for individual measurements. Using fold-change, and by making the assumption that measurements are normally distributed with known variances, we designed a novel statistical test that allows us to detect concomitantly, thus using the same formalism, differentially and similarly expressed genes (http://cds.ihes.fr). Given two sets of gene measurements in different biological conditions, the probabilities of making type I and type II errors in stating that a gene is differentially or similarly expressed from one condition to the other can be calculated. Furthermore, a confidence interval for the fold-change can be delineated. Finally, we demonstrate that the assumption of normality can be relaxed to consider arbitrary distributions numerically. The Concomitant evaluation of Distinctness and Similarity (CDS) statistical test correctly estimates similarities and differences between measurements of gene expression. The implementation, being time and memory efficient, allows the use of the CDS test in high-throughput data analysis such as microarray, dPCR, and NGS experiments. Importantly, the CDS test can be applied to the comparison of single measurements (N = 1) provided the variance (or coefficient of variation) of the signals is known, making CDS a valuable tool also in biomedical analysis where typically a single measurement per subject is available.