We find a new duality for form factors of lightlike Wilson loops in planar N = 4 super-Yang-Mills theory. The duality maps a form factor involving an n-sided lightlike polygonal super-Wilson loop together with m external on-shell states, to the same type of object but with the edges of the Wilson loop and the external states swapping roles. This relation can essentially be seen graphically in Lorentz harmonic chiral (LHC) superspace where it is equivalent to planar graph duality. However there are some crucial subtleties with the cancellation of spurious poles due to the gauge fixing. They are resolved by finding the correct formulation of the Wilson loop and by careful analytic continuation from Minkowski to Euclidean space. We illustrate all of these subtleties explicitly in the simplest non-trivial NMHV-like case. and the identification of coordinates with momenta (1.3) to their supersymmetric analogs. The super-Wilson loop form factors are considered in the planar limit and in the lowestorder perturbative approximation (Born level). The introduction of Grassmann variables (θ i on the Wilson loop contour and η j for the on-shell states) allows us to probe the duality for more complicated configurations of particle helicities. By analogy with the amplitudes, we call the contributions at the lowest level in the Grassmann expansion MHV-like, at the next level NMHV-like, etc. At MHV level we confirm the result of [20]. The NMHV level is much more complicated, the form factor being a non-trivial rational function of the kinematical data. Yet, we show that the duality still works, in a rather simple and suggestive way, by just matching planar Feynman diagrams. This allows us to argue that it should hold for the complete supersymmetric object (at all Grassmann levels) and also beyond the Born approximation.The key to understanding the duality is the appropriate superspace formulation of the Wilson loop and its form factor. In the conventional approach the chiral supersymmetric Wilson loop [4,5,6] is formulated in terms of constrained on-shell super-connections [21,22], which makes the Feynman diagram technique highly inefficient. In this paper we prefer to use the Lorentz harmonic chiral (LHC) superspace approach [23]. It provides an offshell formulation of the chiral N = 4 SYM theory in terms of unconstrained prepotentials, best suited for supersymmetric quantisation. LHC superspace is an alternative to the twistor formulation [24,25], closer in spirit to traditional field theory (see also [26]). The main idea is to consider the interacting theory as a perturbation of the self-dual sector. The twistor formulation has been successfully used to justify the so-called MHV rules for the computation of the amplitude [27], to prove the duality between supersymmetric Wilson loops and amplitudes [5], to compute off-shell correlation functions of the N = 4 stress-tensor multiplet [28]. More recently, the LHC formalism was applied to finding the non-chiral completion of the correlators [29] and to the calculation of form facto...