Detection of myxoma viruses encoding a defective M135R gene from clinical cases of myxomatosis; possible implications for the role of the M135R protein as a virulence factor
© Belsham et al; licensee BioMed Central Ltd. 2010
Received: 23 November 2009
Accepted: 16 January 2010
Published: 16 January 2010
Myxoma virus is a member of the Poxviridae and causes disease in European rabbits. Laboratory confirmation of the clinical disease, which occurs in the autumn of most years in Denmark, has been achieved previously using antigen ELISA and electron microscopy.
An unusually large number of clinically suspected cases of myxomatosis were observed in Denmark during 2007. Myxoma virus DNA was detected, using a new real time PCR assay which targets the M029L gene, in over 70% of the clinical samples submitted for laboratory confirmation. Unexpectedly, further analysis revealed that a high proportion of these viral DNA preparations contained a frame-shift mutation within the M135R gene that has previously been identified as a virulence factor. This frame-shift mutation results in expression of a greatly truncated product. The same frame-shift mutation has also been found recently within an avirulent strain of myxoma virus (6918). However, three other frame-shift mutations found in this strain (in the genes M009L, M036L and M148R) were not shared with the Danish viruses but a single nucleotide deletion in the M138R/M139R intergenic region was a common feature.
It appears that expression of the full-length myxoma virus M135R protein is not required for virulence in rabbits. Hence, the frame-shift mutation in the M135R gene in the nonpathogenic 6918 virus strain is not sufficient to explain the attenuation of this myxoma virus but one/some of the other frame-shift mutations alone or in conjunction with one/some of the thirty two amino acid substitutions must also contribute. The real time PCR assay for myxoma virus is a useful diagnostic tool for laboratory confirmation of suspected cases of myxomatosis.
Myxomatosis is a disease primarily of European rabbits and is caused by infection with myxoma virus which can result in very high levels of mortality (nearly 100%). In consequence, the virus has been deliberately introduced into some countries as a control measure for wild rabbit populations. Although effective initially, co-evolution of virus and the rabbit population has resulted in the generation of attenuated field strains of the virus and also virus-resistant populations of wild rabbits [1–3].
Analysis of suspect myxomatosis samples submitted for laboratory analysis in Denmark in the period 1998-2007
Number of submissions
Number of positivea submissions (% positive)
Myxoma virus is a member of the Leporipoxvirus genus within the Poxviridae (see ). The virus has a dsDNA genome of about 160,000 bp and the complete sequence of the virulent Lausanne strain has been determined . The presence of numerous genes encoding immunomodulatory proteins and the narrow host range of the virus has made it an attractive system for analysing the role of these immunomodulators in pathogenesis [3, 5]. Poxviruses encode a variety of proteins that block host interferon and tumor necrosis factor mediated antiviral responses. Some of these proteins act within the cell, these include the vaccinia virus E3L and K3L proteins (plus their homologues from other viruses) that inhibit the double-stranded RNA activated protein kinase (PKR)-induced phosphorylation of the eukaryotic translation initiation factor 2 (eIF2) (reviewed in ). Other viral proteins act at the cell surface, for example, myxoma virus encodes a protein, called M-T7, which inhibits the interaction of interferon-γ with its cellular receptor . Recently, the product of the myxoma virus M135R gene has been identified as a novel cell surface virulence factor . The M135R protein, which is glycosylated, has some similarity (23% identity) to the vaccinia virus B18R protein which acts as an interferon-α/β receptor. Studies have shown, however, that the M135R protein, which is much shorter than the vaccinia virus B18R product, does not function in this way . However, this same study showed that targeted deletion of the M135R gene (removing 86% of the coding sequence) severely attenuated the virus in European rabbits, with all inoculated animals recovering from the moderate symptoms.
The availability of the myxoma virus genome sequence (, [GenBank:NC_001132]) enabled the design of a sensitive and specific real-time PCR-based assay for this virus. This assay has been used to identify the presence of myxoma virus DNA from within laboratory and clinical samples. Furthermore, sequence analysis of selected regions of these myxoma viruses has been performed to explore the relationship between current Danish viruses and other virulent and attenuated virus strains. Surprisingly, a high proportion of the Danish samples, obtained from clinically diseased animals, contained a frame-shift mutation within the M135R gene, which severely truncates the expressed protein product.
Viral DNA isolation
Cell culture grown samples of myxoma virus, vaccinia virus (WR strain), swinepox and Orf virus or rabbit tissue homogenates were used to prepare samples of viral DNA for analysis. The tissues (including nose, lips, eyelids and kidneys) from dead or euthanized pet rabbits which had exhibited clinical disease and had been provided for laboratory diagnosis, were homogenized in physiological saline (1 ml/g of tissue) and then clarified by centrifugation. The DNA was isolated using the QIAamp DNA Blood Mini kit (Qiagen) using the "Blood and Body Fluid Spin protocol" as described by the manufacturer. Briefly, supernatant samples (200 μl) were incubated with Qiagen protease at 56°C for 10 mins, mixed with an equal volume of ethanol and applied to the QIAamp Spin columns. Bound DNA was eluted in water (50 μl) and then assayed directly. It was found that using 1 μl of the eluted DNA was optimal in the real time PCR assay described below, for some samples higher amounts could produce a lower signal, presumably due to the presence of some inhibitor.
In addition, viral DNA was isolated, as described above, from a sample of Cunivak vaccine (supplied from Impfstoffwerk Dessau-Tornau GmbH, Germany) which contains a live attenuated strain of myxoma virus.
Quantitative real-time PCR detection assay
Oligonucleotide primers and probes used in this work
CCGTCGAC TTA CTCGAGCGCCCCGCCCTGCCACTC
Sequence comparison of the M135R gene from myxoma virus isolates
A standard PCR was performed using the isolated viral DNA preparations from a variety of different myxoma virus samples with the primers M135Rfor and M135Rrev (Table 2) which had been designed to amplify the entire coding sequence of the M135R gene (nt 131699-132235) plus some flanking sequences. The amplicons obtained (653 bp) were purified and sequenced, in both directions, using the same primers on an ABI 3130 automatic sequencer with a BigDye Terminator v3.1 Cycle sequencing kit (Applied Biosystems). The sequences were initially analysed using BioNumerics (Applied-Maths) software and the fasta files were then aligned with the myxoma virus sequence [GenBank:NC_001132] using ClustalW2 .
Expression analysis of M135R proteins
Analysis of M135R protein expression was achieved by making fusion constructs between the coding sequence for chloramphenicol acetyl transferase (CAT) and the M135R coding sequence. The CAT sequence was amplified in a standard PCR using primers CATBamfor and CATXSrev (Table 2) using a derivative of pSV2CAT  as template. These primers introduced BamHI, XhoI and SalI restriction sites plus an in-frame termination codon following the XhoI site. The product (739 bp) was inserted into pCR-XL-TOPO (Invitrogen). From the resulting plasmid, BamHI-XhoI and BamHI-SalI fragments were prepared and ligated into a BamHI-XhoI digested derivative of pGEM3Z (Promega) containing a modified polylinker, (T7-BamHI-SmaI-XhoI-EcoRI-ApaI-SalI, as described , to produce pGEMCATBS and pGEMCATBX respectively. The former plasmid includes a termination codon at the end of the CAT sequence while the latter does not. The myxoma M135R sequences were amplified from various different myxoma virus DNA samples using the primers M135RXFOR and M135RAPAREV (Table 2) which produced a product of 559 nt including the complete coding sequence of the M135R gene flanked by XhoI and ApaI sites. The products were individually inserted into pCR-XL-TOPO and then XhoI-ApaI fragments were isolated and ligated into similarly digested pGEMCATBX (as described above) to produce a fused open reading frame including the CAT and M135R coding sequences. Structures of these plasmids (see below) were confirmed by enzyme digestion and sequence analysis.
These plasmids were assayed by transfection into BHK cells that had been infected with the recombinant vaccinia virus vTF7-3  essentially as described previously  except that FuGene 6 (Roche) was used rather than lipofectin (see ). The vaccinia virus vTF7-3 expresses the T7 RNA polymerase and generates capped mRNAs from the transfected plasmids. Cell extracts were prepared 20 h post transfection and analysed by SDS-PAGE and immunoblotting. Expression of CAT and the CAT-M135R fusion proteins was detected using rabbit anti-CAT antibodies (Sigma) and peroxidase labelled anti-rabbit IgG (Dako) with chemiluminescent reagents (ECL Plus, GE Healthcare).
PCR amplification and sequence analysis of other myxoma virus regions
In addition to the M135R gene, five other regions of the myxoma virus genome were amplified by PCR using primers flanking the insertions and deletions identified within a nonpathogenic myxoma virus strain (6918) by Morales et al. . The primers, listed in Table 2, were designed to amplify a fragment of about 400-500 bp in each case to permit analysis of the M138L/M139R and M153R/M154R intergenic regions plus the coding regions from the M009L, M036L and M148R genes. Amplicons were gel purified and sequenced (by Agowa, Germany) in both directions using the primers used for the PCRs.
Design and use of a real time quantitative PCR assay for detection of myxoma virus
A real time quantitative PCR (qPCR) assay to detect specifically myxoma virus DNA was designed. The targeted region of the myxoma virus (a Leporipoxvirus) genome together with the homologous sequences from vaccinia virus (Orthopoxvirus), swinepox (Suipoxvirus) and Orf virus (Parapoxvirus) are shown in Figure 1 (panel a). The probe is identical to a region of the myxoma virus M029L gene and also to part of the E3L homolog within swinepox virus (SPV032), however the primer sequences selected to amplify the myxoma virus DNA are very different from the corresponding sequences in swinepox virus. The myxoma probe sequence is significantly different from the corresponding region of the vaccinia virus E3L sequence (7 nt differences out of 22) and from the Orf virus E3L homolog (only 13 nt of the Orf virus sequence match the myxoma virus probe).
To characterize the system for the detection of myxoma virus DNA, viral nucleic acid was prepared from cell culture grown samples of myxoma virus, vaccinia virus, swinepox and Orf virus and tested in the qPCR assay. As expected the assay efficiently detected the myxoma virus DNA but did not produce a signal using nucleic acid preparations from cells infected with the other poxviruses tested (see Figure 1, panel b). As controls, analogous qPCR assays for orthopoxviruses (HA gene, as described previously ) and parapoxviruses (B2L gene, encoding the major envelope protein, as described ) respectively were also performed. These assays correctly detected their target viruses within these samples, as expected, but did not detect the myxoma virus DNA (data not shown). The qPCR assay for myxoma virus amplified a fragment of 119 bp (including the primer sequences) and this product could be detected by agarose gel electrophoresis at the end of the assay (data not shown) as well as by the change in fluorescence measured in the assay (Figure 1, panel b). To determine the sensitivity of the myxoma virus assay, a region of the myxoma virus genome (409 bp, containing the entire M029L gene) was amplified by PCR from a viral DNA preparation and inserted into a plasmid vector. The resultant purified plasmid was quantified and a dilution series was then prepared and assayed in the qPCR to produce a standard curve (see Fig. 1, panel c). This showed that the assay could detect as little as 10 copies of the myxoma virus sequence in each sample and produced a linear response up to 107 copies per assay. Thus the assay has a wide range of sensitivity. The specificity of the assay for myxoma virus DNA was confirmed using conventional PCR assays, which targeted other regions of the myxoma virus genome, on DNA preparations from clinical and laboratory samples (see below).
The qPCR assay for myxoma virus has also been applied to 195 tissue samples from dead or euthanized rabbits which had displayed clinical signs of myxomatosis and were submitted by veterinarians to DTU Vet at Lindholm during 2007. From these samples, 142 (73%) proved to be positive (see Table 1) indicating that clinical suspicion had been accurate in most cases. Indeed, most of the clinical samples assayed produced a very high signal (Ct values less than 20). The great increase in the number of samples submitted for laboratory diagnosis is unrelated to the method used for the diagnosis but reflects an exceptional incidence of disease during 2007.
Comparison of myxoma virus isolates by sequence analysis of the M135R gene
Analysis of M135R protein expression from amplified M135R gene sequences
It was anticipated that the construction of plasmids containing the fused CAT-M135R reading frames should allow the expression of the normal M135R coding sequence to be observed and, more importantly, also the prematurely terminated M135R gene product which by itself would only be 40 amino acids in length, ca. 4.5 kDa, and therefore difficult to detect alone. Different plasmids containing either the wt M135R sequence (e.g. as in the myxoma vaccine) or the frame-shift variant (+G) were constructed and then assayed along with a plasmid encoding the normal CAT coding sequence (see Figure 3). The CAT vectors, pGEMCATBS and pGEMCATBX (both lacking any M135R sequences) very efficiently expressed the CAT protein (about 27 kDa) as expected (Figure 3, panel B). Fusing the wt M135R sequence to the CAT coding sequence (as in plasmid pCATM135Rwt) resulted in the expression of a much larger protein being recognized by the anti-CAT antibodies consistent with the expected fusion product of ca. 45 kDa (Figure 3, panel B, lanes 3 and 4). In contrast, the attachment of the mutant M135R sequences (as in pCATM135R+G), containing the insertion of a single G nt, resulted in a small, albeit readily detectable, decrease in the mobility of the CAT protein (Figure 3, panel B, lanes 5 and 6) consistent with the predicted addition of only an extra 40 amino acid residues (ca. 4.5 kDa) to the C-terminus of CAT. In addition, a smaller, weak band, ca. 22 kDa, was also detected from this plasmid, it seems likely that this is a degradation product derived from the fusion protein which may adopt a partially unfolded structure.
PCR and sequence analysis of Danish samples of myxoma virus DNA
Comparison of insertions/deletions within the nonpathogenic Spanish 6918 myxoma virus strain and Danish clinical samples
Country of origin/sample
Strain 6918 (mut)
TGTGGGGG G TG (+G)
A rapid, sensitive and specific quantitative PCR assay has been developed for myxoma virus which complements previously described qPCR assays for orthopoxviruses (e.g. vaccinia virus, ) and parapoxviruses (e.g. Orf virus, ). This assay has been shown to function on tissue culture grown viruses and also directly on material extracted from clinical samples. It removes the requirement for specific antisera (which can be variable) or an electron microscope.
Using this assay for myxoma virus, we tested samples from Danish rabbits with clinically determined myxomatosis and we confirmed the diagnosis of the unusually large number of cases of myxomatosis in Denmark during 2007 (see Table 1). The cases were predominantly located in the eastern part of Denmark around Copenhagen, however a small number of cases were also observed elsewhere (as in most previous years). Unexpectedly, a specific frame-shift mutation was present within the coding sequence of the M135R gene in most of the current Danish samples. The insertion of a single additional G, into a run of 5 G's, was present in 10 out of 13 samples analysed. This was the only sequence difference in the M135R gene between each of these samples and the myxoma virus vaccine (which is also identical in this gene to the pathogenic Lausanne strain, ). The same nucleotide insertion was also detected in subsequent, independent, PCR reactions which used a different primer set. The generation of these amplicons enabled the fusion of the various M135R gene sequences to the coding sequence of the CAT gene (lacking a termination codon). Expression of these fusion gene plasmids within mammalian cells, in transient expression assays, provided clear evidence that the frame-shift mutation, as detected by the DNA sequencing, indeed resulted in the synthesis of a greatly truncated protein from the mutant M135R sequence (Figure 3).
The M135R gene has recently been described as a virulence factor in myxoma virus  since targeted deletion (removing 86% of the coding region) of the gene resulted in severe attenuation of the virus within European rabbits. Experimental inoculation of rabbits with the modified virus gave rise to only moderate symptoms and all animals survived whereas all animals receiving viruses with the intact M135R gene developed severe disease and had to be euthanized. Hence it was surprising that viral DNA obtained from rabbits with clinical disease encoded a greatly truncated version of this gene product comprising only 19 (out of the total 178) amino acids (plus 21 residues translated in a different reading frame). Clearly, the presence of this frame-shift mutation does not severely attenuate the virus within rabbits, since the samples analysed came from animals which had died or were showing significant clinical disease which had lead veterinarians to euthanize the animals. This result is not without precedent. Initial studies, using gene deletion, aimed at determining whether particular poxvirus genes were essential for growth indicated that the vaccinia virus F11L gene was essential for growth , however this gene is split in the modified vaccinia Ankara strain indicating that the intact gene is not essential . Subsequent studies using an improved approach demonstrated that this gene is indeed non-essential and showed that misleading results can be obtained from the use of insertional mutagenesis to test for viability . It is, therefore, possible that gene deletion can give misleading data on the requirement of specific genes for virulence. Thus, we conclude from our observations that expression of the full length M135R gene product is not essential for virulence in European rabbits. It would be interesting to assess by transcription profiling (e.g., as described for vaccinia virus [20, 21]) whether the targeted deletion of the myxoma virus M135R gene has an influence on the expression of neighbouring genes.
Very recently, Morales et al.  reported the complete genome sequence of a nonpathogenic myxoma virus strain (6918) (see ) and found that it encoded thirty-two amino acid substitutions plus six insertions or deletions leading to four frame-shift mutations (plus two intergenic region changes) compared to the parental Lausanne strain of the virus (which was 99.95% identical in sequence). One of these frame-shift mutations, within strain 6918, is within the M135R gene and is identical to the change observed within most of the recent Danish samples analysed. It is interesting to note that analysis of recent Portuguese isolates of myxoma virus found evidence of another frame-shift mutation within a virulent strain of this virus which also occurred at a homopolymeric region .
To analyse the similarity between the Danish samples and the avirulent 6918 strain, five other regions of the myxoma virus genomes from Danish samples with or without the M135R insertion were amplified by PCR and sequenced. These analyses showed that none of the other frame-shift mutations were present in the Danish samples although the deletion of a single T from within the M138L/M139R intergenic region was a common feature (Table 3). Morales et al.  suggested that the frame-shift mutation within the M135R gene may be an important determinant of the attenuation of the 6918 strain since, as described above, deletion of the gene had lead to attenuation . However, the properties of the Danish samples would suggest that some other sequence changes within the 6918 strain are responsible for this phenotype. The critical change(s) could either be amongst the insertions/deletions that are not present in the Danish samples or within the changes leading to the thirty-two different amino acid substitutions.
1) A sensitive and specific qPCR assay has been established for myxoma virus.
2) An unusually high number of clinical cases of myxoma occurred in Denmark in 2007.
3) A high proportion of the Danish samples that were analysed had a frame-shift mutation within the M135R protein coding sequence which results in the expression of a severely truncated product.
4) The presence of this same frame-shift mutation within an attenuated myxoma virus (6918) is not sufficient to explain its attenuation.
5) Analysis of the effect of the frame-shift mutation within the M135R gene, in isolation, on the virulence of myxoma virus in European rabbits is warranted.
We thank Preben Normann, Inge Nielsen, Jiong Li, Kirsten Vestergård, Helle Rasmussen, Janne Holm Hansen and Randi B. Kølle for excellent assistance with these studies.
- Best SM, Kerr PJ: Coevolution of host and virus: the pathogenesis of virulent and attenuated strains of myxoma virus in resistant and susceptible European rabbits. Virology. 2000, 267: 36-48. 10.1006/viro.1999.0104.PubMedView Article
- Best SM, Collins SV, Kerr PJ: Coevolution of host and virus: cellular localization of virus in myxoma virus infection of resistant and susceptible European rabbits. Virology. 2000, 277: 76-91. 10.1006/viro.2000.0505.PubMedView Article
- Stanford MM, Werden SJ, McFadden G: Myxoma virus in the European rabbit: interactions between the virus and its susceptible host. Vet Res. 2007, 38: 299-318. 10.1051/vetres:2006054.PubMedView Article
- Cameron C, Hota-Michell S, Chen L, Barrett J, Cao J-X, Macauley C, Willer D, Evans D, McFadden G: The complete DNA sequence of myxoma virus. Virology. 1999, 264: 298-318. 10.1006/viro.1999.0001.PubMedView Article
- Nash P, Barrett J, Cao JX, Hota-Mitchell S, Lalani AS, Everett H, Xu XM, Robichaud J, Hnatiuk S, Ainslie C, Seet BT, McFadden G: Immunomodulation by viruses: the myxoma virus story. Immunol Rev. 1999, 168: 103-120. 10.1111/j.1600-065X.1999.tb01286.x.PubMedView Article
- García MA, Meurs EF, Esteban M: The dsRNA protein kinase PKR: virus and cell control. Biochimie. 2007, 89: 799-811. 10.1016/j.biochi.2007.03.001.PubMedView Article
- Upton C, Mossman K, McFadden G: Encoding of a homolog of the IFN-gamma receptor by myxoma virus. Science. 1992, 258: 1369-1372. 10.1126/science.1455233.PubMedView Article
- Barrett JW, Sypula J, Wang F, Alston LR, Shao Z, Gao X, Irvine TS, McFadden G: M135R is a novel cell surface virulence factor of myxoma virus. J Virol. 2007, 81: 106-114. 10.1128/JVI.01633-06.PubMedPubMed CentralView Article
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: ClustalW and ClustalX version 2. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.PubMedView Article
- Gorman CM, Moffat LF, Howard BH: Recombinant genomes which express chloramphenicol acetyltransferase in mammalian cells. Mol Cell Biol. 1982, 2: 1044-51.PubMedPubMed CentralView Article
- Belsham GJ, Nielsen I, Normann P, Royall E, Roberts LO: Monocistronic mRNAs containing defective hepatitis C virus-like picornavirus internal ribosome entry site elements in their 5'untranslated regions are efficiently translated in cells by a cap-dependent mechanism. RNA. 2008, 14: 1671-1680. 10.1261/rna.1039708.PubMedPubMed CentralView Article
- Fuerst TR, Niles EG, Studier FW, Moss B: Eukaryotic transient expression system based on recombinant vaccinia virus that synthesizes bacteriophage T7 RNA polymerase. Proc Natl Acad Sci USA. 1986, 83: 8122-8126. 10.1073/pnas.83.21.8122.PubMedPubMed CentralView Article
- Roberts LO, Seamons RA, Belsham GJ: Recognition of picornavirus internal ribosome entry sites within cells; influence of cellular and viral proteins. RNA. 1998, 4: 520-529. 10.1017/S1355838298971989.PubMedPubMed CentralView Article
- Morales M, Ramírez MA, Cano MJ, Párraga M, Castilla J, Pérez-Ordoyo LI, Torres JM, Bárcena J: Genome comparison of a nonpathogenic myxoma virus field strain with its ancestor, the virulent Lausanne strain. J Virol. 2009, 83: 2397-2403. 10.1128/JVI.02189-08.PubMedPubMed CentralView Article
- Kulesh DA, Baker RO, Loveless BM, Norwood D, Zwiers SH, Mucker E, Hartman C, Herrera R, Miller D, Christensen D, Wasieloski LP, Huggins J, Jahrling PB: Smallpox and pan-orthopox virus detection by real-time 3'-minor groove binder TaqMan assays on the Roche Lightcycler and the Cepheid smart cycler platforms. J Clin Micro. 2004, 42: 601-609. 10.1128/JCM.42.2.601-609.2004.View Article
- Nitsche A, Büttner M, Wilhelm S, Pauli G, Meyer H: Real-time PCR detection of parapoxvirus DNA. Clinical Chemistry. 2006, 52: 316-319. 10.1373/clinchem.2005.060335.PubMedView Article
- Bárcena J, Pagès-Manté A, March R, Morales M, Ramírez MA, Sánchez-Vizcaíno JM, Torres JM: Isolation of an attenuated myxoma virus field strain that can confer protection against myxomatosis on contacts of vaccinates. Arch Virol. 2000, 145: 759-771. 10.1007/s007050050669.PubMedView Article
- Kato SEM, Greco FAB, Damaso CRA, Condit RC, Moussatche N: An alternative genetic method to test essential vaccinia virus early genes. J Virol Met. 2004, 115: 31-40. 10.1016/j.jviromet.2003.09.011.View Article
- Antoine G, Scheiflinger F, Dorner F, Falkner FG: The complete genomic sequence of the modified vaccinia Ankara strain: comparison with other orthopoxviruses. Virology. 1998, 44: 365-396. 10.1006/viro.1998.9123.View Article
- Assarsson E, Greenbaum JA, Sundström M, Schaffer L, Hammond JA, Pasquetto V, Oseroff C, Hendrickson RC, Lefkowitz EJ, Tscharke DC, Sidney J, Grey HM, Head SR, Peters B, Sette A: Kinetic analysis of a complete poxvirus transcriptome reveals an immediate-early class of genes. Proc Natl Acad Sci USA. 2008, 105: 2140-2145. 10.1073/pnas.0711573105.PubMedPubMed CentralView Article
- Rubins KH, Hensley LE, Bell GW, Wang C, Lefkowitz EJ, Brown PO, Relman DA: Comparative analysis of viral gene expression programs during poxvirus infection: a transcriptional map of the vaccinia and monkeypox genomes. PLoS ONE. 2008, 3 (e2628):
- Muller A, Silva E, Abrantes J, Esteves PJ, Ferreira PG, Carvalheira JC, Nowotny N, Thompson G: Partial sequencing of recent Portuguese myxoma virus field isolates exhibits a high degree of genetic stability. Vet Microbiol. 2010, 140: 161-166. 10.1016/j.vetmic.2009.07.028.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.