A Multicenter Blinded Analysis Indicates No Association between Chronic Fatigue Syndrome/Myalgic Encephalomyelitis and either Xenotropic Murine Leukemia Virus-Related Virus or Polytropic Murine Leukemia Virus

Study population.A total sample of 293 subjects, 147 patients with CFS and 146 matched healthy control subjects, was recruited from six geographically diverse clinical sites. The sites included the Jen Center of Brigham and Women’s Hospital in Boston, MA; Simmaron Research Institute in Incline Village, NV; Chronic Fatigue & Immune Disorders Research and Treatment Center in Miami, FL; the Levine Clinic in New York City, NY; the Infectious Disease Clinic at Stanford University in Palo Alto, CA; and the Fatigue Consultation Clinic in Salt Lake City, UT. Each of the six clinical sites enrolled between 20 and 30 subjects with matched controls.

Each study subject was assigned a unique study identification (ID) number and a series of linked sample ID numbers. Duplicate aliquots of samples from each subject were labeled with different ID codes (all linked to that subject’s unique study ID) so that laboratory sites would not know whether samples came from the same subject or represented positive controls. Linkage between the study ID and each of the subject-specific sample IDs was retained by the biostatistics team until all laboratory studies were completed and all assay data were reported. All study data and samples received by the coordinating site at Columbia University were deidentified and under code (study ID or sample ID). Personal subject identifiers associated with study ID codes were kept in locked cabinets or on password-protected servers at each study site and were accessible only to key study personnel.

Screening.Potential case and control subjects were screened by each site with a brief telephone questionnaire. Subjects who appeared eligible for the study were scheduled for an evaluation visit at the respective study site.

CFS patients.Clinical site investigators recruited and screened subjects previously diagnosed with CFS for this study. At the study visit, potential cases completed multiple detailed study instruments covering symptoms, past medical history and level of function, and then underwent a physical examination and venipuncture. Based on a study algorithm that used physical examination findings and patient responses on the study instruments, potential case subjects were identified as meeting inclusion criteria for the study if they were between the ages of 18 and 70; met both the 1994 CDC Fukuda criteria for CFS and 2003 Canadian consensus criteria for ME/CFS; reported a viral-like prodrome prior to the onset of CFS, defined as three or more of the clinical features fever, headache, gastrointestinal discomfort/upset, malaise, sore throat, myalgias, arthralgias, and tender lymph nodes; had reduced functional status on at least two of the three subscales of the RAND36 survey (vitality subscale, <35; social functioning subscale, <62.5; role-physical subscale, <50) and <70% on the Karnofsky performance scale of functional impairment; were not pregnant, not <3 months postpartum, and not currently lactating; had none of the exclusionary criteria in the 1994 CDC Fukuda criteria for CFS and 2003 Canadian consensus criteria for ME/CFS case definitions; and had no diagnosis of past or current medical, neurologic, or psychiatric illnesses.

If subjects met the inclusion criteria, 80 ml of venous blood was obtained under each sample ID code and shipped to the coordinating laboratory for processing and division into aliquots. To reduce the potential for diurnal variation, blood samples were drawn between 10 a.m. and 2 p.m. (29). Assays included screens for human immunodeficiency virus infection (HIV), syphilis (rapid plasmin reagin [RPR] test), hypothyroidism (T4 and thyroid-stimulating hormone [TSH] tests), hematocrit, white blood cell count, erythrocyte sedimentation rate, electrolytes, glucose, blood urea nitrogen, creatinine, and transaminases. Potential cases were excluded if screening laboratory values were not within the ranges shown in Table 1 or were HIV or RPR positive.

Healthy controls.Healthy control subjects were recruited at the clinical sites through referral from a case participant or written public advertisement. Control subjects were frequency matched to case subjects by sex and age (within 5 years) as well as by geographic region of residence, season of case subject blood sampling (within 12 weeks), and race/ethnicity (Asian, white, black, Hispanic, or Pacific Islander) (Table 2). Subjects were excluded if they resided in the same household as, or had intimate contact with, a case subject.

Potential control subjects completed the same detailed instruments administered to potential case subjects inquiring about symptoms, medical history, and level of function and underwent a physical examination.

