Effects of transgene insertion loci and copy number on Dnmt3L gene silencing through antisense transgene-derived PIWI-interacting RNAs.

PIWI-interacting RNAs (piRNAs), which are germ cell-specific small RNAs, are essential for spermatogenesis. In fetal mouse germ cells, piRNAs are synthesized from sense and antisense RNAs of transposable element sequences for retrotransposon silencing. In a previous study, we reported that transgenic mice expressing antisense-Dnmt3L under the control of the Miwi2 promoter (Tg-Miwi2P-asDnmt3L) exhibited piRNA-mediated DNMT3L down-regulation. In this study, two transgene integration loci (B3 and E1) were identified on chromosome 18 of the Tg-Miwi2P-asDnmt3L mice; these loci were weak piRNA clusters. Crossbreeding was performed to obtain mice with the transgene cassette inserted into a single locus. DNMT3L was silenced and spermatogenesis was severely impaired in mice with the transgene cassette inserted at the B3 locus (Tg-B mice). In contrast, spermatogenesis in mice bearing the transgene at the E1 locus (Tg-E mice) was normal. The number of piRNAs for Dnmt3L in Tg-B mice was eightfold higher than that in Tg-E mice. Therefore, both gene silencing and impaired spermatogenesis depended on the transgene copy number rather than on the insertion loci. Additionally, the endogenous Dnmt3L promoter was not methylated in Tg mice, suggesting that Dnmt3L silencing was caused by post-transcriptional gene silencing. Based on these data, we discuss a piRNA-dependent gene silencing mechanism against novel gene insertions.

PNLDC1, mouse pre-piRNA Trimmer, is required for meiotic and post-meiotic male germ cell development.

PIWI-interacting RNAs (piRNAs) are germ cell-specific small RNAs essential for retrotransposon gene silencing and male germ cell development. In piRNA biogenesis, the endonuclease MitoPLD/Zucchini cleaves long, single-stranded RNAs to generate 5' termini of precursor piRNAs (pre-piRNAs) that are consecutively loaded into PIWI-family proteins. Subsequently, these pre-piRNAs are trimmed at their 3'-end by an exonuclease called Trimmer. Recently, poly(A)-specific ribonuclease-like domain-containing 1 (PNLDC1) was identified as the pre-piRNA Trimmer in silkworms. However, the function of PNLDC1 in other species remains unknown. Here, we generate Pnldc1 mutant mice and analyze small RNAs in their testes. Our results demonstrate that mouse PNLDC1 functions in the trimming of both embryonic and post-natal pre-piRNAs. In addition, piRNA trimming defects in embryonic and post-natal testes cause impaired DNA methylation and reduced MIWI expression, respectively. Phenotypically, both meiotic and post-meiotic arrests are evident in the same individual Pnldc1 mutant mouse. The former and latter phenotypes are similar to those of MILI and MIWI mutant mice, respectively. Thus, PNLDC1-mediated piRNA trimming is indispensable for the function of piRNAs throughout mouse spermatogenesis.

Mouse GTSF1 is an essential factor for secondary piRNA biogenesis.

The piRNA pathway is a piRNA-guided retrotransposon silencing system which includes processing of retrotransposon transcripts by PIWI-piRNAs in secondary piRNA biogenesis. Although several proteins participate in the piRNA pathway, the ones crucial for the cleavage of target RNAs by PIWI-piRNAs have not been identified. Here, we show that GTSF1, an essential factor for retrotransposon silencing in male germ cells in mice, associates with both MILI and MIWI2, mouse PIWI proteins that function in prospermatogonia. GTSF1 deficiency leads to a severe defect in the production of secondary piRNAs, which are generated from target RNAs of PIWI-piRNAs. Furthermore, in Gtsf1 mutants, a known target RNA of PIWI-piRNAs is left unsliced at the cleavage site, and the generation of secondary piRNAs from this transcript is defective. Our findings indicate that GTSF1 is a crucial factor for the slicing of target RNAs by PIWI-piRNAs and thus affects secondary piRNA biogenesis in prospermatogonia.

GPAT2 is required for piRNA biogenesis, transposon silencing, and maintenance of spermatogonia in mice†.

