I will present work that aim at understanding the failure of analytic hypoellipticity of special differential operators, namely sums of squares of (analytic) vector fields. According to Hörmander, given such an operator P that satisfies the ``bracket condition'' and given a smooth function f, if u is distribution solution to Pu=f then u is also a smooth function. P is then said to be hypoelliptic. But if f is real analytic, then u need not to be analytic but merely smooth Gevrey for some indices, usually to be guessed. Examples were built by Metivier, Matsuzawa, Bove, Baouendi-Goulaouic..., using real variables methods. In this joint work with Paulo Cordaro, by using methods of complex analysis, we show that this failure of analytic hypoellipticity due to the presence of irregular singularity of some holomorphic ODEs the analysis of which defines the best Gevrey indices to be expected. The theory of summability of formal solutions of holomorphic ODEs as developped by Ramis, Malgrange, Braaksma is a fundemental tool here.