SETDB1 deletion causes DNA demethylation and upregulation of multiple zinc-finger genes.

BACKGROUND: SETDB1 (SET domain bifurcated-1) is a histone H3-lysine 9 (H3K9)-specific methyltransferase that mediates heterochromatin formation and repression of target genes. Despite the assumed functional link between DNA methylation and SETDB1-mediated H3K9 trimethylations, several studies have shown that SETDB1 operates autonomously of DNA methylation in a region- and cell-specific manner. This study analyzes SETDB1-null HAP1 cells through a linked methylome and transcriptome analysis, intending to explore genes controlled by SETDB1-involved DNA methylation. METHODS AND RESULTS: We investigated SETDB1-mediated regulation of DNA methylation and gene transcription in human HAP1 cells using reduced-representation bisulfite sequencing (RRBS) and RNA sequencing. While two-thirds of differentially methylated CpGs (DMCs) in genic regions were hypomethylated in SETDB1-null cells, we detected a plethora of C2H2-type zinc-finger protein genes (C2H2-ZFP, 223 of 749) among the DMC-associated genes. Most C2H2-ZFPs with DMCs in their promoters were found hypomethylated in SETDB1-KO cells, while other non-ZFP genes with promoter DMCs were not. These C2H2-ZFPs with DMCs in their promoters were significantly upregulated in SETDB1-KO cells. Similarly, C2H2-ZFP genes were upregulated in SETDB1-null 293T cells, suggesting that SETDB1's function in ZFP gene repression is widespread. There are several C2H2-ZFP gene clusters on chromosome 19, which were selectively hypomethylated in SETDB1-KO cells. CONCLUSIONS: SETDB1 collectively and specifically represses a substantial fraction of the C2H2-ZFP gene family. Through the en-bloc silencing of a set of ZFP genes, SETDB1 may help establish a panel of ZFP proteins that are expressed cell-type specifically and thereby can serve as signature proteins for cellular identity.

Coordinated regulation of microRNA genes in C19MC by SETDB1.

The microRNA (miRNA) gene cluster on chromosome 19, C19MC, is the largest primate-specific miRNA gene cluster. The 46 homologous miRNA genes in C19MC are highly expressed in the placenta, but repressed in other tissues by DNA methylation. Here, we found that the SET domain bifurcated 1(SETDB1), a histone H3-lysine 9 (H3K9)-specific methyltransferase 1, transcriptionally controls C19MC miRNA genes in a coordinated manner in human HAP1 cells. SETDB1 knockout (KO) resulted in the overexpression of C19MC miRNA genes, which was accompanied by a reduction of H3K9 trimethylation (H3K9me3) in the cluster. We found that SETDB1 specifically binds to and modifies the upstream promoter locus of C19MC with H3K9me3, suggesting its role as a C19MC repressor. Overexpression of C19MC miRNA genes was not related to DNA methylation because cytosine methylation levels were not altered in the C19MC of SETDB1 KO cells, indicating that SETDB1 KO does not cause DNA demethylation in the C19MC promoter and body regions. In conclusion, our results suggest that SETDB1 binding and H3K9 methylation at the C19MC promoter and body regions are responsible for the coordinated regulation of miRNA genes in the cluster.

Dnmt1 binds and represses genomic retroelements via DNA methylation in mouse early embryos.

Genome-wide passive DNA demethylation in cleavage-stage mouse embryos is related to the cytoplasmic localization of the maintenance methyltransferase DNMT1. However, recent studies provided evidences of the nuclear localization of DNMT1 and its contribution to the maintenance of methylation levels of imprinted regions and other genomic loci in early embryos. Using the DNA adenine methylase identification method, we identified Dnmt1-binding regions in four- and eight-cell embryos. The unbiased distribution of Dnmt1 peaks in the genic regions (promoters and CpG islands) as well as the absence of a correlation between the Dnmt1 peaks and the expression levels of the peak-associated genes refutes the active participation of Dnmt1 in the transcriptional regulation of genes in the early developmental period. Instead, Dnmt1 was found to associate with genomic retroelements in a greatly biased fashion, particularly with the LINE1 (long interspersed nuclear elements) and ERVK (endogenous retrovirus type K) sequences. Transcriptomic analysis revealed that the transcripts of the Dnmt1-enriched retroelements were overrepresented in Dnmt1 knockdown embryos. Finally, methyl-CpG-binding domain sequencing proved that the Dnmt1-enriched retroelements, which were densely methylated in wild-type embryos, became demethylated in the Dnmt1-depleted embryos. Our results indicate that Dnmt1 is involved in the repression of retroelements through DNA methylation in early mouse development.

Different phases of aging in mouse old skeletal muscle.

With a graying population and increasing longevity, it is essential to identify life transition in later years and discern heterogeneity among older people. Subclassifying the elderly population to inspect the subdivisions for pathophysiological differences is particularly important for the investigation of age-related illnesses. For this purpose, using 24- and 28-month-old mice to represent the "young-old" and "old-old", respectively, we compared their skeletal muscle transcriptomes and found each in a distinct stage: early/gradual (E-aging) and late/accelerated aging phase (L-aging). Principal component analysis showed that the old-old transcriptomes were largely disengaged from the forward transcriptomic trajectory generated in the younger-aged group, indicating a substantial change in gene expression profiles during L-aging. By calculating the transcriptomic distance, it was found that the 28-month group was closer to the two-month group than to the 24-month group. The divergence rate per month for the transcriptomes was the highest in L-aging, twice as fast as the rate in E-aging. Indeed, many of the L-aging genes were significantly altered in transcription, although the changes did not seem random but rather coordinated in a variety of functional gene sets. Of 2,707 genes transcriptionally altered during E-aging, two-thirds were also significantly changed during L-aging, to either downturning or upturning way. The downturn genes were related to mitochondrial function and translational gene sets, while the upturn genes were linked to inflammation-associated gene sets. Our results provide a transcriptomic muscle signature that distinguishes old-old mice from young-old mice. This can help to methodically examine muscle disorders in the elderly.

Intermediate cells of in vitro cellular reprogramming and in vivo tissue regeneration require desmoplakin.

Amphibians and fish show considerable regeneration potential via dedifferentiation of somatic cells into blastemal cells. In terms of dedifferentiation, in vitro cellular reprogramming has been proposed to share common processes with in vivo tissue regeneration, although the details are elusive. Here, we identified the cytoskeletal linker protein desmoplakin (Dsp) as a common factor mediating both reprogramming and regeneration. Our analysis revealed that Dsp expression is elevated in distinct intermediate cells during in vitro reprogramming. Knockdown of Dsp impedes in vitro reprogramming into induced pluripotent stem cells and induced neural stem/progenitor cells as well as in vivo regeneration of zebrafish fins. Notably, reduced Dsp expression impairs formation of the intermediate cells during cellular reprogramming and tissue regeneration. These findings suggest that there is a Dsp-mediated evolutionary link between cellular reprogramming in mammals and tissue regeneration in lower vertebrates and that the intermediate cells may provide alternative approaches for mammalian regenerative therapy.

SETDB1 Overexpression Sets an Intertumoral Transcriptomic Divergence in Non-small Cell Lung Carcinoma.

An increasing volume of evidence suggests that SETDB1 plays a role in the tumorigenesis of various cancers, classifying SETDB1 as an oncoprotein. However, owing to its numerous protein partners and their global-scale effects, the molecular mechanism underlying SETDB1-involved oncogenesis remains ambiguous. In this study, using public transcriptome data of lung adenocarcinoma (ADC) and squamous-cell carcinoma (SCC), we compared tumors with high-level SETDB1 (SH) and those with low-level SETDB1 (comparable with normal samples; SL). The results of principal component analysis revealed a transcriptomic distinction and divergence between the SH and SL samples in both ADCs and SCCs. The results of gene set enrichment analysis indicated that genes involved in the "epithelial-mesenchymal transition," "innate immune response," and "autoimmunity" collections were significantly depleted in SH tumors, whereas those involved in "RNA interference" collections were enriched. Chromatin-modifying genes were highly expressed in SH tumors, and the variance in their expression was incomparably high in SCC-SH, which suggested greater heterogeneity within SCC tumors. DNA methyltransferase genes were also overrepresented in SH samples, and most differentially methylated CpGs (SH/SL) were undermethylated in a highly biased manner in ADCs. We identified interesting molecular signatures associated with the possible roles of SETDB1 in lung cancer. We expect these SETDB1-associated molecular signatures to facilitate the development of biologically relevant targeted therapies for particular types of lung cancer.

