  We compute power-suppressed corrections to the \eta\gamma and
\eta^{\prime}\gamma transition form factors
Q^2F_{\eta(\eta^{\prime})\gamma}(Q^2) arising from the end point regions x \to
0,1 by employing the infrared-renormalon approach. The contribution to the form
factors from the quark and gluon content of the \eta,\eta^{\prime} mesons is
taken into account using for the \eta-\eta^{\prime} mixing the SU_f(3) singlet
\eta_1 and octet \eta_8 basis. The theoretical predictions obtained this way
are compared with the corresponding CLEO data and restrictions on the input
parameters (Gegenbauer coefficients) B_2^q(\eta_1), B_2^g(\eta_1), and
B_2^q(\eta_8) in the distribution amplitudes for the \eta_1,\eta_8 states with
one nonasymptotic term are deduced. Comparison is made with the results from
QCD perturbation theory.
