Construction and validation of a regulatory network for pluripotency and self-renewal of mouse embryonic stem cells.

A 30-node signed and directed network responsible for self-renewal and pluripotency of mouse embryonic stem cells (mESCs) was extracted from several ChIP-Seq and knockdown followed by expression prior studies. The underlying regulatory logic among network components was then learned using the initial network topology and single cell gene expression measurements from mESCs cultured in serum/LIF or serum-free 2i/LIF conditions. Comparing the learned network regulatory logic derived from cells cultured in serum/LIF vs. 2i/LIF revealed differential roles for Nanog, Oct4/Pou5f1, Sox2, Esrrb and Tcf3. Overall, gene expression in the serum/LIF condition was more variable than in the 2i/LIF but mostly consistent across the two conditions. Expression levels for most genes in single cells were bimodal across the entire population and this motivated a Boolean modeling approach. In silico predictions derived from removal of nodes from the Boolean dynamical model were validated with experimental single and combinatorial RNA interference (RNAi) knockdowns of selected network components. Quantitative post-RNAi expression level measurements of remaining network components showed good agreement with the in silico predictions. Computational removal of nodes from the Boolean network model was also used to predict lineage specification outcomes. In summary, data integration, modeling, and targeted experiments were used to improve our understanding of the regulatory topology that controls mESC fate decisions as well as to develop robust directed lineage specification protocols.

Comparison of somatic mutation calling methods in amplicon and whole exome sequence data.

BACKGROUND: High-throughput sequencing is rapidly becoming common practice in clinical diagnosis and cancer research. Many algorithms have been developed for somatic single nucleotide variant (SNV) detection in matched tumor-normal DNA sequencing. Although numerous studies have compared the performance of various algorithms on exome data, there has not yet been a systematic evaluation using PCR-enriched amplicon data with a range of variant allele fractions. The recently developed gold standard variant set for the reference individual NA12878 by the NIST-led "Genome in a Bottle" Consortium (NIST-GIAB) provides a good resource to evaluate admixtures with various SNV fractions. RESULTS: Using the NIST-GIAB gold standard, we compared the performance of five popular somatic SNV calling algorithms (GATK UnifiedGenotyper followed by simple subtraction, MuTect, Strelka, SomaticSniper and VarScan2) for matched tumor-normal amplicon and exome sequencing data. CONCLUSIONS: We demonstrated that the five commonly used somatic SNV calling methods are applicable to both targeted amplicon and exome sequencing data. However, the sensitivities of these methods vary based on the allelic fraction of the mutation in the tumor sample. Our analysis can assist researchers in choosing a somatic SNV calling method suitable for their specific needs.

Effect of Clonal Hematopoiesis Mutations and Canakinumab Treatment on Incidence of Solid Tumors in the CANTOS Randomized Clinical Trial.

Clonal hematopoiesis (CH) is more common in older persons and has been associated with an increased risk of hematological cancers and cardiovascular diseases. The most common CH mutations occur in the DNMT3A and TET2 genes and result in increased pro-inflammatory signaling. The Canakinumab Anti-inflammatory Thrombosis Outcome Study (CANTOS, NCT01327846) evaluated the neutralizing anti-IL-1ß antibody canakinumab in 10,061 randomized patients with a history of myocardial infarction and persistent inflammation; DNA samples were available from 3,923 patients for targeted genomic sequencing. We examined the incidence of non-hematological malignancy by treatment assignment and CH mutations and estimated the cumulative incidence of malignancy events during trial follow-up. Patients with TET2 mutations treated with canakinumab had the lowest incidence of non-hematological malignancy across cancer types. The cumulative incidence of at least one reported malignancy was lower for patients with TET2 mutations treated with canakinumab vs those treated with placebo. These findings support a potential role for canakinumab in cancer prevention and provide evidence of IL-1ß blockade cooperating with CH mutations to modify the disease course.

Aptamer Proteomics for Biomarker Discovery in Heart Failure With Preserved Ejection Fraction: The PARAGON-HF Proteomic Substudy.

BACKGROUND: Prognostic markers and biological pathways linked to detrimental clinical outcomes in heart failure with preserved ejection fraction (HFpEF) remain incompletely defined. METHODS AND RESULTS: We measured serum levels of 4123 unique proteins in 1117 patients with HFpEF enrolled in the PARAGON-HF (Efficacy and Safety of LCZ696 Compared to Valsartan, on Morbidity and Mortality in Heart Failure Patients With Preserved Ejection Fraction) trial using a modified aptamer proteomic assay. Baseline circulating protein concentrations significantly associated with the primary end point and the timing and occurrence of total heart failure hospitalization and cardiovascular death were identified by recurrent events regression, accounting for multiple testing, adjusted for age, sex, treatment, and anticoagulant use, and compared with published analyses in 2515 patients with heart failure with reduced ejection fraction from the PARADIGM-HF (Prospective Comparison of ARNI With ACEI to Determine Impact on Global Mortality and Morbidity in Heart Failure) and ATMOSPHERE (Efficacy and Safety of Aliskiren and Aliskiren/Enalapril Combination on Morbidity-Mortality in Patients With Chronic Heart Failure) clinical trials. We identified 288 proteins that were robustly associated with the risk of heart failure hospitalization and cardiovascular death in patients with HFpEF. The baseline proteins most strongly related to outcomes included B2M (ß-2 microglobulin), TIMP1 (tissue inhibitor of matrix metalloproteinase 1), SERPINA4 (serpin family A member 4), and SVEP1 (sushi, von Willebrand factor type A, EGF, and pentraxin domain containing 1). Overall, the protein-outcome associations in patients with HFpEF did not markedly differ as compared with patients with heart failure with reduced ejection fraction. A proteomic risk score derived in patients with HFpEF was not superior to a previous proteomic score derived in heart failure with reduced ejection fraction nor to clinical risk factors, NT-proBNP (N-terminal pro-B-type natriuretic peptide), or high-sensitivity cardiac troponin. CONCLUSIONS: Numerous serum proteins linked to metabolic, coagulation, and extracellular matrix regulatory pathways were associated with worse HFpEF prognosis in the PARAGON-HF proteomic substudy. Our results demonstrate substantial similarities among serum proteomic risk markers for heart failure hospitalization and cardiovascular death when comparing clinical trial participants with heart failure across the ejection fraction spectrum. REGISTRATION: URL: https://www.clinicaltrials.gov; Unique Identifiers: NCT01920711, NCT01035255, NCT00853658.

Expression2Kinases: mRNA profiling linked to multiple upstream regulatory layers.

MOTIVATION: Genome-wide mRNA profiling provides a snapshot of the global state of cells under different conditions. However, mRNA levels do not provide direct understanding of upstream regulatory mechanisms. Here, we present a new approach called Expression2Kinases (X2K) to identify upstream regulators likely responsible for observed patterns in genome-wide gene expression. By integrating chromatin immuno-precipitation (ChIP)-seq/chip and position weight matrices (PWMs) data, protein-protein interactions and kinase-substrate phosphorylation reactions, we can better identify regulatory mechanisms upstream of genome-wide differences in gene expression. We validated X2K by applying it to recover drug targets of food and drug administration (FDA)-approved drugs from drug perturbations followed by mRNA expression profiling; to map the regulatory landscape of 44 stem cells and their differentiating progeny; to profile upstream regulatory mechanisms of 327 breast cancer tumors; and to detect pathways from profiled hepatic stellate cells and hippocampal neurons. The X2K approach can advance our understanding of cell signaling and unravel drugs mechanisms of action. AVAILABILITY: The software and source code are freely available at: http://www.maayanlab.net/X2K. CONTACT: avi.maayan@mssm.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Effects of IL-1ß inhibition on anemia and clonal hematopoiesis in the randomized CANTOS trial.

Canakinumab, a monoclonal antibody targeting proinflammatory cytokine interleukin-1ß (IL-1ß), improved hemoglobin levels while preventing recurrent cardiovascular events in the Canakinumab Anti-inflammatory Thrombosis Outcomes Study (CANTOS). This cardiovascular (CV) preventive effect was greater in patients with TET2 mutations associated with clonal hematopoiesis (CH). The current proteogenomic analysis aimed to understand the clinical response to canakinumab and underlying proteomic profiles in the context of CH and anemia. The analysis included 4595 patients from the CANTOS study who received either canakinumab or placebo and evaluated multiplexed proteomics (4785 proteins) using SomaScan and targeted deep sequencing for CH mutations. Incident anemia was more common in the presence of CH mutations but reduced by canakinumab treatment. Canakinumab treatment was significantly associated with higher hemoglobin increment in patients with concurrent CH mutations and anemia than patients with CH mutations without anemia or without CH mutations. Compared with those without CH mutations, the presence of CH mutations was associated with proteomic signatures of inflammation and defense response to infection, as well as markers of high-risk CV disease which was further enhanced by the presence of anemia. Canakinumab suppressed hepcidin, proinflammatory cytokines, myeloid activation, and complement pathways, and reversed pathologically deregulated pathways to a greater extent in patients with CH mutations and anemia. These molecular findings provide evidence of the clinical use of IL-1ß blockade and support further study of canakinumab for patients with concurrent anemia and CH mutations. This study was registered at www.clinicaltrials.gov as #NCT01327846.

Clinical characteristics of heart failure with reduced ejection fraction patients with rare pathogenic variants in dilated cardiomyopathy-associated genes: A subgroup analysis of the PARADIGM-HF trial.

AIMS: To evaluate the prevalence of pathogenic variants in genes associated with dilated cardiomyopathy (DCM) in a clinical trial population with heart failure and reduced ejection fraction (HFrEF) and describe the baseline characteristics by variant carrier status. METHODS AND RESULTS: This was a post hoc analysis of the Phase 3 PARADIGM-HF trial. Forty-four genes, divided into three tiers, based on definitive, moderate or limited evidence of association with DCM, were assessed for rare predicted loss-of-function (pLoF) variants, which were prioritized using ClinVar annotations, measures of gene transcriptional output and evolutionary constraint, and pLoF confidence predictions. Prevalence was reported for pLoF variant carriers based on DCM-associated gene tiers. Clinical features were compared between carriers and non-carriers. Of the 1412 HFrEF participants with whole-exome sequence data, 68 (4.8%) had at least one pLoF variant in the 8 tier-1 genes (definitive/strong association with DCM), with Titin being most commonly affected. The prevalence increased to 7.5% when considering all 44 genes. Among patients with idiopathic aetiology, 10.0% (23/229) had tier-1 variants only and 12.6% (29/229) had tier-1, -2 or -3 variants. Compared to non-carriers, tier-1 carriers were younger (4 years; adjusted p-value [padj ] = 4 × 10-3 ), leaner (27.8 kg/m2 vs. 29.4 kg/m2 ; padj = 3.2 × 10-3 ), had lower ejection fraction (27.3% vs. 29.8%; padj = 5.8 × 10-3 ), and less likely to have ischaemic aetiology (37.3% vs. 67.4%; padj = 4 × 10-4 ). CONCLUSION: Deleterious pLoF variants in genes with definitive/strong association with DCM were identified in ∼5% of HFrEF patients from a PARADIGM-HF trial subset, who were younger, had lower ejection fraction and were less likely to have had an ischaemic aetiology.

Network- and attribute-based classifiers can prioritize genes and pathways for autism spectrum disorders and intellectual disability.

