Coexistence in Field Samples of Two Variants of the Infectious Salmon Anemia Virus: A Putative Shift to Pathogenicity

Genetic reassortment plays an important role in the evolution of several segmented RNA viruses and in the epidemiology of their associated diseases. In particular, orthomyxoviruses show rapid fluctuation in the proportion of viral variants coexisting in an infected individual, especially under strong selective pressure. This is particularly relevant in salmon production carried out under confined and stressful conditions where one of the most feared pathogenic agents is the Infectious Salmon Anemia Virus, an orthomyxovirus family member whose biological behavior is only recently beginning to be understood. Pathogenicity of the virus has been mainly associated with deletions of the HPR region in coding segment 6 and the presence or absence of a specific insertion in a key region in coding segment 5. In this study we report, for the first time in Chile, the coexistence of two variants in fully asymptomatic fish. Of five samples analyzed, two were identified as the non-pathogenic variant, HPR0, and two as the highly pathogenic HPR7b variant, though with no clinical signs detectable in the fish. Interestingly, one of the samples unequivocally carried both variants, again without any clinical signs. Considering that in none of the samples the typical insertion in coding segment 5 was detected, it is our impression that this may represent a shift from the non-pathogenic HPR0 variant towards the highly infective HPR7b variant. If this were the case, the transition may be triggered first by deleting the corresponding sequence of the HPR region of segment 6, followed by the putative insertion in segment 5 to generate a virulent strain.     

Novel strategy to evaluate infectious salmon anemia virus variants by high resolution melting

Genetic variability is a key problem in the prevention and therapy of RNA-based virus infections. Infectious Salmon Anemia virus (ISAv) is an RNA virus which aggressively attacks salmon producing farms worldwide and in particular in Chile. Just as with most of the Orthomyxovirus, ISAv displays high variability in its genome which is reflected by a wider infection potential, thus hampering management and prevention of the disease. Although a number of widely validated detection procedures exist, in this case there is a need of a more complex approach to the characterization of virus variability. We have adapted a procedure of High Resolution Melting (HRM) as a fine-tuning technique to fully differentiate viral variants detected in Chile and projected to other infective variants reported elsewhere. Out of the eight viral coding segments, the technique was adapted using natural Chilean variants for two of them, namely segments 5 and 6, recognized as virulence-associated factors. Our work demonstrates the versatility of the technique as well as its superior resolution capacity compared with standard techniques currently in use as key diagnostic tools.     

Functional role of glutamine 28 and arginine 39 in double stranded RNA cleavage by human pancreatic ribonuclease

Human pancreatic ribonuclease (HPR), a member of RNase A superfamily, has a high activity on double stranded (ds) RNA. By virtue of this activity HPR appears to be involved in the host-defense against pathogenic viruses. To delineate the mechanism of dsRNA cleavage by HPR, we have investigated the role of glutamine 28 and arginine 39 of HPR in its activity on dsRNA. A non-basic residue glycine 38, earlier shown to be important for dsRNA cleavage by HPR was also included in the study in the context of glutamine 28 and arginine 39. Nine variants of HPR respectively containing Q28A, Q28L, R39A, G38D, Q28A/R39A, Q28L/R39A, Q28A/G38D, R39A/G38D and Q28A/G38D/R39A mutations were generated and functionally characterized. The far-UV CD-spectral analysis revealed all variants, except R39A, to have structures similar to that of HPR. The catalytic activity of all HPR variants on single stranded RNA substrate was similar to that of HPR, whereas on dsRNA, the catalytic efficiency of all single residue variants, except for the Q28L, was significantly reduced. The dsRNA cleavage activity of R39A/G38D and Q28A/G38D/R39A variants was most drastically reduced to 4% of that of HPR. The variants having reduced dsRNA cleavage activity also had reduction in their dsDNA melting activity and thermal stability. Our results indicate that in HPR both glutamine 28 and arginine 39 are important for the cleavage of dsRNA. Although these residues are not directly involved in catalysis, both arginine 39 and glutamine 28 appear to be facilitating a productive substrate-enzyme interaction during the dsRNA cleavage by HPR.     

Identification of attenuating mutations on the reovirus type 3 S1 double-stranded RNA segment with a rapid sequencing technique.

Reovirus type 3 variants with mutations in the major neutralization domain of the sigma 1 protein have attenuated neurovirulence and restricted neurotropism. We devised a variation of the rapid RNA sequencing technique to facilitate the analysis of double-stranded RNA. We sequenced the S1 double-stranded RNA segment, which encodes the sigma 1 protein, of five attenuated reovirus type 3 variants. Four of the variants have changes in codon 419, and a fifth variant has a change at codon 340, all of which resulted in amino acid substitutions in the sigma 1 protein. We identified two sites on the reovirus type 3 sigma 1 protein that play a critical role in neurovirulence.     

Ultra-Deep Pyrosequencing of Partial Surface Protein Genes from Infectious Salmon Anaemia Virus (ISAV) Suggest Novel Mechanisms Involved in Transition to Virulence

Uncultivable HPR0 strains of infectious salmon anaemia viruses (ISAVs) infecting gills are non-virulent putative precursors of virulent ISAVs (vISAVs) causing systemic disease in farmed Atlantic salmon (Salmo salar). The transition to virulence involves two molecular events, a deletion in the highly polymorphic region (HPR) of the hemagglutinin-esterase (HE) gene and a Q₂₆₆→L₂₆₆ substitution or insertion next to the putative cleavage site (R₂₆₇) in the fusion protein (F). We have performed ultra-deep pyrosequencing (UDPS) of these gene regions from healthy fish positive for HPR0 virus carrying full-length HPR sampled in a screening program, and a vISAV strain from an ISA outbreak at the same farming site three weeks later, and compared the mutant spectra. As the UDPS data shows the presence of both HE genotypes at both sampling times, and the outbreak strain was unlikely to be directly related to the HPR0 strain, this is the first report of a double infection with HPR0s and vISAVs. For F amplicon reads, mutation frequencies generating L₂₆₆ codons in screening samples and Q₂₆₆ codons in outbreak samples were not higher than at any random site. We suggest quasispecies heterogeneity as well as RNA structural properties are linked to transition to virulence. More specifically, a mechanism where selected single point mutations in the full-length HPR alter the RNA structure facilitating single- or sequential deletions in this region is proposed. The data provides stronger support for the deletion hypothesis, as opposed to recombination, as the responsible mechanism for generating the sequence deletions in HE.     

Strain variation, based on the hemagglutinin gene, in Norwegian ISA virus isolates collected from 1987 to 2001: indications of recombination.

Infectious salmon anemia (ISA) is caused by a virus that probably belongs to the Orthomyxoviridae and was first recorded in Norway in 1984. The disease has since spread along the Norwegian coast and has later been found in Canada, Scotland, the Faroe Islands, Chile, and the USA. This study presents sequence variation of the hemagglutinin gene from 37 ISA virus isolates, viz. one isolate from Scotland, one from Canada and 35 from Norway. The hemagglutinin gene contains a highly polymorphic region (HPR), which together with the rest of the gene sequence provides a good tool for studies of epizootics. The gene shows temporal and geographical sequence variation, where certain areas are dominated by distinct groups of isolates. Evidence of transmission of ISA virus isolates within and between regions is given. It is suggested that the hemagglutinin gene from different isolates may recombine. Possible recombination sites are found within the HPR and in the 5'-end flanking region close to the HPR.     

Adaptation of Puumala hantavirus to cell culture is associated with point mutations in the coding region of the L segment and in the noncoding regions of the S segment

We previously developed a model for studies on hantavirus host adaptation and initiated genetic analysis of Puumala virus variants passaged in colonized bank voles and in cultured Vero E6 cells. With the data presented in this paper, the sequence comparison of the wild-type and Vero E6-adapted variants of Puumala virus, strain Kazan, has been completed. The only amino acid substitution that distinguished the two virus variants was found in the L protein, Ser versus Phe at position 2053. Another mutation found in the L segment, the silent transition C1053U, could result from the selection of a variant with altered L RNA folding. Nucleotide substitutions observed in individual L cDNA clones, most of them A-->G and U-->C transitions, suggested that the population of L RNA molecules is represented by quasispecies. The mutation frequency in the L segment quasispecies appeared to be similar to the corresponding values for the S and M quasispecies. Analysis of the cDNA clones with the complete S segment sequences from passage 20 confirmed our earlier conclusion that the cell-adapted genotype of the virus is represented mostly by variants with mutated S segment noncoding regions. However, the spectrum of the S segment quasispecies appeared to be changing, suggesting that, after the initial adaptation (passages 1 to 11), the viral population is still being driven by selection for variants with higher fitness.     

Sequential bottlenecks drive viral evolution in early acute hepatitis C virus infection

Hepatitis C is a pandemic human RNA virus, which commonly causes chronic infection and liver disease. The characterization of viral populations that successfully initiate infection, and also those that drive progression to chronicity is instrumental for understanding pathogenesis and vaccine design. A comprehensive and longitudinal analysis of the viral population was conducted in four subjects followed from very early acute infection to resolution of disease outcome. By means of next generation sequencing (NGS) and standard cloning/Sanger sequencing, genetic diversity and viral variants were quantified over the course of the infection at frequencies as low as 0.1%. Phylogenetic analysis of reassembled viral variants revealed acute infection was dominated by two sequential bottleneck events, irrespective of subsequent chronicity or clearance. The first bottleneck was associated with transmission, with one to two viral variants successfully establishing infection. The second occurred approximately 100 days post-infection, and was characterized by a decline in viral diversity. In the two subjects who developed chronic infection, this second bottleneck was followed by the emergence of a new viral population, which evolved from the founder variants via a selective sweep with fixation in a small number of mutated sites. The diversity at sites with non-synonymous mutation was higher in predicted cytotoxic T cell epitopes, suggesting immune-driven evolution. These results provide the first detailed analysis of early within-host evolution of HCV, indicating strong selective forces limit viral evolution in the acute phase of infection.     

Viral evolution and interferon resistance of hepatitis C virus RNA replication in a cell culture model

Hepatitis C virus (HCV) replicates through an error-prone process that may support the evolution of genetic variants resistant to the host cell antiviral response and interferon (IFN)-based therapy. We evaluated HCV-IFN interactions within a long-term culture system of Huh7 cell lines harboring different variants of an HCV type 1b subgenomic RNA replicon that differed at only two sites within the NS5A-encoding region. A replicon with a K insertion at HCV codon 2040 replicated efficiently and exhibited sequence stability in the absence of host antiviral pressure. In contrast, a replicon with an L2198S point mutation replicated poorly and triggered a cellular response characterized by IFN-beta production and low-level IFN-stimulated gene (ISG) expression. When maintained in long term-culture, the L2198S RNA evolved into a stable high-passage (HP) variant with six additional point mutations throughout the HCV protein-encoding region that enhanced viral replication. The HP RNA transduced Huh7 cells with more than 1,000-fold greater efficiency than its L2198S progenitor or the K2040 sequence. Replication of the HP RNA resisted suppression by IFN-alpha treatment and was associated with virus-directed reduction in host cell expression of ISG56, an antagonist of HCV RNA translation. Accordingly, the HP RNA was retained within polyribosome complexes in vivo that were refractory to IFN-induced disassembly. These results identify ISG56 as a translational control effector of the host response to HCV and provide direct evidence to link this response to viral sequence evolution, ISG regulation, and selection of the IFN-resistant viral phenotype.     

Evolution of the nucleotide sequence of influenza virus RNA segment 7 during drift of the H3N2 subtype

The complete genetic information contained in the influenza virus RNA segment 7 of the A/Bangkok/ $\text{1}\text{79}$ (H3N2) strain has been cloned by in vitro synthesis of the complementary dsDNA and its insertion into plasmid pBR322. The nucleotide sequence of the viral RNA segment has been determined from the cDNA insert. It is 1027 nucleotides long, and contains two open reading frames, as shown for other influenza virus strains. When compared with the previously published sequence for the A/Udorn/72 (H3N2) strain, 15 nucleotide exchanges are observed, most of them silent mutations, and only two causing amino acid changes in each of the M1 and M2 protein sequences.     

Molecular basis of viral persistence: a single amino acid change in the glycoprotein of lymphocytic choriomeningitis virus is associated with suppression of the antiviral cytotoxic T-lymphocyte response and establishment of persistence.

Isolates of lymphocytic choriomeningitis virus (LCMV) that elicit a cytotoxic T-lymphocyte response (CTL+) have been compared with isolates that suppress the CTL response (CTL-) in an effort to map this phenotype. A single amino acid change in the glycoprotein of the LCMV Armstrong (ARM) strain is consistently associated with the CTL- trait and the ability of the virus to persist (P+). The CTL+ P- parental strain spontaneously gives rise to CTL- P+ variants within lymphoid tissues of mice persistently infected from birth. To map the structural basis of the phenotype, the complete RNA sequence of LCMV ARM 53b (CTL+) was compared with that of its variant ARM clone 13 (CTL-). Differences in 5 of 10,600 nucleotides were found. Three changes are noted in the large L RNA segment, and two are noted in the small S RNA segment. Only two of the changes distinguishing CTL+ from CTL- isolates affect amino acid coding: lysine to glutamine at amino acid 1079 of the polymerase protein, and phenylalanine to leucine at amino acid 260 of the envelope glycoprotein (GP). We also analyzed two additional CTL- variants and four spontaneous CTL+ revertants. All three CTL- variants differ from the original CTL+ parental strain at GP amino acid 260, indicating that this amino acid change is consistently associated with the CTL- phenotype. By contrast the other four mutations in LCMV are not associated with the CTL- phenotype. Sequence analysis of the coding regions of four CTL+ revertants of ARM clone 13 did not reveal back mutations at the GP 260 locus. This finding indicates that the GP 260 mutation is necessary but not sufficient for a CTL- P+ phenotype and that the reversion to CTL+ P- is likely either due to secondary mutations in other regions of the viral genome or to quasispecies within the revertant population that make significant contributions to the phenotype.     

Genetic variation occurring on the genome of an in vitro insertion mutant of poliovirus type 1.

An insertion sequence of 72 nucleotides prepared from a polylinker sequence of plasmid pUC18 was introduced at nucleotide position 702 of the 5' noncoding sequence (742 nucleotides long) of the genome of the Sabin strain of poliovirus type 1 by using an infectious cDNA clone of the virus strain. The insertion mutant thus obtained showed a small-plaque phenotype compared with that of the parent virus. Apparent revertants (large-plaque variants) were easily generated from the insertion mutant. Nucleotide sequence analysis was performed on the revertant genomes to determine the mutation(s) by which the plaque size of the parent virus was regained. Some large-plaque variants lacked genomic sequences including all or a part of the insertion sequence. A computer-aided search for secondary structures with respect to the deletion sites detected possible supporting sequences which provided fairly stable secondary structures at the deletion sites. This result was consistent with our supporting sequence-loop model which had been proposed as a new copy-choice model for the generation of genetic rearrangements occurring on single-stranded RNA genomes (S. Kuge, I. Saito, and A. Nomoto, J. Mol. Biol. 192:473-487, 1986). The other large-plaque variants had point mutations at any one of three positions of an AUG existing in the insertion sequence. A small-plaque phenotype was observed when an AUG codon was inserted in frame or out of frame with regard to the initiation site of viral polyprotein synthesis. Our data strongly suggest that an AUG sequence in this genome region is deleterious for efficient poliovirus replication.