Previous Article | Next Article ![]()
Journal of Virology, November 2002, p. 11738-11743, Vol. 76, No. 22
0022-538X/02/$04.00+0 DOI: 10.1128/JVI.76.22.11738-11743.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
pro-pol Nucleotide Sequences
ApoGene Biotechnologie, D-86567 Hilgertshausen, Germany,1 Institut für Tierzucht und Genetik, Veterinärmedizinische Universität Wien,2 Ludwig-Boltzmann-Institut für Immuno-, Zyto- und Molekulargenetische Forschung Wien, A-1210 Vienna, Austria3
Received 2 July 2002/ Accepted 19 August 2002
| ABSTRACT |
|---|
|
|
|---|
1 family, consisting of the subfamilies A, B, and C. The aim of the study was the concise examination of PERV
by the analysis of the retroviral pro-pol sequences. The analysis of 52 pro-pol clones amplified in this study revealed eight PERV
families. In addition to four already-described families (
1,
4,
5,
6), four novel families (
7,
8,
9,
10) were identified. Quantitative analysis of the novel PERV
sequences in selected breeds revealed variations in the endogenous retroviral load. Open reading frames (ORF) in the amplified proviral fragment were only found for PERV
1. In addition, novel ORF-containing PERV
1 clones consisting of hybrid sequences were revealed. Sequence comparison from published full-length PERV
1 clones of the PERV subfamilies A, B, and C resulted in a lack of strict correlation of the classification of pro-pol and env. The results indicated the occurrence of causative recombination events between retroviral genomes. Thus, our study on PERV
provides new data for the evaluation and selection of pigs intended to be used in xenotransplantation. | TEXT |
|---|
|
|
|---|
Endogenous retroviruses (ERV) are copies of exogenous retroviral genomes integrated into the host genome. They are transmitted vertically to the offspring. ERV have been found in multiple copy numbers in all vertebrates examined (6). Most ERV are defective due to deleterious mutations (24). Whereas retrospective examinations of patients exposed to living porcine tissues have not yielded any indications for the infection of the recipients with PERV (10, 21), the potential risk of cross-species infection of PERV has been shown in in vitro and in vivo experiments (4, 8, 28). Xenotransplantation of tissues derived from pigs which are genetically altered for the aim of the reduction of the serious host-versus-graft reactions in humans (7, 14) has been suggested to enhance the risk of infection with PERV (25, 30).
Referring to the envelope (env) sequences, ERV are classified into the retroviral ß (B- or D-type) and
(C-type) genera (29). PERV A, B, and C, which have been observed to be infectious in vitro, are highly homologous in their gag and pro-pol retroviral genes and, therefore, have been classified as PERV
1 (22). Significant differences in the env gene explain their different host tropisms (16, 26). All pig breeds examined contain PERV A and B with a copy number ranging from 10 to 23 and from 7 to 12, respectively. PERV C is only present in some of the breeds with 8 to 15 copies. Differences in the proviral load have been observed in the genome of different animals, and the correlation of an increased proviral load in highly inbred pig breeds has been discussed previously (1, 5, 13, 15-17, 23). Although most of the copies of PERV A and B are defective, several full-length functional PERV A and B from genomic pig loci have been sequenced (11, 19). Recently, additional PERV
clones, including defective sequences as well as mutant full-length copies (PERV E), have been described (17, 22).
Here we precisely examined the PERV
pro-pol nucleotide sequences. We identified the lack of the strict correlation of the classification of pro-pol and env in already-described PERV
1 sequences as well as novel
1 pro-pol clones harboring hybrid sequences and an open reading frame (ORF).
Amplification and analysis of the PERV sequences.
Previous reports have shown a high PERV load in Landrace pigs (16, 17); therefore, the search for additional PERV
families was carried out with genomic DNA of this breed. pro-pol sequences were amplified by PCR with six pairs of degenerate primers (two 5' primers [9, 27] and three 3' primers [9]) and an annealing temperature of 38°C. The 0.5- to 1.2-kb fragments were separated in agarose gels, isolated, and cloned into the pGEM-T Easy vector (Promega, Madison, Wis.). Thirty independent clones were sequenced from the amplification products of each of the six primer pairs. Sequence analysis was carried out with the Amersham DNA sequencing kit (Amersham Pharmacia Biotech, Vienna, Austria) and the ABI PRISM 377 automated DNA sequencer. The clones were sequenced bidirectionally with additional primers annealing within the cloned sequences. Analysis by BLAST search of the GenBank DNA database revealed that 53 of the 180 clones were of retroviral origin, which showed a length of 0.7 to 1.0 kb. Fifty-two clones were classified as PERV
, one clone as PERV ß, and no clone was grouped as a spumavirus. From the 52 PERV
clones amplified in this study, several 0.9-kb clones showed identical sequences, thereby indicating a low rate of polymerase errors in the PCR process as previously described (15).
Phylogenetic analysis.
By using Gene Jockey II (Biosoft, Cambridge, United Kingdom), further comparison of the complete amplified pro-pol fragment of the 53 clones revealed nine families (eight PERV
families and a single clone belonging to a novel PERV ß family) with 1 to 21 members. Amplified clones showing more than 90% identity were pooled to families. Nomenclature of the families was carried out according to Patience et al. (22). In addition to four already-described PERV
families (
1,
4,
5,
6), four novel families (
7,
8,
9,
10) were found. Thirty-eight of the 52 PERV
clones were assigned to two families,
1 (n = 21) and
6 (n = 17).
The phylogenetic examination of the PERV sequences was carried out both for the DNA and for the protein sequences. Partial and full-length sequences were prepared in SeqApp (http://ftp.bio.indiana.edu/soft/molbio/seqapp/) and were aligned by ClustalW (12), with manual adaptation. Most-parsimony, neighbor-joining, and maximum-likelihood phylogenetic trees were generated by using PAUP* (http://www.lms.si.edu/PAUP) and PHYLIP (http://evolution.genetics.washington.edu/phylip.html). The least defective clone was chosen as the representative of the family. Due to the differing sequences in the 5' region of the pro-pol fragments, 510 nucleotides (nt) of the 3' end and the C-terminal 170 amino acids, respectively, were used according to Herniou et al. (9). Figure 1 shows the tree obtained with the nucleotide sequences by using the neighbor-joining method. The phylogenetic relationship of the sequences was confirmed by analogous results which were generated with the different data sets and the various algorithms used (data not shown).
|
1 GenBank sequences.
The intensively investigated mouse leukemia virus-related PERV
1 family includes the three subfamilies PERV A, B, and C, which are defined by their env gene (1, 16). For the classification of the
1 pro-pol sequences amplified in this study, first-step GenBank sequences of full-length PERV
1 clones were compared. Designation of the clones according to their env gene as PERV A, B, and C was carried out by sequence comparison to Y12238, Y12239, and AF038600, respectively. Classification into pro-pol subfamilies was done according to identities in the polymorphic nucleotide positions of the aligned sequences (Fig. 2). Fourteen of 18 examined 928-nt pro-pol fragments strictly correlated to their env gene, whereas two PERV B (AF038601 and AJ133817) and two PERV C (AX052634 and AX052635) pro-pol sequences correlated with the PERV A env gene. The PERV A pro-pol sequences (AF435966, AF435967, AJ279056, and AJ293656) shared only a few common polymorphic nucleotides and showed an increased sequence polymorphism (Fig. 2). Further analysis of five published 928-nt pro-pol sequences which lack the sequence information for their corresponding env gene resulted in classification as
1B (U77599 and X99933) and
1C (AF033259, AX052636, and U77600). One additional sequence (AF274705) harboring a stop codon showed 40-, 42-, and 48-nt mismatches to
1A (AJ293656),
1B (A66552), and
1C (AF038600), respectively, and therefore remained unclassified (data not shown).
|
1 clones amplified in this study.
The classification of the 21 PERV
1 clones amplified in this study was carried out by sequence comparison of the clones to all published PERV
1 sequences. The pro-pol GenBank sequences showing the highest identity were grouped with the respective clones according to identities in the polymorphic nucleotide positions (Fig. 3). Within the examined 928-nt pro-pol fragment, 3 clones were classified as
1A (Fig. 3A) and 14 clones were assigned to
1B (Fig. 3B).
|
1 clones (m317, m1120, m3111, and m3112) did not match to the published PERV
1 sequences (Fig. 3C). Clone m317 (AF511091) showed differences with
1A (AJ293656),
1B (A66552), and
1C (AF038600) in 60, 29, and 61 nt, respectively. Clones m1120 (AF511090), m3111 (AF511096), and m3112 (AF511097) showed hybrid sequences with high homologies to various pro-pol sequences for different regions. Specifically, clone m1120 harboring an ORF was identical to m317 in its 5' region until nt 285 and was identical to AJ293656 downstream of nt 226, with the exception of one A/G nucleotide exchange each. The two clones m3111 and m3112 showed sequence identity to each other 3' of nt 539, whereas high homologies were found in their 5' regions to m319 and to AF274705, respectively. The overlapping region of the clones m3111 and m3112 (nt 539 to 928) showed no identity to other PERV
1 sequences.
In total, sequence analysis revealed nine clones harboring an ORF in the examined pro-pol sequence which were derived from at least six independent loci. Five of these loci were previously undescribed. In detail, PERV
1A included the ORF containing clone m316 (AF511088). Three PERV
1B clones (m319, m313, and m3117) were identical (AF511092), and clone m2112 (AF511095) was closely related. Clone m3119 (AF511098) maintaining an ORF showed relationship to the mutant sequences m314 and m2115. Clone m3120 differed only slightly from several GenBank sequences represented by A66552. In addition, the hybrid clones m1120 (AF511090) and m3111 (AF511096) harbored an ORF.
Determination of the approximate copy number of the PERV
families was carried by Southern blot analysis (BamHI) with six pig breeds (Duroc, Landrace, Large White, Mangaliza, Pietrain, and Turopolja). Ten micrograms of DNA from three individuals per breed was examined. A specific probe hybridizing to ryr1 was used as a quantitative loading control (31). Probes specific for the different PERV
families were amplified by PCR. Genomic mouse DNA was used as a negative control. PERV clones, which were found to be the closest relatives by phylogenetic analysis, were diluted in genomic mouse DNA and were used to monitor cross-hybridization. For the estimation of the copy number, the signal strength was compared to the signals of the probe used diluted in genomic mouse DNA. For PERV
1, the results were consistent with the data already described (Table 1) (16, 17).
|
4 and
5.
Mutant members have been previously described for both PERV
4 and
5 families (22). Eight amplified clones were representatives of PERV
4. Multiple frameshift and premature stop codon mutations were observed in the clones. About 10 copies per haploid genome were determined for PERV
4 (Table 1). One clone belonging to PERV
5 was observed to harbor multiple missense mutations. The sequence was assigned to a branch with human ERV-9, human endogenous retrovirus W (HERV-W), and MSAV (Fig. 1) and was present in a low copy number in all pig breeds examined (Table 1). The amplified clones of PERV
4 and
5 showed high identity to the published sequences AF274708 and AF274709, respectively.
PERV
6.
Mutant full-length PERV E members have been observed in the pig genome (17). Seventeen of our clones corresponded to this family, which was named PERV
6. Compared to clones of PERV
1, the PERV
6 clones showed a higher sequence diversity. Phylogenetic analysis was done as for PERV
1 and revealed at least 11 novel loci (data not shown). None of the clones was found to have an ORF. Analogous to the PERV
1 hybrid sequences, comparison of the 945-nt clones m116 (AF511100) and m2110 (AF511101) revealed sequence identity 3' of nt 480, whereas the 5' region showed differences in 26 nt (data not shown). The copy number of PERV
6 sequences estimated for different pig breeds (Table 1) was similar to that described in previously published studies (17).
Novel PERV families (
7 to
10 and ß5).
Six clones of endogenous retroviral sequences were found not to match to previously described families. Therefore, they represented PERV
7, consisting of two independent clones, and PERV
8,
9,
10, and ß5. None of the clones was found to have an ORF. Southern blot analysis revealed low copy numbers for the novel PERV
families. In addition, breed-specific differences were detected (Table 1).
The extensive search for new PERV
copies carried out by the identification and phylogenetic analysis of more than 50 independent pro-pol sequences resulted in previously undescribed PERV families as well as novel clones harboring an ORF. Occurrence of several independent 0.9-kb PCR clones in PERV
1B,
4, and
6 showing the identical sequence indicated a low incidence for misincorporation of nucleotides during the PCR procedure as previously described (15). The copy numbers of the different PERV
families observed in Landrace pigs (Table 1) correlated well with the number of clones which were amplified in this study from Landrace genomic DNA and subsequently were assigned to the respective families. Previous studies have suggested the correlation of an increased PERV load in highly inbred breeds (13, 17). However, all breeds examined showed a high proviral load for PERV
1 and
6. Recently, breed-specific differences in the chromosomal integration sites of PERV have been described (15) which were confirmed by our results. This indicates that the use of preselected pigs in subsequent breeding concepts might lead to the reduction of the genomic load of putative infectious PERV.
Within the examined 928-nt pro-pol region, five novel PERV
1 loci harboring an ORF (one
1A, two
1B, and two hybrid clones) were identified. This suggests that additional functional PERV
1 proviruses that have not yet been described may exist in the pig genome. In contrast, in addition to the published mutant full-length PERV E (17), 17 defective clones representing at least 11 novel loci were found for PERV
6. Taking the estimated copy number into account, these results indicate a low incidence of the presence of intact PERV
6 genomes.
After in vitro coculture experiments with porcine cells, human-tropic replication-competent particles with hybrid sequences have been observed in human cells which have not been found in the genome of the donor pigs (20). This suggests that coexpression and recombination events of distinct PERV might lead to the production of recombined virus particles with unknown consequences on the potential risk of infection. In this study we showed the presence of such chimeric PERV
1 sequences in the pig genome, where the classification of pro-pol and env did not correlate within the subfamilies. In humans, chimerism between the reverse transcriptase and transmembrane domains has been found both for endogenous retroviruses and infectious retroviruses (3). Recombination events within the diploid retroviral genome have been described to cause this phenomenon (18). We also found hybrid PERV pro-pol sequences within a given animal. This was confirmed by the observation of hybrid PERV pro-pol sequences (5' region including nt 307) in the published clones AF435967 and AJ279056 (Fig. 2). In addition, intragenic recombinants between PERV
1A and
1B have previously been described for the env gene (15).
Thus, the results of our study indicate the presence of sequences deriving from recombination events between distinct PERV
1 sequences in the pig genome. The potential risk of infection can therefore not be ruled out for mutant PERV loci.
Nucleotide sequence accession numbers.
The nucleotide sequence data of the representatives of the PERV families have been submitted to GenBank for PERV
1 (AF511088-AF511099),
6 (AF511100 to AF511110),
7 (AF511111),
8 (AF511112),
9 (AF511113),
10 (AF511114), and ß5 (AF511115).
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
1,3-galactosyltransferase gene in cloned pigs. Nat. Biotechnol. 20:251-255.[CrossRef][Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Mol. Cell. Biol. | Microbiol. Mol. Biol. Rev. |
|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|