Search

Equine hepacivirus (EqHV, Flaviviridae, hepacivirus) is a small, enveloped RNA virus generally causing sub-clinical hepatitis with occasional fatalities. EqHV is reported in equids worldwide, but for Italy data are limited. To address this, a survey study was set up to estimate prevalence at a national level and among different production categories (equestrian; competition; work and meat; reproduction) and national macro-regions (North, Central, South, and Islands). Data obtained testing 1801 horse serum samples by Real-Time RT PCR were compared within the categories and regions. The NS3 fragment of the PCR-positive samples was sequenced by Sanger protocol for phylogenetic and mutational analysis. The tertiary structure of the NS3 protein was also assessed. The estimated national prevalence was 4.27% [1.97–6.59, 95% CI] and no statistical differences were detected among production categories and macro-regions. The phylogenesis confirmed the distribution in Italy of the three known EqHV subtypes, also suggesting a possible fourth sub-type that, however, requires further confirmation. Mutational profiles that could also affect the NS3 binding affinity to the viral RNA were detected. The present paper demonstrates that EqHV should be included in diagnostic protocols when investigating causes of hepatitis, and in quality control protocols for blood derived products due to its parental transmission.

ABSTRACTSpecific HBsAg mutations are known to hamper HBsAg recognition by neutralizing antibodies thus challenging HBV-vaccination efficacy. Nevertheless, information on their impact and spreading over time is limited. Here, we characterize the circulation of vaccine-escape mutations from 2005 to 2019 and their correlation with virological parameters in a large cohort of patients infected with HBV genotype-D (N = 947), dominant in Europe. Overall, 17.7% of patients harbours ≥1 vaccine-escape mutation with the highest prevalence in subgenotype-D3. Notably, complex profiles (characterized by ≥2 vaccine-escape mutations) are revealed in 3.1% of patients with a prevalence rising from 0.4% in 2005–2009 to 3.0% in 2010–2014 and 5.1% in 2015–2019 (P = 0.007) (OR[95%CI]:11.04[1.42–85.58], P = 0.02, by multivariable-analysis). The presence of complex profiles correlates with lower HBsAg-levels (median[IQR]:40[0–2905]IU/mL for complex profiles vs 2078[115–6037]IU/ml and 1881[410–7622]IU/mL for single or no vaccine-escape mutation [P

Abstract Response to ritonavir-boosted-protease inhibitors (PI/r)-based regimen is associated with some Gag mutations among HIV-1 B-clade. There is limited data on Gag mutations and their covariation with mutations in protease among HIV-1 non-B-clades at PI/r-based treatment failure. Thus, we characterized Gag mutations present in isolates from HIV-1 infected individuals treated with a PI/r-regimen (n = 143) and compared them with those obtained from individuals not treated with PI/r (ART-naïve [n = 101] or reverse transcriptase inhibitors (RTI) treated [n = 118]). The most frequent HIV-1 subtypes were CRF02_AG (54.69%), A (13.53%), D (6.35%) and G (4.69%). Eighteen Gag mutations showed a significantly higher prevalence in PI/r-treated isolates compared to ART-naïve (p

(1) Background: SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target to fight COVID-19, and many RdRp inhibitors nucleotide/nucleoside analogs, such as remdesivir, have been identified or are in clinical studies. However, the appearance of resistant mutations could reduce their efficacy. In the present work, we structurally evaluated the impact of RdRp mutations found at baseline in 39 patients treated with remdesivir and associated with a different degree of antiviral response in vivo. (2) Methods: A refined bioinformatics approach was applied to assign SARS-CoV-2 clade and lineage, and to define RdRp mutational profiles. In line with such a method, the same mutations were built and analyzed by combining docking and thermodynamics evaluations with both molecular dynamics and representative pharmacophore models. (3) Results: Clinical studies revealed that patients bearing the most prevalent triple mutant P323L+671S+M899I, which was present in 41% of patients, or the more complex mutational profile P323L+G671S+L838I+D738Y+K91E, which was found with a prevalence of 2.6%, showed a delayed reduced response to remdesivir, as confirmed by the increase in SARS-CoV-2 viral load and by a reduced theoretical binding affinity versus RdRp (ΔGbindWT = −122.70 kcal/mol; ΔGbindP323L+671S+M899I = −84.78 kcal/mol; ΔGbindP323L+G671S+L838I+D738Y+K91E = −96.74 kcal/mol). Combined computational approaches helped to rationalize such clinical observations, offering a mechanistic understanding of the allosteric effects of mutants on the global motions of the viral RNA synthesis machine and in the changes of the interactions patterns of remdesivir during its binding.

In this study, we provided a retrospective overview in order to better define SARS-CoV-2 variants circulating in Italy during the first two years of the pandemic, by characterizing the spike mutational profiles and their association with viral load (expressed as ct values), N-glycosylation pattern, hospitalization and vaccination. Next-generation sequencing (NGS) data were obtained from 607 individuals (among them, 298 vaccinated and/or 199 hospitalized). Different rates of hospitalization were observed over time and among variants of concern (VOCs), both in the overall population and in vaccinated individuals (Alpha: 40.7% and 31.3%, Beta: 0%, Gamma: 36.5% and 44.4%, Delta: 37.8% and 40.2% and Omicron: 11.2% and 7.1%, respectively, both p-values < 0.001). Approximately 32% of VOC-infected individuals showed at least one atypical major spike mutation (intra-prevalence > 90%), with a distribution differing among the strains (22.9% in Alpha, 14.3% in Beta, 41.8% in Gamma, 46.5% in Delta and 15.4% in Omicron, p-value < 0.001). Overall, significantly less atypical variability was observed in vaccinated individuals than unvaccinated individuals; nevertheless, vaccinated people who needed hospitalization showed an increase in atypical variability compared to vaccinated people that did not need hospitalization. Only 5/607 samples showed a different putative N-glycosylation pattern, four within the Delta VOC and one within the Omicron BA.2.52 sublineage. Interestingly, atypical minor mutations (intra-prevalence < 20%) were associated with higher Ct values and a longer duration of infection. Our study reports updated information on the temporal circulation of SARS-CoV-2 variants circulating in Central Italy and their association with hospitalization and vaccination. The results underline how SARS-CoV-2 has changed over time and how the vaccination strategy has contributed to reducing severity and hospitalization for this infection in Italy.

ABSTRACT The process of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genetic diversification is still ongoing and has very recently led to the emergence of a new variant of concern (VOC), defined as Omicron or B.1.1.529. Omicron VOC is the most divergent variant identified so far and has generated immediate concern for its potential capability to increase SARS-CoV-2 transmissibility and, more worryingly, to escape therapeutic and vaccine-induced antibodies. Nevertheless, a clear definition of the Omicron VOC mutational spectrum is still missing. Herein, we provide a comprehensive definition and functional characterization (in terms of infectivity and/or antigenicity) of mutations characterizing the Omicron VOC. In particular, 887,475 SARS-CoV-2 Omicron VOC whole-genome sequences were retrieved from the GISAID database and used to precisely define its specific patterns of mutations across the different viral proteins. In addition, the functional characterization of Omicron VOC spike mutations was finely discussed according to published manuscripts. Lastly, residues characterizing the Omicron VOC and the previous four VOCs (Alpha, Beta, Gamma, and Delta) were mapped on the three-dimensional structure of the SARS-CoV-2 spike protein to assess their localization in the different spike domains. Overall, our study will assist with deciphering the Omicron VOC mutational profile and will shed more light on its clinical implications. This is critical considering that Omicron VOC is currently the predominant variant worldwide. IMPORTANCE The Omicron variant of concern (VOC) has a peculiar spectrum of mutations characterized by the acquisition of mutations or deletions rarely detected in previously identified variants, particularly in the spike glycoprotein. Such mutations, mostly residing in the receptor-binding domain, could play a pivotal role in enhancing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infectivity (by increasing binding affinity for ACE2), jeopardizing spike recognition by therapeutic and vaccine-induced antibodies and causing diagnostic assay failure. To our knowledge, this is one of the first exhaustive descriptions of newly emerged mutations underlying the Omicron VOC and its biological and clinical implications.

Abstract Patients with HBeAg-negative chronic infection (CI) have not been extensively studied because of low viremia. The HBx protein, encoded by HBX, has a key role in viral replication. Here, we analyzed the viral quasispecies at the 5′ end of HBX in CI patients and compared it with that of patients in other clinical stages. Fifty-eight HBeAg-negative patients were included: 16 CI, 19 chronic hepatitis B, 16 hepatocellular carcinoma and 6 liver cirrhosis. Quasispecies complexity and conservation were determined in the region between nucleotides 1255 and 1611. Amino acid changes detected were tested in vitro. CI patients showed higher complexity in terms of mutation frequency and nucleotide diversity and higher quasispecies conservation (p

ABSTRACTIncreasing evidences suggest that HBsAg-production varies across HBV-genotypes. HBsAg C-terminus plays a crucial role for HBsAg-secretion. Here, we evaluate HBsAg-levels in different HBV-genotypes in HBeAg-negative chronic infection, the correlation of specific mutations in HBsAg C-terminus with HBsAg-levels in-vivo, their impact on HBsAg-secretion in-vitro and on structural stability in-silico.HBsAg-levels were investigated in 323 drug-naïve HBeAg-negative patients chronically infected with HBV genotype-D(N = 228), -A(N = 65) and -E(N = 30). Genotype-D was characterized by HBsAg-levels lower than genotype-A and -E (3.3[2.7–3.8]IU/ml; 3.8[3.5–4.2]IU/ml and 3.9[3.7–4.2]IU/ml, P

ABSTRACT Since the beginning of the coronavirus disease 2019 (COVID-19) pandemic caused by it, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been undergoing a genetic diversification leading to the emergence of new variants. Nevertheless, a clear definition of the genetic signatures underlying the circulating variants is still missing. Here, we provide a comprehensive insight into mutational profiles characterizing each SARS-CoV-2 variant, focusing on spike mutations known to modulate viral infectivity and/or antigenicity. We focused on variants and on specific relevant mutations reported by GISAID, Nextstrain, Outbreak.info, Pango, and Stanford database websites that were associated with any clinical/diagnostic impact, according to published manuscripts. Furthermore, 1,223,338 full-length high-quality SARS-CoV-2 genome sequences were retrieved from GISAID and used to accurately define the specific mutational patterns in each variant. Finally, mutations were mapped on the three-dimensional structure of the SARS-CoV-2 spike protein to assess their localization in the different spike domains. Overall, this review sheds light and assists in defining the genetic signatures characterizing the currently circulating variants and their clinical relevance. IMPORTANCE Since the emergence of SARS-CoV-2, several recurrent mutations, particularly in the spike protein, arose during human-to-human transmission or spillover events between humans and animals, generating distinct worrisome variants of concern (VOCs) or of interest (VOIs), designated as such due to their clinical and diagnostic impacts. Characterizing these variants and their related mutations is important in tracking SAR-CoV-2 evolution and understanding the efficacy of vaccines and therapeutics based on monoclonal antibodies, convalescent-phase sera, and direct antivirals. Our study provides a comprehensive survey of the mutational profiles characterizing the important SARS-CoV-2 variants, focusing on spike mutations and highlighting other protein mutations.

SARS-CoV-2 is constantly evolving, leading to new variants. We analysed data from 4400 SARS-CoV-2-positive samples in order to pursue epidemiological variant surveillance and to evaluate their impact on public health in Italy in the period of April–December 2021. The main circulating strain (76.2%) was the Delta variant, followed by the Alpha (13.3%), the Omicron (5.3%), and the Gamma variants (2.9%). The B.1.1 lineages, Eta, Beta, Iota, Mu, and Kappa variants, represented around 1% of cases. There were 48.2% of subjects who had not been vaccinated, and they had a lower median age compared to the vaccinated subjects (47 vs. 61 years). An increasing number of infections in the vaccinated subjects were observed over time, with the highest proportion in November (85.2%). The variants correlated with clinical status; the largest proportion of symptomatic patients (59.6%) was observed with the Delta variant, while subjects harbouring the Gamma variant showed the highest proportion of asymptomatic infection (21.6%), albeit also deaths (5.4%). The Omicron variant was only found in the vaccinated subjects, of which 47% had been hospitalised. The diffusivity and pathogenicity associated with the different SARS-CoV-2 variants are likely to have relevant public health implications, both at the national and international levels. Our study provides data on the rapid changes in the epidemiological landscape of the SARS-CoV-2 variants in Italy.

The performance of diagnostic polymerase chain reaction (PCR) assays can be impacted by SARS-CoV-2 variability as this is dependent on the full complementarity between PCR primers/probes and viral target templates. Here, we investigate the genetic variability of SARS-CoV-2 regions recognized by primers/probes utilized by PCR diagnostic assays based on nucleotide mismatching analysis. We evaluated the genetic variation in the binding regions of 73 primers/probes targeting the Nucleocapsid (N, N = 36), Spike (S, N = 22), and RNA-dependent RNA-polymerase/Helicase (RdRp/Hel, N = 15) of the publicly available PCR-based assays. Over 4.9 million high-quality SARS-CoV-2 genome sequences were retrieved from GISAID and were divided into group-A (all except Omicron, >4.2 million) and group-B (only Omicron, >558 thousand). In group-A sequences, a large range of variability in primers/probes binding regions in most PCR assays was observed. Particularly, 87.7% (64/73) of primers/probes displayed ≥1 mismatch with their viral targets, while 8.2% (6/73) contained ≥2 mismatches and 2.7% (2/73) contained ≥3 mismatches. In group-B sequences, 32.9% (24/73) of primers/probes were characterized by ≥1 mismatch, 13.7% (10/73) by ≥2 mismatches, and 5.5% (4/73) by ≥3 mismatches. The high rate of single and multiple mismatches- found in the target regions of molecular assays used worldwide for SARS-CoV-2 diagnosis reinforces the need to optimize and constantly update these assays according to SARS-CoV-2 genetic evolution and the future emergence of novel variants.

Herein, we report a case of an Italian male infected by Delta sublineage AY.4 harboring an atypical deletion, leading to a N gene target failure (NGTF) by a commercial molecular assay for SARS-CoV-2 diagnosis (AllplexTM SARS-CoV-2 Assay, Seegene). A 59-year-old unvaccinated patient was hospitalized for pulmonary embolism, with first negative results obtained by both molecular and antigen tests. After several days of viral negativity, he presented positive results for E and RdRP/S genes, but negative in N gene. Negativity in N gene was repeatedly confirmed in the following days. Suspecting an infection by the Omicron variant, SARS-CoV-2 genome sequencing was rapidly performed from nasopharyngeal swab by MiSeq and revealed the presence of the Delta sublineage AY.4 variant with an atypical deletion of six nucleotides, leading to G214-G215 deletion in the Nucleocapsid, thus responsible for NGTF. The analysis of GISAID sequences (N = 2,618,373 12 January 2022) showed that G214-G215 deletion is rarely occurring in most circulating Delta lineages and sublineages in the globe and Europe, with an overall prevalence never exceeding 0.2%. Hence, this study highlights the importance to perform SARS-CoV-2 sequencing and to characterize novel mutations/deletions that could jeopardize the proper interpretation of molecular diagnostic tests. Based on these assumptions, the role of deletions in the recently identified Omicron variant deserves further investigation.

Background: If analytical antiretroviral-treatment (ART) interruption (ATI) might significantly impact quantitative or qualitative peripheral-total HIV-DNA is still debated. Methods: Six chronically HIV-1 infected patients enrolled in APACHE-study were analysed for peripheral-total HIV-DNA and residual viremia, major-resistance-mutations (MRMs) and C2-V3-C3 evolution at pre-ATI (T1), during ATI (T2) and at achievement of virological success after ART-resumption (post-ATI, T3). These data were obtained at three comparable time-points in five chronically HIV-1 infected patients on suppressive ART for ≥1 year, enrolled in MODAt-study. Results: At T1, APACHE and MODAt individuals had similar peripheral-total HIV-DNA and residual viremia (p = 0.792 and 0.662, respectively), and no significant changes for these parameters were observed between T1 and T3 in both groups. At T1, 4/6 APACHE and 2/5 MODAt carried HIV-DNA MRMs. MRMs disappeared at T3 in 3/4 APACHE. All disappearing MRMs were characterized by T1 intra-patient prevalence p = 0.04), while no significant changes were found in MODAt. Accordingly, maximum likelihood trees (bootstrap > 70%) and genealogical sorting indices (GSI > 0.50 with p-value < 0.05) showed that T1 C2-V3-C3 DNA sequences were distinct from T2 and T3 viruses in 4/6 APACHE. Virus populations at all three time-points were highly interspersed in MODAt. Conclusions: This pilot study indicates that short ATI does not alter peripheral-total HIV-DNA burden and residual viremia, but in some cases could cause a genetic diversification of peripheral viral reservoir in term of both MRMs rearrangement and viral evolution.

HCV is an important cause of hepatocellular carcinoma (HCC). HCV NS5A domain-1 interacts with cellular proteins inducing pro-oncogenic pathways. Thus, we explore genetic variations in NS5A domain-1 and their association with HCC, by analyzing 188 NS5A sequences from HCV genotype-1b infected DAA-naïve cirrhotic patients: 34 with HCC and 154 without HCC. Specific NS5A mutations significantly correlate with HCC: S3T (8.8% vs. 1.3%, p = 0.01), T122M (8.8% vs. 0.0%, p < 0.001), M133I (20.6% vs. 3.9%, p < 0.001), and Q181E (11.8% vs. 0.6%, p < 0.001). By multivariable analysis, the presence of >1 of them independently correlates with HCC (OR (95%CI): 21.8 (5.7–82.3); p < 0.001). Focusing on HCC-group, the presence of these mutations correlates with higher viremia (median (IQR): 5.7 (5.4–6.2) log IU/mL vs. 5.3 (4.4–5.6) log IU/mL, p = 0.02) and lower ALT (35 (30–71) vs. 83 (48–108) U/L, p = 0.004), suggesting a role in enhancing viral fitness without affecting necroinflammation. Notably, these mutations reside in NS5A regions known to interact with cellular proteins crucial for cell-cycle regulation (p53, p85-PIK3, and β-catenin), and introduce additional phosphorylation sites, a phenomenon known to ameliorate NS5A interaction with cellular proteins. Overall, these results provide a focus for further investigations on molecular bases of HCV-mediated oncogenesis. The role of theseNS5A domain-1 mutations in triggering pro-oncogenic stimuli that can persist also despite achievement of sustained virological response deserves further investigation.

Background: The aim of the study was to investigate the intra-host variability through next-generation-sequencing (NGS) of the NS5A-gene in nosocomial transmission-clusters observed in two Italian hospitals among hepatitis C virus (HCV)-genotype-1b infected patients. Methods: HCV-sequencing was performed by Sanger-sequencing (NS3 + NS5A + NS5B) and by NGS (NS5A, MiSeq-Illumina) in 15 HCV-1b infected patients [five acute with onco-hematologic-disease and 10 (4/6 acute/chronic) with β-thalassemia]. Resistance-associated-substitutions (RAS) were analysed by Geno2pheno-algorithm. Nucleotide-sequence-variability (NSV, at 1%, 2%, 5%, 10% and 15% NGS-cutoffs) and Shannon entropy were estimated. Phylogenetic analysis was performed by Mega6-software and Bayesian-analysis. Results: Phylogenetic analysis showed five transmission-clusters: one involving four HCV-acute onco-hematologic-patients; one involving three HCV-chronic β-thalassemia-patients and three involving both HCV-acute and chronic β-thalassemia-patients. The NS5A-RAS Y93H was found in seven patients, distributed differently among chronic/acute patients involved in the same transmission-clusters, independently from the host-genetic IL-28-polymorphism. The intra-host NSV was higher in chronic-patients versus acute-patients, at all cutoffs analyzed (p < 0.05). Even though Shannon-entropy was higher in chronic-patients, significantly higher values were observed only in chronic β-thalassemia-patients versus acute β-thalassemia-patients (p = 0.01). Conclusions: In nosocomial HCV transmission-clusters, the intra-host HCV quasispecies divergence in patients with acute-infection was very low in comparison to that in chronic-infection. The NS5A-RAS Y93H was often transmitted and distributed differently within the same transmission-clusters, independently from the IL-28-polymorphism.

Chronic HBV + HDV infection is associated with greater risk of liver fibrosis, earlier hepatic decompensation, and liver cirrhosis hepatocellular carcinoma compared to HBV mono-infection. However, to-date no direct anti-HDV drugs are available in clinical practice. Here, we identified conserved and variable regions in HBsAg and HDAg domains in HBV + HDV infection, a critical finding for the design of innovative therapeutic agents. The extent of amino-acid variability was measured by Shannon-Entropy (Sn) in HBsAg genotype-d sequences from 31 HBV + HDV infected and 62 HBV mono-infected patients (comparable for demographics and virological-parameters), and in 47 HDAg genotype-1 sequences. Positions with Sn = 0 were defined as conserved. The percentage of conserved HBsAg-positions was significantly higher in HBV + HDV infection than HBV mono-infection (p = 0.001). Results were confirmed after stratification for HBeAg-status and patients’ age. A Sn = 0 at specific positions in the C-terminus HBsAg were correlated with higher HDV-RNA, suggesting that conservation of these positions can preserve HDV-fitness. Conversely, HDAg was characterized by a lower percentage of conserved-residues than HBsAg (p < 0.001), indicating higher functional plasticity. Furthermore, specific HDAg-mutations were significantly correlated with higher HDV-RNA, suggesting a role in conferring HDV replicative-advantage. Among HDAg-domains, only the virus-assembly signal exhibited a high genetic conservation (75% of conserved-residues). In conclusion, HDV can constrain HBsAg genetic evolution to preserve its fitness. The identification of conserved regions in HDAg poses the basis for designing innovative targets against HDV-infection.