Gather the whole integrand into a RowBox
:
expr = With[{f = Sin[Pi Sqrt[x^2 + y^2]]/Sqrt[x^2 + y^2],
reg = (1 <= x^2 + y^2 <= 4)},
DisplayForm@
InterpretationBox[
RowBox[{UnderscriptBox[StyleBox["∬", FontSize -> 24], reg],
RowBox[{f, RowBox[{"\[DifferentialD]", x}],
RowBox[{"\[DifferentialD]", y}]}]}],
Integrate[f, Element[{x, y}, ImplicitRegion[reg, {x, y}]]]]]
expr // TeXForm