Molecular evolution of the tick-borne encephalitis and Powassan viruses

Issues associated with newly emerging viruses, their genetic diversity, and viral evolution in modern environments are currently attracting growing attention. In this study, a phylogenetic analysis was performed and the evolution rate was evaluated for such pathogenic flaviviruses endemic to Russia as tick-borne encephalitis virus (TBEV) and Powassan virus (PV). The analysis involved 47 nucleotide sequences of the TBEV genome region encoding protein E and 17 sequences of the PV NS5-encoding region. The nucleotide substitution rate was estimated as 1.4 × 10⁻⁴ and 5.4 × 10⁻⁵ substitutions per site per year for the E protein-encoding region of the TBEV genome and for the NS5 genome region of PV, respectively. The ratio of non-synonymous to synonymous nucleotide substitutions (dN/dS) in viral sequences was calculated as 0.049 for TBEV and 0.098 for PV. The highest dN/dS values of 0.201–0.220 were found in the subcluster of Russian and Canadian PV strains, and the lowest value of 0.024 was observed in the cluster of Russian and Chinese strains of the Far Eastern TBEV genotype. Evaluation of time intervals between the events of viral evolution showed that the European subtype of TBEV diverged from the common TBEV ancestor approximately 2750 years ago, while the Siberian and Far Eastern subtypes emerged approximately 2250 years ago. The PV was introduced into its natural foci of the Russian Primorskii krai only approximately 70 years ago; these strains were very close to Canadian PV strains. The pattern of PV evolution in North America was similar to the evolution of the Siberian and Far Eastern TBEV subtypes in Asia. The moments of divergence between major genetic groups of TBEV and PV coincide with historical periods of climate warming and cooling, suggesting that climate change was a key factor in the evolution of flaviviruses in past millennia.     

Multiple transmissions of tick-borne encephalitis virus between Japan and Russia

Tick-borne encephalitis (TBE) is a zoonotic disease causing meningitis, encephalitis, and meningoencephalitis. Tick-borne encephalitis virus (TBEV) is an etiological agent of TBE. From an analysis of five distinct sequences of Japanese TBEV, it has been proposed that Japanese TBEV was transmitted from Russia to Japan on just a single occasion 260-430 years ago. Here thirteen distinct nucleotide sequences encoding the entire region of the envelope protein for Japanese TBEV were analyzed. It is shown, from the phylogenetic analysis, that Japanese TBEV belongs to the Far Eastern subtype, which is known to be highly pathogenic. Japanese TBEV was divided into three groups, and TBEV was inferred to have been transmitted between Japan and Russia at least three times, which were estimated to have occurred several hundred years ago. These results indicate that TBEV has not only been endemic but also transmitted multiple times to Japan.     

Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

Molecular typing of the tick-borne encephalitis virus isolated from patients with different-infection severities in the south of Russia's Far East

The thick-borne encephalitis virus (TBEV), which is widespread in the Eurasian continent, belongs to the Flaviviridae family, Flavirus genus, and comprises the Far Eastern, Siberian and West European subtypes. It was for the first time that the gene part of the E 24 strain envelope glycoprotein of TBEV, which caused infection in residents of the South of Russia's Far East, was analyzed. It was established that the TBEV Far-Eastern subtype causes different-severity disease cases ranging from the focal ones with the lethal outcome to latent infection forms. On the basis of the phylogenetic analysis, the Far-Eastern subtype was shared between 4 sub-clusters, 2 of which constitute a majority of the analyzed TBEV strains.     

Whole-Genome Sequencing of Theileria parva Strains Provides Insight into Parasite Migration and Diversification in the African Continent

The disease caused by the apicomplexan protozoan parasite Theileria parva, known as East Coast fever or Corridor disease, is one of the most serious cattle diseases in Eastern, Central, and Southern Africa. We performed whole-genome sequencing of nine T. parva strains, including one of the vaccine strains (Kiambu 5), field isolates from Zambia, Uganda, Tanzania, or Rwanda, and two buffalo-derived strains. Comparison with the reference Muguga genome sequence revealed 34 814–121 545 single nucleotide polymorphisms (SNPs) that were more abundant in buffalo-derived strains. High-resolution phylogenetic trees were constructed with selected informative SNPs that allowed the investigation of possible complex recombination events among ancestors of the extant strains. We further analysed the dN/dS ratio (non-synonymous substitutions per non-synonymous site divided by synonymous substitutions per synonymous site) for 4011 coding genes to estimate potential selective pressure. Genes under possible positive selection were identified that may, in turn, assist in the identification of immunogenic proteins or vaccine candidates. This study elucidated the phylogeny of T. parva strains based on genome-wide SNPs analysis with prediction of possible past recombination events, providing insight into the migration, diversification, and evolution of this parasite species in the African continent.     

