|Home | About | Journals | Submit | Contact Us | Français|
Theiler's murine encephalitis viruses (TMEV) are divided into two subgroups based on their neurovirulence. Persistent strains resemble Theiler's original viruses (referred to as the TO subgroup), which largely induce a subclinical polioencephalomyelitis during the acute phase of the disease and can persist in the spinal cord of susceptible animals, inducing a chronic demyelinating disease. In contrast, members of the neurovirulent subgroup cause an acute encephalitis characterized by the rapid onset of paralysis and death within days following intracranial inoculation. We report herein the characterization of a novel neurovirulent strain of TMEV, identified using pyrosequencing technology and referred to as NIHE. Complete coverage of the NIHE viral genome was obtained, and it shares <90% nucleotide sequence identity to known TMEV strains irrespective of subgroup, with the greatest sequence variability being observed in genes encoding the leader and capsid proteins. The histopathological analysis of infected brain and spinal cord demonstrate inflammatory lesions and neuronal necrosis during acute infection with no evidence of viral persistence or chronic disease. Intriguingly, genetic analysis indicates the putative expression of the L* protein, considered a hallmark of strains within the persistent subgroup. Thus, the identification and characterization of a novel neurovirulent TMEV strain sharing features previously associated with both subgroups will lead to a deeper understanding of the evolution of TMEV strains and new insights into the determinants of neurovirulence.
Theilioviruses have positive-sense, single-stranded RNA genomes and are members of the Picornaviridae family and the Cardiovirus genus. Several serotypes have been identified, causing disease in rodents, such as Theiler's murine encephalomyelitis virus (TMEV) (36, 37) and rat Theiler-like virus (RTV) (28), and humans, such as Vilyuisk human encephalomyelititis virus (VHEV) (6) and Saffold viruses (7, 13, 44). TMEV is an enteric pathogen that follows the fecal-oral route of infection; however, when inoculated intracranialy (IC), it can cause severe pathology in susceptible mouse strains. As a natural mouse pathogen, TMEV provides a valuable experimental model for understanding host-pathogen interactions. The TMEV serotype is divided into two subgroups based on the pathogenesis of the disease. Members of the TO subgroup, so named because they resemble Theiler's original strains (36), cause a biphasic infection. The early phase is an acute poliomyelitis characterized by viral replication in the neurons of the brain (1, 18). During the late phase of infection, the virus progresses from the gray to the white matter of the spinal cord, where, depending on the genetics of the mouse strain (4), it can persist for the life of the animal within various cell types (2, 21), causing in some instances chronic demyelination (19). In contrast, the second subgroup, consisting of the GDVII and FA strains, are neurovirulent, and their pathology is characterized by an acute and rapidly fatal encephalitis (37). These viruses replicate mainly in the neurons of the brain, and for the few animals that survive primary infection, there does not appear to be viral persistence in the central nervous system (CNS) (18). Both neurovirulent strains were isolated in 1937 by Theiler and Gard; however, full-length sequences and infectious clones have been generated only for the GDVII strain. For the FA strain, only partial genomic sequence covering the 5′-untranslated region (5′UTR) through the end of the P1 region are available (GenBank accession no. U32924 and M80883). Within these 3,262 bp, there is 99.7% identity to the GDVII strain, which has led to speculation that GDVII and FA actually are separate isolates of the same strain due to the time and place of isolation (27).
Here, we report the phenotype and genetic analysis of the complete genome of a novel neurovirulent strain of TMEV that we have designated NIHE based on its isolation at the National Institute of Hygiene and Epidemiology in Hanoi, Vietnam. We demonstrate that the kinetics of disease and the phenotype of infected animals mimic those published for GDVII and FA strains. Moreover, histopathological analyses demonstrate that tissue and cellular tropism in the CNS are consistent with those observed following the inoculation of neurovirulent strains, including the inability to persist past the acute phase in the small percentage of inoculated animals that do survive. Phylogenetic analysis reveals NIHE to be less divergent from the TO strains than GDVII. Interestingly, we present evidence for the presence of an alternate open reading frame (ORF) which is predicted to encode the full-length L* protein, facilitated by a canonical AUG start sequence and previously shown to be present only in the TO subgroup strains (15). These data illustrate the need for continued investigation to establish the genetic and host immune determinants of neurovirulence and the evolution of viral persistence in this family of viruses.
All animal experiments described in the current study were conducted at the Institut Pasteur according to the European Union guidelines for the handling of laboratory animals The protocol was approved by the Institut Pasteur animal care and use committee. Inoculations were performed under anesthesia, and all efforts were made to minimize animal suffering.
Outbred Swiss mice were bred and maintained in the animal facility of the National Institute of Hygiene and Epidemiology in Hanoi, Vietnam. This colony originated from animals imported from France. The initial phenotype was observed in 2008 following the low-number serial passage of brain homogenate from apparently healthy donors. Samples were transferred to the Institut Pasteur in Paris for further characterization. Female 5- to 7-week-old SJL and C57BL/6 mice were purchased from Charles River Laboratories and housed in our SPF animal facility at the Institut Pasteur in Paris. The handling of mice and experimental procedures were performed in accordance with institutional guidelines for animal care. Inoculations were performed in mice via the intracranial injection of 30 to 35 μl of mouse brain homogenates. Specifically, brains from infected animals were prepared as 10% homogenates by mechanical disruption in sterile serum-free Dulbecco's modified essential medium (DMEM). Homogenates were clarified by centrifugation (10,000 × g for 20 min). Control animals received 30 to 35 μl of clarified supernatants from naïve animals at the indicated dilutions. All dilutions were prepared in sterile phosphate-buffered saline (PBS) or DMEM.
Clarified supernatants were generated from the brains of animals demonstrating a disease phenotype. Aliquots of supernatant were combined with the AVL buffer supplied in the QIAamp viral RNA Minikit (Qiagen) and processed according to the manufacturer's specifications. Aliquots were combined, and 454 sequencing was performed by Eurofins MWG Operon (Ebersberg, Germany).
FASTA files containing individual sequencing reads were obtained from the sequencing reactions. The nonredundant protein (nr) and nucleotide (nt) sequence databases were downloaded from NCBI and updated August 2010 prior to query. All sequences were queried against the nr/nt databases using BLAST+ (5). Metagenomic analysis was performed using MEGAN software (version 3.9), which resulted in the grouping of sequences according to phylogeny (11). The de novo assembly of sequences greater than 80 bp was performed using the local alignment algorithm within CodonCode Aligner software (version 3.5), and the resulting consensus sequence was checked for inserted/deleted sequences that can be an artifact of 454 sequencing homopolymeric sequences. Twenty-three sites were manually confirmed and corrected throughout the sequence. The corrected consensus sequence was used in all subsequent comparative analyses.
Infected or control animals were anesthetized and perfused with PBS followed by 4% paraformaldehyde in PBS. Brains and spinal cords were dissected and fixed in 4% paraformaldehyde overnight at 4°C. Tissues were embedded in paraffin by following established procedures, and 5-μm serial sections from the brain and spinal cord were stained with hematoxylin-eosin (HE). Cross-reactivity with TMEV antigens was tested by immunohistochemistry on paraffin sections incubated with a primary rabbit polyclonal antibody raised against Theiler's virus capsid proteins (kindly provided by J. F. Bureau and M. Brahic). Slides were counterstained with hematoxylin.
Nucleotide and associated protein sequences for the following viruses were obtained from GenBank (data concerning three domains were collected: protein sequences coding for the 3D polymerase [3Dpol], P1 [protein 1], and VP1 [viral protein 1]): GDVII (M20562 and NC_001366), FA (U32924 and M80883 combined), DA (M20301), Yale (EU723238), BeAn8386 (M16020), TOB15 (EU718732), VIE415HTR (EU718733), TO4 (U33045 and M80885 combined), WW (M80889 and U33046 combined), Mengo virus (DQ294633), encephalomyocarditis virus (NC_001479), porcine encephalomyocarditis (DQ517424), Saffold virus (FM207487, NC_009448, and EF165067), cardiovirus D (EU681176), human TMEV-like cardiovirus (NC_010810 and EU376394), RTV NGS910 (AB090161), RTV-1 (EU542581 and EU815052), and Vilyuisk (EU723237, M80888, and M94868 combined). Multiple-sequence alignment was performed using the Muscle algorithm implemented within SeaView (version 4.2.4). The 3Dpol data set included 20 taxa and 1,404 characters, the P1 data set had 23 taxa for 2,617 nucleotides, and the VP1 matrix involved 23 taxa for 868 characters.
The phylogenetic reconstructions of three protein-coding sequences (3Dpol, VP1, and V1) from three sets of nucleotides were inferred separately by a maximum-likelihood analysis on matrixes composed of 20, 20, and 23 taxa, respectively. These phylogenetic analyses were implemented by the BEAST software package (version 1.6.1). As this method required a predetermined evolution model estimated on the data set, we specified as one of the models selected by Jmodeltest (30) on each matrix (3Dpol, VP1, and P1). Using the results of Jmodeltest, the three-codon partition model of evolution was used (under the GTR substitution model with a gamma site heterogeneity distribution [GTR+G]), so that each codon position (first, second, and third) could follow its own evolution rate. The Yule speciation process (specified as tree prior) and constant-population-size coalescent models show similar estimates for each data set. The Monte Carlo Markov chain was set to 5,000,000 states for 3Dpol and VP1 data sets and to 10,000,000 for the P1 matrix to reach adequate prior effective sampling size (ESS) values greater than 300 and 2,000, respectively. Trees were annotated with TreeAnnotator (version 1.4.8) and displayed and edited using FigTree (version 1.3.1).
Similarity plots were generated using SimPlot (version 3.5.1) with the following parameters: 200-bp window, 20-bp step, and Kimura 2-parameter distance model (22).
The nucleotide sequence reported in this paper was submitted to GenBank and has been assigned the accession no. HQ652539.
While performing animal experiments involving the serial passage of clarified lysate from mouse brain homogenates, we observed the appearance of a distinct disease phenotype. Following inoculation, several animals displayed a flaccid tail and an abnormal gait, which progressed to unilateral or bilateral limb paralysis 9 to 10 days postinoculation. Following the next passage, 100% of inoculated animals succumbed to infection; continued serial passage has demonstrated consistent lethality. To determine the kinetics of death, groups of C57BL/6 mice (n = 10) were inoculated IC with 10× dilutions (10−1 to 10−6) of infected mouse brain homogenate (IMBH), and survival was monitored (Fig. 1A). Animals receiving a 10−1 dilution had a mean survival time of 7 days, whereas animals receiving 10−2, 10−3, and 10−4 dilutions had mean survival times of 9, 10, and 11.5 days, respectively. The 50% lethal dose (LD50) was determined to be a 10−5 dilution of IMBH. Control animals received a 10−1 dilution of 10% brain homogenate derived from naïve animals with no apparent phenotype (>45 days postinoculation). To quantitatively describe the disease progression in mice, a clinical scoring method was established (Fig. 1B), and disease progression was monitored daily (Fig. 1C). To accurately describe disease progression, only symptomatic animals were included; the number of symptomatic animals per group is indicated in the figure legend.
To determine the identity of the infectious agent, brains from symptomatic mice were pooled and homogenized, and viral RNA was enriched for sequencing. Mass sequencing was performed using 454 GS-FLX technology (Roche). Approximately 700,000 sequencing reads were obtained with an average read length of 330 nucleotides. BLAST analyses were performed, and mouse genomic sequences were excluded. The remaining sequences were blasted against available viral genomes in the NCBI nr/nt databases. Of the identified sequences, 6.9% (~5,100/73,000) showed similarity to members of the Cardiovirus genus, with the greatest sequence identity to TMEV isolates. Sequences were assembled de novo into contigs using a local alignment algorithm (CodonCode Aligner, version 3.5.5). Complete and redundant genomic coverage of a virus (designated NIHE) was obtained with an average of 177 reads per 100-bp segment (Fig. 2). Based on these analyses, a consensus sequence was generated and translated to reveal a full-length polyprotein. The nucleotide and amino acid sequences of NIHE were used to compare percent identities to those of full-length neurovirulent (GDVII) and persistent (BeAn, DA, Yale, TOB15, and Vie415HTR) TMEV strains available in GenBank. Strains with partial sequences were also considered, namely, neurovirulent (FA and ASK-1) and persistent (TO4 and WW) TMEV strains (Tables (Tables11 and and22).
To further characterize the pathogenesis of NIHE-mediated disease, histopathological analysis was performed on brains and spinal cords isolated from SJL and C57BL/6 mice 8 days postinoculation with 10−1 to 10−7 dilutions of IMBH. In the brain, acute inflammatory lesions were observed in several regions, including the cerebral cortex, thalamus, brain stem, and cerebellum, although most of the neuronal destruction and viral antigen was observed in the hippocampus (Fig. 3). In the gray matter, perivascular cuffing was observed and was comprised mainly of mononuclear cells. Many neurons showed characteristics of necrosis and were surrounded by mononuclear cells showing evidence of neuronophagia. Examples of the histopathological findings are shown for the hippocampus (Fig. 3A to C). Infected regions of the spinal cord showed similar, albeit more severe, inflammatory lesions. Perivascular cuffs and neuronal necrosis were observed primarily in the gray matter, but perivascular infiltrates very occasionally were observed in the white matter as well (Fig. 3D). To determine if NIHE is cross-reactive to antibodies known to recognize other TMEV strains (31), immunohistochemistry was performed and representative sections are shown, demonstrating significant viral infection in the brain throughout the hippocampus and thalamus (Fig. 3E and F) as well as gray matter of the spinal cord (data not shown). Numerous foci of infected cells, particularly neurons, are evident in both tissues. Severe neuronal necrosis and TMEV staining were observed, particularly in the CA1 and CA2 pyramidal cell layers of the hippocampus, yet the CA3 and dentate gyrus areas were almost completely spared (Fig. 3A and B).
A shift in the tissue tropism from the gray matter to the white matter of the spinal cord is one of the defining characteristics of the persistent TMEV subgroup. The capacity of NIHE to persist in susceptible SJL animals was examined by the histological analysis of the brain and spinal cord 45 days postinoculation, the time point generally considered the chronic phase of disease (25). Animals were inoculated IC with increasing dilutions of infected mouse brain homogenate. All animals that received a dilution of >10−4 succumbed during the acute phase of disease, 18% of the animals that received a 10−5 dilution died, and all animals that received a ≤10−6 dilution remained asymptomatic (data not shown). Histopathological analysis was performed on brain and spinal cord, and no evidence of inflammatory lesions was present in the tissues, and viral antigens were not detected following anti-TMEV staining (data not shown).
To more fully understand the basis for the neurovirulence of NIHE, as well as to determine its relationship to other TMEV group members, we performed a detailed genetic analysis. All strains with full-length sequence data were aligned and compared to the sequence of NIHE (Fig. 4A). The greatest percent identity is seen in the 5′UTR. Interestingly, the persistent strain DA maintains the highest identity to NIHE through the sequence coding for the L protein (98% identity of amino acid sequence), but it diverges significantly for other genes (Fig. 4A, green line). Regions that demonstrated the lowest percentage of nucleotide identity between NIHE and other TMEV strains were clustered in the L protein and P1 region, particularly the sequences encoding the capsid proteins VP2 and VP1 (Fig. 4A to C). A more focused comparison of the P1 region allows the inclusion of partially sequenced strains and demonstrates additional dissimilarity, specifically within the VP2 and VP1 regions, with no apparent clustering with either subgroup (Fig. 4B). Interestingly, when the same analysis was performed using the GDVII strain as the reference sequence, a separation within the TO subgroup was observed (Fig. 4C). Specifically, the Yale and TOB15 strains (blue and black lines, respectively) maintain a relatively constant degree of sequence identity throughout this region (~90%). This is in contrast to the cluster of TO4/WW, BeAn, DA, and Vie415HTR strains, which differ from GDVII considerably (70 to 80%). NIHE (red line) appears at an intermediate position between these two groups.
To place genetic differences within NIHE in context with other known cardioviruses, a phylogenetic analysis of three proteins or regions was performed: (i) the well-conserved 3Dpol (RNA-dependent RNA polymerase); (ii) the P1 region, which encodes the four capsid proteins; and (iii) the highly variable VP1 gene (Fig. 5). The phylogenetic analysis of the 3Dpol supports the separation of the human viruses, Thera viruses (rat), and TMEV species. NIHE clusters within the TMEV clade; however, it segregates with Vilyuisk and the mildly neurovirulent Vie415HTR strain, while GDVII groups with the TO strains. These groupings are well supported based on the high posterior probability for the nodes (Fig. 5A). Differences in the topography of the P1 region and that of 3Dpol are seen. NIHE segregates with ASK-1 and is more closely related to GDVII and FA than to the TO strains (Fig. 5B). The phylogenetic relationships established for the P1 region are echoed in the VP1 analysis, with NIHE segregating with a group comprised mostly of neurovirulent strains (Fig. 5C).
To identify possible genetic determinants of neurovirulence and persistence in the NIHE genome, we aligned and compared the amino acid sequence for all coding regions of NIHE to the other complete and partial sequences of TMEV strains. The mean percent identity (compared to NIHE) of the amino acid sequence for the indicated viral protein is shown (Fig. 6). Although significant dissimilarity is seen in the VP2 and VP1 proteins, as predicted by the similarity plots (Fig. 4), the lowest percent similarity is seen in the L (88.8%) and L* (87.3%) proteins. A comparative analysis of the amino acid sequence for the carboxy terminus of the L protein, including several distinct domains, is shown in Fig. 6B.
Using modeling software (Swiss-PDB viewer  and RasWin ), we compared the predicted structure of the NIHE VP1-VP3 protomer complex to the structures of BeAn and DA strains (Fig. 7). The calculated structure of NIHE VP1, VP2, and VP3 proteins (blue, green, and red, respectively) was overlaid on the BeAn structure (Fig. 7A) (23). Differences in the amino acid sequence between these two strains are indicated by white-space-filled molecules (40 amino acid changes), and sequences that are unique to NIHE (compared to all published strains) are shown in yellow (14 amino acid changes). The spatial positioning of CD loop I and CD loop II are shown (Fig. 7A, upper right), with the loops delineated by the light-yellow and purple-filled molecules, respectively. Puff A (yellow) and puff B (light-blue) regions of VP2, as well as a protrusion from VP3 known as the knob, are highlighted (Fig. 7A, lower right). To further understand the relationship between NIHE and the other TMEV strains, specifically within these critical loop structures, we compared the amino acid sequences of selected regions (Fig. 7B).
The region that was observed to have the greatest intrastrain variation coded for the L* protein (Fig. 6A). The generation of L* is the result of an additional open reading frame initiated 13 nucleotides downstream of the AUG start site that is responsible for initiating the viral polyprotein (15). This protein is not efficiently produced by the GDVII strain (and, presumably, FA) due to a polymorphism (ACG) in the sequence that initiates translation. Figure 8A depicts the sequence of this region, including the polyprotein start site and both alternative initiation sites located immediately downstream. The putative sequences coding for the L* protein of available strains were aligned, and a midrooted tree was generated (Fig. 8B). An intrastrain alignment was performed, and the difference between putative coding sequences of NIHE and other strains are shown (Fig. 8C).
Our study reports the identification of a novel neurovirulent TMEV strain, the first since Theiler's characterization of the GDVII and FA strains in 1938. In fact, based on the time and location at which these two strains were discovered, as well as the current understanding of quasispecies diversity, it has been proposed that FA is a derivative of GDVII, thus making NIHE only the second fully annotated neurovirulent TMEV. A strain isolated in Japan, referred to as ASK-1 (6), has been proposed to be a member of the neurovirulent subgroup (3, 20); however, no characterization of this strain has been published, and only a partial sequence has been submitted to GenBank (24).
Several characteristics of the NIHE strain led us to classify it as a member of the neurovirulent subgroup: (i) the rapid onset of disease following IC inoculation (Fig. 1); (ii) marked neuronal destruction in the brain and spinal cord, far greater than that observed for persistent strains (Fig. 3); and (iii) a lack of evidence (pathology or viral antigen) that NIHE persists in the CNS of mice surviving low-dose inoculation (data not shown). We do, however, realize the inherent paradox in studying highly virulent TMEV strains, in that it is difficult to distinguish between the inability of the virus to persist and the failure of a productive infection due to low viral inoculation.
The use of next-generation sequencing (NGS) techniques allowed for the complete, redundant sequencing as well as the de novo assembly of the entire NIHE genome (Fig. 2), suggesting the identification of a novel TMEV strain. The degree of identity shared by NIHE and known TMEV strains is <90% at the nucleotide level and <95% in its primary amino acid sequence (Tables (Tables11 and and2).2). This is in contrast to the >99% identity for quasispecies variation seen for the TO4, WW, and DA strains (27) or the laboratory derivatives of S2 BeAn (14) or GDVII substrains (16, 29).
The discovery of NIHE provides new information about the determinants of neurovirulence versus persistence and offers a novel basis of comparison on which to examine the relevance of amino acid differences previously identified between the GDVII and TO strains. Intriguingly, our analysis suggests that although NIHE has a decidedly neurovirulent phenotype, it shares high sequence identity with the TO strains in many important regions, including L, L*, and capsid proteins.
The analysis of the sequence comparison between NIHE and the other TMEV sequences demonstrates clusters of differences in the amino acid sequence (9, 12, 17, 26, 32, 40, 41, 45). Although we present clear evidence that NIHE is neurovirulent, we were surprised to observe the number of instances in which the NIHE sequence segregated with the persistent strains. Here, we restrict our discussion to instances in which NIHE segregates with TO strains within regions important for modulating the host antiviral response, viral persistence, or interaction with the coreceptor sialic acid.
The L protein is produced by all members of the Cardiovirus genus and has been implicated in in vivo persistence because it modulates the host antiviral response (8, 34, 39). Although analysis revealed that the L protein is the second most divergent region across TMEV strains (Fig. 6), there were no unique amino acids predicted by the NIHE sequence. The majority of the sequence variation between NIHE and the other strains was clustered within the Ser/Thr-rich domain (62.5%) and the Theilo domain (79.5%). Interestingly, two positions in the Ser/Thr-rich region (L-55 and L-57) and one in the L-VP4 cleavage site (L-70) have been identified previously as distinguishing between the GDVII and TO strains (27, 35). In all three instances, the NIHE sequence segregates with the TO strains.
A slightly different result emerged when analyzing the L* reading frame. Several nucleotide changes are unique to the NIHE strain, resulting in distinct amino acid sequences at positions L*-123 (Ser), L*-147 (Arg), and L*-153 (Phe) compared to sequences of any of the known TMEV strains (Fig. 8C). Again, the strains encoding L* sequences most closely related to NIHE are those produced by the Yale, BeAn, and TOB15 strains. Perhaps the most notable example of this divergence is the observation that NIHE, like the TO strains, retains the genetic capacity to efficiently express L* by encoding an AUG sequence at the two initiation positions (Fig. 8). It is important to note that some L* transcripts can be produced via the ACG initiation sites, albeit less efficiently (38). What evolutionary advantage this protein provides to a neurovirulent strain is unclear, particularly given that it is not required to infect neurons during the acute phase of disease. Advantages, however, can be imagined if moderate levels of L* are required by the virus during the natural route of transmission or for the crossing of the blood-brain barrier.
Many intraspecies sequence differences mapped to exposed areas of the viral capsid (Fig. 7). Cluster A, a highly variable region within VP1, initially was believed to segregate between GDVII and persistent strains; however, sequence information for several viruses was not then available. The inclusion of recently published strains and our contribution of the NIHE genome further clarifies these observations. Although GDVII and FA segregate with ASK-1 and NIHE at two positions (VP1-139 and VP1-232), NIHE segregates with a subset of TO strains at positions VP1-136, VP1-142, and VP1-180. In addition to these differences, NIHE also encodes two unique amino acids within this domain (VP1-138 and VP1-140), as shown in Fig. 7.
Perhaps the most functionally relevant example of NIHE sequences similar to those of persistent TO strains is observed in VP2 puff B (within the EF loop). Positions VP2-161, VP2-162, VP2-163, and VP2-174 were demonstrated to be the binding site for sialic acid in the persistent strains (43). Although a cellular receptor for TMEV strains currently is unknown, this region has been considered relevant to persistence due to the observation that Ser at position VP2-162 within the GDVII/FA strains prevents VP1-161 (Glu) from forming a hydrogen bond with sialic acid. Alignment with NIHE and ASK-1 indicate that these nonpersistent strains are identical to the TO strains at these positions. VP2-171, VP2-172, and VP2-173 also influence binding to sialic acid by establishing a structural framework for Gly at VP2-174 to bond with sialic acid (24, 42, 43). NIHE is unique in this area, encoding Ser-Gln-Thr, whereas Arg-Gln-Ala is encoded in the GDVII/FA strains and Ser-Arg-Thr is in the majority of the TO strains (Fig. 7B). Interestingly, NIHE is identical to ASK-1, but the consequence on the orientation of Gly at VP2-174 remains unknown. Within this domain, there also exist three unique amino acids within NIHE (VP2-156, VP2-158, and VP2-177). None of these residues has been implicated previously in TMEV pathogenicity, and they may reflect either genetic drift in a nonessential position or compensatory mutations that promote viral entry or replication.
To ascertain information about the relationship of NIHE to other TMEV strains and gain insight into its evolution, we performed the phylogenetic analysis of highly conserved 3Dpol as well as the divergent regions P1 and VP1. The analysis of the 3Dpol gene clearly positions NIHE within the TMEV cluster and demonstrates that NIHE is more closely related to the human Vilyuisk virus and VIE415HTR than the other murine neurovirulent strains (Fig. 5A). The position of NIHE within the topology of the P1 and VP1 phylogenies is well supported and further supports our conclusion of neurovirulence. Furthermore, the addition of the NIHE sequence clarifies the groupings by defining a neurovirulent clade with GDVII, FA, and ASK-1, particularly within the VP1 analysis (Fig. 5C). Evidence of parallel evolutionary paths are predicted for NIHE, ASK-1, and GDVII/FA strains, all of which have obtained neurovirulence, albeit via potentially distinct mutations (Fig. 5B and C). Although it would be difficult to decipher an evolutionary history with such few numbers of sequenced strains, it has been hypothesized that neurovirulent strains have evolved from persistent strains (27).
In summary, our study reports the discovery and characterization of a novel, highly neurovirulent TMEV, the first report of a new strain in this subgroup in more than 70 years. Critically, we characterized the pathogenesis of disease and obtained complete sequencing information without subjecting NIHE to conditions that promote genomic changes (e.g., in vitro expansion, which is known to introduce artifactual compensatory mutations). Perhaps the most significant advance of this work is that the discovery of the NIHE strain helps to illuminate the relevance of sequences putatively implicated in the evolution of neurovirulence versus persistence, data that previously had been based solely on comparison with the GDVII strain. NIHE also provides an opportunity to evaluate new positions as potential determinants for neurovirulence, as NIHE contains unique amino acids at several positions, particularly within the capsid proteins. Further investigation will be required to assess to what extent these positions contribute to the extreme phenotype of neurovirulence observed in the NIHE strain.
This work was supported by grants from the Agence National de la Recherche (ANR) and an EURYI Award from the European Science Foundation (ESF).
We are grateful to Michel Brahic (Stanford University, Stanford, CA) for the anti-TMEV antibody and to Thomas Michels (Université Catholique de Louvain, Brussels, Belgium) for the anti-L* antibody. We also thank Patric Avé (Histotechnologie et Pathologie, Institut Pasteur, Paris, France) for his expert technical assistance in the preparation of the histology slides and Herbert Virgin IV (Washington University School of Medicine, St. Louis, MO) for advice on deep sequencing and data mining. We thank Olivier Schwartz and Felix Rey for their critical reviews of the manuscript.
Published ahead of print on 4 May 2011.