Open Access Research article

Injury-elicited stressors alter endogenous retrovirus expression in lymphocytes depending on cell type and source lymphoid organ

Kang-Hoon Lee, Debora Lim, Tajia Green, David Greenhalgh and Kiho Cho*

Author Affiliations

Department of Surgery, Burn Research, Shriners Hospitals for Children Northern California, University of California, Davis, Sacramento, CA, 95817, USA

For all author emails, please log on.

BMC Immunology 2013, 14:2  doi:10.1186/1471-2172-14-2

The electronic version of this article is the complete one and can be found online at:

Received:19 April 2012
Accepted:11 December 2012
Published:5 January 2013

© 2013 Lee et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



Murine leukemia virus-type endogenous retroviruses (MuLV-ERVs) constitute ~10% of the mouse genome and are associated with various pathophysiologic processes. In this study, we examined whether MuLV-ERVs’ response to burn-elicited stressors is specific for certain lymphocyte populations and/or locations of lymphoid organ.


B- and T-cells, which were sorted from nine lymphoid organs of C57BL/6J mice after burn, were subjected to MuLV-ERV expression analyses. Overall, the post-burn MuLV-ERV expression pattern was dependent on lymphocyte type, time after injury, location of lymphoid organ, and MuLV-ERV type. For instance, the MuLV-ERV expression in T-cells from the thymus and three cervical lymph nodes decreased at 3 hours post-burn while the expression of some MuLV-ERVs was augmented in B-cells derived from the mesenteric lymph node. The MuLV-ERV U3 sequences population of the burn-24 hours group was less diverse in comparison to the no burn and burn-3 hours groups. In addition, it was apparent that at the 24 hours time point, the U3 populations of B-cells from both no burn and burn groups were less heterogeneous than the T-cells’ U3 populations. Using the U3 sequences, some of which were isolated only from specific experimental groups (B- vs. T-cells; no burn vs. burn), as probes, 51 putative MuLV-ERVs, including 16 full-length proviruses, were mapped followed by characterization of their biologic properties.


MuLV-ERVs’ response to burn-elicited stressors may be differentially controlled depending on lymphocyte type, location of lymphoid organ, MuLV-ERV type, and stress duration.

B-cell; T-cell; Burn; Stress; Murine leukemia virus-type endogenous retrovirus


Proviral sequences of the murine leukemia virus-type endogenous retroviruses (MuLV-ERVs) are genomic footprints reflecting the ancient infection of germ-line cells and are presumed to constitute ~10% of the mouse genome [1]. Although most MuLV-ERV genomes are reported to be defective primarily due to mutation and recombination events, some retain intact proviral genome features, such as identical long terminal repeats (LTRs) at both ends, which are formed by retroviruses’ unique replication process, and coding potential for the major structural genes (gag, pol and env) which are essential for virion assembly [2-4]. Expression of MuLV-ERVs is mostly controlled by their U3 promoters located on the 5 LTR, in conjunction with the transcription machinery that is dynamically and uniquely formed within individual cells [5].

Recent studies provide evidence that human endogenous retroviruses (HERVs) are involved in normal physiology as well as various disease processes [6]. The HERV-W envelope protein, called syncytin, is highly expressed in the syncytiotrophoblast layer and plays a critical role in placental differentiation [7]. In addition, the inflammatory properties of syncytin contribute to the degeneration of oligodendrocytes responsible for the development of multiple sclerosis, an autoimmune disease [8].

Patients with severe burn injury often succumb to systemic immune disorder in conjunction with internal organ damage [9]. However, the mechanisms underlying the post-burn pathologic changes have not yet been fully characterized. The results from recent studies suggest that post-burn systemic inflammatory response is associated with the activation and/or depletion of specific lymphocyte subpopulations, such as CD4+CD25+ T regulatory cells, lymphocytes expressing CD25, CD69 and CD71, NK-T-cells, and CD8+ T-cells [10-12]. Sepsis-elicited stressors induce apoptosis in the spleen, resulting in the depletion of B-cells and accompanying immune malfunction [13]. However, the central signaling molecules responsible for the post-burn activation and/or depletion of lymphocytes remain to be elucidated.

Our recent studies demonstrated that burn-elicited stressors differentially alter the expression of MuLV-ERVs in a range of internal organs, including various lymphoid tissues, in a ERV type- and organ-specific manner [14]. Moreover, some MuLV-ERV envelope proteins were able to modulate the mRNA expression of inflammatory cytokines such as IL-6 and IL-1β in macrophage cells [15]. In this study, we hypothesize that the differential post-burn alterations in MuLV-ERV expression, depending on lymphocyte types and/or lymphoid organs, may be associated with systemic pathogeneses, including immune disorder. We examined whether MuLV-ERVs’ response to burn-elicited stressors is specific depending on lymphocyte type and the location of lymphoid organ followed by an investigation of the biologic characteristics of the relevant MuLV-ERVs.


