Three thermostable lactose-hydrolases, namely, two β-glycosidases (bglA and bglB) and one β-galactosidase (bgaA) genes were cloned from the genomic library of Thermus sp. IB-21. The bglA, bglB, and bgaA consisted of 1311bp (436 amino acid residues), 1296bp (431aa), and 1938bp (645 aa) of nucleotides with predicted molecular masses of 49,066, 48,679, and 72,714Da, respectively. These enzymes were overexpressed in Escherichia coli BL21(DE3) using pET21b(+) vector system. The recombinant enzymes were purified to homogeneity by a heat precipitation (70°C, 40min) and a Ni²⁺-affinity chromatography. The molecular masses of the purified enzymes estimated by SDS–PAGE agreed with their predicted values. All the purified enzymes showed their optimal pH at around 5.0–6.0. In contrast, the temperature profiles for activity and thermostability patterns were different for each enzyme. BglB β-glycosidase displayed the best lactose hydrolysis activity of the three enzymes without substrate inhibition up to 200mM lactose at 70°C and pH 7.0. The specific activities (U/mg) of BglA, BglB, and BgaA on 138mM lactose at 70°C and pH 7.0 were 36.8, 160.3, and 8.5, respectively.