Histone Demethylase UTX is an Essential Factor for Zygotic Genome Activation and Regulates Zscan4 Expression in Mouse Embryos

Knockdown of UTX leads to 2-cell retardation

In order to understand the pattern of UTX expression, we collected oocytes, embryos, and different tissues. The qPCR, IF and WB results showed that UTX expression was predominant in the preimplantation embryos, especially in the zygotes and 2-cell stage embryos (Figure 1A, 1B, and 1C).

To explore the role of UTX in early embryonic development, we designed and constructed two short interfering RNAs (siRNAs) specifically targeting Utx. The two siRNAs were mixed in a 1:1 ratio to form Utx-siRNA (Figure S1A). A siRNA without any specificity to Utx or other genes was constructed as a si-control. The results of qPCR analyses revealed that the Utx mRNA expression level decreased by > 90% in the si-Utx-injected embryos compared to that in the si-control-injected embryos (P < 0.01; Figure S1B, S1C). Consistent with the qPCR results, the expression of UTX protein was significantly decreased in the si-Utx-injected embryos as shown using IF and WB (IF, P < 0.01; WB, P < 0.05; Figure 1D, 1E). Furthermore, IF staining showed that the H3K27me3 levels in the si-Utx-injected group were significantly increased (P < 0.05; Figure 1F). Notably, the embryonic development rate from the 4-cell to blastocyst stage was significantly reduced in si-Utx-injected group (Figure 1G). Based on these results, we hypothesized that knockdown of Utx will influenced the occurrence of the ZGA event.

To confirm this hypothesis, we selected 12 ZGA markers and 4 maternal effector genes to detect their expression in 2-cell stage embryos by qPCR. Compared with the si-control group, the expression of Zscan4d and Tcstv1 was significantly down-regulated, whereas the other ZGA-associated genes were not affected (P < 0.05; Figure 1H). The expression of the maternal effector gene Stella was significantly lower than that in the si-control group (P < 0.01; Figure S1D). In addition, IF results showed that the expression of ZSCAN4D protein in the si-Utx-injected group was significantly down-regulated (P < 0.01; Figure 1I). These results indicated that knockdown of UTX impaired embryonic development and resulted in abnormal expression of ZSCAN4D at both mRNA and protein levels.

Overexpression of UTX also results in 2-cell retardation

To determine the effects of UTX overexpression, we constructed an in vitro transcription vector containing Utx N-terminally tagged with Myc tag (Figure S2A, S2B). It is worth noting that the exogenous Myc ectopic expression vector allowed us to track the UTX protein in early embryos without the use of specific antibodies. The qPCR assay showed that Utx mRNA expression in Utx-mRNA-injected group was nearly 6-fold higher than that in the control group at the 2-cell stage (P < 0.01; Figure S2C, S2D). Consistent with the qPCR results, the expression of UTX protein was significantly increased in the Utx-mRNA-injected group as shown using IF and WB (IF, P < 0.001; WB, P < 0.01; Figure 2A, 2B). In addition, IF staining showed that overexpression of UTX markedly reduced the levels of H3K27me3 (P < 0.05; Figure 2A).

The majority of Utx-mRNA-injected embryos arrested at the 2-cell stage, whereas the control group developed to the blastocyst stage (Figure 2C). This phenomenon was consistent with results from the si-Utx-injected group. Compared with the control group, a portion of ZGA-associated genes including MuERVL, Zscan4d, Tcstv1, and Tcstv3 were significantly highly expressed in the Utx-mRNA-injected embryos, as detected by qPCR (MuERVL, Zscan4d, Tcstv1, P < 0.05; Tcstv3, P < 0.01; Figure 2D). The expression of Lin28a, one of the four maternal effector genes, was significantly up-regulated in the Utx mRNA-injected embryos (P < 0.05; Figure S2E). Furthermore, the expression of ZSCAN4D protein was significantly up-regulated in the Utx-mRNA-injected embryos as shown using IF assay, which was consistent with the qPCR result (P < 0.05, Figure 2E). Taken together, these results suggested that overexpression of UTX also impaired embryonic development and led to abnormal expression of ZSCAN4D.

Utx transgenic mice exhibit ZGA defects

Having demonstrated that both knockdown and overexpression of UTX markedly reduced embryonic development rate, we further generated Utx transgenic mice as an in vivo model (referred to as Utx+; Figure S3A, S3B) to explore whether Utx+ mice had similar phenomenon to the in vitro experiment. The body weight of the Utx+ mice showed that neither male nor female mice presented any significant difference compared to that controls, respectively (Figure S3C). When Utx+ males and females were inter-crossed, the number of newborn pups was significantly lower than that from ♂ (Utx+) X ♀ (wild-type) or wild-type inter-crossed groups (Figure 3A). This result indicated that the overexpression of UTX could reduce the fecundity of females. In addition, the number of ovulations in Utx+ females were significantly reduced (P < 0.01; Figure 3B). Furthermore, the size of their ovaries was significantly smaller than that of the controls, and the ovary/body weight ratio also being significantly decreased (P < 0.05; Figure 3C). As expected, the Utx gene significantly highly expressed at the 2-cell embryos as shown by qPCR (P < 0.01; Figure 3D). In the ovaries, the UTX proteins highly expressed in oocytes as shown by immunohistochemistry (IHC) (Figure 3E).

The blastocyst formation rate of the Utx+ embryos was significantly reduced, which was in consistence with the Utx-mRNA-injected embryos (Figure 3F). The qPCR results demonstrated that expression of the ZGA-associated genes Zscan4d and Tcstv3 was significantly up-regulated (P < 0.05; Figure 3G), as was expression of the maternal effector gene Gdf9 (P < 0.05; Figure S3D). These findings suggested that the pattern of gene expression in the Utx+ mice was identical to that in the Utx-mRNA-injected group. Together, the results from Utx-siRNA, Utx-mRNA, and the in vivo Utx+ mice model indicated that Zscan4d represents the target gene of UTX (Figure 3H). To test whether UTX positively regulated Zscan4d in 2-cell embryos, we performed a ChIP-qPCR experiment using an anti-UTX antibody and two specific primer pairs for the promoter (P) and exon 1 (E) of the Zscan4d locus. The results confirmed that UTX does bind to the Zscan4d promoter region (Figure 3I). Furthermore, the UTX occupancy at Zscan4d-P region declined following si-Utx injection (P < 0.01; Figure 3I), supporting the concept that Zscan4d expression can be positively regulated by UTX.

Fine-tuning Zscan4d expression partially rescue si-Utx embryonic development

To further examine the relationship between UTX and Zscan4d, we constructed an in vitro transcription vector wherein Zscan4d was tagged N-terminally with Myc tag (Figure S4A, S4B). The qPCR and IF results showed that the expression of ZSCAN4D was significantly increased in the Zscan4d-mRNA-injected group compared to that in control group (qPCR, P < 0.01; IF, P < 0.001; Figure S4C, S4D).

To explore whether the abnormal embryonic development following UTX modulation is mainly caused by Zscan4d dysregulation, we performed the following rescued experiments in the Utx-knockdown embryos. The Utx-siRNA and Zscan4d-mRNA were simultaneously injected into the zygotes, which showed that the embryonic development significantly improved in the rescued group, especially from 2-cell to 4-cell stage (Figure 4A, Figure S5, Figure S6, and Table S1). However, the blastocyst formation rate did not significantly increase (Figure 4A, Figure S5, Figure S6, and Table S1), and the cell numbers of blastocysts was with no difference between different injected groups (Figure 4B, 4C). To gain further insights into blastocyst lineage segregation, the blastocysts derived from different group were subjected to IF staining of NANOG and CDX2 (Figure 4D). The results showed that the NANOG and CDX2 were exclusively localized to the nuclei of the inner cell mass (ICM) and trophectoderm (TE) in the control group, while the NANOG signals was mislocalized to the ICM in the Utx-siRNA, Utx-mRNA, and Utx+ groups (Figure 4D). These results indicated that fine-tuning expression of Zscan4d partially rescue the impaired embryonic development upon UTX depletion.

