ACE2 and TMPRSS2 variants and expression as candidates to sex and country differences in COVID-19 severity in Italy

Rosanna Asselta, Elvezia Maria Paraboschi, Alberto Mantovani, Stefano Duga, Rosanna Asselta, Elvezia Maria Paraboschi, Alberto Mantovani, Stefano Duga

Abstract

As the outbreak of coronavirus disease 2019 (COVID-19) progresses, prognostic markers for early identification of high-risk individuals are an urgent medical need. Italy has one of the highest numbers of SARS-CoV-2-related deaths and one of the highest mortality rates. Worldwide, a more severe course of COVID-19 is associated with older age, comorbidities, and male sex. Hence, we searched for possible genetic components of COVID-19 severity among Italians by looking at expression levels and variants in ACE2 and TMPRSS2 genes, crucial for viral infection.Exome and SNP-array data from a large Italian cohort were used to compare the rare-variants burden and polymorphisms frequency with Europeans and East Asians. Moreover, we looked into gene expression databases to check for sex-unbalanced expression.While we found no significant evidence that ACE2 is associated with disease severity/sex bias, TMPRSS2 levels and genetic variants proved to be possible candidate disease modulators, prompting for rapid experimental validations on large patient cohorts.

Keywords: ACE2; COVID-19; SARS-CoV-2; TMPRSS2; genetic variants.

Conflict of interest statement

CONFLICTS OF INTEREST: No conflicts of interest to disclose

Figures

