#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

A kinesin Klp10A mediates cell cycle-dependent shuttling of Piwi between nucleus and nuage


Authors: Zsolt G. Venkei aff001;  Charlotte Choi aff002;  Suhua Feng aff003;  Cuie Chen aff001;  Steven E. Jacobsen aff003;  John K. Kim aff002;  Yukiko M. Yamashita aff001;  Charlotte P. Choi aff002
Authors place of work: Life Sciences Institute, University of Michigan, Ann Arbor, Michigan, United States of America aff001;  Department of Biology, Johns Hopkins University, Baltimore, Maryland, United States of America aff002;  Department of Molecular, Cell and Developmental Biology, University of California, Los Angeles, California, United States of America aff003;  Eli and Edyth Broad Center of Regenerative Medicine and Stem Cell Research, University of California, Los Angeles, California, United States of America aff004;  Howard Hughes Medical Institute, University of California, Los Angeles, California, United States of America aff005;  Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, Michigan, United States of America aff006;  Howard Hughes Medical Institute, University of Michigan Ann Arbor, Michigan, United States of America aff007
Published in the journal: A kinesin Klp10A mediates cell cycle-dependent shuttling of Piwi between nucleus and nuage. PLoS Genet 16(3): e32767. doi:10.1371/journal.pgen.1008648
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008648

Summary

The piRNA pathway protects germline genomes from selfish genetic elements (e.g. transposons) through their transcript cleavage in the cytoplasm and/or their transcriptional silencing in the nucleus. Here, we describe a mechanism by which the nuclear and cytoplasmic arms of the piRNA pathway are linked. We find that during mitosis of Drosophila spermatogonia, nuclear Piwi interacts with nuage, the compartment that mediates the cytoplasmic arm of the piRNA pathway. At the end of mitosis, Piwi leaves nuage to return to the nucleus. Dissociation of Piwi from nuage occurs at the depolymerizing microtubules of the central spindle, mediated by a microtubule-depolymerizing kinesin, Klp10A. Depletion of klp10A delays the return of Piwi to the nucleus and affects piRNA production, suggesting the role of nuclear-cytoplasmic communication in piRNA biogenesis. We propose that cell cycle-dependent communication between the nuclear and cytoplasmic arms of the piRNA pathway may play a previously unappreciated role in piRNA regulation.

Keywords:

Transposable elements – Cell cycle and cell division – Drosophila melanogaster – Germ cells – MTS assay – Testes – RNA sequencing – Mitosis

Introduction

Piwi-associated RNAs (piRNAs), a class of endogenous small RNAs found across many organisms, associate with Piwi proteins of the Argonaute family to silence active transposable elements (TEs)[15]. In Drosophila melanogaster, Piwi, Aubergine (Aub), and Argonaute 3 (Ago3) comprise the three Piwi proteins of the piRNA pathway [1, 68]. The piRNA pathway is initiated in the nucleus, where primary piRNA precursors are transcribed from genomic piRNA clusters, translocated to the cytoplasm, and processed into mature piRNAs [1, 9]. In the cytoplasm of the germline, primary piRNAs are amplified in a “ping-pong” cycle that simultaneously cleaves TEs and produces secondary piRNAs [7, 10, 11]. Specifically, primary piRNAs antisense to TEs are loaded into and direct Aub to cleave TE transcripts, which produces secondary piRNAs that are subsequently loaded into Ago3. The secondary piRNAs generated by the ping-pong cycle are then loaded into Piwi and translocated into the nucleus as a piRNA-induced silencing complex (piRISC) [6, 12]. Within the nucleus, Piwi leads to the transcriptional repression of TEs through the deposition of heterochromatic histone marks [1315]. Thus, the cytoplasmic arm of the piRNA pathway intimately interacts with the nuclear arm to affect silencing at the post-transcriptional and transcriptional levels, respectively.

As evidence of the tight coupling between the two arms of the piRNA pathway, ping-pong amplification takes place in nuage–the electron-dense, phase-separated granules that are anchored directly to the cytoplasmic face of the nuclear pore [1618]. The interface between nuage and the nucleus has been shown to be critical for piRNA biogenesis and transposon target recognition and repression. Although nuage has been characterized as a static and long-term platform for piRNA biogenesis, individual nuage components, such as Aub and Ago3, are dynamically shuttled in and out of nuage [1, 6, 1921]. Unlike Aub and Ago3, Piwi shows predominantly nuclear localization [22] and has not been well-characterized within the context of nuage. Although previous studies have identified multiple piRNA pathway components that are required for proper localization of piRNA amplification machinery to nuage, including Spindle-E, Qin, Krimper, and Vasa, the underlying mechanism by which ping-pong-amplified piRNAs translocate into the nucleus remains poorly understood [6, 16, 21, 2325].

Here we show that nuclear and cytoplasmic components of the piRNA pathway interact specifically during mitosis of Drosophila spermatogonia (SGs) to likely facilitate communication between the nuclear and cytoplasmic arms of the pathway. We found that Piwi localizes to nuage specifically during mitosis, and returns to the nucleus upon mitotic exit. Through immunoprecipitation and mass spectrometry, we found that Klp10A, a microtubule (MT)-depolymerizing kinesin motor protein [26], physically interacts with components of the piRNA pathway. In the absence of Klp10A, Piwi remains in nuage for a prolonged period after mitotic exit, thus delaying its return to the nucleus. Cytological observations suggest that dissociation of Piwi from nuage at the end of mitosis occurs at the central spindle while it is being depolymerized in a klp10A-dependent manner. klp10A-depletion leads to increased piRNA production, supporting the physiological role of Klp10A-mediated regulation of Piwi localization. We propose that the interaction of Piwi with nuage during mitosis may represent a previously unappreciated mechanism to regulate piRNA biogenesis.

Results

Klp10A physically interacts with piRNA pathway components in germline stem cells and spermatogonia of Drosophila testis

Drosophila spermatogenesis is supported by asymmetrically dividing germline stem cells (GSCs), which produce a self-renewing GSC and a gonialblast (GB) that initiates the differentiation program. GBs further undergo four rounds of mitotic divisions as spermatogonia (SGs), which then enter meiotic program as spermatocytes (SCs) (Fig 1A). Asymmetric divisions of Drosophila male GSCs, as well as several other systems, are mediated by stereotypical positioning of the centrosomes [2736]. Therefore, stem cell-specific centrosomal components are of significant interest in understanding the mechanisms that regulate asymmetric stem cell divisions.

Klp10A interacts with piRNA pathway proteins in GSCs and SGs.
Fig. 1. Klp10A interacts with piRNA pathway proteins in GSCs and SGs.
A) Drosophila spermatogenesis. The stem cell niche is formed by non-dividing somatic cells (hub, marked by asterisk). Germline stem cells (GSCs) are physically attached to the hub, and divide asymmetrically. Gonialblasts (GBs), the differentiating daughters of GSCs, undergo four rounds of mitotic divisions with incomplete cytokinesis. Resultant 16-cell spermatogonia (SGs) then enter meiotic prophase as spermatocytes (SCs). Klp10A is specifically enriched at the centrosomes of GSCs (indicated by red dots), and at the central spindles of GSCs/GBs/SGs (red lines) [37]. B) Vasa co-immunoprecipitates with Klp10A in the extract from Upd-expressing testes (nos-gal4>UAS-upd). L: lysate input, IP: immunoprecipitant. C) Klp10A co-immunoprecipitates with GFP-Aub in the extract from Upd-expressing testes (nos-gal4>UAS-upd UAS-gfp, vs. nos-gal4>UAS-upd UAS-gfp-aub). D) Klp10A co-immunoprecipites with Vasa in the extract from the Dpp-expressing testes (nos-gal4>UAS-dpp). E) piRNA pathway proteins Aub, Ago3, Piwi and Vasa co-immunoprecipitate with Klp10A in the extract from Dpp-expressing testes (nos-gal4>UAS-dpp).

In a previous study, we showed that Klp10A protein is enriched on the centrosomes specifically in GSCs, but not in differentiating SGs [37] (Fig 1A). In an attempt to isolate proteins that specifically localize to GSC centrosomes, we affinity-purified Klp10A from a GSC-enriched extract, using either specific antibody or anti-GFP antibody combined with expression of Klp10A-GFP (see Methods). To enrich for GSCs, the self-renewal factor Unpaired (Upd) was expressed (nos-gal4>UAS-upd). Upd is a ligand expressed by hub cells to activate the JAK-STAT pathway, which determines GSC self-renewal. It was shown that ectopic expression of Upd in germ cells (nos-gal4>UAS-upd) is sufficient to cause GSC-like tumors in the testis [38, 39]. Anti-Klp10A and anti-GFP immunoprecipitates were analyzed by mass-spectrometry (see Methods). Many MT-associated proteins were enriched in Klp10A pulldown samples (S1S3 Tables), validating the specificity of the Klp10A pull down experiments.

