Role of Direct Sexual Contact in Human Transmission of Monkeypox Virus, Italy.

The 2022 global mpox outbreak was driven by human-to-human transmission, but modes of transmission by sexual relationship versus sexual contact remain unclear. We evaluated sexual transmission of mpox by using monkeypox virus (MPXV) G2R-mRNA as a marker of ongoing viral replication through in vitro experiments. We analyzed clinical samples of 15 MPXV-positive patients in Italy from different biological regions by using the setup method. The presence of MPXV DNA, MPXV G2R-mRNA, or both in all analyzed lesion swab samples, independent of viral load, confirmed a higher infectivity risk from skin lesions. Positivity for MPXV G2R-mRNA in nasopharyngeal swabs was associated with high MPXV load, whereas positive results for MPXV G2R-mRNA were obtained only in the 2 semen samples with the lowest MPXV loads. Our results suggest that close or skin-to-skin contact during sexual intercourse is the main route of sexual transmission and that semen is a minor driver of infection, regardless of MPXV load.

HEV-3 subtypes and strains detected in cases of HEV infection in central Italy from 2015 to 2023.

PURPOSE: HEV is an emerging pathogen in Europe and was previously shown to be hyperendemic in areas of Abruzzo and Lazio, Central Italy. No systematic analysis of the HEV strains responsible for human infections over several years in Central Italy has previously been reported. Aim of the study was the molecular characterization of HEV from autochthonous hepatitis E cases occurred in Abruzzo and Lazio between 2015 and 2023. METHODS: Samples from 118 cases collected as part of virological surveillance in Abruzzo and Lazio from 2015 to 2023 were subjected to HEV sequencing and phylogenetic analysis. RESULTS: The main observed subtype was 3f, followed by 3c and 3e. The annual subtype distribution was quite stable over the observation period, but 3f cases tended to concentrate in winter/early spring whereas 3e cases in summer. Phylogenetic clusters of highly related sequences (a) highlighted unrecognized "point source outbreaks", (b) provided molecular support to temporally and/or geographically linked cases and (c) provided evidence for transmission of identical/highly related strains up to months/years following their first detection. CONCLUSIONS: The data provide an overview of the HEV strains responsible for human infections over eight years in Central Italy. The observed subtype distribution appears to agree better with the subtype distribution reported in Italy in pigs rather than in geographically matched wild boars, suggesting pig and its derivate food was a more frequent source of infection than wild boar in Abruzzo and Lazio. Molecular characterization is essential to recognize "point source outbreaks" and to monitor HEV circulation.

Genomic surveillance of SARS-CoV-2 positive passengers on flights from China to Italy, December 2022.

With numbers of COVID-19 cases having substantially increased at the end of 2022 in China, some countries have started or expanded testing and genomic surveillance of travellers. We report screening results in Italy in late December 2022 of 556 flight passengers in provenance from two Chinese provinces. Among these passengers, 126 (22.7%) tested SARS-CoV-2 positive. Whole genome sequencing of 61 passengers' positive samples revealed Omicron variants, notably sub-lineages BA.5.2.48, BF.7.14 and BQ.1.1, in line with data released from China.

Analysis of hepatitis B virus-mixed genotype infection by ultra deep pyrosequencing in Sudanese patients, 2015-2016.

PURPOSE: The frequency of detection of HBV co-infection with multiple HBV genotypes is influenced by the detection method; usually co-infections are detected by multiplex PCR or hybridization assays, and are rarely confirmed by sequencing and conventional cloning. The objective of this study was to confirm by ultra-deep pyrosequencing (UDPS) mixed HBV infections, previously detected by multiplex-nested PCR. METHODS: Sixteen samples from HBV co-infected Sudanese patients detected by multiplex-nested PCR, were amplified targeting the P/S region and sequenced by UDPS. RESULTS: The only genotypes identified using UDPS were D and E, while A, B, C and F genotypes, previously detected by multiplex-nested PCR, were not detected. Specifically, 10 samples were shown to be mono-infected (D or E); in 3 out of 10 mono-infected D patients, a subtype combination was observed: D1 + D7 in 2 cases and D2 + D6 in 1 case. The remaining 6 subjects were D + E co-infected (harboring different mixtures of D subtypes). CONCLUSIONS: Overall, UDPS is more effective than multiplex-nested PCR for identifying multiple HBV genotypes and subtypes infections.