Demethylation and derepression of genomic retroelements in the skeletal muscles of aged mice.

Changes in DNA methylation influence the aging process and contribute to aging phenotypes, but few studies have been conducted on DNA methylation changes in conjunction with skeletal muscle aging. We explored the DNA methylation changes in a variety of retroelement families throughout aging (at 2, 20, and 28 months of age) in murine skeletal muscles by methyl-binding domain sequencing (MBD-seq). The two following contrasting patterns were observed among the members of each repeat family in superaged mice: (a) hypermethylation in weakly methylated retroelement copies and (b) hypomethylation in copies with relatively stronger methylation levels, representing a pattern of "regression toward the mean" within a single retroelement family. Interestingly, these patterns depended on the sizes of the copies. While the majority of the elements showed a slight increase in methylation, the larger copies (>5 kb) displayed evident demethylation. All these changes were not observed in T cells. RNA sequencing revealed a global derepression of retroelements during the late phase of aging (between 20 and 28 months of age), which temporally coincided with retroelement demethylation. Following this methylation drift trend of "regression toward the mean," aging tended to progressively lose the preexisting methylation differences and local patterns in the genomic regions that had been elaborately established during the early period of development.

Retroelement Insertion in a CRISPR/Cas9 Editing Site in the Early Embryo Intensifies Genetic Mosaicism.

Continued CRISPR/Cas9-mediated editing activity that allows differential and asynchronous modification of alleles in successive cell generations expands allelic complexity. To understand the earliest events during CRISPR/Cas9 editing and the allelic selection among the progeny of subsequent cell divisions, we inspected in detail the genotypes of 4- and 8-cell embryos and embryonic stem cells (ESCs) after microinjection of a CRISPR toolkit into the zygotes. We found a higher editing frequency in 8-cell embryos than in 4-cell embryos, indicating that the CRISPR/Cas9 activity persisted through the 8-cell stage. Analysis of a CRISPR/Cas9 transgenic founder mouse revealed that four different alleles were present in its organs in different combinations and that its germline included three different mutant alleles, as shown by the genotypes of the pups. The indel depth, which measured the extent of indels at the sequence level within single embryos, decreased significantly as the embryos advanced to form ESCs, suggesting that exclusion of fatal indels occurred in the subsequent cell generations. Interestingly, we discovered that the CRISPR sites frequently contained introduced retroelement sequences and that this occurred preferentially with certain classes of retroelements. Therefore, in addition to CRISPR/Cas9's innate mechanism of separate, differential enzymatic modifications of alleles, the frequent retroelement insertions that occur in early mouse embryos during CRISPR/Cas9 editing further expand the allelic diversity and mosaicism in the resulting transgenic founders.

Determination of Oocyte-Manipulation, Zygote-Manipulation, and Genome-Reprogramming Effects on the Transcriptomes of Bovine Blastocysts.

