Publisher Correction: Tet2 and Tet3 in B cells are required to repress CD86 and prevent autoimmunity.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Tet2 and Tet3 in B cells are required to repress CD86 and prevent autoimmunity.

A contribution of epigenetic modifications to B cell tolerance has been proposed but not directly tested. Here we report that deficiency of ten-eleven translocation (Tet) DNA demethylase family members Tet2 and Tet3 in B cells led to hyperactivation of B and T cells, autoantibody production and lupus-like disease in mice. Mechanistically, in the absence of Tet2 and Tet3, downregulation of CD86, which normally occurs following chronic exposure of self-reactive B cells to self-antigen, did not take place. The importance of dysregulated CD86 expression in Tet2- and Tet3-deficient B cells was further demonstrated by the restriction, albeit not complete, on aberrant T and B cell activation following anti-CD86 blockade. Tet2- and Tet3-deficient B cells had decreased accumulation of histone deacetylase 1 (HDAC1) and HDAC2 at the Cd86 locus. Thus, our findings suggest that Tet2- and Tet3-mediated chromatin modification participates in repression of CD86 on chronically stimulated self-reactive B cells, which contributes, at least in part, to preventing autoimmunity.

Highly rigid H3.1/H3.2-H3K9me3 domains set a barrier for cell fate reprogramming in trophoblast stem cells.

The placenta is a highly evolved, specialized organ in mammals. It differs from other organs in that it functions only for fetal maintenance during gestation. Therefore, there must be intrinsic mechanisms that guarantee its unique functions. To address this question, we comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). Our genome-wide, high-throughput analyses revealed that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures with a high degree of histone H3.1/3.2-H3K9me3 accumulation, which we termed TSC-defined highly heterochromatinized domains (THDs). Importantly, depletion of THDs by knockdown of CAF1, an H3.1/3.2 chaperone, resulted in down-regulation of TSC markers, such as Cdx2 and Elf5, and up-regulation of the pluripotent marker Oct3/4, indicating that THDs maintain the trophoblastic nature of TSCs. Furthermore, our nuclear transfer technique revealed that THDs are highly resistant to genomic reprogramming. However, when H3K9me3 was removed, the TSC genome was fully reprogrammed, giving rise to the first TSC cloned offspring. Interestingly, THD-like domains are also present in mouse and human placental cells in vivo, but not in other cell types. Thus, THDs are genomic architectures uniquely developed in placental lineage cells, which serve to protect them from fate reprogramming to stably maintain placental function.

Transient Methionine Deprivation Triggers Histone Modification and Potentiates Differentiation of Induced Pluripotent Stem Cells.

Human induced pluripotent stem cells (iPSCs) require high levels of methionine (Met). Met deprivation results in a rapid decrease in intracellular S-adenosyl-methionine (SAM), poising human iPSCs for differentiation and leading to the apoptosis of undifferentiated cells. Met deprivation triggers rapid metabolic changes, including SAM, followed by reversible epigenetic modifications. Here, we show that short-term Met deprivation impairs the pluripotency network through epigenetic modification in a 3D suspension culture. The trimethylation of lysine 4 on histone H3 (H3K4me3) was drastically affected compared with other histone modifications. Short-term Met deprivation specifically affects the transcription start site (TSS) region of genes, such as those involved in the transforming growth factor ß pathway and cholesterol biosynthetic process, besides key pluripotent genes such as NANOG and POU5F1. The expression levels of these genes decreased, correlating with the loss of H3K4me3 marks. Upon differentiation, Met deprivation triggers the upregulation of various lineage-specific genes, including key definitive endoderm genes, such as GATA6. Upon differentiation, loss of H3K27me3 occurs in many endodermal genes, switching from a bivalent to a monovalent (H3K4me3) state. In conclusion, Met metabolism maintains the pluripotent network with histone marks, and their loss potentiates differentiation.

Codependency and mutual exclusivity for gene community detection from sparse single-cell transcriptome data.

Single-cell RNA-seq (scRNA-seq) can be used to characterize cellular heterogeneity in thousands of cells. The reconstruction of a gene network based on coexpression patterns is a fundamental task in scRNA-seq analyses, and the mutual exclusivity of gene expression can be critical for understanding such heterogeneity. Here, we propose an approach for detecting communities from a genetic network constructed on the basis of coexpression properties. The community-based comparison of multiple coexpression networks enables the identification of functionally related gene clusters that cannot be fully captured through differential gene expression-based analysis. We also developed a novel metric referred to as the exclusively expressed index (EEI) that identifies mutually exclusive gene pairs from sparse scRNA-seq data. EEI quantifies and ranks the exclusive expression levels of all gene pairs from binary expression patterns while maintaining robustness against a low sequencing depth. We applied our methods to glioblastoma scRNA-seq data and found that gene communities were partially conserved after serum stimulation despite a considerable number of differentially expressed genes. We also demonstrate that the identification of mutually exclusive gene sets with EEI can improve the sensitivity of capturing cellular heterogeneity. Our methods complement existing approaches and provide new biological insights, even for a large, sparse dataset, in the single-cell analysis field.

