#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Gtsf1 is essential for proper female sex determination and transposon silencing in the silkworm, Bombyx mori


Authors: Kai Chen aff001;  Ye Yu aff001;  Dehong Yang aff001;  Xu Yang aff001;  Linmeng Tang aff001;  Yujia Liu aff001;  Xingyu Luo aff001;  James R. Walter aff003;  Zulian Liu aff001;  Jun Xu aff001;  Yongping Huang aff001
Authors place of work: Key Laboratory of Insect Developmental and Evolutionary Biology, Center for Excellence in Molecular Plant Sciences, Shanghai Institute of Plant Physiology and Ecology, Chinese Academy of Sciences, Shanghai, China aff001;  CAS Center for Excellence in Biotic Interactions, University of Chinese Academy of Sciences, Beijing, China aff002;  Department of Ecology and Evolutionary Biology, University of Kansas, NV, United States of America aff003
Published in the journal: Gtsf1 is essential for proper female sex determination and transposon silencing in the silkworm, Bombyx mori. PLoS Genet 16(11): e32767. doi:10.1371/journal.pgen.1009194
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1009194

Summary

Sex determination pathways are astoundingly diverse in insects. For instance, the silk moth Bombyx mori uniquely use various components of the piRNA pathway to produce the Fem signal for specification of the female fate. In this study, we identified BmGTSF1 as a novel piRNA factor which participates in B. mori sex determination. We found that BmGtsf1 has a distinct expression pattern compared to Drosophila and mouse. CRISPR/Cas9 induced mutation in BmGtsf1 resulted in partial sex reversal in genotypically female animals by shifting expression of the downstream targets BmMasc and Bmdsx to the male pattern. As levels of Fem piRNAs were substantially reduced in female mutants, we concluded that BmGtsf1 plays a critical role in the biogenesis of the feminizing signal. We also demonstrated that BmGTSF1 physically interacted with BmSIWI, a protein previously reported to be involved in female sex determination, indicating BmGTSF1 function as the cofactor of BmSIWI. BmGtsf1 mutation resulted in piRNA pathway dysregulation, including piRNA biogenesis defects and transposon derepression, suggesting BmGtsf1 is also a piRNA factor in the silkworm. Furthermore, we found that BmGtsf1 mutation leads to gametogenesis defects in both male and female. Our data suggested that BmGtsf1 is a new component involved in the sex determination pathway in B. mori.

Keywords:

Transposable elements – Drosophila melanogaster – Gonads – Ovaries – Immunoprecipitation – Sex determination – Immunostaining – Silkworms

Introduction

The mechanisms of sex determination are highly diverse in different species of insects [14]. In Drosophila melanogaster, the initial signal of sex is the X chromosome dosage present in the early zygote [58]. Individuals with two X chromosomes (XX) will adopt the female fate, while individuals with one X chromosome (XY) will develop to be male [812]. Different from XX/XY systems, the silkworm, B. mori possesses the WZ/ZZ sex determination system in which females are the heterogametic sex [13, 14]. A PIWI-interacting RNA (piRNA) precursor called Feminizer (Fem), which has more than 30 copies distributed in the sex determination region of the W chromosome, has been proposed to act as the female determining factor in the silkworm. piRNAs derived from Fem repress the expression of the Masculinizer gene (Masc) which located on the Z chromosome through piRNA pathway, a mechanism for transposon silencing. When Masc levels are low, the final effector in the pathway, the dsx gene, expresses female-specific variants which instruct female development [15, 16]. In contrast, when Fem is absent in ZZ individuals, Masc is expressed at high levels and activates expression of male-specific variants of dsx resulting in male development. Except for Fem and Masc, there are other two sex determination factors, B. mori insulin-like growth factor II mRNA-binding protein (BmIMP) and P-element Somatic Inhibitor (BmPSI). BmPSI and male-specific BmIMP (BmIMPM) are involved in male-specific Bmdsx splicing through inhibiting the female-specific Bmdsx splicing in males [17, 47].The function of Fem requires assistance of piRNA pathway gene like BmSiwi, deletion of which can cause partial sexual reversal of female [18]. However, except for BmSIWI and its cofactor BmASH2, no any other piRNA pathway factor was shown to be involved in sex determination pathway in the silkworm.

In the past few years, many novel piRNA factors have been uncovered. Several studies have shown that an evolutionarily conserved protein called Gametocyte-specific factor1 (GTSF1) is involved in the piRNA pathway [1922]. GTSF1 is a protein which consists of two U11-48K-like CHHC-type Zinc finger motifs at N-terminal and a C-terminal end without discernable domains [20, 21]. In Drosophila and mouse, GTSF1 is preferentially expressed in germ cells. In Drosophila, GTSF1 (also known as Asterix) is required for female fertility and interacts with PIWI via its C-terminal end. DmGTSF1 is essential for transposon silencing but not for piRNA biogenesis. In absence of GTSF1, localization of PIWI is unaffected, and PIWI can still bind piRNAs, whereas introduction of H3K9me3 marks into transposon loci is impacted [20, 21]. In mouse, GTSF1 is necessary for spermatogenesis and male fertility, and associates with PIWI proteins via its central region instead of its C-terminal tail. Though murine GTSF1 is indispensable for secondary piRNA genesis, it is not needed for primary piRNA biogenesis [22].

Silencing transposable elements (TEs) is the primary function of PIWI-interacting RNA (piRNA) pathway, despite it also plays vital role in the silkworm sex determination. TEs pose a potential threat to genome stability. In particular, their presence in germ cells can cause damage to the germline genome and may result in sterility [2325]. In the germline, piRNA pathway functions as a defense system in silencing TEs, and maintains genomic integrity [2628]. piRNAs are a class of small non-coding RNA with 24–32 nucleotides deeply conserved from yeast to human. piRNAs bind to PIWI proteins and form piRNA induced silencing complex (pi-RISC) to repress TEs expression at the transcriptional and post-transcriptional level [2931]. The biosynthesis of piRNAs consists of two pathways: primary processing pathway and secondary processing pathway (also referred to as the ping-pong cycle) [25, 27]. In Drosophila, three PIWI proteins serve as the key components of piRNA pathway: Piwi, Argonaute (Ago3) and Aubergine (Aub). PIWI and Aub bind to primary piRNAs, while Ago3 binds to secondary piRNAs [3234]. In B. mori, there are only two PIWI proteins: BmSIWI and BmAgo3 [35, 36], indicating discrepant mechanisms of piRNA pathway among species. Beyond PIWI proteins, many other factors are also indispensable for piRNA pathway, and mutational disruption of these piRNA factors can result in spermatogenesis or oogenesis defects in mouse and Drosophila respectively [26, 31, 3743]. However, the research upon piRNA pathway in the silkworm were concentrated on BmN4 cell line, thus, little was known about the function of piRNA pathway in the silkworm gametogenesis.