Figure 1
Figure 1
ACE2 expression levels. All panels show ACE2 mRNA expression levels in human normal lung samples stratified according to sex (or on sex and age). On left panels, data were retrieved for a total of 578 RNAseq experiments from the GTex repository. Expression levels are reported as transcripts per kilobase million (TPM). On the right, data were collected from two different datasets (GSE66499 and GSE19804) from the GEO database. Expression levels are reported as normalized signal intensities. P values were calculated by using either the Kruskal-Wallis or the student t test, using the R software (https://www.r-project.org/).
Figure 2
Figure 2
TMPRSS2 expression levels and eQTLs. (A) Both panels show TMPRSS2 mRNA expression levels in human normal lung samples stratified according to sex. On the left, data were retrieved for a total of 578 RNAseq experiments from the GTex repository. Expression levels are reported as transcripts per kilobase million (TPM). On the right, data were collected for a total of 170 microarray experiments from the GEO database. Expression levels are reported as normalized signal intensities. P values were calculated by using either the Kruskal-Wallis or the student t test. (B) Screenshot from the UCSC Genome browser (http://genome.ucsc.edu/; GRCh37/hg19) highlighting the TMPRSS2 region (coordinates chr21: 42,835,000-42,905,000). The panel shows the following tracks: i) the ruler with the scale at the genomic level; ii) chromosome 21 nucleotide numbering; iii) the UCSC RefSeq track; iv) enhancers (grey and red bars) from GeneHancer database; v) interactions (curved lines) connecting GeneHancer regulatory elements and genes: all curved lines converge towards the androgen-responsive enhancer for the TMPRSS2 gene described by Clinckemalie and colleagues [29].

References

    1. Zhang JJ, Dong X, Cao YY, Yuan YD, Yang YB, Yan YQ, Akdis CA, Gao YD. Clinical characteristics of 140 patients infected with SARS-CoV-2 in Wuhan, China. Allergy. 2020. [Epub ahead of print]. 10.1111/all.14238
    1. Andersen KG, Rambaut A, Lipkin WI, Holmes EC, Garry RF. The proximal origin of SARS-CoV-2. Nat Med. 2020; 26:450–52. 10.1038/s41591-020-0820-9
    1. Oksanen A, Kaakinen M, Latikka R, Savolainen I, Savela N, Koivula A. Regulation and trust: 3-month follow-up study on COVID-19 mortality in 25 european countries. JMIR Public Health Surveill. 2020; 6:e19218. 10.2196/19218
    1. Jin JM, Bai P, He W, Wu F, Liu XF, Han DM, Liu S, Yang JK. Gender differences in patients with COVID-19: focus on severity and mortality. Front Public Health. 2020; 8:152. 10.3389/fpubh.2020.00152
    1. Channappanavar R, Fett C, Mack M, Ten Eyck PP, Meyerholz DK, Perlman S. Sex-based differences in susceptibility to severe acute respiratory syndrome coronavirus infection. J Immunol. 2017; 198:4046–53. 10.4049/jimmunol.1601896
    1. Alghamdi IG, Hussain II, Almalki SS, Alghamdi MS, Alghamdi MM, El-Sheemy MA. The pattern of middle east respiratory syndrome coronavirus in Saudi Arabia: a descriptive epidemiological analysis of data from the saudi ministry of health. Int J Gen Med. 2014; 7:417–23. 10.2147/IJGM.S67061
    1. Klein SL, Flanagan KL. Sex differences in immune responses. Nat Rev Immunol. 2016; 16:626–38. 10.1038/nri.2016.90
    1. Zhou P, Yang XL, Wang XG, Hu B, Zhang L, Zhang W, Si HR, Zhu Y, Li B, Huang CL, Chen HD, Chen J, Luo Y, et al.. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020; 579:270–73. 10.1038/s41586-020-2012-7
    1. Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, Schiergens TS, Herrler G, Wu NH, Nitsche A, Müller MA, Drosten C, Pöhlmann S. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell. 2020; 181:271–80.e8. 10.1016/j.cell.2020.02.052
    1. Kuba K, Imai Y, Rao S, Gao H, Guo F, Guan B, Huan Y, Yang P, Zhang Y, Deng W, Bao L, Zhang B, Liu G, et al.. A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury. Nat Med. 2005; 11:875–9. 10.1038/nm1267
    1. Cheng Z, Zhou J, To KK, Chu H, Li C, Wang D, Yang D, Zheng S, Hao K, Bossé Y, Obeidat M, Brandsma CA, Song YQ, et al.. Identification of TMPRSS2 as a susceptibility gene for severe 2009 pandemic a(H1N1) influenza and a(H7N9) influenza. J Infect Dis. 2015; 212:1214–21. 10.1093/infdis/jiv246
    1. Baena E, Shao Z, Linn DE, Glass K, Hamblen MJ, Fujiwara Y, Kim J, Nguyen M, Zhang X, Godinho FJ, Bronson RT, Mucci LA, Loda M, et al. ETV1 directs androgen metabolism and confers aggressive prostate cancer in targeted mice and patients. Genes Dev. 2013; 27:683–98. 10.1101/gad.211011.112
    1. Tukiainen T, Villani AC, Yen A, Rivas MA, Marshall JL, Satija R, Aguirre M, Gauthier L, Fleharty M, Kirby A, Cummings BB, Castel SE, Karczewski KJ, et al., and GTEx Consortium, Laboratory, Data Analysis & Coordinating Center (LDACC)—Analysis Working Group, Statistical Methods groups—Analysis Working Group, Enhancing GTEx (eGTEx) groups, NIH Common Fund, NIH/NCI, NIH/NHGRI, NIH/NIMH, NIH/NIDA, Biospecimen Collection Source Site—NDRI, Biospecimen Collection Source Site—RPCI, Biospecimen Core Resource—VARI, Brain Bank Repository—University of Miami Brain Endowment Bank, Leidos Biomedical—Project Management, ELSI Study, Genome Browser Data Integration & Visualization—EBI, and Genome Browser Data Integration & Visualization—UCSC Genomics Institute, University of California Santa Cruz. Landscape of X chromosome inactivation across human tissues. Nature. 2017; 550:244–48. 10.1038/nature24265
    1. Cai G. Bulk and single-cell transcriptomics identify tobacco-use disparity in lung gene expression of ACE2, the receptor of 2019-nCov. Preprints. 2020;2020020051. 10.20944/preprints202002.0051.v2
    1. Zhao Y, Zhao Z, Wang Y, Zhou Y, Ma Y, Zuo W. Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCov. bioRxiv. 2020. 10.1101/2020.01.26.919985
    1. Lyon MF. Gene action in the X-chromosome of the mouse (Mus musculus L.). Nature. 1961; 190:372–3. 10.1038/190372a0
    1. Carrel L, Willard HF. X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature. 2005; 434:400–4. 10.1038/nature03479
    1. Cao Y, Li L, Feng Z, Wan S, Huang P, Sun X, Wen F, Huang X, Ning G, Wang W. Comparative genetic analysis of the novel coronavirus (2019-nCoV/SARS-CoV-2) receptor ACE2 in different populations. Cell Discov. 2020; 6:11. 10.1038/s41421-020-0147-1
    1. Do R, Stitziel NO, Won HH, Jørgensen AB, Duga S, Angelica Merlini P, Kiezun A, Farrall M, Goel A, Zuk O, Guella I, Asselta R, Lange LA, et al., and NHLBI Exome Sequencing Project. Exome sequencing identifies rare LDLR and APOA5 alleles conferring risk for myocardial infarction. Nature. 2015; 518:102–6. 10.1038/nature13917
    1. Paraboschi EM, Khera AV, Merlini PA, Gigante L, Peyvandi F, Chaffin M, Menegatti M, Busti F, Girelli D, Martinelli N, Olivieri O, Kathiresan S, Ardissino D, et al.. Rare variants lowering the levels of coagulation factor X are protective against ischemic heart disease. Haematologica. 2019. [Epub ahead of print]. 10.3324/haematol.2019.237750
    1. Chaoxin J, Daili S, Yanxin H, Ruwei G, Chenlong W, Yaobin T. The influence of angiotensin-converting enzyme 2 gene polymorphisms on type 2 diabetes mellitus and coronary heart disease. Eur Rev Med Pharmacol Sci. 2013; 17:2654–9.
    1. Yang M, Zhao J, Xing L, Shi L. The association between angiotensin-converting enzyme 2 polymorphisms and essential hypertension risk: a meta-analysis involving 14,122 patients. J Renin Angiotensin Aldosterone Syst. 2015; 16:1240–4. 10.1177/1470320314549221
    1. Wu YH, Li JY, Wang C, Zhang LM, Qiao H. The ACE2 G8790A polymorphism: involvement in type 2 diabetes mellitus combined with cerebral stroke. J Clin Lab Anal. 2017; 31:e22033. 10.1002/jcla.22033
    1. Kron KJ, Murison A, Zhou S, Huang V, Yamaguchi TN, Shiah YJ, Fraser M, van der Kwast T, Boutros PC, Bristow RG, Lupien M. TMPRSS2-ERG fusion co-opts master transcription factors and activates NOTCH signaling in primary prostate cancer. Nat Genet. 2017; 49:1336–45. 10.1038/ng.3930
    1. Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, Collins RL, Laricchia KM, Ganna A, Birnbaum DP, Gauthier LD, Brand H, Solomonson M, et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. bioRxiv. 2019.
    1. FitzGerald LM, Agalliu I, Johnson K, Miller MA, Kwon EM, Hurtado-Coll A, Fazli L, Rajput AB, Gleave ME, Cox ME, Ostrander EA, Stanford JL, Huntsman DG. Association of TMPRSS2-ERG gene fusion with clinical characteristics and outcomes: results from a population-based study of prostate cancer. BMC Cancer. 2008; 8:230. 10.1186/1471-2407-8-230
    1. Giri VN, Ruth K, Hughes L, Uzzo RG, Chen DY, Boorjian SA, Viterbo R, Rebbeck TR. Racial differences in prediction of time to prostate cancer diagnosis in a prospective screening cohort of high-risk men: effect of TMPRSS2 Met160Val. BJU Int. 2011; 107:466–70. 10.1111/j.1464-410X.2010.09522.x
    1. Kathiresan S, Voight BF, Purcell S, Musunuru K, Ardissino D, Mannucci PM, Anand S, Engert JC, Samani NJ, Schunkert H, Erdmann J, Reilly MP, Rader DJ, et al., and Myocardial Infarction Genetics Consortium, and Wellcome Trust Case Control Consortium. Genome-wide association of early-onset myocardial infarction with single nucleotide polymorphisms and copy number variants. Nat Genet. 2009; 41:334–41. 10.1038/ng.327
    1. Clinckemalie L, Spans L, Dubois V, Laurent M, Helsen C, Joniau S, Claessens F. Androgen regulation of the TMPRSS2 gene and the effect of a SNP in an androgen response element. Mol Endocrinol. 2013; 27:2028–40. 10.1210/me.2013-1098
    1. Jernigan DB, Cox NJ. H7N9: preparing for the unexpected in influenza. Annu Rev Med. 2015; 66:361–71. 10.1146/annurev-med-010714-112311
    1. Ghosh S, Klein RS. Sex drives dimorphic immune responses to viral infections. J Immunol. 2017; 198:1782–90. 10.4049/jimmunol.1601166

Source: PubMed

Sponsorer og samarbeidspartnere

Medisinsk tilstand

Legemiddelintervensjoner

3
Abonnere