Next Article in Journal
Impact of Heavy Metals on Cold Acclimation of Salix viminalis Roots
Next Article in Special Issue
Beyond Amyloid: A Machine Learning-Driven Approach Reveals Properties of Potent GSK-3β Inhibitors Targeting Neurofibrillary Tangles
Previous Article in Journal
E3 Ubiquitin Ligase Uhrf2 Knockout Reveals a Critical Role in Social Behavior and Synaptic Plasticity in the Hippocampus
Previous Article in Special Issue
History of Parkinson’s Disease-Associated Gene, Parkin: Research over a Quarter Century in Quest of Finding the Physiological Substrate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Insights into Endogenous Retrovirus-K Transcripts in Amyotrophic Lateral Sclerosis

by
Laura Moreno-Martinez
1,2,3,4,†,
Sofía Macías-Redondo
5,†,‡,
Mark Strunk
5,
María Isabel Guillén-Antonini
3,
Christian Lunetta
6,7,8,
Claudia Tarlarini
8,
Silvana Penco
8,
Ana Cristina Calvo
1,2,3,4,
Rosario Osta
1,2,3,4 and
Jon Schoorlemmer
3,5,9,*
1
Laboratory of Genetics and Biochemistry (LAGENBIO), Faculty of Veterinary, University of Zaragoza, Miguel Servet 177, 50013 Zaragoza, Spain
2
Centro de Investigación Biomédica en Red en Enfermedades Neurodegenerativas, Instituto de Salud Carlos III (CIBER-CIBERNED-ISCIII), 28029 Madrid, Spain
3
Instituto de Investigación Sanitaria Aragón (IIS Aragón), 50009 Zaragoza, Spain
4
Instituto Agroalimentario de Aragón (IA2), University of Zaragoza-CITA, C/Miguel, Servet 177, 50013 Zaragoza, Spain
5
Instituto Aragonés de Ciencias de la Salud (IACS), Centro de Investigación Biomédica de Aragón (CIBA), 50009 Zaragoza, Spain
6
NEMO (NEuroMuscular Omnicentre) Clinical Center, Fondazione Serena Onlus, 20162 Milan, Italy
7
Neurorehabilitation Department of Milano Institute, Istituti Clinici Scientifici Maugeri IRCCS, 20138 Milan, Italy
8
Medical Genetics Unit, Department of Laboratory Medicine, ASST Grande Ospedale Metropolitano Niguarda, 20162 Milan, Italy
9
ARAID Foundation, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: Certest Biotec, 50840 Zaragoza, Spain.
Int. J. Mol. Sci. 2024, 25(3), 1549; https://doi.org/10.3390/ijms25031549
Submission received: 20 December 2023 / Revised: 12 January 2024 / Accepted: 22 January 2024 / Published: 26 January 2024

Abstract

:
Retroviral reverse transcriptase activity and the increased expression of human endogenous retroviruses (HERVs) are associated with amyotrophic lateral sclerosis (ALS). We were interested in confirming HERVK overexpression in the ALS brain, its use as an accessory diagnostic marker for ALS, and its potential interplay with neuroinflammation. Using qPCR to analyze HERVK expression in peripheral blood mononuclear cells (PBMCs) and in postmortem brain samples from ALS patients, no significant differences were observed between patients and control subjects. By contrast, we report alterations in the expression patterns of specific HERVK copies, especially in the brainstem. Out of 27 HERVK copies sampled, the relative expression of 17 loci was >1.2-fold changed in samples from ALS patients. In particular, the relative expression of two HERVK copies (Chr3-3 and Chr3-5) was significantly different in brainstem samples from ALS patients compared with controls. Further qPCR analysis of inflammation markers in brain samples revealed a significant increase in NLRP3 levels, while TNFA, IL6, and GZMB showed slight decreases. We cannot confirm global HERVK overexpression in ALS, but we can report the ALS-specific overexpression of selected HERVK copies in the ALS brain. Our data are compatible with the requirement for better patient stratification and support the potential importance of particular HERVK copies in ALS.

1. Introduction

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by progressive degeneration and the functional loss of the upper and lower motor neurons (MNs) in the brainstem (including the medulla oblongata), spinal cord, and motor cortex [1]. As no cure is available [2], ALS leads to respiratory muscle dysfunction, paralysis, and death, usually within 3–5 years of diagnosis. While the disease is selective for motor neurons, multiple cell types, including muscle cells, astrocytes, oligodendrocytes, and microglia, are involved in its pathogenesis [3]. In 90–95% of cases, ALS is a sporadic disorder (sALS) with an unknown cause [2]. A small percentage of cases, 5–10%, are considered familial (fALS). Mutations in over 20 genes have been identified in fALS [4]. While apparently disparate, the products of many of these genes are associated with the U1 snRNP and RNA polymerase II (RNAP II) machinery [5], suggesting ALS defects in MN survival and RNA processing pathways. The exact etiology of ALS is still unclear, and a variety of cellular pathways are dysregulated in this disease [6].
ALS is a heterogeneous disease [7,8], and it is as yet unknown whether different phenotypes of ALS represent a clinical continuum or are associated with different pathophysiological mechanisms [9]. Altered inflammatory (auto)immune responses and activities are believed to contribute to the pathogenesis of ALS (reviewed in [10]). The systemic inflammation of both the peripheral and central nervous systems has been described in ALS patients, including the presence of proinflammatory monocytes [11]. In particular, increased mRNA expression of both interleukin 6 (IL6) and NLR family pyrin domain-containing 3 (NLRP3) in monocytes has been associated with progressing ALS [11]. In human ALS tissue, increased levels of several inflammasome components [12], including NLRP3, have been detected in astrocytes [13], and a meta-analysis showed increased blood levels of interleukin 1 beta (IL1B) in ALS [14]. Granzyme A (GZMA) and granzyme B (GZMB) are cell death-inducing enzymes that are released from granules in T-lymphocytes and natural killer cells in inflammatory conditions [15]. Levels of GZMA and GZMB are significantly increased in the sera of ALS patients compared with control groups [16].
ERVs (endogenous retroviruses) are the remnants of ancestral retroviral insertions into the germline. The resulting permanent retroviral integrations, termed proviruses, are transmitted vertically to offspring (reviewed in [17]). Human ERVs (HERVs) are part of the LTR retrotransposon group, exhibiting the typical structure of a provirus: GAG, POL, and ENV genes flanked by two LTRs. Since the initial integration events, HERV elements have become widely distributed and occupy approximately 8% of the human genome [18]. Although many copies express RNA and proteins, no infectious HERVs have been observed to date [19]. In humans, the HERVK (human endogenous retrovirus-K, HML-2 subfamily) group is considered the most recently acquired and is composed of betaretrovirus-like ERV [20]. This group comprises a range of similar but clearly distinguishable elements in the form of over one hundred proviruses and thousands of nearly full-length solitary LTRs. The anomalous expression of LINE-1, endogenous retroviruses (ERVs), and other repetitive elements has been mentioned in the context of degenerative [21,22] and autoimmune diseases, as well as in cancer [23]. Transcripts, proteins, and even retrovirus-like particles originating from HERVK proviruses have been associated with a wide variety of diseases (reviewed in [20,22]).
HERVK overexpression has been reported in serum [24] and the brain and cerebrospinal fluids (CSFs) of ALS patients [25,26], and reverse transcriptase activity can be found in the blood of patients [27,28]. Similarly, an analysis of eleven postmortem brain tissue samples taken from sALS patients yielded 2–3-fold higher levels of HERVK mRNA encoding GAG, POL, and ENV genes in the postmortem cortexes of patients [25]. Over the past decade, multiple studies have shown that the pathogenicity of ALS is significantly impacted by the participation of ERVs [29]. For example, the phenotype of a transgenic mouse model expressing HERVK ENV is reminiscent of ALS [25], and extracellularly released HERVK ENV is toxic to neuronal cells [26]. It was recently shown that HERVK ENV protein is released extracellularly in 11 out of 14 CSF samples taken from ALS patients and causes neurotoxicity both in vivo and according to in vitro neurotoxicity assays via a mechanism that involves receptor binding [26]. HERVK has also been considered a molecular target for ALS therapy in separate clinical trials [30,31]. While these results further support a role for HERVK in ALS pathogenesis in selected cases, the extent of overexpression in postmortem samples remains controversial since additional studies have been unable to confirm the association between elevated cortical HERVK mRNA levels and ALS [26,27,32].

2. Results

2.1. ALS-Associated HERVK Expression in PBMC

We first analyzed the expression of HERVK-encoded POL and ENV transcripts in the peripheral blood mononuclear cells (PBMCs) of ALS patients and compared levels with healthy subjects (non-diseased individuals, NDIs) or subjects suffering from other neuropathies (ONPs). Of the ten ALS patients analyzed (average age, 67 ± 11.16, six males and four females), eight represent cases of familial ALS (fALS) given the presence of one of the four mutations in the ALS susceptibility genes SOD1/TDP43/FUS/C9orf72. Data on the study subjects are detailed in Table 1. Both HERVK ENV and POL expression were readily observed in ALS patients (Figure 1A,B). Relative POL expression was not increased in the ALS group, as relative expression was 0.88 ± 0.21–0.27 compared with 1.0 ± 0.39–0.64 in the control group and 1.01 ± 0.13–0.15 in the ONP group. Relative ENV expression was also not increased in the ALS group, as relative expression was 0.93 ± 0.24–0.33 compared with 1.0 ± 0.33–0.51 in the control group and 1.00 ± 0.08 in the ONP group.
There was a moderate correlation (R2 > 0.50) between the expression of each of these genes in both the ALS samples (Figure 1C) and the combined non-ALS samples (Figure 1D), confirming the reliable detection of expression of the retroviral elements.

2.2. HERVK Expression in ALS Brain

Upregulated HERVK expression was previously described in the frozen brain tissue of ALS patients in [25]. In the absence of increased HERVK levels in our PBMC samples of ALS patients, we repeated our analysis based on RNA obtained from postmortem brain samples (Table 2). We analyzed the presence of HERVK ENV transcripts in brainstem, which represented the region affected in progressive patients: eight NDIs and eight ALS individuals. In addition, disease in this region accounts for most of the worst symptoms of ALS. Simultaneously, we analyzed the cortex samples of six NDI and ten ALS patients to enable a direct comparison with published data [25]. ALS patients were compared with NDIs. Similar to the PBMC results, HERVK ENV and POL expression were readily observed in the ALS patients (Figure 1E–H). No increase in ENV or POL transcripts was observed in ALS compared with the NDIs, neither in the cortex nor in the brainstem.
The presence of contaminating genomic sequences in our RNA preparations was evaluated by comparing cDNA templates treated either with or without reverse transcriptase. The PCR amplification of these templates with specific primers yielded the results depicted in Figure S1. For both genes analyzed, the levels detected were at least 10-fold higher after an RT step, and most samples were much higher, up to 104-fold. These results indicate that the levels of contaminating DNA were low compared with cDNA levels. We conclude that the expression levels we measured were not influenced by contaminating DNA.

2.3. Analysis of Transcribed HERVK Loci

We performed next-generation sequencing (NGS) analysis to identify individual HERVK copies with altered expression in ALS. NGS analysis was performed on amplification products obtained using the specific PCR assay for HERVK ENV, which had been previously used for expression analysis (described in Figure 1 and Table S5). Random-primed cDNAs from the PBMCs of ALS patients (n = 10) and controls (n = 10) were used as templates. Reads (70,694 ± 24,812 per sample; 25,286–136,704; Table S1) were mapped to the human genome according to stringent criteria described in the Section 4. Once assigned to unique genomic locations, reads corresponding to the 27 HERVK ENV loci identified in the PCR simulations in silico were extracted (Table S2). As expected, 99.85–99.99% of mapped reads corresponded to 27 loci identified in PCR simulations in silico (Table S3). Relative frequencies were calculated as the number of reads mapping to an individual HERVK element relative to the total number of mapped reads.
The resulting data from the PBMC samples showed that the majority of reads mapped to a limited number of loci (>5% of total reads/locus), specifically to HERVK copies located on chromosomes 1 (chr1_3); 3 (chr3_2); and, to a minor extent, chromosomes 1 (chr1_3), 5 (chr5_2), 6 (chr6), 7 (chr7_1 and chr7_2), 8 (chr8_2), 12 (chr12), 19 (chr19_2), and 22 (chr22) (Figure 2A). Lower numbers of reads mapped to the remaining 16 loci, with relative frequencies ranging from 0.01 to 1.66%. We found no significant differences between ALS patients and controls in the relative frequency of reads mapping to any of these loci (Figure 2A,B). The average fold change was 1.33 ± 0.22 (Figure 3A). Assuming that the relative frequency after mapping is a measure of relative overexpression, we were unable to find copies of HERVK that were differentially expressed in the PBMCs of ALS patients.
Although we were unable to detect increased overall expression levels of HERVK in ALS-patient-derived postmortem brain samples (Figure 1), we wondered whether specific copies of HERVK might be differentially expressed in these samples. We repeated the analysis described above for PBMCs, including samples obtained from both cortex tissue and the brainstem. The resulting data show that the origin of the reads is different in the brain versus PBMCs (Figure 2A,B). Again, most brain-derived reads mapped to a limited number of loci, in this case, to HERVK copies 3–2 (>48%) and 3–4 (7.5%). In contrast to the PBMC data, we found differences between ALS patients and controls in the relative frequency of reads mapping to several loci to which lesser numbers of reads were mapped. The average relative fold changes were 0.40 ± 1.7 and 0.59 ± 3.58 in the cortex and brainstem, respectively (Figure 3A). To better survey the differences associated with ALS, the relative fold change was calculated, excluding the 10 loci that are unchanged in patients (Table S3). The relative fold changes among the remaining 17 loci were 0.60 ± 1.2 (p = 0.002) and 0.97 ± 4.1 (p = 0.05) in the cortex and brainstem, respectively (Figure 3B).
Reads mapping to a subset of HERVK loci, including copy_chr1-1 (not shown), copy_chr3-3, copy_chr3-5, and copy_chr16-1, differed between ALS patients and NDIs (Figure 3C). The relative frequency (ranging from 0.3 to 2.95%) of reads mapping to these loci showed an increase (copies chr3-5 and chr16-1) in ALS patients compared with controls (Figure 3C) and a greater than two-fold decrease in ALS brainstems compared with NDIs (copy chr3-3). When expressed as log2 fold changes, several copies, including chr2, chr5-1, and chr16-1, displayed differential expression (increased in ALS) to varying extents (−2.42-, −1.87-, and −0.62-fold, respectively). When corrected for multi-locus analysis, significant differential expression (p < 0.05; FDR < 0.05) in two loci was evident (Figure 3C): a decrease in the expression of HERVK_copy_chr 3.3 (1.8-fold) and an increase in the expression of HERVK_copy_chr 3.5 (−3.4-fold). These data indicate that the relative frequency of specific HERVK copies is altered in the brains of ALS patients, even in the absence of overall HERVK overexpression (Figure 2B). Such altered frequency appears to be derived from less abundantly transcribed but potentially ALS-related HERVK loci.

2.4. Inflammation Markers in ALS Brains

Peripheral inflammation is associated with ALS [8], as the induction of IL6 and tumor necrosis factor alpha (TNFα) markers, among others, has been reported in the blood samples of patients [14]. A potential infiltration of activated immune cells and the local activation of microglial cells were analyzed by assessing IL6, TNFA, and GZMB expression levels using qPCR in postmortem tissue from ALS patients. The levels of all three inflammation markers trended lower in ALS patients compared with NDIs (Figure S2). While slight decreases were found for IL6 and TNFA in the brainstem, in the cortex, we observed significant decreases in both GZMB and TNFA levels in ALS patients. TNFA levels were two-fold decreased in ALS patients compared with NDIs.
We also analyzed the expression of NLRP3, an inflammasome component (Figure 4). NLRP3 mRNA levels were increased in postmortem brain tissue from ALS patients compared with NDIs (p = 0.048 for the brainstem; p < 0.001 for all brain samples combined).

3. Discussion

Higher levels of HERVK mRNA encoding GAG, POL, and ENV genes in postmortem brain tissue samples taken from ALS patients have been previously reported [25]. In the hope of improving the detection and diagnosis of ALS, we searched for increased expressions of HERVK in patients with ALS and other neuropathies. We were unable to detect increased levels of the HERVK transcript in PBMCs taken from ALS patients. Although immune alterations have been described in the blood of ALS patients [14], PBMCs are not primarily affected by this disease, foreshadowing the negative results we obtained. Increased levels of HERVK transcripts have been reported in serum [24], but it is unknown whether these reflect circulating DNA or ribonucleotides contained in microvesicles. However, the absence of increased levels of HERVK transcripts in ALS brains was not anticipated. Our results confirm more recent studies in which no difference in HERVK GAG, POL, or ENV levels was observed in brain and spinal cord samples from ALS patients and controls [33,34]. Our failure to detect increased HERVK expression in postmortem ALS brains may reflect the small sample sizes employed or the different cellular composition of the clinical samples analyzed. However, a similar conclusion was reached by other studies analyzing larger numbers of samples [33,34]. Therefore, the particular HERVK elements, or even the HERV family overexpressed, may differ between groups and patients or depend on geographical location. Given the heterogeneity found in ALS, it is as yet unknown whether different phenotypes of ALS represent a clinical continuum or are associated with different pathological mechanisms [35]. In this context, HERVK overexpression may well be restricted to a subgroup of patients, as already mentioned in the initial study reporting increased expression [25]. Such stratification of cases was recently reported for HERVW expression in psychotic spectrum patients [36] and the induction of HERVW ENV expression by SARS-CoV-2 [37]. The disease-associated expression of HERV may, therefore, be more common than previously reported.
Out of three manuscripts closely related to our study, Garson et al. [34] only addressed overall HERVK levels, as opposed to the identification of specific copies. After careful assessment of PCR conditions, no increases in HERVK expression were reported in that study. Two previous studies examined copy-specific HERVK expressions in ALS. These studies analyzed cortex samples [25] and/or spinal cord tissue samples [33] and did not analyze chr3-5 and chr16-1 copies, as they were based on older versions of the genome database.
Using a different approach to ours (Illumina mRNA-Seq sequencing after ribosomal RNA depletion and posterior extraction and counting of HERVK-homologous reads), slight, statistically insignificant differences (<1.15-fold) between control and ALS samples were reported in HERVK copies 7q34 (c7_C), 10p14 (c10_A), and 3q21.2 (K(I)) [25]. The latter corresponds to the copy we refer to as chr3-4 (Table S2). The loci we found differentially expressed in the brainstem were not identified in this previous study [25]. Another recent report describes the transcriptional profiling of HML-2 based on cloning PCR products [33,34]. With respect to HERVK ENV, the results were derived from only a few reads, i.e., 91 for ALS cases and 53 controls. While the study reports that “significant differences in HML-2 loci transcriptional activities were not seen when comparing ALS and controls”, ENV expression from a copy named chr3q13.2_K-3 trends toward a reduction in ALS samples [33]. This is the same copy we call HERVK_copy_chr3-3 (Table S2), whose expression was significantly reduced in brainstem samples from ALS patients (Figure 3C). It is of future interest to repeat these findings in additional samples. Future experiments using our pipeline may be improved by the use of a spiked control to improve normalization and the use of a high-fidelity Taq. A subsequent step will and need be to determine the significance of the altered expression of this locus. As opposed to the increased overall expression of HERVK in ALS clinical samples, we identified specific copies of HERVK that are relatively more abundantly expressed in brain tissue but not in the PBMC samples of patients. HERVK ENV protein was reportedly present in the CSFs of most ALS patients [26], as well as in the sera of ALS individuals [38]. However, the grade of CSF toxicity was not a direct function of the ENV protein concentrations measured. While this may also result from differences in post-translational modifications of these particular ENV proteins and/or the presence of non-ENV components in the CSF, it may very well be related to the specific copy of HERVK expressed.
Our results are reminiscent of a previous analysis of the RNA-seq data of ALS vs. control brain tissue samples [39]. While no global differences were reported for LTR retrotransposon subfamilies, specific copies showed altered expression in ALS. To further evaluate the overexpression of the ALS-related HERVK copies we identified, specific qPCR assays should be designed, validated, and applied to patient samples. In line with this, a recent study that examined the locus-specific expression of HERVs found that individuals with ALS consistently exhibited the upregulation of one specific HERV locus, HML6_3p21.31c, in both the motor cortex and cerebellum when compared with control subjects [40]. In addition, a recent study on the human retrotransposon called long interspersed element-1 (L1) focused on its locus-specific expression in ALS, revealing an overall decrease in intact L1 expression in two brain regions of affected individuals, with distinct clustering based on expression patterns, highlighting the importance of understanding the regulation of specific L1 subsets in these tissues and suggesting avenues for further research [41].
The HERVK copies identified in this study reside in chromosomes 1, 3, and 16, which also harbor the TARDBP (TAR DNA-binding protein or TDP-43), CHMP2B, and FUS genes, which are frequently mutated in familial cases of ALS [42]. As the distance between HERVK elements and ALS susceptibility genes is large, the significance of this localization is presently unclear. Several gene mutations (including TARDBP and FUS) associated with ALS converge on RNA-binding (or -processing) proteins and RNA metabolism [6,39]. As a result, the nucleocytoplasmic shuttling of RNA–protein complexes is altered, resulting in the formation of RNA granules, the loss of RNA binding (including those derived from transposable elements [43]), and defects in RNA splicing [42]. Overexpressed HERVK copies may produce toxic proteins relevant to ALS progression [33]. However, RNA detected in the cytosol or in endosomes is known to activate immune pathways through intracellular receptors [44]. As opposed to increased levels of HERVK RNAs in ALS, pathological pathways in ALS may alter the processing of HERVK RNAs and make them more available to these RNA-detecting pathways [45]. Such a process may depend on the sequence or structure of specific copies. A contribution of HERVK transcripts to ALS can now be envisaged in different ways. HERVK-encoded ERV protein may contribute to ALS-related symptoms such as the degeneration of motor neurons and motor dysfunction, as described for forced (over)expression in mouse brains [25]. In addition, we suggest that the ALS-specific HERVK RNAs identified in this work may contribute to the disease through the activation of innate immune pathways.
Altered inflammatory (auto)immune responses and activities also contribute to the pathogenesis of ALS [10]. Surprisingly, our results do not show increased levels of TNFA and IL6 transcripts in the cortexes or brainstems of postmortem ALS brains. The observed decrease in levels of inflammation markers in ALS brains with respect to control brains could suggest that neuronal loss in the terminal stage of the disease dampens inflammation. The observation of apparently reduced inflammation is based on a very limited set of markers and awaits confirmation via additional methods, such as histological analysis in combination with specific antibody staining in an extended set of samples.
In mouse models of ALS (e.g., SOD1), DAMPs released from injured motor neurons induce microglia to acquire an M1 phenotype, which is associated with the enhanced secretion of the proinflammatory cytokines TNFα, IL1, and IL6 [10,46]. We report increased NLRP3 expression in ALS brains (especially in the cortex; see Figure 4A). This is the first study that describes NLRP3 levels in postmortem brain tissue. The NLRP3 overexpression observed confirms that inflammation in the form of NLRP3 inflammasome activation may be associated with ALS, in accordance with the activation of NLRP3 inflammasome in SOD1G93A mice [47]. At present, we do not know whether other components of this inflammasome or IL1/IL18 levels are increased in the ALS brain samples analyzed.
The NLRP3 inflammasome is activated in response to a variety of agents, including bacterial RNA and toxins, K+ efflux, extracellular ATP, uric acid, and mitochondrial reactive oxygen species [48]. In addition, mutant forms of TDP43 [49], which are genetically associated with ALS (and wild-type forms, to a lesser extent), can activate microglia through the NLRP3 inflammasome [50,51]. As human immunodeficiency virus type-1 single-stranded RNA activates the NLRP3 inflammasome in human microglia [52], we suggest a similar response may occur in response to HERVK-derived RNA. The depletion of TDP43 from the nucleus and aggregation in the neuronal cytoplasm [49] is also associated with the deregulated expression of transposable elements in several model systems [53]. NLRP3 inflammasome activation may depend on dual signals from mutant TDP43 and HERVK derived RNA. Altogether, our results are compatible with the idea that microglial inflammation in the form of inflammasome activation is a common feature of ALS.

4. Materials and Methods

4.1. Participants, Samples, and Approvals

Blood samples from patients (Table 1) were obtained under informed written consent prior to inclusion in this study, as approved by the institutional review boards of the institutional ethical committees of Niguarda Ca’ Granda Hospital, Milan, Italy. This study was conducted according to Declaration of Helsinki principles and according to Directive 2004/23/EC of the European Parliament and of the Council.
Brain samples and data from patients included in this study (patients and non-diseased controls) (Table 2) were collected, processed, and provided by Banco de Tejidos CIEN, Fundación CIEN, Instituto de Salud Carlos III, Madrid, and by Biobanco en Red de la Región de Murcia (BIOBANC-MUR), Murcia, Spain. The latter is integrated into the Spanish Biobanks Network Platform (www.redbiobancos.es (accessed on 4 October 2019)) and registered on the Registro Nacional de Biobancos with registration number B.0000859. Patients and unrelated control individuals were recruited with informed written consent, and samples were processed following standard operating procedures as approved by the local institutional review board (El Comité de Ética de la Investigación de la Comunidad de Aragón (CEICA), approval reference CEICA PI17/0025, Zaragoza, Aragón, Spain) and the Comité Científico del banco de Tejidos de la Fundación CIEN (approval reference CCS17003).
All samples were analyzed by investigators who were blinded to both the clinical conditions and identities of the patients.

4.2. Experimental Design

This study included a random sample of PBMCs extracted from 10 familial ALS patients selected from a previously studied group [54], matched for age and gender as far as possible (Table 1). The “other neuropathies group” (ONP group) included diseases such as myotonic dystrophy type 1 and 2, facioscapulohumeral muscular dystrophy, Becker’s muscular dystrophy, and mild myotonia. For the postmortem brain samples, a total of 11 ALS patients and 7 non-diseased controls (NDIs) were included, and samples were matched according to the anatomical region of the brain studied (Table 2). Sample sizes were based on the availability of tissues.

4.3. RNA Extraction, Preparation of cDNA, and qPCR

RNA extraction and purification from blood samples were described previously in [54]. RNA quality measurements for these samples are not available. For PBMC samples, random-primed cDNA was obtained from 1 μg of total extracted RNA using the High-Capacity cDNA RT kit (Applied Biosystems, Waltham, MA, USA) [54]. RNA extraction and purification from frozen brain samples were carried out as described previously for placentas [55]. No statistical differences in mean RNA quality indicator (RQI) values were encountered between ALS patients and NDI groups (Figure S4). Random-primed cDNA was prepared from 1 µg of these RNA samples using the ThermoScript™ RT-PCR System (11146-016 Thermo Fisher, Waltham, MA, USA), as described in [55]. Using serial dilutions of cDNA produced from samples with low RQI values, we established that qPCR efficiencies (Figure S3) complied with linearity and efficiency thresholds established in [56].
Gene expression of both HERVK and inflammation biomarkers was analyzed via reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) on a ViiA™ 7 System (Applied Biosystems, Waltham, MA, USA) under standard thermal cycling conditions. All technical details concerning qPCR are included in a document (Table S4) according to MIQE guidelines. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and ribosomal protein L19 (RPL19) were used as reference genes. The sequence of the primers used to amplify HERVK sequences has been described before in [25] and is shown in Table S5. Commercial Taqman assays (Thermo Fisher, Waltham, MA, USA) were employed to measure RPL19 (Hs02338565_gH), IL6 (Hs00985639_m1), TNFA (Hs01113624_g1), and NLRP3 (Hs00918082_m1). All other primers are listed in Table S5. Analyses were performed with a SYBR Green qPCR, using the SYBR Green Master Mix (SYBR Premix Ex Taq II (Tli RNase H Plus), RR820A, Takara, Kusatsu-shi, Japan) according to the manufacturer’s protocol. For TaqMan qPCR, the Premix Ex TaqTM (Probe qPCR) Master Mix (RR390A, Takara, Kusatsu-shi, Japan) was used according to the manufacturer’s protocol. Amplification efficiency between 90% and 110% was assessed in standard curves for all primer sets (Figure S3). All reactions were carried out in triplicate, and only measurements with a standard deviation < 0.2 were considered.
The presence of contaminating genomic sequences in our RNA preparations was evaluated by comparing amplification from templates treated in parallel with or without reverse transcriptase in standard RT buffer. To compare the expression between groups, we performed the 2−ΔΔCt method, as previously described [56].

4.4. Next-Generation Sequencing (NGS)

cDNA obtained from ALS patients and NDIs was amplified employing specific HERVK primers [25] (Table S5) using 0.5 U Taq DNA Polymerase (D1806, Sigma-Aldrich, St. Louis, MO, USA). NGS was carried out at the Sequencing and Functional Genomics Core Facility (Servicio Científico Técnico de Secuenciación y Genómica Funcional) of the Aragon Health Sciences Institute (IACS, Zaragoza). Library preparation and sequencing were carried out using kits approved for IonTorrent technology, following procedures recommended by the relevant manufacturers. PCR products were end-repaired and purified with Agencourt AMPure XP paramagnetic beads (Beckman Coulter, Brea, CA, USA), and amplicons were quantified using the Qubit 3 fluorometer (Qubit BR dsDNA Assay, Thermo Fisher, Waltham, MA, USA). Adapters and barcodes (Ion Xpress Barcode Adapters, Thermo Fisher, Waltham, MA, USA) were added via ligation, followed by product purification with AMPure XP beads. Libraries were prepared separately for each sample using the Kapa Lib Prep Kit for Ion Torrent (Roche, Basel, Switzerland). Libraries were quantified via qPCR (Ion Library TaqMan™ Quantitation Kit, Thermo Fisher, Waltham, MA, USA), and equimolar amounts of the samples were pooled before template preparation. Emulsion PCR, template enrichment, and chip loading were carried out manually using the Ion 520&530 Kit-OT2 (Thermo Fisher, Waltham, MA, USA), and samples were sequenced on the Ion Torrent S5XL platform using an Ion 530 chip.
Quality-filtered reads (Table S1) (exported in fastq format) were trimmed to remove primer sequences and subsequently mapped without soft clipping to the human genome (GRCh37/hg19) using TMAP version 5.2.22, applying the setting “random among best hits” for equal best hits. Only reads > 135 bp were analyzed, using (penalty) scoring as follows: match, +1; mismatch, −10; and gap, −10. Accumulated mismatch percentages were <0.75%. Subsequently, reads mapping to selected genomic HERVK loci were extracted and counted. Relevant loci were selected as those HERVK loci identified by BiSearch in a virtual PCR using the HERVK ENV primers employed throughout this study (Tables S2 and S3). As BiSearch uses the GRch 38.92 version of the Homo sapiens genome, locations were transformed into the GRCh37/hg19 version before extraction using the Ensemble assembly converter (https://www.ensembl.org/Homo_sapiens/Tools/AssemblyConverter?db=core (accessed on 11 July 2018)).
The counting reads per locus were recalculated as percentages of reads relative to the total number of reads in the sample (Table S1). Fold change was calculated as the median percentage of reads in a group divided by the percentage in the relevant control group. Relative fold change (the difference between the experimental value and a control value) was defined as the deviation from a 1-fold change (x = fold change; if x > 1, x − 1; if x < 1, (1/x) − 1); if x = 0, relative fold change was also set at 0).
Data were further analyzed using the DESeq2 package [57], which tests for differential expression based on logarithmic fold changes. Separate raw gene-count matrixes were obtained for NDIs and for the ALS samples. The pipeline calculates log2 fold changes and provides the probability values of the loci, correcting p-values for multiple testing (false discovery rate < 0.05). A locus with less than 1 raw count in more than 50% of samples (i.e., HERV.K_copy_chr19.1) was excluded from downstream analysis.

4.5. Statistical Analysis

The SPSS software was used for all analyses (IBM Corp., Armonk, NY, USA, released 2013; IBM SPSS Statistics for Windows, Version 15.0). Normality was assessed with the Shapiro–Wilk test (when n < 50) or Kolmogorov–Smirnov test (when n > 50). For the assessment of the statistical significance of differences, a specific test was performed depending on normality and the number of groups of samples compared. Normally distributed data: Student’s t-test (2 groups) or ANOVA test (>2 groups). Non-normally distributed data: Mann–Whitney U test (2 groups) or Kruskal–Wallis test (>2 groups) were used when appropriate. Graphs were generated using SPSS, GraphPad Prism (version 8), and Microsoft Office Excel 2010.

5. Conclusions

Our observations fail to support the hypothesis that elevated HERVK expression in the brain is generally associated with sporadic ALS. However, we identified specific copies of HERVK, whose expression in the brain may be associated with ALS. As our study is based on a relatively small sample size, the relationship with ALS pathology and potential implications for ALS drug targeting remain to be established. Certainly, our results endorse continued research in this area with respect to patient stratification, the presence of HERVK-derived RNA and proteins in ALS patients and their contribution to this disease, and the potential therapeutic benefits of suppressing HERVK activity.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25031549/s1.

Author Contributions

A.C.C., R.O. and J.S. were responsible for the conception of the study. A.C.C., C.L., C.T. and S.P. recruited the patient cohort. S.M.-R. and M.I.G.-A. prepared RNA and cDNA from brain samples. L.M.-M. and S.M.-R. carried out gene expression analysis. J.S., S.M.-R. and M.S. generated the NGS data. A.C.C., L.M.-M., R.O. and J.S. performed the data analysis (with help from M.S.), interpreted the results, and wrote the manuscript. J.S. supervised the study. All the authors critically revised the manuscript for intellectual and scientific content and approved the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Ayudas de Investigación sobre Medicamentos Huérfanos y Enfermedades Raras 2016 (La Ayuda de Colegios de Farmacéuticos), Fundación Mehuer, Sevilla, Spain. The Osta lab was funded by the Instituto de Salud Carlos III (Grant PI21/00372) and Fondo Europeo de Desarrollo Regional (FEDER), “Una manera de hacer Europa”, from the European Union, Centro de Investigación Biomédica en Red sobre Enfermedades Neurodegenerativas (CIBERNED), Consolidated Groups from Gobierno de Aragón, Fundación FEDER. The Schoorlemmer lab is supported by grants from the Government of Aragon, co-financed by FEDER (“Una manera de hacer Europa, Construyendo Europa desde Aragón”/European Social Funds: B05_17D, B46_20R and B46_23R). S.M.-R. is supported by Ph.D. fellowship C071/2014 from the DGA (Aragón, Spain). BIOBANC-MUR is supported by the “Instituto de Salud Carlos III (proyecto PT13/0010/0018)2” from the “Instituto Murciano de Investigación Biosanitaria Virgen de la Arrixaca, IMIB, (Proyecto AD210/0057)” and by the “Consejería de Sanidad y Política Social de la Comunidad Autónoma de la Región de Murcia”.

Institutional Review Board Statement

The study protocol was approved by the local ethics committee (El Comité de Ética de la Investigación de la Comunidad de Aragón or CEICA) (CP–CI PI17/0025; 15 February 2017) and by the Comité Científico del Banco de Tejídos de la Fundación CIEN (S17003; 15 March 2017).

Informed Consent Statement

Blood samples from patients were obtained with informed written consent prior to inclusion in this study, which was conducted according to Declaration of Helsinki principles, Directive 2004/23/EC of the European Parliament and of the Council, and the institutional ethical committees of Niguarda Ca’ Granda Hospital (approval Nº 636-122015; 23 December 2015).

Data Availability Statement

The datasets generated that form the basis of the analysis presented in the current study are available from the corresponding authors upon request.

Acknowledgments

We acknowledge the use of Servicios Científico Técnicos del CIBA (IACS-Universidad de Zaragoza), specifically, the “Secuenciación y Genómica Funcional” unit. Special thanks go to Montserrat Salas Valero (IACS) for reference management, to Alberto Cebollada Solanas (Unidad de Biocomputación) for statistical analysis, and to Ana Medel-Martínez for control experiments. We are particularly grateful for the generous contribution of the patients and the collaboration of the Biobank Network of the Region of Murcia, BIOBANC-MUR, integrated into the Spanish Biobanks Network Platform and the Banco de Tejídos de la Fundación CIEN.

Conflicts of Interest

The authors declare no conflicts of interest. Neither the funders nor the Biobanks had a role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Mitchell, J.D.; Borasio, G.D. Amyotrophic Lateral Sclerosis. Lancet 2007, 369, 2031–2041. [Google Scholar] [CrossRef]
  2. Andersen, P.M.; Al-Chalabi, A. Clinical Genetics of Amyotrophic Lateral Sclerosis: What Do We Really Know? Nat. Rev. Neurol. 2011, 7, 603–615. [Google Scholar] [CrossRef] [PubMed]
  3. Puentes, F.; Malaspina, A.; van Noort, J.M.; Amor, S. Non-Neuronal Cells in ALS: Role of Glial, Immune Cells and Blood-CNS Barriers. Brain Pathol. 2016, 26, 248–257. [Google Scholar] [CrossRef] [PubMed]
  4. Corcia, P.; Couratier, P.; Blasco, H.; Andres, C.R.; Beltran, S.; Meininger, V.; Vourc’h, P. Genetics of Amyotrophic Lateral Sclerosis. Rev. Neurol. 2017, 173, 254–262. [Google Scholar] [CrossRef] [PubMed]
  5. Chi, B.; O’Connell, J.D.; Yamazaki, T.; Gangopadhyay, J.; Gygi, S.P.; Reed, R. Interactome Analyses Revealed That the U1 SnRNP Machinery Overlaps Extensively with the RNAP II Machinery and Contains Multiple ALS/SMA-Causative Proteins. Sci. Rep. 2018, 8, 8755. [Google Scholar] [CrossRef]
  6. Gao, F.-B.; Almeida, S.; Lopez-Gonzalez, R. Dysregulated Molecular Pathways in Amyotrophic Lateral Sclerosis-Frontotemporal Dementia Spectrum Disorder. EMBO J. 2017, 36, 2931–2950. [Google Scholar] [CrossRef] [PubMed]
  7. Zou, Z.-Y.; Zhou, Z.-R.; Che, C.-H.; Liu, C.-Y.; He, R.-L.; Huang, H.-P. Genetic Epidemiology of Amyotrophic Lateral Sclerosis: A Systematic Review and Meta-Analysis. J. Neurol. Neurosurg. Psychiatry 2017, 88, 540–549. [Google Scholar] [CrossRef]
  8. McCombe, P.A.; Henderson, R.D. The Role of Immune and Inflammatory Mechanisms in ALS. Curr. Mol. Med. 2011, 11, 246–254. [Google Scholar] [CrossRef]
  9. Zou, Z.-Y.; Liu, C.-Y.; Che, C.-H.; Huang, H.-P. Toward Precision Medicine in Amyotrophic Lateral Sclerosis. Ann. Transl. Med. 2016, 4, 27. [Google Scholar] [CrossRef]
  10. Liu, J.; Wang, F. Role of Neuroinflammation in Amyotrophic Lateral Sclerosis: Cellular Mechanisms and Therapeutic Implications. Front. Immunol. 2017, 8, 1005. [Google Scholar] [CrossRef]
  11. Zhao, W.; Beers, D.R.; Hooten, K.G.; Sieglaff, D.H.; Zhang, A.; Kalyana-Sundaram, S.; Traini, C.M.; Halsey, W.S.; Hughes, A.M.; Sathe, G.M.; et al. Characterization of Gene Expression Phenotype in Amyotrophic Lateral Sclerosis Monocytes. JAMA Neurol. 2017, 74, 677–685. [Google Scholar] [CrossRef]
  12. Hornung, V.; Latz, E. Critical Functions of Priming and Lysosomal Damage for NLRP3 Activation. Eur. J. Immunol. 2010, 40, 620–623. [Google Scholar] [CrossRef]
  13. Johann, S.; Heitzer, M.; Kanagaratnam, M.; Goswami, A.; Rizo, T.; Weis, J.; Troost, D.; Beyer, C. NLRP3 Inflammasome Is Expressed by Astrocytes in the SOD1 Mouse Model of ALS and in Human Sporadic ALS Patients. Glia 2015, 63, 2260–2273. [Google Scholar] [CrossRef] [PubMed]
  14. Hu, Y.; Cao, C.; Qin, X.-Y.; Yu, Y.; Yuan, J.; Zhao, Y.; Cheng, Y. Increased Peripheral Blood Inflammatory Cytokine Levels in Amyotrophic Lateral Sclerosis: A Meta-Analysis Study. Sci. Rep. 2017, 7, 9094. [Google Scholar] [CrossRef] [PubMed]
  15. Chowdhury, D.; Lieberman, J. Death by a Thousand Cuts: Granzyme Pathways of Programmed Cell Death. Annu. Rev. Immunol. 2008, 26, 389–420. [Google Scholar] [CrossRef] [PubMed]
  16. Iłżecka, J. Granzymes A and B Levels in Serum of Patients with Amyotrophic Lateral Sclerosis. Clin. Biochem. 2011, 44, 650–653. [Google Scholar] [CrossRef] [PubMed]
  17. Schoorlemmer, J.; Pérez-Palacios, R.; Climent, M.; Guallar, D.; Muniesa, P. Regulation of Mouse Retroelement MuERV-L/MERVL Expression by REX1 and Epigenetic Control of Stem Cell Potency. Front. Oncol. 2014, 4, 14. [Google Scholar] [CrossRef] [PubMed]
  18. de Koning, A.P.J.; Gu, W.; Castoe, T.A.; Batzer, M.A.; Pollock, D.D. Repetitive Elements May Comprise over Two-Thirds of the Human Genome. PLoS Genet. 2011, 7, e1002384. [Google Scholar] [CrossRef]
  19. Dewannieux, M.; Heidmann, T. Endogenous Retroviruses: Acquisition, Amplification and Taming of Genome Invaders. Curr. Opin. Virol. 2013, 3, 646–656. [Google Scholar] [CrossRef]
  20. Subramanian, R.P.; Wildschutte, J.H.; Russo, C.; Coffin, J.M. Identification, Characterization, and Comparative Genomic Distribution of the HERV-K (HML-2) Group of Human Endogenous Retroviruses. Retrovirology 2011, 8, 90. [Google Scholar] [CrossRef]
  21. Kõks, S.; Kõks, G. The Role of Human Endogenous Retroviruses (HERVs) in the Pathologies of the Nervous System. In Molecular-Genetic and Statistical Techniques for Behavioral and Neural Research; Gerlai, R.T., Ed.; Academic Press: Cambridge, MA, USA, 2018; pp. 519–533. ISBN 978-0-12-804078-2. [Google Scholar]
  22. Küry, P.; Nath, A.; Créange, A.; Dolei, A.; Marche, P.; Gold, J.; Giovannoni, G.; Hartung, H.-P.; Perron, H. Human Endogenous Retroviruses in Neurological Diseases. Trends Mol. Med. 2018, 24, 379–394. [Google Scholar] [CrossRef] [PubMed]
  23. Burns, K.H. Repetitive DNA in Disease. Science 2022, 376, 353–354. [Google Scholar] [CrossRef] [PubMed]
  24. Phan, K.; He, Y.; Fu, Y.; Dzamko, N.; Bhatia, S.; Gold, J.; Rowe, D.; Ke, Y.D.; Ittner, L.M.; Hodges, J.R.; et al. Pathological Manifestation of Human Endogenous Retrovirus K in Frontotemporal Dementia. Commun. Med. 2021, 1, 60. [Google Scholar] [CrossRef]
  25. Li, W.; Lee, M.-H.; Henderson, L.; Tyagi, R.; Bachani, M.; Steiner, J.; Campanac, E.; Hoffman, D.A.; von Geldern, G.; Johnson, K.; et al. Human Endogenous Retrovirus-K Contributes to Motor Neuron Disease. Sci. Transl. Med. 2015, 7, 307ra153. [Google Scholar] [CrossRef]
  26. Steiner, J.P.; Bachani, M.; Malik, N.; DeMarino, C.; Li, W.; Sampson, K.; Lee, M.H.; Kowalak, J.; Bhaskar, M.; Doucet-O’Hare, T.; et al. Human Endogenous Retrovirus K Envelope in Spinal Fluid of Amyotrophic Lateral Sclerosis Is Toxic. Ann. Neurol. 2022, 92, 545–561. [Google Scholar] [CrossRef] [PubMed]
  27. Steele, A.J.; Al-Chalabi, A.; Ferrante, K.; Cudkowicz, M.E.; Brown, R.H.J.; Garson, J.A. Detection of Serum Reverse Transcriptase Activity in Patients with ALS and Unaffected Blood Relatives. Neurology 2005, 64, 454–458. [Google Scholar] [CrossRef]
  28. McCormick, A.L.; Brown, R.H.J.; Cudkowicz, M.E.; Al-Chalabi, A.; Garson, J.A. Quantification of Reverse Transcriptase in ALS and Elimination of a Novel Retroviral Candidate. Neurology 2008, 70, 278–283. [Google Scholar] [CrossRef] [PubMed]
  29. Li, W.; Pandya, D.; Pasternack, N.; Garcia-Montojo, M.; Henderson, L.; Kozak, C.A.; Nath, A. Retroviral Elements in Pathophysiology and as Therapeutic Targets for Amyotrophic Lateral Sclerosis. Neurotherapeutics 2022, 19, 1085–1101. [Google Scholar] [CrossRef]
  30. Gold, J.; Rowe, D.B.; Kiernan, M.C.; Vucic, S.; Mathers, S.; van Eijk, R.P.A.; Nath, A.; Garcia Montojo, M.; Norato, G.; Santamaria, U.A.; et al. Safety and Tolerability of Triumeq in Amyotrophic Lateral Sclerosis: The Lighthouse Trial. Amyotroph. Lateral Scler. Front. Degener. 2019, 20, 595–604. [Google Scholar] [CrossRef]
  31. Garcia-Montojo, M.; Fathi, S.; Norato, G.; Smith, B.R.; Rowe, D.B.; Kiernan, M.C.; Vucic, S.; Mathers, S.; van Eijk, R.P.A.; Santamaria, U.; et al. Inhibition of HERV-K (HML-2) in Amyotrophic Lateral Sclerosis Patients on Antiretroviral Therapy. J. Neurol. Sci. 2021, 423, 117358. [Google Scholar] [CrossRef]
  32. Ishihara, T.; Koyama, A.; Hatano, Y.; Takeuchi, R.; Koike, Y.; Kato, T.; Tada, M.; Kakita, A.; Onodera, O. Endogenous Human Retrovirus-K Is Not Increased in the Affected Tissues of Japanese ALS Patients. Neurosci. Res. 2022, 178, 78–82. [Google Scholar] [CrossRef]
  33. Mayer, J.; Harz, C.; Sanchez, L.; Pereira, G.C.; Maldener, E.; Heras, S.R.; Ostrow, L.W.; Ravits, J.; Batra, R.; Meese, E.; et al. Transcriptional Profiling of HERV-K(HML-2) in Amyotrophic Lateral Sclerosis and Potential Implications for Expression of HML-2 Proteins. Mol. Neurodegener. 2018, 13, 39. [Google Scholar] [CrossRef]
  34. Garson, J.A.; Usher, L.; Al-Chalabi, A.; Huggett, J.; Day, E.F.; McCormick, A.L. 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 Neuropathol. Commun. 2019, 7, 45. [Google Scholar] [CrossRef]
  35. Saresella, M.; Piancone, F.; Tortorella, P.; Marventano, I.; Gatti, A.; Caputo, D.; Lunetta, C.; Corbo, M.; Rovaris, M.; Clerici, M. T Helper-17 Activation Dominates the Immunologic Milieu of Both Amyotrophic Lateral Sclerosis and Progressive Multiple Sclerosis. Clin. Immunol. 2013, 148, 79–88. [Google Scholar] [CrossRef]
  36. Tamouza, R.; Meyer, U.; Lucas, A.; Richard, J.R.; Nkam, I.; Pinot, A.; Djonouma, N.; Boukouaci, W.; Charvet, B.; Pierquin, J.; et al. Patients with Psychosis Spectrum Disorders Hospitalized during the COVID-19 Pandemic Unravel Overlooked SARS-CoV-2 Past Infection Clustering with HERV-W ENV Expression and Chronic Inflammation. Transl. Psychiatry 2023, 13, 272. [Google Scholar] [CrossRef] [PubMed]
  37. Charvet, B.; Brunel, J.; Pierquin, J.; Iampietro, M.; Decimo, D.; Queruel, N.; Lucas, A.; Encabo-Berzosa, M.d.M.; Arenaz, I.; Marmolejo, T.P.; et al. SARS-CoV-2 Awakens Ancient Retroviral Genes and the Expression of Proinflammatory HERV-W Envelope Protein in COVID-19 Patients. iScience 2023, 26, 106604. [Google Scholar] [CrossRef]
  38. Garcia-Montojo, M.; Simula, E.R.; Fathi, S.; McMahan, C.; Ghosal, A.; Berry, J.D.; Cudkowicz, M.; Elkahloun, A.; Johnson, K.; Norato, G.; et al. Antibody Response to HML-2 May Be Protective in Amyotrophic Lateral Sclerosis. Ann. Neurol. 2022, 92, 782–792. [Google Scholar] [CrossRef]
  39. Prudencio, M.; Gonzales, P.K.; Cook, C.N.; Gendron, T.F.; Daughrity, L.M.; Song, Y.; Ebbert, M.T.W.; van Blitterswijk, M.; Zhang, Y.-J.; Jansen-West, K.; et al. Repetitive Element Transcripts Are Elevated in the Brain of C9orf72 ALS/FTLD Patients. Hum. Mol. Genet. 2017, 26, 3421–3431. [Google Scholar] [CrossRef]
  40. Jones, A.R.; Iacoangeli, A.; Adey, B.N.; Bowles, H.; Shatunov, A.; Troakes, C.; Garson, J.A.; McCormick, A.L.; Al-Chalabi, A. A HML6 Endogenous Retrovirus on Chromosome 3 Is Upregulated in Amyotrophic Lateral Sclerosis Motor Cortex. Sci. Rep. 2021, 11, 14283. [Google Scholar] [CrossRef] [PubMed]
  41. Pfaff, A.L.; Bubb, V.J.; Quinn, J.P.; Koks, S. Locus Specific Reduction of L1 Expression in the Cortices of Individuals with Amyotrophic Lateral Sclerosis. Mol. Brain 2022, 15, 25. [Google Scholar] [CrossRef] [PubMed]
  42. Taylor, J.P.; Brown, R.H.J.; Cleveland, D.W. Decoding ALS: From Genes to Mechanism. Nature 2016, 539, 197–206. [Google Scholar] [CrossRef] [PubMed]
  43. Li, W.; Jin, Y.; Prazak, L.; Hammell, M.; Dubnau, J. Transposable Elements in TDP-43-Mediated Neurodegenerative Disorders. PLoS ONE 2012, 7, e44099. [Google Scholar] [CrossRef] [PubMed]
  44. Kawai, T.; Akira, S. Toll-like Receptors and Their Crosstalk with Other Innate Receptors in Infection and Immunity. Immunity 2011, 34, 637–650. [Google Scholar] [CrossRef]
  45. Stetson, D.B. Endogenous Retroelements and Autoimmune Disease. Curr. Opin. Immunol. 2012, 24, 692–697. [Google Scholar] [CrossRef] [PubMed]
  46. Cassel, S.L.; Sutterwala, F.S. Sterile Inflammatory Responses Mediated by the NLRP3 Inflammasome. Eur. J. Immunol. 2010, 40, 607–611. [Google Scholar] [CrossRef]
  47. Gugliandolo, A.; Giacoppo, S.; Bramanti, P.; Mazzon, E. NLRP3 Inflammasome Activation in a Transgenic Amyotrophic Lateral Sclerosis Model. Inflammation 2018, 41, 93–103. [Google Scholar] [CrossRef] [PubMed]
  48. He, Y.; Hara, H.; Núñez, G. Mechanism and Regulation of NLRP3 Inflammasome Activation. Trends Biochem. Sci. 2016, 41, 1012–1021. [Google Scholar] [CrossRef] [PubMed]
  49. Cohen, T.J.; Lee, V.M.Y.; Trojanowski, J.Q. TDP-43 Functions and Pathogenic Mechanisms Implicated in TDP-43 Proteinopathies. Trends Mol. Med. 2011, 17, 659–667. [Google Scholar] [CrossRef]
  50. Zhao, W.; Beers, D.R.; Bell, S.; Wang, J.; Wen, S.; Baloh, R.H.; Appel, S.H. TDP-43 Activates Microglia through NF-ΚB and NLRP3 Inflammasome. Exp. Neurol. 2015, 273, 24–35. [Google Scholar] [CrossRef]
  51. Leal-Lasarte, M.M.; Franco, J.M.; Labrador-Garrido, A.; Pozo, D.; Roodveldt, C. Extracellular TDP-43 Aggregates Target MAPK/MAK/MRK Overlapping Kinase (MOK) and Trigger Caspase-3/IL-18 Signaling in Microglia. FASEB J. 2017, 31, 2797–2816. [Google Scholar] [CrossRef]
  52. Rawat, P.; Teodorof-Diedrich, C.; Spector, S.A. Human Immunodeficiency Virus Type-1 Single-Stranded RNA Activates the NLRP3 Inflammasome and Impairs Autophagic Clearance of Damaged Mitochondria in Human Microglia. Glia 2019, 67, 802–824. [Google Scholar] [CrossRef] [PubMed]
  53. Krestel, H.; Meier, J.C. RNA Editing and Retrotransposons in Neurology. Front. Mol. Neurosci. 2018, 11, 163. [Google Scholar] [CrossRef] [PubMed]
  54. Calvo, A.C.; Cibreiro, G.A.; Merino, P.T.; Roy, J.F.; Galiana, A.; Rufián, A.J.; Cano, J.M.; Martín, M.A.; Moreno, L.; Larrodé, P.; et al. Collagen XIX Alpha 1 Improves Prognosis in Amyotrophic Lateral Sclerosis. Aging Dis. 2019, 10, 278–292. [Google Scholar] [CrossRef] [PubMed]
  55. Oros, D.; Strunk, M.; Breton, P.; Paules, C.; Benito, R.; Moreno, E.; Garcés, M.; Godino, J.; Schoorlemmer, J. Altered Gene Expression in Human Placenta after Suspected Preterm Labour. Placenta 2017, 55, 21–28. [Google Scholar] [CrossRef] [PubMed]
  56. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  57. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
Figure 1. HERVK expression in ALS. Expression analysis of human endogenous retrovirus ENV and POL was carried out via qPCR in cDNA samples from ALS patients (ALS), patients with other neuropathies (ONPs), and controls (NDIs). Results were normalized using RPL19 and GAPDH as reference genes, calculated using the 2−ΔΔCt method, and represented as the fold expression compared with the mean expression level in the controls. ENV (A) and POL (B) relative expression levels in PBMC samples from 10 ALS patients were compared with NDI samples (n = 10) and ONP samples (n = 10) (Student’s t-test; p = 0.982). Correlations (C,D) between mRNA expression of HERVK ENV and POL genes in the PBMC samples analyzed in (A,B) are represented, and Pearson’s correlation coefficients are indicated. ENV (E,G) and POL (F,H) levels were analyzed in the cortex and brainstem from ALS patients and control individuals (NDI). Cortex samples: ALS, n = 10; NDI, n = 6; brainstem samples: ALS, n = 8; NDI, n = 8.
Figure 1. HERVK expression in ALS. Expression analysis of human endogenous retrovirus ENV and POL was carried out via qPCR in cDNA samples from ALS patients (ALS), patients with other neuropathies (ONPs), and controls (NDIs). Results were normalized using RPL19 and GAPDH as reference genes, calculated using the 2−ΔΔCt method, and represented as the fold expression compared with the mean expression level in the controls. ENV (A) and POL (B) relative expression levels in PBMC samples from 10 ALS patients were compared with NDI samples (n = 10) and ONP samples (n = 10) (Student’s t-test; p = 0.982). Correlations (C,D) between mRNA expression of HERVK ENV and POL genes in the PBMC samples analyzed in (A,B) are represented, and Pearson’s correlation coefficients are indicated. ENV (E,G) and POL (F,H) levels were analyzed in the cortex and brainstem from ALS patients and control individuals (NDI). Cortex samples: ALS, n = 10; NDI, n = 6; brainstem samples: ALS, n = 8; NDI, n = 8.
Ijms 25 01549 g001
Figure 2. HERVK copies expressed in ALS PBMCs and the brain. A qPCR that specifically amplifies HERVK cDNA (described in Figure 1) was applied to PBMCs and brain samples taken from ALS patients and healthy controls (Table 1 and Table 2). Products were sequenced via next-generation sequencing (NGS), and the reads were mapped to the HERVK copies in the human genome. The relative frequency (%) of individual HERVK copies was calculated as the number of reads mapping to a particular HERVK locus relative to the total number of mapped reads. This relative frequency is a measure of the relative transcription of each individual HERVK copy. The 27 HERVK copies to which mapping was counted are listed on the right; their genomic location can be found in Table S2. Relative frequencies are indicated in pie charts; color codes for each copy are indicated on the right. Analysis was performed on either the PBMCs of ALS patients (n = 10) and healthy controls (n = 10) (A) or on postmortem samples of the brainstems of ALS patients or non-diseased individuals (NDIs) (n = 4 in each group). Both the overall similarity between copies expressed in ALS vs. controls and the differentially expressed individual copies can be discerned. (A) Pie charts indicate the relative frequency (%) of reads mapping to each of the 27 HERVK copies in PBMC. (B) Same for brainstems.
Figure 2. HERVK copies expressed in ALS PBMCs and the brain. A qPCR that specifically amplifies HERVK cDNA (described in Figure 1) was applied to PBMCs and brain samples taken from ALS patients and healthy controls (Table 1 and Table 2). Products were sequenced via next-generation sequencing (NGS), and the reads were mapped to the HERVK copies in the human genome. The relative frequency (%) of individual HERVK copies was calculated as the number of reads mapping to a particular HERVK locus relative to the total number of mapped reads. This relative frequency is a measure of the relative transcription of each individual HERVK copy. The 27 HERVK copies to which mapping was counted are listed on the right; their genomic location can be found in Table S2. Relative frequencies are indicated in pie charts; color codes for each copy are indicated on the right. Analysis was performed on either the PBMCs of ALS patients (n = 10) and healthy controls (n = 10) (A) or on postmortem samples of the brainstems of ALS patients or non-diseased individuals (NDIs) (n = 4 in each group). Both the overall similarity between copies expressed in ALS vs. controls and the differentially expressed individual copies can be discerned. (A) Pie charts indicate the relative frequency (%) of reads mapping to each of the 27 HERVK copies in PBMC. (B) Same for brainstems.
Ijms 25 01549 g002
Figure 3. Altered expression of HERVK copies in ALS brains. (A) Fold change or fold difference between all 27 HERVK loci. The change in each locus (Δfold) was calculated as the mean/average relative frequency (%) of reads mapped to an individual HERVK locus in ALS patients versus NDIs (or the inverse for results < 1). Subsequently, the fold change (ΔΔfold) was calculated for each locus as Δfold-1. The median fold changes over all loci were calculated for PBMCs, the cortex (Cx), and the brainstem (Stem) and are represented in the figure. The corresponding p-values are indicated above. (B) Same analysis as in C for 17 loci showing differences in brains between ALS patients and NDIs. Loci that showed 0.8 < Δfold < 1.2 between ALS patients and controls for both the cortex (Cx) and brainstem were excluded from the analysis. (C) The relative frequency (%) of reads mapped to the individual HERVK copies is indicated. The bottom and top of the box plot are, respectively, the first and third quartiles, and the band inside the box is the second quartile (the median). Lines extending vertically from the boxes indicate variability outside the upper and lower quartiles. Outliers are indicated as individual points. IR: interquartile range. The statistical significance for the difference between ALS and NDI (calculated using the DESeq2 pipeline) is shown in the picture (* adjusted p-value (p-adj.) < 0.05).
Figure 3. Altered expression of HERVK copies in ALS brains. (A) Fold change or fold difference between all 27 HERVK loci. The change in each locus (Δfold) was calculated as the mean/average relative frequency (%) of reads mapped to an individual HERVK locus in ALS patients versus NDIs (or the inverse for results < 1). Subsequently, the fold change (ΔΔfold) was calculated for each locus as Δfold-1. The median fold changes over all loci were calculated for PBMCs, the cortex (Cx), and the brainstem (Stem) and are represented in the figure. The corresponding p-values are indicated above. (B) Same analysis as in C for 17 loci showing differences in brains between ALS patients and NDIs. Loci that showed 0.8 < Δfold < 1.2 between ALS patients and controls for both the cortex (Cx) and brainstem were excluded from the analysis. (C) The relative frequency (%) of reads mapped to the individual HERVK copies is indicated. The bottom and top of the box plot are, respectively, the first and third quartiles, and the band inside the box is the second quartile (the median). Lines extending vertically from the boxes indicate variability outside the upper and lower quartiles. Outliers are indicated as individual points. IR: interquartile range. The statistical significance for the difference between ALS and NDI (calculated using the DESeq2 pipeline) is shown in the picture (* adjusted p-value (p-adj.) < 0.05).
Ijms 25 01549 g003
Figure 4. Inflammation-related gene expression in ALS brain. Expression analysis of the genes indicated was carried out via qPCR on cDNA samples prepared from the brain samples of ALS patients and controls (NDIs). Results were normalized using GAPDH and RPL19 as reference genes, calculated using the 2−ΔΔCt method, and are represented as a percentage of the fold expression compared with the median expression level in the controls. NLRP3 expression was analyzed in both cortex and brainstem (A) and in both tissues together (B). ALS patient samples and control (NDI) samples were prepared from either the cerebral cortex or brainstem. Cortex ALS, n = 10; cortex NDI, n = 3; brainstem ALS, n = 9; brainstem NDI, n = 4. Student’s t-test: * p < 0.05, *** p < 0.001.
Figure 4. Inflammation-related gene expression in ALS brain. Expression analysis of the genes indicated was carried out via qPCR on cDNA samples prepared from the brain samples of ALS patients and controls (NDIs). Results were normalized using GAPDH and RPL19 as reference genes, calculated using the 2−ΔΔCt method, and are represented as a percentage of the fold expression compared with the median expression level in the controls. NLRP3 expression was analyzed in both cortex and brainstem (A) and in both tissues together (B). ALS patient samples and control (NDI) samples were prepared from either the cerebral cortex or brainstem. Cortex ALS, n = 10; cortex NDI, n = 3; brainstem ALS, n = 9; brainstem NDI, n = 4. Student’s t-test: * p < 0.05, *** p < 0.001.
Ijms 25 01549 g004
Table 1. List of PBMC samples analyzed in this study.
Table 1. List of PBMC samples analyzed in this study.
Patient CharacteristicsfALS Patients
(n = 10)
ONP Patients
(n = 10)
NDIs
(n = 10)
Gender (n)6 males
4 females
6 males
4 females
4 males
6 females
Age at illness onset
(mean ± SD)
60.04 ± 11.11
Disease duration, months
(mean ± SD)
34.27 ± 20.73
Age at sampling
(mean ± SD)
66.89 ± 11.6152.37 ± 11.6156.97 ± 10.55
NDI: non-diseased individuals group; ONP: Other neuropathies group; fALS: familial amyotrophic lateral sclerosis group (all the patients, except for two cases, carried mutations in one of the following four genes: SOD1/TDP43/FUS/C9orf72).
Table 2. List of brain samples used in this study.
Table 2. List of brain samples used in this study.
Patient CharacteristicssALS Patients
(n = 10)
NDI
(n = 9)
Brain anatomical region (n)10 frontal cortexes
8 brainstems
6 frontal cortexes
8 brainstems
Gender (n)6 males
4 females
6 males
3 females
Age at sampling (mean ± SD)58.64 ± 8.3759.00 ± 8.73
NDI: Non-diseased individuals group; sALS: sporadic amyotrophic lateral sclerosis group.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Moreno-Martinez, L.; Macías-Redondo, S.; Strunk, M.; Guillén-Antonini, M.I.; Lunetta, C.; Tarlarini, C.; Penco, S.; Calvo, A.C.; Osta, R.; Schoorlemmer, J. New Insights into Endogenous Retrovirus-K Transcripts in Amyotrophic Lateral Sclerosis. Int. J. Mol. Sci. 2024, 25, 1549. https://doi.org/10.3390/ijms25031549

AMA Style

Moreno-Martinez L, Macías-Redondo S, Strunk M, Guillén-Antonini MI, Lunetta C, Tarlarini C, Penco S, Calvo AC, Osta R, Schoorlemmer J. New Insights into Endogenous Retrovirus-K Transcripts in Amyotrophic Lateral Sclerosis. International Journal of Molecular Sciences. 2024; 25(3):1549. https://doi.org/10.3390/ijms25031549

Chicago/Turabian Style

Moreno-Martinez, Laura, Sofía Macías-Redondo, Mark Strunk, María Isabel Guillén-Antonini, Christian Lunetta, Claudia Tarlarini, Silvana Penco, Ana Cristina Calvo, Rosario Osta, and Jon Schoorlemmer. 2024. "New Insights into Endogenous Retrovirus-K Transcripts in Amyotrophic Lateral Sclerosis" International Journal of Molecular Sciences 25, no. 3: 1549. https://doi.org/10.3390/ijms25031549

APA Style

Moreno-Martinez, L., Macías-Redondo, S., Strunk, M., Guillén-Antonini, M. I., Lunetta, C., Tarlarini, C., Penco, S., Calvo, A. C., Osta, R., & Schoorlemmer, J. (2024). New Insights into Endogenous Retrovirus-K Transcripts in Amyotrophic Lateral Sclerosis. International Journal of Molecular Sciences, 25(3), 1549. https://doi.org/10.3390/ijms25031549

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop