An epigenome-wide DNA methylation study of patients with COVID-19

Sarbecovirus disease susceptibility is conserved across viral and host models

Coronaviruses have caused three severe epidemics since the start of the 21st century: SARS, MERS and COVID-19. The severity of the ongoing COVID-19 pandemic and increasing likelihood of future coronavirus outbreaks motivates greater understanding of factors leading to severe coronavirus disease. We screened ten strains from the Collaborative Cross mouse genetic reference panel and identified strains CC006/TauUnc (CC006) and CC044/Unc (CC044) as coronavirus-susceptible and resistant, respectively, as indicated by variable weight loss and lung congestion scores four days post-infection. We generated a genetic mapping population of 755 CC006xCC044 F2 mice and exposed the mice to one of three genetically distinct mouse-adapted coronaviruses: clade 1a SARS-CoV MA15 (n=391), clade 1b SARS-CoV-2 MA10 (n=274), and clade 2 HKU3-CoV MA (n=90). Quantitative trait loci (QTL) mapping in SARS-CoV MA15- and SARS-CoV-2 MA10-infected F2 mice identified genetic loci associated with disease severity. Specifically, we identified seven loci associated with variation in outcome following infection with either virus, including one, HrS43, that is present in both groups. Three of these QTL, including HrS43, were also associated with HKU3-CoV MA outcome. HrS43 overlaps with a QTL previously reported by our lab that is associated with SARS-CoV MA15 outcome in CC011xCC074 F2 mice and is also syntenic with a human chromosomal region associated with severe COVID-19 outcomes in humans GWAS. The results reported here provide: (a) additional support for the involvement of this locus in SARS-CoV MA15 infection, (b) the first conclusive evidence that this locus is associated with susceptibility across the Sarbecovirus subgenus, and (c) demonstration of the relevance of mouse models in the study of coronavirus disease susceptibility in humans.

A methylation clock model of mild SARS-CoV-2 infection provides insight into immune dysregulation

DNA methylation comprises a cumulative record of lifetime exposures superimposed on genetically determined markers. Little is known about methylation dynamics in humans following an acute perturbation, such as infection. We characterized the temporal trajectory of blood epigenetic remodeling in 133 participants in a prospective study of young adults before, during, and after asymptomatic and mildly symptomatic SARS-CoV-2 infection. The differential methylation caused by asymptomatic or mildly symptomatic infections was indistinguishable. While differential gene expression largely returned to baseline levels after the virus became undetectable, some differentially methylated sites persisted for months of follow-up, with a pattern resembling autoimmune or inflammatory disease. We leveraged these responses to construct methylation-based machine learning models that distinguished samples from pre-, during-, and postinfection time periods, and quantitatively predicted the time since infection. The clinical trajectory in the young adults and in a diverse cohort with more severe outcomes was predicted by the similarity of methylation before or early after SARS-CoV-2 infection to the model-defined postinfection state. Unlike the phenomenon of trained immunity, the postacute SARS-CoV-2 epigenetic landscape we identify is antiprotective.

Role of epigenetics in the clinical evolution of COVID-19 disease. Epigenome-wide association study identifies markers of severe outcome

BACKGROUND

COVID-19 has a wide spectrum of clinical manifestations and given its impact on morbidity and mortality, there is an unmet medical need to discover endogenous cellular and molecular biomarkers that predict the expected clinical course of the disease. Recently, epigenetics and especially DNA methylation have been pointed out as a promising tool for outcome prediction in several diseases.

METHODS AND RESULTS

Using the Illumina Infinium Methylation EPIC BeadChip850K, we investigated genome-wide differences in DNA methylation in an Italian Cohort of patients with comorbidities and compared severe (n = 64) and mild (123) prognosis. Results showed that the epigenetic signature, already present at the time of Hospital admission, can significantly predict risk of severe outcomes. Further analyses provided evidence of an association between age acceleration and a severe prognosis after COVID-19 infection. The burden of Stochastic Epigenetic Mutation (SEMs) has been significantly increased in patients with poor prognosis. Results have been replicated in silico considering COVID-19 negative subjects and available previously published datasets.

CONCLUSIONS

