Assignment of the large oligonucleotides of vesicular stomatitis virus to the N, NS, M, G, and L genes and oligonucleotide gene ordering within the L gene.

Analyses of prototype vesicular stomatitis (VSV, Indiana serotype) mRNA-32P-labeled viral RNA duplexes have established the assignments of 65 of the 72 large oligonucleotides that are recovered by two-dimensional electrophoresis of RNase T1 digests of the viral RNA. Fifty of the oligonucleotides are recovered in the L RNA duplex, four each in the N, M, and NS duplexes, and three in the G RNA duplex. Studies of three small defective-particle RNA species indicate that they have only L gene oligonucleotides in addition to three of the seven unassigned oligonucleotides. Some L gene ordering of oligonucleotides can be postulated from the defective-particle RNA sequence analyses. Analyses of naturally occurring alternate isolates of VSV Indiana have established that by comparison to the prototype virus strain, the alternate isolates minimally have genome sequence differences in L, G, N, NS and/or unassigned regions of the genome. Changes in the genome have also been induced by vitro high-level mutagenesis of the prototype virus.     

Comparative genomics of foot-and-mouth disease virus

Here we present complete genome sequences, including a comparative analysis, of 103 isolates of foot-and-mouth disease virus (FMDV) representing all seven serotypes and including the first complete sequences of the SAT1 and SAT3 genomes. The data reveal novel highly conserved genomic regions, indicating functional constraints for variability as well as novel viral genomic motifs with likely biological relevance. Previously undescribed invariant motifs were identified in the 5' and 3' untranslated regions (UTR), as was tolerance for insertions/deletions in the 5' UTR. Fifty-eight percent of the amino acids encoded by FMDV isolates are invariant, suggesting that these residues are critical for virus biology. Novel, conserved sequence motifs with likely functional significance were identified within proteins L(pro), 1B, 1D, and 3C. An analysis of the complete FMDV genomes indicated phylogenetic incongruities between different genomic regions which were suggestive of interserotypic recombination. Additionally, a novel SAT virus lineage containing nonstructural protein-encoding regions distinct from other SAT and Euroasiatic lineages was identified. Insights into viral RNA sequence conservation and variability and genetic diversity in nature will likely impact our understanding of FMDV infections, host range, and transmission.     

Complexity and polyadenylic acid content of visna virus 60-70S RNA.

The genomic complexity of visna virus was measured by quantitative analysis of 18 RNase T1-resistant oligonucleotides from 60-70S RNA. T1-resistant oligonucleotides were separated by two-dimensional polyacrylamide gel electrophoresis. Visna virus had a genomic complexity of 3.6 X 10(6) daltons, very close to the size of a single 30-40S RNA subunit. It was therefore concluded that the visna virus genome is largely polyploid. Visna virus 60-70S RNA polyadenylic acid segment was purified by T1 RNase digestion followed by oligodeoxythymidylic acid-cellulose column chromatography. It contained over 99% AMP and had a size of about 200 nucleotides. The binding capacities on oligodeoxythymidylic acid-cellulose of native 60-70S RNA and purified 30-40S RNA subunits were examined. It was concluded that two out of three intact subunits contain a polyadenylic acid segment.     

RNA synthesis of vesicular stomatitis virus. VIII. Oligonucleotides of the structural genes and mRNA.

The single-stranded RNA genome of vesicular stomatitis virus (VSV, Indiana serotype, San Juan strain) yields approx. 75 RNase T1-resistant oligonucleotides ranging in size from 10 to 50 bases. Each of the five structural genes, isolated as duplex RNA molecules hybridized to complementary mRNA, contains two or more of these large oligonucleotides. One of the oligonucleotides is identified as part of the non-coding region near the 3' end of the genome. Comparison of these results with others indicate that the RNA sequence of VSV is apparently stable in the laboratory but not in the wild. RNase T1-resistant oligonucleotides are also shown for all five VSV mRN species. Whether the mRNA for these digestions are are isolated from duplex RNA molecules or as single-stranded RNA species, the oligonucleotide patterns for each mRNA are virtually identical, indicating that each mRNA is transcribed from contiguous sequences on the genome. Comparison with published oligonucleotide patterns obtained from other isolates of VSV or from VSV deletion mutants indicate that identity and changes in their genome structure can be correlated with specific structural genes.     

Genomic complexities of murine leukemia and sarcoma, reticuloendotheliosis, and visna viruses.

The genetic complexities of several ribodeoxyviruses were measured by quantitative analysis of unique RNase T1-resistant oligonucleotides from 60-70S viral RNAs. Moloney murine leukemia virus was found to have an RNA complexity of 3.5 x 10(6) daltons, whereas Moloney murine sarcoma virus had a significantly smaller genome size of 2.3 x 10(6). Reticuleondotheliosis and visna virus RNAs had complexities of 3.9 x 10(6), respectively. Analysis of RNase A-resistant oligonucleotides of Rous sarcoma virus RNA gave a complexity of 3.6 x 10(6), similar to that previously obtained with RNase T1-resistant oligonucleotides. Since each of these viruses was found to have a unique sequence genomic complexity near the molecular weight of a single 30-40S viral RNA subunit, it was concluded that ribodeoxyvirus genomes are at least largely polyploid.     

Oligonucleotide fingerprints of RNA species obtained from rhabdoviruses belonging to the vesicular stomatitis virus subgroup.

The relationships among the genomes of various rhabdoviruses belonging to the vesicular stomatitis virus subgroup were analyzed by an oligonucleotide fingerprinting technique. Of 10 vesicular stomatitis viruses, Indiana serotype (VSV Indiana), obtained from various sources, either no, few, or many differences were observed in the oligonucleotide fingerprints of the 42S RNA species extracted from standard B virions. Analyses of the oligonucleotides obtained from RNA extracted from three separate preparations of VSV Indiana defective T particles showed that their RNAs contain fewer oligonucleotides than the corresponding B particle RNA species. The fingerprints of RNA obtained from five VSV New Jersey serotype viruses were easily distinguished from those of the VSV Indiana isolates. Three of the VSV New Jersey RNA fingerprints were similar to each other but quite different from those of the other two viruses. The RNA fingerprints of two Chandipura virus isolates (one obtained from India and one from Nigeria) were also unique, whereas the fingerprint of Cocal virus RNA was unlike that of the serologically related VSV Indiana.     

Fixation of mutations in the viral genome during an outbreak of foot-and-mouth disease: heterogeneity and rate variations

Rates of fixation of mutations during the evolution of the foot-and-mouth disease virus (FMDV) C₁ in nature have been estimated by hybridization of viral RNA to cloned cDNAs representing defined FMDV genome segments, and comparison of the selected RNAs by T₁ RNase oligonucleotide fingerprinting. Values ranged from <0.04 × 10⁻² to 4.5 × 10⁻² substitutions per nucleotide per year (s/nt/yr), depending on the time period and the genomic segment considered. Rates for viral structural protein genes were up to sixfold higher than for nonstructural protein genes. Values in excess of 10⁻² s/nt/yr have been measured for the RNA region that encodes VP1–VP3. The nucleotide sequences of the major immunogenic region of capsid protein VP1 have been determined for six new FMDV C₁ isolates, and they are compared with the two previously known sequences of FMDV C₁ (C-S8 and C₁-O). Both oligonucleotide fingerprinting of selected RNA fragments and direct nucleotide sequencing demonstrate that genetic heterogeneity exists among three viruses isolated on the same day, introducing a significant indetermination in the evaluation of fixation rates of mutations. During the FMDV C₁ outbreak, amino acid substitutions did occur that are known to affect the immunological properties of the virus. The proportion of mutations between two viral RNAs does not increase significantly with the time elapsed between the two isolations, suggesting a cocirculation of multiple, related, nonidentical FMDVs (‘evolving quasispecies’) as the mode of evolution of this agent.     

Physical map of the Kirsten sarcoma virus genome as determined by fingerprinting RNase T1-resistant oligonucleotides.

From analysis of the large RNase T1-resistant oligonucleotides of Kirsten sarcoma virus (Ki-SV), a physical map of the virus genome was deduced. Kirsten murine leukemia virus (Ki-MuLV) sequences were detected in T1 oligonucleotides closest to the 3' end of the viral RNA and extended approximately 1,000 nucleotides into the genome. The rat genetic sequences started at this point and extended all the way to the very 5' end of the RNA molecules, where a small stretch of Ki-MuLV sequence was detected. By comparison of the fingerprints of Ki-SV RNA and the RNA of the endogenous rat src genetic sequences, it was found that more than 50% of the T1 oligonucleotides were similar between Ki-SV and the endogenous rat src RNA, suggesting an identical primary nucleotide sequence in over 50% of the viral genomes. The results indicate that Ki-SV arose by recombination between the 5' and 3' ends of Ki-MuLV and a large portion of the homologous sequences of the endogenous rat src RNA.     

Rapid emergence of novel antigenic and genetic variants of equine infectious anemia virus during persistent infection.

Previous results from our laboratory have demonstrated that equine infectious anemia virus displays structural variations in its surface glycoproteins and RNA genome during passage and chronic infections in experimentally infected Shetland ponies (Montelaro et al., J. Biol. Chem. 259:10539-10544, 1984; Payne et al., J. Gen. Virol. 65:1395-1399, 1984). The present study was undertaken to obtain an antigenic and biochemical characterization of equine infectious anemia virus isolates recovered from an experimentally infected pony during sequential disease episodes, each separated by intervals of only 4 to 8 weeks. The virus isolates could be distinguished antigenically by neutralization assays with serum from the infected pony and by Western blot analysis with a monoclonal antibody against the major surface glycoprotein gp90, thus demonstrating that novel antigenic variants of equine infectious anemia virus predominate during each clinical episode. The respective virion glycoproteins displayed different electrophoretic mobilities on sodium dodecyl sulfate-polyacrylamide gels, indicating structural variation. Tryptic peptide and glycopeptide maps of the viral proteins of each virus isolate revealed biochemical alterations involving amino acid sequence and glycosylation patterns in the virion surface glycoproteins gp90 and gp45. In contrast, no structural variation was observed in the internal viral proteins pp15, p26, and p9 from any of the four virus isolates. Oligonucleotide mapping experiments revealed similar but unique RNase T1-resistant oligonucleotide fingerprints of the RNA genomes of each of the virus isolates. Localization of altered oligonucleotides for one virus isolate placed two of three unique oligonucleotides within the predicted env gene region of the genome, perhaps correlating with the structural variation observed in the envelope glycoproteins. Thus these results support the concept that equine infectious anemia virus is indeed capable of relatively rapid genomic variations during replication, some of which result in altered glycoprotein structures and antigenic variants which are responsible for the unique periodic disease nature observed in persistently infected animals. The findings of envelope specific differences in isolates of visna virus and of human T-cell lymphotropic virus III (acquired immune deficiency syndrome-related virus) suggest that this variation may be a common characteristic of the subfamily Lentivirinae.     

Sequence and location of the poly C tract in aphtho- and cardiovirus RNA.

The poly C tract in the RNA of the aphtho- and cardio viruses has been examined in several isolates of foot-and-mouth disease virus (FMDV) and encephalomyocarditis (EMC) virus. The length of the tract is variable, containing 100 to 170 bases in the FMDV isolates and 80 to 250 bases in the EMC virus isolates. Each poly C tract contains c. 10% A and U residues, located at the 5' end, i.e. most of the tract is a continuous run of C residues. The position of the tract on the genome was the same in each of the FMDV isolates, about 400 bases from the 5' end, whereas in the EMC virus isolates it was about 150 bases from the 5' end.     

Characterization of virus produced by a lymphoma induced by inoculation of AKR MCF-247 virus.

We report the characterization of the virus produced by a lymphoid cell line derived from a lymphoma of an AKR mouse after injection of the polytropic AKR virus MCF-247. The virus displays polytropic host range properties and is indistinguishable from MCF-247 as judged by analysis of the large RNase T1-resistant oligonucleotides of the RNA genome. Restriction enzyme analysis of cellular DNA revealed the presence of sequences homologous to MCF-247 genomic RNA. The EcoRI cleavage fragments were characteristic of MCF-247 DNA provirus cleavage products.     

Genetic structure, transforming sequence, and gene product of avian sarcoma virus UR1

We analyzed the genetic structure and gene products of the newly isolated avian sarcoma virus UR1, which recently has been shown to be replication defective and to contain no sequences homologous to the src gene of Rous sarcoma virus. The sizes of the genomic RNAs of UR1 and its associated helper virus, UR1AV, were determined to be 29S and 35S (5.9 and 8.5 kilobases), respectively, by gel electrophoresis and sucrose gradient sedimentation. RNase T1 oligonucleotide mapping of purified viral RNAs indicated that UR1 RNA contains eight unique oligonucleotides in the middle of the genome and shares four 5'-terminal and three 3'-terminal oligonucleotides with UR1AV RNA. The unique sequences of UR1 and Fujinami sarcoma virus were found to be closely related to each other by molecular hybridization of UR1 RNA with DNA complementary to the unique sequence of Fujinami sarcoma virus RNA, but minor differences were found by oligonucleotides fingerprinting. In the regions flanking the unique sequences, UR1 and Fujinami sarcoma viral RNAs contain distinct oligonucleotides, which are shared with oligonucleotides of the respective helper viral RNAs. Cell transformed with UR1 produce a single 29S RNA species which contains a UR1 unique sequence; this species is most likely the mRNA coding for the transforming protein. In UR1-transformed cells, a phosphoprotein fo 150,000 daltons (p150) was detected by immunoprecipitation with antiserum against gag proteins. p150 was associated with a protein kinase activity that was capable of phosphorylating p150 itself, immunoglobulin G of antiserum, and a soluble substrate, alpha-casein. This enzyme transferred phosphate exclusively to tyrosine residues of substrates in vitro, but p 150 labeled in vivo with 32P contained both phosphoserine and phosphotyrosine. The in vitro kinase reaction was not affected by the presence of cyclic AMP or cyclic GMP and strongly preferred Mn2+ over Mg2+. Thus, the properties of UR1 protein are almost identical to those of Fujinami sarcoma virus protein.     

Genome Sequences of SAT 2 Foot-and-Mouth Disease Viruses from Egypt and Palestinian Autonomous Territories (Gaza Strip)

Two foot-and-mouth disease virus (FMDV) genome sequences have been determined for isolates collected from recent field outbreaks in North Africa (Egypt) and the Middle East (Palestinian Autonomous Territories). These data represent the first examples of complete genomic sequences for the FMDV SAT 2 topotype VII, which is thought to be endemic in countries immediately to the south of the Sahara desert. Further studies are now urgently required to provide insights into the epidemiological links between these outbreaks and to define the pathogenicity of this emerging lineage.     

Genetic structure and transforming sequence of avian sarcoma virus UR2.

We have recently shown that a newly isolated avian sarcoma virus, UR2, is defective in replication and contains no sequences homologous to the src gene of Rous sarcoma virus. In this study, we analyzed the genetic structure and transforming sequence of UR2 by oligonucleotide fingerprinting. The sizes of the genomic RNAs of UR2 and its associated helper virus, UR2AV, were determined to be 24S and 35S, respectively, by sucrose gradient sedimentation. The molecular weight of the 24S UR2 genomic RNA was estimated to be 1.1 x 10(6), corresponding to 3,300 nucleotides, by gel electrophoresis under the native and denatured conditions. RNase T1 oligonucleotide mapping indicated that UR2 RNA contains seven unique oligonucleotides in the middle of the genome and shares eight 5'- and six 3'-terminal oligonucleotides with UR2AV RNA. From these data, we estimated that UR2 RNA contains a unique sequence of about 12 kilobases in the middle of the genome, and contains 1.4 and 0.7 kilobases of sequences shared with UR2AV RNA at the 5' and 3' ends, respectively. Partial sequence analysis of the UR2-specific oligonucleotides by RNase A digestion revealed that there are no homologous counterparts to these oligonucleotides in the RNAs of other avian sarcoma and acute leukemia viruses studied to date. UR2-transformed non-virus-producing cells contain a single 24S viral RNA which is most likely the message coding for the transforming protein of UR2. On the basis of the uniqueness of the transforming sequence, we concluded that UR2 is a new member of the defective avian sarcoma viruses.     

