Abstract

Background: It has been suggested that Human endogenous retroviruses (HERVs) are associated with multiple sclerosis (MS) pathogenesis. The objective of this study was to broadly evaluate the expression of HERV core (GAG) and envelope (ENV) genes in diseased brain white matter samples from MS patients compared to normal controls.

Methods: Twenty-eight HERV GAG and 88 ENV gene sequences were retrieved, classified by phylogeny, and grouped into clades. Consensus qPCR primers were designed for each clade, and quantitative PCR was performed on 33 MS and 9 normal control frozen brain samples. MS samples included chronic progressive (n=5), primary progressive (n=4), secondary progressive (n=14), relapsing remitting (n=3) and unclassified confirmed MS cases (n=7). The levels of GAG and ENV RNA within each of the samples were quantitated and normalized using the neuronal reference gene RPL19. Expression differences were analyzed for MS vs control.

Results: Expression of GAG clades 1A, 3B, and 3C mapping to HERV-E and HERV-K were significantly increased compared to controls, while GAG clade 3A expression was decreased. Expression of HERV ENV clades 2, 3A, 3B, mapping to RTVL, HERV-E and HERV-K and MSRV (HERV-W), were significantly increased in the MS group. However, the relative expression differences between the MS and control groups were small, differing less than 1.5-fold.

Conclusion: Expression of GAG and ENV mapping to HERV-E, RTVL and HERV-K10 families were significantly increased in the MS group. However, the relative expression differences between the MS and control groups were small, differing less than 1.5-fold. These results indicate that the expression of HERV GAG and ENV regions do not differ greatly between MS and controls in these frozen brain samples.

Keywords

Human endogenous retroviruses; HERV; qPCR; Multiple sclerosis; Primary progressive multiple sclerosis; GAG; Envelope; Gene expression

Introduction

Multiple sclerosis (MS) is an autoimmune disease characterized by multiple lesions with plaques in the brain and spinal cord. The disease is associated with an inflammatory process that attacks and destroys the myelin sheaths around axons in the brain and spinal cord. The cause of MS is unknown, although the disease is thought to be triggered by environmental factors operating on a predisposing genetic background. Genome-wide association studies (GWAS) have identified numerous susceptibility loci, mainly associated with immunological processes [1]. The human leukocyte antigen (HLA) immune gene-cluster, particularly HLA-DRB1*1501, is the main genetic risk factor for the disease [2,3].

There has been much speculation that an infectious agent may contribute to MS. Based on the epidemiology of MS, including geographic patterns, isolated outbreaks and migration studies, viruses have long been suspected as causative agents [4-6]. There is no direct evidence for the involvement of an exogenous retrovirus in MS, but endogenous retroviruses may be an important factor in MS pathogenesis. The possibility that retroviruses play a role in the pathogenesis of MS has been considered since the 1980s, when Perron showed reverse transcriptase (RT) activity and retroviral particles in leptomeningeal cells taken from a patient with MS [7-10]. Since that time, HERVs have been detected in blood, CSF, and brain tissue from MS patients, as reviewed by Douville RN, Nath A [11]. Human endogenous retroviruses (HERVs) are parts of the human genome, with approximately 98,000 ERV elements and fragments [12-14]. It is believed that HERVs are remnants of prehistoric exogenous retroviruses, integrated in the human genome during evolution through germ line infection [15]. Several transcripts and proteins of genes from the HERV-W family have been detected in the central nervous system and they have been associated with neuroinflammation [16- 18]. Antibodies reactive with GAG or ENV antigens from HERVs can be detected in human sera with elevated levels in patients with autoimmune disease [19].

Previously we performed RNA-seq on frozen human brain specimens and showed statistically significant expression differences between progressive MS and normal controls for some HERV domains [20,21]. Importantly, detection of this signal required using a bioinformatics approach called “comprehensive mapping” which counts all detectable alignments to a given sequencing read, rather than only the most specific match. Furthermore, aggregating sequencing hits by domains suggested that the expression of GAG (core) and ENV (envelope) domains is associated with the MS disease process. The objective of this study was to quantify the expression of selected HERV GAG and ENV domains in additional frozen, cryo preserved brain samples from MS patients in comparison to normal controls. Here, several qPCR experiments were performed to gain an understanding of the overall expression of HERV GAG and ENV genes, divided into multiple HERV clades.