In this study we examined the hitherto unknown functions of the silkworm ortholog of Gtsf1 (BmGtsf1). Of particular interest was our finding that BmGTSF1 plays a critical role in female sex determination. BmGtsf1 mutations resulted in partial sex reversal of female animals, while sexual development of males remain unaffected. Furthermore, BmGTSF1 appears necessary for TE silencing during gametogenesis in both male and female, given that the lack of BmGTSF1 causes sterility most likely due to genomic instability. Subsequently, we found that BmGTSF1 interacted with BmSIWI (Silkworm PIWI protein), and that disruption of BmGTSF1 function led to transposon derepression and dysregulation of piRNA biogenesis. Taken together, our work reaveals that BmGTSF1 is a novel factor in female sex determination and piRNA pathway.

Results

Expression pattern of BmGtsf1 and mutants construction

Previous studies in Drosophila and mouse reported that Gtsf1 is predominantly expressed in the germ line [2022]. In the silkworm we found that BmGtsf1 was ubiquitously expressed at all stages tested with a peak level of expression in gonads (Fig 1A and 1B). High levels of BmGtsf1 mRNA present in gonads from early larval stage to pupal stage suggested an important role of BmGtsf1 in gonadal development. In addition, the expression level of BmGtsf1 in ovaries was higher than that in testes, revealing BmGtsf1 may have more important function in female (Fig 1C). Moreover, we also found that the mRNA level of BmGtsf1 in female was higher than that in male at early embryonic stage (Fig 1D).

Fig. 1. Expression patterns of BmGTSF1.
Expression patterns of BmGTSF1.
(A) qRT-PCR analysis showed relative mRNA levels of BmGtsf1 in eight tissues at three different stages: third day of fifth larval stage (L5D3), wandering stage (W), and day 2 of the prepupal stage (PP2). The mRNA levels were normalized to levels of B. mori ribosomal protein 49 (Bmrp49). Error bars represent ± SD. (B) Western blot analysis of BmGTSF1 in different tissues in wandering stage. (C) The mRNA levels of BmGtsf1 in testis and ovary from day 1 of fifth larval stage (L5D1) to day 6 of pupal stage (P6). Tissues tested were epidermis (EPI), midgut (MG), fat body (FB), silk gland (SG), malpighian tubule (MT), testis (TE), ovary (OV), and wing disc (WD). (D) Relative mRNA level of BmGtsf1 from 9hpo to 30hpo (hours after oviposition).

To explore the physiological function of BmGtsf1, we utilized a previously described binary transgenic CRISPR/Cas9 system [44, 45] to establish somatic mutant lines (ΔBmGtsf1). Briefly, we constructed two transgenic silkworm lines expressing Cas9 and the sgRNA, respectively. The two lines were crossed to obtain somatic BmGtsf1 mutants. We were able to generate lesions in the BmGtsf1 locus in ΔBmGtsf1 animals (S2A Fig and S2B Fig), and the depletion efficiency of RNA and protein levels was further confirmed by qRT-PCR and western blot (S2C Fig and S2D Fig).

BmGtsf1 is involved in female sex determination

Phenotypic analysis of ΔBmGtsf1 animals revealed a critical role of BmGtsf1 in female sex determination. Female mutants (the genotype is WZ) displayed intersexual external genitalia with distinct male-specific structures (Fig 2A, red and green arrow). In addition,female mutant adults developed a male-specific eighth abdominal segment (Fig 2B). Male mutants (the genotype is ZZ) developed normal male genitalia and abdominal segmentation (Fig 2B). In B. mori, female pupae are on average larger and heavier than male pupae. However, both body size and weight of female mutant pupae were reduced and closer to male, while those of male mutant were unaffected (S3 Fig). And we confirmed that the phenotypes were identical among the different individuals.

Fig. 2. Partial sexual reversal of female BmGtsf1 mutants.
Partial sexual reversal of female <i>BmGtsf1</i> mutants.
(A) external genitalia of WT and ΔBmGtsf1 adults. Key to arrows: white, ventral plates; blue, ventral chitin plate; purple, genital papillae; red, clasper (the male specific structure); green, penis (the male specific structure). The red and green arrow showed male specific structures existed in female mutant. (B) Abdominal segment of WT and BmGtsf1 mutant adults. BmGtsf1 female mutant developed male specific eight abdominal segments.

Subsequently, we examined expression levels of BmMasc, the upstream regulator of Bmdsx [17, 45]. Levels of BmMasc transcripts were upregulated by 3.55-fold in the female mutants, but remained unaffected in male mutants (Fig 3A). Disruption of female development of female mutants seems to result from expression of male variants of Bmdsx, the conserved downstream effector gene in the B. mori sex determination pathway [45], as we detected male-specific splicing of Bmdsx in female ΔBmGtsf1 animals (Fig 3B). In addition, we also tested expression of several putative target genes of Bmdsx. The mRNA levels of three male-specific expressed genes in olfactory system, pheromone binding protein 1 (BmPBP1), olfactory receptors 1 (BmOR1), BmOR3, were significantly increased, while expression of two female olfactory system genes, BmOR19, BmOR30, and the yolk gene, vitellogenin (BmVg), were significantly reduced in female mutants (Fig 3C–3H). In contrast, transcripts levels of all these genes were not affected in male mutants (S4 Fig). These results suggested that BmGtsf1 plays a critical role in the B. mori female sex determination pathway, particularly through modulating the expression level of BmMasc, although more evidences were needed to confirm the direct regulation mechanism. Expression of the auxiliary sex determination factors, BmIMP and BmPSI [17, 45], appeared not to be affected (Fig 4A–4D). As levels of Fem piRNA and Masc piRNA were substantially downregulated in female mutants (Fig 4E and 4F), we concluded that BmGtsf1 mutation primarily affects the sex determination pathway at the level of Fem piRNA biogenesis. Through RT-PCR analysis, we found that BmGtsf1 is involved in initial sex determination pathway at early embryonic stage. In the WT, only faint band for male-specific Bmdsx was observed in female. While in the female mutants, male-specific band was enhanced. In addition, the abundance of Fem precursor was down-regulated in female mutant (Fig 4G).

Fig. 3. Alternative splicing pattern of Bmdsx and qRT-PCR analysis of BmMasc and putative target genes of Bmdsx.
Alternative splicing pattern of <i>Bmdsx</i> and qRT-PCR analysis of <i>BmMasc</i> and putative target genes of <i>Bmdsx</i>.
(A) Relative expression of BmMasc in gonads of WT and ΔBmGtsf1 at wandering stage. (B) Splicing patterns of Bmdsx in gonads of WT and BmGtsf1 mutants at wandering stage detected by RT-PCR. (C-H) Relative expression of BmPBP1 (C), BmOR1(D), BmOR3 (E), BmOR19 (F), BmOR30 (G) and BmVg (H) in BmGtsf1 female mutants. RNA extracted from antenna of adults was used for qRT-PCR analyses in Fig 3C–Fig 3G, while the expression of BmVG was detected in fat body at wandering stage. Error bars represent ±SD. The asterisks represent statistically significant differences with p < 0.001.
Fig. 4. BmGtsf1 regulate female sex determination via Fem piRNA.
<i>BmGtsf1</i> regulate female sex determination via <i>Fem</i> piRNA.
(A-C) Relative mRNA levels of BmPSI (A), BmIMP (B) and BmImpM (C, male-specific BmIMP) in gonads of WT and mutants at wandering stage. (D) RT-PCR analysis of BmImpM in gonads of WT and mutants at wandering stage. (E and F) Relative abundance of Fem piRNA (E) and Masc piRNA (F) in ovary and fat body of female mutant at wandering stage, the small RNA U6 was used as the internal reference. (G) Splicing patterns of Bmdsx and expression of Fem precursor in WT and BmGtsf1 mutant embryoes at 24hpo detected by RT-PCR. The asterisks represent statistically significant differences with p < 0.001. Error bars represent ±SD.

