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.

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.

Using gene-network landscape to dissect genotype effects of TCF7L2 genetic variant on diabetes and cardiovascular risk

Vaquero AR, Ferreira NE, Omae SV, Rodrigues MV, Teixeira SK, Krieger JE, Pereira AC. Using gene-network landscape to dissect genotype effects of TCF7L2 genetic variant on diabetes and cardiovascular risk. Physiol Genomics 44: 903-914, 2012. First published August 7, 2012; doi:10.1152/physiolgenomics.00030.2012.-The single nucleotide polymorphism (SNP) within the TCF7L2 gene, rs7903146, is, to date, the most significant genetic marker associated with Type 2 diabetes mellitus (T2DM) risk. Nonetheless, its functional role in disease pathology is poorly understood. The aim of the present study was to investigate, in vascular smooth muscle cells from 92 patients undergoing aortocoronary bypass surgery, the contribution of this SNP in T2DM using expression levels and expression correlation comparison approaches, which were visually represented as gene interaction networks. Initially, the expression levels of 41 genes (seven TCF7L2 splice forms and 40 other T2DM relevant genes) were compared between rs7903146 wild-type (CC) and T2DM-risk (CT + TT) genotype groups. Next, we compared the expression correlation patterns of these 41 genes between groups to observe if the relationships between genes were different. Five TCF7L2 splice forms and nine genes showed significant expression differences between groups. RXR alpha gene was pinpointed as showing the most different expression correlation pattern with other genes. Therefore, T2DM risk alleles appear to be influencing TCF7L2 splice form's expression in vascular smooth muscle cells, and RXR alpha gene is pointed out as a treatment target candidate for risk reduction in individuals with high risk of developing T2DM, especially individuals harboring TCF7L2 risk genotypes.

Differential cell type-specific transcriptional regulation of the CYP1A1 gene

