Skip to main content

Advertisement

Dynamic LTR retrotransposon transcriptome landscape in septic shock patients

Abstract

Background

Sepsis is defined as a life-threatening organ dysfunction caused by a dysregulated host response to infection. Numerous studies have explored the complex and dynamic transcriptome modulations observed in sepsis patients, but a large fraction of the transcriptome remains unexplored. This fraction could provide information to better understand sepsis pathophysiology. Multiple levels of interaction between human endogenous retroviruses (HERV) and the immune response have led us to hypothesize that sepsis is associated with HERV transcription and that HERVs may contribute to a signature among septic patients allowing stratification and personalized management.

Methods

We used a high-density microarray and RT-qPCR to evaluate the HERV and Mammalian Apparent Long Terminal Repeat retrotransposons (MaLR) transcriptome in a pilot study that included 20 selected septic shock patients, stratified on mHLA-DR expression, with samples collected on day 1 and day 3 after inclusion. We validated the results in an unselected, independent cohort that included 100 septic shock patients on day 3 after inclusion. We compared septic shock patients, according to their immune status, to describe the transcriptional HERV/MaLR and conventional gene expression. For differential expression analyses, moderated t tests were performed and Wilcoxon signed-rank tests were used to analyze RT-qPCR results.

Results

We showed that 6.9% of the HERV/MaLR repertoire was transcribed in the whole blood, and septic shock was associated with an early modulation of a few thousand of these loci, in comparison to healthy volunteers. We provided evidence that a subset of HERV/MaLR and conventional genes were differentially expressed in septic shock patients, according to their immune status, using monocyte HLA-DR (mHLA-DR) expression as a proxy. A group of 193 differentially expressed HERV/MaLR probesets, tested in an independent septic shock cohort, identified two groups of patients with different immune status and severity features.

Conclusion

We demonstrated that a large, unexplored part of our genome, which codes for HERV/MaLR, may be linked to the host immune response. The identified set of HERV/MaLR probesets should be evaluated on a large scale to assess the relevance of these loci in the stratification of septic shock patients. This may help to address the heterogeneity of these patients.

Background

Sepsis is defined as a life-threatening organ dysfunction caused by a dysregulated host response to infection, with septic shock being the most severe subtype of sepsis [1]. Despite major progress, due in large part to source control with antibiotics and early initiation of intensive care therapy, sepsis remains a major health issue, with a high worldwide prevalence and high mortality rate. It is characterized by immune dysfunctions with concomitant excessive pro- and anti-inflammatory responses, which can lead to organ failure, immunoparalysis, and secondary infections. The pathophysiological mechanisms of sepsis are not completely understood. The early model of an overwhelming inflammatory reaction failed to capture the complex pathophysiology of the syndrome, as proven by the failure of multiple clinical trials with a variety of anti-inflammatory agents [2]. While an immune dysregulated host response is clearly part of the pathophysiology, the complex interactions between the immune response and other key physiological systems, such as the autonomous nervous system, coagulation, or cellular bioenergetics, are still poorly understood [3, 4]. This translates into an apparent heterogeneity observed in patients, which makes the selection of appropriate therapeutic care a major clinical challenge. Several markers of the immune status have been explored to better stratify patients [5]. There are currently no gold standard or clinical signs to evaluate immunosuppression in the intensive care unit. To date, HLA-DR expression in monocytes, measured by flow cytometry, is a well-accepted and beneficial marker for monitoring immune alterations in critically ill patients (sepsis, trauma, pancreatitis, surgery, and burns). The decreased expression of mHLA-DR has been repeatedly associated with mortality and secondary infections and therefore remains an independent predictor of a poor outcome after sepsis [4, 6,7,8,9,10].

Long terminal repeat (LTR) retrotransposons, which represent 8.3% of the human genome [11] were recently hypothesized to be relevant in the pathophysiology of sepsis [12, 13]. These retrovirus-like sequences consist of approximately 200,000 human endogenous retroviruses (HERVs) and 240,000 mammalian apparent long terminal repeat retrotransposons (MaLR). The expression of HERVs has been observed in both inflammation [14, 15] and immunosuppression [13, 16], and the insertional polymorphism of HERV LTRs at the HLA locus has been associated with several auto-immune diseases [17]. In addition, several lines of evidence support the role of HERVs as contributors to the immune response [18]. For example, loss of DNA methylation triggers cytosolic sensing of double-stranded RNA (dsRNA) and overexpression of HERVs, both of which cause a type I interferon response and apoptosis [19]. Another study in cancer highlighted that overexpression of HERVs triggers pathologic innate immune signaling [20]. Chuong et al. also showed that HERV elements constitute a dynamic reservoir of IFN-inducible enhancers that allow the regulation of essential immune functions [21]. These data led us to hypothesize that sepsis may be associated with HERV transcription, which may modulate the immune response. Moreover, irrespective of causal or consecutive expression, HERVs could be part of an informative molecular signature that may allow a better stratification of septic shock patients and improved patient management.

We have recently shown that several HERVs are expressed and modulated after septic shock and other inflammatory injuries [12]. Using the HERV-V3 microarray, which allows the measurement of HERV/MaLR at the individual locus level [22], we have also shown that HERV/MaLR are tightly modulated in peripheral blood mononuclear cells (PBMCs). We have used an ex vivo endotoxin tolerance model to mimic the septic shock host immune response [13]. The data from this pilot study were used to investigate the HERV/MaLR transcriptome in septic shock patients, who were stratified by monocyte HLA-DR (mHLA-DR) expression on day 3, following intensive care unit (ICU) admission (Fig. 1). We provided a global view of the HERV transcriptome in the whole blood of healthy volunteers and septic shock patients. We also aimed to identify whether a subset of HERV/MaLR was expressed in septic shock patients, according to their immune status, as estimated by mHLA-DR expression. Finally, we highlighted a group of 193 differentially expressed HERV/MaLR probesets, which identified two groups of patients with different immune status and severity features, in an independent septic shock cohort.

Fig. 1
figure1

Schematic representation of the data process analysis. Healthy volunteers and both septic shock patient cohort subsets (IMMUNOSEPSIS and MIP-Rea) are defined on the left. Bioinformatics analysis parameters used for each step are described in the middle. Finally, the three phases of the data process analysis (descriptive phase, discovery phase, and validation phase) are presented. TP time point, D day

Methods

Biological samples

IMMUNOSEPSIS cohort subset

We retrospectively selected patients from a prospective, non-interventional study conducted in Lyon, France, named IMMUNOSEPSIS. A cohort of adult septic shock patients was enrolled from December 2001 to April 2005 in two ICUs from a French university hospital. The clinical study was approved by the regional ethics committee (Comité Consultatif de Protection des Personnes dans la Recherche BioMédicale de Lyon A). The committee waived the need for written informed consent because the study was observational, with a low risk to patients, and no specific procedure, other than routine blood sampling, was required. Oral information and non-opposition to inclusion in the study were mandatory and were systematically obtained before any blood sample was drawn. This was recorded in patients’ clinical files. The cohort has been described in detail elsewhere [23, 24]. In brief, blood samples were collected in PAXgene™ Blood RNA tubes (PreAnalytix) and RNA was extracted and frozen at − 80 °C, as previously described [25]. For the current study, we selected 20 patients for which the mHLA-DR expression level was measured on the monocyte cell surface by flow cytometry from the peripheral whole blood collected in EDTA anticoagulant tubes, as previously described [26, 27]. We chose 10 patients with a normal expression level of mHLA-DR on day 3 (more than 30% of expression on monocytes cell surface) and 10 patients with a low expression level of mHLA-DR on day 3 (less than 30% of expression on monocytes cell surface), for which RNA samples were available on day 1 and day 3 or 4 (to simplify, named day 3, thereafter). The RNA integrity was measured prior to RNA amplification using a Bioanalyser 2100 (Agilent), in accordance with the manufacturer’s instructions (median RNA integrity number 8.5 [min 6.5–max 9.4]). Patient characteristics of the IMMUNOSEPSIS cohort subset are described in Table 1.

Table 1 Patient characteristics of the IMMUNOSEPSIS cohort subset according to the mHLA-DR expression on day 3 (descriptive and discovery phase)

Healthy controls

PAXgene™ Blood RNA tubes (PreAnalytix) from four healthy volunteers (HV), matched in age (median at 52 years [Q1:51–Q3:54]) and sex (75% were male), were independently obtained from EFS (Etablissement Français du Sang) and used immediately. We used the EFS standardized procedures for blood donation and followed provisions of articles R.1243–49 and the French public health code to obtain written non-opposition to the use of donated blood for research purposes from HV. The blood donors’ personal data were deidentified before transfer to our research laboratory. We obtained the favorable notice of the local ethical committee (Comité de Protection des Personnes Sud-Est II, Bâtiment Pinel, 59 Boulevard Pinel, 69,500 Bron) and acceptance from the French ministry of research (Ministère de lʼEnseignement supérieur, de la Recherche et de lʼInnovation, DC-2008-64) for the handling and conservation of these samples. Blood samples were stabilized for at least 4 h at room temperature, after collection, and frozen at − 80 °C, following the manufacturer’s guidelines. Total RNA was extracted using the PAXgene Blood RNA kit (PreAnalytix), in accordance with the manufacturer’s instructions. The RNA quantity and quality were determined using a Nanodrop (Thermo Scientific) Bioanalyser 2100 (Agilent), in accordance with the manufacturer’s instructions. Samples with an RNA integrity number ≤ 6 were excluded due to poor quality RNA.

MIP-Rea cohort subset

We retrospectively selected patients from a prospective, multicenter, non-interventional study, named MIP-Rea, that was conducted in six ICUs in Lyon, France. The study was approved by our institutional ethical review board (Comité d’Ethique des Centres d’investigation Clinique de l’Inter-Région Rhône-Alpes Auvergne – IRB 5044), and consent for an ancillary study was obtained a posteriori. The protocol of this study was submitted to the French CCTIRS and CNIL committees and approved on April 22, 2016, and September 30, 2016, respectively. Further, informed consent was received from patients for inclusion in this specific study. The cohort has been described in detail elsewhere [28]. In brief, blood samples were collected in PAXgene™ Blood RNA tubes (PreAnalytix) and RNA was extracted and frozen at − 80 °C, as previously described [28]. For the current study, we selected only septic shock patients for which RNA samples were available on day 1, day 3, and day 6 (corresponding to a total of 100 septic shock patients). The RNA integrity was measured prior to RNA amplification using a Bioanalyser 2100 (Agilent), in accordance with the manufacturer’s instructions (median RNA integrity number 7.9 [min 6–max 9.6]). Patient characteristics for the MIP-Rea cohort subset are described in Table 2.

Table 2 Patient characteristics of the MIP-Rea cohort subset according to cluster affiliation described in the “Results” section (validation phase)

HERV-V3 processing and analysis

Custom Affymetrix HERV-V3 GeneChip microarray

The HERV-V3 microarray targeted 353,994 loci-elements, which were represented by 4,410,200 probes [22]. The custom HERV GeneChip can discriminate distinct HERV elements through the use of a set of highly informative probesets (located in U3, R, U5 subdomains of solo, 5′ and 3′ individual LTRs and gag/pol/env regions), hereafter named the “HERV prototypes repertoire.” A group of probesets with lower quality annotations (located in the first third and last third of the complete LTR and every 2.5 kb in the region in between the LTRs) was also used and is hereafter referred to as the “HERV/MaLR_Dfam repertoire.” The custom HERV GeneChip also contained probesets that target 1500 genes involved in the immune response. The description of the HERVgDB4 database and the final content of the HERV-V3 microarray are provided in Additional file 1: Table S1 and were previously described in Becker et al. [22].

RNA amplification, labeling, and HERV-V3 microarray hybridization

The cDNA synthesis and amplification steps were performed from 16 ng of RNA, using the Ovation Pico WTA System V2 kit (Nugen), in accordance with the manufacturer’s instructions. Five micrograms of purified, amplified cDNA was fragmented into 50–200 bp fragments and were 3′-labeled using the Encore Biotin Module kit (Nugen), in accordance with the manufacturer’s instructions. The HERV-V3 microarrays were hybridized at 50 °C for 18 h in an oven, with constant stirring (60 rpm). Washing and staining were performed using the protocol supplied by the manufacturer, using the GeneChip fluidic station 450 (Affymetrix). The arrays were scanned using the GeneChip fluorometric scanner 3000 7G (Affymetrix). Images (DAT files) were converted to CEL files using the GCOS software (Affymetrix) [22]. The experimental data generated were deposited in the National Center for Biotechnology Information (NCBI) and are available in the GEO DataSets site, under accession number GSE121352.

Microarray analysis preprocessing

The CEL files were transformed into a matrix, then normalized and adjusted for background noise (RMA normalization). Probes were summarized into probesets with command apt-probeset-summarize (v1.18.0). Microarray preprocessing and statistical analyses were performed using R/Bioconductor (R v3.1.2) [29]. Quality assessment was performed through simpleaffy (v2.42.0) [30] and arrayQualityMetrics (v3.22.1) [31]. For quality control, several criteria were used: RNA quality, images of chips, hybridization spike-in, polyA, amplification and fragmentation, intensity signals (before and after normalization), probeset homogeneity (RLE, NUSE plots), correlation plots (before and after normalization), and Principal Component Analysis. For most criteria, outlier detection was performed by computing the Kolmogorov Smirnof (KS) statistic between each array and the pooled data (default threshold with arrayQualityMetrics library). One array, which corresponded to a HV sample, was removed from the analysis as it did not pass more than four quality controls. This left four HV samples for the analysis. Experiment batch effects were removed using COMBAT [32]. Finally, low expression probesets were filtered to reduce the dataset size and gain statistical power for the analyses. A coefficient of variation (CV) value of 10% was used to determine the intensity threshold (25.5) to filter such probesets with low expression. Probesets under this intensity threshold in more than 90% of all samples (40 samples out of the 44) were removed, which left 120,222 probesets for the analysis. LTR function analysis and RT-qPCR validation not described in the following section are presented in the Additional file 2: Supplementary Methods.

Bioinformatics analysis

For differential expression analyses, moderated t tests were performed (Limma (v3.22.7) R package) [33], and P values were adjusted for multiple testing using the Benjamini-Hochberg procedure [34]. When we compared septic shock patients and HV, we used absolute log2 Fold Change (|log2FC|) higher than 1 and an adjusted P value lower than 0.05. When we compared septic shock patient groups, according to the day of sampling and mHLA-DR, a probeset was considered to show statistically significant differential expression when the |log2FC| was higher than 1 and adjusted P value was lower than 0.1. Graphics were drawn using ggplot2 (v2.2.0) or pheatmap (v1.0.8). The principal component analysis was performed on the HERV expression matrix. The representation on the first two components was carried out using the ggplot2 library. The two vectors drawn on the plot are the median of coordinate differences between day 1 and day 3 for each patient, according to their mHLA-DR status. Finally, the 193 HERV/MaLR probesets that were modulated on day 3, between low and normal mHLA-DR patients from the discovery phase, were tested on the independent cohort using unsupervised hierarchical clustering. The complete clustering method was applied, with Euclidean distance for rows (probesets) and Pearson’s correlation distance for columns (samples).

Results

HERV/MaLR transcriptome landscape in the whole blood

We assessed the HERV/MaLR transcriptome from whole blood in a selected IMMUNOSEPSIS cohort subset [23, 24] that consisted of 20 septic shock patients (samples collected at admission and on day 3) and four HV (named descriptive phase in Fig. 1). A summary of the absolute counts and the relative abundances of transcriptionally active elements of the HERV/MaLR transcriptome are given in Fig. 2a (and extended in Additional file 1: Table S1). Overall, 6.9% (87,912 probesets) of HERV/MaLR were transcriptionally active in the whole blood, which was relatively low compared to the observed percentage of expressed genes (42%) among the set of 1500 immune-related genes (32,310 probesets). All well-defined HERV groups had expressed loci in the whole blood, and we did not observe significant enrichment or depletion of a specific repertoireFootnote 1, classFootnote 2, or HERV groupFootnote 3, between HV and septic shock patients (data not shown). Among the expressed prototype elements, 7.4% belonged to gamma-retroviruses and were mostly found in the HERV-H group (1042 expressed probesets, which corresponded to 19% of the whole group). It should be noted that the HERV-FRD and PRIMA41 groups contained a high absolute number of expressed loci (748 and 732 probesets, respectively). In beta-retroviridae, 2.1% of elements were expressed, mostly from the HML-8 and HML-1 groups (861 probesets, 14% of the whole group; and 216 probesets, 14% of the whole group, respectively). Finally, 2% of the expressed elements belonged to spuma/epsilon-like retrovirus classes, in which the HERV-L group provided the largest count of expressed probesets (1396 probesets, which corresponded to 11% of the whole group). The overall HERV/MaLR transcriptome analysis from the whole blood was also carried out on 100 septic shock patients from the independent and unselected MIP-REA cohort subset [28] (named validation phase in Fig. 1), and the results were similar (Additional file 3: Figure S1). Interestingly, 92% of expressed probesets were common to both cohorts.

Fig. 2
figure2

HERV/MaLR transcriptome landscape in the whole blood (descriptive phase). a Percentages and absolute counts (number presented inside the bars) of expressed loci (probeset intensity above threshold), within individual groups of “HERV_prototypes,” “HERV_Dfam,” and “MaLR_Dfam” repertoires. “HERV_prototypes” were grouped by classes of retroviruses, namely gamma-retrovirus (green), beta-retrovirus (pink), and spuma-epsilon like retrovirus (blue). “HERV and MaLR Dfam” repertoires are depicted as a global homogeneous entity (purple). b Volcano plots derived from the HERV/MaLR differential expression analysis, between septic shock patients and HV on day 1 (left) or day 3 (right). The x-axis represents the log2 fold change values and the y-axis represents the -log10 adjusted P values. Each point represents a given probeset. Statistically and biologically significant probesets (adjusted P value < 0.05, │log2FC│ > 1) are colored (red = downregulated, green = upregulated). c Venn diagrams from the HERV/MaLR differential expression analyses, according to the day of sampling and mHLA-DR expression (compared to HV). Upregulated probesets are in green and downregulated probesets are in red

Beyond the overall transcription analysis, we looked for any difference due to sequence composition (LTRsFootnote 4, versus proviral genes) or function (promoter versus polyA). In Additional file 4: Figure S2, we summarized the HERV/MaLR transcriptomic activity and LTR putative regulatory functions observed in HV and patients. We found enrichment in expressed LTRs when compared to the expressed proviral regions of HERVs. We also found that 22% of LTRs were either putative promoters (Pr) or polyadenylation (pA) signals. About one third of LTRs shifted from a silent to active state (Pr or pA) between samples used in the descriptive phase. Notably, the LTR switch of function from promoter to polyA appeared to be extremely rare. This supported our previous assumption that an LTR is predetermined to act as a promoter or a polyA signal, which is a phenomenon that we named “operational determinism” [13, 35]. Overall, these results confirmed, in septic shock patients, what was previously observed in an endotoxin tolerance model in PBMC in vitro [13].

We then performed a differential expression analysis between the HV and septic shock patients on either day 1 or day 3. The exhaustive dataset of this analysis is presented in Additional file 5: Table S2. Overall, we found 855 differentially expressed HERV/MaLR probesets between combined patient samples (integrating day of sampling and mHLA-DR status) and HV. Volcano plots were depicted at the probeset level for HERV/MaLR (Fig. 2b). There were approximately twice as many downregulated as upregulated HERV/MaLR, when comparing septic shock patients on day 1 with HV (513 upregulated probesets vs 1660 downregulated probesets). The difference between the two groups appeared to be smaller on day 3 (553 upregulated probesets vs 339 downregulated probesets) (Fig. 2b, c. When mHLA-DR expression was added as a variable, the overall number of downregulated HERV/MaLR probesets was more higher in septic shock patients with low mHLA-DR expression versus normal mHLA-DR expression compared to HV on day 1 (931 downregulated probesets vs 298 upregulated probesets) but also on day 3 (323 downregulated probesets vs 12 upregulated probesets) (Fig. 2c and Additional file 5: Table S2). Globally, we confirmed that septic shock is associated with an early specific modulation of the expression of a few thousands loci from the LTR retrotransposon repertoire of the human genome.

The differential expression of HERV/MaLR according to immune status of septic shock patients

To gain an insight into the modulation of genes and HERV/MaLR expression, according to a well-accepted proxy of immune status in patients (i.e., HLA-DR expression on monocytes), we first performed a principal component analysis (PCA) (Fig. 3a). The PCA showed that all septic shock patients were clearly separated from HV, as demonstrated above. The low versus normal mHLA-DR groups appeared to behave differently, as shown by day 1 to day 3 orthogonal progressions for the two sub-populations. These depicted that the patients with normal mHLA-DR on day 3 tended to be closer to the HV group. In addition, when we explored the differentially expressed HERV/MaLR probesets between patients, by day of sampling and by mHLA-DR status, the first observation related to the incapacity to differentiate, on day 1, low versus normal mHLA-DR patients (Fig. 3b). A Venn diagram analysis showed that the majority of HERV/MaLR probesets appeared to be downregulated on day 3 in patients with a low, compared to normal expression of mHLA-DR (166 downregulated and 27 upregulated probesets, a total of 193 HERV/MaLR probesets, which corresponded to 162 HERV/MaLR loci).

Fig. 3
figure3

The differential expression of HERV/MaLR according to the immune status of septic shock patients (discovery phase). a Principal component analysis developed from the HERV expression matrix. Healthy volunteers, patients with normal expression of mHLA-DR, and patients with low expression of mHLA-DR are indicated in gray, pink, and blue, respectively. b Venn diagrams from HERV/MaLR differential expression analyses, according to the day of sampling and mHLA-DR expression (compared between septic shock patient groups). Upregulated probesets are in green and downregulated probesets are in red. c The table shows the number of statistically significant differentially expressed loci (DELs) for HERV/MaLR, differentially expressed genes (DEGs), and the corresponding number of probesets, respectively. Downregulated loci are in red, upregulated loci are in green. For HERV/MaLR loci, the name, number of differentially expressed probesets (between brackets), and chromosomal locations (in italic) are indicated (GRCh38 genome version) for the five most differentially expressed loci. For genes, the official gene symbol and the number of differentially expressed probesets (between brackets) are indicated

The expression of HERV/MaLR elements was investigated and, when we compared low versus normal mHLA-DR expression on day 3, no group appeared to be enriched or depleted (Fig. 3c). Nevertheless, of the 193 differentially expressed HERV/MaLR probesets, according to mHLA-DR expression, 138 were also differentially expressed between septic shock patients on day 3 and HV. Hence, 55 HERV/MaLR probesets showed specific differential expression within septic shock patients, according to mHLA-DR expression level. Differential expression analysis of low versus normal mHLA-DR conditions also identified a downregulation of genes involved in antigen-presentation; TCR activation; differentiation of effector CD8+ T cells, such as CD3G, DPP4, and EOMES; and chemokines and cytokine receptors, such as CCL5 and IL5RA. We found that IL18R1 and CYP1B1 genes, which are involved in cytokine signaling and cell proliferation, migration and survival, HGF, IL1R2, and MERTK genes, which are involved in cell growth, apoptosis, and the TLR innate immune response, were all upregulated in low mHLA-DR samples (Fig. 3c).

We validated the relevance of the loci identified by the microarray using RT-qPCR. We selected 29 HERV/MaLR loci that were significantly differentially expressed in patients on day 3, according to the mHLA-DR status. The HERV/MaLR locus-specific RT-qPCR systems were meticulously designed and validated on genomic DNA matrix to ensure locus specificity (Additional file 6: Figure S3). Additional file 7: Figure S4 illustrates the consistency of 27 out of the 29 microarrays and RT-qPCR profiles. We found three HERV/MaLR profiles that allowed us to determine the immune status of patients on day 3 and sometimes on day 1. Approximately three quarters of the HERV/MaLR loci had a lowered expression in patients with low mHLA-DR, while the remainder conversely showed a significant increase in expression in the same patients. Interestingly, on day 1, some HERV/MaLR loci were able to stratify patients who had normal or low mHLA-DR on day 3, when no difference in mHLA-DR expression was observed on day 1. Consequently, we have identified a link between immune profiles, and the LTR retrotransposon expression.

HERV/MaLR loci can differentiate septic shock patients according to severity status

Finally, we sought to evaluate whether the 193 HERV/MaLR probesets found to be differentially expressed between low and normal mHLA-DR patients on day 3 in the discovery phase were also expressed and able to differentiate septic shock patients in an independent and unstratified cohort subset (MIP-Rea, validation phase). The 193 identified HERV/MaLR probesets were tested on the MIP-Rea cohort, and patients were analyzed by unsupervised hierarchical clustering, according to the expression of these probesets.

Two well-separated groups of septic shock patients were highlighted (Table 2 and Additional file 8: Table S3). The first group, hereafter called cluster 1, consisted of 60 septic shock patients and the second group, called cluster 2, included 40 septic shock patients. Both groups presented significantly different SOFA severity scores on day 3 and day 6 (P value 0.003 and P value < 0.001, respectively), and all hematology criteria measured on day 3 (white cells, neutrophils, lymphocytes, platelets) were also found to be significantly different (P value < 0.001) (Table 2 and extended Additional file 8: Table S3). Even when outcomes were not statistically different, cluster 1 was characterized by 27% of patients presenting with an HAI and 27% that did not survive. Cluster 2 included a similar proportion of patients presenting with an HAI (25%) but with lower mortality (18%). Lastly, when we examined the expression of molecular markers, all were significantly different between the two clusters. The markers were HLA-DRA [36], S100A9 [37], CX3CR1 [38] expression on day 3, and CD74 ratio (D3/D1) [28], which are all linked to the immune system and were previously shown to be prognostic biomarkers (with either mortality or secondary infections) (Table 2 and extended Additional file 8: Table S3). Severity criteria, outcomes, and molecular markers confirmed that cluster 1 contains more severely affected patients.

These observations indicated that the HERV/MaLR probesets identified in the mHLA-DR-stratified patients (discovery phase) helped to differentiate an independent and unselected septic shock cohort (validation phase), based on the septic shock severity, outcomes, and prognostic molecular biomarkers.

Discussion

Over the past 30 years, there has been a considerable decline in the early mortality rate during the first days after septic shock. However, 50 to 70% of the mortality still happens later, after the first week in ICU. As this may reflect a persistent immunosuppressive state [3], there is a need to identify patients who would benefit from immunostimulatory therapies [9, 10]. Although flow cytometric measurement of mHLA-DR is known to be a reliable biomarker associated with death and secondary infections in septic patients [24, 26, 39], the test format makes it difficult to implement in large, multi-centered studies and outside specialized immunology laboratories [27]. Conversely, molecular markers may be easier to introduce, as illustrated by the use of automated tests, with standardized methodologies, for pathogen detection [40]. Several transcriptomic studies have obtained promising candidate signatures to stratify septic patients [41,42,43]. This report suggests that additional loci, identified outside the limited scope of the exome, may help to allow an optimal stratification of patients, taking into account the inter-individual variability and the immune status in the first days after ICU admission.

Firstly, we provided an initial view of the HERV/MaLR transcriptome landscape in the whole blood of HV and septic shock patients. We highlighted that this transcriptome was modified between the two groups. Secondly, the HERV/MaLR transcriptome was differentially expressed between patients, according to their mHLA-DR expression, and we were able to stratify an independent cohort using the identified HERV/MaLR probesets, with a clear difference in severity between the two clusters of patients.

A global view of the HERV transcriptome in the whole blood was obtained using the HERV-V3 chip. We observed that approximately 6.9% of LTR retrotransposons were transcriptionally active, which is similar to the 5.6% we had previously observed, ex vivo, in a PBMCs/endotoxin tolerance model [13]. Of note, 82.4% of the expressed probesets (58,587 probesets) were shared between the whole blood and PBMCs. The differences observed between whole blood and PBMCs may be due to (i) the variability among whole blood samples, (ii) the cell type composition, such as the presence of neutrophils in the whole blood, and (iii) the stimuli released from the endothelial environment. We observed a high proportion of gamma-retroviruses, which included HERV families that potentially code for envelope proteins that comprise an immunosuppressive domain (ISD) [44]. These families included HERV-H, HERV-W, HERV-FRD, HERV-Fc2, and HERV-T [45]. Expression of the two former families was previously observed in PBMCs of healthy subjects [15, 45,46,47,48] and on B cells and monocytes of multiple sclerosis patients [49,50,51]. In addition, we detected the expression of poorly characterized groups, such as MER52A and PRIMA41. The promoter function of MER52A was demonstrated in HEK-293T kidney cells [52] and appears to be negatively controlled by epigenetic factors, such as methylation in stem cells or PBMCs [52]. Similar epigenetic alterations are also seen in sepsis [53, 54]. The LTRC/D of MER41 is a member of the PRIMA41 elements that have been identified as enhancers for adjacent IFN-induced genes and have been shown to be involved in the regulation of essential immune functions [21]. Overall, differential expression analysis between HV and septic shock patients confirmed that septic shock is associated with an early modulation of thousands of HERV/MaLR loci, which is in agreement with (1) in vitro LTR retrotransposons being tightly regulated to endotoxin-induced stress [13] and (2) in vivo preliminary observations in critically ill patients [12].

We showed that individual HERV/MaLR can be differentially expressed between patients, according to mHLA-DR as a proxy of the immune status. We found that 193 probesets, which represented 168 distinct HERV/MaLR loci, were differentially expressed in septic shock patients, depending on their immunosuppression status at day 3. The potential for stratification using these HERV/MaLR loci was validated on an independent cohort, with the identification of a severely affected group that had higher SOFA scores, decreased CX3CR1 expression (molecular marker associated with survival) [38], decreased CD74 and mHLA-DR expression (associated with HAI) [28, 55], and increased S100A9 expression (associated with increased risk of secondary infections and mortality) [56]. As this pilot study was designed as a proof of concept study to show that HERV transcriptome can reflect the immune status, larger studies that involve higher numbers of patients are now warranted. Some of the identified immunosuppression HERV/MaLR loci will thus be evaluated in the REALISM study [57], which aims to determine the incidence, severity, and persistence of innate and adaptive immune alterations in 550 ICU patients (sepsis, severe trauma/burn, and major surgery), compared to 180 age-matched HV. We will study whether HERV/MaLR non-conventional genes can provide added value to conventional genes [28, 38, 41, 42], in regard to the characterization of their immune status and risk prognostication.

Finally, 10% of loci in the 193 HERV/MaLR probesets, identified between septic shock patients with a low vs normal mHLA-DR expression (discovery phase), belonged to the well-annotated region of the HERV/MaLR subpart of the genome (e.g., HERV-H, HERV-L, HERV-E4.1, and MER41). Forty-eight percent and 42% corresponded to the roughly annotated HERV and MaLR repertoires, respectively. The loci were mainly found in LTRs (95%), with internal gag pol env genes being underrepresented (5%) in the whole blood transcriptome. Eighty-five percent of the differentially expressed HERV/MaLR loci were located within or close (≤ 1 kb) to genes, which is significantly higher than the 41% of intragenic expressed LTRs in the PBMC-related endotoxin tolerance model [13] and the 12% of intragenic expressed LTRs in cancer tissues [35]. This bias may reflect either distinct levels of transcriptional control or some putative functional implications. In addition, intragenic LTR orientation was biased towards the antisense representation, i.e., 59% of the 193 identified HERV/MaLR loci were opposed to the transcriptional orientation of the conventional gene. A similar result (65%) was found in PBMCs expressed in the endotoxin tolerance model [13]. This may reflect the capacity of the antisense HERVs to generate dsRNAs [58] that trigger cytosolic sensing and modulate innate immunity, as found in cancer, following IFN-γ exposure [20], or when using methyltransferase inhibitors that lead to a type I interferon response [19]. Conversely, 25% of the 193 identified HERV/MaLR loci presented a similar transcriptional orientation to the associated conventional gene. Although this remains to be demonstrated, they may exert cis-functions, such as promoter or polyadenylation for immunity genes, as observed in oncogenes [59, 60], with 8% and 24% being predicted as a promoter or polyA site, respectively. They may also provide trans-functions, such as being an enhancer reservoir for immunity genes, as recently illustrated by the activation of the AIM2 inflammasome through the transcription factor STAT1 binding to the MER41 loci [21]. Of note, some MER41 solo LTRs were identified and STAT1 is one of the three genes (STAT1, CCR4, and HLA-DRB1/B3) identified in this study as differentially expressed between mHLA-DR-stratified patients but not between patients and HV. Lastly, some loci may encode large proteins or small peptides that drive immunological response [44, 61,62,63], as exemplified by the 170369402HE41 loci that putatively encodes a truncated Env protein, which contains an ISD domain.

Our study had some limitations. Firstly, we used an experimental tool, the HERV-V3 chip, which was validated with the MAQC samples, as previously described [22]. Consequently, we deliberately chose to use a high false discovery rate (FDR at 10%) when we analyzed septic patients, as effectively used in cancer with the previous generation microarray [35, 64]. Nevertheless, the highest confidence should be in those loci that were validated by RT-qPCR. Further studies using RNA-seq could support our results by providing another unbiased view of the HERV transcriptome. Secondly, the unique stratification criteria used in the discovery cohort and the limited number of patients could be considered as restrictive. Even though mHLA-DR expression was not a gold standard for the prediction of immunosuppression status, no gold standard or clinical signs currently exist to accurately evaluate immunosuppression in the ICU. In addition, the decrease in mHLA-DR expression remains an independent predictor of bad outcomes after sepsis [7]. Indeed, HERV/MaLR identified under these conditions were able to define clusters of patients in a larger independent cohort, that were different in terms of severity criteria. To further validate these results in the ICU, the clinical relevance of some HERV/MaLR candidates, together with current conventional markers, will be evaluated using a large unbiased cohort for which the immune status will be objectively defined (REALISM project; NCT01931956). Basically, such association between HERV/MaLR expression and altered states of immunity does not evidence any direct mechanistic connection, which will require further investigations.

Conclusions

The development of immunosuppression shortly after sepsis is now well established, and there is a growing interest in immunostimulatory treatments for these patients. However, we are still lacking the appropriate tools to precisely select the patients who may respond well to this type of treatment. In this study, we used a microarray-based approach to unveil the expression of approximately 6.9% of LTR retrotransposons, which were putatively linked to the immune response, in the whole blood. We identified a panel of 162 HERV/MaLR loci that were differentially expressed in septic shock patients stratified by mHLA-DR expression. We showed that, on an independent and unselected cohort, these loci classified patients according to the severity and distinct immune profiles. The added value of these newly identified HERV/MaLR loci should now be evaluated in a larger cohort of septic patients. If they prove to be robust, they could further serve as a stratification tool prior to immunostimulatory treatment and to monitor drug efficacy, which could contribute to the reduction of mortality in sepsis patients. More generally, we have illustrated the importance of addressing both the exome and repetitive-DNA repertoires [65] to increase our understanding of sepsis pathophysiology.

Availability of data and materials

The microarray expression data were filed in the NCBI Gene Expression Omnibus and are accessible via GEO accession number GSE121352.

All other data related to the study are available upon request to the corresponding author. The sharing of the data was submitted for approval by the REALISM consortium steering committee, as per consortium agreement.

Notes

  1. 1.

    Repertoire: within HERV-V3 microarray, the probesets are grouped into four repertoire: (1) genes, (2) HERV prototype (good annotation), (3) HERV/MaLR Dfam (lower quality annotation), and (4) other probesets (such as lncRNA)

  2. 2.

    Class: each HERV group belongs to different classes: class I (Gamma- and Epsilon-like), II (Beta-like), and III (Spuma-like), according to their pol region and phylogeny similarity

  3. 3.

    Group: also improperly named as “families.” HERV loci belongs to HERV groups based on sequence (nucleotide and Pol amino acid), similarity, degree of fit, and taxonomic markers

  4. 4.

    LTR: repeat flanking sequence of HERV that allows the regulation of their expression may have a promoter or polyadenylation signal function

Abbreviations

AUC:

Area under the curve

D:

Day

DEG:

Differentially expressed genes

DEL:

Differentially expressed loci

DNA:

Desoxyribonucleic acid

FC :

Fold change

FDR:

False discovery rate

h:

Hours

HAI :

Hospital-acquired infection

HERV:

Human endogenous retrovirus

HV:

Healthy volunteer

ICU:

Intensive care unit

IFN:

Interferon

ISD:

Immunosuppressive domain

LTR:

Long terminal repeat

MaLR:

Mammalian apparent LTR retrotransposons

mHLA-DR:

Monocyte HLA-DR

PBMC:

Peripheral blood mononuclear cells

PCA:

Principal component analysis

REALISM:

REAnimation Low Immune Status Markers

RNA:

Ribonucleic acid

RT-qPCR:

Reverse transcription-quantitative polymerase Chain Reaction

SOFA:

Sequential Organ Failure Assessment Score

References

  1. 1.

    Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, Bellomo R, Bernard GR, Chiche JD, Coopersmith CM, et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. 2016;315(8):801–10.

  2. 2.

    Marshall JC. Why have clinical trials in sepsis failed? Trends Mol Med. 2014;20(4):195–203.

  3. 3.

    Hotchkiss RS, Monneret G, Payen D. Sepsis-induced immunosuppression: from cellular dysfunctions to immunotherapy. Nat Rev Immunol. 2013;13(12):862–74.

  4. 4.

    Leentjens J, Kox M, Koch RM, Preijers F, Joosten LA, van der Hoeven JG, Netea MG, Pickkers P. Reversal of immunoparalysis in humans in vivo: a double-blind, placebo-controlled, randomized pilot study. Am J Respir Crit Care Med. 2012;186(9):838–45.

  5. 5.

    Monneret G, Finck ME, Venet F, Debard AL, Bohe J, Bienvenu J, Lepape A. The anti-inflammatory response dominates after septic shock: association of low monocyte HLA-DR expression and high interleukin-10 concentration. Immunol Lett. 2004;95(2):193–8.

  6. 6.

    Monneret G, Venet F. Monocyte HLA-DR in sepsis: shall we stop following the flow? Crit Care. 2014;18(1):102.

  7. 7.

    Tamulyte S, Kopplin J, Brenner T, Weigand MA, Uhle F. Monocyte HLA-DR assessment by a novel point-of-care device is feasible for early identification of ICU patients with complicated courses-a proof-of-principle study. Front Immunol. 2019;10:432.

  8. 8.

    Venet F, Monneret G. Advances in the understanding and treatment of sepsis-induced immunosuppression. Nat Rev Nephrol. 2018;14(2):121–37.

  9. 9.

    Meisel C, Schefold JC, Pschowski R, Baumann T, Hetzger K, Gregor J, Weber-Carstens S, Hasper D, Keh D, Zuckermann H, et al. Granulocyte-macrophage colony-stimulating factor to reverse sepsis-associated immunosuppression: a double-blind, randomized, placebo-controlled multicenter trial. Am J Respir Crit Care Med. 2009;180(7):640–8.

  10. 10.

    Delsing CE, Gresnigt MS, Leentjens J, Preijers F, Frager FA, Kox M, Monneret G, Venet F, Bleeker-Rovers CP, van de Veerdonk FL, et al. Interferon-gamma as adjunctive immunotherapy for invasive fungal infections: a case series. BMC Infect Dis. 2014;14:166.

  11. 11.

    Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al. Initial sequencing and analysis of the human genome. Nature. 2001;409(6822):860–921.

  12. 12.

    Tabone O, Mommert M, Jourdan C, Cerrato E, Legrand M, Lepape A, Allaouchiche B, Rimmele T, Pachot A, Monneret G, et al. Endogenous retroviruses transcriptional modulation after severe infection, Trauma and Burn. Front Immunol. 2018;9:3091.

  13. 13.

    Mommert M, Tabone O, Oriol G, Cerrato E, Guichard A, Naville M, Fournier P, Volff JN, Pachot A, Monneret G, et al. LTR-retrotransposon transcriptome modulation in response to endotoxin-induced stress in PBMCs. BMC Genomics. 2018;19(1):522.

  14. 14.

    Lee KH, Rah H, Green T, Lee YK, Lim D, Nemzek J, Wahl W, Greenhalgh D, Cho K. Divergent and dynamic activity of endogenous retroviruses in burn patients and their inflammatory potential. Exp Mol Pathol. 2014;96(2):178–87.

  15. 15.

    Johnston JB, Silva C, Holden J, Warren KG, Clark AW, Power C. Monocyte activation and differentiation augment human endogenous retrovirus expression: implications for inflammatory brain diseases. Ann Neurol. 2001;50(4):434–42.

  16. 16.

    Lee YJ, Jeong BH, Park JB, Kwon HJ, Kim YS, Kwak IS. The prevalence of human endogenous retroviruses in the plasma of major burn patients. Burns. 2013;39(6):1200–5.

  17. 17.

    Andersson G, Svensson AC, Setterblad N, Rask L. Retroelements in the human MHC class II region. Trends Genet. 1998;14(3):109–14.

  18. 18.

    Hurst TP, Magiorkinis G. Epigenetic control of human endogenous retrovirus expression: focus on regulation of long-terminal repeats (LTRs). Viruses. 2017;9(6):130.

  19. 19.

    Chiappinelli KB, Strissel PL, Desrichard A, Li H, Henke C, Akman B, Hein A, Rote NS, Cope LM, Snyder A, et al. Inhibiting DNA methylation causes an interferon response in cancer via dsRNA including endogenous retroviruses. Cell. 2015;162(5):974–86.

  20. 20.

    Canadas I, Thummalapalli R, Kim JW, Kitajima S, Jenkins RW, Christensen CL, Campisi M, Kuang Y, Zhang Y, Gjini E, et al. Tumor innate immunity primed by specific interferon-stimulated endogenous retroviruses. Nat Med. 2018;24(8):1143–50.

  21. 21.

    Chuong EB, Elde NC, Feschotte C. Regulatory evolution of innate immunity through co-option of endogenous retroviruses. Science. 2016;351(6277):1083–7.

  22. 22.

    Becker J, Perot P, Cheynet V, Oriol G, Mugnier N, Mommert M, Tabone O, Textoris J, Veyrieras JB, Mallet F. A comprehensive hybridization model allows whole HERV transcriptome profiling using high density microarray. BMC Genomics. 2017;18(1):286.

  23. 23.

    Landelle C, Lepape A, Francais A, Tognet E, Thizy H, Voirin N, Timsit JF, Monneret G, Vanhems P. Nosocomial infection after septic shock among intensive care unit patients. Infect Control Hosp Epidemiol. 2008;29(11):1054–65.

  24. 24.

    Landelle C, Lepape A, Voirin N, Tognet E, Venet F, Bohe J, Vanhems P, Monneret G. Low monocyte human leukocyte antigen-DR is independently associated with nosocomial infections after septic shock. Intensive Care Med. 2010;36(11):1859–66.

  25. 25.

    Cazalis MA, Friggeri A, Cave L, Demaret J, Barbalat V, Cerrato E, Lepape A, Pachot A, Monneret G, Venet F. Decreased HLA-DR antigen-associated invariant chain (CD74) mRNA expression predicts mortality after septic shock. Crit Care. 2013;17(6):R287.

  26. 26.

    Monneret G, Lepape A, Voirin N, Bohe J, Venet F, Debard AL, Thizy H, Bienvenu J, Gueyffier F, Vanhems P. Persisting low monocyte human leukocyte antigen-DR expression predicts mortality in septic shock. Intensive Care Med. 2006;32(8):1175–83.

  27. 27.

    Monneret G, Venet F, Meisel C, Schefold JC. Assessment of monocytic HLA-DR expression in ICU patients: analytical issues for multicentric flow cytometry studies. Crit Care. 2010;14(4):432.

  28. 28.

    Peronnet E, Venet F, Maucort-Boulch D, Friggeri A, Cour M, Argaud L, Allaouchiche B, Floccard B, Aubrun F, Rimmele T, et al. Association between mRNA expression of CD74 and IL10 and risk of ICU-acquired infections: a multicenter cohort study. Intensive Care Med. 2017;43(7):1013–20.

  29. 29.

    Huber W, Carey VJ, Gentleman R, Anders S. Orchestrating high-throughput genomic analysis with Bioconductor. Nature Methods. 2015;12(2):115–21.

  30. 30.

    Wilson CL, Miller CJ. Simpleaffy: a BioConductor package for Affymetrix Quality Control and data analysis. Bioinformatics. 2005;21(18):3683–5.

  31. 31.

    Kauffmann A, Gentleman R, Huber W. arrayQualityMetrics--a BioConductor package for quality assessment of microarray data. Bioinformatics. 2009;25(3):415–6.

  32. 32.

    Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics (Oxford). 2007;8(1):118–27.

  33. 33.

    Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.

  34. 34.

    YBaY H. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B (Methodological). 1995;57(1):289–300.

  35. 35.

    Perot P, Mugnier N, Montgiraud C, Gimenez J, Jaillard M, Bonnaud B, Mallet F. Microarray-based sketches of the HERV transcriptome landscape. PLoS One. 2012;7(6):e40194.

  36. 36.

    Xu J, Li J, Xiao K, Zou S, Yan P, Xie X, Xie L. Dynamic changes in human HLA-DRA gene expression and Th cell subsets in sepsis: indications of immunosuppression and associated outcomes. Scand J Immunol. 2019;91(1):e12813.

  37. 37.

    Uhel F, Azzaoui I, Gregoire M, Pangault C, Dulong J, Tadie JM, Gacouin A, Camus C, Cynober L, Fest T, et al. Early expansion of circulating granulocytic myeloid-derived suppressor cells predicts development of nosocomial infections in patients with sepsis. Am J Respir Crit Care Med. 2017;196(3):315–27.

  38. 38.

    Friggeri A, Cazalis MA, Pachot A, Cour M, Argaud L, Allaouchiche B, Floccard B, Schmitt Z, Martin O, Rimmele T, et al. Decreased CX3CR1 messenger RNA expression is an independent molecular biomarker of early and late mortality in critically ill patients. Crit Care. 2016;20(1):204.

  39. 39.

    Pachot A, Monneret G, Brion A, Venet F, Bohe J, Bienvenu J, Mougin B, Lepape A. Messenger RNA expression of major histocompatibility complex class II genes in whole blood from septic shock patients. Crit Care Med. 2005;33(1):31–8 discussion 236-237.

  40. 40.

    Meyers L, Ginocchio CC, Faucett AN, Nolte FS, Gesteland PH, Leber A, Janowiak D, Donovan V, Dien Bard J, Spitzer S, et al. Automated real-time collection of pathogen-specific diagnostic data: syndromic infectious disease epidemiology. JMIR Public Health Surveill. 2018;4(3):e59.

  41. 41.

    Davenport EE, Burnham KL, Radhakrishnan J, Humburg P, Hutton P, Mills TC, Rautanen A, Gordon AC, Garrard C, Hill AV, et al. Genomic landscape of the individual host response and outcomes in sepsis: a prospective cohort study. Lancet Respir Med. 2016;4(4):259–71.

  42. 42.

    Scicluna BP, van Vught LA, Zwinderman AH, Wiewel MA, Davenport EE, Burnham KL, Nurnberg P, Schultz MJ, Horn J, Cremer OL, et al. Classification of patients with sepsis according to blood genomic endotype: a prospective cohort study. Lancet Respir Med. 2017;5(10):816–26.

  43. 43.

    Antcliffe DB, Burnham KL, Al-Beidh F, Santhakumaran S, Brett SJ, Hinds CJ, Ashby D, Knight JC, Gordon AC. Transcriptomic signatures in sepsis and a differential response to steroids. From the VANISH randomized trial. Am J Respir Crit Care Med. 2019;199(8):980–6.

  44. 44.

    Bolze PA, Mommert M, Mallet F. Contribution of syncytins and other endogenous retroviral envelopes to human placenta pathologies. Prog Mol Biol Transl Sci. 2017;145:111–62.

  45. 45.

    de Parseval N, Lazar V, Casella JF, Benit L, Heidmann T. Survey of human genes of retroviral origin: identification and transcriptome of the genes with coding capacity for complete envelope proteins. J Virol. 2003;77(19):10414–22.

  46. 46.

    Balestrieri E, Pica F, Matteucci C, Zenobi R, Sorrentino R, Argaw-Denboba A, Cipriani C, Bucci I, Sinibaldi-Vallebona P. Transcriptional activity of human endogenous retroviruses in human peripheral blood mononuclear cells. Biomed Res Int. 2015;2015:164529.

  47. 47.

    Balestrieri E, Cipriani C, Matteucci C, Capodicasa N, Pilika A, Korca I, Sorrentino R, Argaw-Denboba A, Bucci I, Miele MT, et al. Transcriptional activity of human endogenous retrovirus in Albanian children with autism spectrum disorders. New Microbiol. 2016;39(3):228–31.

  48. 48.

    Kowalczyk MJ, Danczak-Pazdrowska A, Szramka-Pawlak B, Zaba R, Osmola-Mankowska A, Silny W. Human endogenous retroviruses and chosen disease parameters in morphea. Postep Dermatol Alergologii. 2017;34(1):47–51.

  49. 49.

    Antony JM, Deslauriers AM, Bhat RK, Ellestad KK, Power C. Human endogenous retroviruses and multiple sclerosis: innocent bystanders or disease determinants? Biochim Biophys Acta. 2011;1812(2):162–76.

  50. 50.

    Brudek T, Christensen T, Aagaard L, Petersen T, Hansen HJ, Moller-Larsen A. B cells and monocytes from patients with active multiple sclerosis exhibit increased surface expression of both HERV-H Env and HERV-W Env, accompanied by increased seroreactivity. Retrovirology. 2009;6:104.

  51. 51.

    Mameli G, Astone V, Arru G, Marconi S, Lovato L, Serra C, Sotgiu S, Bonetti B, Dolei A. Brains and peripheral blood mononuclear cells of multiple sclerosis (MS) patients hyperexpress MS-associated retrovirus/HERV-W endogenous retrovirus, but not human herpesvirus 6. J General Virol. 2007;88(Pt 1):264–74.

  52. 52.

    Zhang B, Xing X, Li J, Lowdon RF, Zhou Y, Lin N, Zhang B, Sundaram V, Chiappinelli KB, Hagemann IS, et al. Comparative DNA methylome analysis of endometrial carcinoma reveals complex and distinct deregulation of cancer promoters and enhancers. BMC Genomics. 2014;15:868.

  53. 53.

    Saeed S, Quintin J, Kerstens HH, Rao NA, Aghajanirefah A, Matarese F, Cheng SC, Ratter J, Berentsen K, van der Ent MA, et al. Epigenetic programming of monocyte-to-macrophage differentiation and trained innate immunity. Science. 2014;345(6204):1251086.

  54. 54.

    Hassan FI, Didari T, Khan F, Mojtahedzadeh M, Abdollahi M. The role of epigenetic alterations involved in sepsis: an overview. Curr Pharm Des. 2018;24(24):2862–9.

  55. 55.

    Cajander S, Tina E, Backman A, Magnuson A, Stralin K, Soderquist B, Kallman J. Quantitative real-time polymerase chain reaction measurement of HLA-DRA gene expression in whole blood is highly reproducible and shows changes that reflect dynamic shifts in monocyte surface HLA-DR expression during the course of sepsis. PLoS One. 2016;11(5):e0154690.

  56. 56.

    Fontaine M, Planel S, Peronnet E, Turrel-Davin F, Piriou V, Pachot A, Monneret G, Lepape A, Venet F. S100A8/A9 mRNA induction in an ex vivo model of endotoxin tolerance: roles of IL-10 and IFNgamma. PLoS One. 2014;9(6):e100909.

  57. 57.

    Rol ML, Venet F, Rimmele T, Moucadel V, Cortez P, Quemeneur L, Gardiner D, Griffiths A, Pachot A, Textoris J, et al. The REAnimation low immune status markers (REALISM) project: a protocol for broad characterisation and follow-up of injury-induced immunosuppression in intensive care unit (ICU) critically ill patients. BMJ Open. 2017;7(6):e015734.

  58. 58.

    Gogvadze E, Stukacheva E, Buzdin A, Sverdlov E. Human-specific modulation of transcriptional activity provided by endogenous retroviral insertions. J Virol. 2009;83(12):6098–105.

  59. 59.

    Lamprecht B, Walter K, Kreher S, Kumar R, Hummel M, Lenze D, Kochert K, Bouhlel MA, Richter J, Soler E, et al. Derepression of an endogenous long terminal repeat activates the CSF1R proto-oncogene in human lymphoma. Nat Med. 2010;16(5):571–9 571p following 579.

  60. 60.

    Babaian A, Mager DL. Endogenous retroviral promoter exaptation in human cancer. Mob DNA. 2016;7:24.

  61. 61.

    Cherkasova E, Scrivani C, Doh S, Weisman Q, Takahashi Y, Harashima N, Yokoyama H, Srinivasan R, Linehan WM, Lerman MI, et al. Detection of an immunogenic HERV-E envelope with selective expression in clear cell kidney cancer. Cancer Res. 2016;76(8):2177–85.

  62. 62.

    Schiavetti F, Thonnard J Fau - Colau D, Colau D Fau - Boon T, Boon T Fau - Coulie PG, Coulie PG: A human endogenous retroviral sequence encoding an antigen recognized on melanoma by cytolytic T lymphocytes. 2002(0008-5472 (Print)).

  63. 63.

    Mullins CS, Linnebacher M. Human endogenous retroviruses and cancer: causality and therapeutic possibilities. World J Gastroenterol. 2012;18(42):6027–35.

  64. 64.

    Perot P, Mullins CS, Naville M, Bressan C, Huhns M, Gock M, Kuhn F, Volff JN, Trillet-Lenoir V, Linnebacher M, et al. Expression of young HERV-H loci in the course of colorectal carcinoma and correlation with molecular subtypes. Oncotarget. 2015;6(37):40095–111.

  65. 65.

    Slotkin RK. The case for not masking away repetitive DNA. Mob DNA. 2018;9:15.

Download references

Acknowledgements

Collaborators of the REALISM Study Group: For Hospices Civils de Lyon: André Boibieux, Julien Davidson, Laure Fayolle-PivoT, Julie Gatel, Charline Genin, Arnaud Gregoire, Alain Lepape, Anne-Claire Lukaszewicz, Guillaume Marcotte, Marie Matray, Delphine Maucort-Boulch, Guillaume Monneret, Nathalie Panel, Thomas Rimmele, Hélène Vallin, Fabienne Venet. For bioMérieux: Sophie Blein, Karen Brengel-Pesce, Elisabeth Cerrato, Valérie Cheynet, Emmanuelle Gallet-Gorius, Audrey Guichard, Camille Jourdan, Natacha Koenig, François Mallet, Boris Meunier, Virginie Moucadel, Marine Mommert, Guy Oriol, Alexandre Pachot, Claire Schrevel, Olivier Tabone, Julien Textoris, Javier Yugueros Marcos. For BIOASTER: Jérémie Becker, Frédéric Bequet, Yacine Bounab, Florian Brajon, Bertrand Canard, Muriel Collus, Nathalie Garcon, Irène Gorse, Cyril Guyard, Fabien Lavocat, Philippe Leissner, Karen Louis, Maxime Mistretta, Jeanne Moriniere, Yoann MouscaZ, Laura Noailles, Magali Perret, Frédéric Reynier, Cindy Riffaud, Mary-Luz Rol, Nicolas Sapay, Trang Tran, Christophe Vedrine. For Sanofi: Christophe Carre, Pierre Cortez, Aymeric De Monfort, Karine Florin, Laurent Fraisse, Isabelle Fugier, Maïna L’Azou, Sandrine Payrard, Annick PELERAUX, Laurence Quemeneur. For ESPCI Paris: Andrew Griffiths, Stephanie Toetsch. For GSK: Teri Ashton, Peter J. Gough, Scott B. Berger, David Gardiner, Iain Gillespie, Aidan Macnamara, Aparna Raychaudhuri, Rob Smylie, Lionel Tan, Craig Tipple.

Collaborators of the MIPrea Study Group: Christine Alberti-Segui, Bernard Allaouchiche, Laurent Argaud, Frédéric Aubrun, Véronique Barbalat, Thomas Baudry, Julien Bohé, Marie-Angélique Cazalis, Elisabeth Cerrato, Martin Cour, Sylvie De La Salle, Benjamin Delwarde, Bernard Floccard, Arnaud Friggeri, Emmanuelle Gallet-Gorius, Christian Guillaume, Romain Hernu, Audrey Larue-Triolet, Alain Lepape, Guillaume Marcotte, Delphine Maucort-Boulch, Boris Meunier, Guillaume Monneret, Stéphane Morisset, Julie Mouillaux, Alexandre Pachot, Mathieu Page, Nathalie Panel, Estelle Peronnet, Vincent Piriou, Anne Portier, Marion Provent, Thomas Rimmelé, Marie Simon, Julien Textoris, Fabrice Thiolliere, Hélène Vallin, Fabienne Venet.

Funding

The project was funded by a consortium: bioMérieux, SANOFI, GlaxoSmithKline, Ecole Supérieure de Physique Chimie Industrielles de la Ville the Paris – PSL Research University, the University Hospital Hospices Civils de Lyon and the microbiology technological institute BIOASTER. In addition, MM and OT were supported by doctoral grants from bioMérieux and OT was supported by the Association Nationale de la Recherche et de la Technologie (ANRT). bioMérieux and Lyon University Hospital were co-responsible for the study design, collection of data, interpretation of results, decision to publish, and writing of the report. Other partners had no active role in the study.

Author information

All authors were involved in the analysis and interpretation of data as well as drafting the manuscript or revising it critically for important intellectual content. FM, GM, FV, JT, AP, AL, MD, and KBP made substantial contributions to the conception and design of the study. MM, GO, CV, and EC designed and performed the microarray experiments. OT and MM performed the data analyses. MM, OT, JT, and FM performed data interpretations. MM and AG performed the primer pair design, screening and validation, RT-qPCR experiments, and statistical analysis. MM, OT, JT, and FM wrote the paper. All authors read and approved the final manuscript.

Correspondence to Marine Mommert.

Ethics declarations

Ethics approval and consent to participate

Immunosepsis cohort subset

Patients were included from a cohort of adult septic shock patients enrolled from December 2001 to April 2005 in two ICUs from a French university hospital. The clinical study was approved by the regional ethics committee (Comité Consultatif de Protection des Personnes dans la Recherche BioMédicale de Lyon A), which waived the need for written informed consent because the study was observational, with a low risk to the patients, and no specific procedure, other than routine blood sampling, was required. Oral information and non-opposition to inclusion in the study were mandatory and were systematically obtained before any blood sample was drawn. This was recorded in patients’ clinical files.

MIP-Rea cohort subset

The study was approved by our institutional ethical review board (Comité d’Ethique des Centres d’investigation Clinique de l’Inter-Région Rhône-Alpes Auvergne – IRB 5044), and consent for an ancillary study was obtained a posteriori. The protocol of this ancillary retrospective study was submitted to the French CCTIRS and CNIL committees and approved on April 22, 2016, and September 30, 2016, respectively. Further, informed consent was received from patients for inclusion in this specific study.

Consent for publication

Not applicable

Competing interests

MM, AG, GO, EC, VC, AP, KBP, JT, and FM are employees of bioMérieux, an in vitro diagnostic company. All other authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1 : Table S1. The HERV transcriptome in whole blood samples (descriptive phase).

Additional file 2 : Supplementary Methods.

Additional file 3 : Figure S1. Description of HERV/MaLR transcriptome in whole blood of septic shock patients from the MIP-Rea cohort subset (descriptive phase).

Additional file 4 : Figure S2. Genomic, transcriptomic and functional projections of the “HERV_prototype” repertoire in the IMMUNOSEPSIS cohort subset (descriptive phase).

Additional file 5 : Table S2. Differential expression of HERV-V3 chip repertoires induced by immune state changes in whole blood samples (Healthy volunteers and IMMUNOSEPSIS cohort subset, discovery phase).

Additional file 6 : Figure S3. Selection, design and quality criteria for the locus-specific qPCR systems, illustrated with the 060400302-HERV0489uL locus.

Additional file 7 : Figure S4. RT-qPCR validation of 29 HERV/MaLR loci that were identified by microarray and were differentially expressed, according to themHLA-DR expression (discovery phase).

Additional file 8 : Table S3. MIP-Rea patient characteristics according to cluster with the 193 HERV/MaLR that were differentially expressed in the IMMUNOSEPSIS cohort subset on day 3, according to mHLA-DR expression (validation phase).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mommert, M., Tabone, O., Guichard, A. et al. Dynamic LTR retrotransposon transcriptome landscape in septic shock patients. Crit Care 24, 96 (2020). https://doi.org/10.1186/s13054-020-2788-8

Download citation

Keywords

  • HERV transcriptome
  • Whole blood
  • Septic shock patients
  • mHLA-DR expression
  • Immunosuppression