Ultra-Deep Sequencing Characterization of HCV Samples with Equivocal Typing Results Determined with a Commercial Assay.

Hepatitis C virus (HCV) is classified into seven phylogenetically distinct genotypes, which are further subdivided into related subtypes. Accurate assignment of genotype/subtype is mandatory in the era of directly acting antivirals. Several molecular methods are available for HCV genotyping; however, a relevant number of samples with indeterminate, mixed, or unspecified subtype results, or even with misclassified genotypes, may occur. Using NS5B direct (DS) and ultra-deep pyrosequencing (UDPS), we have tested 43 samples, which resulted in genotype 1 unsubtyped (n = 17), mixed infection (n = 17), or indeterminate (n = 9) with the Abbott RealTime HCV Genotype II assay. Genotype 1 was confirmed in 14/17 samples (82%): eight resulted in subtype 1b, and five resulted in subtype 1a with both DS and UDPS, while one was classified as subtype 1e by DS and mixed infection (1e + 1a) by UDPS. Three of seventeen genotype 1 samples resulted in genotype 3h with both sequencing approaches. Only one mixed infection was confirmed by UDPS (4d + 1a), while in 88% of cases a single component of the mixture was detected (five genotype 1a, four genotype 1b, two genotype 3a, two genotype 4m, and two genotype 4d); 44% of indeterminate samples resulted genotype 2c by both DS and UDPS, 22% resulted genotype 3a; one indeterminate sample by Abbott resulted in genotype 4d, one resulted in genotype 6n, and one was classified as subtype 3a by DS, and resulted mixed infection (3a + 3h) by UDPS. The concordance between DS and UDPS was 94%, 88%, and 89% for genotype 1, co-infection, and indeterminate results, respectively. UDPS should be considered very useful to resolve ambiguous HCV genotyping results.

Torquetenovirus Viremia Quantification Using Real-Time PCR Developed on a Fully Automated, Random-Access Platform.

Quantification of Torquetenovirus (TTV) viremia is becoming important for evaluating the status of the immune system in solid organ transplant recipients, monitoring the appearance of post-transplant complications, and controlling the efficacy of maintenance immunosuppressive therapy. Thus, diagnostic approaches able to scale up TTV quantification are needed. Here, we report on the development and validation of a real-time PCR assay for TTV quantification on the Hologic Panther Fusion® System by utilizing its open-access channel. The manual real-time PCR previously developed in our laboratories was optimized to detect TTV DNA on the Hologic Panther Fusion® System. The assay was validated using clinical samples. The automated TTV assay has a limit of detection of 1.6 log copies per ml of serum. Using 112 samples previously tested via manual real-time PCR, the concordance in TTV detection was 93% between the assays. When the TTV levels were compared, the overall agreement between the methods, as assessed using Passing-Bablok linear regression and Bland-Altman analyses, was excellent. In summary, we validated a highly sensitive and accurate method for the diagnostic use of TTV quantification on a fully automated Hologic Panther Fusion® System. This will greatly improve the turnaround time for TTV testing and better support the laboratory diagnosis of this new viral biomarker.

MPXV DNA kinetics in bloodstream and other body fluids samples.

Since spring 2022, the global epidemiology of the monkeypox virus (MPXV) has changed. The unprecedented increase of human clade II MPXV cases worldwide heightened concerns about this emerging zoonotic disease. We analysed the positivity rates, viral loads, infectiousness, and persistence of MPXV DNA for up to 4 months in several biological samples from 89 MPXV-confirmed cases. Our data showed that viral loads and positivity rates were higher during the first two weeks of symptoms for all sample types. Amongst no-skin-samples, respiratory specimens showed higher MPXV DNA levels and median time until viral clearance, suggesting their usefulness in supporting MPXV diagnosis, investigating asymptomatic patients, and monitoring viral shedding. Infectious virus was cultured from respiratory samples, semen, and stools, with high viral loads and collected within the first 10 days. Notably, only one saliva and one semen were found positive for viral DNA after 71 and 31 days from symptoms, respectively. The focus on bloodstream samples showed the best testing sensitivity in plasma, reporting the overall highest MPXV DNA detection rate and viral loads during the 3-week follow-up as compared to serum and whole-blood. The data here presented can be useful for MPXV diagnostics and a better understanding of the potential alternative routes of its onward transmission.

Phylogenetic and Phylodynamic Analysis of Delta Strains Circulating in Italy.