PIWI-interacting RNAs (piRNAs), a subclass of germ cell-specific noncoding small RNAs, are essential for de novo DNA methylation of retrotransposon genes in embryonic testes. PIWIL2/MILI, one of three mouse PIWI family members, is indispensable for piRNA production, DNA methylation of retrotransposons presumably via piRNA, and normal spermatogenesis. In vitro analysis using germline stem cells (GS cells) revealed that glycerol-3-phosphate acyltransferase 2 (GPAT2), which is a mitochondrial outer membrane protein involved in generation of lysophosphatidic acid (LPA) and highly expressed in testes, plays important roles in spermatogenesis. Namely, GPAT2 binds to PIWIL2 and is closely involved in the biogenesis of piRNAs; this process is independent of its enzymatic activity on LPA. However, GS cells recapitulate only a limited phase of spermatogenesis and the biological functions of GPAT2 remain largely unknown. In this study, we generated GPAT2-deficient mice and conducted comprehensive analyses. The deficient mice showed defective piRNA production and subsequent de-silencing of IAP and Line-1 retrotransposons in fetal testes. In addition, apoptosis of pachytene spermatocytes was observed. These abnormalities were all common to the phenotype of PIWIL2-deficient mice, in which piRNA production was impaired. GPAT2-deficient mice exhibited apoptosis in spermatogonia at the neonatal stage, which was not observed in PIWIL2-deficient mice. These data show that GPAT2 plays a critical role in preventing apoptosis in spermatogonia.

Roles of MIWI, MILI and PLD6 in small RNA regulation in mouse growing oocytes.

The mouse PIWI-interacting RNA (piRNA) pathway produces a class of 26-30-nucleotide (nt) small RNAs and is essential for spermatogenesis and retrotransposon repression. In oocytes, however, its regulation and function are poorly understood. In the present study, we investigated the consequences of loss of piRNA-pathway components in growing oocytes. When MILI (or PIWIL2), a PIWI family member, was depleted by gene knockout, almost all piRNAs disappeared. This severe loss of piRNA was accompanied by an increase in transcripts derived from specific retrotransposons, especially IAPs. MIWI (or PIWIL1) depletion had a smaller effect. In oocytes lacking PLD6 (or ZUCCHINI or MITOPLD), a mitochondrial nuclease/phospholipase involved in piRNA biogenesis in male germ cells, the piRNA level was decreased to 50% compared to wild-type, a phenotype much milder than that in males. Since PLD6 is essential for the creation of the 5΄ ends of primary piRNAs in males, the presence of mature piRNA in PLD6-depleted oocytes suggests the presence of compensating enzymes. Furthermore, we identified novel 21-23-nt small RNAs, termed spiRNAs, possessing a 10-nt complementarity with piRNAs, which were produced dependent on MILI and independent of DICER. Our study revealed the differences in the biogenesis and function of the piRNA pathway between sexes.

MVH in piRNA processing and gene silencing of retrotransposons.

VASA is an evolutionarily conserved RNA helicase essential for germ cell development. The mouse PIWI family proteins MILI and MIWI2 are involved in production of Piwi-interacting RNAs (piRNAs) in fetal male germ cells through a ping-pong amplification cycle. Expression of retrotransposons is elevated in MILI- and MIWI2-deficient male germ cells due to defective de novo DNA methylation, which is presumably caused by impaired piRNA expression. Here, we report that essentially the same abnormalities are observed in MVH (mouse VASA homolog)-deficient mice. Comprehensive analysis of piRNAs in MVH-deficient fetal male germ cells showed that MVH plays crucial roles in the early phase of the ping-pong amplification cycle.

piRNAs derived from ancient viral processed pseudogenes as transgenerational sequence-specific immune memory in mammals.

Endogenous bornavirus-like nucleoprotein elements (EBLNs) are sequences within vertebrate genomes derived from reverse transcription and integration of ancient bornaviral nucleoprotein mRNA via the host retrotransposon machinery. While species with EBLNs appear relatively resistant to bornaviral disease, the nature of this association is unclear. We hypothesized that EBLNs could give rise to antiviral interfering RNA in the form of PIWI-interacting RNAs (piRNAs), a class of small RNA known to silence transposons but not exogenous viruses. We found that in both rodents and primates, which acquired their EBLNs independently some 25-40 million years ago, EBLNs are present within piRNA-generating regions of the genome far more often than expected by chance alone (ℙ = 8 × 10(-3)-6 × 10(-8)). Three of the seven human EBLNs fall within annotated piRNA clusters and two marmoset EBLNs give rise to bona fide piRNAs. In both rats and mice, at least two of the five EBLNs give rise to abundant piRNAs in the male gonad. While no EBLNs are syntenic between rodent and primate, some of the piRNA clusters containing EBLNs are; thus we deduce that EBLNs were integrated into existing piRNA clusters. All true piRNAs derived from EBLNs are antisense relative to the proposed ancient bornaviral nucleoprotein mRNA. These observations are consistent with a role for EBLN-derived piRNA-like RNAs in interfering with ancient bornaviral infection. They raise the hypothesis that retrotransposon-dependent virus-to-host gene flow could engender RNA-mediated, sequence-specific antiviral immune memory in metazoans analogous to the CRISPR/Cas system in prokaryotes.

