Kristina Schönfelder1Katharina Breuckmann2Carina Elsner3Ulf Dittmer3David Fistera4Frank Herbstreit5
Joachim Risse4Karsten Schmidt5
Sivagurunathan Sutharsan6Christian Taube6Karl-Heinz Jöckel7Winfried Siffert8Andreas Kribben1
Birte Möhlendick8*- 1Department of Nephrology, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- 2Institute of Diagnostic and Interventional Radiology and Neuroradiology, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- 3Institute for Virology, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- 4Center of Emergency Medicine, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- 5Department of Anesthesiology and Intensive Care Medicine, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- 6Department of Pulmonary Medicine, Ruhrlandklinik, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- 7Institute of Medical Informatics, Biometry and Epidemiology, University of Duisburg-Essen, Essen, Germany
- 8Institute of Pharmacogenetics, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
The transmembrane serine protease 2 (TMPRSS2) is the major host protease that enables entry of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) into host cells by spike (S) protein priming. Single nucleotide polymorphisms (SNPs) in the gene TMPRSS2 have been associated with susceptibility to and severity of H1N1 or H1N9 influenza A virus infections. Functional variants may influence SARS-CoV-2 infection risk and severity of Coronavirus disease 2019 (COVID-19) as well. Therefore, we analyzed the role of SNPs in the gene TMPRSS2 in a German case-control study. We performed genotyping of the SNPs rs2070788, rs383510, and rs12329760 in the gene TMPRSS2 in 239 SARS-CoV-2-positive and 253 SARS-CoV-2-negative patients. We analyzed the association of the SNPs with susceptibility to SARS-CoV-2 infection and severity of COVID-19. SARS-CoV-2-positive and SARS-CoV-2-negative patients did not differ regarding their demographics. The CC genotype of TMPRSS2 rs383510 was associated with a 1.73-fold increased SARS-CoV-2 infection risk, but was not correlated to severity of COVID-19. Neither TMPRSS2 rs2070788 nor rs12329760 polymorphisms were related to SARS-CoV-2 infection risk or severity of COVID-19. In a multivariable analysis (MVA), the rs383510 CC genotype remained an independent predictor for a 2-fold increased SARS-CoV-2 infection risk. In summary, our report appears to be the first showing that the intron variant rs383510 in the gene TMPRSS2 is associated with an increased risk to SARS-CoV-2 infection in a German cohort.
Introduction
The transmembrane serine protease 2 (TMPRSS2) is the major host protease that enables cell entry of several coronaviruses, e.g., HCoV-229E, SARS-CoV, and MERS, respectively (Shulla et al., 2011; Bertram et al., 2013; Shirato et al., 2013). TMPRSS2 is strongly expressed in lung tissue and bronchial transient secretory cells (Lukassen et al., 2020). In a recent study, Hoffmann et al. (2020) demonstrated that SARS-CoV-2 uses the angiotensin-converting enzyme 2 (ACE2) as a receptor for host cell binding, and the viral spike (S) protein is cleaved into S1 and S2 by TMPRSS2 to allow fusion of the viral and cellular membranes.
Single nucleotide polymorphisms (SNPs) in the gene TMPRSS2 (chromosome 21q22.3, reverse strand) have been associated with susceptibility to and severity of H1N1 or H1N9 influenza A virus infections (Cheng et al., 2015). Currently, SNPs altering gene expression or protein function are frequently discussed in the light of Coronavirus disease 2019 (COVID-19) in literature reviews and in in silico analyses utilizing genotyping data from existing databases (Asselta et al., 2020; Barash et al., 2020; Fujikura and Uesaka, 2020; Ghafouri-Fard et al., 2020; Hou et al., 2020; Paniri et al., 2020; Russo et al., 2020; Senapati et al., 2020; Singh et al., 2020; Smatti et al., 2020; Vargas-Alarcón et al., 2020; Zarubin et al., 2020). Three variants, rs2070788, rs383510, and rs12329760, seem to be of particular interest based on preliminary work in the context of other virus infections.
In a study with Chinese patients, a genetic predisposition to severe H1N1 influenza A virus infection was observed for TMPRSS2 rs2070788 GG genotype (OR: 2.11, 95% CI 1.18–3.77, p = 0.01) or G-allele (OR: 1.54, 95% CI 1.14–2.06, p = 0.004) carriers (Cheng et al., 2015). The SNP rs2070788 is a G to A transition located in intron 11–12 (c.1171+587C>T) of TMPRSS2. Minor allele frequencies (MAFs) for the G-allele are higher in Europeans (0.46) than East Asians (0.36) or Africans (0.27; Clarke et al., 2017). In a lung expression quantitative trait loci (eQTL) analysis, rs2070788 was correlated to differential expression of TMPRSS2, with GG genotype carriers showing the highest expression levels (GTEx portal; Irham et al., 2020). The functional basis is still not understood, since rs2070788 is not located in any known regulatory region (ensembl.org, release 102, November 2020).
A second intronic SNP, rs383510 (c.445+1954A>G, intron 5–6), which is in high linkage disequilibrium (LD) with rs2070788 in East Asians (r2/D' > 0.8) but not Europeans (r2 = 0.4, D' = 0.7), is located in a putative enhancer region (Cheng et al., 2015). The MAF of the T-allele is about one and a half times higher for Europeans (0.49) than for East Asians (0.36) or Africans (0.33; Clarke et al., 2017). This variant was shown to be an eQTL in the lung with TT genotype carriers showing the highest expression levels (GTEx portal; Irham et al., 2020). Cheng et al. (2015) observed a significant association with H7N9 influenca A virus infection risk for rs2070788 GG or GA genotype carriers (OR: 1.99, 95% CI 1.14–3.48, p = 0.02) and rs383510 TT or TC genotype carriers (OR: 2.01, 95% CI 1.15–3.54), p = 0.01) in their study comparing the genotype frequencies of Chinese H7N9 influenca A virus infected patients (N = 102) to healthy controls (N = 106; Cheng et al., 2015).
Because of its functional impact, the missense variant rs12329760 (c.478G>A, p.V160M), which alters the amino acid valine at position 160 to methionine, might as well be of interest with regard to COVID-19. The SNP is located in the scavenger receptor cysteine-rich (SRCR) domain of TMPRSS2, which interacts with external pathogens (Paoloni-Giacobino et al., 1997). The amino acid exchange may create a de novo pocket protein (Paniri et al., 2020). It is an eQTL in whole blood with TT genotype carriers showing the highest expression levels (GTEx portal). The MAF of the T-allele in East Asians (0.36) is comparably higher than in Europeans (0.24), South Asians (0.23), or Africans (0.29; Clarke et al., 2017). In prostate cancer, the rs12329760 T-allele is associated with multiple copies of the TMPRSS2-ERG gene fusion responsible for shorter survival and higher tumor recurrence (FitzGerald et al., 2008). In a preliminary small scale, whole-exome sequencing study comprising 131 Italian hospitalized COVID-19 patients, Latini et al. (2020) observed a lower rs12329760 T-allele frequency (p = 0.02) in the patients compared to corresponding European GnomAD database controls.
Based on the above described associations, we investigated the influence of the variants rs2070788, rs383510, and rs12329760 in the gene TMPRSS2 on SARS-CoV-2 infection risk and severity of COVID-19 in a German cohort.
Materials and Methods
Study Participants, Recruitment, and Outcome of the Patients
The study was conducted following the approval of the Ethics Committee of the Medical Faculty of the University of Duisburg-Essen (20-9230-BO) and in cooperation with the West German Biobank (WBE; 20-WBE-088). Written informed consent was obtained from the study patients.
Enrolment started on March 11, 2020, and ended on September 30, 2020. Patients were initially recruited upon presentation with COVID-19 typical symptoms, i.e., fever, cough, and dyspnea or who were admitted to the hospital with already confirmed SARS-CoV-2 infection. We included 239 SARS-CoV-2-positive patients. Patients were classified as SARS-CoV-2-positive with at least one positive real-time reverse transcription PCR (RT-PCR) test result. Follow-up was completed on October 31, 2020, at which time all patients either were discharged from the hospital as “cured” or had a fatal outcome of the disease. We also studied 253 SARS-CoV-2-negative patients, who presented with COVID-19 typical symptoms, but were tested exclusively negative for SARS-CoV-2 by RT-PCR. These patients were hospitalized at the University Hospital Essen or treated as outpatients, due to other medical conditions. Clinical outcome was defined as follows according to the criteria of the European Center of Disease Prevention and Control (ECDC; 2020) - “serious”: hospitalized patients admitted to an intensive care unit and all cases of COVID-19‑related deaths during the hospital stay; “moderate”: outpatients and all other hospitalized patients. In contrast to the ECDC classification, where patients are counted up to three times, every patient only counted once according to the worst clinical outcome observed during the hospital stay in our study.
Genotyping of TMPRSS2
PCR was performed with 2 μl genomic DNA and 30 μl Taq DNA-Polymerase 2x Master Mix Red (Ampliqon, Odense, Denmark) with the following conditions: initial denaturation 95°C for 5 min; 38 cycles with denaturation 95°C for 30 s; annealing at 60°C (rs12329760) or 65°C (rs2070788 and rs383510, respectively) for 30 s and elongation 72°C for 30 s each; and final elongation 72°C for 10 min. The following oligonucleotide primers were used for the respective PCR reactions: rs2070788 5' [BIO]AGTTTCTGCTGATGAGGAGCC 3', 3'GAAGTGCTTAGTGGCAGGCA 5'; rs383510 5' [BIO]ATGGCTGTGCTTGGGAAATAAC 3', 3' CTTATTTCCTGGCCGGACGC 5'; and rs12329760 5' CGCCCGTAGTTCTCGTTCC 3', 3' TTCGCCTCTACGGACCAAAC 5'. PCR products for rs12329760 were digested with Hpy8I (Thermo Scientific, Dreireich, Germany), and restriction fragments were analyzed by agarose gel electrophoresis. For the various genotypes results from restriction fragment length polymorphism (RFLP)-PCR were validated by Sanger sequencing (Supplementary Figure S1). Genotyping of rs2070788 and rs383510 was performed by Pyrosequencing according to the manufacturers’ instructions (Qiagen, Hilden, Germany). In brief, biotinylated PCR amplicons were immobilized on streptavidin-coated sepharose beads (GE Healthcare, Solingen, Germany) with the vacuum tool of the workstation (Qiagen, Hilden, Germany). In the next step, PyroMark denaturation solution (Qiagen, Hilden, Germany) was used to separate the complementary strand from the biotinylated strand. The biotinylated, single-stranded DNA remains immobilized on the vacuum tool. Finally, the single-stranded DNA was released into the sequencing plate containing 0.3 μM sequencing primer (5' TTTATATCCTTCTCAAAAG 3' for rs2070788 and 5' ACGTGGTGGTGCGGGCCT 3' for rs383510, respectively). After incubation at 80°C for 2 min, the hybridized primer and single-stranded template were incubated with the enzymes DNA polymerase, ATP sulfurylase, luciferase, and apyrase, as well as the substrates adenosine 5' phosphosulfate (APS) and luciferin. Addition of dideoxyribonucleotide triphosphates (ddNTPs) was performed sequentially. Each incorporation event is accompanied by the release of pyrophosphate (PPi) in a quantity equimolar to the amount of the incorporated nucleotide. ATP sulfurylase converts PPi to ATP, which drives luciferase-mediated conversion of luciferin to oxyluciferin that generates visible light. The height of each light signal is proportional to the number of nucleotides incorporated and can be visualized in the Pyrogram. Apyrase continuously degrades unincorporated nucleotides and ATP and when degradation is complete another nucleotide is added. Analyses were performed on the PyroMark Q96 MD instrument (Qiagen, Hilden, Germany).
Transmembrane serine protease 2 is located on the reverse strand, and all polymorphisms are reported on the forward strand in this study.
Statistical Analyses
Hardy-Weinberg equilibrium (HWE) was calculated using Pearson’s Χ2 goodness of fit test, and samples were considered as deviant from HWE at a significance level of p < 0.05. LD analysis was performed using HaploView V4.2.
For genetic association, we calculated odds ratio (OR) and 95% CI by Fisher’s exact test using Baptista-Pike method for OR, respectively. Values of p are reported two-sided and values of p < 0.05 were considered significant. Adjustment for multiple comparisons in the genotype-phenotype association analysis was performed using Holm-Bonferroni method (α = 0.0167). MVA was performed for both SARS-COV-2 infection risk and COVID-19 severity to estimate independency of the variables age, sex, and TMPRSS2, and rs2070788, rs383510, and rs12329760 genotypes by logistic regression (likelihood ratio test, backwards).
The eQTL data described in this manuscript were obtained from the Genotype-Tissue Expression (GTEx) Portal (gtexportal.org; GTEx analysis release V8; dbGaP Accession phs000424.v8.p2) on 01/31/2021.
Results
From March 11, 2020 to October 31, 2020, we enrolled and studied 239 SARS-CoV-2-positive and 253 SARS-CoV-2-negative patients to determine associations of the SNPs rs2070788, rs383510, and rs12329760 in the gene TMPRSS2 with susceptibility to SARS-CoV-2 infection and severity of COVID-19. The characteristics and genotypes of SARS-CoV-2-positive and -negative patients are summarized in Table 1. Distribution of sex (p = 0.86) and age (p = 0.10) was similar in both groups.
Table 1. Demographics, transmembrane serine protease 2 (TMPRSS2) genotypes, and outcome of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2)-positive and SARS-CoV-2-negative patients.
Transmembrane serine protease 2 rs2070788 genotypes were compatible with HWE (p = 0.87 and p = 0.73, respectively). The MAF of the G-allele in SARS-CoV-2-positive patients (0.48) was similar to the MAF in SARS-CoV-2-negative patients (0.45). For TMPRSS2 rs12329760, we also did not observe a deviation from HWE for both SARS-CoV-2-negative and -positive patients (p = 0.66 and p = 0.50, respectively). The MAF of the T-allele was higher in SARS-CoV-2-positive patients (0.24 vs. 0.20). Neither the rs2070788 nor rs12329760 variant was correlated with an increased SARS-CoV-2 infection risk or severity of COVID-19 (Table 1).
The observed genotype frequencies for TMPRSS2 rs383510 also were consistent with HWE in SARS-CoV-2-positive patients (p = 0.80). There was a slight deviation from the HWE in SARS-CoV-2-negative patients (p = 0.04). We found 20.1% TT, 48.5% TC, and 31.4% CC genotype carriers in the infected patients. The following genotype frequencies were observed in the non-infected group: 22.5% TT, 56.5% TC, and 21.0% CC. The MAF of the T-allele was higher in the SARS-CoV-2-negative patients (0.51) compared to the MAF in the SARS-CoV-2-positive patients (0.44). Carriers of the CC genotype had a 1.73-fold increased risk to SARS-CoV-2 infection (OR: 1.73, 95% CI 1.15–2.59, p = 0.01). We did not observe an association of any allele to severity of COVID-19.
In this study, we did not observe linkage disequilibrium between the alleles of the three loci. The polymorphisms rs2070788 and rs383510 showed highest LD values (r2 = 0.49, D' = 0.72).
In a multivariable analysis, TMPRSS2 rs383510 remained an independent predictor for SARS-CoV-2 infection risk (OR: 2.00, 95% CI 1.30–3.08, p = 0.002). In a second multivariable analysis, male sex was the only independent predictor for severity of COVID-19 [OR: 2.64, 95% CI (1.44–4.84), p = 0.002].
Discussion and Study Limitations
Currently, the influence of host genetic factors on SARS-CoV-2 infection risk or COVID-19 severity is being investigated in a variety of in silico analyses and with targeted or genome-wide sequencing approaches.
Recently, both blood group A and inborn errors of interferon I immunity were identified as risk factors for severity of COVID-19 in genome-wide association studies (Ellinghaus et al., 2020; Zhang et al., 2020).
So far, mainly in silico studies have been performed to analyze the influence of variants in the gene TMPRSS2 with regard to COVID-19 (Asselta et al., 2020; Barash et al., 2020; Fujikura and Uesaka, 2020; Ghafouri-Fard et al., 2020; Hou et al., 2020; Paniri et al., 2020; Russo et al., 2020; Senapati et al., 2020; Singh et al., 2020; Smatti et al., 2020; Vargas-Alarcón et al., 2020; Zarubin et al., 2020).
In a first small-scale study, Latini et al. (2020) observed a significantly lower frequency of rs12329760 T-allele carriers in Italian SARS-CoV-2-positive patients (N = 131) compared to a European reference group from the GnomAD database. Interestingly, we found a non-significant higher prevalence of T-allele carriers in the SARS-COV-2-positive patients (0.24) compared to the SARS-CoV-2-negative patients (0.20). This would make sense, if one assumes that T-allele carriers have a higher gene expression (GTEx portal), which could enhance TMPRSS2 activity and so may result in a higher infection susceptibility. Our observations would need to be analyzed in additional patients for this purpose.
This is the first case-control study analyzing the role of TMPRSS2 rs2070788 with regard to COVID-19. Although it seems plausible from a biological point of view, we did neither observe differences in allele or genotype frequencies between infected and non-infected patients nor different courses of the disease. The influence of this polymorphism has so far only been shown in a study with Chinese H1N1/H7N9 influenza A virus infected patients (Cheng et al., 2015). It cannot be excluded that this variant is relevant in a larger cohort or in other ethnic groups. Since the allele frequencies differ markedly between East Asians (0.36) and Europeans (0.46; Clarke et al., 2017), the latter hypothesis could well play a role. The variant rs2070788 was in high LD with rs383510 in the study with Chinese individuals by Cheng et al. (2015), which we did not observe in our German cohort. These findings fit the LD database records (Clarke et al., 2017) and again underscore the hypothesis that the variant rs2070788 might have a greater significance in other ethnicities.
Remarkably, we found a significant association to a 1.73-fold increased infection risk for TMPRSS2 rs383510 CC genotype carriers in our German cohort. Since genotyping errors can be excluded we have no explanation for the observation that TMPRSS2 rs383510 genotypes slightly violated HWE in the SARS-CoV-2-negative but not in the SARS-CoV-2-positive group. This could hypothetically be attributed to an unrecognized selection bias in the first-mentioned group. Multivariable logistic regression analysis confirmed the CC genotype as an independent risk factor for SARS-CoV-2 susceptibility. The mechanism behind our finding remains ambiguous or even counterintuitive. Considering the data from the GTEx database, it appears that the C-allele results in decreased gene expression in the lung (GTEx portal, 2020). Similarly, for H7N9 influenza A virus infected patients, it was observed that TT or TC genotype carriers had a higher infection risk, which seems quite plausible when viewed in conjunction with the expression data (Cheng et al., 2015). The GTEx data are based on post-mortem gene expression analysis in 515 individuals. Further functional studies need to be performed to validate the influence of the rs383510 polymorphism on mRNA or protein expression and function especially in the context of COVID-19. As mentioned above, our findings cannot be transferred to populations with different ethnicities and need to be validated in larger cohorts as well.
Regarding a genetic predisposition to a serious course of COVID-19, we did not observe any association with the analyzed TMPRSS2 variants. On the one hand, this could be due to the biological function of TMPRSS2 as a mediator of cell entry for SARS-CoV-2, so that differences in protein function might only influence the infection risk, rather than the course of the disease. On the other hand, the sample size of our cohort with “serious” COVID-19 was relatively small (N = 75). Further analyses are needed to conclusively clarify the influence of TMPRSS2 variants on disease severity.
In conclusion, for the first time, we report here an association of TMPRSS2 rs383510 CC genotype with SARS-CoV-2 infection susceptibility. Neither TMPRSS2 rs2070788 nor rs12329760 polymorphisms were related to SARS-CoV-2 infection risk or severity of COVID-19 in our German cohort. Before universally being acceptable, the role of variants in the gene TMPRSS2 in SARS-CoV-2 infections still needs to be elucidated in larger cohorts encompassing various ethnicities.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the Medical Faculty of the University of Duisburg-Essen (20-9230-BO). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
KrS: resources, methodology, data curation, formal analysis, validation, and writing – review and editing. KB: data curation and investigation. CE: resources and investigation. UD, DF, FH, JR, KaS, SS, and CT: investigation. K-HJ: investigation and validation. WS and AK: conceptualization, validation, supervision, and writing – review and editing. BM: conceptualization, resources, methodology, formal analysis, investigation, supervision, data curation, funding acquisition, project administration, visualization, writing – original draft, and writing – review and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by a grant (to BM) from the Stiftung Universitätsmedizin Essen of the Medical Faculty Essen. The funder of the study had no role in study design, data collection, data analysis, data interpretation, writing of the report, or in decision of submitting the paper for publication. We acknowledge support by the Open Access Publication Fund of the University of Duisburg-Essen.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We kindly thank Iris Manthey, Grit Müller, and Stephanie Büscher for technical support (Institute of Pharmacogenetics). We thank Katharina Jockers, Tanja Lehmann, and Heike Loewendick for their support regarding patient samples from the WBE (University Hospital Essen), Lothar Volbracht, Marc Wichert, Martina Wachsmann, and Petra Kasper for their support with the samples from the central laboratory (University Hospital Essen).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.667231/full#supplementary-material
References
Asselta, R., Paraboschi, E. M., Mantovani, A., and Duga, S. (2020). ACE2 and TMPRSS2 variants and expression as candidates to sex and country differences in COVID-19 severity in Italy. Aging 12, 10087–10098. doi: 10.18632/aging.103415
Barash, A., Machluf, Y., Ariel, I., and Dekel, Y. (2020). The pursuit of COVID-19 biomarkers: putting the spotlight on ACE2 and TMPRSS2 regulatory sequences. Front. Med. 7:582793. doi: 10.3389/fmed.2020.582793
Bertram, S., Dijkman, R., Habjan, M., Heurich, A., Gierer, S., Glowacka, I., et al. (2013). TMPRSS2 activates the human coronavirus 229E for cathepsin-independent host cell entry and is expressed in viral target cells in the respiratory epithelium. J. Virol. 87, 6150–6160. doi: 10.1128/JVI.03372-12
Cheng, Z., Zhou, J., To, K. K., Chu, H., Li, C., Wang, D., et al. (2015). Identification of TMPRSS2 as a susceptibility gene for severe 2009 pandemic A(H1N1) influenza and A(H7N9) influenza. J. Infect. Dis. 212, 1214–1221. doi: 10.1093/infdis/jiv246
Clarke, L., Fairley, S., Zheng-Bradley, X., Streeter, I., Perry, E., Lowy, E., et al. (2017). The international genome sample resource (IGSR): A worldwide collection of genome variation incorporating the 1000 genomes project data. Nucleic Acids Res. 45, D854–D859. doi: 10.1093/nar/gkw829
Ellinghaus, D., Degenhardt, F., Bujanda, L., Buti, M., Albillos, A., Invernizzi, P., et al. (2020). Genomewide association study of severe Covid-19 with respiratory failure. N. Engl. J. Med. 383, 1522–1534. doi: 10.1056/NEJMoa2020283
European Center of Disease Prevention and Control (ECDC) (2020). ECDC surveillance report. Available at: https://covid19-surveillance-report.ecdc.europa.eu/ (Accessed October 31, 2020).
FitzGerald, L. M., Agalliu, I., Johnson, K., Miller, M. A., Kwon, E. M., Hurtado-Coll, A., et al. (2008). Association of TMPRSS2-ERG gene fusion with clinical characteristics and outcomes: results from a population-based study of prostate cancer. BMC Cancer 8:230. doi: 10.1186/1471-2407-8-230
Fujikura, K., and Uesaka, K. (2020). Genetic variations in the human severe acute respiratory syndrome coronavirus receptor ACE2 and serine protease TMPRSS2. J. Clin. Pathol. doi: 10.1136/jclinpath-2020-206867 [Epub ahead of print]
Ghafouri-Fard, S., Noroozi, R., Vafaee, R., Branicki, W., Poṡpiech, E., Pyrc, K., et al. (2020). Effects of host genetic variations on response to, susceptibility and severity of respiratory infections. Biomed. Pharmacother. 128:110296. doi: 10.1016/j.biopha.2020.110296
GTEx portal (2020). The Genotype-Tissue Expression (GTEx) project. Available at: gtexportal.org (Accessed January 31, 2021).
Hoffmann, M., Kleine-Weber, H., Schroeder, S., Krüger, N., Herrler, T., Erichsen, S., et al. (2020). SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181, 271–280. doi: 10.1016/j.cell.2020.02.052
Hou, Y., Zhao, J., Martin, W., Kallianpur, A., Chung, M. K., Jehi, L., et al. (2020). New insights into genetic susceptibility of COVID-19: an ACE2 and TMPRSS2 polymorphism analysis. BMC Med. 18:216. doi: 10.1186/s12916-020-01673-z
Irham, L. M., Chou, W.-H., Calkins, M. J., Adikusuma, W., Hsieh, S.-L., and Chang, W.-C. (2020). Genetic variants that influence SARS-CoV-2 receptor TMPRSS2 expression among population cohorts from multiple continents. Biochem. Biophys. Res. Commun. 529, 263–269. doi: 10.1016/j.bbrc.2020.05.179
Latini, A., Agolini, E., Novelli, A., Borgiani, P., Giannini, R., Gravina, P., et al. (2020). COVID-19 and genetic variants of protein involved in the SARS-CoV-2 entry into the host cells. Gen. Dent. 11:1010. doi: 10.3390/genes11091010
Lukassen, S., Chua, R. L., Trefzer, T., Kahn, N. C., Schneider, M. A., Muley, T., et al. (2020). SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells. EMBO J. 39:e105114. doi: 10.15252/embj.20105114
Paniri, A., Hosseini, M. M., and Akhavan-Niaki, H. (2020). First comprehensive computational analysis of functional consequences of TMPRSS2 SNPs in susceptibility to SARS-CoV-2 among different populations. J. Biomol. Struct. Dyn. 1–18. doi: 10.1080/07391102.2020.1767690 [Epub ahead of print]
Paoloni-Giacobino, A., Chen, H., Peitsch, M. C., Rossier, C., and Antonarakis, S. E. (1997). Cloning of the TMPRSS2 gene, which encodes a novel serine protease with transmembrane, LDLRA, and SRCR domains and maps to 21q22.3. Genomics 44, 309–320. doi: 10.1006/geno.1997.4845
Russo, R., Andolfo, I., Lasorsa, V. A., Iolascon, A., and Capasso, M. (2020). Genetic analysis of the coronavirus SARS-CoV-2 host protease TMPRSS2 in different populations. Front. Genet. 11:872. doi: 10.3389/fgene.2020.00872
Senapati, S., Kumar, S., Singh, A. K., Banerjee, P., and Bhagavatula, S. (2020). Assessment of risk conferred by coding and regulatory variations of TMPRSS2 and CD26 in susceptibility to SARS-CoV-2 infection in human. J. Genet. 99, 1–5. doi: 10.1007/s12041-020-01217-7
Shirato, K., Kawase, M., and Matsuyama, S. (2013). Middle East respiratory syndrome coronavirus infection mediated by the transmembrane serine protease TMPRSS2. J. Virol. 87, 12552–12561. doi: 10.1128/JVI.01890-13
Shulla, A., Heald-Sargent, T., Subramanya, G., Zhao, J., Perlman, S., and Gallagher, T. (2011). A transmembrane serine protease is linked to the severe acute respiratory syndrome coronavirus receptor and activates virus entry. J. Virol. 85, 873–882. doi: 10.1128/JVI.02062-10
Singh, H., Choudhari, R., Nema, V., and Khan, A. A. (2020). ACE2 and TMPRSS2 polymorphisms in various diseases with special reference to its impact on COVID-19 disease. Microb. Pathog. 150:104621. doi: 10.1016/j.micpath.2020.104621
Smatti, M. K., Al-Sarraj, Y. A., Albagha, O., and Yassine, H. M. (2020). Host genetic variants potentially associated with SARS-CoV-2: a multi-population analysis. Front. Genet. 11:578523. doi: 10.3389/fgene.2020.578523
Vargas-Alarcón, G., Posadas-Sánchez, R., and Ramírez-Bello, J. (2020). Variability in genes related to SARS-CoV-2 entry into host cells (ACE2, TMPRSS2, TMPRSS11A, ELANE, and CTSL) and its potential use in association studies. Life Sci. 260:118313. doi: 10.1016/j.lfs.2020.118313
Zarubin, A., Stepanov, V., Markov, A., Kolesnikov, N., Marusin, A., Khitrinskaya, I., et al. (2020). Structural variability, expression profile, and pharmacogenetic properties of TMPRSS2 gene as a potential target for COVID-19 therapy. Gen. Dent. 12:19. doi: 10.3390/genes12010019
Keywords: transmembrane serine protease 2, coronavirus disease 2019, severe acute respiratory syndrome coronavirus type 2, polymorphism, rs383510, rs2070788, rs12329760, genetic association
Citation: Schönfelder K, Breuckmann K, Elsner C, Dittmer U, Fistera D, Herbstreit F, Risse J, Schmidt K, Sutharsan S, Taube C, Jöckel K-H, Siffert W, Kribben A and Möhlendick B (2021) Transmembrane serine protease 2 Polymorphisms and Susceptibility to Severe Acute Respiratory Syndrome Coronavirus Type 2 Infection: A German Case-Control Study. Front. Genet. 12:667231. doi: 10.3389/fgene.2021.667231
Edited by:
Georgios S. Markopoulos, University of Ioannina, GreeceReviewed by:
Byung-Hoon Jeong, Chonbuk National University, South KoreaGeorgios Ntritsos, University of Ioannina, Greece
Copyright © 2021 Schönfelder, Breuckmann, Elsner, Dittmer, Fistera, Herbstreit, Risse, Schmidt, Sutharsan, Taube, Jöckel, Siffert, Kribben and Möhlendick. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Birte Möhlendick, birte.moehlendick@uk-essen.de