GPR146 Deficiency Protects against Hypercholesterolemia and Atherosclerosis.

Although human genetic studies have implicated many susceptible genes associated with plasma lipid levels, their physiological and molecular functions are not fully characterized. Here we demonstrate that orphan G protein-coupled receptor 146 (GPR146) promotes activity of hepatic sterol regulatory element binding protein 2 (SREBP2) through activation of the extracellular signal-regulated kinase (ERK) signaling pathway, thereby regulating hepatic very low-density lipoprotein (VLDL) secretion, and subsequently circulating low-density lipoprotein cholesterol (LDL-C) and triglycerides (TG) levels. Remarkably, GPR146 deficiency reduces plasma cholesterol levels substantially in both wild-type and LDL receptor (LDLR)-deficient mice. Finally, aortic atherosclerotic lesions are reduced by 90% and 70%, respectively, in male and female LDLR-deficient mice upon GPR146 depletion. Taken together, these findings outline a regulatory role for the GPR146/ERK axis in systemic cholesterol metabolism and suggest that GPR146 inhibition could be an effective strategy to reduce plasma cholesterol levels and atherosclerosis.

Genome-wide chromatin state transitions associated with developmental and environmental cues.

Differences in chromatin organization are key to the multiplicity of cell states that arise from a single genetic background, yet the landscapes of in vivo tissues remain largely uncharted. Here, we mapped chromatin genome-wide in a large and diverse collection of human tissues and stem cells. The maps yield unprecedented annotations of functional genomic elements and their regulation across developmental stages, lineages, and cellular environments. They also reveal global features of the epigenome, related to nuclear architecture, that also vary across cellular phenotypes. Specifically, developmental specification is accompanied by progressive chromatin restriction as the default state transitions from dynamic remodeling to generalized compaction. Exposure to serum in vitro triggers a distinct transition that involves de novo establishment of domains with features of constitutive heterochromatin. We describe how these global chromatin state transitions relate to chromosome and nuclear architecture, and discuss their implications for lineage fidelity, cellular senescence, and reprogramming.

Disease-specific induced pluripotent stem cells.

Tissue culture of immortal cell strains from diseased patients is an invaluable resource for medical research but is largely limited to tumor cell lines or transformed derivatives of native tissues. Here we describe the generation of induced pluripotent stem (iPS) cells from patients with a variety of genetic diseases with either Mendelian or complex inheritance; these diseases include adenosine deaminase deficiency-related severe combined immunodeficiency (ADA-SCID), Shwachman-Bodian-Diamond syndrome (SBDS), Gaucher disease (GD) type III, Duchenne (DMD) and Becker muscular dystrophy (BMD), Parkinson disease (PD), Huntington disease (HD), juvenile-onset, type 1 diabetes mellitus (JDM), Down syndrome (DS)/trisomy 21, and the carrier state of Lesch-Nyhan syndrome. Such disease-specific stem cells offer an unprecedented opportunity to recapitulate both normal and pathologic human tissue formation in vitro, thereby enabling disease investigation and drug development.

Inducers of the endothelial cell barrier identified through chemogenomic screening in genome-edited hPSC-endothelial cells.

The blood-retina barrier and blood-brain barrier (BRB/BBB) are selective and semipermeable and are critical for supporting and protecting central nervous system (CNS)-resident cells. Endothelial cells (ECs) within the BRB/BBB are tightly coupled, express high levels of Claudin-5 (CLDN5), a junctional protein that stabilizes ECs, and are important for proper neuronal function. To identify novel CLDN5 regulators (and ultimately EC stabilizers), we generated a CLDN5-P2A-GFP stable cell line from human pluripotent stem cells (hPSCs), directed their differentiation to ECs (CLDN5-GFP hPSC-ECs), and performed flow cytometry-based chemogenomic library screening to measure GFP expression as a surrogate reporter of barrier integrity. Using this approach, we identified 62 unique compounds that activated CLDN5-GFP. Among them were TGF-ß pathway inhibitors, including RepSox. When applied to hPSC-ECs, primary brain ECs, and retinal ECs, RepSox strongly elevated barrier resistance (transendothelial electrical resistance), reduced paracellular permeability (fluorescein isothiocyanate-dextran), and prevented vascular endothelial growth factor A (VEGFA)-induced barrier breakdown in vitro. RepSox also altered vascular patterning in the mouse retina during development when delivered exogenously. To determine the mechanism of action of RepSox, we performed kinome-, transcriptome-, and proteome-profiling and discovered that RepSox inhibited TGF-ß, VEGFA, and inflammatory gene networks. In addition, RepSox not only activated vascular-stabilizing and barrier-establishing Notch and Wnt pathways, but also induced expression of important tight junctions and transporters. Taken together, our data suggest that inhibiting multiple pathways by selected individual small molecules, such as RepSox, may be an effective strategy for the development of better BRB/BBB models and novel EC barrier-inducing therapeutics.

Functional Screening of Candidate Causal Genes for Insulin Resistance in Human Preadipocytes and Adipocytes.

Rationale: Genome-wide association studies have identified genetic loci associated with insulin resistance (IR) but pinpointing the causal genes of a risk locus has been challenging. Objective: To identify candidate causal genes for IR, we screened regional and biologically plausible genes (16 in total) near the top 10 IR-loci in risk-relevant cell types, namely preadipocytes and adipocytes. Methods and Results: We generated 16 human Simpson-Golabi-Behmel syndrome preadipocyte knockout lines each with a single IR-gene knocked out by lentivirus-mediated CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 system. We evaluated each gene knockout by screening IR-relevant phenotypes in the 3 insulin-sensitizing mechanisms, including adipogenesis, lipid metabolism, and insulin signaling. We performed genetic analyses using data on the genotype-tissue expression portal expression quantitative trait loci database and accelerating medicines partnership type 2 diabetes mellitus Knowledge Portal to evaluate whether candidate genes prioritized by our in vitro studies were expression quantitative trait loci genes in human subcutaneous adipose tissue, and whether expression of these genes is associated with risk of IR, type 2 diabetes mellitus, and cardiovascular diseases. We further validated the functions of 3 new adipose IR genes by overexpression-based phenotypic rescue in the Simpson-Golabi-Behmel syndrome preadipocyte knockout lines. Twelve genes, PPARG, IRS-1, FST, PEPD, PDGFC, MAP3K1, GRB14, ARL15, ANKRD55, RSPO3, COBLL1, and LYPLAL1, showed diverse phenotypes in the 3 insulin-sensitizing mechanisms, and the first 7 of these genes could affect all the 3 mechanisms. Five out of 6 expression quantitative trait loci genes are among the top candidate causal genes and the abnormal expression levels of these genes (IRS-1, GRB14, FST, PEPD, and PDGFC) in human subcutaneous adipose tissue could be associated with increased risk of IR, type 2 diabetes mellitus, and cardiovascular disease. Phenotypic rescue by overexpression of the candidate causal genes (FST, PEPD, and PDGFC) in the Simpson-Golabi-Behmel syndrome preadipocyte knockout lines confirmed their function in adipose IR. Conclusions: Twelve genes showed diverse phenotypes indicating differential roles in insulin sensitization, suggesting mechanisms bridging the association of their genomic loci with IR. We prioritized PPARG, IRS-1, GRB14, MAP3K1, FST, PEPD, and PDGFC as top candidate genes. Our work points to novel roles for FST, PEPD, and PDGFC in adipose tissue, with consequences for cardiometabolic diseases.

Generation of hypoimmunogenic human pluripotent stem cells.

Polymorphic HLAs form the primary immune barrier to cell therapy. In addition, innate immune surveillance impacts cell engraftment, yet a strategy to control both, adaptive and innate immunity, is lacking. Here we employed multiplex genome editing to specifically ablate the expression of the highly polymorphic HLA-A/-B/-C and HLA class II in human pluripotent stem cells. Furthermore, to prevent innate immune rejection and further suppress adaptive immune responses, we expressed the immunomodulatory factors PD-L1, HLA-G, and the macrophage "don't-eat me" signal CD47 from the AAVS1 safe harbor locus. Utilizing in vitro and in vivo immunoassays, we found that T cell responses were blunted. Moreover, NK cell killing and macrophage engulfment of our engineered cells were minimal. Our results describe an approach that effectively targets adaptive as well as innate immune responses and may therefore enable cell therapy on a broader scale.

Interactomics Analyses of Wild-Type and Mutant A1CF Reveal Diverged Functions in Regulating Cellular Lipid Metabolism.

