Upregulation of selected HERVW loci in multiple sclerosis
Mobile DNA volume 12, Article number: 18 (2021)
Abstract and Introduction
Human endogenous retrovirus (HERV) are the present day versions of retroviral germline infections that have occured millions of years ago, which occupy about 8 % of the genome . While they are mostly replication deficient, they are known to express RNA and protein  during particular developmental stages, or as a response to aging , inflammation and a wide range of pathologies . A human retrovirus discovered in Multiple Sclerosis (MS) patients , turned out to be the prototype of a novel HERV family referred to as HERVW . The HERVW family consists of 213 elements, 12 out of which are complete proviral copies with intact LTRs . Increased expression of HERVW in peripheral blood mononuclear cells (PBMCs) has been repeatedly associated with MS, and the presence of HERVW protein or elevated RNA transcription has been correlated with disease activity [8,9,10]. While a contribution of HERVW-encoded proteins to brain disease is suggested by their presence in MS-associated brain lesions, expression in peripheral organs may be involved in the disease process through cytokine-induced damage to the blood brain barrier and subsequent infiltration of monocytes. Alterations in peripheral expression may also serve as a useful and practical marker for the diagnostics of this CNS disease. Therefore, we quantified overall HERVW levels and identified individual HERVW loci actually transcribed in PBMCs. Analysis was carried out in patients diagnosed with Clinically Isolated Syndrome (CIS), a precursor to MS, defined by a single episode of neurologic symptoms lasting at least 24 h. CIS is an indicator of future development of MS, as 60 % of the people diagnosed with CIS develop MS . These patients potentially represent the earliest stage of MS routinely available for clinical analysis. We undertook a Next Generation Sequencing (NGS)-based analysis of transcripts amplified from cDNA obtained from patients with CIS and samples from healthy controls. Data presented from this pilot experiment indicate that the relative frequency of specific HERVW copies is altered in PBMC of CIS patients, even in the absence of overall HERVW overexpression. Such altered frequency appears to be derived from less abundantly transcribed but potentially MS-related HERVW loci.
Patients (Table 1)
The local ethics committee (CEICA) approved the study protocol (CP - CI PI14/0021 dated 26/02/2014; modified on 25/10/2017), and patients provided written informed consent (protocolo y información para el paciente v2 de 29/12/2013). Blood samples were collected from MS patients and healthy controls from the Neurology Department of Miguel Servet University Hospital (Zaragoza, Spain). Whole fresh blood was drawn into vacutainer tubes (Becton Dickinson Vacutainer) containing EDTA. Within 24 h, PBMC were isolated as previously described [12, 13].
Expression analysis and PCR
RNA isolation and random-primed cDNA synthesis  was carried out as described before. HERVW ENV levels were determined by triplicate qPCR assays as described [14, 15]. For the identification and localization of transcribed HERVW loci, cDNA was amplified employing the external primers of an established PCR assay for HERVW ENV . Products were purified and subjected to NGS analysis.
Library preparation and sequencing was carried out using the IonTorrent technology workflow on an Ion Torrent S5XL platform using an Ion 530 chip. Resulting reads were mapped to the human reference genome (version hg19) using strict criteria to maximize mapping differences between different HERVW copies. Relative frequencies were calculated as the number of reads mapping to an individual HERVW ENV element relative to the total number of reads. Details in Suppl. Methods.
SPSS software was used for all analyses and graphs (Version 15.0). Normality and statistical significance of differences were assessed using specific tests. Data were further analyzed using the DESeq2 package  to correct p values for multiple testing (False Discovery Rate < 0.05).
More detailed information is available in Suppl. M&M.
We carried out HERVW ENV expression analyses using an optimized assay described by Mameli et al. . No significantly increased expression of HERVW was detected in a small cohort of CIS patients (n = 6) compared to age-matched controls (n = 15)(U-Mann-Whitney p = 0.267) (Fig. 1). Results were not skewed by the use of GAPDH as a reference gene (Fig. 1), as comparison with RPL19 and HSDA reference genes (Table S1 and Suppl Figure 1) showed that there is no statistical difference between the use of either GAPDH or the mean of the three genes (Welch´s t-test; p < 0.05).
In the absence of increased overall expression levels of HERVW in CIS samples, we wondered whether specific copies of HERVW (Table S2) might be differentially expressed. We performed NGS analysis to identify individual HERVW copies with altered expression in PBMC from CIS patients (n = 5) and controls (n = 5). Reads obtained (70,694 ± 24,812 per sample; 25,286–136,704) were mapped to the human genome. Once assigned to unique genomic locations, reads corresponding to 39 HERVW ENV loci were extracted (Table S3 and Table 2). As expected, > 99.85 % of mapped reads correspond to the 39 loci analyzed (data not shown). The resulting data showed that reads obtained from CIS patients mapped to a significant higher number of different HERVW ENV loci (31 ± 13), compared to those obtained from controls (16 ± 5.5) (t-student; p = 0.018) (Fig. 2 A). Over 70 % of the reads mapped to either of two loci: 19q13.2, Xq22.3. Extending the range, reads mapped with high frequency (> 3.6 % of total reads/locus) to a limited number of loci, in particular to HERVW ENV copies located on chromosomes 19q13.2, Xq22.3, 8q21.11, 15q21.3, 12q23.3 and 4q21.22 (Fig. 2 B). We found no significant differences between CIS patients and controls in the relative frequency of reads mapping to these loci (Fig. 2 B; Table 2).
Lower numbers of reads mapped to the remaining 33 loci, with relative frequencies ranging from 0,01 to 1,66 %. We found differences between CIS patients and controls in the relative frequency of reads mapping to several of these. Reads mapping to a subset of HERVW loci, including copy_chr3-1, copy_chr4-3, copy_chr5-1, copy_chr9-1, copy_chr9-2, copy_chr12-4 and copy_chr19-1 differed between CIS patients and controls (Table 3; Fig. 2 C-E). The relative frequency (ranging from 0,3 − 2,95 %) of reads mapping to these loci showed a > 7 fold increase in CIS patients compared to controls (Table 3). When corrected for multi-locus analysis, significant differential expression (p < 0.05; FDR < 0.05) of five of these loci loci was evident (Table 3), with significant increases in expression of HERVW copies 3q11.21, 4q31.1, 9q31.3 and 19p12 and a significant decrease in 12q23.3.
In contrast to the small group of CIS patients analyzed in this study, increased HERVW levels have been associated frequently with MS. Our inability to demonstrate a statistically significant increase of overall HERVW levels in PBMC of CIS patients may be explained by the selection of this particular group or more likely simply by small sample size. However, lack of increased expression is not unprecedented as it was previously reported in a cohort of South African MS patients, although different primers were used for this analysis .
We perfomed NGS analysis to identify individual HERVW copies that show altered expression in PBMC, comparing CIS patients (n = 5) to controls (n = 5). Although more definite answers require future analysis of more subjects, in the CIS patients analyzed more HERVW loci are expressed than in control subjects. A similar increase has been reported previously in MS brain . While previous studies failed to identify MS-specific loci or expression [18, 19], in the CIS patients we found statistically significant overrepresentation of reads corresponding to specific loci (i.e. 3q11.2 and 19p12, see Table 3 for complete list). Locus-specific qPCR assays may first help confirm this finding in a larger patient cohort, and subsequently be evaluated as a potential prognostic assay.
These combined overrepresented loci produce only 1–3 % of total transcripts (Fig. 2 C-E). The combined findings on low levels of overexpression, activation of more loci, and activation of low-expressing HERVW elements in CIS patients suggest that their potential contribution to the pathology may be unrelated to overall high expression levels. None of the copies identified encode full-length ENV protein, as the sequences corresponding to the ENV gene are truncated, lack ATG codons, and/or carry frame shifts and STOP codons (Suppl Figure 2). CIS-associated copies may produce proteins (either or not ENV-related) that are especially active in activation of TLR4 , or RNAs that trigger the native immune system through TLR3 [21, 22]. Although our analysis shows that upregulation of specific HERVW loci in PBMC is associated with CIS, the presence of these transcripts in MS brain is unknown at present. A potential role of these transcripts in proviral protein production and activation of either the peripheral immune system or CNS disease remains to be established.
Availability of data and materials
The datasets used and/or analyzed in this study have been deposited in NCBI's Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE173929 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE173929)
Grandi N, Tramontano E. Human endogenous retroviruses are ancient acquired elements still shaping innate immune responses. Front Immunol. 2018;9:2039.
Bourque G, Burns KH, Gehring M, Gorbunova V, Seluanov A, Hammell M, et al. Ten things you should know about transposable elements. Genome Biol. 2018;19:1–12.
Sedivy JM, Kreiling JA, Neretti N, Cecco M, De, Criscione SW, Hofmann JW, et al. Death by transposition - the enemy within? BioEssays. 2013;35:1035–43.
Cardelli M, Marchegiani F. Good, bad, mobile elements: genome’s most successful “parasites” as emerging players in cell and organismal aging. Curr Pharm Des. 2013;19:1739–52.
Perron H, Garson JA, Bedin F, Beseme F, Paranhos-Baccala G, Komurian-Pradel F, et al. Molecular identification of a novel retrovirus repeatedly isolated from patients with multiple sclerosis. The Collaborative Research Group on Multiple Sclerosis. Proc Natl Acad Sci U S A. 1997;94:7583–8.
Küry P, Nath A, Créange A, Dolei A, Marche P, Gold J, et al. Human Endogenous Retroviruses in Neurological Diseases. Trends Mol. Med. 2018. p. 379–94.
Grandi N, Cadeddu M, Blomberg J, Tramontano E. Contribution of type W human endogenous retroviruses to the human genome: characterization of HERV-W proviral insertions and processed pseudogenes. Retrovirology. 2016;13:67.
Perron H, Germi R, Bernard C, Garcia-Montojo M, Deluen C, Farinelli L, et al. Human endogenous retrovirus type W envelope expression in blood and brain cells provides new insights into multiple sclerosis disease. Mult Scler J. 2012;18:1721–36.
García-Montojo M, de la Hera B, Varadé J, de la Encarnación A, Camacho I, Domínguez-Mozo M, et al. HERV-W polymorphism in chromosome X is associated with multiple sclerosis risk and with differential expression of MSRV. Retrovirology. 2014;11:2.
Morandi E, Tanasescu R, Tarlinton RE, Constantinescu CS, Zhang W, Tench C, et al. The association between human endogenous retroviruses and multiple sclerosis: A systematic review and meta-analysis. PLoS One. 2017;12.
Kuhle J, Disanto G, Dobson R, Adiutori R, Bianchi L, Topping J, et al. Conversion from clinically isolated syndrome to multiple sclerosis: A large multicentre study. Mult Scler. 2015;21:1013–24.
Böyum A. Isolation of mononuclear cells and granulocytes from human blood. Isolation of monuclear cells by one centrifugation, and of granulocytes by combining centrifugation and sedimentation at 1 g. Scand J Clin Lab Invest Suppl. 1968;97:77–89.
Macías-Redondo S, Strunk M, Guillén M, Ara J, Martín J, Alvaréz-Lafuente R, et al. Searching for non-referenced HERVW copies relevant to MS: a splinkerette-PCR method. Submitted. 2020.
Pérez-Palacios R, Macías-Redondo S, Climent M, Contreras-Moreira B, Muniesa P, Schoorlemmer J. In Vivo Chromatin Targets of the Transcription Factor Yin Yang 2 in Trophoblast Stem Cells. PLoS One. 2016;11:e0154268.
Mameli G, Poddighe L, Astone V, Delogu G, Arru G, Sotgiu S, et al. Novel reliable real-time PCR for differential detection of MSRVenv and syncytin-1 in RNA and DNA from patients with multiple sclerosis. J Virol Methods. 2009;161:98–106.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
De Villiers JNP, Treurnicht FK, Warnich L, Carr J, Van Rensburg SJ, Kotze MJ. Analysis of viral and genetic factors in South African patients with multiple sclerosis. Metab Brain Dis. 2006;21:163–9.
Schmitt K, Richter C, Backes C, Meese E, Ruprecht K, Mayer J. Comprehensive Analysis of Human Endogenous Retrovirus Group HERV-W Locus Transcription in Multiple Sclerosis Brain Lesions by High-Throughput Amplicon Sequencing. J Virol. 2013;87:13837–52.
Laufer G, Mayer J, Mueller BF, Mueller-Lantzsch N, Ruprecht K. Analysis of transcribed human endogenous retrovirus W env loci clarifies the origin of multiple sclerosis-associated retrovirus env sequences. Retrovirology. 2009;6:37.
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:7636–44.
Hurst TP, Magiorkinis G. Activation of the innate immune response by endogenous retroviruses. J Gen Virol. 2015;96:1207–18.
Stetson DB. Endogenous retroelements and autoimmune disease. Curr Opin Immunol. 2012;24:692–7.
We acknowledge the use of the Sequencing and Functional Genomics Core Facility (Servicio Científico Técnico de Secuenciación y Genómica Funcional) of the Instituto Aragonés de Ciencias de la Salud Aragon Health Sciences Institute (IACS) and Universidad de Zaragoza. We are particularly grateful to all of the patients for their participation in and commitment to the study. Special thanks go to Carmen Ocabo for her help with informing patients, the collection and the control of clinical samples.
The work described in this manuscript has been financed by the Instituto Carlos III (Government of Spain)(FIS PI13/02518), and by several grants from the Government of Aragon, cofinanced by FEDER (“Una manera de hacer Europa, Construyendo Europa desde Aragón”/European Social Funds): B77 and B05_17D. SMR has been supported by PhD fellowship C071/2014 from the DGA (Aragón, Spain). The funding bodies played no role in the design of the study, in the collection, analysis, and interpretation of data, in writing the manuscript and in the decision to submit the article for publication.
Ethics approval and consent to participate
The local ethics committee (CEICA) approved the study protocol (CP - CI PI14/0021 dated 26/02/2014; modified on 25/10/2017), and patients provided written informed consent (protocolo y información para el paciente v2 de 29/12/2013).
Consent for publication
The authors report no conflict of interest
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this article was revised as the authors noticed the published Table 1 has a distorted headings and columns not in proper order.
The levels of GAPDH, RPL19 and HSDA were determined by qPCR in the samples indicated, as described in the legend to Suppl Fig. 1. The Table lists the Cts obtained for each gene (columns named accordingly, the mean of all three (column “mean”) and the difference between the Cts obtained using either GAPDH (G) or the mean of three reference genes (P) (Column P-G). The difference (mean 2.76; standard deviation 0.29) is statistically constant among samples (Welch´s t-test; p < 0.05).
List of HERVW ENV loci according to the GRCh37.p5 version of the human genome database. Loci identified by unbiased read mapping were verified as HERVW loci by comparison with the Repbase Update library of repeats from the Genetic Information Research Institute (GIRI), using the the RepeatMasker program. Nomenclature used in Tables is indicated, the genomic location of each copy as well as alternative names used in the literature. Sequences corresponding to the coordinates listed were downloaded and compared to primer sequences using the “align” function in SerialCloner (version 2-6-1). Identity with primer sequence is indicated in black, mismatches in red.
Number of reads mapping to the HERVW copies indicated in each of the five CIS or NDI (control) samples. The total number of reads per sample is indicated. The sum of reads in the CIS and NDI groups is indicated.
Comparison of reference genes.
Comparison of HERVW elements.
qPCR primer efficiency standard curve analysis.
About this article
Cite this article
Macías-Redondo, S., Strunk, M., Cebollada-Solanas, A. et al. Upregulation of selected HERVW loci in multiple sclerosis. Mobile DNA 12, 18 (2021). https://doi.org/10.1186/s13100-021-00243-1