The hepatitis delta virus (HDV) exhibits high genetic and evolutionary variability and is classified into eight genotypes (HDV-1 to -8). HDV-1 is the most widespread genotype worldwide and includes several subtypes. It predominates mainly in Europe, the Middle East, North America, and Northern Africa, and is associated with both severe and mild forms of liver disease. In this study, we performed phylogenetic and phylodynamic analyses of HDV strains circulating in Regione Lazio, Italy, to understand when these strains were introduced into the Lazio region and to define their genetic variability in Italy. Fifty HDV RNA positive patient samples were amplified using a nested RT-PCR approach targeting the HDV R0 region and sequenced. A phylogenetic tree of patient-derived sequences and reference sequences representing HDV-1 to -8 was constructed using the GTRGAMMA model in RAxML v8. The results indicated that HDV-1 was the predominant genotype with HDV-1d being the most frequently inferred subtype. HDV-1 sequences clustering with subtypes 1b and 1e were also identified. A phylodynamic analysis of HDV-1 sequences employing a Bayesian birth-death model inferred a clock rate of 3.04 × 10-4 substitutions per site per million years, with a 95% Highest Posterior Density (HPD) interval of 3.45 × 10-5 to 5.72 × 10-4. A Bayesian birth-death analysis with tree calibration based on a sample dating approach indicated multiple original sources of infection (from the late 1950s to late 1980s). Overall, these results suggest that HDV sequences from the native Italian and non-Italian patients analyzed in this study represent multiple lineages introduced across a wide period. A common ancestral origin should be excluded.

Development and validation of a nanoplate-based digital PCR assay for absolute MPXV quantification.

Quantification of mpox virus (MPXV) across different human body anatomical sites can provide insights about the most likely transmission routes, so methods able to release absolute and exact quantitative values of MPXV DNA are crucial. Here, we optimized a new QIAcuity digital PCR (dPCR) protocol for the detection and quantification of MPXV DNA in clinical samples and assessed the performance of the assay by comparing the results obtained in 144 biological samples with those resulting from the use of an in-house real-time PCR (qPCR). Overall, the concordance between the two assays was 95%, with samples identified concordantly as MPXV DNA positive and having a mean number of copies per µl of 1708 (95% CI: 107-2830 copies/µl). The remaining samples gave discordant results, with 5 out of 7 detected with the QIAcuity dPCR assay but not with the in-house qPCR. MPXV DNA levels measured by QIAcuity dPCR were strongly correlated with the Ct values detected by in-house qPCR and with those detected by another dPCR assay previously developed in our laboratories. The QIAcuity dPCR assay may be a robust and easy-to-perform method for MPXV DNA quantification in several biological samples.

TTV and other anelloviruses: The astonishingly wide spread of a viral infection.

The broad family of viruses known as anelloviruses (AV) infects both humans and numerous animal species. They have a tiny, covalently closed single-stranded DNA genome and the astonishing capacity to infect a very high percentage of healthy and ill people with chronic infections that could last a lifetime. AV, and particularly the prototype Torquetenovirus, have established a successful interaction with the host's immune system and the rate at which they replicate is a gauge to measure overall immune function, even though many aspects of their life cycle and pathogenesis are still poorly understood.

Torquetenovirus Loads in Peripheral Blood Predict Both the Humoral and Cell-Mediated Responses to SARS-CoV-2 Elicited by the mRNA Vaccine in Liver Transplant Recipients.

Three years into the COVID-19 pandemic, mass vaccination campaigns have largely controlled the disease burden but have not prevented virus circulation. Unfortunately, many immunocompromised patients have failed to mount protective immune responses after repeated vaccinations, and liver transplant recipients are no exception. Across different solid organ transplant populations, the plasma levels of Torquetenovirus (TTV), an orphan and ubiquitous human virus under control of the immune system, have been shown to predict the antibody response after COVID-19 vaccinations. We show here a single-institution experience with TTV viremia in 134 liver transplant recipients at their first or third dose. We found that TTV viremia before the first and third vaccine doses predicts serum anti-SARS-CoV-2 Spike receptor-binding domain (RBD) IgG levels measured 2-4 weeks after the second or third dose. Pre-vaccine TTV loads were also associated with peripheral blood anti-SARS-CoV-2 cell-mediated immunity but not with serum SARS-CoV-2 neutralizing antibody titers.

Detection of SARS-CoV-2 Variants via Different Diagnostics Assays Based on Single-Nucleotide Polymorphism Analysis.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is characterized by fast evolution with the appearance of several variants. Next-Generation Sequencing (NGS) technology is considered the gold standard for monitoring known and new SARS-CoV-2 variants. However, the complexity of this technology renders this approach impracticable in laboratories located in areas with limited resources. We analyzed the capability of the ThermoFisher TaqPath COVID-19 RT-PCR (TaqPath) and the Seegene Novaplex SARS-CoV-2 Variant assay (Novaplex) to detect Omicron variants; the Allplex VariantII (Allplex) was also evaluated for Delta variants. Sanger sequencing (SaS) was the reference method. The results obtained with n = 355 nasopharyngeal samples were: negative with TaqPath, although positive with other qualitative molecular assays (n = 35); undetermined (n = 40) with both the assays; negative for the ∆69/70 mutation and confirmed as the Delta variant via SaS (n = 100); positive for ∆69/70 and confirmed as Omicron BA.1 via SaS (n = 80); negative for ∆69/70 and typed as Omicron BA.2 via SaS (n = 80). Novaplex typed 27.5% of samples as undetermined with TaqPath, 11.4% of samples as negative with TaqPath, and confirmed 100% of samples were Omicron subtypes. In total, 99/100 samples were confirmed as the Delta variant with Allplex with a positive per cent agreement (PPA) of 98% compared to SaS. As undermined samples with Novaplex showed RdRp median Ct values (Ct = 35.4) statistically higher than those of typed samples (median Ct value = 22.0; p < 0.0001, Mann-Whitney test), the inability to establish SARS-CoV-2 variants was probably linked to the low viral load. No amplification was obtained with SaS among all 35 negative TaqPath samples. Overall, 20% of samples which were typed as negative or undetermined with TaqPath, and among them, twelve were not typed even by SaS, but they were instead correctly identified with Novaplex. Although full-genome sequencing remains the elected method to characterize new strains, our data show the high ability of a SNP-based assay to identify VOCs, also resolving samples typed as undetermined with TaqPath.

Kinetics of viral DNA in body fluids and antibody response in patients with acute Monkeypox virus infection.

We report the follow-up laboratory investigation of three MPXV cases infected in May-June 2022 from diagnosis to disease resolution, monitoring viral shedding in different body fluids and antibody kinetics. Out of 138 non-lesion samples, viral DNA was found in 92.3% saliva, 85.7% semen, 86.2% oropharyngeal swabs, 51.7% plasma, 46.1% stool, and 9.5% urine samples. Viral load quantified by digital PCR widely varied, but tend to be higher in oropharyngeal swabs, saliva, and stool. Replication competent virus was recovered from four out of seventeen samples, including 1 saliva, 1 oropharyngeal swabs, 1 semen, and 1 stool. The analysis of the antibody kinetics revealed that IgM, IgA, and IgG antibodies were detected within two weeks post-symptoms onset for all three patients, with IgG detected early on at day 4-8 and IgM and IgA showing lower titers along the time frame of the study. Antibody levels increased during the second week of illness with IgG reaching high titers.

mRNA- and Adenovirus-Based Vaccines against SARS-CoV-2 in HIV-Positive People.

About two years have passed since the identification of SARS-CoV-2 in China. The rapid spread of this virus all over the world and its high transmissibility and pathogenicity in humans have resulted in a global pandemic. The negative impact of COVID-19 on health, society and the economy at the global level has pushed researchers and pharmaceutical companies to develop effective vaccines to fight SARS-CoV-2. Thanks to this collaborative effort, the first COVID-19 vaccine was developed in less than a year. Since then, several COVID-19 vaccines have been validated for use by the World Health Organization. Among these, mRNA- (BNT162b2 and mRNA1273) and adenovirus-based (ChAdOx1) vaccines were developed through the use of novel technologies. While all three of these vaccines have shown effectiveness against the COVID-19 disease and their immunogenicity was characterized in clinical trials in the general population, data on their efficacy and immunogenicity in people living with HIV (PLWH) are limited. In this review, we provide a description of the characteristics of mRNA- and adenovirus-based vaccines and of the immune response elicited in the general population by vaccination. Then we describe the use of these vaccines and their efficacy and immunogenicity in people living with HIV and we conclude with a discussion regarding some open questions concerning the use of mRNA- and adenovirus-based COVID-19 vaccines in PLWH.

