key: cord-0738443-stln1ees
authors: Qian, Yu; Li, Yi; Liu, Xinxuan; Yuan, Na; Ma, Jinjie; Zheng, Qiwen; Liu, Fan
title: Evidence for CAT gene being functionally involved in the susceptibility of COVID‐19
date: 2021-03-12
journal: FASEB J
DOI: 10.1096/fj.202100008
sha: d13e6601e787e178abc7bf14d97d53a968f81951
doc_id: 738443
cord_uid: stln1ees

Novel coronary pneumonia (COVID‐19) is a respiratory distress syndrome caused by a new type of coronavirus. Understanding the genetic basis of susceptibility and prognosis to COVID‐19 is of great significance to disease prevention, molecular typing, prognosis, and treatment. However, so far, there have been only two genome‐wide association studies (GWASs) on the susceptibility of COVID‐19. Starting with these reported DNA variants, we found the genes regulated by these variants through cis‐eQTL and cis‐meQTL acting. We further did a series of bioinformatics analysis on these potential risk genes. The analysis shows that the genetic variants on EHF regulate the expression of its neighbor CAT gene via cis‐eQTL. There was significant evidence that CAT and the SARS‐CoV‐2‐related S protein binding protein ACE2 interact with each other. Intracellular localization results showed that CAT and ACE2 proteins both exists in the cell membrane and extracellular area and their interaction could have an impact on the cell invasion ability of S protein. In addition, the expression of these three genes showed a significant positive correlation in the lungs. Based on these results, we propose that CAT plays a crucial intermediary role in binding effectiveness of ACE2, thereby affecting the susceptibility to COVID‐19.

these reported DNA variants, we found the genes regulated by these variants through cis-eQTL and cis-meQTL acting. We further did a series of bioinformatics analysis on these potential risk genes. The analysis shows that the genetic variants on EHF regulate the expression of its neighbor CAT gene via cis-eQTL. There was significant evidence that CAT and the SARS-CoV-2-related S protein binding protein ACE2 interact with each other. Intracellular localization results showed that CAT and ACE2 proteins both exists in the cell membrane and extracellular area and their interaction could have an impact on the cell invasion ability of S protein. In addition, the expression of these three genes showed a significant positive correlation in the lungs. Based on these results, we propose that CAT plays a crucial intermediary role in binding effectiveness of ACE2, thereby affecting the susceptibility to COVID-19.

COVID-19, bioinformatics, protein interaction, candidate gene, eQTL This phenomenon of population differences strongly suggests that the susceptibility of COVID-19 is influenced by genetic factors. Understanding the genetic basis of susceptibility to COVID-19 is of great significance to disease prevention, molecular typing, prognosis and treatment.

However, so far, there have been only two genome-wide association studies (GWASs) on COVID-19. 5, 6 One of them by Kachuri et al compared the genotypes of 676 individuals with positive SARS-CoV-2 tests and 1334 with negative tests at the genome-wide level, and found the genetic variation on the EHF gene showing significant association with the tests results. Another GWAS by Ellinghaus et al considered 1610 severe patients as cases and 2205 mild patients or healthy people as controls, and found that genetic variants in the LZTFL1 and ABO genes were significantly associated with the severity of COVID-19. The designs of the two studies were in fact very different in that the first one by Kachuri et al focused on the susceptibility to COVID-19 and the second one by Ellinghaus et al focused on the prognosis of COVID-19. Although these two studies were not very strict in the selection of control samples, the positive findings have provided an important reference for our understanding of COVID-19 susceptibility and prognosis. As for the mechanism by which these positive locus affect the susceptibility and prognosis of COVID-19, it is yet unclear. We next bioinformatically analyzed the potential risk genes associated with COVID-19 susceptibility and prognosis.

We found the genes regulated by these variants through cis-eQTL and cis-meQTL acting. We further performed protein interaction network analysis, intracellular location analysis, and gene expression correlation analysis for these potentially functionally related genes. The results showed that the genetic variants on EHF regulate the expression of its neighbor CAT gene via cis-eQTL acting. There was significant evidence that CAT and the SARS-CoV-2-related S protein binding protein interact with each other. Intracellular localization results showed that CAT protein and ACE2 both exists in the cell membrane and extracellular area. In addition, the expression of these three genes showed a significant positive correlation in the lungs. Based on these results, we propose a hypothesis that CAT plays a crucial intermediary role in binding effectiveness of ACE2, thereby affecting the susceptibility to COVID-19.

