Suppression of a SARS-CoV-2 outbreak in the Italian municipality of Vo'.

On 21 February 2020, a resident of the municipality of Vo', a small town near Padua (Italy), died of pneumonia due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection1. This was the first coronavirus disease 19 (COVID-19)-related death detected in Italy since the detection of SARS-CoV-2 in the Chinese city of Wuhan, Hubei province2. In response, the regional authorities imposed the lockdown of the whole municipality for 14 days3. Here we collected information on the demography, clinical presentation, hospitalization, contact network and the presence of SARS-CoV-2 infection in nasopharyngeal swabs for 85.9% and 71.5% of the population of Vo' at two consecutive time points. From the first survey, which was conducted around the time the town lockdown started, we found a prevalence of infection of 2.6% (95% confidence interval (CI): 2.1-3.3%). From the second survey, which was conducted at the end of the lockdown, we found a prevalence of 1.2% (95% CI: 0.8-1.8%). Notably, 42.5% (95% CI: 31.5-54.6%) of the confirmed SARS-CoV-2 infections detected across the two surveys were asymptomatic (that is, did not have symptoms at the time of swab testing and did not develop symptoms afterwards). The mean serial interval was 7.2 days (95% CI: 5.9-9.6). We found no statistically significant difference in the viral load of symptomatic versus asymptomatic infections (P = 0.62 and 0.74 for E and RdRp genes, respectively, exact Wilcoxon-Mann-Whitney test). This study sheds light on the frequency of asymptomatic SARS-CoV-2 infection, their infectivity (as measured by the viral load) and provides insights into its transmission dynamics and the efficacy of the implemented control measures.

FunTaxIS-lite: a simple and light solution to investigate protein functions in all living organisms.

MOTIVATION: Defining the full domain of protein functions belonging to an organism is a complex challenge that is due to the huge heterogeneity of the taxonomy, where single or small groups of species can bear unique functional characteristics. FunTaxIS-lite provides a solution to this challenge by determining taxon-based constraints on Gene Ontology (GO) terms, which specify the functions that an organism can or cannot perform. The tool employs a set of rules to generate and spread the constraints across both the taxon hierarchy and the GO graph. RESULTS: The taxon-based constraints produced by FunTaxIS-lite extend those provided by the Gene Ontology Consortium by an average of 300%. The implementation of these rules significantly reduces errors in function predictions made by automatic algorithms and can assist in correcting inconsistent protein annotations in databases. AVAILABILITY AND IMPLEMENTATION: FunTaxIS-lite is available on https://www.medcomp.medicina.unipd.it/funtaxis-lite and from https://github.com/MedCompUnipd/FunTaxIS-lite.

ITSoneWB: profiling global taxonomic diversity of eukaryotic communities on Galaxy.

SUMMARY: ITSoneWB (ITSone WorkBench) is a Galaxy-based bioinformatic environment where comprehensive and high-quality reference data are connected with established pipelines and new tools in an automated and easy-to-use service targeted at global taxonomic analysis of eukaryotic communities based on Internal Transcribed Spacer 1 variants high-throughput sequencing. AVAILABILITY AND IMPLEMENTATION: ITSoneWB has been deployed on the INFN-Bari ReCaS cloud facility and is freely available on the web at http://itsonewb.cloud.ba.infn.it/galaxy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

QPARSE: searching for long-looped or multimeric G-quadruplexes potentially distinctive and druggable.

MOTIVATION: G-quadruplexes (G4s) are non-canonical nucleic acid conformations that are widespread in all kingdoms of life and are emerging as important regulators both in RNA and DNA. Recently, two new higher-order architectures have been reported: adjacent interacting G4s and G4s with stable long loops forming stem-loop structures. As there are no specialized tools to identify these conformations, we developed QPARSE. RESULTS: QPARSE can exhaustively search for degenerate potential quadruplex-forming sequences (PQSs) containing bulges and/or mismatches at genomic level, as well as either multimeric or long-looped PQS (MPQS and LLPQS, respectively). While its assessment versus known reference datasets is comparable with the state-of-the-art, what is more interesting is its performance in the identification of MPQS and LLPQS that present algorithms are not designed to search for. We report a comprehensive analysis of MPQS in human gene promoters and the analysis of LLPQS on three experimentally validated case studies from HIV-1, BCL2 and hTERT. AVAILABILITY AND IMPLEMENTATION: QPARSE is freely accessible on the web at http://www.medcomp.medicina.unipd.it/qparse/index or downloadable from github as a python 2.7 program https://github.com/B3rse/qparse. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

NeSSie: a tool for the identification of approximate DNA sequence symmetries.

Motivation: Non-B DNA conformations play an important role in genomic rearrangements, structural three-dimensional organization and gene regulation. Many non-B DNA structures show symmetrical properties as palindromes and mirrors that can form hairpins, cruciform structures or triplexes. A comprehensive tool, capable to perform a fast genome wide search for exact and degenerate symmetrical patterns, is needed for further investigating nucleotide tracts potentially forming non-B DNA structures. Results: We developed NeSSie, an easily customizable C/C++ 64-bit library and tool, based on dynamic programming, to quickly scan for perfect and degenerate DNA palindromes, mirrors and potential triplex forming patterns. In addition, the tool computes linguistic complexity and Shannon entropy measures to verify the repetitive nature of the DNA regions enriched in these motifs. As a case study, the analysis of the Mycobacterium bovis genome is presented. Availability and implementation: http://www.medcomp.medicina.unipd.it/main_site/doku.php? id=nessie and https://github.com/B3rse/nessie. Supplementary information: Supplementary data are available at Bioinformatics online.

G-quadruplex forming sequences in the genome of all known human viruses: A comprehensive guide.

G-quadruplexes are non-canonical nucleic-acid structures that control transcription, replication, and recombination in organisms. G-quadruplexes are present in eukaryotes, prokaryotes, and viruses. In the latter, mounting evidence indicates their key biological activity. Since data on viruses are scattered, we here present a comprehensive analysis of potential quadruplex-forming sequences (PQS) in the genome of all known viruses that can infect humans. We show that occurrence and location of PQSs are features characteristic of each virus class and family. Our statistical analysis proves that their presence within the viral genome is orderly arranged, as indicated by the possibility to correctly assign up to two-thirds of viruses to their exact class based on the PQS classification. For each virus we provide: i) the list of all PQS present in the genome (positive and negative strands), ii) their position in the viral genome, iii) the degree of conservation among strains of each PQS in its genome context, iv) the statistical significance of PQS abundance. This information is accessible from a database to allow the easy navigation of the results: http://www.medcomp.medicina.unipd.it/main_site/doku.php?id=g4virus. The availability of these data will greatly expedite research on G-quadruplex in viruses, with the possibility to accelerate finding therapeutic opportunities to numerous and some fearsome human diseases.

The complete genome sequence analysis of West Nile virus strains isolated in Slovakia (central Europe).

The present study reports the first complete genome sequence analysis of West Nile virus (WNV) strains isolated from brain samples from raptors. The results prove the circulation of closely related WNV lineage II strains in central Europe and genetic analysis revealed seven amino acid substitutions in structural (PrM3, E159 and E231) and in non-structural (NS1109, NS5259, NS5310 and NS5600) proteins. Observed amino acid substitutions Phe3 and Ser231 were common only within the lineage VII Koutango strain isolated from Rhipicephalus guilhoni tick in Senegal. Further research could reveal whether these substitutions influence the biological properties of WNV, including virulence and neuroinvasiveness.

Optimizing PCR primers targeting the bacterial 16S ribosomal RNA gene.

BACKGROUND: Targeted amplicon sequencing of the 16S ribosomal RNA gene is one of the key tools for studying microbial diversity. The accuracy of this approach strongly depends on the choice of primer pairs and, in particular, on the balance between efficiency, specificity and sensitivity in the amplification of the different bacterial 16S sequences contained in a sample. There is thus the need for computational methods to design optimal bacterial 16S primers able to take into account the knowledge provided by the new sequencing technologies. RESULTS: We propose here a computational method for optimizing the choice of primer sets, based on multi-objective optimization, which simultaneously: 1) maximizes efficiency and specificity of target amplification; 2) maximizes the number of different bacterial 16S sequences matched by at least one primer; 3) minimizes the differences in the number of primers matching each bacterial 16S sequence. Our algorithm can be applied to any desired amplicon length without affecting computational performance. The source code of the developed algorithm is released as the mopo16S software tool (Multi-Objective Primer Optimization for 16S experiments) under the GNU General Public License and is available at http://sysbiobig.dei.unipd.it/?q=Software#mopo16S . CONCLUSIONS: Results show that our strategy is able to find better primer pairs than the ones available in the literature according to all three optimization criteria. We also experimentally validated three of the primer pairs identified by our method on multiple bacterial species, belonging to different genera and phyla. Results confirm the predicted efficiency and the ability to maximize the number of different bacterial 16S sequences matched by primers.

Pathway Inspector: a pathway based web application for RNAseq analysis of model and non-model organisms.

SUMMARY: Pathway Inspector is an easy-to-use web application helping researchers to find patterns of expression in complex RNAseq experiments. The tool combines two standard approaches for RNAseq analysis: the identification of differentially expressed genes and a topology-based analysis of enriched pathways. Pathway Inspector is equipped with ad hoc interactive graphical interfaces simplifying the discovery of modulated pathways and the integration of the differentially expressed genes in the corresponding pathway topology. AVAILABILITY AND IMPLEMENTATION: Pathway Inspector is available at the website http://admiral.fmach.it/PI and has been developed in Python, making use of the Django Web Framework. CONTACT: Contact:paolo.fontana@fmach.it

Enhancing protein function prediction with taxonomic constraints--The Argot2.5 web server.

Argot2.5 (Annotation Retrieval of Gene Ontology Terms) is a web server designed to predict protein function. It is an updated version of the previous Argot2 enriched with new features in order to enhance its usability and its overall performance. The algorithmic strategy exploits the grouping of Gene Ontology terms by means of semantic similarity to infer protein function. The tool has been challenged over two independent benchmarks and compared to Argot2, PANNZER, and a baseline method relying on BLAST, proving to obtain a better performance thanks to the contribution of some key interventions in critical steps of the working pipeline. The most effective changes regard: (a) the selection of the input data from sequence similarity searches performed against a clustered version of UniProt databank and a remodeling of the weights given to Pfam hits, (b) the application of taxonomic constraints to filter out annotations that cannot be applied to proteins belonging to the species under investigation. The taxonomic rules are derived from our in-house developed tool, FunTaxIS, that extends those provided by the Gene Ontology consortium. The web server is free for academic users and is available online at http://www.medcomp.medicina.unipd.it/Argot2-5/.

Reprogramming Methods Do Not Affect Gene Expression Profile of Human Induced Pluripotent Stem Cells.

Induced pluripotent stem cells (iPSCs) are pluripotent cells derived from adult somatic cells. After the pioneering work by Yamanaka, who first generated iPSCs by retroviral transduction of four reprogramming factors, several alternative methods to obtain iPSCs have been developed in order to increase the yield and safety of the process. However, the question remains open on whether the different reprogramming methods can influence the pluripotency features of the derived lines. In this study, three different strategies, based on retroviral vectors, episomal vectors, and Sendai virus vectors, were applied to derive iPSCs from human fibroblasts. The reprogramming efficiency of the methods based on episomal and Sendai virus vectors was higher than that of the retroviral vector-based approach. All human iPSC clones derived with the different methods showed the typical features of pluripotent stem cells, including the expression of alkaline phosphatase and stemness maker genes, and could give rise to the three germ layer derivatives upon embryoid bodies assay. Microarray analysis confirmed the presence of typical stem cell gene expression profiles in all iPSC clones and did not identify any significant difference among reprogramming methods. In conclusion, the use of different reprogramming methods is equivalent and does not affect gene expression profile of the derived human iPSCs.

Isolation of infectious Zika virus from saliva and prolonged viral RNA shedding in a traveller returning from the Dominican Republic to Italy, January 2016.

We report the isolation of infectious Zika virus (ZIKV) in cell culture from the saliva of a patient who developed a febrile illness after returning from the Dominican Republic to Italy, in January 2016. The patient had prolonged shedding of viral RNA in saliva and urine, at higher load than in blood, for up to 29 days after symptom onset. Sequencing of ZIKV genome showed relatedness with strains from Latin America.

Infection dynamics in a traveller with persistent shedding of Zika virus RNA in semen for six months after returning from Haiti to Italy, January 2016.

We describe the dynamics of Zika virus (ZIKV) infection in a man in his early 40s who developed fever and rash after returning from Haiti to Italy, in January 2016. Follow-up laboratory testing demonstrated detectable ZIKV RNA in plasma up to day 9 after symptom onset and in urine and saliva up to days 15 and 47, respectively. Notably, persistent shedding of ZIKV RNA was demonstrated in semen, still detectable at 181 days after onset.

keeSeek: searching distant non-existing words in genomes for PCR-based applications.

UNLABELLED: The search for short words that are absent in the genome of one or more organisms (neverwords, also known as nullomers) is attracting growing interest because of the impact they may have in recent molecular biology applications. keeSeek is able to find absent sequences with primer-like features, which can be used as unique labels for exogenously inserted DNA fragments to recover their exact position into the genome using PCR techniques. The main differences with respect to previously developed tools for neverwords generation are (i) calculation of the distance from the reference genome, in terms of number of mismatches, and selection of the most distant sequences that will have a low probability to anneal unspecifically; (ii) application of a series of filters to discard candidates not suitable to be used as PCR primers. KeeSeek has been implemented in C++ and CUDA (Compute Unified Device Architecture) to work in a General-Purpose Computing on Graphics Processing Units (GPGPU) environment. AVAILABILITY AND IMPLEMENTATION: Freely available under the Q Public License at http://www.medcomp.medicina.unipd.it/main_site/doku.php?id=keeseek.

Changes in microRNA expression during disease progression in patients with chronic viral hepatitis.

BACKGROUND & AIMS: MicroRNAs (miRNAs) have been involved in hepatocarcinogenesis, but little is known on their role in the progression of chronic viral hepatitis. Aim of this study was to identify miRNA signatures associated with stages of disease progression in patients with chronic viral hepatitis. METHODS: MiRNA expression profile was investigated in liver biopsies from patients with chronic viral hepatitis and correlated with clinical, virological and histopathological features. Relevant miRNAs were further investigated. RESULTS: Most of the significant changes in miRNA expression were associated with liver fibrosis stages and included the significant up-regulation of a group of miRNAs that were demonstrated to target the master regulators of epithelial-mesenchymal transition ZEB1 and ZEB2 and involved in the preservation of epithelial cell differentiation, but also in cell proliferation and fibrogenesis. In agreement with miRNA data, immunostaining of liver biopsies showed that expression of the epithelial marker E-cadherin was maintained in severe fibrosis/cirrhosis while expression of ZEBs and other markers of epithelial-mesenchymal transition were low or absent. Severe liver fibrosis was also significantly associated with the down-regulation of miRNAs with antiproliferative and tumour suppressor activity. Similar changes in miRNA and target gene expression were demonstrated along with disease progression in a mouse model of carbon tetrachloride (CCl4)-induced liver fibrosis, suggesting that they might represent a general response to liver injury. CONCLUSION: Chronic viral hepatitis progression is associated with the activation of miRNA pathways that promote cell proliferation and fibrogenesis, but preserve the differentiated hepatocyte phenotype.

Reducing bias in RNA sequencing data: a novel approach to compute counts.

BACKGROUND: In the last decade, Next-Generation Sequencing technologies have been extensively applied to quantitative transcriptomics, making RNA sequencing a valuable alternative to microarrays for measuring and comparing gene transcription levels. Although several methods have been proposed to provide an unbiased estimate of transcript abundances through data normalization, all of them are based on an initial count of the total number of reads mapping on each transcript. This procedure, in principle robust to random noise, is actually error-prone if reads are not uniformly distributed along sequences, as happens indeed due to sequencing errors and ambiguity in read mapping. Here we propose a new approach, called maxcounts, to quantify the expression assigned to an exon as the maximum of its per-base counts, and we assess its performance in comparison with the standard approach described above, which considers the total number of reads aligned to an exon. The two measures are compared using multiple data sets and considering several evaluation criteria: independence from gene-specific covariates, such as exon length and GC-content, accuracy and precision in the quantification of true concentrations and robustness of measurements to variations of alignments quality. RESULTS: Both measures show high accuracy and low dependency on GC-content. However, maxcounts expression quantification is less biased towards long exons with respect to the standard approach. Moreover, it shows lower technical variability at low expressions and is more robust to variations in the quality of alignments. CONCLUSIONS: In summary, we confirm that counts computed with the standard approach depend on the length of the feature they are summarized on, and are sensitive to the non-uniform distribution of reads along transcripts. On the opposite, maxcounts are robust to biases due to the non-uniformity distribution of reads and are characterized by a lower technical variability. Hence, we propose maxcounts as an alternative approach for quantitative RNA-sequencing applications.

