key: cord-0750485-4ao3ydt7
authors: Chen, Long; Liu, Yue; Wu, Jiamin; Deng, Chao; Tan, Jianjun; Liu, Huawen; Zhong, Li
title: Lung adenocarcinoma patients have higher risk of SARS-CoV-2 infection
date: 2021-01-10
journal: Aging (Albany NY)
DOI: 10.18632/aging.202375
sha: 39b3443fc64def9586dfe790460f8f22b2af8e42
doc_id: 750485
cord_uid: 4ao3ydt7

Both lung adenocarcinoma and coronavirus disease 2019 would cause pulmonary inflammation. Angiotensin-converting enzyme 2, the functional receptor of SARS-CoV-2, also plays a key role in lung adenocarcinoma. To study the risk of SARS-CoV-2 infection in lung adenocarcinoma patients, mRNA and microRNA profiles were obtained from The Cancer Genome Atlas and Gene Expression Omnibus followed by bioinformatics analysis. A network which regards angiotensin-converting enzyme 2 as the center was structured. In addition, via immunological analysis to explore the essential mechanism of SARS-CoV-2 susceptibility in lung adenocarcinoma. Compared with normal tissue, angiotensin-converting enzyme 2 was increased in lung adenocarcinoma patients. Furthermore, a total of 7 correlated differently expressed mRNAs (ACE2, CXCL9, MMP12, IL6, AZU1, FCN3, HYAL1 and IRAK3) and 5 correlated differently expressed microRNAs (miR-125b-5p, miR-9-5p, miR-130b-5p, miR-381-3p and miR-421) were screened. Interestingly, the most frequent toll-like receptor signaling pathway was enriched by mRNA (interlukin 6) and miRNA (miR-125b-5p) sets simultaneously. In conclusion, it was assumed that miR-125b-5p-ACE2-IL6 axis could alter the risk of SARS-CoV-2 infection in lung adenocarcinoma patients.

In December 2019, there was an outbreak of coronavirus pneumonia (COVID-19) caused by SARS-CoV-2 in Wuhan, Hubei province in China. SARS-CoV-2 could encode spike (S) protein like SARS-CoV. In addition, a same cell entry receptor and manner were used to mediate infection. Unlike other type I fusion proteins, the S protein of coronaviruses is not cleaved in the virusproducing cell [1, 2] . However, in processed coronaviruses, two domains (S1 and S2) with different functions can be defined [3] . The S1 domain induces receptor association, whereas the S2 domain likely undergoes structural rearrangements that mediate membrane fusion. A discrete receptor-binding domain (RBD) of S protein has been defined at residues 318-510 of the S1 domain. This RBD could bind receptor with higher affinity than full S1 domain [4, 5] .

Angiotensin-converting enzyme 2 (ACE2), a terminal carboxypeptidase, catalyzes the conversion of angiotensin II (Ang II) to angiotensin 1-7 (Ang 1-7). Ang II, the major effector molecule of the conventional reninangiotensin system, is implicated in the pathogenesis of cardiovascular disorders like hypertension, atherosclerosis, and myocardial infarction. Whereas ACE2 and its product Ang 1-7 are thought to prevent the detrimental effects of angiotensin II [6] . In respiratory organs, ACE2 is the receptor for SARS-CoV [7, 8] . As a type I transmembrane protein, ACE2 is comprised of a short cytoplasmic domain, a transmembrane domain, and a large ectodomain [9] . The ACE2 ectodomain, includes the first α-helix and N-terminus proximal residues of βsheet, could bind the SARS-CoV S glycoprotein with high affinity [10] . Furthermore, once SARS-CoV binds to ACE2, the abundance on the cell surface, the mRNA expression and enzymatic activity of ACE2 are AGING significantly reduced. ACE2 has also been shown to attenuate inflammation and acute lung injury caused by SARS-CoV infection mainly through inhibition of angiotensin II/NF-κB signaling [11] [12] [13] [14] [15] . The ACE2 dysfunction is implicated in SARS pathogenesis. ACE2 is predominantly localized on the apical surface of welldifferentiated airway epithelia, especially ciliated cells, and also been identified with the pathology of various inflammatory lung diseases like lung cancer [16] [17] [18] . Furthermore, the expression of ACE2 might be influenced by many factors like organ, gender, age, development, lifestyle (cigarette smoking, diet) and some diseases (chronic obstructive pulmonary disease, hypertension, coronary heart disease, chronic kidney disease and diabetes) [19] . Abundant clinical data has verified that the ACE2 level would alter infection risk of SARS-cov-2. For example, COVID-19 patients with comorbidities were more likely to progress to critically ill patients compared with those without comorbidities [20] .

