Background: Thyroid dysfunction has been observed among some patients with COVID-19. It is unclear whether SARS-CoV-2 infection (or its severity) leads to the development of thyroid dysfunction, or vice versa. In this study, we examined the bi-directional causal relationship between host genetic liability to three COVID-19 phenotypes (including SARS-CoV-2 infection, hospitalized and severe COVID-19) and three thyroid dysfunction traits [including hyperthyroidism, hypothyroidism, and autoimmune thyroid disease (AITD)] and three continuous traits of thyroid hormones [including thyroid-stimulating hormone (TSH) and free thyroxine (FT4) within reference range, and TSH in full range]. Methods: Summary statistics from the largest available meta-analyses of human genome-wide association studies were retrieved for the following variables: SARS-CoV-2 infection (n=1,348,701), COVID-19 hospitalization (n=1,557,411), severe COVID-19 (n=1,059,456), hyperthyroidism (n=51,823), hypothyroidism (n=53,423), AITD (n=755,406), TSH within reference range (n=54,288), FT4 within reference range (n=49,269), and TSH in full range (n=119,715). Using a two-sample Mendelian randomization (MR) approach, the inversevariance weighted (IVW) method was adopted as the main MR analysis. Weighted median, contamination mixture, MR-Egger, and MR-PRESSO methods were applied as sensitivity analyses.Results: Host genetic susceptibility to SARS-CoV-2 infection was causally associated with hypothyroidism in the main IVW analysis (per doubling in prevalence of SARS-CoV-2 infection, Odds Ratio (OR)=1.335; 95% CI: 1.167-1.526; p=2.4x10 -5 , surpassing the Bonferroni multiple-testing threshold). Similar causal estimates were observed in the sensitivity analyses (weighted median: