We study the strong coupling behaviour of null polygonal Wilson loops/gluon amplitudes in N = 4 SYM, by using the OPE series and its integrability features. For the hexagon we disentangle the SU (4) matrix structure of the form factors for fermions, organising them in a pattern similar to the Young diagrams used previously for the scalar sector [1,2]. Then, we complete and extend the discussion of [3] by showing, at strong coupling, the appearance of a new effective particle in the series: the fermion-antifermion bound state, the so-called meson. We discuss its interactions in the OPE series with itself by forming (effective) bound states and with the gluons and bound states of them. These lead the OPE series to the known AdS 5 minimal area result for the Wls, described in terms of a set of TBA-like equations. This approach allows us to detect all the one-loop contributions and, once the meson has formed, applies to N = 2 Nekrasov partition function via the parallel meson/instanton (in particular, they share the mechanism by which their bound states emerge and form the TBA node). Finally, to complete the strong coupling analysis, we consider the scalar sector for any polygon, confirming the emergence of a leading contribution from the non-perturbative theory on the sphere S 5In the realm of the supersymmetric gauge theories a special role is played by N = 4 Super Yang-Mills (SYM), with gauge group SU (N c ) and dimensionless coupling constant g Y M . The theory indeed appears at one side of the most known example of AdS/CFT correspondence [4,5,6], i.e. the duality between type IIB superstring theory on AdS 5 × S 5 and N = 4 SYM on the boundary of AdS 5 , the 4d Minkowski spacetime.Of particular interest is the planar limit N c → ∞, in which we also send g Y M → 0 keeping the 't Hooft coupling λ ≡ N c g 2 Y M fixed. The new coupling λ is the only parameter of the planar theory and in literature it is often represented as λ = 16π 2 g 2 . The planar N = 4 SYM shows remarkable connections with 1+1 dimensional integrable models [7,8,9,10,11,12,13,14,15,16,17,18], which have allowed a better comprehension and partial proof of the aforementioned correspondence.Historically, integrability was discovered when dealing with the spectral problem, i.e. the computation of the anomalous dimension of local gauge invariant operators. More recently, it played an important role also in the evaluation of null polygonal Wilson loops (Wls). They have been proposed to be dual to gluon scattering amplitudes [19,20,21], which makes them even more interesting. Currently, this correspondence has been widely tested both at weak and strong coupling, so that it can be considered a well-established fact.In any conformal quantum field theory, the Operator Product Expansion (OPE) technique can be applied, besides to the usual product of local operators, also to the null polygonal Wilson loops [22]. This method has an intrinsic non-perturbative nature and recalls the Form Factor (FF) Infra-Red (IR) spectral series of the correlation functions in ...