Deletion of Tet proteins results in quantitative disparities during ESC differentiation partially attributable to differences in gene expression

Background: The Tet protein family (Tet1, Tet2, and Tet3) regulate DNA methylation through conversion of 5-methylcytosine to 5-hydroxymethylcytosine which can ultimately result in DNA demethylation and play a critical role during early mammalian development and pluripotency¬. While multiple groups have generated knockouts combining loss of different Tet proteins in murine embryonic stem cells (ESCs), differences in genetic background and approaches has made it dicult to directly compare results and discern the direct mechanism by which Tet proteins regulate the transcriptome. To address this concern, we utilized genomic editing in an isogenic pluripotent background which permitted a quantitative, ow-cytometry based measurement of pluripotency in combination with genome-wide assessment of gene expression and DNA methylation changes. Our ultimate goal was to generate a resource of large-scale datasets to permit hypothesis-generating experiments. Results: We demonstrate a quantitative disparity in the differentiation ability among Tet protein deletions, with Tet2 single knockout exhibiting the most severe defect, while loss of Tet1 ¬alone or combinations of Tet genes showed a quantitatively intermediate phenotype. Using a combination of transcriptomic and epigenomic approaches we demonstrate an increase in DNA hypermethylation and a divergence of transcriptional proles in pluripotency among Tet deletions, with loss of Tet2 having the most profound effect in undifferentiated ESCs. Conclusions: We conclude that loss of Tet2 has the most dramatic effect both on the phenotype of ESCs and the transcriptome compared to other genotypes. While loss of Tet proteins increased DNA hypermethylation, especially in gene promoters, these changes in DNA methylation did not correlate with gene expression changes. Thus, while loss of different Tet proteins alters DNA methylation, this change does not appear to be directly responsible for transcriptome changes. Thus, loss of Tet proteins likely regulates the transcriptome epigenetically both through altering 5mC but also through additional mechanisms. Nonetheless, the transcriptome changes in pluripotent Tet2-/- ESCs compared to wild-type implies that the disparities in differentiation can be

Background DNA methylation plays a critical role in regulating gene expression during development and is maintained in mammals through a complex interplay between DNA methyltransferases (DNMT) and a family of proteins termed Ten Eleven Translocation (Tet). DNMTs are the "writer" of DNA methylation, whereas Tet proteins are the "erasers" [1][2][3][4]. All Tet members are capable of oxidizing methylated cytosines (5mC) to 5hydroxymethylated cytosine (5hmC). Further oxidation of the 5hmC mark by Tet proteins yields 5-carboxy cytosine (5caC) and 5-formyl cytosine (5fC; [5]) which are unstable and removed by base excision repair (BER; [6][7][8]), resulting in conversion to an unmodi ed cytosine and is referred to as active DNA demethylation. Tet proteins can also cause DNA demethylation through a passive mechanism. The 5hmC mark is not recognized by maintenance DNMTs during DNA replication and will be absent in the daughter strand, resulting in passive DNA demethylation in a cell-cycle dependent manner. Regulation of DNA methylation is especially important during early embryogenesis where dynamic, rapid changes in the epigenome are required for early differentiation steps.
Tet1 and Tet2 are highly expressed in the inner cell mass (ICM; [9]) from which murine embryonic stem cells (ESCs) are derived, but studies in ESCs on the role of Tet proteins have resulted in con icting reports. In vivo studies of Tet deletion have established that loss of either Tet1 [10] or Tet2 [11] alone had no effect on early development. The combination deletion of Tet1 and Tet2 (double knockout, DKO) exhibited a partially penetrant phenotype in which a fraction of embryos died perinatally, but others were overtly normal [12]. Even though Tet3 is only expressed at low levels in ESCs and early embryos, combined loss of Tet1, Tet2, and Tet3 in ESCs prevents them from contributing to chimeric embryos during blastocyst complementation [13], implying that some amount of Tet protein activity is required for normal embryonic development. Importantly, because of the complex cell:cell interactions within the developing embryo, as well as differences in mouse strain background, it has not been possible to translate these in vivo results into how Tet proteins regulate ESCs and pluripotency. There remains debate about which Tet genes are required for pluripotency and/or differentiation, and how they may interact in combination remains unclear. Part of this likely relates to the differences in how the ESC were isolated, differences between ESC genetic background, as well as the use of classical homologous recombination approaches which induce large genomic deletions [10,12,13]. Genomic editing, which can generate null alleles by inducing small frame-shifts to prevent functional protein production rather than large genomic deletions, has been used at all the Tet genes [14,15], but the resulting ESCs have not been carefully quantitatively characterized in terms of their differentiation and pluripotency. To address this, we utilized a GFP reporter line to permit quantitative tracking of pluripotency in combination with combinatorial genomic-editing to ablate all three Tet genes in an isogenic background. By utilizing both transcriptome and genome-wide DNA methylation analysis we establish that the differences we observe are at least partially related to baseline differences in the epigenome and transcriptome of these cells.

