TIMP metallopeptidase inhibitor 3 (TIMP-3) may contribute to the pathogenesis of venous thromboembolism (VTE). However, few studies have investigated the effect of TIMP-3 on VTE. Therefore, a two-sample Mendelian randomization (MR) analysis was conducted to investigate the association between TIMP-3 levels and VTE. Seven independent single-nucleotide polymorphisms (SNPs) for TIMP-3 levels were obtained from a published genome-wide association study (the KORA Consortium, including 997 Europeans). We obtained outcome datasets for VTE, pulmonary embolism (PE), and deep vein thrombosis (DVT) from the FinnGen Consortium. The primary analytical method used in the MR analysis was the inverse variance weighted (IVW) method. To enhance the robustness of the MR results, some other MR methods including weighted median, MR-Egger, and MR-PRESSO were conducted. Moreover, several sensitivity analyses were performed to identify potential horizontal pleiotropy and heterogeneity. In primary IVW MR analyses, per log increase in genetically predicted TIMP-3 levels were positively associated with the incidence of VTE (odds ratio [OR], 1.03; 95 % confidence interval (CI), 1.01, 1.06; P = 0.010), PE (OR, 1.04; 95 % CI, 1.01, 1.08; P = 0.009), and DVT (OR, 1.06; 95 % CI, 1.02, 1.10; P= 0.003). The results of the weighted median, MR-Egger, and MR-PRESSO were similar to the main findings. No unbalanced pleiotropy or heterogeneity was observed. The study suggests that genetically predicted high levels of TIMP-3 may be associated with an increased risk of VTE. These findings indicate that strategies targeting TIMP-3 may provide a basis for the prevention and treatment of VTE. Further investigation is required to clarify this potential mechanism.