Question: hypergeometricsols bug

I believe I have found a bug or unintended behaviour:

The following second-order DE:

mplinput := ((x^4 + 4*x^3 + 10*x^2 + 4*x + 1)*diff(y(x), [x $ 0]))/(4*x^2*(1 + x)^2*(x^2 + 6*x + 1)) + 2*(x^2 + 5*x + 2)*diff(y(x), [x $ 1])/((1 + x)*(x^2 + 6*x + 1)) + diff(y(x), [x $ 2])=0:

can be solved by

with(DEtools):
hypergeometricsols(mplinput);

in terms of hypergeom functions. However,

dsolve(mplinput, [hypergeometricsols]);

does not find a solution, although it is supposed to use the same algorithm. One can trace this back to the sub-algorithm 'DEtools/2F1: --- find_2f1 ---' which is apparantly used by hypergeometricsols, but not by dsolve.

Another thing is that dsolve(mplinput, [hypergeometricsols]) in fact works in Maple2020, but does not in Maple2021 and Maple2022.