Since the beginning of the pandemic, SARS-CoV-2 has shown a great genomic variability, resulting in the continuous emergence of new variants that has made their global monitoring and study a priority. This work aimed to study the genomic heterogeneity, the temporal origin, the rate of viral evolution and the population dynamics of the main circulating variants (20E.EU1, Alpha and Delta) in Italy, in August 2020-January 2022 period. For phylogenetic analyses, two datasets were evaluated for each variant, the former comprising international genomes and the latter focusing on clusters containing at least 70% of Italian sequences. The international dataset showed 26 (23% Italians, 23% singleton, 54% mixed), 40 (60% mixed, 37.5% Italians, 1 singleton) and 42 (85.7% mixed, 9.5% singleton, 4.8% Italians) clusters with at least one Italian sequence, in 20E.EU1 clade, Alpha and Delta variants, respectively. International clusters presented tMRCAs between 13/06/2020-27/07/2021. R e values showed the highest level between May-June until autumn 2020 in 20E.EU1 clade. The Alpha variant showed an increase in the R e in December 2020, when the highest mean value was estimated. Delta variant presented two peaks: the first between March-May, and the second between June-July 2021. The present work highlighted a different evolutionary dynamic of studied lineages with a high concordance between epidemiological parameters estimation and phylodynamic trends suggesting that the mechanism of replacement of the SARS-CoV-2 variants must be related to a complex of factors involving the transmissibility, as well as the implementation of control measures, and the level of cross-immunization within the population.
