Improvement of experimental testing and network training conditions with genome-wide microarrays for more accurate predictions of drug gene targets

BACKGROUND: Genome-wide microarrays have been useful for predicting chemical-genetic interactions at the gene level. However, interpreting genome-wide microarray results can be overwhelming due to the vast output of gene expression data combined with off-target transcriptional responses many times induced by a drug treatment. This study demonstrates how experimental and computational methods can interact with each other, to arrive at more accurate predictions of drug-induced perturbations. We present a two-stage strategy that links microarray experimental testing and network training conditions to predict gene perturbations for a drug with a known mechanism of action in a well-studied organism. RESULTS: S. cerevisiae cells were treated with the antifungal, fluconazole, and expression profiling was conducted under different biological conditions using Affymetrix genome-wide microarrays. Transcripts were filtered with a formal network-based method, sparse simultaneous equation models and Lasso regression (SSEM-Lasso), under different network training conditions. Gene expression results were evaluated using both gene set and single gene target analyses, and the drug’s transcriptional effects were narrowed first by pathway and then by individual genes. Variables included: (i) Testing conditions – exposure time and concentration and (ii) Network training conditions – training compendium modifications. Two analyses of SSEM-Lasso output – gene set and single gene – were conducted to gain a better understanding of how SSEM-Lasso predicts perturbation targets. CONCLUSIONS: This study demonstrates that genome-wide microarrays can be optimized using a two-stage strategy for a more in-depth understanding of how a cell manifests biological reactions to a drug treatment at the transcription level. Additionally, a more detailed understanding of how the statistical model, SSEM-Lasso, propagates perturbations through a network of gene regulatory interactions is achieved.Published versio

Stable Feature Selection for Biomarker Discovery

Feature selection techniques have been used as the workhorse in biomarker discovery applications for a long time. Surprisingly, the stability of feature selection with respect to sampling variations has long been under-considered. It is only until recently that this issue has received more and more attention. In this article, we review existing stable feature selection methods for biomarker discovery using a generic hierarchal framework. We have two objectives: (1) providing an overview on this new yet fast growing topic for a convenient reference; (2) categorizing existing methods under an expandable framework for future research and development

Aging-associated changes in microRNA expression profile of internal anal sphincter smooth muscle: Role of microRNA-133a.

A comprehensive genomic and proteomic, computational, and physiological approach was employed to examine the (previously unexplored) role of microRNAs (miRNAs) as regulators of internal anal sphincter (IAS) smooth muscle contractile phenotype and basal tone. miRNA profiling, genome-wide expression, validation, and network analyses were employed to assess changes in mRNA and miRNA expression in IAS smooth muscles from young vs. aging rats. Multiple miRNAs, including rno-miR-1, rno-miR-340-5p, rno-miR-185, rno-miR-199a-3p, rno-miR-200c, rno-miR-200b, rno-miR-31, rno-miR-133a, and rno-miR-206, were found to be upregulated in aging IAS. qPCR confirmed the upregulated expression of these miRNAs and downregulation of multiple, predicted targets (Eln, Col3a1, Col1a1, Zeb2, Myocd, Srf, Smad1, Smad2, Rhoa/Rock2, Fn1, Tagln v2, Klf4, and Acta2) involved in regulation of smooth muscle contractility. Subsequent studies demonstrated an aging-associated increase in the expression of miR-133a, corresponding decreases in RhoA, ROCK2, MYOCD, SRF, and SM22α protein expression, RhoA-signaling, and a decrease in basal and agonist [U-46619 (thromboxane A2 analog)]-induced increase in the IAS tone. Moreover, in vitro transfection of miR-133a caused a dose-dependent increase of IAS tone in strips, which was reversed by anti-miR-133a. Last, in vivo perianal injection of anti-miR-133a reversed the loss of IAS tone associated with age. This work establishes the important regulatory effect of miRNA-133a on basal and agonist-stimulated IAS tone. Moreover, reversal of age-associated loss of tone via anti-miR delivery strongly implicates miR dysregulation as a causal factor in the aging-associated decrease in IAS tone and suggests that miR-133a is a feasible therapeutic target in aging-associated rectoanal incontinence

A machine learning pipeline for discriminant pathways identification

Motivation: Identifying the molecular pathways more prone to disruption during a pathological process is a key task in network medicine and, more in general, in systems biology. Results: In this work we propose a pipeline that couples a machine learning solution for molecular profiling with a recent network comparison method. The pipeline can identify changes occurring between specific sub-modules of networks built in a case-control biomarker study, discriminating key groups of genes whose interactions are modified by an underlying condition. The proposal is independent from the classification algorithm used. Three applications on genomewide data are presented regarding children susceptibility to air pollution and two neurodegenerative diseases: Parkinson's and Alzheimer's. Availability: Details about the software used for the experiments discussed in this paper are provided in the Appendix

