Association of Asthma and Rhinitis with Epigenetics of Coronavirus Related Genes

Abstract

Introduction:

Susceptibility factors for coronavirus disease 2019 (COVID-19) include sex and medical conditions such as asthma and rhinitis. DNA methylation (DNAm) is associated with asthma, rhinitis, and several viruses. We examined associations of asthma/rhinitis with DNAm at CpGs located on coronavirus related genes, and if these associations were sex-specific.

Methods:

In total, n = 242 subjects aged 26 years from the Isle of Wight Birth Cohort were included in the study. Linear regressions were used to examine sex specific and non-specific associations of DNAm at CpGs on coronavirus related genes with asthma/rhinitis status. Associations of DNAm with gene expression in blood were assessed for functional relevance of identified CpGs.

Results:

Statistically significant interaction effects of asthma or rhinitis with sex were identified at 40 CpGs for asthma and 27 CpGs for rhinitis. At 21 CpGs, DNAm was associated with asthma, and at 45 CpGs with rhinitis, regardless of sex. Assessment of functional relevance of the identified CpGs indicated a potential of epigenetic regulatory functionality on gene activity at 14 CpGs for asthma and 17 CpGs for rhinitis, and of those 6 CpGs for asthma and 7 CpGs for rhinitis were likely to be sex-specific.

Conclusion:

Subjects with asthma/rhinitis may have altered susceptibility to COVID-19 due to changes in their DNAm associated with these conditions. Sex specificity on association of asthma/rhinitis with DNAm at certain CpGs, and on the association of DNAm at asthma/rhinitis-linked CpGs with gene expression have the potential to explain the reported sex-specificity in COVID-19 morbidity and mortality.

Keywords

COVID-19 asthma rhinitis DNA methylation IOWBC sex-specificity gene expression

Introduction

The novel coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to spread across the world. As of September 8th, 2020, the pandemic has reached 27 236 916 confirmed cases and 891 031 deaths globally.¹ Among those infected, biological sex was found to be an important determinant of COVID-19 severity and mortality.² Males with COVID-19 are at higher risk for worse outcomes and mortality compared to females, independent of age.³ Existing comorbidities and medical conditions of the respiratory tract are risk factors for COVID-19 morbidity and mortality.^4-7 Asthma is a common chronic condition and rhinitis is a common childhood allergic condition that continues in adulthood,⁸ but their role in COVID-19 infections remains unclear. However, recently findings have shown potential associations of pre-existing asthma and rhinitis with SARS-CoV-2 susceptibility and COVID-19 prognosis.^9-12

Centers for Disease Control and Prevention (CDC) states that chronic lung diseases such as moderate to severe asthma increases the likelihood of COVID-19 severity.¹³ Severe asthma defined as asthma with recent use of an oral corticosteroid has been shown to be associated with higher risk of COVID-19 and its related deaths.¹⁴ A case study of a 3-year-old child has been reported with history of allergic rhinitis suffered from COVID-19.¹⁵ However, some studies report low to no prevalence of asthma in COVID-19 patients.^16,17 Underlying reasons for the discrepancy was not fully understood, but 1 possible explanation for asthma being protective against COVID-19 is that for asthmatic subjects on inhaled corticosteroids, such treatment can potentially protect them against COVID-19.

An epigenetic mechanism, DNA methylation (DNAm) has been shown to be associated with asthma and rhinitis in several studies.^18-22 Previous studies have shown association between DNAm and viruses such as human papillomavirus (HPV), Epstein-Barr virus (EBV), Human rhinovirus.^23-25 These viruses modulate epigenetic dysregulation of immune-related gene expression in host cells leading to immune suppression.²⁶

In this study, we examined whether asthma and rhinitis status were associated with DNAm at CpGs located on coronavirus related genes, and whether such associations were sex specific. The identified CpGs have a potential to be further assessed on their association with the risk of COVID-19 and may contribute to a better understanding on the impact of underlying asthma or rhinitis on COVID-19 susceptibility and severity.

Methods

Study population

Isle of Wight whole population birth cohort (IOWBC, n = 1456) was established in 1989/1990 on the Isle of Wight (UK) to study the natural history of allergic diseases.²⁷ Informed consent was provided by the parents at birth and at follow-up assessments at ages 1, 2, 4, 10 years, and by participants at 18, and 26 years of age.

Clinical and demographic variables

Status of asthma and rhinitis at age 26 years were included in the study. Asthma at 26 years was defined using ISAAC questionnaire as “ever had asthma” and “wheezing or whistling in the chest in the last 12 months” and/or “current treatment for asthma.” Rhinitis at 26 years was defined as “In the past 12 months have you had a problem with sneezing or a runny or a blocked nose when you did not have a cold or the flu?.”²⁸ Sex information was based on the questionnaire data and hospital records collected at the time of birth. Active smoking at 26 years was recorded as “yes” if the subject was a current smoker. Body mass index (BMI) was calculated using weight in kilograms divided by height in meters-squared at 26 years. Socio-economic status (SES) was defined based on household income, number of rooms, and maternal education.

Coronavirus related genes and CpGs

Using keywords “Coronavirus” and “Coronavirus silent sweep” in GeneCards (www.genecards.org/),²⁹ genes potentially related to the pathogenesis of SARS-CoV-2 were extracted. A score for each gene (number of times shown to be associated with coronavirus in previous studies) was used to select genes for this study. A scree plot of the scores in descending order was applied to select genes; the minimum score for a gene to be included in our study is such that its score is at a turning point with a large decrease followed by a flat pattern of scores (Rathod, Rathod and Rahimabad, 2021).³⁰ Information regarding the coronavirus associated CpGs on these selected genes, chromosomes, the location relative to CpG islands, and location on genes were identified using the Illumina’s manifestation file.

DNA methylation (DNAm)

DNA was extracted from peripheral whole blood samples using the standard salting out procedure at age 26 years.³¹ DNA concentration was determined by PicoGreen dsDNA quantitation (Molecular Probes, INC. OR, USA) or Qubit (Thermofisher, MA, USA). One microgram of DNA was bisulfite-treated for cytosine to thymine conversion using the EZ 96-DNA Methylation Kit (Zymo Research, Irvine, CA, USA) for each sample. The methylation level for each CpG was assessed using the Infinium MethylationEPIC BeadChips (Illumina, Inc., CA, USA). Arrays were processed using a standard protocol³² with multiple identical control samples assigned to each bisulfite conversion batch to assess assay variability. Using Illumina BeadStation, the beadchips were scanned and the methylation level (beta value) was calculated for each queried CpG locus using the Methylation module of BeadStudio software.

Intensities of DNAm were quantile normalized using the minfi package in R.³³ The intensities were then converted to β values, defined as a ratio of methylated (M) over the sum of methylated and unmethylated (U) probes (β = M/[c + M + U]), where c is a constant to prevent zero in the denominator. Batch effects correction using the ComBat R function was applied to the β values. β values close to 0 or 1 tend to suffer from severe heteroscedasticity. Therefore, base-2 logit transformed β values (denoted as M-values) were used in the analysis.³⁴ In addition, CpG sites were excluded from the analyses if DNAm was missing for at least 1 subject due to low quality or missing at random, or if the minor allele frequency of a probe SNP at that site is >7% (ie, ~⩾10 out of 1456 subjects expected to have the minor allele in the cohort) and the probe SNP was within 10 base pairs of the targeted CpG site.

As whole blood is a mixture of various cell types,³⁵ there is a need to adjust for cell composition to remove the potential confounding effects.³⁶ Cell type proportions were estimated using the minfi package.^37,38 The estimated cell type proportions of CD4⁺ T cells, natural killer cells, neutrophil, B cells, monocytes, and eosinophil cells were included in the analyses as adjusting factors.

Gene expression