Lung adenocarcinoma (LUAD), the most common primary lung cancer which falls under the umbrella of non-small cell lung cancer (NSCLC), has a strong association with smoking. LUAD usually evolves from the mucosal glands and represents about 40% of lung cancers. LUAD usually occurs in the lung periphery, and in many cases, may be found in scars or areas of chronic inflammation [21, 22] . The tumor microenvironment of LUAD may also regulate ACE2 expression level which alters the risk of SARS-CoV-2 infection. Furthermore, there are lots of factors that could cause both ACE2 level alteration and fluctuation of immune-regulation. Thus, through bioinformatics analysis and network structuring, the potential mechanism of COVID-19 in LUAD patient might be explored.

Ang 1-7, the product of ACE2, was confirmed to inhibit the proliferation of human lung cancer cells through interaction with the Mas receptor (MasR) [23] . Ang1-7 not only reduces the size of human lung tumor xenografts in vivo but also markedly decreases their vessel density [24] [25] [26] . Thus, the expression level of ACE2 in LUAD would influence the cancer progression. To compare the ACE2 level between LUAD and normal tissues, sequencing data of a total of 830 samples (347 normal tissues and 483 LUAD tissues) were obtained from the TCGA database ( Figure 1 ). In LUAD patients, a higher expression level of ACE2 was detected while there was no significant difference among various stages of LUAD. Due to the function of combining with S glycoprotein, higher ACE2 level might lead to increased risk of SARS-CoV-2 infection. 

The fluctuation of gene expression often involves complicated regulatory networks. To study the potential mechanism of ACE2 dysregulation in LUAD, the relevant molecular around ACE2 were identified. After screening and comparison, a total of 7 differentially expressed correlative genes (DECGs) and 5 differentially expressed correlative miRNAs (DECMs) were detected respectively. In addition, the DECGs (CXCL9, MMP12, IL6, AZU1, FCN3, HYAL1, IRAK3) were involved in multi defense virus processes and the DECMs (miR-125b-5p, miR-9-5p, miR-130b-5p, miR-381-3p, miR-421) were predicted to target the 3' untranslated regions (UTR) of ACE2 ( Figure 2 and Supplementary Figure 1 ).

Among them, CXCL9, MMP12, miR-9-5p, miR-130b-5p, miR-381-3p and miR-421 were increased while IL6, AZU1, FCN3, HYAL1, IRAK3 and miR-125b-5p were reduced in LUAD patients. Meanwhile, AZU1, FCN3, HYAL1, IRAK3 and miR-125b-5p were positively correlated with ACE2 while CXCL9, MMP12, IL6, miR-9-5p, miR-130b-5p, miR-381-3p, and miR-421 were negatively correlated ( Figure 3) . Notably, only miR-125b-5p has an opposite expression trend with ACE2. Thus, miR-125b-5p was assumed to be the single upstream predicted inhibitor of ACE2. In addition, among the 7 DECGs, IL6 owns the highest mutation frequency (5%) in LUAD patients might indicate that IL6 could mediate carcinogenesis (Supplementary Figure 2) . 

