Analysis of transcribed human endogenous retrovirus W env loci clarifies the origin of multiple sclerosis-associated retrovirus envsequences
© Laufer et al; licensee BioMed Central Ltd. 2009
Received: 14 January 2009
Accepted: 15 April 2009
Published: 15 April 2009
Multiple sclerosis-associated retrovirus (MSRV) RNA sequences have been detected in patients with multiple sclerosis (MS) and are related to the multi-copy human endogenous retrovirus family type W (HERV-W). Only one HERV-W locus (ERVWE1) codes for a complete HERV-W Env protein (Syncytin-1). Syncytin-1 and the putative MSRV Env protein have been involved in the pathogenesis of MS. The origin of MSRV and its precise relation to HERV-W were hitherto unknown.
By mapping HERV-W env cDNA sequences (n = 332) from peripheral blood mononuclear cells of patients with MS and healthy controls onto individual genomic HERV-W env elements, we identified seven transcribed HERV-W env loci in these cells, including ERVWE1. Transcriptional activity of individual HERV-W env elements did not significantly differ between patients with MS and controls. Remarkably, almost 30% of HERV-W env cDNAs were recombined sequences that most likely arose in vitro between transcripts from different HERV-W env elements. Re-analysis of published MSRV env sequences revealed that all of them can be explained as originating from genomic HERV-W env loci or recombinations among them. In particular, a MSRV env clone previously used for the generation of monoclonal antibody 6A2B2, detecting an antigen in MS brain lesions, appears to be derived from a HERV-W env locus on chromosome Xq22.3. This locus harbors a long open reading frame for an N-terminally truncated HERV-W Env protein.
Our data clarify the origin of MSRV env sequences, have important implications for the status of MSRV, and open the possibility that a protein encoded by a HERV-W env element on chromosome Xq22.3 may be expressed in MS brain lesions.
Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease of the central nervous system thought to result from an as yet incompletely understood complex interplay of genetic and environmental factors . Stimulated by the finding that human T cell leukemia virus type 1 (HTLV-1), a human exogenous retrovirus, causes a neurological disease (HTLV-1 associated myelopathy/tropical spastic paraparesis) with similarities to MS , retroviruses have also been searched for in MS. This led to the detection of retrovirus-like particles containing reverse transcriptase activity in the supernatants from cultured cell of patients with MS [3–5]. After molecular characterization of retroviral RNA sequences within such particles, this novel retroviral element was named MS-associated retrovirus (MSRV) [6, 7]. Subsequent probing of human genomic DNA with MSRV sequences revealed endogenous retroviral sequences closely related to MSRV, the human endogenous retrovirus (HERV) family type W (HERV-W) . HERVs are remnants of ancestral germ line infections by active retroviruses, which have thereafter been transmitted in a Mendelian manner. In humans, HERVs comprise approximately 8% of the genome (for review see [9, 10]). They typically consist of an internal region containing gag, pro, pol, and env genes, flanked by two long terminal repeats (LTR). HERV-W is a multicopy family consisting of ~650 elements dispersed in the human genome . About 280 of those elements contain internal sequences. The remaining elements lack internal regions because of recombinational deletion between the two LTRs, leaving a solo LTR . Like almost all HERV families, HERV-W is highly defective due to acquisition of stop-codons, frameshift mutations, and deletions. In addition, many members of the HERV-W family represent processed pseudogenes that were generated through retrotransposition by long interspersed elements (LINE) [11–13]. Notably, no replication-competent HERV-W provirus could be identified so far [8, 14]. Yet, a single HERV-W env locus (ERVWE1, chromosomal location 7q21.2) retained a complete open reading frame (ORF) for a functional envelope (Env) protein, termed Syncytin-1 . Syncytin-1 is highly expressed in the placenta where it likely participates in the fusion of cytotrophoblast cells into the syncytiotrophoblast layer . The ERVWE1 locus therefore appears to have been diverted into a bona fide human gene .
Intriguingly, in addition to the physiological function of Syncytin-1 in placental morphogenesis, several studies have provided evidence in support of a possible pathogenic role of HERV-W Env/Syncytin-1 in MS. HERV-W env/Syncytin-1 RNA levels were found to be higher in autopsied brain tissue from patients with MS than in brain tissue from controls [18–21]. Neuropathological investigations reported increased expression of HERV-W Env/Syncytin-1 protein in astrocytes and microglia in actively demyelinating brain lesions from patients with MS [18, 21, 22]. A potential pathogenic significance of HERV-W Env/Syncytin-1 expression in MS can be inferred from data showing that Syncytin-1 has indirect cytotoxic effects on oligodendrocytes in vitro, and that expression of Syncytin-1 in murine models results in demyelination in vivo [18, 22].
A putative MSRV Env protein was previously reported [7, 23]. The surface (SU) domain of MSRV Env (AF331500), the amino acid sequence of which is 87% identical to Syncytin-1, has been shown to stimulate production of proinflammatory cytokines in human monocytes via engagement of CD14 and toll-like receptor 4. It also triggers maturation of human dendritic cells and confers on them the potential to polarize naive T-cells into Th-1-like effector T-lymphocytes . Production of IFN-γ, IL-6, and IL-12p40 following stimulation with MSRV Env SU of peripheral blood mononuclear cells (PBMC) from patients with MS is significantly higher compared to PBMC from healthy controls . In addition, MSRV Env induces a polyclonal activation of T-cells bearing specific Vβ chains, reminiscent of immunopathogenic effects triggered by superantigens . Altogether, those proinflammatory properties of MSRV Env appear compatible with a potential relevance of MSRV Env in the context of MS, too.
Despite the possible role of MSRV and HERV-W (Syncytin-1) in MS, the exact origin of MSRV and the precise relationship between MSRV and HERV-W have been hitherto unclear. MSRV has been defined by different overlapping cDNA clones that were generated from particle-associated RNA from plasma or supernatants of cultured cells from patients with MS [6, 7, 23]. Sequences of those cloned cDNAs were reported to be similar to HERV-W, but individual HERV-W proviruses from which those sequences may have originated could not be identified so far. To date, there exists no molecular clone containing a complete infectious MSRV genome, and the very nature of MSRV has thus remained uncertain [26, 27]. This has further stimulated a discussion over the relative contributions of MSRV and HERV-W to the pathogenesis of MS [20, 28–30].
By analogy to some well-characterized animal retroviruses (e.g. Jaagsiekte sheep retrovirus), it has been speculated that MSRV could be an exogenous member of an endogenous retrovirus family which consequently might be able to form novel proviral insertions in human genomic DNA [6, 30–32]. If this was the case, it may be assumed that RNA transcripts from such proviral MSRV insertions are specifically detectable in individuals infected with MSRV. However, complicating the identification of such MSRV RNA transcripts, many of the HERV-W elements present in human genomic DNA, which all are very similar to MSRV, may produce RNA transcripts as well. We therefore asked which HERV-W loci are transcriptionally active in human cells, and whether transcripts corresponding to previously published MSRV sequences might be specifically detectable in patients with MS. Given the potential pathogenic role of MSRV/HERV-W Env proteins in MS, we focused our analysis on MSRV/HERV-W env sequences. Published MSRV env sequences were obtained from plasma of patients with MS [7, 23]. Since MSRV sequences found in plasma are likely to originate from PBMC and since MSRV/HERV-W env sequences have previously been detected in human PBMC [18–21] we chose to analyze PBMC in this work.
Herein, we identify transcriptionally active HERV-W env loci in PBMC from patients with MS and healthy controls. We also demonstrate that analysis of transcribed HERV-W env elements is complicated by frequent recombinations that are most likely generated in vitro. Based on these results we show that the published MSRV env and gag sequences can be explained as originating from endogenous HERV-W loci or recombinations among them. By clarifying the origin of MSRV sequences, our data help to settle a longstanding debate, and have important implications for the status of MSRV as well as for the potential role of MSRV/HERV-W Env in the pathogenesis of MS.
Transcriptionally active HERV-W env loci in human PBMC
Characteristics of HERV-W env loci identified in this study as transcribed in human PBMC
HERV-W env locus
location of amplicon in genome
full-length Env ORF
+ (538 aa)
In accordance with previous analyses of the coding capacity of the HERV-W family [14, 15, 40], except for the 7q21.2 HERV-W env locus, none of the transcribed HERV-W env loci disclosed ORFs for full-length Env proteins. Still, a transcriptionally active HERV-W env locus on chromosome Xq22.3 contains an almost complete env ORF, only interrupted by a single premature stop codon in its 5' region (codon 39) followed by several in-frame ATGs. If the longest possible env ORF from this transcribed locus were translated, starting at an in-frame ATG at codon 68, the Xq22.3 HERV-W env locus could give rise to an N-terminally truncated 475 amino acid HERV-W Env protein.
A close inspection of HERV-W env cDNAs reveals a high number of recombined sequences
Ideally, a HERV-W env cDNA sequence is expected to display no nucleotide mismatches to the genomic HERV-W env locus that it originated from. About one third of HERV-W env cDNAs analyzed in this work indeed perfectly matched with genomic DNA sequences. However, the remaining two thirds of HERV-W env cDNAs contained between 1 and 24 nucleotide differences compared to the best matching genomic HERV-W env locus. Although minor nucleotide differences may well be explained by the inaccuracy of Taq polymerase, sequencing errors, or sequence variations (SNPs) in genomic HERV-W env loci, those possibilities seem unlikely to account for the relatively high numbers of nucleotide mismatches observed in some of the cDNA sequences. It has recently been shown that analyses of transcribed HERV sequences are complicated by recombinations between individual HERV transcripts, which most likely arise in vitro during reverse transcription because of template switches of reverse transcriptase and/or through PCR-mediated recombinations . To investigate whether similar recombinations also occurred in the present study, we generated multiple sequence alignments of the 7 transcribed HERV-W env loci and the 332 HERV-W env cDNA sequences. A close inspection of multiple alignments unambiguously demonstrated that a high number of HERV-W env cDNAs, that is, 99 out of 332 (29.8%), represented recombinations between transcripts from different HERV-W env loci. Notably, the alleged breakpoints of recombined sequences appeared to be randomly distributed. Typical examples of recombined sequences are shown in Figure 2.
Recombined HERV-W env cDNA sequences generated in this study
Number of involved loci
Number of cDNA sequences
Sum of involved HERV-W env transcripts
Σ = 332
Σ = 471
Differential transcriptional activity of HERV-W env loci in human PBMC
To compare the transcriptional activity of different HERV-W env loci the PCR efficiencies with which these loci are amplified should be similar. PCR-efficiency mostly relies on the binding of primers to their target sequences. As can be seen in Figure 2, three of the seven amplified loci (7q21.2, 17q12, Xq22.3) perfectly matched to the primers and are thus expected to be amplified with similar efficiencies. The binding regions of three further loci (5q11.2, 6q21, 15q21.3) contained a single nucleotide mismatch in the 5' end of either the forward or the reverse primer. Because the 15q21.3 and the 6q21 loci were the most frequently cloned loci in our study, the single nucleotide mismatches in these two loci appear unlikely to have had a significant negative impact on amplification. In the case of the 5q11.2 locus this possibility cannot be excluded, but seems unlikely given the results for the 15q21.3 and 6q21 loci. One locus (14q21.3) displayed one additional nucleotide and one nucleotide mismatch in the binding region of the forward primer, and it seems possible that those mismatches may have negatively affected its amplification.
Since the finding of recombined HERV-W env cDNA sequences implies that the different HERV-W env loci that contributed to the recombined sequences must have been transcribed, we also estimated relative cloning frequencies based on the 99 recombined sequences, counting individually all HERV-W env loci that were present in the recombined sequences (see also Table 2). The relative cloning frequencies obtained in this evaluation were overall comparable to those of the non-recombined sequences (Figure 4), suggesting that the likelihood of taking part in recombinations correlates with transcript abundance.
Finally, not all transcriptionally active loci were detected as cDNA in every individual. Regarding non-recombined and recombined sequences together, transcripts from the 5q11.2 locus were detected in one, transcripts from the 14q21.3 in three, and transcripts from the 7q21.2 locus in seven of the eight individuals studied. The remaining HERV-W env loci (6q21, 15q21.3, 17q12, Xq22.3) were found to be transcriptionally active in the PBMC of every individual studied.
Similar transcriptional activity of individual HERV-W env loci in PBMC from patients with MS and healthy controls
MSRV sequences can be explained as originating from distinct HERV-W loci and recombinations among them
Having identified transcriptionally active HERV-W env loci in human PBMC, we were interested to know how previously published MSRV env sequences are related to those transcribed HERV-W env loci. Given the high frequency of in vitro recombinations between HERV-W env transcripts, we also wondered whether recombinations may be detectable in MSRV env sequences. To this end, we retrieved published MSRV sequences comprising parts of or the entire env region (n = 5) [7, 23] from the NCBI database and analyzed them by BLAT searches (see Methods for details). Since all published MSRV env sequences are heterogeneous, it is a priori unlikely that all those sequences are derived from a single proviral insertion.
Origin of previously described MSRV sequences
HERV-W source locus/loci
Number of nucleotide mismatches
AF127227 (544 bp)
AF127228 (1932 bp)
AF127229 (2004 bp)
AF123882 (2477 bp)
AF331500 (1629 bp)
AF123881 (1511 bp)
We could also assign the formerly published MSRV gag sequence (AF123881) to a distinct HERV-W element on chromosome 3q26.32. This HERV-W locus has formerly been identified as transcriptionally active in human PBMC  and is identical to a HERV-W gag gene on chromosome 3 (AF156961), previously characterized by Voisset et al. . Although the 3q26.32 HERV-W gag gene is incomplete, it contains the largest HERV-W gag ORF in the human genome, with a putative coding capacity for a 45 kDa HERV-W Gag protein, consisting of a complete matrix domain and a C-terminally truncated capsid, but lacking nucleocapsid .
Notably, the average number of nucleotide mismatches between MSRV env and gag sequences and the respective best matching genomic HERV-W loci (1.97/kb; see Table 3) were in the same range as that observed in our study of transcribed HERV-W env loci in human PBMC (1.53/kb). In summary, our analyses suggest that previously published MSRV sequences originated from genomic HERV-W loci, or recombinations among them.
Relationship between Xq22.3 HERV-W env and MSRV env
In contrast to AF331500, the MSRV env clone AF127228, which likewise originates from the Xq22.3 HERV-W env locus, displays the stop codon at position 39. A DNA fragment comprising amino acids 68 to 446 of the HERV-W env ORF encoded by AF127228 has previously been used to generate the monoclonal anti HERV-W Env antibody 6A2B2 . As shown in Figure 6, except for two amino acid exchanges in its C-terminus, the AF127228 amino acids 68 to 446 sequence is identical to the amino acid sequence of Xq22.3 HERV-W Env, but displays 38 mismatches to the Syncytin-1 amino acid sequence. Nevertheless, the 6A2B2 antibody may cross react with Syncytin-1 [16, 43] and all previous neuropathological studies that reported a dysregulated expression of Syncytin-1 in MS lesions relied on the 6A2B2 antibody [18, 21, 22]. Still, assuming that HERV-W Xq22.3 env may have the potential to code for a HERV-W Env protein, our findings open the intriguing possibility that the protein detected by 6A2B2 in MS lesions could instead have originated from the HERV-W env locus on chromosome Xq22.3.
When studying HERV RNA expression in human diseases, it seems important to clearly dissect from which genomic HERV loci the detected HERV RNA transcripts originate . Consistent with previous findings suggesting that expression of HERV transcripts is a ubiquitious phenomenon occurring in every human tissue [34, 45, 46], we herein show that at least seven HERV-W env loci are transcribed in PBMC from patients with MS and healthy controls. Since the primers used in this investigation only amplify a limited number of genomic HERV-W env elements our study is not exhaustive, and it seems rather likely that more than seven HERV-W env elements are transcriptionally active in human PBMC. Additionally, HERV-W env loci that are transcribed at very low rates could be missed in the cloning procedure unless much higher numbers of clones are generated. Three of the transcribed HERV-W env elements (15q21.3, Xq22.3, 17q12) identified in this study were previously found to be expressed in human PBMC by a cloning and sequencing approach [38, 39]. Assignment of cDNAs to genomic HERV-W env loci in the former investigations was based on rather short sequences (30 bp, excluding primers), containing only few informative nucleotides, that is, nucleotides that are exclusively present in a single genomic HERV-W env locus and thereby allow unambiguous assignment of cDNAs. Usage of a ~600 bp sequence (excluding primers) in the present work resulted in a higher number of informative nucleotides and thus strengthened the accuracy of the assignment. Our finding of ERVWE1 transcripts in human PBMC is consistent with previous observations [19, 20] and corroborates that although ERVWE1 expression is most abundant in placenta this locus is transcribed in non-placental tissues as well .
Several studies have analyzed expression of HERV-W env RNA in PBMC or brain tissue from patients with MS [18–21]. Lack of systematic cloning, sequencing, and assignment of cDNA sequences to genomic HERV-W env loci have impaired the exact identification of transcriptionally active genomic HERV-W env loci responsible for the observed HERV-W env RNA expression in these investigations. Whereas in the present detailed analysis we could identify distinct transcriptionally active HERV-W env loci, we did not observe significant differences in the transcriptional activity of those loci in PBMC from patients with MS versus healthy controls. Although the number of individuals studied was rather small, these data argue against a dysregulated transcription pattern of HERV-W env in PBMC from patients with MS. In contrast, a consistent finding of former investigations was a significantly higher global HERV-W env RNA expression in brain tissue from patients with MS as compared to brain tissue from patients with other neurological diseases or normal brain tissue [18–21]. Using the methodological approach of the present work, it will therefore be interesting to identify the HERV-W env elements underlying upregulated HERV-W env RNA expression in MS brain tissue.
Antony and coworkers addressed this question by designing primers that specifically amplify HERV-W env 7q21.2 (ERVWE1) and the MSRV env clone AF123882 , which, as shown by our analyses, corresponds to a HERV-W env element on chromosome 15q21.3. These authors also employed a pair of degenerate (HERV-Wdeg) env primers that were based on the MSRV env clone AF331500, which, again as shown in this work, corresponds to a recombined cDNA originating from HERV-W env loci on Xq22.3 and 5p12. According to the Antony et al. study, elevated HERV-W env RNA expression in MS brain tissue originates mainly from the HERV-W env elements amplified by the HERV-Wdeg env primers and (somewhat less) from HERV-W env 7q21.2, while HERV-W env 15q21.3 expression was similar in patients and controls . A BLAT-PCR search showed that the HERV-Wdeg env primer pair potentially amplifies at least three genomic HERV-W env loci, among them HERV-W env Xq22.3. It is thus tempting to speculate that HERV-W env Xq22.3 may significantly contribute to increased HERV-W env RNA expression in MS brain tissue. Again, using the methods described herein, this issue could be resolved in a straightforward manner.
Remarkably, we observed a high number (29.8%) of recombined sequences among the analyzed HERV-W env cDNAs. As detailed in a previous study on transcribed HERV-K(HML-2) sequences , those recombinant cDNA sequences very likely resulted from in vitro recombinations that were due to template switches of reverse transcriptase during cDNA synthesis and/or PCR-mediated recombinations. Both of these mechanisms are well-recognized and have been proven experimentally to produce chimeric sequences [41, 47–53]. The percentage of recombined sequences detected in the present study was higher than that in the study on HERV-K(HML-2) in which ~5% of recombined clones were observed . This is most likely explained by the fact that in the HERV-K(HML-2) study only cDNA sequences with more than 17 nucleotide mismatches to the best matching locus were analyzed for recombinations, whereas in the present work all cDNA sequences were scrutinized for recombinations. Altogether, our data indicate that during experimental studies of repetitive elements by RT-PCR, in vitro recombinations are relatively common and almost inevitable complications.
An important result of this investigation is that previously published MSRV env and gag sequences appear to either be derived from transcripts of specific genomic HERV-W elements or to result from recombinations among such transcripts (Table 3). Given the high frequency of in vitro recombinations between transcripts from different HERV loci observed in this and the study by Flockerzi et al. , and given that MSRV clones were generated by methodologically similar approaches, it seems possible that the recombined MSRV env sequences (AF127229, AF331500) have resulted from in vitro recombinations as well.
An alternative explanation is that the recombined MSRV env sequences, and the recombined HERV-W env sequences isolated in this study, originated from novel, recombined, genomic HERV-W insertions. Hypothetically, such insertions could have formed in vivo after recombination of RNA transcripts from different HERV-W env loci through template switches during reverse transcription. Although we cannot formally exclude this possibility, a number of points argue against it. First, all known HERV-W elements are defective and replication-incompetent [8, 14]. Therefore, HERV-W is a priori rather unlikely to have the capacity to form new insertions in human DNA. Second, if there were novel recombined HERV-W loci in human DNA, one would expect to repeatedly observe defined recombined sequences originating from such insertions. However, this is neither the case with the 99 recombined HERV-W env cDNA sequences analyzed in this study nor with the published MSRV env clones. Third, given that about 30% of HERV-W env and 33% (2 of 6) of the investigated MSRV sequences represent recombinants, if all these recombinant MSRV/HERV-W env sequences were derived from novel proviral insertions, formation of such novel insertions would be an astonishingly frequent event. It seems very unlikely that as many recombined HERV-W loci should have been overlooked in previous genome sequencing projects.
Collectively, the most plausible and simplest explanation for the origin of MSRV env and gag sequences seems to be that those sequences originate from RNA transcripts from various endogenous HERV-W loci, or from in vitro recombinations among them. All of the HERV-W loci from which MSRV sequences are derived are defective and except for the 5p12 HERV-W env element, all of those loci resemble processed HERV-W pseudogenes. The human genome sequence was not yet available when MSRV was described, which hampered the identification of the precise origin of MSRV sequences at that time. It was, however, noted that those sequences cannot be attributed to a single replication-competent genome . Nevertheless, the nature of MSRV was subsequently controversial, and it has been speculated that MSRV could be an exogenous, replication-competent retrovirus [6, 30–32]. In contrast, our present data clearly suggest that the published MSRV env and gag RNA sequences are not derived from the genome of a currently replication-competent exogenous retrovirus. In the light of these results and previous observations of an increased prevalence of MSRV pol transcripts in plasma from patients with MS as compared to healthy controls [54, 55], it may similarly be interesting to analyze which HERV-W pol elements those MSRV pol transcripts could be derived from.
Although our findings argue against MSRV being an autonomous retroviral entity, they do by no means rule out that individual HERV-W env loci that correspond to MSRV sequences, or the Syncytin-1 (ERVWE1) gene, could be of relevance in MS. Indeed, we show that two MSRV env clones (AF331500, AF127228), which have been instrumental for the characterization of proinflammatory effects of MSRV Env  and the generation of a monoclonal anti-MSRV/HERV-W Env antibody (6A2B2) , are derived from a HERV-W env locus on chromosome Xq22.3. This locus is highly remarkable as it is interrupted by only a single premature stop at codon position 39 and otherwise harbors a long ORF for a N-terminally truncated 475 amino acid HERV-W Env protein (Figure 6). Bonnaud and colleagues described frameshift insertions/deletions (indels), that is, indels whose length is not a multiple of three, in 33 out of 36 analyzed genomic HERV-W env loci. Interestingly, among the three loci without frameshift indels were the ERVWE1 gene and the Xq22.3 HERV-W env element . We further note that the 475 amino acid Xq22.3 HERV-W env ORF is also present in the orthologous locus in chimpanzees (data not shown). These data may be taken as hints that selective pressure could act on the Xq22.3 HERV-W env locus, raising the possibility that Xq22.3 HERV-W env could exert a biological function. Our finding that the Xq22.3 HERV-W env locus is transcriptionally active in human cells indicates that it fulfills at least one essential prerequisite for a protein expression capacity in vivo.
Neuropathological studies revealed that the 6A2B2 anti-HERV-W Env antibody reacts with an antigen that is strongly expressed by glial cells in MS brain lesions, but not in normal control brain tissue [18, 21, 22]. Because Syncytin-1 has been thought to be the only HERV-W env locus capable of producing a HERV-W Env protein, and because 6A2B2 may crossreact with Syncytin-1 [16, 43], the antigen detected by 6A2B2 in MS brain lesions was considered to be Syncytin-1. However, our analyses show that the protein against which the 6A2B2 antibody was raised is practically identical to the Xq22.3 HERV-W Env protein (Figure 6) . We meanwhile cloned Xq22.3 HERV-W env into a eukaryotic expression vector. Transient transfection of HeLa cells with this clone showed that the Xq22.3 HERV-W env has retained a coding capacity and can produce a HERV-W Env protein in vitro which is detected by the 6A2B2 antibody in immunocytochemistry and immunoblots (C. Crusius, S. Wahl, K. Ruprecht, manuscript in preparation). These data suggest that the antigen recognized by 6A2B2 in MS lesions could likewise originate from the Xq22.3 HERV-W env locus, provided that this locus has a protein expression capacity in vivo. More elaborate studies will be required to clarify the exact nature of the HERV-W Env protein detected in MS lesions. Further characterization of the putative Xq22.3-encoded HERV-W Env protein, especially in comparison to Syncytin-1, will be necessary for such clarification.
In conclusion, we demonstrate that several HERV-W env loci are transcribed in human PBMC, and that analysis of such transcribed HERV-W env elements is complicated by frequent recombinations, which are most likely generated in vitro. Based on these findings, we show that previously reported MSRV env and gag sequences can be explained as originating from (in some instances recombined) transcripts of defective HERV-W elements, arguing against MSRV sequences being derived from an infectious exogenous retrovirus. Our results should help to settle the issue of the nature of MSRV and contribute to the clarification of the roles of MSRV versus HERV-W Env (Syncytin-1) in MS. Indeed, our findings raise the intriguing possibility that a protein encoded by a HERV-W env element on chromosome Xq22.3 could be expressed in MS brain lesions.
Patients with MS and healthy controls
Four patients with MS (3 female, 1 male) and 4 healthy controls (2 female, 2 male) were included in this study. The median age of patients was 34 (range 29–39) and of controls 34.5 (29–41) years. Clinical data of patients with MS were obtained by review of the medical records. All patients with MS had a diagnosis of definite MS according to Poser's criteria . Three patients had relapsing-remitting MS, and one patient had secondary progressive MS. The median expanded disability status scale score of patients with MS was 3.25 (1.5–6.5). One patient was treated with interferon-beta 1a, and two patients were treated with glatiramer acetate by the time of blood collection. None of the patients had been treated with glucocorticosteroids for at least 6 months before blood collection. Participants provided written informed consent, and the study was approved by the ethics committee of the faculty of medicine, Julius-Maximilians University, Würzburg. PBMC samples used in this work were collected and purified with Lymphoprep (Axis Shield, Oslo, Norway) gradient centrifugation as described before . Samples were stored at -80°C prior to the present analysis.
PBMC were thawed and cultured overnight in RPMI 1640 (BioWhittaker) supplemented with 10% FCS and penicillin (100 U/ml) and streptomycin (100 μg/ml) at 37°C in a humidified 5% (v/v) CO2 atmosphere. Total RNA was extracted from PBMC using the RNeasy Mini kit (Qiagen) and eluted in 60 μl of distilled water. RNA concentration and purity were assessed spectrophotometrically. Contaminating DNA was removed using the TURBO DNA-free Kit (Ambion Inc.) following the protocol for rigorous DNAse treatment. In brief, 2 units of TURBO DNase were added to a 50-μl reaction containing 10 μg RNA and incubated for 30 minutes at 37°C. Another 2 units of TURBO DNase were added and the incubation was continued for 30 minutes at 37°C. DNAse was removed using 10 μl of the provided DNAse inactivation reagent. Subsequently, 0.3–0.5 μg of DNase digested cellular RNA was reverse transcribed in a 20-μl reaction using Superscript II (Invitrogen) and 25 μM random hexamer primers (MWG-Biotech AG) according to the protocol of the manufacturer. Negative controls were generated in parallel for each sample by omission of Superscript II from the reaction. PCR primer sequences for amplification of HERV-W env were as follows: forward primer 5'-TTCACTGCCCACACCCAT-3'; reverse primer 5'-GAGGTACCACAGACAAAAAATATTCCT-3'. Conventional PCR was performed in a 50-μl reaction containing 1 μl of cDNA, 0.5 μM of each primer, 200 μM of each dNTP, reaction buffer (10 mM Tris-HCl, 50 mM KCl, 1.5 mM MgCl2), and 0.05 units/μl of Taq DNA Polymerase (D1806, Sigma). Cycling parameters were as follows: 3 minutes at 95°C; 40 cycles of 50 sec at 95°C, 50 sec at 58°C, and 1 minute at 72°C; and 10 minutes at 72°C.
Cloning of HERV-W env transcripts and assignment to proviral HERV-W loci
PCR products were excised from agarose gels, purified (NucleoSpin Extract II, Macherey-Nagel), and ligated into the pGEM-T vector (Promega). Plasmid DNA from randomly selected insert-containing clones was purified with the QIAprep Miniprep kit (Qiagen) and sequenced on an Applied Biosystems 3730x Capillary Sequencer using vector-specific primers (Institut für Immunologie und Genetik, Kaiserlautern, Germany). The quality of chromatograms was assessed by visual inspection. Poor-quality reads (< 0.1% of all sequences) were excluded from the analysis.
Assignment of cDNA sequences to corresponding HERV-W env loci is based on random and thus characteristic nucleotide differences between the various genomic HERV-W env loci. The proviral HERV-W env locus with no or very few nucleotide mismatches to a HERV-W env cDNA sequence can be assumed to represent the origin of this cDNA, if all other alternative loci displayed more nucleotide differences. A detailed discussion of the sequence assignment strategy has recently been provided .
To assign HERV-W env cDNA clones to specific HERV-W env loci in the human genome, HERV-W env cDNA sequences were first analyzed by BLAT searches (http://genome.ucsc.edu/cgi-bin/hgBlat; March 2006 human genome assembly). To further study recombinations between different HERV-W env loci in HERV-W env cDNA sequences, sequences of the seven transcribed HERV-W env loci were retrieved from the human genome sequence (March 2006 assembly) at the Human Genome Browser and multiply aligned with HERV-W env cDNA sequences using Muscle 3.6 . Candidate HERV-W env cDNA sequences were then inspected for recombination events.
Analysis of MSRV sequences
Previously published MSRV env and gag sequences were retrieved from GenBank and analyzed by BLAT searches to identify endogenous HERV-W loci with similarities to MSRV sequences. Alignments of the MSRV sequences with the best matching HERV-W locus were manually inspected for evidence of recombination events. In recombined sequence portions, nucleotide mismatches between MSRV sequences and the best matching HERV-W sequence usually clustered in defined subregions. Presumably recombined subregions were used as probe sequences for another BLAT search to detect their best matching HERV-W locus. Sequences of thus identified HERV-W loci were again retrieved from the Human Genome Browser and aligned with the corresponding MSRV sequences.
This study was supported by grants from HOMFOR to KR and JM. JM and NML are furthermore supported by grants from the Deutsche Forschungsgemeinschaft (DFG).
- Noseworthy JH: Progress in determining the causes and treatment of multiple sclerosis. Nature. 1999, 399: A40-A47.View ArticlePubMedGoogle Scholar
- Gessain A, Barin F, Vernant JC, Gout O, Maurs L, Calendar A, de The G: Antibodies to HTLV-1 in patients with tropical spastic paraparesis. The Lancet. 1985, 2: 407-410. 10.1016/S0140-6736(85)92734-5.View ArticleGoogle Scholar
- Perron H, Geny C, Laurent A, Mouriquand C, Pellat J, Perret J, Seigneurin JM: Leptomeningeal cell line from multiple sclerosis with reverse transcriptase activity and viral particles. Research in Virology. 1989, 140: 551-561. 10.1016/S0923-2516(89)80141-4.View ArticlePubMedGoogle Scholar
- Perron H, Lalande B, Gratacap B, Laurent A, Genoulaz O, Geny C, Mallaret M, Schuller E, Stoebner P, Seigneurin JM: Isolation of retrovirus from patients with multiple sclerosis. Lancet. 1991, 337: 862-863. 10.1016/0140-6736(91)92579-Q.View ArticlePubMedGoogle Scholar
- Haahr S, Sommerlund M, Moller-Larsen A, Nielson R, Hansen HJ: Just another dubious virus in cells from a patient with multiple sclerosis?. Lancet. 1991, 337: 863-864. 10.1016/0140-6736(91)92580-U.View ArticlePubMedGoogle Scholar
- Perron H, Garson JA, Bedin F, Beseme F, Paranhos-Baccala G, Komurian-Pradel F, Mallet F, Tuke PW, Voisset C, Blond JL, et al: Molecular identification of a novel retrovirus repeatedly isolated from patients with multiple sclerosis. Proc Natl Acad Sci USA. 1997, 94: 7583-7588. 10.1073/pnas.94.14.7583.PubMed CentralView ArticlePubMedGoogle Scholar
- Komurian-Pradel F, Paranhos-Baccala G, Bedin F, Ounanian-Paraz A, Sodoyer M, Ott C, Rajoharison A, Garcia E, Mallet F, Mandrand B, et al: Molecular cloning and characterization of MSRV-related sequences associated with retrovirus-like particles. Virology. 1999, 260: 1-9. 10.1006/viro.1999.9792.View ArticlePubMedGoogle Scholar
- Blond J-L, Beseme F, Duret L, Bouton O, Bedin F, Perron H, Mandrand B, Mallet F: Molecular characterization and placental expression of HERV-W, a new human endogenous retrovirus family. J Virol. 1999, 73: 1175-1185.PubMed CentralPubMedGoogle Scholar
- Bannert N, Kurth R: Retroelements and the human genome: new perspectives on an old relation. Proc Natl Acad Sci USA. 2004, 101 (Suppl 2): 14572-14579. 10.1073/pnas.0404838101.PubMed CentralView ArticlePubMedGoogle Scholar
- de Parseval N, Heidmann T: Human endogenous retroviruses: from infectious elements to human genes. Cytogenet Genome Res. 2005, 110: 318-332. 10.1159/000084964.View ArticlePubMedGoogle Scholar
- Pavlicek A, Paces J, Elleder D, Hejnar J: Processed pseudogenes of human endogenous retroviruses generated by LINEs: their integration, stability, and distribution. Genome Res. 2002, 12: 391-399.PubMed CentralView ArticlePubMedGoogle Scholar
- Belshaw R, Katzourakis A, Paces J, Burt A, Tristem M: High copy number in human endogenous retrovirus families is associated with copying mechanisms in addition to reinfection. Mol Biol Evol. 2005, 22 (4): 814-817. 10.1093/molbev/msi088.View ArticlePubMedGoogle Scholar
- Costas J: Characterization of the intragenomic spread of the human endogenous retrovirus family HERV-W. Mol Biol Evol. 2002, 19 (4): 526-533.View ArticlePubMedGoogle Scholar
- Voisset C, Bouton O, Bedin F, Duret L, Mandrand B, Mallet F, Paranhos-Baccala G: Chromosomal distribution and coding capacity of the human endogenous retrovirus HERV-W family. AIDS Res Hum Retroviruses. 2000, 16: 731-740. 10.1089/088922200308738.View ArticlePubMedGoogle Scholar
- de Parseval N, Lazar V, Casella JF, Benit L, Heidmann T: Survey of human genes of retroviral origin: identification and transcriptome of the genes with coding capacity for complete envelope proteins. J Virol. 2003, 77 (19): 10414-10422. 10.1128/JVI.77.19.10414-10422.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Blond JL, Lavillette D, Cheynet V, Bouton O, Oriol G, Chapel-Fernandes S, Mandrand B, Mallet F, Cosset FL: An envelope glycoprotein of the human endogenous retrovirus HERV-W is expressed in the human placenta and fuses cells expressing the type D mammalian retrovirus receptor. J Virol. 2000, 74: 3321-3329. 10.1128/JVI.74.7.3321-3329.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Mallet F, Bouton O, Prudhomme S, Cheynet V, Oriol G, Bonnaud B, Lucotte G, Duret L, Mandrand B: The endogenous retroviral locus ERVWE1 is a bona fide gene involved in hominoid placental physiology. Proc Natl Acad Sci USA. 2004, 101 (6): 1731-1736. 10.1073/pnas.0305763101.PubMed CentralView ArticlePubMedGoogle Scholar
- Antony JM, van Marle G, Opii W, Butterfield DA, Mallet F, Yong VW, Wallace JL, Deacon RM, Warren K, Power C: Human endogenous retrovirus glycoprotein-mediated induction of redox reactants causes oligodendrocyte death and demyelination. Nat Neurosci. 2004, 7 (10): 1088-1095. 10.1038/nn1319.View ArticlePubMedGoogle Scholar
- Antony JM, Izad M, Bar-Or A, Warren KG, Vodjgani M, Mallet F, Power C: Quantitative analysis of human endogenous retrovirus-W env in neuroinflammatory diseases. AIDS Res Hum Retroviruses. 2006, 22 (12): 1253-1259. 10.1089/aid.2006.22.1253.View ArticlePubMedGoogle Scholar
- Antony JM, Zhu Y, Izad M, Warren KG, Vodjgani M, Mallet F, Power C: Comparative expression of human endogenous retrovirus-W genes in multiple sclerosis. AIDS Res Hum Retroviruses. 2007, 23 (10): 1251-1256. 10.1089/aid.2006.0274.View ArticlePubMedGoogle Scholar
- 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-274. 10.1099/vir.0.81890-0.View ArticlePubMedGoogle Scholar
- Antony JM, Ellestad KK, Hammond R, Imaizumi K, Mallet F, Warren KG, Power C: The human endogenous retrovirus envelope glycoprotein, syncytin-1, regulates neuroinflammation and its receptor expression in multiple sclerosis: a role for endoplasmic reticulum chaperones in astrocytes. J Immunol. 2007, 179 (2): 1210-1224.View ArticlePubMedGoogle Scholar
- Perron H, Jouvin-Marche E, Michel M, Ounanian-Paraz A, Camelo S, Dumon A, Jolivet-Reynaud C, Marcel F, Souillet Y, Borel E, et al: Multiple sclerosis retrovirus particles and recombinant envelope trigger an abnormal immune response in vitro, by inducing polyclonal Vbeta16 T-lymphocyte activation. Virology. 2001, 287 (2): 321-332. 10.1006/viro.2001.1045.View ArticlePubMedGoogle Scholar
- Rolland A, Jouvin-Marche E, Viret C, Faure M, Perron H, Marche PN: The envelope protein of a human endogenous retrovirus-W family activates innate immunity through CD14/TLR4 and promotes Th1-like responses. J Immunol. 2006, 176 (12): 7636-7644.View ArticlePubMedGoogle Scholar
- Rolland A, Jouvin-Marche E, Saresella M, Ferrante P, Cavaretta R, Creange A, Marche P, Perron H: Correlation between disease severity and in vitro cytokine production mediated by MSRV (Multiple Sclerosis associated RetroViral element) envelope protein in patients with multiple sclerosis. J Neuroimmunol. 2005, 160 (1–2): 195-203. 10.1016/j.jneuroim.2004.10.019.View ArticlePubMedGoogle Scholar
- Blomberg J, Ushameckis D, Jern P: Evolutionary aspects of human endogenous retroviral sequences (HERVs) and disease. Retroviruses and primate genome evolution. Edited by: Sverdlov ED. 2005, Eurekah.com, 204-238.Google Scholar
- Voisset C, Weiss RA, Griffiths DJ: Human RNA "rumor" viruses: the search for novel human retroviruses in chronic disease. Microbiol Mol Biol Rev. 2008, 72 (1): 157-196. 10.1128/MMBR.00033-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Garson J, Creange A, Dolei A, Ferrante P, Jouvin-Marche E, Marche PN, Rieger F, Ruprecht K, Saresella M, Sotgiu S, et al: MSRV, Syncytin and the role of endogenous retroviral proteins in demyelination. Mult Scler. 2005, 11 (2): 249-250. 10.1191/1352458505ms1160xx.View ArticlePubMedGoogle Scholar
- Dolei A: MSRV/HERV-W/syncytin and its linkage to multiple sclerosis: the usability and the hazard of a human endogenous retrovirus. J Neurovirol. 2005, 11 (2): 232-235. 10.1080/13550280590952899.View ArticlePubMedGoogle Scholar
- Dolei A, Perron H: The multiple sclerosis-associated retrovirus and its HERV-W endogenous family: a biological interface between virology, genetics, and immunology in human physiology and disease. J Neurovirol. 2009, 15 (1): 4-13. 10.1080/13550280802448451.View ArticlePubMedGoogle Scholar
- Serra C, Mameli G, Arru G, Sotgiu S, Rosati G, Dolei A: In vitro modulation of the multiple sclerosis (MS)-associated retrovirus by cytokines: implications for MS pathogenesis. J Neurovirol. 2003, 9 (6): 637-643. 10.1080/714044485.PubMedGoogle Scholar
- Perron H, Perin J-P, Rieger F, Alliel PM: Particle-associated retroviral RNA and tandem RGH/HERV-W copies on human chromosome 7q: possible components of a 'chain-reaction' triggered by infectious agents in multiple sclerosis?. J Neurovirol. 2000, 6: S67-S75.PubMedGoogle Scholar
- Mayer J, Ehlhardt S, Seifert M, Sauter M, Muller-Lantzsch N, Mehraein Y, Zang KD, Meese E: Human endogenous retrovirus HERV-K(HML-2) proviruses with Rec protein coding capacity and transcriptional activity. Virology. 2004, 322 (1): 190-198. 10.1016/j.virol.2004.01.023.View ArticlePubMedGoogle Scholar
- Flockerzi A, Ruggieri A, Frank O, Sauter M, Maldener E, Kopper B, Wullich B, Seifarth W, Muller-Lantzsch N, Leib-Mosch C, et al: Expression patterns of transcribed human endogenous retrovirus HERV-K(HML-2) loci in human tissues and the need for a HERV Transcriptome Project. BMC Genomics. 2008, 9: 354-10.1186/1471-2164-9-354.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruprecht K, Ferreira H, Flockerzi A, Wahl S, Sauter M, Mayer J, Mueller-Lantzsch N: Human endogenous retrovirus family HERV-K(HML-2) RNA transcripts are selectively packaged into retroviral particles produced by the human germ cell tumor line Tera-1 and originate mainly from a provirus on chromosome 22q11.21. J Virol. 2008, 82 (20): 10008-10016. 10.1128/JVI.01016-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ: BLAT – the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-664.PubMed CentralView ArticlePubMedGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110 (1–4): 462-467. 10.1159/000084979.View ArticlePubMedGoogle Scholar
- Nellaker C, Yao Y, Jones-Brando L, Mallet F, Yolken RH, Karlsson H: Transactivation of elements in the human endogenous retrovirus W family by viral infection. Retrovirology. 2006, 3: 44-10.1186/1742-4690-3-44.PubMed CentralView ArticlePubMedGoogle Scholar
- Yao Y, Schroder J, Nellaker C, Bottmer C, Bachmann S, Yolken RH, Karlsson H: Elevated levels of human endogenous retrovirus-W transcripts in blood cells from patients with first episode schizophrenia. Genes Brain Behav. 2007, 7 (1): 103-112.PubMedGoogle Scholar
- Villesen P, Aagaard L, Wiuf C, Pedersen FS: Identification of endogenous retroviral reading frames in the human genome. Retrovirology. 2004, 1: 32-10.1186/1742-4690-1-32.PubMed CentralView ArticlePubMedGoogle Scholar
- Flockerzi A, Maydt J, Frank O, Ruggieri A, Maldener E, Seifarth W, Medstrand P, Lengauer T, Meyerhans A, Leib-Mosch C, et al: Expression pattern analysis of transcribed HERV sequences is complicated by ex vivo recombination. Retrovirology. 2007, 4: 39-10.1186/1742-4690-4-39.PubMed CentralView ArticlePubMedGoogle Scholar
- Pinter A, Kopelman R, Li Z, Kayman SC, Sanders DA: Localization of the labile disulfide bond between SU and TM of the murine leukemia virus envelope protein complex to a highly conserved CWLC motif in SU that resembles the active-site sequence of thiol-disulfide exchange enzymes. J Virol. 1997, 71 (10): 8073-8077.PubMed CentralPubMedGoogle Scholar
- Bonnaud B, Bouton O, Oriol G, Cheynet V, Duret L, Mallet F: Evidence of selection on the domesticated ERVWE1 env retroviral element involved in placentation. Mol Biol Evol. 2004, 21 (10): 1895-1901. 10.1093/molbev/msh206.View ArticlePubMedGoogle Scholar
- Ruprecht K, Mayer J, Sauter M, Roemer K, Mueller-Lantzsch N: Endogenous retroviruses and cancer. Cell Mol Life Sci. 2008, 65: 3366-3382. 10.1007/s00018-008-8496-1.View ArticlePubMedGoogle Scholar
- Seifarth W, Frank O, Zeilfelder U, Spiess B, Greenwood AD, Hehlmann R, Leib-Mosch C: Comprehensive analysis of human endogenous retrovirus transcriptional activity in human tissues with a retrovirus-specific microarray. J Virol. 2005, 79 (1): 341-352. 10.1128/JVI.79.1.341-352.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Stauffer Y, Theiler G, Sperisen P, Lebedev Y, Jongeneel CV: Digital expression profiles of human endogenous retroviral families in normal and cancerous tissues. Cancer Immun. 2004, 4: 2-PubMedGoogle Scholar
- Fagan JB, Pastan I, de Crombrugghe B: Sequence rearrangement and duplication of double stranded fibronectin cDNA probably occurring during cDNA synthesis by AMV reverse transcriptase and Escherichia coli DNA polymerase I. Nucleic Acids Res. 1980, 8 (13): 3055-3064. 10.1093/nar/8.13.3055.PubMed CentralView ArticlePubMedGoogle Scholar
- Bowman RR, Hu WS, Pathak VK: Relative rates of retroviral reverse transcriptase template switching during RNA- and DNA-dependent DNA synthesis. J Virol. 1998, 72 (6): 5198-5206.PubMed CentralPubMedGoogle Scholar
- Luo GX, Taylor J: Template switching by reverse transcriptase during DNA synthesis. J Virol. 1990, 64 (9): 4321-4328.PubMed CentralPubMedGoogle Scholar
- Ouhammouch M, Brody EN: Temperature-dependent template switching during in vitro cDNA synthesis by the AMV-reverse transcriptase. Nucleic Acids Res. 1992, 20 (20): 5443-5450. 10.1093/nar/20.20.5443.PubMed CentralView ArticlePubMedGoogle Scholar
- Svarovskaia ES, Delviks KA, Hwang CK, Pathak VK: Structural determinants of murine leukemia virus reverse transcriptase that affect the frequency of template switching. J Virol. 2000, 74 (15): 7171-7178. 10.1128/JVI.74.15.7171-7178.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Judo MS, Wedel AB, Wilson C: Stimulation and suppression of PCR-mediated recombination. Nucleic Acids Res. 1998, 26 (7): 1819-1825. 10.1093/nar/26.7.1819.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyerhans A, Vartanian JP, Wain-Hobson S: DNA recombination during PCR. Nucleic Acids Res. 1990, 18 (7): 1687-1691. 10.1093/nar/18.7.1687.PubMed CentralView ArticlePubMedGoogle Scholar
- Garson JA, Tuke PW, Giraud G, Paranhos-Baccala G, Perron H: Detection of virion-associated MSRV-RNA in serum of patients with multiple sclerosis. Lancet. 1998, 351: 33-10.1016/S0140-6736(98)24001-3.View ArticlePubMedGoogle Scholar
- Dolei A, Serra C, Mameli G, Pugliatti M, Sechi G, Cirotto MC, Rosati G, Sotgiu S: Multiple sclerosis-associated retrovirus (MSRV) in Sardinian MS patients. Neurology. 2002, 58 (3): 471-473.View ArticlePubMedGoogle Scholar
- Poser CM, Paty DW, Scheinberg L: New diagnostic criteria for multiple sclerosis: guidelines for research protocols. Ann Neurol. 1983, 13: 227-231. 10.1002/ana.410130302.View ArticlePubMedGoogle Scholar
- Gronen F, Ruprecht K, Weissbrich B, Klinker E, Kroner A, Hofstetter H, Rieckmann P: Frequency analysis of HLA-B7-restricted Epstein-Barr virus-specific cytotoxic T lymphocytes in patients with multiple sclerosis and healthy controls. J Neuroimmunol. 2006, 180: 185-192. 10.1016/j.jneuroim.2006.08.008.View ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.