To our surprise, we found that components of the piRNA pathway, such as Vasa, Aub and Piwi, were enriched in the Klp10A pull-downs (S1S3 Tables, Table 1). The fact that many piRNA pathway components were consistently identified in Klp10A-pulldowns led us to speculate this interaction might be of significance. Indeed, immunoprecipitation using anti-Klp10A antibody and GSC-enriched extracts confirmed that Klp10A physically interacts with Vasa (Fig 1B). Similarly, anti-GFP antibody pull down of GFP-Aub confirmed the Aub-Klp10A interaction (Fig 1C). Although it remains elusive whether Klp10A directly binds to these proteins or perhaps through an intermediary such as RNA, the fact that Klp10A co-immunoprecipitates multiple piRNA pathway components indicates that their interaction may be of biological relevance, regardless of whether the interaction is direct or not. Similar interactions between Klp10A and piRNA pathway components (Vasa, Piwi, Aub, Ago3) were observed when extracts from SG tumor was used (nos>dpp, [4042]) (Fig 1D and 1E, S4 and S5 Tables, Table 1), suggesting that the Klp10A interaction with piRNA pathway components is not unique to GSCs. Cytological data confirmed that Klp10A interaction with piRNA pathway components is not limited to GSC but also occurs in SGs (see below). Accordingly, the study diverged from our initial intention to isolate GSC centrosome-specific components. However, we investigated the unexpected role of klp10A in the regulation of the piRNA pathway, as the study revealed an unappreciated mode of regulation of piRNA biogenesis.

Tab. 1. Enrichment of piRNA pathway proteins in GFP-Klp10A and Klp10A pull downs from GSC and SG tumor testes.
Enrichment of piRNA pathway proteins in GFP-Klp10A and Klp10A pull downs from GSC and SG tumor testes.
PSM: number of peptide spectrum matches. Ranking indicates PSM number-based ranking of the protein/all identified proteins of the pull down. piRNA pathway proteins that ranked among top 20% are indicated by bold.

Depletion of klp10A results in alteration of piRNA biogenesis

To explore whether Klp10A interaction with piRNA pathway components was functionally significant, we examined whether klp10A was required for piRNA biogenesis [43]. To characterize the role of klp10A in the piRNA pathway, we first performed deep sequencing of small RNAs from wild type testes or testes with germline-specific klp10A-knockdown (nos-gal4>UAS-klp10ATRiP.HMS00920, validated in our previous study [37], and hereafter referred to as klp10ARNAi) (S1A Fig). We defined piRNAs as reads of 23–29 nucleotides (nt) in length that did not map to microRNAs or ribosomal RNAs. The majority of reads that are 23–29 nucleotides in length start with a Uridine at the 5’ most position, which has shown to be a unique, conserved feature of mature piRNAs [44], whereas reads that are 20–22 nucleotides in length, defined as siRNAs in our study, do not show such enrichment (S2 Fig), validating our analysis.

To assess the global changes in piRNA expression upon klp10A loss, we profiled piRNA reads across the Drosophila transcriptome. klp10ARNAi caused specific upregulation of piRNAs mapping to repetitive elements and piRNA clusters, with minimal differences in piRNAs mapping to other genomic classes, or in other types of small RNAs such as microRNAs (Fig 2A, S3 Fig). These data suggest that klp10A may play a specific role in the piRNA pathway. When we analyzed the abundance of piRNAs mapping to TEs upon loss of klp10A, we observed significant upregulation of piRNA reads both antisense and sense to TEs (Fig 2B–2E) without any noticeable changes in their distribution (Fig 2F, S4 Fig). The total amount of piRNA per TE was upregulated up to ~4 times (Fig 2B), with upregulation of individual piRNA being more prominent (Fig 2C, S4 Fig).

piRNA production is upregulated in <i>klp10A</i><sup><i>RNAi</i></sup> testis.
Fig. 2. piRNA production is upregulated in klp10ARNAi testis.
A) piRNA reads normalized to total library reads at each genomic feature. piRNAs are significantly upregulated in repeat and piRNA cluster regions. Significance was calculated using a Student’s t-test. B) Scatter and bar plots comparing transposon mapped sense and antisense piRNA abundance in klp10ARNAi germ cells vs. wild-type. piRNA abundance was calculated by normalizing piRNA read counts mapping to each transposon to total miRNA hairpin reads. A change of greater than 1.5-fold and FDR of <0.1 were used as cutoffs for differential analysis. Significance for piRNA fold change expression was calculated using Wilcoxon signed-rank test. C) Density of sequenced piRNAs (blue: sense; red: antisense) across HOBO. Each piRNA read was normalized to miRNA hairpin abundance. D) Scatter plot comparing transposon abundance in klp10ARNAi versus wild-type (red: FDR<0.1).

The above data showed that both primary and secondary piRNAs were upregulated, suggesting that the ping-pong cycle itself remains intact in klp10ARNAi. To test this, we calculated the number of piRNAs with a ping-pong signature (i.e. antisense and sense read pairs that have 10-nucleotide complementary base pairing from their 5’ ends) in both wildtype and klp10ARNAi libraries (S5A Fig). We observed a significant bias for piRNA read pairs with 10-nucleotide complementarity in both wildtype and klp10ARNAi, suggesting that the ping-pong pathway indeed remains active upon loss of klp10A. Moreover, we found no significant changes in ping-pong ratios in klp10ARNAi compared to wildtype (S5B Fig), further confirming an intact ping-pong cycle.

To address if piRNA upregulation in klp10ARNAi leads to changes in TE expression, we performed mRNA-sequencing in wildtype and klp10ARNAi germ cells (S6 Fig). We found that most TEs did not show any significant changes in expression levels except for a small subset of TEs, which showed moderate downregulation upon loss of klp10A (Fig 2D). There were no clear common characteristics among these TEs that may explain as to why these TEs exhibits downregulation. Taken together, our results show that klp10A is involved in piRNA production and that Klp10A interaction with piRNA pathway components may have functional significance (see Discussion), prompting us to further examine the underlying mechanisms. It should be noted that klp10ARNAi does not lead to noticeable changes in overall cellular composition of the testes, the process of differentiation [37], or cell cycle distribution (see below). Therefore, it is unlikely that skewed composition of cell types/cell cycle is the cause of observed changes in piRNA levels in klp10ARNAi.

Klp10A colocalizes with piRNA pathway components at the central spindle during telophase in GSCs and SGs

To begin to explore how Klp10A might contribute to piRNA biogenesis, we examined the localization of Klp10A protein and piRNA pathway components throughout the cell cycle of GSCs and SGs. Vasa and Aub showed well-established perinuclear localization in nuage throughout interphase (Fig 3A, 3B, 3E and 3F) [8, 45, 46]. During this period, Klp10A showed centrosome localization in GSCs (Fig 3A, 3E and 3I), or uniform cytoplasmic localization in SGs (Fig 3B, 3F and 3J) as we reported previously [37]. We did not observe clear co-localization between centrosomal Klp10A and nuage, suggesting that the potential role of Klp10A in piRNA pathway is not related to its GSC-specific centrosomal localization.

Klp10A colocalizes with Vasa, Aub and Piwi at the central spindle in GSCs/SGs during mitotic exit.
Fig. 3. Klp10A colocalizes with Vasa, Aub and Piwi at the central spindle in GSCs/SGs during mitotic exit.
A-D”) Localization of Klp10A (red) and GFP-Vasa (green) in GSCs/SGs throughout the cell cycle. DAPI (blue). E-H”) Localization of Klp10A (red) and GFP-Aub (green) in GSCs/SGs throughout the cell cycle, I-L”) Localization of Klp10A (red) and Piwi-GFP (green) in GSCs/SGs throughout the cell cycle. Centrosomal localization of Klp10A is indicated by arrowheads, central spindle localization with arrows. Bars: 5 μm. M-O) Quantification colocalization of GFP-Vasa (M), GFP-Aub (L) and GFP-Piwi (O) with Klp10A in GSCs/SGs during cell cycle. Mender’s colocalization coefficient is calculated based on how much of indicated proteins colocalize with Klp10A. Data points represent single cells, each from different testes. Error bars represent SD.

When cells enter mitosis, nuage became somewhat larger in size, which we refer to as ‘mitotic nuage’ hereafter (Fig 3C, 3G and 3K). Mitotic nuage appears to be the same population described previously as ‘Vasa granule near the mitotic chromosomes’ by Pek and Kai [47]. At this point, Klp10A was observed at the spindle pole (and weakly on the spindle) as previously reported [37], showing little colocalization with Vasa/Aub (Fig 3C, 3G, 3K and 3M–3O).

Although Klp10A did not noticeably colocalize with any of nuage components until anaphase, colocalization became clear during telophase: nuage was associated with Klp10A near the center of the dividing cell (Fig 3D, 3H and 3L–3O), which we confirmed to be bundles of central spindle MTs (S7 Fig). Importantly, Klp10A localization to the central spindle was observed both in GSCs and SGs [37], and colocalization of Klp10A and piRNA pathway components was commonly observed both in GSCs and SGs.

In contrast to Aub and Ago3, which reside in nuage together with Vasa to function in the cytoplasmic arm of the piRNA pathway, Piwi functions in the nucleus to repress TEs at the transcriptional level [15, 16, 43, 4850]. We confirmed that GFP-Piwi localizes to the nucleus during interphase of GSCs and SGs as described previously [22, 51](Fig 3I and 3J). Interestingly, we found that once cells enter mitosis, GFP-Piwi localized to mitotic nuage together with the components of the cytoplasmic arm of the piRNA pathway, such as Aub and Vasa (Fig 3K, see Fig 4 for the confirmation of Piwi’s localization to mitotic nuage). During telophase, GFP-Piwi still colocalized with nuage components at the central spindle (Fig 3L and 3O), after which it returned to the nucleus. Taken together, these data show that piRNA pathway components likely interact with Klp10A at the central spindle during telophase.

