key: cord-1019969-x2s22zdb authors: Kaneko, Syuzo; Takasawa, Ken; Asada, Ken; Shinkai, Norio; Bolatkan, Amina; Yamada, Masayoshi; Takahashi, Satoshi; Machino, Hidenori; Kobayashi, Kazuma; Komatsu, Masaaki; Hamamoto, Ryuji title: Epigenetic Mechanisms Underlying COVID-19 Pathogenesis date: 2021-09-02 journal: Biomedicines DOI: 10.3390/biomedicines9091142 sha: 2cd2c2d39bc3149c0cf926418360295e9789e345 doc_id: 1019969 cord_uid: x2s22zdb In 2019, a novel severe acute respiratory syndrome called coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was reported and was declared a pandemic by the World Health Organization (WHO) in March 2020. With the advancing development of COVID-19 vaccines and their administration globally, it is expected that COVID-19 will converge in the future; however, the situation remains unpredictable because of a series of reports regarding SARS-CoV-2 variants. Currently, there are still few specific effective treatments for COVID-19, as many unanswered questions remain regarding the pathogenic mechanism of COVID-19. Continued elucidation of COVID-19 pathogenic mechanisms is a matter of global importance. In this regard, recent reports have suggested that epigenetics plays an important role; for instance, the expression of angiotensin I converting enzyme 2 (ACE2) receptor, an important factor in human infection with SARS-CoV-2, is epigenetically regulated; further, DNA methylation status is reported to be unique to patients with COVID-19. In this review, we focus on epigenetic mechanisms to provide a new molecular framework for elucidating the pathogenesis of SARS-CoV-2 infection in humans and of COVID-19, along with the possibility of new diagnostic and therapeutic strategies. Coronavirus disease 2019 (COVID- 19) is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Its symptoms are varied and range from mild to severe; common symptoms include headache, loss of sense of smell and taste, nasal congestion and leakage, cough, myalgia, sore throat, fever, diarrhea, and difficulty breathing [1] [2] [3] [4] [5] . Of those with symptoms significant enough to be classified as patients in a study in China, most (81%) had mild to moderate symptoms (up to mild pneumonia), 14% had severe symptoms (dyspnea, hypoxia, or more than lung involvement on imaging), and 5% had severe symptoms (respiratory failure, shock, or multiple organ failure) [6] . COVID-19 spread around the world within a short time and was declared a pandemic by the World Health Organization (WHO) on 11 March 2020. Since the end of 2020, the emergence of SARS-CoV-2 variants with genetic mutations that may affect transmission, severity, and antigenicity has become a major problem. In particular, outbreaks of B.1.1.7 (Alpha variant), first detected in the UK, B.1.351 (Beta variant), first detected in South South Africa, P.1 (Gamma variant), first detected in Japan among returnees from Brazil, and B.1.617.2 (Delta variant), first detected in India, are of global concern [7, 8] . Under these circumstances, the development of vaccines aiming to provide acquired immunity to humans against SARS-CoV-2 has been progressing, and as of July 2021, multiple vaccines with varying methods and manufacturers have been developed, ranging from those that have already been inoculated to those under development (mRNA vaccines, DNA vaccines, virus vector vaccines, inactivated vaccines, recombinant protein vaccines, peptide vaccines, etc.) [9] . In particular, the results of a large placebo-controlled trial of vaccines from Pfizer and Moderna have shown that two repeated doses of the vaccine are highly effective in preventing COVID-19 (more than 90%) [10, 11] . Based on the results of these clinical trials, vaccination with multiple types of COVID-19 vaccines developed by several pharmaceutical companies began worldwide at the end of 2020 or in early 2021 [12] . While hopes for pandemic containment through vaccination are increasing, the emergence of SARS-CoV-2 variants with amino acid mutations in the antigenic determinants of spike proteins, as described above, has been reported worldwide. Further, since there is a possibility that unknown SARS-CoV-2 variants may emerge in the future, it is necessary to continue with basic infection control measures without overconfidence in the effectiveness of vaccines. Currently, drugs such as Remdesivir, Dexamethasone, Baricitinib and Heparin are used as treatments for COVID-19; however, there are still few specific and effective treatments for COVID-19 [13, 14] . This is because there are still many unanswered questions regarding the pathogenic mechanism of COVID-19. Elucidating the pathogenic mechanism of COVID-19 is thus essential to establish effective treatment for patients with severe disease and to develop novel therapeutic agents. In this regard, epigenetics status such as histone modification of host cells is known to be altered upon infection with RNA viruses including coronaviruses, and the importance of epigenetics in the pathogenic mechanism of SARS-CoV-2 infection in humans and COVID-19 was recently pointed out. Epigenetics is commonly defined as the study of heritable phenotypic changes without altering the DNA sequence. The Greek prefix epi-(ἐπι) in epigenetics implies a function "above" or "in addition to" the traditional genetic base [15] . Over the past two decades, epigenetic regulators have been implicated as critical factors in many pathways related to the development and progression of cancer and other diseases, including cell cycle regulation, invasiveness, signaling pathways, chemotherapy resistance, and immune evasion [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] . The three basic systems of epigenetic regulation are DNA methylation of gene regulatory regions; histone protein modifications such as methylation, acetylation, phosphorylation, and sumoylation; and non-coding RNAs [15] . Many techniques for epigenetics analysis have already been developed, and this field is steadily undergoing technological innovation [15, 28, 29] . In this review, we present the latest findings on the importance of epigenetics in the mechanisms of human infection with SARS-CoV-2 and pathogenesis of COVID-19, and discuss future diagnostic and therapeutic strategies for COVID-19 targeting epigenetics. The SARS-CoV-2 genome is composed of single-stranded RNA [30] . SARS-CoV-2 shares 96% genome sequence identity with BatCoV RaTG13 [31] ; 90% identity with Pangolin-CoVs [32] ; 88% genome sequence identity with two bat-derived coronaviruses, bat-SL-CoVZC45 and bat-SL-CoVZXC21; 79% identity with SARS-CoV; and 50% with Middle East respiratory syndrome coronavirus [31, 32] . Although SARS-CoV-2 has highly identical sequences with the above viruses, SARS-CoV-2 only has a functional furin cleavage site at the spike (S) protein [33, 34] . The S protein expressed on the surface of the viral particles is essential to the initial steps of coronavirus infection. The S protein comprises the receptor-binding subunit S1 and the membrane-fusion subunit S2 [35] [36] [37] . The receptor binding S1 consists of two subdomains, an N-terminal domain and a C-terminal domain [38, 39] . πι) in epigenetics implies a function "above" or "in addition to" the traditional genetic base [15] . Over the past two decades, epigenetic regulators have been implicated as critical factors in many pathways related to the development and progression of cancer and other diseases, including cell cycle regulation, invasiveness, signaling pathways, chemotherapy resistance, and immune evasion [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] . The three basic systems of epigenetic regulation are DNA methylation of gene regulatory regions; histone protein modifications such as methylation, acetylation, phosphorylation, and sumoylation; and non-coding RNAs [15] . Many techniques for epigenetics analysis have already been developed, and this field is steadily undergoing technological innovation [15, 28, 29] . In this review, we present the latest findings on the importance of epigenetics in the mechanisms of human infection with SARS-CoV-2 and pathogenesis of COVID-19, and discuss future diagnostic and therapeutic strategies for COVID-19 targeting epigenetics. The SARS-CoV-2 genome is composed of single-stranded RNA [30] . SARS-CoV-2 shares 96% genome sequence identity with BatCoV RaTG13 [31] ; 90% identity with Pangolin-CoVs [32] ; 88% genome sequence identity with two bat-derived coronaviruses, bat-SL-CoVZC45 and bat-SL-CoVZXC21; 79% identity with SARS-CoV; and 50% with Middle East respiratory syndrome coronavirus [31, 32] . Although SARS-CoV-2 has highly identical sequences with the above viruses, SARS-CoV-2 only has a functional furin cleavage site at the spike (S) protein [33, 34] . The S protein expressed on the surface of the viral particles is essential to the initial steps of coronavirus infection. The S protein comprises the receptor-binding subunit S1 and the membrane-fusion subunit S2 [35] [36] [37] . The receptor binding S1 consists of two subdomains, an N-terminal domain and a C-terminal domain [38, 39] . The viral S protein binds to angiotensin I converting enzyme 2 (ACE2) as entry receptor. In addition, the S protein of SARS-CoV-2 cleaves by transmembrane serine protease-2 (TMPRSS2) and furin protease, key factors in its host cell entry [37, 40, 41] . Interestingly, studies have attempted to explain the global pandemic of COVID-19 based on hardness and phylogenetic analysis of the outer shell (M protein) of SARS-CoV-2, and there may be other factors besides the S protein of the virus that characterize SARS-CoV-2 [42] . SARS-CoV-2 is hypothesized to use clathrin-dependent endocytosis [43] [44] [45] , the most common endosomal pathway, to facilitate viral entry. Caveolae-dependent uptake, which is the other pathway, is considered controversial regarding coronavirus entry into host cells and is potentially dependent on cell type [46, 47] . Following genome release into the cytosol, viral genomic RNA is uncoated and translated for the synthesis of nonstructural viral proteins. Importantly, this intrinsic disorganization of viral proteins is an inherent feature and a strategy of viruses to disrupt host nucleocytoplasmic transport to benefit their own replication [48] . Viral RNA is then replicated and subgenomic RNAs are translated in double-membrane vesicles, as previously reported in several viruses, including coronavirus [49] [50] [51] [52] . Multiple components, such as host membrane-derived double-membrane spherules, convoluted membranes, and the endoplasmic reticulum (ER), are used to protect and support the transcribed genomic RNA and are hallmarks of viral infections. Translated structural proteins then translocate into the ER and transit through the ER-to-Golgi intermediate compartment, wherein N-encapsidated, newly produced genomic RNA is packaged in the lumen of secretory vesicles [53, 54] . Viruses are then secreted from the infected cells by exocytosis, and the aforementioned process is repeated in new cells ( Figure 1) . Notably, D614G substitution in the receptor-binding domain of the S protein of SARS-CoV-2 results in the change of a single amino acid, which enhances binding affinity to the ACE2 receptor and increases viral entry into host cells [55] . Thus, to enhance host cell defenses, it is important to understand how viral RNA infects host cells. The viral S protein binds to angiotensin I converting enzyme 2 (ACE2) as entry receptor. In addition, the S protein of SARS-CoV-2 cleaves by transmembrane serine protease-2 (TMPRSS2) and furin protease, key factors in its host cell entry [37, 40, 41] . Interestingly, studies have attempted to explain the global pandemic of COVID-19 based on hardness and phylogenetic analysis of the outer shell (M protein) of SARS-CoV-2, and there may be other factors besides the S protein of the virus that characterize SARS-CoV-2 [42] . SARS-CoV-2 is hypothesized to use clathrin-dependent endocytosis [43] [44] [45] , the most common endosomal pathway, to facilitate viral entry. Caveolae-dependent uptake, which is the other pathway, is considered controversial regarding coronavirus entry into host cells and is potentially dependent on cell type [46, 47] . Following genome release into the cytosol, viral genomic RNA is uncoated and translated for the synthesis of non-structural viral proteins. Importantly, this intrinsic disorganization of viral proteins is an inherent feature and a strategy of viruses to disrupt host nucleocytoplasmic transport to benefit their own replication [48] . Viral RNA is then replicated and subgenomic RNAs are translated in double-membrane vesicles, as previously reported in several viruses, including coronavirus [49] [50] [51] [52] . Multiple components, such as host membrane-derived double-membrane spherules, convoluted membranes, and the endoplasmic reticulum (ER), are used to protect and support the transcribed genomic RNA and are hallmarks of viral infections. Translated structural proteins then translocate into the ER and transit through the ER-to-Golgi intermediate compartment, wherein N-encapsidated, newly produced genomic RNA is packaged in the lumen of secretory vesicles [53, 54] . Viruses are then secreted from the infected cells by exocytosis, and the aforementioned process is repeated in new cells ( Figure 1) . Notably, D614G substitution in the receptor-binding domain of the S protein of SARS-CoV-2 results in the change of a single amino acid, which enhances binding affinity to the ACE2 receptor and increases viral entry into host cells [55] . Thus, to enhance host cell defenses, it is important to understand how viral RNA infects host cells. RNA editing induces substitution in the viral genome and can potentially affect viral infection. RNA epigenetics can be divided into two main categories according to the mechanisms involved. One occurs when RNA modifications are catalyzed by methyltransferases (known as "writer proteins") such as METTL3, and the other occurs when RNA editing is catalyzed by editing enzymes. RNA editing by host deaminases is an innate restriction process used as a defense mechanism against viral infection [56] . RNA sequencing of SARS-CoV-2 isolated from the bronchoalveolar lavage fluid of patients with COVID-19 revealed instances of RNA editing, such as A to G and U to C substitutions, and was accompanied by evidence of adenine to inosine (A-to-I) reactions catalyzed by RNA-specific adenosine deaminase (ADAR). In addition, restriction of viral replication was observed with C to U and G to A substitutions by apolipoprotein B mRNA editing enzyme catalytic polypeptide [57] . Furthermore, using RNA antisense purification and quantitative mass spectrometry in SARS-CoV-2-infected human cells, ADAR was found to interact with SARS-CoV-2 RNA [58] . The authors identified 699 proteins, of which 583 were detected with multiple peptides using liquid chromatography with tandem mass spectrometry. Viruses are known to utilize the host RNA editing machinery, and the resulting substitutions can lead to both proviral and antiviral effects. Recent studies have described the distribution of RNA substitution in the SARS-CoV-2 genome, and have provided insights into the random mutation and the specific substitutions induced by editing enzymes [59, 60] . As ACE2 is an essential intracellular receptor that binds to the S protein encoded by SARS-CoV-2 to cause an infection, understanding the status and mechanism of ACE2 gene expression is essential. Generally, ACE2 is expressed in various human organs, and its organ-and cell-specific expression suggests that it is involved in regulating cardiovascular and renal function and fertility [61] . Interestingly, single-cell RNA sequence datasets revealed that ACE2 was coexpressed with TMPRSS2 within lung type II pneumocytes, ileal absorptive enterocytes, and nasal goblet secretory cells [62] . It is also known that young children exhibit lower ACE2 expression compared to the adults [63] ; the low risk of infection in children may thus be attributed to this age-dependent ACE2 expression [64] . Studies using biochemical approaches have examined the binding of transcription factors to DNA sequences upstream of the transcription start site (TSS) of ACE2. Using the human embryonic kidney cell line HEK293, it has been suggested that there is a functional hepatocyte nuclear factor 1β (HNF1β) binding site in the promoter region and HNF1β in turn promotes ACE2 transcription [65] . Furthermore, HNF1α and HNF1β, which bind to three HNF1-binding motifs that are highly conserved among mammalian species, cooperatively regulate ACE2 activity in insulinoma cells [66] . The ACE2 gene is known to produce various transcripts that are stimulated by interferons and has several different TSSs, wherein the selection of TSSs varies from organ to organ [67] [68] [69] . As the ACE2 gene is located on the X chromosome, a gene dosage effect potentially regulates gene expression. However, new evidence indicates that COVID-19 is a genderbiased disease influenced by myriad variables ranging from biological to social factors, and it is thus difficult to determine the relationship between SARS-CoV-2 infection and sex [70] . DNA methylation is an epigenetic modification reported to be associated with various clinical conditions, such as cancer and asthma [71] [72] [73] . The level of DNA methylation around promoters is thought to regulate the expression of proximal genes by altering the affinity of transcription factor binding. As mentioned above, Pedersen et al. performed promoter deletion analysis for ACE2 and identified two regulatory promoter regions, the distal promoter region and proximal promoter region at −1509 bp to −928 bp and −454 bp to −1 bp from the TSS, respectively [66] . These promoter regions contain several transcription factor binding motifs such as HNF1α, HNF1β, GATA, FOXA, YY1, and C/EBPβ [66, [74] [75] [76] (Figure 2a) . The binding motif of C/EBPβ contains a CpG site, and the hypomethylation of the C/EBPβ binding motif is reported to induce angiotensinogen expression [77] . Therefore, DNA methylation status around the ACE2 promoter may be involved in regulating ACE2 gene expression. Comparative analysis of DNA methylation in various tissues has revealed that several CpG sites are hypomethylated in lung epithelial cells, which exhibit high ACE2 expression [78] (Figure 2b) . Further, Cardenas et al. reported multiple methylation sites in the vicinity of the TSS and ACE2 gene body with variable degrees of methylation dependent on sex and race and a higher level of hypomethylation in females compared to that in males [79] . Although the platform used for the analysis was from a previous generation, Yang et al. reported that patients with asthma showed differentially methylated CpG sites in similar locations [80] . Wang et al. systematically analyzed the aberrant expression of ACE2 and TMPRSS2 in multiple human cancers and found that colorectal cancers with elevated gene expression had reduced DNA methylation levels. Since cancer is considered a risk factor for COVID-19, the outbreak of COVID-19 may require additional care, especially for colorectal cancer patients [81] . to −1 bp from the TSS, respectively [66] . These promoter regions contain several transcription factor binding motifs such as HNF1α, HNF1β, GATA, FOXA, YY1, and C/EBPβ [66, [74] [75] [76] (Figure 2a) . The binding motif of C/EBPβ contains a CpG site, and the hypomethylation of the C/EBPβ binding motif is reported to induce angiotensinogen expression [77] . Therefore, DNA methylation status around the ACE2 promoter may be involved in regulating ACE2 gene expression. Comparative analysis of DNA methylation in various tissues has revealed that several CpG sites are hypomethylated in lung epithelial cells, which exhibit high ACE2 expression [78] (Figure 2b) . Further, Cardenas et al. reported multiple methylation sites in the vicinity of the TSS and ACE2 gene body with variable degrees of methylation dependent on sex and race and a higher level of hypomethylation in females compared to that in males [79] . Although the platform used for the analysis was from a previous generation, Yang et al. reported that patients with asthma showed differentially methylated CpG sites in similar locations [80] . Wang et al. systematically analyzed the aberrant expression of ACE2 and TMPRSS2 in multiple human cancers and found that colorectal cancers with elevated gene expression had reduced DNA methylation levels. Since cancer is considered a risk factor for COVID-19, the outbreak of COVID-19 may require additional care, especially for colorectal cancer patients [81] . The DNA methylation status in the host is reported to change with bacterial and viral infections [82] , and analysis of the DNA methylation state is expected to aid in estimating the infection history and future disease severity [83] . Balnis et al. performed a comprehensive DNA methylation analysis using blood from hospitalized patients (COVID-19 and non-COVID-19) and healthy individuals [84] (Figure 2c ) and found that although global DNA methylation levels did not differ between healthy and hospitalized individuals, 1089 hypo-differentially methylated regions (DMRs) and 416 hyper DMRs were found [84] . In addition, a comparison between patients with and without COVID-19 revealed 47 COVID-19 patient-specific DMRs, 36 of which were inversely correlated with the expression of neighboring genes. Moreover, gene ontology analysis of these 36 DMRs showed that ontologies involving defense response to the virus (27/36 DMRs) and interferon signaling (19/36 DMRs) were enriched [84] . Furthermore, 77 DMRs were identified by comparing patients with mild and severe COVID-19. The results of hierarchical cluster analysis of these 77 DMRs were consistent with the severity assessment indicated by the COVID-GRAM score [85] , indicating that these DMRs can potentially act as biomarkers [84] . Histone modifications are involved in human health conditions, aging, neurological diseases, and cancer development, which have been described in the great reviews [86] [87] [88] . An association between ACE2 gene expression and histone modifications was reported previously [89] . Although core histones in nucleosomes are tightly arranged, histonemodifying enzymes can modify their tails. The counteracting action of various enzymes results in reversible epigenetic modifications. As described in other reviews, the action of enzymes, such as histone acetyltransferases, are balanced by those of histone deacetylases and the action of histone methyltransferases are countered by the action of histone demethylase, to maintain the epigenome [90, 91] . Silent information regulator T1 (SIRT1), a histone deacetylase, is reported to be involved in the transcriptional regulation of ACE2 expression [92] , and this role is thought to be related to the protective role of SIRT1 against cellular stress. Therefore, nonsteroidal anti-inflammatory drugs that inhibit SIRT1 deacetylase activity are expected to exert unexpected anti-infective effects [93] . Notably, administration of atorvastatin to rabbits on a high cholesterol diet has been reported to upregulate ACE2 expression via tissue-specific and promoter-specific histone modifications compared to that in their corresponding controls [94] . Importantly, the elevated ACE2 expression in the lungs of patients with severe COVID-19-related complications may be due to histone modifications of several genes such as HAT1, HDAC2, and lysine demethylase 5 B (KDM5B) [89] . In particular, KDM5B affects chromatin accessibility by removing activated chromatin marks such as the dimethylation and trimethylation of histone H3 (H3K4) from lysine 4, thereby contributing to transcriptional regulation and DNA repair [95] . Inhibition of KDM5B in breast cancer cells has been shown to induce an interferon response, making the cells less susceptible to DNA and RNA viral infection [96] . Therefore, KDM5 demethylase is expected to be a potential target for COVID-19 prevention. Furthermore, diarylheptanoids, also known as diphenylheptanoids have been reported to induce epigenetic silencing of the ACE2 gene mediated by HMGB1 [97] [98] [99] [100] [101] , and this has attracted attention as a possible way to prevent COVID-19 infection. Intriguingly, it has been reported that the expression of TMPRSS2 and ACE2 was decreased by therapies directly targeting androgen receptor (AR) and inhibitors of bromodomain and extra terminal domain (BET) proteins, which are known epigenetic regulators of AR transcriptional activity. Furthermore, these treatments reduced SARS-CoV-2 infection in a cellular model. Therefore, these findings support further research on AR and BET inhibitors as potential treatments for COVID-19 [102, 103] . Histone modifications are not only involved in the transcriptional regulation of ACE2, but also have important pathophysiological functions in COVID-19. Elevated citrullination of histone H3 (Cit-H3) in the serum of COVID-19 patients was recently reported [104] . Citrullination (deamination) of arginine residues in histones leads to generation of un-charged citrullines, facilitates chromatin decondensation, and enhances the accessibility of chromatin [105] [106] [107] (Figure 3 ). Histone citrullination is catalyzed by the enzyme peptidylarginine deiminase, a family of calcium-dependent enzymes that regulate immune activity and are involved in neutrophil extracellular traps (NETs) [106, 107] . Pathologically, NETs are a biological defense response of neutrophils, which are triggered by infection; cell death via NETosis, which is different from cell necrosis and apoptosis, results in the release of neutrophil DNA into the extracellular space to form net-like structures [108] . As there is a positive correlation between Cit-H3 levels and platelet counts, Cit-H3 may contribute to abnormal platelet counts. Furthermore, a recent report suggested that NETs could be markers of COVID-19 severity [104] . Therefore, in the future it will be imperative to investigate the prognostic impact of histone Cit-H3 conversion on SARS-CoV-2 infection and venous thrombosis from both molecular biological and clinical perspectives. model. Therefore, these findings support further research on AR and BET inhibitors as potential treatments for COVID-19 [102, 103] . Histone modifications are not only involved in the transcriptional regulation of ACE2, but also have important pathophysiological functions in COVID-19. Elevated citrullination of histone H3 (Cit-H3) in the serum of COVID-19 patients was recently reported [104] . Citrullination (deamination) of arginine residues in histones leads to generation of uncharged citrullines, facilitates chromatin decondensation, and enhances the accessibility of chromatin [105] [106] [107] (Figure 3 ). Histone citrullination is catalyzed by the enzyme peptidylarginine deiminase, a family of calcium-dependent enzymes that regulate immune activity and are involved in neutrophil extracellular traps (NETs) [106, 107] . Pathologically, NETs are a biological defense response of neutrophils, which are triggered by infection; cell death via NETosis, which is different from cell necrosis and apoptosis, results in the release of neutrophil DNA into the extracellular space to form net-like structures [108] . As there is a positive correlation between Cit-H3 levels and platelet counts, Cit-H3 may contribute to abnormal platelet counts. Furthermore, a recent report suggested that NETs could be markers of COVID-19 severity [104] . Therefore, in the future it will be imperative to investigate the prognostic impact of histone Cit-H3 conversion on SARS-CoV-2 infection and venous thrombosis from both molecular biological and clinical perspectives. Recently, the number of reports regarding the role of microRNAs (miRNAs) in regulating the ACE2 gene has been increasing. miRNAs are non-coding, single-stranded, small RNAs that function by post-transcriptionally inhibiting the expression of target genes in various organisms, from viruses to higher eukaryotes [109] [110] [111] [112] [113] [114] . These small RNAs function by complementarily binding to the 3′-untranslated region (UTR) and sometimes the 5′-UTR or coding regions of their target mRNAs. miRNAs are involved in numerous disease and physiological functions, such as cancer [115] , degenerative neuro disorders [116] , cardiovascular disease [117] , and immunity [118] . miRNAs are also often identified as biomarkers of viral infections, and their contribution to host-pathogen interactions during viral infections is being investigated [119] [120] [121] . Recently, ACE2 mRNA and ACE2 protein levels were reported to be repressed by miR-200c in primary rat cardiomyocytes and human-induced pluripotent stem cell-derived cardiomyocytes, suggesting a direct link between miR-200c and the regulation of ACE2 expression [122] . Many other miRNAs are thought to have the potential to regulate ACE2, but most of them have not been carefully investigated [123, 124] . Furthermore, it has been suggested that TMPRSS2 transcription is regulated by miR-98-5p in two human endothelial cell types, derived from the lung and from the umbilical vein, suggesting strict control by the miRNA-target network [125] . Recently, the number of reports regarding the role of microRNAs (miRNAs) in regulating the ACE2 gene has been increasing. miRNAs are non-coding, single-stranded, small RNAs that function by post-transcriptionally inhibiting the expression of target genes in various organisms, from viruses to higher eukaryotes [109] [110] [111] [112] [113] [114] . These small RNAs function by complementarily binding to the 3 -untranslated region (UTR) and sometimes the 5 -UTR or coding regions of their target mRNAs. miRNAs are involved in numerous disease and physiological functions, such as cancer [115] , degenerative neuro disorders [116] , cardiovascular disease [117] , and immunity [118] . miRNAs are also often identified as biomarkers of viral infections, and their contribution to host-pathogen interactions during viral infections is being investigated [119] [120] [121] . Recently, ACE2 mRNA and ACE2 protein levels were reported to be repressed by miR-200c in primary rat cardiomyocytes and human-induced pluripotent stem cell-derived cardiomyocytes, suggesting a direct link between miR-200c and the regulation of ACE2 expression [122] . Many other miRNAs are thought to have the potential to regulate ACE2, but most of them have not been carefully investigated [123, 124] . Furthermore, it has been suggested that TMPRSS2 transcription is regulated by miR-98-5p in two human endothelial cell types, derived from the lung and from the umbilical vein, suggesting strict control by the miRNA-target network [125] . Long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) have also been implicated in SARS-CoV-2 infection [126] [127] [128] . Like miRNAs, lncRNAs and circRNAs are a class of non-coding RNAs, but lncRNAs are typically 1000-10,000 residues in length, and circRNAs are circular RNAs processed by back splicing [129, 130] . Comparison analysis between normal human bronchial epithelial cells and SARS-CoV-2 infected cells identified several differentially expressed lncRNAs, and the differences of the expression were also confirmed among COVID-19 patients [126] . Involvement of circRNAs in several viral infections such as hepatitis B virus and human papillomavirus were recently reported [128] . For infection of SARS-CoV-2, the two circRNAs, Ppp1r10 and C330019G07RiK, were reported as a part of quintuple competing endogenous RNA networks consisting of miRNA, lncRNA, circRNAs, mRNA, and transcription factor by in silico study [127] . The cutting-edge technologies such as next-generation sequencing (NGS) and CRISPR systems have been contributing to the finding of new insights and development of new diagnostic methods for COVID-19 [131] [132] [133] . The studies using NGS revealed the whole sequences of SARS-CoV-2 at just 4 months from the COVID-19 pandemic [134] , and the sequences were beneficial for understanding the virus genomic variants and pathogenetic mechanisms [131] [132] [133] . These innovative technologies are useful for understanding not only the virological features of SARS-CoV-2 but also the response of patients in a molecular biological aspect. A series of reports have identified factors associated with SARS-CoV-2 infection using genome-wide screening analysis [97, 135, 136] . Wei et al. identified epigenetic regulatory genes involved in diverse biological processes, such as chromatin remodeling, histone modification, intracellular signaling, and RNA regulation, as candidate host genes that affect SARS-CoV-2 infection. For example, HMGB1 appears to play a novel role in epigenetically regulating ACE2 expression, thus increasing susceptibility to SARS-CoV-2 infection (Figure 2a ). HMGB1 encodes a non-histone nuclear DNA-binding protein, belonging to the high mobility group-box superfamily, that regulates transcription and is involved in DNA organization [137, 138] . This protein is also known to act as a danger-associated molecular pattern molecule that amplifies immune responses during tissue injury [139] . Notably, HMGB1 regulates ACE2 expression in a cell-intrinsic manner rather than functioning as a cytokine or alarmin, suggesting a precise mechanism for HMGB1 involvement in SARS-CoV-2 infection [97] . Furthermore, the gene encoding the SWI/SNF chromatin remodeling complex was identified as a SARS-CoV-2 provirus, demonstrating the importance of this complex as a pathogen [97] . Previously, the SWI/SNF complex has been reported to be composed of the ATPase subunits SMARCA2 or SMARCA4, and it catalyzes non-catalytic scaffold core expression via ARID1A, which has no inherent DNA sequence specificity [140, 141] . Its target specificity is thought to be conferred by the recruitment of DNA-binding proteins to target sites on the genome and by sliding nucleosomes to regulate chromatin accessibility and gene expression. However, TMPRSS2, TMPRSS4, and FURIN [142, 143] , which are considered proviral genes for SARS-CoV-2 infection, were not identified in this screening. This discrepancy may be due to technical difficulties, such as the use of CRISPR libraries and variation in gRNA expression levels in the cell type used. In this review, we presented the latest findings on the mechanisms of SARS-CoV-2 infection and COVID-19 pathogenesis with a focus on epigenetics. Since the discovery of the double helix structure of DNA by J.D. Watson and F.H. Crick in 1953 [144] , molecular biology has emerged and various life mechanisms have been elucidated at the genetic level. In addition, with the completion of the Human Genome Project in 2003 and the elucidation of the human whole genome, there has been a growing momentum for the application of genomic information to medicine, known as genomic medicine [15, 27, 28] . U.S. President Barack Obama announced the Precision Medicine Initiative in 2015, leading to a worldwide push for precision medicine, wherein patients are selected for optimal treatment based on genomic mutations [145] . On the contrary, since the latter half of the 20th century, the importance of changes in gene expression or cellular phenotype that are inherited after cell division without changes in DNA sequence has been indicated, and a new research field called epigenetics has developed [15] . In fact, when elucidating the mechanisms of diseases such as cancer, it is difficult to determine the entire pathogenesis based on genetic mutations alone, and the importance of conducting epigenetics research at the same time has been recognized worldwide. Further, considering the ongoing promotion of precision medicine, it is difficult to adequately screen patients based on genetic mutations alone, and epigenetic information must also be used appropriately. In clinical practice, DNA methylation is used for diagnosis [146] , and HDAC inhibitors and DNA methylation inhibitors are used as therapeutic agents [147, 148] . Furthermore, clinical research and clinical trials of new drug types, such as histone methyltransferase inhibitors and histone demethylase inhibitors, are currently underway worldwide for clinical application [149] [150] [151] [152] . Considering these developments, it is important to conduct research focusing on epigenetics to elucidate pathological mechanisms and to develop new diagnostic and therapeutic methods for COVID-19. In this article, we also introduced the observation of DMRs unique to patients with COVID-19, which may serve as biomarkers for predicting the disease severity and the possibility of functioning as a novel therapeutic agent by epigenetically regulating ACE2 expression. Since the involvement of epigenetics factors such as histone deacetylase in ACE2 expression has been indicated, there is a possibility of more effective treatment for COVID-19, for example, by combining HDAC inhibitors, which are already approved by the FDA, with Remdesivir or Dexamethasone. Since COVID-19 was first reported in December 2019, it has only been one year and eight months, and many unanswered questions remain regarding its pathological mechanisms. It is hoped that COVID-19 will be controlled as vaccination progresses, but with SARS-CoV-2 variants being reported consecutively, the situation remains unpredictable. We thus believe that it is of global importance to continue elucidating the pathological mechanisms of COVID-19 and to develop new diagnostic methods and therapeutic strategies to free people worldwide from the threat of COVID-19 and to enable them to lead safe and secure daily lives. Epigenetics is also an important subject of research, and we hope that the clinical importance of epigenetics in the pathogenesis of COVID-19 will be further clarified in the future by conducting analyses using larger clinical samples. Prevalence and characteristics of fever in adult and paediatric patients with coronavirus disease 2019 (COVID-19): A systematic review and meta-analysis of 17,515 patients Prevalence of Headache in Patients with Coronavirus Disease 2019 (COVID-19): A Systematic Review and Meta-Analysis of 14 Prevalence of Olfactory Dysfunction in Coronavirus Disease 2019 (COVID-19): A Meta-analysis of 27,492 Patients Prevalence and Characteristics of Taste Disorders in Cases of COVID-19: A Meta-analysis of 29 Smell and Taste Dysfunction in Patients With COVID-19: A Systematic Review and Meta-analysis Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72314 Cases From the Chinese Center for Disease Control and Prevention SARS-CoV-2 testing and sequencing for international arrivals reveals significant cross border transmission of high risk variants into the United Kingdom Features, Evaluation, and Treatment of Coronavirus (COVID-19) Evolution of the COVID-19 vaccine development landscape Safety and Efficacy of the BNT162b2 mRNA COVID-19 Vaccine Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine The COVID-19 pandemic: A threat to antimicrobial resistance containment Drug treatments for COVID-19: Living systematic review and network meta-analysis Various theranostics and immunization strategies based on nanotechnology against COVID-19 pandemic: An interdisciplinary view Epigenetics Analysis and Integrated Analysis of Multiomics Data, Including Epigenetic Data, Using Artificial Intelligence in the Era of Precision Medicine The oncogenic polycomb histone methyltransferase EZH2 methylates lysine 120 on histone H2B and competes ubiquitination PRMT6 increases cytoplasmic localization of p21CDKN1A in cancer cells through arginine methylation and makes more resistant to cytotoxic agents Oncogenic histone methyltransferase EZH2: A novel prognostic marker with therapeutic potential in endometrial cancer The novel prognostic marker, EHMT2, is involved in cell proliferation via HSPD1 regulation in breast cancer Deregulation of the Histone Lysine-Specific Demethylase 1 Is Involved in Human Hepatocellular Carcinoma Novel prognostic marker PRMT1 regulates cell growth via downregulation of CDKN1A in HCC WHSC1 Promotes Oncogenesis through Regulation of NIMA-related-kinase-7 in Squamous Cell Carcinoma of the Head and Neck Overexpression of LSD1 contributes to human carcinogenesis through chromatin regulation in various cancers Dysregulation of PRMT1 and PRMT6, Type I arginine methyltransferases, is involved in various types of human cancers Critical role of lysine 134 methylation on histone H2AX for gamma-H2AX production and DNA repair Leveraging epigenetics to enhance the efficacy of immunotherapy Analysis of Whole Genome and Epigenome Data Using Machine Learning Technology: Toward the Establishment of Precision Oncology Application of Artificial Intelligence Technology in Oncology: Towards the Establishment of Precision Medicine Chromatin Analysis of FFPE Tissues Using a Dual-Arm Robot with Clinical Potential The Architecture of SARS-CoV-2 Transcriptome A pneumonia outbreak associated with a new coronavirus of probable bat origin Disorder Analysis Suggests That Pangolins Offered a Window for a Silent Spread of an Attenuated SARS-CoV-2 Precursor among Humans Cleavage Site in the Spike Protein of SARS-CoV-2 Is Essential for Infection of Human Lung Cells The proximal origin of SARS-CoV-2 Pre-fusion structure of a human coronavirus spike protein Structure of SARS coronavirus spike receptor-binding domain complexed with receptor Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus Receptor usage and cell entry of porcine epidemic diarrhea coronavirus Possible Transmission Flow of SARS-CoV-2 Based on ACE2 Features Efficient activation of the severe acute respiratory syndrome coronavirus spike protein by the transmembrane protease TMPRSS2 The structural basis of accelerated host cell entry by SARS-CoV-2dagger Computational, Experimental, and Clinical Evidence for a Specific but Peculiar Evolutionary Nature of (COVID-19) SARS-CoV-2 SARS-CoV-2 infects cells following viral entry via clathrin-mediated endocytosis Clathrin-dependent entry of severe acute respiratory syndrome coronavirus into target cells expressing ACE2 with the cytoplasmic tail deleted The clathrin endocytic pathway in viral infection Lipid rafts are involved in SARS-CoV entry into Vero E6 cells SARS coronavirus entry into host cells through a novel clathrinand caveolae-independent endocytic pathway The Role of Protein Disorder in Nuclear Transport and in Its Subversion by Viruses A molecular pore spans the double membrane of the coronavirus replication organelle SARS-coronavirus replication is supported by a reticulovesicular network of modified endoplasmic reticulum Formation of the poliovirus replication complex requires coupled viral translation, vesicle production, and viral RNA synthesis Cellular origin and ultrastructure of membranes induced during poliovirus infection The intracellular cargo receptor ERGIC-53 is required for the production of infectious arenavirus, coronavirus, and filovirus particles Poliovirus infection blocks ERGIC-to-Golgi trafficking and induces microtubule-dependent disruption of the Golgi complex Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus Critical Roles of N(6)-Methyladenosine (m(6)A) in Cancer and Virus Infection Evidence for host-dependent RNA editing in the transcriptome of SARS-CoV-2 The SARS-CoV-2 RNA-protein interactome in infected human cells Point mutation bias in SARS-CoV-2 variants results in increased ability to stimulate inflammatory responses Host Immune Response Driving SARS-CoV-2 Evolution A novel angiotensin-converting enzyme-related carboxypeptidase (ACE2) converts angiotensin I to angiotensin 1-9 SARS-CoV-2 Receptor ACE2 Is an Interferon-Stimulated Gene in Human Airway Epithelial Cells and Is Detected in Specific Cell Subsets across Tissues Nasal Gene Expression of Angiotensin-Converting Enzyme 2 in Children and Adults Epidemiology of COVID-19 Among Children in China Identification of target genes of the transcription factor HNF1beta and HNF1alpha in a human embryonic kidney cell line The transcription factor HNF1alpha induces expression of angiotensin-converting enzyme 2 (ACE2) in pancreatic islets from evolutionarily conserved promoter motifs Interferons and viruses induce a novel truncated ACE2 isoform and not the full-length SARS-CoV-2 receptor Tissue-specific and interferon-inducible expression of nonfunctional ACE2 through endogenous retroelement co-option A novel ACE2 isoform is expressed in human respiratory epithelia and is upregulated in response to interferons and RNA respiratory virus infection A gendered magnifying glass on COVID-19 Epigenetic Changes in Asthma: Role of DNA CpG Methylation. Tubercolosis Respir MMP-9 overexpression is associated with intragenic hypermethylation of MMP9 gene in melanoma Defining Driver DNA Methylation Changes in Human Cancer Epigenetic regulation of ACE2, the receptor of the SARS-CoV-2 virus(1) CCAAT/enhancer-binding protein beta overexpression alleviates myocardial remodelling by regulating angiotensin-converting enzyme-2 expression in diabetes Forkhead Box Transcription Factors of the FOXA Class Are Required for Basal Transcription of Angiotensin-Converting Enzyme 2 Dynamic CCAAT/enhancer binding protein-associated changes of DNA methylation in the angiotensinogen gene DNA methylation analysis of the COVID-19 host cell receptor, angiotensin I converting enzyme 2 gene (ACE2) in the respiratory system reveal age and gender differences DNA methylation architecture of the ACE2 gene in nasal cells of children The nasal methylome and childhood atopic asthma Colorectal Cancer that Highly Express Both ACE2 and TMPRSS2, Suggesting Severe Symptoms to SARS-CoV-2 Infection Hepatitis virus infection affects DNA methylation in mice with humanized livers Bacterial infection remodels the DNA methylation landscape of human dendritic cells Blood DNA methylation and COVID-19 outcomes Development and Validation of a Clinical Risk Score to Predict the Occurrence of Critical Illness in Hospitalized Patients With COVID-19 The hallmarks of aging Epigenetic mechanisms in neurological diseases: Genes, syndromes, and therapies Histone Modifications and Cancer ACE2 Expression Is Increased in the Lungs of Patients With Comorbidities Associated With Severe COVID-19 Erasers of histone acetylation: The histone deacetylase enzymes HDACs and HDAC Inhibitors in Cancer Development and Therapy. Cold Spring Harb Epigenetic regulation of angiotensin-converting enzyme 2 (ACE2) by SIRT1 under conditions of cell energy stress Inhibition of SIRT1 deacetylase and p53 activation uncouples the anti-inflammatory and chemopreventive actions of NSAIDs Tissue specific up regulation of ACE2 in rabbit model of atherosclerosis by atorvastatin: Role of epigenetic histone modifications Histone demethylase KDM5B is a key regulator of genome stability KDM5 histone demethylases repress immune response via suppression of STING Genome-wide CRISPR Screens Reveal Host Factors Critical for SARS-CoV-2 Infection Curcumin inhibits HMGB1 releasing and attenuates concanavalin A-induced hepatitis in mice Curcumin attenuated acute Propionibacterium acnes-induced liver injury through inhibition of HMGB1 expression in mice Curcumin inhibits the lymphangiogenesis of gastric cancer cells by inhibiton of HMGB1/VEGF-D signaling Curcumin Attenuates Pulmonary Inflammation in Lipopolysaccharide Induced Acute Lung Injury in Neonatal Rat Model by Activating Peroxisome Proliferator-Activated Receptor gamma (PPARgamma) Pathway Targeting transcriptional regulation of SARS-CoV-2 entry factors ACE2 and TMPRSS2 Targeting androgen regulation of TMPRSS2 and ACE2 as a therapeutic strategy to combat COVID-19 Neutrophil extracellular traps in COVID-19 Human PAD4 regulates histone arginine methylation levels via demethylimination Histone hypercitrullination mediates chromatin decondensation and neutrophil extracellular trap formation PAD4 is essential for antibacterial innate immunity mediated by neutrophil extracellular traps Neutrophil extracellular traps kill bacteria Identification of virus-encoded microRNAs Ambros, V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14 A microRNA controlling left/right neuronal asymmetry in Caenorhabditis elegans MicroRNAs modulate hematopoietic lineage differentiation The Drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila Identification of Modulated MicroRNAs Associated with Breast Cancer, Diet, and Physical Activity MicroRNAs in Neurodegenerative Diseases MicroRNAs in cardiovascular disease: An introduction for clinicians Contribution of MicroRNAs to autoimmune diseases MicroRNome analysis unravels the molecular basis of SARS infection in bronchoalveolar stem cells Small molecule modifiers of microRNA miR-122 function for the treatment of hepatitis C virus infection and hepatocellular carcinoma Treatment of HCV infection by targeting microRNA MicroRNAs targeting the SARS-CoV-2 entry receptor ACE2 in cardiomyocytes miRDB: An online database for prediction of functional microRNA targets Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data miR-98 Regulates TMPRSS2 Expression in Human Endothelial Cells: Key Implications for COVID-19 Protein Coding and Long Noncoding RNA (lncRNA) Transcriptional Landscape in SARS-CoV-2 Infected Bronchial Epithelial Cells Highlight a Role for Interferon and Inflammatory Response Unravelling host-pathogen interactions: ceRNA network in SARS-CoV-2 infection (COVID-19) Role of circRNAs in viral infection and their significance for diagnosis and treatment (Review) Gene regulation by long non-coding RNAs and its biological functions Regulation of circRNA biogenesis Current and innovative methods for the diagnosis of COVID19 infection (Review) COVID-19 Pandemic: From Molecular Biology, Pathogenesis, Detection, and Treatment to Global Societal Impact Molecular immune pathogenesis and diagnosis of COVID-19 A new coronavirus associated with human respiratory disease in China A genome-wide CRISPR screen identifies host factors that regulate SARS-CoV-2 entry Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells Chromatin conformation regulates the coordination between DNA replication and transcription Substantial histone reduction modulates genomewide nucleosomal occupancy and global transcriptional output A core viral protein binds host nucleosomes to sequester immune danger signals Mechanisms of action and regulation of ATP-dependent chromatin-remodelling complexes SWI/SNF nucleosome remodellers and cancer TMPRSS2 and furin are both essential for proteolytic activation of SARS-CoV-2 in human airway cells TMPRSS2 and TMPRSS4 promote SARS-CoV-2 infection of human small intestinal enterocytes Molecular structure of nucleic acids; a structure for deoxyribose nucleic acid The precision medicine initiative: A new national effort Clinical validation of a novel quantitative assay for the detection of MGMT methylation in glioblastoma patients Zinc-dependent Deacetylase (HDAC) Inhibitors with Different Zinc Binding Groups Cedazuridine/decitabine: From preclinical to clinical development in myeloid malignancies Dysregulation of protein methyltransferases in human cancer: An emerging target class for anticancer therapy Critical roles of non-histone protein lysine methylation in human tumorigenesis The Promise for Histone Methyltransferase Inhibitors for Epigenetic Therapy in Clinical Oncology: A Narrative Review Advances in Epigenetic Cancer Therapeutics The authors show great gratitude to all members of the Hamamoto laboratory. The authors declare no conflict of interest.