Targeted gene silencing in mouse germ cells by insertion of a homologous DNA into a piRNA generating locus.

In germ cells, early embryos, and stem cells of animals, PIWI-interacting RNAs (piRNAs) have an important role in silencing retrotransposons, which are vicious genomic parasites, through transcriptional and post-transcriptional mechanisms. To examine whether the piRNA pathway can be used to silence genes of interest in germ cells, we have generated knock-in mice in which a foreign DNA fragment was inserted into a region generating pachytene piRNAs. The knock-in sequence was transcribed, and the resulting RNA was processed to yield piRNAs in postnatal testes. When reporter genes possessing a sequence complementary to portions of the knock-in sequence were introduced, they were greatly repressed after the time of pachytene piRNA generation. This repression mainly occurred at the post-transcriptional level, as degradation of the reporter RNAs was accelerated. Our results show that the piRNA pathway can be used as a tool for sequence-specific gene silencing in germ cells and support the idea that the piRNA generating regions serve as traps for retrotransposons, enabling the host cell to generate piRNAs against active retrotransposons.

The nuage mediates retrotransposon silencing in mouse primordial ovarian follicles.

Mobilization of endogenous retrotransposons can destabilize the genome, an imminent danger during epigenetic reprogramming of cells in the germline. The P-element-induced wimpy testis (PIWI)-interacting RNA (piRNA) pathway is known to silence retrotransposons in the mouse testes. Several piRNA pathway components localize to the unique, germline structure known as the nuage. In this study, we surveyed mouse ovaries and found, for the first time, transient appearance of nuage-like structures in oocytes of primordial follicles. Mouse vasa homolog (MVH), Piwi-like 2 (PIWIL2/MILI) and tudor domain-containing 9 (TDRD9) are present in these structures, whereas aggregates of germ cell protein with ankyrin repeats, sterile alpha motif and leucine zipper (GASZ) localize separately in the cytoplasm. Retrotransposons are silenced in primordial ovarian follicles, and de-repressed upon reduction of piRNA expression in Mvh, Mili or Gasz mutants. However, these null-mutant females, unlike their male counterparts, are fertile, uncoupling retrotransposon activation from sterility.

HSP90α plays an important role in piRNA biogenesis and retrotransposon repression in mouse.

HSP90, found in all kingdoms of life, is a major chaperone protein regulating many client proteins. We demonstrated that HSP90α, one of two paralogs duplicated in vertebrates, plays an important role in the biogenesis of fetal PIWI-interacting RNAs (piRNA), which act against the transposon activities, in mouse male germ cells. The knockout mutation of Hsp90α resulted in a large reduction in the expression of primary and secondary piRNAs and mislocalization of MIWI2, a PIWI homolog. Whereas the mutation in Fkbp6 encoding a co-chaperone reduced piRNAs of 28-32 nucleotides in length, the Hsp90α mutation reduced piRNAs of 24-32 nucleotides, suggesting the presence of both FKBP6-dependent and -independent actions of HSP90α. Although DNA methylation and mRNA levels of L1 retrotransposon were largely unchanged in the Hsp90α mutant testes, the L1-encoded protein was increased, suggesting the presence of post-transcriptional regulation. This study revealed the specialized function of the HSP90α isofom in the piRNA biogenesis and repression of retrotransposons during the development of male germ cells in mammals.

GPAT2, a mitochondrial outer membrane protein, in piRNA biogenesis in germline stem cells.

piRNA (PIWI-interacting RNA) is a germ cell-specific small RNA in which biogenesis PIWI (P-element wimpy testis) family proteins play crucial roles. MILI (mouse Piwi-like), one of the three mouse PIWI family members, is indispensable for piRNA production, DNA methylation of retrotransposons presumably through the piRNA, and spermatogenesis. The biogenesis of piRNA has been divided into primary and secondary processing pathways; in both of these MILI is involved in mice. To analyze the molecular function of MILI in piRNA biogenesis, we utilized germline stem (GS) cells, which are derived from testicular stem cells and possess a spermatogonial phenotype. We established MILI-null GS cell lines and their revertant, MILI-rescued GS cells, by introducing the Mili gene with Sendai virus vector. Comparison of wild-type, MILI-null, and MILI-rescued GS cells revealed that GS cells were quite useful for analyzing the molecular mechanisms of piRNA production, especially the primary processing pathway. We found that glycerol-3-phosphate acyltransferase 2 (GPAT2), a mitochondrial outer membrane protein for lysophosphatidic acid, bound to MILI using the cells and that gene knockdown of GPAT2 brought about impaired piRNA production in GS cells. GPAT2 is not only one of the MILI bound proteins but also a protein essential for primary piRNA biogenesis.

Locus- and domain-dependent control of DNA methylation at mouse B1 retrotransposons during male germ cell development.

In mammals, germ cells undergo striking dynamic changes in DNA methylation during their development. However, the dynamics and mode of methylation are poorly understood for short interspersed elements (SINEs) dispersed throughout the genome. We investigated the DNA methylation status of mouse B1 SINEs in male germ cells at different developmental stages. B1 elements showed a large locus-to-locus variation in methylation; loci close to RNA polymerase II promoters were hypomethylated, while most others were hypermethylated. Interestingly, a mutation that eliminates Piwi-interacting RNAs (piRNAs), which are involved in methylation of long interspersed elements (LINEs), did not affect the level of B1 methylation, implying a piRNA-independent mechanism. Methylation at B1 loci in SINE-poor genomic domains showed a higher dependency on the de novo DNA methyltransferase DNMT3A but not on DNMT3B, suggesting that DNMT3A plays a major role in methylation of these domains. We also found that many genes specifically expressed in the testis possess B1 elements in their promoters, suggesting the involvement of B1 methylation in transcriptional regulation. Taken altogether, our results not only reveal the dynamics and mode of SINE methylation but also suggest how the DNA methylation profile is created in the germline by a pair of DNA methyltransferases.

Endogenous siRNAs from naturally formed dsRNAs regulate transcripts in mouse oocytes.

RNA interference (RNAi) is a mechanism by which double-stranded RNAs (dsRNAs) suppress specific transcripts in a sequence-dependent manner. dsRNAs are processed by Dicer to 21-24-nucleotide small interfering RNAs (siRNAs) and then incorporated into the argonaute (Ago) proteins. Gene regulation by endogenous siRNAs has been observed only in organisms possessing RNA-dependent RNA polymerase (RdRP). In mammals, where no RdRP activity has been found, biogenesis and function of endogenous siRNAs remain largely unknown. Here we show, using mouse oocytes, that endogenous siRNAs are derived from naturally occurring dsRNAs and have roles in the regulation of gene expression. By means of deep sequencing, we identify a large number of both approximately 25-27-nucleotide Piwi-interacting RNAs (piRNAs) and approximately 21-nucleotide siRNAs corresponding to messenger RNAs or retrotransposons in growing oocytes. piRNAs are bound to Mili and have a role in the regulation of retrotransposons. siRNAs are exclusively mapped to retrotransposons or other genomic regions that produce transcripts capable of forming dsRNA structures. Inverted repeat structures, bidirectional transcription and antisense transcripts from various loci are sources of the dsRNAs. Some precursor transcripts of siRNAs are derived from expressed pseudogenes, indicating that one role of pseudogenes is to adjust the level of the founding source mRNA through RNAi. Loss of Dicer or Ago2 results in decreased levels of siRNAs and increased levels of retrotransposon and protein-coding transcripts complementary to the siRNAs. Thus, the RNAi pathway regulates both protein-coding transcripts and retrotransposons in mouse oocytes. Our results reveal a role for endogenous siRNAs in mammalian oocytes and show that organisms lacking RdRP activity can produce functional endogenous siRNAs from naturally occurring dsRNAs.

Tudor domain containing 7 (Tdrd7) is essential for dynamic ribonucleoprotein (RNP) remodeling of chromatoid bodies during spermatogenesis.