Potential control subjects were excluded if they reported the symptoms or signs of chronic fatigue syndrome as defined by the 1994 CDC Fukuda CFS criteria (21) and 2003 Canadian consensus CFS/ME criteria (20) or if they reported any exclusionary current or past medical, neurologic, or psychiatric illnesses that was not controlled within 5 years of the onset of CFS/ME.

As with potential case subjects, potential control subjects were identified as meeting the inclusion criteria according to a study algorithm based on recorded symptoms, medical history, functional status, and physical examination findings.

If potential control subjects met the inclusion criteria, 80 ml of venous blood was obtained and shipped to the coordinating laboratory, where they were aliquoted and screened. Control subjects who were otherwise eligible were excluded if screening laboratory values were not within the ranges shown in Table 1 or were HIV or RPR positive.

Informed consent.All participants provided informed, written consent in accordance with the protocols approved by the Institutional Review Boards (IRB) for research with human subjects that were associated with each referring institution: Columbia University (AAAI1037), Stanford University (21229), and Brigham and Women’s Hospital (000688). Clinical sites not affiliated with an IRB received approval under the IRB at Columbia University. All participants provided informed, written consent to have their blood tested for HIV in accordance with the requirements of their respective state health departments.

Compensation.All potential case and control subjects were compensated for completing the study visit instruments and questionnaires, undergoing a physical examination, and providing a blood sample.

Instruments.Prior to the study visit, case and control subjects were screened according to a brief screening questionnaire for cases or healthy controls, respectively.

During the study visit, all case and control subjects completed instruments that systematically collected information regarding symptoms, medical history, level of function, quality and intensity of the fatigue, medications, sleep patterns, mood, and pain, including the RAND36 (generic form of the SF-36 instrument) (30), multidimensional fatigue inventory (MFI) (31), symptom inventory, and questionnaires covering demographics, medical history, and comorbid conditions. The physical examination included vital signs and examination of the skin, lymph nodes, head and neck, lungs, heart, abdomen, musculoskeletal system, and nervous system. Additionally, clinical investigators asked subjects to report their level of functioning according to the Karnofsky performance scale adapted for use in CFS (Nancy Klimas, personal communication).

Collection and processing of blood specimens.Fresh EDTA-treated whole-blood specimens were collected from participants on site, deidentified, and shipped overnight at 2 to 4°C to Columbia University for preparation of coded plasma and PBMC aliquots using lymphocyte separation medium. Two distinctly coded aliquots of the same sample were sent in random shipments to each laboratory.

Sample testing.Aliquots of blood specimens received at Columbia University were analyzed for blood chemistries and serologic tests for HIV and syphilis at Columbia University’s Center for Advanced Laboratory Medicine or Quest Diagnostics. Thyroid hormone values were assessed by immunoassay in the CII laboratory (Merck/Millipore). Two distinctly coded aliquots of the same sample were sent in batched shipments for XMRV and pMLV testing to the participating laboratories.

Spiked controls.All cell culture media, cryopreservation reagents, and human samples used as control samples were prescreened for pMLV and XMRV sequences using previously described PCR assays (32). XMRV from 22Rv1 cell culture supernatants and cells were quantified using a previously described PCR assay that can detect a single XMRV copy (32). One-ml coded aliquots of human plasma were spiked with 22Rv1 supernatant containing ~2,000 copies of XMRV RNA. Individual vials of cryopreserved human PBMCs containing 3 million cells were spiked with ~200 22Rv1 cells (2,000 copies/specimen) (33). Spiked plasma samples and PBMC specimens were stored at −80°C until being shipped on dry ice to the testing labs in a blinded fashion with the clinical specimens.

