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.

A white paper on Phospholipid Hydroperoxide Glutathione Peroxidase (GPx4) forty years later

The purification of a protein inhibiting lipid peroxidation led to the discovery of the selenoperoxidase GPx4 forty years ago. Thus, the evidence of the enzymatic activity was reached after identifying the biological effect and unambiguously defined the relationship between the biological function and the enzymatic activity. In the syllogism where GPx4 inhibits lipid peroxidation and its inhibition is lethal, cell death is operated by lipid peroxidation. Based on this rationale, this form of cell death emerged as regulated iron-enforced oxygen toxicity and was named ferroptosis in 2012. In the last decades, we learned that reduction of lipid hydroperoxides is indispensable and, in cooperation with prooxidant systems, controls the critical steady state of lipid peroxidation. This concept defined the GPx4 reaction as both the target for possible anti-cancer therapy and if insufficient, as cause of degenerative diseases. We know the reaction mechanism, but the details of the interaction at the membrane cytosol interface are still poorly defined. We know the gene structure, but the knowledge about expression control is still limited. The same holds true for post-transcriptional modifications. Reverse genetics indicate that GPx4 has a role in inflammation, immunity, and differentiation, but the observations emerging from these studies need a more specifically addressed biochemical evidence. Finally, the role of GPx4 in spermatogenesis disclosed an area unconnected to lipid peroxidation. In its mitochondrial and nuclear form, the peroxidase catalyzes the oxidation of protein thiols in two specific aspects of sperm maturation: stabilization of the mid-piece and chromatin compaction. Thus, although available evidence converges to the notion that GPx4 activity is vital due to the inhibition of lipid peroxidation, it is reasonable to foresee other unknown aspects of the GPx4 reaction to be disclosed.

The glutathione peroxidase family: Discoveries and mechanism

The discoveries leading to our present understanding of the glutathione peroxidases (GPxs) are recalled. The cytosolic GPx, now GPx1, was first described by Mills in 1957 and claimed to depend on selenium by Rotruck et al., in 1972. With the determination of a stoichiometry of one selenium per subunit, GPx1 was established as the first selenoenzyme of vertebrates. In the meantime, the GPxs have grown up to a huge family of enzymes that prevent free radical formation from hydroperoxides and, thus, are antioxidant enzymes, but they are also involved in regulatory processes or synthetic functions. The kinetic mechanism of the selenium-containing GPxs is unusual in neither showing a defined K_M nor any substrate saturation. More recently, the reaction mechanism has been investigated by the density functional theory and nuclear magnetic resonance of model compounds mimicking the reaction cycle. The resulting concept sees a selenolate oxidized to a selenenic acid. This very fast reaction results from a concerted dual attack on the hydroperoxide bond, a nucleophilic one by the selenolate and an electrophilic one by a proton that is unstably bound in the reaction center. Postulated intermediates have been identified either in the native enzymes or in model compounds.

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.

A Circadian Hygiene Education Initiative Covering the Pre-pandemic and Pandemic Period Resulted in Earlier Get-Up Times in Italian University Students: An Ecological Study

The aims of the present study were to obtain sleep quality and sleep timing information in a group of university students and to evaluate the effects of a circadian hygiene education initiative. All students of the University of Padova (approximately 64,000) were contacted by e-mail (major campaigns in October 2019 and October 2020) and directed to an ad hoc website for collection of demographics and sleep quality/timing information. Participants (n = 5,740) received one of two sets of circadian hygiene advice (“A regular life” or “Bright days and dark nights”). Every month, they were then asked how easy it had been to comply and provided with the advice again. At any even month from joining, they completed the sleep quality/timing questionnaires again. Information on academic performance was obtained post hoc, together with representative samples of lecture (n = 5,972) and examination (n = 1,800) timings, plus lecture attendances (n = 25,302). Fifty-two percent of students had poor sleep quality, and 82% showed signs of social jetlag. Those who joined in October 2020, after several months of lockdown and distance learning, had better sleep quality, less social jetlag, and later sleep habits. Over approximately a year, the “Bright days and dark nights” advice resulted in significantly earlier get-up times compared with the “A regular life” advice. Similarly, it also resulted in a trend toward earlier midsleep (i.e., the midpoint, expressed as clock time, between sleep onset and sleep offset) and toward a decrease in the latency between wake-up and get-up time, with no impact on sleep duration. Significant changes in most sleep quality and sleep timing variables (i.e., fewer night awakenings, less social jetlag, and delayed sleep timing during lock-down) were observed in both advice groups over approximately a year, mostly in association with pandemic-related events characterizing 2020. Early chronotype students had better academic performances compared with their later chronotype counterparts. In a multivariate model, sleep quality, chronotype and study subject (science and technology, health and medical, or social and humanities) were independent predictors of academic performance. Taken together, these results underlie the importance of designing circadian-friendly university timetables.

Phylogeography and genomic epidemiology of SARS-CoV-2 in Italy and Europe with newly characterized Italian genomes between February-June 2020

The aims of this study were to characterize new SARS-CoV-2 genomes sampled all over Italy and to reconstruct the origin and the evolutionary dynamics in Italy and Europe between February and June 2020. The cluster analysis showed only small clusters including < 80 Italian isolates, while most of the Italian strains were intermixed in the whole tree. Pure Italian clusters were observed mainly after the lockdown and distancing measures were adopted. Lineage B and B.1 spread between late January and early February 2020, from China to Veneto and Lombardy, respectively. Lineage B.1.1 (20B) most probably evolved within Italy and spread from central to south Italian regions, and to European countries. The lineage B.1.1.1 (20D) developed most probably in other European countries entering Italy only in the second half of March and remained localized in Piedmont until June 2020. In conclusion, within the limitations of phylogeographical reconstruction, the estimated ancestral scenario suggests an important role of China and Italy in the widespread diffusion of the D614G variant in Europe in the early phase of the pandemic and more dispersed exchanges involving several European countries from the second half of March 2020.

Vimentin binds to G-quadruplex repeats found at telomeres and gene promoters

G-quadruplex (G4) structures that can form at guanine-rich genomic sites, including telomeres and gene promoters, are actively involved in genome maintenance, replication, and transcription, through finely tuned interactions with protein networks. In the present study, we identified the intermediate filament protein Vimentin as a binder with nanomolar affinity for those G-rich sequences that give rise to at least two adjacent G4 units, named G4 repeats. This interaction is supported by the N-terminal domains of soluble Vimentin tetramers. The selectivity of Vimentin for G4 repeats versus individual G4s provides an unprecedented result. Based on GO enrichment analysis performed on genes having putative G4 repeats within their core promoters, we suggest that Vimentin recruitment at these sites may contribute to the regulation of gene expression during cell development and migration, possibly by reshaping the local higher-order genome topology, as already reported for lamin B.

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.

Aerobic pyruvate metabolism sensitizes cells to ferroptosis primed by GSH depletion

Ferroptosis is a non-accidental, regulated form of cell death operated by lipid peroxidation under strict control of GPx4 activity. This is consistent with the notion that lipid peroxidation is initiated by radicals produced from decomposition of traces of pre-existing lipid hydroperoxides. The question, therefore, emerges about the formation of these traces of lipid hydroperoxides interacting with Fe²⁺. In the most realistic option, they are produced by oxygen activated species generated during aerobic metabolism. Screening for metabolic sources of superoxide supporting ferroptosis induced by GSH depletion, we failed to detect, in our cell model, a role of respiratory chain. We observed instead that the pyruvate dehydrogenase complex -as other α keto acid dehydrogenases already known as a major source of superoxide in mitochondria- supports ferroptosis. The opposite effect on ferroptosis by silencing either the E1 or the E3 subunit of the pyruvate dehydrogenase complex pointed out the autoxidation of dihydrolipoamide as the source of superoxide. We finally observed that GSH depletion activates superoxide production, seemingly through the inhibition of the specific kinase that inhibits pyruvate dehydrogenase. In summary, this set of data is compatible with a scenario where the more electrophilic status produced by GSH depletion not only activates ferroptosis by preventing GPx4 activity, but also favors the formation of lipid hydroperoxides. In an attractive perspective of tissue homeostasis, it is the activation of energetic metabolism associated to a decreased nucleophilic tone that, besides supporting energy demanding proliferation, also sensitizes cells to a regulated form of death.

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.