Piwi interacts with cytoplasmic piRNA pathway components at the nuage in mitosis, and returns to the nucleus at the end of mitosis in a Klp10A-dependent manner.
Fig. 4. Piwi interacts with cytoplasmic piRNA pathway components at the nuage in mitosis, and returns to the nucleus at the end of mitosis in a Klp10A-dependent manner.
(A-C) GFP-Piwi (green) and Vasa (magenta) localization throughout the GSC cell cycle in wild type testes. The localization patterns were the same in GSCs and SGs. GSCs are encircled by dotted lines. Cartoons show DAPI and MT morphology (shown in S8 Fig) to assess cell cycle stages. GFP-Piwi and Vasa colocalization at the mitotic nuage is indicated by arrowheads. Bars: 5 μm. D) co-immunoprecipitation of Vasa with GFP-Piwi from nos>upd testes after 4h 30min ex vivo MG132 treatment. nos>upd, gfp testes were used as control. L: lysate input, IP: immunoprecipitant. (E-G) GFP-Piwi (green) and Vasa (red) localization throughout the GSC cell cycle in klp10ARNAi testes. Cytoplasmic GFP-Piwi/Vasa granule in post-mitotic interphase cells are highlighted by arrows in panel G. H) co-immunoprecipitation of Vasa with GFP-Piwi from nos>upd testes after klp10ARNAi. nos>upd, gfp and nos>upd, gfp, klp10ARNAi testes were used as control.

Klp10A is required for relocation of Piwi from mitotic nuage to the nucleus at the end of mitosis

Based on the above results that suggest Klp10A interaction with nuage components at the central spindle, we explored how this interaction may impact the piRNA pathway. To this end, we first examined the localization of Piwi and Vasa in control vs. klp10ARNAi GSCs/SGs.

We confirmed the above observation that Piwi localizes to mitotic nuage by examining the colocalization of GFP-Piwi with Vasa in wild type cells. Although their localization was distinct during interphase (i.e. Vasa in nuage, Piwi in the nucleus) (Fig 3, Fig 4A and 4C, S8 Fig), Vasa and Piwi colocalized during mitosis (Fig 4B, S8 Fig, S9 Fig), and Piwi returned to the nucleus at the end of mitosis (Fig 4C, S10 Fig). We further detected physical interaction and colocalization between Vasa and Piwi specifically when cells were blocked in metaphase with MG132 (a proteasome inhibitor that arrests cells in mitosis [52]), or with colcemid (S9 Fig) [31]. These results suggest that Piwi gains access to nuage specifically during mitosis.

Next, we investigated whether klp10ARNAi testes show any defects in the behavior of nuage or Piwi during the cell cycle. Nuage morphology and composition appeared unperturbed in most interphase cells (Fig 4E), and in all mitotic cells until anaphase in both GSCs and SGs of klp10ARNAi testes (Fig 4F, S10A Fig). However, the difference between the control and klp10ARNAi germ cells became clear upon exiting mitosis (‘G1’ in Fig 4C and 4G). In control germ cells, Piwi returned to the nucleus at the end of telophase and became exclusively nuclear in the subsequent interphase as described above (Fig 4C, S10 Fig, S11A Fig). In contrast, in klp10ARNAi germ cells, Piwi persisted in cytoplasmic nuage even after the completion of mitosis (Fig 4G, S10 Fig, S11B and S11C Fig). Piwi-Vasa interaction was detectable by co-immunoprecipitation in klp10ARNAi without being arrested in mitosis (Fig 4H), consistent with persistent colocalization of Piwi with Vasa at nuage in early G1/S cells (Fig 4G, S11B and S11C Fig). Likely reflecting the delayed/incomplete translocation of Piwi from mitotic nuage to the nucleus, we also observed that the amount of nuclear Piwi was slightly but significantly reduced in GSCs/SGs in klp10ARNAi germ cells (S12 Fig). Delayed dissociation of Piwi from mitotic nuage at the end of mitosis was confirmed by live observation of GFP-Piwi and mCherry-Vasa expressed in mitotic SGs (S1 and S2 Videos). Taken together, these results show that 1) Piwi interacts with nuage specifically in mitosis, and 2) klp10A is required for releasing Piwi from mitotic nuage to allow its return to the nucleus at the end of mitosis.

Piwi dissociates from nuage at the central spindle microtubules

How does Piwi dissociate from mitotic nuage at the end of mitosis, and how does klp10A promote this process? A closer examination of live imaging using GFP-Piwi and mCherry-Vasa (S1 Video) revealed that a Piwi-positive compartment (‘Piwi granule’ hereafter) budded off from mitotic nuage and then dispersed near the center of telophase cells at the central spindle (see below for the confirmation of their localization at the central spindle) (Fig 5A, S1 Video). Concomitantly, nuclear Piwi level gradually increased (S1 Video), suggesting that Piwi released from mitotic nuage returned to the nucleus. Morphology of nuage marked by Vasa remained unchanged during this period (S1 Video). In contrast to control SGs, we barely observed the budding off of the Piwi granule in klp10ARNAi SGs (Fig 5B, S2 Video). As a result, Piwi remained in nuage for a prolonged time period even after the completion of mitosis. These results show that klp10A is required for releasing Piwi from mitotic nuage such that Piwi can return to the nucleus. It should be noted that klp10ARNAi germ cells did not exhibit noticeable change in cell cycle progression (see below). Therefore, it is unlikely that prolonged Piwi retention in nuage is an indirect consequence of delayed mitotic exit.

Piwi dissociates from mitotic nuage at the central spindle in a <i>klp10A</i>-dependent manner.
Fig. 5. Piwi dissociates from mitotic nuage at the central spindle in a klp10A-dependent manner.
(A-B) Live imaging of mCherry-Vasa (magenta) and GFP-Piwi (green) in the central region of telophase SGs in wild type (A) and in klp10ARNAi (B) testes. The elapsed time after anaphase-B onset is indicted in min:sec. Arrowheads point to Piwi-GFP signals budding off from the mitotic nuage. Dotted lines indicate the outline of cells. Bars: 1 μm. C-D) Live imaging of mCherry-α-Tubulin (magenta) and GFP-Piwi (green) in telophase SGs in wild type (C) and klp10ARNAi (D) testes. Squares in the first panels indicate the zoomed regions shown by the later panels. Arrowheads in (C) point to a mitotic nuage particle, initially positive for Piwi and gradually losing the signal. Open arrowheads (D) point to a Piwi-positive nuage particle, sliding along and later releasing the MT bundle, without decreasing GFP-Piwi signal. Bars: 5 μm, and 1 μm for zoomed regions. E-F) GFP-Piwi (green) localization in interphase control SGs (E) or mad2 null mutant SGs after 3h colcemid treatment (F). Arrowheads point to Piwi-GFP-positive nuage. Arrows point to interphase nuclei. Bars: 20 μm in left panels, 5 μm in insert panels. G) Frequency (cells/per testis) of interphase germ cells with GFP-Piwi positive nuage 3h after ex vivo colcemid treatment. n = 13–18 testes were scored for each experiment. p values of t-tests are provided.

The above live imaging observations indicated that budding off of the Piwi granule from mitotic nuage occurs near the center of dividing cells (S1 Video). Because Klp10A localizes to the central spindle and nuage components likely interact with Klp10A at the central spindle (Fig 3), we next examined how budding of the Piwi granule occurs at the central spindle. Using GFP-Piwi combined with mCherry-α-Tubulin, we found that the GFP-Piwi signal was dynamically associated with the central spindle MTs and moved back and forth along MTs, gradually decreasing in intensity until it completely dissipated (Fig 5C and S3 Video). In contrast, in klp10ARNAi GSCs/SGs, the majority of Piwi granules did not localize to the central spindle, and even when they are associated with the central spindle, Piwi did not exhibit gradual dispersion as observed in the control (Fig 5D, S4 Video).

These results indicate that budding off and release of Piwi from mitotic nuage occurs along the central spindle MTs. We hypothesized that the interaction of nuage with the central spindle MTs facilitates the dissociation of Piwi from nuage. To more directly test this idea that MTs mediate Piwi dissociation from mitotic nuage, we sought to allow cells to exit mitosis in the absence of MTs. To achieve this, we combined ex vivo colcemid treatment with a mad2 mutation. Colcemid treatment effectively depolymerized MTs (S9F Fig)[31], which would normally cause metaphase arrest due to the spindle assembly checkpoint (S9B and S9D Fig) [31, 53]. However, when this is combined with a mad2 mutation, which inactivates the spindle assembly checkpoint, cells exit mitosis even in the absence of MTs (S13 Fig) (Li et al. 2010). Under this condition, we found that Piwi remained in cytoplasmic nuage after mitotic exit, demonstrating that Piwi’s release from mitotic nuage and return to the nucleus is mediated by MTs (Fig 5E–5G).

klp10A is required for depolymerization of central spindle MTs

The above results led us to hypothesize that the central spindle serves as a platform for Piwi dissociation from mitotic nuage to allow Piwi to return to the nucleus. How does Klp10A, which is a kinesin motor that bends and depolymerizes MTs [26], regulate this process? Because central spindle MTs seemed to facilitate dissociation of Piwi from mitotic nuage (Fig 5), and because Klp10A localizes to the central spindle (S7 Fig), we investigated whether Klp10A may regulate the integrity of the central spindle.

To test this idea, we examined the morphology of the central spindle in control vs. klp10ARNAi testes. First, we found that the frequency of GSCs or SGs that contain central spindle MTs was significantly higher in klp10ARNAi compared to control testes (Fig 6A–6C), suggesting that Klp10A promotes depolymerization of central spindle MTs. In control germ cells, disassembly of the central spindle and completion of S phase, as assessed by pulse labeling with 5-ethynyl-2'-deoxyuridine (EdU), are nearly in synchrony: most of the S phase cells (EdU+) contain central spindle (Fig 6D, arrowhead), whereas all post-S phase cells (EdU-) have resolved central spindles (Fig 6F). In stark contrast, about half of post-S phase cells still contained the central spindle in klp10ARNAi testes (Fig 6E, arrow, Fig 6F), indicating that central spindle disassembly is delayed in the absence of klp10A. Importantly, increased frequency of germ cells with the central spindle is not due to skewed composition of cell cycle stages in klp10ARNAi germ cells, because klp10ARNAi germ cells had comparable frequency of being in S phase (22.6±8.6% EdU+ in control, n = 89 GSC vs. 24.6±15.6% EdU+ in klp10ARNAi, n = 103 GSC, n.s.) or in M phase (0.22±0.03 mitotic GSCs/testis in control, n = 160 vs. 0.194±0.01 mitotic GSCs/testis in klp10ARNAi, n = 222, n.s.).

Klp10A promotes depolymerization of central spindle MTs.
Fig. 6. Klp10A promotes depolymerization of central spindle MTs.
A, B) Apical tip region of wild type (A) and nos>klp10ARNAi (B) testes. Vasa (red), acetylated MTs (acMTs) (green) and DAPI (blue). Arrows indicate acMT bundles of central spindles between interconnected cells. Bars: 5 μm C) Frequency of acMT bundles in GSCs (n = 181 for control, n = 201 for klp10ARNAi) and 4-cell SGs (n = 128 and n = 95). Error bars indicate SD. p values of t-tests are provided. D-E) Apical tip region of a wild type (D) and a klp10ARNAi testis (E) with EdU (green) incorporation after 45min incubation period. Testes were stained for acMTs (red) and Vasa (blue). Arrowheads indicate EdU-positive GSC-GB pairs with acMT bundle. Arrow indicates EdU-negative GSC-GB pairs with acMT bundle in klp10ARNAi. Asterisks indicate the hub, dotted lines indicate GSC-GB pairs. Bars: 5 μm. F) Frequency of acMT bundle in EdU-positive (n = 32 for control, n = 24 for klp10ARNAi) and EdU-negative GSCs (n = 69 for control, n = 79 for klp10ARNAi). Error bars indicate SD. p values of t-tests are provided.

Taken together, these results show that klp10A is required for depolymerization of central spindle MTs. Based on the observation that Piwi dissociates from mitotic nuage along the central spindle as it disassembles, we speculate that Piwi’s dissociation from the central spindle is driven by MT depolymerization, which is facilitated by Klp10A, a MT-depolymerizing kinesin. Understanding how depolymerizing MTs can facilitate dissociation of Piwi from mitotic nuage awaits future investigation.

Discussion

In this study, we reveal an unexpected, dynamic localization of piRNA pathway machinery during the cell cycle of mitotically proliferating germ cells (GSCs and SGs) in the Drosophila testis. Our study is the first to describe a dynamic compositional change in nuage during the cell cycle. Piwi is nuclear in interphase, but associates with nuage specifically during mitosis and then returns to the nucleus at the end of the mitosis. Our study provides several novel insights into the cell biological aspects of piRNA biogenesis.

Our data suggest that this ‘nuage cycle’ in mitotic germ cells of the testis is partly facilitated by disassembly of central spindle MTs: Klp10A regulates depolymerization of central spindle MTs during telophase, and Piwi dissociates from nuage at depolymerizing central spindle MTs. How can depolymerizing MTs facilitate dissociation of Piwi from the nuage compartment? Although the underlying molecular mechanisms remain elusive, the physical interaction of piRNA pathway components and MTs have also been documented [5456]. In addition, MTs have been shown to regulate the assembly, maturation, and disassembly of stress granules, another phase-separated compartment similar to nuage [5760]. Given these earlier studies, our observation might reflect a generalizable regulatory feature of the piRNA pathway. Recently, it was shown that liquid droplets of Tau protein concentrate tubulin dimers to facilitate MT formation, and that MTs within the Tau droplets exhibit liquid-like properties [61]. Polymerized MTs within Tau droplets then deforms Tau droplets [61]. Likewise, tubulin-binding proteins within nuage [5456] may increase tubulin concentration, leading to MT polymerization. It might in turn deform nuage and change biophysical characteristics of phase separation, leading to release of Piwi from nuage. Although our results show the role of Klp10A in dissociating Piwi from nuage at the end of mitotic spindle, underlying mechanism of how Klp10A interacts with piRNA pathway components in facilitating Piwi’s dissociation on the depolymerizing MTs remains unknown. Klp10A protein does not have any noticeable domains other than central motor domain, and the structure does not hint as to how Klp10A might interact with piRNA pathway components. It is possible that Klp10A’s interaction with nuage components are mediated by unidentified linker protein(s). Also, it remains unknown how Piwi associates with nuage at the mitotic entry, which is likely independent of Klp10A. Nuclear envelope breakdown may be sufficient for Piwi to associate with nuage. Alternatively, there may be additional post-translational regulations that enhance Piwi’s affinity to nuage.

We demonstrate that Klp10A is involved in piRNA biogenesis. Loss of klp10A leads to a global upregulation of piRNA expression but only modest depletion of a subset of transposons. The failure of klp10A knockdown to show depletion of a majority of TEs, despite elevated levels of cognate piRNAs, could be because klp10A influences piRNA biogenesis only during limited stages of spermatogenesis, i.e. GSCs and SGs. For example, TE transcripts are not regulated by the piRNA pathway in germ cells of certain stages of spermatogenesis [62], and defects caused by klp10A depletion may be masked by unchanged TE expression in other stages of spermatogenesis. Our cytological data show that klp10A is required for the release of Piwi from the mitotic nuage, such that Piwi can re-enter the nucleus upon completion of mitosis. As a result of delayed return of Piwi to the nucleus, overall nuclear amount of Piwi was reduced in klp10ARNAi. How do these defects result in increased piRNA production? It is possible that the prolonged localization of Piwi to the nuage in klp10ARNAi may enhance the production of pre-piRNAs (or its loading to Piwi), leading to subsequent dysregulation of ping-pong activity to cause upregulation of piRNAs. Alternatively, because Piwi represses the expression of TEs and certain pre-piRNAs [63, 64], reduced nuclear Piwi might lead to derepression of TEs and pre-piRNAs, which, in turn, may result in enhanced piRNA production through intact ping-pong activity in klp10ARNAi SGs.

Expression of piRNAs in male germ cells appears to be distinct from those found in nurse cells during female germline development [65]. Nurse cells are highly polyploid post-mitotic cells and thus do not undergo nuclear envelope breakdown as do SGs [66]. In nurse cells of the Drosophila ovary, Piwi needs to be loaded with piRNAs to translocate into the nucleus [6, 67, 68]. However, Piwi localization to nuage is not visible in wild type cells, likely because Piwi only transiently associates with nuage prior to nuclear translocation. Only when the piRNA pathway is completely compromised, as in the aub ago3 double mutant, can the failure of Piwi to enter the nucleus be detected, leading to its accumulation in nuage [6]. In contrast to this mechanism of piRNA loading and nuclear translocation of Piwi in post-mitotic nurse cells, our finding suggests that mitotically dividing germ cells (SGs) in the testis may utilize mitosis (nuclear envelope breakdown) as a means to load Piwi with piRNAs. Nuclear envelope breakdown may represent a robust and efficient way of loading Piwi with piRNAs in SGs that divide every ~12 hours.

Taken together, we propose that Klp10A-dependent MT depolymerization at the central spindle facilitates Piwi dissociation from nuage to promote its translocation back into the nucleus. Piwi’s interaction with nuage during mitosis might impact piRNA biogenesis, as indicated by piRNA profiling. In summary, the present study reveals a novel cell biological mechanism by which nuclear and cytoplasmic arms of piRNA biogenesis communicate during the cell cycle.

Materials and methods

Fly husbandry and transgenic flies

Flies were raised in standard Bloomington medium at 25°C. The following stocks were used. UAS-gfp-klp10A [69], gfp-vas [70], gfp-piwi (a gift from Katalin Fejes Tóth) [15], mCherry-vas (a gift from Elizabeth Gavis)[71], and the following stocks, obtained from the Bloomington Stock Center: nos-gal4 [72], UAS-upd [73], UAS-gfp (Spana E., 1999, personal communication to FlyBase, ID: FBrf0111645), UAS-gfp-aub [74], UAS- mCherry-α-tubulin [75], UAS-klp10ATRiP.HMS00920 (Flybase: FBrf0214641, FBrf0212437), mad2EY21687 [53], Df(3L)BSC437 (FlyBase: FBrf0204472)

Immunoprecipitation and western blotting

For immunoprecipitation, 150 pairs of testes were lysed in 0.5 ml HEPES based buffer [50mM HEPES pH 7.5, 250mM NaCl, 0.1% Nonidet P40, 0.2% Triton X-100, supplied with cOmplete protease inhibitor mixture (Roche)]. In Klp10A co-immunoprecipitation (co-IP) experiments, the cleared lysates were incubated with rabbit anti-Klp10A serum (1:200, [37]) for 4 hours, extended by 2 extra hours after supplement with Protein A Dynabeads (Life Technologies). The beads were washed four times with lysis buffer and proteins were eluted with SDS-PAGE protein sample buffer. For co-IP experiments with GFP tagged proteins, the cleared lysates were incubated with GFP-Trap magnetic beads (Chromotek) according to manufacturer's instructions. To detect precipitated proteins on western-blots, the following primary antibodies were combined with horseradish peroxidase conjugated secondary antibodies (Abcam): guinea-pig anti-Klp10A (1:10,000, [37]), rabbit anti-GFP (1:3000, ab290, Abcam), rabbit anti-Vasa d-26 (1:3000, Santa Cruz Biotechnology), rabbit anti-Piwi (1:3000 [76]), rabbit anti-Aub (1:3000 [76]), and mouse anti-Ago3 (1:3000 [12]).