Autism spectrum disorders (ASD) are a group of related neurodevelopmental disorders with significant combined prevalence (∼1%) and high heritability. Dozens of individually rare genes and loci associated with high-risk for ASD have been identified, which overlap extensively with genes for intellectual disability (ID). However, studies indicate that there may be hundreds of genes that remain to be identified. The advent of inexpensive massively parallel nucleotide sequencing can reveal the genetic underpinnings of heritable complex diseases, including ASD and ID. However, whole exome sequencing (WES) and whole genome sequencing (WGS) provides an embarrassment of riches, where many candidate variants emerge. It has been argued that genetic variation for ASD and ID will cluster in genes involved in distinct pathways and protein complexes. For this reason, computational methods that prioritize candidate genes based on additional functional information such as protein-protein interactions or association with specific canonical or empirical pathways, or other attributes, can be useful. In this study we applied several supervised learning approaches to prioritize ASD or ID disease gene candidates based on curated lists of known ASD and ID disease genes. We implemented two network-based classifiers and one attribute-based classifier to show that we can rank and classify known, and predict new, genes for these neurodevelopmental disorders. We also show that ID and ASD share common pathways that perturb an overlapping synaptic regulatory subnetwork. We also show that features relating to neuronal phenotypes in mouse knockouts can help in classifying neurodevelopmental genes. Our methods can be applied broadly to other diseases helping in prioritizing newly identified genetic variation that emerge from disease gene discovery based on WES and WGS.

Activation of alternate prosurvival pathways accounts for acquired sunitinib resistance in U87MG glioma xenografts.

Acquired drug resistance represents a major obstacle to using sunitinib for the treatment of solid tumors. Here, we examined the cellular and molecular alterations in tumors that are associated with acquired brain tumor resistance to sunitinib by using an in vivo model. U87MG tumors obtained from nude mice that received sunitinib (40 mg/kg/day) for 30 days were classified into sunitinib-sensitive and -resistant groups based on tumor volume and underwent targeted gene microarray and protein array analyses. The expression of several angiogenesis-associated genes was significantly modulated in sunitinib-treated tumors compared with those in control tumors (p<0.05), whereas no significant differences were observed between sunitinib-sensitive and -resistant tumors (p>0.05). Tumor vasculature based on microvessel density, neurogenin 2 chondroitin sulfate proteoglycan density, and α-smooth muscle actin density was also similar in sunitinib-treatment groups (p>0.05). The moderate increase in unbound sunitinib tumor-to-plasma area-under-the-curve ratio in sunitinib-resistant mice was accompanied by up-regulated ATP-binding cassette G2 expression in tumor. The most profound difference between the sunitinib-sensitive and -resistant groups was found in the expression of several phosphorylated proteins involved in intracellular signaling. In particular, phospholipase C-γ1 phosphorylation in sunitinib-resistant tumors was up-regulated by 2.6-fold compared with that in sunitinib-sensitive tumors (p<0.05). In conclusion, acquired sunitinib resistance in U87MG tumors is not associated with revascularization in tumors, but rather with the activation of alternate prosurvival pathways involved in an escape mechanism facilitating tumor growth and possibly insufficient drug uptake in tumor cells caused by an up-regulated membrane efflux transporter.

Esrrb Regulates Specific Feed-Forward Loops to Transit From Pluripotency Into Early Stages of Differentiation.

Characterization of pluripotent states, in which cells can both self-renew or differentiate, with the irreversible loss of pluripotency, are important research areas in developmental biology. Although microRNAs (miRNAs) have been shown to play a relevant role in cellular differentiation, the role of miRNAs integrated into gene regulatory networks and its dynamic changes during these early stages of embryonic stem cell (ESC) differentiation remain elusive. Here we describe the dynamic transcriptional regulatory circuitry of stem cells that incorporate protein-coding and miRNA genes based on miRNA array expression and quantitative sequencing of short transcripts upon the downregulation of the Estrogen Related Receptor Beta (Esrrb). The data reveals how Esrrb, a key stem cell transcription factor, regulates a specific stem cell miRNA expression program and integrates dynamic changes of feed-forward loops contributing to the early stages of cell differentiation upon its downregulation. Together these findings provide new insights on the architecture of the combined transcriptional post-transcriptional regulatory network in embryonic stem cells.