Animal experiments

The experimental protocols involving mice were approved by the Animal Use and Care Administrative Advisory Committee of the University of California, Davis. Female C57BL/6J mice (12 weeks old) were purchased from Jackson Laboratory-West (Sacramento, CA). The baseline MuLV-ERV expression patterns were examined in nine lymphoid organs (spleen, thymus, and seven lymph nodes [axillary, inguinal, mesenteric and pooled cervicals {mandibular, accessory mandibular, superficial parotid, and deep cervical}]). Three mice were sacrificed for tissue collection by cervical dislocation without any treatment. The burn protocol has been described previously [16]. Briefly, under general anesthesia, an ~18% total body surface area flame burn was generated on the shaved back of mice followed by resuscitation with saline and buprenorphine. Control (no burn) mice were treated the same as burn mice but without the flame burn. The three mice from each experimental group were sacrificed by CO2 inhalation at 3 hours and 24 hours after burn followed by the collection of the same nine lymphoid organs as above without any pooling. The burn experiment was repeated four times.

Cell sorting

The collected tissues were minced, suspended in a trypsin/EDTA solution, and then passed through a cell strainer (BD Biosciences, CA). The cells were washed with PBS containing 0.1% BSA followed by centrifugation at 1,000 ×g for 5 min at 4°C. Following the removal of red blood cells, cells were sorted into B- and T-cells using pan-B and pan-T Dynabeads (Invitrogen, Carlsbad, CA), respectively, according to the protocols provided by the manufacturer.

Total RNA isolation and RT-PCR

Total RNAs were isolated from the sorted cells using an RNeasy Mini Kit (Qiagen, Valencia, CA) following homogenization using a QIAshredder (Qiagen). cDNAs were synthesized using Sensiscript reverse transcriptase (Qiagen), RNase inhibitor (Promega, Madison, WI) and an oligo-dT primer (5-GGC CAC GCG TCG ACT AGT ACT TTT TTT TTT TTT TTT T- 3). A set of primers, ERV-U1 (5-CGG GCG ACT CAG TCT ATC GG-3) and ERV-U2 (5-CAG TAT CAC CAA CTC AAA TC-3), was used to amplify the MuLV-ERV U3 regions [17]. β-actin, serving as a normalization control, was amplified with the primer set: 5-CCA ACT GGG ACG ACA TGG AA-3 and 5-GTA GAT GGG CAC AGT GTG GG-3.

Measurement of relative purity of sorted cell populations

The relative purity of the sorted cell populations was measured by real-time RT-PCR using CD20 (B-cell antigen) primers (5-AAA ACC TCC AGG AAG AGT TTG GTC-3 and 5-CGA TCT CAT TTT CCA CTG GCA AG-3) for B-cells and Thy1 (Thymus cell antigen-1) primers (5-TCC TGC TCT CAG TCT TGC-3 and 5-AGT AGT CGC CCT CAT CCT TG-3) for T-cells. Real-time RT-PCR was performed using the Brilliant-III SYBR Green QPCR kit from Agilent Technologies (La Jolla, CA) and the results were analyzed using the MxPro (v.4.10) program (Agilent Technologies). Cell type specific mRNA levels (CD20 and Thy1) were calculated from each sorted cell population with the level of β-actin as a normalization control, using the same primer set as above. The relative purity of each cell population was calculated with the assumption that there are only B- and T-cells in each sorted cell population and CD20 and Thy1 are exclusively expressed in B-cells and T-cells, respectively.

Densitometry and statistical analysis

The density of the PCR products was measured using KODAK Molecular Imaging Software Version 4.5.1 (Eastman Kodak Company, Rochester, NY). All optical densities were measured based on mass values with the Hi-Lo™ DNA marker (a gift from Minnesota Molecular Inc., Minneapolis, MN). The density of each MuLV-ERV U3 band was normalized to its corresponding β-actin value to obtain the relative intensity of each U3 band. Differences in the relative intensity of the matching MuLV-ERV U3 bands between no burn and burn were calculated as fold changes. Statistical significance (p < 0.05) of fold changes between the matching pairs (no burn and burn) of the individual U3 bands was determined using the data from four burn experiments and Student’s t-test.

Cloning and sequencing

Eighty-eight MuLV-ERV U3 PCR amplicons/bands, which had contrasting densities between individual no burn and burn pairs, were selected for sequence analyses. The PCR amplicons were purified using the QIAquick Gel Extraction kit (Qiagen) and then cloned into the pGEM®-T Easy vector (Promega). Three clones were prepared from each U3 PCR amplicon/band. The plasmid DNAs were prepared using a QIAprep Spin Miniprep kit (Qiagen) for sequencing analysis. Sequencing was performed at a core sequencing facility at the University of California, Davis.

Alignment and phylogenetic analyses

