An approach to adjust for covariates affecting features and/or target in TPOT, based on regressing out the covariates in a manner that avoids 'leakage' during the cross-validation training procedure. Applications of this approach have been demonstrated for toxicogenomics and schizophrenia gene expression data sets.