Appearance frequency modulated gene set enrichment testing

Abstract Background Gene set enrichment testing has helped bridge the gap from an individual gene to a systems biology interpretation of microarray data. Although gene sets are defined a priori based on biological knowledge, current methods for gene set enrichment testing treat all genes equal. It is well-known that some genes, such as those responsible for housekeeping functions, appear in many pathways, whereas other genes are more specialized and play a unique role in a single pathway. Drawing inspiration from the field of information retrieval, we have developed and present here an approach to incorporate gene appearance frequency (in KEGG pathways) into two current methods, Gene Set Enrichment Analysis (GSEA) and logistic regression-based LRpath framework, to generate more reproducible and biologically meaningful results. Results Two breast cancer microarray datasets were analyzed to identify gene sets differentially expressed between histological grade 1 and 3 breast cancer. The correlation of Normalized Enrichment Scores (NES) between gene sets, generated by the original GSEA and GSEA with the appearance frequency of genes incorporated (GSEA-AF), was compared. GSEA-AF resulted in higher correlation between experiments and more overlapping top gene sets. Several cancer related gene sets achieved higher NES in GSEA-AF as well. The same datasets were also analyzed by LRpath and LRpath with the appearance frequency of genes incorporated (LRpath-AF). Two well-studied lung cancer datasets were also analyzed in the same manner to demonstrate the validity of the method, and similar results were obtained. Conclusions We introduce an alternative way to integrate KEGG PATHWAY information into gene set enrichment testing. The performance of GSEA and LRpath can be enhanced with the integration of appearance frequency of genes. We conclude that, generally, gene set analysis methods with the integration of information from KEGG PATHWAY performs better both statistically and biologically.http://deepblue.lib.umich.edu/bitstream/2027.42/112430/1/12859_2010_Article_4457.pd

Bisphenol A‐associated alterations in the expression and epigenetic regulation of genes encoding xenobiotic metabolizing enzymes in human fetal liver

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106679/1/em21823.pd

Intensity-based hierarchical Bayes method improves testing for differentially expressed genes in microarray experiments

BACKGROUND: The small sample sizes often used for microarray experiments result in poor estimates of variance if each gene is considered independently. Yet accurately estimating variability of gene expression measurements in microarray experiments is essential for correctly identifying differentially expressed genes. Several recently developed methods for testing differential expression of genes utilize hierarchical Bayesian models to "pool" information from multiple genes. We have developed a statistical testing procedure that further improves upon current methods by incorporating the well-documented relationship between the absolute gene expression level and the variance of gene expression measurements into the general empirical Bayes framework. RESULTS: We present a novel Bayesian moderated-T, which we show to perform favorably in simulations, with two real, dual-channel microarray experiments and in two controlled single-channel experiments. In simulations, the new method achieved greater power while correctly estimating the true proportion of false positives, and in the analysis of two publicly-available "spike-in" experiments, the new method performed favorably compared to all tested alternatives. We also applied our method to two experimental datasets and discuss the additional biological insights as revealed by our method in contrast to the others. The R-source code for implementing our algorithm is freely available at . CONCLUSION: We use a Bayesian hierarchical normal model to define a novel Intensity-Based Moderated T-statistic (IBMT). The method is completely data-dependent using empirical Bayes philosophy to estimate hyperparameters, and thus does not require specification of any free parameters. IBMT has the strength of balancing two important factors in the analysis of microarray data: the degree of independence of variances relative to the degree of identity (i.e. t-tests vs. equal variance assumption), and the relationship between variance and signal intensity. When this variance-intensity relationship is weak or does not exist, IBMT reduces to a previously described moderated t-statistic. Furthermore, our method may be directly applied to any array platform and experimental design. Together, these properties show IBMT to be a valuable option in the analysis of virtually any microarray experiment

Convergence of genetic influences in comorbidity

Abstract Background Predisposition to complex diseases is explained in part by genetic variation, and complex diseases are frequently comorbid, consistent with pleiotropic genetic variation influencing comorbidity. Genome Wide Association (GWA) studies typically assess association between SNPs and a single-disease phenotype. Fisher meta-analysis combines evidence of association from single-disease GWA studies, assuming that each study is an independent test of the same hypothesis. The Rank Product (RP) method overcomes limitations posed by Fisher assumptions, though RP was not designed for GWA data. Methods We modified RP to accommodate GWA data, and we call it modRP. Using p-values output from GWA studies, we aggregate evidence for association between SNPs and related phenotypes. To assess significance, RP randomly samples the observed ranks to develop the null distribution of the RP statistic, and then places the observed RPs into the null distribution. ModRP eliminates the effect of linkage disequilibrium and controls for differences in power at tested SNPs, to meet RP assumptions in application to GWA data. Results After validating modRP based on both positive and negative control studies, we searched for pleiotropic influences on comorbid substance use disorders in a novel study, and found two SNPs to be significantly associated with comorbid cocaine, opium, and nicotine dependence. Placing these SNPs into biological context, we developed a protein network modeling the interaction of cocaine, nicotine, and opium with these variants. Conclusions ModRP is a novel approach to identifying pleiotropic genetic influences on comorbid complex diseases. It can be used to assess association for related phenotypes where raw data is unavailable or inappropriate for analysis using other approaches. The method is conceptually simple and produces statistically significant, biologically relevant results.http://deepblue.lib.umich.edu/bitstream/2027.42/112931/1/12859_2012_Article_5068.pd

Epigenome‐wide DNA methylation analysis implicates neuronal and inflammatory signaling pathways in adult murine hepatic tumorigenesis following perinatal exposure to bisphenol A

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/133545/1/em22024.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/133545/2/em22024_am.pd

LRpath analysis reveals common pathways dysregulated via DNA methylation across cancer types

