By using Mathematica
Series[MeijerG[{{-(1/2)}, {}}, {{-(1/2), -(1/2), 1/2}, {}}, a x], {x,0, 0}, Assumptions -> a > 0],
I got an asymptotic expansion of MeijerG function at $x\approx 0$ for $a>0$ as $$G_{1,3}^{3,1}\left(a x\left| \begin{array}{c} -\frac{1}{2} \\ -\frac{1}{2},-\frac{1}{2},\frac{1}{2} \\ \end{array} \right.\right)\xrightarrow{} \frac{-\log (a)-\log (x)-2 \gamma }{\sqrt{a} \sqrt{x}}+O\left(\sqrt{x}\right).$$
However, with available asymptotic expression in MeijerG, I could not derive above expression.
Does anyone have an idea of deriving this?
Here is the plot:

We can roundtrip your expression by first finding the
MellinTransformand then instead of applyingInverseMellinTransform, we can just take the first residue. This ends up giving the main series term. [formula]Here's the Mellin transform:
Now in order to find where to take the residue, we need to know the strip in which this transform exists:
It turns out summing residues in the complex plane to the right of the strip will diverge and so we sum to the left. This means we'll start with
s == 1/2: