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.

Skin microbe-dependent TSLP-ILC2 priming axis in early life is co-opted in allergic inflammation.

Although early life colonization of commensal microbes contributes to long-lasting immune imprinting in host tissues, little is known regarding the pathophysiological consequences of postnatal microbial tuning of cutaneous immunity. Here, we show that postnatal exposure to specific skin commensal Staphylococcus lentus (S. lentus) promotes the extent of atopic dermatitis (AD)-like inflammation in adults through priming of group 2 innate lymphoid cells (ILC2s). Early postnatal skin is dynamically populated by discrete subset of primed ILC2s driven by microbiota-dependent induction of thymic stromal lymphopoietin (TSLP) in keratinocytes. Specifically, the indole-3-aldehyde-producing tryptophan metabolic pathway, shared across Staphylococcus species, is involved in TSLP-mediated ILC2 priming. Furthermore, we demonstrate a critical contribution of the early postnatal S. lentus-TSLP-ILC2 priming axis in facilitating AD-like inflammation that is not replicated by later microbial exposure. Thus, our findings highlight the fundamental role of time-dependent neonatal microbial-skin crosstalk in shaping the threshold of innate type 2 immunity co-opted in adulthood.

Parallel shift of DNA methylation and gene expression toward the mean in mouse spleen with aging.

Age-associated DNA-methylation drift (AMD) manifests itself in two ways in mammals: global decrease (hypomethylation) and local increase of DNA methylation (hypermethylation). To comprehend the principle behind this bidirectional AMD, we studied methylation states of spatially clustered CpG dinucleotides in mouse splenic DNA using reduced-representation-bisulfite-sequencing (RRBS). The mean methylation levels of whole CpGs declined with age. Promoter-resident CpGs, generally weakly methylated (<5%) in young mice, became hypermethylated in old mice, whereas CpGs in gene-body and intergenic regions, initially moderately (~33%) and extensively (>80%) methylated, respectively, were hypomethylated in the old. Chromosome-wise analysis of methylation revealed that inter-individual heterogeneities increase with age. The density of nearby CpGs was used to classify individual CpGs, which found hypermethylation in CpG-rich regions and hypomethylation in CpG-poor regions. When genomic regions were grouped by methylation level, high-methylation regions tended to become hypomethylated whereas low-methylation regions tended to become hypermethylated, regardless of genomic structure/function. Data analysis revealed that while methylation level and CpG density were interdependent, methylation level was a better predictor of the AMD pattern representing a shift toward the mean. Further analysis of gene-expression data showed a decrease in the expression of highly-expressed genes and an increase in the expression of lowly-expressed genes with age. This shift towards the mean in gene-expression changes was correlated with that of methylation changes, indicating a potential link between the two age-associated changes. Our findings suggest that age-associated hyper- and hypomethylation events are stochastic and attributed to malfunctioning intrinsic mechanisms for methylation maintenance in low- and high-methylation regions, respectively.

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.

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.

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.

Functional dissection of N-terminal nuclear trafficking signals of SETDB1.

SETDB1 is a histone H3-lysine 9-specific methyltransferase that fulfills epigenetic functions inside the nucleus; however, when overexpressed, SETDB1 majorily localizes in the cytoplasm. SETDB1 has a single nuclear-localization-signal (NLS) motif and two successive nuclear-export-signal (NES1 and NES2) motifs in the N-terminus, suggesting that SETDB1 localization is the consequence of a balance between the two antithetic motifs. Here, we performed a series of motif deletions to characterize their effects on the cellular movement of SETDB1. Given the cytoplasmic localization of GFP-SETDB1 in the whole form, without the NES motifs, GFP-SETDB1 was not nuclear, and 3xNLS addition plus NES removal held the majority of GFP-SETDB1 within the nucleus. The results indicated that the cytoplasmic localization of GFP-SETDB1 is the combined result of weak NLS and robust NESs. In ATF7IP-overexpressing cells, GFP-SETDB1 entered the nucleus only in the presence of the NES1 motif; neither the NES2 nor NLS motif was necessary. Since subcellular fractionation results showed that ATF7IP was nuclear-only, an intermediary protein may interact specifically with the NES1 motif after stimulation by ATF7IP. When GFP-SETDB1 had either NES1 or NES2, it was precipitated (in immunoprecipitation) and colocalized (in immunofluorescence) with ATF7IP, indicating that GFP-SETDB1 interacts with ATF7IP through the NES motifs in the nucleus. The regulated nuclear entry of SETDB1 is assumed to set a tight restriction on its abundance within the nucleus, thereby ensuring balanced nuclear SETDB1 levels.

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.

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.

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.

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.

K120R mutation inactivates p53 by creating an aberrant splice site leading to nonsense-mediated mRNA decay.

The point mutation that substitutes lysine with arginine at position 120 of human p53 has been characterized as a missense mutation. The K120R mutation renders the p53 protein disabled for acetylation and, as a result, defective for apoptotic function, which provides a mechanistic link between the missense mutation and tumorigenesis. However, we noticed the failures of tumorigenesis in mice with the mutation, and of the related studies to notice that it has arbitrarily reflected in amino acid change through a sequence modification (AGA) of the original tumor mutation (AGG) by codon degeneracy. Unlike this modified version, we also discovered a novel splicing site the original mutation, TP53 c.359A>G, may induce. Using a human induced pluripotent stem cell line that was engineered to be homozygous for the original mutation, we here identified that the accidental splicing site generates a defective transcript variant with a frame-shifted premature termination codon which is subjected to nonsense-mediated mRNA decay. The authentic splicing still occurs but in extremely low amounts. Taken together, this mutation causes depletion of cellular p53 via defective mRNA, suggesting a new link to tumorigenesis.

Surveillance of Retroelement Expression and Nucleic-Acid Immunity by Histone Methyltransferase SETDB1.

In human cancers, histone methyltransferase SETDB1 (SET domain, bifurcated 1) is frequently overexpressed but its significance in carcinogenesis remains elusive. A recent study shows that SETDB1 downregulation induces de-repression of retroelements and innate immunity in cancer cells. The possibility of SETDB1 functioning as a surveillant of retroelement expression is discussed in this study: the cytoplasmic presence of retroelement-derived nucleic acids (RdNAs) drives SETDB1 into the nucleus by the RNA-interference route, rendering the corresponding retroelement transcriptionally inert. These RdNAs could, therefore, be signals of genome instability sent out for SETDB1 present in the cytoplasm to maintain genome integrity.

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.

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.

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.

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.