Single-cell transcriptome studies have revealed immune dysfunction in COVID-19 patients, including lymphopenia, T cell exhaustion, and increased levels of pro-inflammatory cytokines, while DNA methylation plays an important role in the regulation of immune response and inflammatory response. The specific cell types of immune responses regulated by DNA methylation in COVID-19 patients will be better understood by exploring the COVID-19 DNA methylation variation at the cell-type level. Here, we developed an analytical pipeline to explore single-cell DNA methylation variations in COVID-19 patients by transferring bulk-tissue-level knowledge to the single-cell level. We discovered that the methylation variations in the whole blood of COVID-19 patients showed significant cell-type specificity with remarkable enrichment in gamma-delta T cells and presented a phenomenon of hypermethylation and low expression. Furthermore, we identified five genes whose methylation variations were associated with several cell types. Among them, S100A9, AHNAK, and CX3CR1 have been reported as potential COVID-19 biomarkers previously, and the others (TRAF3IP3 and LFNG) are closely associated with the immune and virus-related signaling pathways. We propose that they might serve as potential epigenetic biomarkers for COVID-19 and could play roles in important biological processes such as the immune response and antiviral activity.

Postzygotic mutations are acquired in all of the normal tissues throughout an individual’s lifetime and hold clues for identifying mutagenesis causing factors. The process and underlying mechanism of postzygotic mutations in normal tissues is still poorly understood. In this study, we investigated postzygotic mutation spectra in healthy individuals by optimized ultra-deep exome sequencing of time series samples from the same volunteer and samples from different individuals. In cells of blood, sperm, and muscle, we resolved three common types of mutational signature. Two of them are known to represent clock-like mutational processes, and their proportions in mutation profiles associated with polymorphisms of epigenetic regulation genes, suggesting the contribution of personal genetic backgrounds to underlying biological process. Notably, the third signature, characterized by C>T transitions at GpCpN sites, tends to be a feature of diverse normal tissues. Mutations of this type were likely to occur early in embryo development even before the tissue differentiation, as indicated by their relatively high allele frequencies, sharing variants between multiple tissues, and lacking of age-related accumulation. Almost all tumors shown in public datasets did not have this signature detected except for 19.6% of clear cell renal cell carcinoma samples, which featured by activation of the hypoxia-induced signaling pathway. Moreover, in vitro activation of HIF signaling pathway successfully introduced the corresponding mutation profile of this signature in a culture-expanded human embryonic stem cell line. Therefore, embryonic hypoxia may explain this novel signature across multiple normal tissues. Our study suggest hypoxic conditions in the early stage of embryo development may be a crucial factor for the C>T transitions at GpCpN sites and individual genetic background also related to shaping human postzygotic mutation profiles.