An integrative analysis of cancer gene expression studies using Bayesian latent factor modeling

We present an applied study in cancer genomics for integrating data and inferences from laboratory experiments on cancer cell lines with observational data obtained from human breast cancer studies. The biological focus is on improving understanding of transcriptional responses of tumors to changes in the pH level of the cellular microenvironment. The statistical focus is on connecting experimentally defined biomarkers of such responses to clinical outcome in observational studies of breast cancer patients. Our analysis exemplifies a general strategy for accomplishing this kind of integration across contexts. The statistical methodologies employed here draw heavily on Bayesian sparse factor models for identifying, modularizing and correlating with clinical outcome these signatures of aggregate changes in gene expression. By projecting patterns of biological response linked to specific experimental interventions into observational studies where such responses may be evidenced via variation in gene expression across samples, we are able to define biomarkers of clinically relevant physiological states and outcomes that are rooted in the biology of the original experiment. Through this approach we identify microenvironment-related prognostic factors capable of predicting long term survival in two independent breast cancer datasets. These results suggest possible directions for future laboratory studies, as well as indicate the potential for therapeutic advances though targeted disruption of specific pathway components.Comment: Published in at http://dx.doi.org/10.1214/09-AOAS261 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Genetic Variation and Antioxidant Response Gene Expression in the Bronchial Airway Epithelium of Smokers at Risk for Lung Cancer

Prior microarray studies of smokers at high risk for lung cancer have demonstrated that heterogeneity in bronchial airway epithelial cell gene expression response to smoking can serve as an early diagnostic biomarker for lung cancer. As a first step in applying functional genomic analysis to population studies, we have examined the relationship between gene expression variation and genetic variation in a central molecular pathway (NRF2-mediated antioxidant response) associated with smoking exposure and lung cancer. We assessed global gene expression in histologically normal airway epithelial cells obtained at bronchoscopy from smokers who developed lung cancer (SC, n=20), smokers without lung cancer (SNC, n=24), and never smokers (NS, n=8). Functional enrichment analysis showed that the NRF2-mediated, antioxidant response element (ARE)-regulated genes, were significantly lower in SC, when compared with expression levels in SNC. Importantly, we found that the expression of MAFG (a binding partner of NRF2) was correlated with the expression of ARE genes, suggesting MAFG levels may limit target gene induction. Bioinformatically we identified single nucleotide polymorphisms (SNPs) in putative ARE genes and to test the impact of genetic variation, we genotyped these putative regulatory SNPs and other tag SNPs in selected NRF2 pathway genes. Sequencing MAFG locus, we identified 30 novel SNPs and two were associated with either gene expression or lung cancer status among smokers. This work demonstrates an analysis approach that integrates bioinformatics pathway and transcription factor binding site analysis with genotype, gene expression and disease status to identify SNPs that may be associated with individual differences in gene expression and/or cancer status in smokers. These polymorphisms might ultimately contribute to lung cancer risk via their effect on the airway gene expression response to tobacco-smoke exposure.Intramural Research Program of the National Institute of Environmental Health Sciences; National Institutes of Health (Z01 ES100475, U01ES016035, R01CA124640

Hypomethylation of CYP2E1 and DUSP22 Promoters Associated With Disease Activity and Erosive Disease Among Rheumatoid Arthritis Patients.

OBJECTIVE:Epigenetic modifications have previously been associated with rheumatoid arthritis (RA). In this study, we aimed to determine whether differential DNA methylation in peripheral blood cell subpopulations is associated with any of 4 clinical outcomes among RA patients. METHODS:Peripheral blood samples were obtained from 63 patients in the University of California, San Francisco RA cohort (all satisfied the American College of Rheumatology classification criteria; 57 were seropositive for rheumatoid factor and/or anti-cyclic citrullinated protein). Fluorescence-activated cell sorting was used to separate the cells into 4 immune cell subpopulations (CD14+ monocytes, CD19+ B cells, CD4+ naive T cells, and CD4+ memory T cells) per individual, and 229 epigenome-wide DNA methylation profiles were generated using Illumina HumanMethylation450 BeadChips. Differentially methylated positions and regions associated with the Clinical Disease Activity Index score, erosive disease, RA Articular Damage score, Sharp score, medication at time of blood draw, smoking status, and disease duration were identified using robust regression models and empirical Bayes variance estimators. RESULTS:Differential methylation of CpG sites associated with clinical outcomes was observed in all 4 cell types. Hypomethylated regions in the CYP2E1 and DUSP22 gene promoters were associated with active and erosive disease, respectively. Pathway analyses suggested that the biologic mechanisms underlying each clinical outcome are cell type-specific. Evidence of independent effects on DNA methylation from smoking, medication use, and disease duration were also identified. CONCLUSION:Methylation signatures specific to RA clinical outcomes may have utility as biomarkers or predictors of exposure, disease progression, and disease severity