From a total of 260 U3 sequences, 100 unique U3 sequences were identified using Seqman Pro Version 8.1.2 (DNASTAR, Madison, WI). Multiple sequence alignment of the 100 unique U3s was performed using the Clustal W method in MegAlign (DNASTAR), and a phylogenetic tree was generated by the HYPERTREE program [18]. Seven U3 sequences, representing the seven phylogenetic branches, were subjected to a multiple alignment analysis to examine the U3 sequence features of direct repeats (1/1*, 4/4*, 5/5*, and 6/6*), insertion of 190 nucleotides, and unique region (2) based on the protocol by Tomonaga et al.[17,19].

Profiling of transcription regulatory elements

The 100 unique MuLV-ERV U3 sequences were surveyed for putative transcription regulatory elements using the MatInspector program (Genomatix, Munich, Germany) within the vertebrate matrix group with a core similarity higher than 0.95 [20].

Determination of population diversity of U3 sequences within each experimental group

Using the MEGA4 program [21], the population diversity of U3 sequences was measured in each experimental group and compared among different groups in two analyses: a) no burn, burn-3 hours, and burn-24 hours, and b) no burn-B-cells, burn-24 hours-B-cells, no burn-T-cells, and burn-24 hours-T-cells.

In silico mapping of putative MuLV-ERVs and characterization of their biological properties

The C57BL/6J mouse genome database from the National Center for Biotechnology Information (NCBI) was surveyed for putative MuLV-ERVs with a greater than 98% identity using the 100 unique U3 sequences as mining probes and the Advanced BLAST program (NCBI). The coding potential of each MuLV-ERV was determined using the Vector NTI program (Invitrogen). Three polypeptide genes (gag, pol, and env) were noted as intact (+), partial (P), or defective (−) relative to the references (AY219567.2 and AF033811). Each MuLV-ERV was mapped with the chromosomal location and strand orientation. The primer binding site (PBS) of each MuLV-ERV was identified by examination of a stretch of 18 nucleotides immediately downstream of the 5 LTR region [22]. In addition, the annotated gene(s) within 100 Kb upstream and downstream of each MuLV-ERV locus, if any, were mapped.

Results and discussion

Lymphocyte type-specific MuLV-ERV expression patterns in lymphoid organs of normal mice

Prior to the investigation of post-burn alterations in MuLV-ERV expression in B- and T-cells of different lymphoid organs (spleen, thymus, axillary lymph node [LN], inguinal LN, mesenteric LN, and cervical LNs), baseline MuLV-ERV expression levels were examined in B- and T-cells isolated from ~12 week old female C57BL/6J mice. The relative purity of the sorted cell populations was determined by measuring CD20 (for B-cells) and Thy1 (for T-cells) levels by real-time RT-PCR. The relative purity of the sorted B- and T-cell populations was greater than 93%, except for the B-cell population in the thymus (~65%). Throughout the lymphoid organs examined, the MuLV-ERV expression pattern, which was examined by RT-PCR amplification of the 3 U3 regions, was substantially different between B- and T-cells (Figure 1). The expression pattern of an intense ~500 bp and two (~750 bp and ~450 bp) lighter U3 bands was shared by B-cells derived from all lymphoid organs examined. In contrast, four distinct U3 bands, ranging from ~500 bp to ~750 bp, were amplified in T-cells with variable intensity depending on the source lymphoid organs. It was interesting to note that the ~500 bp U3 band was more intense in B-cells from all four LN samples compared to B-cells from the spleen and thymus. The findings from this study suggest that the expression of MuLV-ERVs is differentially regulated depending on the lymphocyte type, source lymphoid organ, and ERV type.

thumbnailFigure 1. MuLV-ERV expression in B- and T-cells derived from lymphoid organs of C57BL/6J mice. MuLV-ERV expression was examined in B- and T-cells isolated from the spleen, thymus, and seven different lymph nodes of female C57BL/6J mice without any treatment. A schematic drawing of a typical MuLV-ERV proviral structure indicates the relative locations of a set of primers used for amplification of the U3 regions. Differential MuLV-ERV expression patterns were observed in B- and T-cells derived from various lymphoid organs. LN (lymph node), SM (size marker).

Burn-elicited alterations in MuLV-ERV expression in lymphocytes depending on lymphocyte type, source lymphoid organ, and ERV type

The effects of burn-elicited stressors on the expression of MuLV-ERVs in the B- and T-cells from the spleen, thymus, and seven lymph nodes were examined at 3 hours and 24 hours after injury by amplifying the 3 U3 regions. At 3 hours post-burn, the MuLV-ERV expression, represented by four distinct U3 bands, was decreased in the T-cells from the thymus and three cervical LNs (mandibular, accessory mandibular, and superficial parotid), but not in the T-cells from the deep cervical LN (Figure 2). In contrast, there was an increase in intensity of one to two U3 bands in the T-cells from the thymus and superficial parotid cervical LN at 24 hours post-burn when there was a decrease in U3 bands in the T-cells of the spleen and mesenteric LN.