In the male germline in mammals, chromatoid bodies, a specialized assembly of cytoplasmic ribonucleoprotein (RNP), are structurally evident during meiosis and haploidgenesis, but their developmental origin and regulation remain elusive. The tudor domain containing proteins constitute a conserved class of chromatoid body components. We show that tudor domain containing 7 (Tdrd7), the deficiency of which causes male sterility and age-related cataract (as well as glaucoma), is essential for haploid spermatid development and defines, in concert with Tdrd6, key biogenesis processes of chromatoid bodies. Single and double knockouts of Tdrd7 and Tdrd6 demonstrated that these spermiogenic tudor genes orchestrate developmental programs for ordered remodeling of chromatoid bodies, including the initial establishment, subsequent RNP fusion with ubiquitous processing bodies/GW bodies and later structural maintenance. Tdrd7 suppresses LINE1 retrotransposons independently of piwi-interacting RNA (piRNA) biogenesis wherein Tdrd1 and Tdrd9 operate, indicating that distinct Tdrd pathways act against retrotransposons in the male germline. Tdrd6, in contrast, does not affect retrotransposons but functions at a later stage of spermiogenesis when chromatoid bodies exhibit aggresome-like properties. Our results delineate that chromatoid bodies assemble as an integrated compartment incorporating both germline and ubiquitous features as spermatogenesis proceeds and that the conserved tudor family genes act as master regulators of this unique RNP remodeling, which is genetically linked to the male germline integrity in mammals.

ADAD2 functions in spermiogenesis and piRNA biogenesis in mice.

BACKGROUND: Adenosine deaminase domain containing 2 (ADAD2) is a testis-specific protein composed of a double-stranded RNA binding domain and a non-catalytic adenosine deaminase domain. A recent study showed that ADAD2 is indispensable for the male reproduction in mice. However, the detailed functions of ADAD2 remain elusive. OBJECTIVES: This study aimed to investigate the cause of male sterility in Adad2 mutant mice and to understand the molecular functions of ADAD2. MATERIALS AND METHODS: Adad2 homozygous mutant mouse lines, Adad2-/- and Adad2Δ/Δ , were generated by CRISPR/Cas9. Western blotting and immunohistochemistry were used to reveal the expression and subcellular localization of ADAD2. Co-immunoprecipitation tandem mass spectrometry was employed to determine the ADAD2-interacting proteins in mouse testes. RNA-sequencing analyses were carried out to analyze the transcriptome and PIWI-interacting RNA (piRNA) populations in wildtype and Adad2 mutant testes. RESULTS: Adad2-/- and Adad2Δ/Δ mice exhibit male-specific sterility because of abnormal spermiogenesis. ADAD2 interacts with multiple RNA-binding proteins involved in piRNA biogenesis, including MILI, MIWI, RNF17, and YTHDC2. ADAD2 co-localizes and forms novel granules with RNF17 in spermatocytes. Ablation of ADAD2 impairs the formation of RNF17 granules, decreases the number of cluster-derived pachytene piRNAs, and increases expression of ping-pong-derived piRNAs. DISCUSSION AND CONCLUSION: In collaboration with RNF17 and other RNA-binding proteins in spermatocytes, ADAD2 directly or indirectly functions in piRNA biogenesis.

A novel class of small RNAs bind to MILI protein in mouse testes.

Small RNAs bound to Argonaute proteins recognize partially or fully complementary nucleic acid targets in diverse gene-silencing processes. A subgroup of the Argonaute proteins--known as the 'Piwi family'--is required for germ- and stem-cell development in invertebrates, and two Piwi members--MILI and MIWI--are essential for spermatogenesis in mouse. Here we describe a new class of small RNAs that bind to MILI in mouse male germ cells, where they accumulate at the onset of meiosis. The sequences of the over 1,000 identified unique molecules share a strong preference for a 5' uridine, but otherwise cannot be readily classified into sequence families. Genomic mapping of these small RNAs reveals a limited number of clusters, suggesting that these RNAs are processed from long primary transcripts. The small RNAs are 26-31 nucleotides (nt) in length--clearly distinct from the 21-23 nt of microRNAs (miRNAs) or short interfering RNAs (siRNAs)--and we refer to them as 'Piwi-interacting RNAs' or piRNAs. Orthologous human chromosomal regions also give rise to small RNAs with the characteristics of piRNAs, but the cloned sequences are distinct. The identification of this new class of small RNAs provides an important starting point to determine the molecular function of Piwi proteins in mammalian spermatogenesis.

