Unsupervised Bayesian linear unmixing of gene expression microarrays

Background: This paper introduces a new constrained model and the corresponding algorithm, called unsupervised Bayesian linear unmixing (uBLU), to identify biological signatures from high dimensional assays like gene expression microarrays. The basis for uBLU is a Bayesian model for the data samples which are represented as an additive mixture of random positive gene signatures, called factors, with random positive mixing coefficients, called factor scores, that specify the relative contribution of each signature to a specific sample. The particularity of the proposed method is that uBLU constrains the factor loadings to be non-negative and the factor scores to be probability distributions over the factors. Furthermore, it also provides estimates of the number of factors. A Gibbs sampling strategy is adopted here to generate random samples according to the posterior distribution of the factors, factor scores, and number of factors. These samples are then used to estimate all the unknown parameters. Results: Firstly, the proposed uBLU method is applied to several simulated datasets with known ground truth and compared with previous factor decomposition methods, such as principal component analysis (PCA), non negative matrix factorization (NMF), Bayesian factor regression modeling (BFRM), and the gradient-based algorithm for general matrix factorization (GB-GMF). Secondly, we illustrate the application of uBLU on a real time-evolving gene expression dataset from a recent viral challenge study in which individuals have been inoculated with influenza A/H3N2/Wisconsin. We show that the uBLU method significantly outperforms the other methods on the simulated and real data sets considered here. Conclusions: The results obtained on synthetic and real data illustrate the accuracy of the proposed uBLU method when compared to other factor decomposition methods from the literature (PCA, NMF, BFRM, and GB-GMF). The uBLU method identifies an inflammatory component closely associated with clinical symptom scores collected during the study. Using a constrained model allows recovery of all the inflammatory genes in a single factor

PKC Signaling Regulates Drug Resistance of the Fungal Pathogen Candida albicans via Circuitry Comprised of Mkc1, Calcineurin, and Hsp90

Fungal pathogens exploit diverse mechanisms to survive exposure to antifungal drugs. This poses concern given the limited number of clinically useful antifungals and the growing population of immunocompromised individuals vulnerable to life-threatening fungal infection. To identify molecules that abrogate resistance to the most widely deployed class of antifungals, the azoles, we conducted a screen of 1,280 pharmacologically active compounds. Three out of seven hits that abolished azole resistance of a resistant mutant of the model yeast Saccharomyces cerevisiae and a clinical isolate of the leading human fungal pathogen Candida albicans were inhibitors of protein kinase C (PKC), which regulates cell wall integrity during growth, morphogenesis, and response to cell wall stress. Pharmacological or genetic impairment of Pkc1 conferred hypersensitivity to multiple drugs that target synthesis of the key cell membrane sterol ergosterol, including azoles, allylamines, and morpholines. Pkc1 enabled survival of cell membrane stress at least in part via the mitogen activated protein kinase (MAPK) cascade in both species, though through distinct downstream effectors. Strikingly, inhibition of Pkc1 phenocopied inhibition of the molecular chaperone Hsp90 or its client protein calcineurin. PKC signaling was required for calcineurin activation in response to drug exposure in S. cerevisiae. In contrast, Pkc1 and calcineurin independently regulate drug resistance via a common target in C. albicans. We identified an additional level of regulatory control in the C. albicans circuitry linking PKC signaling, Hsp90, and calcineurin as genetic reduction of Hsp90 led to depletion of the terminal MAPK, Mkc1. Deletion of C. albicans PKC1 rendered fungistatic ergosterol biosynthesis inhibitors fungicidal and attenuated virulence in a murine model of systemic candidiasis. This work establishes a new role for PKC signaling in drug resistance, novel circuitry through which Hsp90 regulates drug resistance, and that targeting stress response signaling provides a promising strategy for treating life-threatening fungal infections

Expanding the Understanding of Biases in Development of Clinical-Grade Molecular Signatures: A Case Study in Acute Respiratory Viral Infections

The promise of modern personalized medicine is to use molecular and clinical information to better diagnose, manage, and treat disease, on an individual patient basis. These functions are predominantly enabled by molecular signatures, which are computational models for predicting phenotypes and other responses of interest from high-throughput assay data. Data-analytics is a central component of molecular signature development and can jeopardize the entire process if conducted incorrectly. While exploratory data analysis may tolerate suboptimal protocols, clinical-grade molecular signatures are subject to vastly stricter requirements. Closing the gap between standards for exploratory versus clinically successful molecular signatures entails a thorough understanding of possible biases in the data analysis phase and developing strategies to avoid them.Using a recently introduced data-analytic protocol as a case study, we provide an in-depth examination of the poorly studied biases of the data-analytic protocols related to signature multiplicity, biomarker redundancy, data preprocessing, and validation of signature reproducibility. The methodology and results presented in this work are aimed at expanding the understanding of these data-analytic biases that affect development of clinically robust molecular signatures.Several recommendations follow from the current study. First, all molecular signatures of a phenotype should be extracted to the extent possible, in order to provide comprehensive and accurate grounds for understanding disease pathogenesis. Second, redundant genes should generally be removed from final signatures to facilitate reproducibility and decrease manufacturing costs. Third, data preprocessing procedures should be designed so as not to bias biomarker selection. Finally, molecular signatures developed and applied on different phenotypes and populations of patients should be treated with great caution

