Transcription factor-mediated direct cellular reprogramming yields cell-type specific DNA methylation signature.

Direct reprogramming, inducing the conversion of one type of somatic cell into another by the forced expression of defined transcription factors, is a technology with anticipated medical applications. However, due to the many unresolved aspects of the induction mechanisms, it is essential to thoroughly analyze the epigenomic state of the generated cells. Here, we performed comparative genome-wide DNA methylation analyses of mouse embryonic fibroblasts (MEFs) and cells composing organoids formed by intestinal stem cells (ISCs) or induced ISCs (iISCs) that were directly induced from MEFs. We found that the CpG methylation state was similar between cells forming ISC organoids and iISC organoids, while they differed widely from those in MEFs. Moreover, genomic regions that were differentially methylated between ISC organoid- and iISC organoid-forming cells did not significantly affect gene expression. These results demonstrate the accuracy and safety of iISC induction, leading to the medical applications of this technology.

DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins in liver and skeletal muscle.

Each tissue has a dominant set of functional proteins required to mediate tissue-specific functions. Epigenetic modifications, transcription, and translational efficiency control tissue-dominant protein production. However, the coordination of these regulatory mechanisms to achieve such tissue-specific protein production remains unclear. Here, we analyzed the DNA methylome, transcriptome, and proteome in mouse liver and skeletal muscle. We found that DNA hypomethylation at promoter regions is globally associated with liver-dominant or skeletal muscle-dominant functional protein production within each tissue, as well as with genes encoding proteins involved in ubiquitous functions in both tissues. Thus, genes encoding liver-dominant proteins, such as those involved in glycolysis or gluconeogenesis, the urea cycle, complement and coagulation systems, enzymes of tryptophan metabolism, and cytochrome P450-related metabolism, were hypomethylated in the liver, whereas those encoding-skeletal muscle-dominant proteins, such as those involved in sarcomere organization, were hypomethylated in the skeletal muscle. Thus, DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins.

Persistent epigenetic alterations in transcription factors after a sustained virological response in hepatocellular carcinoma.

Background and Aim: The risk of hepatocellular carcinoma (HCC) persists in a condition of sustained virologic response (SVR) after hepatitis C virus (HCV) eradication. Comprehensive molecular analyses were performed to test the hypothesis that epigenetic abnormalities present after an SVR play a role in hepatocarcinogenesis. Methods: Whole-genome methylome and RNA sequencing were performed on HCV, SVR, and healthy liver tissue. Integrated analysis of the sequencing data focused on expression changes in transcription factors and their target genes, commonly found in HCV and SVR. Identified expression changes were validated in demethylated cultured HCC cell lines and an independent validation cohort. Results: The coincidence rates of the differentially methylated regions between the HCV and SVR groups were 91% in the hypomethylated and 71% in the hypermethylated regions in tumorous tissues, and 37% in the hypomethylated and 36% in the hypermethylated regions in non-tumorous tissues. These results indicate that many epigenomic abnormalities persist even after an SVR was achieved. Integrated analysis identified 61 transcription factors and 379 other genes that had methylation abnormalities and gene expression changes in both groups. Validation cohort specified gene expression changes for 14 genes, and gene ontology pathway analysis revealed apoptotic signaling and inflammatory response were associated with these genes. Conclusion: This study demonstrates that DNA methylation abnormalities, retained after HCV eradication, affect the expression of transcription factors and their target genes. These findings suggest that DNA methylation in SVR patients may be functionally important in carcinogenesis, and could serve as biomarkers to predict HCC occurrence.

Pan-cancer methylome analysis for cancer diagnosis and classification of cancer cell of origin.

The accurate and early diagnosis and classification of cancer origin from either tissue or liquid biopsy is crucial for selecting the appropriate treatment and reducing cancer-related mortality. Here, we established the CAncer Cell-of-Origin (CACO) methylation panel using the methylation data of the 28 types of cancer in The Cancer Genome Atlas (7950 patients and 707 normal controls) as well as healthy whole blood samples (95 subjects). We showed that the CACO methylation panel had high diagnostic potential with high sensitivity and specificity in the discovery (maximum AUC = 0.998) and validation (maximum AUC = 1.000) cohorts. Moreover, we confirmed that the CACO methylation panel could identify the cancer cell type of origin using the methylation profile from liquid as well as tissue biopsy, including primary, metastatic, and multiregional cancer samples and cancer of unknown primary, independent of the methylation analysis platform and specimen preparation method. Together, the CACO methylation panel can be a powerful tool for the classification and diagnosis of cancer.

