In animal cells, mitotic spindles are oriented by the dynein/dynactin motor complex, which exerts a pulling force on astral microtubules. Dynein/dynactin localization depends on Mud/NUMA, which is typically recruited to the cortex by Pins/LGN. In Drosophila neuroblasts, the Inscuteable/Baz/Par-6/aPKC complex recruits Pins apically to induce vertical spindle orientation, whereas in epithelial cells Dlg recruits Pins laterally to orient the spindle horizontally. Here we investigate division orientation in the Drosophila imaginal wing disc epithelium. Live imaging reveals that spindle angles vary widely during prometaphase and metaphase, and therefore do not reliably predict division orientation. This finding prompted us to re-examine mutants that have been reported to disrupt division orientation in this tissue. Loss of Mud misorients divisions, but Inscuteable expression and aPKC, dlg and pins mutants have no effect. Furthermore, Mud localizes to the apical-lateral cortex of the wing epithelium independently of both Pins and cell cycle stage. Thus, Pins is not required in the wing disc because there are parallel mechanisms for Mud localization and hence spindle orientation, making it a more robust system than in other epithelia.

Although spindle orientation has been extensively examined in asymmetrically dividing cells, less attention has been given to orientation in symmetrically dividing epithelia. As the tissue develops, most epithelial cell divisions are oriented perpendicular to the plane of the tissue so that both daughter cells lie within the epithelial layer (McCaffrey and Macara, 2011). The orientation of division is determined by the orientation of the mitotic spindle. This orientation depends on a conserved pathway that includes Partner of Inscuteable (Pins; GPR-1/2 in C. elegans, LGN or GPSM2 in vertebrates), which anchors Mushroom body defect (Mud; LIN-5 in C. elegans, NuMA or NUMA1 in vertebrates) to the cortex. This pathway is thought to work in every mitotic cell type and organism.

To ensure that new cells are born within the plane of the tissue, mitotic spindles must be oriented orthogonally, along the plane. This means that the spindle-orienting machinery must be lateral at mitosis to pull the two spindle poles into alignment. To date, studies into the regulation of this localization have focused on Pins/LGN. Work in Madin-Darby canine kidney (MDCK) cells and in the Drosophila imaginal disc has suggested that lateral localization of Pins/LGN is regulated by atypical Protein kinase C (aPKC), which excludes it from the apical cortex (Guilgur et al., 2012; Hao et al., 2010). However, this is not the case in the Drosophila follicle epithelium or chick neuroepithelium, in which spindle orientation is aPKC independent (Bergstralh et al., 2013b; Peyre et al., 2011). In these two tissues, the position of the spindle-orienting machinery is determined by the lateral polarity factor Discs large (Dlg), which provides positional information by interacting directly with Pins/LGN (Bergstralh et al., 2013b; Saadaoui et al., 2014). This interaction is mediated by the C-terminal guanylate kinase (GUK) domain in Dlg, which binds a phosphorylated sequence in Pins/LGN (Johnston et al., 2012; Zhu et al., 2011). Binding is thought to be temporally restricted to mitosis by Lgl [L(2)gl – FlyBase], which binds the GUK domain in interphase and is released at mitosis upon phosphorylation by Aurora A/B (Bell et al., 2015; Carvalho et al., 2015). In agreement with this, Lgl variants that cannot be phosphorylated by Aurora A/B remain cortical in mitosis (Bell et al., 2015; Carvalho et al., 2015).

In the current study, we aimed to determine the kinetics of epithelial cell spindle orientation in the Drosophila imaginal disc. These measurements led to the unexpected finding that, unlike other well-characterized epithelia, spindle orientation in this tissue proceeds through a Pins-independent mechanism.

Spindle angles in the disc vary widely

We used fluorescently tagged Centrosomin (Cnn) and Tubulin (α-Tub) to follow spindle orientation and cell division in the pouch region of live third larval instar wing imaginal discs. Mitotic angles were determined by drawing a line between the two centrosomes and measuring the angle of this line relative to the tissue plane (these angles are labeled αz in figures). Since this method allowed us to track angles prior to spindle formation, we began our measurements 1 min before nuclear envelope breakdown (NEBD) and continued until the appearance of the midbody, which marked the first minute of telophase (Fig. 1A,B, Movie 1).

We next compared angles at different phases of mitosis (Fig. 1C). Because the time between spindle formation and anaphase varied between divisions, we normalized this period and divided it into quintiles. This resulted in sample sizes and times (mean n=50, mean duration=2.3 min) comparable to the period between NEBD and the appearance of a complete spindle (n=50, mean duration=2.3 min) and to the period of anaphase (n=39, mean duration=1.9 min). Centrosome angles are close to random in the first period (the start of prometaphase) and become oriented an average of 6.9 min later (the third quintile). From this point onwards, the distribution of angles is not statistically significant between groups. Thus, spindles are oriented in the disc roughly halfway between NEBD and anaphase.

Our finding that initial spindle angles are nearly random prompted us to investigate the process of spindle orientation. The starting orientation of the mitotic spindle is anticipated by the positions of the two centrosomes at NEBD. To determine how these positions are established, we followed centrosome duplication and movement over time using Abnormal spindle (Asp)-GFP, which marks centrosomes throughout the cell cycle. We found that the behavior of centrosomes in the wing disc is consistent with previous observations in vertebrate pseudostratified epithelia (Spear and Erickson, 2012; Strzyz et al., 2015; Tsai et al., 2010). During interphase, the centrosome is localized at the apical cell surface. As the nucleus undergoes apically directed interkinetic nuclear migration, the centrosome moves towards the nucleus and divides (Fig. 1D). Although the centrosomes sometimes migrate to equivalent apical-basal positions on either side of the nucleus, in other cases one centrosome remains apical while the other moves basally to the opposite side of the nucleus. This orientation of the centrosomes often persists until the spindle is formed. We considered whether the temporal variability between mitoses could be explained by the time it takes for these vertical spindles to orient, but did not find a correlation between initial spindle angle and division time.

Divisions orient along the plane of the tissue

We observed a wide variability in spindle angles during the period between NEBD and anaphase, with 58% of spindles having an angle of >30° when they first form and 50% of spindles exceeding this angle for the first half of this period. To determine whether this variability was an artifact of live imaging ex vivo, we quantified spindle angles in five wing discs from three dissections. The cumulative distribution of these angles agreed with measurements made in live tissue (Fig. S1A).

