Previous Article | Next Article ![]()
Journal of Virology, January 2003, p. 142-149, Vol. 77, No. 1
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.1.142-149.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Division of Cellular and Gene Therapies, Center for Biologics Evaluation and Research, Food and Drug Administration, Bethesda, Maryland 20892,1 Department of Immunology and Microbiology and the Karmanos Cancer Institute, Wayne State University, Detroit, Michigan 48201,2 Institute for Monitoring, Data Management, Biometrics, and Medical Writing, Cologne, Germany3
Received 6 August 2002/ Accepted 24 September 2002
|
|
|---|
|
|
|---|
Southern blot analysis suggested that there may be up to 50 copies of gammaretroviral and up to 32 copies of betaretroviral PERV sequences in the pig genome (15, 35, 36). Although it is not known exactly how many of these loci are capable of encoding infectious PERVs, recent data suggest that this number is small, at least in Large White pigs (5, 19, 32). Of even greater relevance to pig-to-human xenografts is whether the integrated PERVs are transcriptionally expressed. Studies which have examined PERV RNA expression have detected constitutive transcription in many different organs including the kidney, heart, spleen, liver, thymus, lymph nodes, and lung (1, 5, 8, 15, 35). While it has not been determined whether these organs also produce infectious virus, there is evidence for constitutive expression in NIH miniature pigs, as infectious virus was repeatedly and readily isolated from plasma (43). It will be critical to understand tissue-specific expression for the identification and breeding of source animals that do not express replication-competent PERVs.
To better understand the regulation of PERV expression in various tissues and cells, we identified regions of the long terminal repeat (LTR) that control viral transcription. Because the regulatory signals required for retroviral transcription, replication, and integration into the host cellular DNA are located within the LTR, we sequenced and performed functional analysis of the LTRs of the infectious gammaretroviral PERVs that have been identified.
|
|
|---|
Derivation of deleted LTR-SEAP constructs. The NIH.3c LTR-SEAP reporter gene plasmid was used as a base for all of the deleted LTRs included in this study. The NIH.3c LTR was amplified from genomic DNA of PERV-NIH-3o (PERV originally isolated from NIH miniature pig peripheral blood mononuclear cells passaged through HEK 293 cells three times) and then subcloned into the LTR-SEAP reporter plasmid as described above. Deletions of the NIH.3c LTR were derived by PCR with the use of different combinations of primers and then subcloned into the pSEAP2 plasmid (Clontech). All deletions were confirmed by nucleotide sequence analysis.
Transcription initiation site mapping. Total RNA was extracted from cell lines expressing PERV-A, -B, -C, and -NIH. Transcription initiation sites were mapped by using 5' rapid amplification of cDNA ends (Life Technologies, Rockville, Md.). Amplified DNA products were agarose gel purified, cloned into the pCR2.1 TOPO vector, and subjected to DNA sequencing.
Cells used for transfections. Human embryonic kidney cells (HEK 293 cells [ATCC CRL 1573]), human cervical epithelial cells (HeLa cells [ATCC CCL-2]), mink lung fibroblasts (ATCC CCL-64), Mus dunni tail fibroblasts (MDTF; obtained from Leonard Evans, Rocky Mountain Laboratories, Hamilton, Mont.), and swine testis cells (ST-IOWA; obtained from Richard Fister, Tufts University, Boston, Mass.) were maintained in Dulbecco's modified Eagle's medium (Gibco, Grand Island, N.Y.) supplemented with 10% fetal bovine serum (HyClone, Logan, Utah), 2 mM glutamine, 1 mM sodium pyruvate, 100 U of penicillin per ml, and 100 µg of streptomycin per ml. All cells were grown in 5% CO2 at 37°C.
LTR-SEAP transient expression assays. Target cells for transfection of LTR-SEAP plasmids were seeded in 12-well dishes 1 day prior to transfections. The approximate numbers of cells per well were 8 x 104 to 10 x 104 for HEK 293, 2 x 105 for ST-IOWA, and 3 x 104 to 4 x 104 for other cell types. The next day, each well of cells was cotransfected with 0.4 µg of SEAP plasmid and 0.3 µg of cytomegalovirus-ß-galactosidase plasmid (a gift of Edward Max's laboratory, Center for Biologics Evaluation and Research, Bethesda, Md.) by using the Effectene reagent kit (Qiagen, Valencia, Calif.) according to the manufacturer's instructions. All constructs were transfected in duplicate or triplicate wells and assayed by at least two independent transfection experiments. Cells were incubated for an additional 72 h posttransfection, at which time the supernatant and cell lysates were collected for measurement of SEAP activity and ß-galactosidase activity, respectively. P values were calculated with the Student two-tailed paired t test with the software algorithm with Microsoft Excel. P values of <0.05 were considered to represent statistically significant differences.
The SEAP activity in 90 µl of cell-free supernatant was assayed by using the reagents provided in the Great EscAPe SEAP kit (Clontech) according to the manufacturer's instructions for the chemiluminescence assay. The signal was detected with a tube luminometer (TD-20/20; Turner Design Luminometer, Sunnyvale, Calif.). The assay results for SEAP activity were normalized to the values obtained for ß-galactosidase activity in order to control for differences in transfection efficiencies.
Nucleotide sequence accession numbers. The GenBank accession numbers for the LTR sequences of PERV-A, -B, -C, -NIH.3c, and -NIH are AF546883, AF546884, AF546885, AF546886, and AF546887, respectively.
|
|
|---|
![]() View larger version (61K): [in a new window] |
FIG. 1. Sequence alignment of PERV-A, -B, -C, -NIH, and -NIH.3c LTRs. U3, R, and U5 are delineated by bold letters at the start of each region. The R region is further identified with parentheses. Tandem repeated sequences are indicated by brackets. IR sequences are underlined at the 5' and 3' ends of the LTRs. A nonconsensus TATA sequence and reverse CAAT sequence are boxed, and a polyadenylation signal in the R region is underlined. The transcription initiation site is indicated with an arrow. The PERV-A LTR is used as the reference sequence, and nucleotides identical to PERV-A in PERV-B, -C, -NIH, and -NIH.3c are shown as dots. Nucleotide differences are identified by letters, and dashes show missing bases.
|
A common feature among the PERV LTRs is the overall length, between 600 and 800 bp (Fig. 1), which corroborates the size for the PERV-A and PERV-B LTRs reported previously (23). The organization of the PERV LTRs is similar to that of other retroviruses with well-defined U3, R, and U5 regions. A 9-bp inverted repeat (IR) sequence at the 5' end of the LTR is identical for all four PERV LTRs. A perfect IR is present at the 3' end of only the PERV-A and -B LTRs. In comparison, the 3' end IR sequence of PERV-NIH and PERV-C differ by two nucleotide substitutions from the 5' end IR. Similar to other gammaretroviruses is the presence of directly repeated sequences in the U3 region which resemble an enhancer element, a nonconsensus TATA box, and an AATAAA polyadenylation signal in the R region. The initiation site for transcription, which was mapped by using the rapid amplification of cDNA ends technique is located at nucleotides (nt) 508, 468, 361, and 530 for the PERV-A, -B, -C, and -NIH LTRs, respectively. The 5'-ATAAAA-3' sequence, which may be a potential TATA box, was located 32 nt upstream of the identified transcription start sites. We also identified a CAAT box in the reverse orientation 5'-ATTG-3', which is located 66 nt upstream of the start site of transcription. The U5 region for all of the PERV LTRs examined is between 73 to 77 bp long. Both PERV-NIH and PERV-C have binding sites for tRNApro, (5'-TGGGGGCTCGTCCGGGAT-3'), whereas PERV-A and -B have a tRNAgly primer binding site (5'-TGGTGCTTTGGCCGGGGA-3') (data not shown).
PERV LTR repeated sequences. The directly repeated sequence in the PERV LTRs consists of 39 bp. This element is present 3.5 times in the PERV-A LTR and 2.5 times in the PERV-B LTR (Fig. 1 and 2A). A similar sequence with a 5-bp difference is repeated in a similar position 1.5 times in the PERV-C, PERV-NIH, and PERV-NIH.3c LTRs. However, between the 1.5 repeats in the PERV-NIH and PERV-C LTRs, there is a 24-bp sequence that resembles a truncated repeat with an 8-bp difference. Also, the PERV-NIH.3c LTR has two additional copies of a 37-bp repeat compared to a single copy of this sequence in the PERV-NIH LTR. The distance of the PERV repeated sequences from the putative promoter sequences we identified (30 nt) is much smaller than that in other gammaretroviruses (17), which have a region of substantial size (approximately 200 bp) and contain sequences important for both transcriptional regulation and the pathogenicity of some murine leukemia viruses (MLVs) (18, 37, 49, 55).
![]() View larger version (23K): [in a new window] |
FIG. 2. (A) Sequence alignment and identification of transcription factor-binding sites of the repeat elements of PERV-A, -B, -C, -NIH, and -NIH.3c. Each repeat element is indicated in brackets. Dashes indicate missing nucleotides. Identical nucleotides are shown as dots, and differences are indicated by letters. (a), 39-bp repeat element; (b), truncated 39-bp repeat in PERV-C, PERV-NIH, and PERV-NIH.3c; (c), 37-bp sequence repeated three times in PERV-NIH.3c. Binding sites for known transcription factors are underlined and identified by name above the sequence motif, except for PBX1 for PERV-C, -NIH, and -NIH.3c, which is identified below the sequence. Identification of the transcription factor-binding sites was made by using the TRANSFAC database. (B) Enhancer sequence and transcription factor-binding sites for Moloney MLV (Mo-MLV), which represents transcription factor-binding sites for a majority of murine leukemia retroviral enhancers. Double lines indicate overlapping binding sites.
|
Expression analysis of PERV LTRs isolated from different PERV strains. In order to functionally assess the importance of the observed nucleotide differences in the LTR sequences derived from different strains of PERV, each LTR was subcloned into an expression construct, which contained the SEAP gene, as described in Materials and Methods. We also compared an LTR containing 4.5 copies of the PERV-A repeated element (A-4.5), which was isolated from virus-infected HEK 293 cells. The relative transcriptional activity of each LTR was determined in each of five different cell lines: two human cell lines, HEK 293 and HeLa, and three nonhuman cell lines, a porcine fetal swine testis (ST-IOWA), murine MDTF, and mink epithelial (CCL-64). The results were normalized to a positive control containing the simian virus 40 (SV40) early promoter and enhancer and are shown in Fig. 3. No statistical differences were observed in transcriptional activity between any of the LTRs in the ST-IOWA, HeLa, MDTF, or mink cell lines. However, the results with HEK 293 cells were strikingly different from those with the other cell lines. The A-4.5 LTR resulted in an approximately threefold-higher level of SEAP expression than the A-3.5 LTR (P value of 0.03). Likewise, the activity observed for A-4.5 was significantly higher than those for all of the other LTRs examined in HEK 293 cells (all P values were <0.003). The PERV-NIH LTR was consistently the least active in all cell lines tested. Of note, the transcriptional activity of all LTRs relative to the SV40 promoter was significantly lower in the porcine cell line than in cells from other species. Relative SEAP values for the PERV-NIH.3c clone were consistently higher than those observed for the PERV-NIH LTR in HEK 293 cells (data not shown). In order to investigate whether the extra copies of the 37-bp sequence were responsible for the higher SEAP activity, the PERV-NIH.3c LTR was chosen for subsequent analysis.
![]() View larger version (21K): [in a new window] |
FIG. 3. Expression analysis of LTR-SEAP plasmids in different cell lines. Reporter constructs encoding the SEAP cDNAs downstream from the LTRs of different PERV strains were transfected into four different cell lines. A-3.5 and A-4.5 are LTRs derived from the PERV-A class and differ by the number of 39-bp repeats (A-3.5 has 3.5 repeats, and has the sequence represented in the alignment shown in Fig. 1; A-4.5 is identical to that sequence, with the exception of having 4.5 repeats). B and C are LTRs from the PERV-B and PERV-C classes, respectively. NIH is the LTR from the PERV-NIH strain. SV40 represents a positive control plasmid with the SV40 promoter-enhancer element upstream of the SEAP cDNA, and basic represents the negative control plasmid with no promoter-enhancer element upstream of the SEAP cDNA. Results shown are the means plus the standard errors of the mean for each of six separate transfection experiments performed with HEK 293 cells; three experiments with MDTF, mink, and ST-IOWA cells; and two experiments performed with HeLa cells.
|
1F through
3F) had no effect on the NIH.3c LTR activity in HEK 293 cells. However, when an additional 71 bp were deleted (Fig. 4, compare
3F and
4F), we detected a statistically significant twofold decrease in transcriptional activity (P value of 0.008), suggesting the presence of regulatory sequences in this region of the LTR that contribute to the activation of transcription. Removal of the a and b repeat sequences and one copy of the c repeat sequence in
5F did not produce a further decrease in transcriptional activity. The 5' deletion that removed the remaining c repeat sequences and the ATAA box (
6F) resulted in a 50-fold decrease in activity, indicating that these sequences are essential for transcriptional activation. A 3' deletion of 51 bp in the U5 region (
1R and
1F1R) had no effect on transcriptional activity. However, deletion of the remaining U5 sequences and 47 bp of the R region (
1F2R and
4F2R) resulted in a significant increase in NIH.3c LTR activity (P < 0.03 for
1F2R compared with
1F and P < 0.001 for
4F2R compared with
4F). This result indicated the presence of a negative regulatory element or elements in these regions of the PERV-NIH.3c LTR.
![]() View larger version (26K): [in a new window] |
FIG. 4. Schematic depiction of PERV-NIH LTR deletions and relative SEAP activity in HEK 293 and ST-IOWA cells. Reporter constructs carrying various portions of the NIH.3c LTR upstream of the SEAP cDNA were derived as described in Materials and Methods and are schematically depicted. The boxes with gray hatch marks or solid gray represent the repeat elements found in U3. The a box indicates the conserved 39-bp repeat, the b box indicates the 24-bp truncated repeat, and the c box indicates the 37 bp repeated three times in PERV-NIH.3c (shown Fig. 2A). The black box indicates the R region. Each reporter construct was transfected into 293 or ST-IOWA cells in at least two independent experiments containing three replicates each. Shown on the right side are the average SEAP activities for each deletion mutant ± the standard errors normalized to those observed for the NIH.3c parental LTR. ND, not done.
|
1F,
1F2R,
4F, and
4F2R constructs in ST-IOWA cells. A comparison of
1F with
4F produced similar results to those seen in HEK 293 cells, which also indicated the presence of transactivating sequences in the U3 region upstream of the repeat elements that are active in ST-IOWA cells. In contrast, the effect of the 3' deletion in ST-IOWA cells differed from that observed in HEK 293 cells;
1F2R had the same activity as
1F while
4F2R resulted in an approximately threefold decrease in activity relative to
4F. These data suggest that the presence of sequences in the R and U5 regions that upregulate rather than repress transcription. However, results from the
1F2R mutant suggest that additional upstream sequences in the U3 region can compensate for the loss of the
2R sequences. However, the presence of both regions did not have an additive effect on transcription since the SEAP activity of mutant
1F, which contains both the upstream and downstream sequences, was equivalent to that of either
4F or
1F2R. |
|
|---|
4F with
5F). Additional copies of repeated sequences have been detected in the LTRs of other gammaretroviruses, which have led to an increase in replication of their viruses in target cells and an enhancement in virus-induced pathogenicity (3, 20). The enhancer elements of murine gammaretroviruses are usually present as directly repeated sequences of 50 to 120 nt in length (24, 26). Extensive work has demonstrated that retroviral enhancers are important determinants of not only viral replication but also disease specificity and the latency of disease induction (7, 10-13, 17, 22, 24, 27, 28, 50). It has been reported that the PERV-B repeated sequences function as an enhancer element by their ability to potentiate transcription in an orientation-independent manner (39). The direct repeats present in the PERV LTRs differ from the highly conserved sequences present in other gammaretroviral enhancers, including those for the gibbon ape leukemia virus (17), which is closely related by sequence to the PERVs (1, 25). This is the first identification of a mammalian gammaretrovirus that differs in this respect from other known gammaretroviruses. Thus, we predict that the transcription factors that may bind these PERV enhancer elements will differ from those which control the transcriptional activation of other gammaretroviruses (Fig. 2). In transient expression assays, we observed that an increase in the copy number of repeated sequences from 2.5 (PERV-B) to 4.5 (PERV-A-4.5) did not result in a concomitant increase in transcriptional activation in the cell lines we examined, except for HEK 293 cells. This result is consistent with the observation by Scheef et al. (39) that there is a threshold effect for the number of repeats that can increase transcriptional activation in only certain cell types. The statistically significant increase in activity in the LTR-driven expression of the PERV-A-4.5 LTR over the PERV-A-3.5 LTR observed only in HEK 293 cells suggests that the transcription factors necessary to interact with elements found within the repeat sequences are present in HEK 293 cells but may not be present or as active in the other cell lines tested.
We observed that the PERV-A, -B, and -C LTRs had similar transcriptional activities in cells derived from different species and representing different tissue types. Although these results suggest that PERV transcription levels may not vary in different cell types from different species once virus infection has occurred, our deletion mutant analysis, which is discussed further below, indicates that different regions of the LTR may be responsible for the regulation of transcription depending on the cell type. In comparison, the transcriptional activity of the PERV-NIH LTR was consistently lower than that of the other PERV LTRs in all of the cell lines examined. Because the direct repeats of the PERV-NIH and PERV-C are identical, these data support the idea that regulatory sequences outside of this region are responsible for the decreased transcriptional activity detectable for the PERV-NIH LTR. We predict that the region most likely to be responsible for this difference in transcription is the U3 region upstream of the direct repeats, which contains practically all of the nucleotide differences between PERV-NIH and PERV-C (Fig. 1). By a transcription factor-binding site search of the TRANSFAC database, we identified potential protein binding sites that are present in this region of the PERV-C but not in the PERV-NIH LTR, and they may be responsible for this difference. These binding sites include those for the transcription factors SOX5, Ets-1, Evi1, GATA, v-Myb, and CEBP.
It is noteworthy that the transcriptional activity of all of the PERV LTRs tested relative to the SV40 early promoter had significantly lower activity in porcine cells than in cells from other animal species. Perhaps this reflects evolutionary pressure to select for LTRs with reduced transcriptional activity once a retrovirus becomes an endogenous gene (9). Alternatively, it may be that the transcriptional activity observed in ST-IOWA cells may not be representative of primary porcine cells.
A comparison of the deletion mutants in HEK 293 and porcine ST-IOWA cells demonstrated that different regions of the LTR can either activate or repress transcription depending on the cell type. Whether this difference between the HEK 293 and ST-IOWA cells is dependent on species- or cell type-specificity or both requires further analysis. Sequences have been identified for other retroviruses in the R and U5 regions that either activate or repress transcription (6, 14, 34, 38, 48, 53).
A comparison of the
3F and
4F deletion mutants in HEK 293 cells demonstrated the presence of sequences in a 68-bp region upstream of the repeated elements that contribute to transcriptional activation. A TRANSFAC analysis of this region revealed potential protein binding sites for the REBP, AML1, AP1, Sp1, and MZF1 transcription factors, some of which may participate in the transactivation that we observed. Our data also demonstrate that the repeated sequences and ATAA nonconsensus TATA box are essential for transcriptional activation in HEK 293 cells, since deletion of these elements reduced activity by 50-fold.
Our nucleotide sequence analysis of the LTRs of different types of PERV has revealed a structural organization of basal promoter elements that is shared with other gammaretroviruses. In contrast, the repeat elements and putative transcription factor-binding sites that are conserved among the PERV LTRs appear to be unique to these LTRs among other gammaretroviruses. In addition, functional studies of LTR deletion mutants have identified regions with the ability to transactivate or repress transcription, including those that can do both depending on the cell type. This type of information will facilitate the identification of LTR regions that are essential for PERV replication in cells of different types and animal species. Further studies elucidating the cellular factors required for transcriptional activation of PERVs will be critical for the evaluation of the zoonotic risks that may be involved in the use of porcine tissues and organs in human xenotransplantation.
We are grateful to C. Patience, J. Stoye, and Y. Takeuchi for continuous generosity in providing PERV-infected cell lines and plasmid clones without which this study would not have been possible. We also thank Oliver Zill for technical assistance. Finally, we thank Andy Dayton and Xu Lai for critical readings of the manuscript.
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»