Evolutionary Study of the Crassphage Virus at Gene Level

crAss-like viruses are a putative family of bacteriophages recently discovered. The eponym of the clade, crAssphage, is an enteric bacteriophage estimated to be present in at least half of the human population and it constitutes up to 90% of the sequences in some human fecal viral metagenomic datasets. We focused on the evolutionary dynamics of the genes encoded on the crAssphage genome. By investigating the conservation of the genes, a consistent variation in the evolutionary rates across the different functional groups was found. Gene duplications in crAss-like genomes were detected. By exploring the differences among the functional categories of the genes, we confirmed that the genes encoding capsid proteins were the most ubiquitous, despite their overall low sequence conservation. It was possible to identify a core of proteins whose evolutionary trees strongly correlate with each other, suggesting their genetic interaction. This group includes the capsid proteins, which are thus established as extremely suitable for rebuilding the phylogenetic tree of this viral clade. A negative correlation between the ubiquity and the conservation of viral protein sequences was shown. Together, this study provides an in-depth picture of the evolution of different genes in crAss-like viruses.

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.

Copper/Zinc Superoxide Dismutase from the Crocodile Icefish Chionodraco hamatus: Antioxidant Defense at Constant Sub-Zero Temperature

In the present study, we describe the purification and molecular characterization of Cu,Zn superoxide dismutase (SOD) from Chionodraco hamatus, an Antarctic teleost widely distributed in many areas of the Ross Sea that plays a pivotal role in the Antarctic food chain. The primary sequence was obtained using biochemical and molecular biology approaches and compared with Cu,Zn SODs from other organisms. Multiple sequence alignment using the amino acid sequence revealed that Cu,Zn SOD showed considerable sequence similarity with its orthologues from various vertebrate species, but also some specific substitutions directly linked to cold adaptation. Phylogenetic analyses presented the monophyletic status of Antartic Teleostei among the Perciformes, confirming the erratic differentiation of these proteins and concurring with the theory of the "unclock-like" behavior of Cu,Zn SOD evolution. Expression of C. hamatus Cu,Zn SOD at both the mRNA and protein levels were analyzed in various tissues, highlighting the regulation of gene expression related to environmental stress conditions and also animal physiology. The data presented are the first on the antioxidant enzymes of a fish belonging to the Channichthyidae family and represent an important starting point in understanding the antioxidant systems of these organisms that are subject to constant risk of oxidative stress.

Lack of glutathione peroxidase-8 in the ER impacts on lipid composition of HeLa cells microsomal membranes

GPx8 is a glutathione peroxidase homolog inserted in the membranes of endoplasmic reticulum (ER), where it seemingly plays a role in controlling redox status by preventing the spill of H₂O₂. We addressed the impact of GPx8 silencing on the lipidome of microsomal membranes, using stably GPx8-silenced HeLa cells. The two cell lines were clearly separated by Principal Component Analysis (PCA) and Partial Least Square Discriminant analysis (PLS-DA) of lipidome. Considering in detail the individual lipid classes, we observed that unsaturated glycerophospholipids (GPL) decreased, while only in phosphatidylinositols (PI) a substitution of monounsaturated fatty acids (MUFA) for polyunsaturated fatty acids (PUFA) was observed. Among sphingolipids (SL), ceramides (CER) decreased while sphingomyelins (SM) and neutral glycophingolipids (nGSL) increased. Here, in addition, longer chains than in controls in the amide fatty acid were present. The increase up to four folds of the CER (d18:1; c24:0) containing three hexose units, was the most remarkable species increasing in the differential lipidome of siGPx8 cells. Quantitative RT-PCR complied with lipidomic analysis specifically showing an increased expression of: i) acyl-CoA synthetase 5 (ACSL5); ii) CER synthase 2 and 4; iii) CER transporter (CERT); iv) UDP-glucosyl transferase (UDP-GlcT), associated to a decreased expression of UDP-galactosyl transferase (UDP-GalT). A role of the unfolded protein response (UPR) and the spliced form of the transcription factor XBP1 on the transcriptional changes of GPx8 silenced cells was ruled-out. Similarly, also the involvement of Nrf2 and NF-κB. Altogether our results indicate that GPx8-silencing of HeLa yields a membrane depleted by about 24% of polyunsaturated GPL and a corresponding increase of saturated or monounsaturated SM and specific nGSL. This is tentatively interpreted as an adaptive mechanism leading to an increased resistance to radical oxidations. Moreover, the marked shift of fatty acid composition of PI emerges as a possibly relevant issue in respect to the impact of GPx8 on signaling pathways.