Missense Mutations in NKAP Cause a Disorder of Transcriptional Regulation Characterized by Marfanoid Habitus and Cognitive Impairment.

NKAP is a ubiquitously expressed nucleoplasmic protein that is currently known as a transcriptional regulatory molecule via its interaction with HDAC3 and spliceosomal proteins. Here, we report a disorder of transcriptional regulation due to missense mutations in the X chromosome gene, NKAP. These mutations are clustered in the C-terminal region of NKAP where NKAP interacts with HDAC3 and post-catalytic spliceosomal complex proteins. Consistent with a role for the C-terminal region of NKAP in embryogenesis, nkap mutant zebrafish with a C-terminally truncated NKAP demonstrate severe developmental defects. The clinical features of affected individuals are highly conserved and include developmental delay, hypotonia, joint contractures, behavioral abnormalities, Marfanoid habitus, and scoliosis. In affected cases, transcriptome analysis revealed the presence of a unique transcriptome signature, which is characterized by the downregulation of long genes with higher exon numbers. These observations indicate the critical role of NKAP in transcriptional regulation and demonstrate that perturbations of the C-terminal region lead to developmental defects in both humans and zebrafish.

Single-cell transcriptional analysis reveals developmental stage-dependent changes in retinal progenitors in the murine early optic vesicle.

The optic vesicle in the developing embryonic eye contains a multitude of neuroepithelial progenitors that subsequently differentiate into functionally distinct domains of the optic cup, such as the neural retina, pigment epithelium, and optic stalk. To investigate cell-type diversity across early optic vesicles before regionalization of the optic cup, we performed single-cell RNA-sequencing (scRNA-seq) using 7989 cells from the presumptive eye area in mouse embryos at the 12-26-somite stages at five developmental time points. We demonstrated the presence of seven optic vesicle populations. Moreover, the four populations of retinal progenitor cells could be classified according to their stage-dependent time point, and these cells exhibited altered expression of several structural and metabolic key genes, such as Col9a1 and Ckb, just before regionalization of the optic cup. From these data, we provide the first report on stage-dependent transcriptional profiles during initial retinal specification at single-cell resolution and highlight the unexpected developmental heterogeneity of the murine optic vesicle structure.

Identification of Elg1 interaction partners and effects on post-replication chromatin re-formation.

Elg1, the major subunit of a Replication Factor C-like complex, is critical to ensure genomic stability during DNA replication, and is implicated in controlling chromatin structure. We investigated the consequences of Elg1 loss for the dynamics of chromatin re-formation following DNA replication. Measurement of Okazaki fragment length and the micrococcal nuclease sensitivity of newly replicated DNA revealed a defect in nucleosome organization in the absence of Elg1. Using a proteomic approach to identify Elg1 binding partners, we discovered that Elg1 interacts with Rtt106, a histone chaperone implicated in replication-coupled nucleosome assembly that also regulates transcription. A central role for Elg1 is the unloading of PCNA from chromatin following DNA replication, so we examined the relative importance of Rtt106 and PCNA unloading for chromatin reassembly following DNA replication. We find that the major cause of the chromatin organization defects of an ELG1 mutant is PCNA retention on DNA following replication, with Rtt106-Elg1 interaction potentially playing a contributory role.

ARCN1 Mutations Cause a Recognizable Craniofacial Syndrome Due to COPI-Mediated Transport Defects.

Cellular homeostasis is maintained by the highly organized cooperation of intracellular trafficking systems, including COPI, COPII, and clathrin complexes. COPI is a coatomer protein complex responsible for intracellular protein transport between the endoplasmic reticulum and the Golgi apparatus. The importance of such intracellular transport mechanisms is underscored by the various disorders, including skeletal disorders such as cranio-lenticulo-sutural dysplasia and osteogenesis imperfect, caused by mutations in the COPII coatomer complex. In this article, we report a clinically recognizable craniofacial disorder characterized by facial dysmorphisms, severe micrognathia, rhizomelic shortening, microcephalic dwarfism, and mild developmental delay due to loss-of-function heterozygous mutations in ARCN1, which encodes the coatomer subunit delta of COPI. ARCN1 mutant cell lines were revealed to have endoplasmic reticulum stress, suggesting the involvement of ER stress response in the pathogenesis of this disorder. Given that ARCN1 deficiency causes defective type I collagen transport, reduction of collagen secretion represents the likely mechanism underlying the skeletal phenotype that characterizes this condition. Our findings demonstrate the importance of COPI-mediated transport in human development, including skeletogenesis and brain growth.

