Creating subdomains and boundary conditions on imported gmsh mesh

I might be misunderstanding the question, but Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio (that you probably know already) should have everything you need.

The subdomain and boundary markers are those that you set in Physical Surface and Physical Volume.