Insight into the mechanism of ferroptosis inhibition by ferrostatin-1

Ferroptosis is a form of cell death primed by iron and lipid hydroperoxides and prevented by GPx4. Ferrostatin-1 (fer-1) inhibits ferroptosis much more efficiently than phenolic antioxidants. Previous studies on the antioxidant efficiency of fer-1 adopted kinetic tests where a diazo compound generates the hydroperoxyl radical scavenged by the antioxidant. However, this reaction, accounting for a chain breaking effect, is only minimally useful for the description of the inhibition of ferrous iron and lipid hydroperoxide dependent peroxidation. Scavenging lipid hydroperoxyl radicals, indeed, generates lipid hydroperoxides from which ferrous iron initiates a new peroxidative chain reaction. We show that when fer-1 inhibits peroxidation, initiated by iron and traces of lipid hydroperoxides in liposomes, the pattern of oxidized species produced from traces of pre-existing hydroperoxides is practically identical to that observed following exhaustive peroxidation in the absence of the antioxidant. This supported the notion that the anti-ferroptotic activity of fer-1 is actually due to the scavenging of initiating alkoxyl radicals produced, together with other rearrangement products, by ferrous iron from lipid hydroperoxides. Notably, fer-1 is not consumed while inhibiting iron dependent lipid peroxidation. The emerging concept is that it is ferrous iron itself that reduces fer-1 radical. This was supported by electroanalytical evidence that fer-1 forms a complex with iron and further confirmed in cells by fluorescence of calcein, indicating a decrease of labile iron in the presence of fer-1. The notion of such as pseudo-catalytic cycle of the ferrostatin-iron complex was also investigated by means of quantum mechanics calculations, which confirmed the reduction of an alkoxyl radical model by fer-1 and the reduction of fer-1 radical by ferrous iron. In summary, GPx4 and fer-1 in the presence of ferrous iron, produces, by distinct mechanism, the most relevant anti-ferroptotic effect, i.e the disappearance of initiating lipid hydroperoxides.

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.

Molecular characterization of novel mitochondrial peroxiredoxins from the Antarctic emerald rockcod and their gene expression in response to environmental warming

In the present study we describe the molecular characterization of the two paralogous mitochondrial peroxiredoxins from Trematomus bernacchii, a teleost that plays a pivotal role in the Antarctic food chain. The two putative amino acid sequences were compared with orthologs from other fish, highlighting a high percentage of identity and similarity with the respective variant, in particular for the residues that are essential for the characteristic peroxidase activity of these enzymes. The temporal expression of Prdx3 and Prdx5 mRNAs in response to short-term thermal stress showed a general upregulation of prdx3, suggesting that this isoform is the most affected by temperature increase. These data, together with the peculiar differences between the molecular structures of the two mitochondrial Prdxs in T. bernacchii as well as in the tropical species Stegastes partitus, suggest an adaptation that allowed these poikilothermic aquatic vertebrates to colonize very different environments, characterized by different temperature ranges.

Loss of cardiac Wnt/β-catenin signalling in desmoplakin-deficient AC8 zebrafish models is rescuable by genetic and pharmacological intervention

Aims