Phylogenetic and Phylodynamic Analyses of HCV Strains Circulating among Patients Using Injectable Drugs in Central Italy.

Approximately 71 million people worldwide are infected with the hepatitis C virus (HCV). Injectable drug use represents the most common route of transmission in Europe and other developed countries. We studied the molecular characteristics of the HCV infection among mono-infected people who used drugs (PWUD) in Italy. Among 208 PWUD with anti-HCV antibodies, 101 (48.6%) were HCV RNA-positive, the majority (47%) were infected with the HCV genotype (Gt)1a, followed by Gt3a (34.9%), Gt4 (9.1%), Gt1b (4.5%), and Gt2 (4.5%). Bayesian phylogenetic analyses of clustered HCV NS5B sequences from 66 HCV-positive PWUDs with available plasma samples indicated age and neighborhood proximity as the most common characteristics between closely related HCV strains. Population dynamics, as measured by a coalescent Bayesian skyline analysis, revealed an increase in HCV Gt1a infections from the mid-1980s to mid-1990s. While HCV Gt3a infections were first detected in the 1980s, patient numbers with this genotype subtype remained relatively constant. For both Gt1a and Gt3a, Birth-Death Bayesian Skyline analyses produced higher reproduction numbers post 2014. For earlier time intervals, slow growths were observed for both Gt1a and Gt3a with reproduction numbers (Re) of approximately 1. The evolutionary rates for Gt1a and Gt3a were estimated as 2.23 × 10-4 and 3.85 × 10-4, respectively.

Late Relapse and Reinfection in HCV Patients Treated with Direct-Acting Antiviral (DAA) Drugs.

The risk of hepatitis C virus (HCV) recurrence after direct-acting antiviral (DAA) treatment is <0.5%. However, the distinction between HCV RNA late relapse and reinfection still represents a challenge in virological diagnostics. The aim of this study was to employ next-generation sequencing (NGS) to investigate HCV RNA recurrence in patients achieving a sustained virologic response (SVR) at least six months post-treatment. NGS was performed on plasma samples from six HCV-positive patients (Pt1-6) treated with DAA. NGS of HCV NS5B was analyzed before treatment (T0), after HCV RNA rebound (T1), and, for Pt3, after a second rebound (T2). Reinfection was confirmed for Pt5, and for the first rebound observed in Pt3. Conversely, viral relapse was observed when comparing T0 and T1 for Pt6 and T1 and T2 for Pt3. Z-scores were calculated and used to predict whether HCV-positive patient samples at different time points belonged to the same quasispecies population. A low Z-score of <2.58 confirmed that viral quasispecies detected at T0 and T1 were closely related for both Pt1 and Pt2, while the Z-score for Pt4 was suggestive of possible reinfection. NGS data analyses indicate that the Z-score may be a useful parameter for distinguishing late relapse from reinfection.

Clinical Relevance of Torque Teno Virus (TTV) in HIV/HCV Coinfected and HCV Monoinfected Patients Treated with Direct-Acting Antiviral Therapy.

Torque Teno virus (TTV) is a ubiquitous virus that causes chronic infection in humans with unknown clinical consequences. Here, we investigated the influence of TTV infection on HCV direct-acting antiviral (DAA) efficacy in HIV/HCV coinfected and HCV monoinfected patients as controls. Of 92 study patients, 79.3% were TTV DNA positive; untreated patients exhibited a significantly higher proportion of TTV DNA-positivity vs. sustained virological response (SVR) patients (100.0% vs. 65.2%, p < 0.001), while TTV positivity was not significant in DAA failure patients vs. SVR patients despite HIV/HCV coinfection. TTV DNA viral load was higher among HCV monoinfected patients vs. HIV/HCV coinfected, although marginally significant (p = 0.074) and no significant viral load difference was detected between DAA failures and SVR patients, while untreated vs. SVR patients had a significantly higher viral load (19,884, IQR 5977-333,534, vs. 469, IQR 10-4124, p = 0.004). Alpha-genogroup 3 TTV was the most prevalent genetic group, and no specific strain or genogroup was observed in relapser patients. Among HIV/HCV patients with HCV RNA detectable at end of treatment (EOT), TTV DNA was detected in 9/17 treatment responder patients and 3/5 relapser patients, thus, TTV infection does not appear to influence the control HCV viremia after EOT. Levels of IL-6 IL-4, and CD14 were not significantly different between TTV PCR-positive and -negative patients. These results suggest no association between TTV DNA positivity or viral load and HCV DAA failure whether patients were HIV/HCV coinfected or HCV monoinfected.

