Abstract: Given parameters x∈/R<sup>−</sup>∪1 and ν, $\mathrm{Re}(\nu) < 0$, and the space H<em>0 of entire functions in C vanishing at $0$, we consider the family of operators L=c0​⋅δ∘M with constant c0​=ν(1−ν)x/(1−x), δ=zd/dz and integral operator M defined by Mf(z)=∫0​<sup>1</sup>e<sup>−</sup>xz​t<sup>−ν(1−(1−x)t)</sup>f(xz​t<sup>−ν(1−t)</sup>)tdt​,z∈C, for all f∈H0​. Inverting L or M proves equivalent to solve a singular Volterra equation of the first kind. The inversion of operator L on H0​ leads us to derive a new class of linear inversion formulas T=A(x,ν)⋅S⇔S=B(x,ν)⋅T between sequences S=(Sn​)</em>n∈N<sup>∗ and T=(Tn​)<em>n∈N<sup>∗, where the infinite lower-triangular matrix A(x,ν) and its inverse B(x,ν) involve Hypergeometric polynomials F(⋅), namely $$ \left{ \begin{array}{ll} A</em>{n,k}(x,\nu) = \displaystyle (-1)<sup>k\binom{n}{k}F(k-n,-n\nu;-n;x),</sup> B_{n,k}(x,\nu) = \displaystyle (-1)<sup>k\binom{n}{k}F(k-n,k\nu;k;x)</sup> \end{array} \right. $$ for 1⩽k⩽n. Functional relations between the ordinary (resp. exponential) generating functions of the related sequences S and T are also given. These relations finally enable us to derive the integral representation L<sup>−1f(z)</sup>=2iπx1−x​e<sup>z</sup>∫(0+)​<sup>1</sup>t(1−t)e<sup>−xtz​</sup>f(xz(−t)<sup>ν(1−t)<sup>1−ν</sup></sup>)dt,z∈C, for the inverse L<sup>−1 of operator L on H0​, where the integration contour encircles the point 0.