Hydrothermal liquefaction (HTL) of algae is a promising biofuel production technology. However, it is always difficult and time-consuming to identify the best optimal conditions of HTL for different algae by the conventional experimental study. Therefore, machine learning (ML) algorithms were applied to predict and optimize bio-oil production with algae compositions and HTL conditions as inputs, and bio-oil yield (Yield_oil), and the contents of oxygen (O_oil) and nitrogen (N_oil) in bio-oil as outputs. Results indicated that gradient boosting regression (GBR, average test R ∼ 0.90) exhibited better performance than random forest (RF) for both single and multi-target tasks prediction. Furthermore, the model-based interpretation suggested that the relative importance of operating conditions (temperature and residence time) was higher than algae characteristics for the three targets. Moreover, ML-based reverse and forward optimizations were implemented with experimental verifications. The verifications were acceptable, showing great potential of ML-aided HTL for producing desirable bio-oil.