Finding consistent disease subnetworks across microarray datasets

<p>Abstract</p> <p>Background</p> <p>While contemporary methods of microarray analysis are excellent tools for studying individual microarray datasets, they have a tendency to produce different results from different datasets of the same disease. We aim to solve this reproducibility problem by introducing a technique (SNet). SNet provides both quantitative and descriptive analysis of microarray datasets by identifying specific connected portions of pathways that are significant. We term such portions within pathways as “subnetworks”.</p> <p>Results</p> <p>We tested SNet on independent datasets of several diseases, including childhood ALL, DMD and lung cancer. For each of these diseases, we obtained two independent microarray datasets produced by distinct labs on distinct platforms. In each case, our technique consistently produced almost the same list of significant nontrivial subnetworks from two independent sets of microarray data. The gene-level agreement of these significant subnetworks was between 51.18% to 93.01%. In contrast, when the same pairs of microarray datasets were analysed using GSEA, t-test and SAM, this percentage fell between 2.38% to 28.90% for GSEA, 49.60% tp 73.01% for t-test, and 49.96% to 81.25% for SAM. Furthermore, the genes selected using these existing methods did not form subnetworks of substantial size. Thus it is more probable that the subnetworks selected by our technique can provide the researcher with more descriptive information on the portions of the pathway actually affected by the disease.</p> <p>Conclusions</p> <p>These results clearly demonstrate that our technique generates significant subnetworks and genes that are more consistent and reproducible across datasets compared to the other popular methods available (GSEA, t-test and SAM). The large size of subnetworks which we generate indicates that they are generally more biologically significant (less likely to be spurious). In addition, we have chosen two sample subnetworks and validated them with references from biological literature. This shows that our algorithm is capable of generating descriptive biologically conclusions.</p

RhoD Inhibits RhoC-ROCK-Dependent Cell Contraction via PAK6.

RhoA-mediated regulation of myosin-II activity in the actin cortex controls the ability of cells to contract and bleb during a variety of cellular processes, including cell migration and division. Cell contraction and blebbing also frequently occur as part of the cytopathic effect seen during many different viral infections. We now demonstrate that the vaccinia virus protein F11, which localizes to the plasma membrane, is required for ROCK-mediated cell contraction from 2 hr post infection. Curiously, F11-induced cell contraction is dependent on RhoC and not RhoA signaling to ROCK. Moreover, RhoC-driven cell contraction depends on the upstream inhibition of RhoD signaling by F11. This inhibition prevents RhoD from regulating its downstream effector Pak6, alleviating the suppression of RhoC by the kinase. Our observations with vaccinia have now demonstrated that RhoD recruits Pak6 to the plasma membrane to antagonize RhoC signaling during cell contraction and blebbing

The Minimal Autoinhibited Unit of the Guanine Nucleotide Exchange Factor Intersectin

Intersectin-1L is a member of the Dbl homology (DH) domain guanine nucleotide exchange factors (GEF) which control Rho-family GTPase signaling. Intersectin-1L is a GEF that is specific for Cdc42. It plays an important role in endocytosis, and is regulated by several partners including the actin regulator N-WASP. Intact intersectin-1L shows low Cdc42 exchange activity, although the isolated catalytic DH domain shows high activity. This finding suggests that the molecule is autoinhibited. To investigate the mechanism of autoinhibition we have constructed a series of domain deletions. We find that the five SH3 domains of intersectin are important for autoinhibition, with the fifth domain (SH3(E)) being sufficient for the bulk of the autoinhibitory effect. This SH3 domain appears to primarily interact with the DH domain. We have determined the crystal structure of the SH3(E)-DH domain construct, which shows a domain swapped arrangement in which the SH3 from one monomer interacts with the DH domain of the other monomer. Analytical ultracentrifugation and gel filtration, however, show that under biochemical concentrations, the construct is fully monomeric. Thus we propose that the actual autoinhibited structure contains the related intramolecular SH3(E)-DH interaction. We propose a model in which this intramolecular interaction may block or distort the GTPase binding region of the DH domain

The Retrohoming of Linear Group II Intron RNAs in Drosophila melanogaster Occurs by Both DNA Ligase 4–Dependent and –Independent Mechanisms

Mobile group II introns are bacterial retrotransposons that are thought to have invaded early eukaryotes and evolved into introns and retroelements in higher organisms. In bacteria, group II introns typically retrohome via full reverse splicing of an excised intron lariat RNA into a DNA site, where it is reverse transcribed by the intron-encoded protein. Recently, we showed that linear group II intron RNAs, which can result from hydrolytic splicing or debranching of lariat RNAs, can retrohome in eukaryotes by performing only the first step of reverse splicing, ligating their 3′ end to the downstream DNA exon. Reverse transcription then yields an intron cDNA, whose free end is linked to the upstream DNA exon by an error-prone process that yields junctions similar to those formed by non-homologous end joining (NHEJ). Here, by using Drosophila melanogaster NHEJ mutants, we show that linear intron RNA retrohoming occurs by major Lig4-dependent and minor Lig4-independent mechanisms, which appear to be related to classical and alternate NHEJ, respectively. The DNA repair polymerase θ plays a crucial role in both pathways. Surprisingly, however, mutations in Ku70, which functions in capping chromosome ends during NHEJ, have only moderate, possibly indirect effects, suggesting that both Lig4 and the alternate end-joining ligase act in some retrohoming events independently of Ku. Another potential Lig4-independent mechanism, reverse transcriptase template switching from the intron RNA to the upstream exon DNA, occurs in vitro, but gives junctions differing from the majority in vivo. Our results show that group II introns can utilize cellular NHEJ enzymes for retromobility in higher organisms, possibly exploiting mechanisms that contribute to retrotransposition and mitigate DNA damage by resident retrotransposons. Additionally, our results reveal novel activities of group II intron reverse transcriptases, with implications for retrohoming mechanisms and potential biotechnological applications