Deletion of Tet proteins Does Not Disrupt Pluripotency
We utilized a previously generated ESC line in which eGFP was "knocked-in" the Pou5f1 locus, permitting use of eGFP expression as a quantitative marker of pluripotency ( Figures 1A-C). A single clone was used as a parental line for all further experiments and hereafter referred to as wild-type (WT). Published gRNAs [14] to all three Tet genes were transfected into the parental line to generate the following lines: Tet1-/-, Tet2-/-, DKO (Tet1-/-:Tet2-/-), and TKO (Tet1-/-:Tet2-/-:Tet3-/-; Figure 1A). Unless otherwise noted, all experiments were performed with 3 individual clones of WT, Tet1-/and 2 clones of Tet2-/-, DKO, TKO. Generation of indels by genomic editing was con rmed by direct sequencing (data not shown) and loss of protein by Western blot ( Figure 1D). All Tet deletion lines expressed similar levels of eGFP ( Figure 1E) and common pluripotency markers ( Figure 1D). Expression of three lineage markers, Gata6 (endoderm), Cdx2 (trophectoderm), and Brachyury (mesoderm) were variable across lines ( Figure 1F) when measured by RT-qPCR, likely due to clone-to-clone variation given the very low-level expression of these markers.
These data are consistent with published reports that deletion of Tet proteins in any combination does not alter pluripotency per se.

Deletion of Tet Proteins Impairs Differentiation
To determine the effect of Tet deletion on differentiation we performed a leukemia inhibitory factor (LIF) withdrawal assay. LIF maintains ESC pluripotency and its withdrawal promotes differentiation down all three germ layers. All lines were grown without LIF for six days and differentiation assayed by ow cytometry for eGFP expression. All Tet deletion lines, including independent clones for each genotype, were resistant to differentiation as measured by retention of eGFP expression, consistent with prior publications (Figure 2A, B; [10,12,13,16]) . Surprisingly, Tet2-/clones while quantitatively more resistant to differentiation as expected [17] than the other lines, the DKO and TKO lines more closely resembled the Tet1-/-. This suggests that Tet1 deletion is the "dominant" phenotype since the DKO and TKO differentiation mirrors Tet1-/and not Tet2-/-. Given the agreement among clones, we used multiple clones for all subsequent analyses and grouped them as biological replicates, thereby allowing interclonal variations to drive our statistical tests.
To con rm the loss of pluripotency and resistance to differentiation using an independent approach, we performed Western blots for the pluripotency factors Nanog and Oct3/4 and the mesoderm marker Brachyury/T ( Figure 2C). All Tet genotypes retained higher levels of both Nanog and Oct3/4 when compared to the WT ( Figure 2C). Importantly, Oct3/4 protein levels as measured by Western blot closely resemble the eGFP expression, indicating our ow cytometry approach is a quantitative measurement of pluripotency. To identify morphological differences during differentiation, cells were stained for alkaline phosphatase (AP) at day 0 and 6 of LIF withdrawal ( Figure 2E). At D0 no gross morphological changes were observed in any Tet deletion lines ( Figure 2E, top panel). Following six days of LIF withdrawal, WT colonies at D6 appear atter and are negative for AP, Tet2-/colonies were less spread out than WT and remained AP positive whereas the other lines displayed an intermediate phenotype ( Figure 2E, bottom panel). In agreement with several published studies, we found no defects in pluripotency following deletion of Tet proteins but uncovered a quantitative discrepancy between Tet2-/and other genotypes which phenocopied Tet1-/-. This is in contrast to previous studies that suggest TKO cells display the most severe differentiation block [10,12,13,16]. Overall, our use of an isogenic parental line and multiple clones demonstrated a unique quantitative disparity among Tet proteins during ESC differentiation.

