Discovering and exploring the hidden diversity of human gut viruses using highly enriched virome samples.

Viruses are an abundant and crucial component of the human microbiome, but accurately discovering them via metagenomics is still challenging. Currently, the available viral reference genomes poorly represent the diversity in microbiome samples, and expanding such a set of viral references is difficult. As a result, many viruses are still undetectable through metagenomics even when considering the power of de novo metagenomic assembly and binning, as viruses lack universal markers. Here, we describe a novel approach to catalog new viral members of the human gut microbiome and show how the resulting resource improves metagenomic analyses. We retrieved >3,000 viral-like particles (VLP) enriched metagenomic samples (viromes), evaluated the efficiency of the enrichment in each sample to leverage the viromes of highest purity, and applied multiple analysis steps involving assembly and comparison with hundreds of thousands of metagenome-assembled genomes to discover new viral genomes. We reported over 162,000 viral sequences passing quality control from thousands of gut metagenomes and viromes. The great majority of the retrieved viral sequences (~94.4%) were of unknown origin, most had a CRISPR spacer matching host bacteria, and four of them could be detected in >50% of a set of 18,756 gut metagenomes we surveyed. We included the obtained collection of sequences in a new MetaPhlAn 4.1 release, which can quantify reads within a metagenome matching the known and newly uncovered viral diversity. Additionally, we released the viral database for further virome and metagenomic studies of the human microbiome.

Shallow Whole-Genome Sequencing of Aedes japonicus and Aedes koreicus from Italy and an Updated Picture of Their Evolution Based on Mitogenomics and Barcoding.

Aedes japonicus and Aedes koreicus are two invasive mosquitoes native to East Asia that are quickly establishing in temperate regions of Europe. Both species are vectors of arboviruses, but we currently lack a clear understanding of their evolution. Here, we present new short-read, shallow genome sequencing of A. japonicus and A. koreicus individuals from northern Italy, which we used for downstream phylogenetic and barcode analyses. We explored associated microbial DNA and found high occurrences of Delftia bacteria in both samples, but neither Asaia nor Wolbachia. We then assembled complete mitogenomes and used these data to infer divergence times estimating the split of A. japonicus from A. koreicus in the Oligocene, which was more recent than that previously reported using mitochondrial markers. We recover a younger age for most other nodes within Aedini and other Culicidae. COI barcoding and phylogenetic analyses indicate that A. japonicus yaeyamensis, A. japonicus amamiensis, and the two A. koreicus sampled from Europe should be considered as separate species within a monophyletic species complex. Our studies further clarify the evolution of A. japonicus and A. koreicus, and indicate the need to obtain whole-genome data from putative species in order to disentangle their complex patterns of evolution.

Spatial and temporal dynamics of West Nile virus between Africa and Europe.

It is unclear whether West Nile virus (WNV) circulates between Africa and Europe, despite numerous studies supporting an African origin and high transmission in Europe. We integrated genomic data with geographic observations and phylogenetic and phylogeographic inferences to uncover the spatial and temporal viral dynamics of WNV between these two continents. We focused our analysis towards WNV lineages 1 (L1) and 2 (L2), the most spatially widespread and pathogenic WNV lineages. Our study shows a Northern-Western African origin of L1, with back-and-forth exchanges between West Africa and Southern-Western Europe; and a Southern African origin of L2, with one main introduction from South Africa to Europe, and no back introductions observed. We also noticed a potential overlap between L1 and L2 Eastern and Western phylogeography and two Afro-Palearctic bird migratory flyways. Future studies linking avian and mosquito species susceptibility, migratory connectivity patterns, and phylogeographic inference are suggested to elucidate the dynamics of emerging viruses.

First genome-wide data from Italian European beech (Fagus sylvatica L.): Strong and ancient differentiation between Alps and Apennines.

The European beech (Fagus sylvatica L.) is one of the most widespread forest trees in Europe whose distribution and intraspecific diversity has been largely shaped by repeated glacial cycles. Previous studies, mainly based on palaeobotanical evidence and a limited set of chloroplast and nuclear genetic markers, highlighted a complex phylogeographic scenario, with southern and western Europe characterized by a rather heterogeneous genetic structure, as a result of recolonization from different glacial refugia. Despite its ecological and economic importance, the genome of this broad-leaved tree has only recently been assembled, and its intra-species genomic diversity is still largely unexplored. Here, we performed whole-genome resequencing of nine Italian beech individuals sampled from two stands located in the Alpine and Apennine mountain ranges. We investigated patterns of genetic diversity at chloroplast, mitochondrial and nuclear genomes and we used chloroplast genomes to reconstruct a temporally-resolved phylogeny. Results allowed us to test European beech differentiation on a whole-genome level and to accurately date their divergence time. Our results showed comparable, relatively high levels of genomic diversity in the two populations and highlighted a clear differentiation at chloroplast, mitochondrial and nuclear genomes. The molecular clock analysis indicated an ancient split between the Alpine and Apennine populations, occurred between the Günz and the Riss glaciations (approximately 660 kyrs ago), suggesting a long history of separation for the two gene pools. This information has important conservation implications in the context of adaptation to ongoing climate changes.

On the Origin and Propagation of the COVID-19 Outbreak in the Italian Province of Trento, a Tourist Region of Northern Italy.

BACKGROUND: Trentino is an Italian province with a tourism-based economy, bordering the regions of Lombardy and Veneto, where the two earliest and largest outbreaks of COVID-19 occurred in Italy. The earliest cases in Trentino were reported in the first week of March 2020, with most of the cases occurring in the winter sport areas in the Dolomites mountain range. The number of reported cases decreased over the summer months and was followed by a second wave in the autumn and winter of 2020. METHODS: we performed high-coverage Oxford Nanopore sequencing of 253 positive SARS-CoV-2 swabs collected in Trentino between March and December 2020. RESULTS: in this work, we analyzed genome sequences to trace the routes through which the virus entered the area, and assessed whether the autumnal resurgence could be attributed to lineages persisting undetected during summer, or as a consequence of new introductions. CONCLUSIONS: Comparing the draft genomes analyzed with a large selection of European sequences retrieved from GISAID we found that multiple introductions of the virus occurred at the early stage of the epidemics; the two epidemic waves were unrelated; the second wave was due to reintroductions of the virus in summer when traveling restrictions were uplifted.

Epidemiological and Evolutionary Analysis of West Nile Virus Lineage 2 in Italy.

West Nile virus (WNV) is a mosquito-borne virus potentially causing serious illness in humans and other animals. Since 2004, several studies have highlighted the progressive spread of WNV Lineage 2 (L2) in Europe, with Italy being one of the countries with the highest number of cases of West Nile disease reported. In this paper, we give an overview of the epidemiological and genetic features characterising the spread and evolution of WNV L2 in Italy, leveraging data obtained from national surveillance activities between 2011 and 2021, including 46 newly assembled genomes that were analysed under both phylogeographic and phylodynamic frameworks. In addition, to better understand the seasonal patterns of the virus, we used a machine learning model predicting areas at high-risk of WNV spread. Our results show a progressive increase in WNV L2 in Italy, clarifying the dynamics of interregional circulation, with no significant introductions from other countries in recent years. Moreover, the predicting model identified the presence of suitable conditions for the 2022 earlier and wider spread of WNV in Italy, underlining the importance of using quantitative models for early warning detection of WNV outbreaks. Taken together, these findings can be used as a reference to develop new strategies to mitigate the impact of the pathogen on human and other animal health in endemic areas and new regions.

On the correct interpretation of similarity index in codon usage studies: Comparison with four other metrics and implications for Zika and West Nile virus.

Viruses can modify the usage of their synonymous codons to mimic that of their hosts: investigation of codon usage may therefore help increase our understanding of virus biology. The similarity index (SiD or (D(A,B)) is a recent method for estimating the resemblance of the overall codon usage of the hosts to that of a certain virus. Here we report that this index has been systematically misinterpreted in the literature. We first review its mathematical meaning to show that a high value of SiD indicates low codon usage similarity, and not the opposite as frequently reported. We further compared SiD with four other metrics - nCAI, nRCDI, ITE and tAI - using Zika and West Nile arboviruses in relation to various vectors and hosts. We highlight some differences between the five metrics and present a first bioinformatic evidence that Zika virus is poorly adapted to Culex quinquefasciatus, in accordance with the known biology of this virus. We conclude that SiD, if correctly interpreted, is a useful tool for codon studies, but we advocate that the results from this and other methods are context and time dependent and should be taken with care. A script to automatically calculate SiD from two codon usage tables is available at https://github.com/andrea-silverj/SiDcal.