And using the run_harmony function (in pyharmony package, version 1.0.7) and combat function (in Scanpy) methods to deal with batch effect issues if batch effect existing in dataset.