Mass-spectrometry

Klp10A was purified either with anti-Klp10A or anti-GFP nanobodies from GSC-enriched testes (250 pairs of nos-gal4>UAS-upd testes or 1500 pairs of nos-gal4>UAS-upd, UAS-gfp-klp10A testes, respectively) were dissected, and homogenized in lysis buffer (10 mM Tris-HCl pH 7.5; 150 mM NaCl; 0.5 mM EDTA, 0.1% NP40 and protease inhibitor cocktail (cOmplete, Roche)) for 30 min at 4°C. After centrifugation at 13,000 rpm for 10–15 min at 4°C the supernatant was saved as whole cell extract. The protein concentration was measured by absorbance at 562 nm using Pierce BCA Protein Assay Kit (Thermo Fisher Scientific). The whole cell extract of 1 μg total protein mass was incubated with 40 μl packed guinea-pig anti-Klp10A [37] conjugated Protein A-Dynabeads (Thermo Fisher Scientific) for 2h at 4°C, or with 25 μl packed GFP-Trap agarose beads (Chromotek) for 3 hours at 4°C. Beads were washed three times with wash buffer (10 mM Tris-HCl pH 7.5; 150 mM NaCl; 0.5 mM EDTA), and the proteins were eluted in LDS sample loading buffer (1.5x) at 100°C for 5 min and separated on a 10% Bis-Tris Novex mini-gel (Invitrogen) using the MES buffer system. In the anti-Klp10A pull down experiments the gel was stained with Coomassie dye and excised into ten equally sized segments. These segments were analyzed by LC/MS/MS (MS Bioworks, Ann Arbor, MI). The gel digests were analyzed by nano LC/MS/MS with a Waters NanoAcquity HPLC system interfaced to a Thermo Fisher Q Exactive. Peptides were loaded on a trapping column and eluted over a 75 μm analytical column at 350 nL/min; both columns were packed with Luna C18 resin (Phenomenex). The mass spectrometer was operated in data-dependent mode, with MS and MS/MS performed in the Orbitrap at 70,000 FWHM resolution and 17,500 FWHM resolution, respectively. In the anti-GFP pull down experiments after electrophoresis the gel was stained with coomassie and excised into two equally sized segments. The gel digests were analyzed by LC-MS/MS with LTQ Orbitrap XL mass spectrometer (Thermo Fisher Scientific) at Fred Hutchinson Cancer Research Center, Proteomics Resource.

Immunofluorescence staining and microscopy

Fixation and immunofluorescence staining of testes was performed as described previously [37]. The following primary antibodies were used. Mouse anti-Fasciclin III (FasIII; 1:100; 7G10, developed by Goodman C, obtained from Developmental Studies Hybridoma Bank; [77]), rabbit anti-Vasa (1:200; d-26, Santa Cruz Biotechnology), mouse anti-α-tubulin (4.3; 1:50; developed by C. Walsh [78] and obtained from the Developmental Studies Hybridoma Bank), mouse anti-acetylated-Tubulin (1:100, 6-11B-1, Sigma), rabbit anti-phosphorylated (Thr3) histone H3 (PH3) (1:200; clone JY325, Upstate), mouse anti-LaminB (1:200; C20; Santa Cruz Biotechnology), guinea pig anti-Traffic jam (Tj; 1:400; a gift from Dorothea Godt; [79]), rabbit and guinea-pig anti-Klp10A (1:3000/1:1000 respectively; [37]). Alexa Fluor-conjugated secondary antibodies were used (1:200; Thermo Fisher Scientific), EdU was detected by Click-iT Plus EdU Imaging Kit with Alexa Fluor 647 (Thermofisher). Testes were mounted into VECTASHIELD media with 4',6-diamidino-2-phenylindole (DAPI; Vector Labs). Images were captured using a Leica TCS SP8 confocal microscope with a 63×oil-immersion objective (NA = 1.4) and processed by ImageJ software [80]. Colocalization of nuage components with Klp10A or acMTs were quantified by calculating Mander`s coefficient [81, 82] in ImageJ.

Ex vivo treatment of Drosophila testis

To enrich cells in metaphase testes were dissected and transferred to Schneider's insect medium (Sigma) containing MG132 (20 μM final concentration, Sigma-Aldrich) [52] or colcemid (100μM final concentration, Calbiochem) [31]. After 4.5h incubation with MG132 or colcemid at 25°C, GSCs/SGs arrested in metaphase with or without intact bipolar mitotic spindles, respectively. For EdU incorporation the dissected testes were incubated for 45min (10 μM, Thermofisher) in Schneider's insect medium (Gibco) at 25°C.

RNA extraction, library preparations, and sequencing

Total RNAs were extracted by Trizol. 300 pairs of testes were used for each sample, and biological triplicates for each genotype were used. For small RNA sequencing, 2.5 μg of total RNAs were treated with DNase I (Amplification Grade, ThermoFisher) and recovered by RNA Clean & Concentrator-5 columns (Zymo). Libraries were generated from DNaseI treated RNAs using reagents from the TruSeq Small RNA Library Prep Kit (Illumina), while following a previously published protocol from [83] that uses a terminator oligo to deplete Drosophila 2S rRNA from the final sequencing libraries. The libraries were sequenced in one lane on a HiSeq 2500 machine (Illumina) at the UCLA Broad Stem Cell Research Center BioSequencing Core. The terminator oligo complimentary to Drosophila 2S rRNA was TAC AAC CCT CAA CCA TAT GTA GTC CAA GCA /3SpC3/ (Integrated DNA Technologies).

For mRNA sequencing, 2.5 μg of total RNAs were treated with DNase I (Amplification Grade, ThermoFisher) and recovered by RNA Clean & Concentrator-25 columns (Zymo). DNaseI treated RNAs were then treated with Ribo-Zero Gold rRNA Removal Kit (Illumina). Libraries were generated from rRNA depleted RNAs using TruSeq Stranded Total RNA Sample Prep Kit (Illumina). The libraries were sequenced by a HiSeq 4000 machine (Illumina) at the UCLA Broad Stem Cell Research Center BioSequencing Core.

RNA-seq data analysis

Three biological replicates of each condition were used to test for statistical significance and comparative analysis. Raw qseq files were converted to fastq format using custom Python and awk scripts. Small RNAs were clipped from 3’ adapter sequences using Trimmomatic 0.27 [84]. Reads mapping with at most 1 mismatch to rRNA and miRNA hairpin sequences were parsed out by Bowtie2 v2.2.3 [85]. The remaining reads were aligned to the D. melanogaster transcriptome (Dm3) with at most 1 mismatch allowed using Bowtie2 v2.2.3. Abundance estimation was done at annotated genes, piRNA clusters, and TEs using the eXpress pipeline within the piPipes suite [86, 87]. We used Flybase for annotation of protein coding genes and TEs. piRNA cluster annotations were done according to [1]. For piRNA analysis, we selected reads that were 23–29 nucleotides in length to ensure other endo-siRNA species did not influence downstream analysis. miRNA hairpin reads were used to normalize between libraries. Differential expression analysis was done in DESeq2 3.7 using the Wald test and adjusted p-values were corrected using the Benjamin-Hochberg procedure with an FDR threshold of 0.01[88]. Plots were generated using R. Sequencing data is available at GEO: accession number: GSE122596.

Ping-pong activity analysis

piRNA reads were mapped directly to the D. melanogaster transcriptome (Dm3) using Bowtie2 v2.2.3 [85]. BAM files from Bowtie2 were converted to BED files using several pipelines, including BEDTools 2.7 from the piPipes suite, which assigns nucleotide positions of each read across TE transcript bodies [86, 89]. In addition to nucleotide position, reads were categorized as sense and antisense. Sense was defined as piRNAs that were derived from cleavage of the TE mRNAs and antisense was defined as piRNAs that were antisense to annotated TE transcripts. Custom Python scripts were used to calculate the number of 5’ to 5’ complementarity between sense and antisense reads. Ping-pong ratios were calculated for each transposon feature across all libraries. To calculate the ping-pong ratio of each transposon, we used custom Python scripts to calculate the number of piRNAs in which sense piRNAs with an A at the 10-nt position or antisense piRNAs with a U at the 1-nt position showed 10 nt of complementarity from the 5’ end. We then divided the number of such pairs by the total number of piRNA reads. The resulting ratio allowed for quantification of ping-pong activity, without needing to normalize for library size. The calculation of ping-pong ratio was done using custom Python scripts. Plots were generated using R.

mRNA-seq data analysis

Raw qseq files were converted to fastq format using custom Python and awk scripts. Due to the repetitive nature of transposon sequences, reads were aligned to the D. melanogaster genome (Dm3) by STAR-2.6.0a while allowing for up to 100 multi-alignments per read [90]. Feature and abundance estimation were determined using TEtranscripts 2.0.3 [91]. TEtranscripts initially distributes multi-mapped reads evenly among potential matches and optimizes the distribution of multi-mapped reads using the Expectation Maximization approach. Variance was measured based on normalized gene expression counts and we performed principal components analysis and found samples derived from wild-type clustered separately from klp10ARNAi on Principle Component 1 (S6B Fig). Differential analysis was done using DESeq2 3.7 using the Wald test and p-values were corrected using the Benjamin-Hochberg procedure with an FDR threshold of 0.01 [88]. All Plots were generated using R. Pipeline of Analysis is described in (S1A and S1B Fig).

Live imaging

Testes from newly enclosed flies were dissected in Schneider’s Drosophila medium (Gibco) and prepared for live imaging as described previously [92]. The testis tips were placed into a drop of medium in a glass-bottom chamber and were covered by regenerated cellulose membrane (Spectrum Lab). The chamber was mounted on a three-axis computer-controlled piezoelectric stage. An inverted Leica TCS SP8 confocal microscope with a 63× oil immersion objective (NA = 1.4) was used for imaging. Live imaging was performed at ambient room temperature. Images were processed using ImageJ software.

Supporting information

S1 Fig [tif]
Pipelines of small RNA and mRNA sequence analyses.

S2 Fig [tif]
Distribution of 5’ nucleotide identity for reads of 23–29 nucleotides vs. 20–22 nucleotides in length in control and testes.

S3 Fig [a]
Analysis of small RNAs in testes.

S4 Fig [tif]
Characterization of piRNA density across TE transcripts.

S5 Fig [tif]
Ping-pong pathway activity remain unaltered after depletion.

S6 Fig [tif]
Characterization of RNAseq datasets.

S7 Fig [acmts]
Klp10A localization at the central spindle of GSCs/SGs.

S8 Fig [blue]
Identification of cell cycle stage for analysis of Piwi-Vasa colocalization.

S9 Fig [a]
Piwi-Vasa colocalization in mitotic GSCs/SGs.

S10 Fig [green]
GFP-Piwi localization changes during mitotic exit in GSCs/SGs.

S11 Fig [green]
Piwi stays in nuage after mitotic exit in germ cells.

S12 Fig [green]
Nuclear Piwi level is decreased in .

S13 Fig [green]
mutant SGs exit mitosis in the absence of MTs due to defective spindle assembly checkpoint.

S1 Video [green]
Dynamics and composition of mitotic nuage in wild type SGs.

S2 Video [green]
Dynamics and composition of mitotic nuage in a SG upon knockdown of .

S3 Video [green]
Piwi disassembles along the central spindle MTs in a wild type GSC during mitotic exit.

S4 Video [green]
Defective Piwi disassembly in late telophase/G1 SGs upon knockdown of .

S1 Table [xlsx]
Results of anti-GFP pull downs from , vs. , testes.

S2 Table [xlsx]
Results of preimmune IgG vs. anti-Klp10A pull downs from testes.

S3 Table [xlsx]
Results of repeat experiment of preimmune IgG vs. anti-Klp10A pull downs from testes.

S4 Table [xlsx]
Results of anti-Klp10A pull downs from testes.

S5 Table [xlsx]
Results of repeat experiment of anti-Klp10A pull downs from testes.

S1 Dataset [xlsx]
Numerical data associated with each graph.


Zdroje

1. Brennecke J, Aravin AA, Stark A, Dus M, Kellis M, Sachidanandam R, et al. Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell. 2007;128(6):1089–103. Epub 2007/03/10. doi: 10.1016/j.cell.2007.01.043 17346786.

2. Girard A, Sachidanandam R, Hannon GJ, Carmell MA. A germline-specific class of small RNAs binds mammalian Piwi proteins. Nature. 2006;442(7099):199–202. Epub 2006/06/06. doi: 10.1038/nature04917 16751776.

3. Batista PJ, Ruby JG, Claycomb JM, Chiang R, Fahlgren N, Kasschau KD, et al. PRG-1 and 21U-RNAs interact to form the piRNA complex required for fertility in C. elegans. Mol Cell. 2008;31(1):67–78. Epub 2008/06/24. S1097-2765(08)00391-2 [pii] doi: 10.1016/j.molcel.2008.06.002 18571452; PubMed Central PMCID: PMC2570341.

4. Aravin AA, Lagos-Quintana M, Yalcin A, Zavolan M, Marks D, Snyder B, et al. The small RNA profile during Drosophila melanogaster development. Dev Cell. 2003;5(2):337–50. Epub 2003/08/16. doi: 10.1016/s1534-5807(03)00228-4 12919683.

5. Ruby JG, Jan C, Player C, Axtell MJ, Lee W, Nusbaum C, et al. Large-Scale Sequencing Reveals 21U-RNAs and Additional MicroRNAs and Endogenous siRNAs in C. elegans. Cell. 2006;127(6):1193–207. doi: 10.1016/j.cell.2006.10.040 17174894

6. Wang W, Han BW, Tipping C, Ge DT, Zhang Z, Weng Z, et al. Slicing and Binding by Ago3 or Aub Trigger Piwi-Bound piRNA Production by Distinct Mechanisms. Mol Cell. 2015;59(5):819–30. doi: 10.1016/j.molcel.2015.08.007 26340424; PubMed Central PMCID: PMC4560842.

7. Gunawardane LS, Saito K, Nishida KM, Miyoshi K, Kawamura Y, Nagami T, et al. A slicer-mediated mechanism for repeat-associated siRNA 5' end formation in Drosophila. Science. 2007;315(5818):1587–90. doi: 10.1126/science.1140494 17322028.

8. Nishida KM, Saito K, Mori T, Kawamura Y, Nagami-Okada T, Inagaki S, et al. Gene silencing mechanisms mediated by Aubergine piRNA complexes in Drosophila male gonad. RNA. 2007;13(11):1911–22. doi: 10.1261/rna.744307 17872506; PubMed Central PMCID: PMC2040086.

9. Vagin VV, Sigova A, Li C, Seitz H, Gvozdev V, Zamore PD. A distinct small RNA pathway silences selfish genetic elements in the germline. Science. 2006;313(5785):320–4. doi: 10.1126/science.1129333 16809489.

10. 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. Cell. 2009;137(3):509–21. Epub 2009/04/28. doi: 10.1016/j.cell.2009.04.027 19395009; PubMed Central PMCID: PMC2768572.

11. Wang W, Yoshikawa M, Han Bo W, Izumi N, Tomari Y, Weng Z, et al. The Initial Uridine of Primary piRNAs Does Not Create the Tenth Adenine that Is the Hallmark of Secondary piRNAs. Molecular Cell. 2014;56(5):708–16. doi: 10.1016/j.molcel.2014.10.016 25453759

12. Senti KA, Jurczak D, Sachidanandam R, Brennecke J. piRNA-guided slicing of transposon transcripts enforces their transcriptional silencing via specifying the nuclear piRNA repertoire. Genes Dev. 2015;29(16):1747–62. doi: 10.1101/gad.267252.115 26302790; PubMed Central PMCID: PMC4561483.

13. Sienski G, Batki J, Senti KA, Donertas D, Tirian L, Meixner K, et al. Silencio/CG9754 connects the Piwi-piRNA complex to the cellular heterochromatin machinery. Genes Dev. 2015;29(21):2258–71. doi: 10.1101/gad.271908.115 26494711; PubMed Central PMCID: PMC4647559.

14. Klenov MS, Lavrov SA, Korbut AP, Stolyarenko AD, Yakushev EY, Reuter M, et al. Impact of nuclear Piwi elimination on chromatin state in Drosophila melanogaster ovaries. Nucleic Acids Res. 2014;42(10):6208–18. Epub 2014/05/02. doi: 10.1093/nar/gku268 24782529; PubMed Central PMCID: PMC4041442.

15. Le Thomas A, Rogers AK, Webster A, Marinov GK, Liao SE, Perkins EM, et al. Piwi induces piRNA-guided transcriptional silencing and establishment of a repressive chromatin state. Genes Dev. 2013;27(4):390–9. doi: 10.1101/gad.209841.112 23392610; PubMed Central PMCID: PMC3589556.

16. Zhang F, Wang J, Xu J, Zhang Z, Koppetsch BS, Schultz N, et al. UAP56 couples piRNA clusters to the perinuclear transposon silencing machinery. Cell. 2012;151(4):871–84. Epub 2012/11/13. doi: 10.1016/j.cell.2012.09.040 23141543; PubMed Central PMCID: PMC3499805.

17. Lim AK, Kai T. Unique germ-line organelle, nuage, functions to repress selfish genetic elements in Drosophila melanogaster. Proc Natl Acad Sci U S A. 2007;104(16):6714–9. doi: 10.1073/pnas.0701920104 17428915.

18. Malone CD, Brennecke J, Dus M, Stark A, McCombie WR, Sachidanandam R, et al. Specialized piRNA pathways act in germline and somatic tissues of the Drosophila ovary. Cell. 2009;137(3):522–35. Epub 2009/04/28. doi: 10.1016/j.cell.2009.03.040 19395010; PubMed Central PMCID: PMC2882632.

19. Ryazansky SS, Kotov AA, Kibanov MV, Akulenko NV, Korbut AP, Lavrov SA, et al. RNA helicase Spn-E is required to maintain Aub and AGO3 protein levels for piRNA silencing in the germline of Drosophila. Eur J Cell Biol. 2016;95(9):311–22. doi: 10.1016/j.ejcb.2016.06.001 27320195.

20. Snee MJ, Macdonald PM. Live imaging of nuage and polar granules: evidence against a precursor-product relationship and a novel role for Oskar in stabilization of polar granule components. J Cell Sci. 2004;117(Pt 10):2109–20. Epub 2004/04/20. doi: 10.1242/jcs.01059 15090597

21. Andress A, Bei Y, Fonslow BR, Giri R, Wu Y, Yates JR 3rd, et al. Spindle-E cycling between nuage and cytoplasm is controlled by Qin and PIWI proteins. J Cell Biol. 2016;213(2):201–11. Epub 2016/04/20. doi: 10.1083/jcb.201411076 27091448; PubMed Central PMCID: PMC5084268.

22. Gonzalez J, Qi H, Liu N, Lin H. Piwi Is a Key Regulator of Both Somatic and Germline Stem Cells in the Drosophila Testis. Cell Rep. 2015;12(1):150–61. doi: 10.1016/j.celrep.2015.06.004 26119740; PubMed Central PMCID: PMC4497877.

23. 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–84. doi: 10.1016/j.molcel.2011.10.011 22099305; PubMed Central PMCID: PMC3236501.

24. 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–75. doi: 10.1016/j.molcel.2015.07.017 26295961; PubMed Central PMCID: PMC4545750.

25. Sato K, Iwasaki YW, Shibuya A, Carninci P, Tsuchizawa Y, Ishizu H, et al. Krimper Enforces an Antisense Bias on piRNA Pools by Binding AGO3 in the Drosophila Germline. Mol Cell. 2015;59(4):553–63. Epub 2015/07/28. doi: 10.1016/j.molcel.2015.06.024 26212455.

26. Rogers GC, Rogers SL, Schwimmer TA, Ems-McClung SC, Walczak CE, Vale RD, et al. Two mitotic kinesins cooperate to drive sister chromatid separation during anaphase. Nature. 2004;427(6972):364–70. Epub 2003/12/19. doi: 10.1038/nature02256 14681690.

27. Yamashita YM, Jones DL, Fuller MT. Orientation of asymmetric stem cell division by the APC tumor suppressor and centrosome. Science. 2003;301(5639):1547–50. Epub 2003/09/13. doi: 10.1126/science.1087795 12970569.

28. Yamashita YM, Mahowald AP, Perlin JR, Fuller MT. Asymmetric inheritance of mother versus daughter centrosome in stem cell division. Science. 2007;315(5811):518–21. Epub 2007/01/27. doi: 10.1126/science.1134910 17255513; PubMed Central PMCID: PMC2563045.

29. Cheng J, Turkel N, Hemati N, Fuller MT, Hunt AJ, Yamashita YM. Centrosome misorientation reduces stem cell division during ageing. Nature. 2008;456(7222):599–604. Epub 2008/10/17. doi: 10.1038/nature07386 18923395; PubMed Central PMCID: PMC2712891.

30. Inaba M, Venkei ZG, Yamashita YM. The polarity protein Baz forms a platform for the centrosome orientation during asymmetric stem cell division in the Drosophila male germline. Elife. 2015;4. Epub 2015/03/21. doi: 10.7554/eLife.04960 25793442; PubMed Central PMCID: PMC4391501.

31. Venkei ZG, Yamashita YM. The centrosome orientation checkpoint is germline stem cell specific and operates prior to the spindle assembly checkpoint in Drosophila testis. Development. 2015;142(1):62–9. Epub 2014/12/07. doi: 10.1242/dev.117044 25480919.

32. Salzmann V, Chen C, Chiang CY, Tiyaboonchai A, Mayer M, Yamashita YM. Centrosome-dependent asymmetric inheritance of the midbody ring in Drosophila germline stem cell division. Mol Biol Cell. 2014;25(2):267–75. Epub 2013/11/15. doi: 10.1091/mbc.E13-09-0541 24227883; PubMed Central PMCID: PMC3890347.

33. Conduit PT, Raff JW. Cnn dynamics drive centrosome size asymmetry to ensure daughter centriole retention in drosophila neuroblasts. Curr Biol. 2010;20(24):2187–92. Epub 2010/12/15. S0960-9822(10)01521-6 [pii] doi: 10.1016/j.cub.2010.11.055 21145745.

34. Januschke J, Llamazares S, Reina J, Gonzalez C. Drosophila neuroblasts retain the daughter centrosome. Nat Commun. 2011;2:243. Epub 2011/03/17. ncomms1245 [pii] doi: 10.1038/ncomms1245 21407209.

35. Januschke J, Reina J, Llamazares S, Bertran T, Rossi F, Roig J, et al. Centrobin controls mother-daughter centriole asymmetry in Drosophila neuroblasts. Nature cell biology. 2013;15(3):241–8. Epub 2013/01/29. doi: 10.1038/ncb2671 23354166.

36. Wang X, Tsai JW, Imai JH, Lian WN, Vallee RB, Shi SH. Asymmetric centrosome inheritance maintains neural progenitors in the neocortex. Nature. 2009;461(7266):947–55. Epub 2009/10/16. nature08435 [pii] doi: 10.1038/nature08435 19829375.

37. Chen C, Inaba M, Venkei ZG, Yamashita YM. Klp10A, a stem cell centrosome-enriched kinesin, balances asymmetries in Drosophila male germline stem cell division. Elife. 2016;5. Epub 2016/11/26. doi: 10.7554/eLife.20977 27885983; PubMed Central PMCID: PMC5235350.

38. Kiger AA, Jones DL, Schulz C, Rogers MB, Fuller MT. Stem cell self-renewal specified by JAK-STAT activation in response to a support cell cue. Science. 2001;294(5551):2542–5. doi: 10.1126/science.1066707 11752574.

39. Tulina N, Matunis E. Control of stem cell self-renewal in Drosophila spermatogenesis by JAK-STAT signaling. Science. 2001;294(5551):2546–9. doi: 10.1126/science.1066700 11752575.

40. Schulz C, Kiger AA, Tazuke SI, Yamashita YM, Pantalena-Filho LC, Jones DL, et al. A misexpression screen reveals effects of bag-of-marbles and TGF beta class signaling on the Drosophila male germ-line stem cell lineage. Genetics. 2004;167(2):707–23. Epub 2004/07/09. doi: 10.1534/genetics.103.023184 15238523; PubMed Central PMCID: PMC1470893.

41. Kawase E, Wong MD, Ding BC, Xie T. Gbb/Bmp signaling is essential for maintaining germline stem cells and for repressing bam transcription in the Drosophila testis. Development. 2004;131(6):1365–75. doi: 10.1242/dev.01025 14973292.

42. Shivdasani AA, Ingham PW. Regulation of stem cell maintenance and transit amplifying cell proliferation by tgf-Beta signaling in Drosophila spermatogenesis. Curr Biol. 2003;13(23):2065–72. doi: 10.1016/j.cub.2003.10.063 14653996.

43. Toth 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; PubMed Central PMCID: PMC4991928.

44. Aravin A, Gaidatzis D, Pfeffer S, Lagos-Quintana M, Landgraf P, Iovino N, et al. A novel class of small RNAs bind to MILI protein in mouse testes. Nature. 2006;442(7099):203–7. doi: 10.1038/nature04916 16751777

45. Kibanov MV, Egorova KS, Ryazansky SS, Sokolova OA, Kotov AA, Olenkina OM, et al. A novel organelle, the piNG-body, in the nuage of Drosophila male germ cells is associated with piRNA-mediated gene silencing. Mol Biol Cell. 2011;22(18):3410–9. doi: 10.1091/mbc.E11-02-0168 21775629; PubMed Central PMCID: PMC3172265.

46. Nagao A, Mituyama T, Huang H, Chen D, Siomi MC, Siomi H. Biogenesis pathways of piRNAs loaded onto AGO3 in the Drosophila testis. RNA. 2010. Epub 2010/10/29. rna.2270710 [pii] doi: 10.1261/rna.2270710 20980675

47. Pek JW, Kai T. A role for vasa in regulating mitotic chromosome condensation in Drosophila. Current biology: CB. 2011;21(1):39–44. Epub 2010/12/28. doi: 10.1016/j.cub.2010.11.051 21185189.

48. Aravin AA, Sachidanandam R, Bourc'his D, Schaefer C, Pezic D, Toth KF, et al. A piRNA pathway primed by individual transposons is linked to de novo DNA methylation in mice. Mol Cell. 2008;31(6):785–99. Epub 2008/10/17. S1097-2765(08)00619-9 [pii] doi: 10.1016/j.molcel.2008.09.003 18922463; PubMed Central PMCID: PMC2730041.

49. Xiol J, Spinelli P, Laussmann MA, Homolka D, Yang Z, Cora E, et al. RNA Clamping by Vasa Assembles a piRNA Amplifier Complex on Transposon Transcripts. Cell. 2014;157(7):1698–711. doi: 10.1016/j.cell.2014.05.018 24910301.

50. 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.

51. Cox DN, Chao A, Lin H. piwi encodes a nucleoplasmic factor whose activity modulates the number and division rate of germline stem cells. Development. 2000;127(3):503–14. 10631171.

52. Moutinho-Pereira S, Matos I, Maiato H. Drosophila S2 cells as a model system to investigate mitotic spindle dynamics, architecture, and function. Methods Cell Biol. 2010;97:243–57. doi: 10.1016/S0091-679X(10)97014-3 20719275.

53. Li D, Morley G, Whitaker M, Huang JY. Recruitment of Cdc20 to the kinetochore requires BubR1 but not Mad2 in Drosophila melanogaster. Mol Cell Biol. 2010;30(13):3384–95. doi: 10.1128/MCB.00258-10 20421417; PubMed Central PMCID: PMC2897573.

54. Lau NC, Ohsumi T, Borowsky M, Kingston RE, Blower MD. Systematic and single cell analysis of &lt;em&gt;Xenopus&lt;/em&gt; Piwi‐interacting RNAs and Xiwi. The EMBO Journal. 2009;28(19):2945. doi: 10.1038/emboj.2009.237 19713941

55. Rodriguez AJ, Seipel SA, Hamill DR, Romancino DP, Di Carlo M, Suprenant KA, et al. Seawi—a sea urchin piwi/argonaute family member is a component of MT-RNP complexes. RNA. 2005;11(5):646–56. doi: 10.1261/rna.7198205 15840816

