Prominent Renal Expression of a Murine Leukemia Retrovirus in Experimental Systemic Lupus Erythematosus
Jessy J. Alexander*,,
Amit K. Saxena*,
Lihua Bao*,
Alexander Jacob*,
Mark Haas and
Richard J. Quigg*
*Section of Nephrology, The University of Chicago, Chicago, Illinois; and Department of Pathology, The Johns Hopkins University, Baltimore, Maryland.
Correspondence to Dr. Richard J. Quigg, The University of Chicago, 5841 S. Maryland Ave., MC5100, Chicago, IL 60637. Phone: 773-702-0757; Fax: 773-702-4816;E-mail: rquigg{at}medicine.uchicago.edu
ABSTRACT. A role for retroviruses in human systemic lupus erythematosus(SLE) and in mouse lupus models such as the New Zealand Blackand White mice (NZB/W) strain has been postulated. This studycompared the gene profile of nephritic NZB/W kidney with nondiseasedNZW controls. The most highly upregulated gene (5.5-fold) hybridizedwith an expressed sequence tag on a cDNA microarray, which wassequenced and found to correspond with an endogenous murineretrovirus related to the Duplan strain (EDV, L08395). NZB/Wkidney contained the full-length 4.2-kb EDV transcript. By 4wk of age in NZB/W mice, an age preceding renal histologic disease,the EDV transcript was more than threefold increased relativeto NZB or NZW control strains. This upregulated expression tendedto fall with progression of renal histologic disease. By insitu hybridization, the EDV transcript was highly expressedin tubules of NZB/W mice. There was also upregulated expressionof EDV transcript in NZB/W lung and brain, sites of inflammationin this strain, but not in spleen or liver. Thus, using microarrays,the most highly expressed gene in mouse lupus nephritis correspondedto an endogenous retrovirus. This retroviral transcript washighly expressed in the kidneys of lupus mice and tended todecline with advancement of disease. The remarkable upregulationof the EDV transcript only in the setting of active diseasesuggests this transcript is involved in inflammatory disease.
Human systemic lupus erythematosus (SLE) is a multiorgan autoimmunedisease characterized by the loss of tolerance to self antigensand the production of autoantibodies (1). Several spontaneousmurine models of human SLE have been developed, beginning withthe identification in 1959 of New Zealand Black mice (2). Themost extensively studied of these models are the female F1 crossbetween New Zealand Black and White mice (NZB/W) and MRL/Mp-Faslpr(MRL/lpr) mice (2,3). Studies of these mice have contributedgreatly to elucidation of SLE pathogenesis. NZB mice exhibitseveral autoimmune features and mild glomerulonephritis (GN)late in life. NZW mice exhibit only minimal autoimmunity, butthey contribute several genes that are critical to the profoundSLE-like disease phenotype found in NZB/W mice (4,5). MRL/lprmice differ from the lupus-prone congenic MRL/+ strain by thenearly complete absence of the membrane Fas protein, necessaryfor apoptosis, which is due to a retroviral insertion in theFas gene (6). As in human SLE, mice of these strains also developelevated levels of IgG autoantibodies to nuclear antigens, includinganti-double stranded DNA antibodies (Abs) (2).
NZB/W mice develop severe proliferative GN (2). The earliestchanges occur before 20 wk of age and include accumulation ofimmune complexes and proliferation within the mesangial region.Later in the course of disease, immune complexes localize inthe peripheral capillary loops, basement membrane thickeningoccurs, and there is proliferation of intrinsic glomerular cellsleading to obliteration of capillary lumina (7). With the severeglomerular injury, crescent formation occurs. Ultimately, glomerulosclerosisoccurs in the terminal phase of disease. NZB/W mice also havetubular basement membrane immune deposits and develop tubulointerstitialnephritis as in human lupus nephritis (2,8).
Although animals of the various lupus strains are thought todie from renal failure (7), this has come into question fromrecent studies investigating the effects of intercellular adhesionmolecule-1 deficiency in MRL/lpr mice. This prevented mortalityapparently by protecting against pulmonary inflammation whilenot affecting the renal disease (9). In addition to pulmonaryinvolvement (9,10), central nervous system pathology also occursin lupus mice, as is true in human SLE (11,12).
A considerable amount of evidence has been accumulated linkingviruses and particularly retroviruses to the pathogenesis ofSLE. Postulated roles for retroviruses in SLE include as intrinsicor mimicry antigens and/or as stimulants of autoimmunity (1315).Among the best-characterized retroviruses and their productsare the inserted retrotransposon in the Fas gene (6) and theendogenous retroviral envelope glycoprotein, gp70, to whichintrinsic Abs are directed. Immune complexes composed of gp70-anti-gp70deposit in glomeruli of mice with SLE (1619). Murineleukemia viruses (MuLV) have also been implicated in human andmurine autoimmune diseases (13,20).
Given the completion of the first draft of the human genomeand the progress made with mouse genes (21,22), the abilityto determine expression of a diversity of genes in differenttissues at various times and under normal and abnormal conditionsis now possible. Traditional approaches, such as Northern analyses,RT-PCR, and RNAse protection assays, are limited in their capacity.Thus, the technique of massively parallel DNA analysis has beendeveloped over the past decade (23). Here, we used this techniqueto screen for relevant genes in murine lupus nephritis. Withthis approach, we identified a MuLV that was highly expressedin murine lupus and was characterized in detail.
Mice
NZB, NZW, and NZB/W mice were obtained from Jackson Laboratories(Bar Harbor, ME). In initial studies, 20 female NZB/W and 15female NZW mice were obtained at 8 wk of age and followed longitudinally.Blood was obtained monthly for determination of blood urea nitrogen(BUN) levels as measured by a Beckman Autoanalyzer (Fullerton,CA). Groups of three animals were sacrificed every 8 wk startingat 16 wk of age to harvest tissue for RNA isolation and histologicstudies. In subsequent studies, animals along with strain controlswere sacrificed at the indicated times.
Histologic Studies
Kidneys were removed and divided into sections snap frozen forimmunofluorescence (IF) microscopy and fixed in 10% bufferedformalin for light microscopic evaluation. Four-micrometer cryostatsections were processed for direct IF microscopy using FITC-conjugatedAbs to mouse C3 and IgG (Cappel Laboratories, Durham, NC).
For light microscopy, 4-µm sections stained with periodicacid-Schiff were provided as coded slides to a renal pathologist(MH) who was blinded to the origin of each section. For eachslide, the extent of GN was graded semiquantitatively on a scalefrom 0 to 4 according to the schema of Passwell et al. (24).
RNA Isolation
Total RNA was isolated from various tissues by a single-stepguanidinium isothiocyanate-phenol-chloroform extraction (TrizolReagent; Life Technologies, Gaithersburg, MD) (25). For microarraystudies, poly(A+) RNA was purified by oligo-dT affinity chromatography(26).
Microarray Studies
One-microgram poly(A+) RNA from NZW and NZB/W renal corticeswere fluorescently labeled with Cy3 and Cy5, respectively, duringthe reverse transcription reaction to produce cDNA. Labeledtargets were then hybridized with 7854 mouse genes and expressedsequence tags (EST) incorporated into Mouse GEM 1 microarrays(Incyte Genomics, Palo Alto, CA). Expression data were normalizedto internal controls to account for variable intensity of thetwo fluorphores and a balanced differential expression ratioprovided for each gene. In 50 genes that are reasonable candidatesas "housekeeping" genes, the balanced differential expressionratio was 1.0 to 1.3, validating this approach.
Reverse Transcriptase-PCR
cDNA was produced from 5 µg of total RNA by RT using oligo-dTprimers (Superscript Preamplification System, Life Technologies).Subsequent PCR was performed in tubes containing 1/200 of thegenerated cDNA, 50 mM KCl, 20 mM Tris-HCl, 2.5 mM MgCl2, 100µM of deoxynucleotide triphosphate, 0.1 µM of eachprimer, and 2.5 U of Taq polymerase. The primers used were:5'-AGAGAGAGCGCAGAAACGTC-3' and 5'-GGCAATGCACGACCTCTCAA-3', spanningbases 2447 to 2750 of the nucleotide sequence for L08395 (discussedfurther below). Twenty-five cycles of a 1-min denaturation at94°C, 1-min annealing at 60°C, and 1-min extension at72°C were performed. In pilot studies, this input of cDNAand number of cycles was limiting in all instances. To quantifythe PCR product synthesized and to control for the integrityof RNA, PCR was simultaneously performed using primers for glyceraldehyde-3-phosphatedehydrogenase (G3PDH) (27). For these reactions, 1/50 of thecDNA generated by RT from 5 µg of RNA input was subjectedto 30 cycles of PCR, which was also in the linear range of productaccumulation. In all cases, negative controls in which RT wasomitted but samples were otherwise handled identically (i.e.,subjected to PCR) were included.
Products generated by RT-PCR were electrophoresed through agarosegels and stained with ethidium bromide. Gels were photographedunder UV light illumination. Photographs were scanned into TIFF,and band intensity was analyzed using NIH IMAGE software (ScionCorp., Frederick, MD). In each instance where comparisons weremade, RT-PCR was performed contemporaneously and products wereelectrophoresed within the same gel.
Kinetic Real-Time PCR
One µg of total RNA from renal cortex was subjected toRT as described above. Subsequently, kinetic PCR was performedwith the following primers: 5'-ACTGGGAGCCGTATACGATG-3' and 5'-GATAGCTGAGGTGGGTGAGC-3',spanning bases 2548 to 2669 of L08395; the labeled probe was5'-[FAM]CCTGCAGGGAAAACCTGT AA-[TAMRA]3' (bases 2587 to 2606of L08395). Included in the same tube were primers and probefor G3PDH; 5'-GGCAAATTCAACGGCACAGT-3', 5'-AGATGGTGATGGGCTTCCC-3',and 5'-[FAM]AAGGCCGAGAATGGGAAGCTTGTCATC-[TAMRA]3'. Reactionswere performed using a Cepheid Smart Cycler System (Sunnyvale,CA). Concentration curves were determined for both L08395 andG3PDH using serially diluted NZB/W renal cortical RNA. Concentrationsof L08395 and G3PDH RNA were calculated from the respectivecurves. The relative concentration of L08395 was obtained bynormalizing with G3PDH for each sample.
Northern Analyses
Twenty-five micrograms of total RNA from renal cortex was subjectedto electrophoresis through a denaturing agarose gel. RNA wastransferred to a nylon membrane by capillary action and crosslinkedto the membrane by UV irradiation. As probe, the 303-bp PCRproduct derived from kidney by using the primers listed abovewas 32P-labeled with a random primer labeling technique (28).Hybridization and washing of membranes were performed underhigh stringency conditions (28). To control for RNA quantityand quality, the gel was visualized under UV light before transferto examine rRNA bands.
In Situ Hybridization In situ hybridization was performed to determine the distributionof IMAGE EST clone 522713 in kidney. This clone in pBluescriptSK- vector was linearized with EcoRI, and labeled riboprobeswere synthesized with T7 (antisense probe) or T3 RNA polymerase(sense probe) using digoxygenin uridine-triphosphate as substrateaccording to manufacturers instructions (Genius RNA labelingkit; Boehringer Mannheim, Mannheim, Germany). Four-micrometercryostat sections from NZB/W and NZW mice were fixed in 4% paraformaldehydein PBS. Subsequently, sections were treated with proteinaseK (10 µg/ml) for 4 min at 37°C. This reaction wasstopped with 0.1 M glycine in PBS for 10 min at room temperature.Sections were dehydrated sequentially with 70%, 95%, and 100%ethanol, air-dried, and then prehybridized in a humidified chamberfor 1 h in hybridization buffer (50% deionized formamide, 4xSSC, 1x Denhardts solution, 500 µg/ml heat-denaturedherring sperm DNA, 250 µg/ml yeast tRNA, and 10% dextransulfate; all from Sigma Aldrich, St. Louis, MO). Sections wererinsed in 2x SSC and incubated overnight at 42°C with 100µl of hybridization buffer containing 300 ng/ml labeledcRNA probe. The following washes were then performed: 2x SSCfor 1 h at room temperature, 1x SSC for 30 min at room temperature(done twice), 0.5x SSC for 30 min at 42°C, and 0.5x SSCfor 30 min at room temperature. Colorimetric detection of thebound labeled cRNA probe was performed using the Genius NonradioactiveNucleic Acid Detection Kit (Boehringer Mannheim). Slides werewashed in buffer 1 (100 mM Tris-HCl, 150 mM NaCl, pH 7.5), blockedwith 2% normal sheep serum and 0.3% Triton X-100 at room temperatureand incubated overnight at 4°C with anti-digoxigenin Abconjugated with alkaline phosphatase (1:500 dilution). Slideswere washed sequentially in buffer 1 and buffer 2 (100 mM Tris-HCl,100 mM NaCl, 50 mM MgCl2, pH 9.5) and then incubated at roomtemperature with developing solution (45 µl nitrobluetetrazoilum salt [75 mg/ml in dimethylformamide], 35 µlX-phosphate solution [5-bromo-4-chloro-3-indolyl phosphate,mono-(p-toluidium) salt, 50 mg/ml in dimethylformamide], and2.4 mg of levamisole [Sigma] in 10 ml of buffer 2). The colorreaction was stopped with buffer 3 (10 mM Tris-HCl, 1 mM EDTA,pH 8). After dehydration in a series of graded ethanols, slideswere mounted with coverslips and observed under the microscope.Serial sections hybridized with sense probes served as controls.
Sequencing
The nucleotide sequence of IMAGE EST clone 522713 in pBluescriptSK- was determined on an ABI 373A DNA Sequencer using the BigDye Terminator Cycle Sequencing kit (Applied Biosystems, FosterCity, CA). The initial round of sequencing was done with T3and T7 primers flanking the insert and thereafter with primersdesigned from the derived sequence. Sequence comparisons weremade with Genetics Computer Group software (Madison, WI).
Statistics Analyses
Statistical analyses were performed with Minitab software (StateCollege, PA). All data are expressed as mean ± SEM. Two-samplet testing and one-way ANOVA followed by Fishers pairwisecomparisons were used. P values < 0.05 were taken as indicatingstatistical significance.
Longitudinal Analysis of Renal Disease in NZB/W Mice
To examine the course of renal disease occurring over time inNZB/W mice, we followed a cohort of NZB/W mice. At times rangingfrom 16 to 52 wk of age, animals were sacrificed for detailedhistologic analyses. An age-matched group of NZW mice was followedcontemporaneously as a control. As expected, animals developedrenal disease starting with mesangial deposition of immune complexesand complement, along with proliferation of cells in the mesangialregion (Figure 1, A through C). With advancing age, immune complexand complement deposits were present in peripheral capillaryloops of glomeruli, which were associated with endocapillarycell proliferation (Figure 1, D through F). In some animals,extracapillary proliferation (glomerular crescents) and fociof tubulointerstitial inflammation were also present. Animalsthat died spontaneously had renal failure (BUN, 166.2 ±54.6 mg/dl; n = 5) as did one that was studied at 52 wk of age(BUN, 198 mg/dl at sacrifice). In this animal, extensive glomerulosclerosis,fibrocellular and fibrous crescents, and tubulointerstitialinflammation were present (Figure 1G).
Figure 1. Histologic progression of renal disease in a longitudinal study of NZB/W mice. Shown are representative photomicrographs of renal cortex from mice at 24 (A through C), 40 (D through F), and 52 (G) wk of age. At 24 wk, there was predominantly mesangial staining for IgG (A) and C3 (B) by IF, with mesangial matrix expansion and mild mesangial hypercellularity by light microscopy (C). At 40 wk, there was strong staining for IgG (D) and C3 (E) in the glomerular capillary loops and mesangium, with focal, weaker C3 staining within arterioles, tubules, and interstitium (E, arrows). Light microscopy showed global mesangial and endocapillary hypercellularity, with focal mild periglomerular inflammation (F). In one animal with renal failure at 52 wk of age (G), there was advanced glomerulosclerosis, with focal fibrocellular crescents (arrow) and interstitial inflammation (asterisk). Original magnifications: x400 in A through D and F; x200 in E and G.
Consistent with their autoimmune background, NZW mice had amodest extent of mesangial IgG and C3 deposition and cellularproliferation at the later ages studied. At all times points,renal function was normal in these mice (BUN, <30 mg/dl).
Identification of a MuLV Transcript in NZB/W Kidney
To examine gene expression changes that might be relevant tolupus nephritis, our initial studies concentrated on activerenal disease. At 32 wk of age, one of three NZB/W animals hadprominent GN (score = 2), while all three control NZW animalshad normal renal morphology. Renal cortical RNA obtained fromthese two groups were used in a single microarray study in whichthe relative transcript abundance of 7854 mouse genes and ESTwas compared. In this type of analysis, balanced differentialexpression ratios greater than 2 are generally accepted as indicatingsignificant differences in gene expression levels (29). Fromthis experiment, some genes shown to be altered were expected,such as the heightened expression of MHC and complement mRNAin NZB/W kidney (24). The balanced differential expression ratiofor the commonly used "housekeeping" gene, G3PDH, was 1.0; assuch, we used this as control for the remainder of the experimentshere.
The gene that was most upregulated in NZB/W renal cortex comparedwith NZW controls was IMAGE EST clone 522713 (AA087673, derivedfrom C57BL/6 mouse skin) with homology to "mouse DNA with endogenousMuLV" sequence. This gene was 5.5-fold upregulated in the NZB/Wmouse renal cortex compared with the NZW controls. In addition,the absolute expression was among the highest of all genes (signalintensity, 20,554 expression "units"). As the sequence datafrom this clone only contained 174 bases, the entire 1665-bpinsert was sequenced and found to correspond to bases 2394 to4058 of a MuLV known as EDV (endogenous sequence related tothe Duplan strain of murine retrovirus; L08395) (30). This transcriptwas further studied.
Expression of EDV Transcript in Lupus Nephritis
A PCR-based approach was used to amplify a 304-base stretchin open reading frame (ORF) A of the EDV transcript. Initialstudies compared the expression with renal histology in NZB/Wmice over the times when animals had active disease (16 to 52wk of age) using renal tissue obtained from the longitudinalstudy described previously. As shown in Figure 2, expressionof the EDV transcript was elevated in NZB/W renal cortex comparedwith that obtained from NZW animals at all time points studied.As disease progressed, the relative expression of the EDV transcriptdecreased.
Figure 2. EDV (endogenous sequence related to the Duplan strain of murine retrovirus) transcripts are upregulated in New Zealand Black and White (NZB/W) mice compared with NZW mice at all times during the course of disease. Mice were sacrificed at the indicated times to determine expression of the EDV transcript and to score the extent of glomerulonephritis (GN). To normalize results, RT-PCR was also performed for G3PDH in parallel. RT-PCR products from individual mice (n = 14 in both NZB/W and NZW groups) are shown in the top panel. In the lower panel, the expression of EDV transcript relative to G3PDH in individual NZB/W () and NZW () mice is shown. In addition, the average GN scores in the groups of mice are provided as mean values ± SEM.
In the preceding studies, age-matched NZW mice were used ascontrols for NZB/W mice with renal disease. To examine whetherthe NZB background in NZB/W mice contributed to EDV transcriptexpression, additional studies were performed with 26-wk-oldNZB mice. As with NZW mice, the EDV transcript was present inNZB kidneys, but at a much lower expression than in NZB/W kidneys.The mean relative expression of EDV/G3PDH in NZB mice was 0.40± 0.21 (n = 5), with no animal having a relative expressionratio 1, while all NZB/W animals of similar age had relativeexpression ratios between 2.0 and 3.5 (cf, Figure 2), illustratingthat it is the unique combination of NZB and NZW genes thatleads to heightened expression of the EDV transcript in kidney.
As shown in Figure 2, elevated expression of the EDV transcriptoccurred by 16 wk of age in NZB/W mice, a time correspondingto the onset of renal inflammatory disease (16). To evaluatewhether the EDV transcript preceded histologic renal abnormalitiesin NZB/W mice, 8-wk-old NZB/W mice were studied and comparedwith age-matched NZB animals. As shown in Figure 3, the EDVtranscript was elevated in NZB/W animals by 8 wk of age comparedwith control NZB animals (or NZW mice, not shown). Similarly,EDV expression was more than threefold increased in 4-wk-oldNZB/W animals compared with age-matched NZW and NZB mice (P< 0.01). Therefore, increased relative expression of theEDV transcript precedes inflammation in NZB/W mouse kidneys.
Figure 3. EDV transcripts are upregulated in NZB/W mice compared with NZB strain controls at 8 wk of age, which precedes renal histologic disease. Shown are RT-PCR products from individual mice. RT-PCR for G3PDH was performed in parallel.
Additional studies were performed to confirm the EDV transcriptwas increased in the kidneys of lupus mice. The relative abundanceof EDV transcript in NZB/W mice during the evolution of renaldisease was confirmed using kinetic RT-PCR. In this case, primersand the probe spanned bases 2548 to 2669 of the EDV transcript.As shown in Figure 4, data using this quantitative approachconfirmed the results shown previously, in that the EDV transcriptpeaked early in the course of disease and was consistently elevatedcompared with strain control NZW mice at all ages.
Figure 4. Kinetic real-time PCR showing EDV transcripts are upregulated in NZB/W mice () compared with NZW strain controls (). Data shown are mean ± SEM of the relative expression of EDV normalized to G3PDH at the indicated ages.
To confirm the presence of the EDV transcript and to determineits size, Northern analyses were performed with renal tissuefrom NZB/W mice. As described previously (30), the EDV transcriptwas 4.2 kb in size (Figure 5). A smaller transcript of 1.4 kbwas also expressed, which followed intensity of the larger transcriptin individual animals (Figure 5). The identity of this smallertranscript was not pursued, although this corresponds to thesize of ORF A of the EDV (30).
Figure 5. Northern analysis for EDV transcript in 12-wk-old NZB/W mice. A 32P-labeled PCR product from nucleotides 2447 to 2750 of EDV was used as probe. The sizes of the two bands identified by this technique are identified with arrows. Shown in the panel below is the gel stained with ethidium bromide and visualized by UV light before transfer.
Localization in Renal Tissue
Since renal cortex has a number of cell types, localizationof the EDV transcript was accomplished by in situ hybridizationstudies. For these, IMAGE clone 522713 (AA087673) was digoxigenin-labeledand hybridized cRNA probe was identified by an immunohistochemicaltechnique. As shown in Figure 6A, EDV mRNA was present in therenal cortex of 8 wk old NZB/W animals primarily in renal tubules.A similar localization was seen in NZB/W mice at 12 and 24 wkof age (Figs. 6B and C). At the peak of expression, there appearedto be modest glomerular expression of EDV mRNA, but this wasconsiderably less than the tubular expression (Fig. 6B). NZWmice had low expression of EDV by in situ hybridization (Fig. 6D),consistent with previous RNA expression data. Hybridizationof NZB/W kidneys with sense strand as probe was negative (notshown).
Figure 6. In situ hybridization for EDV transcript in renal cortices of NZB/W mice at 8 (A), 12 (B), and 24 (C) wk of age and a 24-wk-old NZW animal (D). Staining is primarily within renal tubules, although glomerular staining was also apparent at the peak of expression in NZB/W mice. In all panels, glomeruli are marked by G. Magnifications: x200 in A, C, and D; x400 in B.
Presence of EDV Transcript in Other Organs in Lupus Mice
Additional studies were done to examine the presence of theEDV transcript in other organs in NZB/W mice. Interestingly,in two organs in which inflammation occurs, brain and lung (912),there was elevated expression of EDV in NZB/W mice comparedwith NZW controls (Figure 7). In contrast, there was no statisticaldifference in liver and spleen expression of EDV between NZB/Wand NZW mice. Notably, the relative expression of the EDV transcriptin these organs was considerably lower than that for kidneyin both strains (note that the y-axis scale in Figure 7 is onetenth of those in Figures 2 and 4). Thus, upregulation of theEDV transcript occurs in organs involved in the inflammatoryprocess in this experimental lupus model.
Figure 7. Expression of EDV transcripts in different tissues from 12-wk-old NZB/W and NZW mice. The expression of EDV transcripts relative to G3PDH by RT-PCR in which cycle number was limiting was measured. Data are mean ± SEM (n = 3 each group). *P < 0.05 versus NZW.
In this study, we used an exploratory approach to identify genesaltered in the NZB/W experimental model of lupus nephritis.After following a cohort of NZB/W mice and their control NZWstrain, we chose animals at 32 wk of age to study for informativegene changes. The gene that was the most remarkably upregulated(5.5-fold) came from an EST (IMAGE clone 522713) derived frommouse skin. The nucleotide sequence of this EST provided cluesto its origin, but we needed to sequence the entire clone toestablish its identity. The 1665-bp insert of this clone correspondedto bases 2394 to 4058 of the endogenous MuLV related to theDuplan retrovirus (EDV). This EDV was originally identifiedin C57BL/6 mouse kidney, in which it was expressed to a muchgreater extent than other organs (30). Consistent with thesefindings is that, in addition to NZB/W mice, the EDV transcriptwas highly expressed in renal cortex of NZB and NZW strainsrelative to other organs.
Microarrays have a number of uses in biologic experimentation.They can determine complex expression patterns, such as thetranscript profile as a cell progresses through the cell cycle(31), or classify samples into biologically meaningful groups,such as different tumor classes and their response to chemotherapy(32,33). The bioinformatic tools to accomplish such tasks arealways expanding and include hierarchical and K-mean clustering,self-organizing maps, and principal-component analysis (34).Another utility for microarrays is to define an organismstranscriptome (35). That is what we have attempted here by screeninglupus renal cortical RNA with over 7800 probes and identifyingone EST as being the most highly expressed relative to the appropriatestrain control.
Although our first clue to the presence of the EDV gene in lupusnephritis came from such a microarray study, we used the traditionalapproaches of conventional and kinetic RT-PCR, Northern blotting,and in situ hybridization to thoroughly study this gene in lupusmice. These results clearly show that the EDV transcript washighly expressed in the kidneys of NZB/W mice. Our careful studiescomparing histopathology with gene expression show that theheightened expression of EDV transcript is at its maximum nearthe very onset of inflammatory disease in kidney and tends todecline with the histologic progression of kidney disease. Aswith EDV transcript in kidneys of C57BL/6 mice (30), the full-lengthtranscript in NZB/W mice is 4.2 kb. A smaller transcript of1.4 kb hybridizing with the probe derived from nucleotides 2447to 2750 was also present. This size corresponds with the lengthof ORF A in EDV, and a 1.4-kb transcript was not apparent inNorthern analyses with a probe derived from the more 5' gagregion (30); it is therefore conceivable that this 1.4-kb transcriptis for this ORF A.
Our studies show that both the NZB and NZW parental strainsto NZB/W contain the EDV transcript. However, at all ages studied,NZB/W mice consistently had threefold or greater expressionof EDV transcript compared with age-matched NZB or NZW mice.These studies were performed when animals were as young as 4wk old, an age preceding significant renal inflammation; therefore,these results indicate that heightened expression of the EDVtranscript cannot be due to inflammation per se. These findingsare consistent with some form of epistatic interaction betweenNZB- and NZW-derived loci, as is known to occur in these andother lupus strains in disease phenomena such as autoantibodiesand sera gp70 levels (3638).
As is true in C57BL/6 mice (30), the EDV transcript was expressedin other tissues of lupus mice, although renal expression wasthe greatest among the organs examined. Both lung and brainfrom NZB/W mice had marked upregulation of the EDV transcriptcompared with tissue from NZW mice; as with kidney, these occurredat a time preceding inflammation. Although the kidney is themost prominently involved organ in NZB/W mice (7,39), lung andbrain are also involved in this inflammatory process (7,10,11).In contrast, although the liver and spleen are clearly participantsin the pathophysiology of SLE as constituents of the mononuclearphagocyte and immune systems, they had low-level expressionof EDV transcript, with no significant differences between NZB/Wand NZB mice. This suggests that the EDV transcript is eithera marker for and/or is involved in inflammation.
An etiologic role for retroviruses in human autoimmune diseasesand their mouse models has been advanced by a number of investigators(1315,40). One of the most thoroughly studied is theenvelope glycoprotein gp70 derived from xenogenic endogenousretroviruses (7,18,41). All lupus-prone strains, including NZB,NZW, and MRL/+ mice, have free gp70 protein in sera, the levelsof which are no different from NZB/W and MRL/lpr mice (42).In the latter, gp70 in serum precedes the onset of inflammatoryrenal disease (42). What is unique about lupus mice is thatthey make pathogenic Abs to gp70, and these resultant immunecomplexes deposit in glomeruli (17,39). In contrast to gp70,we here show that EDV transcript clearly segregates into NZB/Wanimals that develop lupus and NZB and NZW strains that aresimply lupus-prone.
Our data show that the expression of the EDV transcript is upregulatedin the kidneys of lupus mice. We have done all that gene expressionprofiling can accomplish by detailing the time course, relativeamount, and location of EDV expression. Subsequent questionsthat are raised by this technique are whether the transcriptis translated, and if so, whether translated viral protein(s)have a role in SLE. By in situ hybridization, EDV transcriptswere present primarily in renal tubules and to a lesser extentin glomeruli. These anatomic sites in kidney are involved inimmune complex deposition, complement activation, and inflammatorycell accumulation, which lead to endocapillary and extracapillaryglomerular proliferation and tubulointerstitial nephritis (2,8,43).As the disease progresses, fibrosis ensues, leading to the pathologicpicture of sclerosing glomerulonephritis together with tubulointerstitialfibrosis and atrophy (WHO class VI in humans).
The exact pathogenesis of lupus nephritis remains undefined,but it certainly involves a combination of effects from immunecomplexes, complement and inflammatory cells (8,4446).Given the markedly increased expression of EDV transcripts inareas involved in renal pathology, it is conceivable that proteinstranslated from these viral transcripts are directly immunogenicand are involved in immune complex deposition in these sitesor that they interact directly with cells of the immune systemto facilitate inflammation. An alternative explanation to linkincreased expression of EDV transcripts with disease concernstranscription of the EDV gene. Possibilities include that transcriptionof EDV occurs through the actions of cis or trans elements directlyrelevant to the pathogenesis of SLE or that the EDV gene promoterfortuitously drives the expression of pathogenic gene(s) involvedin lupus end-organ disease. Thus, our observations that theEDV transcript is markedly upregulated in lupus nephritis coupledwith future studies designed to determine which of these possibilitiesare correct may provide key insights into the pathogenesis ofthis disease.
Acknowledgments
This work was supported by NIH grant R01DK55357 and by a BiomedicalSciences Grant from the Arthritis Foundation. Dr. Bao was supportedby NIH training grant T32DK07510. We thank Dr. V. Michael Holersfor his critical reading of the manuscript.
Hahn BH: Animal models of systemic lupus erythematosus.In: Dubois Lupus Erythematosus, 5th ed., edited by Wallace DJ, Hahn BH, Baltimore, Williams and Wilkins, 1997,pp 339380
Cohen PL, Eisenberg RA: Lpr and gld: Single gene models of systemic autoimmunity and lymphoproliferative disease. Ann Rev Immunol 9: 243269, 1991[CrossRef][Medline]
Vyse TJ, Morel L, Tanner FJ, Wakeland EK, Kotzin BL: Backcross analysis of genes linked to autoantibody production in New Zealand white mice. J Immunol 157: 27192727, 1996[Abstract]
Mohan C, Morel L, Wakeland EK: Genetic insights into murine lupus.In: Lupus. Molecular and Cellular Pathogenesis,edited by Kammer GM, Tsokos GC, Totowa NJ, Humana Press, 1999,pp 124139
Adachi M, Watanabe-Fukunaga R, Nagata S: Aberrant transcription caused by the insertion of an early transposable element in an intron of the Fas antigen gene of lpr mice. Proc Natl Acad Sci USA 90: 17561760, 1993[Abstract/Free Full Text]
Andrews BS, Eisenberg RA, Theofilopoulos AN, Izui S, Wilson CB, McConahey PJ, Murphy ED, Roths JB, Dixon FJ: Spontaneous murine lupus-like syndromes. Clinical and immunopathological manifestations in several strains. J Exp Med 148: 11981215, 1978[Abstract/Free Full Text]
Zoja C, Liu XH, Donadelli R, Abbate M, Testa D, Corna D, Taraboletti, Vecchi A, Dong QG, Rollins BJ, Bertani T, Remuzzi G: Renal expression of monocyte chemoattractant protein-1 in lupus autoimmune mice. J Am Soc Nephrol 8: 720729, 1997[Abstract]
Lloyd CM, Gonzalo JA, Salant DJ, Just J, Gutierrez-Ramos JC: Intercellular adhesion molecule-1 deficiency prolongs survival and protects against the development of pulmonary inflammation during murine lupus. J Clin Invest 100: 963971, 1997[Medline]
Staszak C, Harbeck RJ: Mononuclear-cell pulmonary vasculitis in NZB/W mice. I. Histopathologic evaluation of spontaneously occurring pulmonary infiltrates. Am J Pathol 120: 99105, 1985[Abstract]
Kier AB: Clinical neurology and brain histopathology in NZB/NZW F1 lupus mice. J Comp Pathol 102: 165177, 1990[Medline]
Vogelweid CM, Johnson GC, Besch-Williford CL, Basler J, Walker SE: Inflammatory central nervous system disease in lupus-prone MRL/lpr mice: Comparative histologic and immunohistochemical findings. J Neuroimmunol 35: 8999, 1991[Medline]
Herrmann M, Hagenhofer M, Kalden JR: Retroviruses and systemic lupus erythematosus. Immunol Rev 152: 145156, 1996[CrossRef][Medline]
Denman AM: Systemic lupus erythematosus is a viral aetiology a credible hypothesis? J Infect 40: 229233, 2000[CrossRef][Medline]
Sekigawa I, Ogasawara H, Kaneko H, Hishikawa T, Hashimoto H: Retroviruses and autoimmunity. Intern Med 40: 8086, 2001[Medline]
Yoshiki T, Mellors RC, Strand M, August JT: The viral envelope glycoprotein of murine leukemia virus and the pathogenesis of immune complex glomerulonephritis of New Zealand mice. J Exp Med 140: 10111027, 1974[Abstract]
Izui S, McConahey PJ, Clark JP, Hang LM, Hara I, Dixon FJ: Retroviral gp70 immune complexes in NZBxNZW F2 mice with murine lupus nephritis. J Exp Med 154: 517528, 1981[Abstract/Free Full Text]
Tucker RM, Vyse TJ, Rozzo S, Roark CL, Izui S, Kotzin BL: Genetic control of glycoprotein 70 autoantigen production and its influence on immune complex levels and nephritis in murine lupus. J Immunol 165: 16651672, 2000[Abstract/Free Full Text]
Haywood ME, Vyse TJ, McDermott A, Thompson EM, Ida A, Walport MJ, Izui S, Morley BJ: Autoantigen glycoprotein 70 expression is regulated by a single locus, which acts as a checkpoint for pathogenic anti-glycoprotein 70 autoantibody production and hence for the corresponding development of severe nephritis, in lupus-prone BXSB mice. J Immunol 167: 17281733, 2001[Abstract/Free Full Text]
Rucheton M, Graafland H, Fanton H, Ursule L, Ferrier P, Larsen CJ: Presence of circulating antibodies against gag-gene MuLV proteins in patients with autoimmune connective tissue disorders. Virology 144: 468480, 1985[CrossRef][Medline]
Venter JC, Adams MD, Myers EW, et al: The sequence of the human genome. Science 291: 13041351, 2001[Abstract/Free Full Text]
Broder S, Venter JC: Sequencing the entire genomes of free-living organisms: The foundation of pharmacology in the new millennium. Annu Rev Pharmacol Toxicol 40: 97132, 2000[CrossRef][Medline]
Brown PO, Botstein D: Exploring the new world of the genome with DNA microarrays. Nat Genet 21: 3337, 1999[CrossRef][Medline]
Passwell J, Schreiner GF, Nonaka M, Beuscher HU, Colten HR: Local extrahepatic expression of complement genes C3, factor B, C2 and C4 is increased in murine lupus nephritis. J Clin Invest 82: 16761684, 1988
Quigg RJ, Lo CF, Alexander JJ, Sneed AE, Moxley G: Molecular characterization of rat Crry: Widespread distribution of two alternative forms of Crry mRNA. Immunogenetics 42: 362367, 1995[Medline]
Quigg RJ, Holers VM, Morgan BP, Sneed AE: Crry and CD59 regulate complement in rat glomerular epithelial cells and are inhibited by the nephritogenic antibody of passive Heymann nephritis. J Immunol 154: 34373443, 1995[Abstract]
Foss RD, Guha-Thakurta N, Conran RM, Gutman P: Effects of fixative and fixation time on the extraction and polymerase chain reaction amplification of RNA from paraffin-embedded tissue. Comparison of two housekeeping gene mRNA controls. Diagn Mol Pathol 3: 148155, 1994[Medline]
Quigg RJ, Sneed AE: Molecular characterization of rat glomerular epithelial cell complement receptors. J Am Soc Nephrol 4: 19001908, 1994
Schena M, Shalon D, Heller R, Chai A, Brown PO, Davis RW: Parallel human genome analysis: Microarray-based expression monitoring of 1000 genes. Proc Natl Acad Sci USA 93: 1061410619, 1996[Abstract/Free Full Text]
Kubo Y, Nakagawa Y, Kakimi K, Matsui H, Higo K, Wang L, Kobayashi H, Hirama T, Ishimoto A: Molecular cloning and characterization of a murine AIDS virus-related endogenous transcript expressed in C57BL/6 mice. J Gen Virol 75: 881888, 1994[Abstract/Free Full Text]
Iyer VR, Eisen MB, Ross DT, Schuler G, Moore T, Lee JF, Trent JM, Staudt LM, Hudson JJ, Boguski MS, Lashkari D, Shalon D, Botstein D, Brown PO: The transcriptional program in the response of human fibroblasts to serum. Science 283: 8387, 1999[Abstract/Free Full Text]
Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, Bloomfield CD, Lander ES: Molecular classification of cancer: Class discovery and class prediction by gene expression monitoring. Science 286: 531537, 1999[Abstract/Free Full Text]
Ratain MJ, Relling MV: Gazing into a crystal ball-cancer therapy in the post-genomic era. Nat Med 7: 283285, 2001[CrossRef][Medline]
Mirnics K: Microarrays in brain research: The good, the bad and the ugly. Nat Rev Neurosci 2: 444447, 2001[CrossRef][Medline]
Burge CB: Chipping away at the transcriptome. Nat Genet 27: 232234, 2001[CrossRef][Medline]
Datta SK, McConahey PJ, Manny N, Theofilopoulos AN, Dixon FJ, Schwartz RS: Genetic studies of autoimmunity and retrovirus expression in crosses of New Zealand black mice. II. The viral envelope glycoprotein gp70. J Exp Med 147: 872881, 1978[Abstract/Free Full Text]
Xie S, Chang S, Yang P, Jacob C, Kaliyaperumal A, Datta SK, Mohan C: Genetic contributions of nonautoimmune SWR mice toward lupus nephritis. J Immunol 167: 71417149, 2001[Abstract/Free Full Text]
Wanstrat A, Wakeland E: The genetics of complex autoimmune diseases: Non-MHC susceptibility genes. Nat Immunol 2: 802809, 2001[CrossRef][Medline]
Theofilopoulos AN, Dixon FJ: Models of systemic lupus erythematosus. Adv Immunol 37: 269390, 1985[Medline]
Krieg AM, Steinberg AD: Retroviruses and autoimmunity. J Autoimmun 3: 137166, 1990[CrossRef][Medline]
Izui S, Elder JH, McConahey PJ, Dixon FJ: Identification of retroviral gp70 and anti-gp70 antibodies involved in circulating immune complexes in NZB X NZW mice. J Exp Med 153: 11511160, 1981[Abstract/Free Full Text]
Izui S, McConahey PJ, Theofilopoulos AN, Dixon FJ: Association of circulating retroviral gp70-anti-gp70 immune complexes with murine systemic lupus erythematosus. J Exp Med 149: 10991116, 1979[Abstract/Free Full Text]
Kashgarian M, Hayslett JP: Renal involvement in systemic lupus erythematosus.In: Renal Pathology with Clinical and Functional Correlations,edited by Tisher CC, Brenner BM, Philadelphia, JB Lipincott, 1989,pp 380408
Clynes R, Dumitru C, Ravetch JV: Uncoupling of immune complex formation and kidney damage in autoimmune glomerulonephritis. Science 279: 10521054, 1998[Abstract/Free Full Text]
Wang Y, Hu Q, Madri JA, Rollins SA, Chodera A, Matis LA: Amelioration of lupus-like autoimmune disease in NZB/W F1 mice after treatment with a blocking monoclonal antibody specific for complement component C5. Proc Natl Acad Sci USA 93: 85638568, 1996[Abstract/Free Full Text]
DAgati VD, Appel GB, Estes D, Knowles DM, Pirani CL: Monoclonal antibody identification of infiltrating mononuclear leukocytes in lupus nephritis. Kidney Int 30: 573581, 1986[Medline]
Received for publication January 21, 2002.
Accepted for publication July 30, 2002.