How do I evaluate the following combination of random variables? Is it martingale? I'm about to analyse the following expression
$$Z_n:=\prod_{k=1}^n \left(\frac{\frac{Y_k}{\prod_{i=1}^k X_i}}{\sum_{j=1}^k \frac{Y_j}{\prod_{i=1}^j X_i}} \right),$$
where $Y_j$ for all $j\in \mathbb{N}$, are independent and $\Gamma(\beta,1)$-distributed random variables and $X_i$, for all $i \in \mathbb{N}$ independent identically BetaDistr.($\alpha,\beta$). I have presumtion that this expression is a martingale and tried to prove the martingale property. Let $\mathcal{F_n}=\{X_i| Y_i: i \le n\}$
$$\mathbb{E}[Z_{n+1}|\mathcal{F}_n]=Z_n\mathbb{E}\left[\frac{\frac{Y_{n+1}}{\prod_{i=1}^{n+1} X_i}}{\sum_{j=1}^{n+1} \frac{Y_j}{\prod_{i=1}^j X_i}} |\mathcal{F}_n \right]=Z_n \frac{1}{\prod_{i=1}^n X_i}\mathbb{E}\left[\frac{\frac{Y_{n+1}}{X_{n+1}}}{\sum_{j=1}^{n} \frac{Y_j}{\prod_{i=1}^j X_i}+\frac{Y_{n+1}}{\prod_{i=1}^{n+1} X_i} } |\mathcal{F}_n \right]=?$$
I could calculate only till this step. Does somebody see how to proceed further? Maybe somebody has already had experience with this combination of random variables and could give me a hint what are another options to evaluate this expression?

With the hint from @Davide Giraudo I started to investigate the expectation $\mathbb{E}\left[\frac{\frac{Y}{X}}{a+b\frac{Y}{X}}\right]$, where $Y$ has $\Gamma(\beta,1)$-distribution and $X$ BetaDistr.($\alpha,\beta$).
Let define $Z=\frac{Y}{X}$. 
In the paper of S. Nadarajah and S. Kotz "On the product and Ration of Gamma and Beta Random Variables" I could find the pdf of $Z=\frac{Y}{X}$ and it's given in paper as
$$f(z)=\frac{\lambda^{\beta}B(\beta+a,b)}{\Gamma(\beta)B(a,b)}z^{\beta-1}F(\beta+a;\beta+a+b;-\lambda z) \text{ for } z>0,$$
where $F(a;b;x)=\sum_{k=0}^{\infty} \frac{(a)_k}{(b)_k}\frac{x^k}{k!}$ is the confluent hypergeometric function, $\Gamma(\beta)$ express Gamma function and $Beta(\alpha,\beta)$ Beta function.
Adjusted to our case: $\beta=\beta$, $\lambda=1$, $a=\alpha, b=\beta$ it results to
$$f(z)=\frac{B(\beta+\alpha,\beta)}{\Gamma(\beta)B(\alpha,\beta)}z^{\beta-1}F(\beta+\alpha;\beta+\alpha+\beta;-z), \text{ for } z>0.$$
Like I understand I have to calculate the following integral
$\mathbb{E}\left[\frac{Z}{a+bZ}\right]=\int_{-\infty}^{\infty} \left( \frac{z}{a+bz}\right)\frac{B(\beta+\alpha,\beta)}{\Gamma(\beta)B(\alpha,\beta)}z^{\beta-1}F(\beta+\alpha;\beta+\alpha+\beta;-z)dz$ 
I'm working on that. I hope I'm on the right way.