Transcriptome dynamics during cholesterol-induced transdifferentiation of human coronary artery smooth muscle cells: A Gene Ontology-centric clustering approach.

Macrophage-like cells derived from vascular smooth muscle cells (SMCs) play critical roles in atherogenesis, and DNA hydroxymethylation was implicated in transdifferentiation. We examined transcriptomes and (hydroxy)methylomes of human coronary artery SMCs during cholesterol-induced transdifferentiation. A unique approach of exhaustive identification of differentially expressed genes followed by Gene Ontology-centric clustering facilitated deeper understanding of multifaceted modulations of genes involved in extracellular matrix organization, angiogenesis, cell migration, hypoxia response, and cholesterol biosynthesis. Intriguingly, type I interferon response was transiently activated, presumably forming an immuno-metabolic circuit with cholesterol metabolism. Neither global nor DEG-proximal changes were evident in (hydroxy)methylation. These results would not only provide a unique data resource for atherosclerosis research but present a potentially useful approach in transcriptome data interpretation.

Hoxa10 mediates positional memory to govern stem cell function in adult skeletal muscle.

Muscle stem cells (satellite cells) are distributed throughout the body and have heterogeneous properties among muscles. However, functional topographical genes in satellite cells of adult muscle remain unidentified. Here, we show that expression of Homeobox-A (Hox-A) cluster genes accompanied with DNA hypermethylation of the Hox-A locus was robustly maintained in both somite-derived muscles and their associated satellite cells in adult mice, which recapitulates their embryonic origin. Somite-derived satellite cells were clearly separated from cells derived from cranial mesoderm in Hoxa10 expression. Hoxa10 inactivation led to genomic instability and mitotic catastrophe in somite-derived satellite cells in mice and human. Satellite cell-specific Hoxa10 ablation in mice resulted in a decline in the regenerative ability of somite-derived muscles, which were unobserved in cranial mesoderm-derived muscles. Thus, our results show that Hox gene expression profiles instill the embryonic history in satellite cells as positional memory, potentially modulating region-specific pathophysiology in adult muscles.

Reconstitution of the oocyte transcriptional network with transcription factors.

During female germline development, oocytes become a highly specialized cell type and form a maternal cytoplasmic store of crucial factors. Oocyte growth is triggered at the transition from primordial to primary follicle and is accompanied by dynamic changes in gene expression1, but the gene regulatory network that controls oocyte growth remains unknown. Here we identify a set of transcription factors that are sufficient to trigger oocyte growth. By investigation of the changes in gene expression and functional screening using an in vitro mouse oocyte development system, we identified eight transcription factors, each of which was essential for the transition from primordial to primary follicle. Notably, enforced expression of these transcription factors swiftly converted pluripotent stem cells into oocyte-like cells that were competent for fertilization and subsequent cleavage. These transcription-factor-induced oocyte-like cells were formed without specification of primordial germ cells, epigenetic reprogramming or meiosis, and demonstrate that oocyte growth and lineage-specific de novo DNA methylation are separable from the preceding epigenetic reprogramming in primordial germ cells. This study identifies a core set of transcription factors for orchestrating oocyte growth, and provides an alternative source of ooplasm, which is a unique material for reproductive biology and medicine.

Identification of genes associated with endometrial cell ageing.

Ageing of the uterine endometrium is a critical factor that affects reproductive success, but the mechanisms associated with uterine ageing are unclear. In this study, we conducted a qualitative examination of age-related changes in endometrial tissues and identified candidate genes as markers for uterine ageing. Gene expression patterns were assessed by two RNA-sequencing experiments using uterine tissues from wild type (WT) C57BL/6 mice. Gene expression data obtained by RNA-sequencing were validated by real-time PCR. Genes expressing the pro-inflammatory cytokines Il17rb and chemokines Cxcl12 and Cxcl14 showed differential expression between aged WT mice and a group of mice composed of 5- and 8-week-old WT (young) animals. Protein expression levels of the above-mentioned genes and of IL8, which functions downstream of IL17RB, were analysed by quantitative immunohistochemistry of unaffected human endometrium tissue samples from patients in their 20s and 40s (10 cases each). In the secretory phase samples, 3,3'- diaminobenzidine staining intensities of IL17RB, CXCL12 and CXCL14 for patients in their 40s were significantly higher than that for patients in their 20s, as detected by a Mann-hitney U test. These results suggest that these genes are candidate markers for endometrial ageing and for prediction of age-related infertility, although confirmation of these findings is needed in larger studies involving fertile and infertile women.