To confirm whether the BmGtsf1 mutation impacts the function of Fem piRNA, we constructed four plasmids: PIZT-Cas9-EGFPsgRNA (control), PIZT-Cas9-U6-Fem (to overexpress Fem piRNA, Fem-OE), PIZT-Cas9-Gtsf1sgRNAs (to knock out Gtsf1, Gtsf1-KO), PIZT-Cas9-Gtsf1sgRNAs-U6-Fem (to knock out Gtsf1 and overexpress Fem piRNA, Gtsf1-ko-Fem-OE) to perform a rescue assay in BmN cells (S5A Fig). BmN cells transfects with PIZT-Cas9-EGFPsgRNA was served as control. In accordance with the results in vivo, BmGtsf1 depletion resulted in Fem piRNA down-regualtion and BmMasc up-regulation in BmN cells. We found that the abundance of Fem piRNA can be rescued through U6 promotor overexpressing. However, the up-regulation of BmMasc caused by BmGtsf1 depletion cannot be rescued, although the level of Fem piRNA was much higher than control (S5B Fig and S5C Fig). These data suggested that BmGtsf1 is also required for the function of Fem piRNA. We next overexpressed BmGtsf1 in vivo (BmGtsf1-OE) using IE1 promoter (S6A Fig), and the results showed that neither the expression of BmMasc nor the splicing of Bmdsx was altered (S6B Fig and S6C Fig). No detectable phenotype was observed in BmGtsf1 overexpressed animals, and all the adults were fertile, indicating BmGtsf1 does not has dosage effect (S6D Fig). In summary, from our data we concluded that BmGtsf1 controls female sex determination through assisting the production and function of Fem piRNAs, which are needed to effectively repress the male determining gene BmMasc.

BmGtsf1 mutation results in gametogenesis defects in both male and female

To investigate the effects of BmGtsf1 mutation on gonadal development, we performed a morphological examination. The testes of the adult male mutants were highly atrophied (Fig 5A, top row), nevertheless the weights of male mutants were unchanged (S3B Fig). Wild-type females contain two ovaries, each of which consists of four ovarioles filled with eggs at various stages of development. In contrast, the ovarioles of female mutants were severely degraded, neither obvious ovariole nor eggs were observed (Fig 5A, bottom row). To further examine the internal structures of the female and male gonads, we stained paraffin-embedded tissue sections with hematoxylin and eosin at the wandering stage. This examination also revealed that the testes development of male mutants were affected. They were much smaller in overall size and the sperm count was greatly reduced (Fig 5B, top row). In B. mori, there are dimorphic sperm: eupyrene sperm bundles and apyrene sperm bundles [46]. In the WT adult males, both kinds of sperm bundles were matured, by comparison, spermatogenesis was arrested in male mutant, and we can only found round spermatogonium (Fig 5C). Similarly, the ovary of mutant was also atrophied, and no clear ovariole was observed (Fig 5B, bottom row). As a result of defective gametogenesis, both male and female mutants were sterile (Fig 5D). Taken together, these data demonstrated that BmGtsf1 deficiency will lead to gametogenesis defects in both male and female silkworms.

Fig. 5. BmGtsf1 mutation results in gametogenesis defects.
<i>BmGtsf1</i> mutation results in gametogenesis defects.
(A) Morphological analysis of gonads in the WT and mutant adults. (B) Internal structure of gonads in the WT and mutant animals at wandering stage using hematoxylin and eosin staining. Scale bar: 200 μm. The green arrow indicates ovariole. (C) Sperm bundles from WT and mutant male adults. C1, eupyrene sperm bundles; C2, apyrene sperm bundles; C3, abnormal sperm bundles in BmGtsf1 mutant. (D) Fertility of WT and BmGtsf1 mutants, n = 20.

BmGTSF1 interacts with BmSIWI

Given the fact that GTSF1 interacts with PIWI in both mouse and Drosophila [2022], we next examined the relationship of the two proteins in the silkworm. Firstly, we detected their localization in the gonads. For clear analysis of immunostaining results, we used a primary antibody recognizing BmVasa, which is a germline specific marker in insects, to differentiate germline cells and somatic cells. In ovary, BmGTSF1 and BmSIWI were present in both germline and surrounding somatic cells (Fig 6A). In testis, the distribution of BmVasa, BmGTSF1 and BmSIWI was almost overlapping, and restricted to spermatogonium (S7 Fig). We also examined the cellular localization of BmSIWI and BmGTSF1 in BmN and Bm12 cells, which were derived from B. mori ovary. The immunostaining results showed that BmSIWI located to the cytoplasm, while BmGTSF1 localized in both nucleus and cytoplasm (Fig 6B), indicating these two proteins have different intracellular localization. To test for a direct interaction between BmGTSF1 and BmSIWI, we fused BmGTSF1 (BmGTSF1-Flag) and BmSiwi (BmSiwi-Flag) C-terminally with a Flag tag and used a monoclonal anti-Flag antibody for immunoprecipitation in BmN cells. DsRed fused C-terminally with Flag tag (DsRed-Flag) was served as control. The western blot followed immunoprecipitation confirmed that the BmGTSF1-Flag coprecipitated with endogenous BmSIWI (Fig 6C), more importantly, the reciprocal immunoprecipitation experiment showed that BmSIWI-Flag also coprecipitated with endogenous BmGTSF1 (Fig 6D). In conclusion, BmGTSF1 and BmSIWI showed the same localization in the gonads of B.mori, and they physically interacted.

Fig. 6. BmGTSF1 interacts with BmSIWI.
BmGTSF1 interacts with BmSIWI.
(A) Localization of BmVasa, BmSIWI and BmGTSF1 in the ovaries of B. mori. A FITC-conjugated secondary antibody was used for fluorescence detection and Hoechst staining (blue) showed the locations of the nuclei. Scale bar, 50μm. The white and red arrow indicates germ line cells and somatic cells respectively. (B) Intracellular localization of BmSIWI and BmGTSF1 in BmN cells and Bm12 cells. (C) Immunoprecipitation followed by western blot showed that endogenous BmSIWI immunoprecipitated with BmGTSF1-Flag rather than DsRed-Flag. (D) Immunoprecipitation followed by western blot showed that endogenous BmGTSF1 immunoprecipitated with BmSIWI-Flag rather than DsRed-Flag.

BmGtsf1 is essential for piRNA induced transposon silencing

Given the well-known functions of Gtsf1 in Drosophila and mouse, we next tested whether BmGtsf1 is also required for piRNA pathway in B.mori. To this end, we analyzed piRNA abundance and transposons expression level in the gonads of WT and mutant by transcriptome analysis and qRT-PCR. In ovaries of ΔBmGtsf1 animals, abundance of 24–27nt long piRNAs was greatly decreased, while the amount of 28–30nt long piRNAs was slightly increased compared to wild-type ovaries (Fig 7A). Moreover, the abundance of 24–30nt long piRNAs was remarkably reduced in testis suggesting that BmGtsf1 is required for piRNA maturation in the gonads (Fig 7B).

