Tet1 controls meiosis by regulating meiotic gene expression

Meiosis is a germ cell-specific cell division process through which haploid gametes are produced for sexual reproduction1. Prior to initiation of meiosis, mouse primordial germ cells (PGCs) undergo a series of epigenetic reprogramming steps2,3, including global erasure of DNA methylation on the 5-position of cytosine (5mC) at CpG4,5. Although several epigenetic regulators, such as Dnmt3l, histone methyltransferases G9a and Prdm9, have been reported to be critical for meiosis6, little is known about how the expression of meiotic genes is regulated and how their expression contributes to normal meiosis. Using a loss of function approach, here we demonstrate that the 5mC-specific dioxygenase Tet1 plays an important role in regulating meiosis in mouse oocytes. Tet1 deficiency significantly reduces female germ cell numbers and fertility. Univalent chromosomes and unresolved DNA double strand breaks are also observed in Tet1-deficient oocytes. Tet1 deficiency does not greatly affect the genome-wide demethylation that takes place in PGCs but leads to defective DNA demethylation and decreased expression of a subset of meiotic genes. Our study thus establishes a function for Tet1 in meiosis and meiotic gene activation in female germ cells.

of DNA methylation 4,5 . Recent demonstration that the Tet family proteins are involved in DNA demethylation prompted us to evaluate the role of the Tet proteins in PGC reprogramming [7][8][9][10] . RT-qPCR analysis demonstrated that Tet1 is preferentially expressed in PGCs, while Tet2 is expressed in both PGCs and somatic cells, and Tet3 is mainly expressed in somatic cells during PGC development (E9.5-E13.5) (Supplementary Fig. 1). To explore a potential role of Tet1 in PGC reprogramming and/or germ cell development, we generated Tet1 gene-trap mice ( Supplementary Fig. 2a-d). Homozygous mutant mice (Tet1 Gt/Gt ) were generated by crossing heterozygous mice. Southern blot analysis, genomic sequencing, and allele-specific PCR confirmed a single site insertion of a tandem-repeated gene-trap cassette into the first intron of Tet1 gene ( Supplementary Fig. 2a-d). β-galactosidase activity analysis revealed that the transgene is almost exclusively expressed in PGCs ( Supplementary Fig.  2e), which is consistent with the Tet1 expression pattern ( Supplementary Fig. 1), supporting a single locus insertion. Western blot analysis demonstrated that the insertion nullified the expression of the full-length Tet1 protein ( Supplementary Fig. 3a). As expected, a fusion protein between the first exon of Tet1 (aa1-621) and βGeo (1303aa) was detected in Tet1 Gt/Gt ES cell. Consistent with loss of Tet1, dot-blot and mass spectrometry analyses revealed about 45% reduction of 5hmC level in E9.5 Tet1 Gt/Gt embryos (Supplementary Fig.  3b-d). RT-qPCR analysis of E11.5 PGCs demonstrated that Tet1 level is less than 5% that of the wild-type level ( Supplementary Fig. 3e). Consistently, Tet1 protein was also not detectable in spreads of Tet1 Gt/Gt PGCs ( Supplementary Fig. 3f). Furthermore, immunostaining revealed loss of the dotted 5hmC staining signal in germ cells of E14.5 Tet1 Gt/Gt genital ridge ( Supplementary Fig. 3g). Collectively, these data indicate that Tet1 expression is effectively abolished and 5hmC level was significantly reduced in Tet1 Gt/Gt PGCs.
Analysis of early backcross generations (N1-2) revealed embryonic lethal for homozygous while heterozygous mice were born in Mendelian ratio (Supplementary Table 1). Although further backcross (N3-6) relieved the embryonic lethal phenotype, the number of viable homozygous mice is still only approximately 1/3 of the numbers expected. Since severity of embryonic lethality is affected by genetic background, we only used later than N6 generation of Tet1 Gt/Gt mice for subsequent analysis. Similar to a recent report 11 , reduced pup numbers were observed when either homozygous male or female animals were crossed with wild-type animals, and the pup numbers are even fewer when homozygous animals were crossed ( Supplementary Fig. 4a). Because male gonad was morphologically normal and no obvious defects in male germ cell development were observed (data not shown), our attention is focused on characterizing the female germ cell phenotypes.
We found that the size of Tet1 Gt/Gt ovary is significantly smaller, with a 30% reduction in the ovary to body weight ratio (Figs. 1a, Supplementary Fig. 4b). Interestingly, ovary agenesis-caused asymmetric ovary size is frequently observed in the Tet1 Gt/Gt animals (Fig.  1a). Both fully-grown oocytes in the ovary and ovulated oocytes after hormonal stimulation are significantly reduced in the Tet1 Gt/Gt animals (Fig. 1b, c). Ovary staining with germ cellspecific markers (Mvh, TRA98, and Msy2) followed by counting revealed a significantly reduction of oocyte number from E16.5-E18.5 (Fig. 1d, Supplementary Fig. 5a), concurrent with a significant increase in apoptotic oocytes ( Supplementary Fig. 5b, c). These results suggest that increased apoptosis is likely one contributing factor for oocyte loss in Tet1 Gt/Gt animals.
Since meiotic prophase takes place in the embryonic ovary and that meiotic defects can cause germ cell apoptosis 12 , we asked whether loss function of Tet1 might lead to a meiotic defect. Immunostaining of oocyte surface spreads with a meiotic marker SYCP3 and a centromere-marker CREST revealed that progression of meiotic prophase is severely impaired in Tet1 Gt/Gt animals. About 50% of oocytes remained at the leptotene stage and no oocyte reached the pachytene stage in E16.5 Tet1 Gt/Gt ovaries (Fig. 2a). Although pachytene stage oocytes were observed in E17.5 and E18.5 Tet1 Gt/Gt ovaries, their percentage is significantly reduced throughout the developmental stages analyzed, indicating a developmental block rather than a developmental delay. This notion is also supported by the fact that Tet1 depletion does not affect the percentage of diplotene stage oocytes at E18.5 (Fig. 2a), and no significant difference in the body size between Tet1 +/Gt and Tet1 Gt/Gt embryos at E16.5 and E18.5 ( Supplementary Fig. 6). The significant decrease of pachytene stage oocytes in Tet1 Gt/Gt ovaries suggests that Tet1 deficiency may cause aberrant synapsis formation. Considering the increased apoptosis in E16.5-E18.5 Tet1 Gt/Gt ovary (Supplemental Fig. 5b, c), severely affected germ cells might be eliminated, explaining the dramatic reduction in germ cell numbers in E18.5 Tet1 Gt/Gt ovaries (Fig. 1d). Thus, meiotic defects are likely the cause of germ cell reduction in Tet1 Gt/Gt ovaries.
In normal meiosis, the axial element (AE) of synaptonemal complex (SC) starts to form at the leptotene stage, AE-alignment and synapsis formation are initiated at the zygotene stage and completed at the pachytene stage ( Supplementary Fig. 7). However, loss function of Tet1 resulted in an increase in unpaired SC at zygotene stage oocytes (Fig. 2b, Supplementary Fig. 7). Despite the presence of Sycp3-positive AE, about 30% of zygonema contained less than 4 SYCP1-positive transverse filaments (TF) in Tet1 Gt/Gt oocytes (Fig.  2b) suggesting a synapsis formation defect. Defects in synapsis formation were also observed in the pachytene stage Tet1 Gt/Gt oocytes. In these oocytes, the great majority of their AE failed to align and remain separated from each other as univalent chromosomes (Fig. 2c, Supplementary Fig. 7). Continuous SYCP3 in short-stretches of SC indicates that those oocytes were at a stage corresponding to pachytene. Quantification of pachytene stage oocytes revealed that 92% of them have 0-4 univalent chromosomes in wild-type oocytes, but that number dropped to 62% in the Tet1 Gt/Gt oocytes, which accompanies the increase of the oocytes with 5 or more univalent chromosomes (Fig. 2c). These results indicate that loss function of Tet1 impaired synapsis formation.
Since about 62% of pachytene stage oocytes are successful for AE pairing, we asked whether they exhibited other defects. At the early phase of meiosis, DNA double strand breaks (DSBs) are introduced in the initiation of homologous recombination, which can be detected by the presence of γH2AX 1 . As expected, immunostaining revealed the presence of γH2AX throughout the nucleoplasm from leptotene to zygotene stages in both wild-type and Tet1 Gt/Gt oocytes ( Supplementary Fig. 8). However, while γH2AX is gradually decreased around pachytene stage and only few foci remained associated with fully synapsed chromosome cores in wild-type oocytes, cloud-like nuclear staining of γH2AX remained in the Tet1 Gt/Gt pachytene and even early diplotene stage oocytes marked by crossover-specific marker MLH1 13 (Fig. 2d, Supplementary Fig. 8). Analysis of γH2AX staining pattern and quantification of the three categories (negative, partially positive, and positive) (Supplementary Fig. 9) clearly demonstrated that Tet1 depletion caused accumulation of γH2AX in pachytene and early diplotene oocytes (Fig. 2e).
The presence of γH2AX in late stage meiotic oocytes was also confirmed by co-staining of γH2AX with late stage meiosis marker Msy2 at E18.5 ovaries (Supplementary Fig. 10).
Importantly, significant increase in γH2AX and cleaved Caspase3 double positive-cells were also observed ( Supplementary Fig. 10b, c), suggesting that increase in apoptotic cell death is likely caused by meiotic defects. Consistent with a DSB repair defect, the DSB repair associated recombinase RAD51 1 remains in pachytene and diplotene oocytes ( Supplementary Fig. 11). The presence of γH2AX and delayed removal of RAD51 in the chromosomes indicate that homologous recombination is impaired in Tet1 Gt/Gt oocytes. Staining with the crossover mark MLH1 indicated that the MLH1 foci numbers are significantly reduced in Tet1 Gt/Gt pachytene and early diplotene oocytes ( Supplementary  Fig. 12), further support a homologous recombination defect. Collectively, the above results support the notion that loss function of Tet1 leads to meiotic defects that include univalent chromosome formation, as well as DSB repair and homologous recombination defects.
Previous studies have shown that establishment of proper pericentric heterochromatin (PCH) structure plays an important role in meiosis 14 . Interestingly, immunostaining revealed specific enrichment of 5hmC at PCH of many prophase meiotic chromosomes and this enrichment is eliminated in Tet1 Gt/Gt oocytes ( Supplementary Fig. 13). Since prophase meiotic PCH possesses specific histone modification pattern 15 , we asked whether loss function of Tet1 can affect PCH structure by affecting PCH histone modification pattern. Immunostaining revealed that loss function of Tet1 neither affected PCH histone modification pattern ( Supplementary Fig. 14), nor centromere clustering or HP1γ localization ( Supplementary Figs. 14, 15). Collectively, these results suggest that the observed meiotic defect is unlikely due to a defect in PCH.
The stage-and cell type-specific Tet1 expression pattern ( Supplementary Fig. 1) and the dual roles of Tet1 in transcription 16,17 suggest that the meiotic defects in Tet1 Gt/Gt germ cells might be due to aberrant transcription in PGCs. Thus, we purified PGCs from female E13.5 embryos and profiled their transcriptome by mRNA sequencing. We choose to focus on this time point because this is the time when female PGCs enter meiotic prophase after epigenetic reprogramming. We used a recently developed Smart-Seq method 18 and generated over 20 million unique reads per sample which allowed us to identify over 13 thousands expressed transcripts in each genotype (Supplementary table 2). Hierarchical clustering and global correlation analysis indicated that the samples were clearly separated by their genotypes with Spearman correlation coefficient of 0.98/0.99 within biological replicates ( Supplementary Fig. S16). Depletion of Tet1 resulted in differential expression of 1,010 genes (FDR<0.05), among which over 80% (899 genes) were down regulated (Fig. 3a,  Supplementary table 3). Gene ontology (GO) analysis revealed that the most significantly enriched pathways of these down-regulated genes are related to cell cycle (p-value <9e-11) and meiosis-related processes (p-value <2e-6) (Fig. 3b, Supplemental Fig. S17a). In contrast, no significant enrichment of pathways or biological processes was identified in the up-regulated gene group. Importantly, genes known to be critical for meiosis are downregulated in Tet1 Gt/Gt PGCs (Fig. 3a, Supplementary table 3). These genes include Stra8, Prdm9, Sycp1, Mael and Sycp3 1, [19][20][21][22] . Notably, this set of meiotic genes remained down regulated even at later developmental stage (Supplemental Fig. S17b, c), consistent with the meiotic defects observed in E16.5 Tet1 Gt/Gt PGCs. The effect of Tet1 on the expression of at least a subset of these genes are direct as ChIP analysis demonstrated that Tet1 occupies Sycp1, Mael, and Sycp3 gene promoters (Fig. 3c). Thus, Tet1 loss directly contributes to aberrant regulation of at least a subset of meiotic genes in PGCs.
To investigate how Tet1 might be involved in activation of these meiotic genes, we performed whole-genome bisulfite sequencing (WGBS) using an ultra-low input, Tn5mCseq method 23 . We generated 945 million and 302 million reads for Tet1 Gt/Gt and wild-type PGCs, respectively. After removing clonal reads due to limited input cells, we obtained 14-16 million CpG sites per genotype at 1.76-2.66x genome coverage (Supplementary table 4), which is over100-fold higher than a previous effort 24 . To our acknowledge, this provides the most comprehensive methylation map in PGCs so far. Consistent with previous findings 24 , we found that PGCs are globally hypomethylated (Fig. 4a), which is verified by immunostaining (Supplemental Fig. S18). Despite no dramatic increase in global DNA methylation, DNA methylation level is generally higher in mutant PGCs, particularly in exons, introns, LTRs, and IAPs (Fig. 4a, p<0.01). With the caveat that 2X genome coverage may not allow robust identification of differentially methylated regions (DMRs) between Tet1 Gt/Gt and wild-type PGCs, we nevertheless performed detailed analysis and identified 4,337 putative DMRs (Supplementary table 5, Supplemental Fig. S19) that are mostly located far away from transcriptional start sites (Supplemental Fig. S20a). These putative DMRs are associated with 5,242 genes, among which 255 exhibited altered expression (Fig.  4b, Supplementary table 6) and are enriched for cell cycle regulation as well as for reproductive and infertility processes (FDR=0.02) (Supplementary Fig. S20b, c). Furthermore, these DMRs are enriched for Tet1 binding in mouse ES cells 17 (P-value <10E-100). To evaluate whether Tet1 binding affects DNA methylation in PGCs, we performed bisulfite sequencing on the three verified Tet1 target genes (Fig. 3c). Consistent with the involvement of Tet1 in DNA demethylation, the methylation level of the Sycp1, Mael, and Sycp3 promoters is increased in the Tet1 Gt/Gt PGCs compared with that in the wild-type PGCs (Fig. 4c). These data indicate that Tet1-mediated demethylation of these genes is likely involved in their activation during PGC development. We note, some down-regulated genes, such as Stra8, showed no obvious change in DNA methylation indicating that they are either regulated indirectly by Tet1 or in a DNA methylation-independent manner.
Taken together, our study provides the first evidence that Tet1 is not responsible for global demethylation in PGCs, rather it plays a specific role in meiotic gene activation at least partly through DNA demethylation. Depletion of Tet1 leads to down-regulation of meiotic genes, which causes defective meiotic prophase including accumulation of non-repaired DSBs, and formation of univalent chromosomes. The meiotic defects cause loss of oocytes and consequent decrease in fertility and small litter size. Previous studies have established that DNA methylation levels of certain meiotic genes are decreased concomitant with genomic reprogramming 25 . Our study has extended this observation by demonstrating that Tet1 mediates locus-specific demethylation and subsequent activation of a subset of meiotic genes, revealing a specific function of Tet1 in germ cell development.

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.   (c) Bisulfite sequencing analysis of the Sycp1, Mael, and Sycp3 gene promoters of Tet1 binding site in wild-type and Tet1 Gt/Gt PGCs. Each CpG is represented by a circle. Open and filled circles represent unmethylated or methylated, respectively. Percentages of DNA methylation are indicated.