thumbnailFigure 2. Post-burn alterations in the expression of MuLV-ERVs in B- and T-cells isolated from lymphoid organs of C57BL/6J mice. Post-burn (3 hours and 24 hours) changes in the expression of MuLV-ERVs were measured in the B- and T-cells isolated from nine different lymphoid organs. The results, which were presented as a fold change (burn/no burn) for the individual MuLV-ERV U3 bands (a, b, and c for B-cells; d, e, f, and g for T-cells), revealed a unique pattern for each experimental group (cell/organ/time). SPL (spleen), THY (thymus), LN (lymph node), AXI (axillary LN), ING (inguinal LN), MES (mesenteric LN), MAN (mandibular cervical LN), ACC (accessory mandibular cervical LN), SUP (superficial parotid cervical LN), DEE (deep cervical LN). A black arrow indicates unmeasured samples due to the low expression level. Error bars indicate standard error of the mean. * P < 0.05.

In B-cells, there was an induction of two (a and b bands) of the three U3 bands in the mesenteric LN at 3 hours post-burn while two bands (a and c) were down regulated in the B-cells of the mandibular cervical LN. However, at 24 hours post-burn, the expression levels of all three U3 bands were increased in the B-cells from the mandibular cervical LN. Also at the 24 hours time point, post-burn increases in one or two U3 bands were also observed in the B-cells of the spleen, thymus, axillary LN, inguinal LN, and/or accessory mandibular cervical LN.

The results from this study suggest that the lymphocyte type (B- or T-cell), source lymphoid organ, and duration after injury play a role in post-burn changes (induction or repression) in the expression of MuLV-ERVs. It is possible that the post-burn differential expression of MuLV-ERVs resulted from changes in the composition of the B- and T-cell sub-populations (e.g. naïve vs. effector cells), which have different transcription environments, instead of (and/or in addition to) the affected MuLV-ERVs’ altered transcriptional potentials. It was interesting to note that the post-burn increase in the MuLV-ERV expression in the B-cells from the mesenteric LN parallels the proliferation of specific B-cell populations in the mesenteric LN and spleen after burn [23].

MuLV-ERV U3 sequences identified in specific lymphocyte type and/or injury (no burn/burn) condition

Three clones from each of the 88 MuLV-ERV U3 amplicons/bands, which showed contrasting expression levels between the individual no burn and burn pairs, were sequenced to identify lymphocyte type- and/or injury (no burn/burn)-specific U3 sequences. From a total of 260 MuLV-ERV U3 sequences cloned from the entire set of MuLV-ERV U3 bands examined, 100 unique U3 sequences, ranging from 346 nucleotides to 615 nucleotides in size, were identified (Additional file 1: Table S1). An initial phylogenetic analysis resulted in seven distinct branches of U3 sequences, primarily segregated by their size (group I [largest U3s] ~ group VII [smallest U3s]) (Figure 3). Group III (440 bp length) had the most number of U3 sequences (40) followed by group VII (346 bp length) (25 sequences) then group I (greater than 603 bp length) (19 sequences). There were only two U3 sequences in groups IV and VI (433 bp and 361 bp in size, respectively). The group II U3 sequences (5) were isolated only from T-cells and the group VI U3 sequences (2) were unique for B-cells. On the other hand, the group IV U3 sequences (2) were identified only from no burn control mice. The other four branches (groups I, III, V, and VII) consisted of U3 sequences isolated from both lymphocyte types derived from both no burn and burn mice.

Additional file 1: Table S1. Sources of the 100 MuLV-ERV U3 sequences analyzed in this study. Yellow and gray highlights indicate unique MuLV-ERV U3 sequences and clones without sequence information, respectively. MAN (mandibular cervical LN), ACC (accessory mandibular cervical LN), SUP (superficial parotid cervical LN), DEE (deep cervical LN) LN (lymph node), SM (size marker).

Format: PDF Size: 66KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 3. Phylogenetic relatedness of 100 unique MuLV-ERV U3 sequences. Phylogenetic relatedness of the 100 unique MuLV-ERV U3 sequences was determined by a multiple alignment analysis. The MuLV-ERV U3 sequences in each branch were analyzed for their presence in a specific lymphocyte type (B- or T-cell) and/or stress condition (no burn or burn). Seven branches were formed, primarily based on the size of the U3 sequences. Number of sequences is in parentheses.

The results obtained from this study suggest that some MuLV-ERVs are selectively expressed depending on lymphocyte type (B- or T-cell) and/or presence of burn-elicited stressors. However, it is possible that an expanded saturation cloning protocol would find that some cell type and/or burn-specific U3 sequences are not exclusively present in one group over the other. In this case, a statistical evaluation may be needed to determine the specificity of the U3 sequences in regard to the lymphocyte type and/or burn-elicited stressors.

Characteristics of seven U3 sequences representing the phylogenetic branches