MORC3, a novel MIWI2 association partner, as an epigenetic regulator of piRNA dependent transposon silencing in male germ cells.

The PIWI (P-element-induced wimpy testis)-interacting-RNA (piRNA) pathway plays a crucial role in the repression of TE (transposable element) expression via de novo DNA methylation in mouse embryonic male germ cells. Various proteins, including MIWI2 are involved in the process. TE silencing is ensured by piRNA-guided MIWI2 that recruits some effector proteins of the DNA methylation machinery to TE regions. However, the molecular mechanism underlying the methylation is complex and has not been fully elucidated. Here, we identified MORC3 as a novel associating partner of MIWI2 and also a nuclear effector of retrotransposon silencing via piRNA-dependent de novo DNA methylation in embryonic testis. Moreover, we show that MORC3 is important for transcription of piRNA precursors and subsequently affects piRNA production. Thus, we provide the first mechanistic insights into the role of this effector protein in the first stage of piRNA biogenesis in embryonic TE silencing mechanism.

Gtsf1/Cue110, a gene encoding a protein with two copies of a CHHC Zn-finger motif, is involved in spermatogenesis and retrotransposon suppression in murine testes.

We recently reported that the Gtsf1/Cue110 gene, a member of the evolutionarily conserved UPF0224 family, is expressed predominantly in male germ cells, and that the GTSF1/CUE110 protein is localized to the cytoplasm of these cells in the adult testis. Here, to analyze the roles of the Gtsf1/Cue110 gene in spermatogenesis, we produced Gtsf1/Cue110-null mice by gene targeting. The Gtsf1/Cue110-null mice grew normally and appeared healthy; however, the males were sterile due to massive apoptotic death of their germ cells after postnatal day 14. In contrast, the null females were fertile. Detailed analyses revealed that the Gtsf1/Cue110-null male meiocytes ceased meiotic progression before the zygotene stage. Thus, the Gtsf1/Cue110 gene is essential for spermatogenesis beyond the early meiotic phase. Furthermore, the loss of the Gtsf1/Cue110 gene caused increased transcription of the long interspersed nucleotide element (Line-1) and the intracisternal A-particle (IAP) retrotransposons, accompanied by demethylation of their promoter regions. These observations indicate that Gtsf1/Cue110 is required for spermatogenesis and involved in retrotransposon suppression in male germ cells.

Associations between PIWI proteins and TDRD1/MTR-1 are critical for integrated subcellular localization in murine male germ cells.

The mouse Piwi family proteins (MILI, MIWI and MIWI2) play pivotal roles in spermatogenesis through transcriptional and post-transcriptional gene regulation. To reveal the molecular functions of these proteins, we investigate the proteins that bind to MILI in adult mouse testes. We found that both MILI and MIWI bind to TDRD1/MTR-1, which is also an essential protein for spermatogenesis. Co-immunoprecipitation assays and subcellular localization of the proteins and mutants thereof revealed a complex formation involving MILI, MIWI and TDRD1/MTR-1. In addition, the subcellular localizations of MILI and TDRD1/MTR-1 were altered, and chromatoid body formation was impaired in the MIWI-null round spermatids. These data suggest that the formation of complexes between MILI, MIWI and TDRD1/MTR-1 is critical for the integrated subcellular localizations of these proteins, and is presumably essential for spermatogenesis.

Relationship between PIWIL4-Mediated H3K4me2 Demethylation and piRNA-Dependent DNA Methylation.

Retrotransposon genes are silenced by DNA methylation because of potential harm due to insertional mutagenesis. DNA methylation of retrotransposon genes is erased and re-established during male germ cell development. Both piRNA-dependent and piRNA-independent mechanisms are active during the re-establishment process, with the piRNA-independent mechanism occurring first. In this study, we analyzed the role of PIWIL4/MIWI2 in the modification of histone H3 and subsequent piRNA-dependent DNA methylation. Dimethylation at H3K4 is highly enriched at piRNA-dependent methylated regions and anti-correlated with de novo DNA methylation during the phase of piRNA-independent DNA methylation. In addition, PIWIL4, which binds the H3K4 demethylases KDM1A and KDM5B, is required for removing H3K4me2 marks. These data show that PIWIL4 plays important roles in histone modification and piRNA-dependent DNA methylation.