The anisotropy of attractive interactions between particles can favor, through a self-assembly process, the formation of linear semi-flexible chains. In the appropriate temperatures and concentration ranges, the growing aspect ratio of the aggregates can induce formation of a nematic phase, as recently experimentally observed in several biologically relevant systems. We present here a numerical study of the isotropic-nematic phase boundary for a model of bifunctional polymerizing hard cylinde...