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.

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.

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.

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.

Endometrial senescence is mediated by interleukin 17 receptor B signaling.

BACKGROUND: We previously identified Il17RB, a member of the IL17 superfamily, as a candidate marker gene for endometrial aging. While IL17RB has been linked to inflammation and malignancies in several organ systems, its function in the endometrium has not been investigated and is thus poorly understood. In the present study, we performed a functional analysis of this receptor with the aim of determining the effects of its age-associated overexpression on the uterine environment. METHODS: We analyzed IL17RB-related signaling pathways and downstream gene expression in an immortalized human endometrial glandular epithelial cell line ("hEM") forced to express the receptor via lentiviral transduction ("IL17RB-hEM"). We also prepared endometrial organoids from human endometrial tissue sourced from hysterectomy patients ("patient-derived EOs") and exposed them to cytokines that are upregulated by IL17RB expression to investigate changes in organoid-forming capacity and senescence markers. We analyzed RNA-seq data (GEO accession number GSE132886) from our previous study to identify the signaling pathways associated with altered IL17RB expression. We also analyzed the effects of the JNK pathway on organoid-forming capacity. RESULTS: Stimulation with interleukin 17B enhanced the NF-κB pathway in IL17RB-hEM, resulting in significantly elevated expression of the genes encoding the senescence associated secretory phenotype (SASP) factors IL6, IL8, and IL1ß. Of these cytokines, IL1ß inhibited endometrial organoid growth. Bioinformatics analysis showed that the JNK signaling pathway was associated with age-related variation in IL17RB expression. When IL17RB-positive cells were cultured in the presence of IL17B, their organoid-forming capacity was slightly but non-significantly lower than in unexposed IL17RB-positive cells, but when IL17B was paired with a JNK inhibitor (SP600125), it was restored to control levels. Further, IL1ß exposure significantly reduced organoid-forming capacity and increased p21 expression in endometrial organoids relative to non-exposure (control), but when IL1ß was paired with SP600125, both indicators were restored to levels comparable to the control condition. CONCLUSIONS: We have revealed an association between IL17RB, whose expression increases in the endometrial glandular epithelium with advancing age, and cellular senescence. Using human endometrial organoids as in vitro model, we found that IL1ß inhibits cell proliferation and leads to endometrial senescence via the JNK pathway.

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.

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.

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.

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.

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.

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.

Real-time spatiotemporal division multiplexing electroholography with a single graphics processing unit utilizing movie features.

We propose a real-time spatiotemporal division multiplexing electroholography utilizing the features of movies. The proposed method spatially divides a 3-D object into plural parts and periodically selects a divided part in each frame, thereby reconstructing a three-dimensional (3-D) movie of the original object. Computer-generated holograms of the selected part are calculated by a single graphics processing unit and sequentially displayed on a spatial light modulator. Visual continuity enables a reconstructed movie of the original 3-D object. The proposed method realized a real-time reconstructed movie of a 3-D object composed of 11,646 points at over 30 frames per second (fps). We also displayed a reconstructed movie of a 3-D object composed of 44,647 points at about 10 fps.

Gene network inference and visualization tools for biologists: application to new human transcriptome datasets.

Gene regulatory networks inferred from RNA abundance data have generated significant interest, but despite this, gene network approaches are used infrequently and often require input from bioinformaticians. We have assembled a suite of tools for analysing regulatory networks, and we illustrate their use with microarray datasets generated in human endothelial cells. We infer a range of regulatory networks, and based on this analysis discuss the strengths and limitations of network inference from RNA abundance data. We welcome contact from researchers interested in using our inference and visualization tools to answer biological questions.

Vasohibin-1 is identified as a master-regulator of endothelial cell apoptosis using gene network analysis.

BACKGROUND: Apoptosis is a critical process in endothelial cell (EC) biology and pathology, which has been extensively studied at protein level. Numerous gene expression studies of EC apoptosis have also been performed, however few attempts have been made to use gene expression data to identify the molecular relationships and master regulators that underlie EC apoptosis. Therefore, we sought to understand these relationships by generating a Bayesian gene regulatory network (GRN) model. RESULTS: ECs were induced to undergo apoptosis using serum withdrawal and followed over a time course in triplicate, using microarrays. When generating the GRN, this EC time course data was supplemented by a library of microarray data from EC treated with siRNAs targeting over 350 signalling molecules.The GRN model proposed Vasohibin-1 (VASH1) as one of the candidate master-regulators of EC apoptosis with numerous downstream mRNAs. To evaluate the role played by VASH1 in EC, we used siRNA to reduce the expression of VASH1. Of 10 mRNAs downstream of VASH1 in the GRN that were examined, 7 were significantly up- or down-regulated in the direction predicted by the GRN.Further supporting an important biological role of VASH1 in EC, targeted reduction of VASH1 mRNA abundance conferred resistance to serum withdrawal-induced EC death. CONCLUSION: We have utilised Bayesian GRN modelling to identify a novel candidate master regulator of EC apoptosis. This study demonstrates how GRN technology can complement traditional methods to hypothesise the regulatory relationships that underlie important biological processes.

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.