Thermostability of β-glucosidase was enhanced by family shuffling, site saturation mutagenesis , and site-directed mutagenesis . Family shuffling was carried out based on β-glucosidase BglC from Thermobifida fusca and β-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 cat towards p -nitrophenyl-β- d -glucopyranoside (pNPG) and cellobiose , respectively. The mutant with enhanced stability would facilitate the recycle of β-glucosidase in the bioconversion of cellulosic biomass.