Epigenome-wide DNA methylation profiling of healthy COVID-19 recoverees reveals a unique signature in circulating immune cells

This article has been Reviewed by the following groups

Read the full article See related articles

Abstract

Background

Epigenetic alterations upon microbial challenge have been described as both a defence strategy and a result of pathogenic manipulation. While most COVID-19 studies focus on inflammatory and immune-mediated responses, little is known about epigenetic modifications in response to SARS-CoV-2 infection.

Methods

Epigenome-wide DNA methylation patterns from COVID-19 convalescents were compared to uninfected controls from before and after the pandemic. Peripheral blood mononuclear cell (PBMC) DNA was extracted from uninfected controls, COVID-19 convalescents and symptom-free individuals with SARS-CoV-2-specific T cell-responses, as well as from PBMCs stimulated in vitro with SARS-CoV-2. Subsequently, the Illumina MethylationEPIC 850K array was performed, and statistical/bioinformatic analyses comprised differential DNA methylation, pathway over-representation and module identification analyses.

Results

Differential DNA methylation patterns distinguished COVID-19 convalescents from uninfected controls, with similar results in an experimental SARS-CoV-2 infection model. A SARS-CoV-2-induced module was identified in vivo , comprising 66 genes of which six ( TP53, INS, HSPA4, SP1, ESR1 and FAS ) were present in corresponding in vitro analyses. Over-representation analyses revealed involvement in Wnt, muscarinic acetylcholine receptor signalling and gonadotropin-releasing hormone receptor pathways. Furthermore, numerous differentially methylated and network genes from both settings interacted with the SARS-CoV-2 interactome.

Conclusions

Altered DNA methylation patterns of COVID-19 convalescents suggest recovery from mild-to-moderate SARS-CoV-2 infection leaves longstanding epigenetic traces. As in vitro SARS-CoV-2 infection corroborated in vivo exposure results, this indicates DNA methylation is involved in immune cell responses to challenge with this virus. Future studies should determine whether this reflects host-induced protective antiviral defence or targeted viral hijacking to evade host defence.

Article activity feed

  1. SciScore for 10.1101/2021.07.05.21260014: (What is this?)

    Please note, not all rigor criteria are appropriate for all manuscripts.

    Table 1: Rigor

    EthicsConsent: For samples from the natural exposure cohort, all participants provided written informed consent, and the present study was approved by the Regional Ethics Committee for Human Research in Linköping (Dnr. 2019-0618).
    IRB: For samples from the natural exposure cohort, all participants provided written informed consent, and the present study was approved by the Regional Ethics Committee for Human Research in Linköping (Dnr. 2019-0618).
    Sex as a biological variablenot detected.
    Randomizationnot detected.
    Blindingnot detected.
    Power Analysisnot detected.
    Cell Line Authenticationnot detected.

    Table 2: Resources

    Antibodies
    SentencesResources
    Cons were defined as neither having any positive circulating IgG-antibody or T cell responses to SARS-CoV-2, while CC19s were defined by the presence of SARS-CoV-2-specific IgG antibodies in plasma using suspension multiplex immunoassay (SMIA), some of which were positive for IgG in saliva, rapid test and in T cell responses as well.
    SARS-CoV-2-specific IgG
    suggested: None
    Stimulation with anti-CD3 antibody at a concentration of 1 µg/ml was used as a positive control for each subject.
    anti-CD3
    suggested: None
    Anti-CD28 antibody (3608-1-50, Mabtech, Sweden) was included at a final concentration of 0.1 μg/ml as a co-stimulator.
    Anti-CD28
    suggested: None
    The beads were resuspended in 100 µl of 1 µg/ml goat anti-human IgG-PE labelled antibody (Southern BioTech,, Birmingham, AL, USA. Cat. #2040-09) in PBS-T + 1% BSA and incubated for 30 min at RT in the dark with rotation at 600 rpm.
    anti-human IgG-PE
    suggested: (SouthernBiotech Cat# 2040-09, RRID:AB_2795648)
    Next, the cells were incubated with mouse-anti-dsRNA antibody (Scions, Code: J2 at 1:100 dilution) for 1.5 h followed by detection using horseradish peroxidase–conjugated goat anti-mouse IgG (heavy plus light chain) (Catalog: 1706516, Bio-Rad Laboratories, Hercules, CA, USA) (1:1000) for 1 h.
    mouse-anti-dsRNA
    suggested: None
    anti-mouse IgG
    suggested: (Bio-Rad Cat# 170-6516, RRID:AB_11125547)
    Experimental Models: Cell Lines
    SentencesResources
    SARS-CoV-2 or mock Vero cell supernatant was added to the PBMC cultures corresponding to a multiplicity of infection of 0.01.
    Vero
    suggested: CLS Cat# 605372/p622_VERO, RRID:CVCL_0059)
    Finally, the SARS CoV-2 infected Vero E6 cells were identified using 3-aminoethylcarbazole (AEC) substrate.
    Vero E6
    suggested: None
    Software and Algorithms
    SentencesResources
    Next, the cells were incubated with mouse-anti-dsRNA antibody (Scions, Code: J2 at 1:100 dilution) for 1.5 h followed by detection using horseradish peroxidase–conjugated goat anti-mouse IgG (heavy plus light chain) (Catalog: 1706516, Bio-Rad Laboratories, Hercules, CA, USA) (1:1000) for 1 h.
    Bio-Rad Laboratories
    suggested: (Bio-Rad Laboratories, RRID:SCR_008426)
    Concentrations of extracted DNA were measured using the Qubit® 4.0 Fluorometer (Thermo Fisher Scientific, Waltham, Massachusetts, U.S), using dsDNA High Sensitivity (HS) Assay Kit and RNA HS Assay Kit.
    Qubit®
    suggested: None
    Thereafter, we performed singular value decomposition (SVD) analyses using the ChAMP package (45
    ChAMP
    suggested: (ChAMP, RRID:SCR_012891)
    Differential DNA methylation analysis in vivo: As we were interested in studying CpGs that were differentially methylated between CC19s and non-infected controls from both before and after the start of the COVID-19 pandemic, we performed differential DNA methylation analyses, using the limma package (version 3.46.0).
    limma
    suggested: (LIMMA, RRID:SCR_010943)
    Pathway over-representation analyses: To make biological sense of the putatively SARS-CoV-2-induced DNA methylation differences, we performed PANTHER pathway over-representation test analyses using the PANTHER database (version 16.0).
    PANTHER
    suggested: (PANTHER, RRID:SCR_004869)
    The significantly enriched pathways were displayed in dot plots generated in R using ggplot2 package (version 3.3.3).
    ggplot2
    suggested: (ggplot2, RRID:SCR_014601)
    The graph clustering algorithm MCODE (48) was used to identify molecular complexes and create a large disease module, which was then fitted to a protein-protein interaction network, and both were analysed and rendered in Cytoscape (version 3.8.0).
    Cytoscape
    suggested: (Cytoscape, RRID:SCR_003032)
    In addition, the inference of modules was performed with two other methods from the MODifieR package (DIAMOnD and WGCNA)(50) to study whether it was possible to condense the module genes to fewer genes of particular interest within the network, for both the in vivo and the in vitro setting.
    DIAMOnD
    suggested: (DIAMOND, RRID:SCR_009457)
    Overlap to SARS-CoV-2 interactome: A publicly available protein-protein interaction (PPI) network of SARS-COV-2 and human genes curated by BioGRID (version 4.4.197) was downloaded from the Network Data Exchange in Cytoscape (version 3.8.0).
    BioGRID
    suggested: (BioGrid Australia, RRID:SCR_006334)

    Results from OddPub: Thank you for sharing your data.


    Results from LimitationRecognizer: We detected the following sentences addressing limitations in the study:
    Although an obvious limitation of the study is the lack of validation of the DNAm findings on a transcriptional level, it serves as a pilot study that generates hypotheses for further studies within the field. Hence, whether the observed DNAm patterns are indeed associated or even causally linked to host protective or host detrimental immune responses still needs to be addressed in future studies. With more well-designed, larger, consecutive sample materials, possibly also in closer proximity to the time of infection with SARS-CoV-2, it will be possible to study the role of DNAm alterations in anti-viral defence and in viral manipulation of the same. An advantage of the investigation of epigenetic modifications in in mild to moderately ill patients, is that we may be able to discern DNAm differences that otherwise would have been masked due to an overriding inflammatory response. These subtle changes may not only be relevant to how a less severe immune response is mounted towards SARS-CoV-2, but also in the case of long-COVID-19. The presentation of longstanding symptoms could be caused by detrimentally changed DNAm patterns, originally triggered as a short-term anti-viral response. This should be explored in detail in further studies since the risk is that these short-term responses may permanently alter and erroneously manifest in the DNA methylome.

    Results from TrialIdentifier: No clinical trial numbers were referenced.


    Results from Barzooka: We did not find any issues relating to the usage of bar graphs.


    Results from JetFighter: We did not find any issues relating to colormaps.


    Results from rtransparent:
    • Thank you for including a conflict of interest statement. Authors are encouraged to include this statement when submitting to a journal.
    • Thank you for including a funding statement. Authors are encouraged to include this statement when submitting to a journal.
    • Thank you for including a protocol registration statement.

    Results from scite Reference Check: We found no unreliable references.


    About SciScore

    SciScore is an automated tool that is designed to assist expert reviewers by finding and presenting formulaic information scattered throughout a paper in a standard, easy to digest format. SciScore checks for the presence and correctness of RRIDs (research resource identifiers), and for rigor criteria such as sex and investigator blinding. For details on the theoretical underpinning of rigor criteria and the tools shown here, including references cited, please follow this link.