Disease-associated XMRV sequences are consistent with laboratory contamination
- Stéphane Hué†1,
- Eleanor R Gray†1,
- Astrid Gall†2,
- Aris Katzourakis3,
- Choon Ping Tan1,
- Charlotte J Houldcroft2,
- Stuart McLaren2,
- Deenan Pillay1,
- Andrew Futreal2,
- Jeremy A Garson1,
- Oliver G Pybus3,
- Paul Kellam1, 2Email author and
- Greg J Towers1Email author
© Hué et al; licensee BioMed Central Ltd. 2010
Received: 2 December 2010
Accepted: 20 December 2010
Published: 20 December 2010
Xenotropic murine leukaemia viruses (MLV-X) are endogenous gammaretroviruses that infect cells from many species, including humans. Xenotropic murine leukaemia virus-related virus (XMRV) is a retrovirus that has been the subject of intense debate since its detection in samples from humans with prostate cancer (PC) and chronic fatigue syndrome (CFS). Controversy has arisen from the failure of some studies to detect XMRV in PC or CFS patients and from inconsistent detection of XMRV in healthy controls.
Here we demonstrate that Taqman PCR primers previously described as XMRV-specific can amplify common murine endogenous viral sequences from mouse suggesting that mouse DNA can contaminate patient samples and confound specific XMRV detection. To consider the provenance of XMRV we sequenced XMRV from the cell line 22Rv1, which is infected with an MLV-X that is indistinguishable from patient derived XMRV. Bayesian phylogenies clearly show that XMRV sequences reportedly derived from unlinked patients form a monophyletic clade with interspersed 22Rv1 clones (posterior probability >0.99). The cell line-derived sequences are ancestral to the patient-derived sequences (posterior probability >0.99). Furthermore, pol sequences apparently amplified from PC patient material (VP29 and VP184) are recombinants of XMRV and Moloney MLV (MoMLV) a virus with an envelope that lacks tropism for human cells. Considering the diversity of XMRV we show that the mean pairwise genetic distance among env and pol 22Rv1-derived sequences exceeds that of patient-associated sequences (Wilcoxon rank sum test: p = 0.005 and p < 0.001 for pol and env, respectively). Thus XMRV sequences acquire diversity in a cell line but not in patient samples. These observations are difficult to reconcile with the hypothesis that published XMRV sequences are related by a process of infectious transmission.
We provide several independent lines of evidence that XMRV detected by sensitive PCR methods in patient samples is the likely result of PCR contamination with mouse DNA and that the described clones of XMRV arose from the tumour cell line 22Rv1, which was probably infected with XMRV during xenografting in mice. We propose that XMRV might not be a genuine human pathogen.
XMRV (Xenotropic murine leukaemia virus-related virus) is a xenotropic murine leukaemia virus (MLV-X) that has been detected in samples from prostate cancer (PC) and chronic fatigue syndrome (CFS) patients[1–6]. This has led to the suggestion that infection with this virus might cause these conditions. MLV-Xs are endogenous gamma retroviruses found in the genomes of mice. They are so named because in vitro they infect cells from a variety of species but were originally found not to infect the inbred strains of mice from which they were derived, due to mutations in the host xenotropic receptor. More recently, murine xenotropic receptor variants have been described which support MLV-X infection revealing a complex evolutionary relationship between MLV-X envelope sequences and their receptors in rodents [7–9]. XMRV has also been detected in 1-6% of healthy human controls in some studies, suggesting that infection may be common in the healthy human population [2, 3, 5]. The association between XMRV and human disease is controversial, with some studies detecting XMRV in up to 67% of patients whilst others have failed to detect XMRV infection [10–18]. Importantly, examination of infected prostate tumours reveals that not all the tumour cells are infected with XMRV suggesting that XMRV insertion is not required for tumourogenesis . XMRV sequences detected in patients are remarkably similar to each other often differing by only a few nucleotides between unlinked patients . This lack of sequence variation appears inconsistent with a retrovirus infecting geographically separated, unconnected individuals. Here, we have examined the specificity of XMRV PCR, XMRV sequence variation and the phylogenetic relationship between XMRV detected in humans and as contaminants in cell culture. We conclude that XMRV in patient samples is likely to be derived from PCR contamination from either mouse DNA or cell lines infected with XMRV, and that XMRV is unlikely to be a human pathogen.
Results and Discussion
Primers reported to be XMRV specific can detect mouse DNA
Screening of genomic DNA from 9 inbred and 3 wild-derived inbred mice using Taqman PCR
Taqman PCR result (Ct FAM)
454 deep sequencing
XMRV - gag-leader
XMRV - int
MLV-X - gag
Human cell lines are commonly contaminated with xenotropic MLVs
Screening human cell line genomic DNA using Taqman PCR
Cancer cell line
Taqman PCR Result (Ct FAM)
Squamous cell carcinoma
Lymphoid neoplasm/Burkitt lymphoma
Haematopoietic neoplasm/Myeloid leukaemia
Phylogenetic analysis of XMRV sequences
Analysis of the genetic diversity and phylogenetic relationships among retroviral sequences, both endogenous  and exogenous  can reveal information about their replication and evolutionary history. We therefore performed extensive evolutionary analysis of published XMRV and related sequences in order to better understand their origin and proliferation. The widely studied prostate cancer line 22Rv1 is reported to produce high levels of a virus closely related to XMRV [25, 26]. We therefore cloned and sequenced gag (n = 16), pol (n = 18) and env (n = 10) PCR products amplified from genomic DNA purified from this cell line. Of these, 13/16 gag sequences, 15/18 pol sequences and 8/10 env sequences were unique. These numbers are consistent with the previously estimated XMRV copy number in the 22Rv1 cell line of 10-20 copies . We analysed these unique cell line sequences together with (i) previously-described full-length endogenous MLV genomes (n = 46), (ii) a previously described XMRV clone from the 22Rv1 cell line (n = 1  GenBank: FN692043) (iii) full-length XMRV sequences reportedly amplified from PC (n = 6), or CFS patient samples (n = 2), (iv) previously reported XMRV pol sequences derived from PC patient material (n = 6), (v) additional C57BL/6 endogenous full-length MLV sequences identified using BLAT (n = 28), and (vi) various other MLV complete genomes (n = 5).
Quite unexpectedly, visual inspection of the 2552 nt XMRV pol sequences  revealed that sequences VP29 and VP184, which were apparently amplified from PC patient material, are recombinants of the 22Rv1 cell line virus and Moloney MLV (MoMLV). A nucleotide BLAST search revealed that the Moloney MLV derived fragment from VP29 (1182 nt) is 100% identical to MoMLV (GenBank AF033811), 11 nucleotides different to the closest known mouse endogenous MLV (GenBank AC153360) and 22 nucleotides different to the XMRV clone derived from the 22Rv1 cell line (FN692043). The recombinant nature of VP29 and VP184 was confirmed by phylogenetic incongruence analysis (Additional File 4; Figure S2). The fact that MoMLV envelope does not have tropism for human cells, that these PCR products were derived from human material, and that the 1182 nt XMRV fragment is identical to common MLV-based plasmids, strongly suggests PCR contamination as the source of the recombinant.
22Rv1 associated XMRV is more diverse than patient derived sequences
Another notable characteristic of the XMRV clade is its asymmetry (B1 asymmetry statistic = 24.47, p < 0.002). This is an expected property of families of endogenous mobile elements . Phylogenetic asymmetry implies that whenever replication occurs, one daughter sequence tends to be inactive whilst the other continues to proliferate. This phenomenon arises naturally when one (or a few) active endogenous viruses in a genome generate inactive copies by re-infection , but is difficult to explain under a hypothesis of host-to-host transmission. Extreme cases of strong selection among genetically diverse variants can cause asymmetry , although in this case, the lack of XMRV genetic diversity is incompatible with this possibility.
Whilst our observations cannot conclusively prove that XMRV is not a human pathogen they appear consistent with the hypothesis that XMRV is not an exogenous virus transmitting among individuals. Instead, multiple lines of evidence suggest that the full length clones of XMRV originated from the 22Rv1 cell line. PCR contamination could arise directly from 22Rv1 cells or from cells inadvertently infected with the 22Rv1 derived virus. We speculate that the 22Rv1 cells became infected with XMRV during their passage through athymic mice . Data in Figure 1 demonstrate that mouse DNA could also contaminate patient samples as a variety of mice encode sequences, with endogenous MLV proviruses, that are detected with PCR protocols that are designed to detect XMRV. It is quite possible therefore that previously published findings are explained by contaminated PCR where the patient samples were contaminated by mouse DNA or DNA from cells infected with MLV-X including that from 22Rv1 cells. A recent study amplified polytropic MLV sequences rather than XMRV from chronic fatigue patient samples  and healthy donors. Unfortunately the MLV sequences described there were too short to carry out a thorough phylogenetic analysis, and we have therefore not included them here. It is difficult to retrospectively establish whether prior studies have contaminated patient samples. Importantly, assay contamination cannot be assessed by detection of murine DNA alone since MLVs contaminate a significant proportion of non-murine cell lines common in laboratories [1, 30]. PCR contamination has previously been found to underlie erroneous association between retroviruses and human disease underlining the difficulties associated with detecting pathogens by PCR [31, 32].
We conclude that future screens for MLV-related sequences use more rigorous PCR containment procedures, such as those used to reliably recover ancient DNA , or manage contamination by controlling for its inevitable frequency, for example by screening equal numbers of controls prepared and stored identically, together with test samples . Positive samples must be sequenced and those that are identical to known endogenous murine sequences, or plasmids present in the host laboratory, should be treated with caution. Whilst true association of XMRV with human disease would be of great medical importance, it is imperative that such an association is rigorously established before it impacts on diagnosis and patient care. We suggest that XMRV as a human virus does not conform to this criterion.
PCR of mouse genomic DNA was performed using primers/probes as previously described (Additional File 1; Table S1) [4, 5, 14]. PCR conditions including buffers (a single batch of 2× Taqman PCR master mix (Applied Biosystems)) and thermocycler conditions were also essentially as described [4, 5, 14]. All human tumour cell line Taqman PCRs were run in a duplex assay using the Taqman® RNase P Control Reagents (VIC) (Applied Biosystems) as an internal control. Cycling conditions were 95°C for 15 secs and annealing/extension at 60°C for 1 minute after an initial denaturation of 10 min. Thresholds were routinely set at default values. Primers are shown in Additional File 1; Table S1.
454 Sequencing of inbred and wild-derived inbred mouse samples
Mouse DNA samples were obtained from the Jackson Labs (Bar Harbor, Maine) except for Balb/c which was obtained from Sigma (D4416). 100-200 ng DNA from each mouse was amplified using Platinum Pfx (Invitrogen) proofreading polymerase and primers EG87 and EG89 (Additional File 1; Table S1). 500 ng of amplified DNAs were sequenced using the Genome Sequencer FLX Instrument and GS FLX Titanium series reagents (Roche/454 Life Sciences) according to the manufacturer's instructions. SFF files were processed using the sfffile and sffinfo commands of the SFF tools, split based on the MIDs, and FASTA files were created for each sample. Reads containing XMRV-specific 24-nt deletion, were identified using a customised python script and their frequency was calculated.
PCR, direct sequencing and sequence analysis
Partial gag, pol and env sequences were amplified from genomic DNA from human tumour cell lines using primers labelled TC in Additional File 1; Table S1. PCR products were purified and subjected to direct sequencing, which obviated PCR error, using an Applied Biosystems 3730×l DNA analyzer. Gag, pol and env sequences from the 22Rv1 cell line were amplified with a single stock of Platinum Pfx (Invitrogen) proofreading polymerase and primers labelled 22Rv1 in Additional File 1; Table S1. PCR product was gel purified and ligated into pZero Blunt and transformed (Invitrogen). Positive plasmid clones derived from individual colonies were then sequenced.
22Rv1 cell line derived gag (1605 nt; n = 11), pol (1635 nt; n = 15) and env (1935; n = 8) unique sequences were manually aligned with 6 full-length XMRV sequences apparently amplified from PC samples (GenBank DQ241301, DQ241302, DQ399707, EF185282, FB579966, NC_007815), 2 full-length XMRV sequences from CFS patient samples (GenBank GQ497343, GQ497344), previously described pol sequences derived from PC patient material (VP29, VP79, VP86, VP88, VP90, VP184; n = 6), one previously reported sequence derived from the 22Rv1 tumour cell line (GenBank FN692043) , 46 nonecotropic endogenous MLV sequences, 28 endogenised MLV full-length sequences identified by BLAT search  of the mouse genome using the VP62 XMRV sequence (GenBank DQ399707) as a query, the DG-75 MLV complete genome sequence (GenBank AF221065), one murine type C retrovirus complete genome (GenBank X94150) and one murine AIDS virus provirus complete sequence (GenBank S80082). The complete genome sequences of AKV and Moloney MLV (GenBank J01998 and AF033811) were added as outgroups. Bayesian phylogenies were reconstructed with the software MrBayes version 3.1.2 , under the General Time Reversible model of nucleotide substitution, with proportion of invariable sites and gamma-distributed rate heterogeneity (GTR+I+G). GTR+I+G parameters were estimated with the program PAUP* version 4b10 using full-length genomes only, and were fixed prior to the phylogenetic reconstruction. The Markov chain Monte Carlo (MCMC) search was set to 3,000,000 iterations, with trees sampled every 1000th generation. Maximum clade credibility trees were selected from the posterior distribution with the program TreeAnnotator version 1.5.2 http://beast.bio.ed.ac.uk/, after discarding a 20% burn in. Trees were edited with the program FigTree version 1.1.2 http://tree.bio.ed.ac.uk/software/figtree/. Phylogenetic reconstructions were also conducted in a gene-specific manner for the gag, pol and env loci following the aformentioned methodology. The posterior probability that the ancestor of the XMRV clade was not in the cell line was estimated by recording the number of times a cell line clone was not basal to the clade on a random sample of 1000 trees extracted from the Bayesian posterior distribution.
Tree shape was assessed using the B1 statistic, accounting for phylogenetic uncertainty by marginalizing B1 across a sample of 1000 posterior trees generated by Mr. Bayes (B1 average = 24.27, 95% highest posterior density interval = 20.98-26.75). The null distribution of B1 was calculated by simulating 1000 phylogenetic trees with 55 taxa (i.e. the size of the XMRV clade). Significance was assessed by counting the number of times B1 values less than the mean B1 value occurred in the simulation (2 in 1000 replicates).
The phylogenies of the gag (783 nt; n = 6), pol (381 nt; n = 9) and env (516 nt; n = 9) MLV-X sequence fragments amplified from tumour cell lines were estimated as described above. The sequences were compared to nonecotropic endogenous MLVs (n = 95), AKV (GenBank Acc. J01998) and Moloney MLV (GenBank Acc. AF033811) sequences, as well as 5 XMRV sequences apparently amplified from PC and 2 from CFS samples. The MCMC search was set to 5,000,000 iterations, with trees sampled every 100th generations.
Following visual inspection of the gene sequences derived from PC patient material (n = 15), a 1185 nt pol gene fragment from patients VP29 and VP184 was used in a nucleotide BLAST search against all available sequences. Both fragments showed maximal identity with two Moloney MLV complete genomes (GenBank AF033811 and J02255; 99% and 100% identity with VP29 and VP184 respectively). Further evidence of recombination between XMRV and Moloney MLV was sought by examining phylogenetic incongruence in two maximum likelihood trees based on i) the 1185 nt pol gene fragment (position 2400 to 3585 of the Moloney MLV AF033811) and ii) the following 1335 nt (position 3586 to 4921). The pol sequences used for the analysis comprised 14 patient-derived sequences (see above), one Moloney MLV sequence (GenBank AF033811), one AKV virus sequence (GenBank J01998) and 39 nonecotropic endogenous MLVs. The trees were reconstructed under the GTR+I+G model of evolution, using PAUP*. The robustness of the topologies was assessed by neighbour joining bootstrapping with 1000 replicates.
Genetic distance analysis
Pairwise nucleotide differences in the gag (1605 nt), pol (1635 nt) and env (1935 nt) of the 22Rv1 and patient-associated sequences were calculated using PAUP*. Genetic distances were estimated i) as the uncorrected number of observed nucleotide substitutions per site and ii) under the GTR+I+G model of evolution. Prior to computation, sequences were screened for APOBEC-3G/F mediated G > A hypermutations, using the Hypermut2.0 algorithm from the Los Alamos HIV Sequence Database, and hypermutations masked. A total of 81 and 5 hypermutations were found in the 22Rv1 and patient-associated sequences respectively. The XMRV/Moloney recombinants VP29 and VP184 were excluded. The null hypothesis that genetic diversity is equal in 22Rv1 clones and patient-derived XMRV sequences was tested using the non-parametric Wilcoxon sum rank test.
MLV-X nucleotide sequences obtained from human tumour cell lines have GenBank accession numbers [FR670581-FR670601] and 22Rv1 derived sequences are [HQ385277-HQ385320].
Acknowledgements and Funding
We thank Robert Silverman for XMRV pol sequences from patients VP29, VP79, VP86, VP88, VP90, VP184, Christine Kozak, Torsten Schaller, Yasuhiro Takeuchi and Robin Weiss for advice, Richard Rance, Vicky Murray and Tarryn Porter for generation of 454 sequencing data, Simon Watson for bioinformatic assistance and John Masters and Robert Kypta for the 22Rv1 cell line. This work was funded by the European Community's Seventh Framework Programme (FP7/2007-2013) under the project 'Collaborative HIV and Anti-HIV Drug Resistance Network (CHAIN)', grant agreement no. 223131 (DP, SH), the National Institute of Health Research UCL/UCLH Comprehensive Biomedical Research Centre (ERG, DP, GJT) Wellcome Trust Senior Fellowships WT076608 and WT090940 (GJT), Wellcome Trust Sanger Institute (AG, CJH, SM, AF, PK), the Medical Research Council (GJT, JAG) and The Royal Society (AK, OGP). The funding sources had no role in study design; in the collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.
- Urisman A, Molinaro RJ, Fischer N, Plummer SJ, Casey G, Klein EA, Malathi K, Magi-Galluzzi C, Tubbs RR, Ganem D, et al: Identification of a novel Gammaretrovirus in prostate tumors of patients homozygous for R462Q RNASEL variant. PLoS Pathog. 2006, 2: e25-10.1371/journal.ppat.0020025.PubMed CentralView ArticlePubMedGoogle Scholar
- Lombardi VC, Ruscetti FW, Das Gupta J, Pfost MA, Hagen KS, Peterson DL, Ruscetti SK, Bagni RK, Petrow-Sadowski C, Gold B, et al: Detection of an Infectious Retrovirus, XMRV, in Blood Cells of Patients with Chronic Fatigue Syndrome. Science. 2009, 326: 585-589. 10.1126/science.1179052.View ArticlePubMedGoogle Scholar
- Fischer N, Hellwinkel O, Schulz C, Chun FK, Huland H, Aepfelbacher M, Schlomm T: Prevalence of human gammaretrovirus XMRV in sporadic prostate cancer. J Clin Virol. 2008, 43: 277-283. 10.1016/j.jcv.2008.04.016.View ArticlePubMedGoogle Scholar
- Dong B, Kim S, Hong S, Das Gupta J, Malathi K, Klein EA, Ganem D, Derisi JL, Chow SA, Silverman RH: An infectious retrovirus susceptible to an IFN antiviral pathway from human prostate tumors. Proc Natl Acad Sci USA. 2007, 104: 1655-1660. 10.1073/pnas.0610291104.PubMed CentralView ArticlePubMedGoogle Scholar
- Schlaberg R, Choe DJ, Brown KR, Thaker HM, Singh IR: XMRV is present in malignant prostatic epithelium and is associated with prostate cancer, especially high-grade tumors. Proc Natl Acad Sci USA. 2009, 106: 16351-16356. 10.1073/pnas.0906922106.PubMed CentralView ArticlePubMedGoogle Scholar
- Arnold RS, Makarova NV, Osunkoya AO, Suppiah S, Scott TA, Johnson NA, Bhosle SM, Liotta D, Hunter E, Marshall FF, et al: XMRV infection in patients with prostate cancer: novel serologic assay and correlation with PCR and FISH. Urology. 2010, 75: 755-761. 10.1016/j.urology.2010.01.038.View ArticlePubMedGoogle Scholar
- Baliji S, Liu Q, Kozak CA: Common Inbred Strains of the Laboratory Mouse That Are Susceptible to Infection by Mouse Xenotropic Gammaretroviruses and the Human-Derived Retrovirus XMRV. J Virol. 2010, 84: 12841-12849. 10.1128/JVI.01863-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan Y, Liu Q, Kozak CA: Six host range variants of the xenotropic/polytropic gammaretroviruses define determinants for entry in the XPR1 cell surface receptor. Retrovirology. 2009, 6: 87-10.1186/1742-4690-6-87.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan Y, Liu Q, Wollenberg K, Martin C, Buckler-White A, Kozak CA: Evolution of functional and sequence variants of the mammalian XPR1 receptor for mouse xenotropic gammaretroviruses and the human-derived retrovirus XMRV. J Virol. 2010, 84: 11970-11980. 10.1128/JVI.01549-10.PubMed CentralView ArticlePubMedGoogle Scholar
- van Kuppeveld FJ, de Jong AS, Lanke KH, Verhaegh GW, Melchers WJ, Swanink CM, Bleijenberg G, Netea MG, Galama JM, van der Meer JW: Prevalence of xenotropic murine leukaemia virus-related virus in patients with chronic fatigue syndrome in the Netherlands: retrospective analysis of samples from an established cohort. BMJ. 2010, 340: c1018-10.1136/bmj.c1018.PubMed CentralView ArticlePubMedGoogle Scholar
- Erlwein O, Kaye S, McClure MO, Weber J, Wills G, Collier D, Wessely S, Cleare A: Failure to detect the novel retrovirus XMRV in chronic fatigue syndrome. PLoS One. 2010, 5: e8519-10.1371/journal.pone.0008519.PubMed CentralView ArticlePubMedGoogle Scholar
- Groom HC, Boucherit VC, Makinson K, Randal E, Baptista S, Hagan S, Gow JW, Mattes FM, Breuer J, Kerr JR, et al: Absence of xenotropic murine leukaemia virus-related virus in UK patients with chronic fatigue syndrome. Retrovirology. 2010, 7: 10-10.1186/1742-4690-7-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Hohn O, Krause H, Barbarotto P, Niederstadt L, Beimforde N, Denner J, Miller K, Kurth R, Bannert N: Lack of evidence for xenotropic murine leukemia virus-related virus(XMRV) in German prostate cancer patients. Retrovirology. 2009, 6: 92-10.1186/1742-4690-6-92.PubMed CentralView ArticlePubMedGoogle Scholar
- McCormick AL, Brown RH, Cudkowicz ME, Al-Chalabi A, Garson JA: Quantification of reverse transcriptase in ALS and elimination of a novel retroviral candidate. Neurology. 2008, 70: 278-283. 10.1212/01.wnl.0000297552.13219.b4.View ArticlePubMedGoogle Scholar
- Cornelissen M, Zorgdrager F, Blom P, Jurriaans S, Repping S, van Leeuwen E, Bakker M, Berkhout B, van der Kuyl AC: Lack of detection of XMRV in seminal plasma from HIV-1 infected men in The Netherlands. PLoS One. 2010, 5: e12040-10.1371/journal.pone.0012040.PubMed CentralView ArticlePubMedGoogle Scholar
- Jeziorski E, Foulongne V, Ludwig C, Louhaem D, Chiocchia G, Segondy M, Rodiere M, Sitbon M, Courgnaud V: No evidence for XMRV association in pediatric idiopathic diseases in France. Retrovirology. 2010, 7: 63-10.1186/1742-4690-7-63.PubMed CentralView ArticlePubMedGoogle Scholar
- Switzer WM, Jia H, Hohn O, Zheng H, Tang S, Shankar A, Bannert N, Simmons G, Hendry RM, Falkenberg VR, et al: Absence of evidence of xenotropic murine leukemia virus-related virus infection in persons with chronic fatigue syndrome and healthy controls in the United States. Retrovirology. 2010, 7: 57-10.1186/1742-4690-7-57.PubMed CentralView ArticlePubMedGoogle Scholar
- Sfanos KS, Sauvageot J, Fedor HL, Dick JD, De Marzo AM, Isaacs WB: A molecular analysis of prokaryotic and viral DNA sequences in prostate tissue from patients with prostate cancer indicates the presence of multiple and diverse microorganisms. Prostate. 2008, 68: 306-320. 10.1002/pros.20680.View ArticlePubMedGoogle Scholar
- Beck JA, Lloyd S, Hafezparast M, Lennon-Pierce M, Eppig JT, Festing MF, Fisher EM: Genealogies of mouse inbred strains. Nat Genet. 2000, 24: 23-25. 10.1038/71641.View ArticlePubMedGoogle Scholar
- Raisch KP, Pizzato M, Sun HY, Takeuchi Y, Cashdollar LW, Grossberg SE: Molecular cloning, complete sequence, and biological characterization of a xenotropic murine leukemia virus constitutively released from the human B-lymphoblastoid cell line DG-75. Virology. 2003, 308: 83-91. 10.1016/S0042-6822(02)00074-0.View ArticlePubMedGoogle Scholar
- Takeuchi Y, McClure MO, Pizzato M: Identification of gammaretroviruses constitutively released from cell lines used for human immunodeficiency virus research. J Virol. 2008, 82: 12585-12588. 10.1128/JVI.01726-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Forbes SA, Bhamra G, Bamford S, Dawson E, Kok C, Clements J, Menzies A, Teague JW, Futreal PA, Stratton MR: The Catalogue of Somatic Mutations in Cancer (COSMIC). Curr Protoc Hum Genet. 2008, Chapter 10: Unit 10 11-PubMedGoogle Scholar
- Katzourakis A, Rambaut A, Pybus OG: The evolutionary dynamics of endogenous retroviruses. Trends Microbiol. 2005, 13: 463-468. 10.1016/j.tim.2005.08.004.View ArticlePubMedGoogle Scholar
- Pybus OG, Rambaut A: Evolutionary analysis of the dynamics of viral infectious disease. Nat Rev Genet. 2009, 10: 540-550. 10.1038/nrg2583.View ArticlePubMedGoogle Scholar
- Knouf EC, Metzger MJ, Mitchell PS, Arroyo JD, Chevillet JR, Tewari M, Miller AD: Multiple integrated copies and high-level production of the human retrovirus XMRV (xenotropic murine leukemia virus-related virus) from 22Rv1 prostate carcinoma cells. J Virol. 2009, 83: 7353-7356. 10.1128/JVI.00546-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Paprotka T, Venkatachari NJ, Chaipan C, Burdick R, Delviks-Frankenberry KA, Hu WS, Pathak VK: Inhibition of xenotropic murine leukemia virus-related virus by APOBEC3 proteins and antiviral drugs. J Virol. 2010, 84: 5719-5729. 10.1128/JVI.00134-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Jern P, Stoye JP, Coffin JM: Role of APOBEC3 in genetic diversity among endogenous murine leukemia viruses. PLoS Genet. 2007, 3: 2014-2022. 10.1371/journal.pgen.0030183.View ArticlePubMedGoogle Scholar
- Grenfell BT, Pybus OG, Gog JR, Wood JL, Daly JM, Mumford JA, Holmes EC: Unifying the epidemiological and evolutionary dynamics of pathogens. Science. 2004, 303: 327-332. 10.1126/science.1090727.View ArticlePubMedGoogle Scholar
- Sramkoski RM, Pretlow TG, Giaconia JM, Pretlow TP, Schwartz S, Sy MS, Marengo SR, Rhim JS, Zhang D, Jacobberger JW: A new human prostate carcinoma cell line, 22Rv1. In Vitro Cell Dev Biol Anim. 1999, 35: 403-409. 10.1007/s11626-999-0115-4.View ArticlePubMedGoogle Scholar
- Lo SC, Pripuzova N, Li B, Komaroff AL, Hung GC, Wang R, Alter HJ: Detection of MLV-related virus gene sequences in blood of patients with chronic fatigue syndrome and healthy blood donors. Proc Natl Acad Sci USA. 2010, 107: 15874-15879. 10.1073/pnas.1006901107.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths DJ, Voisset C, Venables PJ, Weiss RA: Novel endogenous retrovirus in rabbits previously reported as human retrovirus 5. J Virol. 2002, 76: 7094-7102. 10.1128/JVI.76.14.7094-7102.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Voisset C, Weiss RA, Griffiths DJ: Human RNA "rumor" viruses: the search for novel human retroviruses in chronic disease. Microbiol Mol Biol Rev. 2008, 72: 157-196. 10.1128/MMBR.00033-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper A, Poinar HN: Ancient DNA: do it right or not at all. Science. 2000, 289: 1139-10.1126/science.289.5482.1139b.View ArticlePubMedGoogle Scholar
- Weiss RA: A cautionary tale of virus and disease. BMC Biology. 2010, 8: 124-10.1186/1741-7007-8-124.PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.PubMed CentralView ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic analysis using parsimony (* and other methods). 1998, Sunderland, MA: Sinauer Associates, 4Google Scholar
- Kirkpatrick M, Slatkin M: Searching for Evolutionary Patterns in the Shape of a Phylogenetic Tree. Evolution. 1993, 47: 1171-1181. 10.2307/2409983.View ArticleGoogle Scholar
- Rose PP, Korber BT: Detecting hypermutations in viral sequences with an emphasis on G --> A hypermutation. Bioinformatics. 2000, 16: 400-401. 10.1093/bioinformatics/16.4.400.View ArticlePubMedGoogle Scholar
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.