Cytochrome P450 1A1 (CYP1A1) monooxygenase plays an important role in the metabolism of environmental pollutants such as polycyclic aromatic hydrocarbons (PAHs) and halogenated polycyclic aromatic hydrocarbons (HAHs). Oxidation of these compounds converts them to the metabolites that subsequently can be conjugated to hydrophilic endogenous entities e.g. glutathione. Derivates generated in this way are water soluble and can be excreted in bile or urine, which is a defense mechanism. Besides detoxification, metabolism by CYP1A1 may lead to deleterious effects since the highly reactive intermediate metabolites are able to react with DNA and thus cause mutagenic effects, as it is in the case of benzo(a) pyrene (B[a]P). CYP1A1 is normally not expressed or expressed at a very low level in the cells but it is inducible by many PAHs and HAHs e.g. by B[a]P or 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Transcriptional activation of the CYP1A1 gene is mediated by aryl hydrocarbon receptor (AHR), a basic-helix-loop-helix (bHLH) transcription factor. In the absence of a ligand AHR stays predominantly in the cytoplasm. Ligand binding causes translocation of AHR to the nuclear compartment, its heterodimerization with another bHLH protein, the aryl hydrocarbon nuclear translocator (ARNT) and binding of the AHR/ARNT heterodimer to a DNA motif designated dioxin responsive element (DRE). This process leads to the transcriptional activation of the responsive genes containing DREs in their regulatory regions, e.g. that coding for CYP1A1. TCDD is the most potent known agonist of AHR. Since it is not metabolized by the activated enzymes, exposure to this compound leads to a persisting activation of AHR resulting in diverse toxic effects in the organism. To enlighten the molecular mechanisms that mediate the toxicity of xenobiotics like TCDD and related compounds, the AHR-dependent regulation of the CYP1A1 gene was investigated in two cell lines: human cervix carcinoma (HeLa) and mouse hepatoma (Hepa). Study of AHR activation and its consequence concerning expression of the CYP1A1 enzyme confirmed the TCDD-dependent formation of the AHR/ARNT complex on DRE leading to an increase of the CYP1A1 transcription in Hepa cells. In contrast, in HeLa cells formation of the AHR/ARNT heterodimer and binding of a protein complex containing AHR and ARNT to DRE occurred naturally in the absence of TCDD. Moreover, treatment with TCDD did not affect the AHR/ARNT dimer formation and binding of these proteins to DRE in these cells. Even though the constitutive complex on DRE exists in HeLa, transcription of the CYP1A1 gene was not increased. Furthermore, the CYP1A1 level in HeLa cells remained unchanged in the presence of TCDD suggesting repressional mechanism of the AHR complex function which may hinder the TCDD-dependent mechanisms in these cells. Similar to the native, the mouse CYP1A1-driven reporter constructs containing different regulatory elements were not inducible by TCDD in HeLa cells, which supported a presence of cell type specific trans-acting factor in HeLa cells able to repress both the native CYP1A1 and CYP1A1-driven reporter genes rather than species specific differences between CYP1A1 genes of human and rodent origin. The different regulation of the AHR-mediated transcription of CYP1A1 gene in Hepa and HeLa cells was further explored in order to elucidate two aspects of the AHR function: (I) mechanism involved in the activation of AHR in the absence of exogenous ligand and (II) factor that repress function of the exogenous ligand-independent AHR/ARNT complex. Since preliminary studies revealed that the activation of PKA causes an activation of AHR in Hepa cells in the absence of TCDD, the PKA-dependent signalling pathway was the proposed endogenous mechanism leading to the TCDD-independent activation of AHR in HeLa cells. Activation of PKA by forskolin or db-cAMP as well as inhibition of the kinase by H89 in both HeLa and Hepa cells did not lead to alterations in the AHR interaction with ARNT in the absence of TCDD and had no effect on binding of these proteins to DRE. Moreover, the modulators of PKA did not influence the CYP1A1 activity in these cells in the presence and in the absence of TCDD. Thus, an involvement of PKA in the regulation of the CYP1A1 Gen in HeLa cells was not evaluated in the course of this study. Repression of genes by transcription factors bound to their responsive elements in the absence of ligands has been described for nuclear receptors. These receptors interact with protein complex containing histone deacetylase (HDAC), enzyme responsible for the repressional effect. Thus, a participation of histone deacetylase in the transcriptional modulation of CYP1A1 gene by the constitutively DNA-bound AHR/ARNT complex was supposed. Inhibition of the HDAC activity by trichostatin A (TSA) or sodium butyrate (NaBu) led to an increase of the CYP1A1 transcription in the presence but not in the absence of TCDD in Hepa and HeLa cells. Since amount of the AHR and ARNT proteins remained unchanged upon treatment of the cells with TSA or NaBu, the transcriptional upregulation of CYP1A1 gene was not due to an increased expression of the regulatory proteins. These findings strongly suggest an involvement of HDAC in the repression of the CYP1A1 gene. Similar to the native human CYP1A1 also the mouse CYP1A1-driven reporter gene transfected into HeLa cells was repressed by histone deacetylase since the presence of TSA or NaBu led to an increase in the reporter activity. Induction of reporter gene did not require a presence of the promoter or negative regulatory regions of the CYP1A1 gene. A promoter-distal fragment containing three DREs together with surrounding sequences was sufficient to mediate the effects of the HDAC inhibitors suggesting that the AHR/ARNT binding to its specific DNA recognition site may be important for the CYP1A1 repression. Histone deacetylase is recruited to the specific genes by corepressors, proteins that bind to the transcription factors and interact with other members of the HDAC complex. Western blot analyses revealed a presence of HDAC1 and the corepressors mSin3A (mammalian homolog of yeast Sin3) and SMRT (silencing mediator for retinoid and thyroid hormone receptor) in both cell types, while the corepressor NCoR (nuclear receptor corepressor) was expressed exclusively in HeLa cells. Thus the high inducibility of CYP1A1 in Hepa cells may be due to the absence of NCoR in these cells in contrast to the non-responsive HeLa cells, where the presence of NCoR would support repression of the gene by histone deacetylase. This hypothesis was verified in reporter gene experiments where expression constructs coding for the particular members of the HDAC complex were cotransfected in Hepa cells together with the TCDD-inducible reporter constructs containing the CYP1A1 regulatory sequences. An overexpression of NCoR however did not decrease but instead led to a slight increase of the reporter gene activity in the cells. The expected inhibition was observed solely in the case of SMRT that slightly reduced constitutive and TCDD-induced reporter gene activity. A simultaneous expression of NCoR and SMRT shown no further effects and coexpression of HDAC1 with the two corepressors did not alter this situation. Thus, additional factors that are likely involved in the repression of CYP1A1 gene by HDAC complex remained to be identified. Taking together, characterisation of an exogenous ligand independent AHR/ARNT complex on DRE in HeLa cells that repress transcription of the CYP1A1 gene creates a model system enabling investigation of endogenous processes involved in the regulation of AHR function. This study implicates HDAC-mediated repression of CYP1A1 gene that contributes to the xenobiotic-induced expression in a tissue specific manner. Elucidation of these processes gains an insight into mechanisms leading to deleterious effects of TCDD and related compounds.