Ethics statement

All studies adhered to procedures are consistent with the National Research Council Guide for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committee at Inner Mongolia University.

Animals

B6D2F1 (C57BL/6 × DBA/2) female and male mice were purchased at 6-8 weeks of age from Inner Mongolia University (China). Pseudopregnant CD1 mice were used as embryo recipients. All mice were reared in house under specific-pathogen-free conditions and were housed under controlled lighting conditions (light: 08:00~20:00). The animals had free access to food and water. The mice were randomly allocated to each experimental group.

Plasmids and antibodies

The inserts of plasmids pcs2-UTX and pcs2-ZSCAN4D used in this study were PCR-amplified from the cDNA of NIH-3T3 cells and 2-cell mouse embryos, respectively. The PCR products were cloned into the Xho I and Xba I (Thermo, USA) sites downstream of the Myc ectopic tag. The antibodies covered were rabbit anti-UTX (Millipore, USA, for IF and WB), rabbit anti-H3K27me3 (Abcam, USA, for IF), mouse anti-Myc (Santa Cruz, USA, for IF), rabbit anti-ZSCAN4 (Millipore, for IF), rabbit anti-CDX2 (Abcam, for IF), mouse anti-NANOG (Abcam, for IF), rabbit α-tubulin (Proteintech, USA, for WB), Alexa Fluor 594 and Alexa Fluor 488 (Thermo, for IF), anti-UTX (KDM6A) (Abcam, for ChIP), anti-IgG (Abcam, for ChIP).

Transgenic mice generation

Transgenic mice were produced on a B6D2F1 background. For pronuclear microinjection, the pcs2-UTX expression vector was linearized using Sal I and Not I (Thermo). The purified DNA was diluted in 1×TE buffer (10 ng/μl) and microinjected into the well-recognized pronuclei of zygotes. The injected zygotes were transferred into the oviducts of pseudopregnant females. Transgenic mice were identified by PCR analysis using specific primers (Table S2).

Collection of oocytes and fertilized embryos

Superovulation was performed as previously described ¹⁹. The GV stage oocytes were isolated from B6D2F1 female (6-8 weeks) mice 46 h following pregnant mare serum gonadotropin (PMSG, Sansheng, China) injection. The MII stage oocytes were isolated from B6D2F1 female (6-8 weeks) mice 16 h following human chorionic gonadotropin (hCG, Sansheng) injection. Zygotes (24 h after hCG) were obtained from the ampulla of the superovulated female mice after mating with male mice. Then, 2-cell (44 h after hCG), 4-cell (56 h after hCG), 8-cell (67 h after hCG), morula (76 h after hCG), and blastocyst (97 h after hCG) embryos were obtained by culturing zygotes in KSOM medium.

In vitro mRNA transcription, siRNA construction, and microinjection in zygotes

Utx cDNAs and Zscan4d cDNAs were cloned into SP6-driven vectors. To prepare mRNAs for microinjection, Utx and Zscan4d expression vectors were linearized by Sal I and Not I (Thermo), the constructs were then purified with phenol-chloroform extraction and ethanol precipitation. The mRNA was synthesized by in vitro transcription using a mMESSAGE mMACHINE SP6 Ultra Kit (Thermo) according to the manufacturer's instructions. The synthesized mRNA was purified by lithium chloride precipitation and diluted with nuclease-free water. For mRNA injection, zygotes were injected with 100 ng/μl mRNA. The Utx-siRNA-1 and Utx-siRNA-2 respectively located at 585 and 1,518 bp downstream from the start codon (4390771- s75838/s75839, Thermo). The commercially available siRNA without any specificity to the known genes was used as control (4390843, Thermo). The siRNAs were diluted to 20 μM in nuclease-free water and stored at -80°C until use. As previously described ¹² with minor modifications, the zygotes were transferred to the manipulation M2 media and the mRNA or siRNA were microinjected into the zygotes with approximately 10 pl volume using FemtoJet (Eppendorf, German) micromanipulator. Following injection, the zygotes were kept at room temperature for 15 min and then incubated in the incubator.