Using original methylation data and taking advantage of already published datasets, we confirmed in the blood that epigenetics is actively involved in immune response after COVID-19 infection, allowing the identification of a specific signature able to discriminate the disease evolution. Furthermore, the study showed that epigenetic drift and age acceleration are associated with severe prognosis. All these findings prove that host epigenetics undergoes notable and specific rearrangements to respond to COVID-19 infection which can be used for a personalized, timely, and targeted management of COVID-19 patients during the first stages of hospitalization.

Dynamic transcriptome analyses reveal m6A regulated immune non-coding RNAs during dengue disease progression

Dengue infection is one of the most prevalent arthropod-borne viral diseases, which can result in severe complications. Identification of genes and long non-coding RNAs (lncRNAs) involved in dengue infection would help in deciphering potential mechanisms responsible for the disease progression. We comprehensively analyzed the dynamic transcriptome during dengue disease progression and identified critical genes and lncRNAs with expression perturbations. Our findings revealed that the expression of genes (i.e., CCR10 and GNG7) and lncRNAs (i.e., CTBP1-AS and MAFG-AS1) were potentially regulated by m6A RNA methylation. Interestingly, dengue viral proteins prevalently interact with genes or lncRNAs with expression perturbations, which are involved in cell cycle, inflammation signaling pathways and immune response. Dynamically expressed genes and lncRNAs were likely to locate in the central regions of human protein-protein network, which play crucial roles in mediating signaling spread and helping viral replication. Immune microenvironments analysis revealed that plasma cells levels were increased and T cells infiltrations were decreased during dengue disease progression. Dynamically expressed genes and lncRNAs were correlated with immune cell infiltrations. Moreover, network analysis reveals the associations between dengue viral infections and human complex diseases (i.e., digestive diseases and neoplasms). Our comprehensive transcriptome analysis of dengue disease progression identified potential gene and lncRNA biomarkers, providing novel insights for understanding the pathogenesis of and developing effective therapeutic strategies for dengue infection.

Single-cell multiomics revealed the dynamics of antigen presentation, immune response and T cell activation in the COVID-19 positive and recovered individuals

Introduction

Despite numerous efforts to describe COVID-19's immunological landscape, there is still a gap in our understanding of the virus's infections after-effects, especially in the recovered patients. This would be important to understand as we now have huge number of global populations infected by the SARS-CoV-2 as well as variables inclusive of VOCs, reinfections, and vaccination breakthroughs. Furthermore, single-cell transcriptome alone is often insufficient to understand the complex human host immune landscape underlying differential disease severity and clinical outcome.

Methods

By combining single-cell multi-omics (Whole Transcriptome Analysis plus Antibody-seq) and machine learning-based analysis, we aim to better understand the functional aspects of cellular and immunological heterogeneity in the COVID-19 positive, recovered and the healthy individuals.

Results

Based on single-cell transcriptome and surface marker study of 163,197 cells (124,726 cells after data QC) from the 33 individuals (healthy=4, COVID-19 positive=16, and COVID-19 recovered=13), we observed a reduced MHC Class-I-mediated antigen presentation and dysregulated MHC Class-II-mediated antigen presentation in the COVID-19 patients, with restoration of the process in the recovered individuals. B-cell maturation process was also impaired in the positive and the recovered individuals. Importantly, we discovered that a subset of the naive T-cells from the healthy individuals were absent from the recovered individuals, suggesting a post-infection inflammatory stage. Both COVID-19 positive patients and the recovered individuals exhibited a CD40-CD40LG-mediated inflammatory response in the monocytes and T-cell subsets. T-cells, NK-cells, and monocyte-mediated elevation of immunological, stress and antiviral responses were also seen in the COVID-19 positive and the recovered individuals, along with an abnormal T-cell activation, inflammatory response, and faster cellular transition of T cell subtypes in the COVID-19 patients. Importantly, above immune findings were used for a Bayesian network model, which significantly revealed FOS, CXCL8, IL1β, CST3, PSAP, CD45 and CD74 as COVID-19 severity predictors.

Discussion

In conclusion, COVID-19 recovered individuals exhibited a hyper-activated inflammatory response with the loss of B cell maturation, suggesting an impeded post-infection stage, necessitating further research to delineate the dynamic immune response associated with the COVID-19. To our knowledge this is first multi-omic study trying to understand the differential and dynamic immune response underlying the sample subtypes.

Potential links between COVID-19 and periodontitis: a bioinformatic analysis based on GEO datasets

Background

