Genome shuffling methods were explored for Bacillus subtilis strain molecular breeding. Recycling protoplast fusion, recycling transformation and recycling universal transduction were used for genome shuffling in B. subtilis. Four strains with different nutrition-deficiency markers were used as initial strains. After five rounds protoplast fusion, transformation or transduction, the descendant with 4 markers had not been detected, and the rate of descendant with 3 markers were 4.53 x 10(-4), 1.64 x 10(-4), 4.47 x 10(-3), respectively. A computer program was made to simulate the recycling fusion process. Based on simulation result and comparing the genome shuffling result of B. subtilis in this experiment and that of Streptomyces coelicolor reported in references, effective genome shuffling needs a high recombination rate of at least between 10(-3) and 10(-2).