Molecular evolution of hepatitis B virus over 25 years

Determining the longitudinal molecular evolution of hepatitis B virus (HBV) is difficult due to HBV's genomic complexity and the need to study paired samples collected over long periods of time. In this study, serial samples were collected from eight hepatitis B virus e antigen-negative asymptomatic carriers of HBV genotype B in 1979 and 2004, thus providing a 25-year period to document the long-term molecular evolution of HBV. The rate, nature, and distribution of mutations that emerged over 25 years were determined by phylogenetic and linear regression analysis of full-length HBV genome sequences. Nucleotide hypervariability was observed within the polymerase and pre-S/S overlap region and within the core gene. The calculated mean number of nucleotide substitutions/site/year (7.9 x 10(-5)) was slightly higher than published estimates (1.5 x 10(-5) to 5 x 10(-5)). Nucleotide changes in the quasispecies population did not significantly alter the molecular evolutionary rate, based on linear regression analysis of evolutionary distances among serial clone pre-S region sequences. Therefore, the directly amplified or dominant sequence was sufficient to estimate the putative molecular evolutionary rate for these long-term serial samples. On average, the ratio of synonymous (dS) to nonsynonymous (dN) substitutions was highest for the polymerase-coding region and lowest for the core-coding region. The low dS/dN ratios observed within the core suggest that selection occurs within this gene region, possibly as an immune evasion strategy. The results of this study suggest that HBV sequence divergence may occur more rapidly than previously estimated, in a host immune phase-dependent manner.     

Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulata

This is the first study to quantify genomic sequence variation of the major histocompatibility complex (MHC) in wild and ornamental guppies, Poecilia reticulata. We sequenced 196-219 bp of exon 2 MHC class IIB (DAB) in 56 wild Trinidadian guppies and 14 ornamental strain guppies. Each of two natural populations possessed high allelic richness (15-16 alleles), whereas only three or fewer DAB alleles were amplified from ornamental guppies. The disparity in allelic richness between wild and ornamental fish cannot be fully explained by fixation of alleles by inbreeding, nor by the presence of non-amplified sequences (ie null alleles). Rather, we suggest that the same allele is fixed at duplicated MHC DAB loci owing to gene conversion. Alternatively, the number of loci in the ornamental strains has contracted during >100 generations in captivity, a hypothesis consistent with the accordion model of MHC evolution. We furthermore analysed the substitution patterns by making pairwise comparisons of sequence variation at the putative peptide binding region (PBR). The rate of non-synonymous substitutions (dN) only marginally exceeded synonymous substitutions (dS) in PBR codons. Highly diverged sequences showed no evidence for diversifying selection, possibly because synonymous substitutions have accumulated since their divergence. Also, the substitution pattern of similar alleles did not show evidence for diversifying selection, plausibly because advantageous non-synonymous substitutions have not yet accumulated. Intermediately diverged sequences showed the highest relative rate of non-synonymous substitutions, with dN/dS>14 in some pairwise comparisons. Consequently, a curvilinear relationship was observed between the dN/dS ratio and the level of sequence divergence.     

The Population Genetics of dN/dS

Evolutionary pressures on proteins are often quantified by the ratio of substitution rates at non-synonymous and synonymous sites. The dN/dS ratio was originally developed for application to distantly diverged sequences, the differences among which represent substitutions that have fixed along independent lineages. Nevertheless, the dN/dS measure is often applied to sequences sampled from a single population, the differences among which represent segregating polymorphisms. Here, we study the expected dN/dS ratio for samples drawn from a single population under selection, and we find that in this context, dN/dS is relatively insensitive to the selection coefficient. Moreover, the hallmark signature of positive selection over divergent lineages, dN/dS>1, is violated within a population. For population samples, the relationship between selection and dN/dS does not follow a monotonic function, and so it may be impossible to infer selection pressures from dN/dS. These results have significant implications for the interpretation of dN/dS measurements among population-genetic samples.     

Nucleotide sequence of the genome and complete amino acid sequence of a polyprotein of the tick-borne encephalitis virus

