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.

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.

Objectives: We evaluated natural resistance to the new antiretroviral fostemsavir and its potential association with other HIV-1 gp120 polymorphisms.Methods: A total of 1997 HIV-1 B subtype gp120 sequences from the Los Alamos HIV Database were analysed for mutation prevalence at fostemsavir resistance-associated positions and potential association with other gp120 polymorphisms. The role of each fostemsavir resistance-related position and the correlated gp120 mutations, both in protein stability and in reducing the binding affinity between antibody and/or T cell lymphocyte epitopes and the MHC molecules, was estimated.Results: The prevalence of fostemsavir resistance mutations was as follows: L116Q (0.05%), S375H/M/T (0.55%/1.35%/17.73%, the latter being far less relevant in determining resistance), M426L (7.56%), M434I (4.21%) and M475I (1.65%). Additionally, the M426R polymorphism had a prevalence of 16.32%. A significantly higher prevalence in X4 viruses versus R5 viruses was found only for S375M (0.69% versus 3.93%, P = 0.009) and S375T (16.60% versus 22.11%, P = 0.030). Some fostemsavirv resistance positions positively and significantly correlated with specific gp120 polymorphisms: S375T with I371V; S375M with L134W, I154V and I323T; M475I with K322A; and M426R with G167N, K192T and S195N. The topology of the dendrogram suggested the existence of three distinct clusters (bootstrap ≥0.98) involving these fostemsavir resistance mutations and gp120 polymorphisms. Interestingly, all clustered mutations are localized in class I/II-restricted T cell/antibody epitopes, suggesting a potential role in immune HIV escape.Conclusions: A low prevalence of known fostemsavir resistance mutations was found in the HIV-1 B subtype. The detection of novel HIV-1 gp120 polymorphisms potentially relevant for fostemsavir resistance deserves new in-depth in vitro investigations.

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 (Allplex^TM 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.