Abstract Background The relative contribution of epigenetic mechanisms to carcinogenesis is not well understood, including the extent to which epigenetic dysregulation and somatic mutations target similar genes and pathways. We hypothesize that during carcinogenesis, certain pathways or biological gene sets are commonly dysregulated via DNA methylation across cancer types. The ability of our logistic regression-based gene set enrichment method to implicate important biological pathways in high-throughput data is well established. Results We developed a web-based gene set enrichment application called LRpath with clustering functionality that allows for identification and comparison of pathway signatures across multiple studies. Here, we employed LRpath analysis to unravel the commonly altered pathways and other gene sets across ten cancer studies employing DNA methylation data profiled with the Illumina HumanMethylation27 BeadChip. We observed a surprising level of concordance in differential methylation across multiple cancer types. For example, among commonly hypomethylated groups, we identified immune-related functions, peptidase activity, and epidermis/keratinocyte development and differentiation. Commonly hypermethylated groups included homeobox and other DNA-binding genes, nervous system and embryonic development, and voltage-gated potassium channels. For many gene sets, we observed significant overlap in the specific subset of differentially methylated genes. Interestingly, fewer DNA repair genes were differentially methylated than expected by chance. Conclusions Clustering analysis performed with LRpath revealed tightly clustered concepts enriched for differential methylation. Several well-known cancer-related pathways were significantly affected, while others were depleted in differential methylation. We conclude that DNA methylation changes in cancer tend to target a subset of the known cancer pathways affected by genetic aberrations.http://deepblue.lib.umich.edu/bitstream/2027.42/112789/1/12864_2012_Article_4373.pd

Stability of methylation markers in head and neck squamous cell carcinoma

BackgroundAs cancer progresses, methylation patterns change to promote the tumorigenic phenotype. However, stability of methylation markers over time and the extent that biopsy samples are representative of larger tumor specimens are unknown. This information is critical for clinical use of such biomarkers.MethodsNinety‐eight patients with tumor specimens from 2 timepoints were measured for DNA methylation in the promoter regions across 4 genes.ResultsThere were no significant differences in overall methylation of CCNA1 (cyclin A1), NDN (necdin), deleted in colorectal carcinoma (DCC), and cluster of differentiation 1a (CD1A) within paired specimens (p values = .56, .17, .66, and .58, respectively). All genes showed strong correlations between paired specimens across time. Methylation was most consistent for CCNA1 and NDN over time.ConclusionThis report provides the first evidence that methylation markers measured in biopsy samples are representative of gene methylation in later specimens and suggests that biopsy markers could be representative biomarkers for use in defining personalized treatment utilizing epigenetic changes. © 2015 Wiley Periodicals, Head Neck 38: E1325–E1331, 2016Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/137576/1/hed24223.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/137576/2/hed24223_am.pd

A new method to remove hybridization bias for interspecies comparison of global gene expression profiles uncovers an association between mRNA sequence divergence and differential gene expression in Xenopus

The recent sequencing of a large number of Xenopus tropicalis expressed sequences has allowed development of a high-throughput approach to study Xenopus global RNA gene expression. We examined the global gene expression similarities and differences between the historically significant Xenopus laevis model system and the increasingly used X.tropicalis model system and assessed whether an X.tropicalis microarray platform can be used for X.laevis. These closely related species were also used to investigate a more general question: is there an association between mRNA sequence divergence and differences in gene expression levels? We carried out a comprehensive comparison of global gene expression profiles using microarrays of different tissues and developmental stages of X.laevis and X.tropicalis. We (i) show that the X.tropicalis probes provide an efficacious microarray platform for X.laevis, (ii) describe methods to compare interspecies mRNA profiles that correct differences in hybridization efficiency and (iii) show independently of hybridization bias that as mRNA sequence divergence increases between X.laevis and X.tropicalis differences in mRNA expression levels also increase

Gut microbiota induce IGF-1 and promote bone formation and growth

New interventions are needed to improve bone health and reduce the risk for osteoporosis and fracture. Dysbiosis is increasingly linked to metabolic abnormalities, although the effect of the microbiota on skeletal health is poorly understood. Previous studies suggest microbiota are detrimental to bone by increasing resorption. In this report, we show that the gut resident microbiota promote bone formation, as well as resorption, with long-term exposure to microbiota resulting in net skeletal growth. Microbiota induce the hormone insulin-like growth factor 1 (IGF-1), which promotes bone growth and remodeling. Short-chain fatty acids (SCFAs), produced when microbiota ferment fiber, also induce IGF-1, suggesting a mechanism by which microbiota affect bone health. Manipulating the microbiome or its metabolites may afford opportunities to optimize bone health and growth

SEC24A deficiency lowers plasma cholesterol through reduced PCSK9 secretion.

The secretory pathway of eukaryotic cells packages cargo proteins into COPII-coated vesicles for transport from the endoplasmic reticulum (ER) to the Golgi. We now report that complete genetic deficiency for the COPII component SEC24A is compatible with normal survival and development in the mouse, despite the fundamental role of SEC24 in COPII vesicle formation and cargo recruitment. However, these animals exhibit markedly reduced plasma cholesterol, with mutations in Apoe and Ldlr epistatic to Sec24a, suggesting a receptor-mediated lipoprotein clearance mechanism. Consistent with these data, hepatic LDLR levels are up-regulated in SEC24A-deficient cells as a consequence of specific dependence of PCSK9, a negative regulator of LDLR, on SEC24A for efficient exit from the ER. Our findings also identify partial overlap in cargo selectivity between SEC24A and SEC24B, suggesting a previously unappreciated heterogeneity in the recruitment of secretory proteins to the COPII vesicles that extends to soluble as well as trans-membrane cargoes. DOI:http://dx.doi.org/10.7554/eLife.00444.001