Seven representative U3 sequences, which were chosen from the seven main phylogenetic branches, were subjected to an alignment analysis using reference sequences to characterize their sequence features, such as direct repeats, unique sequence, and an insertion of 190 nucleotides. These U3 sequence features were utilized as key parameters for the in silico determination of the MuLV-ERV tropism traits [19]. The U3 sequences of the groups I and II had an insertion of 190 nucleotides, which is consistent with their relatively large size compared to the other five branches (Table 1). There was a well-conserved TATA box in all seven U3 sequences; however, each U3 sequence retained a unique profile in regard to the direct repeats (1/1*, 4/4*, 5/5*, and 6/6*) and the unique sequence (2). The specific features of the individual U3 sequences were translated into the tropism traits of the respective putative MuLV-ERVs [17,19]. However, the sequences features of the group VI U3 sequences were too diverse to predict its tropism trait.

Table 1. Sequence characteristics and putative tropism traits of seven MuLV-ERV U3 sequences representing the seven main phylogenetic branches

Evidence for post-burn increase in the expression of specific MuLV-ERVs

The population diversity of the U3 sequences (a total of 260) was measured in each experimental group and compared among the groups in two different analyses [21]: i) no burn, burn-3 hours, and burn-24 hours, and ii) no burn-B-cells, burn-24 hours-B-cells, no burn-T-cells, and burn-24 hours-T-cells. The results from the first diversity analysis revealed that the U3 population of the burn-24 hours group was substantially less heterogeneous compared to the other two groups (no burn and burn-3 hours) (Figure 4). In contrast, it appears that there was an increase in diversity among the U3 population in the burn-3 hours group. The finding of decreased U3 population diversity within the burn-24 hours group was confirmed by its negative score for the Tajima’s Neutrality test in comparison to the positive values of the other two groups [24]. The results from the second diversity analysis demonstrated that there is a decrease in U3 population diversity in both burn-24 hours-B-cells and burn-24 hours-T-cells groups in comparison to the no burn controls. Furthermore, it was apparent that at the 24 hours time point, the U3 population in B-cells from both no burn and burn were less diverse in comparison to T-cells’ U3 population.

thumbnailFigure 4. Decreased population diversity of MuLV-ERV U3 sequences at 24 hours after burn. The U3 population diversity within each experimental group was measured [21] and the relevant values, indicating different diversity levels, were presented as nucleotide diversity and Tajima statistics scores. There was a marked decrease in the U3 population diversity at 24 hours after burn compared to no burn.

The findings from these diversity analyses indicate that burn-elicited stressors lead to the induction of specific MuLV-ERVs in a time after injury-dependent manner, and the population of MuLV-ERVs expressed in T-cells is more heterogeneous than the population in B-cells.

Differential transcription potentials of MuLV-ERV U3-promoter sequences

The U3 sequence on the 5 LTR, which is presumed to be identical to the U3 of the 3 LTR, serves as an ERV promoter. In this study, the transcription potentials of the 100 unique U3-promoter sequences were surveyed for transcription regulatory elements (TREs). Nine TREs were predominantly present in the U3 sequences isolated from both B- and T-cells of burn mice at the 3 hours and/or 24 hours time points compared to the U3 sequences from no burn mice (Table 2 and Additional file 2: Table S2). Among the nine TREs are BRNF (Brn POU domain factors), AIRE (Autoimmune regulatory element binding factor), and NRF1 (Nuclear respiratory factor 1), which have been reported to be linked to various stress responses [25,26]. Furthermore, it appears that three TREs occurred more often in the U3 sequences of B-cells than T-cells whereas the U3 sequences from T-cells had a higher frequency of 16 other TREs in comparison to B-cells (Table 2 and Additional file 2: Table S2). Among the 16 TREs, MAZF (Myc associated zinc fingers), GATA (GATA binding factors) and MEF2 (myocyte-specific enhancer binding factor) have been reported to be involved in T-cell activation [27-29]. A total of 122 TREs, including the 28 TREs described above, were identified from the 100 unique sequences, and the comprehensive TRE profiles were summarized in Additional file 2: Table S2 and Additional file 3: Table S3.

Additional file 2: Table S2. The profile of transcription regulatory elements in 100 unique MuLV-ERV U3-promoter sequences. The number in the box indicates the occurrence frequency of each transcription regulatory element (TRE). Colored boxes indicate TREs predominantly identified in the U3 sequences derived from specific experimental groups: burn group (green), B-cell group (yellow), and T-cell group (blue).

Format: PDF Size: 88KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 3: Table S3. Abbreviations: transcription regulatory elements.

Format: PDF Size: 13KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Table 2. Differential transcription potentials of MuLV-ERV U3-promoter sequences

It is probable that the TREs (e.g., BRNF, AIRE, NRF1), which were mapped on the U3 sequences almost exclusively derived from burn mice, are closely associated with the modulation of MuLV-ERV expression in conjunction with a transcription environment dynamically formulated by burn-elicited stressors.