We started with the functional annotations of the positive findings of the two GWASs. These included the EHF rs286914 significantly associated with the susceptibility of COVID-19 6 and the LZTFL1 rs11385942 and ABO rs657152 significantly associated with the prognosis of COVID-19. 5 Although one of the GWAS 6 is a preprint, it also laid foundation for follow-up research. The rs286914 is an intronic variant of EHF on chromosome 11 and the A allele of this SNP was associated with an increased risk of positive respiratory virus test (OR = 1.52, Table S1 ). The A allele from the 1000 genomes project showed a large frequency difference between major continents with patterns showing higher frequencies in Africa (AFR 38.20%), America (AME 33.72%), and Europe (EUR 30.82%), but lower frequencies in East Asia (EAS 12.20% ) and South Asia (SAS 15.24%, Figure 1A ). Kachuri GWAS were conducted in UK Biobank participants. As a result, the genomic variants with high frequency in these populations are more likely to be detected. In this case, the variant happened to be rarer in Asians. However, this is merely a reflection of the Winner's curse. We further investigated the existing eQTL 7 and meQTL databases, 8 and found that rs286914 (EHF) is a cis-eQTL of the CAT gene about 140kbp upstream and a cis-meQTL of cg18414381 on the EHF gene in peripheral blood (Tables 1  and 2 , N eQTL = 5311, N meQTL = 4170). Rs286914 also may be the eQTL of CAT in kidney tissue (P = .08, http://mulin lab.org/qtlba se/). This insignificant result might be partially explained by an insufficient sample size (N = 166). Until now, no significant correlation between rs286914 and CAT has been observed in other tissues including the lung. In the future, the sample size needs to be expanded to further verify their correlation in other tissues, especially in lung. This finding suggested a functional role of rs286914 on the expression of CAT while it does not necessarily directly affect the expression of EHF because rs286914 is not a significant cis-eQTL of EHF and cis-meQTL is very common over the genome. 9, 10 We reviewed literature for evidence supporting functional links between these two genes (CAT and EHF) and respiratory diseases. CAT is downregulated in lung cancer 11 and asthma, 12 and promote the apoptosis of non-small cell lung cancer cells by accelerating the degradation of caveolin-1. 13 At the same time, CAT had a protective effect on pulmonary fibrosis 14 and protected lung epithelial cells from hydrogen peroxide-induced apoptosis, 15 supporting that CAT plays an important role in the development of respiratory system diseases. CAT also has a significant anti-inflammatory effect that regulates the production of cytokines in white blood cells, thereby protecting alveolar cells from oxidative damage and inhibiting the replication of SARS-CoV-2, 16 emphasizing its potential utility in the treatment of COVID-19 or other severe inflammation-related diseases in respiratory system. The observation that the EHF variant rs286914 is an eQTL of CAT, together with the lack of literature support for a functional role of EHF in lung diseases including COVID-19, further supports our hypothesis that rs286914 on the EHF may affect COVID-19 via regulatory effects on the expression of CAT.

The GWAS study of Ellinghaus et al found that the individuals carrying the ABO rs657152 A allele showed more Table S1 ). 5 There is no obvious difference in the frequency of rs657152 between continents ( Figure 1B ). In addition, we found that ABO rs657152 is a significant cis-eQTL regulating the expression of GBGT1 in peripheral blood, located about 80 kb upstream of ABO (P = 5.72e-7, Table 1 ). F I G U R E 1 A, Global genotype frequency distribution of rs286914 in EHF; (B) Global genotype frequency distribution of rs657152 in ABO; (C): Global genotype frequency distribution of rs11385942 in LZTFL1; (D) Gene interaction networks created in GeneMania for EHF and CAT. The query genes was analyzed (striped circles) and additionally automatically generated their interacting genes (non-striped circles). The color of the lines connecting the genes denotes the interaction types (supplementary materials). The size of the circle indicates the importance of the gene in the specific interactions while the width of lines indicates the weight of interaction between genes (Table S2 ). E, Correlation of CAT, ACE2, and EHF gene expression among top 6 tissues from the GTEx database (Table S6 ) Meanwhile, rs657152, as a cis-meQTL in peripheral blood, is significantly associated with the methylation level of 23 nearby methylation sites (Table 2 ). This region covers multiple genes (ABO, C9orf7, SURF4, RALGDS, GBGT1 and OBP2B). So far, there is no literature evidence supporting that these genes play an important role in respiratory system diseases. Besides, rs657152 was not observed as eQTL or meQTL in lung. Ellinghaus et al also found that rs11385942, located in the intron region of the LZTFL1, was significantly associated with the prognosis of COVID-19 5 (OR = 2.11, Table S1 ). The risk allele from the 1000 genomes project showed a large frequency difference between major continents with patterns showing higher frequencies in South Asia (SAS 29.55%), but lower frequencies in Europe (EUR 8.05%), Africa (AFR 5.30%), America (AME 4.61%), and East Asia (EAS 0.5%, Figure 1C ). The rs11385942 was not a significant eQTL or meQTL in the peripheral blood and lung. Given the existing literature and databases, we have not sorted out a clear chain of evidence explaining the observed genetic associations (rs657152 and rs11385942) with the prognosis of the COVD-19.

Next, we used GeneMANIA 17 to investigate the protein-protein interactions between EHF and CAT. This analysis identified a functional network consisting of 22 genes. Interestingly, the ACE2 gene, encoding the SARS-CoV-2 S protein receptor, 18 was identified in this network having a protein level interaction with CAT ( Figure 1D , Table S2 ). Noteworthy, among the 163 599 genes from 9 organisms, ACE2 had direct protein-protein interactions with only 6 genes, including CAT ( Figure S1 , Table S3 ). Intracellular localization analysis 19, 20 showed that CAT and ACE2 proteins both exist in the cell membrane and extracellular area which make it clear that the interaction could have an impact on the cell invasion ability of S protein (Table S5) . No direct protein level interaction was found between EHF and ACE2. This result further support our hypothesis that rs286914 (EHF) as an eQTL of CAT may affect the binding efficiency of ACE2 and T A B L E 2 Cis-meQTL of SNPs associated with COVID-19 susceptibility and prognosis in peripheral blood SARS-CoV-2 S protein through CAT, thereby affecting the susceptibility of COVID-19.

There is no protein-protein interaction between the ABO and other genes where the COVID-19 prognosis-associated variants and methylation site are located (GBGT1, C9orf7, SURF4, RALGDS, OBP2B). As illustrated in the network diagram, several genes had the same protein domain, possibly performing same functions ( Figure S2 , Table S4 ). It remains unclear why blood type A individuals had more severe COVID-19 prognosis than those with blood type O. 21 We performed a correlation analysis for the expression of EHF, CAT and ACE2 based on the gene expression data of 49 tissues from the GTEx database (https:// commo nfund.nih.gov/GTEx/) 22 . These 3 genes were simultaneously expressed in 40 different tissues including the lung (Table S6 ). In 30 different tissues including the lung, the expression levels of these 3 genes showed significant correlations for at least one pairs after Bonferroni correction (P < 4.16e-4). The majority of these correlations were positive (80%). The most significant correlation was observed in transverse colon, where the expressions of the three genes were all strongly positively correlated (0.60 < r < 0.73, 2.37e-63 < P < 8.47e-38, Figure 1E ). The 2th significant finding was observed in lung tissue, where the expression of the three genes all showed a significant positive correlation (0.31 < r < 0.49, 1.76e-32 < P < 3.24e-13), which was consistent with that observed in the transverse colon, although the correlation and significance levels were reduced. These correlation patterns indicate that there are intrinsic connections within these three genes and they may regulate the progression of lung-related diseases together. Here, a colocalization analysis would be preferred to test whether the GWAS and eQTL signals are overlapping. However, this analysis could not be done due to the failure to obtain COVID-19 GWAS data. When the GWAS data are available in the future, this analysis can be supplemented to verify the association between CAT and COVID-19.

In summary, EHF rs286914 functionally regulates the expression of CAT via cis-eQTL acting. The expressions of EHF, CAT, and ACE2 are positively correlated with each other in lung. Moreover, there was a significant proteinprotein interaction between CAT and ACE2 while no significant interactions were found between EHF and ACE2. Intracellular location analysis showed that both CAT protein and ACE2 exist in cell membrane and extracellular area and their interaction could have an impact on the cell invasion ability of S protein. These multiple lines of evidence support the hypothesis that EHF may as an intermediary to affect the binding efficiency of ACE2 to SARS-CoV-2 S protein through CAT, thereby affecting the susceptibility of COVID-19.

A novel coronavirus from patients with pneumonia in China

COVID-19: epidemiology, evolution, and cross-disciplinary perspectives

An overview of COVID-19

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

Genomewide association study of severe COVID-19 with respiratory failure

The landscape of host genetic factors involved in infection to common viruses and SARS-CoV-2

Systematic identification of trans eQTLs as putative drivers of known disease associations

Genome-wide identification of DNA methylation QTLs in whole blood highlights pathways for cardiovascular disease

Pancan-meQTL: a database to systematically evaluate the effects of genetic variants on methylation in human cancer

High density methylation QTL analysis in human blood via next-generation sequencing of the methylated genomic DNA fraction

An immunohistochemical analysis of antioxidant and glutathione S-transferase enzyme levels in normal and neoplastic human lung

Nitrotyrosine proteome survey in asthma identifies oxidative mechanism of catalase inactivation

Hydrogen peroxide inhibits non-small cell lung anoikis through the inhibition of caveolin-1 degradation

The role of catalase in pulmonary fibrosis

Mitochondrial localization of catalase provides optimal protection from H2O2-induced cell death in lung epithelial cells

An antioxidant enzyme therapeutic for COVID-19

The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function

Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV

UniProt: the universal protein knowledgebase

UniProt: a hub for protein information

Blood type and outcomes in patients with COVID-19

The genotype-tissue expression (GTEx) project

Program of Chinese Academy of Sciences, Grant No. XDB38010400.

The authors declare that no conflicts of interest exist.

Fan Liu and Yu Qian conceived the study; Yu Qian analyzed the data, prepared the figures and supporting information. All authors took part in writing the article.

Additional supporting information may be found online in the Supporting Information section.