Somatic cell nuclear transfer (scNT) embryos suffer from damage caused by micro-operation (manipulation) and inefficient genome reprograming that hinder their normal development at different levels and in distinct ways. These two effects are inseparable in the nature of the scNT embryo, although methods to separately measure them are needed to improve scNT technology and evaluate incoming reprogramming tools. As an attempt to meet these demands, we made bovine sham nuclear-transfer (shNT) blastocysts, special embryos made with a standard nuclear-transfer procedure at the zygote stage, while retaining an intact genome. We compared their transcriptomes with those of other blastocysts derived by in-vitro fertilization (IVF) or scNT. Correlation analysis revealed a singularity of shNT blastocysts as they separately gathered from the others. Analysis of developmentally important genes revealed that, in shNTs, the stemness-associated differentially expressed genes (DEGs), including OCT4, were mostly underrepresented. Overrepresented epi-driver genes were largely associated with heterochromatin establishment and maintenance. By multilateral comparisons of their transcriptomes, we classified DEGs into three groups: 561 manipulation-associated DEGs (MADs) common to shNTs and scNTs, 764 donor genome-associated DEGs (DADs) specific to scNTs, and 1743 zygote manipulation-associated DEGs (zMADs) specific to shNTs. GO enrichment analysis generated various terms involving "cell-cell adhesion," "translation," and "transcription" for MADs and "cell differentiation" and "embryo implantation" for DADs. Because of the transcriptomic specificity of shNTs, we studied zMADs in detail. GO enrichment analysis with the 854 zMADs underrepresented in shNTs yielded terms related to protein and mitochondria homeostasis, while GO enrichment analysis of 889 shNT-high zMADs yielded terms related to endoplasmic reticulum stress and protein transport. We summarized the DEGs, which, with further investigation, may help improve our understanding of molecular events occurring in cloned embryos and our ability to control clonal reprogramming.

Age-associated bimodal transcriptional drift reduces intergenic disparities in transcription.

This study addressed the question of how well the quantitative transcriptome structure established in early life is maintained and how consistently it appears with increasing age, and if there is age-associated alteration of gene expression (A3GE), how much influence the Huntington's disease (HD) genotype exerts on it. We examined 285 exonic sequences of 175 genes using targeted PCR sequencing in skeletal muscle, brain, and splenic CD4+ T cells of wild-type and HD mice. In contrast to the muscle and brain, T cells exhibited large A3GE, suggesting a strong contribution to functional decline of the organism. This A3GE was markedly intensified in age-matched HD T cells, which exhibited accelerated aging as determined by reduced telomere length. Regression analysis suggested that gene expression levels change at a rate of approximately 3% per month with age. We found a bimodal relationship in A3GE in T cells in that weakly expressed genes in young mice were increasingly transcribed in older animals whereas highly expressed genes in the young were decreasingly expressed with age. This bimodal transcriptional drift in the T cell transcriptome data causes the differences in transcription rate between genes to progressively reduce with age.

Simultaneous Methylation-Level Assessment of Hundreds of CpG Sites by Targeted Bisulfite PCR Sequencing (TBPseq).

Methylated-DNA sequencing technologies are producing vast amounts of methylome data from cancer samples, from which cancer-associated differentially methylated CpG sites (cDMCs) are continuously identified and filed. The inclusion of as many cDMCs as possible helps improve the accuracy of cancer diagnosis and sometimes identify cancer subtypes. However, the lack of an established method for the analysis of 100s of cDMCs practically impedes their robust use in clinical medicine. Here, we tested the availability of targeted bisulfite-PCR-sequencing (TBPseq) technology for the assessment of methylation levels of a myriad of CpGs scattered over the genome. In randomly selected 46 cancer cell lines, multiplexed PCR yielded a variety of amplicons harboring 246 CpGs residing at promoters of 97 cancer-associated genes, all of which were sequenced in the same flow cell. Clustering analysis of the TBPseq-assessed methylation levels of target CpGs showed that the lung and liver cancer cell lines correlated relatively strongly with each other while they weakly correlated with colon cancer cells. CpGs at the LIFR gene promoter, which are known to be hypermethylated in colon cancers, indeed were heavily methylated in the tested colon cancer cells. Moreover, the LIFR promoter hypermethylation was found in colon cancer cells only, but not in biliary tract, liver, lung, and stomach cancers cell lines. A meta-analysis with public cancer methylome data verified the colon cancer specificity of LIFR promoter methylation. These results demonstrate that our TBPseq-based methylation assessment could be considered an effective, accurate, and competitive method to simultaneously examine a large number of target cDMCs and patient samples.

Characterization of X-Chromosome Gene Expression in Bovine Blastocysts Derived by In vitro Fertilization and Somatic Cell Nuclear Transfer.

