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

Gene expression profiles link respiratory viral infection, platelet response to aspirin, and acute myocardial infarction

Background Influenza infection is associated with myocardial infarction (MI), suggesting that respiratory viral infection may induce biologic pathways that contribute to MI. We tested the hypotheses that 1) a validated blood gene expression signature of respiratory viral infection (viral GES) was associated with MI and 2) respiratory viral exposure changes levels of a validated platelet gene expression signature (platelet GES) of platelet function in response to aspirin that is associated with MI. Methods A previously defined viral GES was projected into blood RNA data from 594 patients undergoing elective cardiac catheterization and used to classify patients as having evidence of viral infection or not and tested for association with acute MI using logistic regression. A previously defined platelet GES was projected into blood RNA data from 81 healthy subjects before and after exposure to four respiratory viruses: Respiratory Syncytial Virus (RSV) (n=20), Human Rhinovirus (HRV) (n=20), Influenza A virus subtype H1N1 (H1N1) (n=24), Influenza A Virus subtype H3N2 (H3N2) (n=17). We tested for the change in platelet GES with viral exposure using linear mixed-effects regression and by symptom status. Results In the catheterization cohort, 32 patients had evidence of viral infection based upon the viral GES, of which 25% (8/32) had MI versus 12.2%(69/567) among those without evidence of viral infection (OR 2.3; CI [1.03-5.5], p=0.04). In the infection cohorts, only H1N1 exposure increased platelet GES over time (time course p-value = 1e-04). Conclusions A viral GES of non-specific, respiratory viral infection was associated with acute MI; 18% of the top 49 genes in the viral GES are involved with hemostasis and/or platelet aggregation. Separately, H1N1 exposure, but not exposure to other respiratory viruses, increased a platelet GES previously shown to be associated with MI. Together, these results highlight specific genes and pathways that link viral infection, platelet activation, and MI especially in the case of H1N1 influenza infection

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

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

Systemic Signature of the Lung Response to Respiratory Syncytial Virus Infection

Respiratory Syncytial Virus is a frequent cause of severe bronchiolitis in children. To improve our understanding of systemic host responses to RSV, we compared BALB/c mouse gene expression responses at day 1, 2, and 5 during primary RSV infection in lung, bronchial lymph nodes, and blood. We identified a set of 53 interferon-associated and innate immunity genes that give correlated responses in all three murine tissues. Additionally, we identified blood gene signatures that are indicative of acute infection, secondary immune response, and vaccine-enhanced disease, respectively. Eosinophil-associated ribonucleases were characteristic for the vaccine-enhanced disease blood signature. These results indicate that it may be possible to distinguish protective and unfavorable patient lung responses via blood diagnostics

Uncovering the Genetic Landscape for Multiple Sleep-Wake Traits

Despite decades of research in defining sleep-wake properties in mammals, little is known about the nature or identity of genes that regulate sleep, a fundamental behaviour that in humans occupies about one-third of the entire lifespan. While genome-wide association studies in humans and quantitative trait loci (QTL) analyses in mice have identified candidate genes for an increasing number of complex traits and genetic diseases, the resources and time-consuming process necessary for obtaining detailed quantitative data have made sleep seemingly intractable to similar large-scale genomic approaches. Here we describe analysis of 20 sleep-wake traits from 269 mice from a genetically segregating population that reveals 52 significant QTL representing a minimum of 20 genomic loci. While many (28) QTL affected a particular sleep-wake trait (e.g., amount of wake) across the full 24-hr day, other loci only affected a trait in the light or dark period while some loci had opposite effects on the trait during the light vs. dark. Analysis of a dataset for multiple sleep-wake traits led to previously undetected interactions (including the differential genetic control of number and duration of REM bouts), as well as possible shared genetic regulatory mechanisms for seemingly different unrelated sleep-wake traits (e.g., number of arousals and REM latency). Construction of a Bayesian network for sleep-wake traits and loci led to the identification of sub-networks of linkage not detectable in smaller data sets or limited single-trait analyses. For example, the network analyses revealed a novel chain of causal relationships between the chromosome 17@29cM QTL, total amount of wake, and duration of wake bouts in both light and dark periods that implies a mechanism whereby overall sleep need, mediated by this locus, in turn determines the length of each wake bout. Taken together, the present results reveal a complex genetic landscape underlying multiple sleep-wake traits and emphasize the need for a systems biology approach for elucidating the full extent of the genetic regulatory mechanisms of this complex and universal behavior

Two Genes on A/J Chromosome 18 Are Associated with Susceptibility to Staphylococcus aureus Infection by Combined Microarray and QTL Analyses

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N2 backcross mice (F1 [C18A]×C57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus–challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 β and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies

Temporal Dynamics of Host Molecular Responses Differentiate Symptomatic and Asymptomatic Influenza A Infection

Exposure to influenza viruses is necessary, but not sufficient, for healthy human hosts to develop symptomatic illness. The host response is an important determinant of disease progression. In order to delineate host molecular responses that differentiate symptomatic and asymptomatic Influenza A infection, we inoculated 17 healthy adults with live influenza (H3N2/Wisconsin) and examined changes in host peripheral blood gene expression at 16 timepoints over 132 hours. Here we present distinct transcriptional dynamics of host responses unique to asymptomatic and symptomatic infections. We show that symptomatic hosts invoke, simultaneously, multiple pattern recognition receptors-mediated antiviral and inflammatory responses that may relate to virus-induced oxidative stress. In contrast, asymptomatic subjects tightly regulate these responses and exhibit elevated expression of genes that function in antioxidant responses and cell-mediated responses. We reveal an ab initio molecular signature that strongly correlates to symptomatic clinical disease and biomarkers whose expression patterns best discriminate early from late phases of infection. Our results establish a temporal pattern of host molecular responses that differentiates symptomatic from asymptomatic infections and reveals an asymptomatic host-unique non-passive response signature, suggesting novel putative molecular targets for both prognostic assessment and ameliorative therapeutic intervention in seasonal and pandemic influenza