The distribution of angles varied between discs, and one fixed disc showed a significantly different mean spindle angle from the cumulative average (Fig. S1A). This reveals that fixed-tissue measurements are very sensitive to the proportion of early spindles in the sample, which can give misleading results. Phospho-histone H3, a standard marker for mitotic cells, cannot be used to exclude these spindles because it appears before the spindle has formed (Fig. S1B). Another potential confounding factor is that one of the centrosomes in an early mitotic cell often lies closer to a centrosome within the same z-plane in an adjacent cell than it does to the other centrosome in the same cell, making it difficult to reliably assign spindle angles in the absence of a membrane marker (Fig. S1C).

Because of the variability in spindle orientation during metaphase, we restricted our subsequent analysis of division orientation in the wing disc to measurements of post-metaphase cells with separating chromosomes. In agreement with our live imaging measurements, all but one of the angles we measured in fixed tissue was <30° (Fig. S1A). As a positive control for division misorientation, we examined wing discs mutant for the canonical spindle orientation gene mud (Fig. 2A-C) (Bell et al., 2015; Kraut et al., 1996; Nakajima et al., 2013; Wodarz et al., 2000, 1999). As expected, the distribution of division angles in these discs differed significantly from the control, confirming that Mud is an essential component of the spindle orientation machinery in these cells (Bergstralh et al., 2013b; Nakajima et al., 2013; Saadaoui et al., 2014).

The orientation of division does not require aPKC

In MDCK cell cysts, as in most epithelial tissues examined to date, the polarity kinase aPKC localizes along the apical cell cortex, where it has been proposed to play a key role in spindle orientation by phosphorylating LGN (vertebrate Pins) to exclude it from the apical region (Hao et al., 2010). The same mechanism has been implicated in spindle orientation in the Drosophila wing imaginal disc (Guilgur et al., 2012). A drawback to this model is that aPKC is not apical in this tissue, but is instead concentrated at the uppermost part of the lateral cortex, suggesting that it is not in the appropriate position to regulate spindle orientation by excluding Pins (Georgiou et al., 2008; Guilgur et al., 2012). One possibility is that aPKC moves apically during mitosis. However, we observed that aPKC spreads down the lateral cortex at mitosis but remains absent from the apical cortex (Fig. 2D). Similar observations have been made in the pupal notum, which derives from the same imaginal disc as the wing (Rosa et al., 2015).

These observations prompted us to re-examine the role of aPKC in the wing disc. Clones of the genetic null allele aPKCK06403 do not survive, but wing discs can be isolated from larvae transheterozygous for aPKCK06403 and the temperature-sensitive allele aPKCTS (Guilgur et al., 2012; Rolls et al., 2003). Earlier work suggested that spindles are misoriented in these discs at 25°C and higher temperatures (Guilgur et al., 2012). However, the distribution of division angles at anaphase and telophase was normal in discs isolated from these larvae at both 25°C and 29°C (Fig. 2C). Division orientation was also normal in wing discs transheterozygous for aPKCK06403 and the ‘kinase-dead’ allele aPKCPSU141 (Fig. 2C) (Kim et al., 2009). Consistent with reported results using the aPKCTS allele, extensive apoptotic cell death was observed at the basal surface of these discs (not shown), indicating that aPKC function was compromised (Guilgur et al., 2012). These findings show that aPKC does not regulate spindle orientation in the imaginal wing disc, and are consistent with previous studies in the chick neuroepithelium and the Drosophila notum and follicular epithelium (Bergstralh et al., 2013b; Peyre et al., 2011; Rosa et al., 2015). They contrast, however, with work performed in the zebrafish retinal neuroepithelium, since morpholinos targeting aPKCλ/ζ promote division misorientation in that tissue (Cui et al., 2007; Strzyz et al., 2015).

Ectopically expressed Inscuteable does not reorient divisions in the wing disc

In neuroblasts, apically localized Inscuteable recruits the spindle-orienting machinery to the apical cortex. This provides a pulling force that draws one spindle pole proximal to the apical cortex, thereby aligning the mitotic spindle along the apical-basal axis (reviewed by Bergstralh et al., 2013a). In the Drosophila embryonic ectoderm, optic lobe neuroepithelium and follicular epithelium the ectopic expression of Inscuteable performs the same function, resulting in cell divisions that are reoriented by ∼90° relative to the plane of the tissue (Bergstralh et al., 2015; Egger et al., 2007; Kraut et al., 1996). We therefore examined whether Inscuteable also reorients divisions in the wing imaginal disc. Surprisingly, all measured divisions in Inscuteable-expressing wing discs were aligned within 30° of the plane of the epithelium, showing no detectable difference from wild type (Fig. 2C).

We explored the possibility that the failure of Inscuteable to reorient divisions in the wing disc, as it does in other Drosophila epithelia, could be attributed to aPKC. In neuroblasts, Inscuteable is recruited to the apical cortex by aPKC. Conversely, Inscuteable is itself required for the apical localization of aPKC (Wodarz et al., 2000). We found that this cooperative localization also occurs when Inscuteable is ectopically expressed in the follicular epithelium. In wild-type mitotic follicle cells, aPKC loses its apical enrichment, spreading out around the cortex (Bergstralh et al., 2013b; Morais-de-Sá and Sunkel, 2013). In mitotic follicle cells expressing Inscuteable, aPKC remains enriched at the apical cortex, although some aPKC also spreads laterally (Fig. 2E). Thus, Inscuteable and aPKC are mutually required to localize apically in this epithelial cell type during mitosis.

This raises the question of whether the same mechanism works in the wing disc, in which aPKC is lateral rather than apical at interphase. We used hedgehog-Gal4 to drive Inscuteable expression in the posterior compartment of the wing pouch. During mitosis, Inscuteable is not apical, as it is in the follicle epithelium, but instead localizes at the top of the lateral cortex (Fig. 2F). The localization of neither aPKC nor its partner Bazooka differs from that in the wild type (Fig. 2G, Fig. S2).

Taken together, these results indicate that Inscuteable can stabilize, but not localize, aPKC at the apical cortex of an epithelial cell. They also provide one possible explanation for why Inscuteable does not reorient spindles in the wing disc: Inscuteable cannot facilitate pulling of just one spindle pole proximal to the apical cortex, since it is not localized apically but in a lateral belt.

The Dlg/Pins/Lgl pathway does not regulate division orientation in the wing disc