Population genetic studies highlight a missense variant (G398S) of A1CF that is strongly associated with higher levels of blood triglycerides (TGs) and total cholesterol (TC). Functional analyses suggest that the mutation accelerates the secretion of very low-density lipoprotein (VLDL) from the liver by an unknown mechanism. Here, we used multiomics approaches to interrogate the functional difference between the WT and mutant A1CF. Using metabolomics analyses, we captured the cellular lipid metabolite changes induced by transient expression of the proteins, confirming that the mutant A1CF is able to relieve the TG accumulation induced by WT A1CF. Using a proteomics approach, we obtained the interactomic data of WT and mutant A1CF. Networking analyses show that WT A1CF interacts with three functional protein groups, RNA/mRNA processing, cytosolic translation, and, surprisingly, mitochondrial translation. The mutation diminishes these interactions, especially with the group of mitochondrial translation. Differential analyses show that the WT A1CF-interacting proteins most significantly different from the mutant are those for mitochondrial translation, whereas the most significant interacting proteins with the mutant are those for cytoskeleton and vesicle-mediated transport. RNA-seq analyses validate that the mutant, but not the WT, A1CF increases the expression of the genes responsible for cellular transport processes. On the contrary, WT A1CF affected the expression of mitochondrial matrix proteins and increased cell oxygen consumption. Thus, our studies confirm the previous hypothesis that A1CF plays broader roles in regulating gene expression. The interactions of the mutant A1CF with the vesicle-mediated transport machinery provide mechanistic insight in understanding the increased VLDL secretion in the A1CF mutation carriers.

Asialoglycoprotein receptor 1 is a specific cell-surface marker for isolating hepatocytes derived from human pluripotent stem cells.

Hepatocyte-like cells (HLCs) are derived from human pluripotent stem cells (hPSCs) in vitro, but differentiation protocols commonly give rise to a heterogeneous mixture of cells. This variability confounds the evaluation of in vitro functional assays performed using HLCs. Increased differentiation efficiency and more accurate approximation of the in vivo hepatocyte gene expression profile would improve the utility of hPSCs. Towards this goal, we demonstrate the purification of a subpopulation of functional HLCs using the hepatocyte surface marker asialoglycoprotein receptor 1 (ASGR1). We analyzed the expression profile of ASGR1-positive cells by microarray, and tested their ability to perform mature hepatocyte functions (albumin and urea secretion, cytochrome activity). By these measures, ASGR1-positive HLCs are enriched for the gene expression profile and functional characteristics of primary hepatocytes compared with unsorted HLCs. We have demonstrated that ASGR1-positive sorting isolates a functional subpopulation of HLCs from among the heterogeneous cellular population produced by directed differentiation.

A distant trophoblast-specific enhancer controls HLA-G expression at the maternal-fetal interface.

HLA-G, a nonclassical HLA molecule uniquely expressed in the placenta, is a central component of fetus-induced immune tolerance during pregnancy. The tissue-specific expression of HLA-G, however, remains poorly understood. Here, systematic interrogation of the HLA-G locus using massively parallel reporter assay (MPRA) uncovered a previously unidentified cis-regulatory element 12 kb upstream of HLA-G with enhancer activity, Enhancer L Strikingly, clustered regularly-interspaced short palindromic repeats (CRISPR)/Cas9-mediated deletion of this enhancer resulted in ablation of HLA-G expression in JEG3 cells and in primary human trophoblasts isolated from placenta. RNA-seq analysis demonstrated that Enhancer L specifically controls HLA-G expression. Moreover, DNase-seq and chromatin conformation capture (3C) defined Enhancer L as a cell type-specific enhancer that loops into the HLA-G promoter. Interestingly, MPRA-based saturation mutagenesis of Enhancer L identified motifs for transcription factors of the CEBP and GATA families essential for placentation. These factors associate with Enhancer L and regulate HLA-G expression. Our findings identify long-range chromatin looping mediated by core trophoblast transcription factors as the mechanism controlling tissue-specific HLA-G expression at the maternal-fetal interface. More broadly, these results establish the combination of MPRA and CRISPR/Cas9 deletion as a powerful strategy to investigate human immune gene regulation.

Modeling the Effects of Severe Metabolic Disease by Genome Editing of hPSC-Derived Endothelial Cells Reveals an Inflammatory Phenotype.

