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.

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.

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.

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.

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.

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.

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.

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.

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.

Human Papillomavirus Infections in Cervical Samples From HIV-Positive Women: Evaluation of the Presence of the Nonavalent HPV Genotypes and Genetic Diversity.

Non-nonavalent vaccine (9v) Human papillomavirus (HPV) types have been shown to have high prevalence among HIV-positive women. Here, 1444 cervical samples were tested for HPV DNA positivity. Co-infections of the 9v HPV types with other HPV types were evaluated. The HPV81 L1 and L2 genes were used to investigate the genetic variability of antigenic epitopes. HPV-positive samples were genotyped using the HPVCLART2 assay. The L1 and L2 protein sequences were analyzed using a self-optimized prediction method to predict their secondary structure. Co-occurrence probabilities of the 9v HPV types were calculated. Non9v types represented 49% of the HPV infections; 31.2% of the non9v HPV types were among the low-grade squamous intraepithelial lesion samples, and 27.3% among the high-grade squamous intraepithelial lesion samples, and several genotypes were low risk. The co-occurrence of 9v HPV types with the other genotypes was not correlated with the filogenetic distance. HPV81 showed an amino-acid substitution within the BC loop (N75Q) and the FGb loop (T315N). In the L2 protein, all of the mutations were located outside antigenic sites. The weak cross-protection of the 9v types suggests the relevance of a sustainable and effective screening program, which should be implemented by HPV DNA testing that does not include only high-risk types.

Clinical Characteristics of Acute Hepatitis E and Their Correlation with HEV Genotype 3 Subtypes in Italy.

Genotype 3 (GT3) is responsible for most European autochthonous hepatitis E virus (HEV) infections. This study analyzed circulating genotypes and GT3 subtypes in the Lazio region, Italy, between 2011 and 2019, as well as their pathogenic characteristics. Of the 64 evaluable HEV GT3 patient-derived sequences, identified subtypes included GT3f (n = 36), GT3e (n = 15), GT3c (n = 9), GT3a (n = 1) and three unsubtyped GT3 sequences. GT3c strains were similar to Dutch sequences (96.8-98.1% identity), GT3e strains showed high similarity (96.8%) with a United Kingdom sequence, while the most related sequences to GT3f Italian strains were isolated in France, Belgium and Japan. One sequence was closely related to another Italian strain isolated in raw sewage in 2016. The liver functioning test median values for 56 evaluable GT3 patients were: alanine aminotransferase (ALT), 461 (range 52-4835 U/L); aspartate aminotransferase (AST), 659 (range 64-6588 U/L); and total bilirubin, 3.49 (range 0.4-33 mg/dL). The median HEV RNA viral load for 26 evaluable GT3 patients was 42,240 IU/mL (range 5680-895,490 IU/mL). Of the 37 GT3 patients with available clinical information, no correlation was observed between HEV clinical manifestations and GT3 subtype. HEV symptoms were comparable among GT3c/e/f patients across most analyzed categories except for epigastric pain, which occurred more frequently in patients with HEV GT3e (75%) than in patients with GT3c (50%) or GT3f (19%) (p = 0.01). Additionally, patients with HEV GT3c exhibited significantly higher median international normalized ratio (INR) than patients with GT3e and GT3f (p = 0.033). The severity of GT3 acute hepatitis E was not linked to HEV RNA viral load or to the GT3 subtype.