Gene Sets Net Correlations Analysis (GSNCA): a multivariate differential coexpression test for gene sets

Motivation: To date, Gene Set Analysis (GSA) approaches primarily focus on identifying differentially expressed gene sets (pathways). Methods for identifying differentially coexpressed pathways also exist but are mostly based on aggregated pairwise correlations, or other pairwise measures of coexpression. Instead, we propose Gene Sets Net Correlations Analysis (GSNCA), a multivariate differential coexpression test that accounts for the complete correlation structure between genes.

Results: In GSNCA, weight factors are assigned to genes in proportion to the genes' cross-correlations (intergene correlations). The problem of finding the weight vectors is formulated as an eigenvector problem with a unique solution. GSNCA tests the null hypothesis that for a gene set there is no difference in the weight vectors of the genes between two conditions. In simulation studies and the analyses of experimental data, we demonstrate that GSNCA, indeed, captures changes in the structure of genes' cross-correlations rather than differences in the averaged pairwise correlations. Thus, GSNCA infers differences in coexpression networks, however, bypassing method-dependent steps of network inference. As an additional result from GSNCA, we define hub genes as genes with the largest weights and show that these genes correspond frequently to major and specific pathway regulators, as well as to genes that are most affected by the biological difference between two conditions. In summary, GSNCA is a new approach for the analysis of differentially coexpressed pathways that also evaluates the importance of the genes in the pathways, thus providing unique information that may result in the generation of novel biological hypotheses.

SOCS2 Can Enhance Interleukin-2 (IL-2) and IL-3 Signalling by Accelerating SOCS3 Degradation

Cytokine responses can be regulated by a family of proteins termed suppressors of cytokine signaling (SOCS) which can inhibit the JAK/STAT pathway in a classical negative-feedback manner. While the SOCS are thought to target signaling intermediates for degradation, relatively little is known about how their turnover is regulated. Unlike other SOCS family members, we find that SOCS2 can enhance interleukin-2 (IL-2)- and IL-3-induced STAT phosphorylation following and potentiate proliferation in response to cytokine stimulation. As a clear mechanism for these effects, we demonstrate that expression of SOCS2 results in marked proteasome-dependent reduction of SOCS3 and SOCS1 protein expression. Furthermore, we provide evidence that this degradation is dependent on the presence of an intact SOCS box and that the loss of SOCS3 is enhanced by coexpression of elongin B/C. This suggests that SOCS2 can bind to SOCS3 and elongin B/C to form an E3 ligase complex resulting in the degradation of SOCS3. Therefore, SOCS2 can enhance cytokine responses by accelerating proteasome-dependent turnover of SOCS3, suggesting a mechanism for the gigantism observed in SOCS2 transgenic mice.

Pim2 cooperates with PML-RAR alpha to induce acute myeloid leukemia in a bone marrow transplantation model

Although the potential role of Pim2 as a cooperative oncogene has been well described in lymphoma, its role in leukemia has remained largely unexplored. Here we show that high expression of Pim2 is observed in patients with acute promyelocytic leukemia (APL). To further characterize the cooperative role of Pim2 with promyelocytic leukemia/retinoic acid receptor alpha (PML/RAR alpha), we used a well-established PML-RAR alpha (PR alpha) mouse model. Pim2 coexpression in PR alpha-positive hematopoietic progenitor cells (HPCs) induces leukemia in recipient mice after a short latency. Pim2-PR alpha cells were able to repopulate mice in serial transplantations and to induce disease in all recipients. Neither Pim2 nor PR alpha alone was sufficient to induce leukemia upon transplantation in this model. The disease induced by Pim2 overexpression in PR alpha cells contained a slightly higher fraction of immature myeloid cells, compared with the previously described APL disease induced by PR alpha. However, it also clearly resembled an APL-like phenotype and showed signs of differentiation upon all-trans retinoic acid (ATRA) treatment in vitro. These results support the hypothesis that Pim2, which is also a known target of Flt3-ITD (another gene that cooperates with PML-RAR alpha), cooperates with PR alpha to induce APL-like disease. (Blood. 2010; 115(22): 4507-4516)

Regarding “Co-expression of SNAIL and TWIST determines prognosis in estrogen receptor-positive early breast cancer patients”

We read the interesting research article published by van Nes et al. [1], which described the use of Snail and TWIST together in the prognosis of breast cancer, and in particular in estrogen receptor (ER)-positive breast cancer patients.

Functional analysis of the glycero-manno-heptose 7-phosphate kinase domain from the bifunctional HldE protein, which is involved in ADP-L-glycero-D-manno-heptose biosynthesis

The core oligosaccharide component of the lipopolysaccharide can be subdivided into inner and outer core regions. In Escherichia coli, the inner core consists of two 3-deoxy-d-manno-octulosonic acid and three glycero-manno-heptose residues. The HldE protein participates in the biosynthesis of ADP-glycero-manno-heptose precursors used in the assembly of the inner core. HldE comprises two functional domains: an N-terminal region with homology to the ribokinase superfamily (HldE1 domain) and a C-terminal region with homology to the cytidylyltransferase superfamily (HldE2 domain). We have employed the structure of the E. coli ribokinase as a template to model the HldE1 domain and predict critical amino acids required for enzyme activity. Mutation of these residues renders the protein inactive as determined in vivo by functional complementation analysis. However, these mutations did not affect the secondary or tertiary structure of purified HldE1, as judged by fluorescence spectroscopy and circular dichroism. Furthermore, in vivo coexpression of wild-type, chromosomally encoded HldE and mutant HldE1 proteins with amino acid substitutions in the predicted ATP binding site caused a dominant negative phenotype as revealed by increased bacterial sensitivity to novobiocin. Copurification experiments demonstrated that HldE and HldE1 form a complex in vivo. Gel filtration chromatography resulted in the detection of a dimer as the predominant form of the native HldE1 protein. Altogether, our data support the notions that the HldE functional unit is a dimer and that structural components present in each HldE1 monomer are required for enzymatic activity.