Fig. 7. BmGtsf1 mutation mediated piRNA pathway dysregulation.
<i>BmGtsf1</i> mutation mediated piRNA pathway dysregulation.
(A and B) Size distribution of piRNAs in ovary (A) and testis (B) of WT and BmGtsf1 mutant. (C and D) Fold change of transposons in ovary (C) and testis (D) of BmGtsf1 mutant compared to WT. (E and F) Relative abundance of piRNAs correspond to detected transposons.

In addition, we found that lack of BmGtsf1 activity resulted in derepression of numerous transposons in the gonads of mutant animals (S8A Fig and S8B Fig). More specifically, we examined mRNA levels of seven well annotated transposons [48, 49] by qRT-PCR. The expression of all of these transposons were greatly upregulated (Fig 7C and 7D). In line with the observation, abundance of piRNAs correspond to these transposons were declined (Fig 7E and 7F), except for the piRNAs correspond to Pao and Bm1-Neet (Bm1-related_Neet). In addition, levels of Fem piRNA was severely reduced in ovaries of ΔBmGtsf1 animals (Fig 4E). Fem piRNA depletion was also observed in BmSiwi mutant animals [18], supporting the notion that BmGTSF1 interacts with BmSIWI. Altogether, our data demonstrated that BmGtsf1 is involved in the piRNA pathway. These results also revealed that BmGtsf1 regulate female sex determination via piRNA pathway.

Discussion

The functions of GTSF1 have thus far been examined in three species, Drosophila, mouse and C. elegans [2022, 50]. Despite the fact that GTSF1 is conserved, its functions appear to have diversified in these three species. In both Drosophila and mouse, GTSF1 is necessary for piRNA pathway, however, the mechanisms appear to be quite different. In Drosophila, DmGTSF1 localizes to the nucleus, and is required for Piwi-piRISC-mediated transcriptional transposon silencing but not for piRNA biogenesis [20, 21]. While in mouse, GTSF1 localizes to both cytoplasm and nucleus, and is essential for secondary but not for primary piRNA biogenesis [22]. In C. elegans, GTSF1 does not participate in the piRNA pathway, but instead is involved in the assembly of ERI (enhanced RNAi) complex [50]. In B. mori, as BmGtsf1 mutations caused derepression of TEs and dysregulation of piRNAs, we concluded that BmGTSF1 is not only necessary for transposon silencing but also for piRNA biogenesis. In addition, our co-immunoprecipitation experiments suggested that BmGTSF1 interacts with BmSIWI, such an interaction was also observed in Drosophila and mouse [2022], suggesting a relatively conserved function of GTSF1 in piRNA pathway. However, although these proteins were interacted, their cellular localization was different. We speculated that BmGTSF1 may have another role in nucleus, for example, it might participate in transporting piRNA precursors. There are two facts in accordance with our hypothesis. The one is that BmGTSF1 contains two CHHC-type Zinc finger motifs (S2A Fig), the predicted RNA binding domain. The other is that we observed condensation of BmGTSF1 in nuclues (Fig 6B). The roles of GTSF1 in B. mori and Drosophila seemed quite different, we supposed the different mechanism is related to the fact that transcriptional transposon silencing is absent in the silkworm [51]. In Drosophila and mouse, Gtsf1 is predominantly expressed in germ cells. While in the silkworm, BmGtsf1 was expressed in practically all tissues, consistent with the fact that BmGTSF1 is not only involved in the germline piRNA pathway but also in somatic sex determination.

piRNAs derived from Fem precursor serve as the primary signal for B. mori female sex determination by repressing the male determining gene Masc. Nevertheless, not all of the genes involved in processing and maturation of piRNAs participate in Masc repression. To date, only the core components of the piRNA pathway, BmSiwi, and its cofactor BmAsh2, have been shown to be involved in B. mori sex determination. In contrast, the factors BmAgo3 and BmMael appear not to be required for sex determination [18, 44]. This present study includes another piRNA pathway factor BmGtsf1 in the regulation of female sex determination in B. mori. Compromising the activity of BmGtsf1 results in partial sex reversal of WZ individuals, while no observable effects on sex determination were detected in ZZ male individuals. In accordance with a role in piRNA mediated repression of BmMasc, expression levels of BmMasc were significantly upregulated in WZ mutant individuals. As levels of two auxiliary sex determination factors IMP and PSI were unaffected, this failure to repress BmMasc can be likely attributed to a lack of Fem piRNA activity. Consistent with this interpretation, the abundance of Fem piRNA was found to be significantly low in WZ mutant individuals. Concomitantly, levels of Masc piRNAs derived from BmMasc mRNA through the ping-pong cycle were also drastically reduced. Furthermore, overexpression of Fem piRNAs in BmGtsf1 depletion BmN cells cannot restore repression of BmMasc. Thus, we concluded that BmGtsf1 is involved in female sex determination by assisting the proper processing of functional Fem piRNAs to repress BmMasc. However, we proposed that BmGtsf1 might be an auxiliary factor for female sex determination, since BmGtsf1 mutation caused partial sex reversal of female animals.

piRNA pathway was previously reported to act primarily as a silencing system which prevents TE expression in animal gonads [26, 31, 21, 52]. However, in the silkworm, piRNA pathway is not only required for TE silencing but, as shown in this study, is also involved sex determination. Hence, there seemed to be a somatic piRNA pathway which specifically works in the silkworm sex determination pathway. Consistent with our speculation, we detected that expression of several transposons were barely changed in mutant fat body (S8C Fig and S8D Fig), by contrast, the Fem and Masc piRNA level were significantly downregulated (Fig 4E and 4F). We assumed that BmSiwi was the core component of this somatic piRNA pathway, because BmAsh2 and BmGTSF1, the only two known piRNA factors required for B. mori sex determination, all interact with BmSIWI. Depletion of BmGTSF1; a factor interacting with BmSIWI, showed similar phenotypes as BmSiwi mutants does [18]. Nevertheless, the sex reversal phenotypes were more distinct in BmGtsf1 mutants. We speculated that the more distinct phenotypes were attributed to the extra function of BmGTSF1 in nucleus.

As the piRNA pathway is responsible for genome stability in animal gonads, depletion of piRNA factors often results in sterility attributed to DNA damage caused by TEs derepression [5254]. GTSF1 mutations lead to severely deformed ovaries and female sterility in Drosophila, while the testes appear unaffected [20, 21]. Mouse GTSF1 mutation resulted in much smaller testis and spermatogenesis defects, by contrast, no abnormalities were detected in female [22]. Namely, GTSF1 mutation specifically affects gonadogenesis of the single sex in Drosophila and mouse. Different from the sex-specific role in mouse and Drosophila germ line, we found that BmGtsf1 mutation would result in gametogenesis and gonadogenesis defects in both males and females. We showed that, similar to BmSIWI, BmGTSF1 is required for proper differentiation and maturation of eggs. We assumed that ovarian phenotypes can be attributed to dysregulation of both piRNA pathway and the sex determination pathway. In line with this, the expression of BmVG, a gene which predominantly expressed in female fat body, was drastically down-regulated in the female mutant individuals. An interesting finding was that BmGtsf1 is also required for spermatogenesis. The phenotype appeared to be caused by TEs derepression only, since BmGtsf1 mutation had few effects on male development and sex determination. We concluded that BmGtsf1 is essential for B. mori TE silencing during gametogenesis, although more direct evidence is needed to prove a causal relationship between TEs derepression and defects in gametogenesis.