Arrhythmogenic cardiomyopathy (AC) is an inherited heart disease characterized by life-threatening ventricular arrhythmias and fibro-fatty replacement of the myocardium. More than 60% of AC patients show pathogenic mutations in genes encoding for desmosomal proteins. By focusing our attention on the AC8 form, linked to the junctional protein desmoplakin (DSP), we present here a zebrafish model of DSP deficiency, exploited to identify early changes of cell signalling in the cardiac region.

Methods and results

To obtain an embryonic model of Dsp deficiency, we first confirmed the orthologous correspondence of zebrafish Dsp genes (dspa and dspb) to the human DSP counterpart. Then, we verified their cardiac expression, at embryonic and adult stages, and subsequently we targeted them by antisense morpholino strategy, confirming specific and disruptive effects on desmosomes, like those identified in AC patients. Finally, we exploited our Dsp-deficient models for an in vivo cell signalling screen, using pathway-specific reporter transgenes. Out of nine considered, three pathways (Wnt/β-catenin, TGFβ/Smad3, and Hippo/YAP-TAZ) were significantly altered, with Wnt as the most dramatically affected. Interestingly, under persistent Dsp deficiency, Wnt signalling is rescuable both by a genetic and a pharmacological approach.

Conclusion

Our data point to Wnt/β-catenin as the final common pathway underlying different desmosomal AC forms and support the zebrafish as a suitable model for detecting early signalling pathways involved in the pathogenesis of DSP-associated diseases, possibly responsive to pharmacological or genetic rescue.

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.

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.

Glutathione peroxidase 4-catalyzed reduction of lipid hydroperoxides in membranes: The polar head of membrane phospholipids binds the enzyme and addresses the fatty acid hydroperoxide group toward the redox center

GPx4 is a monomeric glutathione peroxidase, unique in reducing the hydroperoxide group (-OOH) of fatty acids esterified in membrane phospholipids. This reaction inhibits lipid peroxidation and accounts for enzyme's vital role. Here we investigated the interaction of GPx4 with membrane phospholipids. A cationic surface near the GPx4 catalytic center interacts with phospholipid polar heads. Accordingly, SPR analysis indicates cardiolipin as the phospholipid with maximal affinity to GPx4. Consistent with the electrostatic nature of the interaction, KCl increases the K_D. Molecular dynamic (MD) simulation shows that a -OOH posed in the core of the membrane as 13 - or 9 -OOH of tetra-linoleoyl cardiolipin or 15 -OOH stearoyl-arachidonoyl-phosphaphatidylcholine moves to the lipid-water interface. Thereby, the -OOH groups are addressed toward the GPx4 redox center. In this pose, however, the catalytic site facing the membrane would be inaccessible to GSH, but the consecutive redox processes facilitate access of GSH, which further primes undocking of the enzyme, because GSH competes for the binding residues implicated in docking. During the final phase of the catalytic cycle, while GSSG is produced, GPx4 is disconnected from the membrane. The observation that GSH depletion in cells induces GPx4 translocation to the membrane, is in agreement with this concept.

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.

Sox10 contributes to the balance of fate choice in dorsal root ganglion progenitors

The development of functional peripheral ganglia requires a balance of specification of both neuronal and glial components. In the developing dorsal root ganglia (DRGs), these components form from partially-restricted bipotent neuroglial precursors derived from the neural crest. Work in mouse and chick has identified several factors, including Delta/Notch signaling, required for specification of a balance of these components. We have previously shown in zebrafish that the Sry-related HMG domain transcription factor, Sox10, plays an unexpected, but crucial, role in sensory neuron fate specification in vivo. In the same study we described a novel Sox10 mutant allele, sox10baz1, in which sensory neuron numbers are elevated above those of wild-types. Here we investigate the origin of this neurogenic phenotype. We demonstrate that the supernumerary neurons are sensory neurons, and that enteric and sympathetic neurons are almost absent just as in classical sox10 null alleles; peripheral glial development is also severely abrogated in a manner similar to other sox10 mutant alleles. Examination of proliferation and apoptosis in the developing DRG reveals very low levels of both processes in wild-type and sox10baz1, excluding changes in the balance of these as an explanation for the overproduction of sensory neurons. Using chemical inhibition of Delta-Notch-Notch signaling we demonstrate that in embryonic zebrafish, as in mouse and chick, lateral inhibition during the phase of trunk DRG development is required to achieve a balance between glial and neuronal numbers. Importantly, however, we show that this mechanism is insufficient to explain quantitative aspects of the baz1 phenotype. The Sox10(baz1) protein shows a single amino acid substitution in the DNA binding HMG domain; structural analysis indicates that this change is likely to result in reduced flexibility in the HMG domain, consistent with sequence-specific modification of Sox10 binding to DNA. Unlike other Sox10 mutant proteins, Sox10(baz1) retains an ability to drive neurogenin1 transcription. We show that overexpression of neurogenin1 is sufficient to produce supernumerary DRG sensory neurons in a wild-type background, and can rescue the sensory neuron phenotype of sox10 morphants in a manner closely resembling the baz1 phenotype. We conclude that an imbalance of neuronal and glial fate specification results from the Sox10(baz1) protein's unique ability to drive sensory neuron specification whilst failing to drive glial development. The sox10baz1 phenotype reveals for the first time that a Notch-dependent lateral inhibition mechanism is not sufficient to fully explain the balance of neurons and glia in the developing DRGs, and that a second Sox10-dependent mechanism is necessary. Sox10 is thus a key transcription factor in achieving the balance of sensory neuronal and glial fates.

Redox status in a model of cancer stem cells

Reversible oxidation of Cys residues is a crucial element of redox homeostasis and signaling. According to a popular concept in oxidative stress signaling, the oxidation of targets of signals can only take place following an overwhelming of the cellular antioxidant capacity. This concept, however, ignores the activation of feedback mechanisms possibly leading to a paradoxical effect. In a model of cancer stem cells (CSC), stably overexpressing the TAZ oncogene, we observed that the increased formation of oxidants is associated with a globally more reduced state of proteins. Redox proteomics revealed that several proteins, capable of undergoing reversible redox transitions, are indeed more reduced while just few are more oxidized. Among the proteins more oxidized, G6PDH emerges as both more expressed and activated by oxidation. This accounts for the observed more reduced state of the NADPH/NADP⁺ couple. The dynamic redox flux generating this apparently paradoxical effect is rationalized in a computational system biology model highlighting the crucial role of G6PDH activity on the rate of redox transitions eventually leading to the reduction of reversible redox switches.

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.

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.

Peroxiredoxin 6 from the Antarctic emerald rockcod: molecular characterization of its response to warming

In the present study, we describe the purification and molecular characterization of two peroxiredoxins (Prdxs), referred to as Prdx6A and Prdx6B, from Trematomus bernacchii, a teleost widely distributed in many areas of Antarctica, that plays a pivotal role in the Antarctic food chain. The two putative amino acid sequences were compared with Prdx6 orthologs from other fish, highlighting a high percentage of identity and similarity with the respective variant, in particular for the residues that are essential for the characteristic peroxidase and phospholipase activities of these enzymes. Phylogenetic analyses suggest the appearance of the two prdx6 genes through a duplication event before the speciation that led to the differentiation of fish families and that the evolution of the two gene variants seems to proceed together with the evolution of fish orders and families. The temporal expression of Prdx6 mRNA in response to short-term thermal stress showed a general upregulation of prdx6b and inhibition of prdx6a, suggesting that the latter is the variant most affected by temperature increase. The variations of mRNA accumulation are more conspicuous in heart than the liver, probably related to behavioral changes of the specimens in response to elevated temperature. These data, together with the peculiar differences between the molecular structures of the two Prdx6s in T. bernacchii as well as in the tropical species Stegastes partitus, suggest an adaptation that allowed these poikilothermic aquatic vertebrates to colonize very different environments, characterized by different temperature ranges.

Selenocysteine oxidation in glutathione peroxidase catalysis: an MS-supported quantum mechanics study

Glutathione peroxidases (GPxs) are enzymes working with either selenium or sulfur catalysis. They adopted diverse functions ranging from detoxification of H(2)O(2) to redox signaling and differentiation. The relative stability of the selenoenzymes, however, remained enigmatic in view of the postulated involvement of a highly unstable selenenic acid form during catalysis. Nevertheless, density functional theory calculations obtained with a representative active site model verify the mechanistic concept of GPx catalysis and underscore its efficiency. However, they also allow that the selenenic acid, in the absence of the reducing substrate, reacts with a nitrogen in the active site. MS/MS analysis of oxidized rat GPx4 complies with the predicted structure, an 8-membered ring, in which selenium is bound as selenenylamide to the protein backbone. The intermediate can be re-integrated into the canonical GPx cycle by glutathione, whereas, under denaturing conditions, its selenium moiety undergoes β-cleavage with formation of a dehydro-alanine residue. The selenenylamide bypass prevents destruction of the redox center due to over-oxidation of the selenium or its elimination and likely allows fine-tuning of GPx activity or alternate substrate reactions for regulatory purposes.

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.

Understanding mammalian glutathione peroxidase 7 in the light of its homologs

The glutathione peroxidase homologs (GPxs) efficiently reduce hydroperoxides using electrons from glutathione (GSH), thioredoxin (Trx), or protein disulfide isomerase (PDI). Trx is preferentially used by the GPxs of the majority of bacteria, invertebrates, plants, and fungi. GSH or PDI, instead, is preferentially used by vertebrate GPxs that operate by Sec or Cys catalysis, respectively. Mammalian GPx7 and GPx8 are unique homologs that contain a peroxidatic Cys (CP). Being reduced by PDI and located within the endoplasmic reticulum (ER), these enzymes have been involved in oxidative protein folding. Kinetic analysis indicates that oxidation of PDI by recombinant GPx7 occurs at a much faster rate than that of GSH. Nonetheless, activity measurement suggests that, at physiological concentrations, a competition between these two substrates takes place, with the rate of PDI oxidation by GPx7 controlled by the concentration of GSH, whereas the GSSG produced in the competing reaction contributes to the ER redox buffer. A mechanism has been proposed for GPx7 involving two Cys residues, in which an intramolecular disulfide of the CP is formed with an alleged resolving Cys (CR) located in the strongly conserved FPCNQ motif (C86 in humans), a noncanonical position in GPxs. Kinetic measurements and comparison with the other thiol peroxidases containing a functional CR suggest that a resolving function of C86 in the catalytic cycle is very unlikely. We propose that GPx7 is catalytically active as a 1-Cys-GPx, in which CP both reduces H2O2 and oxidizes PDI, and that the CP-C86 disulfide has instead the role of stabilizing the oxidized peroxidase in the absence of the reducing substrate.

Glutathione peroxidase 8 is transcriptionally regulated by HIFα and modulates growth factor signaling in HeLa cells

GPx8 is a mammalian Cys-glutathione peroxidase of the endoplasmic reticulum membrane, involved in protein folding. Its regulation is mostly unknown. We addressed both the functionality of two hypoxia-response elements (HREs) within the promoter, GPx8 HRE1 and GPx8 HRE2, and the GPx8 physiological role. In HeLa cells, treatment with HIFα stabilizers, such as diethyl succinate (DES) or 2-2'-bipyridyl (BP), induces GPx8 expression at both mRNA and protein level. Luciferase activity of pGL3(GPx8wt), containing a fragment of the GPx8 promoter including the two HREs, is also induced by DES/BP or by overexpressing either individual HIFα subunit. Mutating GPx8 HRE1 within pGL3(GPx8wt) resulted in a significantly higher inhibition of luciferase activity than mutating GPx8 HRE2. Electrophoretic mobility-shift assay showed that both HREs exhibit enhanced binding to a nuclear extract from DES/BP-treated cells, with stronger binding by GPx8 HRE1. In DES-treated cells transfected with pGL3(GPx8wt) or mutants thereof, silencing of HIF2α, but not HIF1α, abolishes luciferase activity. Thus GPx8 is a novel HIF target preferentially responding to HIF2α binding at its two novel functional GPx8 HREs, with GPx8 HRE1 playing the major role. Fibroblast growth factor (FGF) treatment increases GPx8 mRNA expression, and reporter gene experiments indicate that induction occurs via HIF. Comparing the effects of depleting GPx8 on the downstream effectors of FGF or insulin signaling revealed that absence of GPx8 results in a 16- or 12-fold increase in phosphorylated ERK1/2 by FGF or insulin treatment, respectively. Furthermore, in GPx8-depleted cells, phosphorylation of AKT by insulin treatment increases 2.5-fold. We suggest that induction of GPx8 expression by HIF slows down proliferative signaling during hypoxia and/or growth stimulation through receptor tyrosine kinases.