In the Drosophila follicle epithelium and the chick neuroepithelium, Dlg determines the positions of the spindle poles by recruiting Pins to the lateral cortex during mitosis (Bergstralh et al., 2013b; Saadaoui et al., 2014). Previous work identified a spindle orientation defect in wing pouches after Dlg (Dlg1) protein was knocked down using UAS-Dlg-shRNA (TRiP.HMS00024) driven by nubbin-Gal4. In our hands, division angles could not be reliably measured in discs from these larvae grown at the standard temperature (25°C) because the tissue was severely disorganized, as expected from the dlg (dlg1) mutant phenotype (Fig. S3A) (Gateff, 1978; Nakajima et al., 2013). We could reduce, but not eliminate, disorganization by allowing the tissue to first develop at 18°C to decrease the efficiency of the Gal4 system, then transferring the larvae to 25°C for the last 24 h. Despite loss of Dlg protein (as measured by immunostaining), misoriented divisions were not observed in the more organized regions of these discs (Fig. 3A,B). This suggests that the spindle phenotype observed in Dlg knockdown discs might be due to a loss of epithelial polarity and organization rather than a direct effect on spindle orientation per se. Furthermore, although Dlg has been proposed to act cooperatively with Scribble to orient spindles in the disc, divisions were also oriented normally in UAS-Scribble-shRNA (TRiP.HMS01490) wing discs (Fig. 3C, Fig. S3B).

The role of Lgl in division orientation in the disc is also unclear, since previous work suggests that Lgl is unlikely to protect the Dlg GUK domain from binding to Pins in this tissue. First, the affinity of Dlg for Pins/LGN (KD=0.33 µM) is over 30 times greater than its affinity for Lgl (KD=10.2 µM) (Zhu et al., 2011, 2014). Thus, Pins would be expected to simply outcompete Lgl for binding. Second, Dlg is restricted to the top (apical) region of the lateral cortex in the wing disc, whereas Lgl extends further down (Fig. 3D). This suggests that cortical localization of Lgl does not require direct interaction with Dlg, although it does not rule out the possibility that Dlg is required to localize Lgl where they overlap.

These observations prompted us to examine the Dlg/Pins/Lgl pathway directly. In the follicle epithelium, dlg18 is thought to disrupt Pins binding without affecting the essential role that Dlg plays in apical-basal polarity (Bergstralh et al., 2013b). This allele is a nonsense mutation that removes the last 43 residues, comprising roughly one-third of the GUK domain (Woods et al., 1996). This amino acid sequence is highly conserved, suggesting functional importance, but does not include the residues that contact phosphorylated binding partners directly, leaving open the possibility that the truncation does not prevent binding (Fig. S3C). We tested this possibility in vitro using vertebrate Dlg4 and LGN (vertebrate Pins). The purified GUK domain lacking the C-terminal 43 residues is unable to bind phosphorylated LGN, its high-affinity target, confirming that the mutation inactivates the phosphoprotein binding activity of the GUK domain (Fig. S3D).

We generated dlg18/dlg18 mitotic clones in the wing disc and observed that spindle angles within the clones fell within 30° of the plane of the epithelium, as in wild type (Fig. 3B,E) (Bergstralh et al., 2013b). We also found normal localization of Lgl in dlg18/dlg18 clones, confirming that localization does not require an interaction between single-phosphorylated Lgl and the Dlg GUK domain (Fig. S3E). To test whether the removal of Lgl from the cortex is necessary for correct division orientation in the disc, we expressed Lgl-ASA (S656A, S664A), a variant of the protein that remains cortical during mitosis because it cannot be phosphorylated by the Aurora kinases (Bell et al., 2015; Carvalho et al., 2015). When this construct was expressed in clones of the null allele lgl4 using the MARCM technique, division angles did not differ from wild type (Fig. 3B,F). Division orientation was also normal in lgl4/lgl334 cells expressing Lgl-ASA (Fig. S3F). Thus, the removal of Lgl from the cortex is not required for correctly aligned divisions in the wing disc. Taken together, these results demonstrate that the Dlg/Pins/Lgl pathway is not required for spindle orientation in this tissue.

Spindle orientation in the wing disc does not require Pins

Pins/LGN is required to orient spindles in every mitotic cell type examined so far, with the exception of the pupal notum (David et al., 2005). However, we have shown that Dlg, Lgl, aPKC and Inscuteable, which are all proposed to exert their spindle-orienting effects through Pins, are not required for division orientation in the wing disc. This raises the question of whether Pins itself is necessary. In agreement with earlier work, we confirmed that Pins is cortically enriched in dividing cells (Fig. 4A, Fig. S4A). To test its functional importance, we generated clones mutant for pinsp62, a 2112 bp deletion that removes the translation start site (Le Borgne and Schweisguth, 2003; Yu et al., 2000). As measured by immunostaining, there is no detectable Pins protein in pinsp62 mutant clones, confirming that it is a null allele (Fig. 4A). The pinsp62 allele randomizes spindle orientation and misorients divisions in the follicle epithelium (Fig. S4B) (Bergstralh et al., 2013b). By contrast, division orientation is normal in pinsp62/pinsp62 wing disc cells (Fig. 4B,C). Thus, the wing imaginal disc epithelium does not require Pins/LGN to orient spindles.

When Dlg is knocked down in the chick neuroepithelium, spindle angles are randomized during metaphase but largely normal at anaphase, suggesting the existence of a correction pathway (Saadaoui et al., 2014). Although Pins is not necessary for spindle orientation in the wing disc, it might still play a redundant role that is compensated for by an anaphase correction mechanism. To address this possibility, we examined spindle orientation throughout mitosis by making time-lapse movies of divisions in wing discs transheterozygous for pinsp62 and the strong allele pins193, a 2658 bp deletion (Parmentier et al., 2000) (Fig. 4C). As we did for the wild type, we normalized progression through mitosis and examined the distribution of centrosome angles at different stages (Fig. 4D). In all of these periods, the distribution of spindle angles in pinsp62/pins193 discs showed no significant deviation from the wild type (Fig. 4E, Fig. S4C). The duration of spindle formation and orientation (appearance of the spindle to anaphase) in pinsp62/pins193 discs was slightly, but not significantly, longer than that of the wild type (control, 12.6±3.3 min; pins mutant, 15.5±4.3 min). Thus, Pins is not required for spindle orientation in the wing disc. This experiment also addresses the possibility that aPKC, Inscuteable, Dlg and Lgl spindles undergo correction, since all four factors are thought to mediate spindle orientation through Pins.

