The superfluid transition of a repulsive Bose gas in the presence of a sinusoidal potential which represents a simple-cubic optical lattice is investigated using quantum Monte Carlo simulations. At the average filling of one particle per well the critical temperature has a nonmonotonic dependence on the interaction strength, with an initial sharp increase and a rapid suppression at strong interactions in the vicinity of the Mott transition. In an optical lattice the positive shift of the transition is strongly enhanced compared to the homogenous gas. By varying the lattice filling we find a crossover from a regime where the optical lattice has the dominant effect to a regime where interactions dominate and the presence of the lattice potential becomes almost irrelevant. The combined effect of interparticle interactions and external potentials plays a fundamental role in determining the quantum-coherence properties of several many-body systems, including He in Vycor or on substrates, paired electrons in superconductors and in Josephson junction arrays, neutrons in the crust of neutron stars [1], and ultracold atoms in optical potentials. However, even the (apparently) simple problem of calculating the superfluid transition temperature T c of a dilute homogeneous Bose gas has challenged theoreticians for decades [2]. Many techniques have been employed, obtaining contradicting results, differing even in the functional dependence of T c on the interaction parameter (the two-body scattering length a) and in the sign of the shift with respect to the ideal gas transition temperature T 0 c (for a review see Ref.[3]). In the weakly interacting limit, the shift of the critical temperature 4,5], where n is the density and the coefficient c ¼ 1.29ð5Þ was determined using Monte Carlo simulations of a classical-field model defined on a discrete lattice [6,7]. Continuous-space quantum Monte Carlo simulations of Bose gases with short-range repulsive interactions have shown that this linear form is valid only in the regime n 1=3 a ≲ 0.01, while at stronger interaction T c reaches a maximum where ΔT c =T 0 c ≃ 6.5% and then decreases for n 1=3 a ≳ 0.2 [8]. This suppression of T c occurs in a regime where universality in terms of the scattering length is lost and other details of the interaction potential become relevant [8][9][10]. In recent years ultracold atomic gases have emerged as the ideal experimental test bed for many-body theories [11]. However, the direct measurement of interactions effects on T c has been hindered by the presence of the harmonic trap. In the presence of confinement the main interactions effect can be predicted by mean-field theory and is due to the broadening of the density profile [12], leading to a suppression of T c . Deviations from the mean-field prediction and effects due to critical correlations have been measured in Refs. [13,14]. A major breakthrough has been achieved recently with the realization of Bose-Einstein condensation in quasiuniform trapping potentials [15]. This setup allows a more ...