Quantitative label-free redox proteomics of reversible cysteine oxidation in red blood cell membranes

Reversible oxidation of cysteine residues is a relevant posttranslational modification of proteins. However, the low activation energy and transitory nature of the redox switch and the intrinsic complexity of the analysis render quite challenging the aim of a rigorous high-throughput screening of the redox status of redox-sensitive cysteine residues. We describe here a quantitative workflow for redox proteomics, where the ratio between the oxidized forms of proteins in the control vs treated samples is determined by a robust label-free approach. We critically present the convenience of the procedure by specifically addressing the following aspects: (i) the accurate ratio, calculated from the whole set of identified peptides rather than just isotope-tagged fragments; (ii) the application of a robust analytical pipeline to frame the most consistent data averaged over the biological variability; (iii) the relevance of using stringent criteria of analysis, even at the cost of losing potentially interesting but statistically uncertain data. The pipeline has been assessed on red blood cell membrane challenged with diamide as a model of a mild oxidative condition. The cluster of identified proteins encompassed components of the cytoskeleton more oxidized. Indirectly, our analysis confirmed the previous observation that oxidized hemoglobin binds to membranes while oxidized peroxiredoxin 2 loses affinity.

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.

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.

Toward more transparent and reproducible omics studies through a common metadata checklist and data publications

Biological processes are fundamentally driven by complex interactions between biomolecules. Integrated high-throughput omics studies enable multifaceted views of cells, organisms, or their communities. With the advent of new post-genomics technologies, omics studies are becoming increasingly prevalent; yet the full impact of these studies can only be realized through data harmonization, sharing, meta-analysis, and integrated research. These essential steps require consistent generation, capture, and distribution of metadata. To ensure transparency, facilitate data harmonization, and maximize reproducibility and usability of life sciences studies, we propose a simple common omics metadata checklist. The proposed checklist is built on the rich ontologies and standards already in use by the life sciences community. The checklist will serve as a common denominator to guide experimental design, capture important parameters, and be used as a standard format for stand-alone data publications. The omics metadata checklist and data publications will create efficient linkages between omics data and knowledge-based life sciences innovation and, importantly, allow for appropriate attribution to data generators and infrastructure science builders in the post-genomics era. We ask that the life sciences community test the proposed omics metadata checklist and data publications and provide feedback for their use and improvement.

Organochalcogen peroxidase mimetics as potential drugs: a long story of a promise still unfulfilled

Organochalcogen compounds have attracted the interest of a multitude of studies to design potential therapeutic agents mimicking the peroxidase activity of selenium-based glutathione peroxidases (GPx's). Starting from the pioneering ebselen, various compounds have been synthesized over the years, which may be traced in three major classes of molecules: cyclic selenenyl amides, diaryl diselenides, and aromatic or aliphatic monoselenides. These compounds share common features and determinants needed to exert an efficient GPx-like activity, such as polarizing groups in close proximity to selenium and steric effects. Nonetheless, the reactivity of selenium, and tellurium as well, poses serious problems for the predictability of the biological effects of these compounds in vivo when used as potential drugs. These molecules, indeed, interfere with thiols of redox-regulated proteins and enzymes, leading to unexpected biological effects. The various chemical aspects of the reaction mechanism of peroxidase mimetics are surveyed here, focusing on experimental evidence and quantum mechanics calculations of organochalcogen representatives of the various classes.