Transcriptional Pro ling of Tet Deletions in Pluripotency and LIF Withdrawal
To delineate the disparities in differentiation among the various genotypes, RNA-seq was performed at D6 of LIF withdrawal. We chose LIF withdrawal because it permits ESCs to differentiate down multiple lineages through loss of a signaling cascade rather than a strong positive differentiation signal which may promote its own set of transcriptome changes. We performed RNA sequencing (RNA-seq) on differentiated (D6) samples and read counts were normalized using ERCC spike-ins to enhance quanti cation [18]. D6 analysis of the top 1000 most variable genes across genotypes revealed a unique transcriptional pro le of Tet2-/cells that was not shared by other Tet deletions or the WT based upon unsupervised hierarchical clustering, consistent with the Tet2-/cells exhibiting a very different cell state after LIF withdrawal compared with the other genotypes ( Figure 3A). Consistent with our eGFP analysis, the Tet1-/-, DKO, and TKO cells most closely resembled each other at the transcriptome level, again indicating that the Tet1-/intermediate differentiation phenotype was dominant over Tet2-/-. To further re ne which genes were differentially expressed, we utilized both a statistical (adjusted p-val <0.05) and fold change (>2-fold) criteria to generate re ned gene lists which likely represented the most differentially expressed as compared to WT cells differentiated in parallel. In terms of the genes which were underexpressed ( Figure 3B) or overexpressed ( Figure 3C) we found that the Tet2-/cells exhibited the most uniquely dysregulated genes compared to WT, consistent with Tet2-/being in a distinct cell state compared to differentiated WT cells. Our results are consistent with the Tet2-/cells at D6 of differentiation being very different from WT, with Tet1-/-, DKO, and TKO cells representing an intermediate defect in differentiation following LIF withdrawal.
Given our results, we wondered whether the loss of Tet genes induced baseline transcriptome changes in the ESCs which were compatible with pluripotency, or alternatively Tet proteins were only required during differentiation as posited in [17]. We performed RNA-seq on undifferentiated ESCs (D0) using the same approach as above. The top 1000 most variable genes are displayed as a heat map in Figure 4A and indicate that Tet2-/cells are transcriptionally distinct from the other genotypes at D0. In contrast to the D6 results, the TKO cells appear to be distinct at D0 than the Tet1-/-, WT, and DKO cells. This was surprising, since Tet3 is generally expressed at low levels in ESCs, although other groups have shown that Tet3 expression is upregulated in DKO cells with minimal effects on global levels of 5hmC [12]. Nonetheless, our results are consistent with the Tet2-/having a very different gene expression pattern compared to the other genotypes. Given the role of Tet proteins in DNA demethylation, it was not surprising that compared to WT cells at D0 our differential expression analysis showed virtually all genes were downregulated ( Figure 4A). Not surprisingly, genes expressed at lower levels in Tet2-/versus WT cells were predominantly unique to this genotype and not shared with the others, again indicating that loss of Tet2 induced a unique change on the transcriptome in ESCs compared to the other genotypes at D0.
Prior work from others has demonstrated that loss of Tet2 in ESCs delays differentiation because of a failure to activate enhancers critical to differentiation programs [17]. Thus, loss of Tet proteins could prevent proper differentiation either through baseline transcriptome differences, an inability to activate differentiation-critical gene expression programs, or perhaps both. To address this, gene lists for signi cantly downregulated genes at D0 were compared for overlap between different knockouts ( Figure  4B). Tet2-/had the most downregulated genes (299) as compared to the WT and showed little overlap with any other Tet genotype besides TKO (60 shared). Taken together, we conclude that Tet proteins predominantly function to activate gene expression and Tet2-/induces unique transcriptome changes not shared by other Tet deletions at D0. This would be consistent with baseline transcriptome differences in the different genotypes at D0. Collectively, our data along with the literature [17] would suggest loss of Tet proteins prevents differentiation through multiple mechanisms.
To obtain a more global view of transcriptome differences between the various genotypes at both D0 and D6, we performed a principal component analysis (PCA; Figure 4C). Consistent with our prior analysis, differentiated (D6) WT cells segregated into a unique quadrant, indicating they were highly distinct from the other cells. D6 DKO, TKO, and Tet1-/clustered together and were distinct from the D6 Tet2-/-, again con rming, using an alternative analytic approach, that D6 Tet2-/cells were different than other genotypes. Surprisingly, the D0 cells, all of which display a pluripotent phenotype, were less clustered than anticipated. The D0 WT cells were the most distinct, and again the D0 DKO, TKO, and Tet1-/formed a cluster. The D0 Tet2-/cells were distinct from the other cells again, indicating that at baseline their transcriptome is distinct from the other genotypes.
To better delineate the pathways which may be altered in Tet1-/and Tet2-/as compared to WT cells we rst performed Gene Ontology analysis using PANTHER [19] to assign differentially expressed downregulated transcripts to their respective biological process ( Figure 5A). Because Tet1-/had far fewer down-regulated transcripts (24) then Tet2-/-(421), fewer overall pathways were identi ed overall. Importantly, there was rough agreement in general between the processes effected by loss of either Tet1 or Tet2. Thus, while loss of Tet2 had a more dramatic effect on the overall transcriptome then loss of Tet1, there was general agreement between the biological processes impacted by loss of either Tet1 or Tet2.
Tet Deletion Alters DNA Methylation which does not correlate with gene expression changes.
We next hypothesized that changes in DNA methylation in pluripotent ESCs explains their transcriptome differences, at least partially. To quantitate DNA methylation differences, we performed reduced representation bisul te sequencing (RRBS) on D0 samples. RRBS is highly quantitative and provides robust coverage of CpG islands which are abundant in mammalian gene promoters. Differentially methylated regions (DMRs) were de ned as having a change of at least 25% compared to WT. Overall, we observed more hypermethylated than hypomethylated regions in all Tet genotypes, consistent with their role in DNA demethylation ( Figure 5B; [3,4]). It should be noted that the Tet2-/cells displayed the most hypermethylation overall compared to wild-type, whereas the Tet1-/-, DKO, and TKO all showed a similar phenotype. DMRs were found at consistent ratios in all genotypes across promoters, introns, exons, and intergenic elements (data not shown), consistent with a genome-wide change in DNA methylation. Importantly, when we looked at a small window around TSS (+/-2kb) corresponding to gene promoters, we observed a larger predominance of DNA hypermethylation rather than hypomethylation ( Figure 5C), consistent with prior literature that Tet proteins bind to a large fraction of promoters [20,21] To determine if changes in DNA methylation correlated with changes in mRNA, we identi ed DMRs within promoters for each genotype and then mapped the fold change in the gene ( Figure 6). Surprisingly, changes in DNA methylation at gene promoters did not correlate with transcriptional changes in the D0 RNA-seq. This is consistent with prior literature [10,12,13,20,21]. Thus, even though we utilized genomic editing in an isogenic background with multiple clones, we were unable to mechanistically link changes in DNA methylation to altered gene expression.

