The presence of dyes in the aquatic environment is a global problem because of the negative consequences on the ecosystem quality. Adsorption process is used as a low cost and effective method for removing hazardous dyes. This study investigated removal of Bromophenol Blue from aqueous solutions by Fomes fomentarius and Trametes versicolor. The statistical relationships of the parameters in the biosorption process by considering three controllable factors including initial pH value (2, 7, 12), initial dye concentration (25, 50, 75 mg L-1 ) and contact time (30, 150, 270 minutes) were investigated at three different levels and the biosorption process was optimized using Taguchi orthogonal array (L9 OA) experimental design. As a result of experimental studies, the optimal dye removal conditions were determined as the contact time of 270 minutes, pH: 2 and initial Bromophenol Blue concentration of 50 mg L-1 for F. fomentarius. The optimal dye removal conditions were determined as the contact time of 150 minutes, pH: 2 and initial Bromophenol Blue concentration of 50 mg L-1 for T. versicolor. For F. fomentarius and T. versicolor, the removal efficiencies at optimum conditions were 78.34% and 69.28%, respectively. Biosorption isotherms were modeled using the Langmuir and Freundlich models. The biosorption of Bromephenol Blue on F. fomentarius and T. versicolor fitted better in the Freundlich and Langmuir models by non-linearized equations, respectively. Both biosorbents can be promising for dyes removal from aqueous environment. [ABSTRACT FROM AUTHOR]