Mud can localize independently of Pins

In the canonical pathway, Mud/NuMA orients spindles by exerting a pulling force on astral microtubules. This mechanism is likely to be conserved in the wing disc, since this tissue requires both Mud and centrosomes to orient divisions (Bell et al., 2015; Nakajima et al., 2013; Poulton et al., 2014). However, we have shown that the cortical anchor Pins is dispensable. We therefore examined the localization of Mud in the wing disc. Surprisingly, we observed that Mud localizes to discrete cortical foci, both during interphase and mitosis (Fig. 5A,B). This finding agrees with recently published work (Bosveld et al., 2016). To demonstrate antibody specificity, we used hedgehog-Gal4 to drive expression of UAS-Mud-shRNA (TRiP.HMS01458) in the posterior compartment, which abrogated immunoreactivity (Fig. S4C).

In the disc, Mud localization is not affected in pinsp62 mutant clones or in aPKCK06403/aPKCTS transheterozygous mutant wing discs (Fig. 5C, Fig. S4E). When viewed along the apical-basal axis, Mud localizes to the apical region of the lateral cortex, where it only partially overlaps with septate junctions (marked by Dlg). These junctions are proximal to spindle poles in this tissue (Nakajima et al., 2013) (Fig. 5D). Thus, Mud is in the appropriate position to orient spindles in the wing disc by a mechanism that is independent of both cell cycle and Pins.

The possibility that Mud can localize without Pins is not without precedent. In cultured vertebrate cells NuMA can associate with the membrane independently of LGN. This depends on the dephosphorylation of NuMA at a conserved Cdk1 target sequence at anaphase (Kiyomitsu and Cheeseman, 2013; Kotak et al., 2013; Seldin et al., 2013; Zheng et al., 2014). Given that the relevant Cdk1 target sequence is not conserved in Drosophila (not shown), and that Mud localizes throughout the cell cycle in the wing disc, this mechanism is unlikely to be at work in this tissue. We also explored the possibility that Mud interacts directly with Dishevelled (Dsh), as it does in sensory organ precursor (pI) cells (Ségalen et al., 2010). However, we did not observe defective spindle orientation in dsh1/dsh1 mutant wing discs (Fig. S4F,G). Thus, the nature of the mechanism that localizes Mud in the wing imaginal disc remains to be identified.

Pins/LGN (GPR-1/2 in nematodes) is required to orient spindles in almost every instance of mitotic spindle orientation in C. elegans, Drosophila and vertebrate cells studied to date. This includes the symmetrically dividing epithelial cells of the Drosophila follicle epithelium and chick neuroepithelium, suggesting that Pins/LGN is broadly required in epithelial tissues (Bergstralh et al., 2013b; Saadaoui et al., 2014). In support of this, we and two other groups have observed cortical localization of Pins in mitotic cells of the Drosophila imaginal disc (Dewey et al., 2015; Guilgur et al., 2012). Thus, it seemed likely that Pins retains its spindle orienting function in this tissue. Surprisingly, when we tested the role of Pins directly, we found that it is dispensable for spindle orientation in this tissue. This probably continues into the pupal notum, which derives from the same tissue as the larval imaginal disc (David et al., 2005).

These observations contrast with those of other studies, in which both direct evidence (derived from fixed pins mutant tissue) and indirect evidence (derived from fixed tissue lacking the function of putative Pins-regulatory factors) has been used to demonstrate a crucial role for Pins in orienting mitotic spindles in the wing disc. This contradiction is explained by a technical consideration. By imaging spindles in live tissue, we determined that mitotic spindle angles in the wing disc can vary widely over time and are not reliable predictors of division orientation until anaphase. Spindle angle measurements are thus sensitive to the stage of mitosis, which was not accounted for in the earlier work. By restricting our measurements to post-metaphase cell division angles, we determined that neither Pins nor its putative regulatory factors are required to orient divisions in the wing disc.

Pins-independent spindle orientation pathways have been described in three other cases, but none of these seems to be related to spindle orientation in the imaginal wing disc. In sensory organ precursor pI cells in the pupal notum, Mud is recruited to one side of the cell by Pins and to the other side by Dsh (Ségalen et al., 2010). This is unlikely to be relevant, since we have shown that Dsh is dispensable for horizontal orientation in the imaginal disc epithelium. In asymmetrically dividing mouse skin progenitors that divide perpendicularly relative to the plane of the tissue, LGN (vertebrate Pins) cooperates with mInsc (the mammalian ortholog of Inscuteable) to orient spindles along the apical-basal axis (Williams et al., 2014). Surprisingly, LGN is not required for horizontal orientation (Williams et al., 2011). However, this orientation is also independent of NuMA (vertebrate Mud). By contrast, we and others have shown that spindle orientation in the wing disc requires Mud (Nakajima et al., 2013). As mentioned earlier, vertebrate NuMA also has an LGN-independent activity during anaphase (Kiyomitsu and Cheeseman, 2013; Kotak et al., 2013; Seldin et al., 2013; Zheng et al., 2014). This function relies on the dephosphorylation of its C-terminal Cdk1 site, which allows the C-terminal region of the protein to interact with the plasma membrane. However, neither the Cdk1 site nor the C-terminal plasma membrane-binding domain is conserved in Drosophila, and the cortical localization of Mud in the wing disc is not cell cycle regulated.

Nevertheless, our data clearly demonstrate that the cortical localization of Mud in the wing disc does not require Pins. In fact, direct examination revealed that even in the absence of Pins, Mud is enriched in cortical foci throughout the cell cycle. This result explains the finding that Pins is not required to orient spindles, but raises further questions about the mechanism that localizes Mud. After submission of our manuscript, another group published that Mud foci in the Drosophila notum and wing disc correspond to tricellular junctions (Bosveld et al., 2016). These specialized structures, characterized by distinct protein components including Gliotactin and Anakonda (Bark beetle – FlyBase), form in epithelial tissues with mature septate junctions (Byri et al., 2015; Schulte et al., 2006). Since they are located at lateral cell-cell contacts, they fulfill the requirement for the location of the cortical pulling force that drives horizontal spindle orientation in epithelial tissues.