Using gene expression profiles from peripheral blood to identify asymptomatic responses to acute respiratory viral infections

<p>Abstract</p> <p>Background</p> <p>A recent study reported that gene expression profiles from peripheral blood samples of healthy subjects prior to viral inoculation were indistinguishable from profiles of subjects who received viral challenge but remained asymptomatic and uninfected. If true, this implies that the host immune response does not have a molecular signature. Given the high sensitivity of microarray technology, we were intrigued by this result and hypothesize that it was an artifact of data analysis.</p> <p>Findings</p> <p>Using acute respiratory viral challenge microarray data, we developed a molecular signature that for the first time allowed for an accurate differentiation between uninfected subjects prior to viral inoculation and subjects who remained asymptomatic after the viral challenge.</p> <p>Conclusions</p> <p>Our findings suggest that molecular signatures can be used to characterize immune responses to viruses and may improve our understanding of susceptibility to viral infection with possible implications for vaccine development.</p

Multidrug-Resistant Tuberculosis Treatment Outcomes in Karakalpakstan, Uzbekistan: Treatment Complexity and XDR-TB among Treatment Failures

BACKGROUND: A pilot programme to treat multidrug-resistant TB (MDR-TB) was implemented in Karakalpakstan, Uzbekistan in 2003. This region has particularly high levels of MDR-TB, with 13% and 40% among new and previously treated cases, respectively. METHODOLOGY: This study describes the treatment process and outcomes for the first cohort of patients enrolled in the programme, between October 2003 and January 2005. Confirmed MDR-TB cases were treated with an individualised, second-line drug regimen based on drug susceptibility test results, while suspected MDR-TB cases were treated with a standardised regimen pending susceptibility results. PRINCIPAL FINDINGS: Of 108 MDR-TB patients, 87 were started on treatment during the study period. Of these, 33 (38%) were infected with strains resistant to at least one second-line drug at baseline, but none had initial ofloxacin resistance. Treatment was successful for 54 (62%) patients, with 13 (15%) dying during treatment, 12 (14%) defaulting and 8 (8%) failing treatment. Poor clinical condition and baseline second-line resistance contributed to treatment failure or death. Treatment regimens were changed in 71 (82%) patients due to severe adverse events or drug resistance. Adverse events were most commonly attributed to cycloserine, ethionamide and p-aminosalicylic acid. Extensively drug resistant TB (XDR-TB) was found among 4 of the 6 patients who failed treatment and were still alive in November 2006. CONCLUSIONS: While acceptable treatment success was achieved, the complexity of treatment and the development of XDR-TB among treatment failures are important issues to be addressed when considering scaling up MDR-TB treatment

Caveolin-1 protects B6129 mice against Helicobacter pylori gastritis.

Caveolin-1 (Cav1) is a scaffold protein and pathogen receptor in the mucosa of the gastrointestinal tract. Chronic infection of gastric epithelial cells by Helicobacter pylori (H. pylori) is a major risk factor for human gastric cancer (GC) where Cav1 is frequently down-regulated. However, the function of Cav1 in H. pylori infection and pathogenesis of GC remained unknown. We show here that Cav1-deficient mice, infected for 11 months with the CagA-delivery deficient H. pylori strain SS1, developed more severe gastritis and tissue damage, including loss of parietal cells and foveolar hyperplasia, and displayed lower colonisation of the gastric mucosa than wild-type B6129 littermates. Cav1-null mice showed enhanced infiltration of macrophages and B-cells and secretion of chemokines (RANTES) but had reduced levels of CD25+ regulatory T-cells. Cav1-deficient human GC cells (AGS), infected with the CagA-delivery proficient H. pylori strain G27, were more sensitive to CagA-related cytoskeletal stress morphologies ("humming bird") compared to AGS cells stably transfected with Cav1 (AGS/Cav1). Infection of AGS/Cav1 cells triggered the recruitment of p120 RhoGTPase-activating protein/deleted in liver cancer-1 (p120RhoGAP/DLC1) to Cav1 and counteracted CagA-induced cytoskeletal rearrangements. In human GC cell lines (MKN45, N87) and mouse stomach tissue, H. pylori down-regulated endogenous expression of Cav1 independently of CagA. Mechanistically, H. pylori activated sterol-responsive element-binding protein-1 (SREBP1) to repress transcription of the human Cav1 gene from sterol-responsive elements (SREs) in the proximal Cav1 promoter. These data suggested a protective role of Cav1 against H. pylori-induced inflammation and tissue damage. We propose that H. pylori exploits down-regulation of Cav1 to subvert the host's immune response and to promote signalling of its virulence factors in host cells

Aberrant Cell Cycle and Apoptotic Changes Characterise Severe Influenza A Infection – A Meta-Analysis of Genomic Signatures in Circulating Leukocytes