RNA extraction and quantitative real-time PCR

Total RNA was extracted from each pool of embryos (n = 3 pools of 20 embryos) using the Pico-Pure RNA Isolation Kit (Thermo). The reverse transcription-qPCR (RT-qPCR) was transcribed from total RNA using the SuperScript III First-Strand Synthesis System (Thermo) according to the manufacturer's instructions. All primers were listed in Table S2. The quantitative real-time PCR (qPCR) was performed using a TB Green Premix Ex Taq (Takara, Japan), and the signals were detected with ABI7500 Real-Time PCR System (Thermo). Analyses of the relative gene expression was measured using the 2-ΔΔCt method.

Immunofluorescence staining and immunohistochemistry

As previously described ¹², the immunofluorescence (IF) staining of embryos conducted with minor modifications. The embryos were rinsed three times in phosphate-buffered saline (PBS) with 0.3% (m/v) bovine serum albumin (BSA), fixed with 4% paraformaldehyde for 1 h and then permeabilized with 0.5% Triton X-100 for 15 min at room temperature. The embryos were blocked with 3% (m/v) BSA in PBS for 1 h and then incubated with the primary antibodies for 1 h at 37°C. After washing three times with PBS containing 0.3% BSA, the embryos were incubated with the secondary antibody for 1 h at 37°C. The nuclei were stained with 10 μg/ml DAPI (4',6-diamidino-2-phenylindole) (Sigma, USA). The embryos were then mounted on glass slides and examined with a Confocal Laser-Scanning Microscope (A1R, Nikon, Japan). The steps of immunohistochemistry (IHC) were as follows. The ovaries from 8-week-old mice were fixed in 4% paraformaldehyde (Sigma) overnight at 4 °C, then dehydrated through a series of graded ethanol solutions and xylene, and embedded in paraffin (Amresco, USA). The ovaries were sectioned at a thickness of 5 μm, incubated with 1:250 anti-Myc antibody after deparaffinization at 4°C by overnight. The sections were stained with hematoxylin (Amresco) and imaged under an optical microscope (80i, Nikon, Japan).

Western blot

Whole cell lysates from 500 embryos were lysed in 2× Laemmli buffer (0.5 M Tris HCl, pH 6.8, 0.4% SDS, 2% glycerol, 0.5% β-mercaptoethanol) and Bromophenol Blue (Bio-Rad, USA). The lysates were heated at 100 °C for 10 min, then cooled rapidly for 1 min. Proteins were separated by using SDS polyacrylamide gel electrophoresis (90 V for 30 min, 110 V 100 min,) and blotted onto nitrocellulose membranes (Biosharp, China). The membranes were blocked in TBST (TBS containing 0.1% Tween 20) containing 5% non-fat milk at room temperature 1 h, followed by overnight incubation with antibody at 4°C. After being probed with the primary antibodies, the membranes were washed in TBST, incubated with an HRP-linked secondary antibody for 1 h at 37°C, and washed three times with TBST. The bound antibodies were detected with SuperSignal West Femto Substrate Trial Kit (Thermo). The intensity of the protein bands was calculated by ImageJ software (NIH, Bethesda, USA).

ChIP-qPCR analysis

The ChIP assays were performed according to the protocol of ultra-low-input chromatin immunoprecipitation (ULI-ChIP) ⁴¹. Samples of 2-cell embryos of si-Control and si-Utx-injected group were harvested for ChIP. 500 embryos were used for each reaction and three replicates were performed. 2 μg of either anti-UTX (KDM6A) or anti-IgG was used for each immunoprecipitation (IP) reaction. ChIP-qPCR were detected using a TB Green Premix Ex Taq (Takara), and the signals were detected with ABI7500 Real-Time PCR System (Thermo). Primers used for ChIP-qPCR were listed in Table S2.

Statistical Analysis

Independent t-tests were performed using SPSS software version 22.0 (SPSS Inc., Chicago, IL, USA) to compare difference between two groups. A P-value < 0.05 was considered significant.