Study population in NHANES
This study included data from six continuous cycles of the NHANES database (2007–2008, 2009–2010, 2011–2012, 2013–2014, 2015–2016, 2017–2018). We included 17,907 female participants aged ≥ 20 years with the following exclusion criteria: (1) Being pregnant (n=274), (2) Lack of informative data on diagnosis of hypertension (n=3), (3) insufficient hysterectomy data (n=2646), (4) lack of necessary demographic information and health-related data (n=2356). Finally, 12,628 participants with complete NHANES data were included in this work. Fig. 1 displays the participant selection procedure.
Exposure and Outcomes
Information about hysterectomy status was obtained in the form of a questionnaire. Hysterectomy status was defined using self-reported history by asking the following question in the reproductive health section:“Have you had a hysterectomy?”In line with American Heart Association/American College of Cardiology (AHA/ACC) 2017 guidelines and previous studies [30,31], The diagnosis of hypertension was determined by a questionnaire about hypertension and the average of four consecutive blood pressure measurements. Participants were diagnosed with hypertension if they had self-reported hypertension, DBP ≥ 90 mmHg, SBP ≥ 140 mmHg, or previous or current medication for hypertension.
Covariates
Covariates in the present work included demographic features (age, education, race, marital status, body mass index (BMI), poverty/income ratio (PIR), drinking and smoking), medical history (diabetes, oophorectomy, hormone therapy), laboratory tests (total cholesterol (TC), total glyceride (TG), high-density lipoprotein cholesterol (HDL), uric acid (UA), creatinine (Cr), glycosylated hemoglobin (GHb)) and dietary data (sodium intake).
PIR was computed through the division of family (or individual) income by poverty criteria applicable to survey year, with the greater PIR indicating superior family income status. The PIR for the not poor was defined as ≥ 1, and for the poor, it was defined as < 1. The diagnosis of diabetes was made according to a self-reported diabetes history, glycosylated hemoglobin ≥ 6.5%, fasting blood glucose ≥ 126 mg/dl, insulin or antihyperglycemic agents medication [32]. Dietary sodium intake was defined as mean sodium intake of two 24-h recalls.
Statistical analysis for NHANES analysis
When performing the NHANES analyses, sampling weights provided by NHANES were utilized for weighting. Continuous variables were represented by mean±standard deviation (SD), whereas categorical variables were represented by numbers (percentages). Relationship of hysterectomy status with hypertension risk was evaluated using multivariate logistic regression models. Model 1 was not adjusted for all covariates. Model 2 was adjusted for age, race, educational level, marital status, PIR. In Model 3, BMI, smoking, alcohol, diabetes, oophorectomy, hormone therapy, levels of TC, HDL, TG, UA, Cr, GHb and sodium intake was adjusted based on Model 2. For determining whether the relation of hysterectomy status with hypertension risk was different among different subgroups, we carried out the independent stratification. Additionally, we utilized Wald test for calculating p-value for interaction. The restricted cubic splines (RCS) analysis was constructed for exploring association nonlinearity and depicting the overall trends. R version 4.1.3 was adopted in statistical analysis. p-value < 0.05 (two-sided) revealed significant difference.
Assumptions and data sources of two‑sample MR
Hysterectomy exposure data were obtained based on MRC Integrated Epidemiology Unit (MRC-IEU) Consortium, involving 462,933 Europeans (46,411 patients and 416,522 control participants) in total (ID:ukb-b-3700). Summary data of hypertension were acquired based on a large UK Biobank-based GWAS [33]. Researchers used ICD-10 codes, UK Biobank self-reported disease diagnosis to identify individuals with a history of hypertension. Through selection, this GWAS enrolled 129,909 cases and 354,689 controls in total (ID:ebi-a-GCST90038604).
In addition, we also included data on SBP (ID:ieu-b-38) and DBP (ID:ieu-b-39) as outcomes, which were acquired based on the meta-analysis of GWAS involving 757,601 subjects from UK Biobank and International Council on Blood Pressure (ICBP) associations [34]. We obtained all the above GWAS data from Integrated Epidemiology Unit (IEU) OpenGWAS database (https://gwas. mrcieu. ac. uk/). Data download details can be observed from Additional file 1: Table S1.
Selection of genetic instruments
IVs used for MR analysis must satisfy three requirements: (1) the IVs should be convincingly correlated to hysterectomy; (2) IVs are not related to all confounders for exposure-outcome relation; and (3) are just related to outcomes via the interested exposure rather than via additional routes.
To identify genetic variants for the causality of hysterectomy status with risk of developing several outcomes (hypertension, increased SBP, and DBP), we set genome-wide significance level as p < 5 × 10-8 for screening for genetic variants closely related to exposure. Afterward, the linkage disequilibrium (LD) clumping test was conducted for identifying independent SNP (r2 < 0.001; 10,000 kb). After excluding palindrome SNPs, effect alleles were coordinated in outcome and exposure datasets. We then examined SNPs via LDtrait Tool (https://ldlink.nih.gov/?tab=ldtrait) for removing those closely related to potential confounders such as BMI, waist circumference, diabetes, dyslipidemia [15, 35]. Subsequently, to further assess the strength of each IV, F-statistic was determined for the IVs in the exposure and excluded SNPs with F < 10 to ensure that the IVs had adequate validity and instrumental strength [36].
Statistical analysis for two‑sample MR
In two-sample MR study, random-effects inverse-variance weighting (IVW) was used to be the primary method for assessing causal relationships of genetically predicted hysterectomy status with hypertension risk and increased SBP, and DBP. Fixed-effects estimates are not appropriate when there is excessive result heterogeneity; consequently, random-effects IVW model was adopted in our analysis, accounting for heterogeneity when assessing causality [37]. In addition, we used four complementary MR analysis approaches, such as simple mode, MR‒Egger, weighted mode, and weighted median approaches, to validate IVW results. MR‒Egger regression intercept was adopted for assessing potential horizontal multiplicity effects. Horizontal pleiotropy was considered not to be present at p-value >0.05 [36]. Cochrane's Q test was utilized for testing for possible SNP heterogeneity. Analyses with p-values > 0.05 did not reveal any obvious heterogeneity. Heterogeneity does not invalidate causality estimates in MR analyses because the random effects IVW approach balances potential total heterogeneity to some extent [38].
We also performed MR pleiotropy residual sum and outlier (MR-PRESSO) test in determining existence of outlier IVs. In addition, residual sum of squares accounts for a heterogeneity measure, which equals Cochran's Q statistic. In the presence of outlier IVs, we utilized the MR-PRESSO outlier-corrected test for obtaining corrected causal effects through removing outliers, and then assessed distortion of causal estimates prior to and following removing outliers with MR-PRESSO distortion test. Typically, MR-PRESSO test is useful if few genetic variants had heterogeneous ratio estimates, because they were excluded and therefore did not impact total estimate [39, 40]. Additionally, leave one-out analyses were carried out for assessing whether single-sensitive SNP affected IVW test. The results were also validated using funnel plots and scatter plots for further validation. Two-sample MR package (version 0.5.6) in R (version 4.1.3) was employed for analysis.
Fig. 1 shows the research flowchart in the present study.