TET2-Driven Clonal Hematopoiesis and Response to Canakinumab: An Exploratory Analysis of the CANTOS Randomized Clinical Trial.

Importance: Clonal hematopoiesis of indeterminate potential (CHIP) is associated with increased risk of atherosclerotic cardiovascular disease, and mouse experiments suggest that CHIP related to Tet2 loss of function in myeloid cells accelerates atherosclerosis via augmented interleukin (IL) 1ß signaling. Objective: To assess whether individuals with CHIP have greater cardiovascular event reduction in response to IL-1ß neutralization in the Canankinumab Anti-inflammatory Thrombosis Outcomes Trial (CANTOS). Design, Setting, and Participants: This randomized clinical trial took place from April 2011 to June 2017 at more than 1000 clinical sites in 39 countries. Targeted deep sequencing of genes previously associated with CHIP in a subset of trial participants using genomic DNA prepared from baseline peripheral blood samples were analyzed. All participants had prior myocardial infarction and elevated high-sensitivity C-reactive protein level above 0.20 mg/dL. Analysis took place between June 2017 and December 2021. Interventions: Canakinumab, an anti-IL-1ß antibody, given at doses of 50, 150, and 300 mg once every 3 months. Main Outcomes and Measures: Major adverse cardiovascular events (MACE). Results: A total of 338 patients (8.6%) were identified in this subset with evidence for clonal hematopoiesis. As expected, the incidence of CHIP increased with age; the mean (SD) age of patients with CHIP was 66.3 (9.2) years and 61.5 (9.6) years in patients without CHIP. Unlike other populations that were not preselected for elevated C-reactive protein, in the CANTOS population variants in TET2 were more common than DNMT3A (119 variants in 103 patients vs 86 variants in 85 patients). Placebo-treated patients with CHIP showed a nonsignificant increase in the rate of MACE compared with patients without CHIP using a Cox proportional hazard model (hazard ratio, 1.32 [95% CI, 0.86-2.04]; P = .21). Exploratory analyses of placebo-treated patients with a somatic variant in either TET2 or DNMT3A (n = 58) showed an equivocal risk for MACE (hazard ratio, 1.65 [95% CI, 0.97-2.80]; P = .06). Patients with CHIP due to somatic variants in TET2 also had reduced risk for MACE while taking canakinumab (hazard ratio, 0.38 [95% CI, 0.15-0.96]) with equivocal difference compared with others (P for interaction = .14). Conclusions and Relevance: These results are consistent with observations of increased risk for cardiovascular events in patients with CHIP and raise the possibility that those with TET2 variants may respond better to canakinumab than those without CHIP. Future studies are required to further substantiate this hypothesis. Trial Registration: ClinicalTrials.gov Identifier: NCT01327846.

ChEA: transcription factor regulation inferred from integrating genome-wide ChIP-X experiments.

MOTIVATION: Experiments such as ChIP-chip, ChIP-seq, ChIP-PET and DamID (the four methods referred herein as ChIP-X) are used to profile the binding of transcription factors to DNA at a genome-wide scale. Such experiments provide hundreds to thousands of potential binding sites for a given transcription factor in proximity to gene coding regions. RESULTS: In order to integrate data from such studies and utilize it for further biological discovery, we collected interactions from such experiments to construct a mammalian ChIP-X database. The database contains 189,933 interactions, manually extracted from 87 publications, describing the binding of 92 transcription factors to 31,932 target genes. We used the database to analyze mRNA expression data where we perform gene-list enrichment analysis using the ChIP-X database as the prior biological knowledge gene-list library. The system is delivered as a web-based interactive application called ChIP Enrichment Analysis (ChEA). With ChEA, users can input lists of mammalian gene symbols for which the program computes over-representation of transcription factor targets from the ChIP-X database. The ChEA database allowed us to reconstruct an initial network of transcription factors connected based on shared overlapping targets and binding site proximity. To demonstrate the utility of ChEA we present three case studies. We show how by combining the Connectivity Map (CMAP) with ChEA, we can rank pairs of compounds to be used to target specific transcription factor activity in cancer cells. AVAILABILITY: The ChEA software and ChIP-X database is freely available online at: http://amp.pharm.mssm.edu/lib/chea.jsp.