To understand the involved pathways and functions of the DECGs and DECMs in LUAD, enrichment analysis was performed (Supplementary Figure 3) . A total of 41 pathways and 129 GO functions were enriched based on DECGs and DECMs. Interestingly, both 2 DECGs (CXCL9, IL6) and 4 DECMs (miR-125b-5p, miR-9-5p, miR-130b-5p, and miR-421) were enriched in the toll-like receptor signaling pathway (Supplementary Figure 4) . In this pathway, CXCL9 takes part in chemotactic effects of T cells while IL6 participates in proinflammatory effects. Then according to the relationships between DECGs, DECMs and virusassociated processes, an interaction network was constructed ( Figure 4) . Notably, not only ACE2 but IL6 were predicted as the targets of miR-125b-5p. Similarly, both ACE2 and CXCL9 were regarded as the predicted target of miR-9-5p. However, due to the expression trend of CXCL9 and miR-9-5p were not opposite with ACE2, it was suggested that both CXCL9 and miR-9-5p may not play a decisive role in ACE2 expression.

In conclusion, the reduced miR-125b-5p might be the primary inhibitor of ACE2 in LUAD. Once ACE2 was dysregulated, IL6 in toll-like receptor pathway might activate the immune system as a downstream effector. No matter the defense of SARS-CoV-2 or pneumonia in LUAD, the altered immunoreaction was the primary cause. Thus, the ability of ACE2 and IL6 to regulate immune system was evaluated based on the correlation between gene characters (mRNA expression, copy number and methylation) and immune elements (lymphocyte, immune-inhibitor, immune-stimulator, MHC molecule, chemokine, and chemokine receptor) ( Figure 5 ). Interestingly, CXCL9 could mediate AGING immunoreaction as a chemokine which also influenced by ACE2 and IL6 (expression level and methylation). Therefore, it was suggested that even CXCL9 might not participate in regulating ACE2 level, it would like to mediate immune-regulation of ACE2 and IL6.

Once infected with SARS-CoV-2, the pathogenic T cells were activated rapidly followed by immune response and virus eliminating. However, the fatal cause of viral pneumonia was the uncontrolled inflammation. In addition, previous study has demonstrated that COVID-19 might be more likely to cause an endothelial damage in the vessel system. ACE2, not only as the functional receptor of coronavirus, also emerged as a potent negative regulator of the RAS which maintaining homeostasis of blood pressure and inflammatory responses [27] . Altered activation of the RAS is often attributed to the pathogenesis of many diseases such as hypertension, myocardial infarction and inflammatory lung disease [28, 29] . Compared with normal COVID-19 patients, whether the vessel system of LUAD has higher ACE2 level is unclear. In murine acute respiratory distress syndrome models, it was confirmed that the lack of ACE2 expression in the lung resulted in attenuated vascular permeability, enhanced lung edema, neutrophil infiltration, and further deteriorated lung function [30] . Thus, though higher ACE2 level would increase the risk of SARS-CoV-2 infection, the normal level of ACE2 would also improve prognosis of patients with acute lung injury.

The RAS contained ACE2 participate in maintaining homeostasis of inflammatory responses. However, once SARS-CoV-2 binds to ACE2, the abundance on the cell surface, mRNA expression and the enzymatic activity of ACE2 are significantly reduced which lead to the inflammatory storm. Coincidentally, IL6 was the main factor which could also induce an inflammatory storm. Once infected with SARS-CoV-2, the activated pathogenic T cells could produce granulocytemacrophage colony-stimulating factors (GM-CSF) and IL6. The GM-CSF could further activate CD14+CD16+ inflammatory monocytes. In addition, more cytokines include IL6 were produced via positive feedback loop. A large concentration of immune cells and tissue fluid in the lungs can block gas exchange between alveoli and capillaries, leading to acute respiratory distress syndrome. Once a cytokine storm forms, the immune system would kill lots of normal cells in the lung while eliminating the virus, severely damaging the lung's ventilation function. In a word, similar to ACE2, IL6 play key role in defense of SARS-CoV-2 with normal expression level. Once deregulated, IL6 will become a lethal factor for patients.

Conventionally, the LUAD patients have lower immunity than healthy persons. The increased ACE2 level in the lung might lead to susceptibility to infection. Nevertheless, the pathology of COVID-19 is complicated, except miRNA introduced in present study, transcription factor, lncRNA and ceRNA may participate in mechanism of COVID. The result of our study may provide a reference for further research.

Dataset GSE74190 was acquired from GEO Dataset. All experiments were approved by the local ethics committee. Expression profiles of 723 human miRNAs were investigated in selected cancerous cells and normal cells populations derived from 36 LUAD and 44 adjacent normal tissues. ALL cells derived from 82 snap-frozen surgical specimens. The sequencing platform is Agilent-019118 Human miRNA Microarray V1 G4470A [miRBase release 9.1 miRNA ID version].

The potential target genes of differentially expressed miRNAs in LUAD were acquired from the widely used online databases TargetScanHuman 7.2 [31] and DIANA-microT [32] . To reduce false positives, the predicted target genes which appeared at both databases were accepted. Following this, the list of predicted target genes of individual miRNAs was imported to DIANA-mirPath, a miRNA pathway analysis web server. Canonical pathways significantly affected by individual differentially expressed miRNAs were made a contrast and then we obtained the common pathways of which elements were also as the targets of deregulation miRNAs. All these canonical pathways were identified from the Kyoto Encyclopaedia of Genes and Genomes (KEGG) databases.

Total RNA from tissues was isolated using Trizol extractions (Invitrogen). The RNA quantity was assessed by NanoDrop®ND-1000 spectrophotometer and RNA 6000 NanoChips with the Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA, USA). 100 ng of total RNA was amplified using the Ambion® WT Expression Kit (4411973, Life Technologies). Small RNA libraries were prepared using 1μg of total RNA according to the TruSeq Small RNA Sample Preparation Guide (Illumina, San Diego, CA, USA). To generate count data, the raw sequences were compared to human mature miRNA sequences (from miRBase version 17) and non-coding RNA sequences (Rfam version 10) by MEGABLAST.

Background deletion, quantile normalization, and probe assembly were performed. Different expression mRNAs between normal vs. LUAD samples were detected by the R package DESeq [33] . P-values were adjusted for multiple comparisons via using the Benjamini-Hochberg procedure [34] . MiRNA with an adjusted pvalue of < 0.05 and logFC ≥ 2.0 was considered as differentially expressed. Gene and miRNA enrichment analyses were performed with DAVID version 6.7 and DIANA-mirPATH v.3 respectively. The enriched biological GO and pathway terms were identified [35] . The interaction network was drawn by Cytoscape. Some other databases used are listed in Table 1 .

Results are presented as mean values ± standard error of the mean (SEM). Unless mentioned otherwise, the statistical comparison between groups was performed by using t-test, a maximum of three comparisons were performed per panel, and robustness of statistical significance was verified after correction for multiple testing. Probability was considered to be significant at p < 0.05.

Publicly available datasets were analyzed in this study. All data was available.

coronavirus pneumonia-19; ACE2: angiotensin-converting enzyme 2; Ang II: angiotensin II; Ang 1-7: angiotensin 1-7; RAS: renin-angiotensin system; LUAD: lung adenocarcinoma; NSCLC: non-small cell lung cancer; KEGG: Kyoto Encyclopaedia of Genes and Genomes; SEM: standard error of the mean; MasR: Mas receptor; DECG: differentially expressed correlative gene; DECM: differentially expressed correlative miRNA; RAS: renin-angiotensin system; GM-CSF: granulocytemacrophage colony-stimulating factors.

All authors contributed equally to this work. LC, YL, HL and LZ contributed to the conception and design of the work. LC analyzed data and wrote this manuscript. CD, JT and LZ contributed to interpretation of data for the work. All authors contributed to final approval of the version to be published. All authors agree to be accountable for all aspects of the work ACKNOWLEDGMENTS This study was supported by the sequencing and clinical datasets from TCGA and GEO databases. Thanks for guidance from Pro. Zhong and Dr. Liu.

The authors declare that they have no conflicts of interest. 

The SARS-CoV S glycoprotein: expression and functional characterization

Retroviruses pseudotyped with the severe acute respiratory syndrome coronavirus spike protein efficiently infect cells expressing angiotensinconverting enzyme 2

Coronavirus spike proteins in viral entry and pathogenesis

A 193-amino acid fragment of the SARS coronavirus S protein efficiently binds angiotensin-converting enzyme 2

Amino acids 270 to 510 of the severe acute respiratory syndrome coronavirus spike protein are required for interaction with receptor

Angiotensin-converting enzyme 2 and angiotensin 1-7: novel therapeutic targets

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury

A human homolog of angiotensinconverting enzyme. Cloning and functional expression as a captopril-insensitive carboxypeptidase

Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2

Ectodomain shedding of angiotensin converting enzyme 2 in human airway epithelia

Loss of angiotensin-converting enzyme 2 enhances TGFβ/smad-mediated renal fibrosis and NF-κB-driven renal inflammation in a mouse model of obstructive nephropathy

Angiotensin-converting enzyme 2/angiotensin-(1-7)/Mas axis protects against lung fibrosis by inhibiting the MAPK/NF-κB pathway

Angiotensinconverting enzyme 2 protects from lethal avian influenza a H5N1 infections

Angiotensin-converting enzyme 2 is a potential therapeutic target for EGFRmutant lung adenocarcinoma

Angiotensin-converting enzyme 2 attenuates the metastasis of non-small cell lung cancer through inhibition of epithelial-mesenchymal transition

Pulmonary angiotensin-converting enzyme 2 (ACE2) and inflammatory lung disease

Physiological and pathological regulation of ACE2, the SARS-CoV-2 receptor

and China Medical Treatment Expert Group for Covid-19. Clinical characteristics of coronavirus disease 2019 in China

Adenosquamous carcinoma of the lung

ACE2: a new target for cardiovascular disease therapeutics

Characterization and significance of ACE2 and Mas receptor in human colon adenocarcinoma

Inhibition of human lung cancer cell growth by angiotensin

Angiotensin-(1-7) inhibits growth of human lung adenocarcinoma xenografts in nude mice through a reduction in cyclooxygenase-2

A century old reninangiotensin system still grows with endless possibilities: AT1 receptor signaling cascades in cardiovascular physiopathology

Serum and pulmonary angiotensin converting enzyme as a marker of acute lung

Role of the renin-angiotensin system in the systemic microvascular inflammation of alveolar hypoxia

Angiotensin-converting enzyme 2 protects from severe acute lung failure

Prediction of mammalian microRNA targets

A combined computational-experimental approach predicts human microRNA targets

Differential expression analysis for sequence count data

Mathematical formulae for the prediction of the residual beta cell function during the first two years of disease in children and adolescents with insulin-dependent diabetes mellitus

Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources

Supplementary Figure 1. The correlated mRNAs (A) and miRNAs (B) of ACE2 in LUAD. Due to the number of correlated genes is too large, only top 50 mRNAs were list. The red point represents positive correlation while the blue ones represents negative correlation. (C) The sketch map of predicted target sites of miRNAs in ACE2 3'-UTR

Integrated plot of clinical data and 7 mutation of DECGs in 586 LUAD samples. From top to bottom panels indicate: mutation spectrum, overall survival status, overall survival (months), tissue source site, American Joint Committee on Cancer tumor stage code, patients' vital status, sex, tissue prospective collection indicator, tissue retrospective collection indicator, mutation symbol of DECGs, heatmap of DECGs expression and methylation