HPV16 E7-Dependent Transformation Activates NHE1 through a PKA-RhoA-Iinduced Inhibition of p38alpha

Background: Neoplastic transformation originates from a large number of different genetic alterations. Despite this genetic variability, a common phenotype to transformed cells is cellular alkalinization. We have previously shown in human keratinocytes and a cell line in which transformation can be turned on and followed by the inducible expression of the E7 oncogene of human papillomavirus type 16 (HPV16), that intracellular alkalinization is an early and essential physiological event driven by the up-regulation of the Na/H-+(+) exchanger isoform 1 (NHE1) and is necessary for the development of other transformed phenotypes and the in vivo tumor formation in nude mice.Methodology: Here, we utilize these model systems to elucidate the dynamic sequence of alterations of the upstream signal transduction systems leading to the transformation-dependent activation of NHE1.Principal Findings: We observe that a down-regulation of p38 MAPK activity is a fundamental step in the ability of the oncogene to transform the cell. Further, using pharmacological agents and transient transfections with dominant interfering, constitutively active, phosphorylation negative mutants and siRNA strategy to modify specific upstream signal transduction components that link HPV16 E7 oncogenic signals to up-regulation of the NHE1, we demonstrate that the stimulation of NHE1 activity is driven by an early rise in cellular cAMP resulting in the down-stream inhibition of p38 MAPK via the PKA-dependent phosphorylation of the small G-protein, RhoA, and its subsequent inhibition.Conclusions: All together these data significantly improve our knowledge concerning the basic cellular alterations involved in oncogene-driven neoplastic transformation

Structural basis of RhoA activation by leukemia-associated RhoGEF

textSmall GTPases of the Rho family are regulators of cytoskeletal organization, neuronal morphogenesis, and transcription and are implicated in the development of cancer. The activity of Rho GTPases is dependent on their binding to GDP (inactive) or GTP (active) and is tightly regulated by accessory proteins. RhoGEFs (Rho guanine nucleotide exchange factors) activate Rho GTPases by stabilizing their nucleotide-free form via a DH/PH (Dbl homology/Pleckstrin homology) domains module. Leukemiaassociated RhoGEF (LARG) belongs to a subfamily of RhoGEFs, the RH-RhoGEFs, that also contain an RH (Regulator of G protein signaling homology) domain N-terminal to the DH/PH domains and specifically activate RhoA, and not the two other RhoGTPases, Cdc42 and Rac1. RH-RhoGEFs are coupled to G protein-coupled receptor (GPCR) activation because their RH domains interact with Gα12/13 proteins, which in turn activates their GEF activity. Although the LARG DH domain is sufficient for catalysis the PH domain contributes to nucleotide exchange. To better understand how the LARG PH and RH domains contribute to its activity, and to elucidate the structural determinants of RhoA-specificity, structures of the LARG DH/PH domains alone and in complex with RhoA were determined by x-ray crystallography at 2.1 and 3.2 Å resolution, respectively. To verify the structural findings, mutants were generated and assessed using a fluorescence assay. A novel N-terminal subdomain of the DH domain was discovered, which seems to be important for activity and might be used as a switch to regulate LARG activity. The sequence of this N-terminal extension is conserved throughout the Lbc family of RhoA specific RhoGEFs. PH domain-assisted nucleotide exchange in LARG is dependent on the structural integrity of the junction between the DH and PH domains and not on direct contacts of the PH domain with RhoA. A hydrophobic patch on the PH domain has been discovered, which might be a protein-docking site that is used for regulation of GEF activity by other proteins (e.g Gα13) or other domains within LARG (e.g. RH domain). Fluorescence assays of fragments including the RH domain showed an inhibitory effect of the RH domain in vitro, the structural basis of which will be investigated in future experiments.Biochemistr

Bacterial expression and purification of active hematopoietic cell kinase

Src family kinases (SFKs) are traditionally purified from eukaryotic expression systems. These expression systems can be costly, yield heterogeneously phosphorylated protein samples and present difficulties when metabolic labeling is required for structural studies. Therefore, many attempts have been made to develop bacterial purification systems for SFKs. So far, high-yield bacterial expression systems have only been achieved for SFK kinase domains or for inactive mutants of constructs containing the regulatory SH3 and SH2 domains, but not for their active forms. Herein described is a bacterial expression system for the wild type, active SFK Hck containing SH3, SH2 and kinase domains. Hck plays an important role in phagocyte function as well as the etiology of chronic myeloid leukemia as Hck is an interaction partner of Bcr-Abl. Structural studies of Hck are essential to fully understand the signaling processes involved in host defense and leukemogenesis. Successful bacterial expression of Hck was possible by a dual strategy: (1) co-expression with YopH phosphatase in order to control host toxicity, and (2) expression in a bacterial strain that is RNase E deficient, which dramatically increased overall expression levels. The expressed Hck construct is unphosphorylated and appears to be in an open conformation. Bacterially expressed Hck is capable of autophosphorylation, phosphorylates substrate at rates comparable to insect cell expressed Hck, and can be inhibited by staurosporine and Csk.Romana Kristelly, Theresa Wenli Qiu, Natalie J. Gunn, Denis B. Scanlon, Terrence D. Mulher