Base-resolution methylomes of gliomas bearing histone H3.3 mutations reveal a G34 mutant-specific signature shared with bone tumors.

Two recurrent mutations, K27M and G34R/V, in H3F3A, encoding non-canonical histone H3.3, are reported in pediatric and young adult gliomas, whereas G34W mutation is prevalent in bone tumors. In contrast to K27M mutation, it remains elusive how G34 mutations affect the epigenome. Here we performed whole-genome bisulfite sequencing of four G34R-mutated gliomas and the G34V-mutated glioma cell line KNS-42 for comparison with gliomas harboring K27M and no mutations in H3F3A and with G34W-mutated bone tumors. G34R-mutated gliomas exhibited lower global methylation levels, similar CpG island (CGI) methylation levels, and compromised hypermethylation of telomere-proximal CGIs, compared to the other two glioma subgroups. Hypermethylated regions specific to G34R-mutated gliomas were enriched for CGIs, including those of OLIG1, OLIG2, and canonical histone genes in the HIST1 cluster. They were notably hypermethylated in osteosarcomas with, but not without, G34W mutation. Independent component analysis revealed that G34 mutation-specific components shared a significant similarity between glioma and osteosarcoma, suggesting that G34 mutations exert characteristic methylomic effects regardless of the tumor tissue-of-origin. CRISPR/Cas9-mediated disruption of G34V-allele in KNS-42 cells led to demethylation of a subset of CGIs hypermethylated in G34R-mutated gliomas. These findings will provide a basis for elucidating epigenomic roles of G34 oncohistone in tumorigenesis.

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.

ATF7-Dependent Epigenetic Changes Are Required for the Intergenerational Effect of a Paternal Low-Protein Diet.

Paternal dietary conditions may contribute to metabolic disorders in offspring. We have analyzed the role of the stress-dependent epigenetic regulator cyclic AMP-dependent transcription factor 7 (ATF7) in paternal low-protein diet (pLPD)-induced gene expression changes in mouse liver. Atf7+/- mutations cause an offspring phenotype similar to that caused by pLPD, and the effect of pLPD almost vanished when paternal Atf7+/- mice were used. ATF7 binds to the promoter regions of ∼2,300 genes, including cholesterol biosynthesis-related and tRNA genes in testicular germ cells (TGCs). LPD induces ATF7 phosphorylation by p38 via reactive oxygen species (ROS) in TGCs. This leads to the release of ATF7 and a decrease in histone H3K9 dimethylation (H3K9me2) on its target genes. These epigenetic changes are maintained and induce expression of some tRNA fragments in spermatozoa. These results indicate that LPD-induced and ATF7-dependent epigenetic changes in TGCs play an important role in paternal diet-induced metabolic reprograming in offspring.

Base-Resolution Methylome of Retinal Pigment Epithelial Cells Used in the First Trial of Human Induced Pluripotent Stem Cell-Based Autologous Transplantation.

The first-in-human trial of induced pluripotent stem cell (iPSC)-based autologous transplantation was successfully performed on a female patient with age-related macular degeneration. Here we delineated the base-resolution methylome of the iPSC-derived retinal pigment epithelium (iRPE) used in this trial. The methylome of iRPE closely resembled that of native RPE (nRPE), although partially methylated domains (PMDs) emerged in iRPE but not nRPE. Most differentially methylated regions between iRPE and nRPE appeared to originate from (de)methylation errors during differentiation, whereas errors at reprogramming resulted in aberrant genomic imprinting and X chromosome reactivation. Moreover, non-CpG methylation was prominent in nRPE but not iRPE. Intriguingly, xenotransplantation to mouse remodeled the iRPE methylome to demethylate a subset of suppressed genes and accumulate non-CpG methylation, but failed to resolve PMDs and hypermethylated CpG islands. Although the impacts of these alterations remain elusive, our findings should provide a useful guide for methylome analyses of other iPSC-derived cells.

Highly efficient single-stranded DNA ligation technique improves low-input whole-genome bisulfite sequencing by post-bisulfite adaptor tagging.

Whole-genome bisulfite sequencing (WGBS) is the current gold standard of methylome analysis. Post-bisulfite adaptor tagging (PBAT) is an increasingly popular WGBS protocol because of high sensitivity and low bias. PBAT originally relied on two rounds of random priming for adaptor-tagging of single-stranded DNA (ssDNA) to attain high efficiency but at a cost of library insert length. To overcome this limitation, we developed terminal deoxyribonucleotidyl transferase (TdT)-assisted adenylate connector-mediated ssDNA (TACS) ligation as an alternative to random priming. In this method, TdT attaches adenylates to the 3'-end of input ssDNA, which are then utilized by RNA ligase as an efficient connector to the ssDNA adaptor. A protocol that uses TACS ligation instead of the second random priming step substantially increased the lengths of PBAT library fragments. Moreover, we devised a dual-library strategy that splits the input DNA to prepare two libraries with reciprocal adaptor polarity, combining them prior to sequencing. This strategy ensured an ideal base-color balance to eliminate the need for DNA spike-in for color compensation, further improving the throughput and quality of WGBS. Adopting the above strategies to the HiSeq X Ten and NovaSeq 6000 platforms, we established a cost-effective, high-quality WGBS, which should accelerate various methylome analyses.

Mapping of histone-binding sites in histone replacement-completed spermatozoa.

The majority of histones are replaced by protamines during spermatogenesis, but small amounts are retained in mammalian spermatozoa. Since nucleosomes in spermatozoa influence epigenetic inheritance, it is important to know how histones are distributed in the sperm genome. Conflicting data, which may result from different conditions used for micrococcal nuclease (MNase) digestion, have been reported: retention of nucleosomes at either gene promoter regions or within distal gene-poor regions. Here, we find that the swim-up sperm used in many studies contain about 10% population of sperm which have not yet completed the histone-to-protamine replacement. We develop a method to purify histone replacement-completed sperm (HRCS) and to completely solubilize histones from cross-linked HRCS without MNase digestion. Our results indicate that histones are retained at specific promoter regions in HRCS. This method allows the study of epigenetic status in mature sperm.

The AP-1 transcription factor JunB is required for Th17 cell differentiation.

Interleukin (IL)-17-producing T helper (Th17) cells are crucial for host defense against extracellular microbes and pathogenesis of autoimmune diseases. Here we show that the AP-1 transcription factor JunB is required for Th17 cell development. Junb-deficient CD4+ T cells are able to develop in vitro into various helper T subsets except Th17. The RNA-seq transcriptome analysis reveals that JunB is crucial for the Th17-specific gene expression program. Junb-deficient mice are completely resistant to experimental autoimmune encephalomyelitis, a Th17-mediated inflammatory disease, and naive T helper cells from such mice fail to differentiate into Th17 cells. JunB appears to activate Th17 signature genes by forming a heterodimer with BATF, another AP-1 factor essential for Th17 differentiation. The mechanism whereby JunB controls Th17 cell development likely involves activation of the genes for the Th17 lineage-specifying orphan receptors RORγt and RORα and reduced expression of Foxp3, a transcription factor known to antagonize RORγt function.

Enhanced Stomatal Conductance by a Spontaneous Arabidopsis Tetraploid, Me-0, Results from Increased Stomatal Size and Greater Stomatal Aperture.

The rate of gas exchange in plants is regulated mainly by stomatal size and density. Generally, higher densities of smaller stomata are advantageous for gas exchange; however, it is unclear what the effect of an extraordinary change in stomatal size might have on a plant's gas-exchange capacity. We investigated the stomatal responses to CO2 concentration changes among 374 Arabidopsis (Arabidopsis thaliana) ecotypes and discovered that Mechtshausen (Me-0), a natural tetraploid ecotype, has significantly larger stomata and can achieve a high stomatal conductance. We surmised that the cause of the increased stomatal conductance is tetraploidization; however, the stomatal conductance of another tetraploid accession, tetraploid Columbia (Col), was not as high as that in Me-0. One difference between these two accessions was the size of their stomatal apertures. Analyses of abscisic acid sensitivity, ion balance, and gene expression profiles suggested that physiological or genetic factors restrict the stomatal opening in tetraploid Col but not in Me-0. Our results show that Me-0 overcomes the handicap of stomatal opening that is typical for tetraploids and achieves higher stomatal conductance compared with the closely related tetraploid Col on account of larger stomatal apertures. This study provides evidence for whether larger stomatal size in tetraploids of higher plants can improve stomatal conductance.