Influenza A infection is a global disease that has been responsible for four pandemics over the last one hundred years. However, it remains poorly understood as to why some infected individuals succumb to life threatening complications whilst others recover and are relatively unaffected. Using gene-expression analysis of circulating leukocytes, here we show that the progression towards severe influenza A infection is characterised by an abnormal transcriptional reprogramming of cell cycle and apoptosis pathways. In severely infected humans, leukocyte gene-expression profiles display opposing cell cycle activities; an increased aberrant DNA replication in the G1/S phase yet delayed progression in the G2/M phase. In mild infection, cell cycle perturbations are fewer and are integrated with an efficient apoptotic program. Importantly, the loss of integration between cell cycle perturbations and apoptosis marks the transition from a mild viral illness to a severe, life threatening infection. Our findings suggest that circulating immune cells may play a significant role in the evolution of the host response. Further study may reveal alternative host response factors previously unrecognized in the current disease model of influenza

A Simple Method for Combining Genetic Mapping Data from Multiple Crosses and Experimental Designs

Over the past decade many linkage studies have defined chromosomal intervals containing polymorphisms that modulate a variety of traits. Many phenotypes are now associated with enough mapping data that meta-analysis could help refine locations of known QTLs and detect many novel QTLs.We describe a simple approach to combining QTL mapping results for multiple studies and demonstrate its utility using two hippocampus weight loci. Using data taken from two populations, a recombinant inbred strain set and an advanced intercross population we demonstrate considerable improvements in significance and resolution for both loci. 1-LOD support intervals were improved 51% for Hipp1a and 37% for Hipp9a. We first generate locus-wise permuted P-values for association with the phenotype from multiple maps, which can be done using a permutation method appropriate to each population. These results are then assigned to defined physical positions by interpolation between markers with known physical and genetic positions. We then use Fisher's combination test to combine position-by-position probabilities among experiments. Finally, we calculate genome-wide combined P-values by generating locus-specific P-values for each permuted map for each experiment. These permuted maps are then sampled with replacement and combined. The distribution of best locus-specific P-values for each combined map is the null distribution of genome-wide adjusted P-values.Our approach is applicable to a wide variety of segregating and non-segregating mapping populations, facilitates rapid refinement of physical QTL position, is complementary to other QTL fine mapping methods, and provides an appropriate genome-wide criterion of significance for combined mapping results

Genetic susceptibility to aspergillosis in allogeneic stem-cell transplantation

Invasive aspergillosis (IA) is a major threat to positive outcomes for allogeneic stem-cell transplantation (allo-SCT) patients. Despite presenting similar degrees of immunosuppression, not all individuals at-risk ultimately develop infection. Therefore, the traditional view of neutropenia as a key risk factor for aspergillosis needs to be accommodated within new conceptual advances on host immunity and its relationship to infection. Polymorphisms in innate immune genes, such as those encoding TLRs, cytokines and cytokine receptors, have recently been associated with susceptibility to IA in allo-SCT recipients. This suggests that understanding host-pathogen interactions at the level of host genetic susceptibility will allow the formulation of new targeted and patient-tailored antifungal therapeutics, including improved donor screening.Fundação para a Ciência e a Tecnologia (FCT) - SFRH/BD/65962/2009, SFRH/BPD/46292/2008Specific Targeted Research Projects MANASP (LSHE-CT-2006), contract number 037899 (FP6), Italian Project PRIN2007KLCKP8_004

Integrating Factor Analysis and a Transgenic Mouse Model to Reveal a Peripheral Blood Predictor of Breast Tumors

Abstract Background Transgenic mouse tumor models have the advantage of facilitating controlled in vivo oncogenic perturbations in a common genetic background. This provides an idealized context for generating transcriptome-based diagnostic models while minimizing the inherent noisiness of high-throughput technologies. However, the question remains whether models developed in such a setting are suitable prototypes for useful human diagnostics. We show that latent factor modeling of the peripheral blood transcriptome in a mouse model of breast cancer provides the basis for using computational methods to link a mouse model to a prototype human diagnostic based on a common underlying biological response to the presence of a tumor. Methods We used gene expression data from mouse peripheral blood cell (PBC) samples to identify significantly differentially expressed genes using supervised classification and sparse ANOVA. We employed these transcriptome data as the starting point for developing a breast tumor predictor from human peripheral blood mononuclear cells (PBMCs) by using a factor modeling approach. Results The predictor distinguished breast cancer patients from healthy individuals in a cohort of patients independent from that used to build the factors and train the model with 89% sensitivity, 100% specificity and an area under the curve (AUC) of 0.97 using Youden's J-statistic to objectively select the model's classification threshold. Both permutation testing of the model and evaluating the model strategy by swapping the training and validation sets highlight its stability. Conclusions We describe a human breast tumor predictor based on the gene expression of mouse PBCs. This strategy overcomes many of the limitations of earlier studies by using the model system to reduce noise and identify transcripts associated with the presence of a breast tumor over other potentially confounding factors. Our results serve as a proof-of-concept for using an animal model to develop a blood-based diagnostic, and it establishes an experimental framework for identifying predictors of solid tumors, not only in the context of breast cancer, but also in other types of cancer.</p