Thermostability of beta-glucosidase was enhanced by family shuffling, site saturation mutagenesis, and site-directed mutagenesis. Family shuffling was carried out based on beta-glucosidase BglC from Thermobifida fusca and beta-glucosidase BglB from Paebibacillus polymxyxa with the help of synthetic primers. High-throughput screening revealed mutants with higher thermostability than both parental enzymes. The most: thermostable mutant VM2 containing three key amino acid changes in L444Y/G447S/A433V had a 144-fold increase in half-life of inactivation as compared to the parental enzyme BglC. The mutant VM2 showed 28% and 94% increase in k(car) towards p-nitrophenyl-beta-D-glucopyranoside (pNPG) and cellobiose, respectively. The mutant with enhanced stability would facilitate the recycle of beta-glucosidase in the bioconversion of cellulosic biomass. (C) 2010 Elsevier Ltd. All rights reserved.