Anti-HEV IgG Avidity Testing: Utility for Diagnosing Acute and Resolved Genotype 3 Infections.

European Association of the Study of the Liver (EASL) guidelines specify HEV RNA, as well as anti-HEV IgG and IgM as positive markers for acute HEV infection. HEV RNA assay sensitivity limitations may lead to false negative test results in patients with low levels of viremia. Moreover, anti-HEV IgM positivity is not a reliable indicator for distinguishing between acute and resolved infections given the ability of this antibody to persist several months after a resolved infection. Our study aims were to assess HEV IgG avidity for diagnosing acute and resolved infections, regardless of the anti-HEV IgM serostatus, and examine assay reliability when evaluating different genotype 3 (GT3) HEV subtypes. Patient serum samples (n = 104) were tested for HEV IgG avidity by utilizing the DIA.PRO kit on a DSX automated instrument. Among patients identified with acute HEV infections, 32 were infected with GT3: GT3c (n = 5), GT3e (n = 8), 3f (n = 17) and GT3-unsubtyped (n = 2). Avidity sensitivity was 91.2% and specificity was 100%. For patients with long-lasting anti-HEV IgM persistence, an Avidity Index >70% was observed. Thus, the DIA.PRO avidity assay may be utilized to distinguish between recently acquired and resolved HEV GT3 infections. However, for equivocal results (Avidity Index > 40-70%), HEV RNA molecular testing will be required to confirm a recent infection.

Genetic variability in E6, E7 and L1 protein of HPV81 from HIV-1 positive women in Italy.

The genetic variability of E6, E7 and L1 of HPV81 from HIV-1 positive women carrying multiple HPV infections was investigated by clonal analysis for E6 and E7. The range of maximal divergence from the prototype was 0.6%-2.6% for E6 and 1.0%-3.1% for E7. Compared to prototype HPV81, 13 and 10 mutations were identified in E6 and E7, respectively. In the pRB binding domain of E7, all HPV81 clones showed D21, as reported for prototype HPV81 and for HPV16 and 18, while G22 is reported in HPV6 and 11. In the CR3 region, CxxC motif was conserved in all but one clone. The L1 sequence of a single clone from 5 study patients was also established. The range of similarity with prototype HPV81 was 97.8%-99.2%, with 25 polymorphic sites. Two substitutions (R492K and T493S) were observed in 5/5, one (T287N) in 4/5 patients. Among L1 immune-related regions, BC loop presented T56N in 1/5, while FGb loop presented T287N in 4/5 patients. Our data indicate the presence of polymorphisms in all 3 HPV81 genes analyzed, with a certain degree of intra-patient diversity. The importance of polymorphisms on HPV81 persistence and pathogenicity needs to be addressed in longitudinal studies involving larger patient numbers.

Origin of HAV strains responsible for 2016-2017 outbreak among MSM: Viral phylodynamics in Lazio region.

In Europe HAV infection occurs mainly among specific risk groups, such as consumers of specific food. Sexual transmission of HAV has been demonstrated, particularly among Men-Who-Have-Sex-With-Men (MSM), causing MSM-specific outbreaksin Europe. Here we report a molecular epidemiologic and phylodynamic analysis on HAV sequences in Lazio (central Italy)to identify genetic background and the phylogenetic relations, and test the HAV infection dynamics during a large outbreak through phylodynamic model.Among all HAV sequences found during 2013-2018 in Lazio, low genetic diversity was observed in HAV population in 2016 and 2017, along with high frequenciesVRD_521_2016and RIVM-HAV16-090, suggesting a large expansion event of viral population. The initial expansion of both VRD_521_2016 and RIVM-HAV16-090 clusters dated back to 2012 (95% HPD:2006-2015). During the2016-2017outbreak in Lazio region, the Re peaked around mid-2016, with a value of 1.73 (95% HPD: 1.03-2.37), consistent with incidence trend of AHA cases in Lazio between 2016 and mid-2017. This study showed the magnitude of HAV outbreak in Lazio during 2016-2017, demonstrating the epidemic continuity to MSM-specific outbreak in Europe. The HAV dataset is available on interactive phylodynamic platform https://nextstrain.org to real-time update of future outbreaks.