Haptoglobin promoter polymorphism rs5472 as a prognostic biomarker for peptide vaccine efficacy in castration-resistant prostate cancer patients.

Personalized peptide vaccination (PPV) is an attractive approach to cancer immunotherapy with strong immune-boosting effects conferring significant clinical benefit. However, as with most therapeutic agents, there is a difference in clinical efficacy among patients receiving PPV. Therefore, a useful biomarker is urgently needed for prognosticating clinical outcomes to preselect patients who would benefit the most from PPV. In this retrospective study, to detect a molecular prognosticator of clinical outcomes for PPV, we analyzed whole-genome gene expression profiles of peripheral blood mononuclear cells (PBMCs) in castration-resistant prostate cancer (CRPC) patients before administration of PPV. Cox regression analysis revealed that mRNA expression of myeloperoxidase, haptoglobin, and neutrophil elastase was significantly associated with overall survival (OS) among vaccinated CRPC patients (adjusted P < 0.01). By promoter sequence analysis of these three genes, we found that rs5472 of haptoglobin (HP), an acute-phase plasma glycoprotein, was strongly correlated to OS of vaccinated CRPC patients (P = 0.0047, hazard ratio 0.47; 95 % confidence interval 0.28-0.80). Furthermore, both HP mRNA expression in PBMCs and protein level in plasma of CRPC patients before administration of PPV exhibited rs5472 dependence (P < 0.001 for mRNA expression and P < 0.05 for protein level). Our findings suggest that rs5472 may play an important role in the immune response to PPV via regulation of HP. Thus, we concluded that rs5472 is a potential prognostic biomarker for PPV.

Changepoint detection in base-resolution methylome data reveals a robust signature of methylated domain landscape.

BACKGROUND: Base-resolution methylome data generated by whole-genome bisulfite sequencing (WGBS) is often used to segment the genome into domains with distinct methylation levels. However, most segmentation methods include many parameters to be carefully tuned and/or fail to exploit the unsurpassed resolution of the data. Furthermore, there is no simple method that displays the composition of the domains to grasp global trends in each methylome. RESULTS: We propose to use changepoint detection for domain demarcation based on base-resolution methylome data. While the proposed method segments the methylome in a largely comparable manner to conventional approaches, it has only a single parameter to be tuned. Furthermore, it fully exploits the base-resolution of the data to enable simultaneous detection of methylation changes in even contrasting size ranges, such as focal hypermethylation and global hypomethylation in cancer methylomes. We also propose a simple plot termed methylated domain landscape (MDL) that globally displays the size, the methylation level and the number of the domains thus defined, thereby enabling one to intuitively grasp trends in each methylome. Since the pattern of MDL often reflects cell lineages and is largely unaffected by data size, it can serve as a novel signature of methylome. CONCLUSIONS: Changepoint detection in base-resolution methylome data followed by MDL plotting provides a novel method for methylome characterization and will facilitate global comparison among various WGBS data differing in size and even species origin.

Real-time time-division color electroholography using a single GPU and a USB module for synchronizing reference light.

We propose real-time time-division color electroholography using a single graphics processing unit (GPU) and a simple synchronization system of reference light. To facilitate real-time time-division color electroholography, we developed a light emitting diode (LED) controller with a universal serial bus (USB) module and the drive circuit for reference light. A one-chip RGB LED connected to a personal computer via an LED controller was used as the reference light. A single GPU calculates three computer-generated holograms (CGHs) suitable for red, green, and blue colors in each frame of a three-dimensional (3D) movie. After CGH calculation using a single GPU, the CPU can synchronize the CGH display with the color switching of the one-chip RGB LED via the LED controller. Consequently, we succeeded in real-time time-division color electroholography for a 3D object consisting of around 1000 points per color when an NVIDIA GeForce GTX TITAN was used as the GPU. Furthermore, we implemented the proposed method in various GPUs. The experimental results showed that the proposed method was effective for various GPUs.