In conclusion, our data provided a first functional analysis of BmGTSF1, the Drosophila ortholog of GTSF1 in the silkworm. We showed that BmGTSF1 is a new addition to the piRNA processing factors required not only in canonical TE silencing but also female sex determination pathway in B. mori (Fig 8). Although the piRNA pathway is involved in both TE silencing and female sex determination in the silkworm, not all of the factors regulating piRNA activity participate in both pathways. For future studies, it will be necessary to identify new piRNA factors and test whether or not they are also required in B. mori sex determination. Such studies will contribute to a better understanding and distinction of these two vitally important roles of the piRNA pathway in silkworm development.

Fig. 8. Proposed model for functions of BmGTSF1 in B. mori.
Proposed model for functions of BmGTSF1 in <i>B</i>. <i>mori</i>.
Although BmGTSF1 interacts with BmSiwi, the function of BmGTSF1 is different in male and female. In female, BmGTSF1 is reqiured not only in canonical TE silencing but also sex determination, both of which are vital for oogenesis. It was reported that another factor BmAsh2, which interactes with BmSiwi, is also involved in sex determination, although whether BmAsh2 interactes with BmGTSF1 is unkonwn. While in male, BmGTSF1 is only essential for transposon silencing, which ensures proper spermatogenesis.

Materials and methods

Silkworm strain and cell line

The multivoltine, nondiapausing silkworm strain, Nistari, was used in this study. The larvae were reared on fresh mulberry leaves under standard conditions [55]. The cell line BmN and Bm12, which are derived from silkworm ovary, were used for transfection and immunofluorescent staining. The cell line was cultured at 27°C in TC100 insect medium with 10% fetal calf serum (FBS) [56].

Information of BmGtsf1 gene

