python - How to do a Spatial Join after opening a GDB file with Fiona?

After I successfully open up an ESRI GDB file/folder with Fiona, what open-source Python package should I use to do a spatial join?

I was told Shapely is not the appropriate package for this.

I am completely new to GIS.