Absence of circularly permuted and largely redundant sequences in the genome of visna virus.

A previous study of the infectivity of visna virus proviral DNA suggested that the genetic information of the virus is distributed over at least two of the RNA subunits. Because the genetic complexity of visna virus corresponds to the size of one subunit, this result may imply that sequence redundancies exist within each subunit. In the present article we have examined this question by constructing a map of the large RNase T1-resistant oligonucleotides of the viral genome. Our principal results are as follows: (i) all 36S RNA subunits have the same genetic content regardless of their polyadenylic acid [poly(A)] content; (ii) the poly(A) tract is present at the 3' end of the molecule; and (iii) the recoveries of 19 large RNase T1-resistant oligonucleotides from poly(A)-tagged RNA fragments of various sizes demonstrate that the oligonucleotides are organized in the same linear order within all subunits. Our results, therefore, exclude the existence of large sequence redundancies in the genome of visna virus.     

Biochemical analysis of murine leukemia viruses isolated from radiation-induced leukemias of strain BALB/c.

Murine leukemia viruses isolated from radiation-induced BALB/c leukemias were characterized with respect to viral proteins and RNA. Analysis by sodium dodecyl sulfate-polyacrylamide gel electrophoresis of the viral structural proteins revealed that for p12, p15, p30, and gp70, three of four electrophoretic variants of each could be detected. There was no correlation found between any of these mobilities and N- or B-tropism of the viruses. Proteins of all xenotropic viral isolates were identical in their gel electrophoretic profiles. The similar phenotypes of multiple viral clones from individual leukemias and of isolates grown in different cells suggest that the polymorphism of ecotropic viruses was generated in vivo rather than during in vitro virus growth. By two-dimensional fingerprinting of RNase T1-resistant oligonucleotides from 70S viral DNA, the previously reported association of N- and B-tropism with two distinct oligonucleotides was confirmed. The presence of two other oligonucleotides was correlated with positive and negative phenotypes of the virus-coded GIX cell surface antigen. The RNAs of two B-tropic isolates with distinctive p15 and p12 phenotypes differed from the RNA of a prototype N-tropic virus by the absence of three oligonucleotides mapping in the 5' portion (gag region) of the prototype RNA. In addition, one small-plaque B-tropic virus displayed extensive changes in the RNA sequences associated with the env region of the prototype.     

A universal protocol to generate consensus level genome sequences for foot-and-mouth disease virus and other positive-sense polyadenylated RNA viruses using the Illumina MiSeq

Background

Next-Generation Sequencing (NGS) is revolutionizing molecular epidemiology by providing new approaches to undertake whole genome sequencing (WGS) in diagnostic settings for a variety of human and veterinary pathogens. Previous sequencing protocols have been subject to biases such as those encountered during PCR amplification and cell culture, or are restricted by the need for large quantities of starting material. We describe here a simple and robust methodology for the generation of whole genome sequences on the Illumina MiSeq. This protocol is specific for foot-and-mouth disease virus (FMDV) or other polyadenylated RNA viruses and circumvents both the use of PCR and the requirement for large amounts of initial template.

Results

The protocol was successfully validated using five FMDV positive clinical samples from the 2001 epidemic in the United Kingdom, as well as a panel of representative viruses from all seven serotypes. In addition, this protocol was successfully used to recover 94% of an FMDV genome that had previously been identified as cell culture negative. Genome sequences from three other non-FMDV polyadenylated RNA viruses (EMCV, ERAV, VESV) were also obtained with minor protocol amendments. We calculated that a minimum coverage depth of 22 reads was required to produce an accurate consensus sequence for FMDV O. This was achieved in 5 FMDV/O/UKG isolates and the type O FMDV from the serotype panel with the exception of the 5′ genomic termini and area immediately flanking the poly(C) region.

Conclusions

We have developed a universal WGS method for FMDV and other polyadenylated RNA viruses. This method works successfully from a limited quantity of starting material and eliminates the requirement for genome-specific PCR amplification. This protocol has the potential to generate consensus-level sequences within a routine high-throughput diagnostic environment.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-828) contains supplementary material, which is available to authorized users.     

Foot-and-mouth disease virus localisation on follicular dendritic cells and sustained induction of neutralising antibodies is dependent on binding to complement receptors (CR2/CR1)

Previous studies have shown after the resolution of acute infection and viraemia, foot-and-mouth disease virus (FMDV) capsid proteins and/or genome are localised in the light zone of germinal centres of lymphoid tissue in cattle and African buffalo. The pattern of staining for FMDV proteins was consistent with the virus binding to follicular dendritic cells (FDCs). We have now demonstrated a similar pattern of FMDV protein staining in mouse spleens after acute infection and showed FMDV proteins are colocalised with FDCs. Blocking antigen binding to complement receptor type 2 and 1 (CR2/CR1) prior to infection with FMDV significantly reduced the detection of viral proteins on FDCs and FMDV genomic RNA in spleen samples. Blocking the receptors prior to infection also significantly reduced neutralising antibody titres, through significant reduction in their avidity to the FMDV capsid. Therefore, the binding of FMDV to FDCs and sustained induction of neutralising antibody responses are dependent on FMDV binding to CR2/CR1 in mice.     

Genetic and phenotypic variation of foot-and-mouth disease virus during serial passages in a natural host

Foot-and-mouth disease virus (FMDV), like other RNA viruses, exhibits high mutation rates during replication that have been suggested to be of adaptive value. However, even though genetic variation in RNA viruses and, more specifically, FMDV has been extensively examined during virus replication in a wide variety of in vitro cell cultures, very little is known regarding the generation and effects of genetic variability of virus replication in the natural host under experimental conditions and no genetic data are available regarding the effects of serial passage in natural hosts. Here, we present the results of 20 serial contact transmissions of the highly pathogenic, pig-adapted O Taiwan 97 (O Tw97) isolate of FMDV in swine. We examined the virus genomic consensus sequences for a total of 37 full-length viral genomes recovered from 20 in vivo passages. The characteristics and distributions of changes in the sequences during the series of pig infections were analyzed in comparison to the O Tw97 genomes recovered from serially infected BHK-21 cell cultures. Unexpectedly, a significant reduction of virulence upon pig passages was observed, and finally, interruption of the viral transmission chain occurred after the14th pig passage (T14). Virus was, however, isolated from the tonsils and nasal swabs of the asymptomatic T15 pigs at 26 days postcontact, consistent with a natural establishment of the carrier state previously described only for ruminants. Surprisingly, the region encoding the capsid protein VP1 (1D) did not show amino acid changes during in vivo passages. These data demonstrate that contact transmission of FMDV O Tw97 in pigs mimics the fitness loss induced by the bottleneck effect, which was previously observed by others during plaque-to-plaque FMDV passage in vitro, suggesting that unknown mechanisms of virulence recovery might be necessary during the evolution and perpetuation of FMDV in nature.     

The genome of Uukuniemi virus consists of three unique RNA segments

The three RNA species isolated from virions of Uukuniemi virus, a proposed member of the newly defined Bunyaviridae family, have been characterized by analysis of ³²P-labeled ribonuclease T1 oligonucleotides separated on two-dimensional polyacrylamide gels. Each RNA species contains unique oligonucleotides not present in the two others, indicating that the genome of this virus is segmented. Each segment appears to contain a unique primary sequence with little or no overlapping among the segments. The complexities of the RNA segments as calculated from the radioactivity in unique oligonucleotides of defined lengths are about 8000 (L RNA), 3500 (M) and 1900 (S) nucleotides. Since these values are similar to the molecular weights determined by other methods, each size class of RNA corresponds to a single molecular species. The presence of a 5′ terminal pppAp … structure in each RNA segment confirms indications from electron microscopy that the apparently circular RNA segments are not covalently closed. The absence of either a 5′ terminal “cap” or 3′ terminal poly(A) supports the concept that Uukuniemi virus is a negative strand virus.