In this paper, we find integral presentations of incomplete gamma functions. Using these results we construct bivariate gamma distributions and prove some properties of such distributions.