Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2021 Mar 2:10:e67569.
doi: 10.7554/eLife.67569.

Association of Toll-like receptor 7 variants with life-threatening COVID-19 disease in males: findings from a nested case-control study

Chiara Fallerini #  1   2 , Sergio Daga #  1   2 , Stefania Mantovani #  3 , Elisa Benetti  2 , Nicola Picchiotti  4   5 , Daniela Francisci  6   7 , Francesco Paciosi  6   7 , Elisabetta Schiaroli  6 , Margherita Baldassarri  1   2 , Francesca Fava  1   2   8 , Maria Palmieri  1   2 , Serena Ludovisi  3   9 , Francesco Castelli  10 , Eugenia Quiros-Roldan  10 , Massimo Vaghi  11 , Stefano Rusconi  12   13 , Matteo Siano  12 , Maria Bandini  14 , Ottavia Spiga  5   15 , Katia Capitani  1   16 , Simone Furini  2 , Francesca Mari  1   2   8 ; GEN-COVID Multicenter Study  1 ; Alessandra Renieri  1   2   8 , Mario U Mondelli  3   9 , Elisa Frullanti  1   2
Collaborators, Affiliations

Association of Toll-like receptor 7 variants with life-threatening COVID-19 disease in males: findings from a nested case-control study

Chiara Fallerini et al. Elife. .

Abstract

Background: Recently, loss-of-function variants in TLR7 were identified in two families in which COVID-19 segregates like an X-linked recessive disorder environmentally conditioned by SARS-CoV-2. We investigated whether the two families represent the tip of the iceberg of a subset of COVID-19 male patients.

Methods: This is a nested case-control study in which we compared male participants with extreme phenotype selected from the Italian GEN-COVID cohort of SARS-CoV-2-infected participants (<60 y, 79 severe cases versus 77 control cases). We applied the LASSO Logistic Regression analysis, considering only rare variants on young male subsets with extreme phenotype, picking up TLR7 as the most important susceptibility gene.

Results: Overall, we found TLR7 deleterious variants in 2.1% of severely affected males and in none of the asymptomatic participants. The functional gene expression profile analysis demonstrated a reduction in TLR7-related gene expression in patients compared with controls demonstrating an impairment in type I and II IFN responses.

Conclusions: Young males with TLR7 loss-of-function variants and severe COVID-19 represent a subset of male patients contributing to disease susceptibility in up to 2% of severe COVID-19.

Funding: Funded by private donors for the Host Genetics Research Project, the Intesa San Paolo for 2020 charity fund, and the Host Genetics Initiative.

Clinical trial number: NCT04549831.

Keywords: COVID-19; LASSO Logistic Regression Analysis; TLR7; genetics; genomics; human; medicine.

PubMed Disclaimer

Conflict of interest statement

CF, SD, SM, EB, NP, DF, FP, ES, MB, FF, MP, SL, FC, EQ, MV, SR, MS, MB, OS, KC, SF, FM, AR, MM, EF No competing interests declared

Figures

Figure 1.
Figure 1.. Rare TLR7 variants and association with COVID-19.
LASSO logistic regression on boolean representation of rare variants of all genes of the X chromosome is presented. TLR7 is picked up by LASSO logistic regression as one of the most important genes on the X chr (Panel A). The LASSO logistic regression model provides an embedded feature selection method within the binary classification tasks (male patients with life-threatening COVID-19 vs infected asymptomatic male participants). The upward histograms (positive weights) reflect a susceptible behavior of the features to the target COVID-19, whereas the downward histograms (negative weights) a protective action. Panel B represents the cross-validation accuracy score for the grid of LASSO regularization parameters; the error bar is given by the standard deviation of the score within the 10 folds; the red circle (1.26) corresponds to the parameter chosen for the fitting procedure. Performances are evaluated through the confusion matrix of the aggregated predictions in the 10 folds of the cross-validation (Panel C) and with the boxplot (Panel D) of accuracy (60% average value), precision (59%), sensitivity (75%), specificity (43%), and ROC-AUC score (68%). The box extends from the Q1 to Q3 quartile, with a line at the median (Q2) and a triangle for the average.
Figure 2.
Figure 2.. Gene expression profile analysis in peripheral blood mononuclear cells (PBMCs) and in HEK293 cells transfected with the functional variants after stimulation with a TLR7 agonist for 4 hr.
(A) 5 × 105 PBMCs from COVID-19 patients and six unaffected male and female controls were stimulated for 4 hr with the TLR7 agonist imiquimod at 5 μg/mL or cell culture medium. Quantitative PCR assay was performed and the 2-ΔΔCt calculated using HPRT1 as housekeeping gene. Fold change in mRNA expression of TLR7 and type 1 IFN-related genes ISG15, IRF7, IFN-ɑ and IFN-γ induced by TLR7 agonist imiquimod was compared with cell culture medium. Ctl indicates healthy controls (white bar); C1, the asymptomatic mutated control (diagonal lines bar); P2, P5, cases with neutral variants (vertical lines bar); P1, P3, P8, P7 cases with functional variants (gray bar) (as in Table 2). (B) Histograms of intracellularly expressed TLR7 protein in HEK293 cells transfected with the different TLR7 plasmids. (C) Gene expression profile analysis of IFN-ɑ in transfected cells after stimulation with the TLR7 agonist imiquimod. WT indicates cells transfected with WT TLR7 plasmid. Quantitative PCR assay was performed and the 2-ΔΔCt calculated using HPRT1 as housekeeping gene. Fold change in mRNA expression induced by imiquimod was compared with cell culture medium. Error bars show standard deviation. p values were calculated for the reduction using an unpaired t test: *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001.
Figure 3.
Figure 3.. Segregation analysis.
Fold change in mRNA expression following Imiquimod stimulation of TLR7 itself and its main effectors, IRF7, ISG15, IFN-alpha, and IFN-gamma is shown in Panel A. Gray columns represent individuals harboring the TLR7 variant and black columns are severely affected SARS-CoV-2 cases. Pedigree (Panel B) and respective segregation of TLR7 variant and COVID-19 status (Panel C) are also shown. Squares represent male family members; circles, females. Individuals infected by SARS-CoV-2 are indicated by a virus cartoon close to the individual symbol (formula image).
See this image and copyright information in PMC

Comment in

References

    1. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nature Methods. 2010;7:248–249. doi: 10.1038/nmeth0410-248. - DOI - PMC - PubMed
    1. Bastard P, Rosen LB, Zhang Q, Michailidis E, Hoffmann HH, Zhang Y, Dorgham K, Philippot Q, Rosain J, Béziat V, Manry J, Shaw E, Haljasmägi L, Peterson P, Lorenzo L, Bizien L, Trouillet-Assant S, Dobbs K, de Jesus AA, Belot A, Kallaste A, Catherinot E, Tandjaoui-Lambiotte Y, Le Pen J, Kerner G, Bigio B, Seeleuthner Y, Yang R, Bolze A, Spaan AN, Delmonte OM, Abers MS, Aiuti A, Casari G, Lampasona V, Piemonti L, Ciceri F, Bilguvar K, Lifton RP, Vasse M, Smadja DM, Migaud M, Hadjadj J, Terrier B, Duffy D, Quintana-Murci L, van de Beek D, Roussel L, Vinh DC, Tangye SG, Haerynck F, Dalmau D, Martinez-Picado J, Brodin P, Nussenzweig MC, Boisson-Dupuis S, Rodríguez-Gallego C, Vogt G, Mogensen TH, Oler AJ, Gu J, Burbelo PD, Cohen JI, Biondi A, Bettini LR, D'Angio M, Bonfanti P, Rossignol P, Mayaux J, Rieux-Laucat F, Husebye ES, Fusco F, Ursini MV, Imberti L, Sottini A, Paghera S, Quiros-Roldan E, Rossi C, Castagnoli R, Montagna D, Licari A, Marseglia GL, Duval X, Ghosn J, Tsang JS, Goldbach-Mansky R, Kisand K, Lionakis MS, Puel A, Zhang SY, Holland SM, Gorochov G, Jouanguy E, Rice CM, Cobat A, Notarangelo LD, Abel L, Su HC, Casanova JL, HGID Lab, NIAID-USUHS Immune Response to COVID Group, COVID Clinicians, COVID-STORM Clinicians, Imagine COVID Group, French COVID Cohort Study Group, Milieu Intérieur Consortium, CoV-Contact Cohort, Amsterdam UMC Covid-19 Biobank, COVID Human Genetic Effort, Autoantibodies against type I IFNs in patients with life-threatening COVID-19. Science. 2020;370:eabd4585. doi: 10.1126/science.abd4585. - DOI - PMC - PubMed
    1. Chehadeh W, Alkhabbaz M. Differential TLR7-mediated expression of proinflammatory and antiviral cytokines in response to laboratory and clinical Enterovirus strains. Virus Research. 2013;174:88–94. doi: 10.1016/j.virusres.2013.03.006. - DOI - PubMed
    1. Daga S, Fallerini C, Baldassarri M, Fava F, Valentino F, Doddato G, Benetti E, Furini S, Giliberti A, Tita R, Amitrano S, Bruttini M, Meloni I, Pinto AM, Raimondi F, Stella A, Biscarini F, Picchiotti N, Gori M, Pinoli P, Ceri S, Sanarico M, Crawley FP, Birolo G, Renieri A, Mari F, Frullanti E, GEN-COVID Multicenter Study Employing a systematic approach to biobanking and analyzing clinical and genetic data for advancing COVID-19 research. European Journal of Human Genetics. 2021;1:1–15. doi: 10.1038/s41431-020-00793-7. - DOI - PMC - PubMed
    1. Garg S, Kim L, Whitaker M, O'Halloran A, Cummings C, Holstein R, Prill M, Chai SJ, Kirley PD, Alden NB, Kawasaki B, Yousey-Hindes K, Niccolai L, Anderson EJ, Openo KP, Weigel A, Monroe ML, Ryan P, Henderson J, Kim S, Como-Sabetti K, Lynfield R, Sosin D, Torres S, Muse A, Bennett NM, Billing L, Sutton M, West N, Schaffner W, Talbot HK, Aquino C, George A, Budd A, Brammer L, Langley G, Hall AJ, Fry A. Hospitalization rates and characteristics of patients hospitalized with Laboratory-Confirmed coronavirus disease 2019 - COVID-NET, 14 states, march 1-30, 2020. MMWR. Morbidity and Mortality Weekly Report. 2020;69:458–464. doi: 10.15585/mmwr.mm6915e3. - DOI - PMC - PubMed

Publication types

Associated data