Aligning BOLD data across sessions in functional output space

Is there a way in fmriprep to align BOLD data from multiple runs and across sessions with the output being in func space (i.e. not registering to a standard or anatomical space?
Ideally, this process would involve multiple steps: 1. calculating the motion transforms within run using mcflirt and 2. use flirt to calculate the transform between each run’s boldref and a single template (let’s say from the very first run) and 3. reslicing in one step. Is anything like this possible in fmriprep?

Hi some attemps have been made, but nothing is ready out of the box in fmriprep right now: