Does sum of random variables follows gaussian distribution although each have different mean and variances?
You can add random variables with different means and variances. For example, look at two continuous uncorrelated normal-distributed random variables $x_1$ and $x_2$ with means and variances $\mu_1, \sigma_1^2$ and $\mu_2, \sigma_2^3$
Then $x_3 = x_1 + x_2 $ is also a normal-distributed random variable with $\mu_3 = \mu_1 + \mu_2$ and $ \sigma_3^2 = \sigma_1^2 + \sigma_2^2 $. For uncorrelated random variables the means and variances just sum.
The distribution of a sum of two continuous random variables is the convolution of the individual  distributions. That's actually the core of the central limit theorem. If you convolve a Gaussian with a Gaussian you get another Gaussian (see for example https://jeremy9959.net/Math-5800-Spring-2020/notebooks/convolution_of_gaussians.html) . If you keep convolving reasonably looking functions with each other eventually it will approach a Gaussian.
A simple example is  summing a few uniformly distributed random variables. The first one has a rectangular distribution. The second one is triangular. The third one already shows a nice round belly in the middle and once you are up to summing 10 or 12 it looks very much Gaussian. That's in fact one method how gaussian random number generators are implemented.
For the Central Limit Theorem, there is no requirement for random variables to have the same mean or variance.