We study the production of primordial black hole (PBH) binaries and the resulting merger rate, accounting for an extended PBH mass function and the possibility of a clustered spatial distribution. Under the hypothesis that the gravitational wave events observed by LIGO were caused by PBH mergers, we show that it is possible to satisfy all present constraints on the PBH abundance, and find the viable parameter range for the lognormal PBH mass function. The non-observation of a gravitational wave background allows us to derive constraints on the fraction of dark matter in PBHs, which are stronger than any other current constraint in the PBH mass range 0.5 − 16M . We show that the predicted gravitational wave background can be observed by the coming runs of LIGO, and its non-observation would indicate that the observed events are not of primordial origin. As the PBH mergers convert matter into radiation, they may have interesting cosmological implications, for example in the context of relieving the tension between high and low redshift measurements of the Hubble constant. However, we find that these effects are suppressed as, after recombination, no more that 1% of dark matter can be converted into gravitational waves.