To better understand X-chromosome reactivation (XCR) during early development, we analyzed transcriptomic data obtained from bovine male and female blastocysts derived by in-vitro fertilization (IVF) or somatic-cell nuclear transfer (SCNT). We found that X-linked genes were upregulated by almost two-fold in female compared with male IVF blastocysts. The upregulation of X-linked genes in female IVFs indicated a transcriptional dimorphism between the sexes, because the mean autosomal gene expression levels were relatively constant, regardless of sex. X-linked genes were expressed equivalently in the inner-cell mass and the trophectoderm parts of female blastocysts, indicating no imprinted inactivation of paternal X in the trophectoderm. All these features of X-linked gene expression observed in IVFs were also detected in SCNT blastocysts, although to a lesser extent. A heatmap of X-linked gene expression revealed that the initial resemblance of X-linked gene expression patterns between male and female donor cells turned sexually divergent in host SCNTs, ultimately resembling the patterns of male and female IVFs. Additionally, we found that sham SCNT blastocysts, which underwent the same nuclear-transfer procedures, but retained their embryonic genome, closely mimicked IVFs for X-linked gene expression, which indicated that the embryo manipulation procedure itself does not interfere with XCR in SCNT blastocysts. Our findings indicated that female SCNTs have less efficient XCR, suggesting that clonal reprogramming of X chromosomes is incomplete and occurs variably among blastocysts, and even among cells in a single blastocyst.

Age-associated chromatin relaxation is enhanced in Huntington's disease mice.

Expansion of polyglutamine stretch in the huntingtin (HTT) protein is a major cause of Huntington's disease (HD). The polyglutamine part in HTT interacts with various proteins implicated in epigenetic regulation of genes, suggesting that mutant HTT may disturb the integrity of the epigenetic system. Here, we used a PCRseq-based method to examine expression profile of 395 exonic segments from 260 "epi-driver" genes in splenic T lymphocytes from aged HD mice. We identified 67 exonic segments differentially expressed between young and aged HD mice, most of them upregulated in the aged. Polycomb-repressive complex (PRC)-regulated genes (PRGs) were markedly upregulated in aged HD mice, consistent with downregulation of PRC genes. Epi-driver gene categories of lysine-methylation, lysine-demethylation, arginine-methylation, and PRG showed differential age-associated changes between HD and control. Analyzing the pattern of change in epi-driver gene expressions hinted at an enhanced shift in HD chromatin to a more accessible state with age, which was experimentally demonstrated by DNase-I-hypersensitivity sequencing showing increased chromatin accessibility in HD cells compared to control. We suggest the global change can potentially relieve chromatin-induced repression of many genes, and the unintended expressions of some detrimental proteins could alter T cell function to a greater degree in aged HD mice.

The HIST1 Locus Escapes Reprogramming in Cloned Bovine Embryos.

Epigenetic reprogramming is necessary in somatic cell nuclear transfer (SCNT) embryos in order to erase the differentiation-associated epigenetic marks of donor cells. However, such epigenetic memories often persist throughout the course of clonal development, thus decreasing cloning efficiency. Here, we explored reprogramming-refractory regions in bovine SCNT blastocyst transcriptomes. We observed that histone genes residing in the 1.5 Mb spanning the cow HIST1 cluster were coordinately downregulated in SCNT blastocysts. In contrast, both the nonhistone genes of this cluster, and histone genes elsewhere remained unaffected. This indicated that the downregulation was specific to HIST1 histone genes. We found that, after trichostatin A treatment, HIST1 histone genes were derepressed, and DNA methylation at their promoters was decreased to the level of in vitro fertilization embryos. Therefore, our results indicate that the reduced expression of HIST1 histone genes is a consequence of poor epigenetic reprogramming in SCNT blastocysts.

Differences in the gene expression profiles of slow- and fast-forming preinduced pluripotent stem cell colonies.

Induced pluripotent stem cells (iPSCs) are generated through a gradual process in which somatic cells undergo a number of stochastic events. In this study, we examined whether two different doxycycline-inducible iPSCs, slow-forming 4F2A-iPSCs and fast-forming NGFP-iPSCs, have equivalent levels of pluripotency. Multiplex reverse-transcriptase PCR generated gene expression profiles (GEPs) of 13 pluripotency genes in single initially formed-iPSC (if-iPSC) colonies of NGFP and 4F2A group. Assessment of GEP difference using a weighted root mean square deviation (wRMSD) indicates that 4F2A if-iPSCs are more closely related to mESCs than NGFP if-iPSCs. Consistently, Nanog and Sox2 genes were more frequently derepressed in 4F2A if-iPSC group. We further examined 20 genes that are implicated in reprogramming. They were, overall, more highly expressed in NGFP if-iPSCs, differing from the pluripotency genes being more expressed in 4F2A if-iPSCs. wRMSD analysis for these reprogramming-related genes confirmed that the 4F2A if-iPSC colonies were less deviated from mESCs than the NGFP if-iPSC colonies. Our findings suggest that more important in attaining a better reprogramming is the mode of action by the given reprogramming factors, rather than the total activity of them exerting to the cells, as the thin-but-long-lasting mode of action in 4F2A if-iPSCs is shown to be more effective than its full-but-short-lasting mode in NGFP if-iPSCs.

Transcriptomic Features of Bovine Blastocysts Derived by Somatic Cell Nuclear Transfer.

Reprogramming incompletely occurs in most somatic cell nuclear transfer (SCNT) embryos, which results in misregulation of developmentally important genes and subsequent embryonic malfunction and lethality. Here we examined transcriptome profiles in single bovine blastocysts derived by in vitro fertilization (IVF) and SCNT. Different types of donor cells, cumulus cell and ear-skin fibroblast, were used to derive cSCNT and fSCNT blastocysts, respectively. SCNT blastocysts expressed 13,606 genes on average, similar to IVF (13,542). Correlation analysis found that both cSCNT and fSCNT blastocyst groups had transcriptomic features distinctive from the IVF group, with the cSCNT transcriptomes closer to the IVF ones than the fSCNT. Gene expression analysis identified 56 underrepresented and 78 overrepresented differentially expressed genes in both SCNT groups. A 400-kb locus harboring zinc-finger protein family genes in chromosome 18 were found coordinately down-regulated in fSCNT blastocysts, showing a feature of reprogramming-resistant regions. Probing into different categories of genes important for blastocyst development revealed that genes involved in trophectoderm development frequently were underrepresented, and those encoding epigenetic modifiers tended to be overrepresented in SCNT blastocysts. Our effort to identify reprogramming-resistant, differentially expressed genes can help map reprogramming error-prone loci onto the genome and elucidate how to handle the stochastic events of reprogramming to improve cloning efficiency.

Rapid generation of secondary fibroblasts through teratoma formation.

The use of secondary or reprogrammable cells in the production of induced pluripotent stem cells (iPSCs) circumvents random infection by various viral particles and random, uncontrollable integrations of the viral genomes into different genomic loci. We have developed a convenient method for repeatedly producing genetically identical secondary fibroblasts via teratoma formation using pre-existing iPSCs. The iPSCs used in this study carried doxycycline (Dox)-inducible transgenes for four transcription factors in their genome. Teratoma-derived primary cells (TOFs) were obtained in a huge amount during the culture of teratomas and showed good ability to form iPSCs similar to that of regular secondary fibroblasts. Immunohistochemistry analysis demonstrated the potential of TOF-derived iPSCs to differentiate into all three germ layers. The gene expression profiles of these TOFs and their iPSCs closely mimicked those of regular embryonic fibroblasts and embryonic stem cells/iPSCs, respectively. The possibility that the iPSCs were derived from a small part of pluripotent cells lurking in the TOF population was precluded by the observation of doxycycline-dependent and PluriSin (a compound selectively eliminating pluripotent cells)- independent formations of iPSCs. Our results showed that the TOFs retained the capability to mediate cellular reprogramming, similar to that of regular secondary fibroblasts.