We use the semi-analytic method of \citet{Esma91} to calculate the abundances of Helium and Deuterium produced during Big Bang nucleosynthesis assuming the fine structure constant and the Higgs vacuum expectation value may vary in time. We analyze the dependence on the fundamental constants of the nucleon mass, nuclear binding energies and cross sections involved in the calculation of the abundances. Unlike previous works, we do not assume the chiral limit of QCD. Rather, we take into account the quark masses and consider the one-pion exchange potential, within perturbation theory, for the proton-neutron scattering. However, we do not consider the time variation of the strong interactions scale but attribute the changes in the quark masses to the temporal variation of the Higgs vacuum expectation value. Using the observational data of the helium and deuterium, we put constraints on the variation of the fundamental constants between the time of nucleosynthesis and the present time.

Variation of constants in the very early universe can generate inflation. We consider a scenario where the strong coupling constant was changing in time and where the gluon condensate underwent a phase transition ending the inflation.

We consider a scenario where the strong coupling constant was changing in the early universe. We attribute this change to a variation in the colour charge within a Bekenstein-like model. Treating the vacuum gluon condensate $<G^2>$ as a free paramete r, we could generate inflation with the required properties to solve the fluctuation and other standard cosmology problems. A possible approach to end the inflation is suggested.

We propose to generalize Bekenstein model for the time variation of the fine structure "constant" $\alpha_{em}$ to QCD strong coupling constant $\alpha_S$. We find that, except for a ``fine tuned'' choice of the free parameters, the extension can not be performed trivially without being in conflict with experimental constraints and this rules out $\alpha_S$ variability. This is due largely to the huge numerical value of the QCD vacuum gluon condensate when compared to the matter density of the universe.