Elevation of insulin-like growth factor binding protein-2 level in Pallister-Killian syndrome: implications for the postnatal growth retardation phenotype.

Pallister-Killian syndrome (PKS) is a multi-system developmental disorder caused by tetrasomy 12p that exhibits tissue-limited mosaicism. Probands with PKS often demonstrate a unique growth profile consisting of macrosomia at birth with deceleration of growth postnatally. We have previously demonstrated that cultured skin fibroblasts from PKS probands have significantly elevated expression of insulin-like growth factor binding protein-2 (IGFBP2). To further evaluate the role of IGFBP2 in PKS, the amount of IGFBP2 secreted from cultured skin fibroblast cell lines and serum IGFBP2 levels were measured in probands with PKS. Approximately 60% of PKS fibroblast cell lines secreted higher levels of IGFBP2 compared to control fibroblasts, although the remaining 40% of PKS samples produced comparable level of IGFBP2 to that of control fibroblasts. Serum IGFBP2 levels were also measured in PKS probands and were elevated in 40% of PKS probands. PKS probands with elevated IGFBP2 manifested with severe postnatal growth retardation. IGFBPs are the family of related proteins that bind IGFs with high affinity and are typically thought to attenuate IGF action. We suggest that elevated IGFBP2 levels might play a role in the growth retardation phenotype of PKS.

Genome-wide analysis of murine renal distal convoluted tubular cells for the target genes of mineralocorticoid receptor.

BACKGROUND AND OBJECTIVE: Mineralocorticoid receptor (MR) is a member of nuclear receptor family proteins and contributes to fluid homeostasis in the kidney. Although aldosterone-MR pathway induces several gene expressions in the kidney, it is often unclear whether the gene expressions are accompanied by direct regulations of MR through its binding to the regulatory region of each gene. The purpose of this study is to identify the direct target genes of MR in a murine distal convoluted tubular epithelial cell-line (mDCT). METHODS: We analyzed the DNA samples of mDCT cells overexpressing 3xFLAG-hMR after treatment with 10(-7)M aldosterone for 1h by chromatin immunoprecipitation with deep-sequence (ChIP-seq) and mRNA of the cell-line with treatment of 10(-7)M aldosterone for 3h by microarray. RESULTS: 3xFLAG-hMR overexpressed in mDCT cells accumulated in the nucleus in response to 10(-9)M aldosterone. Twenty-five genes were indicated as the candidate target genes of MR by ChIP-seq and microarray analyses. Five genes, Sgk1, Fkbp5, Rasl12, Tns1 and Tsc22d3 (Gilz), were validated as the direct target genes of MR by quantitative RT-qPCR and ChIP-qPCR. MR binding regions adjacent to Ctgf and Serpine1 were also validated. CONCLUSIONS: We, for the first time, captured the genome-wide distribution of MR in mDCT cells and, furthermore, identified five MR target genes in the cell-line. These results will contribute to further studies on the mechanisms of kidney diseases.

Cell cycle-dependent gene networks for cell proliferation activated by nuclear CK2α complexes.

Nuclear expression of protein kinase CK2α is reportedly elevated in human carcinomas, but mechanisms underlying its variable localization in cells are poorly understood. This study demonstrates a functional connection between nuclear CK2 and gene expression in relation to cell proliferation. Growth stimulation of quiescent human normal fibroblasts and phospho-proteomic analysis identified a pool of CK2α that is highly phosphorylated at serine 7. Phosphorylated CK2α translocates into the nucleus, and this phosphorylation appears essential for nuclear localization and catalytic activity. Protein signatures associated with nuclear CK2 complexes reveal enrichment of apparently unique transcription factors and chromatin remodelers during progression through the G1 phase of the cell cycle. Chromatin immunoprecipitation-sequencing profiling demonstrated recruitment of CK2α to active gene loci, more abundantly in late G1 phase than in early G1, notably at transcriptional start sites of core histone genes, growth stimulus-associated genes, and ribosomal RNAs. Our findings reveal that nuclear CK2α complexes may be essential to facilitate progression of the cell cycle, by activating histone genes and triggering ribosomal biogenesis, specified in association with nuclear and nucleolar transcriptional regulators.