It is unclear why the wing imaginal disc has evolved a Pins-independent mechanism to orient mitoses. One possibility is that this is related to the fact that spindle orientation is essential to maintain cells in the epithelial layer in this tissue. Misplaced daughter cells in the disc undergo apoptosis and are extruded basally, whereas other epithelia can compensate for misaligned divisions by simply reintegrating the misplaced cells (Bergstralh et al., 2015; Nakajima et al., 2013). These other epithelia have immature septate junctions, and the lateral adhesion proteins that drive reintegration, such as Neuroglian and Fasciclin 2, localize along the entire length of the lateral cortex. In the wing disc, by contrast, the lateral adhesion proteins are tightly restricted to the septate junctions in the apical region of the lateral cortex. This means that they are not in the correct position to adhere to cells that have been basally displaced by misoriented divisions and therefore cannot drive their reintegration. The presence of mature septate junctions in the wing disc also means that this tissue has tricellular junctions, unlike the epithelia in which reintegration occurs. We suggest that in order to compensate for its inability to reintegrate misplaced cells, the wing disc has taken advantage of its tricellular junctions to provide a robust backup mechanism for localizing Mud to the lateral cortex, and thus for spindle orientation.

Drosophila mutants

The following mutant alleles and transgenic constructs have been described previously: mud3 and mud4 (Yu et al., 2006), pins193 (Parmentier et al., 2000), pinsp62 (Yu et al., 2000), dlg18 (Woods and Bryant, 1989), aPKCk06403 (Wodarz et al., 2000), apkcts (aPKCTS) (Guilgur et al., 2012), lgl4 (Ohshiro et al., 2000), dsh1 (Perrimon and Mahowald, 1987), UAS-Lgl-ASA-GFP (Bell et al., 2015), UAS-Inscuteable (Kraut et al., 1996), nubbin-Gal4 (Thompson and Cohen, 2006) and hedgehog-Gal4 (Tanimoto et al., 2000). dlg18 FRT19A (a gift from Floris Bosveld, Institut Curie, Paris, France) and FRT82B pinsp62 were described previously (Bergstralh et al., 2013b). The following background stocks were used to generate mitotic clones, which were induced by heat shock at 37°C for multiple periods of 2 h: RFP-nls, hsflp, FRT19A, hsflp;FRT40A RFP-nls, hsflp;; FRT82B RFP-nls and hsflp;; FRT82B GFP-nls. Mosaic analysis with a repressible cell marker (MARCM; after the method of Lee and Luo, 1999) was carried out using GFP-mCD8 (under the control of an actin promoter) as the marker. The background stock was generated by Aram Sayadian (Gurdon Institute, Cambridge, UK). We thank the Transgenic RNAi Project at Harvard Medical School [NIH/NIGMS R01-GM084947] for providing UAS-Dlg-shRNA (HMS00024), UAS-Scribble-shRNA (HMS01490) and UAS-Mud-shRNA (HMS01458).

Fluorescent marker stocks

We used the following fluorescent markers: Dlg::YFP (Bergstralh et al., 2013b), Ubi-Asp-GFP (Rujano et al., 2013), Ubi-Cnn-RFP (Basto et al., 2008), Ubi-Cnn-GFP (Conduit et al., 2010), Ubi-α-Tub-RFP (Basto et al., 2008) and Ubi-α-Tub84B-GFP (Rebollo et al., 2004).

Reagents

The following antibodies were used: rabbit anti-Centrosomin (gift from J. Raff; Lucas and Raff, 2007), rabbit anti-Inscuteable (gift from J. Knoblich; Kraut et al., 1996), rabbit anti-phospho-histone H3 (Ser10, Cell Signaling, #9701, lot #13), rabbit anti-Bazooka (gift from A. Wodarz; Wodarz et al., 2000), rabbit anti-Mud (gift from R. Basto; Rujano et al., 2013), rabbit anti-Pins (gift from F. Matsuzaki; Izumi et al., 2006), rabbit anti-aPKC and anti-Lgl (Santa Cruz, sc-27509, dN-16, lot #H3107), mouse anti-Dlg (Developmental Studies Hybridoma Bank, clone 4F3, 6/5/14), and mouse FITC-conjugated anti-α-Tubulin (Sigma, clone DM1A, lot #114M4817V). Cy5- and Texas Red-conjugated secondary antibodies were purchased from Jackson ImmunoResearch. Phalloidin was purchased from Invitrogen and Vectashield with DAPI was purchased from Vector Labs. Primary and secondary antibodies were used at a dilution of 1:150.

Imaging

Immunofluorescence and fixed cell imaging were performed as previously described (Bergstralh et al., 2013b). Live imaging was performed using a Leica SP5 (63×/1.4 HCX PL Apo CS oil). z-stacks of planes spaced 0.5 µm apart were taken at 1 min intervals. Wing discs were dissected and imaged in 0.8% agarose in Schneider's medium (Sigma) containing 10 µg/ml insulin (Sigma). Images were collected with Leica Application Suite AF and processed (Gaussian blur) using ImageJ (NIH).

Spindle angle measurements

Centrosome angles were calculated using ImageJ. Angles were determined by drawing a first line connecting the two spindle poles and a second line along the apical surface of the tissue, then measuring the angle between them. These measurements frequently required correction in the xy plane, such that both spindle poles were apparent in a single z-plane. Statistical analyses were performed using Prism (GraphPad). To prevent biased distribution of fixed images, we counted all angles in each wing disc examined. For analysis of live images, we counted all complete divisions within the 1-h window imaged. Images were analyzed by four independent researchers. As described, centrosome angles were measured from the minute preceding NEBD until the first minute of telophase. NEBD was negatively marked by tubulin, which is clearly excluded from the nucleus prior to NEBD. Telophase was marked by the midbody, which is distinguished by its size, morphology and position.

Test for correlation

We used the Pearson product-moment correlation to test for correlation between initial spindle angle and the duration of mitosis. We obtained an r value of 0.2275, which is less than the critical value of 0.360 for a sample size of 22 (degrees of freedom=20).

Isothermal titration calorimetry measurements

ITC measurement was performed on an ITC200 microcalorimeter (MicroCal) at 25°C. Protein and peptide samples were dissolved in buffer containing 50 mM Tris pH 8.0, 100 mM NaCl and 1 mM EDTA. The protein concentrations used in the cell (GK mutant) and in the syringe (phospho-LGN peptide) were 0.05 mM and 0.48 mM, respectively. The titration was carried out at time intervals of 2 min to ensure that the titration peak returned to the baseline. The titration data were analyzed using Origin 7.0 software (MicroCal).

We thank Floris Bosveld and Yohanns Bellaïche for helpful discussions over the course of our study; Jürgen Knoblich, Barry Thompson, Rui Martinho and their labs for fly stocks; and Iwo Kuciński, Joanna Kosałka, Judy Sayers and the D.S.J. lab for technical help and criticism.

Author contributions

D.T.B. and D.S.J. conceived the project. D.T.B., J.Z. and D.S.J. designed the experiments. D.T.B., H.E.L., I.K., N.S.D. and S.C. performed the Drosophila experiments. J.Z. and R.Z. performed the biochemistry. D.T.B. and D.S.J. wrote the manuscript.

Funding

This work was supported by a Wellcome Trust Principal Fellowship to D.S.J. [080007] and by core support from the Wellcome Trust [092096] and Cancer Research UK [A14492]. D.T.B. was supported by a Marie Curie Fellowship (funded by the European Commission) and the Wellcome Trust. H.E.L. was supported by a Herchel Smith Studentship, University of Cambridge. Deposited in PMC for immediate release.

Basto
,
R.
,
Brunk
,
K.
,
Vinadogrova
,
T.
,
Peel
,
N.
,
Franz
,
A.
,
Khodjakov
,
A.
and
Raff
,
J. W.
(
2008
).
Centrosome amplification can initiate tumorigenesis in flies
.
Cell
133
,
1032
-
1042
.
Bell
,
G. P.
,
Fletcher
,
G. C.
,
Brain
,
R.
and
Thompson
,
B. J.
(
2015
).
Aurora kinases phosphorylate lgl to induce mitotic spindle orientation in Drosophila epithelia
.
Curr. Biol.
25
,
61
-
68
.
Bergstralh
,
D. T.
,
Haack
,
T.
and
St Johnston
,
D.
(
2013a
).
Epithelial polarity and spindle orientation: intersecting pathways
.
Philos. Trans. R. Soc. B Biol. Sci.
368
,
20130291
.
Bergstralh
,
D. T.
,
Lovegrove
,
H. E.
and
St Johnston
,
D.
(
2013b
).
Discs large links spindle orientation to apical-Basal polarity in Drosophila epithelia
.
Curr. Biol.
23
,
1707
-
1712
.
Bergstralh
,
D. T.
,
Lovegrove
,
H. E.
and
St Johnston
,
D.
(
2015
).
Lateral adhesion drives reintegration of misplaced cells into epithelial monolayers
.
Nat. Cell Biol.
17
,
1497
-
1503
.
Bosveld
,
F.
,
Markova
,
O.
,
Guirao
,
B.
,
Martin
,
C.
,
Wang
,
Z.
,
Pierre
,
A.
,
Balakireva
,
M.
,
Gaugue
,
I.
,
Ainslie
,
A.
,
Christophorou
,
N.
, et al. 
(
2016
).
Epithelial tricellular junctions act as interphase cell shape sensors to orient mitosis
.
Nature
530
,
495
-
498
.
Byri
,
S.
,
Misra
,
T.
,
Syed
,
Z. A.
,
Bätz
,
T.
,
Shah
,
J.
,
Boril
,
L.
,
Glashauser
,
J.
,
Aegerter-Wilmsen
,
T.
,
Matzat
,
T.
,
Moussian
,
B.
, et al. 
(
2015
).
The triple-repeat protein Anakonda controls epithelial tricellular junction formation in Drosophila
.
Dev. Cell
33
,
535
-
548
.
Carvalho
,
C. A.
,
Moreira
,
S.
,
Ventura
,
G.
,
Sunkel
,
C. E.
and
Morais-de-Sá
,
E.
(
2015
).
Aurora A triggers Lgl cortical release during symmetric division to control planar spindle orientation
.
Curr. Biol.
25
,
53
-
60
.
Conduit
,
P. T.
,
Brunk
,
K.
,
Dobbelaere
,
J.
,
Dix
,
C. I.
,
Lucas
,
E. P.
and
Raff
,
J. W.
(
2010
).
Centrioles regulate centrosome size by controlling the rate of Cnn incorporation into the PCM
.
Curr. Biol.
20
,
2178
-
2186
.
Cui
,
S.
,
Otten
,
C.
,
Rohr
,
S.
,
Abdelilah-Seyfried
,
S.
and
Link
,
B. A.
(
2007
).
Analysis of aPKClambda and aPKCzeta reveals multiple and redundant functions during vertebrate retinogenesis
.
Mol. Cell. Neurosci.
34
,
431
-
444
.
David
,
N. B.
,
Martin
,
C. A.
,
Segalen
,
M.
,
Rosenfeld
,
F.
,
Schweisguth
,
F.
and
Bellaïche
,
Y.
(
2005
).
Drosophila Ric-8 regulates Galphai cortical localization to promote Galphai-dependent planar orientation of the mitotic spindle during asymmetric cell division
.
Nat. Cell Biol.
7
,
1083
-
1090
.
Dewey
,
E. B.
,
Sanchez
,
D.
and
Johnston
,
C. A.
(
2015
).
Warts phosphorylates Mud to promote pins-mediated mitotic spindle orientation in Drosophila, independent of Yorkie
.
Curr. Biol.
25
,
2751
-
2762
.
Egger
,
B.
,
Boone
,
J. Q.
,
Stevens
,
N. R.
,
Brand
,
A. H.
and
Doe
,
C. Q.
(
2007
).
Regulation of spindle orientation and neural stem cell fate in the Drosophila optic lobe
.
Neural Dev.
2
,
1
.
Gateff
,
E.
(
1978
).
Malignant neoplasms of genetic origin in Drosophila melanogaster
.
Science
200
,
1448
-
1459
.
Georgiou
,
M.
,
Marinari
,
E.
,
Burden
,
J.
and
Baum
,
B.
(
2008
).
Cdc42, Par6, and aPKC regulate Arp2/3-mediated endocytosis to control local adherens junction stability
.
Curr. Biol.
18
,
1631
-
1638
.
Guilgur
,
L. G.
,
Prudencio
,
P.
,
Ferreira
,
T.
,
Pimenta-Marques
,
A. R.
and
Martinho
,
R. G.
(
2012
).
Drosophila aPKC is required for mitotic spindle orientation during symmetric division of epithelial cells
.
Development
139
,
503
-
513
.
Hao
,
Y.
,
Du
,
Q.
,
Chen
,
X.
,
Zheng
,
Z.
,
Balsbaugh
,
J. L.
,
Maitra
,
S.
,
Shabanowitz
,
J.
,
Hunt
,
D. F.
and
Macara
,
I. G.
(
2010
).
Par3 controls epithelial spindle orientation by aPKC-mediated phosphorylation of apical Pins
.
Curr. Biol.
20
,
1809
-
1818
.
Izumi
,
Y.
,
Ohta
,
N.
,
Hisata
,
K.
,
Raabe
,
T.
and
Matsuzaki
,
F.
(
2006
).
Drosophila Pins-binding protein Mud regulates spindle-polarity coupling and centrosome organization
.
Nat. Cell Biol.
8
,
586
-
593
.
Johnston
,
C. A.
,
Doe
,
C. Q.
and
Prehoda
,
K. E.
(
2012
).
Structure of an enzyme-derived phosphoprotein recognition domain
.
PLoS ONE
7
,
e36014
.
Kim
,
S.
,
Gailite
,
I.
,
Moussian
,
B.
,
Luschnig
,
S.
,
Goette
,
M.
,
Fricke
,
K.
,
Honemann-Capito
,
M.
,
Grubmüller
,
H.
and
Wodarz
,
A.
(
2009
).
Kinase-activity-independent functions of atypical protein kinase C in Drosophila
.
J. Cell Sci.
122
,
3759
-
3771
.
Kiyomitsu
,
T.
and
Cheeseman
,
I. M.
(
2013
).
Cortical dynein and asymmetric membrane Elongation coordinately position the spindle in Anaphase
.
Cell
154
,
391
-
402
.
Kotak
,
S.
,
Busso
,
C.
and
Gönczy
,
P.
(
2013
).
NuMA phosphorylation by CDK1 couples mitotic progression with cortical dynein function
.
EMBO J.
32
,
2517
-
2529
.
Kraut
,
R.
,
Chia
,
W.
,
Jan
,
L. Y.
,
Jan
,
Y. N.
and
Knoblich
,
J. A.
(
1996
).
Role of inscuteable in orienting asymmetric cell divisions in Drosophila
.
Nature
383
,
50
-
55
.
Le Borgne
,
R.
and
Schweisguth
,
F.
(
2003
).
Unequal segregation of Neuralized biases Notch activation during asymmetric cell division
.
Dev. Cell
5
,
139
-
148
.
Lee
,
T.
and
Luo
,
L.
(
1999
).
Mosaic analysis with a repressible cell marker for studies of gene function in neuronal morphogenesis
.
Neuron
22
,
451
-
461
.
Lucas
,
E. P.
and
Raff
,
J. W.
(
2007
).
Maintaining the proper connection between the centrioles and the pericentriolar matrix requires Drosophila centrosomin
.
J. Cell Biol.
178
,
725
-
732
.
McCaffrey
,
L. M.
and
Macara
,
I. G.
(
2011
).
Epithelial organization, cell polarity and tumorigenesis
.
Trends Cell Biol.
21
,
727
-
735
.
Morais-de-Sá
,
E.
and
Sunkel
,
C.
(
2013
).
Adherens junctions determine the apical position of the midbody during follicular epithelial cell division
.
EMBO Rep.
14
,
696
-
703
.
Nakajima
,
Y.-I.
,
Meyer
,
E. J.
,
Kroesen
,
A.
,
McKinney
,
S. A.
and
Gibson
,
M. C.
(
2013
).
Epithelial junctions maintain tissue architecture by directing planar spindle orientation
.
Nature
500
,
359
-
362
.
Ohshiro
,
T.
,
Yagami
,
T.
,
Zhang
,
C.
and
Matsuzaki
,
F.
(
2000
).
Role of cortical tumour-suppressor proteins in asymmetric division of Drosophila neuroblast
.
Nature
408
,
593
-
596
.
Parmentier
,
M. L.
,
Woods
,
D.
,
Greig
,
S.
,
Phan
,
P. G.
,
Radovic
,
A.
,
Bryant
,
P.
and
O'Kane
,
C. J.
(
2000
).
Rapsynoid/partner of inscuteable controls asymmetric division of larval neuroblasts in Drosophila
.
J. Neurosci.
20
,
RC84
.
Perrimon
,
N.
and
Mahowald
,
A. P.
(
1987
).
Multiple functions of segment polarity genes in Drosophila
.
Dev. Biol.
119
,
587
-
600
.
Peyre
,
E.
,
Jaouen
,
F.
,
Saadaoui
,
M.
,
Haren
,
L.
,
Merdes
,
A.
,
Durbec
,
P.
and
Morin
,
X.
(
2011
).
A lateral belt of cortical LGN and NuMA guides mitotic spindle movements and planar division in neuroepithelial cells
.
J. Cell Biol.
193
,
141
-
154
.
Poulton
,
J. S.
,
Cuningham
,
J. C.
and
Peifer
,
M.
(
2014
).
Acentrosomal Drosophila epithelial cells exhibit abnormal cell division, leading to cell death and compensatory proliferation
.
Dev. Cell
30
,
731
-
745
.
Rebollo
,
E.
,
Llamazares
,
S.
,
Reina
,
J.
and
Gonzalez
,
C.
(
2004
).
Contribution of noncentrosomal microtubules to spindle assembly in Drosophila spermatocytes
.
PLoS Biol.
2
,
e8
.
Rolls
,
M. M.
,
Albertson
,
R.
,
Shih
,
H.-P.
,
Lee
,
C.-Y.
and
Doe
,
C. Q.
(
2003
).
Drosophila aPKC regulates cell polarity and cell proliferation in neuroblasts and epithelia
.
J. Cell Biol.
163
,
1089
-
1098
.
Rosa
,
A.
,
Vlassaks
,
E.
,
Pichaud
,
F.
and
Baum
,
B.
(
2015
).
Ect2/Pbl acts via Rho and polarity proteins to direct the assembly of an isotropic actomyosin cortex upon mitotic entry
.
Dev. Cell
32
,
604
-
616
.
Rujano
,
M. A.
,
Sanchez-Pulido
,
L.
,
Pennetier
,
C.
,
le Dez
,
G.
and
Basto
,
R.
(
2013
).
The microcephaly protein Asp regulates neuroepithelium morphogenesis by controlling the spatial distribution of myosin II
.
Nat. Cell Biol.
15
,
1294
-
1306
.
Saadaoui
,
M.
,
Machicoane
,
M.
,
di Pietro
,
F.
,
Etoc
,
F.
,
Echard
,
A.
and
Morin
,
X.
(
2014
).
Dlg1 controls planar spindle orientation in the neuroepithelium through direct interaction with LGN
.
J. Cell Biol.
206
,
707
-
717
.
Schulte
,
J.
,
Charish
,
K.
,
Que
,
J.
,
Ravn
,
S.
,
MacKinnon
,
C.
and
Auld
,
V. J.
(
2006
).
Gliotactin and Discs large form a protein complex at the tricellular junction of polarized epithelial cells in Drosophila
.
J. Cell Sci.
119
,
4391
-
4401
.
Ségalen
,
M.
,
Johnston
,
C. A.
,
Martin
,
C. A.
,
Dumortier
,
J. G.
,
Prehoda
,
K. E.
,
David
,
N. B.
,
Doe
,
C. Q.
and
Bellaïche
,
Y.
(
2010
).
The Fz-Dsh planar cell polarity pathway induces oriented cell division via Mud/NuMA in Drosophila and zebrafish
.
Dev. Cell
19
,
740
-
752
.
Seldin
,
L.
,
Poulson
,
N. D.
,
Foote
,
H. P.
and
Lechler
,
T.
(
2013
).
NuMA localization, stability, and function in spindle orientation involve 4.1 and Cdk1 interactions
.
Mol. Biol. Cell
24
,
3651
-
3662
.
Spear
,
P. C.
and
Erickson
,
C. A.
(
2012
).
Apical movement during interkinetic nuclear migration is a two-step process
.
Dev. Biol.
370
,
33
-
41
.
Strzyz
,
P. J.
,
Lee
,
H. O.
,
Sidhaye
,
J.
,
Weber
,
I. P.
,
Leung
,
L. C.
and
Norden
,
C.
(
2015
).
Interkinetic nuclear migration is centrosome independent and ensures apical cell division to maintain tissue integrity
.
Dev. Cell
32
,
203
-
219
.
Tanimoto
,
H.
,
Itoh
,
S.
,
ten Dijke
,
P.
and
Tabata
,
T.
(
2000
).
Hedgehog creates a gradient of DPP activity in Drosophila wing imaginal discs
.
Mol. Cell
5
,
59
-
71
.
Thompson
,
B. J.
and
Cohen
,
S. M.
(
2006
).
The Hippo pathway regulates the bantam microRNA to control cell proliferation and apoptosis in Drosophila
.
Cell
126
,
767
-
774
.
Tsai
,
J.-W.
,
Lian
,
W.-N.
,
Kemal
,
S.
,
Kriegstein
,
A. R.
and
Vallee
,
R. B.
(
2010
).
Kinesin 3 and cytoplasmic dynein mediate interkinetic nuclear migration in neural stem cells
.
Nat. Neurosci.
13
,
1463
-
1471
.
Williams
,
S. E.
,
Beronja
,
S.
,
Pasolli
,
H. A.
and
Fuchs
,
E.
(
2011
).
Asymmetric cell divisions promote Notch-dependent epidermal differentiation
.
Nature
470
,
353
-
358
.
Williams
,
S. E.
,
Ratliff
,
L. A.
,
Postiglione
,
M. P.
,
Knoblich
,
J. A.
and
Fuchs
,
E.
(
2014
).
Par3–mInsc and Gαi3 cooperate to promote oriented epidermal cell divisions through LGN
.
Nat. Cell Biol.
16
,
758
-
769
.
Wodarz
,
A.
,
Ramrath
,
A.
,
Kuchinke
,
U.
and
Knust
,
E.
(
1999
).
Bazooka provides an apical cue for Inscuteable localization in Drosophila neuroblasts
.
Nature
402
,
544
-
547
.
Wodarz
,
A.
,
Ramrath
,
A.
,
Grimm
,
A.
and
Knust
,
E.
(
2000
).
Drosophila atypical protein kinase C associates with Bazooka and controls polarity of epithelia and neuroblasts
.
J. Cell Biol.
150
,
1361
-
1374
.
Woods
,
D. F.
and
Bryant
,
P. J.
(
1989
).
Molecular cloning of the lethal(1)discs large-1 oncogene of Drosophila
.
Dev. Biol.
134
,
222
-
235
.
Woods
,
D. F.
,
Hough
,
C.
,
Peel
,
D.
,
Callaini
,
G.
and
Bryant
,
P. J.
(
1996
).
Dlg protein is required for junction structure, cell polarity, and proliferation control in Drosophila epithelia
.
J. Cell Biol.
134
,
1469
-
1482
.
Yu
,
F.
,
Morin
,
X.
,
Cai
,
Y.
,
Yang
,
X.
and
Chia
,
W.
(
2000
).
Analysis of partner of inscuteable, a novel player of Drosophila asymmetric divisions, reveals two distinct steps in inscuteable apical localization
.
Cell
100
,
399
-
409
.
Yu
,
J. X.
,
Guan
,
Z.
and
Nash
,
H. A.
(
2006
).
The mushroom body defect gene product is an essential component of the meiosis II spindle apparatus in Drosophila oocytes
.
Genetics
173
,
243
-
253
.
Zheng
,
Z.
,
Wan
,
Q.
,
Meixiong
,
G.
and
Du
,
Q.
(
2014
).
Cell cycle-regulated membrane binding of NuMA contributes to efficient anaphase chromosome separation
.
Mol. Biol. Cell
25
,
606
-
619
.
Zhu
,
J.
,
Shang
,
Y.
,
Xia
,
C.
,
Wang
,
W.
,
Wen
,
W.
and
Zhang
,
M.
(
2011
).
Guanylate kinase domains of the MAGUK family scaffold proteins as specific phospho-protein-binding modules
.
EMBO J.
30
,
4986
-
4997
.
Zhu
,
J.
,
Shang
,
Y.
,
Wan
,
Q.
,
Xia
,
Y.
,
Chen
,
J.
,
Du
,
Q.
and
Zhang
,
M.
(
2014
).
Phosphorylation-dependent interaction between tumor suppressors Dlg and Lgl
.
Cell Res.
24
,
451
-
463
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information