Quantitative analysis of human endogenous retrovirus-K transcripts in postmortem premotor cortex fails to confirm elevated expression of HERV-K RNA in amyotrophic lateral sclerosis
Acta Neuropathologica Communicationsvolume 7, Article number: 45 (2019)
The Letter to the Editor to this article has been published in Acta Neuropathologica Communications 2019 7:101
Over the past two decades a number of studies have demonstrated activity of the retroviral enzyme reverse transcriptase in the serum of patients with sporadic amyotrophic lateral sclerosis (ALS). Known human exogenous retroviruses such as HIV-1 have been eliminated as possible sources of this activity and investigators have therefore considered the possibility that human endogenous retroviruses (HERVs) might be involved. HERV-K (HML-2) is the most recent retroviral candidate to be proposed following the observation of elevated HERV-K expression in cortical and spinal neurons of ALS patients and the demonstration of HERV-K envelope protein neurotoxicity in vitro and in transgenic mice. This retroviral hypothesis is an attractive one, not least because it raises the possibility that ALS might become treatable using antiretroviral drugs. In the present study we have attempted independent confirmation of the observation that HERV-K RNA levels are elevated in ALS brain. Total RNA was extracted from the postmortem premotor cortex of 34 patients with ALS and 23 controls. Quantitative real-time reverse transcription PCR (RT-qPCR) was performed according to the MIQE guidelines using HERV-K gag, pol and env primer sets. Data was analysed by the 2-∆∆Ct method with normalisation against two reference genes, GAPDH and XPNPEP1. Geometric mean HERV-K RNA expression levels in the premotor cortex of ALS patients were not found to be different from the expression levels in non-ALS controls. Our findings do not confirm the recently reported association between elevated cortical HERV-K RNA levels and ALS, and thus raise doubts about the role of this endogenous retrovirus in ALS pathogenesis. The results of this study may have implications for ongoing clinical trials aiming to suppress HERV-K activity with antiretroviral drugs.
Amyotrophic lateral sclerosis (ALS), also known as motor neuron disease, is a fatal neurodegenerative disease characterised by loss of motor neurons from the brain and spinal cord. In Europe and the USA the incidence is about 2 cases per 100,000 person-years and survival post diagnosis is typically 3 to 5 years. In about 5–10% a family history of ALS may be obtained in a first or second degree relative. Variants causative of or associated with ALS have been identified in at least 25 different genes, and in people with or without a family history, but the aetiology of most apparently sporadic cases remains uncertain .
Numerous non-genetic environmental risk factors for apparently sporadic ALS have been proposed including electromagnetic fields, heavy metals, pesticides, smoking, dietary factors, physical trauma and viral infection . Enteroviruses and herpesviruses have been considered [5, 8] but more recent attention has focused on retroviruses. Retroviruses are known to cause motor neuron disease in mice , and the human retroviruses HIV-1 and HTLV-1 are both capable of causing ALS-like syndromes [23, 32], which in some cases have been shown to respond well to antiretroviral therapy [1, 22].
In a series of studies looking directly for evidence of retroviral involvement in ALS we have been able to exclude known exogenous retroviruses whilst repeatedly demonstrating an increased prevalence and raised levels of reverse transcriptase activity (a generic retroviral marker) in the serum of patients with ALS [3, 4, 25, 29]. The increased prevalence of reverse transcriptase activity in ALS was subsequently confirmed independently by another group . In one of our studies  reverse transcriptase activity was detected more frequently in the serum of unaffected blood relatives of ALS patients than in unrelated controls and spouses. This raised the possibility that the observed reverse transcriptase activity might be associated with an inherited endogenous retrovirus.
Human endogenous retroviruses (HERVs), which constitute around 8% of the human genome, are thought to be the relics of retroviral germline infections that occurred millions of years ago . Most HERVs are considered inactive due to the accumulation of mutations and deletions but there is increasing evidence that some of them may be capable of expressing full-length RNA transcripts, proteins and even retroviral particles. The most recently integrated HERVs such as HERV-K (HML-2) are thought to be the most intact and potentially biologically active  and it was therefore of great interest when elevated HERV-K RNA levels were reported in the cerebral cortex of patients with ALS [9, 19]. In addition to increased cortical HERV-K RNA, Li et al.  reported that HERV-K envelope protein was selectively expressed in cortical and spinal neurons of ALS patients and that the envelope protein was neurotoxic in stem-cell derived human neurons and in a transgenic mouse ALS model, strongly suggesting that HERV-K contributes to motor neuron disease.
In the present study we have attempted independent confirmation of elevated HERV-K transcripts in the cortex of ALS patients by using exactly the same GAPDH (Glyceraldehyde 3-phosphate dehydrogenase)-normalised RT-qPCR methods for quantification of HERV-K gag, pol and env transcripts as those described previously . We have also repeated the quantification of cortical HERV-K transcripts using an alternative validated reference gene, XPNPEP1 (X-prolyl aminopeptidase 1), as recommended by the MIQE guidelines . Additionally, we investigated HERV-W, which has previously been associated with multiple sclerosis, schizophrenia and chronic inflammatory demyelinating polyneuropathy , in the same ALS and control samples by estimating its env RNA expression using a similar RT-qPCR method, to exclude it as a possible cause of the previously observed raised levels of serum reverse transcriptase in ALS.
Materials and methods
Frozen postmortem brain material was obtained from the Medical Research Council (MRC) Neurodegenerative Disease Brain Bank Network, Institute of Psychiatry, Psychology and Neuroscience, Kings College London, UK. Premotor cortex (i.e. part of the frontal lobe just anterior to the primary motor cortex) from 34 patients with sporadic ALS and 23 non-ALS controls was analysed. Details including age, gender, diagnosis, postmortem delay and RNA integrity number (RIN) are presented in Additional file 1: Table S1. All ALS patients had their ALS clinical diagnosis confirmed by neuropathological examination of the brain postmortem.
RNA extraction and quality control
60 mg pieces of frozen (− 80 °C) premotor cortex were homogenised (TissueRuptor II, Qiagen Ltd., Crawley, UK) on dry ice in QIAzol lysis buffer (Qiagen) and extracted using the RNeasy Lipid Tissue kit (Qiagen) with on-column DNase treatment (RNase-free DNase Set, Qiagen) according to the manufacturer’s instructions. RNA concentration was measured by Qubit™ RNA BR assay (Thermo Fisher Scientific Inc. Waltham, MA, USA) and RNA quality established by Agilent RNA 6000 Nano assay (Agilent Technologies, Inc. Santa Clara, CA, USA).
Reverse transcription of 1 μg of total RNA in a 20 μl reaction was performed using the Invitrogen SuperScript III First-Strand Synthesis Supermix for qRT-PCR (Thermo Fisher) according to manufacturer’s instructions. RNase H digestion was carried out at 37 °C for 20 mins and cDNA stored at − 20 °C. Negative control reverse transcription reactions without adding RNA were included in each batch of cDNA syntheses.
Real-time PCR was performed in an Applied Biosystems QuantStudio™ 5 thermocycler (Thermo Fisher), 96 well format, using Fast SYBR Green Master Mix (Thermo Fisher) in a 20 μl reaction with 2 μL cDNA and final concentrations of 0.25 μM forward primer and 0.25 μM reverse primer or 1X primer pool for XPNPEP1. Primer sequences for amplification of HERV-K gag, pol and env, and GAPDH were as described previously  and detailed in Table 1. Thermal cycling parameters were as follows: DNA polymerase activation, 95 °C for 20 s, then 45 cycles of denaturation, 95 °C for 1 s and anneal/extension 60 °C for 20 s. Cycling was followed by melt curve analysis. Automatic baseline settings were used with a manual threshold setting of 0.2 in all experiments. Each sample was analysed in duplicate and any sample with a Cq (quantification cycle) standard deviation of > 0.2 was repeated. All 96 well plates contained samples from both ALS patients and control individuals. No template controls (NTCs) were run in each experiment and every sample was analysed with and without reverse transcription in order to identify any residual genomic DNA contamination in the extracted RNA. Relative HERV-K RNA expression levels (i.e. relative with respect to the geometric mean of the non-ALS controls) were calculated using the 2-∆∆Ct method  by normalisation against both GAPDH and XPNPEP1 reference transcripts. RT-qPCR procedures were performed in accordance with the MIQE guidelines  (Additional file 1: Table S2) and all experiments were conducted ‘blind’ with the identity of each sample being hidden from the investigator. RT-qPCR experiments were conducted simultaneously on ALS cases and controls so as to eliminate such confounders as batch effects.
Reference gene validation
A panel of 9 candidate reference gene primer sets were evaluated in order to establish the one with the most stable level of expression in ALS and non-ALS control premotor cortex samples. Briefly, RT-qPCR was used to quantify the RNA expression levels of GAPDH, ACTB, CYC1, SDHA, UBC, RPL13A, XPNPEP1, EIF4A2 and YWHAZ in 5 ALS samples and 5 non-ALS control samples in triplicate. All primer sets apart from XPNPEP1 were obtained from Primerdesign Ltd., Camberley, UK. The evaluation was conducted using qBase+ software, version 3.1 (Biogazelle, Zwijnaarde, Belgium), which utilises the geNorm selection algorithm , along with additional statistical tools. Further verification was performed using RefFinder software  which exploits the computational programs Normfinder , BestKeeper , geNorm and the ΔCt method  to comprehensively rank and compare candidate reference genes.
Amplicons generated by the HERV-K gag, HERV-K pol, HERV-K env, HERV-W env, GAPDH and XPNPEP1 PCRs were subjected to Sanger dideoxy sequencing (Eurofins GTAC Biotech, Germany) and sequences analysed by BLAST (Basic Local Alignment Search Tool) in order to check the specificity of each assay. In each case amplicon size was confirmed by agarose gel electrophoresis.
The statistical significance of differences between groups was assessed by 2-tailed Mann-Whitney U test (GraphPad Prism 7 software). p-values of < 0.05 were considered significant. Linear regression p-values were calculated using Microsoft Excel Add-In Daniel’s XL toolbox v6.22.
Matching of ALS patients with non-ALS control group
Due to limited availability of suitable postmortem material it was not always possible to obtain a perfect match of all clinical parameters between the ALS patient group and the non-ALS control group. The mean age at death was 66.9 years for the ALS patients and 73.5 years for the controls (p = 0.03). There was also a difference in gender distribution between the ALS and control groups. The ALS group was 29% female whereas the control group was 48% female. There was no significant difference in mean postmortem delay between the patients and controls; for the ALS group it was 45.1 h and for the control group 41.7 h (p = 0.64). Finally, there was a small difference between the mean RNA integrity value (RIN) of the ALS samples, 6.53, and that of the controls, 6.05 (p = 0.01). The potential effect of such imperfectly matched parameters on HERV-K expression comparisons between ALS patients and controls is considered below.
RT-qPCR performance characteristics
Agarose gel electrophoresis confirmed that the amplicons generated by each of the six different RT-qPCR assays (HERV-K gag, HERV-K pol, HERV-K env, HERV-W env, GAPDH, XPNPEP1) were single bands of the expected size. Melt curve analysis also revealed single dominant peaks for all assays. Amplification specificity was further confirmed by Sanger sequencing of the amplicons which demonstrated that the nucleotide sequence was as expected for all six RT-qPCR assays (Basic Local Alignment Search Tool, BLASTn, analysis). PCR efficiencies for each assay were estimated from slopes of serial dilution standard curves using the formula E = 10^(− 1/slope) – 1 and all fell within the range 99% ± 6%. For each assay the R2 value of the standard curve was > 0.99. The no template controls were negative in all experimental runs and the no reverse transcriptase controls revealed that residual genomic DNA in RNA extracts typically contributed less than 1%, but always less than 3%, of the total signal in all samples.
Reference gene validation
qBase+ and RefFinder ranking of the nine candidate reference genes revealed that XPNPEP1 and GAPDH had the most stable expression in ALS and control samples (Additional file 1: Table S3). XPNPEP1 was therefore selected as the additional ‘validated’ reference gene to use in addition to GAPDH for the normalisation of HERV-K and HERV-W RNA expression levels.
HERV-K RNA expression in ALS and controls
Figure 1 and Table 2 summarise the relative expression levels of HERV-K gag, pol and env RNA in the 34 ALS cases and 23 non-ALS controls investigated. When the GAPDH reference gene was used for normalisation, the geometric mean expression levels for HERV-K gag and HERV-K pol were very slightly higher in the ALS patients than in the control group but the converse was true for HERV-K env. None of these marginal differences in HERV-K expression levels approached statistical significance. When the XPNPEP1 reference gene was used for normalisation the geometric mean expression levels for all three HERV-K genes were slightly lower in ALS cases than in controls. Once again these small differences in geometric mean expression levels failed to reach statistical significance. However, the expression levels of HERV-K gag, pol and env were correlated with each other whether the data was normalised by GAPDH or XPNPEP1 reference genes (Fig. 2).
As perfect matching between the ALS patients and non-ALS controls of parameters such as age, gender, postmortem delay and RIN values was not always possible due to limited availability of suitable postmortem material, these variables were also examined for possible correlation with HERV-K expression. Neither age nor gender was correlated with HERV-K gag, pol or env expression levels whether data was normalised by GAPDH or XPNPEP1 reference genes (Additional file 1: Figure S1 and Figure S2). A trend towards reduced HERV-K RNA levels with increased postmortem delay was observed (Additional file 1: Figure S3) for HERV-K pol and env, whether normalisation was to GAPDH or XPNPEP1 reference genes (p ≤ 0.01). RIN values were not correlated with HERV-K gag, pol or env RNA levels when GAPDH normalisation was used but there was a slight negative correlation between high RIN values and XPNPEP1-normalised relative HERV-K gag and env RNA levels (gag, p = 0.04 and env, p = 0.03, Additional file 1: Figure S4).
HERV-W env RNA expression in ALS and controls
Figure 3 shows the relative HERV-W env RNA expression levels in ALS cases and non-ALS controls. When GAPDH was employed for normalisation, the geometric mean HERV-W env RNA level for the ALS cases was 0.87 and for the controls 1.00 (p = 0.26). However, when XPNPEP1 was used for normalisation the geometric mean HERV-W env RNA level in ALS was lower at 0.75 than in controls at 1.00 (p = 0.04).
HERV-W env RNA levels showed no correlation with age or gender whether normalised by GAPDH or XPNPEP1 (Additional file 1: Figure S5 and Figure S6). Postmortem delay was not significantly correlated with HERV-W env RNA levels when XPNPEP1 normalisation was used but there was a downward trend of HERV-W env RNA level with increased PMD when GAPDH normalisation was used (p = 0.04) (Additional file 1: Figure S7). With XPNPEP1 normalisation, high RIN values were correlated with low relative expression of HERV-W env RNA (p = 0.001) but there was no correlation when GAPDH normalisation was used (Additional file 1: Figure S8).
Correlation between GAPDH-normalised and XPNPEP1-normalised expression data
GAPDH-normalised relative expression levels and XPNPEP1-normalised relative expression levels correlated well with each other for HERV-K gag, pol and env, and also for HERV-W env transcripts (Additional file 1: Figure S9).
Overexpression of HERV-K has recently been proposed as a possible causative factor in patients with ALS . Although this is an attractive hypothesis, not least because it raises the possibility that ALS might become treatable using antiretroviral drugs  or antibodies , the contribution of HERV-K to the disease process has yet to be conclusively proven. The magnitude of the difference reported previously  between the mean HERV-K RNA expression level in ALS patients and non-ALS controls was less than threefold for gag and env, and less than two fold for pol. Although such relatively modest differences can in principle be resolved by RT-qPCR due to the high technical precision of the method, it is essential to demonstrate their reproducibility in other patient cohorts by independent testing such as that undertaken here. In this study we have therefore attempted confirmation of the observation that GAPDH-normalised HERV-K RNA levels are elevated in the cerebral cortex of ALS patients with respect to non-ALS controls. Postmortem premotor cortex samples from 34 patients with ALS and 23 non-ALS controls were tested using the same RT-qPCR methods (including the same reverse transcription method, the same DNase treatment method, the same PCR reagents and the same thermal cycling parameters), the same 2-∆∆Ct data analysis method and identical primer sets to those used previously . Our observations were concordant with previous findings  regarding the good correlation between the relative expression levels of the three HERV-K transcripts gag, pol and env, suggesting that the entire viral genome was expressed. However, in contrast to the findings of Li et al. , we were unable to demonstrate any difference in HERV-K gag, pol or env RNA levels between ALS patients and controls, whether the data were normalised by GAPDH or XPNPEP1 reference genes.
The reasons for this conflict are uncertain but several possibilities exist. Firstly, our ALS patients and non-ALS controls were not perfectly matched for age, gender, postmortem delay or RIN values and the differences between the groups were statistically significant for age, gender and RIN. However, none of these three parameters was found to be correlated with GAPDH-normalised HERV-K RNA expression and so their imperfect matching does not appear to explain why our results differ from those published previously . Secondly, it is conceivable that differences between the USA ALS cohort  and our UK ALS cohort might explain the discrepant results. We consider this unlikely because the diagnosis of ALS was confirmed by neuropathological examination in both cohorts and tissue samples in the two cohorts were derived from similar neuroanatomical regions. Nevertheless, it is theoretically possible that there is a small subset of ALS patients who do have elevated cortical HERV-K expression and that members of that subset were present in the ALS cohort studied by Li et al. but not in our cohort. However, in comparison with our ALS samples, the mean postmortem delay was shorter and the mean RIN value was higher in the samples analysed previously . Although we found no significant correlation between GAPDH-normalised HERV-K RNA expression and RIN we did observe a trend towards lower GAPDH-normalised HERV-K RNA levels (pol and env only) with increasing postmortem delay. Nevertheless, we consider it unlikely that this difference in mean postmortem delay explains our failure to confirm the previously published findings  because the HERV-K RNA RT-qPCR quantitative data is not ‘absolute’ but ‘relative’ with respect to the geometric mean level of the non-ALS controls which had postmortem delays not significantly different from the ALS cases. In agreement with a previous report on human postmortem brain RNA quality  we found no correlation between RIN value and postmortem delay (Additional file 1: Figure S10). Thirdly, it is conceivable that the discrepancy may be related to differences in the proportion of controls with cancer. Twenty five percent of the control group used by Li et al. had cancer as against 47% in the present study. Since upregulation of HERV-K expression occurs in various types of tumour tissue  it could be argued that there is a remote possibility that HERV-K expression could be elevated in the brains of the controls who had cancer, even without cerebral metastases. However, this explanation can be dismissed because our conclusion that patients with ALS do not have significantly higher levels of cerebral HERV-K RNA expression than non-ALS controls remains unchanged on reanalysing the data following exclusion of those with cancer from the control group (Additional file 1: Figure S11). Finally, the relatively low number of ALS cases studied by Li et al. (n = 11) may have resulted in less statistically robust findings than those generated by our larger analysis of 34 cases.
For ensuring accurate quantification, the importance of selecting appropriate reference genes for normalisation in RNA expression studies has been stressed by many authors, as has the advantage of using more than one reference gene [7, 15]. In the present study we evaluated a panel of nine candidate reference genes including GAPDH which had been used previously . This evaluation revealed that XPNPEP1 and GAPDH had the most stable expression in ALS and non-ALS control material, and XPNPEP1 was therefore chosen as the additional reference gene. It is noteworthy that Durrenberger et al.  in an evaluation of 12 candidate reference genes for use on human CNS postmortem tissue in various neurological diseases including ALS, also identified XPNPEP1 as the most stable and suitable candidate. GAPDH-normalised relative expression levels and XPNPEP1-normalised relative expression levels were generally well correlated. Our conclusion that HERV-K RNA levels are not elevated in the premotor cortex of ALS patients with respect to non-ALS controls was the same whichever reference gene was used for data normalisation.
In the study by Li et al. , cerebral cortical expression of a number of other HERVs, in addition to HERV-K, was measured. They examined HERV-E, HERV-R and HERV-P by RT-qPCR but did not observe significant ALS-associated elevation of transcripts in any of these. We considered that HERV-W should also be investigated because it had been implicated in a number of other neurological conditions including multiple sclerosis, schizophrenia and chronic inflammatory demyelinating polyneuropathy . Although our study does not show any association between HERV-W RNA expression level and ALS when data is GAPDH-normalised, a small negative correlation (p = 0.04) between HERV-W RNA expression and ALS was observed with XPNPEP1 normalisation due in part to four control samples showing relatively high HERV-W env expression when normalised to this reference gene. We speculate however that this small negative correlation may be unreliable due to our ALS samples having a slightly higher mean RIN value than controls, in conjunction with high RIN values being significantly correlated with lower XPNPEP1-normalised relative HERV-W expression.
Another study investigating the potential role of HERV-K in ALS has just been published by Mayer et al. . Using an alternative RT-qPCR method they examined levels of GAPDH-normalised HERV-K gag RNA in frozen postmortem neural tissues from a total of 108 ALS and control samples obtained from brain banks in the USA. As in our study, Mayer and colleagues were unable to confirm the previous findings  and concluded that levels of HERV-K transcripts were not significantly different between ALS and controls. Furthermore, they were unable to demonstrate any significant difference between ALS and controls in the pattern of transcriptionally active HERV-K loci by sequencing-based transcriptional profiling, or indeed any evidence of full length HERV-K envelope protein in either ALS or control tissues. These transcriptional profiling results also conflict with the findings of Li et al.  who observed three HERV-K loci to be transcribed at higher levels in ALS cases than in controls. However, it is important to note that the failure of our study and that of Mayer et al. to confirm elevated levels of HERV-K RNA in ALS post-mortem brain tissue does not diminish the significance of the report by Li et al.  that expression of HERV-K in human neurons in vitro causes toxicity, or that expression of HERV-K env in transgenic mice causes degeneration of motor neurons and progressive motor dysfunction.
In conclusion, our observations and those recently published by Mayer and colleagues  fail to confirm the findings of Li et al.  and provide no support for the hypothesis that elevated HERV-K expression in the cerebral cortex is associated with sporadic amyotrophic lateral sclerosis. These conflicting results may have significant implications for ALS clinical trials aiming to suppress HERV-K activity with antiretroviral drugs, and suggest that further research in this area is required to discover the source of the increased serum reverse transcriptase activity seen in this disease. Future studies that we propose include the use of next generation sequencing and custom-made microarrays to undertake a broad screening of the expression profiles of all human endogenous retrovirus families in ALS and controls.
Amyotrophic lateral sclerosis
quantification cycle (also referred to as Ct or threshold cycle)
- GAPDH :
Glyceraldehyde 3-phosphate dehydrogenase
human endogenous retrovirus-K (HML-2 subfamily)
human endogenous retrovirus-W
- MIQE guidelines:
Minimum Information for Publication of Quantitative Real-Time PCR Experiments
RNA integrity number
quantitative real-time reverse transcription polymerase chain reaction
- XPNPEP1 :
Alfahad T, Nath A (2013) Retroviruses and amyotrophic lateral sclerosis. Antivir Res 99:180–187.
Andersen CL, Jensen JL, Orntoft TF (2004) Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res 64:5245–5250.
Andrews WD, Al Chalabi A, Garson JA (1997) Lack of evidence for HTLV tax-rex DNA in motor neuron disease. J Neurol Sci 153:86–90.
Andrews WD, Tuke PW, Al-Chalabi A, Gaudin P, Ijaz S, Parton MJ, Garson JA (2000) Detection of reverse transcriptase activity in the serum of patients with motor neuron disease. J Med Virol 61:527–532.
Berger MM, Kopp N, Vital C, Redl B, Aymard M, Lina B (2000) Detection and cellular localization of enterovirus RNA sequences in spinal cord of patients with ALS. Neurology 54:20–25.
Brown RH, Al-Chalabi A (2017) Amyotrophic lateral sclerosis. N Engl J Med 377:162–172.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT (2009) The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 55:611–622.
Cermelli C, Vinceti M, Beretti F et al (2003) Risk of sporadic amyotrophic lateral sclerosis associated with seropositivity for herpesviruses and echovirus-7. Eur J Epidemiol 18:123–127.
Douville R, Liu J, Rothstein J, Nath A (2011) Identification of active loci of a human endogenous retrovirus in neurons of patients with amyotrophic lateral sclerosis. Ann Neurol 69:141–151.
Durrenberger PF, Fernando FS, Magliozzi R et al (2012) Selection of novel reference genes for use in the human central nervous system: a BrainNet Europe study. Acta Neuropathol 124:893–903.
Durrenberger PF, Fernando S, Kashefi SN et al (2010) Effects of antemortem and postmortem variables on human brain mRNA quality: a BrainNet Europe study. J Neuropathol Exp Neurol 69:70–81.
Gardner MB (1991) Retroviral leukemia and lower motor neuron disease in wild mice: natural history, pathogenesis, and genetic resistance. Adv Neurol 56:473–479.
GeNeuro (2018) Novel pHERV-K Antibodies for ALS in collaboration with NINDS. Press release. http://www.geneuro.com/data/news/GeNeuro-PR-NIH-ALS-License-Press-release-ENG.pdf Accessed 17 Dec 2018.
Gröger V, Cynis H (2018) Human endogenous retroviruses and their putative role in the development of autoimmune disorders such as multiple sclerosis. Front Microbiol 9(265). https://doi.org/10.3389/fmicb.2018.00265.
Huggett J, Dheda K, Bustin S, Zumla A (2005) Real-time RT-PCR normalisation; strategies and considerations. Genes Immun 6:279–284.
Ingre C, Roos PM, Piehl F, Kamel F, Fang F (2015) Risk factors for amyotrophic lateral sclerosis. Clin Epidemiol 7:181–193.
Küry P, Nath A, Créange A, Dolei A, Marche P, Gold J, Giovannoni G, Hartung HP, Perron H (2018) Human endogenous retroviruses in neurological diseases. Trends Mol Med 24:379–394.
Levet S, Medina J, Joanou J, Demolder A, Queruel N, Réant K, Normand M, Seffals M, Dimier J, Germi R, Piofczyk T, Portoukalian J, Touraine JL, Perron H (2017) An ancestral retroviral protein identified as a therapeutic target in type-1 diabetes. JCI Insight 2:e94387. https://doi.org/10.1172/jci.insight.94387.
Li W, Lee MH, Henderson L, Tyagi R, Bachani M, Steiner J, Campanac E, Hoffman DA, von Geldern G, Johnson K, Maric D, Morris HD, Lentz M, Pak K, Mammen A, Ostrow L, Rothstein J, Nath A (2015) Human endogenous retrovirus-K contributes to motor neuron disease. Sci Transl Med 7:307ra153.
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2-∆∆Ct method. Methods 25:402–408.
MacGowan DJ, Scelsa SN, Imperato TE, Liu KN, Baron P, Polsky B (2007) A controlled study of reverse transcriptase in serum and CSF of HIV-negative patients with ALS. Neurology 68:1944–1946.
MacGowan DJ, Scelsa SN, Waldron M (2001) An ALS-like syndrome with new HIV infection and complete response to antiretroviral therapy. Neurology 57:1094–1097.
Matsuzaki T, Nakagawa M, Nagai M et al (2000) HTLV-Iassociated myelopathy (HAM)/tropical spastic paraparesis (TSP) with amyotrophic lateral sclerosis-like manifestations. J Neurovirol 6:544–548.
Mayer J, Harz C, Sanchez L, Pereira GC, Maldener E, Heras SR, Ostrow LW, Ravits J, Batra R, Meese E, García-Pérez JL, Goodier JL (2018) Transcriptional profiling of HERV-K(HML-2) in amyotrophic lateral sclerosis and potential implications for expression of HML-2 proteins. Mol Neurodegener 13:39.
McCormick AL, Brown RH, Cudkowicz ME, Al-Chalabi A, Garson JA (2008) Quantification of reverse transcriptase in ALS and elimination of a novel retroviral candidate. Neurology 70:278–283.
Nath A, National Institute of Neurological Disorders and Stroke (2015) HERV-K suppression using antiretroviral therapy in volunteers with amyotrophic lateral sclerosis (ALS). ClinicalTrials.gov identifier: NCT02437110 https://clinicaltrials.gov/ct2/show/NCT02437110.
Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP (2004) Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper--excel-based tool using pair-wise correlations. Biotechnol Lett 26:509–515.
Silver N, Best S, Jiang J, Thein SL (2006) Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol Biol 7:33.
Steele AJ, Al-Chalabi A, Ferrante K, Cudkowicz ME, Brown RH Jr, Garson JA (2005) Detection of serum reverse transcriptase activity in patients with ALS and unaffected blood relatives. Neurology 64:454–458.
Subramanian RP, Wildschutte JH, Russo C, Coffin JM (2011) Identification, characterization, and comparative genomic distribution of the HERV-K (HML-2) group of human endogenous retroviruses. Retrovirology 8:90.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3:research0034.1.
Verma A, Berger JR (2006) ALS syndrome in patients with HIV-1 infection. J Neurol Sci 240:59–64.
Xie F, Xiao P, Chen D, Xu L, Zhang B (2012) miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol Biol 80:75–84.
We thank Dr. Ashley Jones for technical assistance.
This study was supported by a generous grant from The Amyotrophic Lateral Sclerosis Association, Washington, DC 20005, USA. This is in part an EU Joint Programme - Neurodegenerative Disease Research (JPND) project. The project is supported through the following funding organisations under the aegis of JPND - www.jpnd.eu (United Kingdom, Medical Research Council (MR/L501529/1; MR/R024804/1) and Economic and Social Research Council (ES/L008238/1)) and through the Motor Neurone Disease Association. This study represents independent research part funded by the National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London. The work leading up to this publication was funded by the European Community’s Horizon 2020 Programme (H2020-PHC-2014-two-stage; grant agreement number 633413). The funding bodies had no part in the design of the study and collection, analysis, and interpretation of data or in writing the manuscript.
Availability of data and materials
The datasets used and/or analysed during the current study are available from corresponding authors on reasonable request.
Ethics approval and consent to participate
Full ethical approval for this study was obtained from the relevant University Research Ethics Committees. Consent was obtained for autopsy in accordance with the regulations of the Medical Research Council (MRC) Neurodegenerative Disease Brain Bank Network, Institute of Psychiatry, Psychology and Neuroscience, Kings College London, UK.
Consent for publication
AAC reports consultancies for Biogen Idec, Mitsubishi-Tanabe Pharma, Cytokinetics Inc., Treeway, Chronos Therapeutics Inc., GSK and OrionPharma in the past 5 years, and was Chief Investigator for clinical trials from Cytokinetics Inc. and OrionPharma. AAC is a collaborator in a Phase 2 clinical trial of Triumeq, an antiretroviral drug, in ALS. None of the other co-authors report any competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figures S1.–Figure S11. and Tables S1–Table S3. (PDF 874 kb)