Lacking cytokine production in ES cells and ES-cell-derived vascular cells stimulated by TNF-alpha is rescued by HDAC inhibitor trichostatin A

Inflammation and TNF-alpha signaling play a central role in most of the pathological conditions where cell transplantation could be applied. As shown by initial experiments, embryonic stem (ES) cells and ES-cell derived vascular cells express very low levels of TNF-alpha receptor I (TNFRp55) and thus do not induce cytokine expression in response to TNF-alpha stimulation. Transient transfection analysis of wild-type or deletion variants of the TNFRp55 gene promoter showed a strong activity for a 250-bp fragment in the upstream region of the gene. This activity was abolished by mutations targeting the Sp1/Sp3 or AP1 binding sites. Moreover, treatment with trichostatin A (TSA) led to a pronounced increase in TNFRp55 mRNA and promoter activity. Overexpression of Sp1 or c-fos further enhanced the TSA-induced luciferase activity, and this response was attenuated by Sp3 or c-jun coexpression. Additional experiments revealed that TSA did not affect the Sp1/Sp3 ratio but caused transcriptional activation of the c-fos gene. Thus, we provide the first evidence that ES and ES-cell-derived vascular cells lack cytokine expression in response to TNF-alpha stimulation due to low levels of c-fos and transcriptional activation of Sp1 that can be regulated by inhibition of histone deacetylase activity.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Nuclear LYRIC/AEG-1 interacts with PLZF and relieves PLZF-mediated repression

LYRIC/AEG-1 and its altered expression have been linked to carcinogenesis in prostate, brain and melanoma as well as promoting chemoresistance and metastasis in breast cancer. LYRIC/AEG-1 function remains unclear, although LYRIC/AEG-1 is activated by oncogenic HA-RAS, through binding of c-myc to its promoter, which in turn regulates the key components of the PI3-kinase and nuclear factor-kappaB pathways. We have identified the transcriptional repressor PLZF as an interacting protein of LYRIC/AEG through a yeast two-hybrid screen. PLZF regulates the expression of genes involved in cell growth and apoptosis including c-myc. Coexpression of LYRIC/AEG-1 with PLZF leads to a reduction in PLZF-mediated repression by reducing PLZF binding to promoters. We have confirmed that nuclear LYRIC/AEG-1 and PLZF interact in mammalian cells via the N- and C termini of LYRIC/AEG-1 and a region C terminal to the RD2 domain of PLZF. Both proteins colocalize to nuclear bodies containing histone deacetylases, which are known to promote PLZF-mediated repression. Our data suggest one mechanism for cells with altered LYRIC/AEG-1 expression to evade apoptosis and increase cell growth during tumourigenesis through the regulation of PLZF repression.

Co-expression of NPY and VIP in human dental pulp

Introduction: The regulation of pulpal haemodynamics in health and disease involves sympathetic and parasympathetic mechanisms in which both neuropeptide Y (NPY; a sympathetic vasoconstrictor) and vasoactive intestinal polypeptide (VIP; a parasympathetic vasodilator) may play potential pathophysiological roles. We have previously investigated the levels of NPY or VIP present in human dental pulp tissue and shown that their expression is up-regulated in caries induced pulpal inflammation. Objectives: The aim of this study was to investigate the potential correlation between NPY and VIP levels measured in the same dental pulp samples using radioimmunoassay (RIA). Methods: Pulp tissue was obtained from extracted teeth, classified as follows; healthy (n=22), moderately carious (n=20) and grossly carious (n=26). Samples were processed for RIA by boiling in acetic acid as previously described. The levels of NPY and VIP, measured by RIA, were expressed as ng/gram of pulp tissue. The nature of the relationship between NPY and VIP levels in human pulp tissue was tested by calculating Pearson's product moment correlation coefficient using the linear regression test. Results: Calculation of Pearson product moment correlation coefficient showed a significant negative correlation between NPY and VIP levels in pulp tissue samples from non-carious teeth (p = 0.02, r = -.48). This negative correlation in non-carious teeth changed to a significant positive correlation in carious teeth when the levels of NPY and VIP were compared (p = 0.03, r= 0.311). Conclusions: In non-carious teeth, the negative correlation between NPY and VIP levels is in keeping with the previously described modulatory influence of cholinergic nerves on sympathetic function which may be perturbed as caries develops.

Data for the co-expression and purification of human recombinant CaMKK2 in complex with calmodulin in Escherichia coli

Calcium/calmodulin-dependent kinase kinase 2 (CaMKK2) has been implicated in a range of conditions and pathologies from prostate to hepatic cancer. Here, we describe the expression in Escherichia coli and the purification protocol for the following constructs: full-length CaMKK2 in complex with CaM, CaMKK2 'apo', CaMKK2 (165-501) in complex with CaM, and the CaMKK2 F267G mutant. The protocols described have been optimized for maximum yield and purity with minimal purification steps required and the proteins subsequently used to develop a fluorescence-based assay for drug binding to the kinase, "Using the fluorescent properties of STO-609 as a tool to assist structure-function analyses of recombinant CaMKK2"