HSF1 phosphorylation establishes an active chromatin state via the TRRAP-TIP60 complex and promotes tumorigenesis.

Transcriptional regulation by RNA polymerase II is associated with changes in chromatin structure. Activated and promoter-bound heat shock transcription factor 1 (HSF1) recruits transcriptional co-activators, including histone-modifying enzymes; however, the mechanisms underlying chromatin opening remain unclear. Here, we demonstrate that HSF1 recruits the TRRAP-TIP60 acetyltransferase complex in HSP72 promoter during heat shock in a manner dependent on phosphorylation of HSF1-S419. TRIM33, a bromodomain-containing ubiquitin ligase, is then recruited to the promoter by interactions with HSF1 and a TIP60-mediated acetylation mark, and cooperates with the related factor TRIM24 for mono-ubiquitination of histone H2B on K120. These changes in histone modifications are triggered by phosphorylation of HSF1-S419 via PLK1, and stabilize the HSF1-transcription complex in HSP72 promoter. Furthermore, HSF1-S419 phosphorylation is constitutively enhanced in and promotes proliferation of melanoma cells. Our results provide mechanisms for HSF1 phosphorylation-dependent establishment of an active chromatin status, which is important for tumorigenesis.

MAB21L1 modulates gene expression and DNA metabolic processes in the lens placode.

Mutations in human MAB21L1 cause aberrations in lens ectoderm morphogenesis and lead to congenital cerebellar, ocular, craniofacial and genital (COFG) syndrome. Murine Mab21l1-null mutations cause severe cell-autonomous defects in lens formation, leading to microphthalmia; therefore, Mab21l1-null mice are used as a mouse model for COFG syndrome. In this study, we investigated the early-onset single-cell-level phenotypes of murine Mab21l1-null lens ectoderms using electron microscopy and single-cell RNA sequencing (scRNA-seq). Electron microscopy and immunohistochemical analyses indicated endoplasmic reticulum stress at the 24- to 26-somite stage in Mab21l1-null lens placodes. scRNA-seq analysis revealed that 131 genes were downregulated and 148 were upregulated in Mab21l1-null lens ectoderms relative to the wild type. We successfully identified 21 lens-specific genes that were downregulated in Mab21l1-null cells, including three key genes involved in lens formation: Pitx3, Maf and Sfrp2. Moreover, gene ontology analysis of the 279 differentially expressed genes indicated enrichment in housekeeping genes associated with DNA/nucleotide metabolism prior to cell death. These findings suggest that MAB21L1 acts as a nuclear factor that modulates not only lens-specific gene expression but also DNA/nucleotide metabolic processes during lens placode formation.

Small extracellular vesicles derived from interferon-γ pre-conditioned mesenchymal stromal cells effectively treat liver fibrosis.

Mesenchymal stromal cells (MSCs) are used for ameliorating liver fibrosis and aiding liver regeneration after cirrhosis; Here, we analyzed the therapeutic potential of small extracellular vesicles (sEVs) derived from interferon-γ (IFN-γ) pre-conditioned MSCs (γ-sEVs). γ-sEVs effectively induced anti-inflammatory macrophages with high motility and phagocytic abilities in vitro, while not preventing hepatic stellate cell (HSC; the major source of collagen fiber) activation in vitro. The proteome analysis of MSC-derived sEVs revealed anti-inflammatory macrophage inducible proteins (e.g., annexin-A1, lactotransferrin, and aminopeptidase N) upon IFN-γ stimulation. Furthermore, by enabling CX3CR1+ macrophage accumulation in the damaged area, γ-sEVs ameliorated inflammation and fibrosis in the cirrhosis mouse model more effectively than sEVs. Single cell RNA-Seq analysis revealed diverse effects, such as induction of anti-inflammatory macrophages and regulatory T cells, in the cirrhotic liver after γ-sEV administration. Overall, IFN-γ pre-conditioning altered sEVs resulted in efficient tissue repair indicating a new therapeutic strategy.

Expression of the peroxisome proliferator activated receptor gamma gene is repressed by DNA methylation in visceral adipose tissue of mouse models of diabetes.

BACKGROUND: Adipose tissues serve not only as a store for energy in the form of lipid, but also as endocrine tissues that regulates metabolic activities of the organism by secreting various kinds of hormones. Peroxisome proliferator activated receptor gamma (PPARgamma) is a key regulator of adipocyte differentiation that induces the expression of adipocyte-specific genes in preadipocytes and mediates their differentiation into adipocytes. Furthermore, PPARgamma has an important role to maintain the physiological function of mature adipocyte by controlling expressions of various genes properly. Therefore, any reduction in amount and activity of PPARgamma is linked to the pathogenesis of metabolic syndrome. RESULTS: In this study, we investigated the contribution of epigenetic transcriptional regulatory mechanisms, such as DNA methylation, to the expression of the PPARgamma gene, and further evaluated the contribution of such epigenetic regulatory mechanisms to the pathogenesis of metabolic syndrome. In 3T3-L1 preadipocytes, the promoter of the PPARgamma2 gene was hypermethylated, but was progressively demethylated upon induction of differentiation, which was accompanied by an increase of mRNA expression. Moreover, treatment of cells with 5'-aza-cytideine, an inhibitor of DNA methylation, increased expression of the PPARgamma gene in a dose-dependent manner. Methylation in vitro of a PPARgamma promoter-driven reporter construct also repressed the transcription of a downstream reporter gene. These results suggest that the expression of the PPARgamma gene is inhibited by methylation of its promoter. We next compared the methylation status of the PPARgamma promoters in adipocytes from wild-type (WT) mice with those from two diabetic mouse models: +Leprdb/+Leprdb and diet-induced obesity mice. Interestingly, we found increased methylation of the PPARgamma promoter in visceral adipose tissues (VAT) of the mouse models of diabetes, compared to that observed in wild-type mice. We observed a concomitant decrease in the level of PPARgamma mRNA in the diabetic mice compared to the WT mice. CONCLUSION: We conclude that the expression of PPARgamma gene is regulated by DNA methylation of its promoter region and propose that reduced expression of PPARgamma owing to DNA methylation in adipocytes of the VAT may contribute to the pathogenesis of metabolic syndrome.

Age-Dependent Ribosomal DNA Variations in Mice.

The rRNA gene, which consists of tandem repetitive arrays (ribosomal DNA [rDNA] repeat), is one of the most unstable regions in the genome. The rDNA repeat in the budding yeast Saccharomyces cerevisiae is known to become unstable as the cell ages. However, it is unclear how the rDNA repeat changes in aging mammalian cells. Using quantitative single-cell analyses, we identified age-dependent alterations in rDNA copy number and levels of methylation in mice. The degree of methylation and copy number of rDNA from bone marrow cells of 2-year-old mice were increased by comparison to levels in 4-week-old mice in two mouse strains, BALB/cA and C57BL/6. Moreover, the level of pre-rRNA transcripts was reduced in older BALB/cA mice. We also identified many sequence variations in the rDNA. Among them, three mutations were unique to old mice, and two of them were found in the conserved region in budding yeast. We established yeast strains with the old-mouse-specific mutations and found that they shortened the life span of the cells. Our findings suggest that rDNA is also fragile in mammalian cells and that alterations within this region have a profound effect on cellular function.

Mosaic ratio quantification of isochromosome 12p in Pallister-Killian syndrome using droplet digital PCR.

BACKGROUND: Pallister-Killian syndrome (PKS) is a prototypic mosaic aneuploidy syndrome caused by mosaic supernumerary marker isochromosome 12p. Cells possessing the isochromosome 12p rapidly diminish after birth in the peripheral blood, often necessitating a skin biopsy for diagnosis. Therefore, a genomic testing that is capable of detecting low percent mosaic isochromosome 12p is preferred for the diagnosis of PKS. METHODS: The utility of the droplet digital PCR system in quantifying the mosaic ratio of isochromosome 12p in PKS was evaluated. RESULTS: Droplet digital PCR was able to precisely quantify isochromosome 12p mosaic ratio, and copy number measured by droplet digital PCR was correlated well with that of fluorescence in situ hybridization analysis. CONCLUSION: Droplet digital PCR should be considered as an effective tool for both clinical and research analytics to precisely quantify mosaic genomic copy number alterations or mosaic mutations.

PPARγ-induced PARylation promotes local DNA demethylation by production of 5-hydroxymethylcytosine.

Recent studies have shown that DNA demethylation goes through the conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) by Tet proteins. However, it is still unclear how the target regions for demethylation are distinguished within their genomic context. Here we show that the nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ) has the ability to direct local demethylation around its binding sites, the PPAR response elements (PPREs), during adipocyte differentiation. PPARγ is a key regulator of the differentiation process that forms a PPARγ co-activator complex on PPREs and activates the expression of adipocyte-specific genes. The complex is poly(ADP-ribosyl)ated (PARylated) on PPREs, and Tet proteins catalyse the conversion of 5mC to 5hmC locally by their ability to bind to the PAR polymer, thereby inducing region-specific demethylation. Our study demonstrates that a sequence-dependent transcription factor complex can, through its post-translational modification, serve for Tet proteins as a landmark to identify sites of DNA demethylation.