Glycogen-rich pleomorphic xanthoastrocytoma with clear-cell features: confirmatory report of a rare variant with implications for differential diagnosis

Central nervous system space-occupying lesions with clear-cell features encompass a nosologically heterogeneous array, ranging from reactive histiocytic proliferations to neuroepithelial or meningothelial neoplasms of various grades and to metastases. In the face of such differential diagnostic breadth, recognizing cytoplasmic lucency as part of the morphological spectrum of some low grade gliomas will directly have an impact on patient care. We describe a prevailing clear-cell change in an epileptogenic left temporal pleomorphic xanthoastrocytoma surgically resected from a 36-year-old man. Mostly subarachnoid and focally calcified, the tumor was composed of fascicles of moderately atypical spindle cells with optically lucent cytoplasm that tended to intermingle with a desmoplastic mesh of reticulin fibers. Immunohistochemically, coexpression of S100 protein, vimentin, GFAP, and CD34 was noted. Conversely, neither punctate staining for EMA nor positivity for CD68 was seen. Mitotic activity was absent, and the MIB1 labeling index was 2-3% on average. Diastase-sensitive PAS-positive granula indicated clear-cell change to proceed from glycogen storage. Electron microscopy showed tumor cell cytoplasm to be largely obliterated by non-lysosomal-bound pools of glycogen, while hardly any fat vacuole was encountered. Neither ependymal-derived organelles nor annular lamellae suggesting oligodendroglial differentiation were detected. The latter differential diagnosis was further invalidated by lack of codeletion of chromosomal regions 1p36 and 19q13 on molecular genetic testing. By significantly interfering with pattern recognition as an implicit approach in histopathology, clear-cell change in pleomorphic xanthoastrocytoma is likely to suspend its status as a "classic", and to prompt more deductive differential diagnostic strategies to exclude look-alikes, especially clear-cell ependymoma and oligodendroglioma.

A Petrov-Galerkin method for a singularly perturbed ordinary differential equation with non-smooth data

In this paper, a singularly perturbed ordinary differential equation with non-smooth data is considered. The numerical method is generated by means of a Petrov-Galerkin finite element method with the piecewise-exponential test function and the piecewise-linear trial function. At the discontinuous point of the coefficient, a special technique is used. The method is shown to be first-order accurate and singular perturbation parameter uniform convergence. Finally, numerical results are presented, which are in agreement with theoretical results.

Some Remarks on the Approximation of Transitional Density in Stochastic Differential Equations

Aijt-Sahalia (2002) introduced a method to estimate transitional probability densities of di®usion processes by means of Hermite expansions with coe±cients determined by means of Taylor series. This note describes a numerical procedure to ¯nd these coe±cients based on the calculation of moments. One advantage of this procedure is that it can be used e®ectively when the mathematical operations required to ¯nd closed-form expressions for these coe±cients are otherwise infeasible.

A broadcast authentication and integrity augmentation for trusted differential GPS in marine navigation

In this paper we propose an efficient authentication and integrity scheme to support DGPS corrections using the RTCM protocol, such that the identified vulnerabilities in DGPS are mitigated. The proposed scheme is based on the TESLA broadcast protocol with modifications that make it suitable for the bandwidth and processor constrained environment of marine DGPS.

Organizational ambidexterity and corporate entrepreneurship : the differential effects on venturing, innovation and renewal processes

Most corporate entrepreneurship studies have focused on either innovation, venturing or strategic renewal making comparison between the antecedents of all three aspects of corporate entrepreneurship difficult. Moreover, studies on corporate entrepreneurship hardly address organizational antecedents, while simultaneously managing and organizing CE and mainstream activities has been seen as a major challenge for incumbent firms. Our findings show that organizational ambidexterity has strong and differential effects on venturing, innovation and renewal. We find, for example, that innovation is affected by horizontal integration, while strategic renewal is significantly influenced by integration on top management team level.