Nucleic acid testing (NAT).The CDC received plasma specimens and performed several NAT assays as described in detail elsewhere (18, 27, 34-36). Briefly, 900 µl of plasma was ultracentrifuged at 45,000 rpm, and RNA was extracted from the pellet using a QIAamp viral RNA minikit (Qiagen). One negative-control human plasma specimen was included with every eleven test plasma specimens during the extraction process. Quantitative real-time RT-PCR assays (qRT-PCR) for generic pMLV/XMRV pro (protease) and gag detection were performed on RNA extracts. The AccessQuick RT-PCR system (Promega) with avian myeloblastosis virus (AMV) RT (Promega), and Tfl DNA polymerase (Promega) and an AgPath one-step RT-PCR kit (ABI/Ambion) with the ArrayScript RT and AmpliTaq Gold DNA polymerase were used for cDNA synthesis and amplification in the pro and gag qRT-PCR assays, respectively. A third PCR was done using the primers XPOLOF and XPOLOR, followed by a nested PCR with the primers XPOLIF and XPOLIR for the generic detection of MLV/XMRV 216-bp pol sequences (18, 27, 34-36). For this reaction, cDNA synthesis and amplification of RNA was done using AMV RT (Promega) and a RobustI RT-PCR kit (Finnzymes). Each PCR experiment included 20 water-only reactions to control for contamination. Positive bands of the correct size in each PCR assay were excised, purified, and sequenced. Detection of viral RNA and not DNA was verified by repeat qRT-PCR and nested-PCR testing of positive specimens with and without RT. To test for mouse DNA contamination, a quantitative real-time PCR assay for mouse intracisternal A particle (IAP) sequences was performed on samples that were PCR positive for XMRV/pMLV sequences (27). Specimens were considered positive if XMRV/pMLV RNA sequences were detected in at least two of three PCR assays and had undetectable IAP sequences or XMRV/pMLV DNA.

The FDA laboratory of Lo and colleagues performed nested RT-PCR and PCR assays as previously described, with some modifications (11). Briefly, 500 µl of plasma was lysed by mixing with 2 ml viral lysis buffer AVL (Qiagen) containing 5 µg carrier RNA, and nucleic acids were purified using a QIAamp viral RNA kit (Qiagen). cDNA synthesis and the first-round PCR were carried out with recombinant Thermus thermophilus (rTth) DNA polymerase (Applied Biosystems) in triplicate using two-thirds of the total nucleic acids extracted from each specimen. For quality assurance, one negative-control plasma sample and one positive-control plasma sample (spiked with ~20 copies/ml of XMRV RNA from a 22Rv1 culture) were included with every 22 test plasma samples during the nucleic acid isolation process. During RT-PCR the following negative and positive controls were tested in triplicate: RNA diluent negative control, 2 copies of XMRV RNA per reaction, 6 copies of XMRV RNA per reaction, and 20 copies of XMRV RNA per reaction. Each assay run was considered valid if one of three XMRV-spiked positive-control plasma samples containing 20 copies/ml tested positive and at least two of three positive controls containing 6 copies of XMRV RNA per reaction tested positive. All negative controls always tested negative. Genomic DNA was extracted from PBMCs using a DNeasy blood and tissue kit (Qiagen). Nested PCR for MLV gag was performed in six replicates using 50 to 60 ng DNA and the 419F/1154R outer primer pair and GAG-I–F/GAG-I–R inner primer set using conditions as previously described (11). Internal positive controls consisted of normal human PBMC spiked with 22Rv1 cells (containing ~10 copies of XMRV genomes per cell). Each test run was considered valid if 2.5 copies of XMRV from 22Rv1 cells could be detected in all 3 triplicates by nested PCR. Any product of approximately the correct size (~746 bp for the outer primer pair and ~410 bp for the internal primer pairs) was excised, purified, and sequenced. PCR for detection of mouse mitochondrial DNA (mtDNA) was performed on all positive samples using a seminested PCR targeting the mouse-specific D-loop sequence as described previously (32). The mtDNA assay detected 3 fg of both mouse spleen DNA and EL4 mouse cell DNA used to spike into 50 ng of human DNA run in parallel as the positive controls. Detection of murine IAP sequences was also performed on all samples using an IAP PCR assay with 45 amplification cycles and a 72°C annealing temperature (37). The IAP PCR assay detected 100 fg of both mouse spleen DNA and EL4 mouse cell DNA used to spike 50 ng of human DNA and run in parallel as the positive controls (11).

PCR assays for the Mikovits/Ruscetti/Hanson group were performed in the Hanson laboratory at Cornell University. RNA and DNA extracts of cultured PBMC (see “Virus culture assays” below) were received on dry ice and stored at −80°C until use. Nucleic acid concentrations were measured on a NanoDrop spectrophotometer (Thermo, Fisher), and RNA was diluted in RNase-free water to 750 ng/14 µl. Fourteen microliters of each RNA sample was used for cDNA synthesis with a SuperScript VILO cDNA synthesis kit (Invitrogen). Precautions against amplicon and mouse DNA contamination, nested PCR with the Gag-O/Gag-I primers, and single-round gagL PCR with 419F/711R primers was performed as described by Lee et al. (38). Any bands of about the expected size were excised, purified with a PureLink quick gel extraction kit (Invitrogen), sequenced, and analyzed using BLASTn. Adequate amounts of RNA were available except for three gag-negative samples, for which less than 750 ng of RNA was available. IAP PCR with sensitivity for detection of 10 fg of mouse DNA was performed on any gag-positive samples to exclude specimen and tissue culture contamination (38, 39). DNA was not available for IAP PCR for six gag-positive cultures. All gag-positive cDNA samples were negative for mouse IAP sequences. PCR assays were considered positive if XMRV/MLV-like gag sequences were detected after amplification with either set of primers. In order for a specimen to be considered positive, both of the duplicate aliquots that were received were required to test positive with at least one set of primers.

Virus culture assays.Virus culture was performed on PBMC specimens by the Mikovits/Ruscetti/Hanson group at the Ruscetti laboratory. Briefly, 3 million PBMCs were stimulated with 200 U/ml interleukin 2 (IL-2) (rather than 20 U, which we observed to result in considerable cell death following thawing cells) and 2 µg phytohemagglutinin in RPMI containing 10% fetal bovine serum (FBS), 1% l-glutamine, and 1% penicillin-streptomycin in a T25 flask for 48 h. To reduce potential for contamination there was no coculture with LNCaP cells (32). Within the first 2 weeks, 10 to 15 million proliferating cells were then treated with 5 µM 5-aza-2-deoxycytidine for 4 to 5 days in order to inhibit DNA methylation, which could silence integrated retroviruses. The cells in two T25 flasks were pelleted and resuspended in Trizol (Invitrogen) for RNA isolation. The contents of two additional flasks were pelleted for DNA isolation with standard Trizol methods.

Serologic assays.At the Ruscetti lab, the Mikovits/Ruscetti/Hanson group performed a flow cytometry-based serologic assay modified slightly from what was previously reported (10). One million murine BaF3ER cells or BaF3ER cells expressing the recombinant spleen focus-forming virus (SFFV) envelope (Env) were incubated with 1:10, 1:100, and 1:1,000 dilutions of patient sera or control antisera for 20 min at room temperature. SFFV and XMRV are highly similar in the Env N terminus and would thus be recognized by cross-reactive antibodies in serologic assays. The 7C10 rat monoclonal antibody recognizes the N terminus of the XMRV/MLV Env and was used as a positive-control serum. XMRV Env blocks the binding of patient plasma to cell surface XMRV. Cells were washed, incubated with phycoerythrin (PE)-conjugated anti-human IgG (1:200) in FBS-PBS medium for 20 min at 4°C, and washed in cold PBS. Cells were then examined using a FACScan (Becton, Dickinson). Nonspecific reactivity to the BaF3ER cells was determined by incubating BaF3ER cells with patient sera diluted 1:10. Instrument settings were adjusted so that the center of the resulting negative histogram was in the center of the first log. Using the BaF3ER cells as a reference, live cells were determined using forward versus side scatter with the resulting plot set as negative and copied and pasted to the other samples in that group. This comparison was used to determine if the plasma was positive based on the level of fluorescence intensity of viable cells. Reactive human plasma does not always titrate with the 1:100 or 1:1,000 dilutions sometimes showing higher anti-human IgG binding to the BaF3ER-SFFV Env cells than the 1:10 dilution. Plasma specimens were determined to be positive if two of three dilutions reacted to the BaF3ER-SFFV Env cells at levels 3-fold that of the BaF3ER control cells. Data analysis was done using FlowJo software (Tree Star, Ashland, OR).

Statistical analysis.Individual test results reported for the two duplicate, coded sample aliquots received at each of the three laboratory sites were compared for concordance. If both aliquots of a set of duplicate samples were classified as leading to positive results in any single laboratory (based on each individual laboratory’s criteria for defining positive and negative results for individual samples), the result for that sample in that laboratory was defined as positive. Discordance between test results for the duplicate aliquots examined at any one laboratory was defined as a negative result for that sample and for that laboratory. The linkage of the pair of samples within each set of duplicate sample aliquots was performed by the biostatistics team after all assays were complete and all laboratory data had been received. Serologic data were examined for group differences by an exact Mantel-Haenszel test, stratified by clinical site, with a nominal alpha level of 0.05.