Selecton 2007: advanced models for detecting positive and purifying selection using a Bayesian inference approach

Biologically significant sites in a protein may be identified by contrasting the rates of synonymous (Ks) and non-synonymous (Ka) substitutions. This enables the inference of site-specific positive Darwinian selection and purifying selection. We present here Selecton version 2.2 (http://selecton.bioinfo.tau.ac.il), a web server which automatically calculates the ratio between Ka and Ks (ω) at each site of the protein. This ratio is graphically displayed on each site using a color-coding scheme, indicating either positive selection, purifying selection or lack of selection. Selecton implements an assembly of different evolutionary models, which allow for statistical testing of the hypothesis that a protein has undergone positive selection. Specifically, the recently developed mechanistic-empirical model is introduced, which takes into account the physicochemical properties of amino acids. Advanced options were introduced to allow maximal fine tuning of the server to the user's specific needs, including calculation of statistical support of the ω values, an advanced graphic display of the protein's 3-dimensional structure, use of different genetic codes and inputting of a pre-built phylogenetic tree. Selecton version 2.2 is an effective, user-friendly and freely available web server which implements up-to-date methods for computing site-specific selection forces, and the visualization of these forces on the protein's sequence and structure

Mapping of Tilapia Lake Virus entry pathways with inhibitors reveals dependence on dynamin activity and cholesterol but not endosomal acidification

Tilapia Lake Virus (TiLV) is an emerging virus lethal to tilapia, which threatens the global tilapia aquaculture with severe implications for food security. TiLV possesses similar features to orthomyxoviruses but is classified in the sole and the monotypic genus Tilapinevirus of the family Amnoonviridae. TiLV enveloped virions encapsidate a genome comprising ten segments of single-stranded, negative RNA. Remarkably, nine of TiLV’s ten major proteins lack sequence homology to any known viral or cellular proteins. The mode of TiLV entry into tilapia cells is not known. Following the measurement of the entry window of TiLV (∼3 h), we applied a panel of inhibitors of known regulators of endocytic functions to map the molecular requirements for TiLV entry. We identified productive entry by quantification of TiLV nucleoprotein expression and the generation of infectious particles. Inhibition of dynamin activity with dynasore or dynole, or depletion of cholesterol with methyl-β-cyclodextrin, strongly inhibited TiLV protein synthesis and infectious virion production. Moreover, inhibition of actin cytoskeleton polymerization with latrunculin A or microtubule polymerization with nocodazole within the entry window resulted in partial inhibition of TiLV infection. In contrast, inhibitors of endosomal acidification (NH4Cl, bafilomycin A1, or chloroquine), an inhibitor of clathrin-coated pit assembly (pitstop 2), and erlotinib—an inhibitor of the endocytic Cyclin G-associated kinase (GAK), did not affect TiLV entry. Altogether, these results suggest that TiLV enters via dynamin-mediated endocytosis in a cholesterol-, cytoskeleton-dependent manner, and clathrin-, pH-independent manner. Thus, despite being an orthomyx-like virus, when compared to the prototypical orthomyxovirus (influenza A virus), TiLV shows a distinct set of requirements for entry into cells

Synonymous site conservation in the HIV-1 genome

Background: Synonymous or silent mutations are usually thought to evolve neutrally. However, accumulating recent evidence has demonstrated that silent mutations may destabilize RNA structures or disrupt cis regulatory motifs superimposed on coding sequences. Such observations suggest the existence of stretches of codon sites that are evolutionary conserved at both DNA-RNA and protein levels. Such stretches may point to functionally important regions within protein coding sequences not necessarily reflecting functional constraints on the amino-acid sequence. The HIV-1 genome is highly compact, and often harbors overlapping functional elements at the protein, RNA, and DNA levels. This superimposition of functions leads to complex selective forces acting on all levels of the genome and proteome. Considering the constraints on HIV-1 to maintain such a highly compact genome, we hypothesized that stretches of synonymous conservation would be common within its genome. Results: We used a combined computational-experimental approach to detect and characterize regions exhibiting strong purifying selection against synonymous substitutions along the HIV-1 genome. Our methodology is based on advanced probabilistic evolutionary models that explicitly account for synonymous rate variation among sites and rate dependencies among adjacent sites. These models are combined with a randomization procedure to automatically identify the most statistically significant regions of conserved synonymous sites along the genome. Using this procedure we identified 21 conserved regions. Twelve of these are mapped to regions within overlapping genes, seven correlate with known functional elements, while the functions of the remaining four are yet unknown. Among these four regions, we chose the one that deviates most from synonymous rate homogeneity for in-depth computational and experimental characterization. In our assays aiming to quantify viral fitness in both early and late stages of the replication cycle, no differences were observed between the mutated and the wild type virus following the introduction of synonymous mutations. Conclusions: The contradiction between the inferred purifying selective forces and the lack of effect of these mutations on viral replication may be explained by the fact that the phenotype was measured in single-cycle infection assays in cell culture. Such a system does not account for the complexity of HIV-1 infections in vivo, which involves multiple infection cycles and interaction with the host immune system

Detection of Tilapia Lake Virus (TiLV) in Clinical Samples by Culturing and Nested RT-PCR

Tilapia are an important group of farmed fish that serve as a significant protein source worldwide. In recent years, substantial mortality of wild tilapia has been observed in the Sea of Galilee and in commercial ponds in Israel and Ecuador. We have identified the etiological agent of these mass die-offs as a novel orthomyxo-like virus and named it tilapia lake virus (TiLV). Here, we provide the conditions for efficient isolation, culturing, and quantification of the virus, including the use of susceptible fish cell lines. Moreover, we describe a sensitive nested reverse transcription-PCR (RT-PCR) assay allowing the rapid detection of TiLV in fish organs. This assay revealed, for the first time to our knowledge, the presence of TiLV in diseased Colombian tilapia, indicating a wider distribution of this emerging pathogen and stressing the risk that TiLV poses for the global tilapia industry. Overall, the described procedures should provide the tilapia aquaculture industry with important tools for the detection and containment of this pathogen

Characterization of a Novel Orthomyxo-like Virus Causing Mass Die-Offs of Tilapia

Tilapia are an important global food source due to their omnivorous diet, tolerance for high-density aquaculture, and relative disease resistance. Since 2009, tilapia aquaculture has been threatened by mass die-offs in farmed fish in Israel and Ecuador. Here we report evidence implicating a novel orthomyxo-like virus in these outbreaks. The tilapia lake virus (TiLV) has a 10-segment, negative-sense RNA genome. The largest segment, segment 1, contains an open reading frame with weak sequence homology to the influenza C virus PB1 subunit. The other nine segments showed no homology to other viruses but have conserved, complementary sequences at their 5′ and 3′ termini, consistent with the genome organization found in other orthomyxoviruses. In situ hybridization indicates TiLV replication and transcription at sites of pathology in the liver and central nervous system of tilapia with disease

The Gag Cleavage Product, p12, is a Functional Constituent of the Murine Leukemia Virus Pre-Integration Complex

The p12 protein is a cleavage product of the Gag precursor of the murine leukemia virus (MLV). Specific mutations in p12 have been described that affect early stages of infection, rendering the virus replication-defective. Such mutants showed normal generation of genomic DNA but no formation of circular forms, which are markers of nuclear entry by the viral DNA. This suggested that p12 may function in early stages of infection but the precise mechanism of p12 action is not known. To address the function and follow the intracellular localization of the wt p12 protein, we generated tagged p12 proteins in the context of a replication-competent virus, which allowed for the detection of p12 at early stages of infection by immunofluorescence. p12 was found to be distributed to discrete puncta, indicative of macromolecular complexes. These complexes were localized to the cytoplasm early after infection, and thereafter accumulated adjacent to mitotic chromosomes. This chromosomal accumulation was impaired for p12 proteins with a mutation that rendered the virus integration-defective. Immunofluorescence demonstrated that intracellular p12 complexes co-localized with capsid, a known constituent of the MLV pre-integration complex (PIC), and immunofluorescence combined with fluorescent in situ hybridization (FISH) revealed co-localization of the p12 proteins with the incoming reverse transcribed viral DNA. Interactions of p12 with the capsid and with the viral DNA were also demonstrated by co-immunoprecipitation. These results imply that p12 proteins are components of the MLV PIC. Furthermore, a large excess of wt PICs did not rescue the defect in integration of PICs derived from mutant p12 particles, demonstrating that p12 exerts its function as part of this complex. Altogether, these results imply that p12 proteins are constituent of the MLV PIC and function in directing the PIC from the cytoplasm towards integration

Evolutionary Modeling of Rate Shifts Reveals Specificity Determinants in HIV-1 Subtypes

A hallmark of the human immunodeficiency virus 1 (HIV-1) is its rapid rate of evolution within and among its various subtypes. Two complementary hypotheses are suggested to explain the sequence variability among HIV-1 subtypes. The first suggests that the functional constraints at each site remain the same across all subtypes, and the differences among subtypes are a direct reflection of random substitutions, which have occurred during the time elapsed since their divergence. The alternative hypothesis suggests that the functional constraints themselves have evolved, and thus sequence differences among subtypes in some sites reflect shifts in function. To determine the contribution of each of these two alternatives to HIV-1 subtype evolution, we have developed a novel Bayesian method for testing and detecting site-specific rate shifts. The RAte Shift EstimatoR (RASER) method determines whether or not site-specific functional shifts characterize the evolution of a protein and, if so, points to the specific sites and lineages in which these shifts have most likely occurred. Applying RASER to a dataset composed of large samples of HIV-1 sequences from different group M subtypes, we reveal rampant evolutionary shifts throughout the HIV-1 proteome. Most of these rate shifts have occurred during the divergence of the major subtypes, establishing that subtype divergence occurred together with functional diversification. We report further evidence for the emergence of a new sub-subtype, characterized by abundant rate-shifting sites. When focusing on the rate-shifting sites detected, we find that many are associated with known function relating to viral life cycle and drug resistance. Finally, we discuss mechanisms of covariation of rate-shifting sites

Oncolytic Virotherapy: The Cancer Cell Side

Cell autonomous immunity genes mediate the multiple stages of anti-viral defenses, including recognition of invading pathogens, inhibition of viral replication, reprogramming of cellular metabolism, programmed-cell-death, paracrine induction of antiviral state, and activation of immunostimulatory inflammation. In tumor development and/or immunotherapy settings, selective pressure applied by the immune system results in tumor immunoediting, a reduction in the immunostimulatory potential of the cancer cell. This editing process comprises the reduced expression and/or function of cell autonomous immunity genes, allowing for immune-evasion of the tumor while concomitantly attenuating anti-viral defenses. Combined with the oncogene-enhanced anabolic nature of cancer-cell metabolism, this attenuation of antiviral defenses contributes to viral replication and to the selectivity of oncolytic viruses (OVs) towards malignant cells. Here, we review the manners by which oncogene-mediated transformation and tumor immunoediting combine to alter the intracellular milieu of tumor cells, for the benefit of OV replication. We also explore the functional connection between oncogenic signaling and epigenetic silencing, and the way by which restriction of such silencing results in immune activation. Together, the picture that emerges is one in which OVs and epigenetic modifiers are part of a growing therapeutic toolbox that employs activation of anti-tumor immunity for cancer therapy