BmGtsf1 is located on chromosome11 according to SilkBase (http://silkbase.ab.a.u-tokyo.ac.jp/cgi-bin/index.cgi), and the gene ID of BmGtsf1 in NCBI (https://www.ncbi.nlm.nih.gov/gene/101745848) is 101745848.

Analysis of mutants genotype

The genomic DNA of the mutants was using standard SDS lysis-phenol treatment, incubated with proteinase K, treated with RNase treatment, and purified. 150ng of the genomic DNA was used as template, and the primer pair, Fem-seq-F and Fem-seq-R (S1 Table), which flank Fem precursor was used (S1 Fig). To confirm the genotype of early embryo, the genome was extracted according to the manufacturer's instructions for TRIzol reagent.

RNA isolation and cDNA synthesis

Total RNA was isolated from the tissues of three individual mutants or WT animals using the TRIzol reagent according to the manufacturer's instructions, followed by treatment with DNase I and purification. The quality of total RNA was confirmed by agarose gel electrophoresis. An aliquot of 1 μg of the total RNA was used to synthesize cDNA using ReverAid First Strand cDNA Synthesis Kit. For piRNA detecting, the total RNA was firstly incubated with piRNA specific stem-loop primer at 65°C for 5min, followed by normal cDNA synthesis steps according to the manufacturer's instructions.

Quantitative real-time PCR (qRT-PCR)

qRT-PCR analyses were performed using SYBR Green Real-time PCR Master Mix to analyze mRNA levels of selected genes and piRNAs abundance. mRNA levels were normalized to levels of B. mori ribosomal protein 49 (Bmrp49). The relative expression levels of piRNAs were measured by using the stem-loop method [57, 58], and the small RNA U6 was used as the internal reference. Three independent biological replicates of each qRT-PCR analysis were performed. All the primers used for qRT-PCR were listed in S1Table.

Western blot analysis

For western blot analysis, the following primary antibodies were used: polyclonal rabbit anti-BmGTSF1 (1:5000), polyclonal rabbit anti-BmSIWI (1:2000), monoclonal mouse anti-β-actin (1:5000), monoclonal mouse anti-Flag (1:5000). The peptides corresponding to the full length of BmGTSF1 was synthesized and used for generating polyclonal rabbit anti-BmGTSF1 primary antibody. For the secondary antibody, horseradish peroxidase-conjugated goat anti-rabbit and anti-mouse IgG (1:5000) were used. The ECL Plus Western blotting Detection Kit was used to detect the protein signal.

Paraffin sectioning and hematoxylin-eosin staining

Gonads of mutants and WT animals were dissected from the wandering stage larvae and prefixed with Qurnah's fixative (anhydrous ethanol: acetic acid: chloroform, 6:1:3 (v/v/v)) for 24h. Samples were dehydrated using anhydrous ethanol for three times (1h per time), followed by clearing for three times (10min per time) using xylene. Tissues were embedded in paraffin overnight, and were cut into cross-sections (5 μm) with a Leica RM2235 microtome. The sections were rehydrated and then stained using a mixture of hematoxylin and eosin to visualize morphology. The stained sections were photographed by using an Olympus BX53 microscope.

Gonads immunostaining

The gonads paraffin-embedded sections were firstly rehydrated by steps as follow: xylene for twice, 10min per time, 95% ethanol for twice, 5min per time, 80% ethanol for 5min, 70% ethanol for 5min, H2O for twice, 3min per time, followed by PBS (phosphate buffered saline) for 5min. Sections were treated with 0.1% trisodium citrate containing 0.1% Triton X-100 for 10min at room temperature for antigen retrieval. The samples were washed using PBS for twice, followed by blocking with 1% BSA (bovine serum albumin) for 1 hour at room temperature. The sections were then incubated with primary antibodies for 48 hours at 4°C. Primary antibodies used for immunostaining were as follow: polyclonal rabbit anti-BmGTSF1 (1:500), polyclonal rabbit anti-BmSIWI (1:200), polyclonal rabbit anti-BmVasa (1:200). The samples were washed with PBS for three times and incubated with FITC-AffiniPure Goat Anti-Rabbit IgG (H+L) (1:200) secondary antibody for 2 hours at room temperature. The samples were washed for once and treated with Hoechst33258 for 10min at room temperature. The samples were washed for four times and analyzed with an Olympus BX53 microscope.

Immunostaining of BmN and Bm12 cells

BmN and Bm12 cells were seed on cover glass put in culture dish. 12 hours after cell passaging, cells were fixed (Immunol Staining Fix Solution) for 1 hour at room temperature. The samples were washed twice for 5min and then permeabilized with 0.1% Triton X-100 in 0.1% trisodium citrate for15 minutes. After washing with PBS for 5min, the samples were blocked using 1% BSA for 2 hours at room temperature. The cells were incubated with primary antibodies for 48 hours at 4°C. Following three washing steps, cells were incubated with secondary antibody for 2 hours at room temperature. The antibodies used in cell lines staining were the same as gonads immunostaining. The samples were washed for once and treated with Hoechst33258 for 10min at room temperature. After washing for 4 times for 10min, the results were analyzed using Olympus Fluoview FV1000 Confocal Microscope.

Immunoprecipitation assay

Coding sequences of BmSIWI, BmGTSF1 and DsRed tagged with 3×FLAG were cloned into pIZT/V5-His plasmid. The plasmids were transfected into the BmN cell line respectively using Effectene Transfection Reagent according to the manufacturer's instructions. Crude proteins from one 60mm dish cells were collected using lysis buffer (50 mM Tris HCl, pH 7.4, 150 mM NaCl, 1 mM EDTA, 1% TRITON X-100) 72 hours after transfection, and then used for immunoprecipitation. Immunoprecipitation was performed using ANTI-FLAG M2 Affinity Gel. Immunoprecipitates were washed ten times with TBS (Tris Buffered Saline, 50 mM Tris HCl,150 mM NaCl, pH 7.4) at 4°C, 5min per time. The immunoprecipitation results were detected using western blot.

CRISPR/Cas9 system mediated mutant construction

A binary transgenic CRISPR/Cas9 system was applied to construct BmGtsf1 mutant. The plasmids construction and the silkworm germline transformation were described previously [44, 45]. Briefly, two silkworm lines were constructed, one line expressed Cas9 nuclease under the control of the B. mori nanos (nos) promoter, the other line expressed two guide RNAs activated by U6 promoter [44, 45]. The mixture of transformation plasmids and helper plasmids were microinjected into preblastoderm G0 embryos for germline transformation.

Mutagenesis analysis

Genomic DNA was extracted from the somatic mutants using standard SDS lysis-phenol treatment, incubated with proteinase K, treated with RNase treatment, and purified. To perform gene-specific PCR amplification, 100 ng of the genomic DNA was used as template, and the primer pair, BmGtsf1-seq-F and BmGtsf1-seq-R (S1 Table), which flank both targets was used. The PCR products were extracted and cloned into the pJET-1.2 vector, and 20 clones from 5 individual were sequenced. The results of sequencing are listed in S2 Table.

RNA-seq analysis

Total RNA from the testes and ovaries of wandering stage animals was extracted from three individual animals of ΔBmGtsf1 and WT, and mixed together. The total RNA was used for both mRNA sequencing and small RNA sequencing analysis. For mRNA sequencing, the total RNA was firstly enriched and then fragmented, and used for cDNA synthesis and library construction. The library was sequenced using BGISEQ-500 technology, raw data was qualified, filtered, and mapped (SOAPnuke, HISAT, Bowtie2) to the reference silkworm genome database (http://sgp.dna.affrc.go.jp/KAIKObase/), mRNA abundance was measured by fragment per kilobase of exon per million fragments mapped (FPKM). The data was then mapped to 1,811 silkworm transposons [59, 60] and analyzed the transposon expression. For small RNA sequencing, RNAs ranging in size between 18 and 32 nucleotides (nt) were gel-purified and used for library construction and sequencing (BGISEQ-500 technology). The generated reads were filtered and mapped to the silkworm genome (http://sgp.dna.affrc.go.jp/KAIKObase/). The mapped reads from 24 to 30nt in length were further mapped to 1811 silkworm transposons [59, 60] allowing no mismatches. These perfect mapping reads were defined as piRNAs. The relative abundance of piRNAs were measured by reads per million (RPM) by normalizing to the total number of perfect genome mapped sRNAs in each library [61].

Fluorescent staining of sperm bundles

Testes of mutant and WT animals were dissected from adult animals. Next, the testes were avulsed and the contents were fixed for 1 hour at room temperature. Samples were then washed three times using PBS, followed by staining with TRITC Phalloidin for 1 hour. Nuclei were stained with Hoechst33258 for 10 min at room temperature. Samples were washed three times with PBS, dropped onto a microslide, and immediately analyzed by a fluorescence microscopy (Olympus, BX53).

Statistics

Experimental data were analyzed with the Student's t-test. The asterisks indicate different p value, *, p<0.05; **, p<0.01; ***, p<0.001.

Supporting information

S1 Table [xlsx]
Primers used in this work.

S2 Table [xlsx]
Summary of the mutagenesis analysis.

S1 Fig [tif]
Analysis of mutants’ genotype.

S2 Fig [a]
CRISPR/Csa9 system mediated mutant construction.

S3 Fig [a]
Body size of WT and mutant pupa.

S4 Fig [a]
qRT-PCR analysis of male or female specific genes in male mutants.

S5 Fig [a]
is essential for the function of piRNA.

S6 Fig [a]
No obvious phenotype was detected in overexpression lines.

S7 Fig [a]
Immunostaining of BmVasa, BmSIWI and BmGTSF1 in . testis.

S8 Fig [a]
RNA-seq analysis reveals up-regulated transposons in gonads and qRT-PCR analysis reveals identical transposon levels in fat body.


Zdroje

1. Bachtrog D, Mank J, Peichel C, Kirkpatrick M, Otto S, Ashman T, et al. Sex Determination: Why So Many Ways of Doing It? PLoS biology. 2014;12(7):e1001899. doi: 10.1371/journal.pbio.1001899 24983465

2. Bopp D, Saccone G, Beye M. Sex determination in insects: variations on a common theme. Sex Dev. 2014;8(1–3):20–28. doi: 10.1159/000356458 24335049

3. Gempe T and Beye M. Function and evolution of sex determination mechanisms, genes and pathways in insects. Bioessays. 2011;33(1):52–60. doi: 10.1002/bies.201000043 21110346

4. Salz HK. Sex determination in insects: a binary decision based on alternative splicing. Curr Opin Genet Dev. 2011;21(4):395–400. doi: 10.1016/j.gde.2011.03.001 21474300

5. Salz HK and Erickson JW. Sex determination in Drosophila: The view from the top. Fly (Austin). 2010;4(1):60–70. doi: 10.4161/fly.4.1.11277 20160499

6. Erickson JW, Quintero JJ. Indirect effects of ploidy suggest X chromosome dose, not the X: A ratio, signals sex in Drosophila. PLoS Biol. 2007;5(12):e332. doi: 10.1371/journal.pbio.0050332 18162044

7. González AN, Lu H, Erickson JW. A shared enhancer controls a temporal switch between promoters during Drosophila primary sex determination. Proc Natl Acad Sci U S A. 2008;105(47):18436–18441. doi: 10.1073/pnas.0805993105 19011108

8. Harrison DA. Sex determination: controlling the master. Curr Biol. 2007;17(9):R328–330. doi: 10.1016/j.cub.2007.03.012 17470347

9. Siera SG, Cline TW. Sexual back talk with evolutionary implications: stimulation of the Drosophila sex-determination gene sex-lethal by its target transformer. Genetics. 2008;180(4):1963–1981. doi: 10.1534/genetics.108.093898 18845845

10. Evans DS, Cline TW. Drosophila switch gene Sex-lethal can bypass its switch-gene target transformer to regulate aspects of female behavior. Proc Natl Acad Sci U S A. 2013;110(47):E4474–E4481. doi: 10.1073/pnas.1319063110 24191002

11. Johnson ML, Nagengast AA, Salz HK. PPS, a large multidomain protein, functions with sex-lethal to regulate alternative splicing in Drosophila. PLoS Genet. 2010;6(3):e1000872. doi: 10.1371/journal.pgen.1000872 20221253

12. Moschall R, Gaik M, Medenbach J. Promiscuity in post-transcriptional control of gene expression: Drosophila sex-lethal and its regulatory partnerships. FEBS Lett. 2017;591(11):1471–1488. doi: 10.1002/1873-3468.12652 28391641

13. Traut W, Sahara K, Marec F. Sex chromosomes and sex determination in Lepidoptera. Sex Dev. 2007;1(6):332–346. doi: 10.1159/000111765 18391545

14. Abe H, Fujii T, Tanaka N, Yokoyama T, Kakehashi H, Ajimura M, et al. Identification of the female-determining region of the W chromosome in Bombyx mori. Genetica. 2008;133(3):269–282. doi: 10.1007/s10709-007-9210-1 17901928

15. Kiuchi T, Koga H, Kawamoto M, Shoji K, Sakai H, Arai Y, et al. A single female-specific piRNA is the primary determiner of sex in the silkworm. Nature. 2014;509(7502):633–636. doi: 10.1038/nature13315 24828047

16. Katsuma S, Kawamoto M, Kiuchi T. Guardian small RNAs and sex determination. RNA Biol. 2014;11(10):1238–1242. doi: 10.1080/15476286.2014.996060 25588029

17. Sakai H, Sumitani M, Chikami Y, Yahata K, Uchino K, Kiuchi T, et al. Transgenic Expression of the piRNA-Resistant Masculinizer Gene Induces Female-Specific Lethality and Partial Female-to-Male Sex Reversal in the Silkworm, Bombyx mori. PLoS Genet. 2016;12(8):e1006203. doi: 10.1371/journal.pgen.1006203 27579676

18. Li Z, You L, Yan D, James AA, Huang Y, Tan A. Bombyx mori histone methyltransferase BmAsh2 is essential for silkworm piRNA-mediated sex determination. PLoS Genet. 2017;14(2):e1007245. doi: 10.1371/journal.pgen.1007245 29474354

19. Muerdter F, Guzzardo PM, Gillis J, Luo Y, Yu Y, Chen C, et al. A genome-wide RNAi screen draws a genetic framework for transposon control and primary piRNA biogenesis in Drosophila. Mol Cell. 2013;50(5):736–748. doi: 10.1016/j.molcel.2013.04.006 23665228

20. Dönertas D, Sienski G, Brennecke J. Drosophila Gtsf1 is an essential component of the Piwi-mediated transcriptional silencing complex. Genes Dev. 2013;27(15):1693–1705. doi: 10.1101/gad.221150.113 23913922

21. Ohtani H, Iwasaki YW, Shibuya A, Siomi H, Siomi MC, Saito K. DmGTSF1 is necessary for Piwi-piRISC-mediated transcriptional transposon silencing in the Drosophila ovary. Genes Dev. 2013;27(15):1656–1661. doi: 10.1101/gad.221515.113 23913921

22. Yoshimura T, Watanabe T, Kuramochi-Miyagawa S, Takemoto N, Shiromoto Y, Kudo A, et al. Mouse GTSF1 is an essential factor for secondary piRNA biogenesis. EMBO Rep. 2018;19(4):e42054. doi: 10.15252/embr.201642054 29437694

23. Wang SH, Elgin SC. Drosophila Piwi functions downstream of piRNA production mediating a chromatin-based transposon silencing mechanism in female germ line. Proc Natl Acad Sci U S A. 2011;108(52):21164–211649. doi: 10.1073/pnas.1107892109 22160707

24. Tóth KF, Pezic D, Stuwe E, Webster A. The piRNA Pathway Guards the Germline Genome Against Transposable Elements. Adv Exp Med Biol. 2016;886:51–77. doi: 10.1007/978-94-017-7417-8_4 26659487

25. Iwasaki YW, Siomi MC, Siomi H. PIWI-Interacting RNA: Its Biogenesis and Functions. Annu Rev Biochem. 2015;84:405–433. doi: 10.1146/annurev-biochem-060614-034258 25747396

26. Sienski G, Dönertas D, Brennecke J. Transcriptional silencing of transposons by Piwi and maelstrom and its impact on chromatin state and gene expression. Cell. 2012;151(5):964–980. doi: 10.1016/j.cell.2012.10.040 23159368

27. Ozata DM, Gainetdinov I, Zoch A, O'Carroll D, Zamore PD. PIWI-interacting RNAs: small RNAs with big functions. PIWI-interacting RNAs: small RNAs with big functions. Nat Rev Genet. 2019;20(2):89–108. doi: 10.1038/s41576-018-0073-3 30446728

28. Siomi MC, Sato K, Pezic D, Aravin AA. PIWI-interacting small RNAs: the vanguard of genome defence. Nat Rev Mol Cell Biol. 2011;12(4):246–458. doi: 10.1038/nrm3089 21427766

29. Czech B, Munafò M, Ciabrelli F, Eastwood EL, Fabry MH, Kneuss E, et al. piRNA-Guided Genome Defense: From Biogenesis to Silencing. Annu Rev Genet. 2018;52:131–157. doi: 10.1146/annurev-genet-120417-031441 30476449

30. Théron E, Maupetit-Mehouas S, Pouchin P, Baudet L, Brasset E, Vaury C. The interplay between the Argonaute proteins Piwi and Aub within Drosophila germarium is critical for oogenesis, piRNA biogenesis and TE silencing. Nucleic Acids Res. 2018;46(19):10052–10065. doi: 10.1093/nar/gky695 30113668

31. Nishimura T, Nagamori I, Nakatani T, Izumi N, Tomari Y, Kuramochi-Miyagawa S, et al. PNLDC1, mouse pre-piRNA Trimmer, is required for meiotic and post-meiotic male germ cell development. EMBO Rep. 2018;19(3):e44957. doi: 10.15252/embr.201744957 29444933

32. Huang X, Fejes Tóth K, Aravin AA. piRNA Biogenesis in Drosophila melanogaster. Trends Genet. 2017;33(11):882–894. doi: 10.1016/j.tig.2017.09.002 28964526

33. Li C, Vagin VV, Lee S, Xu J, Ma S, Xi H, et al. Collapse of germline piRNAs in the absence of Argonaute3 reveals somatic piRNAs in flies. 2009;137(3):509–521. doi: 10.1016/j.cell.2009.04.027 19395009

34. Rozhkov NV, Hammell M, Hannon GJ. Multiple roles for Piwi in silencing Drosophila transposons. Genes Dev. 2013;27(4):400–412. doi: 10.1101/gad.209767.112 23392609

35. Zhang Z, Xu J, Koppetsch BS, Wang J, Tipping C, Ma S, et al. Heterotypic piRNA Ping-Pong requires qin, a protein with both E3 ligase and Tudor domains. Mol Cell. 2011;44(4):572–584. doi: 10.1016/j.molcel.2011.10.011 22099305

36. Matsumoto N, Nishimasu H, Sakakibara K, Nishida KM, Hirano T, Ishitani R, et al. Crystal Structure of Silkworm PIWI-Clade Argonaute Siwi Bound to piRNA. Cell. 2016;167(2):484–497. doi: 10.1016/j.cell.2016.09.002 27693359

37. Zhang Y, Guo R, Cui Y, Zhu Z, Zhang Y, Wu H, et al. An essential role for PNLDC1 in piRNA 3′ end trimming and male fertility in mice. Cell research. 2017;27(11):1392–1396. doi: 10.1038/cr.2017.125 28994417

38. Czech B, Hannon GJ. One Loop to Rule Them All: The Ping-Pong Cycle and piRNA-Guided Silencing. Trends Biochem Sci. 2016;41(4):324–337. doi: 10.1016/j.tibs.2015.12.008 26810602

39. Webster A, Li S, Hur JK, Wachsmuth M, Bois JS, Perkins EM, et al. Aub and Ago3 Are Recruited to Nuage through Two Mechanisms to Form a Ping-Pong Complex Assembled by Krimper. Mol Cell. 2015;59(4):564–575. doi: 10.1016/j.molcel.2015.07.017 26295961

40. Nishida KM, Iwasaki YW, Murota Y, Nagao A, Mannen T, Kato Y, et al. Respective functions of two distinct Siwi complexes assembled during PIWI-interacting RNA biogenesis in Bombyx germ cells. Cell Rep. 2015;10(2):193–203. doi: 10.1016/j.celrep.2014.12.013 25558067

41. Nishida KM, Sakakibara K, Iwasaki YW, Yamada H, Murakami R, Murota Y, et al. Hierarchical roles of mitochondrial Papi and Zucchini in Bombyx germline piRNA biogenesis. Nature. 2018;555(7695):260–264. doi: 10.1038/nature25788 29489748

42. Ge DT, Wang W, Tipping C, Gainetdinov I, Weng Z, Zamore PD. The RNA-Binding ATPase, Armitage, Couples piRNA Amplification in Nuage to Phased piRNA Production on Mitochondria. Mol Cell. 2019;74(5)982–995.e6. doi: 10.1016/j.molcel.2019.04.006 31076285

43. Horwich MD, Li C, Matranga C, Vagin V, Farley G, Wang P, et al. The Drosophila RNA methyltransferase, DmHen1, modifies germline piRNAs and single-stranded siRNAs in RISC. Curr Biol. 2007;17(14):1265–1272. doi: 10.1016/j.cub.2007.06.030 17604629

44. Chen K, Chen S, Xu J, Yu Y, Liu Z, Tan A, et al. Maelstrom regulates spermatogenesis of the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2019;109:43–51. doi: 10.1016/j.ibmb.2019.03.012 30970276

45. Xu J, Chen S, Zeng B, James AA, Tan A, Huang Y. Bombyx mori P-element Somatic Inhibitor (BmPSI) Is a Key Auxiliary Factor for Silkworm Male Sex Determination. PLoS Genet. 2017;13(1):e1006576. doi: 10.1371/journal.pgen.1006576 28103247

46. Hiroki S, Hiroyuki O, Kodai Y, Hiroki G, Takaaki D, Toshinobu Y, et al. Dimorphic Sperm Formation by Sex-lethal. Proc Natl Acad Sci U S A. 2019;116(21):10412–10417. doi: 10.1073/pnas.1820101116 31036645

47. Suzuki MG, Imanishi S, Dohmae N, Asanuma M, Matsumoto S. Identification of a male-specific RNA binding protein that regulates sex-specific splicing of Bmdsx by increasing RNA binding activity of BmPSI. Mol Cell Biol. 2010;30:5776–5786. doi: 10.1128/MCB.00444-10 20956562

48. Hara K, Fujii T, Suzuki Y, Sugano S, Shimada T, Katsuma S, et al. Altered expression of testis-specific genes, piRNAs, and transposons in the silkworm ovary masculinized by a W chromosome mutation. BMC Genomics. 2012;13:119. doi: 10.1186/1471-2164-13-119 22452797

49. Kawaoka S, Arai Y, Kadota K, Suzuki Y, Hara K, Sugano S, et al. Zygotic amplification of secondary piRNAs during silkworm embryogenesis. RNA. 2011;17(7):1401–1407. doi: 10.1261/rna.2709411 21628432

50. Almeida MV, Dietz S, Redl S, Karaulanov E, Hildebrandt A, Renz C, et al. GTSF-1 is required for formation of a functional RNA-dependent RNA Polymerase complex in Caenorhabditis elegans. EMBO J. 2018;37(12):e99325. doi: 10.15252/embj.201899325 29769402

51. Sakakibara K, Siomi MC. The PIWI-Interacting RNA Molecular Pathway: Insights From Cultured Silkworm Germline Cells. Bioessays. 2018;40(1). doi: 10.1002/bies.201700068 29164638

52. Ding D, Liu J, Dong K, Midic U, Hess RA, Xie H, et al. PNLDC1 is essential for piRNA 3' end trimming and transposon silencing during spermatogenesis in mice. Nat Commun. 2017;8(1):819. doi: 10.1038/s41467-017-00854-4 29018194

53. Aravin AA and Hannon GJ. Small RNA silencing pathways in germ and stem cells. Cold Spring Harb Symp Quant Biol. 2008;73:283–290. doi: 10.1101/sqb.2008.73.058 19270082

54. Khurana JS and Theurkauf W. piRNAs, transposon silencing, and Drosophila germline development. J Cell Biol. 2010;191(5):905–913. doi: 10.1083/jcb.201006034 21115802

55. Tan A, Fu G, Jin L, Guo Q, Li Z, Niu B, et al. Transgene-based, female-specific lethality system for genetic sexing of the silkworm, Bombyx mori. Proc Natl Acad Sci U S A. 2013;110(17):6766–6770. doi: 10.1073/pnas.1221700110 23569267

56. Yu Y, Liu XJ, Ma X, Zhang ZJ, Wang TC, Sun F, et al. A palmitoyltransferase Approximated gene Bm-app regulates wing development in Bombyx mori. Insect Sci. 2020;27(1):2–13. doi: 10.1111/1744-7917.12629 29943911

57. Ling L, Ge X, Li Z, Zeng B, Xu J, Chen X, et al. MiR-2 family targets awd and fng to regulate wing morphogenesis in Bombyx mori. RNA Biol, 2015;12(7):742–748. doi: 10.1080/15476286.2015.1048957 26037405

58. Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005;33(20):e179. doi: 10.1093/nar/gni178 16314309

59. Osanai-Futahashi M, Suetsugu Y, Mita K, Fujiwara H. Genome-wide screening and characterization of transposable elements and their distribution analysis in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2008;38(12):1046–1057. doi: 10.1016/j.ibmb.2008.05.012 19280695

60. Izumi N, Shoji K, Sakaguchi Y, Honda S, Kirino Y, Suzuki T, et al. Identification and Functional Analysis of the Pre-piRNA 3' Trimmer in Silkworms. Cell. 2016;164(5):962–973. doi: 10.1016/j.cell.2016.01.008 26919431

61. Kawaoka S, Kadota K, Arai Y, Suzuki Y, Fujii T, Abe H, et al. The silkworm W chromosome is a source of female-enriched piRNAs. RNA. 2011;17(12):2144–2151. doi: 10.1261/rna.027565.111 22020973


Článek vyšel v časopise

PLOS Genetics


2020 Číslo 11
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#