In silico cloning and characterization of putative MuLV-ERVs harboring the unique U3 sequences

To investigate the biological properties of the putative MuLV-ERVs harboring the U3 sequences identified in this study, the C57BL/6J mouse genome was surveyed to map the proviral loci using the 100 unique U3 sequences as mining probes. From the in silico mining experiment, 49 ERVs were found with two identical LTR sequences at both ends and two, MuERV-49 and MuERV-26, retained slightly different LTR sequences due to a few point mutations. Forty-nine of the 51 putative MuLV-ERVs identified from this survey, including 15 full-length proviruses retaining a set of intact coding potentials for gag, pol and env polypeptides, have been reported in previous studies [14,30-33]. The two newly identified MuLV-ERVs are: one full-length provirus (8,759 nucleotides) with intact coding potentials and one defective provirus (6,846 nucleotides) (Table 3). An analysis of the primer binding site revealed that both MuLV-ERVs had the tRNAGln priming sequence. All 51 putative MuLV-ERVs were mapped using the U3 probes derived from both no burn and burn mice, except for BM-a-2.11e [31], which was identified using a U3 probe unique for B-cells of burn mice. On the other hand, 30 putative MuLV-ERV loci were retrieved using U3 probes exclusively isolated from T-cells. Future studies focusing on the biological properties of these putative MuLV-ERVs would be of interest: 1) effects (e.g., cytotoxicity) of their replication on immune cells in conjunction with infection followed by random integration into the genome and 2) functions (e.g., regulation of inflammatory mediators) of their gene products (gag, pol, and env polypeptides).

Table 3. Genomic map and characteristics of putative MuLV-ERVs harboring the unique U3 sequences

Identification of genes neighboring the MuLV-ERV loci

The U3 promoter and the enhancer sequences embedded in MuLV-ERVs can affect the transcriptional activities of the neighboring genes via the promotion/repression of transcription, and alternative splicing and polyadenylation [34]. In this study, we identified the genes located within 100 Kb upstream and downstream of each MuLV-ERV integration site. Of a total of 131 genes identified, nine genes had MuLV-ERVs integrated into their introns (Additional file 4: Table S4). It will be interesting to examine how the expression patterns (depending on lymphocyte type or burn-elicited stress) of these proviruses are correlated with the expression of the relevant genes, for example, the relationship between BM-a-2.11e MuLV-ERV, which was identified using a U3 probe uniquely identified from the B-cells of burn mice, and Camk2b (calcium/calmodulin-dependent protein kinase II-beta) locus. The Camk2b expression is associated with T-cell activation and cytotoxic T-cell activity [35,36]. It is possible that the burn-elicited stressors may affect the promoter and/or enhancer activities of certain MuLV-ERVs paralleling modulation of the expression of adjacent genes in a specific lymphocyte type.

Additional file 4: Table S4. Annotated genes neighboring 100 Kb upstream and downstream of the individual MuLV-ERV integration loci. Chr (Chromosome number), Ori (Orientation). Gray highlight indicates the annotated genes harboring an ERV in the intron. *MuLV-ERV newly identified in this study. Please refer to Table 3 for the MuLV-ERV names.

Format: PDF Size: 120KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data


It is expected that burn injury elicits a range of stress signals that are directly linked to pathologic changes in the immune system and other organs [37]. The majority of the underlying mechanisms regarding how these stressors contribute to post-burn pathogenesis still remain to be characterized. The findings from this study suggest that certain MuLV-ERVs are differentially modulated depending on the lymphocyte type, location of lymphoid organ, and duration after injury in conjunction with the burn-elicited stressors. As a result of activation of certain MuLV-ERVs following burn injury, the gene products and/or newly produced virus particles of the activated ERVs may participate in specific burn pathogeneses, such as immune disorder and cytotoxicity, in a cell type-specific manner. In fact, our laboratory has reported that some MuLV-ERVs are activated in the liver of mice after burn injury which resulted in evident hepatic damage [33]. Further investigations are needed to elucidate the link between lymphocyte type-specific post-burn MuLV-ERV activation and systemic inflammatory disorder as well as the other aspects of burn pathogeneses.


HERV: Human endogenous retroviruses; LN: Lymph node; LTR: Long terminal repeat; MuLV-ERVs: Murine leukemia virus-type endogenous retroviruses; PBS: Primer binding site; TREs: Transcription regulatory elements.

Competing interests

There is no conflicting financial interest.

Authors’ contributions

This study was conceived and managed by KC. DGG participated in scientific discussions. KHL, DL, and TG performed the experiments. KHL generated the figures and drafted the manuscript. All authors read and approved the final manuscript.


This study was supported by grants from Shriners of North America (No. 86800 to KC and No. 84302 to KHL [postdoctoral fellowship]) and the National Institutes of Health (R01 GM071360 to KC).


  1. Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, et al.: Initial sequencing and comparative analysis of the mouse genome.

    Nature 2002, 420:520-562. PubMed Abstract | Publisher Full Text OpenURL

  2. Mayer J, Sauter M, Racz A, Scherer D, Mueller-Lantzsch N, Meese E: An almost-intact human endogenous retrovirus K on human chromosome 7.

    Nat Genet 1999, 21:257-258. PubMed Abstract | Publisher Full Text OpenURL

  3. Mueller-Lantzsch N, Sauter M, Weiskircher A, Kramer K, Best B, Buck M, Grasser F: Human endogenous retroviral element K 10 (HERV-K 10) encodes a full-length gag homologous 73-kDa protein and a functional protease.

    AIDS Res Hum Retroviruses 1993, 9:343-350. PubMed Abstract | Publisher Full Text OpenURL

  4. Vogt VM: Retroviral Virions and Genomes. In Retroviruses. Edited by Coffin JM, Hughes SH, Varmus HE. Cold Spring Harbor (NY): Cold Spring Harbor Laboratory Press; 1997:27-69. OpenURL

  5. Lavie L, Kitova M, Maldener E, Meese E, Mayer J: CpG methylation directly regulates transcriptional activity of the human endogenous retrovirus family HERV-K (HML-2).

    J Virol 2005, 79:876-883. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Moyes DL, Martin A, Sawcer S, Temperton N, Worthington J, Griffiths DJ, Venables PJ: The distribution of the endogenous retroviruses HERV-K113 and HERV-K115 in health and disease.

    Genomics 2005, 86:337-341. PubMed Abstract | Publisher Full Text OpenURL

  7. Frendo JL, Olivier D, Cheynet V, Blond JL, Bouton O, Vidaud M, Rabreau M, Evain-Brion D, Mallet F: Direct involvement of HERV-W Env glycoprotein in human trophoblast cell fusion and differentiation.

    Mol Cell Biol 2003, 23:3566-3574. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Mameli G, Astone V, Arru G, Marconi S, Lovato L, Serra C, Sotgiu S, Bonetti B, Dolei A: Brains and peripheral blood mononuclear cells of multiple sclerosis (MS) patients hyperexpress MS-associated retrovirus/HERV-W endogenous retrovirus, but not Human herpesvirus 6.

    J Gen Virol 2007, 88:264. PubMed Abstract | Publisher Full Text OpenURL

  9. Purcell EM, Dolan SM, Kriynovich S, Mannick JA, Lederer JA: Burn injury induces an early activation response by lymph node CD4+ T-cells.

    Shock 2006, 25:135. PubMed Abstract | Publisher Full Text OpenURL

  10. Schneider DF, Glenn CH, Faunce DE: Innate lymphocyte subsets and their immunoregulatory roles in burn injury and sepsis.

    J Burn Care Res 2007, 28:365-379. PubMed Abstract | Publisher Full Text OpenURL

  11. Maldonado MD, Venturoli A, Franco A, Nunez-Roldan A: Specific changes in peripheral blood lymphocyte phenotype from burn patients. Probable origin of the thermal injury-related lymphocytopenia.

    Burns 1991, 17:188-192. PubMed Abstract | Publisher Full Text OpenURL

  12. Barlow Y: T lymphocytes and immunosuppression in the burned patient: a review.

    Burns 1994, 20:487-490. PubMed Abstract | Publisher Full Text OpenURL

  13. Hotchkiss RS, Tinsley KW, Swanson PE, Schmieg RE, Hui JJ, Chang KC, Osborne DF, Freeman BD, Cobb JP, Buchman TG: Sepsis-induced apoptosis causes progressive profound depletion of B and CD4+ T lymphocytes in humans 1.

    J Immunol 2001, 166:6952-6963. PubMed Abstract | Publisher Full Text OpenURL

  14. Lee YK, Chew A, Phan H, Greenhalgh DG, Cho K: Genome-wide expression profiles of endogenous retroviruses in lymphoid tissues and their biological properties.

    Virology 2008, 373:263-273. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Lee YK, Chew A, Greenhalgh DG, Cho K: Tropism, cytotoxicity, and inflammatory properties of two envelope genes of murine leukemia virus type-endogenous retroviruses of C57BL/6J mice.

    Mediators Inflamm 2011, 2011:509604. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Cho K, Adamson LK, Greenhalgh DG: Induction of murine AIDS virus-related sequences after burn injury.

    J Surg Res 2002, 104:53-62. PubMed Abstract | Publisher Full Text OpenURL

  17. Tomonaga K, Coffin JM: Structures of endogenous nonecotropic murine leukemia virus (MLV) long terminal repeats in wild mice: implication for evolution of MLVs.

    J Virol 1999, 73:4327-4340. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Bingham J, Sudarsanam S: Visualizing large hierarchical clusters in hyperbolic space.

    Bioinformatics 2000, 16:660-661. PubMed Abstract | Publisher Full Text OpenURL

  19. Tomonaga K, Coffin JM: Structure and distribution of endogenous nonecotropic murine leukemia viruses in wild mice.

    J Virol 1998, 72:8289-8300. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Quandt K, Frech K, Karas H, Wingender E, Werner T: MatInd and MatInspector: new fast and versatile tools for detection of consensus matches in nucleotide sequence data.

    Nucleic Acids Res 1995, 23:4878-4884. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0.

    Mol Biol Evol 2007, 24:1596-1599. PubMed Abstract | Publisher Full Text OpenURL

  22. Nikbakht KN, Ou CY, Boone LR, Glover PL, Yang WK: Nucleotide sequence analysis of endogenous murine leukemia virus-related proviral clones reveals primer-binding sites for glutamine tRNA.

    J Virol 1985, 54:889-893. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Kotadia BK, Ravindranath TM, Choudhry MA, Haque F, Al-Ghoul W, Sayeed MM: Effects of pentoxyfylline on mesenteric lymph node T-cells in a rat model of thermal injury.

    Shock 2003, 20:517-520. PubMed Abstract | Publisher Full Text OpenURL

  24. Tajima F: DNA polymorphism in a subdivided population: the expected number of segregating sites in the two-subpopulation model.

    Genetics 1989, 123:229-240. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Kovacs KJ, Sawchenko PE: Sequence of stress-induced alterations in indices of synaptic and transcriptional activation in parvocellular neurosecretory neurons.

    J Neurosci 1996, 16:262-273. PubMed Abstract | Publisher Full Text OpenURL

  26. Suliman HB, Carraway MS, Welty-Wolf KE, Whorton AR, Piantadosi CA: Lipopolysaccharide stimulates mitochondrial biogenesis via activation of nuclear respiratory factor-1.

    J Biol Chem 2003, 278:41510-41518. PubMed Abstract | Publisher Full Text OpenURL

  27. Duncan DD, Stupakoff A, Hedrick SM, Marcu KB, Siu G: A Myc-associated zinc finger protein binding site is one of four important functional regions in the CD4 promoter.

    Mol Cell Biol 1995, 15:3179-3186. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Ting CN, Olson MC, Barton KP, Leiden JM: Transcription factor GATA-3 is required for development of the T-cell lineage.

    Nature 1996, 384:474-478. PubMed Abstract | Publisher Full Text OpenURL

  29. Volker JL, Rameh LE, Zhu Q, DeCaprio J, Hansen U: Mitogenic stimulation of resting T-cells causes rapid phosphorylation of the transcription factor LSF and increased DNA-binding activity.

    Genes Dev 1997, 11:1435-1446. PubMed Abstract | Publisher Full Text OpenURL

  30. Kwon DN, Lee YK, Greenhalgh DG, Cho K: Lipopolysaccharide stress induces cell-type specific production of murine leukemia virus type-endogenous retroviral virions in primary lymphoid cells.

    J Gen Virol 2011, 92:292-300. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Kwon DN, Greenhalgh DG, Cho K: Cloning and characterization of endogenous retroviruses associated with postinjury stress signals in lymphoid tissues.

    Shock 2009, 32:80. PubMed Abstract | Publisher Full Text OpenURL

  32. Stoye JP, Kaushik N, Jeremiah S, Best S: Genetic map of the region surrounding the retrovirus restriction locus, Fv1, on mouse chromosome 4.

    Mamm Genome 1995, 6:31-36. PubMed Abstract | Publisher Full Text OpenURL

  33. Lee YK, Chew A, Fitzsimon L, Thomas R, Greenhalgh D, Cho K: Genome-wide changes in expression profile of murine endogenous retroviruses (MuERVs) in distant organs after burn injury.

    BMC Genomics 2007, 8:440. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  34. Selten G, Cuypers HT, Berns A: Proviral activation of the putative oncogene Pim-1 in MuLV induced T-cell lymphomas.

    EMBO J 1985, 4:1793. PubMed Abstract | PubMed Central Full Text OpenURL

  35. Mack KD, Von Goetz M, Lin M, Venegas M, Barnhart J, Lu Y, Lamar B, Stull R, Silvin C, Owings P, et al.: Functional identification of kinases essential for T-cell activation through a genetic suppression screen.

    Immunol Lett 2005, 96:129-145. PubMed Abstract | Publisher Full Text OpenURL

  36. Liu X, Yao M, Li N, Wang C, Zheng Y, Cao X: CaMKII promotes TLR-triggered proinflammatory cytokine and type I interferon production by directly binding and activating TAK1 and IRF3 in macrophages.

    Blood 2008, 112:4961-4970. PubMed Abstract | Publisher Full Text OpenURL

  37. Deitch EA: Role of the gut lymphatic system in multiple organ failure.

    Curr Opin Crit Care 2001, 7:92. PubMed Abstract | Publisher Full Text OpenURL