Materials and Methods

a. Brain specimens

Cryo preserved white matter from MS plaques (N=33) and normal white matter control brain specimens from patients without any brain diseases (N=9) were obtained from the Rocky Mountain MS Center Tissue Bank (Westminster, CO, USA) and the UCLA Human Brain and Spinal Fluid Resource Center (Los Angeles, CA, USA). The MS samples included; chronic progressive MS (CPMS, n=5), primary progressive MS (PPMS, n=4), secondary progressive (SPMS, n=14), relapsing remitting MS (RRMS, n=3). Confirmed MS cases with no clinical subtype (unclassified MS, n=7) were also included in the study for a total of 33 MS samples. These specimens, as well as accompanying clinical and pathologic information (Table 1), were de-identified before shipment. The research plan was submitted to the University of Utah Health Sciences Center IRB for review. Since this research involved only de-identified post-mortem material, it was found to be exempt from IRB review and oversight (IRB #00028658).

Group	Specimen	PMI	Age	Sex	Collection Year	Clinical History/MS Type	Specimen Location
Normal Controls	5	2	55	F	1989	NA	cerebrum
	202*	4	57	M	1983	post-surgical infection	cerebrum
	214*	4	56	M	1981	heart disease	cerebrum
	3236	16	67	M	2002	esophageal cancer	frontal WM
	3276*	19	54	M	2002	heart disease	frontal WM
	33*	5	69	M	1983	lung disease	cerebrum
	3371*	16	52	M	2002	lung cancer	frontal WM
	3348*	9	76	F	2002	heart disease, diabetes	frontal WM
	3540	11	68	M	2003	lung cancer	frontal WM
MS	159	NA	53	M	1983	chronic progressive	WM
	164	NA	NA	M	1983	Unclassified	WM
	216	5	71	F	2003	Unclassified	WM
	2096	NA	NA	M	NA	secondary progressive	periventricular WM
	2485*	9	69	M	1997	primary progressive	frontal WM
	2696*	21	86	F	1998	primary progressive	periventricular WM
	2742	17	60	M	1998	primary progressive	periventricular WM, cerebellum
	2743	23	59	F	1998	secondary progressive	subcortical WM
	2758	NA	60	M	1999	secondary progressive	periventricular WM
	2778	39	61	F	1998	secondary progressive	WM, temporal cortex
	2800	9	64	F	1998	secondary progressive	WM, frontal cortex
	2946*	NA	59	M	1999	secondary progressive	periventricular WM, parietal lobe
	2966	NA	55	F	2000	secondary progressive	periventricular WM
	3010	NA	49	F	2000	secondary progressive	frontal WM, frontal lobe
	3056	10	61	M	2000	relapsing remitting	WM, frontal cortex
	3093	18	68	M	2000	Unclassified	periventricular WM, basal ganglia
	3161*	20	51	F	2001	secondary progressive	WM, frontal-parietal
	3250	17	75	M	2002	secondary progressive	periventricular WM
	3289	29	54	F	2001	Unclassified	periventricular WM
	3413	10	38	F	2002	relapsing remitting	periventricular WM
	3422	12	62	M	2002	relapsing remitting	periventricular WM, cerebellum
	3474	9	70	F	2002	chronic progressive	periventricular WM
	3502	16	78	M	2003	secondary progressive	periventricular WM,
	3509*	11	74	F	2003	primary progressive	periventricular WM
	3831	21	82	F	2004	chronic progressive	occipital WM
	3835	29	65	M	2004	chronic progressive	periventricular WM
	3867	13	75	M	2004	secondary progressive	periventricular WM
	3891	25	53	M	2005	chronic progressive	occipital WM
	3894	32	71	M	2004	secondary progressive	parietal WM
	3928	10	53	F	2004	secondary progressive	occipital WM
	4201	14	75	F	2006	Unclassified	occipital WM
	4212	19	50	F	2006	Unclassified	frontal WM
	4218	15	63	F	2006	Unclassified	frontal WM

Table 1: Characteristics of the frozen brain samples used in this study.
*Specimen used in the preceding RNA sequencing study [21].
NA=information not available from the participating brain banks
PMI=post-mortem interval (i.e., time elapsed between death and collection of the tissue specimen), rounded to the nearest hour

b. HERV database and Primer design

The retroviral gene catalog (RVGC) was compiled by aligning the protein domains in the Gypsy 2.0 database with the human genome (NCBI build 37.p13) [22]. The RVGC contains human genomic sequences with detectable protein homology (BLASTX) to entries in GyDB 2.0 domains database: core (GAG), envelope (ENV), reverse transcriptase (RT), integrase (INT) SCAN-A2, KRAB-A, dUTPase, RNase H, protease (AP) and chromo domain (CHR) (Supplementary Data). RVCG entries were named with arbitrary unique codes of the form: _ U (Table S1).

Deep sequencing showed that some HERV domains including GAG and ENV were overexpressed in MS samples [20,21]. Based on these results, the 28 GAG and 88 ENV gene sequences in the RVGC were aligned using Clustal Omega [23]. The direction of each sequence was determined by finding the longest open reading frame of the sequence. Before alignment all domain sequences were rectified to the sense orientation if necessary. GAG and ENV phylogenies were constructed using neighbor-joining (Figures 1 and 2). The major nodes within the phylogenies were used to define clades within the GAG and ENV domains (Figures 1 and 2).

In order to amplify multiple similar GAG and ENV loci spread across the human genome, primers were designed to be clade-specific, rather than locus specific. For this purpose, clade consensus sequences from the Clustal Omega multiple alignments were used to derive qPCR primer pairs, using Primer3 [24]. These primer pairs were 20-24 bp in length, with annealing temperatures of 58-60°C, with no predicted low-energy self-dimers, predicted to give 100-120bp amplicons. The primer sequences were aligned to the human genome using BLAST to predict possible amplicons [25]. Primers predicted to specifically bind the clade-specific consensus sequence were selected to minimize their hybridization to the other clades. Several primer pairs were designed for each of the GAG and ENV clades. The primer sequences used in the study are specified in Table 2.

Domain	Primer Pair	Forward sequence (5’-3’)	Reverse Sequence (5’-3’)
GAG	Clade 1A	GCAACATCTTGGAGCCTTG	ATCTGCCTTAGAGCCTGGGA
	Clade 1B	ATCTTGGAGCCTTGCCACAA	ATCTGCCTTAGAGCCTGGGA
	Clade 1C	AACATCTTGGAGCCTTGCCA	ATCTGCCTTAGAGCCTGGGA
	Clade 2	GCCTTCACATATTCTGTAATT	GCAAGGTTGCAAGATGCAGCTC
	Clade 3A	CAAATGCCCCTGAGAGAGCA	TTCCAGTTGAGAAGGTCGGC
	Clade 3B	CCACCAGATAACCACACCCC	TGCTCTCTCAGGGGCATTTG
	Clade 3C	TTAAGAGGACAGGCAGCAGC	GGGGTGTGGTTATCTGGTGG
ENV	Clade 1A	TGGGAAAACAGAATGGCCCT	AAGGCCCTTGTTATGCTCCC
	Clade 1B	TCCATGGACTGACAGTGGGA	AAAGCAATTCCAGCAGCCAC
	Clade 2	AGACCCACCAGTCAAATGGC	TTTGGTGTGATCCCACCTGG
	Clade 3A	GCGGTAAGTCTTTGCAAGGC	CCAGAATGGCTCCAGACATGT
	Clade 3B	CTGGACAATCAGGACCCACC	GGCCTTGCAAAGCCTTTGTT
	MSRV	CTTCCAGAATTGAAGCTGTAAAGC	GGGTTGTGCAGTTGAGATTTCC
Control	RPL19	ATGTATCACAGCCTGTACCTG	TTCTTGGTCTCTTCCTCCTTG
	RPL13	CCTGGAGGAGAAGAGGAAAGAGA	TTGAGGACCTCTGTGTATTTGTCAA
	UBC	ATTTGGGTCGCGGTTCTTG	TGCCTTGACATTCTCGATGGT

Table 2: Primers designed for this study

c. RNA extraction and cDNA synthesis

Study specimens from the frozen blocks were obtained using 6 mm sterile circular skin biopsy punch blades (VWR, Radnor, PA, USA). The tissue specimens (each approximately 100 mg) were homogenized by finely chopping them with a sterile blade in a sterile weigh boat, all on dry ice, before placing them into RNA-lysis buffer. The RNA was extracted using the RN easy Lipid Tissue Mini Kit (Qiagen, Hilden, Germany and Germantown, MD, USA) and quantified using an RNA High Sensitivity (HS) Assay Kit (Life Technologies/Thermo Fisher Scientific, Carlsbad, CA, USA) on a Qubit 2.0 Fluorometer instrument (Invitrogen/Thermo Fisher Scientific, Carlsbad, CA, USA). Extracted RNA was DNase treated using a TURBO DNA-free Kit (Life Technologies/Thermo Fisher Scientific, Carlsbad, CA, USA). Nucleic acids in the DNase treated specimens were quantified by using a Qubitds DNA HS Assay Kit. A total of 500 ng of RNA was used for the reverse transcription (RT) reaction using Omniscript RT (Qiagen, Hilden, Germany and Germantown, MD, USA), using random hexamers according to the manufacturer’s instructions. The resulting cDNA was aliquoted and stored at –20°C.

Several reference genes, previously used in studies of gene expression in neuronal tissue, were considered for use in this study: ubiquitin C (UBC), RNA polymerase II polypeptide (RP II), ribosomal large subunit protein 19 (RPL19), β-2-microglobulin (B2M), and glyceraldehyde-3- phosphate dehydrogenase (GAPDH) [26,27]. RNA expression levels for each of the candidate reference genes were determined by qPCR of the MS and control specimens. The geNorm application was used to determine the best reference gene [28]. RPL19 was selected as the reference gene due to the least variation in expression levels among all the specimens.

Five primer pairs were designed from consensus sequences for each of the 3 GAG and 3 ENV clades, for a total of 30 primer pairs. Melting curves from qPCR reaction were examined to determine suitable primer pairs to carry forward in the study. Specifically, 12 primer pairs demonstrating single homogenous amplification products were used for the study. Standard curves for each primer pair including reference gene candidates are shown in supporting data.

d. qPCR and Data analysis

Each reaction was run in duplicate with no-template, and no-RT controls. Average Ct value was considered for analysis according to the MIQE guidelines [29]. Control and MS samples (5 µl each) were pooled together and serially diluted 2-fold. This dilution series was used to obtain Ct standard curves for each gene of interest (Supplementary Data). Calibration curves with slopes in the range of .95-1.29 and R² <0.95 were considered for analysis. For each plate, PCR amplification efficiency was determined from the standard curve using the equation: PCR efficiency=(2^-1/slope-1) × 100% (Supplementary Figures). Gene expression values were calculated from the standard curve for each control and MS specimen, normalized by RPL19 expression for each sample. The median expression values for each primer pair was derived for the control and MS groups. The expression change (fold change) was derived by dividing the median expression value of the MS group from the median of the control group. The Mann-Whitney U-test was used to detect differences among the groups [30]. We considered the difference to be statistically significant when p< 0.05.

Results

a. Frozen Brain samples

Characteristics of the study samples are shown in Table 1. Thirty-three MS samples were compared to 9 normal controls. The number of control specimens was lower than MS samples due to limited availability of normal controls from the brain bank sources. The post-mortem interval (PMI), i.e., the time between death and collection of the brain sample, range was 2-39 hours. The PMI was significantly longer in the MS group (µ=17 hrs) than in the controls (µ=9.5 hrs, p=0.0002). Age, sex, and year of collection were not significantly different between the groups.

b. HERV GAG and ENV Phylogenies

Figure 1 shows the related GAG sequences classified by phylogeny into three different clades. According to Gypsy database annotation these clades include GAGs which belong to family or subfamily HERV-K10, K-HERV and HERV-E (specified in Table S1). Phylogeny using a neighbor joining tree of the ENV sequences also revealed three separate clades (Figure 2). Env sequence of HERV-W family; MSRV/HERV-W (AF331500.1), was also included in the ENV phylogeny [31]. The ENV sequences belonged to families including HERV-K10, K-HERV, HERV-E, Retrovirus-like Element-Isoleucine (RTVL-Ia), Mouse Mammary Tumor Virus (MMTV), and the D-type retrovirus Squirrel Monkey-like Retrovirus (SMRV-H) characterized from human B lymphocytes [32,33] (Table S1).

Figure 1: GAG sequence phylogeny-Twenty-eight non redundant GAG sequences were pulled from the RVGC (see Methods). Phylogeny was constructed using a neighbor-joining algorithm from the distance matrix. GAG sequences were classified into three different clades outlined by the boxes. GAG sequences over expressed in MS compared to controls and as determined by RNA-seq are displayed in red font.

Figure 2: ENV sequence phylogeny-Eighty-eight non redundant ENV sequences were pulled from the RVGC (see Methods). Phylogeny was constructed using a neighbor-joining algorithm from the distance matrix. ENV sequences were classified into three different clades outlined by the boxes. ENV sequences over expressed in MS compared to controls and as determined by RNA-seq are displayed in red font.

c. Expression analysis

Twelve clade-specific primer pairs were tested for the GAG and ENV domains. The MSRV-ENV primer sequences were kindly provided by H. Perron [32]. Each primer pair was tested by performing qPCR with on the pooled RNA extractions (from all control and MS samples). In silico analysis of these primer pairs was performed using Mega BLAST. This showed a heterogeneous population of predicted amplicons that mapped to multiple different human chromosomal locations specified in Table S1.

RPL19 was selected as the reference gene for normalization of expression because it had the lowest expression variation among the samples. Normalized gene expression levels for each of the targeted HERV GAG domains are shown in Figure 3. The data is summarized as relative fold change (log₂ ) for each primer pair. Expression of GAG clade 1A, 3B, and 3C were significantly increased while GAG clade 3A was decreased in the MS samples (N=33) compared to controls. The pattern of GAG expression for the MS subtypes among each of the clades was similar to that for the MS group as a whole (Figure S1). The magnitude of the significant changes in GAG expression was small with fold-changes only 0.74-1.25.

Figure 3: GAG domain expression in MS vs control-Differences in normalized expression for each primer pair are shown as Log2 (MS median/control median). Significance in expression differences was determined using the nonparametric Mann-Whitney test comparing expression values for the MS patients (N=33) to the controls (N=9). p-values for the statistically significant comparisons are shown on the chart.

Expression differences between the MS and control groups among the ENV clades studied were also small, with fold-changes less than 1.5 (Figure 4). No differences were seen in the expression of ENV clade 1. The expression of ENV clade 2 (fold-change 1.33), clades 3A and 3B (foldchange 1.2-1.27), and ENV MSRV/HERV-W (fold-change 1.15) were all increased in the MS group compared to controls. The pattern of ENV expression for the MS subtypes among each of the clades was similar to that for the MS group as a whole (Figure S2).

Figure 4: ENV domain expression in MS vs control-Differences in normalized expression for each primer pair are shown as Log₂ (MS median/control median). Significance in expression differences was determined using the nonparametric Mann-Whitney test comparing expression values for the MS patients (N=33) to the controls (N=9). p value for statistically significant groups is mentioned.

Discussion

This study compared HERV GAG and ENV domain expression measured by qPCR in frozen brain tissue specimens collected from individuals with MS compared to neurologically normal controls. While this study showed several significant differences in HERV expression, the magnitude of these differences was relatively small (less than 1.5-fold). These results indicate that the expression of HERV GAG and ENV regions do not differ greatly between different types of MS and controls in these frozen brain samples.

There are approximately 40,000 HERV elements distributed among thirty five human HERV families classified by Repbase [14,34]. Transcriptional activity of representative members of 20 HERV families in several different normal human tissues has been studied by microarray [34]. Our attempts to clone and amplify specific HERV loci of interest from the cDNA derived from these frozen brain samples was not successful, likely due to the highly homologous nature of the HERV sequences we were trying to study (data not shown). For this reason, a more unconventional approach was employed utilizing alignments among related HERVs (clades).

Specifically, this study used an alignment approach, employing HERV phylogenies (Figures 1 and 2) to design consensus primer sequences (Table 2) amplifying multiple different HERV ENV and GAG regions. These regions, designated as clades and subclades in this study, were derived from 28 GAG and 88 ENV domain sequences distributed among many different HERV families (Table S1). The approach we took does not formally exclude the possibility that some very specific HERV loci might very elevated in MS brain tissue. Rather, we were interested in getting a broad overview of HERV expression in MS brain tissue, complementing our previous sequencing data.

HERVs constitute about 8% of the human genome [35]. They are scattered widely throughout the human genome, and most cannot produce virus particles. HERVs have recently been shown to act as tissue-specific enhancers of gene transcription [36], and some human proteins are the result of HERV “domestication.” One example is syncytin, an important protein that allows for development of the placenta [37,38]. HERVs have been proposed as etiological cofactors in chronic diseases such as cancer, autoimmune disorders, and neurological disease [39]. Two complete ENV ORFs from the HERV-W family are thought to be associated with MS: MSRV, and ERVWE1 (syncytin-1) [10].

Our group recently reported an RNA-seq study showing that some HERVs and other retro elements are significantly over expressed in a smaller set of frozen MS and control brain tissue samples. That study found over expression of GAG sequences mapping to GAG Clade 1 (HERV-K10, K-HERV) and ENV sequences mapping to ENV clades 2 ( RTVL-Ia) and 3 (K-HERV) (Figures 1 and 2). However, these findings were partially confirmed by the present qPCR study where expression of ENV clades 3A and 3B (HERV-K and HERV-E) were increased in MS. Other groups have reported increased levels of ENV MSRV (HERV-W) in blood, spinal fluid, and brain samples in MS patients [17]. We did confirm that finding here, but the magnitude of the MSRV expression increase was relatively small – only 1.15 fold-change).

A number of assumptions are inherent to this study. First, HERV GAG and ENV regions were chosen for evaluation because multiple reports show that these are the retrovirus genes most likely to be involved in immunopathogenesis [17,40-43]. Second, we assume that the collapse of GAG and ENV genes from the Gypsy 2.0 database into the RVGC is sufficiently representative of all GAG and ENV HERV domains. Third, we assume that the primer pairs designed from the HERV GAG and ENV phylogenies and divided into related sequences (clades), providing a reasonable representation of these domains’ RNA within the samples. Finally, we assume that these postmortem samples, typically taken several hours after death, accurately reflect the premortem RNA content of the tissue. There is no doubt that cellular RNA is not static after death, and this clearly a limitation of our study, but fresh brain samples from neurologically-normal or diseased people are simply not available for comparison. Our MS samples had significantly longer PMIs than the controls, and this could have decreased the expression measured in the MS group more than in the controls, leading to an underestimation of expression differences between the groups. However, any PMI effect on expression should be canceled by RPL19 normalization which, presumably, deteriorates at a similar rate to the HERV domains we measured. All these assumptions can be criticized, but the goal of the study was to provide a broad view of HERV expression in MS white matter samples.

Conclusion

Expression of GAG and ENV mapping to HERV-E, RTVL and HERV-K10 families were significantly increased in the MS group. However, the relative expression differences between the MS and control groups were small, differing less than 1.5-fold. These results indicate that the expression of HERV GAG and ENV regions do not differ greatly between MS and controls in these frozen brain samples. This work only partially confirms the sequencing data previous derived from similar frozen MS brain specimens.

Financial Disclosure

This work was funded by grant R21NS077023 from the NIH/NINDS, “Deep Sequencing for the Detection of Viral Sequences in Primary Progressive Multiple Sclerosis Brains.” Some follow up work was funded by the National MS Society, grant #RG4807A4.

References

1. Sawcer S, Franklin RJ, Ban M (2014) Multiple sclerosis genetics. Lancet Neurol 13: 700-709. [Ref.]
2. Cree BA (2014) Multiple sclerosis genetics. Handb clin neurol 122: 193-209. [Ref.]
3. Alcina A, Abad-Grau Mdel M, Fedetz M, Izquierdo G, Lucas M, et al. (2012) Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations. PloS one 7: e29819. [Ref.]
4. Kurtzke J (2005) Epidemiology and etiology of multiple sclerosis. Phys Med Rehabil Clin N Am 16: 327-349. [Ref.]
5. Meinl E (1999) Concepts of viral pathogesis of MS. Current Opinion in Neurology 12: 303-307. [Ref.]
6. Murray J (2002) Infection as a cause of multiple sclerosis. BMJ 325: 1128. [Ref.]
7. Garson JA, Tuke PW, Giraud P, Paranhos-Baccala G, Perron H (1998) Detection of virion-associated MSRV-RNA in serum of patients with multiple sclerosis. Lancet 351: 33. [Ref.]
8. Perron H, Geny C, Laurent A, Mouriquand C, Pellat J, et al. (1989) Leptomeningeal cell line from multiple sclerosis with reverse transcriptase activity and viral particles. Res Virol 140: 551-561. [Ref.]
9. Perron H, Suh M, Lalande B, Gratacap B, Laurent A, et al. (1993) Herpes simplex virus ICP0 and ICP4 immediate early proteins strongly enhance expression of a retrovirus harboured by a leptomeningeal cell line from a patient with multiple sclerosis. J Gen Virol 74: 65-72. [Ref.]
10. Dolei A, Garson JA, Arru G, Clerici M, Germi R, et al. (2014) Multiple sclerosis-associated retrovirus and related human endogenous retrovirus-W in patients with multiple sclerosis. Journal of neuroimmunology 266: 87- 88. [Ref.]
11. Douville RN, Nath A (2014) Human endogenous retroviruses and the nervous system. Handb Clin Neurol 123: 465-485. [Ref.]
12. Paces J, Pavlicek A, Paces V (2002) HERVd: database of human endogenous retroviruses. Nucleic acids research 30: 205-206. [Ref.]
13. Paces J, Pavlicek A, Zika R, Kapitonov VV, Jurka J, et al. (2004) HERVd: the Human Endogenous RetroViruses Database: update. Nucleic acids research 32: D50. [Ref.]
14. Kim TH, Jeon YJ, Yi JM, Kim DS, Huh JW, et al. (2004) The distribution and expression of HERV families in the human genome. Mol Cells 18: 87-93. [Ref.]
15. Belshaw R, Katzourakis A, Paces J, Burt A, Tristem M (2005) High copy number in human endogenous retrovirus families is associated with copying mechanisms in addition to reinfection. Mol Biol Evol 22: 814-817. [Ref.]
16. Antony JM, Zhu Y, Izad M, Warren KG, Vodjgani M, et al. (2007) Comparative expression of human endogenous retrovirus-W genes in multiple sclerosis. AIDS Res Hum Retroviruses 23: 1251-1256.
17. Perron H, Lang A (2010) The human endogenous retrovirus link between genes and environment in multiple sclerosis and in multifactorial diseases associating neuroinflammation. Clin Rev Allergy Immunol 39: 51-61. [Ref.]
18. Perron H, Lazarini F, Ruprecht K, Pechoux-Longin C, Seilhean D, et al. (2005) Human endogenous retrovirus (HERV)-W ENV and GAG proteins: physiological expression in human brain and pathophysiological modulation in multiple sclerosis lesions. J Neurovirol 11: 23-33. [Ref.]
19. Vogetseder W, Dumfahrt A, Mayersbach P, Schonitzer D, Dierich MP (1993) Antibodies in human sera recognizing a recombinant outer membrane protein encoded by the envelope gene of the human endogenous retrovirus K. AIDS Res Hum Retroviruses 9: 687-694. [Ref.]
20. Nath A, Küry P, Olival GSd, Dolei A, Karlsson H, et al. (2015) First international workshop on human endogenous retroviruses and diseases, HERVs & disease 2015. Mob DNA 6: 20. [Ref.]
21. Kriesel JD, Bhetariya PJ, Chan BK, Wilson T, Fischer KF (2014) Enrichment of Retroviral Sequences in Brain Tissue from Patients with Progressive Multiple Sclerosis. ECTRIMS Online Library: 64391. [Ref.]
22. Llorens C, Futami R, Covelli L, Dominguez-Escriba L, Viu JM, et al. (2011) The Gypsy Database (GyDB) of mobile genetic elements: release 2.0. Nucleic Acids Res 39: D70-D74. [Ref.]
23. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, et al. (2011) Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7: 539. [Ref.]
24. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, et al. (2012) Primer3--new capabilities and interfaces. Nucleic acids Res 40: e115. [Ref.]
25. Morgulis A, Coulouris G, Raytselis Y, Madden TL, Agarwala R, et al. (2008) Database indexing for production MegaBLAST searches. Bioinformatics 24: 1757-1764. [Ref.]
26. Penna I, Vella S, Gigoni A, Russo C, Cancedda R, et al. (2011) Selection of candidate housekeeping genes for normalization in human postmortem brain samples. Int J Mol Sci 12: 5461-5470. [Ref.]
27. Coulson DT, Brockbank S, Quinn JG, Murphy S, Ravid R, et al. (2008) Identification of valid reference genes for the normalization of RT qPCR gene expression data in human brain tissue. BMC mol biol 9: 46. [Ref.]
28. Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, et al. (2009) A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol 10: R64. [Ref.]
29. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, et al. (2009) The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 55: 611-622. [Ref.]
30. Lowry R (2016) VassarStats: Website for Statistical Computation: Vassar College open source web application for statistical testing. [Ref.]
31. Mameli G, Poddighe L, Astone V, Delogu G, Arru G, et al. (2009) 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 161: 98-106. [Ref.]
32. Oda T, Ikeda S, Watanabe S, Hatsushika M, Akiyama K, et al. (1988) Molecular cloning, complete nucleotide sequence, and gene structure of the provirus genome of a retrovirus produced in a human lymphoblastoid cell line. Virology 167: 468-476. [Ref.]
33. Mason AL, Gilady SY, Mackey JR (2011) Mouse Mammary Tumor Virus in Human Breast Cancer: Red Herring or Smoking Gun? Am J Pathol 179:1588-1590. [Ref.]
34. Seifarth W, Frank O, Zeilfelder U, Spiess B, Greenwood AD, et al. (2005) Comprehensive analysis of human endogenous retrovirus transcriptional activity in human tissues with a retrovirus-specific microarray. J virol 79: 341-352. [Ref.]
35. Dolei A (2006) Endogenous retroviruses and human disease. Expert Rev Clin Immunol 2: 149-167. [Ref.]
36. Pavlicev M, Hiratsuka K, Swaggart KA, Dunn C, Muglia L (2015) Detecting endogenous retrovirus-driven tissue-specific gene transcription. Genome biology and evolution 7: 1082-1097. [Ref.]
37. Potgens AJ, Drewlo S, Kokozidou M, Kaufmann P (2004) Syncytin: the major regulator of trophoblast fusion? Recent developments and hypotheses on its action. Human Reprod update 10: 487-496. [Ref.]
38. Mi S, Lee X, Li X, Veldman GM, Finnerty H, et al. (2000) Syncytin is a captive retroviral envelope protein involved in human placental morphogenesis. Nature 403: 785-789. [Ref.]
39. Antony JM, Deslauriers AM, Bhat RK, Ellestad KK, Power C (2011) Human endogenous retroviruses and multiple sclerosis: innocent bystanders or disease determinants? Biochim biophysic acta 1812: 162-176. [Ref.]
40. Christensen T, Sorensen PD, Hansen HJ, Moller-Larsen A (2003) Antibodies against a human endogenous retrovirus and the preponderance of env splice variants in multiple sclerosis patients. Mult Scler 9: 6-15. [Ref.]
41. Clerici M, Fusi ML, Caputo D, Guerini FR, Trabattoni D, et al. (1999) Immune responses to antigens of human endogenous retroviruses in patients with acute or stable multiple sclerosis. Journal of neuroimmunology 99: 173-182. [Ref.]
42. Perron H, Bernard C, Bertrand JB, Lang AB, Popa I, et al. (2009) Endogenous retroviral genes, Herpesviruses and gender in Multiple Sclerosis. J Neurol Sci 286: 65-72. [Ref.]
43. Perron H, Germi R, Bernard C, Garcia-Montojo M, Deluen C, et al. (2012) Human endogenous retrovirus type W envelope expression in blood and brain cells provides new insights into multiple sclerosis disease. Mult Scler 18: 1721-1736. [Ref.]

Download Supplementary Data Here

 

Download Provisional PDF Here

 

Article Information

Aritcle Type: Research Article

Citation: Bhetariya PJ, Kriesel JD, Fischer KF (2017) Analysis of Human Endogenous Retrovirus Expression in Multiple Sclerosis Plaques. J Emerg Dis Virol 3(2): doi http://dx.doi.org/10.16966/2473-1846.133

Copyright: © 2017 Bhetariya PJ, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Publication history: 

Received date: 07 Jun 2017

Accepted date: 18 Jul 2017

Published date: 24 Jul 2017