Discussion
In this study we used CRISPR/Cas9 to generate Tet1-/-, Tet2-/-, DKO, and TKO in an isogenic ESC line. In accordance with published studies, we did not observe any perturbation of pluripotency following Tet deletion. We did however uncover a quantitative disparity across genotypes during differentiation. All genotypes were resistant to differentiation, consistent with published reports [10,12,13,16]. Importantly, these experiments revealed three key ndings. First, while loss of Tet proteins caused a block in differentiation, the loss of Tet2 induced a near complete block with a >90% retention of eGFP expression. Second, loss of Tet1 caused only a partial reduction in differentiation. Third, is that in both DKO and TKO cells their differentiation potential phenocopied loss of Tet1 rather than loss of Tet2, implying that the Tet1-/phenotype is dominant. This quantitative difference has not previously been described within the literature and was most likely di cult to conclusively observe because of the lack of isogenic backgrounds and our use of at least two, independent clones for all genotypes minimizes the chance that differences are simply related to clone-to-clone variation. In addition, by using LIF withdrawal rather than a potent differentiation agent such as retinoic acid we were able to observe a more subtle but nonetheless quantitative difference among genotypes. Importantly, many labs have published using TKO embryonic cells representing the most dramatic phenotype [13,22]. Our data would argue this may not be the case, since the Tet2-/-ESCs exhibit a substantial, quantitative reduction in their ability to differentiate. Importantly, while our study focused on the role of Tet proteins in vitro, these ndings are inconsistent with the observed in vivo phenotypes. As mentioned earlier, loss of either Tet1 or Tet2 appears to be welltolerated by embryos, and in fact double knockouts (Tet1-/-;Tet2-/-) are able to generate healthy adults, although not at the expected Mendelian ratios [10][11][12]. Only triple knock-out (Tet1-/-;Tet2-/-;Tet3-/-) animals exhibit early developmental defects, illustrating that the in vivo and in vitro results of ESCs are not fully compatible [13]. Given the complex environment in vivo where both external signals including cell-cell interactions and diffusible signaling factors play critical roles in organismal development, it is perhaps not surprising there is discordance between our in vitro results and in vivo studies. Nonetheless, while our studies may not lend insights into the in vivo role of Tet proteins, they do provide some key insights into how they operate to regulate the transcriptome of pluripotent cells.
A second key area of disparity with the literature surrounds seminal studies by Meelad Dawlaty and Rudy Jaenisch [10,12,13], in which ESCs de cient in different combinations of Tet proteins were derived. In these studies, classical homologous recombination was used to generate animals with combinatorial loss of Tets, and then ESCs were derived from blastocysts. While in vivo these cells contributed poorly to chimeras during blastocyst complementation, they were able to form embryoid bodies in vitro, consistent with our ndings that TKO cells have an impaired but not complete block in differentiation. This disparity between our groups results and the Jaenisch ndings could be for two different reasons. The rst may be technical, given the signi cant differences between genomic editing in a single, parental ESC line as opposed to classical homologous recombination, breeding adult animals, followed by derivation of ESCs. The second is our use of an Oct4-EGFP reporter, originally developed by the Jaenisch group [15,23], which permits quantitative measurement of pluripotency and a direct comparison among the different genotypes. Irrespective of the differences, it is clear from our work and others [17] that loss of Tet2 has the most dramatic phenotype on ESC differentiation, but further studies will be required to explain the discordance between our work and the studies from others.
The second part of this study used genome-wide approaches to quantitate transcriptome and epigenomic differences. To our surprise given their phenotypic similarity at D0, RNA-seq showed substantial gene expression differences between Tet2-/and other genotypes. These changes did not affect pluripotency but imply the differences during LIF withdrawal are at least partly attributable to baseline disparities in the transcriptomes between genotypes. Importantly, we cannot exclude that Tet2 is also required for remodeling chromatin during differentiation [17], which warrants further investigation to tease apart these possibilities. From the D0 results we conclude that Tet2-/cells induce a unique transcriptional program as compared to loss of the other Tets. Given the profound block in differentiation we saw in the Tet2-/-ESCs, we were surprised to nd that additional deletion of Tet1-/in DKO or TKO cells more closely phenocopied Tet1-/cells, both in terms of their differentiation potential and their baseline transcriptome. This phenotype requires further investigation, since in terms of baseline transcriptome changes loss of Tet1 alone had modest effects. This implies that the combinatorial loss of both Tet1 and Tet2 that is critical to the intermediate phenotype. Importantly, given the close correlation between the baseline D0 transcriptome of Tet1-/-, DKO, and TKO cells, it is likely that the intermediate differentiation phenotype again is at least partially attributable to the role of different Tet proteins in chromatin remodeling during differentiation. Nonetheless, this pattern of differentiation block among the various genotypes could only have been uncovered using a single, parental, isogenic line which easily permitted the quantitative measurement of pluripotency.
Since Tet proteins canonically function in DNA demethylation, we hypothesized that changes in D0 RNAseq would be secondary to alterations in DNA methylation. Overall, we observed a majority of DMRs in the Tet deletion genotypes were hypermethylated compared to WT. Tet2-/had the largest number of hypermethylated DMRs even though this genotype had the second fewest of total DMRs. Nonetheless, our data is consistent with the work of others, in that there was little correlation between the changes in DNA methylation and transcriptome changes we observed [9,10,12,13,20]. This would imply that the observed transcriptome changes are separate from the role of Tet proteins in active DNA demethylation. It may be that our transcriptome differences are secondary to loss of 5hmC, which we did not directly measure. Alternatively, it may be the result of Tet proteins modulating gene expression independent of their oxidative capacity. In either case, this important question remains unresolved and requires further study.

