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.

Prevalence and virulence potential of Aeromonas spp. isolated from human diarrheal samples in North East Italy

IMPORTANCE

In this work, we demonstrate the epidemiologic relevance of the Aeromonas genus as the cause of infective diarrhea in North East Italy, both in children and adult subjects, with the significative presence of highly pathogenic strains. Aeromonas strains possess a heterogeneous armamentarium of pathogenicity factors that allows the microbe to affect a wide range of human intestinal epithelial cell processes that justify the ability to induce diarrhea through different mechanisms and cause diseases of variable severity, as observed for other gastrointestinal pathogens. However, it remains to be determined whether specific genotype(s) are associated with clinical pictures of different severity to implement the diagnostic and therapeutic approaches for this relevant enteric pathogen.

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.

Role of the Extracytoplasmic Function Sigma Factor SigE in the Stringent Response of Mycobacterium tuberculosis

Bacteria respond to nutrient starvation implementing the stringent response, a stress signaling system resulting in metabolic remodeling leading to decreased growth rate and energy requirements. A well-characterized model of stringent response in Mycobacterium tuberculosis is the one induced by growth in low phosphate. The extracytoplasmic function (ECF) sigma factor SigE was previously suggested as having a key role in the activation of stringent response. In this study, we challenge this hypothesis by analyzing the temporal dynamics of the transcriptional response of a sigE mutant and its wild-type parental strain to low phosphate using RNA sequencing. We found that both strains responded to low phosphate with a typical stringent response trait, including the downregulation of genes encoding ribosomal proteins and RNA polymerase. We also observed transcriptional changes that support the occurring of an energetics imbalance, compensated by a reduced activity of the electron transport chain, decreased export of protons, and a remodeling of central metabolism. The most striking difference between the two strains was the induction in the sigE mutant of several stress-related genes, in particular, the genes encoding the ECF sigma factor SigH and the transcriptional regulator WhiB6. Since both proteins respond to redox unbalances, their induction suggests that the sigE mutant is not able to maintain redox homeostasis in response to the energetics imbalance induced by low phosphate. In conclusion, our data suggest that SigE is not directly involved in initiating stringent response but in protecting the cell from stress consequent to the low phosphate exposure and activation of stringent response. IMPORTANCE Mycobacterium tuberculosis can enter a dormant state enabling it to establish latent infections and to become tolerant to antibacterial drugs. Dormant bacteria's physiology and the mechanism(s) used by bacteria to enter dormancy during infection are still unknown due to the lack of reliable animal models. However, several in vitro models, mimicking conditions encountered during infection, can reproduce different aspects of dormancy (growth arrest, metabolic slowdown, drug tolerance). The stringent response, a stress response program enabling bacteria to cope with nutrient starvation, is one of them. In this study, we provide evidence suggesting that the sigma factor SigE is not directly involved in the activation of stringent response as previously hypothesized, but it is important to help the bacteria to handle the metabolic stress related to the adaptation to low phosphate and activation of stringent response, thus giving an important contribution to our understanding of the mechanism behind stringent response development.

Impact of antigen test target failure and testing strategies on the transmission of SARS-CoV-2 variants

Population testing remains central to COVID-19 control and surveillance, with countries increasingly using antigen tests rather than molecular tests. Here we describe a SARS-CoV-2 variant that escapes N antigen tests due to multiple disruptive amino-acid substitutions in the N protein. By fitting a multistrain compartmental model to genomic and epidemiological data, we show that widespread antigen testing in the Italian region of Veneto favored the undetected spread of the antigen-escape variant compared to the rest of Italy. We highlight novel limitations of widespread antigen testing in the absence of molecular testing for diagnostic or confirmatory purposes. Notably, we find that genomic surveillance systems which rely on antigen population testing to identify samples for sequencing will bias detection of escape antigen test variants. Together, these findings highlight the importance of retaining molecular testing for surveillance purposes, including in contexts where the use of antigen tests is widespread.

Human Virus Genomes Are Enriched in Conserved Adenine/Thymine/Uracil Multiple Tracts That Pause Polymerase Progression

The DNA secondary structures that deviate from the classic Watson and Crick base pairing are increasingly being reported to form transiently in the cell and regulate specific cellular mechanisms. Human viruses are cell parasites that have evolved mechanisms shared with the host cell to support their own replication and spreading. Contrary to human host cells, viruses display a diverse array of nucleic acid types, which include DNA or RNA in single-stranded or double-stranded conformations. This heterogeneity improves the possible occurrence of non-canonical nucleic acid structures. We have previously shown that human virus genomes are enriched in G-rich sequences that fold in four-stranded nucleic acid secondary structures, the G-quadruplexes.Here, by extensive bioinformatics analysis on all available genomes, we showed that human viruses are enriched in highly conserved multiple A (and T or U) tracts, with such an array that they could in principle form quadruplex structures. By circular dichroism, NMR, and Taq polymerase stop assays, we proved that, while A/T/U-quadruplexes do not form, these tracts still display biological significance, as they invariably trigger polymerase pausing within two bases from the A/T/U tract. “A” bases display the strongest effect. Most of the identified A-tracts are in the coding strand, both at the DNA and RNA levels, suggesting their possible relevance during viral translation. This study expands on the presence and mechanism of nucleic acid secondary structures in human viruses and provides a new direction for antiviral research.

Longitudinal analysis of T cell receptor repertoires reveals shared patterns of antigen-specific response to SARS-CoV-2 infection

T cells play a prominent role in orchestrating the immune response to viral diseases, but their role in the clinical presentation and subsequent immunity to SARS-CoV-2 infection remains poorly understood. As part of a population-based survey of the municipality of Vo', Italy, conducted after the initial SARS-CoV-2 outbreak, we sampled the T cell receptor (TCR) repertoires of the population 2 months after the initial PCR survey and followed up positive cases 9 and 15 months later. At 2 months, we found that 97.0% (98 of 101) of cases had elevated levels of TCRs associated with SARS-CoV-2. T cell frequency (depth) was increased in individuals with more severe disease. Both depth and diversity (breadth) of the TCR repertoire were positively associated with neutralizing antibody titers, driven mostly by CD4+ T cells directed against spike protein. At the later time points, detection of these TCRs remained high, with 90.7% (78 of 96) and 86.2% (25 of 29) of individuals having detectable signal at 9 and 15 months, respectively. Forty-three individuals were vaccinated by month 15 and showed a significant increase in TCRs directed against spike protein. Taken together, these results demonstrate the central role of T cells in mounting an immune defense against SARS-CoV-2 that persists out to 15 months.

Prospective epidemiological, molecular, and genetic characterization of a novel coronavirus disease in the Val Venosta/Vinschgau: the CHRIS COVID-19 study protocol

The COVID-19 pandemic has been threatening the healthcare and socioeconomic systems of entire nations. While population-based surveys to assess the distribution of SARS-CoV-2 infection have become a priority, pre-existing longitudinal studies are ideally suited to assess the determinants of COVID-19 onset and severity.The Cooperative Health Research In South Tyrol (CHRIS) study completed the baseline recruitment of 13,393 adults from the Venosta/Vinschgau rural district in 2018, collecting extensive phenotypic and biomarker data, metabolomic data, densely imputed genotype and whole-exome sequencing data.Based on CHRIS, we designed a prospective study, called CHRIS COVID-19, aimed at: 1) estimating the incidence of SARS-CoV-2 infections; 2) screening for and investigating the determinants of incident infection among CHRIS participants and their household members; 3) monitoring the immune response of infected participants prospectively.An online screening questionnaire was sent to all CHRIS participants and their household members. A random sample of 1450 participants representative of the district population was invited to assess active (nasopharyngeal swab) or past (serum antibody test) infections. We prospectively invited for complete SARS-CoV-2 testing all questionnaire completers gauged as possible cases of past infection and their household members. In positive tested individuals, antibody response is monitored quarterly for one year. Untested and negative participants receive the screening questionnaire every four weeks until gauged as possible incident cases or till the study end.Originated from a collaboration between researchers and community stakeholders, the CHRIS COVID-19 study aims at generating knowledge about the epidemiological, molecular, and genetic characterization of COVID-19 and its long-term sequelae.

Disease Severity and Prognosis of SARS-CoV-2 Infection in Hospitalized Patients Is Not Associated With Viral Load in Nasopharyngeal Swab

Background: The impact of viral burden on severity and prognosis of patients hospitalized for Coronavirus Disease 2019 (COVID-19) is still a matter of debate due to controversial results. Herein, we sought to assess viral load in the nasopharyngeal swab and its association with severity score indexes and prognostic parameters. Methods: We included 127 symptomatic patients and 21 asymptomatic subjects with a diagnosis of SARS-CoV-2 infection obtained by reverse transcription polymerase chain reaction and presence of cycle threshold. According to the level of care needed during hospitalization, the population was categorized as high-intensity (HIMC, n = 76) or low intensity medical care setting (LIMC, n = 51). Results: Viral load did not differ among asymptomatic, LIMC, and HIMC SARS-CoV-2 positive patients [4.4 (2.9-5.3) vs. 4.8 (3.6-6.1) vs. 4.6 (3.9-5.7) log10 copies/ml, respectively; p = 0.31]. Similar results were observed when asymptomatic individuals were compared to hospitalized patients [4.4 (2.9-5.3) vs. 4.68 (3.8-5.9) log10 copies/ml; p = 0.13]. When the study population was divided in High (HVL, n = 64) and Low Viral Load (LVL, n = 63) group no differences were observed in disease severity at diagnosis. Furthermore, LVL and HVL groups did not differ with regard to duration of hospital stay, number of bacterial co-infections, need for high-intensity medical care and number of deaths. The viral load was not an independent risk factor for HIMC in an adjusted multivariate regression model (OR: 1.59; 95% CI: 0.46-5.55, p = 0.46). Conclusions: Viral load at diagnosis is similar in asymptomatic and hospitalized patients and is not associated with either worse outcomes during hospitalization. SARS CoV-2 viral load might not be the right tool to assist clinicians in risk-stratifying hospitalized patients.

SARS-CoV-2 antibody dynamics and transmission from community-wide serological testing in the Italian municipality of Vo'

In February and March 2020, two mass swab testing campaigns were conducted in Vo', Italy. In May 2020, we tested 86% of the Vo' population with three immuno-assays detecting antibodies against the spike and nucleocapsid antigens, a neutralisation assay and Polymerase Chain Reaction (PCR). Subjects testing positive to PCR in February/March or a serological assay in May were tested again in November. Here we report on the results of the analysis of the May and November surveys. We estimate a seroprevalence of 3.5% (95% Credible Interval (CrI): 2.8-4.3%) in May. In November, 98.8% (95% Confidence Interval (CI): 93.7-100.0%) of sera which tested positive in May still reacted against at least one antigen; 18.6% (95% CI: 11.0-28.5%) showed an increase of antibody or neutralisation reactivity from May. Analysis of the serostatus of the members of 1,118 households indicates a 26.0% (95% CrI: 17.2-36.9%) Susceptible-Infectious Transmission Probability. Contact tracing had limited impact on epidemic suppression.

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.

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) infection¹. 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 province². In response, the regional authorities imposed the lockdown of the whole municipality for 14 days³. 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.

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) infection¹. 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 province². In response, the regional authorities imposed the lockdown of the whole municipality for 14 days³. 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.

Human West Nile Virus Lineage 2 Infection: Epidemiological, Clinical, and Virological Findings

West Nile virus (WNV) lineage 2 is expanding and causing large outbreaks in Europe. In this study, we analyzed the epidemiological, clinical, and virological features of WNV lineage 2 infection during the large outbreak that occurred in northern Italy in 2018. The study population included 86 patients with neuroinvasive disease (WNND), 307 with fever (WNF), and 34 blood donors. Phylogenetic analysis of WNV full genome sequences from patients’ samples showed that the virus belonged to the widespread central/southern European clade of WNV lineage 2 and was circulating in the area at least since 2014. The incidence of WNND and WNF progressively increased with age and was higher in males than in females. Among WNND patients, the case fatality rate was 22%. About 70% of blood donors reported symptoms during follow-up. Within the first week after symptom onset, WNV RNA was detectable in the blood or urine of 80% of patients, while 20% and 40% of WNND and WNF patients, respectively, were WNV IgM-seronegative. In CSF samples of WNND patients, WNV RNA was typically detectable when WNV IgM antibodies were absent. Blunted or no WNV IgM response and high WNV IgG levels were observed in seven patients with previous flavivirus immunity.

Measles Virus Infection and Immunity in a Suboptimal Vaccination Coverage Setting

Despite efforts to improve surveillance and vaccination coverage, measles virus (MeV) continues to cause outbreaks also in high-income countries. As the reference laboratory of the Veneto Region, Italy, we analyzed changes in population immunity, described measles outbreaks, investigated MeV genetic diversity, and evaluated cross-protection of measles vaccination against MeV epidemic strains. Like most European areas, the Veneto Region has suboptimal measles vaccination coverage and is facing a growing public mistrust of vaccination. A progressive decline of measles vaccine uptake was observed during the last decade in the Veneto Region, leading to immunity gaps in children and young adults. Measles outbreaks were caused by the same MeV genotype B3, D4, and D8 strains that were circulating in other European countries. Eleven cases of measles were observed in immunized subjects. These cases were not associated with particular MeV genotypes nor with mutations in epitopes recognized by neutralizing antibodies. Accordingly, sera from fully vaccinated subjects cross-neutralized epidemic MeV strains, including the genotypes B3, D4, and D8, with the same high efficiency demonstrated against the vaccine strain. In fully vaccinated subjects, high MeV IgG antibody titers persisted up to 30 years following vaccination. These results support the use of the current measles-containing vaccines and strategies to strengthen vaccination.

The CAFA challenge reports improved protein function prediction and new functional annotations for hundreds of genes through experimental screens

BACKGROUND

The Critical Assessment of Functional Annotation (CAFA) is an ongoing, global, community-driven effort to evaluate and improve the computational annotation of protein function.

RESULTS

Here, we report on the results of the third CAFA challenge, CAFA3, that featured an expanded analysis over the previous CAFA rounds, both in terms of volume of data analyzed and the types of analysis performed. In a novel and major new development, computational predictions and assessment goals drove some of the experimental assays, resulting in new functional annotations for more than 1000 genes. Specifically, we performed experimental whole-genome mutation screening in Candida albicans and Pseudomonas aureginosa genomes, which provided us with genome-wide experimental data for genes associated with biofilm formation and motility. We further performed targeted assays on selected genes in Drosophila melanogaster, which we suspected of being involved in long-term memory.

CONCLUSION

We conclude that while predictions of the molecular function and biological process annotations have slightly improved over time, those of the cellular component have not. Term-centric prediction of experimental annotations remains equally challenging; although the performance of the top methods is significantly better than the expectations set by baseline methods in C. albicans and D. melanogaster, it leaves considerable room and need for improvement. Finally, we report that the CAFA community now involves a broad range of participants with expertise in bioinformatics, biological experimentation, biocuration, and bio-ontologies, working together to improve functional annotation, computational function prediction, and our ability to manage big data in the era of large experimental screens.

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 (PrM₃, E₁₅₉ and E₂₃₁) and in non-structural (NS1₁₀₉, NS5₂₅₉, NS5₃₁₀ and NS5₆₀₀) proteins. Observed amino acid substitutions Phe₃ and Ser₂₃₁ 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.

Transcriptome and Cell Physiological Analyses in Different Rice Cultivars Provide New Insights Into Adaptive and Salinity Stress Responses

Salinity tolerance has been extensively investigated in recent years due to its agricultural importance. Several features, such as the regulation of ionic transporters and metabolic adjustments, have been identified as salt tolerance hallmarks. Nevertheless, due to the complexity of the trait, the results achieved to date have met with limited success in improving the salt tolerance of rice plants when tested in the field, thus suggesting that a better understanding of the tolerance mechanisms is still required. In this work, differences between two varieties of rice with contrasting salt sensitivities were revealed by the imaging of photosynthetic parameters, ion content analysis and a transcriptomic approach. The transcriptomic analysis conducted on tolerant plants supported the setting up of an adaptive program consisting of sodium distribution preferentially limited to the roots and older leaves, and in the activation of regulatory mechanisms of photosynthesis in the new leaves. As a result, plants resumed grow even under prolonged saline stress. In contrast, in the sensitive variety, RNA-seq analysis revealed a misleading response, ending in senescence and cell death. The physiological response at the cellular level was investigated by measuring the intracellular profile of H₂O₂ in the roots, using a fluorescent probe. In the roots of tolerant plants, a quick response was observed with an increase in H₂O₂ production within 5 min after salt treatment. The expression analysis of some of the genes involved in perception, signal transduction and salt stress response confirmed their early induction in the roots of tolerant plants compared to sensitive ones. By inhibiting the synthesis of apoplastic H₂O₂, a reduction in the expression of these genes was detected. Our results indicate that quick H₂O₂ signaling in the roots is part of a coordinated response that leads to adaptation instead of senescence in salt-treated rice plants.

Mapping and characterization of G-quadruplexes in Mycobacterium tuberculosis gene promoter regions

Mycobacterium tuberculosis is the causative agent of tuberculosis (TB), one of the top 10 causes of death worldwide in 2015. The recent emergence of strains resistant to all current drugs urges the development of compounds with new mechanisms of action. G-quadruplexes are nucleic acids secondary structures that may form in G-rich regions to epigenetically regulate cellular functions. Here we implemented a computational tool to scan the presence of putative G-quadruplex forming sequences in the genome of Mycobacterium tuberculosis and analyse their association to transcription start sites. We found that the most stable G-quadruplexes were in the promoter region of genes belonging to definite functional categories. Actual G-quadruplex folding of four selected sequences was assessed by biophysical and biomolecular techniques: all molecules formed stable G-quadruplexes, which were further stabilized by two G-quadruplex ligands. These compounds inhibited Mycobacterium tuberculosis growth with minimal inhibitory concentrations in the low micromolar range. These data support formation of Mycobacterium tuberculosis G-quadruplexes in vivo and their potential regulation of gene transcription, and prompt the use of G4 ligands to develop original antitubercular agents.

Reconstructing the recent West Nile virus lineage 2 epidemic in Europe and Italy using discrete and continuous phylogeography

West Nile virus lineage 2 (WNV-2) was mainly confined to sub-Saharan Africa until the early 2000s, when it was identified for the first time in Central Europe causing outbreaks of human and animal infection. The aim of this study was to reconstruct the origin and dispersion of WNV-2 in Central Europe and Italy on a phylodynamic and phylogeographical basis. To this aim, discrete and continuous space phylogeographical models were applied to a total of 33 newly characterised full-length viral genomes obtained from mosquitoes, birds and humans in Northern Italy in the years 2013–2015 aligned with 64 complete sequences isolated mainly in Europe. The European isolates segregated into two highly significant clades: a small one including three sequences and a large clade including the majority of isolates obtained in Central Europe since 2004. Discrete phylogeographical analysis showed that the most probable location of the root of the largest European clade was in Hungary a mean 12.78 years ago. The European clade bifurcated into two highly supported subclades: one including most of the Central/East European isolates and the other encompassing all of the isolates obtained in Greece. The continuous space phylogeographical analysis of the Italian clade showed that WNV-2 entered Italy in about 2008, probably by crossing the Adriatic sea and reaching a central area of the Po Valley. The epidemic then spread simultaneously eastward, to reach the region of the Po delta in 2013, and westward to the border area between Lombardy and Piedmont in 2014; later, the western strain changed direction southward, and reached the central area of the Po valley once again in 2015. Over a period of about seven years, the virus spread all over an area of northern Italy by following the Po river and its main tributaries.

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.

An expanded evaluation of protein function prediction methods shows an improvement in accuracy

BACKGROUND

A major bottleneck in our understanding of the molecular underpinnings of life is the assignment of function to proteins. While molecular experiments provide the most reliable annotation of proteins, their relatively low throughput and restricted purview have led to an increasing role for computational function prediction. However, assessing methods for protein function prediction and tracking progress in the field remain challenging.

RESULTS

We conducted the second critical assessment of functional annotation (CAFA), a timed challenge to assess computational methods that automatically assign protein function. We evaluated 126 methods from 56 research groups for their ability to predict biological functions using Gene Ontology and gene-disease associations using Human Phenotype Ontology on a set of 3681 proteins from 18 species. CAFA2 featured expanded analysis compared with CAFA1, with regards to data set size, variety, and assessment metrics. To review progress in the field, the analysis compared the best methods from CAFA1 to those of CAFA2.

CONCLUSIONS

The top-performing methods in CAFA2 outperformed those from CAFA1. This increased accuracy can be attributed to a combination of the growing number of experimental annotations and improved methods for function prediction. The assessment also revealed that the definition of top-performing algorithms is ontology specific, that different performance metrics can be used to probe the nature of accurate predictions, and the relative diversity of predictions in the biological process and human phenotype ontologies. While there was methodological improvement between CAFA1 and CAFA2, the interpretation of results and usefulness of individual methods remain context-dependent.

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.

Zika virus: from pathogenesis to disease control

Zika virus is a mosquito-borne flavivirus discovered in Uganda in 1947. The virus has emerged in recent years and spread in the Pacific Area and the Americas, where it has caused large human outbreaks. The factors involved in the virus's emergence are still unknown, but probably include its introduction in naïve environments characterised by the presence of high densities of competent Aedes spp. mosquitoes and susceptible human hosts in urban areas. Unique features of Zika virus infection are sexual and transplacental transmission and associated neurological morbidities, i.e. Guillain-Barré syndrome and fetal microcephaly. Diagnosis relies on the detection of viral nucleic acids in biological samples, while detection of a specific antibody response may be inconclusive because of the broad cross-reactivity of antibodies among flaviviruses. Experimental studies have clarified some mechanisms of Zika virus pathogenesis and have identified potential targets for antiviral drugs. In animal models, the virus can infect and efficiently replicate in the placenta and in the brain, and induce fetal demise or neural damage, recapitulating human diseases. These animal models have been used to evaluate candidate vaccines and promising results have been obtained.

Eliciting the Functional Taxonomy from protein annotations and taxa

The advances of omics technologies have triggered the production of an enormous volume of data coming from thousands of species. Meanwhile, joint international efforts like the Gene Ontology (GO) consortium have worked to provide functional information for a vast amount of proteins. With these data available, we have developed FunTaxIS, a tool that is the first attempt to infer functional taxonomy (i.e. how functions are distributed over taxa) combining functional and taxonomic information. FunTaxIS is able to define a taxon specific functional space by exploiting annotation frequencies in order to establish if a function can or cannot be used to annotate a certain species. The tool generates constraints between GO terms and taxa and then propagates these relations over the taxonomic tree and the GO graph. Since these constraints nearly cover the whole taxonomy, it is possible to obtain the mapping of a function over the taxonomy. FunTaxIS can be used to make functional comparative analyses among taxa, to detect improper associations between taxa and functions, and to discover how functional knowledge is either distributed or missing. A benchmark test set based on six different model species has been devised to get useful insights on the generated taxonomic rules.

Third generation sequencing technologies applied to diagnostic microbiology: benefits and challenges in applications and data analysis

INTRODUCTION

The diagnosis of infectious diseases is among the most successful areas of application of new generation sequencing technologies. The field has seen the development of numerous experimental and analytical approaches for the detection and the fine description of pathogenic and non-pathogenic microorganisms.

AREAS COVERED

Without claiming to be exhaustive with respect to all applications and methods developed over the years, this review focuses on the advantages and the issues brought by the new technologies, with an eye in particular to third generation sequencing methods. Both experimental procedures and algorithmic strategies are presented, following the most relevant publications which have led to progress in our ability of detecting infectious agents. Expert commentary: The technical advance brought by third generation sequencing platforms has the potential to significantly expand the range of diagnostic tools that will be available to clinicians. Nonetheless, the implementation of these technologies in clinical practice is still far from being actionable and will temporally follow the path undertaken by second generation methods, which still require the setup of standardized pipelines in both wet and dry laboratory procedures.

Full Genome Sequence-Based Comparative Study of Wild-Type and Vaccine Strains of Infectious Laryngotracheitis Virus from Italy

Infectious laryngotracheitis (ILT) is an acute and highly contagious respiratory disease of chickens caused by an alphaherpesvirus, infectious laryngotracheitis virus (ILTV). Recently, full genome sequences of wild-type and vaccine strains have been determined worldwide, but none was from Europe. The aim of this study was to determine and analyse the complete genome sequences of five ILTV strains. Sequences were also compared to reveal the similarity of strains across time and to discriminate between wild-type and vaccine strains. Genomes of three ILTV field isolates from outbreaks occurred in Italy in 1980, 2007 and 2011, and two commercial chicken embryo origin (CEO) vaccines were sequenced using the 454 Life Sciences technology. The comparison with the Serva genome showed that 35 open reading frames (ORFs) differed across the five genomes. Overall, 54 single nucleotide polymorphisms (SNPs) and 27 amino acid differences in 19 ORFs and two insertions in the UL52 and ORFC genes were identified. Similarity among the field strains and between the field and the vaccine strains ranged from 99.96% to 99.99%. Phylogenetic analysis revealed a close relationship among them, as well. This study generated data on genomic variation among Italian ILTV strains revealing that, even though the genetic variability of the genome is well conserved across time and between wild-type and vaccine strains, some mutations may help in differentiating among them and may be involved in ILTV virulence/attenuation. The results of this study can contribute to the understanding of the molecular bases of ILTV pathogenicity and provide genetic markers to differentiate between wild-type and vaccine strains.

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/.

Phylogenetic characterization of Central/Southern European lineage 2 West Nile virus: analysis of human outbreaks in Italy and Greece, 2013-2014

In recent years, West Nile virus (WNV) lineage 2 has been spreading and causing disease outbreaks in humans and animals in Europe. In order to characterize viral diversity, we performed full-length genome sequencing of WNV lineage 2 from human samples collected during outbreaks in Italy and Greece in 2013 and 2014. Phylogenetic analysis showed that these WNV lineage 2 genomes belonged to a monophyletic clade derived from a single introduction into Europe of the prototype Hungarian strain. Correlation of phylogenetic data with geospatial information showed geographical clustering of WNV genome sequences both in Italy and in Greece, indicating that the virus had evolved and diverged during its dispersal in Europe, leading to the emergence of novel genotypes, as it adapted to local ecological niches. These genotypes carried divergent conserved amino acid substitutions, which might have been relevant for viral adaptation, as suggested by selection pressure analysis and in silico and experimental modelling of sequence changes. In conclusion, the results of this study provide further information on WNV lineage 2 transmission dynamics in Europe, and emphasize the need for WNV surveillance activities to monitor viral evolution and diversity.

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.

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

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.

Oxidative damage in the progression of chronic liver disease to hepatocellular carcinoma: An intricate pathway

The histo-pathologic and molecular mechanisms leading to initiation and progression of hepatocellular carcinoma (HCC) are still ill-defined; however, there is increasing evidence that the gradual accumulation of mutations, genetic and epigenetic changes which occur in preneoplastic hepatocytes results in the development of dysplastic foci, nodules, and finally, overt HCC. As well as many other neoplasias, liver cancer is considered an “inflammatory cancer”, arising from a context of inflammation, and characterized by inflammation-related mechanisms that favor tumor cell survival, proliferation, and invasion. Molecular mechanisms that link inflammation and neoplasia have been widely investigated, and it has been well established that inflammatory cells recruited at these sites with ongoing inflammatory activity release chemokines that enhance the production of reactive oxygen species. The latter, in turn, probably have a major pathogenic role in the continuum starting from hepatitis followed by chronic inflammation, and ultimately leading to cancer. The relationship amongst chronic liver injury, free radical production, and development of HCC is explored in the present review, particularly in the light of the complex network that involves oxidative DNA damage, cytokine synthesis, telomere dysfunction, and microRNA regulation.

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.

Large human outbreak of West Nile virus infection in north-eastern Italy in 2012

Human cases of West Nile virus (WNV) disease have been reported in Italy since 2008. So far, most cases have been identified in north-eastern Italy, where, in 2012, the largest outbreak of WNV infection ever recorded in Italy occurred. Most cases of the 2012 outbreak were identified in the Veneto region, where a special surveillance plan for West Nile fever was in place. In this outbreak, 25 cases of West Nile neuroinvasive disease and 17 cases of fever were confirmed. In addition, 14 WNV RNA-positive blood donors were identified by screening of blood and organ donations and two cases of asymptomatic infection were diagnosed by active surveillance of subjects at risk of WNV exposure. Two cases of death due to WNND were reported. Molecular testing demonstrated the presence of WNV lineage 1 in all WNV RNA-positive patients and, in 15 cases, infection by the novel Livenza strain was ascertained. Surveillance in other Italian regions notified one case of neuroinvasive disease in the south of Italy and two cases in Sardinia. Integrated surveillance for WNV infection remains a public health priority in Italy and vector control activities have been strengthened in areas of WNV circulation.

Genomic comparative analysis and gene function prediction in infectious diseases: application to the investigation of a meningitis outbreak

Background

Next generation sequencing (NGS) is being increasingly used for the detection and characterization of pathogens during outbreaks. This technology allows rapid sequencing of pathogen full genomes, useful not only for accurate genotyping and molecular epidemiology, but also for identification of drug resistance and virulence traits.

Methods

In this study, an approach based on whole genome sequencing by NGS, comparative genomics, and gene function prediction was set up and retrospectively applied for the investigation of two N. meningitidis serogroup C isolates collected from a cluster of meningococcal disease, characterized by a high fatality rate.

Results

According to conventional molecular typing methods, all the isolates had the same typing results and were classified as outbreak isolates within the same N. meningitidis sequence type ST-11, while full genome sequencing demonstrated subtle genetic differences between the isolates. Looking for these specific regions by means of 9 PCR and cycle sequencing assays in other 7 isolates allowed distinguishing outbreak cases from unrelated cases. Comparative genomics and gene function prediction analyses between outbreak isolates and a set of reference N. meningitidis genomes led to the identification of differences in gene content that could be relevant for pathogenesis. Most genetic changes occurred in the capsule locus and were consistent with recombination and horizontal acquisition of a set of genes involved in capsule biosynthesis.

Conclusions

This study showed the added value given by whole genome sequencing by NGS over conventional sequence-based typing methods in the investigation of an outbreak. Routine application of this technology in clinical microbiology will significantly improve methods for molecular epidemiology and surveillance of infectious disease and provide a bulk of data useful to improve our understanding of pathogens biology.

Genome sequencing of West Nile Virus from human cases in Greece, 2012

A West Nile Virus (WNV) lineage 2 strain, named Nea Santa-Greece-2010, has been demonstrated to be responsible for the large outbreaks of neuroinvasive disease (WNND) that have been occurring in Greece since 2010, based on sequence similarities of viral isolates identified between 2010–2012. However, knowledge on the evolution of this strain is scarce because only partial WNV genome sequences are available from Greece. The aim of this study was to get the complete genome sequence of WNV from patients with infection. To this aim, plasma and urine samples collected during the 2012 Greek outbreak were retrospectively investigated. Full WNV genome sequence was obtained from a patient with WNND. The genome had 99.7% sequence identity to Nea Santa, higher than to other related WNV lineage 2 strains, and five amino acid changes apparently not relevant for viral pathogenicity or fitness. In addition, infection by WNV lineage 2 was confirmed in additional nine patients with WNND; in three of them the infection with WNV Nea Santa was demonstrated by sequencing. In conclusion, this study characterized for the first time a WNV full genome from a patient with WNND from Greece, demonstrated the persistence of the Nea Santa strain, and suggested that the virus might have locally evolved.

Whole genome sequencing and phylogenetic analysis of West Nile virus lineage 1 and lineage 2 from human cases of infection, Italy, August 2013

A human outbreak of West Nile virus (WNV) infection caused by WNV lineage 2 is ongoing in northern Italy. Analysis of six WNV genome sequences obtained from clinical specimens demonstrated similarities with strains circulating in central Europe and Greece and the presence of unique amino acid changes that identify a new viral strain. In addition, WNV lineage 1 Livenza, responsible for a large outbreak in north-eastern Italy in 2012, was fully sequenced from a blood donor during this 2013 outbreak.

Next-generation sequencing technologies in diagnostic virology

The data deluge produced by next-generation sequencing (NGS) technologies is an appealing feature for clinical virologists that are involved in the diagnosis of emerging viral infections, molecular epidemiology of viral pathogens, drug-resistance testing, and also like to do some basic and clinical research. Indeed, NGS platforms are being implemented in many clinical and research laboratories, as the costs of these platforms are progressively decreasing. We provide here some suggestions for virologists who are planning to implement a NGS platform in their clinical laboratory and an overview on the potential applications of these technologies in diagnostic virology.

A large-scale evaluation of computational protein function prediction

A report on the results of the first large-scale community-based critical assessment of protein function annotation (CAFA) experiment.

Automated annotation of protein function is challenging. As the number of sequenced genomes rapidly grows, the overwhelming majority of protein products can only be annotated computationally. If computational predictions are to be relied upon, it is crucial that the accuracy of these methods be high. Here we report the results from the first large-scale community-based critical assessment of protein function annotation (CAFA) experiment. Fifty-four methods representing the state of the art for protein function prediction were evaluated on a target set of 866 proteins from 11 organisms. Two findings stand out: (i) today's best protein function prediction algorithms substantially outperform widely used first-generation methods, with large gains on all types of targets; and (ii) although the top methods perform well enough to guide experiments, there is considerable need for improvement of currently available tools.