The mode of inheritance of total milk yield and its genetic parameters were investigated in Churra dairy sheep through segregation analyses using a Monte Carlo Markov Chains (MCMC) method. Data which consisted of 7,126 lactations belonging to 5,154 ewes were collected between 1999 and 2002 from 15 Spanish Churra dairy flocks. A postulated major gene was assumed to be additive and priors used for variance components were uniform. Based on 50 000 Gibbs samples from ten replicates chains of 100,000 cycles, the estimated marginal posterior means±posterior standard deviations of variance components of milk yield were 23.17±18.42, 65.20±25.05, 120.40±42.12 and 420.83±40.26 for major gene variance ( 2 G σ ), polygenic variance ( 2 u σ ), permanent environmental variance ( 2 pe σ ) and error variance ( 2 e σ ), respectively. The results of this study showed the postulated major locus was not significant, and the 95% highest posterior density regions (HPDs 95% ) of most major gene parameters included 0, and particularly for the major gene variance. The estimated transmission probabilities for the 95% highest posterior density regions (HPDs 95% ) were overlapped. These results indicated that segregation of a major gene was unlikely and that the mode of inheritance of total milk yield in Churra dairy sheep is purely polygenic. Based on 50,000 Gibbs samples from ten replicates chains of 100,000 cycles, the estimated polygenic heritability and repeatability were h 2 = 0.20±0.05 and r = 0.34±0.06, respectively.