CONCLUSIONS
Tet proteins are required for proper embryonic development [10,12,13] and in ESC by interacting with the pluripotency-associated transcription factor Nanog [24]. While the in vivo roles during early embryogenesis has been well characterized, there remains con icting literature about the role either individually or combinatorially for the Tet proteins in ESCs. To address this gap, we utilized genomic editing to delete different Tet genes individually and in combination. Key ndings are that loss of Tet2 by itself displayed the most signi cant, quantitative block in differentiation during LIF withdrawal. What is most surprising is that Tet1-/-, DKO, and TKO cells phenocopied each other during differentiation and displayed very similar transcriptome changes both prior and following LIF withdrawal. Consistent with other publications however, we were unable to correlate changes in gene expression with altered DNA methylation. Thus, while transcriptome differences are at best only partially attributable to the DNA demethylase activity of Tet proteins. Nonetheless, the baseline changes in the transcriptome induced by loss of Tet2 in particular are likely at least partially responsible for the inability of these cells to differentiate. Given the isogenic background we utilized and the use of multiple clones for each genotype, our work represents a valuable resource for investigators interested in pluripotency, epigenetics, and the role of Tet proteins in regulating gene expression.
RNA Isolation and RT-qPCR RNA was isolated using TRIzol (Thermo Fisher 15596026) according to the manufacturer's instructions.
Further puri cation of RNA was done using the Qiagen RNeasy Mini Kit (Qiagen 7404) and converted to cDNA with the iScript cDNA Synthesis Kit (Bio-Rad 4106228). Approximately 20ng of cDNA was used for each reaction. Primers used are listed in Table 1. All primers, unless previously published, were designed using mm9.

LIF Withdrawal and Flow Cytometry
Approximately 5,000 ESC were plated on one well of a six-well gelatin-coated plate in ESC media with LIF and replaced the following day with ESC media lacking LIF. Media was changed every day and on day 6 the cells were analyzed using the BD LSR II ow cytometer and analyzed with FlowJo.

Statistical Analysis
Statistically signi cant differences were measured using a two-tailed Student's t-test and a p-value <.05.
Bar graphs represent the mean of all experiments and errors bars are standard error of the mean (SEM).

Next-Generation Sequencing Library Preparation
RNA-seq: 1 µg total RNA was obtained, as described above; ERCC RNA (Thermo 4456740) was then added to each sample prior to Poly-A mRNA selection (NEB E7490). RNA-seq libraries were made using the NEBNext Ultra RNA Library Prep Kit for Illumina (E7530). The following numbers of clones were used for both RNA-seq and RRBS analysis: 3 WT, 3 Tet1-/-, 2 Tet2-/-, 2 DKO, and 2 TKO. RRBS: Libraries were made from 100 ng gDNA using Premium Reduced Representation Bisul te Sequencing Kit (Diagenode C02030033). All libraries were run as paired-end (76 total cycles) on an Illumina NextSeq 500.
The R package methylKit v1.4.1 [31] was used for further analysis. The aligned les from Bismark were utilized to extract the methylation calls occurring at only CpG dinucleotides with minimum of 10 read coverage. Differentially Methylated Regions (DMRs) were identi ed for each comparison between treatment and wild type. DMRs were selected based on q-value < 0.05 and those that meet the minimum percent methylation difference cut-off of 25%. Reported DMRs were then annotated using HOMER v4.10 software [32] and the distribution of DMRs in different genomic elements were plotted. Promoters were de ned as -1kb to +100bp around TSS of RefSeq gene. Intergenic partitions were de ned as genomic regions that did not overlap with promoters, exons and introns.

RNA-seq Analysis
The raw RNA sequence reads were mapped to mouse reference genome build mm9 using STAR v2.5.1 [33] using default parameters and including normalization using the ERCC spike-ins. Quality control matrices were con rmed with a FastQC program. Differential Expression analysis was done using DESeq package in R [34]. Differentially expressed genes were called as signi cant at Benjamin-Hochberg adjusted p-value < 0.05 and fold change of 2. Principal Component Analysis (PCA) was done in R to see the variance between the samples. Heatmaps were done in R using heatmap.2 in R with Pearson method for distance and average agglomeration for clustering. Top 1000 variable genes with log2 read counts were plotted in heatmap after removing duplicates and miRNA genes.

List Of Abbreviations
Tet: Ten Eleven Translocation