We have cloned and sequenced RNA encoding all virion and nonstructural proteins of tick-borne encephalitis virus (TBEV). Its length is 10,477 bases with a single open reading frame (nucleotides 127-10,363) encoding 3412 amino acids. The 5'- and 3'-noncoding regions have stem- and- loop structure. The polyprotein precursor is proteolytically cleaved, apparently, by a mechanism resembling that proposed for the expression of polyproteins of other flaviviruses, such as yellow fever, West Nile and Kunjin viruses. The deduced TBEV gene order is 5'-C-preM (M)-E-NS1-NS2A-NS2B-NS3-NS4A-NS4B-NS5++ +-3'. The genome and the polyprotein of TBEV and other flaviviruses appears to be structurally similar, although these flaviviruses are transmitted to and from their vertebrate hosts by different carriers, such as ticks or mosquitoes. Analysis of sequence homologies of polyproteins of flaviviruses suggests that TBEV is more closely related to yellow fever virus than to other serological subgroups of flaviviruses (West Nile or Dengue viruses). The hydrophobic profiles of the flaviviruses are highly conservative. Nonstructural proteins NS2A, NS2B, NS4A and NS4B are extremely hydrophobic, suggesting that they are likely to be associated with cellular membranes. Proteins E, NS1, NS3 and NS5 are the most conservative and may be involved in general enzymatic activities related to viral replication and virion assembly.     

Characterization of a siberian virus isolated from a patient with progressive chronic tick-borne encephalitis

A strain of Tick-borne encephalitis virus designated Zausaev (Za) was isolated in Siberia from a patient who died of a progressive (2-year) form of tick-borne encephalitis 10 years after being bitten by a tick. The complete genomic sequence of this virus was determined, and an attempt was made to correlate the sequence with the biological characteristics of the virus. Phylogenetic analysis demonstrated that this virus belongs to the Siberian subtype of Tick-borne encephalitis virus. Comparison of Za virus with two related viruses, a Far Eastern isolate, Sofjin, and a Siberian isolate, Vasilchenko, revealed differences among the three viruses in pathogenicity for Syrian hamsters, cytopathogenicity for PS cells, plaque morphology, and the electrophoretic profiles of virus-specific nonstructural proteins. Comparative amino acid alignments revealed 10 individual amino acid substitutions in the Za virus polyprotein sequence that were different from those of other tick-borne flaviviruses. Notably, the dimeric form of the Za virus NS1 protein migrated in polyacrylamide gels as a heterogeneous group of molecules with a significantly higher electrophoretic mobility than those of the Sofjin and Vasilchenko viruses. Two amino acid substitutions, T(277)-->V and E(279)-->G, within the NS1 dimerization domain are probably responsible for the altered oligomerization of Za virus NS1. These studies suggest that the patient from whom Za virus was isolated died due to increased pathogenicity of the latent virus following spontaneous mutagenesis.     

A Sliding Window-Based Method to Detect Selective Constraints in Protein-Coding Genes and Its Application to RNA Viruses

Here we present a new sliding window-based method specially designed to detect selective constraints in specific regions of a multiple protein-coding sequence alignment. In contrast to previous window-based procedures, our method is based on a nonarbitrary statistical approach to find the appropriate codon-window size to test deviations of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions from the expectation. The probabilities of dN and dS are obtained from simulated data and used to detect significant deviations of dN and dS in a specific window region of the real sequence alignment. The nonsynonymous-to-synonymous rate ratio (w = dN/dS) was used to highlight selective constraints in any window wherein dS or dN was significantly different from the expectation. In these significant windows, w and its variance [V(w)] were calculated and used to test the neutral hypothesis. Computer simulations showed that the method is accurate even for highly divergent sequences. The main advantages of the new method are that it (i) uses a statistically appropriate window size to detect different selective patterns, (ii) is computationally less intensive than maximum likelihood methods, and (iii) detects saturation of synonymous sites, which can give deviations from neutrality. Hence, it allows the analysis of highly divergent sequences and the test of different alternative hypothesis as well. The application of the method to different human immunodeficiency virus type 1 and to foot-and-mouth disease virus genes confirms the action of positive selection on previously described regions as well as on new regions.     

The caddisfly fauna (Trichoptera) of Russia: on the centenary of exploration

A brief overview of the history of faunistic studies in Russia during the past century is presented. The existence of 643 species of 148 genera (28 families) is reported according to the results of the long-term investigation of the Russian fauna of the Trichoptera. There are 4 regional faunas of Trichoptera in Russia: Boreal European, Caucasian, Siberian, and Far Eastern. The tasks and prospects of the future studies are considered.     

Lineage‐specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage‐specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage‐specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage‐specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.     

Differential rates of evolution for the ZFY-related zinc finger genes,Zfy, Zfx,and Zfa in the mouse genus Mus

A comparative study of the last exon of the zinc finger genes Zfx, Zfy, and Zfa from species of mice in the genus Mus was conducted to assess the extent of gene-specific and chromosome-specific effects on the evolutionary patterns among related X-, Y-, and autosomal-linked genes. Phylogenetic analyses of 29 sequences from Zfx, Zfa, and Zfy from 10 taxa were performed to infer relatedness among the zinc finger loci, and codon-based maximum likelihood analyses were conducted to assess evolutionary pattern among genes. Five models of nucleotide sequence evolution were applied and compared using a likelihood ratio test. Estimates of nonsynonymous to synonymous changes (dN/dS) for these genes suggest that amino acid substitutions are occurring at a more rapid rate across the autosomal- and Y-specific lineages compared to the X-specific lineage, with the Y-specific lineage showing the highest rate under certain models. The data suggest the action of gene-specific effects on evolutionary pattern. In particular, Zfa and Zfy genes, both with presumed restricted expression, appear less functionally constrained relative to ubiquitously expressed Zfx. Slightly elevated dN/dS for Zfy genes in comparison to Zfa also suggest Y-specific effects.     

Uncorrected nucleotide bias in mtDNA can mimic the effects of positive Darwinian selection

The relative rates of nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions have been widely used to infer the action of natural selection from comparative sequence data. It is known, however, that mutational and repair biases can affect rates of evolution at both synonymous and nonsynonymous sites. More importantly, it is also known that synonymous sites are particularly prone to the effects of nucleotide bias. This means that nucleotide biases may affect the calculated ratio of substitution rates at synonymous and nonsynonymous sites. Using a large data set of animal mitochondrial sequences, we demonstrate that this is, in fact, the case. Highly biased nucleotide sequences are characterized by significantly elevated dN/dS ratios, but only when the nucleotide frequencies are not taken into account. When the analysis is repeated taking the nucleotide frequencies at each codon position into account, such elevated ratios disappear. These results suggest that the recently reported differences in dN/dS ratios between vertebrate and invertebrate mitochondrial sequences could be explained by variations in mitochondrial nucleotide frequencies rather than the effects of positive Darwinian selection.     

A tick-borne Langat virus mutant that is temperature sensitive and host range restricted in neuroblastoma cells and lacks neuroinvasiveness for immunodeficient mice

Langat virus (LGT), the naturally attenuated member of the tick-borne encephalitis virus (TBEV) complex, was tested extensively in clinical trials as a live TBEV vaccine and was found to induce a protective, durable immune response; however, it retained a low residual neuroinvasiveness in mice and humans. In order to ablate or reduce this property, LGT mutants that produced a small plaque size or temperature-sensitive (ts) phenotype in Vero cells were generated using 5-fluorouracil. One of these ts mutants, clone E5-104, exhibited a more than 10(3)-fold reduction in replication at the permissive temperature in both mouse and human neuroblastoma cells and lacked detectable neuroinvasiveness for highly sensitive immunodeficient mice. The E5-104 mutant possessed five amino acid substitutions in the structural protein E and one change in each of the nonstructural proteins NS3 and NS5. Using reverse genetics, we demonstrated that a Lys(46)-->Glu substitution in NS3 as well as a single Lys(315)-->Glu change in E significantly impaired the growth of LGT in neuroblastoma cells and reduced its peripheral neurovirulence for SCID mice. This study and our previous experience with chimeric flaviviruses indicated that a decrease in viral replication in neuroblastoma cells might serve as a predictor of in vivo attenuation of the neurotropic flaviviruses. The combination of seven mutations identified in the nonneuroinvasive E5-104 mutant provided a useful foundation for further development of a live attenuated TBEV vaccine. An evaluation of the complete sequence of virus recovered from brain of SCID mice inoculated with LGT mutants identified sites in the LGT genome that promoted neurovirulence/neuroinvasiveness.     

Replication enhancer elements within the open reading frame of tick-borne encephalitis virus and their evolution within the Flavivirus genus

We provide experimental evidence of a replication enhancer element (REE) within the capsid gene of tick-borne encephalitis virus (TBEV, genus Flavivirus). Thermodynamic and phylogenetic analyses predicted that the REE folds as a long stable stem-loop (designated SL6), conserved among all tick-borne flaviviruses (TBFV). Homologous sequences and potential base pairing were found in the corresponding regions of mosquito-borne flaviviruses, but not in more genetically distant flaviviruses. To investigate the role of SL6, nucleotide substitutions were introduced which changed a conserved hexanucleotide motif, the conformation of the terminal loop and the base-paired dsRNA stacking. Substitutions were made within a TBEV reverse genetic system and recovered mutants were compared for plaque morphology, single-step replication kinetics and cytopathic effect. The greatest phenotypic changes were observed in mutants with a destabilized stem. Point mutations in the conserved hexanucleotide motif of the terminal loop caused moderate virus attenuation. However, all mutants eventually reached the titre of wild-type virus late post-infection. Thus, although not essential for growth in tissue culture, the SL6 REE acts to up-regulate virus replication. We hypothesize that this modulatory role may be important for TBEV survival in nature, where the virus circulates by non-viraemic transmission between infected and non-infected ticks, during co-feeding on local rodents.     

Apparent Selection Intensity for the Cytochrome Oxidase Subunit I Gene Varies with Mode of Reproduction in Echinoderms

When most amino acid substitutions in protein-coding genes are slightly deleterious rather than selectively neutral, life history differences can potentially modify the effective population size or the selective regime, resulting in altered ratios of non-synonymous to synonymous substitutions among taxa. We studied substitution patterns for the mitochondrial cytochrome oxidase subunit I (COI) gene in a sea star genus (Leptasterias spp.) with an obligate brood-protecting mode of reproduction and small-scale population genetic subdivision, and compared the results to available COI sequences in nine other genera of echinoderms with pelagic larvae: three sea stars, five sea urchins and one brittle star. We predicted that this life history difference would be associated with differences in the ratio of non-synonymous (dN) to synonymous (dS) substitution rates. Leptasterias had a significantly greater dN/dS ratio (both between species and within species), a significantly smaller transition/transversion rate ratio, and a significantly lower average nucleotide diversity within species, than did the non-brooding genera. Other explanations for the results, such as altered mutation rates or selective sweeps, were not supported by the data analysis. These findings highlight the potential influence of reproductive traits and other life history factors on patterns of nucleotide substitution within and between species.     

Non-Hemagglutinating Flaviviruses: Molecular Mechanisms for the Emergence of New Strains via Adaptation to European Ticks

Tick-borne encephalitis virus (TBEV) causes human epidemics across Eurasia. Clinical manifestations range from inapparent infections and fevers to fatal encephalitis but the factors that determine disease severity are currently undefined. TBEV is characteristically a hemagglutinating (HA) virus; the ability to agglutinate erythrocytes tentatively reflects virion receptor/fusion activity. However, for the past few years many atypical HA-deficient strains have been isolated from patients and also from the natural European host tick, Ixodes persulcatus. By analysing the sequences of HA-deficient strains we have identified 3 unique amino acid substitutions (D67G, E122G or D277A) in the envelope protein, each of which increases the net charge and hydrophobicity of the virion surface. Therefore, we genetically engineered virus mutants each containing one of these 3 substitutions; they all exhibited HA-deficiency. Unexpectedly, each genetically modified non-HA virus demonstrated increased TBEV reproduction in feeding Ixodes ricinus, not the recognised tick host for these strains. Moreover, virus transmission efficiency between infected and uninfected ticks co-feeding on mice was also intensified by each substitution. Retrospectively, the mutation D67G was identified in viruses isolated from patients with encephalitis. We propose that the emergence of atypical Siberian HA-deficient TBEV strains in Europe is linked to their molecular adaptation to local ticks. This process appears to be driven by the selection of single mutations that change the virion surface thus enhancing receptor/fusion function essential for TBEV entry into the unfamiliar tick species. As the consequence of this adaptive mutagenesis, some of these mutations also appear to enhance the ability of TBEV to cross the human blood-brain barrier, a likely explanation for fatal encephalitis. Future research will reveal if these emerging Siberian TBEV strains continue to disperse westwards across Europe by adaptation to the indigenous tick species and if they are associated with severe forms of TBE.