Comparative analysis of algorithms for whole-genome assembly of pyrosequencing data.

Next-generation sequencing technologies have fostered an unprecedented proliferation of high-throughput sequencing projects and a concomitant development of novel algorithms for the assembly of short reads. In this context, an important issue is the need of a careful assessment of the accuracy of the assembly process. Here, we review the efficiency of a panel of assemblers, specifically designed to handle data from GS FLX 454 platform, on three bacterial data sets with different characteristics in terms of reads coverage and repeats content. Our aim is to investigate their strengths and weaknesses in the reconstruction of the reference genomes. In our benchmarking, we assess assemblers' performance, quantifying and characterizing assembly gaps and errors, and evaluating their ability to solve complex genomic regions containing repeats. The final goal of this analysis is to highlight pros and cons of each method, in order to provide the final user with general criteria for the right choice of the appropriate assembly strategy, depending on the specific needs. A further aspect we have explored is the relationship between coverage of a sequencing project and quality of the obtained results. The final outcome suggests that, for a good tradeoff between costs and results, the planned genome coverage of an experiment should not exceed 20-30 ×.

Immune System Deficiencies Do Not Alter SARS-CoV-2 Evolutionary Rate but Favour the Emergence of Mutations by Extending Viral Persistence.

During the COVID-19 pandemic, immunosuppressed patients showed prolonged SARS-CoV-2 infections, with several studies reporting the accumulation of mutations in the viral genome. The weakened immune system present in these individuals, along with the effect of antiviral therapies, are thought to create a favourable environment for intra-host viral evolution and have been linked to the emergence of new viral variants which strongly challenged containment measures and some therapeutic treatments. To assess whether impaired immunity could lead to the increased instability of viral genomes, longitudinal nasopharyngeal swabs were collected from eight immunocompromised patients and fourteen non-immunocompromised subjects, all undergoing SARS-CoV-2 infection. Intra-host viral evolution was compared between the two groups through deep sequencing, exploiting a probe-based enrichment method to minimise the possibility of artefactual mutations commonly generated in amplicon-based methods, which heavily rely on PCR amplification. Although, as expected, immunocompromised patients experienced significantly longer infections, the acquisition of novel intra-host viral mutations was similar between the two groups. Moreover, a thorough analysis of viral quasispecies showed that the variability of viral populations in the two groups is comparable not only at the consensus level, but also when considering low-frequency mutations. This study suggests that a compromised immune system alone does not affect SARS-CoV-2 within-host genomic variability.

A functional and structural comparative analysis of large tumor antigens reveals evolution of different importin α-dependent nuclear localization signals.

Nucleocytoplasmic transport regulates the passage of proteins between the nucleus and cytoplasm. In the best characterized pathway, importin (IMP) α bridges cargoes bearing basic, classical nuclear localization signals (cNLSs) to IMPß1, which mediates transport through the nuclear pore complex. IMPα recognizes three types of cNLSs via two binding sites: the major binding site accommodates monopartite cNLSs, the minor binding site recognizes atypical cNLSs, while bipartite cNLSs simultaneously interact with both major and minor sites. Despite the growing knowledge regarding IMPα-cNLS interactions, our understanding of the evolution of cNLSs is limited. We combined bioinformatic, biochemical, functional, and structural approaches to study this phenomenon, using polyomaviruses (PyVs) large tumor antigens (LTAs) as a model. We characterized functional cNLSs from all human (H)PyV LTAs, located between the LXCXE motif and origin binding domain. Surprisingly, the prototypical SV40 monopartite NLS is not well conserved; HPyV LTA NLSs are extremely heterogenous in terms of structural organization, IMPα isoform binding, and nuclear targeting abilities, thus influencing the nuclear accumulation properties of full-length proteins. While several LTAs possess bipartite cNLSs, merkel cell PyV contains a hybrid bipartite cNLS whose upstream stretch of basic amino acids can function as an atypical cNLS, specifically binding to the IMPα minor site upon deletion of the downstream amino acids after viral integration in the host genome. Therefore, duplication of a monopartite cNLS and subsequent accumulation of point mutations, optimizing interaction with distinct IMPα binding sites, led to the evolution of bipartite and atypical NLSs binding at the minor site.