A recent increase in HIV diagnoses among men-having-sex-with-men (MSM) has been shown by surveillance data from Europe and Italy, and new approaches to inferring viral population dynamics from heterochronously sampled gene sequences have been developed. The aim of this study was to reconstruct the epidemiological history of HIV-1 subtype B in a homogeneous group of Italian MSM using a coalescent-based Bayesian framework. A total of 125 HIV-1 subtype B pol sequences were analyzed using Bayesian methods and a relaxed molecular clock to reconstruct their dated phylogeny and estimate population dynamics. At least 10 epidemiological clusters of 3-9 isolates were identified: half including the largest clades originated in the early 1990s and the other half radiated from 1999. Demographic analysis showed that the HIV epidemic grew in accordance with a logistic model characterized by a rapid exponential increase in the effective number of infections (r = 1.54 year) starting from the early 1980s and reaching a plateau 10 years later. Our data suggest that the HIV B epidemic entered our MSM population through multiple transmission chains about 20 years later than in other Western European country. Epidemiological clusters originating in the early 2000s suggest a recent re-emergence of HIV in Italian MSM.