56. Sato K, Nishida KM, Shibuya A, Siomi MC, Siomi H. Maelstrom coordinates microtubule organization during Drosophila oogenesis through interaction with components of the MTOC. Genes Dev. 2011;25(22):2361–73. doi: 10.1101/gad.174110.111 22085963; PubMed Central PMCID: PMC3222902.

57. Ivanov PA, Chudinova EM, Nadezhdina ES. Disruption of microtubules inhibits cytoplasmic ribonucleoprotein stress granule formation. Experimental Cell Research. 2003;290(2):227–33. doi: 10.1016/s0014-4827(03)00290-8 14567982

58. Loschi M, Leishman CC, Berardone N, Boccaccio GL. Dynein and kinesin regulate stress-granule and P-body dynamics. Journal of cell science. 2009;122(Pt 21):3973–82. Epub 2009/10/13. doi: 10.1242/jcs.051383 19825938.

59. Chernov KG, Barbet A, Hamon L, Ovchinnikov LP, Curmi PA, Pastré D. Role of Microtubules in Stress Granule Assembly: MICROTUBULE DYNAMICAL INSTABILITY FAVORS THE FORMATION OF MICROMETRIC STRESS GRANULES IN CELLS. Journal of Biological Chemistry. 2009;284(52):36569–80. doi: 10.1074/jbc.M109.042879 19843517

60. Shao J, Gao F, Zhang B, Zhao M, Zhou Y, He J, et al. Aggregation of SND1 in Stress Granules is Associated with the Microtubule Cytoskeleton During Heat Shock Stimulus. The Anatomical Record. 2017;300(12):2192–9. doi: 10.1002/ar.23642 28758359

61. Hernández-Vega A, Braun M, Scharrel L, Jahnel M, Wegmann S, Hyman BT, et al. Local Nucleation of Microtubule Bundles through Tubulin Concentration into a Condensed Tau Phase. Cell Rep. 2017;20(10):2304–12. doi: 10.1016/j.celrep.2017.08.042 28877466.

62. Vourekas A, Zheng Q, Alexiou P, Maragkakis M, Kirino Y, Gregory BD, et al. Mili and Miwi target RNA repertoire reveals piRNA biogenesis and function of Miwi in spermiogenesis. Nat Struct Mol Biol. 2012;19(8):773–81. Epub 2012/07/31. doi: 10.1038/nsmb.2347 22842725; PubMed Central PMCID: PMC3414646.

63. Chang TH, Mattei E, Gainetdinov I, Colpan C, Weng Z, Zamore PD. Maelstrom Represses Canonical Polymerase II Transcription within Bi-directional piRNA Clusters in Drosophila melanogaster. Molecular Cell. 2019;73(2):291–303.e6. doi: 10.1016/j.molcel.2018.10.038 30527661

64. Sato K, Siomi MC. Two distinct transcriptional controls triggered by nuclear Piwi-piRISCs in the Drosophila piRNA pathway. Current Opinion in Structural Biology. 2018;53:69–76. doi: 10.1016/j.sbi.2018.06.005 29990672

65. Malone CD, Hannon GJ. Molecular Evolution of piRNA and Transposon Control Pathways in Drosophila. Cold Spring Harbor Symposia on Quantitative Biology. 2009;74:225–34. doi: 10.1101/sqb.2009.74.052 20453205

66. Spradling A. Developmental Genetics of Oogenesis. Martinez-Arias MBaA, editor: New York: Cold Spring Harbor Laboratory Press.; 1993.

67. Handler D, Olivieri D, Novatchkova M, Gruber FS, Meixner K, Mechtler K, et al. A systematic analysis of Drosophila TUDOR domain-containing proteins identifies Vreteno and the Tdrd12 family as essential primary piRNA pathway factors. The EMBO journal. 2011;30(19):3977–93. Epub 2011/08/25. doi: 10.1038/emboj.2011.308 21863019; PubMed Central PMCID: PMC3209783.

68. Olivieri D, Sykora MM, Sachidanandam R, Mechtler K, Brennecke J. An in vivo RNAi assay identifies major genetic and cellular requirements for primary piRNA biogenesis in Drosophila. EMBO J. 2010;29(19):3301–17. doi: 10.1038/emboj.2010.212 20818334; PubMed Central PMCID: PMC2957214.

69. Inaba M, Buszczak M, Yamashita YM. Nanotubes mediate niche-stem-cell signalling in the Drosophila testis. Nature. 2015;523(7560):329–32. Epub 2015/07/02. doi: 10.1038/nature14602 26131929; PubMed Central PMCID: PMC4586072.

70. Sano H, Nakamura A, Kobayashi S. Identification of a transcriptional regulatory region for germline-specific expression of vasa gene in Drosophila melanogaster. Mech Dev. 2002;112(1–2):129–39. doi: 10.1016/s0925-4773(01)00654-2 11850184.

71. Lerit DA, Gavis ER. Transport of germ plasm on astral microtubules directs germ cell development in Drosophila. Current biology: CB. 2011;21(6):439–48. Epub 2011/03/08. doi: 10.1016/j.cub.2011.01.073 21376599; PubMed Central PMCID: PMC3062663.

72. Van Doren M, Williamson AL, Lehmann R. Regulation of zygotic gene expression in Drosophila primordial germ cells. Curr Biol. 1998;8(4):243–6. doi: 10.1016/s0960-9822(98)70091-0 9501989.

73. Zeidler MP, Perrimon N, Strutt DI. Polarity determination in the Drosophila eye: a novel role for unpaired and JAK/STAT signaling. Genes Dev. 1999;13(10):1342–53. Epub 1999/05/27. doi: 10.1101/gad.13.10.1342 10346822; PubMed Central PMCID: PMC316719.

74. Harris AN, Macdonald PM. Aubergine encodes a Drosophila polar granule component required for pole cell formation and related to eIF2C. Development. 2001;128(14):2823–32. 11526087.

75. Rusan NM, Peifer M. A role for a novel centrosome cycle in asymmetric cell division. J Cell Biol. 2007;177(1):13–20. doi: 10.1083/jcb.200612140 17403931.

76. Klattenhoff C, Xi H, Li C, Lee S, Xu J, Khurana JS, et al. The Drosophila HP1 homolog Rhino is required for transposon silencing and piRNA production by dual-strand clusters. Cell. 2009;138(6):1137–49. doi: 10.1016/j.cell.2009.07.014 19732946; PubMed Central PMCID: PMC2770713.

77. Patel NH, Snow PM, Goodman CS. Characterization and cloning of fasciclin III: a glycoprotein expressed on a subset of neurons and axon pathways in Drosophila. Cell. 1987;48(6):975–88. doi: 10.1016/0092-8674(87)90706-9 3548998.

78. Walsh C. Synthesis and assembly of the cytoskeleton of Naegleria gruberi flagellates. J Cell Biol. 1984;98(2):449–56. doi: 10.1083/jcb.98.2.449 6363422; PubMed Central PMCID: PMC2113116.

79. Li MA, Alls JD, Avancini RM, Koo K, Godt D. The large Maf factor Traffic Jam controls gonad morphogenesis in Drosophila. Nat Cell Biol. 2003;5(11):994–1000. doi: 10.1038/ncb1058 14578908.

80. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676–82. doi: 10.1038/nmeth.2019 22743772; PubMed Central PMCID: PMC3855844.

81. Dunn KW, Kamocka MM, McDonald JH. A practical guide to evaluating colocalization in biological microscopy. American Journal of Physiology-Cell Physiology. 2011;300(4):C723–C42. doi: 10.1152/ajpcell.00462.2010 21209361

82. Adler J, Parmryd I. Quantifying colocalization by correlation: The Pearson correlation coefficient is superior to the Mander's overlap coefficient. Cytometry Part A. 2010;77A(8):733–42. doi: 10.1002/cyto.a.20896 20653013

83. Wickersheim ML, Blumenstiel JP. Terminator oligo blocking efficiently eliminates rRNA from Drosophila small RNA sequencing libraries. Biotechniques. 2013;55(5):269–72. doi: 10.2144/000114102 24215643

84. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. doi: 10.1093/bioinformatics/btu170 24695404

85. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nature Methods. 2012;9:357. doi: 10.1038/nmeth.1923 https://www.nature.com/articles/nmeth.1923#supplementary-information. 22388286

86. Han BW, Wang W, Zamore PD, Weng Z. piPipes: a set of pipelines for piRNA and transposon analysis via small RNA-seq, RNA-seq, degradome- and CAGE-seq, ChIP-seq and genomic DNA sequencing. Bioinformatics. 2015;31(4):593–5. doi: 10.1093/bioinformatics/btu647 25342065

87. Roberts A, Pachter L. Streaming fragment assignment for real-time analysis of sequencing experiments. Nature Methods. 2012;10:71. doi: 10.1038/nmeth.2251 https://www.nature.com/articles/nmeth.2251#supplementary-information. 23160280

88. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology. 2014;15(12):550. doi: 10.1186/s13059-014-0550-8 25516281

89. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2. doi: 10.1093/bioinformatics/btq033 20110278

90. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21. doi: 10.1093/bioinformatics/bts635 23104886

91. Jin Y, Tam OH, Paniagua E, Hammell M. TEtranscripts: a package for including transposable elements in differential expression analysis of RNA-seq datasets. Bioinformatics. 2015;31(22):3593–9. doi: 10.1093/bioinformatics/btv422 26206304

92. Cheng J, Hunt AJ. Time-lapse live imaging of stem cells in Drosophila testis. Curr Protoc Stem Cell Biol. 2009;Chapter 2:Unit 2E doi: 10.1002/9780470151808.sc02e02s11 19885824; PubMed Central PMCID: PMC2778854.


Článek vyšel v časopise

PLOS Genetics


2020 Číslo 3
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
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#