An Esrrb and Nanog Cell Fate Regulatory Module Controlled by Feed Forward Loop Interactions.

Cell fate decisions during development are governed by multi-factorial regulatory mechanisms including chromatin remodeling, DNA methylation, binding of transcription factors to specific loci, RNA transcription and protein synthesis. However, the mechanisms by which such regulatory "dimensions" coordinate cell fate decisions are currently poorly understood. Here we quantified the multi-dimensional molecular changes that occur in mouse embryonic stem cells (mESCs) upon depletion of Estrogen related receptor beta (Esrrb), a key pluripotency regulator. Comparative analyses of expression changes subsequent to depletion of Esrrb or Nanog, indicated that a system of interlocked feed-forward loops involving both factors, plays a central part in regulating the timing of mESC fate decisions. Taken together, our meta-analyses support a hierarchical model in which pluripotency is maintained by an Oct4-Sox2 regulatory module, while the timing of differentiation is regulated by a Nanog-Esrrb module.

Identification and developmental expression of Dec2 in zebrafish.

The involvement of Dec2, a member of the basic helix-loop-helix (bHLH) family, in cellular differentiation, hypoxia response, and circadian regulation has been investigated. Here we report the previously unknown spatiotemporal expression of Dec2 in zebrafish embryogenesis. Dec2 is dynamically expressed in zebrafish pineal gland, tract of the postoptic commissure, brain, notochord, heart, common cardinal vein (CCV), axial vein, pronephric duct, swim bladder, and early somites during embryogenesis, which implies that Dec2 is involved in zebrafish central nervous system development, cardiogenesis, and internal organs and somites formation. The embryonic expression patterns of zebrafish Dec2 and its homolog Dec1 partially overlap, but are distinct from each other. The Dec2 expression level was lower than that of Dec1 during zebrafish embryogenesis. Although Dec1 also contributed to zebrafish somites formation, cardiogenesis, and internal organs and central nervous system development, the two Dec genes were not likely to be simply redundant during zebrafish embryogenesis. Our results imply that Dec2, like its homolog Dec1, is involved in zebrafish cardiogenesis, central nervous system development, and internal organs and somites formation with distinct developmental roles.

Selection of GalNAc-conjugated siRNAs with limited off-target-driven rat hepatotoxicity.

Small interfering RNAs (siRNAs) conjugated to a trivalent N-acetylgalactosamine (GalNAc) ligand are being evaluated in investigational clinical studies for a variety of indications. The typical development candidate selection process includes evaluation of the most active compounds for toxicity in rats at pharmacologically exaggerated doses. The subset of GalNAc-siRNAs that show rat hepatotoxicity is not advanced to clinical development. Potential mechanisms of hepatotoxicity can be associated with the intracellular accumulation of oligonucleotides and their metabolites, RNA interference (RNAi)-mediated hybridization-based off-target effects, and/or perturbation of endogenous RNAi pathways. Here we show that rodent hepatotoxicity observed at supratherapeutic exposures can be largely attributed to RNAi-mediated off-target effects, but not chemical modifications or the perturbation of RNAi pathways. Furthermore, these off-target effects can be mitigated by modulating seed-pairing using a thermally destabilizing chemical modification, which significantly improves the safety profile of a GalNAc-siRNA in rat and may minimize the occurrence of hepatotoxic siRNAs across species.

Reversal of siRNA-mediated gene silencing in vivo.

We report rapid, potent reversal of GalNAc-siRNA-mediated RNA interference (RNAi) activity in vivo with short, synthetic, high-affinity oligonucleotides complementary to the siRNA guide strand. We found that 9-mers with five locked nucleic acids (LNAs) have the highest potency across several targets. Our modular, sequence-specific approach, named REVERSIR, may enhance the therapeutic profile of any long-acting GalNAc-siRNA (short interfering RNA) conjugate by enabling control of RNAi pharmacology.

A Randomized Controlled Trial Comparing Suture-Fixation Mucopexy and Doppler-Guided Hemorrhoidal Artery Ligation in Patients with Grade III Hemorrhoids.

Background. We aimed to evaluate the effectiveness of a suture-fixation mucopexy procedure by comparing with Doppler-guided hemorrhoidal artery ligation (DGHAL) in the management of patients with grade III hemorrhoids. Methods. This was a randomized controlled trial. One hundred patients with grade III hemorrhoids were randomly assigned to receive suture-fixation mucopexy (n = 50) or DGHAL (n = 50). Outcome assessments were performed at 2 weeks, 12 months, and 24 months. Assessments included resolution of clinical symptoms, postoperative complications, duration of hospitalization, and total costs. Results. At 2 weeks, one (2%) patient in suture-fixation group and four (8%) patients in DGHAL group had persistent prolapsing hemorrhoids. Postoperative bleeding was observed in two patients (4%) in suture-fixation group and one patient in DGHAL group. There was no significant difference in short-term recurrence between groups. Postoperative complications and duration of hospitalization were comparable between the two groups. Rates of recurrence of prolapse or bleeding at 12 months did not differ between groups. However, recurrence of prolapse at 24 months was significantly more common in DGHAL group (19.0% versus 2.3%, p = 0.030). Conclusions. Compared with DGHAL, the suture-fixation mucopexy technique had comparable short-term outcomes and favorable long-term outcomes.

ESCAPE: database for integrating high-content published data collected from human and mouse embryonic stem cells.

High content studies that profile mouse and human embryonic stem cells (m/hESCs) using various genome-wide technologies such as transcriptomics and proteomics are constantly being published. However, efforts to integrate such data to obtain a global view of the molecular circuitry in m/hESCs are lagging behind. Here, we present an m/hESC-centered database called Embryonic Stem Cell Atlas from Pluripotency Evidence integrating data from many recent diverse high-throughput studies including chromatin immunoprecipitation followed by deep sequencing, genome-wide inhibitory RNA screens, gene expression microarrays or RNA-seq after knockdown (KD) or overexpression of critical factors, immunoprecipitation followed by mass spectrometry proteomics and phosphoproteomics. The database provides web-based interactive search and visualization tools that can be used to build subnetworks and to identify known and novel regulatory interactions across various regulatory layers. The web-interface also includes tools to predict the effects of combinatorial KDs by additive effects controlled by sliders, or through simulation software implemented in MATLAB. Overall, the Embryonic Stem Cell Atlas from Pluripotency Evidence database is a comprehensive resource for the stem cell systems biology community. Database URL: http://www.maayanlab.net/ESCAPE

Oct4 links multiple epigenetic pathways to the pluripotency network.

Oct4 is a well-known transcription factor that plays fundamental roles in stem cell self-renewal, pluripotency, and somatic cell reprogramming. However, limited information is available on Oct4-associated protein complexes and their intrinsic protein-protein interactions that dictate Oct4's critical regulatory activities. Here we employed an improved affinity purification approach combined with mass spectrometry to purify Oct4 protein complexes in mouse embryonic stem cells (mESCs), and discovered many novel Oct4 partners important for self-renewal and pluripotency of mESCs. Notably, we found that Oct4 is associated with multiple chromatin-modifying complexes with documented as well as newly proved functional significance in stem cell maintenance and somatic cell reprogramming. Our study establishes a solid biochemical basis for genetic and epigenetic regulation of stem cell pluripotency and provides a framework for exploring alternative factor-based reprogramming strategies.