Genomic epidemiology and phylogenetics applied to the study of SARS-CoV-2 pandemic.

The study of characteristics, prevalence and patterns of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is significant to monitor and define the status of the pandemic, helping to design and evaluate control strategies. In this setting, the continuous emergence of new variants and their dynamic of replacement underline the importance of implementing genomic epidemiology and phylogenetic methods for the molecular monitoring and surveillance of this new virus. The current profile of the pandemic can change rapidly when new variants emerge and spread, impacting epidemiology and public health in terms of prevention and treatment and making it necessary to develop new molecules and formulate vaccines. In this paper, we reviewed and synthesized the main studies on molecular genomics and phylogeny of SARS-CoV-2 during the pandemic, and highlighted their contributions to our understanding of this new emergent pathogen.

Epidemiological and Clinical Features of SARS-CoV-2 Variants Circulating between April-December 2021 in Italy.

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.

From Clinical Specimen to Whole Genome Sequencing of A(H3N2) Influenza Viruses: A Fast and Reliable High-Throughput Protocol.

(1) Background: Over the last few years, there has been growing interest in the whole genome sequencing (WGS) of rapidly mutating pathogens, such as influenza viruses (IVs), which has led us to carry out in-depth studies on viral evolution in both research and diagnostic settings. We aimed at describing and determining the validity of a WGS protocol that can obtain the complete genome sequence of A(H3N2) IVs directly from clinical specimens. (2) Methods: RNA was extracted from 80 A(H3N2)-positive respiratory specimens. A one-step RT-PCR assay, based on the use of a single set of specific primers, was used to retro-transcribe and amplify the entire IV type A genome in a single reaction, thus avoiding additional enrichment approaches and host genome removal treatments. Purified DNA was quantified; genomic libraries were prepared and sequenced by using Illumina MiSeq platform. The obtained reads were evaluated for sequence quality and read-pair length. (3) Results: All of the study specimens were successfully amplified, and the purified DNA concentration proved to be suitable for NGS (at least 0.2 ng/µL). An acceptable coverage depth for all eight genes of influenza A(H3N2) virus was obtained for 90% (72/80) of the clinical samples with viral loads >105 genome copies/mL. The mean depth of sequencing ranged from 105 to 200 reads per position, with the majority of the mean depth values being above 103 reads per position. The total turnaround time per set of 20 samples was four working days, including sequence analysis. (4) Conclusions: This fast and reliable high-throughput sequencing protocol should be used for influenza surveillance and outbreak investigation.

Evidence of SARS-CoV-2 Antibodies and RNA on Autopsy Cases in the Pre-Pandemic Period in Milan (Italy).

In this study, we analyzed blood samples obtained from 169 cadavers subjected to an autopsy from 1 October 2019 to 27 March 2020. The presence of anti-severe acute respiratory syndrome coronavirus 2 (anti-SARS-CoV-2) antibodies was searched by lateral flow immunochromatographic assay (LFIA) and ELISA tests and the SARS-CoV-2 RNA was tested in blood and available lung tissues by real-time PCR (RT-PCR) and droplet digital PCR (ddPCR). Five cases resulted in positives at the serological screening for anti-SARS-CoV-2. Three results were weakly positive for IgM while only one showed strong reactivity for IgG antibodies. The fifth subject (who died in December 2019) resulted positive for the ELISA test. The detection of SARS-CoV-2 RNA resulted in positive only in the blood and lung tissues of such cases. These data suggest that cadaveric blood may be a suitable substrate for the assessment of SARS-CoV-2 infection; moreover, they extend the observations of sporadic cases of SARS-CoV-2 infection in North Italy prior to the first confirmed cases.

Sofosbuvir Selects for Drug-Resistant Amino Acid Variants in the Zika Virus RNA-Dependent RNA-Polymerase Complex In Vitro.

The nucleotide analog sofosbuvir, licensed for the treatment of hepatitis C, recently revealed activity against the Zika virus (ZIKV) in vitro and in animal models. However, the ZIKV genetic barrier to sofosbuvir has not yet been characterized. In this study, in vitro selection experiments were performed in infected human hepatoma cell lines. Increasing drug pressure significantly delayed viral breakthrough (p = 0.029). A double mutant in the NS5 gene (V360L/V607I) emerged in 3 independent experiments at 40-80 µM sofosbuvir resulting in a 3.9 ± 0.9-fold half- maximal inhibitory concentration (IC50) shift with respect to the wild type (WT) virus. A triple mutant (C269Y/V360L/V607I), detected in one experiment at 80 µM, conferred a 6.8-fold IC50 shift with respect to the WT. Molecular dynamics simulations confirmed that the double mutant V360L/V607I impacts the binding mode of sofosbuvir, supporting its role in sofosbuvir resistance. Due to the distance from the catalytic site and to the lack of reliable structural data, the contribution of C269Y was not investigated in silico. By a combination of sequence analysis, phenotypic susceptibility testing, and molecular modeling, we characterized a double ZIKV NS5 mutant with decreased sofosbuvir susceptibility. These data add important information to the profile of sofosbuvir as a possible lead for anti-ZIKV drug development.

Time and Mode of Epidemic HCV-2 Subtypes Spreading in Europe: Phylodynamics in Italy and Albania.

Hepatitis C virus (HCV) genotype 2 causes about 10% of global infections and has the most variable circulation profile in Europe. The history of "endemic" HCV-2 subtypes has been satisfactorily reconstructed, instead there is little information about the recent spread of the "epidemic" subtypes, including HCV-2c. To investigate the origin and dispersion pathways of HCV-2c, 245 newly characterized Italian and Albanian HCV-2 NS5B sequences were aligned with 247 publicly available sequences and included in phylogeographic and phylodynamic analyses using the Bayesian framework. Our findings show that HCV-2c was the most prevalent subtype in Italy and Albania. The phylogeographic analysis suggested an African origin of HCV-2c before it reached Italy about in the 1940s. Phylodynamic analysis revealed an exponential increase in the effective number of infections and Re in Italy between the 1940s and 1960s, and in Albania between the 1990s and the early 2000s. It seems very likely that HCV-2c reached Italy from Africa at the time of the second Italian colonization but did not reach Albania until the period of dramatic migration to Italy in the 1990s. This study contributes to reconstructing the history of the spread of epidemic HCV-2 subtypes to Europe.

Evolutionary Dynamics of the Lineage 2 West Nile Virus That Caused the Largest European Epidemic: Italy 2011-2018.

Lineage 2 West Nile virus (WNV) caused a vast epidemic in Europe in 2018, with the highest incidence being recorded in Italy. To reconstruct the evolutionary dynamics and epidemiological history of the virus in Italy, 53 envelope gene and 26 complete genome sequences obtained from human and animal samples were characterised by means of next-generation sequencing. Phylogenetic analysis revealed two Italian strains originating between 2010 and 2012: clade A, which apparently became extinct in 2013-2014, and clade B, which was responsible for the 2018 epidemic. The mean genetic distances in clade B increased over time and with the distance between sampling locations. Bayesian birth-death and coalescent skyline plots of the clade B showed that the effective number of infections and the effective reproduction number (Re) increased between 2015 and 2018. Our data suggest that WNV-2 entered Italy in 2011 as a result of one or a few penetration events. Clade B differentiated mainly as a result of genetic drift and purifying selection, leading to the appearance of multiple locally circulating sub-clades for different times. Phylodynamic analysis showed a current expansion of the infection among reservoir birds and/or vectors.