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
Listed in
- Evaluated articles (ScreenIT)
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
-
-
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
Ethics Consent: 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 variable not detected. Randomization not detected. Blinding not detected. Power Analysis not detected. Cell Line Authentication not detected. Table 2: Resources
Antibodies Sentences Resources Cons were defined as neither having any positive … 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
Ethics Consent: 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 variable not detected. Randomization not detected. Blinding not detected. Power Analysis not detected. Cell Line Authentication not detected. Table 2: Resources
Antibodies Sentences Resources 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 IgGsuggested: NoneStimulation with anti-CD3 antibody at a concentration of 1 µg/ml was used as a positive control for each subject. anti-CD3suggested: NoneAnti-CD28 antibody (3608-1-50, Mabtech, Sweden) was included at a final concentration of 0.1 μg/ml as a co-stimulator. Anti-CD28suggested: NoneThe 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-PEsuggested: (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-dsRNAsuggested: Noneanti-mouse IgGsuggested: (Bio-Rad Cat# 170-6516, RRID:AB_11125547)Experimental Models: Cell Lines Sentences Resources SARS-CoV-2 or mock Vero cell supernatant was added to the PBMC cultures corresponding to a multiplicity of infection of 0.01. Verosuggested: 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 E6suggested: NoneSoftware and Algorithms Sentences Resources 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 Laboratoriessuggested: (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: NoneThereafter, we performed singular value decomposition (SVD) analyses using the ChAMP package (45 ChAMPsuggested: (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). limmasuggested: (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). PANTHERsuggested: (PANTHER, RRID:SCR_004869)The significantly enriched pathways were displayed in dot plots generated in R using ggplot2 package (version 3.3.3). ggplot2suggested: (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). Cytoscapesuggested: (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. DIAMOnDsuggested: (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). BioGRIDsuggested: (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.
-