The R2021b release of the System Identification ToolboxTM for MATLAB contains new features that enable the use of machine learning techniques for nonlinear system identification. With this release it is possible to build nonlinear ARX models with regression tree ensemble and Gaussian process regression mapping functions. The release contains several other enhancements including, but not limited to, (a) online state estimation using the extended Kalman filter and the unscented Kalman filter with code generation capability; (b) improved handling of initial conditions for transfer functions and polynomial models; (c) a new architecture of nonlinear black-box models that streamlines regressor handling, reduces memory footprint and improves numerical accuracy; and (d) easy incorporation of identification apps in teaching tools and interactive examples by leveraging the Live Editor tasks of MATLAB. Copyright (C) 2021 The Authors.