Gene expression levels at 26 years from peripheral blood samples of subjects in IOWBC were measured using paired-end (2 × 75 bp) RNA sequencing with the Illumina Tru-Seq Stranded mRNA Library Preparation Kit with IDT for Illumina Unique Dual Index (UDI) barcode primers following manufacturer’s recommendations. All samples were sequenced twice using the identical protocol and for each sample the output from both runs were combined. FASTQC were run to assess the quality of the FASTQ files (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Reads were mapped against Human Genome (GRch37 version 75) using HISAT2 (v2.1.0) aligner.³⁹ The alignment files, produced in the Sequence Alignment Map (SAM) format, were converted into the Binary Alignment Map (BAM) format using SAMtools (v1.3.1).⁴⁰ HTseq (v0.11.1) was used to count the number of reads mapped to each gene in the same reference genome used for alignment.⁴¹ Normalized read count FPKM (Fragments Per Kilobase of transcript per Million mapped reads) were calculated using the countToFPKM package (https://github.com/AAlhendi1707/countToFPKM) and their log transformed values were used for data analysis.

Statistical analysis

Linear regressions were applied to evaluate the association of DNAm at each CpG site on the coronavirus-related genes (dependent variable) with asthma status (asthma-free as reference) at age 26 years. The same analytical model was applied to the status of rhinitis (rhinitis-free as reference). As previous studies have shown sex differences in asthma prevalence^42-44 as well as COVID-19 severity and mortality,^3,45-47 interaction effects of asthma/rhinitis, and sex on DNAm were assessed. Subjects with asthma free and males are the reference groups.

Association between DNAm and gene expression

To assess the biological relevance of the identified CpGs showing an association between DNAm and asthma or rhinitis status, we further assessed their association with gene expression. Linear regressions were applied to evaluate the association of DNAm levels with gene expression for the CpG mapped genes as well as the nearby genes on the same chromosome. When selecting nearby genes, we considered 500 000 base pairs (bps) (250 000 bps up stream and 250 000 bps downstream of a CpG site) following suggestion in a recent study.¹⁸ Gene expression was used as dependent variable, DNAm and sex (male = reference group) as independent variables in the model. Additionally, to assess sex-specificity of DNAm and expression association, an interaction term of DNAm × sex was also included in the model as we previously found this association to be different in both sexes (Patel, Solatikia and Zhang, 2021).⁵² Statistical significance was considered at α = .05 level. All the analyses were performed using SAS 9.4.

Sensitivity analysis

Of the CpGs that were found to be associated with asthma and rhinitis (interaction and main effects, they were further assessed using the same analytical model as above with covariates BMI, active smoking, and SES at age 26 years included in the model. We would examine whether the interaction effects of asthma or rhinitis with sex were changed in terms of direction of associations after including these covariates in the model.

Results

Following the elbow method of a scree plot on gene scores of >2.8, 66 genes were selected shown to be associated with coronavirus in the literature. Additionally, 14 genes shown to be associated with coronavirus silent sweep infection were also included. Thus, in total 80 genes were selected for the purpose of this study (Supplemental Table 1). Of the 80 genes, we had data for 78 genes on autosomes. In total, 1555 CpGs mapped to 78 genes were extracted from Illumina manifestation file. Of these, 4 CpGs with “ch” prefix were excluded. After preprocessing and exclusion of CpG sites with missingness, low quality, or close to probe SNPs, DNAm at 677 CpGs were available at 26 years in IOWBC (n = 242). Among these 242 subjects, 41 (females = 26, 17.3%) had asthma and the remaining without asthma (females = 124, 82.7%), and 120 (females = 78, 52%) had rhinitis and others without rhinitis (females = 72, 48%). We also examined the difference in asthma prevalence between males and females (Supplemental Tables 4(A) and 4(B)). In the complete cohort, the difference in prevalence between males and females was marginally significant (P-value = .08) for asthma, and was statistically significant for rhinitis (P-value < .001). For both health conditions, prevalence was higher in females. For the subsamples included in the current study, the prevalence in females was also higher, but the difference between the 2 sexes was not statistically significant.

We assessed the association of asthma/rhinitis and DNAm at each of the 677 CpGs, and whether such an association was sex-specific using linear regressions. Two separate models were tested for. Model 1 assessed interaction of asthma and sex while model 2 assessed interaction of rhinitis and sex. Of the 677 CpGs, statistically significant interaction effects with sex were observed at 40 CpGs with asthma (26 genes) (Supplemental Table 2(A)) and 27 CpGs (17 genes) with rhinitis (Supplemental Table 3(A)). A potential sex-specific association of asthma/rhinitis with DNAm was shown at 36 CpGs for asthma (Figure 1a and Supplemental Table 2(A)) and 25 CpGs for rhinitis (Figure 1b and Supplemental Table 3(A)). CpG sites cg05927401 and cg25866075 were identified in both the analyses of asthma and of rhinitis, and at each of these 2 CpGs, the association of asthma/rhinitis with DNAm was sex-specific.

Figure 1.

(a) Effects of asthma on 40 CpGs of coronavirus genes, stratified by sex-specific effects. At 36 CpGs, opposite effects of asthma on DNAm were observed in males and females. (The Y-axis is for regression coefficients of asthma status. Negative values indicate that on average DNAm was lower in asthmatic subjects compared to subjects without asthma. The X-axis lists the CpG with their mapped gene names.) (b) Effects of rhinitis on 27 CpGs of coronavirus genes, stratified by sex-specific effects. At 25 CpGs, opposite effects of rhinitis on DNAm were observed in males and females. (The Y-axis is for regression coefficients of rhinitis status. Negative values indicate that on average DNAm was lower in rhinitis subjects compared to subjects without rhinitis. The X-axis lists the CpG with their mapped gene names.)

For CpGs that we did not observe significant interaction effects between sex and asthma/rhinitis status, we examined the main effects of asthma/rhinitis on DNAm at these CpGs. In total, 637 CpGs for asthma and 650 CpGs for rhinitis were included in the analyses. Statistically significant associations of DNAm were observed at 21 CpGs with asthma (17 genes) (Supplemental Table 2(B)) and at 45 CpGs with rhinitis (31 genes) (Supplemental Table 3(B)). One CpG site, cg16628205, was common between the 21 and 45 identified CpGs.

Results from sensitivity analyses indicated factors such as BMI, smoke exposure, and SES did not substantially affect the results. Of the 40 identified CpGs showing asthma and sex interaction effects, statistical significance remained at 31 CpGs (~80%) with the same direction of associations as in the model without these covariates (Supplemental Table 5(A)). Of the 21 CpGs with statistically significant main effects of asthma, after including the covariates, at 18 CpGs (~90%), statistical significance was intact, and the directions of association were the same as those when covariates were not included (Supplemental Table 5(B)). For rhinitis, similar results were obtained; interaction effects at 21 (80%) of the 27 identified CpGs were intact (Supplemental Table 6(A)), and main effects at 38 (~84%) of the 45 CpGs were not affected (Supplemental Table 6(B)).

To assess the biological relevance of 61 (=40 with sex-interactions +21 with main effects only) and 72 (=27 with sex-interactions +45 with main effects only) CpGs that showed association with asthma and rhinitis respectively, we evaluated the association of DNAm at these CpGs with the expression of genes that the CpGs were mapped to as well as their nearby genes by considering a window of 500 000 base pairs (bps) of the CpG site (250 000 bps up and down stream). Additionally, we also assessed whether such associations were sex-specific. We did not have gene expression data for gene CRP (from the study on asthma and rhinitis) and TMPRSS2 (from the study on rhinitis). For asthma, statistically significant interaction effects were observed at 6 of the 61 CpG sites on their association with expression of 5 genes (Table 1), and at 8 CpG sites, DNAm was associated with expression of 8 genes regardless of sex (Table 2). For rhinitis, statistically significant interaction effects were observed at 6 CpG sites on their association with expression of 7 genes (Table 3), and at 11 CpG sites, DNAm was associated with expression of 8 genes regardless of sex (Table 4). In particular, an increase in DNAm was associated with downregulation of 5 genes in males and 4 genes in females for asthma, indicated by negative regression coefficients for each sex. While for rhinitis, an increase in DNAm was associated with downregulation of 4 genes in males and 5 genes in females. For the sex non-specific associations, an increase in DNAm was associated with downregulation of 7 genes for asthma and 6 genes for rhinitis. For the remaining genes, an upregulating relationship was identified. The flowchart of the study is shown in figure 2.

Table 1.

Asthma: sex-specific association of DNAm at 6 CpGs with expression of their mapped genes as well as nearby genes (male = reference group).

CpG site	Gene	Int.	Est.	P value*	Chr.	CpG islands	Gene location
cg12455187	CCL5	−1.75	1.68	3.57 × 10⁻⁵	17	—	TSS1500
cg01281450	HLA_B	0.69	−0.26	.03	12	—	3′UTR
cg15096505	IL10	2.12	−3.56	.01	1	—	Body; Intron
cg17744604	IL10	1.34	−1.62	.04	1	—	TSS1500
cg19307543	MYOM2	−1.21	−0.37	.02	8	—	Body; Intron
cg06582411	TFR2	0.79	−0.77	.03	7	N_Shore	Body; Intron

Common genes between the findings related to asthma and those on rhinitis (Table 1 and 3) are in bold font.

Abbreviations: Chr., Chromosome; Est., main estimate; Int., interaction coefficient of the asthma with the sex.

P-values for the interaction effects of CpG × sex.

Table 2.

Asthma: sex non-specific association of DNAm at 8 CpGs with expression of their mapped genes as well as nearby genes (male = reference group).

CpG site	Gene	Est.	P value*	Chr.	CpG islands	Gene location
cg08596817	IFITM3	−0.81	.0002	11	S_Shore	TSS200
cg01281450	IFNG	−1.39	8.37 × 10⁻⁵	12	—	3′UTR
cg05495029	IRF3	−0.32	.01	19	Island	TSS200
cg17498313	MAPK14	0.37	.02	6	Island	TSS200
cg01439248	MYOM2	−1.46	.04	8	S_Shore	Body; Intron
cg08643950	MYOM2	−1.44	.03	8	N_Shore	Body
cg13986130	CEACAM1	−0.36	.03	1	Island	Body; Exon
cg13986130	MX1	−0.55	.01	1	Island	—
cg14530382	SH2D3A	0.35	1.63 × 10⁻⁶	19	Island	Body; Exon

Common genes between the findings related to asthma and those on rhinitis (Table 2 and 4) are in bold font.

Abbreviations: Chr., chromosome; Est., main estimate.

P-values for the main effects of CpGs.

Table 3.

Rhinitis: sex-specific association of DNAm at 6 CpGs with expression of their mapped genes as well as nearby genes (male = reference group).

CpG site	Gene	Int.	Est.	P value*	Chr.	CpG islands	Gene location
cg19776453	CEACAM1	2.07	−2.72	0.02	19	—	TSS1500
cg19776453	MX1	2.44	−2.71	0.04	19	—	TSS1500
cg04038163	CXCL9	−1.42	0.81	0.01	4	—	Body; Intron
cg16582889	GPT	2.97	−1.08	0.02	8	N_Shore	TSS1500
cg13902357	HLA_B	0.47	−0.25	0.01	6	Island	Body
cg07821354	SH2D3A	−0.74	0.59	0.04	19	Island	TSS1500
cg09369818	TFR2	−0.54	0.42	0.002	7	S_Shore	Body; Exon

Common genes between the findings related to asthma and those on rhinitis (Table 1 and 3) are in bold font.Abbreviations: Chr., chromosome; Est., main estimate; Int., interaction coefficient of the rhinitis with the sex.

P-values for the interaction effects of CpG × sex.

Table 4.

Rhinitis: sex non-specific association of DNAm at 11 CpGs with expression of their mapped genes as well as nearby genes (male = reference group).

CpG site	Gene	Est.	P value*	Chr.	CpG islands	Gene location
cg22316575	ADA	−0.40	.04	20		3′UTR
cg05626117	CLEC4G	0.74	.0005	19	S_Shore	TSS1500
cg01064373	HLA_B	0.46	.002	6	N_Shore	Body; Exon
cg01064373	HLA_C	1.04	3.37 × 10⁻⁶	6	N_Shore	Body; Exon
cg11867651	HLA_B	0.12	6.1 × 10⁻⁶	6	N_Shore	Body; Intron
cg11867651	HLA_C	0.36	0	6	N_Shore	Body; Intron
cg18511546	HLA_B	0.44	0.0004	6	N_Shore	Body; Intron
cg18511546	HLA_C	1.09	2 × 10⁻⁹	6	N_Shore	Body, Intron
cg06490385	IRF3	−0.21	0.03	19	Island	TSS200
cg05759948	NLRP12	0.48	0.03	10	—	TSS1500
cg05563672	MX1	−0.37	0.02	21	Island	TSS200
cg01360627	TNF	−0.46	4.4 × 10⁻⁶	6	N_Shelf	Body; Exon
cg15989608	TNF	−0.60	6.3 × 10⁻⁶	6	N_Shelf	3′UTR
cg26736341	TNF	−0.61	1 × 10⁻⁵	6	N_Shelf	3′UTR

Common genes between the findings related to asthma and those on rhinitis (Table 2 and 4) are in bold font. Abbreviations: Chr., chromosome; Est., main estimate.

P-values for the main effects of CpGs.

Figure 2.

Flowchart of the study.

Discussion

In this study, we assessed the association of asthma/rhinitis with DNAm at CpGs on genes associated with coronavirus in a sex-specific and non-specific manner. We identified statistically significant interaction effects between sex and asthma/rhinitis at 40 CpGs (26 genes) for asthma and 27 CpGs (17 genes) for rhinitis. In addition, at 21 CpGs (17 genes) with asthma and at 45 CpGs (31 genes) with rhinitis, asthma/rhinitis was associated with DNAm regardless of sex. Our sensitivity analysis indicated that the associations of DNAm at these CpGs with asthma/rhinitis were not substantially impacted by potential confounding factors such as BMI, smoking exposure, or SES. Assessment of the biological relevance on the identified CpGs indicated a potential of epigenetic regulatory functionality on gene activity at some of the CpGs and such epigenetic functionality was likely to be sex specific. In addition, for most genes showing association between DNAm and gene expression, the association was stronger in males (indicated by larger effect sizes for males).

Current findings acknowledge the sex difference in COVID-19 susceptibility.^3,45-47 The prevalence of asthma and rhinitis was different by sex in our complete cohort, which was consistently observed in the study samples (although the difference was statistically insignificant), with females having a higher prevalence. This has been shown in earlier studies (Patel, Solatikia and Zhang, 2021).^48-52 The disproportionality in asthma and rhinitis between males and females possibly contributed to the observed interaction effects of such health conditions with sex on DNAm of genes linked to coronavirus in the present study. Although further investigations are certainly needed to help us better understand such associations, with a focus on CpGs located on genes shown to be associated with coronavirus, our results of sex-specificity on the association of asthma/rhinitis with DNAm and on the association of DNAm with gene expression support sex-specificity in COVID-19 susceptibility and severity.³ Our findings have a potential to help explain the disproportionality of COVID-19 between males and females. In patients with asthma and rhinitis, the DNAm might have been altered and thereby influence the risk of coronavirus susceptibility or severity of manifestations. This suggests that gene expression regulated by epigenetics, for example, DNAm, may play a critical role in a host’s susceptibility of coronavirus.²⁶ Further epidemiological and experimental investigations are warranted to examine the epigenetic functionality of the identified CpGs and assess their role in the occurrence of diseases linked to coronavirus.

At 3 genes, HLA-B, IRF3, and MX1, for both asthma and rhinitis, DNA methylation was associated with expression of genes. All these 3 genes have been shown to play a role in pathogenesis of SARS-CoV-2 infection. First, HLA-B was shown to be associated with the severity of SARS-CoV-2 infection.⁵³ Alleles of HLA-B gene such as HLA-B*07:03, HLA-B*46:01 have been shown to be associated with SARS-CoV-1 susceptibility.⁵⁴ As HLA-B*46:01 is expected to have few possible binding peptides for SARS-CoV-2, individuals with this allele may be vulnerable to COVID-19. While HLA-B*15:03 allele has highly conserved SARS-CoV-2 peptides making individuals with this allele to have cross-protective T-cell dependent immunity.⁵³ Other studies have shown HLA-B alleles HLA-B^*07:03,⁵⁵ HLA-B^*46:01,⁵⁶ and HLA-B*52:01:01:02⁵⁷ to be significantly associated with severity of COVID-19. A significant correlation of HLA-B*27:07 was observed in a group of 99 Italian patients affected by severe or extremely severe course of COVID-19.⁵⁸ Interestingly, DNAm has a quite strong potential of sex-specific regulatory functionality on HLA-B gene (indicated by rather small P-value 6.1 × 10⁻⁶), which deserves further investigation regarding its potential to explain infection severity differentiation between sex. Second, IRF3 inactivates interferon signaling which aids viral replication of SARS-CoV-2 infection.^59-61 Gene expression of IRF3 has been significantly shown to be upregulated in COVID-19 patients by approximately 250-fold compared to uninfected patients.⁶² Another study showed SARS-CoV-2 inhibits activation of IRF3 thereby slowing IFNβ production upon infection.⁶³ Third, the gene MX1, an interferon responsive gene target, was found to be highly upregulated in SARS-CoV-2 infected normal human bronchial epithelial cells suggesting an activation of antiviral interferon innate response.⁶⁴ MX1 antiviral gene expression was significantly higher in COVID-19 patients when compared to non-COVID-19 patients, has been shown to be triggered by SARS-CoV-2, thereby playing a crucial role in inhibiting viral replication.⁶⁵ Minor alleles of 5 SNPs near MX1 were found to be correlated with decreased risk to develop severe COVID-19 and high level of MX1 expression in blood.⁶⁶ MX1 has been suggested as a potential candidate for drug targets in COVID-19 treatment.^65-67 All these 3 genes have been previously demonstrated to be associated with SARS-CoV-2 pathogenesis. Our study from a different angle examined these genes in terms of their role in the connection between asthma or rhinitis and the susceptibility of COVID-19. We showed that DNAm at CpGs on these genes was associated with asthma or rhinitis, which was further associated with expression of genes. Our findings indicated that asthma or rhinitis might have resulted in DNAm changes on these genes and the potential regulatory functionality (shown by their association with gene expressions) of these changes in methylation levels at CpGs of these genes may influence gene expression, thereby influencing the COVID-19 susceptibility and severity among people with asthma and/or rhinitis. We hope our epigenetic findings will help us better understand the pathogenesis of COVID-19 from the epigenetic perspective.

The strength of this study is the availability of DNAm and expression data of genes, shown in the literature to be associated with coronavirus, enabling us to examine the association of DNAm at CpGs on those genes with asthma/rhinitis among adults, and how such association was sex-specific. To our knowledge this is the first study to examine SARS-CoV-2 related genes from an epigenetic perspective.

We aimed to provide a comprehensive epigenetic insight on genes related to coronavirus. This study has some limitations worth mentioning. First, in our data, DNA was extracted from peripheral blood cells. Coronaviruses affect different cell types, primarily, cells of the respiratory tract.⁶⁸ In our analysis, we extracted candidate genes potentially related to SARS-CoV-2 infection, and did not have information on which tissue or cells these genes were identified in. It has been shown that DNAm of blood cells has concordance with that of the respiratory system cells,⁶⁹ although some differences exist between the 2 and tissues specifically related to the respiratory system are ideal. For CpGs with DNAm associated with gene expressions, our findings supported their biological relevance and indicated that those CpGs at least had a potential to serve as markers. However, for genes not showing any association with expression of genes, studies focus on different tissues are desired to examine epigenetic perspectives of asthma and rhinitis on coronavirus related genes. Second, the candidate genes included in our study were not randomly selected or untargeted (eg, as in genome-wide studies). In such situations, forcing multiple testing would potentially increase type II error and is not recommended.⁷⁰ However, we believe it is important to further test the identified CpGs in independent cohorts to assess their reproducibility and generalizability. Third, we were not able to assess the impact of the identified CpGs on the risk of COVID-19 incidence due to lack of COVID-19 data in those subjects. Nevertheless, our findings showed that subjects with asthma or rhinitis had different DNAm compared to people without the diseases on certain genes related to coronavirus, which suggests a possibility that subjects with asthma or rhinitis may be at an altered risk of susceptibility to the SARS-CoV-2 infection due to changes in their DNAm. An assessment on the direction of risk, increased/decreased/intact, requires information on infection status, and certainly deserves further investigation.

Supplemental Material

sj-docx-1-gae-10.1177_25168657211039224 – Supplemental material for Association of Asthma and Rhinitis with Epigenetics of Coronavirus Related Genes

Supplemental material, sj-docx-1-gae-10.1177_25168657211039224 for Association of Asthma and Rhinitis with Epigenetics of Coronavirus Related Genes by Aniruddha Rathod, Rutu Rathod, Hongmei Zhang, Parnian Kheirkhah Rahimabad, Wilfried Karmaus and Hasan Arshad in Epigenetics Insights

Footnotes

Acknowledgements

The authors are thankful to the nurses and staff at the David Hide Asthma and Research Allergy Centre, Isle of Wight, UK, for their help in recruitment and sample collections, and are thankful to all the cohort participants.

Funding:

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by the National Institutes of Health (NIH) research fund R01AI121226 (MPI: H. Zhang and J.W. Holloway) and by NIH research fund R01HL132321 and R01AI091905 (PI: W. Karmaus). The Isle of Wight Birth Cohort assessments were supported by NIH research fund R01HL082925 (PI: H. Arshad), Asthma UK (Grant no. 364. S.H. Arshad), and the David Hide Asthma and Allergy Research Trust.

Declaration of conflicting interests:

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Author Contributions

HZ, WK, and AR conceived the study. HZ, RR, and AR designed the analytical plan. SHA provided data. PR, WK, and HZ searched the genes related to coronavirus. AR and RR did data management and performed the statistical analyses. AR, RR, and HZ drafted the manuscript. All authors reviewed and approved the final manuscript.

Availability of Data and Materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethics Approval and Consent to Participate

The IoW birth cohort study was approved by the Isle of Wight, Portsmouth and Hampshire Local Research Ethics Committee (now known as the National Research Ethics Service, NRES Committee South Central—Hampshire A) (06/Q1701/34) DNA methylation related studies have been approved by the University of Memphis Institutional Review Board in Memphis, U.S. (FWA00006815).

ORCID iDs

Aniruddha Rathod

Hongmei Zhang

Supplemental Material

Supplemental material for this article is available online.

References

World Health Organization. WHO coronavirus disease (COVID-19) dashboard 2020. Accessed June 25, 2020. https://covid19.who.int/

Márquez

Trowbridge

Kuchel

Banchereau

Ucar

. The lethal sex gap: COVID-19. Immun Ageing. 2020;17:13.

Jin

J-M

Bai

, et al. Gender differences in patients with COVID-19: focus on severity and mortality. Front Public Health. 2020;8:152.

Fang

Karakiulakis

Roth

. Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection? Lancet Respir Med. 2020;8:e21.

Wang

, et al. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan, China. JAMA. 2020;323:1061-1069.

Huang

Wang

, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395:497-506.

Yang

Zheng

Gou

, et al. Prevalence of comorbidities and its effects in patients infected with SARS-CoV-2: a systematic review and meta-analysis. Int J Infect Dis. 2020;94:91-95.

Zhang

Kaushal

Soto-Ramírez

, et al. Acquisition, remission, and persistence of eczema, asthma, and rhinitis in children. Clin Exp Allergy. 2018;48:568-576.

Hosoki

Chakraborty

Sur

. Molecular mechanisms and epidemiology of COVID-19 from an allergist’s perspective. J Allergy Clin Immunol. 2020;146:285-299.

10.

Zimmermann

Curtis

. Coronavirus infections in children including COVID-19: an overview of the epidemiology, clinical features, diagnosis, treatment and prevention options in children. Pediatr Infect Dis J. 2020;39:355-368.

11.

Stawicki

Jeanmonod

Miller

, et al. The 2019–2020 novel coronavirus (severe acute respiratory syndrome coronavirus 2) pandemic: a Joint American College of Academic International Medicine-World Academic Council of Emergency Medicine Multidisciplinary COVID-19 Working Group Consensus Paper. J Glob Infect Dis. 2020;12:47-93.

12.

Guo

Pan

, et al. New insights of emerging SARS-CoV-2: epidemiology, etiology, clinical features, clinical treatment, and prevention. Front Cell Dev Biol. 2020;8:410.

13.

Prevention CfDCa. COVID-19: People with moderate to severe asthma 2021. https://www.cdc.gov/coronavirus/2019-ncov/need-extra-precautions/asthma.html. Accessed on June 4, 2021.

14.

Williamson

Walker

Bhaskaran

, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584:430-436.

15.

Dong

Cao

, et al. Eleven faces of coronavirus disease 2019. Allergy. 2020;75:1699-1709.

16.

Chhiba

Patel

THT

, et al. Prevalence and characterization of asthma in hospitalized and nonhospitalized patients with COVID-19. J Allergy Clin Immunol. 2020;146:307-314.e4.

17.

Wang

Xie

. The association between COVID-19 and asthma: a systematic review and meta-analysis. Clin Exp Allergy. 2020;50:1274-1277.

18.

Reese

den Dekker

, et al. Epigenome-wide meta-analysis of DNA methylation and childhood asthma. J Allergy Clin Immunol. 2019;143:2062-2074.

19.

Jiang

Yang

, et al. Nasal DNA methylation profiling of asthma and rhinitis. J Allergy Clin Immunol. 2020;145:1655-1663.

20.

Gao

Huang

Wang

Cai

. Identification of DNA methylation module in seasonal allergic rhinitis. Int J Pediatr Otorhinolaryngol. 2019;117:163-166.

21.

Meng

Wang

Zhang

. Recent developments and highlights in allergic rhinitis. Allergy. 2019;74:2320-2328.

22.

Rathod

Duan

Zhang

, et al. Interweaving between genetic and epigenetic studies on childhood asthma. Epigenet Insights. 2020;13:2516865720923395.

23.

Clarke

Gradissimo

Schiffman

, et al. Human Papillomavirus DNA methylation as a biomarker for cervical precancer: consistency across 12 genotypes and potential impact on management of HPV-positive women. Clin Cancer Res. 2018;24:2194-2202.

24.

Lam

WKJ

Jiang

Chan

KCA

, et al. Methylation analysis of plasma DNA informs etiologies of Epstein-Barr virus-associated diseases. Nat Commun. 2019;10:3256.

25.

Pech

Weckmann

König

, et al. Rhinovirus infections change DNA methylation and mRNA expression in children with asthma. PLoS One. 2018;13:e0205275.

26.

Kuss-Duerkop

Westrich

Pyeon

. DNA tumor virus regulation of host DNA methylation and its implications for immune evasion and oncogenesis. Viruses. 2018;10(2):82.

27.

Arshad

Holloway

Karmaus

, et al. Cohort profile: the Isle of Wight Whole Population Birth Cohort (IOWBC). Int J Epidemiol. 2018;47:1043-1044i.

28.

Ziyab

Karmaus

Zhang

, et al. Association of filaggrin variants with asthma and rhinitis: is eczema or allergic sensitization status an effect modifier? Int Arch Allergy Immunol. 2014;164:308-318.

29.

Stelzer

Rosen

Plaschkes

, et al. The GeneCards suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinform. 2016;54:1.30.1-1.30.33.

30.

Rathod

Rahimabad

Duan

Zhang

Arshad

Karmaus

. Methylation of Host Genes Associated with Coronavirus Infection from Birth to 26 Years. Genes. 2021;12(8):1198. https://doi.org/10.3390/genes12081198

31.

Miller

Dykes

Polesky

. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16:1215.

32.

Bibikova

Fan

. GoldenGate assay for DNA methylation profiling. Methods Mol Biol. 2009;507:149-163.

33.

Aryee

Jaffe

Corrada-Bravo

, et al. Minfi: a flexible and comprehensive bioconductor package for the analysis of infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363-1369.

34.

Zhang

Huang

, et al. Comparison of beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587.

35.

Reinius

Acevedo

Joerink

, et al. Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS One. 2012;7:e41361.

36.

Koestler

Christensen

Karagas

, et al. Blood-based profiles of DNA methylation predict the underlying distribution of cell types: a validation analysis. Epigenetics. 2013;8:816-826.

37.

Jaffe

Irizarry

. Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. 2014;15:R31.

38.

Houseman

Accomando

Koestler

, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.

39.

Kim

Langmead

Salzberg

. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357-360.

40.

Anders

Pyl

Huber

. HTSeq – a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166-169.

41.

Handsaker

Wysoker

, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078-2079.

42.

Shah

Newcomb

. Sex bias in asthma prevalence and pathogenesis. Front Immunol. 2018;9:2997.

43.

Pignataro

Bonini

Forgione

Melandri

Usmani

. Asthma and gender: the female lung. Pharmacol Res. 2017;119:384-390.

44.

Zhang

Zein

. Novel insights on sex-related differences in asthma. Curr Allergy Asthma Rep. 2019;19:44.

45.

Gausman

Langer

. Sex and gender disparities in the COVID-19 pandemic. J Womens Health. 2020;29:465-466.

46.

Gebhard

Regitz-Zagrosek

Neuhauser

Morgan

Klein

. Impact of sex and gender on COVID-19 outcomes in Europe. Biol Sex Differ. 2020;11:29.

47.

Huang

Wang

Liang

Huang

, et al. COVID-19 patients’ clinical characteristics, discharge rate, and fatality rate of meta-analysis. J Med Virol. 2020;92:577-583.

48.

Postma

. Gender differences in asthma development and progression. Gend Med. 2007;S133-S146.

49.

Leynaert

Sunyer

Garcia-Esteban

, et al. Gender differences in prevalence, diagnosis and incidence of allergic and non-allergic asthma: a population-based cohort. Thorax. 2012;67:625-631.

50.

Cazzoletti

Ferrari

Olivieri

, et al. The gender, age and risk factor distribution differs in self-reported allergic and non-allergic rhinitis: a cross-sectional population-based study. Allergy Asthma Clin Immunol. 2015;11:36.

51.

Fröhlich

Pinart

Keller

, et al. Is there a sex-shift in prevalence of allergic rhinitis and comorbid asthma from childhood to adulthood? A meta-analysis. Clin Transl Allergy. 2017;7:44.

52.

Patel

Solatikia

Zhang

, et al. Sex-specific associations of asthma acquisition with changes in DNA methylation during adolescence. Clin Exp Allergy. 2021;51:318-328.

53.

Nguyen

David

Maden

, et al. Human leukocyte antigen susceptibility map for severe acute respiratory syndrome coronavirus 2. J Virol. 2020;94:e00510-20.

54.

Sanchez-Mazas

. HLA studies in the context of coronavirus outbreaks. Swiss Med Wkly. 2020;150:w20248.

55.

Lau

, et al. Association of human-leukocyte-antigen class I (B*0703) and class II (DRB1*0301) genotypes with susceptibility and resistance to the development of severe acute respiratory syndrome. J Infect Dis. 2004;190:515-518.

56.

Lin

Tseng

Trejaut

, et al. Association of HLA class I with severe acute respiratory syndrome coronavirus infection. BMC Med Genet. 2003;4:9.

57.

Khor

S-S

Omae

Nishida

, et al. HLA-A*11:01:01:01, HLA-c*12:02:02:01-HLA-B*52:01:02:02, age and Sex are associated with severity of Japanese COVID-19 with respiratory failure. Front Immunol. 2021;12:658570.

58.

Novelli

Andreani

Biancolella

Liberatoscioli

Passarelli

Colona

, et al. HLA allele frequencies and susceptibility to COVID-19 in a group of 99 Italian patients. HLA. 2020;96:610-614.

59.

Yoshimoto

. The proteins of severe acute respiratory syndrome coronavirus-2 (SARS CoV-2 or n-cOV19), the cause of COVID-19. Protein J. 2020;39:198-216.

60.

Tay

Poh

Rénia

MacAry

LFP

. The trinity of COVID-19: immunity, inflammation and intervention. Nat Rev Immunol. 2020;20:363-374.

61.

Tang

Comish

Kang

. The hallmarks of COVID-19 disease. PLoS Pathog. 2020;16:e1008536.

62.

Salman

Waheed

Ali-Abdulsahib

Atwan

. Low type I interferon response in COVID-19 patients: Interferon response may be a potential treatment for COVID-19. Biomed Rep. 2021;14:43.

63.

Glanz

Chakravarty

Varghese

, et al. Transcriptional and non-transcriptional activation, posttranslational modifications, and antiviral functions of interferon regulatory factor 3 and viral antagonism by the SARS-Coronavirus. Viruses. 2021;13:575.

64.

Vishnubalaji

Shaath

Alajez

. 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. Genes. 2020;11(7):760.

65.

Bizzotto

Sanchis

Abbate

, et al. SARS-CoV-2 infection boosts MX1 antiviral effector in COVID-19 patients. iScience. 2020;23:101585.

66.

Andolfo

Russo

Lasorsa

, et al. Common variants at 21q22.3 locus influence MX1 and TMPRSS2 gene expression and susceptibility to severe COVID-19. iScience. 2021;24:102322.

67.

Prasad

Khatoon

Rashid

, et al. Targeting hub genes and pathways of innate immune response in COVID-19: a network biology perspective. Int J Biol Macromol. 2020;163:1-8.

68.

Hui

KPY

Cheung

Perera

RAPM

, et al. Tropism, replication competence, and innate immune responses of the coronavirus SARS-CoV-2 in human respiratory tract and conjunctiva: an analysis in ex-vivo and in-vitro cultures. Lancet Respir Med. 2020;8:687-695.

69.

Stueve

Shi

, et al. Epigenome-wide analysis of DNA methylation in lung tissue shows concordance with blood studies and identifies tobacco smoke-inducible enhancers. Hum Mol Genet. 2017;26:3014-3027.

70.

Rothman

. No adjustments are needed for multiple comparisons. Epidemiology. 1990;1:43-46.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.07 MB