extfieldsofstudy: Biology Medicine

text:
Acute Versus Chronic Loss of Mammalian Azi1/Cep131 Results in Distinct Ciliary Phenotypes

Defects in cilium and centrosome function result in a spectrum of clinically-related disorders, known as ciliopathies. However, the complex molecular composition of these structures confounds functional dissection of what any individual gene product is doing under normal and disease conditions. As part of an siRNA screen for genes involved in mammalian ciliogenesis, we and others have identified the conserved centrosomal protein Azi1/Cep131 as required for cilia formation, supporting previous Danio rerio and Drosophila melanogaster mutant studies. Acute loss of Azi1 by knock-down in mouse fibroblasts leads to a robust reduction in ciliogenesis, which we rescue by expressing siRNA-resistant Azi1-GFP. Localisation studies show Azi1 localises to centriolar satellites, and traffics along microtubules becoming enriched around the basal body. Azi1 also localises to the transition zone, a structure important for regulating traffic into the ciliary compartment. To study the requirement of Azi1 during development and tissue homeostasis, Azi1 null mice were generated (Azi1Gt/Gt). Surprisingly, Azi1Gt/Gt MEFs have no discernible ciliary phenotype and moreover are resistant to Azi1 siRNA knock-down, demonstrating that a compensation mechanism exists to allow ciliogenesis to proceed despite the lack of Azi1. Cilia throughout Azi1 null mice are functionally normal, as embryonic patterning and adult homeostasis are grossly unaffected. However, in the highly specialised sperm flagella, the loss of Azi1 is not compensated, leading to striking microtubule-based trafficking defects in both the manchette and the flagella, resulting in male infertility. Our analysis of Azi1 knock-down (acute loss) versus gene deletion (chronic loss) suggests that Azi1 plays a conserved, but non-essential trafficking role in ciliogenesis. Importantly, our in vivo analysis reveals Azi1 mediates novel trafficking functions necessary for flagellogenesis. Our study highlights the importance of both acute removal of a protein, in addition to mouse knock-out studies, when functionally characterising candidates for human disease.

Introduction
Centrosomes are conserved animal organelles which function as the major microtubule organising centre (MTOC), and are required for diverse processes including formation of cilia and flagella, intracellular trafficking events, cell polarity and division. Structurally, the centrosome consists of a pair of cylindrical centrioles surrounded by a proteinaceous matrix of pericentriolar material (PCM) [1]. Importantly, centrioles replicate only once per cell cycle and are essential for the formation of cilia, key signalling organelles during development and homeostasis. In post-mitotic cells, the centrosome moves to the apical surface where the mother centriole docks with the cell membrane to become a basal body and a template for the axonemal microtubules of the primary cilium. Cilia assembly and function requires diverse trafficking events, including intraflagellar transport (IFT), coordinated by the basal body 'hub', which regulates traffic in and out of the ciliary compartment. Enrichment of key signalling receptors and downstream effectors in cilia allows these structures to function as effective signalling organelles, with unique protein and lipid composition [2]. The transition zone, a highly specialised structure just distal to the basal body, is thought to act as an additional ciliary gate, controlling traffic into and out of the cilium [3]. Many of these aspects of ciliogenesis are highly conserved [4].
Despite these common features, cilia are also structurally and functionally diverse. Cilia play important sensory roles, acting as transducers of developmental signalling pathways, detecting fluid flow, as well as highly specialised sensory receptors [5]. Some cilia are motile, involved in generating fluid flow in the embryonic node, airways, oviduct and brain, as well as in the propulsion of sperm. How the core ciliary assembly programme is modified and elaborated on to account for these species-and cell-specific variations is not well understood [4].
Mutations in conserved cilial and centrosomal genes have been identified in a growing spectrum of clinical disorders, termed ciliopathies, with both distinct and overlapping clinical features including polydactyly, skeletal defects, situs inversus, infertility and neuropathology [6,7]. Proteomic and genetic studies in several organisms estimate the molecular composition of cilia/centrosomes to include hundreds to thousands of putative components, many of them unknown [8,9]. Functional dissection of the role and requirement of many of these ciliopathy candidates in cilia formation and function are often performed using cell culture [10,11,12] and zebrafish knock-down models [13]. Mouse mutant models are analysed less often as these are more costly in time and resources to produce. Given the phenotypic complexities of clinical features in ciliopathies [14], what is the best way to understand the underlying molecular mechanisms for candidate genes in relation to human disease?
Recently, centriolar satellites have reported to be the site of localisation of many ciliopathy proteins, and are involved in their ciliary targeting, including OFD1 (oral-facia1-digital syndrome 1), BBS4 (Bardet-Biedl Syndrome 4) and CEP290 [15]. Conserved among vertebrates, but not present in arthropods, centriolar satellites are electron dense, multi-protein complexes enriched in the area surrounding the centrosome/basal body [16,17,18]. These are dynamic structures trafficking along microtubules towards the centrosome utilising dynein motors [16,18]. Centriolar satellites have been shown to regulate ciliogenesis and centriole biogenesis, in part by regulating trafficking of proteins to and/or sequestering of proteins away from the centrosome/basal body [19,20,21,22]. Centriolar satellites are defined by pericentriolar material 1 (PCM1) which is a key scaffolding component of centriolar satellites and to date, all centriolar satellite-localised proteins have been shown to interact with PCM1 [17]. However, what the functional significance of vertebrate-specific centriolar satellites to mammalian development is and how they affect the function of highly conserved components in ciliogenesis and centriole biogenesis is unknown.
Here, we address the role and requirement of 5-azacytidine induced gene 1 (Azi1)/Cep131 (MGI:107440), a highly conserved centrosomal protein, in ciliogenesis. We screened a subset of cilia-enriched orthologous candidates from Drosophila melanogaster studies [23] by RNAi to identify genes involved in mammalian ciliogenesis, and identified Azi1/Cep131. This finding agrees with a previous siRNA screen, which showed a role for the human orthologue, AZI1, in ciliogenesis [10]. Both Drosophila melanogaster mutants and Danio rerio morphants of Azi1 (dila/CG1625 and cep131, respectively) phenocopy mutations of known ciliary genes [24,25], suggesting Azi1 plays a conserved function in ciliogenesis. AZI1 was described as a centrosomal protein (Centrosomal protein 131: Cep131) in a large scale proteomics screen and this localisation was recently refined to the centriolar satellites [26,27,28]. The mouse protein is highly expressed in the testes in germ cells during the period of flagellar formation [29]. More recently, additional roles for AZI1 include involvement in genome stability and centriole duplication. Knock-down of AZI1 leads to an increase in double-stranded DNA breaks, indicated by cH2AX staining, as well as a slight increase in cells with extra centrioles [28,30]. However, little is known of the in vivo role of mouse Azi1 and its requirement for development.
Here we utilise knock-down, localisation and live-imaging techniques, to further investigate the role of Azi1 in mammalian ciliogenesis at the cellular level. To determine the requirement for Azi1 in mouse development, we generated Azi1 null mutant mice and focused on the in vivo role of Azi1 in ciliogenesis and genome stability. Our analysis of Azi1 knock-down (acute loss) versus gene deletion (chronic loss) suggests that Azi1 plays a conserved, but non-critical trafficking role in ciliogenesis. Importantly, our in vivo analysis reveals Azi1 mediates novel trafficking events necessary for spermiogenesis and male fertility.

Results
Functional cell-based screening of putative ciliary candidates from Drosophila identifies Dila orthologue Azi1 as required for mammalian ciliogenesis Using a set of forty orthologous putative ciliary genes identified as highly expressed in ciliated cell types in Drosophila melanogaster [23], we carried out an siRNA screen in a mouse fibroblast cell line to identify genes involved in mammalian ciliogenesis. Cilia formation was assayed as the percentage of cells with a cilium, marked by anti-Arl13b, a ciliary membrane marker [31], and antiacetylated a-Tubulin, a ciliary axoneme marker, using an automated imaging and image analysis system.
We identified Azi1/Cep131 as the top hit for genes involved in cilia formation, with at least two of four siRNAs giving a significant reduction in ciliogenesis across three independent assays (data not shown). This observation is supported by the study of Graser et al. (2007), who found a reduction in ciliogenesis on AZI1 knock-down in human hTERT-RPE1 cells [10]. To exclude off-target effects of the siRNAs, we co-transfected a different pool of four siRNAs, specifically targeting the 39 untranslated region (UTR) of only Azi1, along with a plasmid encoding either GFP or Azi1-GFP. The Azi1-GFP plasmid lacks the 39 UTR of Azi1 and so is resistant to these siRNA. Transfection of Azi1 39 UTR siRNA leads to a reduction in Azi1 protein to 10% of wild type levels, which can be partially rescued by co-transfection with Azi1-GFP ( Figure 1A). Azi1 knock-down leads to a 50% reduction of transfected cells with cilia ( Figure 1B, D and F), similar to that seen with a positive control siRNA targeting Ift88, a gene essential for ciliogenesis [32,33]. Importantly, co-expression of Azi1-GFP rescues the phenotype back to control levels demonstrating that the ciliary phenotype observed upon addition of Azi1 siRNA is not due to off-target effects of the siRNA ( Figure 1C, E and F). We conclude that Azi1 is involved in mammalian cilia formation.

Author Summary
Cilia are slender projections from the surface of most mammalian cells and have sensory and sometimes motile functions. They are essential for mammalian development and defects in cilia lead to a group of human diseases, termed ciliopathies, with variable symptoms including embryonic lethality, lung and kidney defects, blindness and infertility. Cilia are complex structures composed of hundreds of components, whose individual functions are poorly understood. We screened for mammalian genes important in building cilia, and identified Azi1/Cep131, a gene previously shown to be required for cilia formation and function in fish and flies. We show that if we acutely reduce levels of Azi1 in mouse cells, fewer cells form cilia, but if we generate cells chronically lacking all Azi1, cilia form normally. In addition, mice without any Azi1 are healthy and viable, confirming normal cilia function. However, in these mice, the highly specialised ciliary structure of the sperm tail does not form, resulting in male infertility. We suggest Azi1 has conserved trafficking roles in both primary cilia and the specialised sperm flagella. Abruptly removing Azi1 results in instability causing the existing cilia network to collapse, whereas chronic deletion of Azi1 allows the system to re-equilibrate, and cilia to form normally. Azi1 knock-down leads to reduced ciliogenesis. Mouse ShhLIGHT-II fibroblast cells were transfected with siRNA (a non-targeting control siRNA (Ctrl), two positive control siRNAs targeting Ift88 (Ift88 #3 and #4) or a pool of four siRNAs targeting the 39 UTR of Azi1 (Azi1 39UTR)), along with plasmids encoding either GFP or Azi1-GFP (which lacks the 39UTR of Azi1). (A) Western blot probed with anti-Azi1 antibody (SF91) shows endogenous Azi1 levels are reduced upon Azi1 39UTR siRNA addition (lower arrow). Tagged Azi1-GFP protein is also detected (upper arrow). The blot was reprobed with anti-a Tubulin as a loading control. Below is quantification of total Azi1 levels (endogenous plus overexpressed), relative to a-Tubulin. (B-E) Cells were stained with anti-Arl13b to mark the cilia and GFP Booster (Chromotek) to enhance the GFP signal. Transfected cells are highlighted with arrowheads; closed arrowheads highlight cells with cilia, open arrowheads highlight cells without cilia. (F) Addition of Azi1 39UTR siRNA significantly reduced the percentage of cells with cilia, and this reduction is rescued to wild type levels by co-transfection with Azi1-GFP. Shown is the mean +/2 SEM of two technical and two biological duplicates (*** P,0.001, chi-squared test). Scale bars represent 10 mm. doi:10.1371/journal.pgen.1003928.g001 Azi1 traffics along microtubules towards the centrosome/ciliary base, where it localises to the transition zone AZI1 was originally identified as a centrosomal protein (Cep131 [26]). We have spatially refined the localisation of mouse endogenous and GFP-tagged Azi1 to centriolar satellites, marked by anti-PCM1, which confirms recent human AZI1 immunofluorescence reports ( Figure S1A and G) [28,34]. We identified a further pool of human and mouse Azi1 at the transition zone (Figures 2A, B and S1B-C), indicated by co-staining with antipolyglutamylated tubulin, which stains the ciliary axoneme and basal body but, importantly, is absent from the transition zone [35]. The transition zone is an area at the base of the cilia involved in regulating traffic into the cilium [3]. Co-staining with anti-NPHP1, a marker of the transition zone [36,37], confirms this localisation ( Figures 2C, D and S1D). This is consistent with the observation that in D. melanogaster ciliated sensory neurons Azi1 homologue dila localises distal to the basal body at the putative transition zone [24].
Recently, CEP290 has also been reported to localise to both centriolar satellites and the transition zone [38,39,40], raising the possibility this could be a general trend for centriolar satellite proteins. To test this, we investigated the localisation of PCM1, the core component of centriolar satellites at the transition zone. Indeed, PCM1 localises to the transition zone of most, but not all cilia, as indicated by the gap between the basal body and axoneme on anti-polyglutamylated tubulin staining ( Figure 2E and F). Previous reports had shown OFD1 and PCM1 to similarly localise to the distal portion of basal bodies [15]. Interestingly, the putative functional orthologue of OFD1, UNC, is also found at the putative transition zone of Drosophila mechanosensory neurons [24,41]. Together, this suggests that docking at the transition zone may be a conserved feature of components of mammalian centriolar satellites.
We used live imaging of Azi1-GFP to address the dynamics of Azi1 localisation. In interphase cells, centriolar satellites have been proposed to function in dynein motor-dependent, microtubulebased trafficking of proteins to the centrosome [19,40], and it has been shown recently that the pericentriolar localisation of AZI1 is microtubule dependent [28]. Similarly, trafficking of cargo associated with IFT motors along microtubules into the ciliary compartment is selectively regulated in part by the transition zone [42]. To examine Azi1 trafficking more directly we imaged Azi1-GFP movement in live mouse NIH-3T3 cells. Azi1-GFP was observed to traffic along microtubules, co-labelled with Map4-RFP ( Figure 2G). Azi1-GFP traffics with a dynamic saltatory motion, with periods of fast movement, for an average distance of 3.4 mm (range 1.3-10.2 mm) interspersed with sometimes long stationary periods. Azi1-GFP was observed to move both towards and away from the centrosome, similar to observations of PCM1-GFP [18] (Movie S1). Although average speeds varied according to direction, 1.860.2 mm/s (mean 6 SEM) towards the minus end of microtubules at the centrosome, and 1.060.3 mm/s away from the centrosome, they were consistent with speeds observed previously for microtubule-based motors in vivo [43,44]. This suggests Azi1 is involved in microtubule-based trafficking to and from the centrosome/basal body.
Higher levels of endogenous AZI1 staining are observed around basal bodies and surrounding centriolar satellites of ciliated cells compared to non-ciliated cells ( Figure 2H and I). However, total levels of AZI1 do not change upon serum starvation ( Figure 2J) indicating that under these conditions to induce ciliogenesis, there is a redistribution of AZI1 towards the basal body area.
It has been proposed that centriolar satellites act as proteinaceous scaffolds to physically restrict access of proteins to the centrosome/cilium complex [22]. Disruption of core components results in dissolution or dispersal of centriolar satellites, and relocalisation of associated centriolar satellite proteins to the centrosome/basal body [15,20,22,40]. Despite its localisation, Azi1 siRNA knock-down in mouse cells does not alter centriolar satellite integrity as shown by Pcm1 localisation, consistent with recent reports for human AZI1 [28]. This suggests Azi1 is not required for mammalian centriolar satellite integrity nor retention of Pcm1 to these structures ( Figure S1E-H).

Azi1 is dispensable for mouse embryonic development
Given the high conservation of Azi1 among ciliates (Table S1), including arthropods which lack centriolar satellites, together with phenotypic mutant data from diverse organisms [24,25,45,46], we predicted Azi1 would have a central role in mammalian cilia biology in vivo and thus generated mouse mutants null for Azi1. Azi1 Gt(CCOG35)Wtsi embryonic stem (ES) cells, which have a gene trap inserted into intron 2 of Azi1 ( Figure 3A), were used to generate Azi1 +/Gt(CCOG35)Wtsi mice (referred to as Azi1 Gt/+ throughout this manuscript). Azi1 Gt/Gt mice are born at sub-Mendelian ratios, with approximately two thirds of the expected numbers of Azi1 Gt/Gt mice remaining after weaning (see Table 1, P = 0.0025). A significant reduction in Azi1 Gt/Gt numbers was observed at embryonic day 11.5-13.5 (E11.5-13.5), suggesting roughly a third of mutants are lost before mid-gestation, although further work is needed to determine exactly when this loss occurs (Table 1). Azi1 Gt/Gt mice that are born appear morphologically normal, and are the same weight as wild type littermates ( Figure S5G and H). Viable Azi1 mutant mice showed none of the gross abnormalities associated with cilia dysfunction in mice, including failure to thrive, hydrocephalus, situs inversus, and chronic airway infections.
Azi1 has several coiled-coil domains as well as a predicted t-SNARE domain (IPR010989) implicated in membrane fusion events during vesicular transport ( Figure S2B). The gene trap is predicted to truncate Azi1 after the initial 69 amino acids such that any remaining Azi1 trapped protein in Azi1 Gt/Gt mice will lack all the predicted domains in the more highly conserved C terminus ( Figure S2, Table S1), and hence unlikely to be functional. To confirm that the gene trap eliminated gene expression, we examined Azi1 mRNA expression levels across the gene trap insertion site ( Figure 3B and C). No expression across the insertion site was detected in kidneys, ovaries or testes from Azi1 Gt/Gt mice by RT-PCR or qRT-PCR ( Figure 3B and C), whereas robust expression of the gene trap was detected ( Figure 3B). X-Gal staining of E11.5 Azi1 Gt/+ embryos further confirmed expression of the gene trap b-geo gene, and showed Azi1 expression is ubiquitous during development, with higher expression in tissues with high levels of cilia-dependent developmental signalling such as the limbs, eyes, somite derivatives and brain ( Figure 3D).
We detected some low level expression of the 39 end of Azi1 in Azi1 Gt/Gt mice; when quantified by qRT-PCR this was less than 2% of wild type levels ( Figure S2C-E). Importantly, no Azi1 protein was detected in Azi1 Gt/Gt mice when probing with an antibody raised against the C terminal of Azi1 ( Figure 3E and S2B), despite a single strong band of the expected size in Azi1 +/+ and a band of roughly 50% intensity in Azi1 Gt/+ samples. Furthermore, anti-Azi1 immunofluorescence analysis of Azi1 Gt/Gt mutant mouse embryonic fibroblasts (MEFs) or multiciliated airway epithelial cells detected no signal, despite clear centrosomal/basal body localisation of Azi1 in littermate controls ( Figure 3F-H). This suggests that any low level transcription  The lack of wild type transcription across the gene trap insertion site was confirmed by qPCR on testes cDNA. (D) LacZ staining of E11.5 Azi1 Gt/+ embryos demonstrates the gene trap-bGeo cassette is expressed. Azi1 expression is ubiquitous during development with higher expression in the limbs, somite derivatives, eyes and brain. (E) No full length or truncated Azi1 protein expression was detected in Azi1 Gt/Gt mutant testes. A strong band was detected in wild type and at lower levels in Azi1 +/Gt testes extract using an antibody raised to the C-terminal of Azi1 (SF91) (See Figure S3A). Densitometry was used to give Azi1 levels relative to a-Tubulin, normalised to 1 in wild type (lower of panel D). Azi1 Gt/Gt mice are null for Azi1. Cilia, centrioles and centriolar satellites of Azi1 null MEFs are grossly normal Primary cilia are required for mammalian development, so the fact most Azi1 null mice survive without any patterning defects suggests Azi1 is not required for mammalian ciliogenesis in vivo. In contrast, transient siRNA knock-down of Azi1 leads to a two-fold reduction in ciliogenesis ( Figure 1). On careful examination of cilia formation in Azi1 null primary MEFs we found normal numbers of cilia, marked by anti-Arl13b ( Figure 4A, F and K). Moreover, cilia compartmentalisation appeared normal, with correct distributions of the ciliary membrane protein Arl13b, the IFT-B protein Ift88, and components of the transition zone Nphp1 and Mks1 ( Figure 4A-D, F-I). Post-translational modifications of tubulin also appeared normal: acetylated a-Tubulin is present at the basal body and axoneme, polyglutamylated tubulin is present at the axoneme and basal body but absent from the transition zone, and c-Tubulin is present at the basal body ( Figure 4A-J). Localisation of Rab8, a small GTPase involved in trafficking to the cilium [21], whose localisation has been shown to be centriolar satellitedependent [40], was also localised normally in Azi1 Gt/Gt MEFs ( Figure 4E and J). Together, this analysis suggests that Azi1 is dispensable for primary cilia formation and compartmentalisation both in vivo and in primary cells.

Genotype
Azi1  we examined formation of centriolar satellites in the absence of Azi1. Similar to Azi1 knock-down, Pcm1 is correctly localised in Azi1 Gt/Gt MEFs ( Figure S3A-C), as are centriolar satellite components Cep72 and Cep290, involved in modulating localisation and composition of centriolar satellites [22,40] (data not shown). This confirms Azi1 is not required for centriolar satellite formation, or localisation of key structural components like Pcm1, nor regulatory components like Cep290 or Cep72. As centriolar satellites are proposed to have a role in regulating centrosome and centriole biogenesis, and AZI1 knock-down in human cells leads to increased Centrin2-positive foci [17,19,28,47], we examined the numbers of centrioles (marked by Centrin2-GFP and anti-Centrin3) and centrosomes (marked by anti-c Tubulin) in cells lacking Azi1 ( Figure 4A, F, L-R). We found both centrosome and centriole numbers were normal in Azi1 Gt/Gt MEFs ( Figure 4A, F and N). Multiciliated cells have the ability to assemble hundreds of centrioles through two parallel pathways [48,49], in both of which, protein-rich fibrous granules, akin to centriolar satellites, marked by PCM1, are found surrounding the elongating centrioles [18]. After assembly in the cytoplasm, these centrioles move apically to dock at the plasma membrane and each extend a ciliary axoneme. Ultrastructural analysis of motile multiciliated epithelia of adult trachea revealed normal numbers and docking of basal bodies and appendage formation in Azi1 null mice ( Figure S3D and G). Confirming the immunofluorescent analysis in Azi1 mutant cells, transition zone ultrastructure appeared morphologically normal in Azi1 null mice ( Figure S3E and H).
In summary, and in contrast with the acute Azi1 knock-down in both human and mouse cells ( Figure 1) [10,28], cilia and centriole structure appears grossly normal in Azi1 null cells.

Azi1 null MEFs have compensated for the loss of Azi1
The difference in ciliary phenotypes observed with acute Azi1 knock-down versus its chronic absence in Azi1 Gt/Gt MEFs is intriguing. To rule out differences between the cell line used for the screen and primary cells, we co-transfected wild-type MEFs with siRNA against Azi1 along with plasmids encoding either GFP or Azi1-GFP and examined cilia formation. Similar to the results obtained in the embryonic fibroblast cell line ( Figure 1), Azi1 knock-down in primary MEFs led to a robust reduction in ciliogenesis which was rescued by the siRNA-resistant Azi1-GFP ( Figure 5A-D and I). To rule out any residual Azi1 function in our mutant Azi1 Gt/Gt cells, and to further eliminate off-target effects, we transfected Azi1 Gt/Gt MEFs with Azi1 siRNA. While our positive control Ift88 siRNA gave a robust reduction in ciliogenesis in Azi1 Gt/Gt MEFs, ciliogenesis was unaffected upon Azi1 knockdown ( Figure 5E-I), demonstrating that Azi1 null MEFs have compensated for the loss of Azi1. Discrepancies between the phenotypic severity observed with siRNA knock-down versus genetic deletion has previously been attributed to the acute nature of knock-down, allowing less time for compensation to occur [50,51]. We conclude that although Azi1 is involved in cilia formation in mouse, compensation during embryogenesis in the absence of Azi1 allows ciliogenesis to proceed normally in most tissues of Azi1 null mice.
A global DNA damage response is not observed in Azi1 null mice A role for cilia/centrosomal proteins in genome stability and DNA damage response pathways has been proposed [52], and was recently reported for AZI1 [28,30]. Unlike the reported AZI1 knock-down [28,30], we saw no increase in the number or intensity of cH2AX foci in Azi1 null MEFs ( Figure S4A, F and K and data not shown). To test whether Azi1 mutants were more susceptible to DNA damaging agents, we challenged MEFs with either hydroxyurea (HU) or ionising radiation and examined cH2AX foci. No significant difference was observed between genotypes at lower concentrations of HU or upon challenge with ionising radiation, although Azi1 null MEFs were more susceptible to high doses of HU (5 mM) ( Figure S4B-E, and G-K). As Staples et al. (2012) reported an increase in micronuclei in AZI1 depleted U2OS cells [28], we examined Azi1 Gt/Gt mice using a peripheral blood micronuclei assay -a highly sensitive method for detecting in vivo DNA damage [53]. Whilst positive control Mcph1 Gt/Gt mice show a marked increase in micronucleated erythrocytes as previously documented (http://www.sanger.ac.uk/mouseportal/ phenotyping/MBGX/micronuclei/) [54], Azi1 null mice show no such increase, suggesting there is no elevation in DNA damage in the peripheral blood of Azi1 null mice ( Figure S4L and M). Taken together, there is no gross evidence for chromosomal instability in Azi1 null animals, neither in vivo nor in primary cell culture. Once again it is possible that the differences observed between siRNA knock-down and genetic null of Azi1 are due to compensation for Azi1 loss in the genetic null, as seen for the ciliogenic phenotypes observed.

Azi1 function is not compensated in the sperm flagella resulting in male infertility
We examined Azi1 null mice carefully for adult-onset ciliopathic phenotypes to determine whether cilia formation and function was normal in all tissues. Azi1 mutant males display complete infertility with no evidence of pregnancy or pups born from more than 25 plugged dams. In contrast, Azi1 Gt/Gt female mice have normal litter numbers and sizes ( Figure 6A). Azi1 Gt/Gt males have reduced testes weight (Azi1 Gt/Gt : 185 mg+/223.9 mg vs. Azi1 Gt/+ : 127 mg+/29.7 mg, P,0.05, Student's t-test) corresponding to a drastic reduction in sperm density (less than 2% of wild type) ( Figure 6B). Male infertility is a common symptom of ciliopathies, often in conjunction with airway dysfunction and late-onset phenotypes including retinal degeneration, kidney cysts and obesity in mouse mutants of cilia genes [55,56,57]. Immunofluorescent and ultrastructural analysis of postnatal multiciliated airway epithelium revealed mutant cilia to be morphologically normal, consistent with the lack of chronic airway infections in these mice ( Figure 3G-H and S3D-I). In aged cohorts of littermates, no signs of retinal degeneration were observed in Azi1 null eyes either by ophthalmoscopic examination (data not shown) or histologically at six months (n = 7, Figure S5A-C). No cysts were observed in Azi1 null kidneys aged 6 months or older (n = 6, Figure S5D-F), nor was any obesity observed in aged Azi1 Gt/Gt mice (n = 8 at 3 months, n = 5 at 6 months, Figure S5G and H). In support of our cellular analysis, these in vivo studies of Azi1 null mice demonstrate that Azi1 is not required for cilia formation or function in general but is required for formation and function of the specialised cilia derivative, the sperm flagella.

Spermiogenesis arrests at
Step 9 in Azi1 null tubules, with severe flagellar defects and teratozoospermia To determine when spermatogenesis is disrupted in Azi1 mutant mice, we examined periodic acid-Schiff (PAS) stained sections of testes. Azi1 Gt/Gt mutant testes show a significant reduction in tubule lumen size (Lumen diameter: 15.1+/21.2 mm in Azi1 Gt/Gt , compared to 36.7+/21.4 mm in Azi1 +/+ , mean +/2 SEM, P,0.001, Student's t-test, n = 3), with a drastic reduction in the number of sperm flagella visible in the lumens of Azi1 Gt/Gt testes ( Figure 6C-L). Anti-acetylated a-Tubulin, which marks the  axonemes, reveals that whilst control tubule lumens are filled with sperm flagella, almost none were detected in mutant tubules, and any flagella seen appeared shorter and morphologically abnormal, suggesting Azi1 is necessary for flagella formation ( Figure 6M and N). Light microscopy revealed the pre-spermiogenic stages of spermatogenesis, up to Step 7-8 spermatids, appeared to be normal ( Figure 6E-L, O-P). However, from Step 9, spermatid morphology is highly abnormal, with mutant elongating spermatids mislocalised and misorientated within the tubule. In addition to axonemal defects, mutant elongating spermatids also exhibit teratozoospermia, where sperm heads show an abnormal club-shaped nuclear morphology, as opposed to the normal hook-shaped head in wild type elongate spermatids ( Figure 6E-L, P). Very few spermatids reach maturity and are successfully passed into the epididymis ( Figure 6B, S6A and B). Importantly these phenotypes were observable from the first wave of spermatogenesis ( Figure S6C-J). Whereas Azi1 Gt/Gt tubules resemble Azi1 +/+ tubules at postnatal day 20 (P20) and P25 ( Figure  S6C-F), by P30, defects such as lack of flagella and misorientated spermatids displaying teratozoospermia become apparent. These results confirm observations in the adult mutant tubules that defects in spermiogenesis arise as the spermatids begin to elongate (P25-30) ( Figure S6G-J).
To investigate whether the reduction in mature sperm was due to an increase in spermatid death in Azi1 Gt/Gt tubules, we analysed levels of activated Caspase 3a and TUNEL staining, both indicators of apoptosis, in adult mutant tubules. Azi1 mutant tubules showed increases in both activated Caspase 3a staining and TUNEL positive cells ( Figure S6K-P). The restricted spatial distribution of activated Caspase 3a positive cells, combined with the cell counts showing a reduction in elongating spermatids ( Figure 6O) indicate Azi1 Gt/Gt spermatids undergo apoptosis from Step 9 onwards.
The germline is particularly sensitive to DNA damage and given the suggested role for AZI1 in genome stability we considered whether Azi1 null male infertility is due to an increase in DNA damage. Infertility due to defects in DNA damage response pathways generally presents as an early arrest in spermiogenesis, with spermatocytes not progressing through meiosis [58,59,60]. This is in contrast to the relatively late arrest in post-meiotic spermatogenesis seen in Azi1 null mice, reminiscent of other ciliopathic mouse models [55,56,57,61,62]. To affirm that the arrest in spermatogenesis in Azi1 mutant mice is not due to genome instability, we stained Azi1 null testes with anti-cH2AX. In Azi1 mutant tubules, we observed anti-cH2AX staining comparable to controls, emphasising that there is no increase in DNA double-stranded breaks in Azi1 null testes ( Figure S4N and O). Together with the previous in vitro and in vivo data ( Figure S4), this shows there is no increase in DNA damage in the absence of Azi1 under physiological conditions.

Loss of Azi1 results in multiple microtubule-dependent trafficking defects in mutant spermatids
Cauda or caput epididymides from Azi1 mutant mice only contained debris and degenerating sperm, none of which were motile ( Figure 6B, 7A and B, Movie S2 and S3, and data not shown). Even in the testes, flagella were rarely observed by ultrastructure analysis of mutant tubules, which exhibited drastically reduced lumen diameters filled with vacuolar cells and proteinaceous cell debris, in contrast to the open, flagella filled wild type tubule lumens ( Figure 7C and D). The rare flagella that remained showed evidence of abnormal trafficking, with swollen flagellar lumens and ectopic mistrafficked outer dense fibres, occasionally associated with microtubules ( Figure 7E and F, S7). Remaining sperm were morphologically abnormal with flagella of severely reduced length ( Figure 7G-I). The processes involved in extension of the mammalian sperm flagella are not well understood, but are thought to involve IFT-mediated trafficking [63,64,65]. Truncated Azi1 mutant axonemes exhibit abnormal post-translational modifications of microtubules and irregular distributions of IFT trains, as marked by anti-Ift88 ( Figure 7G and H). Extended axonemal structures are rarely detected by longitudinal TEM sections of mutant spermatids, instead occasional swollen shortened flagellar remnants with build-up of IFT cargo can be seen ( Figure S7). Together these data suggest disruptions in regulated IFT account for the failure of flagellar axoneme extension.
However, some of the abnormalities in Azi1 mutant spermatid development, such as teratozoospermia cannot easily be explained by defects in IFT. These defects are first observed at step 9, when spermatids undergo a series of complex morphological changes, including nuclear remodelling and formation of the transient microtubule structure of the manchette [66,67]. This microtubular ''sleeve'' structure surrounds the head and is assembled concurrently with the elongation and condensation of the spermatid nucleus, as well as growth of the centrosome-derived axoneme [68]. As in IFT, motor-driven trafficking along this track of microtubules delivers cargo from the Golgi-derived acrosome toward the centrosome and nascent sperm tail, in a process of intramanchette transport (IMT) [67]. Abnormal club shaped nuclei were previously observed in mutants with defects in manchette formation and function [69,70]. Ultrastructural analysis revealed the manchette to be present in Azi1 mutant spermatids but it often appears kinked, and is occasionally misnucleated away from the spermatid head ( Figure 8A and B, S8A and B). Late stage spermatids exhibited abnormal nuclear morphologies, consistent with the histological analyses ( Figure 6 and S6), often with detachment of the acrosome from the nucleus ( Figure S8C-E). Formation of the sperm tail involves the migration and modification of a peripheral pair of centrioles to the caudal pole of the nucleus, opposite the acrosome, where they become lodged to form the neck or head-tail coupling apparatus (HTCA). In early spermatid differentiation, modifications to the proximal centriole, lodging of the centrioles into the implantation fossa of the nuclear membrane and formation of the centriolar adjunct appear grossly normal in Azi1 mutant spermatids [71] (Figure 8C and D). A range of HTCA phenotypes are observed in later stage mutant spermatids, including misalignment of HTCA with the nucleus and/or displaced implantation fossa [72] (Figure 8E and F, S8F). Together, these results suggest Azi1 may be required for maturation and functional integrity of the HTCA.
To understand these spermiogenic phenotypes observed in Azi1 mutants, we examined Azi1 localisation during sperm development. Azi1 is found at the Golgi-derived acrosome (Step 10-12 spermatids: Figure 8G), then in the centrosome-containing HTCA at the flagellar base in later stage spermatids ( Figure 8I). Importantly, no Azi1 staining is detected in any stage of mutant spermatids, confirming specificity of this Azi1 localisation type MEFs show a reduction in the percentage of cells with cilia upon treatment with Azi1 siRNA, which is rescued by overexpressing Azi1-GFP. Azi1 Gt/Gt MEFs show normal cilia numbers when transfected with non-targeting siRNA, and cilia numbers remain unchanged upon treatment with Azi1 siRNA in Azi1 Gt/Gt MEFs. Shown is mean +/2 standard deviation (n = 3). doi:10.1371/journal.pgen.1003928.g005 Azi1/Cep131 in Mammalian Cilia and Flagella  Figure 8H and J). This dynamic stage-specific redistribution of Azi1 is consistent with Azi1 undergoing IMT, although we failed to detect Azi1 specifically in the manchette, possibly due to dispersal of Azi1 below the detection threshold during IMT. This redistribution of Azi1 to the HTCA is consistent with a role in the maturation and functional integrity of this structure, and is reminiscent of enrichment of Azi1 around basal bodies upon ciliation in somatic cells.
We next examined localisation of an IMT cargo, Hook1, a coiled-coil protein implicated in vesicular transport which is mobilised progressively from the acrosome to HTCA during spermiogenesis [73,74]. Post-translational modification of microtubules in the manchette are proposed to determine trafficking events by motors [74]. Manchette microtubules in wild type spermatids, although stabilised, are not labelled by the usual microtubule stabilising modifications, such as acetylation [75] ( Figure 8K-M). However, strong anti-acetylated a-Tubulin staining is observed in Hook1-positive Azi1 mutant manchettes at step 9-12 ( Figure 8N and O). Subsequently, Hook1 is prematurely lost from Step 14-15 mutant manchettes ( Figure 8M and P). These results suggest Azi1 mutant spermatids exhibit altered microtubule dynamics and IMT cargo localisation. While intramanchette transport (IMT) is essential for both normal sperm head morphology and flagella formation [66,73,76], mouse mutations of components trafficked by IMT, like Hook1 [68,69,73] and RIM-BP3 [70], do not block extension of the flagellar axoneme completely.
These results suggest that loss of Azi1 disrupts the microtubulebased trafficking of both flagellar-directed IMT, and intraflagellar transport, resulting in both abnormal sperm head morphology together with the lack of flagella in Azi1 mutant spermatids.

Discussion
Genomic and proteomic studies on the centrosome/cilium complex (the ''ciliome'') have estimated that thousands of gene products are involved in the formation and function of this important structure [8,9]. Although many of the components are conserved across ciliates, it has generally been assumed that extrapolation from diverse genetic model organisms such as photosensitive, motile flagella of Chlamydomonas and ciliated mechanosensory chordotonal organs of Drosophila can be confidently extended to vertebrate cilia. In particular, there has been a growing trend of morpholino-based knock-down in the D. rerio model to functionally analyse putative ciliopathy candidate genes [13,25]. Several recent mammalian cell culture-based knock-down screens have identified novel components involved in regulating cilia biology [10,11,12]. Facing a flood of next-generation sequencing detecting human disease variants, the daunting challenge is how best to investigate the requirement and function of novel and poorly characterised genes during development and disease.
Although Azi1/Cep131 was predicted to be essential for cilia from studies in fly and zebrafish [24,25], we present here detailed analysis of both acute loss (by siRNA) and chronic absence (by genetic mutation), and show that Azi1 has a conserved but non-essential role in mammalian ciliogenesis, and is essential for the formation of sperm flagella (Figure 9).
Acute versus chronic lack of Azi1: Lessons for functional characterisation of putative human disease genes We and others have shown transient knock-down of mammalian Azi1 leads to a reduction in ciliogenesis [10]. Importantly, we show this Azi1 knock-down ciliogenesis defect is rescued by overexpressing Azi1-GFP, emphasising the phenotype is ''ontarget'' (Figure 1). Surprisingly, following genetic deletion of all gene function, cilia develop and function normally in vivo and in primary cells (Figure 3, 4, S3 and S5). Ruling out any subdetectable ''leakiness'' of our Azi1 gene-trap, transfection of Azi1 siRNA into Azi1 Gt/Gt MEFs does not affect ciliogenesis ( Figure 5), suggesting that compensation for the lack of Azi1 has occurred. It has been previously suggested that acute knock-down of proteins can give a more severe phenotype than long-term deletions due to compensation in vivo [50,51]. We present the first demonstration of a lack of phenotype in null cells treated with siRNA against the gene of interest, proving this functional compensation exists.
In the absence of Azi1, mammalian cilia develop and function properly in all mouse tissues except for the developing sperm ( Figure 3, 4 S3 and S5). Azi1 is not essential for ciliogenesis and any involvement it may have in cilia biology can be compensated for in most tissues, aside from the modified cilia of the sperm flagella. Azi1 null mice are born at sub-Mendelian ratios and a third of mutants appear to be lost before mid-gestation ( Table 1), suggesting that a proportion of mutant embryos may fail to compensate for the loss of Azi1 and die earlier in development. Several ciliopathy mouse models are also born at sub-Mendelian ratios, suggesting stochastic events may affect the requirement for ciliary proteins during embryogenesis [55,56,57,77].
Our study highlights the importance of functional follow-up studies of siRNA data, and cautions against direct extrapolation of siRNA phenotypes to the genetic in vivo phenotype. On the other hand, acute knock-down by siRNA can expose roles for genes that would otherwise be overlooked due to compensation or redundancy in long-term deletion studies in vivo.
Azi1 null male mice are infertile, suggesting the loss of Azi1 cannot be compensated for during spermiogenesis. Azi1 null mice exhibit post-meiotic defects in spermatogenesis with misorientation and abnormal morphology of elongating and elongated spermatids, including teratozoospermia, from Step 9 of spermatid development onwards (Figure 6 and S6). Sperm flagella are mostly absent and any remaining axonemes are truncated and immotile, with swollen flagella lumens and evidence of mistrafficking of IFT components and cargo (Figure 7 and S7), suggesting Azi1 is required for IFT during the formation of the sperm flagella.
Additional roles for ciliary proteins in spermatid development outside the axoneme have been suggested, with parallels drawn between intraflagellar transport (IFT) and intramanchette transport (IMT) [67]. The abnormal club-shaped sperm head morphology is similar to mutants with defective manchettes [68,69,70]. Azi1 mutant manchettes are structurally abnormal, and exhibit altered microtubule post-translational modifications, step 10 onwards. (M and N) Azi1 +/+ and Azi1 Gt/Gt testes sections stained with anti-acetylated a-Tubulin, marking the flagella. Almost no flagella are seen in Azi1 Gt/Gt tubules. Those remaining appear short and morphologically abnormal (Inset, N). (O) Spermatid numbers at each step of development as a percentage of total cells in the tubule. (P) Representative images of spermatids at each step of development taken from (C-K), with a schematic of spermatid morphology, highlighting the acrosome in pink. Spermatid numbers are reduced in Azi1 Gt/Gt tubules from step 9, and morphological abnormalities in the spermatid head are seen from step 12, with a club-shaped morphology. Subtle morphological abnormalities are visible at step 10 with many Azi1 Gt/Gt spermatids lacking the tapered end of the sperm head in. Scale bars represent 50 mm (C-L) or 100 mm (M and N). doi:10.1371/journal.pgen.1003928.g006 Azi1/Cep131 in Mammalian Cilia and Flagella which are implicated in motor-selection for cargo delivery [74] ( Figure 8 and S8). Altered expression of Hook1, an effector of IMT implicated in vesicular transport [74], is also observed in these mutant structures, suggesting the process of IMT is disrupted in Azi1 mutants. Many motors and associated proteins, some linked to IFT, have been localised to the manchette and the testicular phenotype of Ift88 orpk/orpk hypomorphic mice phenocopies the Azi1-dependent trafficking defects [67,74]. Although the role of centriolar satellite or transition zone proteins in spermiogenesis is not well understood, it is tempting to draw parallels between the microtubule-based trafficking of centriolar satellites towards the centrosome/basal body of primary cilia, and the trafficking of proteins along the transient microtubular structure of the manchette to the highly specialised motile sperm flagellum. The dynamic relocalisation of Azi1 from the acrosome to the HTCA during spermatid development suggests Azi1 undergoes IMT, as this is the main transport route between these structures (Figure 8). Interestingly, components of the dynactin complex, which localise to centriolar satellites and is required for their peri-centrosomal localisation, have also been localised to the manchette [19,47,74,78].
It is a conundrum for the cilia field as to why broadly expressed ''core'' cilia genes have clinically restricted pheno-types when mutated. This ''tissue-specific'' requirement is demonstrated by the limited phenotypes of late-onset ciliopathies, and many centriolar satellite and transition zone-specific mouse mutants display tissue specific ciliary phenotypes [56,61,62,77,79,80]. In the case of the transition zone complex components, Nphp4 and Nphp1 null mutant mice display male infertility without kidney phenotypes, in contrast to the human disease in which patients have severe nephronophthisis [61,62,81]. This discrepancy could reflect differences in species and/or functional nature of the mutation, but the severity in humans could also reflect the mutational load in other components of the complex [82,83]. Importantly primary cilia formation and function appears grossly normal in these mutants, like in our Azi1 null mice, emphasising that compensation for function of key cilia genes is likely to be a recurrent theme given  Figure S7A), or is misnucleated away from the spermatid head (B, black arrows, also see Figure S7B). (C-F) Centrioles and the head tail coupling apparatus (HTCA) in early (C and D) or later stage (E and F) elongating spermatids from control (Azi1 +/+ or Azi1 Gt/+ ) (C and E) or Azi1 Gt/Gt testes (D and F). The centrioles implant normally at the early stages of Azi1 Gt/Gt spermatid elongation, and some normal looking later stage HTCA can be found (Fig. S8D), although in many cases late stage HTCA are misaligned, with the implantation fossa or centrioles off-centre (F, brackets and arrow). (G-J) Elongating spermatids prepared by testes squashes from Azi1 +/+ or Azi1 Gt/Gt animals were stained with anti-Azi1. Azi1 localises to the acrosome in Step 10-12 spermatids (G), then to punctae concentrated around centrosomes within the HTCA in control Step 13-14 spermatids (I). Importantly no Azi1 was detected in Azi1 Gt/Gt spermatids, demonstrating the staining is specific.

Evolutionary function and divergence of the requirement for Azi1
Azi1/Cep131 is a conserved protein found in all ciliates except for nematodes, leading to the suggestion that Azi1 may be involved in cilia motility (Table S1) [46]. Alternatively, Azi1/Cep131 could be required to build a canonical nine-triplet centriole, as nematodes build specialized, non-canonical centrioles. However, an Azi1/Cep131 orthologue is present in the Toxoplasma gondii genome which also build non-canonical centrioles [84]. Furthermore, knock-down of Azi1/Cep131 in planaria does not affect centriole formation, but affects cilia motility [45]. While deletion mutations of dila in D. melanogaster [24] lead to defects in specialised ''motile'' mechanosenory Type I neurons and sperm, cilia do form with defects in mechanosensory neuronal cilia morphology characteristic of IFT mutants, suggesting defective trafficking. Similarly, Azi1 D. rerio morphants phenocopy IFT morphants with cilia forming but displaying tissue-specific reductions in length [25]. Interestingly, in trypanosomes TzAZI1 localises to the flagellar pocket, a dynamic endo-exocytic organelle implicated in membrane trafficking surrounding the flagellar base, and TzAZI1 RNAi affects flagellar function, as opposed to its formation or maintenance [46]. These studies suggest a conserved trafficking function for Azi1 in regulating ciliary-bound cargo.
We confirmed mammalian Azi1 localises to centriolar satellites and provide the first direct observation of Azi1 trafficking along microtubules, both towards and away from the centrosome ( Figure 2G, Supplementary Movie 1), similar to the movement observed for PCM1-GFP [18]. The transport away from the centrosome may involve kinesin motors, and it has been shown that CEP290, which binds AZI1, interacts with both Dynactin components and the kinesin motor KIF3a [28,79], supporting the theory that these proteins could undergo bidirectional trafficking along microtubules.
Centriolar satellites are thought to spatially restrict centrosomal access of proteins involved in basal body maturation and ciliogenesis [22]. While centriolar satellites and the key centriolar satellite protein, Pcm1, are not found in Drosophila, tight transcriptional control of dila mRNA limits expression to just before the onset of ciliogenesis where the protein can localise to the PCM [24]. Although the regulatory mechanisms are different, in both mammals and flies, it appears Azi1/dila recruitment is involved in centrosome to basal body maturation. In D. melanogaster, mutations in two other coiled-coil proteins YURI, conserved only among the Drosophila genus [85,86], and UNC, also insect-specific but for which centriolar satellite protein OFD1 is proposed to be a functional orthologue [87] partially phenocopy dila mutants and genetically interact with dila [24]. These proteins are involved in the proper maturation and anchoring of the sperm basal bodies to the nuclear membrane [24,86,87]. Formation of Drosophila sperm flagella axoneme is unusual in that it is IFTindependent, forming instead in cytoplasmic cysts. Sperm axonemes do form in dila mutant flies, although the sperm display defective HTCA formation. This is reminiscent of the defects we observe in compromised integrity of HTCA in Azi1 null sperm (Figure 8 and S8), although the requirement for IFT in mammalian sperm axonemal formation may explain the more severe IFT-based flagellar phenotypes observed in the mouse mutant spermatids.
Interestingly, similar to Ma et al. (2011) who showed transition zone localisation for DILA and UNC in D. melanogaster [24], we show Azi1 and PCM1 also localise to transition zone of primary In mammalian somatic cells, Azi1/Cep131 is a highly conserved and ubiquitously expressed coiled-coil protein found associated with centriolar satellites. During ciliogenesis it becomes enriched during basal body maturation and is found at the transition zone, or ciliary gate, suggesting it plays a role in membrane docking, axoneme extension and/or cargo delivery to nascent cilia. It is one component of a complex and dynamic network required to build and maintain cilia. (B) Transient Azi1/Cep131 knock-down (acute depletion) abruptly reduces levels of this protein and destabilises the entire ciliogenic network, resulting in lack of cilia in somatic cells. (C) Knock-out of Azi1/Cep131 (chronic depletion) in somatic cells allows the system to re-equilibrate the ciliogenic network through compensation and continue to build functional cilia. (D) Knock-out of Azi1/Cep131 (chronic depletion) in developing sperm cannot be compensated for and in its absence, the network never assembles properly. This results in defects in basal body maturation affecting axoneme extension and IFT-dependent cargo delivery necessary to build a functional flagellum. This Azi1/Cep131 null sperm phenotype is similar to Azi1 mutations in other model organisms. doi:10.1371/journal.pgen.1003928.g009 cilium. While this manuscript was in preparation, OFD1 was also shown to localise to the transition zone [41]. Transition zone and centriolar satellite localisation has previously been described for CEP290, which was recently shown to interact with AZI1 [28,38,39]. Our data supports that redistribution of centriolar satellite proteins to the transition zone during ciliogenesis is of functional significance. It is possible that centriolar satelliteassociated cargo is transported to the centrosome as it matures into a basal body with a functional transition zone ready to be trafficked into the cilium. Alternatively, Azi1/CEP290 could also have a role in the gating function at the transition zone, regulating protein content in the cilium, accounting for the IFT-like phenotypes reported in dila Drosophila mutants [24] and Cep290 Chlamydomonas mutants [42]. Given that centriolar satellites are scaffolds for controlling activity of ciliopathy-associated proteins [15,22], it will be important to define the composition of these specific sub-complexes that move to the transition zone, and ask whether they are part of the compensation mechanism observed in Azi1 null animals. If true, we propose that these interactions among the centriolar satellite proteins could extend to multi-allelic mutational load, including AZI1, in a subset of human ciliopathies with diverse clinical presentations beyond male infertility.

Cell culture and transfection
ShhLIGHT-II (ATCC, genetically modified NIH-3T3) and NIH-3T3 cells were maintained in DMEM (Life Technologies), hTERT-RPE cells were maintained in DMEM-F12 (Life Technologies), all supplemented with 1.5 g/L sodium bicarbonate (Sigma), 10% foetal calf serum, 5610 8 U/L penicillin and 11 mM streptomycin at 37uC, 5% CO 2 . Early passage MEFs were maintained in Optimem (Life Technologies) plus 0.5 mM betamercaptoethanol (Sigma), 10% foetal calf serum, 5610 8 U/L penicillin and 11 mM streptomycin at 37uC, 5% CO 2 , 3% O 2 . To induce ciliogenesis, serum was removed for 48 hours. ShhLIGHT-II cells were co-transfected with 25 nM Dharmacon OnTarget Plus siRNA and 1 mg/mL plasmid DNA with Dharmafect Duo (Dharmacon), serum was removed after 24 h and samples were analysed 72 hours after transfection. MEFs were co-transfected with 50 nM Dharmacon OnTarget Plus siRNA and 1.6 mg/mL plasmid DNA using the Invitrogen Neon, according to the manufacturer's protocol. For DNA damage assays, MEFs were challenged with hydroxyurea (Sigma) at given concentrations for 18 hours. Alternatively, MEFs were irradiated in culture medium at 2Gy/ minute using a Faxitron CellRad cabinet X-ray system (Faxitron Bioptics), cultured for 3 hours and then fixed and analysed.

Generation of plasmids
Azi1 was PCR amplified from mouse cDNA and TA cloned into pcDNA6.2-C-Em-GFP-GW-TOPO cDNA plasmid (Life Technologies). Centrin2 was PCR amplified from mouse cDNA, adding EcoR1 and Sal1 sites. This was restriction enzyme cloned into pEGFP-N1 (Clontech).

Western blot and densitometry
Cells or testes were homogenised in cell lysis buffer (Cell Signaling Technology) plus 1 mM phenylmethylsulfonyl fluoride (PMSF) (Thermo Scientific) and Complete Protease Inhibitor Cocktail (Roche), sonicated for 2630 seconds. Testes extracts were concentrated with Amicon Ultra-0.5 mL 30 kDa centrifugal filters (Millipore). Samples were separated on Novex 3-8% Tris Acetate gels (Life Technologies) then transferred to Hybond nitrocellulose membrane (GE Healthcare), which were blocked in 5% milk. Membranes were incubated in primary antibodies (Table  S4) washed, incubated with horse radish peroxidase (HRP)conjugated secondaries (Table S5) and developed with Amersham ECL-plus western blotting detection system. Densitometry was performed using ImageJ.

Immunocytochemistry
Cells were fixed in 4% paraformaldehyde/phosphate buffered solution (PFA/PBS) for 10 minutes at room temperature, or 4% PFA/PHEM (120 mM PIPES, 140 mM HEPES, 20 mM EGTA, 16 mM MgSO 4 , pH7) for 10 minutes at 37uC. Alternatively, preextraction was performed for 30 seconds on ice in 0.1 M PIPES pH 6.8, 2 mM EGTA and 1 mM MgSO 4 , then cells were fixed in ice cold methanol on ice for 10 minutes. Cells were blocked in 10% donkey serum/0.1% Triton-X in Tris buffered solution (TBS). Primary antibodies were added (Table S4), cells were washed then incubated in Alexa Fluor-conjugated secondaries (Table S5) and slides were mounted using Prolong Gold (Life Technologies).

Histology and immunohistochemistry
Tissue samples were collected, kidneys were fixed in 4% PFA/ PBS, testes were fixed in Bouin's fixative, and eyes were fixed in Davidson's fixative according to standard protocols. Tissues were serially dehydrated and embedded in paraffin. Microtome sections of 8 mm thickness were examined histologically via haematoxylin and eosin (H&E) or periodic acid-Schiff (PAS) staining.
For immunofluorescent analysis, paraffin sections were dewaxed and re-hydrated via ethanol series. Antigen retrieval was performed by boiling the sections for 15 minutes in the microwave in citrate buffer. Sections were blocked in 10% donkey serum/ 0.1% Triton-X in PBS and primary antibodies were diluted in 1% donkey serum/PBS (Table S4). Slides were washed and incubated in Alexafluor conjugated secondary antibodies (Table S5), washed and mounted in ProLong Gold (Life technologies).
For immunohistochemistry, the same procedure was used, with the addition of one step after the re-hydration. Slides were immersed in 3% H 2 O 2 in PBS for 20 minutes to block endogenous peroxidases. Slides were incubated in primary antibody, washed, then incubated in biotin-conjugated secondary antibody (Vector laboratories). This was detected using the Vector ABC kit and Vector NovaRed peroxidase substrate kit.
For TUNEL, after dewaxing, sections were incubated in 0.25% Triton-X and labelling was performed using Click-IT TUNEL staining kit following manufacturer's instructions (Life Technologies).

Sperm preparation
Testes, cauda and caput epididymides were dissected into M2 media (Invitrogen). For live imaging, sperm were imaged in M2 media or 1% methyl cellulose (Sigma), in capillary tubes (Vitrotubes Mountain Leaks) sealed with Cristaseal (Hawskley). Sperm counts were performed on sperm from the cauda epididymides, diluted in H 2 O using a haemocytometer, only counting intact sperm (with both head and tail). For fixed samples, either sperm spreads or testes squashes were prepared. For sperm spreads, testes were placed through a 100 mm nylon mesh (BD Biosciences). Sperm from the caput epididymus and testes were then placed on a 20-40% Percoll gradient (GE Healthcare) and spun at 3,000 g for 10 minutes. The sperm were spread on Poly-D-lysine slides (BD Biosciences) and fixed 4% PFA/PBS. Cauda sperm was spread directly on Poly-D-Lysine slides and fixed with 4% PFA/PBS. In all cases, sperm were permeabilised with 0.4% Triton-X in PBS and immunofluorescence was performed as described. For testes squashes, tubules were dissected as described [88], placed on a slide and squashed with a coverslip. Slides were snap frozen in liquid nitrogen, coverslips removed and samples fixed and permeabilised by 10 minutes in 220uC methanol, 30 seconds in acetone and then 15 minutes in 4% PFA/PBS. Immunofluorescence was then performed as described.

Nasal brushing
Nasal brushing was performed as described [89]. Cells were fixed for 30 minutes on ice in 4% PFA, cytospun onto slides, fixed with 220uC methanol for 10 minutes then immunofluorescence was performed as described.

Transmission electron microscopy
Samples were dissected into PBS. Samples were fixed in 2% PFA/2.5% gluteraldehyde/0.1 M Sodium Cacodylate Buffer pH7.4 + 0.04% CaCl 2 . Testes capsules were removed prior to immersion in fix. After 30 minutes at room temperatures, samples were cut into 1 mm cubes and fixed overnight or longer at 4uC. Tissue was rinsed in 0.1 M sodium cacodylate buffer, post-fixed in 1% OsO 4 (Agar Scientific) for one hour and dehydrated in sequential steps of acetone prior to impregnation in increasing concentrations of resin (TAAB Lab Equipment) in acetone followed by 100%, placed in moulds and polymerised at 60uC for 24 hours.
Ultrathin sections of 70 nm were subsequently cut using a diamond knife on a Leica EM UC7 ultramicrotome. Sections were stretched with chloroform to eliminate compression and mounted on Pioloform filmed copper grids prior to staining with 1% aqueous uranyl acetate and lead citrate (Leica). They were viewed on a Philips CM100 Compustage Transmission Electron Microscope with images collected using an AMT CCD camera (Deben).

Micronuclei assay
Mice were tail tipped and blood was collected using a microhematocrit capillary tube with heparin coating (Globe Scientific) into heparin (Sigma). This was fixed in super-chilled methanol. Saline was added and then blood was pelleted at 600 g for 5 minutes. Pellets were treated with RNAseA and anti-CD71primary antibody (Lifespan Biosciences). Cells were incubated on ice for 30 minutes then at room temperature for 30 minutes. Propidium iodide was added and flow cytometry was performed on a FACScalibur (BD Biosciences). Data was analysed using FlowJo software (v7.6.1, Tree Star) as described by [53].

Imaging and image analysis
The initial siRNA screen imaging was carried out on the Olympus Scan R microscope, imaging 16 frames per well. Image analysis, including identification and counting of cells and cilia was performed using the Olympus Scan R Analysis Software.
Confocal images were captured with a Nikon A1R confocal microscope, comprising a Nikon Eclipse TiE inverted microscope and four laser modules: 405 (laser diode), 457, 488, 514 (multiline Argon) 561 (diode-pumped solid-state) and 638 nm (laser diode). For live imaging of Azi1-GFP and MAP4-RFP, a Zeiss Axiovert 200 fluorescence microscope was used equipped with 1006/1.4 plan apochromat objective (Carl Zeiss, Welwyn, UK), Lambda LS 300W Xenon source with liquid light guide and 10-position excitation, neutral density and emission filterwheels (Sutter Instrument, Novato, CA), ASI PZ2000 3-axis XYZ stage with integrated piezo Z-drive (Applied Scientific Instrumentation, Eugene, OR) and a Photometrics Coolsnap HQ 2 CCD camera (Roper Scientific, Tucson, AZ). Brightfield images were captured with a Coolsnap HQ CCD camera (Photometrics Ltd, Tucson, AZ) on a Zeiss Axioplan II fluorescence microscope with Planneofluar objectives (Carl Zeiss, Welwyn Garden City, UK). Colour additive filters (Andover Corporation, Salem, NH) installed in a motorised emission filter wheel (Prior Scientific Instruments, Cambridge, UK) were used sequentially to collect red, green and blue images that were then superimposed to form a colour image.
For live imaging of sperm, a Qimaging Retiga camera running at 30 frames per second (bin 262 half frame) captured image sequences with a 56 objective at zoom 5 on a Nikon AZ100 macroscope. Still figures show all time points superimposed into one image to depict the movement, or lack thereof, during the movie.
Apart from the initial screen, image analysis including intensity profiling was performed in ImageJ.

Statistics
Throughout P,0.05 is considered significant. Statistics were carried out in Microsoft Excel or GraphPad Prism (La Jolla, CA) and the test used is specified in the text/figure legend. Figure S1 Mouse and human AZI1 localise to centriolar satellites and the transition zone, but are not required for centriolar satellite formation. Mouse and human AZI1 localises to centriolar satellites and the transition zone in NIH-3T3 cells. (A) Co-staining with anti-Azi1 SF91 and anti-Pcm1 shows mouse Azi1 localises to centriolar satellites in NIH-3T3 cells. Enlargements to the right highlight the co-localisation of Pcm1 (middle) and Azi1 assessed in the testes by anti-cH2AX staining. As DNA doublestranded breaks occur during meiosis, cH2AX staining occurs normally in certain cell types such as in foci in some spermatogonia and preleptotene to zygotene spermatocytes, in the XY body of pachytene spermatocytes and in the nucleus of round spermatids [90]. Comparing stage matched tubules (assessed by serial PAS sections), the Azi1 Gt/Gt testes showed no increase in cH2AX staining. In both Azi1 +/+ and Azi1 Gt/Gt tubules, anti-cH2AX stains the XY body (white arrows) of pachytene spermatocytes (PS) and, less intensely, round spermatids (RS). It is absent from elongating spermatids (ES). Scale bars represent 50 mm (A-J) or 25 mm (N and O). (TIF) Figure S5 No late-onset ciliopathy phenotypes are observed in Azi1 Gt/Gt mice. (A-C) Eyes of seven mutant mice at 6 months of age were directly examined by opthalmoscope and no retinal degeneration was observed (data not shown). This was confirmed histologically on H&E stained wax sections in four mutants. Shown are representative sections of the retina from wild type (A) or two Azi1 null mice (B and C). (D-F) Kidneys were taken from seven mutants aged 6 months or more, and no kidney cysts were seen. This was confirmed by H&E stained sections for four mutants. (E and F) Representative kidney sections from 6 months old Azi1 Gt/Gt mice, and Azi1 Gt/+ littermate, showing no cysts in the mutant kidneys. (G and H) The weight of male (G) and female (H) Azi1 Gt/Gt mice is no different to Azi1 +/+ and Azi1 Gt/+ mice at 1-9 months old. Scale bars represent 50 mm.


extfieldsofstudy: Medicine Biology

text:
Cleavage Events and Sperm Dynamics in Chick Intrauterine Embryos

This study was undertaken to elucidate detailed event of early embryogenesis in chicken embryos using a noninvasive egg retrieval technique before oviposition. White Leghorn intrauterine eggs were retrieved from 95 cyclic hens aged up to 54-56 weeks and morphogenetic observation was made under both bright field and fluorescent image in a time course manner. Differing from mammals, asymmetric cleavage to yield preblastodermal cells was observed throughout early embryogenesis. The first two divisions occurred synchronously and four polarized preblastodermal cells resulted after cruciform cleavage. Then, asynchronous cleavage continued in a radial manner and overall cell size in the initial cleavage region was smaller than that in the distal area. Numerous sperms were visible, regardless of zygotic nuclei formation. Condensed sperm heads were present mainly in the perivitelline space and cytoplasm, and rarely in the yolk region, while decondensed sperm heads were only visible in the yolk. In conclusion, apparent differences in sperm dynamics and early cleavage events compared with mammalian embryos were detected in chick embryo development, which demonstrated polarized cleavage with penetrating supernumerary sperm into multiple regions.

Introduction
Avian models have tremendous value as ex vivo-model systems for both basic and clinical purposes, enabling monitoring of cell differentiation, transformation, and organogenesis under specific conditions. Nevertheless, limited work has been conducted due to technical difficulties in egg retrieval before oviposition. Furthermore, avian embryos demonstrate discoidal meroblastic cleavage with a large amount of yolk and a small germinal disc [1,2], which hinders monitoring early embryo development. Indeed, very little information on early development before oviposition has been reported [3][4][5][6] in comparison with that available after laying of stage X [7]. In this study, we employed a non-surgical intrauterine egg collection by abdominal rubbing [7], which contributes to overcoming current technical limitation.
Lots of information on cell-fate determination occurring in early embryogenesis was given in a variety of invertebrate and vertebrate species [8][9][10]. Differing from mammals, polyspermic penetration was physiologically occurred in avian eggs, but detailed observation has not been reported to date. In this study, we employed a non-invasive egg retrieval technique with comparative classifying of egg shall formation and embryogenesis for monitoring details of sperm penetration and early cleavage events.

Experimental animals
White Leghorn (WL) hens (54-56 weeks old) were used for the collection of intrauterine eggs. We managed chickens according to our standard operation protocol. Relevant experimental procedures for the study were approved by the Institutional Animal Care and Use Committee, Seoul National University before undertaking experiments (SNU-070823-5).

Collection of intrauterine eggs from hens
Intrauterine eggs retrieved from WL hens were harvested by an abdominal massage technique slightly modified from Eyal-Giladi and Kochav [7]. Briefly, the abdomen of hens was pushed gently until exposure of the shell gland, and the surface of the shell gland expanded when an egg was located there for eggshell formation. After expansion of the surface of the shell gland, massage was used to move the egg gently toward the cloaca until the intrauterine egg was released ( Figure 1A).

Analysis of cleavage stages in the intrauterine embryos
Intrauterine embryos were separated from the egg using sterilized paper [11] and the shell membrane and albumen were detached from the yolk. A piece of square-type filter paper (Whatman, Maidstone, Kent, UK) with the hole at the center was placed over the germinal disc. After cutting around the paper containing the intrauterine embryo, it was gently turned over and transferred to saline buffer to further remove the yolk and the vitelline membrane for embryo collection [12]. Collected embryos were fixed with 4% (w/v) paraformaldehyde in 1× phosphate-buffered saline (PBS) and the fixed embryos were classified according to the cleavage stages proposed by Eyal-Giladi and Kochav [7]. Unfertilized and abnormal embryos were identified by the morphological criteria of cleavage furrows.
Photographs of the dorsal part of intrauterine embryos were taken with a stereoscopic zoom microscope (SMZ1000; Nikon Corporation, Tokyo, Japan) and EG&K stage I-II embryo was cultured in Chamlide incubator system (Live Cell Instrument, Seoul, Korea) at 41.5 °C with 5% of CO 2 gas for live cell imaging. Shell membrane-and albumen-detached eggs were put into the 25 ml plastic cup (40025; SPL Life Sciences, Pocheon, Korea) with 10 ml of albumen on the bottom and the surface area on the top was covered with 3 ml of albumen. For retaining embryo viability, all procedures were undertaken less than five minutes in the heated room (more than 30 °C). Timelapse images were taken by multi-purpose zoom confocal microscope (AZ100; Nikon Corporation, Tokyo, Japan).

Phalloidin and DAPI staining of intrauterine embryos
After fixation with 4% paraformaldehyde, the intrauterine embryos were washed in PBS three times and incubated in 0.1% (v/v) Triton X-100 in PBS (PBST). The fixed embryos , and III stages were designated, which were equal to embryonic EG&K stages I-IV, V-VI, and VII-X, respectively. The phase I stage represented as an egg with a yellowish and soft eggshell, phase II was an egg with a light yellow-colored, flexible eggshell, and phase III was an egg with stiffened and calciumdeposited eggshell with a milky-white color. doi: 10.1371/journal.pone.0080631.g001 were incubated with Alexa Fluor 488 phalloidin (A12379; Invitrogen, Carlsbad, CA, USA) diluted 1:40 in PBST overnight at room temperature. After overnight incubation, the embryos were washed three times in PBS and mounted with Prolong Gold antifade reagent with 4',6-diamidino-2-phenylindole (DAPI) (P36931; Invitrogen). The stained embryos were observed under a fluorescence microscope (Ti-U; Nikon Corporation). In addition, the intrauterine embryos were embedded with paraffin and sectioned (12 µm) using a microtome and after being mounted with Prolong Gold antifade reagent with DAPI, the embryonic nuclei were evaluated under a fluorescence microscope.

Statistical analysis
Statistical analyses were performed using the Student t test in SAS version 9.3 software (SAS Institute, Cary, NC). The significance levels between control and treatment groups were analyzed using the general linear model (PROC-GLM) in SAS software. Differences between treatments were deemed to be significant when P was less than 0.05.

Retrieval of intrauterine eggs
The general procedure for the noninvasive collection of intrauterine eggs by abdominal massage is shown in Figure 1. This procedure resulted in minimal stress to the hens, which continued to lay eggs from the second day after harvest. Ninety-five WL hens at 54-56-week-old were provided for egg retrieval, and intrauterine eggs were retrieved from all hens. Among the 95 collected embryos, 38 were of EG&K stage I, 26 of stage II, 11 of stage III, 13 of stage IV, and 7 of stage V. In total, 67.4% of the harvested intrauterine embryos were classified as early EG&K stages I-II. Intrauterine eggs can be divided into three categories based on morphological characteristics ( Figure 1B): eggs with a yellowish soft eggshell membrane of EG&K stages I-V, eggs with a light yellowish flexible eggshell of EG&K stages V-VII, and eggs with a milkywhite stiffened eggshell of EG&K stages VII-X. Eggshell formation advanced gradually in the shell gland. The calciumdeposited eggshell was well formed during EG&K stages V-VI (8 h in the shell gland), hardening of the eggshell was observed at EG&K stage VII, and eggshell formation was complete by EG&K stages IX-X. Overall times to retrieve each stage were expected to be 0-8, 8-12, and 12-20 h after entering into shell gland for phases I, II and III, respectively.

Morphogenesis of cleavage furrows in intrauterine embryos
Of the 38 EG&K stage I embryos collected from the shell glands, five were undergoing the first cleavage ( Figure 2A). The first cleavage furrow was observed in the central region, while a few showed the initiation of cleavage in the peripheral area. Six of the 38 underwent synchronous cleavage up to the third cleavage, perpendicular to the previous cleavage furrow. The fourth cleavage separates central and peripheral cells (schematic diagram; Figure 2B). Distinguishable from the main cleavage furrows formed in a cruciform manner, peripheral cleavage furrows were formed at the embryo boundary until EG&K stage V (Figure 2A, C). The peripheral furrows disappeared gradually after EG&K stage V and became invisible. During cell divisions between EG&K stages I and V, cell size decreased gradually and was approximately tenfold smaller (from 250-300 to 15-40 μm) at EG&K stage V than that of the first cleavage stage ( Table 1). As shown in Table 1, preblastodermal cells, indicating completely closed cells detached from the yolk, were detected from EG&K stage III, but the size varied due to rapid cleavage after EG&K stage II. The subgerminal cavity was initially formed with completely closed cells beginning at EG&K stage III. At EG&K stage IV, the central cells began to form cell layers, and three to six cell layers were detected at EG&K stage V; at this stage, preblastodermal cells were observed in both the central and the peripheral regions ( Figures 2C, 3A).
To further examine cell division, time-lapse live-imaging of the cleaving embryo (EG&K stage I-II) was taken ( Figure 4). Cleavage of two laterally closed cells at the central region, which were indicated as '1' and '2' in the first panel of Figure  4A, was monitored during 4 hours of culture. Asymmetric division with asynchronous cleavage was notable in the observation of two cells. The cell surface area of the cell number '1' was 11258.92 μm 2 at onset of culture and those of its daughter cells were 6855.68 and 3711.55 μm 2 at 58 minutes after culture, that indicated asymmetric division in each of the two cells ( Figure 4B left). In terms of cleavage duration, the second division in one of daughter cells of the cell number '1' completed at 144 minutes after the onset of culture, while that in the other daughter cell completed at 204 minutes after the onset of culture, that indicated asynchronous division ( Figure  4B left). The cell number '2' also showed asymmetric division during culture ( Figure 4B right).
To trace the division direction of open cells, time-lapse liveimaging of the total three cleaving embryos (EG&K stage I) was taken and the one representative embryo is shown in

Localization of F-actin to the cleavage furrows and division patterns in intrauterine embryos
Nuclear and F actin staining respectively with DAPI and phalloidin was conducted to examine the cleavage pattern of intrauterine embryos. Strong F actin staining was detected in the main cleavage furrow and in the peripheral area of EG&K stage I embryos ( Figure 6). Subsequently, F-actin was detected strongly in the second and third cleavage furrows. The newly developed cleavage furrows appeared not to be initiated from the dorsal surface, but rather from deeper  cytoplasmic regions underneath the surface ( Figure 6A, B). During this early stage, F-actin-stained cleavage furrows from the center did not reach the peripheral area of the embryos ( Figure 6). F-actin-stained peripheral cleavage furrows were formed in an irregular (linear, dot-shaped, circular) manner ( Figure 6C). From EG&K stage I, the dividing cells in the center became closed first (Figure 6), whereas the peripheral cells were still open before stage IV. Closed cells in the peripheral area were detected primarily in stage IV, and the majority of cells were completely closed in stage V ( Figure 6H, I). Doublestaining with phalloidin and DAPI clearly showed cell division patterns in the intrauterine embryos in EG&K stages II-V ( Figure 6).

Embryonic and supernumerary sperm nuclei in the intrauterine embryos
Three types of nuclei were observed in the intrauterine embryos according to their morphology, size and position: embryonic (zygotic) nuclei, condensed supernumerary sperm nuclei, and decondensed supernumerary sperm nuclei. Condensed sperm nuclei were mainly present in the dorsal surface and cytoplasm, and rarely in the yolk region underneath the cytoplasm (Figures 3, 7) with a linear shape ( Figure 8B), whereas the decondensed sperm nuclei were spread in the peripheral yolk region and yolk region underneath the cytoplasm (Figures 3, 7) with an irregular shape and smaller size compared to embryonic nuclei ( Figure 8B). Also, the three-dimensional depth coding image showed that the decondensed sperm nucleus was located under the cytoplasm, while embryonic nuclei were in the cytoplasm ( Figure 8A). Less than ten to several thousand condensed and decondensed supernumerary sperm nuclei were detected in the cleavage stages of intrauterine embryos. In particular, the numbers of condensed supernumerary sperm on the dorsal side of EG&K stages I-III embryos ranged from 1 to 10 to more than 1000 per embryo (Table 2). However, late EG&K stage embryos contained very low numbers of supernumerary sperm nuclei. It was obvious that observed sperms were penetrated because the vitelline membrane of all embryos was removed before staining. In the yolk on the ventral side, only decondensed sperm heads were observed in the majority of embryos. In a few embryos, a few condensed sperm heads were also observed on the ventral side as well as decondensed sperm heads.
To examine the spatial distribution of supernumerary sperm nuclei, condensed and decondensed sperm nuclei on the dorsal and ventral side of EG&K stage I-II embryos were counted respectively (Figure 9). On the dorsal side, condensed sperm nuclei and embryonic nuclei were detectable while decondensed sperm nuclei were present on the ventral side ( Figure 9A2, 9B2). Also, the mean number of condensed sperm nuclei was significantly higher on the periphery region than center region ( Figure 9A3). The number of condensed and decondensed sperm nuclei per 1 mm 2 of cell surface area was shown in Figure 9A4 and 9B4.

Discussion
The finding of this study clearly demonstrated different aspects of sperm penetration and embryo cleavage between birds (chicken) and mammals. There was a unique, radiating progress of preblastoderm furrowing which showed different furrowing status between the dorsal and the ventral surfaces. Interestingly, different status of spermatozoa penetrated into egg preblastoderm was detected and uneven distribution of condensed and decondensed sperm heads were detected in central (furrowing-completed, cleavage-initiated region) and peripheral (furrowing-incomplete, cleavage-progressing region) parts of the preblastoderm. Although it was not certain whether supernumerary sperm move from the center toward the periphery, it was obvious that they were abundant in the periphery than the center. To clarify the exact function of supernumerary sperm on cleavages, what components of sperm contribute to embryos should be identified in further studies. In the yolk on the ventral side, decondensed sperm nuclei were mainly detected, which might imply either the presence of decondensation factor in the yolk or the entry of sperm into the yolk area only through the preblastodermal region. In any case, this is the unique phenomenon in chick embryos, which is not seen in the mammals.
In this study, we used a modified noninvasive collection method (abdominal massage) for retrieving intrauterine EG&K stage embryos, which was originally reported by Eyal-Giladi and Kochav [7]. Based on this original technique, we provided the detailed information for the classification, which reflected egg shell formation and a compatible comparison was possible between the newly suggested classification and the "conventional" EG&K classification. There has been no classification reflecting both eggshell formation and embryo development. Combining of EG&K classification with eggshell formation, formation of area pellucida begins from EG&K stage VII, thus this stage was the first lineage differentiation in chicken. Calcium-deposited eggshell was formed from EG&K stage V and eggshell hardening was observed from EG&K stage VII. There seems to be a close correlation between eggshell formation and formation of area pellucida. By employing this modified classification, it will be feasible to identify and to collect embryos at various intrauterine stages.
We found a significant difference in the dynamics of the sperm that had penetrated into oocytes and in early cleavage. Polyspermic fertilization, with large numbers of decondensed or condensed sperm in an oocyte, was observed. Differing from mammals, many unfertilized supernumerary sperm heads were observed in the yolk area, as well as in the cytoplasm. Several sperm heads in the yolk were undergoing decondensation. The sperm tract from the extracellular space into the yolk was unknown, whether it was direct penetration into the yolk or penetration via the cytoplasm. Decondensed sperm may pass through the cytoplasm during the initial stage of egg development.
Asymmetric cleavage was initiated as early as from the first cleavage, which triggered radiation-oriented progress from central to peripheral part. Central cells in a cleaving embryo seemed to divide very rapidly while peripheral cells, including open cells, divided relatively very slowly. The peripheral furrowing could be readily distinguished from the central one by their length and origin. The peripheral cleavage furrows formed from the peripheral edge of embryo, elongated toward the center, and were more easily visible from the ventral side; however, they were not detectable after EG&K stage IV. This furrowing-type cleavage yielded lots of differences compared with the cleavage of mammalian embryos. In mammals, asymmetric, polarized cleavage signs the initiation of differentiation, while in chick, each part of the preblastoderm  . Mitotic nuclei stained with DAPI were observed before cleavage furrow formation (arrows), and the furrow formed after detection of mitotic nuclei (arrowheads). New cleavage furrows developing between two daughter nuclei were observed from the ventral, rather than the dorsal side, showing completion of diakinesis before cytokinesis. The order of cleavage furrow formation was indicated in Arabic numerals (B). (C, D) Cleavage of EG&K stage I-II embryos were monitored. Multinuclear preblastodermal cells having two daughter nuclei were detected, while due to vigorous proliferation, the size of preblastodermal cells in the cleavage initiation region was smaller than that of the cells in the peripheral region at the initial cleavage stages. Decondensed sperm heads were visible on the ventral side of the embryos (C') and condensed sperm heads were visible on the dorsal side (C") and formation of the large number of cleavage furrows before cytokinesis was visible primarily in the peripheral region (C). Formation of cleavage furrows with mitotic nuclei in stage II was visible (arrowheads in D). (E, F) EG&K stage II-III embryos had many decondensed sperm heads, considered to be penetrated sperm, in the peripheral yolk part. being separated was still connected to each other at the ventral side even after initial furrowing. So, it is difficult to simply reflect the knowledge from the mammals and to further justify the signs of initial differentiation.
Preblastodermal cell divides rapidly. Bellairs et al. [4] stated that the open cells mitotically divide into two daughter cells:  [7] may be an inducible factor for dividing central cells vertically to create two or more layers.
Polyspermy or supernumerary sperm are not common in mammals, whereas they are consistently found in avian species [13]. Chick embryos begin normal development after numerous sperm penetrate the oocyte cell membrane, suggesting that supernumerary sperm may be important to ensure karyogamy [14]. Considering the small area of the germinal disc in relation to the entire ovum of the chicken, polyspermy or supernumerary sperm are necessary to ensure fertilization [13]. Previous reports have shown that low sperm penetration reduces the fertilization rate in chickens [15,16]. Co-localization of decondensed supernumerary sperm in peripheral small cleavage furrows suggested that decondensation of sperm nuclei is a prerequisite for the shortlived supernumerary sperm-associated peripheral cleavage furrows. We found that decondensed sperm were located mainly on the ventral side of the embryos, specifically underneath the cytoplasm, whereas condensed sperm were located mainly on the dorsal side. This might indicate different role of intracytoplasmic, decondensation factors in development of chicken embryos, compared with mammalian embryos.  Table 2. Approximate number of condensed sperm heads on the * dorsal side of EG&K stage I-III embryos after penetration.