- Research article
- Open Access
Kinesin-2 and kinesin-9 have atypical functions during ciliogenesis in the male gametophyte of Marsilea vestita
BMC Cell Biologyvolume 17, Article number: 29 (2016)
Spermatogenesis in the semi-aquatic fern, Marsilea vestita, is a rapid, synchronous process that is initiated when dry microspores are placed in water. Development is post-transcriptionally driven and can be divided into two phases. The first phase consists of nine mitotic division cycles that produce 7 sterile cells and 32 spermatids. During the second phase, each spermatid differentiates into a corkscrew-shaped motile spermatozoid with ~140 cilia.
Analysis of the transcriptome from the male gametophyte of Marsilea revealed that one kinesin-2 (MvKinesin-2) and two kinesin-9 s (MvKinesin-9A and MvKinesin-9B) are present during spermatid differentiation and ciliogenesis. RNAi knockdowns show that MvKinesin-2 is required for mitosis and cytokinesis in spermatogenous cells. Without MvKinesin-2, most spermatozoids contain two or more coiled microtubule ribbons with attached cilia and very large cell bodies. MvKinesin-9A is required for the correct placement of basal bodies along the organelle coil. Knockdowns of MvKinesin-9A have basal bodies and cilia that are irregularly positioned. Spermatozoid swimming behavior in MvKinesin-2 and -9A knockdowns is altered because of defects in axonemal placement or ciliogenesis. MvKinesin-2 knockdowns only quiver in place while MvKinesin-9A knockdowns swim erratically compared to controls. In contrast, spermatozoids produced after the silencing of MvKinesin-9B exhibit normal morphology and swimming behavior, though development is slower than normal for these gametes.
Our results show that MvKinesin-2 and MvKinesin-9A are required for ciliogenesis and motility in the Marsilea male gametophyte; however, these kinesins display atypical roles during these processes. MvKinesin-2 is required for cytokinesis, a role not typically associated with this protein, as well as for ciliogenesis during rapid development and MvKinesin-9A is needed for the correct orientation of basal bodies. Our results are the first to investigate the kinesin-linked mechanisms that regulate ciliogenesis in a land plant.
Spermatogenesis in the semi-aquatic water fern, Marsilea vestita, is a rapid, synchronous process that produces motile gametes in only 11 h . Similar to other rapidly developing systems, spermatogenesis in male gametophytes of Marsilea is controlled at a post-transcriptional level. Virtually all the RNA required for development is present in the dry spore at the time of its rehydration and almost no additional transcription is needed for spermiogenesis to reach completion [2–4]. Rapid development is dependent on the unmasking of pre-mRNAs , which are stored in dry spores within nuclear speckles , and then, the processing of these pre-mRNAs for translation . Spermatogenesis begins when microspores are exposed to water. Shortly after hydration, the gametophyte initiates a developmental program that culminates with the production of 32 corkscrew-shaped spermatozoids, each with ~140 cilia [8–11]. Development can be divided into two phases. The first phase consists of a series of nine mitotic division cycles that produce 32 spermatids that are surrounded by seven sterile cells. All divisions are complete approximately 5 h after microspore hydration. During the second phase, only the spermatids undergo drastic morphological changes and differentiate into motile, corkscrew-shaped spermatozoids [4, 12].
This unusual shaping of the gamete is achieved through the elongation and coiling of the nucleus and mitochondria along a coiled ribbon of crosslinked microtubules. Basal bodies, already formed de novo, are placed in two rows at regular intervals along the dorsal face of the microtubule ribbon to become the sites of ciliogenesis [8, 11]. At first, basal bodies are oriented so cilia diverge away from each other and are parallel to the plasma membrane of each spermatid. Near the end of spermiogenesis, the basal bodies rotate 90° so that the ciliary axonemes protrude vertically from the microtubule ribbon and nuclear coil . At 9.5 h of development, an extension of cytoplasm begins to grow around the anterior end of each spermatid and eventually fuses together to surround each cell. This creates an internal, but extracellular space that contains the microtubule ribbon and organelle coil plus all of the cilia. Upon release from the microspore, each spermatozoid breaks free from the surrounding cytoplasmic extension and leaves behind a thin vesicle-like structure . The ciliary axonemes have the typical 9 + 2 architecture  found in motile organisms and spermatozoids are able to swim towards the megaspore for fertilization.
We are interested in the processes that regulate spermatid differentiation and ciliogenesis during male gametophyte development in Marsilea. Proteins important for ciliary assembly and function are moved to the distal ends of forming axonemes by intraflagellar transport (IFT) involving members of the kinesin-2 family. Heterotrimeric kinesin-2 consists of a kinesin-2α, a kinesin-2β, and a kinesin associated protein (KAP) that regulates cargo binding . In Chlamydomonas, heterotrimeric kinesin-2 is necessary for IFT [15–18]. Kinesin-2 can also function as a homodimer of two kinesin-2γ subunits. Kinesin-2γ, also known as OSM-3 or KIF17, has a distinct role in assembling sensory cilia that is separate from heterotrimeric kinesin-2 [19–24]. Although the most common mechanism for axonemal assembly is dependent on IFT function, Plasmodium falciparum and the sperm flagella of Drosophila are able to build motile cilia and flagella using IFT-independent mechanisms. In this case, cilia are assembled in the cytoplasm [25–27]. Ciliogenesis in Marsilea does not occur in this way; instead, cilia are assembled in growing membrane extensions from basal bodies that are positioned along a microtubule and organelle coil. Therefore, we anticipated that ciliogenesis in Marsilea is reliant on IFT-dependent mechanisms of transport and assembly. In support of this idea, we found a variety of transcripts encoding IFT-associated proteins (kinesin-2, dynein-1b, IFT-A, and IFT-B subcomplex proteins) in the assembled transcriptome from the male gametophyte of this organism .
Phylogenetic analysis of kinesin-9 shows the existence of two subfamilies, kinesin-9A and kinesin-9B. Kinesin-9A includes KIF9 and kinesin-like protein 1 (KLP1), while kinesin-9B includes the KIF6 protein [28, 29]. Kinesin-9A is localized to the central pair of microtubules in the Chlamydomonas axoneme  and is necessary for motility in both Chlamydomonas  and Trypanosoma brucei . Kinesin-9A appears to be important for motility by regulating the activity of flagellar dynein  and by interacting with Hydin, which is necessary for motility in algae , trypanosomes , and mice . Less is known about the function of kinesin-9B; however, in T. brucei, KIF9B localizes to basal bodies and the flagellum where it is necessary for the construction of the paraflagellar rod .
Although the processes that regulate the structure and function of cilia are highly conserved, not all organisms make ciliated cells. This is most pronounced during the evolution and adaption of land plants from green algae. Conifers and angiosperms never make any ciliated cells, while lower plants known as embryophytes, including the ferns, mosses, liverworts, and certain members of the gymnosperms (e.g., Ginkgo biloba and the cycads) only produce cilia in their male gametes. In conjunction with the loss of cilia in the ‘higher’ organisms, the reduction , or complete absence, of proteins important for ciliogenesis and motility is also observed. For example, comparative analyses of the kinesin family has shown that kinesin-2, kinesin-9, and the more recently identified, kinesin-‘orphan’ III (also referred to as kinesin-16) and kinesin-17, are only found in organisms that are ciliated at some point during the life cycle [28, 36]. Members of these kinesin families can be found in ciliated plants such as Chlamydomonas, the moss Physcomitrella, and the water fern Marsilea, except kinesin-17, which is only present in Chlamydomonas [28, 36–38]. The genome of Arabidopsis, a flowering plant that is never makes ciliated cells, contains none of these kinesins . Due to the reduction of cilia and IFT proteins in land plants, the majority of research on ciliogenesis in plants has been conducted in Chlamydomonas. It is somewhat ironic that the conserved 9 + 2 microtubule organization of motile axonemes was first observed in the spermatozoid of a fern .
Here, we are using spermatogenesis in Marsilea as a model for ciliogenesis. This provides us the unique ability to study the construction, organization and motility of ciliary apparatus produced in gametes of a land plant. Moreover, this gametohyte provides insights on the mechanisms that evolved to regulate de novo ciliogenesis in specialized cells of otherwise nonmotile organisms. Transcripts that encode members of the kinesin-2 and kinesin-9 families were selected as targets for our studies on the regulation of ciliogenesis in spermatids of Marsilea. mRNAs that encode kinesin-2 and kinesin-9 increase in abundance during the stage of development associated with spermatid shaping and ciliogenesis . These late rises in transcript abundance led us to suspect that kinesin-2 and kinesin-9 may play critical roles for spermatid differentiation and ciliogenesis. Unlike many other systems, the male gametophyte of Marsilea has only one kinesin-2 that is apparent in its transcriptome. It is most similar to the kinesin-2 found in Physcomitrella and is divergent from the typical heterotrimeric kinesin-2 associated with IFT. In this study, we show that MvKinesin-2 is required for two separate events in this gametophyte. It is necessary for cytokinesis in spermatogenous cells and is also important for regulating the length of cilia during later phases of spermatid maturation. The Marsilea gametophyte has two transcripts that encode members of the kinesin-9 family; one is similar to the kinesin-9A and the other is most like kinesin-9B. In the gametophyte, we show that MvKinesin-9A is involved in the proper positioning of basal bodies that are required for ciliogenesis and it is necessary for motility. MvKinesin-9B is needed for the timely differentiation of motile spermatozoids in the rapidly developing gametophyte.
Characterization of the kinesin-2 and kinesin-9 family in Marsilea
After conducting RNAseq with poly(A+)-RNA isolates obtained from male gametophytes of Marsilea during different intervals of development, we were able to assemble a transcriptome de novo , and calculate changes in transcript abundance that occur during different phases of spermatogenesis. We searched our combined reference transcriptome for members of the kinesin-2 and the kinesin-9 families using Physcomitrella and Chlamydomonas kinesins-2 and kinesin-9 sequences (Additional file 1). We used these kinesins as the base of our search because the entire kinesin family in Physcomitrella has recently been identified  and because Physcomitrella and Chlamydomonas are two of the more closely related ciliated species to Marsilea with an annotated kinesin family. We then verified our results by comparing them to the 56 kinesin-like transcripts that were previously identified in our transcriptome . Our search identified one candidate kinesin-2 (MvKinesin-2) transcript and two candidate kinesin-9 (MvKinesin-9A and MvKinesin-9B) transcripts (Additional file 1).
Unlike other species, there is only one kinesin-2 present in Physcomitrella  and Marsilea. The translated kinesin-2 sequence from Marsilea (MvKinesin-2) shows that MvKinesin-2 recognizes members of the kinesin-2 family from other ciliated eukaryotes. Reciprocal Blastp analysis shows that MvKinesin-2 recognizes CrFLA8 (XP_001697037) and CrFLA10 (XP_001701510) with an e-value of 0.0. Phylogenetic analysis of the motor domain of established kinesin-2 s in plants and animals (Additional file 2) demonstrates that the single kinesin-2 motor in Marsilea cannot be placed in the standard kinesin-2α (-2A, shown in red), -2β (-2B, shown in orange), or -2γ (-2C, shown in green) subgroups (Fig. 1a). This is similar to our findings for both FLA10 and FLA8 from Chlamydomonas and for kinesin-2 from Physcomitrella. MvKinesin-2 is most similar to the kinesin-2 in Physcomitrella (Fig. 1a). This is not surprising since the only ciliated cells in Marsilea and Physcomitrella are their male gametes. It is unclear whether this kinesin functions as a homodimer, a heterotrimer with yet unidentified partners, or as a unique single protein.
The male gametophyte of Marsilea has two transcripts that encode kinesin-9-like proteins. Phylogenetic analysis of the motor domain of established kinesin-9 s in plants and animals (Fig. 1b, Additional file 2) demonstrates that the kinesin-9 s in Marsilea can be separated into the kinesin-9A (red) and kinesin-9B (blue) subfamilies (Fig. 1b). This provides further evidence for the existence of two distinct kinesin-9 subfamilies. Reciprocal Blastp analysis shows that translated kinesin-9A (MvKinesin9A) and kinesin-9B (MvKinesin9B) sequences in Marsilea recognize other members of the kinesin-9 family. MvKinesin9A is similar to CrKLP1 (XP_001701617) and TbKIF9A (XP_846252) with e-values of 4e-112 and 8e-50, respectively. MvKinesin9B recognizes CrKIF9B (Cre01.g036800.t1.1) and TbKIF9B (XP_846346) with e-values of 0.0 and 6e-89, respectively.
MvKinesin-2, MvKinesin9A, and MvKinesin9B have typical kinesin motor and ATP-binding domain architectures including the P-loop, switch I, and switch II motifs (Additional file 3). Outside of this highly conserved motor domain, MvKineisn-2, MvKinesin9A, and MvKinesin9B show a similarity to counterpart kinesins in Physcomitrella and Chlamydomonas, however, large insertions and gaps exist in the non-motor regions for each of these alignments (Additional file 3).
RT-PCR using poly(A+)-RNA isolated at 1–2 h, 3–5 h, and 6–8 h of development shows that MvKinesin-2, MvKinesin-9A, and MvKinesin-9B transcripts increase in abundance from the 1–2 h to the 6–8 h time intervals of gametophyte development (Fig. 1c). This pattern of transcript availability and abundance suggests that kinesin-2 and kinesin-9 transcripts are unmasked, spliced and polyadenylated prior to the time interval associated with ciliogenesis during gametophyte development . It is unclear why kinesin-9B transcripts decrease during 3–5 h time interval. We have seen in our transcriptome analysis that the middle time point of development (3–5 h post hydration) appears to be a transition time between the early stage of development, marked by cell division, and the later portion of development, which is dedicated to spermatid differentiation and ciliogenesis. GO terms that are enriched during this transition time includea variety of nucleases, and proteasome and ubiquitin ligase components. It is possible that kinesin-9B transcripts and proteins are destroyed during this time period of development, and then additional kinesin-9B transcripts are unmasked, processed and translated for subsequent kinesin-9B functions during ciliogenesis.
Kinesin-2 and kinesin-9 are involved in different aspects of spermatogenesis
To determine the function of these kinesins during spermatogenesis, we performed RNAi knockdown experiments, treating gametophytes at the time of spore hydration with dsRNA. To do this we transcribed and constructed dsRNA from unique regions of each kinesin (Fig. 2a–b). It was important to use unique regions to make dsRNA in order to prevent broad silencing of off-target sequences and to ensure that only one kinesin was being silenced. The effectiveness of RNAi was measured using RT-PCR for each transcript after knockdown at 8 h of development. In each case, the presence of the transcript could not be detected after its individual knockdown, while the presence of other kinesin and centrin mRNAs remained unchanged (Fig. 2c–f). Microspores were grown for 8 h, fixed, embedded in methacrylate, and sectioned. The sections were then stained with toluidine blue (TBO) and examined with brightfield microscopy to observe broad morphological changes in development that might have resulted from the knockdown.
After 8 h of development in untreated gametophytes, all cell divisions are complete and spermatids can be clearly distinguished from jacket cells by their size, location within the microspore, shape, and staining pattern. Jacket cells are located towards the outside of the microspore near the microspore wall and contain starch filled plastids. Jacket cells do not contain very much mRNA or proteins and therefore do not stain intensely with TBO. During this stage of development, each spermatid begins to differentiate into a motile spermatid by forming a coiled nucleus and microtubule ribbon, creating a slightly boomerang shaped cell (Fig. 2g) .
Spermatogenous cells can clearly be distinguished from jacket cells in MvKinesin-2 knockdowns; however, in the majority of gametophytes (79/107; 74 %) spermatogenous cells are larger than normal and do not have a consistent size and shape compared to controls. This suggests that one or more cell division cycles have been skipped or failed to reach completion (Fig. 2h) . In a minority (28/107; 26 %) of these knockdowns, spermatogenous cells are the correct size and shape and appear to be morphologically normal (Fig. 2i).
In MvKinesin-9A knockdowns, spermatogenous cells were the same size and shape as controls, suggesting normal progression through all of the cell division cycles (Fig. 2j). However, a distinct phenocopy was observed with knockdowns of kinesin-9B. In these gametophytes, spermatogenous cells appear larger than controls and the cells are rounded. This suggests that nuclear elongation failed to occur. Thus, while cell divisions have proceeded normally, differentiation is aberrant (Fig. 2k).
Spermatid differentiation is incomplete after knockdown of MvKinesin-2, MvKinesin-9A, and MvKinesin-9B
In order to assess altered patterns of spermatid maturation, we fixed and sectioned gametophytes that had undergone silencing of MvKinesin-2, MvKinesin-9A, or MvKinesin-9B. We stained the methacrylate sections of fixed gametophytes with DAPI to observe nuclear elongation and labeled them with anti-centrin antibodies to observe the presence and distribution of basal bodies in the spermatids . During normal spermatid differentiation, the basal bodies become situated at regular intervals along the microtubule ribbon and the nuclear coil where they serve as templates for the growth of ciliary axonemes (Fig. 3a). By observing the distribution centrin protein in kinesin knockdowns, we are able to make conclusions about the formation and localization of basal bodies in these cells [3, 5, 7].
In each knockdown, the appearance and placement of basal bodies is altered. At 8 h of development, anti-centrin staining in the majority (85/114; 75 %) of MvKinesin-2 knockdowns is diffuse and there are few, if any, aggregates of centrin that resemble typical basal body staining in spermatogenous cells. Some larger centrin aggregates that possibly resemble the blepharoplast can be found in spermatogenous cells (Fig. 3b). DAPI staining in these cells is also distinctly different from controls and it is apparent that the nuclear coil is not properly formed at this stage (Fig. 3b). Just as a minority of microspores treated with dsRNA complementary to MvKinesin-2 appeared to proceed through development normally when observed with TBO, a similar percentage (29/114; 25 %) also have normal centrin and DAPI staining (Fig. 3c). These microspores were either not affected by RNAi or it is possible that there are two separate phenocopies of MvKinesin-2 knockdown during spermatogenesis in Marsilea.
In MvKinesin-9A knockdowns, we observed aggregates of centrin staining that resemble basal bodies. The presumptive basal bodies become localized in areas of the spermatids where they normally are not seen. Later in development, these basal bodies become positioned at irregular intervals around the cell periphery; they are not restricted to the anterior (coiled) portion of the spermatid, as they would be in normal cells. DAPI staining in these cells shows a normally elongated nucleus (Fig. 3d). Therefore, the elongation stage of spermatid differentiation is occurring, but there are anomalies in the localization of basal bodies in these cells.
MvKinesin-9B knockdowns show large aggregates of centrin protein (Fig. 3e) that resemble blepharoplast particles. The blepharoplast is a cytoplasmic particle that forms during the last spermatogenous cell division. During the last mitotic division, the blepharoplast functions like a centrosome at the spindle pole, though it lacks any organized centrioles. As the spermatids are being formed, the blepharoplast disappears and then reforms to serve as a site for de novo basal body assembly . In this knockdown, no further maturation of the blepharoplast (for basal body formation) was observed. Spermatogenous cells in MvKinesin-9B knockdowns exhibit nuclei that are round in shape (Fig. 3e), while normal spermatids undergo nuclear elongation and coiling . This pattern of centrin immunolabeling and DAPI nuclear staining is characteristic of arrest at a stage of development preceding basal body formation in spermatids, and suggests that either spermatid maturation has stopped early in gamete differentiation or that the differentiation process has been slowed substantially.
Silencing of MvKinesin-2, MvKinesin-9A, and MvKinesin-9B cause defects in ciliogenesis and motility
To determine whether the abnormalities observed with these kinesin knockdowns during spermatogenesis had any impact on ciliogenesis or on motility in released gametes, we observed the kinesin knockdowns at 11 h of development. At 11 h, the normal process of spermatogenesis reaches completion, and 32 motile spermatozoids, each with ~140 cilia, break free from their enclosing microspore walls. First, the opaque outer exine wall of the microspore begins to thin, revealing a much thinner and translucent intine wall (Fig. 4a–1). Then two clusters of 16 spermatozoids emerge from the microspore (Fig. 4a–2). The gametes are each contained within a thin wall-like structure built from an extension of the cytoplasm that originally surrounded the spermatid. Each motile spermatozoid rapidly spins in place, presumably in order to break free from its surrounding wall (Additional file 4). The shape of each spermatozoid resembles a corkscrew with cilia extending vertically from the edges of the microtubule and nuclear coil (Fig. 4a–3). Once released, the spermatozoids swim in a shallow helical path, rotating as they swim. By borrowing terms from aerodynamics, we have termed this rotational behavior “rolling” (Fig. 4b, Additional file 5). This pattern is repeated until the spermatozoid dies or reaches the egg cell in the female gametophyte (produced by the megaspore) for fertilization.
After the silencing of MvKinesin-2, for the majority of gametophytes (20/27; 74 %), the emergence of spermatozoids from the spore wall differed considerably from controls. Although the exine wall thinned as usual (Fig. 4c–1), the spermatozoids exited from the microspore as individuals, rather than as groups of 16 cells (Fig. 4c–2, Additional file 6). The spermatozoids that emerged this way were much larger than controls and often contained two or more coiled microtubule ribbons with attached cilia. Also, these cells lacked the characteristic extension of cytoplasm that unusually surrounds each spermatid (Fig. 4c–3). Due to their large size and multiple coils (Fig. 4c–3), these cells were unable to swim and could only quiver in place or yaw along their y-axes (Fig. 4b, Additional file 7). We call this phenocopy ‘Monster’ because each cell has multiple coils and cilia. It is likely that these cells failed to complete cytokinesis during the last spermatogenous cell division, but nevertheless, they were able to form microtubule coils in association with elongated nuclei and subsequently, each coil was associated with it complement of cilia.
In a minority (7/27; 26 %) of MvKinesin-2 knockdowns, the spermatids emerged from the microspore as normal; however, these spermatozoids exhibited abnormally long ciliary axonemes (Fig. 4d–3). The spermatozoids with long cilia were unable to swim in a shallow helical path. These cells could only roll and essentially, swim in place (Fig. 4b, Additional file 8). We refer to this phenocopy as ‘Rapunzel’ for its long cilia. Monster and Rapunzel are two different phenocopies caused by MvKinesin-2 knockdowns. It is possible that MvKinesin-2 participates in two distinct processes during spermatogenesis; the first occurs during cytokinesis leading to spermatid formation, and the second occurs during ciliary axonemal formation and affects ciliary length control. An alternate explanation is that the Monster and Rapunzel phenocopies are two manifestations of the same biological event caused by knockdown of kinesin-2. In this scenario, the Monster phenocopy is manifested earlier during spermatogenesis and Rapunzel later, in ciliogenesis. Because Monster spermatids were released from the microspores as groups of cells and Rapunzel spermatids as individuals, we suspect that each of these phenocopies occurs in different microspores and that all the spermatozoids emerging from one microspore display the same defects in development.
The thinning of the exine wall and the emergence of the spermatids from the microspore is normal in kinesin-9A knockdowns (Fig. 4e–1, 2). Once the spermatozoids emerge (Fig. 4e–3), their swimming behavior is affected by the knockdown (Additional file 9). These gametes roll, yaw, and pitch in place (Fig. 4b), effectively spinning without vectorial movement. Directional swimming for these cells occurs infrequently; they travel only short linear distances before switching back to spinning behavior (Additional file 9). Based on the movies, it is difficult to determine why the kinesin-9A knockdowns are unable to swim properly.
We have named our kinesin-9B knockdowns ‘Late Bloomer’ because the cells emerge slowly after the seemingly normal thinning of the exine (Fig. 4f–1, 2). The cells finally emerge ~16 h after microspore hydration, in contrast to 11 h for controls, and for the MvKinesin-2, and MvKinesin-9A knockdowns. The Late Bloomer spermatozoids that emerge are morphologically normal and have the same shape as control gametes (Fig. 4f–3). These cells also exhibit a helical swimming behavior similar to that of normal cells (Additional file 10).
MvKinesin-9A knockdowns have irregularly positioned cilia
To determine why kinesin-9A knockdowns were unable to swim normally when compared to controls, we used a variety of methods to observe the arrangement of cilia on the fully developed spermatozoids. Firstly, we fixed fully emerged spermatozoids quickly in 2 % PFA on the surface of a microscope slide. We then observed these fixed cells with DIC microscopy. In controls, the corkscrew shape of the spermatozoid cell body, the nuclear coil, and the ciliary axonemes can clearly be seen (Fig. 5a). Normal spermatozoids have the cilia placed in two regular rows along the dorsal flank of the anterior, coiled cell body [11, 43]. Kinesin-9A knockdowns show cilia that are irregularly positioned on the spermatozoid cell body and are not restricted to the anterior portion of the cell (Fig. 5c). We named this phenocopy ‘Porcupine’ for the cilia that emerge from all points of the cell. None of the other kinesin silencing treatments exhibited irregularly placed cilia similar to Porcupine in spite of the fact that MvKinesin-2 knockdowns (Monster) generally have cell bodies that are misshapen and larger than normal (Fig. 5b) and Kinesin-9B knockdowns (Late Bloomer) appear to be structurally normal even though they emerge at 16 h instead of 11 h after spore hydration (Fig. 5d).
Cilia are irregularly positioned in MvKinesin-9A knockdowns due to the improper positioning of basal bodies in these cells
The abnormal distribution of the cilia and aberrant swimming patterns observed in the Porcupine phenocopy appears to be caused by the improper positioning of basal bodies along the microtubule ribbon in these cells (Fig. 3d). In order to test this hypothesis, we sectioned fixed microspores at 10.5 h after spore hydration and stained the sections with TBO, DAPI, and an anti-centrin antibody. Sections were visualized with phase contrast and fluorescence microscopy.
At 10.5 h of development, the normal spermatozoids are fully developed, but have yet to emerge entirely from the microspore wall. Using TBO staining at this time point, the organelle coil can easily be observed in each nearly-mature spermatozoid. Cilia emerging along this coil are also visible. The cilia are wrapped around the cell body coil and are seen in transverse-section in some of the gametes as a white area adjacent to the organelle coil. We refer to this area as the ciliary band (Fig. 6a). Anti-centrin staining in control sections form aggregates that resemble basal bodies. In control spermatids, the basal bodies are restricted to areas on the dorsal side of the coiled nucleus (Fig. 6b). Overlaying the fluorescence images with phase contrast images of the same section, it appears that the centrin aggregates associate uniformly along the edges of the organelle coil and are often within areas of the ciliary band (Fig. 6c).
In comparison with untreated gametophytes, in MvKinesin-9A knockdowns stained with TBO, it is difficult to visualize the ciliary band. In cells that possess a visible ciliary band, the band is irregularly shaped compared to controls (Fig. 6g). This suggests anomalies in the orientation and localization of the basal bodies in the spermatid. Centrin staining confirms this, showing aggregates that resemble basal bodies present throughout the cytoplasm of the nearly-mature spermatozoids. Aggregates of centrin are present at the dorsal face of the spermatozoid along the edges of the organelle coil as in controls, but they are also abundant along the posterior portion of many spermatozoids (Fig. 6h, i). The failure of these cells to position the axonemes correctly during differentiation is a likely explanation for the Porcupine phenocopy and the odd swimming patterns observed in these knockdowns.
In contrast to this phenocopy, MvKinesin-2 knockdowns visualized with TBO have large, irregularly shaped cells. We suspect this is the consequence of arrested, aberrant or blocked cell division cycles. Coils are visible only in some of the maturing spermatids (Fig. 6d). Anti-centrin antibody staining in these knockdowns appears as a diffuse cloud of fluorescence in spermatogenous cells. Aggregates that resemble basal bodies are present in some cells; however, they do not appear localized in any reliable pattern that resembles basal body distributions in normal cells (Fig. 6e, f). This pattern is reminiscent of the Monster phenocopy where the released gametes are larger than normal and possess multiple coils containing cilia.
In kinesin-9B knockdowns, the spermatogenous cells resemble to control cells, though the coil and cilia are not as apparent (Fig. 6j). Anti-centrin antibody staining forms aggregates that resemble basal bodies and are positioned at regular intervals along the anterior portion of the nuclear coil (Fig. 6k, i). These cells do not look like they have developed to the same stage as the controls though clearly, they have proceeded into the later stages of differentiation. This provides further evidence that knockdown of kinesin-9B (Late Bloomer) causes a general slowing of differentiation.
MvKinesin-2 has atypical functions in cytokinesis and ciliogenesis during male gametophyte development in Marsilea
Phylogenetic analysis of MvKinesin-2 shows that this kinesin is atypical when compared to most other members of the kinesin-2 family. Firstly, MvKinesin-2 cannot be placed in the usual kinesin-2α, β, or γ subgroups, though it most resembles the kinesin-2 isoform found in Physcomitrella. This was unsurprising since both FLA10 and FLA8 in Chlamydomonas  and the kinesin-2 in Physcomitrella also group independently (Fig. 1a). Similarly, many non-animal kinesin-2 s also cannot be placed in the established kinesin-2α, β, or γ subgroups (Fig. 1a). Secondly, we found only a single kinesin-2 transcript in our transcriptome. This was surprising because kinesin-2 is best known to function in heterotrimeric complexes. A single kinesin-2 motor can also be found in Physcomitrella and in certain protists [35–37]. Although the reason is unclear, organisms that are only ciliated for a short time during the life cycle show a general reduction in cilia-associated proteins . The single kinesin-2 motor found in Marsilea may function as a homodimer, similarly to C. elegans OSM-3, though evidence to support this existence of homodimeric kinesin-2 outside of sensory cilia is limited . From this analysis, we were unable to discern whether in gametophytes MvKinesin-2 functions as a homodimer, with yet unidentified partners, or if it acts alone. It is possible that MvKinesin-2 and single kinesin-2 proteins found in other organisms function separately from the more well studied kinesin-2 motors.
Functional silencing assays showed that MvKinesin-2 is required for cytokinesis during the cell divisions that ultimately produce 32 spermatids in the gametophyte (Monster phenocopy) and for regulating the length of cilia (Rapunzel phenocopy). However, we cannot eliminate the possibility that the Monster and Rapunzel phenocopies represent two different manifestations of the same biological event where the role for MvKinesin-2 in cytokinesis precedes its involvement in ciliogenesis. Monster spermatids formed can still produce motile cilia, but these abnormally large gametes contain multiple coils with attached cilia and are unable to swim and can only quiver in place (Additional file 7). In addition, these cells lack the extension of cytoplasm that usually surrounds the anterior portion of each spermatid (Fig. 4c–3). Although not typically associated with kinesin-2, a role for this motor in cytokinesis has been reported [45–47]. How exactly kinesin-2 participates in cell division has yet to be established as mutants are able to pass through mitosis without visible consequence [48–50]. We suspect that the cytokinesis defect phenocopy observed in Marsilea is related to the function of kinesin-2 in the intracellular transport of Golgi-derived vesicles and membrane dynamics [45, 51–53]. During plant cell cytokinesis, the coordinated efforts of a variety of motor proteins are necessary to deliver Golgi-derived vesicles along phragmoplast microtubules to build the cell plate . Heretofore, kinesin-2 has not been implicated in this process and a recent study in Physcomitrella showed that kinesin-2 was not expressed in caulonemal cells undergoing mitosis . Since both Marsilea and Physcomitrella only produce ciliated cells in their gametophytes, it is likely that the expression of kinesin-2 is restricted in land plants to the formation of male gametes. It is possible that this motor has evolved or retained a role in cytokinesis and membrane dynamics specifically for spermatogenesis.
In a minority of cells, MvKinesin-2 knockdowns produced abnormally long ciliary axonemes (Rapunzel phenocopy) while maintaining normal cell size and shape. This result suggests a function for MvKinesin-2 during ciliogenesis that is distinct from its role in cytokinesis. Rapunzel spermatozoids were unable to swim vectorially and could only rotate in place. Heterotrimeric kinesin-2 is responsible for anterograde transport during ciliogenesis. Mutations in this kinesin isoform typically result in no change in ciliary length or the absence of full-length cilia due to suppressed IFT activity [15, 17, 21, 22, 56]. Therefore, the production of long cilia was unexpected. One possible explanation is that other proteins are responsible for modulating IFT dynamics and ciliary length during male gametophyte development in Marsilea. Evidence in support of this hypothesis can be found from the cephalic male (CEM) cilia of C. elegans. Here loss of KLP-11 (a member of the heterotrimeric kinesin-2 complex) results in long cilia due to the complex interaction of KLP-11, OSM-3, and the kinesin-3 motor, KLP-6, during IFT . In Chlamydomonas kinesin-13 is important for regulating axoneme assembly and length [58, 59] and KCBP, a plant-specific member of the kinesin-14 family, is necessary for the function of flagella . It is possible that kinesin-13, KCBP, other plant-specific, or orphan kinesins are responsible for modulating IFT dynamics during ciliogenesis in Marsilea. Additional functional studies on these kinesins are necessary to investigate this possibility.
Kinesin-9 is required to orient basal bodies along the microtubule ribbon and for spermatid differentiation during spermiogenesis
Our analysis of the kinesin-9 family echoes previous findings that there are two kinesin-9 subgroups [28, 29]. In the Marsilea gametophyte transcriptome, a transcript that encodes one member of each kinesin-9 subgroup is present. With silencing treatments, we found that kinesin-9A is necessary for the organization and placement of basal bodies that become the sites of ciliogenesis during spermatid maturation. Since the basal bodies are not properly placed along the microtubule ribbon, and, as a consequence, the cilia are misplaced in the cell body, the spermatozoids cannot maintain normal swimming patterns and they are only able to spin and flip in place. In both Chlamydomonas and T. brucei, kinesin-9A is also important for motility; however, this is achieved through the interaction of kinesin-9A with central pair microtubules in Chlamydomonas [29, 31] and not through regulating the orientation of basal bodies. The typical phenotype observed in kinesin-9A mutants of a reduction in ciliary beating [29, 31] was not observed in Marsilea kinesin-9A knockdowns. Instead, cilia were unable to maintain normal swimming and the gamete moved in a highly disorganized pattern (Additional file 9).
We also found that knockdowns of kinesin-9B produced a general slowing of differentiation as fully mature and functional spermatozoids were released 5 h later than controls. Early stages of development were not altered after kinesin-9B knockdown. This result implies that kinesin-9B is needed for the differentiation of spermatids into motile cells. The specific role for kinesin-9B in differentiation is not known, but we suspect that there is sufficient redundancy among kinesin isoforms to allow development to proceed, albeit slowly, in the absence of this motor protein. Functional studies on kinesin-9B in other systems have begun to appear in the literature. In T. brucei, kinesin-9B is responsible for building the paraflagellar rod (PFR) . The PFR is unique to kinetoplastic protozoa like T. brucei so, in a fashion similar to the spermatozoid of Marsilea, it is not surprising in other eukaryotes that kinesin-9B would perform functions beyond PFR assembly during ciliogenesis.
Among the 56 kinesin family members whose mRNAs are detectable in the male gametophyte of Marsilea , we have found that kinesin-2, kinesin-9A, and kinesin-9B are each required for specific events during spermatid formation or maturation and for ciliary assembly or positioning during rapid development. mRNAs that encode these kinesins increase in abundance during the developmental stage associated with spermatid differentiation and ciliogenesis. RNAi knockdowns of each of these kinesins show distinct and separate phenocopies during spermatogenesis and spermiogenesis. Functional analyses showed that kinesin-2 is involved in cell division (Monster phenocopy) and in regulating the length of motile cilia (Rapunzel phenocopy). Knockdown experiments showed that kinesin-9A is required for the proper positioning of basal bodies (Porcupine phenocopy) during spermatid differentiation leading to spermatozoids that are unable to power directional swimming. In contrast, kinesin-9B (Late Bloomer phenocopy) does not appear to be required for the construction of motile spermatozoids in Marsilea.
Heterotrimeric kinesin-2 is the main motor associated with anterograde IFT in motile cilia and mutations typically produce cells that lack cilia due to resultant problems with IFT [15–18]. Therefore, the involvement of kinesin-2 in cytokinesis and in regulating ciliary length in the Marsilea male gametophyte is atypical. Although best known for its functions during ciliogenesis, kinesin-2 has also been implicated in mitosis and in intracellular membrane trafficking in a number of cell types [45, 51–53] thus potentially explaining the Monster phenocopy observed during gametophyte development. Also kinesin-2 is not the only motor implicated in anterograde IFT. In the C. elegans cephalic male cilia (CEM) KLP-6, a member of the kinesin-3 family, works together with kinesin-2 during IFT to regulate ciliary length . It is possible that the Marsilea kinesin-2 interacts with other motors during ciliogenesis and in IFT thus explaining the Rapunzel phenocopy. The Porcupine phenocopy observed in Marsilea after knockdown of kinesin-9A was also different than currently established roles for this motor. The typical phenotype observed in kinesin-9A mutants of a reduction in ciliary beating due to the interaction of kinesin-9A with central pair microtubules [29–31] In Marsilea, we attribute disorganized ciliary beating and the inability to maintain normal swimming patterns on the mislocalization of basal bodies in mature spermatozoids. The reason for these contradictions is not currently understood; however, significant differences in the structure of axoneme in fern spermatozoids including those from Marsilea [61–63] may be to blame.
Sporophytes of Marsilea vestita
Our original collection of sporocarps of M. vestita was obtained from Dr. Peter K. Hepler (University of Massachusetts), who had initially collected sporocarps in the 1970s from Lake Lagunita, which is located on the Stanford University campus in California. Dry sporocarps remain viable for extended periods [1, 42]. In our annual effort to raise new sporophytes, approximately 20–30 dry sporocarps were scarified and immersed in water in an aerated aquarium to activate development in microspores and megaspores. After 1 or 2 weeks, small sporophytes were transferred to water-saturated soil in growth flats and raised indoors until they were approximately 1–2 cm in height. At this stage they were transferred to pots and placed into artificial ponds located at the University of Maryland Research Greenhouse Complex, where they were allowed to grow for 6–8 months as semi-submersed plants . By this time, the plants were actively making new sporocarps. Thereafter, we allowed the ponds to dry down gradually over the course of 2–3 months and during this period, the sporophytes turned yellow, and the spores contained in sporangia within the sporocarps underwent their natural desiccation process. Sporocarps were then harvested from the dried plant material and stored at room temperature. Successive plantings over the years allowed us to create a stockpile of sporocarps that contain viable, desiccated spores. A set of sporophyte voucher specimens is being prepared for storage at the Norton-Brown Herbarium at the University of Maryland (College Park, MD).
Marsilea vestita microspores were isolated from dry sporocarps as previously described [1, 42]. Four milligrams of microspores were then cultured in 1 ml commercial spring water (Deer Park) for 1 h at 20 °C with rotational shaking. Gametophytes were then transferred to 50 ml flasks containing a total of 10 ml spring water and cultured with rotational shaking at 20 °C [42, 64].
Identifying Kinesin-2 and Kinesin-9 Sequences in Marsilea
Kinesin-2 and kinesin-9 sequences were identified by searching the entire gametophyte transcriptome with kinesin-2 and kinesin-9 sequences from Physcomitrella (Phypa_425592, Phypa_425498, Phypa_458410, Phypa_428375) and Chlamydomonas (CrFLA8, CrFLA10, CrKLP1, CrKIF9b). Stand-alone BLAST was used for the analysis. A maximum e-value of 1E-100 was used to yield the best possible matches and avoid identifying kinesins unrelated to kinesin-2 and -9 in the search (Additional file 1).
MvKinesin-2 and MvKinesin-9 Phylogenetics
Sequences for previously analyzed kinesin-2 and kinesin-9 family members were obtained using NCBI , Phytozome , and Cosmoss  databases. The motor domain from each sequence was identified using the NCBI Conserved Domain Database  and isolated for analysis. Phylogenetic trees for kinesin-2 and kinesin-9 were built by importing the sequence files into MEGA  and generating a multiple sequence alignment using the ClustalW default parameters. Maximum likelihood (ML) phylogenetic trees were built using the GARLI (Genetic Algorithm for Rapid Likelihood Inference) web service [70, 71]. Trees were generated using a fast ML stepwise-addition algorithm and 50 attachment braches were evaluated for each taxon. The parameters used were as follows: rate matrix-LG, state frequencies-estimate, proportion of invariable sites-estimate, rate heterogeneity model-gamma distribution, number of rate categories-4. Bootstrap analysis was done at 500 replicates. Phylogenetic trees were then viewed using FigTree (Fig. 1a–b) .
Poly(A) + RNA Isolation
Four milligrams of microspores were cultured in spring water for the desired amount of time. RNA was isolated from the developing gametophytes using the Magnetic mRNA Isolation Kit from NEB (#S1550S, New England Biolabs) according to manufacturer instructions.
Primers for MvKinesin-2, MvKinesin-9A, MvKinesin-9B, and MvCentrin were generated using sequence data from the gametophyte transcriptome (Additional file 11).
Reverse Transcription Polymerase Chain Reaction (RT-PCR)
Reverse transcription reactions were performed using poly(A+)-RNA isolated at 1–2 h, 3–5 h, and 6–8 h of gametophyte development. AMV RT enzyme (New England Biolabs) was used to catalyze the reaction according to manufacture instructions and 5 ng of polyA + RNA was used for each reaction. PCR was carried out using 10 μl of cDNA from each RT reaction and amplified using Taq polymerase (New England Biolabs) for 30 cycles. RT-PCR products were run on a 1 % TAE agarose gel and visualized using ethidium bromide and UV light.
dsRNA was constructed and RNAi was performed as described previously . Spectrophotometry and gel electrophoresis was used to determine the quality and quantity of dsRNA used for RNAi (Fig. 2b). dsRNA construct sequences can be found in Additional file 11. RNAi was performed by adding 200 μg dsRNA to the microspores with 1 ml of spring water in a 2 ml microcentrifuge tube [42, 64]. After 1 h, gametophytes were transferred to 50 ml flasks containing a total of 10 ml spring water and cultured with rotational shaking at 20 °C.
Fixing, embedding, and sectioning microspores
Staining and Immunofluorescence
Toluidine Blue O (TBO) staining was performed on fixed and sectioned microspores and observed with bright field microscopy as previously described . Immunolabeling of sectioned microspores was preformed as previously described [74, 75]. The primary antibody used was mouse anti-centrin clone 20H5 (Milipore 04-1624) diluted to 1:200 in PBS. The secondary antibody was Alexa Fluor 594-conjugated goat anti-mouse (Molecular probes cat# A11005) diluted to 1:1000 in PBST. DAPI was added at 2.5 μg/ml during the last 10 min of staining . Incident light fluorescence imaging was conducted on a Zeiss Axio microscope with standard Fluorescein, TexasRed, and UV filter sets. For all images, thousands of microspores are treated, sectioned, and viewed. About 100 microspores were imaged for each treatment. Percentages were calculated by counting the incidence of each phenocopy.
Analyzing Spermatozoid Swimming Behavior
Gametophytes were allowed to develop for 10.5–16 h and viewed by placing several drops of culture media containing emerging and fully developed spermatozoids in a shallow dish with a coverslip for a bottom. Spermatozoids were analyzed using differential interference contrast (DIC) microscopy with a Zeiss Axiovert 135 TV microscope. Movies were taken with a PIKE F032B monochrome camera (Allied Vision Technology) using StreamPix software at 30 frames per second. Movie files were converted to .mov format. Still frames from movies were captured and converted into tiff format for analysis.
Visualizing Whole Spermatozoids
Spermatozoids were viewed by transferring a drop of the culture media at 11 h to a glass microscope slide. A drop of 4 % PFA was added to fix the spermatozoids rapidly. A coverslip was then placed on top of the drop containing spermatozoids, culture media, and PFA. To prevent crushing fully developed spermatozoids on the slide, the edges of the coverslip were coated in a thin layer of Vaseline. Spermatozoids were then visualized using DIC microscopy.
ATP, adenosine triphosphate; B, blepharoplast; BB, basal bodies; Blastp, basic local alignment sequence tool, protein; CEM, cephalic male cilia; Cr, Chlamydomonas reinhardtii; DAPI, 4’, 6-diamidino-2-phenylindole; DIC, differential interference contrast microscopy; dsRNA, double stranded RNA; e-value, expectancy value; FLA8, flagellar assembly protein 8; FLA10, flagellar assembly protein 10; GARLI, genetic algorithm for rapid likelihood inference; IFT, intraflagellar transport; J, jacket cells; KAP, kinesin associated protein; KCBP, kinesin calmodulin binding protein; KIF, kinesin family member; KLP, kinesin-like protein; mRNA, messenger RNA; ML, Maximum likelihood; Mv, Marsilea vestita; OSM-3, osmotic avoidance abnormal protein 3; PFA, paraformaldehyde; PFR, paraflagellar rod; Poly(A), polyadenylated; RNA, ribonucleic acid; RNAi, RNA interference; RNAseq, RNA sequencing; RT-PCR, reverse transcriptase polymerase chain reaction; SP, spermatogenous cells; Tb, Trypanosoma brucei; TBO, toluidine blue-o
Hepler PK. The blepharoplast of Marsilea: Its De Novo Formation and Spindle Association. J Cell Sci. 1976;21:361–90.
Hart PE, Wolniak SM. Spermiogenesis in Marsilea vestita: A temporal correlation between centrin expression and blepharoplast differentiation. Cell Motil Cytoskeleton. 1998;41:39–48.
Klink VP, Wolniak SM. Changes in the abundance and distribution of conserved centrosomal, cytoskeletal and ciliary proteins during spermiogenesis in Marsilea vestita. Cell Motil Cytoskeleton. 2003;56:57–73.
Wolniak SM, Boothby TC, van der Weele CM. Posttranscriptional control over rapid development and ciliogenesis in Marsilea. Methods Cell Biol. 2015;127:403–44.
Deeb F, Wolniak SM, van der Weele CM. Spermidine is a morphogenetic determinant for cell fate specification in the male gametophyte of the Water Fern Marsilea vestita. Plant Cell. 2010;22:3678–91.
Boothby TC, Wolniak SM. Masked mRNA is stored with aggregated nuclear speckles and its asymmetric redistribution requires a homolog of mago nashi. BMC Cell Biol. 2011;12:45.
Boothby TC, Zipper RS, van der Weele CM, Wolniak SM. Removal of retained introns regulates translation in the rapidly developing gametophyte of Marsilea vestita. Dev Cell. 2013;24:517–29.
Sharp LW. Spermatogenesis in Marsilia. Bot Gaz. 1914;58:419–31.
Mizukami I, Gall J. Centriole replication. II. Sperm formation in the fern, Marsilea and the cycad, Zamia. J Cell Biol. 1966;29:97–111.
Myles DG. An ultrastructural study of the spermatozoid of the fern, Marsilea vestita. J Cell Sci. 1975;17:633–45.
Myles DG, Hepler PK. Spermiogenesis in the fern Marsilea: microtubules, nuclear shaping, and cytomorphogenesis. J Cell Sci. 1977;23:57–83.
Wolniak SM, van der Weele CM, Deeb F, Boothby T, Klink VP. Extremes in rapid cellular morphogenesis: post-transcriptional regulation of spermatogenesis in Marsilea vestita. Protoplasma. 2011;248:457–73.
Myles DG, Southworth D, Hepler PK. Cell surface topography during Marsilea spermiogenesis: flagellar reorientation and membrane particle arrays. Protoplasma. 1987;93:405–17.
Scholey JM. Kinesin-2: A family of heterotrimeric and homodimeric motors with diverse intracellular transport functions. Annu Rev Cell Dev Biol. 2013;29:443–69.
Huang B, Rifkin MR, Luck DJL. Temperature-sensitive mutations affecting flagellar assembly and function in Chlamydomonas reinhardtii. J Cell Biol. 1977;72:67–85.
Walther Z, Vashishtha M, Hall JL. The Chlamydomonas FLA10 gene encodes a novel kinesin-homologous protein. J Cell Biol. 1994;126:175–88.
Kozminski KG, Beech PL, Rosembaum JL. The Chlamydomonas kinesin-like protein FLA10 is involved in motility associated with the flagellar membrane. J Cell Biol. 1995;131:1517–27.
Cole DG, Diener DR, Himelblau AL, Beech PL, Fuster JC, Rosembaum JL. Chlamydomonas kinesin-II-dependent intraflagellar transport (IFT): IFT particles contain proteins required for ciliary assembly in Caenorhabditis elegans sensory neurons. J Cell Biol. 1998;141:993–1008.
Signor D, Wedaman KP, Rose LS, Scholey JM. Two heteromeric kinesin complexes in chemosensory neurons and sensory cilia of Caenorhabditis elegans. Mol Biol Cell. 1999;2:345–60.
Setou M, Nakagawa T, Seog DH, Hirokawa N. Kinesin superfamily motor protein KIF17 and mLin-10 in NMDA receptor-containing vesicle transport. Science. 2000;288:1796–802.
Snow JJ, Ou G, Gunnarson AL, Walker MR, Zhou HM, Brust-Mascher I, Scholey JM. Two anterograde intraflagellar transport motors cooperate to build sensory cilia on C. elegans neurons. Nat Cell Biol. 2004;11:1109–13.
Evans JE, Snow JJ, Gunnarson AL, Ou G, Stahlberg H, McDonald KL, Scholey JM. Functional modulation of IFT kinesins extends the sensory repertoire of ciliated neurons in Caenorhabditis elegans. J Cell Biol. 2006;172:663–9.
Imanishi M, Endres NF, Gennerich A, Vale RD. Autoinhibition regulates the motility of the C. elegans intraflagellar transport motor OSM-3. J Cell Biol. 2006;174:931–7.
Pan X, Ou G, Civelekoglu-Scholey G, Blacque OE, Endres NF, Tao L, Mogilner A, Leroux MR, Vale RD, Scholey JM. Mechanisms of transport of IFT particles in C. elegans cilia by the concerted actions of kinesin-II and OSM-3 motors. J Cell Biol. 2006;174:1035–45.
Han YG, Kwok BH, Kernan MJ. Intraflagellar transport is required for Drosophila to differentiation sensory cilia but not sperm. Curr Biol. 2003;13:1679–86.
Sarpal R, Todi SV, Sivan-Loukianova E, Shirolikar S, Subramanian S, Raff EC, Erickson JW, Ray K, Eberl DF. Drosophila KAP interacts with kinesin II motor subunit KLP64D to assemble chordotonal sensory cilia, but not sperm tails. Curr Biol. 2003;13:1687–96.
Briggs LJ, Davidge JA, Wickstead B, Ginger ML, Gull K. More than one way to build a flagellum: comparative genomics of parasitic protozoa. Curr Biol. 2004;14:611–2.
Wickstead B, Gull K. A “holistic” kinesin phylogeny reveals new kinesin families and predicts protein functions. Mol Biol Cell. 2006;17:1734–43.
Demonchy R, Blisnick T, Deprez C, Toutirais G, Loussert C, Marande W, et al. Kinesin 9 family members perform separate functions in the trypanosome flagellum. J Cell Biol. 2009;187:615–22.
Bernstein M, Beech PL, Katz SG, Rosenbaum JL. A New Kinesin-like Protein (Klp1) Localized to a Single Microtubule of the Chlamydomonas Flagellum. J Cell Biol. 1994;125:1313–26.
Yokoyama R, OToole E, Ghosh S, Mitchell DR. Regulation of flagellar dynein activity by a central pair kinesin. Proc Natl Acad Sci U S A. 2004;101:17398–403.
Lechtreck KF, Witman GB. Chlamydomonas reinhardtii hydin is a central pair protein required for flagellar motility. J Cell Biol. 2007;176:473–82.
Dawe HR, Shaw MK, Farr H, Gull K. The hydrocephalus inducing gene product, Hydin, positions axonemal central pair microtubules. BMC Biol. 2007;5:33.
Lechtreck KF, Delmotte P, Robinson ML, Sanderson MJ, Witman GB. Mutations in Hydin impair ciliary motility in mice. J Cell Biol. 2008;180:633–43.
Marande W, Kohl L. Flagellar kinesins in protists. Future Microbiol. 2011;6:231–46.
Wickstead B, Gull K, Richards TA. Patterns of kinesin evolution reveal a complex ancestral eukaryote with a multifunctional cytoskeleton. BMC Evol Biol. 2010;10:110.
Shen Z, Collatos AR, Bibeau JP, Furt F, Vidali L. Phylogenetic analysis of the Kineisn superfamily from Physcomitrella. Front Plant Sci. 2012;3:230.
Tomei EJ, Wolniak SM. Transcriptome analysis reveals a diverse family of kinesins essential for spermatogenesis in the fern Marsilea. Cytoskeleton. 2016;73:145–59.
Reddy ASN, Day IS. Kinesins in the Arabidopsis genome: A comparative analysis among eukaryotes. BMC Genomics. 2001;2:2.
Manton I, Clarke B. Demonstration of compound cilia in a fern spermatozoid with the electron microscope. J Exp Bot. 1951;2:125–8.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Klink VP, Wolniak SM. Centrin is necessary for the formation of the motile apparatus in spermatids of Marsilea. Mol Biol Cell. 2001;12:761–76.
Myles DG, Hepler PK. Shaping of the sperm nucleus in Marsilea: a distinction between factors responsible for shape generation and shape determination. Dev Biol. 1982;90:238–52.
Awan A, Bernstein M, Hamasaki T, Satir P. Cloning and characterization of Kin5, a novel Tetrahymena ciliary kinesin II. Cell Motil Cytoskeleton. 2004;58:1–9.
Fan J, Beck KA. A role for the spectrin superfamily member Syne-1 and kinesin II in cytokinesis. J Cell Sci. 2004;117:619–29.
Haraguchi K, Hayashi T, Jimbo T, Yamamoto T, Akiyama T. Role of the kinesin-2 family protein, KIF3, during mitosis. J Biol Chem. 2006;281:4094–9.
Keil R, Kiessling C, Hatzfeld M. Targeting of p0071 to the midbody depends on KIF3. J Cell Sci. 2009;122:1174–83.
Morris RL, Scholey JM. Heterotrimeric kinesin-II is required for the assembly of motile 9 + 2 ciliary axonemes on sea urchin embryos. J Cell Biol. 1997;138:1009–22.
Matsuura K, Lefebvre PA, Kamiya R, Hirono M. Kinesin-II is not essential for mitosis and cell growth in Chlamydomonas. Cell Motil Cytoskeleton. 2002;52:195–201.
Miller MS, Esparza JM, Lippa AM, Lux 3rd FG, Cole DB, Dutcher SK. Mutant kinesin-2 motor subunits increase chromosome loss. Mol Biol Cell. 2005;16:3810–20.
Le Bot N, Antony C, Whitle J, Karsenti E, Vemos I. Role of xklp3, a subunit of the Xenopus kinesin II heterotrimeric complex, in membrane transport between endoplasmic reticulum and the Golgi apparatus. J Cell Biol. 1998;143:1559–73.
Stauber T, Simpson JC, Pepperkok R, Vernos I. A role for kinesin-2 in COP1-depenent recycling between the ER and the Golgi complex. Curr Biol. 2006;16:2245–51.
Nekrasova OE, Amargo EV, Smith WO, Chen J, Kreitzer GE, Green KJ. Desmosomal cadherins utilize distinct kinesins for assembly into desmosomes. J Cell Biol. 2011;195:1185–203.
Lipka E, Müller S. Potential roles for kinesins at the cortical division site. Front Plant Sci. 2012;3:158.
Miki T, Naito H, Hishina M, Goshima G. Endogenous localizome identifies 43 mitotic kinesins in a plant cell. Proc Natl Acad Sci U S A. 2014;111:E1053–61.
Mukhopadhyay S, Lu Y, Qin H, Lanjuin A, Shaham S, Sengupta P. Distinct IFT mechanisms contribute to the generation of ciliary structural diversity in C. elegans. EMBO J. 2007;26:2966–80.
Morsci NS, Barr MM. Kinesin-3 KLP-6 regulates intraflagellar transport in male-specific cilia of Caenorhabditis elegans. Curr Biol. 2011;21:1239–44.
Piao T, Luo M, Wang L, Guo Y, Li D, Li P, Snell WJ, Pan J. A microtubule depolymerizing kinesin functions during both flagellar disassembly and flagellar assembly in Chlamydomonas. Proc Natl Acad Sci U S A. 2009;106:4713–18.
Wang L, Piao T, Cao M, Qin T, Huang L, Deng H, Mao T, Pan J. Flagellar regeneration requires cytoplasmic microtubule depolymerization and kinesin-13. J Cell Sci. 2013;126:1531–40.
Dymek EE, Goduti D, Kramer T, Smith EF. A kinesin-like calmodulin-binding protein in Chlamydomonas: evidence for a role in cell division and flagellar functions. J Cell Sci. 2006;119:3107–16.
Wolniak SM, Cande WZ. Physiological requirements for ciliary reactivation of bracken fern spermatozoids. J Cell Sci. 1980;43:195–207.
Hyams JS. Binding of tetrahymena dynein to axonemes of Marsilea vestita lacking the outer dynein arm. J Cell Sci. 1985;73:299–310.
Hyams JS, Campbell CJ. Widespread absence of outer dynein arms in the spermatozoid of lower plants. Cell Biol Int Rep. 1985;9:841–8.
Tsai CW, Wolniak SM. Cell cycle arrest allows centrin translation but not basal body formation during spermiogenesis in Marsilea. J Cell Sci. 2001;114:4265–72.
NCBI Database. http://www.ncbi.nlm.nih.gov. Accessed Feb 2016.
Phytozome Database, http://phytozome.jgi.doe.gov/pz/portal.html. Accessed Feb 2016.
Cosmoss Database. http://www.cosmoss.org. Accessed Feb 2016.
NCBI Conserved Domain Database. http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml. Accessed Feb 2016.
MEGA Software. http://www.megasoftware.net. Accessed Feb 2016.
Genetic Algorithm for Rapid Likelihood Inference (GARLI) Web Service. http://www.molecularevolution.org. Accessed Feb 2016.
Bazinet AL, Zwickl DJ, Cummings MP. A Gateway for phylogenetic analysis powered by grid computing featuring GARLI 2.0. Syst Biol. 2014;63:812–8.
FigTree Software. http://tree.bio.ed.ac.uk/software/figtree/. Accessed Feb 2016.
O’Brien TP, McCully ME. The study of plant structure. Principles and selected methods. Melbourne: Termarcarphi Pty Ltd; 1981.
Baskin TI, Wilson JE. Inhibitors of protein kinase and phosphatases alter root morphology and disorganize cortical microtubules. Plant Physiol. 1997;113:493–502.
Van der Weele CM, Tsai CW, Wolniak SM. Mago Nashi is essential for spermatogenesis in Marsilea. Mol Cell Biol. 2007;18:3711–22.
Notredame C, Higgins DG, Heringa J. T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;8:205–17.
Di Tommasso D, Moretti S, Xenarious I, Orobitg M, Montanyola A, Chang JM, Taly JF, Notredame C. T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension. Nucleic Acids Res. 2011;39:13–7.
We gratefully acknowledge project support from the Maryland Agricultural Experiment Station and comments, suggestions, and early-morning help from our colleagues, Dr. Corine van der Weele, Dr. Stephen M. Mount, Dr. Thomas C. Boothby, and Mr. Richard Zipper.
Grant support for this work from the National Science Foundation (MCB 0842525 to SMW) and project support from the Maryland Agricultural Experiment Station.
Availability of data and material
The datasets supporting the conclusions of Fig. 1a and b are available in the TreeBase repository [http://purl.org/phylo/treebase/phylows/study/TB2:S19463]. The datasets supporting the conclusions of the remainder of this article are included within the article and its additional files.
EJT performed experiments, participated in the design of experiments, and interpretation of the data, and writing of the manuscript. SMW participated in the design of experiments, the interpretation of data, and the writing of the manuscript. EJT and SMW both read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Identifying kinesin-2 and kinesin-9 sequences in Marsilea. (A) Protein sequences from members of the kinesin-2 and kinesin-9 family in Physcomitrella patens and Chlamydomonas reinhardtii were used to search the Marsilea transcriptome. (B) Stand alone blast used to search the transcriptome for kinesin-2 and kinesin-9 sequences. Results show that there is one kinesin-2 mRNA sequence (KT986235) and two kinesin-9 mRNA sequences (KT986258, KT986259) in the male gametophyte of Marsilea. (PDF 3038 kb)
Building a kinesin-2 and a kinesin-9 phylogenetic tree. Source of kinesin-2 and kinesin-9 motor domains used for multiple sequence alignment (MSA) and to build phylogenetic trees. (PDF 1418 kb)
Domain architecture and multiple sequence alignments for (A) kinesin-2, (B) kinesin-9A, and (C) kinesin-9B. Bold regions correspond to the motor domains of each respective kinesin. Regions that are important for motor domain function, including P-loop (yellow), switch I (blue), and switch II (green), are represented in each motor. Regions that are within the motor domain are highlighted and are very highly conserved. Non-motor regions are somewhat conserved, however, large gaps and insertions exist in each alignment. Alignments generated using T-coffee [76, 77]. (PDF 1034 kb)
Untreated, control spermatozoids emerge from the microspore. (MOV 4129 kb)
Untreated, control spermatozoids exhibit normal swimming behavior by spinning in place. (MOV 3274 kb)
MvKinesin-2 knockdowns emerging from the microspore. (MOV 5165 kb)
MvKinesin-2 knockdowns swimming in place exhibiting multiple coils with attached cilia on one cell body (Monster). (MOV 4106 kb)
MvKinesin-2 knockdowns swimming in place exhibiting long cilia (Rapunzel). (MOV 4527 kb)
MvKinesin-9A knockdowns swimming erratically. (MOV 4190 kb)
MvKinesin-9B knockdowns swimming in a similar manner to control spermatozoids (Late Bloomer). (MOV 4043 kb)
Primers and dsRNA constructs. Primers for MvKinesin-2, MvKinesin-9A, MvKinesin-9B, and MvCentrin were generated from the gametophyte transcriptome for (A) RT-PCR and (B) to make dsRNA. (C) Full dsRNA constructs for RNAi. (PDF 564 kb)