Human Endogenous Retrovirus (HERV) Transcriptome Is Dynamically Modulated during SARS-CoV-2 Infection and Allows Discrimination of COVID-19 Clinical Stages

SARS-CoV-2 infection is known to trigger an important inflammatory response, which has a major role in COVID-19 pathogenesis. In infectious and inflammatory contexts, the modulation of human endogenous retroviruses (HERV) has been broadly reported, being able to further sustain innate immune responses due to the expression of immunogenic viral transcripts, including double-stranded DNA (dsRNA), and eventually, immunogenic proteins. To gain insights on this poorly characterized interplay, we performed a high-throughput expression analysis of similar to 3,300 specific HERV loci in the peripheral blood mononuclear cells (PBMCs) of 10 healthy controls and 16 individuals being either convalescent after the infection (6) or retesting positive after convalescence (10) (Gene Expression Onmibus [GEO] data set ). Results showed that the exposure to SARS-CoV-2 infection modulates HERV expression according to the disease stage and reflecting COVID-19 immune signatures. The differential expression analysis between healthy control (HC) and COVID-19 patients allowed us to identify a total of 282 differentially expressed HERV loci (deHERV) in the individuals exposed to SARS-CoV-2 infection, independently from the clinical form. In addition, 278 and 60 deHERV loci that were specifically modulated in individuals convalescent after COVID19 infection (C) and patients that retested positive to SARS-CoV-2 after convalescence (RTP) as individually compared to HC, respectively, as well as 164 deHERV loci between C and RTP patients were identified. The identified HERV loci belonged to 36 different HERV groups, including members of all three classes. The present study provides an exhaustive picture of the HERV transcriptome in PBMCs and its dynamic variation in the presence of COVID-19, revealing specific modulation patterns according to the infection stage that can be relevant to the disease clinical manifestation and outcome. IMPORTANCE We report here the first high-throughput analysis of HERV loci expression along SARS-CoV-2 infection, as performed with peripheral blood mononuclear cells (PBMCs). Such cells are not directly infected by the virus but have a crucial role in the plethora of inflammatory and immune events that constitute a major hallmark of COVID-19 pathogenesis. Results provide a novel and exhaustive picture of HERV expression in PBMCs, revealing specific modulation patterns according to the disease condition and the concomitant immune activation. To our knowledge, this is the first set of deHERVs whose expression is dynamically modulated across COVID-19 stages, confirming a tight interplay between HERV and cellular immunity and revealing specific transcriptional signatures that can have an impact on the disease clinical manifestation and outcome.