Bayesian Gaussian Copula Factor Models for Mixed Data.

Gaussian factor models have proven widely useful for parsimoniously characterizing dependence in multivariate data. There is a rich literature on their extension to mixed categorical and continuous variables, using latent Gaussian variables or through generalized latent trait models acommodating measurements in the exponential family. However, when generalizing to non-Gaussian measured variables the latent variables typically influence both the dependence structure and the form of the marginal distributions, complicating interpretation and introducing artifacts. To address this problem we propose a novel class of Bayesian Gaussian copula factor models which decouple the latent factors from the marginal distributions. A semiparametric specification for the marginals based on the extended rank likelihood yields straightforward implementation and substantial computational gains. We provide new theoretical and empirical justifications for using this likelihood in Bayesian inference. We propose new default priors for the factor loadings and develop efficient parameter-expanded Gibbs sampling for posterior computation. The methods are evaluated through simulations and applied to a dataset in political science. The models in this paper are implemented in the R package bfa.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

Bayesian inference of the number of factors in gene-expression analysis: application to human virus challenge studies.

BACKGROUND: Nonparametric Bayesian techniques have been developed recently to extend the sophistication of factor models, allowing one to infer the number of appropriate factors from the observed data. We consider such techniques for sparse factor analysis, with application to gene-expression data from three virus challenge studies. Particular attention is placed on employing the Beta Process (BP), the Indian Buffet Process (IBP), and related sparseness-promoting techniques to infer a proper number of factors. The posterior density function on the model parameters is computed using Gibbs sampling and variational Bayesian (VB) analysis. RESULTS: Time-evolving gene-expression data are considered for respiratory syncytial virus (RSV), Rhino virus, and influenza, using blood samples from healthy human subjects. These data were acquired in three challenge studies, each executed after receiving institutional review board (IRB) approval from Duke University. Comparisons are made between several alternative means of per-forming nonparametric factor analysis on these data, with comparisons as well to sparse-PCA and Penalized Matrix Decomposition (PMD), closely related non-Bayesian approaches. CONCLUSIONS: Applying the Beta Process to the factor scores, or to the singular values of a pseudo-SVD construction, the proposed algorithms infer the number of factors in gene-expression data. For real data the "true" number of factors is unknown; in our simulations we consider a range of noise variances, and the proposed Bayesian models inferred the number of factors accurately relative to other methods in the literature, such as sparse-PCA and PMD. We have also identified a "pan-viral" factor of importance for each of the three viruses considered in this study. We have identified a set of genes associated with this pan-viral factor, of interest for early detection of such viruses based upon the host response, as quantified via gene-expression data.

Growth factor erv1-like modulates Drp1 to preserve mitochondrial dynamics and function in mouse embryonic stem cells.

The relationship of mitochondrial dynamics and function to pluripotency are rather poorly understood aspects of stem cell biology. Here we show that growth factor erv1-like (Gfer) is involved in preserving mouse embryonic stem cell (ESC) mitochondrial morphology and function. Knockdown (KD) of Gfer in ESCs leads to decreased pluripotency marker expression, embryoid body (EB) formation, cell survival, and loss of mitochondrial function. Mitochondria in Gfer-KD ESCs undergo excessive fragmentation and mitophagy, whereas those in ESCs overexpressing Gfer appear elongated. Levels of the mitochondrial fission GTPase dynamin-related protein 1 (Drp1) are highly elevated in Gfer-KD ESCs and decreased in Gfer-overexpressing cells. Treatment with a specific inhibitor of Drp1 rescues mitochondrial function and apoptosis, whereas expression of Drp1-dominant negative resulted in the restoration of pluripotency marker expression in Gfer-KD ESCs. Altogether, our data reveal a novel prosurvival role for Gfer in maintaining mitochondrial fission-fusion dynamics in pluripotent ESCs.

MYC activity mitigates response to rapamycin in prostate cancer through eukaryotic initiation factor 4E-binding protein 1-mediated inhibition of autophagy.

Loss of PTEN and activation of phosphoinositide 3-kinase are commonly observed in advanced prostate cancer. Inhibition of mammalian target of rapamycin (mTOR), a downstream target of phosphoinositide 3-kinase signaling, results in cell cycle arrest and apoptosis in multiple in vitro and in vivo models of prostate cancer. However, single-agent use of mTOR inhibition has limited clinical success, and the identification of molecular events mitigating tumor response to mTOR inhibition remains a critical question. Here, using genetically engineered human prostate epithelial cells (PrEC), we show that MYC, a frequent target of genetic gain in prostate cancers, abrogates sensitivity to rapamycin by decreasing rapamycin-induced cytostasis and autophagy. Analysis of MYC and the mTOR pathway in human prostate tumors and PrEC showed selective increased expression of eukaryotic initiation factor 4E-binding protein 1 (4EBP1) with gain in MYC copy number or forced MYC expression, respectively. We have also found that MYC binds to regulatory regions of the 4EBP1 gene. Suppression of 4EBP1 expression resulted in resensitization of MYC-expressing PrEC to rapamycin and increased autophagy. Taken together, our findings suggest that MYC expression abrogates sensitivity to rapamycin through increased expression of 4EBP1 and reduced autophagy.

Implant healing in experimental animal models of diabetes.

Diabetes mellitus is becoming increasingly prevalent worldwide. Additionally, there is an increasing number of patients receiving implantable devices such as glucose sensors and orthopedic implants. Thus, it is likely that the number of diabetic patients receiving these devices will also increase. Even though implantable medical devices are considered biocompatible by the Food and Drug Administration, the adverse tissue healing that occurs adjacent to these foreign objects is a leading cause of their failure. This foreign body response leads to fibrosis, encapsulation of the device, and a reduction or cessation of device performance. A second adverse event is microbial infection of implanted devices, which can lead to persistent local and systemic infections and also exacerbates the fibrotic response. Nearly half of all nosocomial infections are associated with the presence of an indwelling medical device. Events associated with both the foreign body response and implant infection can necessitate device removal and may lead to amputation, which is associated with significant morbidity and cost. Diabetes mellitus is generally indicated as a risk factor for the infection of a variety of implants such as prosthetic joints, pacemakers, implantable cardioverter defibrillators, penile implants, and urinary catheters. Implant infection rates in diabetic patients vary depending upon the implant and the microorganism, however, for example, diabetes was found to be a significant variable associated with a nearly 7.2% infection rate for implantable cardioverter defibrillators by the microorganism Candida albicans. While research has elucidated many of the altered mechanisms of diabetic cutaneous wound healing, the internal healing adjacent to indwelling medical devices in a diabetic model has rarely been studied. Understanding this healing process is crucial to facilitating improved device design. The purpose of this article is to summarize the physiologic factors that influence wound healing and infection in diabetic patients, to review research concerning diabetes and biomedical implants and device infection, and to critically analyze which diabetic animal model might be advantageous for assessing internal healing adjacent to implanted devices.

Recovery from an acute infection in C. elegans requires the GATA transcription factor ELT-2.

The mechanisms involved in the recognition of microbial pathogens and activation of the immune system have been extensively studied. However, the mechanisms involved in the recovery phase of an infection are incompletely characterized at both the cellular and physiological levels. Here, we establish a Caenorhabditis elegans-Salmonella enterica model of acute infection and antibiotic treatment for studying biological changes during the resolution phase of an infection. Using whole genome expression profiles of acutely infected animals, we found that genes that are markers of innate immunity are down-regulated upon recovery, while genes involved in xenobiotic detoxification, redox regulation, and cellular homeostasis are up-regulated. In silico analyses demonstrated that genes altered during recovery from infection were transcriptionally regulated by conserved transcription factors, including GATA/ELT-2, FOXO/DAF-16, and Nrf/SKN-1. Finally, we found that recovery from an acute bacterial infection is dependent on ELT-2 activity.

Modeling economy-wide vs sectoral climate policies using combined aggregate-sectoral models

Economic analyses of climate change policies frequently focus on reductions of energy-related carbon dioxide emissions via market-based, economy-wide policies. The current course of environment and energy policy debate in the United States, however, suggests an alternative outcome: sector-based and/or inefficiently designed policies. This paper uses a collection of specialized, sector-based models in conjunction with a computable general equilibrium model of the economy to examine and compare these policies at an aggregate level. We examine the relative cost of different policies designed to achieve the same quantity of emission reductions. We find that excluding a limited number of sectors from an economy-wide policy does not significantly raise costs. Focusing policy solely on the electricity and transportation sectors doubles costs, however, and using non-market policies can raise cost by a factor of ten. These results are driven in part by, and are sensitive to, our modeling of pre-existing tax distortions. Copyright © 2006 by the IAEE. All rights reserved.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.