2019 Coronavirus disease (COVID-19) is an infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The COVID-19 pandemic has already had a serious influence on human existence, causing a huge public health concern for countries all around the world. Because SARS-CoV-2 infection can be spread by contact with the oral cavity, the link between oral illness and COVID-19 is gaining traction. Through bioinformatics approaches, we explored the possible molecular mechanisms linking the COVID-19 and periodontitis to provide the basis and direction for future research.

Methods

Transcriptomic data from blood samples of patients with COVID-19 and periodontitis was downloaded from the Gene Expression Omnibus database. The shared differentially expressed genes were identified. The analysis of Gene Ontology, Kyoto Encyclopedia of Genesand Genomes pathway, and protein–protein interaction network was conducted for the shared differentially expressed genes. Top 5 hub genes were selected through Maximal Clique Centrality algorithm. Then mRNA-miRNA network of the hub genes was established based on miRDB database, miRTarbase database and Targetscan database. The Least absolute shrinkage and selection operator regression analysis was used to discover possible biomarkers, which were then investigated in relation to immune-related genes.

Results

Fifty-six shared genes were identified through differential expression analysis in COVID-19 and periodontitis. The function of these genes was enriched in regulation of hormone secretion, regulation of secretion by cell. Myozenin 2 was identified through Least absolute shrinkage and selection operator regression Analysis, which was down-regulated in both COVID-19 and periodontitis. There was a positive correlation between Myozenin 2 and the biomarker of activated B cell, memory B cell, effector memory CD4 T cell, Type 17 helper cell, T follicular helper cell and Type 2 helper cell.

Conclusion

By bioinformatics analysis, Myozenin 2 is predicted to correlate to the pathogenesis and immune infiltrating of COVID-19 and periodontitis. However, more clinical and experimental researches are needed to validate the function of Myozenin 2.

HLA-C dysregulation as a possible mechanism of immune evasion in SARS-CoV-2 and other RNA-virus infections

One of the mechanisms by which viruses can evade the host’s immune system is to modify the host’s DNA methylation pattern. This work aims to investigate the DNA methylation and gene expression profile of COVID-19 patients, divided into symptomatic and asymptomatic, and healthy controls, focusing on genes involved in the immune response. In this study, changes in the methylome of COVID-19 patients’ upper airways cells, the first barrier against respiratory infections and the first cells presenting viral antigens, are shown for the first time. Our results showed alterations in the methylation pattern of genes encoding proteins implicated in the response against pathogens, in particular the HLA-C gene, also important for the T-cell mediated memory response. HLA-C expression significantly decreases in COVID-19 patients, especially in those with a more severe prognosis and without other possibly confounding co-morbidities. Moreover, our bionformatic analysis revealed that the identified methylation alteration overlaps with enhancers regulating HLA-C expression, suggesting an additional mechanism exploited by SARS-CoV-2 to inhibit this fundamental player in the host’s immune response. HLA-C could therefore represent both a prognostic marker and an excellent therapeutic target, also suggesting a preventive intervention that conjugate a virus-specific antigenic stimulation with an adjuvant increasing the T-cell mediated memory response.

A Review: Highlighting the Links between Epigenetics, COVID-19 Infection, and Vitamin D

The highly transmittable and infectious COVID-19 remains a major threat worldwide, with the elderly and comorbid individuals being the most vulnerable. While vaccines are currently available, therapeutic drugs will help ease the viral outbreak and prevent serious health outcomes. Epigenetic modifications regulate gene expression through changes in chromatin structure and have been linked to viral pathophysiology. Since epigenetic modifications contribute to the life cycle of the virus and host immune responses to infection, epigenetic drugs are promising treatment targets to ameliorate COVID-19. Deficiency of the multifunctional secosteroid hormone vitamin D is a global health threat. Vitamin D and its receptor function to regulate genes involved in immunity, apoptosis, proliferation, differentiation, and inflammation. Amassed evidence also indicates the biological relations of vitamin D with reduced disease risk, while its receptor can be modulated by epigenetic mechanisms. The immunomodulatory effects of vitamin D suggest a role for vitamin D as a COVID-19 therapeutic agent. Therefore, this review highlights the epigenetic effects on COVID-19 and vitamin D while also proposing a role for vitamin D in COVID-19 infections.

EWAS of post-COVID-19 patients shows methylation differences in the immune-response associated gene, IFI44L, three months after COVID-19 infection