The kinase AKT2 (PKB) is an important mediator of insulin signaling, for which loss-of-function knockout (KO) mutants lead to early onset diabetes mellitus, and dominant active mutations lead to early development of obesity and endothelial cell (EC) dysfunction. To model EC dysfunction, we used edited human pluripotent stem cells (hPSCs) that carried either a homozygous deletion of AKT2 (AKT2 KO) or a dominant active mutation (AKT2 E17K), which, along with the parental wild type (WT), were differentiated into ECs. Profiling of EC lines indicated an increase in proinflammatory and a reduction in anti-inflammatory fatty acids, an increase in inflammatory chemokines in cell supernatants, increased expression of proinflammatory genes, and increased binding to the EC monolayer in a functional leukocyte adhesion assay for both AKT2 KO and AKT2 E17K. Collectively, these findings suggest that vascular endothelial inflammation that results from dysregulated insulin signaling (homeostasis) may contribute to coronary artery disease, and that either downregulation or upregulation of the insulin pathway may lead to inflammation of endothelial cells. This suggests that the standard of care for patients must be expanded from control of metabolic parameters to include control of inflammation, such that endothelial dysfunction and cardiovascular disorders can ultimately be prevented.

FPLD2 LMNA mutation R482W dysregulates iPSC-derived adipocyte function and lipid metabolism.

Lipodystrophies are disorders that directly affect lipid metabolism and storage. Familial partial lipodystrophy type 2 (FPLD2) is caused by an autosomal dominant mutation in the LMNA gene. FPLD2 is characterized by abnormal adipose tissue distribution. This leads to metabolic deficiencies, such as insulin-resistant diabetes mellitus and hypertriglyceridemia. Here we have derived iPSC lines from two individuals diagnosed with FPLD2, and differentiated these cells into adipocytes. Adipogenesis and certain adipocyte functions are impaired in FPLD2-adipocytes. Consistent with the lipodystrophic phenotype, FPLD2-adipocytes appear to accumulate markers of autophagy and catabolize triglycerides at higher levels than control adipocytes. These data are suggestive of a mechanism causing the lack of adipose tissue in FPLD2 patients.

Genome-Edited Human Pluripotent Stem Cell-Derived Macrophages as a Model of Reverse Cholesterol Transport--Brief Report.

OBJECTIVE: To create isogenic human pluripotent stem cell-derived macrophages with and without ABCA1 expression as a model for reverse cholesterol transport. APPROACH AND RESULTS: The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 (Cas9) genome-editing system was used to introduce frameshift mutations into the coding sequence of ATP-binding cassette, subfamily A, member 1. Individual human pluripotent stem cell clones with deleterious mutations were identified, expanded, and differentiated into mature macrophages with a cytokine-based, feeder-free differentiation protocol. Wild-type cells demonstrated effective cholesterol efflux to apoAI acceptor, whereas ABCA1(-/-) cells displayed significantly reduced efflux ability and increased expression of proinflammatory cytokines. CONCLUSIONS: Human pluripotent stem cell-derived macrophages capable of reverse cholesterol transport can be rapidly generated and genetically edited with CRISPR/Cas9. Introduction of homozygous frameshift mutations results in loss of ABCA1 expression in differentiated macrophages and subsequent reduction of cholesterol efflux capability. This facile genome-editing approach and differentiation protocol pave the way for future studies of the molecular determinants of reverse cholesterol transport and other macrophage properties.

Identification of adult nephron progenitors capable of kidney regeneration in zebrafish.

Loss of kidney function underlies many renal diseases. Mammals can partly repair their nephrons (the functional units of the kidney), but cannot form new ones. By contrast, fish add nephrons throughout their lifespan and regenerate nephrons de novo after injury, providing a model for understanding how mammalian renal regeneration may be therapeutically activated. Here we trace the source of new nephrons in the adult zebrafish to small cellular aggregates containing nephron progenitors. Transplantation of single aggregates comprising 10-30 cells is sufficient to engraft adults and generate multiple nephrons. Serial transplantation experiments to test self-renewal revealed that nephron progenitors are long-lived and possess significant replicative potential, consistent with stem-cell activity. Transplantation of mixed nephron progenitors tagged with either green or red fluorescent proteins yielded some mosaic nephrons, indicating that multiple nephron progenitors contribute to a single nephron. Consistent with this, live imaging of nephron formation in transparent larvae showed that nephrogenic aggregates form by the coalescence of multiple cells and then differentiate into nephrons. Taken together, these data demonstrate that the zebrafish kidney probably contains self-renewing nephron stem/progenitor cells. The identification of these cells paves the way to isolating or engineering the equivalent cells in mammals and developing novel renal regenerative therapies.

Adipocyte insulin receptor activity maintains adipose tissue mass and lifespan.

Type 2 diabetes follows a well-defined progressive pathogenesis, beginning with insulin resistance in metabolic tissues such as the adipose. Intracellular signaling downstream of insulin receptor activation regulates critical metabolic functions of adipose tissue, including glucose uptake, lipogenesis, lipolysis and adipokine secretion. Previous studies have used the aP2 promoter to drive Cre recombinase expression in adipose tissue. Insulin receptor (IR) knockout mice created using this aP2-Cre strategy (FIRKO mice) were protected from obesity and glucose intolerance. Later studies demonstrated the promiscuity of the aP2 promoter, casting doubts upon the tissue specificity of aP2-Cre models. It is our goal to use the increased precision of the Adipoq promoter to investigate adipocyte-specific IR function. Towards this end we generated an adipocyte-specific IR knockout (AIRKO) mouse using an Adipoq-driven Cre recombinase. Here we report AIRKO mice are less insulin sensitive throughout life, and less glucose tolerant than wild-type (WT) littermates at the age of 16 weeks. In contrast to WT littermates, the insulin sensitivity of AIRKO mice is unaffected by age or dietary regimen. At any age, AIRKO mice are comparably insulin resistant to old or obese WT mice and have a significantly reduced lifespan. Similar results were obtained when these phenotypes were re-examined in FIRKO mice. We also found that the AIRKO mouse is protected from high-fat diet-induced weight gain, corresponding with a 90% reduction in tissue weight of major adipose depots compared to WT littermates. Adipose tissue mass reduction is accompanied by hepatomegaly and increased hepatic steatosis. These data indicate that adipocyte IR function is crucial to systemic energy metabolism and has profound effects on adiposity, hepatic homeostasis and lifespan.

Permanent alteration of PCSK9 with in vivo CRISPR-Cas9 genome editing.

RATIONALE: Individuals with naturally occurring loss-of-function proprotein convertase subtilisin/kexin type 9 (PCSK9) mutations experience reduced low-density lipoprotein cholesterol levels and protection against cardiovascular disease. OBJECTIVE: The goal of this study was to assess whether genome editing using a clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated system can efficiently introduce loss-of-function mutations into the endogenous PCSK9 gene in vivo. METHODS AND RESULTS: We used adenovirus to express CRISPR-associated 9 and a CRISPR guide RNA targeting Pcsk9 in mouse liver, where the gene is specifically expressed. We found that <3 to 4 days of administration of the virus, the mutagenesis rate of Pcsk9 in the liver was as high as >50%. This resulted in decreased plasma PCSK9 levels, increased hepatic low-density lipoprotein receptor levels, and decreased plasma cholesterol levels (by 35-40%). No off-target mutagenesis was detected in 10 selected sites. CONCLUSIONS: Genome editing with the CRISPR-CRISPR-associated 9 system disrupts the Pcsk9 gene in vivo with high efficiency and reduces blood cholesterol levels in mice. This approach may have therapeutic potential for the prevention of cardiovascular disease in humans.

From noncoding variant to phenotype via SORT1 at the 1p13 cholesterol locus.

Recent genome-wide association studies (GWASs) have identified a locus on chromosome 1p13 strongly associated with both plasma low-density lipoprotein cholesterol (LDL-C) and myocardial infarction (MI) in humans. Here we show through a series of studies in human cohorts and human-derived hepatocytes that a common noncoding polymorphism at the 1p13 locus, rs12740374, creates a C/EBP (CCAAT/enhancer binding protein) transcription factor binding site and alters the hepatic expression of the SORT1 gene. With small interfering RNA (siRNA) knockdown and viral overexpression in mouse liver, we demonstrate that Sort1 alters plasma LDL-C and very low-density lipoprotein (VLDL) particle levels by modulating hepatic VLDL secretion. Thus, we provide functional evidence for a novel regulatory pathway for lipoprotein metabolism and suggest that modulation of this pathway may alter risk for MI in humans. We also demonstrate that common noncoding DNA variants identified by GWASs can directly contribute to clinical phenotypes.

Rewirable gene regulatory networks in the preimplantation embryonic development of three mammalian species.

Mammalian preimplantation embryonic development (PED) is thought to be governed by highly conserved processes. While it had been suggested that some plasticity of conserved signaling networks exists among different mammalian species, it was not known to what extent modulation of the genomes and the regulatory proteins could "rewire" the gene regulatory networks (GRN) that control PED. We therefore generated global transcriptional profiles from three mammalian species (human, mouse, and bovine) at representative stages of PED, including: zygote, two-cell, four-cell, eight-cell, 16-cell, morula and blastocyst. Coexpression network analysis suggested that 40.2% orthologous gene triplets exhibited different expression patterns among these species. Combining the expression data with genomic sequences and the ChIP-seq data of 16 transcription regulators, we observed two classes of genomic changes that contributed to interspecies expression difference, including single nucleotide mutations leading to turnover of transcription factor binding sites, and insertion of cis-regulatory modules (CRMs) by transposons. About 10% of transposons are estimated to carry CRMs, which may drive species-specific gene expression. The two classes of genomic changes act in concert to drive mouse-specific expression of MTF2, which links POU5F1/NANOG to NOTCH signaling. We reconstructed the transition of the GRN structures as a function of time during PED. A comparison of the GRN transition processes among the three species suggested that in the bovine system, POU5F1's interacting partner SOX2 may be replaced by HMGB1 (a TF sharing the same DNA binding domain with SOX2), resulting in rewiring of GRN by a trans change.

Chromatin regulator SMARCAL1 modulates cellular lipid metabolism.

Biallelic mutations of the chromatin regulator SMARCAL1 cause Schimke Immunoosseous Dysplasia (SIOD), characterized by severe growth defects and premature mortality. Atherosclerosis and hyperlipidemia are common among SIOD patients, yet their onset and progression are poorly understood. Using an integrative approach involving proteomics, mouse models, and population genetics, we investigated SMARCAL1's role. We found that SmarcAL1 interacts with angiopoietin-like 3 (Angptl3), a key regulator of lipoprotein metabolism. In vitro and in vivo analyses demonstrate SmarcAL1's vital role in maintaining cellular lipid homeostasis. The observed translocation of SmarcAL1 to cytoplasmic peroxisomes suggests a potential regulatory role in lipid metabolism through gene expression. SmarcAL1 gene inactivation reduces the expression of key genes in cellular lipid catabolism. Population genetics investigations highlight significant associations between SMARCAL1 genetic variations and body mass index, along with lipid-related traits. This study underscores SMARCAL1's pivotal role in cellular lipid metabolism, likely contributing to the observed lipid phenotypes in SIOD patients.

Human embryonic stem cells: emerging technologies and practical applications.

Human embryonic stem cells possess the unique ability to differentiate into any adult cell type. Recent advances in the understanding of stem cell biology make new applications possible for stem cell based technology. Of note, it is now possible to reprogram terminally differentiated human somatic cells into pluripotent cells that are functionally equivalent to embryonic stem cells. These induced pluripotent cells may become the substrate for future disease models and cell-based therapies. In addition, novel techniques for genetic manipulation have increased the ease with which genes can be modified into stem cells. In this review, we describe these novel technologies as well as developments in the understanding of basic biology of stem cell pluripotency and differentiation.

Genomewide analysis of PRC1 and PRC2 occupancy identifies two classes of bivalent domains.

In embryonic stem (ES) cells, bivalent chromatin domains with overlapping repressive (H3 lysine 27 tri-methylation) and activating (H3 lysine 4 tri-methylation) histone modifications mark the promoters of more than 2,000 genes. To gain insight into the structure and function of bivalent domains, we mapped key histone modifications and subunits of Polycomb-repressive complexes 1 and 2 (PRC1 and PRC2) genomewide in human and mouse ES cells by chromatin immunoprecipitation, followed by ultra high-throughput sequencing. We find that bivalent domains can be segregated into two classes -- the first occupied by both PRC2 and PRC1 (PRC1-positive) and the second specifically bound by PRC2 (PRC2-only). PRC1-positive bivalent domains appear functionally distinct as they more efficiently retain lysine 27 tri-methylation upon differentiation, show stringent conservation of chromatin state, and associate with an overwhelming number of developmental regulator gene promoters. We also used computational genomics to search for sequence determinants of Polycomb binding. This analysis revealed that the genomewide locations of PRC2 and PRC1 can be largely predicted from the locations, sizes, and underlying motif contents of CpG islands. We propose that large CpG islands depleted of activating motifs confer epigenetic memory by recruiting the full repertoire of Polycomb complexes in pluripotent cells.