I can see the kummer transformation of the confluent hypergeometric function of first kind throught the integral representation. However, I failed to see that for the second kind. More specificially, if $$U(a,b,z)=\frac{1}{\Gamma(a)}\int_0^\infty e^{-z t}t^{a-1}(1+t)^{b-a-1}dt,$$ then why $$U(a,b,z)=z^{1-b}U(1+a-b,2-b,z)?$$
Any references would be appreciated.
Term by term expansion. We reduce the original integral to obviously symmetric Pochhammer symbols.
We need several facts.
${\displaystyle \intop_{t=0}^{\infty}}e^{-z\cdot t}\cdot t^{\left(w-1\right)}=\Gamma\left(w\right)\cdot z^{-w}$
$\alpha\in\mathbb{C}$
$\left(\alpha\right)_{n}=\frac{\Gamma\left(\alpha+n\right)}{\Gamma\left(\alpha\right)}$
$\left(\begin{array}{c} \alpha\\ n \end{array}\right)=\frac{\Gamma\left(\alpha+1\right)}{\Gamma\left(n+1\right)\cdot\Gamma\left(\alpha-n+1\right)} $
Letting $m,n\in Z$ and $s\in\mathbb{C}$ and using http://mathworld.wolfram.com/BinomialCoefficient.html $\frac{\Gamma\left(s-m+1\right)}{\Gamma\left(s-l+1\right)}=\left(-1\right)^{l-m}\frac{\Gamma\left(l-s\right)}{\Gamma\left(m-s\right)}$
Starting from the given integral: $${\displaystyle \intop_{t=0}^{\infty}}\frac{1}{\Gamma\left(a\right)}e^{-z\cdot t}\cdot t^{a-1}\cdot\left(1+t\right)^{\left(b-a-1\right)}dt$$
We select the $n^{th}$term of the binomial to work with implicit summation at the end. This assumes the the binomial summation and integral can be interchanged.
$\left[t^{n}\right]\left(1+t\right)^{\left(b-a-1\right)}=\frac{\Gamma\left(\left(b-a-1\right)+1\right)}{\Gamma\left(n+1\right)\cdot\Gamma\left(\left(b-a-1\right)-n+1\right)}=\frac{\Gamma\left(b-a\right)}{\Gamma\left(n+1\right)\cdot\Gamma\left(\left(b-a\right)-n\right)} $ [Gamma-reciprocal]
and the term integral is: $\frac{1}{\Gamma\left(a\right)}\cdot\frac{\Gamma\left(b-a\right)}{\Gamma\left(n+1\right)\cdot\Gamma\left(\left(b-a\right)-n\right)}\cdot{\displaystyle \intop_{t=0}^{\infty}}t^{\left(a-1\right)+n}\cdot e^{-zt}dt $
$=\frac{1}{\Gamma\left(a\right)}\cdot\frac{\Gamma\left(b-a\right)}{\Gamma\left(n+1\right)\cdot\Gamma\left(\left(b-a\right)-n\right)}\cdot\Gamma\left(a+n\right)\cdot z^{-\left(a+n\right)} $ Parsing $\left(\frac{\Gamma\left(a+n\right)}{\Gamma\left(a\right)}\right)\cdot\left(\frac{\Gamma\left(b-a\right)}{\Gamma\left(\left(b-a\right)-n\right)}\right)\cdot\frac{z^{-n}z^{-a}}{\Gamma\left(n+1\right)} $
Using [Gamma-reciprocal]:Taking $s=b-a,m=1,l=n+1$
$=\left(a\right)_{n}\cdot\left(\frac{\Gamma\left(b-a\right)}{\Gamma\left(\left(b-a\right)-n\right)}\right)\cdot\frac{z^{-n}z^{-a}}{\Gamma\left(n+1\right)}=\left(a\right)_{n}\cdot\left(\frac{\Gamma\left(\left(n+1\right)-\left(b-a\right)\right)}{\Gamma\left(1-\left(b-a\right)\right)}\right)\cdot\left(-1\right)^{n}\cdot\frac{z^{-n}z^{-a}}{\Gamma\left(n+1\right)} $
$=\left(a\right)_{n}\cdot\left(a-b+1\right)_{n}\cdot\left(-1\right)^{n}\cdot\frac{z^{-n}z^{-a}}{\Gamma\left(n+1\right)} $
Which is obviously symetric with: $a'=a-b+1,b'=2-b$
And the power $z^{1-b'}\cdot z^{-a}\rightarrow z^{1-b'}z^{-\left(a'-b'+1\right)}\rightarrow z^{-a'}$