Although substantial progress has been made in managing COVID-19, it is still difficult to predict a patient’s prognosis. We explored the epigenetic signatures of COVID-19 in peripheral blood using data from an ongoing prospective observational study of COVID-19 called the Norwegian Corona Cohort Study. A series of EWASs were performed to compare the DNA methylation profiles between COVID-19 cases and controls three months post-infection. We also investigated differences associated with severity and long-COVID. Three CpGs—cg22399236, cg03607951, and cg09829636—were significantly hypomethylated (FDR < 0.05) in COVID-19 positive individuals. cg03607951 is located in IFI44L which is involved in innate response to viral infection and several systemic autoimmune diseases. cg09829636 is located in ANKRD9, a gene implicated in a wide variety of cellular processes, including the degradation of IMPDH2. The link between ANKRD9 and IMPDH2 is striking given that IMPDHs are considered therapeutic targets for COVID-19. Furthermore, gene ontology analyses revealed pathways involved in response to viruses. The lack of significant differences associated with severity and long-COVID may be real or reflect limitations in sample size. Our findings support the involvement of interferon responsive genes in the pathophysiology of COVID-19 and indicate a possible link to systemic autoimmune diseases.

Epigenetic rewiring of pathways related to odour perception in immune cells exposed to SARS-CoV-2 in vivo and in vitro

A majority of SARS-CoV-2 recoverees develop only mild-to-moderate symptoms, while some remain completely asymptomatic. Although viruses, including SARS-CoV-2, may evade host immune responses by epigenetic mechanisms including DNA methylation, little is known about whether these modifications are important in defence against and healthy recovery from COVID-19 in the host. To this end, 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. 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. Altered DNA methylation patterns of COVID-19 convalescents suggest recovery from mild-to-moderate SARS-CoV-2 infection leaves longstanding epigenetic traces. Both in vitro and in vivo exposure caused epigenetic modulation of pathways thataffect odour perception. Future studies should determine whether this reflects host-induced protective antiviral defense or targeted viral hijacking to evade host defence.

Longitudinal Study of DNA Methylation and Epigenetic Clocks Prior to and Following Test-Confirmed COVID-19 and mRNA Vaccination

The host epigenetic landscape rapidly changes during SARS-CoV-2 infection, and evidence suggest that severe COVID-19 is associated with durable scars to the epigenome. Specifically, aberrant DNA methylation changes in immune cells and alterations to epigenetic clocks in blood relate to severe COVID-19. However, a longitudinal assessment of DNA methylation states and epigenetic clocks in blood from healthy individuals prior to and following test-confirmed non-hospitalized COVID-19 has not been performed. Moreover, the impact of mRNA COVID-19 vaccines upon the host epigenome remains understudied. Here, we first examined DNA methylation states in the blood of 21 participants prior to and following test-confirmed COVID-19 diagnosis at a median time frame of 8.35 weeks; 756 CpGs were identified as differentially methylated following COVID-19 diagnosis in blood at an FDR adjusted p-value < 0.05. These CpGs were enriched in the gene body, and the northern and southern shelf regions of genes involved in metabolic pathways. Integrative analysis revealed overlap among genes identified in transcriptional SARS-CoV-2 infection datasets. Principal component-based epigenetic clock estimates of PhenoAge and GrimAge significantly increased in people over 50 following infection by an average of 2.1 and 0.84 years. In contrast, PCPhenoAge significantly decreased in people fewer than 50 following infection by an average of 2.06 years. This observed divergence in epigenetic clocks following COVID-19 was related to age and immune cell-type compositional changes in CD4+ T cells, B cells, granulocytes, plasmablasts, exhausted T cells, and naïve T cells. Complementary longitudinal epigenetic clock analyses of 36 participants prior to and following Pfizer and Moderna mRNA-based COVID-19 vaccination revealed that vaccination significantly reduced principal component-based Horvath epigenetic clock estimates in people over 50 by an average of 3.91 years for those who received Moderna. This reduction in epigenetic clock estimates was significantly related to chronological age and immune cell-type compositional changes in B cells and plasmablasts pre- and post-vaccination. These findings suggest the potential utility of epigenetic clocks as a biomarker of COVID-19 vaccine responses. Future research will need to unravel the significance and durability of short-term changes in epigenetic age related to COVID-19 exposure and mRNA vaccination.